Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinematics and 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measuring the streaming motion in the Milky Way disc with Gaia EDR3+

We map the 3D kinematics of the Galactic disc out to 3.5 kpc from the Sun, and within 0.75 kpc from the mid-plane of the Galaxy. To this end, we combine high-quality astrometry from Gaia EDR3, with heliocentric line-of-sight velocities from Gaia DR2, and spectroscopic surveys including APOGEE, GALAH, and LAMOST. We construct an axisymmetric model for the mean velocity field, and subtract this on a star-by-star basis to obtain the residual velocity field in the Galactocentric components (V Φ , V R , V z ), and V los . The velocity residuals are quantified using the power spectrum, and we find that the peak power (A/[km s –1 ]) in the mid-plane (|z| < 0.25 kpc) is (A Φ , A R , A Z , A los ) = (4.2,8.5,2.6,4.6), at 0.25 < |z|/[kpc] < 0.5, is (A Φ , A R , A Z , A los ) = (4.0,7.9,3.6,5.3), and at 0.5 < |z|/[kpc] < 0.75, is (A Φ , A R , A Z , A los ) = (1.9,6.9,5.2,6.4). Our results provide a sophisticated measurement of the streaming motion in the disc and in the individual components. We find that streaming is most significant in VR, and at all heights (|Z|) probed, but is also non-negligible in other components. Additionally, we find that patterns in velocity field overlap spatially with models for spiral arms in the Galaxy. Our simulations show that phase-mixing of disrupting spiral arms can generate such residuals in the velocity field, where the radial component is dominant, just as in real data. Furthermore, we also find that with time evolution, both the amplitude and physical scale of the residual motion decrease.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The stellar mass of the Gaia -Sausage/Enceladus accretion remnant

The Gaia-Sausage/Enceladus (GS/E) structure is an accretion remnant that comprises a large fraction of the Milky Way’s stellar halo. We study GS/E using high-purity samples of kinematically selected stars from APOGEE DR16 and Gaia. Employing a novel framework to account for kinematic selection biases using distribution functions, we fit density profiles to these GS/E samples and measure their masses. We find that GS/E has a shallow density profile in the inner Galaxy, with a break between 15 and 25 kpc beyond which the profile steepens. We also find that GS/E is triaxial, with axis ratios 1:0.55:0.45 (nearly prolate), and the major axis is oriented about 80° from the Sun–Galactic centre line and 16° above the plane. We measure a stellar mass for GS/E of $1.45\, ^{+0.92}_{-0.51}\, \mathrm{(stat.)}\, ^{+0.13}_{-0.37} \mathrm{(sys.)}\ \times 10^{8}$ M ⊙ . Our mass estimate is lower than others in the literature, a finding we attribute to the excellent purity of the samples we work with. We also fit a density profile to the entire Milky Way stellar halo, finding a mass in the range of 6.7–8.4 × 10 8 M ⊙ , and implying that GS/E could make up as little as 15–25 per cent of the mass of the Milky Way stellar halo. Our lower stellar mass combined with standard stellar mass-to-halo mass relations implies that GS/E constituted a minor 1:8 mass-ratio merger at the time of its accretion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A comparison of the baryonic Tully–Fisher relation in MaNGA and IllustrisTNG

We compare an observed baryonic Tully–Fisher relation (BTFR) from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) and HI-MaNGA surveys to a simulated BTFR from the cosmological magnetohydrodynamical simulation, IllustrisTNG. To do so, we calibrate the BTFR of the local universe using 377 galaxies from the MaNGA and HI-MaNGA surveys, and perform mock 21 cm observations of matching galaxies from IllustrisTNG. The mock observations are used to ensure that the comparison with the observed galaxies is fair since it has identical measurement algorithms, observational limitations, biases, and uncertainties. For comparison, we also calculate the BTFR for the simulation without mock observations and demonstrate how mock observations are necessary to fairly and consistently compare between observational and theoretical data. We report a MaNGA BTFR of log$_{10} (M_{ \rm Bary}/M_\odot)= (2.97 \pm 0.18)$ log$_{10} V_{ \rm Rot} + (4.04 \pm 0.41)\, \log _{10}{M_{\odot }}$ and an IllustrisTNG BTFR of log$_{10} (M_{ \rm Bary}/M_\odot) = (2.94 \pm 0.23$) log$_{10} V_{ \rm Rot} + (4.15 \pm 0.44)\, \log _{10}{M_{\odot }}$. Thus, MaNGA and IllustrisTNG produce BTFRs that agree within uncertainties, demonstrating that IllustrisTNG has created a galaxy population that obeys the observed relationship between mass and rotation velocity in the observed universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Galaxy Zoo: kinematics of strongly and weakly barred galaxies

