Fast evolution strategies

Evolution strategies are a class of general optimisation algorithms which are applicable to functions that are multimodal, nondifferentiable, or even discontinuous. Although recombination operators have been introduced into evolution strategies, the primary search operator is still mutation. Classical evolution strategies rely on Gaussian mutations. A new mutation operator based on the Cauchy distribution is proposed in this paper. It is shown empirically that the new evolution strategy based on Cauchy mutation outperforms the classical evolution strategy on most of the 23 benchmark problems tested in this paper. The paper also shows empirically that changing the order of mutating the objective variables and mutating the strategy parameters does not alter the previous conclusion significantly, and that Cauchy mutations with different scaling parameters still outperform the Gaussian mutation with self-adaptation. However, the advantage of Cauchy mutations disappears when recombination is used in evolution strategies. It is argued that the search step size plays an important role in determining evolution strategies' performance. The large step size of recombination plays a similar role as Cauchy mutation.

Evolution of a mate recognition system after hybridization between two Drosophila species

I investigated the genetic relationship between male and female components of the mate recognition system and how this relationship influenced the subsequent evolution of the two traits, in a series of replicate populations of interspecific hybrids. Thirty populations of hybrids between Drosophila serrata and Drosophila birchii were established and maintained for 24 generations. At the fifth generation after hybridization, the mating success of hybrid individuals with the D. serrata parent was determined. The genetic correlation between male and female components of the male recognition system, as a consequence of pleiotropy or tight physical linkage, was found to be significant but low (r = 0.388). This result suggested that pleiotropy may play only a minor role in the evolution of mate recognition in this system. At the twenty-fourth generation after hybridization, the mating success of the hybrids was again determined. The evolution of male and female components was investigated by analyzing the direction of evolution of each hybrid line with respect to its initial position in relation to the genetic regression. Male and female components appeared to converge on a single equilibrium point, rather than evolving along trajectories with slope equal to the genetic regression, toward a line of equilibria.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Thermal and mechanical properties of radiation crosslinked poly(tetrafluoroethylene-co-perfluoromethyl vinyl ether)

The gamma-radiolysis of poly(tetrafluoroethylene-co-perfluoromethyl vinyl ether) (TFE/PMVE) was investigated using chemical and mechanical analyses. The polymer was found to form an insoluble network with a dose of gelation of 15.8 kGy. Tensile and glass transition temperature measurements indicated the predominance of crosslinking, with optimal elastomeric properties reached in the dose range of 120 to 200 kGy. Photoacoustic FTIR spectroscopy CPAS) showed the formation of new carboxylic acid end groups on irradiation. These new end groups were shown to decrease the thermal oxidative stability of the crosslinked network as determined by thermal gravimetric analysis. Electron spin resonance (ESR) studies of the polymer at 77 K indicated the presence of radical precursors. A G-value of 1.1 was determined for radical production at 77 K. Comparison of radical concentrations for a copolymer with a different mole ratio of PMVE, indicated that the PMVE units contribute to scission reactions. (C) 1998 Elsevier Science Ltd. All rights reserved.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Genetic structure and evolution of species in the mangrove genus Avicennia (Avicenniaceae) in the Indo-West Pacific

