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.

The subtype of hot algol semidetached binaries dubbed Double Periodic Variables (DPVs) are characterized by a photometric cycle longer than the orbital one, whose nature has been related to a magnetic dynamo in the donor component controlling the mass transfer rate. We aim to understand the morphologic changes observed in the light curve of OGLE-BLG-ECL-157529 that are linked to the long cycle. In particular, we want to explain the changes in relative depth of primary and secondary eclipses. We analyze $I$ and $V$-band OGLE photometric times series spanning 18.5 years and model the orbital light curve. We find that OGLE-BLG-ECL-157529 is a new eclipsing Galactic DPV of orbital period 24\fd8, and that its long cycle length decreases in amplitude and length during the time baseline. We show that the changes of the orbital light curve can be reproduced considering an accretion disk of variable thickness and radius, surrounding the hottest stellar component. Our models indicate changes in the temperatures of hot spot and bright spot during the long cycle, and also in the position of the bright spot. This, along with the changes in disk radius might indicate a variable mass transfer in this 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 report the discovery and analysis of the planetary microlensing event OGLE-2017-BLG-0406, which was observed both from the ground and by the ${\it Spitzer}$ satellite in a solar orbit. At high magnification, the anomaly in the light curve was densely observed by ground-based-survey and follow-up groups, and it was found to be explained by a planetary lens with a planet/host mass ratio of $q=7.0 \times 10^{-4}$ from the light-curve modeling. The ground-only and ${\it Spitzer}$-"only" data each provide very strong one-dimensional (1-D) constraints on the 2-D microlens parallax vector $\bf{\pi_{\rm E}}$. When combined, these yield a precise measurement of $\bf{\pi_{\rm E}}$, and so of the masses of the host $M_{\rm host}=0.56\pm0.07\,M_\odot$ and planet $M_{\rm planet} = 0.41 \pm 0.05\,M_{\rm Jup}$. The system lies at a distance $D_{\rm L}=5.2 \pm 0.5 \ {\rm kpc}$ from the Sun toward the Galactic bulge, and the host is more likely to be a disk population star according to the kinematics of the lens. The projected separation of the planet from the host is $a_{\perp} = 3.5 \pm 0.3 \ {\rm au}$, i.e., just over twice the snow line. The Galactic-disk kinematics are established in part from a precise measurement of the source proper motion based on OGLE-IV data. By contrast, the ${\it Gaia}$ proper-motion measurement of the source suffers from a catastrophic $10\,\sigma$ error.

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.

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

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.

The number and properties of observed gravitational microlensing events depend on the distribution and kinematics of stars and other compact objects along the line of sight. In particular, precise measurements of the microlensing optical depth and event rate toward the Galactic bulge enable strict tests of competing models of the Milky Way. Previous estimates, based on samples of up to a few hundred events, gave larger values than expected from the Galactic models and were difficult to reconcile with other constraints on the Galactic structure. Here we used long-term photometric observations of the Galactic bulge by the Optical Gravitational Lensing Experiment (OGLE) to select a homogeneous sample of 8000 gravitational microlensing events. We created the largest and most accurate microlensing optical depth and event rate maps of the Galactic bulge. The new maps ease the tension between the previous measurements and Galactic models. They are consistent with some earlier calculations based on bright stars and are systematically ~30% smaller than the other estimates based on "all-source" samples of microlensing events. The difference is caused by the careful estimation of the source star population. The new maps agree well with predictions based on the Besancon model of the Galaxy. Apart from testing the Milky Way models, our maps may have numerous other applications, such as the measurement of the initial mass function or constraining the dark matter content in the Milky Way center. The new maps will also inform the planning of future space-based microlensing experiments by revising the expected number of events.

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.

The Milky Way is a barred spiral galaxy, with physical properties inferred from various tracers informed by the extrapolation of structures seen in other galaxies. However, the distances of these tracers are measured indirectly and are model-dependent. We constructed a map of the Milky Way in three-dimensions based on the positions and distances of thousands of classical Cepheid variable stars. This map shows the structure of our Galaxy's young stellar population, and allows us to constrain the warped shape of the Milky Way's disk. A simple model of star formation in the spiral arms reproduces the observed distribution of Cepheids.

