Search NASA⌕ Search

SEARCH · Search NASA

Results for “orbital dynamics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Mott localization in the van der Waals crystal CrI 3 : A GGA+DMFT study

Using the generalized gradient approximation plus dynamical mean-field theory (GGA+DMFT) we confirm the importance of multi-orbital dynamical correlations in determining the paramagnetic insulating state of CrI 3 . Furthermore, while the ferromagnetic phase reveals weak electronic correlation effects due to strong spin-orbital polarization, the Mott insulating state of paramagnetic CrI 3 crystal is shown to be driven by the interplay between orbital-dependent one-electron lineshape and multi-orbital electronic interactions. To probe the paramagnetic Mott insulating state we performed x-ray absorption spectroscopy (XAS) measurements for the two structural phases of CrI 3 . Our study is relevant to understanding the orbital-selective electronic structure reconstruction of Mott insulators and should be applicable to other van der Waals bonded materials from bulk to the ultrathin limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

MEGASIM: Lifetimes and Resonances of Earth Trojan Asteroids—The Death of Primordial ETAs?

Abstract We present an analysis of lifetimes and resonances of Earth Trojan Asteroids (ETAs) in the MEGASIM data set. Trojan asteroids co-orbit the Sun with a planet, but remain bound to the Lagrange points, L4 (60° leading the planet) or L5 (60° trailing). In the circular three-body approximation, the stability of a Trojan asteroid depends on the ratio of the host planet mass and the central mass. For the inner planets, the range of stability becomes increasingly small, so perturbations from the planets have made primordial Trojans rare. To date, there have been just two ETAs (2010 TK 7 and 2020 XL 5 ), several Mars Trojans, and a Venus Trojan discovered. The estimated lifetimes of the known inner system Trojans are shorter than a million years, suggesting they are interlopers rather than members of a stable and long-lasting population. With the largest ETA n -body simulation to date, we are able to track their survival across a wide initialized parameter space. We find that the remaining fraction of ETAs over time is well fit with a stretched exponential function that, when extrapolated beyond our simulation run time, predicts zero ETAs by 2.33 Gyr. We also show correlations between ETA ejections and the periods of the Milankovitch cycles. Though Earth’s orbital dynamics dominate the instabilities of ETAs, we provide evidence that ETA ejections are linked to resonances found in the variation of the orbital elements of many if not all of the planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Capture and migration of Jupiter and Saturn in mean motion resonance in a gaseous protoplanetary disc

Here, we study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets’ orbits, Δa SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inwards or outwards). We also provide an assessment of the planet’s orbital dynamics at different epochs of Saturn’s growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Δa SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the time-scales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low-probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Planetary waves can activate resonant drag instabilities in 3D dusty gaseous discs

Resonant drag instabilities (RDIs) in protoplanetary discs are driven by the aerodynamic back-reaction of dust on gas and occur when the relative dust-gas motion resonate with a wave mode intrinsic to the gas fluid. Axisymmetric models indicate that the RDI generates filamentary perturbations, leading to grain clumping and planetesimal formation. Motivated by these findings, we investigate the dust–gas interaction in a non-axisymmetric inviscid protoplanetary disc with an embedded low-mass planet ($M_{\mathrm{p}}\in [0.3, 3] \, \rm M_{\oplus }$, here $\, \rm M_{\oplus }$ is the Earth mass). We conduct global 3D high-resolution two-fluid simulations, with the dust being parametrized by the Stokes number $\mathrm{St}\in [0.01,0.5]$. We find that planetary waves (PWs; also known as Rossby waves), which propagate along the downstream separatrices of the horseshoe region, resonate with the streaming motion and trigger the RDI. The consequent development of a global-scale filamentary dust distribution does not sensitively depend on the Stokes number, nor does it depend on the fast dust settling that takes place in an inviscid disc. The rapid onset of this instability, which is comparable to the dynamical orbital time-scale, suppresses the formation of asymmetric structures in the dust in the vicinity of the planet (such as dust voids and filaments). Additionally, we find that the dust feedback enables buoyancy resonances in an otherwise non-buoyant (globally isothermal) disc. Therefore, our results provide the first numerical evidence of RDIs generation driven by planetary waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Optical Properties of the Infinite-Layer La 1-x Sr x NiO 2 and Hidden Hund's Physics

