Multi-wavelength Observations Research Articles

Multiwavelength Campaign Observations of a Young Solar-type Star, EK Draconis. I. Discovery of Prominence Eruptions Associated with Superflares

Young solar-type stars frequently produce superflares, serving as a unique window into the young Sun-Earth environments. Large solar flares are closely linked to coronal mass ejections (CMEs) associated with filament/prominence eruptions, but observational evidence for stellar superflares remains scarce. Here, we present a 12-day, multiwavelength campaign observation of young solar-type star EK Draconis (G1.5V, 50–120 Myr age) utilizing the Transiting Exoplanet Survey Satellite, the Neutron star Interior Composition ExploreR, and the Seimei telescope. The star has previously exhibited blueshifted Hα absorptions as evidence for a filament eruption associated with a superflare. Our simultaneous optical and X-ray observations identified three superflares of 1.5 × 1033–1.2 × 1034 erg. We report the first discovery of two prominence eruptions on a solar-type star, observed as blueshifted Hα emissions at speeds of 690 and 430 km s−1 and masses of 1.1 × 1019 and 3.2 × 1017 g, respectively. The faster, massive event shows a candidate of post-flare X-ray dimming with the amplitude of up to ∼10%. Several observational aspects consistently point to the occurrence of a fast CME associated with this event. The comparative analysis of the estimated length scales of flare loops, prominences, possible dimming region, and starspots provides the overall picture of the eruptive phenomena. Furthermore, the energy partition of the observed superflares in the optical and X-ray bands is consistent with flares from the Sun, M-dwarfs, and close binaries, yielding the unified empirical relations. These discoveries provide profound implications of the impact of these eruptive events on early Venus, Earth, and Mars and young exoplanets.

Onset and Evolution of Solar Flares: Application of 2D and 3D Models of Magnetic Reconnection

Contemporary multi-wavelength observations have revealed various important features during solar flares. On one hand, these observations support the two-dimensional (2D) “standard flare model,” while on other hand, they also call for exploration of three-dimensional (3D) magnetic field topologies involved in flares. Traditionally, the formation of parallel ribbons on both sides of the polarity inversion line (PIL) and associated overlying loop arcades have been recognized as the most prominent features of eruptive flares, forming the basis for the development of the standard model that provides a 2D description of flare-associated phenomena. However, the actual flare occurs within a more complicated 3D magnetic structure. Thus, despite its general applicability, the standard model has limited or no scope in explaining some features that exclusively require a 3D description. In this context, the observations of “circular ribbon flares” stand out, where one of the ribbons presents an almost fully closed quasi-circular or quasi-ellipsoidal shape, providing evidence of the involvement of a typical fan-spine magnetic configuration. In this article, we discuss observational features vis-à-vis theoretical understanding of solar flares in view of 2D and 3D models of magnetic reconnection. We highlight a few complex cases of circular ribbon flares exhibiting parallel ribbons, a coronal jet, and/or an erupting magnetic flux rope. Exploring various 3D topologies also enables us to probe similarities between the circumstances that govern the onset of jets, confined flares and CME-producing eruptive flares.

Exploring Magnetic Coupling of Solar Atmosphere Through Frequency Modulations of 3-min Slow Magnetoacoustic Waves

Coronal fan loops rooted in sunspot umbra show outward propagating waves with subsonic phase speed and a period around 3 min. However, their source region in the lower atmosphere is still ambiguous. We performed multi-wavelength observations of a clean fan loop system rooted in a sunspot utilizing data from the Interface Region Imaging Spectrograph (IRIS) and the Solar Dynamics Observatory (SDO). We explored the less studied property of frequency modulation of these 3-min waves from the photosphere to the corona, and identified periodicities in the ranges of 14-20 min, and 24-35 min. Based on our findings, we interpret that the 3-min slow waves observed in the coronal fan loops are driven by 3-min oscillations observed at the photospheric footpoints of these fan loops in the umbral region. We also investigated possible connections between the 3-min and 5-min oscillations observed at the photosphere, which are still poorly understood. Our results provide clear evidence of magnetic coupling of the solar umbral atmosphere through the propagation of 3-min waves along the fan loops at different atmospheric heights.

Multiwavelength observation of 1A 0535+262=HD 245770 from 2010 to 2021

