Monte Carlo simulation of Ising chains showing several phase transitions

A Monte Carlo simulation of Ising chains with competing short-range and infiniterange interactions has been carried out. Results show that whenever the system does not enter a metastable state, variation of temperature brings about phase transitions in the Ising chain. These phase transitions, except for two sets of interaction strengths, are generally of higher order and involve changes in the long-range order while the short-range order remains unaffected.

A dynamic renormalization group study of active nematics

We carry out a systematic construction of the coarse-grained dynamical equation of motion for the orientational order parameter for a two-dimensional active nematic, that is a nonequilibrium steady state with uniaxial, apolar orientational order. Using the dynamical renormalization group, we show that the leading nonlinearities in this equation are marginally irrelevant. We discover a special limit of parameters in which the equation of motion for the angle field bears a close relation to the 2d stochastic Burgers equation. We find nevertheless that, unlike for the Burgers problem, the nonlinearity is marginally irrelevant even in this special limit, as a result of a hidden fluctuation-dissipation relation. 2d active nematics therefore have quasi-long-range order, just like their equilibrium counterparts.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of Delta C-p of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Turbulence-induced melting of a nonequilibrium vortex crystal in a forced thin fluid film

To gain a better understanding of recent experiments on the turbulence-induced melting of a periodic array of vortices in a thin fluid film, we perform a direct numerical simulation of the two-dimensional Navier-Stokes equations forced such that, at low Reynolds numbers, the steady state of the film is a square lattice of vortices. We find that as we increase the Reynolds number, this lattice undergoes a series of nonequilibrium phase transitions, first to a crystal with a different reciprocal lattice and then to a sequence of crystals that oscillate in time. Initially, the temporal oscillations are periodic; this periodic behaviour becoming more and more complicated with increasing Reynolds number until the film enters a spatially disordered nonequilibrium statistical steady state that is turbulent. We study this sequence of transitions using fluid-dynamics measures, such as the Okubo-Weiss parameter that distinguishes between vortical and extensional regions in the flow, ideas from nonlinear dynamics, e.g. Poincare maps, and theoretical methods that have been developed to study the melting of an equilibrium crystal or the freezing of a liquid and that lead to a natural set of order parameters for the crystalline phases and spatial autocorrelation functions that characterize short- and long-range order in the turbulent and crystalline phases, respectively.

Multiple conformational states in crystals and in solution in alpha gamma hybrid peptides. Fragility of the C-12 helix in short sequences

The conformational properties of foldamers generated from alpha gamma hybrid peptide sequences have been probed in the model sequence Boc-Aib-Gpn-Aib-Gpn-NHMe. The choice of alpha-aminoisobutyryl (Aib) and gabapentin (Gpn) residues greatly restricts sterically accessible coil formational space. This model sequence was anticipated to be a short segment of the alpha gamma C-12 helix, stabilized by three successive 4 -> 1 hydrogen bonds, corresponding to a backbone-expanded analogue of the alpha polypeptide 3(10)-helix. Unexpectedly, three distinct crystalline polymorphs were characterized in the solid state by X-ray diffraction. In one form, two successive C-12 hydrogen bonds were obtained at the N-terminus, while a novel C-17 hydrogen-bonded gamma alpha gamma turn was observed at the C-terminus. In the other two polymorphs, isolated C-9 and C-7 hydrogen-bonded turns were observed at Gpn (2) and Gpn (4). Isolated C-12 and C-9 turns were also crystallographically established in the peptides Boc-Aib-Gpn-Aib-OMe and Boc-Gpn-Aib-NHMe, respectively. Selective line broadening of NH resonances and the observation of medium range NH(i)<-> NH(i+2) NOEs established the presence of conformational heterogeneity for the tetrapeptide in CDCl3 solution. The NMR results are consistent with the limited population of the continuous C-12 helix conformation. Lengthening of the (alpha gamma)(n) sequences in the nonapeptides Boc-Aib-Gpn-Aib-Gpn-Aib-Gpn-Aib-Gpn-Xxx (Xxx = Aib, Leu) resulted in the observation of all of the sequential NOEs characteristic of an alpha gamma C-12 helix. These results establish that conformational fragility is manifested in short hybrid alpha gamma sequences despite the choice of conformationally constrained residues, while stable helices are formed on chain extension.

