Showing 1–50 of 57 results for author: Campbell, D

A survey for variable young stars with small telescopes: VIII -- Properties of 1687 Gaia selected members in 21 nearby clusters

Authors: Dirk Froebrich, Aleks Scholz, Justyn Campbell-White, Siegfried Vanaverbeke, Carys Herbert, Jochen Eislöffel, Thomas Urtly, Timothy P. Long, Ivan L. Walton, Klaas Wiersema, Nick J. Quinn, Tony Rodda, Juan-Luis González-Carballo, Mario Morales Aimar, Rafael Castillo García, Francisco C. Soldán Alfaro, Faustino García de la Cuesta, Domenico Licchelli, Alex Escartin Perez, José Luis Salto González, Marc Deldem, Stephen R. L. Futcher, Tim Nelson, Shawn Dvorak, Dawid Moździerski , et al. (38 additional authors not shown)

Abstract: The Hunting Outbursting Young Stars (HOYS) project performs long-term, optical, multi-filter, high cadence monitoring of 25 nearby young clusters and star forming regions. Utilising Gaia DR3 data we have identified about 17000 potential young stellar members in 45 coherent astrometric groups in these fields. Twenty one of them are clear young groups or clusters of stars within one kiloparsec and t… ▽ More The Hunting Outbursting Young Stars (HOYS) project performs long-term, optical, multi-filter, high cadence monitoring of 25 nearby young clusters and star forming regions. Utilising Gaia DR3 data we have identified about 17000 potential young stellar members in 45 coherent astrometric groups in these fields. Twenty one of them are clear young groups or clusters of stars within one kiloparsec and they contain 9143 Gaia selected potential members. The cluster distances, proper motions and membership numbers are determined. We analyse long term (about 7yr) V, R, and I-band light curves from HOYS for 1687 of the potential cluster members. One quarter of the stars are variable in all three optical filters, and two thirds of these have light curves that are symmetric around the mean. Light curves affected by obscuration from circumstellar materials are more common than those affected by accretion bursts, by a factor of 2-4. The variability fraction in the clusters ranges from 10 to almost 100 percent, and correlates positively with the fraction of stars with detectable inner disks, indicating that a lot of variability is driven by the disk. About one in six variables shows detectable periodicity, mostly caused by magnetic spots. Two thirds of the periodic variables with disk excess emission are slow rotators, and amongst the stars without disk excess two thirds are fast rotators - in agreement with rotation being slowed down by the presence of a disk. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: accepted for publication in MNRAS, 1 table, 9 figures

An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations

Authors: Nicholas Van Alfen, Duncan Campbell, Jonathan Blazek, C. Danielle Leonard, Francois Lanusse, Andrew Hearin, Rachel Mandelbaum, The LSST Dark Energy Science Collaboration

Abstract: We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce n… ▽ More We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the IllustrisTNG suite. We create an HOD mock with TNG-like correlations using our method, constraining the associated IA model parameters, with the $χ^2_{\rm dof}$ between our model's correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation--position and orientation--orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to $0.1~h^{-1}{\rm Mpc}$. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime. △ Less

Submitted 3 June, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 19 pages, 13 figures, 3 tables, accepted for publication by The Open Journal of Astrophysics, code available at https://github.com/astropy/halotools

Journal ref: Open Journal of Astrophysics, 7 (2024)

A survey for variable young stars with small telescopes: VI -- Analysis of the outbursting Be stars NSW284, Gaia19eyy, and VES263

Authors: Dirk Froebrich, Lynne A. Hillenbrand, Carys Herbert, Kishalay De, Jochen Eislöffel, Justyn Campbell-White, Ruhee Kahar, Franz-Josef Hambsch, Thomas Urtly, Adam Popowicz, Krzysztof Bernacki, Andrzej Malcher, Slawomir Lasota, Jerzy Fiolka, Piotr Jozwik-Wabik, Franky Dubois, Ludwig Logie, Steve Rau, Mark Phillips, George Fleming, Rafael Gonzalez Farfán, Francisco C. Soldán Alfaro, Tim Nelson, Stephen R. L. Futcher, Samantha M. Rolfe , et al. (22 additional authors not shown)

Abstract: This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterised as rare, truly regularly outbursting Be stars. We inte… ▽ More This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterised as rare, truly regularly outbursting Be stars. We interpret the photometric data within a framework for modelling light curve morphology, and find that the models correctly predict the burst shapes, including their larger amplitudes and later peaks towards longer wavelengths. We are thus able to infer the start and end times of mass loading into the circumstellar disks of these stars. The disk sizes are typically 3-6 times the areas of the central star. The disk temperatures are ~40%, and the disk luminosities are ~10% of those of the central Be star, respectively. The available spectroscopy is consistent with inside-out evolution of the disk. Higher excitation lines have larger velocity widths in their double-horned shaped emission profiles. Our observations and analysis support the decretion disk model for outbursting Be stars. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted for publication by MNRAS. 20 pages, 11 figures, 2 tables

The Swansong of the Galactic Center Source X7: An Extreme Example of Tidal Evolution near the Supermassive Black Hole

Authors: Anna Ciurlo, Randall D. Campbell, Mark R. Morris, Tuan Do, Andrea M. Ghez, Eric E. Becklin, Rory O. Bentley, Devin S. Chu, Abhimat K. Gautam, Yash A. Gursahani, Aurelien Hees, Kelly Kosmo O'Neil, Jessica R. Lu, Gregory D. Martinez, Smadar Naoz, Shoko Sakai, Rainer Schoedel

Abstract: We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital p… ▽ More We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital parameters and quantitative characterization of the evolution of its morphology and mass. We find that the leading edge of X7 appears to be on a mildly eccentric (e~0.3), relatively short-period (170 years) orbit and is headed towards periapse passage, estimated to occur in ~2036. Furthermore, our kinematic measurements rule out the earlier suggestion that X7 is associated with the stellar source S0-73 or with any other point source that has overlapped with X7 during our monitoring period. Over the course of our observations, X7 has (1) become more elongated, with a current length-to-width ratio of 9, (2) maintained a very consistent long-axis orientation (position angle of 50 deg), (3) inverted its radial velocity differential from tip to tail from -50 to +80 km/sec, and (4) sustained its total brightness (12.8 Lp magnitudes at the leading edge) and color temperature (425 K), which suggest a constant mass of ~50 MEarth. We present a simple model showing that these results are compatible with the expected effect of tidal forces exerted on it by the central black hole and we propose that X7 is the gas and dust recently ejected from a grazing collision in a binary system. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: accepted by ApJ, 20 pages, 16 figures

Galaxies on graph neural networks: towards robust synthetic galaxy catalogs with deep generative models

Authors: Yesukhei Jagvaral, Francois Lanusse, Sukhdeep Singh, Rachel Mandelbaum, Siamak Ravanbakhsh, Duncan Campbell

Abstract: The future astronomical imaging surveys are set to provide precise constraints on cosmological parameters, such as dark energy. However, production of synthetic data for these surveys, to test and validate analysis methods, suffers from a very high computational cost. In particular, generating mock galaxy catalogs at sufficiently large volume and high resolution will soon become computationally un… ▽ More The future astronomical imaging surveys are set to provide precise constraints on cosmological parameters, such as dark energy. However, production of synthetic data for these surveys, to test and validate analysis methods, suffers from a very high computational cost. In particular, generating mock galaxy catalogs at sufficiently large volume and high resolution will soon become computationally unreachable. In this paper, we address this problem with a Deep Generative Model to create robust mock galaxy catalogs that may be used to test and develop the analysis pipelines of future weak lensing surveys. We build our model on a custom built Graph Convolutional Networks, by placing each galaxy on a graph node and then connecting the graphs within each gravitationally bound system. We train our model on a cosmological simulation with realistic galaxy populations to capture the 2D and 3D orientations of galaxies. The samples from the model exhibit comparable statistical properties to those in the simulations. To the best of our knowledge, this is the first instance of a generative model on graphs in an astrophysical/cosmological context. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Comments: Accepted as extended abstract at ICML 2022 Workshop on Machine Learning for Astrophysics. Condensed version of arXiv:2204.07077

Galaxies and Halos on Graph Neural Networks: Deep Generative Modeling Scalar and Vector Quantities for Intrinsic Alignment

Authors: Yesukhei Jagvaral, François Lanusse, Sukhdeep Singh, Rachel Mandelbaum, Siamak Ravanbakhsh, Duncan Campbell

Abstract: In order to prepare for the upcoming wide-field cosmological surveys, large simulations of the Universe with realistic galaxy populations are required. In particular, the tendency of galaxies to naturally align towards overdensities, an effect called intrinsic alignments (IA), can be a major source of systematics in the weak lensing analysis. As the details of galaxy formation and evolution releva… ▽ More In order to prepare for the upcoming wide-field cosmological surveys, large simulations of the Universe with realistic galaxy populations are required. In particular, the tendency of galaxies to naturally align towards overdensities, an effect called intrinsic alignments (IA), can be a major source of systematics in the weak lensing analysis. As the details of galaxy formation and evolution relevant to IA cannot be simulated in practice on such volumes, we propose as an alternative a Deep Generative Model. This model is trained on the IllustrisTNG-100 simulation and is capable of sampling the orientations of a population of galaxies so as to recover the correct alignments. In our approach, we model the cosmic web as a set of graphs, where the graphs are constructed for each halo, and galaxy orientations as a signal on those graphs. The generative model is implemented on a Generative Adversarial Network architecture and uses specifically designed Graph-Convolutional Networks sensitive to the relative 3D positions of the vertices. Given (sub)halo masses and tidal fields, the model is able to learn and predict scalar features such as galaxy and dark matter subhalo shapes; and more importantly, vector features such as the 3D orientation of the major axis of the ellipsoid and the complex 2D ellipticities. For correlations of 3D orientations the model is in good quantitative agreement with the measured values from the simulation, except for at very small and transition scales. For correlations of 2D ellipticities, the model is in good quantitative agreement with the measured values from the simulation on all scales. Additionally, the model is able to capture the dependence of IA on mass, morphological type and central/satellite type. △ Less

