MOA-2010-BLG-073Lb

MOA-2010-BLG-073L: An M-Dwarf with a Substellar Companion at the Planet/Brown Dwarf Boundary

R.A. Street, J.-Y. Choi, Y. Tsapras, C. Han, K. Furusawa, M. Hundertmark, A. Gould, T. Sumi, I.A. Bond, D. Wouters, R. Zellem, A. Udalski
and
(The RoboNet Collaboration)
C. Snodgrass, K. Horne, M. Dominik, P. Browne, N. Kains, D.M. Bramich, D. Bajek, I.A. Steele,
S. Ipatov, and
(The MOA Collaboration)
F. Abe, D.P. Bennett, C.S. Botzler, P. Chote, M. Freeman, A. Fukui, P. Harris, Y. Itow, C.H. Ling, K. Masuda, Y. Matsubara, N. Miyake, Y. Muraki, T. Nagayama, S. Nishimaya, K. Ohnishi, N. Rattenbury, To. Saito, D.J. Sullivan, D. Suzuki, W.L. Sweatman, P.J. Tristram, K. Wada, P.C.M. Yock
and
(The OGLE Collaboration)
M. K. Szymański, M. Kubiak, G. Pietrzyński, I. Soszyński, R. Poleski, K. Ulaczyk, Ł. Wyrzykowski and
(The FUN Collaboration)
J. Yee, S. Dong, I.-G. Shin, C.-U. Lee, J. Skowron, L. Andrade De Almeida, D.L. DePoy, B.S. Gaudi, L.-W. Hung, F. Jablonski, S. Kaspi, N. Klein, K.-H. Hwang, J.-R. Koo, D. Maoz, J.A. Muñoz, R.W. Pogge, D. Polishhook, A. Shporer, J. McCormick, G. Christie, T. Natusch, B. Allen, J.Drummond, D. Moorhouse, G. Thornley, M. Knowler, M. Bos, G. Bolt,
and
(The PLANET Collaboration)
J.-P. Beaulieu, M.D. Albrow, V. Batista, S. Brillant, J.A.R. Caldwell, A. Cassan, A. Cole, E. Corrales, Ch. Coutures, S. Dieters, D. Dominis Prester, J. Donatowicz, P. Fouqué, E. Bachelet, J. Greenhill, S.R. Kane, D. Kubas, J.-B. Marquette, R. Martin, J. Menzies, K.R. Pollard, K.C. Sahu, J. Wambsganss, A. Williams, M. Zub
and
(MiNDSTEp)
K.A. Alsubai, V. Bozza, M.J. Burgdorf, S. Calchi Novati, P. Dodds, S. Dreizler, F. Finet, T. Gerner, S. Hardis, K. Harpsøe, F. Hessman, T.C. Hinse, U.G. Jørgensen, E. Kerins, C. Liebig, L. Mancini, M. Mathiasen, M.T. Penny, S. Proft, S. Rahvar, D. Ricci, G. Scarpetta, S. Schäfer, F. Schönebeck, J. Southworth, J. Surdej
LCOGT, 6740 Cortona Drive, Suite 102, Goleta, CA 93117, USA. Department of Physics, Institute for Astrophysics, Chungbuk National University, Cheongju, 361-763, Korea School of Mathematical Sciences, Queen Mary, University of London, Mile End Road, London E1 4NS, UK Solar-Terrestrial Environment Laboratory, Nagoya University, Nagoya, 464-8601, Japan SUPA/St Andrews, Dept. of Physics and Astronomy, North Haugh, St. Andrews, Fife, KY16 9SS, UK Dept. of Astronomy, Ohio State University, McPherson Laboratory, 140 W. 18th Avenue, Columbus, Ohio 43210-1173, USA Dept. of Earth and Space Science, Graduate School of Science, Osaka University, 1-1 Machikaneyama-cho, Toyonaka, Osaka 560-0043, Japan Institute of Information and Mathematical Sciences, Massey University, Private Bag 102-904, North Shore Mail Centre, Auckland, New Zealand UPMC-CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98bis boulevard Arago, F-75014, Paris, France Lunar and Planetary Laboratory, Department of Planetary Sciences, University of Arizona, 1629 E. University Blvd., Tucson AZ 85721-0092, USA Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warszawa, Poland Max Planck Institute for Solar System Research, Max-Planck-Str. 2, 37191 Katlenburg-Lindau, Germany European Southern Observatory, Karl-Schwarzschild-Straße 2, 85748 Garching bei München, Germany Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead, Wirral., CH41 1LD, UK Qatar Foundation, P.O. Box 5825, Doha, Qatar Dept. of Physics, University of Notre Dame, Notre Dame, IN 46556, USA Dept. of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand School of Chemical and Physical Sciences, Victoria University, Wellington, New Zealand Dept. of Physics, Konan University, Nishiokamoto 8-9-1, Kobe 658-8501, Japan Dept. of Physics and Astrophysics, Faculty of Science, Nagoya University, Nagoya 464-8602, Japan Extrasolar Planet Detection Project Office, National Astronomical Observatory of Japan (NAOJ), Osawa 2-12-1, Mitaka, Tokyo 181-8588, Japan Nagano National College of Technology, Nagano, 381-8550, Japan Tokyo Metropolitan College of Industrial Technology, Tokyo, 116-8523, Japan Mt. John Observatory, P.O. Box 56, Lake Tekapo 8770, New Zealand Universidad de Concepción, Departamento de Astronomia, Casilla 160–C, Concepción, Chile Institute for Advanced Study, Einstein Drive, Princeton, New Jersey 08540, USA Korea Astronomy & Space Science Institute, Daejeon, 305-348, Korea Divisao de Astrofisica, Instituto Nacional de Pesquisas Espeaciais, Avenida dos Astronauntas, 1758 Sao José dos Campos, 12227-010 SP, Brazil Dept. of Physics and Astronomy, Texas A&M University College Station, TX 77843-4242, USA School of Physics and Astronomy and Wise Observatory, Tel-Aviv University, Tel-Aviv 69978, Israel Departamento de Astronomiá y Astrofísica, Universidad de Valencia, E-46100 Burjassot, Valencia, Spain Farm Cove Observatory, Farm Cove, Pakuranga, Auckland 2010, New Zealand Auckland Observatory, Auckland, New Zealand Vintage Lane Observatory, Blenheim, New Zealand Possum Observatory, Patutahi, New Zealand Kumeu Observatory, Kumeu, West Auckland, New Zealand Molehill Astronomical Observatory, North Shore City, New Zealand Craigie, Perth, Western Australia, Australia University of Canterbury, Dept. of Physics and Astronomy, Private Bag 4800, Christchurch 8020, New Zealand European Southern Observatory, Casilla 19001, Vitacura 19, Santiago, Chile McDonald Observatory, 16120 St Hwy Spur 78 #2, Fort Davis, Tx 79734, USA Institut d’Astrophysique de Paris, UMR7095 CNRS-Université Pierre & Marie Curie, 98 bis boulevard Arago, 75014 Paris, France School of Math and Physics, University of Tasmania, Private Bag 37, GPO Hobart, Tasmania 7001, Australia Physics Department, Faculty of Arts and Sciences, University of Rijeka, Omladinska 14, 51000 Rijeka, Croatia Technical University of Vienna, Department of Computing, Wiedner Hauptstrasse 10, Vienna, Austria Université de Toulouse, UPS-OMP, IRAP, Toulouse, France CNRS, IRAP, 14 Avenue Edouard Belin, F-31400 Toulouse, France NASA Exoplanet Science Institute, Caltech, MS 100-22, 770 S. Wilson Ave., Pasadena, CA 91125, USA Perth Observatory, Walnut Road, Bickley, Perth 6076, Australia South African Astronomical Observatory, P.O. Box 9, Observatory 7935, South Africa Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg (ZAH), Mönchhofstr. 12-14, 69120 Heidelberg, Germany Universitá degli Studi di Salerno, Dipartimento di Fisica “E.R. Caianiello”, Via S. Allende, 84081, Baronissi (SA), Italy SOFIA Science Center, NASA Ames Research Center, Mail Stop N211-3, Moffett Field, CA 94035, USA Istituto Internazionale per gli Alti Studi Scientifici (IIASS), Vietri Sul Mare (SA) Italy Institut für Astrophysik, Georg-August-Universität, Friedrich-Hund-Platz1, 37077 Göttingen, Germany Institut d’Astrophysique et de Géophysique, Allée du 6 Août 17, Sart Tilman, Bât. B5c, 4000 Liége, Belgium Niels Bohr Institute, University of Copenhagen, Juliane Maries vej 30, 2100 Copenhagen, Denmark Centre for Star and Planet Formation, Geological Museum, Øster Voldgade 5, 1350 Copenhagen, Denmark Jodrell Bank Centre for Astrophysics, University of Manchester, Oxford Road, Manchester, M13 9PL, UK Max Planck Institute for Astronomy, Königstuhl 17, 619117 Heidelberg, Germany Dept. of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran INFN, Gruppo Collegato di Salerno, Sezione di Napoli, Italy Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK Corresponding author Royal Society University Research Fellow rstreet@lcogt.net
Abstract