Context. 1A 0535+262 is a high-mass X-ray binary that went into a giant X-ray outburst in 2020. During this event, the X-ray luminosity reached the highest value measured over the last 30 years. Aims. Our aim is to study the long-term variability of 1A 0535+262 before and after the 2020 major X-ray outburst and to uncover the mechanism that led to the X-ray outburst. Methods. We used the long-term photometric light curve and the equivalent widths of the Hα and He I λ6678 lines to monitor the state of the Be star’s circumstellar disk. The Hα line profiles show evidence for V/R variability, which we revealed by fitting the Hα spectral line profiles with two Gaussian functions. In addition, we divided our data into four periods according to the intensity of the X-ray, optical, and infrared emission. Results. The Hα line profiles show single-peaked profiles in most cases. This is consistent with the previously reported orbital inclination of i = 37° ±2°. Unlike the Hα lines, the He I λ6678 lines show a maximal intensity in October 2020, which is one month before the giant X-ray outburst in 2020. Based on the behavior of the equivalent widths of the Hα and He I λ6678 lines, and the V-band magnitude, we find two mass ejection processes from the Be star to the Be disk on MJD 55820 and MJD 56600. The V/R quasi-period is about two years during 2011–2015, which is different from 1994 to 1995. Furthermore, the periods I → II → III → IV in the (B − V) color index versus V-band magnitude diagram constitute a cycle. From the behavior of the V/R ratio of Hα lines, and the variability of the V band, we believe that the precession of the density perturbation inside the disk is retrograde.

Multiline Stokes Synthesis of Ellerman Bombs: Obtaining Seamless Information from Photosphere to Chromosphere

Magnetic reconnection in the lower atmosphere is a critical process in determining the chromospheric dynamics, such as Ellerman bombs and UV bursts. Because the heating of the atmosphere significantly depends on the ionization degree and plasma β, which varies with height, it is essential to diagnose the height at which the magnetic reconnection takes place. Multiwavelength spectropolarimetry is a powerful solution to fulfill this requirement. We verify the diagnostic capabilities and usefulness of near-infrared multiwavelength spectropolarimetric observations for understanding magnetic reconnection phenomena by synthesizing the Stokes vector from a realistic magnetohydrodynamic simulation. The analysis considers two magnetic reconnection regions occurring at different heights. In the case of magnetic reconnection at low altitude, both red- and blueshifted components originating from reconnection bidirectional flow are identified in the photospheric lines, Fe i 8468 Å, K i 7664 Å, and K i 7698 Å. In the case of magnetic reconnection at high altitudes, chromospheric lines, Ca ii 8498 Å and 8542 Å, show emission due to the heating that occurs at the upper part of the formation layer. These results suggest that multiwavelength spectropolarimetric observations are capable of distinguishing the height where magnetic reconnection occurs.

A shock flash breaking out of a dusty red supergiant.

Shock-breakout emission is light that arises when a shockwave, generated by the core-collapse explosion of a massive star, passes through its outer envelope. Hitherto, the earliest detection of such a signal was at several hours after the explosion1, although a few others had been reported2-7. The temporal evolution of early light curves should provide insights into the shock propagation, including explosion asymmetry and environment in the vicinity, but this has been hampered by the lack of multiwavelength observations. Here we report the instant multiband observations of a type II supernova (SN 2023ixf) in the galaxy M101 (at a distance of 6.85 ± 0.15 Mpc; ref. 8), beginning at about 1.4 h after the explosion. The exploding star was a red supergiant with a radius of about 440 solar radii. The light curves evolved rapidly, on timescales of 1-2 h, and appeared unusually fainter and redder than predicted by the models9-11 within the first few hours, which we attribute to an optically thick dust shell before it was disrupted by the shockwave. We infer that the breakout and perhaps the distribution of the surrounding dust were not spherically symmetric.

The Triggering Process of an X-class Solar Flare on a Small Quadrupolar Active Region

The occurrence of X-class solar flares and their potential impact on space weather often receive greater attention than other flares. But predicting when and where an X-class flare will occur is still a challenge. With the multiwavelength observation from the Solar Dynamics Observatory and FengYun-3E satellite, we investigate the triggering of a GOES X1.0 flare occurring in the NOAA active region (AR) 12887. Our results show that this unique X-class flare is bred in a relatively small but complex quadrupolar AR. Before the X-class flare, two filaments (F1 and F2) exist below a null-point topology of the quadrupolar AR. Magnetic field extrapolation and observation reveal that F1 and F2 correspond to two magnetic flux ropes with the same chirality and that their adjacent feet are respectively rooted at nonconjugated opposite polarities. Interestingly, these two polarities collide rapidly, accompanied by photospheric magnetic flux emergence, cancellation, and shear motion in the AR center. Above this site, F1 and F2 subsequently intersect with each other and merge into a longer filament (F3) via a tether-cutting-like reconnection process. As a result, F3 rises and erupts upward, eventually leading to a coronal mass ejection and the X-class flare. This observation suggests that the rapid collision of nonconjugated opposite polarities provides a key condition for the triggering of this X-class flare and also provides a featured case for a flare trigger mechanism and space weather forecasting.

