Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled-cluster methods”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Positive Anharmonicities: The Oxonide Anion as an Example

An accurate ab initio quartic force field for the ozonide anion has been determined at the singles and doubles coupled-cluster level of theory that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), using the augmented valence triple-zeta correlation consistent one-particle basis set. Convergence of the harmonic frequencies with respect to the one-particle basis set has been demonstrated by computing quadratic force fields at the CCSD(T) level using augmented valence double-zeta and quadruple-zeta basis sets. Fundamental vibrational frequencies have been determined via second-order ro-vibrational perturbation theory and also using exact variational methods. Agreement is very good, and they both show that the antisymmetric O-O stretch, v 3, possesses a positive anharmonic correction (that is, the fundamental frequency is predicted to be higher in energy than the harmonic frequency). Comparison of the O_3 and O3 quartic force fields shows that the positive anharmonic correction is the result of a particularly large F3333 symmetry internal coordinate force constant. The reasonableness of this force constant has been tested by computing both F33 and F3333 at the CCSD, CCSD(T), and CASPT2 levels of theory. A discussion of known positive anharmonicities for stretching modes is presented.

Lee, Timothy J.↗
DOE OSTI2023

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗
DOE OSTI2022

Systematics of E 2 strength in the s d shell with the valence-space in-medium similarity renormalization group

Recent developments in ab initio nuclear theory demonstrate promising results in medium- to heavy-mass nuclei. A particular challenge for many of the many-body methodologies, however, is an accurate treatment of the electric-quadrupole, E2, strength associated with collectivity. The valence-space in-medium similarity renormalization group (VS-IMSRG) is a particularly powerful method for accessing medium- and high-mass nuclei but has been found to underpredict E2 strengths. The purpose of this work is to evaluate the isospin dependence of this underprediction. We perform a systematic comparison of VS-IMSRG calculations with available literature. We make use of isoscalar and isovector contributions to the E2 matrix elements to assess isoscalar and isovector contributions to the missing strength. It is found that the E2 strength is consistent throughout T z =|12|, T z =|1|, T z =|32|, and T z =2 pairs within the sd shell. Furthermore, no isovector contribution to the deficiency is identified. A comparison with toy-models and coupled-cluster calculations is used to discuss potential origins of the missing strength, which arises from missing many-particle, many-hole excitations out of the model space. The absence of any significant isovector contribution to the missing E2 strength indicates that the E2 strength discrepancy, and therefore any correction, is largely independent of the isospin of the nuclei in question.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

Exploring Parameter Redundancy in the Unitary Coupled-Cluster Ansätze for Hybrid Variational Quantum Computing

One of the commonly used chemical-inspired approaches in variational quantum computing is the unitary coupled-cluster (UCC) ansatze. Despite being a systematic way of approaching the exact limit, the number of parameters in the standard UCC ansatze exhibits unfavorable scaling with respect to the system size, hindering its practical use on near-term quantum devices. Efforts have been taken to propose some variants of UCC ansatze with better scaling. In this paper we explore the parameter redundancy in the preparation of unitary coupled-cluster singles and doubles (UCCSD) ansatze employing spin-adapted formulation, small amplitude filtration, and entropy-based orbital selection approaches. Numerical results of using our approach on some small molecules have exhibited a significant cost reduction in the number of parameters to be optimized and in the time to convergence compared with conventional UCCSD-VQE simulations. Further, we also discuss the potential application of some machine learning techniques in further exploring the parameter redundancy, providing a possible direction for future studies.

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

