Search NASA⌕ Search

SEARCH · Search NASA

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 289 records · Page 16

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↗

CASM — A software package for first-principles based study of multicomponent crystalline solids

CASM is a software package that enables first-principles based studies of crystalline materials. It has been designed to treat coupled chemical, mechanical, vibrational, and magnetic degrees of freedom to determine ground state and finite temperature properties of crystals. The symmetry of the underlying parent crystal structure is used to enumerate perturbations of the parent crystal structure and generate derivative structures which can be input to first-principles calculations in order to explore the ground state energy landscape. CASM algorithmically constructs cluster expansions that fully couple discrete and continuous degrees of freedom, and generates highly efficient code to evaluate the cluster expansion basis functions. Widely used machine learning methods are integrated for fitting expansion coefficients to first-principle calculations. The fully parameterized cluster expansions can be combined with (kinetic) Monte Carlo methods to calculate finite temperature thermodynamic and kinetic properties. CASM Alloy Manager identifies distinct parent crystal structures in an alloy system, creating individual projects for each, and integrating the results. The integrated infrastructure facilitates the linkage between first-principles statistical mechanics predictions with higher length scale computational methods, such as phase field simulations. In conclusion, CASM projects are designed to be easy to share in repositories, to re-use, and to extend to include additional chemical species or types of degrees of freedom.

36 MATERIALS SCIENCE↗

Multi-scale modeling of fuel fragmentation and microstructural evolution

An overview of NEAMS program activities to model fuel fragmentation and microstructual evolution to a group of industry attendees at a meeting hosted by EPRI. Topics include atomistic methods to determine grain boundary strength, phase-field fracture modeling to determine critical pressure at grain boundaries, and coupled phase-field-cluster dynamics to model microstructural evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Clustering-based biased Monte Carlo Approach to Protein Titration Curve Prediction

We develop and implement a novel approach to computing the ensemble averages in systems characterized by pair-wise interactions between the entities. Methods involving full enumeration of the configuration space result in exponential complexity. Sampling methods such as Markov Chain Monte Carlo (MCMC) algorithms have been proposed to tackle the exponential complexity of these problems. In certain scenarios where significant energetic coupling exists between the entities, the accuracy of the such algorithms can be diminished. We propose a strategy to improve the accuracy of the MCMC runs by taking advantage of the cluster structure in the interaction energy matrix. We propose two different schemes for performing the biased MCMC runs on the partitioned systems and show that they are valid MCMC schemes. We then apply these algorithms to the problem of computing the protonation fractions and hence the titration curves of titratable protein residues that constitute a given protein. We leverage both synthesized and real-world systems and show the improved performance of our biased MCMC methods when compared to the regular MCMC method.

Visweswara Sathanur, Arun↗

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↗

Magnetic phase diagram of a two-orbital model for bilayer nickelates with varying doping

Motivated by the recently discovered high-T c bilayer nickelate superconductor La 3⁢ Ni 2 ⁢O 7 , we comprehensively research a bilayer 2×2×2 cluster for different electronic densities n by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying n. Based on the spin structure factor S(q), we have obtained a rich magnetic phase diagram in the plane defined by n and U/W, at fixed Hund coupling, where U is the Hubbard strength and W the bandwidth. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. At half-filling, n=2 (two electrons per Ni site, corresponding to N=16 electrons), the canonical superexchange interaction leads to a robust G-AFM state (π,π,π) with antiferromagnetic couplings both in-plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the “half-empty” and “half-full” mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At n=1.5 (or N=12), density corresponding to La 3 ⁢Ni 2 ⁢O 7 , we obtained the “Stripe 2” ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the z axis) in the 2×2×2 cluster. In addition, we obtained a much stronger AFM coupling along the z axis than the magnetic coupling in the xy plane. The random-phase approximation calculations with varying n give very similar results as Lanczos, even though both techniques are based on quite different procedures. Additionally, a state with q/π=(0.6,0.6,1) close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to n=1.25, possibly responsible for the E-phase SDW recently observed in experiments. In conclusion, our predictions can be tested by chemically doping La 3 ⁢Ni 2 ⁢O 7 .

36 MATERIALS SCIENCE↗

Iron(III)–Oxo Cluster Chemistry with Dimethylarsinate Ligands: Structures, Magnetic Properties, and Computational Studies

