SEARCH · Search NASA

Results for “COUPLING”

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.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Functional group activation and coupling agent migration induced by plasma treatment in adhesive for enhanced toughness of metal-composite joints

Commercial adhesive paste was additionally activated using low-power oxygen plasma with different plasma exposure durations to investigate its interfacial bonding contribution in adhesively-bonded Aural-5/CFRP-PA6 double cantilever beam (DCB) joints with plasma-treated adherends. The plasma-treated adhesive showed enhanced functional peak intensities in Fourier transfer infrared spectroscopy (FTIR) as plasma exposure time was increased, and some degree of oxidation was also detected via X-ray photoelectron spectroscopy (XPS). A small molecule, dicyandiamde (coupling agent) in adhesive was migrated to adhesive layer, where double concentration of coupling agent was detected after 10 minutes of plasma treatment. As a result, characterization of the adhesive/CFRP-PA6 interface after DCB fracture showed additional chemical bonds formed through amine-initiated epoxy ring opening polymerization, as well as amide bonds and ester bonds. In conclusion, these additional bonds at the adhesive/CFRP-PA6 interfaces led to significantly increased failure extensions and fracture energies of the joints as the plasma exposure time on the adhesive increased, compared to the joints with plasma-treated adherends but non-plasma-treated adhesive as commonly seen in the literature.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Quantifying Phospholipids in Organic Samples Using a Hydrophilic Interaction Liquid Chromatography–Inductively Coupled Plasma High-Resolution Mass Spectrometry (HILIC-ICP-HRMS) Method

Here, in this study, a novel method using hydrophilic interaction liquid chromatography (HILIC) coupled with inductively coupled plasma high-resolution mass spectrometry (ICP-HRMS) was introduced for the quantification of phospholipids in oil samples. The method employed a bridged ethyl hybrid (BEH) stationary phase HILIC column with a tetrahydrofuran (THF)/water mobile phase, enhancing the solubility and detection of phospholipids. During the study, a gradient/matrix effect on ICP-HRMS sensitivity was observed and successfully compensated for experimentally, ensuring reliable quantification results. This approach has proven effective for a wide range of different oil samples including vegetable oils, animal fats, and phospholipid supplements. Notably, this method allowed the direct quantification of phospholipids in oil samples, bypassing the need for prior sample preparation methods, such as solid phase extraction (SPE), thereby streamlining the analytical process. The precision, accuracy, and reduced need for extensive sample preparation offered by this method mark a significant advancement in lipids analysis. Its robustness and broad applicability have substantial implications for industries such as food and renewable energy production, where both efficient and accurate lipid identification and quantification are crucial.

09 BIOMASS FUELS↗

DOE OSTI2026

Single-Particle Isotopic Analysis Using the Liquid Sampling-Atmospheric Pressure Glow Discharge Microplasma Coupled to an Orbitrap Mass Spectrometer

Isotope ratio (IR) determinations on individual particles can provide valuable information relative to the sample, including processing and dating, which could be valuable to geological and nuclear material analysis communities. In comparison to “bulk” measurements, in which all particle information is lost and homogenized within a sample, “particle” measurements allow for fingerprinting and can provide unique insights related to the sample’s nano- and micro- compositions. Furthermore, the complex nature of particles, with potential isobaric interferences from either the matrix or other elements within the same particle, poses a significant analytical challenge for conventional mass spectrometry platforms employed for elemental analysis due to their limited mass resolution. Presented here is the first demonstration of an ultrahigh resolution method for single particle (SP) isotopic analysis using the liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma ionization source coupled to an Orbitrap mass spectrometer and FTMS Booster X2T acquisition/processing system. For proof of concept, well-characterized CeO 2 microparticles with a nominal diameter of 1.0 μm were analyzed at a mass resolution of ∼330,000. Important instrumental parameters were optimized to enable the detection of single particles. The 142 Ce/ 140 Ce ratio of individual particles and the population average were determined and compared to an SP inductively coupled plasma time-of-flight mass spectrometry (ICP-TOF-MS) analysis. The isotopic accuracy and precision were determined by two methods and found to be in good agreement with the SP-ICP-TOF-MS method, with the linear regression slope method providing a more accurate ratio, showing a ∼1.4% relative difference from the ratio determined from a particle digest. The encouraging performance of the method was further supported by a determined detection limit of 2.9 fg ( 142 Ce). The developed method is anticipated to overcome many of the challenges posed by isobaric interferences encountered in conventional mass spectrometry techniques when analyzing real-world samples for particle isotopic composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗

DOE OSTI2025

Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems

