The unprecedented astrometry from Gaia DR2 provides us with an opportunity to study in detail molecular clouds in the solar neighbourhood. Extracting the wealth of information in these data remains a challenge, however. We have further improved our Gaussian Processes-based, three-dimensional dust mapping technique to allow us to study molecular clouds in more detail. These improvements include a significantly better scaling of the computational cost with the number of stars, and taking into account distance uncertainties to individual stars. Using Gaia DR2 astrometry together with 2MASS and WISE photometry for 30 000 stars, we infer the distribution of dust out to 600 pc in the direction of the Orion A molecular cloud. We identify a bubble-like structure in front of Orion A, centred at a distance of about 350 pc from the Sun. The main Orion A structure is visible at slightly larger distances, and we clearly see a tail extending over 100 pc that is curved and slightly inclined to the line-of-sight. The location of our foreground structure coincides with 5-10 Myr old stellar populations, suggesting a star formation episode that predates that of the Orion Nebula Cluster itself. We identify also the main structure of the Orion B molecular cloud, and in addition discover a background component to this at a distance of about 460 pc from the Sun. Finally, we associate our dust components at different distances with the plane-of-the-sky magnetic field orientation as mapped by Planck. This provides valuable information for modelling the magnetic field in 3D around star forming regions.

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.

We present a mock stellar catalog, matching in volume, depth and data model the content of the planned Gaia early data release 3 (Gaia EDR3). We have generated our catalog (GeDR3mock) using galaxia, a tool to sample stars from an underlying Milky Way (MW) model or from N-body data. We used an updated Besan\c{c}on Galactic model together with the latest PARSEC stellar evolutionary tracks, now also including white dwarfs. We added the Magellanic clouds and realistic open clusters with internal rotation. We empirically modelled uncertainties based on Gaia DR2 (GDR2) and scaled them according to the longer baseline in Gaia EDR3. The apparent magnitudes were reddened according to a new selection of 3D extinction maps. To help with the Gaia selection function we provide all-sky magnitude limit maps in G and BP for a few relevant GDR2 subsets together with the routines to produce these maps for user-defined subsets. We supplement the catalog with photometry and extinctions in non-Gaia bands. The catalog is available in the Virtual Observatory and can be queried just like the actual Gaia EDR3 will be. We highlight a few capabilities of the Astronomy Data Query Language (ADQL) with educative catalog queries. We use the data extracted from those queries to compare GeDR3mock to GDR2, which emphasises the importance of adding observational noise to the mock data. Since the underlying truth, e.g. stellar parameters, is know in GeDR3mock, it can be used to construct priors as well as mock data tests for parameter estimation. All code, models and data used to produce GeDR3mock are linked and contained in galaxia_wrap, a python package, representing a fast galactic forward model, able to project MW models and N-body data into realistic Gaia observables.

The discovery of the second interstellar object 2I/Borisov on 2019 August 30 raises the question of whether it was ejected recently from a nearby stellar system. Here we compute the asymptotic incoming trajectory of 2I/Borisov, based on both recent and pre-discovery data extending back to December 2018, using a range of force models that account for cometary outgassing. From Gaia DR2 astrometry and radial velocities, we trace back in time the Galactic orbits of 7.4 million stars to look for close encounters with 2I/Borisov. The closest encounter we find took place 910kyr ago with the M0V star Ross 573, at a separation of 0.068pc (90% confidence interval of 0.053-0.09pc) with a relative velocity of 23km/s. This encounter is nine times closer than the closest past encounter identified for the first interstellar object 1I/'Oumuamua. Ejection of 2I/Borisov via a three-body encounter in a binary or planetary system is possible, although such a large ejection velocity is unlikely to be obtained and Ross 573 shows no signs of binarity. We also identify and discuss some other recent close encounters, recognizing that if 2I/Borisov is more than about 10Myr old, our search would be unlikely to find its parent system.

The Pioneer 10 and 11 and Voyager 1 and 2 spacecraft, launched in the 1970s, are heading out of the solar system. Using the astrometric and radial velocity data from the second Gaia data release, we integrate the trajectories of 7.4 million stars, and the spacecraft, through a Galactic potential in order to identify those stars the spacecraft will pass closest to. The closest encounters for all spacecraft take place at separations between 0.2 and 0.5 pc within the next million years. The closest encounter will be by Pioneer 10 with the K8 dwarf HIP 117795, at 0.23 pc in 90 kyr at a high relative velocity of 291 km/s.

We construct a supervised classifier based on Gaussian Mixture Models to probabilistically classify objects in Gaia data release 2 (GDR2) using only photometric and astrometric data in that release. The model is trained empirically to classify objects into three classes -- star, quasar, galaxy -- for G<=14.5 mag down to the Gaia magnitude limit of G=21.0 mag. Galaxies and quasars are identified for the training set by a cross-match to objects with spectroscopic classifications from the Sloan Digital Sky Survey. Stars are defined directly from GDR2. When allowing for the expectation that quasars are 500 times rarer than stars, and galaxies 7500 times rarer than stars (the class imbalance problem), samples classified with a threshold probability of 0.5 are predicted to have purities of 0.43 for quasars and 0.28 for galaxies, and completenesses of 0.58 and 0.72 respectively. The purities can be increased up to 0.60 by adopting a higher threshold. Not accounting for this expected low frequency of extragalactic objects (the class prior) would give both erroneously optimistic performance predictions and severely impure samples. Applying our model to all 1.20 billion objects in GDR2 with the required features, we classify 2.3 million objects as quasars and 0.37 million objects as galaxies (with individual probabilities above 0.5). The small number of galaxies is due to the strong bias of the satellite detection algorithm and on-ground data selection against extended objects. We infer the true number of quasars and galaxies -- as these classes are defined by our training set -- to be 690,000 and 110,000 respectively (+/- 50%). The aim of this work is to see how well extragalactic objects can be classified using only GDR2 data. Better classifications should be possible with the low resolution spectroscopy (BP/RP) planned for GDR3.

The first detected interstellar object 'Oumuamua that passed within 0.25au of the Sun on 2017 September 9 was presumably ejected from a stellar system. We use its newly determined non-Keplerian trajectory together with the reconstructed Galactic orbits of 7 million stars from Gaia DR2 to identify past close encounters. Such an "encounter" could reveal the home system from which 'Oumuamua was ejected. The closest encounter, at 0.60pc (0.53-0.67pc, 90% confidence interval), was with the M2.5 dwarf HIP 3757 at a relative velocity of 24.7km/s, 1Myr ago. A more distant encounter (1.6pc) but with a lower encounter (ejection) velocity of 10.7km/s was with the G5 dwarf HD 292249, 3.8Myr ago. Two more stars have encounter distances and velocities intermediate to these. The encounter parameters are similar across six different non-gravitational trajectories for 'Oumuamua. Ejection of 'Oumuamua by scattering from a giant planet in one of the systems is plausible, but requires a rather unlikely configuration to achieve the high velocities found. A binary star system is more likely to produce the observed velocities. None of the four home candidates have published exoplanets or are known to be binaries. Given that the 7 million stars in Gaia DR2 with 6D phase space information is just a small fraction of all stars for which we can eventually reconstruct orbits, it is a priori unlikely that our current search would find 'Oumuamua's home star system. As 'Oumuamua is expected to pass within 1pc of about 20 stars and brown dwarfs every Myr, the plausibility of a home system depends also on an appropriate (low) encounter velocity.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

For the vast majority of stars in the second Gaia data release, reliable distances cannot be obtained by inverting the parallax. A correct inference procedure must instead be used to account for the nonlinearity of the transformation and the asymmetry of the resulting probability distribution. Here we infer distances to essentially all 1.33 billion stars with parallaxes published in the second \gaia\ data release. This is done using a weak distance prior that varies smoothly as a function of Galactic longitude and latitude according to a Galaxy model. The irreducible uncertainty in the distance estimate is characterized by the lower and upper bounds of an asymmetric confidence interval. Although more precise distances can be estimated for a subset of the stars using additional data (such as photometry), our goal is to provide purely geometric distance estimates, independent of assumptions about the physical properties of, or interstellar extinction towards, individual stars. We analyse the characteristics of the catalogue and validate it using clusters. The catalogue can be queried on the Gaia archive using ADQL at http://gea.esac.esa.int/archive/ and downloaded from http://www.mpia.de/~calj/gdr2_distances.html .

Large stellar surveys are sensitive to interstellar dust through the effects of reddening. Using extinctions measured from photometry and spectroscopy, together with three-dimensional (3D) positions of individual stars, it is possible to construct a three-dimensional dust map. We present the first continuous map of the dust distribution in the Galactic disk out to 7 kpc within 100 pc of the Galactic midplane, using red clump and giant stars from SDSS APOGEE DR14. We use a non-parametric method based on Gaussian Processes to map the dust density, which is the local property of the ISM rather than an integrated quantity. This method models the dust correlation between points in 3D space and can capture arbitrary variations, unconstrained by a pre-specified functional form. This produces a continuous map without line-of-sight artefacts. Our resulting map traces some features of the local Galactic spiral arms, even though the model contains no prior suggestion of spiral arms, nor any underlying model for the Galactic structure. This is the first time that such evident arm structures have been captured by a dust density map in the Milky Way. Our resulting map also traces some of the known giant molecular clouds in the Galaxy and puts some constraints on their distances, some of which were hitherto relatively uncertain.

Passing stars may play an important role in the evolution of our solar system. We search for close stellar encounters to the Sun among all 7.2 million stars in Gaia-DR2 that have six-dimensional phase space data. We characterize encounters by integrating their orbits through a Galactic potential and propagating the correlated uncertainties via a Monte Carlo resampling. After filtering to remove spurious data, we find 694 stars that have median (over uncertainties) closest encounter distances within 5 pc, all occurring within 15 Myr from now. 26 of these have at least a 50% chance of coming closer than 1 pc (and 7 within 0.5 pc), all but one of which are newly discovered here. We further confirm some and refute several other previously-identified encounters, confirming suspicions about their data. The closest encounter in the sample is Gl 710, which has a 95% probability of coming closer than 0.08 pc (17 000 AU). Taking mass estimates from Gaia astrometry and multiband photometry for essentially all encounters, we find that Gl 710 also has the largest impulse on the Oort cloud. Using a Galaxy model, we compute the completeness of the Gaia-DR2 encountering sample as a function of perihelion time and distance. Only 15% of encounters within 5 pc occurring within +/- 5 Myr of now have been identified, mostly due to the lack of radial velocities for faint and/or cool stars. Accounting for the incompleteness, we infer the present rate of encounters within 1 pc to be 19.7 +/- 2.2 per Myr, a quantity expected to scale quadratically with the encounter distance out to at least several pc. Spuriously large parallaxes in our sample from imperfect filtering would tend to inflate both the number of encounters found and this inferred rate. The magnitude of this effect is hard to quantify.