Solution structure of O-glycosylated C-terminal leucine zipper domain of human salivary mucin (MUC7)

Solution structures of a 23 residue glycopeptide II (KIS* RFLLYMKNLLNRIIDDMVEQ, where * denotes the glycan Gal-beta-(1-3)-alpha-GalNAc) and its deglycosylated counterpart I derived from the C-terminal leucine zipper domain of low molecular weight human salivary mucin (MUC7) were studied using CD, NMR spectroscopy and molecular modeling. The peptide I was synthesized using the Fmoc chemistry following the conventional procedure and the glycopeptide II was synthesized incorporating the O-glycosylated building block (N alpha-Fmoc-Ser-[Ac-4,-beta-D-Gal-(1,3)-Ac(2)alpha-D-GalN(3)]-OPfp) at the appropriate position in stepwise assembly of peptide chain. Solution structures of these glycosylated and nonglycosylated peptides were studied in water and in the presence of 50% of an organic cosolvent, trifluoroethanol (TFE) using circular dichroism (CD), and in 50% TFE using two-dimensional proton nuclear magnetic resonance (2D H-1 NMR) spectroscopy. CD spectra in aqueous medium indicate that the apopeptide I adapts, mostly, a beta-sheet conformation whereas the glycopeptide II assumes helical structure. This transition in the secondary structure, upon glycosylation, demonstrates that the carbohydrate moiety exerts significant effect on the peptide backbone conformation. However, in 50% TFE both the peptides show pronounced helical structure. Sequential and medium range NOEs, C alpha H chemical shift perturbations, (3)J(NH:C alpha H) couplings and deuterium exchange rates of the amide proton resonances in water containing 50% TFE indicate that the peptide I adapts alpha-helical structure from Ile2-Val21 and the glycopeptide II adapts alpha-helical structure from Ser3-Glu22. The observation of continuous stretch of helix in both the peptides as observed by both NMR and CD spectroscopy strongly suggests that the C-terminal domain of MUC7 with heptad repeats of leucines or methionine residues may be stabilized by dimeric leucine zipper motif. The results reported herein may be invaluable in understanding the aggregation (or dimerization) of MUC7 glycoprotein which would eventually have implications in determining its structure-function relationship.

Solution conformation of a tetradecapeptide stabilized by two di-n-propyl glycine residues

The solution conformation of a designed tetradecapeptide Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (Dpg-14) containing two di-n-propyl glycine (Dpg) residues has been investigated by H-1 NMR and circular dichroism in organic solvents. The peptide aggregates formed at a concentration of 3 mm in the apolar solvent CDCl3 were broken by the addition of 12% v/v of the more polar solvent DMSO-d(6). Successive NiH <-> Ni+1H NOEs observed over the entire length of the sequence in this solvent mixture together with the observation of several characteristic medium-range NOEs support a major population of continuous helical conformations for Dpg-14. Majority of the observed coupling constants ((3)(alpha)(JNHC)(H)) also support phi values in the helical conformation. Circular dichroism spectra recorded in methanol and propan-2-ol give further support in favor of helical conformation for Dpg-14 and the stability of the helix at higher temperature. Copyright (C) 2010 European Peptide Society and John Wiley & Sons, Ltd.

Rapidly solidified Al–Cr alloys: Crystalline and quasicrystalline phases

Rapidly solidified Al–Cr alloys up to 20 at. % Cr were studied to delineate the extent of crystalline and quasicrystalline phase formation in these alloys in comparison with as-cast alloys by using transmission electron microscopy and x-ray diffraction technique. The icosahedral quasicrystals are observed from 7 to 15 at. % Cr alloys, while equilibrium ?–Al11Cr2 phase is completely absent. Both rapid solidification and subsequent thermal decomposition studies indicate that the main competing phase is ?–Al2Cr up to 15 at. % Cr. Beyond this composition ?–Al4Cr is the dominant phase together with a small amount of ?4–Al7Cr3. We have shown that the electron diffraction patterns of Al–Cr quasicrystals are often associated with a diffuse intensity distribution, indicative of short-range order. The change in quasilattice constant with composition suggests the existence of structural vacancies. Further, a sudden change from coarse to ultrafine quasicrystalline grain structure in Al-7 at. % Cr alloy points to a change in nucleation mechanism from heterogeneous to homogeneous mode during the rapid solidification.

