On the intrinsic bias in detecting caustic crossings between Galactic halo and self-lensing events in the Magellanic Clouds

In this paper, we investigate the intrinsic bias in detecting caustic crossings between the Galactic halo and self-lensing gravitational microlensing events in the Magellanic Clouds. For this, we determine the region for optimal caustic-crossing detection in the parameter space of the physical binary separations, ℓ, and the total binary lens mass, M , and find that the optimal regions for both populations of events are similar to each other. In particular, if the Galactic halo is composed of lenses with the claimed average mass of 〈 M 〉∼0.5 M_⊙, the optimal binary separation range of Galactic halo events of 3.5 au≲ℓ≲14 au matches well with that of a Magellanic Cloud self-lensing event caused by a binary lens with a total mass of M ∼1 M_⊙; well within the mass range of the most probable lens population of stars in the Magellanic Clouds. Therefore, our computation implies that if the binary fractions and the distributions of binary separations of the two populations of lenses are not significantly different from each other, there is no strong detection bias against Galactic halo caustic-crossing events.     

The origin of the magnetic fields in white dwarfs

Magnetic white dwarfs with fields in excess of ∼10⁶ G (the high field magnetic white dwarfs; HFMWDs) constitute about ∼10 per cent of all white dwarfs and show a mass distribution with a mean mass of  ∼0.93 M_⊙  compared to  ∼0.56 M_⊙  for all white dwarfs. We investigate two possible explanations for these observations. First, that the initial–final mass relationship (IFMR) is influenced by the presence of a magnetic field and that the observed HFMWDs originate from stars on the main sequence that are recognized as magnetic (the chemically peculiar A and B stars). Secondly, that the IFMR is essentially unaffected by the presence of a magnetic field, and that the observed HFMWDs have progenitors that are not restricted to these groups of stars. Our calculations argue against the former hypothesis and support the latter. The HFMWDs have a higher than average mass because on the average they have more massive progenitors and not because the IFMR is significantly affected by the magnetic field. A requirement of our model is that ∼40 per cent of main-sequence stars more massive than  ∼4.5 M_⊙  must either have magnetic fields in the range of ∼10–100 G, which is below the current level of detection, or generate fields during subsequent stellar evolution towards the white dwarf phase. In the former case, the magnetic fields of the HFMWDs could be fossil remnants from the main-sequence phase consistent with the approximate magnetic flux conservation.     

Constraints on Jupiters from observations of Galactic bulge microlensing events during 2000

We present observations of eight Galactic bulge microlensing events taken with the 1.0-m Jacobus Kapteyn Telescope (JKT) on La Palma during 2000 June and July. The JKT observing schedule was optimized using a prioritizing algorithm to automatically update the target list. For most of these events we have sampled the light curves at times where no information was available from the OGLE alert team. We assume a point-source point-lens (PSPL) model and perform a maximum likelihood fit to both our data and the OGLE data to constrain the event parameters of the fit. We then refit the data assuming a binary lens and proceed to calculate the probability of detecting planets with mass ratio   q = 10⁻³  . We have seen no clear signatures of planetary deviations on any of the eight events and we quantify constraints on the presence of planetary companions to the lensing stars. For two well-observed events, 2000BUL31 and 2000BUL33, our detection probabilities peak at ∼30 and ∼20 per cent respectively for   q = 10⁻³  and   a ∼ R _E  for a  Δχ²  threshold value of 60.     

The potential for Earth-mass planet formation around brown dwarfs

Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the substellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low-mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarf mass 0.05 M_⊙, our models predict a maximum planetary mass of  ∼5   M_⊕  , orbiting with semimajor axis ∼ 1 au. However, we note that the predictions for the mass–semimajor axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the substellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasize the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.     

A source of high-velocity white dwarfs

