Search NASA⌕ Search

SEARCH · Search NASA

Results for “ORBITAL CALCULATION”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A hybrid framework of first principles molecular orbital calculations and a three-dimensional integral equation theory for molecular liquids: Multi-center molecular Ornstein–Zernike self-consistent field approach

In this study, we reported the development of a new quantum mechanics/molecular mechanics (QM/MM)-type framework to describe chemical processes in solution by combining standard molecular-orbital calculations with a three-dimensional formalism of integral equation theory for molecular liquids (multi-center molecular Ornstein–Zernike (MC-MOZ) method). The theoretical procedure is very similar to the 3D-reference interaction site model self-consistent field (RISM-SCF) approach. Since the MC-MOZ method is highly parallelized for computation, the present approach has the potential to be one of the most efficient procedures to treat chemical processes in solution. Benchmark tests to check the validity of this approach were performed for two solute (solute water and formaldehyde) systems and a simple S{sub N}2 reaction (Cl{sup −} + CH{sub 3}Cl → ClCH{sub 3} + Cl{sup −}) in aqueous solution. The results for solute molecular properties and solvation structures obtained by the present approach were in reasonable agreement with those obtained by other hybrid frameworks and experiments. In particular, the results of the proposed approach are in excellent agreements with those of 3D-RISM-SCF.

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

A new generation of effective core potentials from correlated and spin–orbit calculations: Selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with 4d, 5d, 6s, and 6p valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin–orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and SO splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin–orbit effects. The subsequent quantum Monte Carlo calculations point out the importance of accurate trial wave functions that, in some cases (mid-series transition elements), require treatment well beyond a single-reference.

Wang, Guangming (ORCID:0000000231327925)↗
DOE OSTI2016

Improvements to an ion orbit loss calculation in the tokamak edge

An existing model of collisionless particle, momentum, and energy ion orbit loss from the edge region of a diverted tokamak plasma has been extended. The extended ion orbit loss calculation now treats losses of both thermal ions and fast neutral beam injection ions and includes realistic flux surface and magnetic field representations, particles returning to the plasma from the scrape off layer, and treatment of x-transport and x-loss. More realistic flux surface geometry allows the intrinsic rotation calculation to predict a peaking in the profile closer to the separatrix, which is consistent with experiment; and particle tracking calculations reveal a new mechanism of “x-transport pumping,” which predicts larger ion losses when coupling conventional ion orbit loss and x-loss mechanisms, though still dominated by conventional ion orbit loss. Finally, sensitivity to these ion orbit loss model enhancements is illustrated by fluid predictions of neoclassical rotation velocities and radial electric field profiles, with and without the enhancements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Composite Correlated Molecular Orbital Theory Calculations of Ring Strain for Use in Predicting Polymerization Reactions

Strained ring systems play an important role in synthesis and can be characterized by the ring strain energy (RSE). The RSE of 3, 4, 5, and 6 membered saturated and unsaturated ring systems containing N, O, P, and S heteroatoms and H, F, SiMe 3 , and SO 2 Me substituents were calculated at the G3(MP2) composite correlated molecular orbital theory level using up to 5 models to predict the RSE. Generally, the RSE decreased as ring size increased with a substantial decrease from 4 to 5 membered rings. Replacement of a ring CH 2 with P or S reduced the RSE, consistent with less angle strain. The RSE for unsaturated systems were generally greater than for saturated systems due to increased angle strain. No general trends were found with respect to substituent effects. In conclusion, the RSE values suggest that 3-pyrroline and 2-pyrroline and their derivatives may be able to support ring opening metathesis polymerization and warrant further study.

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

General framework for calculating spin–orbit couplings using spinless one-particle density matrices: Theory and application to the equation-of-motion coupled-cluster wave functions

