Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass diffusivity”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Divergence of Velocity Fields in Electrochemical Systems

The passage of current through electrochemical systems results in the development of concentration gradients in the electrolytic phase that can be modeled using concentrated solution theory. Application of this theory requires knowledge of three concentration-dependent transport coefficients, which are often taken to be conductivity, diffusion coefficient, and the cation transference number with respect to the solvent velocity. The governing diffusion equation for molar concentration contains two additional terms - the thermodynamic factor which is related to activity of the electrolytic species and the solvent velocity. The main advance in this paper is the derivation of an expression for the divergence of the solvent velocity. Solving this equation requires knowledge of the partial molar volume of the electrolyte. Analogous expressions are derived for the mass average and molar average velocities. These velocities occur naturally in the diffusion equation if concentration is expressed as weight fraction and mole fraction of the electrolytic species, respectively.

diffusion↗
DOE OSTI2026

Beyond Mass and Multiscale Environments: What Shapes Low Surface Brightness Galaxies? Evidence from MaNGA

The origin of low surface brightness (LSB) galaxies remains a key open question in galaxy formation, reflecting the balance between internal mechanisms and environmental influence. Using MaNGA integral-field spectroscopy, we investigate whether LSB and high surface brightness (HSB) galaxies of comparable stellar mass (9 < logM * < 10) occupy distinct environments or differ primarily through internal evolution. Our late-type sample comprises 113 central and 29 satellite LSB galaxies, and 374 central and 142 satellite HSB galaxies. We characterize environments on scales from ∼100 kpc to 10 Mpc, analyzing radial profiles of stellar mass surface density (Σ * ), star formation activity, and gas-phase metallicity. Central LSB and HSB galaxies inhabit similarly low-density large-scale (>200 kpc) environments, but LSB galaxies are more isolated on small scales (∼100 kpc). Even after matching in stellar mass and environment, LSB galaxies show systematically lower Σ * , Σ SFR , and metallicities, often hosting diffuse, weakly star-forming bulges embedded in extended disks. These results indicate that LSB structure and star formation are not primarily governed by a large-scale environment or halo mass. While secondary halo properties such as spin, concentration, or gas accretion history are often invoked, their environmental dependence appears weak. Instead, LSB–HSB differences for centrals likely reflect divergent assembly or interaction histories and internal processes—such as angular momentum-driven disk evolution or inefficient gas conversion—largely decoupled from a large-scale environment. Nonetheless, the environment still influences the observed star formation and chemical differences between central and satellite LSB galaxies.