We investigate whether the recently observed population of high-velocity white dwarfs can be derived from a population of binaries residing initially within the thin disc of the Galaxy. In particular, we consider binaries where the primary is sufficiently massive to explode as a Type II supernova. A large fraction of such binaries are broken up when the primary then explodes as a supernova, owing to the combined effects of the mass loss from the primary and the kick received by the neutron star on its formation. For binaries where the primary evolves to fill its Roche lobe, mass transfer from the primary leads to the onset of a common envelope phase during which the secondary and the core of the primary spiral together as the envelope is ejected. Such binaries are the progenitors of X-ray binaries if they are not broken up when the primary explodes. For those systems that are broken up, a large number of the secondaries receive kick velocities ∼100–200 km s⁻¹ and subsequently evolve into white dwarfs. We compute trajectories within the Galactic potential for this population of stars and relate the birth rate of these stars over the entire Galaxy to those seen locally with high velocities relative to the local standard of rest (LSR) . We show that for a reasonable set of assumptions concerning the Galactic supernova rate and the binary population, our model produces a local number density of high-velocity white dwarfs compatible with that inferred from observations. We therefore propose that a population of white dwarfs originating in the thin disc may make a significant contribution to the observed population of high-velocity white dwarfs.     

Magnetic fields and rotation in white dwarfs and neutron stars

Recent spectropolarimetric observations of Ap and Bp stars with improved sensitivity have suggested that most Ap and Bp stars are magnetic with dipolar fields of at least a few hundred gauss. These new estimates suggest that the range of magnetic fluxes found for the majority of magnetic white dwarfs is similar to that of main-sequence Ap–Bp stars, thus strengthening the empirical evidence for an evolutionary link between magnetism on the main sequence and magnetism in white dwarfs. We draw parallels between the magnetic white dwarfs and the magnetic neutron stars and argue that the observed range of magnetic fields in isolated neutron stars  ( B_p ∼ 10¹¹–10¹⁵ G)  could also be explained if their mainly O-type progenitors have effective dipolar fields in the range of a few gauss to a few kilogauss, assuming approximate magnetic flux conservation with the upper limit being consistent with the recent measurement of a field of   B_p ∼ 1100 G  for θ Orion C.
In the magnetic field–rotation diagram, the magnetic white dwarfs can be divided into three groups of different origin: a significant group of strongly magnetized slow rotators  ( P _rot∼ 50 –100 yr)  that have originated from single-star evolution, a group of strongly magnetized fast rotators  ( P _rot∼ 700 s)  , typified by EUVE J0317–853, that have originated from a merger, and a group of modest rotators ( P _rot∼ hours–days) of mixed origin (single-star and CV-type binary evolution). We propose that the neutron stars may similarly divide into distinct classes at birth , and suggest that the magnetars may be the counterparts of the slowly rotating high-field magnetic white dwarfs.     

Simulating Gaia performances on white dwarfs

One of the most promising space missions of the European Space Agency is the astrometric satellite Gaia , which will provide very precise astrometry and multicolour photometry, for all 1.3 billion objects to   V ∼ 20  , and radial velocities with accuracies of a few km s⁻¹ for most stars brighter than   V ∼ 17  . Consequently, full homogeneous six-dimensional phase-space information for a huge number of stars will become available. Our Monte Carlo simulator has been used to estimate the number of white dwarfs potentially observable by Gaia . From this we assess the white dwarf luminosity functions that Gaia will obtain and discuss in depth the scientific returns of Gaia in the specific field of white dwarf populations. Scientifically attainable goals include, among others, a reliable determination of the age of the Galactic disc, a better knowledge of the halo of the Milky Way and the reconstruction of the star formation history of the Galactic disc. Our results also demonstrate the potential impact of a mission such as Gaia within the context of current understanding of white dwarf cooling theory.     

Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.     

Distribution of caustic-crossing intervals for Galactic binary-lens microlensing events

