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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

Quantum dynamics of charge state in silicon field evaporation

The charge state of an ion field-evaporating from a silicon-atom cluster is analyzed using time-dependent density functional theory coupled to molecular dynamics. The final charge state of the ion is shown to increase gradually with increasing external electrostatic field in agreement with the average charge state of silicon ions detected experimentally. When field evaporation is triggered by laser-induced electronic excitations the charge state also increases with increasing intensity of the laser pulse. At the evaporation threshold, the charge state of the evaporating ion does not depend on the electrostatic field due to the strong contribution of laser excitations to the ionization process both at low and high laser energies. A neutral silicon atom escaping the cluster due to its high initial kinetic energy is shown to be eventually ionized by external electrostatic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Phase diagram of magnetic shape memory alloy Ni 50 Mn $50–x$ In $x$ , 0 < $x$ , 25 from first principles, via spin cluster expansion and phonon vibrational entropies

The metamagnetic shape memory Heusler alloy Ni 50 Mn $50–x$ In $x$ exhibits a rich phase diagram featuring competing magnetic states, coupled magnetic–structural phase transitions, and strong compositional sensitivity. Existing first-principles approaches struggletocapturetheintertwinedchemical, magnetic, andvibrationaleffectsinthesealloys, necessitating a more integrated modeling framework. We develop a spin cluster expansion (spin-CE) framework augmented by a quasi-harmonic phonon model to capture both configurational (chemical and magnetic) and vibrational contributions to the free energy of Ni 50 Mn $50–x$ In $x$ over the full compositional range 0 ≤ x ≤25. The spin-CE includes both chemical clusters and composition-dependent Ising spin interactions, with parameters fit to a first-principles density functional theory (DFT) dataset. Using this approach, we predict the complete magnetostructural phase diagram and transformation temperatures of Ni 50 Mn $50–x$ In $x$ across the composition space. We find that vibrational entropy alone is insufficient to reproduce the martensitic transformation in the magnetic shape memory alloy regime, highlighting the essential role of magnetism. Incorporating both magnetic and vibrational contributions allows us to reproduce all experimentally known phases, including the disappearance of the stable martensite phase at a critical In concentration and the Curie temperature of the austenite phase. The method also captures the transition with increasing In in martensite from antiferromagnetic to ferromagnetic order and predicts re-entrant ferromagnetism, though the latter occurs at higher In content than reported experimentally. We discuss possible sources of this discrepancy and highlight the broader applicability of the method to other magnetostructurally complex materials, where it may offer mechanistic insight and predictive design capabilities.

Cluster expansion↗
DOE OSTI2015

The hyperfine structure in the rotational spectra of D{sub 2}{sup 17}O and HD{sup 17}O: Confirmation of the absolute nuclear magnetic shielding scale for oxygen