The second Gaia data release (GDR2) provides precise five-parameter astrometric data (positions, proper motions and parallaxes) for an unprecedented amount of sources (more than $1.3$ billion, mostly stars). The use of this wealth of astrometric data comes with a specific challenge: how does one properly infer from these data the astrophysical parameters of interest? The main - but not only - focus of this paper is the issue of the estimation of distances from parallaxes, possibly combined with other information. We start with a critical review of the methods traditionally used to obtain distances from parallaxes and their shortcomings. Then we provide guidelines on how to use parallaxes more efficiently to estimate distances by using Bayesian methods. In particular also we show that negative parallaxes, or parallaxes with relatively larger uncertainties still contain valuable information. Finally, we provide examples that show more generally how to use astrometric data for parameter estimation, including the combination of proper motions and parallaxes and the handling of covariances in the uncertainties. The paper contains examples based on simulated Gaia data to illustrate the problems and the solutions proposed. Furthermore, the developments and methods proposed in the paper are linked to a set of tutorials included in the Gaia archive documentation that provide practical examples and a good starting point for the application of the recommendations to actual problems. In all cases the source code for the analysis methods is provided. Our main recommendation is to always treat the derivation of (astro-) physical parameters from astrometric data, in particular when parallaxes are involved, as an inference problem which should preferably be handled with a full Bayesian approach.

We present the second Gaia data release, Gaia DR2, consisting of astrometry, photometry, radial velocities, and information on astrophysical parameters and variability, for sources brighter than magnitude 21. In addition epoch astrometry and photometry are provided for a modest sample of minor planets in the solar system. A summary of the contents of Gaia DR2 is presented, accompanied by a discussion on the differences with respect to Gaia DR1 and an overview of the main limitations which are still present in the survey. Recommendations are made on the responsible use of Gaia DR2 results. Gaia DR2 contains celestial positions and the apparent brightness in G for approximately 1.7 billion sources. For 1.3 billion of those sources, parallaxes and proper motions are in addition available. The sample of sources for which variability information is provided is expanded to 0.5 million stars. This data release contains four new elements: broad-band colour information in the form of the apparent brightness in the $G_\mathrm{BP}$ (330--680 nm) and $G_\mathrm{RP}$ (630--1050 nm) bands is available for 1.4 billion sources; median radial velocities for some 7 million sources are presented; for between 77 and 161 million sources estimates are provided of the stellar effective temperature, extinction, reddening, and radius and luminosity; and for a pre-selected list of 14000 minor planets in the solar system epoch astrometry and photometry are presented. Finally, Gaia DR2 also represents a new materialisation of the celestial reference frame in the optical, the Gaia-CRF2, which is the first optical reference frame based solely on extragalactic sources. There are notable changes in the photometric system and the catalogue source list with respect to Gaia DR1, and we stress the need to consider the two data releases as independent.

The second Gaia data release (Gaia-DR2) contains, beyond the astrometry, three-band photometry for 1.38 billion sources. We have used these three broad bands to infer stellar effective temperatures, Teff, for all sources brighter than G=17 mag with Teff in the range 3000-10 000 K (161 million sources). Using in addition the parallaxes, we infer the line-of-sight extinction, A_G, and the reddening, E[BP-RP], for 88 million sources. Together with a bolometric correction we derive luminosity and radius for 77 million sources. These quantities as well as their estimated uncertainties are part of Gaia-DR2. Here we describe the procedures by which these quantities were obtained, including the underlying assumptions, comparison with literature estimates, and the limitations of our results. Typical accuracies are of order 324 K (Teff), 0.46 mag (A_G), 0.23 mag (E[BP-RP]), 15% (luminosity), and 10% (radius). Being based on only a small number of observable quantities and limited training data, our results are necessarily subject to some extreme assumptions that can lead to strong systematics in some cases (not included in the aforementioned accuracy estimates). One aspect is the non-negativity contraint of our estimates, in particular extinction. Yet in several regions of parameter space our results show very good performance, for example for red clump stars and solar analogues. Large uncertainties render the extinctions less useful at the individual star level, but they show good performance for ensemble estimates. We identify regimes in which our parameters should and should not be used and we define a "clean" sample. Despite the limitations, this is the largest catalogue of uniformly-inferred stellar parameters to date. More precise and detailed astrophysical parameters based on the full BP/RP spectrophotometry are planned as part of the third Gaia data release.

We present a mock catalog of Milky Way stars, matching in volume and depth the content of the Gaia data release 2 (GDR2). We generated our catalog using Galaxia, a tool to sample stars from a Besancon Galactic model, together with a realistic 3D dust extinction map. The catalog mimicks the complete GDR2 data model and contains most of the entries in the Gaia source catalog: 5-parameter astrometry, 3-band photometry, radial velocities, stellar parameters, and associated scaled nominal uncertainty estimates. In addition, we supplemented the catalog with extinctions and photometry for non-Gaia bands. This catalog can be used to prepare GDR2 queries in a realistic runtime environment, and it can serve as a Galactic model against which to compare the actual GDR2 data in the space of observables. The catalog is hosted through the virtual observatory GAVO's Heidelberg data center service and thus can be queried using ADQL as for GDR2 data.

We present a map of the three-dimensional (3D) distribution of dust in the Orion complex. Orion is the closest site of high-mass star formation, making it an excellent laboratory for studying the interstellar medium and star formation. We use data from the Gaia-TGAS catalogue combined with photometry from 2MASS and WISE to get the distances and extinctions of individual stars in the vicinity of the Orion complex. We adopt the non-parametric method of Rezaei Kh. et al. (2017) to infer the probability distribution function of the dust densities at arbitrary points throughout the region. We map the dust distribution towards different parts of the Orion complex and find the distance and depth of the cloud compatible with other recent works which show the ability of the method to be applicable on the local molecular clouds to map their 3D dust distribution. We also demonstrate the danger of only using colours of stars to derive their extinctions without considering further physical constraints like the colour-magnitude diagram (CMD).

I report on close encounters of stars to the Sun found in the first Gaia data release (GDR1). Combining Gaia astrometry with radial velocities of around 320 000 stars drawn from various catalogues, I integrate orbits in a Galactic potential to identify those stars which pass within a few parsecs. Such encounters could influence the solar system, for example through gravitational perturbations of the Oort cloud. 16 stars are found to come within 2 pc (although a few of these have dubious data). This is fewer than were found in a similar study based on Hipparcos data, even though the present study has many more candidates. This is partly because I reject stars with large radial velocity uncertainties (>10 km/s), and partly because of missing stars in GDR1 (especially at the bright end). The closest encounter found is Gl 710, a K dwarf long-known to come close to the Sun in about 1.3 Myr. The Gaia astrometry predict a much closer passage than pre-Gaia estimates, however: just 16 000 AU (90% confidence interval: 10 000 to 21 000 AU), which will bring this star well within the Oort cloud. Using a simple model for the spatial, velocity, and luminosity distributions of stars, together with an approximation of the observational selection function, I model the incompleteness of this Gaia-based search as a function of the time and distance of closest approach. Applying this to a subset of the observed encounters (excluding duplicates and stars with implausibly large velocities), I estimate the rate of stellar encounters within 5 pc averaged over the past and future 5 Myr to be 545 +/- 59 per Myr. Assuming a quadratic scaling of the rate within some encounter distance (which my model predicts), this corresponds to 87 +/- 9 per My within 2 pc. A more accurate analysis and assessment will be possible with future Gaia data releases.

We present H$\alpha$ spectroscopic observations and detailed modelling of the Balmer filaments in the supernova remnant Tycho. We used Galaxy H$\alpha$ Fabry-P\'erot Spectrometer on the William Herschel Telescope with a 3.4'$\times$3.4' field-of-view, 0.2" pixel scale and $\sigma_\rm{instr}=8.1$ km/s resolution at 1" seeing for $\sim10$ hours, resulting in 82 spatial-spectral bins that resolve the narrow H$\alpha$ line in the entire Tycho's northeastern rim. For the first time, we can mitigate artificial line broadening from unresolved differential motion, and probe H$\alpha$ emission parameters in varying shock and ambient medium conditions. Broad H$\alpha$ line remains unresolved within spectral coverage of 392 km/s. We employed Bayesian inference to obtain reliable parameter confidence intervals, and quantify the evidence for models with multiple line components. The median H$\alpha$ narrow-line full-width at half-maximum of all bins and models is $W_\rm{NL}=(54.8\pm1.8)$ km/s at the $95\%$ confidence level, varying within [35, 72] km/s between bins and clearly broadened compared to the intrinsic (thermal) $\approx20$ km/s. Possible line splits are accounted for, significant in $\approx18\%$ of the filament, and presumably due to remaining projection effects. We also find wide-spread evidence for intermediate-line emission of a broad-neutral precursor, with median $W_\rm{IL}=(180\pm14)$ km/s ($95\%$ confidence). Finally, we present a measurement of the remnant's systemic velocity, $V_\rm{LSR}=-34$ km/s, and map differential line-of-sight motions. Our results confirm the existence and interplay of shock precursors in Tycho's remnant. In particular, we show that suprathermal narrow-line emission is near-universal in Tycho and that, in absence of an alternative explanation, collisionless supernova remnant shocks constitute a viable acceleration source for Galactic TeV Cosmic-Ray protons.

We present wide-field, spatially and highly resolved spectroscopic observations of Balmer filaments in the northeastern rim of Tycho's supernova remnant in order to investigate the signal of cosmic-ray (CR) acceleration. The spectra of Balmer-dominated shocks (BDSs) have characteristic narrow (FWHM $\sim$ 10 kms$^{-1}$) and broad (FWHM $\sim$ 1000 kms$^{-1}$) H$\alpha$ components. CRs affect the H$\alpha$-line parameters: heating the cold neutrals in the interstellar medium results in broadening of the narrow H$\alpha$-line width beyond 20 kms$^{-1}$, but also in reduction of the broad H$\alpha$-line width due to energy being removed from the protons in the post-shock region. For the first time we show that the width of the narrow H$\alpha$ line, much larger than 20 kms$^{-1}$, is not a resolution or geometric effect nor a spurious result of a neglected intermediate (FWHM $\sim$ 100 kms$^{-1}$) component resulting from hydrogen atoms undergoing charge exchange with warm protons in the broad-neutral precursor. Moreover, we show that a narrow line width $\gg$ 20 kms$^{-1}$ extends across the entire NE rim, implying CR acceleration is ubiquitous, and making it possible to relate its strength to locally varying shock conditions. Finally, we find several locations along the rim, where spectra are significantly better explained (based on Bayesian evidence) by inclusion of the intermediate component, with a width of 180 kms$^{-1}$ on average.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Pulsating stars, such as Cepheids, Miras, and RR Lyrae stars, are important distance indicators and calibrators of the "cosmic distance ladder", and yet their period-luminosity-metallicity (PLZ) relations are still constrained using simple statistical methods that cannot take full advantage of available data. To enable optimal usage of data provided by the Gaia mission, we present a probabilistic approach that simultaneously constrains parameters of PLZ relations and uncertainties in Gaia parallax measurements. We demonstrate this approach by constraining PLZ relations of type $ab$ RR Lyrae stars in near-infrared W1 and W2 bands, using Tycho-Gaia Astrometric Solution (TGAS) parallax measurements for a sample of $\approx100$ type $ab$ RR Lyrae stars located within 2.5 kpc of the Sun. The fitted PLZ relations are consistent with previous studies, and in combination with other data, deliver distances precise to 6% (once various sources of uncertainty are taken into account). To a precision of 0.05 mas ($1\sigma$), we do not find a statistically significant offset in TGAS parallaxes for this sample of distant RR Lyrae stars (median parallax of 0.8 mas and distance of 1.4 kpc). With only minor modifications, our probabilistic approach can be used to constrain PLZ relations of other pulsating stars, and we intend to apply it to Cepheid and Mira stars in the near future.

