BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

C. A. L. Bailer-Jones

C. A. L. Bailer-Jones appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.EPastro-ph.IMastro-phMachine Learningphysics.data-anphysics.geo-phPopulations and EvolutionApplicationsastro-ph.COBiological Physicsastro-ph.HENeural and Evolutionary Computingphysics.ao-ph

Catalog footprint

What is connected

37works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

37 published item(s)

preprint2022arXiv

Astrophysical Parameters from Gaia DR2, 2MASS & AllWISE

Stellar physical and dynamical properties are essential knowledge to understanding the structure, formation, and evolution of our Galaxy. We produced an all-sky uniformly derived catalog of stellar astrophysical parameters (APs; age, mass, temperature, bolometric luminosity, distance, dust extinction) to give insight into the physical properties of Milky-Way stars. Exploiting the power of multi-wavelength and multi-survey observations from Gaia DR2 parallaxes and integrated photometry along with 2MASS and AllWISE photometry, we introduce an all-sky uniformly derived catalog of stellar astrophysical parameters, including dust extinction (A0) and average grain size (R0) along the line of sight, for 123,097,070 stars. In contrast with previous works, we do not use a Galactic model as prior in our analysis. We validate our results against other literature (e.g., benchmark stars, interferometry, Bayestar, StarHorse). The limited optical information in the Gaia photometric bands or the lack of ultraviolet or spectroscopic information renders the chemistry inference prior dominated. We demonstrate that Gaia parallaxes bring sufficient leverage to explore the detailed structures of the interstellar medium in our Milky Way. In Gaia DR3, we will obtain the dispersed optical light information to break through some limitations of this analysis, allowing us to infer stellar chemistry in particular. Gaia promises us data to construct the most detailed view of the chemo-dynamics of field star populations in our Galaxy. Our catalog is available from GAVO at http://dc.g-vo.org/tableinfo/gdr2ap.main (soon Gaia Archive and VizieR)

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stars that approach within one parsec of the Sun: New and more accurate encounters identified in Gaia Data Release 3

Close encounters of stars to the Sun could affect life on Earth through gravitational perturbations of comets in the Oort cloud or exposure to ionizing radiation. By integrating orbits through the Galactic potential, I identify which of 33 million stars in Gaia DR3 with complete phase space information come close to the Sun. 61 stars formally approach within 1 pc, although there is high confidence in only 42 (two thirds) of these, the rest being spurious measurements or (in) binary systems. Most of the stars will encounter within the past or future 6 Myr; earlier/later encounters are less common due to the magnitude limit of the Gaia radial velocities (RVs). Several close encountering stars are identified for the first time, and the encounter times, distances, and velocities of previously known close encounters are determined more precisely on account of the significantly improved precision of Gaia DR3 over earlier releases. The K7 dwarf Gl 710 remains the closest known encounter, with an estimated (median) encounter distance of 0.0636 pc (90% confidence interval 0.0595-0.0678 pc) to take place in 1.3 Myr. The new second closest encounter took place 2.8 Myr ago: this was the G3 dwarf HD 7977, now 76 pc away, which approached within less than 0.05 pc of the Sun with a probability of one third. The apparent close encounter of the white dwarf UPM J0812-352 is probably spurious due to an incorrect RV in Gaia DR3.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Estimating distances from parallaxes. V: Geometric and photogeometric distances to 1.47 billion stars in Gaia Early Data Release 3

