Computational algorithms inspired by biological processes and evolution

In recent times computational algorithms inspired by biological processes and evolution are gaining much popularity for solving science and engineering problems. These algorithms are broadly classified into evolutionary computation and swarm intelligence algorithms, which are derived based on the analogy of natural evolution and biological activities. These include genetic algorithms, genetic programming, differential evolution, particle swarm optimization, ant colony optimization, artificial neural networks, etc. The algorithms being random-search techniques, use some heuristics to guide the search towards optimal solution and speed-up the convergence to obtain the global optimal solutions. The bio-inspired methods have several attractive features and advantages compared to conventional optimization solvers. They also facilitate the advantage of simulation and optimization environment simultaneously to solve hard-to-define (in simple expressions), real-world problems. These biologically inspired methods have provided novel ways of problem-solving for practical problems in traffic routing, networking, games, industry, robotics, economics, mechanical, chemical, electrical, civil, water resources and others fields. This article discusses the key features and development of bio-inspired computational algorithms, and their scope for application in science and engineering fields.

Evolution of electronic structure with dimensionality in divalent nickelates

We investigate the evolution of electronic structure with dimensionality (d) of Ni-O-Ni connectivity in divalent nickelates, NiO (3-d), La2NiO4, Pr2NiO4 (2-d), Y2BaNiO5 (1-d) and Lu2BaNi5 (0-d), by analyzing the valence band and the Ni 2p core-level photoemission spectra in conjunction with detailed many-body calculations including full multiplet interactions. Experimental results exhibit a reduction in the intensity of correlation-induced satellite features with decreasing dimensionality. The calculations based on the cluster model, but evaluating both Ni 3d and O 2p related photoemission processes on the same footing, provide a consistent description of both valence-band and core-level spectra in terms of various interaction strengths. While the correlation-induced satellite features in NiO is dominated by poorly screened d(8) states as described in the existing literature, we find that the satellite features in the nickelates with lower dimensional Ni-O-Ni connectivity are in fact dominated by the over-screened d(10)L(2) states. It is found that the changing electronic structure with the dimensionality is primarily driven by two factors: (i) a suppression of the nonlocal contribution to screening; and (ii) a systematic decrease of the charge-transfer energy Delta driven by changes in the Madelung potential. [S0163-1829(99)09619-8].

Evolution of sex ratios in social hymenoptera: consequences of finite brood size

Evolutionarily stable sex ratios are determined for social hymenoptera under local mate competition (LMC) and when the brood size is finite. LMC is modelled by the parameter d. Of the reproductive progeny from a single foundress nest, a fraction d disperses (outbreeding), while (1-d) mate amongst themselves (sibmating). When the brood size is finite, d is taken to be the probability of an offspring dispersing, and similarly, r, the proportion of male offspring, the probability of a haploid egg being laid. Under the joint influence of these two stochastic processes, there is a nonzero probability that some females remain unmated in the nest. As a result, the optimal proportion of males (corresponding to the evolutionarily stable strategy, ESS) is higher than that obtained when the brood size is infinite. When the queen controls the sex ration, the ESS becomes more female biased under increased inbreeding (lower d), However, the ESS under worker control shows an unexpected pattern, including an increase in the proportion of males with increased inbreeding. This effect is traced to the complex interaction between inbreeding and local mate competition.

The Evolution of Guanylyl Cyclases as Multidomain Proteins: Conserved Features of Kinase-Cyclase Domain Fusions

Guanylyl cyclases (GCs) are enzymes that generate cyclic GMP and regulate different physiologic and developmental processes in a number of organisms. GCs possess sequence similarity to class III adenylyl cyclases (ACs) and are present as either membrane-bound receptor GCs or cytosolic soluble GCs. We sought to determine the evolution of GCs using a large-scale bioinformatic analysis and found multiple lineage-specific expansions of GC genes in the genomes of many eukaryotes. Moreover, a few GC-like proteins were identified in prokaryotes, which come fused to a number of different domains, suggesting allosteric regulation of nucleotide cyclase activity Eukaryotic receptor GCs are associated with a kinase homology domain (KHD), and phylogenetic analysis of these proteins suggest coevolution of the KHD and the associated cyclase domain as well as a conservation of the sequence and the size of the linker region between the KHD and the associated cyclase domain. Finally, we also report the existence of mimiviral proteins that contain putative active kinase domains associated with a cyclase domain, which could suggest early evolution of the fusion of these two important domains involved in signa transduction.

