Electron solvation in water clusters following charge transfer from iodide

The dynamics following charge transfer to solvent from iodide to a water cluster are studied using time-resolved photoelectron imaging of I-(H2O)n and I-(D2O)n clusters with n< or =28. The results show spontaneous conversion, on a time scale of approximately 1 ps, from water cluster anions with surface-bound electrons to structures in which the excess electron is more strongly bound and possibly more internalized within the solvent network. The resulting dynamics provide valuable insight into the electron solvation dynamics in water clusters and the relative stabilities between recently observed isomers of water cluster anions.     

Ionic and electronic structure of metal clusters

Structure and electron dynamics of sodium clusters are investigated within the local-density approximation for the electrons. We compare results from detailed ionic structure with those from a structure averaged jellium model and find that the dominance of the electron cloud overlays most of the differences in the background. Ionic structure is indispensable, however, to compute the surface energy of clusters and to provide an unprejudiced picture of cluster fission. For all cases, we compute the resonance spectra associated with electron dynamics. In particular, the very strong deformations during fission deliver unusual resonance modes with a broad spectral fragmentation.     

Excess electron relaxation dynamics at water/air interfaces

We have performed mixed quantum-classical molecular dynamics simulations of the relaxation of a ground state excess electron at interfaces of different phases of water with air. The investigated systems included ambient water/air, supercooled water/air, Ih ice/air, and amorphous solid water/air interfaces. The present work explores the possible connections of the examined interfacial systems to finite size cluster anions and the three-dimensional infinite, fully hydrated electron. Localization site analyses indicate that in the absence of nuclear relaxation the electron localizes in a shallow potential trap on the interface in all examined systems in a diffuse, surface-bound (SB) state. With relaxation, the weakly bound electron undergoes an ultrafast localization and stabilization on the surface with the concomitant collapse of its radius. In the case of the ambient liquid interface the electron slowly (on the 10 ps time scale) diffuses into the bulk to form an interior-bound state. In each other case, the excess electron persists on the interface in SB states. The relaxation dynamics occur through distinct SB structures which are easily distinguishable by their energetics, geometries, and interactions with the surrounding water bath. The systems exhibiting the most stable SB excess electron states (supercooled water/air and Ih ice/air interfaces) are identified by their characteristic hydrogen-bonding motifs which are found to contain double acceptor-type water molecules in the close vicinity of the electron. These surface states correlate reasonably with those extrapolated to infinite size from simulated water cluster anions.     

Ab initio molecular dynamics simulation of a medium-sized water cluster anion: from an interior to a surface-located excess electron via a delocalized state

We present a computational study of the structure and dynamics of an excess electron in a medium-sized water cluster aimed at addressing the question of interior vs exterior solvation. Ab initio Born-Oppenheimer molecular dynamics simulations were performed within the DFT framework, employing a hybrid Gaussian and plane-wave formalism together with the PBE exchange-correlation functional and norm-conserving pseudopotentials. Analysis of a 15-ps trajectory allowed us to reach the following conclusions: (i) the excess electron is predominantly located at the cluster surface (even if it is initially placed in the interior), (ii) the computed electron binding energies correlate with the electron localization rather than with its bulk vs surface location, and (iii) a dynamical interconversion between two different H-bond patterns around the electron occurs. The computed electron binding energies and the most relevant features of the IR spectrum are in a very good agreement with results of previous experimental studies.     

Excess electron localization sites in neutral water clusters

We present approximate pseudopotential quantum-mechanical calculations of the excess electron states of equilibrated neutral water clusters sampled by classical molecular dynamics simulations. The internal energy of the clusters are representative of those present at temperatures of 200 and 300 K. Correlated electronic structure calculations are used to validate the pseudopotential for this purpose. We find that the neutral clusters support localized, bound excess electron ground states in about 50% of the configurations for the smallest cluster size studied (n = 20), and in almost all configurations for larger clusters (n > 66). The state is always exterior to the molecular frame, forming typically a diffuse surface state. Both cluster size and temperature dependence of energetic and structural properties of the clusters and the electron distribution are explored. We show that the stabilization of the electron is strongly correlated with the preexisting instantaneous dipole moment of the neutral clusters, and its ground state energy is reflected in the electronic radius. The findings are consistent with electron attachment via an initial surface state. The hypothetical spectral dynamics following such attachment is also discussed.     

Electron hydration dynamics in water clusters: A direct ab initio molecular dynamics approach