Standard implementations of nonrelativistic excited-state calculations compute only one component of spin multiplets (i.e., Ms = 0 triplets); however, matrix elements for all components are necessary for deriving spin-dependent experimental observables. Wigner–Eckart’s theorem allows one to circumvent explicit calculations of all multiplet components. We generate all other spin–orbit matrix elements by applying Wigner–Eckart’s theorem to a reduced one-particle transition density matrix computed for a single multiplet component. In addition to computational efficiency, this approach also resolves the phase issue arising within Born–Oppenheimer’s separation of nuclear and electronic degrees of freedom. A general formalism and its application to the calculation of spin–orbit couplings using equation-of-motion coupled-cluster wave functions are presented. The two-electron contributions are included via the mean-field spin–orbit treatment. Intrinsic issues of constructing spin–orbit mean-field operators for open-shell references are discussed, and a resolution is proposed. The method is benchmarked by using several radicals and diradicals. The merits of the approach are illustrated by a calculation of the barrier for spin inversion in a high-spin tris(pyrrolylmethyl)amine Fe(II) complex.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Isolation and Study of Ruthenium–Cobalt Oxo Cubanes Bearing a High-Valent, Terminal Ru V –Oxo with Significant Oxyl Radical Character

High-valent Ru V -oxo intermediates have long been proposed in catalytic oxidation chemistry, but investigations into their electronic and chemical properties have been limited due to their reactive nature and rarity. The incorporation of Ru into the [Co 3 O 4 ] subcluster via the single-step assembly reaction of Co II (OAc) 2 (H 2 O) 4 (OAc = acetate), perruthenate (RuO 4 - ), and pyridine (py) yielded an unprecedented Ru(O)Co 3 (μ 3 -O) 4 (OAc) 4 (py) 3 cubane featuring an isolable, yet reactive, RuV-oxo moiety. EPR, ENDOR, and DFT studies reveal a valence-localized [Ru V (S = 1/2)Co III 3(S = 0)O 4 ] configuration and non-negligible covalency in the cubane core. Significant oxyl radical character in the Ru V -oxo unit is experimentally demonstrated by radical coupling reactions between the oxo cubane and both 2,4,6-tri-tert-butylphenoxyl and trityl radicals. The oxo cubane oxidizes organic substrates and, notably, reacts with water to form an isolable μ-oxo bis-cubane complex [(py) 3 (OAc) 4 Co 3 (μ 3 -O) 4 Ru]-O-[RuCo 3 (μ 3 -O) 4 (OAc) 4 (py) 3 ]. Redox activity of the Ru V -oxo fragment is easily tuned by the electron-donating ability of the distal pyridyl ligand set at the Co sites demonstrating strong electronic communication throughout the entire cubane cluster. Natural bond orbital calculations reveal cooperative orbital interactions of the [Co 3 O 4 ] unit in supporting the Ru V -oxo moiety via a strong π-electron donation.

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

Orbital-resolved calculations of two-center interferences in linear triatomic molecules

In this report we perform ab initio calculations of high-harmonic spectroscopy (HHS) of two-center interference phenomena in oriented carbon-dichalcogen molecules, using time-dependent density functional theory (TDDFT). We show that by resolving the harmonic response into contributions from individual Kohn-Sham orbitals, we can extract target-specific characteristics for both the spectral amplitude and phase. We also discuss that this extraction is predicated on a careful analysis and normalization of the harmonic spectrum. Finally, we present field-free scattering calculations that mimic the recollision step in high-order-harmonic generation and allow the calculation of recombination dipole matrix elements without explicitly calculating the scattering states of a molecule. We show that the orbital-resolved TDDFT HHS results and results based on our field-free scattering calculations are in very good agreement with each other.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Speckle Interferometry of Red Dwarf Stars

We report high-resolution optical speckle observations of 336 M dwarfs, which results in 113 measurements of the relative position of 80 systems and 256 other stars with no indications of duplicity. These are the first measurements for two of the systems. We also present the earliest measurements of relative position for 17 others. We include orbits for six of the systems, two revised and four reported for the first time. For one of the systems with a new orbit, G 161-7, we determine masses of 0.156 ± 0.011 and 0.1175±0.0079 M{sub ⊙} for the A and B components, respectively. All six of these new calculated orbits have short periods (between five and 38 years) and hold the promise of deriving accurate masses in the near future. For many other pairs we can establish their nature as physical or chance alignment, depending on their relative motion. Of the 80 systems, 32 have calculated orbits, 25 others are physical pairs, four are optical pairs, and 19 are currently unknown.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Method to integrate full particle orbit in toroidal plasmas