We investigate the optical properties of the normal state of the infinite-layer La 1-x Sr x NiO 2 using density functional theory plus dynamical mean-field theory. We find a correlated metal which exhibits substantial transfer of spectral weight to high energies relative to the density functional theory. The correlations are not due to Mott physics, which would suppress the charge fluctuations and the integrated optical spectral weight as we approach a putative insulating state. Instead, we find the unusual situation, that the integrated optical spectral weight decreases with doping and increases with increasing temperature. We contrast this with the coherent component of the optical conductivity, which decreases with increasing temperature as a result of a coherence-incoherence crossover. Our studies reveal that the effective crystal field splitting is dynamical and increases strongly at low frequency. This leads to a picture of a Hund's metallic state, where dynamical orbital fluctuations are visible at intermediate energies, while at low energies a Fermi surface with primarily d x 2 -y 2 character emerges. The infinite-layer nickelates are thus in an intermediate position between the iron based high temperature superconductors where multiorbital Hund's physics dominates and a one-band system such as the cuprates. To capture this physics we propose a low-energy two-band model with atom centered e g states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

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

Precise Dynamical Masses and Orbital Fits for β Pic b and β Pic c

We present a comprehensive orbital analysis to the exoplanets β Pictoris b and c that resolves previously reported tensions between the dynamical and evolutionary mass constraints on β Pic b. We use the Markov Chain Monte Carlo orbit code orvara to fit 15 years of radial velocities and relative astrometry (including recent GRAVITY measurements), absolute astrometry from Hipparcos and Gaia, and a single relative radial velocity measurement between β Pic A and b. We measure model-independent masses of 9.3{sub −2.5}{sup +2.6} M {sub Jup} for β Pic b and 8.3 ± 1.0 M {sub Jup} for β Pic c. These masses are robust to modest changes to the input data selection. We find a well-constrained eccentricity of 0.119 ± 0.008 for β Pic b, and an eccentricity of 0.21{sub −0.09}{sup +0.16} for β Pic c, with the two orbital planes aligned to within ∼05. Both planets’ masses are within ∼1σ of the predictions of hot-start evolutionary models and exclude cold starts. We validate our approach on N-body synthetic data integrated using REBOUND. We show that orvara can account for three-body effects in the β Pic system down to a level ∼5 times smaller than the GRAVITY uncertainties. Systematics in the masses and orbital parameters from orvara’s approximate treatment of multiplanet orbits are a factor of ∼5 smaller than the uncertainties we derive here. Future GRAVITY observations will improve the constraints on β Pic c’s mass and (especially) eccentricity, but improved constraints on the mass of β Pic b will likely require years of additional radial velocity monitoring and improved precision from future Gaia data releases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Solvated Nuclear–Electronic Orbital Structure and Dynamics

Nonadiabatic dynamical processes such as proton-coupled electron transfer and excited state intramolecular proton transfer have been the subject of much research. One of the promising theoretical methods to describe these processes is the nuclear–electronic orbital (NEO) approach. This approach inherently accounts for nuclear quantum effects within quantum chemistry calculations, and it has recently been extended to directly simulate nonadiabatic processes with the development of real-time NEO methods. These processes can also be significantly dependent on the surrounding chemical environment, however, and capturing the effects of the environment is often necessary for analyzing experimentally relevant systems. This work couples the NEO density functional theory and real-time time-dependent density functional theory approaches with solvation through the polarizable continuum model. Here, the effects of this coupling are investigated for ground state properties, solvent-dependent vibrational frequencies, and direct excited state intramolecular proton transfer dynamics.

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

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

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

Subhalos in Galaxy Clusters: Coherent Accretion and Internal Orbits

Subhalo dynamics in galaxy cluster host halos govern the observed distribution and properties of cluster member galaxies. We use the IllustrisTNG simulation to investigate the accretion and orbits of subhalos found in cluster-size halos. We find that the median change in the major axis direction of cluster-size host halos is approximately 80° between a ∼ 0.1 and the present day. We identify coherent regions in the angular distribution of subhalo accretion, and ∼68% of accreted subhalos enter their host halo through ∼38% of the surface area at the virial radius. The majority of galaxy clusters in the sample have ∼2 such coherent regions. We further measure angular orbits of subhalos with respect to the host major axis and use a clustering algorithm to identify distinct orbit modes with varying oscillation timescales. The orbit modes correlate with subhalo accretion conditions. Subhalos in orbit modes with shorter oscillations tend to have lower peak masses and accretion directions somewhat more aligned with the major axis. One orbit mode, exhibiting the least oscillatory behavior, largely consists of subhalos that accrete near the plane perpendicular to the host halo major axis. Our findings are consistent with expectations from inflow from major filament structures and internal dynamical friction: most subhalos accrete through coherent regions, and more massive subhalos experience fewer orbits after accretion. Our work offers a unique quantification of subhalo dynamics that can be connected to how the intracluster medium strips and quenches cluster galaxies.

