Accurateab initiopotential for argon dimer including highly repulsive region

An algorithm for very high pressure molecular dynamics simulations

We describe a method to run simulations of ground or excited state dynamics under extremely high pressures. The method is based on the introduction of a fictitious ideal gas that exerts the required pressure on a molecular sample and is therefore called XP-GAS (eXtreme Pressure by Gas Atoms in a Sphere). The algorithm is most suitable for approximately spherical clusters of molecules described by quantum chemistry methods, Molecular Mechanics or mixed QM/MM approaches. We compare the results obtained by the algorithm here presented and by the XP-PCM approach, based on a continuum description of the environment. As a test case, we study the conformational dynamics of 1,3-butadiene either as an isolated molecule ("naked" butadiene) or embedded in a cluster of argon atoms, under pressures up to 15 GPa. Overall, our results show that the XP-GAS QM/MM simulation method is in good agreement with the XP-PCM QM/Continuum model (Cammi model) in describing the effect of the pressure on static properties as the equilibrium geometry of butadiene in the ground state. Furthermore, the comparison of XP-GAS simulations with naked butadiene and butadiene in argon shows the importance, for XP-GAS and related methods, of a realistic representation of the medium in modelling pressure effects.

A nuclear configuration interaction approach to study nuclear spin effects: an application to ortho- and para-3 He2 @C60

We introduce a non-orthogonal configuration interaction approach to investigate nuclear quantum effects on energies and densities of confined fermionic nuclei. The Hamiltonian employed draws parallels between confined systems and many-electron atoms, where effective non-Coulombic potentials represent the interactions of the trapped particles. One advantage of this method is its generality, as it offers the potential to study the nuclear quantum effects of various confined species affected by effective isotropic or anisotropic potentials. As a first application, we analyze the quantum states of two 3 He atoms encapsulated in C60 . At the Hartree-Fock level, we observe the breaking of spin and spatial symmetries. To ensure wavefunctions with the correct symmetries, we mix the broken-symmetry Hartree-Fock states within the non-orthogonal configuration interaction expansion. Our proposed approach predicts singly and triply degenerate ground states for the singlet (para-3 He2 @C60 ) and triplet (ortho-3 He2 @C60 ) nuclear spin configurations, respectively. The ortho-3 He2 @C60 ground state is 5.69 cm-1 higher in energy than the para-3 He2 @C60 ground state. The nuclear densities obtained for these states exhibit the icosahedral symmetry of the C60 embedding potential. Importantly, our calculated energies for the lowest 85 states are in close agreement with perturbation theory results based on a harmonic oscillator plus rigid rotor model of 3 He2 @C60 .

Generalized perturbative singles corrections to the random phase approximation method: Impact on noncovalent interaction energies of closed- and open-shell dimers

The post-Kohn-Sham (KS) random phase approximation (RPA) method may provide a poor description of interaction energies of weakly bonded molecules due to inherent density errors in approximate KS functionals. To overcome these errors, we develop a generalized formalism to incorporate perturbative singles (pS) corrections to the RPA method using orbital rotations as a perturbation parameter. The pS schemes differ in the choice of orbital-rotation gradient and Hessian. We propose a pS scheme termed RPA singles (RPAS)[Hartree-Fock (HF)] that uses the RPA orbital-rotation gradient and time-dependent HF Hessian. This correction reduces the errors in noncovalent interaction energies of closed- and open-shell dimers. For the open-shell dimers, the RPAS(HF) method leads to a consistent error reduction by 50% or more compared to the RPA method for the cases of hydrogen-bonding, metal-solvent, carbene-solvent, and dispersion interactions. We also find that the pS corrections are more important in error reduction compared to higher-order exchange corrections to the RPA method. Overall, for open shells, the RPAS(HF)-corrected RPA method provides chemical accuracy for noncovalent interactions and is more reliable than other perturbative schemes and dispersion-corrected density functional approximations, highlighting its importance as a reliable beyond-RPA correction.

Molecules in Environments: Toward Systematic Quantum Embedding of Electrons and Drude Oscillators

