Asteroid Terrestrial-impact Last Alert System Research Articles

A Comparative Study of TESS Light-curve Extraction Methods Applied to Blazars

Blazars are characterized by largely aperiodic variability on timescales ranging from minutes to decades across the electromagnetic spectrum. The Transiting Exoplanet Survey Satellite (TESS) mission provides continuous sampling of blazar variability on timescales ranging from tens of minutes to 27 days for a single sector observation. Proper removal of the background, thermal ramping, and onboard systematic effects are crucial to the extraction of a reliable blazar light curve. Multiple publicly available procedures have been created to correct for these effects. Using ground-based observations from the Zwicky Transit Facility (ZTF) and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as “ground truth” observations, we compare six different methods (regression, cotrending basis vectors, pixel-level decorrelation, eleanor, quaver, and simple differential photometry) to each other, and to our “ground truth” observations, to identify which methods properly correct light curves of a sample of 11 bright blazars, including the suspected neutrino source TXS 0506+056. In addition to comparing the resulting light curves, we compare the slopes of the power spectral densities, perform least-square fitting to simultaneous ZTF data, and quantify other statistical qualities generated from the light curves of each method. We find that only three of the six methods compared (simple differential photometry, eleanor, and quaver) produce TESS light curves consistent with the ground-based ZTF and ATLAS observations.

Main-belt and Trojan asteroid phase curves from the ATLAS survey

ABSTRACT Sparse and serendipitous asteroid photometry obtained by wide field surveys such as the Asteroid Terrestrial-impact Last Alert System (ATLAS) is a valuable resource for studying the properties of large numbers of small Solar System bodies. We have gathered a large data base of ATLAS photometry in wideband optical cyan and orange filters, consisting of 9.6 × 107 observations of 4.5 × 105 main belt asteroids and Jupiter Trojans. We conduct a phase curve analysis of these asteroids considering each apparition separately, allowing us to accurately reject outlying observations and to remove apparitions and asteroids not suitable for phase curve determination. We obtain a data set of absolute magnitudes and phase parameters for over 100 000 selected asteroids observed by ATLAS, $\sim 66\, 000$ of which had sufficient measurements to derive colours in the ATLAS filters. To demonstrate the power of our data set we consider the properties of the Nysa–Polana complex, for which the ATLAS colours and phase parameters trace the S-like and C-like compositions amongst family members. We also compare the properties of the leading and trailing groups of Jupiter Trojans, finding no significant differences in their phase parameters or colours as measured by ATLAS, supporting the consensus that these groups were captured from a common source population during planetary migration. Furthermore, we identify ∼9000 asteroids that exhibit large shifts in derived absolute magnitude between apparitions, indicating that these objects have both elongated shapes and spin axes with obliquity ∼90 degrees.

Discovery of the magnetic cataclysmic variable XMM J152737.4-205305.9 with a deep eclipse-like feature