Establishing the impact of luminous AGN with multi-wavelength observations and simulations

AbstractCosmological simulations fail to reproduce realistic galaxy populations without energy injection from active galactic nuclei (AGN) into the interstellar medium (ISM) and circumgalactic medium (CGM); a process called ‘AGN feedback’. Consequently, observational work searches for evidence that luminous AGN impact their host galaxies. Here, we review some of this work. Multi-phase AGN outflows are common, some with potential for significant impact. Additionally, multiple feedback channels can be observed simultaneously; e.g., radio jets from ‘radio quiet’ quasars can inject turbulence on ISM scales, and displace CGM-scale molecular gas. However, caution must be taken comparing outflows to simulations (e.g., kinetic coupling efficiencies) to infer feedback potential, due to a lack of comparable predictions. Furthermore, some work claims limited evidence for feedback because AGN live in gas-rich, star-forming galaxies. However, simulations do not predict instantaneous, global impact on molecular gas or star formation. The impact is expected to be cumulative, over multiple episodes.

Constraining the magnetic field geometry of the millisecond pulsar PSR J0030+0451 from joint radio, thermal X-ray, andγ-ray emission

Context.With the advent of multi-wavelength electromagnetic observations of neutron stars – spanning many decades in photon energies – from radio wavelengths up to X-rays andγ-rays, it has become possible to significantly constrain the geometry and the location of the associated emission regions.Aims.In this work, we use results from the modelling of thermal X-ray observations of PSR J0030+0451 from the Neutron Star Interior Composition Explorer (NICER) mission and phase-aligned radio andγ-ray pulse profiles to constrain the geometry of an off-centred dipole that is able to reproduce the light curves in these respective bands simultaneously.Methods.To this aim, we deduced a configuration with a simple dipole off-centred from the location of the centre of the thermal X-ray hot spots. We show that the geometry is compatible with independent constraints from radio andγ-ray pulsations only, leading to a fixed magnetic obliquity ofα ≈ 75° and a line-of-sight inclination angle ofζ ≈ 54°.Results.We demonstrate that an off-centred dipole cannot be rejected by accounting for the thermal X-ray pulse profiles. Moreover, the crescent shape of one spot is interpreted as the consequence of a small-scale surface dipole on top of the large-scale off-centred dipole.

The Close AGN Reference Survey (CARS): An Interplay between Radio Jets and AGN Radiation in the Radio-quiet AGN HE0040-1105

We present a case study of HE 0040-1105, an unobscured radio-quiet active galactic nucleus (AGN) at a high accretion rate of λ Edd = 0.19 ± 0.04. This particular AGN hosts an ionized gas outflow with the largest spatial offset from its nucleus compared to all other AGNs in the Close AGN Reference Survey. By combining multiwavelength observations from the Very Large Telescope/MUSE, Hubble Space Telescope/Wide Field Camera 3, Very Large Array, and European VLBI Network, we probe the ionization conditions, gas kinematics, and radio emission from host galaxy scales to the central few parsecs. We detect four kinematically distinct components, one of which is a spatially unresolved AGN-driven outflow located within the central 500 pc, where it locally dominates the interstellar medium conditions. Its velocity is too low to escape the host galaxy’s gravitational potential, and may be re-accreted onto the central black hole via chaotic cold accretion. We detect compact radio emission in HE 0040-1105 within the region covered by the outflow, varying on a timescale of ∼20 yr. We show that neither AGN coronal emission nor star formation processes wholly explain the radio morphology/spectrum. The spatial alignment between the outflowing ionized gas and the radio continuum emission on 100 pc scales is consistent with a weak jet morphology rather than diffuse radio emission produced by AGN winds. >90% of the outflowing ionized gas emission originates from the central 100 pc, within which the ionizing luminosity of the outflow is comparable to the mechanical power of the radio jet. Although radio jets might primarily drive the outflow in HE 0040-1105, radiation pressure from the AGN may contribute to this process.

