Search

The Einstein Probe (EP) achieved its first detection and localization of a bright X-ray flare, EP240219a, on 2024 February 19, during its commissioning phase. Subsequent targeted searches triggered by the EP240219a alert identified a faint, untriggered gamma-ray burst (GRB) in the archived data of Fermi Gamma-ray Burst Monitor (GBM), Swift Burst Alert Telescope (BAT), and Insight-HXMT/HE. The EP Wide-field X-ray Telescope (WXT) light curve reveals a long duration of approximately 160 s with a slow decay, whereas the Fermi/GBM light curve shows a total duration of approximately 70 s. The peak in the Fermi/GBM light curve occurs slightly later with respect to the peak seen in the EP/WXT light curve. Our spectral analysis shows that a single cutoff power-law (PL) model effectively describes the joint EP/WXT–Fermi/GBM spectra in general, indicating coherent broad emission typical of GRBs. The model yielded a photon index of ∼–1.70 ± 0.05 and a peak energy of ∼257 ± 134 keV. After detection of GRB 240219A, long-term observations identified several candidates in optical and radio wavelengths, none of which was confirmed as the afterglow counterpart during subsequent optical and near-infrared follow-ups. The analysis of GRB 240219A classifies it as an X-ray-rich GRB (XRR) with a high peak energy, presenting both challenges and opportunities for studying the physical origins of X-ray flashes, XRRs, and classical GRBs. Furthermore, linking the cutoff PL component to nonthermal synchrotron radiation suggests that the burst is driven by a Poynting flux-dominated outflow.

Previous investigations of ionospheric electrodynamical responses to solar flares primarily focused on the main phases (MPs) of solar flares. Typical solar irradiance models for driving global ionosphere models do not include the extreme ultraviolet (EUV) late phase (ELP) of flares, which was recently observed with new high-quality solar EUV spectra. Thus, it is still unclear how ionospheric electrodynamics respond to the flare ELP. Here, we analyzed the ionospheric electrodynamical response to the MP and ELP of the X9.3 flare on 2017 September 6, using observations from ground magnetometers, along with simulation results from an ionosphere–thermosphere coupled model. Observations indicated an intensification of the dayside eastward equatorial electrojet (EEJ) by approximately 12 nT at the ELP peak as compared to the quiet day reference. Additionally, the dayside eastward electric field increased due to the ELP, which is different from the reduction of dayside electric fields during MP. The upward E × B plasma drifts decreased by 2.5 m s ^–1 during MP but increased by 0.75 m s ^–1 during the ELP. Altitude-dependent responses of ionospheric conductivities to the ELP modulated the relative contribution of the E- and F-region wind dynamo to zonal electric fields, resulting in an overall increase in the daytime eastward electric fields. Furthermore, combined effects of electric fields and conductivities enhancements contributed to EEJ intensification during the ELP. This study enhances our understanding of how solar flares with ELP change global ionospheric electric fields and currents.

The first source catalog of the Large High Altitude Air Shower Observatory (LHAASO) reported the detection of a very high energy gamma-ray source, 1LHAASO J1219+2915. This Letter presents a further detailed study of the spectral and temporal behavior of this pointlike source. The best-fit position of the TeV source (R.A. = 185.°05 ± 0.°04, decl. = 29.°25 ± 0.°03) is compatible with NGC 4278 within ∼0.°03. Variation analysis shows an indication of variability on a timescale of a few months in the TeV band, which is consistent with low-frequency observations. Based on these observations, we report the detection of TeV γ -ray emissions from this low-luminosity active galactic nucleus. The observation by LHAASO's Water Cherenkov Detector Array during the active period has a significance level of 8.8 σ with a best-fit photon spectral index Γ = 2.56 ± 0.14 and a flux f _1–10 TeV = (7.0 ± 1.1 _sta ± 0.35 _syst ) × 10 ^−13 photons cm ^−2 s ^−1 , or approximately 5% of the Crab Nebula. The discovery of VHE gamma-ray emission from NGC 4278 indicates that compact, weak radio jets can efficiently accelerate particles and emit TeV photons.

The progenitors of gamma-ray bursts (GRBs) have long been an unresolved issue. GRB 230307A stands out as an exceptionally bright event, belonging to the long-duration GRBs but also exhibiting a late-emission component reminiscent of a kilonova. Together with the similar events GRBs 060614 and 211211A, they make up a new subgroup of GRBs with intriguing progenitors. If such long-duration merger-type GRBs originated from the coalescence of a white dwarf (WD) with a neutron star (NS) or a black hole (BH), as proposed in the recent literature, then the larger tidal disruption radius of the WD, together with a nonnegligible residual orbital eccentricity, would make repeated partial tidal disruptions inevitable. This may modulate the mass accretion and jet launching process at the NS or BH, resulting in a quasiperiodic modulation (QPM) in the light curve of the GRB, with a period equal to the orbital period. The detection of potential QPMs during the early episode of prompt emission of these three GRBs supports this scenario, and the relatively slow QPM (> 1 s) suggests that the lighter object cannot be an NS. We propose that the progenitor system of GRBs 230307A, 060614, and 211211A consist of a WD of mass 1.3 M _⊙ , 0.9 M _⊙ , and 1.4 M _⊙ , respectively, and an NS (or BH). After several cycles of modulations, the WD is completely destroyed, and the accretion of the remaining debris dominates the extended emission episode.

Gamma-ray bursts (GRBs) are believed to launch relativistic jets, which generate prompt emission by internal processes, and produce long-lasting afterglows by driving external shocks into the surrounding medium. However, how the jet powers the external shock is poorly known. The unprecedented observations of the keV–MeV emission with Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor and the TeV emission with LHAASO of the brightest-of-all-time GRB 221009A offer a great opportunity to study the prompt-to-afterglow transition and the impact of jet on the early dynamics of external shock. In this Letter, we find that the cumulative light curve of keV–MeV emission could well fit the rising stage of the TeV light curve of GRB 221009A, with a time delay, ${4.45}_{-0.26}^{+0.26}$ s, of TeV emission. Moreover, both the rapid increase in the initial stage and the excess from about T _ref + 260 s to 270 s in the TeV light curve are tracking the light-curve bumps in the prompt keV–MeV emission. The close relation between the keV–MeV and TeV emission reveals the continuous energy injection into the external shock. Assuming an energy injection rate exactly following the keV–MeV flux of GRB 221009A, including the very early precursor, we build a continuous energy injection model where the jet Lorentz factor is derived from the TeV time delay, and the TeV data are well fitted, with the TeV excesses interpreted by inverse-Compton scatterings of the inner-coming prompt emission by the energetic electrons in external shock.

