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.

Shape and chemically anisotropic colloidal particles: A theoretical study of their structure, phase behavior, and slow dynamics

A detailed non-equilibrium state diagram of shape-anisotropic particle fluids is constructed. The effects of particle shape are explored using Naive Mode Coupling Theory (NMCT), and a single particle Non-linear Langevin Equation (NLE) theory. The dynamical behavior of non-ergodic fluids are discussed. We employ a rotationally frozen approach to NMCT in order to determine a transition to center of mass (translational) localization. Both ideal and kinetic glass transitions are found to be highly shape dependent, and uniformly increase with particle dimensionality. The glass transition volume fraction of quasi 1- and 2- dimensional particles fall monotonically with the number of sites (aspect ratio), while 3-dimensional particles display a non-monotonic dependence of glassy vitrification on the number of sites. Introducing interparticle attractions results in a far more complex state diagram. The ideal non-ergodic boundary shows a glass-fluid-gel re-entrance previously predicted for spherical particle fluids. The non-ergodic region of the state diagram presents qualitatively different dynamics in different regimes. They are qualified by the different behaviors of the NLE dynamic free energy. The caging dominated, repulsive glass regime is characterized by long localization lengths and barrier locations, dictated by repulsive hard core interactions, while the bonding dominated gel region has short localization lengths (commensurate with the attraction range), and barrier locations. There exists a small region of the state diagram which is qualified by both glassy and gel localization lengths in the dynamic free energy. A much larger (high volume fraction, and high attraction strength) region of phase space is characterized by short gel-like localization lengths, and long barrier locations. The region is called the attractive glass and represents a 2-step relaxation process whereby a particle first breaks attractive physical bonds, and then escapes its topological cage. The dynamic fragility of fluids are highly particle shape dependent. It increases with particle dimensionality and falls with aspect ratio for quasi 1- and 2- dimentional particles. An ultralocal limit analysis of the NLE theory predicts universalities in the behavior of relaxation times, and elastic moduli. The equlibrium phase diagram of chemically anisotropic Janus spheres and Janus rods are calculated employing a mean field Random Phase Approximation. The calculations for Janus rods are corroborated by the full liquid state Reference Interaction Site Model theory. The Janus particles consist of attractive and repulsive regions. Both rods and spheres display rich phase behavior. The phase diagrams of these systems display fluid, macrophase separated, attraction driven microphase separated, repulsion driven microphase separated and crystalline regimes. Macrophase separation is predicted in highly attractive low volume fraction systems. Attraction driven microphase separation is charaterized by long length scale divergences, where the ordering length scale determines the microphase ordered structures. The ordering length scale of repulsion driven microphase separation is determined by the repulsive range. At the high volume fractions, particles forgo the enthalpic considerations of attractions and repulsions to satisfy hard core constraints and maximize vibrational entropy. This results in site length scale ordering in rods, and the sphere length scale ordering in Janus spheres, i.e., crystallization. A change in the Janus balance of both rods and spheres results in quantitative changes in spinodal temperatures and the position of phase boundaries. However, a change in the block sequence of Janus rods causes qualitative changes in the type of microphase ordered state, and induces prominent features (such as the Lifshitz point) in the phase diagrams of these systems. A detailed study of the number of nearest neighbors in Janus rod systems reflect a deep connection between this local measure of structure, and the structure factor which represents the most global measure of order.

Actin Dynamics in Muscle Cells

