A singular thermodynamically consistent temperature at the origin of the anomalous behavior of liquid water

Do Specific Ion Effects on Collective Relaxation Arise from Perturbation of Hydrogen-Bonding Network Structure?

The change in the transport properties (i.e., water diffusivity, shear viscosity, etc.) when adding salts to water has been used to classify ions as either being chaotropic or kosmotropic, a terminology based on the presumption that this phenomenon arises from respective breakdown or enhancement of the hydrogen-bonding network structure. Recent quasi-elastic neutron scattering measurements of the collective structural relaxation time, τC, in aqueous salt solutions were interpreted as confirming this proposed origin of ion effects on the dynamics of water. However, we find similar changes in τC in the same salt solutions based on molecular dynamics (MD) simulations using a coarse-grained water model in which no hydrogen bonding exists, challenging this conventional interpretation of mobility change resulting from the addition of salts to water. A thorough understanding of specific ion effects should be useful in diverse material manufacturing and biomedical applications, where these effects are prevalent, but poorly understood.

Local and global expansivity in water

The supra-molecular structure of a liquid is strongly connected to its dynamics, which in turn control macroscopic properties such as viscosity. Consequently, detailed knowledge about how this structure changes with temperature is essential to understand the thermal evolution of the dynamics ranging from the liquid to the glass. Here, we combine infrared spectroscopy (IR) measurements of the hydrogen (H) bond stretching vibration of water with molecular dynamics simulations and employ a quantitative analysis to extract the inter-molecular H-bond length in a wide temperature range of the liquid. The extracted expansivity of this H-bond differs strongly from that of the average nearest neighbor distance of oxygen atoms obtained through a common conversion of mass density. However, both properties can be connected through a simple model based on a random loose packing of spheres with a variable coordination number, which demonstrates the relevance of supra-molecular arrangement. Furthermore, the exclusion of the expansivity of the inter-molecular H-bonds reveals that the most compact molecular arrangement is formed in the range of ∼316-331K (i.e., above the density maximum) close to the temperature of several pressure-related anomalies, which indicates a characteristic point in the supra-molecular arrangement. These results confirm our earlier approach to deduce inter-molecular H-bond lengths via IR in polyalcohols [Gabriel et al. J. Chem. Phys. 154, 024503 (2021)] quantitatively and open a new alley to investigate the role of inter-molecular expansion as a precursor of molecular fluctuations on a bond-specific level.

Prediction of water anomalous properties by introducing the two-state theory in SAFT

Water is one of the most abundant substances on earth, but it is still not entirely understood. It shows unusual behavior, and its properties present characteristic extrema unlike any other fluid. This unusual behavior has been linked to the two-state theory of water, which proposes that water forms different clusters, one with a high density and one with a low density, which may even form two distinct phases at low temperatures. Models incorporating the two-state theory manage to capture the unusual extrema of water, unlike traditional equations of state, which fail. In this work, we have derived the framework to incorporate the two-state theory of water into the Statistical-Associating-Fluid-Theory (SAFT). More specifically, we have assumed that water is an ideal solution of high density water molecules and low density water molecules that are in chemical equilibrium. Using this assumption, we have generalized the association term SAFT to allow for the simultaneous existence of the two water types, which have the same physical parameters but different association properties. We have incorporated the newly derived association term in the context of the Perturbed Chain-SAFT (PC-SAFT). The new model is referred to as PC-SAFT-Two-State (PC-SAFT-TS). Using PC-SAFT-TS, we have succeeded in predicting the characteristic extrema of water, such as its density and speed of sound maximum, etc., without loss of accuracy compared to the original PC-SAFT. This new framework is readily extended to mixtures, and PC-SAFT-TS manages to capture the solubility minimum of hydrocarbons in water in a straightforward manner.

Light-to-Heat Conversion of Optically Trapped Hot Brownian Particles

Anisotropic hybrid nanostructures stand out as promising therapeutic agents in photothermal conversion-based treatments. Accordingly, understanding local heat generation mediated by light-to-heat conversion of absorbing multicomponent nanoparticles at the single-particle level has forthwith become a subject of broad and current interest. Nonetheless, evaluating reliable temperature profiles around a single trapped nanoparticle is challenging from all of the experimental, computational, and fundamental viewpoints. Committed to filling this gap, the heat generation of an anisotropic hybrid nanostructure is explored by means of two different experimental approaches from which the local temperature is measured in a direct or indirect way, all in the context of hot Brownian motion theory. The results were compared with analytical results supported by the numerical computation of the wavelength-dependent absorption efficiencies in the discrete dipole approximation for scattering calculations, which has been extended to inhomogeneous nanostructures. Overall, we provide a consistent and comprehensive view of the heat generation in optical traps of highly absorbing particles from the viewpoint of the hot Brownian motion theory.

