Beyond Local Structures in Critical Supercooled Water through Unsupervised Learning.

The presence of a second critical point in water has been a topic of intense investigation for the last few decades. The molecular origins underlying this phenomenon are typically rationalized in terms of the competition between local high-density (HD) and low-density (LD) structures. Their identification often requires designing parameters that are subject to human intervention. Herein, we use unsupervised learning to discover structures in atomistic simulations of water close to the liquid-liquid critical point (LLCP). Encoding the information on the environment using local descriptors, we do not find evidence for two distinct thermodynamic structures. In contrast, when we deploy nonlocal descriptors that probe instead heterogeneities on the nanometer length scale, this leads to the emergence of LD and HD domains rationalizing the microscopic origins of the density fluctuations close to criticality.

Exploring the Mechanisms behind Non-aromatic Fluorescence with the Density Functional Tight Binding Method.

Recent experimental findings reveal nonconventional fluorescence emission in biological systems devoid of conjugated bonds or aromatic compounds, termed non-aromatic fluorescence (NAF). This phenomenon is exclusive to aggregated or solid states and remains absent in monomeric solutions. Previous studies focused on small model systems in vacuum show that the carbonyl stretching mode along with strong interaction of short hydrogen bonds (SHBs) remains the primary vibrational mode explaining NAF in these systems. In order to simulate larger model systems taking into account the effects of the surrounding environment, in this work we propose using the density functional tight-binding (DFTB) method in combination with non-adiabatic molecular dynamics (NAMD) and the mixed quantum/molecular mechanics (QM/MM) approach. We investigate the mechanism behind NAF in the crystal structure of l-pyroglutamine-ammonium, comparing it with the related nonfluorescent amino acid l-glutamine. Our results extend our previous findings to more realistic systems, demonstrating the efficiency and robustness of the proposed DFTB method in the context of NAMD in biological systems. Furthermore, due to its inherent low computational cost, this method allows for a better sampling of the nonradiative events at the conical intersection which is crucial for a complete understanding of this phenomenon. Beyond contributing to the ongoing exploration of NAF, this work paves the way for future application of this method in more complex biological systems such as amyloid aggregates, biomaterials, and non-aromatic proteins.

Photoinactivation of Enterococcus faecalis Biofilm: In Vitro Antimicrobial Effect of Photoexcited Rutin-Gallium(III) Complex via Visible Blue Light.

INTRODUCTION: Endodontic infection is a common problem that can result in tooth loss if not effectively treated. This study focused on investigating the use of rutin-gallium (Ga)(III) complex-mediated antimicrobial photodynamic therapy (aPDT) for the photoinactivation of Enterococcus faecalis biofilm. METHODS: The minimum biofilm eradication concentration of the rutin-Ga(III) complex and the minimum biofilm eradication dose of light-emitting diode against E. faecalis were evaluated. The antimicrobial effect of rutin-Ga(III) complex-mediated aPDT against E. faecalis was assessed. Additionally, the expression of genes associated with E. faecalis virulence, such as ace, gelE, and esp, as well as the production of reactive oxygen species within the cells were evaluated. RESULTS: The minimum biofilm eradication concentration of the rutin-Ga(III) complex was determined to be 25 µmol/L, whereas the minimum biofilm eradication dose of light-emitting diode irradiation was defined as 5 minutes with an energy density of 300-420 J/cm2. Rutin-Ga(III) complex-mediated aPDT demonstrated a significant dose-dependent reduction in the growth of E. faecalis biofilms. Moreover, aPDT led to increased intracellular reactive oxygen species generation in treated E. faecalis cells. Furthermore, the messenger RNA levels of ace, gelE, and esp genes were significantly down-regulated in E. faecalis treated with rutin-Ga(III) complex-mediated aPDT (P < .05). CONCLUSIONS: Rutin-Ga(III) complex-mediated aPDT effectively reduces E. faecalis biofilm growth by disrupting biofilm structure and down-regulating virulence genes. These findings highlight the potential of aPDT with the rutin-Ga(III) complex as an adjuvant therapeutic approach against E. faecalis biofilms.

Lessons Learned from Multiobjective Automatic Optimizations of Classical Three-Site Rigid Water Models Using Microscopic and Macroscopic Target Experimental Observables.