In every cell, actin is a key component involved in migration, cytokinesis, endocytosis and generation of contraction. In non-muscle cells, actin filaments are very dynamic and regulated by an array of proteins that interact with actin filaments and/or monomeric actin. Interestingly, in non-muscle cells the barbed ends of the filaments are the predominant assembly place, whereas in muscle cells actin dynamics was reported to predominate at the pointed ends of thin filaments. The actin-based thin filament pointed (slow growing) ends extend towards the middle of the sarcomere's M-line where they interact with the thick filaments to generate contraction. The actin filaments in muscle cells are organized into a nearly crystalline array and are believed to be significantly less dynamic than the ones in other cell types. However, the exact mechanisms of the sarcomere assembly and turnover are largely unknown. Interestingly, although sarcomeric actin structures are believed to be relatively non-dynamic, many proteins promoting actin dynamics are expressed also in muscle cells (e.g ADF/cofilin, cyclase-associated protein and twinfilin). Thus, it is possible that the muscle-specific isoforms of these proteins promote actin dynamics differently from their non-muscle counterparts, or that actin filaments in muscle cells are more dynamic than previously thought. To study protein dynamics in live muscle cells, I used primary cell cultures of rat cardiomyocytes. My studies revealed that a subset of actin filaments in cardiomyocyte sarcomeres displays rapid turnover. Importantly, I discovered that the turnover of actin filaments depends on contractility of the cardiomyocytes and that the contractility-induced actin dynamics plays an important role in sarcomere maturation. Together with previous studies those findings suggest that sarcomeres undergo two types of actin dynamics: (1) contractility-dependent turnover of whole filaments and (2) regulatory pointed end monomer exchange to maintain correct thin filament length. Studies involving an actin polymerization inhibitor suggest that the dynamic actin filament pool identified here is composed of filaments that do not contribute to contractility. Additionally, I provided evidence that ADF/cofilins, together with myosin-induced contractility, are required to disassemble non-productive filaments in developing cardiomyocytes. In addition, during these studies we learned that isoforms of actin monomer binding protein twinfilin, Twf-1 and Twf-2a localise to myofibrils in cardiomyocytes and may thus contribute to actin dynamics in myofibrils. Finally, in collaboration with Roberto Dominguez s laboratory we characterized a new actin nucleator in muscle cells - leiomodin (Lmod). Lmod localises towards actin filament pointed ends and its depletion by siRNA leads to severe sarcomere abnormalities in cardiomyocytes. The actin filament nucleation activity of Lmod is enhanced by interactions with tropomyosin. We also revealed that Lmod expression correlates with the maturation of myofibrils, and that it associates with sarcomeres only at relatively late stages of myofibrillogenesis. Thus, Lmod is unlikely to play an important role in myofibril formation, but rather might be involved in the second step of the filament arrangement and/or maintenance through its ability to promote tropomyosin-induced actin filament nucleation occurring at the filament pointed ends. The results of these studies provide valuable new information about the molecular mechanisms underlying muscle sarcomere assembly and turnover. These data offer important clues to understanding certain physiological and pathological behaviours of muscle cells. Better understanding of the processes occurring in muscles might help to find strategies for determining, diagnosis, prognosis and therapy in heart and skeletal muscles diseases.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Application of static charge transfer within an ionic-liquid force field and its effect on structure and dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3-dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are +/- 0.5 e. While the three-dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly-with the positive result that for ionic charges of +/- 0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial-charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic-liquid models in which the ionic charges are assumed to be +/- 1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled-charge sets may be of benefit in the simulation of these systems-especially when looking at properties beyond liquid structure-thus providing on alternative to computationally expensive polarisable force fields.

Simplification of structured population dynamics in complex ecological communities

Goldstone's idea of slow dynamics resulting from spontaneously broken symmetries is applied to Hubbell's neutral hypothesis of community dynamics, to efficiently simplify stage-structured multi-species models-introducing the quasi-neutral approximation (QNA). Rather than assuming population-dynamical neutrality in the QNA, deviations from ideal neutrality, thought to be small, drive dynamics. The QNA is systematically derived to first and second order in a two-scale singular perturbation expansion. The total reproductive value of species, as computed from the effective life-history parameters resulting from the non-linear interactions with the surrounding community, emerges as the new dynamic variables in this aggregated description. Using a simple stage-structured community-assembly model, the QNA is demonstrated to accurately reproduce population dynamics in large, complex communities. Further, the utility of the QNA in building intuition for management problems is illustrated by estimating the responses of a fish stock to harvesting and variations in fecundity.

