Molecular Polarizabilities in Aqueous Systems from First-Principles

The heterogeneity of aqueous solutions: the current situation in the context of experiment and theory

The advancement of experimental methods has provided new information about the structure and structural fluctuations of water. Despite the appearance of numerous models, which aim to describe a wide range of thermodynamic and electrical characteristics of water, there is a deficit in systemic understanding of structuring in aqueous solutions. A particular challenge is the fact that even pure water is a heterogeneous, multicomponent system composed of molecular and supramolecular structures. The possibility of the existence of such structures and their nature are of fundamental importance for various fields of science. However, great difficulties arise in modeling relatively large supramolecular structures (e.g. extended hydration shells), where the bonds between molecules are characterized by low energy. Generally, such structures may be non-equilibrium but relatively long-lived. Evidently, the short times of water microstructure exchanges do not mean short lifetimes of macrostructures, just as the instability of individual parts does not mean the instability of the entire structure. To explain this paradox, we review the data from experimental and theoretical research. Today, only some of the experimental results on the lifetime of water structures have been confirmed by modeling, so there is not a complete theoretical picture of the structure of water yet. We propose a new hierarchical water macrostructure model to resolve the issue of the stability of water structures. In this model, the structure of water is presented as consisting of many hierarchically related levels (the stratification model). The stratification mechanism is associated with symmetry breaking at the formation of the next level, even with minimal changes in the properties of the previous level. Such a hierarchical relationship can determine the unique physico-chemical properties of water systems and, in the future, provide a complete description of them.

Nonlinear Effects of Hydrophobic Confinement on the Electronic Structure and Dielectric Response of Water

Fundamental studies of the dielectrics of confined water are critical to understand the ion transport across biological and synthetic nanochannels. The relevance of these fundamental studies, however, surmounts the difficulty of probing water's dielectric constant as a function of a fine variation in confinement. In this work, we explore the computational efficiency of machine learning potentials to derive the confinement effects on the dielectric constant, polarization, and dipole moment of water. Our simulations predict an enhancement of the axial dielectric constant of water under extreme confinement, arising from either the formation of ferroelectric structures of ordered water or larger dipole fluctuations facilitated by the disruption of water's H-bond network. Our study highlights the impact of hydrophobic nanoconfinement on the dielectric constant and on the ionic and electronic structure of water molecules, pointing to the importance of geometric flexibility and electronic polarizability to properly model confinement effects on water.

Raman Spectra of Electrified Si-Water Interfaces: First-Principles Simulations

We investigate the Raman spectra of liquid water in contact with a semiconductor surface using first-principles molecular dynamics simulations. We focus on a hydrogenated silicon-water interface and compute the Raman spectra from time correlation functions of the polarizability. We establish a relationship between Raman spectral signatures and structural properties of the liquid at the interface, and we identify the vibrational impacts of an applied electric field. We show that negative bias leads to a reduction of the number of hydrogen bonds (HBs) formed between the surface and the topmost water layer and an enhancement of the HB interactions between water molecules. Instead, positive bias leads to an enhancement of both the HB interactions between water and the surface and between water molecules, creating a semi-ordered interfacial layer. Our work provides molecular-level insights into electrified semiconductor/water interfaces and the identification of specific structural features through Raman spectroscopy.

Accurate and Efficient Prediction of Post-Hartree-Fock Polarizabilities of Condensed-Phase Systems

To accurately and efficiently predict the molecular response properties (such as polarizability) at post-Hartree-Fock levels for condensed-phase systems under periodic boundary conditions (PBC) is still an unaccomplished and ongoing task. We demonstrate that static isotropic polarizabilities can be cost-effectively predicted at post-Hartree-Fock levels by combining the linear-scaling generalized energy-based fragmentation (GEBF) and information-theoretic approach (ITA) quantities. In PBC-GEBF, the total molecular polarizability of an extended system is obtained as a linear combination of the corresponding quantities of a series of small embedded subsystems of several monomers. Here, we show that in the PBC-GEBF-ITA framework, one can obtain the molecular polarizabilities and establish linear relations to ITA quantities. Once these relations are established for smaller subsystems, one can predict the polarizabilities of larger subsystems directly from the molecular wavefunction (or electron density) via ITA quantities. Alternatively, one can determine the total molecular polarizability via a linear combination equation in PBC-GEBF. We have corroborated that this newly proposed PBC-GEBF-ITA protocol is much more efficient than the original PBC-GEBF approach but is not much less accurate and that this conclusion holds for both many-body perturbation theory and the coupled cluster calculations. Good efficiency and transferability of the PBC-GEBF-ITA protocol are demonstrated for periodic systems with several hundred atoms in a unit cell.

Knowledge-Based Atomic Polarizabilities Used to Model Circular Dichroism Spectra of Proteins

We present a model of circular dichroism for proteins, which is mainly based on both the classical theory of optical activity and a series of effective atomic polarizabilities. Such polarizabilities are extracted from the analysis of a set of synchrotron radiation circular dichroism spectra and their corresponding three-dimensional structures from the Protein Data Bank. Each modeled spectrum is obtained from the protein atomic coordinates and the identification of its secondary structure elements. The resulting spectra are in good agreement with additional experimental data and also with the predictions of some other models. Among them, only our approach is able to describe the effect of d-amino acids. Moreover, our model is also utilized to evaluate protein reconstructions as well as structural changes.

Highly Heterogeneous Polarization and Solvation of Gold Nanoparticles in Aqueous Electrolytes

The performance of gold nanoparticles (NPs) in applications depends critically on the structure of the NP-solvent interface, at which the electrostatic surface polarization is one of the key characteristics that affects hydration, ionic adsorption, and electrochemical reactions. Here, we demonstrate significant effects of explicit metal polarizability on the solvation and electrostatic properties of bare gold NPs in aqueous electrolyte solutions of sodium salts of various anions (Cl^-, BF₄^-, PF₆^-, nitrophenolate, and 3- and 4-valent hexacyanoferrate), using classical molecular dynamics simulations with a polarizable core-shell model for the gold atoms. We find considerable spatial heterogeneity of the polarization and electrostatic potentials on the NP surface, mediated by a highly facet-dependent structuring of the interfacial water molecules. Moreover, ion-specific, facet-dependent ion adsorption leads to considerable alterations of the interfacial polarization. Compared to nonpolarizable NPs, surface polarization modifies water local dipole densities only slightly but has substantial effects on the electrostatic surface potentials and leads to significant lateral redistributions of ions on the NP surface. Besides, interfacial polarization effects cancel out in the far field for monovalent ions but not for polyvalent ions, as anticipated from continuum "image-charge" concepts. Far-field effective Debye-Hückel surface potentials change accordingly in a valence-specific fashion. Hence, the explicit charge response of metal NPs is crucial for the accurate description and interpretation of interfacial electrostatics (e.g., for charge transfer and interfacial polarization in catalysis and electrochemistry).