Showing 1–50 of 127 results for author: Naylor, T

Revisiting the accretion disc spectra of Dwarf Novae and Novalike variables: implications for the standard disc model

Authors: Gabriella Zsidi, C. J. Nixon, T. Naylor, J. E. Pringle

Abstract: Accretion discs are fundamental to much of astronomy. They can occur around stars both young and old, around compact objects they provide a window into the extremes of physics, and around supermassive black holes in galaxy centres they generate spectacular luminosities that can outshine the entire galaxy. However, our understanding of the inner workings of accretion discs remains far from complete… ▽ More Accretion discs are fundamental to much of astronomy. They can occur around stars both young and old, around compact objects they provide a window into the extremes of physics, and around supermassive black holes in galaxy centres they generate spectacular luminosities that can outshine the entire galaxy. However, our understanding of the inner workings of accretion discs remains far from complete. Here we revisit a conundrum in the observations of some of the simplest accreting systems; the Cataclysmic Variables (CVs). The high-accretion-rate states of (non-magnetic) CVs can be divided into the short-lived outbursts ($\sim$ a week) typical of dwarf novae (DNe) and the long-lived (and sometimes perpetual) high states of nova-like (NL) CVs. Since both sorts of high-state occur in approximately steady-state accretion discs with similar properties and accretors, we would expect them to display similar spectral energy distributions. However, previous analyses based on UV spectra from the {\it International Ultraviolet Explorer} have shown that their spectral energy distributions are different. We perform a re-analysis of the data using up to date calibrations and distance (and thus dereddening) estimates to test whether this difference persists and whether it is statistically significant over the sample. We find that it does persist and it is statistically significant. We propose routes to investigating this discrepancy further and discuss the implications this has for other accreting systems, such as X-ray binaries, Active Galactic Nuclei and protoplanetary discs. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: Accepted for publication in MNRAS. 12 pages, 8 figures

The magnetically quiet solar surface dominates HARPS-N solar RVs during low activity

Authors: Ben S. Lakeland, Tim Naylor, Raphaëlle Haywood, Nadège Meunier, Federica Rescigno, Shweta Dalal, Annelies Mortier, Samantha J. Thompson, Andrew Collier Cameron, Xavier Dumusque, Mercedes López-Morales, Francesco Pepe, Ken Rice, Alessandro Sozzetti, Stéphane Udry, Eric Ford, Adriano Ghedina, Marcello Lodi

Abstract: Using images from the Helioseismic and Magnetic Imager aboard the \textit{Solar Dynamics Observatory} (SDO/HMI), we extract the radial-velocity (RV) signal arising from the suppression of convective blue-shift and from bright faculae and dark sunspots transiting the rotating solar disc. We remove these rotationally modulated magnetic-activity contributions from simultaneous radial velocities obser… ▽ More Using images from the Helioseismic and Magnetic Imager aboard the \textit{Solar Dynamics Observatory} (SDO/HMI), we extract the radial-velocity (RV) signal arising from the suppression of convective blue-shift and from bright faculae and dark sunspots transiting the rotating solar disc. We remove these rotationally modulated magnetic-activity contributions from simultaneous radial velocities observed by the HARPS-N solar feed to produce a radial-velocity time series arising from the magnetically quiet solar surface (the 'inactive-region radial velocities'). We find that the level of variability in the inactive-region radial velocities remains constant over the almost 7 year baseline and shows no correlation with well-known activity indicators. With an RMS of roughly 1 m/s, the inactive-region radial-velocity time series dominates the total RV variability budget during the decline of solar cycle 24. Finally, we compare the variability amplitude and timescale of the inactive-region radial velocities with simulations of supergranulation. We find consistency between the inactive-region radial-velocity and simulated time series, indicating that supergranulation is a significant contribution to the overall solar radial velocity variability, and may be the main source of variability towards solar minimum. This work highlights supergranulation as a key barrier to detecting Earth twins. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 12 pages, 11 figures, accepted to MNRAS

Observational Bias and Young Massive Cluster Characterisation II. Can Gaia accurately observe young clusters and associations?

Authors: Anne S. M. Buckner, Tim Naylor, Clare L. Dobbs, Steven Rieder, Thomas J. R. Bending

Abstract: Observations of clusters suffer from issues such as completeness, projection effects, resolving individual stars and extinction. As such, how accurate are measurements and conclusions are likely to be? Here, we take cluster simulations (Westerlund2- and Orion- type), synthetically observe them to obtain luminosities, accounting for extinction and the inherent limits of Gaia, then place them within… ▽ More Observations of clusters suffer from issues such as completeness, projection effects, resolving individual stars and extinction. As such, how accurate are measurements and conclusions are likely to be? Here, we take cluster simulations (Westerlund2- and Orion- type), synthetically observe them to obtain luminosities, accounting for extinction and the inherent limits of Gaia, then place them within the real Gaia DR3 catalogue. We then attempt to rediscover the clusters at distances of between 500pc and 4300pc. We show the spatial and kinematic criteria which are best able to pick out the simulated clusters, maximising completeness and minimising contamination. We then compare the properties of the 'observed' clusters with the original simulations. We looked at the degree of clustering, the identification of clusters and subclusters within the datasets, and whether the clusters are expanding or contracting. Even with a high level of incompleteness (e.g. $<2\%$ stellar members identified), similar qualitative conclusions tend to be reached compared to the original dataset, but most quantitative conclusions are likely to be inaccurate. Accurate determination of the number, stellar membership and kinematic properties of subclusters, are the most problematic to correctly determine, particularly at larger distances due to the disappearance of cluster substructure as the data become more incomplete, but also at smaller distances where the misidentification of asterisms as true structure can be problematic. Unsurprisingly, we tend to obtain better quantitative agreement of properties for our more massive Westerlund2-type cluster. We also make optical style images of the clusters over our range of distances. △ Less

Submitted 13 November, 2023; v1 submitted 31 October, 2023; originally announced October 2023.

Comments: 19 pages, 10 figures, 6 tables. Accepted for publication in MNRAS

A trio of month long flares in the nova-like variable V704 And

Authors: Gabriella Zsidi, C. J. Nixon, T. Naylor, J. E. Pringle, K. L. Page

Abstract: We present the discovery of an unusual set of flares in the nova-like variable V704 And. Using data from AAVSO, ASAS-SN, and ZTF, of the nova-like variable V704 And, we have discovered a trio of brightening events that occurred during the high state. These events elevate the optical brightness of the source from $\sim13.5$ magnitude to $\sim12.5$ magnitude. The events last for roughly a month, and… ▽ More We present the discovery of an unusual set of flares in the nova-like variable V704 And. Using data from AAVSO, ASAS-SN, and ZTF, of the nova-like variable V704 And, we have discovered a trio of brightening events that occurred during the high state. These events elevate the optical brightness of the source from $\sim13.5$ magnitude to $\sim12.5$ magnitude. The events last for roughly a month, and exhibit the unusual shape of a slow rise and faster decay. Just after the third event we obtained data from regular monitoring with Swift, although by this time the flares had ceased and the source returned to its pre-flare level of activity in the high-state. The Swift observations confirm that during the high-state the source is detectable in the X-rays, and provide simultaneous UV and optical fluxes. As the source is already in the high-state prior to the flares, and thus the disc is expected to already be in the high-viscosity state, we conclude that the driver of the variations must be changes in the mass transfer rate from the companion star and we discuss possible mechanisms for such short-timescale mass transfer variations to occur. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 5 pages + appendix. Accepted for publication in A&A Letters

Journal ref: A&A 679, L3 (2023)

Classification of Chandra X-ray Sources in Cygnus OB2

Authors: Vinay L. Kashyap, Mario G. Guarcello, Nicholas J. Wright, Jeremy J. Drake, Ettore Flaccomio, Tom L. Aldcroft, Juan F. Albacete Colombo, Kevin Briggs, Francesco Damiani, Janet E. Drew, Eduardo L. Martin, Giusi Micela, Tim Naylor, Salvatore Sciortino

Abstract: We have devised a predominantly Naive Bayes method to classify the optical/IR matches to X-ray sources detected by Chandra in the Cygnus OB2 association into foreground, member, and background objects. We employ a variety of X-ray, optical, and infrared characteristics to construct likelihoods using training sets defined by well-measured sources. Combinations of optical photometry from SDSS (riz)… ▽ More We have devised a predominantly Naive Bayes method to classify the optical/IR matches to X-ray sources detected by Chandra in the Cygnus OB2 association into foreground, member, and background objects. We employ a variety of X-ray, optical, and infrared characteristics to construct likelihoods using training sets defined by well-measured sources. Combinations of optical photometry from SDSS (riz) and IPHAS (riHa), IR magnitudes from UKIDSS and 2MASS (JHK), X-ray quantiles and hardness ratios, and estimates of extinction Av are used to compute the relative probabilities that a given source belongs to one of the classes. We use Principal Component Analysis of photometric magnitude combinations to isolate the best axes for classification. We incorporate measurement errors into the classification. We evaluate the accuracy of the classification by inspection and reclassify a number of sources based on IR magnitudes, presence of disks, and X-ray spectral hardness. We also consider systematic errors due to extinction. We find that about 6100 objects are association members, 1400 are background, and 500 are foreground objects. The overall classification accuracy is 95%. △ Less

Submitted 11 June, 2023; originally announced June 2023.

Comments: 27 pages, 23 figures, 6 tables; accepted for publication in ApJS. Full Table 3 is in Zenodo at https://zenodo.org/record/8025756

MSC Class: 85A35 (Primary); 85A15; 62H30; 62P99 (Secondary)

A Naive Bayes Classifier for identifying Class II YSOs

Authors: Andrew J. Wilson, Ben S. Lakeland, Tom J. Wilson, Tim Naylor

Abstract: A naive Bayes classifier for identifying Class II YSOs has been constructed and applied to a region of the Northern Galactic Plane containing 8 million sources with good quality Gaia EDR3 parallaxes. The classifier uses the five features: Gaia $G$-band variability, WISE mid-infrared excess, UKIDSS and 2MASS near-infrared excess, IGAPS H$α$ excess and overluminosity with respect to the main sequenc… ▽ More A naive Bayes classifier for identifying Class II YSOs has been constructed and applied to a region of the Northern Galactic Plane containing 8 million sources with good quality Gaia EDR3 parallaxes. The classifier uses the five features: Gaia $G$-band variability, WISE mid-infrared excess, UKIDSS and 2MASS near-infrared excess, IGAPS H$α$ excess and overluminosity with respect to the main sequence. A list of candidate Class II YSOs is obtained by choosing a posterior threshold appropriate to the task at hand, balancing the competing demands of completeness and purity. At a threshold posterior greater than 0.5 our classifier identifies 6504 candidate Class II YSOs. At this threshold we find a false positive rate around 0.02 per cent and a true positive rate of approximately 87 per cent for identifying Class II YSOs. The ROC curve rises rapidly to almost one with an area under the curve around 0.998 or better, indicating the classifier is efficient at identifying candidate Class II YSOs. Our map of these candidates shows what are potentially three previously undiscovered clusters or associations. When comparing our results to published catalogues from other young star classifiers, we find between one quarter and three quarters of high probability candidates are unique to each classifier, telling us no single classifier is finding all young stars. △ Less

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

