Going beyond the frozen core approximation: Development of coordinate-dependent pseudopotentials and application to Na 2+

Solvent Control of Chemical Identity Can Change Photodissociation into Photoisomerization

In solution-phase chemistry, the solvent is often considered to be merely a medium that allows reacting solutes to encounter each other. In this work, however, we show that moderate locally specific solute-solvent interactions can affect not only the nature of the solute but also the types of reactive chemistry. We use quantum simulation methods to explore how solvent participation in solute chemical identity alters reactions involving the breaking of chemical bonds. In particular, we explore the photoexcitation dynamics of Na₂⁺ dissolved in liquid tetrahydrofuran. In the gas phase, excitation of Na₂⁺ directly leads to dissociation, but in solution, photoexcitation leads to an isomerization reaction involving rearrangement of the first-shell solvent molecules; this isomerization must go to completion before the solute can dissociate. Despite the complexity, the solution-phase reaction dynamics can be captured by a two-dimensional energy surface where one dimension involves only the isomerization of the first-shell solvent molecules.

Hydrated Electrons in High-Concentration Electrolytes Interact with Multiple Cations: A Simulation Study

Experimental studies have demonstrated that the hydrated electron's absorption spectrum undergoes a concentration-dependent blue-shift in the presence of electrolytes such as NaCl. The blue-shift increases roughly linearly at low salt concentration but saturates as the solubility limit of the salt is approached. Previous attempts to understand the origin of the concentration-dependent spectral shift using molecular simulation have only examined the interaction between the hydrated electron and a single sodium cation, and these simulations predicted a spectral blue-shift that was an order of magnitude larger than that seen experimentally. Thus, in this paper, we first explore the reasons for the exaggerated spectral blue-shift when a simulated hydrated electron interacts with a single Na⁺. We find that the issue arises from nonpairwise additivity of the Na⁺-e^- and H₂O-e^- pseudopotentials used in the simulation. This effect arises because the solvating water molecules donate charge into the empty orbitals of Na⁺, lowering the effective charge of the cation and thus reducing the excess electron-cation interaction. Careful analysis shows, however, that although this nonpairwise additivity changes the energetics of the electron-Na⁺ interaction, the forces between the electron, Na⁺, and water are unaffected, so that mixed quantum/classical (MQC) simulations produce the correct structure and dynamics. With this in hand, we then use MQC simulations to explore the behavior of the hydrated electron as an explicit function of NaCl salt concentration. We find that the simulations correctly reproduce the observed experimental spectral shifting behavior. The reason for the spectral shift is that as the electrolyte concentration increases, the average number of cations simultaneously interacting in contact pairs with the hydrated electron increases from 1.0 at low concentrations to ∼2.5 near the saturation limit. As the number of cations that interact with the electron increases, the cation/electron interactions becomes slightly weaker, so that the corresponding Na⁺-e^- distance increases with increasing salt concentration. We also find that the dielectric constant of the solution plays little role in the observed spectroscopy, so that the electrolyte-dependent spectral shifts of the hydrated electron are directly related to the concentration-dependent number of closely interacting cations.

To be or not to be in a cavity: the hydrated electron dilemma

The hydrated electron-the species that results from the addition of a single excess electron to liquid water-has been the focus of much interest both because of its role in radiation chemistry and other chemical reactions, and because it provides for a deceptively simple system that can serve as a means to confront the predictions of quantum molecular dynamics simulations with experiment. Despite all this interest, there is still considerable debate over the molecular structure of the hydrated electron: does it occupy a cavity, have a significant number of interior water molecules, or have a structure somewhere in between? The reason for all this debate is that different computer simulations have produced each of these different structures, yet the predicted properties for these different structures are still in reasonable agreement with experiment. In this Feature Article, we explore the reasons underlying why different structures are produced when different pseudopotentials are used in quantum simulations of the hydrated electron. We also show that essentially all the different models for the hydrated electron, including those from fully ab initio calculations, have relatively little direct overlap of the electron's wave function with the nearby water molecules. Thus, a non-cavity hydrated electron is better thought of as an "inverse plum pudding" model, with interior waters that locally expel the surrounding electron's charge density. Finally, we also explore the agreement between different hydrated electron models and certain key experiments, such as resonance Raman spectroscopy and the temperature dependence and degree of homogeneous broadening of the optical absorption spectrum, in order to distinguish between the different simulated structures. Taken together, we conclude that the hydrated electron likely has a significant number of interior water molecules.

Comment on "Going beyond the frozen core approximation: development of coordinate-dependent pseudopotentials and application to Na2(+)" [J. Chem. Phys. 138, 054110 (2013)]

The new coordinate-dependent pseudopotential for Na2(+) by Kahros and Schwartz [J. Chem. Phys. 138, 054110 (2013)] is assessed and compared to the pseudopotential approach by Fuentealba et al. [Chem. Phys. Lett. 89, 418 (1982)] which incorporates the coordinate-dependent core-polarization potential by Müller and Meyer [J. Chem. Phys. 80, 3311 (1984)]. In contrast to the latter approach, the one by Kahros and Schwartz does not reproduce the accurately known experimental data and∕or high level theoretical results for Na2(+). The treatment of core polarization by Kahros and Schwartz neglects the dynamic polarization of atomic cores which is much more important for Na2(+) than the static one. On the other hand, the Kahros and Schwartz method heavily overestimates frozen-core corrections at the Hartree-Fock level by compounding them with artifacts of a superposition of non-norm-conserving pseudopotentials.

Response to "Comment on 'Going beyond the frozen core approximation: development of coordinate-dependent pseudopotentials and application to Na2(+)'" [J. Chem. Phys. 139, 147101 (2013)]

Stoll, Fuentealba, and Szentpály (SFS) argue that the coordinate-dependent pseudopotential we developed for the sodium dimer cation molecule is inferior to other potentials that have been presented in the literature for this molecule. The goal of our work, however, was to present a novel method for the development of rigorous coordinate-dependent pseudopotentials. Our method is designed to reproduce all-electron Hartree-Fock calculations without the inclusion of adjustable parameters. Moreover, our method starts from the superposition of unoptimized, non-norm-conserved atomic potentials, so that when complete, the resulting norm-conserving potential can reproduce an all-electron Hartree-Fock calculation without the inclusion of adjustable parameters. We chose the sodium dimer cation system as a proof of principle for our method, and showed that our method does indeed allow a one-electron calculation to correctly reproduce the all-electron Hartree-Fock calculation from bonding to the dissociation limit. Our purpose in developing this method is to use such potentials in condensed-phase mixed quantum/classical molecular dynamics simulations, where inclusion of valence polarization effects is unimportant or can be added on after the fact. Thus we do not claim that our method provides a potential that is superior to potentials that have been specifically constructed to go beyond the static exchange approximation and/or include valence polarization effects-such potentials are beyond the scope of our work. We also note that although we made a numerical error in the application of our method to Na2(+) in our original work [A. Kahros and B. J. Schwartz, J. Chem. Phys. 138, 054110 (2013)] that led to an overestimation of the magnitude of core polarization effects for this particular molecule, out method does work as derived for this molecule and the error does not affect the significance of our method or its general applicability.