Source author record

Siming Liu

Siming Liu 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

42works
13topics
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

42 published item(s)

preprint2026arXiv

EduBench: A Comprehensive Benchmarking Dataset for Evaluating Large Language Models in Diverse Educational Scenarios

As large language models continue to advance, their application in educational contexts remains underexplored and under-optimized. In this paper, we address this gap by introducing the first diverse benchmark tailored for educational scenarios, incorporating synthetic data containing 9 major scenarios and over 4,000 distinct educational contexts. To enable comprehensive assessment, we propose a set of multi-dimensional evaluation metrics that cover 12 critical aspects relevant to both teachers and students. We further apply human annotation to ensure the effectiveness of the model-generated evaluation responses. Additionally, we succeed to train a relatively small-scale model on our constructed dataset and demonstrate that it can achieve performance comparable to state-of-the-art large models (e.g., Deepseek V3, Qwen Max) on the test set. Overall, this work provides a practical foundation for the development and evaluation of education-oriented language models. Code and data are released at https://github.com/ybai-nlp/EduBench.

preprint2022arXiv

A Three-component Model for Cosmic-ray Spectrum and Dipole Anisotropy

Using a three-component, multi-scale diffusion model, we show that the cosmic-ray (CR) proton and helium spectra and the dipole anisotropy can be explained with reasonable parameters. The model includes a nearby source associated with the supernova remnant (SNR) that gave rise to the Geminga pulsar, a source at the Galactic center, and a component associated with the Galactic disk. The CR flux below TeV is dominated by the disk component. The center source with a continuous injection of CRs starting about 18 Myr ago is needed to explain the anisotropy above 100 TeV. With the assumption of universal CR spectra injected by all SNRs, the nearby source can produce a TeV spectral bump observed at Earth via slow diffusion across the interstellar magnetic field, which needs to have an angle $ θ\approx 5^{\circ} $ between the field line and the line of sight toward the source, and have weak magnetic turbulence with the Alfvén Mach number $ M_{\text{A}}\approx 0.1 $. Considering the modulation of the Galactic-scale anisotropy by this magnetic field, in a quasi-local approach the field may be directed at a right ascension about $ -90^{\circ} $ and a declination about $ -7.4^{\circ} $ in the equatorial coordinate system.

preprint2022arXiv

Observational Signatures of Tearing Instability in the Current Sheet of a Solar Flare

Magnetic reconnection is a fundamental physical process converting magnetic energy into not only plasma energy but also particle energy in various astrophysical phenomena. In this letter, we show a unique dataset of a solar flare where various plasmoids were formed by a continually stretched current sheet. EUV images captured reconnection inflows, outflows, and particularly the recurring plasma blobs (plasmoids). X-ray images reveal nonthermal emission sources at the lower end of the current sheet, presumably as large plasmoids with a sufficiently amount of energetic electrons trapped in. In the radio domain, an upward slowly drifting pulsation structure, followed by a rare pair of oppositely drifting structures, was observed. These structures are supposed to map the evolution of the primary and the secondary plasmoids formed in the current sheet. Our results on plasmoids at different locations and scales shed important light on the dynamics, plasma heating, particle acceleration, and transport processes in the turbulent current sheet and provide observational evidence for the cascading magnetic reconnection process.

preprint2022arXiv

The Large High Altitude Air Shower Observatory (LHAASO) Science Book (2021 Edition)

Since the science white paper of the Large High Altitude Air Shower Observatory (LHAASO) published on arXiv in 2019 [e-Print: 1905.02773 (astro-ph.HE)], LHAASO has completed the transition from a project to an operational gamma-ray astronomical observatory LHAASO is a new generation multi-component facility located in Daocheng, Sichuan province of China, at an altitude of 4410 meters. It aims at measuring with unprecedented sensitivity the spectrum, composition, and anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$~eV, and acting simultaneously as a wide aperture (one stereoradiant) continuously operating gamma-ray telescope in the energy range between 10$^{11}$ and $10^{15}$~eV with the designed sensitivity of 1.3\% of the Crab Unit (CU) above 100 TeV. LHAASO's capability of measuring simultaneously different shower components (electrons, muons, and Cherenkov/fluorescence light), will allow it to investigate the origin, acceleration, and propagation of CR through measurement of the energy spectrum, elemental composition, and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO will play a key role in CR physics and gamma-ray astronomy for a general and comprehensive exploration of the high energy universe and will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. The LHAASO Collaboration organized an editorial working group and finished all editorial work of this science book, to summarize the instrumental features and outline the prospects of scientific researches with the LHAASO experiment.

preprint2022arXiv

Wave emission of non-thermal electron beams generated by magnetic reconnection

Magnetic reconnection in solar flares can efficiently generate non-thermal electron beams. The energetic electrons can, in turn, cause radio waves through microscopic plasma instabilities as they propagate through the ambient plasma along the magnetic field lines. We aim at investigating the wave emission caused by fast moving electron beams (FEBs) with characteristic non-thermal electron velocity distribution functions (EVDFs) generated by kinetic magnetic reconnection: two-streaming EVDFs along the separatrices and in the diffusion region, and perpendicular crescent-shaped EVDFs closer to the diffusion region. For this purpose, we utilized 2.5D fully kinetic Particle-In-Cell (PIC) code simulations in this study. We found that:(1) the two-streaming EVDFs plus the background ions are unstable to electron/ion (streaming) instabilities which cause ion acoustic waves and Langmuir waves due to the net current. This can lead to multiple harmonic plasma emission in the diffusion region and the separatrices of reconnection. (2) The perpendicular crescent-shaped EVDFs can cause multiple harmonic electromagnetic electron cyclotron waves through the electron cyclotron maser instabilities in the diffusion region of reconnection. Our results are applicable to diagnose the plasma parameters which control magnetic reconnection in solar flares by means of radio waves observations.

preprint2021arXiv

