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.

Energy funneling and macromolecular conformational dynamics: a 2D Raman correlation study of PEG melting

Cascading energy landscapes through funneling has been postulated as a mechanistic route for achieving the lowest energy configuration of a macromolecular system (such as proteins and polymers). In particular, understanding the molecular mechanism for the melting and crystallization of polymers is a challenging fundamental question. The structural modifications that lead to the melting of poly(ethylene glycol) (PEG) are investigated here. Specific Raman bands corresponding to different configurations of the PEG chain have been identified, and the molecular structural dynamics of PEG melting have been addressed using a combination of Raman spectroscopy, 2D Raman correlation and density functional theory (DFT) calculations. The melting dynamics of PEG have been unambiguously explained along the C-O bond rotation coordinate.

Diffusion on a rugged energy landscape with spatial correlations

Rugged energy landscapes find wide applications in diverse fields ranging from astrophysics to protein folding. We study the dependence of diffusion coefficient (D) of a Brownian particle on the distribution width (epsilon) of randomness in a Gaussian random landscape by simulations and theoretical analysis. We first show that the elegant expression of Zwanzig Proc. Natl. Acad. Sci. U.S.A. 85, 2029 (1988)] for D(epsilon) can be reproduced exactly by using the Rosenfeld diffusion-entropy scaling relation. Our simulations show that Zwanzig's expression overestimates D in an uncorrelated Gaussian random lattice - differing by almost an order of magnitude at moderately high ruggedness. The disparity originates from the presence of ``three-site traps'' (TST) on the landscape - which are formed by the presence of deep minima flanked by high barriers on either side. Using mean first passage time formalism, we derive a general expression for the effective diffusion coefficient in the presence of TST, that quantitatively reproduces the simulation results and which reduces to Zwanzig's form only in the limit of infinite spatial correlation. We construct a continuous Gaussian field with inherent correlation to establish the effect of spatial correlation on random walk. The presence of TSTs at large ruggedness (epsilon >> k(B)T) gives rise to an apparent breakdown of ergodicity of the type often encountered in glassy liquids. (C) 2014 AIP Publishing LLC.

A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

Hierarchical sampling for metastable conformers determines biomolecular recognition: the case of malectin and diglucosylated N-glycan interactions

Structural information over the entire course of binding interactions based on the analyses of energy landscapes is described, which provides a framework to understand the events involved during biomolecular recognition. Conformational dynamics of malectin's exquisite selectivity for diglucosylated N-glycan (Dig-N-glycan), a highly flexible oligosaccharide comprising of numerous dihedral torsion angles, are described as an example. For this purpose, a novel approach based on hierarchical sampling for acquiring metastable molecular conformations constituting low-energy minima for understanding the structural features involved in a biologic recognition is proposed. For this purpose, four variants of principal component analysis were employed recursively in both Cartesian space and dihedral angles space that are characterized by free energy landscapes to select the most stable conformational substates. Subsequently, k-means clustering algorithm was implemented for geometric separation of the major native state to acquire a final ensemble of metastable conformers. A comparison of malectin complexes was then performed to characterize their conformational properties. Analyses of stereochemical metrics and other concerted binding events revealed surface complementarity, cooperative and bidentate hydrogen bonds, water-mediated hydrogen bonds, carbohydrate-aromatic interactions including CH-pi and stacking interactions involved in this recognition. Additionally, a striking structural transition from loop to beta-strands in malectin CRD upon specific binding to Dig-N-glycan is observed. The interplay of the above-mentioned binding events in malectin and Dig-N-glycan supports an extended conformational selection model as the underlying binding mechanism.

Majorana fermions in superconducting wires: Effects of long-range hopping, broken time-reversal symmetry, and potential landscapes

