Fragmentation dynamics of ionized neon clusters (Ne(n), n=3-14) embedded in helium nanodroplets

We report a theoretical study of the nonadiabatic fragmentation dynamics of ionized neon clusters embedded in helium nanodroplets for cluster sizes up to n=14 atoms. The dynamics of the neon atoms is modeled using the molecular dynamics with quantum transitions method of Tully [J. Chem. Phys. 93, 1061 (1990)] with the nuclei treated classically and transitions between electronic states quantum mechanically. The potential-energy surfaces are derived from a diatomics-in-molecules model to which induced dipole-induced dipole interactions are added. The effect of the spin-orbit interaction is also discussed. The helium environment is modeled by a friction force acting on charged atoms whose speed exceeds the critical Landau velocity. The dependence of the fragment size distribution on the friction strength and on the initial nanodroplet size is investigated. By comparing with the available experimental data obtained for Ne3+ and Ne4+, a reasonable value for the friction coefficient, the only parameter of the model, is deduced. This value is then used to predict the effect of the helium environment on the dissociation dynamics of larger neon clusters, n=5-14. The results show stabilization of larger fragments than in the gas phase, but fragmentation is not completely caged. In addition, two types of dynamics are characterized for Ne4+: fast and explosive, therefore leaving no time for friction to cool down the process when dynamics starts on one of the highest electronic states, and slower, therefore leading to some stabilization by helium when it starts on one of the lowest electronic states.     

Modelization of the fragmentation dynamics of krypton clusters (Kr(n),n=2-11) following electron impact ionization

We present the first prediction for the fragmentation dynamics following electron impact ionization of neutral krypton clusters from 2 to 11 atoms. Fragment proportions and parent ion lifetimes are deduced from a molecular dynamics with quantum transitions study in which the nuclei are treated classically and the transitions between electronic states quantum mechanically. The potential-energy surfaces are derived from a diatomics-in-molecules model to which induced dipole-induced dipole and spin-orbit interactions are added. The results show surprisingly fast and extensive fragmentation for clusters of such a heavy atom, although not as extensive as in the case of neon clusters studied previously [D. Bonhommeau et al., J. Chem. Phys. 123, 54316 (2005)]. The parent ion lifetimes range from 2.8 to 0.7 ps, and the most abundant fragment is Kr(2) (+) for all studied sizes, followed by Kr(+) for sizes smaller than 7 atoms and by Kr(3) (+) for larger sizes. Trimer and larger fragments are found to originate from the lower electronic states of parent ions. The comparison with preliminary results from experiments on size-selected neutral clusters conducted by Steinbach et al. (private communication) reveal a good agreement on the extensive character of the fragmentation. It is checked that the additional internal energy brought by the helium scattering technique used for size selection does not affect the fragment proportions. In addition, the existence of long-lived trajectories is revealed, and they are found to be more and more important for larger cluster sizes and to favor the stabilization of larger fragments. The implications of this work for microsecond-scale dynamics of ionized rare-gas clusters are discussed. In particular, given the extent of fragmentation of the parent clusters and the fast kinetics of the whole process, the small cluster ions that exhibit a monomer loss in the microsecond time window must originate from much larger neutral precursors. The decay rate of the II(12)(u) state of the ionic dimer Kr(2) (+) by spin-orbit coupling is found to be of the order of 3 ps, in contrast to the expected tens of microseconds, but only reasonably faster than the corresponding state of HeNe(+). Finally, the spin-orbit interaction strongly affects both the Kr(+)Kr(2) (+) ratio and some of the characteristic times of the dynamics, especially for smaller sizes, but not the overall dependence of the fragment proportions as a function of cluster size.     

Fragmentation dynamics of argon clusters (Ar(n), n = 2 to 11) following electron-impact ionization: modeling and comparison with experiment