The development of accurate water models is of primary importance for molecular simulations. Despite their intrinsic approximations, three-site rigid water models are still ubiquitously used to simulate a variety of molecular systems. Automatic optimization approaches have been recently used to iteratively refine three-site water models to fit macroscopic (average) thermodynamic properties, providing state-of-the-art three-site models that still present some deviations from the liquid water properties. Here, we show the results obtained by automatically optimizing three-site rigid water models to fit a combination of microscopic and macroscopic experimental observables. We use Swarm-CG, a multiobjective particle-swarm-optimization algorithm, for training the models to reproduce the experimental radial distribution functions of liquid water at various temperatures (rich in microscopic-level information on, e.g., the local orientation and interactions of the water molecules). We systematically analyze the agreement of these models with experimental observables and the effect of adding macroscopic information to the training set. Our results demonstrate how adding microscopic-rich information in the training of water models allows one to achieve state-of-the-art accuracy in an efficient way. Limitations in the approach and in the approximated description of water in these three-site models are also discussed, providing a demonstrative case useful for the optimization of approximated molecular models, in general.

The carbonyl-lock mechanism underlying non-aromatic fluorescence in biological matter.

Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics.

Evaluation of anti-biofilm effect of antimicrobial sonodynamic therapy-based periodontal ligament stem cell-derived exosome-loaded kojic acid on Enterococcus faecalis biofilm.

Introduction. Antimicrobial sonodynamic therapy (aSDT) is an approach that uses ultrasound waves (UWs) and a sonosensitizer to generate reactive oxygen species (ROS) to damage microbial cells in biofilms. Using nano-carriers, such as exosomes (Exos), to deliver the sonosensitizer can potentially enhance the effectiveness of aSDT.Hypothesis/Gap Statement. aSDT can downregulate the expression of gelE and sprE genes, increasing the production of endogenous ROS and degradation of pre-formed Enterococcus faecalis biofilms.Aim. This study investigated the anti-biofilm effect of aSDT-based periodontal ligament stem cell-derived exosome-loaded kojic acid (KA@PDL-Exo) on pre-formed E. faecalis biofilms in root canals.Methodology. Following the isolation and characterization of PDL-Exo, KA@PDL-Exo was prepared and confirmed. The minimal biofilm inhibitory concentration (MBIC) of KA, PDL-Exo, KA@PDL-Exo and sodium hypochlorite (NaOCl) was determined, and their anti-biofilm effects were assessed with and without UWs. The binding affinity of KA with GelE and SprE proteins was evaluated using in silico molecular docking. Additionally, the study measured the generation of endogenous ROS and evaluated changes in the gene expression levels of gelE and sprE.Results. The results revealed a dose-dependent decrease in the viability of E. faecalis cells within biofilms. KA@PDL-Exo was the most effective, with an MBIC of 62.5 µg ml-1, while NaOCl, KA and PDL-Exo had MBIC values of 125, 250 and 500 µg ml-1, respectively. The use of KA@PDL-Exo-mediated aSDT resulted in a significant reduction of the E. faecalis biofilm (3.22±0.36 log10 c.f.u. ml-1; P<0.05). The molecular docking analysis revealed docking scores of -5.3 and -5.2 kcal mol-1 for GelE-KA an SprE-KA, respectively. The findings observed the most significant reduction in gene expression of gelE and sprE in the KA@PDL-Exo group, with a decrease of 7.9- and 9.3-fold, respectively, compared to the control group (P<0.05).Conclusion. The KA@PDL-Exo-mediated aSDT was able to significantly reduce the E. faecalis load in pre-formed biofilms, decrease the expression of gelE and srpE mRNA, and increase the generation of endogenous ROS. These findings imply that KA@PDL-Exo-mediated aSDT could be a promising anti-biofilm strategy that requires additional in vitro and in vivo investigations.

ZundEig: The Structure of the Proton in Liquid Water from Unsupervised Learning.

