Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermophysical properties of FLiBe using moment tensor potentials

Fluoride salts are prospective materials for applications in some next-generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (66.6% LiF – 33.3% BeF2) predicted by the ab initio (DFT-D3) method. Here, we then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi/Fumi.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A systematic approach to generating accurate neural network potentials: the case of carbon

Abstract Availability of affordable and widely applicable interatomic potentials is the key needed to unlock the riches of modern materials modeling. Artificial neural network-based approaches for generating potentials are promising; however, neural network training requires large amounts of data, sampled adequately from an often unknown potential energy surface. Here we propose a self-consistent approach that is based on crystal structure prediction formalism and is guided by unsupervised data analysis, to construct an accurate, inexpensive, and transferable artificial neural network potential. Using this approach, we construct an interatomic potential for carbon and demonstrate its ability to reproduce first principles results on elastic and vibrational properties for diamond, graphite, and graphene, as well as energy ordering and structural properties of a wide range of crystalline and amorphous phases.

Chemistry↗
DOE OSTI2020

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF–Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born–Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. In this study, using tetrahydrofuran (THF)–water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D 8 THF to examine the molecular origins of the isotope-dependent loss of the THF–water miscibility gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S–S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S–S and C–S bond cleavage, which involve the lowest excited electronic states S 1 and S 2 . Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S 1 –S 6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S–S and C–S bond dissociation coordinates. Finally, we present calculations of spin–orbit coupling constants at the Franck–Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S 0 –S 6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck–Condon region and multiconfigurational character toward the bond-dissociation limits. In conclusion, we compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Rishi, Varun↗
DOE OSTI2026

Unsupervised atomic data mining via multi-kernel graph autoencoders for machine learning force fields

Constructing a chemically diverse dataset while avoiding sampling bias is critical to training efficient and generalizable force fields. However, in computational chemistry and materials science, many common dataset generation techniques are prone to oversampling regions of the potential energy surface. Furthermore, these regions can be difficult to identify and isolate from each other or may not align well with human intuition, making it challenging to systematically remove bias in the dataset. While traditional clustering and pruning (down-sampling) approaches can be useful for this, they can often lead to information loss or a failure to properly identify distinct regions of the potential energy surface due to difficulties associated with the high dimensionality of atomic descriptors. In this work, we introduce the Multi-kernel Edge Attention-based Graph Autoencoder (MEAGraph) model, an unsupervised approach for analyzing atomic datasets. MEAGraph combines multiple linear kernel transformations with attention-based message passing to capture geometric sensitivity and enable effective dataset pruning without relying on labels or extensive training. Demonstrated applications on niobium, tantalum, and iron datasets show that MEAGraph efficiently groups similar atomic environments, allowing for the use of basic pruning techniques for removing sampling bias. This approach provides an effective method for representation learning and clustering that can be used for data analysis, outlier detection, and dataset optimization.

Materials science↗
DOE OSTI2024

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Ultrafast solvent-to-solute proton transfer mediated by intermolecular coherent vibrations

Ultrafast photoinduced excited-state proton transfer (ESPT) plays a crucial role in protecting biomolecules and functional materials from photodamage. However, the influence of solute-solvent interactions on ESPT dynamics remains under active investigation. Here, we present an ultrafast spectroscopic study of ESPT in the photobase 2-(2´-pyridyl)benzimidazole (PBI) in methanol. Ultrafast absorption spectroscopy, supported by quantum chemical calculations, reveals three distinct kinetic steps: (1) a 2.2 ps solvent-to-solute proton transfer, (2) subsequent nonradiative relaxation to the ground state within 31 ps, producing a vibrationally hot ensemble with substantial excess kinetic energy, and (3) equilibration as this energy dissipates into the surrounding solvent bath over 186 ps. Femtosecond-resolved dynamics exhibit oscillatory signals indicative of coherent wavepacket motion on the S 1 potential energy surface. A phase flip in the excited-state absorption maximum confirms this assignment. Fourier analysis resolves two dominant periods (∼117 fs and ∼340 fs), corresponding to in-plane and out-of-plane vibrational modes coupled between PBI and the hydrogen-bonded methanol molecule. The rapid dephasing ( < 300 fs) suggests that the nuclear wavefunction evolves on an anharmonic potential energy surface while traversing the ESPT reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Tracking ultrafast non-adiabatic dissociation dynamics of the deuterated water dication molecule