ABSTRACT We study the bar pattern speeds and corotation radii of 225 barred galaxies, using integral field unit data from MaNGA and the Tremaine–Weinberg method. Our sample, which is divided between strongly and weakly barred galaxies identified via Galaxy Zoo, is the largest that this method has been applied to. We find lower pattern speeds for strongly barred galaxies than for weakly barred galaxies. As simulations show that the pattern speed decreases as the bar exchanges angular momentum with its host, these results suggest that strong bars are more evolved than weak bars. Interestingly, the corotation radius is not different between weakly and strongly barred galaxies, despite being proportional to bar length. We also find that the corotation radius is significantly different between quenching and star-forming galaxies. Additionally, we find that strongly barred galaxies have significantly lower values for $\mathcal {R}$, the ratio between the corotation radius and the bar radius, than weakly barred galaxies, despite a big overlap in both distributions. This ratio classifies bars into ultrafast bars ($\mathcal {R} \lt $ 1.0; 11 per cent of our sample), fast bars (1.0 $\lt \mathcal {R} \lt $ 1.4; 27 per cent), and slow bars ($\mathcal {R} \gt $ 1.4; 62 per cent). Simulations show that $\mathcal {R}$ is correlated with the bar formation mechanism, so our results suggest that strong bars are more likely to be formed by different mechanisms than weak bars. Finally, we find a lower fraction of ultrafast bars than most other studies, which decreases the recently claimed tension with Lambda cold dark matter. However, the median value of $\mathcal {R}$ is still lower than what is predicted by simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The non-explosive stellar merging origin of the ultra-massive carbon-rich white dwarfs

We have investigated the origin of a sub-class of carbon-polluted white dwarfs (DQ) originally identified as the “hot DQ” white dwarfs. These objects are relatively hot ($10\, 000\lesssim T_{\rm eff}\lesssim 25\, 000$ K), have markedly higher carbon abundance (C-enriched), are more massive (M ≳ 0.8 M⊙) than ordinary DQs (M ∼ 0.6 M⊙), and display high space velocities. Hence, despite their young appearance their kinematic properties are those of an old white dwarf population. The way out of this dilemma is to assume that they formed via the merging of two white dwarfs. In this paper, we examine the observed characteristics of this population of “C-enriched” DQ white dwarfs and confirm that nearly half of the 63 known objects have kinematic properties consistent with those of the Galactic thick disc or halo. We have also conducted population synthesis studies and found that the merging hypothesis is indeed compatible with observations. Studies of this sub-class of white dwarfs have important implications for our understanding of Type Ia Supernovae (SNeIa), commonly used to determine the expansion history of the Universe, since the same formation channel applies to both kinds of objects. Hence, probing the properties of these white dwarfs that failed to explode may yield important constraints to the modelling of the mechanisms leading to a thermonuclear runaway.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Robust clustering of the local Milky Way stellar kinematic substructures with Gaia eDR3

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm HDBSCAN on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the HDBSCAN member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Probabilistic inference of the structure and orbit of Milky Way satellites with semi-analytic modelling