We study the gas inflows toward the Galactic Central Molecular Zone (CMZ) based on the gas morphological and kinematic features from the Milky Way Imaging Scroll Painting in the region of l = 1.°2–19.°0 and ∣ b ∣ ≲ 3.°0. We find that the near dust lane appears to extend to l ∼ 15°, in which the end of the large-scale gas structure intersects with the 3 kpc ring at a distance of ∼5 kpc. Intriguingly, many filamentary molecular clouds (MCs), together with the bow-like/ballistic-like clouds and continuous CO features with notable velocity gradient, are finely outlined along the long structure. These MCs also have relatively large velocity dispersions, indicating the shocked gas generated by local continuous accretion and thus the enhanced turbulence along the entire gas structure. We suggest that the ∼3.1–3.6 kpc-long CO structure originates from the accretion molecular gas driven by the Galactic bar. The gas near the bar end at the 3 kpc ring region becomes an important reservoir for the large-scale accreting flows inward to the CMZ through the bar channel. The inclination angle of the bar is estimated to be ϕ _bar = 23° ± 3°, while the pattern speed of the bar is Ω _bar ≲ 32.5 ± 2.5 km s ^−1 kpc ^−1 . The total mass of the whole near gas lane is about 1.3 ± 0.4 × 10 ^7 M _⊙ according to the calculated X _CO ∼ 1.0 ± 0.4 × 10 ^20 cm ^−2 (K km s ^−1 ) ^−1 from the large-scale ^12 CO and ^13 CO data and the complementary H i data. We revisit the gas inflow rate as a mean value of 1.1 ± 0.3 M _⊙ yr ^−1 , which seems to be comparable to the outflow's rate of the Galactic nuclear winds after applying the updated lower X-factor above.

The Imaging X-ray Polarimetry Explorer measured with high significance the X-ray polarization of the brightest Z-source, Sco X-1, resulting in the nominal 2–8 keV energy band in a polarization degree of 1.0% ± 0.2% and a polarization angle of 8° ± 6° at a 90% confidence level. This observation was strictly simultaneous with observations performed by NICER, NuSTAR, and Insight-HXMT, which allowed for a precise characterization of its broadband spectrum from soft to hard X-rays. The source has been observed mainly in its soft state, with short periods of flaring. We also observed low-frequency quasiperiodic oscillations. From a spectropolarimetric analysis, we associate a polarization to the accretion disk at

Investigating the extinction law in regions of high dust extinction, such as the Pillars of Creation within the M16 region, is crucial for understanding the densest parts of the interstellar medium (ISM). In this study, we utilize observations from the Near-Infrared Camera and the Mid-Infrared Instrument onboard the James Webb Space Telescope (JWST) to analyze the color-excess ratios E (F090W − λ )/ E (F090W − F200W) across a wavelength range of 0.9–7.7 μ m. Our method involves performing linear regression on color–color diagrams to derive these ratios. The enhanced detection capabilities of JWST data allow us to probe the distinct extinction law to the densest regions in M16 corresponding to an extinction depth up to A _V ∼ 60 mag. Remarkably, the resultant color-excess ratio curve exhibits a flatter profile than predicted by typical dust extinction models with R _V = 5.5 for dense ISM environments. Moreover, we observe that the mid-infrared extinction law diverges from the near-infrared power law, showing a tendency for the slope to flatten as the wavelength increases. These findings have significant implications for our understanding of the dust properties in dense interstellar environments.

In this work, we report a study on the relationship between flux and intensity for molecular clouds. Our analysis is established on high-quality CO images from the Milky Way Imaging Scroll Painting project. The new flux–intensity relation characterizes the flux variation of molecular clouds above specific intensity levels. We found that the flux–intensity relation exhibits two prominent features. First, the flux–intensity relation generally follows exponential shapes; second, hierarchical structures of molecular clouds are imprinted on flux–intensity relations. Specifically, ^12 CO flux–intensity relations are composed of one or more exponential segments, and for molecular clouds with segmented flux–intensity relations, the edge and the flux of the high-temperature component are strikingly consistent with ^13 CO emission. Further analysis shows that a similar relationship also exists between ^13 CO flux–intensity relations and C ^18 O emission. The mean brightness temperature of molecular clouds is tightly associated with the decay rate of flux, the break temperature of exponential segments, and, to a certain extent, the flux fraction of the high-temperature component. Broadly, the flux–intensity relation of a molecular tracer, either in optically thick or in optically thin cases, has the capability to outline the silhouette of internal structures of molecular clouds, proving to be a potent tool for probing structures of molecular clouds.

Magnetic reconnection plays a crucial role in the energy release process for different kinds of solar eruptions and activities. The rapid solar eruption requires a fast reconnection model. Plasmoid instability in the reconnecting current sheets is one of the most acceptable fast reconnection mechanisms for explaining the explosive events in the magnetohydrodynamics (MHD) scale, which is also a potential bridge between the macroscopic MHD reconnection process and microscale dissipations. Plenty of high-resolution observations indicate that the plasmoid-like structures exist in the high-temperature solar corona, but such evidences are very rare in the lower solar atmosphere with partially ionized plasmas. Utilizing joint observations from the Goode Solar Telescope and the Solar Dynamics Observatory, we discovered a small-scale eruptive phenomenon in NOAA Active Region 13085, characterized by clear reconnection cusp structures, supported by nonlinear force-free field extrapolation results. The plasmoid-like structures with a size of about 150 km were observed to be ejected downward from the current sheet at a maximum velocity of 24 km s ^−1 in the H α line wing images, followed by enhanced emissions at around the postflare loop region in multiple wavelengths. Our 2.5D high-resolution MHD simulations further reproduced such a phenomenon and revealed reconnection fine structures. These results provide comprehensive evidences for the plasmoid-mediated reconnection in partially ionized plasmas, and suggest a unified reconnection model for solar flares with different length scales from the lower chromosphere to the corona.

In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow.