We present a comprehensive study of two of the most experimentally relevant extensions of Kitaev's spinless model of a one-dimensional p-wave superconductor: those involving (i) longer-range hopping and superconductivity and (ii) inhomogeneous potentials. We commence with a pedagogical review of the spinless model and, as a means of characterizing topological phases exhibited by the systems studied here, we introduce bulk topological invariants as well as those derived from an explicit consideration of boundary modes. In time-reversal symmetric systems, we find that the longer range hopping leads to topological phases characterized by multiple Majorana modes. In particular, we investigate a spin model that respects a duality and maps to a fermionic model with multiple Majorana modes; we highlight the connection between these topological phases and the broken symmetry phases in the original spin model. In the presence of time-reversal symmetry breaking terms, we show that the topological phase diagram is characterized by an extended gapless regime. For the case of inhomogeneous potentials, we explore phase diagrams of periodic, quasiperiodic, and disordered systems. We present a detailed mapping between normal state localization properties of such systems and the topological phases of the corresponding superconducting systems. This powerful tool allows us to leverage the analyses of Hofstadter's butterfly and the vast literature on Anderson localization to the question of Majorana modes in superconducting quasiperiodic and disordered systems, respectively. We briefly touch upon the synergistic effects that can be expected in cases where long-range hopping and disorder are both present.

On the measurement of the surface energy budget over a land surface during the summer monsoon

The measurement of surface energy balance over a land surface in an open area in Bangalore is reported. Measurements of all variables needed to calculate the surface energy balance on time scales longer than a week are made. Components of radiative fluxes are measured while sensible and latent heat fluxes are based on the bulk method using measurements made at two levels on a micrometeorological tower of 10 m height. The bulk flux formulation is verified by comparing its fluxes with direct fluxes using sonic anemometer data sampled at 10 Hz. Soil temperature is measured at 4 depths. Data have been continuously collected for over 6 months covering pre-monsoon and monsoon periods during the year 2006. The study first addresses the issue of getting the fluxes accurately. It is shown that water vapour measurements are the most crucial. A bias of 0.25% in relative humidity, which is well above the normal accuracy assumed the manufacturers but achievable in the field using a combination of laboratory calibration and field intercomparisons, results in about 20 W m(-2) change in the latent heat flux on the seasonal time scale. When seen on the seasonal time scale, the net longwave radiation is the largest energy loss term at the experimental site. The seasonal variation in the energy sink term is small compared to that in the energy source term.

Determination of diffusion parameters and activation energy of diffusion in V3Si phase with A15 crystal structure

Diffusion such is the integrated diffusion coefficient of the phase, the tracer diffusion coefficient of species at different temperatures and the activation energy for diffusion, are determined in V3Si phase with A15 crystal structure. The tracer diffusion coefficient of Si Was found to be negligible compared to the tracer diffusion coefficient of V. The calculated diffusion parameters will help to validate the theoretical analysis of defect structure of the phase, which plays an important role in the superconductivity.

Phase Equilibria in the System $Al_{2}O_{3}-CaO-CoO$ and Gibbs Energy of Formation of $Ca_{3}CoAl_{4}O_{10}$

The floating-zone method with different growth ambiences has been used to selectively obtain hexagonal or orthorhombic DyMnO3 single crystals. The crystals were characterized by x-ray powder diffraction of ground specimens and a structure refinement as well as electron diffraction. We report magnetic susceptibility, magnetization and specific heat studies of this multiferroic compound in both the hexagonal and the orthorhombic structure. The hexagonal DyMnO3 shows magnetic ordering of Mn3+ (S = 2) spins on a triangular Mn lattice at T-N(Mn) = 57 K characterized by a cusp in the specific heat. This transition is not apparent in the magnetic susceptibility due to the frustration on the Mn triangular lattice and the dominating paramagnetic susceptibility of the Dy3+ (S = 9/2) spins. At T-N(Dy) = 3 K, a partial antiferromagnetic order of Dy moments has been observed. In comparison, the magnetic data for orthorhombic DyMnO3 display three transitions. The data broadly agree with results from earlier neutron diffraction experiments, which allows for the following assignment: a transition from an incommensurate antiferromagnetic ordering of Mn3+ spins at T-N(Mn) = 39 K, a lock-in transition at Tlock-in = 16 K and a second antiferromagnetic transition at T-N(Dy) = 5 K due to the ordering of Dy moments. Both the hexagonal and the orthorhombic crystals show magnetic anisotropy and complex magnetic properties due to 4f-4f and 4f-3d couplings.

Effect of spatial dispersion of the dielectric on the binding energy of D- ion in Si and Ge

Using a multivalley effective mass theory, we obtain the binding energy of a D- ion in Si and Ge taking into account the spatial variation of the host dielectric function. We find that on comparison with experimental results the effect of spatial dispersion is important in the estimation of binding energy for the D- formed by As in Si and Ge. The effect is less significant for the case of D- formed by P and Sb donors.