The structure of the excess proton in liquid water has been the subject of lively debate on both experimental and theoretical fronts for the last century. Fluctuations of the proton are typically interpreted in terms of limiting states referred to as the Eigen and Zundel species. Here, we put these ideas under the microscope, taking advantage of recent advances in unsupervised learning that use local atomic descriptors to characterize environments of acidic water combined with advanced clustering techniques. Our agnostic approach leads to the observation of only one charged cluster and two neutral ones. We demonstrate that the charged cluster involving the excess proton is best seen as an ionic topological defect in water's hydrogen bond network, forming a single local minimum on the global free-energy landscape. This charged defect is a highly fluxional moiety, where the idealized Eigen and Zundel species are neither limiting configurations nor distinct thermodynamic states. Instead, the ionic defect enhances the presence of neutral water defects through strong interactions with the network. We dub the combination of the charged and neutral defect clusters as ZundEig, demonstrating that the fluctuations between these local environments provide a general framework for rationalizing more descriptive notions of the proton in the existing literature.

Zirconium Coordination Chemistry and Its Role in Optimizing Hydroxymate Chelation: Insights from Molecular Dynamics.

In the past decade, there has been a growth in using Zirconium-89 (89Zr) as a radionuclide in nuclear medicine for cancer diagnostic imaging and drug discovery processes. Although one of the most popular chelators for 89Zr, desferrioxamine (DFO) is typically presented as a hexadentate ligand, our work suggests a different scenario. The coordination structure of the Zr4+-DFO complex has primarily been informed by DFT-based calculations, which typically ignore temperature and therefore entropic and dynamic solvent effects. In this work, free energy calculations using molecular dynamics simulations, where the conformational fluctuations of both the ligand and the solvent are explicitly included, are used to compare the binding of Zr4+ cations with two different chelators, DFO and 4HMS, the latter of which is an octadentate ligand that has been recently proposed as a better chelator due to the presence of four hydroxymate groups. We find that thermally induced disorder leads to an open hexadentate chelate structure of the Zr4+-DFO complex, leaving the Zr4+ metal exposed to the solvent. A stable coordination of Zr4+ with 4HMS, however, is formed by involving both hydroxamate groups and water molecules in a more closely packed structure.

Free energy decompositions illuminate synergistic effects in interfacial binding thermodynamics of mixed surfactant systems.

Recent vibrational sum frequency generation spectroscopic experiments [Sengupta et al., J. Phys. Chem. Lett. 13, 11391-11397 (2022)] demonstrated synergistic interfacial adsorption effects between the anionic dodecyl sulfate (DS-) and the polar, but charge-neutral hexaethylene glycol monododecyl ether (C12E6), surfactants. In this study, the interfacial adsorption thermodynamics underlying these synergistic effects are analyzed through free energy decompositions. A general decomposition method utilizing alchemical intermediate states is outlined. Combining free energy decompositions with the potential distribution theorem illuminates the statistical interpretations of correlated effects between different system components. This approach allows for the identification of the physical effects leading to synergistic adsorption thermodynamics of DS- binding to the air-C12E6-water interface. The binding properties are found to result from a combination of effects predominantly including energetic van der Waals stabilization between DS- and C12E6, as well as competing energetic and entropic effects due to changes in the interfacial water structure as a result of introducing a C12E6 monolayer into the bare air-water interface.

Examining the origins of observed terahertz modes from an optically pumped atomistic model protein in aqueous solution.

The microscopic origins of terahertz (THz) vibrational modes in biological systems are an active and open area of current research. Recent experiments [Phys Rev X. 8, 031061 (2018)] have revealed the presence of a pronounced mode at ∼0.3 THz in fluorophore-decorated bovine serum albumin (BSA) protein in aqueous solution under nonequilibrium conditions induced by optical pumping. This result was heuristically interpreted as a collective elastic fluctuation originating from the activation of a low-frequency phonon mode. In this work, we show that the sub-THz spectroscopic response emerges in a statistically significant manner (>2σ) from such collective behavior, illustrating how photoexcitation can alter specific THz vibrational modes. We revisit the theoretical analysis with proof-of-concept molecular dynamics that introduce optical excitations into the simulations. Using information theory techniques, we show that these excitations can give rise to a multiscale response involving two optically excited chromophores (tryptophans), other amino acids in the protein, ions, and water. Our results motivate new experiments and fully nonequilibrium simulations to probe these phenomena, as well as the refinement of atomistic models of Fröhlich condensates that are fundamentally determined by nonlinear interactions in biology.

Solvation thermodynamics from cavity shapes of amino acids.

