The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) $V-$band light curves for sources brighter than V$\lesssim17$ mag across the whole sky. We produced V-band light curves for a total of ${\sim}61.5$ million sources and systematically searched these sources for variability. We identified ${\sim} 426,000$ variables, including ${\sim} 219,000$ new discoveries. Most (${\sim}74\%$) of our discoveries are in the Southern hemisphere. Here we use spectroscopic information from LAMOST, GALAH, RAVE, and APOGEE to study the physical and chemical properties of these variables. We find that metal-poor eclipsing binaries have orbital periods that are shorter than metal-rich systems at fixed temperature. We identified rotational variables on the main-sequence, red giant branch and the red clump. A substantial fraction (${\gtrsim}80\%$) of the rotating giants have large $v_{\rm rot}$ or large NUV excesses also indicative of fast rotation. The rotational variables have unusual abundances suggestive of analysis problems. Semi-regular variables tend to be lower metallicity ($\rm [Fe/H]{\sim}-0.5$) than most giant stars. We find that the APOGEE DR16 temperatures of oxygen-rich semi-regular variables are strongly correlated with the $W_{RP}-W_{JK}$ color index for $\rm T_{eff}\lesssim3800$ K. Using abundance measurements from APOGEE DR16, we find evidence for Mg, O and N enrichment in the semi-regular variables. We find that the Aluminum abundances of the semi-regular variables are strongly correlated with the pulsation period, where the variables with $\rm P\gtrsim 60$ days are significantly depleted in Al.

We characterize ${\sim} 71,200$ W UMa type (EW) contact binaries, including ${\sim} 12,600$ new discoveries, using ASAS-SN $V$-band all-sky light curves along with archival data from Gaia, 2MASS, AllWISE, LAMOST, GALAH, RAVE, and APOGEE. There is a clean break in the EW period-luminosity relation at $\rm \log (\rm P/d){\simeq}-0.30$, separating the longer period early-type (A sub-type) EW binaries from the shorter period, late-type (W sub-type) systems. The two populations are even more cleanly separated in the space of period and effective temperature, by $\rm T_{eff}=6710\,K-1760\,K\,\log(P/0.5\,d)$. Early-type and late-type EW binaries follow opposite trends in $\rm T_{eff}$ with orbital period. For longer periods, early-type EW binaries are cooler, while late-type systems are hotter. We derive period-luminosity relationships (PLRs) in the $W_{JK}$, $V$, Gaia DR2 $G$, $J$, $H$, $K_s$ and $W_1$ bands for the late-type and early-type EW binaries separated both by period and effective temperature, and by period alone. The dichotomy of contact binaries is almost certainly related to the Kraft break and the related changes in envelope structure, winds and angular momentum loss.

The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 ($\mu_{\rm cl}$(N,E) = (+0.36+-0.10, +1.42+-0.10) mas/yr) as well as the source and show that the lens is not a cluster member. Even though this particular event does not probe the distribution of planets in globular clusters, other potential cluster lens events can be verified using our methodology. Additionally, we find that microlens parallax measured using OGLE photometry is consistent with the value found based on the light curve displacement between Earth and Spitzer.

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.

The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) light curves for sources brighter than V$\lesssim17$ mag across the whole sky. As part of our effort to characterize the variability of all the stellar sources visible in ASAS-SN, we have produced ${\sim}30.1$ million V-band light curves for sources in the southern hemisphere using the APASS DR9 catalog as our input source list. We have systematically searched these sources for variability using a pipeline based on random forest classifiers. We have identified ${\sim} 220,000$ variables, including ${\sim} 88,300$ new discoveries. In particular, we have discovered ${\sim}48,000$ red pulsating variables, ${\sim}23,000$ eclipsing binaries, ${\sim}2,200$ $\delta$-Scuti variables and ${\sim}10,200$ rotational variables. The light curves and characteristics of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). The pre-computed ASAS-SN V-band light curves for all the ${\sim}30.1$ million sources are available through the ASAS-SN photometry database (https://asas-sn.osu.edu/photometry). This effort will be extended to provide ASAS-SN light curves for sources in the northern hemisphere and for V$\lesssim17$ mag sources across the whole sky that are not included in APASS DR9.

We explore the synergy between photometric and spectroscopic surveys by searching for periodic variable stars among the targets observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) using photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN). We identified 1924 periodic variables among more than 258000 APOGEE targets; 465 are new discoveries. We homogeneously classified 430 eclipsing and ellipsoidal binaries, 139 classical pulsators (Cepheids, RR Lyrae and delta Scuti), 719 long period variables (pulsating red giants) and 636 rotational variables. The search was performed using both visual inspection and machine learning techniques. The light curves were also modeled with the damped random walk stochastic process. We find that the median [Fe/H] of variable objects is lower by 0.3 dex than that of the overall APOGEE sample. Eclipsing binaries and ellipsoidal variables are shifted to a lower median [Fe/H] by 0.2 dex. Eclipsing binaries and rotational variables exhibit significantly broader spectral lines than the rest of the sample. We make ASAS-SN light curves for all the APOGEE stars publicly available and provide parameters for the variable objects.

We present the analysis of stellar binary microlensing event OGLE-2015-BLG-0060 based on observations obtained from 13 different telescopes. Intensive coverage of the anomalous parts of the light curve was achieved by automated follow-up observations from the robotic telescopes of the Las Cumbres Observatory. We show that, for the first time, all main features of an anomalous microlensing event are well covered by follow-up data, allowing us to estimate the physical parameters of the lens. The strong detection of second-order effects in the event light curve necessitates the inclusion of longer-baseline survey data in order to constrain the parallax vector. We find that the event was most likely caused by a stellar binary-lens with masses $M_{\star1} = 0.87 \pm 0.12 M_{\odot}$ and $M_{\star2} = 0.77 \pm 0.11 M_{\odot}$. The distance to the lensing system is 6.41 $\pm 0.14$ kpc and the projected separation between the two components is 13.85 $\pm 0.16$ AU. Alternative interpretations are also considered.

We analyze the binary microlensing event MOA-2016-BLG-231, which was observed from the ground and from Spitzer. The lens is composed of very low-mass brown dwarfs (BDs) with $M_1 = 21^{+12}_{-5} \ M_J$ and $M_2 = 9^{+5}_{-2}\ M_J$, and it is located in the Galactic disk $D_{\rm L} = 2.85^{+0.88}_{-0.50}\ {\rm kpc}$. This is the fifth binary brown dwarf discovered by microlensing, and the BD binary is moving counter to the orbital motion of disk stars. Constraints on the lens physical properties come from late time, non-caustic-crossing features of the Spitzer light curve. Thus, MOA-2016-BLG-231 shows how Spitzer plays a crucial role in resolving the nature of BDs in binary BD events with short timescale ($\lesssim 10$ days).

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.

We analyze the gravitational binary-lensing event OGLE-2016-BLG-0156, for which the lensing light curve displays pronounced deviations induced by microlens-parallax effects. The light curve exhibits 3 distinctive widely-separated peaks and we find that the multiple-peak feature provides a very tight constraint on the microlens-parallax effect, enabling us to precisely measure the microlens parallax $\pi_{\rm E}$. All the peaks are densely and continuously covered from high-cadence survey observations using globally located telescopes and the analysis of the peaks leads to the precise measurement of the angular Einstein radius $\theta_{\rm E}$. From the combination of the measured $\pi_{\rm E}$ and $\theta_{\rm E}$, we determine the physical parameters of the lens. It is found that the lens is a binary composed of two M dwarfs with masses $M_1=0.18\pm 0.01\ M_\odot$ and $M_2=0.16\pm 0.01\ M_\odot$ located at a distance $D_{\rm L}= 1.35\pm 0.09\ {\rm kpc}$. According to the estimated lens mass and distance, the flux from the lens comprises an important fraction, $\sim 25\%$, of the blended flux. The bright nature of the lens combined with the high relative lens-source motion, $\mu=6.94\pm 0.50\ {\rm mas}\ {\rm yr}^{-1}$, suggests that the lens can be directly observed from future high-resolution follow-up observations.

We report on the discovery and follow-up of a peculiar transient, OGLE17aaj, which occurred in the nucleus of a weakly active galaxy. We investigate whether it can be interpreted as a new candidate for a tidal disruption event (TDE). We present the OGLE-IV light curve that covers the slow 60-day-long rise to maximum along with photometric, spectroscopic, and X-ray follow-up during the first year. OGLE17aaj is a nuclear transient exhibiting some properties similar to previously found TDEs, including a long rise time, lack of colour-temperature evolution, and high black-body temperature. On the other hand, its narrow emission lines and slow post-peak evolution are different from previously observed TDEs. Its spectrum and light-curve evolution is similar to F01004-2237 and AT 2017bgt. Signatures of historical low-level nuclear variability suggest that OGLE17aaj may instead be related to a new type of accretion event in active super-massive black holes.