Detection of caustic crossings of binary-lens gravitational microlensing events is important because by detecting them one can obtain useful information about both the lens and the source star. In this paper, we compute the distribution of the intervals between two successive caustic crossings, f ( t _cc), for Galactic bulge binary-lens events to investigate the observational strategy for the optimal detection and resolution of caustic crossings. From this computation, we find that the distribution is highly skewed towards short t _cc and peaks at t _cc∼1.5 d. For the maximal detection of caustic crossings, therefore, prompt initiation of follow-up observations for intensive monitoring of events will be important. We estimate that, under the strategy of the current follow-up observations with a second caustic-crossing preparation time of ∼2 d, the fraction of events with resolvable caustic crossing is ∼80 per cent. We find that if the follow-up observations can be initiated within 1 d after the first caustic crossing by adopting more aggressive observational strategies, the detection rate can be improved to ∼90 per cent.     

Substellar companions and isolated planetary-mass objects from protostellar disc fragmentation

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M _disc/ M _*= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M_⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.     

Repeating microlensing events in the OGLE data

Microlensing events are usually selected among single-peaked non-repeating light curves in order to avoid confusion with variable stars. However, a microlensing event may exhibit a second microlensing brightening episode when the source or/and the lens is a binary system. A careful analysis of these repeating events provides an independent way to study the statistics of wide binary stars and to detect extrasolar planets. Previous theoretical studies predicted that 0.5–2 per cent of events should repeat due to wide binary lenses. We present a systematic search for such events in about 4000 light curves of microlensing candidates detected by the Optical Gravitational Lensing Experiment (OGLE) towards the Galactic bulge from 1992 to 2007. The search reveals a total of 19 repeating candidates, with six clearly due to a wide binary lens. As a by-product, we find that 64 events (∼2 per cent of the total OGLE-III sample) have been misclassified as microlensing; these misclassified events are mostly nova or other types of eruptive stars. The number and importance of repeating events will increase considerably when the next-generation wide-field microlensing experiments become fully operational in the future.     

L dwarfs in the Hyades

We present the results of a proper motion survey of the Hyades to search for brown dwarfs, based on UKIRT Deep Sky Survey (UKIDSS) and Two-Micron All Sky Survey (2MASS) data. This survey covers  ∼275 deg²  to a depth of   K ∼ 15  mag, equivalent to a mass of  ∼0.05 M_⊙  assuming a cluster age of 625 Myr. The discovery of 12 L dwarf Hyades members is reported. These members are also brown dwarfs, with masses between  0.05 < M < 0.075 M_⊙  . A high proportion of these L dwarfs appear to be photometric binaries.     

Accretion and dynamical interactions in small-N star-forming clusters: N= 5

We present results from high-resolution hydrodynamical simulations that explore the effects of small-scale clustering in star-forming regions. A large ensemble of small- N clusters with five stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data.
Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses that is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed initial mass function, in both the stellar and substellar regimes. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities ≥2 km s⁻¹, which might deprive low-mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long-period (uncircularized) binaries.
The results obtained support a mechanism in which a significant fraction of brown dwarfs form under similar circumstances as those of normal stars but are ejected from the common envelope of unstable multiple systems before their masses exceed the hydrogen burning limit. We predict that many close binary stars should have wide brown dwarf companions. Brown dwarfs, and, in general, very low-mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or substellar primaries being scarce. Where such binaries exist, they are expected either to be close enough (semimajor axis ∼10 au) to survive strong interactions with more massive binaries or to be born in very small molecular cloud cores.     

Something about the structure of the Galaxy