Coupled-cluster theory (CC) is a highly accurate and versatile method for simulating complex interactions within quantum systems. The extension of CC theory to model mixed electron-photon processes with quantum electrodynamics (QED) has improved our capability to predict cavity-modified chemistry, a field where photons are used as cost-effective and eco-friendly alternatives to catalyze/inhibit chemical reactions. However, calculations with CC methods, even without incorporating QED effects, are often prohibitively expensive. Simulations of larger systems require scalable infrastructures that exist for traditional CC methods but not for QED-CC methods. As such, we present a GPU-enabled, high-performance, open-source implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) within the ExaChem quantum chemistry software package. ExaChem relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure: a parallel heterogeneous tensor library designed to achieve scalable performance on modern heterogeneous supercomputing platforms. Furthermore, we discuss theoretical foundations, algorithmic details, and numerical benchmarks to showcase the larger systems that ExaChem can simulate and how the integration of photonic degrees-of-freedom alters their ground-state properties.

Basis sets↗

DOE OSTI2026

Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Local Pair Natural Orbital-Based Coupled-Cluster Theory through Full Quadruples (DLPNO–CCSDTQ)

In this work, we implement a local pair natural orbitalbased coupled-cluster method through the full treatment of quadruple excitations (CCSDTQ). The domain-based local pair natural orbital (DLPNO) approach, which has successfully been applied to lower levels of coupled-cluster theory, is utilized in our algorithm, and thus our algorithm is called DLPNO-CCSDTQ. For simplicity in the working equations and in the implementation, we t 1 -dress the twoelectron integrals as well as Fock matrix elements. Our method can recover CCSDTQ-CCSDT and CCSDTQ-CCSDT(Q) energy differences on the order of 0.01−0.05 kcal mol −1 , even at a loose quadruples natural orbital (QNO) occupation number cutoff of 3.33 × 10 −6 . To highlight the capabilities of our code and its potential future applications, we showcase computations that would be intractable with canonical CCSDTQ, such as the benzene dimer, (H 2 O) 17 , and adamantane. With sufficient computing resources, computations up to 15 heavy atoms (40 atoms overall) may be feasible for fully bonded 3D systems.

Cluster chemistry↗

DOE OSTI2026

One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory

In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗

DOE OSTI2025

Synthesis of Sandwich-α-Diimine Ligands via Copper-Mediated Stille Coupling

While sandwich diimines have found use in transition metal catalysis, the general method for their preparation has been lacking. Here, this approach utilizes a copper-mediated Stille-type coupling between aryl stannanes and α-diimines bearing 8-bromonaphthylimino group in commercial-grade N,N-dimethylformamide solvent and employs CuI or thiophene-2-carboxylate to promote the reactions. No other additives or bases are required. Electron-poor, electron-rich, heterocyclic, and hindered aryl stannanes can be employed as coupling reagents. Ligand scaffolds derived from structurally diverse α-diketones can be used.

Catalysts↗

DOE OSTI2025

Activation, Dehydrogenation, and Carbon–Carbon Coupling of Methane by Iridium Cations Studied by Infrared Multiple Photon Dissociation Spectroscopy and Density Functional Theory

Products resulting from the sequential activation of one, two, three, and four methane molecules by atomic iridium cations were characterized by gas-phase infrared multiple photon dissociation spectroscopy and density functional theory (DFT) calculations. Iridium cations were generated using a laser ablation source and reacted with methane in a linear radiofrequency ion trap before mass analysis and spectroscopic interrogation in a Fourier transform ion cyclotron resonance mass spectrometer coupled to the free-electron laser for intracavity experiments (FELICE) beamline. Product ions were irradiated using infrared light over the 250–1500 cm –1 range. Comparisons between the experimental and DFT-calculated spectra enabled structural determination of the products formed. The observed products include HIrCH + , Ir(CH 2 ) 2 + , H s Ir(C 3 H 5 ) + , and Ir(CH 3 ) s (C 3 H 5 ) + , where the subscript s denotes a syn orientation of the two ligands. Furthermore, formation of the latter two products provides evidence for efficient C–H bond activation and subsequent C–C coupling on the atomic iridium cation.

Chemical reactions↗

DOE OSTI2025

Coupling of Charge Regulation and Geometry in Soft Ionizable Molecular Assemblies

The size, shape, and charge of structures, such as proteins and amphiphile assemblies, respond in an interconnected manner to solution ionic conditions. Here, we analyze assemblies of an amphiphile (C 16 K 2 ), with two ionizable amino acids [lysine (K)] coupled to a 16-carbon alkyl tail, via small-angle X-ray scattering (SAXS), nonlinear Poisson–Boltzmann theory (nl-PB), and hybrid Monte Carlo-molecular dynamics (MC-MD) simulations. SAXS revealed structural transitions from spherical micelles to cylindrical micelles to bilayers with increasing pH. By combining SAXS-determined structural information and nl-PB, we derived the molecular degree of ionization as a function of pH. The back-calculated titration curves matched the experimental data over an extended pH range, without adjustable parameters. Similarly, the SAXS data on the evolution of spherical micelle structure with ionic strength were combined with nl-PB and MC-MD to derive the bare and effective charges. MC-MD, which considered finite ion sizes, showed that bare and effective charges saturate quickly with increasing salt concentration. Furthermore, the calculated effective charges closely matched results from Zeta-potential measurements. The presented approach has advantages over customary methods for charge regulation, such as the Henderson–Hasselbalch (HH) or Hill models, where molecular ionization/deionization in assemblies is described by effective pKs that are distinct from the pK for isolated molecules. However, these models lack a physical explanation for these pK shifts. By contrast, our approach of combining structural details with an electrostatic model and simulations provides a more intuitive understanding of structure-charge coupling and a framework for understanding charge regulation in many synthetic and biological systems.

