Source author record

S. Mao

S. Mao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

MOA-2019-BLG-008Lb: a new microlensing detection of an object at the planet/brown dwarf boundary

We report on the observations, analysis and interpretation of the microlensing event MOA-2019- BLG-008. The observed anomaly in the photometric light curve is best described through a binary lens model. In this model, the source did not cross caustics and no finite source effects were observed. Therefore the angular Einstein ring radius cannot be measured from the light curve alone. However, the large event duration, t E about 80 days, allows a precise measurement of the microlensing parallax. In addition to the constraints on the angular radius and the apparent brightness I s of the source, we employ the Besancon and GalMod galactic models to estimate the physical properties of the lens. We find excellent agreement between the predictions of the two Galactic models: the companion is likely a resident of the brown dwarf desert with a mass Mp about 30 MJup and the host is a main sequence dwarf star. The lens lies along the line of sight to the Galactic Bulge, at a distance of less then4 kpc. We estimate that in about 10 years, the lens and source will be separated by 55 mas, and it will be possible to confirm the exact nature of the lensing system by using high-resolution imaging from ground or space-based observatories.

preprint2022arXiv

Precision measurement of a brown dwarf mass in a binary system in the microlensing event OGLE-2019-BLG-0033/MOA-2019-BLG-035

Context. Brown dwarfs are poorly understood transition objects between stars and planets, with several competing mechanisms having been proposed for their formation. Mass measurements are generally difficult for isolated objects but also for brown dwarfs orbiting low-mass stars, which are often too faint for spectroscopic follow-up. Aims. Microlensing provides an alternative tool for the discovery and investigation of such faint systems. Here we present the analysis of the microlensing event OGLE-2019-BLG-0033/MOA-2019-BLG-035, which is due to a binary system composed of a brown dwarf orbiting a red dwarf. Methods. Thanks to extensive ground observations and the availability of space observations from Spitzer, it has been possible to obtain accurate estimates of all microlensing parameters, including parallax, source radius and orbital motion of the binary lens. Results. After accurate modeling, we find that the lens is composed of a red dwarf with mass $M_1 = 0.149 \pm 0.010M_\odot$ and a brown dwarf with mass $M_2 = 0.0463 \pm 0.0031M_\odot$, at a projected separation of $a_\perp = 0.585$ au. The system has a peculiar velocity that is typical of old metal-poor populations in the thick disk. Percent precision in the mass measurement of brown dwarfs has been achieved only in a few microlensing events up to now, but will likely become common with the Roman space telescope.

preprint2019arXiv

Full orbital solution for the binary system in the northern Galactic disc microlensing event Gaia16aye

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

preprint2019arXiv