Role of the triplet state in the green emission peak of polyfluorene films: A time evolution study

The blue emission of ethyl-hexyl substituted polyfluorene (PF2/6) films is accompanied by a low energy green emission peak around 500 nm in inert atmosphere. The intensity of this 500 nm peak is large in electroluminescence (EL) compared to photoluminescence (PL)measurements. Furthermore, the green emission intensity reduces dramatically in the presence of molecular oxygen. To understand this, we have modeled various nonradiative processes by time dependent quantum many body methods. These are (i) intersystem crossing to study conversion of excited singlets to triplets leading to a phosphorescence emission, (ii) electron-hole recombination (e-hR) process in the presence of a paramagnetic impurity to follow the yield of triplets in a polyene system doped with paramagnetic metal atom, and (iii) quenching of excited triplet states in the presence of oxygen molecules to understand the low intensity of EL emission in ambient atmosphere, when compared with that in nitrogen atmosphere. We have employed the Pariser-Parr-Pople Hamiltonian to model the molecules and have invoked electron-electron repulsions beyond zero differential approximation while treating interactions between the organic molecule and the rest of the system. Our time evolution methods show that there is a large cross section for triplet formation in the e-hR process in the presence of paramagnetic impurity with degenerate orbitals. The triplet yield through e-hR process far exceeds that in the intersystem crossing pathway, clearly pointing to the large intensity of the 500 nm peak in EL compared to PL measurements. We have also modeled the triplet quenching process by a paramagnetic oxygen molecule which shows a sizable quenching cross section especially for systems with large sizes. These studies show that the most probable origin of the experimentally observed low energy EL emission is the triplets.

Evolution of domain combinations in protein kinases and its implications for functional diversity

Protein kinases phosphorylating Ser/Thr/Tyr residues in several cellular proteins exert tight control over their biological functions. They constitute the largest protein family in most eukaryotic species. Protein kinases classified based on sequence similarity in their catalytic domains, cluster into subfamilies, which share gross functional properties. Many protein kinases are associated or tethered covalently to domains that serve as adapter or regulatory modules,naiding substrate recruitment, specificity, and also serve as scaffolds. Hence the modular organisation of the protein kinases serves as guidelines to their functional and molecular properties. Analysis of genomic repertoires of protein kinases in eukaryotes have revealed wide spectrum of domain organisation across various subfamilies of kinases. Occurrence of organism-specific novel domain combinations suggests functional diversity achieved by protein kinases in order to regulate variety of biological processes. In addition, domain architecture of protein kinases revealed existence of hybrid protein kinase subfamilies and their emerging roles in the signaling of eukaryotic organisms. In this review we discuss the repertoire of non-kinase domains tethered to multi-domain kinases in the metazoans. Similarities and differences in the domain architectures of protein kinases in these organisms indicate conserved and unique features that are critical to functional specialization. (C) 2009 Elsevier Ltd. All rights reserved.

The Evolution of Tribomaterial During Sliding: A Brief Introduction

This brief introductory article summarizes key findings from experiments and from computer simulations concerning the dramatic changes that commonly occur adjacent to sliding interfaces. We conclude that a wide range of observed features depends on a few basic processes (plastic deformation, interactions with the environment (including the counterface) and mechanical mixing) and that sliding leads to flow patterns similar to those expected in fluid flow.

Effect of Strain Rate on Evolution of the Deformation Microstructure and Texture in Polycrystalline Copper and Nickel