Electron attachment dynamics of excess electron in water cluster (H2O)n (n = 2 and 3) have been investigated by means of full-dimensional direct ab initio molecular dynamics (MD) method at the MP26-311++G(d,p) level. It was found that the hydrogen bond breaking due to the excess electron is an important process in the first stage of electron capture in water trimer. Time scale of electron localization and hydrogen bond breaking were determined by the direct ab initio MD simulation. The initial process of hydration in water cluster is clearly visualized in the present study. In n = 3, an excess electron is first trapped around the cyclic water trimer with a triangular form, where the excess electron is equivalently distributed on the three water molecules at time zero. After 50 fs, the excess electron is concentrated into two water molecules, while the potential energy of the system decreases by -1.5 kcal/mol from the vertical point. After 100 fs, the excess electron is localized in one of the water molecules and the potential energy decreases by -5.3 kcal/mol, but the triangular form still remained. After that, one of the hydrogen bonds in the triangular form is gradually broken by the excess electron, while the structure becomes linear at 100-300 fs after electron capture. The time scale of hydrogen bond breaking due to the excess electron is calculated to be about 300 fs. Finally, a dipole bound state is formed by the linear form of three water molecules. In the case of n = 2, the dipole bound anion is formed directly. The mechanism of electron hydration dynamics was discussed on the basis of theoretical results.     

Reaction dynamics following electron capture of chlorofluorocarbon adsorbed on water cluster: a direct density functional theory molecular dynamics study

The electron capture dynamics of halocarbon and its water complex have been investigated by means of the full dimensional direct density functional theory molecular dynamics method in order to shed light on the mechanism of electron capture of a halocarbon adsorbed on the ice surface. The CF(2)Cl(2) molecule and a cyclic water trimer (H(2)O)(3) were used as halocarbon and water cluster, respectively. The dynamics calculation of CF(2)Cl(2) showed that both C-Cl bonds are largely elongated after the electron capture, while one of the Cl atoms is dissociated from CF(2)Cl(2) (-) as a Cl(-) ion. Almost all total available energy was transferred into the internal modes of the parent CF(2)Cl radical on the product state, while the relative translational energy of Cl(-) was significantly low due to the elongation of two C-Cl bonds. In the case of a halocarbon-water cluster system, the geometry optimization of neutral complex CF(2)Cl(2)(H(2)O)(3) showed that one of the Cl atoms interacts with n orbital of water molecules of trimer and the other Cl atom existed as a dangling Cl atom. After the electron capture, only one C-Cl bond (dangling Cl atom) was rapidly elongated, whereas the other C-Cl bond is silent during the reaction. The dangling Cl atom was directly dissociated from CF(2)Cl(2) (-)(H(2)O)(3) as Cl(-). The fast Cl(-) ion was generated from CF(2)Cl(2) (-)(H(2)O)(3) on the water cluster. The mechanism of the electron capture of halocarbon on water ice was discussed on the basis of the theoretical results.     

Molecular dynamics simulations of the Ag+ or Na+ cation with an excess electron in bulk water

The properties of an excess electron interacting with a monovalent cation in bulk water are studied by molecular dynamics simulations. Sodium and silver cations are chosen as prototypical cases because of their very different redox properties. In both cases, mixed quantum classical molecular dynamics simulations reproduce the experimental UV-Vis spectra. In the case of silver, we observe a highly polarized neutral atom, corresponding to a dipolar excitonic state. For sodium a contact cation/electron pair is observed. Free energy curves along the cation electron coordinate are calculated using quantum Umbrella Sampling technique. The relative stability of the different chemical species is discussed.     

Dissociative electron capture of halocarbon caused by the internal electron transfer from water trimer anion