Quantitative structural refinement of mn K edge XANES in LaMnO3 and CaMnO3 perovskites

The special magnetotransport properties of hole doped manganese perovskites originate from a complex interplay among structural, magnetic and electronic degree of freedom. In this picture the local atomic structure around Mn ions plays a special role and this is the reason why short range order techniques like X-ray absorption spectroscopy (XAS) have been deeply exploited for studying these compounds. The analysis of near edge region features (XANES) of XAS spectra can provide very fine details of the local structure around Mn, complementary to the EXAFS, so contributing to the full understanding of the peculiar physical properties of these materials. Nevertheless the XANES analysis is complicated by the large amount of structural and electronic details involved making difficult the quantitative interpretation.This work exploits the recently developed MXAN code to achieve a full structural refinement of the Mn K edge XANES of LaMnO3 and CaMnO3 compounds; they are the end compounds of the doped manganite series LaxCa1-xMnO3, in which the Mn ions are present only in one charge state as Mn3+ and Mn4+ respectively. The good agreement between the results derived from the analysis of near edge and extended region of the XAS spectra demonstrates that a quantitative picture of the local structure call be obtained from structural refinement of Mn K edge XANES data in these crystalline compounds. The XANES analysis offers, in addition.. the possibility to directly achieve information on the topology of local atomic structure around the absorber not directly achievable from EXAFS.

Local structure and magneto-transport in Sr2FeMoO6 oxides

Double perovskite oxides Sr2FeMoO6 have attracted a great interest for their peculiar magneto-transport properties, and, ill particular, for the large values of low-field magneto-resistance (MR) which remains elevated even at room temperature, thanks to their high Curie temperature (T-c > 400 K). These properties are strongly influenced by chemical cation disorder, that is by the relative arrangement of Fe and Mo on their sublattices: the regular alternation of Fe and Mo enhances the M R and saturation magnetization. On the contrary the disorder generally depresses the magnetization and worsen the MR response. In this work the X-ray absorption fine structure (XAFS) technique has been employed in order to probe the cation order from a local point of view. XAFS spectra were collected at the Fe and Mo K edges on Sr2FeMoO6 samples with different degree of long-range chemical order. The XAFS results prove that a high degree of short-range cation order is preserved, despite the different long-range order: the Fe-Mo correlations are always preferred over the Fe-Fe and Mo-Mo ones in the perfectly ordered as well as in highly disordered samples.

Ordering near the percolation threshold in models of two-dimensional interacting bosons with quenched dilution

Randomly diluted quantum boson and spin models in two dimensions combine the physics of classical percolation with the well-known dimensionality dependence of ordering in quantum lattice models. This combination is rather subtle for models that order in two dimensions but have no true order in one dimension, as the percolation cluster near threshold is a fractal of dimension between 1 and 2: two experimentally relevant examples are the O(2) quantum rotor and the Heisenberg antiferromagnet. We study two analytic descriptions of the O(2) quantum rotor near the percolation threshold. First a spin-wave expansion is shown to predict long-ranged order, but there are statistically rare points on the cluster that violate the standard assumptions of spin-wave theory. A real-space renormalization group (RSRG) approach is then used to understand how these rare points modify ordering of the O(2) rotor. A new class of fixed points of the RSRG equations for disordered one-dimensional bosons is identified and shown to support the existence of long-range order on the percolation backbone in two dimensions. These results are relevant to experiments on bosons in optical lattices and superconducting arrays, and also (qualitatively) for the diluted Heisenberg antiferromagnet La-2(Zn,Mg)(x)Cu1-xO4.

First-principles analysis of electron correlation, spin ordering and phonons in the normal state of FeSe1-x

We present first-principles density-functional-theory-based calculations to determine the effects of the strength of on-site electron correlation, magnetic ordering, pressure and Se vacancies on phonon frequencies and electronic structure of FeSe1-x. The theoretical equilibrium structure (lattice parameters) of FeSe depends sensitively on the value of the Hubbard parameter U of on-site correlation and magnetic ordering. Our results suggest that there is a competition between different antiferromagnetic states due to comparable magnetic exchange couplings between first- and second-neighbor Fe sites. As a result, a short range order of stripe antiferromagnetic type is shown to be relevant to the normal state of FeSe at low temperature. We show that there is a strong spin-phonon coupling in FeSe (comparable to its superconducting transition temperature) as reflected in large changes in the frequencies of certain phonons with different magnetic ordering, which is used to explain the observed hardening of a Raman-active phonon at temperatures (similar to 100 K) where magnetic ordering sets in. The symmetry of the stripe antiferromagnetic phase permits an induced stress with orthorhombic symmetry, leading to orthorhombic strain as a secondary order parameter at the temperature of magnetic ordering. The presence of Se vacancies in FeSe gives rise to a large peak in the density of states near the Fermi energy, which could enhance the superconducting transition temperature within the BCS-like picture.

