Dissociation of GaN2+ and AlN2+ in APT: Electronic structure and stability in strong DC field.

We investigate from a theoretical point of view the stability of AlN2+ and GaN2+ dications produced under high static electric fields like those reached in Atom Probe Tomography (APT) experiments. By means of quantum chemical calculations of the electronic structure of these molecules, we show that their stability is governed by two independent processes. On the one hand, the spin-orbit coupling allows some molecular excited states to dissociate by inter-system crossing. On the other hand, the action of the electric field lowers the potential energy barrier, which ensures the dication stability in standard conditions. We present a detailed example of field emission dynamics in the specific case of the 11Δ states for a parabolic tip, which captures the essentials of the process by means of a simplified model. We show that the dissociation dynamics of AlN2+ and GaN2+ is completely different despite the strong resemblance of their electronic structure.

Dissociation of GaN2+ and AlN2+ in APT: Analysis of experimental measurements.

The use of a tip-shaped sample for the atom probe tomography technique offers the unique opportunity to analyze the dynamics of molecular ions in strong DC fields. We investigate here the stability of AlN2+ and GaN2+ dications emitted from an Al0.25Ga0.75N sample in a joint theoretical and experimental study. Despite the strong chemical resemblance of these two molecules, we observe only stable AlN2+, while GaN2+ can only be observed as a transient species. We simulate the emission dynamics of these ions on field-perturbed potential energy surfaces obtained from quantum chemical calculations. We show that the dissociation is governed by two independent processes. For all bound states, a mechanical dissociation is induced by the distortion of the potential energy surface in the close vicinity of the emitting tip. In the specific case of GaN2+, the relatively small electric dipole of the dication in its ground 13Σ- and excited 11Δ states induces a weak coupling with the electric field so that the mechanical dissociation into Ga+ + N+ lasts for sufficient time to be observed. By contrast, the AlN2+ mechanical dissociation leads to Al2+ + N which cannot be observed as a correlated event. For some deeply bound singlet excited states, the spin-orbit coupling with lower energy triplet states gives another chance of dissociation by system inter-system crossing with specific patterns observed experimentally in a correlated time of flight map.

Electronic structure and stability of the SiO2+ dications produced in tomographic atom probe experiments.

The molecular electronic states of the SiO2+ dication have been investigated in a joint theoretical and experimental analysis. The use of a tip-shaped sample for tomographic atom probe analysis offers the unique opportunity to produce and to analyze the lifetime of some excited states of this dication. The perturbation brought by the large electric field of the polarized tip along the ion trajectory is analyzed by means of molecular dynamics simulation. For the typical electric fields used in the experiment, the lowest energy triplet states spontaneously dissociate, while the lowest energy singlet states do not. We show that the emission process leads to the formation of some excited singlet state, which dissociates by means of spin-orbit coupling with lower-energy triplet states to produce specific patterns associated with Si+ + O+ and Si2+ + O dissociation channels. These patterns are recorded and observed experimentally in a correlated time-of-flight map.

Role of a Neighbor Ion in the Fragmentation Dynamics of Covalent Molecules.

Fragmentation of molecular nitrogen dimers (N_{2})_{2} induced by collision with low energy 90 keV Ar^{9+} ions is studied to evidence the influence of a molecular environment on the fragmentation dynamics of N_{2} cations. Following the capture of three or four electrons from the dimer, the three-body N_{2}^{+}+N^{m+}+N^{n+} [with (m,n)=(1,1) or (1, 2)] fragmentation channels provide clean experimental cases where molecular fragmentation may occur in the presence of a neighbor molecular cation. The effect of the environment on the fragmentation dynamics within the dimer is investigated through the comparison of the kinetic energy release (KER) spectra for these three-body channels and for isolated N_{2}^{(m+n)+} monomer cations. The corresponding KER spectra exhibit energy shifts of the order of 10 eV, attributed to the deformation of the N^{m+}+N^{n+} potential energy curves in the presence of the neighboring N_{2}^{+} cation. The KER structures remain unchanged, indicating that the primary collision process is not significantly affected by the presence of a neighbor molecule.

Spin-orbit coupling in the dissociative excitation of alkali atoms at the surface of rare gas clusters: A theoretical study.