The fragmentation dynamics of argon clusters ionized by electron impact is investigated for initial cluster sizes up to n = 11 atoms. The dynamics of the argon atoms is modeled using a mixed quantum-classical method in which the nuclei are treated classically and the transitions between electronic states quantum mechanically. The potential-energy surfaces are derived from a diatomics-in-molecules model with the addition of the induced dipole-induced dipole and spin-orbit interactions. The results show extensive and fast fragmentation. The dimer is the most abundant ionic fragment, with a proportion increasing from 66% for n = 2 to a maximum of 95% for n = 6 and then decreasing down to 67% for n = 11. The next abundant fragment is the monomer for n < 7 and the trimer otherwise. The parent ion dissociation lifetimes are all in the range of 1 ps. Long-lived trajectories appear for initial cluster sizes of seven and higher, and favor the formation of the larger fragments (trimers and tetramers). Our results show quantitative agreement with available experimental results concerning the extensive character of the fragmentation: Ar+ and Ar2(+) are the only ionic fragments for sizes up to five atoms; their overall proportion is in quantitative agreement for all the studied sizes; Ar2(+) is the main fragment for all sizes; stable Ar3(+) fragments only appear for n > or = 5, and their proportion increases smoothly with cluster size from there. However, the individual ionic monomer and dimer fragment proportions differ. The experimental ones exhibit oscillations with initial cluster size, with a slight tendency to decrease on average for the monomer. In contrast our results show a monotonic, systematic evolution, similar to what was found in our earlier studies on neon and krypton clusters. Several hypotheses are discussed in order to find the origin of this discrepancy. Finally, the metastable II(1/2)u and II(1/2)g states of Ar2(+) are found to decay with a lifetime of 3.5 and 0.1 ps, respectively, due to spin-orbit coupling. The difference with the commonly accepted microsecond range value for rare-gas dimer ions could originate from the role of autoionizing states in the formation of the parent ions.     

Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations

The non-relativistic quantum dynamics of nuclei and electrons is solved within the framework of quantum hydrodynamics using the adiabatic representation of the electronic states. An on-the-fly trajectory-based nonadiabatic molecular dynamics algorithm is derived, which is also able to capture nuclear quantum effects that are missing in the traditional trajectory surface hopping approach based on the independent trajectory approximation. The use of correlated trajectories produces quantum dynamics, which is in principle exact and computationally very efficient. The method is first tested on a series of model potentials and then applied to study the molecular collision of H with H(2) using on-the-fly TDDFT potential energy surfaces and nonadiabatic coupling vectors.     

Fragmentation of ionized doped helium nanodroplets: theoretical evidence for a dopant ejection mechanism

We report a theoretical study of the effect induced by a helium nanodroplet environment on the fragmentation dynamics of a dopant. The dopant is an ionized neon cluster Ne(n) (+) (n=4-6) surrounded by a helium nanodroplet composed of 100 atoms. A newly designed mixed quantum/classical approach is used to take into account both the large helium cluster zero-point energy due to the light mass of the helium atoms and all the nonadiabatic couplings between the Ne(n) (+) potential-energy surfaces. The results reveal that the intermediate ionic dopant can be ejected from the droplet, possibly with some helium atoms still attached, thereby reducing the cooling power of the droplet. Energy relaxation by helium atom evaporation and dissociation, the other mechanism which has been used in most interpretations of doped helium cluster dynamics, also exhibits new features. The kinetic energy distribution of the neutral monomer fragments can be fitted to the sum of two Boltzmann distributions, one with a low kinetic energy and the other with a higher kinetic energy. This indicates that cooling by helium atom evaporation is more efficient than was believed so far, as suggested by recent experiments. The results also reveal the predominance of Ne(2) (+) and He(q)Ne(2) (+) fragments and the absence of bare Ne(+) fragments, in agreement with available experimental data (obtained for larger helium nanodroplets). Moreover, the abundance in fragments with a trimeric neon core is found to increase with the increase in dopant size. Most of the fragmentation is achieved within 10 ps and the only subsequent dynamical process is the relaxation of hot intermediate He(q)Ne(2) (+) species to Ne(2) (+) by helium atom evaporation. The dependence of the ionic fragment distribution on the parent ion electronic state reached by ionization is also investigated. It reveals that He(q)Ne(+) fragments are produced only from the highest electronic state, whereas He(q)Ne(2) (+) fragments originate from all the electronic states. Surprisingly, the highest electronic states also lead to fragments that still contain the original ionic dopant species. A mechanism is conjectured to explain this fragmentation inhibition.     