It is vital to integrate full particle orbit accurately when studying charged particle dynamics in electromagnetic waves with frequency higher than cyclotron frequency. We have derived a form of the Boris scheme using magnetic coordinates, which can be used effectively to integrate the cyclotron orbit in toroidal geometry over a long period of time. The proposed method has been verified by a full particle orbit simulation in toroidal geometry without high frequency waves. The full particle orbit calculation recovers guiding center banana orbit. This approach has better numeric properties than the conventional Runge-Kutta method for conserving particle energy and magnetic moment. The toroidal precession frequency is discovered to match that from guiding center simulation. Many other important phenomena in the presence of an electric field, such as E×B drift, Ware pinch effect and neoclassical polarization drift are also verified by the full orbit simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Guided ion beam and theoretical studies of the bond energy of SmS +

Previous work has shown that atomic samarium cations react with carbonyl sulfide to form SmS + + CO in an exothermic and barrierless process. To characterize this reaction further, the bond energy of SmS + is determined in the present study using guided ion beam tandem mass spectrometry. Reactions of SmS + with Xe, CO, and O 2 are examined. Results for collision-induced dissociation processes with all three molecules along with the endothermicity of the SmS + + CO → Sm + + COS exchange reaction are combined to yield D 0 (Sm + –S) = 3.37 ± 0.20 eV. The CO and O 2 reactions also yield a SmSO + product, with measured endothermicities that indicate D 0 (SSm + –O) = 3.73 ± 0.16 eV and D 0 (OSm + –S) = 1.38 ± 0.27 eV. The SmS + bond energy is compared with theoretical values characterized at several levels of theory, including CCSD(T) complete basis set extrapolations using all-electron basis sets. Multireference configuration interaction calculations with explicit spin-orbit calculations along with composite thermochemistry using the Feller-Peterson-Dixon method and all-electron basis sets were also explored for SmS+, and for comparison, SmO, SmO + , and EuO

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

The Separation and Hα Contrasts of Massive Accreting Planets in the Gaps of Transitional Disks: Predicted Hα Protoplanet Yields for Adaptive Optics Surveys

We present a massive accreting gap planet model that ensures large gaps in transitional disks are kept dust free by the scattering action of three coplanar quasi-circular planets in a 1:2:4 mean motion resonance (MMR). This model uses the constraint of the observed gap size, and the dust-free nature of the gap, to determine within ∼10% the possible orbits for three massive planets in an MMR. Calculated orbits are consistent with the observed orbits and Hα emission (the brightest line to observe these planets) for LkCa 15 b, PDS 70 b, and PDS 70 c within observational errors. Moreover, the model suggests that the scarcity of detected Hα planets is likely a selection effect of the current limitations of non-coronagraphic, low (<10%) Strehl, Hα imaging with adaptive optics (AO) systems used in past Hα surveys. We predict that as higher Strehl AO systems (with high-performance custom coronagraphs; like the 6.5 m Magellan Telescope MagAO-X system) are utilized at Hα, the number of detected gap planets will substantially increase by more than tenfold. For example, we show that >25 ± 5 new Hα “gap planets” are potentially discoverable by a survey of the best 19 transitional disks with MagAO-X. Detections of these accreting protoplanets will significantly improve our understanding of planet formation, planet growth and accretion, solar system architectures, and planet–disk interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Spin-orbit couplings within the equation-of-motion coupled-cluster framework: Theory, implementation, and benchmark calculations

We present a formalism and an implementation for calculating spin-orbit couplings (SOCs) within the EOM-CCSD (equation-of-motion coupled-cluster with single and double substitutions) approach. The following variants of EOM-CCSD are considered: EOM-CCSD for excitation energies (EOM-EE-CCSD), EOM-CCSD with spin-flip (EOM-SF-CCSD), EOM-CCSD for ionization potentials (EOM-IP-CCSD) and electron attachment (EOM-EA-CCSD). We employ a perturbative approach in which the SOCs are computed as matrix elements of the respective part of the Breit-Pauli Hamiltonian using zeroth-order non-relativistic wave functions. We follow the expectation-value approach rather than the response-theory formulation for property calculations. Both the full two-electron treatment and the mean-field approximation (a partial account of the two-electron contributions) have been implemented and benchmarked using several small molecules containing elements up to the fourth row of the periodic table. The benchmark results show the excellent performance of the perturbative treatment and the mean-field approximation. When used with an appropriate basis set, the errors with respect to experiment are below 5% for the considered examples. The findings regarding basis-set requirements are in agreement with previous studies. The impact of different correlation treatment in zeroth-order wave functions is analyzed. Overall, the EOM-IP-CCSD, EOM-EA-CCSD, EOM-EE-CCSD, and EOM-SF-CCSD wave functions yield SOCs that agree well with each other (and with the experimental values when available). Using an EOM-CCSD approach that provides a more balanced description of the target states yields more accurate results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Tidal radii and destruction rates of globular clusters in the Milky Way due to bulge-bar and disk shocking