The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) light curves for sources brighter than V$\lesssim17$ mag across the whole sky. The Transiting Exoplanet Survey Satellite (TESS) has started to produce high-quality light curves with a baseline of at least 27 days, eventually for most of the sky. The combination of ASAS-SN and TESS light curves probes both long and short term variability in great detail, especially towards the TESS continuous viewing zones (CVZ) at the ecliptic poles. We have produced ${\sim}1.3$ million V-band light curves covering a total of ${\sim}1000 \, \rm deg^2$ towards the southern TESS CVZ and have systematically searched these sources for variability. We have identified ${\sim} 11,700$ variables, including ${\sim} 7,000$ new discoveries. The light curves and characteristics of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). We also introduce an online resource to obtain pre-computed ASAS-SN V-band light curves (https://asas-sn.osu.edu/photometry) starting with the light curves of the ${\sim}1.3$ million sources studied in this work. This effort will be extended to provide ASAS-SN light curves for ${\sim}50\;$million sources over the entire sky.

We present an analysis of the microlensing event OGLE-2015-BLG-0232. This event is challenging to characterize for two reasons. First, the light curve is not well sampled during the caustic crossing due to the proximity of the full Moon impacting the photometry quality. Moreover, the source brightness is difficult to estimate because this event is blended with a nearby K dwarf star. We found that the light curve deviations are likely due to a close brown dwarf companion (i.e., s = 0.55 and q = 0.06), but the exact nature of the lens is still unknown. We finally discuss the potential of follow-up observations to estimate the lens mass and distance in the future.

We present the analysis of the caustic-crossing binary microlensing event OGLE-2017-BLG-0039. Thanks to the very long duration of the event, with an event time scale $t_{\rm E}\sim 130$ days, the microlens parallax is precisely measured despite its small value of $\pie\sim 0.06$. The analysis of the well-resolved caustic crossings during both the source star's entrance and exit of the caustic yields the angular Einstein radius $\thetae\sim 0.6$~mas. The measured $\pie$ and $\thetae$ indicate that the lens is a binary composed of two stars with masses $\sim 1.0~M_\odot$ and $\sim 0.15~M_\odot$, and it is located at a distance of $\sim 6$ kpc. From the color and brightness of the lens estimated from the determined lens mass and distance, it is expected that $\sim 2/3$ of the $I$-band blended flux comes from the lens. Therefore, the event is a rare case of a bright lens event for which high-resolution follow-up observations can confirm the nature of the lens.

We report the identification of the bright (V${\sim}13.3$ mag) star FY Sct as a long period detached eclipsing binary using a combined ASAS-SN and ASAS light curve spanning 2000-2018. The orbital period is P${\sim}2.57$ years and the primary eclipse lasts ${\sim}73$ days. The eclipse profile is suggestive of a disk eclipsing binary rather than a stellar component. We also detect ${\sim}0.4$ mag pulsations with a period of P${\sim78}$ d. The next eclipse begins on September, 28, 2018. Further photometric and spectroscopic observations are encouraged, particularly when the system is in eclipse.

The All-Sky Automated Survey for Supernovae (ASAS-SN) is the first optical survey to routinely monitor the whole sky with a cadence of $\sim2-3$ days down to V$\lesssim17$ mag. ASAS-SN has monitored the whole sky since 2014, collecting $\sim100-500$ epochs of observations per field. The V-band light curves for candidate variables identified during the search for supernovae are classified using a random forest classifier and visually verified. We present a catalog of 66,179 bright, new variable stars discovered during our search for supernovae, including 27,479 periodic variables and 38,700 irregular variables. V-band light curves for the ASAS-SN variables are available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). The database will begin to include the light curves of known variable stars in the near future along with the results for a systematic, all-sky variability survey.

The variable stars in the VSX catalog are derived from a multitude of inhomogeneous data sources and classification tools. This inhomogeneity complicates our understanding of variable star types, statistics, and properties, and it directly affects attempts to build training sets for current (and next) generation all-sky, time-domain surveys. We homogeneously analyze the ASAS-SN V-band light curves of ${\sim}412,000$ variables from the VSX catalog. The variables are classified using an updated random forest classifier with an $F_1$ score of 99.4% and refinement criteria for individual classifications. We have derived periods for ${\sim}44,000$ variables in the VSX catalog that lack a period, and have reclassified ${\sim} 17,000$ sources into new broad variability groups with high confidence. We have also reclassified ${\sim} 94,000$ known variables with miscellaneous/generic classifications. The light curves, classifications, and a range of properties of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). We also include the V-band light curves for a set of ${\sim}4,000$ rare variables and transient sources, including cataclysmic variables, symbiotic binaries and flare stars.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

We report the discovery and the analysis of the short timescale binary-lens microlensing event, MOA-2015-BLG-337. The lens system could be a planetary system with a very low mass host, around the brown dwarf/planetary mass boundary, or a brown dwarf binary. We found two competing models that explain the observed light curves with companion/host mass ratios of q~0.01 and ~0.17, respectively. From the measurement of finite source effects in the best-fit planetary model, we find a relatively small angular Einstein radius of theta_E ~ 0.03 mas which favors a low mass lens. We conduct a Bayesian analysis to obtain the probability distribution of the lens properties. The results for the planetary models strongly depend on the minimum mass, M_min, in the assumed mass function. In summary, there are two solutions of the lens system: (1) a brown dwarf/planetary mass boundary object orbited by a super-Neptune (the planetary model with M_min=0.001 M_sun) and (2) a brown dwarf binary (the binary model). If the planetary models is correct, this system can be one of a new class of planetary system, having a low host mass and also a planetary mass ratio (q <0.03) between the host and its companion. The discovery of the event is important for the study of planetary formation in very low mass objects. In addition, it is important to consider all viable solutions in these kinds of ambiguous events in order for the future comprehensive statistical analyses of planetary/binary microlensing events.

The Gaia Data Release 2 (DR2): we summarise the processing and results of the identification of variable source candidates of RR Lyrae stars, Cepheids, long period variables (LPVs), rotation modulation (BY Dra-type) stars, delta Scuti & SX Phoenicis stars, and short-timescale variables. In this release we aim to provide useful but not necessarily complete samples of candidates. The processed Gaia data consist of the G, BP, and RP photometry during the first 22 months of operations as well as positions and parallaxes. Various methods from classical statistics, data mining and time series analysis were applied and tailored to the specific properties of Gaia data, as well as various visualisation tools. The DR2 variability release contains: 228'904 RR Lyrae stars, 11'438 Cepheids, 151'761 LPVs, 147'535 stars with rotation modulation, 8'882 delta Scuti & SX Phoenicis stars, and 3'018 short-timescale variables. These results are distributed over a classification and various Specific Object Studies (SOS) tables in the Gaia archive, along with the three-band time series and associated statistics for the underlying 550'737 unique sources. We estimate that about half of them are newly identified variables. The variability type completeness varies strongly as function of sky position due to the non-uniform sky coverage and intermediate calibration level of this data. The probabilistic and automated nature of this work implies certain completeness and contamination rates which are quantified so that users can anticipate their effects. This means that even well-known variable sources can be missed or misidentified in the published data. The DR2 variability release only represents a small subset of the processed data. Future releases will include more variable sources and data products; however, DR2 shows the (already) very high quality of the data and great promise for variability studies.

We present the analysis of the binary-lens microlensing event OGLE-2017-BLG-0537. The light curve of the event exhibits two strong caustic-crossing spikes among which the second caustic crossing was resolved by high-cadence surveys. It is found that the lens components with a mass ratio $\sim 0.5$ are separated in projection by $\sim 1.3\thetae$, where $\thetae$ is the angular Einstein radius. Analysis of the caustic-crossing part yields $\thetae=1.77\pm 0.16$~mas and a lens-source relative proper motion of $\mu =12.4\pm 1.1~{\rm mas}~{\rm yr}^{-1}$. The measured $\mu$ is the third highest value among the events with measured proper motions and $\sim 3$ times higher than the value of typical Galactic bulge events, making the event a strong candidate for follow-up observations to directly image the lens by separating it from the source. From the angular Einstein radius combined with the microlens parallax, it is estimated that the lens is composed of two main-sequence stars with masses $M_1\sim 0.4~M_\odot$ and $M_2\sim 0.2~M_\odot$ located at a distance of $D_{\rm L}\sim 1.2$~kpc. However, the physical lens parameters are not very secure due to the weak microlens-parallax signal, and thus we cross check the parameters by conducting a Bayesian analysis based on the measured Einstein radius and event timescale combined with the blending constraint. From this, we find that the physical parameters estimated from the Bayesian analysis are consistent with those based on the measured microlens parallax. Resolving the lens from the source can be done in about 5 years from high-resolution follow-up observations and this will provide a rare opportunity to test and refine the microlensing model.

Microlensing is a powerful and unique technique to probe isolated objects in the Galaxy. To study the characteristics of these interesting objects based on the microlensing method, measurement of the microlens parallax is required to determine the properties of the lens. Of the various methods to measure microlens parallax, the most robust way is to make simultaneous ground- and space-based observations, i.e., by measuring the space-based microlens parallax. However, space-based campaigns usually require "expensive" resources. Gould & Yee (2012) proposed an idea called the "cheap space-based microlens parallax" that can measure the lens-parallax using only $2$ or $3$ space-based observations of high-magnification events. This cost-effective observation strategy to measure microlens parallaxes could be used by space-borne telescopes to build a complete sample for studying isolated objects. This would enable a direct measurement of the mass function including both extremely low-mass objects and high-mass stellar remnants. However, to adopt this idea requires a test to check how it would work in actual situations. Thus, we present the first practical test of the idea using the high-magnification microlensing event OGLE-2016-BLG-1045, for which a subset of Spitzer observations fortuitously duplicate the prescription of Gould & Yee (2012). From the test, we confirm that the measurement of the lens-parallax adopting this idea has sufficient accuracy to determine the physical properties of the isolated lens.

We show that dense OGLE and KMTNet $I$-band survey data require four bodies (sources plus lenses) to explain the microlensing light curve of OGLE-2015-BLG-1459. However, these can equally well consist of three lenses and one source (3L1S), two lenses and two sources (2L2S) or one lens and three sources (1L3S). In the 3L1S and 2L2S interpretations, the host is a brown dwarf and the dominant companion is a Neptune-class planet, with the third body (in the 3L1S case) being a Mars-class object that could have been a moon of the planet. In the 1L3S solution, the light curve anomalies are explained by a tight (five stellar radii) low-luminosity binary source that is offset from the principal source of the event by $\sim 0.17\,\au$. These degeneracies are resolved in favor of the 1L3S solution by color effects derived from comparison to MOA data, which are taken in a slightly different ($R/I$) passband. To enable current and future ($WFIRST$) surveys to routinely characterize exomoons and distinguish among such exotic systems requires an observing strategy that includes both a cadence faster than 9 min$^{-1}$ and observations in a second band on a similar timescale.

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.