Comments: 38 pages, 28 figures, 15 tables. Accepted for publication in MNRAS

Journal ref: MNRAS vol 521, p354-388 (2023)

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

Authors: Shoko Jin, Scott C. Trager, Gavin B. Dalton, J. Alfonso L. Aguerri, J. E. Drew, Jesús Falcón-Barroso, Boris T. Gänsicke, Vanessa Hill, Angela Iovino, Matthew M. Pieri, Bianca M. Poggianti, D. J. B. Smith, Antonella Vallenari, Don Carlos Abrams, David S. Aguado, Teresa Antoja, Alfonso Aragón-Salamanca, Yago Ascasibar, Carine Babusiaux, Marc Balcells, R. Barrena, Giuseppina Battaglia, Vasily Belokurov, Thomas Bensby, Piercarlo Bonifacio , et al. (190 additional authors not shown)

Abstract: WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrogr… ▽ More WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator. △ Less

Submitted 31 October, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

ABORAS: polarimetric, 10cm/s RV observations of the Sun as a star

Authors: Casper Farret Jentink, Annelies Mortier, Frans Snik, Patrick Dorval, Samantha J. Thompson, Ramon Navarro, Tim Naylor

Abstract: We present a description of A dual-Beam pOlarimetric Robotic Aperture for the Sun (ABORAS), to serve as a Solar input with a dedicated Stokes V polarimeter for the HARPS3 high-resolution spectrograph. ABORAS has three main science drivers: trying to understand the physics behind stellar variability, tracking the long-term stability of HARPS3, and serve as a benchmark for Earth-sized exoplanet dete… ▽ More We present a description of A dual-Beam pOlarimetric Robotic Aperture for the Sun (ABORAS), to serve as a Solar input with a dedicated Stokes V polarimeter for the HARPS3 high-resolution spectrograph. ABORAS has three main science drivers: trying to understand the physics behind stellar variability, tracking the long-term stability of HARPS3, and serve as a benchmark for Earth-sized exoplanet detection with HARPS3 by injecting an Earth RV signal into the data. By design, ABORAS will (together with the HARPS3 instrument) be able to measure 10cm/s variations in RV of the integrated Solar disk and detect integrated magnetic field levels at sub 1 Gauss level through circularly polarized light. △ Less

Submitted 18 July, 2022; v1 submitted 11 July, 2022; originally announced July 2022.

Comments: 8 pages, 5 figures, SPIE Proceedings pre-print, draft version 2

Towards an understanding of YSO variability: A multi-wavelength analysis of bursting, dipping, and symmetrically varying light curves of disc-bearing YSOs

Authors: Ben S. Lakeland, Tim Naylor

Abstract: Using simultaneous optical and infrared light curves of disc-bearing young stars in NGC 2264, we perform the first multi-wavelength structure function study of YSOs. We find that dippers have larger variability amplitudes than bursters and symmetric variables at all timescales longer than a few hours. By analysing optical-infrared colour time-series, we also find that the variability in the burste… ▽ More Using simultaneous optical and infrared light curves of disc-bearing young stars in NGC 2264, we perform the first multi-wavelength structure function study of YSOs. We find that dippers have larger variability amplitudes than bursters and symmetric variables at all timescales longer than a few hours. By analysing optical-infrared colour time-series, we also find that the variability in the bursters is systematically less chromatic at all timescales than the other variability types. We propose a model of YSO variability in which symmetric, bursting, and dipping behaviour is observed in systems viewed at low, intermediate, and high inclinations, respectively. We argue that the relatively short thermal timescale for the disc can explain the fact that the infrared light curves for bursters are more symmetric than their optical counterparts. △ Less

Submitted 26 May, 2022; originally announced May 2022.

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

Observational Bias and Young Massive Cluster Characterisation I. 2D Perspective Effects

Authors: Anne S. M. Buckner, Kong You Liow, Clare L. Dobbs, Tim Naylor, Steven Rieder

Abstract: Understanding the formation and evolution of high mass star clusters requires comparisons between theoretical and observational data to be made. Unfortunately, while the full phase space of simulated regions is available, often only partial 2D spatial and kinematic data is available for observed regions. This raises the question as to whether cluster parameters determined from 2D data alone are re… ▽ More Understanding the formation and evolution of high mass star clusters requires comparisons between theoretical and observational data to be made. Unfortunately, while the full phase space of simulated regions is available, often only partial 2D spatial and kinematic data is available for observed regions. This raises the question as to whether cluster parameters determined from 2D data alone are reliable and representative of clusters real parameters and the impact of line-of-sight orientation. In this paper we derive parameters for a simulated cluster formed from a cloud-cloud collision with the full 6D phase space, and compare them with those derived from three different 2D line-of-sight orientations for the cluster. We show the same qualitative conclusions can be reached when viewing clusters in 2D versus 3D, but that drawing quantitative conclusions when viewing in 2D is likely to be inaccurate. The greatest divergence occurs in the perceived kinematics of the cluster, which in some orientations appears to be expanding when the cluster is actually contracting. Increases in the cluster density compounds pre-existing perspective issues, reducing the relative accuracy and consistency of properties derived from different orientations. This is particularly problematic for determination of the number, and membership, of subclusters present in the cluster. We find the fraction of subclusters correctly identified in 2D decreases as the cluster evolves, reaching less than 3.4% at the evolutionary end point for our cluster. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: 16 pages, 10 figures, 1 table. Accepted for publication in MNRAS

Spatial Statistics in Star Forming Regions: Is Star Formation Driven By Column Density Alone?

Authors: Brendan Retter, Jennifer Hatchell, Tim Naylor

Abstract: Star formation is known to occur more readily where more raw materials are available. This is often expressed by a 'Kennicutt-Schmidt' relation where the surface density of Young Stellar Objects (YSOs) is proportional to column density to some power, $μ$. The aim of this work was to determine if column density alone is sufficient to explain the locations of Class 0/I YSOs within Serpens South, Ser… ▽ More Star formation is known to occur more readily where more raw materials are available. This is often expressed by a 'Kennicutt-Schmidt' relation where the surface density of Young Stellar Objects (YSOs) is proportional to column density to some power, $μ$. The aim of this work was to determine if column density alone is sufficient to explain the locations of Class 0/I YSOs within Serpens South, Serpens Core, Ophiuchus, NGC1333 and IC348, or if there is clumping or avoidance that would point to additional influences on the star formation. Using the O-ring test as a summary statistic, 95 per cent confidence envelopes were produced for different values of $μ$ from probability models made using the Herschel column density maps. The YSOs were tested against four distribution models: the best-estimate of $μ$ for the region, $μ=0$ above a minimum column density threshold and zero probability elsewhere, $μ=1$, and the power-law that best represents the five regions as a collective, $μ=2.05 \pm 0.20$. Results showed that $μ=2.05$ model was consistent with the majority of regions and, for those regions, the spatial distribution of YSOs at a given column density is consistent with being random. Serpens South and NGC1333 rejected the $μ= 2.05$ model on small scales of $\sim 0.15 \mathrm{pc}$ which implies that small-scale interactions may be necessary to improve the model. On scales above 0.15 pc, the positions of YSOs in all five regions can be well described using column density alone. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: 19 pages, 13 figures. MNRAS accepted

High resolution H-alpha imaging of the Northern Galactic Plane, and the IGAPS images database

Authors: R. Greimel, J. E. Drew, M. Monguió, R. P. Ashley, G. Barentsen, J. Eislöffel, A. Mampaso, R. A. H. Morris, T. Naylor, C. Roe, L. Sabin, B. Stecklum, N. J. Wright, P. J. Groot, M. J. Irwin, M. J. Barlow, C. Fariña, A. Fernández-Martín, Q. A. Parker, S. Phillipps, S. Scaringi, A. A. Zijlstra

Abstract: The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. These capture the entire northern Galactic plane within the Galactic coordinate range, -5<b<+5 deg. and 30<l<215 deg. From the beginning, the incorporation of narrowband H-alpha imaging has been a unique and dist… ▽ More The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. These capture the entire northern Galactic plane within the Galactic coordinate range, -5<b<+5 deg. and 30<l<215 deg. From the beginning, the incorporation of narrowband H-alpha imaging has been a unique and distinctive feature of this effort. Alongside a focused discussion of the nature and application of the H-alpha data, we present the IGAPS world-accessible database of images for all 5 survey filters, i, r, g, U-RGO and narrowband H-alpha, observed on a pixel scale of 0.33 arcsec and at an effective (median) angular resolution of 1.1 to 1.3 arcsec. The background, noise, and sensitivity characteristics of the narrowband H-alpha filter images are outlined. Typical noise levels in this band correspond to a surface brightness at full one-arcsec resolution of around 2e-16 erg/cm2/s/arcsec2. Illustrative applications of the H-alpha data to planetary nebulae and Herbig-Haro objects are outlined and, as part of a discussion of mosaicking technique, we present a very large background-subtracted narrowband mosaic of the supernova remnant, Simeis 147. Finally we lay out a method that exploits the database via an automated selection of bright ionized diffuse interstellar emission targets for the coming generation of wide-field massive-multiplex spectrographs. Two examples of the diffuse H-alpha maps output from this selection process are presented and compared with previously published data. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: 21 pages, 14 main-text figures, 11 appendix figures. Images database and other supplementary items mentioned in the paper are available from http://www.igapsimages.org

Journal ref: A&A 655, A49 (2021)

The JCMT Transient Survey: Four Year Summary of Monitoring the Submillimeter Variability of Protostars

Authors: Yong-Hee Lee, Doug Johnstone, Jeong-Eun Lee, Gregory Herczeg, Steve Mairs, Carlos Contreras-Peña, Jennifer Hatchell, Tim Naylor, Graham S. Bell, Tyler L. Bourke, Colton Broughton, Logan Francis, Aashish Gupta, Daniel Harsono, Sheng-Yuan Liu, Geumsook Park, Spencer Plovie, Gerald H. Moriarty-Schieven, Aleks Scholz, Tanvi Sharma, Paula Stella Teixeira, Yao-Te Wang, Yuri Aikawa, Geoffrey C. Bower, Huei-Ru Vivien Chen , et al. (27 additional authors not shown)