Semi-analytic modelling furnishes an efficient avenue for characterizing dark matter haloes associated with satellites of Milky Way-like systems, as it easily accounts for uncertainties arising from halo-to-halo variance, the orbital disruption of satellites, baryonic feedback, and the stellar-to-halo mass (SMHM) relation. We use the SatGen semi-analytic satellite generator, which incorporates both empirical models of the galaxy–halo connection as well as analytic prescriptions for the orbital evolution of these satellites after accretion onto a host to create large samples of Milky Way-like systems and their satellites. By selecting satellites in the sample that match observed properties of a particular dwarf galaxy, we can infer arbitrary properties of the satellite galaxy within the cold dark matter paradigm. For the Milky Way’s classical dwarfs, we provide inferred values (with associated uncertainties) for the maximum circular velocity v max and the radius r max at which it occurs, varying over two choices of baryonic feedback model and two prescriptions for the SMHM relation. While simple empirical scaling relations can recover the median inferred value for v max and r max , this approach provides realistic correlated uncertainties and aids interpretability. We also demonstrate how the internal properties of a satellite’s dark matter profile correlate with its orbit, and we show that it is difficult to reproduce observations of the Fornax dwarf without strong baryonic feedback. Furthermore, the technique developed in this work is flexible in its application of observational data and can leverage arbitrary information about the satellite galaxies to make inferences about their dark matter haloes and population statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Exploring the interaction between the MW and LMC with a large sample of blue horizontal branch stars from the DESI survey

The Large Magellanic Cloud (LMC) is a Milky Way (MW) satellite that is massive enough to gravitationally attract the MW disc and inner halo, causing significant motion of the inner MW with respect to the outer halo. In this work, we probe this interaction by constructing a sample of 9866 blue horizontal branch (BHB) stars with radial velocities from the DESI spectroscopic survey out to 120 kpc from the Galactic centre. This is the largest spectroscopic set of BHB stars in the literature to date, and it contains four times more stars with Galactocentric distances beyond 50 kpc than previous BHB catalogues. Using the DESI BHB sample combined with SDSS BHBs, we measure the bulk radial velocity of stars in the outer halo and observe that the velocity in the Southern Galactic hemisphere is different by 3.7σ from the North. Modelling the projected velocity field shows that its dipole component is directed at a point 22 deg away from the LMC along its orbit, which we interpret as the travel direction of the inner MW. The velocity field includes a monopole term that is –24 km s –1 ⁠, which we refer to as compression velocity. This velocity is significantly larger than predicted by the current models of the MW and LMC interaction. This work uses DESI data from its first 2 yr of observations, but we expect that with upcoming DESI data releases, the sample of BHB stars will increase and our ability to measure the MW–LMC interaction will improve significantly.

Galaxy: evolution↗
DOE OSTI2025

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗
DOE OSTI2025

Scaling relations in the phase-space structure of dark matter haloes

We present new scaling relations for the isotropic phase-space distribution functions (DFs) and energy distributions of simulated dark matter haloes. These relations are inspired by those for the singular isothermal sphere with density profile $\rho (r)\propto r^{-2}$, for which the DF satisfies $f(E) \propto r_{\mathrm{max}}^{-2}(E)$ and the energy distribution satisfies $\mathrm{ d}M/\mathrm{ d}E \propto r_{\mathrm{max}}(E)$, with $r_{\mathrm{max}}(E)$ being the radius where the gravitational potential equals energy $E$. For the simulated haloes, we find $f(E)\propto r_{\mathrm{max}}^{-2.08}(E)$ and $\mathrm{ d}M/\mathrm{ d}E \propto r_{\mathrm{max}}(E)$ across broad energy ranges. In addition, the proportionality coefficients depend on the gravitational constant and the parameters of the best-fitting Navarro–Frenk–White density profile. These scaling relations are satisfied by haloes over a wide mass range and provide an efficient method to approximate their DFs and energy distributions. Understanding the origin of these relations may shed more light on halo formation.

galaxies: haloes↗
DOE OSTI2024

Bayesian inference of structured latent spaces from neural population activity with the orthogonal stochastic linear mixing model

