The homogeneous crystallization of water at low temperature is believed to occur through the direct nucleation of cubic (Ic) and hexagonal (Ih) ices. Here, we provide evidence from molecular ...simulations that the nucleation of ice proceeds through the formation of a new metastable phase, which we name Ice 0. We find that Ice 0 is structurally similar to the supercooled liquid, and that on growth it gradually converts into a stacking of Ice Ic and Ih. We suggest that this mechanism provides a thermodynamic explanation for the location and pressure dependence of the homogeneous nucleation temperature, and that Ice 0 controls the homogeneous nucleation of low-pressure ices, acting as a precursor to crystallization in accordance with Ostwald's step rule of phases. Our findings show that metastable crystalline phases of water may play roles that have been largely overlooked.
Some liquids, if cooled rapidly enough to avoid crystallization, can be frozen into a nonergodic glassy state. The tendency for a material to form a glass when quenched is called “glass-forming ...ability,” and it is of key significance both fundamentally and for materials science applications. Here, we consider liquids with competing orderings, where an increase in the glass-forming ability is signaled by a depression of the melting temperature towards its minimum at triple or eutectic points. With simulations of two model systems where glass-forming ability can be tuned by an external parameter, we are able to interpolate between crystal-forming and glass-forming behavior. We find that the enhancement of the glass-forming ability is caused by an increase in the structural difference between liquid and crystal: stronger competition in orderings towards the melting point minimum makes a liquid structure more disordered (more complex). This increase in the liquid-crystal structure difference can be described by a single adimensional parameter, i.e., the interface energy cost scaled by the thermal energy, which we call the “thermodynamic interface penalty.” Our finding may provide a general physical principle for not only controlling the glass-forming ability but also the emergence of glassy behavior of various systems with competing orderings, including orderings of structural, magnetic, electronic, charge, and dipolar origin.
We introduce a sequence-dependent parametrization for a coarse-grained DNA model T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, J. Chem. Phys. 134, 085101 (2011) originally designed to reproduce ...the properties of DNA molecules with average sequences. The new parametrization introduces sequence-dependent stacking and base-pairing interaction strengths chosen to reproduce the melting temperatures of short duplexes. By developing a histogram reweighting technique, we are able to fit our parameters to the melting temperatures of thousands of sequences. To demonstrate the flexibility of the model, we study the effects of sequence on: (a) the heterogeneous stacking transition of single strands, (b) the tendency of a duplex to fray at its melting point, (c) the effects of stacking strength in the loop on the melting temperature of hairpins, (d) the force-extension properties of single strands, and (e) the structure of a kissing-loop complex. Where possible, we compare our results with experimental data and find a good agreement. A simulation code called oxDNA, implementing our model, is available as a free software.
DNA is acquiring a primary role in material development, self-assembling by design into complex supramolecular aggregates, the building block of a new-materials world. Using DNA nanoconstructs to ...translate sophisticated theoretical intuitions into experimental realizations by closely matching idealized models of colloidal particles is a much less explored avenue. Here we experimentally show that an appropriate selection of competing interactions enciphered in multiple DNA sequences results into the successful design of a one-pot DNA hydrogel that melts both on heating and on cooling. The relaxation time, measured by light scattering, slows down dramatically in a limited window of temperatures. The phase diagram displays a peculiar re-entrant shape, the hallmark of the competition between different bonding patterns. Our study shows that it is possible to rationally design biocompatible bulk materials with unconventional phase diagrams and tuneable properties by encoding into DNA sequences both the particle shape and the physics of the collective response.
Although the thermodynamics of DNA hybridization is generally well established, the kinetics of this classic transition is less well understood. Providing such understanding has new urgency because ...DNA nanotechnology often depends critically on binding rates. Here, we explore DNA oligomer hybridization kinetics using a coarse-grained model. Strand association proceeds through a complex set of intermediate states, with successful binding events initiated by a few metastable base-pairing interactions, followed by zippering of the remaining bonds. But despite reasonably strong interstrand interactions, initial contacts frequently dissociate because typical configurations in which they form differ from typical states of similar enthalpy in the double-stranded equilibrium ensemble. Initial contacts must be stabilized by two or three base pairs before full zippering is likely, resulting in negative effective activation enthalpies. Non-Arrhenius behavior arises because the number of base pairs required for nucleation increases with temperature. In addition, we observe two alternative pathways-pseudoknot and inchworm internal displacement-through which misaligned duplexes can rearrange to form duplexes. These pathways accelerate hybridization. Our results explain why experimentally observed association rates of GC-rich oligomers are higher than rates of AT- rich equivalents, and more generally demonstrate how association rates can be modulated by sequence choice.
We present a new, nucleotide-level model for RNA, oxRNA, based on the coarse-graining methodology recently developed for the oxDNA model of DNA. The model is designed to reproduce structural, ...mechanical, and thermodynamic properties of RNA, and the coarse-graining level aims to retain the relevant physics for RNA hybridization and the structure of single- and double-stranded RNA. In order to explore its strengths and weaknesses, we test the model in a range of nanotechnological and biological settings. Applications explored include the folding thermodynamics of a pseudoknot, the formation of a kissing loop complex, the structure of a hexagonal RNA nanoring, and the unzipping of a hairpin motif. We argue that the model can be used for efficient simulations of the structure of systems with thousands of base pairs, and for the assembly of systems of up to hundreds of base pairs. The source code implementing the model is released for public use.