Shen, Mengting [Xiamen University (China); SDSS Co↗
DOE OSTI2022

Enhancing reactant selectivity for Ni/Mg reforming catalysts using silicalite-1 shells: A modeling study

The ability of a zeolite shell to enhance the selective conversion of hydrocarbons through diffusional limitations was investigated using a multi-scale model of a fixed-bed reactor. The impact of shell thickness and molecule/pore size on the catalytic performance of silicalite-1 zeolite encapsulated nickel catalyst pellets for steam reforming of C1-C7 hydrocarbons is reported. A reaction–diffusion model using kinetic expressions established in literature was employed. The model was verified through comparison with reported experimental results for steam reforming data over a temperature range of 748 – 1113 K and pressure of 1 – 10 bar. Comparisons are also made against experimental data for steam reforming in the presence of a zeolite shell. Evaluation of the Weisz-Prater criterion for both the core and encapsulated catalyst confirmed mass transfer limitation induced by the utilization of a zeolite shell. The model was used to suggest an optimal thickness that balances diffusional limitations imposed by the zeolite layer on methane versus that of the heavier hydrocarbons. The optimum thickness varied as a function of hydrocarbon size and shape which determined the diffusion rates. For toluene and heptane, a 50 nm thick shell was sufficient to wholly prevent reaction. Hydrocarbons like propane and butane required a shell 7.5 and 5 μm thick. Increasing the gas-hourly-space-velocity from 10,000 to 60,000 h –1 caused a decrease in the optimum shell thickness. Furthermore, this approach can be modified for application to other mixed hydrocarbon systems to predict optimal catalyst design.

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

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{&#x221E;}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Elucidating the Transition of 3D Morphological Evolution of Binary Alloys in Molten Salts with Metal Ion Additives

Molten salts serve as effective high-temperature heat transfer fluids and thermal storage media used in a wide range of energy generation and storage facilities, including concentrated solar power plants, molten salt reactors and high-temperature batteries. However, at the salt–metal interfaces, a complex interplay of charge-transfer reactions involving various metal ions, generated either as fission products or through corrosion of structural materials, takes place. Simultaneously, there is a mass transport of ions or atoms within the molten salt and the parent alloys. The precise physical and chemical mechanisms leading to the diverse morphological changes in these materials remain unclear. Here, to address this knowledge gap, this work employed a combination of synchrotron X-ray nanotomography and electron microscopy to study the morphological and chemical evolution of Ni-20Cr in molten KCl-MgCl 2 , while considering the influence of metal ions (Ni 2+ , Ce 3+ , and Eu 3+ ) and variations in salt composition. Our research suggests that the interplay between interfacial diffusivity and reactivity determines the morphological evolution. The summary of the associated mass transport and reaction processes presented in this work is a step forward toward achieving a fundamental comprehension of the interactions between molten salts and alloys. Overall, the findings offer valuable insights for predicting the diverse chemical and structural alterations experienced by alloys in molten salt environments, thus aiding in the development of protective strategies for future applications involving molten salts.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

A horseshoe-shaped ring of diffuse emission detected at 1.4 GHz

We identify a horseshoe-shaped ring (HSR) of diffuse emission in J1407+0453 from the Faint Images of Radio Sky at Twenty-cm (FIRST) survey using the Very Large Array telescope at 1.4 GHz. An optical galaxy, SDSS J140709.01+045302.1, is present near the limb of the HSR of J1407+0453, with a spectroscopic redshift ofz = 0.13360. The total extent of the source, including the diffuse emission of J1407+0453, is 65 arcsec (with a physical extent of 160 kpc), whereas the diameter of the HSR is approximately 10 arcsec (25 kpc). The flux density of the HSR is ∼47 mJy at 1400 MHz, whereas the flux densities of the whole diffuse emission of J1407+0453 at 1400 MHz and 150 MHz are 172 mJy and 763 mJy, respectively. We measure the radio luminosity of HSR J1407+0453 as 1.94 × 10 24 W Hz −1 , with a spectral index,α 150 1400 = −0.67. The black hole mass of J1407+0453 is 5.8 × 10 8 M ⊙ . We compare the HSR of diffuse emission of J1407+0453 with other discovered diffuse circular sources. The possible formation scenarios for J1407+0453 are discussed, so as to understand the nature of the source. We present a spectral index map of source J147+0453 to study the spectral properties of the source.

Astronomy & Astrophysics↗
DOE OSTI2021

Adapting confocal Raman microscopy for in situ studies of redox transformations at electrode-electrolyte interfaces

Confocal Raman microscopy was applied to quantify redox species present within the diffusion layer adjacent to an electrode surface under potentiostatic control. A glass microscope coverslip with a thin indium tin oxide (ITO) coating served as both the working electrode and optical window for a microscope-stage mountable spectroelectrochemical cell. A high numerical aperture objective mounted in an inverted microscope frame just below the stage brought excitation radiation through the coverslip window and to a tight focus a few micrometers above the ITO film surface. Species diffusing into the confocal probe volume defined by the excitation beam focus and the collected light region were detected, identified and quantified based on their Raman scattering frequencies and intensities. Here, in measurements that interrogated the interconversion of ferrocyanide and ferricyanide ions as a function of applied voltage, least-squares regression analysis of spectral datasets predicted the formal potential and relative surface concentrations of the ions in good agreement with the expected Nernstian response. Preliminary studies of methyl viologen reduction at an ITO film/Nafion membrane interface were conducted and show the possibility for estimation of mass transport coefficients of redox species within ionic polymer materials.

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

Random Walks and Sticky Surfaces: Single-Molecule Measurements of Solute Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Nanoporous anodic aluminum oxide (AAO) membranes are now being explored for use in advanced chemical separations, including in the dehydration of biofuels such as ethanol. Optimization of membrane performance requires an in-depth understanding of how solvent mixtures and solutes behave under nanoconfinement. In this work, the diffusion of rhodamine B (RhB) dye through 10 nm and 20 nm AAO nanopores filled with a series of ethanol/water mixtures (0 - 33% water) is explored by fluorescence correlation spectroscopy (FCS). RhB was found to diffuse through the pores by two distinct mechanisms with mean diffusion coefficients, $D_f$ and $D_s$, reflecting fast and slow diffusive motions, respectively, with values that differ by nearly two orders of magnitude. Further, both $D_f$ and $D_s$ increased with pore size and were significantly smaller than $D_b$, the RhB diffusion coefficient in bulk liquid. Mean $D_f$ values follow a composition-dependent trend that closely mimics the viscosity dependence of $D_b$. Additional slowing of fast RhB diffusion is attributed to both hydrodynamic drag and electrostatic interactions with the nanopore surface. The mean $D_s$ values exhibit a different trend with increasing water content, revealing an increase in $D_s$ and a decrease in the contributions of slow diffusion to the observed dynamics. The fluorescence time transient data used in the analysis show that the slow diffusion process is strongly hindered and likely involves frequent adsorption of RhB to the pore surfaces. These results provide new insights into the detailed molecular-level mechanisms of mass transport in nanoporous AAO membranes.

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

Methane mass transfer in mesoporous silica saturated with liquid hydrocarbons

Mass transfer across gas/liquid interfaces plays a central role in many industrial applications. In particular, gas dissolution and diffusion in liquid hydrocarbon mixtures, confined in nanometer-sized pores, is an essential mechanism during enhanced oil recovery (EOR) from unconventional formations. In this work, we have measured methane (C 1 ) diffusion in n-decane (C 10 ), n-hexadecane (C 16 ), and mixtures of C 10 + C 16 in a mesoporous material with an average pore size of 4 nm at 50 °C and ~ 8 MPa of gas pressure. A key conclusion of this work is that the diffusivities, measured in the bulk phase, for the relevant binary systems, are sufficient to predict the diffusion behavior of the corresponding multicomponent systems in the porous medium by using Wilke’s equation for the evaluation of effective component diffusivities, combined with an accurate equation of state (EOS) representation of the (bulk) phase behavior. This observation can facilitate accurate prediction of recovery processes in unconventional formations and, potentially, guide other applications that entail gas-liquid interface mass transfer in mesoporous materials.

42 ENGINEERING↗
DOE OSTI2020

Mesoporous silica-encapsulated gold core–shell nanoparticles for active solvent-free benzyl alcohol oxidation

Silica-encapsulated gold core@shell nanoparticles (Au@SiO2 CSNPs) were synthesized via a bottom-up procedure and used to catalyze the selective oxidation of benzyl alcohol. The pore size, morphology, crystallinity and composition of Au@SiO2 was evaluated using non-local density functional theory, transmission electron microscopy, high-energy x-ray diffraction and inductively coupled plasma-mass spectroscopy, respectively. The nanoparticles exhibit a mesoporous shell with an average thickness of 25.5 Å which can enhance selectivity via preferential transport of the desired product (i.e., benzaldehyde) relative to larger, undesired products (i.e, benzoic acid/benzyl benzoate). GC-FID analysis revealed the addition of potassium carbonate to the solvent-free oxidation of benzyl alcohol increased conversion from 17.3 to 60.4% while decreasing selectivity from 98.7 to 75.0%. Under equivalent conditions, a bare gold nanoparticle control catalyst deposited on a silica support with a similar gold surface area took 6 times as long to reach the same conversion, achieving only 49.4% selectivity. These results suggest that the pore size distribution within the inert silica shell of Au@SiO2 CSNPs inhibits the formation of undesired products to facilitate the selective oxidation of benzaldehyde despite a basic environment, which reduces selectivity under typical conditions. The CSNPs demonstrated a much lower activation energy than the Au-SiO2 control catalyst, 37 ± 1.9 kJ/mol and 72 ± 7.1 kJ/mol, respectively. Thiele modulus analysis indicates the CSNP pore structure does not create a mass transport limitation due to the nano-scale path of diffusion through the pore structure to the active surface. The lower activation energy and mesopore distribution together suggest the Au@SiO2 catalyst demonstrates higher activity through beneficial in-pore orientation, both reducing competitive adsorption and promoting a single, lower activation energy mechanistic pathway.

benzyl alcohol oxidation, catalysis, confinement↗
DOE OSTI2024

DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations

DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. Here, the flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Interacting Dark Matter and the Origin of Ultradiffuse Galaxies NGC1052-DF2 and -DF4

Observations of ultradiffuse galaxies NGC 1052-DF2 and -DF4 show they may contain little dark matter, challenging our understanding of galaxy formation. Using controlled N-body simulations, we explore the possibility that their properties can be reproduced through tidal stripping from the elliptical galaxy NGC 1052, in both cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. To explain the dark matter deficiency, we find that a CDM halo must have a very low concentration so that it can lose sufficient inner mass in the tidal field. In contrast, SIDM favors a higher and more reasonable concentration as core formation enhances tidal mass loss. Final stellar distributions in our SIDM benchmarks are more diffuse than the CDM one, and hence the former are in better agreement with the data. We further show that a cored CDM halo model modified by strong baryonic feedback is unlikely to reproduce the observations. Our results indicate that SIDM is more favorable for the formation of dark-matter-deficient galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. In conclusion, our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Degradation Mechanism Due to Water Ingress Effect on the Top Contact of Cu(In,Ga)Se2 Solar Cells

The impact of moisture ingress on the surface of copper indium gallium diselenide (CIGS) solar cells was studied. While industry-scale modules are encapsulated in specialized polymers and glass, over time, the glass can break and the encapsulant can degrade. During such conditions, water can potentially degrade the interior layers and decrease performance. The first layer the water will come in contact with is the transparent conductive oxide (TCO) layer. To simulate the impact of this moisture ingress, complete devices were immersed in deionized water. To identify the potential sources of degradation, a common window layer for CIGS devices—a bilayer of intrinsic zinc oxide (i-ZnO) and conductive indium tin oxide (ITO)—was deposited. The thin films were then analyzed both pre and post water soaking. To determine the extent of ingress, dynamic secondary ion mass spectroscopy (SIMS) was performed on completed devices to analyze impurity diffusion (predominantly sodium and potassium) in the devices. The results were compared to device measurements, and indicated a degradation of device efficiency (mostly fill factor, contrary to previous studies), potentially due to a modification of the alkali profile.

14 SOLAR ENERGY↗
DOE OSTI2021

The Population of H ii Regions in NGC 7020

The galaxy NGC 7020 displays an exotic hexagonal ringlike central structure with conspicuous ansae located at two opposite vertices and a tenuous external ring populated by H ii regions. Inside and around the hexagonal structure, Hα emission is also present at the inner disk. To characterize the population of the H ii regions, as well as their ionizing clusters, we imaged NGC 7020 with narrowband Hα and nearby continuum filters attached to GMOS-S installed on the Gemini South telescope. We found 202 H ii regions or complexes of H ii regions evenly distributed between the outer ring and the central disk The nucleus and ansae also present Hα emission. The equivalent width of the Hα line (W {sub Hα}) is systematically greater at the regions of the outer ring relative to those of the inner disk. We discuss the influence of the metallicity gradient of the disk and the upper limit of the masses of the initial mass function on W {sub Hα}, and we conclude that the data are still consistent with the occurrence of a younger burst of H ii region formation in the outer ring. The central regions present more massive clusters, M ≥ 10{sup 6} M {sub ⊙}, than those of the outer ring (M ≤ 10{sup 6} M {sub ⊙}). Three clusters within 5″ of the nucleus present masses higher than 10{sup 8} M {sub ⊙}. The presence of diffuse Hα emission in the inner 5″ suggests gas flows in the nuclear region.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Heterogeneous oxidation of amorphous organic aerosol surrogates by O 3 , NO 3 , and OH at typical tropospheric temperatures

Typical tropospheric temperatures render possible phase states of amorphous organic aerosol (OA) particles of solid, semisolid, and liquid. This will affect the multiphase oxidation kinetics involving the organic condensed-phase and gaseous oxidants and radicals. To quantify this effect, we determined the reactive uptake coefficients (γ) of O 3 , NO 3 , and OH by substrate films composed of single and binary OA surrogate species under dry conditions for temperatures from 213 to 313 K. A temperature-controlled coated-wall flow reactor coupled to a chemical ionization mass spectrometer was applied to determine γ with consideration of gas diffusion transport limitation and gas flow entrance effects, which can impact heterogeneous reaction kinetics. The phase state of the organic substrates was probed via the poke-flow technique, allowing the estimation of the substrates' glass transition temperatures. γ values for O 3 and OH uptake to a canola oil substrate, NO 3 uptake to a levoglucosan and a levoglucosan / xylitol substrate, and OH uptake to a glucose and glucose / 1,2,6-hexanetriol substrate have been determined as a function of temperature. We observed the greatest changes in γ with temperature for substrates that experienced the largest changes in viscosity as a result of a solid-to-liquid phase transition. Organic substrates that maintain a semisolid or solid phase state and as such a relatively higher viscosity do not display large variations in heterogeneous reactivity. From 213 to 293 K, γ values of O 3 with canola oil, of NO 3 with a levoglucosan / xylitol mixture, and of OH with a glucose / 1,2,6-hexanetriol mixture and canola oil, increase by about a factor of 34, 3, 2, and 5, respectively, due to a solid-to-liquid phase transition of the substrate. These results demonstrate that the surface and bulk lifetime of the OA surrogate species can significantly increase due to the slowed heterogeneous kinetics when OA species are solid or highly viscous in the middle and upper troposphere. This experimental study will further our understanding of the chemical evolution of OA particles with subsequent important consequences for source apportionment, air quality, and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Summary of Graphite Data Stored within NDMAS

The Graphite Technology Development Project provides data to support the design of graphite core components within specific reactor service conditions of the next generation of high-temperature, gas-cooled nuclear reactors. Physical, mechanical, and thermal properties of nuclear grade graphite were characterized for specimens that were unirradiated, irradiated, and irradiated under various stress conditions. The material properties include diameter, length, mass, density, compressive strength, tensile strength, flexural strength, modulus, resistivity, thermal diffusivity, and thermal expansion coefficient. Baseline graphite specimens are unirradiated from different grades (2114, IG-110, NBG-17, NBG-18, and PCEA) and different types used in different characterization tests (compressive, flexural, tensile, one-inch cylinder, and quarter-inch cylinder). The Advanced Graphite Creep (AGC) irradiation specimens are from a much larger number of grades and cylinder types (creep, piggyback, and pencil). For baseline graphite, characterization data for 7,756 specimens extracted from thirteen graphite billets were captured to the NDMAS database. For AGC experiments, four irradiation campaigns have been completed: AGC-1, AGC-2, AGC-3, and AGC-4. The ongoing HDG-1 (High-Dose Graphite) experiment, which began irradiation with Cycle 168B on August 26, 2020, includes specimens previously irradiated in AGC-2 in addition to the specimens originally destined for AGC-5.. Besides the characterization data, the AGC data includes irradiation monitoring and physics data representing the irradiation conditions of AGC specimens. Currently, all data for AGC-1, AGC-2, and AGC-3 have been captured to the NDMAS database. Only AGC-4 pre-irradiation and irradiation monitoring data have been captured, and HDG-1 pre-irradiation data are in process of being captured for unirradiated specimens. To date, a total of 38,149 material property records have been captured into NDMAS database for baseline and AGC specimens. All characterization data are qualified for use according to their perspective data verification reports. For the AGC irradiation campaigns, a total of 41,502 qualified physics calculation records were added for AGC-1, AGC-2, and AGC-3. Finally, a total of 173,698,505 AGC irradiation monitoring data records (thermocouple temperature, gas flow rate, gas pressure, gas moisture, applied load, and specimen displacement) have been captured to NDMAS; the majority of those records (~94%) are qualified data records.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Milky Way mass with K giants and BHB stars using LAMOST, SDSS/SEGUE, and Gaia : 3D spherical Jeans equation and tracer mass estimator

ABSTRACT We measure the enclosed Milky Way mass profile to Galactocentric distances of ∼70 and ∼50 kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Sloan Digital Sky Survey/Sloan Extension for Galactic Understanding and Exploration (SDSS/SEGUE) K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and Gaia DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to ∼70 and ∼50 kpc, the Jeans method yields total enclosed masses of 4.3 ± 0.95 (random) ±0.6 (systematic) × 1011 M⊙ and 4.1 ± 1.2 (random) ±0.6 (systematic) × 1011 M⊙ for the KG and BHB stars, respectively. For the KG and BHB samples, we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) ±0.083 (systematic) × 1012 M⊙ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) ±0.15 (systematic) × 1012 M⊙, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