Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computers orbits”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

The spin and orbital contributions to the total magnetic moments of free Fe, Co, and Ni clusters

We present size dependent spin and orbital magnetic moments of cobalt (Co{sub n}{sup +}, 8 ≤ n ≤ 22), iron (Fe{sub n}{sup +}, 7 ≤ n ≤ 17), and nickel cluster (Ni{sub n}{sup +}, 7 ≤ n ≤ 17) cations as obtained by X-ray magnetic circular dichroism (XMCD) spectroscopy of isolated clusters in the gas phase. The spin and orbital magnetic moments range between the corresponding atomic and bulk values in all three cases. We compare our findings to previous XMCD data, Stern-Gerlach data, and computational results. We discuss the application of scaling laws to the size dependent evolution of the spin and orbital magnetic moments per atom in the clusters. We find a spin scaling law “per cluster diameter,” ∼n{sup −1/3}, that interpolates between known atomic and bulk values. In remarkable contrast, the orbital moments do likewise only if the atomic asymptote is exempt. A concept of “primary” and “secondary” (induced) orbital moments is invoked for interpretation.

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

Merger of a Neutron Star with a Newtonian Black Hole

Newtonian smooth particle hydro simulations are presented of the merger of a 1.4 solar mass neutron star with a black hole of equal mass. The initial state of the system is modeled with a stiff polytrope orbiting a point mass. Dynamical instability sets in when the orbital separation is equal to about three stellar radii. The ensuing mass transfer occurs on the dynamical timescale. No accretion torus is formed. At the end of the computation a corona of large extent shrouds an apparently stable binary system of a 0.25 solar mass star orbiting a 2.3 solar mass black hole.

Lee, William H.↗
DOE OSTI2022

3D ion gyro-orbit heat load predictions for NSTX-U

Abstract High power tokamaks operate with divertor heat loads capable of destroying the plasma facing components (PFCs). High fidelity heat load predictions are necessary to ascertain the PFC state for design and during operation. Typical heat flux calculations are 2D, time invariant, and assume that power flows directly along the magnetic field lines (the optical approximation). These assumptions neglect the complex 3D geometries employed to protect the PFCs, the time varying nature of the plasma and PFC thermal state, and the helical trajectories of ions with finite Larmor radii (the gyro-orbit approximation). An integrated software framework, the heat flux engineering analysis toolkit (HEAT), was developed to generate time varying optical heat loads applied to real engineering computer aided design (CAD) (Looby et al 2022 Fusion Sci. Technol. 78 10–27). Recently, an ion-gyro orbit module has been added to HEAT. This module calculates the helical trajectories of ions as they gyrate about the magnetic field lines using kinetic theory macro-particles to accelerate the calculation. First, the new gyro-orbit module will be presented. Next, a comparison to existing research is performed. Finally, an analysis of the gyro-orbit heat loads for NSTX-U is presented for diverted discharges using the engineering CAD models utilized for PFC fabrication. Including these gyro-orbit effects can enhance the PFC performance by ‘smearing’ out the magnetic shadows associated with the castellated fish-scaled geometry. Simultaneously, the helical trajectories can degrade performance when they load narrow regions on edges and corners with high heat fluxes. Analysis of the trade-offs between these competing effects is included, and regions for further investigation are identified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Facing the Challenges of Borderline Oxidation State Assignments Using State-of-the-Art Computational Methods.

The oxidation state (OS) of metals and ligands in inorganic complexes may be defined by carefully curated rules, such as from IUPAC, or by computational procedures such as the effective oxidation state (EOS) or localized orbital bonding analysis (LOBA). Such definitions typically agree for systems with simple ionic bonding and innocent ligands but may disagree as the boundary between ionic and covalent bonds is approached, or as the role of ligand noninnocence becomes nontrivial, or high oxidation states of metals are supported by heavy dative bonding, and so on. This work systematically compares IUPAC, EOS, and LOBA across a series of complexes where OS assignment is challenging. These systems include high-valent transition metal oxides, transition metal complexes with noninnocent ligands such as dithiolate and nitrosyl, metal sulfur dioxide adducts, and two transition metal carbene complexes. The differences in OS assignment by the three methods are carefully discussed, demonstrating the synergy between EOS and LOBA. In addition, a clarity index for LOBA OS assignments is introduced that provides an indication of whether or not its predictions are close to the ionic-covalent boundary.

Gimferrer, Martí↗
DOE OSTI2026

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Spin–Orbit Coupling via Four-Component Multireference Methods: Benchmarking on p-Block Elements and Tentative Recommendations