We report the identification and subsequent examination of a polar-type cataclysmic variable named XMM J152737.4-205305.9 newly discovered with the X-ray Multi-Mirror Mission ( This discovery was made by matching the data archive with the cataclysmic variable candidate catalog provided by Data Release 3 (DR3). The use of X-ray photometry led to the identification of two distinct dips that exhibit a recurring pattern with a precise period of 112.4(1) minutes in two observations made one year apart. The data obtained from photometry provided by the Zwicky Transient Facility (ZTF) and the Asteroid Terrestrial-impact Last Alert System (ATLAS) survey consistently indicate the presence of mass-accretion states that differ by up to 2 mag. Following the optical data, the Spectrum Roentgen Gamma(SRG)/ ero All Sky Survey observed the system at two different X-ray levels, which may imply different accretion states. Following these observations, the low-resolution spectrum obtained using SALT spectroscopy exposes the prominent hydrogen Balmer and helium emission lines, strongly supporting the categorization of this system as a polar-type magnetic cataclysmic variable. The observations conducted at various X-ray levels reveal a consistent pattern of a deep dip-like feature with a width of $ 9.1$ min. This feature implies the presence of an eclipse in both observations. According to data, the object is located at a distance of $\,pc, and its X-ray luminosity lies within the X range.

GW190425: Pan-STARRS and ATLAS coverage of the skymap and limits on optical emission associated with FRB 20190425A

ABSTRACT GW190425 is the second of two binary neutron star (BNS) merger events to be significantly detected by the Laser Interferometer Gravitational Wave (GW) Observatory (LIGO), Virgo and the Kamioka Gravitational Wave (KAGRA) detector network. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real-time searching in 2019. Here, we summarize Asteroid Terrestrial-Impact Last Alert System (ATLAS) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) wide-field optical coverage of the skymap beginning within 1 and 3 h, respectively, of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 20190425A has been suggested, given their spatial and temporal coincidences. The smaller sky localization area of FRB 20190425A and its dispersion measure led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 ± 10 Mpc. Our optical imaging covered the galaxy 6.0 h after GW190425 was detected and 3.5 h after the FRB 20190425A. No optical emission was detected and further imaging at +1.2 and +13.2 d also revealed no emission. If the FRB 20190425A and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 20190425A from a BNS merger involves a supramassive magnetized neutron star spinning down by dipole emission on the time-scale of hours. We show that magnetar-enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB–GW link for this source.

Dust Production Rates in Jupiter-family Comets: A Two Year Study with ATLAS Photometry

Jupiter-family comets (JFCs) exhibit a wide range of activity levels and mass loss over their orbits. We analyzed high-cadence observations of 42 active JFCs with the wide-field Asteroid Terrestrial-impact Last Alert System (ATLAS) in 2020–2021. We measured the dust production rates of the JFCs using the Af ρ parameter and its variation as a function of heliocentric distance. There is a tendency for our JFC sample to exhibit a maximum Af ρ after perihelion, with 254P/McNaught and P/2020 WJ5 (Lemmon) having their maximum Af ρ over a year after perihelion. On average, the rate of change of activity postperihelion was shallower than preperihelion. We also estimated the maximum mass-loss rate for 17 of the JFCs in our sample, finding 4P/Faye to be the most active. We present a subset of comets whose measured Af ρ have been interpolated and extrapolated to a common distance of 2 au preperihelion and postperihelion. From these measurements we found no correlation of intrinsic activity with current perihelion distance. For three of the JFCs in our sample, 6P/d’Arrest, 156P/Russell–LINEAR, and 254P/McNaught, there was no visible coma but a constant absolute magnitude, which we attribute to a probable detection of the nucleus. We derived upper limits for the nuclear radii of ≤2.1 ± 0.3 km, ≤2.0 ± 0.2 km, and ≤4.0 ± 0.8 km, respectively. Finally, we found that 4P/Faye, 108P/Ciffreo, 132P/Helin–Roman–Alu 2, 141P/Machholz 2, and 398P/Boattini experienced outbursts between 2020 and 2022.

Linking Sky-plane Observations of Moving Objects

The Asteroid Terrestrial-impact Last Alert System (ATLAS) observes the visible sky every night in search of dangerous asteroids. With four (soon five) sites ATLAS is facing new challenges for scheduling observations and linking detections to identify moving asteroids. Flexibility in coping with diverse observation sites and times of detections that can be linked is critical, as is optimization of observing time for coverage versus depth. We present new algorithms to fit orbits rapidly to sky-plane observations, and to test and link sets of detections to find the ones which belong to moving objects. The PUMA algorithm for fitting orbits to angular positions on the sky executes in about a millisecond, orders of magnitude faster than the methods currently in use by the community, without sacrifice in accuracy. The PUMA software should be generally useful to anyone who needs to test many sets of detections for consistency with a real orbit. The PUMALINK algorithm to find linkages among sets of detections has similarities to other approaches, notably HelioLinC, but it functions well at asteroid ranges of a small fraction of an astronomical unit. PUMALINK is fast enough to test 10 million possible tracklets against one another in a half hour of computer time. Candidate linkages are checked by the PUMA library to test that the detections correspond to a real orbit, even at close range, and the false alarm rate is manageable. Sky surveys that produce large numbers of detections from large numbers of exposures may find the PUMALINK software helpful. We present the results of tests of PUMALINK on three data sets which illustrate PUMALINK’s effectiveness and economy: 2 weeks of all ATLAS detections over the sky, 2 weeks of special ATLAS opposition observations with long exposure time, and 2 weeks of simulated LSST asteroid observations. Detection probabilities of linkages must be traded against false alarm rate, but a representative choice for PUMALINK might be 90% detection probability for real objects while keeping the false alarm rate below 10% for a 100:1 population of false:real. Although optimization of the tradeoffs between detection probability, execution time, and false alarm rate is application specific and beyond the scope of this paper, we provide guidance on methods to distinguish false alarms from correct linkages of real objects.