We analyze the role of the spin-orbit (SO) coupling in the dissociative dynamics of excited alkali atoms at the surface of small rare gas clusters. The electronic structure of the whole system is deduced from a one-electron model based on core polarization pseudo-potentials. It allows us to obtain in the same footing the energy, forces, and non-adiabatic couplings used to simulate the dynamics by means of a surface hopping method. The fine structure state population is analyzed by considering the relative magnitude of the SO coupling ξ, with respect to the spin-free potential energy. We identify three regimes of ξ-values leading to different evolution of adiabatic state population after excitation of the system in the uppermost state of the lowest np (2)P shell. For sufficiently small ξ, the final population of the J=12 atomic states, P12, grows up linearly from P12=13 at ξ = 0 after a diabatic dynamics. For large values of ξ, we observe a rather adiabatic dynamics with P12 decreasing as ξ increases. For intermediate values of ξ, the coupling is extremely efficient and a complete transfer of population is observed for the set of parameters associated to NaAr3 and NaAr4 clusters.

Spectroscopy and Dynamics of K Atoms on Argon Clusters.

We present a combined experimental and simulation study of the 4s → 4p photoexcitation of the K atom trapped at the surface of ArN clusters made of a few hundred Ar atoms. Our experimental method based on photoelectron spectroscopy allows us to firmly establish that one single K atom is trapped at the surface of the cluster. The absorption spectrum is characterized by the splitting of the atomic absorption line into two broad bands, a Π band associated with p orbitals parallel to the cluster surface and a Σ band associated with the perpendicular orientation. The spectrum is consistent with observations reported for K atoms trapped on lighter inert gas clusters, but the splitting between the Π and Σ bands is significantly larger. We show that a large amount of K atoms are transiently stuck and eventually lost by the Ar cluster, in contrast with previous observations reported for alkaline earth metal systems. The excitation in the Σ band leads systematically to the ejection of the K atom from the Ar cluster. On the contrary, excitation in the Π band leads to the formation of a bound state. In this case, the analysis of the experimental photoelectron spectrum by means of nonadiabatic molecular dynamics simulation shows that the relaxation drives the system toward a basin where the coordination of the K atom is 2.2 Ar atoms on the average, in a poorly structured surface.

Potential energy curves and spin-orbit coupling of light alkali-heavy rare gas molecules.

The potential energy curves of the X, A, and B states of alkali-rare gas diatomic molecules, MKr and MXe, are investigated for M = Li, Na, K. The molecular spin-orbit coefficients a(R)=<(2)Π(½)|H(SO)|(2)Π(½)> and b(R)=<(2)Π(-½)|H(SO)|(2)Σ(½)> are calculated as a function the interatomic distance R. We show that a(R) increases and b(R) decreases as R decreases. This effect becomes less and less important as the mass of the alkali increases. A comparison of the rovibrational properties deduced from our calculations with experimental measurements recorded for NaKr and NaXe shows the quality of the calculations.

Spectroscopic properties of alkali atoms embedded in Ar matrix.

We present a theoretical investigation of visible absorption and related luminescence of alkali atoms (Li, Na, and K) embedded in Ar matrix. We used a model based on core polarization pseudopotentials, which allows us to determine accurately the gas-to-matrix shifts of various trapping sites. The remarkable agreement between our calculated results and the experimental spectra recorded by several authors allows us to establish a clear assignment of the observed spectra, which are made of contributions from crystalline sites on the one hand, and of grain boundary sites on the other hand. Our study reveals remarkably large Stokes shifts, up to 9000 cm(-1), which could be observed experimentally to identify definitely the trapping sites.

Nonadiabatic molecular dynamics of photoexcited Li2(+) Ne(n) clusters.

We investigate the relaxation of photoexcited Li(2)(+) chromophores solvated in Ne(n) clusters (n = 2-22) by means of molecular dynamics with surface hopping. The simplicity of the electronic structure of these ideal systems is exploited to design an accurate and computationally efficient model. These systems present two series of conical intersections between the states correlated with the Li+Li(2s) and Li+Li(2p) dissociation limits of the Li(2)(+) molecule. Frank-Condon transition from the ground state to one of the three lowest excited states, hereafter indexed by ascending energy from 1 to 3, quickly drives the system toward the first series of conical intersections, which have a tremendous influence on the issue of the dynamics. The states 1 and 2, which originate in the Frank-Condon area from the degenerated nondissociative 1(2)Π(u) states of the bare Li(2)(+) molecule, relax mainly to Li+Li(2s) with a complete atomization of the clusters in the whole range of size n investigated here. The third state, which originates in the Frank-Condon area from the dissociative 1(2)Σ(u)(+) state of the bare Li(2)(+) molecule, exhibits a richer relaxation dynamics. Contrary to intuition, excitation into state 3 leads to less molecular dissociation, though the amount of energy deposited in the cluster by the excitation process is larger than for excitation into state 1 and 2. This extra amount of energy allows the system to reach the second series of conical intersections so that approximately 20% of the clusters are stabilized in the 2(2)Σ(g)(+) state potential well for cluster sizes n larger than 6.