Abstract: We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources.… ▽ More We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted for publication in the Astrophysical Journal

Statistical Fitting of Evolution Models to Rotation Rates of Sun-Like Stars

Authors: Angela A. Breimann, Sean P. Matt, Tim Naylor

Abstract: We apply for the first time a two-dimensional fitting statistic, $τ^2$, to rotational evolution models (REMs) of stars (0.1 to 1.3 $M_{\odot}$) on the period-mass plane. The $τ^2$ statistic simultaneously considers all cluster rotation data to return a goodness of fit, allowing for data-driven improvement of REMs. We construct data sets for Upper Sco, the Pleiades and Praesepe, to which we tune ou… ▽ More We apply for the first time a two-dimensional fitting statistic, $τ^2$, to rotational evolution models (REMs) of stars (0.1 to 1.3 $M_{\odot}$) on the period-mass plane. The $τ^2$ statistic simultaneously considers all cluster rotation data to return a goodness of fit, allowing for data-driven improvement of REMs. We construct data sets for Upper Sco, the Pleiades and Praesepe, to which we tune our REMs. We use consistently determined stellar masses (calculated by matching $K_\textrm{s}$ magnitudes to isochrones) and literature rotation periods. As a first demonstration of the $τ^2$ statistic, we find the best-fitting gyrochronology age for Praesepe, which is in good agreement with the literature. We then systematically vary three parameters which determine the dependence of our stellar wind torque law on Rossby number in the saturated and unsaturated regimes, and the location of the transition between the two. By minimising $τ^2$, we find best-fit values for each parameter. These values vary slightly between clusters, mass determinations and initial conditions, highlighting the precision of $τ^2$ and its potential for constraining REMs, gyrochronology, and understanding of stellar physics. Our resulting REMs, which implement the best-possible fitting form of a broken power-law torque, are statistically improved on previous REMs using similar formulations, but still do not simultaneously describe the observed rotation distributions of the lowest masses, which have both slow and fast rotators by the Praesepe age, and the shape of the converged sequence for higher masses. Further complexity in the REMs is thus required to accurately describe the data. △ Less

Submitted 24 March, 2021; originally announced March 2021.

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

Young Faithful: The Eruptions of EC 53 as It Cycles through Filling and Draining the Inner Disk

Authors: Yong-Hee Lee, Doug Johnstone, Jeong-Eun Lee, Gregory Herczeg, Steve Mairs, Watson Varricatt, Klaus W. Hodapp, Tim Naylor, Carlos Contreras Peña, Giseon Baek, Martin Haas, Rolf Chini, The JCMT Transient Team

Abstract: While young stellar objects sometimes undergo bursts of accretion, these bursts usually occur sporadically, making them challenging to study observationally and to explain theoretically. We build a schematic description of cyclical bursts of the young stellar object EC 53 using near-IR and sub-mm monitoring obtained over six cycles, each lasting $\approx530$ days. EC 53 brightens over $0.12$ yr by… ▽ More While young stellar objects sometimes undergo bursts of accretion, these bursts usually occur sporadically, making them challenging to study observationally and to explain theoretically. We build a schematic description of cyclical bursts of the young stellar object EC 53 using near-IR and sub-mm monitoring obtained over six cycles, each lasting $\approx530$ days. EC 53 brightens over $0.12$ yr by $0.3$ mag at 850 $μ$m, $2$ mag at 3.35 $μ$m, and $1.5$ mag at near-IR wavelengths, to a maximum luminosity consistent with an accretion rate of $\sim8\times10^{-6}$ M$_\odot$ yr$^{-1}$. The emission then decays with an e-folding timescale of $\approx0.74$ yr until the accretion rate is $\sim1\times10^{-6}$ M$_\odot$ yr$^{-1}$. The next eruption then occurs, likely triggered by the buildup of $\sim5\times10^{-6}$ M$_\odot$ of mass in the inner disk, enough that it becomes unstable and drains onto the star. Just before outburst, when the disk is almost replenished, the near-IR colors become redder, indicating an increase in the geometrical height of the disk by this mass buildup. The reddening disappears soon after the initial burst, as much of the mass is drained from the disk. We quantify physical parameters related to the accretion process in EC 53 by assuming an $α$-disk formulation, constrained by the observed disk properties and accretion rate. While we can only speculate about the possible trigger for these faithful eruptions, we hope that our quantified schematic will motivate theorists to test the hypothesized mechanisms that could cause the cyclical buildup and draining of mass in the inner disk. △ Less

Submitted 17 September, 2020; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: Accepted for publication in the Astrophysical Journal

Submillimetre Transient Science in the Next Decade: EAO Submillimetre Futures White Paper Series, 2019

Authors: Steve Mairs, Gregory Herczeg, Doug Johnstone, Jeong-Eun Lee, Simon Coude, Alexandra J. Tetarenko, Jenny Hatchell, Aleks Scholz, Bhavana Lalchand, Wen-Ping Chen, Carlos Contreras Pena, Tim Naylor, Kevin Lacaille, Peter Scicluna

Abstract: This white paper gives a brief summary of the time domain science that has been performed with the JCMT in recent years and highlights the opportunities for continuing work in this field over the next decade. The main focus of this document is the JCMT Transient Survey, a large program initiated in 2015 to measure the frequency and amplitude of variability events associated with protostars in near… ▽ More This white paper gives a brief summary of the time domain science that has been performed with the JCMT in recent years and highlights the opportunities for continuing work in this field over the next decade. The main focus of this document is the JCMT Transient Survey, a large program initiated in 2015 to measure the frequency and amplitude of variability events associated with protostars in nearby star-forming regions. After summarising the major accomplishments so far, an outline is given for extensions to the current survey, featuring a discussion on what will be possible with the new 850 micron camera that is expected to be installed in late 2022. We also discuss possible applications of submillimetre monitoring to active galactic nuclei, X-ray binaries, asymptotic giant branch stars, and flare stars. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: 11 pages, 4 figures; Part of the 2019 EAO Submillimetre Futures Paper Series

Characterising the i-band variability of YSOs over six orders of magnitude in timescale

Authors: Darryl J. Sergison, Tim Naylor, S. P. Littlefair, Cameron P. M. Bell, C. D. H. Williams

Abstract: We present an $i$-band photometric study of over 800 young stellar objects in the OB association Cep OB3b, which samples timescales from 1 minute to ten years. Using structure functions we show that on all timescales ($τ$) there is a monotonic decrease in variability from Class I to Class II through the transition disc (TD) systems to Class III, i.e. the more evolved systems are less variable. The… ▽ More We present an $i$-band photometric study of over 800 young stellar objects in the OB association Cep OB3b, which samples timescales from 1 minute to ten years. Using structure functions we show that on all timescales ($τ$) there is a monotonic decrease in variability from Class I to Class II through the transition disc (TD) systems to Class III, i.e. the more evolved systems are less variable. The Class Is show an approximately power-law increase ($τ^{0.8}$) in variability from timescales of a few minutes to ten years. The Class II, TDs and Class III systems show a qualitatively different behaviour with most showing a power-law increase in variability up to a timescale corresponding to the rotational period of the star, with little additional variability beyond that timescale. However, about a third of the Class IIs show lower overall variability, but their variability is still increasing at 10 years. This behaviour can be explained if all Class IIs have two primary components to their variability. The first is an underlying roughly power-law variability spectrum, which evidence from the infrared suggests is driven by accretion rate changes. The second component is approximately sinusoidal and results from the rotation of the star. We suggest that the systems with dominant longer-timescale variability have a smaller rotational modulation either because they are seen at low inclinations or have more complex magnetic field geometries. We derive a new way of calculating structure functions for large simulated datasets (the "fast structure function"), based on fast Fourier transforms. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: 23 pages and 22 figures. Accepted for publication in MNRAS. At time of submission the electronic tables are only available from https:// doi.org/10.24378/exe.2124

The young stellar content of the giant HII regions M8, G333.6-0.2, and NGC6357 with VLT/KMOS

Authors: M. C. Ramírez-Tannus, J. Poorta, A. Bik, L. Kaper, A. de Koter, J. De Ridder, H. Beuther, W. Brandner, B. Davies, M. Gennaro, D. Guo, T. Henning, H. Linz, T. Naylor, A. Pasquali, O. H. Ramírez-Agudelo, H. Sana

Abstract: Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations… ▽ More Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A. The arXiv version includes appendix C, which is an online only figure in A&A

Journal ref: A&A 633, A155 (2020)

Exploring the M-dwarf Luminosity--Temperature--Radius Relationships using Gaia DR2

Authors: Sam Morrell, Tim Naylor

Abstract: There is growing evidence that M-dwarf stars suffer radius inflation when compared to theoretical models, suggesting that models are missing some key physics required to completely describe stars at effective temperatures $(T_{\rm SED})$ less than about 4000K. The advent of Gaia DR2 distances finally makes available large datasets to determine the nature and extent of this effect. We employ an all… ▽ More There is growing evidence that M-dwarf stars suffer radius inflation when compared to theoretical models, suggesting that models are missing some key physics required to completely describe stars at effective temperatures $(T_{\rm SED})$ less than about 4000K. The advent of Gaia DR2 distances finally makes available large datasets to determine the nature and extent of this effect. We employ an all-sky sample, comprising of $>$15\,000 stars, to determine empirical relationships between luminosity, temperature and radius. This is accomplished using only geometric distances and multiwave-band photometry, by utilising a modified spectral energy distribution fitting method. The radii we measure show an inflation of $3 - 7\%$ compared to models, but no more than a $1 - 2\%$ intrinsic spread in the inflated sequence. We show that we are currently able to determine M-dwarf radii to an accuracy of $2.4\%$ using our method. However, we determine that this is limited by the precision of metallicity measurements, which contribute $1.7\%$ to the measured radius scatter. We also present evidence that stellar magnetism is currently unable to explain radius inflation in M-dwarfs. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS 3 August 2019. 19 pages, 19 figures. Online supplementary material currently available at https://github.com/sammorrell/mn19-supplementary-material and https://doi.org/10.24378/exe.1683

A Volume Limited Sample of Cataclysmic Variables from $\mathit{Gaia}$ DR2: Space Density and Population Properties

Authors: A. F. Pala, B. T. Gänsicke, E. Breedt, C. Knigge, J. J. Hermes, N. P. Gentile Fusillo, M. A. Hollands, T. Naylor, I. Pelisoli, M. R. Schreiber, S. Toonen, A. Aungwerojwit, E. Cukanovaite, E. Dennihy, C. J. Manser, M. L. Pretorius, S. Scaringi, O. Toloza