The Event Horizon Telescope observed the horizon-scale synchrotron emission region around the Galactic center supermassive black hole, Sagittarius A* (Sgr A*), in 2017. These observations revealed a bright, thick ring morphology with a diameter of 51.8 ± 2.3 μ as and modest azimuthal brightness asymmetry, consistent with the expected appearance of a black hole with mass M ≈ 4 × 10 ^6 M _⊙ . From these observations, we present the first resolved linear and circular polarimetric images of Sgr A*. The linear polarization images demonstrate that the emission ring is highly polarized, exhibiting a prominent spiral electric vector polarization angle pattern with a peak fractional polarization of ∼40% in the western portion of the ring. The circular polarization images feature a modestly (∼5%–10%) polarized dipole structure along the emission ring, with negative circular polarization in the western region and positive circular polarization in the eastern region, although our methods exhibit stronger disagreement than for linear polarization. We analyze the data using multiple independent imaging and modeling methods, each of which is validated using a standardized suite of synthetic data sets. While the detailed spatial distribution of the linear polarization along the ring remains uncertain owing to the intrinsic variability of the source, the spiraling polarization structure is robust to methodological choices. The degree and orientation of the linear polarization provide stringent constraints for the black hole and its surrounding magnetic fields, which we discuss in an accompanying publication.

The early afterglow of a gamma-ray burst (GRB) can provide critical information on the jet and progenitor of the GRB. The extreme brightness of GRB 221009A allows us to probe its early afterglow in unprecedented detail. In this Letter, we report comprehensive observation results of the early afterglow of GRB 221009A (from T _0 +660 s to T _0 +1860 s, where T _0 is the Insight-HXMT/HE trigger time) in X-ray/gamma-ray energy band (from 20 keV to 20 MeV) by Insight-HXMT High Energy X-ray Telescope, GECAM-C, and Fermi/Gamma-ray Burst Monitor. We find that the spectrum of the early afterglow in 20 keV–20 MeV can be well described by a cutoff power law with an extra power law that dominates the low- and high-energy bands, respectively. The cutoff power law E _peak is ∼30 keV, and the power-law photon index is ∼1.8 throughout the early afterglow phase. By fitting the light curves in different energy bands, we find that a significant achromatic break (from keV to TeV) is required at T _0 + ${1246}_{-26}^{+27}$ s (i.e., 1021 s since the afterglow starting time T _AG = T _0 +225 s), providing compelling evidence of a jet break. Interestingly, both the pre-break and post-break decay slopes vary with energy, and these two slopes become closer in the lower energy band, making the break less identifiable. Intriguingly, the spectrum of the early afterglow experienced a slight hardening before the break and a softening after the break. These results provide new insights into the physics of this remarkable GRB.

Based on the observations by the Parker Solar Probe (PSP) during its encounter phases of approaching the Sun, I. C. Jebaraj et al. found that fundamental–harmonic (F-H) pairs constitute a majority of interplanetary (IP) type III radio bursts. In the present Letter, spectral characteristics of the IP F-H pairs are identified and analyzed further. The observations were made with the Radio Frequency Spectrometer (RFS) experiment on the PSP spacecraft in its encounter phase from the first to the ninth orbit as it traveled from 0.17 to 0.074 au from the Sun. The result shows that the occurrence rate of F-H pairs rises significantly with the rise in the number of IP type III radio bursts detected by the PSP or the enhancement in the time resolution of the RFS instrument. In particular, we compare the relationship between F and H spectral characteristics, such as the frequency-drift rate, emission intensity, relative bandwidth, duration, and fine structure. The results will be helpful for us to understand the physics underlying the generation and evolution of the IP F-H pairs as well as other IP type III radio bursts.

Massive stars mainly form in close binaries, where their mutual interactions can profoundly alter their evolutionary paths. Evolved binaries consisting of a massive OB-type main-sequence star with a stripped helium star or a compact companion represent a crucial stage in the evolution toward double compact objects, whose mergers are (potentially) detectable via gravitational waves. The recent detection of X-ray-quiet OB+black hole binaries and OB+stripped helium star binaries has set the stage for discovering more of these systems in the near future. In this work, based on 3670 detailed binary-evolution models and using empirical distributions of initial binary parameters, we compute the expected population of such evolved massive binaries in coeval stellar populations, including stars in star clusters and in galaxies with starburst activities, for ages up to 100 Myr. Our results are vividly illustrated in an animation that shows the evolution of these binaries in the color–magnitude diagram over time. We find that the number of OB+black hole binaries peaks around 10 Myr, and OB+neutron star binaries are most abundant at approximately 20 Myr. Both black holes and neutron stars can potentially be found in populations with ages up to 90 Myr. Additionally, we analyze the properties of such binaries at specific ages. We find that OB+helium stars and OB+black hole binaries are likely to be identifiable as single-lined spectroscopic binaries. Our research serves as a guide for future observational efforts to discover such binaries in young star clusters and starburst environments.

The heliosphere is permeated with highly structured solar wind originating from the Sun. One of the primary science objectives of Parker Solar Probe (PSP) is to determine the structures and dynamics of the plasma and magnetic fields at the sources of the solar wind. However, establishing the connection between in situ measurements and structures and dynamics in the solar atmosphere is challenging: most of the magnetic footpoint mapping techniques have significant uncertainties in the source localization of a plasma parcel observed in situ, and the PSP plasma measurements suffer from a limited field of view. Therefore, it lacks a universal tool to self-contextualize the in situ measurements. Here we develop a novel time series visualization method named Gaussianity Scalogram. Utilizing this method, by analyzing the magnetic magnitude data from both PSP and Ulysses, we successfully identify in situ structures that are possible remnants of solar atmospheric and magnetic structures spanning more than 7 orders of magnitude, from years to seconds, including polar and midlatitude coronal holes, as well as structures compatible with supergranulation, “jetlets” and “picoflares.” Furthermore, computer simulations of Alfvénic turbulence successfully reproduce the Gaussianization of magnetic magnitude, supporting the observed distribution. Building upon these discoveries, the Gaussianity Scalogram can help future studies to reveal the fractal-like fine structures in the solar wind time series from both PSP and a decades-old data archive.

