Source author record

Edison Liang

Edison Liang 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

23works
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

23 published item(s)

preprint2022arXiv

Research Opportunities in Plasma Astrophysics

Major scientific questions and research opportunities are described on 10 unprioritized plasma astrophysics topics: (1) magnetic reconnection, (2) collisionless shocks and particle acceleration, (3) waves and turbulence, (4) magnetic dynamos, (5) interface and shear instabilities, (6) angular momentum transport, (7) dusty plasmas, (8) radiative hydrodynamics, (9) relativistic, pair-dominated and strongly magnetized plasmas, (10) jets and outflows. Note that this is a conference report from a Workshop on Opportunities in Plasma Astrophysics (WOPA, https://w3.pppl.gov/conferences/2010/WOPA/) in January 2010, that attracted broad representation from the community and was supported by the U.S. Department of Energy, National Aeronautics and Space Administration, National Science Foundation, American Physical Society's Topical Group for Plasma Astrophysics and Division of Plasma Physics, and Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas. Although there has been much planning and many developments in both science and infrastructure since the report was written, most of the motivation, priorities, problems and technical challenges discussed therein remain unaddressed and are relevant at the time of posting.

preprint2020arXiv

Fermi-type particle acceleration from magnetic reconnection at the termination shock of a relativistic striped wind

An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell (PIC) simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultra-relativistic regime with bulk Lorentz factor up to γ_{0}=10^{6}. Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. We find that the resulting energy spectra crucially depend on λ/d_{e}. When λ/d_{e} is large (λ\gtrsim40d_{e}) , the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime with upstream magnetization parameter σ_{0}=10. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as γ_{0} and σ_{0} increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor γ_{0}\approx10^{6} and magnetization parameter σ_{0}=10, which can explain the recent observations of high-energy gamma-rays from pulsar wind nebulae (PWNe).

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.

preprint2020arXiv

Studying particle acceleration from driven magnetic reconnection at the termination shock of a relativistic striped wind using particle-in-cell simulations

A rotating pulsar creates a surrounding pulsar wind nebula (PWN) by steadily releasing an energetic wind into the interior of the expanding shockwave of supernova remnant or interstellar medium. At the termination shock of a PWN, the Poynting-flux-dominated relativistic striped wind is compressed. Magnetic reconnection is driven by the compression and converts magnetic energy into particle kinetic energy and accelerating particles to high energies. We carrying out particle-in-cell (PIC) simulations to study the shock structure as well as the energy conversion and particle acceleration mechanism. By analyzing particle trajectories, we find that many particles are accelerated by Fermi-type mechanism. The maximum energy for electrons and positrons can reach hundreds of TeV.

preprint2020arXiv

Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

preprint2019arXiv

Modeling hydrodynamics, magnetic fields and synthetic radiographs for high-energy-density plasma flows in shock-shear targets

Three-dimensional FLASH radiation-magnetohydrodynamics (radiation-MHD) modeling is carried out to study the hydrodynamics and magnetic fields in the shock-shear derived platform. Simulations indicate that fields of tens of Tesla can be generated via Biermann battery effect due to vortices and mix in the counter-propagating shock-induced shear layer. Synthetic proton radiography simulations using MPRAD and synthetic X-ray image simulations using SPECT3D are carried out to predict the observable features in the diagnostics. Quantifying the effects of magnetic fields in inertial confinement fusion (ICF) and high-energy-density (HED) plasmas represents frontier research that has far-reaching implications in basic and applied sciences.

preprint2019arXiv

MPRAD: A Monte Carlo and ray-tracing code for the proton radiography in high-energy-density plasma experiments

Proton radiography is used in various high-energy-density (HED) plasma experiments. In this paper, we describe a Monte Carlo and ray-tracing simulation tool called MPRAD that can be used for modeling the deflection of proton beams in arbitrary three dimensional electromagnetic fields, as well as the diffusion of the proton beams by Coulomb scattering and stopping power. The Coulomb scattering and stopping power models in cold matter and fully ionized plasma are combined using interpolation. We discuss the application of MPRAD in a few setups relevant to HED plasma experiments where the plasma density can play a role in diffusing the proton beams and affecting the prediction and interpretation of the proton images. It is shown how the diffusion due to plasma density can affect the resolution and dynamical range of the proton radiography.

preprint2014arXiv

Effects of dust feedback on vortices in protoplanetary disks

We carried out two-dimensional high-resolution simulations to study the effect of dust feedback on the evolution of vortices induced by massive planets in protoplanetary disks. Various initial dust to gas disk surface density ratios ($0.001$ -- $0.01$) and dust particle sizes (Stokes number $4\times10^{-4}$ -- $0.16$) are considered. We found that while dust particles migrate inwards, vortices are very effective in collecting them. When dust density becomes comparable to gas density within the vortex, a dynamical instability is excited and it alters the coherent vorticity pattern and destroys the vortex. This dust feedback effect is stronger with higher initial dust/gas density ratio and larger dust grain. Consequently, we found that the disk vortex lifetime can be reduced up to a factor of 10. We discuss the implications of our findings on the survivability of vortices in protoplanetary disks and planet formation.

preprint2014arXiv

Ultra-Intense Gamma-Rays Created Using the Texas Petawatt Laser

In a series of experiments at the Texas Petawatt Laser (TPW) in Austin, TX, we have used attenuation spectrometers, dosimeters, and a new Forward Compton Electron Spectrometer (FCES) to measure and characterize the angular distribution, fluence, and energy spectrum of the X-rays and gamma rays produced by the TPW striking multi-millimeter thick gold targets. Our results represent the first such measurements at laser intensities > 10 21 W*cm-2 and pulse durations < 150 fs. We obtain a maximum yield of X-ray and gamma ray energy with respect to laser energy of 4% and a mean yield of 2%. We futher obtain a Full Width Half Maximum (FWHM) of the gamma distribution of 37°. We were able to characterize the gamma-ray spectrum from 3 MeV to 90 MeV using a Forward Compton Electron Spectrometer, with an energy resolution of 0.5 MeV and mean kT of ~ 6 MeV.. We were able to characterize the spectrum from 1 to 7 MeV using a Filter Stack Spectrometer, measuring a mean gamma-ray temperature for the spectrum from 3 to 7 MeV of 2.1 MeV.

preprint2013arXiv

Hot Electron and Pair Production from the Texas Petawatt Laser Irradiating Thick Gold Targets

We present data for relativistic hot electron production by the Texas Petawatt Laser irradiating solid Au targets with thickness between 1 and 4 mm. The experiment was performed at the short focus target chamber TC1 in July 2011, with laser energies around 50 J. We measured hot electron spectra out to 50 MeV which show a narrow peak around 10 - 20 MeV plus high energy exponential tail. The hot electron spectral shape differs from those reported for other PW lasers. We did not observe direct evidence of positron production above background.

preprint2013arXiv

Magnetic Field Generation and Particle Energization at Relativistic Shear Boundaries in Collisionless Electron-Positron Plasmas

Using 2.5-dimensional Particle-in-Cell simulations, we study the kinetic physics of relativistic shear flow boundary in collisionless electron-positron (e+e-) plasmas. We find efficient magnetic field generation and particle energization at the shear boundary, driven by streaming instabilities across the shear interface and sustained by the shear flow. Nonthermal, anisotropic high-energy particles are accelerated across field lines to produce a power-law tail, truncated at energies below the shear Lorentz factor. These results have important implications for the dissipation and radiation of jets in blazars, gamma-ray bursts and other relativistic outflows.

preprint2013arXiv

Relativistic Positron-Electron-Ion Shear Flows and Application to Gamma-Ray Bursts

We present Particle-in-Cell simulation results of relativistic shear flows for hybrid positron-electron-ion plasmas and compare to those for pure e+e- and pure e-ion plasmas. Among the three types of relativistic shear flows, we find that only hybrid shear flow is able to energize the electrons to form a high-energy spectral peak plus a hard power-law tail. Such electron spectra are needed to model the observational properties of gamma-ray bursts.

preprint2012arXiv

Electrostatic field acceleration of laser-driven ion bunch by using double layer thin foils

Monoenergetic ion bunch generation and acceleration from double layer thin foil target irradiated by intense linearly polarized (LP) laser pulse is investigated using two-dimensional (2D) particle-in-cell (PIC) simulations. The low-Z ions in the front layer of the target are accelerated by the laser-driven hot electrons and penetrate through the high-Z ion layer to generate a quasi-monoenergetic ion bunch, and this bunch will continue to be accelerated by the quasi-stable electrostatic sheath field which is formed by the immobile high-Z ions and the hot electrons. This mechanism offers possibility to generate monoenergetic ion bunch without ultrahigh-contrast and ultrahigh gradient laser pulses in beam generation experiments, which is confirmed by our simulations.

preprint2012arXiv

Time-dependent simulations of emission from FSRQ PKS1510-089: multiwavelength variability of external Compton and SSC models

[abridged] We present results of modeling the SED and multiwavelength variability of the bright FSRQ PKS1510-089 with our time-dependent multizone Monte Carlo/Fokker-Planck code (Chen et al. 2001). As primary source of seed photons for inverse Compton scattering, we consider radiation from the broad line region (BLR), from the molecular torus, and the local synchrotron radiation (SSC). Different scenarios are assessed by comparing simulated light curves and SEDs with one of the best flares by PKS1510-089, in March 2009. The time-dependence of our code and its correct handling of light travel time effects allow us to fully take into account the effect of the finite size of the active region, and in turn to fully exploit the information carried by time resolved observed SEDs, increasingly available since the launch of Fermi. We confirm that the spectrum adopted for the external radiation has an important impact on the modeling of the SED, in particular for the lower energy end of the Compton component, observed in the X-ray band, which in turn is one of the most critical bands to assess the differences between EC and SSC emission. In the context of the scenario presented here, where the flaring is caused by the increase of the number of relativistic electrons ascribed to the effect of the interaction of a portion of the jet (blob) with a shock, we can not firmly discriminate the three main scenarios for gamma-ray emission. However, results show clearly the differences produced by a more realistic treatment of the emitting source in the shape of SEDs and their time variability over relevant, observable time-scales, and demonstrate the crucial importance of time-dependent multi-zone models to advance our understanding of the physics of these sources, by taking full advantage of the wealth of information offered by the high quality data of current multiwavelength campaigns.

preprint2011arXiv

Numerical Modeling of Multi-wavelength Spectra of M87 Core Emission

Spectral fits to M87 core data from radio to hard x-ray are generated via a specially selected software suite, comprised of the HARM GRMHD accretion disk model and a 2D Monte Carlo radiation transport code. By determining appropriate parameter changes necessary to fit x-ray quiescent and flaring behavior of M87's core, we assess the reasonableness of various flaring mechanisms. This shows that an accretion disk model of M87's core out to 28 GM/c^2 can describe the inner emissions. High spin rates show GRMHD-driven polar outflow generation, without citing an external jet model. Our results favor accretion rate changes as the dominant mechanism of x-ray flux and index changes, with variations in density of approximately 20% necessary to scale between the average x-ray spectrum and flaring or quiescent spectra. The best fit parameters are black hole spin a/M > 0.8 and maximum accretion flow density n <= 3x10^7 cm^-3, equivalent to horizon accretion rates between m_dot = M_dot/M_dot_Edd ~ 2x10^-6 and 1x10^-5 (with M_dot_Edd defined assuming a radiative efficiency eta = 0.1). These results demonstrate that the immediate surroundings of M87's core are appropriate to explain observed x-ray variability.

preprint2011arXiv

Scaling of Laser Produced Pair Plasmas And Comparison With Experimental Data

We report simulation results of pair production by ultra-intense lasers irradiating a gold target using the GEANT4 Monte-Carlo code. Certain experimental features of the positron and electron energy spectra are reproduced, as well as trends with regard to target thickness and hot electron temperature Te. For Te in the range 5-10 MeV, the optimal target thickness for pair production is found to be about 3 mm. Further Monte-Carlo simulations may aid in the optimization of laser-driven positron sources.

preprint2011arXiv

Time dependent simulations of multiwavelength variability of the blazar Mrk 421 with a Monte Carlo multi-zone code

(abridged) We present a new time-dependent multi-zone radiative transfer code and its application to study the SSC emission of Mrk 421. The code couples Fokker-Planck and Monte Carlo methods, in a 2D geometry. For the first time all the light travel time effects (LCTE) are fully considered, along with a proper treatment of Compton cooling, which depends on them. We study a set of simple scenarios where the variability is produced by injection of relativistic electrons as a `shock front' crosses the emission region. We consider emission from two components, with the second one either being pre-existing and co-spatial and participating in the evolution of the active region, or spatially separated and independent, only diluting the observed variability. Temporal and spectral results of the simulation are compared to the multiwavelength observations of Mrk 421 in March 2001. We find parameters that can adequately fit the observed SEDs and multiwavelength light curves and correlations. There remain however a few open issues, most notably: i) systematic soft intra-band X-ray lags. ii) The quadratic correlation between the TeV and X-ray flux during the flare decay has not been reproduced. These features are among those affected by the spatial extent and geometry of the source. The difficulty of producing hard X-ray lags is exacerbated by a bias towards soft lags caused by the combination of energy dependent radiative cooling time-scales and LCTE. About the second emission component, our results strongly favor the scenario where it is co-spatial and it participates in the flare evolution, suggesting that different phases of activity may occur in the same region. The cases presented in this paper represent only an initial study, and despite their limited scope they make a strong case for the need of true time-dependent and multi-zone modeling.

preprint2011arXiv

X-ray Time Lags in TeV Blazars

We use Monte Carlo/Fokker-Planck simulations to study the X-ray time lags. Our results show that soft lags will be observed as long as the decay of the flare is dominated by radiative cooling, even when acceleration and cooling timescales are similar. Hard lags can be produced in presence of a competitive achromatic particle energy loss mechanism if the acceleration process operates on a timescale such that particles are slowly moved towards higher energy while the flare evolves. In this type of scenario, the γ -ray/X-ray quadratic relation is also reproduced.

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.

preprint2009arXiv

Electron heating and acceleration by magnetic reconnection in hot accretion flows

Both analytical and numerical works show that magnetic reconnection must occur in hot accretion flows. This process will effectively heat and accelerate electrons. In this paper we use the numerical hybrid simulation of magnetic reconnection plus test-electron method to investigate the electron acceleration and heating due to magnetic reconnection in hot accretion flows. We consider fiducial values of density, temperature, and magnetic parameter $β_e$ (defined as the ratio of the electron pressure to the magnetic pressure) of the accretion flow as $n_{0} \sim 10^{6} {\rm cm^{-3}}$, $T_{e}^0\sim 2\times 10^9 {\rm K}$, and $β_e=1$. We find that electrons are heated to a higher temperature $T_{e}=5\times 10^9$K, and a fraction $η\sim 8%$ of electrons are accelerated into a broken power-law distribution, $dN(γ)\propto γ^{-p}$, with $p\approx 1.5$ and 4 below and above $\sim 1$ MeV, respectively. We also investigate the effect of varying $β$ and $n_0$. We find that when $β_e$ is smaller or $n_0$ is larger, i.e, the magnetic field is stronger, $T_e$, $η$, and $p$ all become larger.

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.

preprint2006arXiv

Sustained Comoving Acceleration of Overdense Plasmas by Colliding Laser Pulses

We review recent PIC simulation results which show that double-sided irradiation of a thin over- dense plasma slab by ultra-intense laser pulses from both sides can lead to sustained comoving acceleration of surface electrons to energies much higher than the conventional ponderomotive limit. The acceleration stops only when the electrons drift transversely out of the laser beam. We show latest 2.5D results of parameter studies based on finite laser spot size and discuss future laser experiments that can be used to test these results.