galev evolutionary synthesis models – I. Code, input physics and web interface

galev (GALaxy EVolution) evolutionary synthesis models describe the evolution of stellar populations in general, of star clusters as well as of galaxies, both in terms of resolved stellar populations and of integrated light properties over cosmological time-scales of ≥13 Gyr from the onset of star formation shortly after the big bang until today.
For galaxies, galev includes a simultaneous treatment of the chemical evolution of the gas and the spectral evolution of the stellar content, allowing for what we call a chemically consistent treatment: we use input physics (stellar evolutionary tracks, stellar yields and model atmospheres) for a large range of metallicities and consistently account for the increasing initial abundances of successive stellar generations.
Here we present the latest version of the galev evolutionary synthesis models that are now interactively available at http://www.galev.org . We review the currently used input physics, and also give details on how this physics is implemented in practice. We explain how to use the interactive web interface to generate models for user-defined parameters and also give a range of applications that can be studied using galev , ranging from star clusters, undisturbed galaxies of various types E–Sd to starburst and dwarf galaxies, both in the local and the high-redshift Universe.     

The spatial evolution of stellar structures in the Large Magellanic Cloud

We present an analysis of the spatial distribution of various stellar populations within the Large Magellanic Cloud (LMC). We combine mid-infrared selected young stellar objects, optically selected samples with mean ages between ∼9 and ∼1000 Myr and existing stellar cluster catalogues to investigate how stellar structures form and evolve within the LMC. For the analysis we use Fractured Minimum Spanning Trees, the statistical Q parameter and the two-point correlation function. Restricting our analysis to young massive (OB) stars, we confirm our results obtained for M33, namely that the luminosity function of the groups is well described by a power law with index −2, and that there is no characteristic length-scale of star-forming regions. We find that stars in the LMC are born with a large amount of substructure, consistent with a two-dimensional fractal distribution with dimension     and evolve towards a uniform distribution on a time-scale of ∼175 Myr. This is comparable to the crossing time of the galaxy, and we suggest that stellar structure, regardless of spatial scale, will be eliminated in a crossing time. This may explain the smooth distribution of stars in massive/dense young clusters in the Galaxy, while other, less massive, clusters still display large amounts of structure at similar ages. By comparing the stellar and star cluster distributions and evolving time-scales, we show that infant mortality of clusters (or 'popping clusters') has a negligible influence on the galactic structure. Finally, we quantify the influence of the elongation, differential extinction and contamination of a population on the measured Q value.     

Fast identification of transits from light-curves

We present an algorithm that allows fast and efficient detection of transits, including planetary transits, from light-curves. The method is based on building an ensemble of fiducial models and compressing the data using the moped compression algorithm. We describe the method and demonstrate its efficiency by finding planet-like transits in simulated Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) light-curves. We show that our method is independent of the size of the search space of transit parameters. In large sets of light-curves, we achieve speed-up factors of the order of 10³ times over an optimized adaptive search in the χ² space. We discuss how the algorithm can be used in forthcoming large surveys like Pan-STARRS and the Large Synoptic Survey Telescope (LSST), and how it may be optimized for future space missions like Kepler and COROT where most of the processing must be done on board.     

Planetary transit observations at the University Observatory Jena: TrES‐2

