Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orbital motion”

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 109 records · Page 6

Models of bars − II. Exponential profiles

ABSTRACT We present a new model for galactic bars with exponentially falling major axis luminosity profiles and Gaussian cross-sections. This is based on the linear superposition of Gaussian potential–density pairs with an exponential weight function, using an extension of the method originally introduced by Long & Murali. We compute the density, potential, and forces, using Gaussian quadrature. These quantities are given as explicit functions of position. There are three independent scaled bar parameters that can be varied continuously to produce bespoke bars of a given mass and shape. We categorize the effective potential by splitting a reduced parameter space into six regions. Unusually, we find bars with three stable Lagrange points on the major axis are possible. Our model reveals a variety of unexpected orbital structure, including a bifurcating x1 orbit coexisting with a stable x4 orbit. Propeller orbits are found to play a dominant role in the orbital structure, and we find striking similarities between our bar configuration and the model of Kaufmann & Contopoulos. We find a candidate orbital family, sired from the propeller orbits, that may be responsible for the observed high-velocity peaks in the Milky Way’s bar. As a cross-check, we inspect, for the first time, the proper motions of stars in the high-velocity peaks, which also match our suggested orbital family well. This work adds to the increasing body of evidence that real galactic bars may be supported at least partly by propeller orbits rather than solely by elliptical-like orbits of the x1 family.

McGough, D. P.↗

Real-time time-dependent density functional theory in Quantum Espresso and related codes

A set of subroutines that interact with an existing computer program, Quantum Espresso (QE), that takes a set of electronic orbitals from QE and propagates them forward in time in response to a number of perturbations. These perturbations include the motion of ions and a wide spectrum of applied electromagnetic fields (e.g., static magnetic fields, slowly-varying electric fields, x-rays, etc.). From the time evolution of the electronic orbitals comes time-varying values of physical observables like the electronic charge and current densities and atomic forces. These are further post-processed to yield properties of interest like stopping powers, conductivities, and the dynamic structure factor. The electronic orbitals are described in a plane wave basis and the equations of motion are carried out using a unitary or approximately unitary time propagation algorithm. Related work will include subroutines for interfacing this capability, as well as others that already exist within QE, with a fork of QE that is more extensible and modular.SAND2019-12611 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Baczewski, Andrew↗

Effects of turbulent diffusion and back-reaction on the dust distribution around two resonant planets

In evolved and dusty circumstellar discs, two planets with masses comparable to Jupiter and Saturn that migrate outwards while maintaining an orbital resonance can produce distinctive features in the dust distribution. Dust accumulates at the outer edge of the common gas gap, which behaves as a dust trap, where the local dust concentration is significantly enhanced by the planets’ outward motion. Concurrently, an expanding cavity forms in the dust distribution inside the planets’ orbits, because dust does not filter through the common gaseous gap and grain depletion in the region continues via inward drifting. There is no cavity in the gas distribution because gas can filter through the gap, although ongoing gas accretion on the planets can reduce the gas density in the inner disc. Such behaviour was demonstrated by means of simulations neglecting the effects of dust diffusion due to turbulence and of dust backreaction on the gas. Both effects may alter the formation of the dust peak at the gap outer edge and of the inner dust cavity, by letting grains filter through the dust trap. We performed high-resolution hydrodynamical simulations of the coupled evolution of gas and dust species, the latter treated as pressureless fluids, in the presence of two giant planets. We show that diffusion and backreaction can change some morphological aspects of the dust distribution but do not alter some main features, such as the outer peak and the expanding inner cavity. Furthermore, these findings are confirmed for different parametrizations of gas viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗

The Milky Way’s stellar streams and globular clusters do not align in a Vast Polar Structure

