Variational principles in evolution

For a one-locus selection model, Svirezhev introduced an integral variational principle by defining a Lagrangian which remained stationary on the trajectory followed by the population undergoing selection. It is shown here (i) that this principle can be extended to multiple loci in some simple cases and (ii) that the Lagrangian is defined by a straightforward generalization of the one-locus case, but (iii) that in two-locus or more general models there is no straightforward extension of this principle if linkage and epistasis are present. The population trajectories can be constructed as trajectories of steepest ascent in a Riemannian metric space. A general method is formulated to find the metric tensor and the surface-in the metric space on which the trajectories, which characterize the variations in the gene structure of the population, lie. The local optimality principle holds good in such a space. In the special case when all possible linkage disequilibria are zero, the phase point of the n-locus genetic system moves on the surface of the product space of n higher dimensional unit spheres in a certain Riemannian metric space of gene frequencies so that the rate of change of mean fitness is maximum along the trajectory. In the two-locus case the corresponding surface is a hyper-torus.

First-principles DFT plus GW study of oxygen vacancies in rutile TiO2

We perform first-principles calculations of the quasiparticle defect states, charge transition levels, and formation energies of oxygen vacancies in rutile titanium dioxide. The calculations are done within the recently developed combined DFT + GW formalism, including the necessary electrostatic corrections for the supercells with charged defects. We find the oxygen vacancy to be a negative U defect, where U is the defect electron addition energy. For Fermi level values below similar to 2.8 eV (relative to the valence-band maximum), we find the +2 charge state of the vacancy to be the most stable, while above 2.8 eV we find that the neutral charge state is the most stable.

First principles study of magnetism in divalent Eu perovskites

A comparative first principles study has been carried out for EuLiH3 (ELH) and EuTiO3 (ETO) using the generalized gradient approximation +U approach. While ELH exhibits ferromagnetic ground state for all volumes, the magnetic ground state of ETO has the tendency to switch from G-type antiferromagnetic to a ferromagnetic state with change in volume. The marked difference in magnetic behavior and magnitude of the nearest neighbors exchange interaction of both the compounds are shown to be related to the difference in their respective electronic structure near the Fermi level. The Ti 3d states are shown to play predominant role in weakening the strength of the exchange interaction in ETO.

Cholinergic Muscarinic Receptors in Human Fetal Brain: Ontogeny of [3H]Quinuclidinyl Benzilate Binding Sites in Corpus Striatum, Brainstem, and Cerebellum

