Microlensing of Sub-parsec Massive Binary Black Holes in Lensed QSOs: Light Curves and Size-Wavelength Relation

Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circum-binary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared-optical-UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius vs. wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius-wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Large Synoptic Survey telescope (LSST) and Euclid.Comment: 18 pages, 17 figures, accepted for publication in the Astrophysical Journa

Object DUO 2: A New Binary Lens Candidate

We present the light curve of an unusual variable object, DUO 2, detected during the search for microlensing events by the DUO project. The star remained stable for more than 150 days before it brightened by more than two magnitudes in 6 days in the B and R bands. The light curves are achromatic during the variability. We consider possible explanations of the photometric behavior, with particular emphasis on the binary lens interpretation of the event. The masses of the lenses are quite small, with the companion possibly in the range of a brown dwarf or even a few times of Jupiter. We report evidence of blending of the source by a companion through the first detection of shift in the light centroid among all the microlensing experiments. This shift sets a lower limit of $0.3^{\prime\prime}$ on the separation between the stars. The best lens model obtained requires moderate blending, which was what motivated us to check the centroid shift that was subsequently found. The best lens model predicts a separation of $1^{\prime\prime}$ between the two blended stars. This prediction was recently tested using two CCD images taken under good seeing conditions. Both images show two components. Their separation and position angle are in good agreement with our model.Comment: uuencoded, compressed PostScript, 4 pages, 4 figures (in text). Accepted for publication in Astronomy and Astrophysics Letter

A New Photometric Model of the Galactic Bar using Red Clump Giants

We present a study of the luminosity density distribution of the Galactic bar using number counts of red clump giants (RCGs) from the OGLE-III survey. The data were recently published by Nataf et al. (2013) for 9019 fields towards the bulge and have $2.94\times 10^6$ RC stars over a viewing area of $90.25 \,\textrm{deg}^2$. The data include the number counts, mean distance modulus ($\mu$), dispersion in $\mu$ and full error matrix, from which we fit the data with several tri-axial parametric models. We use the Markov Chain Monte Carlo (MCMC) method to explore the parameter space and find that the best-fit model is the $E_3$ model, with the distance to the GC is 8.13 kpc, the ratio of semi-major and semi-minor bar axis scale lengths in the Galactic plane $x_{0},y_{0}$, and vertical bar scale length $z_0$, is $x_0:y_0:z_0 \approx 1.00:0.43:0.40$ (close to being prolate). The scale length of the stellar density profile along the bar's major axis is $\sim$ 0.67 kpc and has an angle of $29.4^\circ$, slightly larger than the value obtained from a similar study based on OGLE-II data. The number of estimated RC stars within the field of view is $2.78 \times 10^6$, which is systematically lower than the observed value. We subtract the smooth parametric model from the observed counts and find that the residuals are consistent with the presence of an X-shaped structure in the Galactic centre, the excess to the estimated mass content is $\sim 5.8$. We estimate the total mass of the bar is $\sim 1.8 \times 10^{10} M_\odot$. Our results can be used as a key ingredient to construct new density models of the Milky Way and will have implications on the predictions of the optical depth to gravitational microlensing and the patterns of hydrodynamical gas flow in the Milky Way.Comment: 15 pages, 6 figures, 4 tables. MNRAS accepte

The Structure and Clustering of Lyman Break Galaxies

The number density and clustering properties of Lyman-break galaxies (LBGs) are consistent with them being the central galaxies of the most massive dark halos present at z~3. This conclusion holds in all currently popular hierarchical models for structure formation, and is almost independent of the global cosmological parameters. We examine whether the sizes, luminosities, kinematics and star-formation rates of LBGs are also consistent with this identification. Simple formation models tuned to give good fits to low redshift galaxies can predict the distribution of these quantities in the LBG population. The LBGs should be small (with typical half-light radii of 0.6-2 kpc/h), should inhabit haloes of moderately high circular velocity (180-290 km/s) but have low stellar velocity dispersions (70-120 km/s) and should have substantial star formation rates (15-100 Msun/yr). The numbers here refer to the predicted median values in the LBG sample of Adelberger et al. (1998); the first assumes an Omega=1 universe and the second a flat universe with Omega=0.3. For either cosmology these predictions are consistent with the current (rather limited) observational data. Following the work of Kennicutt (1998) we assume stars to form more rapidly in gas of higher surface density. This predicts that LBG samples should preferentially contain objects with low angular momentum, and so small size, for their mass. In contrast, samples of damped Lyman alpha systems (DLSs), should be biased towards objects with large angular momentum. Bright LBGs and DLSs may therefore form distinct populations, with very different sizes and star formation rates, LBGs being smaller and more metal-rich than DLSs of similar mass and redshift.Comment: 27 pages, 9 figures, MNRAS submitte

Mercury deposition in southern New Hampshire, 2006–2009

The atmospheric deposition of mercury (Hg) occurs via several mechanisms including wet, dry, and occult processes. In an effort to understand the atmospheric cycling and seasonal depositional characteristics of Hg, event-based wet deposition samples and reactive gaseous Hg (RGM) measurements were collected for approximately 3 years at Thompson Farm (TF), a near-coastal rural site in Durham, NH, part of the University of New Hampshire AIRMAP Observing Network. Total aqueous mercury exhibited seasonal patterns in Hg wet deposition at TF. The lowest Hg wet deposition was measured in the winter with an average total seasonal deposition of 1.56 μg m−2compared to the summer average of 4.71 μg m−2. Inter-annual differences in total wet deposition are generally linked with precipitation volume, with the greatest deposition occurring in the wettest year. Relationships between surface level RGM and Hg wet deposition were also investigated based on continuous RGM measurements at TF from November 2006 to September 2009. No correlations were observed between RGM mixing ratios and Hg wet deposition, however the ineffective scavenging of RGM during winter precipitation events was evidenced by the less frequent depletion of RGM below the detection level. Seasonal dry deposition of reactive gaseous Hg (RGM) was estimated using an order-of-magnitude approach. RGM mixing ratios and dry deposition estimates were greatest during the winter and spring. The seasonal ratios of Hg wet deposition to RGM dry deposition vary by up to a factor of 80

Solvable senescence model with positive mutations

We build upon our previous analytical results for the Penna model of senescence to include positive mutations. We investigate whether a small but non-zero positive mutation rate gives qualitatively different results to the traditional Penna model in which no positive mutations are considered. We find that the high-lifespan tail of the distribution is radically changed in structure, but that there is not much effect on the bulk of the population. Th e mortality plateau that we found previously for a stochastic generalization of the Penna model is stable to a small positive mutation rate.Comment: 3 figure

Discrete Razumikhin-type technique and stability of the Euler-Maruyama method to stochastic functional differential equations

A discrete stochastic Razumikhin-type theorem is established to investigate whether the Euler--Maruyama (EM) scheme can reproduce the moment exponential stability of exact solutions of stochastic functional differential equations (SFDEs). In addition, the Chebyshev inequality and the Borel-Cantelli lemma are applied to show the almost sure stability of the EM approximate solutions of SFDEs. To show our idea clearly, these results are used to discuss stability of numerical solutions of two classes of special SFDEs, including stochastic delay differential equations (SDDEs) with variable delay and stochastically perturbed equations