Source author record

Tim Naylor

Tim Naylor appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2025arXiv

Using Stellar Spectral Energy Distributions to Measure Exoplanet Parameters

The ability to make accurate determinations of planetary parameters is inextricably linked to measuring physical parameters of the host star, in particular the stellar radius. In this paper we fit the stellar spectral energy distributions of exoplanet hosts to measure their radii, making use of only archival photometry, the $Gaia$ parallaxes and $Gaia$ extinction maps. Using the extinction maps frees us of the degeneracy between temperature and extinction which has plagued this method in the past. The resulting radii have typical random uncertainties of about 2 per cent. We perform a quantitative study of systematic uncertainties affecting the methodology and find they are similar to, or smaller than, the random ones. We discuss how the stellar parameters can be used to derive the properties of both transiting exoplanets, and those where only a radial-velocity curve is available. We then explore in detail the improvements the method makes possible for the parameters of the PanCET sample of transiting planets. For this sample we find the best literature measurements of the planetary radii have mean uncertainties about 40 per cent larger than those presented here, with the new measurements achieving precisions of 2 per cent in radius and 10 per cent in mass. In contrast to much recent work, these transiting exoplanets parameters are derived without using theoretical models of stellar interiors, freeing them of the assumptions those models contain, and any priors for stellar age. As the data used are available for the whole sky, the method can be used for self-consistent measurements of the planetary parameters of a very large fraction of known exoplanets.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2019arXiv

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

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.

preprint2016arXiv

HARPS3 for a Roboticized Isaac Newton Telescope

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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/.

preprint2015arXiv

Isochronal age scale of young moving groups in the solar neighbourhood

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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).

preprint2014arXiv

The Spatial Structure of Young Stellar Clusters. I. Subclusters

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.

preprint2013arXiv

A lithium depletion boundary age of 22 Myr for NGC 1960

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.

preprint2013arXiv

Ages of young stars

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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

preprint2013arXiv

No evidence for intense, cold accretion onto YSOs from measurements of Li in T-Tauri stars

We have used medium resolution spectra to search for evidence that proto-stellar objects accrete at high rates during their early 'assembly phase'. Models predict that depleted lithium and reduced luminosity in T-Tauri stars are key signatures of 'cold' high-rate accretion occurring early in a star's evolution. We found no evidence in 168 stars in NGC 2264 and the Orion Nebula Cluster for strong lithium depletion through analysis of veiling corrected 6708 angstrom lithium spectral line strengths. This suggests that 'cold' accretion at high rates (M_dot > 5 x 10-4 M_sol yr-1) occurs in the assembly phase of fewer than 0.5 per cent of 0.3 < M < 1.9 M_sol stars. We also find that the dispersion in the strength of the 6708 angstrom lithium line might imply an age spread that is similar in magnitude to the apparent age spread implied by the luminosity dispersion seen in colour magnitude diagrams. Evidence for weak lithium depletion (< 10 per cent in equivalent width) that is correlated with luminosity is also apparent, but we are unable to determine whether age spreads or accretion at rates less than 5 x 10-4 M_sol yr-1 are responsible.

preprint2013arXiv

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

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.

preprint2013arXiv

Pre-main-sequence isochrones -- II. Revising star and planet formation timescales

We have derived ages for 13 young (<30 Myr) star-forming regions and find they are up to a factor two older than the ages typically adopted in the literature. This result has wide-ranging implications, including that circumstellar discs survive longer (~10-12 Myr) and that the average Class I lifetime is greater (~1 Myr) than currently believed. For each star-forming region we derived two ages from colour-magnitude diagrams. First we fitted models of the evolution between the zero-age main-sequence and terminal-age main-sequence to derive a homogeneous set of main-sequence ages, distances and reddenings with statistically meaningful uncertainties. Our second age for each star-forming region was derived by fitting pre-main-sequence stars to new semi-empirical model isochrones. For the first time (for a set of clusters younger than 50 Myr) we find broad agreement between these two ages, and since these are derived from two distinct mass regimes that rely on different aspects of stellar physics, it gives us confidence in the new age scale. This agreement is largely due to our adoption of empirical colour-Teff relations and bolometric corrections for pre-main-sequence stars cooler than 4000 K. The revised ages for the star-forming regions in our sample are: ~2 Myr for NGC 6611 (Eagle Nebula; M 16), IC 5146 (Cocoon Nebula), NGC 6530 (Lagoon Nebula; M 8), and NGC 2244 (Rosette Nebula); ~6 Myr for σ Ori, Cep OB3b, and IC 348; ~10 Myr for λ Ori (Collinder 69); ~11 Myr for NGC 2169; ~12 Myr for NGC 2362; ~13 Myr for NGC 7160; ~14 Myr for χ Per (NGC 884); and ~20 Myr for NGC 1960 (M 36).

