1
|
Malosso C, Manko N, Izzo MG, Baroni S, Hassanali A. Evidence of ferroelectric features in low-density supercooled water from ab initio deep neural-network simulations. Proc Natl Acad Sci U S A 2024; 121:e2407295121. [PMID: 39083416 PMCID: PMC11317578 DOI: 10.1073/pnas.2407295121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2024] [Accepted: 06/30/2024] [Indexed: 08/02/2024] Open
Abstract
Over the last decade, an increasing body of evidence has emerged, supporting the existence of a metastable liquid-liquid critical point in supercooled water whereby two distinct liquid phases of different densities coexist. Analyzing long molecular dynamics simulations performed using deep neural-network force fields trained to accurate quantum mechanical data, we demonstrate that the low-density liquid phase displays a strong propensity toward spontaneous polarization, as witnessed by large and long-lived collective dipole fluctuations. Our findings suggest that the dynamical stability of the low-density phase, and hence the transition from high-density to low-density liquid, is triggered by a collective process involving an accumulation of rotational angular jumps, which could ignite large dipole fluctuations. This dynamical transition involves subtle changes in the electronic polarizability of water molecules which affects their rotational mobility within the two phases. These findings hold the potential for catalyzing activity in the search for dielectric-based probes of the putative second critical point.
Collapse
Affiliation(s)
- Cesare Malosso
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
| | - Natalia Manko
- Condensed Matter and Statistical Physics (CMSP), The Abdus Salam Centre for Theoretical Physics, Trieste34151, Italy
| | - Maria Grazia Izzo
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
| | - Stefano Baroni
- Scuola Internazionale Superiore di Studi Avanzati, Trieste34136, Italy
- Consiglio Nazionale delle Ricerche-Istituto Officina dei Materiali, Scuola Internazionale Superiore di Studi Avanzati Unit, Trieste34136, Italy
| | - Ali Hassanali
- Condensed Matter and Statistical Physics (CMSP), The Abdus Salam Centre for Theoretical Physics, Trieste34151, Italy
| |
Collapse
|
2
|
Lee SC, Z Y. Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity. J Chem Phys 2024; 161:044502. [PMID: 39037136 DOI: 10.1063/5.0215370] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2024] [Accepted: 07/02/2024] [Indexed: 07/23/2024] Open
Abstract
Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl2 from Tm + 85 K to Tm + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl2. We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. This result may contribute to new understandings of the structural relaxation of other network liquids.
Collapse
Affiliation(s)
- Shao-Chun Lee
- Department of Nuclear, Plasma, and Radiological Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Y Z
- Department of Nuclear, Plasma, and Radiological Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Department of Nuclear Engineering and Radiological Science, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Materials Science and Engineering, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
3
|
Lupi L, Gallo P. Mode coupling behavior and fragile to strong transition of trehalose in a binary mixture with water upon supercooling. J Chem Phys 2024; 160:244501. [PMID: 38912627 DOI: 10.1063/5.0218369] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2024] [Accepted: 06/06/2024] [Indexed: 06/25/2024] Open
Abstract
We perform molecular dynamics simulations of a binary mixture of water and trehalose with the TIP4P/Ice water model. We analyze the slow dynamics of trehalose molecules in the mildly supercooled region for concentrations of 3.66 and 18.57 wt. %. We previously studied the dynamics of water in the same mixtures. Supercooled TIP4P/Ice water solvating trehalose molecules was found to follow the Mode Coupling Theory (MCT) and to undergo a transition from a fragile to a strong behavior for both concentrations. Here, we show that also the dynamics of trehalose molecules follows the MCT and displays a fragile to strong crossover (FSC). The results show that trehalose in binary mixtures with water shares with it the dynamical behavior typical of glass forming liquids. Moreover, the FSC for trehalose structural relaxation times is found to occur at temperatures close to those previously obtained for water in the same solutions, showing that the dynamics of the solute is strongly coupled to that of the solvent. We also perform a MCT test showing that the trehalose dynamics obeys the MCT time-temperature superposition principle and that the exponents derived from the theory and the ones obtained from fitting procedure of the relaxation times are comparable, confirming that trehalose molecules in supercooled water solutions follow the MCT of glassy dynamics. Moreover, as predicted by the theory, trehalose particles have MCT parameters comparable to those of water in the same mixtures. This is an important result, given that MCT was originally formulated for monoatomic particles.
Collapse
Affiliation(s)
- Laura Lupi
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - Paola Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| |
Collapse
|
4
|
Iorio A, Perin L, Gallo P. Structure and slow dynamics of protein hydration water with cryopreserving DMSO and trehalose upon cooling. J Chem Phys 2024; 160:244502. [PMID: 38912631 DOI: 10.1063/5.0205569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Accepted: 06/05/2024] [Indexed: 06/25/2024] Open
Abstract
We study, through molecular dynamics simulations, three aqueous solutions with one lysozyme protein and three different concentrations of trehalose and dimethyl sulfoxide (DMSO). We analyze the structural and dynamical properties of the protein hydration water upon cooling. We find that trehalose plays a major role in modifying the structure of the network of HBs between water molecules in the hydration layer of the protein. The dynamics of hydration water presents, in addition to the α-relaxation, typical of glass formers, a slower long-time relaxation process, which greatly slows down the dynamics of water, particularly in the systems with trehalose, where it becomes dominant at low temperatures. In all the solutions, we observe, from the behavior of the α-relaxation times, a shift of the Mode Coupling Theory crossover temperature and the fragile-to-strong crossover temperature toward higher values with respect to bulk water. We also observe a strong-to-strong crossover from the temperature behavior of the long-relaxation times. In the aqueous solution with only DMSO, the transition shifts to a lower temperature than in the case with only lysozyme reported in the literature. We observe that the addition of trehalose to the mixture has the opposite effect of restoring the original location of the strong-to-strong crossover. In all the solutions analyzed in this work, the observed temperature of the protein dynamical transition is slightly shifted at lower temperatures than that of the strong-to-strong crossover, but their relative order is the same, showing a correlation between the motion of the protein and that of the hydration water.
Collapse
Affiliation(s)
- Antonio Iorio
- Dipartimento di Fisica, Università Roma Tre, Via della Vasca Navale 84, I-00146 Roma, Italy
| | - Leonardo Perin
- Dipartimento di Fisica, Università Roma Tre, Via della Vasca Navale 84, I-00146 Roma, Italy
| | - Paola Gallo
- Dipartimento di Fisica, Università Roma Tre, Via della Vasca Navale 84, I-00146 Roma, Italy
| |
Collapse
|
5
|
Tavagnacco L, Zanatta M, Buratti E, Bertoldo M, Chiessi E, Appel M, Natali F, Orecchini A, Zaccarelli E. Water slowing down drives the occurrence of the low temperature dynamical transition in microgels. Chem Sci 2024; 15:9249-9257. [PMID: 38903230 PMCID: PMC11186305 DOI: 10.1039/d4sc02650k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/05/2024] [Indexed: 06/22/2024] Open
Abstract
The protein dynamical transition marks an increase in atomic mobility and the onset of anharmonic motions at a critical temperature (T d), which is considered relevant for protein functionality. This phenomenon is ubiquitous, regardless of protein composition, structure and biological function and typically occurs at large protein content, to avoid water crystallization. Recently, a dynamical transition has also been reported in non-biological macromolecules, such as poly(N-isopropyl acrylamide) (PNIPAM) microgels, bearing many similarities to proteins. While the generality of this phenomenon is well-established, the role of water in the transition remains a subject of debate. In this study, we use atomistic molecular dynamics (MD) simulations and elastic incoherent neutron scattering (EINS) experiments with selective deuteration to investigate the microscopic origin of the dynamical transition and distinguish water and PNIPAM roles. While a standard analysis of EINS experiments would suggest that the dynamical transition occurs in PNIPAM and water at a similar temperature, simulations reveal a different perspective, also qualitatively supported by experiments. From room temperature down to about 180 K, PNIPAM exhibits only modest changes of dynamics, while water, being mainly hydration water under the probed extreme confinement, significantly slows down and undergoes a mode-coupling transition from diffusive to activated. Our findings therefore challenge the traditional view of the dynamical transition, demonstrating that it occurs in proximity of the water mode-coupling transition, shedding light on the intricate interplay between polymer and water dynamics.
Collapse
Affiliation(s)
- Letizia Tavagnacco
- CNR Institute of Complex Systems, Uos Sapienza Piazzale Aldo Moro 2 00185 Rome Italy
- Department of Physics, Sapienza University of Rome Piazzale Aldo Moro 2 00185 Rome Italy
| | - Marco Zanatta
- Department of Physics, University of Trento Via Sommarive 14 38123 Trento Italy
| | - Elena Buratti
- Department of Chemical, Pharmaceutical and Agricultural Sciences, University of Ferrara Via L. Borsari 46 44121 Ferrara Italy
| | - Monica Bertoldo
- Department of Chemical, Pharmaceutical and Agricultural Sciences, University of Ferrara Via L. Borsari 46 44121 Ferrara Italy
| | - Ester Chiessi
- Department of Chemical Science and Technologies, University of Rome Tor Vergata Via della Ricerca Scientifica I 00133 Rome Italy
| | - Markus Appel
- Institut Laue-Langevin 71 avenue des Martyrs, CS 20156 38042 Grenoble Cedex 9 France
| | - Francesca Natali
- CNR-IOM, Operative Group Grenoble (OGG), Institut Laue Langevin F-38042 Grenoble France
| | - Andrea Orecchini
- Dipartimento di Fisica e Geologia, Università di Perugia Via Alessandro Pascoli 06123 Perugia Italy
- CNR-IOM c/o Dipartimento di Fisica e Geologia, Università di Perugia Via Alessandro Pascoli 06123 Perugia Italy
| | - Emanuela Zaccarelli
- CNR Institute of Complex Systems, Uos Sapienza Piazzale Aldo Moro 2 00185 Rome Italy
- Department of Physics, Sapienza University of Rome Piazzale Aldo Moro 2 00185 Rome Italy
| |
Collapse
|
6
|
Kumar S, Sarkar S, Bagchi B. Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum. J Chem Phys 2024; 160:224501. [PMID: 38856065 DOI: 10.1063/5.0210699] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Accepted: 05/24/2024] [Indexed: 06/11/2024] Open
Abstract
In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.
Collapse
Affiliation(s)
- Shubham Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Sarmistha Sarkar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Biman Bagchi
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
7
|
Voulgarakis NK. Multilayered noise model for transport in complex environments. Phys Rev E 2023; 108:064105. [PMID: 38243501 DOI: 10.1103/physreve.108.064105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Accepted: 11/09/2023] [Indexed: 01/21/2024]
Abstract
Transport in complex fluidic environments often exhibits transient subdiffusive dynamics accompanied by non-Gaussian probability density profiles featuring a nonmonotonic non-Gaussian parameter. Such properties cannot be adequately explained by the original theory of Brownian motion. Based on an extension of kinetic theory, this study introduces a chain of hierarchically coupled random walks approach that effectively captures all these intriguing characteristics. If the environment consists of a series of independent white noise sources, then the problem can be expressed as a system of hierarchically coupled Ornstein-Uhlenbech equations. Due to the linearity of the system, the most essential transport properties have a closed analytical form.
Collapse
Affiliation(s)
- Nikolaos K Voulgarakis
- Department of Mathematics and Statistics, Washington State University, Pullman, Washington 99164, USA
| |
Collapse
|
8
|
Lupi L, Gallo P. Glassy dynamics of water in TIP4P/Ice aqueous solutions of trehalose in comparison with the bulk phase. J Chem Phys 2023; 159:154504. [PMID: 37850697 DOI: 10.1063/5.0168933] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 09/25/2023] [Indexed: 10/19/2023] Open
Abstract
We perform molecular dynamics simulations of TIP4P/Ice water in solution with trehalose for 3.65 and 18.57 wt. % concentrations and of bulk TIP4P/Ice water at ambient pressure, to characterize the structure and dynamics of water in a sugar aqueous solution in the supercooled region. We find here that TIP4P/Ice water in solution with trehalose molecules follows the Mode Coupling Theory and undergoes a fragile to strong transition up to the highest concentration investigated, similar to the bulk. Moreover, we perform a Mode Coupling Theory test, showing that the Time Temperature Superposition principle holds for both bulk TIP4P/Ice water and for TIP4P/Ice water in the solutions and we calculate the exponents of the theory. The direct comparison of the dynamical results for bulk water and water in the solutions shows upon cooling along the isobar a fastening of water dynamics for lower temperatures, T < 240 K. We found that the counter-intuitive behavior for the low temperature solutions can be explained with the diffusion anomaly of water leading us to the conclusion that the fastening observed below T = 240 K in water dynamics is only fictitious, due to the fact that the density of water molecules in the solutions is higher than the density of the bulk at the same temperature and pressure. This result should be taken into account in experimental investigations which are often carried out at constant pressure.
Collapse
Affiliation(s)
- Laura Lupi
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - Paola Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| |
Collapse
|
9
|
Wang J, Xie SJ. The influence of force fields on the structure and dynamics of water confined in ZIF-8 from atomistic simulations. Phys Chem Chem Phys 2023; 25:23100-23110. [PMID: 37602670 DOI: 10.1039/d3cp02075d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/22/2023]
Abstract
The complexity of modeling flexible crystals, such as ZIF-8, mainly stems from the handling of intramolecular interactions. Numerous force fields have been proposed in the literature to describe the interactions between atoms in ZIF-8. We employ seven force fields to examine the structure and dynamic behavior of water molecules confined in ZIF-8, with the aim of investigating the impact of force fields on simulation results. Various structural characterization methods consistently indicate that the choice of different force fields has quantitative effects but no qualitative effects on the structural characteristics of confined water. Additionally, the force fields do not impact the qualitative description of the diffusion mechanism. Both mean-square displacement and van Hove autocorrelation function reveal two characteristic movements of water molecules diffusing in ZIF-8: a short-time intra-cavity hopping process and a long-time inter-cavity hopping process. However, the framework flexibility is found to play a crucial role in determining the order of spatial arrangement and local structure, self-diffusion coefficient and reorientational dynamics of confined water. Specifically, the DREIDING force field gives rise to an unrealistic stiff framework, enhancing the order of spatial arrangement and diminishing the local ordered structure of confined water. Meanwhile, it results in much slower translational and reorientational dynamics. Hence, the general DREIDING force field cannot be considered for providing a quantitative description of the water structure and dynamics.
Collapse
Affiliation(s)
- Jing Wang
- Center for Membrane Separation and Water Science & Technology, Zhejiang University of Technology, Hangzhou 310014, P. R. China.
| | - Shi-Jie Xie
- Center for Membrane Separation and Water Science & Technology, Zhejiang University of Technology, Hangzhou 310014, P. R. China.
| |
Collapse
|
10
|
Loubet NA, Verde AR, Lockhart JA, Appignanesi GA. Turning an energy-based defect detector into a multi-molecule structural indicator for water. J Chem Phys 2023; 159:064512. [PMID: 37578063 DOI: 10.1063/5.0159060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2023] [Accepted: 07/27/2023] [Indexed: 08/15/2023] Open
Abstract
Recent studies have provided conclusive evidence for the existence of a liquid-liquid critical point in numerical models of water. Such a scenario implies the competition between two local molecular arrangements of different densities: a high-density liquid (HDL) and a low-density liquid (LDL). Within this context, the development of accurate structural indicators to properly characterize the two interconverting local structures is demanded. In a previous study, we introduced a reliable energy-based structural descriptor that properly discriminates water molecules into tetrahedrally arranged molecules (T molecules) and distorted molecules (D molecules). The latter constitute defects in terms of hydrogen bond (HB) coordination and have been shown to represent a minority component, even at high temperatures above the melting point. In addition, the D molecules tend to form high-quality HBs with three T molecules and to be surrounded by T and D molecules at further distances. Thus, it became evident that, while the LDL state might consist of a virtually pure T state, the HDL state would comprise mixed molecular arrangements including the D molecules. Such a need to abandon the single-molecule description requires the investigation of the degree of structural information to be incorporated in order to build an appropriate multi-molecule indicator. Hence, in this work, we shall study the effect of the local structural constraints on the water molecules in order to discriminate the different molecular arrangements into two disjoint classes. This will enable us to build a multi-molecule structural indicator for water whose performance will then be investigated within the water's supercooled regime.
Collapse
Affiliation(s)
- Nicolás A Loubet
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Alejandro R Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Jano A Lockhart
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
11
|
Li R, Zhu P, Xu Y, Lu H, Rong J. Molecular Dynamics Simulation of Droplet Impact on a Hydrophobic 3D Elastic Surface. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023. [PMID: 37450274 DOI: 10.1021/acs.langmuir.3c01519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/18/2023]
Abstract
The phenomenon of droplets impacting elastic surfaces is common in nature and in many engineering applications. It has been shown that droplet impact on an elastic surface drastically reduces droplet contact time and hinders droplet spreading. However, most of the current studies are based on experiments, and the analysis of the influence mechanism of the elastic substrate on the dynamic behavior of droplets is not complete. In addition, the simulations of droplet impact on elastic substrates are mainly focused on 2D elastic films or vibrating rigid substrates, ignoring the effect of 3D elastic substrate deformation on the droplet dynamic behavior. Therefore, in this paper, we propose to model the droplet impact on a 3D hydrophobic elastic substrate using the molecular dynamics method. We find that droplet pancake rebound can substantially reduce the droplet contact time. Moreover, we record the conditions required for the pancake rebound of the droplet. Furthermore, we investigated the effects of the elastic modulus of the substrate and the initial velocity of the droplet on the droplet contact time, contact area, and spreading factor. This study further elucidates the influence mechanism of the elastic substrate on the dynamic behavior of the droplet and provides theoretical guidance for regulating the dynamic behavior of the droplet in related fields.
Collapse
Affiliation(s)
- Rao Li
- School of Mechanical Engineering, Beijing Institute of Technology, Beijing 10081, PR China
| | - Pengzhe Zhu
- School of Mechanical Engineering, Beijing Institute of Technology, Beijing 10081, PR China
| | - Yimeng Xu
- School of Mechanical Engineering, Beijing Institute of Technology, Beijing 10081, PR China
| | - Hongsheng Lu
- School of Mechanical Engineering, Beijing Institute of Technology, Beijing 10081, PR China
| | - Jiacheng Rong
- School of Mechanical Engineering, Beijing Institute of Technology, Beijing 10081, PR China
| |
Collapse
|
12
|
Shikata K, Kikutsuji T, Yasoshima N, Kim K, Matubayasi N. Revealing the hidden dynamics of confined water in acrylate polymers: Insights from hydrogen-bond lifetime analysis. J Chem Phys 2023; 158:2887576. [PMID: 37125720 DOI: 10.1063/5.0148753] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 04/13/2023] [Indexed: 05/02/2023] Open
Abstract
Polymers contain functional groups that participate in hydrogen bond (H-bond) with water molecules, establishing a robust H-bond network that influences bulk properties. This study utilized molecular dynamics (MD) simulations to examine the H-bonding dynamics of water molecules confined within three poly(meth)acrylates: poly(2-methoxyethyl acrylate) (PMEA), poly(2-hydroxyethyl methacrylate) (PHEMA), and poly(1-methoxymethyl acrylate) (PMC1A). Results showed that H-bonding dynamics significantly slowed as the water content decreased. Additionally, the diffusion of water molecules and its correlation with H-bond breakage were analyzed. Our findings suggest that when the H-bonds between water molecules and the methoxy oxygen of PMEA are disrupted, those water molecules persist in close proximity and do not diffuse on a picosecond time scale. In contrast, the water molecules H-bonded with the hydroxy oxygen of PHEMA and the methoxy oxygen of PMC1A diffuse concomitantly with the breakage of H-bonds. These results provide an in-depth understanding of the impact of polymer functional groups on H-bonding dynamics.
Collapse
Affiliation(s)
- Kokoro Shikata
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Takuma Kikutsuji
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuhiro Yasoshima
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Department of Information and Computer Engineering, National Institute of Technology, Toyota College, 2-1 Eiseicho, Toyota, Aichi 471-8525, Japan
| | - Kang Kim
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
13
|
Malik S, Karmakar S, Debnath A. Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers. J Chem Phys 2023; 158:114503. [PMID: 36948835 DOI: 10.1063/5.0138681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023] Open
Abstract
The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.
Collapse
Affiliation(s)
- Sheeba Malik
- Department of Chemistry, IIT Jodhpur, Jodhpur, Rajasthan, India
| | - Smarajit Karmakar
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad, India
| | - Ananya Debnath
- Department of Chemistry, IIT Jodhpur, Jodhpur, Rajasthan, India
| |
Collapse
|
14
|
Verde AR, Alarcón LM, Appignanesi GA. Correlations between defect propensity and dynamical heterogeneities in supercooled water. J Chem Phys 2023; 158:114502. [PMID: 36948825 DOI: 10.1063/5.0139118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023] Open
Abstract
A salient feature of supercooled liquids consists in the dramatic dynamical slowdown they undergo as temperature decreases while no significant structural change is evident. These systems also present dynamical heterogeneities (DH): certain molecules, spatially arranged in clusters, relax various orders of magnitude faster than the others. However, again, no static quantity (such as structural or energetic measures) shows strong direct correlations with such fast-moving molecules. In turn, the dynamic propensity approach, an indirect measure that quantifies the tendency of the molecules to move in a given structural configuration, has revealed that dynamical constraints, indeed, originate from the initial structure. Nevertheless, this approach is not able to elicit which structural quantity is, in fact, responsible for such a behavior. In an effort to remove dynamics from its definition in favor of a static quantity, an energy-based propensity has also been developed for supercooled water, but it could only find positive correlations between the lowest-energy and the least-mobile molecules, while no correlations could be found for those more relevant mobile molecules involved in the DH clusters responsible for the system's structural relaxation. Thus, in this work, we shall define a defect propensity measure based on a recently introduced structural index that accurately characterizes water structural defects. We shall show that this defect propensity measure provides positive correlations with dynamic propensity, being also able to account for the fast-moving molecules responsible for the structural relaxation. Moreover, time dependent correlations will show that defect propensity represents an appropriate early-time predictor of the long-time dynamical heterogeneity.
Collapse
Affiliation(s)
- Alejandro R Verde
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Laureano M Alarcón
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000 Bahía Blanca, Argentina
| |
Collapse
|
15
|
Krishnan G, Harbola U. Quantum uncertainty effects in the dynamics of supercooled liquids: A molecular dynamics study. Phys Rev E 2022; 106:064604. [PMID: 36671091 DOI: 10.1103/physreve.106.064604] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2022] [Accepted: 11/20/2022] [Indexed: 06/17/2023]
Abstract
Dynamics of density fluctuations in quantum supercooled liquids is analyzed using molecular dynamics simulations. In contrast to the classical case, the uncertainty in the particle position (delocalization of quantum particle in space) leads to significant differences in the dynamics of quantum liquids, both in the short- and long-time limits. The effect of uncertainty is found to be significant for length scales smaller than the uncertainty itself, and diminishes as the length scale grows. The dynamic heterogeneity of the system at short times is enhanced due to uncertainty. In the intermediate (β-relaxation) time regime, the heterogeneity tends to get suppressed due to quantum uncertainty. The probability distribution of particle displacements shows highly nonclassical behavior with double-peak structure at short timescales.
Collapse
Affiliation(s)
- Gopika Krishnan
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| | - Upendra Harbola
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
16
|
Gomez A, Piskulich ZA, Thompson WH, Laage D. Water Diffusion Proceeds via a Hydrogen-Bond Jump Exchange Mechanism. J Phys Chem Lett 2022; 13:4660-4666. [PMID: 35604934 DOI: 10.1021/acs.jpclett.2c00825] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The self-diffusion of water molecules plays a key part in a broad range of essential processes in biochemistry, medical imaging, material science, and engineering. However, its molecular mechanism and the role played by the water hydrogen-bond network rearrangements are not known. Here we combine molecular dynamics simulations and analytic modeling to determine the molecular mechanism of water diffusion. We establish a quantitative connection between the water diffusion coefficient and hydrogen-bond jump exchanges, and identify the features that determine the underlying energetic barrier. We thus provide a unified framework to understand the coupling between translational, rotational, and hydrogen-bond dynamics in liquid water. It explains why these different dynamics do not necessarily exhibit identical temperature dependences although they all result from the same hydrogen-bond exchange events. The consequences for the understanding of water diffusion in supercooled conditions and for water transport in complex aqueous systems, including ionic, biological, and confined solutions, are discussed.
Collapse
Affiliation(s)
- Axel Gomez
- PASTEUR, Department of Chemistry, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Damien Laage
- PASTEUR, Department of Chemistry, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
17
|
Roy P, Menon S, Sengupta N. Dynamical Manifestations of Supercooling in Amyloid Hydration. J Phys Chem B 2021; 126:44-53. [PMID: 34941279 DOI: 10.1021/acs.jpcb.1c07724] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The effect of extreme temperature on amyloidogenic species remains sparsely explored. In a recent study (J. Phys. Chem. Lett., 2019, 10, (10)), we employed exhaustive molecular dynamics simulations to explore the cold thermal response of a putative small amyloid oligomer and to elicit the role of solvent modulation. Herein, we investigate the dynamical response of the hydration waters of the oligomer within the supercooled states. Using NMR-based formalism, we delineate the entropic response in terms of the side-chain conformational entropy that corroborates the weakening of the hydrophobic core with lowering of temperature. The translational dynamics of the protein and hydration waters reveal the coupling of protein dynamical fluctuations with solvent dynamics under supercooled conditions. Probing the translational motion as a space-time correlation indicates glassy dynamics exhibited by hydration waters in the supercooled regime. Caging of the water molecules with lowering of temperature and the resultant hopping dynamics are reflected in the longer β-relaxation timescales of translational motion. Furthermore, we utilized mode-coupling theory (MCT) and derived the ideal glass transition temperature from translational and rotational dynamics, around ∼196 and 209 K, respectively. Interestingly, rotational motion in the supercooled regime deviates from the MCT law, exhibits Arrhenius motion, and marks a fragile-to-strong crossover at 227 K. The low-frequency vibrational modes also coincide with the dynamical transition. This exposition lends dynamical insights into the hydration coupling of an amyloid aggregate under cryogenic conditions.
Collapse
Affiliation(s)
- Priti Roy
- Department of Biological Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal, India 741246
| | - Sneha Menon
- Tata Institute of Fundamental Research Hyderabad, Telangana 500046, India
| | - Neelanjana Sengupta
- Department of Biological Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal, India 741246
| |
Collapse
|
18
|
Gallo P, Bachler J, Bove LE, Böhmer R, Camisasca G, Coronas LE, Corti HR, de Almeida Ribeiro I, de Koning M, Franzese G, Fuentes-Landete V, Gainaru C, Loerting T, de Oca JMM, Poole PH, Rovere M, Sciortino F, Tonauer CM, Appignanesi GA. Advances in the study of supercooled water. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:143. [PMID: 34825973 DOI: 10.1140/epje/s10189-021-00139-1] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 10/17/2021] [Indexed: 06/13/2023]
Abstract
In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.
Collapse
Affiliation(s)
- Paola Gallo
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy.
| | - Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Livia E Bove
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
- Sorbonne Université, CNRS UMR 7590, IMPMC, 75005, Paris, France
| | - Roland Böhmer
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Gaia Camisasca
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Luis E Coronas
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Horacio R Corti
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina
| | - Ingrid de Almeida Ribeiro
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859, Campinas, São Paulo, Brazil
| | - Maurice de Koning
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859, Campinas, São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861, Campinas, São Paulo, Brazil
| | - Giancarlo Franzese
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Violeta Fuentes-Landete
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Catalin Gainaru
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | | | - Peter H Poole
- Department of Physics, St. Francis Xavier University, Antigonish, NS, B2G 2W5, Canada
| | - Mauro Rovere
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
| | - Christina M Tonauer
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
| |
Collapse
|
19
|
Corti HR, Appignanesi GA, Barbosa MC, Bordin JR, Calero C, Camisasca G, Elola MD, Franzese G, Gallo P, Hassanali A, Huang K, Laria D, Menéndez CA, de Oca JMM, Longinotti MP, Rodriguez J, Rovere M, Scherlis D, Szleifer I. Structure and dynamics of nanoconfined water and aqueous solutions. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:136. [PMID: 34779954 DOI: 10.1140/epje/s10189-021-00136-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Accepted: 10/06/2021] [Indexed: 06/13/2023]
Abstract
This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.
Collapse
Affiliation(s)
- Horacio R Corti
- Departmento de Física de la Materia Condensada & Instituto de Nanociencia y Nanotecnología (CNEA-CONICET), Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina.
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, 8000, Bahía Blanca, Argentina
| | - Marcia C Barbosa
- Institute of Physics, Federal University of Rio Grande do Sul, 91501-970, Porto Alegre, Brazil
| | - J Rafael Bordin
- Department of Physics, Institute of Physics and Mathematics, 96050-500, Pelotas, RS, Brazil
| | - Carles Calero
- Secció de Física Estadística i Interdisciplinària - Departament de Física de la Matèria Condensada, Universitat de Barcelona & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, 08028, Barcelona, Spain
| | - Gaia Camisasca
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, 00146, Roma, Italy
| | - M Dolores Elola
- Departmento de Física de la Materia Condensada & Instituto de Nanociencia y Nanotecnología (CNEA-CONICET), Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina
| | - Giancarlo Franzese
- Secció de Física Estadística i Interdisciplinària - Departament de Física de la Matèria Condensada, Universitat de Barcelona & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, 08028, Barcelona, Spain
| | - Paola Gallo
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, 00146, Roma, Italy
| | - Ali Hassanali
- Condensed Matter and Statistical Physics Section (CMSP), The International Center for Theoretical Physics (ICTP), Trieste, Italy
| | - Kai Huang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong, China
| | - Daniel Laria
- Departmento de Física de la Materia Condensada & Instituto de Nanociencia y Nanotecnología (CNEA-CONICET), Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina
- Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Cintia A Menéndez
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, 8000, Bahía Blanca, Argentina
| | - Joan M Montes de Oca
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, 8000, Bahía Blanca, Argentina
| | - M Paula Longinotti
- Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Javier Rodriguez
- Departmento de Física de la Materia Condensada & Instituto de Nanociencia y Nanotecnología (CNEA-CONICET), Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina
- Escuela de Ciencia y Tecnología, Universidad Nacional de General San Martín, San Martín, Buenos Aires, Argentina
| | - Mauro Rovere
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, 00146, Roma, Italy
| | - Damián Scherlis
- Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Igal Szleifer
- Biomedical Engineering Department, Northwestern University, Evanston, USA
| |
Collapse
|
20
|
Tetrahedral structure of supercooled water at ambient pressure and its influence on dynamic relaxation: Comparative study of water models. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117269] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|
21
|
Lupi L, Vázquez Ramírez B, Gallo P. Dynamical crossover and its connection to the Widom line in supercooled TIP4P/Ice water. J Chem Phys 2021; 155:054502. [PMID: 34364341 DOI: 10.1063/5.0059190] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We perform molecular dynamics simulations with the TIP4P/Ice water model to characterize the relationship between dynamics and thermodynamics of liquid water in the supercooled region. We calculate the relevant properties of the phase diagram, and we find that TIP4P/Ice presents a retracing line of density maxima, similar to what was previously found for atomistic water models and models of other tetrahedral liquids. For this model, a liquid-liquid critical point between a high-density liquid and a low-density liquid was recently found. We compute the lines of the maxima of isothermal compressibility and the minima of the coefficient of thermal expansion in the one phase region, and we show that these lines point to the liquid-liquid critical point while collapsing on the Widom line. This line is the line of the maxima of correlation length that emanates from a second order critical point in the one phase region. Supercooled water was found to follow mode coupling theory and to undergo a transition from a fragile to a strong behavior right at the crossing of the Widom line. We find here that this phenomenology also happens for TIP4P/Ice. Our results appear, therefore, to be a general characteristic of supercooled water, which does not depend on the interaction potential used, and they reinforce the idea that the dynamical crossover from a region where the relaxation mechanism is dominated by cage relaxation to a region where cages are frozen and hopping dominates is correlated in water to a phase transition between a high-density liquid and a low-density liquid.
Collapse
Affiliation(s)
- Laura Lupi
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - Benjamín Vázquez Ramírez
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - Paola Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| |
Collapse
|
22
|
Kikutsuji T, Kim K, Matubayasi N. Transition pathway of hydrogen bond switching in supercooled water analyzed by the Markov state model. J Chem Phys 2021; 154:234501. [PMID: 34241244 DOI: 10.1063/5.0055531] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
In this work, we examine hydrogen-bond (H-bond) switching by employing the Markov State Model (MSM). During the H-bond switching, a water hydrogen initially H-bonded with water oxygen becomes H-bonded to a different water oxygen. MSM analysis was applied to trajectories generated from molecular dynamics simulations of the TIP4P/2005 model from a room-temperature state to a supercooled state. We defined four basis states to characterize the configuration between two water molecules: H-bonded ("H"), unbound ("U"), weakly H-bonded ("w"), and alternative H-bonded ("a") states. A 16 × 16 MSM matrix was constructed, describing the transition probability between states composed of three water molecules. The mean first-passage time of the H-bond switching was estimated by calculating the total flux from the HU to UH states. It is demonstrated that the temperature dependence of the mean first-passage time is in accordance with that of the H-bond lifetime determined from the H-bond correlation function. Furthermore, the flux for the H-bond switching is decomposed into individual pathways that are characterized by different forms of H-bond configurations of trimers. The dominant pathway of the H-bond switching is found to be a direct one without passing through such intermediate states as "w" and "a," the existence of which becomes evident in supercooled water. The pathway through "w" indicates a large reorientation of the donor molecule. In contrast, the pathway through "a" utilizes the tetrahedral H-bond network, which is revealed by the further decomposition based on the H-bond number of the acceptor molecule.
Collapse
Affiliation(s)
- Takuma Kikutsuji
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kang Kim
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
23
|
Camisasca G, De Marzio M, Gallo P. Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water. J Chem Phys 2021; 153:224503. [PMID: 33317300 DOI: 10.1063/5.0033526] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.
Collapse
Affiliation(s)
- Gaia Camisasca
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - Margherita De Marzio
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - Paola Gallo
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| |
Collapse
|
24
|
Sepulveda-Medina PI, Tyagi M, Wang C, Vogt BD. Water dynamics within nanostructured amphiphilic statistical copolymers from quasielastic neutron scattering. J Chem Phys 2021; 154:154903. [PMID: 33887940 DOI: 10.1063/5.0045341] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Understanding the properties of water under either soft or hard confinement has been an area of great interest, but nanostructured amphiphilic polymers that provide a secondary confinement have garnered significantly less attention. Here, a series of statistical copolymers of 2-hydroxyethyl acrylate (HEA) and 2-(N-ethylperfluorooctane sulfonamido)ethyl methacrylate (FOSM) are swollen to equilibrium in water to form nanostructured physically cross-linked hydrogels to probe the effect of soft confinement on the dynamics of water. Changing the composition of the copolymer from 10 to 21 mol. % FOSM decreases the average size of the assembled FOSM cross-link, but also the spacing between the cross-links in the hydrogels with the mean distance between the FOSM aggregates decreasing from 3.9 to 2.7 nm. The dynamics of water within the hydrogels were assessed with quasielastic neutron scattering. These hydrogels exhibit superior performance for inhibition of water crystallization on supercooling in comparison to analogous hydrogels with different hydrophilic copolymer chemistries. Despite the lower water crystallinity, the self-diffusion coefficient for these hydrogels from the copolymers of HEA and FOSM decreases precipitously below 260 K, which is a counter to the nearly temperature invariant water dynamics reported previously with an analogous hydrogel [Wiener et al., J. Phys. Chem. B 120, 5543 (2016)] that exhibits nearly temperature invariant dynamics to 220 K. These results point to chemistry dependent dynamics of water that is confined within amphiphilic hydrogels, where the interactions of water with the hydrophilic segments can qualitatively alter the temperature dependent dynamics of water in the supercooled state.
Collapse
Affiliation(s)
| | - Madhusudan Tyagi
- Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Chao Wang
- School of Polymer Science and Polymer Engineering, University of Akron, Akron, Ohio 44325, USA
| | - Bryan D Vogt
- School of Polymer Science and Polymer Engineering, University of Akron, Akron, Ohio 44325, USA
| |
Collapse
|
25
|
Shen Z, Ma J, Carrillo JMY, Chen WR, Sumpter BG, Wang Y. Spatiotemporal mapping of mesoscopic liquid dynamics. Phys Rev E 2021; 103:022609. [PMID: 33736070 DOI: 10.1103/physreve.103.022609] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 01/29/2021] [Indexed: 11/07/2022]
Abstract
The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.
Collapse
Affiliation(s)
- Zhiqiang Shen
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Jihong Ma
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Jan-Michael Y Carrillo
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Wei-Ren Chen
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Bobby G Sumpter
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Yangyang Wang
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
26
|
Arbe A, Alvarez F, Colmenero J. Insight into the Structure and Dynamics of Polymers by Neutron Scattering Combined with Atomistic Molecular Dynamics Simulations. Polymers (Basel) 2020; 12:E3067. [PMID: 33371357 PMCID: PMC7767341 DOI: 10.3390/polym12123067] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2020] [Revised: 12/04/2020] [Accepted: 12/11/2020] [Indexed: 11/23/2022] Open
Abstract
Combining neutron scattering and fully atomistic molecular dynamics simulations allows unraveling structural and dynamical features of polymer melts at different length scales, mainly in the intermolecular and monomeric range. Here we present the methodology developed by us and the results of its application during the last years in a variety of polymers. This methodology is based on two pillars: (i) both techniques cover approximately the same length and time scales and (ii) the classical van Hove formalism allows easily calculating the magnitudes measured by neutron scattering from the simulated atomic trajectories. By direct comparison with experimental results, the simulated cell is validated. Thereafter, the information of the simulations can be exploited, calculating magnitudes that are experimentally inaccessible or extending the parameters range beyond the experimental capabilities. We show how detailed microscopic insight on structural features and dynamical processes of various kinds has been gained in polymeric systems with different degrees of complexity, and how intriguing questions as the collective behavior at intermediate length scales have been faced.
Collapse
Affiliation(s)
- Arantxa Arbe
- Centro de Física de Materiales (CSIC, UPV/EHU) and Materials Physics Center MPC, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain; (A.A.); (F.A.)
| | - Fernando Alvarez
- Centro de Física de Materiales (CSIC, UPV/EHU) and Materials Physics Center MPC, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain; (A.A.); (F.A.)
- Departamento de Polímeros y Materiales Avanzados, Física, Química y Tecnología (UPV/EHU), Apartado 1072, E-20080 San Sebastián, Spain
| | - Juan Colmenero
- Centro de Física de Materiales (CSIC, UPV/EHU) and Materials Physics Center MPC, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain; (A.A.); (F.A.)
- Departamento de Polímeros y Materiales Avanzados, Física, Química y Tecnología (UPV/EHU), Apartado 1072, E-20080 San Sebastián, Spain
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, E-20018 San Sebastián, Spain
| |
Collapse
|
27
|
Tenuzzo L, Camisasca G, Gallo P. Protein-Water and Water-Water Long-Time Relaxations in Protein Hydration Water upon Cooling-A Close Look through Density Correlation Functions. Molecules 2020; 25:molecules25194570. [PMID: 33036320 PMCID: PMC7583983 DOI: 10.3390/molecules25194570] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Revised: 09/29/2020] [Accepted: 09/30/2020] [Indexed: 11/16/2022] Open
Abstract
We report results on the translational dynamics of the hydration water of the lysozyme protein upon cooling obtained by means of molecular dynamics simulations. The self van Hove functions and the mean square displacements of hydration water show two different temperature activated relaxation mechanisms, determining two dynamic regimes where transient trapping of the molecules is followed by hopping phenomena to allow to the structural relaxations. The two caging and hopping regimes are different in their nature. The low-temperature hopping regime has a time scale of tenths of nanoseconds and a length scale on the order of 2–3 water shells. This is connected to the nearest-neighbours cage effect and restricted to the supercooling, it is absent at high temperature and it is the mechanism to escape from the cage also present in bulk water. The second hopping regime is active at high temperatures, on the nanoseconds time scale and over distances of nanometers. This regime is connected to water displacements driven by the protein motion and it is observed very clearly at high temperatures and for temperatures higher than the protein dynamical transition. Below this temperature, the suppression of protein fluctuations largely increases the time-scale of the protein-related hopping phenomena at least over 100 ns. These protein-related hopping phenomena permit the detection of translational motions of hydration water molecules longly persistent in the hydration shell of the protein.
Collapse
|
28
|
Prasad M, English NJ, Nath Chakraborty S. Relaxation dynamics and power spectra of liquid water: a molecular dynamics simulation study. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1733117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Mahabir Prasad
- Department of Chemistry, Sikkim University, Gangtok, India
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | | |
Collapse
|
29
|
Jung G, Caraglio M, Schrack L, Franosch T. Dynamical properties of densely packed confined hard-sphere fluids. Phys Rev E 2020; 102:012612. [PMID: 32795038 DOI: 10.1103/physreve.102.012612] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Accepted: 06/28/2020] [Indexed: 11/07/2022]
Abstract
Numerical solutions of the mode-coupling theory (MCT) equations for a hard-sphere fluid confined between two parallel hard walls are elaborated. The governing equations feature multiple parallel relaxation channels which significantly complicate their numerical integration. We investigate the intermediate scattering functions and the susceptibility spectra close to structural arrest and compare to an asymptotic analysis of the MCT equations. We corroborate that the data converge in the β-scaling regime to two asymptotic power laws, viz. the critical decay and the von Schweidler law. The numerical results reveal a nonmonotonic dependence of the power-law exponents on the slab width and a nontrivial kink in the low-frequency susceptibility spectra. We also find qualitative agreement of these theoretical results to event-driven molecular dynamics simulations of polydisperse hard-sphere systems. In particular, the nontrivial dependence of the dynamical properties on the slab width is well reproduced.
Collapse
Affiliation(s)
- Gerhard Jung
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21A, A-6020 Innsbruck, Austria
| | - Michele Caraglio
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21A, A-6020 Innsbruck, Austria
| | - Lukas Schrack
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21A, A-6020 Innsbruck, Austria
| | - Thomas Franosch
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21A, A-6020 Innsbruck, Austria
| |
Collapse
|
30
|
McMillin PJ, Alegrete M, Peric M, Luchko T. Electron Paramagnetic Resonance Measurements of Four Nitroxide Probes in Supercooled Water Explained by Molecular Dynamics Simulations. J Phys Chem B 2020; 124:3962-3972. [PMID: 32301326 DOI: 10.1021/acs.jpcb.0c00684] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Electron paramagnetic resonance (EPR) measurements of the rotational diffusion of small nitroxide probes have been demonstrated to be a powerful technique for experimentally investigating the properties of supercooled liquids, such as water. However, since only the rotational diffusion of the probe molecules is measured and EPR measurements are indirect, it is not clear what the relationship between the behavior of water and the probe molecule is. To address this, we have performed molecular dynamics simulations of four nitroxide probes in TIP4P-Ew and OPC water models to directly compare with EPR experiments and to determine the behavior of the water and the underlying microscopic coupling between the water and the probes. In all, 200 ns simulations were run for 23 temperatures between 253 and 283 K for all four probes with each water model for an aggregate of 36.8 μs of simulation time. Simulations for both water models systematically underestimated the rotational diffusion coefficients for both water and probes, though OPC simulations were generally in better agreement with the experiments than TIP4P-Ew simulations. Despite this, when the temperature dependence of the data was fit to a power law, fit parameters for TIP4P-Ew were generally in better agreement with the experiments than OPC. For probe molecules, the singular temperature was found to be T0 = 226.5 ± 0.4 K from experiments, T0 = 208 ± 2 K for OPC water, and T0 = 215 ± 2 K for TIP4P-Ew water. While for water molecules, the singular temperature was found to be T0 = 220.3 ± 0.2 K from experiments, T0 = 208 ± 2 K for OPC water, and T0 = 220 ± 1 K for TIP4P-Ew water. Systematic underestimation of the rotational diffusion coefficients was most pronounced at lower temperatures and was clearly observed in changes to the Arrhenius activation energy. Above the maximum density temperature of Tρmax = 277 K, an activation energy of EA ≈ 16.7 kJ/mol was observed for the probes from experiments, while OPC had EA ≈ 15.2 kJ/mol and TIP4P-Ew had EA ≈ 14.6 kJ/mol. Below the maximum density temperature, the activation energy jumped to EA ≈ 32.5 kJ/mol for experiments but only EA ≈ 23 kJ/mol for OPC and EA ≈ 22 kJ/mol for TIP4P-Ew. In all cases, we saw good agreement between the behavior of the probe molecules and water. To understand why, we calculated the average number of hydrogen bonds between the probe molecules and water. From this, we were able to explain the rotational diffusion times for all of the probes. These results show that current molecular models are sufficient to capture physical phenomena observed with EPR and to help elucidate why the probes provide accurate insights into the behavior of supercooled water.
Collapse
Affiliation(s)
- Patrick J McMillin
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| | - Matthew Alegrete
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| | - Miroslav Peric
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| | - Tyler Luchko
- Department of Physics and Astronomy, Center for Biological Physics, California State University, Northridge, Northridge, California 91330, United States
| |
Collapse
|
31
|
Ortiz de Urbina J, Sesé G. Orientational dynamics in methanol: Influence of temperature and hydrogen bonding. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2019.112374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
32
|
Srivastava A, Karmakar S, Debnath A. Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling. SOFT MATTER 2019; 15:9805-9815. [PMID: 31746927 DOI: 10.1039/c9sm01725a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.
Collapse
Affiliation(s)
- Abhinav Srivastava
- Department of Chemistry, Indian Institute of Technology Jodhpur, Jodhpur 342037, India.
| | | | | |
Collapse
|
33
|
Consistency of geometrical definitions of hydrogen bonds based on the two-dimensional potential of mean force with respect to the time correlation in liquid water over a wide range of temperatures. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111603] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
34
|
Structure, Dynamics, and Wettability of Water at Metal Interfaces. Sci Rep 2019; 9:14805. [PMID: 31616006 PMCID: PMC6794260 DOI: 10.1038/s41598-019-51323-5] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2019] [Accepted: 09/27/2019] [Indexed: 11/24/2022] Open
Abstract
The water/metal interface often governs important chemophysical processes in various technologies. Therefore, from scientific and engineering perspectives, the detailed molecular-level elucidation of the water/metal interface is of high priority, but the related research is limited. In experiments, the surface-science techniques, which can provide full structural details of the surface, are not easy to directly apply to the interfacial systems under ambient conditions, and the well-defined facets cannot be entirely free from contamination at the contact with water. To answer long-standing debates regarding the wettability, structure, and dynamics of water at metal interfaces, we here develop reliable first-principles-based multiscale simulations. Using the state-of-the-art simulations, we find that the clean metal surfaces are actually superhydrophilic and yield zero contact angles. Furthermore, we disclose an inadequacy of widespread ice-like bilayer model of the water adlayers on metal surfaces from both averaged structural and dynamic points of view. Our findings on the nature of water on metal surfaces provide new molecular level perspectives on the tuning and design of water/metal interfaces that are at the heart of many energy applications.
Collapse
|
35
|
Voloshin V, Smolin N, Geiger A, Winter R, Medvedev NN. Dynamics of TMAO and urea in the hydration shell of the protein SNase. Phys Chem Chem Phys 2019; 21:19469-19479. [PMID: 31461098 DOI: 10.1039/c9cp03184g] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Using all-atom molecular dynamics simulations of aqueous solutions of the globular protein SNase, the dynamic behavior of water molecules and cosolvents (trimethylamine-N-oxide (TMAO) and urea) in the hydration shell of the protein was studied for different solvent compositions. TMAO is a potent protein-stabilizing osmolyte, whereas urea is known to destabilize proteins. For molecules that are initially located in successive narrow layers at a given distance from the protein, the mean displacements and the distribution of displacements for short time intervals are calculated. For molecules that are initially located in solvation shells of a given thickness around the protein, the characteristic residence times in these shells are determined to characterize the dynamic behavior of the solvent molecules as a function of the distance to the protein. A combined consideration of these characteristics allows to reveal additional features of the dynamics of the cosolvents. It is shown that TMAO molecules leave the nearest vicinity of the protein faster than urea molecules, despite the fact that the mobility of TMAO molecules, measured by their mean displacements, is lower than that of urea. Moreover, we show that the rate of release of TMAO molecules from the hydration shell is lower in ternary (TMAO + urea + H2O) solvent mixtures than in the binary ones. This is consistent with a recent observation that the fraction of TMAO near the protein decreases in the presence of urea. From the analysis of the decay of the number of particles initially located in the region of the first peak of the distribution function of solvent molecules around the protein, we estimated that about 20 water molecules and 6-7 urea molecules stay near the protein for more than 1000 ps.
Collapse
Affiliation(s)
- Vladimir Voloshin
- Institute of Chemical Kinetics and Combustion, SB RAS, 630090 Novosibirsk, Russia.
| | - Nikolai Smolin
- Department of Cell and Molecular Physiology, Loyola University Chicago, Maywood, Illinois 60153, USA
| | - Alfons Geiger
- Physikalische Chemie, Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44221 Dortmund, Germany.
| | - Roland Winter
- Physikalische Chemie, Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44221 Dortmund, Germany.
| | - Nikolai N Medvedev
- Institute of Chemical Kinetics and Combustion, SB RAS, 630090 Novosibirsk, Russia. and Novosibirsk State University, 630090 Novosibirsk, Russia
| |
Collapse
|
36
|
Iorio A, Camisasca G, Rovere M, Gallo P. Characterization of hydration water in supercooled water-trehalose solutions: The role of the hydrogen bonds network. J Chem Phys 2019; 151:044507. [PMID: 31370561 DOI: 10.1063/1.5108579] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The structural and dynamical properties of hydration water in aqueous solutions of trehalose are studied with molecular dynamics simulation. We simulate the systems in the supercooled region to investigate how the interaction with the trehalose molecules modifies the hydrogen bond network, the structural relaxation, and the diffusion properties of hydration water. The analysis is performed by considering the radial distribution functions, the residence time of water molecules in the hydration shell, the two body excess entropy, and the hydrogen bond water-water and water-trehalose correlations of the hydration water. The study of the two body excess entropy shows the presence of a fragile to strong crossover in supercooled hydration water also found in the relaxation time of the water-water hydrogen bond correlation function, and this is in agreement with predictions of the mode coupling theory and of previous studies of the oxygen-oxygen density correlators [A. Iorio et al., J. Mol. Liq. 282, 617 (2019); Sci. China: Phys., Mech. Astron. 62, 107011 (2019)]. The water-trehalose hydrogen bond correlation function instead evidences a strong to strong crossover in the relaxation time, and this crossover is related to a trehalose dynamical transition. This signals the role that the strong interplay between the soluted molecules and the surrounding solvent has in determining the dynamical transition common to both components of the system that happens upon cooling and that is similar to the well known protein dynamical transition. We connect our results with the cryoprotecting role of trehalose molecules.
Collapse
Affiliation(s)
- A Iorio
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - G Camisasca
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - M Rovere
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - P Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| |
Collapse
|
37
|
Song S, Park SJ, Kim M, Kim JS, Sung BJ, Lee S, Kim JH, Sung J. Transport dynamics of complex fluids. Proc Natl Acad Sci U S A 2019; 116:12733-12742. [PMID: 31175151 PMCID: PMC6600932 DOI: 10.1073/pnas.1900239116] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Thermal motion in complex fluids is a complicated stochastic process but ubiquitously exhibits initial ballistic, intermediate subdiffusive, and long-time diffusive motion, unless interrupted. Despite its relevance to numerous dynamical processes of interest in modern science, a unified, quantitative understanding of thermal motion in complex fluids remains a challenging problem. Here, we present a transport equation and its solutions, which yield a unified quantitative explanation of the mean-square displacement (MSD), the non-Gaussian parameter (NGP), and the displacement distribution of complex fluids. In our approach, the environment-coupled diffusion kernel and its time correlation function (TCF) are the essential quantities that determine transport dynamics and characterize mobility fluctuation of complex fluids; their time profiles are directly extractable from a model-free analysis of the MSD and NGP or, with greater computational expense, from the two-point and four-point velocity autocorrelation functions. We construct a general, explicit model of the diffusion kernel, comprising one unbound-mode and multiple bound-mode components, which provides an excellent approximate description of transport dynamics of various complex fluidic systems such as supercooled water, colloidal beads diffusing on lipid tubes, and dense hard disk fluid. We also introduce the concepts of intrinsic disorder and extrinsic disorder that have distinct effects on transport dynamics and different dependencies on temperature and density. This work presents an unexplored direction for quantitative understanding of transport and transport-coupled processes in complex disordered media.
Collapse
Affiliation(s)
- Sanggeun Song
- Creative Research Initiative Center for Chemical Dynamics in Living Cells, Chung-Ang University, 06974 Seoul, Republic of Korea
- Department of Chemistry, Chung-Ang University, 06974 Seoul, Republic of Korea
- National Institute of Innovative Functional Imaging, Chung-Ang University, 06974 Seoul, Republic of Korea
| | - Seong Jun Park
- Creative Research Initiative Center for Chemical Dynamics in Living Cells, Chung-Ang University, 06974 Seoul, Republic of Korea
- Department of Chemistry, Chung-Ang University, 06974 Seoul, Republic of Korea
- National Institute of Innovative Functional Imaging, Chung-Ang University, 06974 Seoul, Republic of Korea
| | - Minjung Kim
- Department of Chemistry, College of Natural Sciences, Seoul National University, 08826 Seoul, Republic of Korea
| | - Jun Soo Kim
- Department of Chemistry and Nanoscience, Ewha Womans University, 03760 Seoul, Republic of Korea
| | - Bong June Sung
- Department of Chemistry, Sogang University, 04107 Seoul, Republic of Korea
| | - Sangyoub Lee
- Department of Chemistry, College of Natural Sciences, Seoul National University, 08826 Seoul, Republic of Korea
| | - Ji-Hyun Kim
- Creative Research Initiative Center for Chemical Dynamics in Living Cells, Chung-Ang University, 06974 Seoul, Republic of Korea;
| | - Jaeyoung Sung
- Creative Research Initiative Center for Chemical Dynamics in Living Cells, Chung-Ang University, 06974 Seoul, Republic of Korea;
- Department of Chemistry, Chung-Ang University, 06974 Seoul, Republic of Korea
- National Institute of Innovative Functional Imaging, Chung-Ang University, 06974 Seoul, Republic of Korea
| |
Collapse
|
38
|
Camisasca G, Galamba N, Wikfeldt KT, Pettersson LGM. Translational and rotational dynamics of high and low density TIP4P/2005 water. J Chem Phys 2019; 150:224507. [PMID: 31202216 DOI: 10.1063/1.5079956] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We use molecular dynamics simulations using TIP4P/2005 to investigate the self- and distinct-van Hove functions for different local environments of water, classified using the local structure index as an order parameter. The orientational dynamics were studied through the calculation of the time-correlation functions of different-order Legendre polynomials in the OH-bond unit vector. We found that the translational and orientational dynamics are slower for molecules in a low-density local environment and correspondingly the mobility is enhanced upon increasing the local density, consistent with some previous works, but opposite to a recent study on the van Hove function. From the analysis of the distinct dynamics, we find that the second and fourth peaks of the radial distribution function, previously identified as low density-like arrangements, show long persistence in time. The analysis of the time-dependent interparticle distance between the central molecule and the first coordination shell shows that particle identity persists longer than distinct van Hove correlations. The motion of two first-nearest-neighbor molecules thus remains coupled even when this correlation function has been completely decayed. With respect to the orientational dynamics, we show that correlation functions of molecules in a low-density environment decay exponentially, while molecules in a local high-density environment exhibit bi-exponential decay, indicating that dynamic heterogeneity of water is associated with the heterogeneity among high-density and between high-density and low-density species. This bi-exponential behavior is associated with the existence of interstitial waters and the collapse of the second coordination sphere in high-density arrangements, but not with H-bond strength.
Collapse
Affiliation(s)
- Gaia Camisasca
- Department of Physics, Stockholm University, 106 91 Stockholm, Sweden
| | - Nuno Galamba
- Centre of Chemistry and Biochemistry and Biosystems and Integrative Sciences Institute, Faculty of Sciences of the University of Lisbon, C8, Campo Grande, 1749-016 Lisbon, Portugal
| | | | | |
Collapse
|
39
|
Kikutsuji T, Kim K, Matubayasi N. Diffusion dynamics of supercooled water modeled with the cage-jump motion and hydrogen-bond rearrangement. J Chem Phys 2019; 150:204502. [DOI: 10.1063/1.5095978] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Takuma Kikutsuji
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kang Kim
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
40
|
Iorio A, Camisasca G, Gallo P. Slow dynamics of hydration water and the trehalose dynamical transition. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.02.088] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
41
|
M R, Ayappa KG. Enhancing the Dynamics of Water Confined between Graphene Oxide Surfaces with Janus Interfaces: A Molecular Dynamics Study. J Phys Chem B 2019; 123:2978-2993. [DOI: 10.1021/acs.jpcb.8b12341] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
|
42
|
Tavagnacco L, Chiessi E, Zanatta M, Orecchini A, Zaccarelli E. Water-Polymer Coupling Induces a Dynamical Transition in Microgels. J Phys Chem Lett 2019; 10:870-876. [PMID: 30735054 PMCID: PMC6416711 DOI: 10.1021/acs.jpclett.9b00190] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2019] [Accepted: 02/08/2019] [Indexed: 06/09/2023]
Abstract
The long debated protein dynamical transition was recently found also in nonbiological macromolecules, such as poly- N-isopropylacrylamide (PNIPAM) microgels. Here, by using atomistic molecular dynamics simulations, we report a description of the molecular origin of the dynamical transition in these systems. We show that PNIPAM and water dynamics below the dynamical transition temperature T d are dominated by methyl group rotations and hydrogen bonding, respectively. By comparing with bulk water, we unambiguously identify PNIPAM-water hydrogen bonding as mainly responsible for the occurrence of the transition. The observed phenomenology thus crucially depends on the water-macromolecule coupling, being relevant to a wide class of hydrated systems, independently from the biological function.
Collapse
Affiliation(s)
- Letizia Tavagnacco
- CNR-ISC
and Department of Physics, Sapienza University
of Rome, Piazzale A.
Moro 2, 00185 Rome, Italy
| | - Ester Chiessi
- Department
of Chemical Sciences and Technologies, University
of Rome Tor Vergata, Via della Ricerca Scientica I, 00133 Rome, Italy
| | - Marco Zanatta
- Department
of Computer Science, University of Verona, Strada Le Grazie 15, 37138 Verona, Italy
| | - Andrea Orecchini
- Department
of Physics and Geology, University of Perugia, Via A. Pascoli, 06123 Perugia, Italy
- CNR-IOM
c/o Department of Physics and Geology, University
of Perugia, Via A. Pascoli, 06123 Perugia, Italy
| | - Emanuela Zaccarelli
- CNR-ISC
and Department of Physics, Sapienza University
of Rome, Piazzale A.
Moro 2, 00185 Rome, Italy
| |
Collapse
|
43
|
Shi R, Russo J, Tanaka H. Common microscopic structural origin for water's thermodynamic and dynamic anomalies. J Chem Phys 2018; 149:224502. [PMID: 30553247 DOI: 10.1063/1.5055908] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water displays a vast array of unique properties, known as water's anomalies, whose origin remains subject to hot debate. Our aim in this article is to provide a unified microscopic physical picture of water's anomalies in terms of locally favored structures, encompassing both thermodynamic and dynamic anomalies, which are often attributed to different origins. We first identify locally favored structures via a microscopic structural descriptor that measures local translational order and provide direct evidence that they have a hierarchical impact on the anomalies. At each state point, the strength of thermodynamic anomalies is directly proportional to the amount of locally favored structures, while the dynamic properties of each molecule depend on the local structure surrounding both itself and its nearest neighbors. To incorporate this, we develop a novel hierarchical two-state model. We show by extensive simulations of two popular water models that both thermodynamic and kinetic anomalies can be almost perfectly explained by the temperature and pressure dependence of these local and non-local versions of the same structural descriptor, respectively. Moreover, our scenario makes three unique predictions in supercooled water, setting it apart from other scenarios: (1) Presence of an "Arrhenius-to-Arrhenius" crossover upon cooling, as the origin of the apparent "fragile-to-strong" transition; (2) maximum of dynamic heterogeneity around 20 K below the Widom line and far above the glass transition; (3) violation of the Stokes-Einstein-Debye relation at ∼2T g, rather than 1.2T g typical of normal glass-formers. These predictions are verified by recent measurement of water's diffusion at very low temperatures (point 1) and discoveries from our extensive simulations (points 2-3). We suggest that the same scenario may generally apply to water-like anomalies in liquids tending to form locally favored structures, including not only other important tetrahedral liquids such as silicon, germanium, and silica, but also metallic and chalcogenide liquids.
Collapse
Affiliation(s)
- Rui Shi
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - John Russo
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
44
|
|
45
|
Tsimpanogiannis IN, Moultos OA, Franco LFM, Spera MBDM, Erdős M, Economou IG. Self-diffusion coefficient of bulk and confined water: a critical review of classical molecular simulation studies. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1511903] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Luís F. M. Franco
- School of Chemical Engineering, University of Campinas, Campinas, Brazil
| | | | - Máté Erdős
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Ioannis G. Economou
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
46
|
Zanatta M, Tavagnacco L, Buratti E, Bertoldo M, Natali F, Chiessi E, Orecchini A, Zaccarelli E. Evidence of a low-temperature dynamical transition in concentrated microgels. SCIENCE ADVANCES 2018; 4:eaat5895. [PMID: 30276264 PMCID: PMC6162076 DOI: 10.1126/sciadv.aat5895] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2018] [Accepted: 08/21/2018] [Indexed: 05/07/2023]
Abstract
A low-temperature dynamical transition has been reported in several proteins. We provide the first observation of a "protein-like" dynamical transition in nonbiological aqueous environments. To this aim, we exploit the popular colloidal system of poly-N-isopropylacrylamide (PNIPAM) microgels, extending their investigation to unprecedentedly high concentrations. Owing to the heterogeneous architecture of the microgels, water crystallization is avoided in concentrated samples, allowing us to monitor atomic dynamics at low temperatures. By elastic incoherent neutron scattering and molecular dynamics simulations, we find that a dynamical transition occurs at a temperature T d ~ 250 K, independently from PNIPAM mass fraction. However, the transition is smeared out on approaching dry conditions. The quantitative agreement between experiments and simulations provides evidence that the transition occurs simultaneously for PNIPAM and water dynamics. The similarity of these results with hydrated protein powders suggests that the dynamical transition is a generic feature in complex macromolecular systems, independently from their biological function.
Collapse
Affiliation(s)
- Marco Zanatta
- Department of Computer Science, University of Verona, Strada le Grazie 15, 37134 Verona, Italy
| | - Letizia Tavagnacco
- Consiglio Nazionale delle Ricerche (CNR)–Istituto dei Sistemi Complessi, Sede Sapienza, and Department of Physics, Sapienza University of Rome, Piazzale A. Moro 2, 00185 Rome, Italy
| | - Elena Buratti
- CNR–Istituto per i Processi Chimico-Fisici, Sede Secondaria di Pisa, Area della Ricerca, via G. Moruzzi 1, 56124 Pisa, Italy
| | - Monica Bertoldo
- CNR–Istituto per i Processi Chimico-Fisici, Sede Secondaria di Pisa, Area della Ricerca, via G. Moruzzi 1, 56124 Pisa, Italy
- Corresponding author. (M.B.); (A.O.); (E.Z.)
| | - Francesca Natali
- CNR–Istituto Officina dei Materiali (IOM), Operative Group in Grenoble, c/o Institut Laue Langevin, 6 rue Jules Horowitz, BP 156, 38042 Grenoble Cedex 9, France
| | - Ester Chiessi
- Department of Chemical Sciences and Technologies, University of Rome Tor Vergata, Via della Ricerca Scientifica I, 00133 Rome, Italy
| | - Andrea Orecchini
- Department of Physics and Geology, University of Perugia, Via A. Pascoli, 06123 Perugia, Italy
- CNR-IOM c/o Department of Physics and Geology, University of Perugia, Via A. Pascoli, 06123 Perugia, Italy
- Corresponding author. (M.B.); (A.O.); (E.Z.)
| | - Emanuela Zaccarelli
- Consiglio Nazionale delle Ricerche (CNR)–Istituto dei Sistemi Complessi, Sede Sapienza, and Department of Physics, Sapienza University of Rome, Piazzale A. Moro 2, 00185 Rome, Italy
- Corresponding author. (M.B.); (A.O.); (E.Z.)
| |
Collapse
|
47
|
Kikutsuji T, Kim K, Matubayasi N. How do hydrogen bonds break in supercooled water?: Detecting pathways not going through saddle point of two-dimensional potential of mean force. J Chem Phys 2018; 148:244501. [DOI: 10.1063/1.5033419] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Takuma Kikutsuji
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kang Kim
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
48
|
Stachura SS, Malajczuk CJ, Mancera RL. Molecular dynamics simulations of a DMSO/water mixture using the AMBER force field. J Mol Model 2018; 24:174. [DOI: 10.1007/s00894-018-3720-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2018] [Accepted: 06/12/2018] [Indexed: 12/23/2022]
|
49
|
Camisasca G, De Marzio M, Rovere M, Gallo P. High density liquid structure enhancement in glass forming aqueous solution of LiCl. J Chem Phys 2018; 148:222829. [DOI: 10.1063/1.5024375] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- G. Camisasca
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. De Marzio
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. Rovere
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| | - P. Gallo
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| |
Collapse
|
50
|
Espinosa JR, Navarro C, Sanz E, Valeriani C, Vega C. On the time required to freeze water. J Chem Phys 2018; 145:211922. [PMID: 28799362 DOI: 10.1063/1.4965427] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avrami's expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude.
Collapse
Affiliation(s)
- J R Espinosa
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Navarro
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Valeriani
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|