Instability and dewetting of ultrathin solid viscoelastic films on homogeneous and heterogeneous substrates

Instability and dewetting engendered by the van der Waals force in soft thin (<100 nm) linear viscoelastic solid (e. g., elastomeric gel) films on uniform and patterned surfaces are explored. Linear stability analysis shows that, although the elasticity of the film controls the onset of instability and the corresponding critical wavelength, the dominant length-scale remains invariant with the elastic modulus of the film. The unstable modes are found to be long-wave, for which a nonlinear long-wave analysis and simulations are performed to uncover the dynamics and morphology of dewetting. The stored elastic energy slows down the temporal growth of instability significantly. The simulations also show that a thermodynamically stable film with zero-frequency elasticity can be made unstable in the presence of physico-chemical defects on the substrate and can follow an entirely different pathway with far fewer holes as compared to the viscous films. Further, the elastic restoring force can retard the growth of a depression adjacent to the hole-rim and thus suppress the formation of satellite holes bordering the primary holes. These findings are in contrast to the dewetting of viscoelastic liquid films where nonzero frequency elasticity accelerates the film rupture and promotes the secondary instabilities. Thus, the zero-frequency elasticity can play a major role in imposing a better-defined long-range order to the dewetted structures by arresting the secondary instabilities. (C) 2011 American Institute of Physics. doi: 10.1063/1.3554748]

Conjectures about phase turbulence in the complex Ginzburg-Landau equation

In the complex Ginzburg-Landau equation, we consider possible ''phase turbulent'' regimes, where asymptotic correlations are controlled by phase fluctuations rather than by topological defects. Conjecturing that the decay of such correlations is governed by the Kardar-Parisi-Zhang (KPZ) model of growing interfaces, we derive the following results: (1) A scaling ansatz implies that equal-time spatial correlations in 1d, 2d, and 3d decay like e(-Ax2 zeta), where A is a nonuniversal constant, and zeta=1/2 in 1d. (2) Temporal correlations decay as exp(-t(2 beta)h(t/L(z))), with the scaling law <(beta)over bar> = <(zeta)over bar>/z, where z = 3/2, 1.58..., and 1.66..., for d = 1,2, and 3 respectively. The scaling function h(y) approaches a constant as y --> 0, and behaves like y(2(beta-<(beta)over bar>)), for large y. If in 3d the associated KPZ model turns out to be in its weak-coupling (''smooth'') phase, then, instead of the above behavior, the CGLE exhibits rotating long-range order whose connected correlations decay like 1/x in space or 1/t(1/2) in time. (3) For system sizes, L, and times t respectively less than a crossover length, L(c), and time, t(c), correlations are governed by the free-field or Edwards-Wilkinson (EW) equation, rather than the KPZ model. In 1d, we find that L(c) is large: L(c) similar to 35,000; for L < L(c) we show numerical evidence for stretched exponential decay of temporal correlations with an exponent consistent with the EW value beta(EW)= 1/4.

Origin of a trend in ECMWF wind analysis and its objective removal

The suitability of the European Centre for Medium Range Weather Forecasting (ECMWF) operational wind analysis for the period 1980-1991 for studying interannual variability is examined. The changes in the model and the analysis procedure are shown to give rise to a systematic and significant trend in the large scale circulation features. A new method of removing the systematic errors at all levels is presented using multivariate EOF analysis. Objectively detrended analysis of the three-dimensional wind field agrees well with independent Florida State University (FSU) wind analysis at the surface. It is shown that the interannual variations in the detrended surface analysis agree well in amplitude as well as spatial patterns with those of the FSU analysis. Therefore, the detrended analyses at other levels as well are expected to be useful for studies of variability and predictability at interannual time scales. It is demonstrated that this trend in the wind field is due to the shift in the climatologies from the period 1980-1985 to the period 1986-1991.