Detection of stellar spots from the observations of caustic-crossing binary-lens gravitational microlensing events

Recently, Heyrovský & Sasselov investigated the sensitivity of single-lens gravitational microlensing event light curves to spots and found that, during source transit, spots can cause deviations in amplification larger than 2 per cent, and thus be detectable. In this paper, we explore the feasibility of spot detection from the observations of binary-lens microlensing events instead of single-lens events. For this we investigate the sensitivity of binary-lens event light curves to spots and compare it with that of single-lens events. From this investigation, we find that during caustic crossings the fractional amplification deviations of light curves from those of spotless source events are equivalent to those of single-lens events, implying that spots can also be detected with a similar photometric precision to that required for spot detection by observing single-lens events. We discuss the relative advantages of observing binary-lens events over the observations of single-lens events in detecting stellar spots.     

A systematic fitting scheme for caustic-crossing microlensing events

We outline a method for fitting binary-lens caustic-crossing microlensing events based on the alternative model parametrization proposed and detailed by Cassan. As an illustration of our methodology, we present an analysis of OGLE-2007-BLG-472, a double-peaked Galactic microlensing event with a source crossing the whole caustic structure in less than three days. In order to identify all possible models we conduct an extensive search of the parameter space, followed by a refinement of the parameters with a Markov Chain Monte Carlo algorithm. We find a number of low-  χ²  regions in the parameter space, which lead to several distinct competitive best models. We examine the parameters for each of them, and estimate their physical properties. We find that our fitting strategy locates several minima that are difficult to find with other modelling strategies and is therefore a more appropriate method to fit this type of event.     

On the baseline flux determination of microlensing events detectable with the difference image analysis method

To improve photometric precision by removing the blending effect, a newly developed technique of difference image analysis (DIA) has been adopted by several gravitational microlensing experiment groups. However, the principal problem of the DIA method is that, by its nature, it has difficulties in measuring the baseline flux F ₀ of a source star, causing a degeneracy problem in determining the lensing parameters of an event. Therefore, it is often believed that the DIA method is not as powerful as the classical method based on PSF photometry for determining the Einstein time-scales t _E of events.
In this paper, we demonstrate that the degeneracy problem in microlensing events, detectable from searches using the DIA method, is not as serious as is often thought. This is because a substantial fraction of events will be high amplification events for which the deviations of the amplification curves, constructed with the wrong baseline fluxes from their corresponding best-fit standard amplification curves, will be considerable, even for a small amount of the fractional baseline flux deviation Δ F ₀ F ₀. With a model luminosity function of source stars and under realistic observational conditions, we find that ∼30 per cent of detectable Galactic bulge events are expected to have high amplifications and their baseline fluxes can be determined with uncertainties Δ F ₀ F ₀≤0.5.     

The brown dwarf mass function from pixel microlensing

We argue that gravitational microlensing is a feasible technique for measuring the mass function of brown dwarf stars in distant galaxies. Microlensing surveys of the bulge of M31, and of M87 in the Virgo cluster, may provide enough events to differentiate the behaviour of the mass function of lenses below the hydrogen-burning limit (although we find that M87 is a more favourable target). Such objects may provide a significant supply of baryonic dark matter, an interesting possibility for the study of galactic dynamics. Furthermore, these systems have different metallicities from the solar neighbourhood, which may affect the mass function. These considerations are relevant in the context of star formation studies.     

Variation of spot-induced anomalies in caustic-crossing binary microlensing event light curves

We investigate the pattern of anomalies in the light curves of caustic-crossing binary microlensing events induced by spot(s) on the lensed source star. To this end, we perform simulations of events with various models of spots. From these simulations we find that the spot-induced anomalies take various forms depending on the physical state of spots, which is characterized by the surface brightness contrast, the size, the number, the umbra/penumbra structure, the shape and the orientation with respect to the sweeping caustic. We also examine the feasibility of distinguishing the two possibly degenerate types of anomalies caused by a spot and a transiting planet and find that in many cases the degeneracy can be separated from the characteristic multiple deviation feature in the spot-induced anomaly pattern caused by the multiplicity of spots.     

Predicting the second caustic crossing in binary microlensing events

We fit binary lens models to the data covering the initial part of real microlensing events in an attempt to predict the time of the second caustic crossing. We use approximations during the initial search through the parameter space for light curves that roughly match the observed ones. Exact methods for calculating the lens magnification of an extended source are used when we refine our best initial models. Our calculations show that the reliable prediction of the second crossing can only be made very late, when the light curve has risen appreciably after the minimum between the two caustic crossings. The best observational strategy is therefore to sample as frequently as possible once the light curve starts to rise after the minimum.     

On the astrometric behaviour of binary microlensing events