Allozyme variation in species of the mangrove genus Avicennia was screened in 25 populations collected from 22 locations in the Indo-West Pacific and eastern North America using 11 loci. Several fixed gene differences supported the specific status of Avicennia alba, A. integra, A. marina, and A. rumphiana from the Indo-West Pacific, and A. germinans from the Atlantic-East Pacific. The three varieties of A. marina, var. marina, var. eucalyptifolia, and var. australasica, had higher genetic similarities (Nei's I) and no fixed gene differences, confirming their conspecific status. Strong genetic structuring was observed in A. marina, with sharp changes in gene frequencies at the geographical margins of varietal distributions. The occurrence of alleles found otherwise in only one variety, in only immediately adjacent populations of another variety, provided evidence of introgession between varieties. The varieties appear to have diverged recently in the Pleistocene and are apparently not of ancient Cretaceous origin, as suggested earlier. Despite evidence of high degrees of outcrossing, gene flow among populations was relatively low (N(e)m less than or equal to 1-2), except where populations were geographically continuous, questioning assumptions that these widespread mangrove species achieve high levels of long-distance dispersal.

Crystal structure of human T cell leukemia virus type 1 gp21 ectodomain crystallized as a maltose-binding protein chimera reveals structural evolution of retroviral transmembrane proteins

Retroviral entry into cells depends on envelope glycoproteins, whereby receptor binding to the surface-exposed subunit triggers membrane fusion by the transmembrane protein (TM) subunit. We determined the crystal structure at 2.5-Angstrom resolution of the ectodomain of gp21, the TM from human T cell leukemia virus type 1. The gp21 fragment was crystallized as a maltose-binding protein chimera, and the maltose-binding protein domain was used to solve the initial phases by the method of molecular replacement. The structure of gp21 comprises an N-terminal trimeric coiled coil, an adjacent disulfide-bonded loop that stabilizes a chain reversal, and a C-terminal sequence structurally distinct from HIV type 1/simian immunodeficiency virus gp41 that packs against the coil in an extended antiparallel fashion. Comparison of the gp21 structure with the structures of other retroviral TMs contrasts the conserved nature of the coiled coil-forming region and adjacent disulfide-bonded loop with the variable nature of the C-terminal ectodomain segment. The structure points to these features having evolved to enable the dual roles of retroviral TMs: conserved fusion function and an ability to anchor diverse surface-exposed subunit structures to the virion envelope and infected cell surface. The structure of gp21 implies that the N-terminal fusion peptide is in close proximity to the C-terminal transmembrane domain and likely represents a postfusion conformation.

Some new examples for nonuniqueness of the evolution problem of harmonic maps

We find some new examples to show nonuniquence for the heat flow of harmonic maps where weak solutions satisfy the same monotonicity property.

A numerical study of pore-fluid, thermal and mass flow iu fluid-saturated porous rock basins

We present a numerical methodology for the study of convective pore-fluid, thermal and mass flow in fluid-saturated porous rock basins. lit particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore-fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore-fluid flow and mass transport in fluid-saturated hydrothermal basins; (2) Convective pore-fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence off the distribution patterns of convective pore;fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore-fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

Rapid evolution towards equal sex ratios in a system with heterogamety

The equal sex ratios found in many species with heterogametic sex determination may be a consequence of selection for equality or the result of the Mendelian segregation of the two sex chromosomes. A lack of genetic variation in sex ratio in species with heterogamety has been the major obstacle in distinguishing between these two hypotheses. We overcome this obstacle by generating hybrids between two species of Drosophila. The resulting hybrid lines had biased sex ratios, allowing us to observe the evolution of sex ratio in replicate populations. Sex ratio converged towards 1:1 after 16 generations of natural selection. These changes in sex ratio were not due to differences in viability between the sexes and the loci underlying the variation in sex ratio were not sex-linked. Equal sex ratios may therefore be the result of natural selection as Fisher predicted.

Thermal acclimation of locomotor performance in tadpoles of the frog Limnodynastes peronii

Previous analyses of thermal acclimation of locomotor performance in amphibians have only examined the adult life history stage and indicate that the locomotor system is unable to undergo acclimatory changes to temperature. In this study, we examined the ability of tadpoles of the striped marsh frog (Limnodynastes peronii) to acclimate their locomotor system by exposing them to either 10 degrees C or 24 degrees C for 6 weeks and testing their burst swimming performance at 10, 24, and 34 degrees C. At the test temperature of 10 degrees C, maximum velocity (U-max) of the 10 degrees C-acclimated tadpoles was 47% greater and maximum acceleration (A(max)) 53% greater than the 24 degrees C-acclimated animals. At 24 degrees C, U-max was 16% greater in the 10 degrees C-acclimation group, while there was no significant difference in A(max) or the time taken to reach U-max (T-U-max). At 34 degrees C, there was no difference between the acclimation groups in either U-max or A(max), however T-U-max was 36% faster in the 24 degrees C-acclimation group. This is the first study to report an amphibian (larva or adult) possessing the capacity to compensate for cool temperatures by thermal acclimation of locomotor performance. To determine whether acclimation period affected the magnitude of the acclimatory response, we also acclimated tadpoles of L. peronii to 10 degrees C for 8 months and compared their swimming performance with tadpoles acclimated to 10 degrees C for 6 weeks. At the test temperatures of 24 degrees C and 34 degrees C, U-max and A(max) were significantly slower in the tadpoles acclimated to 10 degrees C for 8 months. At 10 degrees C, T-U-max was 40% faster in the 8-month group, while there were no differences in either U-max or A(max). Although locomotor performance was enhanced at 10 degrees C by a longer acclimation period, this was at the expense of performance at higher temperatures.

Electrochemistry of macrocyclic cobalt(III/II) hexaamines: Electrocatalytic hydrogen evolution in aqueous solution

The macrocyclic cobalt hexaamines [Co(trans-diammac)](3+) and [Co(cis-diammac)](3+) (diammac = 6,13-dimethyl-1,4,8,11-tetraazacyclotetradecane-6,13-diamine) are capable of reducing the overpotential for hydrogen evolution on a mercury cathode in aqueous solution. Protons are reduced in a catalytic process involving reoxidation of the Co-II species to its parent Co-III complex. The cycle is robust at neutral pH with no decomposition of catalyst. The stability of the [Co(trans-diammac)](2+) and [Co(cis-diammac)](2+) complexes depends on the pH of the solution and the coordinating properties of the supporting electrolyte. Electrochemical studies indicate that the adsorbed Co-II complex on the surface of mercury is the active catalyst for the reduction of protons to dihydrogen.

Evolution of the genetic covariance between male and female components of mate recognition: An experimental test

The evolution of a positive genetic correlation between male and female components of mate recognition systems will result as a consequence of assortative mating and, in particular, is central to a number of theories of sexual selection. Although the existence of such genetic correlations has been investigated in a number of taxa, it has yet to be shown that such correlations evolve and whether they may evolve as rapidly as suggested by sexual selection models. In this study, I used a hybridization experiment to disrupt natural mate recognition systems and then observed the subsequent evolutionary dynamics of the genetic correlation between male and female components for 56 generations in hybrids between Drosophila serrata and Drosophila birchii. The genetic correlation between male and female components evolved from 0.388 at generation 5 to 1.017 at generation 37 and then declined to -0.040 after a further 19 generations. These results indicated that the genetic basis of the mate recognition system in the hybrid populations evolved rapidly. The initial rapid increase in the genetic correlation was consistent with the classic assumption that male and female components will coevolve under sexual selection. The subsequent decline in genetic correlation may be attributable to the fixation of major genes or, alternatively, may be a result of a cyclic evolutionary change in mate recognition.

Thermal and radiation curing of phenylethynyl terminated macromers

The thermal and gamma-irradiation induced curing of two phenylethynyl terminated composite resin systems, DFB/BPF and PETI5A, was investigated. Thermal curing of these matrix resin samples was performed at a temperature of 360 degrees C, gamma irradiation of the samples was conducted at 300 degrees C at a dose rate of 2.2 kGy h(-1). The reaction and subsequent loss of ethynyl groups in the resins for both cure methods was demonstrated by observing the decrease of the 2215 cm(-1) peak in the Raman spectra of the resins. Fully cured resin samples were found to have glass transition temperatures of 244-246 degrees C and 278-280 degrees C for DFB/BPF and PETI5A respectively. Similar relationships between T-g and fractional conversion were observed in both resins. The apparent polymerization rate, R-p, for thermal cure at 360 degrees C, was found to be 4.79 x 10(-2)% s(-1) in PETI5A and 3.22 x 10(-2)% s(-1) in DFB/BPF. Catastrophic degradation under nitrogen was observed to commence near 450 degrees C and 530 degrees C, with 5% weight losses occurring at 455 degrees C and 540 degrees C for DFB/BPF and PETI5A respectively. Gamma radiation induced cure at 300 degrees C was shown to be feasible, with full cure being reached with doses of 40 kGy for DFB/BPF and 100 kGy for PETI5A.