Liquid Water: A Single Approach to Its Two Continuous Phase Transitions

In this paper, we show that both continuous phase transitions of liquid water, the liquid-gas and the liquid-liquid, can be articulated within a single thermodynamic analytical formalism. This result follows from a combination of the two-liquid model (TLM), recently confirmed for water, with the idea of a thermal-dependent excluded volume, v_e, concept introduced by van der Waals, in his famous state equation. Starting from the fundamentals of thermodynamics, it will be shown that the TLM naturally leads to the idea of an extensive thermal-dependent v_e that acts as a parameter of the sample thermodynamic potentials. This procedure effectively separates the thermodynamics of the system into two parts: the first concerns the clusters' thermodynamics, taken as wandering particles, and the second concerns the thermal behavior of its internal structure (geometry and number of particles). From this result, we demonstrate that the condition of mechanical instability leads to not one but two critical points, each happening in one of the above-described parts of the system.

The Hydrophobic Effect Studied by Using Interacting Colloidal Suspensions

Interactions between nanoparticles (NPs) determine their self-organization and dynamic processes. In these systems, a quantitative description of the interparticle forces is complicated by the presence of the hydrophobic effect (HE), treatable only qualitatively, and due to the competition between the hydrophobic and hydrophilic forces. Recently, instead, a sort of crossover of HE from hydrophilic to hydrophobic has been experimentally observed on a local scale, by increasing the temperature, in pure confined water and studying the occurrence of this crossover in different water-methanol solutions. Starting from these results, we then considered the idea of studying this process in different nanoparticle solutions. By using photon correlation spectroscopy (PCS) experiments on dendrimer with OH terminal groups (dissolved in water and methanol, respectively), we show the existence of this hydrophobic-hydrophilic crossover with a well defined temperature and nanoparticle volume fraction dependence. In this frame, we have used the mode coupling theory extended model to evaluate the measured time-dependent density correlation functions (ISFs). In this context we will, therefore, show how the measured spectra are strongly dependent on the specificity of the interactions between the particles in solution. The observed transition demonstrates that just the HE, depending sensitively on the system thermodynamics, determines the hydrophobic and hydrophilic interaction properties of the studied nanostructures surface.

Exploring the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes

The breakdown of Stokes-Einstein relation in liquid water is one of the many anomalies that take place upon cooling and indicates the decoupling of diffusion and viscosity. It is hypothesized that these anomalies manifest due to the appearance of nanometer-scale spatial fluctuations, which become increasingly pronounced in the supercooled regime. Here, we explore the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes. We capture the diffusive dynamics of the probes using dynamic light scattering and target dynamics at different length scales by varying the probe size, from ≈100 nm silica spheres to molecular-sized polyhydroxylated fullerenes (≈1 nm). We find that all the studied probes, independent of size, display similar diffusive dynamics with an Arrhenius activation energy of ≈23 kJ mol^-1. Analysis of the diffusion coefficient further indicates that the probes, independent of their size, experience similar dynamic environment, which coincides with the macroscopic viscosity, while single water molecules effectively experience a comparatively lower viscosity. Finally, we conclude that our results indicate that the Stokes-Einstein relation is preserved for diffusion of probes in supercooled water T ≥ 260 K with size as small as ≈1 nm.

From Critical Point to Critical Point: The Two-States Model Describes Liquid Water Self-Diffusion from 623 to 126 K

Liquid’s behaviour, when close to critical points, is of extreme importance both for fundamental research and industrial applications. A detailed knowledge of the structural–dynamical correlations in their proximity is still today a target to reach. Liquid water anomalies are ascribed to the presence of a second liquid–liquid critical point, which seems to be located in the very deep supercooled regime, even below 200 K and at pressure around 2 kbar. In this work, the thermal behaviour of the self-diffusion coefficient for liquid water is analyzed, in terms of a two-states model, for the first time in a very wide thermal region (126 K < T < 623 K), including those of the two critical points. Further, the corresponding configurational entropy and isobaric-specific heat have been evaluated within the same interval. The two liquid states correspond to high and low-density water local structures that play a primary role on water dynamical behavior over 500 K.