Free Energy Gap Dependence of the Electron-Transfer Rate from the Inverted to the Normal Region

numerical study of the free energy gap (FEG) dependence of the electron-transfer rate in polar solvents is presented. This study is based on the generalized multidimensional hybrid model, which not only includes the solvent polarization and the molecular vibration modes, but also the biphasic polar response of the solvent. The free energy gap dependence is found to be sensitive to several factors, including the solvent relaxation rate, the electronic coupling between the surfaces, the frequency of the high-frequency quantum vibrational mode, and the magnitude of the solvent reorganization energy. It is shown that in some cases solvent relaxation can play an important role even in the Marcus normal regime. The minimal hybrid model involves a large number of parameters, giving rise to a diverse non-Marcus FEG behavior which is often determined collectively by these parameters. The model gives the linear free energy gap dependence of the logarithmic rate over a substantial range of FEG, spanning from the normal to the inverted regime. However, even for favorable values of the relevant parameters, a linear free energy gap dependence of the rate could be obtained only over a range of 5000-6000 cm(-1) (compared to the experimentally observed range of 10000 cm(-1) reported by Benniston et al.). The present work suggests several extensions/generalizations of the hybrid model which might be necessary to fully understand the observed free energy gap dependence.

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.

Long range resonance energy transfer from a dye molecule to graphene has (distance)-4 dependence

In our previous report on resonance energy transfer from a dye molecule to graphene [J. Chem. Phys.129, 054703 (2008)], we had derived an expression for the rate of energy transfer from a dye to graphene. An integral in the expression for the rate was evaluated approximately. We found a Yuwaka-type dependence of the rate on the distance. We now present an exact evaluation of the integral involved, leading to very interesting results. For short distances (z < 20 A), the present rate and the previous rate are in good agreement. For larger distances, the rate is found to have a z(-4) dependence on the distance, exactly. Thus we predict that for the case of pyrene on graphene, it is possible to observe fluorescence quenching up to a distance of 300 A. This is in sharp contrast to the traditional fluorescence resonance energy transfer where the quenching is observable only up to 100 A.

Constraining the low-energy pion electromagnetic form factor with space-like and phase of time-like data

The Taylor coefficients c and d of the EM form factor of the pion are constrained using analyticity, knowledge of the phase of the form factor in the time-like region, 4m(pi)(2) <= t <= t(in) and its value at one space-like point, using as input the (g - 2) of the muon. This is achieved using the technique of Lagrange multipliers, which gives a transparent expression for the corresponding bounds. We present a detailed study of the sensitivity of the bounds to the choice of time-like phase and errors present in the space-like data, taken from recent experiments. We find that our results constrain c stringently. We compare our results with those in the literature and find agreement with the chiral perturbation-theory results for c. We obtain d similar to O(10) GeV-6 when c is set to the chiral perturbation-theory values.

UVB irradiation alters the activities and kinetic properties of the enzymes of energy metabolism in rat lens during aging

Ultraviolet (UV) radiation is one of the major risk factors of cataract (loss of eye-lens transparency). The influence of UVB radiation (300 nm, 100 mu W cm(-2)) on the activity and apparent kinetic constants (K-m and V-max) of rat lens hexokinase (HK;EC2.7.1.1), phosphofructokinase (PFK;EC2.7.1.11), isocitrate dehydrogenase (ICDH;EC1.1.1.41) and malate dehydrogenase (MDH;EC1.1.1.37) of energy metabolism has been investigated by irradiating the lens homogenate of three-and 12-month-old rats. In the three-month-old group specific activities of HK and PFK are reduced by 56 and 43 %, respectively, and there is no change in ICDH and MDH activities after a 24 h exposure. On the other hand, in the 12-month-old group the decreases are 72, 71, 24 and 16 % for HK, PFK. ICDH and MDH, respectively. UVB irradiation increases the apparent K-m of HK and PFK (in both age groups), whereas the K-m of ICDH and MDH is not altered. While the decrease in V-max of these enzymes due to UVB exposure is only marginal in three-month-old rats, it is more pronounced (significant) in 12-month-old rats. A similar decrease in enzyme activities of HK and PFK is also observe upon UVB exposure of the intact rat lens. The photoinduced changes in energy metabolism may in turn have a bearing on lens transparency, particularly at an older age.