Ultrafast Dynamics of Chlorine−Water and Bromine−Water Radical Complexes Following Electron Photodetachment in Their Anionic Precursors

Prediction of Aqueous Reduction Potentials of X•, ChH•, and XO• Radicals with X = Halogen and Ch = Chalcogen

The aqueous electron affinity and aqueous reduction potentials for F•, Cl•, Br•, I•, OH•, SH•, SeH•, TeH•, ClO•, BrO•, and IO• were calculated using electronic structure methods for explicit cluster models coupled with a self-consistent reaction field (SMD) to treat the aqueous solvent. Calculations were conducted using MP2 and correlated molecular orbital theory up to the CCSD(T)-F12b level for water tetramer clusters and MP2 for octamer cluster. Inclusion of explicit waters was found to be important for accurately predicting the redox potentials in a number of cases. The calculated reduction potentials for X• and ChH• were predicted to within ∼0.1 V of the reported literature values. Fluorine is anomalous due to abstraction of a hydrogen from one of the surrounding water molecules to form a hydroxyl radical and hydrogen fluoride, so its redox potential was calculated using only an implicit model. Larger deviations from experiment were predicted for ClO• and BrO•. These deviations are due to the free energy of solvation of the anion being too negative, as found in the pKa calculations, and that for the neutral being too positive with the current approach.

Solvent Dynamics of Aqueous Halides before and after Photoionization

Electron transfer reactions can be strongly influenced by solvent dynamics. We study the photoionization of halides in water as a model system for such reactions. There are no internal nuclear degrees of freedom in the solute, allowing the dynamics of the solvent to be uniquely identified. We simulate the equilibrium solvent dynamics for Cl^-, Br^-, I^-, and their respective neutral atoms in water, comparing quantum mechanical/molecular mechanical (QM/MM) and classical molecular dynamics (MD) methods. On the basis of the obtained configurations, we calculate the extended X-ray absorption fine structure (EXAFS) spectra rigorously based on the MD snapshots and compare them in detail with other theoretical and experimental results available in the literature. We find our EXAFS spectra based on QM/MM MD simulations in good agreement with their experimental counterparts for the ions. Classical MD simulations for the ions lead to EXAFS spectra that agree equally well with the experiment when it comes to the oscillatory period of the signal, even though they differ from the QM/MM radial distribution functions extracted from the MD. The amplitude is, however, considerably overestimated. This suggests that to judge the reliability of theoretical simulation methods or to elucidate fine details of the atomistic dynamics of the solvent based on EXAFS spectra, the amplitude as well as the oscillatory period need to be considered. If simulations fail qualitatively, as does the classical MD for the aqueous neutral halogen atoms, the resulting EXAFS will also be strongly affected in both oscillatory period and amplitude. The good reliability of QM/MM-based EXAFS simulations, together with clear qualitative differences in the EXAFS spectra found between halides and their atomic counterparts, suggests that a combined theory and experimental EXAFS approach is suitable for elucidating the nonequilibrium solvent dynamics in the photoionization of halides and possibly also for electron transfer reactions in more complex systems.

Aqueous Microsolvation of Mercury Halide Species†

The effects of aqueous solvation on the thermochemistry of reactions between mercury and small halogen molecules has been investigated by the microsolvation approach using ab initio and density functional theory (DFT) calculations. The structures, vibrational frequencies, and binding energies of 1, 2, and 3 water molecules with mercury-halide (HgBr2, HgBrCl, HgCl2, HgBr, and HgCl) and related mercury and halogen species (Br2, BrCl, Cl2, Cl, Hg, and Br) have been computed with second order Møller-Plesset perturbation theory (MP2) and the B3LYP density functional method. Accurate incremental water binding energies have been obtained at the complete basis set (CBS) limit using sequences of correlation consistent basis sets, including higher order correlation effects estimated from coupled cluster calculations. The resulting energetics were used to calculate the influence of water molecules on the thermochemistry of a number of reactions between mercury and small halogen-containing molecules. In general, the presence of water favors the formation of oxidized mercury halide species.

Femtosecond dynamics of Cu(CD3OD)