Submitted 22 July, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 15 pages, 7 figures (+2 figures in Appendix), matches published version at MNRAS

Optical and NIR data and modelling of nova V5668 Sgr

Authors: L. Takeda, M. Diaz, R. D. Campbell, J. E. Lyke, S. S. Lawrence, J. D. Linford, K. V. Sokolovsky

Abstract: We present HST optical images, Keck-OSIRIS NIR IFS data cubes and Keck-NIRC2 NIR images of nova V5668 Sgr from 2016 to 2019. The observations indicate enhanced emission at the polar caps and equatorial torus for low ionization lines, and enhanced high ionization emission lines only at the polar caps. The radial velocities are compatible with a homogeneous expansion velocity of v=590 km s$^{-1}$ an… ▽ More We present HST optical images, Keck-OSIRIS NIR IFS data cubes and Keck-NIRC2 NIR images of nova V5668 Sgr from 2016 to 2019. The observations indicate enhanced emission at the polar caps and equatorial torus for low ionization lines, and enhanced high ionization emission lines only at the polar caps. The radial velocities are compatible with a homogeneous expansion velocity of v=590 km s$^{-1}$ and a system inclination angle of 24$^o$. These values were used to estimate an expansion parallax distance of 1200 $\pm$ 400 pc. The NIRC2 data indicate the presence of dust in 2016 and 2017, but no dust emission could be detected in 2019. The observational data were used for assembling 3D photoionization models of the ejecta. The model results indicate that the central source has a temperature of $1.88\times10^{5}$ K and a luminosity of $1.6\times10^{35}$ erg s$^{-1}$ in August of 2017 (2.4 years post eruption), and that the shell has a mass of $6.3\times10^{-5}$ M$_{\odot}$. The models also suggest an anisotropy of the ionizing flux, possibly by the contribution from a luminous accretion disc. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Validating Synthetic Galaxy Catalogs for Dark Energy Science in the LSST Era

Authors: Eve Kovacs, Yao-Yuan Mao, Michel Aguena, Anita Bahmanyar, Adam Broussard, James Butler, Duncan Campbell, Chihway Chang, Shenming Fu, Katrin Heitmann, Danila Korytov, François Lanusse, Patricia Larsen, Rachel Mandelbaum, Christopher B. Morrison, Constantin Payerne, Marina Ricci, Eli Rykoff, F. Javier Sánchez, Ignacio Sevilla-Noarbe, Melanie Simet, Chun-Hao To, Vinu Vikraman, Rongpu Zhou, Camille Avestruz , et al. (14 additional authors not shown)

Abstract: Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation proc… ▽ More Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation procedure that ensures that the simulated galaxy properties faithfully capture observations and delivers an assessment of the level of realism attained by the catalog. We present here a suite of validation tests that have been developed by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). We discuss how the inclusion of each test is driven by the scientific targets for static ground-based dark energy science and by the availability of suitable validation data. The validation criteria that are used to assess the performance of a catalog are flexible and depend on the science goals. We illustrate the utility of this suite by showing examples for the validation of cosmoDC2, the extragalactic catalog recently released for the LSST DESC second Data Challenge. △ Less

Submitted 13 January, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 46 pages, 33 figures

Probabilistic model for dynamic galaxy decomposition

Authors: Yesukhei Jagvaral, Duncan Campbell, Rachel Mandelbaum, Markus Michael Rau

Abstract: In the era of precision cosmology and ever-improving cosmological simulations, a better understanding of different galaxy components such as bulges and discs will give us new insight into galactic formation and evolution. Based on the fact that the stellar populations of the constituent components of galaxies differ by their dynamical properties, we develop two simple models for galaxy decompositi… ▽ More In the era of precision cosmology and ever-improving cosmological simulations, a better understanding of different galaxy components such as bulges and discs will give us new insight into galactic formation and evolution. Based on the fact that the stellar populations of the constituent components of galaxies differ by their dynamical properties, we develop two simple models for galaxy decomposition using the IllustrisTNG cosmological hydrodynamical simulation. The first model uses a single dynamical parameter and can distinguish 4 components: thin disc, thick disc, counter-rotating disc and bulge. The second model uses one more dynamical parameter, was defined in a probabilistic manner, and distinguishes two components: bulge and disc. The number fraction of disc-dominated galaxies at a given stellar mass obtained by our models agrees well with observations for masses exceeding $ \log_{10}(M_*/M_\odot)=10$. Sérsic indices and half-mass radii for the bulge components agree well with those for real galaxies. The mode of the distribution of Sérsic indices for the disc components is at the expected value of $n=1$. However, disc half-mass radii are smaller than those for real galaxies, in accordance with previous findings that the IllustrisTNG simulation produces undersized discs. △ Less

Submitted 5 May, 2021; originally announced May 2021.

Comments: 15 pages, 12 figures

Spin state and moment of inertia of Venus

Authors: Jean-Luc Margot, Donald B. Campbell, Jon D. Giorgini, Joseph S. Jao, Lawrence G. Snedeker, Frank D. Ghigo, Amber Bonsall

