Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled cluster method”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

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↗
DOE OSTI2016

Communication: Variation after response in quantum Monte Carlo

We present a new method for modeling electronically excited states that overcomes a key failing of linear response theory by allowing the underlying ground state ansatz to relax in the presence of an excitation. The method is variational, has a cost similar to ground state variational Monte Carlo, and admits both open and periodic boundary conditions. We present preliminary numerical results showing that, when paired with the Jastrow antisymmetric geminal power ansatz, the variation-after-response formalism delivers accuracies for valence and charge transfer single excitations on par with equation of motion coupled cluster, while surpassing coupled cluster’s accuracy for excitations with significant doubly excited character.

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

Dynamics of Water Molecules and Ions in Concentrated Lithium Chloride Solutions Probed with Ultrafast 2D IR Spectroscopy

Water and ion dynamics in concentrated LiCl solutions were studied using ultrafast 2D IR spectroscopy with the methyl thiocyanate (MeSCN) CN stretch as the vibrational probe. The IR absorption spectrum of MeSCN has two peaks, one peak for water associated with the nitrogen lone pair of MeSCN (W) and the other peak corresponding to Li + associated with the lone pair (L). To extract the spectral diffusion (structural dynamics) of W and L species, we developed a method that isolates the peak of interest by subtracting the 2D Gaussian proxies of multiple interfering peaks. Center line slope data (normalized frequency–frequency correlation function) for 2D bands from the W and L are fit with triexponential functions. The fastest component (1.1–1.6 ps) is assigned to local hydrogen bond length fluctuations. The intermediate timescale (~4.0 ps) corresponds to the hydrogen bond network rearrangement. The slowest component decays in ~40 ps and corresponds to ion pair and ion cluster dynamics. The very similar W and L spectral diffusion indicates that the motions of the water and ions are strongly coupled. Orientational relaxations of the W and L species were extracted using a new method to eliminate the effects of overlapping peaks. The results show that MeSCN bound to water undergoes orientational relaxation significantly faster than MeSCN bound to Li + . The orientational and spectral diffusion results are compared. In conclusion, a Stark coupling model is used to extract the root mean square average electric field caused by the ion clouds along the CN moiety as a function of concentration.

14 SOLAR ENERGY↗
DOE OSTI2020

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.↗
DOE OSTI2016

Multiscale modelling of evolving foams

We present a set of multi-scale interlinked algorithms to model the dynamics of evolving foams. These algorithms couple the key effects of macroscopic bubble rearrangement, thin film drainage, and membrane rupture. For each of the mechanisms, we construct consistent and accurate algorithms, and couple them together to work across the wide range of space and time scales that occur in foam dynamics. These algorithms include second order finite difference projection methods for computing incompressible fluid flow on the macroscale, second order finite element methods to solve thin film drainage equations in the lamellae and Plateau borders, multiphase Voronoi Implicit Interface Methods to track interconnected membrane boundaries and capture topological changes, and Lagrangian particle methods for conservative liquid redistribution during rearrangement and rupture. We derive a full set of numerical approximations that are coupled via interface jump conditions and flux boundary conditions, and show convergence for the individual mechanisms. We demonstrate our approach by computing a variety of foam dynamics, including coupled evolution of three-dimensional bubble clusters attached to an anchored membrane and collapse of a foam cluster.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

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↗
DOE OSTI2017

PLATO (Parallel Load Assignment Tool): A Parallel Workload Partitioner for Particle-based Methods