We report the femtosecond nuclear dynamics of Cu(CD3OD) van der Waals clusters, investigated using photodetachment-photoionization spectroscopy. Photodetachment of an electron from Cu-(CD3OD) with a 150 fs, 398 nm laser pulse produces a vibrationally excited neutral complex that undergoes ligand reorientation and dissociation. The dynamics of Cu(CD3OD) on the neutral surface is interrogated by delayed femtosecond resonant two-photon ionization. Analysis of the resulting time-dependent signals indicates that the nascent Cu(CD3OD) complex dissociates on two distinct time scales of 3 and 30 ps. To understand the origins of the observed time scales, complimentary studies were performed. These included measurement of the photoelectron spectrum of Cu-(CD3OD) as well as a series of calculations of the structure and the electronic and vibrational energies of the anion and neutral complexes. Based on the comparisons of the experimental and calculated results for Cu(CD3OD) with those obtained from earlier studies of Cu(H2O), we conclude that the 3 ps time scale reflects the energy transfer from the rotation of CD3OD in the complex to the dissociation coordinate, while the 30 ps time scale reflects the energy transfer from the excited methyl torsion states to the dissociation coordinate.

Microsolvation of F- in Water

Microsolvation of F- in water is studied by ionization and double ionization spectra of (H2O)1-3F- calculated by ab initio methods. It is shown that the presence of the fluorine electrons introduces many-body properties in the spectra which cannot be reproduced by the presence of a negative point charge. The increase of the solvation shell increases the complexity in particular of the double ionization spectra. Ionization and double ionization energies slowly increase with continued solvation, and many-body effects in the inner valence spectra become more prominent.

Femtosecond dynamics of Cu(H2O)2

The ultrafast relaxation dynamics of Cu(H(2)O)(2) is investigated using femtosecond photodetachment-photoionization spectroscopy. In addition, stationary points on the Cu(H(2)O)(2) anion, neutral, and cation potential energy surfaces are characterized by ab initio electronic structure calculations. Electron photodetachment from Cu(-)(H(2)O)(2) initiates the dynamics on the ground-state potential energy surface of neutral Cu(H(2)O)(2). The resulting Cu(H(2)O)(2) complexes experience large-amplitude H(2)O reorientation and dissociation. The time evolution of the Cu(H(2)O)(2) fragmentation products is monitored by time-resolved resonant multiphoton ionization. The parent ion, Cu(+)(H(2)O)(2), is not detected above background levels. The rise to a maximum of the Cu(+) signal from Cu(-)(H(2)O)(2), and the decay of the Cu(+)(H(2)O) signal from Cu(-)(H(2)O)(2) have similar tau approximately 10 ps time dependences to the corresponding signals from Cu(-)(H(2)O), but display clear differences at very short and long times. The experimental observations can be understood in terms of the following picture. Prompt dissociation of H(2)O from nascent Cu(H(2)O)(2) gives rise to a vibrationally excited Cu(H(2)O) complex, which dissociates to Cu+H(2)O due to coupling of H(2)O internal rotation to the dissociation coordinate. This prompt dissociation removes all intra-H(2)O vibrational excitation from the intermediate Cu(H(2)O) fragment, which quenches the long time vibrational predissociation to Cu+H(2)O previously observed in analogous experiments on Cu(-)(H(2)O).

Femtosecond study of Cu(H2O) dynamics

The short-time nuclear dynamics of Cu(H(2)O) is investigated using femtosecond photodetachment-photoionization spectroscopy and time-dependent quantum wave packet calculations. The Cu(H(2)O) dynamics is initiated in the electronic ground state of the complex by electron photodetachment from the Cu(-)(H(2)O) complex, where hydrogen atoms are oriented toward Cu. Several time-resolved resonant multiphoton ionization schemes are used to probe the ensuing reorientation and dissociation. Immediately following photodetachment, the neutral complex is far from its minimum energy geometry and possesses an internal energy comparable to the Cu-H(2)O dissociation energy and undergoes both large-amplitude H(2)O motion and dissociation. Dissociation is observed to occur on three distinct time scales: 0.6, 8, and 100 ps. These results are compared to the results of time-dependent J=0 wave packet calculations, propagating the initial anion vibrational wave functions on the ground-state potential of the neutral complex. An excellent agreement is obtained between the experimental results and the ionization signals derived from the calculated probability amplitudes. Related experiments and calculations are carried out on the Cu(D(2)O) complex, with results very similar to those of Cu(H(2)O).