A program has been initiated to develop Fe III /oxo cluster chemistry with the ‘pseudo-carboxylate’ ligand dimethylarsinate (Me 2 AsO 2 - ) for comparison with the well investigated Fe III /oxo/carboxylate cluster area. The synthesis and characterization of three polynuclear Fe III complexes are reported, [Fe 12 O 4 (O 2 C t Bu) 8 (O 2 AsMe 2 ) 17 (H 2 O) 3 ]Cl 3 (1), Na 2 [Fe 12 Na 2 O 4 (O 2 AsMe 2 ) 20 (NO 3 ) 6 (Me 2 AsO 2 H) 2 (H 2 O) 4 ](NO 3 ) 6 (2) and [Fe 3 (O 2 AsMe 2 ) 6 (Me 2 AsO 2 H) 2 (hqn) 2 ](NO 3 ) (3), where hqn is 8-hydroxyquinoline. The Fe 12 core of 1 is a type never previously encountered in Fe III carboxylate chemistry, consisting of two Fe 6 units each of which comprises two {Fe 3 (μ 3 -O 2- )} units bridged by three Me 2 AsO 2 - groups and linked into an Fe 12 loop structure by two anti-anti η 1 :η 1 :μ Me 2 AsO 2 - groups, a bridging mode extremely rare with carboxylates. 2 also consists of two Fe 6 units, differing in their ligation from those in 1, and this time linked together into a linear structure by a central {Na 2 (NO 3 ) 2 } bridging unit. 3 is a linear Fe 3 complex with no monoatomic bridges between Fe III ions, a very rare situation in Fe III chemistry with any ligands, and unprecedented in Fe carboxylate chemistry. The distinct differences observed in arsinate vs carboxylate ligation modes are rationalized largely based on the greater basicity of the former vs the latter. Variable-temperature dc and ac magnetic susceptibility data reveal all Fe 2 pairwise interactions to be antiferromagnetic. For 1 and 2, the different J ij couplings were estimated by use of a magnetostructural correlation for high nuclearity Fe III -oxo clusters and by density functional theory calculations using broken symmetry methods, allowing identification of their relative spin vector alignments and thus rationalization of their S = 0 ground states. The J ij values were then used as input values to give excellent fits of the experimental χM T vs T data. For 3, the fits of the experimental χM T vs T data to the Van Vleck equation or with PHI gave a very weak J 12 = -0.8(1) cm -1 (H = –2JŠ i ·Š j convention) between adjacent Fe III ions, and an S = 5/2 ground state. Furthermore, these initial Fe III arsinate complexes also provide structural parameters that help validate literature assignments of arsinate binding modes to iron oxide/hydroxide minerals as part of environmental concerns of using arsenic-containing herbicides in agriculture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Synthesis of Co‐Doped MoS 2 Monolayers with Enhanced Valley Splitting

Abstract Internal magnetic moments induced by magnetic dopants in MoS 2 monolayers are shown to serve as a new means to engineer valley Zeeman splitting (VZS). Specifically, successful synthesis of monolayer MoS 2 doped with the magnetic element Co is reported, and the magnitude of the valley splitting is engineered by manipulating the dopant concentration. Valley splittings of 3.9, 5.2, and 6.15 meV at 7 T in Co‐doped MoS 2 with Co concentrations of 0.8%, 1.7%, and 2.5%, respectively, are achieved as revealed by polarization‐resolved photoluminescence (PL) spectroscopy. Atomic‐resolution electron microscopy studies clearly identify the magnetic sites of Co substitution in the MoS 2 lattice, forming two distinct types of configurations, namely isolated single dopants and tridopant clusters. Density functional theory (DFT) and model calculations reveal that the observed enhanced VZS arises from an internal magnetic field induced by the tridopant clusters, which couples to the spin, atomic orbital, and valley magnetic moment of carriers from the conduction and valence bands. The present study demonstrates a new method to control the valley pseudospin via magnetic dopants in layered semiconducting materials, paving the way toward magneto‐optical and spintronic devices.

Zhou, Jiadong↗

Helium droplet infrared spectroscopy of the butyl radicals