Within current electronic structure theory methods, fully relativistic four-component (4c) approaches based on the Dirac Hamiltonian treat spin–orbit coupling with the most rigor. The spin treatment arises naturally from the formulation and does not need to be included ad hoc. Spin–orbit splittings can provide insightful benchmark criteria for the assessment of 4c methods; however, there have not been extensive studies in this respect. Spin–orbit splittings of the p-block elements B–I were computed using the 4c-CASSCF, 4c-CASPT2, and 4c-MR-CISD+Q methods, as recently implemented in BAGEL, with uncontracted Dunning basis sets. Comparison with experiment reveals that the four-component methods yield good results, with most of the computed splittings falling within 15% of the experimental values. A large basis set is needed to obtain accurate splittings of the light elements B–F, while splittings of heavier elements show little basis dependence. The 4c-MR-CISD+Q method gave the best splittings for light elements, while 4c-CASSCF showed the best splittings for elements beyond fluorine. In conclusion, the 4c-CASPT2 method gave the best splittings for group 13 atoms.

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

Gas-Phase Complexes of Americium and Lanthanides with a Bis-triazinyl Pyridine: Reactivity and Bonding of Archetypes for F-Element Separations

Bis-triazinyl pyridines (BTPs) exhibit solution selectivity for trivalent americium over lanthanides (Ln), the origins of which remain uncertain. Here in this paper, electrospray ionization was used to generate gas-phase complexes [ML 3 ] 3+ , where M = La, Lu, or Am and L is EtBTP 2,6-bis(5,6-diethyl-1,2,4-triazin-3-yl)-pyridine. Collision-induced dissociation (CID) of [ML 3 ] 3+ in the presence of H 2 O yielded a protonated ligand [L(H)] + and hydroxide [ML 2 (OH)] 2+ or hydrate [ML(L–H)(H 2 O)] 2+ , where (L–H) - is a deprotonated ligand. Although solution affinities indicate stronger binding of BTPs toward Am 3+ versus Ln 3+ , the observed CID process is contrastingly more facile for M = Am versus Ln. To understand the disparity, density functional theory was employed to compute potential energy surfaces for two possible CID processes, for M = La and Am. In accordance with the CID results, both the rate determining transition state barrier and the net energy are lower for [AmL 3 ] 3+ versus [LaL 3 ] 3+ and for both product isomers, [ML 2 (OH)] 2+ and [ML(L–H)(H 2 O)] 2+ . More facile removal of a ligand from [AmL 3 ] 3+ by CID does not necessarily contradict stronger Am 3+ –L binding, as inferred from solution behavior. In particular, the formation of new bonds in the products can distort kinetics and thermodynamics expected for simple bond cleavage reactions. In addition to correctly predicting the seemingly anomalous CID behavior, the computational results indicate greater participation of Am 5f versus La 4f orbitals in metal–ligand bonding.

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

Near and Above Ionization Electronic Excitations with Non-Hermitian Real-Time Time-Dependent Density Functional Theory

We present a real-time time-dependent density functional theory (RT-TDDFT) prescription for capturing near and post-ionization excitations based on non-Hermitian von Neumann density matrix propagation with atom-centered basis sets, tuned range-separated DFT, and a phenomenological imaginary molecular orbital-based absorbing potential to mimic coupling to the continuum. The computed extreme ultraviolet absorption spectra for acetylene (C2H2), water (H2O), and Freon 12 (CF2Cl2) agree well with electron energy loss spectroscopy (EELS) data over the range 0 to 50 eV. The absorbing potential removes spurious high energy finite basis artifacts, yielding correct bound to bound transitions, metastable (autoionizing) resonance states, and consistent overall absorption shapes.

absorption↗
DOE OSTI2014

Relativistic theory of nuclear spin-rotation tensor with kinetically balanced rotational London orbitals

Both kinetically balanced (KB) and kinetically unbalanced (KU) rotational London orbitals (RLO) are proposed to resolve the slow basis set convergence in relativistic calculations of nuclear spin-rotation (NSR) coupling tensors of molecules containing heavy elements [Y. Xiao and W. Liu, J. Chem. Phys. 138, 134104 (2013)]. While they perform rather similarly, the KB-RLO Ansatz is clearly preferred as it ensures the correct nonrelativistic limit even with a finite basis. Moreover, it gives rise to the same “direct relativistic mapping” between nuclear magnetic resonance shielding and NSR coupling tensors as that without using the London orbitals [Y. Xiao, Y. Zhang, and W. Liu, J. Chem. Theory Comput. 10, 600 (2014)].

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