Hydrophilic and Hydrophobic Effects on the Structure and Themodynamic Properties of Confined Water: Water in Solutions

NMR spectroscopy is used in the temperature range 180–350 K to study the local order and transport properties of pure liquid water (bulk and confined) and its solutions with glycerol and methanol at different molar fractions. We focused our interest on the hydrophobic effects (HE), i.e., the competition between hydrophilic and hydrophobic interactions. Nowadays, compared to hydrophilicity, little is known about hydrophobicity. Therefore, the main purpose of this study is to gain new information about hydrophobicity. As the liquid water properties are dominated by polymorphism (two coexisting liquid phases of high and low density) due to hydrogen bond interactions (HB), creating (especially in the supercooled regime) the tetrahedral networking, we focused our interest to the HE of these structures. We measured the relaxation times (T1 and T2) and the self-diffusion (DS). From these times, we took advantage of the NMR property to follow the behaviors of each molecular component (the hydrophilic and hydrophobic groups) separately. In contrast, DS is studied in terms of the Adam–Gibbs model by obtaining the configurational entropy (Sconf) and the specific heat contributions (CP,conf). We find that, for the HE, all of the studied quantities behave differently. For water–glycerol, the HB interaction is dominant for all conditions; water–methanol, two different T-regions above and below 265 K are observable, dominated by hydrophobicity and hydrophilicity, respectively. Below this temperature, where the LDL phase and the HB network develops and grows, with the times and CP,conf change behaviors leading to maxima and minima. Above it, the HB becomes weak and less stable, the HDL dominates, and hydrophobicity determines the solution.

The effect of multiple freeze-thaw cycles on the viscoelastic properties and microstructure of bovine superficial digital flexor tendon

The most common injuries of the human musculoskeletal system are related to soft tissue structures such as tendons or ligaments. To repair torn structures, surgical intervention and application of a biological or synthetic graft may be required. A typical procedure for the processing, storage, and distribution of soft tissue grafts involves at least two freezing/thawing (F/T) cycles. Even though repeated F/T cycles decrease the mechanical performance and change the structure of tendons, it is unclear whether there exists a maximum number of F/T cycles above which tendons should not be approved for use as a tissue allograft. To fill this research gap, we present an ex vivo study on the effects of repetitive F/T cycles on the biomechanical stability of bovine superficial digital flexor tendon tissue. Using mechanical testing supported with scanning electron microscopy imaging, we show that multiple F/T cycles affect the viscoelastic and structural properties of tissue by significantly reducing its tensile modulus after the 3^rd or 4^th F/T cycle (depending on the strain range), stress drop during relaxation after the 8^th F/T cycle (regardless the strain values), mechanical hysteresis after the 10^th F/T cycle, and by causing a significant decrease in collagen fibril diameter. Our results provide a deeper insight into understanding the mechanisms responsible for tissue damage during multiple F/T cycles, and thus, may be useful for the future optimization of tissue storage protocols.

Machine learning for reparameterization of four-site water models: TIP4P-BG and TIP4P-BGT

Parameterizing an effective water model is a challenging issue because of the difficulty in maintaining a comprehensive balance among the diverse physical properties of water with a limited number of parameters. The advancement in machine learning provides a promising path to search for a reliable set of parameters. Based on the TIP4P water model, hence, about 6000 molecular dynamics (MD) simulations for pure water at 1 atm and in the range of 273-373 K are conducted here as the training data. The back-propagation (BP) neural network is then utilized to construct an efficient mapping between the model parameters and four crucial physical properties of water, including the density, vaporization enthalpy, self-diffusion coefficient and viscosity. Without additional time-consuming MD simulations, this mapping operation could result in sufficient and accurate data for high-population genetic algorithm (GA) to optimize the model parameters as much as possible. Based on the proposed parameterizing strategy, TIP4P-BG (a conventional four-site water model) and TIP4P-BGT (an advanced model with temperature-dependent parameters) are established. Both the water models exhibit excellent performance with a reasonable balance among the four crucial physical properties. The relevant mean absolute percentage errors are 3.53% and 3.08%, respectively. Further calculations on the temperature of maximum density, isothermal compressibility, thermal expansion coefficient, radial distribution function and surface tension are also performed and the resulting values are in good agreement with the experimental values. Through this water modeling example, the potential of the proposed data-driven machine learning procedure has been demonstrated for parameterizing a MD-based material model.

