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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Parallel multiscale simulations of a brain aneurysm

Cardiovascular pathologies, such as a brain aneurysm, are affected by the global blood circulation as well as by the local microrheology. Hence, developing computational models for such cases requires the coupling of disparate spatial and temporal scales often governed by diverse mathematical descriptions, e.g., by partial differential equations (continuum) and ordinary differential equations for discrete particles (atomistic). However, interfacing atomistic-based with continuum-based domain discretizations is a challenging problem that requires both mathematical and computational advances. We present here a hybrid methodology that enabled us to perform the first multiscale simulations of platelet depositions on the wall of a brain aneurysm. The large scale flow features in the intracranial network are accurately resolved by using the high-order spectral element Navier–Stokes solver NεκTαr. The blood rheology inside the aneurysm is modeled using a coarse-grained stochastic molecular dynamics approach (the dissipative particle dynamics method) implemented in the parallel code LAMMPS. The continuum and atomistic domains overlap with interface conditions provided by effective forces computed adaptively to ensure continuity of states across the interface boundary. A two-way interaction is allowed with the time-evolving boundary of the (deposited) platelet clusters tracked by an immersed boundary method. The corresponding heterogeneous solvers (NεκTαr and LAMMPS) are linked together by a computational multilevel message passing interface that facilitates modularity and high parallel efficiency. Results of multiscale simulations of clot formation inside the aneurysm in a patient-specific arterial tree are presented. We also discuss the computational challenges involved and present scalability results of our coupled solver on up to 300 K computer processors. Validation of such coupled atomistic-continuum models is a main open issue that has to be addressed in future work.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Unified ab initio approach to bound and unbound states: No-core shell model with continuum and its application to 7 He

In this study, we introduce a unified approach to nuclear bound and continuum states based on the coupling of the no-core shell model (NCSM), a bound-state technique, with the no-core shell model/resonating group method (NCSM/RGM), a nuclear scattering technique. This new ab initio method, no-core shell model with continuum (NCSMC), leads to convergence properties superior to either NCSM or NCSM/RGM while providing a balanced approach to different classes of states. In the NCSMC, the ansatz for the many-nucleon wave function includes (i) a square-integrable A-nucleon component expanded in a complete harmonic oscillator basis and (ii) a binary-cluster component with asymptotic boundary conditions that can properly describe weakly bound states, resonances, and scattering. The Schrödinger equation is transformed into a system of coupled-channel integral-differential equations that we solve using a modified microscopic R-matrix formalism within a Lagrange mesh basis. We demonstrate the usefulness of the approach by investigating the unbound 7 He nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Manipulating magnetic moments by superconducting currents

Here, we show that the interaction between a superconducting order parameter and the magnetic moment of an atomic cluster in a two-dimensional s-wave superconductor with Rashba spin-orbit coupling generates magnetic anisotropy that can be stronger or comparable to the magnetic anisotropy due to the crystal field and the shape of the cluster. Transport current through the superconductor produces the effective magnetic field acting on the cluster’s magnetic moment. The direction of the effective field depends on the direction of the current, thus allowing one to manipulate the magnetic moment by the superconducting current. Due to the large density of the superconducting current this method of magnetization reversal can be more advantageous at low temperatures than the spin-transfer torque method that requires a large spin-polarized current through a normal metal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Efficient Step-Merged Quantum Imaginary Time Evolution Algorithm for Quantum Chemistry

In this work, we develop a resource-efficient step-merged quantum imaginary time evolution approach (smQITE) to solve for the ground state of a Hamiltonian on quantum computers. This heuristic method features a fixed shallow quantum circuit depth along the state evolution path. We use this algorithm to determine the binding energy curves of a set of molecules, including H 2 , H 4 , H 6 , LiH, HF, H 2 O, and BeH 2 , and find highly accurate results. The required quantum resources of smQITE calculations can be further reduced by adopting the circuit form of the variational quantum eigensolver (VQE) technique, such as the unitary coupled cluster ansatz. We demonstrate that smQITE achieves a similar computational accuracy as VQE at the same fixed-circuit ansatz, without requiring a generally complicated high-dimensional nonconvex optimization. Finally, smQITE calculations are carried out on Rigetti quantum processing units, demonstrating that the approach is readily applicable on current noisy intermediate-scale quantum devices.

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

Structural and Electronic Effects on the Properties of Fe 2 (dobdc) upon Oxidation with N 2 O