Analysis of temperature time-series: embedding dynamics into the MDS method

Global warming and the associated climate changes are being the subject of intensive research due to their major impact on social, economic and health aspects of the human life. Surface temperature time-series characterise Earth as a slow dynamics spatiotemporal system, evidencing long memory behaviour, typical of fractional order systems. Such phenomena are difficult to model and analyse, demanding for alternative approaches. This paper studies the complex correlations between global temperature time-series using the Multidimensional scaling (MDS) approach. MDS provides a graphical representation of the pattern of climatic similarities between regions around the globe. The similarities are quantified through two mathematical indices that correlate the monthly average temperatures observed in meteorological stations, over a given period of time. Furthermore, time dynamics is analysed by performing the MDS analysis over slices sampling the time series. MDS generates maps describing the stations’ locus in the perspective that, if they are perceived to be similar to each other, then they are placed on the map forming clusters. We show that MDS provides an intuitive and useful visual representation of the complex relationships that are present among temperature time-series, which are not perceived on traditional geographic maps. Moreover, MDS avoids sensitivity to the irregular distribution density of the meteorological stations.

Comments on 'Balance and the slow quasimanifold: some explicit results'

The concept of slow vortical dynamics and its role in theoretical understanding is central to geophysical fluid dynamics. It leads, for example, to “potential vorticity thinking” (Hoskins et al. 1985). Mathematically, one imagines an invariant manifold within the phase space of solutions, called the slow manifold (Leith 1980; Lorenz 1980), to which the dynamics are constrained. Whether this slow manifold truly exists has been a major subject of inquiry over the past 20 years. It has become clear that an exact slow manifold is an exceptional case, restricted to steady or perhaps temporally periodic flows (Warn 1997). Thus the concept of a “fuzzy slow manifold” (Warn and Ménard 1986) has been suggested. The idea is that nearly slow dynamics will occur in a stochastic layer about the putative slow manifold. The natural question then is, how thick is this layer? In a recent paper, Ford et al. (2000) argue that Lighthill emission—the spontaneous emission of freely propagating acoustic waves by unsteady vortical flows—is applicable to the problem of balance, with the Mach number Ma replaced by the Froude number F, and that it is a fundamental mechanism for this fuzziness. They consider the rotating shallow-water equations and find emission of inertia–gravity waves at O(F2). This is rather surprising at first sight, because several studies of balanced dynamics with the rotating shallow-water equations have gone beyond second order in F, and found only an exponentially small unbalanced component (Warn and Ménard 1986; Lorenz and Krishnamurthy 1987; Bokhove and Shepherd 1996; Wirosoetisno and Shepherd 2000). We have no technical objection to the analysis of Ford et al. (2000), but wish to point out that it depends crucially on R 1, where R is the Rossby number. This condition requires the ratio of the characteristic length scale of the flow L to the Rossby deformation radius LR to go to zero in the limit F → 0. This is the low Froude number scaling of Charney (1963), which, while originally designed for the Tropics, has been argued to be also relevant to mesoscale dynamics (Riley et al. 1981). If L/LR is fixed, however, then F → 0 implies R → 0, which is the standard quasigeostrophic scaling of Charney (1948; see, e.g., Pedlosky 1987). In this limit there is reason to expect the fuzziness of the slow manifold to be “exponentially thin,” and balance to be much more accurate than is consistent with (algebraic) Lighthill emission.

Fast dynamics in the optical storage with Langmuir-Blodgett films of a diazocrown ether molecule

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Managing cattle grazing intensity: effects on soil organic matter and soil nitrogen