The inconsistency between experiments in the measurements of the local Universe expansion rate, the Hubble constant, suggests unknown systematics in the existing experiments or new physics. Gravitational-wave standard sirens, a method to independently provide direct measurements of the Hubble constant, have the potential to address this tension. Before that, it is critical to ensure there are no substantial systematics in the standard siren method. A significant systematic has been identified when the viewing angle of the gravitational-wave sources, the compact binary coalescences, was inferred inaccurately from electromagnetic observations of the sources. Such a systematic has led to a more than 10% discrepancy in the standard siren Hubble constant measurements with the observations of binary neutron star merger, GW170817. In this Letter, we develop a new formalism to infer and mitigate this systematic. We demonstrate that the systematic uncertainty of the Hubble constant measurements can be reduced to a level smaller than their statistical uncertainty with 5, 10, and 20 binary neutron star merger observations. We show that our formalism successfully reduces the systematics even if the shape of the biased viewing angle distribution does not follow precisely the model we choose. Our formalism ensures unbiased standard siren Hubble constant measurements when the binary viewing angles are inferred from electromagnetic observations.

The largest geomagnetic storm in two decades occurred in 2024 May with a minimum D _st of −412 nT. We examine its solar and interplanetary origins by combining multipoint imaging and in situ observations. The source active region, NOAA AR 13664, exhibited extraordinary activity and produced successive halo eruptions, which were responsible for two complex ejecta observed at the Earth. In situ measurements from STEREO A, which was 12.°6 apart, allow us to compare the “geo-effectiveness” at the Earth and STEREO A. We obtain key findings concerning the formation of solar superstorms and how mesoscale variations of coronal mass ejections affect geo-effectiveness: (1) the 2024 May storm supports the hypothesis that solar superstorms are “perfect storms” in nature, i.e., a combination of circumstances resulting in an event of an unusual magnitude; (2) the first complex ejecta, which caused the geomagnetic superstorm, shows considerable differences in the magnetic field and associated “geo-effectiveness” between the Earth and STEREO A, despite a mesoscale separation; and (3) two contrasting cases of complex ejecta are found in terms of the geo-effectiveness at the Earth, which is largely due to different magnetic field configurations within the same active region.

We present cosmological constraints from the sample of Type Ia supernovae (SNe Ia) discovered and measured during the full 5 yr of the Dark Energy Survey (DES) SN program. In contrast to most previous cosmological samples, in which SNe are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being an SN Ia, we find 1635 DES SNe in the redshift range 0.10 < z < 1.13 that pass quality selection criteria sufficient to constrain cosmological parameters. This quintuples the number of high-quality z > 0.5 SNe compared to the previous leading compilation of Pantheon+ and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints, we combine the DES SN data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning 0.025 < z < 0.10. Using SN data alone and including systematic uncertainties, we find Ω _M = 0.352 ± 0.017 in flat ΛCDM. SN data alone now require acceleration ( q _0 < 0 in ΛCDM) with over 5 σ confidence. We find $({{\rm{\Omega }}}_{{\rm{M}}},w)=({0.264}_{-0.096}^{+0.074},-{0.80}_{-0.16}^{+0.14})$ in flat w CDM. For flat w _0 w _a CDM, we find $({{\rm{\Omega }}}_{{\rm{M}}},{w}_{0},{w}_{a})=({0.495}_{-0.043}^{+0.033},-{0.36}_{-0.30}^{+0.36},-{8.8}_{-4.5}^{+3.7})$ , consistent with a constant equation of state to within ∼2 σ . Including Planck cosmic microwave background, Sloan Digital Sky Survey baryon acoustic oscillation, and DES 3 × 2pt data gives (Ω _M , w ) = (0.321 ± 0.007, −0.941 ± 0.026). In all cases, dark energy is consistent with a cosmological constant to within ∼2 σ . Systematic errors on cosmological parameters are subdominant compared to statistical errors; these results thus pave the way for future photometrically classified SN analyses.

Some X-ray binaries containing an energetic pulsar in orbit around a normal star accelerate particles to high energies in the shock cone formed where the pulsar and stellar winds collide. The magnetic field geometry in the acceleration region in such binaries is unknown. We performed the first measurement of the polarization of the X-ray synchrotron emission from a gamma-ray-emitting binary system. We observed PSR B1259–63 with the Imaging X-ray Polarimetry Explorer during an X-ray bright phase following the periastron passage in 2024 June. X-ray polarization is detected with a polarization degree of 8.3% ± 1.5% at a significance of 5.3 σ . The X-ray polarization angle is aligned with the axis of the shock cone at the time of the observation. This indicates that the predominant component of the magnetic field in the acceleration region is oriented perpendicular to the shock-cone axis.

Fast radio bursts (FRBs) are millisecond-duration flashes with unknown origins. Their formation rate is crucial for unveiling physical origins. However, the luminosity and formation rate are degenerate when directly fitting the redshift distribution of FRBs. In contrast to previous forward-fitting methods, we use Lynden-Bell’s c ^− method to derive the luminosity function and formation rate of FRBs without any assumptions. Using the nonrepeating FRBs from the first Canadian Hydrogen Intensity Mapping Experiment FRB catalog, we find a relatively strong luminosity evolution, and luminosity function can be fitted by a broken power-law model with a break at 1.33 × 10 ^41 erg s ^−1 . The formation rate declines rapidly as (1 + z ) ^−4.9±0.3 with a local rate of 1.13 × 10 ^4 Gpc ^−3 yr ^−1 . This monotonic decrease is similar to the rate of short gamma-ray bursts. After comparing this function with the star formation rate and stellar mass density, we conclude that the old populations, including neutron stars and black holes, are closely related to the origins of FRBs.

This Letter reports the first observation of the onset of fully developed turbulence in the solar corona. Long time series of white-light coronal images, acquired by Metis aboard Solar Orbiter at 2 minutes cadence and spanning about 10 hr, were studied to gain insight into the statistical properties of fluctuations in the density of the coronal plasma in the time domain. From pixel-by-pixel spectral frequency analysis in the whole Metis field of view, the scaling exponents of plasma fluctuations were derived. The results show that, over timescales ranging from 1 to 10 hr and corresponding to the photospheric mesogranulation-driven dynamics, the density spectra become shallower moving away from the Sun, resembling a Kolmogorov-like spectrum at 3 R _⊙ . According to the latest observation and interpretive work, the observed 5/3 scaling law for density fluctuations is indicative of the onset of fully developed turbulence in the corona. Metis observation-based evidence for a Kolmogorov turbulent form of the fluctuating density spectrum casts light on the evolution of 2D turbulence in the early stages of its upward transport from the low corona.