Radial Profiles of Non-thermal Emission from Supernova Remnant RX J1713.7-3946

Supernova remnant RX J1713.7-3946 (also named as G347.3-0.5) has exhibited largest surface brightness, detailed spectral and shell-type morphology, it is one of the brightest TeV sources. The recent H.E.S.S. observation of RX J1713.7-3946 revealed \textbf{a} broken power-law spectrum of GeV-TeV gamma-ray spectrum and more extended gamma-ray spatial radial profile than that in the X-ray band. Based on the diffusion shock acceleration model, we solve spherically symmetric hydrodynamic equations and transport equations of particles, and investigate multi-band non-thermal emission of RX J1713.7-3946 and radial profiles of its surface brightness for two selected zones in the leptonic scenario for the $γ$-ray emission. We found (1) the diffusion coefficient has a weak energy-dependent, and the Kolmogorov type is favored; (2) the magnetic field strength could vary linearly or nonlinearly with radius for different surrounding environments because of possible turbulence in shock downstream region, and a compressional amplification is likely to exist at the shock front; (3) the non-thermal photons from radio to X-ray bands are dominated by synchrotron emission from relativistic electrons, if the GeV-TeV gamma-rays are produced by inverse Compton scattering from these electrons interacting with the background photons, then the X-ray and gamma-ray radial profiles can be reproduced except for the more extended $γ$-ray emission.

preprint2020arXiv

Crowd Scene Analysis by Output Encoding

Crowd scene analysis receives growing attention due to its wide applications. Grasping the accurate crowd location (rather than merely crowd count) is important for spatially identifying high-risk regions in congested scenes. In this paper, we propose a Compressed Sensing based Output Encoding (CSOE) scheme, which casts detecting pixel coordinates of small objects into a task of signal regression in encoding signal space. CSOE helps to boost localization performance in circumstances where targets are highly crowded without huge scale variation. In addition, proper receptive field sizes are crucial for crowd analysis due to human size variations. We create Multiple Dilated Convolution Branches (MDCB) that offers a set of different receptive field sizes, to improve localization accuracy when objects sizes change drastically in an image. Also, we develop an Adaptive Receptive Field Weighting (ARFW) module, which further deals with scale variation issue by adaptively emphasizing informative channels that have proper receptive field size. Experiments demonstrate the effectiveness of the proposed method, which achieves state-of-the-art performance across four mainstream datasets, especially achieves excellent results in highly crowded scenes. More importantly, experiments support our insights that it is crucial to tackle target size variation issue in crowd analysis task, and casting crowd localization as regression in encoding signal space is quite effective for crowd analysis.

preprint2020arXiv

From period to quasi-period to chaos: A continuous spectrum of orbits of charged particles trapped in a dipole magnetic field

Via evaluation of the Lyapunov exponent, we report the discovery of three prominent sets of phase space regimes of quasi-periodic orbits of charged particles trapped in a dipole magnetic field. Besides the low energy regime that has been studied extensively and covers more than 10% in each dimension of the phase space of trapped orbits, there are two sets of high energy regimes, the largest of which covers more than 4% in each dimension of the phase space of trapped orbits. Particles in these high energy orbits may be observed in space and be realized in plasma experiments on the Earth. It is well-known that there are quasi-periodic orbits around stable periodic orbits in Hamiltonian systems with 2 degrees of freedom and these quasi-periodic orbits are stable as well. Since periodic orbits appear to have a negligible measure in the phase space, they are difficult to realize in nature. Quasi-periodic orbits, on the other hand, may occupy a finite volume in the 4 dimensional (4D) phase space and be readily detectable. A chaotic orbit has at least one positive Lyapunov exponent. The Lyapunov exponents of quasi-periodic orbits, on the other hand, should be zero. Via calculation of the Lyapunov exponent of orbits of trapped charged particles in a dipole magnetic field, we scanned the corresponding phase space and found several prominent regimes of quasi-periodic orbits associated with stable periodic orbits in the equatorial plane. These regimes appear to be connected to some small regimes of quasi-periodic orbits associated with stable periodic orbits in the Meridian plane. Our numerical results also show a continuous spectrum of these orbits from stable periodic, to quasi-periodic with vanishing Lyapunov exponents, and eventually to chaotic ones with at least one positive Lyapunov exponent and there are unstable periodic orbits with a positive maximum Lyapunov exponent.

preprint2020arXiv

Hadronic vs leptonic models for $γ$-ray emission from VER J2227+608

Recent observations of VER J2227+608 reveal a broken power $γ$-ray spectrum with the spectral index increasing from $\sim 1.8$ in the GeV energy range to $\sim 2.3$ in the TeV range. Such a spectral break can be attributed to radiative energy loss of energetic electrons in the leptonic scenario for the $γ$-ray emission, which, in combination with characteristic age of the nearby pulsar, can be used to constrain magnetic field in the emission region. We show that the radio and X-ray observations can also be explained in such a scenario. In the hadronic scenario, the spectral break can be attributed to diffusion of energetic ions in a turbulent medium and detailed spectral measurement can be used to constrain the diffusion coefficient. These two models, however, predict drastically different spectra above 100 TeV, which will be uncovered with future high-resolution observations, such as LHAASO and/or CTA.

preprint2020arXiv

Is PSR J0855$-$4644 responsible for the 1.4 TeV electron spectral bump hinted by DAMPE?

DAMPE observation on the cosmic ray electron spectrum hints a narrow excess at $\sim$ 1.4 TeV. Although the excess can be ascribed to dark matter particles, pulsars and pulsar wind nebulae are believed to be a more natural astrophysical origin: electrons injected from nearby pulsars at their early ages can form a bump-like feature in the spectrum due to radiative energy losses. In this paper, with a survey of nearby pulsars, we find 4 pulsars that may have notable contributions to $\sim$ 1.4 TeV cosmic ray electrons. Among them, PSR J0855$-$4644 has a spin down luminosity more than 50 times higher than others and presumably dominates the electron fluxes from them. X-ray observations on the inner compact part (which may represent a tunnel for the transport of electrons from the pulsar) of PWN G267.0$-$01.0 are then used to constrain the spectral index of high energy electrons injected by the pulsar. We show that high-energy electrons released by PSR J0855$-$4644 could indeed reproduce the 1.4 TeV spectral feature hinted by the DAMPE with reasonable parameters.