ionization↗

DOE OSTI2024

Multiconfigurational Electronic Structure of Nickel Cross-Coupling Catalysts Revealed by X-ray Absorption Spectroscopy

Ni II 2,2'-bipyridine complexes are commonly invoked intermediates in metallaphotoredox cross-coupling reactions. Despite their ubiquity, design principles targeting improved catalytic performance remain underdetermined. A series of Ni( R bpy)( R 'Ar)Cl (R = MeOOC, t-Bu, R' = CH 3 , CF 3 ) complexes were proposed to have multiconfigurational electronic structures on the basis of multiconfigurational/multireference calculations, with significant mixing of Ni → bpy metal-to-ligand charge transfer (MLCT) configurations into the ground-state wave function. Here, Ni K-edge and L 2,3 -edge X-ray absorption spectroscopies provide experimental support for the highly covalent and multiconfigurational electronic structures of these complexes. The pre-edge intensity in the K-edge spectrum reflects highly covalent Ni–aryl bonding. The L 3 -edge spectral shape is dependent on ligand functionalization, and a feature reflecting the MLCT character is assigned using prior ab initio and new semiempirical calculations. The results suggest the push/pull effects of the aryl/bpy ligands moderate the changes in electron density on Ni during the multiredox cross-coupling reaction cycle.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Binding the Power of Cycloaddition and Cross-Coupling in a Single Mechanism: An Unexpected Bending Journey to Radical Chemistry of Butadiynyl with Conjugated Dienes

What if an experiment could combine the power of cycloaddition and cross-coupling with the in situ formation of an aromatic molecule in a single collision? Crossed molecular beam experiments augmented with electronic structure and statistical calculations provided compelling evidence on a novel radical route involving 1,3-butadiynyl (HCCCC; X 2 Σ + ) radicals synthesizing (substituted) arylacetylenes in the gas phase upon reactions with 1,3-butadiene (CH 2 CHCHCH 2 ; X 1 Ag) and 2-methyl-1,3-butadiene (isoprene; CH 2 C(CH 3 )CHCH 2 ; X 1 A'). Furthermore, this elegant mechanism de facto merges two previously disconnected concepts of cross-coupling and cycloaddition–aromatization in a single collision event via the formation of two new C(sp 2 )–C(sp 2 ) bonds and bending the 180° moiety of the linear 1,3-butadiynyl radical out of the ordinary by 60° to 120°. In addition to its importance to fundamental organic chemistry, this unconventional mechanism links two previously separated routes of gas-phase molecular mass growth processes of polyacetylenes and polycyclic aromatic hydrocarbons (PAHs), respectively, in low-temperature environments such as in cold molecular clouds like the Taurus Molecular Cloud (TMC-1) and in hydrocarbon-rich atmospheres of planets and their moons such as Titan, what revises the established understanding of low-temperature molecular mass growth processes in the Universe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Raman Fingerprints of Phase Transitions and Ferroic Couplings in van der Waals Multiferroic CuCrP 2 S 6

CuCrP 2 S 6 (CCPS), a type-II multiferroic material, exhibits unique phase transitions involving ferroelectric, antiferroelectric, and antiferromagnetic ordering. Here, in this study, we conduct a comprehensive investigation on the intricate phase transitions and their multiferroic couplings in CCPS across a wide temperature range from 4 to 345 K through Raman spectroscopic measurements down to 5 cm –1 . We first assign the observed Raman modes with the support of theoretical calculations and angle-resolved polarized Raman measurements. We further present clear signatures of phase transitions from the analyses of temperature-dependent Raman spectral parameters. Particularly, two low-frequency soft modes are observed at 36.1 cm –1 and 70.5 cm –1 below 145 K, indicating the antiferroelectric to quasi-antiferroelectric transition. Moreover, phonon mode hardening is observed when the temperature increases from 4 to 65 K, suggesting negative thermal expansion (NTE) and strong magnetoelastic coupling below 65 K. These findings advance the understanding of vdW multiferroic CCPS, paving the way for future design and engineering of multiferroicity in cutting-edge technologies, such as spintronics and quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