Investigating Possible Correlations between Gamma-Ray and Optical Lightcurves for TeV-Detected Northern Blazars over 8 Years of Observations

Blazars are a subclass of active galactic nuclei (AGN) having relativistic jets aligned within a few degrees of our line-of-sight and form the majority of the AGN detected in the TeV regime. The Fermi-Large Area Telescope (LAT) is a pair-conversion telescope, sensitive to photons having energies between 20 MeV and 2 TeV, and is capable of scanning the entire gamma-ray sky every three hours. Despite the remarkable success of the Fermi mission, many questions still remain unanswered, such as the site of gamma-ray production and the emission mechanisms involved. The Asteroid Terrestrial-impact Last Alert System (ATLAS) is a high cadence all sky survey system optimized to be efficient for finding potentially dangerous asteroids, as well as in tracking and searching for highly variable and transient sources, such as AGN. In this study, we investigate possible correlations between the Fermi-LAT observations in the 100 MeV–300 GeV energy band and the ATLAS optical data in the R-band, centered at 679 nm, for a sample of 18 TeV-detected northern blazars over 8 years of observations between 2015 and 2022. Under the assumption that the optical and gamma-ray flares are produced by the same outburst propagating down the jet, the strong correlations found for some sources suggest a single-zone leptonic model of emission.

The value-added catalogue of ASAS-SN eclipsing binaries – II. Properties of extra-physics systems

ABSTRACT Detached eclipsing binaries are the primary tool used to measure precise masses and radii of stars. In our previous paper estimating the parameters of more than 30 000 detached eclipsing binaries, we identified 766 eclipsing binaries with additional features in their All-Sky Automated Survey for Supernovae (ASAS-SN) and Transiting Exoplanet Survey Satellite light curves. Here, we characterize these ‘extra-physics’ systems, identifying eclipsing binaries with spotted stars, pulsating components, and candidate triple/quadruple systems. We use the Gaia, Asteroid Terrestrial-impact Last Alert System (ATLAS), Zwicky Transient Facility (ZTF), and ASAS-SN variable star catalogues to search for possible blends. We use MIST (MESA Isochrones and Stellar Tracks) isochrones and evolutionary tracks to identify systems with main-sequence, subgiant, and giant primaries and highlight systems in sparsely populated regions of the colour–magnitude diagram. We find that the orbital period distribution of spotted binaries is divided by evolutionary state and find 68 with X-ray detections. For the candidate triple/quadruples and pulsating systems, we calculate the extra orbital/pulsational period and identify systems with resonances. Finally, we highlight a number of exotic systems, including eclipsing cataclysmic variables, subdwarfs, and binaries with discs.

SCAT uncovers ATLAS’s first tidal disruption event ATLAS18mlw: a faint and fast TDE in a quiescent Balmer strong Galaxy