The brain produces diverse functions, from perceiving sounds to producing arm reaches, through the collective activity of populations of many neurons. Determining if and how the features of these exogenous variables (e.g., sound frequency, reach angle) are reflected in population neural activity is important for understanding how the brain operates. Often, high-dimensional neural population activity is confined to low-dimensional latent spaces. However, many current methods fail to extract latent spaces that are clearly structured by exogenous variables. This has contributed to a debate about whether or not brains should be thought of as dynamical systems or representational systems. Here, we developed a new latent process Bayesian regression framework, the orthogonal stochastic linear mixing model (OSLMM) which introduces an orthogonality constraint amongst time-varying mixture coefficients, and provide Markov chain Monte Carlo inference procedures. We demonstrate superior performance of OSLMM on latent trajectory recovery in synthetic experiments and show superior computational efficiency and prediction performance on several real-world benchmark data sets. We primarily focus on demonstrating the utility of OSLMM in two neural data sets: μ ECoG recordings from rat auditory cortex during presentation of pure tones and multi-single unit recordings form monkey motor cortex during complex arm reaching. We show that OSLMM achieves superior or comparable predictive accuracy of neural data and decoding of external variables (e.g., reach velocity). Most importantly, in both experimental contexts, we demonstrate that OSLMM latent trajectories directly reflect features of the sounds and reaches, demonstrating that neural dynamics are structured by neural representations. Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale biological time-series datasets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Spatial Decorrelation of Young Stars and Dense Gas as a Probe of the Star Formation–Feedback Cycle in Galaxies

The spatial decorrelation of dense molecular gas and young stars observed on ≲1 kpc scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. In this work, we explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC 300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC 300, and we show that this agreement is due to different aspects of feedback, including H 2 dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affects the correlation on ≲100 pc scales, while supernovae play a significant role on ≳ 100 pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per free fall time, ϵ ff , which provides a strong observational constraint on ϵ ff when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

Abstract We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard–Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in a collisionless system, also known as the Case–van Kampen modes, is eliminated. Furthermore, because the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest for N -body simulations as well as observations of gravitational systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The AGORA High-resolution Galaxy Simulations Comparison Project. V. Satellite Galaxy Populations in a Cosmological Zoom-in Simulation of a Milky Way–Mass Halo

We analyze and compare the satellite halo populations at z ~ 2 in the high-resolution cosmological zoom-in simulations of a 10 12 M ⊙ target halo (z = 0 mass) carried out on eight widely used astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Arepo-t, and Gizmo) for the AGORA High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near z = 2 for each code (hereafter "z ~ 2") at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the CosmoRun simulations and the DMO simulations, we discover that each CosmoRun halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at z ~ 2, the number of satellite galaxies is significantly fewer than that of dark matter halos in all participating AGORA simulations and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called "missing satellite problem" is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in AGORA and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution (≲100 proper pc at z = 2). We also compare other properties such as the stellar mass–halo mass relation and the mass–metallicity relation. Our work highlights the value of comparison studies such as AGORA, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

MaNGA DynPop – I. Quality-assessed stellar dynamical modelling from integral-field spectroscopy of 10K nearby galaxies: a catalogue of masses, mass-to-light ratios, density profiles, and dark matter

ABSTRACT This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here, we present a catalogue of dynamically determined quantities for over 10 000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically aligned JAMcyl or spherically aligned JAMsph. We also make four assumptions for the models’ dark versus luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results.

Astronomy & Astrophysics↗
DOE OSTI2022

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Theoretical Approach for the Fast Estimation of the Turbulent Kinematic Viscosity for Internal Flows

An analytical expression for turbulent kinematic viscosity (νt), based solely on the hydraulic Reynolds number (Re), was derived and evaluated. The analytical expression is valid for the fast estimation of νt for internal, isotropic, fully developed flows. The expression was compared with experimental and simulation data for air, water, and liquid sodium, and was shown to provide reasonable values for 2100 ≤ Re ≤ 3.6 × 106 and Prandtl number (Pr) range of 0.0107 ≤ Pr ≤ 5.65. In addition, new expressions suitable for the central portion of internal flows, away from the wall, were derived for the turbulent Reynolds number (ReT), showing its relationship to Re, as well as to the ratio of νt and the molecular kinematic viscosity (ν).

42 ENGINEERING↗