Dielectric constant and low-frequency infrared spectra for liquid water and ice Ih within the E3B model

Terahertz Spectroscopy of Aqueous Solutions of Sodium Halides (NaX): Self- and Cross-Correlation Contributions of Ions and Hydration Shell Water for X- = F-, Cl-, Br-, and I. J Phys Chem B 2023;127:9323-9335. [PMID: 37871257 DOI: 10.1021/acs.jpcb.3c05228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

We investigated the terahertz (THz) absorption spectra of aqueous sodium halide solutions through molecular dynamics simulations using polarizable models of both water and ions. Specifically, we have considered aqueous solutions (∼1 M) of NaF, NaCl, NaBr, and NaI and calculated the difference THz spectrum of these solutions by subtracting the corresponding pure water contribution from the total THz spectrum of an ionic solution. The difference absorption spectrum of a given solution is then dissected into contributions from ion and ion-water correlations and also modifications of water-water correlations in the presence of the ions. The different components are further decomposed into induced dipole and permanent charge/dipole components and also into self- and cross-correlation components. The ion-water cross-correlation components are subsequently decomposed into contributions coming from different solvation shells through radially resolved calculations of such ion-water cross-correlations. Through all of these dissections, we could investigate the origin of different parts of the difference THz spectra of the sodium halide solutions studied here. It is found that while features below or around 100 cm-1 and also around 200 cm-1 arise mainly from ion and ion-water motion, that at the librational region above 600 cm-1 primarily originates from changes in water librational motion influenced by the ions. The variations of intensities of different components are also linked to the size and charge density of the anions in the solutions.

Incorporating Polarization and Charge Transfer into a Point-Charge Model for Water Using Machine Learning

Rigid nonpolarizable water models with fixed point charges have been widely employed in molecular dynamics simulations due to their efficiency and reasonable accuracy for the potential energy surface. However, the dipole moment surface of water is not necessarily well-described by the same fixed charges, leading to failure in reproducing dipole-related properties. Here, we developed a machine-learning model trained against electronic structure data to assign point charges for water, and the resulting dipole moment surface significantly improved the predictions of the dielectric constant and the low-frequency IR spectrum of liquid water. Our analysis reveals that within our atom-centered point-charge description of the dipole moment surface, the intermolecular charge transfer is the major source of the peak intensity at 200 cm^-1, whereas the intramolecular polarization controls the enhancement of the dielectric constant. The effects of exact Hartree-Fock exchange in the hybrid density functional on these properties are also discussed.

Isochores and Heat Capacity of Liquid Water in Terms of the Ion–Molecular Model

Thermodynamics of liquid water in terms of a non-standard approach—the ion–molecular model—is considered. Water is represented as a dense gas of neutral H2O molecules and single charged H3O+ and OH− ions. The molecules and ions perform thermal collisional motion and interconvert due to ion exchange. The energy-rich process—vibrations of an ion in a hydration shell of molecular dipoles—well known to spectroscopists with its dielectric response at 180 cm−1 (5 THz), is suggested to be key for water dynamics. Taking into account this ion–molecular oscillator, we compose an equation of state of liquid water to obtain analytical expressions for the isochores and heat capacity.

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree-Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

Superheating and Homogeneous Melting Dynamics of Bulk Ice

Homogeneous melting of crystals is a complex multistep process involving the formation of transient states at temperatures considerably higher than the melting point. The nature and persistence of these metastable structures are intimately connected to the melting process, and a precise definition of the temporal boundaries of these phenomena is not yet available. We set up a specifically designed experiment to probe by transient infrared absorption spectroscopy the entire dynamics, ranging from tens of picoseconds to microseconds, of superheating and melting of an ice crystal. In spite of a large excess of energy provided, only about 30% of the micrometric crystal liquefies in the first 20-25 ns because of the long persistence of the superheated metastable phase that extends for more than 100 ns. This behavior is ascribed to the population of low-energy states that trap a large amount of energy, favoring the formation of a metastable, likely plastic, ice phase.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

Non-equilibrium thermodynamics and collective vibrational modes of liquid water in an inhomogeneous electric field

In this experiment liquid water is subject to an inhomogeneous electric field (∇(2)Ea≈ 10(10) V m(2)) using a high voltage (20 kV) point-plane electrode system. Using interferometry it was found that the application of a strong electric field gradient to water generates local changes in the refractive index of the liquid, polarizes the surface and creates a downward moving electro-convective jet. A maximum temperature difference of 1 °C is measured in the immediate vicinity of the point electrode. Raman spectroscopy performed on water reveals an enhancement of the vibrational collective modes (3250 cm(-1)) as well as an increase in the local mode (3490 cm(-1)) energy. This bimodal enhancement indicates that the spectral changes are not due to temperature changes. The intense field gradient thus establishes an excited subpopulation of vibrational oscillators far from thermal equilibrium. Delocalization of the collective vibrational mode spatially expands this excited population beyond the microscale. Hindered rotational freedom due to electric field pinning of molecular dipoles retards the heat flow and generates a chemical potential gradient. These changes are responsible for the observed changes in the refractive index and temperature. It is demonstrated that polar liquids can thus support local non-equilibrium thermodynamic transient states critical to biochemical and environmental processes.

Heavy snow: IR spectroscopy of isotope mixed crystalline water ice

Temperature and isotopic dependence of simulated and experimental spectra shed light on the vibrational modes of crystalline water ice.