Stellar distances constitute a foundational pillar of astrophysics. The publication of 1.47 billion stellar parallaxes from Gaia is a major contribution to this. Yet despite Gaia's precision, the majority of these stars are so distant or faint that their fractional parallax uncertainties are large, thereby precluding a simple inversion of parallax to provide a distance. Here we take a probabilistic approach to estimating stellar distances that uses a prior constructed from a three-dimensional model of our Galaxy. This model includes interstellar extinction and Gaia's variable magnitude limit. We infer two types of distance. The first, geometric, uses the parallax together with a direction-dependent prior on distance. The second, photogeometric, additionally uses the colour and apparent magnitude of a star, by exploiting the fact that stars of a given colour have a restricted range of probable absolute magnitudes (plus extinction). Tests on simulated data and external validations show that the photogeometric estimates generally have higher accuracy and precision for stars with poor parallaxes. We provide a catalogue of 1.47 billion geometric and 1.35 billion photogeometric distances together with asymmetric uncertainty measures. Our estimates are quantiles of a posterior probability distribution, so they transform invariably and can therefore also be used directly in the distance modulus (5log10(r)-5). The catalogue may be downloaded or queried using ADQL at various sites (see http://www.mpia.de/homes/calj/gedr3_distances.html) where it can also be cross-matched with the Gaia catalogue.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Three-dimensional dust density structure of the Orion, Cygnus X,Taurus, and Perseus star-forming regions

Interstellar dust affects astronomical observations through absorption and reddening, yet this extinction is also a powerful tool for studying interstellar matter in galaxies. 3D reconstructions of dust extinction and density in the Milky Way have suffered from artefacts such as the fingers-of-god effect and negative densities, and have been limited by large computational costs. Here we aim to overcome these issues with a novel algorithm that derives the 3D extinction density of dust in the Milky Way using a latent variable Gaussian Process and variational inference. Our model maintains non-negative density and hence monotonically non-decreasing extinction along all lines-of-sight, and performs inference within a reasonable computational time. Using extinctions for hundreds of thousands of stars computed from optical and near-IR photometry, together with distances based on Gaia parallaxes, we use our algorithm to infer the structure of the Orion, Taurus, Perseus, and Cygnus X star-forming regions. A number of features superimposed in 2D extinction maps are clearly deblended in 3D dust extinction density maps. We find a large filament on the edge of Orion, identify filament that in the Taurus and Perseus regions, and show that Cygnus X is located at 1300-1500pc. We compute dust masses of the regions and find these to be slightly higher than previous estimates, likely a consequence of our input data recovering the highest column densities more effectively. Comparing our predicted extinctions to Planck data, we find that known relationships between density and dust processing, where high extinction lines of sight have the most processed grains, hold up in resolved observations when density is included, and that they exist at smaller scales than previously suggested. This can be used to study the changes in size or composition of dust as they are processed in molecular clouds. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evidence of a dynamically evolving Galactic warp

In a cosmological setting, the disc of a galaxy is expected to continuously experience gravitational torques and perturbations from a variety of sources, which can cause the disc to wobble, flare and warp. Specifically, the study of galactic warps and their dynamical nature can potentially reveal key information on the formation history of galaxies and the mass distribution of their halos. Our Milky Way presents a unique case study for galactic warps, thanks to the detailed knowledge of its stellar distribution and kinematics. Using a simple model of how the warp's orientation is changing with time, we here measure the precession rate of the Milky Way's warp using 12 million giant stars from Gaia Data Release 2, finding that it is precessing at $10.86 \pm 0.03_{stat} \pm 3.20_{sys}$ km/s/kpc in the direction of Galactic rotation, about one third the angular velocity at the Sun's position in the Galaxy. The direction and magnitude of the warp's precession rate favour the scenario that the warp is the results of a recent or ongoing encounter with a satellite galaxy, rather than the relic of the ancient assembly history of the Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dust in three dimensions in the Galactic Plane

We present three dimensional maps in monochromatic extinction $A_{\rm 0}$ and the extinction parameter $R_0$ within a few degrees of the Galactic plane. These are inferred using photometry from the Pan-STARRS1 and Spitzer Glimpse surveys of nearly $20$ million stars located in the region $l = 0-250^{\circ}$ and from $b = -4.5^{\circ}$ to $b=4.5^{\circ}$. Given the available stellar number density, we use an angular resolution of $7' \times 7'$ and steps of $1{\rm mag}$ in distance modulus. We simultaneously estimate distance modulus and effective temperature $T_{\rm eff}$ alongside the other parameters for stars individually using the method of \citet{Hanson2014} before combining these estimates to a complete map. The full maps are available via the MNRAS website.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Estimating distances from parallaxes

Astrometric surveys such as Gaia and LSST will measure parallaxes for hundreds of millions of stars. Yet they will not measure a single distance. Rather, a distance must be estimated from a parallax. In this didactic article, I show that doing this is not trivial once the fractional parallax error is larger than about 20%, which will be the case for about 80% of stars in the Gaia catalogue. Estimating distances is an inference problem in which the use of prior assumptions is unavoidable. I investigate the properties and performance of various priors and examine their implications. A supposed uninformative uniform prior in distance is shown to give very poor distance estimates (large bias and variance). Any prior with a sharp cut-off at some distance has similar problems. The choice of prior depends on the information one has available - and is willing to use - concerning, for example, the survey and the Galaxy. I demonstrate that a simple prior which decreases asymptotically to zero at infinite distance has good performance, accommodates non-positive parallaxes, and does not require a bias correction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finding the imprints of stellar encounters in long period comets

The solar system's Oort cloud can be perturbed by the Galactic tide and by individual passing stars. These perturbations can inject Oort cloud objects into the inner parts of the solar system, where they may be observed as the long-period comets (periods longer than 200 years). Using dynamical simulations of the Oort cloud under the perturbing effects of the tide and 61 known stellar encounters, we investigate the link between long-period comets and encounters. We find that past encounters were responsible for injecting at least 5% of the currently known long-period comets. This is a lower limit due to the incompleteness of known encounters. Although the Galactic tide seems to play the dominant role in producing the observed long-period comets, the non-uniform longitude distribution of the cometary perihelia suggests the existence of strong -- but as yet unidentified -- stellar encounters or other impulses. The strongest individual future and past encounters are probably HIP 89825 (Gliese 710) and HIP 14473, which contribute at most 8% and 6% to the total flux of long-period comets, respectively. Our results show that the strength of an encounter can be approximated well by a simple proxy, which will be convenient for quickly identifying significant encounters in large data sets. Our analysis also indicates a smaller population of the Oort cloud than is usually assumed, which would bring the mass of the solar nebula into line with planet formation theories.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Obliquity and precession as pacemakers of Pleistocene deglaciations

The Milankovitch theory states that the orbital eccentricity, precession, and obliquity of the Earth influence our climate by modulating the summer insolation at high latitudes in the northern hemisphere. Despite considerable success of this theory in explaining climate change over the Pleistocene epoch (2.6 to 0.01 Myr ago), it is inconclusive with regard to which combination of orbital elements paced the 100 kyr glacial-interglacial cycles over the late Pleistocene. Here we explore the role of the orbital elements in pacing the Pleistocene deglaciations by modeling ice-volume variations in a Bayesian approach. When comparing models, this approach takes into account the uncertainties in the data as well as the different degrees of model complexity. We find that the Earth's obliquity (axial tilt) plays a dominant role in pacing the glacial cycles over the whole Pleistocene, while precession only becomes important in pacing major deglaciations after the transition of the dominant period from 41 kyr to 100 kyr (the mid-Pleistocene transition). We also find that geomagnetic field and orbital inclination variations are unlikely to have paced the Pleistocene deglaciations. We estimate that the mid-Pleistocene transition took place over a 220 kyr interval centered on a time 715 kyr ago, although the data permit a range of 600--1000 kyr. This transition, occurring within just two 100\,kyr cycles, indicates a relatively rapid change in the climate response to insolation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Close encounters of the stellar kind

Stars which pass close to the Sun can perturb the Oort cloud, injecting comets into the inner solar system where they may collide with the Earth. Using van Leeuwen's re-reduction of the Hipparcos data complemented by the original Hipparcos and Tycho-2 catalogues, along with recent radial velocity surveys, I integrate the orbits of over 50 000 stars through the Galaxy to look for close encounters. The search uses a Monte Carlo simulation over the covariance of the data in order to properly characterize the uncertainties in the times, distances, and speeds of the encounters. I show that modelling stellar encounters by assuming instead a linear relative motion produces, for many encounters, inaccurate and biased results. I find 42, 14, and 4 stars which have encounter distances below 2, 1, and 0.5 pc respectively, although some of these stars have questionable data. Of the 14 stars coming within 1 pc, 5 were found by at least one of three previous studies (which found a total of 7 coming within 1 pc). The closest encounter appears to be Hip 85605, a K or M star, which has a 90% probability of coming between 0.04 and 0.20 pc between 240 and 470 kyr from now (90% Bayesian confidence interval). However, its astrometry may be incorrect, in which case the closest encounter found is the K7 dwarf GL 710, which has a 90% probability of coming within 0.10-0.44 pc in about 1.3 Myr. A larger perturbation may have been caused by gamma Microscopii, a G6 giant with a mass of about 2.5 Msol, which came within 0.35-1.34 pc (90% confidence interval) around 3.8 Myr ago.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exploring the role of the Sun's motion in terrestrial comet impacts

The cratering record on the Earth and Moon shows that our planet has been exposed to high velocity impacts for much or all of its existence. Some of these craters were produced by the impact of long period comets (LPCs). These probably originated in the Oort cloud, and were put into their present orbits through gravitational perturbations arising from the Galactic tide and stellar encounters, both of which are modulated by the solar motion about the Galaxy. Here we construct dynamicalmodels of these mechanisms in order to predict the time-varying impact rate of LPCs and the angular distribution of their perihelia (which is observed to be non-uniform). Comparing the predictions of these dynamical models with other models, we conclude that cometary impacts induced by the solar motion contribute only a small fraction of terrestrial impact craters over the past 250 Myr. Over this time scale the apparent cratering rate is dominated by a secular increase towards the present, which might be the result of the disruption of a large asteroid. Our dynamical models, together with the solar apex motion, predict a non-uniform angular distribution of the perihelia, without needing to invoke the existence of a massive body in the outer Oort cloud. Our results are reasonably robust to changes in the parameters of the Galaxy model, Oort cloud, and stellar encounters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

3D Galactic dust extinction mapping with multi-band photometry

We present a method to simultaneously infer the interstellar extinction parameters $A_0$ and $R_0$, stellar effective temperature $T_{\rm eff}$, and distance modulus $μ$ in a Bayesian framework. Using multi-band photometry from SDSS and UKIDSS, we train a forward model to emulate the colour-change due to physical properties of stars and the interstellar medium for temperatures from 4000 to 9000 K and extinctions from 0 to 5 mag. We introduce a Hertzsprung-Russel diagram prior to account for physical constraints on the distribution of stars in the temperature-absolute magnitude plane. This allows us to infer distances probabilistically. Influences of colour information, priors and model parameters are explored. Residual mean absolute errors (MAEs) on a set of objects for extinction and temperature are 0.2 mag and 300 K, respectively, for $R_0$ fixed to 3.1. For variable $R_0$, we obtain MAEs of 0.37 mag, 412.9 K and 0.74 for $A_0$, $T_{\rm eff}$ and $R_0$, respectively. Distance moduli are accurate to approximately 2 mag. Quantifying the precisions of individual parameter estimates with $68\%$ confidence interval of the posterior distribution, we obtain 0.05 mag, 66 K, 2 mag and 0.07 for $A_0$, $T_{\rm eff}$, $μ$ and $R_0$, respectively, although we find that these underestimate the accuracy of the model. We produce two-dimensional maps in extinction and $R_0$ that are compared to previous work. Furthermore we incorporate the inferred distance information to compute fully probabilistic distance profiles for individual lines of sight. The individual stellar AP estimates, combined with inferred 3D information will make possible many Galactic science and modelling applications. Adapting our method to work with other surveys, such as Pan-STARRS and Gaia, will allow us to probe other regions of the Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A history of the gamma-ray burst flux at the Earth from Galactic globular clusters

Nearby gamma-ray bursts (GRBs) are likely to have represented a significant threat to life on the Earth. Recent observations suggest that a significant source of such bursts is compact binary mergers in globular clusters. This link between globular clusters and GRBs offers the possibility to find time intervals in the past with higher probabilities of a nearby burst, by tracing globular cluster orbits back in time. Here we show that the expected flux from such bursts is not flat over the past 550 Myr but rather exhibits three broad peaks, at 70, 180 and 340 Myr ago. The main source for nearby GRBs for all three time intervals is the globular cluster 47 Tuc, a consequence of its large mass and high stellar encounter rate, as well as the fact that it is one of the globular clusters which comes quite close to the Sun. Mass extinction events indeed coincide with all three time intervals found in this study, although a chance coincidence is quite likely. Nevertheless, the identified time intervals can be used as a guide to search for specific signatures of GRBs in the geological record around these times.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Assessing the influence of the solar orbit on terrestrial biodiversity

The terrestrial fossil record shows a significant variation in the extinction and origination rates of species during the past half billion years. Numerous studies have claimed an association between this variation and the motion of the Sun around the Galaxy, invoking the modulation of cosmic rays, gamma rays and comet impact frequency as a cause of this biodiversity variation. However, some of these studies exhibit methodological problems, or were based on coarse assumptions (such as a strict periodicity of the solar orbit). Here we investigate this link in more detail, using a model of the Galaxy to reconstruct the solar orbit and thus a predictive model of the temporal variation of the extinction rate due to astronomical mechanisms. We compare these predictions as well as those of various reference models with paleontological data. Our approach involves Bayesian model comparison, which takes into account the uncertainties in the paleontological data as well as the distribution of solar orbits consistent with the uncertainties in the astronomical data. We find that various versions of the orbital model are not favored beyond simpler reference models. In particular, the distribution of mass extinction events can be explained just as well by a uniform random distribution as by any other model tested. Although our negative results on the orbital model are robust to changes in the Galaxy model, the Sun's coordinates and the errors in the data, we also find that it would be very difficult to positively identify the orbital model even if it were the true one. (In contrast, we do find evidence against simpler periodic models.) Thus while we cannot rule out there being some connection between solar motion and biodiversity variations on the Earth, we conclude that it is difficult to give convincing positive conclusions of such a connection using current data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evidence for periodicities in the extinction record? Response to Melott & Bambach [arXiv:1307.1884]

In a recent paper, we applied time series analysis methods to study the possible influence of the solar motion through the Galaxy on terrestrial extinction (Feng & Bailer-Jones 2013). We drew conclusions about the relative probabilities of how well different models explain the geological extinction record. We found no strong evidence for either a periodicity in the extinction record or for a link between the solar orbit and this extinction record. In response to this, Melott & Bambach (2013) dedicate the majority of their article to criticizing our analysis. Their main objection is to evidence-based model comparison. Here we deal with their criticisms one by one, and show that they are largely misplaced, leaving our conclusions unaffected. The bottom line is that although one can get a periodic model at some period to fit the extinction data, there are other models which explain the data better.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Gaia astrophysical parameters inference system (Apsis). Pre-launch description

The Gaia satellite will survey the entire celestial sphere down to 20th magnitude, obtaining astrometry, photometry, and low resolution spectrophotometry on one billion astronomical sources, plus radial velocities for over one hundred million stars. Its main objective is to take a census of the stellar content of our Galaxy, with the goal of revealing its formation and evolution. Gaia's unique feature is the measurement of parallaxes and proper motions with hitherto unparalleled accuracy for many objects. As a survey, the physical properties of most of these objects are unknown. Here we describe the data analysis system put together by the Gaia consortium to classify these objects and to infer their astrophysical properties using the satellite's data. This system covers single stars, (unresolved) binary stars, quasars, and galaxies, all covering a wide parameter space. Multiple methods are used for many types of stars, producing multiple results for the end user according to different models and assumptions. Prior to its application to real Gaia data the accuracy of these methods cannot be assessed definitively. But as an example of the current performance, we can attain internal accuracies (RMS residuals) on F,G,K,M dwarfs and giants at G=15 (V=15-17) for a wide range of metallicites and interstellar extinctions of around 100K in effective temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]), and 0.25dex in surface gravity (logg). The accuracy is a strong function of the parameters themselves, varying by a factor of more than two up or down over this parameter range. After its launch in November 2013, Gaia will nominally observe for five years, during which the system we describe will continue to evolve in light of experience with the real data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