Abstract: Fundamental properties of the planet Venus, such as its internal mass distribution and variations in length of day, have remained unknown. We used Earth-based observations of radar speckles tied to the rotation of Venus obtained in 2006-2020 to measure its spin axis orientation, spin precession rate, moment of inertia, and length-of-day variations. Venus is tilted by 2.6392 $\pm$ 0.0008 degrees (… ▽ More Fundamental properties of the planet Venus, such as its internal mass distribution and variations in length of day, have remained unknown. We used Earth-based observations of radar speckles tied to the rotation of Venus obtained in 2006-2020 to measure its spin axis orientation, spin precession rate, moment of inertia, and length-of-day variations. Venus is tilted by 2.6392 $\pm$ 0.0008 degrees ($1σ$) with respect to its orbital plane. The spin axis precesses at a rate of 44.58 $\pm$ 3.3 arcseconds per year ($1σ$), which gives a normalized moment of inertia of 0.337 $\pm$ 0.024 and yields a rough estimate of the size of the core. The average sidereal day on Venus in the 2006-2020 interval is 243.0226 $\pm$ 0.0013 Earth days ($1σ$). The spin period of the solid planet exhibits variations of 61 ppm ($\sim$20 minutes) with a possible diurnal or semidiurnal forcing. The length-of-day variations imply that changes in atmospheric angular momentum of at least $\sim$4% are transferred to the solid planet. △ Less

Submitted 30 March, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: 20 pages, 7 figures, supplementary information. Submitted to Nature Astronomy on October 14, 2020

The Future Of The Arecibo Observatory: The Next Generation Arecibo Telescope

Authors: D. Anish Roshi, N. Aponte, E. Araya, H. Arce, L. A. Baker, W. Baan, T. M. Becker, J. K. Breakall, R. G. Brown, C. G. M. Brum, M. Busch, D. B. Campbell, T. Cohen, F. Cordova, J. S. Deneva, M. Devogele, T. Dolch, F. O. Fernandez-Rodriguez, T. Ghosh, P. F. Goldsmith, L. I. Gurvits, M. Haynes, C. Heiles, J. W. T. Hessel, D. Hickson , et al. (49 additional authors not shown)

Abstract: The Arecibo Observatory (AO) is a multidisciplinary research and education facility that is recognized worldwide as a leading facility in astronomy, planetary, and atmospheric and space sciences. AO's cornerstone research instrument was the 305-m William E. Gordon telescope. On December 1, 2020, the 305-m telescope collapsed and was irreparably damaged. In the three weeks following the collapse, A… ▽ More The Arecibo Observatory (AO) is a multidisciplinary research and education facility that is recognized worldwide as a leading facility in astronomy, planetary, and atmospheric and space sciences. AO's cornerstone research instrument was the 305-m William E. Gordon telescope. On December 1, 2020, the 305-m telescope collapsed and was irreparably damaged. In the three weeks following the collapse, AO's scientific and engineering staff and the AO users community initiated extensive discussions on the future of the observatory. The community is in overwhelming agreement that there is a need to build an enhanced, next-generation radar-radio telescope at the AO site. From these discussions, we established the set of science requirements the new facility should enable. These requirements can be summarized briefly as: 5 MW of continuous wave transmitter power at 2 - 6 GHz, 10 MW of peak transmitter power at 430 MHz (also at 220MHz under consideration), zenith angle coverage 0 to 48 deg, frequency coverage 0.2 to 30 GHz and increased Field-of-View. These requirements determine the unique specifications of the proposed new instrument. The telescope design concept we suggest consists of a compact array of fixed dishes on a tiltable, plate-like structure with a collecting area equivalent to a 300m dish. This concept, referred to as the Next Generation Arecibo Telescope (NGAT), meets all of the desired specifications and provides significant new science capabilities to all three research groups at AO. This whitepaper presents a sample of the wide variety of the science that can be achieved with the NGAT, the details of the telescope design concept and the need for the new telescope to be located at the AO site. We also discuss other AO science activities that interlock with the NGAT in the white paper. △ Less

Submitted 1 April, 2021; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 82 pages (executive summary 10 pages), 21 figures, Arecibo observatory white paper (Updated with the complete author list and minor edits)

Upper Limit on Brackett-gamma Emission from the Immediate Accretion Flow onto the Galactic Black Hole

Authors: Anna Ciurlo, Mark R. Morris, Randall D. Campbell, Andrea M. Ghez, Tuan Do, Devin S. Chu

Abstract: We present the first observational constraint on the Br-gamma recombination line emission associated with the supermassive black hole at the center of our Galaxy, known as Sgr A*. By combining 13 years of data with the Adaptive Optics fed integral field spectrograph OSIRIS at the W. M. Keck Observatory obtained as part of the Galactic Center Orbits Initiative, we extract the near-infrared spectrum… ▽ More We present the first observational constraint on the Br-gamma recombination line emission associated with the supermassive black hole at the center of our Galaxy, known as Sgr A*. By combining 13 years of data with the Adaptive Optics fed integral field spectrograph OSIRIS at the W. M. Keck Observatory obtained as part of the Galactic Center Orbits Initiative, we extract the near-infrared spectrum within ~0.2'' of the black hole and we derive an upper limit on the Br-gamma flux. The aperture was set to match the size of the disk-like structure that was recently reported based on millimeter-wave ALMA observations of the hydrogen recombination line, H30-alpha. Our stringent upper limit is at least a factor of 80 (and up to a factor of 245) below what would be expected from the ALMA measurements and strongly constrains possible interpretation of emission from this highly under-luminous supermassive black hole. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: Accepted by ApJ

A small actively-controlled high-resolution spectrograph based on off-the-shelf components

Authors: H. R. A. Jones, W. E. Martin, G. Anglada-Escudé, R. Errmann, D. A. Campbell, C. Baker, C. Boonsri, P. Choochalerm

Abstract: We present the design and testing of a prototype in-plane echelle spectrograph based on an actively controlled fibre-fed double-pass design. This system aims to be small and efficient with the minimum number of optical surfaces - currently a collimator/camera lens, cross-dispersing prism, grating and a reflector to send light to the detector. It is built from catalogue optical components and has d… ▽ More We present the design and testing of a prototype in-plane echelle spectrograph based on an actively controlled fibre-fed double-pass design. This system aims to be small and efficient with the minimum number of optical surfaces - currently a collimator/camera lens, cross-dispersing prism, grating and a reflector to send light to the detector. It is built from catalogue optical components and has dimensions of approximately 20x30 cm. It works in the optical regime with a resolution of >70,000. The spectrograph is fed by a bifurcated fibre with one fibre to a telescope and the other used to provide simultaneous Thorium Argon light illumination for wavelength calibration. The positions of the arc lines on the detector are processed in real time and commercial auto-guiding software is used to treat the positions of the arc lines as guide stars. The guiding software sends any required adjustments to mechanical piezo-electric actuators which move the mirror sending light to the camera removing any drift in the position of the arc lines. The current configuration using an sCMOS detector provides a precision of 3.5 milli-pixels equivalent to 4 m/s in a standard laboratory environment. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 21 pages, 9 figures, 2 tables, accepted PASP

Void Galaxies Follow a Distinct Evolutionary Path in the Environmental COntext Catalog

Authors: Jonathan Florez, Andreas A. Berlind, Sheila J. Kannappan, David V. Stark, Kathleen D. Eckert, Victor F. Calderon, Amanda J. Moffett, Duncan Campbell, Manodeep Sinha

Abstract: We measure the environmental dependence, where environment is defined by the distance to the third nearest neighbor, of multiple galaxy properties inside the Environmental COntext (ECO) catalog. We focus primarily on void galaxies, which we define as the $10 \%$ of galaxies having the lowest local density. We compare the properties of void and non-void galaxies: baryonic mass, color, fractional st… ▽ More We measure the environmental dependence, where environment is defined by the distance to the third nearest neighbor, of multiple galaxy properties inside the Environmental COntext (ECO) catalog. We focus primarily on void galaxies, which we define as the $10 \%$ of galaxies having the lowest local density. We compare the properties of void and non-void galaxies: baryonic mass, color, fractional stellar mass growth rate (FSMGR), morphology, and gas-to-stellar-mass ratio (estimated from a combination of HI data and photometric gas fractions calibrated with the RESOLVE survey). Our void galaxies typically have lower baryonic masses than galaxies in denser environments, and they display the properties expected of a lower mass population: they have more late-types, are bluer, have higher FSMGR, and are more gas rich. We control for baryonic mass and investigate the extent to which void galaxies are different at fixed mass. Void galaxies are bluer, more gas-rich, and more star forming at fixed mass than non-void galaxies, which is a possible signature of galaxy assembly bias. Furthermore, we show that these trends persist even at fixed mass and morphology, and we find that voids host a distinct population of early-types that are bluer and more star-forming than the typical red and quenched early-types. In addition to these empirical observational results, we also present theoretical results from mock catalogs with built-in galaxy assembly bias. We show that a simple matching of galaxy properties to (sub)halo properties, such as mass and age, can recover the observed environmental trends in ECO galaxies. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 18 pages, 11 figures. Accepted for publication in ApJ

A population of dust-enshrouded objects orbiting the Galactic black hole

Authors: Anna Ciurlo, Randall D. Campbell, Mark R. Morris, Tuan Do, Andrea M. Ghez, Aurelien Hees, Breann N. Sitarski, Kelly Kosmo O'Neil, Devin S. Chu, Gregory D. Martinez, Smadar Naoz, Alexander P. Stephan

Abstract: The central 0.1 parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A*, a cluster of young, massive stars (the S stars) and various gaseous features. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the… ▽ More The central 0.1 parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A*, a cluster of young, massive stars (the S stars) and various gaseous features. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the order of 100 astronomical units, except at periapse, where the tidal interaction with the black hole stretches them along the orbit) and they show both thermal dust emission and line emission from ionized gas. G1 and G2 have generated attention because they appear to be tidally interacting with the supermassive Galactic black hole, possibly enhancing its accretion activity. No broad consensus has yet been reached concerning their nature: the G objects show the characteristics of gas and dust clouds but display the dynamical properties of stellar-mass objects. Here we report observations of four additional G objects, all lying within 0.04 parsecs of the black hole and forming a class that is probably unique to this environment. The widely varying orbits derived for the six G objects demonstrate that they were commonly but separately formed. △ Less

Submitted 22 January, 2020; originally announced January 2020.

Journal ref: Nature volume 577, pages337-340(2020)

Generating Synthetic Cosmological Data with GalSampler

Authors: Andrew Hearin, Danila Korytov, Eve Kovacs, Andrew Benson, Han Aung, Christopher Bradshaw, Duncan Campbell

Abstract: As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-anal… ▽ More As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: 12 pages, 5 figures

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

Authors: Danila Korytov, Andrew Hearin, Eve Kovacs, Patricia Larsen, Esteban Rangel, Joseph Hollowed, Andrew J. Benson, Katrin Heitmann, Yao-Yuan Mao, Anita Bahmanyar, Chihway Chang, Duncan Campbell, Joseph Derose, Hal Finkel, Nicholas Frontiere, Eric Gawiser, Salman Habib, Benjamin Joachimi, François Lanusse, Nan Li, Rachel Mandelbaum, Christopher Morrison, Jeffrey A. Newman, Adrian Pope, Eli Rykoff , et al. (5 additional authors not shown)

