Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spherical 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 487 records · Page 27

The ellipticity parameterization for an NFW profile: An overlooked angular structure in strong lens modeling

Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is a Navarro Frenk White (NFW) profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because an unmodeled azimuthal structure has been shown to be able to bias lens model results, we explored the degree to which this azimuthal structure that was introduced can affect the model accuracy. We constructed input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculated the value of the Hubble parameter H 0 one would determine from the lensing fit. We found that the fits to the NFWp input return H 0 values that are systematically biased by about 3% lower than the NFWm counterparts. We explored whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, ξ 2 . We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias was found beyond the MST due to the azimuthal structure exterior to the Einstein radius. We recommend modelers use an NFWm prescription in the future, such that the azimuthal structure can be introduced explicitly rather than implicitly.

79 ASTRONOMY AND ASTROPHYSICS↗

First-principles thermal conductivity of warm-dense deuterium plasmas for inertial confinement fusion applications

Thermal conductivity (κ) of both the ablator materials and deuterium-tritium (DT) fuel plays an important role in understanding and designing inertial confinement fusion (ICF) implosions. The extensively used Spitzer model for thermal conduction in ideal plasmas breaks down for high-density, low-temperature shells that are compressed by shocks and spherical convergence in imploding targets. A variety of thermal-conductivity models have been proposed for ICF hydrodynamic simulations of such coupled and degenerate plasmas. The accuracy of these κ models for DT plasmas has recently been tested against first-principles calculations using the quantum molecular-dynamics (QMD) method; although mainly for high densities (ρ >100 g/cm^3), large discrepancies in κ have been identified for the peak-compression conditions in ICF. To cover the wide range of density-temperature conditions undergone by ICF imploding fuel shells, we have performed QMD calculations of κ for a variety of deuterium densities of ρ = 1.0 to 673.518 g/cm^3, at temperatures varying from T = 5 × 10^3 K to T = 8 × 10^6 K. The resulting κQMD of deuterium is fitted with a polynomial function of the coupling and degeneracy parameters and θ, which can then be used in hydrodynamic simulation codes. Compared with the “hybrid” Spitzer-Lee-More model currently adopted in our hydrocode LILAC, the hydrosimulations using the fitted κQMD have shown up to 20% variations in predicting target performance for different ICF implosions on OMEGA and direct-drive–ignition designs for the National Ignition Facility (NIF). The lower the adiabat of an imploding shell, the more variations in predicting target performance using κQMD. Moreover, the use of κQMD also modifies the shock conditions and the density-temperature profiles of the imploding shell at early implosion stage, which predominantly affects the final target performance. This is in contrast to the previous speculation that κQMD changes mainly the inside ablation process during the hot-spot formation of an ICF implosion.

Hu, S. X.↗

The Thomas theorem and the Efimov States within a generalized Bohr model

Using a generalized Bohr model and the hyper-spherical formalism for a three-body system, we derive the universal energy functions. We apply our model to 4He atom where only the Coulomb potential is dominant and find that its binding energy is well reproduced (less than 5.5% difference). Later, we focus on the equal mass three-body systems and derive the Thomas theorem assuming a simple interaction depending on the range of the potential. We discuss the conditions for which an unbound two-body system produces a bound three-body system and apply our model to 4He and triton atoms as well as to the triton nucleus. Using their scattering lengths and effective ranges, we are able to reproduce the two-body or the three-body binding energies (less than 5% difference) with only one parameter fitted. Prediction for excited (Efimov) levels are also given and in particular we demonstrate that for some hyper-angles two equal minima appear which indicate a phase (shape) transition similar to the Landau's theory of phase transition. We suggest that the observed excited levels in two different experiments for the triton nucleus are indeed Efimov levels and there may be more surprises.

74 ATOMIC AND MOLECULAR PHYSICS↗

HST-STIS Spectra of Saturn's Rings and Implications for Their Reddening Agent