What will Gaia tell us about the Galactic disk?

Gaia will provide parallaxes and proper motions with accuracy ranging from 10 to 1000 microarcsecond on up to one billion stars. Most of these will be disk stars: for an unreddened K giant at 6 kpc, it will measure the distance accurate to 15% and the transverse velocity to an accuracy of about 1 km/s. Gaia will observe tracers of Galactic structure across the whole HR diagram, including Cepheids, RR Lyrae, white dwarfs, F dwarfs and HB stars. Onboard low resolution spectrophotometry will permit -- in addition to a Teff estimate -- dwarf/giant discrimination, metallicity measurement and extinction determination. For the first time, then, Gaia will provide us with a 3D spatial/properties map and at least a 2D velocity map of these tracers (RVs will be obtained too for brighter stars.) This will be a goldmine of information from which to learn about the origin and evolution of the Galactic disk. I briefly review the Gaia mission, and then show how the expected astrometric accuracies translate into distance and velocity accuracies and statistics. I examine the impact Gaia should have on a few scientific areas relevant to the Galactic disk. I discuss how a better determination of the spiral arm locations and pattern speed, plus a better reconstruction of the Sun's orbit over the past billion years (from integration through the Gaia-measured gravitational potential) will allow us to assess the possible role of spiral arm crossings in ice ages and mass extinctions on the Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Bayesian method for the analysis of deterministic and stochastic time series

