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.     

The lens and source of the optical Einstein ring gravitational lens ER 0047–2808

We perform a detailed analysis of the optical gravitational lens ER 0047–2808 imaged with the Wide Field Planetary Camera 2 on the Hubble Space Telescope . Using software specifically designed for the analysis of resolved gravitational lens systems, we focus on how the image alone can constrain the mass distribution in the lens galaxy. We find that the data are of sufficient quality to strongly constrain the lens model with no a priori assumptions about the source. Using a variety of mass models, we find statistically acceptable results for elliptical isothermal-like models with an Einstein radius of 1.17 arcsec. An elliptical power-law model  (Σ∝ R ^−β)  for the surface mass density favours a slope slightly steeper than isothermal with  β= 1.08 ± 0.03  . Other models including a constant mass-to-light ratio (M/L), pure Navarro, Frenk & White halo and (surprisingly) an isothermal sphere with external shear are ruled out by the data. We find the galaxy light profile can only be fit with a Sérsic plus point-source model. The resulting total  M/L _B   contained within the images is  4.7  h ₆₅± 0.3  . In addition, we find the luminous matter is aligned with the total mass distribution within a few degrees. This is the first time a resolved optical gravitational lens image has been quantitatively reproduced using a non-parametric source.
The source, reconstructed by the software, is revealed to have two bright regions, with an unresolved component inside the caustic and a resolved component straddling a fold caustic. The angular size of the entire source is ∼0.1 arcsec and its (unlensed) Lyα flux is  3 × 10⁻¹⁷ erg s⁻¹ cm⁻²  .     

Analytic relations for magnifications and time delays in gravitational lenses with fold and cusp configurations

Gravitational lensing provides a unique and powerful probe of the mass distributions of distant galaxies. Four-image lens systems with fold and cusp configurations have two or three bright images near a critical point. Within the framework of singularity theory, we derive analytic relations that are satisfied for a light source that lies a small but finite distance from the astroid caustic of a four-image lens. Using a perturbative expansion of the image positions, we show that the time delay between the close pair of images in a fold lens scales with the cube of the image separation, with a constant of proportionality that depends on a particular third derivative of the lens potential. We also apply our formalism to cusp lenses, where we develop perturbative expressions for the image positions, magnifications and time delays of the images in a cusp triplet. Some of these results were derived previously for a source asymptotically close to a cusp point, but using a simplified form of the lens equation whose validity may be in doubt for sources that lie at astrophysically relevant distances from the caustic. Along with the work of Keeton, Gaudi & Petters, this paper demonstrates that perturbation theory plays an important role in theoretical lensing studies.     

Suzaku observations of Markarian 335: evidence for a distributed reflector

We report on a 151-ks net exposure Suzaku observation of the narrow-line Seyfert 1 galaxy Mrk 335. The 0.5–40 keV spectrum contains a broad Fe line, a strong soft excess below about 2 keV and a Compton hump around 20–30 keV. We find that a model consisting of a power law and two reflectors provides the best fit to the time-averaged spectrum. In this model, an ionized, heavily blurred, inner reflector produces most of the soft excess, while an almost neutral outer reflector (outside ∼ 40 r _g) produces most of the Fe line emission. The spectral variability of the observation is characterized by spectral hardening at very low count rates. In terms of our power-law + two-reflector model it seems like this hardening is mainly caused by pivoting of the power law. The rms spectrum of the entire observation has the curved shape commonly observed in active galactic nuclei, although the shape is significantly flatter when an interval which does not contain any deep dip in the light curve is considered. We also examine a previous 133-ks XMM–Newton observation of Mrk 335. We find that the XMM–Newton spectrum can be fitted with a similar two-reflector model as the Suzaku data and we confirm that the rms spectrum of the observation is flat. The flat rms spectra, as well as the high-energy data from the Suzaku PIN detector, disfavour an absorption origin for the soft excess in Mrk 335.     

Galaxy density profiles and shapes – II. Selection biases in strong lensing surveys

Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio and image separation distribution, with a flexible treatment of magnification bias to mimic different survey strategies. We present our results for two families of density profiles: cusped and deprojected Sérsic models. While we use fixed lens and source redshifts for most of the analysis, we show that the results are applicable to other redshift combinations, and we also explore the physics of how our results change for very different redshifts. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sérsic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.     