We present an analysis of the anomalous microlensing event, MOA-2010-BLG-073, announced by the Microlensing Observations in Astrophysics survey on 2010-03-18. This event was remarkable because the source was previously known to be photometrically variable. Analyzing the pre-event source lightcurve, we demonstrate that it is an irregular variable over time scales 200 d. Its dereddened color, , is 1.2210.051 mag and from our lens model we derive a source radius of 14.71.3 , suggesting that it is a red giant star. We initially explored a number of purely microlensing models for the event but found a residual gradient in the data taken prior to and after the event. This is likely to be due to the variability of the source rather than part of the lensing event, so we incorporated a slope parameter in our model in order to derive the true parameters of the lensing system. We find that the lensing system has a mass ratio of =0.06540.0006. The Einstein crossing time of the event, =44.30.1 d, was sufficiently long that the lightcurve exhibited parallax effects. In addition, the source trajectory relative to the large caustic structure allowed the orbital motion of the lens system to be detected. Combining the parallax with the Einstein radius, we were able to derive the distance to the lens, =2.80.4 kpc, and the masses of the lensing objects. The primary of the lens is an M-dwarf with =0.160.03  while the companion has =11.02.0 , putting it in the boundary zone between planets and brown dwarfs.

Gravitational lensing: micro, microlensing, MOA-2010-BLG-073, MOA-2010-BLG-073L, brown dwarf, low mass binary

1 Introduction

The mass function of individual compact objects (brown dwarfs and planets) in the Galaxy remains poorly understood, particularly at the low-mass end. Brown dwarfs (BD) are commonly defined as objects with masses between the deuterium and hydrogen burning limits (DBL and HBL, respectively) but these can be hard to detect, being intrinsically faint and fading further as they cool over time. Below the DBL, the mass function for individual objects is even more poorly measured. Unbound, free-floating objects of planetary mass have been discovered via direct imaging of clusters (for example in  Orionis (Béjar et al., 2012)) and in the field (e.g. the 6–25  object reported by Kirkpatrick et al. (2006)). Sumi et al. (2011) reported a population of planets which are either unbound or at very wide separations, discovered when their gravity caused short timescale microlensing events.

At least partially as a result of these poor constraints the origin of low-mass compact objects remain unclear. Although traditionally thought of as separate classes of objects, planets and brown dwarfs form a continuous scale of mass and are best distinguished by the circumstances of their formation (Burrows et al., 2001; Chabrier et al., 2005; Sahlmann et al., 2010; Chabrier et al., 2011). Planets form in disks of material orbiting a protostellar object and may subsequently migrate to different period orbits. Brown dwarfs on the other hand, are considered to be the extreme low mass end of the star formation process by fragmentation of locally over dense cores caused by turbulence in a cloud (Chabrier et al., 2011), which can themselves form protoplanetary disks (Klein et al., 2003; Scholz et al., 2006). This mechanism may also produce objects of a few Jupiter masses. For a recent review, see Luhman (2012).

The mass function of free-floating low-mass objects is likely to be different from those bound to stars. Marcy & Butler (2000) identified a paucity of BDs orbiting close (3 AU) to their host stars, a region where Jovian-mass planets are commonly found. This “brown dwarf desert” may represent the gap between the largest objects that can form in protoplanetary disks and the smallest objects that can concurrently collapse/condense next to a star.

Two different theories have been proposed to explain the formation of Jovian planets in disks (for a review, see Zhou et al. (2012)). The core accretion model predicts planets form from protoplanetary cores, growing up to tens of Jupiter masses (Mordasini et al., 2009; Baraffe et al., 2010, e.g.) but predicts few giant planets around M-dwarfs. Higher mass stars are thought to have disks with enhanced surface densities which allow the cores to grow more rapidly (Laughlin et al., 2004), as do disks with a high fraction of dust, leading to enhanced planet formation around high metallicity stars (Ida & Lin, 2004). Alternatively, the model of planet formation via gravitational instabilities in the disk (Boss, 2006, e.g.) tends to favor the formation of more massive planets, in generally wider orbits.

A number of lines of evidence support the core accretion theory. There is a well-established correlation of increasing planet frequency with stellar metallicity (Santos, Israelian & Mayor, 2001; Fischer & Valenti, 2005; Maldonado et al., 2012). The results of radial velocity surveys imply there is a derth of M-dwarf stars with massive, close-in planets. In part, this reflects an observational bias against these faint objects but the sample is sufficiently large that a real statistical trend is emerging (Cumming et al., 2008; Johnson et al., 2010; Bonfils et al., 2011), for companions with  d. Recent spectroscopic and Kepler results have confirmed the prediction of a rapid increase in frequency for planets with small radii (down to 2) and  d for all spectral types, and found that these small planets are several times more common around stars of late spectral type (Bonfils et al., 2011; Howard et al, 2012).

However, the core accretion model has difficulty forming massive planets at large orbital radii and such systems have been discovered, for example HR 8799 (Marois et al., 2008). Furthermore, a number of planets have been found orbiting M-dwarf hosts at larger orbital separations, for example Dong et al. (2006); Forveille et al. (2011); Batista et al. (2011). These systems may instead form through gravitational instability in the disk, which can account for companions up to several Jupiter masses around M-dwarfs (Boss, 2006).

To better understand the formation mechanisms of heavy substellar companions in bound systems we need to trace the distributions of the physical and orbital properties (such as mass ratio, orbital separation, occurrence frequency) of a significant number of systems. Yet relatively few bound brown dwarf companions have been reported, despite their being easy to detect at close orbital separations (the “brown dwarf desert”).

Microlensing offers a complementary window onto BD and planet formation by probing for cooler companions of all masses in orbital radii between 0.2 – 10 AU, separations which are difficult or time consuming to explore by other methods (Shin et al., 2012a). It can probe the companion mass function down to M- and brown dwarf hosts, and is sensitive to companions from nearly equal mass down to terrestrial masses.

Sixteen systems have been published to date111Listed on exoplanet.eu, and thanks to large-scale galactic lensing surveys and efficient follow-up, each season’s Bulge observing campaign is now producing a regular yield of new discoveries (e.g. Bachelet et al. (2012); Yee et al. (2012); Miyake et al. (2012)). Of these 16 companions, 3 are giant planets orbiting M-dwarf stars: OGLE-2005-BLG-071Lb, a 3.8  planet (Dong et al., 2009), MOA-2009-BLG-0387Lb, with =2.5  planet (Batista et al., 2011) and MOA-2011-BLG-293Lb, which hosts a 2.4  companion (Yee et al., 2012).