Here, we applied reaction microscopy to elucidate fast non-adiabatic dissociation dynamics of deuterated water molecules after direct photo-double ionization at 61 eV with synchrotron radiation. For the very rare D + + O + + D breakup channel, the particle momenta, angular, and energy distributions of electrons and ions, measured in coincidence, reveal distinct electronic dication states and their dissociation pathways via spin–orbit coupling and charge transfer at crossings and seams on the potential energy surfaces. Notably, we could distinguish between direct and fast sequential dissociation scenarios. For the latter case, our measurements reveal the geometry and orientation of the deuterated water molecule with respect to the polarization vector that leads to this rare 3-body molecular breakup channel. Aided by multi-reference configuration-interaction calculations, the dissociation dynamics could be traced on the relevant potential energy surfaces and particularly their crossings and seams. This approach also unraveled the ultrafast time scales governing these processes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Pt-Assisted Carbon Remediation of Mo2C Materials for CO Disproportionation

Using the CO disproportionation (Boudouard) reaction as a probe reaction, an in-depth analysis of temperature-programmed pulse response data shows that the addition of Pt to Mo2C mitigates deactivation of Mo active sites by acting as a carbon collector. CO2 production on Mo2C and Pt/Mo2C materials is dependent on both the activation energy and the CO surface concentration. Detailed plane-wave density functional theory calculations of the CO adsorption and disproportion reactions on Mo2C-supported Pt nanoparticles (NPs) are reported. The Mo2C was modeled by the ß-Mo2C (100) surface, and the Pt/Mo2C interface was modeled by the addition of 12 Pt atoms to the Mo2C (100) surface (12Pt@Mo2C). The potential energy surfaces of the Boudouard reaction were calculated on pure Mo2C, 12Pt@Mo2C, and Pt (111) surfaces. CO dissociation readily occurs on the Mo2C (100) surface, but not on the Pt (111) surface, with the former being exothermic and the latter being endothermic. At the Pt/Mo2C interface, CO dissociation is still exothermic, but with a larger energy barrier. The Boudouard reaction takes place on the Mo2C region, where CO2 is formed from a surface O atom dissociated from one CO molecule in reaction with another CO molecule, leaving one C atom on the surface. C adsorption is preferential on the Pt site in comparison to the Mo site. The supported Pt domains can collect the remaining C atoms, facilitating further CO2 formation on the active Mo sites. A Bader charge analysis shows that the surface metal-carbon bond is a mixture of covalent and ionic bonds, whereas the surface metal-oxygen bond is ionic. Electron localization function (ELF) and partial charge density calculations agree well with the Bader charge analysis. These computational results are consistent with experimental observations of the interaction of CO with Mo2C nanotube supported Pt domains in the transient regime under far from equilibrium conditions. The Boudouard reaction is an important side reaction, and the unexpected role found for Pt as a carbon collector, with Mo serving as a disproportionation site, provides a unique vantage point for understanding carbon and coke formation on catalytic materials.

Pt-assisted carbon remediation, Pt/Mo2C interface,↗
DOE OSTI2021

Improvements to the macroscopic-microscopic approach of nuclear fission

The well-established macroscopic-microscopic (mac-mic) description of nuclear fission enables the prediction of fission-fragment yields for a broad range of fissioning systems. In this work, we present several key enhancements to this approach. We improve upon the microscopic sector of nuclear potential-energy surfaces by magnifying the resolution of the Lipkin-Nogami equations and strengthening the Strutinsky procedure, thus reducing spurious effects from the continuum. We further present a novel deterministic method for calculating fission dynamics under the assumption of strongly damped nucleonic motion. Our technique directly determines the evolution of the scissioned shape distribution according to the number of random-walk steps rather than the statistical accumulation of fission events. We show that our new technique is equivalent to the Metropolis random walk pioneered over the past decade by Randrup and colleagues. It further improves upon it because we remove the need for altering the nuclear landscape via a biased potential. With our final improvement, we calculate fission fragments mass and charge distributions using particle number projection, which affords the simultaneous calculation of both mass- and charge-yield distributions. Fission fragments are thus calculated from the quantum-mechanical A -body states of the potential-energy surface rather than from the collective mass asymmetry variable α g of the finite-range liquid-drop model used in past work. We highlight the success of our enhancements by predicting the odd-even staggering and the charge polarization for the neutron-induced fission of 233 U and 235 U .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Pairing a Global Optimization Algorithm with EXAFS to Characterize Lanthanide Structure in Solution

Ensemble-average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. Here, to accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO 3 ) 3 dissolved in acetonitrile with a terpyridine ligand (terpyNO 2 ), and Nd(NO 3 ) 3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO 3 ) 3 (terpyNO 2 ) and Nd(NO 3 )(acetonitrile) 3 , with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Measurement of coherent vibrational dynamics with X-ray Transient Absorption Spectroscopy simultaneously at the Carbon K- and Chlorine L 2,3 - edges

