Free energy landscape of a dense hard-sphere system

The topography of the free energy landscape in phase space of a dense hard-sphere system characterized by a discretized free energy functional of the Ramakishnan-Yussouff form is investigated numerically using a specially devised Monte Carlo procedure. We locate a considerable number of glassy local minima of the free energy and analyze the distributions of the free energy at a minimum and an appropriately defined phase-space "distance" between different minima. We find evidence for the existence of pairs of closely related glassy minima("two-level systems"). We also investigate the way the system makes transitions as it moves from the basin of attraction of a minimum to that of another one after a start under nonequilibrium conditions. This allows us to determine the effective height of free energy barriers that separate a glassy minimum from the others. The dependence of the height of free energy barriers on the density is investigated in detail. The general appearance of the free energy landscape resembles that of a putting green: relatively deep minima separated by a fairly flat structure. We discuss the connection of our results with the Vogel-Fulcher law and relate our observations to other work on the glass transition.

Identity for harmonic oscillator brackets

An identity satisfied by the harmonic oscillator (Talmi-Moshinsky) brackets is derived from two equivalent methods for evaluating an integral often encountered in cluster model studies.

Generalization of Dynkin's Identity and Some Applications

Let X(t) be a right continuous temporally homogeneous Markov pro- cess, Tt the corresponding semigroup and A the weak infinitesimal genera- tor. Let g(t) be absolutely continuous and r a stopping time satisfying E.( S f I g(t) I dt) < oo and E.( f " I g'(t) I dt) < oo Then for f e 9iJ(A) with f(X(t)) right continuous the identity Exg(r)f(X(z)) - g(O)f(x) = E( 5 " g'(s)f(X(s)) ds) + E.( 5 " g(s)Af(X(s)) ds) is a simple generalization of Dynkin's identity (g(t) 1). With further restrictions on f and r the following identity is obtained as a corollary: Ex(f(X(z))) = f(x) + k! Ex~rkAkf(X(z))) + n-1E + (n ) )!.E,(so un-1Anf(X(u)) du). These identities are applied to processes with stationary independent increments to obtain a number of new and known results relating the moments of stopping times to the moments of the stopped processes.

Energy landscape view of phase transitions and slow dynamics in thermotropic liquid crystals

Thermotropic liquid crystals are known to display rich phase behavior on temperature variation. Although the nematic phase is orientationally ordered but translationally disordered, a smectic phase is characterized by the appearance of a partial translational order in addition to a further increase in orientational order. In an attempt to understand the interplay between orientational and translational order in the mesophases that thermotropic liquid crystals typically exhibit upon cooling from the high-temperature isotropic phase, we investigate the potential energy landscapes of a family of model liquid crystalline systems. The configurations of the system corresponding to the local potential energy minima, known as the inherent structures, are determined from computer simulations across the mesophases. We find that the depth of the potential energy minima explored by the system along an isochor grows through the nematic phase as temperature drops in contrast to its insensitivity to temperature in the isotropic and smectic phases. The onset of the growth of the orientational order in the parent phase is found to induce a translational order, resulting in a smectic-like layer in the underlying inherent structures; the inherent structures, surprisingly, never seem to sustain orientational order alone if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. The Arrhenius temperature dependence of the orientational relaxation time breaks down near the isotropic-nematic transition. We find that this breakdown occurs at a temperature below which the system explores increasingly deeper potential energy minima.

A matrix identity and perturbation expressions

An identity expressing formally the diagonal and off-diagonal elements of an inverse of a matrix is deduced employing operator techniques. Several well-known perturbation expressions for the self-energy are deduced as special cases. A new approximation and other applications following from the above formalism are briefly indicated through illustrations from a perturbed harmonic oscillator, chemisorption approximations and Kelly's result in the problem of electron correlation.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Both Surface Proteins (VP4 and VP7) of an Asymptomatic Neonatal Rotavirus Strain (1321) Have High Levels of Sequence Identity with the Homologous Proteins of a Serotype 10 Bovine Rotavirus

The nucleotide sequence of genes 4 and 9, encoding the outer capsid proteins VP4 and VP7 of a serotype 10 tissue culture-adapted strain, 1321, representative of asymptomatic neonatal rotaviruses isolated from neonates in Bangalore, India, were determined. Comparison of nucleotide and deduced amino acid sequences of 1321 VP4 and VP7 with previously published sequences of various serotypes revealed that both genes were highly homologous to the respective genes of serotype 10 bovine rotavirus, B223. The VP4 of 1321 represents a new human P serotype and the 1321 and related strains represent the first description of neonatal rotaviruses that appear to derive both surface proteins from an animal rotavirus.