Ultraviolet-visible absorption of small silver clusters in neon: Ag(n) (n = 1-9)

We present optical absorption and fluorescence spectra in the UV-visible range of size selected neutral Ag(n) clusters (n = 1-9) in solid neon. Rich and detailed optical spectra are found with linewidths as small as 50 meV. These spectra are compared to time dependent density functional theory implemented in the TURBOMOLE package. Excellent agreement between theory and experiment is achieved in particular for the dominant spectroscopic features at photon energies below 4.5 eV. This allows a clear attribution of the observed electronic transitions to specific isomers. Optical transitions associated to the s-electrons are concentrated in the energy range between 3 and 4 eV and well separated from transitions of the d-electrons. This is in contrast to the other coinage metals (Au and Cu) which show a strong coupling of the d-electrons.     

Optical absorption of small copper clusters in neon: Cu(n), (n = 1-9)

We present optical absorption spectra in the UV-visible range (1.6 eV < ?ω < 5.5 eV) of mass selected neutral copper clusters Cu(n)(n = 1-9) embedded in a solid neon matrix at 7 K. The atom and the dimer have already been measured in neon matrices, while the absorption spectra for sizes between Cu(3) and Cu(9) are entirely (n = 6-9) or in great part new. They show a higher complexity and a larger number of transitions distributed over the whole energy range compared to similar sizes of silver clusters. The experimental spectra are compared to the time dependent density functional theory (TD-DFT) implemented in the TURBOMOLE package. The analysis indicates that for energies larger than 3 eV the transitions are mainly issued from d-type states; however, the TD-DFT scheme does not reproduce well the detailed structure of the absorption spectra. Below 3 eV the agreement for transitions issued from s-type states is better.     

Modeling the fragmentation dynamics of ionic clusters inside helium nanodroplets: the case of He100Ne4+

We present simulation results on the effect of a helium nanodroplet environment on the fragmentation dynamics of embedded molecular systems. The helium atoms are treated explicitly, with zero-point effects taken into account through an effective helium-helium interaction potential. The ionized neon tetramer is used as a model molecular system because, like all the small rare-gas clusters, it fragments extensively upon ionization. All the nonadiabatic effects between electronic states of the ionized neon cluster are taken into account. The results reveal a predominance of Ne2+ and HepNe2+ fragments and the absence of bare Ne+ fragments, in agreement with available experimental data. The neutral monomer fragments exhibit a rather wide kinetic energy distribution that can be fitted to the sum of two Boltzmann distributions, one with a low kinetic energy and the other with a higher kinetic energy. This indicates that cooling by helium atom evaporation is more efficient than was believed so far, as suggested by recent experimental results. Purely classical calculations are shown to strongly overestimate the amount of cage effect (cooling), clearly indicating the need to take into account zero-point effects.     

On-the-fly, electric-field-driven, coupled electron-nuclear dynamics

An on-the-fly, electric field driven, coupled electron-nuclear dynamics approach is developed and applied to model the photodissociation of water in the A((1)B1) excited state. In this method, a quantum propagator evolves the photon-induced electronic dynamics in the ultrafast time scale, and a quasi-classical surface hopping approach describes the nuclear dynamics in the slower time scale. In addition, strong system-field interactions are explicitly included in the electronic propagator. This theoretical development enables us to study rapid photon-induced bond dissociation dynamics and demonstrates the partial breakdown of electronic coherence as well as electronic population trapping in the excited state when the molecular vibrations detune the system with respect to the applied field. The method offers a practical way to use on-the-fly dynamics for modeling light-molecule interactions that lead to interesting photochemical events.     