I introduce a general, Bayesian method for modelling univariate time series data assumed to be drawn from a continuous, stochastic process. The method accommodates arbitrary temporal sampling, and takes into account measurement uncertainties for arbitrary error models (not just Gaussian) on both the time and signal variables. Any model for the deterministic component of the variation of the signal with time is supported, as is any model of the stochastic component on the signal and time variables. Models illustrated here are constant and sinusoidal models for the signal mean combined with a Gaussian stochastic component, as well as a purely stochastic model, the Ornstein-Uhlenbeck process. The posterior probability distribution over model parameters is determined via Monte Carlo sampling. Models are compared using the "cross-validation likelihood", in which the posterior-averaged likelihood for different partitions of the data are combined. In principle this is more robust to changes in the prior than is the evidence (the prior-averaged likelihood). The method is demonstrated by applying it to the light curves of 11 ultra cool dwarf stars, claimed by a previous study to show statistically significant variability. This is reassessed here by calculating the cross-validation likelihood for various time series models, including a null hypothesis of no variability beyond the error bars. 10 of 11 light curves are confirmed as being significantly variable, and one of these seems to be periodic, with two plausible periods identified. Another object is best described by the Ornstein-Uhlenbeck process, a conclusion which is obviously limited to the set of models actually tested.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The expected performance of stellar parametrization with Gaia spectrophotometry

Gaia will obtain astrometry and spectrophotometry for essentially all sources in the sky down to a broad band magnitude limit of G=20, an expected yield of 10^9 stars. Its main scientific objective is to reveal the formation and evolution of our Galaxy through chemo-dynamical analysis. In addition to inferring positions, parallaxes and proper motions from the astrometry, we must also infer the astrophysical parameters of the stars from the spectrophotometry, the BP/RP spectrum. Here we investigate the performance of three different algorithms (SVM, ILIUM, Aeneas) for estimating the effective temperature, line-of-sight interstellar extinction, metallicity and surface gravity of A-M stars over a wide range of these parameters and over the full magnitude range Gaia will observe (G=6-20mag). One of the algorithms, Aeneas, infers the posterior probability density function over all parameters, and can optionally take into account the parallax and the Hertzsprung-Russell diagram to improve the estimates. For all algorithms the accuracy of estimation depends on G and on the value of the parameters themselves, so a broad summary of performance is only approximate. For stars at G=15 with less than two magnitudes extinction, we expect to be able to estimate Teff to within 1%, logg to 0.1-0.2dex, and [Fe/H] (for FGKM stars) to 0.1-0.2dex, just using the BP/RP spectrum (mean absolute error statistics are quoted). Performance degrades at larger extinctions, but not always by a large amount. Extinction can be estimated to an accuracy of 0.05-0.2mag for stars across the full parameter range with a priori unknown extinction between 0 and 10mag. Performance degrades at fainter magnitudes, but even at G=19 we can estimate logg to better than 0.2dex for all spectral types, and [Fe/H] to within 0.35dex for FGKM stars, for extinctions below 1mag.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A semi-empirical library of galaxy spectra for Gaia classification based on SDSS data and PEGASE models