preprint2020arXiv

Origin of Strong Linear Polarization from Fast Radio Bursts

The detection of almost 100% linearly polarized emission from the fast radio burst source FRB 121102 implies coherent emission of relativistic electrons moving in perpendicular to the ambient magnetic field. The origin of such a particle distribution is very intriguing. Given that FRB 121102 is likely driven by a neutron star, we explored orbits of charged particles trapped in a dipole magnetic field (the Störmer problem). Most previous studies focused on particles with relatively low energies so that the guiding center approximation may be applied. High energy particles usually have chaotic orbits except those on a periodic orbit or near stable periodic orbits. Via evaluation of the maximum Lyapunov exponent of orbits of particles launched from the equatorial plane with an axial velocity (the angular velocity sets the length and energy scales of the system), we found prominent regions of quasi-periodic orbits around stable periodic orbits in the equatorial plane at high-energies. Particles in these orbits oscillate around the equatorial plane and their radial distance from the dipole can vary by a factor of/,2. Relativistic electrons in such orbits may be responsible for the almost 100% polarized emission from FRB 121102.

preprint2020arXiv

Self-similar Blast Wave for A Two-component Fluid with Variable Adiabatic Index

We propose a self-similar (SS) solution to hydrodynamic non-relativistic flow behind a spherical strong blast wave (BW) passing through a homogeneous plasma with efficient relativistic particle acceleration at the shock front. The flow is described by an ideal two-fluid model with a relativistic component so that the post-shock gas has an effective SS adiabatic index $ γ$ varying from $ 5/3 $ to $ 4/3 $. This solution is calculated numerically and compared with the standard Sedov solution. We find that the BW center in our solution is dominated by the relativistic component with $ γ=4/3 $ for the divergence of expansion there, and the relativistic component dominates the interior for a moderate acceleration efficiency at the shock front. The overall efficiency of relativistic particle acceleration can be enhanced by a factor of $ 2 $ due to the slower adiabatic energy loss rate of the relativistic component during expansion. Tendency of the dominance by the relativistic component may be common in expanding astrophysical two-fluid systems such as supernova remnants, lobes of radio galaxies.

preprint2019arXiv

Discovery of a spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099: a TeV halo or supernova remnant ?

Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we search for possible GeV counterparts of three TeV halo candidates with Fermi Large Area Telescopes. We detect a new spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099, which is also detected by HESS (HESS J1912+101). We find that the size of the GeV source is significantly larger than that of the TeV emission measured by HESS, and a spatial template characteristic of a PWN can fit the GeV data. We suggest that the GeV source is an extended PWN powered by the central middle-aged pulsar PSR J1913+1011. This discovery favors the TeV halo scenario for the TeV source 2HWC J1912+099 (HESS J1912+101), although the possible shell-like morphology measured by HESS challenges this interpretation. Alternatively, the TeV emission could be dominated by a supernova remnant (SNR) via the hadronic process. Future multi-wavelength observations of the source and more precise measurements of the spatial profile of the TeV emission will be useful to distinguish between the two scenarios.

preprint2019arXiv

VER J2227+608: A Hadronic PeVatron Pulsar Wind Nebula?

We report the detection of GeV $γ$-ray emission from the very-high-energy (VHE) $γ$-ray source VER J2227+608 associated with the "tail" region of SNR G106.3+2.7. The GeV $γ$-ray emission is extended and spatially coincident with molecular clouds traced by CO emission. The broadband GeV to TeV emission of VER J2227+608 can be well fitted by a single power-law function with an index of 1.90$\pm$0.04, without obvious indication of spectral cutoff toward high energies. The pure leptonic model for the $γ$-ray emission can be marginally ruled out by the X-ray and TeV data. In the hadronic model, the low energy content of CRs and the hard $γ$-ray spectrum, in combination with the center-bright source structure, suggest that VER J2227+608 may be powered by the PWN instead of shocks of the SNR. And the cutoff energy of the proton distribution needs to be higher than $\sim$ 400 TeV, which makes it an attractive PeVatron candidate. Future observations by the upcoming Large High Altitude Air Shower Observatory (LHAASO) and the Cherenkov Telescope Array in the north (CTA-North) could distinguish these models and constrain the maximum energy of cosmic rays in supernova remnants.

preprint2016arXiv

Energy Flows in the Jupiter-Io System

With the laws of mass conservation, momentum conservation and energy conservation, incorporating the processes of neutral gas ionization and ion diffusion, we develop a self-consistent model for the bright ribbon --- the most prominent feature in Io's plasma torus. The model parameters are well constrained by earlier {\it in situ} observations with the Galileo and Voyager spacescrafts. Our model calculation indicates that the total power dissipated inside the torus is 3.6 times bigger than the total power transported to Jovian ionosphere via Birkeland current. The power dissipation inside the torus is relatively uniform. Most of the power transportation associated with the Birkeland current, however, is localized near the flux tube of Io. With a height-intergrated conductivity of 0.15$\,$mho in Jovian ionosphere, consistent with earlier aeronomy models, the model gives a reasonable fit to the recent observations of the FUV Io tail on Jupiter. Extra mass loading near Io is required in the model. This excess of mass injection into the plasma torus may explain the Io-correlated energy source revealed by earlier EUV observations with the Voyager spacecrafts. Further investigation of the model implications is warranted.

preprint2016arXiv

On the Power-Law Distributions of X-ray Fluxes from Solar Flares Observed with GOES