Pixel lensing as a way to detect extrasolar planets in M31

We study the possibility to detect extrasolar planets in M31 through pixel-lensing observations. Using a Monte Carlo approach, we select the physical parameters of the binary lens system, a star hosting a planet, and we calculate the pixel-lensing light curve taking into account the finite source effects. Indeed, their inclusion is crucial since the sources in M31 microlensing events are mainly giant stars. Light curves with detectable planetary features are selected by looking for significant deviations from the corresponding Paczyński shapes. We find that the time-scale of planetary deviations in light curves increase (up to 3–4 d) as the source size increases. This means that only few exposures per day, depending also on the required accuracy, may be sufficient to reveal in the light curve a planetary companion. Although the mean planet mass for the selected events is about     , even small mass planets  ( M _P < 20 M_⊕)  can cause significant deviations, at least in the observations with large telescopes. However, even in the former case, the probability to find detectable planetary features in pixel-lensing light curves is at most a few per cent of the detectable events, and therefore many events have to be collected in order to detect an extrasolar planet in M31. Our analysis also supports the claim that the anomaly found in the candidate event PA-99-N2 towards M31 can be explained by a companion object orbiting the lens star.     

Intragroup diffuse light in compact groups of galaxies – II. HCG 15, 35 and 51

This continuing study of intragroup light in compact groups of galaxies aims to establish new constraints to models of formation and evolution of galaxy groups, specially of compact groups, which are a key part in the evolution of larger structures, such as clusters. In this paper we present three additional groups (HCG 15, 35 and 51) using deep wide-field B - and R -band images observed with the LAICA camera at the 3.5-m telescope at the Calar Alto observatory (CAHA). This instrument provides us with very stable flat-fielding, a mandatory condition for reliably measuring intragroup diffuse light. The images were analysed with the ov_wav package, a wavelet technique that allows us to uncover the intragroup component in an unprecedented way. We have detected that 19, 15 and 26 per cent of the total light of HCG 15, 35 and 51, respectively, are in the diffuse component, with colours that are compatible with old stellar populations and with mean surface brightness that can be as low as 28.4 B mag arcsec⁻². Dynamical masses, crossing times and mass-to-light ratios were recalculated using the new group parameters. Also tidal features were analysed using the wavelet technique.     

Stellar disruption by a supermassive black hole: is the light curve really proportional to t−5/3?

In this paper, we revisit the arguments for the basis of the time evolution of the flares expected to arise when a star is disrupted by a supermassive black hole. We present a simple analytic model relating the light curve to the internal density structure of the star. We thus show that the standard light curve proportional to   t ^−5/3  only holds at late times. Close to the peak luminosity the light curve is shallower, deviating more strongly from   t ^−5/3  for more centrally concentrated (e.g. solar type) stars. We test our model numerically by simulating the tidal disruption of several stellar models, described by simple polytropic spheres with index γ. The simulations agree with the analytical model given two considerations. First, the stars are somewhat inflated on reaching pericentre because of the effective reduction of gravity in the tidal field of the black hole. This is well described by a homologous expansion by a factor which becomes smaller as the polytropic index becomes larger. Secondly, for large polytropic indices wings appear in the tails of the energy distribution, indicating that some material is pushed further away from parabolic orbits by shocks in the tidal tails. In all our simulations, the   t ^−5/3  light curve is achieved only at late stages. In particular, we predict that for solar-type stars, this happens only after the luminosity has dropped by at least 2 mag from the peak. We discuss our results in the light of recent observations of flares in otherwise quiescent galaxies and note the dependence of these results on further parameters, such as the star/hole mass ratio and the stellar orbit.     

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.     

Inferring dark halo structure from observed scaling laws of late-type galaxies and LSBs