Here we present the newly discovered system, MOA-2010-BLG-073L (catalog MOA-2010-BLG-073L), an M-dwarf star with a companion whose mass is close to the deuterium burning limit of 12.6 . The discovery and follow-up observations are described in § 2 and we discuss the variability of the source star in § 3. We describe our analysis in § 4, from which we derive the physical properties of the lens in § 5. Finally, we discuss our findings in § 6.

2 Observations

The microlensing event, MOA-2010-BLG-073 (catalog MOA-2010-BLG-073), was first announced by the Microlensing Observations in Astrophysics222www.phys.canterbury.ac.nz/moa (MOA, Bond et al. (2001); Sumi et al. (2003) on the 1.8 m telescope at Mt. John Observatory, New Zealand) survey on 2010-03-18. A background source star in the Galactic Bulge, =18:10:11.342, =-26:31:22.544 (J2000.0), previously having a mean baseline magnitude of 16.5 mag, was discovered to be rising smoothly in brightness consistent with a point-source, point-lens (PSPL) microlensing event. However, on 2010-05-03 the event was found to show an anomalous brightening of 0.5 mag and an alert was issued (K. Furusawa, private comm.).

Microlensing follow-up teams worldwide – RoboNet-II333robonet.lcogt.net (Tsapras et al., 2009), FUN444www.astronomy.ohio-state.edu/microfun (Gould et al., 2006), PLANET555planet.iap.fr (Beaulieu et al., 2006) and MiNDSTEp666www.mindstep-science.org (Dominik et al., 2010) – responded to provide intensive coverage of the event for the duration of the anomaly (2 days), and monitored the event as it returned to baseline, over the course of the next 2 months.

In addition to the MOA data, taken with the wide band “MOA-Red” filter (corresponding to + bandpasses), the event was observed from several other sites in New Zealand. The 0.41 m telescope at Auckland Observatory, the 0.35 m at Kumeu Observatory and the 0.41 m Possum Observatory all used -band filters while the 0.304 m Molehill Astronomical Observatory (MAO), the 0.35 m telescope at Farm Cove Observatory (FCO) and the 0.4 m telescope at Vintage Lane Observatory (VLO) all observed it unfiltered. The event was then picked up from three sites in Australia, firstly in -band by the 1 m Canopus telescope in Tasmania followed by the 2 m Faulkes Telescope South (FTS), where an SDSS- filter was used. The 0.6 m telescope in Perth also observed in band. Of the observing sites around longitude zero, the event was imaged from the 1 m telescope at the South African Astronomical Observatory (SAAO) using an I band filter and in , and by the 1.4 m Infra Red Survey Facility (IRSF), also at SAAO. The 2 m Liverpool Telescope (LT) observed the event in SDSS- from the Canary Islands.

As darkness fell in the Americas, a number of Chilean telescopes picked up the observing baton: the SMARTS 1.3 m at the Cerro Tololo Interamerican Observatory (CTIO) obtained data in , and bands with the andicam camera, and the Danish 1.54 m used an filter. Though in the midst of commissioning the new OGLE-IV camera at the time, the 1.3 m Warsaw telescope also covered the event in -band777ogle.astrouw.edu.pl. The 1 m Mt. Lemmon Telescope in Arizona imaged the event in the -band and in the extreme west, the 2 m Faulkes Telescope North (FTN), Hawai’i used an SDSS-i filter to complete the 24-hour coverage of the event from the Pacific. Table 1 summarizes the data obtained, which are plotted in Figure 1.

The high density of Galactic Bulge star fields and the consequent degree of overlap (or blending) in stellar point-spread functions (PSF) has long since made difference image analysis (DIA) the photometry method of choice among microlensing teams. Both MOA and OGLE make their photometry available to the community, automatically reducing their data with their custom pipelines described respectively in Bond et al. (2001) and Udalski et al. (2003). The RoboNet data (from FTN, FTS and the LT) were reduced with the project’s automated data reduction pipeline, which is based around the DanDIA package (Bramich, 2008). This software was also later used to reduce data from Canopus, the Perth 0.6 m, the SAAO 1 m and the -band data from CTIO, while the DIAPL package was used to process the images from the Danish telescope. The PLANET team released their photometry (produced by the WISIS pipeline) in real time via their website, and the Pysis DIA pipeline (Albrow et al., 2009) was used for later re-reduction of these data sets.

Telescope & Filter N frames N frames
aperture [m] total used [mag]
MOA 1.8 / 0.7027 0.6118 1747 1726 1.305 0.005
OGLE 1.3 0.6098 0.5103 47 42 2.188 0.005
Auckland 0.41 0.7027 0.6118 136 136 0.910 0.000
Canopus 1.0 0.6098 0.5103 162 159 1.310 0.005
CTIO 1.3 0.7817 0.7048 19 18 0.603 0.000
CTIO 1.3 0.6098 0.5103 162 162 1.010 0.003
CTIO 1.3 0.4145 0.3206 586 575 1.340 0.016
Danish 1.54 0.6098 0.5103 498 491 1.130 0.014
Farm Cove 0.4 Unfiltered - 0.5611 225 225 0.975 0.000
FTN 2.0 SDSS- 0.6098 0.5103 159 158 1.055 0.011
FTS 2.0 SDSS- 0.6098 0.5103 129 129 1.125 0.006
IRSF 1.4 0.4836 0.3844 4 4 1.000 0.000
IRSF 1.4 0.4145 0.3206 4 4 1.000 0.000
IRSF 1.4 0.3550 0.2684 4 4 1.000 0.000
Kumeu 0.35 0.7027 0.6118 272 272 0.772 0.000
Lemmon 1.0 0.6098 0.5103 116 105 1.290 0.020
LT 2.0 SDSS- 0.6098 0.5103 167 167 1.155 0.006
MAO 0.304 Unfiltered - 0.5611 238 238 1.025 0.018
Perth 0.6 0.6098 0.5103 66 66 1.095 0.008
Possum 0.7027 0.6118 15 15 1.030 0.009
SAAO 1.0 0.6098 0.5103 30 30 1.050 0.011
Vintage Lane 0.4 Unfiltered - 0.5611 124 124 0.995 0.000
aafootnotetext: Includes only OGLE-IV data taken during event MOA-2010-BLG-073.bbfootnotetext: For unfiltered or very broadband observations we adopted a limb-darkening parameter value which was the average of that for R and I bands: .ccfootnotetext: , are defined by Eqns. 3, 4, 5 and , by Eqn. 6.
Table 1: Summary of observations.
Figure 1: Top panel: The complete light curve of MOA-2010-BLG-073 with the data sets from the various observatories overlaid with our best fitting model. The inset shows the anomaly in greater detail. Lower panels: The residuals of the fit for each of the best-fitting models of each class.

3 Variability of the Source Star

MOA-2010-BLG-073 was present in the fields of the OGLE-II and OGLE-III surveys so the source star’s I-band photometric record extends from 1998 to 2006 (Fig. 2). OGLE-IV was in the commissioning phase when this event took place. From this excellent baseline it was immediately clear that the source is variable over many-month time scales. This raised the possibility that shorter-term variability might obfuscate the microlensing signal, making it difficult to determine its properties.

To investigate this possibility, we performed a search for periodicities in the baseline OGLE-II and OGLE-III data, excluding the lensing event, using the ANOVA algorithm (Schwarzenberg-Czerny, 1996). Due to the seasonal gaps in the baseline, we analyzed the OGLE-II data in yearly subsets as it is the best sampled, searching for periods between =0.5–200 d. As Figure 3 demonstrates, there are no significant or persistent periodicities, other than the expected integer multiples of the 1-day sampling alias.

