First principles multielectron mixed quantum/classical simulations in the condensed phase. II. The charge-transfer-to-solvent states of sodium anions in liquid tetrahydrofuran

Gas-phase atomic anions lack bound electronic excited states, yet in solution many of these anions exhibit intense absorption bands due to the presence of excited states, referred to as charge-transfer-to-solvent (CTTS) states that are bound only by the presence of the solvent. CTTS spectra thus serve as delicate probes of solute-solvent interactions, but the fact that they are created by the interactions of a solute with many solvent molecules makes them a challenge to describe theoretically. In this paper, we use mixed quantum/classical molecular dynamics with the two-electron Fourier-grid (2EFG) electronic structure method presented in the previous paper [W. J. Glover, R. E. Larsen, and B. J. Schwartz, J. Chem. Phys. 132, 144101 (2010)] to simulate the CTTS states of a sodium anion in liquid tetrahydrofuran, Na(-)/THF. Since our 2EFG method is based on configuration interaction with single and double excitations in a grid basis, it allows for an exact treatment of the two valence electrons of the sodium anion. To simulate Na(-)/THF, we first develop a new electron-THF pseudopotential, and we verify the accuracy of this potential by reproducing the experimental absorption spectrum of an excess electron in liquid THF with near quantitative accuracy. We also are able to reproduce the CTTS spectrum of Na(-)/THF and find that the CTTS states of Na(-) exhibit a Rydberg-like progression due to the pre-existing long-range solvent polarization around the anion. We also find that the CTTS states are highly mixed with the disjoint electronic states supported by naturally occurring solvent cavities that exist in liquid THF. This mixing explains why the solvated electrons that are ejected following CTTS excitation appear with their equilibrium absorption spectrum. The mixing of the CTTS and solvent-cavity states also explains why the recombination of the electron and its geminate Na(0) partner occurs on slower time scales when photoexciting in the blue compared to in the red of the CTTS band: blue excitation accesses CTTS states whose charge densities lies further from the position of the anion, whereas red excitation accesses CTTS states that lie primarily within the anion's first solvation shell. Finally, we see that the radial character of the CTTS states near the Na(+) core matches that of Na(0), explaining why the spectrum of this species appears instantly after photoexciting Na(-).