A landscape approach to conservation of birds

Landscape ecology as a discipline in science is rather young. However its principles appear promising in outlining conservation strategies including a wide range of organisms, particularly birds. Birds due to their mobility use a variety of environmental resources, especially habitats. However, currently these habitats are only available in patches over most of the tropical world. Further whatever is left is under constant human pressure. This paper, therefore, addresses this problem and suggests means of dealing with it using the landscape approach as outlined by landscape ecology. The landscape approach starts with the realization that patches of habitats are open and interact with one another. Corridors of trees along roads, hedgerows and canals in a landscape can aid in the movement of species. Hence the landscape approach considers patches of habitats as interacting elements in the large matrix of the landscape. The landscape approach also integrates concepts. It puts together often debated issues such as whether to preserve maximum species diversity, to maximize representativeness, or to preserve only the valuable species. Based on a case study of the Uttara Kannada district in Karnataka, these oft-opposing views and complications can be dealt with practically and synthesized into a conservation strategy far the diverse avifauna of the Western Chats.

Assortative mating and the adaptive landscape

The ability of a population to shift from one adaptive peak to another was examined for a two-locus model with different degrees of assortative mating, selection, and linkage. As expected, if the proportion of the population that mates assortatively increases, so does its ability to shift to a new peak. Assortative mating affects this process by allowing the mean fitness of a population to increase monotonically as it passes through intermediate gene frequencies on the way to a new, higher, homozygotic peak. Similarly, if the height of the new peak increases or selection against intermediates becomes less severe, the population becomes more likely to shift to a new peak. Close linkage also helps the shift to a new adaptive peak and acts similarly to assortative mating, but it is not necessary for such a shift as was previously thought. When a population shifts to a new peak, the number of generations required is significantly less than that needed to return to the original peak when that happens. The short period of time required may be an explanation for rapid changes in the geological record. Under extremely high degrees of assortative mating, the shift takes longer, presumably because of the difficulty of breaking up less favored allelic combinations.

Synthon polymorphism and pseudopolymorphism in co-crystals. The 4,4 `-bipyridine-4-hydroxybenzoic acid structural landscape

Co-crystals of 4,4'-bipyridine and 4-hydroxybenzoic acid (1 : 2) show synthon polymorphism with the former being more stable. A 2 : 1 co-crystal is pseudopolymorphic within the same structural landscape with the structural roles of the two bipyridine N-atoms being distinct, as evidenced by mimicry by 4-phenylpyridine.

Polymorphs, Pseudopolymorphs, and Co-Crystals of Orcinol: Exploring the Structural Landscape with High Throughput Crystallography

An extensive search of the structural landscape of orcinol, 5-methyl-1,3-dihydroxybenzene, has been carried out with high throughput techniques. Polymorphs, pseudopolymorphs (solvates), and co-crystals are described. Several packing modes driven by O-H center dot center dot center dot N hydrogen bonds were identified for the orcinol N-base co-crystals and their hydrates. In these several structural variations, the OH group conformations in the orcinol molecule were found to depend on the choice of co-formers and the crystallization conditions employed. The structural landscape of a molecule is properly described by a sufficiently large number of related crystal structures, and high throughput crystallization followed by rapid structure determinations enables one to access these structures efficiently. Any understanding of this landscape would enable the crystal engineer to reasonably anticipate crystal structures of benzene-1,3-diol co-crystals with N-bases.

Levitation effect and its relationship with the underlying potential energy landscape

Positions of potential energy minima for spherical monatomic sorbates in zeolite NaY have been identified for different sizes of the sorbate. It is found that for small sorbates (sigma less than or equal to 4.96 Angstrom) there are only six adsorption sites per alpha-cage, which are located close to the inner surface of the alpha-cage. For larger sorbates, additional sites of comparable energies appear close to the 12-ring window which forms the bottleneck for intercage diffusion. Minimum energy paths between these sites have been computed. These suggest that the barriers for both intracage and intercage site-to-site migrations are comparable and decrease with increase in sorbate size. Earlier simulation studies on the diffusion of monatomic sorbates in zeolites indicated that there is a dramatic change in the nature of dependence of D on sorbate size around 4.96 Angstrom, for zeolite NaY. Therefore, the present results suggest that the dependence of D on sorbate size and the changes in the potential energy landscape are correlated. The sorbate-zeolite system is characterized by a flatter potential energy landscape when the sorbate size is large. (C) 1999 American Institute of Physics. [S0021-9606(99)51110-0].