Figure 2: All available OGLE-II (1998–2001) and OGLE-III (June 2001–2006) lightcurve data for the source star of MOA-2010-BLG-073, taken prior to the event. The instrument upgrade to OGLE-III occurred around HJD=2452000.0. The apparent offset in magnitude at this time may be due in part to the difficulties of accurately calibrating photometry between the two wide field surveys.
Figure 3: Periodograms of the OGLE-II (1998–2001) seasons of data, analyzed separately with the ANOVA algorithm. No significant periodicities were found in the range =0.5–200 d, aside from the expected aliases.

We then combined the OGLE-II and III data sets in order to search for periods up to =4000 d. Figure 2 indicates a slight (0.04 mag) magnitude offset between the OGLE-II and III data. This can occur as a residual of OGLE’s photometric calibration between the two surveys but it might also be the result of the intrinsic stellar variation. Therefore, we performed a search for periods between 0.5 and 4000 d based on the combined data both with and without this offset (estimated visually). In both cases the periodogram was dominated by the window function; the only significant power was found in the peaks marking multiples of the 1-day alias, plus one peak at extreme low frequency corresponding to the finite length of the data set. We conclude that this star is an irregular, long-term variable, most likely as a result of pulsations.

However, there remained the possibility that the star could be irregularly variable on time scales comparable to that of the lensing event. To test this possibility, we binned the OGLE-II light curve on a range of time scales between 2–1200 d. To each bin we fitted two functions: one of constant brightness, and one with a linear slope and calculated the weighted root-mean-square (RMS) photometric scatter around each function per bin. We plot the average and maximum RMS (calculated over all bins in the lightcurve) against the width of the bins in time in Figure 4. On time scales shorter than 200 d the RMS scatter in the binned light curve is reasonably constant implying no significant short term variability. The longer term trend becomes clearly evident in the constant brightness curves for time scales longer than 400 d. We note that for bin widths between 800–1200 d, this curve has an RMS actually exceeding that of the whole light curve; this is because these bins were sufficiently wide that the first bin included the majority of the data, and the most variable sections of the light curve. As the bin widths became longer, they included more data points from the relatively stable section towards the end of the OGLE-II data set, and the RMS drops. The deviation of the OGLE lightcurve from a constant brightness exceeds 3 for timescales longer than 750 d. However, it deviates from a linear slope by 1.6 for timescales less than 1000 d, so we represent this variation as a gradient in the lightcurve over the duration of the event.

Figure 4: The OGLE-II photometry was binned using a series of bins of varying width in time. The data in each bin was fitted with two functions a) constant brightness b) with a linear slope, and the RMS around these curves was averaged over all the bins in the lightcurve. The average and maximum RMS for each binning is plotted here against the width of the bin, indicating variability over different timescales.

The most intuitive way to account for the variation of the source was to measure the gradient of the lightcurve taken at baseline immediately before and after the event. Unfortuantely only one of the available datasets covered these periods. Fitting a straight-line model (via a non-linear least-squares Marquardt-Levenberg algorithm) to the MOA 2009 and 2011 season data, we measured a slope of 0.018 mag/yr. However, the RMS scatter in the residuals of this fit were 0.016 mag, making it difficult to properly determine the slope. Additionally, in order to remove this slope from the other datasets, it would be necessary to transform the fluxes measured by each telescope on to the same scale as the MOA data. We attempted this via a linear regression approach but found that significant residuals remained. These contributed to overall higher  values when the corrected data were fit with binary lensing models. We therefore adopted the alternative method of incorporating the slope as an additional parameter in our lensing model which we fit to the original, uncorrected data and we describe this approach in the following sections.

4 Analysis

In the analysis of this event, we used the established modeling software developed by S. Dong and C. Han (Dong et al., 2006; Shin et al., 2012b).

4.1 Initial Parameters

As a starting point for our analysis, we needed approximate values for the three parameters of the standard model for a PSPL event (not yet including the slope; this is discussed in § 4.8): , the time of peak magnification occurring at the closest projected separation between the lens and source,  and the Einstein radius crossing time, . Following standard convention, all distances are quoted in units of the angular Einstein radius, , of the lens.

To estimate these parameters, we combined all available data sets into a single light curve using the following scaling to take account of the varying degrees of PSF blending from different instruments:

(1)

where is the measured flux of the target at time from data set , is the lensing magnification at that time, is the flux of the source star and represents the flux of all stars blended with the source in the data set. A regression fit was used to measure  and  for each data set, producing an aligned light curve. Although the resulting parameter estimates are somewhat different from their ‘true’ values due to the existence of the anomalous deviation on the light curve, they provided a starting point in parameter space.

We note that two additional parameters can contribute to a PSPL model. They are the lens parallax parameters  and  that account for the light curve deviation caused by the motion of the Earth in its orbit over the course of the event. The vector microlens parallax, , where is the Einstein radius projected onto the observer plane and:

(2)

where represents the direction of lens motion relative to the source as a counter-clockwise angle, north through east. However, this is generally significant only for long (months) time scale events, and was included at a later stage (see Section 4.6).

4.2 Finite Source

The sharp spike feature in the lightcurve is indicative of the source closely approaching or crossing a caustic. In these circumstances, it cannot be approximated by a point light source and must be treated as a disk of finite angular radius, , with wavelength-dependent limb-darkening. This is addressed within our software using the ray-shooting approach (Kayser et al., 1986): the path of light rays is traced from the image plane back to the source, taking into account the bending of the trajectory according to the lens equation. If a ray is found to “land” within the radius of the source, its intensity is computed taking limb-darkening into account. We derived this from the linear limb darkening law:

(3)

where is the intensity of the source at radius from the center, relative to the central intensity in the same wavelength, , scaled by the coefficient . While more accurate limb darkening models are available, they are not commonly used in microlensing analyses due to the complexity introduced by combining data from many sources (Bachelet et al. (2012) discussed this in more detail). The values of for each passband were calculated from the Kurucz ATLAS9 stellar atmosphere models presented by Kurucz (1979) using the method of Heyrovský (2007). However, within the microlensing community and software, Equation 3 more commonly follows the formalism derived by Albrow et al. (1999):

(4)

where is the total flux from the source in a given passband and is the angle between the line of sight to the observer and the normal to the stellar surface. The limb darkening coefficient, is related to by:

(5)

The values of and applied for each dataset are presented in Table 1. The lensing magnification is then computed as the ratio of the number of rays reaching the source plane relative to the number in the image plane. This approach is only required while the source is close to the caustic. At larger separations, the software employs a semi-analytic hexadecapole approximation to the finite source calculation to improve computation speeds (Pejcha & Heyrovský, 2009; Gould, 2008).

4.3 Standard Binary Model Grid Search

To model the light curve of a binary lens event, we introduced three additional parameters: , the ratio of the masses of the two bodies composing the lens where  is the more massive component, , the projected separation of those masses and , the angle of the trajectory of the lensed source star, relative to the lens’ binary axis. The frame of reference was defined to be at rest with respect to the Earth at time , which we took to be the time of caustic crossing at HJD=2455321.0, estimated from the easily identifiable feature in the light curve (following the notation of Skowron et al. 2011).

With seven variables in the model (, , , , , , ), a number of different lens/source configurations may produce similar light curves, so it was necessary to thoroughly explore a large area of parameter space in order to ensure all possible solutions are identified. We therefore constructed a grid of models, spanning set ranges in the values of the three variables upon which the overall  of the fit depended most sensitively, ,  and . Each node in this grid took fixed values of (,,) and used a Markov Chain Monte Carlo (MCMC) approach (Dong et al., 2006) to find the best fitting model by optimizing the other parameters. To improve efficiency, a magnification map is generated by rayshooting for each point in the grid from which the model light curves used to compute the  are drawn. The grid covered the following range: in steps of 0.012, in steps of 0.05 and in steps of 0.6.