Alpha–proton differential flow ( V _α _p ) of coronal mass ejections (CMEs) and solar wind from the Sun to 1 au and beyond could influence the instantaneous correspondence of absolute abundances of alpha particles (He ^2+ /H ^+ ) between the solar corona and interplanetary space as the abundance of a coronal source can vary with time. Previous studies based on Ulysses and Helios showed that V _α _p is negligible within CMEs from 5 to 0.3 au, similar to slow solar wind (

The systematic properties are largely unknown for the black hole X-ray binary Swift J151857.0–572147 newly discovered in the 2024 outburst. The nature of a black hole can be completely defined by specifying the mass and dimensionless spin parameter. Therefore, accurate measurement of the two fundamental parameters is important for understanding the nature of black holes. The joint spectral fitting of a reflection component with simultaneous observations from Insight-HXMT, NICER, and NuSTAR reveals for the first time a black hole dimensionless spin of ${0.84}_{-0.26}^{+0.17}$ and an inclination angle of ${21.1}_{-3.6}^{+4.5}$ degrees for this system. Monitoring of the soft state by NICER results in disk flux and temperature following ${F}_{\mathrm{disk}}\propto {T}_{\mathrm{in}}^{3.83\pm 0.17}$ . For the standard thin disk, ${L}_{\mathrm{disk}}\approx 4\pi {R}_{\mathrm{in}}^{2}\sigma {T}_{\mathrm{in}}^{4}$ , so the relationship between the flux and temperature of the disk we measured indicates that the inner radius of the disk is stable and the disk is in the innermost stable circular orbit. With an empirical relation built previously between the black hole outburst profile and the intrinsic power output, the source distance is estimated as 5.8 ± 2.5 kpc according to the outburst profile and peak flux observed by Insight-HXMT and NICER. Finally, a black hole mass of 3.67 ± 1.79–8.07 ± 4.20 M _⊙ can be inferred from a joint diagnostic of the aforementioned parameters measured for this system. This system is also consistent with most black hole X-ray binaries with high spin and a mass in the range of 5–20 M _⊙ .

The periodic coronal rain and in-phase radiative intensity pulsations have been observed in multiple wavelengths in recent years. However, due to the lack of three-dimensional coronal magnetic fields and thermodynamic data in observations, it remains challenging to quantify the coronal heating rate that drives the mass cycles. In this work, based on the MURaM code, we conduct a three-dimensional radiative magnetohydrodynamic simulation spanning from the convective zone to the corona, where the solar atmosphere is heated self-consistently through dissipation resulting from magnetoconvection. For the first time, we model the periodic coronal rain in an active region. With a high spatial resolution, the simulation well resembles the observational features across different extreme-ultraviolet wavelengths. These include the realistic interweaving coronal loops, periodic coronal rain, and periodic intensity pulsations, with two periods of 3.0 hr and 3.7 hr identified within one loop system. Moreover, the simulation allows for a detailed three-dimensional depiction of coronal rain on small scales, revealing adjacent shower-like rain clumps ∼500 km in width and showcasing their multithermal internal structures. We further reveal that these periodic variations essentially reflect the cyclic energy evolution of the coronal loop under thermal nonequilibrium state. Importantly, as the driver of the mass circulation, the self-consistent coronal heating rate is considerably complex in time and space, with hour-level variations in 1 order of magnitude, minute-level bursts, and varying asymmetry reaching ten times between footpoints. This provides an instructive template for the ad hoc heating function and further enhances our understanding of the coronal heating process.

We present James Webb Space Telescope (JWST) NIRCam and NIRSpec observations of a Type Ic supernova (SN Ic) and its host galaxy (JADES-GS+53.13533-27.81457) at z = 2.83. This SN (named SN 2023adta) was identified in deep JWST/NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) program. Follow-up observations with JWST/NIRSpec provided a spectroscopic redshift of z = 2.83 and the classification as an SN Ic-BL. The light curve of SN 2023adta matches well with other stripped-envelope SNe, and we find a high peak luminosity, M _V = −19.0 ± 0.2 mag, based on the distribution of best-fit SNe. The broad absorption features in its spectrum are consistent with other SNe Ic-BL 1–3 weeks after peak brightness. We measure a Ca ii near-IR triplet expansion velocity of 29,000 ± 2000 km s ^−1 . The host galaxy of SN 2023adta is irregular, and modeling of its spectral energy distribution indicates a metallicity of $Z={0.35}_{-0.08}^{+0.16}{Z}_{\odot }$ . This environment is consistent with the population of low- z SNe Ic-BL that prefer lower metallicities relative to other stripped-envelope SNe and track long-duration γ -ray burst environments. We do not identify any γ -ray bursts that are coincident with SN 2023adta. Given the rarity of SNe Ic-BL in the local Universe, the detection of an SN Ic-BL at z = 2.83 could indicate that their rates are enhanced at high redshift.

The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high angular resolution (∼0.″3, 870 au) map of the magnetic field projected on the plane of the sky ( B _POS ) toward the high-mass star-forming region G333.46−0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array at 1.2 mm as part of the Magnetic fields in Massive star-forming Regions survey. The B _POS morphology found in this region is consistent with a canonical “hourglass” with an embedded flattened envelope in a perpendicular direction, which suggests a dynamically important field. This region is fragmented into two protostars that appear to be gravitationally bound in a stable binary system with a separation of ∼1740 au. Interestingly, by analyzing H ^13 CO ^+ ( J = 3–2) line emission, we find no velocity gradient over the extent of the continuum, which is consistent with a strong field. We model the B _POS , obtaining a marginally supercritical mass-to-flux ratio of 1.43, suggesting an initially strongly magnetized environment. Based on the Davis–Chandrasekhar–Fermi method, the magnetic field strength toward G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows toward the central region of G333 suggests strong magnetic braking, consistent with a highly magnetized environment. Our study shows that despite being a strong regulator, the magnetic energy fails to prevent the process of fragmentation, as revealed by the formation of the two protostars in the central region.

