Source author record

Yajie Yuan

Yajie Yuan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Magnetar bursts due to Alfvén wave nonlinear breakout

The most common form of magnetar activity is short X-ray bursts, with durations from milliseconds to seconds, and luminosities ranging from $10^{36}$ to $10^{43}\ {\rm erg}\,{\rm s}^{-1}$. Recently, an X-ray burst from the galactic magnetar SGR 1935+2154 was detected to be coincident with two fast radio burst (FRB) like events from the same source, providing evidence that FRBs may be linked to magnetar bursts. Using fully 3D force-free electrodynamics simulations, we show that such magnetar bursts may be produced by Alfvén waves launched from localized magnetar quakes: a wave packet propagates to the outer magnetosphere, becomes nonlinear, and escapes the magnetosphere, forming an ultra-relativistic ejecta. The ejecta pushes open the magnetospheric field lines, creating current sheets behind it. Magnetic reconnection can happen at these current sheets, leading to plasma energization and X-ray emission. The angular size of the ejecta can be compact, $\lesssim 1$ sr if the quake launching region is small, $\lesssim 0.01$ sr at the stellar surface. We discuss implications for the FRBs and the coincident X-ray burst from SGR 1935+2154.

preprint2020arXiv

A Numerical Model for the Multi-wavelength Lightcurves of PSR J0030+0451

Recent modeling of Neutron Star Interior Composition Explorer(NICER) observations of the millisecond pulsar PSR J0030+0451 suggests that the magnetic field of the pulsar is non-dipolar. We construct a magnetic field configuration where foot points of the open field lines closely resemble the hotspot configuration from NICER observations. Using this magnetic field as input, we perform force-free simulations of the magnetosphere of PSR J0030+0451, showing the three-dimensional structure of its plasma-filled magnetosphere. Making simple and physically motivated assumptions about the emitting regions, we are able to construct the multi-wavelength lightcurves that qualitatively agree with the corresponding observations. The agreement suggests that multipole magnetic structures are the key to modeling this type of pulsars, and can be used to constrain the magnetic inclination angle and the location of radio emission.

preprint2020arXiv

Alfvén wave mode conversion in pulsar magnetospheres

The radio emission anomaly coincident with the 2016 glitch of the Vela pulsar may be caused by a star quake that launches Alfvén waves into the magnetosphere, disturbing the original radio emitting region. To quantify the lifetime of the Alfvén waves, we investigate a possible energy loss mechanism, the conversion of Alfvénwaves into fast magnetosonic waves. Using axisymmetric force-free simulations, we follow the propagation of Alfvén waves launched from the stellar surface with small amplitude into the closed zone of a force-free dipolar pulsar magnetosphere. We observe mode conversion happening in the ideal force-free regime. The conversion efficiency during the first passage of the Alfvén wave through the equator can be large, for waves that reach large amplitudes as they travel away from the star, or propagate on the field lines passing close to the Y-point. However, the conversion efficiency is reduced due to dephasing on subsequent passages and considerable Alfvén power on the closed field lines remains. Thus while some leakage into the fast mode happens, we need detailed understanding of the original quenching in order to say whether mode conversion alone can lead to reactivation of the pulsar on a short timescale.

preprint2020arXiv

Plasmoid ejection by Alfvén waves and the fast radio bursts from SGR 1935+2154

Using numerical simulations we show that low-amplitude Alfvén waves from a magnetar quake propagate to the outer magnetosphere and convert to "plasmoids" (closed magnetic loops) which accelerate from the star, driving blast waves into the magnetar wind. Quickly after its formation, the plasmoid becomes a thin relativistic pancake. It pushes out the magnetospheric field lines, and they gradually reconnect behind the pancake, generating a variable wind far stronger than the normal spindown wind of the magnetar. Repeating ejections drive blast waves in the amplified wind. We suggest that these ejections generate the simultaneous X-ray and radio bursts detected from SGR 1935+2154. A modest energy budget of the magnetospheric perturbation $\sim 10^{40}$ erg is sufficient to produce the observed bursts. Our simulation predicts a narrow (a few ms) X-ray spike from the magnetosphere, arriving almost simultaneously with the radio burst emitted far outside the magnetosphere. This timing is caused by the extreme relativistic motion of the ejecta.

preprint2020arXiv

Relativistic Alfvén Waves Entering Charge Starvation in the Magnetospheres of Neutron Stars

Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current $j_{\rm A}$. We derive an analytic expression for the evolution of the wave vector $\boldsymbol{k}$ and the growth of $j_{\rm A}$. In the strongly sheared regime, $j_{\rm A}$ may exceed the maximum current $j_{0}$ that can be supported by the background $e^{\pm}$ plasma. We investigate these "charge-starved" waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave continues to propagate successfully even when $κ\equiv j_{\rm A}/j_{0} \gg 1$. It sustains $j_{\rm A}$ by compressing and advecting the plasma along the magnetic field lines with particle Lorentz factors $\sim κ^{1/2}$. The simulations show how plasma instabilities lead to gradual dissipation of the wave energy, giving a dissipation power $L_{\rm diss}\sim 10^{35}(κ/100)^{1/2} (B_w/10^{11}\,{\rm G})\,\mathrm{erg/s}$, where $B_w$ is the wave amplitude. Our results imply that dissipation due to charge starvation is not sufficient to power observed fast radio bursts (FRBs), in contrast to recent proposals.

preprint2019arXiv

Physics of Pair Producing Gaps in Black Hole Magnetospheres II -- General Relativity