Analysis and comparison of CVS-ADC approaches up to third order for the calculation of core-excited states

The extended second order algebraic-diagrammatic construction (ADC(2)-x) scheme for the polarization operator in combination with core-valence separation (CVS) approximation is well known to be a powerful quantum chemical method for the calculation of core-excited states and the description of X-ray absorption spectra. For the first time, the implementation and results of the third order approach CVS-ADC(3) are reported. Therefore, the CVS approximation has been applied to the ADC(3) working equations and the resulting terms have been implemented efficiently in the adcman program. By treating the α and β spins separately from each other, the unrestricted variant CVS-UADC(3) for the treatment of open-shell systems has been implemented as well. The performance and accuracy of the CVS-ADC(3) method are demonstrated with respect to a set of small and middle-sized organic molecules. Therefore, the results obtained at the CVS-ADC(3) level are compared with CVS-ADC(2)-x values as well as experimental data by calculating complete basis set limits. The influence of basis sets is further investigated by employing a large set of different basis sets. Besides the accuracy of core-excitation energies and oscillator strengths, the importance of cartesian basis functions and the treatment of orbital relaxation effects are analyzed in this work as well as computational timings. It turns out that at the CVS-ADC(3) level, the results are not further improved compared to CVS-ADC(2)-x and experimental data, because the fortuitous error compensation inherent in the CVS-ADC(2)-x approach is broken. While CVS-ADC(3) overestimates the core excitation energies on average by 0.61% ± 0.31%, CVS-ADC(2)-x provides an averaged underestimation of −0.22% ± 0.12%. Eventually, the best agreement with experiments can be achieved using the CVS-ADC(2)-x method in combination with a diffuse cartesian basis set at least at the triple-ζ level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Tuned range separated hybrid functionals for solvated low bandgap oligomers

The description of charge transfer excitations has long been a challenge to time dependent density functional theory. The recently developed concept of “optimally tuned range separated hybrid (OT-RSH) functionals” has proven to describe charge transfer excitations accurately in many cases. However, describing solvated or embedded systems is yet a challenge. This challenge is not only computational but also conceptual, because the tuning requires identifying a specific orbital, typically the highest occupied one of the molecule under study. For solvated molecules, this orbital may be delocalized over the solvent. We here demonstrate that one way of overcoming this problem is to use a locally projected self-consistent field diagonalization on an absolutely localized molecular orbital expansion. We employ this approach to determine ionization energies and the optical gap of solvated oligothiophenes, i.e., paradigm low gap systems that are of relevance in organic electronics. Dioxane solvent molecules are explicitly represented in our calculations, and the ambiguities of straightforward parameter tuning in solution are elucidated. We show that a consistent estimate of the optimal range separated parameter (ω) at the limit of bulk solvation can be obtained by gradually extending the solvated system. In particular, ω is influenced by the solvent beyond the first coordination sphere. For determining ionization energies, a considerable number of solvent molecules on the first solvation shell must be taken into account. We demonstrate that accurately calculating optical gaps of solvated systems using OT-RSH can be done in three steps: (i) including the chemical environment when determining the range-separation parameter, (ii) taking into account the screening due to the solvent, and (iii) using realistic molecular geometries.

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

Multi-reference approach to the calculation of photoelectron spectra including spin-orbit coupling

X-ray photoelectron spectra provide a wealth of information on the electronic structure. The extraction of molecular details requires adequate theoretical methods, which in case of transition metal complexes has to account for effects due to the multi-configurational and spin-mixed nature of the many-electron wave function. Here, the restricted active space self-consistent field method including spin-orbit coupling is used to cope with this challenge and to calculate valence- and core-level photoelectron spectra. The intensities are estimated within the frameworks of the Dyson orbital formalism and the sudden approximation. Thereby, we utilize an efficient computational algorithm that is based on a biorthonormal basis transformation. The approach is applied to the valence photoionization of the gas phase water molecule and to the core ionization spectrum of the [Fe(H{sub 2}O){sub 6}]{sup 2+} complex. The results show good agreement with the experimental data obtained in this work, whereas the sudden approximation demonstrates distinct deviations from experiments.

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

Low-lying excited states by constrained DFT