We obtained HST-STIS spectra of Saturn's main rings in May 2011, using the G230L (and G430L) gratings, with final averaged radial resolution of 160 (and 330) km/pixel. The dataset filled a previous 200-330nm "spectral gap" between Cassini and ground-based spectra. The data provide radial profiles as a function of wavelength, but our most basic product at this point is a set of very low-noise spectra, radially averaged over broad regions of the rings (A, B, C, and Cassini Division). The raw spectra required special processing to remove artifacts due to extended-source grating scatter. We have modeled the spectra using a new particle surface model, which corrects for on-surface shadowing due to the likely very rough ring particle surfaces, and avoids overestimation of intra-mixed "neutral absorber". We correct for non-classical layer effects and finite ring optical depth, and relate our observed reflectivities to the spherical albedos of individual smooth particles. We model these smooth particle albedos using standard Hapke theory for regolith grain mixtures that are either homogeneous and "intramixed" (nonicy absorbers dispersed in water ice regolith grains) or heterogeneous "intimate" mixtures. As candidates for the nonicy contaminants we have considered amorphous carbon, aromatic-rich and aliphatic-rich organic tholins, silicates, hematite and iron metal. For the A and B rings, we find that iron metal (including a new theoretical estimate of the refractive indices of nanometer-sized grains of iron) is not spectrally steep enough in the 200-300nm range, and that aliphatic-rich tholins are either too steep at short wavelengths or too flat at long wavelengths. However, less than 1% by mass of aromatic-rich tholins provides a very good fit across the entire spectral range with no gratuitous "neutral absorber" needed, and a minimum of additional free parameters. The best fits require forward-scattering regolith grains. For the C Ring and Cassini Division, additional absorbers are needed (updated results will be given).

Saturn's Rings↗

Interstellar Dust Models

A viable interstellar dust model - characterized by the composition, morphology, and size distribution of the dust grains and by the abundance of the different elements locked up in the dust - should fit all observational constraints arising primarily from the interactions of the dust with incident radiation or the ambient gas. As a minimum, these should include the average interstellar extinction, the infrared emission from the diffuse interstellar medium (ISM), and the observed interstellar abundances of the various refractory elements. The last constraint has been largely ignored, resulting in dust models that require more elements to be in the dust phase than available in the ISM. In this talk I will describe the most recent advances towards the construction of a comprehensive dust model made by Zubko, Dwek, and Arendt, who, for the first time, included the interstellar abundances as explicit constraints in the construction of interstellar dust models. The results showed the existence of many distinct models that satisfy the basic set of observational constraints, including bare spherical silicate and graphite particles, PAHs, as well as spherical composite particles containing silicate, organic refractories, water ice, and voids. Recently, a new interstellar dust constituent has emerged, consisting of metallic needles. These needles constitute a very small fraction of the interstellar dust abundance, and their existence is primarily manifested in the 4 to 8 micron wavelength region, where they dominate the interstellar extinction. Preliminary studies show that these models may be distinguished by their X-ray halos, which are produced primarily by small angle scattering off large dust particles along the line of sight to bright X-ray sources, and probe dust properties largely inaccessible at other wavelengths.

Dwek, Eli↗

The Far-Infrared Emission Line and Continuum Spectrum of the Seyfert Galaxy NGC 1068

We report on the analysis of the first complete far-infrared spectrum (43-197 microns) of the Seyfert 2 galaxy NGC 1068 as observed with the Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory (ISO). In addition to the 7 expected ionic fine structure emission lines, the OH rotational lines at 79, 119 and 163 microns were all detected in emission, which is unique among galaxies with full LWS spectra, where the 119 micron line, where detected, is always in absorption. The observed line intensities were modelled together with IS0 Short Wavelength Spectrometer (SWS) and optical and ultraviolet line intensities from the literature, considering two independent emission components: the AGN component and the starburst component in the circumnuclear ring of approximately 3kpc in size. Using the UV to mid-IR emission line spectrum to constrain the nuclear ionizing continuum, we have confirmed previous results: a canonical power-law ionizing spectrum is a poorer fit than one with a deep absorption trough, while the presence of a big blue bump is ruled out. Based on the instantaneous starburst age of 5 Myr constrained by the Br gamma equivalent width in the starburst ring, and starburst synthesis models of the mid- and far-infrared fine-structure line emission, a low ionization parameter (U=10(exp -3.5)) and low densities (n=100 cm (exp -3)) are derived. Combining the AGN and starburst components, we succeed in modeling the overall UV to far-IR atomic spectrum of SGC 1068, reproducing the line fluxes to within a factor 2.0 on average with a standard deviation of 1.4. The OH 119 micron emission indicates that the line is collisionally excited, and arises in a warm and dense region. The OH emission has been modeled using spherically symmetric, non-local, non-LTE radiative transfer models. The models indicate that the bulk of the emission arises from the nuclear region, although some extended contribution from the starburst is not ruled out. The OH abundance in the nuclear region is expected to be approximately 10(exp -5), characteristic of X-ray dominated regions.