Dissociative electron capture dynamics of halocarbon absorbed on water cluster anion, caused by internal electron transfer from the water trimer anion to the halocarbon, have been investigated by means of the direct density functional theory (DFT)-molecular dynamics (MD) method. The CF(2)Cl(2) molecule and a water trimer anion e(-)(H(2)O)(3) were used as a halocarbon and a trapped electron, respectively. First, the structure of trapped electron state, expressed by e(-)(H(2)O)(3)-CF(2)Cl(2), was fully optimized. The excess electron was trapped by a dipole moment of water trimer. Next, initial geometries were randomly generated around the equilibrium point of the trapped electron state, and then trajectories were run. The direct DFT-MD calculations showed that the spin density distribution of excess electron is gradually changed from the water cluster (trapped electron state) to CF(2)Cl(2) as a function of time. Immediately, the Cl(-) ion was dissociated from CF(2)Cl(2)(-) adsorbed on the water cluster. The reaction was schematically expressed by e(-)(H(2)O)(3)-CF(2)Cl(2)-->[(H(2)O)(3)-->-CF(2)CL(2)](-) --> (H(2O)(3) + CF(2)CL + CI(-) (I) where [(H(2)O)(3)-CF(2)Cl(2)](-) indicates a transient intermediate state in which the excess electron is widely distributed on both the water cluster and CF(2)Cl(2). The mechanism of the electron capture of halocarbon from the trapped electron in water ice was discussed on the basis of the theoretical results. Also, the dynamics feature was compared with those of the direct electron capture reactions of CF(2)Cl(2) and CF(2)Cl(2)-(H(2)O)(3), i.e. e(-) + CF(2)Cl(2), and e(-) + CF(2)Cl(2)-(H(2)O)(3), investigated in our previous paper [Tachikawa and Abe, J. Chem. Phys., 2007, 126, 194310].     

Analysis of localization sites for an excess electron in neutral methanol clusters using approximate pseudopotential quantum-mechanical calculations

We have used a recently developed electron-methanol molecule pseudopotential in approximate quantum mechanical calculations to evaluate and statistically analyze the physical properties of an excess electron in the field of equilibrated neutral methanol clusters ((CH(3)OH)(n), n=50-500). The methanol clusters were generated in classical molecular dynamics simulations at nominal 100 and 200 K temperatures. Topological analysis of the neutral clusters indicates that methyl groups cover the surface of the clusters almost exclusively, while the associated hydroxyl groups point inside. Since the initial neutral clusters are lacking polarity on the surface and compact inside, the excess electron can barely attach to these structures. Nevertheless, most of the investigated cluster configurations do support weakly stabilized cluster anion states. We find that similarly to water clusters, the pre-existing instantaneous dipole moment of the neutral clusters binds the electron. The localizing electrons occupy diffuse, weakly bound surface states that largely engulf the cluster although their centers are located outside the cluster molecular frame. The initial localization of the excess electron is reflected in its larger radius compared to water due to the lack of free OH hydrogens on the cluster surface. The stabilization of the excess electron increases, while the radius decreases monotonically as the clusters grow in size. Stable, interior bound states of the excess electron are not observed to form neither in finite size methanol clusters nor in the equilibrium bulk.     

How important are temperature effects for cluster polarizabilities?

State-of-the-art first-principle all-electron density functional theory calculations on small sodium clusters are performed to study the temperature dependency of their polarizabilities. For this purpose Born-Oppenheimer molecular dynamics simulations with more than 100,000 time steps (>200 ps) are recorded employing gradient corrected functionals in combination with a double-zeta valence polarization basis set. For each cluster 18 trajectories between 50 and 900 K are collected. The cluster polarizabilities are then calculated along these trajectories employing a triple-zeta valence polarization basis set augmented with field-induced polarization functions. The analysis of these calculations shows that the temperature dependency of the sodium cluster polarizabilities varies strongly with cluster size. For several clusters characteristic changes in the polarizability per atom as a function of temperature are observed. It is shown that the inclusion of finite temperature effects resolves the long-standing mismatch between calculated and measured sodium cluster polarizabilities.     

A first principles molecular dynamics study of excess electron and lithium atom solvation in water-ammonia mixed clusters: structural, spectral, and dynamical behaviors of [(H2O)5NH3]- and Li(H2O)5NH3 at finite temperature

First principles molecular dynamics simulations are carried out to investigate the solvation of an excess electron and a lithium atom in mixed water-ammonia cluster (H(2)O)(5)NH(3) at a finite temperature of 150 K. Both [(H(2)O)(5)NH(3)](-) and Li(H(2)O)(5)NH(3) clusters are seen to display substantial hydrogen bond dynamics due to thermal motion leading to many different isomeric structures. Also, the structures of these two clusters are found to be very different from each other and also very different from the corresponding neutral cluster without any excess electron or the metal atom. Spontaneous ionization of Li atom occurs in the case of Li(H(2)O)(5)NH(3). The spatial distribution of the singly occupied molecular orbital shows where and how the excess (or free) electron is primarily localized in these clusters. The populations of single acceptor (A), double acceptor (AA), and free (NIL) type water and ammonia molecules are found to be significantly high. The dangling hydrogens of these type of water or ammonia molecules are found to primarily capture the free electron. It is also found that the free electron binding motifs evolve with time due to thermal fluctuations and the vertical detachment energy of [(H(2)O)(5)NH(3)](-) and vertical ionization energy of Li(H(2)O)(5)NH(3) also change with time along the simulation trajectories. Assignments of the observed peaks in the vibrational power spectra are done and we found a one to one correlation between the time-averaged populations of water and ammonia molecules at different H-bonding sites with the various peaks of power spectra. The frequency-time correlation functions of OH stretch vibrational frequencies of these clusters are also calculated and their decay profiles are analyzed in terms of the dynamics of hydrogen bonded and dangling OH modes. It is found that the hydrogen bond lifetimes in these clusters are almost five to six times longer than that of pure liquid water at room temperature.     

Quantum-classical simulation of electron localization in negatively charged methanol clusters

A series of quantum molecular dynamics simulations have been performed to investigate the energetic, structural, dynamic, and spectroscopic properties of methanol cluster anions, [(CH(3)OH)(n)](-), (n = 50-500). Consistent with the inference from photo-electron imaging experiments, we find two main localization modes of the excess electron in equilibrated methanol clusters at ～200 K. The two different localization patterns have strikingly different physical properties, consistent with experimental observations, and are manifest in comparable cluster sizes to those observed. Smaller clusters (n ≤ 128) tend to localize the electron in very weakly bound, diffuse electronic states on the surface of the cluster, while in larger ones the electron is stabilized in solvent cavities, in compact interior-bound states. The interior states exhibit properties that largely resemble and smoothly extrapolate to those simulated for a solvated electron in bulk methanol. The surface electronic states of methanol cluster anions are significantly more weakly bound than the surface states of the anionic water clusters. The key source of the difference is the lack of stabilizing free hydroxyl groups on a relaxed methanol cluster surface. We also provide a mechanistic picture that illustrates the essential role of the interactions of the excess electron with the hydroxyl groups in the dynamic process of the transition of the electron from surface-bound states to interior-bound states.     

Thermal effects on energetics and dynamics in water cluster anions (H2O)n(-)

The electron binding energies and relaxation dynamics of water cluster anions (H(2)O)(n)(-) (11 ≤ n ≤ 80) formed in co-expansions with neon were investigated using one-photon and time-resolved photoelectron imaging. Unlike previous experiments with argon, water cluster anions exhibit only one isomer class, the tightly bound isomer I with approximately the same binding energy as clusters formed in argon. This result, along with a decrease in the internal conversion lifetime of excited (H(2)O)(n)(-) (25 ≤ n ≤ 40), indicates that clusters are vibrationally warmer when formed in neon. Over the ranges studied, the vertical detachment energies and lifetimes appear to converge to previously reported values.     

Ab initio molecular dynamics studies of ionic dissolution and precipitation of sodium chloride and silver chloride in water clusters, NaCl(H2O)n and AgCl(H2O)n, n = 6, 10, and 14

An ab initio molecular dynamics method was used to compare the ionic dissolution of soluble sodium chloride (NaCl) in water clusters with the highly insoluble silver chloride (AgCl). The investigations focused on the solvation structures, dynamics, and energetics of the contact ion pair (CIP) and of the solvent-separated ion pair (SSIP) in NaCl(H(2)O)(n) and AgCl(H(2)O)(n) with cluster sizes of n = 6, 10 and 14. We found that the minimum cluster size required to stabilize the SSIP configuration in NaCl(H(2)O)(n) is temperature-dependent. For n = 6, both configurations are present as two distinct local minima on the free-energy profile at 100 K, whereas SSIP is unstable at 300 K. Both configurations, separated by a low barrier (<10 kJ mol(-1)), are identifiable on the free energy profiles of NaCl(H(2)O)(n) for n = 10 and 14 at 300 K, with the Na(+)/Cl(-) pairs being internally solvated in the water cluster and the SSIP configuration being slightly higher in energy (<5 kJ mol(-1)). In agreement with the low bulk solubility of AgCl, no SSIP minimum is observed on the free-energy profiles of finite AgCl(H(2)O)(n) clusters. The AgCl interaction is more covalent in nature, and is less affected by the water solvent. Unlike NaCl, AgCl is mainly solvated on the surface in finite water clusters, and ionic dissolution requires a significant reorganization of the solvent structure.     

Dynamics of electron solvation in I(-)(CH3OH)n clusters (4 ≤ n ≤ 11)

The dynamics of electron solvation following excitation of the charge-transfer-to-solvent precursor state in iodide-doped methanol clusters, I(-)(CH(3)OH)(n = 4-11), are studied with time-resolved photoelectron imaging. This excitation produces a I···(CH(3)OH)(n)(-) cluster that is unstable with respect to electron autodetachment and whose autodetachment lifetime increases monotonically from ~800 fs to 85 ps as n increases from 4 to 11. The vertical detachment energy (VDE) and width of the excited state feature in the photoelectron spectrum show complex time dependence during the lifetime of this state. The VDE decreases over the first 100-400 fs, then rises exponentially to a maximum with a ~1 ps time constant, and finally decreases by as much as 180 meV with timescales of 3-20 ps. The early dynamics are associated with electron transfer from the iodide to the methanol cluster, while the longer-time changes in VDE are attributed to solvent reordering, possibly in conjunction with ejection of neutral iodine from the cluster. Changes in the observed width of the spectrum largely follow those of the VDEs; the dynamics of both are attributed to the major rearrangement of the solvent cluster during relaxation. The relaxation dynamics are interpreted as a reorientation of at least one methanol molecule and the disruption and formation of the solvent network in order to accommodate the excess charge.     

Statistical mechanically averaged molecular properties of liquid water calculated using the combined coupled cluster/molecular dynamics method

Liquid water is investigated theoretically using combined molecular dynamics (MD) simulations and accurate electronic structure methods. The statistical mechanically averaged molecular properties of liquid water are calculated using the combined coupled cluster/molecular mechanics (CC/MM) method for a large number of configurations generated from MD simulations. The method includes electron correlation effects at the coupled cluster singles and doubles level and the use of a large correlation consistent basis set. A polarizable force field has been used for the molecular dynamics part in both the CC/MM method and in the MD simulation. We describe how the methodology can be optimized with respect to computational costs while maintaining the quality of the results. Using the optimized method we study the energetic properties including the heat of vaporization and electronic excitation energies as well as electric dipole and quadrupole moments, the frequency dependent electric (dipole) polarizability, and electric-field-induced second harmonic generation first and second hyperpolarizabilities. Comparisons with experiments are performed where reliable data are available. Furthermore, we discuss the important issue on how to compare the calculated microscopic nonlocal properties to the experimental macroscopic measurements.     

Electron attachment to molecules in a cluster environment

Low-energy dissociative electron attachment (DEA) to the CF(2)Cl(2) and CF(3)Cl molecules in a water cluster environment is investigated theoretically. Calculations are performed for the water trimer and water hexamer. It is shown that the DEA cross section is strongly enhanced when the attaching molecule is embedded in a water cluster, and that this cross section grows as the number of water molecules in the cluster increases. This growth is explained by a trapping effect that is due to multiple scattering by water molecules while the electron is trapped in the cluster environment. The trapping increases the resonance lifetime and the negative ion survival probability. This confirms qualitatively existing experiments on electron attachment to the CF(2)Cl(2) molecule placed on the surface of H(2)O ice. The DEA cross sections are shown to be very sensitive to the position of the attaching molecule within the cluster and the orientation of the electron beam relative to the cluster.     

Charge-transfer-to-solvent-driven dissolution dynamics of I- (H2O)2-5 upon excitation: excited-state ab initio molecular dynamics simulations

In contrast to the extensive theoretical investigation of the solvation phenomena, the dissolution phenomena have hardly been investigated theoretically. Upon the excitation of hydrated halides, which are important substances in atmospheric chemistry, an excess electron transfers from the anionic precursor (halide anion) to the solvent and is stabilized by the water cluster. This results in the dissociation of hydrated halides into halide radicals and electron-water clusters. Here we demonstrate the charge-transfer-to-solvent (CTTS)-driven femtosecond-scale dissolution dynamics for I-(H2O)n=2-5 clusters using excited state (ES) ab initio molecular dynamics (AIMD) simulations employing the complete-active-space self-consistent-field (CASSCF) method. This study shows that after the iodine radical is released from I-(H2O)n=2-5, a simple population decay is observed for small clusters (2 相似文献   

Unique Dynamics of Water Confined on the Surface of Metal Oxide Clusters

Water confined on metal oxide surface plays significant roles in heterogeneous catalysis. Heteropolyacid, a 1.2 nm-metal oxide cluster with well-defined structure, is applied as a model to understand the dynamics of water on its surface. The surface water strongly associates with heteropolyacid cluster and form the so-called ‘pseudoliquid phase’ where catalytic reactions are conducted. Broadband dielectric spectroscopy and differential scanning calorimetry have been applied to probe the dynamics of water in this pseudoliquid phase. A supercooling phase transition of water below its normal melting temperature and a dipolar glassy relaxation behaviour due to the hindered dynamics of water have been observed. The rich dynamic behavior on the surface of such well-defined metal clusters provide new perspectives to understand the properties of surface water and their relation to catalytic performance of heteropolyacid.