We present a non-parametric model for inferring the three-dimensional (3D) distribution of dust density in the Milky Way. Our approach uses the extinction measured towards stars at different locations in the Galaxy at approximately known distances. Each extinction measurement is proportional to the integrated dust density along its line-of-sight. Making simple assumptions about the spatial correlation of the dust density, we can infer the most probable 3D distribution of dust across the entire observed region, including along sight lines which were not observed. This is possible because our model employs a Gaussian Process to connect all lines-of-sight. We demonstrate the capability of our model to capture detailed dust density variations using mock data as well as simulated data from the Gaia Universe Model Snapshot. We then apply our method to a sample of giant stars observed by APOGEE and Kepler to construct a 3D dust map over a small region of the Galaxy. Due to our smoothness constraint and its isotropy, we provide one of the first maps which does not show the "fingers of god" effect.

We infer distances and their asymmetric uncertainties for two million stars using the parallaxes published in the Gaia DR1 (GDR1) catalogue. We do this with two distance priors: A minimalist, isotropic prior assuming an exponentially decreasing space density with increasing distance, and an anisotropic prior derived from the observability of stars in a Milky Way model. We validate our results by comparing our distance estimates for 105 Cepheids which have more precise, independently estimated distances. For this sample we find that the Milky Way prior performs better (the RMS of the scaled residuals is 0.40) than the exponentially decreasing space density prior (RMS is 0.57), although for distances beyond 2 kpc the Milky Way prior performs worse, with a bias in the scaled residuals of -0.36 (vs. -0.07 for the exponentially decreasing space density prior). We do not attempt to include the photometric data in GDR1 due to the lack of reliable colour information. Our distance catalogue is available at http://www.mpia.de/homes/calj/tgas distances/main.html as well as at CDS. This should only be used to give individual distances. Combining data or testing models should be done with the original parallaxes, and attention paid to correlated and systematic uncertainties.

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.

Estimating a distance by inverting a parallax is only valid in the absence of noise. As most stars in the Gaia catalogue will have non-negligible fractional parallax errors, we must treat distance estimation as a constrained inference problem. Here we investigate the performance of various priors for estimating distances, using a simulated Gaia catalogue of one billion stars. We use three minimalist, isotropic priors, as well an anisotropic prior derived from the observability of stars in a Milky Way model. The two priors that assume a uniform distribution of stars--either in distance or in space density---give poor results: The root mean square fractional distance error, f_RMS, grows far in excess of 100% once the fractional parallax error, f_true, is larger than 0.1. A prior assuming an exponentially decreasing space density with increasing distance performs well once its single scale length parameter has been set to an appropriate value: f_RMS is roughly equal to f_true for f_true < 0.4, yet does not increase further as f_true increases up to to 1.0. The Milky Way prior performs well except towards the Galactic centre, due to a mismatch with the (simulated) data. Such mismatches will be inevitable (and remain unknown) in real applications, and can produce large errors. We therefore suggest to adopt the simpler exponentially decreasing space density prior, which is also less time-consuming to compute. Including Gaia photometry improves the distance estimation significantly for both the Milky Way and exponentially decreasing space density prior, yet doing so requires additional assumptions about the physical nature of stars.

We present a machine learning package for the classification of periodic variable stars. Our package is intended to be general: it can classify any single band optical light curve comprising at least a few tens of observations covering durations from weeks to years, with arbitrary time sampling. We use light curves of periodic variable stars taken from OGLE and EROS-2 to train the model. To make our classifier relatively survey-independent, it is trained on 16 features extracted from the light curves (e.g. period, skewness, Fourier amplitude ratio). The model classifies light curves into one of seven superclasses - Delta Scuti, RR Lyrae, Cepheid, Type II Cepheid, eclipsing binary, long-period variable, non-variable - as well as subclasses of these, such as ab, c, d, and e types for RR Lyraes. When trained to give only superclasses, our model achieves 0.98 for both recall and precision as measured on an independent validation dataset (on a scale of 0 to 1). When trained to give subclasses, it achieves 0.81 for both recall and precision. In order to assess classification performance of the subclass model, we applied it to the MACHO, LINEAR, and ASAS periodic variables, which gave recall/precision of 0.92/0.98, 0.89/0.96, and 0.84/0.88, respectively. We also applied the subclass model to Hipparcos periodic variable stars of many other variability types that do not exist in our training set, in order to examine how much those types degrade the classification performance of our target classes. In addition, we investigate how the performance varies with the number of data points and duration of observations. We find that recall and precision do not vary significantly if the number of data points is larger than 80 and the duration is more than a few weeks. The classifier software of the subclass model is available from the GitHub repository (https://goo.gl/xmFO6Q).

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.

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.

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.

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.

(Abridged) We aim at identifying the least massive population of the solar metallicity, young (120 Myr), nearby (133.5 pc) Pleiades star cluster with the ultimate goal of understanding the physical properties of intermediate-age, free-floating, low-mass brown dwarfs and giant planetary-mass objects, and deriving the cluster substellar mass function across the deuterium-burning mass limit at ~0.012 Msol. We performed a deep photometric and astrometric J- and H-band survey covering an area of ~0.8 deg^2. The images with completeness and limiting magnitudes of J,H ~ 20.2 and ~ 21.5 mag were acquired ~9 yr apart (proper motion precision of +/-6 mas/yr). J- and H-band data were complemented with Z, K, and mid-infrared magnitudes up to 4.6 micron coming from UKIDSS, WISE, and follow-up observations of our own. Pleiades member candidates were selected to have proper motions compatible with that of the cluster, and colors following the known Pleiades sequence in the interval J = 15.5-8.8 mag, and Z_UKIDSS - J > 2.3 mag or Z nondetections for J > 18.8 mag. We found a neat sequence of astrometric and photometric Pleiades substellar member candidates in the intervals J = 15.5-21.2 mag and ~0.072-0.008 Msol. The faintest objects show very red near- and mid-infrared colors exceeding those of field high-gravity dwarfs by >0.5 mag. The Pleiades photometric sequence does not show any color turn-over because of the presence of photospheric methane absorption down to J = 20.3 mag, which is about 1 mag fainter than predicted by the color-computed models. Pleiades brown dwarfs have a proper motion dispersion of 6.4-7.5 mas/yr and are dynamically relaxed at the age of the cluster. The Pleiades mass function extends down to the deuterium burning-mass threshold, with a slope fairly similar to that of other young star clusters and stellar associations.

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.

The EPOCH (EROS-2 periodic variable star classification using machine learning) project aims to detect periodic variable stars in the EROS-2 light curve database. In this paper, we present the first result of the classification of periodic variable stars in the EROS-2 LMC database. To classify these variables, we first built a training set by compiling known variables in the Large Magellanic Cloud area from the OGLE and MACHO surveys. We crossmatched these variables with the EROS-2 sources and extracted 22 variability features from 28 392 light curves of the corresponding EROS-2 sources. We then used the random forest method to classify the EROS-2 sources in the training set. We designed the model to separate not only $\delta$ Scuti stars, RR Lyraes, Cepheids, eclipsing binaries, and long-period variables, the superclasses, but also their subclasses, such as RRab, RRc, RRd, and RRe for RR Lyraes, and similarly for the other variable types. The model trained using only the superclasses shows 99% recall and precision, while the model trained on all subclasses shows 87% recall and precision. We applied the trained model to the entire EROS-2 LMC database, which contains about 29 million sources, and found 117 234 periodic variable candidates. Out of these 117 234 periodic variables, 55 285 have not been discovered by either OGLE or MACHO variability studies. This set comprises 1 906 $\delta$ Scuti stars, 6 607 RR Lyraes, 638 Cepheids, 178 Type II Cepheids, 34 562 eclipsing binaries, and 11 394 long-period variables. A catalog of these EROS-2 LMC periodic variable stars will be available online at http://stardb.yonsei.ac.kr and at the CDS website (http://vizier.u-strasbg.fr/viz-bin/VizieR).

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 $\mu$ 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}$, $\mu$ 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.

We present two epochs of MPG/ESO 2.2m GROND simultaneous 6-band ($r'i'z'JHK$) photometric monitoring of the closest known L/T transition brown dwarf binary WISE J104915.57-531906.1AB. We report here the first resolved variability monitoring of both the T0.5 and L7.5 components. We obtained 4 hours of focused observations on the night of UT 2013-04-22, as well as 4 hours of defocused (unresolved) observations on the night of UT 2013-04-16. We note a number of robust trends in our light curves. The $r'$ and $i'$ light curves appear to be anticorrelated with $z'$ and $H$ for the T0.5 component and in the unresolved lightcurve. In the defocused dataset, $J$ appears correlated with $z'$ and $H$ and anticorrelated with $r'$ and $i'$, while in the focused dataset we measure no variability for $J$ at the level of our photometric precision, likely due to evolving weather phenomena. In our focused T0.5 component lightcurve, the $K$ band lightcurve displays a significant phase offset relative to both $H$ and $z'$. We argue that the measured phase offsets are correlated with atmospheric pressure probed at each band, as estimated from 1D atmospheric models. We also report low-amplitude variability in $i'$ and $z'$ intrinsic to the L7.5 component.

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.

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.

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.

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.

The optical light curves of many quasars show variations of tenths of a magnitude or more on time scales of months to years. This variation often cannot be described well by a simple deterministic model. We perform a Bayesian comparison of over 20 deterministic and stochastic models on 6304 QSO light curves in SDSS Stripe 82. We include the damped random walk (or Ornstein-Uhlenbeck [OU] process), a particular type of stochastic model which recent studies have focused on. Further models we consider are single and double sinusoids, multiple OU processes, higher order continuous autoregressive processes, and composite models. We find that only 29 out of 6304 QSO lightcurves are described significantly better by a deterministic model than a stochastic one. The OU process is an adequate description of the vast majority of cases (6023). Indeed, the OU process is the best single model for 3462 light curves, with the composite OU process/sinusoid model being the best in 1706 cases. The latter model is the dominant one for brighter/bluer QSOs. Furthermore, a non-negligible fraction of QSO lightcurves show evidence that not only the mean is stochastic but the variance is stochastic, too. Our results confirm earlier work that QSO light curves can be described with a stochastic model, but place this on a firmer footing, and further show that the OU process is preferred over several other stochastic and deterministic models. Of course, there may well exist yet better (deterministic or stochastic) models which have not been considered here.

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.

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.

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.

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.

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)

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.

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.

