Rakesh S. Singh

Anisotropy induced crossover from weakly to strongly first order melting of two dimensional solids

Melting and freezing transitions in two dimensional (2D) systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems: the melting of 2D molecular solids is still largely unexplored. In order to understand the role of anisotropy as well as multiple energy and length scales present in molecular systems, here we report computer simulation studies of melting of 2D molecular systems. We computed a limited portion of the solid-liquid phase diagram. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram consisting of isotropic liquid and two crystalline phases-honeycomb and oblique. The nature of the transition depends on the relative strength of the anisotropic interaction and a strongly first order melting turns into a weakly first order transition on increasing the strength of the isotropic interaction. This crossover can be attributed to an increase in stiffness of the solid phase free energy minimum on increasing the strength of the anisotropic interaction. The defects involved in melting of molecular systems are quite different from those known for the atomic systems.

Correlation between thermodynamic anomalies and pathways of ice nucleation in supercooled water.

The well-known classical nucleation theory (CNT) for the free energy barrier towards formation of a nucleus of critical size of the new stable phase within the parent metastable phase fails to take into account the influence of other metastable phases having density/order intermediate between the parent metastable phase and the final stable phase. This lacuna can be more serious than capillary approximation or spherical shape assumption made in CNT. This issue is particularly significant in ice nucleation because liquid water shows rich phase diagram consisting of two (high and low density) liquid phases in supercooled state. The explanations of thermodynamic and dynamic anomalies of supercooled water often invoke the possible influence of a liquid-liquid transition between two metastable liquid phases. To investigate both the role of thermodynamic anomalies and presence of distinct metastable liquid phases in supercooled water on ice nucleation, we employ density functional theoretical approach to find nucleation free energy barrier in different regions of phase diagram. The theory makes a number of striking predictions, such as a dramatic lowering of nucleation barrier due to presence of a metastable intermediate phase and crossover in the dependence of free energy barrier on temperature near liquid-liquid critical point. These predictions can be tested by computer simulations as well as by controlled experiments.

Density and bond-orientational relaxations in supercooled water

ABSTRACT Recent computational studies have reported evidence of a metastable liquid–liquid phase transition (LLPT) in molecular models of water under deeply supercooled conditions. A competing hypothesis suggests, however, that non-equilibrium artefacts associated with coarsening of the stable crystal phase have been mistaken for an LLPT in these models. Such artefacts are posited to arise due to a separation of time scales in which density fluctuations in the supercooled liquid relax orders of magnitude faster than those associated with bond-orientational order. Here, we use molecular simulation to investigate the relaxation of density and bond-orientational fluctuations in three molecular models of water (ST2, TIP5P and TIP4P/2005) in the vicinity of their reported LLPT. For each model, we find that density is the slowly relaxing variable under such conditions. We also observe similar behaviour in the coarse-grained mW model of water. Our findings, therefore, challenge the key physical assumption underlying the competing hypothesis.

Orientational order as the origin of the long-range hydrophobic effect

The long range attractive force between two hydrophobic surfaces immersed in water is observed to decrease exponentially with their separation-this distance-dependence of effective force is known as the hydrophobic force law (HFL). We explore the microscopic origin of HFL by studying distance-dependent attraction between two parallel rods immersed in 2D Mercedes Benz model of water. This model is found to exhibit a well-defined HFL. Although the phenomenon is conventionally explained by density-dependent theories, we identify orientation, rather than density, as the relevant order parameter. The range of density variation is noticeably shorter than that of orientational heterogeneity. The latter is comparable to the observed distances of hydrophobic force. At large separation, attraction between the rods arises primarily from a destructive interference among the inwardly propagating oppositely oriented heterogeneity generated in water by the two rods. As the rods are brought closer, the interference increases leading to a decrease in heterogeneity and concomitant decrease in free energy of the system, giving rise to the effective attraction. We notice formation of hexagonal ice-like structures at the onset of attractive region which suggests that metastable free energy minimum may play a role in the origin of HFL.