Liquid dibromomethane under pressure: a computational study

Molecular dynamics simulations have been performed on liquid dibromomethane at thermodynamic states corresponding to temperature in the range 268-328 K and pressure varying from 1 bar to 3000 bar. The interaction model is a simple effective two-body pair potential with atom-atom Coulomb and Lennard-Jones interactions and molecules are rigid. Thermodynamic properties have been studied, including the isobaric thermal expansion coefficient, the isothermal compressibility, the heat capacities and the speed of sound. The simulation results exhibit a crossing of the isotherms of the isobaric thermal expansion coefficient at about 800 bar in very good agreement with the prediction of an isothermal fluctuation equation of state predicting such a crossing in the pressure range 650-900 bar, though experimental results up to 1000 bar do not find any crossing.

Some Aspects of the Liquid Water Thermodynamic Behavior: From The Stable to the Deep Supercooled Regime

Liquid water is considered to be a peculiar example of glass forming materials because of the possibility of giving rise to amorphous phases with different densities and of the thermodynamic anomalies that characterize its supercooled liquid phase. In the present work, literature data on the density of bulk liquid water are analyzed in a wide temperature-pressure range, also including the glass phases. A careful data analysis, which was performed on different density isobars, made in terms of thermodynamic response functions, like the thermal expansion αP and the specific heat differences CP−CV, proves, exclusively from the experimental data, the thermodynamic consistence of the liquid-liquid transition hypothesis. The study confirms that supercooled bulk water is a mixture of two liquid “phases”, namely the high density (HDL) and the low density (LDL) liquids that characterize different regions of the water phase diagram. Furthermore, the CP−CV isobars behaviors clearly support the existence of both a liquid–liquid transition and of a liquid–liquid critical point.

Thermal perturbation of NMR properties in small polar and non-polar molecules

Water is an important constituent in an abundant number of chemical systems; however, its presence complicates the analysis of in situ¹H MAS NMR investigations due to water’s ease of solidification and vaporization, the large changes in mobility, affinity for hydrogen bonding interactions, etc., that are reflected by dramatic changes in temperature-dependent chemical shielding. To understand the evolution of the signatures of water and other small molecules in complex environments, this work explores the thermally-perturbed NMR properties of water in detail by in situ MAS NMR over a wide temperature range. Our results substantially extend the previously published temperature-dependent ¹H and ¹⁷O chemical shifts, linewidths, and spin-lattice relaxation times over a much wider range of temperatures and with significantly enhanced thermal resolution. The following major results are obtained: Hydrogen bonding is clearly shown to weaken at elevated temperatures in both ¹H and ¹⁷O spectra, reflected by an increase in chemical shielding. At low temperatures, transient tetrahedral domains of H-bonding networks are evidenced and the observation of the transition between solid ice and liquid is made with quantitative considerations to the phase change. The ¹H chemical shift properties in other small polar and non-polar molecules have also been described over a range of temperatures, showing the dramatic effect hydrogen bonding perturbation on polar species. Gas phase species are observed and chemical exchange between gas and liquid phases is shown to play an important role on the observed NMR shifts. The results disclosed herein lay the foundation for a clear interpretation of complex systems during the increasingly popular in situ NMR characterization at elevated temperatures and pressures for studying chemical systems.

The Proton Density of States in Confined Water (H2O)

The hydrogen density of states (DOS) in confined water has been probed by inelastic neutron scattering spectra in a wide range of its P–T phase diagram. The liquid–liquid transition and the dynamical crossover from the fragile (super-Arrhenius) to strong (Arrhenius) glass forming behavior have been studied, by taking into account the system polymorphism in both the liquid and amorphous solid phases. The interest is focused in the low energy region of the DOS (E<10 meV) and the data are discussed in terms of the energy landscape (local minima of the potential energy) approach. In this latest research, we consider a unit scale energy (EC) linked to the water local order governed by the hydrogen bonding (HB). All the measured spectra, scaled according to such energy, evidence a universal power law behavior with different exponents (γ) in the strong and fragile glass forming regions, respectively. In the first case, the DOS data obey the Debye squared-frequency law, whereas, in the second one, we obtain a value predicted in terms of the mode-coupling theory (MCT) (γ≃1.6).

Some considerations on the water polymorphism and the liquid-liquid transition by the density behavior in the liquid phase