We report the discovery of a giant planet in the OGLE-2017-BLG-1522 microlensing event. The planetary perturbations were clearly identified by high-cadence survey experiments despite the relatively short event timescale of $t_{\rm E} \sim 7.5$ days. The Einstein radius is unusually small, $\theta_{\rm E} = 0.065\,$mas, implying that the lens system either has very low mass or lies much closer to the microlensed source than the Sun, or both. A Bayesian analysis yields component masses $(M_{\rm host}, M_{\rm planet})=(46_{-25}^{+79}, 0.75_{-0.40}^{+1.26})~M_{\rm J}$ and source-lens distance $D_{\rm LS} = 0.99_{-0.54}^{+0.91}~{\rm kpc}$, implying that this is a brown-dwarf/Jupiter system that probably lies in the Galactic bulge, a location that is also consistent with the relatively low lens-source relative proper motion $\mu = 3.2 \pm 0.5~{\rm mas}~{\rm yr^{-1}}$. The projected companion-host separation is $0.59_{-0.11}^{+0.12}~{\rm AU}$, indicating that the planet is placed beyond the snow line of the host, i.e., $a_{sl} \sim 0.12~{\rm AU}$. Planet formation scenarios combined with the small companion-host mass ratio $q \sim 0.016$ and separation suggest that the companion could be the first discovery of a giant planet that formed in a protoplanetary disk around a brown dwarf host.

We analyze the binary gravitational microlensing event OGLE-2017-BLG-1130 (mass ratio q~0.45), the first published case in which the binary anomaly was only detected by the Spitzer Space Telescope. This event provides strong evidence that some binary signals can be missed by observations from the ground alone but detected by Spitzer. We therefore invert the normal procedure, first finding the lens parameters by fitting the space-based data and then measuring the microlensing parallax using ground-based observations. We also show that the normal four-fold space-based degeneracy in the single-lens case can become a weak eight-fold degeneracy in binary-lens events. Although this degeneracy is resolved in event OGLE-2017-BLG-1130, it might persist in other events.

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, via the microlensing method, of a new very-low-mass binary system. By combining measurements from Earth and from the Spitzer telescope in Earth-trailing orbit, we are able to measure the microlensing parallax of the event, and find that the lens likely consists of an $(12.0 \pm 0.6) M_{\rm J}$ + $(15.7 \pm 1.5) M_{\rm J}$ super-Jupiter / brown-dwarf pair. The binary is located at a distance of $(3.08 \pm 0.18)$ kpc in the Galactic Plane, and the components have a projected separation of $(0.43 \pm 0.03)$ AU. Two alternative solutions with much lower likelihoods are also discussed, an 8- and 6-$M_{\rm J}$ model and a 90- and 70-$M_{\rm J}$ model. Although disfavored at the 3-$\sigma$ and 5-$\sigma$ levels, these alternatives cannot be rejected entirely. We show how the more-massive of these models could be tested with future direct imaging.

Current microlensing surveys are sensitive to free-floating planets down to Earth-mass objects. All published microlensing events attributed to unbound planets were identified based on their short timescale (below two days), but lacked an angular Einstein radius measurement (and hence lacked a significant constraint on the lens mass). Here, we present the discovery of a Neptune-mass free-floating planet candidate in the ultrashort ($t_{\rm E}=0.320\pm0.003$ days) microlensing event OGLE-2016-BLG-1540. The event exhibited strong finite-source effects, which allowed us to measure its angular Einstein radius of $\theta_{\rm E}=9.2\pm0.5\,\mu$as. There remains, however, a degeneracy between the lens mass and distance. The combination of the source proper motion and source-lens relative proper motion measurements favors a Neptune-mass lens located in the Galactic disk. However, we cannot rule out that the lens is a Saturn-mass object belonging to the bulge population. We exclude stellar companions up to 15 au.

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).

We report the discovery of a new gravitationally lensed quasar (double) from the Optical Gravitational Lensing Experiment (OGLE) identified inside the $\sim$670 sq. deg area encompassing the Magellanic Clouds. The source was selected as one of $\sim$60 "red W1-W2" mid-IR objects from WISE and having a significant amount of variability in OGLE for both two (or more) nearby sources. This is the first detection of a gravitational lens, where the discovery is made "the other way around", meaning we first measured the time delay between the two lensed quasar images of $-132<t_{\rm AB}<-76$ days (90% CL), with the median $t_{\rm AB}\approx-102$ days (in the observer frame), and where the fainter image B lags image A. The system consists of the two quasar images separated by 1.5" on the sky, with $I\approx20.0$ mag and $I\approx19.6$ mag, respectively, and a lensing galaxy that becomes detectable as $I \approx 21.5$ mag source, 1.0" from image A, after subtracting the two lensed images. Both quasar images show clear AGN broad emission lines at $z=2.16$ in the NTT spectra. The SED fitting with the fixed source redshift provided the estimate of the lensing galaxy redshift of $z \approx 0.9 \pm 0.2$ (90% CL), while its type is more likely to be elliptical (the SED-inferred and lens-model stellar mass is more likely present in ellipticals) than spiral (preferred redshift by the lens model).

We present a collection of classical, type II, and anomalous Cepheids detected in the OGLE fields toward the Galactic center. The sample contains 87 classical Cepheids pulsating in one, two or three radial modes, 924 type II Cepheids divided into BL Her, W Vir, peculiar W Vir, and RV Tau stars, and 20 anomalous Cepheids - first such objects found in the Galactic bulge. Additionally, we upgrade the OGLE Collection of RR Lyr stars in the Galactic bulge by adding 828 newly identified variables. For all Cepheids and RR Lyr stars, we publish time-series VI photometry obtained during the OGLE-IV project, from 2010 through 2017. We discuss basic properties of our classical pulsators: their spatial distribution, light curve morphology, period-luminosity relations, and position in the Petersen diagram. We present the most interesting individual objects in our collection: a type II Cepheid with additional eclipsing modulation, W Vir stars with the period doubling effect and the RVb phenomenon, a mode-switching RR Lyr star, and a triple-mode anomalous RRd star.

We report the discovery of a sub-Jupiter mass planet orbiting beyond the snow line of an M-dwarf most likely in the Galactic disk as part of the joint Spitzer and ground-based monitoring of microlensing planetary anomalies toward the Galactic bulge. The microlensing parameters are strongly constrained by the light curve modeling and in particular by the Spitzer-based measurement of the microlens parallax, $\pi_\mathrm{E}$. However, in contrast to many planetary microlensing events, there are no caustic crossings, so the angular Einstein radius, $\theta_\mathrm{E}$ has only an upper limit based on the light curve modeling alone. Additionally, the analysis leads us to identify 8 degenerate configurations: the four-fold microlensing parallax degeneracy being doubled by a degeneracy in the caustic structure present at the level of the ground-based solutions. To pinpoint the physical parameters, and at the same time to break the parallax degeneracy, we make use of a series of arguments: the $\chi^2$ hierarchy, the Rich argument, and a prior Galactic model. The preferred configuration is for a host at $D_L=3.73_{-0.67}^{+0.66}~\mathrm{kpc}$ with mass $M_\mathrm{L}=0.30_{-0.12}^{+0.15}~\mathrm{M_\odot}$, orbited by a Saturn-like planet with $M_\mathrm{planet}=0.43_{-0.17}^{+0.21}~\mathrm{M_\mathrm{Jup}}$ at projected separation $a_\perp = 1.70_{-0.39}^{+0.38}~\mathrm{au}$, about 2.1 times beyond the system snow line. Therefore, it adds to the growing population of sub-Jupiter planets orbiting near or beyond the snow line of M-dwarfs discovered by microlensing. Based on the rules of the real-time protocol for the selection of events to be followed up with Spitzer, this planet will not enter the sample for measuring the Galactic distribution of planets.

The Magellanic System (MS) encompasses the nearest neighbors of the Milky Way, the Large (LMC) and Small (SMC) Magellanic Clouds, and the Magellanic Bridge (MBR). This system contains a diverse sample of star clusters. Their parameters, such as the spatial distribution, chemical composition and age distribution yield important information about the formation scenario of the whole Magellanic System. Using deep photometric maps compiled in the fourth phase of the Optical Gravitational Lensing Experiment (OGLE-IV) we present the most complete catalog of star clusters in the Magellanic System ever constructed from homogeneous, long time-scale photometric data. In this second paper of the series, we show the collection of star clusters found in the area of about 360 square degrees in the MBR and in the outer regions of the SMC. Our sample contains 198 visually identified star cluster candidates, 75 of which were not listed in any of the previously published catalogs. The new discoveries are mainly young small open clusters or clusters similar to associations.

We report the discovery of OGLE-2016-BLG-1190Lb, which is likely to be the first Spitzer microlensing planet in the Galactic bulge/bar, an assignation that can be confirmed by two epochs of high-resolution imaging of the combined source-lens baseline object. The planet's mass M_p= 13.4+-0.9 M_J places it right at the deuterium burning limit, i.e., the conventional boundary between "planets" and "brown dwarfs". Its existence raises the question of whether such objects are really "planets" (formed within the disks of their hosts) or "failed stars" (low mass objects formed by gas fragmentation). This question may ultimately be addressed by comparing disk and bulge/bar planets, which is a goal of the Spitzer microlens program. The host is a G dwarf M_host = 0.89+-0.07 M_sun and the planet has a semi-major axis a~2.0 AU. We use Kepler K2 Campaign 9 microlensing data to break the lens-mass degeneracy that generically impacts parallax solutions from Earth-Spitzer observations alone, which is the first successful application of this approach. The microlensing data, derived primarily from near-continuous, ultra-dense survey observations from OGLE, MOA, and three KMTNet telescopes, contain more orbital information than for any previous microlensing planet, but not quite enough to accurately specify the full orbit. However, these data do permit the first rigorous test of microlensing orbital-motion measurements, which are typically derived from data taken over <1% of an orbital period.