We analyse a sample of 507 evolved (OH/IR) stars in the region (10°>ℓ>−45°), (| b |<3°). We derive average ages for subsets of this sample, and use those sets as beacons for the evolution of the Galaxy. In the bulge, the oldest OH/IR stars in the plane are 7.5 Gyr (1.3 M_⊙), and in the disc 2.7 Gyr (2.3 M_⊙). The vertical distribution of almost all AGB stars in the disc is found to be nearly exponential, with scaleheight increasing from 100 pc for stars ≲1 Gyr old to 500 pc for stars ≳5 Gyr old. There may be a small, disjunct population of OH/IR stars. The radial distribution of AGB stars is dictated by the metallicity gradient. Unequivocal morphological evidence is presented for the existence of a central bar, but parameters can be constrained only for a given spatial-density model. Using a variety of indicators, we identify the radii of the inner ultraharmonic (2.5 kpc) and corotation resonance (3.5 kpc). We show that the 3-kpc arm is likely to be an inner ring, as observed in other barred galaxies, by identifying a group of evolved stars that is connected to the 3-kpc H  i filament. Also, using several observed features, we argue that an inner-Lindblad resonance exists, at ∼1–1.5 kpc. The compositions of OH/IR populations within 1 kpc of the Galactic Centre give insight into the bar-driven evolution of the inner regions. We suggest that the bar is ∼8 Gyr old, relatively weak (SAB), and may be in a final stage of its existence.     

Spatial distribution of stars and brown dwarfs in σ Orionis

I have re-visited the spatial distribution of stars and high-mass brown dwarfs in the σ Orionis (σ Ori) cluster (∼3 Ma, ∼360 pc). The input was a catalogue of 340 cluster members and candidates at separations less than 30 arcmin to σ Ori AB. Of them, 70 per cent have features of extreme youth. I fitted the normalized cumulative number of objects counting from the cluster centre to several power-law, exponential and King radial distributions. The cluster seems to have two components: a dense core that extends from the centre to   r ≈  20 arcmin and a rarified halo at larger separations. The radial distribution in the core follows a power law proportional to r ¹, which corresponds to a volume density proportional to   r ⁻²  . This is consistent with the collapse of an isothermal spherical molecular cloud. The stars more massive than  3.7 M_⊙  concentrate, however, towards the cluster centre, where there is also an apparent deficit of very low mass objects  ( M < 0.16 M_⊙)  . Last, I demonstrated through Monte Carlo simulations that the cluster is azimuthally asymmetric, with a filamentary overdensity of objects that runs from the cluster centre to the Horsehead Nebula.     

Contribution to the search for binaries among Am stars – VIII. New spectroscopic orbits of eight systems and statistical study of a sample of 91 Am stars

Red clump giant (RCG) stars can be used as distance indicators to trace the mass distribution of the Galactic bar. We use RCG stars from 44 bulge fields from the OGLE-II microlensing collaboration data base to constrain analytic triaxial models for the Galactic bar. We find the bar major-axis is oriented at an angle of 24°–27° to the Sun–Galactic Centre line-of-sight. The ratio of semimajor and semiminor bar axis scalelengths in the Galactic plane   x ₀, y ₀  , and vertical bar scalelength z ₀, is   x ₀ :  y ₀ :  z ₀= 10 : 3.5 : 2.6  , suggesting a slightly more prolate bar structure than the working model of Gerhard which gives the scalelength ratios as   x ₀ :  y ₀ :  z ₀= 10 : 4 : 3  .     

Feast and Famine: regulation of black hole growth in low-redshift galaxies

We analyse the observed distribution of Eddington ratios  ( L / L _Edd)  as a function of supermassive black hole mass for a large sample of nearby galaxies drawn from the Sloan Digital Sky Survey. We demonstrate that there are two distinct regimes of black hole growth in nearby galaxies. The first is associated with galaxies with significant star formation [   M _*/star formation rate (SFR) ∼  a Hubble time] in their central kiloparsec regions, and is characterized by a broad lognormal distribution of accretion rates peaked at a few per cent of the Eddington limit. In this regime, the Eddington ratio distribution is independent of the mass of the black hole and shows little dependence on the central stellar population of the galaxy. The second regime is associated with galaxies with old central stellar populations (   M _*/SFR ≫  a Hubble time), and is characterized by a power-law distribution function of Eddington ratios. In this regime, the time-averaged mass accretion rate on to black holes is proportional to the mass of stars in the galaxy bulge, with a constant of proportionality that depends on the mean stellar age of the stars. This result is once again independent of black hole mass. We show that both the slope of the power law and the decrease in the accretion rate on to black holes in old galaxies are consistent with population synthesis model predictions of the decline in stellar mass loss rates as a function of mean stellar age. Our results lead to a very simple picture of black hole growth in the local Universe. If the supply of cold gas in a galaxy bulge is plentiful, the black hole regulates its own growth at a rate that does not further depend on the properties of the interstellar medium. Once the gas runs out, black hole growth is regulated by the rate at which evolved stars lose their mass.     

