Separatrix divergence of stellar streams in galactic potentials

ABSTRACT Flattened axisymmetric galactic potentials are known to host minor orbit families surrounding orbits with commensurable frequencies. The behaviour of orbits that belong to these orbit families is fundamentally different than that of typical orbits with non-commensurable frequencies. We investigate the evolution of stellar streams on orbits near the boundaries between orbit families (separatrices) in a flattened axisymmetric potential. We demonstrate that the separatrix divides these streams into two groups of stars that belong to two different orbit families, and that as a result, these streams diffuse more rapidly than streams that evolve elsewhere in the potential. We utilize Hamiltonian perturbation theory to estimate both the time-scale of this effect and the likelihood of a stream evolving close enough to a separatrix to be affected by it. We analyse two prior reports of stream-fanning in simulations with triaxial potentials, and conclude that at least one of them is caused by separatrix divergence. These results lay the foundation for a method of mapping the orbit families of galactic potentials using the morphology of stellar streams. Comparing these predictions with the currently known distribution of streams in the Milky Way presents a new way of constraining the shape of our Galaxy’s potential and distribution of dark matter.

Strong chemical tagging with APOGEE: 21 candidate star clusters that have dissolved across the Milky Way disc

ABSTRACT Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the 8D chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high-quality sample of 182 538 giant stars, we detect 21 candidate clusters with more than 15 members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in log g versus Teff space. Each group’s members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely bound groups.

Systematic errors in strong gravitational lensing reconstructions, a numerical simulation perspective

ABSTRACT We present the analysis of a sample of 24 SLACS-like galaxy–galaxy strong gravitational lens systems with a background source and deflectors from the Illustris-1 simulation. We study the degeneracy between the complex mass distribution of the lenses, substructures, the surface brightness distribution of the sources, and the time delays. Using a novel inference framework based on Approximate Bayesian Computation, we find that for all the considered lens systems, an elliptical and cored power-law mass density distribution provides a good fit to the data. However, the presence of cores in the simulated lenses affects most reconstructions in the form of a Source Position Transformation. The latter leads to a systematic underestimation of the source sizes by 50 per cent on average, and a fractional error in H0 of around $25_{-19}^{+37}$ per cent. The analysis of a control sample of 24 lens systems, for which we have perfect knowledge about the shape of the lensing potential, leads to a fractional error on H0 of $12_{-3}^{+6}$ per cent. We find no degeneracy between complexity in the lensing potential and the inferred amount of substructures. We recover an average total projected mass fraction in substructures of fsub < 1.7–2.0 × 10−3 at the 68 per cent confidence level in agreement with zero and the fact that all substructures had been removed from the simulation. Our work highlights the need for higher resolution simulations to quantify the lensing effect of more realistic galactic potentials better, and that additional observational constraint may be required to break existing degeneracies.

Hypervelocity stars in the Gaia era

Context. Young massive stars in the halo are assumed to be runaway stars from the Galactic disk. Possible ejection scenarios are binary supernova ejections (BSE) or dynamical ejections from star clusters (DE). Hypervelocity stars (HVSs) are extreme runaway stars that are potentially unbound from the Galaxy. Powerful acceleration mechanisms such as the tidal disruption of a binary system by a supermassive black hole (SMBH) are required to produce them. Therefore, HVSs are believed to originate in the Galactic center (GC), the only place known to host an SMBH. Aims. The second Gaia data release (DR2) offers the opportunity of studying HVSs in an unprecedented manner. We revisit some of the most interesting high-velocity stars, that is, 15 stars (11 candidate HVSs and 4 radial velocity outliers) for which proper motions with the Hubble Space Telescope. were obtained in the pre-Gaia era, to unravel their origin. Methods. By carrying out kinematic analyses based on revised spectrophotometric distances and proper motions from Gaia DR2, kinematic properties were obtained that help constrain the spatial origins of these stars. Results. Stars that were previously considered (un)bound remain (un)bound in Galactic potentials favored by Gaia DR2 astrometry. For nine stars (five candidate HVSs plus all four radial velocity outliers), the GC can be ruled out as spatial origin at least at 2σ confidence level, suggesting that a large portion of the known HVSs are disk runaway stars launched close to or beyond Galactic escape velocities. The fastest star in the sample, HVS 3, is confirmed to originate in the Large Magellanic Cloud. Conclusions. Because the ejection velocities of five of our non-GC stars are close to or above the upper limits predicted for BSE and DE, another powerful dynamical ejection mechanism (e.g., involving massive perturbers such as intermediate-mass black holes) is likely to operate in addition to the three classical scenarios mentioned above.

NGC 6705 a young α-enhanced open cluster from OCCASO data

Context. The stellar [α/Fe] abundance is sometimes used as a proxy for stellar age, following standard chemical evolution models for the Galaxy, as seen by different observational results. Aim. In this work, we aim to show that the open cluster NGC 6705/M 11 has a significant α-enhancement [α/Fe] > 0.1 dex, despite its young age (~300 Myr), challenging the current paradigm. Methods. We used high resolution (R > 65 000) high signal-to-noise (~70) spectra of eight red clump stars, acquired within the OCCASO survey. We determined very accurate chemical abundances of several α elements, using an equivalent width methodology (Si, Ca and Ti), and spectral synthesis fits (Mg and O). Results. We obtain [Si/Fe] = 0.13 ± 0.05, [Mg/Fe] = 0.14 ± 0.07, [O/Fe] = 0.17 ± 0.07, [Ca/Fe] = 0.06 ± 0.05, and [Ti/Fe] = 0.03 ± 0.03. Our results place these clusters within the group of young [α/Fe]-enhanced field stars recently found by several authors in the literature. The ages of our stars have an uncertainty of around 50 Myr, much more precise than for field stars. By integrating the cluster’s orbit in several non-axisymmetric Galactic potentials, we establish the M 11’s most likely birth radius as lying between 6.8–7.5 kpc from the Galactic centre, not far from its current position. Conclusions. With the robust open cluster age scale, our results prove that a moderate [α/Fe]-enhancement is no guarantee for a star to be old, and that not all α-enhanced stars can be explained with an evolved blue straggler scenario. Based on our orbit calculations, we further argue against a Galactic bar origin of M 11.

Phantom: A Smoothed Particle Hydrodynamics and Magnetohydrodynamics Code for Astrophysics

We present Phantom, a fast, parallel, modular, and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary, and high energy astrophysics, and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, dust–gas mixtures, H2 chemistry, physical viscosity, external forces including numerous galactic potentials, Lense–Thirring precession, Poynting–Robertson drag, and stochastic turbulent driving. Phantom is hereby made publicly available.

Hypervelocity star candidates in Gaia DR1/TGAS

Hypervelocity stars (HVSs) are characterized by a total velocity in excess of the Galactic escape speed, and with trajectories consistent with coming from the Galactic Centre. We apply a novel data mining routine, an artificial neural network, to discover HVSs in the TGAS subset of the first data release of the Gaia satellite, using only the astrometry of the stars. We find 80 stars with a predicted probability >90% of being HVSs, and we retrieved radial velocities for 47 of those. We discover 14 objects with a total velocity in the Galactic rest frame >400 km s−1, and 5 of these have a probability >50% of being unbound from the Milky Way. Tracing back orbits in different Galactic potentials, we discover 1 HVS candidate, 5 bound HVS candidates, and 5 runaway star candidates with remarkably high velocities, between 400 and 780 km s−1. We wait for future Gaia releases to confirm the goodness of our sample and to increase the number of HVS candidates.