Aims:This paper is the third in a series implementing a classification system for Gaia observations of unresolved galaxies. The system makes use of template galaxy spectra in order to determine spectral classes and estimate intrinsic astrophysical parameters. In previous work we used synthetic galaxy spectra produced by PEGASE.2 code to simulate Gaia observations and to test the performance of Support Vector Machine (SVM) classifiers and parametrizers. Here we produce a semi-empirical library of galaxy spectra by fitting SDSS spectra with the previously produced synthetic libraries. We present (1) the semi-empirical library of galaxy spectra, (2) a comparison between the observed and synthetic spectra, and (3) first results of claassification and parametrization experiments with simulated Gaia spectrophotometry of this library. Methods: We use chi2-fitting to fit SDSS galaxy spectra with the synthetic library in order to construct a semi-empirical library of galaxy spectra in which (1) the real spectra are extended by the synthetic ones in order to cover the full wavelength range of Gaia, and (2) astrophysical parameters are assigned to the SDSS spectra by the best fitting synthetic spectrum. The SVM models were trained with and applied to semi-empirical spectra. Tests were performed for the classification of spectral types and the estimation of the most significant galaxy parameters (in particular redshift, mass to light ratio and star formation history). Results: We produce a semi-empirical library of 33670 galaxy spectra covering the wavelength range 250 to 1050 nm at a sampling of 1 nm or less. Using the results of the fitting of the SDSS spectra with our synthetic library, we investigate the range of the input model parameters that produces spectra which are in good agreement with observations. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bayesian time series analysis of terrestrial impact cratering

Giant impacts by comets and asteroids have probably had an important influence on terrestrial biological evolution. We know of around 180 high velocity impact craters on the Earth with ages up to 2400Myr and diameters up to 300km. Some studies have identified a periodicity in their age distribution, with periods ranging from 13 to 50Myr. It has further been claimed that such periods may be causally linked to a periodic motion of the solar system through the Galactic plane. However, many of these studies suffer from methodological problems, for example misinterpretation of p-values, overestimation of significance in the periodogram or a failure to consider plausible alternative models. Here I develop a Bayesian method for this problem in which impacts are treated as a stochastic phenomenon. Models for the time variation of the impact probability are defined and the evidence for them in the geological record is compared using Bayes factors. This probabilistic approach obviates the need for ad hoc statistics, and also makes explicit use of the age uncertainties. I find strong evidence for a monotonic decrease in the recorded impact rate going back in time over the past 250Myr for craters larger than 5km. The same is found for the past 150Myr when craters with upper age limits are included. This is consistent with a crater preservation/discovery bias modulating an otherwise constant impact rate. The set of craters larger than 35km (so less affected by erosion and infilling) and younger than 400Myr are best explained by a constant impact probability model. A periodic variation in the cratering rate is strongly disfavoured in all data sets. There is also no evidence for a periodicity superimposed on a constant rate or trend, although this more complex signal would be harder to distinguish.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Photometric Classification Server for Pan-STARRS1

The Pan-STARRS1 survey is obtaining multi-epoch imaging in 5 bands (gps rps ips zps yps) over the entire sky North of declination -30deg. We describe here the implementation of the Photometric Classification Server (PCS) for Pan-STARRS1. PCS will allow the automatic classification of objects into star/galaxy/quasar classes based on colors, the measurement of photometric redshifts for extragalactic objects, and constrain stellar parameters for stellar objects, working at the catalog level. We present tests of the system based on high signal-to-noise photometry derived from the Medium Deep Fields of Pan-STARRS1, using available spectroscopic surveys as training and/or verification sets. We show that the Pan-STARRS1 photometry delivers classifications and photometric redshifts as good as the Sloan Digital Sky Survey (SDSS) photometry to the same magnitude limits. In particular, our preliminary results, based on this relatively limited dataset down to the SDSS spectroscopic limits and therefore potentially improvable, show that stars are correctly classified as such in 85% of cases, galaxies in 97% and QSOs in 84%. False positives are less than 1% for galaxies, ~19% for stars and ~28% QSOs. Moreover, photometric redshifts for 1000 luminous red galaxies up to redshift 0.5 are determined to 2.4% precision with just 0.4% catastrophic outliers and small (-0.5%) residual bias. PCS will create a value added catalog with classifications and photometric redshifts for eventually many millions sources.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The substellar population of Sigma Orionis: A deep wide survey

We present a deep I,Z photometric survey covering a total area of 1.12 deg^{2} of the Sigma Orionis cluster (Icompl=22 and Zcompl=21.5mag). From I, I-Z color-magnitude diagrams we have selected 153 candidates that fit the previously known sequence of the cluster. Using J-band photometry, we find that 124 of the 151 candidates follow the previously known infrared photometric sequence of the cluster and are probably members. We have studied the spatial distribution of these candidates and found that there are objects located at distances greater than 30 arcmin to the north and west of Sigma Orionis that probably belong to different populations of the Orion's Belt. For the 102 bona fide Sigma Orionis cluster member candidates, we find that the radial surface density can be represented by a decreasing exponential function (sigma = sigma_0 e^{-r/r_0}) with a central density of sigma_0=0.23+/-0.03 object/arcmin^{2} and a characteristic radius of r_0=9.5+/-0.7 arcmin. From a statistical comparison with Monte Carlo simulations, we conclude that the spatial distribution of the cluster member candidates is compatible with a Poissonian distribution and, hence, they are not mainly forming aggregations or sub-clustering. Using near-infrared JHK-band data from 2MASS and UKIDSS and mid-infrared data from IRAC/Spitzer, we find that 5-9 % of the brown dwarf candidates in the Sigma Orionis cluster have K-band excesses and 31+/-7 % of them show mid-infrared excesses at wavelengths longer than 5.8 microns, which are probably related to the presence of disks. We have also calculated the initial mass spectrum (dN/dm) of Sigma Orionis from very low mass stars (0.10 Msol) to the deuterium-burning mass limit (0.012-0.013 Msol). This is a rising function toward lower masses and can be represented by a power-law distribution (dN/dm = m^{-alpha}) with an exponent alpha of 0.7+/-0.3 for an age of 3 Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bayesian inference of stellar parameters and interstellar extinction using parallaxes and multiband photometry