Another channel to detect close-in binary companions via gravitational microlensing

Gaudi & Gould showed that close companions of remote binary systems can be efficiently detected by using gravitational microlensing via the deviations in the lensing light curves induced by the existence of the lens companions. In this paper, we introduce another channel to detect faint close-in binary companions by using microlensing. This method utilizes a caustic-crossing binary lens event with a source also composed of binary stars, where the companion is a faint star. Detection of the companion is possible because the flux of the companion can be highly amplified when it crosses the lens caustic. The detection is facilitated since the companion is more amplified than the primary because it, in general, has a smaller size than the primary, and thus experiences less finite source effect. The method is an extension of the previous one suggested to detect close-in giant planets by Graff & Gaudi and Lewis & Ibata and further developed by Ashton & Lewis. From the simulations of realistic Galactic bulge events, we find that companions of K-type main-sequence or brighter stars can be efficiently detected from the current type of microlensing follow-up observations by using the proposed method. We also find that compared with the method of detecting lens companions for which the efficiency drops significantly for binaries with separations ≲0.2 of the angular Einstein ring radius, θ _E, the proposed method has an important advantage of being able to detect companions with substantially smaller separations down to ∼     .     

The birthrate of magnetars

Magnetars, neutron stars with ultrastrong magnetic fields  ( B ∼ 10¹⁴−10¹⁵G)  , manifest their exotic nature in the form of soft gamma-ray repeaters and anomalous X-ray pulsars. This study estimates the birthrate of magnetars to be ∼0.22 per century with a Galactic population comprising ∼17 objects. A population synthesis was carried out based on the five anomalous X-ray pulsars detected in the ROSAT All Sky Survey by comparing their number to that of massive OB stars in a well-defined volume. Additionally, the group of seven X-ray dim isolated neutron stars detected in the same survey were found to have a birthrate of ∼2 per century with a Galactic population of ∼22 000 objects.     

Inversion of stellar statistics equation for the Galactic bulge

A method based on Lucy's iterative algorithm is developed to invert the equation of stellar statistics for the Galactic bulge and is then applied to the K -band star counts from the Two-Micron Galactic Survey in a number of off-plane regions (10°>| b |>2°, | l |<15°). The top end of the K -band luminosity function is derived and the morphology of the stellar density function is fitted to triaxial ellipsoids, assuming a non-variable luminosity function within the bulge. The results, which have already been outlined by López-Corredoira et al., are shown in this paper with a full explanation of the steps of the inversion: the luminosity function shows a sharp decrease brighter than M _K =−8.0  mag when compared with the disc population; the bulge fits triaxial ellipsoids with the major axis in the Galactic plane at an angle with the line of sight to the Galactic centre of 12° in the first quadrant; the axial ratios are 1:0.54:0.33, and the distance of the Sun from the centre of the triaxial ellipsoid is 7860 pc. The major–minor axial ratio of the ellipsoids is found not to be constant, the best fit to the gradient being K _z =(8.4±1.7)×exp(− t /(2000±920) pc), where t is the distance along the major axis of the ellipsoid in parsecs. However, the interpretation of this is controversial. An eccentricity of the true density-ellipsoid gradient and a population gradient are two possible explanations. The best fit for the stellar density, for 1300 pc< t <3000 pc, is calculated for both cases, assuming an ellipsoidal distribution with constant axial ratios, and when K _z is allowed to vary. From these, the total number of bulge stars is ∼3×10¹⁰ or ∼4×10¹⁰, respectively.