We present microlensing planet OGLE-2017-BLG-0173Lb, with planet-host mass ratio either $q\simeq 2.5\times 10^{-5}$ or $q\simeq 6.5\times 10^{-5}$, the lowest or among the lowest ever detected. The planetary perturbation is strongly detected, $\Delta\chi^2\sim 10,000$, because it arises from a bright (therefore, large) source passing over and enveloping the planetary caustic: a so-called "Hollywood" event. The factor $\sim 2.5$ offset in $q$ arises because of a previously unrecognized discrete degeneracy between Hollywood events in which the caustic is fully enveloped and those in which only one flank is enveloped, which we dub "Cannae" and "von Schlieffen", respectively. This degeneracy is "accidental" in that it arises from gaps in the data. Nevertheless, the fact that it appears in a $\Delta\chi^2=10,000$ planetary anomaly is striking. We present a simple formalism to estimate the sensitivity of other Hollywood events to planets and show that they can lead to detections close to, but perhaps not quite reaching, the Earth/Sun mass ratio of $3\times 10^{-6}$. This formalism also enables an analytic understanding of the factor $\sim 2.5$ offset in $q$ between the Cannae and von Schlieffen solutions. The Bayesian estimates for the host-mass, system distance, and planet-host projected separation are $M=0.39^{+0.40}_{-0.24}\,M_\odot$, $D_L=4.8^{+1.5}_{-1.8}\,\kpc$, and $a_\perp=3.8\pm 1.6\,\au$. The two estimates of the planet mass are $m_p=3.3^{+3.8}_{-2.1}\,M_\oplus$ and $m_p=8^{+11}_{-6}\,M_\oplus$. The measured lens-source relative proper motion $\mu=6\,\masyr$ will permit imaging of the lens in about 15 years or at first light on adaptive-optics imagers on next-generation telescopes. These will allow to measure the host mass but probably cannot resolve the planet-host mass-ratio degeneracy.

We analyze an ensemble of microlensing events from the 2015 Spitzer microlensing campaign, all of which were densely monitored by ground-based high-cadence survey teams. The simultaneous observations from Spitzer and the ground yield measurements of the microlensing parallax vector $\pi_{\rm E}$, from which compact constraints on the microlens properties are derived, including $\lesssim$25\% uncertainties on the lens mass and distance. With the current sample, we demonstrate that the majority of microlenses are indeed in the mass range of M dwarfs. The planet sensitivities of all 41 events in the sample are calculated, from which we provide constraints on the planet distribution function. In particular, assuming a planet distribution function that is uniform in $\log{q}$, where $q$ is the planet-to-star mass ratio, we find a $95\%$ upper limit on the fraction of stars that host typical microlensing planets of 49\%, which is consistent with previous studies. Based on this planet-free sample, we develop the methodology to statistically study the Galactic distribution of planets using microlensing parallax measurements. Under the assumption that the planet distributions are the same in the bulge as in the disk, we predict that $\sim$1/3 of all planet detections from the microlensing campaigns with Spitzer should be in the bulge. This prediction will be tested with a much larger sample, and deviations from it can be used to constrain the abundance of planets in the bulge relative to the disk.

We present our study of OGLE-2014-SN-073, one of the brightest Type II SN ever discovered, with an unusually broad lightcurve combined with high ejecta velocities. From our hydrodynamical modelling we infer a remarkable ejecta mass of $60^{+42}_{-16}$~M$_\odot$, and a relatively high explosion energy of $12.4^{+13.0}_{-5.9} \times10^{51}$~erg. We show that this object belongs, with a very small number of other hydrogen-rich SNe, to an energy regime that is not explained by standard core-collapse (CC) neutrino-driven explosions. We compare the quantities inferred by the hydrodynamical modelling with the expectations of various exploding scenarios, trying to explain the high energy and luminosity released. We find some qualitative similarities with pair-instabilities SNe, although a prompt injection of energy by a magnetar seems also a viable alternative to explain such extreme event.

We present the result of microlensing event MOA-2016-BLG-290, which received observations from the two-wheel Kepler (K2), Spitzer, as well as ground-based observatories. A joint analysis of data from K2 and the ground leads to two degenerate solutions of the lens mass and distance. This degeneracy is effectively broken once the (partial) Spitzer light curve is included. Altogether, the lens is found to be an extremely low-mass star located in the Galactic bulge. MOA-2016-BLG-290 is the first microlensing event for which we have signals from three well-separated ($\sim1$ AU) locations. It demonstrates the power of two-satellite microlensing experiment in reducing the ambiguity of lens properties, as pointed out independently by S. Refsdal and A. Gould several decades ago.

We present an analysis of microlensing event OGLE-2016-BLG-0693, based on the survey-only microlensing observations by the OGLE and KMTNet groups. In order to analyze the light curve, we consider the effects of parallax, orbital motion, and baseline slope, and also refine the result using a Galactic model prior. From the microlensing analysis, we find that the event is a binary composed of a low-mass brown dwarf 49+-20 M_J companion and a K- or G-dwarf host, which lies at a distance 5.0+-0.6 kpc toward the Galactic bulge. The projected separation between the brown dwarf and its host star is less than 5 AU, and thus it is likely that the brown dwarf companion is located in the brown dwarf desert.

The reduction of the \emph{K2}'s Campaign 9 (K2C9) microlensing data is challenging mostly because of the very crowded field and the unstable pointing of the spacecraft. In this work, we present the first method that can extract microlensing signals from this K2C9 data product. The raw light curves and the astrometric solutions are first derived, using the techniques from Soares-Furtado et al. and Huang et al. for \emph{K2} dense field photometry. We then minimize and remove the systematic effect by performing simultaneous modeling with the microlensing signal. We also derive precise $(K_p-I)$ vs. $(V-I)$ color-color relations that can predict the microlensing source flux in the \emph{Kepler} bandpass. By implementing the color-color relation in the light curve modeling, we show that the microlensing parameters can be better constrained. In the end, we use two example microlensing events, OGLE-2016-BLG-0980 and OGLE-2016-BLG-0940, to test our method.

We report the discovery of a planet-mass companion to the microlens OGLE-2016-BLG-0263L. Unlike most low-mass companions that were detected through perturbations to the smooth and symmetric light curves produced by the primary, the companion was discovered through the channel of a repeating event, in which the companion itself produced its own single-mass light curve after the event produced by the primary had ended. Thanks to the continuous coverage of the second peak by high-cadence surveys, the possibility of the repeating nature due to source binarity is excluded with a $96\%$ confidence level. The mass of the companion estimated by a Bayesian analysis is $M_{\rm p}=4.1_{-2.5}^{+6.5}\ M_{\rm J}$. The projected primary-companion separation is $a_\perp = 6.5^{+1.3}_{-1.9}$ au. The ratio of the separation to the snow-line distance of $a_\perp/a_{\rm sl}\sim 15.4$ corresponds to the region beyond Neptune, the outermost planet of the solar system. We discuss the importance of high-cadence surveys in expanding the range of microlensing detections of low-mass companions and future space-based microlensing surveys.

Inspired by the discovery of the red nova V1309 Sco (Nova Scorpii 2008) and the fact that its progenitor was a binary system with a rapidly decreasing orbital period, we have searched for period changes in OGLE binary stars. We have selected a sample of 22 462 short-period (P_orb < 4d) eclipsing binary stars observed toward the Galactic bulge by the OGLE-III survey in years 2001-2009. This dataset was extended with photometry from OGLE-II (1997-2000) and the first six years of OGLE-IV (2010-2015). For some stars, the data were supplemented with OGLE-I photometry (1992-1995). After close inspection of the whole sample we have found 56 systems with realistic period decrease and 52 systems with realistic period increase. We have also recognized 35 systems with cyclic period variations. The highest negative period change rate of -1.943 x10^-4 d/y has been detected in detached eclipsing binary OGLE-BLG-ECL-139622 with P_orb = 2.817 d, while all other found systems are contact binaries with orbital periods mostly shorter than 1.0 d. For 22 our systems with decreasing orbital period the absolute rate is higher than the value reported recently for eclipsing binary KIC 9832227. Interestingly, there is an excess of systems with high negative period change rate over systems with positive rate. We cannot exclude the possibility that some of the contact binaries with relatively long orbital period and high negative period change rate will merge in the future. However, our results rather point to the presence of tertiary companions in the observed systems and/or spot activity on the surface of the binary components.

We present the analysis of the binary gravitational microlensing event MOA-2015-BLG-020. The event has a fairly long timescale (about 63 days) and thus the light curve deviates significantly from the lensing model that is based on the rectilinear lens-source relative motion. This enables us to measure the microlensing parallax through the annual parallax effect. The microlensing parallax parameters constrained by the ground-based data are confirmed by the Spitzer observations through the satellite parallax method. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined. It is found that the binary lens is composed of two dwarf stars with masses $M_1 = 0.606 \pm 0.028M_\odot$ and $M_2 = 0.125 \pm 0.006M_\odot$ in the Galactic disk. Assuming the source star is at the same distance as the bulge red clump stars, we find the lens is at a distance $D_L = 2.44 \pm 0.10 kpc$. In the end, we provide a summary and short discussion of all published microlensing events in which the annual parallax effect is confirmed by other independent observations.

Gravitational microlensing is the only method capable of exploring the entire population of free-floating planets down to Mars-mass objects, because the microlensing signal does not depend on the brightness of the lensing object. A characteristic timescale of microlensing events depends on the mass of the lens: the less massive the lens, the shorter the microlensing event. A previous analysis of 474 microlensing events found an excess of very short events (1-2 days) - more than known stellar populations would suggest - indicating the existence of a large population of unbound or wide-orbit Jupiter-mass planets (reported to be almost twice as common as main-sequence stars). These results, however, do not match predictions of planet formation theories and are in conflict with surveys of young clusters. Here we report the analysis of a six times larger sample of microlensing events discovered during the years 2010-2015. Although our survey has very high sensitivity (detection efficiency) to short-timescale (1--2 days) microlensing events, we found no excess of events with timescales in this range, with a 95% upper limit on the frequency of Jupiter-mass free-floating or wide-orbit planets of 0.25 planet per main-sequence star. We detected a few possible ultrashort-timescale events (with timescales of less than 0.5 day), which may indicate the existence of Earth- and super-Earth-mass free-floating planets, as predicted by planet-formation theories. [abridged]