We develop a quantum embedding method that enables accurate and efficient treatment of interactions between molecules and an environment, while explicitly including many-body correlations. The molecule is composed of classical nuclei and quantum electrons, whereas the environment is modeled via charged quantum harmonic oscillators. We construct a general Hamiltonian and introduce a variational Ansatz for the correlated ground state of the fully interacting molecule-environment system. This wave function is optimized via the variational Monte Carlo method and the ground state energy is subsequently estimated through the diffusion Monte Carlo method. The proposed scheme allows an explicit many-body treatment of electrostatic, polarization, and dispersion interactions between the molecule and the environment. We study solvation energies and excitation energies of benzene derivatives, obtaining excellent agreement with explicit ab initio calculations and experiments.

The ab initio potential energy curves of atom pairs and transport properties of high-temperature vapors of Cu and Si and their mixtures with He, Ar, and Xe gases

The potential energy curves (PECs) for the homonuclear He-He, Ar-Ar, Cu-Cu, and Si-Si dimers, as well as heteronuclear Cu-He, Cu-Ar, Cu-Xe, Si-He, Si-Ar, and Si-Xe dimers, are obtained in quantum Monte Carlo (QMC) calculations. It is shown that the QMC method provides the PECs with an accuracy comparable with that of the state-of-the-art coupled cluster singles and doubles with perturbative triples corrections [CCSD(T)] calculations. The QMC data are approximated by the Morse long range (MLR) and (12-6) Lennard-Jones (LJ) potentials. The MLR and LJ potentials are used to calculate the deflection angles in binary collisions of corresponding atom pairs and transport coefficients of Cu and Si vapors and their mixtures with He, Ar, and Xe gases in the range of temperature from 100 K to 10 000 K. It is shown that the use of the LJ potentials introduces significant errors in the transport coefficients of high-temperature vapors and gas mixtures. The mixtures with heavy noble gases demonstrate anomalous behavior when the viscosity and thermal conductivity can be larger than that of the corresponding pure substances. In the mixtures with helium, the thermal diffusion factor is found to be unusually large. The calculated viscosity and diffusivity are used to determine parameters of the variable hard sphere and variable soft sphere molecular models as well as parameters of the power-law approximations for the transport coefficients. The results obtained in the present work include all information required for kinetic or continuum simulations of dilute Cu and Si vapors and their mixtures with He, Ar, and Xe gases.

Rotational spectroscopy of the argon dimer by time-resolved Coulomb explosion imaging of rotational wave packets

We report time-domain rotational spectroscopy of the argon dimer, Ar₂, by implementing time-resolved Coulomb explosion imaging of rotational wave packets. The rotational wave packets are created in Ar₂ with a linearly polarized, nonresonant, ultrashort laser pulse, and their spatiotemporal evolution is fully characterized by measuring angular distribution of the fragmented Ar⁺ promptly ejected from Ar₂²⁺ generated by the more intense probe pulse. The pump-probe measurements have been carried out up to a delay time of 16 ns. The alignment parameters, derived from the observed images, exhibit periodic oscillation lasting for more than 15 ns. The pure rotational spectrum of Ar₂ is obtained by Fourier transformation of the time traces of the alignment parameters. The frequency resolution in the spectrum is about 90 MHz, the highest ever achieved for Ar₂. The rotational constant and the centrifugal distortion constant are determined with much improved precision than the previous experimental results: B₀ = 1.72713 ± 0.00009 GHz and D₀ = 0.0310 ± 0.0005 MHz. The present B₀ value does not match within the quoted experimental uncertainty with that from the VUV spectroscopy, so far accepted as an experimental reference to assess theories. The present improved constants would stand as new references to calibrate state-of-the-art theoretical investigations and an indispensable experimental source for the construction of an accurate empirical intermolecular potential.

Theoretical modeling of molecules in weakly interacting environments: trifluoride anions in argon

The properties of molecules can be affected by the presence of a host environment. Even in inert rare gas matrices such effects are observable, as for instance in matrix isolation spectroscopy. In this work we study the trifluoride anion in cryogenic argon environments. To investigate the structure and vibrational properties of the guest-host systems, a potential energy surface of compound F-3-argon structures is determined from ab initio calculations with the CCSD(T)-F12b approach. Argon environments are probed with minima hopping optimizations of extended trifluoride-argon clusters. The vibrations of F-3 within the optimized environments are examined with anharmonic vibrational analyses. Among the three identified structural surroundings for the trifluoride, two are characterized by relatively favorable guest-host and host-host interactions as well as vibrational zero-point energies. A striking dependence of the trifluoride properties on the particular argon environment reveals the delicate influence of the host atoms on the guest molecule. Very good agreement with measured data suggests that in experiment F-3 occupies a double-vacancy site.