Abstract: This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `… ▽ More This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data. △ Less

Submitted 27 July, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 27 pages, 17 figures, submitted to APJS

Astro2020 Activities and Projects White Paper: Arecibo Observatory in the Next Decade

Authors: D. Anish Roshi, L. D. Anderson, E. Araya, D. Balser, W. Brisken, C. Brum, D. Campbell, S. Chatterjee, E. Churchwell, J. Condon, J. Cordes, F. Cordova, Y. Fernandez, J. Gago, T. Ghosh, P. F. Goldsmith, C. Heiles, D. Hickson, B. Jeffs, K. M. Jones, J. Lautenbach, B. M. Lewis, R. S. Lynch, P. K. Manoharan, S. Marshall , et al. (20 additional authors not shown)

Abstract: The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a… ▽ More The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a universal backend and (d) augment the VLBI facility by integrating the 12m telescope for phase referencing. These upgrades to the Arecibo telescope are critical to keep the national facility in the forefront of research in radio astronomy while maintaining its dominance in radar studies of near-Earth asteroids, planets and satellites. In the next decade, the Arecibo telescope will play a synergistic role with the upcoming facilities such as ngVLA, SKA and the now commissioned FAST telescope. Further, the observatory will be actively engaged in mentoring and training programs for students from a diverse background. △ Less

Submitted 13 July, 2019; originally announced July 2019.

Comments: Submitted as an Activities and Projects White Paper for the Astro2020 decadal survey (13 pages, 3 figures, 1 table)

Efficient Follow-Up of Exoplanet Transits Using Small Telescopes

Authors: Peter Beck, Luke Robson, Mark Gallaway, Hugh R. A. Jones, David Campbell

Abstract: This paper is to introduce an online tool for the prediction of exoplanet transit light curves. Small telescopes can readily capture exoplanet transits under good weather conditions when the combination of a bright star and a large transiting exoplanet results in a significant depth of transit. However, in reality there are many considerations that need to be made in order to obtain useful measure… ▽ More This paper is to introduce an online tool for the prediction of exoplanet transit light curves. Small telescopes can readily capture exoplanet transits under good weather conditions when the combination of a bright star and a large transiting exoplanet results in a significant depth of transit. However, in reality there are many considerations that need to be made in order to obtain useful measurements. This paper and accompanying website layout a procedure based on time series differential photometry that has been successfully employed using 0.4m aperture telescopes to predict the expected precision for a whole light curve. This enables robust planning to decide whether the observation of a particular exoplanet transit should be attempted and in particular to be able to readily see when it should not to be attempted. This may result in a significant increase in the number of transit observations captured by non-specialists. The technique and website are also appropriate for planning a variety of variable star observations where a prediction of the light curve can be made. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 11 pages, 5 figures, to be published in PASP, comments welcome

How to Optimally Constrain Galaxy Assembly Bias: Supplement Projected Correlation Functions with Count-in-cells Statistics

Authors: Kuan Wang, Yao-Yuan Mao, Andrew R. Zentner, Frank C. van den Bosch, Johannes U. Lange, Chad M. Schafer, Antonia S. Villarreal, Andrew P. Hearin, Duncan Campbell

Abstract: Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real U… ▽ More Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $ΔΣ(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+ΔΣ(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $ΔΣ(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection. △ Less

Submitted 6 September, 2019; v1 submitted 22 March, 2019; originally announced March 2019.

Comments: Figures 3 and 4 show the main results. Published in Monthly Notices of the Royal Astronomical Society

Characterizing and Improving the Data Reduction Pipeline for the Keck OSIRIS Integral Field Spectrograph

Authors: Kelly E. Lockhart, Tuan Do, James E. Larkin, Anna Boehle, Randy D. Campbell, Samantha Chappell, Devin Chu, Anna Ciurlo, Maren Cosens, Michael P. Fitzgerald, Andrea Ghez, Jessica R. Lu, Jim E. Lyke, Etsuko Mieda, Alexander R. Rudy, Andrey Vayner, Gregory Walth, Shelley A. Wright

Abstract: OSIRIS is a near-infrared (1.0--2.4 $μ$m) integral field spectrograph operating behind the adaptive optics system at Keck Observatory, and is one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics (LGS AO) papers on Keck. The complexity of its raw data format neces… ▽ More OSIRIS is a near-infrared (1.0--2.4 $μ$m) integral field spectrograph operating behind the adaptive optics system at Keck Observatory, and is one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics (LGS AO) papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP including characterization of the flux assignment algorithm; exploration of spatial rippling in the reduced data cubes; and improvements to several calibration files, including the rectification matrix, the bad pixel mask, and the wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design. △ Less

Submitted 5 December, 2018; originally announced December 2018.

Comments: 18 pages, 16 figures; accepted for publication in AJ

GD358: three decades of observations for the in-depth asteroseismology of a DBV star

Authors: Agnes Bischoff-Kim, J. L. Provencal, P. A. Bradley, M. H. Montgomery, H. L. Shipman, Samuel T. Harrold, B. Howard, W. Strickland, D. Chandler, D. Campbell, A. Arredondo, R. Linn, D. P. Russell, D. Doyle, A. Brickhouse, D. Peters, S. -L. Kim, X. J. Jiang, Y-N. Mao, A. V. Kusakin, A. V. Sergeev, M. Andreev, S. Velichko, R. Janulis, E. Pakstiene , et al. (16 additional authors not shown)

Abstract: We report on the analysis of 34 years of photometric observations of the pulsating helium atmosphere white dwarf GD358. The complete data set includes archival data from 1982-2006, and 1195.2 hours of new observations from 2007- 2016. From this data set, we extract 15 frequencies representing g-mode pulsation modes, adding 4 modes to the 11 modes known previously. We present evidence that these 15… ▽ More We report on the analysis of 34 years of photometric observations of the pulsating helium atmosphere white dwarf GD358. The complete data set includes archival data from 1982-2006, and 1195.2 hours of new observations from 2007- 2016. From this data set, we extract 15 frequencies representing g-mode pulsation modes, adding 4 modes to the 11 modes known previously. We present evidence that these 15 modes are ell = 1 modes, 13 of which belong to a consecutive sequence in radial overtone k. We perform a detailed asteroseismic analysis using models that include parameterized, complex carbon and oxygen core composition profiles to fit the periods. Recent spectroscopic analyses place GD358 near the red edge of the DBV instability strip, at 24,000 plus or minus 500 K and a log g of 7.8 plus or minus 0.08 dex. The surface gravity translates to a mass range of 0.455 to 0.540 solar masses. Our best fit model has a temperature of 23,650 K and a mass of 0.5706 solar masses. That is slightly more massive than suggested by most the recent spectroscopy. We find a pure helium layer mass of 10^-5.50, consistent with the result of previous studies and the outward diffusion of helium over time. △ Less

Submitted 27 October, 2018; originally announced October 2018.

Comments: 94 pages, 14 figures, 7 tables, Accepted for publication in the ApJ

A Potential Progenitor for the Type Ic Supernova 2017ein

Authors: Charles D. Kilpatrick, Tyler Takaro, Ryan J. Foley, Camille N. Leibler, Yen-Chen Pan, Randall D. Campbell, Wynn V. Jacobson-Galan, Hilton A. Lewis, James E. Lyke, Claire E. Max, Sophia A. Medallon, Armin Rest

Abstract: We report the first detection of a credible progenitor system for a Type Ic supernova (SN Ic), SN 2017ein. We present spectra and photometry of the SN, finding it to be similar to carbon-rich, low-luminosity SNe Ic. Using a post-explosion Keck adaptive optics image, we precisely determine the position of SN 2017ein in pre-explosion \hst\ images, finding a single source coincident with the SN posit… ▽ More We report the first detection of a credible progenitor system for a Type Ic supernova (SN Ic), SN 2017ein. We present spectra and photometry of the SN, finding it to be similar to carbon-rich, low-luminosity SNe Ic. Using a post-explosion Keck adaptive optics image, we precisely determine the position of SN 2017ein in pre-explosion \hst\ images, finding a single source coincident with the SN position. This source is marginally extended, and is consistent with being a stellar cluster. However, under the assumption that the emission of this source is dominated by a single point source, we perform point-spread function photometry, and correcting for line-of-sight reddening, we find it to have $M_{\rm F555W} = -7.5\pm0.2$ mag and $m_{\rm F555W}-m_{\rm F814W}$=$-0.67\pm0.14$ mag. This source is bluer than the main sequence and brighter than almost all Wolf-Rayet stars, however it is similar to some WC+O- and B-star binary systems. Under the assumption that the source is dominated by a single star, we find that it had an initial mass of $55\substack{+20-15} M_{\odot}$. We also examined binary star models to look for systems that match the overall photometry of the pre-explosion source and found that the best-fitting model is a $80$+$48 M_{\odot}$ close binary system in which the $80 M_{\odot}$ star is stripped and explodes as a lower mass star. Late-time photometry after the SN has faded will be necessary to cleanly separate the progenitor star emission from the additional coincident emission. △ Less

Submitted 8 August, 2018; originally announced August 2018.

Comments: 13 pages, 8 figures, published in MNRAS

CCAT-prime: a novel telescope for submillimeter astronomy

Authors: Stephen C. Parshley, Jörg Kronshage, James Blair, Terry Herter, Mike Nolta, Gordon J. Stacey, Andrew Bazarko, Frank Bertoldi, Ricardo Bustos, Donald B. Campbell, Scott Chapman, Nicholas Cothard, Mark Devlin, Jens Erler, Michel Fich, Patricio A. Gallardo, Riccardo Giovanelli, Urs Graf, Scott Gramke, Martha P. Haynes, Richard Hills, Michele Limon, Jeffrey G. Mangum, Jeff McMahon, Michael D. Niemack , et al. (6 additional authors not shown)

Abstract: The CCAT-prime telescope is a 6-meter aperture, crossed-Dragone telescope, designed for millimeter and sub-millimeter wavelength observations. It will be located at an altitude of 5600 meters, just below the summit of Cerro Chajnantor in the high Atacama region of Chile. The telescope's unobscured optics deliver a field of view of almost 8 degrees over a large, flat focal plane, enabling it to acc… ▽ More The CCAT-prime telescope is a 6-meter aperture, crossed-Dragone telescope, designed for millimeter and sub-millimeter wavelength observations. It will be located at an altitude of 5600 meters, just below the summit of Cerro Chajnantor in the high Atacama region of Chile. The telescope's unobscured optics deliver a field of view of almost 8 degrees over a large, flat focal plane, enabling it to accommodate current and future instrumentation fielding >100k diffraction-limited beams for wavelengths less than a millimeter. The mount is a novel design with the aluminum-tiled mirrors nested inside the telescope structure. The elevation housing has an integrated shutter that can enclose the mirrors, protecting them from inclement weather. The telescope is designed to co-host multiple instruments over its nominal 15 year lifetime. It will be operated remotely, requiring minimum maintenance and on-site activities due to the harsh working conditions on the mountain. The design utilizes nickel-iron alloy (Invar) and carbon-fiber-reinforced polymer (CFRP) materials in the mirror support structure, achieving a relatively temperature-insensitive mount. We discuss requirements, specifications, critical design elements, and the expected performance of the CCAT-prime telescope. The telescope is being built by CCAT Observatory, Inc., a corporation formed by an international partnership of universities. More information about CCAT and the CCAT-prime telescope can be found at www.ccatobservatory.org. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: Event: SPIE Astronomical Telescope + Instrumentation, 2018, Austin, Texas, USA; Proceedings Volume 10700, Ground-based and Airborne Telescopes VII; 107005X (2018)

HOD modelling of high redshift galaxies using the BLUETIDES simulation

Authors: Aklant Kumar Bhowmick, Duncan Campbell, Tiziana DiMatteo, Yu Feng

Abstract: We construct halo occupation distribution (HOD) models of high redshift ($z \gtrsim 7.5$) galaxies with $M_{*}>10^8~M_{\odot}/h$ using the BlueTides hydrodynamic simulation suite, with a particular emphasis on modelling the small scale / 1-halo clustering ($0.01\lesssim r \lesssim 1~ h^{-1}\rm{Mpc}$). Similar to low redshift studies, we find that the central and satellite mean HODs (… ▽ More We construct halo occupation distribution (HOD) models of high redshift ($z \gtrsim 7.5$) galaxies with $M_{*}>10^8~M_{\odot}/h$ using the BlueTides hydrodynamic simulation suite, with a particular emphasis on modelling the small scale / 1-halo clustering ($0.01\lesssim r \lesssim 1~ h^{-1}\rm{Mpc}$). Similar to low redshift studies, we find that the central and satellite mean HODs ($\left<N_{\mathrm{cen}}\right>$ and $\left<N_{\mathrm{sat}}\right>$) can be modeled by a smoothed step function and a power law respectively. The number density of satellite galaxies is however significantly suppressed compared to low redshift (satellite fractions drop from $\sim 50 \%$ at $z=0$ to $\lesssim 10 \%$ at $z=7.5$). The mean number of satellites, $\left<N_{\mathrm{sat}}\right> < 1$ for halo masses below $3 \times 10^{11} M_{\odot}/h$ (a rare halo at these redshifts). For the radial number density profiles, satellites with $10^8 \lesssim M^* \lesssim 10^{9} M_{\odot}/h$ in halos with $M_H \gtrsim 3 \times10^{11} M_{\odot}/h$ are consistent with NFW (with concentrations $c_{\mathrm{sat}} \sim 10-40$). Within halos of mass $M_H\lesssim 3 \times 10^{11} M_{\odot}/h$ satellites exhibit a power law profile with slope -3. Because these halos dominate the small scale clustering, the resulting 1-halo term is steeper than predicted using standard NFW profiles. Using this power-law profile for satellites, we can successfully reproduce the small-scale clustering exhibited by BlueTides galaxies using HOD modelling. We predict the highest probability of detecting satellites at $z>7.5$ is around centrals of $M^*\sim 3 \times 10^{10} M_{\odot}/h$ (with $M^{*}\gtrsim$ a few $10^{7} M_{\odot}/h$ ). This should be achievable with the James Webb Space Telescope (JWST). △ Less

Submitted 27 June, 2018; originally announced June 2018.

Comments: 14

DESCQA: An Automated Validation Framework for Synthetic Sky Catalogs

Authors: Yao-Yuan Mao, Eve Kovacs, Katrin Heitmann, Thomas D. Uram, Andrew J. Benson, Duncan Campbell, Sofía A. Cora, Joseph DeRose, Tiziana Di Matteo, Salman Habib, Andrew P. Hearin, J. Bryce Kalmbach, K. Simon Krughoff, François Lanusse, Zarija Lukić, Rachel Mandelbaum, Jeffrey A. Newman, Nelson Padilla, Enrique Paillas, Adrian Pope, Paul M. Ricker, Andrés N. Ruiz, Ananth Tenneti, Cristian Vega-Martínez, Risa H. Wechsler , et al. (2 additional authors not shown)

Abstract: The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline developm… ▽ More The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs. △ Less

Submitted 8 February, 2018; v1 submitted 27 September, 2017; originally announced September 2017.

Comments: 16+4 pages, 9 figures, 2+1 tables. Published in the ApJ Supplement. Code repository at https://github.com/LSSTDESC/descqa

Journal ref: ApJS (2018), 234, 36

The Galaxy Clustering Crisis in Abundance Matching

Authors: Duncan Campbell, Frank C. van den Bosch, Nikhil Padmanabhan, Yao-Yuan Mao, Andrew R. Zentner, Johannes U. Lange, Fangzhou Jiang, Antonia Villarreal

Abstract: Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxi… ▽ More Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxies produced in these models compared to those based on $V_{\rm peak}$. Despite success in predicting clustering, a simple $V_{\rm peak}$ based SHAM model results in predictions for galaxy growth that are at odds with observations. We evaluate three possible remedies that could "save" mass-based SHAM: (1) SHAM models require a significant population of "orphan" galaxies as a result of artificial disruption/merging of sub-haloes in modern high resolution dark matter simulations; (2) satellites must grow significantly after their accretion; and (3) stellar mass is significantly affected by halo assembly history. No solution is entirely satisfactory. However, regardless of the particulars, we show that popular SHAM models based on $M_{\rm peak}$ cannot be complete physical models as presented. Either $V_{\rm peak}$ truly is a better predictor of stellar mass at $z\sim 0$ and it remains to be seen how the correlation between stellar mass and $V_{\rm peak}$ comes about, or SHAM models are missing vital component(s) that significantly affect galaxy clustering. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: 25 pages, 22 figures, submitted to MNRAS, comments welcome

Brightest galaxies as halo centre tracers in SDSS DR7

Authors: Johannes U. Lange, Frank C. van den Bosch, Andrew Hearin, Duncan Campbell, Andrew R. Zentner, Antonia Villarreal, Yao-Yuan Mao

Abstract: Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight veloc… ▽ More Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight velocity and spatial offsets between brightest galaxies and their neighbours. We compare those to detailed mock catalogues, constructed from high-resolution, dark-matter-only $N$-body simulations, in which it is assumed that satellite galaxies trace dark matter subhaloes. This allows us to place constraints on the fraction $f_{\rm BNC}$ of haloes in which the brightest galaxy is not the central. Compared to previous studies we explicitly take into account the unrelaxed state of the host haloes, velocity offsets of halo cores and correlations between $f_{\rm BNC}$ and the satellite occupation. We find that $f_{\rm BNC}$ strongly decreases with the luminosity of the brightest galaxy and increases with the mass of the host halo. Overall, in the halo mass range $10^{13} - 10^{14.5} h^{-1} M_\odot$ we find $f_{\rm BNC} \sim 30\%$, in good agreement with a previous study by Skibba et al. We discuss the implications of these findings for studies inferring the galaxy--halo connection from satellite kinematics, models of the conditional luminosity function and galaxy formation in general. △ Less

Submitted 28 November, 2017; v1 submitted 14 May, 2017; originally announced May 2017.

Comments: 24 pages, 15 figures. Accepted for publication in MNRAS

The Immitigable Nature of Assembly Bias: The Impact of Halo Definition on Assembly Bias

Authors: Antonia S. Villarreal, Andrew R. Zentner, Yao-Yuan Mao, Chris W. Purcell, Frank C. van den Bosch, Benedikt Diemer, Johannes U. Lange, Kuan Wang, Duncan Campbell

Abstract: Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $Δ$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo m… ▽ More Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $Δ$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $Δ$. For conventional halo definitions ($Δ\sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $Δ\sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $Δ$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $Δ\gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass. △ Less

Submitted 16 September, 2017; v1 submitted 11 May, 2017; originally announced May 2017.

Comments: 19 pages, 13 figures. Updated to published version. Main result summarized in Figure 10

Journal ref: MNRAS, 472, 1088-1105 (2017)

Quenching and ram pressure stripping of simulated Milky Way satellite galaxies

Authors: Christine M. Simpson, Robert J. J. Grand, Facundo A. Gómez, Federico Marinacci, Rüdiger Pakmor, Volker Springel, David J. R. Campbell, Carlos S. Frenk

Abstract: We present predictions for the quenching of star formation in satellite galaxies of the Local Group from a suite of 30 cosmological zoom simulations of Milky Way-like host galaxies. The Auriga simulations resolve satellites down to the luminosity of the classical dwarf spheroidal galaxies of the Milky Way. We find strong mass-dependent and distance-dependent quenching signals, where dwarf systems… ▽ More We present predictions for the quenching of star formation in satellite galaxies of the Local Group from a suite of 30 cosmological zoom simulations of Milky Way-like host galaxies. The Auriga simulations resolve satellites down to the luminosity of the classical dwarf spheroidal galaxies of the Milky Way. We find strong mass-dependent and distance-dependent quenching signals, where dwarf systems beyond 600 kpc are only strongly quenched below a stellar mass of $10^7$ M$_\odot$. Ram pressure stripping appears to be the dominant quenching mechanism and 50% of quenched systems cease star formation within 1 Gyr of first infall. We demonstrate that systems within a host galaxy's $R_{200}$ radius are comprised of two populations: (i) a first infall population that has entered the host halo within the past few Gyrs and (ii) a population of returning `backsplash' systems that have had a much more extended interaction with the host. Backsplash galaxies that do not return to the host galaxy by redshift zero exhibit quenching properties similar to galaxies within $R_{200}$ and are distinct from other external systems. The simulated quenching trend with stellar mass has some tension with observations, but our simulations are able reproduce the range of quenching times measured from resolved stellar populations of Local Group dwarf galaxies. △ Less

