Review Special Issues

Biofilms at work: Bio-, phyto- and rhizoremediation approaches for soils contaminated with polychlorinated biphenyls

  • Organohalide contaminants such as polychlorinated biphenyls (PCBs) have been released into the environment for decades due to anthropogenic activities, but are also naturally produced in small amounts through volcanic eruptions and geochemical processes. Although toxic to humans and other organisms, the natural production of these compounds has resulted in the evolution of naturally occurring organohalide-respiring bacteria that possess the enzymes necessary to degrade PCB compounds to non-toxic products. The efficiency of PCB degradation can be improved by facilitating the formation of organohalide-respiring biofilms. During biofilm colonization on a surface or interface, bacteria are encased in an extracellular polymeric substance (EPS) or “slime,” which allows them to share nutrients and remain protected from environmental stresses. Effective bioremediation of PCBs involves facilitation of biofilm growth to promote cooperation between bacteria, which can be further enhanced by the presence of certain plant species. This review aims to give an overview of biofilm processes involved in the detoxification of PCBs including anaerobic and aerobic PCB degradation by bacteria as well as the ability of plants to stimulate microbial activity and degradation (rhizoremediation and phytoremediation).

    Citation: Merily Horwat, Meggie Tice, Birthe V. Kjellerup. Biofilms at work: Bio-, phyto- and rhizoremediation approaches for soils contaminated with polychlorinated biphenyls[J]. AIMS Bioengineering, 2015, 2(4): 324-334. doi: 10.3934/bioeng.2015.4.324

    Related Papers:

    [1] Martha Garlick, James Powell, David Eyre, Thomas Robbins . Mathematically modeling PCR: An asymptotic approximation with potential for optimization. Mathematical Biosciences and Engineering, 2010, 7(2): 363-384. doi: 10.3934/mbe.2010.7.363
    [2] Ellina Grigorieva, Evgenii Khailov, Andrei Korobeinikov . Parametrization of the attainable set for a nonlinear control model of a biochemical process. Mathematical Biosciences and Engineering, 2013, 10(4): 1067-1094. doi: 10.3934/mbe.2013.10.1067
    [3] Paula Mercurio, Di Liu . Identifying transition states of chemical kinetic systems using network embedding techniques. Mathematical Biosciences and Engineering, 2021, 18(1): 868-887. doi: 10.3934/mbe.2021046
    [4] Xiaoxuan Pei, Kewen Li, Yongming Li . A survey of adaptive optimal control theory. Mathematical Biosciences and Engineering, 2022, 19(12): 12058-12072. doi: 10.3934/mbe.2022561
    [5] Huangjing Yu, Yuhao Wang, Heming Jia, Laith Abualigah . Modified prairie dog optimization algorithm for global optimization and constrained engineering problems. Mathematical Biosciences and Engineering, 2023, 20(11): 19086-19132. doi: 10.3934/mbe.2023844
    [6] Yancong Zhou, Chenheng Xu, Yongqiang Chen, Shanshan Li . C4 olefin production conditions optimizing based on a hybrid model. Mathematical Biosciences and Engineering, 2023, 20(7): 12433-12453. doi: 10.3934/mbe.2023553
    [7] Qing Wu, Chunjiang Zhang, Mengya Zhang, Fajun Yang, Liang Gao . A modified comprehensive learning particle swarm optimizer and its application in cylindricity error evaluation problem. Mathematical Biosciences and Engineering, 2019, 16(3): 1190-1209. doi: 10.3934/mbe.2019057
    [8] Rong Zheng, Heming Jia, Laith Abualigah, Shuang Wang, Di Wu . An improved remora optimization algorithm with autonomous foraging mechanism for global optimization problems. Mathematical Biosciences and Engineering, 2022, 19(4): 3994-4037. doi: 10.3934/mbe.2022184
    [9] Liping Wu, Zhongyi Xiang . A study of integrated pest management models with instantaneous and non-instantaneous impulse effects. Mathematical Biosciences and Engineering, 2024, 21(2): 3063-3094. doi: 10.3934/mbe.2024136
    [10] Andriy A. Avramenko, Igor V. Shevchuk . Renormalization group analysis of heat transfer in the presence of endothermic and exothermic chemical reactions. Mathematical Biosciences and Engineering, 2019, 16(4): 2049-2062. doi: 10.3934/mbe.2019100
  • Organohalide contaminants such as polychlorinated biphenyls (PCBs) have been released into the environment for decades due to anthropogenic activities, but are also naturally produced in small amounts through volcanic eruptions and geochemical processes. Although toxic to humans and other organisms, the natural production of these compounds has resulted in the evolution of naturally occurring organohalide-respiring bacteria that possess the enzymes necessary to degrade PCB compounds to non-toxic products. The efficiency of PCB degradation can be improved by facilitating the formation of organohalide-respiring biofilms. During biofilm colonization on a surface or interface, bacteria are encased in an extracellular polymeric substance (EPS) or “slime,” which allows them to share nutrients and remain protected from environmental stresses. Effective bioremediation of PCBs involves facilitation of biofilm growth to promote cooperation between bacteria, which can be further enhanced by the presence of certain plant species. This review aims to give an overview of biofilm processes involved in the detoxification of PCBs including anaerobic and aerobic PCB degradation by bacteria as well as the ability of plants to stimulate microbial activity and degradation (rhizoremediation and phytoremediation).


    It is prevalent to solve the optimization problems in our life, which is generally done by selecting a set of variables to reach an optimal value while satisfying a set of constraints. For instance, antenna arrays are a typical global optimization problem that is important for detecting and processing signals from different directions. The goal of antenna array synthesis is to generate a physical layout of the array in which the radiation pattern is as close as possible to the desired pattern. Unfortunately, most optimization problems are always nonlinear, non-differentiable, and multipolar. Numerous optimization methods have been proposed for solving these optimization problems. Classical optimization methods, including gradient methods and alternating projection methods, are constantly being proposed. Nevertheless, it is known that classical optimization methods require a fixed point that is fairly near to the final value, otherwise they are likely to get trapped in local optimal. Moreover, the quality of the solution strongly depends on the estimation of the initial values. However, the goodness of the final solution scheme of the algorithm mostly depends on the setting of the initial scheme. If the initial solution falls in the worse region of the solution, then the search will only be conducted to find the best solution in the worse solution region. Considering the shortcomings of these techniques, evolutionary algorithms are the best choice for solving large-scale optimization problems since evolutionary algorithms are population-based search algorithms with high convergence rates and have been successfully used for large-scale optimization problems extensively [1,2]. Therefore, different kinds of evolutionary algorithms, such as particle swarm optimization algorithms (PSO) [3,4], ant colony optimization (ACO), biogeography-based optimization (BBO) [5,6,7,8], animal migration optimization (AMO) [9,10], cuckoo search algorithm (CS) [11,12,13], and artificial bee colony (ABC) [14], have been adopted to handle antenna arrays, especially on large-scale optimization problems.

    Recently, chemical reaction optimization (CRO) [15] is a population-based heuristic evolutionary algorithm that simulates molecular collisions in chemical reactions. The algorithm simulates the intermolecular interactions of molecules from a high-energy state into a low-energy state to explore the global optimal solution in the search space. To achieve efficient performance and extensive applications of the algorithm, different variations of the CRO algorithm were proposed successively to enhance the performance of the traditional CRO algorithm. In addition, there are many works that focus on the applications of CRO algorithms; for instance, Lam and Li [16] developed a real-coded CRO, called RCCRO for addressing continuous optimization problems. The experimental results indicate that the average performance of RCCRO is superior to other algorithms. In [17], the authors developed a double molecular structure-based chemical reaction optimization (DMSCRO) method and the CRO scheme is employed to outline the directed acyclic graph jobs scheduling in heterogeneous computing systems. Xu et al. [18] developed a hybrid chemical reaction optimization framework to solve the DAG-based task scheduling problem. It combines the novel heuristic approaches and a proposed selection strategy. Dam et al. [19] proposed a new method that integrates the chemical reaction optimization scheme with the unified tabu search heuristic method for solving the capacitated vehicle routing problem.

    Lam and Li [20] solved the peer to peer live streaming with chemical reaction optimization and experimental results demonstrate that CRO beats many commonly-used techniques for containing population transitions in the system applications. In [21], Lam and Li used chemical reaction optimization to resolve the grid scheduling problem. Experimental results indicate that the proposed algorithm surpasses many existing evolutionary methods for the majority of cases tested. Truong et al. [22] introduced a chemical reaction optimization algorithm (CROG) coupled with a greedy strategy to tackle the quadratic allocation problem. The algorithm also proposes a design method for the update operator and a parameter tuning method. Unfortunately, it is difficult for these approaches to balance the exploitation and exploration capability of CRO algorithms simultaneously. Therefore, it is necessary to develop new optimization strategies for chemical reaction optimization to further enhance the performance of the algorithms.

    To address these issues, we propose a new realistic encoded chemical reaction optimization algorithm that can boost the performance of traditional chemical reaction algorithms. First, inspired by the search mechanism in differential evolution algorithms, we propose a modified chemical reaction optimization (MCRO) model based on the optimal solution in the population. Then, to further improve the performance of the algorithm, the perturbation frequency of individuals is manipulated by increasing the modification rate of individuals and balance the exploitation and exploration ability of the algorithm. Moreover, to further accelerate the convergence of the algorithm, we assemble the stochastic initialization, opposition learning methods, and generalized opposition learning methods to obtain a new population initialization method. Finally, to demonstrate the superiority of the proposed algorithm, the algorithm is compared with several benchmark functions and applied on two antenna arrays. Based on the experimental results, it can be observed that our proposed algorithm is superior to other state-of-the-art algorithms from several perspectives.

    Chemical Reaction Optimization (CRO) is an evolutionary algorithm first published by Lam and Li in [22,23], which simulates the collision of molecules in a chemical reaction. For each molecule, its molecular structure is similar to the optimization problem and the molecule denotes the solution to the problem under consideration. Potential energy ($ \chi $) is the objective function value and kinetic energy ($ \eta $) can be interpreted as the tolerance of the system, which accepts a solution that is worse than the current one. This approach can assist the solution to be partially replaced locally. Assuming that the numerator is $ \omega $, $ \chi $ can be calculated as follows: $ \chi_{\omega} = f(\omega) $.

    In a standard CRO, the algorithm simulates four different types of chemical reactions, comprising on-wall ineffective collision, decomposition, inter-molecular ineffective collision, and synthesis. They are mainly employed to control the redistribution of energy between molecules and buffers in the solution. Among them, single operations include wall ineffective collisions and disintegration, and the rest are multiple operations.

    When a molecule strikes the wall of a container and then bounces away, the on-wall ineffective collision reaction is invoked. The process will have some change. The molecular structure $ \omega $ and $ \chi $ will change to a generated molecular $ \omega^{\prime} $ in the neighborhood of $ \omega $ if the condition is satisfied under Eq 1.

    $ χω+ηωχω
    $
    (1)

    After an invalid collision at the wall, a certain fraction of the $ \eta $ of the transformed molecule is withdrawn to the buffer. Then, the $ \eta $ of the new molecule $ \omega^{\prime} $ can be as follows:

    $ ηω=(χω+ηωηω)×a
    $
    (2)

    where $ a \in[\eta^{LR}, 1] $, $ \eta^{LR} $ is a statistical parameter that limits the maximum percentage of $ \eta $ that can be lost at each time. As long as there is enough $ \chi $ at the beginning, it is possible for a molecule with lower $ \chi $ to transform into a molecule with higher $ \chi $, corresponding to a poorer solution. Following the collision, the molecule has fewer $ \eta $.

    The decomposition collision reaction results from a molecule hitting a wall, which then breaks one molecule into several molecules. Assume that $ \omega $ produces $ \omega_{1}^{\prime} $ and $ \omega_{2}^{\prime} $, i.e., $ \omega \rightarrow \omega_{1}^{\prime}+\omega_{2}^{\prime} $. Any mechanism, which can generate $ \omega_{1}^{\prime} $ and $ \omega_{2}^{\prime} $ from $ \omega $, is allowed. The algorithm considers two cases for the decomposition reaction: first, each molecule needs enough energy to complete the decomposition. The $ \chi $ of the resulting molecule is denoted by the following:

    $ ηω1=(χω+ηωχω1χω2)×qηω2=(χω+ηωχω1χω2)×(1q)
    $
    (3)

    Based on the above analysis, the following conditions need to be satisfied for molecular decomposition:

    $ χω+ηωχω1+χω2
    $
    (4)

    where $ q $ is a random number that takes values in the range [0, 1]. Meanwhile, the algorithm considers another condition:

    $ ηω1=(χω+ηωχω1χω2+ buffer )×k1×k2ηω2=(χω+ηωχω1χω2+ buffer )×k3×k4
    $
    (5)
    $ χω+ηω+ buffer χω1+χω2
    $
    (6)

    where $ k_{1}, k_{2}, k_{3} $ and $ k_{4} $ are random numbers which take values in the range of [0, 1]. The buffer can be defined as below:

    $  buffer =χω+ηωχω1χω2+bufferηω1ηω2
    $
    (7)

    Null collisions between molecules happen when two or more molecules collide with each other and then separate. The energy conservation condition can be expressed in the following way:

    $ χω+ηω+ηω1+ηω2χω1+χω2
    $
    (8)

    Then the $ \eta $ of two transformed molecules can be described as follows:

    $ ηω0=(χω+ηω+ηω1+ηω2χω1χω2)×qηω2=(χω+ηω+ηω1+ηω2χω1χω2)×(1q)
    $
    (9)

    where $ q $ is a random number that takes values in the range [0, 1].

    The synthesis is the opposite of decomposition. Synthesis occurs when a molecule collides and bond together. The reaction conditions that allow the synthesis reaction are description below:

    $ χω+ηω+ηω1+ηω2χω
    $
    (10)

    Then the $ \eta $ for the new molecule is considered by the following method.

    $ ηω=χω+ηω+ηω1+ηω2χω
    $
    (11)

    To address the continuous optimization problem more efficiently, Lam and Li proposed a real coded version of CRO, namely RCCRO. Subsequently, a similar modified neighborhood search approach was used. Among them, the Gaussian distribution model was proposed to generate perturbations to search continuous space neighborhood. Their update method can be characterized here below:

    $ ω(i)=ω(i)+N(0,σ2)
    $
    (12)

    where $ \sigma $ holds the step size and determines the efficiency of the algorithm. Because different problems have different characteristics, various $ \sigma $ should be adopted. In the algorithm, the neighborhood method will be employed for one-wall and intermolecular ineffective collision. The framework of the algorithm is shown in Figure 1. In summary, the CRO algorithm starts with the initialization. Then if it is an inter-molecular ineffective collision, the single molecule reaction is employed. Otherwise we will carry out the bimolecule reaction. In terms of the single molecule reaction, if the composition condition is satisfied, then it will conduct the decomposition, otherwise, it will conduct the on-wall ineffective collision. In terms of the bimolecular reaction, if the synthesis condition is met, the synthesis reaction is conducted, otherwise the inter-molecular ineffective collision reaction is conducted. After these steps, new minimum points are checked. If the stopping criterion is not matched, then the above steps will be repeated. Otherwise, the global minimum point is obtained.

    Figure 1.  The framework of the CRO algorithm.

    For evolutionary algorithms, population initialization is a crucial component. The goodness of population initialization directly influences the final convergence speed of the algorithm and the quality of the optimal solution. In fact, most algorithms use random initialization to generate the initial population in the absence of a priori information. Moreover, it is difficult to guarantee the quality of initialized population solutions by individual population initialization methods. Therefore, this paper proposes a new initialization method to generate the initial population by fusing random initialization, opposition-based learning initialization [24] and generalized opposition learning initialization [25]. This method can enhance the performance of CRO by preventing premature convergence to local optimal. Table 4 summarized the framework of this method. Figure 2 depicted the difference between the random population initiation and the opposition-based population initialization discussed above.

    Figure 2.  The difference between the random population initiation and above mentioned opposition based population initialization.
    Procedure A new initialization approach
    Begin
      Randomly initialize each individual in population x.
      Set the individual counter i=1, j=1
      for $ i=1 \rightarrow NP $ do
        for $ j=1 \rightarrow D $ do
          $ o x_{i, j}=x_{\min, j}+x_{\max, j}-x_{i, j} $
        end for
      end for
      Set the individual counter $ i=1 $, $ j=1 $
      for $ i=1 \rightarrow NP $ do
        for $ j=1 \rightarrow D $ do
          $ oo x_{i, j}=x_{\min, j}+x_{\max, j}-x_{i, j} $
        end for
      end for
      The fitness of each individual was assessed and the most suitable individual for NP from $ \{X(N P) \cup O X(N P) \cup O O X(N P)\} $ was selected as the initial population.
    End

     | Show Table
    DownLoad: CSV

    Differential evolution algorithm (DE) [26] is an efficient stochastic search method which was originally proposed mainly for solving continuous optimization problems. Its simple structure and ease of completion allowed the algorithm to be quickly applied to different optimizations and achieve impressive performance. The algorithm also has the identical evolutionary process as other evolutionary algorithms. In the algorithm, we can see that the algorithm consists of three main processes: mutation, crossover and selection. At the beginning of the algorithm, a mid-autumn is generated using a random initialization, followed by the operations of mutation, crossover and selection to generate a new population. It is worth noting that several different variation operations have been proposed in the DE algorithm, such as:

    $  DE/best /1:Vi=Xbest +F(Xr1Xr2)
    $
    (13)

    where $ i \in\{1, \cdots, N P\} $ and $ r_{1} $ and $ r_{2} $ are random integer indices selected from $ \{1, \cdots, N P\} $ that are different from each other. $ F $ is a scaling factor or amplification factor. $ X_{\text{best }} $, the base vector to be perturbed, is the best individual in the current population and indicates the best information in the population.

    In fact, the information about the best individuals in the current population is very useful for the whole population, and if the information is used wisely it is possible to accelerate the convergence rate of the population. For the strategy "DE/best/1", the best solution can explore the movement of the current population. Based on DE algorithm and the property of CRO, we propose a novel search mechanism to improve CRO:

    $ vij=xbestj+φij(xijxbj)
    $
    (14)

    where $ \varphi_{i j} $ denotes the scaling factor that controls the speed of optimization of the algorithm. It is calculated at the beginning of each generation. If $ rand < rand $, then $ \varphi_{i j} $ is equal to $ 5 \times ({rand-rand }){\times}(rand-rand) $, otherwise, $ \varphi_{i j} $ equals to $ randn $.

    We can note that "$ rand $" denotes a random number generated from the range of $ [-1, 1], i \in\{1, 2, \cdots, NP\} $, and $ j \in\{1, 2, \cdots, D\} $ are randomly chosen indexes, "$ randn $" denotes a number generated through a Gaussian distribution with mean "0" and standard deviation "1".

    To control the frequency of perturbation, it makes the algorithm can solve the difficult problem better. A control parameter, the modification rate (MR), is proposed, which serves mainly to generate a candidate position $ v_{i j} $ from the current memory of $ x_{i j} $ as described below:

    $ vij={xbest +φij(xijxbj), if rijMRxij otherwise 
    $
    (15)

    where $ \mathrm{r}_{\mathrm{ij}} $ denotes a random number generated from the range of $ [-1, 1] $. This new rule is similar with the binomial crossover of the DE algorithm, which can enhance the potential diversity of the population.

    For evolutionary algorithms, the setting of bounds is crucial. During the search process, the individuals in the part of MCRO algorithm can easily move to the boundary of the whole search space, which renders many of the solutions invalid. Therefore, it is important to restrict the individuals around the boundary to the search space once again. In other words, if most of the individuals are on the boundary, it is difficult for the algorithm to find the optimal solution to the problem, which will necessarily make the algorithm fall into local optimal solutions and make the whole population very easy to lose the population diversity. To address this problem, we propose to use the following rule for the setting of off-boundary individuals:

    $ xi={xmin+mod((xminxi),(xixmin))xmaxmod((xixmax),(xmaxxmin))
    $
    (16)

    In this section, the time complexity analysis of the proposed MCRO algorithm for the benchmark function optimization is analyzed. The time mainly depends on the searching loop. At first, the time complexity of the initialization method is $ O(NP\times D) $, where $ NP $ is the population size and $ D $ is the number of variables. Then for each iteration in the loop, the time is mainly spent in evaluating the fitness of each individual, which also costs $ O(NP \times D) $. Therefore, the time complexity of MCRO is $ O((I+1)(NP\times D)) $, where $ I $ is the number of iterations for the optimization process.

    To validate the effectiveness of the MCRO algorithm proposed in this paper, we first applied the algorithm to solve nine standard global optimization algorithms that are widely used. To ensure the fairness of the algorithm comparison, we did not modify any of the problems, which are tabulated in Table 1. As can be observed from the table, the first five problems are unimodal global optimization problems, where f05 is a noisy quadratic function. The remaining problems are all multimodal global optimization problems. It is worth noting that the local minima of these problems are increasing as the size of the problem increases, which is the reason why the difficulty of solving the problem increases with the size of the problem.

    Table 1.  Benchmark functions based in our experimental study.
    Test Function D Range optimum
    $ f_{01}=\sum_{i=1}^{n} x_{i}^{2} $ 30 [-100, 100] 0
    $ f_{02}=\sum_{i=1}^{n}\left|x_{i}\right|+\prod_{i=1}^{n}\left|x_{i}\right| $ 30 [-10, 10] 0
    $ f_{03}=\max _{i}\left\{x_{i} \mid, 1 \leq i \leq D\right\} $ 30 [-100, 100] 0
    $ f_{04}=\sum_{i=1}^{D}\left(\left\lfloor x_{i}+0.5\right\rfloor^{2}\right. $ 30 [-100, 100] 0
    $ f_{05}=\sum_{i=1}^{D} i x_{i}^{4}+random[0, 1) $ 30 [-1.28, 1.28] 0
    $ f_{06}=\frac{1}{400} \sum_{i-1}^{D} x_{i}^{2}-\prod_{i=1}^{D} \cos \left(\frac{x_{i}}{\sqrt{i}}\right)+1 $ 30 [-32, 32] 0
    $ f_{07}=-20 \exp \left(-0.2 \sqrt{\frac{1}{D} \sum_{i-1}^{D} x_{i}^{2}}\right)-\exp \left(\frac{1}{D} \sum_{i=1}^{D} \cos 2 \pi x_{i}\right)+20+e $ 30 [-600, 600] 0
    $ f_{08}=\frac{\pi}{D}\left\{10 \sin ^{2}\left(\pi y_{i}\right)+\sum_{i=1}^{D-1}\left(y_{i}-1\right)^{2}\left[1+10 \sin ^{2}\left(\pi y_{i}+1\right)\right]\right. $
    $ \left.+(y D-1)^{2}+\sum_{i=1}^{D} u\left(x_{i}, 10, 100, 4\right)\right\} $
    $ y_{i}=1+\frac{x_{i}+1}{4} \quad u\left(x_{i}, a, k, m\right)=\left\{k(xia)mxi>a0a<xi<ak(xia)mxi<a
    \right. $
    30 [-50, 50] 0
    $ f_{09}=0.1\left\{10 \sin ^{2}\left(\pi y_{i}\right)+\sum_{i=1}^{D-1}\left(y_{i}-1\right)^{2}\left[1+10 \sin ^{2}\left(\pi y_{i}+1\right)\right]+(y D-1)^{2}\right\} $
    $ +\sum_{i=1}^{D} u\left(x_{i}, 10, 100, 4\right) $
    30 [-50, 50] 0

     | Show Table
    DownLoad: CSV

    To illustrate the superiority of the MCRO algorithm proposed in our paper, we compared the MCRO algorithm with the standard CRO, the Differential Evolution algorithm (DE), the Adaptive Differential Evolution (jDE) [27] and the Modified Artificial Bee Colony algorithm (MABC) [28]. Across the experiments, the algorithms were computed with a maximum number of functions of 150, 000 and 50 independent experiments were conducted for each problem, and their average values were finally compared as the final results. To compare fairly, the parameter of different algorithm can be set as follows: for the CRO algorithm, the population size NP is 20, KELossRate is 0.1, Molecoll is 0.2, InitialKE is 1000, and the value of decThres and synThres is 150000 and 0, respectively. The value of MR is 0.2. For the DE algorithm, the population size is 100, F is 0.3 and CR is 0.7. The parameter of the jDE is the same with the paper[27]. For modified artificial bee colony algorithm (MABC), the parameters of this algorithm are identical to the paper [28]. The results are listed in Table 2. From Table 2, we can observe that the results of MCRO are better than the other algorithms while Figure 3 shows the 2D plot of f07.

    Table 2.  Performance comparisons of DE, jDE, MABC, CRO and MCRO.
    F DE jDE MABC CRO MCRO
    f01 6.1695e-14$ \pm $ 4.7200e-14 2.0141e-28$ \pm $ 4.1612e-28 1.9244e-39$ \pm $4.3030e-39 6.7348e-007$ \pm $2.1237e-007 5.5207e-71 $ \pm $1.0949e-70
    f02 4.2283e-07 $ \pm $3.0823e-07 1.6421e-17$ \pm $ 1.0644e-17 6.8088e-18$ \pm $ 1.5132e-17 0.00213.6637e-004 3.1288e-43$ \pm $3.9422e-43
    f03 0.0458 $ \pm $0.0809 0.8695$ \pm $ 0.7707 13.3649 $ \pm $ 3.6487 0.0097$ \pm $0.0035 2.7781e-10 $ \pm $ 2.2871e-10
    f04 0$ \pm $ 0 0$ \pm $0 0 $ \pm $ 0 0$ \pm $0 0$ \pm $ 0
    f05 0.0096$ \pm $0.0021 0.0063$ \pm $0.0014 0.0311 $ \pm $ 0.0083 0.0036$ \pm $0.0013 0.0083$ \pm $ 0.0023
    f06 1.3524e-13 $ \pm $7.7685e-14 0$ \pm $0 2.9016e-08$ \pm $ 6.4883e-08 4.7687e-07$ \pm $1.5030e-007 0$ \pm $0
    f07 7.8443e-08$ \pm $ 3.1063e-08 7.9936e-15$ \pm $0 2.8747e-06 $ \pm $ 6.4281e-06 0.0024$ \pm $4.2876e-004 1.6520e-14$ \pm $ 5.9448e-15
    f 08 3.5919e-15 $ \pm $2.2490e-15 4.6115e-30 $ \pm $5.3009e-30 2.3631e-07 $ \pm $5.2841e-07 0.1371$ \pm $0.3257 1.6222e-32 $ \pm $ 7.0698e-34
    f 09 4.7192e-14 $ \pm $3.9138e-14 5.5873e-29$ \pm $7.6744e-29 1.4035e-09 $ \pm $3.1382e-09 7.5513e-006$ \pm $5.3363e-006 7.0479e-31 $ \pm $ 7.6213e-31

     | Show Table
    DownLoad: CSV
    Figure 3.  2D Ackley function.

    To further demonstrate the effectiveness of the MCRO algorithm proposed in this paper, the proposed algorithm is applied to two different antenna array problems, including reconfigurable Antenna Array and linear antenna. {To enable a fair comparison of the algorithms}, the same number of function calculations is necessary. Therefore, in this paper, we use the same number of function evaluations for each algorithm as the final termination condition. As a result, we can evaluate the performance of the algorithms by the results of different algorithms and finally illustrate the efficiency of the proposed algorithm.

    In this section, we applied our proposed MCRO to address the reconfigurable dual-beam antenna array. In this problem, when the phase distribution of the array is suitably modified, the amplitude distribution can produce pencil-shaped or sector power patterns. In addition, to reduce the effect of coupling between elements, an additional term is incorporated in the objective function [29]. Therefore, we can optimize the amplitude-excitation dynamic range (ARD) to minimize the mutual coupling problem [30,31]. Then, the objective function to be optimized in this paper is defined as below.

    $ Ec(P)=3i=1(P(p)i,dP(i)i)2+4i=1(P(s)i,dP(s)i)2+ADR
    $
    (17)

    where ADR denotes the amplitude-dynamic ratio. Based on the above description, it can be noticed that the algorithm can motivate the difference between the amplitudes by optimizing the ADR. Meanwhile, we can reveal that the difference between ADR and amplitude is linear, and minimizing the ADR is accompanied by a smaller difference between amplitudes. Therefore, in this paper, we optimize the reconfigurable antenna-array by optimizing the effect of coupling.

    In order to efficiently design the above reconfigurable antenna array problem, we use our algorithm MCRO algorithm to tackle the problem. We compare the MCRO algorithm with some other algorithms including G3-GA, DE, jDE, MABC and MCRO. In the MCRO algorithm, the population size is 20 and the number of function computations is 20000, and the value of MR is 0.7. The rest of the parameters are the same as those of the previous global optimization algorithm. The experimental results are tabulated in Tables 3 and 4, mainly to determine the amplitude and phase excitation patterns for the dual-beam optimization. From these tables, we can see that the MCRO algorithm is superior to other algorithms, especially for reconfigurable antenna array designs with coupling effects. For the problems without coupling effects, we can find that most of these algorithms can find the optimal solutions. For problems with coupling effects, we can find that MCRO gives the best ADR and fitness values. These results show that our proposed algorithm is superior to other evolutionary algorithms. Figures 4 and 5 show the optimized excitation pattern and the dual-beam pattern, respectively. Figure 5 shows the design parameters that satisfy both pencil and sector beams. Figures 6 and 7 show the same case as Figures 4 and 5.

    Table 3.  The results of Experimental I and Experiment I with ADR.
    Element Number Experiment I Experiment I with ADR
    Amplitude Phase[deg.] Amplitude Phase[deg.]
    1/20 0.129861 -10.6276 0.149815 -0.80199
    2/19 0.174584 -25.4484 0.149813 -38.653
    3/18 0.255019 -36.7229 0.154374 -29.1362
    4/17 0.331274 -64.5232 0.256851 -71.1288
    5/16 0.367296 87.84077 0.307724 -72.5951
    6/15 0.556305 108.9342 0.412374 -97.438
    7/14 0.652548 -80.0141 0.483299 97.30298
    8/13 0.737711 -87.0023 0.564111 53.28282
    9/12 0.831561 -13.7764 0.616299 -106.146
    10/11 0.868279 43.79851 0.645976 88.67032
    ADR 6.68 4.31
    Fitness value 0.16 0.04

     | Show Table
    DownLoad: CSV
    Table 4.  Performance comparison of MCRO with G3-GA, DE, jDE, MABC and MCRO.
    Exp-I without ADR Exp-I with ADR
    fitness ADR fitness
    G3-GA 0.16 4.4137 0.1028
    DE 0.16 4.3470 0.04
    jDE 0.36 4.40 0.10
    MABC 0.16 4.31 0.06
    MCRO 0.16 4.31 0.04

     | Show Table
    DownLoad: CSV
    Figure 4.  Amplitude and phase excitation without coupling effect.
    Figure 5.  Dual-beam array patterns without coupling effect.
    Figure 6.  Dual-beam array patterns with coupling effect.
    Figure 7.  Dual-beam array patterns with coupling effect.

    In this section, we applied our proposed MCRO to address the linear antenna. To design of a linear antenna array [32], let us assume that we have 2N isotropic radiators placed symmetrically along x-axis. The objective function of side lobe suppression is:

    $ f1=i1Δϕiϕiiϕii|AF(ϕ)2|dϕ
    $
    (18)

    And for null control:

    $ f2=k|AF(ϕk)|2dϕ
    $
    (19)

    {where $ \Delta\phi_{i} $ is the bandwidth to be restrained} as $ \phi_{u i}-phi_{l i}, \phi_{k} $ is the null direction. To optimize both objectives simultaneously, we sum Eqs 18 and 19 and complete their sum to make the objective function that our algorithm needs to optimize.

    In order to elaborate the efficiency of the proposed MCRO algorithm, we came to test the reduction of the level of the side lobes in the linear array. Concurrently, we ran each algorithm independently for 25 times in each problem. The termination condition for each run of each algorithm was $ 5 \times 10^{4} $ number of function evaluations. All other parameters of the algorithms were consistent with the previous experiments. Meanwhile, the same number of function evaluations is applied to the other algorithms to construct a fair experiment.

    In this section, we design an array of 26 elements that has minimum SLL in bands [$ 0^{\circ} $, $ 82^{\circ} $] and [$ 98^{\circ} $, $ 180^{\circ} $] and null direction in $ 20^{\circ} $. Tables 5 and 6 summarize the results of the eight different evolutionary algorithms. From these tables, it can be viewed that our proposed MCRO algorithm is to be superior to other algorithms. Also, Figure 8 summarizes the array patterns of MCRO algorithm with other algorithms. From Figure 8, it can be seen that MCRO can suppress the generation of side lobes to the maximum extent. Moreover, it also obtains the minimum gain in the $ 20^{circ} $ null.

    Table 5.  Geometry of the 26 element linear array, normalized numbers with respect to $ \lambda/2 $ and null at $ 20^{o} $.
    MCRO 0.301 0.939 1.98 2.012 3.367 3.541 5.284 7.049 4.590 8.021 6.093 10.881 9.411
    MABC 0.330 0.979 1.821 2.039 3.267 2.897 4.151 6.574 7.749 4.880 5.799 9.544 10.909
    DE 0.472 1.049 2.090 2.058 2.897 3.799 4.825 6.318 5.373 7.475 8.61 10.231 11.558
    jDE 0.338 1.252 1.527 1.527 2.743 3.077 4.341 6.834 8.475 5.132 10.25 7.954 11.614

     | Show Table
    DownLoad: CSV
    Table 6.  Mean final objective function value, standard deviation, best, worst, median, and the Rank.
    Algorithm MCRO DE jDE MABC
    Mean 0.01194 0.02763 0.02356 0.0177
    Best 0.00921 0.01984 0.01736 0.0131
    Median 0.01219 0.02915 0.02319 0.0167
    Worst 0.01488 0.03396 0.02982 0.0224
    std 0.00209 0.00368 0.00443 0.0039
    Rank 1 4 3 2

     | Show Table
    DownLoad: CSV
    Figure 8.  13-element array for minimum SLL [0°, 82°] and [98°, 180°] and NULL 20°.

    In this paper, a new evolutionary search strategy is employed to enhance the performance of the algorithm. Moreover, we present a new population initialization by fusing stochastic initialization, opposition learning method, and generalized opposition learning method. It is worth noting that our proposed algorithm is simple in structure and easy to implement. To demonstrate the superiority of our proposed MCRO algorithm, we used nine benchmark functions as test functions. The results show that the proposed algorithm significantly outperforms the other algorithms. In addition, the algorithm was tested on the synthesis problem of two antenna arrays. We analyzed and evaluated the experimental results and compared them with other evolutionary optimization algorithms. In terms of the algorithm's performance, our algorithm is significantly better than the other algorithms. In the future, we would like to apply our proposed algorithm MCRO to solve other antenna arrays, such as circular antenna arrays. More quantifiable experiments on those arrays are expected. Besides, we will focus on exploring other strategies such as multiobjective mechanism to enhance the efficiency of the algorithm. Meanwhile, other evaluation criteria will be computed to estimate the performance of those extended algorithms.

    The work described in this paper was substantially supported by the National Natural Science Foundation of China under Grant No. 62076109.

    The authors declare there is no conflicts of interest.

    [1] Sowers KR, May HD (2013) In Situ Treatment of PCBs by Anaerobic Microbial Dechlorination in Aquatic Sediment: Are We There Yet? Curr Opin Biotechnol 24: 482-488. doi: 10.1016/j.copbio.2012.10.004
    [2] U.S. Environmental Protection Agency (2014) Polychlorinated Biphenyls (PCBs). Available at: http://www.epa.gov/epawaste/hazard/tsd/pcbs/.
    [3] World Health Organization Regional Office for Europe (1996) Air Quality Guidelines for Europe. 2nd Ed. WHO Regional Publications.
    [4] U.S. Department of Health and Human Services Agency for Toxic Substances and Disease Registry (2014) ATSDR Case Studies in Environmental Medicine Polychlorinated Biphenyls (PCBs) Toxicity.
    [5] Khairy MA, Weinstein MP, Lohmann R (2014) Trophodynamic Behavior of Hydrophobic Organic Contaminants in the Aquatic Food Web of a Tidal River. Environ Sci Technol 48: 12533-12542. doi: 10.1021/es502886n
    [6] Wakeman TH, Themelis NJ (2001) A Basin-Wide Approach to Dredged Material Management in New York/New Jersey Harbor. J Hazard Mater 85: 1-13. doi: 10.1016/S0304-3894(01)00218-7
    [7] Kinney CA, Furlong ET, Zaugg SD, et al. (2006) Survey of Organic Wastewater Contaminants in Biosolids Destined for Land Application. Environ Sci Technol 40: 7207-7215. doi: 10.1021/es0603406
    [8] Colborn T, Vom Saal FS, Soto AM (1993) Developmental Effects of Endocrine-Disrupting Chemicals in Wildlife and Humans. Environ Health Perspect 101(5): 378-384.
    [9] U.S. Environmental Protection Agency (2012) Superfund Site Information. Available at: http://www.epa.gov/superfund/sites/cursites/.
    [10] Payne RB, Fagervold SK, May HD, et al. (2013) Remediation of polychlorinated biphenyl impacted sediment by concurrent bioaugmentation with anaerobic halorespiring and aerobic degrading bactera. Environ Sci Technol 47(8): 3807-3815.
    [11] Johansen HR, Alexander J, Rossland OJ, et al. (1996) PCDDs, PCDFs, and PCBs in human blood in relation to consumption of crabs from a contaminated fjord area in Norway. Environ Health Perspect 104: 756-764. doi: 10.1289/ehp.96104756
    [12] Yao M, Li Z, Zhang X, et al. (2014) Polychlorinated biphenyls in the centralized wastewater Tteatment plant in a chemical industry zone: source, distribution, and removal. J Chem 2014.
    [13] Fuoco R, Giannarelli S, Onor M, et al. (2012) A snow/firn four-century record of polycyclic aromatic hydrocarbons (PAHs) and polychlorobiphenyls (PCBs) at Talos Dome (Antarctica). Microchem J 105: 133-141. doi: 10.1016/j.microc.2012.05.018
    [14] Fuoco R, Colombini M, Ceccarini A, et al. (1996) Polychlorobiphenyls in Antarctica. Microchem J 54(4): 384-390.
    [15] Kjellerup BV, Paul P, Ghosh U, et al. (2012) Spatial Distribution of PCB Dechlorinating Bacteria and Activities in Contaminated Soil. Appl Environ Soil Sci 2012.
    [16] Stoodley P, Sauer K, Davies DG, et al. (2002) Biofilms as Complex Differentiated Communities. Annu Rev Microbiol 56: 187-209. doi: 10.1146/annurev.micro.56.012302.160705
    [17] Leys D, Adrian L, Smidt H (2013) Organohalide respiration: microbes breathing chlorinated molecules. Philos Trans R Soc London Ser B, Biol Sci 368: 2012316.
    [18] Mercier A, Joulian C, Michel C, et al. (2014) Evaluation of Three Activated Carbons for Combined Adsorption and Biodegradation of PCBs in Aquatic Sediment. Water Res 59(0): 304-315.
    [19] Meagher RB (2000) Phytoremediation of toxic elemental and organic pollutants. Curr Opin Plant Biol 3: 153-162. doi: 10.1016/S1369-5266(99)00054-0
    [20] Mackova M, Macek T, Ocenaskova J, et al. (1997) Biodegradation of polychlorinated biphenyls by plant cells. Int Biodeterior Biodegradation 39(4): 317-325.
    [21] Passatore L, Rossetti S, Juwarkar AA, et al. (2014) Phytoremediation and bioremediation of polychlorinated biphenyls (PCBs): State of knowledge and research perspectives. J Hazard Mater 278: 189-202. doi: 10.1016/j.jhazmat.2014.05.051
    [22] Zanaroli G, Negroni A, Häggblom MM, et al. (2015) Microbial dehalogenation of organohalides in marine and estuarine environments. Curr Opin Biotechnol 33: 287-295. doi: 10.1016/j.copbio.2015.03.013
    [23] Cardon ZG, Gage DJ (2006) Resource Exchange in the Rhizosphere: Molecular Tools and the Microbial Perspective. Annu Rev Ecol Evol Syst 37(1): 459-488.
    [24] Van Aken B, Correa PA, Schnoor JL (2010) Phytoremediation of Polychlorinated Biphenyls: New Trends and Promises. Environ Sci Technol 44(8): 2767-2776.
    [25] Leigh MB, et al. (2006) Polychlorinated Biphenyl (PCB)-Degrading Bacteria Associated with Trees in a PCB-Contaminated Site. Appl Environ Microbiol 72(4): 2331-2342.
    [26] Erickson MD (1997) Analytical Chemistry of PCBs, 2 Ed. CRC Press; 90 p.
    [27] Pieper DH, Seeger M (2008) Bacterial metabolism of polychlorinated biphenyls. J Mol Microbiol Biotechnol 15: 121-138. doi: 10.1159/000121325
    [28] Cooper M, et al. (2015) Anaerobic Microbial Transformation of Halogenated Aromatics and Fate Prediction Using Electron Density Modelling. Environ Sci Technol Technol 49: 6018-6028. doi: 10.1021/acs.est.5b00303
    [29] Bedard DL (2008) A Case Study for Microbial Biodegradation: Anaerobic Bacterial Reductive Dechlorination of Polychlorinated Biphenyls-From Sediment to Defined Medium. Annu Rev Microbiol 62: 253-270. doi: 10.1146/annurev.micro.62.081307.162733
    [30] Vasilyeva GK, Strijakova ER (2007) Bioremediation of soils and sediments contaminated by polychlorinated biphenyls. Microbiology 76(6): 639-653.
    [31] Donlan RM (2002) Biofilms: microbial life on surfaces. Emerg Infect Dis 8(9): 881-890.
    [32] Borja J, Taleon DM, Auresenia J, et al. (2005) Polychlorinated biphenyls and their biodegradation. Process Biochem 40(6): 1999-2013.
    [33] Liu J, Schnoor JL (2008) Uptake and Translocation of Lesser-Chlorinated Polychlorinated Biphenyls (PCBs) in Whole Hybrid Poplar Plants after Hydroponic Exposure. Chemosphere 73(10): 1608-1616.
    [34] Pilon-Smits E (2005) Phytoremediation. Annu Rev Plant Biol 56: 15-39. doi: 10.1146/annurev.arplant.56.032604.144214
    [35] Rezek J, Macek T, Mackova M, et al. (2008) Hydroxy-PCBs, Methoxy-PCBs and Hydroxy-Methoxy-PCBs: Metabolites of Polychlorinated Biphenyls Formed In Vitro by Tobacco Cells. Environ Sci Technol 42(15): 5746-5751.
    [36] Moeckel C, Thomas GO, Barber JL, et al. (2008) Uptake and storage of PCBs by plant cuticles. Environ Sci Technol 42: 100-105. doi: 10.1021/es070764f
    [37] Wang W, Gorsuch JW, Hughes J (1997) Plants for Environmental Studies (CRC Press), p 484.
    [38] O'Connor GA, Kiehl D, Eiceman GA, et al. (1990) Plant uptake of sludge-borne PCBs. J Environ Qual 19: 113-118.
    [39] Mackova M, et al. (2009) Phyto/rhizoremediation studies using long-term PCB-contaminated soil. Environ Sci Pollut Res 16: 817-829. doi: 10.1007/s11356-009-0240-3
    [40] Campanella BF, Bock C, Schröder P (2002) Phytoremediation to increase the degradation of PCBs and PCDD/Fs. Environ Sci Pollut Res 9(1): 73-85.
    [41] Pivetz BE (2001) Phytoremediation of Contaminated Soil and Ground Water at Hazardous Waste Sites. U.S. Environmental Protection Agency, Office of Research and Development, Office of Solid Waste and Emergency Response.
    [42] Kuiper I, Lagendijk EL, Bloemberg GV, et al. (2004) Rhizoremediation: A Beneficial Plant-Microbe Interaction. Mol Plant-Microbe Interact 17(1): 6-15.
    [43] Hall-Stoodley L, Costerton JW, Stoodley P (2004) Bacterial biofilms: from the natural environment to infectious diseases. Nat Rev Microbiol 2(2): 95-108.
    [44] Macek T, Mackova M, Kas J (2000) Exploitation of plants for the removal of organics in environmental remediation. Biotechnol Adv 18: 23-34. doi: 10.1016/S0734-9750(99)00034-8
    [45] Liang Y, Meggo R, Hu D, et al. (2015) Microbial community analysis of switchgrass planted and unplanted soil microcosms displaying PCB dechlorination. Appl Microbiol Biotechnol 99: 6515-6526. doi: 10.1007/s00253-015-6545-x
    [46] Field J, Sierra-Alvarez R (2008) Microbial transformation and degradation of polychlorinated biphenyls. Environ Pollut 155: 1-12. doi: 10.1016/j.envpol.2007.10.016
    [47] Fagervold SK, Watts JEM, May HD, et al. (2005) Sequential Reductive Dechlorination of Meta-Chlorinated Polychlorinated Biphenyl Congeners in Sediment Microcosms by Two Different Chloroflexi Phylotypes. Appl Environ Microbiol 71(12): 8085-8090.
    [48] Vasilyeva GK, Strijakova ER, Shea PJ (2006) Proceedings of the NATO Advanced Research Workshop on Viable Methods of Soil and Water Pollution Monitoring, Protection and Remediation. Springer; p 309-322.
  • This article has been cited by:

    1. Huaijun Deng, Linna Liu, Jianyin Fang, Boyang Qu, Quanzhen Huang, A novel improved whale optimization algorithm for optimization problems with multi-strategy and hybrid algorithm, 2023, 205, 03784754, 794, 10.1016/j.matcom.2022.10.023
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(8453) PDF downloads(1355) Cited by(5)

Figures and Tables

Figures(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog