Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled cluster”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2023

SPT-CL J2215−3537: A Massive Starburst at the Center of the Most Distant Relaxed Galaxy Cluster

Abstract We present the discovery of the most distant, dynamically relaxed cool core cluster, SPT-CL J2215−3537 (SPT2215), and its central brightest cluster galaxy (BCG) at z = 1.16. Using new X-ray observations, we demonstrate that SPT2215 harbors a strong cool core with a central cooling time of 200 Myr (at 10 kpc) and a maximal intracluster medium cooling rate of 1900 ± 400 M ⊙ yr −1 . This prodigious cooling may be responsible for fueling the extended, star-forming filaments observed in Hubble Space Telescope imaging. Based on new spectrophotometric data, we detect bright [O ii ] emission in the BCG, implying an unobscured star formation rate (SFR) of 320 − 140 + 230 M ⊙ yr −1 . The detection of a weak radio source (2.0 ± 0.8 mJy at 0.8 GHz) suggests ongoing feedback from an active galactic nucleus (AGN), though the implied jet power is less than half the cooling luminosity of the hot gas, consistent with cooling overpowering heating. The extreme cooling and SFR of SPT2215 are rare among known cool core clusters, and it is even more remarkable that we observe these at such high redshift, when most clusters are still dynamically disturbed. The high mass of this cluster, coupled with the fact that it is dynamically relaxed with a highly isolated BCG, suggests that it is an exceptionally rare system that must have formed very rapidly in the early universe. Combined with the high SFR, SPT2215 may be a high- z analog of the Phoenix cluster, potentially providing insight into the limits of AGN feedback and star formation in the most massive galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Quantum mechanical calculations to chemical accuracy

The accuracy of current molecular-structure calculations is illustrated with examples of quantum mechanical solutions for chemical problems. Two approaches are considered: (1) the coupled-cluster singles and doubles (CCSD) with a perturbational estimate of the contribution of connected triple excitations, or CCDS(T); and (2) the multireference configuration-interaction (MRCI) approach to the correlation problem. The MRCI approach gains greater applicability by means of size-extensive modifications such as the averaged-coupled pair functional approach. The examples of solutions to chemical problems include those for C-H bond energies, the vibrational frequencies of O3, identifying the ground state of Al2 and Si2, and the Lewis-Rayleigh afterglow and the Hermann IR system of N2. Accurate molecular-wave functions can be derived from a combination of basis-set saturation studies and full configuration-interaction calculations.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1994

Cr2 Revisited

Abstract Cr2 is studied at the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples [CCSD(T)]. Both spin restricted, RCCSD(T), and spin unrestricted, UCCSD(T), approaches are used. The UCCSD(T) bond length is much longer than the experimental value, while the RCCSD(T) bond length is shorter. This difference is due to the inability of the spin restricted Hartree-Fock reference to describe the antiferromagnetic coupling at larger r values. In contrast with the results obtained using density functional theory in conjunction with older functionals, the results obtained with recently developed functionals are in agreement with the UCCSD(T) results.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2022

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Theoretical Study of the First Transition Row Oxides and Sulfides

The first transition row oxides and sulfides are studied using several different levels of theory. The calculations show the bonding mechanism in the sulfides and oxides to be very similar. For the oxides, accurate experimental data allow the theoretical methods to be calibrated. The same level of theory is used to study the sulfides where there is far less experimental information. For ScO through MnO and CuO the coupled cluster singles and doubles technique including a perturbational estimate of the unliked triple excitations [CCSD(T)] yields spectroscopic constants ((tau)e, (omega)e, and D0) in good agreement with experiment. The triple excitations are found to be very important in achieving this accuracy. For FeO to NiO, the single determinant self-consistent-field (SCF) approach yields pi orbitals that are localized on the metal or oxygen. This appears to cause problems for the single reference techniques; this is discussed in detail for NiO. The complete-active-space SCF/internally contracted averaged coupled pair functional approach (CASSCF/ICACPF) works well for FeO to NiO. The calculation of accurate dipole moments is found to be very difficult.