We present a detailed elemental abundance study of 90 F and G dwarf, turn-off and subgiant stars in the Galactic bulge. Based on high-resolution spectra acquired during gravitational microlensing events, stellar ages and abundances for 11 elements (Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Zn, Y and Ba) have been determined. We find that the Galactic bulge has a wide metallicity distribution with significant peaks at [Fe/H]=-1.09, -0.63, -0.20, +0.12, +0.41. We also find a high fraction of intermediate-age to young stars: at [Fe/H]>0 more than 35 % are younger than 8 Gyr. For [Fe/H]<-0.5 most stars are 10 Gyr or older. We have also identified several episodes when significant star formation in the bulge happened: 3, 6, 8, and 12 Gyr ago. We further find that the "knee" in the alpha-element abundance trends of the sub-solar metallicity bulge is located at about 0.1 dex higher [Fe/H] than in the local thick disk. The Galactic bulge has complex age and abundance properties that appear to be tightly connected to the main Galactic stellar populations. In particular, the peaks in the metallicity distribution, the star formation episodes, and the abundance trends, show similarities with the properties of the Galactic thin and thick disks. At the same time there are additional components not seen outside the bulge region, and that most likely can be associated with the Galactic bar. For instance, the star formation rate appears to have been slightly faster in the bulge than in the local thick disk, which most likely is an indication of the denser stellar environment closer to the Galactic centre. Our results strengthen the observational evidence that support the idea of a secular origin for the Galactic bulge, formed out of the other main Galactic stellar populations present in the central regions of our Galaxy.

We present a three-dimensional analysis of a sample of 22 859 type $ab$ RR Lyrae stars in the Magellanic System from the OGLE-IV Collection of RR Lyrae stars. The distance to each object was calculated based on its photometric metallicity and a theoretical relation between color, absolute magnitude and metallicity. The LMC RR Lyrae distribution is very regular and does not show any substructures. We demonstrate that the bar found in previous studies may be an overdensity caused by blending and crowding effects. The halo is asymmetrical with a higher stellar density in its north-eastern area, which is also located closer to us. Triaxial ellipsoids were fitted to surfaces of a constant number density. Ellipsoids farther from the LMC center are less elongated and slightly rotated toward the SMC. The inclination and position angle change significantly with the $a$ axis size. The median axis ratio is $1:1.23:1.45$. The RR Lyrae distribution in the SMC has a very regular, ellipsoidal shape and does not show any substructures or asymmetries. All triaxial ellipsoids fitted to surfaces of a constant number density have virtually the same shape (axis ratio) and are elongated along the line of sight. The median axis ratio is $1:1.10:2.13$. The inclination angle is very small and thus the position angle is not well defined. We present the distribution of RR Lyrae stars in the Magellanic Bridge area, showing that the Magellanic Clouds' halos overlap. A comparison of the distributions of RR Lyrae stars and Classical Cepheids shows that the former are significantly more spread and distributed regularly, while the latter are very clumped and form several distinct substructures.

More than a century ago, Henrietta Leavitt discovered the first Cepheids in the Magellanic Clouds together with the famous period-luminosity relationship revealed by these stars, which soon after revolutionized our view of the Universe. Over the years, the number of known Cepheids in these galaxies has steadily increased with the breakthrough in the last two decades thanks to the new generation of large-scale long-term sky variability surveys. Here we present the final upgrade of the OGLE Collection of Cepheids in the Magellanic System which already contained the vast majority of known Cepheids. The updated collection now comprises 9649 classical and 262 anomalous Cepheids. Type-II Cepheids will be updated shortly. Thanks to high completeness of the OGLE survey the sample of classical Cepheids includes virtually all stars of this type in the Magellanic Clouds. Thus, the OGLE survey concludes the work started by Leavitt. Additionally, the OGLE sample of RR Lyrae stars in the Magellanic System has been updated. It now counts 46 443 variables. A collection of seven anomalous Cepheids in the halo of our Galaxy detected in front of the Magellanic Clouds is also presented. OGLE photometric data are available to the astronomical community from the OGLE Internet Archive. The time-series photometry of all pulsating stars in the OGLE Collection has been supplemented with new observations.

Regular intrinsic brightness variations observed in many stars are caused by pulsations. These pulsations provide information on the global and structural parameters of the star. The pulsation periods range from seconds to years, depending on the compactness of the star and properties of the matter that forms its outer layers. Here, we report the discovery of more than a dozen of previously unknown short-period variable stars: blue large-amplitude pulsators. These objects show very regular brightness variations with periods in the range of 20-40 min and amplitudes of 0.2-0.4 mag in the optical passbands. The phased light curves have a characteristic sawtooth shape, similar to the shape of classical Cepheids and RR Lyrae-type stars pulsating in the fundamental mode. The objects are significantly bluer than main sequence stars observed in the same fields, which indicates that all of them are hot stars. Follow-up spectroscopy confirms a high surface temperature of about 30,000 K. Temperature and colour changes over the cycle prove the pulsational nature of the variables. However, large-amplitude pulsations at such short periods are not observed in any known type of stars, including hot objects. Long-term photometric observations show that the variable stars are very stable over time. Derived rates of period change are of the order of 10^-7 per year and, in most cases, they are positive. According to pulsation theory, such large-amplitude oscillations may occur in evolved low-mass stars that have inflated helium-enriched envelopes. The evolutionary path that could lead to such stellar configurations remains unknown.

The microlens parallax is a crucial observable for conclusively identifying the nature of lens systems in microlensing events containing or composed of faint (even dark) astronomical objects such as planets, neutron stars, brown dwarfs, and black holes. With the commencement of a new era of microlensing in collaboration with space-based observations, the microlens parallax can be routinely measured. In addition, space-based observations can provide opportunities to verify the microlens parallax measured from ground-only observations and to find a unique solution of the lensing lightcurve analysis. However, since most space-based observations cannot cover the full lightcurves of lensing events, it is also necessary to verify the reliability of the information extracted from fragmentary space-based lightcurves. We conduct a test based on the microlensing event OGLE-2016-BLG-0168 created by a binary lens system consisting of almost equal mass M-dwarf stars to demonstrate that it is possible to verify the microlens parallax and to resolve degeneracies by using the space-based lightcurve even though the observations are fragmentary. Since space-based observatories will frequently produce fragmentary lightcurves due to their short observing windows, the methodology of this test will be useful for next-generation microlensing experiments that combine space-based and ground-based collaboration.

We report the discovery of a binary composed of two brown dwarfs, based on the analysis of the microlensing event OGLE-2016-BLG-1469. Thanks to detection of both finite-source and microlens-parallax effects, we are able to measure both the masses $M_1\sim 0.05\ M_\odot$, $M_2\sim 0.01\ M_\odot$, and distance $D_{\rm L} \sim 4.5$ kpc, as well as the projected separation $a_\perp \sim 0.33$ au. This is the third brown-dwarf binary detected using the microlensing method, demonstrating the usefulness of microlensing in detecting field brown-dwarf binaries with separations less than 1 au.

We report the analysis of the first resolved caustic-crossing binary-source microlensing event OGLE-2016-BLG-1003. The event is densely covered by the round-the-clock observations of three surveys. The light curve is characterized by two nested caustic-crossing features, which is unusual for typical caustic-crossing perturbations. From the modeling of the light curve, we find that the anomaly is produced by a binary source passing over a caustic formed by a binary lens. The result proves the importance of high-cadence and continuous observations, and the capability of second-generation microlensing experiments to identify such complex perturbations that are previously unknown. However, the result also raises the issues of the limitations of current analysis techniques for understanding lens systems beyond two masses and of determining the appropriate multiband observing strategy of survey experiments.

We report discovery of the lowest mass ratio exoplanet to be found by the microlensing method in the light curve of the event OGLE~2016--BLG--1195. This planet revealed itself as a small deviation from a microlensing single lens profile from an examination of the survey data soon after the planetary signal. The duration of the planetary signal is $\sim 2.5\,$hours. The measured ratio of the planet mass to its host star is $q = 4.2\pm 0.7 \times10^{-5}$. We further estimate that the lens system is likely to comprise a cold $\sim$3 Earth mass planet in a $\sim\,$2 AU wide orbit around a 0.2 Solar mass star at an overall distance of 7.1 kpc.

We analyze the single microlensing event OGLE-2015-BLG-1482 simultaneously observed from two ground-based surveys and from \textit{Spitzer}. The \textit{Spitzer} data exhibit finite-source effects due to the passage of the lens close to or directly over the surface of the source star as seen from \textit{Spitzer}. Such finite-source effects generally yield measurements of the angular Einstein radius, which when combined with the microlens parallax derived from a comparison between the ground-based and the \textit{Spitzer} light curves, yields the lens mass and lens-source relative parallax. From this analysis, we find that the lens of OGLE-2015-BLG-1482 is a very low-mass star with the mass $0.10 \pm 0.02 \ M_\odot$ or a brown dwarf with the mass $55\pm 9 \ M_{J}$, which are respectively located at $D_{\rm LS} = 0.80 \pm 0.19\ \textrm{kpc}$ and $ D_{\rm LS} = 0.54 \pm 0.08\ \textrm{kpc}$, and thus it is the first isolated low-mass microlens that has been decisively located in the Galactic bulge. The fundamental reason for the degeneracy is that the finite-source effect is seen only in a single data point from \textit{Spitzer} and this single data point gives rise to two solutions for $\rho$. Because the $\rho$ degeneracy can be resolved only by relatively high cadence observations around the peak, while the \textit{Spitzer} cadence is typically $\sim 1\,{\rm day}^{-1}$, we expect that events for which the finite-source effect is seen only in the \textit{Spitzer} data may frequently exhibit this $\rho$ degeneracy. For OGLE-2015-BLG-1482, the relative proper motion of the lens and source for the low-mass star is $\mu_{\rm rel} = 9.0 \pm 1.9\ \textrm{mas yr$^{-1}$}$, while for the brown dwarf it is $5.5 \pm 0.5\ \textrm{mas yr$^{-1}$}$. Hence, the degeneracy can be resolved within $\sim 10\ \rm yrs$ from direct lens imaging by using next-generation instruments with high spatial resolution.

We report the discovery of a high mass-ratio planet $q=0.012$, i.e., 13 times higher than the Jupiter/Sun ratio. The host mass is not presently measured but can be determined or strongly constrained from adaptive optics imaging. The planet was discovered in a small archival study of high-magnification events in pure-survey microlensing data, which was unbiased by the presence of anomalies. The fact that it was previously unnoticed may indicate that more such planets lie in archival data and could be discovered by similar systematic study. In order to understand the transition from predominantly survey+followup to predominately survey-only planet detections, we conduct the first analysis of these detections in the observational $(s,q)$ plane. Here $s$ is projected separation in units of the Einstein radius. We find some evidence that survey+followup is relatively more sensitive to planets near the Einstein ring, but that there is no statistical difference in sensitivity by mass ratio.

