We have used photometric data on almost 91 000 fundamental-mode RR Lyrae stars (type RRab) detected by the OGLE survey to investigate properties of old populations in the Milky Way. Based on their metallicity distributions, we demonstrate that the Galaxy is built from three distinct old components: halo, bulge, and disk. The distributions reach their maxima at approximately [Fe/H]_J95 = -1.2, -1.0, and -0.6 dex on the Jurcsik's metallicity scale, respectively. We find that, very likely, the entire halo is formed from infalling dwarf galaxies. It is evident that halo stars penetrate the inner regions of the Galactic bulge. We estimate that about one-third of all RR Lyr stars within the bulge area belong in fact to the halo population. The whole old bulge is dominated by two populations, A and B, represented by a double sequence in the period-amplitude (Bailey) diagram. The boundary in iron abundance between the halo and the disk population is at about [Fe/H]_J95 = -0.8 dex. Using Gaia DR2 for RRab stars in the disk area, we show that the observed dispersion of proper motions along the Galactic latitude decreases smoothly with the increasing metal content excluding a bump around [Fe/H]_J95 = -1.0 dex.

We analyze KMT-2019-BLG-1339, a microlensing event with an obvious but incompletely resolved brief anomaly feature around the peak of the light curve. Although the origin of the anomaly is identified to be a companion to the lens with a low mass ratio $q$, the interpretation is subject to two different degeneracy types. The first type is the ambiguity in $\rho$, representing the angular source radius scaled to the angular radius of the Einstein ring, $\theta_{\rm E}$, and the other is the $s\leftrightarrow s^{-1}$ degeneracy. The former type, `finite-source degeneracy', causes ambiguities in both $s$ and $q$, while the latter induces an ambiguity only in $s$. Here $s$ denotes the separation (in units of $\theta_{\rm E}$) in projection between the lens components. We estimate that the lens components have masses $(M_1, M_2)\sim (0.27^{+0.36}_{-0.15}~M_\odot, 11^{+16}_{-7}~M_{\rm J})$ and $\sim (0.48^{+0.40}_{-0.28}~M_\odot, 1.3^{+1.1}_{-0.7}~M_{\rm J})$ according to the two solutions subject to the finite-source degeneracy, indicating that the lens comprises an M dwarf and a companion with a mass around the planet/brown dwarf boundary or a Jovian-mass planet. It is possible to lift the finite-source degeneracy by conducting future observations utilizing a high resolution instrument because the relative lens-source proper motion predicted by the solutions are widely different.

We report the discovery of a planet in the microlensing event OGLE-2018-BLG-1269, with planet-host mass ratio $q \sim 6\times10^{-4}$, i.e., $0.6$ times smaller than the Jupiter/Sun mass ratio. Combined with the $Gaia$ parallax and proper motion, a strong one-dimensional constraint on the microlens parallax vector allows us to significantly reduce the uncertainties of lens physical parameters. A Bayesian analysis that ignores any information about light from the host yields that the planet is a cold giant $(M_{2} = 0.69_{-0.22}^{+0.44}\,M_{\rm J})$ orbiting a Sun-like star $(M_{1} = 1.13_{-0.35}^{+0.72}\,M_{\odot})$ at a distance of $D_{\rm L} = 2.56_{-0.62}^{+0.92}\,{\rm kpc}$. The projected planet-host separation is $a_{\perp} = 4.61_{-1.17}^{+1.70}\,{\rm au}$. Using {\it Gaia} astrometry, we show that the blended light lies $\lesssim 12\,$mas from the host and therefore must be either the host star or a stellar companion to the host. An isochrone analysis favors the former possibility at $>99.6\%$. The host is therefore a subgiant. For host metallicities in the range of $0.0 \leq {\rm [Fe/H]} \leq +0.3$, the host and planet masses are then in the range of $1.16 \leq M_{1}/M_{\odot} \leq 1.38$ and $0.74 \leq M_{2}/M_{\rm J} \leq 0.89$, respectively. Low host metallicities are excluded. The brightness and proximity of the lens make the event a strong candidate for spectroscopic followup both to test the microlensing solution and to further characterize the system.

Searches for gravitational microlensing events are traditionally concentrated on the central regions of the Galactic bulge but many microlensing events are expected to occur in the Galactic plane, far from the Galactic Center. Owing to the difficulty in conducting high-cadence observations of the Galactic plane over its vast area, which are necessary for the detection of microlensing events, their global properties were hitherto unknown. Here, we present results of the first comprehensive search for microlensing events in the Galactic plane. We searched an area of almost 3000 square degrees along the Galactic plane (|b|<7, 0<l<50, 190<l<360 deg) observed by the Optical Gravitational Lensing Experiment (OGLE) during 2013-2019 and detected 630 events. We demonstrate that the mean Einstein timescales of Galactic plane microlensing events are on average three times longer than those of Galactic bulge events, with little dependence on the Galactic longitude. We also measure the microlensing optical depth and event rate as a function of Galactic longitude and demonstrate that they exponentially decrease with the angular distance from the Galactic Center (with the characteristic angular scale length of 32 deg). The average optical depth decreases from $0.5\times 10^{-6}$ at l=10 deg to $1.5\times 10^{-8}$ in the Galactic anticenter. We also find that the optical depth in the longitude range 240<l<330 deg is asymmetric about the Galactic equator, which we interpret as a signature of the Galactic warp.

We present the most extensive and detailed reddening maps of the Magellanic Clouds (MCs) derived from the color properties of Red Clump (RC) stars. The analysis is based on the deep photometric maps from the fourth phase of the Optical Gravitational Lensing Experiment (OGLE-IV), covering approximately 670 deg2 of the sky in the Magellanic System region. The resulting maps provide reddening information for 180 deg2 in the Large Magellanic Cloud (LMC) and 75 deg2 in the Small Magellanic Cloud (SMC), with a resolution of 1.7x1.7 arcmin in the central parts of the MCs, decreasing to approximately 27x27 arcmin in the outskirts. The mean reddening is E(V-I) = 0.100 +- 0.043 mag in the LMC and E(V-I) = 0.049 +- 0.024 mag in the SMC. We refine methods of calculating the RC color to obtain the highest possible accuracy in determining the intrinsic color of RC stars. We show that the intrinsic RC color in the LMC is $(V-I)_0 = 0.876 +- 0.013 mag and there is no color gradient in the verifiable range, i.e. from 3 up to 9 deg from the galaxy center. In the case of the SMC (V-I)_0 = 0.856 +- 0.016 mag and is constant between 1 and 4 deg from the galaxy center, followed by a slight decline beyond 4 deg from the center of the SMC. The reddening map of the MCs is available on-line both in the downloadable form and as an interactive interface.

High-cadence observations of the Galactic bulge by the microlensing surveys led to the discovery of a handful of extremely short-timescale microlensing events that can be attributed to free-floating or wide-orbit planets. Here, we report the discovery of another strong free-floating planet candidate, which was found from the analysis of the gravitational microlensing event OGLE-2019-BLG-0551. The light curve of the event is characterized by a very short duration (<3 d) and a very small amplitude (< 0.1 mag). From modeling of the light curve, we find that the Einstein timescale, tE = 0.381 +/- 0.017 d, is much shorter, and the angular Einstein radius, thetaE = 4.35 +/- 0.34 uas, is much smaller than those of typical lensing events produced by stellar-mass lenses (tE ~ 20 d, thetaE ~ 0.3 mas), indicating that the lens is very likely to be a planetary-mass object. We conduct an extensive search for possible signatures of a companion star in the light curve of the event, finding no significant evidence for the putative host star. For the first time, we also demonstrate that the angular Einstein radius of the lens does not depend on blending in the low-magnification events with strong finite source effects.

We report the analysis of the microlensing event OGLE-2018-BLG-0677. A small feature in the light curve of the event leads to the discovery that the lens is a star-planet system. Although there are two degenerate solutions that could not be distinguished for this event, both lead to a similar planet-host mass ratio. We perform a Bayesian analysis based on a Galactic model to obtain the properties of the system and find that the planet corresponds to a super-Earth/sub-Neptune with a mass $M_{\mathrm{planet}} = {3.96}^{+5.88}_{-2.66}\mathrm{M_\oplus}$. The host star has a mass $ M_{\mathrm{host}} = {0.12}^{+0.14}_{-0.08}\mathrm{M_\odot}$. The projected separation for the inner and outer solutions are ${0.63}^{+0.20}_{-0.17}$~AU and ${0.72}^{+0.23}_{-0.19}$~AU respectively. At $\Delta\chi^2=\chi^2({\rm 1L1S})-\chi^2({\rm 2L1S})=46$, this is by far the lowest $\Delta\chi^2$ for any securely-detected microlensing planet to date, a feature that is closely connected to the fact that it is detected primarily via a "dip" rather than a "bump".

We present an upgrade of the OGLE Collection of RR Lyrae stars in the Galactic bulge and disk. The size of our sample has been doubled and reached 78 350 RR Lyr variables, of which 56 508 are fundamental-mode pulsators (RRab stars), 21 321 pulsate solely in the first-overtone (RRc stars), 458 are classical double-mode pulsators (RRd stars), and 63 are anomalous RRd variables (including five triple-mode pulsators). For all the newly identified RR Lyr stars, we publish time-series photometry obtained during the OGLE Galaxy Variability Survey. We present the spatial distribution of RR Lyr stars on the sky, provide a list of globular clusters hosting RR Lyr variables, and discuss the Petersen diagram for multimode pulsators. We find new RRd stars belonging to a compact group in the Petersen diagram (with period ratios P_1O/P_F = 0.74 and fundamental-mode periods P_F = 0.44 d) and we show that their spatial distribution is roughly spherically symmetrical around the Milky Way center.

We present an updated three dimensional map of the Milky Way based on a sample of 2431 classical Cepheid variable stars, supplemented with about 200 newly detected classical Cepheids from the OGLE survey. The new objects were discovered as a result of a dedicated observing campaign of the ~280 square degree extension of the OGLE footprint of the Galactic disk during 2018-2019 observing seasons. These regions cover the main part of the northern Galactic warp that has been deficient in Cepheids so far. We use direct distances to the sample of over 2390 classical Cepheids to model the distribution of the young stellar population in the Milky Way and recalculate the parameters of the Galactic disk warp. Our data show that its northern part is very prominent and its amplitude is ~10% larger than that of the southern part. By combining Gaia astrometric data with the Galactic rotation curve and distances to Cepheids from our sample, we construct a map of the vertical component of the velocity vector for all Cepheids in the Milky Way disk. We find large-scale vertical motions with amplitudes of 10-20 km/s, such that Cepheids located in the northern warp exhibit large positive vertical velocity (toward the north Galactic pole), whereas those in the southern warp - negative vertical velocity (toward the south Galactic pole).

With $q=(4.09\pm 0.27) \times 10^{-5}$, KMT-2019-BLG-0842Lb is the second cold planet with planet-host mass ratio $q$ below that of Uranus relative to the Sun ($q=4.37\times 10^{-5}$), the first having been discovered only in 2018. Hence, such cold, sub-Uranus-ratio planets are probably not rare. We discuss possible improvements in the search for these planets. The Bayesian estimates for the host mass, planet mass, system distance, and planet-host projected separation are $M_{\rm host}=0.76\pm 0.40 M_\odot$, $M_{\rm planet}=10.3\pm 5.5 M_\oplus$, $D_L = 3.3\pm 1.3$ kpc, and $a_\perp = 3.3\pm 1.4$ AU, respectively. The blended light is plausibly due to the host star, an identification that could be tested by high resolution images taken during the next (2020) bulge season.

Short-timescale microlensing events are likely to be produced by substellar brown dwarfs (BDs), but it is difficult to securely identify BD lenses based on only event timescales $t_{\rm E}$ because short-timescale events can also be produced by stellar lenses with high relative lens-source proper motions. In this paper, we report three strong candidate BD-lens events found from the search for lensing events not only with short timescales ($t_{\rm E} \lesssim 6~{\rm days}$) but also with very small angular Einstein radii ($\theta_{\rm E}\lesssim 0.05~{\rm mas}$) among the events that have been found in the 2016--2019 observing seasons. These events include MOA-2017-BLG-147, MOA-2017-BLG-241, and MOA-2019-BLG-256, in which the first two events are produced by single lenses and the last event is produced by a binary lens. From the Bayesian analysis conducted with the combined $t_{\rm E}$ and $\theta_{\rm E}$ constraint, it is estimated that the lens masses of the individual events are $0.051^{+0.100}_{-0.027}~M_\odot$, $0.044^{+0.090}_{-0.023}~M_\odot$, and $0.046^{+0.067}_{-0.023}~M_\odot/0.038^{+0.056}_{-0.019}~M_\odot$ and the probability of the lens mass smaller than the lower limit of stars is $\sim 80\%$ for all events. We point out that routine lens mass measurements of short time-scale lensing events require survey-mode space-based observations.

We report the discovery of a very young high-mass X-ray binary (HMXB) system associated with the supernova remnant (SNR) MCSNRJ0513-6724 in the Large Magellanic Cloud (LMC), using XMM-Newton X-ray observations. The HMXB is located at the geometrical centre of extended soft X-ray emission, which we confirm as an SNR. The HMXB spectrum is consistent with an absorbed power law with spectral index ~1.6 and a luminosity of 7x10^{33} ergs/s (0.2--12 keV). Tentative X-ray pulsations are observed with a periodicity of 4.4 s and the OGLE I-band light curve of the optical counterpart from more than 17.5 years reveals a period of 2.2324\pm0.0003 d, which we interpret as the orbital period of the binary system. The X-ray spectrum of the SNR is consistent with non-equilibrium shock models as expected for young/less evolved SNRs. From the derived ionisation time scale we estimate the age of the SNR to be <6 kyr. The association of the HMXB with the SNR makes it the youngest HMXB, in the earliest evolutionary stage known to date. A HMXB as young as this can switch on as an accreting pulsar only when the spin period has reached a critical value. Under this assumption, we obtain an upper limit to the magnetic field of < 5x10^{11} G. This implies several interesting possibilities including magnetic field burial, possibly by an episode of post-supernova hyper-critical accretion. Since these fields are expected to diffuse out on a timescale of 10^{3}-10^{4} years, the discovery of a very young HMXB can provide us the unique opportunity to observe the evolution of the observable magnetic field for the first time in X-ray binaries.

We report the discovery of a planet in a binary that was discovered from the analysis of the microlensing event OGLE-2018-BLG-1700. We identify the triple nature of the lens from the fact that the complex anomaly pattern can be decomposed into two parts produced by two binary-lens events, in which one binary pair has a very low mass ratio of $\sim 0.01$ between the lens components and the other pair has a mass ratio of $\sim 0.3$. We find two sets of degenerate solutions, in which one solution has a projected separation between the primary and its stellar companion less than the angular Einstein radius $\thetae$ (close solution), while the other solution has a separation greater than $\thetae$ (wide solution). From the Bayesian analysis with the constraints of the event time scale and angular Einstein radius together with the location of the source lying in the far disk behind the bulge, we find that the planet is a super-Jupiter with a mass of $4.4^{+3.0}_{-2.0}~M_{\rm J}$ and the stellar binary components are early and late M-type dwarfs with masses $0.42^{+0.29}_{-0.19}~M_\odot$ and $0.12^{+0.08}_{-0.05}~M_\odot$, respectively, and the planetary system is located at a distance of $D_{\rm L}=7.6^{+1.2}_{-0.9}~{\rm kpc}$. The planet is a circumstellar planet according to the wide solution, while it is a circumbinary planet according to the close solution. The projected primary-planet separation is $2.8^{+3.2}_{-2.5}~{\rm au}$ commonly for the close and wide solutions, but the primary-secondary binary separation of the close solution, $0.75^{+0.87}_{-0.66}~{\rm au}$, is widely different from the separation, $10.5^{+12.1}_{-9.2}~{\rm au}$, of the wide solution.

We present here a new major part of the OGLE Collection of Variable Stars - OGLE Collection of Galactic Cepheids. The new dataset was extracted from the Galaxy Variability Survey images - a dedicated large-scale survey of the Galactic disk and outer bulge conducted by the OGLE project since 2013. The OGLE collection contains 2721 Cepheids of all types - classical, type II and anomalous. It more than doubles the number of known Galactic classical Cepheids. Due to the long-term monitoring and large number of epochs the selected sample is very pure, generally free from contaminating stars of other types often mimicking Cepheids. Its completeness is high at 90% level for classical Cepheids - tested using recent samples of Galactic Cepheids: ASAS-SN, ATLAS, Gaia DR2 and Wise catalog of variable stars. Our comparisons indicate that the completeness of the two latter datasets, Gaia DR2 and Wise catalog, is very low, at < 10% level in the magnitude range of the OGLE GVS survey (10.8 < I < 19.5 mag). Both these samples are severely contaminated by non-Cepheids (the purity is 67% and 56%, respectively). We also present several interesting objects found in the new OGLE Collection - multi-mode pulsators, a first Galactic candidate for eclipsing system containing Cepheid, a binary Cepheid candidate.

We present the final release of the OGLE collection of classical pulsators (Cepheids and RR Lyr stars) in the Large and Small Magellanic Clouds. The sky coverage has been increased from 670 to 765 square degrees compared to the previous edition of the collection. We also add some Cepheids and RR Lyr stars found by the Gaia team and reclassify three Cepheids. Ultimately, our collection consists of 9650 classical Cepheids, 343 type II Cepheids, 278 anomalous Cepheids, and 47828 RR Lyr stars inside and toward the Magellanic System.

We report the discovery of the microlensing planet OGLE-2018-BLG-0740Lb. The planet is detected with a very strong signal of $\Delta\chi^2\sim 4630$, but the interpretation of the signal suffers from two types of degeneracies. One type is caused by the previously known close/wide degeneracy, and the other is caused by an ambiguity between two solutions, in which one solution requires to incorporate finite-source effects, while the other solution is consistent with a point-source interpretation. Although difficult to be firmly resolved based on only the photometric data, the degeneracy is resolved in strong favor of the point-source solution with the additional external information obtained from astrometric and spectroscopic observations. The small astrometric offset between the source and baseline object supports that the blend is the lens and this interpretation is further secured by the consistency of the spectroscopic distance estimate of the blend with the lensing parameters of the point-source solution. The estimated mass of the host is $1.0\pm 0.1~M_\odot$ and the mass of the planet is $4.5\pm 0.6~M_{\rm J}$ (close solution) or $4.8\pm 0.6~M_{\rm J}$ (wide solution) and the lens is located at a distance of $3.2\pm 0.5$~kpc. The bright nature of the lens, with $I\sim 17.1$ ($V\sim 18.2$), combined with its dominance of the observed flux suggest that radial-velocity (RV) follow-up observations of the lens can be done using high-resolution spectrometers mounted on large telescopes, e.g., VLT/ESPRESSO, and this can potentially not only measure the period and eccentricity of the planet but also probe for close-in planets. We estimate that the expected RV amplitude would be $\sim 60\sin i ~{\rm m~s}^{-1}$.

We report the discovery of a $Spitzer$ microlensing planet OGLE-2018-BLG-0596Lb, with preferred planet-host mass ratio $q \sim 2\times10^{-4}$. The planetary signal, which is characterized by a short $(\sim 1~{\rm day})$ "bump" on the rising side of the lensing light curve, was densely covered by ground-based surveys. We find that the signal can be explained by a bright source that fully envelops the planetary caustic, i.e., a "Hollywood" geometry. Combined with the source proper motion measured from $Gaia$, the $Spitzer$ satellite parallax measurement makes it possible to precisely constrain the lens physical parameters. The preferred solution, in which the planet perturbs the minor image due to lensing by the host, yields a Uranus-mass planet with a mass of $M_{\rm p} = 13.9\pm1.6~M_{\oplus}$ orbiting a mid M-dwarf with a mass of $M_{\rm h} = 0.23\pm0.03~M_{\odot}$. There is also a second possible solution that is substantially disfavored but cannot be ruled out, for which the planet perturbs the major image. The latter solution yields $M_{\rm p} = 1.2\pm0.2~M_{\oplus}$ and $M_{\rm h} = 0.15\pm0.02~M_{\odot}$. By combining the microlensing and $Gaia$ data together with a Galactic model, we find in either case that the lens lies on the near side of the Galactic bulge at a distance $D_{\rm L} \sim 6\pm1~{\rm kpc}$. Future adaptive optics observations may decisively resolve the major image/minor image degeneracy.

We report the discovery of the planet OGLE-2018-BLG-0532Lb, with very obvious signatures in the light curve that lead to an estimate of the planet-host mass ratio $q=M_{\rm planet}/M_{\rm host}\simeq 1\times10^{-4}$. Although there are no obvious systematic residuals to this double-lens/single-source (2L1S) fit, we find that $\chi^2$ can be significantly improved by adding either a third lens (3L1S, $\Delta\chi^2=81$) or second source (2L2S, $\Delta\chi^2=65$) to the lens-source geometry. After thorough investigation, we conclude that we cannot decisively distinguish between these two scenarios and therefore focus on the robustly-detected planet. However, given the possible presence of a second planet, we investigate to what degree and with what probability such additional planets may affect seemingly single-planet light curves. Our best estimates for the properties of the lens star and the secure planet are: a host mass $M\sim 0.25\,M_\odot$, system distance $D_L\sim 1\,$kpc and planet mass $m_{p,1}= 8\,M_\oplus$ with projected separation $a_{1,\perp}=1.4\,$au. However, there is a relatively bright $I=18.6$ (and also relatively blue) star projected within $<50\,$mas of the lens, and if future high-resolution images show that this is coincident with the lens, then it is possible that it is the lens, in which case, the lens would be both more massive and more distant than the best-estimated values above.

We present a detailed analysis of Magellanic Bridge Cepheid sample constructed using the OGLE Collection of Variable Stars. Our updated Bridge sample contains 10 classical and 13 anomalous Cepheids. We calculate their individual distances using optical period--Wesenheit relations and construct three-dimensional maps. Classical Cepheids on-sky locations match very well neutral hydrogen and young stars distributions, thus they add to the overall Bridge young population. In three dimensions, eight out of ten classical Cepheids form a bridge-like connection between the Magellanic Clouds. The other two are located slightly farther and may constitute the Counter Bridge. We estimate ages of our Cepheids to be less than 300 Myr for five up to eight out of ten, depending on whether the rotation is included. This is in agreement with a scenario where these stars were formed in-situ after the last encounter of the Magellanic Clouds. Cepheids' proper motions reveal that they are moving away from both Large and Small Magellanic Cloud. Anomalous Cepheids are more spread than classical Cepheids in both two and three dimensions. Even though, they form a rather smooth connection between the Clouds. However, this connection does not seem to be bridge-like, as there are many outliers around both Magellanic Clouds.

We use the extended and updated OGLE Collection of Variable Stars to thoroughly analyze distribution of RR Lyrae stars in the Magellanic Bridge. We use photometric metallicities to derive absolute Wesenheit magnitude and then individual distance of each RR Lyr star. We confirm results from our earlier study showing that RR Lyr stars are present in between the Magellanic Clouds, though their three-dimensional distribution rather resembles two extended overlapping structures than a strict bridge-like connection. The contours do connect in the southern parts of the Bridge, albeit on a too low level to state that there exists an evident connection. To test the sample numerically, we use multi-Gaussian fitting and conclude that there is no additional population or overdensity located in the Bridge. We also try to reproduce results on putative RR Lyr Magellanic Bridge stream by selecting RR Lyr candidates from Gaia Data Release 1. We show that we are not able to obtain the evident connection of the Clouds without many spurious sources in the sample, as the cuts are not able to remove artifacts and not eliminate the evident connection at the same time. Moreover, for the first time we present the Gaia Data Release 2 RR Lyr stars in the Magellanic Bridge area and show that their distribution matches our results.

In this work, we present the analysis of the binary microlensing event OGLE-2018-BLG-0022 that is detected toward the Galactic bulge field. The dense and continuous coverage with the high-quality photometry data from ground-based observations combined with the space-based {\it Spitzer} observations of this long time-scale event enables us to uniquely determine the masses $M_1=0.40 \pm 0.05~M_\odot$ and $M_2=0.13\pm 0.01~M_\odot$ of the individual lens components. Because the lens-source relative parallax and the vector lens-source relative proper motion are unambiguously determined, we can likewise unambiguously predict the astrometric offset between the light centroid of the magnified images (as observed by the {\it Gaia} satellite) and the true position of the source. This prediction can be tested when the individual-epoch {\it Gaia} astrometric measurements are released.