Submitted 20 March, 2018; v1 submitted 8 May, 2017; originally announced May 2017.

Comments: 22 pages, 17 figures, accepted for publication in MNRAS

Magnetic field formation in the Milky Way-like disk galaxies of the Auriga project

Authors: Ruediger Pakmor, Facundo A. Gomez, Robert J. J. Grand, Federico Marinacci, Christine M. Simpson, Volker Springel, David J. R. Campbell, Carlos S. Frenk, Thomas Guillet, Christoph Pfrommer, Simon D. M. White

Abstract: The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogethe… ▽ More The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{μG}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution. △ Less

Submitted 2 May, 2017; v1 submitted 24 January, 2017; originally announced January 2017.

Comments: 16 pages, 16 figures, accepted by MNRAS, section 4.1 (turbulent dynamo) extended substantially, now includes magnetic and kinetic power spectra

Thermal Properties and an Improved Shape Model for Near-Earth Asteroid (162421) 2000 ET70

Authors: Sean E. Marshall, Ellen S. Howell, Christopher Magri, Ronald J. Vervack Jr., Donald B. Campbell, Yanga R. Fernández, Michael C. Nolan, Jenna L. Crowell, Michael D. Hicks, Kenneth J. Lawrence, Patrick A. Taylor

Abstract: We present thermal properties and an improved shape model for potentially hazardous asteroid (162421) 2000 ET70. In addition to the radar data from 2000 ET70's apparition in 2012, our model incorporates optical lightcurves and infrared spectra that were not included in the analysis of Naidu et al. (2013, Icarus 226, 323-335). We confirm the general "clenched fist" appearance of the Naidu et al. mo… ▽ More We present thermal properties and an improved shape model for potentially hazardous asteroid (162421) 2000 ET70. In addition to the radar data from 2000 ET70's apparition in 2012, our model incorporates optical lightcurves and infrared spectra that were not included in the analysis of Naidu et al. (2013, Icarus 226, 323-335). We confirm the general "clenched fist" appearance of the Naidu et al. model, but compared to their model, our best-fit model is about 10% longer along its long principal axis, nearly identical along the intermediate axis, and about 25% shorter along the short axis. We find the asteroid's dimensions to be 2.9 km $\times$ 2.2 km $\times$ 1.5 km (with relative uncertainties of about 10%, 15%, and 25%, respectively). With the available data, 2000 ET70's period and pole position are degenerate with each other. The radar and lightcurve data together constrain the pole direction to fall along an arc that is about twenty-three degrees long and eight degrees wide. Infrared spectra from the NASA InfraRed Telescope Facility (IRTF) provide an additional constraint on the pole. Thermophysical modeling, using our SHERMAN software, shows that only a subset of the pole directions, about twelve degrees of that arc, are compatible with the infrared data. Using all of the available data, we find that 2000 ET70 has a sidereal rotation period of 8.944 hours ($\pm$ 0.009 h) and a north pole direction of ecliptic coordinates $(52^{\circ}, -60^{\circ}) \pm 6^{\circ}$. The infrared data, acquired over several dates, require that the thermal properties (albedo, thermal inertia, surface roughness) must change across the asteroid's surface. By incorporating the detailed shape model and spin state into our thermal modeling, the multiple ground-based observations at different viewing geometries have allowed us to constrain the levels of the variations in the surface properties. △ Less