Guided by theoretical predictions, the hyperfine structures of the rotational spectra of mono- and bideuterated-water containing {sup 17}O have been experimentally investigated. To reach sub-Doppler resolution, required to resolve the hyperfine structure due to deuterium quadrupole coupling as well as to spin-rotation (SR) and dipolar spin-spin couplings, the Lamb-dip technique has been employed. The experimental investigation and in particular, the spectral analysis have been supported by high-level quantum-chemical computations employing coupled-cluster techniques and, for the first time, a complete experimental determination of the hyperfine parameters involved was possible. The experimentally determined {sup 17}O spin-rotation constants of D{sub 2}{sup 17}O and HD{sup 17}O were used to derive the paramagnetic part of the corresponding nuclear magnetic shielding constants. Together with the computed diamagnetic contributions as well as the vibrational and temperature corrections, the latter constants have been employed to confirm the oxygen nuclear magnetic shielding scale, recently established on the basis of spin-rotation data for H{sub 2}{sup 17}O [Puzzarini et al., J. Chem. Phys. 131, 234304 (2009)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An Orthogonal Recursive Bisection (ORB) Based Time Advancement Algorithm for CFD-DEM Solvers

The time integration of the granular phase in coupled computational fluid dynamics (CFD) – discrete element method (DEM) simulations presents a unique computational challenge brought about by the large variations in particle collisional time scales. Particles in the dilute regions of the computational domain can be advanced with large time steps while dense regions require much smaller time increments. However, the time step size in most solvers is globally set as the limit for accuracy and stability imposed by the collisions and is typically orders of magnitude less than that required away from collisions. This work addresses this precise issue and provides a strategy to avoid the use of a global conservative small time step size for the entire set of particles.A novel time stepping algorithm for CFD-DEM solvers using a partitioning approach using orthogonal recursive bisection (ORB) that allows for variable time steps among particles is described and its computational performance is compared against baseline explicit methods, typically used in several CFD-DEM solvers. ORB has advantages of being relatively quick and easy to update incrementally and has the required heuristic behavior (i.e., it will split the region in half with a cluster on each side) when groups of particles are well separated (clustered). The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for subsets of particles that are present at the leaves of the ORB, thereby resulting in substantial reduction of computational cost. The parallel implementation of this method where a ``knapsack” algorithm is used in tandem with ORB for effective load-balancing is also presented, where a best possible partitioning is obtained based on number of particles and local time-stepping costs. The algorithm is tested against benchmark problems with varying particle distributions that include fluidized bed and riser flow scenarios. Preliminary results indicate that the approach is 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

adaptive timestepping↗
NASA NTRS1985

The missing bulge globular clusters in M31 - New optical candidates

A new method to attack the question of the 'missing' globular clusters in the bulge of M31 is used. Image-processing techniques were used on 13 videocamera fields to obtain an accurate photometric census of stellar objects in M31's bulge down to a limiting B magnitude of 21. This luminosity distribution is compared with the Bahcall-Soneira model of galactic foreground stars. A statistically significant excess of bright images in the luminosity range of globular clusters at M31's distance is found. If the optical candidates considered prove to be globular clusters, they would double the number of known globular clusters in the surveyed region. The colors of a subsample of the candidates are the same as those of the known globular clusters. It is concluded that the previously observed flattening away from a de Vaucouleurs law in the radial distribution of M31 may be an observational selection effect. As an offshoot of this analysis, no evidence is found for very luminous stars in the inner bulge of M31. The lack of such stars indicates that there has not been active star formation (with a normal IMF) in the recent past. Coupled with the existence of many planetary nebulae in the bulge, this may strengthen the case for a galactic wind in M31's bulge.

Wirth, A.↗
DOE OSTI2018

Performance tuning of CEED software and first wave apps

The goal of this milestone was the performance tuning of the CEED software and first wave apps. In this milestone, the CEED team developed optimization techniques and tuned for performance the CEED software. Specifically, the focus was on the following: (1) Fast finite element operator storage and evaluation using partial assembly/matrix-free algorithms that take advantage of tensor-product element structure; (2) Architecture optimizations, including the development of performant discretization libraries targeting heterogeneous systems with multi-core CPUs, GPUs, and/or many-core CPUs in single node, cluster, or large-scale ECP supercomputer configurations; and (3) Global kernels for finite element operators, including the development of algorithms and libraries for fast gather-scatter exchanges for action assembly on unstructured graphs and for components needed for global linear solvers. A main part of this milestone was the performance tuning of the CEED first-wave ECP applications. This included the ExaSMR application – Coupled Monte Carlo Neutronics and Fluid Flow Simulation of Small Modular Reactors (ORNL), and the MARBL application – Next-Gen Multi-physics Simulation Code (LLNL). The artifacts delivered include performance improvements in CEED’s 1st wave of applications, and tuned CEED software for various architectures through a number of backends, freely available in the CEED’s repository on GitHub. See the CEED website, http://ceed.exascaleproject.org and the CEED GitHub organization, http://github.com/ceed for more details. In addition to details and results from the above R&D efforts, in this document we are also reporting on other project-wide activities performed in Q4 of FY18 including: CEED’s second annual meeting, a successful minisymposium at the premier international conference on high-order methods (ICOSAHOM), making the CEED milestone reports publicly available, new Laghos release, and other outreach efforts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

SOPPA and CCSD vibrational corrections to NMR indirect spin-spin coupling constants of small hydrocarbons

We present zero-point vibrational corrections to the indirect nuclear spin-spin coupling constants in ethyne, ethene, cyclopropene and allene. The calculations have been carried out both at the level of the second order polarization propagator approximation (SOPPA) employing a new implementation in the DALTON program, at the density functional theory level with the B3LYP functional employing also the Dalton program and at the level of coupled cluster singles and doubles (CCSD) theory employing the implementation in the CFOUR program. Specialized coupling constant basis sets, aug-cc-pVTZ-J, have been employed in the calculations. We find that on average the SOPPA results for both the equilibrium geometry values and the zero-point vibrational corrections are in better agreement with the CCSD results than the corresponding B3LYP results. Furthermore we observed that the vibrational corrections are in the order of 5 Hz for the one-bond carbon-hydrogen couplings and about 1 Hz or smaller for the other couplings apart from the one-bond carbon-carbon coupling (11 Hz) and the two-bond carbon-hydrogen coupling (4 Hz) in ethyne. However, not for all couplings lead the inclusion of zero-point vibrational corrections to better agreement with experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multiscale Physics of Atomic Nuclei from First Principles

Atomic nuclei exhibit multiple energy scales ranging from hundreds of MeV in binding energies to fractions of an MeV for low-lying collective excitations. As the limits of nuclear binding are approached near the neutron and proton drip lines, traditional shell structure starts to melt with an onset of deformation and an emergence of coexisting shapes. It is a long-standing challenge to describe this multiscale physics starting from nuclear forces with roots in quantum chromodynamics. Here, we achieve this within a unified and nonperturbative quantum many-body framework that captures both short- and long-range correlations starting from modern nucleon-nucleon and three-nucleon forces from chiral effective field theory. The short-range (dynamic) correlations which account for the bulk of the binding energy are included within a symmetry-breaking framework, while long-range (static) correlations (and fine details about the collective structure) are included by employing symmetry projection techniques. Our calculations accurately reproduce—within theoretical error bars—available experimental data for low-lying collective states and the electromagnetic quadrupole transitions in 20−30 Ne. In addition, we reveal coexisting spherical and deformed shapes in 30 Ne, which indicates the breakdown of the magic neutron number 𝑁 = 20 as the key nucleus 28 O is approached, and we predict that the drip line nuclei 32,34 Ne are strongly deformed and collective. By developing reduced-order models for symmetry-projected states, we perform a global sensitivity analysis and find that the subleading singlet 𝑆-wave contact and a pion-nucleon coupling strongly impact nuclear deformation in chiral effective field theory. The techniques developed in this work clarify how microscopic nuclear forces generate the multiscale physics of nuclei spanning collective phenomena as well as short-range correlations and allow one to capture emergent and dynamical phenomena in finite fermion systems such as atom clusters, molecules, and atomic nuclei.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Statistical Tests for Diagnosing Fission Source Convergence and Undersampling in Monte Carlo Criticality Calculations [Slides]

There is a very strong need for statistical testing to determine fission source convergence in Monte Carlo criticality calculations. Automation of such tests will greatly streamline and support the work carried out by NCS practitioners. Recent R&D work has shown that no single statistical test for convergence is sufficiently reliable, robust, and “guaranteed.” However, a combination of several standard statistical tests for the similarity of distributions, coupled with a high-fidelity estimate of the fission-matrix source is sufficiently robust, reliable, and repeatable that convergence can be “guaranteed.” During the course of the EG-AMCT studies, a number of statistical metrics and tests were proposed for diagnosing clustering and undersampling. None of these was robust and reliable enough for practical use in production codes. However, the expert group efforts came close. Some recent R&D work stemming from those past efforts has been very successful and promising. This Sub-Group will provide international input and collaboration on the development and implementation of statistical tests for convergence, with the primary goal of having the MC codes automatically detect convergence (or lack thereof). Newly proposed statistical tests to detect undersampling (after convergence) will also be reviewed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Driven similarity renormalization group: Third-order multireference perturbation theory

Here, a third-order multireference perturbation theory based on the driven similarity renormalization group (DSRG-MRPT3) approach is presented. The DSRG-MRPT3 method has several appealing features: (a) it is intruder free, (b) it is size consistent, (c) it leads to a non-iterative algorithm with O(N 6 ) scaling, and (d) it includes reference relaxation effects. The DSRG-MRPT3 scheme is benchmarked on the potential energy curves of F 2 , H 2 O 2 , C 2 H 6 , and N 2 along the F–F, O–O, C–C, and N–N bond dissociation coordinates, respectively. The nonparallelism errors of DSRG-MRPT3 are consistent with those of complete active space third-order perturbation theory and multireference configuration interaction with singles and doubles and show significant improvements over those obtained from DSRG second-order multireference perturbation theory. Our efficient implementation of the DSRG-MRPT3 based on factorized electron repulsion integrals enables studies of medium-sized open-shell organic compounds. This point is demonstrated with computations of the singlet-triplet splitting (Δ ST = E T –E S ) of 9,10-anthracyne. At the DSRG-MRPT3 level of theory, our best estimate of the adiabatic Δ ST is 3.9 kcal mol –1 , a value that is within 0.1 kcal mol –1 from multireference coupled cluster results.

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

Vibronic coupling in asymmetric bichromophores: Theory and application to diphenylmethane-d{sub 5}

A theoretical model based on Fulton and Gouterman dimer Hamiltonian [J. Chem. Phys. 35, 1059 (1961)] is used to understand the jet cooled spectra of partly deuterated diphenylmethane (DPM-d{sub 5}), reported in adjoining paper by Zwier and co-workers [“Vibronic coupling in asymmetric bichromophores: Experimental investigation of diphenylmethane-d5,” J. Chem. Phys. 141, 064316 (2014)]. The model is extended to include multiple vibrational modes, vibrational asymmetries between the monomers, and inter-monomer modes. Time dependent density functional theory and equation of motion coupled cluster electronic structure calculations are used to obtain parameters for the model Hamiltonian. The monomer units for DPM-d{sub 5} are toluene and partially deuterated toluene-d{sub 5}. The comparisons of modeled and experimental spectra show that the model captures the vibronic effects on a multitude of observables, including anomalously intense vibrational states, transition dipole moment directions, and localization/delocalization phenomena of both the vibronic and electronic wavefunctions. Following the success of the model, a qualitative picture of vibronic coupling in DPM-d{sub 5} is built using only three active vibrations. This illustrates the ability of the model to not just qualitatively predict but also phenomenologically explain the vibronic effects observed in spectra of complex bichromophores.

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

Formation and detection of metastable formic acid in a supersonic expansion: High resolution infrared spectroscopy of the jet-cooled cis -HCOOH conformer

High-resolution direct absorption infrared spectra of metastable cis-formic acid (HCOOH) trapped in a cis-well resonance behind a 15 kcal/mol barrier are reported for the first time, with the energetically unstable conformer produced in a supersonic slit plasma expansion of trans-formic acid/H 2 mixtures. We present a detailed high-resolution rovibrational analysis for cis-formic acid species in the OH stretch (ν 1 ) fundamental, providing first precision vibrational band origin, rotational constants, and term values, which in conjunction with ab initio calculations at the couple-cluster with single, double, and perturbative triple [CCSD(T)]/ANOn (n = 0, 1, 2) level support the experimental assignments and establish critical points on the potential energy surface for internal rotor trans-to-cis isomerization. Relative intensities for a- and b-type transitions observed in the spectra permit the transition dipole moment components to be determined in the body fixed frame and prove to be in good agreement with ab initio CCSD(T) theoretical estimates but in poor agreement with simple bond-dipole predictions. As a result, the observed signal dependence on H 2 in the discharge suggests the presence of a novel H atom radical chemical mechanism for strongly endothermic “up-hill” internal rotor isomerization between trans- and cis-formic acid conformers.

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

Nonadiabatic Hydrogen Dissociation on Copper Nanoclusters

Copper surfaces display high catalytic selectivity but have poor hydrogen dissociation kinetics; therefore, we consider icosahedral Cu 13 nanoclusters to understand how nanoscale structure might improve catalytic prospects. We find that the spin state is a surprisingly important design consideration. Cu 13 clusters have large magnetic moments due to finite size and symmetry effects and exhibit magnetization-dependent catalytic behavior. The most favorable transition state for hydrogen dissociation has a lower activation energy than that on single-crystal copper surfaces but requires a magnetization switch from 5 to 3 μ B . Without this switch, the activation energy is higher than that on single-crystal surfaces. Weak spin–orbit coupling hinders this switch, decreasing the kinetic rate of hydrogen dissociation by a factor of 16. Here, we consider strategies to facilitate magnetization switches through optical excitations, substitution, charge states, and co-catalysts; these considerations demonstrate how control of magnetic properties could improve catalytic performance.

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

Bond energies of ThO{sup +} and ThC{sup +}: A guided ion beam and quantum chemical investigation of the reactions of thorium cation with O{sub 2} and CO

Kinetic energy dependent reactions of Th{sup +} with O{sub 2} and CO are studied using a guided ion beam tandem mass spectrometer. The formation of ThO{sup +} in the reaction of Th{sup +} with O{sub 2} is observed to be exothermic and barrierless with a reaction efficiency at low energies of k/k{sub LGS} = 1.21 ± 0.24 similar to the efficiency observed in ion cyclotron resonance experiments. Formation of ThO{sup +} and ThC{sup +} in the reaction of Th{sup +} with CO is endothermic in both cases. The kinetic energy dependent cross sections for formation of these product ions were evaluated to determine 0 K bond dissociation energies (BDEs) of D{sub 0}(Th{sup +}–O) = 8.57 ± 0.14 eV and D{sub 0}(Th{sup +}–C) = 4.82 ± 0.29 eV. The present value of D{sub 0} (Th{sup +}–O) is within experimental uncertainty of previously reported experimental values, whereas this is the first report of D{sub 0} (Th{sup +}–C). Both BDEs are observed to be larger than those of their transition metal congeners, TiL{sup +}, ZrL{sup +}, and HfL{sup +} (L = O and C), believed to be a result of lanthanide contraction. Additionally, the reactions were explored by quantum chemical calculations, including a full Feller-Peterson-Dixon composite approach with correlation contributions up to coupled-cluster singles and doubles with iterative triples and quadruples (CCSDTQ) for ThC, ThC{sup +}, ThO, and ThO{sup +}, as well as more approximate CCSD with perturbative (triples) [CCSD(T)] calculations where a semi-empirical model was used to estimate spin-orbit energy contributions. Finally, the ThO{sup +} BDE is compared to other actinide (An) oxide cation BDEs and a simple model utilizing An{sup +} promotion energies to the reactive state is used to estimate AnO{sup +} and AnC{sup +} BDEs. For AnO{sup +}, this model yields predictions that are typically within experimental uncertainty and performs better than density functional theory calculations presented previously.

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

Role of Interactions and Correlations on Collective Dynamics of Molecular Motors Along Parallel Filaments

Cytoskeletal motors known as motor proteins are molecules that drive cellular transport along several parallel cytoskeletal filaments and support many biological processes. Experimental evidences suggest that they interact with the nearest molecules of their filament while performing any mechanical work. These interactions modify the microscopic level properties of motor proteins. In this work, a new version of two-channel totally asymmetric simple exclusion process, that incorporates the intra-channel interactions in a thermodynamically consistent way, is proposed. As the existing approaches for multi-channel systems deviate from analyzing the combined effect of inter and intra-channel interactions, a new approach known as modified vertical cluster mean field is developed. The approach along with Monte Carlo simulations successfully encounters some correlations and computes the complex dynamic properties of the system. Role of symmetry of interactions and inter-channel coupling is observed on the phase diagrams, maximal particle current and its corresponding optimal interaction strength. Surprisingly, for all values of coupling rate and most of the interaction splittings, the optimal interaction strength corresponding to maximal current belongs to the case of weak repulsive interactions. Moreover, for weak interaction splittings and with an increase in the coupling rate, the optimal interaction strength tends towards the known experimental results. The effect of coupling as well as interaction energy is also measured for correlations. They are found to be short-range and weaker for repulsive and weak attractive interactions while they are long-range and stronger for large attractions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Cosmic variance and the inhomogeneous UV luminosity function of galaxies during reionization

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.

79 ASTRONOMY AND ASTROPHYSICS↗