Monte Carlo simulations of the structures and optical absorption spectra of Na atoms in Ar clusters, surfaces, and solids

Trapping and thermal migration of the first- and second-row atoms in Ar, Kr and Xe crystals

Trapping and temperature-induced migration (TIM) of the first- and second-row atoms A from H to Ne in the face-centered cubic rare gas RG = Ar, Kr and Xe crystals are investigated within the classical crystal model parameterized by the empirically modified pairwise potentials. New ab initio coupled cluster A-RG potentials computed in a uniform way for all the atoms A are used to represent the atom-crystal interactions. Absolute and relative stabilities of the substitutional and interstitial trapping sites, their structures, interstitial migration pathways, related activation energies and rough estimates of the TIM rates are obtained. The isotropic model, which neglects non-zero atomic electronic orbital momentum, reveals that migration of interstitial atoms along the network of conjugated fcc octahedral voids is the generic case for atomic mobility. Anisotropic interactions with a crystal inherent to P-state atoms B, C, O and F are accounted for using the non-relativistic diatomics-in-molecule method. Depending on its sign, interaction anisotropy can alter the structures of interstitial trapping sites and transition states remarkably. This, in turn, can dramatically affect the TIM rates. Comparison with reliable experimental data available for oxygen and hydrogen indicates a systematic overestimation of the measured activation energies, by 30% at worst. A comprehensive literature review accomplished for other atoms reveals a lack of information on the TIM processes and rates, though makes it possible to verify a part of the present results on the trapping site energies and structures.

Pairwise Model Potential and DFT Study of Li+Nen Clusters (n = 1-20): The Structural, Electronic, and Thermodynamic Properties

The structural properties, relative stabilities, electronic, and thermodynamic properties, of Li+Nen (n = 1-20) clusters have been studied based on a pairwise model and density functional theory (DFT) methods. In the pairwise method, the potential energy surface considered interactions between Li+Ne, Ne - Ne, and many-body term. For the DFT calculations, the B3LYP functional combined with the 6-311 + + G (2d,2p) basis sets has been employed. In both methods, the Li+Ne6 cluster demonstrated high stability with an octahedral structure, where the Li+ cation was surrounded by Ne atoms. Thus, the octahedral Li+Ne6 structure was considered to be the core for larger cluster sizes. Relative stabilities were assessed based on binding energies, second-order differences of energies, transition dipole moment, and HOMO-LUMO energy gaps. Furthermore, thermodynamic properties were calculated, revealing that the formation process of Li+Nen clusters is endothermic and nonspontaneous.

Modeling of the thermal migration mechanisms of atomic oxygen in Ar, Kr, and Xe crystals

Accommodation and migration of the ground-state (2s²2p^{4 3}P) oxygen atom in the ideal Ar, Kr, and Xe rare gas crystals are investigated using the classical model. The model accounts for anisotropy of interaction between guest and host atoms, spin-orbit coupling, and lattice relaxation. Interstitial and substitutional accommodations are found to be the only thermodynamically stable sites for trapping atomic oxygen. Mixing of electronic states coupled to lattice distortions justifies that its long-range thermal migration follows the adiabatic ground-state potential energy surface. Search for the migration paths reveals a common direct mechanism for interstitial diffusion. Substitutional atoms are activated by the point lattice defects, whereas the direct guest-host exchange meets a higher activation barrier. These three low-energy migration mechanisms provide plausible interpretation for multiple migration activation thresholds observed in Kr and Xe free-standing crystals, confirmed by reasonable agreement between calculated and measured activation energies. An important effect of interaction anisotropy and a minor role of spin-orbit coupling are emphasized.

Stable axially symmetric atomic impurity in an fcc solid-Ba in rare gases

Closed-shell metal atoms in rare gas solids tend to occupy highly symmetric polyhedral crystal sites, as follows from the generic triplet Jahn-Teller splitting of the S → P excitation bands and complies with the isotropic nature of the dispersion forces. Atypical 2 + 1 Jahn-Teller splitting inherent to axially symmetric sites observed recently for Ba atoms has been therefore interpreted as the defect accommodation. By modeling the structure, stability, and spectra of the Ba atom in the face-centered cubic rare gas crystals, we identify thermodynamically stable crystal site of axial C_3v symmetry that explains experimental observations. We also demonstrate the dramatic effect of the interaction anisotropy on the trapping site structure and stability for an excited P-state atom. Our results provide strong evidence for stable axially symmetric accommodation of isotropic impurity in a close-packed lattice.