A key first step in many large-scale high-performance scientific computing applications is the distribution of the computational domain across all participating processors of a parallel computing platform such as a shared memory machine, cluster or a tightly-coupled supercomputer. The primary goal of such a partitioning step is to ensure that the computational work performed by each process remains load-balanced throughout the execution of the parallel application. In this project, a general-purpose parallel partitioning tool, code-named PLATO, was designed and implemented. PLATO assumes that the workload is a set of points in an n-dimensional metric space. These points are abstractions of computational units that may represent point-like computational units such as celestial objects in astrophysics, atoms in materials sciences and molecules in biology/chemistry, or may represent more abstract computational units such as grid points in finite difference methods, grid elements in finite element methods or n-dimensional points in the feature space of a host of other applications. PLATO supports several key and useful capabilities, namely: (a) multiple load partitioning strategies, (b) support for arbitrary user-defined processor topology, (c) pre-fetching data from neighboring processors based on user-defined parameters, and (d) support for arbitrary number and types of attributes associated with each point in the computational domain. PLATO is implemented in the C programming language using Message Passing Interface (MPI) for inter-processor communications and its parallel performance tested on up to 48 processors of a shared memory machine. This report includes a detailed performance study of PLATO using parallel speedup and iso-efficiency as the two main metrics of parallel performance. In addition, a user manual and DOXYGEN generated code documentation have been included for completeness.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

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↗
DOE OSTI2024

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↗
DOE OSTI2015

Theoretical investigation of HNgNH{sub 3}{sup +} ions (Ng = He, Ne, Ar, Kr, and Xe)

The equilibrium geometries, harmonic frequencies, and dissociation energies of HNgNH{sub 3}{sup +} ions (Ng = He, Ne, Ar, Kr, and Xe) were investigated using the following method: Becke-3-parameter-Lee-Yang-Parr (B3LYP), Boese-Matrin for Kinetics (BMK), second-order Møller-Plesset perturbation theory (MP2), and coupled-cluster with single and double excitations as well as perturbative inclusion of triples (CCSD(T)). The results indicate that HHeNH{sub 3}{sup +}, HArNH{sub 3}{sup +}, HKrNH{sub 3}{sup +}, and HXeNH{sub 3}{sup +} ions are metastable species that are protected from decomposition by high energy barriers, whereas the HNeNH{sub 3}{sup +} ion is unstable because of its relatively small energy barrier for decomposition. The bonding nature of noble-gas atoms in HNgNH{sub 3}{sup +} was also analyzed using the atoms in molecules approach, natural energy decomposition analysis, and natural bond orbital analysis.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

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↗
DOE OSTI2026

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↗
DOE OSTI2020

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↗
DOE OSTI2014

Thermal decomposition of 1,3,3-trinitroazetidine (TNAZ): A density functional theory and ab initio study

Density functional theory and ab initio methods are employed to investigate decomposition pathways of 1,3,3-trinitroazetidine initiated by unimolecular loss of NO{sub 2} or HONO. Geometry optimizations are performed using M06/cc-pVTZ and coupled-cluster (CC) theory with single, double, and perturbative triple excitations, CCSD(T), is used to calculate accurate single-point energies for those geometries. The CCSD(T)/cc-pVTZ energies for NO{sub 2} elimination by N–N and C–N bond fission are, including zero-point energy (ZPE) corrections, 43.21 kcal/mol and 50.46 kcal/mol, respectively. The decomposition initiated by trans-HONO elimination can occur by a concerted H-atom and nitramine NO{sub 2} group elimination or by a concerted H-atom and nitroalkyl NO{sub 2} group elimination via barriers (at the CCSD(T)/cc-pVTZ level with ZPE corrections) of 47.00 kcal/mol and 48.27 kcal/mol, respectively. Thus, at the CCSD(T)/cc-pVTZ level, the ordering of these four decomposition steps from energetically most favored to least favored is: NO{sub 2} elimination by N–N bond fission, HONO elimination involving the nitramine NO{sub 2} group, HONO elimination involving a nitroalkyl NO{sub 2} group, and finally NO{sub 2} elimination by C–N bond fission.

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

Hierarchical Inorganic Assemblies for Artificial Photosynthesis

