The Galli group webpage

Water and simple aqueous solutions

First principles simulations of liquid water are an ongoing challenge.
Accounting for the properties of water requires the ability to describe different types of bonds: intra-molecular covalent bonds, inter-molecular hydrogen bonds and van der Waals interactions. None of the known density functionals or hybrid functionals is capable of giving an equally accurate description of all these bonds. Additional difficult tasks include accounting for proton quantum effects, and efficiently sampling the complex potential energy surface of water as a function of time.
In spite of these challenges, we are making progress in understanding solvation properties of simple ions, unraveling interactions between water and surfaces, and in predicting spectroscopic properties of aqueous systems and of water and ice under pressure. We mostly use ab initio molecular dynamics (MD) and the Qbox code, but also classical simulations to investigate ice nucleation processes and some of the properties of water yet inaccessible to ab initio MD.
Determining how the structure of water is modified by the presence of salts is instrumental to understanding the solvation of biomolecules and in general, the role played by salts in biochemical processes. Yet the extent of hydrogen bonding disruption induced by salts remains controversial. We performed extensive first-principles simulations of solutions of a simple salt (NaCl) and found that while the cation does not significantly change the structure of water beyond the first solvation shell, the anion has a further reaching effect, modifying the hydrogen-bond network even outside its second solvation shell. We found that a distinctive fingerprint of hydrogen bonding modification is the change in polarizability of water molecules. Molecular dipole moments are instead insensitive probes of long-range modifications induced by Na+ and Cl- ions. Though noticeable, the long-range effect of Cl- is expected to be too weak to affect solubility of large biomolecules. We obtained results using both PBE and hybrid functionals.
  • "Local and Global Effects of Dissolved Sodium Chloride on the Structure of Water ", Alex P. Gaiduk and Giulia Galli, J. Phys. Chem. Lett. 8, pp 1496–1502 (2017)
  • "Structural and electronic properties of aqueous NaCl solutions from ab initio molecular dynamics simulations with hybrid density functionals", Alex P. Gaiduk, Cui Zhang, François Gygi, and Giulia Galli, Chem. Phys. Lett. 604, 89 (2014)

We investigated the photoelectron spectrum of a simple aqueous solution of NaCl, by combining ab initio calculations and experiments. Measurements were conducted on microjets, and first-principles calculations were performed using hybrid functionals and many-body perturbation theory at the G0W0 level, starting with wave functions computed in ab initio molecular dynamics simulations. We found excellent agreement between theory and experiments for the positions of both the solute and solvent excitation energies on an absolute energy scale and for peak intensities. The best comparison was obtained using wave functions obtained with dielectric-dependent self-consistent and range-separated hybrid functionals. Our computational protocol opens the way to accurate, predictive calculations of the electronic properties of electrolytes, of interest to a variety of energy problems.

Our ab initio simulations showed that the study of the electronic properties of anions in water requires the use of hybrid functionals, and of many body perturbation theory, in order to properly describe the correct alignment of their energy levels with those of water.
  • "Electronic Energy Levels and Band Alignment for Aqueous Phenol and Phenolate from First Principles", Daniel Opalka, Tuan Anh Pham, Giulia Galli and Michiel Sprik, J. Phys. Chem. B 119, 9651 (2015)
  • "Electronic Structure of Aqueous Sulfuric Acid from First-Principles Simulations with Hybrid Functionals", Quan Wan, Leonardo Spanu, Francois Gygi, and Giulia Galli, J. Phys. Chem. Lett. 5 (15), 2562 (2014)
  • "The ionization potential of aqueous hydroxide computed using many-body perturbation theory", Daniel Opalka, Tuan Anh Pham, Michiel Sprik, and Giulia Galli, J. Chem. Phys. 5, 2562 (2014)
  • "Structural and electronic properties of aqueous NaCl solutions from ab initio molecular dynamics simulations with hybrid density functionals", Alex P. Gaiduk, Cui Zhang, Francois Gygi, and Giulia Galli, Chem. Phys. Lett. 604, 89 (2014)
  • "Electronic structure of the solvated chloride anion from first principles molecular dynamics", C.Zhang, T.A.Pham, F.Gygi and G.Galli, J. Chem. Phys- Commun. 138, 181102 (2013)
  • "Solvation Properties of Microhydrated Sulfate Anion Clusters: Insights From Ab Initio Calculations", Quan Wan, Leonardo Spanu and Giulia Galli, J. Phys. Chem. B 116, 9460 (2012)

