Abstract We present pyLIMASS, a novel algorithm for estimating the physical properties of the lensing system in microlensing events. The main idea of pyLIMASS is to combine all available information ...regarding the microlensing event, defined as observables , and to estimate the parameter distributions of the system, such as the lens mass and distance. The algorithm is based on isochrones for the stars model and combines the observables using a Gaussian mixture approach. After describing the mathematical formalism and its implementation, we discuss the algorithm’s performance on simulated and published events. Generally, the pyLIMASS estimations are in good agreement (i.e., within 1 σ ) with the results of the selected published events, making it an effective tool to estimate the lens properties and their distribution. The applicability of the method was tested by using a catalog of realistically simulated events that could be observed by the future Galactic Bulge Time Domain Survey of the Nancy Grace Roman Space Telescope. By solely using constraints from the Roman lightcurves and images, pyLIMASS estimates the masses of the lens of the Roman catalog with a median precision of 20% with almost no bias.
Abstract
Galactic science encompasses a wide range of subjects in the study of the Milky Way and Magellanic Clouds, from young stellar objects to X-ray binaries. Mapping these populations, and ...exploring transient phenomena within them, are among the primary science goals of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time. While early versions of the survey strategy dedicated relatively few visits to the Galactic Plane region, more recent strategies under consideration envision a higher cadence within selected regions of high scientific interest. The range of galactic science presents a challenge in evaluating which strategies deliver the highest scientific returns. Here we present metrics designed to evaluate Rubin survey strategy simulations, based on the cadence of observations they deliver within regions of interest to different topics in galactic science, using variability categories defined by timescale. We also compare the fractions of exposures obtained in each filter with those recommended for the different science goals. We find that the
baseline
_
v2.x
simulations deliver observations of the high-priority regions at sufficiently high cadence to reliably detect variability on timescales >10 days or more. Follow-up observations may be necessary to properly characterize variability, especially transients, on shorter timescales. Combining the regions of interest for all the science cases considered, we identify those areas of the Galactic Plane and Magellanic Clouds of highest priority. We recommend that these refined survey footprints be used in future simulations to explore rolling cadence scenarios, and to optimize the sequence of observations in different bandpasses.
We present a general framework for matching the point-spread function (PSF), photometric scaling and sky background between two images, a subject which is commonly referred to as difference image ...analysis (DIA). We introduce the new concept of a spatially varying photometric scale factor which will be important for DIA applied to wide-field imaging data in order to adapt to transparency and airmass variations across the field-of-view. Furthermore, we demonstrate how to separately control the degree of spatial variation of each kernel basis function, the photometric scale factor and the differential sky background. We discuss the common choices for kernel basis functions within our framework, and we introduce the mixed-resolution delta basis functions to address the problem of the size of the least-squares problem to be solved when using delta basis functions. We validate and demonstrate our algorithm on simulated and real data. We also describe a number of useful optimizations that may be capitalized on during the construction of the least-squares matrix and which have not been reported previously. We pay special attention to presenting a clear notation for the DIA equations which are set out in a way that will hopefully encourage developers to tackle the implementation of DIA software.
We report the analysis of additional multiband photometry and spectroscopy and new adaptive optics (AO) imaging of the nearby planetary microlensing event TCP J05074264+2447555 (Kojima-1), which was ...discovered toward the Galactic anticenter in 2017 (Nucita et al.). We confirm the planetary nature of the light-curve anomaly around the peak while finding no additional planetary feature in this event. We also confirm the presence of apparent blending flux and the absence of significant parallax signal reported in the literature. The AO image reveals no contaminating sources, making it most likely that the blending flux comes from the lens star. The measured multiband lens flux, combined with a constraint from the microlensing model, allows us to narrow down the previously unresolved mass and distance of the lens system. We find that the primary lens is a dwarf on the K/M boundary (0.581 0.033 M ) located at 505 47 pc, and the companion (Kojima-1Lb) is a Neptune-mass planet (20.0 2.0 M⊕) with a semimajor axis of au. This orbit is a few times smaller than those of typical microlensing planets and is comparable to the snow-line location at young ages. We calculate that the a priori detection probability of Kojima-1Lb is only ∼35%, which may imply that Neptunes are common around the snow line, as recently suggested by the transit and radial velocity techniques. The host star is the brightest among the microlensing planetary systems (Ks = 13.7), offering a great opportunity to spectroscopically characterize this system, even with current facilities.
Abstract During the last 25 yr, hundreds of binary stars and planets have been discovered toward the Galactic bulge by microlensing surveys. Thanks to a new generation of large-sky surveys, it is now ...possible to regularly detect microlensing events across the entire sky. The OMEGA Key Projet at the Las Cumbres Observatory carries out automated follow-up observations of microlensing events alerted by these surveys with the aim of identifying and characterizing exoplanets as well as stellar remnants. In this study, we present the analysis of the binary lens event Gaia20bof. By automatically requesting additional observations, the OMEGA Key Project obtained dense time coverage of an anomaly near the peak of the event, allowing characterization of the lensing system. The observed anomaly in the lightcurve is due to a binary lens. However, several models can explain the observations. Spectroscopic observations indicate that the source is located at ≤2.0 kpc, in agreement with the parallax measurements from Gaia. While the models are currently degenerate, future observations, especially the Gaia astrometric time series as well as high-resolution imaging, will provide extra constraints to distinguish between them.
As the
Kepler
mission has done for hot exoplanets, the ESA
Euclid
and NASA
Roman
missions have the potential to create a breakthrough in our understanding of the demographics of cool exoplanets, ...including unbound, or free-floating, planets (FFPs).
Roman
will dedicate part of its core survey program to the detection of cool exoplanets via microlensing, while
Euclid
may undertake a microlensing program as an ancillary science goal. In this study, we demonstrate the complementarity of the two missions and propose two joint surveys to better constrain the mass and distance of microlensing events. We first demonstrate that an early brief
Euclid
survey (~7 h) of the
Roman
microlensing fields will allow the measurement of at least 30% of the events’ relative proper motions
µ
rel
and 42% of the lens magnitudes. This survey would place strong constraints on the mass and distance on thousands of microlensing events observed by
Roman
just after the first year of observation. Then, we study the potential of simultaneous observations by
Roman
and
Euclid
to enable the measurement of the microlensing parallax for the shortest microlensing events and, ultimately, obtain a direct measurement of the masses, distances, and transverse motions of FFPs. Using detailed simulations of the joint detection yield we show that within one year
Roman-Euclid
observations will be at least an order of magnitude more sensitive than current ground-based measurements. The recent tentative detection of an excess of short-duration events by the OGLE survey is consistent with a scenario of up to ten Earth-mass FFPs per Galactic star. For such a scenario a joint
Roman-Euclid
campaign should detect around 130 FFP events within a year, including 110 with measured parallax that strongly constrain the FFP mass, and around 30 FFP events with direct mass and distance measurements. The ability of the joint survey to completely break the microlens mass-distance-velocity degeneracy for a significant subset of events provides a unique opportunity to verify unambiguously the FFP hypothesis or else place abundance limits for FFPs between Earth and Jupiter masses that are up to two orders of magnitude stronger than provided by ground-based surveys. Finally, we study the capabilities of the joint survey to enhance the detection and characterization of exomoons, and find that it could lead to the detection of the first exomoon.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
7.
The Thresher: Lucky imaging without the waste Hitchcock, J A; Bramich, D M; Foreman-Mackey, D ...
Monthly notices of the Royal Astronomical Society,
03/2022, Volume:
511, Issue:
4
Journal Article
Peer reviewed
Open access
ABSTRACT
In traditional lucky imaging (TLI), many consecutive images of the same scene are taken with a high frame-rate camera, and all but the sharpest images are discarded before constructing the ...final shift-and-add image. Here, we present an alternative image analysis pipeline – The Thresher – for these kinds of data, based on online multi-frame blind deconvolution. It makes use of all available data to obtain the best estimate of the astronomical scene in the context of reasonable computational limits; it does not require prior estimates of the point-spread functions in the images, or knowledge of point sources in the scene that could provide such estimates. Most importantly, the scene it aims to return is the optimum of a justified scalar objective based on the likelihood function. Because it uses the full set of images in the stack, The Thresher outperforms TLI in signal-to-noise ratio; as it accounts for the individual-frame PSFs, it does this without loss of angular resolution. We demonstrate the effectiveness of our algorithm on both simulated data and real Electron-Multiplying CCD images obtained at the Danish 1.54-m telescope (hosted by ESO, La Silla). We also explore the current limitations of the algorithm, and find that for the choice of image model presented here, non-linearities in flux are introduced into the returned scene. Ongoing development of the software can be viewed at https://github.com/jah1994/TheThresher.
We present 17 transit light curves of the ultrashort period planetary system WASP-103, a strong candidate for the detection of tidally-induced orbital decay. We use these to establish a ...high-precision reference epoch for transit timing studies. The time of the reference transit mid-point is now measured to an accuracy of 4.8 s, versus 67.4 s in the discovery paper, aiding future searches for orbital decay. With the help of published spectroscopic measurements and theoretical stellar models, we determine the physical properties of the system to high precision and present a detailed error budget for these calculations. The planet has a Roche lobe filling factor of 0.58, leading to a significant asphericity; we correct its measured mass and mean density for this phenomenon. A high-resolution Lucky Imaging observation shows no evidence for faint stars close enough to contaminate the point spread function of WASP-103. Our data were obtained in the Bessell RI and the SDSS griz passbands and yield a larger planet radius at bluer optical wavelengths, to a confidence level of 7.3σ. Interpreting this as an effect of Rayleigh scattering in the planetary atmosphere leads to a measurement of the planetary mass which is too small by a factor of 5, implying that Rayleigh scattering is not the main cause of the variation of radius with wavelength.
ABSTRACT
Transits in the planetary system WASP-4 were recently found to occur 80 s earlier than expected in observations from the TESS satellite. We present 22 new times of mid-transit that confirm ...the existence of transit timing variations, and are well fitted by a quadratic ephemeris with period decay dP/dt = −9.2 ± 1.1 ms yr−1. We rule out instrumental issues, stellar activity, and the Applegate mechanism as possible causes. The light-time effect is also not favoured due to the non-detection of changes in the systemic velocity. Orbital decay and apsidal precession are plausible but unproven. WASP-4 b is only the third hot Jupiter known to show transit timing variations to high confidence. We discuss a variety of observations of this and other planetary systems that would be useful in improving our understanding of WASP-4 in particular and orbital decay in general.
Planet population synthesis models predict an abundance of planets with semimajor axes between 1 and 10 au, yet they lie at the edge of the detection limits of most planet finding techniques. ...Discovering these planets and studying their distribution is critical to understanding the physical processes that drive planet formation. ROME/REA is a gravitational microlensing project whose main science driver is to discover exoplanets in the cold outer regions of planetary systems. To achieve this, it uses a novel approach combining a multiband survey with reactive follow-up observations, exploiting the unique capabilities of the Las Cumbres Observatory global network of robotic telescopes combined with a Target and Observation Manager system. We present the main science objectives and a technical overview of the project, including initial results.
Full text
Available for:
BFBNIB, NMLJ, NUK, PNG, UL, UM, UPUK