galaxy clusters↗
DOE OSTI2023

Novel Phases and Dynamics of Multi-orbital Mott Insulators

The goal of this project is to gain theoretical insight into magnetism in the presence of spin-orbit coupling (SOC) for multi-orbital Mott insulators. Inspired by the exact solution of the Kitaev model that harbors a quantum spin liquid with novel excitations and struck by its relevance for materials with anisotropic orbital interactions, the PI will explore various new classes of 4d and 5d transition metal oxides. Starting with all electron Hamiltonians, the PI will derive minimal magnetic models. The aim will be to understand the role played by orbital frustration, even in the absence of any geometric frustration, in creating orbitally ordered as well as spin-orbital liquid phases. Orbital frustration arises primarily from the directional or anisotropic nature of d-orbitals in contrast to the isotropic nature of the spin degree of freedom. These models will be investigated by a variety of theoretical and numerical methods, including exact diagonalization, mean field theories, density matrix renormalization group and quantum Monte Carlo methods. Testable predictions for three experiments: nuclear magnetic resonance (NMR) and resonant x-ray scattering (RXS) to probe orbital ordering and pump-probe experiments to probe quantum dynamics will test the validity of these models for materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Spin-orbit coupling induced demagnetization in Ni: Ab initio nonadiabatic molecular dynamics perspective

Spin-orbit coupling (SOC), which can induce spin flip during the relaxation of photoexcited charge carrier, plays a crucial role in spin dynamics. Here in this work, we have used time-domain ab initio nonadiabatic molecular dynamics (NAMD) method to study the SOC induced ultrafast demagnetization in Ni at 300K. The spin-diabatic representation using spin-polarized Kohn-Sham (KS) basis sets and spin-adiabatic representation using spinor basis sets have been applied, and both of them achieve demagnetization in Ni with a timescale around 100fs. The spin-diabatic representation suggests a picture that the electron-phonon coupling (EPC) provides direct energy relaxation channel among the same-spin states, while the SOC can induce spin flip. After photoexcitation, it is found the spin-minority electrons relax to the same-spin states rather than the opposite-spin states, since EPC is larger than SOC by one order of magnitude. By contrast, for the spin-majority electrons, spin flip occurs since there are no empty same-spin states as electron acceptor above the Fermi level. The different relaxation pathways for spin-majority and spin-minority electrons induce the demagnetization. The spin-adiabatic representation provides an Elliott-Yafet spin-phonon scattering picture. The SOC induced reduction of magnetic moment in Ni may induce magnon to drive further demagnetization. The ab initio NAMD simulation provides a critical angle to understand how the SOC and EPC affect demagnetization process in Ni.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum Sensing and Imaging of Spin–Orbit-Torque-Driven Spin Dynamics in the Non-Collinear Antiferromagnet Mn 3 Sn

Novel non-collinear antiferromagnets with spontaneous time-reversal symmetry breaking, non-trivial band topology, and unconventional transport properties have received immense research interest over the past decade due to their rich physics and enormous promise in technological applications. One of the central focuses in this emerging field is exploring the relationship between the microscopic magnetic structure and exotic material properties. Nanoscale imaging of both spin–orbit-torque-induced deterministic magnetic switching and chiral spin rotation in non-collinear antiferromagnet Mn 3 Sn films using nitrogen-vacancy (NV) centers are reported. Direct evidence of the off-resonance dipole–dipole coupling between the spin dynamics in Mn 3 Sn and proximate NV centers is also demonstrated by NV relaxometry measurements. These results demonstrate the unique capabilities of NV centers in accessing the local information of the magnetic order and dynamics in these emergent quantum materials and suggest new opportunities for investigating the interplay between topology and magnetism in a broad range of topological magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