This is the second paper in a series where we examine the physics of pair producing gaps in low-luminosity accreting supermassive black hole systems. In this paper, we carry out time-dependent self-consistent fully general relativistic 1D PIC simulations of the gap, including full inverse Compton scattering and photon tracking. Similar to the previous paper, we find a highly time-dependent solution where a macroscopic vacuum gap can open quasi-periodically, producing bursts of $e^\pm$ pairs and high energy radiation. We present the light curve, particle and photon spectra from this process. Using an empirical scaling relation, we rescale the parameters to the inferred values at the base of the jet in M87, and find that the observed TeV flares could potentially be explained by this model under certain parameter assumptions.

preprint2016arXiv

Kinetic simulations of the lowest-order unstable mode of relativistic magnetostatic equilibria

We present the results of particle-in-cell numerical pair plasma simulations of relativistic 2D magnetostatic equilibria known as the 'ABC' fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of 'X-point collapse' first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration time scale. We describe in detail the evolving structure of the dynamical current layers, and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

preprint2016arXiv

Kinetic study of radiation-reaction-limited particle acceleration during the relaxation of unstable force-free equilibria

Many powerful and variable gamma-ray sources, including pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, seem capable of accelerating particles to gamma-ray emitting energies efficiently over very short time scales. These are likely due to rapid dissipation of electromagnetic energy in a highly magnetized, relativistic plasma. In order to understand the generic features of such processes, we have investigated simple models based on relaxation of unstable force-free magnetostatic equilibria. In this work, we make the connection between the corresponding plasma dynamics and the expected radiation signal, using 2D particle-in-cell simulations that self-consistently include synchrotron radiation reaction. We focus on the lowest order unstable force-free equilibrium in a 2D periodic box. We find that rapid variability, with modest apparent radiation efficiency as perceived by a fixed observer, can be produced during the evolution of the instability. The "flares" are accompanied by an increased polarization degree in the high energy band, with rapid variation in the polarization angle. Furthermore, the separation between the acceleration sites and the synchrotron radiation sites for the highest energy particles facilitates acceleration beyond the synchrotron radiation reaction limit. We also discuss the dynamical consequences of radiation reaction, and some astrophysical applications of this model. Our current simulations with numerically tractable parameters are not yet able to reproduce the most dramatic gamma-ray flares, e.g., from Crab Nebula. Higher magnetization studies are promising and will be carried out in the future.

preprint2015arXiv

On the Implications of Recent Observations of the Inner Knot in the Crab Nebula

Recent observations of the Crab Nebula (Rudy et al 2015) have maintained its reputation for high energy astrophysical enlightenment and its use as a testbed for theories of the behaviour of magnetized, relativistic plasma. In particular, new observations of the inner knot located 0.65" SE from the pulsar confirm that it is compact, elongated transversely to the symmetry axis and curved concave towards the pulsar. 60 percent polarization has been measured along the symmetry axis (Moran et al 2013). The knot does not appear to be involved in the gamma ray flares. The new observations both reinforce the interpretation of the knot as dissipation of the pulsar wind at a strong shock and challenge the details of existing models of this process. In particular, it is argued that the compactness, high polarization and curvature are difficult to reconcile with simple relativistic shock models. Alternative possibilities include deflection of the outflow ahead of the shock and spatial variation in which the knot is interpreted as a caustic. Some future observations are proposed and new theoretical investigations are suggested.

preprint2015arXiv

Spontaneous decay of periodic magnetostatic equilibria

In order to understand the conditions which lead a highly magnetized, relativistic plasma to become unstable, and in such cases how the plasma evolves, we study a prototypical class of magnetostatic equilibria where the magnetic field satisfies $\nabla \times\mathbf B = α\mathbf B$, where $α$ is spatially uniform, on a periodic domain. Using numerical solutions we show that generic examples of such equilibria are unstable to ideal modes (including incompressible ones) which are marked by exponential growth in the linear phase. We characterize the unstable mode, showing how it can be understood in terms of merging magnetic and current structures, and explicitly demonstrate its instability using the energy principle. Following the nonlinear evolution of these solutions, we find that they rapidly develop regions with relativistic velocities and electric fields of comparable magnitude to the magnetic field, liberating magnetic energy on dynamical timescales and eventually settling into a configuration with the largest allowable wavelength. These properties make such solutions a promising setting for exploring the mechanisms behind extreme cosmic sources of gamma rays.

preprint2014arXiv

Cosmic Ray Origins: An Introduction

Physicists have pondered the origin of cosmic rays for over a hundred years. However the last few years have seen an upsurge in the observation, progress in the theory and a genuine increase in the importance attached to the topic due to its intimate connection to the indirect detection of evidence for dark matter. The intent of this talk is to set the stage for the meeting by reviewing some of the basic features of the entire cosmic ray spectrum from GeV to ZeV energy and some of the models that have been developed. The connection will also be made to recent developments in understanding general astrophysical particle acceleration in pulsar wind nebulae, relativistic jets and gamma ray bursts. The prospects for future discoveries, which may elucidate the origin of cosmic rays, are bright.

preprint2013arXiv

Chandra, Keck and VLA Observations of the Crab Nebula during the 2011-April Gamma-ray Flare

We present results from our analysis of Chandra X-ray Observatory, W. M. Keck Observatory, and Karl G. Jansky Very Large Array (VLA) images of the Crab Nebula that were contemporaneous with the gamma-ray flare of 2011 April. Despite hints in the X-ray data, we find no evidence for statistically significant variations that pinpoint the specific location of the flares within the Nebula. The Keck observations extend this conclusion to the "inner knot", i.e., the feature within an arcsecond of the pulsar. The VLA observations support this conclusion. We also discuss theoretical implications of the gamma-ray flares and suggest that the most dramatic gamma-ray flares are due to radiation-reaction-limited synchrotron emission associated with sudden, dissipative changes in the current system sustained by the central pulsar.