Bauschlicher, Charles W., Jt.↗
DOE OSTI2024

Emulating ab initio computations of infinite nucleonic matter

We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

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

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Threshold Ionization Spectroscopy and Theoretical Calculations of LnO (Ln = La and Ce)

Mass-analyzed threshold ionization (MATI) spectroscopy was used to measure the vibronic spectra of LnO (Ln = La and Ce). Single-reference coupled cluster and relativistic multireference configuration calculations were carried out to compare with the measured vibronic energies. The spectrum of LaO displays a single vibronic band system, while that of CeO shows multiple ones. The ionization energies of LaO and CeO are measured as 5.2446 (6) and 5.3332(6) eV, respectively, which are a hundredfold improvement over the literature values. The vibrational energies of the neutral molecule and corresponding ion reveal the charge effect on the metal-oxygen bond of both species. The single band system in the spectrum of LaO arises from the transition of the ground state of the neutral molecule with the La(6s 1 )O(2p 6 ) valence configuration to the ground state of the singly charged ion with the La(6s o )O(2p 6 ) configuration. Here, the multiple band systems in the spectrum of CeO are attributed to the spin-orbit coupling for the Ce(4f 1 6s 1 )O(2p 6 ) configuration of the neutral molecule and an excited state for the Ce(4f 1 )O(2p 6 ) configuration of the ion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Excited-State Dynamics during Primary C–I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

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

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗
NASA NTRS1995

The Anharmonic Force Field of Ethylene, C2H4, by Means of Accurate Ab Initio Calculations

The quartic force field of ethylene, C2H4, has been calculated ab initio using augmented coupled cluster, CCSD(T), methods and correlation consistent basis sets of spdf quality. For the C-12 isotopomers C2H4, C2H3D, H2CCD2, cis-C2H2D2, trans-C2H2D2, C2HD3, and C2D4, all fundamentals could be reproduced to better than 10 per centimeter, except for three cases of severe Fermi type 1 resonance. The problem with these three bands is identified as a systematic overestimate of the Kiij Fermi resonance constants by a factor of two or more; if this is corrected for, the predicted fundamentals come into excellent agreement with experiment. No such systematic overestimate is seen for Fermi type 2 resonances. Our computed harmonic frequencies suggest a thorough revision of the accepted experimentally derived values. Our computed and empirically corrected re geometry differs substantially from experimentally derived values: both the predicted rz geometry and the ground-state rotational constants are, however, in excellent agreement with experiment, suggesting revision of the older values. Anharmonicity constants agree well with experiment for stretches, but differ substantially for stretch-bend interaction constants, due to equality constraints in the experimental analysis that do not hold. Improved criteria for detecting Fermi and Coriolis resonances are proposed and found to work well, contrary to the established method based on harmonic frequency differences that fails to detect several important resonances for C2H4 and its isotopomers. Surprisingly good results are obtained with a small spd basis at the CCSD(T) level. The well-documented strong basis set effect on the v8 out-of-plane motion is present to a much lesser extent when correlation-optimized polarization functions are used. Complete sets of anharmonic, rovibrational coupling, and centrifugal distortion constants for the isotopomers are available as supplementary material to the paper.

Martin, Jan M. L.↗
DOE OSTI2024

Elucidating Anomalous Intensity Ratios in Chlorine L-Edge X-ray Absorption Spectroscopy: Multiplet Effects and Core Rydberg Transitions

A relativistic core–valence-separated equation-of-motion coupled cluster (CVS-EOM-CC) study of chlorine L 2,3 -edge X-ray absorption near-edge structure (XANES) spectra using CH 3 Cl and CH 2 ICl as representative molecules is reported. The nearly identical intensity for the main features in the L 2 - and L 3 -edge XANES spectra is attributed to multiplet effects and the overlap between core–valence and core Rydberg transitions. The multiplet effects originating from the interaction between the core hole and the C–Cl σ* orbitals account for around half of the deviation of the L 3 and L 2 intensity ratio from the 2:1 ratio of the numbers of 2p 3/2 and 2p 1/2 electrons. The 2p 3/2 → 4s core Rydberg transitions are shown to overlap with the 2p 1/2 → σ* transitions and contribute to the other half of the intensity anomaly. Here, we demonstrate that triple excitations in CVS-EOM-CC calculations play important roles in accurate simulation of the overlap between the 2p 1/2 → σ* and 2p 3/2 → 4s transitions.

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