Astrometric surveys provide the opportunity to measure the absolute magnitudes of large numbers of stars, but only if the individual line-of-sight extinctions are known. Unfortunately, extinction is highly degenerate with stellar effective temperature when estimated from broad band optical/infrared photometry. To address this problem, I introduce a Bayesian method for estimating the intrinsic parameters of a star and its line-of-sight extinction. It uses both photometry and parallaxes in a self-consistent manner in order to provide a non-parametric posterior probability distribution over the parameters. The method makes explicit use of domain knowledge by employing the Hertzsprung--Russell Diagram (HRD) to constrain solutions and to ensure that they respect stellar physics. I first demonstrate this method by using it to estimate effective temperature and extinction from BVJHK data for a set of artificially reddened Hipparcos stars, for which accurate effective temperatures have been estimated from high resolution spectroscopy. Using just the four colours, we see the expected strong degeneracy (positive correlation) between the temperature and extinction. Introducing the parallax, apparent magnitude and the HRD reduces this degeneracy and improves both the precision (reduces the error bars) and the accuracy of the parameter estimates, the latter by about 35%. The resulting accuracy is about 200K in temperature and 0.2mag in extinction. I then apply the method to estimate these parameters and absolute magnitudes for some 47000 F,G,K Hipparcos stars which have been cross-matched with 2MASS. The method can easily be extended to incorporate the estimation of other parameters, in particular metallicity and surface gravity, making it particularly suitable for the analysis of the 10^9 stars from Gaia.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The ILIUM forward modelling algorithm for multivariate parameter estimation and its application to derive stellar parameters from Gaia spectrophotometry

I introduce an algorithm for estimating parameters from multidimensional data based on forward modelling. In contrast to many machine learning approaches it avoids fitting an inverse model and the problems associated with this. The algorithm makes explicit use of the sensitivities of the data to the parameters, with the goal of better treating parameters which only have a weak impact on the data. The forward modelling approach provides uncertainty (full covariance) estimates in the predicted parameters as well as a goodness-of-fit for observations. I demonstrate the algorithm, ILIUM, with the estimation of stellar astrophysical parameters (APs) from simulations of the low resolution spectrophotometry to be obtained by Gaia. The AP accuracy is competitive with that obtained by a support vector machine. For example, for zero extinction stars covering a wide range of metallicity, surface gravity and temperature, ILIUM can estimate Teff to an accuracy of 0.3% at G=15 and to 4% for (lower signal-to-noise ratio) spectra at G=20. [Fe/H] and logg can be estimated to accuracies of 0.1-0.4dex for stars with G<=18.5. If extinction varies a priori over a wide range (Av=0-10mag), then Teff and Av can be estimated quite accurately (3-4% and 0.1-0.2mag respectively at G=15), but there is a strong and ubiquitous degeneracy in these parameters which limits our ability to estimate either accurately at faint magnitudes. Using the forward model we can map these degeneracies (in advance), and thus provide a complete probability distribution over solutions. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The substellar mass function in the central region of the open cluster Praesepe from deep LBT observations

Studies of the mass function (MF) of open clusters of different ages allow us to probe the efficiency with which brown dwarfs (BDs) are evaporated from clusters to populate the field. Surveys in old clusters (age > 100 Myr) do not suffer so severely from several problems encountered in young clusters, such as intra-cluster extinction and large uncertainties in BD models. Here we present the results of a deep photometric survey to study the MF of the old open cluster Praesepe (age 590 Myr and distance 190 pc), down to a 5 sigma detection limit at i~25.6 mag (~40M_Jup). We identify 62 cluster member candidates, of which 40 are substellar, from comparison with predictions from a dusty atmosphere model. The MF rises from the substellar boundary until ~60M_Jup and then declines. This is quite different from the form inferred for other open clusters older than 50 Myr, but seems to be similar to those found in very young open cluster, whose MFs peak at ~10M_Jup. Either Praesepe really does have a different MF from other clusters or they had similar initial MFs but have differed in their dynamical evolution. We further have identified six foreground T dwarf candidates towards Praesepe, which require follow-up spectroscopy to confirm their nature.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Constraint on brown dwarf formation via ejection: radial variation of the stellar and substellar mass function of the young open cluster IC2391

We present the stellar and substellar mass function of the open cluster IC2391, plus its radial dependence, and use this to put constraints on the formation mechanism of brown dwarfs. Our multiband optical and infrared photometric survey with spectroscopic follow-up covers 11 square degrees, making it the largest survey of this cluster to date. We observe a radial variation in the mass function over the range 0.072 to 0.3Msol, but no significant variation in the mass function below the substellar boundary at the three cluster radius intervals analyzed. This lack of radial variation for low masses is what we would expect with the ejection scenario for brown dwarf formation, although considering that IC2391 has an age about three times older than its crossing time, we expect that brown dwarfs with a velocity greater than the escape velocity have already escaped the cluster. Alternatively, the variation in the mass function of the stellar objects could be an indication that they have undergone mass segregation via dynamical evolution. We also observe a significant variation across the cluster in the colour of the (background) field star locus in colour-magnitude diagrams and conclude that this is due to variable background extinction in the Galactic plane. From our preliminary spectroscopic follow-up to confirm brown dwarf status and cluster membership, we find that all candidates are M dwarfs (in either the field or the cluster), demonstrating the efficiency of our photometric selection method in avoiding contaminants (e.g. red giants). About half of our photometric candidates for which we have spectra are spectroscopically-confirmed as cluster members; two are new spectroscopically-confirmed brown dwarf members of IC2391.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Brown dwarfs and very low mass stars in the Praesepe open cluster: a dynamically unevolved mass function?

[Abridged] In this paper, we present the results of a photometric survey to identify low mass and brown dwarf members of the old open cluster Praesepe (age of 590[+150][-120]Myr and distance of 190[+6.0][-5.8]pc) and use this to infer its mass function which we compare with that of other clusters. We have performed an optical (Ic-band) and near-infrared (J and Ks-band) photometric survey of Praesepe with a spatial coverage of 3.1deg^2. With 5sigma detection limits of Ic=23.4 and J=20.0, our survey is sensitive to objects with masses from about 0.6 to 0.05Msol. The mass function of Praesepe rises from 0.6Msol down to 0.1Msol and then turns-over at ~0.1Msol. The rise observed is in agreement with the mass function derived by previous studies, including a survey based on proper motion and photometry. Comparing our mass function with that for another open cluster with a similar age, the Hyades (age ~ 600Myr), we see a significant difference. Possible reasons are that dynamical evaporation has not influenced the Hyades and Praesepe in the same way, or that the clusters did not have the same initial mass function, or that dynamical interactions have modified the evolution of one or both clusters. Although a difference in the binary fractions of the clusters could cause the observed (i.e. system) mass functions to differ, measurements in the literature give no evidence for a significant difference in the binary fractions of the two clusters. Of our cluster candidates, six have masses predicted to be equal to or below the stellar/substellar boundary at 0.072Msol.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Candidate free-floating super-Jupiters in the young sigma Orionis open cluster