Submitted 10 March, 2017; v1 submitted 14 October, 2016; originally announced October 2016.

Comments: Revised and submitted to Icarus; 22 pages, 13 figures; 51 supplementary figures (in a 25MB PDF) can be found at http://astro.cornell.edu/~seanm/2000et70/

The Auriga Project: the properties and formation mechanisms of disc galaxies across cosmic time

Authors: Robert J. J. Grand, Facundo A. Gómez, Federico Marinacci, Ruediger Pakmor, Volker Springel, David J. R. Campbell, Carlos S. Frenk, Adrian Jenkins, Simon D. M. White

Abstract: We introduce a suite of thirty cosmological magneto-hydrodynamical zoom simulations of the formation of galaxies in isolated Milky Way mass dark haloes. These were carried out with the moving mesh code \textlcsc{AREPO}, together with a comprehensive model for galaxy formation physics, including AGN feedback and magnetic fields, which produces realistic galaxy populations in large cosmological simu… ▽ More We introduce a suite of thirty cosmological magneto-hydrodynamical zoom simulations of the formation of galaxies in isolated Milky Way mass dark haloes. These were carried out with the moving mesh code \textlcsc{AREPO}, together with a comprehensive model for galaxy formation physics, including AGN feedback and magnetic fields, which produces realistic galaxy populations in large cosmological simulations. We demonstrate that our simulations reproduce a wide range of present-day observables, in particular, two component disc dominated galaxies with appropriate stellar masses, sizes, rotation curves, star formation rates and metallicities. We investigate the driving mechanisms that set present-day disc sizes/scale lengths, and find that they are related to the angular momentum of halo material. We show that the largest discs are produced by quiescent mergers that inspiral into the galaxy and deposit high angular momentum material into the pre-existing disc, simultaneously increasing the spin of dark matter and gas in the halo. More violent mergers and strong AGN feedback play roles in limiting disc size by destroying pre-existing discs and by suppressing gas accretion onto the outer disc, respectively. The most important factor that leads to compact discs, however, is simply a low angular momentum for the halo. In these cases, AGN feedback plays an important role in limiting central star formation and the formation of a massive bulge. △ Less

Submitted 11 January, 2017; v1 submitted 4 October, 2016; originally announced October 2016.

Comments: 32 pages, 25 figures. MNRAS accepted after significant improvements. Information, images and movies of the Auriga project can be found at http://auriga.h-its.org

Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection

Authors: ChangHoon Hahn, Mohammadjavad Vakili, Kilian Walsh, Andrew P. Hearin, David W. Hogg, Duncan Campbell

Abstract: Standard approaches to Bayesian parameter inference in large scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be correct. Likelihood free inferences such as Approximate Bayesian Computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a fo… ▽ More Standard approaches to Bayesian parameter inference in large scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be correct. Likelihood free inferences such as Approximate Bayesian Computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a forward generative model of the data and a metric for measuring the distance between the model and data. In this work, we demonstrate that ABC is feasible for LSS parameter inference by using it to constrain parameters of the halo occupation distribution (HOD) model for populating dark matter halos with galaxies. Using specific implementation of ABC supplemented with Population Monte Carlo importance sampling, a generative forward model using HOD, and a distance metric based on galaxy number density, two-point correlation function, and galaxy group multiplicity function, we constrain the HOD parameters of mock observation generated from selected "true" HOD parameters. The parameter constraints we obtain from ABC are consistent with the "true" HOD parameters, demonstrating that ABC can be reliably used for parameter inference in LSS. Furthermore, we compare our ABC constraints to constraints we obtain using a pseudo-likelihood function of Gaussian form with MCMC and find consistent HOD parameter constraints. Ultimately our results suggest that ABC can and should be applied in parameter inference for LSS analyses. △ Less

Submitted 10 April, 2017; v1 submitted 6 July, 2016; originally announced July 2016.

Comments: 16 pages, 10 figures

Forward Modeling of Large-Scale Structure: An open-source approach with Halotools

Authors: Andrew Hearin, Duncan Campbell, Erik Tollerud, Peter Behroozi, Benedikt Diemer, Nathan J. Goldbaum, Elise Jennings, Alexie Leauthaud, Yao-Yuan Mao, Surhud More, John Parejko, Manodeep Sinha, Brigitta Sipocz, Andrew Zentner

Abstract: We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo model… ▽ More We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python. △ Less

Submitted 22 September, 2017; v1 submitted 13 June, 2016; originally announced June 2016.

Comments: Revisions match version accepted for publication in AAS

Spiral-induced velocity and metallicity patterns in a cosmological zoom simulation of a Milky Way-sized galaxy

Authors: Robert J. J. Grand, Volker Springel, Daisuke Kawata, Ivan Minchev, Patricia Sánchez-Blázquez, Facundo A. Gómez, Federico Marinacci, Rüdiger Pakmor, David J. R. Campbell

Abstract: We use a high resolution cosmological zoom simulation of a Milky Way-sized halo to study the observable features in velocity and metallicity space associated with the dynamical influence of spiral arms. For the first time, we demonstrate that spiral arms, that form in a disc in a fully cosmological environment with realistic galaxy formation physics, drive large-scale systematic streaming motions.… ▽ More We use a high resolution cosmological zoom simulation of a Milky Way-sized halo to study the observable features in velocity and metallicity space associated with the dynamical influence of spiral arms. For the first time, we demonstrate that spiral arms, that form in a disc in a fully cosmological environment with realistic galaxy formation physics, drive large-scale systematic streaming motions. In particular, on the trailing edge of the spiral arms the peculiar galacto-centric radial and azimuthal velocity field is directed radially outward and azimuthally backward, whereas it is radially inward and azimuthally forward on the leading edge. Owing to the negative radial metallicity gradient, this systematic motion drives, at a given radius, an azimuthal variation in the residual metallicity that is characterised by a metal rich trailing edge and a metal poor leading edge. We show that these signatures are theoretically observable in external galaxies with Integral Field Unit instruments such as VLT/MUSE, and if detected, would provide evidence for large-scale systematic radial migration driven by spiral arms. △ Less

Submitted 16 May, 2016; v1 submitted 4 April, 2016; originally announced April 2016.

Comments: 5 pages, 3 figures. Accepted to MNRAS after minor corrections

Knowing the unknowns: uncertainties in simple estimators of galactic dynamical masses

Authors: David J. R. Campbell, Carlos S. Frenk, Adrian Jenkins, Vincent R. Eke, Julio F. Navarro, Till Sawala, Matthieu Schaller, Azadeh Fattahi, Kyle A. Oman, Tom Theuns

Abstract: The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make us… ▽ More The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make use of the APOSTLE cosmological hydrodynamical simulations of the Local Group to test the validity and accuracy of such mass estimators for both dispersion and rotation-supported galaxies, for field and satellite galaxies, and for galaxies of varying masses, shapes, and velocity dispersion anisotropies. We find that the mass estimators of Walker et al. and Wolf et al. are able to recover the masses of dispersion-dominated systems with little systematic bias, but with a 1-sigma scatter of 25 and 23 percent, respectively. The error on the estimated mass is dominated by the impact of the 3D shape of the stellar mass distribution, which is difficult to constrain observationally. This intrinsic scatter becomes the dominant source of uncertainty in the masses estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky Way, where the observational errors in their sizes and velocity dispersions are small. Such scatter may also affect the inner density slopes of dSphs derived from multiple stellar populations, relaxing the significance with which Navarro-Frenk-White profiles may be excluded, depending on the degree to which the relevant properties of the different stellar populations are correlated. Finally, we derive a new optimal mass estimator that removes the residual biases and achieves a statistically significant reduction in the scatter to 20 percent overall for dispersion-dominated galaxies, allowing more precise and accurate mass estimates. △ Less

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

Comments: Accepted for publication in MNRAS. 26 pages, 21 figures, 3 tables

Introducing Decorated HODs: modeling assembly bias in the galaxy-halo connection

Authors: Andrew P. Hearin, Andrew R. Zentner, Frank C. van den Bosch, Duncan Campbell, Erik Tollerud