Power-law frequency distributions of the peak flux of solar flare X-ray emission have been studied extensively and attributed to a system of self-organized criticality (SOC). In this paper, we first show that, so long as the shape of the normalized light curve is not correlated with the peak flux, the flux histogram of solar flares also follows a power-law distribution with the same spectral index as the power-law frequency distribution of the peak flux, which may partially explain why power-law distributions are ubiquitous in the Universe. We then show that the spectral indexes of the histograms of soft X-ray fluxes observed by GOES satellites in two different energy channels are different: the higher energy channel has a harder distribution than the lower energy channel, which challenges the universal power-law distribution predicted by SOC models and implies a very soft distribution of thermal energy content of plasmas probed by the GOES. The temperature ($T$) distribution, on the other hand, approaches a power-law distribution with an index of 2 for high values of $T$. Application of SOC models to statistical properties of solar flares needs to be revisited.

preprint2015arXiv

Limitation of the Least Square Method in the Evaluation of Dimension of Fractal Brownian Motions

With the standard deviation for the logarithm of the re-scaled range $\langle |F(t+τ)-F(t)|\rangle$ of simulated fractal Brownian motions $F(t)$ given in a previous paper \cite{q14}, the method of least squares is adopted to determine the slope, $S$, and intercept, $I$, of the log$(\langle |F(t+τ)-F(t)|\rangle)$ vs $\rm{log}(τ)$ plot to investigate the limitation of this procedure. It is found that the reduced $χ^2$ of the fitting decreases with the increase of the Hurst index, $H$ (the expectation value of $S$), which may be attributed to the correlation among the re-scaled ranges. Similarly, it is found that the errors of the fitting parameters $S$ and $I$ are usually smaller than their corresponding standard deviations. These results show the limitation of using the simple least square method to determine the dimension of a fractal time series. Nevertheless, they may be used to reinterpret the fitting results of the least square method to determine the dimension of fractal Brownian motions more self-consistently. The currency exchange rate between Euro and Dollar is used as an example to demonstrate this procedure and a fractal dimension of 1.511 is obtained for spans greater than 30 transactions.

preprint2014arXiv

Fermi Large Area Telescope Detection of Supernova Remnant RCW 86

Using $5.4$ year Fermi-LAT data, we report the detection of GeV $γ$-ray emission from the shell-type supernova remnant RCW 86 (G315.4-2.3) with a significance of $\sim5.1σ$. The data slightly favors an extended emission of this supernova remnant. The spectral index of RCW 86 is found to be very hard, $Γ\sim1.4$, in the $0.4$ to $300$ GeV range. A one zone leptonic model can well fit the multi-wavelength data from radio to very high energy $γ$-rays. The very hard GeV $γ$-ray spectrum and the inferred low gas density seem to disfavor the hadronic origin of the $γ$-rays. The $γ$-ray behavior of RCW 86 is very similar to several other TeV shell-type supernova remnants, e.g., RX J1713.7-3946, RX J0852.0-4622, SN 1006 and HESS J1731-347.

preprint2014arXiv

Fermi Large Area Telescope observations of the supernova remnant HESS J1731-347

Context: HESS J1731-347 has been identified as one of the few TeV-bright shell-type supernova remnants (SNRs). These remnants are dominated by nonthermal emission, and the nature of TeV emission has been continuously debated for nearly a decade. Aims: We carry out the detailed modeling of the radio to gamma-ray spectrum of HESS J1731-347 to constrain the magnetic field and energetic particles sources, which we compare with those of the other TeV-bright shell-type SNRs explored before. Methods: Four years of data from Fermi Large Area Telescope (LAT) observations for regions around this remnant are analyzed, leading to no detection correlated with the source discovered in the TeV band. The Markov Chain Monte Carlo method is used to constrain parameters of one-zone models for the overall emission spectrum. Results: Based on the 99.9% upper limits of fluxes in the GeV range, one-zone hadronic models with an energetic proton spectral slope greater than 1.8 can be ruled out, which favors a leptonic origin for the gamma-ray emission, making this remnant a sibling of the brightest TeV SNR RX J1713.7-3946, the Vela Junior SNR RX J0852.0-4622, and RCW 86. The best-fit leptonic model has an electron spectral slope of 1.8 and a magnetic field of about 30 muG, which is at least a factor of 2 higher than those of RX J1713.7-3946 and RX J0852.0-4622, posing a challenge to the distance estimate and/or the energy equipartition between energetic electrons and the magnetic field of this source. A measurement of the shock speed will address this challenge and has implications on the magnetic field evolution and electron acceleration driven by shocks of SNRs.

preprint2014arXiv

The structure of TeV-bright shell-type supernova remnants

Aims. Two-dimensional MHD simulations are used to model the emission properties of TeV-bright shell-type supernova remnants (SNRs) and to explore their nature. Methods. In the leptonic scenario for the TeV emission, the $γ$-ray emission is produced via Inverse Compton scattering of background soft photons by high-energy electrons accelerated by the shocks of the SNRs. The TeV emissivity is proportional to the magnetic field energy density and MHD simulations can be used to model the TeV structure of such remnants directly. 2D MHD simulations for SNRs are then performed under the assumption that the ambient interstellar medium is turbulent with the magnetic field and density fluctuations following a Kolmogorov-like power-law spectrum. Results. (1) As expected, these simulations confirm early 1D and 2D modelings of these sources, namely the hydrodynamical evolution of the shock waves and amplification of magnetic field by Rayleigh-Taylor convective flows and by shocks propagating in a turbulent medium; (2) We reproduce rather complex morphological structure for $γ$-rays, suggesting intrinsic variations of the source morphology not related to the structure of the progenitor and environment; (3)The observed radial profile of several remnants are well reproduced with an ambient medium density of $0.1-1$ cm$^{-3}$. An even lower ambient density leads to a sharper drop of the TeV brightness with radius than what is observed near the outer edge of these remnants. Conclusions. In a turbulent background medium, we can reproduce the observed characteristics of several shell-type TeV SNRs with reasonable parameters except for a higher ambient density than that inferred from X-ray observations.

