Decay rates of inner-valence excitations in noble gas atoms

A Fano - algebraic diagrammatic construction - Stieltjes method has been recently developed for ab initio calculations of nonradiative decay rates [V. Averbukh and L. S. Cederbaum, J. Chem. Phys. 123, 204107 (2005)] of singly ionized states. In the present work this method is generalized for the case of electronic decay of excited states. The decay widths of autoionizing inner-valence-excited states of Ne, Ar, and Kr are calculated. Apart from the lowest excitation of Kr, they are found to be in good to excellent agreement with the experimental values. Comparison with the other theoretical studies shows that in many cases the new method performs better than the previously available techniques.     

Interatomic decay of inner-valence-excited states in clusters

In an isolated atom, excitation of an inner valence electron above the outer valence subshell leads to creation of an autoionizing state. Recently, it has been demonstrated experimentally that in a cluster, the inner-valence-excited states can decay also by an interatomic mechanism which has been called resonant interatomic Coulombic decay (RICD). Here we show that RICD is indeed the leading but not the only possible interatomic decay mode of the inner-valence excitations in clusters. Using Ne (2s-->3p) excitation in MgNe cluster as an example, we explore the possible decay mechanisms and draw conclusions on their relative importance and on the nature of the corresponding decay products.     

Interatomic Coulombic decay in mixed NeKr clusters

We report the occurrence of interatomic Coulombic decay (ICD) in mixed NeKr clusters. A well-defined feature ranging from 9 to 12 eV in kinetic energy is observed in coincidence with the Ne 2s photoelectrons. It derives from an ICD process, in which an initial Ne 2s vacancy is filled by a Ne 2p electron and an electron is emitted from a 4p level on a neighboring Kr atom. We have studied the dependence of the effect on photon energy, cluster composition, and cluster size. Interestingly, the ICD electron energy increases slightly and grows a shoulder on going from 2% to 5% Kr in the coexpansion process, which we interpret in terms of surface versus bulk effects.     

Low-lying electronic states of HBr2+

The present study describes the characterization of energy and structure of HBr(2+) in its low-lying electronic states, achieved through an extension of a new empirical method [Chem. Phys. Lett. 379, 139 (2003)] recently introduced to evaluate the interatomic interaction in the HX(2+) (X=F,Cl,Br,I) molecular dications. The method is based on identification of the main components of the interaction and their evaluation through some simple correlation formulas. Potential energy curves, given in a simple, natural, and analytical form, made possible the calculations of some important properties, such as double-photoionization energy thresholds, vibrational spacing, average lifetime, and Franck-Condon factors. The predictions, compared with data available in the literature, are of great interest for the analysis and interpretation of some new experimental results.     

Ab initio interatomic decay widths of excited states by applying Stieltjes imaging to Lanczos pseudospectra

Electronically excited states of atoms and molecules in an environment may decay in interatomic processes by transferring excess energy to neighboring species and ionizing them. The corresponding interatomic decay width is the most important characteristic of the decay allowing to calculate its efficiency and the final states' distribution. In this paper we present calculations of interatomic widths by the Fano-Stieltjes method applied to Lanczos pseudospectra, which has been previously shown to provide accurate autoionization widths in atoms and molecules. The use of Lanczos pseudospectra allows one to avoid the full diagonalization bottleneck and makes the method applicable to larger systems. We apply the present method to the calculation of interatomic decay widths in NeMg, NeAr and HCN[middle dot]Mg(n), n = 1, 2 clusters. The results are compared with widths obtained analytically and by other ab initio methods where available.     

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.     

Monte Carlo simulations of critical cluster sizes and nucleation rates of water

We have calculated the critical cluster sizes and homogeneous nucleation rates of water at temperatures and vapor densities corresponding to experiments by Wolk and Strey [J. Phys. Chem B 105, 11683 (2001)]. The calculations have been done with an expanded version of a Monte Carlo method originally developed by Vehkamaki and Ford [J. Chem. Phys. 112, 4193 (2000)]. Their method calculates the statistical growth and decay probabilities of molecular clusters. We have derived a connection between these probabilities and kinetic condensation and evaporation rates, and introduce a new way for the calculation of the work of formation of clusters. Three different interaction potential models of water have been used in the simulations. These include the unpolarizable SPC/E [J. Phys. Chem. 91, 6269 (1987)] and TIP4P [J. Chem. Phys. 79, 926 (1983)] models and a polarizable model by Guillot and Guissani [J. Chem. Phys. 114, 6720 (2001)]. We show that TIP4P produces critical cluster sizes and a temperature and vapor density dependence for the nucleation rate that agree well with the experimental data, although the magnitude of nucleation rate is constantly overestimated by a factor of 2 x 10(4). Guissani and Guillot's model is somewhat less successful, but both the TIP4P and Guillot and Guissani models are able to reproduce a much better experimental temperature dependency of the nucleation rate than the classical nucleation theory. Using SPC/E results in dramatically too small critical clusters and high nucleation rates. The water models give different average binding energies for clusters. We show that stronger binding between cluster molecules suppresses the decay probability of a cluster, while the growth probability is not affected. This explains the differences in results from different water models.     