O(N) Stochastic Evaluation of Many-Body van der Waals Energies in Large Complex Systems

We propose a new strategy to solve the key equations of the many-body dispersion (MBD) model by Tkatchenko, DiStasio Jr., and Ambrosetti. Our approach overcomes the original O(N3) computational complexity that limits its applicability to large molecular systems within the context of O(N) density functional theory. First, to generate the required frequency-dependent screened polarizabilities, we introduce an efficient solution to the Dyson-like self-consistent screening equations. The scheme reduces the number of variables and, coupled to a direct inversion of the iterative subspace extrapolation, exhibits linear-scaling performances. Second, we apply a stochastic Lanczos trace estimator resolution to the equations evaluating the many-body interaction energy of coupled quantum harmonic oscillators. While scaling linearly, it also enables communication-free pleasingly parallel implementations. As the resulting O(N) stochastic massively parallel MBD approach is found to exhibit minimal memory requirements, it opens up the possibility of computing accurate many-body van der Waals interactions of millions-atoms' complex materials and solvated biosystems with computational times in the range of minutes.

Solvation of Isoelectronic Halide and Alkali Metal Ions by Argon Atoms

This work systematically examines the interactions of alkali metal cations and their isoelectronic halide counterparts with up to six solvating Ar atoms (M⁺Ar_n and X^-Ar_n, where M = Li, Na, K, and Rb; X = H, F, Cl, and Br; and n = 1-6) via full geometry optimizations with the MP2 method and robust, correlation-consistent quadruple-ζ (QZ) basis sets. 116 unique M⁺Ar_n and X^-Ar_n stationary points have been characterized on the MP2/QZ potential energy surface. To the best of our knowledge, approximately two dozen of these stationary points have been reported here for the first time. Some of these new structures are either the lowest-energy stationary point for a particular cluster or energetically competitive with it. The CCSD(T) method was employed to perform additional single-point energy computations upon all MP2/QZ-optimized structures using the same basis set. CCSD(T)/QZ results indicate that internally solvated structures with the ion at/near the geometric center of the cluster have appreciably higher energies than those placing the ion on the periphery. While this study extends the prior investigations of M⁺Ar_n clusters found within the literature, it notably provides one of the first thorough characterizations of and comparisons to the corresponding negatively charged X^-Ar_n clusters.

Thermophysical properties of low-density neon gas from highly accurate first-principles calculations and dielectric-constant gas thermometry measurements

New interatomic potential energy and interaction-induced polarizability curves for two ground-state neon atoms were developed and used to predict the second density, acoustic, and dielectric virial coefficients and the dilute gas shear viscosity and thermal conductivity of neon at temperatures up to 5000 K. The potential energy curve is based on supermolecular coupled-cluster (CC) calculations at very high levels up to CC with single, double, triple, quadruple, and perturbative pentuple excitations [CCSDTQ(P)]. Scalar and spin-orbit relativistic effects, the diagonal Born-Oppenheimer correction, and retardation of the dispersion interactions were taken into account. The interaction-induced polarizability curve, which in this work is only needed for the calculation of the second dielectric virial coefficient, is based on supermolecular calculations at levels up to CCSDT and includes a correction for scalar relativistic effects. In addition to these first-principles calculations, highly accurate dielectric-constant gas thermometry (DCGT) datasets measured at temperatures from 24.5 to 200 K were analyzed to obtain the difference between the second density and dielectric virial coefficients with previously unattained accuracy. The agreement of the DCGT values with the ones resulting from the first-principles calculations is, despite some small systematic deviations, very satisfactory. Apart from this combination of two virial coefficients, the calculated thermophysical property values of this work are significantly more accurate than any available experimental data.

The development of a full range analytical interatomic potential

A chronological account is given to the development of a full range interatomic potential. Starting with a simple phenomenological model, the terms in the model are gradually modified, so that they can carry some definite physical meaning. To gain insight, a systematic, order by order interaction potential theory is developed. Conversely, this theory suggests the functional form for the potential model. At present, we have a simple interaction model that is capable of describing the van der Waals potentials of many systems from R = 0 to R→∞.