According to common physical chemistry wisdom, the solvent cavities hosting a solute are tightly sewn around it, practically coinciding with its van der Waals surface. Solvation entropy is primarily determined by the surface and the volume of the cavity while enthalpy is determined by the solute-solvent interaction. In this work, we challenge this picture, demonstrating by molecular dynamics simulations that the cavities surrounding the 20 amino acids deviate significantly from the molecular surface. Strikingly, the shape of the cavity alone can be used to predict the solvation free energy, entropy, enthalpy, and hydrophobicity. Solute-solvent interactions involving the different chemical moieties of the amino acid, determine indirectly the cavity shape, and the properties of the branches but do not have to be taken explicitly into account in the prediction model.

Effect of Quantum Delocalization on Temperature Dependent Double Proton Transfer in Molecular Crystals of Terephthalic Acid.

Double proton transfers (DPTs) are important for several physical processes, both in molecules and in the condensed phase. While these have been widely studied in biological systems, their study in crystalline environments is rare. In this work, using path integral molecular dynamics simulations, we have studied temperature dependent DPT in molecular crystals of terephthalic acid (TPA). In accordance with experimental reports, we find evidence for a double proton transfer induced order-to-disorder transition that is sensitive to the inclusion of nuclear quantum effects. Our simulations show that in addition to the presence of L and R tautomers of terepthalic acid, there are a small but non-negligible concentration of positive and negatively charged pairs of TPA molecules. At the onset of the transition at low temperatures, DPT likely occurs through a tunneling mechanism while at room temperature, likely involving the dominance of activated hopping. Through an analysis of the electronic structure of the system using Wannier functions, we show that the H atom shuttling between the donor and acceptor O atoms involves a proton.

Water Self-Dissociation is Insensitive to Nanoscale Environments.

Nanoconfinement effects on water dissociation and reactivity remain controversial, despite their importance to understand the aqueous chemistry at interfaces, pores, or aerosols. The pKw in confined environments has been assessed from experiments and simulations in a few specific cases, leading to dissimilar conclusions. Here, with the use of carefully designed ab initio simulations, we demonstrate that the energetics of bulk water dissociation is conserved intact to unexpectedly small length-scales, down to aggregates of only a dozen molecules or pores of widths below 2 nm. The reason is that most of the free-energy involved in water autoionization comes from breaking the O-H covalent bond, which has a comparable barrier in the bulk liquid, in a small droplet of nanometer size, or in a nanopore in the absence of strong interfacial interactions. Thus, dissociation free-energy profiles in nanoscopic aggregates or in 2D slabs of 1 nm width reproduce the behavior corresponding to the bulk liquid, regardless of whether the corresponding nanophase is delimited by a solid or a gas interface. The present work provides a definite and fundamental description of the mechanism and thermodynamics of water dissociation at different scales with broader implications on reactivity and self-ionization at the air-liquid interface.

Do Machine-Learning Atomic Descriptors and Order Parameters Tell the Same Story? The Case of Liquid Water.

Machine-learning (ML) has become a key workhorse in molecular simulations. Building an ML model in this context involves encoding the information on chemical environments using local atomic descriptors. In this work, we focus on the Smooth Overlap of Atomic Positions (SOAP) and their application in studying the properties of liquid water both in the bulk and at the hydrophobic air-water interface. By using a statistical test aimed at assessing the relative information content of different distance measures defined on the same data space, we investigate if these descriptors provide the same information as some of the common order parameters that are used to characterize local water structure such as hydrogen bonding, density, or tetrahedrality to name a few. Our analysis suggests that the ML description and the standard order parameters of the local water structure are not equivalent. In particular, a combination of these order parameters probing local water environments can predict SOAP similarity only approximately, and vice versa, the environments that are similar according to SOAP are not necessarily similar according to the standard order parameters. We also elucidate the role of some of the metaparameters in the SOAP definition in encoding chemical information.

Nonreciprocal forces enable cold-to-hot heat transfer between nanoparticles.