We determined the equilibrium density and compressibility of water and ice from first-principles molecular dynamics simulations using gradient-corrected (PBE) and hybrid (PBE0) functionals and the Qbox code. Both functionals predicted the density of ice to be larger than that of water, by 15 (PBE) and 35% (PBE0). The PBE0 functional yielded a lower density of both ice and water with respect to PBE, leading to better agreement with experiment for ice but not for liquid water. Approximate inclusion of dispersion interactions on computed molecular-dynamics trajectories led to a substantial improvement of the PBE0 results for the density of liquid water, which, however, resulted to be slightly lower than that of ice.
  • "Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals", Alex P. Gaiduk, François Gygi, and Giulia Galli, J. Phys. Chem. Lett. 6, 2902 (2015)

The use of hybrid functionals - in particular PBE0 - improves the description of the structural, diffusion and vibrational properties of liquid water, with respect to that provided by semi-local density functionals (e.g. PBE). Such an improved description stems from two effects: a more accurate account, at the PBE0 level of theory, of the vibrational properties of the monomer and dimer, and a structural model for the liquid with a smaller number of hydrogen bonds and oxygen coordination than those obtained with semilocal functionals. Van der Waals (vdW) functionals may also improve the description of water and ice properties, with respect to PBE, depending on the choice of the local exchange functional.

Ab initio infrared spectra of water confined between non polar surfaces show that electronic charge fluctuations at the interface occur even in the case of highly hydrophobic substrates. These results indicate that hydrophobic surfaces may not be regarded as hard, inert walls. Electronic fluctuations are responsible for specific features present in IR signals and for important differences between IR spectra and vibrational densities of states.
  • "Electronic effects in the IR spectrum of water under confinement", D. Donadio, G. Cicero, E. Schwegler, M. Sharma and G.Galli, J. Phys. Chem. B., 113, 4170 (2009)
  • "Probing properties of water under confinement: Infrared spectra", M.Sharma, D.Donadio, E.Schwegler and G.Galli, Nano Lett., 8, 2959 (2008)

The model of water confined between non polar substrates obtained from ab initio simulations points at a highly inhomogeneous, mobile interfacial layer. This layer is composed of a region of about 2.5 Å of zero particle density (and depleted electronic density), followed by a region with particle density higher than in the bulk, with molecules exhibiting a smaller dipole moment than in bulk water. In confined salty water simulated from first principles, simple ions such as Na+ and Cl- preferentially reside at the interface.
  • "Probing the Structure of Salt Water under Confinement with First-Principles Molecular Dynamics and Theoretical X-ray Absorption Spectroscopy", Heather J. Kulik, Eric Schwegler and Giulia Galli, J. Phys. Chem. Lett., 3, 2653-2658 (2012)
  • "Water confined in nanotubes and between graphene sheets: a first principle study", G. Cicero, J. C. Grossman, E. Schwegler, F. Gygi and G. Galli, J. Phys. Chem. B., 113, 4170 (2009)
  • "Water confined in carbon nanotubes: Magnetic response and proton chemical shieldings", P. Huang, E. Schwegler and G. Galli, J. Phys. Chem. C, 113, 8696 (2009)

First principles molecular dynamics simulations and ab initio calculations of X-ray absorption spectra were used to investigate the aqueous solvation of cations in MgCl2, CaCl2, and NaCl solutions and of salty water confined in nanotubes. The results helped rationalize and interpret several experimental findings. In particular spectral signatures of single donor, double donor and acceptor HBs were identified, and assigned to specific configurations of water molecules in the solvation shells of the mono- and di-valent cations.
  • "Probing the Structure of Salt Water under Confinement with First-Principles Molecular Dynamics and Theoretical X-ray Absorption Spectroscopy", Heather J. Kulik, Eric Schwegler and Giulia Galli, J. Phys. Chem. Lett., 3, 2653-2658 (2012)
  • "Local Effects in the X-ray Absorption Spectrum of Salt Water", H. J. Kulik, N. Marzari, A. A. Correa, D. Prendergast, E. Schwegler and G. Galli, J. Phys. Chem. B, 114, 9594 (2010)

The calculation of the static dielectric screening in water and ice allowed us to identify intra- and inter-molecular screening modes (as linear combinations of the eigenvectors of the dielectric matrix). The results indicate that in ice and water screening has strong intermolecular components and cannot be simply derived from perturbed values of gas phase molecular polarizabilities. The work on dielectric screening led to find efficient ways to represent large dielectric matrices, that are now used in many body perturbation theory calculations (BSE and GW) of absorption and emission spectra of molecules and nanostructures.


Top