preprint2013arXiv

Energy Partition between Energetic Electrons and Turbulent Magnetic Field in Supernova Remnant RX J1713.7-3946

Current observations of supernova remnant (SNR) RX J1713.7-3946 favor the leptonic scenario for the TeV emission, where the radio to X-ray emission is produced via the synchrotron process and the $γ$-ray emission is produced via the inverse Comptonization of soft background photons, and the electron distribution can be inferred from the observed $γ$-ray spectrum with a spectral inversion method. It is shown that the observed correlation between the X-ray and $γ$-ray brightness of SNR RX J1713.7-3946 can be readily explained with the assumption that the energy density of energetic electrons is proportional to that of the magnetic field in such a scenario. A 2D magnetohydrodynamic simulation is then carried out to model the overall emission spectrum. It is found that the total energy of electrons above $\sim 1$ GeV is equal to that of the magnetic field. This is the first piece of observational evidence for energy equipartition between energetic electrons and magnetic field in the downstream of strong collision-less astrophysical shocks of SNRs.

preprint2013arXiv

Error Assessment in Modeling with Fractal Brownian Motions

To model a given time series $F(t)$ with fractal Brownian motions (fBms), it is necessary to have appropriate error assessment for related quantities. Usually the fractal dimension $D$ is derived from the Hurst exponent $H$ via the relation $D=2-H$, and the Hurst exponent can be evaluated by analyzing the dependence of the rescaled range $\langle|F(t+τ)-F(t)|\rangle$ on the time span $τ$. For fBms, the error of the rescaled range not only depends on data sampling but also varies with $H$ due to the presence of long term memory. This error for a given time series then can not be assessed without knowing the fractal dimension. We carry out extensive numerical simulations to explore the error of rescaled range of fBms and find that for $0<H<0.5$, $|F(t+τ)-F(t)|$ can be treated as independent for time spans without overlap; for $0.5<H<1$, the long term memory makes $|F(t+τ)-F(t)|$ correlated and an approximate method is given to evaluate the error of $\langle|F(t+τ)-F(t)|\rangle$. The error and fractal dimension can then be determined self-consistently in the modeling of a time series with fBms.

preprint2013arXiv

Impulsive Thermal X-ray Emission from a Low-lying Coronal Loop

Understanding the relationship among different emission components plays an essential role in the study of particle acceleration and energy conversion in solar flares. In flares where gradual and impulsive emission components can be readily identified the impulsive emission has been attributed to non-thermal particles. We carry out detailed analysis of Hα and X-ray observations of a GOES class B microflare loop on the solar disk. The impulsive hard X-ray emission, however, is found to be consistent with a hot, quasi-thermal origin, and there is little evidence of emission from chromospheric footpoints, which challenges conventional models of flares and reveals a class of microflares associated with dense loops. Hα observations indicate that the loop lies very low in the solar corona or even in the chromosphere and both emission and absorption materials evolve during the flare. The enhanced Hα emission may very well originate from the photosphere when the low-lying flare loop heats up the underlying chromosphere and reduces the corresponding Hα opacity. These observations may be compared with detailed modeling of flare loops with the internal kink instability, where the mode remains confined in space without apparent change in the global field shape, to uncover the underlying physical processes and to probe the structure of solar atmosphere.

preprint2012arXiv

Stochastic Electron Acceleration in SNR RX J1713.7-3946

Stochastic acceleration of charged particles due to their interactions with plasma waves may be responsible for producing superthermal particles in a variety of astrophysical systems. This process can be described as a diffusion process in the energy space with the Fokker-Planck equation. In this paper, a time-dependent numerical code is used to solve the reduced Fokker-Planck equation involving only time and energy variables with general forms of the diffusion coefficients. We also propose a self-similar model for particle acceleration in Sedov explosions and use the TeV SNR RX J1713.7-3946 as an example to demonstrate the model characteristics. Markov Chain Monte Carlo method is utilized to constrain model parameters with observations.

preprint2012arXiv

Unified model for the gamma-ray emission of supernova remnants

Shocks of supernova remnants (SNRs) are important (and perhaps the dominant) agents for production of the Galactic cosmic rays. Recent $γ$-ray observations of several SNRs have made this case more compelling. However, these broadband high-energy measurements also reveal a variety of spectral shape demanding more comprehensive modeling of emissions from SNRs. According to the locally observed fluxes of cosmic ray protons and electrons, the electron-to-proton number ratio is known to be about 1%. Assuming such a ratio is universal for all SNRs and identical spectral shape for all kinds of accelerated particles, we propose a unified model that ascribes the distinct $γ$-ray spectra of different SNRs to variations of the medium density and the spectral difference between cosmic ray electrons and protons observed at Earth to transport effects. For low density environments, the $γ$-ray emission is inverse-Compton dominated. For high density environments like systems of high-energy particles interacting with molecular clouds, the $γ$-ray emission is $π^0$-decay dominated. The model predicts a hadronic origin of $γ$-ray emission from very old remnants interacting mostly with molecular clouds and a leptonic origin for intermediate age remnants whose shocks propagate in a low density environment created by their progenitors via e.g., strong stellar winds. These results can be regarded as evidence in support of the SNR-origin of the Galactic cosmic rays.

preprint2011arXiv

A statistical model for the $γ$-ray variability of the Crab nebula