Super active region NOAA 13664 produced 12 X-class flares (including the largest one so far, an occulted X8.7 flare, in solar cycle 25) during 2024 May 8–15, and 11 of them are identified as white-light flares. Here we present various features of these X-class white-light flares observed by the White-light Solar Telescope (WST) on board the Advanced Space-based Solar Observatory and the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory. It is found that both the white-light emissions at WST 3600 Å (Balmer continuum) and HMI 6173 Å (Paschen continuum) show up in different regions of the sunspot group in these flares, including outside the sunspots and within the penumbra and umbra of the sunspots. They exhibit a point-, ribbon-, loop-, or ejecta-like shape, which can come from flare ribbons (or footpoints), flare loops, and plasma ejecta depending on the perspective view. The white-light duration and relative enhancement are measured and both parameters for 3600 Å emission have greater values than those for 6173 Å emission. It is also found that these white-light emissions are cospatial well with the hard X-ray (HXR) sources in the on-disk flares but have some offsets with the HXR emissions in the off-limb flares. In addition, it is interesting that the 3600 and 6173 Å emissions show different correlations with the peak HXR fluxes, with the former one more sensitive to the HXR emission. All these greatly help us understand the white-light flares of a large magnitude from a super active region on the Sun and also provide important insights into superflares on Sun-like stars.

Based on observations from the Parker Solar Probe in the near-Sun solar wind, this study identifies an ion-scale wave event characterized by two distinct frequency bands. The lower-band waves exhibit right-hand polarization, while the upper-band waves have left-hand polarization. Alongside these waves, there are clear indications of the existence of both proton core and beam components, with the perpendicular temperature being higher than the parallel temperature in the measured proton velocity distribution functions (VDFs). Utilizing the plasma parameters derived from typical proton VDFs, instability analyses are conducted to investigate the mode nature of the observed waves and their excitation mechanism. The lower-band waves are identified as sunward ion cyclotron waves (ICWs), generated through the proton beam cyclotron instability; the upper-band waves are recognized as antisunward ICWs, induced by the proton core cyclotron instability. This study provides the first direct observational evidence confirming the presence of counterpropagating ICWs and proton cyclotron instability in the solar wind.

We present the James Webb Space Telescope (JWST) discovery of SN 2023adsy, a transient object located in a host galaxy JADES-GS+53.13485−27.82088 with a host spectroscopic redshift of 2.903 ± 0.007. The transient was identified in deep (JWST)/NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) program. Photometric and spectroscopic follow-up with NIRCam and NIRSpec, respectively, confirm the redshift and yield UV-NIR light-curve, NIR color, and spectroscopic information all consistent with a Type Ia classification. Despite its classification as a likely SN Ia, SN 2023adsy is both fairly red ( c ∼ 0.9) despite a host galaxy with low extinction and has a high Ca ii velocity (19,000 ± 2000 km s ^−1 ) compared to the general population of SNe Ia. While these characteristics are consistent with some Ca-rich SNe Ia, particularly SN 2016hnk, SN 2023adsy is intrinsically brighter than the low- z Ca-rich population. Although such an object is too red for any low- z cosmological sample, we apply a fiducial standardization approach to SN 2023adsy and find that the SN 2023adsy luminosity distance measurement is in excellent agreement (≲1 σ ) with ΛCDM. Therefore unlike low- z Ca-rich SNe Ia, SN 2023adsy is standardizable and gives no indication that SN Ia standardized luminosities change significantly with redshift. A larger sample of distant SNe Ia is required to determine if SN Ia population characteristics at high z truly diverge from their low- z counterparts and to confirm that standardized luminosities nevertheless remain constant with redshift.