ABSTRACT There is increasing evidence that a substantial fraction of Milky Way satellite galaxies align in a rotationally supported plane of satellites, a rare configuration in cosmological simulations of galaxy formation. It has been suggested that other Milky Way substructures (namely young halo globular clusters and stellar/gaseous streams) similarly tend to align with this plane, accordingly dubbed the Vast Polar Structure (VPOS). Using systemic proper motions inferred from Gaia data, we find that globular cluster orbital poles are not clustered in the VPOS direction, though the population with the highest VPOS membership fraction is the young halo clusters (∼30 per cent). We additionally provide a current census of stellar streams, including new streams discovered using the Dark Energy Survey and Gaia data sets, and find that stellar stream normals are also not clustered in the direction of the VPOS normal. We also find that, based on orbit modelling, there is a likely association between NGC 3201 and the Gjöll stellar stream and that, based on its orbital pole, NGC 4147 is likely not a Sagittarius globular cluster. That the Milky Way’s accreted globular clusters and streams do not align in the same planar configuration as its satellites suggests that the plane of satellites is either a particularly stable orbital configuration or a population of recently accreted satellites. Neither of these explanations is particularly likely in light of other recent studies, leaving the plane of satellites problem as one of the more consequential open problems in galaxy formation and cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Meteoroid Stream of Comet C/1961 T1 (Seki) and Its Relation to the December ρ-Virginids and γ-Sagittariids

We answer the question of whether long-period comet C/1961 T1 (Seki) is the parent body of a meteor shower observed in Earth’s atmosphere. We model discrete parts of a meteoroid stream associated with this comet and follow the orbital evolution of its particles via a numerical integration of their orbits. The gravitational perturbations of major planets and the Poynting–Robertson effect, which influences the motion of the particles, are taken into account. We predict two showers which are related to C/1961 T1. One of these showers was identified with the December ρ-Virginids, #502. This identification confirms an earlier suggested relationship between the shower and the comet. The second predicted shower was identified with the γ-Sagittariids, #657. The radiant areas of both showers are symmetrical in respect to the apex of Earth’s heliocentric motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Pegasus IV: Discovery and Spectroscopic Confirmation of an Ultra-faint Dwarf Galaxy in the Constellation Pegasus

We report the discovery of Pegasus IV, an ultra-faint dwarf galaxy found in archival data from the Dark Energy Camera processed by the DECam Local Volume Exploration Survey. Pegasus IV is a compact, ultra-faint stellar system ( pc; M $_{V}$ = −4.25 ± 0.2 mag) located at a heliocentric distance of . Based on spectra of seven nonvariable member stars observed with Magellan/IMACS, we confidently resolve Pegasus IV’s velocity dispersion, measuring km s$^{−1}$ (after excluding three velocity outliers); this implies a mass-to-light ratio of for the system. From the five stars with the highest signal-to-noise spectra, we also measure a systemic metallicity of [Fe/H] = dex, making Pegasus IV one of the most metal-poor ultra-faint dwarfs. We tentatively resolve a nonzero metallicity dispersion for the system. These measurements provide strong evidence that Pegasus IV is a dark-matter-dominated dwarf galaxy, rather than a star cluster. We measure Pegasus IV’s proper motion using data from Gaia Early Data Release 3, finding (μ $_{α*}$, μ $_{δ}$) = (0.33 ± 0.07, −0.21 ± 0.08) mas yr$^{−1}$. When combined with our measured systemic velocity, this proper motion suggests that Pegasus IV is on an elliptical, retrograde orbit, and is currently near its orbital apocenter. Lastly, we identify three potential RR Lyrae variable stars within Pegasus IV, including one candidate member located more than 10 half-light radii away from the system’s centroid. The discovery of yet another ultra-faint dwarf galaxy strongly suggests that the census of Milky Way satellites is still incomplete, even within 100 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗

A Comprehensive Analysis of NGC 2158 in the Gaia Era: Photometric Parameters, Apex, and Orbit