Chromium-doped germanium clusters CrGen (n = 1-5): geometry, electronic structure, and topology of chemical bonding

The structure and properties of small neutral and cationic CrGen(0,+) clusters, with n from 1 to 5, were investigated using quantum chemical calculations at the CASSCF/CASPT2 and DFT/B3LYP levels. Smaller clusters prefer planar geometries, whereas the lowest-lying electronic states of the neutral CrGe4, CrGe5, and cationic CrGe5+ forms exhibit nonplanar geometries. Most of the clusters considered prefer structures with high-spin ground state and large magnetic moments. Relative to the values obtained for the pure Gen clusters, fragmentation energies of doped CrGen clusters are smaller when n is 3 and 4 and larger when n = 5. The averaged binding energy tends to increase with the increasing number of Ge atoms. For n = 5, the binding energies for Ge5, CrGe5, and CrGe5+ are similar to each other, amounting to approximately 2.5 eV. The Cr atom acts as a general electron donor in neutral CrGen clusters. Electron localization function (ELF) analyses suggest that the chemical bonding in chromium-doped germanium clusters differs from that of their pure or Li-doped counterparts and allow the origin of the inherent high-spin ground state to be understood. The differential DeltaELF picture, obtained in separating both alpha and beta electron components, is consistent with that derived from spin density calculations. For CrGen, n = 2 and 3, a small amount of d-pi back-donation is anticipated within the framework of the proposed bonding model.     

Structural and electronic properties of Si n, Si n+, and AlSi n-1 (n=2-13) clusters: theoretical investigation based on ab initio molecular orbital theory

The geometric and electronic structures of Si(n), Si(n) (+), and AlSi(n-1) clusters (2< or =n< or =13) have been investigated using the ab initio molecular orbital theory under the density functional theory formalism. The hybrid exchange-correlation energy function (B3LYP) and a standard split-valence basis set with polarization functions [6-31G(d)] were employed for this purpose. Relative stabilities of these clusters have been analyzed based on their binding energies, second difference in energy (Delta (2)E) and fragmentation behavior. The equilibrium geometry of the neutral and charged Si(n) clusters show similar structural growth. However, significant differences have been observed in the electronic structure leading to their different stability pattern. While for neutral clusters, the Si(10) is magic, the extra stability of the Si(11) (+) cluster over the Si(10) (+) and Si(12) (+) bears evidence for the magic behavior of the Si(11) (+) cluster, which is in excellent agreement with the recent experimental observations. Similarly for AlSi(n-1) clusters, which is isoelectronic with Si(n) (+) clusters show extra stability of the AlSi(10) cluster suggesting the influence of the electronic structures for different stabilities between neutral and charged clusters. The ground state geometries of the AlSi(n-1) clusters show that the impurity Al atom prefers to substitute for the Si atom, that has the highest coordination number in the host Si(n) cluster. The fragmentation behavior of all these clusters show that while small clusters prefers to evaporate monomer, the larger ones dissociate into two stable clusters of smaller size.     

Hydride transfer in liver alcohol dehydrogenase: quantum dynamics, kinetic isotope effects, and role of enzyme motion.