Free-floating substellar candidates with estimated theoretical masses of as low as ~5 Jupiter masses have been found in the ~3 Myr old sigma Orionis open cluster. As the overlap with the planetary mass domain increases, the question of how these objects form becomes important. The determination of their number density and whether a mass cut-off limit exists is crucial to understanding their formation. We propose to search for objects of yet lower masses in the cluster and determine the shape of the mass function at low mass. Using new- and (re-analysed) published IZJHKs[3.6]-[8.0]-band data of an area of 840 arcmin2, we performed a search for LT-type cluster member candidates in the magnitude range J=19.5-21.5 mag, based on their expected magnitudes and colours. Besides recovering the T type object S Ori 70 and two other known objects, we find three new cluster member candidates, S Ori 72-74, with J=21 mag and within 12 arcmin of the cluster centre. They have theoretical masses of 4 (-2,+3) M_Jup and are among the least massive free-floating objects detected by direct imaging outside the Solar System. The photometry in archival Spitzer [3.6]-[5.8]-band images infers that S Ori 72 is an L/T transition candidate and S Ori 73 a T-type candidate, following the expected cluster sequence in the mid-infrared. Finally, the L-type candidate S Ori 74 with lower quality photometry is located at 11.8 arcsec (~4250 AU) of a stellar member of sigma Orionis and could be a companion. After contaminant correction in the area complete to J=21.1 mag, we estimate that there remain between zero and two cluster members in the mass interval 6-4 M_Jup. Our result suggests a possible turnover in the substellar mass spectrum below ~6 Jupiter masses, which could be investigated further by wider and deeper photometric surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Polarisation of very-low-mass stars and brown dwarfs

Ultra-cool dwarfs of the L spectral type (Teff=1400-2200K) are known to have dusty atmospheres. Asymmetries of the dwarf surface may arise from rotationally-induced flattening and dust-cloud coverage, and may result in non-zero linear polarisation through dust scattering. We aim to study the heterogeneity of ultra-cool dwarfs' atmospheres and the grain-size effects on the polarisation degree in a sample of nine late M, L and early T dwarfs. We obtain linear polarimetric imaging measurements using FORS1 at the Very Large Telescope, in the Bessel I filter, and for a subset in the Bessel R and the Gunn z filters. We measure a polarisation degree of (0.31+/-0.06)% for LHS102BC. We fail to detect linear polarisation in the rest of our sample, with upper-limits on the polarisation degree of each object of 0.09% to 0.76% (95% CL). For those targets we do not find evidence of large-scale cloud horizontal structure in our data. Together with previous surveys, our results set the fraction of ultra-cool dwarfs with detected linear polarisation to (30+10-6)% (1-sigma). For three brown dwarfs, our observations indicate polarisation degrees different (at the 3-sigma level) than previously reported, giving hints of possible variations. Our results fail to correlate with the current model predictions for ultra-cool dwarf polarisation for a flattening-induced polarisation, or with the variability studies for a polarisation induced by an hetereneous cloud cover. This stresses the intricacy of each of those tasks, but may as well proceed from complex and dynamic atmospheric processes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The evidence for and against astronomical impacts on climate change and mass extinctions: A review

Numerous studies over the past 30 years have suggested there is a causal connection between the motion of the Sun through the Galaxy and terrestrial mass extinctions or climate change. Proposed mechanisms include comet impacts (via perturbation of the Oort cloud), cosmic rays and supernovae, the effects of which are modulated by the passage of the Sun through the Galactic midplane or spiral arms. Supposed periodicities in the fossil record, impact cratering dates or climate proxies over the Phanerozoic (past 545 Myr) are frequently cited as evidence in support of these hypotheses. This remains a controversial subject, with many refutations and replies having been published. Here I review both the mechanisms and the evidence for and against the relevance of astronomical phenomena to climate change and evolution. This necessarily includes a critical assessment of time series analysis techniques and hypothesis testing. Some of the studies have suffered from flaws in methodology, in particular drawing incorrect conclusions based on ruling out a null hypothesis. I conclude that there is little evidence for intrinsic periodicities in biodiversity, impact cratering or climate on timescales of tens to hundreds of Myr. Furthermore, Galactic midplane and spiral arm crossings seem to have little or no impact on biological or climate variation above background level. (truncated)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Towards a library of synthetic galaxy spectra and preliminary results of classification and parametrization of unresolved galaxies for Gaia - II

This paper is the second in a series, implementing a classification system for Gaia observations of unresolved galaxies. Our goals are to determine spectral classes and estimate intrinsic astrophysical parameters via synthetic templates. Here we describe (1) a new extended library of synthetic galaxy spectra, (2) its comparison with various observations, and (3) first results of classification and parametrization experiments using simulated Gaia spectrophotometry of this library. Using the PEGASE.2 code, based on galaxy evolution models that take account of metallicity evolution, extinction correction, and emission lines (with stellar spectra based on the BaSeL library), we improved our first library and extended it to cover the domain of most of the SDSS catalogue. We produce an extended library of 28885 synthetic galaxy spectra at zero redshift covering four general Hubble types of galaxies, over the wavelength range between 250 and 1050 nm at a sampling of 1 nm or less. The library is also produced for 4 random values of redshift in the range of 0-0.2. It is computed on a random grid of four key astrophysical parameters (infall timescale and 3 parameters defining the SFR) and, depending on the galaxy type, on two values of the age of the galaxy. The synthetic library was compared and found to be in good agreement with various observations. The first results from the SVM classifiers and parametrizers are promising, indicating that Hubble types can be reliably predicted and several parameters estimated with low bias and variance.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Finding rare objects and building pure samples: Probabilistic quasar classification from low resolution Gaia spectra