We present an investigation of NGC 2158 using Gaia DR2 data. We identified 800 most likely cluster members with membership probability higher than 90%. The mean proper motions of this object are determined as (μ {sub x} = −0.203 ± 0.003, μ {sub y} = −1.99 ± 0.004) mas yr{sup −1}. The limiting radius, log(age), and distance of the cluster are obtained as 23.5 arcmin, 9.38 ± 0.04 Gyr, and 4.69 ± 0.22 kpc, respectively. The overall mass function slope (0.93 ± 0.14) is flatter than the Salpeter value (1.35) within the mass range of 1.17–1.44 M {sub ⊙}. This cluster also shows the mass-segregation effect and our study demonstrates that NGC 2158 is a dynamically relaxed open cluster. Using the apex diagram (AD), the apex coordinates of the cluster are obtained in different ways and are examined using the (μ {sub U},μ {sub T}) diagram. The best value of the apex coordinates is determined as A = 8724 ± 160, D = −3661 ± 530. We also determined the orbit of the cluster and found that NGC 2158 moves almost in the solar antapex direction. The resulting spatial velocity of NGC 2158 is 51 km s{sup −1}. A significant oscillation along the Z-coordinate up to 529 pc is detected. Various scenarios regarding the origin of this cluster are also discussed.

42 ENGINEERING↗

A KELT–TESS Eclipsing Binary in a Young Triple System Associated with the Local “Stellar String” Theia 301

HD 54236 is a nearby, wide common-proper-motion visual pair that has been previously identified as likely being very young by virtue of strong X-ray emission and lithium absorption. Here, we report the discovery that the brighter member of the wide pair, HD 54236A, is itself an eclipsing binary (EB), comprising two near-equal solar-mass stars on a 2.4 days orbit. It represents a potentially valuable opportunity to expand the number of benchmark-grade EBs at young stellar ages. Using new observations of Ca ii H and K emission and lithium absorption in the wide K-dwarf companion, HD 54236B, we obtain a robust age estimate of 225 ± 50 Myr for the system. This age estimate and Gaia proper motions show HD 54236 is associated with Theia 301, a newly discovered local “stellar string,” which itself may be related to the AB Dor moving group through shared stellar members. Applying this age estimate to AB Dor itself alleviates reported tension between observation and theory that arises for the luminosity of AB Dor C when younger age estimates are used.

79 ASTRONOMY AND ASTROPHYSICS↗

Skyrmion soliton motion on periodic substrates by atomistic and particle-based simulations

Abstract We compare the dynamical behavior of magnetic skyrmions interacting with square and triangular defect arrays just above commensuration using both atomistic and particle-based models. Under applied drives, the initial motion is a kink traveling through the pinned skyrmion lattice. For the square defect array, both models agree well and show a regime in which the soliton motion is locked along 45°. The atomistic model also produces locking of a soliton along 30°, which is absent in the particle-based model. For the triangular defect array, the atomistic model exhibits 30° soliton motion over a wide region of external current values. In contrast, the particle-based model gives 45° soliton motion over a small range of external driving force values. The difference arises because the nondeforming particle model facilitates meandering skyrmion orbits while the deformable atomistic model enables stronger skyrmion-skyrmion interactions that reduce the meandering. Our results indicate that soliton motion through pinned skyrmion lattices on a periodic substrate is a robust effect.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spectroscopic Analysis of the Bulge Globular Cluster ESO 456-SC38

It has been suggested that the oldest stellar populations in the Milky Way are tightly bound and confined to the central regions of the Galaxy. This is one of the reasons why a handful of globular clusters located in the bulge region are thought to be remnants of the primeval formation stages of the Milky Way. The globular cluster, ESO 456-SC38 (Djorgovski 2), is one such cluster; it has a blue horizontal branch, is projected very close to the center of the Galaxy, and has an orbit confining it to the bulge/bar region. The first α abundances of seven stars in this heavily reddened cluster are presented using Apache Point Observatory Galaxy Evolution Experiment Data Release 16. A significant spread in the abundances of N, C, Na, and Al indicates the presence of multiple stellar populations in this cluster. Using Gaia DR2 proper motions and radial velocities from the Bulge Radial Velocity Assay for RR Lyrae stars, we confirm that RR Lyrae stars belong to this globular cluster.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of state-interaction and spinor-representation calculations of spin-orbit coupling within exact two-component coupled-cluster theories