The ontogeny of muscarinic receptors was studied in human fetal striatum, brainstem, and cerebellum to investigate general principles of synaptogenesis as well as the physiological balance between various chemical synapses during development in a given region of the brain. [3H]Quinuclidinyl benzilate ([-'H]QNB) binding was assayed in total particulate fraction (TPF) from various parts of brain. In the corpus striatum, QNB binding sites are present at 16 weeks of gestation (average concentration 180 fmol/mg protein of TPF), slowly increase up to 24 weeks (average concentration 217 fmol/mg protein), and rapidly increase during the third trimester to 480 fmol/mg protein of TPF. In contrast, dopaminergic receptors exist as two subpopulations. one with low affinity and the other with high affinity up to the 24th week of gestation; all of them acquire the highaffinity characteristic during the third trimester. In brainstem, the muscarinic receptors show maximum concentration by 16 weeks of age (360 fmolimg protein of TPF). Subsequently the muscarinic receptor concentration shows a gradual decline in the brainstem. In cerebellum, except for a slight increase at 24 weeks (average concentration 90 fmol/mg protein of TPF), the receptor concentration remained nearly constant at about 60-70 fmolimg protein of TPF throughout fetal life. This study demonstrates that the ontogeny of muscarinic receptors varies among the different regions, and the patterns observed suggest that receptor formation occurs principally in the third trimester. Also noteworthy is the finding that the QNB binding sites decreased in all regions of the human brain during adult life. Key Words: Cholinergic muscarinic receptors-Quinuclidinyl benzilate-Corpus striaturn-Brainstem-Cerebellum. Ravikumar B. V. and Sastry P. S. Cholinergic muscarinic receptors in human fetal brain: Ontogeny of [3H]quinuclidinyl benzilate binding sites in corpus striatum, brainstem, and cerebellum. J. Neurochem. 45, 1948- 1950 (1985).

Adsorption isitherms for neutral organic compounds-A hierarchy in modelling: Part IV. Charge of maximum adsorption-Statistical mechanical basis

The charge at which adsorption of orgamc compounds attains a maximum ( \sigma MAX M) at an electrochenucal interface is analysed using several multi-state models in a hierarchical manner The analysis is based on statistical mechamcal results for the following models (A) two-state site parity, (B) two-state muhl-slte, and (C) three-state site parity The coulombic interactions due to permanent and reduced dipole effects (using mean field approximation), electrostatic field effects and specific substrate interactions have been taken into account. The simplest model in the hierarchy (two-state site parity) yields the exphcit dependence of ( \sigma MAX M) on the permanent dipole moment, polarizability of the solvent and the adsorbate, lattice spacing, effective coordination number, etc Other models in the baerarchy bring to hght the influence of the solvent structure and the role of substrate interactions, etc As a result of this approach, the "composition" of oM.x m terms of the fundamental molecular constants becomes clear. With a view to use these molecular results to maxamum advantage, the derived results for ( \sigma MAX M) have been converted into those involving experimentally observable parameters lake Co, C 1, E N, etc Wherever possible, some of the earlier phenomenologlcal relations reported for ( \sigma MAX M), notably by Parsons, Damaskm and Frumkln, and Trasattl, are shown to have a certain molecular basis, vlz a simple two-state sate panty model.As a corollary to the hxerarcbacal modelling, \sigma MAX M and the potential corresponding to at (Emax) are shown to be constants independent of 0max or Corg for all models The lmphcatlon of our analysis f o r OmMa x with respect to that predicted by the generalized surface layer equation (which postulates Om~ and Ema x varlaUon with 0) is discussed in detail Finally we discuss an passing o M. and the electrosorptlon valency an this context.

Gauge theory of a group of diffeomorphisms. I. General principles

Any (N+M)-parameter Lie group G with an N-parameter subgroup H can be realized as a global group of diffeomorphisms on an M-dimensional base space B, with representations in terms of transformation laws of fields on B belonging to linear representations of H. The gauged generalization of the global diffeomorphisms consists of general diffeomorphisms (or coordinate transformations) on a base space together with a local action of H on the fields. The particular applications of the scheme to space-time symmetries is discussed in terms of Lagrangians, field equations, currents, and source identities. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Acyclic Edge Coloring of Graphs with Maximum Degree 4

An acyclic edge coloring of a graph is a proper edge coloring such that there are no bichromatic cycles. The acyclic chromatic index of a graph is the minimum number k such that there is an acyclic edge coloring using k colors and is denoted by a'(G). It was conjectured by Alon, Sudakov, and Zaks that for any simple and finite graph G, a'(G) <= Delta+2, where Delta=Delta(G) denotes the maximum degree of G. We prove the conjecture for connected graphs with Delta(G)<= 4, with the additional restriction that m <= 2n-1, where n is the number of vertices and m is the number of edges in G. Note that for any graph G, m <= 2n, when Delta(G)<= 4. It follows that for any graph G if Delta(G)<= 4, then a'(G) <= 7.

On the origin of the inflectional instability of a laminar separation bubble

This is an experimental and theoretical Study of a laminar separation bubble and the associated linear stability mechanisms. Experiments were performed over a flat plate kept in a wind tunnel, with an imposed pressure gradient typical of an aerofoil that would involve a laminar separation bubble. The separation bubble was characterized by measurement of surface-pressure distribution and streamwise velocity using hot-wire anemometry. Single component hot-wire anemometry was also used for a detailed study of the transition dynamics. It was foundthat the so-called dead-air region in the front portion of the bubble corresponded to a region of small disturbance amplitudes, with the amplitude reaching a maximum value close to the reattachment point. An exponential growth rate of the disturbance was seen in the region upstream of the mean maximum height of the bubble, and this was indicative of a linear instability mechanism at work. An infinitesimal disturbance was impulsively introduced into the boundary layer upstream of separation location, and the wave packet was tracked (in an ensemble-averaged sense) while it was getting advected downstream. The disturbance was found to be convective in nature. Linear stability analyses (both the Orr-Sommerfeld and Rayleigh calculations) were performed for mean velocity profiles, starting from an attached adverse-pressure-gradient boundary layer all the way up to the front portion of the separation-bubble region (i.e. up to the end of the dead-air region in which linear evolution of the disturbance could be expected). The conclusion from the present work is that the primary instability mechanism in a separation bubble is inflectional in nature, and its origin can be traced back to upstream of the separation location. In other words, the inviscid inflectional instability of the separated shear layer should be logically seen as an extension of the instability of the upstream attached adverse-pressure-gradient boundary layer. This modifies the traditional view that pegs the origin of the instability in a separation bubble to the detached shear layer Outside the bubble, with its associated Kelvin-Helmholtz mechanism. We contendthat only when the separated shear layer has moved considerably away from the wall (and this happens near the maximum-height location of the mean bubble), a description by the Kelvin-Helmholtz instability paradigm, with its associated scaling principles, Could become relevant. We also propose a new scaling for the most amplified frequency for a wall-bounded shear layer in terms of the inflection-point height and the vorticity thickness and show it to be universal.

Theoretical Maximum Limits on Power-Handling Capacity of HVDC Cables

In the case of an ac cable, power transmission is limited by the length of the cable due to the capacitive reactive current component. It is well known that high-voltage direct current (HVDC) cables do not have such limitations. However, insulation-related thermal problems pose a limitation on the power capability of HVDC cables. The author presents a viable theoretical development, a logical extension to Whitehead's theory on thermal limitations of the insulation. The computation of the maximum power-carrying capability of HVDC cables subject to limits on the maximum operable temperature of the insulation is presented. The limitation on the power-carrying capability is closely associated with the electrothermal insulation failure. The effect of environmental interaction by way of external thermal resistance, an important aspect, is also considered in the formulations. The Lagrange multiplier method has been used to handle the ensuing optimization problem. The theory is based on an accepted theory of thermal breakdown in insulation and is an important and a coherent extension of great significance.

On the axiomatic approach to the maximum entropy principle of inference

Recent axiomatic derivations of the maximum entropy principle from consistency conditions are critically examined. We show that proper application of consistency conditions alone allows a wider class of functionals, essentially of the form ∝ dx p(x)[p(x)/g(x)] s , for some real numbers, to be used for inductive inference and the commonly used form − ∝ dx p(x)ln[p(x)/g(x)] is only a particular case. The role of the prior densityg(x) is clarified. It is possible to regard it as a geometric factor, describing the coordinate system used and it does not represent information of the same kind as obtained by measurements on the system in the form of expectation values.

Approximation algorithm for maximum independent set in planar traingle-free graphs

The maximum independent set problem is NP-complete even when restricted to planar graphs, cubic planar graphs or triangle free graphs. The problem of finding an absolute approximation still remains NP-complete. Various polynomial time approximation algorithms, that guarantee a fixed worst case ratio between the independent set size obtained to the maximum independent set size, in planar graphs have been proposed. We present in this paper a simple and efficient, O(|V|) algorithm that guarantees a ratio 1/2, for planar triangle free graphs. The algorithm differs completely from other approaches, in that, it collects groups of independent vertices at a time. Certain bounds we obtain in this paper relate to some interesting questions in the theory of extremal graphs.

On the Maximum Cloud Zone and the ITCZ over Indian, Longitudes during the Southwest Monsoon

An investigation is presented of the daily variation of the maximum cloud zone (MCZ) and the 7W mb trough in the Northern Hemisphere over the Indian longitudes 70–90°E during April–October for 1973–77. It is found that during June–September there are two favorable locations for a MCZ over these longitudes–on a majority of days the MCZ is present in the monsoon zone north of 15°N, and often a secondary MCZ occurs in the equatorial region (0–10°N). The monsoon MCZ gets established by northward movement of the MCZ occurring over the equatorial Indian ocean in April and May. The secondary MCZ appears intermittently, and is characterized by long spells of persistence only when the monsoon MCZ is absent. In each of the seasons studied, the MCZ temporarily disappeared from the mean summer monsoon location (15–28°N) about four weeks after it was established near the beginning of July. It is reestablished by the northward movement of the secondary MCZ, which becomes active during the absence of the monsoon MCZ, in a manner strikingly similar to that observed in the spring to summer transition. A break in monsoon conditions prevails just prior to the temporary disappearance of the monsoon MCZ. Thus we conclude that the monsoon MCZ cannot survive for longer than a month without reestablishment by the secondary MCZ. Possible underlying mechanisms are also discussed.

Optimal simultaneous maximum a posterioriestimation of states, noise statistics and parameters I. Algorithm

The simultaneous state and parameter estimation problem for a linear discrete-time system with unknown noise statistics is treated as a large-scale optimization problem. The a posterioriprobability density function is maximized directly with respect to the states and parameters subject to the constraint of the system dynamics. The resulting optimization problem is too large for any of the standard non-linear programming techniques and hence an hierarchical optimization approach is proposed. It turns out that the states can be computed at the first levelfor given noise and system parameters. These, in turn, are to be modified at the second level.The states are to be computed from a large system of linear equations and two solution methods are considered for solving these equations, limiting the horizon to a suitable length. The resulting algorithm is a filter-smoother, suitable for off-line as well as on-line state estimation for given noise and system parameters. The second level problem is split up into two, one for modifying the noise statistics and the other for modifying the system parameters. An adaptive relaxation technique is proposed for modifying the noise statistics and a modified Gauss-Newton technique is used to adjust the system parameters.

Optimal simultaneous maximum a posterioriestimation of states, noise statistics and parameters II. Numerical performance

A very general and numerically quite robust algorithm has been proposed by Sastry and Gauvrit (1980) for system identification. The present paper takes it up and examines its performance on a real test example. The example considered is the lateral dynamics of an aircraft. This is used as a vehicle for demonstrating the performance of various aspects of the algorithm in several possible modes.

Study of Structural Stability and Electronic Structure of Nonstoichiometric CdS Nano Clusters from First Principles

Structural stability of small sized nonstoichiometric CdS nano clusters between zincblende and wurtzite structures has been investigated using first-principles density functional calculations. Our study shows that the relative stability of these two structures depends sensitively on whether the surface is S-terminated or Cd-terminated. The associated band gap also exhibits non-monotonic behavior as a function of cluster size. Our findings may shed light on contradictory reports of experimentally observed structures of CdS nano clusters found in the literature.