Spinoglio, Luigi↗

Using Cassini CIRS Data to Constrain Enceladus' Libration State

Given the non-spherical shape of Enceladus, the satellite may experience gravitational torques that will cause it to physically librate as it orbits Saturn. Physical libration would produce a diurnal oscillation in the longitude of Enceladus' tidal bulge, which could have a profound effect on the diurnal stresses experienced by the surface of the satellite. Although Cassini ISS has placed an observational upper limit on Enceladus' libration amplitude, small amplitude librations may have geologically significant consequences. For example, a physical libration will affect heat production along the tiger stripes as produced by tidal shear heating. We have modeled the expected power emitted along the tiger stripes for various types of physical libration and have quantified which types of physical libration best reproduce the observed power f1ux as detailed in Cassini CIRS data. We find that including a physical libration does allow better fits to the observations and we have identified regions of the libration phase space that where these fits are optimized. A physical libration has important implications for tidal dissipation within Ence1adus and if identified may provide an additional constraint on its interior mass distribution.

Hurford, T. A.↗

Chain trajectories, domain shapes, and terminal boundaries in block copolymers

The packing geometry of macromolecules in complex mesophases is of key importance to self-organization in synthetic and biological soft materials. While approximate or heuristic models rely on often-untested assumptions about how flexible molecules “fit in” to distinct locations of complex assemblies, physical assemblies derive from ensembles of fluctuating conformations, obscuring the connection between mesophase geometry and the underlying arrangements. Here, we present an approach to extract and analyze features of molecular packing in diblock block copolymer (BCP) melts, a prototypical soft matter system, based on the statistical description of chain conformations in self-consistent field (SCF) theory. We show how average BCP chain trajectories in ordered morphologies can be analyzed from the SCF-derived orientational order parameter of chain segments. We use these extracted trajectories to analyze the features of local packing geometry, including chain bending and tilt, as well as the terminal boundaries that delineate distinct domains in ordered BCP morphologies. In conclusion, we illustrate this analysis by focusing on measurable features of packing frustration in 2D (columnar) and 3D (spherical and bicontinuous) morphologies, notably establishing an explicit link between chain conformations in complex morphologies and their medial geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite element formulation of fluctuating hydrodynamics for fluids filled with rigid particles using boundary fitted meshes

In this paper, we present a finite element implementation of fluctuating hydrodynamics with a moving boundary fitted mesh for treating the suspended particles. The thermal fluctuations are incorporated into the continuum equations using the Landau and Lifshitz approach [1]. The proposed implementation fulfills the fluctuation–dissipation theorem exactly at the discrete level. Since we restrict the equations to the creeping flow case, this takes the form of a relation between the diffusion coefficient matrix and friction matrix both at the particle and nodal level of the finite elements. Brownian motion of arbitrarily shaped particles in complex confinements can be considered within the present formulation. A multi-step time integration scheme is developed to correctly capture the drift term required in the stochastic differential equation (SDE) describing the evolution of the positions of the particles. The proposed approach is validated by simulating the Brownian motion of a sphere between two parallel plates and the motion of a spherical particle in a cylindrical cavity. The time integration algorithm and the fluctuating hydrodynamics implementation are then applied to study the diffusion and the equilibrium probability distribution of a confined circle under an external harmonic potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using Cassini CIRS Data to Constrain Enceladus' Libration State

