When a solute molecule is placed in solution, is it acceptable to presume that its electronic structure is essentially the same as that in the gas phase? In this paper, we address this question from a simulation perspective for the case of the sodium dimer cation (Na2+) molecule in both liquid Ar and liquid tetrahydrofuran (THF). In previous work, we showed that, when local specific interactions between a solute and solvent are energetically on the order of a hydrogen bond, the solvent can become part of the chemical identity of the solute. Here, using mixed quantum/classical molecular dynamics simulations, we see that, for the Na2+ molecule, solute-solvent interactions lead to two stable, chemically distinct coordination states (Na(THF)4-Na(THF)5+ and Na(THF)5-Na(THF)5+) that are not only stable themselves as gas-phase molecules but that also have a completely new electronic structure with important implications for the excited-state photodissociation of this molecule in the condensed phase. Furthermore, we show through a set of comparative classical simulations that treating the solute's bonding electron explicitly quantum mechanically is necessary to understand both the ground-state dynamics and chemical identity of this simple diatomic molecule; even use of the quantum-derived potential of mean force is insufficient to describe the behavior of the molecule classically. Finally, we calculate the results of a proposed transient hole-burning experiment that could be used to spectroscopically disentangle the presence of the different coordination states.
Read full abstract