preprint2013arXiv

The MYStIX InfraRed-Excess Source Catalog

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.

preprint2013arXiv

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

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.

preprint2012arXiv

MYStIX First Results: Spatial Structures of Massive Young Stellar Clusters

Observations of the spatial distributions of young stars in star-forming regions can be linked to the theory of clustered star formation using spatial statistical methods. The MYStIX project provides rich samples of young stars from the nearest high-mass star-forming regions. Maps of stellar surface density reveal diverse structure and subclustering. Young stellar clusters and subclusters are fit with isothermal spheres and ellipsoids using the Bayesian Information Criterion to estimate the number of subclusters. Clustering is also investigated using Cartwright and Whitworth's Q statistic and the inhomogeneous two-point correlation function. Mass segregation is detected in several cases, in both centrally concentrated and fractally structured star clusters, but a few clusters are not mass segregated.

preprint2012arXiv

Pre-main-sequence isochrones -- I. The Pleiades benchmark

We present a critical assessment of commonly used pre-main-sequence isochrones by comparing their predictions to a set of well-calibrated colour-magnitude diagrams of the Pleiades in the wavelength range 0.4 to 2.5 microns. Our analysis shows that for temperatures less than 4000 K the models systematically overestimate the flux by a factor two at 0.5 microns, though this decreases with wavelength, becoming negligible at 2.2 microns. In optical colours this will result in the ages for stars younger than 10 Myr being underestimated by factors between two and three. We show that using observations of standard stars to transform the data into a standard system can introduce significant errors in the positioning of pre-main-sequences in colour-magnitude diagrams. Therefore we have compared the models to the data in the natural photometric system in which the observations were taken. Thus we have constructed and tested a model of the system responses for the Wide-Field Camera on the Isaac Newton Telescope. As a benchmark test for the development of pre-main-sequence models we provide both our system responses and the Pleiades sequence.

preprint2012arXiv

Spitzer Imaging of the Nearby Rich Young Cluster, Cep OB3b

We map the full extent of a rich massive young cluster in the Cep OB3b association with the IRAC and MIPS instruments aboard the {\it Spitzer} Space Telescope and the ACIS instrument aboard the $\it{Chandra}$ X-Ray Observatory. At 700 pc, it is revealed to be the second nearest large ($>1000$ member), young ($< 5$ Myr) cluster known. In contrast to the nearest large cluster, the Orion Nebula Cluster, Cep OB3b is only lightly obscured and is mostly located in a large cavity carved out of the surrounding molecular cloud. Our infrared and X-ray datasets, as well as visible photometry from the literature, are used to take a census of the young stars in Cep OB3b. We find that the young stars within the cluster are concentrated in two sub-clusters; an eastern sub-cluster, near the Cep B molecular clump, and a western sub-cluster, near the Cep F molecular clump. Using our census of young stars, we examine the fraction of young stars with infrared excesses indicative of circumstellar disks. We create a map of the disk fraction throughout the cluster and find that it is spatially variable. Due to these spatial variations, the two sub-clusters exhibit substantially different average disk fractions from each other: $32% \pm 4%$ and $50% \pm 6%$. We discuss whether the discrepant disk fractions are due to the photodestruction of disks by the high mass members of the cluster or whether they result from differences in the ages of the sub-clusters. We conclude that the discrepant disk fractions are most likely due to differences in the ages.

preprint2011arXiv

Accretion-induced luminosity spreads in young clusters: evidence from stellar rotation