Given the non-spherical shape of Enceladus, the satellite may experience gravitational torques that will cause it to physically librate as it orbits Saturn. Physical Libration would produce a diurnal oscillation in the longitude of Enceladus' tidal bulge, which could have a profound effect on the diurnal stresses experienced by the surface of the satellite. Although Cassini ISS has placed an observational upper limit on Enceladus' libration amplitude, small amplitude librations may have geologically significant consequences. For example, a physical libration will affect heat production. along the tiger stripes as produced by tidal shear heating. We have modeled the expected power en-litted along the tiger stripes for various types of physical libration and have quantified which types of physical libration best reproduce the observed power flux as detailed in Cassini CIRS data. We find that including a physical libration does allow better fits to the observations and we have identified regions of the libration phase space that where these fits are optimized. A physical libration has important implications for tidal dissipation within Enceladus and if identified may provide an additional constraint on its interior mass distribution.

Hurford, T. A.↗

Toward a NASA Deep Space Optical Communications System

As discussed at SpaceOps in 2016, we expect the data rates from deep space missions to increase approximately one order of magnitude per decade for the next 50 years. The first order of magnitude improvement will come from existing plans for radio frequency (RF) communications including enhancements to both spacecraft and Deep Space Network (DSN) facilities. The next two orders of magnitude are predicted to come from the introduction of deep space optical communications. Studies indicate that optical receive apertures of between 8m12m are desired. The large cost of dedicated receive telescopes makes this method unrealistic – at least in the near-term. The cost of large optical ground terminals is driven primarily by the cost of the optics and by the cost of a stable structure for the telescope. We propose a novel hybrid design in which existing DSN 34m beam waveguide (BWG) radio antennas can be modified to include an 8m equivalent optical primary. By utilizing a low-cost segmented spherical mirror optical design, pioneered by the optical astronomical community, and by exploiting the already existing extremely stable large radio aperture structures in the DSN, we can minimize both of these cost drivers for implementing large optical communications ground terminals. Two collocated hybrid RF/optical antennas could be arrayed to synthesize the performance of an 11.3m receive aperture to support more capable or more distant space missions or used separately to communicate with two optical spacecraft simultaneously. NASA is in the midst of building six new 34m BWG antennas in the DSN. The final two are planned to be built at the DSN Goldstone, California and Canberra complexes. We are now investigating building these last two antennas as RF/optical hybrids. By delaying their operational dates by two years, we would be able to add the 8m optical receive capability for these two antennas while fitting within existing budgetary constraints. This paper describes the hybrid antenna design, the technical challenges being addressed, and plan for using this concept, together with ongoing work on optical flight terminals, to infuse operation optical communications into deep space missions.

Cornwell, Donald M.↗

High-degree Gravity Models from GRAIL Primary Mission Data

We have analyzed Ka‒band range rate (KBRR) and Deep Space Network (DSN) data from the Gravity Recovery and Interior Laboratory (GRAIL) primary mission (1 March to 29 May 2012) to derive gravity models of the Moon to degree 420, 540, and 660 in spherical harmonics. For these models, GRGM420A, GRGM540A, and GRGM660PRIM, a Kaula constraint was applied only beyond degree 330. Variance‒component estimation (VCE) was used to adjust the a priori weights and obtain a calibrated error covariance. The global root‒mean‒square error in the gravity anomalies computed from the error covariance to 320×320 is 0.77 mGal, compared to 29.0 mGal with the pre‒GRAIL model derived with the SELENE mission data, SGM150J, only to 140×140. The global correlations with the Lunar Orbiter Laser Altimeter‒derived topography are larger than 0.985 between l = 120 and 330. The free‒air gravity anomalies, especially over the lunar farside, display a dramatic increase in detail compared to the pre‒GRAIL models (SGM150J and LP150Q) and, through degree 320, are free of the orbit‒track‒related artifacts present in the earlier models. For GRAIL, we obtain an a posteriori fit to the S‒band DSN data of 0.13 mm/s. The a posteriori fits to the KBRR data range from 0.08 to 1.5 micrometers/s for GRGM420A and from 0.03 to 0.06 micrometers/s for GRGM660PRIM. Using the GRAIL data, we obtain solutions for the degree 2 Love numbers, k20=0.024615+/-0.0000914, k21=0.023915+/-0.0000132, and k22=0.024852+/-0.0000167, and a preliminary solution for the k30 Love number of k30=0.00734+/-0.0015, where the Love number error sigmas are those obtained with VCE.