Here, we report electronic, vibrational, and magnetic properties, together with their structural dependences, for the metal–organic framework Fe 2 (dobdc) (dobdc 4 – = 2,5-dioxido-1,4-benzenedicarboxylate) and its derivatives, Fe 2 (O) 2 (dobdc) and Fe 2 (OH) 2 (dobdc)—species arising in the previously proposed mechanism for the oxidation of ethane to ethanol using N 2 O as an oxidant. Magnetic susceptibility measurements reported for Fe 2 (dobdc) in an earlier study and reported in the current study for Fe II 0.26 [Fe III (OH)] 1.74 (dobdc)(DMF) 0.15 (THF) 0.22 , which is more simply referred to as Fe 2 (OH) 2 (dobdc), were used to confirm the computational results. Theory was also compared to experiment for infrared spectra and X-ray powder diffraction structures. Structural and magnetic properties were computed by using Kohn-Sham density functional theory both with periodic boundary conditions and with cluster models. In addition, we studied the effects of different treatments of the exchange interactions on the magnetic coupling parameters by comparing several approaches to the exchange-correlation functional: generalized gradient approximation (GGA), GGA with empirical Coulomb and exchange integrals for 3d electrons (the GGA+U method), nonseparable gradient approximation with empirical Coulomb and exchange integrals for 3d electrons (GAM+U), hybrid GGA, meta-GGA, and hybrid meta- GGA. We found the coupling between the metal centers to be ferromagnetic in the case of Fe 2 (dobdc) and antiferromagnetic in the cases of Fe 2 (O) 2 (dobdc) and Fe 2 (OH) 2 (dobdc). The shift in magnetic coupling behavior correlates with the changing electronic structure of the framework, which derives from both structural and electronic changes that occur upon metal oxidation and addition of the charge-balancing oxo and hydroxo ligands.

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

Investigating the Ground-State Rotamers of n-Propylperoxy Radical

The n-propylperoxy radical has been described as a molecule of critical importance to studies of low temperature combustion. Ab initio methods were used to study this three-carbon alkylperoxy radical, normal propylperoxy. Reliable CCSD(T)/ANO0 geometries were predicted for the molecule's five rotamers. For each rotamer, energetic predictions were made using basis sets as large as the cc-pV5Z in conjunction with coupled cluster levels of theory up to CCSDT(Q). Along with the extrapolations, corrections for relativistic effects, zero-point vibrational energies, and diagonal Born--Oppenheimer corrections were used to further refine energies. The results indicate that the lowest conformer is the gauche-gauche (GG) rotamer followed by the gauche-trans (0.12 kcal mol^-1 above GG), trans-gauche (0.44 kcal mol^-1), gauche'-gauche (0.47 kcal mol^-1), and trans-trans (0.57 kcal mol^-1). Fundamental vibrational frequencies were obtained using second-order vibrational perturbation theory (VPT2). This is the first time anharmonic frequencies have been computed for this system. The most intense IR features include all but one of the C-H stretches. The O-O fundamental (1063 cm^-1 for the GG structure) also has a significant IR intensity, 19.6 km mol^-1. The anharmonicity effects on the potential energy surface were also used to compute vibrationally averaged r_g,0 K bond lengths, accounting for zero-point vibrations present within the molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Here we developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

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

Integrated profiling of three dimensional cell culture models and 3D microscopy

Motivation: Our goal is to develop a screening platform for quantitative profiling of colony organizations in 3D cell culture models. The 3D cell culture models, which are also imaged in 3D, are functional assays that mimic the in vivo characteristics of the tissue architecture more faithfully than the 2D cultures. However, they also introduce significant computational challenges, with the main barriers being the effects of growth conditions, fixations and inherent complexities in segmentation that need to be resolved in the 3D volume. Results: A segmentation strategy has been developed to delineate each nucleus in a colony that overcomes (i) the effects of growth conditions, (ii) variations in chromatin distribution and (iii) ambiguities formed by perceptual boundaries from adjacent nuclei. The strategy uses a cascade of geometric filters that are insensitive to spatial nonuniformity and partitions a clump of nuclei based on the grouping of points of maximum curvature at the interface of two neighboring nuclei. These points of maximum curvature are clustered together based on their coplanarity and proximity to define dissecting planes that separate the touching nuclei. The proposed curvature-based partitioning method is validated with both synthetic and real data, and is shown to have a superior performance against previous techniques. Validation and sensitivity analysis are coupled with the experimental design that includes a non-transformed cell line and three tumorigenic cell lines, which covers a wide range of phenotypic diversity in breast cancer. Colony profiling, derived from nuclear segmentation, reveals distinct indices for the morphogenesis of each cell line.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1972

Automatic classification of soils and vegetation with ERTS-1 data

Preliminary results of a test of a computerized analysis method using ERTS 1 data are presented. The method consisted of a four-spectral-band supervised, maximum likelihood, Gaussian classifier with training statistics derived through a combination of clustering and manual methods. The multivariate analysis method leads to the assignment of each resolution element of the data to one of a preselected set of discrete classes. The data frame was an area over the Texas-Oklahoma border including Lake Texoma. The study suggests that multispectral scanner data coupled with machine processing shows promise for earth surface cover surveys. Futhermore, the processing time is short and consequently the costs are low; a full frame can be analyzed completely within 48 hours.