The bulk liquid water density data (ρ) are studied in a very large temperature pressure range including also the glass phases. A thorough analysis of their isobars, together with the suggestions of recent thermodynamical studies, gives evidence of two crossovers at T^* and P^* above which the hydrogen bond interaction is unable to arrange the tetrahedral network that is at the basis of the liquid polymorphism giving rise to the low density liquid (LDL). The curvatures of these isobars, as a function of T, are completely different: concave below P^* (where maxima are) and convex above. In both the cases, a continuity between liquid and glass is observed with P^* as the border of the density evolution toward the two different polymorphic glasses (low and high density amorphous). The experimental data of the densities of these two glasses also show a markedly different pressure dependence. Here, on the basis of these observations in bulk water and by considering a recent study on the growth of the LDL phase, by decreasing temperature, we discuss the water liquid-liquid transition and evaluate the isothermal compressibility inside the deep supercooled regime. Such a quantity shows an additional maximum that is pressure dependent that under ambient conditions agrees with a recent X-ray experiment. In particular, the present analysis suggests the presence of a liquid-liquid critical point located at about 180 MPa and 197 K.

In Vivo Water Dynamics in Shewanella oneidensis Bacteria at High Pressure

Following observations of survival of microbes and other life forms in deep subsurface environments it is necessary to understand their biological functioning under high pressure conditions. Key aspects of biochemical reactions and transport processes within cells are determined by the intracellular water dynamics. We studied water diffusion and rotational relaxation in live Shewanella oneidensis bacteria at pressures up to 500 MPa using quasi-elastic neutron scattering (QENS). The intracellular diffusion exhibits a significantly greater slowdown (by −10–30%) and an increase in rotational relaxation times (+10–40%) compared with water dynamics in the aqueous solutions used to resuspend the bacterial samples. Those results indicate both a pressure-induced viscosity increase and slowdown in ionic/macromolecular transport properties within the cells affecting the rates of metabolic and other biological processes. Our new data support emerging models for intracellular organisation with nanoscale water channels threading between macromolecular regions within a dynamically organized structure rather than a homogenous gel-like cytoplasm.

Pressure response of the THz spectrum of bulk liquid water revealed by intermolecular instantaneous normal mode analysis

The radial distribution functions of liquid water are known to change significantly their shape upon hydrostatic compression from ambient conditions deep into the kbar pressure regime. It has been shown that despite their eye-catching changes, the fundamental locally tetrahedral fourfold H-bonding pattern that characterizes ambient water is preserved up to about 10 kbar (1 GPa), which is the stability limit of liquid water at 300 K. The observed increase in coordination number comes from pushing water molecules into the first coordination sphere without establishing an H-bond, resulting in roughly two such additional interstitial molecules at 10 kbar. THz spectroscopy has been firmly established as a powerful experimental technique to analyze H-bonding in aqueous solutions given that it directly probes the far-infrared lineshape and thus the prominent H-bond network mode around 180 cm^-1. We, therefore, set out to assess pressure effects on the THz response of liquid water at 10 kbar in comparison to the 1 bar (0.1 MPa) reference, both at 300 K, with the aim to trace back the related lineshape changes to the structural level. To this end, we employ the instantaneous normal mode approximation to rigorously separate the H-bonding peak from the large background arising from the pronounced librational tail. By exactly decomposing the total molecular dynamics into hindered translations, hindered rotations, and intramolecular vibrations, we find that the H-bonding peak arises from translation-translation and translation-rotation correlations, which are successively decomposed down to the level of distinct local H-bond environments. Our utmost detailed analysis based on molecular pair classifications unveils that H-bonded double-donor water pairs contribute most to the THz response around 180 cm^-1, whereas interstitial waters are negligible. Moreover, short double-donor H-bonds have their peak maximum significantly shifted toward higher frequencies with respect to such long H-bonds. In conjunction with an increasing relative population of these short H-bonds versus the long ones (while the population of other water pair classes is essentially pressure insensitive), this explains not only the blue-shift of the H-bonding peak by about 20-30 cm^-1 in total from 1 bar to 10 kbar but also the filling of the shallow local minimum of the THz lineshape located in between the network peak and the red-wing of the librational band at 1 bar. Based on the changing populations as a function of pressure, we are also able to roughly estimate the pressure-dependence of the H-bond network mode and find that its pressure response and thus the blue-shifting are most pronounced at low kbar pressures.

The Role of Hydrogen Bonding in the Folding/Unfolding Process of Hydrated Lysozyme: A Review of Recent NMR and FTIR Results