The quantum dynamics of the hydride transfer reaction catalyzed by liver alcohol dehydrogenase (LADH) are studied with real-time dynamical simulations including the motion of the entire solvated enzyme. The electronic quantum effects are incorporated with an empirical valence bond potential, and the nuclear quantum effects of the transferring hydrogen are incorporated with a mixed quantum/classical molecular dynamics method in which the transferring hydrogen nucleus is represented by a three-dimensional vibrational wave function. The equilibrium transition state theory rate constants are determined from the adiabatic quantum free energy profiles, which include the free energy of the zero point motion for the transferring nucleus. The nonequilibrium dynamical effects are determined by calculating the transmission coefficients with a reactive flux scheme based on real-time molecular dynamics with quantum transitions (MDQT) surface hopping trajectories. The values of nearly unity for these transmission coefficients imply that nonequilibrium dynamical effects such as barrier recrossings are not dominant for this reaction. The calculated deuterium and tritium kinetic isotope effects for the overall rate agree with experimental results. These simulations elucidate the fundamental nature of the nuclear quantum effects and provide evidence of hydrogen tunneling in the direction along the donor-acceptor axis. An analysis of the geometrical parameters during the equilibrium and nonequilibrium simulations provides insight into the relation between specific enzyme motions and enzyme activity. The donor-acceptor distance, the catalytic zinc-substrate oxygen distance, and the coenzyme (NAD(+)/NADH) ring angles are found to strongly impact the activation free energy barrier, while the donor-acceptor distance and one of the coenzyme ring angles are found to be correlated to the degree of barrier recrossing. The distance between VAL-203 and the reactive center is found to significantly impact the activation free energy but not the degree of barrier recrossing. This result indicates that the experimentally observed effect of mutating VAL-203 on the enzyme activity is due to the alteration of the equilibrium free energy difference between the transition state and the reactant rather than nonequilibrium dynamical factors. The promoting motion of VAL-203 is characterized in terms of steric interactions involving THR-178 and the coenzyme.     

Geometries and electronic properties of the neutral and charged rare earth Yb-doped Si(n) (n = 1-6) clusters: a relativistic density functional investigation

The neutral and charged YbSi(n) (n = 1-6) clusters considering different spin configurations have been systematically investigated by using the relativistic density functional theory with generalized gradient approximation. The total bonding energies, equilibrium geometries, Mulliken populations (MP), Hirshfeld charges (HC), fragmentation energies, and highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gaps are calculated and discussed. The optimized geometries indicate that the most stable YbSi(n) (n = 1-6) clusters keep basically the analogous frameworks as the low-lying Si(n)(+1) clusters, while the charged species deviate from their neutral counterparts, and that the doped Yb tends to occupy the substitutional site of the neutral and charged YbSi(n) isomers. The relative stabilities are investigated in terms of the calculated fragmentation energies, exhibiting enhanced stabilities for the remarkably stable neutral and charged YbSi2 and YbSi5 clusters. Furthermore, the calculated MP and HC values show that the charges of the neutral and charged YbSi(n) clusters transfer from the Yb atom to Si(n) atoms and the Yb atom acts as an electron donor, and that the f orbitals of the Yb atom in the neutral and charged YbSi(n) clusters behave as core without involvement in chemical bonding. The calculated HOMO-LUMO gaps indicate that the YbSi2 and YbSi4+ clusters have stronger chemical stabilities. Comparisons of the Yb-doped Si(n) (n = 1-6) with available theoretical results of transition-metal-doped silicon clusters are made. The growth pattern is investigated also.     

One-Electron Pseudo-Potential Determination of Stable Isomers of the Li*Ar n Excited Clusters: Absorption Spectra

We present pseudo-potential calculations of geometrical structures of stable isomers of LiAr _n clusters with both an electronic ground state and excited states of the lithium atom. The Li atom is perturbed by argon atoms in LiAr _n clusters. Its electronic structure obtained as the eigenfunctions of a single-electron operator describing the electron in the field of a Li⁺Ar _n core, the Li⁺ and Ar atoms are replaced by pseudo-potentials. These pseudo-potentials include core-polarization operators to account for the polarization and correlation of the inert core with the valence Lithium electron [J Chem Phys 116, 1839 1]. The geometry optimization of the ground and excited states of LiAr _n (n = 1–12) clusters is carried out via the Basin-Hopping method of Wales et al. [J Phys Chem 101, 5111 2; J Chem Phys 285, 1368 3]. The geometries of the ground and ionic states of LiAr _n clusters were used to determine the energy of the high excited states of the neutral LiAr _n clusters. The variation of the excited state energies of LiAr _n clusters as a function of the number of argon atoms shows an approximate Rydberg character, corresponding to the picture of an excited electron surrounding an ionic cluster core, is already reached for the 3s state. The result of optical transitions calculations shows that the absorption spectral features are sensitive to isomer structure. It is clearly the case for transitions close to the 2p levels of Li which are distorted by the cluster environment.     