Studies of the mass function (MF) of open clusters of different ages allow us to probe the efficiency with which brown dwarfs evaporate from clusters to populate the field. Surveys of older clusters (age >100 Myr) are not affected so severely by several problems encountered in young clusters, such as intra-cluster extinction and large uncertainties in brown dwarf models. We present the results of a deep photometric survey to study the MF of the central region of the old open cluster Praesepe (age~590$^{+150}_{-120}$ Myr, distance ~190$^{+6.0}_{-5.8}$ pc), down to the substellar regime. We performed an optical ($riz$ and $Y$-band) photometric survey of Praesepe using the Large Binocular Telescope Camera covering an area of ~0.59 deg$^2$ in the cluster centre from $i$~19.0 mag (~100 M_Jup) down to a 5$\sigma$ detection limit at $i$~25.6mag (~40 M_Jup). The survey is approximately 95% complete at $i=23.8$ mag and $z=22.0$ mag (~55 M_Jup). We identify 59 cluster member candidates, of which 37 are substellar, by comparing with the predictions of a dusty atmosphere model. The MF of those candidates rises from the substellar boundary until ~67 M_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 clusters, the MFs of which peak at ~10 M_Jup. Either Praesepe really does have a different MF from other clusters or they had similar initial MFs but a different dynamical evolution. Since most of the candidates are faint, we lack astrometric or spectroscopic follow-ups to test their memberships. However, the contaminations by field dwarfs, galaxies, or giants are found to have little effect on the shape of MF and therefore the MF of `real' cluster members should have similar characteristics.

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.

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.

An efficient separation between dwarfs and giants in surveys of bright stars is important, especially for studies in which distances are estimated through photometric parallax relations. We use the available spectroscopic log g estimates from the second RAVE data release (DR2) to assign each star a probability for being a dwarf or subgiant/giant based on mixture model fits to the log g distribution in different color bins. We further attempt to use these stars as a labeled training set in order to classify stars which lack log g estimates into dwarfs and giants with a SVM algorithm. We assess the performance of this classification against different choices of the input feature vector. In particular, we use different combinations of reduced proper motions, 2MASS JHK, DENIS IJK and USNO-B B2R2 apparent magnitudes. Our study shows that -- for our color ranges -- the infrared bands alone provide no relevant information to separate dwarfs and giants. Even when optical bands and reduced proper motions are added, the fraction of true giants classified as dwarfs (the contamination) remains above 20%. Using only the dwarfs with available spectroscopic log g and distance estimates (the latter from Breddels et al. 2010), we then repeat the stream search by Klement, Fuchs & Rix (2008, KFR08), which assumed all stars were dwarfs and claimed the discovery of a new stellar stream at V = -160 km/s in a sample of 7015 stars from RAVE DR1. Our re-analysis of the pure DR2 dwarf sample exhibits an overdensity of 5 stars at the phase-space position of the KFR08 stream, with a metallicity distribution that appears inconsistent with that of stars at comparably low rotational velocities. Compared to several smooth Milky Way models, the mean standardized deviation of the KFR08 stream is only marginal at 1.6$\pm$0.4... (abbreviated)

We investigate the performance of some common machine learning techniques in identifying BHB stars from photometric data. To train the machine learning algorithms, we use previously published spectroscopic identifications of BHB stars from SDSS data. We investigate the performance of three different techniques, namely k nearest neighbour classification, kernel density estimation and a support vector machine (SVM). We discuss the performance of the methods in terms of both completeness and contamination. We discuss the prospect of trading off these values, achieving lower contamination at the expense of lower completeness, by adjusting probability thresholds for the classification. We also discuss the role of prior probabilities in the classification performance, and we assess via simulations the reliability of the dataset used for training. Overall it seems that no-prior gives the best completeness, but adopting a prior lowers the contamination. We find that the SVM generally delivers the lowest contamination for a given level of completeness, and so is our method of choice. Finally, we classify a large sample of SDSS DR7 photometry using the SVM trained on the spectroscopic sample. We identify 27,074 probable BHB stars out of a sample of 294,652 stars. We derive photometric parallaxes and demonstrate that our results are reasonable by comparing to known distances for a selection of globular clusters. We attach our classifications, including probabilities, as an electronic table, so that they can be used either directly as a BHB star catalogue, or as priors to a spectroscopic or other classification method. We also provide our final models so that they can be directly applied to new data.

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)

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

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.

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.

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.

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)

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.

We present some of the strategies being developed to classify and parameterize objects obtained with spectra from the Sloan Digital Sky Survey (SDSS) and the RAdial Velocity Experiment (RAVE) and present some results. We estimate stellar atmospheric parameters (effective temperature, gravity, and metallicity) from spectral and photometric data and use these to analyse Galactic populations. We demonstrate this through the selection of a sample of candidate Blue Horizontal-Branch and RR Lyrae stars selected from SDSS/SEGUE.

Using the Wide Field Imager (WFI) at the ESO 2.2m telescope at La Silla and the CPAPIR camera at the CTIO 1.5m telescope at Cerro Tololo, we have performed an extensive, multiband photometric survey of the open cluster IC2391 (D~146pc, age~50Myr, solar metallicity). Here we present the results from our photometric survey and from a spectroscopic follow-up of the central part of the survey.

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)

Using effective temperature and metallicity derived from SDSS spectra for ~60,000 F and G type main sequence stars (0.2<g-r<0.6), we develop polynomial models for estimating these parameters from the SDSS u-g and g-r colors. We apply this method to SDSS photometric data for about 2 million F/G stars and measure the unbiased metallicity distribution for a complete volume-limited sample of stars at distances between 500 pc and 8 kpc. The metallicity distribution can be exquisitely modeled using two components with a spatially varying number ratio, that correspond to disk and halo. The two components also possess the kinematics expected for disk and halo stars. The metallicity of the halo component is spatially invariant, while the median disk metallicity smoothly decreases with distance from the Galactic plane from -0.6 at 500 pc to -0.8 beyond several kpc. The absence of a correlation between metallicity and kinematics for disk stars is in a conflict with the traditional decomposition in terms of thin and thick disks. We detect coherent substructures in the kinematics--metallicity space, such as the Monoceros stream, which rotates faster than the LSR, and has a median metallicity of [Fe/H]=-0.96, with an rms scatter of only ~0.15 dex. We extrapolate our results to the performance expected from the Large Synoptic Survey Telescope (LSST) and estimate that the LSST will obtain metallicity measurements accurate to 0.2 dex or better, with proper motion measurements accurate to ~0.2 mas/yr, for about 200 million F/G dwarf stars within a distance limit of ~100 kpc (g<23.5). [abridged]

The halo of the Milky Way provides unique elemental abundance and kinematic information on the first objects to form in the Universe, which can be used to tightly constrain models of galaxy formation and evolution. Although the halo was once considered a single component, evidence for its dichotomy has slowly emerged in recent years from inspection of small samples of halo objects. Here we show that the halo is indeed clearly divisible into two broadly overlapping structural components -- an inner and an outer halo -- that exhibit different spatial density profiles, stellar orbits and stellar metallicities (abundances of elements heavier than helium). The inner halo has a modest net prograde rotation, whereas the outer halo exhibits a net retrograde rotation and a peak metallicity one-third that of the inner halo. These properties indicate that the individual halo components probably formed in fundamentally different ways, through successive dissipational (inner) and dissipationless (outer) mergers and tidal disruption of proto-Galactic clumps.

Gaia is an ambitious space astrometry mission of ESA with a main objective to map the sky in astrometry and photometry down to a magnitude 20 by the end of the next decade. While the mission is built and operated by ESA and an industrial consortium, the data processing is entrusted to a consortium formed by the scientific community, which was formed in 2006 and formally selected by ESA one year later. The satellite will downlink around 100 TB of raw telemetry data over a mission duration of 5 years from which a very complex iterative processing will lead to the final science output: astrometry with a final accuracy of a few tens of microarcseconds, epoch photometry in wide and narrow bands, radial velocity and spectra for the stars brighter than 17 mag. We discuss the general principles and main difficulties of this very large data processing and present the organisation of the European Consortium responsible for its design and implementation.

I outline a method for estimating astrophysical parameters (APs) from multidimensional data. It is a supervised method based on matching observed data (e.g. a spectrum) to a grid of pre-labelled templates. However, unlike standard machine learning methods such as ANNs, SVMs or k-nn, this algorithm explicitly uses domain information to better weight each data dimension in the estimation. Specifically, it uses the sensitivity of each measured variable to each AP to perform a local, iterative interpolation of the grid. It avoids both the non-uniqueness problem of global regression as well as the grid resolution limitation of nearest neighbours.

Models of brown dwarf atmospheres suggest they exhibit complex physical behaviour. Observations have shown that they are indeed dynamic, displaying small photometric variations over timescales of hours. Here I report results of infrared (0.95-1.64 micron) spectrophotometric monitoring of four field L and T dwarfs spanning timescales of 0.1-5.5 hrs, the goal being to learn more about the physical nature of this variability. Spectra are analysed differentially with respect to a simultaneously observed reference source in order to remove Earth-atmospheric variations. The variability amplitude detected is typically 2-10%, depending on the source and wavelength. I analyse the data for correlated variations between spectral indices. This approach is more robust than single band or chisq analyses, because it does not assume an amplitude for the (often uncertain) noise level (although the significance test still assumes a shape for the noise power spectrum). Three of the four targets show significant evidence for correlated variability. Some of this can be associated with specific features including Fe, FeH, VO and KI, and there is good evidence for intrinsic variability in water and possibly also methan. Yet some of this variability covers a broader spectral range which would be consistent with dust opacity variations. The underlying common cause is plausibly localized temperature or composition fluctuations caused by convection. Looking at the high signal-to-noise ratio stacked spectra we see many previously identified spectral features of L and T dwarfs, such as KI, NaI, FeH, water and methane. In particular we may have detected methane absorption at 1.3-1.4 micron in the L5 dwarf SDSS 0539-0059.

We report high-resolution spectroscopy of 125 field stars previously observed as part of the Sloan Digital Sky Survey and its program for Galactic studies, the Sloan Extension for Galactic Understanding and Exploration (SEGUE). These spectra are used to measure radial velocities and to derive atmospheric parameters, which we compare with those reported by the SEGUE Stellar Parameter Pipeline (SSPP). The SSPP obtains estimates of these quantities based on SDSS ugriz photometry and low-resolution (R = 2000) spectroscopy. For F- and G-type stars observed with high signal-to-noise ratios (S/N), we empirically determine the typical random uncertainties in the radial velocities, effective temperatures, surface gravities, and metallicities delivered by the SSPP to be 2.4 km/s, 130 K (2.2%), 0.21 dex, and 0.11 dex, respectively, with systematic uncertainties of a similar magnitude in the effective temperatures and metallicities. We estimate random errors for lower S/N spectra based on numerical simulations.

We validate the performance and accuracy of the current SEGUE (Sloan Extension for Galactic Understanding and Exploration) Stellar Parameter Pipeline (SSPP), which determines stellar atmospheric parameters (effective temperature, surface gravity, and metallicity) by comparing derived overall metallicities and radial velocities from selected likely members of three globular clusters (M 13, M 15, and M 2) and two open clusters (NGC 2420 and M 67) to the literature values. Spectroscopic and photometric data obtained during the course of the original Sloan Digital Sky Survey (SDSS-I) and its first extension (SDSS-II/SEGUE) are used to determine stellar radial velocities and atmospheric parameter estimates for stars in these clusters. Based on the scatter in the metallicities derived for the members of each cluster, we quantify the typical uncertainty of the SSPP values, sigma([Fe/H]) = 0.13 dex for stars in the range of 4500 K < Teff < 7500 K and 2.0 < log g < 5.0, at least over the metallicity interval spanned by the clusters studied (-2.3 < [Fe/H] < 0). The surface gravities and effective temperatures derived by the SSPP are also compared with those estimated from the comparison of the color-magnitude diagrams with stellar evolution models; we find satisfactory agreement. At present, the SSPP underestimates [Fe/H] for near-solar-metallicity stars, represented by members of M 67 in this study, by about 0.3 dex.

We describe the development and implementation of the SEGUE (Sloan Extension for Galactic Exploration and Understanding) Stellar Parameter Pipeline (SSPP). The SSPP derives, using multiple techniques, radial velocities and the fundamental stellar atmospheric parameters (effective temperature, surface gravity, and metallicity) for AFGK-type stars, based on medium-resolution spectroscopy and $ugriz$ photometry obtained during the course of the original Sloan Digital Sky Survey (SDSS-I) and its Galactic extension (SDSS-II/SEGUE). The SSPP also provides spectral classification for a much wider range of stars, including stars with temperatures outside of the window where atmospheric parameters can be estimated with the current approaches. This is Paper I in a series of papers on the SSPP; it provides an overview of the SSPP, and initial tests of its performance using multiple data sets. Random and systematic errors are critically examined for the current version of the SSPP, which has been used for the sixth public data release of the SDSS (DR-6).

(Abridged) The nature of S Ori 70, a faint mid-T type object found towards the direction of the young sigma Orionis cluster, is still under debate. We intend to disentangle whether it is a field brown dwarf or a 3-Myr old planetary-mass member of the cluster. We report on near-infrared JHK_s and mid-infrared [3.6] and [4.5] IRAC/Spitzer photometry recently obtained for S Ori 70. The new near-infrared images (taken 3.82 yr after the discovery data) have allowed us to derive a very small proper motion (11.0 +/- 5.9 mas/yr) for this object, which is consistent with the proper motion of the cluster within 1.5 sigma the astrometric uncertainty. The colors (H-K_s), (J-K_s) and K_s-[3.6] appear discrepant when compared to T4-T7 dwarfs in the field. This behavior could be ascribed to a low-gravity atmosphere or alternatively to an atmosphere with a metallicity significantly different than solar. Taking into account the small proper motion of S Ori 70 and its new near- and mid-infrared colors, a low-gravity atmosphere remains as the most likely explanation to account for our observations. This supports S Ori 70's membership in sigma Orionis, with an estimated mass in the interval 2-7 Mjup, in agreement with our previous derivation.

Often the relation between the variables constituting a multivariate data space might be characterized by one or more of the terms: ``nonlinear'', ``branched'', ``disconnected'', ``bended'', ``curved'', ``heterogeneous'', or, more general, ``complex''. In these cases, simple principal component analysis (PCA) as a tool for dimension reduction can fail badly. Of the many alternative approaches proposed so far, local approximations of PCA are among the most promising. This paper will give a short review of localized versions of PCA, focusing on local principal curves and local partitioning algorithms. Furthermore we discuss projections other than the local principal components. When performing local dimension reduction for regression or classification problems it is important to focus not only on the manifold structure of the covariates, but also on the response variable(s). Local principal components only achieve the former, whereas localized regression approaches concentrate on the latter. Local projection directions derived from the partial least squares (PLS) algorithm offer an interesting trade-off between these two objectives. We apply these methods to several real data sets. In particular, we consider simulated astrophysical data from the future Galactic survey mission Gaia.

Aims:The Gaia astrometric survey mission will, as a consequence of its scanning law, obtain low resolution optical (330-1000 nm) spectrophotometry of several million unresolved galaxies brighter than V=22. We present the first steps in a project to design and implement a classification system for these data. The goal is both to determine morphological classes and to estimate intrinsic astrophysical parameters via synthetic templates. Here we describe (1) a new library of synthetic galaxy spectra, and (2) first results of classification and parametrization experiments using simulated Gaia spectrophotometry of this library. Methods:We have created a large grid of synthetic galaxy spectra using the PEGASE.2 code, which is based on galaxy evolution models that take into account metallicity evolution, extinction correction, emission lines (with stellar spectra based on the BaSeL library). Our classification and regression models are Support Vector Machines (SVMs), which are kernel-based nonlinear estimators. Results:We produce a basic library of about 4000 zero redshift galaxy spectra covering the main Hubble types over wavelength range 250 to 1050 nm at a sampling of 1 nm or less. It is computed on a regular grid of four key astrophysical parameters for each type and for intermediate random values of the same parameters. An extended library reproduces this at a series of redshifts. Initial 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. Comparing the colours of our synthetic library with Sloan Digital Sky Survey (SDSS) spectra we find good agreement over the full range of Hubble types and parameters.

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.

We present techniques for the estimation of stellar atmospheric parameters (Teff,logg,[Fe/H]) for stars from the SDSS/SEGUE survey. The atmospheric parameters are derived from the observed medium-resolution (R=2000) stellar spectra using non-linear regression models trained either on (1) pre-classified observed data or (2) synthetic stellar spectra. In the first case we use our models to automate and generalize parametrization produced by a preliminary version of the SDSS/SEGUE Spectroscopic Parameter Pipeline (SSPP). In the second case we directly model the mapping between synthetic spectra (derived from Kurucz model atmospheres) and the atmospheric parameters, independently of any intermediate estimates. After training, we apply our models to various samples of SDSS spectra to derive atmospheric parameters, and compare our results with those obtained previously by the SSPP for the same samples. We obtain consistency between the two approaches, with RMS deviations of 150K in Teff, 0.35dex in logg, and 0.22dex in [Fe/H]. The models are applied to pre-processed spectra, either via Principal Components Analysis or a Wavelength Range Selection method, which employs a subset of the full 3850-9000A spectral range. This is both for computational reasons, and because it delivers higher accuracy. From an analysis of cluster candidates with available SDSS spectroscopy (M15, M13, M2, and NGC2420), we find evidence for small systematic offsets in Teff and/or logg for the atmospheric parameter estimates from the model trained on real data with the SSPP. Thus, this model turns out to derive more precise, but less accurate, atmospheric parameters than the model trained on synthetic data.

Gaia is ESA's ambitious space astrometry mission the main objective of which is to astrometrically and spectro-photometrically map 1000 Million celestial objects (mostly in our galaxy) with unprecedented accuracy. The announcement of opportunity for the data processing will be issued by ESA late in 2006. The Gaia Data Processing and Analysis Consortium (DPAC) has been formed recently and is preparing an answer. The satellite will downlink close to 100 TB of raw telemetry data over 5 years. To achieve its required accuracy of a few 10s of Microarcsecond astrometry, a highly involved processing of this data is required. In addition to the main astrometric instrument Gaia will host a Radial Velocity instrument, two low-resolution dispersers for multi-color photometry and two Star Mappers. Gaia is a flying Giga Pixel camera. The various instruments each require relatively complex processing while at the same time being interdependent. We describe the overall composition of the DPAC and the envisaged overall architecture of the Gaia data processing system. We shall delve further into the core processing - one of the nine, so-called, coordination units comprising the Gaia processing system.

Gaia is a fully-approved all-sky astrometric and photometric survey due for launch in 2011. It will measure accurate parallaxes and proper motions for everything brighter than G=20 (ca. 10^9 stars). Its primary objective is to study the composition, origin and evolution of our Galaxy from the 3D structure, 3D velocities, abundances and ages of its stars. In some respects it can be considered as a cosmological survey at redshift zero. Several other upcoming space-based surveys, in particular JWST and Herschel, will study star and galaxy formation in the early (high-redshift) universe. In this paper I briefly describe these missions, as well as SIM and Jasmine, and explain why they need to observe from space. I then discuss some Galactic science contributions of Gaia concerning dark matter, the search for substructure, stellar populations and the mass--luminosity relation. The Gaia data are complex and require the development of novel analysis methods; here I summarize the principle of the astrometric processing. In the last two sections I outline how the Gaia data can be exploited in connection with other observational and theoretical work in order to build up a more comprehensive picture of galactic evolution.

We present a search for low-mass brown dwarfs in the Pleiades open cluster. The identification of Pleiades members fainter and cooler than those currently known allows us to constrain evolutionary models for L dwarfs and to extend the study of the cluster mass function to lower masses. We conducted a 1.8 deg^2 near-infrared J-band survey at the 3.5m Calar Alto Telescope, with completeness J~19.0. The detected sources were correlated with those of previously available optical I-band images (completeness I~22). Using a J versus I-J colour-magnitude diagram, we identified 18 faint red L-type candidates, with magnitudes 17.4<J<19.7 and colours I-J>3.2. If Pleiades members, their masses would span ~0.040-0.020 M_Sol. We performed follow-up HKs-band imaging to further confirm their cluster membership by photometry and proper motion. Out of 11 IJ candidates with proper motion measurements, we find six cluster members, two non-members and three whose membership is uncertain and depends on the intrinsic velocity dispersion of Pleiades brown dwarfs. This dispersion (>4 mas yr^-1) is at least four times that of cluster stars with masses >1 M_Sol. Five of the seven other IJ candidates are discarded because their J-Ks colours are bluer than those of confirmed members. The J versus I-J sequence of the L-type candidates at J>18 is not as red as theoretical models predict; it rather follows the field L-dwarf sequence translated to the cluster distance. This sequence overlapping, also observed in the J versus J-H and J-K diagrams, suggests that Pleiades and field L dwarfs may have similar spectral energy distributions and luminosities, and thus possibly similar radii. Also, we find alpha=0.5+-0.2 for a power-law approximation dN/dM propor. M^-alpha of the survey mass spectrum in the mass range 0.5-0.026 M_Sol.

We present simultaneous UV-G-R-I monitoring of 19 M dwarfs that revealed a huge flare on the M9 dwarf 2MASSW J1707183+643933 with an amplitude in the UV of at least 6 magnitudes. This is one of the strongest detections ever of an optical flare on an M star and one of the first in an ultracool dwarf (UCD, spectral types later than about M7). Four intermediate strength flares (Delta m_UV < 4 mag) were found in this and three other targets. For the whole sample we deduce a flare probability of 0.013 (rate of 0.018/hr), and 0.049 (0.090/hr) for 2M1707+64 alone. Deviations of the flare emission from a blackbody is consistent with strong Halpha line emission. We also confirm the previously found rotation period for 2M1707+64 (Rockenfeller, Bailer-Jones & Mundt (2006), http://arxiv.org/abs/astro-ph/0511614/) and determine it more precisely to be 3.619 +/- 0.015 hr.

The European Gaia astrometry mission is due for launch in 2011. Gaia will rely on the proven principles of ESA's Hipparcos mission to create an all-sky survey of about one billion stars throughout our Galaxy and beyond, by observing all objects down to 20th magnitude. Through its massive measurement of stellar distances, motions and multi-colour photometry it will provide fundamental data necessary for unravelling the structure, formation and evolution of the Galaxy. This paper presents the design and performance of the broad- and medium-band set of photometric filters adopted as the baseline for Gaia. The nineteen selected passbands (extending from the ultraviolet to the far-red), the criteria, and the methodology on which this choice has been based are discussed in detail. We analyse the photometric capabilities for characterizing the luminosity, temperature, gravity and chemical composition of stars. We also discuss the automatic determination of these physical parameters for the large number of observations involved, for objects located throughout the entire Hertzsprung-Russell diagram. Finally, the capability of the photometric system to deal with the main Gaia science case is outlined.

We present simultaneous, multiband monitoring of 19 field M dwarfs (spectral types M2-M9). Significant variability was found in seven objects in at least one out of the three channels I, R and G. Two objects, LHS370 (M5V) and 2M1707+64 (M9V), show periods of 5.9+-2.0 and 3.65+-0.1 hours respectively. On account of the agreement with the typical values of v*sin(i) published for M dwarfs (Mohanty & Basri 2003), we claim these to be the objects' rotation periods. In the case of 2M1707+64 we furthermore find evidence for the presence of magnetically-induced cool spots as the reason for the observed variability (at a temperature contrast of 4-8 percent, with a projected surface coverage factor of less than 0.075) and can rule out the scenario of dust clouds (as represented by the COND or DUSTY models). Comparing the frequency of occurrence of variability in this and various L dwarf samples published over the past few years, we find that variability is more common in field L dwarfs than in field M dwarfs. This implies a change in the physical nature and/or extent of surface features when moving from M to L dwarfs.

We present results of an optical and near-infrared (IR) 1.8 deg^2 survey in the Pleiades open cluster to search for substellar objects. From optical I-band images from the CFHT and J-band images from the 3.5 m CAHA Telescope, we identify 18 faint and very red L brown dwarf candidates, with I> 20.9 and I-J> 3.2. The follow-up observations of nine objects in the H- and Ks-bands confirm that eight belong to the IR sequence of the cluster and the proper motion measurements of seven candidates confirm that they are Pleiades members. A preliminary estimation of the substellar mass spectrum dN/dM in the form of a power law M^-alpha provides alpha=0.57+-0.14. We extrapolate this function to estimate the number of very low-mass brown dwarfs and planetary mass objects that could be present in the cluster down to 1 M_Jup. Sensitive searches combining far red and near infrared observations may unveal these objects in a near future.

We present results from the most recent (2002-2004) observing campaigns of the eclipsing system KH 15D, in addition to re-reduced data obtained at Van Vleck Observatory (VVO) between 1995 and 2000. Phasing nine years of photometric data shows substantial evolution in the width and depth of the eclipses. The most recent data indicate that the eclipses are now approximately 24 days in length, or half the orbital period. These results are interpreted and discussed in the context of the recent models for this system put forward by Winn et al. and Chiang & Murray-Clay. A periodogram of the entire data set yields a highly significant peak at 48.37 +/- 0.01 days, which is in accord with the spectroscopic period of 48.38 +/- 0.01 days determined by Johnson et al. Another significant peak, at 9.6 days, was found in the periodogram of the out-of-eclipse data at two different epochs. We interpret this as the rotation period of the visible star and argue that it may be tidally locked in pseudosynchronism with its orbital motion. If so, application of Hut's theory implies that the eccentricity of the orbit is e = 0.65 +/- 0.01. Analysis of the UVES/VLT spectra obtained by Hamilton et al. shows that the v sin(i) of the visible star in this system is 6.9 +/- 0.3 km/sec. Using this value of v sin(i) and the measured rotation period of the star, we calculate the lower limit on the radius to be R = (1.3 +/- 0.1), R_Sun, which concurs with the value obtained by Hamilton et al. from its luminosity and effective temperature. Here we assume that i = 90 degrees since it is likely that the spin and orbital angular momenta vectors are nearly aligned.

We have employed feedforward neural networks trained on synthetic spectra in the range 3800 to 5600 AA with resolutions of 2-3 AA to determine metallicities from spectra of about 1000 main-sequence turn-off, subgiant and red giant stars in the globular clusters M55 and omega Cen. The overall metallicity accuracies are of the order of 0.15 to 0.2 dex. In addition, we tested how well the stellar parameters logg and Teff can be retrieved from such data without additional colour or photometric information. We find overall uncertainties of 0.3 to 0.4 dex for logg and 140 to 190 K for Teff. In order to obtain some measure of uncertainty for the determined values of [Fe/H], logg and Teff, we applied the bootstrap method for the first time to neural networks for this kind of parametrization problem. The distribution of metallicities for stars in omega Cen clearly shows a large spread in agreement with the well known multiple stellar populations in this cluster.

Large, deep surveys must typically rely on multiband photometry rather than spectroscopy for determining the astrophysical properties (APs) of stars. Yet designing an optimal photometric system for a wide range of objects is complex, because it must trade off conflicting scientific requirements. I present a new method for designing photometric systems and apply it to the Gaia Galactic survey satellite, which will observe one billion stars brighter than V=20. The principle is to optimally sample stellar spectra in order to best determine the APs (e.g. Teff, [Fe/H], interstellar extinction). By considering a filter system as a set of free parameters (central wavelengths, FWHM etc.), it may be designed by optimizing a figure-of-merit (FoM) with respect to these parameters. The FoM is a measure of how well the filter system can vectorially 'separate' stars in the data space to avoid AP degeneracies. The resulting systems show some interesting features, in particular broad, overlapping filters, which may be desirable from a multivariate classification perspective. These systems are competitive with others proposed for Gaia.

Gaia will observe more than one billion objects brighter than V=20, including stars, asteroids, galaxies and quasars. As Gaia performs real time detection (i.e. without an input catalogue) the intrinsic properties of most of these objects will not be known a priori. An integral part of the Gaia data processing is therefore to classify everything observed. This will be based primarily on multiband photometry provided by Gaia, but should also make optimal use of the high resolution spectroscopy (for brighter stars) and the parallaxes. In addition to a broad classification, we can also determine fundamental stellar parameters, in particular effective temperature, metallicity and the line-of-sight interstellar extinction. Such information will be essential for fully exploiting the astrometric part of the Gaia catalogue for stellar population studies. However, extracting this information is a significant challenge, and will need to make use of appropriate multidimensional data analysis techniques. I outline some of the problems and the strategies being developed to tackle them.

Gaia is an all sky, high precision astrometric and photometric satellite of the European Space Agency (ESA) due for launch in 2010-2011. Its primary mission is to study the composition, formation and evolution of our Galaxy. Gaia will measure parallaxes and proper motions of every object in the sky brighter than V=20, amounting to a billion stars, galaxies, quasars and solar system objects. It will achieve an astrometric accuracy of 10muas at V=15 - corresponding to a distance accuracy of 1% at 1kpc. With Gaia, tens of millions of stars will have their distances measured to a few percent or better. This is an improvement over Hipparcos by several orders of magnitude in the number of objects, accuracy and limiting magnitude. Gaia will also measure radial velocities for source brighter than V~17. To characterize the objects, each object is observed in 15 medium and broad photometric bands with an onboard CCD camera. With these capabilities, Gaia will make significant advances in a wide range of astrophysical topics. These include a detailed kinematical map of stellar populations, stellar structure and evolution, the discovery and characterization of thousands of exoplanetary systems and General Relativity on large scales. I give an overview of the mission, its operating principles and its expected scientific contributions. For the latter I provide a quick look in five areas on increasing scale size in the universe: the solar system, exosolar planets, stellar clusters and associations, Galactic structure and extragalactic astronomy.

Over the past few years monitoring programs have shown ultra cool dwarfs (UCDs) to be photometrically variable. Of the 60 sources monitored in the field and some 120 monitored in clusters, about 40% show variability in both cases. For mid to late M dwarfs in young (<100 Myr) clusters, this variability is generally periodic with amplitudes of up to a few tenths of a magnitude and periods of between a few hours and several days. For older field dwarfs (covering late M, L and T types) this variability is often nonperiodic with smaller amplitudes (up to 0.1 mag in I) and timescales of order hours. The former may be attributed to the rotational modulation of magnetically-induced photospheric spots, as seen in higher mass T Tauri stars. The nonperiodic variability, on the other hand, may be caused by a rapid evolution of surface features (which 'mask' the otherwise observable rotational modulation). This could be related to the formation and dissipation of inhomogeneities in dust clouds in the photospheres of UCDs.

The variability and rotation of ultra cool dwarfs (UCDs) provide important information on their atmospheres and evolution. As part of an ongoing program to investigate this, projected rotation speeds, vsini, are presented for 16 field UCDs (M9V-L7.5V). This doubles the number of L dwarfs for which vsini has been measured. All targets are found to have vsini between 10 and 40 km/s, confirming that L dwarfs are rapid rotators. From the random distribution of the rotation axes, i, one-sided confidence intervals are placed on the rotation periods of individual object. These are compared with published period data obtained from photometric monitoring. From this, the period of 31 hrs for the L0.5 dwarf 2M0746+2000 published by Gelino et al. (2002) may be ruled out as the rotation period. The period of 11.2 +/- 0.8 hrs for the L1.5 dwarf 2M1145+2317 obtained by Bailer-Jones & Mundt (2001) is plausibly the rotation period. The inclination of the rotation axis is i = 62-90 deg with an expectation value of 76 deg. Alternatively the data set a lower limit on the radius of 0.1 Rsol. Similarly, the period of 2.7 +/- 0.1 hrs for 2M1334+1940 is the likely rotation period; the inclination is i = 27-44 deg (mean i = 34 deg). Where no variability or period was detected by the monitoring programs the likely reason is low contrast surface features. However, in three cases variability but no period was detected, despite sensitivity to the likely rotation periods inferred from vsini. This reinforces the `masking hypothesis', the idea that the evolution of photospheric features obscures regular modulation of the light curve. As has been previously discussed, a likely candidate for such features is inhomogeneous dust clouds.