Mapping out the  for each node in this grid, we found a number of local minima. Visual inspection of these models overlaid on the light curve demonstrated that some more closely followed the data than others. Our first pass analysis included substantial baseline photometry before and after the event. This was not well fit by the models due to the variability of the source and hence the  map gave a distorted view of regions in parameter space that best match the event. For this reason, we proceeded by repeating the grid search using just data taken during 2010. This produced two clear minima in  of which one model stood out as by far the best match to the data. We then conducted a refined grid search over this restricted region of parameter space, taking smaller incremental steps.

4.4 Optimized Standard Binary Model

The refined grid search produced a reasonable model, fitting the majority of the data from all telescopes. This was used as a guide to identifying likely outlying data points, for which the quality of the reduction was then checked. A handful of data points were removed at this stage. However, this model included only fixed values for ,  and . To properly determine the standard binary model for this event, our next step was to allow the seven parameters (, , , , , , ) to be optimized during the MCMC fitting process, which used the grid search results as its starting point. At this point we included the extended baseline data from MOA for seasons 2009 and 2011, as these fall within the period for which the source’s variability can be approximated with a straight line; we address this in Section 4.8.

4.5 Normalization of Photometric Errors

When fitting microlensing events, the reduced  of the fit on a per data set basis, , typically produces a range of values both less than and exceeding the expected unity value. This can occur as different groups have slightly different ways of estimating photometric errors, but can lead to over- or under-emphasis being placed on particular data sets during the modeling process.

A common technique to address this issue is to arrive at a complete model for the event and then use this model to renormalize the original photometric errors of each data set, , according to the expression:

(6)

We first conducted the sequence of models described in the following sections in order to find the best model for the event. We then set the coefficients , such that the  relative to that model equaled unity; the adopted values are given in Table 1. We then repeated our MCMC fitting process, starting with the Standard binary model and systematically adding parameters in to determine the extent of improvement in the model  in each case. We compare models in Tables 2 & 3, and we plot the residuals () in Figure 1. In the following sections, the  values given are those post-renormalization.

4.6 Parallax

Given that the event’s 44 d0.12 yr, it was necessary to include parallax in our model. Using the parameters of the Standard Binary Model as a starting point, we allowed our fitting process to optimize for ,  also. We found that this significantly reduced the  of the overall fit to 6064.773. By default, this procedure explored models with positive projected separations at closest approach of the lens and source, that is , which we define as the source’s trajectory passing the caustic at positive values of in the lens plane (see Figure 6). For the standard model case the symmetry with respect to the binary axis of the caustic means that the  solutions are identical. However once parallax was included, this was no longer the case, so we also explored  solutions (this degeneracy is further discussed in Park et al. (2004)). The parameters of the  model were taken as a starting point for the fit, except that the sign of  was reversed and the  value became . We found the best fitting  model to be slightly less favored, with =6099.487.

4.7 Lens Orbital Motion

The mass ratio and projected separation determined from this model put this event close to the boundary between close and intermediate/resonant caustic structure. In this regime, small changes in the projected separation of the lensing bodies due to their orbital motion can effectively change the shape of the caustic (see Fig. 6) while the event is underway, sometimes causing detectable deviations in the lightcurve. To explore this possibility, we included additional parameters in our model to describe the change in projected binary separation,  and the rate of change of the angle of the projected binary axis, . Again we found that the  model gave the best fit, with =5544.071, compared with 5690.003 when . This type of orbital motion is classified as “separational” in the schema put forward by Gaudi (2009); Penny et al. (2010), and is detected in this event as the source happens to cross the cusp of the caustic in the position where the caustic changes most rapidly.

4.8 Sloping Baseline

While taking these second-order effects into account significantly improved the fit to the data, the overall  remained rather high. Visual inspection of the light curve still showed a gradient, especially in the 2010 baseline before, relative to after, the event. Based on our analysis in § 3, this trend is likely to be part of the longer-term variability of the source and not associated with the microlens. In order to determine the true lens/source characteristics though, this trend had to be taken into account.

The OGLE-II and III data demonstrate that the source variation over the 150 d time scale of the event can be approximated by a straight-line, rather than a higher-order function. We therefore introduced a ‘slope’ parameter to our model, representing the linear rate of change in magnitude during the event. This further improved the , and the best-fitting model was once again the  solution with =4782.367.

We note that there exist degeneracies between the slope parameter and those for parallax and orbital motion as they can be used to fit similar residuals in the lightcurve. To test for possible degeneracies, we also fit a standard model plus the slope parameter alone and found that =5334.378. The value for the slope from this model, -0.01600.0004, was consistent with that derived from the model including parallax and orbital motion, -0.01530.0004.

4.9 Second Order Effects

With the slope parameter included we had accounted for all the physical effects which we expected to be present in the lightcurve. Having found that the residuals showed no further variation at a level detectable above the photometric noise, we did not attempt to include second order effects such as xallarap etc.

4.10 Final Model

All our models are compared in Table 2 and the parameters of the best-fitting models are presented in Table 3 and Figure 1. We plot the  for each link in the MCMC chain for all parameters against one another in the best-fitting model in Figure 5. This plot was used as a diagnostic throughtout the fitting process, as any correlations between parameters display distinct trends as the chain moves towards the minimum. The caustic structure changed during the course of event, so in Figure 6 we show the structure at two distinct times; the first at the time of the first caustic crossing during the anomaly, and the second at the time of closest approach.

Model
Standard 6331.950 6331.950
Standard+Parallax 6064.773 6099.487
Standard+Parallax+Orbital Motion 5544.071 5690.003
Standard+Parallax+Orbital Motion with Slope 4782.367 4802.606
Standard+Slope 5334.378 5334.378
Table 2: The for the best-fitting model in each class, comparing  and  solutions.
Parameter Standard Standard Standard Standard+Parallax Standard+Parallax
(Units) +Slope +Parallax +Orbital Motion +Orbital Motion+Slope
6331.950 5334.378 6064.773 5544.071 4782.367
1549.583 552.011 1282.406 761.704 0.0
 (HJD) 5344.32 5344.38 5344.47 5344.83 5344.69
0.01 0.01 0.02 0.03 0.02
0.4089 0.403 0.403 0.381 0.386
0.0009 0.001 0.001 0.001 0.001
 (d) 43.82 44.84 43.49 43.4 44.3
0.08 0.09 0.09 0.1 0.1
0.7692 0.7725 0.7717 0.7792 0.7750
0.0005 0.0005 0.0006 0.0007 0.0007
0.0705 0.0677 0.0695 0.0683 0.0654
0.0005 0.0005 0.0006 0.0006 0.0006
0.180 0.171 0.198 0.297 0.221
0.003 0.003 0.003 0.006 0.007
0.01963 0.01912 0.01931 0.0163 0.0165
0.00006 0.00007 0.00008 0.0001 0.0001
0.18 0.96 0.37
0.02 0.04 0.05
-0.124 0.09 0.01
0.007 0.01 0.01
 () 0.53 0.49
0.02 0.02
 () -1.21 -0.37
0.06 0.08
Slope (mag/yr) -0.0160 -0.0153
0.0004 0.0004
aafootnotetext: All timestamps are abbreviated to HJD= HJD-2450000.0.bbfootnotetext: Improvement in  relative to that of the best fitting model.
Table 3: The best fitting parameters.
Figure 5: Chi squared contours plotted as a function of the parameters fitted in the MCMC fit for the best model. The red, orange, green, light blue, dark blue and purple colors indicate the regions with 1–6 (respectively) from the best-fit solution. Inset: Close-up of the coutours for the parallax parameters.
Figure 6: Maps of the caustic structure during the anomaly and at the time of closest approach. The red line and arrow indicate the trajectory of the source in a reference frame centered on the barycenter of the lensing system, while the red dots indicate the position of the source at these times. The green and blue dots indicate the positions of (largest) and at both times (radii of dots not to scale).

