Global structure and kinematics of stellar haloes in cosmological hydrodynamic simulations

We use the Galaxies–Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r≲ 30 kpc. In the present study, we show that this component was primarily formed in a protodisc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, which is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axial ratio of ∼0.6 at radii of ≲30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ versus accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a ‘dual stellar halo’

The chemo-dynamical groups of Galactic globular clusters

Galaxie

Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations

We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic ΛCDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with V 20 deg. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. (2015) that preserves the phase-space distribution of the model stars. These publicly available catalogues contain five-parameter astrometry, radial velocities, multiband photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that (i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A - and B -dwarf stars up to radii of ∼15 kpc, and (ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to (i) develop and test analysis methods for the Gaia DR2 data, (ii) gauge the limitations and biases of the data, and (iii) interpret the data in the light of theoretical predictions from realistic ab initio simulations of galaxy formation in the ΛCDM cosmological model

Stellar splashback: the edge of the intracluster light

Galaxie

Space Telescope and Optical Reverberation Mapping Project. VII. Understanding the Ultraviolet Anomaly in NGC 5548 with X-Ray Spectroscopy

During the Space Telescope and Optical Reverberation Mapping Project observations of NGC 5548, the continuum and emission-line variability became decorrelated during the second half of the six-month-long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with the Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all three observations: the UV emission-line flux decrease, the soft-excess increase, and the emission-line anomaly

Axi-symmetric models of auroral current systems in Jupiter's magnetosphere with predictions for the Juno mission

We develop two related models of magnetosphere-ionosphere coupling in the jovian system by combining previous models defined at ionospheric heights with magnetospheric magnetic models that allow system parameters to be extended appropriately into the magnetosphere. The key feature of the combined models is thus that they allow direct connection to be made between observations in the magnetosphere, particularly of the azimuthal field produced by the magnetosphere-ionosphere coupling currents and the plasma angular velocity, and the auroral response in the ionosphere. The two models are intended to reflect typical steady-state sub-corotation conditions in the jovian magnetosphere, and transient super-corotation produced by sudden major solar wind-induced compressions, respectively. The key simplification of the models is that of axi-symmetry of the field, flow, and currents about the magnetic axis, limiting their validity to radial distances within ~30 RJ of the planet, though the magnetic axis is appropriately tilted relative to the planetary spin axis and rotates with the planet. The first exploration of the jovian polar magnetosphere is planned to be undertaken in 2016–2017 during the NASA New Frontiers Juno mission, with observations of the polar field, plasma, and UV emissions as a major goal. Evaluation of the models along Juno planning orbits thus produces predictive results that may aid in science mission planning. It is shown in particular that the low-altitude near-periapsis polar passes will generally occur underneath the corresponding auroral acceleration regions, thus allowing brief examination of the auroral primaries over intervals of ~1–3 min for the main oval and ~10 s for narrower polar arc structures, while the "lagging" field deflections produced by the auroral current systems on these passes will be ~0.1°, associated with azimuthal fields above the ionosphere of a few hundred nT

The progenitors of the Milky Way stellar halo: big bricks favoured over little bricks

We present a census of blue horizontal branch (BHB) and blue straggler (BS) stars belonging to dwarf galaxies and globular clusters, and compare these counts to that of the Milky Way stellar halo. We find, in agreement with earlier studies, that the ratio of BS-to-BHB stars in these satellite populations is dependent on stellar mass. Dwarf galaxies show an increasing BS-to-BHB ratio with luminosity. In contrast, globular clusters display the reverse trend, with NBS/NBHB (≲ 1) decreasing with luminosity. The faintest (L < 105 Lȯ) dwarfs have similar numbers of BS and BHB stars (NBS/NBHB ∼ 1), whereas more-massive dwarfs tend to be dominated by BS stars (NBS/NBHB ∼ 2–40). We find that the BS-to-BHB ratio in the stellar halo is relatively high (NBS/NBHB ∼ 5–6), and thus inconsistent with the low ratios found in both ultra-faint dwarfs and globular clusters. Our results favour more-massive dwarfs as the dominant ‘building blocks’ of the stellar halo, in good agreement with current predictions from Λ cold dark matter models

The Milky Way stellar halo out to 40 kpc: squashed, broken but smooth

We introduce a new maximum likelihood method to model the density profile of Blue Horizontal Branch and Blue Straggler stars and apply it to the Sloan Digital Sky Survey Data Release 8 (DR8) photometric catalogue. There are a large number (20,000) of these tracers available over an impressive 14,000 square degrees in both Northern and Southern Galactic hemispheres, and they provide a robust measurement of the shape of the Milky Way stellar halo. After masking out stars in the vicinity of the Virgo Overdensity and the Sagittarius stream, the data are consistent with a smooth, oblate stellar halo with a density that follows a broken power-law. The best fitting model has an inner slope = 2.3 and an outer slope = 4.6, together with a break radius occurring at 27 kpc and a constant halo flattening (that is, ratio of minor axis to major axis) of q = 0.6. Although a broken power-law describes the density fall-off most adequately, it is also well fit by an Einasto profile. There is no strong evidence for variations in flattening with radius, or for triaxiality of the stellar halo.Comment: 13 pages, 12 figures. Accepted for publication by MNRA

Mass estimators in the Gaia era

Forthcoming astrometric missions such as the Gaia satellite will bring to the fore the problem of estimating the enclosed mass from a set of positions, radial velocities and proper motions of tracer stars. Here, we show how to construct the tracer mass estimator when the proper motion data are available in addition to the usual line-of-sight velocity data. Notably, the mass estimators do not require any assumption on the anisotropy, as it is naturally incorporated through the different components of proper motions. In addition, the separate treatment of the proper motions and the line-of-sight velocities is desirable because they are observationally independent and thus the propagation of the combined uncertainties is rather straightforward. The extension to projected data is also sketched, together with a possible application of measuring the masses of Galactic globular clusters.Comment: submitted to MNRAS, including 7 journal pages, 4 figures, 3 table

Accurate Mass Estimators for Navarro-Frenk-White Halos

We consider the problem of estimating the virial mass of a dark halo from the positions and velocities of a tracer population. Although a number of general tools are available, more progress can be made if we are able to specify the functional form of the halo potential (although not its normalization). Here, we consider the particular case of the cosmologically motivated Navarro-Frenk-White (NFW) halo and develop two simple estimators. We demonstrate their effectiveness against numerical simulations and use them to provide new mass estimates of Carina, Fornax, Sculptor, and Sextans dSphs.Comment: Ap J (Letters), submitte