We present an analysis of the rotation of young stars in the associations Cepheus OB3b, NGC 2264, NGC 2362 and the Orion Nebula Cluster (ONC). We discover a correlation between rotation rate and position in a colour-magnitude diagram (CMD) such that stars which lie above an empirically determined median pre-main sequence rotate more rapidly than stars which lie below this sequence. The same correlation is seen, with a high degree of statistical significance, in each association studied here. If position within the CMD is interpreted as being due to genuine age spreads within a cluster, then the stars above the median pre-main sequence would be the youngest stars. This would in turn imply that the most rapidly rotating stars in an association are the youngest, and hence those with the largest moments of inertia and highest likelihood of ongoing accretion. Such a result does not fit naturally into the existing picture of angular momentum evolution in young stars, where the stars are braked effectively by their accretion discs until the disc disperses. Instead, we argue that, for a given association of young stars, position within the CMD is not primarily a function of age, but of accretion history. We show that this hypothesis could explain the correlation we observe between rotation rate and position within the CMD.

preprint2011arXiv

No wide spread of stellar ages in the Orion Nebula Cluster

The wide luminosity dispersion seen for stars at a given effective temperature in the H-R diagrams of young clusters and star forming regions is often interpreted as due to significant (~10 Myr) spreads in stellar contraction age. In the scenario where most stars are born with circumstellar discs, and that disc signatures decay monotonically (on average) over timescales of only a few Myr, then any such age spread should lead to clear differences in the age distributions of stars with and without discs. We have investigated large samples of stars in the Orion Nebula Cluster (ONC) using three methods to diagnose disc presence from infrared measurements. We find no significant difference in the mean ages or age distributions of stars with and without discs, consistent with expectations for a coeval population. Using a simple quantitative model we show that any real age spread must be smaller than the median disc lifetime. For a log-normal age distribution, there is an upper limit of <0.14 dex (at 99% confidence) to any real age dispersion, compared to the ~=0.4 dex implied by the H-R diagram. If the mean age of the ONC is 2.5 Myr, this would mean at least 95% of its low-mass stellar population has ages between 1.3--4.8 Myr. We suggest that the observed luminosity dispersion is caused by a combination of observational uncertainties and physical mechanisms that disorder the conventional relationship between luminosity and age for pre main-sequence stars. This means that individual stellar ages from the H-R diagram are unreliable and cannot be used to directly infer a star formation history. Irrespective of what causes the wide luminosity dispersion, the finding that any real age dispersion is less than the median disc lifetime argues strongly against star formation scenarios for the ONC lasting longer than a few Myr.

preprint2009arXiv

Are pre-main-sequence stars older than we thought?

We fit the colour-magnitude diagrams of stars between the zero-age main-sequence and terminal-age main sequence in young clusters and associations. The ages we derive are a factor 1.5 to 2 longer than the commonly used ages for these regions, which are derived from the positions of pre-main-sequence stars in colour-magnitude diagrams. From an examination of the uncertainties in the main-sequence and pre-main-sequence models, we conclude that the longer age scale is probably the correct one, which implies we must revise upwards the commonly used ages for young clusters and associations. Such a revision would explain the discrepancy between the observational lifetimes of proto-planetary discs and theoretical calculations of the time to form planets. It would also explain the absence of clusters with ages between 5 and 30Myr. We use the $τ^2$ statistic to fit the main-sequence data, but find that we must make significant modifications if we are to fit sequences which have vertical segments in the colour-magnitude diagram. We present this modification along with improvements to methods of calculating the goodness-of-fit statistic and parameter uncertainties. Software implementing the methods described in this paper is available from http://www.astro.ex.ac.uk/people/timn/tau-squared/

preprint2009arXiv

Are pre-MS stars older than we thought?

We present a consistent age ordering for young clusters and groups determined using the contraction of stars through their pre-main-sequence phase. We compare these with ages derived from the evolution of the upper main-sequence stars, and find the upper MS ages are older by a factor 1.5 to 2. We show that increasing the binary fraction and number of equal-mass binaries amongst the O-stars compared to the rest of the MS cannot remove this discrepancy.

preprint2009arXiv

Pre-main-sequence variability across the radiative-convective gap

We use I band imaging to perform a variability survey of the 13 Myr-old cluster h Per. We find a significant fraction of the cluster members to be variable. Most importantly, we find that variable members lie almost entirely on the convective side of the gap in the cluster sequence between fully convective stars and those which have a radiative core. This result is consistent with a scenario in which the magnetic field changes topology when the star changes from being fully convective, to one containing a radiative core. When the star is convective the magnetic field appears dominated by large-scale structures, resulting in global-size spots that drive the observed variability. For those stars with radiative cores we observe a marked absence of variability due to spots, which suggests a switch to a magnetic field dominated by smaller-scale structures, resulting in many smaller spots and thus less apparent variability. This implies that wide field variability surveys may only be sensitive to fully convective stars. On the one hand this reduces the chances of picking out young groups (since the convective stars are the lower mass and therefore fainter objects), but conversely the absolute magnitude of the head of the convective sequence provides a straightforward measure of age for those groups which are discovered.