5 Physical Parameters

The purpose of this model is to ultimately arrive at the physical parameters of the lens and source, which can be achieved using the known relations between these and the lensing parameters obtained from the modeling.

Chiefly of interest is the mass of the lensing system, , which can be determined explicitly for events where parallax is measurable provided the angular extent of the Einstein radius, , is known from:

(7)

where is the Einstein radius projected from the source onto the observer’s plane. The model parameter represents the angular size of the source star  in units of the angular Einstein radius . We derive this from the crossing time taken for the source to travel behind the lens, where is the relative source-lens proper motion. can then be written as:

(8)

These parameters also yield the distance to the lens,

(9)

which in turn yields the projected separation between the lens components:

(10)

and the relative proper motion between lens and source, when combined with :

(11)

The appreciable lens orbital motion during this event also allows us to test whether the companion object is bound to the primary lensing mass, via the ratio of its kinetic to potential energy:

(12)

where relates the two lens orbital parameters: =(/) + () and where the masses are in units of , the distances in AU and time measured in years.

However, these expressions include two key terms which are as yet, unknown: , the angular source radius and , the distance to the source. In order to extract the physical characteristics of the lens, we therefore turned our attention to the characteristics of the source.

5.1 Source Star

Long-exposure , images were acquired by the CTIO 1.3 m at several epochs which enabled us to plot the color-magnitude diagram for the field including the source star (Fig. 7). By observing the event at different levels of lensing magnification, these data can be incorporated into the model which yields the source and blended light fluxes, ,  for those data, and hence the instrumental magnitudes and colors of the source and blend. But we note that these uncalibrated fluxes also suffer from the high degree of extinction along the line of sight to the Galactic Bulge. To calculate the dereddened color, , and magnitude, , of the source, we needed to calibrate the instrumental fluxes  and  relative to a standard candle.

Fig. 7 clearly shows a locus of stars centered at =15.8210.05 mag, =-0.3500.05 mag. This consists of a clump of red giant stars, for which stellar theory predicts a stable absolute luminosity, varying only slightly with age and chemical composition. Their frequent occurrence makes these objects useful as standard candles. Stanek et al. (1998) established photometric calibrations for red clump magnitudes which were later refined by Alves et al. (2002) using Hipparcos data. Most recently, Nataf et al. (2012) were able to measure the dereddened apparent magnitude of the red clump stars at the Galactocentric distance, = 14.443. By mapping the distances, , to red clump stars in the Galactic Bar as a function of Galactic longitude, they found an apparent viewing angle on the Bar of =40,

(13)

where Nataf et al. (2012) measured to be 8.20 kpc. For the field of MOA-2010-BLG-073 (Galactic coordinates: =4.81030, =-3.50131), we derive =7.48 kpc, and we assumed that the source star lies behind the same amount of dust as the Red Clump stars and at the same distance. Scaling the dereddened apparent magnitude of the red clump stars, , appropriately for the slightly closer distance of the stars in this field, = + , where:

(14)

We found = 0.20 mag, and so the distance modulus to the red clump and the source in this field is = 14.24 mag. Bensby et al. (2011) determined the intrinsic =1.09 mag for red clump stars, so we were able to derive their absolute magnitude of =0.97 mag. Combining these results with the measured and between the source and red clump in the CTIO data, we then derived the dereddened color,  and magnitude  of the source, summarized in Table 4.

Bessell & Brett (1988) provided a relationship between and color indices, and Kervella et al. (2004) related to angular radius for giant and dwarf stars. Having thus determined  and , we derived the physical parameters of the lensing system, which are also summarized in Table 4. The color and large source radius of 14.71.3  implies this star is a K-type giant, which is consistent with the observed photometric variability. Jorissen et al. (1997) found that red giant stars with spectral types later than early-K are all variable, with amplitudes increasing from microvariability to several magnitudes towards cooler temperatures and timescales from days to years. Kiss et al. (2006) notes that irregular photometric variability may be caused by large convection cells, or may actually be the result of a number of simultaneous periodic pulsation modes, and many examples have been identified from time-domain surveys (Wray et al., 2004; Woźniak et al., 2004; Eyer & Blake, 2005; Ciechanowska et al., 2006, e.g.). We note that the star was detected by 2mass (Skrutskie et al., 2006) as source 2MASS~J18101138-2631226 (catalog 2MASS~J18101138-2631226) with colors =0.760.078 mag and =0.2840.079 mag. Although the 2mass field is crowded the star’s PSF is distinct and their photometry for it has the best-quality AAA flag. These colors are consistent with a giant star and with our derived value for =12.345 mag, when we take into account =0.24 mag from the vvv survey (Gonzalez et al., 2012).

Figure 7: The instrumental color-magnitude diagram for the field of view including the lensing source star. The position of the source is marked with a blue star, relative to the center of the Red Giant Clump highlighted with a red diamond. The green circle indicates the color and magnitude of light blended in the CTIO photometry.
Parameter Units Value
as 9.1430.792
mas 0.5570.09
14.71.3
0.160.03
11.02.0
0.170.03
kpc 2.80.4
AU 1.210.16
KE/PE 0.079
Proper motion mas yr 4.600.4
mag 15.5540.007
mag 15.3350.007
mag -0.220.01
mag 15.8210.05
mag -0.3500.05
mag 14.443
mag 1.09
mag 13.976
mag 15.197
mag 1.2210.051
mag 2.852
mag 12.345
(2mass) mag 13.6860.053
(2mass) mag 12.9260.057
(2mass) mag 12.6420.054
Table 4: The physical parameters of the lensing system and source star, derived from the model including slope, parallax and orbital motion, plus color information.

Nataf et al. (2012) explain that their value of our viewing angle of the Galactic Bar is a “soft upper bound” because the distance along the plane to the greatest density of stars along the line-of-sight to a triaxial Bar structure is less than the distance to the structure’s major axis on the far side and greater on the near side. As the physical parameters derived for MOA-2012-BLG-073 are somewhat dependent on the value of our viewing angle of the Galactic Bar and Nataf et al. (2012) quote consistent results with values as low as 25, we explored the potential impact of this on our results. A reduced viewing angle would produce a smaller distance to the source, changing it’s dereddened magnitude and color. The resulting increase in source radius produces a corresponding reduction in the value of and increases in the lens masses. However we found that the physical parameter values do not change by more than the errors quoted in Table 4, implying that this is not the dominant source of uncertainty. Finally, we computed the physical parameters derived from the best-fitting model for comparison and found that the masses derived changed by .

6 Discussion

The Earth’s movement during this relatively long time scale (=44.3 d) microlensing event resulted in a gradual shift in our perspective on the lensing system, breaking the symmetry of the caustic. Meanwhile, the change in projected separation of the lensing objects modified the shape of the caustic just as the source’s trajectory happened to pass close by. If not for these subtle variations, it can be seen from Figure 6 that a source trajectory  would produce exactly the same light curve as a  trajectory. As it is, for this event, the  solution best explains our observations, though the difference in  relative to the corresponding  solution is very small (=20.2) compared with the  of both fits.

The measurable parallax signature enables us to determine the masses of the lensing bodies. The primary lensing object has =0.160.03 , making it an M-dwarf star. The companion’s mass is =11.02.0 . This places it the brown dwarf desert, though we note that this traditionally refers to close-in companions, and since microlensing and direct imaging measure only the projected separation, we know only their minimum orbital semi-major radii. Regardless, it is clear that MOA-2010-BLG-073L b is close to the mass threshold for deuterium burning (0.012 =12.6 ) quoted as the nominal boundary between planets and brown dwarfs established by the IAU (Chabrier et al., 2005). So what kind of object is it?