Effect of the range of particle cohesion on the phase behavior and thermodynamic properties of fluids

Molecular simulations are performed for the (m + 1, m) potential to systematically investigate the effect of changing the range of particle cohesion on both vapor-liquid equilibria and thermodynamic properties of fluids. The results are reported for m = 4-11, which represent a progressive narrowing of the potential energy well. The conventional Lennard-Jones potential is used as a reference point for normal fluid behavior. Small values of m result in a broadening of the phase envelope compared with the Lennard-Jones potential, whereas a contraction is observed in other cases. The critical properties are reported, and a relationship between the critical temperature and the Boyle temperature is determined. The low values of the critical compressibility factor when m < 6 reflect the behavior observed for real fluids such as n-alkanes. The results for supercritical thermodynamic properties are much more varied. Properties such as pressure, potential energy, isochoric thermal pressure coefficient, and thermal expansion coefficient vary consistently with m, whereas other properties such as the Joule-Thomson coefficient exhibit much more nuanced behavior. Maximum and minimum values are reported for both the isochoric heat capacity and isothermal compressibility. A minimum in the speed of sound is also observed.

Conformal Analytical Potential for All the Rare Gas Dimers over the Full Range of Internuclear Distances

An analytical model for the potential between two rare gas atoms at distances between R=0 to R→∞ is assumed to be conformal with the previously published potential for He_{2} [J. Chem. Phys. 142, 131102 (2015)JCPSA60021-960610.1063/1.4916740]. The potential curves of the rare gas dimers all have the same shape and only depend on the well parameters D_{e} and R_{e}. The potentials and the vibrational levels for the 11 homonuclear and heteronuclear dimers for which recent ab initio calculations are available agree, within several percent, with the ab initio results. For the other rare gas dimers, the new potential provides the first realistic estimates for the potentials.

Dispersion Energy from Local Polarizability Density

A simple nonlocal functional for calculation of dispersion energies is proposed. Compared to a similar formula used earlier, we introduced a regularization to remove its singularities and used a dynamic polarizability density similar to those in the so-called van der Waals density functionals. The performance of the new functional is tested on dispersion energies for a set of representative dimers, and it is found that it is significantly more accurate than published nonlocal functionals.

Do Semilocal Density-Functional Approximations Recover Dispersion Energies at Small Intermonomer Separations?

The methods that add dispersion energies to interaction energies computed using density-functional theory (DFT), known as DFT+D methods, taper off the dispersion energies at distances near van der Waals minima and smaller based on an assumption that DFT starts to reproduce the dispersion energies there. We show that this assumption is not correct as the alleged contribution behaves unphysically and originates to a large extent from nonexchange-correlation terms. Thus, dispersion functions correct DFT in this region for deficiencies unrelated to dispersion interactions.

Determination of Interatomic Potentials of He_{2}, Ne_{2}, Ar_{2}, and H_{2} by Wave Function Imaging

We report on a direct method to measure the interatomic potential energy curve of diatomic systems. A cold target recoil ion momentum spectroscopy reaction microscope was used to measure the squares of the vibrational wave functions of H_{2}, He_{2}, Ne_{2}, and Ar_{2}. The Schrödinger equation relates the curvature of the wave function to the potential V(R) and therefore offers a simple but elegant way to extract the shape of the potential.

Low-temperature chemistry using the R-matrix method

Techniques for producing cold and ultracold molecules are enabling the study of chemical reactions and scattering at the quantum scattering limit, with only a few partial waves contributing to the incident channel, leading to the observation and even full control of state-to-state collisions in this regime. A new R-matrix formalism is presented for tackling problems involving low- and ultra-low energy collisions. This general formalism is particularly appropriate for slow collisions occurring on potential energy surfaces with deep wells. The many resonance states make such systems hard to treat theoretically but offer the best prospects for novel physics: resonances are already being widely used to control diatomic systems and should provide the route to steering ultracold reactions. Our R-matrix-based formalism builds on the progress made in variational calculations of molecular spectra by using these methods to provide wavefunctions for the whole system at short internuclear distances, (a regime known as the inner region). These wavefunctions are used to construct collision energy-dependent R-matrices which can then be propagated to give cross sections at each collision energy. The method is formulated for ultracold collision systems with differing numbers of atoms.