A statistical scenario is proposed to explain the $γ$-ray variability and flares of the Crab nebula, which were observed recently by the Fermi/LAT. In this scenario electrons are accelerated in a series of knots, whose sizes follow a power-law distribution. These knots presumably move outwards from the pulsar and have a distribution in the Doppler boost factor. The maximal electron energy is assumed to be proportional to the size of the knot. Fluctuations at the highest energy end of the overall electron distribution will result in variable $γ$-ray emission via the synchrotron process in the $\sim 100$ MeV range. Since highly boosted larger knots are rarer than smaller knots, the model predicts that the variability of the synchrotron emission increases with the photon energy. We realize such a scenario with a Monte-Carlo simulation and find that the model can reproduce both the two $γ$-ray flares over a period of $\sim$ year and the monthly scale $γ$-ray flux fluctuations as observed by the Fermi/LAT. The observed $γ$-ray spectra in both the steady and flaring states are also well reproduced.

preprint2011arXiv

A statistical model to explain the gamma-ray variability and flares of the Crab nebula

Recently the AGILE and Fermi/LAT detectors uncovered giant $γ$-ray flares from the Crab nebula. The duration of these flares is a few days. The Fermi/LAT data with monthly time binning further showed significant variability of the synchrotron tail of the emission, while the inverse Compton component was stable. The simultaneous or follow-up observations in X-ray, optical, infrared and radio bands did not find significant flux variation. Based on these observations, we propose that the $γ$-ray variability and flares are due to statistical fluctuations of knots that can accelerate electrons to $\sim$PeV energies. The maximum achievable energy of electrons is adopted to be proportional to the size of the knot, which is assumed to follow a power-law distribution. Thus the low energy electron flux will be stable due to the large number of small knots, while the high energy electron flux may experience large fluctuations. Monte Carlo realization of such a picture can reproduce the observational data quite well given proper model parameters.

preprint2011arXiv

Derivation of the Electron Distribution in SNR RX J1713.7-3946 via a Spectral Inversion Method

We show that the radio, X-ray and gamma-ray spectrum of the supernova remnant RX J1713.7-3946 can be accounted for with the simplest emission model, where all of these emissions are attributed to a population of relativistic electrons interacting with the cosmic microwave background radiation, IR interstellar photons, and a background magnetic field. With a spectral inversion method (Johns & Lin 1992), the parent electron distribution and its uncertainties are derived from the observed photon spectrum. These results are independent on the model of particle acceleration and strongly support the leptonic scenario for the TeV emission.

preprint2011arXiv

Modeling the Multi-Wavelength Emission of Shell-Type Supernova Remnant RX J1713.7-3946

Emission mechanisms of the shell-type supernova remnant (SNR) RX J1713.7-3946 are studied with multi-wavelength observational data from radio, X-ray, GeV $γ$-ray to TeV $γ$-ray band. A Markov Chain Monte Carlo method is employed to explore the high-dimensional model parameter space systematically. Three scenarios for the $γ$-ray emission are investigated: the leptonic, the hadronic and a hybrid one. Thermal emission from the background plasma is also included to constrain the gas density, assuming ionization equilibrium, and a 2$σ$ upper limit of about 0.03 cm$^{-3}$ is obtained as far as thermal energies account for a significant fraction of the dissipated kinetic energy of the SNR shock. Although systematic errors dominate the $χ^2$ of the spectral fit of all models, we find that 1) the leptonic model has the best constrained model parameters, whose values can be easily accommodated with a typical supernova, but gives relatively poor fit to the TeV $γ$-ray data; 2) The hybrid scenario has one more parameter than the leptonic one and improves the overall spectral fit significantly; 3) The hadronic one, which has three more parameters than the leptonic model, gives the best fit to the overall spectrum with relatively not-well-constrained model parameters and very hard spectra of accelerated particles. The uncertainties of the model parameters decrease significantly if the spectral indices of accelerated electrons and protons are the same. The hybrid and hadronic models also require an energy input into high-energy protons, which seems to be too high compared with typical values of a supernova explosion. Further investigations are required to reconcile these observations with SNR theories.

preprint2010arXiv

General Relativistic Magnetohydrodynamic and Monte Carlo Modeling of Sagittarius A*

We present results of models of the physical space and parameters of the accretion disk of Sagittarius A*, as well as simulations of its emergent spectrum. This begins with HARM, a 2D general relativistic magneto-hydrodynamic (GRMHD) model, specifically set up to evolve the space around a black hole. Data from HARM are then fed into a 2D Monte-Carlo (MC) code which generates and tracks emitted photons, allowing for absorption and scattering before they escape the volume.

preprint2010arXiv

Relationship between Hard and Soft X-ray Emission Components of a Solar Flare

X-ray observations of solar flares routinely reveal an impulsive high-energy and a gradual low-energy emission component, whose relationship is one of the key issues of solar flare study. The gradual and impulsive emission components are believed to be associated with, respectively, the thermal and nonthermal components identified in spectral fitting. In this paper, a prominent about 50 second hard X-ray (HXR) pulse of a simple GOES class C7.5 flare on 20 February 2002 is used to study the association between high energy, non-thermal and impulsive evolution, and low energy, thermal and gradual evolution. We use regularized methods to obtain time derivatives of photon fluxes to quantify the time evolution as a function of photon energy, obtaining a break energy between impulsive and gradual behavior. These break energies are consistent with a constant value of about 11 keV in agreement with those found spectroscopically between thermal and non-thermal components, but the relative errors of the former are greater than 15% and much greater than the a few percent errors found from the spectral fitting. These errors only weakly depend on assuming an underlying spectral model for the photons, pointing to the current data being inadequate to reduce the uncertainties rather than there being a problem associated with an assumed model. The time derivative method is used to test for the presence of a 'pivot energy' in this flare. Although these pivot energies are marginally consistent with a constant value of about 9 keV, its values in the HXR rise phase appear to be lower than those in the decay phase.

preprint2009arXiv

Elementary Energy Release Events in Flaring Loops: Effects of Chromospheric Evaporation on X-rays