We calculate orbits, tidal radii, and bulge-bar and disk shocking destruction rates for 63 globular clusters in our Galaxy. Orbits are integrated in both an axisymmetric and a nonaxisymmetric Galactic potential that includes a bar and a three-dimensional model for the spiral arms. With the use of a Monte Carlo scheme, we consider in our simulations observational uncertainties in the kinematical data of the clusters. In the analysis of destruction rates due to the bulge-bar, we consider the rigorous treatment of using the real Galactic cluster orbit instead of the usual linear trajectory employed in previous studies. We compare results in both treatments. We find that the theoretical tidal radius computed in the nonaxisymmetric Galactic potential compares better with the observed tidal radius than that obtained in the axisymmetric potential. In both Galactic potentials, bulge-shocking destruction rates computed with a linear trajectory of a cluster at its perigalacticons give a good approximation of the result obtained with the real trajectory of the cluster. Bulge-shocking destruction rates for clusters with perigalacticons in the inner Galactic region are smaller in the nonaxisymmetric potential than those in the axisymmetric potential. For the majority of clusters with high orbital eccentricities (e > 0.5), their total bulge+disk destruction rates are smaller in the nonaxisymmetric potential.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Chemical Bonding Governs Complex Magnetism in MnPt 5 P

Subtle changes in chemical bonds may result in dramatic revolutions in magnetic properties in solid-state materials. MnPt 5 P, a derivative of the rare-earth-free ferromagnetic MnPt 5 As, was discovered and is presented in this work. Specifically, MnPt 5 P was synthesized, and its crystal structure and chemical composition were characterized by X-ray diffraction as well as energy-dispersive X-ray spectroscopy. Accordingly, MnPt 5 P crystallizes in the layered tetragonal structure with the space group P4/mmm (No. 123), in which the face-shared Mn@Pt 12 polyhedral layers are separated by P layers. In contrast to the ferromagnetism observed in MnPt 5 As, the magnetic properties measurements on MnPt 5 P show antiferromagnetic ordering occurs at ~188 K with a strong magnetic anisotropy in and out of the ab-plane. Moreover, a spin-flop transition appears when a high magnetic field is applied. An A-type antiferromagnetic structure was obtained from the analysis of powder neutron diffraction (PND) patterns collected at 150 and 9 K. Calculated electronic structures imply that hybridization of Mn-3d and Pt-5d orbitals is critical for both the structural stability and observed magnetic properties. Semiempirical molecular orbitals calculations on both MnPt 5 P and MnPt 5 As indicate that the lack of 4p character on the P atoms at the highest occupied molecular orbital (HOMO) in MnPt 5 P may cause the different magnetic behavior in MnPt 5 P compared to MnPt 5 As. The discovery of MnPt 5 P, along with our previously reported MnPt5As, parametrizes the end points of a tunable system to study the chemical bonding which tunes the magnetic ordering from ferromagnetism to antiferromagnetism with the strong spin–orbit coupling (SOC) effect.

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

An under-approximation of entropy for elemental multiconfigurational ground state electronic structures

A combinatorial approach has been applied to the allowable permutations of quantum electronic configurations under the constraints of Hund's rule for established ground state configurations toward an under-approximation of electronic structure entropy. Combined with a previously reported over-approximation, the approximations are used in conjunction in an attempt to bracket the upper and lower entropy limits for multiconfigurational ground state electronic structure entropy and compared to known standard molar entropies for the elements. This formality has been used for the application of a classical statistical mechanics methodology to be applied to the discrete sets of quantum mechanical states of Pu in order to calculate orbital occupancies in Pu's multiconfigurational ground state. Without consideration of the relative energies of various possible electronic configurations contributing to the multiconfigurational ground state, the calculations are performed under a general energy degeneracy assumption weighted to the number of permutations for specific configurations. The number of configurations assumed to significantly contribute is gradually constrained in order to approach a low-order approximation of orbital occupancies in Pu that are then compared to experimental and other calculated results from the literature.

Beaux, II, Miles F. (ORCID:000000032192626X)↗