Thermodynamics and equilibrium structure of Ne38 cluster: quantum mechanics versus classical

The equilibrium properties of classical Lennard-Jones (LJ38) versus quantum Ne38 Lennard-Jones clusters are investigated. The quantum simulations use both the path-integral Monte Carlo (PIMC) and the recently developed variational-Gaussian wave packet Monte Carlo (VGW-MC) methods. The PIMC and the classical MC simulations are implemented in the parallel tempering framework. The classical heat capacity Cv(T) curve agrees well with that of Neirotti et al. [J. Chem. Phys. 112, 10340 (2000)], although a much larger confining sphere is used in the present work. The classical Cv(T) shows a peak at about 6 K, interpreted as a solid-liquid transition, and a shoulder at approximately 4 K, attributed to a solid-solid transition involving structures from the global octahedral (Oh) minimum and the main icosahedral (C5v) minimum. The VGW method is used to locate and characterize the low energy states of Ne38, which are then further refined by PIMC calculations. Unlike the classical case, the ground state of Ne38 is a liquidlike structure. Among the several liquidlike states with energies below the two symmetric states (Oh and C5v), the lowest two exhibit strong delocalization over basins associated with at least two classical local minima. Because the symmetric structures do not play an essential role in the thermodynamics of Ne38, the quantum heat capacity is a featureless curve indicative of the absence of any structural transformations. Good agreement between the two methods, VGW and PIMC, is obtained. The present results are also consistent with the predictions by Calvo et al. [J. Chem. Phys. 114, 7312 (2001)] based on the quantum superposition method within the harmonic approximation. However, because of its approximate nature, the latter method leads to an incorrect assignment of the Ne38 ground state as well as to a significant underestimation of the heat capacity.     

Simulation of vibrational dephasing of I(2) in solid Kr using the semiclassical Liouville method

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.     

On the interatomic Coulombic decay in the Ne dimer

The interatomic Coulombic decay (ICD) in the Ne dimer is discussed in view of the recent experimental results. The ICD electron spectrum and the kinetic energy release of the Ne+ fragments resulting after Coulomb explosion of Ne2 (2+) are computed and compared to the measured ones. A very good agreement is found, confirming the dynamics predicted for this decay mechanism. The effect of the temperature on the electron spectrum is briefly investigated.     

Ground state, growth, and electronic properties of small lanthanum clusters

The DMol cluster method based on density-functional theory has been employed to study the structural stability and electronic structure of La(n) (n=2-14) clusters. The ground states have been found out for lanthanum clusters. The Jahn-Teller effect plays an important role in this process because there are many isomers near the ground state. The magnetism is not sensitive to interatomic spacing when the change of interatomic spacing is in a small range. Lanthanum clusters grow in an icosahedral pattern. The results of the mean binding energy, of the second derivative of binding energy, and of the formation energy show strong odd-even alternation and that 7- and 13-atom clusters are magic. Further, the HOMO-LUMO gap, the mean nearest bond lengths, and the mean magnetic moments suggest that the convergence to bulk is slow and it shows an oscillatory behavior for small lanthanum clusters.     

Electronic excited-state mixing in NeCl2

Ab initio calculations that explicitly include spin-orbit interactions are reported for the NeCl2 system of electronic states. A surprising curve crossing is observed for the C2v, T-shaped geometry. Away from the C2v geometry, the states mix, as expected. On the basis of these new results we propose a new mechanism for electronic energy transfer from highly vibrationally excited levels of the B electronic state of the chlorine molecule. It is proposed that as long as vibrational predissociation of NeCl2 proceeds by direct coupling of the initial state to the continuum states the Ne atom does not sample geometries that efficiently quench the Cl2 B electronic state. However, when the vibrational dynamics changes to the intramolecular vibrational relaxation regime the Ne atom becomes quite effective at coupling the Cl2 B3Pi0u+ state with a 3Pi2g state.     

An excited state paired interacting orbital method

A new method for analyzing and visualizing the molecular excited states, named "excited state paired interacting orbital (EPIO)," is proposed. The method is based both on the paired interacting orbital (PIO) proposed by Fujimoto and Fukui [J. Chem. Phys. 60, 572 (1974)] and the natural transition orbital (NTO) by Martin [J. Chem. Phys. 118, 4775 (2003)]. Within the PIO method, orbital interactions between the two fragmented molecules are represented practically only by a few pairs of fragment orbitals. The NTO method is a means of finding a compact orbital representation for the electronic transitions in the excited states. With the method, electronic transitions are expressed by a few particle-hole orbital pairs and a clear picture on the electronic transitions is obtained. EPIO method is designed to have both properties of the preceding two methods: electronic transitions in composite molecular systems can be expressed with a few pairs of EPIOs which are constructed with fragmented molecular orbitals (MOs). Excited state characters, such as charge transfer and local excitations, are analyzed by using EPIOs with their generation probabilities. Thus, the present method gives us clear information on the composition of MOs which play an important role in the molecular excitation processes, e.g., optical processes.