Despite the suspected binarity for a significant fraction of Galactic lenses, the current photometric surveys detected binary microlensing events only for a small fraction of the total events. The detection efficiency is especially low for non-caustic crossing events, which comprise the majority of the binary lensing events, as a result of the absence of distinctive features in their light curves combined with small deviations from the standard light curve of a single point-mass event. In addition, even if they are detected, it will be difficult to determine the solution of the binary lens parameters owing to the severe degeneracy problem. In this paper, we investigate the properties of binary lensing event expected when they are astrometrically observed by using high-precision interferometers. For this, we construct vector field maps of excess centroid shifts, which represent the deviations of the binary lensing centroid shifts from those of a single lensing event as a function of source position. From the analysis of the maps, we find that the excess centroid shifts are substantial in a considerably large area around caustics. In addition, they have characteristic sizes and directions depending strongly on the source positions with respect to the caustics and the resulting trajectories of the light centroid (astrometric trajectories) have distinctive features, which can be distinguished from the deviations caused by other reasons. We classify the types of the deviations and investigate where they occur. Because of the strong dependence of the centroid shifts on the lens system geometry combined with the distinctive features in the observed astrometric trajectories, astrometric binary lensing observations will provide an important tool that can probe the properties of the Galactic binary lens population.     

The microlensing rate and mass function versus dynamics of the Galactic bar

With the steady increase of the sample size of observed microlenses towards the central regions of the Galaxy, the main source of the uncertainty in the lens mass will shift from the simple Poisson noise to the intrinsic non-uniqueness of our dynamical models of the inner Galaxy, particularly the Galactic bar. We use a set of simple self-consistent bar models to investigate how the microlensing event rate varies as a function of axis ratio, bar angle and velocity distribution. The non-uniqueness of the velocity distribution of the bar model adds a significant uncertainty (by about a factor of 1.5) to any prediction of the lens mass. Kinematic data and self-consistent models are critical to lift the non-uniqueness. We discuss the implications of these results for the interpretation of microlensing observations of the Galactic bulge. In particular we show that Freeman bar models scaled to the mass of the Galactic bulge/bar imply a typical lens mass of around 0.8 M⊙, a factor of 3–5 times larger than the value from other models.     

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.     

Giant cluster arcs as a constraint on the scattering cross-section of dark matter

We carry out ray tracing through five high-resolution simulations of a galaxy cluster, to study how its ability to produce giant gravitationally lensed arcs is influenced by the collision cross-section of its dark matter. In three cases typical dark matter particles in the cluster core undergo between 1 and 100 collisions per Hubble time; two more explore the long ('collisionless') and short ('fluid') mean free path limits. We study the size and shape distributions of arcs and compute the cross-section for producing 'extreme' arcs of various sizes. Even a few collisions per particle modifies the core structure enough to destroy the ability of the cluster to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to unrealistically large masses before it regains the ability to produce giant arcs. None of our models with self-interacting dark matter (except the 'fluid' limit) is able to produce radial arcs; even the case with the smallest scattering cross-section must be scaled to the upper limit of observed cluster masses before it produces radial arcs. Apparently the elastic collision cross-section of dark matter in clusters must be very small, below 0.1 cm² g⁻¹, to be compatible with the observed ability of clusters to produce both radial arcs and giant arcs.     

Effects of cluster galaxies on arc statistics

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics.
We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N -body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White density profile.
We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in the presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12 arcsec, while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4 arcsec.     

The non-Gaussian tail of cosmic-shear statistics

Owing to gravitational instability, an initially Gaussian density field develops non-Gaussian features as the Universe evolves. The most prominent non-Gaussian features are massive haloes, visible as clusters of galaxies. The distortion of high-redshift galaxy images because of the tidal gravitational field of the large-scale matter distribution, called cosmic shear, can be used to investigate the statistical properties of the large‐scale structure (LSS) . In particular, non-Gaussian properties of the LSS will lead to a non-Gaussian distribution of cosmic-shear statistic. The aperture mass ( M _ap) statistics, recently introduced as a measure for cosmic shear, is particularly well suited for measuring these non-Gaussian properties. In this paper we calculate the highly non-Gaussian tail of the aperture mass probability distribution, assuming Press–Schechter theory for the halo abundance and the 'universal' density profile of haloes as obtained from numerical simulations. We find that for values of M _ap much larger than its dispersion, this probability distribution is closely approximated by an exponential, rather than a Gaussian. We determine the amplitude and shape of this exponential for various cosmological models and aperture sizes, and show that wide-field imaging surveys can be used to distinguish between some of the currently most popular cosmogonies. Our study here is complementary to earlier cosmic-shear investigations, which focused more on two- and three-point statistical properties.