Multiwavelength Study of Dark Globule DC 314.8–5.1: Point-source Identification and Diffuse Emission Characterization

We present an analysis of multiwavelength observations of the dark globule DC 314.8–5.1, using data from the Gaia optical, Two Micron All Star Survey near-infrared, and Wide-field Infrared Survey Explorer mid-infrared surveys, dedicated imaging with the Spitzer Space Telescope, and X-ray data obtained with the Swift X-Ray Telescope (XRT). The main goal was to identify possible pre-main-sequence stars (PMSs) and young stellar objects (YSOs) associated with the globule. For this, we studied the infrared colors of all point sources within the boundaries of the cloud. After removing sources with nonstellar spectra, we investigated the Gaia parallaxes for the YSO candidates and found that none are physically related to DC 314.8–5.1. In addition, we searched for X-ray emission from PMSs with Swift-XRT, and found no 0.5–10 keV emission down to a luminosity level ≲1031 erg s−1, typical of a PMS with mass ≥2 M ⊙. Our detailed inspection therefore supports a very young, “prestellar core” evolutionary stage for the cloud. Based on archival Planck and IRAS data, we moreover identify the presence of hot dust, with temperatures ≳100 K, in addition to the dominant dust component at 14 K, originating with the associated reflection nebula.

Coronal diagnostics of solar type III radio bursts using LOFAR and PSP observations

Context. Solar type III radio bursts are common phenomena, recognized as the result of accelerated electron beams propagating through the solar corona. These bursts are of particular interest as they provide valuable information about the magnetic field and plasma conditions in the corona, which are difficult to measure directly. Aims. This study aims to investigate the ambiguous source and the underlying physical processes of the type III radio bursts that occurred on April 3, 2019, through the utilization of multi-wavelength observations from the Low-Frequency Array (LOFAR) radio telescope and the Parker Solar Probe (PSP) space mission, as well as incorporating results from a Potential Field Source Surface (PFSS) and magnetohydrodynamic (MHD) models. The primary goal is to identify the spatial and temporal characteristics of the radio sources, as well as the plasma conditions along their trajectories. Methods. We applied data preprocessing techniques to combine high- and low-frequency observations from LOFAR and PSP between 2.6 kHz and 80 MHz. We then extracted information on the frequency drift and speed of the accelerated electron beams from the dynamic spectra. Additionally, we used LOFAR interferometric observations to image the sources of the radio emission at multiple frequencies and determine their locations and kinematics in the corona. Lastly, we analyzed the plasma parameters and magnetic field along the trajectories of the radio sources using PFSS and MHD model results. Results. We present several notable findings related to type III radio bursts. Firstly, through our automated implementation, we were able to effectively identify and characterize 9 type III radio bursts in the LOFAR-PSP combined dynamic spectrum and 16 type III bursts in the LOFAR dynamic spectrum. We found that the frequency drift for the detected type III bursts in the combined spectrum ranges between 0.24 and 4 MHz s−1, while the speeds of the electron beams range between 0.013 and 0.12 C. Secondly, our imaging observations show that the electrons responsible for these bursts originate from the same source and within a short time frame of fewer than 30 min. Finally, our analysis provides informative insights into the physical conditions along the path of the electron beams. For instance, we found that the plasma density obtained from the magnetohydrodynamic algorithm outside a sphere (MAS) model is significantly lower than the expected theoretical density.

Multiwavelength Observations of a B-class Flare Using XSM, AIA, and XRT

We present multiwavelength observations by Chandrayaan-2/Solar X-ray Monitor, Solar Dynamics Observatory/Atmospheric Imaging Assembly, and Hinode/X-Ray Telescope (XRT) of a B-class flare observed on 2021 February 25, originating from an active region (AR 12804) near the northwest limb. The microflare lasts for ∼30 minutes and is composed of hot loops reaching temperatures of 10 MK. We report excellent agreement (within 20%) for the average effective temperatures obtained at the flare peak from all the three instruments, which have different temperature sensitivities. The XRT filter combination of Be-thin and Be-med provides an excellent opportunity to measure the high temperatures in such microflare events. The elemental abundances during the evolution of the microflare are also studied and observed to drop toward photospheric values at the flare peak time, compared to coronal values during the rise and decay phase. This is consistent with previous XSM studies.

Hunting down white dwarf–main sequence binaries using multiwavelength observations