Solid-solid collapse transition in a two dimensional model molecular system

Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real difficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises from the interplay between different energy and length scales involved in molecular systems, often mediated through a solvent. In this work we employ Monte-Carlo simulation to study the collapse transition in a model molecular system interacting via both isotropic as well as anisotropic interactions having different length and energy scales. The model we use is known as Mercedes-Benz (MB), which, for a specific set of parameters, sustains two solid phases: honeycomb and oblique. In order to study the temperature induced collapse transition, we start with a metastable honeycomb solid and induce transition by increasing temperature. High density oblique solid so formed has two characteristic length scales corresponding to isotropic and anisotropic parts of interaction potential. Contrary to the common belief and classical nucleation theory, interestingly, we find linear strip-like nucleating clusters having significantly different order and average coordination number than the bulk stable phase. In the early stage of growth, the cluster grows as a linear strip, followed by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate balance between two types of interactions, which enables the dominance of stabilizing energy over destabilizing surface energy. The nucleus of stable oblique phase is wetted by intermediate order particles, which minimizes the surface free energy. In the case of pressure induced transition at low temperature the collapsed state is a disordered solid. The disordered solid phase has diverse local quasi-stable structures along with oblique-solid like domains.

Microscopic Origin of Hysteresis in Water Sorption on Protein Matrices

Despite the importance of water sorption isotherms for a fundamental understanding of protein-water interactions, the microscopic origin of hysteresis between the adsorption and desorption branches is not well understood. Using our recently developed simulation technique, we compute the water sorption isotherms of two proteins, lysozyme and Trp-cage, a miniprotein. We explicitly compare protein-water interactions in adsorption and desorption processes, by analyzing local hydration in terms of hydrogen bonding, water density, and solvent-accessible surface area. We find that significant differences in hydration behavior between adsorption and desorption manifest themselves at the individual amino acid level, in particular around polar or charged residues. We confirm this observation by demonstrating that Trp-cages hysteresis can be significantly reduced by mutating charged residues to alanine, a neutral and nonpolar amino acid.

Sensitivity of nucleation phenomena on range of interaction potential.

Theoretical and computational investigations of nucleation have been plagued by the sensitivity of the phase diagram to the range of the interaction potential. As the surface tension depends strongly on the range of interaction potential and as the classical nucleation theory (CNT) predicts the free energy barrier to be directly proportional to the cube of the surface tension, one expects a strong sensitivity of nucleation barrier to the range of the potential; however, CNT leaves many aspects unexplored. We find for gas-liquid nucleation in Lennard-Jones system that on increasing the range of interaction the kinetic spinodal (KS) (where the mechanism of nucleation changes from activated to barrierless) shifts deeper into the metastable region. Therefore the system remains metastable for larger value of supersaturation and this allows one to explore the high metastable region without encountering the KS. On increasing the range of interaction, both the critical cluster size and pre-critical minima in the free energy surface of kth largest cluster, at respective kinetic spinodals, shift towards smaller cluster size. In order to separate surface tension contribution to the increase in the barrier from other non-trivial factors, we introduce a new scaling form for surface tension and use it to capture both the temperature and the interaction range dependence of surface tension. Surprisingly, we find only a weak non-trivial contribution from other factors to the free energy barrier of nucleation.

Effects of disulfide bridges and backbone connectivity on water sorption by protein matrices

Understanding the water sorption behavior of protein powders is important in applications such as the preservation of protein-based pharmaceuticals. Most globular proteins exhibit a characteristic sigmoidal water adsorption isotherm at ambient conditions. However, it is not well understood how water sorption behavior is influenced by intrinsic factors that are related to structural properties of proteins. We investigate computationally how structural constraints on proteins influence the water sorption isotherms of amorphous protein powders. Specifically, we study the effects of non-local disulfide linkages and backbone connectivity using pheromone ER-23 and lysozyme as model proteins. We find that non-local disulfide linkages can significantly restrict structural changes during hydration and dehydration, and this in turn greatly reduces the extent of hysteresis between the adsorption and desorption branches. Upon removing the backbone connectivity by breaking all peptide bonds in lysozyme, we find that the hysteresis shifts towards the lower humidity regime, and the water uptake capacity is significantly enhanced. We attribute these changes to the higher aggregation propensity of the constraint-free amino acids in dehydrated condition, and the formation of a spanning water network at high hydration levels.