Reaction dynamics of S( 3 P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S( 3 P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H 2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C 3 H 4 S (loss of ethene or H 2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Lang, Jinxin↗
DOE OSTI2025

Toward a Balanced Description of Ground and Excited States with Transcorrelated F12 Methods

By correlating only the 1-particle states occupied in the reference determinant, the conventional design for the single-reference R12/F12 explicitly correlated methods biases them toward the ground-state description, thereby making the treatment of response properties of the ground state, and energies and other properties of excited states less robust. While the use of multireference methods and/or extensions of the standard SP-projected geminals can achieve a more balanced description of ground and excited states, here we show that the same goals can be achieved by extending the action of F12 correlators to the occupied and valence unoccupied 1-particle states only. This design choice reflects the strong dependence of the optimal correlation length scale of the F12 ansatz on the orbital energies/structure, and helps to avoid the unphysical raising of the ground-state energy if the F12 geminals are used to correlate pairs of all 1-particle states. The improved F12 geminal design is incorporated into the unitary transcorrelation framework to produce a unitary 2-body Hamiltonian that incorporates the short-range dynamical correlation physics for ground and low-energy excited states in a balanced manner. This explicitly correlated effective Hamiltonian reduces the basis set requirement on the correlation-consistent basis cardinal number by 1 or more over the uncorrelated counterpart for the description of ground-state coupled-cluster singles and doubles (CCSD) energies, vertical excitation energies, and harmonic vibrational frequencies of equation-of-motion CCSD low-energy excited states.

Hamiltonians↗
DOE OSTI2022

On the prospects of optical cycling in diatomic cations: effects of transition metals, spin–orbit couplings, and multiple bonds

Molecules with optical cycling centers (OCCs) are highly desirable in the context of fundamental studies as well as applications (e.g., quantum computing) because they can be effectively cooled to very low temperatures by repeated absorption and emission (hence, cycling). Charged species offer additional advantages for experimental control and manipulation. Here, we present a systematic computational study of a series of diatomic radical-cations made of a d-block metal and a p-block ligand, that are isoelectronic (in their valence shell) to the successfully laser-cooled neutral molecules. Using high-level electronic structure methods, we characterize state and transition properties of low-lying electronic states and compute Franck-Condon factors. The computed branching ratios and radiative lifetimes reveal that the electronic transitions analogous to those successfully used in the laser cooling of neutral molecules are less than optimal in the cations. We propose alternative transitions suitable for optical cycling and highlight trends that could assist future designs of OCCs in charged or neutral molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

A classical model for three-body interactions in aqueous ionic systems

We present a classical induction model to evaluate the three-body ion–water–water (I–W–W) and water–water–water (W–W–W) interactions in aqueous ionic systems. The classical description of the induction energy is based on electrostatic distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole on the O and H atoms of water. The monatomic ions were described by a point charge and a dipole–dipole polarizability, while for the polyatomic ions, distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole were used. The accuracy of the classical model is benchmarked against an accurate dataset of 936 (I–W–W) and 2184 (W–W–W) three-body terms for 13 different monatomic and polyatomic cation and anion systems. The classical model shows excellent agreement with the reference second order Moller–Plesset and coupled-cluster single double and perturbative triple [CCSD(T)] three-body energies. The Root-Mean-Square-Errors (RMSEs) for monatomic cations, monatomic anions, and polyatomic ions were 0.29, 0.25, and 0.12 kcal/mol, respectively. The corresponding RMSE for 1744 CCSD(T)/aVTZ three-body (W–W–W) energies, used to train MB-pol, was 0.12 kcal/mol. The accuracy of the proposed classical model demonstrates that the three-body term for aqueous ionic systems can be accurately modeled classically. This approach provides a fast, efficient, and as-accurate path toward modeling the three-body term in aqueous ionic systems that is fully transferable across systems with different ions without the need to fit to tens of thousands of ab initio calculations for each ion to extend existing many-body force fields to interactions between water and ions.

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

High-resolution CH stretch spectroscopy of jet-cooled cyclopentyl radical: First insights into equilibrium structure, out-of-plane puckering, and IVR dynamics

First, high-resolution sub-Doppler infrared spectroscopic results for cyclopentyl radical (C 5 H 9 ) are reported on the α-CH stretch fundamental with suppression of spectral congestion achieved by adiabatic cooling to T rot ≈ 19(4) K in a slit jet expansion. Surprisingly, cyclopentyl radical exhibits a rotationally assignable infrared spectrum, despite 3N – 6 = 36 vibrational modes and an upper vibrational state density (ρ ≈ 40–90 #/cm –1 ) in the critical regime (ρ ≈ 100 #/cm –1 ) necessary for onset of intramolecular vibrational relaxation (IVR) dynamics. Such high-resolution data for cyclopentyl radical permit detailed fits to a rigid-rotor asymmetric top Hamiltonian, initial structural information for ground and vibrationally excited states, and opportunities for detailed comparison with theoretical predictions. Specifically, high level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples (CCSD(T))/ANO0, 1 level are used to calculate an out-of-plane bending potential, which reveals a C 2 symmetry double minimum 1D energy surface over a C 2v transition state. Here, the inversion barrier [V barrier ≈ 3.7(1) kcal/mol] is much larger than the effective moment of inertia for out-of-plane bending, resulting in localization of the cyclopentyl wavefunction near its C 2 symmetry equilibrium geometry and tunneling splittings for the ground state too small (<1 MHz) to be resolved under sub-Doppler slit jet conditions. The persistence of fully resolved high-resolution infrared spectroscopy for such large cyclic polyatomic radicals at high vibrational state densities suggests a “deceleration” of IVR for a cycloalkane ring topology, much as low frequency torsion/methyl rotation degrees of freedom have demonstrated a corresponding “acceleration” of IVR processes in linear hydrocarbons.

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

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗
DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

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

“Best” Iterative Coupled-Cluster Triples Model? More Evidence for 3CC

To follow up on the unexpectedly good performance of several coupled-cluster models with approximate inclusion of 3-body clusters we performed a more complete assessment of the 3CC method for accurate computational thermochemistry in the standard HEAT framework. New spin-integrated implementation of the 3CC method applicable to closed- and open-shell systems utilizes a new automated toolchain for derivation, optimization, and evaluation of operator algebra in many-body electronic structure. We found that with a double-ζ basis set the 3CC correlation energies and their atomization energy contributions are almost always more accurate (with respect to the CCSDTQ reference) than the CCSDT model as well as the standard CCSD(T) model. The mean absolute errors in cc-pVDZ {3CC, CCSDT, and CCSD(T)} electronic (per valence electron) and atomization energies relative to the CCSDTQ reference for the HEAT data set, were {24, 70, 122} μE h /e and {0.46, 2.00, 2.58} kJ/mol, respectively. The mean absolute errors in the complete-basis-set limit {3CC, CCSDT, and CCSD(T)} atomization energies relative to the HEAT model reference, were {0.52, 2.00, and 1.07} kJ/mol, The significant and systematic reduction of the error by the 3CC method and its lower cost than CCSDT suggests it as a viable candidate for post- CCSD(T) thermochemistry applications, as well as the preferred alternative to CCSDT in general.

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

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

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

The many-body expansion for metals. I. The alkaline earth metals Be, Mg, and Ca

We examine the many-body expansion (MBE) for alkaline earth metal clusters, Be n , Mg n , Ca n ( n = 4, 5, 6), at the Møller–Plesset second order perturbation theory, coupled-cluster singles and doubles with perturbative triples, multi-reference perturbation theory, and multi-reference configuration interaction levels of theory. The magnitude of each term in the MBE is evaluated for several geometrical configurations. We find that the behavior of the MBE for these clusters depends strongly on the geometrical arrangement and, to a lesser extent, on the level of theory used. Another factor that affects the MBE is the in situ (ground or excited) electronic state of the individual atoms in the cluster. For most geometries, the three-body term is the largest, followed by a steady decrease in absolute energy for subsequent terms. Though these systems exhibit non-negligible multi-reference effects, there was little qualitative difference in the MBE when employing single vs multi-reference methods. Useful insights into the connectivity and stability of these clusters have been drawn from the respective potential energy surfaces and quasi-atomic orbitals for the various dimers, trimers, and tetramers. Through these analyses, we investigate the similarities and differences in the binding energies of different-sized clusters for these metals.

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

Converging high-level coupled-cluster energetics via adaptive selection of excitation manifolds driven by moment expansions

A novel approach to rapidly converging high-level coupled-cluster (CC) energetics in an automated fashion is proposed. The key idea is an adaptive selection of excitation manifolds defining higher--than--two-body components of the cluster operator inspired by CC(P;Q) moment expansions. Further, the usefulness of the resulting methodology is illustrated by molecular examples where the goal is to recover the electronic energies obtained using the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT) when the noniterative triples corrections to CCSD fail.

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

Direct reactions with the AT-TPC

Direct reactions are crucial tools for accessing properties of the atomic nucleus. Fundamental and exotic phenomena such as collective modes, pairing, weakbinding effects and evolution of single-particles energies can be investigated in peripheral collisions between a heavy nucleus and a light target. The necessity of using inverse kinematics to reveal how these structural properties change with isospin imbalance renders direct reactions a challenging technique when using the missing mass method. In this scenario, Active Target Time Projection Chambers (AT-TPC) have demonstrated an outstanding performance in enabling these types of reactions even under conditions of very low beam intensities. The AT-TPC of the Facility for Rare Isotope Beams (FRIB) is a next generation multipurpose Active Target. When operated inside a solenoidal magnet, direct reactions benefit from the measurement of the magnetic rigidity that enables particle identification and the determination of the excitation energy with high resolution without the need of auxiliary detectors. Additionally, the AT-TPC can be coupled to a magnetic spectrometer improving even further its spectroscopic investigation capability. In this contribution, we discuss inelastic scattering and transfer reaction data obtained via the AT-TPC and compare them to theory. In particular, we present the results for the 14 C(p,p′) and 12 Be (p,d) 11 Be reactions. For 14 C, we compare the experimental excitation energy of the first 1 – excited state with coupled-cluster calculationsbased on nuclear interactions from chiral effective field theory and with available shell-model predictions. For 12 Be, we determine the theoretical spectroscopic factors of the 12 Be (p,d) 11 Be transfer reaction in the shell modeland compare them to the experimental excitation spectrum from a qualitative standpoint.

active target↗
NASA NTRS1990

Vibrations in small Mg clusters

The equilibrium geometries, binding energies, and harmonic frequencies of Mg3 and Mg4 have been determined using large atomic natural orbital basis sets in conjunction with high levels of electron correlation. The correlation treatments comprise multireference configuration interaction (MRCI), singles and doubles coupled-cluster (CCSD) theory and the CCSD(T) extension that includes a perturbational estimate of connected triple excitations. As with previous studies of small Be clusters, the CCSD(T) method is found to reproduce the MRCI results with a remarkable degree of accuracy and at a fraction of the computational cost. Using the CCSD(T) method, full quartic force fields for Mg3 and Mg4 have been determined and anharmonic analyses have been performed using second-order perturbation theory. Vibration-rotation interaction constants and centrifugal distortion constants have also been determined.

Lee, Timothy J.↗