X-ray Transient Absorption Spectroscopy (XTAS) is a powerful probe for ultrafast molecular dynamics. The evolution of XTAS signal is controlled by the shapes of potential energy surfaces of the associated core-excited states, which are difficult to directly measure. Here, we study the vibrational dynamics of Raman activated CCl 4 with XTAS targeting the C 1s and Cl 2p electrons. The totally symmetric stretching mode leads to concerted elongation or contraction in bond lengths, which in turn induce an experimentally measurable red or blue shift in the X-ray absorption energies associated with inner-shell electron excitations to the valence antibonding levels. The ratios between slopes of different core-excited potential energy surfaces (CEPESs) thereby extracted agree very well with Restricted Open-Shell Kohn-Sham calculations. The other, asymmetric, modes do not measurably contribute to the XTAS signal. The results highlight the ability of XTAS to reveal coherent nuclear dynamics involving < 0.01 Å atomic displacements and also provide direct measurement of forces on CEPESs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

SpectrumSDT: A program for parallel calculation of coupled rotational-vibrational energies and lifetimes of bound states and scattering resonances in triatomic systems

In this work, we present SpectrumSDT – a program for calculations of energies and lifetimes of bound rotational-vibrational states below and scattering resonances above the dissociation threshold on a global potential energy surface of a triatomic system, which may include stable molecules, weekly-bound van-der-Waals complexes, and unbound atom + diatom scattering systems. Large-amplitude vibrational motion is treated explicitly using hyper-spherical coordinates. Three options for the rotational-vibrational interaction are supported: uncoupled (symmetric top rotor), partially coupled (to include interaction between several nearest states only) and full-coupled (vibrating asymmetric-top rotor). In addition to energies and lifetimes, SpectrumSDT is able to integrate ro-vibrational wave functions over the user-defined regions of potential energy surface, which helps to classify these states. In this release of the code, SpectrumSDT is limited to ABA-type molecules with wave functions that do not extend into the regions near Eckart singularities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

First observation of a shape isomer and a low-lying strongly-coupled prolate band in neutron-deficient semi-magic 187 Pb

Prompt and delayed gamma-ray spectroscopy of the neutron-deficient, semi-magic isotope 187Pb has been performed using the recoil-decay and isomer-decay tagging techniques at the Argonne Gas-Filled Analyzer. A new 5.15(15)-mu s isomeric state at only 308 keV above the spherical 3/2 - ground state is identified and classified as a shape isomer. A strongly-coupled band is observed on top of the isomer, which is nearly identical to the one built on the prolate 7/2 - [514] Nilsson state in the isotone 185 Hg. Based on this similarity and on the result of the potential-energy surface calculations, the new isomer in 187 Pb is proposed to originate from the same configuration. The retarded character of the 308-keV (7/2( - )→ 3/2$_{gs}^{-}$ transition with a deduced B(E2) = 5.6(2) x 10 -4 W.u. can be well explained by the significant difference between the prolate parent and spherical daughter configurations, leading to the shape isomerism. The excitation energy of the isomer is surprisingly low, being roughly half of the excitation energies of the known 0 + intruder bandheads in the neighboring 186,188 Pb isotopes. The combined results of the present work and the previous alpha-decay and laser spectroscopy studies present evidence for triple shape coexistence at low energy in the negative-parity configurations of 187 Pb, which is well reproduced by the potential-energy surface calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Observation of Slow Eigen-Zundel Interconversion in H + (H 2 O) 6 Clusters upon Isomer-Selective Vibrational Excitation and Buffer Gas Cooling in a Cryogenic Ion Trap

The formation of isomers when trapping floppy cluster ions in a temperature-controlled ion trap is a generally observed phenomenon. This involves collisional quenching of the ions initially formed at high temperature by buffer gas cooling until their internal energies fall below the barriers in the potential energy surface that separate them. Here we explore the kinetics at play in the case of the two isomers adopted by the H + (H 2 O) 6 cluster ion that differ in the proton accommodation motif. One of these is most like the Eigen cation with a tricoordinated hydronium motif (denoted E), and the other is most like the Zundel ion with the proton equally shared between two water molecules (denoted Z). After initial cooling to about 20 K in the radiofrequency (Paul) trap, the relative populations of these two spectroscopically distinct isomers are abruptly changed through isomer-selective photoexcitation of bands in the OH stretching region with a pulsed (~6 ns) infrared laser while the ions are in the trap. We then monitor the relaxation of the vibrationally excited clusters and reformation of the two cold isomers by recording infrared photodissociation spectra with a second IR laser as a function of delay time from the initial excitation. The latter spectra are obtained after ejecting the trapped ions into a time-of-flight photofragmentation mass spectrometer, thus enabling long (~0.1 s) delay times. Excitation of the Z isomer is observed to display long-lived vibrationally excited states that are collisionally cooled on a ms time scale, some of which quench into the E isomer. These excited E species then display spontaneous interconversion to the Z form on a ~10 ms time scale. Furthermore, these qualitative observations set the stage for a series of experimental measurements that can provide quantitative benchmarks for theoretical simulations of cluster dynamics and the potential energy surfaces that underlie them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