Showing 1–50 of 93 results for author: Smith, M C

Polarisation-insensitive state preparation for trapped-ion hyperfine qubits

Authors: A. D. Leu, M. C. Smith, M. F. Gely, D. M. Lucas

Abstract: Quantum state preparation for trapped-ion qubits often relies on high-quality circularly-polarised light, which may be difficult to achieve with chip-based integrated optics technology. We propose and implement a hybrid optical/microwave scheme for intermediate-field hyperfine qubits which instead relies on frequency selectivity. Experimentally, we achieve $99.94\%$ fidelity for linearly-polarised… ▽ More Quantum state preparation for trapped-ion qubits often relies on high-quality circularly-polarised light, which may be difficult to achieve with chip-based integrated optics technology. We propose and implement a hybrid optical/microwave scheme for intermediate-field hyperfine qubits which instead relies on frequency selectivity. Experimentally, we achieve $99.94\%$ fidelity for linearly-polarised ($σ^+$/$σ^-$) light, using $^{43}$Ca$^+$ at 28.8 mT. We find that the fidelity remains above $99.8\%$ for a mixture of all polarisations ($σ^+$/$σ^-$/$π$). We calculate that the method is capable of $99.99\%$ fidelity in $^{43}$Ca$^+$, and even higher fidelities in heavier ions such as $^\text{137}$Ba$^\text{+}$. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Focusing of quantum gate interactions using dynamical decoupling

Authors: M. C. Smith, A. D. Leu, M. F. Gely, D. M. Lucas

Abstract: In 1995, Cirac and Zoller proposed the first concrete implementation of a small-scale quantum computer, using laser beams focused to micron spot sizes to address individual trapped ions in a linear crystal. Here we propose a method to focus entangling gate interactions, but driven by microwave fields, to micron-sized zones, corresponding to $10^{-5}$ microwave wavelengths. We demonstrate the abili… ▽ More In 1995, Cirac and Zoller proposed the first concrete implementation of a small-scale quantum computer, using laser beams focused to micron spot sizes to address individual trapped ions in a linear crystal. Here we propose a method to focus entangling gate interactions, but driven by microwave fields, to micron-sized zones, corresponding to $10^{-5}$ microwave wavelengths. We demonstrate the ability to suppress the spin-dependent force using a single ion, and find the required interaction introduces $3.7(4)\times 10^{-4}$ error per emulated gate in a single-qubit benchmarking sequence. We model the scheme for a 17-qubit ion crystal, and find that any pair of ions should be addressable with an average crosstalk error of $\sim 10^{-5}$. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Fast, high-fidelity addressed single-qubit gates using efficient composite pulse sequences

Authors: A. D. Leu, M. F. Gely, M. A. Weber, M. C. Smith, D. P. Nadlinger, D. M. Lucas