Extensive cattle grazing is the dominant land use in northern Australia. It has been suggested that grazing intensity and rainfall have profound effects on the dynamics of soil nutrients in northern Australia’s semi-arid rangelands. Previous studies have found positive, neutral and negative effects of grazing pressure on soil nutrients. These inconsistencies could be due to short-term experiments that do not capture the slow dynamics of some soil nutrients and the effects of interannual variability in rainfall. In a long-term cattle grazing trial in northern Australia on Brown Sodosol–Yellow Kandosol complex, we analysed soil organic matter and mineral nitrogen in surface soils (0–10 cm depth) 11, 12 and 16 years after trial establishment on experimental plots representing moderate stocking (stocked at the long-term carrying capacity for the region) and heavy stocking (stocked at twice the long-term carrying capacity). Higher soil organic matter was found under heavy stocking, although grazing treatment had little effect on mineral and total soil nitrogen. Interannual variability had a large effect on soil mineral nitrogen, but not on soil organic matter, suggesting that soil nitrogen levels observed in this soil complex may be affected by other indirect pathways, such as climate. The effect of interannual variability in rainfall and the effects of other soil types need to be explored further.

Time evolution of resistance in response to magnetic field: Evidence of glassy transport in La0.85Sr0.15CoO3

We demonstrate the distinct glassy transport phenomena associated with the phase separated and spin-glass-like phases of La0.85Sr0.15CoO3, prepared under different heat-treatment conditions. The low-temperature annealed (phase-separated) sample, exhibits a small change in resistance, with evolution of time, as compared to the high-temperature annealed (spin glass) one. However, the resistance change as a function of time, in both cases, is well described by a stretched exponential fit, signifying the slow dynamics. Moreover, the ultraviolet spectroscopy study evidences a relatively higher density of states in the vicinity of EF for low-temperature annealed sample and this correctly points to its less semiconducting behavior.

Evidence of collective relaxation behavior in the reentrant phase of oxygen-rich LaMnO3

The reentrant low temperature phase of the perovskite manganite LaMnO3+delta (delta=0.22) has been investigated with ac susceptibility and dc magnetization studies. A critical examination of the memory effects in ac susceptibility leads us to the conclusion that the slow dynamics in the system is a consequence of collective relaxation processes resulting from interactions between ferromagnetic clusters, whose presence was indicated in earlier studies. Here, we postulate that the collective behavior is due to the existence of long-range (dipolar) interactions between the large ferromagnetic `superspins'. This is also confirmed by an abnormally large microscopic spin-flip time (similar to 10(-9) s) compared to a canonical spin glass. (C) 2013 Elsevier B.V. All rights reserved.

Memory effect in SrRu(1-x)O3 (0.01 < x < 0.07)

SrRuO3 is widely known to be an itinerant ferromagnet with a T-C similar to 160 K. It is well known that glassy materials exhibit time dependent phenomena such as memory effect due to their generic slow dynamics. However, for the first time, we have observed memory effect in SrRu(1-x)O3 (0.01

Evolution of ferromagnetism from a frustrated state in LixNi(2-x)O2 (0.67 < x < 0.98)

Two-dimensional triangular-lattice antiferromagnetic systems continue to be an interesting area in condensed matter physics and LiNiO2 is one such among them. Here we present a detailed experimental magnetic study of the quasi-stoichiometric LixNi2-xO2 system (0.67slow dynamics (magnetic relaxation, magnetic memory effect etc). By tuning the Li deficiency in a controlled manner, an increase in the ordering temperature is observed. Most strikingly, with the Li deficiency the nature of the magnetic ground state is changed from spin glass to ferromagnetic, with two intermediate states-namely cluster glass and re-entrant spin glass. The critical behaviour of the re-entrant spin glass is also studied here. The critical exponents (beta, gamma and delta) are extracted from the modified Arrot plot, Kouvel-Fisher method, and critical isotherm analysis. The critical exponents match with the long-range mean-field model. The values of the critical exponents are confirmed by the Widom scaling law: delta = 1 + gamma beta(-1). Furthermore, the universality class of the scaling relations is verified, where the scaled m and scaled h collapse into two branches. Finally, based on our observations, a phase diagram is constructed.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.