Predicting vapor liquid equilibria using density functional theory: A case study of argon

Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase.

Simple Modifications of the SCAN Meta-Generalized Gradient Approximation Functional

We analyzed various possibilities to improve upon the SCAN meta-generalized gradient approximation density functional obeying all known properties of the exact functional that can be satisfied at this level of approximation. We examined the necessity of locally satisfying a strongly tightened lower bound for the exchange energy density in single-orbital regions, the nature of the error cancellation between the exchange and correlation parts in two-electron regions, and the effect of the fourth-order term in the gradient expansion of the correlation energy density. We have concluded that the functional can be modified to separately reproduce the exchange and correlation energies of the helium atom by locally releasing the strongly tightened lower bound for the exchange energy density in single-orbital regions, but this leads to an unbalanced improvement in the single-orbital electron densities. Therefore, we decided to keep the F_X ≤ 1.174 exact condition for any single-orbital density, where F_X is the exchange enhancement factor. However, we observed a general improvement in the single-orbital electron densities by revising the correlation functional form to follow the second-order gradient expansion in a wider range. Our new revSCAN functional provides more-accurate atomization energies for the systems with multireference character, compared to the SCAN functional. The nonlocal VV10 dispersion-corrected revSCAN functional yields more-accurate noncovalent interaction energies than the VV10-corrected SCAN functional. Furthermore, its global hybrid version with 25% of exact exchange, called revSCAN0, generally performs better than the similar SCAN0 for reaction barrier heights. Here, we also analyzed the possibility of the construction of a local hybrid from the SCAN exchange and a specific locally bounded nonconventional exact exchange energy density. We predict compatibility problems since this nonconventional exact exchange energy density does not really obey the strongly tightened lower bound for the exchange energy density in single-orbital regions.

Ab initio calculations of many-body interactions for compressed solid argon

An investigation on many-body effects of solid argon at high pressure was conducted based on a many-body expansion of interaction energy. The three- and four-body terms in the expansion were calculated using the coupled-cluster method with single, double, and noniterative triple theory and incremental method, in which the configurations of argon trimers and tetramers were chosen as the same as those in the actual lattice. The four-body interactions in compressed solid argon were estimated for the first time, and the three-body interaction ab initio calculations were extended to a small distance. It shows that the four-body contribution is repulsive at high densities and effectively cancels the three-body lattice energy. The dimer potential plus three-body interaction can well reproduce the measurements of equation of state at pressure approximately lower than ∼60 GPa, when including the four-body effects extends the agreement up to the maximum experimental pressure of 114 GPa.

Strongly Constrained and Appropriately Normed Semilocal Density Functional

The ground-state energy, electron density, and related properties of ordinary matter can be computed efficiently when the exchange-correlation energy as a functional of the density is approximated semilocally. We propose the first meta-generalized-gradient approximation (meta-GGA) that is fully constrained, obeying all 17 known exact constraints that a meta-GGA can. It is also exact or nearly exact for a set of "appropriate norms," including rare-gas atoms and nonbonded interactions. This strongly constrained and appropriately normed meta-GGA achieves remarkable accuracy for systems where the exact exchange-correlation hole is localized near its electron, and especially for lattice constants and weak interactions.

Ab initio study of the structures and electronic states of small neutral and ionic DABCO--Ar(n) clusters

In the present theoretical work, we investigated the stationary points (minima and transition states) on the ground state potential energy surfaces of neutral and ionic 1,4-diazabicyclo[2.2.2]octane (DABCO)--Ar(n)⁰,⁺¹ (n = 1-4) clusters. As established in our systematic work on DABCO--Ar cluster (Mathivon et al., J Chem Phys 139:164306, 2013), the (R)MP2/aug-cc-pVDZ level is accurate enough for validating the prediction of stable forms. For n = 1 and 2, further computations at the MP2/aug-cc-pVTZ level confirm these assumptions. We show that some of the already known isomers of these heteroclusters derived using lower levels of theory are not realistic. More interestingly, our work reveals that DABCO is subject to slight deformations when binding to a small number of Ar atoms. Moreover, we computed the potential energy surfaces of the lowest singlet electronic states of DABCO--Ar(n)(n = 1-3) and of DABCO⁺--Ar(n)(n = 1-3), and the transition moments for the Sp(p = 1-3) ← S0 neutral transitions. These electronic states are found to be Rydberg in nature. The shape of their potentials is mainly repulsive with slight stabilization in the S2 potentials. Finally, the effects of microsolvation of DABCO in Ar clusters in ground and electronic excited states are discussed. The photophysical and photochemical dynamics of these electronic states may be complex.