Quantum fluctuation of electronic wave-packet dynamics coupled with classical nuclear motions

An ab initio electronic wave-packet dynamics coupled with the simultaneous classical dynamics of nuclear motions in a molecule is studied. We first survey the dynamical equations of motion for the individual components. Reflecting the nonadiabatic dynamics that electrons can respond to nuclear motions only with a finite speed, the equations of motion for nuclei include a force arising from the kinematic (nuclear momentum) coupling from electron cloud. To materialize these quantum effects in the actual ab initio calculations, we study practical implementation of relevant electronic matrix elements that are related to the derivatives with respect to the nuclear coordinates. Applications of the present scheme are performed in terms of the configuration state functions (CSF) using the canonical molecular orbitals as basis functions without transformation to particular diabatic basis. In the CSF representation, the nonadiabatic interaction due to the kinematic coupling is anticipated to be rather small, and instead it should be well taken into account through the off-diagonal elements of the electronic Hamiltonian matrix. Therefore it is expected that the nonadiabatic dynamics based on this CSF basis neglecting the kinematic coupling may work. To verify this anticipation and to quantify the actual effects of the kinematic coupling, we compare the dynamics with and without the kinematic-coupling terms using the same CSF set. Applications up to the fifth electronically excited states in a nonadiabatic collision between H(2) and B(+) shows that the overall behaviors of these two calculations are surprisingly similar to each other in an average sense except for a fast fluctuation reflecting the electronic time scale. However, at the same time, qualitative differences in the collision events are sometimes observed. Therefore it turns out after all that the kinematic-coupling terms cannot be neglected in the CSF-basis representation. The present applications also demonstrate that the nonadiabatic electronic wave-packet dynamics within ab initio quantum chemical calculation is feasible.     

A new approach to decoherence and momentum rescaling in the surface hopping algorithm

As originally proposed, the fewest switches surface hopping (FSSH) algorithm does not allow for decoherence between wavefunction amplitudes on different adiabatic surfaces. In this paper, we propose an inexpensive correction to standard FSSH dynamics wherein we explicitly model the decoherence of nuclear wave packets on distinct electronic surfaces. Our augmented fewest switches surface hopping approach is conceptually simple and, thus far, it has allowed us to capture several key features of the exact quantum results. Two points in particular merit attention. First, we obtain the correct branching ratios when a quantum particle passes through more than one region of nonadiabatic coupling. Second, our formalism provides a new and natural approach for rescaling nuclear momenta after a surface hop. Both of these features should become increasingly important as surface hopping schemes are applied to higher-dimensional problems.     

Structure and stability of oxygen adsorption on Si(n) (n = 5-10) clusters

The structures, binding energies, and electronic properties of one oxygen atom (O) and two oxygen atoms (2O) adsorption on silicon clusters Si(n) with n ranging from 5 to 10 are studied systematically by ab initio calculations. Twelve stable structures are obtained, two of which are in agreement with those reported in previous literature and the others are new structures that have not been proposed before. Further investigations on the fragmentations of Si(n)O and Si(n)O2 (n = 5-10) clusters indicate that the pathways Si(n)O --> Si(n-1) + SiO and Si(n)O2 --> Si(n-2) + Si2O2 are most favorable from thermodynamic viewpoint. Among the studied silicon oxide clusters, Si8O, Si9O, Si5O2 and Si8O2 correspond to large adsorption energies of silicon clusters with respect to O or 2O, while Si8O, with the smallest dissociation energy, has a tendency to separate into Si7 + SiO. Using the recently developed quasi-atomic minimal-basis-orbital method, we have also calculated the unsaturated valences of the neutral Si(n) clusters. Our calculation results show that the Si atoms which have the largest unsaturated valences are more attractive to O atom. Placing O atom right around the Si atoms with the largest unsaturated valences usually leads to stable structures of the silicon oxide clusters.