The biological activity of proteins depends on their three-dimensional structure, known as the native state. The main force driving the correct folding mechanism is the hydrophobic effect and when this folding kinetics is altered, aggregation phenomena intervene causing the occurrence of illnesses such as Alzheimer and Parkinson’s diseases. The other important effect is performed by water molecules and by their ability to form a complex network of hydrogen bonds whose dynamics influence the mobility of protein amino acids. In this work, we review the recent results obtained by means of spectroscopic techniques, such as Fourier Transform Infrared (FTIR) and Nuclear Magnetic Resonance (NMR) spectroscopies, on hydrated lysozyme. In particular, we explore the Energy Landscape from the thermal region of configurational stability up to that of the irreversible denaturation. The importance of the coupling between the solute and the solvent will be highlighted as well as the different behaviors of hydrophilic and hydrophobic moieties of protein amino acid residues.

Activity of Supercooled Water on the Ice Curve and Other Thermodynamic Properties of Liquid Water up to the Boiling Point at Standard Pressure

A simple model for thermodynamic properties of water from subzero temperatures up to 373 K was derived at ambient pressure. The heat capacity of supercooled water was assessed as lambda transition. The obtained properties for supercooled water such as heat capacity, vapor pressure, density and thermal expansion are in excellent agreement with literature data. Activity of water on ice curve, independent of used electrolyte and Debye-Hückel constant applied in modeling, is also calculated. Thus, the ice curve activity of supercooled water can be used as a universal basis for thermodynamic modeling of aqueous solutions, precipitating hydrated and anhydrous solids. A simple model for heat capacity, density and thermal expansion of ice are also derived from 170 K up to melting point.

The evaluation of the hydrophilic-hydrophobic interactions and their effect in water-methanol solutions: A study in terms of the thermodynamic state functions in the frame of the transition state theory

Aqueous solutions of amphiphilic molecules are characterized by the competition between hydrophilic and hydrophobic interactions. These interactions have a different energetic dependence with the temperature. Whereas hydrophilic interactions have been well characterized, a complete theory for the hydrophobic ones is still lacking as well as the comprehension of the effect that the solvent exerts on the solute and vice versa. In this paper from the measured relaxation time, we evaluated the thermodynamic state functions of water-methanol solutions in the frame of the transition state theory. In particular we study the behavior of the Gibbs free energy, enthalpy and entropy of water, methanol and some of their solutions as a function of both temperature and water molar fraction. Our results indicate that the temperature of about 280 K represents a crossover between two regions dominated by hydrophobicity (high T) and hydrophilicity (low T).

Temperature-dependent optoacoustic response and transient through zero Grüneisen parameter in optically contrasted media

Non-invasive optoacoustic mapping of temperature in tissues with low blood content can be enabled by administering external contrast agents. Some important clinical applications of such approach include temperature mapping during thermal therapies in a prostate or a mammary gland. However, the technique would require a calibration that establishes functional relationship between the measured normalized optoacoustic response and local tissue temperature. In this work, we investigate how a key calibration parameter - the temperature of zero optoacoustic response (T₀ ) - behaves in different environments simulating biological tissues augmented with either dissolved or particulate (nanoparticles) contrast agents. The observed behavior of T₀ in ionic and molecular solutions suggests that in-vivo temperature mapping is feasible for contrast agents of this type, but requires knowledge of local concentrations. Oppositely, particulate contrast agents (plasmonic or carbon nanoparticles) demonstrated concentration-independent thermal behavior of optoacoustic response with T₀ defined by the thermoelastic properties of the local environment.

Unveiling Molecular Changes in Water by Small Luminescent Nanoparticles

Nowadays a large variety of applications are based on solid nanoparticles dispersed in liquids-so called nanofluids. The interaction between the fluid and the nanoparticles plays a decisive role in the physical properties of the nanofluid. A novel approach based on the nonradiative energy transfer between two small luminescent nanocrystals (GdVO₄ :Nd³⁺ and GdVO₄ :Yb³⁺ ) dispersed in water is used in this work to investigate how temperature affects both the processes of interaction between nanoparticles and the effect of the fluid on the nanoparticles. From a systematic analysis of the effect of temperature on the GdVO₄ :Nd³⁺ → GdVO₄ :Yb³⁺ interparticle energy transfer, it can be concluded that a dramatic increase in the energy transfer efficiency occurs for temperatures above 45 °C. This change is properly explained by taking into account a crossover existing in diverse water properties that occurs at about this temperature. The obtained results allow elucidation on the molecular arrangement of water molecules below and above this crossover temperature. In addition, it is observed that an energy transfer process is produced as a result of interparticle collisions that induce irreversible ion exchange between the interacting nanoparticles.