Structure and photoabsorption properties of cationic alkali dimers solvated in neon clusters.

We present a theoretical investigation of the structure and optical absorption of M(2)(+) alkali dimers (M=Li,Na,K) solvated in Ne(n) clusters for n=1 to a few tens Ne atoms. For all these alkali, the lowest-energy isomers are obtained by aggregation of the first Ne atoms at the extremity of the alkali molecule. This particular geometry, common to other M(2)(+)-rare gas clusters, is intimately related to the shape of the electronic density of the X (2)Σ(g)(+) ground state of the bare M(2)(+) molecules. The structure of the first solvation shell presents equilateral Ne(3) and capped pentagonal Ne(6) motifs, which are characteristic of pure rare gas clusters. The size and geometry of the complete solvation shell depend on the alkali and were obtained at n=22 with a D(4h) symmetry for Li and at n=27 with a D(5h) symmetry for Na. For K, our study suggests that the closure of the first solvation shell occurs well beyond n=36. We show that the atomic arrangement of these clusters has a profound influence on their optical absorption spectrum. In particular, the XΣ transition from the X (2)Σ(g)(+) ground state to the first excited (2)Σ(u)(+) state is strongly blueshifted in the Frank-Condon area.

Stepwise hydration and evaporation of adenosine monophosphate nucleotide anions: a multiscale theoretical study.

The structure and finite-temperature properties of hydrated nucleotide anion adenosine 5'-monophosphate (AMP) have been theoretically investigated with a variety of methods. Using a polarizable version of the Amber force field and replica-exchange molecular dynamics simulations, putative lowest-energy structures have been located for the AMP(-)(H(2)O)(n) cluster anions with n = 0-20. The hydration energies obtained with the molecular mechanics potential slightly overestimate experimental measurements. However, closer values are found after reoptimizing the structures locally at more sophisticated levels, namely semi-empirical (PM6) and density-functional theory (B3LYP/6-31+G*). Upon heating the complexes, various indicators such as the heat capacity, number of hydrogen bonds or surface area provide evidence that the water cluster melts below 200 K but remains bonded to the AMP anion. The sequential loss of water molecules after sudden heating has been studied using a statistical approach in which unimolecular evaporation is described using the orbiting transition state version of phase space theory, together with anharmonic densities of vibrational states. The evaporation rates are calibrated based on the results of molecular dynamics trajectories at high internal energy. Our results indicate that between 4 and 10 water molecules are lost from AMP(-)(H(2)O)(20) after one second depending on the initial heating in the 250-350 K range, with a concomitant cooling of the remaining cluster by 75-150 K.

Accurate evaporation rates of pure and doped water clusters in vacuum: A statistico-dynamical approach.

Unimolecular evaporation of selected pure (H(2)O)(n) and heterogeneous (H(2)O)(n-1)X(+) water clusters containing a single hydronium or ammonium impurity is investigated in the framework of phase space theory (PST) in its orbiting transition state version. Using the many-body polarizable Kozack-Jordan potential and its extensions for X(+)=H(3)O(+) and NH(4) (+), the thermal evaporation of clusters containing 21 and 50 molecules is simulated at several total energies. Numerous molecular dynamics (MD) trajectories at high internal energies provide estimates of the decay rate constant, as well as the kinetic energy and angular momentum released upon dissociation. Additional Monte Carlo simulations are carried out to determine the anharmonic densities of vibrational states, which combined with suitable forms for the rotational densities of states provide expressions for the energy-resolved differential rates. Successful comparison between the MD results and the independent predictions of PST for the distributions of kinetic energy and angular momentum released shows that the latter statistical approach is quantitative. Using MD data as a reference, the absolute evaporation rates are calculated from PST over broad energy and temperature ranges. Based on these results, the presence of an ionic impurity is generally found to decrease the rate, however the effect is much more significant in the 21-molecule clusters. Our calculations also suggest that due to backbendings in the microcanonical densities of states the variations of the evaporation rates may not be strictly increasing with energy or temperature.