Abstract: The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models… ▽ More The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models designed to account for assembly bias. Decorated HODs minimally expand the parameter space and maximize the independence between traditional and novel HOD parameters. We use decorated HODs to quantify the influence of assembly bias on clustering and lensing statistics. For SDSS-like samples, the impact of assembly bias on galaxy clustering can be as large as a factor of two on r ~ 200 kpc scales and ~15% in the linear regime. Assembly bias can either enhance or diminish clustering on large scales, but generally increases clustering on scales r <~ 1 Mpc. We performed our calculations with Halotools, an open-source, community-driven python package for studying the galaxy-halo connection (http://halotools.readthedocs.org). We conclude by describing the use of decorated HODs to treat assembly bias in otherwise conventional likelihood analyses. △ Less

Submitted 9 December, 2015; originally announced December 2015.

Comments: Figure 1 provides a cartoon visualization of the analytical formalism. Figures 4 & 5 show the characteristic signature that assembly bias imprints on the clustering and lensing of SDSS-like galaxy samples. 21 pages, 8 figures, submitted to MNRAS

Vertical disc heating in Milky Way-sized galaxies in a cosmological context

Authors: Robert J. J. Grand, Volker Springel, Facundo A. Gómez, Federico Marinacci, Rüdiger Pakmor, David J. R. Campbell, Adrian Jenkins

Abstract: Vertically extended, high velocity dispersion stellar distributions appear to be a ubiquitous feature of disc galaxies, and both internal and external mechanisms have been proposed to be the major driver of their formation. However, it is unclear to what extent each mechanism can generate such a distribution, which is likely to depend on the assembly history of the galaxy. To this end, we perform… ▽ More Vertically extended, high velocity dispersion stellar distributions appear to be a ubiquitous feature of disc galaxies, and both internal and external mechanisms have been proposed to be the major driver of their formation. However, it is unclear to what extent each mechanism can generate such a distribution, which is likely to depend on the assembly history of the galaxy. To this end, we perform 16 high resolution cosmological-zoom simulations of Milky Way-sized galaxies using the state-of-the-art cosmological magneto-hydrodynamical code \textlcsc{AREPO}, and analyse the evolution of the vertical kinematics of the stellar disc in connection with various heating mechanisms. We find that the bar is the dominant heating mechanism in most cases, whereas spiral arms, radial migration, and adiabatic heating from mid-plane density growth are all sub-dominant. The strongest source, though less prevalent than bars, originates from external perturbations from satellites/sub-halos of masses log$_{10} (M/\rm M_{\odot}) \gtrsim 10$. However, in many simulations the orbits of newborn star particles become cooler with time, such that they dominate the shape of the age-velocity dispersion relation and overall vertical disc structure unless a strong external perturbation takes place. △ Less

Submitted 2 April, 2016; v1 submitted 7 December, 2015; originally announced December 2015.

Comments: 24 pages, 15 figures (new figs. 6&7). Accepted in MNRAS. Much improved over originally submitted version

On the Segregation of Dark Matter Substructure

Authors: Frank C. van den Bosch, Fangzhou Jiang, Duncan Campbell, Peter Behroozi

Abstract: We present the first comprehensive analysis of the segregation of dark matter subhaloes in their host haloes. Using numerical simulations, we examine the segregation of twelve different subhalo properties with respect to both orbital energy and halo-centric radius (in real space as well as in projection). Subhaloes are strongly segregated by accretion redshift, which is an outcome of the inside-ou… ▽ More We present the first comprehensive analysis of the segregation of dark matter subhaloes in their host haloes. Using numerical simulations, we examine the segregation of twelve different subhalo properties with respect to both orbital energy and halo-centric radius (in real space as well as in projection). Subhaloes are strongly segregated by accretion redshift, which is an outcome of the inside-out assembly of their host haloes. Since subhaloes that were accreted earlier have experienced more tidal stripping, subhaloes that have lost a larger fraction of their mass at infall are on more bound orbits. Subhaloes are also strongly segregated in their masses and maximum circular velocities at accretion. We demonstrate that part of this segregation is already imprinted in the infall conditions. For massive subhaloes it is subsequently boosted by dynamical friction, but only during their first radial orbit. The impact of these two effects is counterbalanced, though, by the fact that subhaloes with larger accretion masses are accreted later. Because of tidal stripping, subhaloes reveal little to no segregation by present-day mass or maximum circular velocity, while the corresponding torques cause subhaloes on more bound orbits to have smaller spin. There is a weak tendency for subhaloes that formed earlier to be segregated towards the center of their host halo, which is an indirect consequence of the fact that (sub)halo formation time is correlated with other, strongly segregated properties. We discuss the implications of our results for the segregation of satellite galaxies in galaxy groups and clusters. △ Less

Submitted 6 October, 2015; originally announced October 2015.

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

Assessing Colour-dependent Occupation Statistics Inferred from Galaxy Group Catalogues

Authors: Duncan Campbell, Frank C van den Bosch, Andrew Hearin, Nikhil Padmanabhan, Andreas Berlind, H. J. Mo, Jeremy Tinker, Xiaohu Yang

Abstract: We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most col… ▽ More We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent. △ Less

Submitted 18 May, 2015; originally announced May 2015.

Comments: 27 pages, 18 figures, accepted for publication in MNRAS

A new methodology to test galaxy formation models using the dependence of clustering on stellar mass

Authors: David J. R. Campbell, Carlton M. Baugh, Peter D. Mitchell, John C. Helly, Violeta Gonzalez-Perez, Cedric G. Lacey, Claudia del P. Lagos, Vimal Simha, Daniel J. Farrow

Abstract: We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different… ▽ More We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger timescales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from GAMA, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from SDSS and VIPERS, depending on the GALFORM model used. △ Less

Submitted 10 June, 2015; v1 submitted 11 December, 2014; originally announced December 2014.

Comments: Accepted for publication in MNRAS; 20 pages, 9 figure, 2 tables

The first planet detected in the WTS: an inflated hot-Jupiter in a 3.35 day orbit around a late F-star [ERRATUM]

Authors: M. Cappetta, R. P. Saglia, J. L. Birkby, J. Koppenhoefer, D. J. Pinfield, S. T. Hodgkin, P. Cruz, G. Kovacs, B. Sipocz, D. Barrado, B. Nefs, Y. V. Pavlenko, L. Fossati, C. del Burgo, E. L. Martin, I. Snellen, J. Barnes, D. A. Campbell, S. Catalan, M. C. Galvez-Ortiz, N. Goulding, C. Haswell, O. Ivanyuk, H. Jones, M. Kuznetsov , et al. (13 additional authors not shown)

Abstract: We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope (UKIRT) in August 2007. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J ~ 16 were constructed for ~60000 stars and searched for periodic transit signals. For one of the most p… ▽ More We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope (UKIRT) in August 2007. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J ~ 16 were constructed for ~60000 stars and searched for periodic transit signals. For one of the most promising transiting candidates, high-resolution spectra taken at the Hobby-Eberly Telescope (HET) allowed us to estimate the spectroscopic parameters of the host star, a late-F main sequence dwarf (V=16.13) with possibly slightly subsolar metallicity, and to measure its radial velocity variations. The combined analysis of the light curves and spectroscopic data resulted in an orbital period of the substellar companion of 3.35 days, a planetary mass of 4.01 +- 0.35 Mj and a planetary radius of 1.49+0.16-0.18 Rj. WTS-1b has one of the largest radius anomalies among the known hot Jupiters in the mass range 3-5 Mj. The high irradiation from the host star ranks the planet in the pM class. △ Less

Submitted 13 October, 2014; originally announced October 2014.

Comments: 16 pages, 10 figures

Coming of Age in the Dark Sector: How Dark Matter Haloes grow their Gravitational Potential Wells

Authors: Frank C. van den Bosch, Fangzhou Jiang, Andrew Hearin, Duncan Campbell, Douglas Watson, Nikhil Padmanabhan

Abstract: We present a detailed study of how dark matter haloes assemble their mass and grow their (central) potential well. We characterize these via their mass accretion histories (MAHs) and potential well growth histories (PWGHs), which we extract from the Bolshoi simulation and from semi-analytical merger trees supplemented with a method to compute the maximum circular velocity, Vmax, of progenitor halo… ▽ More We present a detailed study of how dark matter haloes assemble their mass and grow their (central) potential well. We characterize these via their mass accretion histories (MAHs) and potential well growth histories (PWGHs), which we extract from the Bolshoi simulation and from semi-analytical merger trees supplemented with a method to compute the maximum circular velocity, Vmax, of progenitor haloes. The results of both methods are in excellent agreement, both in terms of the average and the scatter. We show that the MAH and PWGH are tightly correlated, and that growth of the central potential precedes the assembly of mass; the maximum circular velocity is already half the present day value by the time the halo has accreted only 2 percent of its final mass. Finally, we demonstrate that MAHs have a universal form, which we use to develop a new and improved universal model that can be used to compute the average or median MAH and PWGH for a halo of any mass in any LCDM cosmology, without having to run a numerical simulation or a set of halo merger trees. △ Less

Submitted 9 September, 2014; originally announced September 2014.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

Keck-I MOSFIRE spectroscopy of compact star-forming galaxies at z$\gtrsim$2: High velocity dispersions in progenitors of compact quiescent galaxies

Authors: G. Barro, J. R. Trump, D. C. Koo, A. Dekel, S. A. Kassin, D. D. Kocevski, S. M. Faber, A. van der Wel, Y. Guo, P. G. Perez-Gonzalez, E. Toloba, J. J. Fang, C. Pacifici, R. Simons, R D. Campbell, D. Ceverino, S. L. Finkelstein, B. Goodrich, M. Kassis, A. M. Koekemoer, N. P. Konidaris, R. C. Livermore, J. E. Lyke, B. Mobasher, H. Nayyeri , et al. (6 additional authors not shown)

Abstract: We present Keck-I MOSFIRE near-infrared spectroscopy for a sample of 13 compact star-forming galaxies (SFGs) at redshift $2\leq z \leq2.5$ with star formation rates of SFR$\sim$100M$_{\odot}$ y$^{-1}$ and masses of log(M/M$_{\odot}$)$\sim10.8$. Their high integrated gas velocity dispersions of $σ_{\rm{int}}$=230$^{+40}_{-30}$ km s$^{-1}$, as measured from emission lines of H$_α$ and [OIII], and th… ▽ More We present Keck-I MOSFIRE near-infrared spectroscopy for a sample of 13 compact star-forming galaxies (SFGs) at redshift $2\leq z \leq2.5$ with star formation rates of SFR$\sim$100M$_{\odot}$ y$^{-1}$ and masses of log(M/M$_{\odot}$)$\sim10.8$. Their high integrated gas velocity dispersions of $σ_{\rm{int}}$=230$^{+40}_{-30}$ km s$^{-1}$, as measured from emission lines of H$_α$ and [OIII], and the resultant M$_{\star}-σ_{\rm{int}}$ relation and M$_{\star}$$-$M$_{\rm{dyn}}$ all match well to those of compact quiescent galaxies at $z\sim2$, as measured from stellar absorption lines. Since log(M$_{\star}$/M$_{\rm{dyn}}$)$=-0.06\pm0.2$ dex, these compact SFGs appear to be dynamically relaxed and more evolved, i.e., more depleted in gas and dark matter ($<$13$^{+17}_{-13}$\%) than their non-compact SFG counterparts at the same epoch. Without infusion of external gas, depletion timescales are short, less than $\sim$300 Myr. This discovery adds another link to our new dynamical chain of evidence that compact SFGs at $z\gtrsim2$ are already losing gas to become the immediate progenitors of compact quiescent galaxies by $z\sim2$. △ Less

Submitted 27 May, 2014; originally announced May 2014.

Comments: 12 pages, 7 figures, submitted to ApJ

WTS-2 b: a hot Jupiter orbiting near its tidal destruction radius around a K-dwarf

Authors: J. L. Birkby, M. Cappetta, P. Cruz, J. Koppenhoefer, O. Ivanyuk, A. J. Mustill, S. T. Hodgkin, D. J. Pinfield, B. Sipőcz, G. Kovács, R. Saglia, Y. Pavlenko, D. Barrado, A. Bayo, D. Campbell, S. Catalan, L. Fossati, M. -C. Gálvez-Ortiz, M. Kenworthy, J. Lillo-Box E. L. Martín, D. Mislis, E. J. W. de Mooij, S. V. Nefs, I. A. G. Snellen, H. Stoev , et al. (12 additional authors not shown)

Abstract: We report the discovery of WTS-2 b, an unusually close-in 1.02-day hot Jupiter (Mp=1.12MJ, Rp=1.363RJ) orbiting a K2V star, which has a possible gravitationally-bound M-dwarf companion at 0.6 arcsec separation contributing ~20 percent of the total flux in the observed J-band light curve. The planet is only 1.5 times the separation from its host star at which it would be destroyed by Roche lobe ove… ▽ More We report the discovery of WTS-2 b, an unusually close-in 1.02-day hot Jupiter (Mp=1.12MJ, Rp=1.363RJ) orbiting a K2V star, which has a possible gravitationally-bound M-dwarf companion at 0.6 arcsec separation contributing ~20 percent of the total flux in the observed J-band light curve. The planet is only 1.5 times the separation from its host star at which it would be destroyed by Roche lobe overflow, and has a predicted remaining lifetime of just ~40 Myr, assuming a tidal dissipation quality factor of Q'*=10^6. Q'* is a key factor in determining how frictional processes within a host star affect the orbital evolution of its companion giant planets, but it is currently poorly constrained by observations. We calculate that the orbital decay of WTS-2 b would correspond to a shift in its transit arrival time of T_shift~17 seconds after 15 years assuming Q'*=10^6. A shift less than this would place a direct observational constraint on the lower limit of Q'* in this system. We also report a correction to the previously published expected T_shift for WASP-18 b, finding that T_shift=356 seconds after 10 years for Q'*=10^6, which is much larger than the estimated 28 seconds quoted in WASP-18 b discovery paper. We attempted to constrain Q'* via a study of the entire population of known transiting hot Jupiters, but our results were inconclusive, requiring a more detailed treatment of transit survey sensitivities at long periods. We conclude that the most informative and straight-forward constraints on Q'* will be obtained by direct observational measurements of the shift in transit arrival times in individual hot Jupiter systems. We show that this is achievable across the mass spectrum of exoplanet host stars within a decade, and will directly probe the effects of stellar interior structure on tidal dissipation. △ Less

Submitted 21 February, 2014; originally announced February 2014.

Comments: 21 pages, 16 figures, 6 tables, accepted for publication in MNRAS

Efficiency Measurements and Installation of a New Grating for the OSIRIS Spectrograph at Keck Observatory

Authors: Etsuko Mieda, Shelley A. Wright, James E. Larkin, James R. Graham, Sean M. Adkins, James E. Lyke, Randy D. Campbell, Jerome Maire, Tuan Do, Jacob Gordon

Abstract: OSIRIS is a near-infrared integral field spectrograph operating behind the adaptive optics system at W. M. Keck Observatory. While OSIRIS has been a scientifically productive instrument to date, its sensitivity has been limited by a grating efficiency that is less than half of what was expected. The spatially averaged efficiency of the old grating, weighted by error, is measured to be 39.5 +/- 0.8… ▽ More OSIRIS is a near-infrared integral field spectrograph operating behind the adaptive optics system at W. M. Keck Observatory. While OSIRIS has been a scientifically productive instrument to date, its sensitivity has been limited by a grating efficiency that is less than half of what was expected. The spatially averaged efficiency of the old grating, weighted by error, is measured to be 39.5 +/- 0.8 % at λ = 1.310 μm, with large field dependent variation of 11.7 % due to efficiency variation across the grating surface. Working with a new vendor, we developed a more efficient and uniform grating with a weighted average efficiency at λ = 1.310 μm of 78.0 +/- 1.6 %, with field variation of only 2.2 %. This is close to double the average efficiency and five times less variation across the field. The new grating was installed in December 2012, and on- sky OSIRIS throughput shows an average factor of 1.83 improvement in sensitivity between 1 and 2.4 microns. We present the development history, testing, and implementation of this new near-infrared grating for OSIRIS and report the comparison with the predecessors. The higher sensitivities are already having a large impact on scientific studies with OSIRIS. △ Less

Submitted 17 January, 2014; originally announced January 2014.

Journal ref: PASP. 126 (2014) 250-263

A new model of galaxy formation: How sensitive are predicted galaxy luminosities to the choice of SPS model?

Authors: V. Gonzalez-Perez, C. G. Lacey, C. M. Baugh, C. D. P. Lagos, J. Helly, D. J. R. Campbell, P. D. Mitchell

Abstract: We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the WMAP7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven diff… ▽ More We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the WMAP7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven different SPS models. In each case we obtain the rest-frame luminosity function in the far-ultra-violet, optical and near-infrared (NIR) wavelength ranges. We find that both the predicted rest-frame ultra-violet and optical luminosity function are insensitive to the choice of SPS model. However, we find that the predicted evolution of the rest-frame NIR luminosity function depends strongly on the treatment of the thermally pulsating asymptotic giant branch (TP-AGB) stellar phase in the SPS models, with differences larger than a factor of 2 for model galaxies brighter than $M_{\rm AB}(K)-5$log$h<-22$ ($\sim$L$_*$ for $0\leq z\leq 1.5$). We have also explored the predicted number counts of galaxies, finding remarkable agreement between the results with different choices of SPS model, except when selecting galaxies with very red optical-NIR colours. The predicted number counts of these extremely red galaxies appear to be more affected by the treatment of star formation in disks than by the treatment of TP-AGB stars in the SPS models. △ Less

Submitted 6 February, 2014; v1 submitted 26 September, 2013; originally announced September 2013.

Comments: 22 pages, 19 figures, MNRAS. Monthly Notices of the Royal Astronomical Society 2014

Phase-Resolved Infrared Spectroscopy and Photometry of V1500 Cygni, and a Search for Similar Old Classical Novae

Authors: Thomas E. Harrison, Randy D. Campbell, James E. Lyke

Abstract: We present phase-resolved near-infrared photometry and spectroscopy of the classical nova V1500 Cyg to explore whether cyclotron emission is present in this system. While the spectroscopy do not indicate the presence of discrete cyclotron harmonic emission, the light curves suggest that a sizable fraction of its near-infrared fluxes are due to this component. The light curves of V1500 Cyg appear t… ▽ More We present phase-resolved near-infrared photometry and spectroscopy of the classical nova V1500 Cyg to explore whether cyclotron emission is present in this system. While the spectroscopy do not indicate the presence of discrete cyclotron harmonic emission, the light curves suggest that a sizable fraction of its near-infrared fluxes are due to this component. The light curves of V1500 Cyg appear to remain dominated by emission from the heated face of the secondary star in this system. We have used infrared spectroscopy and photometry to search for other potential magnetic systems amongst old classical novae. We have found that the infrared light curves of V1974 Cyg superficially resemble those of V1500 Cyg, suggesting a highly irradiated companion. The old novae V446 Her and QV Vul have light curves with large amplitude variations like those seen in polars, suggesting they might have magnetic primaries. We extract photometry for seventy nine old novae from the 2MASS Point Source Catalog and use those data to derive the mean, un-reddened infrared colors of quiescent novae. We also extract WISE data for these objects and find that forty five of them were detected. Surprisingly, a number of these systems were detected in the WISE 22 mum band. While two of those objects produced significant dust shells (V705 Cas and V445 Pup), the others did not. It appears that line emission from their ionized ejected shells is the most likely explanation for those detections. △ Less

Submitted 6 June, 2013; originally announced June 2013.

Comments: Accepted in the Astronomical Journal

The first planet detected in the WTS: an inflated hot-Jupiter in a 3.35 day orbit around a late F-star

Authors: M. Cappetta, R. P. Saglia, J. L. Birkby, J. Koppenhoefer, D. J. Pinfield, S. T. Hodgkin, P. Cruz, G. Kovács, B. Sipöcz, D. Barrado, B. Nefs, Y. V. Pavlenko, L. Fossati, C. del Burgo, E. L. Martín, I. Snellen, J. Barnes, A. M. Bayo, D. A. Campbell, S. Catalan, M. C. Gálvez-Ortiz, N. Goulding, C. Haswell, O. Ivanyuk, H. Jones , et al. (14 additional authors not shown)

Abstract: We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J=16 were constructed for 60000 stars and searched for periodic transit signals. For one of the most promising transiting candid… ▽ More We report the discovery of WTS-1b, the first extrasolar planet found by the WFCAM Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Telescope. Light curves comprising almost 1200 epochs with a photometric precision of better than 1 per cent to J=16 were constructed for 60000 stars and searched for periodic transit signals. For one of the most promising transiting candidates, high-resolution spectra taken at the Hobby-Eberly Telescope allowed us to estimate the spectroscopic parameters of the host star, a late-F main sequence dwarf (V=16.13) with possibly slightly subsolar metallicity, and to measure its radial velocity variations. The combined analysis of the light curves and spectroscopic data resulted in an orbital period of the substellar companion of 3.35 days, a planetary mass of 4.01+-0.35 Mj and a planetary radius of 1.49+-0.17 Rj. WTS-1b has one of the largest radius anomalies among the known hot Jupiters in the mass range 3-5 Mj. △ Less

Submitted 3 October, 2012; originally announced October 2012.

Comments: 15 pages, 10 figures, 9 tables, accepted for publication in MNRAS