The evolution of crystallographic texture in polycrystalline copper and nickel has been studied. The deformation texture evolution in these two materials over seven orders of magnitude of strain rate from 3 x 10(-4) to similar to 2.0 x 10(+3) s(-1) show little dependence on the stacking fault energy (SFE) and the amount of deformation. Higher strain rate deformation in nickel leads to weakerh < 101 > texture because of extensive microband formation and grain fragmentation. This behavior, in turn, causes less plastic spin and hence retards texture evolution. Copper maintains the stable end < 101 > component over large strain rates (from 3 x 10(-4) to 10(+2) s(-1)) because of its higher strain-hardening rate that resists formation of deformation heterogeneities. At higher strain rates of the order of 2 x 10(+3) s(-1), the adiabatic temperature rise assists in continuous dynamic recrystallization that leads to an increase in the volume fraction of the < 101 > component. Thus, strain-hardening behavior plays a significant role in the texture evolution of face-centered cubic materials. In addition, factors governing the onset of restoration mechanisms like purity and melting point govern texture evolution at high strain rates. SFE may play a secondary role by governing the propensity of cross slip that in turn helps in the activation of restoration processes.

Texture evolution in commercially pure titanium after warm equal channel angular extrusion

Texture development in commercially pure titanium during equal channel angular extrusion (ECAE) through Routes A, Be and C has been studied up to three passes at 400 C. Textures were measured using X-ray diffraction, while the microstructural analyses were performed using electron back-scattered diffraction as well as transmission electron microscopy. Occurrences of dynamic restoration processes (recovery and recrystallization) were clearly noticed at all levels of deformations. Finally, the textures were simulated using a viscoplastic polycrystal self-consistent (VPSC) model. Simulations were performed incorporating basal, prismatic and pyramidal slip systems as well as tensile and compressive twinning. The simulated textures corroborate well with experimental textures in spite of the occurrence of dynamic restoration processes. (C) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.

Determination of solid fraction evolution during solidification of aluminum alloy in presence of linear electromagnetic stirring

In many industrial casting processes, knowledge of the solid fraction evolution during the solidification process is a key factor in determining the process parameters such as cooling rate, stirring intensity and in estimating the total solidification time. In the present work, a new method of estimating solid fraction is presented, which is based on calorimetric principles. In this method, the cooling curve data at each point in the melt, along with the thermal boundary conditions, are used to perform energy balance in the mould, from which solid fraction generation during any time interval can be estimated. This method is applied to the case of a rheocasting process, in which Al-Si alloy (A356 alloy) is solidified by stirring in a cylindrical mould placed in the annulus of a linear electromagnetic stirrer. The metal in the mould is simultaneously cooled and stirred to produce a cylindrical billet with non-dendritic globular microstructure. Temperature is measured at key locations in the mould to assess the various heat exchange processes prevalent in the mould and to monitor the solidification rate. The results obtained by energy balance method are compared with those by the conventional procedure of calculating solid fraction using the Schiel equation.

Genesis of return stroke current evolution at the wavefront

The channel dynamics at the wavefront is quite complex and is basically responsible for the evolution of return stroke current. The physical processes that actually contribute to the current evolution are not very clearly known. The enhancement of channel conductance at the wavefront is necessary for the current evolution and hence, return stroke. With regard to this, several questions arise like: (i) what causes the enhancement of this conductance, (ii) as the channel core temperature and electrical conductance are closely related, does one support the other and (iii) is the increase in core temperature on the nascent section of the channel is the result of free burning arc of the wavefront just below. These questions are investigated in detail in this work with appropriate transient thermal analysis and a macroscopic physical model for the lightning return stroke. Results clearly indicate that the contribution from the thermal field of the wavefront region to the adjacent nascent channel section is negligible as compared to the field enhancement brought in by the same. In other words, the whole process of return stroke evolution is dependent on the local heat generation at the nascent section caused by the enhancement of the electric field due to the arrival of the wavefront.

Microstructure and texture evolution during beta extrusion of boron modified Ti-6A1-4V alloy

Boron addition to conventional titanium alloys below the eutectic limit refines the cast microstructure and improves mechanical properties. The present work explores the influence of hypoeutectic boron addition on the microstructure and texture evolution in Ti-6Al-4V alloy under beta extrusion. The beta extruded microstructure of Ti-6Al-4V is characterized by shear bands parallel to the extrusion direction. In contrast, the extruded Ti-6Al-4V-0.1B alloy shows a regular beta worked microstructure consisting of fine prior beta grains and acicular alpha-lamellae with no signs of the microstructural instability. Crystallographic texture after extrusion was almost identical for the two alloys indicating the similarity in their transformation behavior, which is attributed to complete dynamic recrystallization during beta processing. Microstructural features as well as crystallographic texture indicate dominant grain boundary related deformation processes for the boron modified alloy that leads to homogeneous deformation without instability formation. The absence of shear bands has significant technological importance as far as the secondary processing of boron added alloys in (alpha + beta)-phase field are concerned. (C) 2012 Elsevier B.V. All rights reserved.