OGLE-2013-BLG-0911Lb: A Secondary on the Brown-Dwarf Planet Boundary around an M-dwarf

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q~0.03 which differs from that reported in Shvartzvald+2016. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s~0.15 or s~7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M-dwarf orbited by a massive Jupiter companion at very close (M_{host}=0.30^{+0.08}_{-0.06} M_{Sun}, M_{comp}=10.1^{+2.9}_{-2.2} M_{Jup}, a_{exp}=0.40^{+0.05}_{-0.04} au) or wide (M_{host}=0.28^{+0.10}_{-0.08} M_{Sun}, M_{comp}=9.9^{+3.8}_{-3.5}M_{Jup}, a_{exp}=18.0^{+3.2}_{-3.2} au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q=0.03 which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M_{Jup}. It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M_{Jup} in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M-dwarfs and even BDs which is challenging for other exoplanet detection methods.

preprint2019arXiv

Rest-frame UV Properties of Luminous Strong Gravitationally Lensed Ly$α$ Emitters from the BELLS GALLERY Survey

We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lya emitters (LAEs) at $2.36<z<2.82$ selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as L(Lya) = 30x10$^{42}$ erg s-1 and M(UV) = -23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lya rest-frame equivalent widths, EW(Lya), ranging from 16Åto 50Åand Lya escape fractions, fesc(Lya), from 10% to 40%. Large EW(Lya) and fesc(Lya) are found predominantly in LAEs showing weak low-ionization ISM absorption (EW < 1Å) and narrow Lya profiles (< 300 km s-1 FWHM) with their peak close (< 80 km s-1) to their systemic redshifts, suggestive of less scatter from low HI column densities that favours the escape of Lya photons. We infer stellar metallicities of Z/Zsun ~ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines NV1240Åand CIV1549Åwith those from stellar synthesis models. We also find a trend between M(UV) and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEs show star formation rates up to 180 Msun yr-1, yet they appear relatively blue ($β$(UV) ~ -1.8 to -2.0) showing evidence of little dust attenuation (E(B-V) = 0.10-0.14). These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

preprint2016arXiv

Campaign 9 of the $K2$ Mission: Observational Parameters, Scientific Drivers, and Community Involvement for a Simultaneous Space- and Ground-based Microlensing Survey

$K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $π_{\rm E}$ for $\gtrsim$127 microlensing events. These will include several that are planetary in nature as well as many short-timescale microlensing events, which are potentially indicative of free-floating planets (FFPs). These satellite parallax measurements will in turn allow for the direct measurement of the masses of and distances to the lensing systems. In this white paper we provide an overview of the $K2$C9 space- and ground-based microlensing survey. Specifically, we detail the demographic questions that can be addressed by this program, including the frequency of FFPs and the Galactic distribution of exoplanets, the observational parameters of $K2$C9, and the array of resources dedicated to concurrent observations. Finally, we outline the avenues through which the larger community can become involved, and generally encourage participation in $K2$C9, which constitutes an important pathfinding mission and community exercise in anticipation of $WFIRST$.

preprint2016arXiv

First simultaneous microlensing observations by two space telescopes: $Spitzer$ & $Swift$ reveal a brown dwarf in event OGLE-2015-BLG-1319

Simultaneous observations of microlensing events from multiple locations allow for the breaking of degeneracies between the physical properties of the lensing system, specifically by exploring different regions of the lens plane and by directly measuring the "microlens parallax". We report the discovery of a 30-55$M_J$ brown dwarf orbiting a K dwarf in microlensing event OGLE-2015-BLG-1319. The system is located at a distance of $\sim$5 kpc toward the Galactic bulge. The event was observed by several ground-based groups as well as by $Spitzer$ and $Swift$, allowing the measurement of the physical properties. However, the event is still subject to an 8-fold degeneracy, in particular the well-known close-wide degeneracy, and thus the projected separation between the two lens components is either $\sim$0.25 AU or $\sim$45 AU. This is the first microlensing event observed by $Swift$, with the UVOT camera. We study the region of microlensing parameter space to which $Swift$ is sensitive, finding that while for this event $Swift$ could not measure the microlens parallax with respect to ground-based observations, it can be important for other events. Specifically, for detecting nearby brown dwarfs and free-floating planets in high magnification events.

preprint2016arXiv

Mass Measurements of Isolated Objects from Space-based Microlensing

We report on the mass and distance measurements of two single-lens events from the 2015 \emph{Spitzer} microlensing campaign. With both finite-source effect and microlens parallax measurements, we find that the lens of OGLE-2015-BLG-1268 is very likely a brown dwarf. Assuming that the source star lies behind the same amount of dust as the Bulge red clump, we find the lens is a $45\pm7$ $M_{\rm J}$ brown dwarf at $5.9\pm1.0$ kpc. The lens of of the second event, OGLE-2015-BLG-0763, is a $0.50\pm0.04$ $M_\odot$ star at $6.9\pm1.0$ kpc. We show that the probability to definitively measure the mass of isolated microlenses is dramatically increased once simultaneous ground- and space-based observations are conducted.

preprint2016arXiv

OGLE-2015-BLG-0479LA,B: Binary Gravitational Microlens Characterized by Simultaneous Ground-based and Space-based Observation

We present a combined analysis of the observations of the gravitational microlensing event OGLE-2015-BLG-0479 taken both from the ground and by the {\it Spitzer Space Telescope}. The light curves seen from the ground and from space exhibit a time offset of $\sim 13$ days between the caustic spikes, indicating that the relative lens-source positions seen from the two places are displaced by parallax effects. From modeling the light curves, we measure the space-based microlens parallax. Combined with the angular Einstein radius measured by analyzing the caustic crossings, we determine the mass and distance of the lens. We find that the lens is a binary composed of two G-type stars with masses $\sim 1.0\ M_\odot$ and $\sim 0.9\ M_\odot$ located at a distance $\sim 3$ kpc. In addition, we are able to constrain the complete orbital parameters of the lens thanks to the precise measurement of the microlens parallax derived from the joint analysis. In contrast to the binary event OGLE-2014-BLG-1050, which was also observed by {\it Spitzer}, we find that the interpretation of OGLE-2015-BLG-0479 does not suffer from the degeneracy between $(\pm,\pm)$ and $(\pm,\mp)$ solutions, confirming that the four-fold parallax degeneracy in single-lens events collapses into the two-fold degeneracy for the general case of binary-lens events. The location of the blend in the color-magnitude diagram is consistent with the lens properties, suggesting that the blend is the lens itself. The blend is bright enough for spectroscopy and thus this possibility can be checked from future follow-up observations.

preprint2016arXiv

Statistical analysis of $m/n$ = 2/1 locked and quasi-stationary modes with rotating precursors at DIII-D

A database has been developed to study the evolution, the nonlinear effects on equilibria, and the disruptivity of locked and quasi-stationary modes with poloidal and toroidal mode numbers $m=2$ and $n=1$ at DIII-D. The analysis of 22,500 discharges shows that more than 18% of disruptions are due to locked or quasi-stationary modes with rotating precursors (not including born locked modes). A parameter formulated by the plasma internal inductance $l_i$ divided by the safety factor at 95% of the poloidal flux, $q_{95}$, is found to exhibit predictive capability over whether a locked mode will cause a disruption or not, and does so up to hundreds of milliseconds before the disruption. Within 20 ms of the disruption, the shortest distance between the island separatrix and the unperturbed last closed flux surface, referred to as $d_{edge}$, performs comparably to $l_i/q_{95}$ in its ability to discriminate disruptive locked modes. Out of all parameters considered, $d_{edge}$ also correlates best with the duration of the locked mode. Within 50 ms of a locked mode disruption, average behavior includes exponential growth of the $n=1$ perturbed field, which might be due to the 2/1 locked mode. Even assuming the aforementioned 2/1 growth, disruptivity following a locked mode shows little dependence on island width up to 20 ms before the disruption. Separately, greater deceleration of the rotating precursor is observed when the wall torque is large. At locking, modes are often observed to align with a residual error field. Timescales associated with the mode evolution are also studied and dictate the response times necessary for disruption avoidance and mitigation.

preprint2015arXiv

Spitzer Microlens Measurement of a Massive Remnant in a Well-Separated Binary

We report the detection and mass measurement of a binary lens OGLE-2015-BLG-1285La,b, with the more massive component having $M_1>1.35\,M_\odot$ (80% probability). A main-sequence star in this mass range is ruled out by limits on blue light, meaning that a primary in this mass range must be a neutron star or black hole. The system has a projected separation $r_\perp= 6.1\pm 0.4\,{\rm AU}$ and lies in the Galactic bulge. These measurements are based on the "microlens parallax" effect, i.e., comparing the microlensing light curve as seen from $Spitzer$, which lay at $1.25\,{\rm AU}$ projected from Earth, to the light curves from four ground-based surveys, three in the optical and one in the near infrared. Future adaptive optics imaging of the companion by 30m class telescopes will yield a much more accurate measurement of the primary mass. This discovery both opens the path and defines the challenges to detecting and characterizing black holes and neutron stars in wide binaries, with either dark or luminous companions. In particular, we discuss lessons that can be applied to future $Spitzer$ and $Kepler$ K2 microlensing parallax observations.

preprint2015arXiv

Spitzer Parallax of OGLE-2015-BLG-0966: A Cold Neptune in the Galactic Disk

We report the detection of a Cold Neptune m_planet=21+/-2MEarth orbiting a 0.38MSol M dwarf lying 2.5-3.3 kpc toward the Galactic center as part of a campaign combining ground-based and Spitzer observations to measure the Galactic distribution of planets. This is the first time that the complex real-time protocols described by Yee et al. (2015), which aim to maximize planet sensitivity while maintaining sample integrity, have been carried out in practice. Multiple survey and follow-up teams successfully combined their efforts within the framework of these protocols to detect this planet. This is the second planet in the Spitzer Galactic distribution sample. Both are in the near-to-mid disk and clearly not in the Galactic bulge.

preprint2014arXiv

ExELS: an exoplanet legacy science proposal for the ESA Euclid mission. II. Hot exoplanets and sub-stellar systems

The Exoplanet Euclid Legacy Survey (ExELS) proposes to determine the frequency of cold exoplanets down to Earth mass from host separations of ~1 AU out to the free-floating regime by detecting microlensing events in Galactic Bulge. We show that ExELS can also detect large numbers of hot, transiting exoplanets in the same population. The combined microlensing+transit survey would allow the first self-consistent estimate of the relative frequencies of hot and cold sub-stellar companions, reducing biases in comparing "near-field" radial velocity and transiting exoplanets with "far-field" microlensing exoplanets. The age of the Bulge and its spread in metallicity further allows ExELS to better constrain both the variation of companion frequency with metallicity and statistically explore the strength of star-planet tides. We conservatively estimate that ExELS will detect ~4100 sub-stellar objects, with sensitivity typically reaching down to Neptune-mass planets. Of these, ~600 will be detectable in both Euclid's VIS (optical) channel and NISP H-band imager, with ~90% of detections being hot Jupiters. Likely scenarios predict a range of 2900-7000 for VIS and 400-1600 for H-band. Twice as many can be expected in VIS if the cadence can be increased to match the 20-minute H-band cadence. The separation of planets from brown dwarfs via Doppler boosting or ellipsoidal variability will be possible in a handful of cases. Radial velocity confirmation should be possible in some cases, using 30-metre-class telescopes. We expect secondary eclipses, and reflection and emission from planets to be detectable in up to ~100 systems in both VIS and NISP-H. Transits of ~500 planetary-radius companions will be characterised with two-colour photometry and ~40 with four-colour photometry (VIS,YJH), and the albedo of (and emission from) a large sample of hot Jupiters in the H-band can be explored statistically.

preprint2013arXiv

ExELS: an exoplanet legacy science proposal for the ESA Euclid mission I. Cold exoplanets

Euclid is the second M-class mission of the ESA Cosmic Vision programme, with the principal science goal of studying dark energy. Euclid is also expected to undertake additional Legacy Science programmes. One proposal is the Exoplanet Euclid Legacy Survey (ExELS) which will be the first survey able to measure the abundance of exoplanets down to Earth mass for host separations from ~1AU out to the free-floating (unbound) regime. The cold and free-floating exoplanet regimes represent a crucial discovery space for testing planet formation theories. ExELS will use the gravitational microlensing technique and will detect 1000 microlensing events per month over 1.6 deg^2 of the Galactic bulge. We assess how many of these events will have detectable planetary signatures using a detailed multi-wavelength microlensing simulator (MABuLS) which incorporates the Besancon Galactic model with 3D extinction. MABuLS is the first theoretical simulation of microlensing to treat the effects of point spread function (PSF) blending self-consistently with the underlying Galactic model. We use MABuLS, together with current numerical models for the Euclid PSFs, to explore a number of designs and de-scope options for ExELS, including the exoplanet yield as a function of filter choice and slewing time, and the effect of systematic photometry errors. Using conservative extrapolations of current empirical exoplanet mass functions determined from ground-based microlensing and radial velocity surveys, ExELS can expect to detect a few hundred cold exoplanets around mainly G, K and M-type stellar hosts, including ~45 Earth-mass planets and ~6 Mars-mass planets for an observing programme totalling 10 months. ExELS will be capable of measuring the cold exoplanet mass function down to Earth mass or below, with orbital separations from ~1AU to the free-floating regime. (Abridged)

preprint2013arXiv

Photometric properties and luminosity function of nearby massive early-type galaxies

We perform photometric analyses for a bright early-type galaxy (ETG) sample with 2949 galaxies ($M_{\rm r}<-22.5$ mag) in the redshift range of 0.05 to 0.15, drawn from the SDSS DR7 with morphological classification from Galaxy Zoo 1. We measure the Petrosian and isophotal magnitudes, as well as the corresponding half-light radius for each galaxy. We find that for brightest galaxies ($M_{\rm r}<-23$ mag), our Petrosian magnitudes, and isophotal magnitudes to 25 ${\rm mag/arcsec^2}$ and 1\% of the sky brightness are on average 0.16 mag, 0.20 mag, and 0.26 mag brighter than the SDSS Petrosian values, respectively. In the first case the underestimations are caused by overestimations in the sky background by the SDSS PHOTO algorithm, while the latter two are also due to deeper photometry. Similarly, the typical half-light radii ($r_{50}$) measured by the SDSS algorithm are smaller than our measurements. As a result, the bright-end of the $r$-band luminosity function is found to decline more slowly than previous works. Our measured luminosity densities at the bright end are more than one order of magnitude higher than those of Blanton et al. (2003), and the stellar mass densities at $M_{\ast}\sim 5\times10^{11} M_{\odot}$ and $M_{\ast}\sim 10^{12} M_{\odot}$ are a few tenths and a factor of few higher than those of Bernardi et al. (2010). These results may significantly alleviate the tension in the assembly of massive galaxies between observations and predictions of the hierarchical structure formation model.

preprint2012arXiv

Influence of baryons on the orbital structure of dark matter haloes

We explore the dynamical signatures imprinted by baryons on dark matter haloes during the formation process using the OverWhelmingly Large Simulations (OWLS), a set of state-of-the-art high-resolution cosmological hydrodynamical simulations. We present a detailed study of the effects of the implemented feedback prescriptions on the orbits of dark matter particles, stellar particles and subhaloes, analysing runs with no feedback, with stellar feedback and with feedback from supermassive black holes. We focus on the central regions (0.25 r_{200}) of haloes with virial masses ~ 6 x 10^{13} (~ 7 x 10^{11}) Msun/h at z = 0(2). We also investigate how the orbital content (relative fractions of the different orbital types) of these haloes depends on several key parameters such as their mass, redshift and dynamical state. The results of spectral analyses of the orbital content of these simulations are compared, and the change in fraction of box, tube and irregular orbits is quantified. Box orbits are found to dominate the orbital structure of dark matter haloes in cosmological simulations. There is a strong anticorrelation between the fraction of box orbits and the central baryon fraction. While radiative cooling acts to reduce the fraction of box orbits, strong feedback implementations result in a similar orbital distribution to that of the dark matter only case. The orbital content described by the stellar particles is found to be remarkably similar to that drawn from the orbits of dark matter particles, suggesting that either they have forgotten their dynamical history, or that subhaloes bringing in stars are not biased significantly with respect to the main distribution. The orbital content of the subhaloes is in broad agreement with that seen in the outer regions of the particle distributions.

preprint2012arXiv

Moderate Galaxy-Galaxy Lensing

We study moderate gravitational lensing where a background galaxy is magnified substantially, but not multiply imaged, by an intervening galaxy. We focus on the case where both the lens and source are elliptical galaxies. The signatures of moderate lensing include isophotal distortions and systematic shifts in the fundamental plane and Kormendy relation, which can potentially be used to statistically determine the galaxy mass profiles. These effects are illustrated using Monte Carlo simulations of galaxy pairs where the foreground galaxy is modelled as a singular isothermal sphere model and observational parameters appropriate for the Large Synoptic Survey Telescope (LSST). The range in radius probed by moderate lensing will be larger than that by strong lensing, and is in the interesting regime where the density slope may be changing.

preprint2012arXiv

Molecular Gas in Infrared Ultraluminous QSO Hosts

We report CO detections in 17 out of 19 infrared ultraluminous QSO (IR QSO) hosts observed with the IRAM 30m telescope. The cold molecular gas reservoir in these objects is in a range of 0.2--2.1$\times 10^{10}M_\odot$ (adopting a CO-to-${\rm H_2}$ conversion factor $α_{\rm CO}=0.8 M_\odot {\rm (K km s^{-1} pc^2)^{-1}}$). We find that the molecular gas properties of IR QSOs, such as the molecular gas mass, star formation efficiency ($L_{\rm FIR}/L^\prime_{\rm CO}$) and the CO (1-0) line widths, are indistinguishable from those of local ultraluminous infrared galaxies (ULIRGs). A comparison of low- and high-redshift CO detected QSOs reveals a tight correlation between L$_{\rm FIR}$ and $L^\prime_{\rm CO(1-0)}$ for all QSOs. This suggests that, similar to ULIRGs, the far-infrared emissions of all QSOs are mainly from dust heated by star formation rather than by active galactic nuclei (AGNs), confirming similar findings from mid-infrared spectroscopic observations by {\it Spitzer}. A correlation between the AGN-associated bolometric luminosities and the CO line luminosities suggests that star formation and AGNs draw from the same reservoir of gas and there is a link between star formation on $\sim$ kpc scale and the central black hole accretion process on much smaller scales.

preprint2012arXiv

OGLE-2008-BLG-510: first automated real-time detection of a weak microlensing anomaly - brown dwarf or stellar binary?

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

preprint2010arXiv

EUCLID : Dark Universe Probe and Microlensing planet Hunter

There is a remarkable synergy between requirements for Dark Energy probes by cosmic shear measurements and planet hunting by microlensing. Employing weak and strong gravitational lensing to trace and detect the distribution of matter on cosmic and Galactic scales, but as well as to the very small scales of exoplanets is a unique meeting point from cosmology to exoplanets. It will use gravity as the tool to explore the full range of masses not accessible by any other means. EUCLID is a 1.2m telescope with optical and IR wide field imagers and slitless spectroscopy, proposed to ESA Cosmic Vision to probe for Dark Energy, Baryonic acoustic oscillation, galaxy evolution, and an exoplanet hunt via microlensing. A 3 months microlensing program will already efficiently probe for planets down to the mass of Mars at the snow line, for free floating terrestrial or gaseous planets and habitable super Earth. A 12+ months survey would give a census on habitable Earth planets around solar like stars. This is the perfect complement to the statistics that will be provided by the KEPLER satellite, and these missions combined will provide a full census of extrasolar planets from hot, warm, habitable, frozen to free floating.

preprint2010arXiv

Made-to measure galaxy models - I Methodology

We re-derive the made-to-measure method of Syer and Tremaine 1996, for modelling stellar systems and individual galaxies, and demonstrate how extensions to the made-to-measure method may be implemented and used. We illustrate the enhanced made-to-measure method by determining the mass-to-light ratio of a galaxy modelled as a Plummer sphere. From the standard galactic observables of surface brightness and line-of-sight velocity dispersion together with the h_4 Gauss-Hermite coefficient of the line-of-sight velocity distribution, we successfully recover the true mass-to-light ratio of our toy galaxy. Using kinematic data from Kleyna et al 2002, we then estimate the mass-to-light ratio of the dwarf spheroidal galaxy Draco achieving a V-band value of 539 \pm 136 M_{\odot} / L_{\odot}. We describe the main aspects of creating a made-to-measure galaxy model and show how the key modelling parameters may be determined.

preprint2010arXiv

Satellites in the field and lens galaxies: SDSS/COSMOS vs. SLACS/CLASS

The incidence of sub-galactic level substructures is an important quantity, as it is a generic prediction of high-resolution Cold Dark Matter (CDM) models which is susceptible to observational test. Confrontation of theory with observations is currently in an uncertain state. In particular, gravitational lens systems appear to show evidence for flux ratio anomalies, which are expected from CDM substructures although not necessarily in the same range of radius as observed. However, the current small samples of lenses suggest that the lens galaxies in these systems are unusually often accompanied by luminous galaxies. Here we investigate a large sample of unlensed elliptical galaxies from the COSMOS survey, and determine the fraction of objects with satellites, in excess of background counts, as a function of satellite brightness and separation from the primary object. We find that the incidence of luminous satellites within 20 kpc is typically a few tenths of one percent for satellites of a few tenths of the primary flux, comparable to what is observed for the wider but shallower SDSS survey. Although the environments of lenses in the SLACS survey are compatible with this observation, the CLASS radio survey lenses are significantly in excess of this.

preprint2010arXiv

Substructure Lensing: Effects of Galaxies, Globular Clusters & Satellite Streams

Lensing flux-ratio anomalies have been frequently observed and taken as evidence for the presence of abundant dark matter substructures in lensing galaxies, as predicted by the cold dark matter (CDM) model of cosmogony. In previous work, we examined the cusp-caustic relations of the multiple images of background quasars lensed by galaxy-scale dark matter haloes, using a suite of high-resolution N-body simulations (the Aquarius simulations). In this work, we extend our previous calculations to incorporate both the baryonic and diffuse dark components in lensing haloes. We include in each lensing simulation: (1) a satellite galaxy population derived from a semi-analytic model applied to the Aquarius haloes, (2) an empirical Milky-Way globular cluster population and (3) satellite streams (diffuse dark component) identified in the simulations. Accounting for these extra components, we confirm our earlier conclusion that the abundance of intrinsic substructures (dark or bright, bound or diffuse) is not sufficient to explain the observed frequency of cusp-caustic violations in the CLASS survey. We conclude that the observed effect could be the result of the small number statistics of CLASS, or intergalactic haloes along the line of sight acting as additional sources of lensing flux anomalies. Another possibility is that this discrepancy signals a failure of the CDM model.

preprint2009arXiv

Effects of Dark Matter Substructures on Gravitational Lensing: Results from the Aquarius Simulations

We use high-resolution Aquarius simulations of Milky Way-sized haloes in the LCDM cosmology to study the effects of dark matter substructures on gravitational lensing. Each halo is resolved with ~ 10^8 particles (at a mass resolution ~ 10^3-4 M_sun/h) within its virial radius. Subhaloes with masses larger than 10^5 M_sun/h are well resolved, an improvement of at least two orders of magnitude over previous lensing studies. We incorporate a baryonic component modelled as a Hernquist profile and account for the response of the dark matter via adiabatic contraction. We focus on the "anomalous" flux ratio problem, in particular on the violation of the cusp-caustic relation due to substructures. We find that subhaloes with masses less than ~ 10^8 M_sun/h play an important role in causing flux anomalies; such low mass subhaloes have been unresolved in previous studies. There is large scatter in the predicted flux ratios between different haloes and between different projections of the same halo. In some cases, the frequency of predicted anomalous flux ratios is comparable to that observed for the radio lenses, although in most cases it is not. The probability for the simulations to reproduce the observed violations of the cusp lenses is about 0.001. We therefore conclude that the amount of substructure in the central regions of the Aquarius haloes is insufficient to explain the observed frequency of violations of the cusp-caustic relation. These conclusions are based purely on our dark matter simulations which ignore the effect of baryons on subhalo survivability.

preprint2009arXiv

Impacts of Source Properties on Strong Lensing by Rich Galaxy Clusters

We use a high-resolution $N$-body simulation to investigate the influence of background galaxy properties, including redshift, size, shape and clustering, on the efficiency of forming giant arcs by gravitational lensing of rich galaxy clusters. Two large sets of ray-tracing simulations are carried out for 10 massive clusters at two redshifts, i.e. $z_{\rm l} \sim 0.2$ and 0.3. The virial mass ($M_{\rm vir}$) of the simulated lens clusters at $z\sim0.2$ ranges from $6.8\times10^{14} h^{-1} {M_{\odot}}$ to $1.1\times 10^{15} h^{-1} M_{\odot}$. The information of background galaxies brighter than 25 magnitude in the $I$-band is taken from Cosmological Evolution Survey (COSMOS) imaging data. Around $1.7\times 10^5$ strong lensing realizations with these images as background galaxies have been performed for each set. We find that the efficiency for forming giant arcs for $z_{\rm l}=0.2$ clusters is broadly consistent with observations. The efficiency of producing giant arcs by rich clusters is weakly dependent on the source size and clustering. Our principal finding is that a small proportion ($\sim 1/3$) of galaxies with elongated shapes (e.g. ellipticity $ε=1-b/a>0.5$) can boost the number of giant arcs substantially. Compared with recent studies where a uniform ellipticity distribution from 0 to 0.5 is used for the sources, the adoption of directly observed shape distribution increases the number of giant arcs by a factor of $\sim2$. Our results indicate that it is necessary to account for source information and survey parameters (such as point-spread-function, seeing) to make correct predictions of giant arcs and further to constrain the cosmological parameters.(abridged)

preprint2009arXiv

Properties of Long Gamma-Ray Burst Host Galaxies in Cosmological Simulation

We use galaxy catalogues constructed by combining high-resolution N-body simulations with semi-analytic models of galaxy formation to study the properties of Long Gamma-Ray Burst (LGRB) host galaxies. We assume that LGRBs originate from the death of massive young stars and analyse how results are affected by different metallicity constraints on the progenitor stars. As expected, the host sample with no metallicity restriction on the progenitor stars provides a perfect tracer of the cosmic star formation history. When LGRBs are required to be generated by low-metallicity stars, they trace a decreasing fraction of the cosmic star formation rate at lower redshift, as a consequence of the global increase in metallicity. We study the properties of host galaxies up to high redshift (~9), finding that they typically have low-metallicity (Z<0.5 Z_sun) and that they are small (M<10^9 M_sun), bluer and younger than the average galaxy population, in agreement with observational data. They are also less clustered than typical L_* galaxies in the Universe, and their descendents are massive, red and reside in groups of galaxies with halo mass between 10^{13} M_sun to 10^{14} M_sun.

preprint2006arXiv

Luminous Infrared Galaxies in the Local Universe

We study the morphology and star formation properties of 159 local luminous infrared galaxy (LIRG) using multi-color images from Data Release 2 (DR2) of the Sloan Digital Sky Survey (SDSS). The LIRGs are selected from a cross-correlation analysis between the IRAS survey and SDSS. They are all brighter than 15.9 mag in the r-band and below redshift ~ 0.1, and so can be reliably classified morphologically. We find that the fractions of interacting/merging and spiral galaxies are ~ 48% and ~ 40% respectively. Our results complement and confirm the decline (increase) in the fraction of spiral (interacting/merging) galaxies from z ~1 to z ~ 0.1, as found by Melbourne, Koo & Le Floc'h (2005). About 75% of spiral galaxies in the local LIRGs are barred, indicating that bars may play an important role in triggering star formation rates > 20 M_{sun}/yr in the local universe. Compared with high redshift LIRGs, local LIRGs have lower specific star formation rates, smaller cold gas fractions and a narrower range of stellar masses. Local LIRGs appear to be either merging galaxies forming intermediate mass ellipticals or spiral galaxies undergoing high star formation activities regulated by bars.

preprint2006arXiv

Smoothing Algorithms and High-order Singularities in Gravitational Lensing

We propose a new smoothing method for obtaining surface densities from discrete particle positions from numerical simulations. This is an essential step for many applications in gravitational lensing. This method is based on the ``scatter'' interpretation of the discrete density field in the Smoothed Particle Hydrodynamics. We use Monte Carlo simulations of uniform density fields and one isothermal ellipsoid to empirically derive the noise properties, and best smoothing parameters (such as the number of nearest neighbors used). A cluster from high-resolution simulations is then used to assess the reality of high-order singularities such as swallowtails and butterflies in caustics, which are important for the interpretation of substructures in gravitational lenses. We also compare our method with the Delaunay tesselation field estimator using the galaxy studied by Bradac et al. (2004), and find good agreements. We show that higher order singularities are not only connected with bound subhaloes but also with the satellite streams. However, the presence of high-order singularities are sensitive to not only the fluctuation amplitude of the surface density, but also the detailed form of the underlying smooth lensing potential (such as ellipticity and external shear).

preprint2005arXiv

Interferometric Visibility and Closure Phase of Microlensing Events with Finite Source Size

Interferometers from the ground and space will be able to resolve the two images in a microlensing event. This will at least partially lift the inherent degeneracy between physical parameters in microlensing events. To increase the signal-to-noise ratio, intrinsically bright events with large magnifications will be preferentially selected as targets. These events may be influenced by finite source size effects both photometrically and astrometrically. Using observed finite source size events as examples, we show that the fringe visibility can be affected by 5% - 10%, and the closure phase by a few degrees: readily detectable by ground and space interferometers. Such detections will offer unique information about the lens-source trajectory relative to the baseline of the interferometers. Combined with photometric finite source size effects, interferometry offers a way to measure the angular sizes of the source and the Einstein radius accurately. Limb-darkening changes the visibility by a small amount compared with a source with uniform surface brightness, marginally detectable with ground-based instruments. We discuss the implications of our results for the plans to make interferometric observations of future microlensing events.

preprint2002arXiv

Galaxy Formation in Preheated Intergalactic Media

We outline a scenario of galaxy formation in which the gas in galaxy-forming regions was preheated to high entropy by vigorous energy feedback associated with the formation of stars in old ellipticals and bulges and with AGN activity. Such preheating likely occurred at redshifts z ~ 2-3, and can produce the entropy excess observed today in low-mass clusters of galaxies without destroying the bulk of the Lyman alpha forest. Subsequent galaxy formation is affected by the preheating, because the gas no longer follows the dark matter on galaxy scales. The hot gas around galaxy haloes has very shallow profiles and emits only weakly in the X-ray. Cooling in a preheated halo is not inside-out, because the cooling efficiency does not change significantly with radius. Only part of the gas in a protogalaxy region can cool and be accreted into the final galaxy halo. The accreted gas is likely in diffuse clouds and so does not lose angular momentum to the dark matter. Cluster ellipticals are produced by mergers of stellar systems formed prior to the preheating, while large galaxy disks form in low-density environments where gas accretion can continue to the present time.

preprint2000arXiv

Analytic Time Delays and H_0 Estimates for Gravitational Lenses

We study gravitational lens time delays for a general family of lensing potentials, which includes the popular singular isothermal elliptical potential and singular isothermal elliptical density distribution but allows general angular structure. Using a novel approach, we show that the time delay can be cast in a very simple form, depending only on the observed image positions. Including an external shear changes the time delay proportional to the shear strength, and varying the radial profile of the potential changes the time delay approximately linearly. These analytic results can be used to obtain simple estimates of the time delay and the Hubble constant in observed gravitational lenses. The naive estimates for four of five time delay lenses show surprising agreement with each other and with local measurements of H_0; the complicated Q 0957+561 system is the only outlier. The agreement suggests that it is reasonable to use simple isothermal lens models to infer H_0, although it is still important to check this conclusion by examining detailed models and by measuring more lensing time delays.