We report on observations of several transit events of the transiting planet TrES‐2 obtained with the Cassegrain‐Teleskop‐Kamera at the University Observatory Jena. Between March 2007 and November 2008 ten different transits and almost a complete orbital period were observed. Overall, in 40 nights of observation 4291 exposures (in total 71.52 h of observation) of the TrES‐2 parent star were taken. With the transit timings for TrES‐2 from the 34 events published by the TrES‐network, the Transit Light Curve project and the Exoplanet Transit Database plus our own ten transits, we find that the orbital period is P = (2.470614 ± 0.000001) d, a slight change by ∼0.6 s compared to the previously published period. We present new ephemeris for this transiting planet. Furthermore, we found a second dip after the transit which could either be due to a blended variable star or occultation of a second star or even an additional object in the system. Our observations will be useful for future investigations of timing variations caused by additional perturbing planets and/or stellar spots and/or moons (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Results from the Wide Angle Search for Planets Prototype (WASP0) – II. Stellar variability in the Pegasus field

Recent wide field photometric surveys, which target a specific field for long durations, are ideal for studying both long- and short-period stellar variability. Here, we report on 75 variable stars detected during the observations of a field in Pegasus using the Wide Angle Search for Planets Prototype (WASP0) instrument, 73 of which are new discoveries. The variables detected include 16 δ Scuti stars, 34 eclipsing binaries, 3 BY Draconis stars and 4 RR Lyraes. We estimate that the fraction of stars in the field brighter than   V ∼ 13.5  exhibiting variable behaviour with an amplitude greater than 0.6 per cent rms is ∼0.4 per cent. These results are compared with other wide field stellar variability surveys, and implications for detecting transits due to extra-solar planets are discussed.     

A revised orbital ephemeris for HAT-P-9b

We present here three transit observations of HAT-P-9b taken on 14 February 2010 and 05 April 2010 UT from the University of Arizona’s 1.55 m Kuiper telescope on Mt. Bigelow. Our two light curves were obtained in the Arizona-I filter for all our observations, and underwent the same reduction process. Both transits occurred approximately 24 min earlier than expected from the ephemeris of Shporer et al. (2009). However, due to the large time span between our observed transits and those of Shporer et al. (2009), a 6.5 s shift downwards in orbital period from the value of Shporer et al. (2009) is sufficient to explain all available transit data. We find a new period of 3.922814 ± 0.000002 days for HAT-P-9b, an order of magnitude more precise than previous measurements, with no evidence for significant nonlinearities in the transit period.     

Upper limits on the hot Jupiter fraction in the field of NGC 7789

We describe a method of estimating the abundance of short-period extra-solar planets based on the results of a photometric survey for planetary transits. We apply the method to a 21-night survey with the 2.5-m Isaac Newton Telescope of ∼32 000 stars in a ∼0.5 × 0.5 deg² field including the open cluster NGC 7789. From the colour–magnitude diagram, we estimate the mass and radius of each star by comparison with the cluster main sequence. We search for injected synthetic transits throughout the light curve of each star in order to determine their recovery rate, and thus calculate the expected number of transit detections and false alarms in the survey. We take proper account of the photometric accuracy, time sampling of the observations and criteria (signal-to-noise ratio and number of transits) adopted for transit detection. Assuming that none of the transit candidates found in the survey will be confirmed as real planets, we place conservative upper limits on the abundance of planets as a function of planet radius, orbital period and spectral type.     

Planetary transit observations at the University Observatory Jena: XO‐1b and TrES‐1

We report on observations of transit events of the transiting planets XO‐1b and TrES‐1 with a 25 cm telescope of the University Observatory Jena. With the transit timings for XO‐1b from all 50 available XO, SuperWASP, Transit Light Curve (TLC)‐Project‐ and Exoplanet Transit Database (ETD)‐data, including our own I ‐band photometry obtained in March 2007, we find that the orbital period is P = (3.941501 ± 0.000001) d, a slight change by ∼3 s compared to the previously published period. We present new ephemeris for this transiting planet. Furthermore, we present new R ‐band photometry of two transits of TrES‐1. With the help of all available transit times from literature this allows us to refine the estimate of the orbital period: P = (3.0300722 ± 0.0000002) d. Our observations will be useful for future investigations of timing variations caused by additional perturbing planets and/or stellar spots and/or moons (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Black holes and core expansion in massive star clusters

In this study we present the results from realistic N -body modelling of massive star clusters in the Magellanic Clouds. We have computed eight simulations with   N ∼ 10⁵  particles; six of these were evolved for at least a Hubble time. The aim of this modelling is to examine in detail the possibility of large-scale core expansion in massive star clusters, and search for a viable dynamical origin for the radius–age trend observed for such objects in the Magellanic Clouds. We identify two physical processes which can lead to significant and prolonged cluster core expansion – mass-loss due to rapid stellar evolution in a primordially mass-segregated cluster, and heating due to a retained population of stellar mass black holes, formed in the supernova explosions of the most massive cluster stars. These two processes operate over different time-scales and during different periods of a cluster's life. The former occurs only at early times and cannot drive core expansion for longer than a few hundred Myr, while the latter typically does not begin until several hundred Myr have passed, but can result in core expansion lasting for many Gyr. We investigate the behaviour of each of these expansion mechanisms under different circumstances – in clusters with varying degrees of primordial mass segregation, and in clusters with varying black hole retention fractions. In combination, the two processes can lead to a wide variety of evolutionary paths on the radius–age plane, which fully cover the observed cluster distribution and hence define a dynamical origin for the radius–age trend in the Magellanic Clouds. We discuss in some detail the implications of core expansion for various aspects of globular cluster research, as well as the possibility of observationally inferring the presence of a significant population of stellar mass black holes in a cluster.     

Periastron precession measurements in transiting extrasolar planetary systems at the level of general relativity

Transiting exoplanetary systems are surpassingly important among the planetary systems since they provide the widest spectrum of information for both the planet and the host star. If a transiting planet is on an eccentric orbit, the duration of transits T _D is sensitive to the orientation of the orbital ellipse relative to the line of sight. The precession of the orbit results in a systematic variation in both the duration of individual transit events and the observed period between successive transits,   P _obs  . The periastron of the ellipse slowly precesses due to general relativity and possibly the presence of other planets in the system. This secular precession can be detected through the long-term change in   P _obs  (transit timing variations, TTV) or in T _D (transit duration variations, TDV). We estimate the corresponding precession measurement precision for repeated future observations of the known eccentric transiting exoplanetary systems (XO-3b, HD 147506b, GJ 436b and HD 17156b) using existing or planned space-borne instruments. The TDV measurement improves the precession detection sensitivity by orders of magnitude over the TTV measurement. We find that TDV measurements over a approximately 4 yr period can typically detect the precession rate to a precision well exceeding the level predicted by general relativity.     

Open clusters and the galactic disk

It is textbook knowledge that open clusters are conspicuous members of the thin disk of our Galaxy, but their role as contributors to the stellar population of the disk was regarded as minor. Starting from a homogenous stellar sky survey, the ASCC‐2.5, we revisited the population of open clusters in the solar neighbourhood from scratch. In the course of this enterprise we detected 130 formerly unknown open clusters, constructed volume‐ and magnitude‐limited samples of clusters, re‐determined distances, motions, sizes, ages, luminosities and masses of 650 open clusters. We derived the present‐day luminosity and mass functions of open clusters (not the stellar mass function in open clusters), the cluster initial mass function CIMF and the formation rate of open clusters. We find that open clusters contributed around 40 percent to the stellar content of the disk during the history of our Galaxy. Hence, open clusters are important building blocks of the Galactic disk (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Front‐ vs. back‐illuminated CCD cameras for photometric surveys: a noise budget analysis

Exoplanetary transit and stellar oscillation surveys require a very high precision photometry. The instrumental noise has therefore to be minimized. First, we perform a semi‐analytical model of different noise sources. We show that the noise due the CCD electrodes can be overcome using a Gaussian PSF (Point Spread Function) of full width half maximum larger than 1.6 pixels. We also find that for a PSF size of a few pixels, the photometric aperture has to be at least 2.5 times larger than the PSF full width half maximum. Then, we compare a front‐ with a back‐illuminated CCD through a Monte‐Carlo simulation. Both cameras give the same results for a PSF full width half maximum larger than 1.5 pixels. All these simulations are applied to the A STEP (Antarctica Search for Transiting Extrasolar Planets) project. As a result, we choose a front‐illuminated camera for A STEP because of its better resolution and lower price, and we will use a PSF larger than 1.6 pixels. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Detectability of exoplanetary transits from radial velocity surveys

Of the known transiting extrasolar planets, a few have been detected through photometric follow-up observations of radial velocity planets. Perhaps the best known of these is the transiting exoplanet HD 209458b. For hot Jupiters (periods less than ∼5 d), the a priori information that 10 per cent of these planets will transit their parent star due to the geometric transit probability leads to an estimate of the expected transit yields from radial velocity surveys. The radial velocity information can be used to construct an effective photometric follow-up strategy which will provide optimal detection of possible transits. Since the planet-harbouring stars are already known in this case, one is only limited by the photometric precision achievable by the chosen telescope/instrument. The radial velocity modelling code presented here automatically produces a transit ephemeris for each planet data set fitted by the program. Since the transit duration is brief compared with the fitted period, we calculate the maximum window for obtaining photometric transit observations after the radial velocity data have been obtained, generalizing for eccentric orbits. We discuss a typically employed survey strategy which may contribute to a possible radial velocity bias against detection of the very hot Jupiters which have dominated the transit discoveries. Finally, we describe how these methods can be applied to current and future radial velocity surveys.     

The Kepler Mission: Astrophysics and Eclipsing Binaries

The Kepler Mission is a photometric space mission that will continuously observe a single 100 square degree field of view (FOV) of the sky of more than 100,000 stars in the Cygnus-Lyra region for four or more years with a precision of 14 parts per million (ppm) for a 6.5 hour integration including shot noise for a twelfth magnitude star. The primary goal of the mission is to detect Earth-size planets in the habitable zone of solar-like stars. In the process, many eclipsing binaries (EB) will also be detected. Prior to launch, the stellar characteristics will have been determined for all the stars in the FOV with K<14.5. As part of the verification process, stars with transits (about 5%) will need to have follow-up radial velocity observations performed to determine the component masses and thereby separate grazing eclipses caused by stellar companions from transits caused by planets. The result will be a rich database on EBs. The community will have access to the archive for uses such as for EB modeling of the high-precision light curves. A guest observer program is also planned for objects not already on the target list.     

Collaborative research of open star clusters

Preliminary results on observations of open clusters are presented. The project has been initiated in the framework of the Uzbek-Taiwan and Taiwan-Baltic collaboration, mainly to upgrade and make use of facilities at Maidanak Observatory. We present detailed, multiwavelength studies of the young cluster NGC 6823 and the associated complex nebulosity, to diagnose the young stellar population and star formation history in the region. In addition, 7 compact open clusters have been monitored for stellar variability. We show how observations like these could feasibly be used to look for exoplanet transit events. We also expect to join the Whole-Earth Telescope effort in future campaigns for asteroseismology.     

The period–luminosity relation for type II Cepheids in globular clusters

We report the result of our near-infrared observations ( JHK _s) for type II Cepheids (including possible RV Tau stars) in galactic globular clusters. We detected variations of 46 variables in 26 clusters (10 new discoveries in seven clusters) and present their light curves. Their periods range from 1.2 d to over 80 d. They show a well-defined period–luminosity relation at each wavelength. Two type II Cepheids in NGC 6441 also obey the relation if we assume the horizontal branch stars in NGC 6441 are as bright as those in metal-poor globular clusters in spite of the high metallicity of the cluster. This result supports the high luminosity which has been suggested for the RR Lyr variables in this cluster. The period–luminosity relation can be reproduced using the pulsation equation     assuming that all the stars have the same mass. Cluster RR Lyr variables were found to lie on an extrapolation of the period–luminosity relation. These results provide important constraints on the parameters of the variable stars.
Using Two Micron All-Sky Survey (2MASS) data, we show that the type II Cepheids in the Large Magellanic Cloud (LMC) fit our period–luminosity relation within the expected scatter at the shorter periods. However, at long periods (   P > 40  d, i.e. in the RV Tau star range) the LMC field variables are brighter by about one magnitude than those of similar periods in galactic globular clusters. The long-period cluster stars also differ from both these LMC stars and galactic field RV Tau stars in a colour–colour diagram. The reasons for these differences are discussed.     

A robust statistical estimation of the basic parameters of single stellar populations – I. Method

The colour–magnitude diagrams of resolved single stellar populations, such as open and globular clusters, have provided the best natural laboratories to test stellar evolution theory. Whilst a variety of techniques have been used to infer the basic properties of these simple populations, systematic uncertainties arise from the purely geometrical degeneracy produced by the similar shape of isochrones of different ages and metallicities. Here we present an objective and robust statistical technique which lifts this degeneracy to a great extent through the use of a key observable: the number of stars along the isochrone. Through extensive Monte Carlo simulations we show that, for instance, we can infer the four main parameters (age, metallicity, distance and reddening) in an objective way, along with robust confidence intervals and their full covariance matrix. We show that systematic uncertainties due to field contamination, unresolved binaries, initial or present-day stellar mass function are either negligible or well under control. This technique provides, for the first time, a proper way to infer with unprecedented accuracy the fundamental properties of simple stellar populations, in an easy-to-implement algorithm.     

BVI CCD photometry of the distant open star clusters Berkeley 81, Berkeley 99, NGC 6603 and NGC 7044

We present CCD observations for the distant northern open star clusters Berkeley 81, Berkeley 99, NGC 6603 and NGC 7044 in B V I photometric passbands. A total of 9 900 stars have been observed in fields of about 6 × 6 arcmin² of the sky around the clusters. Colour–magnitude diagrams in V , ( B  −  V ) and V , ( V  −  I ) have been generated down to V  = 22 mag and, for the first time, such diagrams have been produced for the clusters Berkeley 81 and Berkeley 99. The data serve as a base for the study of mass functions and for comparison with theoretical models. Analysis of the radial distribution of stellar surface density indicates that the radius values for Berkeley 81, Berkeley 99, NGC 6603 and NGC 7044 are 2.7, 2.8, 2.8 and 2.2 arcmin respectively. By fitting the latest convective core overshooting isochrones to the colour–magnitude diagram and using its morphological features, reddenings, distances and ages of the star clusters have been determined. Broad but well-defined main sequences with stellar evolutionary effects in the brighter stars are clearly visible in colour–magnitude diagrams of all the clusters under study. Some blue stragglers along with well-developed giant branches and red giant clumps are also clearly seen in all of them. The clusters studied here are located at a distance of ∼ 3 kpc, except for Berkeley 99 which is located at a distance of 4.9 kpc. Their linear sizes lie between 3.8 and 8.0 pc; E ( B  −  V ) values range from 0.3 to 1.0 mag, while their ages are between 0.5 and 3.2 Gyr. Thus the star clusters studied here are of intermediate and high age but are compact and distant objects.     

A search for planets transiting the M-dwarf debris disc host, AU Microscopii

We present high-cadence, high-precision multiband photometry of the young, M1Ve, debris disc star, AU Microscopii. The data were obtained in three continuum filters spanning a wavelength range from 4500 to 6600 Å, plus Hα, over 28 nights in 2005. The light curves show intrinsic stellar variability due to star-spots with an amplitude in the blue band of 0.051 mag and a period of 4.847 d. In addition, three large flares were detected in the data which all occur near the minimum brightness of the star. We remove the intrinsic stellar variability and combine the light curves of all the filters in order to search for transits by possible planetary companions orbiting in the plane of the nearly edge-on debris disc. The combined final light curve has a sampling of 0.35 min and a standard deviation of 6.8 mmag. We performed Monte Carlo simulations by adding fake transits to the observed light curve and find with 95 per cent significance that there are no Jupiter mass planets orbiting in the plane of the debris disc on circular orbits with periods,   P ≤ 5  d. In addition, there are no young Neptune like planets (with radii 2.5 times smaller than the young Jupiter) on circular orbits with periods,   P ≤ 3  d.     

The Monitor project: searching for occultations in young open clusters

The Monitor project is a photometric monitoring survey of nine young (1–200 Myr) clusters in the solar neighbourhood to search for eclipses by very low mass stars and brown dwarfs and for planetary transits in the light curves of cluster members. It began in the autumn of 2004 and uses several 2- to 4-m telescopes worldwide. We aim to calibrate the relation between age, mass, radius and where possible luminosity, from the K dwarf to the planet regime, in an age range where constraints on evolutionary models are currently very scarce. Any detection of an exoplanet in one of our youngest targets (≲10 Myr) would also provide important constraints on planet formation and migration time-scales and their relation to protoplanetary disc lifetimes. Finally, we will use the light curves of cluster members to study rotation and flaring in low-mass pre-main-sequence stars.
The present paper details the motivation, science goals and observing strategy of the survey. We present a method to estimate the sensitivity and number of detections expected in each cluster, using a simple semi-analytic approach which takes into account the characteristics of the cluster and photometric observations, using (tunable) best-guess assumptions for the incidence and parameter distribution of putative companions, and we incorporate the limits imposed by radial velocity follow-up from medium and large telescopes. We use these calculations to show that the survey as a whole can be expected to detect over 100 young low and very low mass eclipsing binaries, and ∼3 transiting planets with radial velocity signatures detectable with currently available facilities.