Abstract: We present the first volume-limited sample of cataclysmic variables (CVs), selected using the accurate parallaxes provided by the second data release (DR2) of the ESA $\mathit{Gaia}$ space mission. The sample is composed of 42 CVs within $150\,$pc, including two new systems discovered using the $\mathit{Gaia}$ data, and is $(77 \pm 10)\,$per cent complete. We use this sample to study the intrinsic… ▽ More We present the first volume-limited sample of cataclysmic variables (CVs), selected using the accurate parallaxes provided by the second data release (DR2) of the ESA $\mathit{Gaia}$ space mission. The sample is composed of 42 CVs within $150\,$pc, including two new systems discovered using the $\mathit{Gaia}$ data, and is $(77 \pm 10)\,$per cent complete. We use this sample to study the intrinsic properties of the Galactic CV population. In particular, the CV space density we derive, $ρ=(4.8^{+0.6}_{-0.8}) \times10^{-6}\,\mathrm{pc}^{-3}$, is lower than predicted by most binary population synthesis studies. We also find a low fraction of period bounce CVs, seven per cent, and an average white dwarf mass of $\langle M_\mathrm{WD} \rangle = (0.83 \pm 0.17)\,\mathrm{M}_\odot$. Both findings confirm previous results, ruling out the presence of observational biases affecting these measurements, as has been suggested in the past. The observed fraction of period bounce CVs falls well below theoretical predictions, by at least a factor of five, and remains one of the open problems in the current understanding of CV evolution. Conversely, the average white dwarf mass supports the presence of additional mechanisms of angular momentum loss that have been accounted for in the latest evolutionary models. The fraction of magnetic CVs in the $150\,$pc sample is remarkably high at $36\,$per cent. This is in striking contrast with the absence of magnetic white dwarfs in the detached population of CV progenitors, and underlines that the evolution of magnetic systems has to be included in the next generation of population models. △ Less

Submitted 6 May, 2020; v1 submitted 30 July, 2019; originally announced July 2019.

Comments: Published in MNRAS. Definitive version. A typo in Table 6 (space density of novalikes) has been corrected

Spatial Statistics in Star Forming Regions: Testing the Limits of Randomness

Authors: Brendan Retter, Jennifer Hatchell, Tim Naylor