ABSTRACT We present the discovery that ATLAS18mlw was a tidal disruption event (TDE) in the galaxy WISEA J073544.83+663717.3, at a luminosity distance of 334 Mpc. Initially discovered by the Asteroid Terrestrial Impact Last Alert System (ATLAS) on 2018 March 17.3, the TDE nature of the transient was uncovered only recently with the re-reduction of a SuperNova Integral Field Spectrograph (SNIFS) spectrum. This spectrum, taken by the Spectral Classification of Astronomical Transients (SCAT) survey, shows a strong blue continuum and a broad H α emission line. Here, we present roughly 6 yr of optical survey photometry beginning before the TDE to constrain active galactic nucleus activity, optical spectroscopy of the transient, and a detailed study of the host galaxy properties through analysis of archival photometry and a host spectrum. ATLAS18mlw was detected in ground-based light curves for roughly 2 months. From a blackbody fit to the transient spectrum and bolometric correction of the optical light curve, we conclude that ATLAS18mlw is best explained by a low-luminosity TDE with a peak luminosity of log(L [erg s−1]) = 43.5 ± 0.2. The TDE classification is further supported by the quiescent Balmer strong nature of the host galaxy. We also calculated the TDE decline rate from the bolometric light curve and find ΔL40 = −0.7 ± 0.2 dex, making ATLAS18mlw a member of the growing class of ‘faint and fast’ TDEs with low peak luminosities and fast decline rates.

What can Gaussian processes really tell us about supernova light curves? Consequences for Type II(b) morphologies and genealogies

ABSTRACT Machine learning has become widely used in astronomy. Gaussian process (GP) regression in particular has been employed a number of times to fit or resample supernova (SN) light curves, however by their nature typical GP models are not suited to fit SN photometric data and they will be prone to overfitting. Recently GP resampling was used in the context of studying the morphologies of Type II and IIb SNe and they were found to be clearly distinct with respect to four parameters: the rise time (trise), the magnitude difference between 40 and 30 d post-explosion (Δm40–30), the earliest maximum (post-peak) of the first derivative (dm1), and minimum of the second derivative (dm2). Here we take a close look at GP regression and its limitations in the context of SN light curves in general, and we also discuss the uncertainties on these specific parameters, finding that dm1 and dm2 cannot give reliable astrophysical information. We do reproduce the clustering in trise–Δm40–30 space, although it is not as clear cut as previously presented. The best strategy to accurately populate the trise–Δm40–30 space will be to use an expanded sample of high-quality light curves [such as those in the Asteroid Terrestrial-impact Last Alert System (ATLAS) transient survey] and analytical fitting methods. Finally, using the bpass fiducial models, we predict that future photometric studies will reveal clear clustering of the Type IIb and II light curve morphologies with a distinct continuum of transitional events.

The ‘Giraffe’: discovery of a stripped red giant in an interacting binary with an ∼2 M⊙ lower giant

ABSTRACT We report the discovery of a stripped giant + lower giant Algol-type binary, 2M04123153+6738486 (2M0412), identified during a search for non-interacting compact object – star binaries. 2M0412 is an evolved (Teff, giant ≃ 4000 K), luminous (Lgiant ≃ 150 L⊙) red giant in a circular P = 81.2 d binary. 2M0412 is a known variable star previously classified as a semiregular variable. The cross-correlation functions of follow-up Keck/HIRES and LBT/Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) spectra show an radial velocity variable second component with implied mass ratio q = Mgiant/Mcomp ≃ 0.20 ± 0.01. The All-Sky Automated Survey for SuperNovae (ASAS-SN), Asteroid Terrestrial-impact Last Alert System (ATLAS), Transiting Exoplanet Survey Telescope ( TESS), and Zwicky Transient Facility (ZTF) light curves show that the giant is a Roche lobe filling ellipsoidal variable with an inclination of 49.4 ± 0.3°, and a giant mass of Mgiant = 0.38 ± 0.01 M⊙ for a distance of ≃ 3.7 kpc. The mass of the giant indicates that its envelope has been partially stripped. The giant companion on the lower red giant branch has a mass of Mcomp = 1.91 ± 0.03 M⊙ with Teff, comp ≃ 5000 K, Lcomp ≃ 60 L⊙, and Rcomp ≃ 11 R⊙. The lower giant contributes ${\sim }35{{\ \rm per\ cent}}$ of the flux in the V band. We also identify an orbital phase dependent, broad $\rm H\,\alpha$ emission line which could indicate ongoing accretion from the stripped red giant on to the companion.

