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

Quantum Orbital Minimization Method for Excited States Calculation on a Quantum Computer

Herein we propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H 2 , LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A new method to compute periodic orbits in general symplectic maps

The search of high-order periodic orbits has been typically restricted to problems with symmetries that help to reduce the dimension of the search space. Well-known examples include reversible maps with symmetry lines. The present work proposes a new method to compute high-order periodic orbits in twist maps without the use of symmetries. The method is a combination of the parameterization method in Fourier space and a Newton–Gauss multiple shooting scheme. The parameterization method has been successfully used in the past to compute quasi-periodic invariant circles. However, this is the first time that this method is used in the context of periodic orbits. Numerical examples are presented showing the accuracy and efficiency of the proposed method. Furthermore, the method is also applied to verify the renormalization prediction of the residues’ convergence at criticality (extensively studied in reversible maps) in the relatively unexplored case of maps without symmetries.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

SPECTROSCOPIC ORBITS FOR 15 LATE-TYPE STARS

Spectroscopic orbital elements are determined for 15 stars with periods from 8 to 6528 days with six orbits computed for the first time. Improved astrometric orbits are computed for two stars and one new orbit is derived. Visual orbits were previously determined for four stars, four stars are members of multiple systems, and five stars have Hipparcos “G” designations or have been resolved by speckle interferometry. For the nine binaries with previous spectroscopic orbits, we determine improved or comparable elements. For HD 28271 and HD 200790, our spectroscopic results support the conclusions of previous authors that the large values of their mass functions and lack of detectable secondary spectrum argue for the secondary in each case being a pair of low-mass dwarfs. The orbits given here may be useful in combination with future interferometric and Gaia satellite observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A New Orbit for the Eclipsing Binary V577 Oph

Pulsating stars in eclipsing binary systems are unique objects for providing constraints on stellar models. To fully leverage the information available from the binary system, full orbital radial velocity curves must be obtained. We report 23 radial velocities for components of the eclipsing binary V577 Oph, whose primary star is a δ Sct variable. The velocities cover a nearly complete orbit and a time base of 20 years. We computed orbital elements for the binary and compared them to the ephemeris computed by Creevey et al. The comparison shows marginally different results. In particular, a change in the systemic velocity by −2 km s{sup −1} is suggested by our results. We compare this systemic velocity difference to that expected due to reflex motion of the binary in response to the third body in the system. The systemic velocity difference is consistent with reflex motion, given our mass determination for the eclipsing binary and the orbital parameters determined by Volkov and Volkova for the three-body orbit. We see no evidence for the third body in our spectra, but we do see strong interstellar Na D lines that are consistent in strength with the direction and expected distance of V577 Oph.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets

Natural orbitals are often used to achieve a more compact representation of correlated wave-functions. Using natural orbitals computed as eigenstates of the virtual–virtual block of the state density matrix instead of the canonical Hartree–Fock orbitals results in smaller errors when the same fraction of virtual space is frozen. This strategy, termed frozen natural orbital (FNO) approach, is effective in reducing the cost of regular coupled-cluster (CC) calculations and some multistate methods, such as EOM-IP-CC (equation-of-motion CC for ionization potentials). This contribution extends the FNO approach to the EOM-SF-CC ansatz (EOM-CC with spin-flip). In contrast to EOM-IP-CCSD, EOM-SF-CCSD relies on high-spin open-shell references. Using FNOs computed for an open-shell reference leads to an erratic behavior of the EOM-SF-CC energies and properties due to an inconsistent truncation of the α and β orbital spaces. A general solution to problems arising in the EOM-CC calculations utilizing open-shell references, termed OSFNO (open-shell FNO), is proposed. Furthermore, by means of singular value decomposition (SVD) of the overlap matrix of the α and β orbitals, the OSFNO algorithm identifies the corresponding orbitals and determines virtual orbitals corresponding to the singly occupied space. This is followed by SVD of the singlet part of the state density matrix in the remaining virtual orbital subspace. The so-computed FNOs preserve the spin purity of the open-shell orbital subspace to the extent allowed by the original reference, thus facilitating a safe truncation of the virtual space. The performance of OSFNO is benchmarked for selected diradicals and triradicals.

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

THE SUN WAS NOT BORN IN M67