Landgrebe, D. A.↗
NASA NTRS2001

The Transition from Hydrogen Bonding to Ionization in (HCI)n(NH3)n and (HCI)n(H2O)n Clusters: Consequences for Anharmonic Vibrational Spectroscopy

Anharmonic vibrational frequencies and intensities are calculated for 1:1 and 2:2 (HCl)(sub n)(NH3)(sub n) and (HCl)(sub n)(H2O)(sub n) complexes, employing the correlation-corrected vibrational self-consistent field method with ab initio potential surfaces at the MP2/TZP computational level. In this method, the anharmonic coupling between all vibrational modes is included, which is found to be important for the systems studied. For the 4:4 (HCl)(sub n)(H2O)(sub n) complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Just as the (HCl)(sub n)(NH3)(sub n) structure switches from hydrogen-bonded to ionic for n=2, the (HCl)(sub n)(H2O)(sub n) switches to ionic structure for n=4. For (HCl)2(H2O)2, the lowest energy structure corresponds to the hydrogen-bonded form. However, configurations of the ionic form are separated from this minimum by a barrier of less than an O-H stretching quantum. This suggests the possibility of experiments on ionization dynamics using infrared excitation of the hydrogen-bonded form. The strong cooperative effects on the hydrogen bonding, and concomitant transition to ionic bonding, makes an accurate estimate of the large anharmonicity crucial for understanding the infrared spectra of these systems. The anharmonicity is typically of the order of several hundred wave numbers for the proton stretching motions involved in hydrogen or ionic bonding, and can also be quite large for the intramolecular modes. In addition, the large cooperative effects in the 2:2 and higher order (HCl(sub n)(H2O)(sub n) complexes may have interesting implications for solvation of hydrogen halides at ice surfaces.

Chaban, Galina M.↗
DOE OSTI2017

Coupled Cluster and Density Functional Studies of Atomic Fluorine Chemisorption on Coronene as Model Systems for Graphene Fluorination

Potential energy curves (PECs) of atomic fluorine adsorption on coronene as a model for graphene or nanocarbon surfaces have been computed. The PECs were obtained by scanning the fluorine atom distance to one of the center carbon atoms of the coronene molecule as model system in a “top” position from 4.0 to1.0 Å in intervals of 0.1 Å using a variety of quantum chemical methods. Various density functional theory (DFT) functionals, such as B3LYP, PBE, PBE0, CAM-B3LYP, and LC-ωPBE; approximate DFT methods such as several levels of the density-functional tight-binding (DFTB) method, as well as ab initio wave function theory methods such as MP2, CCSD, CCSD(T), and G2MS extrapolations, were used to evaluate energies for B3LYP/cc-pVDZ PEC geometries. G2MS is an approximation to the highly accurate CCSD(T)/cc-pVTZ level of theory in our work. We found that fluorine is chemically adsorbed on coronene with a binding energy of 22.9 and 23.3 kcal/mol at the B3LYP/cc-pVDZ and G2MS levels of theory, respectively, and 18.3 and 19.3 kcal/mol after counterpoise correction. Here, we found that pure DFT functionals and their DFTB approximations fail to predict the correct dissociation limit due to the DFT-inherent self-interaction error and various limitations in the DFTB approximation itself.

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

U.S. stock market interaction network as learned by the Boltzmann machine

Here, we study historical dynamics of joint equilibrium distribution of stock returns in the U.S. stock market using the Boltzmann distribution model being parametrized by external fields and pairwise couplings. Within Boltzmann learning framework for statistical inference, we analyze historical behavior of the parameters inferred using exact and approximate learning algorithms. Since the model and inference methods require use of binary variables, effect of this mapping of continuous returns to the discrete domain is studied. The presented results show that binarization preserves the correlation structure of the market. Properties of distributions of external fields and couplings as well as the market interaction network and industry sector clustering structure are studied for different historical dates and moving window sizes. We demonstrate that the observed positive heavy tail in distribution of couplings is related to the sparse clustering structure of the market. We also show that discrepancies between the model’s parameters might be used as a precursor of financial instabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Coupled cluster Green's function: Past, Present, and Future

Coupled cluster Green's function (CCGF) approach has drawn much attention in recent years for targeting the molecular and material electronic structure problems from a many-body perspective in a systematically improvable way. Here, we will present a brief review of the history of how the Green's function method evolved with the wavefunction, early and recent development of CCGF theory, and more recently scalable CCGF software development. We will highlight some of the recent applications of CCGF approach and propose some potential applications that would emerge in the near future.

Peng, Bo↗
DOE OSTI2014

Electron correlation in solids via density embedding theory

Density matrix embedding theory [G. Knizia and G. K.-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)] and density embedding theory [I. W. Bulik, G. E. Scuseria, and J. Dukelsky, Phys. Rev. B 89, 035140 (2014)] have recently been introduced for model lattice Hamiltonians and molecular systems. In the present work, the formalism is extended to the ab initio description of infinite systems. An appropriate definition of the impurity Hamiltonian for such systems is presented and demonstrated in cases of 1, 2, and 3 dimensions, using coupled cluster theory as the impurity solver. Additionally, we discuss the challenges related to disentanglement of fragment and bath states. The current approach yields results comparable to coupled cluster calculations of infinite systems even when using a single unit cell as the fragment. The theory is formulated in the basis of Wannier functions but it does not require separate localization of unoccupied bands. The embedding scheme presented here is a promising way of employing highly accurate electronic structure methods for extended systems at a fraction of their original computational cost.

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