With the elementary energy release events introduced in a previous paper (Liu & Fletcher 2009) we model the chromospheric evaporation in flaring loops. The thick-target hard X-ray (HXR) emission produced by electrons escaping from the acceleration region dominates the impulsive phase and the thin-target emission from the acceleration region dominates the low-energy thermal component in the gradual phase, as observed in early impulsive flares. Quantitative details depend on properties of the thermal background, which leads to variations in the correlation between HXR flux and spectral index. For lower temperature and/or higher density of the background electrons, the HXRs both rise and decay more quickly with a plateau near the peak. The plateau is less prominent at higher energies. Given the complexity of transport of mass, momentum, and energy along loops in the impulsive phase, we propose a strategy to apply this single-zone energy release and electron acceleration model to observations of flares associated with single loops so that the energy release, electron acceleration, and evaporation processes may be studied quantitatively.

preprint2009arXiv

MC Simulations of the Broadband Spectra of Sagittarius A* through the use of GRMHD

We present results of simulations of the spectrum of the accretion flow onto the supermassive black hole in our Galactic Centre, Sagittarius A*, generated with a coupling of Monte-Carlo (MC) radiation and general relativistic magnetohydrodynamic (GRMHD) codes. In our modeling, we use the 2D HARM GRMHD code to first model the physical parameters of the disk, then feed its results into our 2D MC photon transport code. We will discuss results obtained which fit radio, IR, and Chandra-obtained flaring or quiescent x-ray data points, as well as the validity of the amount of scaling of input parameters (density, temperature, and magnetic field) required to fit these points. HARM output will be used to suggest whether the scaling is within reasonable limits.

preprint2009arXiv

Polarized Emission of Sagittarius A*

We explore the parameter space of the two temperature pseudo-Newtonian Keplerian accretion flow model for the millimeter and shorter wavelength emission from Sagittarius A*. A general relativistic ray-tracing code is used to treat the radiative transfer of polarized synchrotron emission from the flow. The synchrotron self-Comptonization and bremsstrahlung emission components are also included. It is shown that the model can readily account for the millimeter to sub-millimeter emission characteristics with an accretion rate of ~6x10^17g.s^-1 and an inclination angle of ~40 deg. However, the corresponding model predicted near-infrared and X-ray fluxes are more than one order of magnitude lower than the observed 'quiescent' state values. While the extended quiescent-state X-ray emission has been attributed to thermal emission from the large-scale accretion flow, the NIR emission and flares are likely dominated by emission regions either within the last stable orbit of a Schwarzschild black hole or associated with outflows. With the viscous parameter derived from numerical simulations, there is still a degeneracy between the electron heating rate and the magnetic parameter. A fully general relativistic treatment with the black hole spin incorporated will resolve these issues.

preprint2008arXiv

Stochastic Electron Acceleration in Shell-Type Supernova Remnants

We study the stochastic electron acceleration by fast mode waves in the turbulent downstream of weakly magnetized collisionless astrophysical shocks. The acceleration is most efficient in a dissipative layer, and the model characteristics are determined by the shock speed, density, magnetic field, and turbulence decay length. The model explains observations of shell-type supernova remnants RX J1713.7-3946 and J0852.0-4622 and can be tested by observations in hard X-rays with the HXMT and NuSTAR or gamma-rays with the GLAST.

preprint2006arXiv

Correlation between Flux and Spectral Index during Flares in Sagittarius A*

Flares in Sagittarius A* are produced by hot plasmas within a few Schwarzschild radii of the supermassive black hole at the Galactic center. The recent detection of a correlation between the spectral index and flux during a near infrared (NIR) flare provides a means to conduct detailed investigations of the plasma heating and radiation processes. We study the evolution of the electron distribution function under the influence of a turbulent magnetic field in a hot collisionless plasma. The magnetic field, presumably generated through instabilities in the accretion flow, can both heat the plasma via resonant wave-particle coupling and cool the electrons via radiation. The electron distribution can generally be approximated as relativistic Maxwellian. To account for the observed correlation, we find that the magnetic field needs to be anti-correlated with the electron ''temperature''. NIR and X-ray light curves are produced for a cooling and a heating phase. The model predicts simultaneous flare activity in the NIR and X-ray bands, which can be compared with observations. These results can be applied to MHD simulations to study the radiative characteristics of collisionless plasmas, especially accretion flows in low-luminosity AGNs.

preprint2006arXiv

Probing the Density in the Galactic Center Region: Wind-Blown Bubbles and High-Energy Proton Constraints

Recent observations of the Galactic center in high-energy gamma-rays (above 0.1TeV) have opened up new ways to study this region, from understanding the emission source of these high-energy photons to constraining the environment in which they are formed. We present a revised theoretical density model of the inner 5pc surrounding Sgr A* based on the fact that the underlying structure of this region is dominated by the winds from the Wolf-Rayet stars orbiting Sgr A*. An ideal probe and application of this density structure is this high energy gamma-ray emission. We assume a proton-scattering model for the production of these gamma-rays and then determine first whether such a model is consistent with the observations and second whether we can use these observations to further constrain the density distribution in the Galactic center.

preprint2006arXiv

Stochastic Acceleration in the Galactic Center HESS Source

Stochastic acceleration of electrons interacting resonantly with a turbulent magnetic field in a small accretion torus appears to be the likely mechanism responsible for much of Sagittarius A*'s millimeter and shorter wavelength spectrum. The longer wavelength radiation is produced at larger radii by electrons either diffusing from smaller scales or accelerated in situ. An important prediction of this model is the ejection of a significant flux of relativistic protons from a magnetic-field-dominated acceleration site into the wind-shocked medium surrounding the black hole. Recently, several air Cerenkov telescopes, notably HESS, have detected TeV emission from the Galactic center, with characteristics hinting at a p-p-induced pion decay process for the γ-ray emission. Given (1) the size of this acceleration region measured in the radio band and (2) the wind-injected ISM mapped with Chandra using the diffuse X-rays, it is feasible to test the idea that protons accelerated within \~20 Schwarzschild radii of the black hole produce the TeV emission farther out. We show a fraction of TeV protons scattering about once within ~3 pc of Sagittarius A* and the proton power (~10^37 erg s^-1) produced in concert with the 7 mm radio emission matches the TeV luminosity well. This model explains why the TeV source does not vary on a timescale of a year or less. The particle cascade generated by the p-p scatterings also produces bremsstrahlung, inverse Compton, and synchrotron emission at longer wavelengths from secondary particles. We compare these with current measurements and demonstrate that GLAST will detect this source during its one-year all-sky survey.