Primary 1↗

Crustal Magnetic Spectra from Correlated Sources on Mars

The spectral method for distinguishing crustal from core-source magnetic fields has been re-examined, modified and applied to both a comprehensive geomagnetic field model and an altitude normalized magnetic map of Mars. These observational spectra are fairly fitted by theoretical forms expected from certain elementary classes of magnetic sources. For Earth we found fields from a core radius 3512 +/- 64 km and a crust represented by a shell of random dipolar sources at radius 6367 +/- 14 km. For Mars we found only a field from a crust represented in the same way, but 46 +/- 10 km below the planetary menu radius. More realistic theoretical spectra, allowing for crustal thickness, oblateness and magnetization by a planet centered dipole, were derived and discussed, as were spectral effects of laterally correlated sources. The main effect of laterally correlated sources is to soften the spectrum at high degrees. We tend to over-estimate source shell depth when this is omitted. To include this effect simply size and magnetization distribution functions for extended sources are recast as a characteristic diameter and mean square magnetization amplitude for an ensemble of vertically magnetized spherical caps on a shell. For small caps, and at moderated degrees, the practical derivatives of the log-theoretical spectrum with respect to amplitude, shell radius, and cap diameter are approximately proportional to l, n, and -n **2, respectively. Separation of diameter from amplitude from layer thickness. Results from applications to observational spectra are discussed, noting that there are now several fine field models for Mars; moreover, the terrestrial magnetic spectrum at high degrees, as revised, updated with high precision Oerstad data, and upgrade with high resolution Champ data, appears softer than before.

Voorhies, Coerte V.↗

Strong evidence that the galactic bulge is shining in gamma rays

There is growing evidence that the Galactic Center Excess identified in the Fermi-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the Fermi bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the GALPROP team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an admixture formation scenario, in which some of the relevant binaries are primordial and the rest are formed dynamically, is preferred over a primordial-only formation scenario at 7.6σ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site.

79 ASTRONOMY AND ASTROPHYSICS↗

SuperCIC: Enhanced Winding Current Density for Hybrid Windings of Tokamaks

The design is presented for a new approach to cable-in-conduit (SuperCIC) for use in the high-field windings of tokamaks. Two layers of high-field superconductor wires are cabled onto a thin-wall perforated center tube. An overwrap is applied and the cable is inserted as a loose fit into a sheath tube. The sheath tube is drawn down onto the cable to compress the wires onto the center tube and immobilize them. The SuperCIC is then co-wound with a high-strength armor extrusion, which is kerf-cut so that the co-winding onto a coil mandrel can be made without deformation within the armor or the CIC. SuperCIC facilitates hybrid windings, in which subwindings of NbTi are used where B < 7 T, Nb 3 Sn where B < 14 T, and Bi-2212where B > 14 T to minimize superconductor cost. Demount-able splices are used to interconnect layers. A conceptual design has been prepared for a tokamak with aspect ratio A ~ 2.0, magnetic field of 17.4Tat the inner windings and 6.7 T at the plasma,and winding current density ~140 A/mm 2 -sufficient for optimizing the fusion power density in a compact spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling CLPX IOP3 Radiometric Data by Means of the Dense Media Theory: Preliminary Results for the LSOS Test Site