Exploiting the machinery of Constrained Density Functional Theory (CDFT), we propose a variational method for calculating low-lying excited states of molecular systems. We dub this method eXcited CDFT (XCDFT). Excited states are obtained by self-consistently constraining a user-defined population of electrons, N c , in the virtual space of a reference set of occupied orbitals. By imposing this population to be N c = 1.0, we computed the first excited state of 15 molecules from a test set. Our results show that XCDFT achieves an accuracy in the predicted excitation energy only slightly worse than linear-response time-dependent DFT (TDDFT), but without incurring into problems of variational collapse typical of the more commonly adopted ΔSCF method. In addition, we selected a few challenging processes to test the limits of applicability of XCDFT. We find that in contrast to TDDFT, XCDFT is capable of reproducing energy surfaces featuring conical intersections (azobenzene and H3) with correct topology and correct overall energetics also away from the intersection. Venturing to condensed-phase systems, XCDFT reproduces the TDDFT solvatochromic shift of benzaldehyde when it is embedded by a cluster of water molecules. Finally, we find XCDFT to be a competitive method among single-reference methods for computations of excited states in terms of time to solution, rate of convergence, and accuracy of the result.

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

Survival of the most transferable at the top of Jacob’s ladder: Defining and testing the ω B97M(2) double hybrid density functional

A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resu|ting functional, ωB97M(2) is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 datq points clearly surpasses that of all of the tested density funstionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.

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

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

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

STARS DO NOT EAT THEIR YOUNG MIGRATING PLANETS: EMPIRICAL CONSTRAINTS ON PLANET MIGRATION HALTING MECHANISMS

The discovery of ''hot Jupiters'' very close to their parent stars confirmed that Jovian planets migrate inward via several potential mechanisms. We present empirical constraints on planet migration halting mechanisms. We compute model density functions of close-in exoplanets in the orbital semi-major axis-stellar mass plane to represent planet migration that is halted via several mechanisms, including the interior 1:2 resonance with the magnetospheric disk truncation radius, the interior 1:2 resonance with the dust sublimation radius, and several scenarios for tidal halting. The models differ in the predicted power-law dependence of the exoplanet orbital semi-major axis as a function of stellar mass, and thus we also include a power-law model with the exponent as a free parameter. We use a Bayesian analysis to assess the model success in reproducing empirical distributions of confirmed exoplanets and Kepler candidates that orbit interior to 0.1 AU. Our results confirm a correlation of the halting distance with stellar mass. Tidal halting provides the best fit to the empirical distribution of confirmed Jovian exoplanets at a statistically robust level, consistent with the Kozai mechanism and the spin-orbit misalignment of a substantial fraction of hot Jupiters. We can rule out migration halting at the interior 1:2 resonances with the magnetospheric disk truncation radius and the interior 1:2 resonance with the dust disk sublimation radius, a uniform random distribution, and a distribution with no dependence on stellar mass. Note that our results do not rule out Type-II migration, but rather eliminate the role of a circumstellar disk in stopping exoplanet migration. For Kepler candidates, which have a more restricted range in stellar mass compared to confirmed planets, we are unable to discern between the tidal dissipation and magnetospheric disk truncation braking mechanisms at a statistically significant level. The power-law model favors exponents in the range of 0.38-0.9. This is larger than that predicted for tidal halting (0.23-0.33), which suggests that additional physics may be missing in the tidal halting theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Direct probe of dark energy through gravitational lensing effect

We show that gravitational lensing can provide a direct method to probe the nature of dark energy at astrophysical scales. For lensing system as an isolated astrophysical object, we derive the dark energy contribution to gravitational potential as a repulsive power-law term, containing a generic equation of state parameter w . We find that it generates w -dependent and position-dependent modification to the conventional light orbital equation of w =−1. With post-Newtonian approximation, we compute its direct effect for an isolated lensing system at astrophysical scales and find that the dark energy force can deflect the path of incident light rays. We demonstrate that the dark-energy-induced deflection angle Δα{sub DE}∝ M {sup (1+1/3} {sup w} {sup )} (with 1+1/3 w > 0), which increases with the lensing mass M and consistently approaches zero in the limit M → 0. This effect is distinctive because dark energy tends to diffuse the rays and generates concave lensing effect . This is in contrast to the conventional convex lensing effect caused by both visible and dark matter. Measuring such concave lensing effect can directly probe the existence and nature of dark energy. We estimate this effect and show that the current gravitational lensing experiments are sensitive to the direct probe of dark energy at astrophysical scales. For the special case w =−1, our independent study favors the previous works that the cosmological constant can affect light bending, but our prediction qualitatively and quantitatively differ from the literature, including our consistent realization of Δα{sub DE} → 0 (under 0 M → ) at the leading order.

79 ASTRONOMY AND ASTROPHYSICS↗