We study the heat transfer between two nanoparticles held at different temperatures that interact through nonreciprocal forces, by combining molecular dynamics simulations with stochastic thermodynamics. Our simulations reveal that it is possible to construct nano refrigerators that generate a net heat transfer from a cold to a hot reservoir at the expense of power exerted by the nonreciprocal forces. Applying concepts from stochastic thermodynamics to a minimal underdamped Langevin model, we derive exact analytical expressions predictions for the fluctuations of work, heat, and efficiency, which reproduce thermodynamic quantities extracted from the molecular dynamics simulations. The theory only involves a single unknown parameter, namely an effective friction coefficient, which we estimate fitting the results of the molecular dynamics simulation to our theoretical predictions. Using this framework, we also establish design principles which identify the minimal amount of entropy production that is needed to achieve a certain amount of uncertainty in the power fluctuations of our nano refrigerator. Taken together, our results shed light on how the direction and fluctuations of heat flows in natural and artificial nano machines can be accurately quantified and controlled by using nonreciprocal forces.

The collective burst mechanism of angular jumps in liquid water.

Understanding the microscopic origins of collective reorientational motions in aqueous systems requires techniques that allow us to reach beyond our chemical imagination. Herein, we elucidate a mechanism using a protocol that automatically detects abrupt motions in reorientational dynamics, showing that large angular jumps in liquid water involve highly cooperative orchestrated motions. Our automatized detection of angular fluctuations, unravels a heterogeneity in the type of angular jumps occurring concertedly in the system. We show that large orientational motions require a highly collective dynamical process involving correlated motion of many water molecules in the hydrogen-bond network that form spatially connected clusters going beyond the local angular jump mechanism. This phenomenon is rooted in the collective fluctuations of the network topology which results in the creation of defects in waves on the THz timescale. The mechanism we propose involves a cascade of hydrogen-bond fluctuations underlying angular jumps and provides new insights into the current localized picture of angular jumps, and its wide use in the interpretations of numerous spectroscopies as well in reorientational dynamics of water near biological and inorganic systems. The role of finite size effects, as well as of the chosen water model, on the collective reorientation is also elucidated.

Ex-vivo effects of propolis quantum dots-nisin-nanoquercetin-mediated photodynamic therapy on Streptococcus mutans biofilms and white spot lesions.

BACKGROUND: White spot lesions (WSLs) remain one of the most critical adverse sequelae of fixed orthodontic treatment, despite materials and techniques advances in orthodontics. WSLs seem to be a multi-factorial interaction including increased microbial plaque due to intrabuccal appliances that limit the oral-cleansing mechanism and change in the oral microbiome during fixed appliance wear. The aim of this study was to investigate the synergistic effect of propolis quantum dots (PQD), nisin (Nis), and quercetin nanoparticles (nQCT)-mediated photodynamic therapy (PQD-Nis-nQCT-mediated aPDT) in the eradication of Streptococcus mutans biofilms and the remineralization of WSLs ex-vivo. MATERIALS AND METHODS: The cytotoxicity of PQD-Nis-nQCT composite on human gingival fibroblasts was evaluated using neutral red. Intracellular reactive oxygen species (ROS) generation following PQD-Nis-nQCT-mediated aPDT was measured. Enamel slabs were prepared and demineralized using a demineralization solution containing S. mutans. Demineralized enamel slabs were divided into 9 groups (n = 10) and treated in the following groups: 1) Artificial saliva (negative control), 2) 2% neutral sodium fluoride gel (NSF; positive control or treatment control, 3) PQD, 4) Nis, 5) nQCT, 6) Nis-nQCT, 7) PQD-Nis-nQCT 8) Blue laser irradiation (light), 9) PQD-Nis-nQCT with irradiation (PQD-Nis-nQCT-mediated aPDT). Then, the surface changes, microhardness, and surface topography of the demineralized slabs were examined following each treatment using DIAGNOdent Pen reading, digital hardness tester, and SEM, respectively. After the determination of minimum biofilm eradication concentration (MBEC) of PQD, Nis, and nQCT by microtiter plate assay, the synergistic antimicrobial effects of PQD and Nis-nQCT were determined via evaluation of fractional biofilm eradication concentration (FBEC) index. The anti-biofilm effects of each treatment on S. mutans were assessed using a colorimetric assay. The virulence­associated gtfB gene expression was assessed following PQD-Nis-nQCT-mediated aPDT by quantitative real­time PCR. RESULTS: PQD-Nis-nQCT at 2048 µg/mL had no significant cell cytotoxicity on human gingival fibroblasts compared to the control group (P > 0.05). A significantly increased (7.6 fold) in intracellular ROS was observed following PQD-Nis-nQCT-mediated aPDT (13.9 ± 1.41) when compared to the control (1.83 ± 0.13). Following each treatment, the microhardness of the demineralized enamel surface significantly increased except for the artificial saliva (negative) and blue laser irradiation groups. The highest change in microhardness improvement was detected in the PQD-Nis-nQCT-mediated aPDT group (P < 0.05). Also, DIAGNODent Pen reading revealed the highest significant improved change in the level of mineralization degree in the PQD-Nis-nQCT-mediated aPDT group. Nis and blue light irradiation groups, like the artificial saliva-treated demineralized enamel slabs (control group), did not lead to remineralization (P > 0.05). Also, the PQD-Nis-nQCT-mediated aPDT treatment results obtained from SEM revealed that remineralization of demineralized enamel slabs in that group has significantly improved compared to the others. Light-activated nQCT, PQD, Nis-nQCT, and PQD-Nis-nQCT composite significantly reduced pre-formed biofilms of S. mutans compared with unactivated forms of test materials. The relative expression level of the virulence gtfB gene was significantly decreased (7.53-fold) in the presence of PQD-Nis-nQCT-mediated aPDT (P < 0.05). CONCLUSION: PQD-Nis-nQCT-mediated aPDT can be used for the eradication of S. mutans biofilms and remineralization of WSLs. The found in vitro efficacy should be tested further through clinical studies.