Using the most recent proper-motion determination of the old, solar-metallicity, Galactic open cluster M67 in orbital computations in a non-axisymmetric model of the Milky Way, including a bar and three-dimensional spiral arms, we explore the possibility that the Sun once belonged to this cluster. We have performed Monte Carlo numerical simulations to generate the present-day orbital conditions of the Sun and M67, and all the parameters in the Galactic model. We compute 3.5 Multiplication-Sign 10{sup 5} pairs of orbits Sun-M67 looking for close encounters in the past with a minimum distance approach within the tidal radius of M67. In these encounters we find that the relative velocity between the Sun and M67 is larger than 20 km s{sup -1}. If the Sun had been ejected from M67 with this high velocity by means of a three-body encounter, this interaction would have either destroyed an initial circumstellar disk around the Sun or dispersed its already formed planets. We also find a very low probability, much lower than 10{sup -7}, that the Sun was ejected from M67 by an encounter of this cluster with a giant molecular cloud. This study also excludes the possibility that the Sun and M67 were born in the same molecular cloud. Our dynamical results convincingly demonstrate that M67 could not have been the birth cluster of our solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

INFRARED SPECTROSCOPY OF SYMBIOTIC STARS. X. ORBITS FOR THREE S-TYPE SYSTEMS: V1044 CENTAURI, HEN 3-1213, AND SS 73-96

Employing new infrared radial velocities, we have computed orbits of the cool giants in three southern S-type symbiotic systems. The orbit for V1044 Cen, an M5.5 giant, has a period of 985 days and a modest eccentricity of 0.16. Hen 3-1213 is a K4 giant, yellow symbiotic with an orbital period of 533 days and a similar eccentricity of 0.18. For the M2 giant SS 73-96 the orbital period is 828 days, and this system has a somewhat larger eccentricity of 0.26. Measurement of the H i Paschen δ emission lines, which may at least partially reflect the motion of the secondary in SS 73-96, results in a mass ratio of 2.4 for the M giant relative to the presumed white dwarf. The estimated orbital inclinations of V1044 Cen and Hen 3-1213 are low, about 40°. However, for SS 73-96 the predicted inclination is 90°, and so an ephemeris for eclipses of the secondary or the hot nebula surrounding it is provided. A search of the orbital velocity residuals of V1044 Cen and SS 73-96 for pulsation periods produced no realistic or convincing period for either star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics

Abstract The weight function W ( X ) is a diagnostic sensitivity to phase-space variables X that relates the measured signal C to the distribution function F ( X ) through the equation C = ∫ W ( X ) F ( X ) d X . In the present work, an algorithm to calculate W for a diagnostic that measures 3 MeV protons produced in d(d,p)t fusion reactions between a fast ion and a thermal deuteron is developed. The emitted protons escape the tokamak on curved orbits and are detected. These curved orbits constitute effective diagnostic ‘sightlines’. The presented algorithm accounts for the complications associated with these curved sightlines. An initial calculation of time-reversed proton orbits computes effective solid angles and sightlines for the relevant range of incident proton velocity vectors. These precomputed orbits are then used within the framework of FIDASIM (2020 Plasma Phys. Cont. Fusion 62 105008) to calculate the reactivity averaged over the thermal distribution of the ‘target’ deuterons and the probability that a fast ion of specified energy and pitch has a gyroangle that is consistent with the kinematic equations along each of the sightlines. Comparisons with analytic formulas and with independent calculations for the Mega Amp Spherical Tokamak 3 MeV proton diagnostic verify the algorithm.

Physics↗
DOE OSTI2020

An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Here we present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Proper Motions, Orbits, and Tidal Influences of Milky Way Dwarf Spheroidal Galaxies

We combine Gaia early data release 3 astrometry with accurate photometry and utilize a probabilistic mixture model to measure the systemic proper motion of 52 dwarf spheroidal (dSph) satellite galaxies of the Milky Way (MW). For the 46 dSphs with literature line-of-sight velocities we compute orbits in both a MW and a combined MW + Large Magellanic Cloud (LMC) potential and identify Car II, Car III, Hor I, Hyi I, Phx II, and Ret II as likely LMC satellites. 40% of our dSph sample has a >25% change in pericenter and/or apocenter with the MW + LMC potential. For these orbits, we use a Monte Carlo sample for the observational uncertainties for each dSph and the uncertainties in the MW and LMC potentials. We predict that Ant II, Boo III, Cra II, Gru II, and Tuc III should be tidally disrupting by comparing each dSph's average density relative to the MW density at its pericenter. dSphs with large ellipticity (CVn I, Her, Tuc V, UMa I, UMa II, UMi, Wil 1) show a preference for their orbital direction to align with their major axis even for dSphs with large pericenters. We compare the dSph radial orbital phase to subhalos in MW-like N-body simulations and infer that there is not an excess of satellites near their pericenter. With projections of future Gaia data releases, we find that dSph's orbital precision will be limited by uncertainties in the distance and/or MW potential rather than in proper motion precision. Finally, we provide our membership catalogs to enable community follow-up.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