Discovery of superslow rotating asteroids with ATLAS and ZTF photometry

ABSTRACT We present here the discovery of a new class of superslow rotating asteroids (Prot ≳1000 h) in data extracted from the Asteroid Terrestrial-impact Last Alert System (ATLAS) and Zwicky Transient Facility (ZTF) all-sky surveys. Of the 39 rotation periods we report here, 32 have periods longer than any previously reported unambiguous rotation periods currently in the Asteroid Light Curve Data base. In our sample, seven objects have a rotation period >4000 h and the longest period we report here is 4812 h (∼200 d). We do not observe any correlation between taxonomy, albedo, or orbital properties with superslow rotating status. The most plausible mechanism for the creation of these very slow rotators is if their rotations were slowed by YORP spin-down. Superslow rotating asteroids may be common, with at least 0.4 per cent of the main-belt asteroid population with a size range between 2 and 20 km in diameter rotating with periods longer than 1000 h.

CONFIRMATION OF SEVEN FAINT ATLAS VARIABLE STAR CANDIDATES

A pilot survey conducted at the Lookout Observatory has confirmed seven faint (V ~ 13 to 17) variables in the region of Kepler-76b that were recently discovered by the Asteroid Terrestrial-Impact Last Alert System (ATLAS). The ATLAS survey identified 315,000 probably variables within its wide-field survey in 2018. The faintness (down to r ~ 18) and small amplitudes (down to 0.02 mag) included in these candidates makes external validation difficult. Our confirmation of seven such variable stars gives credibility to the ATLAS list. Further, the agreement between various surveys and LO data validates the use of our new survey for variable star and exoplanet research.

Asteroid models reconstructed from ATLAS photometry

Context. The Asteroid Terrestrial-impact Last Alert System (ATLAS) is an all-sky survey primarily aimed at detecting potentially hazardous near-Earth asteroids. Apart from the astrometry of asteroids, it also produces their photometric measurements that contain information about asteroid rotation and their shape. Aims. To increase the current number of asteroids with a known shape and spin state, we reconstructed asteroid models from ATLAS photometry that was available for approximately 180 000 asteroids observed between 2015 and 2018. Methods. We made use of the light-curve inversion method implemented in the Asteroids@home project to process ATLAS photometry for roughly 100 000 asteroids with more than a hundred individual brightness measurements. By scanning the period and pole parameter space, we selected those best-fit models that were, according to our setup, a unique solution for the inverse problem. Results. We derived ~2750 unique models, 950 of them were already reconstructed from other data and published. The remaining 1800 models are new. About half of them are only partial models, with an unconstrained pole ecliptic longitude. Together with the shape and spin, we also determined for each modeled asteroid its color index from the cyan and orange filter used by the ATLAS survey. We also show the correlations between the color index, albedo, and slope of the phase-angle function. Conclusions. The current analysis is the first inversion of ATLAS asteroid photometry, and it is the first step in exploiting the huge scientific potential that ATLAS photometry has. ATLAS continues to observe, and in the future, this data, together with other independent photometric measurements, can be inverted to produce more refined asteroid models.

Design and Operation of the ATLAS Transient Science Server

The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5 m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above ) every two nights, exposing four times per night, typically reaching magnitude per exposure when the moon is illuminated and magnitude per exposure in dark skies. Construction is underway of two further units to be sited in Chile and South Africa which will result in an all-sky daily cadence from 2021. Initially designed for detecting potentially hazardous near earth objects, the ATLAS data enable a range of astrophysical time domain science. To extract transients from the data stream requires a computing system to process the data, assimilate detections in time and space and associate them with known astrophysical sources. Here we describe the hardware and software infrastructure to produce a stream of clean, real, astrophysical transients in real time. This involves machine learning and boosted decision tree algorithms to identify extragalactic and Galactic transients. Typically we detect 10–15 supernova candidates per night which we immediately announce publicly. The ATLAS discoveries not only enable rapid follow-up of interesting sources but will provide complete statistical samples within the local volume of 100 Mpc. A simple comparison of the detected supernova rate within 100 Mpc, with no corrections for completeness, is already significantly higher (factor 1.5 to 2) than the current accepted rates.