Here a benchmark study of state-interaction and spinor-representation calculations of spin-orbit coupling using the exact two-component Hamiltonians with atomic mean-field integrals (the X2CAMF schemes) at the equation-of-motion coupled-cluster singles and doubles level is reported. We adopt a version of the X2CAMF scheme with spin-orbit integrals correct to first order for the state-interaction calculations and a version correct to infinite order for the spinor-representation calculations. The differences between the state-interaction calculations with minimum active spaces to account for (quasi-)degeneracy and the spinor-representation calculations thus correspond to higher-order spin-orbit contributions. This state-interaction approach with scalar-relativistic effects accurately included in the reference functions and in the spin-orbit integrals is shown to exhibit robust performance for elements across the periodic table. On the other hand, the more rigorous spinor representation shows more rapid convergence with respect to the number of correlated electrons and is the preferred choice for accurate calculations for heavy elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From synchrotrons for XFELs: the soft x-ray near-edge spectrum of the ESCA molecule

A predictive understanding of soft x-ray near-edge absorption spectra of small molecules is an enduring theoretical challenge and of current interest for x-ray probes of molecular dynamics. We report the experimental absorption spectrum for the electron spectroscopy for chemical analysis (ESCA) molecule (ethyl trifluoroacetate) near the carbon 1s absorption edge between 285-300 eV. Furthermore, the ESCA molecule with four chemically distinct carbon sites has previously served as a theoretical benchmark for photoelectron spectra and now for photoabsorption spectra. We report a simple edge-specific approach for systematically expanding standard basis sets to properly describe diffuse Rydberg orbitals and the importance of triple excitations in equation-of-motion coupled-cluster calculations of the energy interval between valence and Rydberg excitations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Discovery and Spectroscopic Confirmation of Aquarius III: A Low-mass Milky Way Satellite Galaxy

Abstract We present the discovery of Aquarius III, an ultra-faint Milky Way satellite galaxy identified in the second data release of the DECam Local Volume Exploration survey. Based on deeper follow-up imaging with DECam, we find that Aquarius III is a low-luminosity ( M V = − 2.5 − 0.5 + 0.3 ; L V = 850 − 260 + 380 L ⊙ ), extended ( r 1 / 2 = 41 − 8 + 9 pc) stellar system located in the outer halo ( D ⊙ = 85 ± 4 kpc). From medium-resolution Keck/DEIMOS spectroscopy, we identify 11 member stars and measure a mean heliocentric radial velocity of v sys = − 13.1 − 0.9 + 1.0 km s − 1 for the system and place an upper limit of σ v < 3.5 km s −1 ( σ v < 1.6 km s −1 ) on its velocity dispersion at the 95% (68%) credible level. Based on calcium-triplet metallicities of the six brightest red giant members, we find that Aquarius III is very metal-poor ([Fe/H]= − 2.61 ± 0.21) with a statistically significant metallicity spread ( σ [ Fe / H ] = 0.46 − 0.14 + 0.26 dex). We interpret this metallicity spread as strong evidence that the system is a dwarf galaxy as opposed to a star cluster. Combining our velocity measurement with Gaia proper motions, we find that Aquarius III is currently situated near its orbital pericenter in the outer halo ( r peri = 78 ± 7 kpc) and that it is plausibly on first infall onto the Milky Way. This orbital history likely precludes significant tidal disruption from the Galactic disk, notably unlike other satellites with comparably low velocity dispersion limits in the literature. Thus, if further velocity measurements confirm that its velocity dispersion is truly below σ v ≲ 2 km s −1 , Aquarius III may serve as a useful laboratory for probing galaxy formation physics in low-mass halos.

79 ASTRONOMY AND ASTROPHYSICS↗

HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b: Four Warm Giant Planets Transiting K Dwarfs