The capabilities of the Dense Media Radiative Transfer model using the Quasi Crystalline Approximation with Coherent Potential (QCA-CP) to reproduce measured radiometric data were tested using the University of Tokyo Ground Based Microwave Radiometer (GBMR-7) during the third Intensive Observation Period (IOP3) of the NASA Cold-land Processes Field Experiment (CLPX). The data were collected at the Local-Scale Observation Site (LSOS), a 0.8-ha study site consisting of two open meadows separated by trees. Intensive measurements were also made of snow depth and temperature, density, and grain size profiles. A DMRT model is needed to describe radiative transfer in a medium such as snow because the assumption of independent scattering used in classical radiative transfer theory (CRT) is not valid. Validation of the DMRT approach requires a relationship between measured snow grain size and the DMRT approximation of snow grain radius as spherical particles with a mean radius of the log-normal particle-size distribution. This relationship is very important for a better understanding of snow modelling and for practical applications. DMRT simulations were compared with observations of microwave brightnesses at 18.7, 36.5 and 89-GHz (V and H polarizations) collected on February-1 9-25, 2003. Observation angles ranged from 30\deg to 70\deg. Model inputs included measured snow parameters except mean grain size. The average snow temperature, fractional volume and depth were held constant, together with the ice and soil permittivities. The minimum and maximum measured mean grain sizes were used to test the capabilities of the DMRT to reproduce the brightnesses as upper and lower limits. The sensitivity to the largest and smallest measured grain size in the three classes of minimum, medium and maximum observed grain sizes was also investigated. DMRT particle sizes yielding a best-fit to the experimental data for each date were computed. Results show that the measured brightnesses fall within the range of simulated brightnesses using the smallest and largest measured grain size values. The DMRT best-fit radii are comparable to the average radii for the medium observed grain sizes.

Tedesco, Marco↗

Temperature-Dependent Hole Transfer from Photoexcited Quantum Dots to Molecular Species: Evidence for Trap-Mediated Transfer

The effect of temperature on the rate of hole transfer from photoexcited quantum dots (QDs) is investigated by measuring the driving force dependence of the charge transfer rate for different sized QDs across a range of temperatures from 78 to 300 K. Spherical CdSe/CdS core/shell QDs were used with a series of ferrocene-derived molecular hole acceptors with an 800 meV range in electrochemical potential. Time-resolved photoluminescence measurements and photoluminescence quantum yield measurements in an integrating sphere were both performed from 78 to 300 K to obtain temperature-dependent rates for a series of driving forces as dictated by the nature of the molecular acceptor. For both QD sizes studied and all ligands, the Arrhenius plot of hole transfer exhibited an activated (linear) regime at higher temperatures and a temperature-independent regime at low temperatures. The extracted activation energies in the high-temperature regime were consistent across all ligands for a given QD size. This observation is not consistent with direct charge transfer from the QD valence band to the ferrocene acceptor. Instead, a model in which charge transfer is mediated by a shallow and reversible trap more accurately fits the experimental results. Implications for this observed trap-mediated transfer are discussed including as a strategy to more efficiently extract charge from QDs.

cadmium selenide↗

Contribution of Gaia Sausage to the Galactic Stellar Halo Revealed by K Giants and Blue Horizontal Branch Stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, Sloan Digital Sky Survey, and Gaia

We explore the contribution of the Gaia Sausage to the stellar halo of the Milky Way by making use of a Gaussian mixture model (GMM) and applying it to halo star samples of Large Sky Area Multi-Object Fiber Spectroscopic Telescope K giants, Sloan Extension for Galactic Understanding and Exploration K giants, and Sloan Digital Sky Survey blue horizontal branch stars. The GMM divides the stellar halo into two parts, of which one represents a more metal-rich and highly radially biased component associated with an ancient, head-on collision referred to as the Gaia Sausage, and the other one is a more metal-poor and isotropic halo. A symmetric bimodal Gaussian is used to describe the distribution of spherical velocity of the Gaia Sausage, and we find that the mean absolute radial velocity of the two lobes decreases with the Galactocentric radius. We find that the Gaia Sausage contributes about 41%–74% of the inner (Galactocentric radius r gc < 30 kpc) stellar halo. The fraction of stars of the Gaia Sausage starts to decline beyond r gc ~ 25–30 kpc, and the outer halo is found to be significantly less influenced by the Gaia Sausage than the inner halo. After the removal of halo substructures found by integrals of motion, the contribution of the Gaia Sausage falls slightly within r gc ~ 25 kpc but is still as high as 30%–63%. Finally, we select several possible Sausage-related substructures consisting of stars on highly eccentric orbits. The GMM/Sausage component agrees well with the selected substructure stars in their chemodynamical properties, which increases our confidence in the reliability of the GMM fits.

79 ASTRONOMY AND ASTROPHYSICS↗