Dynamical properties of water-methanol solutions

We study the relaxation times tα in the water-methanol system. We examine new data and data from the literature in the large temperature range 163 < T < 335 K obtained using different experimental techniques and focus on how tα affects the hydrogen bond structure of the system and the hydrophobicity of the alcohol methyl group. We examine the relaxation times at a fixed temperature as a function of the water molar fraction XW and observe two opposite behaviors in their curvature when the system moves from high to low T regimes. This behavior differs from that of an ideal solution in that it has excess values located at different molar fractions (XW = 0.5 for high T and 0.75 in the deep supercooled regime). We analyze the data and find that above a crossover temperature T ∼ 223 K, hydrophobicity plays a significant role and below it the water tetrahedral network dominates. This temperature is coincident with the fragile-to-strong dynamical crossover observed in confined water and supports the liquid-liquid phase transition hypothesis. At the same time, the reported data suggest that this crossover temperature (identified as the Widom line temperature) also depends on the alcohol concentration.

The interplay between dynamic heterogeneities and structure of bulk liquid water: A molecular dynamics simulation study

In order to study the interplay between dynamical heterogeneities and structural properties of bulk liquid water in the temperature range 130-350 K, thus including the supercooled regime, we use the explicit trend of the distribution functions of some molecular properties, namely, the rotational relaxation constants, the atomic mean-square displacements, the relaxation of the cross correlation functions between the linear and squared displacements of H and O atoms of each molecule, the tetrahedral order parameter q and, finally, the number of nearest neighbors (NNs) and of hydrogen bonds (HBs) per molecule. Two different potentials are considered: TIP4P-Ew and a model developed in this laboratory for the study of nanoconfined water. The results are similar for the dynamical properties, but are markedly different for the structural characteristics. In particular, for temperatures higher than that of the dynamic crossover between "fragile" (at higher temperatures) and "strong" (at lower temperatures) liquid behaviors detected around 207 K, the rotational relaxation of supercooled water appears to be remarkably homogeneous. However, the structural parameters (number of NNs and of HBs, as well as q) do not show homogeneous distributions, and these distributions are different for the two water models. Another dynamic crossover between "fragile" (at lower temperatures) and "strong" (at higher temperatures) liquid behaviors, corresponding to the one found experimentally at T(∗) ∼ 315 ± 5 K, was spotted at T(∗) ∼ 283 K and T(∗) ∼ 276 K for the TIP4P-Ew and the model developed in this laboratory, respectively. It was detected from the trend of Arrhenius plots of dynamic quantities and from the onset of a further heterogeneity in the rotational relaxation. To our best knowledge, it is the first time that this dynamical crossover is detected in computer simulations of bulk water. On the basis of the simulation results, the possible mechanisms of the two crossovers at molecular level are discussed.

Energy landscape in protein folding and unfolding

We use (1)H NMR to probe the energy landscape in the protein folding and unfolding process. Using the scheme ⇄ reversible unfolded (intermediate) → irreversible unfolded (denatured) state, we study the thermal denaturation of hydrated lysozyme that occurs when the temperature is increased. Using thermal cycles in the range 295 < T < 365 K and following different trajectories along the protein energy surface, we observe that the hydrophilic (the amide NH) and hydrophobic (methyl CH3 and methine CH) peptide groups evolve and exhibit different behaviors. We also discuss the role of water and hydrogen bonding in the protein configurational stability.

Understanding n-Octane Behavior near Graphene with Scaled Solvent-Solute Attractions

