A 12.4-day periodicity in a close binary system after a supernova.

Neutron stars and stellar-mass black holes are the remnants of massive star explosions1. Most massive stars reside in close binary systems2, and the interplay between the companion star and the newly formed compact object has been theoretically explored3, but signatures for binarity or evidence for the formation of a compact object during a supernova explosion are still lacking. Here we report a stripped-envelope supernova, SN 2022jli, which shows 12.4-day periodic undulations during the declining light curve. Narrow Hα emission is detected in late-time spectra with concordant periodic velocity shifts, probably arising from hydrogen gas stripped from a companion and accreted onto the compact remnant. A new Fermi-LAT γ-ray source is temporally and positionally consistent with SN 2022jli. The observed properties of SN 2022jli, including periodic undulations in the optical light curve, coherent Hα emission shifting and evidence for association with a γ-ray source, point to the explosion of a massive star in a binary system leaving behind a bound compact remnant. Mass accretion from the companion star onto the compact object powers the light curve of the supernova and generates the γ-ray emission.

A rotating white dwarf shows different compositions on its opposite faces.

White dwarfs, the extremely dense remnants left behind by most stars after their death, are characterized by a mass comparable to that of the Sun compressed into the size of an Earth-like planet. In the resulting strong gravity, heavy elements sink towards the centre and the upper layer of the atmosphere contains only the lightest element present, usually hydrogen or helium1,2. Several mechanisms compete with gravitational settling to change a white dwarf's surface composition as it cools3, and the fraction of white dwarfs with helium atmospheres is known to increase by a factor of about 2.5 below a temperature of about 30,000 kelvin4-8; therefore, some white dwarfs that appear to have hydrogen-dominated atmospheres above 30,000 kelvin are bound to transition to be helium-dominated as they cool below it. Here we report observations of ZTF J203349.8+322901.1, a transitioning white dwarf with two faces: one side of its atmosphere is dominated by hydrogen and the other one by helium. This peculiar nature is probably caused by the presence of a small magnetic field, which creates an inhomogeneity in temperature, pressure or mixing strength over the surface9-11. ZTF J203349.8+322901.1 might be the most extreme member of a class of magnetic, transitioning white dwarfs-together with GD 323 (ref. 12), a white dwarf that shows similar but much more subtle variations. This class of white dwarfs could help shed light on the physical mechanisms behind the spectral evolution of white dwarfs.

An infrared transient from a star engulfing a planet.

Planets with short orbital periods (roughly under 10 days) are common around stars like the Sun1,2. Stars expand as they evolve and thus we expect their close planetary companions to be engulfed, possibly powering luminous mass ejections from the host star3-5. However, this phase has never been directly observed. Here we report observations of ZTF SLRN-2020, a short-lived optical outburst in the Galactic disk accompanied by bright and long-lived infrared emission. The resulting light curve and spectra share striking similarities with those of red novae6,7-a class of eruptions now confirmed8 to arise from mergers of binary stars. Its exceptionally low optical luminosity (approximately 1035 erg s-1) and radiated energy (approximately 6.5 × 1041 erg) point to the engulfment of a planet of fewer than roughly ten Jupiter masses by its Sun-like host star. We estimate the Galactic rate of such subluminous red novae to be roughly between 0.1 and several per year. Future Galactic plane surveys should routinely identify these, showing the demographics of planetary engulfment and the ultimate fate of planets in the inner Solar System.

A dense 0.1-solar-mass star in a 51-minute-orbital-period eclipsing binary.

Of more than a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes1-9. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution before interacting with the white dwarf10-14, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs13-18. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute-orbital-period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater owing to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modelling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs.

Candidate Tidal Disruption Event AT2019fdr Coincident with a High-Energy Neutrino.

The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time x-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources.

A 62-minute orbital period black widow binary in a wide hierarchical triple.

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original 'black widow', the eight-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20)1, high-energy emission originating from the pulsar2 is irradiating and may eventually destroy3 a low-mass companion. These systems are not only physical laboratories that reveal the interesting results of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars4, allowing for robust tests of the neutron star equation of state. Here we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate, the optical flux of which varies by a factor of more than ten. ZTF J1406+1222 pushes the boundaries of evolutionary models5, falling below the 80-minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low-metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic Centre, making it a probe of formation channels, neutron star kick physics6 and binary evolution.

General relativistic orbital decay in a seven-minute-orbital-period eclipsing binary system.

General relativity1 predicts that short-orbital-period binaries emit considerable amounts of gravitational radiation. The upcoming Laser Interferometer Space Antenna2 (LISA) is expected to detect tens of thousands of such systems3 but few have been identified4, of which only one5 is eclipsing-the double-white-dwarf binary SDSS J065133.338+284423.37, which has an orbital period of 12.75 minutes. Here we report the discovery of an eclipsing double-white-dwarf binary system, ZTF J153932.16+502738.8, with an orbital period of 6.91 minutes. This system has an orbit so compact that the entire binary could fit within the diameter of the planet Saturn. The system exhibits a deep eclipse, and a double-lined spectroscopic nature. We see rapid orbital decay, consistent with that expected from general relativity. ZTF J153932.16+502738.8 is a strong source of gravitational radiation close to the peak of LISA's sensitivity, and we expect it to be detected within the first week of LISA observations, once LISA launches in approximately 2034.

Bringing the visible universe into focus with Robo-AO.