Butyl radicals (n-, s-, i-, and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo-tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700 to 3125 cm -1 , allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800–3000 cm -1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photoluminescence, photophysics, and photochemistry of the V B - defect in hexagonal boron nitride

Extensive photochemical and spectroscopic properties of the V B - defect in hexagonal boron nitride are calculated, concluding that the observed photoemission associated with recently observed optically detected magnetic resonance is most likely of (1) 3 E'' → (1) 3 A' 2 origin. Rapid intersystem crossing from the defect's triplet to singlet manifolds explains the observed short excited-state lifetime and very low quantum yield. New experimental results reveal smaller intrinsic spectral bandwidths than previously recognized, interpreted in terms of spectral narrowing and zero-phonon-line shifting induced by the Jahn-Teller effect. Different types of computational methods are applied to map out the complex triplet and singlet defect manifolds, including the doubly ionized formulation of the equation-of-motion coupled-cluster theory that is designed to deal with the open-shell nature of defect states, and mixed quantum-mechanics/molecular-mechanics schemes enabling 5763-atom simulations. Finally, two other energetically feasible spectral assignments from among the singlet and triplet manifolds are considered, but ruled out based on inappropriate photochemical properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Rethinking the magnetic properties of lepidocrocite: A density functional theory and cluster expansion study

The iron oxyhydroxide lepidocrocite (g-FeOOH) is an abundant mineral critical to a number of chemical and technological applications. Of particular interest is the ground state and nite temperature magnetic order, and the subsequent impact this has upon crystal properties. The magnetic properties, investigated in this work are governed primarily through superexchange interactions, and have been calculated using density functional theory and cluster expansion methods. Quantication of these exchange terms has facilitated the determination of the ground state magneto-crystalline structure and subsequent calculation of its lattice constants, elastic moduli, cohesive enthalpy, and electronic density of states. Further, using a collinear magnetic conguration model, the magnetic heat capacity versus temperature has been studied and the Neel temperature obtained.

Pope, Daniel J.↗

Sequential surface synthesis of dispersed sub-nanometer iridium on titanium nitride for acidic water oxidation

Maximizing iridium utilization while maintaining high oxygen evolution reaction (OER) performance remains a persistent challenge in acidic water electrolysis. Immobilizing Ir on conductive, acid-stable supports is promising, yet simultaneously achieving sub-nanometer size, high area coverage, and strong electronic coupling is difficult. Here, we report a sequential surface synthesis on titanium nitride (TiN) that yields uniformly distributed sub-nanometer Ir arrays (∼0.7 nm). Our method uses ethylenediaminetetraacetic acid (EDTA) as a temporal scaffold: it chemisorbs to TiN to install dense chelating sites, captures Ir 3+ ions, and confines Ir cluster growth. A subsequent thermal treatment at 500 °C in a reducing atmosphere removes the ligand shell, while preserving ultrasmall particle size and establishing direct Ir–TiN electronic coupling. The optimized catalyst exhibits mixed Ir 0 /Ir x+ coordination with low charge-transfer resistance (R ct = 19.2 Ω), delivering a mass activity of 342 A g Ir −1 at 1.54 V in acidic electrolyte. In situ X-ray absorption spectroscopy reveals irreversible surface oxidation as the primary stability-limiting factor. This stepwise strategy provides a general framework for supported catalysts that maximize precious metal utilization via sub-nanometer dispersion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extrapolation of nuclear structure observables with artificial neural networks

Calculations of nuclei are often carried out in finite model spaces. Thus, finite-size corrections enter, and it is necessary to extrapolate the computed observables to infinite model spaces. In this work, we employ extrapolation methods based on artificial neural networks for observables such as the ground-state energy and the point-proton radius. We extrapolate results from no-core shell model and coupled-cluster calculations to very large model spaces and estimate uncertainties. Training the network on different data typically yields extrapolation results that cluster around distinct values. We show that a preprocessing of input data, and the inclusion of correlations among the input data, reduces the problem of multiple solutions and yields more stable extrapolated results and consistent uncertainty estimates. We perform extrapolations for ground-state energies and radii in 4 He, 6 Li, and 16 O, and compare the predictions from neural networks with results from infrared extrapolations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

State Preparation in Quantum Algorithms for Fragment-Based Quantum Chemistry

State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space–unitary coupled cluster (LAS–UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. Here, in this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C–C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗

Core excitations with excited state mean field and perturbation theory

In this work, we test the efficacy of excited state mean field theory and its excited-state-specific perturbation theory on the prediction of K-edge positions and x-ray peak separations. We find that the mean field theory is surprisingly accurate, even though it contains no accounting of differential electron correlation effects. In the perturbation theory, we test multiple core-valence separation schemes and find that, with the mean field theory already so accurate, electron-counting biases in one popular separation scheme become a dominant error when predicting K-edges. Happily, these appear to be relatively easy to correct for, leading to a perturbation theory for K-edge positions that is lower scaling and more accurate than coupled cluster theory and competitive in accuracy with recent high-accuracy results from restricted open-shell Kohn-Sham theory. For peak separations, our preliminary data show excited state mean field theory to be exceptionally accurate, but more extensive testing will be needed to see how it and its perturbation theory compare to coupled cluster peak separations more broadly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