We present the analysis of the event OGLE-2017-BLG-1186 from the 2017 Spitzer microlensing campaign. This is a remarkable microlensing event because its source is photometrically bright and variable, which makes it possible to perform an asteroseismic analysis using ground-based data. We find that the source star is an oscillating red giant with average timescale of $\sim 9$ d. The asteroseismic analysis also provides us source properties including the source angular size ($\sim 27~\mu{\rm as}$) and distance ($\sim 11.5$ kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian Processes (GPs) in handling the correlated noise caused by the variable source. We find that the parameters from the GP model are generally more loosely constrained than those from the traditional $\chi^2$ minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a $\sim 0.045~M_{\odot}$ brown dwarf at distance $\sim 9.0$ kpc, or a $\sim 0.073~M_{\odot}$ ultracool dwarf at distance $\sim 9.8$ kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a $\sim 35$ per cent probability for the lens to be a bulge object and $\sim 65$ per cent to be a background disc object.

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 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 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 a multiplanetary system found from the analysis of microlensing event OGLE-2018-BLG-1011, for which the light curve exhibits a double-bump anomaly around the peak. We find that the anomaly cannot be fully explained by the binary-lens or binary-source interpretations and its description requires the introduction of an additional lens component. The 3L1S (3 lens components and a single source) modeling yields three sets of solutions, in which one set of solutions indicates that the lens is a planetary system in a binary, while the other two sets imply that the lens is a multiplanetary system. By investigating the fits of the individual models to the detailed light curve structure, we find that the multiple-planet solution with planet-to-host mass ratios $\sim 9.5\times 10^{-3}$ and $\sim 15\times 10^{-3}$ are favored over the other solutions. From the Bayesian analysis, we find that the lens is composed of two planets with masses $1.8^{+3..4}_{-1.1}~M_{\rm J}$ and $2.8^{+5.1}_{-1.7}~M_{\rm J}$ around a host with a mass $0.18^{+0.33}_{-0.10}~M_\odot$ and located at a distance $7.1^{+1.1}_{-1.5}~{\rm kpc}$. The estimated distance indicates that the lens is the farthest system among the known multiplanetary systems. The projected planet-host separations are $a_{\perp,2}=1.8^{+2.1}_{-1.5}~{\rm au}$ ($0.8^{+0.9}_{-0.6}~{\rm au}$) and $a_{\perp,3}=0.8^{+0.9}_{-0.6}~{\rm au}$, where the values of $a_{\perp,2}$ in and out the parenthesis are the separations corresponding to the two degenerate solutions, indicating that both planets are located beyond the snow line of the host, as with the other four multiplanetary systems previously found by microlensing.

We present the discovery and statistical analysis of $12\;660$ spotted variable stars toward and inside the Galactic bulge from over two-decade-long Optical Gravitational Lensing Experiment (OGLE) data. We devise a new method of dereddening of individual stars toward the Galactic bulge where strong and highly nonuniform extinction is present. In effect, $11\;812$ stars were classified as giants and $848$ as dwarfs. Well defined correlations between the luminosity, variability amplitude and rotation period were found for the giants. Rapidly rotating dwarfs with periods $P \leq 2$ d show I-band amplitudes lower than 0.2 mag which is substantially less than the amplitudes of up to 0.8 mag observed in giants and slowly rotating dwarfs. We also notice that amplitudes of stars brighter than $I_0 \approx 16$ mag do not exceed 0.3-0.4 mag. We divide the stars into three groups characterized by correlation between light and color variations. The positive correlation is characteristic for stars that are cooler when fainter, which results from the variable coverage of the stellar surface with spots similar to the sunspots. The variability of stars that are cooler when brighter (negative correlation) can be characterized by chemical spots with overabundance of heavy elements inside and variable line-blanketing effect, which is observed in chemically peculiar stars. The null correlation may results from very high level of the magnetic activity with rapidly variable magnetic fields. This division is readily visible on the color-magnitude diagram (CMD), which suggests that it may depend on the radius of the stars. We detect 79 flaring objects and discuss briefly their properties. Among others, we find that relative brightening during flares is correlated with brightness amplitude.

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 report the mass and distance measurements of two single-lens events from the 2017 Spitzer microlensing campaign. The ground-based observations yield the detection of finite-source effects, and the microlens parallaxes are derived from the joint analysis of ground-based observations and Spitzer observations. We find that the lens of OGLE-2017-BLG-1254 is a $0.60 \pm 0.03 M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.11~\text{kpc}$, where $D_{\rm LS}$ is the distance between the lens and the source. The second event, OGLE-2017-BLG-1161, is subject to the known satellite parallax degeneracy, and thus is either a $0.51^{+0.12}_{-0.10} M_{\odot}$ star with $D_{\rm LS} = 0.40 \pm 0.12~\text{kpc}$ or a $0.38^{+0.13}_{-0.12} M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.19~\text{kpc}$. Both of the lenses are therefore isolated stars in the Galactic bulge. By comparing the mass and distance distributions of the eight published Spitzer finite-source events with the expectations from a Galactic model, we find that the Spitzer sample is in agreement with the probability of finite-source effects occurrence in single lens events.

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.

Planet formation theories predict the existence of free-floating planets that have been ejected from their parent systems. Although they emit little or no light, they can be detected during gravitational microlensing events. Microlensing events caused by rogue planets are characterized by very short timescales $t_{\rm E}$ (typically below two days) and small angular Einstein radii $\theta_{\rm E}$ (up to several uas). Here we present the discovery and characterization of two ultra-short microlensing events identified in data from the Optical Gravitational Lensing Experiment (OGLE) survey, which may have been caused by free-floating or wide-orbit planets. OGLE-2012-BLG-1323 is one of the shortest events discovered thus far ($t_{\rm E}$=0.155 +/- 0.005 d, $\theta_{\rm E}$=2.37 +/- 0.10 uas) and was caused by an Earth-mass object in the Galactic disk or a Neptune-mass planet in the Galactic bulge. OGLE-2017-BLG-0560 ($t_{\rm E}$=0.905 +/- 0.005 d, $\theta_{\rm E}$=38.7 +/- 1.6 uas) was caused by a Jupiter-mass planet in the Galactic disk or a brown dwarf in the bulge. We rule out stellar companions up to a distance of 6.0 and 3.9 au, respectively. We suggest that the lensing objects, whether located on very wide orbits or free-floating, may originate from the same physical mechanism. Although the sample of ultrashort microlensing events is small, these detections are consistent with low-mass wide-orbit or unbound planets being more common than stars in the Milky Way.

The kinematics of isolated brown dwarfs in the Galaxy, beyond the solar neighborhood, is virtually unknown. Microlensing has the potential to probe this hidden population, as it can measure both the mass and five of the six phase-space coordinates (all except the radial velocity) even of a dark isolated lens. However, the measurements of both the microlens parallax and finite-source effects are needed in order to recover the full information. Here, we combine $Spitzer$ satellite parallax measurement with the ground-based light curve, which exhibits strong finite-source effects, of event OGLE-2017-BLG-0896. We find two degenerate solutions for the lens (due to the known satellite-parallax degeneracy), which are consistent with each other except for their proper motion. The lens is an isolated brown dwarf with a mass of either $18\pm1M_J$ or $20\pm1M_J$. This is the lowest isolated-object mass measurement to date, only $\sim$45\% more massive than the theoretical deuterium-fusion boundary at solar metallicity, which is the common definition of a free-floating planet. The brown dwarf is located at either $3.9\pm0.1$ kpc or $4.1\pm0.1$ kpc toward the Galactic bulge, but with proper motion in the opposite direction of disk stars, with one solution suggesting it is moving within the Galactic plane. While it is possibly a halo brown dwarf, it might also represent a different, unknown population.

The nature of type II Cepheids and anomalous Cepheids is still not well known and their evolutionary channels leave many unanswered questions. We use complete collection of classical pulsating stars in the Magellanic Clouds discovered by the OGLE project, to compare their spatial distributions, which are one of the characteristic features directly related to the star formation history. In this analysis we use 9649 classical Cepheids, 262 anomalous Cepheids, 338 type II Cepheids and 46 443 RR Lyr stars from both Magellanic Clouds. We compute three-dimensional Kolmogorov-Smirnov tests for every possible pair of type II and anomalous Cepheids with classical Cepheids, and RR Lyr stars. We confirm that BL Her stars are as old as RR Lyr variable stars - their spatial distributions are similar, and they create a vast halo around both galaxies. We discover that spatial distribution of W Vir stars has attributes characteristic for both young and old stellar populations. Hence, it seems that these similarities are related to the concentration of these stars in the center of the Large Magellanic Cloud, and the lack of a vast halo. This leads to the conclusion that W Vir variables could be a mixture of old and intermediate-age stars. Our analysis of the three-dimensional distributions of anomalous Cepheids shows that they differ significantly from classical Cepheids. Statistical tests of anomalous Cepheids distributions with RR Lyr distributions do not give unambiguous results. We consider that these two distributions can be similar through the vast halos they create. This similarity would confirm anomalous Cepheids evolution scenario that assumes coalescence of a binary system.

We have used deep Optical Gravitational Lensing Experiment (OGLE-IV) images ($V \lesssim 23$ mag, $I \lesssim 23$ mag at $3\sigma$) of the Magellanic System, encompassing an area of $\sim$670 deg$^2$, to perform a search for high-$z$ quasar candidates. We combined the optical OGLE data with the mid-IR Wide-field Infrared Survey Explorer (WISE) 3.4/4.6/12 $\mu$m data, and devised a multi-color selection procedure. We have identified 33 promising sources and then spectroscopically observed the two most variable ones. We report the discovery of two high-$z$ quasars, OGLE J015531-752807 at a redshift $z=5.09$ and OGLE J005907-645016 at a redshift of $z=4.98$. The variability amplitude of both quasars at the rest-frame wavelength $\sim$1300\AA\ is much larger ($\sim$0.4 mag) than other quasars ($<0.15$ mag) at the same rest-frame wavelength but lower redshifts ($2<z<4$). To verify if there exist an increased variability amplitude in high-$z$ population of quasars, simply a larger sample of such sources with a decade long (or longer) light curves is necessary, which will be enabled by the Large Synoptic Survey Telescope (LSST) providing light curves for sources 3-4 mag fainter than OGLE.

We analyze the combined Spitzer and ground-based data for OGLE-2017-BLG-1140 and show that the event was generated by a Jupiter-class $(m_p\simeq 1.6\,M_{\rm jup})$ planet orbiting a mid-late M dwarf $(M\simeq 0.2\,M_\odot)$ that lies $D_{LS}\simeq 1.0\,\mathrm{kpc}$ in the foreground of the microlensed, Galactic-bar, source star. The planet-host projected separation is $a_\perp \simeq 1.0\,\mathrm{au}$, i.e., well-beyond the snow line. By measuring the source proper motion ${\mathbf{\mu}}_s$ from ongoing, long-term OGLE imaging, and combining this with the lens-source relative proper motion ${\mathbf{\mu}}_\mathrm{rel}$ derived from the microlensing solution, we show that the lens proper motion ${\mathbf{\mu}}_l={\mathbf{\mu}}_\mathrm{rel} + {\mathbf{\mu}}_s$ is consistent with the lens lying in the Galactic disk, although a bulge lens is not ruled out. We show that while the Spitzer and ground-based data are comparably well fitted by planetary (i.e., binary-lens, 2L1S) models and by binary-source (1L2S) models, the combination of Spitzer and ground-based data decisively favor the planetary model. This is a new channel to resolve the 2L1S/1L2S degeneracy, which can be difficult to break in some cases.

We report the discovery of microlensing planet OGLE-2017-BLG-0373Lb. We show that while the planet-host system has an unambiguous microlens topology, there are two geometries within this topology that fit the data equally well, which leads to a factor 2.5 difference in planet-host mass ratio, i.e., $q=1.5\times 10^{-3}$ vs. $q=0.6\times 10^{-3}$. We show that this is an "accidental degeneracy" in the sense that it is due to a gap in the data. We dub it "the caustic-chirality degeneracy". We trace the mathematical origins of this degeneracy, which should enable similar degenerate solutions to be easily located in the future. A Bayesian estimate, based on a Galactic model, yields a host mass $M=0.25^{+0.30}_{-0.15} M_\odot$ at a distance $D_L=5.9^{+1.3}_{-1.95}$ kpc. The lens-source relative proper motion is relatively fast, $\mu=9$ mas/yr, which implies that the host mass and distance can be determined by high-resolution imaging after about 10 years. The same observations could in principle resolve the discrete degeneracy in $q$, but this will be more challenging.

We report the discovery of a planetary system in which a super-earth orbits a late M-dwarf host. The planetary system was found from the analysis of the microlensing event OGLE-2017-BLG-0482, wherein the planet signal appears as a short-term anomaly to the smooth lensing light curve produced by the host. Despite its weak signal and short duration, the planetary signal was firmly detected from the dense and continuous coverage by three microlensing surveys. We find a planet/host mass ratio of $q\sim 1.4\times 10^{-4}$. We measure the microlens parallax $\pi_{\rm E}$ from the long-term deviation in the observed lensing light curve, but the angular Einstein radius $\theta_{\rm E}$ cannot be measured because the source trajectory did not cross the planet-induced caustic. Using the measured event timescale and the microlens parallax, we find that the masses of the planet and the host are $M_{\rm p}=9.0_{-4.5}^{+9.0}\ M_\oplus$ and $M_{\rm host}=0.20_{-0.10}^{+0.20}\ M_\odot$, respectively, and the projected separation between them is $a_\perp=1.8_{-0.7}^{+0.6}$ au. The estimated distance to the lens is $D_{\rm L}=5.8_{-2.1}^{+1.8}$ kpc. The discovery of the planetary system demonstrates that microlensing provides an important method to detect low-mass planets orbiting low-mass stars.

Mass measurements of gravitational microlenses require one to determine the microlens parallax $\pie$, but precise $\pie$ measurement, in many cases, is hampered due to the subtlety of the microlens-parallax signal combined with the difficulty of distinguishing the signal from those induced by other higher-order effects. In this work, we present the analysis of the binary-lens event OGLE-2017-BLG-0329, for which $\pie$ is measured with a dramatically improved precision using additional data from space-based $Spitzer$ observations. We find that while the parallax model based on the ground-based data cannot be distinguished from a zero-$\pie$ model at 2$\sigma$ level, the addition of the $Spitzer$ data enables us to identify 2 classes of solutions, each composed of a pair of solutions according to the well-known ecliptic degeneracy. It is found that the space-based data reduce the measurement uncertainties of the north and east components of the microlens-parallax vector $\pivec_{\rm E}$ by factors $\sim 18$ and $\sim 4$, respectively. With the measured microlens parallax combined with the angular Einstein radius measured from the resolved caustic crossings, we find that the lens is composed of a binary with components masses of either $(M_1,M_2)\sim (1.1,0.8)\ M_\odot$ or $\sim (0.4,0.3)\ M_\odot$ according to the two solution classes. The first solution is significantly favored but the second cannot be securely ruled out based on the microlensing data alone. However, the degeneracy can be resolved from adaptive optics observations taken $\sim 10$ years after the event.

We report the discovery of a cold Super-Earth planet (m_p=4.4 +/- 0.5 M_Earth) orbiting a low-mass (M=0.23 +/- 0.03 M_Sun) M dwarf at projected separation a_perp = 1.18 +/- 0.10 AU, i.e., about 1.9 times the snow line. The system is quite nearby for a microlensing planet, D_Lens = 0.86 +/- 0.09 kpc. Indeed, it was the large lens-source relative parallax pi_rel=1.0 mas (combined with the low mass M) that gave rise to the large, and thus well-measured, "microlens parallax" that enabled these precise measurements. OGLE-2017-BLG-1434Lb is the eighth microlensing planet with planet-host mass ratio q < 1 * 10^-4. We apply a new planet-detection sensitivity method, which is a variant of "V/V_max", to seven of these eight planets to derive the mass-ratio function in this regime. We find dN/d(ln q) ~ q^p, with p = 1.05 (+0.78,-0.68), which confirms the "turnover" in the mass function found by Suzuki et al. relative to the power law of opposite sign n = -0.93 +/- 0.13 at higher mass ratios q >~ 2 * 10^-4. We combine our result with that of Suzuki et al. to obtain p = 0.73 (+0.42,-0.34).