preprint2005arXiv

Stochastic Acceleration of $^3$He and $^4$He in Solar Flares by Parallel Propagating Plasma Waves: General Results

We study the acceleration in solar flares of $^3$He and $^4$He from a thermal background by parallel propagating plasma waves with a general broken power-law spectrum that takes into account the turbulence generation processes at large scales and the thermal damping effects at small scales. The exact dispersion relation for a cold plasma is used to describe the relevant wave modes. Because low-energy $α$-particles only interact with small scale waves in the $^4$He-cyclotron branch, where the wave frequencies are below the $α$-particle gyro-frequency, their pitch angle averaged acceleration time is at least one order of magnitude longer than that of $^3$He ions, which mostly resonate with relatively higher frequency waves in the proton-cyclotron (PC) branch. The $α$-particle acceleration rate starts to approach that of $^3$He beyond a few tens of keV nucleon$^{-1}$, where $α$-particles can also interact with long wavelength waves in the PC branch. However, the $^4$He acceleration rate is always smaller than that of $^3$He. Consequently, the acceleration of $^4$He is suppressed significantly at low energies, and the spectrum of the accelerated $α$-particles is always softer than that of $^3$He. The model gives reasonable account of the observed low-energy $^3$He and $^4$He fluxes and spectra in the impulsive solar energetic particle events observed with the {\it Advanced Composition Explorer}. We explore the model parameter space to show how observations may be used to constrain the model.

preprint2004arXiv

Stochastic Acceleration of 3He and 4He by Parallel Propagating Plasma Waves

Stochastic acceleration of $^3$He and $^4$He from a thermal background by parallel propagating turbulent plasma waves with a single power-law spectrum of the wavenumber is studied. In the model, both ions interact with several resonant modes. When one of these modes dominates, the acceleration rate is reduced considerably. At low energies, this happens for $^4$He, but not for $^3$He where contributions from the two stronger modes are comparable so that acceleration of $^3$He is very efficient. As a result, the acceleration of $^4$He is suppressed by a barrier below $\sim 100$ keV nucleon$^{-1}$ and there is a prominent quasi-thermal component in the $^4$He spectra, while almost all the injected $^3$He ions are accelerated to high energies. This accounts for the large enrichment of $^3$He at high energies observed in impulsive solar energetic particle events. With reasonable plasma parameters this also provides a good fit to the spectra of both ions. Beyond $\sim 1$ MeV nucleon$^{-1}$, the spectrum of $^3$He is softer than that of $^4$He, which is consistent with the observed decrease of the $^3$He to $^4$He ratio with energy. This study also indicates that the acceleration, Coulomb losses and diffusive escape of the particles from the acceleration site {\it all} play important roles in shaping the ion spectra. This can explain the varied spectral shapes observed recently by the {\it Advanced Composition Explorer}.

preprint2004arXiv

Stochastic Acceleration of Electrons and Protons. I. Acceleration by Parallel Propagating Waves

Stochastic acceleration of electrons and protons by waves propagating parallel to the large scale magnetic fields of magnetized plasmas is studied with emphasis on the feasibility of accelerating particles from a thermal background to relativistic energies and with the aim of determining the relative acceleration of the two species in one source. In general, the stochastic acceleration by these waves results in two distinct components in the particle distributions, a quasi-thermal and a hard nonthermal, with the nonthermal one being more prominent in hotter plasmas and/or with higher level turbulence. This can explain many of the observed features of solar flares. Regarding the proton to electron ratio, we include the effects of $^4$He in the dispersion relation and the damping of the turbulent waves by the thermal background plasma. The relative acceleration of protons and electrons is very sensitive to the ratio of the electron plasma frequency to its gyro-frequency. Protons are preferentially accelerated in weakly magnetized plasmas. The formalism developed here is applicable to the acceleration of other ion species and to other astrophysical systems.

preprint2001arXiv

Accretion Processes in the Nucleus of M31

The hypothesized supermassive black hole in the nucleus of M31 (which we shall hereafter call M31*) has many features in common with Sgr A* at the Galactic Center, yet they differ in several significant and important ways. Though M31* is probably ten times heavier, its radio luminosity at 3.6 cm is only one third that of Sgr A*. At the same time, M31* is apparently thousands of times more luminous in X-rays than its Galactic Center counterpart. Thus, a comparative study of these objects can be valuable in helping us to understand the underlying physical basis for their activity. We show here that the accretion model being developed for Sgr A* comprises two branches of solutions, distinguished by the relative importance of cooling compared to compressional heating at the radius $r_C$ where the ambient gas is captured by the black hole. For typical conditions in the ISM, the initial temperature ($T[r_C]\sim 10^6-10^7$ K) sits on the unstable branch of the cooling function. Depending on the actual value of $T(r_C)$ and the accretion rate, the plasma settles either onto a hot branch (attaining a temperature as high as $10^{10} K$ or so at small radii) or a cold branch, in which $T$ drops to $\sim 10^4$ K. Sgr A* is presumably a `hot' black hole. We show here that the VLA, UV and {\it Chandra} observations of M31* reveal it to be a member of the `cold' black hole family. We discuss several predicted features in the spectrum of M31* that may be testable by future multi-wavelength observations, including the presence of a prominent UV spike (from hydrogen line emission) that would be absent on the hot branch.