Designing a photometric system to best fulfil a set of scientific goals is a complex task, demanding a compromise between conflicting requirements and subject to various constraints. A specific example is the determination of stellar astrophysical parameters (APs) - effective temperature, metallicity etc. - across a wide range of stellar types. I present a novel approach to this problem which makes minimal assumptions about the required filter system. By considering a filter system as a set of free parameters it may be designed by optimizing some figure-of-merit (FoM) with respect to these parameters. In the example considered, the FoM is a measure of how well the filter system can `separate' stars with different APs. This separation is vectorial in nature, in the sense that the local directions of AP variance are preferably mutually orthogonal to avoid AP degeneracy. The optimization is carried out with an evolutionary algorithm, which uses principles of evolutionary biology to search the parameter space. This model, HFD (Heuristic Filter Design), is applied to the design of photometric systems for the Gaia space astrometry mission. The optimized systems show a number of interesting features, not least the persistence of broad, overlapping filters. These HFD systems perform as least as well as other proposed systems for Gaia, although inadequacies remain in all. The principles underlying HFD are quite generic and may be applied to filter design for numerous other projects, such as the search for specific types of objects or photometric redshift determination.

We present the results of an extensive search for periodic and irregular variable pre-main sequence (PMS) stars in the young (2-4 Myr) open cluster NGC 2664, based on photometric monitoring using the Wide Field Imager (WFI) on the 2.2m telescope on La Silla (Chile). In total, about 10600 stars with I_c magnitudes between 9.8mag and 21mag have been monitored in our 34x33 arcmin field. Time series data were obtained in the I_c band in 44 nights between Dec. 2000 and March 2001; altogether we obtained 88 data points per star. Using two different time series analysis techniques (Scargle and CLEAN) we found 543 periodically variable stars with periods between 0.2 days and 15 days. Also, 484 irregular variable stars were identified using a chi-squared-test. In addition we have carried out nearly simultaneous observations in V, R_c and a narrow-band Halpha filter. The photometric data enable us to reject background and foreground stars from our sample of variable stars according to their location in the I_c vs. (R_c-I_c) colour-magnitude and (R_c-Halpha) vs. (R_c-I_c) colour-colour diagram. We could in the end identify 405 periodically variable and 184 irregular variable PMS stars as cluster members using these two different tests. In addition 35 PMS stars for which no significant variablilty were detected could be identified as members using an Halpha emission index criterion. This yields to a total of 624 PMS stars in NGC 2264, of which only 182 were known before. Most of the newly found PMS stars are fainter than I_c = 16 and of late spectral type (>M2). (abridged)

We present high resolution (R $\sim$ 44,000) UVES spectra of the eclipsing pre-main sequence star KH 15D covering the wavelength range 4780 to 6810 {\AA} obtained at three phases: out of eclipse, near minimum light and during egress. The system evidently acts like a natural coronagraph, enhancing the contrast relative to the continuum of hydrogen and forbidden emission lines during eclipse. At maximum light the H$\alpha$ equivalent width was $\sim$2 {\AA} and the profile showed broad wings and a deep central absorption. During egress the equivalent width was much higher ($\sim$70 {\AA}) and the broad wings, which extend to $\pm$ 300 km/s, were prominent. During eclipse totality the equivalent width was less than during egress ($\sim$40 {\AA}) and the high velocity wings were much weaker. H$\beta$ showed a somewhat different behavior, revealing only the blue-shifted portion of the high velocity component during eclipse and egress. [OI] $\lambda\lambda$6300, 6363 lines are easily seen both out of eclipse and when the photosphere is obscured and exhibit little or no flux variation with eclipse phase. Our interpretation is that KH 15D, although clearly a weak-line T Tauri star by the usual criteria, is still accreting matter from a circumstellar disk, and has a well-collimated bipolar jet. As the knife-edge of the occulting matter passes across the close stellar environment it is evidently revealing structure in the magnetosphere of this pre-main sequence star with unprecedented spatial resolution. We also show that there is only a small, perhaps marginally significant, change in the velocity of the K7 star between the maximum light and egress phases probed here.

VLT/FORS spectroscopy and 2MASS near-infrared photometry, together with previously known data, have been used to establish the membership and the properties of a sample of low-mass candidate members of the sigma Orionis cluster with masses spanning from 1 Msun down to about 0.013 Msun (i.e., deuterium-burning mass limit). We have observed K-band infrared excess and remarkably intense H(alpha) emission in various cluster members, which, in addition to the previously detected forbidden emision lines and the presence of LiI in absorption at 6708 A, have allowed us to tentatively classify sigma Orionis members as classical or weak-line TTauri stars and substellar analogs. Variability of the H(alpha) line has been investigated and detected in some objects. Based on the K-band infrared excesses and the intensity of H(alpha) emission, we estimate that the minimum disk frequency of the sigma Orionis low-mass population is in the range 5-12%.

We have assessed how well stellar parameters (T_eff, logg and [Fe/H]) can be retrieved from low-resolution dispersed images to be obtained by the DIVA satellite. Although DIVA is primarily an all-sky astrometric mission, it will also obtain spectrophotometric information for about 13 million stars (operational limiting magnitude V ~ 13.5 mag). Constructional studies foresee a grating system yielding a dispersion of ~200nm/mm on the focal plane (first spectral order). For astrometric reasons there will be no cross dispersion which results in the overlapping of the first to third diffraction orders. The one-dimensional, position related intensity function is called a DISPI (DISPersed Intensity). We simulated DISPIS from synthetic spectra (...) for a limited range of metallicites i.e. our results are for [Fe/H] in the range -0.3 to 1 dex. We show that there is no need to deconvolve these low resolution signals in order to obtain basic stellar parameters. Using neural network methods and by including simulated data of DIVA's UV telescope, we can determine T_eff to an average accuracy of about 2% for DISPIS from stars with 2000 K < T_eff < 20000 K and visual magnitudes of V=13 mag (end of mission data). logg can be determined for all temperatures with an accuracy better than 0.25 dex for magnitudes brighter than V=12 mag. For low temperature stars with 2000 K < T_eff < 5000 K and for metallicities in the range -0.3 to +1 dex a determination of [Fe/H] is possible (to better than 0.2 dex) for these magnitudes. Additionally we examined the effects of extinction E(B-V) on DISPIS and found that it can be determined to better than 0.07 mag for magnitudes brighter than V=14 mag if the UV information is included.

Recent monitoring programs of ultra cool field M and L dwarfs (low mass stars or brown dwarfs) have uncovered low amplitude photometric I-band variations which may be associated with an inhomogeneous distribution of photospheric condensates. Further evidence hints that this distribution may evolve on very short timescales, specifically of order a rotation period or less. In an attempt to study this behaviour in more detail, we have carried out a pilot program to monitor three L dwarfs in the near infrared where these objects are significantly brighter than at shorter wavelengths. We present a robust data analysis method for improving the precision and reliability of infrared photometry. No significant variability was detected in either the J or Km bands in 2M1439 and SDSS1203 above a peak-to-peak amplitude of 0.04 mag (0.08 mag for 2M1112). The main limiting factor in achieving lower detection limits is suspected to be second order extinction effects in the Earth's atmosphere, on account of the very different colours of the target and reference stars. Suggestions are given for overcoming such effects which should improve the sensitivity and reliability of infrared variability searches.

Results of an international campaign to photometrically monitor the unique pre-main sequence eclipsing object KH 15D are reported. An updated ephemeris for the eclipse is derived that incorporates a slightly revised period of 48.36 d. There is some evidence that the orbital period is actually twice that value, with two eclipses occurring per cycle. The extraordinary depth (~3.5 mag) and duration (~18 days) of the eclipse indicate that it is caused by circumstellar matter, presumably the inner portion of a disk. The eclipse has continued to lengthen with time and the central brightness reversals are not as extreme as they once were. V-R and V-I colors indicate that the system is slightly bluer near minimum light. Ingress and egress are remarkably well modeled by the passage of a knife-edge across a limb-darkened star. Possible models for the system are briefly discussed.

I present the results of a program to spectrophotometrically monitor the L1.5 dwarf 2MASSW J1145572+231730 to identify the cause of photospheric variability in ultra cool dwarfs. Plausible candidates are magnetically-induced star spots and inhomogeneous photospheric dust clouds. Based on the atmospheric models and synthetic spectra of Allard et al. (2001), the expected signatures of these phenomena in the 0.5-2.5 um wavelength region are presented and discussed. Near infrared spectra of 2M1145 were obtained along with a nearby reference star observed simultaneously in the spectrograph slit. No convincing variability over a 54 hour interval is found in any one of several colour indices designed to be most sensitive to dust- and spot-related variability. Nonetheless, a significant correlation between the variability of two colour indices is found. This is slightly more consistent with the dust-related variability model than the cool spot one considered. Based on the theoretically predicted signatures and the median errors in the colour indices (0.03-0.05 magnitudes), upper limits are placed on the coverage of possible spots and clouds. Assuming the L dwarf to be best modelled by a dusty atmosphere at 1900K, coherent clear clouds are limited to a coverage of 10-15% of the projected surface area and 200K cooler spots to a 20% coverage. A larger coverage of many small features varying incoherently cannot be ruled out with this method. A lower effective temperature restricts coherent clear clouds to be much smaller; a higher temperature allows both clouds and spots to be larger. These upper limits are consistent with the two separate variability detections in the I-band reported by Bailer-Jones & Mundt (2001).

The GAIA Galactic survey satellite will obtain photometry in 15 filters of over 10^9 stars in our Galaxy across a very wide range of stellar types. No other planned survey will provide so much photometric information on so many stars. I examine the problem of how to determine fundamental physical parameters (Teff, log g, [Fe/H] etc.) from these data. Given the size, multidimensionality and diversity of this dataset, this is a challenging task beyond any encountered so far in large-scale stellar parametrization. I describe the problems faced (initial object identification, interstellar extinction, multiplicity, missing data etc.) and present a framework in which they can can be addressed. A probabilistic approach is advocated on the grounds that it can take advantage of additional information (e.g. priors and data uncertainties) in a consistent and useful manner, as well as give meaningful results in the presence of poor or degenerate data. Furthermore, I suggest an approach to parametrization which can use the other information GAIA will acquire, in particular the parallax, which has not previously been available for large-scale multidimensional parametrization. Several of the problems identified and ideas suggested will be relevant to other large surveys, such as SDSS, DIVA, FAME, VISTA and LSST, as well as stellar parametrization in a virtual observatory.

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.

We have determined new rotation periods for 404 stars in the Orion Nebula Cluster using the Wide Field Imager attached to the MPG/ESO 2.2 m telescope on La Silla, Chile. Mass estimates are available for 335 of these and most have M < 0.3 M_sun. We confirm the existence of a bimodal period distribution for the higher mass stars in our sample and show that the median rotation rate decreases with increasing mass for stars in the range 0.1 < M <0.4 M_sun. While the spread in angular momentum (J) at any given mass is more than a factor of 10, the majority of lower mass stars in the ONC rotate at rates approaching 30% of their critical break-up velocity, as opposed to 5-10% for solar-like stars. This is a consequence of both a small increase in observed specific angular momentum (j=J/M) and a larger decrease in the critical value of j with decreasing mass. Perhaps the most striking fact, however, is that j varies by so little - less than a factor of two - over the interval 0.1-1.0 M_sun. The distribution of rotation rates with mass in the ONC (age ~ 1 My) is similar in nature to what is found in the Pleiades (age ~ 100 My). These observations provide a significant new guide and test for models of stellar angular momentum evolution during the proto-stellar and pre-main sequence phases.

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.

Based on near-infrared polarimetric observations we constrain the degree of linear polarization of the afterglow light of GRB 000301C to less than 30% 1.8 days after the burst.

We present Ulysses and NEAR data from the detection of the short or intermediate duration (2 s) gamma-ray burst GRB000301C (2000 March 1.41 UT). The gamma-ray burst (GRB) was localised by the Inter Planetary Network (IPN) and RXTE to an area of 50 arcmin^2. A fading optical counterpart was subsequently discovered with the Nordic Optical Telescope (NOT) about 42h after the burst. The GRB lies at the border between the long-soft and the short-hard classes of GRBs. If GRB000301C belongs to the latter class, this would be the first detection of an afterglow to a short-hard burst. We present UBRI and JHK photometry from the time of the discovery until 11 days after the burst. Finally, we present spectroscopic observations of the optical afterglow obtained with the ESO VLT Antu telescope 4 and 5 days after the burst. The optical light curve is consistent with being achromatic from 2 to 11 days after the burst and exhibits a break. A broken power-law fit yields a shallow pre-break decay power-law slope of a_1=-0.72+-0.06, a break time of t_b=4.39+-0.26 days after the burst, and a post-break slope of a_2=-2.29+-0.17, which is best explained by a sideways expanding jet in an ambient medium of constant mean density. In the optical spectrum we find absorption features that are consistent with FeII, CIV, CII, SiII and Ly-a at a redshift of 2.0404+-0.0008. We find evidence for a curved shape of the spectral energy distribution of the observed afterglow. It is best fitted with a power-law spectral distribution with index b ~ -0.7 reddened by an SMC-like extinction law with A_V~0.1 mag. Based on the Ly-a absorption line we estimate the HI column density to be log(N(HI))=21.2+-0.5. This is the first direct indication of a connection between GRB host galaxies and Damped Ly-a Absorbers.

Artificial neural networks are algorithms which have been developed to tackle a range of computational problems. These range from modelling brain function to making predictions of time-dependent phenomena to solving hard (NP-complete) problems. In this introduction we describe a single, yet very important, type of network known as a feedforward network. This network is a mathematical model which can be trained to learn an arbitrarily complex relationship between a data and a parameter domain, so can be used to solve interpolation and classification problems. We discuss the structure, training and interpretation of these networks, and their implementation, taking the classification of stellar spectra as an example.

Current and future large astronomical surveys will yield multiparameter databases on millions or even billions of objects. The scientific exploitation of these will require powerful, robust, and automated classification tools tailored to the specific survey. Partly motivated by this, the past five to ten years has seen a significant increase in the amount of work focused on automated classification and its application to astronomical data. In this article, I review this work and assess the current status of automated stellar classification, with particular regard to its potential application to large astronomical surveys. I examine both the strengths and weaknesses of the various techniques and how they have been applied to different classification and parametrization problems. I finish with a brief look at the developments still required in order to apply a stellar classifier to a large survey.

Photometric I band light curves of 21 ultra cool M and L dwarfs are presented. Variability with amplitudes of 0.01 to 0.055 magnitudes (RMS) with typical timescales of an hour to several hours are discovered in half of these objects. Periodic variability is discovered in a few cases, but interestingly several variable objects show no significant periods, even though the observations were almost certainly sensitive to the expected rotation periods. It is argued that in these cases the variability is due to the evolution of the surface features on timescales of a few hours. This is supported in the case of 2M1145 for which no common period is found in two separate light curves. It is speculated that these features are photospheric dust clouds, with their evolution possibly driven by rotation and turbulence. An alternative possibility is magnetically-induced surface features. However, chromospheric activity undergoes a sharp decrease between M7 and L1, whereas a greater occurrence of variability is observed in objects later than M9, lending support to the dust interpretation.

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)

Large scale, deep survey missions such as GAIA will collect enormous amounts of data on a significant fraction of the stellar content of our Galaxy. These missions will require a careful optimisation of their observational systems in order to maximise their scientific return, and will require reliable and automated techniques for parametrizing the very large number of stars detected. To address these two problems, I investigate the precision to which the three principal stellar parameters (Teff, logg, [M/H]) can be determined as a function of spectral resolution and signal-to-noise (SNR) ratio, using a large grid of synthetic spectra. The parametrization technique is a neural network, which is shown to provide an accurate three-dimensional physical parametrization of stellar spectra across a wide range of parameters. It is found that even at low resolution (50-100 AA FWHM) and SNR (5-10 per resolution element), Teff and \met can be determined to 1% and 0.2 dex respectively across a large range of temperatures (4000-30000 K) and metallicities (-3.0 to +1.0 dex), and that logg is measurable to +/- 0.2 dex for stars earlier than solar. The accuracy of the results is probably limited by the finite parameter sampling of the data grid. The ability of medium band filter systems (with 10-15 filters) for determining stellar parameters is also investigated. Although easier to implement in a unpointed survey, it is found that they are only competitive at higher SNRs (> 50).

The ongoing development of large infrared array detectors has enabled wide field, deep surveys to be undertaken. There are, however, a number of challenges in building an infrared instrument which has both excellent optical quality and high sensitivity over a wide field. We discuss these problems in the context of building a wide field imaging camera for the 3.5m telescope at Calar Alto with the new 2K*2K HgCdTe HAWAII-2 focal plane array. Our final design is a prime focus camera with a 15' field-of-view, called Omega 2000. To achieve excellent optical quality over the whole field, we have had to dispense with the reimaging optics and cold Lyot stop. We show that creative baffling schemes, including the use of undersized baffles, can compensate for the lost K band sensitivity. A moving baffle will be employed in Omega 2000 to allow full transmission in the non-thermal J and H bands.

We have undertaken a differential photometric I band search for variability in three Pleiades brown dwarfs and three very low mass field L dwarfs. Analysis of the resultant time series for the Pleiades brown dwarfs (Teide 1, Calar 3, Roque 11) gives no evidence for variability (above 99% confidence) greater than 0.05 mags in any of these objects on timescales between 25 minutes and 27 hours. Time series of two of the L dwarfs, 2MASSW J0913032+184150 and 2MASSW J1146345+223053, are also consistent with no intrinsic variability (no significant variability above 0.07 mags for 30 mins < t < 126 hours, and above 0.025 mags for 30 mins < t < 75 hours, for the two objects respectively). However, the L dwarf 2MASSW J1145572+231730 shows evidence for variability of amplitude 0.04 magnitudes with a tentative period of about 7 hours. It is therefore interesting that Kirkpatrick et al. (1999) report an Halpha emission line of 4.2AA equivalent width for this star: If our detection is confirmed as a rotation period, it would support the link between rotation and Halpha emission in very low mass stars.

We describe a method for the extraction of spectra from high dispersion objective prism plates. Our method is a catalogue driven plate solution approach, making use of the Right Ascension and Declination coordinates for the target objects. In contrast to existing methods of photographic plate reduction, we digitize the entire plate and extract spectra off-line. This approach has the advantages that it can be applied to CCD objective prism images, and spectra can be re-extracted (or additional spectra extracted) without having to re-scan the plate. After a brief initial interactive period, the subsequent reduction procedure is completely automatic, resulting in fully-reduced, wavelength justified spectra. We also discuss a method of removing stellar continua using a combination of non-linear filtering algorithms. The method described is used to extract over 12,000 spectra from a set of 92 objective prism plates. These spectra are used in an associated project to develop automated spectral classifiers based on neural networks.

We investigate the application of neural networks to the automation of MK spectral classification. The data set for this project consists of a set of over 5000 optical (3800-5200 AA) spectra obtained from objective prism plates from the Michigan Spectral Survey. These spectra, along with their two-dimensional MK classifications listed in the Michigan Henry Draper Catalogue, were used to develop supervised neural network classifiers. We show that neural networks can give accurate spectral type classifications (sig_68 = 0.82 subtypes, sig_rms = 1.09 subtypes) across the full range of spectral types present in the data set (B2-M7). We show also that the networks yield correct luminosity classes for over 95% of both dwarfs and giants with a high degree of confidence. Stellar spectra generally contain a large amount of redundant information. We investigate the application of Principal Components Analysis (PCA) to the optimal compression of spectra. We show that PCA can compress the spectra by a factor of over 30 while retaining essentially all of the useful information in the data set. Furthermore, it is shown that this compression optimally removes noise and can be used to identify unusual spectra.

We present a technique which employs artificial neural networks to produce physical parameters for stellar spectra. A neural network is trained on a set of synthetic optical stellar spectra to give physical parameters (e.g. T_eff, log g, [M/H]). The network is then used to produce physical parameters for real, observed spectra. Our neural networks are trained on a set of 155 synthetic spectra, generated using the SPECTRUM program written by Gray (Gray & Corbally 1994, Gray & Arlt 1996). Once trained, the neural network is used to yield T_eff for over 5000 B-K spectra extracted from a set of photographic objective prism plates (Bailer-Jones, Irwin & von Hippel 1997a). Using the MK classifications for these spectra assigned by Houk (1975, 1978, 1982, 1988) we have produced a temperature calibration of the MK system based on this set of 5000 spectra. It is demonstrated through the metallicity dependence of the derived temperature calibration that the neural networks are sensitive to the metallicity signature in the real spectra. With further work it is likely that neural networks will be able to yield reliable metallicity measurements for stellar spectra.