We report the discovery of four transiting giant planets around K dwarfs. The planets HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b have masses of 0.369{sub −0.021}{sup +0.031} M{sub J}, 0.243{sub −0.030}{sup +0.022} M{sub J}, 0.353{sub −0.027}{sup +0.038} M{sub J}, and 0.1254±0.0039 M{sub J}, respectively, and radii of 1.117±0.014 R{sub J}, 0.800±0.015 R{sub J}, 0.765±0.013 R{sub J}, and 0.7224±0.0032 R{sub J}, respectively. The planets orbit close to their host stars with orbital periods of 3.9228 days, 3.1317 days, 4.1480 days, and 7.3279 days, respectively. The hosts are main-sequence K dwarfs with masses of 0.674{sub −0.012}{sup +0.016} M{sub ⊙}, 0.7279 ±0.0066 M{sub ⊙}, 0.7133±0.0075 M{sub ⊙}, and 0.7311±0.0028, and with V-band magnitudes of V=14.829±0.010, 14.35±0.11, 14.998±0.040 and 12.469±0.010. The super-Neptune HATS-72b (a.k.a. WASP-191b and TOI 294.01) was independently identified as a transiting planet candidate by the HATSouth, WASP, and TESS surveys, and we present a combined analysis of all of the data gathered by each of these projects (and their follow-up programs). An exceptionally precise mass is measured for HATS-72b thanks to high-precision radial velocity (RV) measurements obtained with VLT/ESPRESSO, FEROS, HARPS, and Magellan/PFS. We also incorporate TESS observations of the warm Saturn–hosting systems HATS-47 (a.k.a. TOI 1073.01), HATS-48A, and HATS-49. HATS-47 was independently identified as a candidate by the TESS team, while the other two systems were not previously identified from the TESS data. The RV orbital variations are measured for these systems using Magellan/PFS. HATS-48A has a resolved [[Formula]] neighbor in Gaia DR2, which is a common-proper-motion binary star companion to HATS-48A with a mass of 0.22 M{sub ⊙} and a current projected physical separation of ∼1400 au.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

HD 143811 AB b: A Directly Imaged Planet Orbiting a Spectroscopic Binary in Sco-Cen

We present confirmation of HD 143811 AB b, a substellar companion to spectroscopic binary HD 143811 AB through direct imaging with the Gemini Planet Imager (GPI) and Keck NIRC2. HD 143811 AB was observed as a part of the GPI Exoplanet Survey in 2016 and 2019 and is a member of the Sco-Cen star formation region. The exoplanet is detected ∼430 mas from the host star by GPI. With two GPI epochs and one from Keck/NIRC2 in 2022, we confirm through common proper motion analysis that the object is bound to its host star. We derive an orbit with a semimajor axis of $64^{+32}_{-14}$ au and eccentricity $0.23^{+0.24}_{-0.16}$. Spectral analysis of the GPI H-band spectrum and NIRC2 L′ photometry provides additional proof that this object is a substellar companion. We compare the spectrum of HD 143811 AB b to PHOENIX stellar models and Exo-Radioactive-Convective Equilibrium Model (REM) exoplanet atmosphere models and find that Exo-REM models provide the best fits to the data. From the Exo-REM models, we derive an effective temperature of $1042^{+178}_{-132}$ K for the planet and translate the derived luminosity of the planet to a mass of 5.6 ± 1.1 M Jup assuming hot-start evolutionary models. HD 143811 AB b is the first directly imaged planet around a binary that is not on an ultrawide orbit. Future characterization of this object will shed light on the formation of planets around binary star systems.

Astronomy and AstroPhysics↗

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of the Impact of Flow on MHD Perturbations with the NIMROD Code (Final Report)

This work will enable improved understanding of tokamak edge flows and impurity-species dynamics and their impact on MHD stability as it relates to ELMs, RMPs and QH-mode by leveraging computational developments of the NIMROD code. Experimentally, it is established that the tokamak flows and impurities species can have a substantial effect on MHD stability. This is particularly true for the edge plasma where, for example, large flow shear is correlated with the occurrence of the Quiescent H-mode (QH- mode) state as opposed to operation with edge localized modes (ELMs). Understanding plasma flows is challenging because of the many physical effects that come into play: particle orbits loss, neutrals dynamics, and interaction between multiple ion species with decoupled motion. We will develop a model that incorporates multiple ion species into the MHD framework that includes self-consistent magnetic-field evolution. With this model we will understand how the presence of impurities impacts the tokamak edge-pedestal flows. Finally, we will quantify the impact of the new multiple species modeling on 3D QH-mode simulations. This work will have broader impacts on all MHD studies when multiple species are present (as is the case in modern tokamaks), as well as contributing to the edge and transport communities who also are interested in the studies of multiple-ion-species flows in the edge pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