We employ molecular dynamics simulations of n-octane near a layered graphene surface to study the related phenomena of solvation, density fluctuations, wettability, and structure and dynamics of n-octane molecules in the inhomogeneous interfacial environment. That solvation in bulk n-octane displays a lengthscale-dependent crossover similar to that of hydrophobic solvation in water is known. Here we show that, near an extended graphene interface having attractive interactions with n-octane, lengthscale-dependent solvation is similar to that in the bulk and displays a small to large crossover. However, as the n-octane-graphene interactions are reduced to make the surface increasingly solvophobic, the crossover behavior is modulated and essentially absent near the most solvophobic surfaces, similar to that in water near hydrophobic interfaces. We show that the macroscopic measure of wettability, namely, the contact angle, characterizes n-octane-graphene coupling over a limited range of attractions. In contrast, molecular measures such as the free energy of cavity formation or the local compressibility in the interfacial region provide an effective measure of this coupling over a broader range of attractions. Finally, as n-octane-graphene attractions are increased, the n-octane liquid displays a wetting transition and corresponding change from sigmoidal to layered density profile. Analysis of the local structure shows that n-octane molecules prefer approximately linear conformations and surface-parallel orientations near the graphene surface, and their translational dynamics slow down with increasing n-octane-graphene attractions. Our study highlights molecular scale behavior of n-octane molecules that is relevant to understanding nanoparticle-solvent coupling in composite materials with enhanced mechanical or thermal properties.

Some considerations on the transport properties of water-glycerol suspensions

We study the self-diffusion coefficient and viscosity of a water-glycerol mixture for several glycerol molar fractions as a function of temperature well inside the metastable supercooled regime. We perform NMR experiments and verify that the system has at different concentration a fragile-to-strong crossover accompanied by the violation of the Stokes-Einstein relation. We observe that the crossover temperature depends on the water amount. Studying the fractional representation of the Stokes-Einstein relation, we find that in these systems dynamical arrest does not exhibit criticality and the transport parameters have a universal behavior.

Density and structural anomalies in soft-repulsive dimeric fluids

We performed a simulation study of the fluid structure of dimeric particles interacting via a core-softened potential and shed light on their anomalous behaviours upon varying both geometrical and interaction parameters.

The structure of water around the compressibility minimum

Here we present diffraction data that yield the oxygen-oxygen pair distribution function, g(OO)(r) over the range 254.2-365.9 K. The running O-O coordination number, which represents the integral of the pair distribution function as a function of radial distance, is found to exhibit an isosbestic point at 3.30(5) Å. The probability of finding an oxygen atom surrounding another oxygen at this distance is therefore shown to be independent of temperature and corresponds to an O-O coordination number of 4.3(2). Moreover, the experimental data also show a continuous transition associated with the second peak position in g(OO)(r) concomitant with the compressibility minimum at 319 K.

Thermodynamic properties of bulk and confined water

The thermodynamic response functions of water display anomalous behaviors. We study these anomalous behaviors in bulk and confined water. We use nuclear magnetic resonance (NMR) to examine the configurational specific heat and the transport parameters in both the thermal stable and the metastable supercooled phases. The data we obtain suggest that there is a behavior common to both phases: that the dynamics of water exhibit two singular temperatures belonging to the supercooled and the stable phase, respectively. One is the dynamic fragile-to-strong crossover temperature (T(L) ≃ 225 K). The second, T* ∼ 315 ± 5 K, is a special locus of the isothermal compressibility K(T)(T, P) and the thermal expansion coefficient α(P)(T, P) in the P-T plane. In the case of water confined inside a protein, we observe that these two temperatures mark, respectively, the onset of protein flexibility from its low temperature glass state (T(L)) and the onset of the unfolding process (T*).

The thermodynamical response functions and the origin of the anomalous behavior of liquid water

The density maximum of water dominates the thermodynamics of the system under ambient conditions, is strongly P-dependent, and disappears at a crossover pressure P(cross) approximately 1.8 kbar. We study this variable across a wide area of the T-P phase diagram. We consider old and new data of both the isothermal compressibility K(T)(T, P), the pressure constant specific heat C(P)(T) and the coefficient of thermal expansion alpha(P) (T, P). We observe that K(T)(T) shows a minimum at T* approximately 315 +/- 5 K for all of the studied pressures, whereas, at the same temperature, C(P)(T) has the minimal variation as a function of P in the interval 1 bar-4 kbar. We find the behavior of alpha(P) also to be surprising: all the alpha(P)(T) curves measured at different P cross at T*. The experimental data show a "singular and universal expansivity point" at T* approximately 315 K and alpha(P)(T*) = 0.44 10(-3) K(-1). Unlike other water singularities, we find this temperature to be thermodynamically consistent in the relationship connecting the three response functions. By considering also the P-T behavior of the self-diffusion coefficient D(S) and of the NMR proton chemical shift delta we have the information that at T* the water local order points out, with decreasing T, the crossover from a normal fluid to the anomalous and complex liquid characterized by the many anomalies.