Figure 8: The ratio of planet mass to host star mass plotted against semi-major axis for all exoplanets for which these parameters are available (source: exoplanet.eu Schneider et al. (2011)). The true orbital semi-major axis is plotted where it has been measured, otherwise the projected (minimum) separation has been indicated with an arrow, where the base of the arrow marks the measured value. Very few objects have been found between =-1 – -2. This corresponds to the brown dwarf desert. The location of MOA-2010-BLG-073L, highlighted in cyan and bold, lies between the brown dwarf and planetary regimes.

No further orbital or metallicity information is available for this event, which might have shed light on its evolutionary history. Theoretical isochrones predict that a star of this low mass will not have lost a significant amount of material over its lifetime, so we can say that MOA-2010-BLG-073L b formed as a high mass-ratio binary. Models of protoplanetary disks have the expectation that disk mass, , will scale linearly with star mass (Williams & Cieza, 2011), % at young ages, but this would limit – and hence – to 1–10  in the case of MOA-2010-BLG-073L. So it seems questionable whether such a massive companion could have formed in a protostellar disk, via either core accretion or gravitational instability.

Bonnell et al. (2008); Kroupa & Bouvier (2003) discuss a number of mechanisms which can produce an M-dwarf/Brown Dwarf binary following gravitational fragmentation in a molecular cloud:

  1. Embryo rejection model the nascent binary was ejected from a dynamically unstable multiple protostellar system, leading to the loss of its accretion envelope before the secondary component could acquire enough mass to become a star.

  2. Collision model the binary was prematurely ejected from a larger protostellar system by the close passage of another star.

  3. Photo-evaporation model the accretion envelope around the binary was photo-evaporated by the nearby presence of a massive star in the birth cluster before the secondary could accrete enough mass to become a star.

  4. “Star-like” model the object formed as a normal stellar binary with low-mass components.

The embryo rejection model predicts (Bate, Bonnell & Bromm, 2002) that the maximum separation of binaries surviving this process is  AU. We cannot rule out this scenario as we only measure the projected separation of the lens, which is nevertheless below . Since MOA-2010-BLG-073L is a field object we have no information regarding the proximity of other stars during its birth, so the collisional and photo-evapouration models are equally plausible. However, we note that Whitworth & Stamatellos (2006) derived a minimum mass for primary fragmentation components of 0.004 4.2 . As this threshold is below the mass of MOA-2010-BLG-073Lb the simplest explanation is that the companion is an extremely low-mass product of star formation. However, we note that the mass ratio of this system would match that of an object in the brown dwarf desert if the host were a more massive star. Recent surveys (Dieterich et al. (2012, e.g.); Evans et al. (2012, e.g.), though restricted to massive companions with  AU) hint that the brown dwarf desert may extend to M-dwarfs and beyond 3 AU, which would make MOA-2010-BLG-073L a rare example of its type.

RAS is grateful to the Chungbuk University group and especially to C. Han and J.-Y. Choi for their advice and hospitality in Korea where much of this work was completed. RAS would also like to express appreciation for Y. Tsapras, K. Horne, M. Hundertmark, S. Dong and P. Fouqué for many useful discussions. KA, DMB, MD, KH, MH, CL, CS, RAS and YT would like to thank the Qatar Foundation for support from QNRF grant NPRP-09-476-1-078. Work by C. Han was supported by the Creative Research Initiative Program (2009-0081561) of the National Research Foundation of Korea. CS received funding from the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 268421. The OGLE project has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013) / ERC grant agreement no. 246678 to AU. This publication makes use of data products from the Exoplanet Encyclopeida and the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. MOA acknowledge funding JSPS20340052, JSPS22403003 and JSPS19340058. TCH gratefully acknowledges financial support from the Korea Research Council for Fundamental Science and Technology (KRCF) through the Young Research Scientist Fellowship Program. CUL and TCH acknowledge financial support from KASI (Korea Astronomy and Space Science Institute) grant number 2012-1-410-02. DR (boursier FRIA), FF (boursier ARC) and J. Surdej acknowledge support from the Communauté francaise de Belgique - Actions de recherche concertées - Académie universitaire Wallonie-Europe. A. Gould acknowledges support from NSF AST-1103471. B.S. Gaudi, A. Gould and R.W. Pogge acknowledge support from NASA grant NNX12AB99G. Work by J.C. Yee is supported by a National Science Foundation Graduate Research Fellowship under Grant No. 2009068160. S.D. is supported through a Ralph E. and Doris M. Hansmann Membership at the IAS and NSF grant AST-0807444. Facilities: FTN, FTN, Liverpool: 2m, CTIO: 1.3m, MtCanopus: 1m, SAAO: 1m, Danish 1.54m Telescope