An accurate model potential for alkali neon systems.

We present a detailed investigation of the ground and lowest excited states of M-Ne dimers, for M=Li, Na, and K. We show that the potential energy curves of these Van der Waals dimers can be obtained accurately by considering the alkali neon systems as one-electron systems. Following previous authors, the model describes the evolution of the alkali valence electron in the combined potentials of the alkali and neon cores by means of core polarization pseudopotentials. The key parameter for an accurate model is the M(+)-Ne potential energy curve, which was obtained by means of ab initio CCSD(T) calculation using a large basis set. For each MNe dimer, a systematic comparison with ab initio computation of the potential energy curve for the X, A, and B states shows the remarkable accuracy of the model. The vibrational analysis and the comparison with existing experimental data strengthens this conclusion and allows for a precise assignment of the vibrational levels.

Structure, stability, and infrared spectroscopy of (H2O)nNH4(+) clusters: a theoretical study at zero and finite temperature.

The combined effects of size and temperature on the stable structures of water clusters doped with one ammonium molecule have been investigated theoretically using an empirical potential and density-functional theory (DFT) calculations. Global optimization with Monte Carlo methods has been performed using an explicit intermolecular potential based on the Kozack-Jordan polarizable model. Putative lowest-energy structures based on this empirical potential are reported. Our results indicate a high propensity for the NH(4)(+) impurity to be fully solvated by water molecules. Clathratelike patterns are formed for clusters containing more than 11 molecules. Local reoptimizations of candidate structures carried out at the DFT level with the B3LYP hybrid functional and the 6-311++G(d,p) basis set confirm the general trends obtained with the intermolecular potential. However, some reorderings between isomers often due to zero-point energy corrections are found in small clusters, leading to stable geometries in agreement with other first-principles studies. Temperature effects have been assessed using a simple harmonic superposition approximation for selected cluster sizes and using dedicated Monte Carlo simulations for (H(2)O)(20)NH(4)(+). The clusters are found to melt near 200 K, and possibly isomerize already below 50 K. The free energy barrier for core/surface isomerization of the impurity in the 21-molecule cluster is estimated to be only a few kcal/mol at 150 K. The vibrational spectroscopic signatures of the clusters obtained from the electronic structure calculations show the usual four O-H stretching bands. As the cluster size increases, the double acceptor-single donor band near 3700 cm(-1) increasingly dominates over the three other bands. While we do not find conclusive evidence for a O-H stretching spectroscopic signature of the ammonium impurity to be in the core or at the surface in the 20-molecule cluster, a possible signature via the N-H stretching bands is suggested near 2800-2900 cm(-1). In the larger (H(2)O)(49)NH(4)(+) cluster, the impurity is slightly more stable at the surface.

Solvation of Na(2) (+) in Ar(n) clusters. I. Structures and spectroscopic properties.

We present a theoretical study of Na(2) (+) solvation in an argon matrix Ar(n) for n=1 to a few tens. We use a model based on an explicit description of valence electron interaction with Na(+) and Ar cores by means of core polarization pseudopotential. The electronic structure determination is thus reduced to a one-electron problem, which can be handled efficiently. We investigate the ground state geometry and related optical absorption of Na(2) (+)Ar(n) clusters. For nA (2)Sigma(u) (+)), which reveals the confinement of the excited A (2)Sigma(u) (+) state. The Na(2) (+) energy spectrum is so strongly perturbed that the A (2)Sigma(u) (+) state becomes higher than the B (2)Pi(u) (+) states. The closure of the first solvation shell is observed at n=17. Above this size, the second solvation shell develops. Its structure is dominated by a pentagonal organization around the Na(2) (+) molecular axis. The optical transitions vary smoothly with n and the A (2)Sigma(u) (+) and B (2)Pi(u) states are no longer inverted, though the first optical transition remains strongly blueshifted.

Experimental demonstration of a pure third-order optical parametric downconversion process.

We report what is, to the best of our knowledge, the first experiment on third-order optical difference frequency generation. The nonlinear medium is a phase-matched KTiOPO4 crystal pumped at 532 and 1665 nm in the picosecond regime. The proposed experimental configuration leads to an efficient and pure cubic process, with generation of 4.5 microJ/pulse at 1474 nm. This result opens the way to new quantum experiments on three-photon Greenberger-Horne-Zeilinger entanglement.