We re-examine the Fall & Efstathiou scenario for galaxy formation, including the dark halo gravitational reaction to the formation of the baryon disc, as well as continuous variations in the intrinsic halo density profile. The recently published rotation curves of low surface brightness (LSB) and dwarf galaxies together with previously known scaling relations provide sufficient information on the present-day structure of late-type disc galaxies to invert the problem. By requiring that the models reproduce all the observational restrictions we can fully constrain the initial conditions of galaxy formation, with a minimum of assumptions, in particular without the need to specify a cold dark matter (CDM) halo profile. This allows one to solve for all the initial conditions, in terms of the halo density profile, the baryon fraction and the total angular momentum. We find that a unique initial halo shape is sufficient to accurately reproduce the rotation curves of both LSB and normal late-type spiral galaxies. This unique halo profile differs substantially from that found in standard CDM models. A galactic baryon fraction of 0.065 is found. The initial value of the dimensionless angular momentum is seen to be the principal discriminator between the galaxy classes we examine. The present-day scalings between structural parameters are seen to originate in the initial conditions.     

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.     

NICMOS images of JVAS/CLASS gravitational lens systems

We present Hubble Space Telescope ( HST ) infrared images of four gravitational lens systems from the JVAS/CLASS gravitational lens survey and compare the new infrared HST pictures with previously published WFPC2 HST optical images and radio maps. Apart from the wealth of information that we get from the flux ratios and accurate positions and separations of the components of the lens systems, which we can use as inputs for better constraints on the lens models, we are able to discriminate between reddening and optical/radio microlensing as the possible cause of differences observed in the flux ratios of the components across the three wavelength bands. Substantial reddening has been known to be present in the lens system B1600+434 and has been further confirmed by the present infrared data. In the two systems B0712+472 and B1030+074 microlensing has been pinpointed as the main cause of the flux ratio discrepancy both in the optical/infrared and in the radio, the radio possibly caused by the substructure revealed in the lensing galaxies. In B0218+357, however, the results are still not conclusive. If we are actually seeing the two 'true' components of the lens system then the flux ratio differences are attributed to a combination of microlensing and reddening or are alternatively the result of some variability in at least one of the images. Otherwise the second 'true' component of B0218+357 may be completely absorbed by a molecular cloud and the anomalous flux density ratios and large difference in separation between the optical/infrared and radio that we see can be explained by emission either from a foreground object or from part of the lensing galaxy.     

The flux-dependent amplitude of broadband noise variability in X-ray binaries and active galaxies

Standard shot-noise models, which seek to explain the broadband noise variability that characterizes the X-ray light curves of X-ray binaries and active galaxies, predict that the power spectrum of the X-ray light curve is stationary (i.e. constant amplitude and shape) on short time-scales. We show that the broadband noise power spectra of the black hole candidate Cyg X-1 and the accreting millisecond pulsar SAX J1808.4−3658 are intrinsically non-stationary, in that rms variability scales linearly with flux. Flux-selected power spectra confirm that this effect is due to changes in power-spectral amplitude and not shape. The light curves of three Seyfert galaxies are also consistent with a linear relationship between rms variability and flux, suggesting that it is an intrinsic feature of the broadband noise variability in compact accreting systems over more than six decades of central object mass. The rms variability responds to flux variations on all measured time-scales, raising fundamental difficulties for shot-noise models which seek to explain this result by invoking variations in the shot parameters. We suggest that models should be explored where the longest time-scale variations are fundamental and precede the variations on shorter time-scales. Possible models which can explain the linear rms-flux relation include the fractal break-up of large coronal flares, or the propagation of fluctuations in mass accretion rate through the accretion disc. The linear relationship between rms variability and flux in Cyg X-1 and SAX J1808.4−3658 is offset on the flux axis, suggesting the presence of a second, constant-flux component to the light curve which contributes ∼25 per cent of the total flux. The spectrum of this constant component is similar to the total spectrum, suggesting that it may correspond to quiet, non-varying regions in the X-ray emitting corona.     

Modelling the first probable two-plane lens system B2114+022: reproducing two compact radio cores A and D

We test possible lensing scenarios of the JVAS system B2114+022, in which two galaxies at different redshifts ('G1' at z ₁=0.3157 and 'G2' at z ₂=0.5883) are found within 2 arcsec of quadruple radio sources. For our investigation, we use possible lensing constraints derived from a wealth of data on the radio sources obtained with VLA, MERLIN, VLBA and EVN as well as HST imaging data on the two galaxies, which were presented recently in Augusto et al. In the present study, we focus on reproducing the widest separated, observationally similar radio components A and D as lensed images. We first treat G2 (which is the more distant one from the geometric centre) as a shear term, and then consider two-plane lensing explicitly including G2's potential at the z ₂ plane as the first case of two-plane lens modelling. Our modelling results not only support the hypothesis that the system includes gravitationally lensed images of a higher-redshift extragalactic object, but they also show that the explicit inclusion of G2's potential at the second lens plane is necessary in order to fit the data with astrophysically plausible galaxy parameters. Finally, we illustrate a natural consequence of a two-plane lens system, namely the prediction of distortion as well as shift and stretching of G2's isophotes by G1's potential, which can in principle be measured by subtracting out G1's light distribution in an image of high signal-to-noise ratio and good angular resolution, especially a multicolour one.     