Abstract: Observational studies of star formation reveal spatial distributions of Young Stellar Objects (YSOs) that are `snapshots' of an ongoing star formation process. Using methods from spatial statistics it is possible to test the likelihood that a given distribution process could produce the observed patterns of YSOs. The aim of this paper is to determine the usefulness of the spatial statistics tests… ▽ More Observational studies of star formation reveal spatial distributions of Young Stellar Objects (YSOs) that are `snapshots' of an ongoing star formation process. Using methods from spatial statistics it is possible to test the likelihood that a given distribution process could produce the observed patterns of YSOs. The aim of this paper is to determine the usefulness of the spatial statistics tests Diggle's G function (G), the `free-space' function (F), Ripley's K and O-ring for application to astrophysical data. The spatial statistics tests were applied to simulated data containing 2D Gaussian clusters projected on random distributions of stars. The number of stars within the Gaussian cluster and number of background stars were varied to determine the tests' ability to reject complete spatial randomness (CSR) with changing signal-to-noise. The best performing test was O-ring optimised with overlapping logarithmic bins, closely followed by Ripley's K. The O-ring test is equivalent to the 2-point correlation function. Both F and G (and the minimum spanning tree, of which G is a subset) performed significantly less well, requiring a cluster with a factor of two higher signal-to-noise in order to reject CSR consistently. We demonstrate the tests on example astrophysical datasets drawn from the Spitzer catalogue. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: 15 pages, 14 figures, 2 tables, accepted for publication in MNRAS

Determining the recurrence timescale of long-lasting YSO outbursts

Authors: C. Contreras Peña, T. Naylor, S. Morrell

Abstract: We have determined the rate of large accretion events in class I and II young stellar objects (YSOs) by comparing the all-sky digitised photographic plate surveys provided by SuperCOSMOS with the latest data release from Gaia (DR2). The long mean baseline of 55 years along with a large sample of class II YSOs ($\simeq$15,000) allows us to study approximately 1 million YSO-years. We find 139 object… ▽ More We have determined the rate of large accretion events in class I and II young stellar objects (YSOs) by comparing the all-sky digitised photographic plate surveys provided by SuperCOSMOS with the latest data release from Gaia (DR2). The long mean baseline of 55 years along with a large sample of class II YSOs ($\simeq$15,000) allows us to study approximately 1 million YSO-years. We find 139 objects with $ΔR\geq1$~mag, most of which are found at amplitudes between 1 and 3 mag. The majority of YSOs in this group show irregular variability or long-lasting fading events, which is best explained as hot spots due to accretion or by variable extinction. There is a tail of YSOs at $ΔR\geq3$~mag and they seem to represent a different population. Surprisingly many objects in this group show high-amplitude irregular variability over timescales shorter than 10 years, in contrast with the view that high-amplitude objects always have long outbursts. However, we find 6 objects that are consistent with undergoing large, long lasting accretion events, 3 of them previously unknown. This yields an outburst recurrence timescale of 112 kyr, with a 68\% confidence interval [74 to 180] kyr. This represents the first robust determination of the outburst rate in class II YSOs and shows that YSOs in their planet-forming stage do in fact undergo large accretion events, and with timescales of $\simeq$100,000 years. In addition, we find that outbursts in the class II stage are $\simeq$10 times less frequent than during the class I stage. △ Less

Submitted 8 April, 2019; originally announced April 2019.

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

Gaia 17bpi: An FU Ori Type Outburst

Authors: Lynne A. Hillenbrand, Carlos Contreras Peña, Sam Morrell, Tim Naylor, Michael A. Kuhn, Roc M. Cutri, Luisa M. Rebull, Simon Hodgkin, Dirk Froebrich, Amy K. Mainzer

Abstract: We report on the source Gaia 17bpi and identify it as a new, ongoing FU Ori type outburst, associated with a young stellar object. The optical lightcurve from Gaia exhibited a 3.5 mag rise with the source appearing to plateau in mid/late 2018. Mid-infrared observations from NEOWISE also show a $>$3 mag rise that occurred in two stages, with the second one coincident with the optical brightening, a… ▽ More We report on the source Gaia 17bpi and identify it as a new, ongoing FU Ori type outburst, associated with a young stellar object. The optical lightcurve from Gaia exhibited a 3.5 mag rise with the source appearing to plateau in mid/late 2018. Mid-infrared observations from NEOWISE also show a $>$3 mag rise that occurred in two stages, with the second one coincident with the optical brightening, and the first one preceding the optical brightening by $\sim$1.5 years. We model the outburst as having started between October and December of 2014. This wavelength-dependent aspect of young star accretion-driven outbursts has never been documented before. Both the mid-infrared and the optical colors of the object become bluer as the outburst proceeds. Optical spectroscopic characteristics in the outburst phase include: a GK-type absorption spectrum, strong wind/outflow in e.g. Mgb, NaD, H$α$, KI, OI, and CaII profiles, and detection of LiI 6707 Å. The infrared spectrum in the outburst phase is similar to that of an M-type spectrum, notably exhibiting prominent $H_2O$ and $^{12}$CO (2-0) bandhead absorption in the K-band, and likely HeI wind in the Y-band. The new FU Ori source Gaia 17bpi is associated with a little-studied dark cloud in the galactic plane, located at a distance of 1.27 kpc. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: Accepted to ApJ

The Diverse Science Return from a Wide-Area Survey of the Galactic Plane

Authors: R. A. Street, M. B. Lund, S. Khakpash, M. Donachie, W. A. Dawson, N. Golovich, L. Wyrzykowski, P. Szkody, T. Naylor, M. Penny, N. Rattenbury, M. Dall'Ora, W. I. Clarkson, D. Bennett, J. Pepper, M. Rabus, M. P. G. Hundertmark, Y. Tsapras, R. Di Stefano, S. Ridgway, M. Liu, E. Lin

Abstract: The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would gre… ▽ More The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would greatly benefit from an LSST survey that included the Galactic Plane, Magellanic Clouds and Bulge at a cadence of 2-3 d. Particular highlights include measuring the mass spectrum of black holes, and mapping the population of exoplanets in the Galaxy in relation to variations in star forming environments. But the same survey data will provide a goldmine for a wide range of science, and we explore possible survey strategies which maximize the scientific return for a number of fields including young stellar objects, cataclysmic variables and Neptune Trojans. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: White Paper submitted in response to the 2018 call for LSST cadence proposals

A Contaminant-Free Catalogue of Gaia DR2-WISE Galactic Plane Matches: Including the Effects of Crowding in the Cross-Matching of Photometric Catalogues

Authors: Tom J Wilson, Tim Naylor

Abstract: Faint, hidden contaminants in the point-spread functions (PSFs) of stars cause shifts to their measured positions. Wilson & Naylor (2017) showed failing to account for these shifts can lead to a drastic decrease in the number of returned catalogue matches in crowded fields. Here we highlight the effect these perturbations have on cross-matching, for matches between Gaia DR2 and WISE stars in a cro… ▽ More Faint, hidden contaminants in the point-spread functions (PSFs) of stars cause shifts to their measured positions. Wilson & Naylor (2017) showed failing to account for these shifts can lead to a drastic decrease in the number of returned catalogue matches in crowded fields. Here we highlight the effect these perturbations have on cross-matching, for matches between Gaia DR2 and WISE stars in a crowded Galactic plane region. Applying the uncertainties as quoted to Gaussian-based astrometric uncertainty functions (AUFs) can lead, in dense Galactic fields, to only matching 55% of the counterparts. We describe the construction of empirical descriptions for AUFs, building on the cross-matching method of Wilson & Naylor (2018), utilising the magnitudes of both catalogues to discriminate between true and false counterparts. We apply the improved cross-matching method to the Galactic plane |b| < 10. We provide the most likely counterpart matches and their respective probabilities. We also analyse several cases to verify the robustness of the results, highlighting some important caveats and considerations. Finally, we discuss the effect PSF resolution has by comparing the intra-catalogue nearest neighbour separation distributions of a sample of likely contaminated WISE objects and their corresponding Spitzer counterpart. We show that some WISE contaminants are resolved in Spitzer, with smaller intra-catalogue separations. We have highlighted the effect contaminant stars have on WISE, but it is important for all photometric catalogues, playing an important role in the next generation of surveys, such as LSST. △ Less

Submitted 31 August, 2018; originally announced September 2018.

Comments: 22 pages, 16 figures, 3 tables. Accepted for publication in MNRAS. Dataset of Gaia-WISE matches described in paper is available from the University of Exeter ORE repository in FITS format, and will appear on the CDS once the paper is assigned Volume and Page numbers

The Massive Star-Forming Regions Omnibus X-Ray Catalog, Second Installment

Authors: Leisa K. Townsley, Patrick S. Broos, Gordon P. Garmire, Gemma E. Anderson, Eric D. Feigelson, Tim Naylor, Matthew S. Povich

Abstract: We present the second installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC2), a compilation of X-ray point sources detected in Chandra/ACIS observations of 16 Galactic MSFRs and surrounding fields. MOXC2 includes 13 ACIS mosaics, three containing a pair of unrelated MSFRs at different distances, with a total catalog of 18,396 point sources. The MSFRs sampled range ov… ▽ More We present the second installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC2), a compilation of X-ray point sources detected in Chandra/ACIS observations of 16 Galactic MSFRs and surrounding fields. MOXC2 includes 13 ACIS mosaics, three containing a pair of unrelated MSFRs at different distances, with a total catalog of 18,396 point sources. The MSFRs sampled range over distances of 1.3 kpc to 6 kpc and populations varying from single massive protostars to the most massive Young Massive Cluster known in the Galaxy. By carefully detecting and removing X-ray point sources down to the faintest statistically-significant limit, we facilitate the study of the remaining unresolved X-ray emission. Through comparison with mid-infrared images that trace photon-dominated regions and ionization fronts, we see that the unresolved X-ray emission is due primarily to hot plasmas threading these MSFRs, the result of feedback from the winds and supernovae of massive stars. The 16 MSFRs studied in MOXC2 more than double the MOXC1 sample, broadening the parameter space of ACIS MSFR explorations and expanding Chandra's substantial contribution to contemporary star formation science. △ Less

Submitted 13 February, 2018; originally announced February 2018.

Comments: 60 pages, 24 figures. Accepted by ApJS. Figures have been significantly downgraded to meet arXiv requirements. See https://doi.org/10.5281/zenodo.1067748 for higher-resolution versions of this paper and for supplemental data products

Improving Catalogue Matching By Supplementing Astrometry with Additional Photometric Information

Authors: Tom J. Wilson, Tim Naylor

Abstract: The matching of sources between photometric catalogues can lead to cases where objects of differing brightness are incorrectly assumed to be detections of the same source. The rejection of unphysical matches can be achieved through the inclusion of information about the sources' magnitudes. The method described here uses the additional photometric information from both catalogues in the process of… ▽ More The matching of sources between photometric catalogues can lead to cases where objects of differing brightness are incorrectly assumed to be detections of the same source. The rejection of unphysical matches can be achieved through the inclusion of information about the sources' magnitudes. The method described here uses the additional photometric information from both catalogues in the process of accepting or rejecting counterparts, providing approximately a factor 10 improvement in Bayes' factor with its inclusion. When folding in the photometric information we avoid using prior astrophysical knowledge. Additionally, the method allows for the possibility of no counterparts to sources as well as the possibility that sources overlap multiple potential counterparts. We formally describe the probability of two sources being the same astrometric object, allowing systematic effects of astrometric perturbation (by, e.g., contaminant objects) to be accounted for. We apply the method to two cases. First, we test IPHAS-Gaia matches to compare the resulting matches in two catalogues of similar wavelength coverage but differing dynamical ranges. Second, we apply the method to matches between IPHAS and 2MASS and show that the method holds when considering two catalogues with approximately equal astrometric precision. We discuss the importance of including the magnitude information in each case. Additionally, we discuss extending the method to multiple catalogue matches through an iterative matching process. The method allows for the selection of high-quality matches by providing an overall probability for each pairing, giving the flexibility to choose stars known to be good matches. △ Less

Submitted 31 August, 2018; v1 submitted 13 October, 2017; originally announced October 2017.

Comments: 19 pages, 8 figures, 3 tables; accepted for publication in MNRAS. Updated equations 26 and 27 to fix minor typographical errors

The Effect of Unresolved Contaminant Stars on the Cross-Matching of Photometric Catalogues

Authors: Tom J. Wilson, Tim Naylor

Abstract: A fundamental process in astrophysics is the matching of two photometric catalogues. It is crucial that the correct objects be paired, and that their photometry does not suffer from any spurious additional flux. We compare the positions of sources in WISE, IPHAS, 2MASS, and APASS with Gaia DR1 astrometric positions. We find that the separations are described by a combination of a Gaussian distribu… ▽ More A fundamental process in astrophysics is the matching of two photometric catalogues. It is crucial that the correct objects be paired, and that their photometry does not suffer from any spurious additional flux. We compare the positions of sources in WISE, IPHAS, 2MASS, and APASS with Gaia DR1 astrometric positions. We find that the separations are described by a combination of a Gaussian distribution, wider than naively assumed based on their quoted uncertainties, and a large wing, which some authors ascribe to proper motions. We show that this is caused by flux contamination from blended stars not treated separately. We provide linear fits between the quoted Gaussian uncertainty and the core fit to the separation distributions. We show that at least one in three of the stars in the faint half of a given catalogue will suffer from flux contamination above the 1% level when the density of catalogue objects per PSF area is above approximately 0.005. This has important implications for the creation of composite catalogues. It is important for any closest neighbour matches as there will be a given fraction of matches that are flux contaminated, while some matches will be missed due to significant astrometric perturbation by faint contaminants. In the case of probability-based matching, this contamination affects the probability density function of matches as a function of distance. This effect results in up to 50% fewer counterparts being returned as matches, assuming Gaussian astrometric uncertainties for WISE-Gaia matching in crowded Galactic plane regions, compared with a closest neighbour match. △ Less

Submitted 29 March, 2017; v1 submitted 10 March, 2017; originally announced March 2017.

Comments: 10 pages, 6 figures; accepted for publication in the Monthly Notices of the Royal Astronomical Society

HARPS3 for a Roboticized Isaac Newton Telescope

Authors: Samantha J. Thompson, Didier Queloz, Isabelle Baraffe, Martyn Brake, Andrey Dolgopolov, Martin Fisher, Michel Fleury, Joost Geelhoed, Richard Hall, Jonay I. Gonzalez Hernandez, Rik ter Horst, Jan Kragt, Ramon Navarro, Tim Naylor, Francesco Pepe, Nikolai Piskunov, Rafael Rebolo, Louis Sander, Damien Segransan, Eugene Seneta, David Sing, Ignas Snellen, Frans Snik, Julien Spronck, Eric Stempels , et al. (3 additional authors not shown)

Abstract: We present a description of a new instrument development, HARPS3, planned to be installed on an upgraded and roboticized Isaac Newton Telescope by end-2018. HARPS3 will be a high resolution (R = 115,000) echelle spectrograph with a wavelength range from 380-690 nm. It is being built as part of the Terra Hunting Experiment - a future 10 year radial velocity measurement programme to discover Earth-l… ▽ More We present a description of a new instrument development, HARPS3, planned to be installed on an upgraded and roboticized Isaac Newton Telescope by end-2018. HARPS3 will be a high resolution (R = 115,000) echelle spectrograph with a wavelength range from 380-690 nm. It is being built as part of the Terra Hunting Experiment - a future 10 year radial velocity measurement programme to discover Earth-like exoplanets. The instrument design is based on the successful HARPS spectrograph on the 3.6m ESO telescope and HARPS-N on the TNG telescope. The main changes to the design in HARPS3 will be: a customised fibre adapter at the Cassegrain focus providing a stabilised beam feed and on-sky fibre diameter ~ 1.4 arcsec, the implementation of a new continuous flow cryostat to keep the CCD temperature very stable, detailed characterisation of the HARPS3 CCD to map the effective pixel positions and thus provide an improved accuracy wavelength solution, an optimised integrated polarimeter and the instrument integrated into a robotic operation. The robotic operation will optimise our programme which requires our target stars to be measured on a nightly basis. We present an overview of the entire project, including a description of our anticipated robotic operation. △ Less

Submitted 16 August, 2016; originally announced August 2016.

Comments: 13 pages, 8 figures, SPIE conference proceedings

Journal ref: Proc. SPIE 9908, Ground-based and Airborne Instrumentation for Astronomy VI, 99086F (August 9, 2016)

V346 Normae: First post-outburst observations of an FU Orionis star

Authors: Stefan Kraus, Alessio Caratti o Garatti, Rebeca Garcia-Lopez, Alexander Kreplin, Alicia Aarnio, John D. Monnier, Tim Naylor, Gerd Weigelt

Abstract: During their formation phase stars gain most of their mass in violent episodic accretion events, such as observed in FU Orionis (FUor) and EXor stars. V346 Normae is a well-studied FUor that underwent a strong outburst beginning in $\sim1980$. Here, we report photometric and spectroscopic observations which show that the visual/near-infrared brightness has decreased dramatically between the 1990s… ▽ More During their formation phase stars gain most of their mass in violent episodic accretion events, such as observed in FU Orionis (FUor) and EXor stars. V346 Normae is a well-studied FUor that underwent a strong outburst beginning in $\sim1980$. Here, we report photometric and spectroscopic observations which show that the visual/near-infrared brightness has decreased dramatically between the 1990s and 2010 ($ΔR\approx10.9^{\rm m}$, $ΔJ\approx7.8^{\rm m}$, $ΔK\approx5.8^{\rm m}$). The spectral properties of this fading event cannot be explained with variable extinction alone, but indicate a drop in accretion rate by 2-3 orders of magnitude, marking the first time that a member of the FUor class has been observed to switch to a very low accretion phase. Remarkably, in the last few years (2011-2015) V346 Nor has brightened again at all near-infrared wavelengths, indicating the onset of a new outburst event. The observed behaviour might be consistent with the clustered luminosity bursts that have been predicted by recent gravitational instability and fragmentation models for the early stages of protostellar evolution. Given V346 Nor's unique characteristics (concerning outburst duration, repetition frequency, and spectroscopic diagnostics), our results also highlight the need for revisiting the FUor/EXor classification scheme. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Comments: 5 pages, 2 figures, accepted at MNRAS Letters

Photoevaporation and close encounters: how the environment around Cygnus OB2 affects the evolution of protoplanetary disks

Authors: M. G. Guarcello, J. J. Drake, N. J. Wright, J. F. Albacete-Colombo, C. Clarke, B. Ercolano, E. Flaccomio, V. Kashyap, G. Micela, T. Naylor, N. Schneider, S. Sciortino, J. S. Vink

Abstract: In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedbac… ▽ More In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy. △ Less

Submitted 5 May, 2016; originally announced May 2016.

Comments: Accepted for publication in ApJS as part of the special issue on the Chandra Cygnus OB2 Legacy Project

Rapid Circumstellar Disk Evolution and an Accelerating Star Formation Rate in the Infrared Dark Cloud M17 SWex

Authors: Matthew S. Povich, Leisa K. Townsley, Thomas P. Robitaille, Patrick S. Broos, Wesley T. Orbin, Robert R. King, Tim Naylor, Barbara A. Whitney

Abstract: We present a catalog of 840 X-ray sources and first results from a 100 ks Chandra X-ray Observatory imaging study of the filamentary infrared dark cloud G014.225$-$00.506, which forms the central regions of a larger cloud complex known as the M17 southwest extension (M17 SWex). In addition to the rich population of protostars and young stellar objects with dusty circumstellar disks revealed by Spi… ▽ More We present a catalog of 840 X-ray sources and first results from a 100 ks Chandra X-ray Observatory imaging study of the filamentary infrared dark cloud G014.225$-$00.506, which forms the central regions of a larger cloud complex known as the M17 southwest extension (M17 SWex). In addition to the rich population of protostars and young stellar objects with dusty circumstellar disks revealed by Spitzer Space Telescope archival data, we discover a population of X-ray-emitting, intermediate-mass pre--main-sequence stars (IMPS) that lack infrared excess emission from circumstellar disks. We model the infrared spectral energy distributions of this source population to measure its mass function and place new constraints on the inner dust disk destruction timescales for 2-8 $M_{\odot}$ stars. We also place a lower limit on the star formation rate (SFR) and find that it is quite high ($\dot{M}\ge 0.007~M_{\odot}$ yr$^{-1}$), equivalent to several Orion Nebula Clusters in G14.225$-$0.506 alone, and likely accelerating. The cloud complex has not produced a population of massive, O-type stars commensurate with its SFR. This absence of very massive (${\ge}20~M_{\odot}$) stars suggests that either (1) M17 SWex is an example of a distributed mode of star formation that will produce a large OB association dominated by intermediate-mass stars but relatively few massive clusters, or (2) the massive cores are still in the process of accreting sufficient mass to form massive clusters hosting O stars. △ Less

Submitted 21 April, 2016; originally announced April 2016.

Comments: 29 pages, 9 figures, accepted to ApJ

The age of Taurus - environmental effects on disc lifetimes

Authors: Jon M. Rees, T. Wilson, C. P. M. Bell, R. D. Jeffries, Tim Naylor

Abstract: Using semi-empirical isochrones, we find the age of the Taurus star-forming region to be 3-4 Myr. Comparing the disc fraction in Taurus to young massive clusters suggests discs survive longer in this low density environment. We also present a method of photometrically de-reddening young stars using $iZJH$ data. Using semi-empirical isochrones, we find the age of the Taurus star-forming region to be 3-4 Myr. Comparing the disc fraction in Taurus to young massive clusters suggests discs survive longer in this low density environment. We also present a method of photometrically de-reddening young stars using $iZJH$ data. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: 2 pages, 1 figure, To appear in "Young Stars and Planets Near the Sun", Proceedings of IAU Symposium No. 314 (Cambridge University Press), J.H. Kastner, B. Stelzer, S.A. Metchev, eds

A self-consistent, absolute isochronal age scale for young moving groups in the solar neighbourhood

Authors: Cameron P. M. Bell, Eric E. Mamajek, Tim Naylor

Abstract: We present a self-consistent, absolute isochronal age scale for young (< 200 Myr), nearby (< 100 pc) moving groups in the solar neighbourhood based on homogeneous fitting of semi-empirical pre-main-sequence model isochrones using the tau^2 maximum-likelihood fitting statistic of Naylor & Jeffries in the M_V, V-J colour-magnitude diagram. The final adopted ages for the groups are: 149+51-19 Myr for… ▽ More We present a self-consistent, absolute isochronal age scale for young (< 200 Myr), nearby (< 100 pc) moving groups in the solar neighbourhood based on homogeneous fitting of semi-empirical pre-main-sequence model isochrones using the tau^2 maximum-likelihood fitting statistic of Naylor & Jeffries in the M_V, V-J colour-magnitude diagram. The final adopted ages for the groups are: 149+51-19 Myr for the AB Dor moving group, 24+/-3 Myr for the β Pic moving group (BPMG), 45+11-7 Myr for the Carina association, 42+6-4 Myr for the Columba association, 11+/-3 Myr for the η Cha cluster, 45+/-4 Myr for the Tucana-Horologium moving group (Tuc-Hor), 10+/-3 Myr for the TW Hya association, and 22+4-3 Myr for the 32 Ori group. At this stage we are uncomfortable assigning a final, unambiguous age to the Argus association as our membership list for the association appears to suffer from a high level of contamination, and therefore it remains unclear whether these stars represent a single population of coeval stars. Our isochronal ages for both the BPMG and Tuc-Hor are consistent with recent lithium depletion boundary (LDB) ages, which unlike isochronal ages, are relatively insensitive to the choice of low-mass evolutionary models. This consistency between the isochronal and LDB ages instills confidence that our self-consistent, absolute age scale for young, nearby moving groups is robust, and hence we suggest that these ages be adopted for future studies of these groups. Software implementing the methods described in this study is available from http: //www.astro.ex.ac.uk/people/timn/tau-squared/. △ Less

Submitted 24 August, 2015; originally announced August 2015.

Comments: Accepted for publication in MNRAS, 24 pages, 12 figures and 2 tables

Isochronal age scale of young moving groups in the solar neighbourhood

Authors: Cameron P. M. Bell, Eric E. Mamajek, Tim Naylor

Abstract: We present a self-consistent, absolute isochronal age scale for young (< 200 Myr), nearby (< 100 pc) moving groups, which is consistent with recent lithium depletion boundary ages for both the Beta Pic and Tucana-Horologium moving groups. This age scale was derived using a set of semi-empirical pre-main-sequence model isochrones that incorporate an empirical colour-Teff relation and bolometric cor… ▽ More We present a self-consistent, absolute isochronal age scale for young (< 200 Myr), nearby (< 100 pc) moving groups, which is consistent with recent lithium depletion boundary ages for both the Beta Pic and Tucana-Horologium moving groups. This age scale was derived using a set of semi-empirical pre-main-sequence model isochrones that incorporate an empirical colour-Teff relation and bolometric corrections based on the observed colours of Pleiades members, with theoretical corrections for the dependence on log g. Absolute ages for young, nearby groups are vital as these regions play a crucial role in our understanding of the early evolution of low- and intermediate-mass stars, as well as providing ideal targets for direct imaging and other measurements of dusty debris discs, substellar objects and, of course, extrasolar planets. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: 8 pages, Proceedings IAU Symposium No. 314, Young Stars & Planets Near the Sun, 2015, J. H. Kastner, B. Stelzer, & S. A. Metchev, eds

Optical and infrared counterparts of the X-ray sources detected in the Chandra Cygnus OB2 Legacy Survey

Authors: M. G. Guarcello, J. J. Drake, N. J. Wright, T. Naylor, E. Flaccomio, V. L. Kashyap, D. Garcia-Alvarez

Abstract: The young massive OB association Cygnus OB2, in the Cygnus X complex, is the closest (1400 pc) star forming region to the Sun hosting thousands of young low mass stars and up to 1000 OB stars, among which are some of the most massive stars known in our Galaxy. This region holds great importance for several fields of modern astrophysics, such as the study of the physical properties of massive and y… ▽ More The young massive OB association Cygnus OB2, in the Cygnus X complex, is the closest (1400 pc) star forming region to the Sun hosting thousands of young low mass stars and up to 1000 OB stars, among which are some of the most massive stars known in our Galaxy. This region holds great importance for several fields of modern astrophysics, such as the study of the physical properties of massive and young low-mass stars and the feedback provided by massive stars on star and planet formation process. Cygnus OB2 has been recently observed with Chandra/ACIS-I as part of the 1.08Msec Chandra Cygnus OB2 Legacy Project. This survey detected 7924 X-ray sources in a square degree area centered on Cyg OB2. Since a proper classification and study of the observed X-ray sources also requires the analysis of their optical and infrared counterparts, we combined a large and deep set of optical and infrared catalogs available for this region with our new X-ray catalog. In this paper we describe the matching procedure and present the combined catalog containing 5703 sources. We also briefly discuss the nature of the X-ray sources with optical and infrared counterparts using their position in the color-magnitude and color-color diagrams. △ Less

Submitted 15 January, 2015; originally announced January 2015.

Comments: Matched catalog can be provided by the authors; Accepted for publication by ApJ

Pre-main-sequence isochrones -- III. The Cluster Collaboration isochrone server

Authors: Cameron P. M. Bell, Jon M. Rees, Tim Naylor, N. J. Mayne, R. D. Jeffries, Eric E. Mamajek, John Rowe

Abstract: We present an isochrone server for semi-empirical pre-main-sequence model isochrones in the following systems: Johnson-Cousins, Sloan Digital Sky Survey, Two-Micron All-Sky Survey, Isaac Newton Telescope (INT) Wide-Field Camera, and INT Photometric H$α$ Survey (IPHAS)/UV-Excess Survey (UVEX). The server can be accessed via the Cluster Collaboration webpage {http://www.astro.ex.ac.uk/people/timn/is… ▽ More We present an isochrone server for semi-empirical pre-main-sequence model isochrones in the following systems: Johnson-Cousins, Sloan Digital Sky Survey, Two-Micron All-Sky Survey, Isaac Newton Telescope (INT) Wide-Field Camera, and INT Photometric H$α$ Survey (IPHAS)/UV-Excess Survey (UVEX). The server can be accessed via the Cluster Collaboration webpage {http://www.astro.ex.ac.uk/people/timn/isochrones/}. To achieve this we have used the observed colours of member stars in young clusters with well-established age, distance and reddening to create fiducial loci in the colour-magnitude diagram. These empirical sequences have been used to quantify the discrepancy between the models and data arising from uncertainties in both the interior and atmospheric models, resulting in tables of semi-empirical bolometric corrections (BCs) in the various photometric systems. The model isochrones made available through the server are based on existing stellar interior models coupled with our newly derived semi-empirical BCs. As part of this analysis we also present new cluster parameters for both the Pleiades and Praesepe, yielding ages of $135^{+20}_{-11}$ and $665^{+14}_{-7}\,\rm{Myr}$ as well as distances of $132 \pm 2$ and $184 \pm 2\,\rm{pc}$ respectively (statistical uncertainty only). △ Less

Submitted 18 September, 2014; originally announced September 2014.

Comments: Accepted for publication in MNRAS; 18 pages, 8 figures and 7 tables. Tables 3, 4, 6 and 7 (for which only samples are given in the paper) are available in full at the Cluster Collaboration homepage {http://www.astro.ex.ac.uk/people/timn/Catalogues/paper3.html}

The Second Data Release of the INT Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS DR2)

Authors: Geert Barentsen, H. J. Farnhill, J. E. Drew, E. A. González-Solares, R. Greimel, M. J. Irwin, B. Miszalski, C. Ruhland, P. Groot, A. Mampaso, S. E. Sale, A. A. Henden, A. Aungwerojwit, M. J. Barlow, P. J. Carter, R. L. M. Corradi, J. J. Drake, J. Eislöffel, J. Fabregat, B. T. Gänsicke, N. P. Gentile Fusillo, S. Greiss, A. S. Hales, S. Hodgkin, L. Huckvale , et al. (28 additional authors not shown)

Abstract: The INT/WFC Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS) is a 1800 square degrees imaging survey covering Galactic latitudes |b| < 5 deg and longitudes l = 30 to 215 deg in the r, i and H-alpha filters using the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) in La Palma. We present the first quality-controlled and globally-calibrated source catalogue der… ▽ More The INT/WFC Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS) is a 1800 square degrees imaging survey covering Galactic latitudes |b| < 5 deg and longitudes l = 30 to 215 deg in the r, i and H-alpha filters using the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) in La Palma. We present the first quality-controlled and globally-calibrated source catalogue derived from the survey, providing single-epoch photometry for 219 million unique sources across 92% of the footprint. The observations were carried out between 2003 and 2012 at a median seeing of 1.1 arcsec (sampled at 0.33 arcsec/pixel) and to a mean 5σ-depth of 21.2 (r), 20.0 (i) and 20.3 (H-alpha) in the Vega magnitude system. We explain the data reduction and quality control procedures, describe and test the global re-calibration, and detail the construction of the new catalogue. We show that the new calibration is accurate to 0.03 mag (rms) and recommend a series of quality criteria to select the most reliable data from the catalogue. Finally, we demonstrate the ability of the catalogue's unique (r-Halpha, r-i) diagram to (1) characterise stellar populations and extinction regimes towards different Galactic sightlines and (2) select H-alpha emission-line objects. IPHAS is the first survey to offer comprehensive CCD photometry of point sources across the Galactic Plane at visible wavelengths, providing the much-needed counterpart to recent infrared surveys. △ Less

Submitted 12 August, 2014; v1 submitted 18 June, 2014; originally announced June 2014.

Comments: Accepted for publication in MNRAS. The catalogue is available in full from the survey website at http://www.iphas.org and has been submitted to Vizier

The Spatial Structure of Young Stellar Clusters. I. Subclusters

Authors: Michael A. Kuhn, Eric D. Feigelson, Konstantin V. Getman, Adrian J. Baddeley, Patrick S. Broos, Alison Sills, Matthew R. Bate, Matthew S. Povich, Kevin L. Luhman, Heather A. Busk, Tim Naylor, Robert R. King

Abstract: The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of… ▽ More The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid. △ Less

Submitted 17 March, 2014; originally announced March 2014.

Comments: Accepted for publication in The Astrophysical Journal; 67 pages, 9 figures, and 4 tables. Supplemental materials from the online-version of this article, including machine-readable tables and source code in R are available at the public MYStIX website (http://astro.psu.edu/mystix/)

Age Gradients in the Stellar Populations of Massive Star Forming Regions Based on a New Stellar Chronometer

Authors: Konstantin V. Getman, Eric D. Feigelson, Michael A. Kuhn, Patrick S. Broos, Leisa K. Townsley, Tim Naylor, Matthew S. Povich, Kevin L. Luhman, Gordon P. Garmire

Abstract: A major impediment to understanding star formation in massive star forming regions (MSFRs) is the absence of a reliable stellar chronometer to unravel their complex star formation histories. We present a new estimation of stellar ages using a new method that employs near-infrared (NIR) and X-ray photometry, AgeJX. Stellar masses are derived from X-ray luminosities using the Lx - Mass relation from… ▽ More A major impediment to understanding star formation in massive star forming regions (MSFRs) is the absence of a reliable stellar chronometer to unravel their complex star formation histories. We present a new estimation of stellar ages using a new method that employs near-infrared (NIR) and X-ray photometry, AgeJX. Stellar masses are derived from X-ray luminosities using the Lx - Mass relation from the Taurus cloud. J-band luminosities are compared to mass-dependent pre-main-sequence evolutionary models to estimate ages. AgeJX is sensitive to a wide range of evolutionary stages, from disk-bearing stars embedded in a cloud to widely dispersed older pre-main sequence stars. The MYStIX (Massive Young Star-Forming Complex Study in Infrared and X-ray) project characterizes 20 OB-dominated MSFRs using X-ray, mid-infrared, and NIR catalogs. The AgeJX method has been applied to 5525 out of 31,784 MYStIX Probable Complex Members. We provide a homogeneous set of median ages for over a hundred subclusters in 15 MSFRs; median subcluster ages range between 0.5 Myr and 5 Myr. The important science result is the discovery of age gradients across MYStIX regions. The wide MSFR age distribution appears as spatially segregated structures with different ages. The AgeJX ages are youngest in obscured locations in molecular clouds, intermediate in revealed stellar clusters, and oldest in distributed populations. The NIR color index J-H, a surrogate measure of extinction, can serve as an approximate age predictor for young embedded clusters. △ Less

Submitted 11 March, 2014; originally announced March 2014.

Comments: Accepted for publication in ApJ; 89 pages, 23 figures, 2 Tables; High quality version is at http://astro.psu.edu/mystix

The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+)

Authors: J. E. Drew, E. Gonzalez-Solares, R. Greimel, M. J. Irwin, A. Kupcu Yoldas, J. Lewis, G. Barentsen, J. Eisloeffel, H. J. Farnhill, W. E. Martin, J. R. Walsh, N. A. Walton, M. Mohr-Smith, R. Raddi, S. E. Sale, N. J. Wright, P. Groot, M. J. Barlow, R. L. M. Corradi, J. J. Drake, J. Fabregat, D. J. Frew, B. T. Gaensicke, C. Knigge, A. Mampaso , et al. (10 additional authors not shown)

Abstract: The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+) is surveying the southern Milky Way in u, g, r, i and Halpha at 1 arcsec angular resolution. Its footprint spans the Galactic latitude range -5 < b < +5 at all longitudes south of the celestial equator. Extensions around the Galactic Centre to Galactic latitudes +/-10 bring in much of the Galactic Bulge. This ESO p… ▽ More The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+) is surveying the southern Milky Way in u, g, r, i and Halpha at 1 arcsec angular resolution. Its footprint spans the Galactic latitude range -5 < b < +5 at all longitudes south of the celestial equator. Extensions around the Galactic Centre to Galactic latitudes +/-10 bring in much of the Galactic Bulge. This ESO public survey, begun on 28th December 2011, reaches down to 20th magnitude (10-sigma) and will provide single-epoch digital optical photometry for around 300 million stars. The observing strategy and data pipelining is described, and an appraisal of the segmented narrowband Halpha filter in use is presented. Using model atmospheres and library spectra, we compute main-sequence (u - g), (g - r), (r - i) and (r - Halpha) stellar colours in the Vega system. We report on a preliminary validation of the photometry using test data obtained from two pointings overlapping the Sloan Digital Sky Survey. An example of the (u - g, g - r) and (r - Halpha, r - i) diagrams for a full VPHAS+ survey field is given. Attention is drawn to the opportunities for studies of compact nebulae and nebular morphologies that arise from the image quality being achieved. The value of the u band as the means to identify planetary-nebula central stars is demonstrated by the discovery of the central star of NGC 2899 in survey data. Thanks to its excellent imaging performance, the VST/OmegaCam combination used by this survey is a perfect vehicle for automated searches for reddened early-type stars, and will allow the discovery and analysis of compact binaries, white dwarfs and transient sources. △ Less

Submitted 2 March, 2014; v1 submitted 27 February, 2014; originally announced February 2014.

Comments: 24 pages, 23 figures (supplementary tables not included, see www.vphas.eu for publisher's version). Accepted for publication by MNRAS, 24th February 2014. Replacing original arxiv upload to correct corrupt figures

An Anomalous Extinction Law in the Cep OB3b Young Cluster: Evidence for dust processing during gas dispersal

Authors: Thomas S. Allen, Jakub J. Prchlik, S. Thomas Megeath, Robert A. Gutermuth, Judith L. Pipher, Tim Naylor, Rob D. Jeffries

Abstract: We determine the extinction law through Cep OB3b, a young cluster of 3000 stars undergoing gas dispersal. The extinction is measured toward 76 background K giants identified with MMT/Hectospec spectra. Color excess ratios were determined toward each of the giants using $V$ and $R$ photometry from the literature, $g$,$r$,$i$ and $z$ photometry from SDSS and $J$, $H$, and $K_{s}$ photometry from 2MA… ▽ More We determine the extinction law through Cep OB3b, a young cluster of 3000 stars undergoing gas dispersal. The extinction is measured toward 76 background K giants identified with MMT/Hectospec spectra. Color excess ratios were determined toward each of the giants using $V$ and $R$ photometry from the literature, $g$,$r$,$i$ and $z$ photometry from SDSS and $J$, $H$, and $K_{s}$ photometry from 2MASS. These color excess ratios were the used to construct the extinction law through the dusty material associated with Cep OB3b. The extinction law through Cep OB3b is intermediate between the $R_{V} = 3.1$ and $R_{V} = 5$ laws commonly used for the diffuse atomic ISM and dense molecular clouds, respectively. The dependence of the extinction law on line-of-sight $A_{V}$ is investigated and we find the extinction law becomes shallower for regions with $A_{V} > 2.5$ magnitudes. We speculate that the intermediate dust law results from dust processing during the dispersal of the molecular cloud by the cluster. △ Less

Submitted 6 February, 2014; originally announced February 2014.

Comments: 31 pages, 10 Figures, 3 Tables, accepted for publication in ApJ

Ages of young stars

Authors: David R. Soderblom, Lynne A. Hillenbrand, Rob D. Jeffries, Eric E. Mamajek, Tim Naylor

Abstract: Determining the sequence of events in the formation of stars and planetary systems and their time-scales is essential for understanding those processes, yet establishing ages is fundamentally difficult because we lack direct indicators. In this review we discuss the age challenge for young stars, specifically those less than ~100 Myr old. Most age determination methods that we discuss are primaril… ▽ More Determining the sequence of events in the formation of stars and planetary systems and their time-scales is essential for understanding those processes, yet establishing ages is fundamentally difficult because we lack direct indicators. In this review we discuss the age challenge for young stars, specifically those less than ~100 Myr old. Most age determination methods that we discuss are primarily applicable to groups of stars but can be used to estimate the age of individual objects. A reliable age scale is established above 20 Myr from measurement of the Lithium Depletion Boundary (LDB) in young clusters, and consistency is shown between these ages and those from the upper main sequence and the main sequence turn-off -- if modest core convection and rotation is included in the models of higher-mass stars. Other available methods for age estimation include the kinematics of young groups, placing stars in Hertzsprung-Russell diagrams, pulsations and seismology, surface gravity measurement, rotation and activity, and lithium abundance. We review each of these methods and present known strengths and weaknesses. Below ~20 Myr, both model-dependent and observational uncertainties grow, the situation is confused by the possibility of age spreads, and no reliable absolute ages yet exist. The lack of absolute age calibration below 20 Myr should be borne in mind when considering the lifetimes of protostellar phases and circumstellar material. △ Less

Submitted 27 November, 2013; originally announced November 2013.

Comments: Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning

The dependence of stellar age distributions on GMC environment

Authors: C. L. Dobbs, J. E. Pringle, T. Naylor

Abstract: In this Letter, we analyse the distributions of stellar ages in Giant Molecular Clouds (GMCs) in spiral arms, inter-arm spurs, and at large galactic radii, where the spiral arms are relatively weak. We use the results of numerical simulations of galaxies, which follow the evolution of GMCs and include star particles where star formation events occur. We find that GMCs in spiral arms tend to have p… ▽ More In this Letter, we analyse the distributions of stellar ages in Giant Molecular Clouds (GMCs) in spiral arms, inter-arm spurs, and at large galactic radii, where the spiral arms are relatively weak. We use the results of numerical simulations of galaxies, which follow the evolution of GMCs and include star particles where star formation events occur. We find that GMCs in spiral arms tend to have predominantly young (< 10 Myr) stars. By contrast, clouds which are the remainders of spiral arm GMAs that have been sheared into inter-arm GMCs, contain fewer young (< 10 Myr) stars, and more ~20 Myr stars. We also show that clouds which form in the absence of spiral arms, due to local gravitational and thermal instabilities, contain preferentially young stars. We propose the age distributions of stars in GMCs will be a useful diagnostic to test different cloud evolution scenarios, the origin of spiral arms, and the success of numerical models of galactic star formation. We discuss the implications of our results in the context of Galactic and extragalactic molecular clouds. △ Less

Submitted 24 September, 2013; originally announced September 2013.

Comments: 5 pages, 2 figures, accepted for publication in MNRAS Letters

Identifying Young Stars in Massive Star-Forming Regions for the MYStIX Project

Authors: Patrick S. Broos, Konstantin V. Getman, Matthew S. Povich, Eric D. Feigelson, Leisa K. Townsley, Tim Naylor, Michael A. Kuhn, R. R. King, Heather A. Busk

Abstract: The Massive Young star-forming Complex Study in Infrared and X-rays (MYStIX) project requires samples of young stars that are likely members of 20 nearby Galactic massive star-forming regions. Membership is inferred from statistical classification of X-ray sources, from detection of a robust infrared excess that is best explained by circumstellar dust in a disk or infalling envelope, and from publ… ▽ More The Massive Young star-forming Complex Study in Infrared and X-rays (MYStIX) project requires samples of young stars that are likely members of 20 nearby Galactic massive star-forming regions. Membership is inferred from statistical classification of X-ray sources, from detection of a robust infrared excess that is best explained by circumstellar dust in a disk or infalling envelope, and from published spectral types that are unlikely to be found among field stars. We present the MYStIX membership lists here, and describe in detail the statistical classification of X-ray sources via a \Naive Bayes Classi △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 47 pages, 7 figures. Accepted for publication of Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of all MYStIX papers are available at http://www.astro.psu.edu/mystix

The MYStIX InfraRed-Excess Source Catalog

Authors: Matthew S. Povich, Michael A. Kuhn, Konstantin V. Getman, Heather A. Busk, Eric D. Feigelson, Patrick S. Broos, Leisa K. Townsley, Robert R. King, Tim Naylor

Abstract: The MYStIX project (Massive Young Star-Forming Complex Study in Infrared and X-rays) provides a comparative study of 20 Galactic massive star-forming complexes (d = 0.4 to 3.6 kpc). Probable stellar members in each target complex are identified using X-ray and/or infrared data via two pathways: (1) X-ray detections of young/massive stars with coronal activity/strong winds; or (2) infrared excess (… ▽ More The MYStIX project (Massive Young Star-Forming Complex Study in Infrared and X-rays) provides a comparative study of 20 Galactic massive star-forming complexes (d = 0.4 to 3.6 kpc). Probable stellar members in each target complex are identified using X-ray and/or infrared data via two pathways: (1) X-ray detections of young/massive stars with coronal activity/strong winds; or (2) infrared excess (IRE) selection of young stellar objects (YSOs) with circumstellar disks and/or protostellar envelopes. We present the methodology for the second pathway, using Spitzer/IRAC, 2MASS, and UKIRT imaging and photometry. Although IRE selection of YSOs is welltrodden territory, MYStIX presents unique challenges. We combine IR spectral energy distribution (SED) fitting with IR color cuts and spatial clustering analysis to identify IRE sources and isolate probable YSO members in each MYStIX target field from the myriad types of contaminating sources that can resemble YSOs: extragalactic sources, evolved stars, nebular knots, and even unassociated foreground/background YSOs. Applying our methodology consistently across 18 of the target complexes, we produce the MYStIX IRE Source (MIRES) Catalog comprising 20,719 sources, including 8686 probable stellar members of the MYStIX target complexes. We also classify the SEDs of 9365 IR counterparts to MYStIX X-ray sources to assist the first pathway, the identification of X-ray detected stellar members. The MIRES catalog provides a foundation for follow-up studies of diverse phenomena related to massive star cluster formation, including protostellar outflows, circumstellar disks, and sequential star formation triggered by massive star feedback processes. △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 67 pages, 11 figures. Accepted for publication in Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of all MYStIX papers are available at http://www.astro.psu.edu/mystix

Bayesian Matching for X-ray and Infrared Sources in the MYStIX Project

Authors: Tim Naylor, Patrick S. Broos, Eric D. Feigelson

Abstract: Identifying the infrared counterparts of X-ray sources in Galactic Plane fields such as those of the MYStIX project presents particular difficulties due to the high density of infrared sources. This high stellar density makes it inevitable that a large fraction of X-ray positions will have a faint field star close to them, which standard matching techniques may incorrectly take to be the counterpa… ▽ More Identifying the infrared counterparts of X-ray sources in Galactic Plane fields such as those of the MYStIX project presents particular difficulties due to the high density of infrared sources. This high stellar density makes it inevitable that a large fraction of X-ray positions will have a faint field star close to them, which standard matching techniques may incorrectly take to be the counterpart. Instead we use the infrared data to create a model of both the field star and counterpart magnitude distributions, which we then combine with a Bayesian technique to yield a probability that any star is the counterpart of an X-ray source. In our more crowded fields, between 10 and 20% of counterparts that would be identified on the grounds of being the closest star to X-ray position within a 99% confidence error circle are instead identified by the Bayesian technique as field stars. These stars are preferentially concentrated at faint magnitudes. Equally importantly the technique also gives a probability that the true counterpart to the X-ray source falls beneath the magnitude limit of the infrared catalog. In deriving our method, we place it in the context of other procedures for matching astronomical catalogs. △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 12 pages, 6 figures. Accepted for publication in Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of all MYStIX papers are available at http://www.astro.psu.edu/mystix

The Wide-Field Near Infrared Data: Optimal Photometry in Crowded Fields

Authors: R. R. King, Tim Naylor, Patrick S. Broos, Konstantin V. Getman, Eric D. Feigelson

Abstract: We present JHK infrared data from the UK Infrared Telescope for a subset of the regions of the MYStIX (Massive Young Star-Forming Complex Study in Infrared and X-ray) survey. Some of the data were obtained specifically for the MYStIX project, and some as part of the UKIRT Infrared Deep Sky Survey's Galactic Plane Survey. In most of these fields crowding is a significant issue for aperture photomet… ▽ More We present JHK infrared data from the UK Infrared Telescope for a subset of the regions of the MYStIX (Massive Young Star-Forming Complex Study in Infrared and X-ray) survey. Some of the data were obtained specifically for the MYStIX project, and some as part of the UKIRT Infrared Deep Sky Survey's Galactic Plane Survey. In most of these fields crowding is a significant issue for aperture photometry, and so we have re-extracted the photometry from the processed images using an optimal extraction technique, and we describe how we adapt the optimal technique to mitigate the effects of crowding. △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 13 pages, 5 figures. Accepted for publication in Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of this and other MYStIX papers are available at http://www.astro.psu.edu/mystix

Overview of the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Project

Authors: Eric D. Feigelson, Leisa K. Townsley, Patrick S. Broos, Heather A. Busk, Konstantin V. Getman, Robert King, Michael A. Kuhn, Tim Naylor, Matthew Povich, Adrian Baddeley, Matthew Bate, Remy Indebetouw, Kevin Luhman, Mark McCaughrean, Julian Pittard, Ralph Pudritz, Alison Sills, Yong Song, James Wadsley

Abstract: MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined… ▽ More MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula. △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 74 pages, 10 figures. To appear in Astrophysical Journal Supplements (with 6 other MYStIX papers). A version with complete figures is available at http://www.astro.psu.edu/mystix

A lithium depletion boundary age of 22 Myr for NGC 1960

Authors: R. D. Jeffries, Tim Naylor, N. J. Mayne, Cameron P. M. Bell, S. P. Littlefair

Abstract: We present a deep Cousins RI photometric survey of the open cluster NGC 1960, complete to R_C \simeq 22, I_C \simeq 21, that is used to select a sample of very low-mass cluster candidates. Gemini spectroscopy of a subset of these is used to confirm membership and locate the age-dependent "lithium depletion boundary" (LDB) --the luminosity at which lithium remains unburned in its low-mass stars. Th… ▽ More We present a deep Cousins RI photometric survey of the open cluster NGC 1960, complete to R_C \simeq 22, I_C \simeq 21, that is used to select a sample of very low-mass cluster candidates. Gemini spectroscopy of a subset of these is used to confirm membership and locate the age-dependent "lithium depletion boundary" (LDB) --the luminosity at which lithium remains unburned in its low-mass stars. The LDB implies a cluster age of 22 +/-4 Myr and is quite insensitive to choice of evolutionary model. NGC 1960 is the youngest cluster for which a LDB age has been estimated and possesses a well populated upper main sequence and a rich low-mass pre-main sequence. The LDB age determined here agrees well with precise age estimates made for the same cluster based on isochrone fits to its high- and low-mass populations. The concordance between these three age estimation techniques, that rely on different facets of stellar astrophysics at very different masses, is an important step towards calibrating the absolute ages of young open clusters and lends confidence to ages determined using any one of them. △ Less

Submitted 26 June, 2013; originally announced June 2013.

Comments: Accepted for publication in MNRAS