The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. The timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR J0437–4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition Explorer (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ( M _p / R ) of the star. Direct measurements of the mass M _p from pulsar timing therefore substantially improve constraints on the radius R and the equation of state. Here we use observations spanning 26 yr from Murriyang, the 64 m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass M _p = 1.418 ± 0.044 M _⊙ , distance D = 156.96 ± 0.11 pc, and orbital inclination angle i = 137.°506 ± 0.°016, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modeled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar or excess low-frequency noise such as a gravitational-wave background.

The gas-phase metallicity of galaxies is regulated by multiple astrophysical processes, which makes it a crucial diagnostic of galaxy formation and evolution. Beyond the fundamental mass–metallicity relation, a debate about the secondary galaxy property to predict the metallicity of galaxies arises. Motivated by this, we systematically examine the relationship between gas-phase metallicity and other galaxy properties, i.e., the star formation rate (SFR) and galaxy size, in addition to stellar mass in both observation and simulation. We utilize the data from the Mapping Nearby Galaxies at Apache Point Observatory survey and the TNG50 simulations. We find that the combination of ${M}_{* }/{R}_{{\rm{e}}}^{\beta }$ with β ∼ 0.6–1 is in much stronger correlation to the metallicity than stellar mass alone, regardless of whether the SFR is included or not, in both observation and simulation. This indicates that galaxy size plays a more important role in determining gas-phase metallicity of galaxies than SFR. In addition, The Next Generation simulation predicts that the SFR, although being a subdominant role, becomes increasingly important in the high- z universe. Finally, we speculate that the SFR modulates metallicity on the temporal dimension, synchronized with time-varying gas inflows, and galaxy size regulates metallicity on the spatial dimension by affecting the gravitational potential and the mass-loading factor.

To date, the dayside thermal structure of ultra-hot Jupiters (UHJs) is generally considered to be inverted, but their nightside thermal structure has been less explored. Here, we explore the impact of nightside thermal emission on high-resolution infrared transmission spectroscopy, which should not be neglected, especially for UHJs. We present a general equation for the high-resolution transmission spectrum that includes planetary nightside thermal emission. This provides a new way to infer the thermal structure of the planetary nightside with high-resolution transmission spectroscopy. Using the cross-correlation technique, we find evidence for the presence of an H _2 O emission signature on the UHJ WASP-33 b during the transit, indicating an inverted temperature structure on its nightside. Such a result suggests a stronger heat transport through the circulation than currently expected. An alternative explanation is that the rotating visible hemisphere during transit leads to the potential contribution of the limb and dayside atmospheres to the detected emission signature. In the future, the combination of high-resolution, full-phase-curve spectroscopic observations and general circulation models will hopefully solve this puzzle and provide a complete picture of the three-dimensional nature of the chemistry, circulation, and thermal structure of UHJs.

We present early-phase good-cadence (hour-to-day) simultaneous multiband ( ugi and vrz bands) imaging of the nearby supernova SN 2024ggi, which exploded in the nearby galaxy NGC 3621. A quick follow-up was conducted within less than a day after the explosion and continued for ∼23 days. The uvg band light curves display a rapid rise (∼1.4 mag day ^−1 ) to maximum in ∼4 days and absolute magnitude M _g ∼ −17.75 mag. The postpeak decay rate in redder bands is ∼0.01 mag day ^−1 . Different colors (e.g., u − g and v − r ) of SN 2024ggi are slightly redder than SN 2023ixf. A significant rise (∼12.5 kK) in blackbody temperature (optical) was noticed within ∼2 days after the explosion, which successively decreased, indicating shock break out inside a dense circumstellar medium surrounding the progenitor. Using semianalytical modeling, the ejecta mass and progenitor radius were estimated as 1.2 M _☉ and ∼550 R _☉ . The archival deep images ( g , r , i , and z bands) from the Dark Energy Camera Legacy Survey were examined, and a possible progenitor was detected in each band (∼22–22.5 mag) and had a mass range of 14–17 M _☉ .

Wolf–Rayet stars in close binary systems can be tidally spun up by their companions, potentially leaving behind fast-spinning, highly magnetized neutron stars, known as “magnetars,” after core collapse. These newborn magnetars can transfer rotational energy into heating and accelerating the ejecta, producing hydrogen-poor superluminous supernovae (SLSNe). In this Letter, we propose that the magnetar wind of the newborn magnetar could significantly evaporate its companion star, typically a main-sequence or helium star, if the binary system is not disrupted by the abrupt mass loss and supernova (SN) kick. The subsequent heating and acceleration of the evaporated star material along with the SN ejecta by the magnetar wind can produce a postpeak bump in the SLSN lightcurve. Our model can reproduce the primary peaks and postpeak bumps of four example observed multiband SLSN lightcurves, revealing that the mass of the evaporated material could be ∼0.4–0.6 M _⊙ if the material is hydrogen-rich. We propose that the magnetar could induce strongly enhanced evaporation from its companion star near the pericenter if the orbit of the post-SN binary is highly eccentric, ultimately generating multiple postpeak bumps in the SLSN lightcurves. This “magnetar–star binary engine” model may offer a possible explanation for the evolution of polarization, along with the origin and velocity broadening of late-time hydrogen or helium broad spectral features observed in some bumpy SLSNe. The diversity in the lightcurves and spectra of SLSNe may be attributed to the wide variety of companion stars and post-SN binary systems.

Highly r -process-enhanced (RPE) stars are rare and usually metal poor ([Fe/H] < −1.0), and they mainly populate the Milky Way halo and dwarf galaxies. This study presents the discovery of a relatively bright ( V = 12.72), highly RPE ( r -II) star ([Eu/Fe] = +1.32, [Ba/Eu] = −0.95), LAMOST J020623.21+494127.9. This star was selected from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope medium-resolution ( R ∼ 7500) spectroscopic survey; follow-up high-resolution ( R ∼ 25,000) observations were conducted with the High Optical Resolution Spectrograph installed on the Gran Telescopio Canarias. The stellar parameters ( T _eff = 4130 K, $\mathrm{log}\,{\rm{g}}$ = 1.52, [Fe/H] = −0.54, ξ = 1.80 km s ^−1 ) have been inferred taking into account nonlocal thermodynamic equilibrium effects. The abundances of [Ce/Fe], [Pr/Fe], and [Nd/Fe] are +0.19, +0.65, and +0.64, respectively, relatively low compared to the Solar r -process pattern normalized to Eu. This star has a high metallicity ([Fe/H] = −0.54) compared to most other highly RPE stars and has the highest measured abundance ratio of Eu to H ([Eu/H] = +0.78). It is classified as a thin-disk star based on its kinematics and does not appear to belong to any known stream or dwarf galaxy.

Mrk 501 is a prototypical high-synchrotron-peaked blazar and serves as one of the primary targets for the Imaging X-ray Polarimetry Explorer (IXPE). In this study, we report X-ray polarization measurements of Mrk 501 based on six IXPE observations. The detection of X-ray polarization at a confidence level exceeding 99% is achieved in four out of the six observations conducted across the entire energy range (2–8 keV) of IXPE. The maximum polarization degree (Π _X ) is measured to be 15.8% ± 2.8%, accompanied by a polarization angle ( ψ _X ) of 98.°0 ± 5.°1 at a confidence level of 5.6 σ . During the remaining two observations, only an upper limit of Π _X < 12% could be derived at the 99% confidence level. No temporal variability in polarization is observed throughout all six IXPE observations for Mrk 501. A discernible trend of energy-dependent variation in the polarization degree is detected in optical spectropolarimetry; however, no analogous indication is observed in Π _X . The chromatic behavior of Π and the consistent values of ψ across different frequencies from X-ray to radio bands, along with the agreement between ψ and jet position angle, strongly support the interpretation of the energy-stratified model with shock-accelerated particles in the jet of Mrk 501. Additionally, the possibility of the presence of a global helical magnetic field in the jet of Mrk 501 is discussed.

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce very high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy γ -ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the γ -ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The γ -ray luminosity above 0.5 TeV of both sources is (2–3) × 10 ^35 erg s ^−1 . This exceeds by more than a factor of 2 the luminosity of HESS J1646−458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the γ -ray emission from each source is extended with a significance of >3 σ and a Gaussian width of about 30 pc. For 30 Dor C, a connection between the γ -ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the γ -ray signal from R136 are discussed.

Type IIb supernovae are an important subclass of stripped-envelope supernovae (SNe), which show H lines only at early times. Their progenitors are believed to contain a low-mass H envelope before explosion. This work reports the discovery of a progenitor candidate in preexplosion Hubble Space Telescope images for the Type IIb SN 2017gkk. With detailed analysis of its spectral energy distribution and local environment, we suggest that the progenitor is most likely a yellow supergiant with significant circumstellar extinction and has an initial mass of about 16 M _⊙ , effective temperature log( T _eff /K) = 3.72 ± 0.08, and luminosity log( L / L _⊙ ) = 5.17 ± 0.04. This progenitor is not massive enough to strip envelope through stellar wind, and it supports an interacting binary progenitor channel and adds to the growing list of direct progenitor detections for Type IIb SNe. Future late-time observations will confirm whether this progenitor candidate has disappeared and reveal the putative binary companion that has survived the explosion.

We present the first X-ray spectropolarimetric results for Cygnus X-1 in its soft state from a campaign of five IXPE observations conducted during 2023 May–June. Companion multiwavelength data during the campaign are likewise shown. The 2–8 keV X-rays exhibit a net polarization degree PD = 1.99% ± 0.13% (68% confidence). The polarization signal is found to increase with energy across the Imaging X-ray Polarimetry Explorer’s (IXPE) 2–8 keV bandpass. The polarized X-rays exhibit an energy-independent polarization angle of PA = −25.°7 ± 1.°8 east of north (68% confidence). This is consistent with being aligned to Cyg X-1’s au-scale compact radio jet and its parsec-scale radio lobes. In comparison to earlier hard-state observations, the soft state exhibits a factor of 2 lower polarization degree but a similar trend with energy and a similar (also energy-independent) position angle. When scaling by the natural unit of the disk temperature, we find the appearance of a consistent trend line in the polarization degree between the soft and hard states. Our favored polarimetric model indicates that Cyg X-1’s spin is likely high ( a _* ≳ 0.96). The substantial X-ray polarization in Cyg X-1's soft state is most readily explained as resulting from a large portion of X-rays emitted from the disk returning and reflecting off the disk surface, generating a high polarization degree and a polarization direction parallel to the black hole spin axis and radio jet. In IXPE’s bandpass, the polarization signal is dominated by the returning reflection emission. This constitutes polarimetric evidence for strong gravitational lensing of X-rays close to the black hole.

We present a detailed analysis of the progenitor and its local environment for the recently discovered Type II supernova (SN) 2024ggi at a distance of about 6.7 Mpc, by utilizing the pre-explosion images from the Hubble Space Telescope and Spitzer Space Telescope. The progenitor is identified as a red bright variable star, with absolute F814W-band magnitudes being −6.2 mag in 1995 to −7.2 mag in 2003, respectively, consistent with that of a normal red supergiant star. Combining with the historical mid-infrared light curves, a pulsational period of about 379 days can be inferred for the progenitor star. Fitting its spectral energy distribution with stellar spectral models yields the stellar parameters of temperature, radius, and bolometric luminosity as ${T}_{* }={3290}_{-27}^{+19}$ K, ${R}_{* }={887}_{-51}^{+60}$ R _⊙ , and log( L / L _⊙ ) $\,=\,{4.92}_{-0.04}^{+0.05}$ , respectively. The above parameters indicate that the progenitor of SN 2024ggi is consistent with the stellar evolutionary track of a solar-metallicity massive star with an initial mass of ${13}_{-1}^{+1}$ M _⊙ . Moreover, our analysis indicates a relatively low mass-loss rate (i.e.,

Fast X-ray Transients (FXTs) are extragalactic bursts of soft X-rays first identified ≳10 yr ago. Since then, nearly 40 events have been discovered, although almost all of these have been recovered from archival Chandra and XMM-Newton data. To date, optical sky surveys and follow-up searches have not revealed any multiwavelength counterparts. The Einstein Probe, launched in 2024 January, has started surveying the sky in the soft X-ray regime (0.5–4 keV) and will rapidly increase the sample of FXTs discovered in real time. Here we report the first discovery of both an optical and radio counterpart to a distant FXT, the fourth source publicly released by the Einstein Probe. We discovered a fast-fading optical transient within the 3′ localization radius of EP 240315a with the all-sky optical survey ATLAS, and our follow-up Gemini spectrum provides a redshift, z = 4.859 ± 0.002. Furthermore, we uncovered a radio counterpart in the S band (3.0 GHz) with the MeerKAT radio interferometer. The optical (rest-frame UV) and radio luminosities indicate that the FXT most likely originates from either a long gamma-ray burst or a relativistic tidal disruption event. This may be a fortuitous early mission detection by the Einstein Probe or may signpost a mode of discovery for high-redshift, high-energy transients through soft X-ray surveys, combined with locating multiwavelength counterparts.

On 2022 March 28 two successive coronal mass ejections (CMEs) were observed by multiple spacecraft and resulted in a magnetic cloud (MC) at 1 au. We investigate the propagation and interaction properties of the two CMEs correlated with the MC using coordinated multipoint remote sensing and in situ observations from Solar Orbiter, STEREO-A, Solar and Heliospheric Observatory, and Wind. The first CME was triggered by a filament eruption with a high inclination angle. Roughly 9 hr later, the second CME originating from the same active region erupted with a smaller tilt angle and faster speed compared to the first one. The second CME overtook the preceding CME and formed a merged front at approximately 75 R _☉ , which developed into a complex ejecta at 1 au. The descending speed and low proton temperature inside the complex ejecta suggest that the two CMEs have fully merged before reaching 1 au, leading them to begin expanding rather than compressing against each other. The complex ejecta appears to have the magnetic field and plasma signatures of an MC, although there is a discontinuity in the magnetic field implying previous interactions. The cross section of the complex ejecta, reconstructed from in situ data using a Grad–Shafranov technique, exhibits a right-handed flux rope structure. These results highlight that an MC-like complex ejecta lacking interaction features could arise from the complete merging of two CMEs.

We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) band 6 data for the outbursting massive protostar M17 MIR. The ALMA CO J = 2–1 data reveal a collimated and bipolar north–south outflow from M17 MIR. The blueshifted outflow exhibits four CO knots (N1 to N4) along the outflow axis, while the redshifted outflow appears as a single knot (S1). The extremely high velocity (EHV) emissions of N1 and S1 are jetlike and contain subknots along the outflow axis. Assuming the nearest EHV subknots trace the ejecta from the accretion outbursts in the past decades, a tangential ejection velocity of ∼421 km s ^−1 is derived for M17 MIR. Assuming the same velocity, the dynamical times of the multiple ejecta, traced by the four blueshifted CO knots, range from 20 to 364 yr. The four blueshifted CO knots imply four clustered accretion outbursts with a duration of tens of years in the past few hundred years. The intervals between the four clustered accretion outbursts are also about tens of years. These properties of the four clustered accretion outbursts are in line with the disk gravitational instability and fragmentation model. The episodic accretion history of M17 MIR traced by episodic outflow suggests that a massive star can form from a lower-mass protostar via frequent episodic accretion events triggered by disk gravitational instability and fragmentation. The first detection of the knotty outflow from an outbursting massive protostar suggests that mass ejections accompanied with accretion events could serve as an effective diagnostic tool for the episodic accretion histories of massive protostars.