In the process of analyzing an observed light curve, one often confronts various scenarios that can mimic the planetary signals causing difficulties in the accurate interpretation of the lens system. In this paper, we present the analysis of the microlensing event OGLE-2016-BLG-0733. The light curve of the event shows a long-term asymmetric perturbation that would appear to be due to a planet. From the detailed modeling of the lensing light curve, however, we find that the perturbation originates from the binarity of the source rather than the lens. This result demonstrates that binary sources with roughly equal-luminosity components can mimic long-term perturbations induced by planets with projected separations near the Einstein ring. The result also represents the importance of the consideration of various interpretations in planet-like perturbations and of high-cadence observations for ensuring the unambiguous detection of the planet.

We present a collection of 450 598 eclipsing and ellipsoidal binary systems detected in the OGLE fields toward the Galactic bulge. The collection consists of binary systems of all types: detached, semi-detached, and contact eclipsing binaries, RS CVn stars, cataclysmic variables, HW Vir binaries, double periodic variables, and even planetary transits. For all stars we provide the I- and V-band time-series photometry obtained during the OGLE-II, OGLE-III, and OGLE-IV surveys. We discuss methods used to identify binary systems in the OGLE data and present several objects of particular interest.

We present the collection of eclipsing binaries in the Large and Small Magellanic Clouds, based on the OGLE survey. It contains 48605 systems, 40204 belonging to the LMC and 8401 to the SMC. Out of the total number of the presented binaries, 16374 are the new ones. We present the time-series photometry obtained for the selected objects during the fourth phase of the OGLE project. The catalog has been created using a two step machine learning procedure based on the Random Forest algorithm.

LSST is expected to yield ~10^7 light curves over the course of its mission, which will require a concerted effort in automated classification. Stochastic processes provide one means of quantitatively describing variability with the potential advantage over simple light curve statistics that the parameters may be physically meaningful. Here, we survey a large sample of periodic, quasi-periodic, and stochastic OGLE-III variables using the damped random walk (DRW, CARMA(1,0)) and quasi-periodic oscillation (QPO, CARMA(2,1)) stochastic process models. The QPO model is described by an amplitude, a period, and a coherence time-scale, while the DRW has only an amplitude and a time-scale. We find that the periodic and quasi-periodic stellar variables are generally better described by a QPO than a DRW, while quasars are better described by the DRW model. There are ambiguities in interpreting the QPO coherence time due to non-sinusoidal light curve shapes, signal-to-noise, error mischaracterizations, and cadence. Higher-order implementations of the QPO model that better capture light curve shapes are necessary for the coherence time to have its implied physical meaning. Independent of physical meaning, the extra parameter of the QPO model successfully distinguishes most of the classes of periodic and quasi-periodic variables we consider.

We present a comparison of the Gaia DR1 samples of pulsating variable stars - Cepheids and RR Lyrae type - with the OGLE Collection of Variable Stars aiming at the characterization of the Gaia mission performance in the stellar variability domain. Out of 575 Cepheids and 2322 RR Lyrae candidates from the Gaia DR1 samples located in the OGLE footprint in the sky, 559 Cepheids and 2302 RR Lyrae stars are genuine pulsators of these types. The number of misclassified stars is low indicating reliable performance of the Gaia data pipeline. The completeness of the Gaia DR1 samples of Cepheids and RR Lyrae stars is at the level of 60-75% as compared to the OGLE Collection dataset. This level of completeness is moderate and may limit the applicability of the Gaia data in many projects.

In an era of extensive photometric observations, the catalogs of RR Lyr type variable stars number tens of thousands of objects. The relation between the iron abundance [Fe/H] and the Fourier parameters of the stars light curve allows us to investigate mean metallicities and metallicity gradients in various stellar environments, independently of time-consuming spectroscopic observations. In this paper we use almost 6500 $V$- and $I$-band light curves of fundamental mode RR Lyr stars from the OGLE-IV survey to provide a relation between the $V$- and $I$-band phase parameter $\varphi_{31}$ used to estimate [Fe/H]. The relation depends on metallicity, which limits its applicability. We apply this relation to metallicity formulae developed for the Johnson $V$- and the Kepler $Kp$-band to obtain the relation between [Fe/H] and $\varphi_{31}$ for the $I$-band photometry. Last, we apply the new relation of Nemec to the OGLE-IV fundamental mode RR Lyr stars data and construct a metallicity map of the Magellanic Clouds. Median [Fe/H] is $-1.39\pm0.44$ dex for the LMC and $-1.77\pm0.48$ dex for the SMC, on the Jurcsik metallicity scale. We also find a metallicity gradient within the LMC with a slope of $-0.029\pm0.002$ dex/kpc in the inner 5 kpc and $-0.030 \pm0.003$ dex/kpc beyond 8 kpc, and no gradient in-between ($-0.019\pm0.002$ dex/kpc integrally). We do not observe a metallicity gradient in the SMC, although we show that the metal-rich RRab stars are more concentrated toward the SMC center than the metal-poor.

In this paper, we present the analysis of the binary gravitational microlensing event OGLE-2015-BLG-0196. The event lasted for almost a year and the light curve exhibited significant deviations from the lensing model based on the rectilinear lens-source relative motion, enabling us to measure the microlens parallax. The ground-based microlens parallax is confirmed by the data obtained from space-based microlens observations using the {\it Spitzer} telescope. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined up to the two-fold degeneracy: $u_0<0$ and $u_0>0$ solutions caused by the well-known "ecliptic" degeneracy. It is found that the binary lens is composed of two M dwarf stars with similar masses $M_1=0.38\pm 0.04\ M_\odot$ ($0.50\pm 0.05\ M_\odot)$ and $M_2=0.38\pm 0.04\ M_\odot$ ($0.55\pm 0.06\ M_\odot$) and the distance to the lens is $D_{\rm L}=2.77\pm 0.23$ kpc ($3.30\pm 0.29$ kpc). Here the physical parameters out and in the parenthesis are for the $u_0<0$ and $u_0>0$ solutions, respectively.

We present the discovery and first three months of follow-up observations of a currently on-going unusual transient detected by the OGLE-IV survey, located in the centre of a galaxy at redshift z=0.1655. The long rise to absolute magnitude of -20.5 mag, slow decline, very broad He and H spectral features make OGLE16aaa similar to other optical/UV Tidal Disruption Events (TDEs). Weak narrow emission lines in the spectrum and archival photometric observations suggest the host galaxy is a weak-line Active Galactic Nucleus (AGN), which has been accreting at higher rate in the past. OGLE16aaa, along with SDSS J0748, seems to form a sub-class of TDEs by weakly or recently active super-massive black holes (SMBHs). This class might bridge the TDEs by quiescent SMBHs and flares observed as "changing-look QSOs", if we interpret the latter as TDEs. If this picture is true, the previously applied requirement for identifying a flare as a TDE that it had to come from an inactive nucleus, could be leading to observational bias in TDE selection, thus affecting TDE-rate estimations.

The Magellanic System (MS), consisting of the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC) and the Magellanic Bridge (MBR), contains diverse sample of star clusters. Their spatial distribution, ages and chemical abundances may provide important information about the history of formation of the whole System. We use deep photometric maps derived from the images collected during the fourth phase of The Optical Gravitational Lensing Experiment (OGLE-IV) to construct the most complete catalog of star clusters in the Large Magellanic Cloud using the homogeneous photometric data. In this paper we present the collection of star clusters found in the area of about 225 square degrees in the outer regions of the LMC. Our sample contains 679 visually identified star cluster candidates, 226 of which were not listed in any of the previously published catalogs. The new clusters are mainly young small open clusters or clusters similar to associations.

Cataclysmic variable stars (CVs) are close binary systems consisting of a white dwarf (primary) that is accreting matter from a low-mass companion star (secondary). From time to time such systems undergo large-amplitude brightenings. The most spectacular eruptions, over $10^4$ times in brightness, occur in classical novae and are caused by a thermonuclear runaway on the surface of the white dwarf. Such eruptions are thought to recur on timescales of $10^4-10^6$. In between, the system's properties depend primarily on the mass-transfer rate: if it is lower than a $10^{-9} M_{\odot}$/year, the accretion becomes unstable and the matter is dumped onto the white dwarf during quasi-periodic dwarf nova outbursts. The hibernation hypothesis predicts that nova eruptions strongly affect the mass-transfer rate, keeping it high for centuries after the event. Subsequently, the mass-transfer rate should significantly decrease for $10^3-10^6$ years, starting the hibernation phase. After that the nova awakes again - with accretion returning to the pre-eruption level and leading to a new nova explosion. The hibernation model predicts cyclical evolution of CVs through phases of high and low mass-transfer. The theory gained some support from the discovery of ancient nova shells around dwarf novae Z Cam and AT Cnc, but direct evidence for considerable mass-transfer changes prior, during and after nova eruptions has not hitherto been found. Here we report long-term observations of the classical nova V1213 Cen (Nova Cen 2009) covering its pre- and post-eruption phases. Within the six years before the explosion, the system revealed dwarf nova outbursts indicative of a low mass-transfer rate. The post-nova is two orders of magnitude brighter than the pre-nova at minimum light with no trace of dwarf nova behavior, implying that the mass-transfer rate increased considerably as a result of the nova explosion.

