Search NASA⌕ Search

SEARCH · Search NASA

Results for “concentration profile fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Ink-sprayed transparent ceramic multilayer planar waveguides for ribbon lasers

Multilayer transparent ceramic planar waveguides were fabricated using a newly developed particle-loaded ink spraying method, characterized, and lased. Waveguides with five Yb:YAG gain regions and four Lu:YAG index-matched regions were produced for application as a “ribbon laser,” as well as a homogeneous Yb:YAG waveguide. Elemental mapping revealed that the waveguide thicknesses varied from 40 to 69 μm, comprised of layers of less than 10 µm. The Yb concentration profiles were fit to a diffusion profile. The homogeneous and multilayer waveguides were tested in a laser cavity and produced similar slope efficiencies, with a maximum of 31% achieved. In conclusion, their beam profiles provided evidence of the ribbon structure affecting the preferred mode.

Osborne, Ross [Lawrence Livermore National Laborat↗

Estimation of Fission Product Transport Parameters for Cesium in the AGR-3/4 TRISO Fuel Experiment

A one-dimensional (1D) finite-element model of fission product transport in the AGR-3/4 experiment has been developed using the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and implemented in the fuel performance code, BISON. The model resolves capsule-specific geometries, materials, and temperature histories and simulates radial migration of fission products from the fuel compact through the inner ring, outer ring, and into the sink ring. Model parameters governing diffusion and sorption were estimated for key fission products – cesium (Cs), and europium (Eu) – by simultaneously fitting modeled isotopic concentration profiles and total ring inventories to a post-irradiation experimental measurement. These data include gamma scanning, liquid scintillation for Sr-90, radial deconsolidation leach-burn-leach analysis, tomographic reconstructions, and destructive physical sampling. A mortar-based interfacial sorption framework was implemented to enforce physically consistent mass transfer and flux conservation across gas gaps. Two classes of parameter sets were derived: a least-squares best-fit, and a safety-oriented conservative-fit, what applies strong penalties for underprediction of sink inventories. Across all twelve capsules, the model successfully reproduces the dominant radial transport trends for Cs, Sr, with decreasing concentrations from the compact outward through successive rings. Cs behavior is captured most consistently, while strontium predictions reveal systematic trade-offs between compact accuracy and conservative sink-ring bounding. The results demonstrate that sink ring weighted calibration provides conservative, safety relevant bounds on low temperature fission product transport, but at the cost of underpredicting compact inventories for Sr isotopes. These discrepancies highlight the need for additional physics, including fast-slow diffusion model, incorporating trapping mechanism in the transport behavior. Overall, this work establishes a robust, capsule-specific modeling framework for AGR-3/4 fission product transport and provides a defensible basis for parameter selection in source-term and fuel performance analyses for high temperature gas-cooled reactors.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

($\mathrm{INVITED}$)Counter-ion effect on the diffusion behavior of $\mathrm{Y}$b, $\mathrm{L}$u, and $\mathrm{N}$d ions in $\mathrm{YAG}$ transparent ceramics

The ability to fabricate additively manufactured laser waveguides with sharp dopant concentration interfaces is limited by diffusion of the dopants at the temperatures required to fully densify the material. Compositional analysis of bilayer samples, where each layer was either undoped YAG or YAG doped with Yb, Lu, or Nd, were fabricated such that all combinations were available for testing. Samples were fabricated at both 1750°C and 1850°C to determine the diffusion behavior of each dopant alone and also in the presence of a second dopant. It was found that the experimental concentration profiles exhibited both intragranular (bulk) and grain boundary contributions, and thus fitting to a complementary error function equation required the use of two diffusion coefficients. Nd always diffused further along grain boundaries than the other dopants, due in part to its small segregation coefficient in YAG. It is shown that the presence of Nd as a counter dopant inhibits intragranular diffusion of other dopants while enhancing their grain boundary diffusion. All the observed trends were attributed to a combination of intragranular lattice strain due to: dopant ions replacing yttrium substitutionally, the relative driving forces for the segregation of dopant ions to grain boundaries, and the ability of one dopant to affect the diffusion of a different dopant in the other direction to maintain charge neutrality.

36 MATERIALS SCIENCE↗

Resolving atomic transport through iron oxide under irradiation using isotopic tracers

Structural materials in nuclear environments are exposed to complex and harsh conditions, including high temperatures, corrosive solutions, and irradiation. Material stability and likewise degradation are fundamentally controlled by atomic transport mechanisms. Here, we investigate the influence of radiation-induced, non-equilibrium point defect populations on mass transport within model oxide systems using isotopic tracers and atom probe tomography. An isotopically labelled hematite (a-Fe2O3) film was synthesized with 16O- and 18O-rich layers (>30nm and 10 nm, respectively) sequentially deposited via molecular beam epitaxy. One region of the film was proton irradiated to ~0.1 dpa at 450°C, and regions outside the beam spot provided a non-irradiated control of thermal annealing only. 3D APT was used to directly visualize the elemental and O-isotopic redistribution following both thermal and combined thermal and irradiation conditions. Complementary scanning transmission electron microscopy (STEM) was also employed to study the microstructural evolution. Error function fits of the 18O-concentration profiles enabled quantification of diffusion coefficients.

isotopic tracers, diffusion, proton irradiation, h↗

Understanding Multi-Ion Transport Mechanisms in Bipolar Membranes

Bipolar membranes (BPMs) have the potential to become critical components in electrochemical devices for a variety of electrolysis and electrosynthesis applications. Because they can operate under large pH gradients, BPMs enable favorable environments for electrocatalysis at the individual electrodes. Critical to the implementation of BPMs in these devices is understanding the kinetics of water dissociation that occurs within the BPM as well as the co- and counter-ion crossover through the BPM, which both present significant obstacles to developing efficient and stable BPM-electrolyzers. In this study, a continuum model of multi-ion transport in a BPM is developed and fit to experimental data. Specifically, concentration profiles are determined for all ionic species, and the importance of a water-dissociation catalyst is demonstrated. The model describes internal concentration polarization and co- and counter-ion crossover in BPMs, determining the mode of transport for ions within the BPM and revealing the significance of salt-ion crossover when operated with pH gradients relevant to electrolysis and electrosynthesis. Lastly, a sensitivity analysis reveals that the performance and lifetime of BPMs can be improved substantially by using of thinner dissociation catalysts, managing water transport, modulating the thickness of the individual layers in the BPM to control salt-ion crossover, and increasing the ion-exchange capacity of the ion-exchange layers in order to amplify the water-dissociation kinetics at the interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Current–Induced Ion Concentration Polarization at a Perfect Ion–Exchange Patch in an Infinite Insulating Wall

This research examines theoretically the ion concentration polarization, ion fluxes, and electrostatic fields near an ion exchange patch in the wall of an electrified fluidic channel. These phenomena are important in related microfluidic ion-preconcentration systems. Under an electric field, counterions enter the ion-exchange patch at one side and leave at the other, with salt depletion occurring near the entrance and accumulation near the exit. Furthermore, the high patch conductivity and the concentration profiles lead to local electric field perturbations that may facilitate preconcentration. This study includes analytical expressions of ion concentrations and electrochemical potentials at small to moderate electric fields, as well as numerical simulations. Additionally, a simple matrix of polynomial coefficients (obtained via fitting of numerical data) enables analytical calculation of the 2-dimensional concentration profiles at all electric fields within the range investigated in the numerical simulations. This is possible because a single dimensionless parameter controls this problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The solubility of titanite in silicate melt determined from growth and dissolution experiments

The solubility of titanite (CaTiSiO 5 ) in Si-rich melts was measured experimentally through growth experiments at 800–1000 °C, 0.5–1.0 GPa, log fO 2 ~ CCO–0.8, t = 72–168 h, and H 2 O = 0 to 4 wt.%, and in dissolution experiments at 925–1300 °C, 0.8 GPa, t = 18–118 h, and H 2 O = 1–10 wt.% in a piston-cylinder apparatus. Run product glasses in growth experiments were homogeneous, and iron loss suppressed ilmenite crystallization. Saturation concentrations in dissolution experiments were estimated by fitting measured diffusion profiles. Titanite solubility increases with increasing temperature and melt composition parameter \(M = {\text{ molar }}\left( {{\text{Na }} + {\text{ K }} + {\text{ 2Ca}}} \right)/\left( {{\text{Al }} \times {\text{ Si}}} \right)\) . Multiple linear regression of glass composition data from growth and dissolution experiments ( n = 29) plus 39 experiments from the LEPR database (Hirschmann et al. 2008 ) yielded the titanite solubility equation (adj. r 2 = 0.95): \(\left( {TiO_{2} } \right)^{melt} \left( {wt.\% } \right) = 0.978 \times M + 0.0048 \times T\left( K \right){-}5.90\) . This model correctly predicted undersaturation in 95% of 2344 experiments from the LEPR database that did not have titanite or rutile. Application to natural rocks yields saturation temperatures that are similar to independent temperature estimates. Finally, this equation should be useful for constraining the temperatures of titanite-saturated melts, for determining whether titanite saturation in magmatic source regions is likely, and for determining when titanite can crystallize and begin to exert an influence on melt geochemistry.

58 GEOSCIENCES↗

Quantifying Northern High Latitude Gross Primary Productivity (GPP) Using Carbonyl Sulfide (OCS)

The northern high latitude (NHL, 40°N to 90°N) is where the second peak region of gross primary productivity (GPP) other than the tropics. The summer NHL GPP is about 80% of the tropical peak, but both regions are still highly uncertain. Carbonyl sulfide (OCS) provides an important proxy for photosynthetic carbon uptake. Here we optimize the OCS plant uptake fluxes across the NHL by fitting atmospheric concentration simulation with the GEOS-CHEM global transport model to the aircraft profiles acquired over Alaska during NASA's Carbon in Arctic Reservoirs Vulnerability Experiment (2012–2015). We use the empirical biome-specific linear relationship between OCS plant uptake flux and GPP to derive the six plant uptake OCS fluxes from different GPP data. Such GPP-based fluxes are used to drive the concentration simulations. We evaluate the simulations against the independent observations at two ground sites of Alaska. The optimized OCS fluxes suggest the NHL plant uptake OCS flux of -247 Gg S year -1 , about 25% stronger than the ensemble mean of the six GPP-based OCS fluxes. GPP-based OCS fluxes systematically underestimate the peak growing season across the NHL, while a subset of models predict early start of season in Alaska, consistent with previous studies of net ecosystem exchange. The OCS optimized GPP of 34 PgC yr -1 for NHL is also about 25% more than the ensembles mean from six GPP data. Further work is needed to fully understand the environmental and biotic drivers and quantify their rate of photosynthetic carbon uptake in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗

The Odd Dark Matter Halos of Isolated Gas-rich Ultradiffuse Galaxies

Abstract We analyze circular velocity profiles of seven ultradiffuse galaxies (UDGs) that are isolated and gas-rich. Assuming that the dark matter halos of these UDGs have a Navarro–Frenk–White (NFW) density profile or a Read density profile (which allows for constant-density cores), the inferred halo concentrations are systematically lower than the cosmological median, even as low as −0.6 dex (about 5 σ away) in some cases. Alternatively, similar fits can be obtained with a density profile that scales roughly as 1/ r 2 for radii larger than a few kiloparsecs. Both solutions require the radius where the halo circular velocity peaks ( R max ) to be much larger than the median expectation. Surprisingly, we find an overabundance of such large- R max halos in the IllustrisTNG dark-matter-only simulations compared to the Gaussian expectation. These halos form late and have higher spins compared to median halos of similar masses. The inner densities of the most extreme among these late-forming halos are higher than their NFW counterparts, leading to a ∼1/ r 2 density profile. However, the two well-resolved UDGs in our sample strongly prefer lower dark matter densities in the center than the simulated ones. Comparing to IllustrisTNG hydrodynamical simulations, we also find a tension in getting both low enough circular velocities and high enough halo mass to accommodate the measurements. Our results indicate that the gas-rich UDGs present a significant challenge for galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗

Data for Catalyst and Reactor Design Considerations for Selective Production of Acids by Oxidative Cleavage of Alkenes with H2O2

Oxidative cleavage of alkenes and unsaturated fatty acids with hydrogen peroxide gives an efficient and sustainable process to obtain mono- and di-acids for polymers and lubricants with fewer safety risks and less environmental impact than processes that utilize ozone or other inorganic oxidizers (e.g., permanganate, dichromate, etc.). Guided by insight into the mechanisms for competing reaction pathways (i.e., epoxidation of alkene on W–(η2-O2) complexes vs. H2O2 decomposition) and the apparent kinetics derived from kinetic experiments, here, we postulate that W-based heterogeneous catalysts can provide high performance and stable operations at low H2O2 concentrations. Semi-batch reactors with continuous introduction of H2O2 solutions offer the means to maintain low H2O2 concentrations while providing sufficient quantities of H2O2 to satisfy the reaction stoichiometry. We derived simple kinetic model equations for the epoxidation, ring-opening, oxidative cleavage, and oxidation steps and fit theses equations to batch experimental data to obtain kinetic parameters. This kinetic model describes the concentration profiles of reactant, oxidant, and products well as shown by agreement with experimental data. Further predictions of the optimal H2O2 feed rate for semi-batch operation utilized by the proposed rate expressions and the reactor design equations suggest that low H2O2 feed rate increases selectivity towards oxidative cleavage products and selective use of H2O2 for oxidative cleavage pathway. Comparisons of oxidative cleavage of 4-octene in batch and semi-batch reactors show that semi-batch reactors with optimized molar feed rates of H2O2 increased oxidative cleavage product selectivities (76% to 99%; with an increase in butyric acid selectivity from 1% to 55%) and H2O2 selectivity (3% to 30%). In addition, semi-batch reaction conditions used avoid H2O2-mediated dissolution of W-atoms from the catalyst. Analysis of these findings suggest that solid oxide catalysts will be effective for continuous oxidative cleavage reactions if deployed within fixed-bed reactors that allow for distributed introduction of reactants and therefore low in situ concentrations of H2O2.

Catalysis↗

Catalyst and reactor design considerations for selective production of acids by oxidative cleavage of alkenes and unsaturated fatty acids with H 2 O 2

Oxidative cleavage of alkenes and unsaturated fatty acids with hydrogen peroxide gives an efficient and sustainable process to obtain mono- and di-acids for polymers and lubricants with fewer safety risks and less environmental impact than processes that utilize ozone or other inorganic oxidizers (e.g., permanganate, dichromate, etc.). Guided by insight into the mechanisms for competing reaction pathways (i.e., epoxidation of alkene on W–(η 2 -O 2 ) complexes vs. H 2 O 2 decomposition) and the apparent kinetics derived from kinetic experiments, here, we postulate that W-based heterogeneous catalysts can provide high performance and stable operations at low H 2 O 2 concentrations. Semi-batch reactors with continuous introduction of H 2 O 2 solutions offer the means to maintain low H 2 O 2 concentrations while providing sufficient quantities of H 2 O 2 to satisfy the reaction stoichiometry. We derived simple kinetic model equations for the epoxidation, ring-opening, oxidative cleavage, and oxidation steps and fit theses equations to batch experimental data to obtain kinetic parameters. Furthermore, this kinetic model describes the concentration profiles of reactant, oxidant, and products well as shown by agreement with experimental data. Further predictions of the optimal H 2 O 2 feed rate for semi-batch operation utilized by the proposed rate expressions and the reactor design equations suggest that low H 2 O 2 feed rate increases selectivity towards oxidative cleavage products and selective use of H 2 O 2 for oxidative cleavage pathway. Comparisons of oxidative cleavage of 4-octene in batch and semi-batch reactors show that semibatch reactors with optimized molar feed rates of H 2 O 2 increased oxidative cleavage product selectivities (76% to 99%; with an increase in butyric acid selectivity from 1% to 55%) and H 2 O 2 selectivity (3% to 30%). In addition, semibatch reaction conditions used avoid H 2 O 2 -mediated dissolution of W-atoms from the catalyst. Analysis of these findings suggest that solid oxide catalysts will be effective for continuous oxidative cleavage reactions if deployed within fixed-bed reactors that allow for distributed introduction of reactants and therefore low in situ concentrations of H 2 O 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fine Root and Soil Organic Carbon Depth Distributions are Inversely Related Across Fertility and Rainfall Gradients in Lowland Tropical Forests

Humid tropical forests contain some of the largest soil organic carbon (SOC) stocks on Earth. Much of this SOC is in subsoil, yet variation in the distribution of SOC through the soil profile remains poorly characterized across tropical forests. We used a correlative approach to quantify relationships among depth distributions of SOC, fine root biomass, nutrients and texture to 1 m depths across 43 lowland tropical forests in Panama. The sites span rainfall and soil fertility gradients, and these are largely uncorrelated for these sites. We used fitted b parameters to characterize depth distributions, where b is a numerical index based on an asymptotic relationship, such that larger b values indicate greater concentrations of root biomass or SOC at depth in the profile. Root b values ranged from 0.82 to 0.95 and were best predicted by soil pH and extractable potassium (K) stocks. For example, the three most acidic (pH < 4) and Kpoor (< 20 g K m -2 ) soils contained 76 ± 5% of fine root biomass from 0 to 10 cm depth, while the three least acidic (pH > 6.0) and most K-rich (> 50 g K m -2 ) soils contained only 41 ± 9% of fine root biomass at this depth. Root b and SOC b values were inversely related, such that a large fine root biomass in surface soils corresponded to large SOC stocks in subsoils (50–100 cm). SOC b values were best predicted by soil pH and base cation stocks, with the three most base-poor soils containing 34 ± 8% of SOC from 50 to 100 cm depth, and the three most base-rich soils containing just 9 ± 2% of SOC at this depth. Nutrient depth distributions were not related to Root b or SOC b values. Here, these data show that large surface root biomass stocks are associated with large subsoil C stocks in strongly weathered tropical soils. Further studies are required to evaluate why this occurs, and whether changes in surface root biomass, as may occur with global change, could in turn influence SOC storage in tropical forest subsoils.

54 ENVIRONMENTAL SCIENCES↗

New Determination of Fundamental Properties of Palomar 5 Using Deep DESI Imaging Data

The legacy imaging surveys for the Dark Energy Spectroscopic Instrument project provides multiple-color photometric data, which are about 2 mag deeper than those from the SDSS. In this study, we redetermine the fundamental properties for an old halo globular cluster of Palomar 5 based on these new imaging data, including structure parameters, stellar population parameters, and luminosity and mass functions. These characteristics, together with its tidal tails, are key for dynamical studies of the cluster and constraining the mass model of the Milky Way. By fitting the King model to the radial surface density profile of Palomar 5, we derive the core radius of [[Formula]], tidal radius of [[Formula]], and concentration parameter of c = 0.78 ± 0.04. We apply a Bayesian analysis method to derive the stellar population properties and get an age of 11.508 ± 0.027 Gyr, metallicity of [Fe/H] = −1.798 ± 0.014, reddening of E(B − V) = 0.0552 ± 0.0005, and distance modulus of (m−M){sub 0} = 16.835±0.006. The main-sequence luminosity and mass functions for both the cluster center and tidal tails are investigated. The luminosity and mass functions at different distances from the cluster center suggest that there is obvious spatial mass segregation. Many faint low-mass stars have been evaporated at the cluster center, and the tidal tails are enhanced by low-mass stars. Both the concentration and relaxation times suggest that Palomar 5 is a totally relaxed system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Substructure at High Speed. II. The Local Escape Velocity and Milky Way Mass with Gaia eDR3

Abstract Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib & Lin, we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the “tail” of the velocity distribution and the consistency of the data with different underlying models. In this paper, we apply this method to the Gaia eDR3 data release and find that a model with two components is preferred, although results from a single-component fit are also consistent. Based on a fit to retrograde data with two bound components to account for the relaxed halo and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be v esc = 445 − 8 + 25 km s −1 . A fit with a single component to the same data gives v esc = 472 − 12 + 17 km s −1 . Assuming a Navarro−Frenck−White dark matter profile, we find a Milky Way concentration of c 200 = 19 − 7 + 11 and a mass of M 200 = 4.6 − 0.8 + 1.5 × 10 11 M ⊙ , which is considerably lighter than previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the halo occupation distribution model for emission line galaxies

We study the modelling of the halo occupation distribution (HOD) for the eBOSS DR16 emission line galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both Navarro–Frenk–White and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, w p , captures most of the effects from the PDF and satellites profile. The quadrupole, ξ 2 , captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. Here, we fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