Search

Several computational studies on different water models reported evidence of a phase transition in supercooled conditions between two liquid states of water differing in density: the high-density liquid (HDL) and the low-density liquid (LDL). Yet, conclusive experimental evidence of the existence of a phase transition between the two liquid water phases could not be obtained due to fast crystallization in the region where the phase transition should occur. For the same reason, the investigation of possible transition mechanisms between the two phases is committed to computational investigations. In this work, we simulate an out-of-equilibrium temperature-induced transition from the LDL to the HDL-like state in the TIP4P/2005 water model. To structurally characterize the system relaxation, we use the node total communicability (NTC) we recently proposed as an effective order parameter to discriminate the two liquid phases differing in density. We find that the relaxation process is compatible with a spinodal-like scenario. We observe the formation of HDL-like domains in the LDL phase and we characterize their fluctuating behavior and subsequent coarsening and stabilization. Furthermore, we find that the formation of stable HDL-like domains is favored in the regions where the early formation of small patches of highly connected HDL-like molecules (i.e., with very high NTC values) is observed. Besides characterizing the LDL- to HDL-like relaxation from a structural point of view, these results also show that the NTC order parameter can serve as an early-time predictor of the regions from which the transition process initiates. [ABSTRACT FROM AUTHOR]

In this paper, we apply a theoretical model for fluid state thermodynamics to investigate simulated water in supercooled conditions. This model, which we recently proposed and applied to sub- and super-critical fluid water [Zanetti-Polzi et al., J. Chem. Phys. 156(4), 44506 (2022)], is based on a combination of the moment-generating functions of the enthalpy and volume fluctuations as provided by two gamma distributions and provides the free energy of the system as well as other relevant thermodynamic quantities. The application we make here provides a thermodynamic description of supercooled water fully consistent with that expected by crossing the liquid–liquid Widom line, indicating the presence of two distinct liquid states. In particular, the present model accurately reproduces the Widom line temperatures estimated with other two-state models and well describes the heat capacity anomalies. Differently from previous models, according to our description, a cluster of molecules that extends beyond the first hydration shell is necessary to discriminate between the statistical fluctuation regimes typical of the two liquid states. [ABSTRACT FROM AUTHOR]

Estimating the kinetics of electron transfer (ET) processes in biologically relevant systems using theoretical‐computational methods remains a formidable task. This challenge arises from the inherent complexity of these systems, which makes it impractical to apply a fully quantum‐mechanical treatment. Hybrid quantum mechanical/classical mechanical computational approaches have been devised to enable the explicit simulation of electron transfer kinetics. This concept article focuses on a specific theoretical‐computational method employed in this context, namely the Perturbed Matrix Method (PMM), which has the merit of being able to include large‐scale conformational effects in the ET kinetics and potential multiple, alternative, ET channels. We describe its underlying physical principles, examine its advantages and limitations, and offer insights into its applications. Examples of the approach are discussed in the context of estimating photo‐induced electron transfer kinetics in proteins. The non‐exponential behavior observed in the presented case studies mainly arises from an active coupling with the environment fluctuations, but also partly stems from the presence of branching ET pathways. [ABSTRACT FROM AUTHOR]

We propose in this paper a theoretical model for fluid state thermodynamics based on modeling the fluctuation distributions and, hence, the corresponding moment generating functions providing the free energy of the system. Using the relatively simple and physically coherent gamma model for the fluctuation distributions, we obtain a complete theoretical equation of state, also giving insight into the statistical/molecular organization and phase or pseudo-phase transitions occurring under the sub- and super-critical conditions, respectively. Application to sub- and super-critical fluid water and a comparison with the experimental data show that this model provides an accurate description of fluid water thermodynamics, except close to the critical point region where limited but significant deviations from the experimental data occur. We obtain quantitative evidence of the correspondence between the sub- and super-critical thermodynamic behaviors, with the super-critical water pseudo-liquid and pseudo-gas phases being the evolution of the sub-critical water liquid and gas phases, respectively. Remarkably, according to our model, we find that for fluid water the minimal subsystem corresponding to either the liquid-like or the gas-like condition includes an infinite number of molecules in the sub-critical regime (providing the expected singularities due to macroscopic phase transitions) but only five molecules in the super-critical regime (coinciding with the minimal possible hydrogen-bonding cluster), thus suggesting that the super-critical regime be characterized by the coexistence of nanoscopic subsystems in either the pseudo-liquid or the pseudo-gas phase with each subsystem fluctuating between forming and disrupting the minimal hydrogen-bonding network. [ABSTRACT FROM AUTHOR]