We report the discovery of a new subclass of double-mode RR Lyrae stars in the Large and Small Magellanic Clouds. The sample of 22 pulsating stars have been extracted from the latest edition of the OGLE collection of RR Lyrae variables in the Magellanic System. The stars pulsating simultaneously in the fundamental (F) and first-overtone (1O) modes have distinctly different properties than regular double-mode RR Lyrae variables (RRd stars). The $P_{1O}/P_F$ period ratios of our anomalous RRd stars are within a range 0.725-0.738, while "classical" double-mode RR Lyrae variables have period ratios in the range 0.742-0.748. In contrast to the typical RRd stars, in the majority of the anomalous pulsators the F-mode amplitudes are higher than the 1O-mode amplitudes. The light curves associated with the F-mode in the anomalous RRd stars show different morphology than the light curves of, both, regular RRd stars and single-mode RRab stars. Most of the anomalous double-mode stars show long-term modulations of the amplitudes (Blazhko-like effect). Translating the period ratios into the abundance parameter, Z, we find for our stars Z in the range (0.002,0.005) - an order of magnitude higher values than typical for RR Lyrae stars. The mass range of the RRd stars inferred from the $W_I$ vs. $P_F$ diagram is (0.55-0.75) of the solar mass. These parameters cannot be accounted for with single star evolution assuming a Reimers-like mass loss. Much greater mass loss caused by interaction with other stars is postulated. We blame the peculiar pulsation properties of our stars to the parametric resonance instability of the 1O-mode to excitation of the F- and 2O-modes as with the inferred parameters of the stars $2\omega_{\rm 1O}\approx\omega_{\rm F}+\omega_{\rm 2O}$.

We use the fundamental-mode RR Lyr-type variable stars (RRab) from OGLE-IV survey to draw a 3D picture of the central part of the tidally disrupted Sagittarius Dwarf Spheroidal (Sgr dSph) galaxy. We estimate the line-of-sight thickness of the Sgr dSph stream to be FWHM_cen=2.42 kpc. Based on OGLE-IV observations collected in seasons 2011-2014 we conduct a comprehensive study of stellar variability in the field of the globular cluster M54 (NGC 6715) residing in the core of this dwarf galaxy. Among the total number of 268 detected variable stars we report the identification of 174 RR Lyr stars, four Type II Cepheids, 51 semi-regular variable red giants, three SX Phe-type stars, 18 eclipsing binary systems. Eighty-three variable stars are new discoveries. The distance to the cluster determined from RRab stars is d_M54 = 26.7 +/-0.03(stat) +/-1.3(sys) kpc. From the location of RRab stars in the period-amplitude (Bailey) diagram we confirm the presence of two old populations, both in the cluster and the Sgr dSph stream.

We present the largest collection of RR Lyrae stars in the Magellanic System and in its foreground. The sample consists of 45 451 RR Lyr stars, of which 39 082 were detected toward the Large Magellanic Cloud and 6369 toward the Small Magellanic Cloud. We provide long-term time-series photometric measurements collected during the fourth phase of the Optical Gravitational Lensing Experiment (OGLE-IV). We discuss several potential astrophysical applications of our collection: investigation of the structure of the Magellanic Clouds and the Galactic halo, studies of the globular clusters in the Magellanic System, analysis of double-mode RR Lyr stars, and search for RR Lyr stars in eclipsing binary systems.

We analyzed a sample of 9418 fundamental-mode and first-overtone Classical Cepheids from the OGLE-IV Collection of Classical Cepheids. The distance to each Cepheid was calculated using the period-luminosity relation for the Wesenheit magnitude, fitted to our data. The classical Cepheids in the LMC are situated mainly in the bar and in the northern arm. The eastern part of the LMC is closer to us and the plane fit to the whole LMC sample yields the inclination $i=24.2\pm0.7$ deg and position angle ${\rm P.A.}=151.4\pm1.7$ deg. We redefined the LMC bar by extending it in the western direction and found no offset from the plane of the LMC contrary to previous studies. On the other hand, we found that the northern arm is offset from a plane by about $-0.5$ kpc, which was not observed before. The age distribution of the LMC Cepheids shows one maximum at about 100 Myr. We demonstrate that the SMC has a non-planar structure and can be described as an extended ellipsoid. We identified two large ellipsoidal off-axis structures in the SMC. The northern one is located closer to us and is younger, while the south-western is farther and older. The age distribution of the SMC Cepheids is bimodal with one maximum at 110 Myr, and another one at 220 Myr. Younger stars are located in the closer part of this galaxy while older ones are more distant. We classified nine Cepheids from our sample as Magellanic Bridge objects. These Cepheids show a large spread in three-dimensions although five of them form a connection between the Clouds. The closest one is closer than any of the LMC Cepheids, while the farthest one -- farther than any SMC Cepheid. All but one Cepheids in the Magellanic Bridge are younger than 300 Myr. The oldest one can be associated with the SMC Wing.

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

We report the discovery of an extrasolar planet detected from the combined data of a microlensing event OGLE-2015-BLG-0051/KMT-2015-BLG-0048 acquired by two microlensing surveys. Despite that the short planetary signal occurred in the very early Bulge season during which the lensing event could be seen for just about an hour, the signal was continuously and densely covered. From the Bayesian analysis using models of the mass function, matter and velocity distributions combined with the information of the angular Einstein radius, it is found that the host of the planet is located in the Galactic bulge. The planet has a mass $0.72_{-0.07}^{+0.65}\ M_{\rm J}$ and it is orbiting a low-mass M-dwarf host with a projected separation $d_\perp=0.73 \pm 0.08$ AU. The discovery of the planet demonstrates the capability of the current high-cadence microlensing lensing surveys in detecting and characterizing planets.

We report the characterization of a massive (m_p=3.9 +- 1.4 M_jup) microlensing planet (OGLE-2015-BLG-0954Lb) orbiting an M dwarf host (M=0.33 +- 0.12 M_sun) at a distance toward the Galactic bulge of 0.6 (+0.4,-0.2) kpc, which is extremely nearby by microlensing standards. The planet-host projected separation is a_perp ~ 1.2 AU. The characterization was made possible by the wide-field (4 sq. deg.) high cadence (Gamma = 6/hr) monitoring of the Korea Microlensing Telescope Network (KMTNet), which had two of its three telescopes in commissioning operations at the time of the planetary anomaly. The source crossing time t_* = 16 min is among the shortest ever published. The high-cadence, wide-field observations that are the hallmark of KMTNet are the only way to routinely capture such short crossings. High-cadence resolution of short caustic crossings will preferentially lead to mass and distance measurements for the lens. This is because the short crossing time typically implies a nearby lens, which enables the measurement of additional effects (bright lens and/or microlens parallax). When combined with the measured crossing time, these effects can yield planet/host masses and distance.

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

In this paper, we demonstrate the severity of the degeneracy between the microlens-parallax and lens-orbital effects by presenting the analysis of the gravitational binary-lens event OGLE-2015-BLG-0768. Despite the obvious deviation from the model based on the the linear observer motion and the static binary, it is found that the residual can be almost equally well explained by either the parallactic motion of the Earth or the rotation of the binary lens axis, resulting in the severe degeneracy between the two effects. We show that the degeneracy can be readily resolved with the additional data provided by space-based microlens parallax observations. Enabling to distinguish between the two higher-order effects, space-based microlens parallax observations will make it possible not only to accurately determine the physical lens parameters but also to further constrain the orbital parameters of binary lenses.

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

We present a statistical assessment of both, observed and reported, photometric uncertainties in the OGLE-IV Galactic bulge microlensing survey data. This dataset is widely used for the detection of variable stars, transient objects, discovery of microlensing events, and characterization of the exo-planetary systems. Large collections of RR Lyrae stars and Cepheids discovered by the OGLE project toward the Galactic bulge provide light curves based on this dataset. We describe the method of analysis, and provide the procedure, which can be used to update preliminary photometric uncertainties, provided with the light curves, to the ones reflecting the actual observed scatter at a given magnitude and for a given CCD detector of the OGLE-IV camera.This is of key importance for data modeling, in particular, for the correct estimation of the goodness of fit.

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

Most stellar remnants so far have been found in binary systems, where they interact with matter from their companions. Isolated neutron stars and black holes are difficult to find as they are dark, yet they are predicted to exist in our Galaxy in vast numbers. We explored the OGLE-III database of 150 million objects observed in years 2001-2009 and found 59 microlensing events exhibiting a parallax effect due to the Earth's motion around the Sun. Combining parallax and brightness measurements from microlensing light curves with expected proper motions in the Milky Way, we identified 13 microlensing events which are consistent with having a white dwarf, neutron star or a black hole lens and we estimated their masses and distances. The most massive of our black hole candidates has 9.3 MSun and is at a distance of 2.4 kpc. The distribution of masses of our candidates indicates a continuum in mass distribution with no mass gap between neutron stars and black holes. We also present predictions on how such events will be observed by the astrometric Gaia mission.

Spitzer microlensing parallax observations of OGLE-2015-BLG-1212 decisively breaks a degeneracy between planetary and binary solutions that is somewhat ambiguous when only ground-based data are considered. Only eight viable models survive out of an initial set of 32 local minima in the parameter space. These models clearly indicate that the lens is a stellar binary system possibly located within the bulge of our Galaxy, ruling out the planetary alternative. We argue that several types of discrete degeneracies can be broken via such space-based parallax observations.

This work describes the analysis of 64 early-type, massive contact or near-contact eclipsing systems from the Large Magellanic Cloud discovered by the OGLE-III survey. It presents the determination of the period-luminosity-colour relation followed by these objects, that is different from the one previously known for late-type W UMa stars. The relation for massive stars has a significantly steeper dependence on the colour, which is related to a much higher bolometric correction, however it is shallower in the period term. This leads to the conclusion, that the relation for the total population of main sequence contact binaries is non-linear. When studied separately, genuinely-contact and near-contact systems follow two slightly different relations.

We present a sample of the most interesting classical Cepheids selected from the OGLE collection of classical Cepheids in the Magellanic System. The main selection criterion for this sample was the presence of non-standard, unique pulsational properties. The sample contains the first known double-mode Cepheid pulsating in the second- and third-overtone modes and a large number of objects with non-radial modes excited. We also found Cepheids revealing Blazhko-like light curve modulation, objects ceasing pulsations or showing atypical shapes of their light curves. Additionally, the status of several triple mode classical Cepheids is updated based on OGLE-IV photometry extending the time baseline to 15 years.