Abstract: We use electronic microwave control methods to implement addressed single-qubit gates with high speed and fidelity, for $^{43}\text{Ca}^{+}$ hyperfine "atomic clock" qubits in a cryogenic (100K) surface trap. For a single qubit, we benchmark an error of $1.5$ $\times$ $10^{-6}$ per Clifford gate (implemented using $600~\text{ns}$ $π/2$-pulses). For two qubits in the same trap zone (ion separation… ▽ More We use electronic microwave control methods to implement addressed single-qubit gates with high speed and fidelity, for $^{43}\text{Ca}^{+}$ hyperfine "atomic clock" qubits in a cryogenic (100K) surface trap. For a single qubit, we benchmark an error of $1.5$ $\times$ $10^{-6}$ per Clifford gate (implemented using $600~\text{ns}$ $π/2$-pulses). For two qubits in the same trap zone (ion separation $5~μ\text{m}$), we use a spatial microwave field gradient, combined with an efficient 4-pulse scheme, to implement independent addressed gates. Parallel randomized benchmarking on both qubits yields an average error $3.4$ $\times$ $10^{-5}$ per addressed $π/2$-gate. The scheme scales theoretically to larger numbers of qubits in a single register. △ Less

Submitted 20 September, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Journal ref: Phys. Rev. Lett. 131, 120601 (2023)

Arkenstone I: A Novel Method for Robustly Capturing High Specific Energy Outflows In Cosmological Simulations

Authors: Matthew C. Smith, Drummond B. Fielding, Greg L. Bryan, Chang-Goo Kim, Eve C. Ostriker, Rachel S. Somerville, Jonathan Stern, Kung-Yi Su, Rainer Weinberger, Chia-Yu Hu, John C. Forbes, Lars Hernquist, Blakesley Burkhart, Yuan Li

Abstract: Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arep… ▽ More Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model. △ Less

Submitted 7 March, 2024; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

Journal ref: MNRAS, 527, 1216 (2024)

Bounded Disturbance Amplification for Mass Chains with Passive Interconnection

Authors: Kaoru Yamamoto, Malcolm C. Smith

Abstract: This paper introduces the problem of passive control of a chain of N identical masses in which there is an identical passive connection between neighbouring masses and a similar connection to a movable point. The problem arises in the design of multi-storey buildings which are subjected to earthquake disturbances, but applies in other situations, for example vehicle platoons. The paper studies the… ▽ More This paper introduces the problem of passive control of a chain of N identical masses in which there is an identical passive connection between neighbouring masses and a similar connection to a movable point. The problem arises in the design of multi-storey buildings which are subjected to earthquake disturbances, but applies in other situations, for example vehicle platoons. The paper studies the scalar transfer functions from the disturbance to a given intermass displacement. It is shown that these transfer functions can be conveniently represented in the form of complex iterative maps and that these maps provide a method to establish boundedness in N of the H-infinity norm of these transfer functions for certain choices of interconnection impedance. △ Less

Submitted 24 December, 2022; originally announced December 2022.

MSC Class: 93A14; 93A15

Journal ref: IEEE Transactions on Automatic Control (2015) vol. 61 no. 6, pp. 1565 - 1574

Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods

Authors: Chia-Yu Hu, Matthew C. Smith, Romain Teyssier, Greg L. Bryan, Robbert Verbeke, Andrew Emerick, Rachel S. Somerville, Blakesley Burkhart, Yuan Li, John C. Forbes, Tjitske Starkenburg

Abstract: We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid mo… ▽ More We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density-temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger supernova-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: gas in Lagrangian codes collapses to much higher densities than in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed. △ Less

Submitted 3 May, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: accepted version by ApJ (including a new simulation in Appendix B suggested by the referee)

Two can play at that game: constraining the role of supernova and AGN feedback in dwarf galaxies with cosmological zoom-in simulations

Authors: Sophie Koudmani, Debora Sijacki, Matthew C. Smith

Abstract: There is growing observational evidence for dwarf galaxies hosting active galactic nuclei (AGN), including hints of AGN-driven outflows in dwarfs. However, in the common theoretical model of galaxy formation, efficient supernova (SN) feedback is the tool of choice for regulating star formation in the low-mass regime. In this paper, we present a suite of high-resolution cosmological dwarf zoom-in s… ▽ More There is growing observational evidence for dwarf galaxies hosting active galactic nuclei (AGN), including hints of AGN-driven outflows in dwarfs. However, in the common theoretical model of galaxy formation, efficient supernova (SN) feedback is the tool of choice for regulating star formation in the low-mass regime. In this paper, we present a suite of high-resolution cosmological dwarf zoom-in simulations relaxing the assumption of strong SN feedback, with the goal to determine whether more moderate SN feedback in combination with an efficient AGN could be a suitable alternative. Importantly, we find that there are sufficient amounts of gas to power brief Eddington-limited accretion episodes in dwarfs. This leads to a variety of outcomes depending on the AGN accretion model: from no additional suppression to moderate regulation of star formation to catastrophic quenching. Efficient AGN can drive powerful outflows, depleting the gas reservoir of their hosts via ejective feedback and then maintaining a quiescent state through heating the circumgalactic medium. Moderate AGN outflows can be as efficient as the strong SN feedback commonly employed, leading to star formation regulation and HI gas masses in agreement with observations of field dwarfs. All efficient AGN set-ups are associated with overmassive black holes (BHs) compared to the (heavily extrapolated) observed BH mass - stellar mass scaling relations, with future direct observational constraints in this mass regime being crucially needed. Efficient AGN activity is mostly restricted to high redshifts, with hot, accelerated outflows and high X-ray luminosities being the clearest tell-tale signs for future observational campaigns. △ Less

Submitted 1 September, 2022; v1 submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS; 30 pages, 14 figures (including 3 appendices)

Searching Extra-tidal Features around the Globular Cluster Whiting 1

Authors: Jundan Nie, Hao Tian, Jing Li, Chao Liu, Martin C. Smith, Baitian Tang, Julio A. Carballo-Bello, Jun Ma, Haijun Tian, Jiaxin Wang, Zhenyu Wu, Xiyan Peng, Jiali Wang, Tianmeng Zhang, Xu Zhou, Zhimin Zhou, Hu Zou

Abstract: Whiting 1 is a faint and young globular cluster in the halo of the Milky Way, and was suggested to have originated in the Sagittarius spherical dwarf galaxy (Sgr dSph). In this paper, we use the deep DESI Legacy Imaging Surveys to explore tentative spatial connection between Whiting 1 and the Sgr dSph. We redetermine the fundamental parameters of Whiting 1 and use the best-fitting isochrone (age… ▽ More Whiting 1 is a faint and young globular cluster in the halo of the Milky Way, and was suggested to have originated in the Sagittarius spherical dwarf galaxy (Sgr dSph). In this paper, we use the deep DESI Legacy Imaging Surveys to explore tentative spatial connection between Whiting 1 and the Sgr dSph. We redetermine the fundamental parameters of Whiting 1 and use the best-fitting isochrone (age $τ$=6.5 Gyr, metalicity Z=0.005 and $\rm d_{\odot}$=26.9 kpc) to construct a theoretical matched filter for the extra-tidal features searching. Without any smooth technique to the matched filter density map, we detect a round-shape feature with possible leading and trailing tails on either side of the cluster. This raw image is not totally new compared to old discoveries, but confirms that no more large-scale features can be detected under a depth of r<=22.5 mag. In our results, the whole feature stretches 0.1-0.2 degree along the orbit of Whiting 1, which gives a much larger area than the cluster core. The tails on both sides of the cluster align along the orbital direction of the Sgr dSph as well as the cluster itself, which implies that these debris are probably stripped remnants of Whiting 1 by the Milky Way. △ Less

Submitted 25 April, 2022; v1 submitted 23 April, 2022; originally announced April 2022.

Comments: 10 pages, 9 figures, accepted for the publication of ApJ

The evolution of scientific literature as metastable knowledge states

Authors: Sai Dileep Koneru, David Rench McCauley, Michael C. Smith, David Guarrera, Jenn Robinson, Sarah Rajtmajer

Abstract: The problem of identifying common concepts in the sciences and deciding when new ideas have emerged is an open one. Metascience researchers have sought to formalize principles underlying stages in the life-cycle of scientific research, determine how knowledge is transferred between scientists and stakeholders, and understand how new ideas are generated and take hold. Here, we model the state of sc… ▽ More The problem of identifying common concepts in the sciences and deciding when new ideas have emerged is an open one. Metascience researchers have sought to formalize principles underlying stages in the life-cycle of scientific research, determine how knowledge is transferred between scientists and stakeholders, and understand how new ideas are generated and take hold. Here, we model the state of scientific knowledge immediately preceding new directions of research as a metastable state and the creation of new concepts as combinatorial innovation. We find that, through the combined use of natural language clustering and citation graph analysis, we can predict the evolution of ideas over time and thus connect a single scientific article to past and future concepts in a way that goes beyond traditional citation and reference connections. △ Less

Submitted 11 September, 2022; v1 submitted 25 February, 2022; originally announced February 2022.

Towards Aggregated Asynchronous Checkpointing

Authors: Mikaila J. Gossman, Bogdan Nicolae, Jon C. Calhoun, Franck Cappello, Melissa C. Smith

Abstract: High-Performance Computing (HPC) applications need to checkpoint massive amounts of data at scale. Multi-level asynchronous checkpoint runtimes like VELOC (Very Low Overhead Checkpoint Strategy) are gaining popularity among application scientists for their ability to leverage fast node-local storage and flush independently to stable, external storage (e.g., parallel file systems) in the background… ▽ More High-Performance Computing (HPC) applications need to checkpoint massive amounts of data at scale. Multi-level asynchronous checkpoint runtimes like VELOC (Very Low Overhead Checkpoint Strategy) are gaining popularity among application scientists for their ability to leverage fast node-local storage and flush independently to stable, external storage (e.g., parallel file systems) in the background. Currently, VELOC adopts a one-file-per-process flush strategy, which results in a large number of files being written to external storage, thereby overwhelming metadata servers and making it difficult to transfer and access checkpoints as a whole. This paper discusses the viability and challenges of designing aggregation techniques for asynchronous multi-level checkpointing. To this end we implement and study two aggregation strategies, their limitations, and propose a new aggregation strategy specifically for asynchronous multi-level checkpointing. △ Less

Submitted 4 December, 2021; originally announced December 2021.

Comments: Accepted submission to the SuperCheck Workshop at the SuperComputing Conference held in St. Louis, MO. November 14-19, 2021(SC'21)

Speech Enhancement Using Multi-Stage Self-Attentive Temporal Convolutional Networks

Authors: Ju Lin, Adriaan J. van Wijngaarden, Kuang-Ching Wang, Melissa C. Smith

Abstract: Multi-stage learning is an effective technique to invoke multiple deep-learning modules sequentially. This paper applies multi-stage learning to speech enhancement by using a multi-stage structure, where each stage comprises a self-attention (SA) block followed by stacks of temporal convolutional network (TCN) blocks with doubling dilation factors. Each stage generates a prediction that is refined… ▽ More Multi-stage learning is an effective technique to invoke multiple deep-learning modules sequentially. This paper applies multi-stage learning to speech enhancement by using a multi-stage structure, where each stage comprises a self-attention (SA) block followed by stacks of temporal convolutional network (TCN) blocks with doubling dilation factors. Each stage generates a prediction that is refined in a subsequent stage. A fusion block is inserted at the input of later stages to re-inject original information. The resulting multi-stage speech enhancement system, in short, multi-stage SA-TCN, is compared with state-of-the-art deep-learning speech enhancement methods using the LibriSpeech and VCTK data sets. The multi-stage SA-TCN system's hyper-parameters are fine-tuned, and the impact of the SA block, the fusion block and the number of stages are determined. The use of a multi-stage SA-TCN system as a front-end for automatic speech recognition systems is investigated as well. It is shown that the multi-stage SA-TCN systems perform well relative to other state-of-the-art systems in terms of speech enhancement and speech recognition scores. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: Preprint

A LAMOST BHB Catalog and Kinematics Therein I: Catalog and Halo Properties

Authors: John J. Vickers, Zhao-Yu Li, Martin C. Smith, Juntai Shen

Abstract: We collect a sample of stars observed both in LAMOST and Gaia which have colors implying a temperature hotter than 7000 K. We train a machine learning algorithm on LAMOST spectroscopic data which has been tagged with stellar classifications and metallicities, and use this machine to construct a catalog of Blue Horizontal Branch stars (BHBs) with metallicity information. Another machine is trained… ▽ More We collect a sample of stars observed both in LAMOST and Gaia which have colors implying a temperature hotter than 7000 K. We train a machine learning algorithm on LAMOST spectroscopic data which has been tagged with stellar classifications and metallicities, and use this machine to construct a catalog of Blue Horizontal Branch stars (BHBs) with metallicity information. Another machine is trained using Gaia parallaxes to predict absolute magnitudes for these stars. The final catalog of 13,693 BHBs is thought to be about 86\% pure, with $σ_{[Fe/H]}\sim$0.35 dex and $σ_{G}\sim$0.31 mag. These values are confirmed via comparison to globular clusters, although a covariance error seems to affect our magnitude and abundance estimates. We analyze a subset of this catalog in the Galactic Halo. We find that BHB populations in the outer halo appear redder, which could imply a younger population, and that the metallicity gradient is relatively flat around [Fe/H] = -1.9 dex over our sample footprint. We find that our metal rich BHB stars are on more radial velocity dispersion dominated orbits ($β\sim 0.70$) at all radii than our metal poor BHB stars ($β\sim 0.62$). △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 13 Pages, 7 Figures, accepted to apj

The sensitivity of stellar feedback to IMF averaging versus IMF sampling in galaxy formation simulations

Authors: Matthew C. Smith

Abstract: Galaxy formation simulations frequently use Initial Mass Function (IMF) averaged feedback prescriptions, where star particles are assumed to represent single stellar populations that fully sample the IMF. This approximation breaks down at high mass resolution, where stochastic variations in stellar populations become important. We discuss various schemes to populate star particles with stellar mas… ▽ More Galaxy formation simulations frequently use Initial Mass Function (IMF) averaged feedback prescriptions, where star particles are assumed to represent single stellar populations that fully sample the IMF. This approximation breaks down at high mass resolution, where stochastic variations in stellar populations become important. We discuss various schemes to populate star particles with stellar masses explicitly sampled from the IMF. We use Monte Carlo numerical experiments to examine the ability of the schemes to reproduce an input IMF in an unbiased manner while conserving mass. We present our preferred scheme which can easily be added to pre-existing star formation prescriptions. We then carry out a series of high resolution isolated simulations of dwarf galaxies with supernovae, photoionization and photoelectric heating to compare the differences between using IMF averaged feedback and explicitly sampling the IMF. We find that if supernovae are the only form of feedback, triggering individual supernovae from IMF averaged rates gives identical results to IMF sampling. However, we find that photoionization is more effective at regulating star formation when IMF averaged rates are used, creating more, smaller H II regions than the rare, bright sources produced by IMF sampling. We note that the increased efficiency of the IMF averaged feedback versus IMF sampling is not necessarily a general trend and may be reversed depending on feedback channel, resolution and other details. However, IMF sampling is always the more physically motivated approach. We conservatively suggest that it should be used for star particles less massive than $\sim500\,\mathrm{M_\odot}$. △ Less

Submitted 20 April, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 23 pages, 14 figures, published MNRAS. Minor alterations to reflect published version

Journal ref: MNRAS, 502, 5417 (2021)

A Framework for Multiphase Galactic Wind Launching using TIGRESS

Authors: Chang-Goo Kim, Eve C. Ostriker, Drummond B. Fielding, Matthew C. Smith, Greg L. Bryan, Rachel S. Somerville, John C. Forbes, Shy Genel, Lars Hernquist

Abstract: Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations using the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing ga… ▽ More Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations using the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing gas is characterized by two distinct thermal phases, cool (T<10^4 K) and hot (T>10^6 K), with most mass carried by the cool phase and most energy and newly-injected metals carried by the hot phase. Both components have a broad distribution of outflow velocity, and especially for cool gas this implies a varying fraction of escaping material depending on the halo potential. Informed by the TIGRESS results, we develop straightforward analytic formulae for the joint probability density functions (PDFs) of mass, momentum, energy, and metal loading as distributions in outflow velocity and sound speed. The model PDFs have only two parameters, SFR surface density Σ_SFR and the metallicity of the ISM, and fully capture the behavior of the original TIGRESS simulation PDFs over Σ_SFR~(10^{-4},1)M_sun/kpc^2/yr. Employing PDFs from resolved simulations will enable galaxy formation subgrid model implementations with wind velocity and temperature (as well as total loading factors) that are based on theoretical predictions rather than empirical tuning. This is a critical step to incorporate advances from TIGRESS and other high-resolution simulations in future cosmological hydrodynamics and semi-analytic galaxy formation models. We release a python package to prototype our model and to ease its implementation. △ Less

Submitted 18 October, 2020; originally announced October 2020.

Comments: ApJL accepted. For associated python package, see https://twind.readthedocs.io/en/latest/

LAMOST Observations in 15 \textit{K}2 Campaigns: I. Low resolution spectra from LAMOST DR6

Authors: Jiangtao Wang, Jian-Ning Fu, Weikai Zong, M. C. Smith, Peter De Cat, Jianrong Shi, Ali Luo, Haotong Zhang, A. Frasca, C. J. Corbally, J. Molenda-Żakowicz, G. Catanzaro, R. O. Gray, Jiaxin Wang, Yang Pan

Abstract: The LAMOST-\textit{K}2 (L\textit{K}2) project, initiated in 2015, aims to collect low-resolution spectra of targets in the \textit{K}2 campaigns, similar to LAMOST-\textit{Kepler} project. By the end of 2018, a total of 126 L\textit{K}2 plates had been observed by LAMOST. After cross-matching the catalog of the LAMOST data release 6 (DR6) with that of the \textit{K}2 approved targets, we found 160… ▽ More The LAMOST-\textit{K}2 (L\textit{K}2) project, initiated in 2015, aims to collect low-resolution spectra of targets in the \textit{K}2 campaigns, similar to LAMOST-\textit{Kepler} project. By the end of 2018, a total of 126 L\textit{K}2 plates had been observed by LAMOST. After cross-matching the catalog of the LAMOST data release 6 (DR6) with that of the \textit{K}2 approved targets, we found 160,619 usable spectra of 84,012 objects, most of which had been observed more than once. The effective temperature, surface gravity, metallicity, and radial velocity from 129,974 spectra for 70,895 objects are derived through the LAMOST Stellar Parameter Pipeline (LASP). The internal uncertainties were estimated to be 81 K, 0.15 dex, 0.09 dex and 5 kms$^{-1}$, respectively, when derived from a spectrum with a signal-to-noise ratio in the $g$ band (SNR$_g$) of 10. These estimates are based on results for targets with multiple visits. The external accuracies were assessed by comparing the parameters of targets in common with the APOGEE and GAIA surveys, for which we generally found linear relationships. A final calibration is provided, combining external and internal uncertainties for giants and dwarfs, separately. We foresee that these spectroscopic data will be used widely in different research fields, especially in combination with \textit{K}2 photometry. △ Less

Submitted 14 October, 2020; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: 31 pages, 9 figures, 6 tables, accepted by ApJS

Efficient early stellar feedback can suppress galactic outflows by reducing supernova clustering

Authors: Matthew C. Smith, Greg L. Bryan, Rachel S. Somerville, Chia-Yu Hu, Romain Teyssier, Blakesley Burkhart, Lars Hernquist

Abstract: We present a novel set of stellar feedback models, implemented in the moving-mesh code Arepo, designed for galaxy formation simulations with near-parsec (or better) resolution. These include explicit sampling of stars from the IMF, allowing feedback to be linked to individual massive stars, an improved method for the modelling of H II regions, photoelectric heating from a spatially varying FUV fie… ▽ More We present a novel set of stellar feedback models, implemented in the moving-mesh code Arepo, designed for galaxy formation simulations with near-parsec (or better) resolution. These include explicit sampling of stars from the IMF, allowing feedback to be linked to individual massive stars, an improved method for the modelling of H II regions, photoelectric heating from a spatially varying FUV field and supernova feedback. We perform a suite of 32 simulations of isolated $M_\mathrm{vir} = 10^{10}\,\mathrm{M_\odot}$ galaxies with a baryonic mass resolution of $20\,\mathrm{M_\odot}$ in order to study the non-linear coupling of the different feedback channels. We find that photoionization and supernova feedback are both independently capable of regulating star formation to the same level, while photoelectric heating is inefficient. Photoionization produces a considerably smoother star formation history than supernovae. When all feedback channels are combined, the additional suppression of star formation rates is minor. However, outflow rates are substantially reduced relative to the supernova only simulations. We show that this is directly caused by a suppression of supernova clustering by the photoionization feedback, disrupting star forming clouds prior to the first supernovae. We demonstrate that our results are robust to variations of our star formation prescription, feedback models and the gas fraction of the disk. Our results also imply that the burstiness of star formation and the mass loading of outflows may be overestimated if the adopted star particle mass is considerably larger than the mass of individual stars because this imposes a minimum cluster size. △ Less

Submitted 20 August, 2021; v1 submitted 23 September, 2020; originally announced September 2020.

Comments: Published in MNRAS, 37 pages (including appendices), 23 figures. Updated for consistency with published version

Journal ref: 2021, MNRAS, Volume 506, Issue 3, pp.3882-3915

First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

Authors: Bhawna Motwani, Shy Genel, Greg L. Bryan, Chang-Goo Kim, Viraj Pandya, Rachel S. Somerville, Matthew C. Smith, Eve C. Ostriker, Dylan Nelson, Annalisa Pillepich, John C. Forbes, Francesco Belfiore, Rüdiger Pakmor, Lars Hernquist

Abstract: Physical and chemical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of 8 resolved ISM properties in star-… ▽ More Physical and chemical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space. △ Less

Submitted 22 October, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: Updated intro, minor text and abstract modifications. More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

The Splash without a merger

Authors: João A. S. Amarante, Leandro Beraldo e Silva, Victor P. Debattista, Martin C. Smith

Abstract: The Milky Way's progenitor experienced several merger events which left their imprints on the stellar halo, including the Gaia-Sausage/Enceladus. Recently, it has been proposed that this event perturbed the proto-disk and gave rise to a metal rich ([Fe/H] $>-1$) low angular momentum ($v_φ < 100$ km/s) stellar population. These stars have dynamical and chemical properties different from the accrete… ▽ More The Milky Way's progenitor experienced several merger events which left their imprints on the stellar halo, including the Gaia-Sausage/Enceladus. Recently, it has been proposed that this event perturbed the proto-disk and gave rise to a metal rich ([Fe/H] $>-1$) low angular momentum ($v_φ < 100$ km/s) stellar population. These stars have dynamical and chemical properties different from the accreted stellar halo, but are continuous with the canonical thick disk. In this letter, we use a hydrodynamical simulation of an isolated galaxy which develops clumps that produce a bimodal thin$+$thick disk chemistry to explore whether it forms such a population. We demonstrate clump scattering forms a metal-rich, low angular momentum population, without the need for a major merger. We show that, in the simulation, these stars have chemistry, kinematics and density distribution in good agreement with those in the Milky Way. △ Less

Submitted 21 February, 2020; v1 submitted 29 December, 2019; originally announced December 2019.

Comments: Accepted for publication in ApJL

The tale of the tail -- disentangling the high transverse velocity stars in Gaia DR2

Authors: João A. S. Amarante, Martin C. Smith, Corrado Boeche

Abstract: Although the stellar halo accounts for just ~1% of the total stellar mass of the Milky Way, the kinematics of halo stars encode valuable information about the origins and evolution of our Galaxy. It has been shown that the high transverse velocity stars in Gaia DR2 reveal a double sequence in the Hertzsprung-Russell (HR) diagram, indicating a bifurcation in the local stellar halo within 1 kpc. We… ▽ More Although the stellar halo accounts for just ~1% of the total stellar mass of the Milky Way, the kinematics of halo stars encode valuable information about the origins and evolution of our Galaxy. It has been shown that the high transverse velocity stars in Gaia DR2 reveal a double sequence in the Hertzsprung-Russell (HR) diagram, indicating a bifurcation in the local stellar halo within 1 kpc. We fit these stars by updating the popular Besançon/Galaxia model, incorporating the latest observational results for the stellar halo and an improved kinematic description for the thick-disk from Schönrich & Binney (2012). We are able to obtain a good match to the Gaia data and provide new constraints on the properties of the Galactic disc and stellar halo. In particular, we show that the kinematically defined thick disc contribution to this high velocity tail is ~13%. We look in greater detail using chemistry from LAMOST DR5, identifying a population of retrograde stars with thick-disc chemistry. Our thick disc kinematic model cannot account for this population and so we conclude there is likely to be a contribution from heated or accreted stars in the Solar Neighbourhood. We also investigate proposed dynamical substructures in this sample, concluding that they are probably due to resonant orbits rather than accreted populations. Finally we provide new insights on the nature of the two sequences and their relation with past accretion events and the primordial Galactic disc. △ Less

Submitted 11 January, 2020; v1 submitted 29 December, 2019; originally announced December 2019.

Comments: Accepted for publication in MNRAS

Report number: MN-19-3099-MJ.R2

Young stars as tracers of a barred-spiral Milky Way

Authors: Alex R. Pettitt, Sarah E. Ragan, Martin C. Smith

Abstract: Identifying the structure of our Galaxy has always been fraught with difficulties, and while modern surveys continue to make progress building a map of the Milky Way, there is still much to understand. The arm and bar features are important drivers in shaping the interstellar medium, but their exact nature and influence still require attention. We present results of smoothed particle hydrodynamic… ▽ More Identifying the structure of our Galaxy has always been fraught with difficulties, and while modern surveys continue to make progress building a map of the Milky Way, there is still much to understand. The arm and bar features are important drivers in shaping the interstellar medium, but their exact nature and influence still require attention. We present results of smoothed particle hydrodynamic simulations of gas in the Milky Way including star formation, stellar feedback, and ISM cooling, when exposed to different arm and bar features, with the aim of better understanding how well newly formed stars trace out the underlying structure of the Galaxy. The bar is given a faster pattern speed than the arms, resulting in a complex, time-dependent morphology and star formation. Inter-arm branches and spurs are easily influenced by the bar, especially in the two-armed spiral models where there is a wide region of resonance overlap in the disc. As the bar over-takes the spiral arms it induces small boosts in star formation and enhances spiral features, which occur at regularly spaced beat-like intervals. The locations of star formation events are similar to those seen in observational data, and do not show a perfect 1:1 correspondence with the underlying spiral potential, though arm tangencies are generally well traced by young stars. Stellar velocity fields from the newly formed stars are compared to data from Gaia DR2, showing that the spiral and bar features can reproduce many of the non-axisymmetric features seen in the data. A simple analytical model is used to show many of these feature are a natural response of gas to rigidly rotating spiral and bar potentials. △ Less

Submitted 10 November, 2019; originally announced November 2019.

Comments: 21 pages, 26 figures, accepted for publication in MNRAS. Movies of simulations included in this work can be found at: https://zenodo.org/record/3239994

Dynamical Relics of the Ancient Galactic Halo

Authors: Zhen Yuan, G. C. Myeong, Timothy C. Beers, N. Wyn Evans, Young Sun Lee, Projjwal Banerjee, Dmitrii Gudin, Kohei Hattori, Haining Li, Tadafumi Matsuno, Vinicius M. Placco, M. C. Smith, Devin D. Whitten, Gang Zhao

Abstract: We search for dynamical substructures in the LAMOST DR3 very metal-poor (VMP) star catalog. After cross-matching with Gaia DR2, there are 3300 VMP stars with available high-quality astrometric information that have halo-like kinematics. We apply a method based on self-organizing maps to find groups clustered in the 4D space of orbital energy and angular momentum. We identify 57 dynamically tagged… ▽ More We search for dynamical substructures in the LAMOST DR3 very metal-poor (VMP) star catalog. After cross-matching with Gaia DR2, there are 3300 VMP stars with available high-quality astrometric information that have halo-like kinematics. We apply a method based on self-organizing maps to find groups clustered in the 4D space of orbital energy and angular momentum. We identify 57 dynamically tagged groups, which we label DTG-1 to DTG-57. Most of them belong to existing substructures in the nearby halo, such as the $Gaia$ Sausage or Sequoia. The stream identified by Helmi et al. is recovered, but the two disjoint portions of the substructure have distinct dynamical properties. The very retrograde substructure Rg5 found previously by Myeong et al. is also retrieved. We report 6 new DTGs with highly retrograde orbits, 2 with very prograde orbits, and 12 with polar orbits. By mapping other datasets (APOGEE halo stars, and catalogs of r-process-enhanced and CEMP stars) onto the trained neuron map, we can associate stars with detailed chemical abundances to the DTGs, and look for associations with chemically peculiar stars. The highly eccentric $Gaia$ Sausage groups contain representatives both of debris from the satellite itself (which is $α$-poor) and the Splashed Disk, sent up into eccentric halo orbits from the encounter (and is $α$-rich). The new prograde substructures also appear to be associated with the Splashed Disk. The DTGs belonging to the $Gaia$ Sausage host two relatively metal-rich $r$-II stars and six CEMP stars in different sub-classes, consistent with the idea that the $Gaia$ Sausage progenitor is a massive dwarf galaxy. Rg5 is dynamically associated with two highly $r$-process-enhanced stars with [Fe/H] $\sim -$3. This finding indicates that its progenitor might be an ultra-faint dwarf galaxy that has experienced $r$-process enrichment from neutron star mergers. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: ApJ, submitted

The Biggest Splash

Authors: Vasily Belokurov, Jason L. Sanders, Azadeh Fattahi, Martin C. Smith, Alis J. Deason, N. Wyn Evans, Robert J. J. Grand

Abstract: Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of element… ▽ More Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics and stellar ages in the plane of azimuthal velocity v_phi and metallicity [Fe/H], we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the alpha-rich (or "thick") disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the Milky Way's proto-disc prior to the massive ancient accretion event which drastically altered their orbits. We can not, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the Milky Way's formation. △ Less

Submitted 10 September, 2019; originally announced September 2019.

Comments: Submitted to MNRAS. Comments are welcome!

Chemo-dynamical properties of the Anticenter Stream: a surviving disc fossil from a past satellite interaction

Authors: Chervin F. P. Laporte, Vasily Belokurov, Sergey E. Koposov, Martin C. Smith, Vanessa Hill

Abstract: Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly… ▽ More Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars ($\sim 10\,\rm{Gyr}$), whereas Monoceros has a pronounced tail of younger stars ($6-10\, \rm{Gyr}$) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 6 pages, 6 Figure, submitted to MNRAS, comments welcome

Cellular State Transformations using Generative Adversarial Networks

Authors: Colin Targonski, Benjamin T. Shealy, Melissa C. Smith, F. Alex Feltus

Abstract: We introduce a novel method to unite deep learning with biology by which generative adversarial networks (GANs) generate transcriptome perturbations and reveal condition-defining gene expression patterns. We find that a generator conditioned to perturb any input gene expression profile simulates a realistic transition between source and target RNA expression states. The perturbed samples follow a… ▽ More We introduce a novel method to unite deep learning with biology by which generative adversarial networks (GANs) generate transcriptome perturbations and reveal condition-defining gene expression patterns. We find that a generator conditioned to perturb any input gene expression profile simulates a realistic transition between source and target RNA expression states. The perturbed samples follow a similar distribution to original samples from the dataset, also suggesting these are biologically meaningful perturbations. Finally, we show that it is possible to identify the genes most positively and negatively perturbed by the generator and that the enriched biological function of the perturbed genes are realistic. We call the framework the Transcriptome State Perturbation Generator (TSPG), which is open source software available at https://github.com/ctargon/TSPG. △ Less

Submitted 28 June, 2019; originally announced July 2019.

Comments: 11 pages, 5 figures

Principles of lossless adjustable one-ports

Authors: Tryphon T. Georgiou, Faryar Jabbari, Malcolm C. Smith

Abstract: This paper explores the possibility to construct two-terminal mechanical devices (one-ports) which are lossless and adjustable. To be lossless, the device must be passive (i.e. not requiring a power supply) and non-dissipative. To be adjustable, a parameter of the device should be freely variable in real time as a control input. For the simplest lossless one ports, the spring and inerter, the ques… ▽ More This paper explores the possibility to construct two-terminal mechanical devices (one-ports) which are lossless and adjustable. To be lossless, the device must be passive (i.e. not requiring a power supply) and non-dissipative. To be adjustable, a parameter of the device should be freely variable in real time as a control input. For the simplest lossless one ports, the spring and inerter, the question is whether the stiffness and inertance may be varied freely in a lossless manner. We will show that the typical laws which have been proposed for adjustable springs and inerters are necessarily active and that it is not straightforward to modify them to achieve losslessness, or indeed passivity. By means of a physical construction using a lever with moveable fulcrum we will derive device laws for adjustable springs and inerters which satisfy a formal definition of losslessness. We further provide a construction method which does not require a power supply for physically realisable translational and rotary springs and inerters. The analogous questions for lossless adjustable electrical devices are examined. △ Less

Submitted 15 May, 2019; originally announced May 2019.

Comments: 10 pages, 8 figures

MSC Class: 93A99; 70Q99; 70E99

The Pisces Plume and the Magellanic Wake

Authors: Vasily Belokurov, Alis J. Deason, Denis Erkal, Sergey E. Koposov, Julio A. Carballo-Bello, Martin C Smith, Prashin Jethwa, Camila Navarrete

Abstract: Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream.… ▽ More Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream. Using follow-up VLT FORS2 spectroscopy, we have measured the velocity distribution of the Pisces candidate member stars and have shown it to be as broad as that of the Galactic halo but offset to negative velocities. Using a suite of numerical simulations, we demonstrate that the structure has many properties in common with the predicted behaviour of the Magellanic wake, i.e. the Galactic halo overdensity induced by the in-fall of the Magellanic Clouds. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: Submitted to MNRAS Letters. Comments welcome!

Revealing the Complicated Story of the Cetus Stream with StarGO

Authors: Zhen Yuan, M. C. Smith, Xiang-Xiang Xue, Jing Li, Chao Liu, Yue Wang, Lu Li, Jiang Chang

Abstract: We use a novel cluster identification tool StarGO to explore the metal poor ([Fe/H] $<$ -1.5) outer stellar halo (d $>$ 15 kpc) of the Milky Way using data from Gaia, LAMOST and SDSS. Our method is built using an unsupervised learning algorithm, a self-organizing map, which trains a 2-D neural network to learn the topological structures of a data set from an n-D input space. Using a 4-D space of a… ▽ More We use a novel cluster identification tool StarGO to explore the metal poor ([Fe/H] $<$ -1.5) outer stellar halo (d $>$ 15 kpc) of the Milky Way using data from Gaia, LAMOST and SDSS. Our method is built using an unsupervised learning algorithm, a self-organizing map, which trains a 2-D neural network to learn the topological structures of a data set from an n-D input space. Using a 4-D space of angular momentum and orbital energy, we identify three distinct groups corresponding to the Sagittarius, Orphan, and Cetus Streams. For the first time we are able to discover a northern counterpart to the Cetus stream. We test the robustness of this new detection using mock data and find that the significance is more than 5-sigma. We also find that the existing southern counterpart bifurcates into two clumps with different radial velocities. By exploiting the visualization power of StarGO, we attach MW globular clusters to the same trained neural network. The Sagittarius stream is found to have five related clusters, confirming recent literature studies, and the Cetus stream has one associated cluster, NGC 5824. This latter association has previously been postulated, but can only now be truly confirmed thanks to the high-precision Gaia proper motions and large numbers of stellar spectra from LAMOST. The large metallicity dispersion of the stream indicates that the progenitor cannot be a globular cluster. Given the mean metallicity of the stream, we propose that the stream is the result of a merger of a low-mass dwarf galaxy that hosted a large nuclear star cluster (NGC 5824). △ Less

Submitted 28 June, 2019; v1 submitted 14 February, 2019; originally announced February 2019.

Comments: 11 pages, 6 figures, accepted to ApJ

Fast and energetic AGN-driven outflows in simulated dwarf galaxies

Authors: Sophie Koudmani, Debora Sijacki, Martin A. Bourne, Matthew C. Smith

Abstract: The systematic analysis of optical large-scale surveys has revealed a population of dwarf galaxies hosting AGN, which have been confirmed by X-ray follow-up observations. Recently, the MaNGA survey identified six dwarf galaxies that appear to have an AGN that is preventing on-going star formation. It is therefore timely to study the physical properties of dwarf galaxies, in particular whether the… ▽ More The systematic analysis of optical large-scale surveys has revealed a population of dwarf galaxies hosting AGN, which have been confirmed by X-ray follow-up observations. Recently, the MaNGA survey identified six dwarf galaxies that appear to have an AGN that is preventing on-going star formation. It is therefore timely to study the physical properties of dwarf galaxies, in particular whether the presence of an AGN can affect their evolution. Using the moving mesh code AREPO, we have investigated different models of AGN activity, ranging from simple energy-driven spherical winds to collimated, mass-loaded, bipolar outflows in high resolution simulations of isolated dwarf galaxies hosting an active black hole. Our simulations also include a novel implementation of star formation and mechanical supernova (SN) feedback. We find that AGN outflows have a small but systematic effect on the central star formation rates (SFRs) for all set-ups explored, while substantial effects on the global SFR are only obtained with strong SNe and a sustained high-luminosity AGN with an isotropic wind. This suggests that AGN feedback in dwarf galaxies is unlikely to directly regulate their global SFRs. There is, however, a significant effect on outflow properties, which are notably enhanced by the AGN to much higher outflow temperatures and velocities, in agreement with kinematic signatures from the MaNGA survey. This indicates that AGN may play an indirect role in regulating the baryon cycle in dwarf galaxies by hindering cosmic gas inflows. △ Less

Submitted 8 January, 2019; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: Accepted for publication in MNRAS, 22 pages, 11 figures

Cosmological simulations of dwarfs: the need for ISM physics beyond SN feedback alone

Authors: Matthew C. Smith, Debora Sijacki, Sijing Shen

Abstract: The dominant feedback mechanism in low mass haloes is usually assumed to take the form of massive stars exploding as supernovae (SNe). We perform very high resolution cosmological zoom-in simulations of five dwarf galaxies to z = 4 with our mechanical SN feedback model. This delivers the correct amount of momentum corresponding to the stage of the SN remnant evolution resolved, and has been shown… ▽ More The dominant feedback mechanism in low mass haloes is usually assumed to take the form of massive stars exploding as supernovae (SNe). We perform very high resolution cosmological zoom-in simulations of five dwarf galaxies to z = 4 with our mechanical SN feedback model. This delivers the correct amount of momentum corresponding to the stage of the SN remnant evolution resolved, and has been shown to lead to realistic dwarf properties in isolated simulations. We find that in 4 out of our 5 simulated cosmological dwarfs, SN feedback has insufficient impact resulting in excessive stellar masses, extremely compact sizes and central super-solar stellar metallicities. The failure of the SN feedback in our dwarfs is physical in nature within our model and is the result of the build up of very dense gas in the early universe due to mergers and cosmic inflows prior to the first SN occurring. We demonstrate that our results are insensitive to resolution (provided that it is high enough), details of the (spatially uniform) UV background and reasonable alterations within our star formation prescription. We therefore conclude that the ability of SNe to regulate dwarf galaxy properties is dependent on other physical processes, such as turbulent pressure support, clustering and runaway of SN progenitors and other sources of stellar feedback. △ Less

Submitted 2 April, 2019; v1 submitted 11 July, 2018; originally announced July 2018.

Comments: 18 pages, 13 figures, accepted MNRAS

Journal ref: 2019 MNRAS 485, 3317

StarGO: A New Method to Identify the Galactic Origins of Halo Stars

Authors: Zhen Yuan, Jiang Chang, Projjwal Banerjee, Jiaxin Han, Xi Kang, M. C. Smith

Abstract: We develop a new method StarGO (Stars' Galactic Origin) to identify the galactic origins of halo stars using their kinematics. Our method is based on self-organizing map (SOM), which is one of the most popular unsupervised learning algorithms. StarGO combines SOM with a novel adaptive group identification algorithm with essentially no free parameters. In order to evaluate our model, we build a syn… ▽ More We develop a new method StarGO (Stars' Galactic Origin) to identify the galactic origins of halo stars using their kinematics. Our method is based on self-organizing map (SOM), which is one of the most popular unsupervised learning algorithms. StarGO combines SOM with a novel adaptive group identification algorithm with essentially no free parameters. In order to evaluate our model, we build a synthetic stellar halo from mergers of nine satellites in the Milky Way. We construct the mock catalogue by extracting a heliocentric volume of 10 kpc from our simulations and assigning expected observational uncertainties corresponding to bright stars from Gaia DR2 and LAMOST DR5. We compare the results from StarGO against that from a Friends-of-Friends (FoF) based method in the space of orbital energy and angular momentum. We show that StarGO is able to systematically identify more satellites and achieve higher number fraction of identified stars for most of the satellites within the extracted volume. When applied to data from Gaia DR2, StarGO will enable us to reveal the origins of the inner stellar halo in unprecedented detail. △ Less

Submitted 17 June, 2018; originally announced June 2018.

Comments: 11 pages, 7 figures, Accepted for publication in ApJ

The Lives of Stars: Insights From the TGAS-RAVE-LAMOST Dataset

Authors: John J. Vickers, Martin C. Smith

Abstract: In this paper we investigate how the chemical and kinematic properties of stars vary as a function of age. Using data from a variety of photometric, astrometric and spectroscopic surveys, we calculate the ages, phase space information and orbits for $\sim$125,000 stars covering a wide range of stellar parameters. We find indications that the inner regions of the disk reached high levels of enric… ▽ More In this paper we investigate how the chemical and kinematic properties of stars vary as a function of age. Using data from a variety of photometric, astrometric and spectroscopic surveys, we calculate the ages, phase space information and orbits for $\sim$125,000 stars covering a wide range of stellar parameters. We find indications that the inner regions of the disk reached high levels of enrichment early, while the outer regions were more substantially enriched in intermediate and recent epochs. We consider these enrichment histories through comparison of the ages of stars, their metallicities, and kinematic properties, such as their angular momentum in the solar neighborhood (which is a proxy for orbital radius). We calculate rates at which the velocity dispersions evolve, investigate the Oort constants for different aged populations (finding a slightly negative $\partial V_{C} / \partial R$ and $\partial V_{R} / \partial R$ for all ages, being most negative for the oldest stars), as well as examine the behavior of the velocity vertex deviation angle as a function of age (which we find to fall from $\sim$15 degrees for the 2 Gyr aged population to $\sim$6 degrees at around 6.5 Gyr of age, after which it remains unchanged). We find evidence for stellar churning, and find that the churned stars have a slightly younger age distribution than the rest of the data. △ Less

Submitted 6 May, 2018; originally announced May 2018.

Comments: 18 Pages, 14 Figures, Accepted ApJ

LAMOST DR1: Stellar parameters and chemical abundances with SP_Ace

Authors: C. Boeche, M. C. Smith, E. K. Grebel, J. Zhong, J. L. Hou, L. Chen, D. Stello

Abstract: We present a new analysis of the LAMOST DR1 survey spectral database performed with the code SP_Ace, which provides the derived stellar parameters T$_{\rm eff}$, log (g), [Fe/H], and [$α$/Fe] for 1,097,231 stellar objects. We tested the reliability of our results by comparing them to reference results from high spectral resolution surveys. The expected errors can be summarized as $\sim$120 K in T… ▽ More We present a new analysis of the LAMOST DR1 survey spectral database performed with the code SP_Ace, which provides the derived stellar parameters T$_{\rm eff}$, log (g), [Fe/H], and [$α$/Fe] for 1,097,231 stellar objects. We tested the reliability of our results by comparing them to reference results from high spectral resolution surveys. The expected errors can be summarized as $\sim$120 K in T$_{\rm eff}$, $\sim$0.2 in log (g), $\sim$0.15 dex in [Fe/H], and $\sim$0.1 dex in [$α$/Fe] for spectra with S/N$>$40, with some differences between dwarf and giant stars. SP_Ace provides error estimations consistent with the discrepancies observed between derived and reference parameters. Some systematic errors are identified and discussed. The resulting catalog is publicly available at the LAMOST and CDS websites. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 23 pages. Accepted for publication in AJ. Catalogue now available from http://dr1.lamost.org/doc/vac and will appear on the CDS website upon publication

The velocity ellipsoid in the Galactic disc using Gaia DR1

Authors: Borja Anguiano, Steven R. Majewski, Kenneth C. Freeman, Arik W. Mitschang, Martin C. Smith

Abstract: The stellar velocity ellipsoid of the solar neighbour is re-examined using intermediate-old mono-abundances stellar groups with high quality chemistry data together with parallaxes and proper motions from Gaia DR1. We find the average velocity dispersion values for the three space velocity components for the thin and thick disc of (σ_{U},σ_{V},σ_{W})_{thin} = (33 \pm 4, 28 \pm 2, 23 \pm 2) and (σ_… ▽ More The stellar velocity ellipsoid of the solar neighbour is re-examined using intermediate-old mono-abundances stellar groups with high quality chemistry data together with parallaxes and proper motions from Gaia DR1. We find the average velocity dispersion values for the three space velocity components for the thin and thick disc of (σ_{U},σ_{V},σ_{W})_{thin} = (33 \pm 4, 28 \pm 2, 23 \pm 2) and (σ_{U},σ_{V},σ_{W})_{thick} = (57 \pm 6, 38 \pm 5, 37 \pm 4) km s^{-1}, respectively. The mean values of the ratio between the semi-axes of the velocity ellipsoid for the thin disc are found to be, σ_{V}/σ_{U} = 0.70 \pm 0.13 and σ_{W}/σ_{U} is 0.64 \pm 0.08, while for the thick disc σ_{V}/σ_{U} = 0.67 \pm 0.11 and σ_{W}/σ_{U} is 0.66 \pm 0.11. Inputting these dispersions into the linear Strömberg relation for the thin disc groups, we find the Sun's velocity with respect to the LSR in Galactic rotation to be V_{\sun} = 13.9 \pm 3.4 km s^{-1}. A relation is found between the vertex deviation and the chemical abundances for the thin disc, ranging from -5 to +40^{\circ} as iron-abundance increases. For the thick disc we find a vertex deviation of l_{uv} \sim -15^{\circ}. The tilt angle (l_{uw}) in the U-W plane for the thin disc groups ranges from -10 to +15^\circ, but there is no evident relation between l_{uw} and the mean abundances. However we find a weak relation for l_{uw} as a function of iron abundances and α-elements for most of the groups in the thick disc, where the tilt angle decreases from -5 to -20^\circ when [Fe/H] decreases and [α/Fe] increases. The velocity anisotropy parameter is independent of the chemical group abundances and its value is nearly constant for both discs (β\sim 0.5), suggesting that the combined disc is dynamically relaxed. △ Less

Submitted 23 October, 2017; originally announced October 2017.

Comments: 13 pages, 13 figures, accepted for publication in MNRAS

Supernova feedback in numerical simulations of galaxy formation: separating physics from numerics

Authors: Matthew C. Smith, Debora Sijacki, Sijing Shen

Abstract: While feedback from massive stars exploding as supernovae (SNe) is thought to be one of the key ingredients regulating galaxy formation, theoretically it is still unclear how the available energy couples to the interstellar medium and how galactic scale outflows are launched. We present a novel implementation of six sub-grid SN feedback schemes in the moving-mesh code Arepo, including injections o… ▽ More While feedback from massive stars exploding as supernovae (SNe) is thought to be one of the key ingredients regulating galaxy formation, theoretically it is still unclear how the available energy couples to the interstellar medium and how galactic scale outflows are launched. We present a novel implementation of six sub-grid SN feedback schemes in the moving-mesh code Arepo, including injections of thermal and/or kinetic energy, two parametrizations of delayed cooling feedback and a `mechanical' feedback scheme that injects the correct amount of momentum depending on the relevant scale of the SN remnant resolved. All schemes make use of individually time-resolved SN events. Adopting isolated disk galaxy setups at different resolutions, with the highest resolution runs reasonably resolving the Sedov-Taylor phase of the SN, we aim to find a physically motivated scheme with as few tunable parameters as possible. As expected, simple injections of energy overcool at all but the highest resolution. Our delayed cooling schemes result in overstrong feedback, destroying the disk. The mechanical feedback scheme is efficient at suppressing star formation, agrees well with the Kennicutt-Schmidt relation and leads to converged star formation rates and galaxy morphologies with increasing resolution without fine tuning any parameters. However, we find it difficult to produce outflows with high enough mass loading factors at all but the highest resolution, indicating either that we have oversimplified the evolution of unresolved SN remnants, require other stellar feedback processes to be included, require a better star formation prescription or most likely some combination of these issues. △ Less

Submitted 26 April, 2018; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: 32 pages, 22 figures. Accepted for publication in MNRAS

Generalized Value Iteration Networks: Life Beyond Lattices

Authors: Sufeng Niu, Siheng Chen, Hanyu Guo, Colin Targonski, Melissa C. Smith, Jelena Kovačević

Abstract: In this paper, we introduce a generalized value iteration network (GVIN), which is an end-to-end neural network planning module. GVIN emulates the value iteration algorithm by using a novel graph convolution operator, which enables GVIN to learn and plan on irregular spatial graphs. We propose three novel differentiable kernels as graph convolution operators and show that the embedding based kerne… ▽ More In this paper, we introduce a generalized value iteration network (GVIN), which is an end-to-end neural network planning module. GVIN emulates the value iteration algorithm by using a novel graph convolution operator, which enables GVIN to learn and plan on irregular spatial graphs. We propose three novel differentiable kernels as graph convolution operators and show that the embedding based kernel achieves the best performance. We further propose episodic Q-learning, an improvement upon traditional n-step Q-learning that stabilizes training for networks that contain a planning module. Lastly, we evaluate GVIN on planning problems in 2D mazes, irregular graphs, and real-world street networks, showing that GVIN generalizes well for both arbitrary graphs and unseen graphs of larger scale and outperforms a naive generalization of VIN (discretizing a spatial graph into a 2D image). △ Less

Submitted 26 October, 2017; v1 submitted 7 June, 2017; originally announced June 2017.

Comments: 14 pages, conference

Identification of the white dwarf companion to millisecond pulsar J2317+1439

Authors: S. Dai, M. C. Smith, S. Wang, S. Okamoto, R. X. Xu, Y. L. Yue, J. F. Liu

Abstract: We report identification of the optical counterpart to the companion of the millisecond pulsar J2317+1439. At the timing position of the pulsar, we find an object with $g=22.96\pm0.05$, $r=22.86\pm0.04$ and $i=22.82\pm0.05$. The magnitudes and colors of the object are consistent with it being a white dwarf. By comparing with white dwarf cooling models, we estimate that it has a mass of… ▽ More We report identification of the optical counterpart to the companion of the millisecond pulsar J2317+1439. At the timing position of the pulsar, we find an object with $g=22.96\pm0.05$, $r=22.86\pm0.04$ and $i=22.82\pm0.05$. The magnitudes and colors of the object are consistent with it being a white dwarf. By comparing with white dwarf cooling models, we estimate that it has a mass of $0.39^{+0.13}_{-0.10}$ M$_{\odot}$, an effective temperature of $8077^{+550}_{-470}$ K and a cooling age of $10.9\pm0.3$ Gyr. Combining our results with published constraints on the orbital parameters obtained through pulsar timing, we estimate the pulsar mass to be $3.4^{+1.4}_{-1.1}$ M$_{\odot}$. Although the constraint on the pulsar mass is still weak, there is a significant possibility that the pulsar could be more massive than two solar mass. △ Less

Submitted 7 May, 2017; originally announced May 2017.

Comments: 7 pages, 6 figures, accepted for publication in ApJ

Red Runaways II: Low mass Hills stars in SDSS Stripe 82

Authors: Yanqiong Zhang, Martin C. Smith, Jeffrey L. Carlin

Abstract: Stars ejected from the Galactic centre can be used to place important constraints on the Milky Way potential. Since existing hypervelocity stars are too distant to accurately determine orbits, we have conducted a search for nearby candidates using full three-dimensional velocities. Since the efficacy of such studies are often hampered by deficiencies in proper motion catalogs, we have chosen to ut… ▽ More Stars ejected from the Galactic centre can be used to place important constraints on the Milky Way potential. Since existing hypervelocity stars are too distant to accurately determine orbits, we have conducted a search for nearby candidates using full three-dimensional velocities. Since the efficacy of such studies are often hampered by deficiencies in proper motion catalogs, we have chosen to utilize the reliable, high-precision SDSS Stripe 82 proper motion catalog. Although we do not find any candidates which have velocities in excess of the escape speed, we identify 226 stars on orbits that are consistent with Galactic centre ejection. This number is significantly larger than what we would expect for halo stars on radial orbits and cannot be explained by disk or bulge contamination. If we restrict ourselves to metal-rich stars, we find 29 candidates with [Fe/H] > -0.8 dex and 10 with [Fe/H] > -0.6 dex. Their metallicities are more consistent with what we expect for bulge ejecta, and so we believe these candidates are especially deserving of further study. We have supplemented this sample using our own radial velocities, developing an algorithm to use proper motions for optimizing candidate selection. This technique provides considerable improvement on the blind spectroscopic sample of SDSS, being able to identify candidates with an efficiency around 20 times better than a blind search. △ Less

Submitted 13 August, 2016; originally announced August 2016.

Comments: 13 pages, accepted for publication in ApJ

The Eruption of the Candidate Young Star ASASSN-15qi

Authors: Gregory J. Herczeg, Subo Dong, Benjamin J. Shappee, Ping Chen, Lynne A. Hillenbrand, Jessy Jose, Christopher S. Kochanek, Jose L. Prieto, K. Z. Stanek, Kyle Kaplan, Thomas W. -S. Holoien, Steve Mairs, Doug Johnstone, Michael Gully-Santiago, Zhaohuan Zhu, Martin C. Smith, David Bersier, Gijs D. Mulders, Alexei V. Filippenko, Kazuya Ayani, Joseph Brimacombe, Jonathan S. Brown, Michael Connelley, Jussi Harmanen, Ryosuke Ito , et al. (5 additional authors not shown)

Abstract: Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band… ▽ More Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling. △ Less

Submitted 15 August, 2016; v1 submitted 21 July, 2016; originally announced July 2016.

Comments: 23 pages; accepted by ApJ

The Gaia-ESO survey: Metal-rich bananas in the bulge

Authors: Angus A. Williams, N. W. Evans, Matthew Molloy, Georges Kordopatis, M. C. Smith, J. Shen, G. Gilmore, S. Randich, T. Bensby, P. Francois, S. E Koposov, A. Recio-Blanco, A. Bayo, G. Carraro, A. Casey, T. Costado, E. Franciosini, A. Hourihane, P. de Laverny, J. Lewis, K. Lind, L. Magrini, L. Monaco, L. Morbidelli, G. G Sacco , et al. (3 additional authors not shown)

Abstract: We analyse the kinematics of $\sim 2000$ giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region $-10^\circ \lesssim \ell \lesssim 10^\circ$ and $-11^\circ \lesssim b \lesssim -3^\circ$. We find distinct kinematic trends in the metal rich ($\mathrm{[M/H]}>0$) and metal poor ($\mathrm{[M/H]}<0$) stars in the data. The velocity dispersion of the metal-ric… ▽ More We analyse the kinematics of $\sim 2000$ giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region $-10^\circ \lesssim \ell \lesssim 10^\circ$ and $-11^\circ \lesssim b \lesssim -3^\circ$. We find distinct kinematic trends in the metal rich ($\mathrm{[M/H]}>0$) and metal poor ($\mathrm{[M/H]}<0$) stars in the data. The velocity dispersion of the metal-rich stars drops steeply with latitude, compared to a flat profile in the metal-poor stars, as has been seen previously. We argue that the metal-rich stars in this region are mostly on orbits that support the boxy-peanut shape of the bulge, which naturally explains the drop in their velocity dispersion profile with latitude. The metal rich stars also exhibit peaky features in their line-of-sight velocity histograms, particularly along the minor axis of the bulge. We propose that these features are due to stars on resonant orbits supporting the boxy-peanut bulge. This conjecture is strengthened through the comparison of the minor axis data with the velocity histograms of resonant orbits generated in simulations of buckled bars. The 'banana' or 2:1:2 orbits provide strongly bimodal histograms with narrow velocity peaks that resemble the Gaia-ESO metal-rich data. △ Less

Submitted 31 May, 2016; originally announced May 2016.

Comments: 6 pages, 8 figures, accepted to ApJ Letters

Kinematically Detected Halo Streams

Authors: Martin C. Smith

Abstract: Clues to the origins and evolution of our Galaxy can be found in the kinematics of stars around us. Remnants of accreted satellite galaxies produce over- densities in velocity-space, which can remain coherent for much longer than spatial over-densities. This chapter reviews a number of studies that have hunted for these accretion relics, both in the nearby solar-neighborhood and the more-distant s… ▽ More Clues to the origins and evolution of our Galaxy can be found in the kinematics of stars around us. Remnants of accreted satellite galaxies produce over- densities in velocity-space, which can remain coherent for much longer than spatial over-densities. This chapter reviews a number of studies that have hunted for these accretion relics, both in the nearby solar-neighborhood and the more-distant stellar halo. Many observational surveys have driven this field forwards, from early work with the Hipparcos mission, to contemporary surveys like RAVE & SDSS. This active field continues to flourish, providing many new discoveries, and will be revolutionised as the Gaia mission delivers precise proper motions for a billion stars in our Galaxy. △ Less

Submitted 7 March, 2016; originally announced March 2016.

Comments: 27 pages, 10 figures. Chapter from Springer ASSL Volume entitled "Tidal Streams in the Local Group and Beyond". Affluent readers may wish to purchase the full volume here: http://link.springer.com/book/10.1007/978-3-319-19336-6

Journal ref: From "Tidal Streams in the Local Group and Beyond", Astrophysics and Space Science Library, Volume 420. ISBN 978-3-319-19335-9. Springer International Publishing Switzerland, 2016, p. 113

Selecting M-giants with infra-red photometry: Distances, metallicities and the Sagittarius stream

Authors: Jing Li, Martin C. Smith, Jing Zhong, Jinliang Hou, Jeffrey L. Carlin, Heidi Jo Newberg, Chao Liu, Li Chen, Linlin Li, Zhengyi Shao, Emma Small, Hao Tian

Abstract: Using a spectroscopically confirmed sample of M-giants, M-dwarfs and quasars from the LAMOST survey, we assess how well WISE $\&$ 2MASS color-cuts can be used to select M-giant stars. The WISE bands are very efficient at separating M-giants from M-dwarfs and we present a simple classification that can produce a clean and relatively complete sample of M-giants. We derive a new photometric relation… ▽ More Using a spectroscopically confirmed sample of M-giants, M-dwarfs and quasars from the LAMOST survey, we assess how well WISE $\&$ 2MASS color-cuts can be used to select M-giant stars. The WISE bands are very efficient at separating M-giants from M-dwarfs and we present a simple classification that can produce a clean and relatively complete sample of M-giants. We derive a new photometric relation to estimate the metallicity for M-giants, calibrated using data from the APOGEE survey. We find a strong correlation between the $(W1-W2)$ color and $\rm [M/H]$, where almost all of the scatter is due to photometric uncertainties. We show that previous photometric distance relations, which are mostly based on stellar models, may be biased and devise a new empirical distance relation, investigating trends with metallicity and star formation history. Given these relations, we investigate the properties of M-giants in the Sagittarius stream. The offset in the orbital plane between the leading and trailing tails is reproduced and, by identifying distant M-giants in the direction of the Galactic anti-center, we confirm that the previously detected debris in the outer halo is the apocenter of the trailing tail. We also find tentative evidence supporting an existing overdensity near the leading tail in the Northern Galactic hemisphere, possibly an extension to the trailing tail (so-called Branch C). We have measured the metallicity distribution along the stream, finding a clear metallicity offset between the leading and trailing tails, in agreement with models for the stream formation. We include an online table of M-giants to facilitate further studies. △ Less

Submitted 14 April, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: 15 pages, 13 figures, accepted by ApJ. Minor issue fixed with Fig 12 and Section 5.3

Journal ref: ApJ 823 59 (2016)

LAMOST 1: A Disrupted Satellite in the Constellation Draco

Authors: John J. Vickers, Martin C. Smith, Yonghui Hou, Yuefei Wang, Yong Zhang

Abstract: Using LAMOST spectroscopic data, we find a strong signal of a comoving group of stars in the constellation of Draco. The group, observed near the apocenter of its orbit, is 2.6 kpc from the Sun with a metallicity of -0.64 dex. The system is observed as a streaming population of unknown provenance with mass of about 2.1E4 solar masses and an absolute V band magnitude of about -3.6. Its high metalli… ▽ More Using LAMOST spectroscopic data, we find a strong signal of a comoving group of stars in the constellation of Draco. The group, observed near the apocenter of its orbit, is 2.6 kpc from the Sun with a metallicity of -0.64 dex. The system is observed as a streaming population of unknown provenance with mass of about 2.1E4 solar masses and an absolute V band magnitude of about -3.6. Its high metallicity, diffuse physical structure, and eccentric orbit may indicate that the progenitor satellite was a globular cluster rather than a dwarf galaxy or an open cluster. △ Less

Submitted 16 December, 2015; originally announced December 2015.

Comments: 6 pages, 4 Figures, 1 Table, Accepted to ApJL

An extended view of the Pisces Overdensity from the SCUSS survey

Authors: J. D. Nie, M. C. Smith, V. Belokurov, X. H. Fan, Z. Fan, M. J. Irwin, Z. J. Jiang, Y. P. Jing, S. E. Koposov, M. Lesser, J. Ma, S. Y. Shen, J. L. Wang, Z. Y. Wu, T. M. Zhang, X. Zhou, Z. M. Zhou, H. Zou

Abstract: SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this adva… ▽ More SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system. △ Less

Submitted 7 September, 2015; v1 submitted 5 August, 2015; originally announced August 2015.

Comments: 8 pages, 4 figures, accepted by ApJ

Red Runaways: Hypervelocity Stars, Hills Ejecta and Other Outliers in the F-M Star Regime

Authors: John J. Vickers, Martin C. Smith, Eva K. Grebel

Abstract: In this paper we analyze a sample of metal-rich (>-0.8 dex) main sequence stars in the extended solar neighborhood, investigating kinematic outliers from the background population. The data, which are taken from the Sloan Digital Sky Survey, are kinematically profiled as a function of distance from the Galactic plane using full six dimensional phase space information. Each star is examined in the… ▽ More In this paper we analyze a sample of metal-rich (>-0.8 dex) main sequence stars in the extended solar neighborhood, investigating kinematic outliers from the background population. The data, which are taken from the Sloan Digital Sky Survey, are kinematically profiled as a function of distance from the Galactic plane using full six dimensional phase space information. Each star is examined in the context of these kinematic profiles and likelihoods are assigned to quantify whether a star matches the underlying profile. Since some of these stars are likely to have been ejected from the disc, we trace back their orbits in order to determine potential ejection radii. We find that objects with low probability (i.e. `outliers') are typically more metal poor, faster and, most importantly, have a tendency to originate from the inner Galaxy compared to the underlying population. We also compose a sample of stars with velocities exceeding the local escape velocity. Although we do not discount that our sample could be contaminated by objects with spurious proper motions, a number of stars appear to have been ejected from the disc with exceptionally high velocities. Some of these are consistent with being ejected from the spiral arms and hence are a rich resource for further study. Finally we look at objects whose orbits are consistent with them being ejected at high speeds from the Galactic center. Of these objects we find that one, J135855.65+552538.19, is inconsistent with halo, bulge and disk kinematics and could plausibly have been ejected from the Galactic nucleus via a Hills mechanism. △ Less

Submitted 9 June, 2015; originally announced June 2015.

Comments: 17 Pages, 12 Figures, Accepted to AJ

Estimation of distances to stars with stellar parameters from LAMOST

Authors: Jeffrey L. Carlin, Chao Liu, Heidi Jo Newberg, Timothy C. Beers, Li Chen, Licai Deng, Puragra Guhathakurta, Jinliang Hou, Yonghui Hou, Sebastien Lepine, Guangwei Li, A-Li Luo, Martin C. Smith, Yue Wu, Ming Yang, Brian Yanny, Haotong Zhang, Zheng Zheng

Abstract: We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star's absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Mul… ▽ More We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star's absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) survey. Because LAMOST obtains roughly 3000 stellar spectra simultaneously within each ~5-degree diameter "plate" that is observed, we can use the stellar parameters of the observed stars to account for the stellar luminosity function and target selection effects. This removes biasing assumptions about the underlying populations, both due to predictions of the luminosity function from stellar evolution modeling, and from Galactic models of stellar populations along each line of sight. Using calibration data of stars with known distances and stellar parameters, we show that our method recovers distances for most stars within ~20%, but with some systematic overestimation of distances to halo giants. We apply our code to the LAMOST database, and show that the current precision of LAMOST stellar parameters permits measurements of distances with ~40% error bars. This precision should improve as the LAMOST data pipelines continue to be refined. △ Less

Submitted 20 May, 2015; originally announced May 2015.

Comments: 11 pages, 12 figures; accepted for publication in AJ

Resonant Orbits and the High Velocity Peaks Towards the Bulge

Authors: Matthew Molloy, Martin C. Smith, N. Wyn Evans, Juntai Shen

Abstract: We extract the resonant orbits from an N-body bar that is a good representation of the Milky Way, using the method recently introduced by Molloy et al. (2015). By decomposing the bar into its constituent orbit families, we show that they are intimately connected to the boxy-peanut shape of the density. We highlight the imprint due solely to resonant orbits on the kinematic landscape towards the Ga… ▽ More We extract the resonant orbits from an N-body bar that is a good representation of the Milky Way, using the method recently introduced by Molloy et al. (2015). By decomposing the bar into its constituent orbit families, we show that they are intimately connected to the boxy-peanut shape of the density. We highlight the imprint due solely to resonant orbits on the kinematic landscape towards the Galactic centre. The resonant orbits are shown to have distinct kinematic features and may be used to explain the cold velocity peak seen in the APOGEE commissioning data (Nidever at al., 2012). We show that high velocity peaks are a natural consequence of the motions of stars in the 2:1 orbit family and that stars on other higher order resonances can contribute to the peaks. The locations of the peaks vary with bar angle and, with the tacit assumption that the observed peaks are due to the 2:1 family, we find that the locations of the high velocity peaks correspond to bar angles in the range 10 < theta_bar < 25 (deg). However, some important questions about the nature of the peaks remain, such as their apparent absence in other surveys of the Bulge and the deviations from symmetry between equivalent fields in the north and south. We show that the absence of a peak in surveys at higher latitudes is likely due to the combination of a less prominent peak and a lower number density of bar supporting orbits at these latitudes. △ Less

Submitted 19 August, 2015; v1 submitted 16 May, 2015; originally announced May 2015.

Comments: 7 Figures, 1 Table, Now includes figures & discussion of higher order resonances, Minor revisions to text throughout, Conclusions unchanged

Kinematics of the X-shaped Milky Way Bulge: Expectations from a Self-consistent N-body Model

Authors: Yujing Qin, Juntai Shen, Zhao-Yu Li, Shude Mao, Martin C. Smith, R. Michael Rich, Andrea Kunder, Chao Liu

Abstract: We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bu… ▽ More We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along l=0 degree, where the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage. △ Less

Submitted 20 April, 2015; v1 submitted 23 March, 2015; originally announced March 2015.

Comments: Minor changes to match the ApJ accepted version; 17 pages; emulateapj format. The electronic tables of our model result are available upon request

Gravitational Microlensing by Neutron Stars and Radio Pulsars: Event Rates, Timescale Distributions, and Mass Measurements

Authors: Shi Dai, Martin C. Smith, Mengxiang Lin, Youling Yue, George Hobbs, Renxin Xu

Abstract: We investigate properties of Galactic microlensing events in which a stellar object is lensed by a neutron star. For an all-sky photometric microlensing survey, we determine the number of lensing events caused by $\sim10^{5}$ potentially-observable radio pulsars to be $\sim0.2\ \rm{yr^{-1}}$ for $10^{10}$ background stellar sources. We expect a few detectable events per year for the same number of… ▽ More We investigate properties of Galactic microlensing events in which a stellar object is lensed by a neutron star. For an all-sky photometric microlensing survey, we determine the number of lensing events caused by $\sim10^{5}$ potentially-observable radio pulsars to be $\sim0.2\ \rm{yr^{-1}}$ for $10^{10}$ background stellar sources. We expect a few detectable events per year for the same number of background sources from an astrometric microlensing survey. We show that such a study could lead to precise measurements of radio pulsar masses. For instance, if a pulsar distance could be constrained through radio observations, then its mass would be determined with a precision of $\sim10\%$. We also investigate the time-scale distributions for neutron star events, finding that they are much shorter than had been previously thought. For photometric events towards the Galactic centre that last $\sim15$ days, around $7\%$ will have a neutron star lens. This fraction drops rapidly for longer time-scales. Away from the bulge region we find that neutron stars will contribute $\sim40\%$ of the events that last less than $\sim10$ days. These results are in contrast to earlier work which found that the maximum fraction of neutron star events would occur on time-scales of hundreds of days. △ Less

Submitted 2 July, 2015; v1 submitted 9 February, 2015; originally announced February 2015.

Comments: 10 pages, accepted for publication in ApJ. v2 updated to reflect change of title in proof stage

Journal ref: 2015 ApJ, 802, 120

An Investigation of the Absolute Proper Motions of the SCUSS Catalog

Authors: Xiyan Peng, Zhaoxiang Qi, Zhenyu Wu, Jun Ma, Cuihua Du, Xu Zhou, Yong Yu, Zhenghong Tang, Zhaoji Jiang, Hu Zou, Zhou Fan, Xiaohui Fan, Martin C. Smith, Linhua Jiang, Yipeng Jing, Mario G. Lattanzi, Brian J. Mclean, Michael Lesser, Jundan Nie, Shiyin Shen, Jiali Wang, Tianmeng Zhang, Zhimin Zhou, Songhu Wang

Abstract: Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-depend… ▽ More Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for μα cos δ and μδ, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for μα cos δ and μδ. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent. △ Less

Submitted 9 February, 2015; originally announced February 2015.

Comments: 8 pages, 10 figures ,PASP accepted

Nearby Low-Mass Hypervelocity Stars

Authors: Yanqiong Zhang, Martin C. Smith, Jeffrey L. Carlin

Abstract: Hypervelocity stars are those that have speeds exceeding the escape speed and are hence unbound from the Milky Way. We investigate a sample of low-mass hypervelocity candidates obtained using data from the high-precision SDSS Stripe 82 catalogue, which we have combined with spectroscopy from the 200-inch Hale Telescope at Palomar Observatory. We find four good candidates, but without metallicities… ▽ More Hypervelocity stars are those that have speeds exceeding the escape speed and are hence unbound from the Milky Way. We investigate a sample of low-mass hypervelocity candidates obtained using data from the high-precision SDSS Stripe 82 catalogue, which we have combined with spectroscopy from the 200-inch Hale Telescope at Palomar Observatory. We find four good candidates, but without metallicities it is difficult to pin-down their distances and therefore total velocities. Our best candidate has a significant likelihood that it is escaping the Milky Way for a wide-range of metallicities. △ Less

Submitted 30 January, 2015; originally announced January 2015.

Comments: 5 pages; Contribution to proceedings for "The Milky Way Unravelled by Gaia" conference, Barcelona, Dec 2015