Long-term optical and X-ray variability of the Seyfert galaxy Markarian 79

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than 5 yr. We find that on short to medium time-scales (days to a few tens of days) the 2–10 keV X-ray and optical u - and V -band fluxes are significantly correlated, with a delay between the bands consistent with 0 d. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.     

The anatomy of a quadruply imaged gravitational lens system

The key to using a strong gravitational lens system to measure the Hubble constant is to obtain an accurate model of the lens potential. In this paper, we investigate the properties of gravitational lens B1608+656, a quadruply imaged lens system with an extended source intensity distribution. Our analysis is valid for generic quadruply lensed systems. Limit curves and isophotal separatrices are defined for such systems, and we show that the isophotal separatrices must intersect at the critical curves and the satellite isophotes must be tangent to the limit curves. The most recent model of B1608+656 by Koopmans et al. satisfies these criteria for some, but not all, of the isophotal separatrices within the observational uncertainty. We study a non-parametric method of potential reconstruction proposed by Blandford, Surpi & Kundic and demonstrate that although the method works in principle and elucidates image formation, the initial potential only converges to the true model when it is within ∼1 per cent of the true model.     

Variation of the broad X-ray iron line in MCG–6-30-15 during a flare

We report results on the broad iron emission line of the Seyfert galaxy MCG–6-30-15, obtained from the second long ASCA observation in 1997. The time-averaged profile of the broad line is very similar to that seen with ASCA in 1994, so confirming the detailed model fit then obtained. A bright flare is seen in the light curve, during which the continuum was soft. At that time the emission line peaks around 5 keV and most of its emission is shifted below 6 keV with no component detected at 6.4 keV  (  This can be interpreted as the result of an extraordinarily large gravitational redshift owing to a dominant flare occurring very close to the black hole at a radius of  ≲5  .     

Separation of the visible and dark matter in the Einstein ring LBG J213512.73−010143

We model the mass distribution in the recently discovered Einstein ring LBG J213512.73−010143 (the 'Cosmic Eye') using archival Hubble Space Telescope imaging. We reconstruct the mass density profile of the z = 0.73 lens and the surface brightness distribution of the z = 3.07 source and find that the observed ring is best fitted with a dual-component lens model consisting of a baryonic Sersic component nested within a dark matter halo. The dark matter halo has an inner slope of 1.42^+0.24_−0.22, consistent with cold dark matter simulations after allowing for baryon contraction. The baryonic component has a mass-to-light ratio of  1.71^+0.28_−0.38 M_⊙/L _{B ⊙}  which when evolved to the present day is in agreement with local ellipticals. Within the Einstein radius of 0.77 arcsec (5.6 kpc), the baryons account for 46 ± 11 per cent of the projected lens mass. External shear from a nearby foreground cluster is accurately predicted by the model. The reconstructed surface brightness distribution in the source plane clearly shows two peaks. Through a generalization of our lens inversion method, we conclude that the redshifts of both peaks are consistent with each other, suggesting that we are seeing structure within a single galaxy.     

B0712+472: a new radio four-image gravitational lens

A new four-image gravitational lens system, B0712+472, has been discovered during the Cosmic Lens All-Sky Survey. This system consists of four flat-spectrum radio images that are also seen on a Hubble Space Telescope ( HST ) image, together with the lensing galaxy. We present MERLIN, VLA and VLBA maps and WHT spectra of the system as well as the HST images. The light distribution of the lensing galaxy is highly elongated and so too is the mass distribution deduced from modelling. We suggest a redshift of ∼1.33 for the lensed object; the lens redshift will require further investigation. The discovery of this new system further increases the ratio of four-image to two-image lens systems currently known, exacerbating problems of required ellipticity of matter distributions in lensing galaxies.