We develop and demonstrate a probabilistic method for classifying rare objects in surveys with the particular goal of building very pure samples. It works by modifying the output probabilities from a classifier so as to accommodate our expectation (priors) concerning the relative frequencies of different classes of objects. We demonstrate our method using the Discrete Source Classifier, a supervised classifier currently based on Support Vector Machines, which we are developing in preparation for the Gaia data analysis. DSC classifies objects using their very low resolution optical spectra. We look in detail at the problem of quasar classification, because identification of a pure quasar sample is necessary to define the Gaia astrometric reference frame. By varying a posterior probability threshold in DSC we can trade off sample completeness and contamination. We show, using our simulated data, that it is possible to achieve a pure sample of quasars (upper limit on contamination of 1 in 40,000) with a completeness of 65% at magnitudes of G=18.5, and 50% at G=20.0, even when quasars have a frequency of only 1 in every 2000 objects. The star sample completeness is simultaneously 99% with a contamination of 0.7%. Including parallax and proper motion in the classifier barely changes the results. We further show that not accounting for class priors in the target population leads to serious misclassifications and poor predictions for sample completeness and contamination. (Truncated)

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The substellar mass function in sigma Orionis. II. Optical, near-infrared and IRAC/Spitzer photometry of young cluster brown dwarfs and planetary-mass objects

We investigate the mass function in the substellar domain down to a few Jupiter masses in the young sigma Orionis open cluster (3+/-2 Ma, d = 360^+70_-60 pc). We have performed a deep IJ-band search, covering an area of 790 arcmin^2 close to the cluster centre. This survey was complemented with an infrared follow-up in the HKs- and Spitzer 3.6-8.0 mum-bands. Using colour-magnitude diagrams, we have selected 49 candidate cluster members in the magnitude interval 16.1 mag < I < 23.0 mag. Accounting for flux excesses at 8.0 mum and previously known spectral features of youth, 30 objects are bona fide cluster members. Four are first identified from our optical-near infrared data. Eleven have most probable masses below the deuterium burning limit and are classified as planetary-mass object candidates. The slope of the substellar mass spectrum (Delta N / Delta M = a M^-alpha) in the mass interval 0.11 Msol M < 0.006 Msol is alpha = +0.6+/-0.2. Any opacity mass-limit, if these objects form via fragmentation, may lie below 0.006 Msol. The frequency of sigma Orionis brown dwarfs with circumsubstellar discs is 47+/-15 %. The continuity in the mass function and in the frequency of discs suggests that very low-mass stars and substellar objects, even below the deuterium-burning mass limit, may share the same formation mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Optical spectroscopy of isolated planetary mass objects in the sigma Orionis cluster

We have obtained low resolution optical spectra of 15 isolated planetary mass objects (IPMOs) in the sigma Orionis cluster, and derived spectral types by comparison with nearby M and L dwarfs. The spectral types are in the range late M - mid L, in agreement with our expectations based on colors and magnitudes for bona fide members. Therefore, most of these objects have masses below the deuterium burning limit. About 2/3 show H(alpha) in emission at our spectral resolution. From our spectroscopic and photometric data, we infer that three IPMOs in this sample may be binaries with components of similar masses. These results confirm that the substellar mass function of the sigma Orionis cluster, in the form dN/dM, keeps rising in the planetary domain.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The Substellar Mass Function in sigma Orionis

We combine results from imaging searches for substellar objects in the sigma Orionis cluster and follow-up photometric and spectroscopic observations to derive a census of the brown dwarf population in a region of 847 arcmin^2. We identify 64 very low-mass cluster member candidates in this region. We have available three color (IZJ) photometry for all of them, spectra for 9 objects, and K photometry for 27% of our sample. These data provide a well defined sequence in the I vs I-J, I-K color magnitude diagrams, and indicate that the cluster is affected by little reddening despite its young age (~5 Myr). Using state-of-the-art evolutionary models, we derive a mass function from the low-mass stars (0.2 Msol) across the complete brown dwarf domain (0.075 Msol to 0.013 Msol), and into the realm of free-floating planetary-mass objects (<0.013 Msol). We find that the mass spectrum (dN/dm ~ m^{-alpha}) increases toward lower masses with an exponent alpha = 0.8+/-0.4. Our results suggest that planetary-mass isolated objects could be as common as brown dwarfs; both kinds of objects together would be as numerous as stars in the cluster. If the distribution of stellar and substellar masses in sigma Orionis is representative of the Galactic disk, older and much lower luminosity free-floating planetary-mass objects with masses down to about 0.005 Msol should be abundant in the solar vicinity, with a density similar to M-type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Variability in ultra cool dwarfs: evidence for the evolution of surface features

We present photometric light curves for a sample of 21 ultra cool M and L dwarfs in the field and in the young open clusters sigma Orionis and the Pleiades. Evidence for variability with RMS amplitudes (in the I band) of 0.01 to 0.055 magnitudes on timescales of 0.4 to 100 hours is discovered in half of these objects. Power spectral analysis using the CLEAN algorithm was performed to search for evidence of periodic variability. Some objects show strong periodicities at around a few hours, which could be due to rotational modulation of the light curve by surface features. However, several objects do not have any significant periodicities to explain their variability. The vsini values of a similar population of objects makes it very likely that our time sampling was sensitive to the expected range of rotation periods, and simulations show that we would have detected these if they were caused by long-lived surface features. We argue that this absence of periodicity is due to the evolution of the brightness, and presumably also the physical size, of surface features on timescales of a few to a few tens of hours. This is supported in the case of 2M1145 for which two light curves have been obtained one year apart and show no common periodicity. The surface features could plausibly be photospheric dust clouds or magnetically-induced spots. The recently observed decline in chromospheric activity for late type M and L dwarfs hints towards the former explanation for at least our later-type objects. Furthermore, our sample suggests that variability to be more common in objects later than M9, indicating that the variability may be related to dust formation. One light curve shows a brief, but significant, dip, which could be a short-lived feature or possibly an eclipse by a companion. (abridged)

0 citations0 reviewsNo code linkedOpen access