Beyond Gaia: Asteroseismic Distances of M Giants Using Ground-based Transient Surveys

Abstract Evolved stars near the tip of the red giant branch show solar-like oscillations with periods spanning hours to months and amplitudes ranging from ∼1 mmag to ∼100 mmag. The systematic detection of the resulting photometric variations with ground-based telescopes would enable the application of asteroseismology to a much larger and more distant sample of stars than is currently accessible with space-based telescopes such as Kepler or the ongoing Transiting Exoplanet Survey Satellite mission. We present an asteroseismic analysis of 493 M giants using data from two ground-based surveys: the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the All-Sky Automated Survey for Supernovae (ASAS-SN). By comparing the extracted frequencies with constraints from Kepler, the Sloan Digital Sky Survey Apache Point Observatory Galaxy Evolution Experiment, and Gaia we demonstrate that ground-based transient surveys allow accurate distance measurements to oscillating M giants with a precision of ∼15%. Using stellar population synthesis models we predict that ATLAS and ASAS-SN can provide asteroseismic distances to ∼2 × 106 galactic M giants out to typical distances of 20–50 kpc, vastly improving the reach of Gaia and providing critical constraints for Galactic archeology and galactic dynamics.

Investigating Taxonomic Diversity within Asteroid Families through ATLAS Dual-band Photometry

Abstract We present here the c–o colors for identified Flora, Vesta, Nysa–Polana, Themis, and Koronis family members within the historic data set (2015–2018) of the Asteroid Terrestrial-impact Last Alert System (ATLAS). The Themis and Koronis families are known to be relatively pure C- and S-type Bus-DeMeo taxonomic families, respectively, and the extracted color data from the ATLAS broadband c- and o-filters of these two families are used to demonstrate that the ATLAS c–o color is a sufficient parameter to distinguish between the C- and S-type taxonomies. The Vesta and Nysa–Polana families are known to display a mixture of taxonomies possibly due to Vesta’s differentiated parent body origin and Nysa–Polana actually consisting of two nested families with differing taxonomies. Our data show that the Flora family also displays a large degree of taxonomic mixing and the data reveal a substantial H-magnitude dependence on color. We propose and exclude several interpretations for the observed taxonomic mix. Additionally, we extract rotation periods of all of the targets reported here and find good agreement with targets that have previously reported periods.

The New EXor Outburst of ESO-Hα 99 Observed by Gaia ATLAS and TESS

Abstract We report photometry and spectroscopy of the outburst of the young stellar object ESO-Hα99. The outburst was first noticed in Gaia alert Gaia18dvc and later by the Asteroid Terrestrial-impact Last Alert System (ATLAS). We have established the outburst light curve with archival ATLAS orange filter photometry, Gaia data, new V-band photometry, and J, H, and K s photometry from the Infrared Imaging System (IRIS) and the United Kingdom Infrared Telescope (UKIRT). The brightness has fluctuated several times near the light curve maximum. The Transiting Exoplanet Survey Satellite (TESS) satellite observed ESO-Hα 99 with high cadence during one of these minor minima and found brightness fluctuations on timescales of days and hours. Imaging with UKIRT shows the outline of an outflow cavity, and we find one knot of emission, now named MHO 1520, on the symmetry axis of this nebula, indicating recent collimated outflow activity from ESO-Hα 99. Its pre-outburst SED shows a flat far-infrared spectrum, confirming its early evolutionary state and its similarity to other deeply embedded objects in the broader EXor class. The pre-outburst luminosity is 34 L ⊙, a much higher luminosity than typical EXors, indicating that ESO-Hα 99 may be a star of intermediate mass. Infrared and optical spectroscopy show a rich emission-line spectrum, including H i lines, strong red Ca ii emission, as well as infrared CO bandhead emission, all characteristic EXors in the broadest sense. Comparison of the present spectra with an optical spectrum obtained in 1993, presumably in the quiescent state of the object, shows that during the present outburst the continuum component of the spectrum has increased notably more than the emission lines. The Hα equivalent width during the outburst is down to one-half of its 1993 level, and shock-excited emission lines are much less prominent.