ABSTRACT Identification of white dwarfs (WD) with main-sequence (MS) companions and characterization of their properties can put important constraints on our understanding of binary stellar evolution and guide the theoretical predictions for a wide range of interesting transient events relevant for, e.g. Rubin Observatory (LSST), ZTF, and LISA. In this study, we combine ultraviolet (UV) and optical colour–magnitude diagrams (CMDs) to identify unresolved WD–MS binaries. In particular, we combine high-precision astrometric and photometric data in the optical from Gaia-DR3 and UV data from GALEX GR6/7 to identify 93 WD–MS candidates within 100 pc. Of these, 80 are newly identified. Using the Virtual Observatory SED Analyzer (VOSA), we fit the spectral energy distributions (SEDs) of all our candidates and derive stellar parameters, such as effective temperature, bolometric luminosity, and radius for both companions. We find that our identification method helps identify hotter and smaller WD companions (majority with ≥10 000 K and ≤0.02 R⊙) relative to the WDs identified by past surveys. We infer that these WDs are relatively more massive (median ∼ $0.76\, \mathrm{M}_\odot$). We find that most of the MS companions in our binaries are of the K and M spectral types.

Gemini near-infrared spectroscopy of high-redshift Fermi blazars: jetted black holes in the early universe were overly massive

ABSTRACT Jetted active galactic nuclei (AGNs) are the principal extragalactic γ-ray sources. Fermi-detected high-redshift (z > 3) blazars are jetted AGNs thought to be powered by massive, rapidly spinning supermassive black holes (SMBHs) in the early universe (<2 Gyr). They provide a laboratory to study early black hole (BH) growth and super-Eddington accretion – possibly responsible for the more rapid formation of jetted BHs. However, previous virial BH masses of z > 3 blazars were based on C iv λ1549 in the observed optical, but C iv λ1549 is known to be biased by strong outflows. We present new Gemini/GNIRS near-infrared spectroscopy for a sample of nine z > 3 Fermi γ-ray blazars with available multiwavelength observations that maximally sample the spectral energy distributions (SEDs). We estimate virial BH masses based on the better calibrated broad H β and/or Mg ii λ2800. We compare the new virial BH masses against independent mass estimates from SED modelling. Our work represents the first step in campaigning for more robust virial BH masses and Eddington ratios for high-redshift Fermi blazars. Our new results confirm that high-redshift Fermi blazars indeed host overly massive SMBHs as suggested by previous work, which may pose a theoretical challenge for models of the rapid early growth of jetted SMBHs.

Circumnuclear Multiphase Gas in the Circinus Galaxy. V. The Origin of the X-Ray Polarization in the Circinus Galaxy

The Imaging X-ray Polarimetry Explorer (IXPE) detected X-ray polarization in the nearest Seyfert 2 galaxy, the Circinus galaxy, for the first time. To reproduce the IXPE results, we computed the degree of polarization based on two types of radiative hydrodynamic simulations: a parsec-scale three-dimensional model and a sub-parsec-scale axisymmetric model with a higher spatial resolution. In a series of papers, we confirmed that these models naturally explain the multiwavelength observations of the Circinus galaxy from radio to X-rays. We used a Monte Carlo Simulation for Astrophysics and Cosmology code to compute the linear polarization of continuum emission. We found that the degree of polarization based on the parsec-scale radiation-driven fountain model was smaller than that observed with the IXPE. The degree of polarization based on the sub-parsec-scale model depends on the hydrogen number density of the disk (d), and the degree of polarization obtained from our simulation is consistent with that observed with the IXPE in the case of . We investigate where the photons are Compton scattered and imply that the origin of the X-ray polarization in the Circinus galaxy is the outflow inside 0.01 pc. In this case, the degree of polarization may change over a timescale of approximately 10 yr.

Off-axis Afterglow Closure Relations and Fermi-LAT Detected Gamma-Ray Bursts