CAS without SCF—Why to use CASCI and where to get the orbitals

The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree–Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Lastly, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.

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

Toward Full Configuration Interaction for Transition-Metal Complexes

In this work, an efficacious approximation to full configuration interaction (FCI) is adapted to calculate singlet-triplet gaps for transition-metal complexes. This strategy, incremental FCI (iFCI), uses a many-body expansion to systematically add correlation to a simple reference wave function and therefore achieves greatly reduced computational costs compared to FCI. iFCI through the 3-body expansion is demonstrated on four model transition-metal complexes involving the metals Zn, V, and Cu. Screening techniques to increase the computational efficiency of iFCI are proposed and tested, showing reduction in the number of 3-body terms by more than 90% with controlled errors. The largest complex treated by iFCI has 142 valence electrons, all of which are correlated among the full set of 444 active orbitals. Computed spin gaps approach experimental results for the four complexes, though room for improvement remains.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

RADIAL VELOCITIES OF GALACTIC O-TYPE STARS. II. SINGLE-LINED SPECTROSCOPIC BINARIES

We report on new radial velocity measurements of massive stars that are either suspected binaries or lacking prior observations. This is part of a survey to identify and characterize spectroscopic binaries among O-type stars with the goal of comparing the binary fraction of field and runaway stars with those in clusters and associations. We present orbits for HDE 308813, HD 152147, HD 164536, BD-16 Degree-Sign 4826, and HDE 229232, Galactic O-type stars exhibiting single-lined spectroscopic variation. By fitting model spectra to our observed spectra, we obtain estimates for effective temperature, surface gravity, and rotational velocity. We compute orbital periods and velocity semiamplitudes for each system and note the lack of photometric variation for any system. These binaries probably appear single-lined because the companions are faint and because their orbital Doppler shifts are small compared to the width of the rotationally broadened lines of the primary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Fabrication of DNA/Hydroxyapatite nanocomposites by simulated body fluid for gene delivery

The hydroxyapatite (HA) formation on the surface of DNA molecules in simulated body fluid (SBF) was examined. The osteoconductivity is estimated using SBF having ion concentrations approximately equal to those of human blood plasma. After immersion for 4 weeks in SBF at 36.5 °C, the HA crystallites possessing 1-14 micrometer in diameter grew on the surface of DNA molecules. The leaf flake-like and spherical shapes morphologies were observed through scanning electron microscopy analysis. Original peaks of both of DNA and HA were characterized by fourier transform infrared spectroscopy. The Ca/P ratio (1.1-1.5) in HA was estimated by energy dispersive X-ray analysis. After biomineralization, the calculated weight ratio of DNA/HA was 18/82 by thermogravimetry/differential thermal analysis. The molecular orbital computer simulation has been used to probe the interaction of DNA with two charge-balancing ions, CaOH{sup +} and CaH{sub 2}PO{sub 4}{sup +}. The adsorption enthalpy of the two ions on DNA having negative value was the evidence for the interface in mineralization of HA in SBF.

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

The ezSpectra suite: An easy‐to‐use toolkit for spectroscopy modeling

Abstract A molecule's spectrum encodes information about its structure and electronic properties. It is a unique fingerprint that can serve as a molecular ID. Quantum chemistry calculations provide key ingredients for interpreting spectra, but modeling the spectra rarely ends there; it requires additional steps that entail combined treatments of electronic and nuclear degrees of freedom and account for specifics of the experimental setup (light energy, polarization, averaging over molecular orientations, temperature, etc.). This Software Focus article describes the ezSpectra suite, which currently comprises two stand‐alone open‐source codes: ezFCF and ezDyson . ezFCF calculates Franck–Condon factors, which yield vibrational progressions for polyatomic molecules, within the double‐harmonic approximation. ezDyson calculates absolute cross‐sections for photodetachment/photoionization processes and photoelectron angular distributions using Dyson orbitals computed by a quantum chemistry program. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy Software > Simulation Methods

Gozem, Samer↗
DOE OSTI2021

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