Luminescence of Atomic Barium in Rare Gas Matrices: A Two-Dimensional Excitation/Emission Spectroscopy Study

A detailed characterization is made of the distinct sites occupied by atomic barium isolated in the three rare gas hosts Ar, Kr, and Xe in excitation scans extracted from the recorded total 6s6p ¹P₁ → (6s)²¹S₀ fluorescence. Extensive use has been made of two-dimensional excitation/emission (2D-EE) spectroscopy to achieve a comprehensive characterization for the wide variety of sites present in the Ba/RG matrix systems. The 2D-EE technique has proved to be a very powerful method to probe the effects of strong intersite reabsorption when extensive spectral overlap occurs between emission and resonance 6s6p ¹P₁ ← (6s)²¹S₀ absorption of barium atoms occupying multiple sites. Two-dimensional excitation/emission scans have also been used in this study to monitor the effects of sample annealing and thereby identify the thermally stable sites of isolation. Sites of the same type occupied by atomic barium in the three host solids are identified in resolved excitation spectra and are associated on the basis of the observed matrix shift versus host polarizability. Following site associations, the photophysical properties of each matrix site were characterized revealing that the Stokes shift was greatest in the blue site, smallest for the violet site, and intermediate for the green site. The emission temperature dependences and excited state lifetimes were recorded, indicating that measured radiative lifetimes of 4-5 ns were in good agreement with the gas phase value of 8.4 ns when corrected for the effective field of the solids. The only exception to this was the blue site in Ba/Xe, where a nonradiative quenching channel exists even at 9.8 K that competes effectively with the nanosecond fluorescence. An unusual, asymmetric 2 + 1 excitation band has been recorded for atomic barium in the three rare gas hosts in addition to the threefold split, Jahn-Teller bands typically observed for P ← S absorptions of matrix-isolated metal atoms. Possible assignments of the sites responsible for these band shapes are made on the basis of recent spectral simulations obtained from molecular dynamics calculations on the Ba/Xe system.

An investigation of the sites occupied by atomic barium in solid xenon-A 2D-EE luminescence spectroscopy and molecular dynamics study

A detailed characterisation of the luminescence recorded for the 6p ¹P₁-6s ¹S₀ transition of atomic barium isolated in annealed solid xenon has been undertaken using two-dimensional excitation-emission (2D-EE) spectroscopy. In the excitation spectra extracted from the 2D-EE scans, two dominant thermally stable sites were identified, consisting of a classic, three-fold split Jahn-Teller band, labeled the blue site, and an unusual asymmetric 2 + 1 split band, the violet site. A much weaker band has also been identified, whose emission is strongly overlapped by the violet site. The temperature dependence of the luminescence for these sites was monitored revealing that the blue site has a non-radiative channel competing effectively with the fluorescence even at 9.8 K. By contrast, the fluorescence decay time of the violet site was recorded to be 4.3 ns and independent of temperature up to 24 K. The nature of the dominant thermally stable trapping sites was investigated theoretically with Diatomics-in-Molecule (DIM) molecular dynamics simulations. The DIM model was parameterized with ab initio multi-reference configuration interaction calculations for the lowest energy excited states of the Ba⋅Xe pair. The simulated absorption spectra are compared with the experimental results obtained from site-resolved excitation spectroscopy. The simulations allow us to assign the experimental blue feature spectrum to a tetra-vacancy trapping site in the bulk xenon fcc crystal-a site often observed when trapping other metal atoms in rare gas matrices. By contrast, the violet site is assigned to a specific 5-atom vacancy trapping site located at a grain boundary.

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.

Absorption Spectroscopy, a Tool for Probing Local Structures and the Onset of Large-Amplitude Motions in Small KAr(n) Clusters at Increasing Temperatures