References

  • Albrow et al. (1999) Albrow, M.D., Beaulieu, J.-P., Caldwell, J.A.R. et al. 1999, ApJ, 522, 1022.
  • Albrow et al. (2009) Albrow, M. D., Horne, K., Bramich, D. M. et al. 2009, MNRAS, 397, 2099
  • Alves et al. (2002) Alves, D. R., Rejkuba, M., Minniti, D.,Cook, K. H. 2002, ApJ, 573, L51
  • Bachelet et al. (2012) Bachelet, E., Shin, I.-G., Han.C et al. 2012, ApJ, 754, 73.
  • Beaulieu et al. (2006) Beaulieu, J.-P., Bennett, D.P.,Fouqué, P. et al. 2006, Nature, 439, 437.
  • Baraffe et al. (2010) Baraffe, I., Chabrier, G.,Barman,T. 2010, Rep. Prog. Phys. 73, 016901.
  • Bate, Bonnell & Bromm (2002) Bate, M.R., Bonnell, I.A., Bromm, V. 2002, MNRAS, 332, L65.
  • Batista et al. (2011) Batista, V., Gould, A., Dieters, S. etal. 2011, A&A, 529, 102.
  • Béjar et al. (2012) Béjar, V.J.S., Zapatero Osorio,M.R., Rebolo, R. et al. 2012, ApJ, 743, 64.
  • Bensby et al. (2011) Bensby, T., Adén, D., Meléndez, J. et al. 2011, A&A, 533, 134.
  • Bessell & Brett (1988) Bessell, M.S. & Brett, J.M.1988, PASP, 100, 1134.
  • Bond et al. (2001) Bond, I., Abe, F., Dodd, R. J. et al. 2001, MNRAS, 327, 868.
  • Bonnell et al. (2008) Bonnell, I.A., Clark, P., Bate, M.R.2008, MNRAS, 389, 1556.
  • Bonfils et al. (2011) Bonfils, X., Delfosse, X., Udry, S. et al.2011, A&A, submitted.
  • Boss (2006) Boss, A.P. 2006, ApJ, 643, 501.
  • Bramich (2008) Bramich, D.M. 2008, MNRAS, 386, L77.
  • Burrows et al. (2001) Burrows, A., Hubbard, J.I.,Lunine, J.I. et al. 2001, Rev. Mod. Phys., 73, 719.
  • Chabrier et al. (2005) Chabrier, G., Baraffe, I., Allard,F. et al. 2005, ASP Conf. Series ”Resolved Stellar Populations”, Cancun.
  • Chabrier et al. (2011) Chabrier, G., Leconte, J.,Baraffe, I. 2011, Proc. IAU Symp., 276, 171.
  • Ciechanowska et al. (2006) Ciechanowska, A. et al. 2006, Acta Astron., 56, 219.
  • Cumming et al. (2008) Cumming, A., Butler, R.P., Marcy, G.W. et al. 2008, PASP, 120, 531.
  • Díaz et al. (2012) Díaz, R.F., Santerne, A., Sahlmann, J. et al. 2012, A&A, 538, A113.
  • Dieterich et al. (2012) Dieterich, S.B., Henry, T.J., Golimowski, D.A. et al. 2012, AJ, 144, 64.
  • Dominik et al. (2010) Dominik, M., Jørgensen, U.G., Rattenbury, N.J. et al. 2010, Ast. Nac., 331, 671.
  • Dong et al. (2006) Dong, S., DePoy, D. L., Gaudi, B. S. et al. 2006, ApJ, 642, 842.
  • Dong et al. (2009) Dong, S., Gould, A., Udalski, A. et al. 2009, ApJ, 695, 970.
  • Evans et al. (2012) Evans, T.M., Ireland, M.J., Kraus, A.L. et al. 2012, ApJ, 744, 120.
  • Eyer & Blake (2005) Eyer, L. & Blake, C. 2005, MNRAS, 358, 30.
  • Fischer & Valenti (2005) Fischer, D.A. & Valenti, J. 2005, ApJ, 622, 1102.
  • Forveille et al. (2011) Forveille, T., Bonfils, X., Lo Curto, G. et al. 2011, A&A, 526, A141.
  • Gaudi (2009) Gaudi, B.S. 2009, Oral presentation at 13 Microlensing Workshop, Paris.
  • Gonzalez et al. (2012) Gonzalez, O.A., Rejkuba, M., Zoccali, M. 2012, A&A, 543, A13.
  • Gould et al. (2006) Gould, A., Udalski, A., An, D. et al. 2006, ApJ, 644, L37.
  • Gould (2008) Gould, A. 2008, ApJ, 681, 1598.
  • Heyrovský (2007) Heyrovský, D. 2007, ApJ, 656, 483.
  • Howard et al (2012) Howard, A.W., Marcy, G.W., Bryson, S.T. et al. 2012, ApJ, 201, 15.
  • Ida & Lin (2004) Ida, S. & Lin, D.N.C. 2004, ApJ, 616, 567.
  • Ida & Lin (2005) Ida, S. & Lin, D.N.C. 2005, ApJ, 626, 1045.
  • Johnson et al. (2007) Johnson, J.A., Butler, R.P., Marcy, G.W. et al. 2007, ApJ, 670, 833.
  • Johnson et al. (2010) Johnson, J.A., Aller, K.M., Howard, A.W. et al. 2010, PASP, 122, 905.
  • Jorissen et al. (1997) Jorissen, A., Mowlavi, N., Sterken, C., Manifroid, J. 1997, A&A, 324, 578.
  • Kervella et al. (2004) Kervella, P., Bersier, D., Mourard, D. et al. 2004, A&A, 428, 587.
  • Kayser et al. (1986) Kayser, R., Refsdal, S. & Stabell, R.1986, A&A, 166, 36.
  • Kirkpatrick et al. (2006) Kirkpatrick, J. D., Barman, T. S., Burgasser, A. et al. 2006, ApJ, 639, 1120.
  • Kiss et al. (2006) Kiss, L.L., Szabó, Gy.M., Bedding, T.R. 2006, MNRAS, 372, 1721.
  • Klein et al. (2003) Klein, R., Apai, D., Pascucci, I. et al. 2003, ApJ, 593, L57.
  • Kroupa & Bouvier (2003) Kroupa, P. & Bouvier, J. 2003,MNRAS, 346, 369.
  • Kurucz (1979) Kurucz, R.L. 1979, ApJS, 40, 1.
  • Laughlin et al. (2004) Laughlin, G., Bodenheimer, P. & Adams,F.C. 2004, ApJ, 612, L73.
  • Luhman (2012) Luhman, K.L. 2012, ARA&A, 50, 65.
  • Marcy & Butler (2000) Marcy, G.W. & Butler, R.P. 2000, PASP, 112, 137.
  • Maldonado et al. (2012) Maldonado, J., Eiroa, C., Villaver, E. et al. 2012, A&A, 541, 40.
  • Marois et al. (2008) Marois, C., MacIntosh, B., Barman, T., 2008, Science, 322, 1348
  • Mordasini et al. (2009) Mordasini, C., Alibert, Y., Benz, W. 2009,A&A, 501, 1139.
  • Miyake et al. (2012) Miyake, N., Udalski, A., Sumi, T. et al. 2012, ApJ, 752, 82.
  • Nataf et al. (2012) Nataf, D.M., Gould, A., Fouqué, P. et al. 2012, ApJ,submitted.
  • Park et al. (2004) Park, B.-G., DePoy, D.L., Gaudi, B.S. et al. 2004, ApJ, 609, 166.
  • Pejcha & Heyrovský (2009) Pejcha, O. & Heyrovský, D. 2009, ApJ, 690, 1772.
  • Penny et al. (2010) Penny, M.T., Mao, S. & Kerins, E. 2010, MNRAS, 412, 607.
  • Santos, Israelian & Mayor (2001) Santos, N.C., Israelian, G. & Mayor, M. 2001, A&A, 373, 1019.
  • Sahlmann et al. (2010) Sahlmann, J., Ségransan, D., Queloz, D. et al. 2010, A&A, 525, A95.
  • Schneider et al. (2011) Schneider, J., Dedieu, C., Le Sidaner, P. et al. 2011, A&A, 532, A79.
  • Scholz et al. (2006) Scholz, A., Jayawardhana, R. & Wood, K. 2006, ApJ, 645, 1498.
  • Schwarzenberg-Czerny (1996) Schwarzenberg-Czerny, A. 1996,ApJ, 460, L107–L110.
  • Shin et al. (2012a) Shin, I.-G., Han, C., Gould, A. et al. 2012, ApJ, in prep.
  • Shin et al. (2012b) Shin, I.-G., Han, C., Choi, J.-Y. et al. 2012, ApJ, 746, 127.
  • Skowron et al. (2011) Skowron, J., Udalski, A., Gould, A. et al. 2011, ApJ, 738, 87.
  • Skrutskie et al. (2006) Skrutskie, M.F., Cutri, R.M., Stiening, R. et al. 2006, AJ, 131, 1163.
  • Stanek et al. (1998) Stanek, K.Z. & Garnavich, P.M. 1998, ApJ, 503, 131.
  • Sumi et al. (2003) Sumi, T., Bond, I. A., Dodd, R. J. 2003, ApJ, 591, 204.
  • Sumi et al. (2011) Sumi, T., Kamiya, K., Bennett, D.P. et al. 2011,Nature, 473, 349.
  • Tsapras et al. (2009) Tsapras, Y, Street, R., Horne, K. et al. 2009, Ast. Nac., 330,4.
  • Udalski et al. (2003) Udalski, A. 2003, Acta Astron.,53, 291.
  • Whitworth & Stamatellos (2006) Whitworth, A.P. & Stamatellos, D. 2006, A&A, 458, 817.
  • Williams & Cieza (2011) Williams, J.P. & Cieza, L.A. 2011, ARA&A, 49, 67.
  • Woźniak et al. (2004) Woźniak, P.R., Williams, S.J., Vestrand, W.T., Gupta, V. 2004, AJ, 128, 2965.
  • Wray et al. (2004) Wray, J.J., Eyer, L., Paczyński, B. 2004,MNRAS, 349, 1059.
  • Yee et al. (2012) Yee, J. C., Shvartzvald, Y., Gal-Yam, A. et al. 2012, ApJ, 755, 102.
  • Zhou et al. (2012) Zhou, J.-L., Xie, J-W, Liu, H-G et al. 2012, Res. in Astron. & Astrophy., 12, 1081.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
138518
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description