No Arabic abstract
Photoelectron spectroscopy experiments in ionic solutions reveal important electronic structure information, in which the interaction between hydrated ions and water solvent can be inferred. Based on many-body perturbation theory with GW approximation, we theoretically compute the quasiparticle electronic structure of chloride anion solution, which is modeled by path-integral $ab$ $initio$ molecular dynamics simulation by taking account the nuclear quantum effects (NQEs). The electronic levels of hydrated anion as well as water are determined and compared to the recent experimental photoelectron spectra. It is found that NQEs improve the agreement between theoretical prediction and experiment because NQEs effectively weaken the hybridization of the between the $rm Cl^-$ anion and water. Our study indicates that NQEs plays a small but non-negligible role in predicting the electronic structure of the aqueous solvation of ions of the Hofmeister series.
Path-integral ab initio molecular dynamics (PI-AIMD) calculations have been employed to probe the nature of chloride ion solvation in aqueous solution. Nuclear quantum effects (NQEs) are shown to weaken hydrogen bonding between the chloride anion and the solvation shell of water molecules. As a consequence, the disruptive effect of the anion on the solvent water structure is significantly reduced compared to what is found in the absence of NQEs. The chloride hydration structure obtained from PI-AIMD agrees well with information extracted from neutron scattering data. Inparticular, the observed satellite peak in the hydrogen-chloride-hydrogen triple angular distribution serves as a clear signature of NQEs. The present results suggest that NQEs are likely to play acrucial role in determining the structure of saline solutions.
We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQE) in molecular simulations: the adaptive Quantum Thermal Bath (adQTB). In this method, zero point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation-dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and we report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.
We study the solvation and electrostatic properties of bare gold (Au) nanoparticles (NPs) of $1$-$2$ nm in size in aqueous electrolyte solutions of sodium salts of various anions with large physicochemical diversity (Cl$^-$, BF$_4$$^-$, PF$_6$$^-$, Nip$^-$(nitrophenolate), 3- and 4-valent hexacyanoferrate (HCF)) using nonpolarizable, classical molecular dynamics computer simulations. We find a substantial facet selectivity in the adsorption structure and spatial distribution of the ions at the Au-NPs: while sodium and some of the anions (e.g., Cl$^-$, HCF$^{3-}$) adsorb more at the `edgy (100) and (110) facets of the NPs, where the water hydration structure is more disordered, other ions (e.g., BF$_4$$^-$, PF$_6$$^-$, Nip$^-$) prefer to adsorb strongly on the extended and rather flat (111) facets. In particular, Nip$^-$, which features an aromatic ring in its chemical structure, adsorbs strongly and perturbs the first water monolayer structure on the NP (111) facets substantially. Moreover, we calculate adsorptions, radially-resolved electrostatic potentials, as well as the far-field effective electrostatic surface charges and potentials by mapping the long-range decay of the calculated electrostatic potential distribution onto the standard Debye-Huckel form. We show how the extrapolation of these values to other ionic strengths can be performed by an analytical Adsorption-Grahame relation between effective surface charge and potential. We find for all salts negative effective surface potentials in the range from $-10$ mV for NaCl down to about $-80$ mV for NaNip, consistent with typical experimental ranges for the zeta-potential. We discuss how these values depend on the surface definition and compare them to the explicitly calculated electrostatic potentials near the NP surface, which are highly oscillatory in the $pm 0.5$ V range.
This paper presents a combined theoretical and experimental investigation of aqueous near-neutral electrolytes based on chloride salts for rechargeable zinc-air batteries (ZABs). The resilience of near-neutral chloride electrolytes in air could extend ZAB lifetime, but theory-based simulations predict that such electrolytes are vulnerable to other challenges including pH instability and the unwanted precipitation of mixed zinc hydroxide chloride products. In this work, we combine theory-based simulations with experimental methods such as full cell cycling, operando pH measurements, ex-situ XRD, SEM, and EDS characterization to investigate the performance of ZABs with aqueous chloride electrolytes. The experimental characterization of near-neutral ZAB cells observes the predicted pH instability and confirms the composition of the final discharge products. Steps to promote greater pH stability and control the precipitation of discharge products are proposed.
We discuss the analytic and diagrammatic structure of ionization potential (IP) and electron affinity (EA) equation-of-motion coupled-cluster (EOM-CC) theory, in order to put it on equal footing with the prevalent $GW$ approximation. The comparison is most straightforward for the time-ordered one-particle Greens function, and we show that the Greens function calculated by EOM-CC with single and double excitations (EOM-CCSD) includes fewer ring diagrams at higher order than does the $GW$ approximation, due to the formers unbalanced treatment of time-ordering. However, the EOM-CCSD Greens function contains a large number of vertex corrections, including ladder diagrams, mixed ring-ladder diagrams, and exchange diagrams. By including triple excitations, the EOM-CCSDT Greens function includes all diagrams contained in the $GW$ approximation, along with many high-order vertex corrections. In the same language, we discuss a number of common approximations to the EOM-CCSD equations, many of which can be classified as elimination of diagrams. Finally, we present numerical results by calculating the principal charged excitations energies of the molecules contained in the so-called $GW$100 test set [J. Chem. Theory Comput. 2015, 11, 5665-5687]. We argue that (in molecules) exchange is as important as screening, advocating for a Hartree-Fock reference and second-order exchange in the self-energy.