Artificial photosynthesis is an attractive approach for renewable fuel generation because it offers the prospect of a technology suitable for deployment on highly abundant, non-arable land. Recent leaps forward in the development of efficient and durable light absorbers and catalysts for oxygen evolution and the growing attention to catalysts for carbon dioxide activation brings into focus the tasks of hierarchically integrating the components into assemblies for closing of the photosynthetic cycle. A particular challenge is the efficient coupling of the multi-electron processes of CO 2 reduction and H 2 O oxidation. Among the most important requirements for a complete integrated system are catalytic rates that match the solar flux, efficient charge transport between the various components, and scalability of the photosynthetic assembly on the unprecedented scale of terawatts in order to have impact on fuel consumption. To address these challenges, we have developed in this paper a heterogeneous inorganic materials approach with molecularly precise control of light absorption and charge transport pathways. Oxo-bridged heterobinuclear units with metal-to-metal charge-transfer transitions absorbing deep in the visible act as single photon, single charge transfer pumps for driving multi-electron catalysts. A photodeposition method has been introduced for the spatially directed assembly of nanoparticle catalysts for selective coupling to the donor or acceptor metal of the light absorber. For CO 2 reduction, a Cu oxide cluster is coupled to the Zr center of a ZrOCo light absorber, while coupling of an Ir nanoparticle catalyst for water oxidation to the Co donor affords closing of the photosynthetic cycle of CO 2 conversion by H 2 O to CO and O 2 . Optical, vibrational, and X-ray spectroscopy provide detailed structural knowledge of the polynuclear assemblies. Time resolved visible and rapid-scan FT-IR studies reveal charge transfer mechanisms and transient surface intermediates under photocatalytic conditions for guiding performance improvements. Separation of the water oxidation and carbon dioxide reduction half reactions by a membrane is essential for efficient photoreduction of CO 2 by H 2 O to liquid fuel products. A concept of a macroscale artificial photosystem consisting of arrays of Co oxide–silica core–shell nanotubes is introduced in which each tube operates as a complete, independent photosynthetic unit with built-in membrane separation. The ultrathin amorphous silica shell with embedded molecular wires functions as a proton conducting, molecule impermeable membrane. Photoelectrochemical and transient optical measurements confirm tight control of charge transport through the membrane by the orbital energetics of the wire molecules. Finally, hierarchical arrangement of the components is accomplished by a combination of photodeposition, controlled anchoring, and atomic layer deposition methods.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Oscillator strengths, first-order properties, and nuclear gradients for local ADC(2)

We describe theory and implementation of oscillator strengths, orbital-relaxed first-order properties, and nuclear gradients for the local algebraic diagrammatic construction scheme through second order. The formalism is derived via time-dependent linear response theory based on a second-order unitary coupled cluster model. The implementation presented here is a modification of our previously developed algorithms for Laplace transform based local time-dependent coupled cluster linear response (CC2LR); the local approximations thus are state specific and adaptive. The symmetry of the Jacobian leads to considerable simplifications relative to the local CC2LR method; as a result, a gradient evaluation is about four times less expensive. Test calculations show that in geometry optimizations, usually very similar geometries are obtained as with the local CC2LR method (provided that a second-order method is applicable). As an exemplary application, we performed geometry optimizations on the low-lying singlet states of chlorophyllide a.

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

Accurate and efficient open-source implementation of domain-based local pair natural orbital (DLPNO) coupled-cluster theory using a t1-transformed Hamiltonian

We present an efficient, open-source formulation for coupled-cluster theory through perturbative triples with domain-based local pair natural orbitals [DLPNO-CCSD(T)]. Similar to the implementation of the DLPNO-CCSD(T) method found in the ORCA package, the most expensive integral generation and contraction steps associated with the CCSD(T) method are linear-scaling. In this work, we show that the t1-transformed Hamiltonian allows for a less complex algorithm when evaluating the local CCSD(T) energy without compromising efficiency or accuracy. Our algorithm yields sub-kJ mol−1 deviations for relative energies when compared with canonical CCSD(T), with typical errors being on the order of 0.1 kcal mol−1, using our TightPNO parameters. We extensively tested and optimized our algorithm and parameters for non-covalent interactions, which have been the most difficult interaction to model for orbital (PNO)-based methods historically. To highlight the capabilities of our code, we tested it on large water clusters, as well as insulin (787 atoms).

Chemistry↗