Photoabsorption spectra of KArn (n = 1-10) are simulated at temperatures ranging between 5 and 25 K. The calculations associate a Monte Carlo (MC) method to sample cluster geometries at temperature T, with a one-electron ab initio model to calculate the ground-state and excited-state energies of the cluster. The latter model replaces the K(+) core electrons and all the electrons of the Ar atoms by appropriate pseudopotentials, complemented by core polarization potentials. It also provides the necessary oscillator strengths to simulate the spectra. Global optimization by basin-hopping is used in combination with MC simulation at low temperature (5 K) to identify the most stable isomer and remarkable isomers of ground-state KArn clusters, which are stable with respect to deformations of the order of those expected with Zero Point Energy motions. The absorption spectra calculated for each of these isomers at 5 K suggest that absorption spectroscopy can probe sensitively the local environment of K atom: surface location of K with respect to a close-packed Ar moiety, number of Ar atom in close vicinity, and local symmetry about K. Simulation at increasing temperatures, up to the evaporation limit of K out of the cluster, shows the onset of large amplitude motions above 20 K, when the K atom experiences a variety of local environments.

Infrared-active spin-orbit transitions of halogen atom dopants in solid parahydrogen: the role of trapping site geometry

We present theoretical calculations of the (2)P(1/2) ← (2)P(3/2) spin-orbit transition of Cl dopants embedded as substitutional impurities in solid parahydrogen (pH2) matrices. In the lower-energy (2)P(3/2) spin-orbit level, the Cl atom's electron density distribution is anisotropic, and slightly distorts the geometry of the atom's trapping site. This distortion leads to a blue shift in the spin-orbit transition energy; the blue shift is enhanced when we account for the large-amplitude zero point motions of the pH2 molecules surrounding the Cl dopant. We also show that the intensity of the transition depends on the geometry of the trapping site. In the gas phase, the (2)P(1/2) ← (2)P(3/2) atomic transition is electric dipole forbidden. However, when the Cl atom resides in trapping sites that mimic the hexagonal close packed morphology of pure solid pH2, the transition becomes electric dipole allowed through interaction-induced transition dipole moments. These transition dipole moments originate in the anisotropic electron density distribution of the lower-energy (2)P(3/2) spin-orbit level.

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)Π(½)|Ĥ(SO)|(2)Π(½)> and b(R)=<(2)Π(-½)|Ĥ(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.

Rydberg states of small NaArn* clusters

The 4s and 5s Rydberg excited states of NaAr(n)* clusters are investigated using a pseudopotential quantum-classical method. While NaAr(n) clusters in their ground state are known to be weakly bound van der Waals complexes with Na lying at the surface of the argon cluster, isomers in 4s or 5s electronically excited states of small NaAr(n)* clusters (n< or =10) are found to be stable versus dissociation. The relationship between electronic excitation and cluster geometry is analyzed as a function of cluster size. For both 4s and 5s states, the stable exciplex isomers essentially appear as sodium-centered structures with similar topologies, converging towards those of the related NaAr(n)+ positive ions when the excitation level is increased. This is consistent with a Rydberg-type picture for the electronically excited cluster, described by a central sodium ion solvated by an argon shell, and an outer diffuse electron orbiting around this NaAr(n)+ cluster core.

The OH radical-H2O molecular interaction potential

The OH radical is one of the most important oxidants in the atmosphere due to its high reactivity. The study of hydrogen-bonded complexes of OH with the water molecules is a topic of significant current interest. In this work, we present the development of a new analytical functional form for the interaction potential between the rigid OH radical and H(2)O molecules. To do this we fit a selected functional form to a set of high level ab initio data. Since there is a low-lying excited state for the H(2)O.OH complex, the impact of the excited state on the chemical behavior of the OH radical can be very important. We perform a potential energy surface scan using the CCSD(T)/aug-cc-pVTZ level of electronic structure theory for both excited and ground states. To model the physics of the unpaired electron in the OH radical, we develop a tensor polarizability generalization of the Thole-type all-atom polarizable rigid potential for the OH radical, which effectively describes the interaction of OH with H(2)O for both ground and excited states. The stationary points of (H(2)O)(n)OH clusters were identified as a benchmark of the potential.

Path-integral Monte Carlo simulation of the recombination of two Al atoms embedded in parahydrogen

We report the use of path-integral Monte Carlo (PIMC) simulations in the study of the stability against recombination of two Al atoms trapped in solid parahydrogen (pH2) at 4 K. The many-body interactions involving open-shell Al atoms are described with a pairwise additive Hamiltonian model. To estimate the lifetime against recombination, we use PIMC simulations to define an effective potential averaged over the position of the pH2 molecules, followed by a transition-state treatment. Different initial embedding sites are explored. If the initial substitution sites are within a distance of approximately 13 bohrs, the Al atoms will significantly distort the lattice structure to allow recombination, with an accompanying release of energy during the process. For substitution distances longer than approximately 14 bohrs, the dispersion of Al atoms is shown to be metastable, with lifetimes varying from approximately 30 min to several days. The electronic anisotropy is a factor that helps to stabilize the dispersion.

Moment Analysis Method As Applied to the 2S → 2P Transition in Cryogenic Alkali Metal/Rare Gas Matrices

The moment analysis method (MA) has been tested for the case of 2S --> 2P ([core]ns1 --> [core]np1) transitions of alkali metal atoms (M) doped into cryogenic rare gas (Rg) matrices using theoretically validated simulations. Theoretical/computational M/Rg system models are constructed with precisely defined parameters that closely mimic known M/Rg systems. Monte Carlo (MC) techniques are then employed to generate simulated absorption and magnetic circular dichroism (MCD) spectra of the 2S --> 2P M/Rg transition to which the MA method can be applied with the goal of seeing how effective the MA method is in re-extracting the M/Rg system parameters from these known simulated systems. The MA method is summarized in general, and an assessment is made of the use of the MA method in the rigid shift approximation typically used to evaluate M/Rg systems. The MC-MCD simulation technique is summarized, and validating evidence is presented. The simulation results and the assumptions used in applying MA to M/Rg systems are evaluated. The simulation results on Na/Ar demonstrate that the MA method does successfully re-extract the 2P spin-orbit coupling constant and Landé g-factor values initially used to build the simulations. However, assigning physical significance to the cubic and noncubic Jahn-Teller (JT) vibrational mode parameters in cryogenic M/Rg systems is not supported.

Simple DFT model of clusters embedded in rare gas matrix: Trapping sites and spectroscopic properties of Na embedded in Ar

We present a theoretical model to study the dynamics of metallic clusters embedded in a rare gas matrix. We describe the active electrons of the embedded cluster using time dependent density functional theory, while the surrounding matrix is described in terms of classical molecular dynamics of polarizable atoms. The coupling between the cluster and the rare gas atoms is deduced from the work of Gross and Spiegelmann [J. Chem. Phys. 108, 4148 (1998)] and reformulated explicitly in a simple and efficient density functional form. The electron rare gas interaction takes the form of an averaged dipole fluctuation term, which retains the van der Waals long range interaction, and a short range repulsive pseudopotential, which accounts for the Pauli repulsion of the electron by the rare gas atom. We applied our model to Na clusters embedded in Ar matrix. For the latter we developed an efficient local pseudopotential, which allows studying systems containing more than 10(3) Ar atoms. We show that large systems are indeed necessary to account properly for long range polarization of the matrix, that competes with the matrix confinement effect. We focus our study on Na(2), Na(4), and Na(8). For each system, we have determined the geometry of the most favorable trapping site by means of damped molecular dynamics. We present the effect of matrix embedding on the optical absorption spectrum. For Na(2), the trapping site can be unambiguously identified by comparison of the absorption spectrum with experiment. For Na(4) the spectrum of the embedded cluster is significantly different from the free cluster spectrum, while for Na(8) differences are less pronounced.

Modeling of the Three-Body Effects in the Neutral Trimers in the Quartet State by ab initio Calculations. H3, Na3, and Na2B

The Na2B, Na3, and H3 trimers in the lowest quartet states were studied by ab initio methods, using both the supermolecular approach and the intermolecular Møller-Plesset perturbation theory. Partitioning of the nonadditive contribution into the orientational two-body part and the genuine three-body part was proposed. The lowest quartet state of the Na3 trimer and all the three lowest quartet states of the Na2B trimer are bound, and the forms of these clusters are essentially determined by two-body forces. In the case of the Na2B trimer the orientational two-body nonadditivity proved to be crucial. In addition, in the title metal trimers, in the region of the van der Waals minima, the genuine nonadditivity is very important, and amounts to 30% in Na2B and up to 70% in Na3. The leading nonadditive term is the triple-exchange Heitler-London exchange term. For triangular arrangements it considerably enhances the total stabilization. The single-exchange term and the SCF deformation play only a secondary role. The dispersion nonadditivity is negligible. The isotropic part of the basis set superposition error (BSSE) is large and must be corrected by the counterpoise method. The anisotropic contribution to BSSE is practically negligible.