Quasi-degenerate perturbation theory using matrix product states

In this work, we generalize the recently proposed matrix product state perturbation theory (MPSPT) for calculating energies of excited states using quasi-degenerate (QD) perturbation theory. Our formulation uses the Kirtman-Certain-Hirschfelder canonical Van Vleck perturbation theory, which gives Hermitian effective Hamiltonians at each order, and also allows one to make use of Wigner’s 2n + 1 rule. Further, our formulation satisfies Granovsky’s requirement of model space invariance which is important for obtaining smooth potential energy curves. Thus, when we use MPSPT with the Dyall Hamiltonian, we obtain a model space invariant version of quasi-degenerate n-electron valence state perturbation theory (NEVPT), a property that the usual formulation of QD-NEVPT2 based on a multipartitioning technique lacked. We use our method on the benchmark problems of bond breaking of LiF which shows ionic to covalent curve crossing and the twist around the double bond of ethylene where significant valence-Rydberg mixing occurs in the excited states. In accordance with our previous work, we find that multi-reference linearized coupled cluster theory is more accurate than other multi-reference theories of similar cost.

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

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Classical Optimal Control for Energy Minimization Based On Diffeomorphic Modulation under Observable-Response-Preserving Homotopy

In this work, we introduce the so-called “Classical Optimal Control Optimization” (COCO) method for global energy minimization based on the implementation of the diffeomorphic modulation under observable-response-preserving homotopy (DMORPH) gradient algorithm. A probe particle with time-dependent mass $m(t;β)$ and dipole $μ(r,t;β)$ is evolved classically on the potential energy surface $V(r)$ coupled to an electric field $E(t;β)$, as described by the time-dependent density of states represented on a grid, or otherwise as a linear combination of Gaussians generated by the k-means clustering algorithm. Control parameters $β$ defining $m(t;β), μ(r,t;β),$ and $E(t;β)$ are optimized by following the gradients of the energy with respect to $β,$ adapting them to steer the particle toward the global minimum energy configuration. We find that the resulting COCO algorithm is capable of resolving near-degenerate states separated by large energy barriers and successfully locates the global minima of golf potentials on flat and rugged surfaces, previously explored for testing quantum annealing methodologies and the quantum optimal control optimization (QuOCO) method. Preliminary results show successful energy minimization of multidimensional Lennard-Jones clusters. Beyond the analysis of energy minimization in the specific model systems investigated, we anticipate COCO should be valuable for solving minimization problems in general, including optimization of parameters in applications to machine learning and molecular structure determination.

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

Formation and Combustion of Unconfined Drop Clusters in Microgravity

Single-drop and droplet array studies have become common methods to isolate and investigate the effects of any of the complexities that enter into the drop combustion process. Microgravity environments are required to allow larger drops to be studied while minimizing or eliminating the confounding effects of buoyancy. Based on the results from current isolated drop, drop array, and spray studies funded through the Microgravity Science and Applications Division, it has become clear that even with the effects of buoyancy removed, the extrapolation of results from droplet array studies to spray flames is difficult. The problem occurs because even the simplest spray systems introduce complexities of multi-disperse drop sizes and drop-drop interactions, coupled with more complicated fluid dynamics. Not only do these features make the interpretation of experimental data difficult, they also make the problem very difficult to analyze computationally. Group combustion models, in which the interaction between droplets is treated on a statistical manner, have become a popular method to investigate the behavior of large numbers of interacting droplets, particularly through the work of Ryan et al. and Bellan and co-workers. While these models idealize the actual spray systems to a point where they can be treated computationally, the experimental analogy to these models is difficult to achieve because it requires the formation and Combustion of drop clusters without the effects of buoyancy. Therefore, even though these models have provided useful and insightful information, the verification of the results by direct comparison with experimental data is still lacking.

Liu, S.↗