The most intriguing hypothesis explaining many water anomalies is a metastable liquid–liquid phase transition (LLPT) at high pressure and low temperatures, experimentally hidden by homogeneous nucleation. Recent infrared spectroscopic experiments showed that upon addition of hydrazinium trifluoroacetate to water, the supercooled ionic solution undergoes a sharp, reversible LLPT at ambient pressure, possible offspring of that in pure water. Here, we calculate the temperature-dependent signature of the OH-stretching band, reporting on the low/high density phase of water, in neat water and in the same experimentally investigated ionic solution. The comparison between the infrared signature of the pure liquid and that of the ionic solution can be achieved only computationally, providing insight into the nature of the experimentally observed phase transition and allowing us to investigate the effects of ionic compounds on the high to low density supercooled liquid water transition. We show that the experimentally observed crossover behavior in the ionic solution can be reproduced only if the phase transition between the low- and high-density liquid states of water is coupled to a mixing–unmixing transition between the water component and the ions: at low temperatures, water and ions are separated and the water component is a low density liquid. At high temperatures, water and ions get mixed and the water component is a high-density liquid. The separation at low temperatures into ion-rich and ion-poor regions allows unveiling the polyamorphic nature of liquid water, leading to a crossover behavior resembling that observed in supercooled neat water under high pressure. [ABSTRACT FROM AUTHOR]

Infrared spectroscopy is a widely used technique to characterize protein structures and protein mediated processes. While the amide I band provides information on proteins' secondary structure, amino acid side chains are used as infrared probes for the investigation of protein reactions and local properties. In this paper, we use a hybrid quantum mechanical/classical molecular dynamical approach based on the perturbed matrix method to compute the infrared band due to the C=O stretching mode of amide-containing side chains. We calculate, at first, the infrared band of zwitterionic glutamine in water and obtain results in very good agreement with the experimental data. Then, we compute the signal arising from glutamine side chains in a microcrystal of the yeast prion Sup35-derived peptide, GNNQQNY, with a fibrillar structure. The infrared bands obtained by selective isotopic labeling of the two glutamine residues, Q4 and Q5, of each peptide were experimentally used to investigate the local hydration in the fibrillar microcrystal. The experimental spectra of the two glutamine residues, which experience different hydration environments, feature different spectral signals that are well reproduced by the corresponding calculated spectra. In addition, the analysis of the simulated spectra clarifies the molecular origin of the experimentally observed spectroscopic differences that arise from the different local electric field experienced by the two glutamine residues, which is, in turn, determined by a different hydrogen bonding pattern. [ABSTRACT FROM AUTHOR]

Liquid water, besides being fundamental for life on Earth, has long fascinated scientists due to several anomalies. Different hypotheses have been put forward to explain these peculiarities. The most accredited one foresees the presence in the supercooled region of two phases at different densities: the low-density liquid phase and the high-density liquid phase. In our previous work [Faccio et al. (2022), J. Mol. Liq. , 355 , 118922], we showed that it is possible to identify these two forms in water networks through a computational approach based on molecular dynamics simulation and on the calculation of the total communicability of the associated graph, in which the nodes correspond to water molecules and the edges represent the connections (interactions) between molecules. In this article, we present a more in-depth investigation of the application of graph-theory based approaches to the analysis of the structure of water networks. In particular, we investigate different connectivity and centrality measures and we report on the use of a variety of global metrics aimed at giving a topological and geometrical characterization of liquid water. [ABSTRACT FROM AUTHOR]