The angular resolution of ground-based optical telescopes is limited by the degrading effects of the turbulent atmosphere. In the absence of an atmosphere, the angular resolution of a typical telescope is limited only by diffraction, i.e., the wavelength of interest, λ, divided by the size of its primary mirror's aperture, D. For example, the Hubble Space Telescope (HST), with a 2.4-m primary mirror, has an angular resolution at visible wavelengths of ~0.04 arc seconds. The atmosphere is composed of air at slightly different temperatures, and therefore different indices of refraction, constantly mixing. Light waves are bent as they pass through the inhomogeneous atmosphere. When a telescope on the ground focuses these light waves, instantaneous images appear fragmented, changing as a function of time. As a result, long-exposure images acquired using ground-based telescopes--even telescopes with four times the diameter of HST--appear blurry and have an angular resolution of roughly 0.5 to 1.5 arc seconds at best. Astronomical adaptive-optics systems compensate for the effects of atmospheric turbulence. First, the shape of the incoming non-planar wave is determined using measurements of a nearby bright star by a wavefront sensor. Next, an element in the optical system, such as a deformable mirror, is commanded to correct the shape of the incoming light wave. Additional corrections are made at a rate sufficient to keep up with the dynamically changing atmosphere through which the telescope looks, ultimately producing diffraction-limited images. The fidelity of the wavefront sensor measurement is based upon how well the incoming light is spatially and temporally sampled. Finer sampling requires brighter reference objects. While the brightest stars can serve as reference objects for imaging targets from several to tens of arc seconds away in the best conditions, most interesting astronomical targets do not have sufficiently bright stars nearby. One solution is to focus a high-power laser beam in the direction of the astronomical target to create an artificial reference of known shape, also known as a 'laser guide star'. The Robo-AO laser adaptive optics system, employs a 10-W ultraviolet laser focused at a distance of 10 km to generate a laser guide star. Wavefront sensor measurements of the laser guide star drive the adaptive optics correction resulting in diffraction-limited images that have an angular resolution of ~0.1 arc seconds on a 1.5-m telescope.

Study of a MEMS-based Shack-Hartmann wavefront sensor with adjustable pupil sampling for astronomical adaptive optics.

We introduce a Shack-Hartmann wavefront sensor for adaptive optics that enables dynamic control of the spatial sampling of an incoming wavefront using a segmented mirror microelectrical mechanical systems (MEMS) device. Unlike a conventional lenslet array, subapertures are defined by either segments or groups of segments of a mirror array, with the ability to change spatial pupil sampling arbitrarily by redefining the segment grouping. Control over the spatial sampling of the wavefront allows for the minimization of wavefront reconstruction error for different intensities of guide source and different atmospheric conditions, which in turn maximizes an adaptive optics system's delivered Strehl ratio. Requirements for the MEMS devices needed in this Shack-Hartmann wavefront sensor are also presented.

Large-scale wave-front reconstruction for adaptive optics systems by use of a recursive filtering algorithm.

We propose a new recursive filtering algorithm for wave-front reconstruction in a large-scale adaptive optics system. An embedding step is used in this recursive filtering algorithm to permit fast methods to be used for wave-front reconstruction on an annular aperture. This embedding step can be used alone with a direct residual error updating procedure or used with the preconditioned conjugate-gradient method as a preconditioning step. We derive the Hudgin and Fried filters for spectral-domain filtering, using the eigenvalue decomposition method. Using Monte Carlo simulations, we compare the performance of discrete Fourier transform domain filtering, discrete cosine transform domain filtering, multigrid, and alternative-direction-implicit methods in the embedding step of the recursive filtering algorithm. We also simulate the performance of this recursive filtering in a closed-loop adaptive optics system.

Comparison of measurements of the outer scale of turbulence by three different techniques.

We have made simultaneous and nearly simultaneous measurements of L0, the outer scale of turbulence, at the Palomar Observatory by using three techniques: angle-of-arrival covariance measurements with the Generalized Seeing Monitor (GSM), differential-image-motion measurements with the adaptive-optics system on the Hale 5-m telescope, and fringe speed measurements with the Palomar Testbed Interferometer (PTI). The three techniques give consistent results, an outer scale of approximately 10-20 m, despite the fact that the spatial scales of the three instruments vary from 1 m for the GSM to 100 m for the PTI.

Fast wave-front reconstruction by solving the Sylvester equation with the alternating direction implicit method.

Large degree-of-freedom real-time adaptive optics (AO) control requires reconstruction algorithms that are computationally efficient and readily parallelized for hardware implementation. In particular, we find the wave-front reconstruction for the Hudgin and Fried geometry can be cast into a form of the well-known Sylvester equation using the Kronecker product properties of matrices. We derive the filters and inverse filtering formulas for wave-front reconstruction in two-dimensional (2-D) Discrete Cosine Transform (DCT) domain for these two geometries using the Hadamard product concept of matrices and the principle of separable variables. We introduce a recursive filtering (RF) method for the wave-front reconstruction on an annular aperture, in which, an imbedding step is used to convert an annular-aperture wave-front reconstruction into a squareaperture wave-front reconstruction, and then solving the Hudgin geometry problem on the square aperture. We apply the Alternating Direction Implicit (ADI) method to this imbedding step of the RF algorithm, to efficiently solve the annular-aperture wave-front reconstruction problem at cost of order of the number of degrees of freedom, O(n). Moreover, the ADI method is better suited for parallel implementation and we describe a practical real-time implementation for AO systems of order 3,000 actuators.