preprint2009arXiv

Rotation of young stars in Cepheus OB3b

We present a photometric study of I-band variability in the young association Cepheus OB3b. The study is sensitive to periodic variability on timescales of less than a day, to more than 20 days. After rejection of contaminating objects using V, I, R and narrowband H-alpha photometry, we find 475 objects with measured rotation periods, which are very likely pre-main-sequence members of the Cep OB3b star forming region. We revise the distance and age to Cep OB3b, putting it on the self-consistent age and distance ladder of Mayne & Naylor (2008). This yields a distance modulus of 8.8 +/- 0.2 mags, corresponding to a distance of 580 +/- 60 pc, and an age of 4-5Myrs. The rotation period distribution confirms the general picture of rotational evolution in young stars, exhibiting both the correlation between accretion (determined in this case through narrowband H-alpha photometry) and rotation expected from disc locking, and the dependence of rotation upon mass that is seen in other star forming regions. However, this mass dependence is much weaker in our data than found in other studies. Comparison to the similarly aged NGC 2362 shows that the low-mass stars in Cep OB3b are rotating much more slowly. This points to a possible link between star forming environment and rotation properties. Such a link would call into question models of stellar angular momentum evolution, which assume that the rotational period distributions of young clusters and associations can be assembled into an evolutionary sequence, thus ignoring environmental effects.

preprint2004arXiv

Contamination and exclusion in the sigma Orionis young group

We present radial velocities for 38 low-mass candidate members of the sigma Orionis young group. We have measured their radial velocities by cross-correlation of high resolution (R~6000) AF2/WYFFOS spectra of the gravity sensitive NaI doublet at 8183, 8195Angstroms. The total sample contained 117 objects of which 54 have sufficient signal-to-noise to detect NaI at an equivalent width of 3Angstroms, however we only detect NaI in 38 of these. This implies that very low-mass members of this young group display weaker NaI absorption than similarly aged objects in the Upper Scorpius OB association. We develop a technique to assess membership using radial velocities with a range of uncertainties that does not bias the selection when large uncertainties are present. The resulting membership probabilities are used to assess the issue of exclusion in photometric selections, and we find that very few members are likely to be excluded by such techniques. We also assess the level of contamination in the expected pre-main sequence region of colour-magnitude space brighter than I = 17. We find that contamination by non-members in the expected PMS region of the colour-magnitude diagram is small. We conclude that although radial velocity alone is insufficient to confirm membership, high signal-to-noise observations of the NaI doublet provide the opportunity to use the strength of NaI absorption in concert with radial velocities to asses membership down to the lowest masses, where Lithium absorption no longer distinguishes youth.

preprint2003arXiv

Evidence for high accretion-rates in Weak-Line T Tauri stars?

We have discovered T Tauri stars which show startling spectral variability between observations seperated by 20 years. In spectra published by Bouvier & Appenzeller (1992) these objects showed very weak H-alpha emission, broad CaII absorption and so called ``composite spectra'', where the spectral type inferred from the blue region is earlier than that inferred from the red. We present here new spectroscopy which shows that all four stars now exhibit strong H-alpha emission, narrow CaII emission and a spectral type which is consistent at all wavelengths. We propose a scheme to understand these changes whereby the composite spectra of these stars can be explained by a period of active accretion onto the central, young star. In this scheme the composite spectrum consists of a contribution from the stellar photosphere and a contribution from a hot, optically thick, accretion component. The optically thick nature of the accretion flow explains the weakness of the H-alpha emission during this phase. Within this scheme, the change to a single spectral type at all wavelengths and emergence of strong H-alpha emission are consistent with the accretion columns becoming optically thin, as the accretion rate drops. There is a strong analogy here with the dwarf novae class of interacting binaries, which show similar behaviour during the decline from outbursts of high mass-transfer rate. The most important consequence of this interpretation is that these objects bring into question the association of Weak-Line T Tauri stars (WTTs) with non-accreting or discless objects. In light of this result we consider the justification for this paradigm.