Evolution of microstructure and texture in friction stir processed Al-Mg-Mn alloy

Friction stir processing was carried out on the Al-Mg-Mn alloy to achieve ultrafine grained microstructure. The evolution of microstructure and micro-texture was studied in different regions of the deformed sample, namely nugget zone, thermo-mechanically affected zone (TMAZ) and base metal. The average grain sizes of the nugget zone, TMAZ and base metal are 1.5 mu m +/- 0.5 mu m, 15 mu m +/- 8 mu m, and 80 mu m +/- 10 mu m, respectively. The TMAZ exhibits excessive deformation banding structure and sub-grain formation. The orientation gradient within the sub-grain is dependent on grain size, orientation, and distance from nugget zone. The microstructure was partitioned based on the grain orientation spread and grain size values to separate the recrystallized fraction from the deformed region in order to understand the micromechanism of grain refinement. The texture of both deformed and recrystallized regions are similar in nature. Microstructure and texture analysis suggest that the restoration processes are different in different regions of the processed sample. The transition region between nugget zone and TMAZ exhibits large elongated grains surrounded by fine equiaxed grains of different orientation which indicate the process of discontinuous dynamic recrystallization. Within the nugget zone, similar texture between deformed and recrystallized grain fraction suggests that the restoration mechanism is a continuous process.

Chromium isotope variations (delta Cr-53/52) in mantle-derived sources and their weathering products: Implications for environmental studies and the evolution of delta Cr-53/52 in the Earth's mantle over geologic time

Here we report chromium isotope compositions, expressed as delta Cr-53/ 52 in per mil (&) relative to NIST 979, measured in selected Cr-rich minerals and rocks formed by the primary magmatic as well as the secondary metamorphic and weathering processes. The main objectives of this study were: (i) to further constrain the isotope composition of the Earth's mantle Cr inventory and its possible variation during geological history, based on the analysis of globally distributed and stratigraphically constrained mantle-derived chromites; and (ii) to investigate the magnitude and systematics of Cr isotope fractionation during oxidative weathering and secondary alteration (i. e., hydration, serpentinization) of the magmatic Cr sources. Specifically, we analyzed delta Cr-53/ 52 in a set of globally distributed mantle-derived chromites (FeMgCr2O4, n = 30) collected from various locations in Europe, Asia, Africa and South America, and our results confirm that a chromite-hosted Earth's mantle Cr inventory is uniform at - 0.079 +/- 0.129& (2SD), which we named here as a ` canonical' mantle d 53/ 52 Cr signature. Furthermore our dataset of stratigraphically constrained chromites, whose crystallization ages cover most of the Earth's geological history, indicate that the bulk Cr isotope composition of the chromite-hosted mantle inventory has remained uniform, within about +/- 0.100&, since at least the Early Archean times (similar to 3500 million years ago, Ma). To investigate the systematics of Cr isotope fractionation associated with alteration processes we analyzed a number of secondary Cr-rich minerals and variably altered ultramafic rocks (i. e., serpentinized harzburgites, lherzolites) that revealed large positive delta Cr-53/ 52 anomalies that are systematically shifted to higher values with an increasing degree of alteration and serpentinization. The degree of aqueous alteration and serpentinization was quantified by the abundances of fluid-mobile (Rb, K) elements, and by the Loss On Ignition (LOI) parameter, which determines the amount of structurally bound water (OH/ H2O) present in secondary hydrated minerals like serpentine. Overall, we observed that altered ultramafic rocks that yielded the highest LOI values, and the lowest amounts of fluid mobile elements, also yielded the heaviest delta Cr-53/ 52 signatures. Therefore, we conclude that secondary alteration (i.e., hydration, serpentinization) of ultramafic rocks in near-surface oxidative environments tend to shift the bulk Cr isotope composition of the weathered products to isotopically heavier values, pointing to a dynamic redox cycling of Cr in the Earth's crustal and near-surface environments. Hence, if validated by future