Electronic structure of Li 1,2,3 +,0,– and nature of the bonding in Li 2,3 +,0,–

Abstract The current study of the small lithium molecules Li 2 +,0,− and Li 3 +,0,− focuses on the nature of the bonding in these molecules as well as their structures and energetics (bond energies, ionization energies, and electron affinities). Valence CASSCF (2s,2p) calculations incorporate nondynamical electron correlation in the calculations, while the corresponding multireference configuration interaction and coupled cluster calculations incorporate dynamical electron correlation. Treatment of nondynamical correlation is critical for properly describing the Li 2,3 +,0,− molecules as well as the Li − anion with dynamical correlation, in general, only fine‐tuning the predictions. All lithium molecules and ions are bound, with the Li 3 + and Li 2 + ions being the most strongly bound, followed by Li 3 − , Li 2 , Li 2 − and Li 3 . The minimum energy structures of Li 3 +,0,− are, respectively, an equilateral triangle, an isosceles triangle, and a linear structure. The results of SCGVB calculations are analyzed to obtain insights into the nature of the bonding in these molecules. An important finding of this work is that interstitial orbitals, a concept first put forward by McAdon and Goddard in 1985, play an essential role in the bonding of all lithium molecules considered here except for Li 2 . The interstitial orbitals found in the Li 3 +,0 molecules likely give rise to the non‐nuclear attractors/maxima observed in these molecules.

Chemistry↗
DOE OSTI2024

Ab initio exploration of low–lying electronic states of linear and bent MNX + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) and their origins

High-level multireference and coupled cluster quantum calculations were employed to analyze low-lying electronic states of linear-MNX + and side-bonded-M[NX] + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) species. Their full potential energy curves (PECs), dissociation energies (D e s), geometric parameters, excitation energies (T e s), and harmonic vibrational frequencies (ω e s) are reported. The first three chemically bound electronic states of MNX + and M[NX] + are 3 Σ – , 1 Δ, 1 Σ + and 3 A", 1 A', 1 A", respectively. The 3 Σ – , 1 Δ, 1 Σ + of MNX + originate from the M + ( 2 D) + NX( 2 Π) fragments, whereas the 3 A", 1 A', 1A" states of M[NX] + dissociate to M + ( 2 S) + NX( 2 Π) as a result of avoided crossings. The MNX + and M[NX] + are real minima on the potential energy surface and their interconversions are possible. The M 2+ NX – /M 2+ [NX] – ionic structure is an accurate representation for their low-lying electronic states. The D e s of MNX + species were found to depend on the dipole moment (μ) of the corresponding NX ligands and a linear relationship between these two parameters was observed.

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

Evaluation of Density-Functional Tight-Binding Methods for Simulation of Protic Molecular Ion Pairs

In this work, we benchmark the accuracy of the density-functional tight-binding (DFTB) method, namely the long-range corrected second-order (LC-DFTB2) and third-order (DFTB3) models, for predicting energetics of imidazolium-based ionic liquid (IL) ion pairs. We compare the DFTB models against popular density functionals such as LC-ωPBE and B3LYP, using ab initio domain-based local pair-natural orbital coupled cluster (DLPNO-CC) energies as reference. Calculations were carried out in the gas phase, as well as in aqueous solution using implicit solvent methods. We find that the LC-DFTB2 model shows excellent performance in the gas phase and agrees well with reference energies in implicit solvent, often outperforming DFTB3 predictions for complexation energetics. Our study identifies a range of opportunities for use of the LC-DFTB method and quantifies its sensitivity to protonation states and the types of chemical interactions between ion pairs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