Observation of Strong Synergy in the Interfacial Water Response of Binary Ionic and Nonionic Surfactant Mixtures.

Interfacial vibrational footprints of the binary mixture of sodium dodecyl sulfate (SDS) and hexaethylene glycol monododecyl ether (C12E6) were probed using heterodyne detected vibrational sum frequency generation (HDVSFG). Our results show that in the presence of C12E6 at CMC (70 µM) the effect of SDS on the orientation of interfacial water molecules is enhanced 10 times compared to just pure surfactants. The experimental results contest the traditional Langmuir adsorption model predictions. This is also evidenced by our molecular dynamics simulations that show a remarkable restructuring and enhanced orientation of the interfacial water molecules upon DS- adsorption to the C12E6 surface. The simulations show that the adsorption free energy of DS- ions to a water surface covered with C12E6 is an enthalpy-driven process and more attractive by ∼10 kBT compared to the adsorption energy of DS- to the surface of pure water.

Non-Aromatic Fluorescence in Biological Matter: The Exception or the Rule?

While in the vast majority of cases fluorescence in biological matter has been attributed to aromatic or conjugated groups, peptides associated with neurodegenerative diseases, such as Alzheimer's, Parkinson's, or Huntington's, have been recently shown to display an intrinsic visible fluorescence even in the absence of aromatic residues. This has called the attention of researchers from many different fields, trying to understand the origin of this peculiar behavior and, at the same time, motivating the search for novel strategies to control the optical properties of new biophotonic materials. Today, after nearly 15 years of its discovery, there is a growing consensus about the mechanism underlying this phenomenon, namely, that electronic interactions between non-optically active molecules can result in supramolecular assemblies that are fluorescent. Despite this progress, many aspects of this phenomenon remain uncharted territory. In this Perspective, we lay down the state-of-the-art in the field highlighting the open questions from both experimental and theoretical fronts in this fascinating emerging area of non-aromatic fluorescence.

Charge Gradients around Dendritic Voids Cause Nanoscale Inhomogeneities in Liquid Water.

Water is the matrix of life and serves as a solvent for numerous physical and chemical processes. The origins of the nature of inhomogeneities that exist in liquid water and the time scales over which they occur remains an open question. Here, we report femtosecond elastic second harmonic scattering (fs-ESHS) of liquid water in comparison to an isotropic liquid (CCl4) and show that water is indeed a nonuniform liquid. The coherent fs-ESHS intensity was interpreted, using molecular dynamics simulations, as arising from charge density fluctuations with enhanced nanoscale polarizabilities around transient voids having an average lifetime of 300 fs. Although voids were also present in CCl4, they were not characterized by hydrogen bond defects and did not show strong polarizability fluctuations, leading to fs-ESHS of an isotropic liquid. The voids increased in number at higher temperatures above room temperature, in agreement with the fs-ESHS results.