We present here a nearly complete census of classical Cepheids in the Magellanic System. The sample extends the set of Cepheids published in the past by the Optical Gravitational Lensing Experiment (OGLE) to the outer regions of the Large (LMC) and Small Magellanic Cloud (SMC). The entire collection consists of 9535 Cepheids of which 4620 belong to the LMC and 4915 are members of the SMC. We provide the I- and V-band time-series photometry of the published Cepheids, their finding charts, and basic observational parameters. Based on this unique OGLE sample of Cepheids we present updated period--luminosity relations for fundamental, first, and second mode of pulsations in the I- and V-bands and for the W_I extinction-free Wesenheit index. We also show the distribution of classical Cepheids in the Magellanic System. The OGLE collection contains several classical Cepheids in the Magellanic Bridge - the region of interaction between the Magellanic Clouds. The discovery of classical Cepheids and their estimated ages confirm the presence of young stellar population between these galaxies.

We present a census of eclipsing binary systems with classical Cepheid as a component. All such systems known were found in the OGLE collection of classical Cepheids in the Magellanic System. We extend the list of potential candidates adding four new objects found in the OGLE-IV photometric data. One of the new Cepheids in the eclipsing system, OGLE-SMC-CEP-3235, revealed only one eclipse during 15 years of the OGLE photometric monitoring. However, it additionally shows very well pronounced light-time effect indicating that the binarity is real. We also search for the light-time effect in other known eclipsing Cepheids and we clearly detect it in OGLE-LMC-CEP-1812. We discuss application of this tool for the search for Cepheids in non-eclipsing binary systems.

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

We present a detailed study of the classical Cepheid in the double-lined, highly eccentric eclipsing binary system OGLE-LMC562.05.9009. The Cepheid is a fundamental mode pulsator with a period of 2.988 days. The orbital period of the system is 1550 days. Using spectroscopic data from three 4-8-m telescopes and photometry spanning 22 years, we were able to derive the dynamical masses and radii of both stars with exquisite accuracy. Both stars in the system are very similar in mass, radius and color, but the companion is a stable, non-pulsating star. The Cepheid is slightly more massive and bigger (M_1 = 3.70 +/- 0.03M_sun, R_1 = 28.6 +/- 0.2R_sun) than its companion (M_2 = 3.60 +/- 0.03M_sun, R_2 = 26.6 +/- 0.2R_sun). Within the observational uncertainties both stars have the same effective temperature of 6030 +/- 150K. Evolutionary tracks place both stars inside the classical Cepheid instability strip, but it is likely that future improved temperature estimates will move the stable giant companion just beyond the red edge of the instability strip. Within current observational and theoretical uncertainties, both stars fit on a 205 Myr isochrone arguing for their common age. From our model, we determine a value of the projection factor of p = 1.37 +/- 0.07 for the Cepheid in the OGLE-LMC562.05.9009 system. This is the second Cepheid for which we could measure its p-factor with high precision directly from the analysis of an eclipsing binary system, which represents an important contribution towards a better calibration of Baade-Wesselink methods of distance determination for Cepheids.

Based on data from the ongoing OGLE Galaxy Variability Survey (OGLE GVS) we have verified observed properties of stars detected by the near-infrared VVV survey in a direction near the Galactic plane at longitude l~-27 deg and recently tentatively classified as classical Cepheids belonging to a, hence claimed, dwarf galaxy at a distance of about 90 kpc from the Galactic Center. Three of four stars are detected in the OGLE GVS I-band images. We show that two of the objects are not variable at all and the third one with a period of 5.695 d and a nearly sinusoidal light curve of an amplitude of 0.5 mag cannot be a classical Cepheid and is very likely a spotted object. These results together with a very unusual shape of the Ks-band light curve of the fourth star indicate that very likely none of them is a Cepheid and, thus, there is no evidence for a background dwarf galaxy. Our observations show that a great care must be taken when classifying objects by their low-amplitude close-to-sinusoidal near-infrared light curves, especially with a small number of measurements. We also provide a sample of high-amplitude spotted stars with periods of a few days that can mimick pulsations and even eclipses.

We develop a new photometry algorithm that is optimized for $Spitzer$ time series in crowded fields and that is particularly adapted to faint and/or heavily blended targets. We apply this to the 170 targets from the 2015 $Spitzer$ microlensing campaign and present the results of three variants of this algorithm in an online catalog. We present detailed accounts of the application of this algorithm to two difficult cases, one very faint and the other very crowded. Several of $Spitzer$'s instrumental characteristics that drive the specific features of this algorithm are shared by $Kepler$ and $WFIRST$, implying that these features may prove to be a useful starting point for algorithms designed for microlensing campaigns by these other missions.

There are by now ten published detections of fast radio bursts (FRBs), single bright GHz-band millisecond pulses of unknown origin. Proposed explanations cover a broad range from exotic processes at cosmological distances to atmospheric and terrestrial sources. Loeb et al. have previously suggested that FRB sources could be nearby flare stars, and pointed out the presence of a W-UMa-type contact binary within the beam of one out of three FRB fields that they examined. Using time-domain optical photometry and spectroscopy, we now find possible flare stars in additional FRB fields, with one to three such cases among eight FRB fields studied. We evaluate the chance probabilities of these possible associations to be in the range 0.1% to 9%, depending on the input assumptions. Further, we re-analyze the probability that two FRBs recently discovered 3 years apart within the same radio beam are unrelated. Contrary to other claims, we conclude with 99% confidence that the two events are from the same repeating source. The different dispersion measures between the two bursts then rule out a cosmological origin for the dispersion measure, but are consistent with the flare-star scenario with a varying plasma blanket between bursts. Finally, we review some theoretical objections that have been raised against a local flare-star FRB origin, and show that they are incorrect.

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

We present a collection of 250 anomalous Cepheids (ACs) discovered in the OGLE-IV fields toward the Large (LMC) and Small Magellanic Cloud (SMC). The LMC sample is an extension of the OGLE-III Catalog of ACs published in 2008, while the SMC sample contains the first known bona fide ACs in this galaxy. The total sample is composed of 141 ACs in the LMC and 109 ACs in the SMC. All these stars pulsate in single modes: fundamental (174 objects) or first overtone (76 objects). Additionally, we report the discovery of four ACs located in the foreground of the Magellanic Clouds. These are the first fundamental-mode ACs known in the Galactic field. We demonstrate that the coefficients phi_21 and phi_31 determined by the Fourier light curve decomposition are useful discriminators between classical Cepheids and ACs, at least in the LMC and in the field of the Milky Way. In the SMC, the light curve shapes and mean magnitudes of short-period classical Cepheids make them similar to ACs, which is a source of difficulties in the discrimination of both classes of pulsators. The presence of unidentified ACs in the catalogs of classical Cepheids may be partly responsible for the observed non-linearity of the period-luminosity relation observed for short-period Cepheids in the SMC. We compare spatial distributions of ACs, classical Cepheids and RR Lyr stars. We show that the distribution of ACs resembles that of old stars (RR Lyr variables), although in the LMC there are visible structures typical for young population (classical Cepheids): the bar and spiral arms. This may suggest that ACs are a mixture of relatively young stars and mergers of very old stars.

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

We present a sample of 242 ultra-short-period (P < 0.22 d) eclipsing and ellipsoidal binary stars identified in the OGLE fields toward the Galactic bulge. Based on the light curve morphology, we divide the sample into candidates for contact binaries and non-contact binaries. In the latter group we distinguish binary systems consisting of a cool main-sequence star and a B-type subdwarf (HW Vir stars) and candidates for cataclysmic variables, including five eclipsing dwarf novae. One of the detected eclipsing binary systems - OGLE-BLG-ECL-000066 - with the orbital period below 0.1 d, likely consists of M dwarfs in a nearly contact configuration. If confirmed, this would be the shortest-period M-dwarf binary system currently known. We discuss possible evolutionary mechanisms that could lead to the orbital period below 0.1 d in an M-dwarf binary.

We present and study the largest and the most comprehensive catalog of microlensing events ever constructed. The sample of standard microlensing events comprises 3718 unique events from years 2001--2009, with 1409 not detected before in real-time by the Early Warning System of the Optical Gravitational Lensing Experiment (OGLE). The search pipeline makes use of Machine Learning algorithms in order to help find rare phenomena among 150 million objects and derive the detection efficiency. Applications of the catalog can be numerous, from analyzing individual events to large statistical studies for the Galactic mass and kinematics distributions and planetary abundances. We derive the maps of the mean Einstein ring crossing time of events spanning 31 sq. deg. toward of the Galactic Center and compare the observed distributions with the most recent models. We find good agreement within the observed region and we see the signature of the tilt of the bar in the microlensing data. However, the asymmetry of the mean time-scales seems to rise more steeply than predictions, indicating either a somewhat different orientation of the bar or a larger bar width. The map for the events with sources in the Galactic bulge shows a dependence of the mean time-scale on the Galactic latitude, signaling an increasing contribution from disk lenses closer to the plane, related with the height of the disk. Our data present a perfect set for comparing and enhancing new models of the central parts of the Milky Way and creating the 3D picture of the Galaxy.

We report the discovery of two distinct types of ellipsoidal binary systems occupying, so called, sequence E on the period-luminosity (P-L) diagram. We propose that steeper P-L relation is composed of giant-dwarf binaries, while the other consists of giant-giant binary systems. Analysis is based on a sample of 5334 objects, which we select from the OGLE-III survey data toward the the Large Magellanic Cloud. We show that one of the components of ellipsoidal binaries is typically either a Red Clump or a Red Giant Branch star, which leads to clear separation split of the sequence E at P = 40 d. In its short-period part, we identify two subsequences corresponding to the two types of binary systems (E2 and E3), while in the longer-period part the two groups merge forming a single subsequence E1. We extract a group of 271 ellipsoidal systems with eccentric orbits, from our sample. We present the period-luminosity relation they follow.

We present a large sample of eclipsing binary stars detected in the Small Magellanic Cloud fields covering about 14 square degrees that have been monitored for eight years during the third phase of the OGLE survey. This is the largest set of such variables containing 6138 objects, of which 777 are contact and 5361 non-contact binaries. The estimated completeness of this sample is around 82%. We analyze the statistical properties of the sample and present selected interesting objects: 32 systems having eccentric orbit with visible apsidal motion, one Transient Eclipsing Binary, ten RS CVn type stars, 22 still unexplained Double-Periodic Variable stars, and 15 candidates for doubly eclipsing quadruple systems. Based on the OGLE-III proper motions, we classified 47 binaries from our sample as foreground Galactic stars. We also list candidates suitable for the SMC distance determination.