Gamma-ray bursts (GRBs) are one of the most promising transient events for studying multiwavelength observations in extreme conditions. Observation of GeV photons from bursts would provide crucial information on GRB physics, including the off-axis emission. The Second Gamma-ray Burst Catalog (2FLGC) has been announced by the Fermi Large Area Telescope (Fermi-LAT) Collaboration. This catalog includes 29 bursts with photon energy higher than 10 GeV. While the synchrotron forward-shock model has well explained the afterglow data of GRBs, photon energies greater than 10 GeV are very difficult to interpret within this framework. To study the spectral and temporal indices of those bursts described in 2FLGC, Fraija et al. () proposed the closure relations (CRs) of the synchrotron self-Compton (SSC) emission of GRBs emitted from an on-axis jet, which decelerates in stellar-wind and constant-density medium. In this paper, we extend the CRs of the SSC afterglow from an on-axis scenario to an off-axis scenario, including the synchrotron afterglow radiation that seems to be off-axis. In order to investigate the spectral and temporal index evolution of those bursts reported in 2FLGC, we consider hydrodynamical evolution with energy injection in the adiabatic and radiative regimes for an electron distribution with a spectral index of 1 < p < 2 and 2 < p. The results show that the most likely scenario for synchrotron emission corresponds to the stellar wind whether or not there is energy injection and that the most likely scenario for SSC emission corresponds to the constant density when there is no energy injection and to the stellar wind when there is energy injection.

Internal Activities in a Solar Filament and Heating in Its Threads

Filaments are one of the most common features in the solar atmosphere and are of significance in solar, stellar, and laboratory plasma physics. Using data from the Chinese Hα Solar Explorer, the Solar Upper Transition Region Imager, and the Solar Dynamics Observatory, we report on multiwavelength imaging and spectral observations of the activation of a small filament. The filament activation produces several localized dynamic brightenings, which are probably produced by internal reconnections of the braided magnetic fields in the filament. The filament expands during the activation, and its threads reconnect with the ambient magnetic fields, which leads to the formation of hot arcades or loops overlying the filament. The thermal energy of each of these localized brightenings is estimated in the order of 1025–1027 erg, and the total energy is estimated to be ∼1.77 × 1028 erg. Our observations demonstrate that the internal magnetic reconnections in the filament can lead to localized heating in the filament threads and prompt external reconnections with ambient corona structures and thus could contribute to the energy and mass transferring into the corona.

Multiwavelength observations of the lensed quasar PKS 1830-211 during the 2019 γ-ray flare

ABSTRACT PKS 1830 -211 is a γ-ray emitting, high-redshift (z =2.507 ± 0.002), lensed flat-spectrum radio quasar. During the period 2019 mid-February to mid-April, this source underwent a series of strong γ-ray flares that were detected by both AGILE-GRID (Gamma-Ray Imaging Detector) and Fermi Large Area Telescope (Fermi-LAT), reaching a maximum γ-ray flux of $F_{\rm E\gt 100\, MeV}\approx 2.3\times 10^{-5}$ photons cm−2 s−1. Here, we report on a coordinated campaign from both on-ground [Medicina, Owens Valley Radio Observatory (OVRO), Rapid Eye Mount (REM), and Sardinia Radio Telescope (SRT)] and orbiting facilities (AGILE, Fermi, INTEGRAL, NuSTAR, Swift, and Chandra), with the aim of investigating the multiwavelength properties of PKS 1830-211 through nearly simultaneous observations presented here for the first time. We find a possible break in the radio spectra in different epochs above 15 GHz, and a clear maximum of the 15 GHz data approximately 110 d after the γ-ray main activity periods. The spectral energy distribution shows a very pronounced Compton dominance (&amp;gt; 200) which challenges the canonical one-component emission model. Therefore, we propose that the cooled electrons of the first component are re-accelerated to a second component by, for example, kink or tearing instability during the γ-ray flaring periods. We also note that PKS 1830-211 could be a promising candidate for future observations with both Compton satellites [e.g. enhanced ASTROGAM (e-ASTROGAM)] and Cherenkov arrays [Cherenkov Telescope Array Observatory (CTAO)] which will help, thanks to their improved sensitivity, in extending the data availability in energy bands currently uncovered.

Radio Emission and Electric Gaps in Pulsar Magnetospheres

The origin of pulsar radio emission is one of the old puzzles in theoretical astrophysics. In this Letter, we present a global kinetic plasma simulation that shows from first principles how and where radio emission can be produced in pulsar magnetospheres. We observe the self-consistent formation of electric gaps that periodically ignite electron-positron discharge. The gaps form above the polar cap and in the bulk return current. Discharge of the gaps excites electromagnetic modes, which share several features with the radio emission of real pulsars. We also observe the excitation of plasma waves and charge bunches by beam instabilities in the outer magnetosphere. Our numerical experiment demonstrates that global kinetic models can provide deep insight into the emission physics of pulsars and may help interpret their multiwavelength observations.