The EREBOS project: Investigating the effect of substellar and low-mass stellar companions on late stellar evolution

Eclipsing post-common-envelope binaries are highly important for resolving the poorly understood, very short-lived common-envelope phase of stellar evolution. Most hot subdwarfs (sdO/Bs) are the bare helium-burning cores of red giants that have lost almost all of their hydrogen envelope. This mass loss is often triggered by common-envelope interactions with close stellar or even substellar companions. Cool companions to hot subdwarf stars such as late-type stars and brown dwarfs are detectable from characteristic light-curve variations – reflection effects and often eclipses. In the recently published catalog of eclipsing binaries in the Galactic Bulge and in the Asteroid Terrestrial-impact Last Alert System (ATLAS) survey, we discovered 125 new eclipsing systems showing a reflection effect seen by visual inspection of the light curves and using a machine-learning algorithm, in addition to the 36 systems previously discovered by the Optical Gravitational Lesing Experiment (OGLE) team. The Eclipsing Reflection Effect Binaries from Optical Surveys (EREBOS) project aims at analyzing all newly discovered eclipsing binaries of the HW Vir type (hot subdwarf + close, cool companion) based on a spectroscopic and photometric follow up to derive the mass distribution of the companions, constrain the fraction of substellar companions, and determine the minimum mass needed to strip off the red-giant envelope. To constrain the nature of the primary we derived the absolute magnitude and the reduced proper motion of all our targets with the help of the parallaxes and proper motions measured by the Gaia mission and compared those to the Gaia white-dwarf candidate catalog. It was possible to derive the nature of a subset of our targets, for which observed spectra are available, by measuring the atmospheric parameter of the primary, confirming that less than 10% of our systems are not sdO/Bs with cool companions but are white dwarfs or central stars of planetary nebula. This large sample of eclipsing hot subdwarfs with cool companions allowed us to derive a significant period distribution for hot subdwarfs with cool companions for the first time showing that the period distribution is much broader than previously thought and is ideally suited to finding the lowest-mass companions to hot subdwarf stars. The comparison with related binary populations shows that the period distribution of HW Vir systems is very similar to WD+dM systems and central stars of planetary nebula with cool companions. In the future, several new photometric surveys will be carried out, which will further increase the sample of this project, providing the potential to test many aspects of common-envelope theory and binary evolution.

PS18kh: A New Tidal Disruption Event with a Non-axisymmetric Accretion Disk

Abstract We present the discovery of PS18kh, a tidal disruption event discovered at the center of SDSS J075654.53+341543.6 (d ≃ 322 Mpc) by the Pan-STARRS Survey for Transients. Our data set includes pre-discovery survey data from Pan-STARRS, the All-sky Automated Survey for Supernovae, and the Asteroid Terrestrial-impact Last Alert System as well as high-cadence, multiwavelength follow-up data from ground-based telescopes and Swift, spanning from 56 days before peak light until 75 days after. The optical/UV emission from PS18kh is well-fit as a blackbody with temperatures ranging from T ≃ 12,000 K to T ≃ 25,000 K and it peaked at a luminosity of L ≃ 8.8 × 1043 erg s−1. PS18kh radiated E = (3.45 ± 0.22) × 1050 erg over the period of observation, with (1.42 ± 0.20) × 1050 erg being released during the rise to peak. Spectra of PS18kh show a changing, boxy/double-peaked Hα emission feature, which becomes more prominent over time. We use models of non-axisymmetric accretion disks to describe the profile of the Hα line and its evolution. We find that at early times the high accretion rate leads the disk to emit a wind which modifies the shape of the line profile and makes it bell-shaped. At late times, the wind becomes optically thin, allowing the non-axisymmetric perturbations to show up in the line profile. The line-emitting portion of the disk extends from r in ∼ 60r g to an outer radius of r out ∼ 1400r g and the perturbations can be represented either as an eccentricity in the outer rings of the disk or as a spiral arm in the inner disk.