Steering the nuclear motion in singly ionized argon dimers with mutually detuned laser pulses

We demonstrate that the vibrational nuclear motion of singly ionized argon dimers can be controlled with two ultrashort laser pulses of different wavelengths. In particular, we observe a striking "gap" in the pump-probe-delay-dependent kinetic-energy release spectrum only if the probe-pulse wavelength exceeds the pump-pulse wavelength. This "frustrated dissociation effect" is reproduced by our two-state quantum mechanical model, validating its interpretation as a pump-pulse-initiated population transfer between dipole-coupled Born-Oppenheimer electronic states of the dissociating Ar(2)(+) molecular ion. Our numerical results also reproduce the measured collapse and fractional revival of the oscillating Ar(2)(+) nuclear wave packet, and, for single-pulse dissociation, the decrease of the kinetic-energy release with increasing laser wavelength.

Dispersion energy evaluated by using locally projected occupied and excited molecular orbitals for molecular interaction

The dispersion terms are evaluated with the perturbation theory based on the locally projected molecular orbitals. A series of model systems, including some of the S22 set, is examined, and the calculated binding energies are compared with the published results. The basis set dependence is also examined. The dispersion energy correction is evaluated by taking into account the double excitations only of the dispersion type electron configurations and is added to the 3rd order single excitation perturbation energy, which is a good approximation to the counterpoise (CP) corrected Hartree-Fock (HF) binding energy. The procedure is the approximate "CP corrected HF + D" method. It ensures that the evaluated binding energy is approximately free of the basis set superposition error without the CP procedure. If the augmented basis functions are used, the evaluated binding energies for the predominantly dispersion-bound systems, such as rare gas dimers and halogen bonded clusters, agree with those of the reference calculations within 1 kcal mol(-1) (4 kJ mol(-1)). The limitation of the present method is also discussed.

Multiple recapture of electrons in multiple ionization of the argon dimer by a strong laser field

We observe multiply frustrated tunneling ionization-induced dissociation of the argon dimers by intense linearly polarized ultrashort laser pulses. By measuring the kinetic energy release and angular distribution of the Coulomb explosion of up to eightfold ionized argon dimers, we can trace the recapture of up to two electrons to Rydberg states of the highly charged compound at the end of the laser pulse. Upon dissociation of the dimer, the Rydberg electron prefers to localize at the atomic ion with the higher charge state. We probe the electron recapture dynamics by a time-delayed weak pulse.

Dynamics of relaxation and fragmentation in size-selected icosahedral Ar(n)[NO(-)(v = 1)] clusters

We study the vibrational relaxation and solvation dynamics in size-selected icosahedral Ar(n)(NO(-)) at 300 K, where NO(-)(X(3)Σ(-)) is in v = 1 and n = 1-12, using a classical dynamics method and an interaction model consisting of detailed host-guest and host-host interactions. Two relaxation time scales are found: (i) the short-time (<200 ps), in which rate is nearly independent of cluster size, and (ii) the ns scale, in which a slow energy transfer process occurs between NO(-) vibration and argon modes at a rate (~10(8) s(-1)) decreasing slightly from n = 12 to 6 and rapidly from n = 5 to 1 (~10(6) s(-1)). In Ar(12)(NO(-)), less than one-quarter of the host atoms sampled evaporate, nearly 60% of evaporation occurring within 200 ps caused by rapid energy transfer from NO(-) at short time. The fraction of evaporation decreases nearly exponentially with increasing evaporation time, but ~16% of evaporation still occurs on a time scale longer than 1 ns. Evaporation from one hemisphere of Ar(12)(NO(-)) dominates the rest. Final cluster sizes commonly produced from the fragmentation of Ar(12)(NO(-)) are n = 6-11 (evaporation of 6-1 atoms) and n = 12 (no evaporation).