Mechanical active galactic nucleus (AGN) feedback plays a key role in massive galaxies, galaxy groups and clusters. However, the energy content of AGN jets that mediate this feedback process is still far from clear. Here we present a preliminary study of radial elongations $\tau$ of a large sample of X-ray cavities, which are apparently produced by mechanical AGN feedback. All the cavities in our sample are elongated along the angular (type-I) or jet directions (type-II), or nearly circular (type-III). The observed value of $\tau$ roughly decreases as the cavities rise buoyantly, confirming the same trend found in hydrodynamic simulations. For young cavities, both type-I and II cavities exist, and the latter dominates. Assuming a spheroidal cavity shape, we derive an analytical relation between the intrinsic radial elongation $\bar{\tau}$ and the inclination-angle-dependent value of $\tau$, showing that projection effect makes cavities appear more circular, but does not change type-I cavities into type-II ones, or vice versa. We summarize radial elongations of young cavities in hydrodynamic simulations, and find that $\bar{\tau}$ increases with the kinetic fraction of AGN jets. While mild jets always produce type-II cavities, thermal-energy-dominated strong jets produce type-I cavities, and kinetic-energy-dominated strong jets produce type-II cavities. The existence of type-I young cavities indicates that some AGN jets are strong and dominated by thermal energy (or cosmic rays). If most jets are dominated by non-kinetic energies, our results suggest that they must be long-duration mild jets. However, if most jets are strong, they must be dominated by the kinetic energy.

Magnetic reconnection in the relativistic and trans-relativistic regimes is able to accelerate particles to hard power law energy spectra $f \propto \gamma^{-p}$ (approaching $p=1$). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during magnetic reconnection in the trans-relativistic regime of a proton-electron plasma. While earlier work in this parameter regime has focused on the effects of electric field parallel to the local magnetic field on the particle injection (from thermal energy to the lower energy bound of the power-law spectrum), here we examine the roles of both parallel and perpendicular electric fields to gain a more complete understanding on the injection process and further development of a power-law spectrum. We show that the parallel electric field does contribute significantly to particle injection, and is more important in the initial phase of magnetic reconnection. However, as the simulation proceeds, the acceleration by the perpendicular electric field becomes more important for particle injection and completely dominates the acceleration responsible for the high-energy power-law spectrum. This holds robustly, in particular for longer reconnection times and larger systems, i.e. in simulations that are more indicative of the processes in astrophysical sources.

Magnetic reconnection in strongly magnetized astrophysical plasma environments is believed to be the primary process for fast energy release and particle energization. Currently there is strong interest in relativistic magnetic reconnection, in that it may provide a new explanation for high-energy particle acceleration and radiation in strongly magnetized astrophysical systems. We review recent advances in particle acceleration and reconnection physics in the magnetically-dominated regime. More discussion is focused on the physics of particle acceleration, power-law formation as well as the reconnection rate problem. In addition, we provide an outlook for studying reconnection acceleration mechanisms and kinetic physics in the next step.

The Milky Way (MW) mass $M_{\rm vir}$ is a fundamental quantity in astronomy. Although it has been measured extensively, it is still uncertain to more than a factor of two due to limited number or spatial coverage of kinematic tracers. Here we use a novel method to constrain $M_{\rm vir}$ based on the properties of the MW corona. We build a hydrostatic corona model with non-thermal pressure support and a physically-motivated density profile, and derive the temperature distribution, which depends on $M_{\rm vir}$. While the temperature profile decreases substantially with radius, the X-ray-emission-weighted average temperature $T_{\rm em}$ is quite uniform toward different sight lines, consistent with X-ray observations. Using available measurements of $T_{\rm em}$, we find that $M_{\rm vir}$=$1.60^{+1.35}_{-0.41}$ $\times$ $10^{12}$ $M_{\rm sun}$ assuming an Navarro-Frenk-White (NFW) total matter distribution. This estimate is independent of the uncertain total corona mass and the gas temperature at very large radii, and is on the high mass side of current measurements. Adopting a $\beta$ model for the corona density distribution or a total matter distribution contributed by an NFW dark matter distribution and a central cold baryonic matter distribution leads to similar estimates of $M_{\rm vir}$. Non-thermal pressure support, which likely exists in the corona, leads to higher values of $M_{\rm vir}$.

In the standard model of solar flares, a large-scale reconnection current sheet is postulated as the central engine for powering the flare energy release and accelerating particles. However, where and how the energy release and particle acceleration occur remain unclear due to the lack of measurements for the magnetic properties of the current sheet. Here we report the measurement of spatially-resolved magnetic field and flare-accelerated relativistic electrons along a current-sheet feature in a solar flare. The measured magnetic field profile shows a local maximum where the reconnecting field lines of opposite polarities closely approach each other, known as the reconnection X point. The measurements also reveal a local minimum near the bottom of the current sheet above the flare loop-top, referred to as a "magnetic bottle". This spatial structure agrees with theoretical predictions and numerical modeling results. A strong reconnection electric field of ~4000 V/m is inferred near the X point. This location, however, shows a local depletion of microwave-emitting relativistic electrons. These electrons concentrate instead at or near the magnetic bottle structure, where more than 99% of them reside at each instant. Our observations suggest that the loop-top magnetic bottle is likely the primary site for accelerating and/or confining the relativistic electrons.

Active galactic nucleus (AGN) jets are believed to be important in solving the cooling flow problem in the intracluster medium (ICM), while the detailed mechanism is still in debate. Here we present a systematic study on the energy coupling efficiency $\eta_{\rm cp}$, the fraction of AGN jet energy transferred to the ICM. We first estimate the values of $\eta_{\rm cp}$ analytically in two extreme cases, which are further confirmed and extended with a parameter study of spherical outbursts in a uniform medium using hydrodynamic simulations. We find that $\eta_{\rm cp}$ increases from $\sim 0.4$ for a weak isobaric injection to $\gtrsim 0.8$ for a powerful point injection. For any given outburst energy, we find two characteristic outburst powers that separate these two extreme cases. We then investigate the energy coupling efficiency of AGN jet outbursts in a realistic ICM with hydrodynamic simulations, finding that jet outbursts are intrinsically different from spherical outbursts. For both powerful and weak jet outbursts, $\eta_{\rm cp}$ is typically around $0.7-0.9$, partly due to the non-spherical nature of jet outbursts, which produce backflows emanating from the hotspots, significantly enhancing the ejecta-ICM interaction. While for powerful outbursts a dominant fraction of the energy transferred from the jet to the ICM is dissipated by shocks, shock dissipation only accounts for $\lesssim 30\%$ of the injected jet energy for weak outbursts. While both powerful and weak outbursts could efficiently heat cooling flows, powerful thermal-energy-dominated jets are most effective in delaying the onset of the central cooling catastrophe.

The Fermi bubbles are two giant bubbles in gamma rays lying above and below the Galactic center (GC). Despite numerous studies on the bubbles, their origin and emission mechanism remain elusive. Here we use a suite of hydrodynamic simulations to study the scenario where the cosmic rays (CRs) in the bubbles are mainly accelerated at the forward shocks driven by a pair of opposing jets from Sgr A*. We find that an active galactic nucleus (AGN) jet event happened $5-6$ Myr ago can naturally reproduce the bilobular morphology of the bubbles, and the postshock gas temperature in the bubbles is heated to $\sim0.4$ keV, consistent with recent X-ray observations. The forward shocks compress the hot halo gas, and at low latitudes, the compressed gas shows an X-shaped structure, naturally explaining the biconical X-ray structure in the ROSAT 1.5 keV map in both morphology and X-ray surface brightness. CR acceleration is most efficient in the head regions of the bubbles during the first 2 Myrs. The opposing jets release a total energy of $\sim 10^{55}$ erg with an Eddington ratio of $\sim 10^{-3}$, which falls well in the range of the hot accretion flow mode for black holes. Our simulations further show that the forward shocks driven by spherical winds at the GC typically produce bubbles with much wider bases than observed, and could not reproduce the biconical X-ray structure at low latitudes. This suggests that starburst or AGN winds are unlikely the origin of the bubbles in the shock scenario.

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.

Theoretical and observational arguments suggest that there is a large amount of hot ($\sim 10^6$ K), diffuse gas residing in the Milky Way's halo, while its total mass and spatial distribution are still unclear. In this work, we present a general model for the gas density distribution in the Galactic halo, and investigate the gas evolution under radiative cooling with a series of 2D hydrodynamic simulations. We find that the mass inflow rate in the developed cooling flow increases with gas metallicity and the total gas mass in the halo. For a fixed halo gas mass, the spatial gas distribution affects the onset time of the cooling catastrophe, which starts earlier when the gas distribution is more centrally-peaked, but does not substantially affect the final mass inflow rate. The gravity from the Galactic bulge and disk affects gas properties in inner regions, but has little effect on the final inflow rate either. We confirm our results by investigating cooling flows in several density models adopted from the literature, including the Navarro-Frenk-White (NFW) model, the cored-NFW model, the Maller & Bullock model, and the $\beta$ model. Typical mass inflow rates in our simulations range from $\sim 5 M_{\odot}$ yr$^{-1}$ to $\sim 60 M_{\odot}$ yr$^{-1}$, and are much higher than the observed star formation rate in our Galaxy, suggesting that stellar and active galactic nucleus feedback processes may play important roles in the evolution of the Milky Way (MW) and MW-type galaxies.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

We derive the axion-nucleon interaction Lagrangian in heavy baryon chiral perturbation theory up to next-to-next-to-leading order. The effective axion-nucleon coupling is calculated to a few percent accuracy.

It is well known that collisionless systems are dissipation free from the perspective of particle collision and thus conserve entropy. On the other hand, processes such as magnetic reconnection and turbulence appear to convert large-scale magnetic energy into heat. In this paper, we investigate the energization and heating of collisionless plasma. The dissipation process is discussed in terms of fluid entropy in both isotropic and gyrotropic forms. Evolution equations for the entropy are derived and they reveal mechanisms that lead to changes in fluid entropy. These equations are verified by a collisionless particle-in-cell simulation of multiple reconnecting current sheets. In addition to previous findings regarding the pressure tensor, we emphasize the role of heat flux in the dissipation process.

Nonthermal loop-top sources in solar flares are the most prominent observational signature that suggests energy release and particle acceleration in the solar corona. Although several scenarios for particle acceleration have been proposed, the origin of the loop-top sources remains unclear. Here we present a model that combines a large-scale magnetohydrodynamic simulation of a two-ribbon flare with a particle acceleration and transport model for investigating electron acceleration by a fast-mode termination shock at the looptop. Our model provides spatially resolved electron distribution that evolves in response to the dynamic flare geometry. We find a concave-downward magnetic structure located below the flare termination shock, induced by the fast reconnection downflows. It acts as a magnetic trap to confine the electrons at the looptop for an extended period of time. The electrons are energized significantly as they cross the shock front, and eventually build up a power-law energy spectrum extending to hundreds of keV. We suggest that this particle acceleration and transport scenario driven by a flare termination shock is a viable interpretation for the observed nonthermal loop-top sources.

During magnetically dominated relativistic reconnection, inflowing plasma depletes the initial relativistic pressure at the x-line and collisionless plasma heating inside the diffusion region is insufficient to overcome this loss. The resulting pressure drop causes a collapse at the x-line, essentially a localization mechanism of the diffusion region necessary for fast reconnection. The extension of this low-pressure region further explains the bursty nature of anti-parallel reconnection because a once opened outflow exhaust can also collapse, which repeatedly triggers secondary tearing islands. However, a stable single x-line reconnection can be achieved when an external guide field exists, since the reconnecting magnetic field component rotates out of the reconnection plane at outflows, providing additional magnetic pressure to sustain the opened exhausts.

Enhancement of minor ions such as $^3$He and heavy ions in flare-associated solar energetic particle (SEP) events remains one of the major puzzles in heliophysics. In this work, we use 3D hybrid simulations (kinetic protons and fluid electrons) to investigate particle energization in a turbulent low-beta environment similar to solar flares. It is shown that in this regime the injected large-amplitude Alfv\'en waves develop into compressible and anisotropic turbulence, which efficiently heats thermal ions of different species. We find that temperature increase of heavy ions is inversely proportional to the charge to mass ratio, which is consistent with observations of impulsive SEP events. Further analysis reveals that ions are energized by interacting with nearly perpendicular magnetosonic waves near proton inertial scale.

Solar flare termination shocks have been suggested as one of the promising drivers for particle acceleration in solar flares, yet observational evidence remains rare. By utilizing radio dynamic spectroscopic imaging of decimetric stochastic spike bursts in an eruptive flare, Chen et al. found that the bursts form a dynamic surface-like feature located at the ending points of fast plasma downflows above the looptop, interpreted as a flare termination shock. One piece of observational evidence that strongly supports the termination shock interpretation is the occasional split of the emission band into two finer lanes in frequency, similar to the split-band feature seen in fast-coronal-shock-driven type II radio bursts. Here we perform spatially, spectrally, and temporally resolved analysis of the split-band feature of the flare termination shock event. We find that the ensemble of the radio centroids from the two split-band lanes each outlines a nearly co-spatial surface. The high-frequency lane is located slightly below its low frequency counterpart by ~0.8 Mm, which strongly supports the shock upstream-downstream interpretation. Under this scenario, the density compression ratio across the shock front can be inferred from the frequency split, which implies a shock with a Mach number of up to 2.0. Further, the spatiotemporal evolution of the density compression along the shock front agrees favorably with results from magnetohydrodynamics simulations. We conclude that the detailed variations of the shock compression ratio may be due to the impact of dynamic plasma structures in the reconnection outflows, which results in distortion of the shock front.

While observations have suggested that power-law electron energy spectra are a common outcome of strong energy release during magnetic reconnection, e.g., in solar flares, kinetic simulations have not been able to provide definite evidence of power-laws in energy spectra of non-relativistic reconnection. By means of 3D large-scale fully kinetic simulations, we study the formation of power-law electron energy spectra in non-relativistic low-$\beta$ reconnection. We find that both the global spectrum integrated over the entire domain and local spectra within individual regions of the reconnection layer have power-law tails with a spectral index $p \sim 4$ in the 3D simulation, which persist throughout the non-linear reconnection phase until saturation. In contrast, the spectrum in the 2D simulation rapidly evolves and quickly becomes soft. We show that 3D effects such as self-generated turbulence and chaotic magnetic field lines enable the transport of high-energy electrons across the reconnection layer and allow them to access several main acceleration regions. This leads to a sustained and nearly constant acceleration rate for electrons at different energies. We construct a model that explains the observed power-law spectral index in terms of the dynamical balance between particle acceleration and escape from main acceleration regions, which are defined based upon a threshold for the curvature drift acceleration term. This result could be important for explaining the formation of power-law energy spectrum in solar flares.

We present numerical modelling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation with spatial diffusion both along and across the magnetic field. We show that the location on the shock where the high-energy particle intensity is the largest, depends on the energy of the particles and on time. The acceleration of particles to more than 100 MeV mainly occurs in the shock-streamer interaction region, due to perpendicular shock geometry and the trapping effect of closed magnetic fields. A comparison of the particle spectra to that in a radial magnetic field shows that the intensity at 100 MeV (200 MeV) is enhanced by more than one order (two orders) of magnitude. This indicates that the streamer-like magnetic field can be an important factor in producing large solar energetic particle events. We also show that the energy spectrum integrated over the simulation domain consists of two different power laws. Further analysis suggests that it may be a mixture of two distinct populations accelerated in the streamer and open field regions, where the acceleration rate differs substantially. Our calculations also show that the particle spectra are affected considerably by a number of parameters, such as the streamer tilt angle, particle spatial diffusion coefficient, and shock compression ratio. While the low-energy spectra agree well with standard diffusive shock acceleration theory, the break energy ranges from $\sim$1 MeV to $\sim$90 MeV and the high-energy spectra can extend to $\sim$1 GeV with a slope of $\sim$2-3.

While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, since it opens up the possiblity of modeling the energetic spectra without resolving microscopic diffusion regions.

By means of fully kinetic particle-in-cell simulations, we study whether the proton-to-electron mass ratio $m_i/m_e$ influences the energy spectrum and underlying acceleration mechanism during magnetic reconnection. While kinetic simulations are essential for studying particle acceleration during magnetic reconnection, a reduced $m_i/m_e$ is often used to alleviate the demanding computing resources, which leads to artificial scale separation between electron and proton scales. Recent kinetic simulations with high-mass-ratio have suggested new regimes of reconnection, as electron pressure anisotropy develops in the exhaust region and supports extended current layers. In this work, we study whether different $m_i/m_e$ changes the particle acceleration processes by performing a series of simulations with different mass ratio ($m_i/m_e=25-400$) and guide-field strength in a low-$\beta$ plasma. We find that mass ratio does not strongly influence reconnection rate, magnetic energy conversion, ion internal energy gain, plasma energization processes, ion energy spectra, and the acceleration mechanisms for high-energy ions. Simulations with different mass ratios are different in electron acceleration processes, including electron internal energy gain, electron energy spectrum and the acceleration efficiencies for high-energy electrons. We find that high-energy electron acceleration becomes less efficient when the mass ratio gets larger because the \textit{Fermi}-like mechanism associated with particle curvature drift becomes less efficient. These results indicate that when particle curvature drift dominates high-energy particle acceleration, the further the particle kinetic scales are from the magnetic field curvature scales ($\sim d_i$), the weaker the acceleration will be.

Magnetic reconnection is a fundamental physical process in many laboratory, space, and astrophysical plasma contexts. Solar flares serve as an outstanding laboratory to study the magnetic reconnection and the associated energy release and conversion processes under plasma conditions difficult to reproduce in the laboratory, and with considerable spatiotemporal details not possible elsewhere in astrophysics. Here we emphasize the unique power of remote-sensing observations of solar flares at radio wavelengths. In particular, we discuss the transformative technique of broadband radio dynamic imaging spectroscopy in making significant contributions to addressing several outstanding challenges in magnetic reconnection, including the capability of pinpointing magnetic reconnection sites, measuring the time-evolving reconnecting magnetic fields, and deriving the spatially and temporally resolved distribution function of flare-accelerated electrons.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

In eruptive solar flares, termination shocks (TSs), formed when high-speed reconnection outflows collide with closed dense flaring loops, are believed to be one of the possible candidates for plasma heating and particle acceleration. In this work, we perform resistive magnetohydrodynamic simulations in a classic Kopp-Pneuman flare configuration to study the formation and evolution of TSs, and analyze in detail the dynamic features of TSs and variations of the shock strength in space and time. This research focuses on the fast reconnection phase when plasmoids form and produce small-scale structures inside the flare current sheet. It is found that the TS emerges once the downward outflow colliding with closed magnetic loops becomes super-magnetosonic, and immediately becomes highly dynamical. The morphology of a TS can be flat, oblique, or curved depending on the detailed interactions between the outflows/plasmoids and the highly dynamic plasma in the looptop region. The TS becomes weaker when a plasmoid is crossing through, or may even be destroyed by well developed plasmoids and then re-constructed above the plasmoids. We also perform detailed statistical analysis on important physical quantities along and across the shock front. The density and temperature ratios range from 1 to 3 across the TS front, and the pressure ratio typically has larger values up to 10. We show that weak guide fields do not strongly affect the Mach number and compression ratios, and the TS length becomes slightly larger in the case with thermal conduction.

In solar flares and other astrophysical systems, a major challenge for solving particle acceleration problem associated with magnetic reconnection is the enormous scale separation between kinetic scales and observed reconnection scale. Because of this, it has been difficult to draw any definite conclusions by just using kinetic simulations. Particle acceleration model that solves energetic particle transport equation can capture the main acceleration physics found in kinetic simulations, and thus provide a practical way to make observable predictions and directly compare model results with observations. Here we study compression particle acceleration in magnetic reconnection by solving Parker (diffusion-advection) transport equation using velocity and magnetic fields from two-dimensional high-Lundquist-number magnetohydrodynamics (MHD) simulations of a low-$\beta$ reconnection layer. We show that the compressible reconnection layer can give significant particle acceleration, leading to the formation of power-law particle energy distributions. We analyze the acceleration rate and find that the acceleration in the reconnection layer is a mixture of first order and second order Fermi processes. When including a guide field, we find the spectrum becomes steeper and both the power-law cutoff energy and maximum particle energy decrease as plasma becomes less compressible. This model produces 2D particle distribution that one can use to generate radiation map and directly compare with solar flare observations. This provides a framework to explain particle acceleration at large-scale astrophysical reconnection sites, such as solar flares.

Magnetic flux ropes are commonly observed throughout the heliosphere, and recent studies suggest that interacting flux ropes are associated with some energetic particle events. In this work, we carry out 2D particle-in-cell (PIC) simulations to study the coalescence of two magnetic flux ropes (or magnetic islands), and the subsequent plasma energization processes. The simulations are initialized with two magnetic islands embedded in a reconnecting current sheet. The two islands collide and eventually merge into a single island. Particles are accelerated during this process as the magnetic energy is released and converted to the plasma energy, including bulk kinetic energy increase by the ideal electric field, and thermal energy increase by the fluid compression and the non-ideal electric field. We find that contributions from these different energization mechanisms are all important and comparable with each other. Fluid shear and a non-gyrotropic pressure tensor also contribute to the energy conversion process. For simulations with different box sizes ranging from $L_x \sim 25$--$100 d_i$ and ion-to-electron mass ratios $m_i / m_e = 25$, 100 and 400, we find that the general evolution is qualitatively the same for all runs, and the energization depends only weakly on either the system size or the mass ratio. The results may help us understand plasma energization in solar and heliospheric environments.

Magnetic reconnection, the central engine that powers explosive phenomena throughout the Universe, is also perceived as one of the principal mechanisms for accelerating particles to high energies. Although various signatures of magnetic reconnection have been frequently reported, observational evidence that links particle acceleration directly to the reconnection site has been rare, especially for space plasma environments currently inaccessible to $\textit{in situ}$ measurements. Here we utilize broadband radio dynamic imaging spectroscopy available from the Karl G. Jansky Very Large Array to observe decimetric type III radio bursts in a solar jet with high angular ($\sim$20$''$), spectral ($\sim$1 %), and temporal resolution (50 milliseconds). These observations allow us to derive detailed trajectories of semi-relativistic (tens of keV) electron beams in the low solar corona with unprecedentedly high angular precision ($<0''.65$). We found that each group of electron beams, which corresponds to a cluster of type III bursts with 1-2-second duration, diverges from an extremely compact region ($\sim$600 km$^2$) in the low solar corona. The beam-diverging sites are located behind the erupting jet spire and above the closed arcades, coinciding with the presumed location of magnetic reconnection in the jet eruption picture supported by extreme ultraviolet/X-ray data and magnetic modeling. We interpret each beam-diverging site as a reconnection null point where multitudes of magnetic flux tubes join and reconnect. Our data suggest that the null points likely consist of a high level of density inhomogeneities possibly down to 10-km scales. These results, at least in the present case, strongly favor a reconnection-driven electron acceleration scenario.

Relativistic magnetic reconnection events can widely exist in magnetized plasmas in astrophysical systems. During this process, oppositely directed magnetic field lines reconnect and release magnetic energy, efficiently accelerating nonthermal particles. However, so far there is little clear observational signatures of relativistic magnetic reconnection events in astrophysical systems. Blazars are relativistic magnetized plasma outflows from supermassive black holes. Their multi-wavelength flares may be powered by relativistic magnetic reconnection. The highly variable radiation and polarization signatures are well covered by multi-wavelength observation campaigns, making them ideal targets to examine the magnetic reconnection model. Recent observations have found that several blazar flares are accompanied by optical polarization angle swings which may be of as large amplitude as $> 180^{\circ}$, challenging existing theoretical models. In this paper, we present integrated particle-in-cell (PIC) and polarized radiation transfer simulations of magnetic reconnection events. We find that plasmoid coalescences in the reconnection layer can give rise to highly variable light curves, low and fluctuating polarization degree, and rotating polarization angle. In particular, large-amplitude polarization angle swings, similar to those observed during blazar flares, can be a unique signature of relativistic magnetic reconnection events.

This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for a given accretion rate are adopted and their interaction with the gas in the host galaxy is calculated. Physical processes such as star formation and SNe feedback are taken into account. Consistent with observation, we find the AGN spends most of the time in the low-luminosity regime. AGN feedback overall suppresses the star formation; but depending on location in the galaxy and time, it can also enhance it. The light curve of specific star formation rate is not synchronous with the AGN light curve. We find that wind usually plays a dominant role in controlling the AGN luminosity and star formation, but radiation also cannot be neglected.

Galaxy clusters host a large reservoir of diffuse plasma with radially-varying temperature profiles. The efficiency of thermal conduction in the intracluster medium (ICM) is complicated by the existence of turbulence and magnetic fields, and has received a lot of attention in the literature. Previous studies suggest that the magnetothermal instability developed in outer regions of galaxy clusters would drive magnetic field lines preferentially radial, resulting in efficient conduction along the radial direction. Using a series of spherically-symmetric simulations, here we investigate the impact of thermal conduction on the observed temperature distributions in outer regions of three massive clusters, and find that thermal conduction substantially modifies the ICM temperature profile. Within 3 Gyr, the gas temperature at a representative radius of $0.3r_{500}$ typically decreases by ~10 - 20% and the average temperature slope between $0.3r_{500}$ and $r_{500}$ drops by ~ 30 - 40%, indicating that the observed ICM would not stay in a long-term equilibrium state in the presence of thermal conduction. However, X-ray observations show that the outer regions of massive clusters have remarkably similar radially-declining temperature profiles, suggesting that they should be quite stable. Our study thus suggests that the effective conductivity along the radial direction must be suppressed below the Spitzer value by a factor of 10 or more, unless additional heating sources offset conductive cooling and maintain the observed temperature distributions. Our study provides a smoking-gun evidence for the suppression of parallel conduction along magnetic field lines in low-collisionality plasmas by kinetic mirror or whistler instabilities.

Recent Chandra X-ray observations of many galaxy clusters find evidence for hot metal-rich outflows preferentially aligned with the large-scale axes of X-ray cavities with typical outflow masses of around $10^{9} $ - $10^{10} M_{\odot}$. Here we perform a suite of three hydrodynamic simulations to investigate whether AGN jets could drive these metal-rich outflows in a representative cluster. By using both the tracer variable and virtual particle methods, and additionally following the gas metallicity evolution, we show that metal-rich gas initially located in central regions can indeed be uplifted by the AGN bubble to large distances, a phenomenon called Darwin drift in fluid mechanics, and forming a filamentary trailing outflow extending beyond $100$ kpc behind the bubble. The gas entrained in the trailing outflow is entirely outflowing with an average outflow rate of nearly $100M_{\odot}$/yr during the first $100$ Myr, and at later times, a growing lower part flows back towards the cluster center due to gravity. The outflow mass rises up to about $10^{10} M_{\odot}$ with entrained iron mass of about $10^{6} - 10^{7}M_{\odot}$, consistent with observations and predictions from the drift model. By the end of our simulation ($\sim 800$ Myr after the AGN event), several $10^{9}M_{\odot}$ of the uplifted high-metallicity gas still remains at large altitudes, potentially contributing to the enrichment of the bulk ICM and the broadening of central metallicity peaks observed in cool core clusters.

We demonstrate the dragging of the magnetic field by the super-Alfvenic shear flows out of the reconnection plane can strongly localize the reconnection x-line in collisionless plasmas, reversing the current direction at the x-line. Reconnection with this new morphology, which is impossible in resistive-magnetohydrodynamic (MHD), is enabled by electron inertia. Surprisingly, the quasi-steady reconnection rate remains of order 0.1 even though the aspect ratio of the local x-line geometry is larger than unity. We explain this by examining the transport of the reconnected magnetic flux and the opening angle made by the upstream magnetic field, concluding that the reconnection rate is still limited by the constraint imposed at the inflow region. This study further suggests the nearly universal fast rate value of order 0.1 cannot be explained by the physics of tearing modes, nor can it be explained by a universal localization mechanism.

We investigate the effects of AGN feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. The inner boundary of the simulation is chosen so that the Bondi radius is resolved. Compared to previous works, the two accretion modes, namely hot and cold, which correspond to different accretion rates and have different radiation and wind outputs, are carefully discriminated and the feedback effects by radiation and wind in each mode are taken into account. The most updated AGN physics, including the descriptions of radiation and wind from the hot accretion flows and wind from cold accretion disks, are adopted. Physical processes like star formation, Type Ia and Type II supernovae are taken into account. We study the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. We compare our simulation results with observations and find general consistency. Comparisons with previous simulation works find significant differences, indicating the importance of AGN physics. The respective roles of radiation and wind feedbacks are examined and it is found that they are different for different problems of interest such as AGN luminosity and star formation. We find that it is hard to neglect any of them, so we suggest to use the names of "cold feedback mode" and "hot feedback mode" to replace the currently used ones.

Particle acceleration in space and astrophysical reconnection sites is an important unsolved problem in studies of magnetic reconnection. Earlier kinetic simulations have identified several acceleration mechanisms that are associated with particle drift motions. Here, we show that, for sufficient large systems, the energization processes due to particle drift motions can be described as fluid compression and shear, and that the shear energization is proportional to the pressure anisotropy of energetic particles. By analyzing results from fully kinetic simulations, we show that the compression energization dominates the acceleration of high-energy particles in reconnection with a weak guide field, and the compression and shear effects are comparable when the guide field is 50\% of the reconnecting component. Spatial distributions of those energization effects reveal that reconnection exhausts, contracting islands, and island merging regions are the three most important regions for compression and shear acceleration. This study connects particle energization by particle guiding-center drift motions with that due to background fluid motions, as in the energetic particle transport theory. It provides foundations for building particle transport models for large-scale reconnection acceleration such as those in solar flares.

The unexpectedly high flux of cosmic ray positrons detected at Earth may originate from nearby astrophysical sources, dark matter, or unknown processes of cosmic-ray secondary production. We report the detection, using the HighAltitude Water Cherenkov Observatory (HAWC), of extended tera-electron volt gamma-ray emission coincident with the locations of two nearby middle-aged pulsars (Geminga and PSR B0656+14). The HAWC observations demonstrate that these pulsars are indeed local sources of accelerated leptons, but the measured tera-electron volt emission profile constrains the diffusion of particles away from these sources to be much slower than previously assumed. We demonstrate that the leptons emitted by these objects are therefore unlikely to be the origin of the excess positrons, which may have a more exotic origin.

Recent observations have shown that coronal shocks driven by coronal mass ejections can develop and accelerate particles within several solar radii in large solar energetic particle (SEP) events. Motivated by this, we present an SEP acceleration study including the process that a fast shock propagates through a streamer-like magnetic field with both closed and open field lines in the low corona region. The acceleration of protons is modeled by numerically solving the Parker transport equation with spatial diffusion both along and across the magnetic field. We show that particles can be sufficiently accelerated to up to several hundred MeV within 2-3 solar radii. When the shock propagates through a streamer-like magnetic field, particles are more efficiently accelerated compared to the case with a simple radial magnetic field, mainly due to perpendicular shock geometry and the natural trapping effect of closed magnetic fields. Our results suggest that the coronal magnetic field configuration is an important factor for producing large SEP events. We further show that the coronal magnetic field configuration strongly influences the distribution of energetic particles, leading to different locations of source regions along the shock front where most of high energy particles are concentrated. This work may have strong implications to SEP observations. The upcoming Parker Solar Probe will provide in situ observations for the distribution of energetic particles in the coronal shock region, and test the results of the study.

The cooling flow problem is one of the central problems in galaxy clusters, and active galactic nucleus (AGN) feedback is considered to play a key role in offsetting cooling. However, how AGN jets heat and suppress cooling flows remains highly debated. Using an idealized simulation of a cool-core cluster, we study the development of central cooling catastrophe and how a subsequent powerful AGN jet event averts cooling flows, with a focus on complex gasdynamical processes involved. We find that the jet drives a bow shock, which reverses cooling inflows and overheats inner cool core regions. The shocked gas moves outward in a rarefaction wave, which rarefies the dense core and adiabatically transports a significant fraction of heated energy to outer regions. As the rarefaction wave propagates away, inflows resume in the cluster core, but a trailing outflow is uplifted by the AGN bubble, preventing gas accumulation and catastrophic cooling in central regions. Inflows and trailing outflows constitute meridional circulations in the cluster core. At later times, trailing outflows fall back to the cluster centre, triggering central cooling catastrophe and potentially a new generation of AGN feedback. We thus envisage a picture of cool cluster cores going through cycles of cooling-induced contraction and AGN-induced expansion. This picture naturally predicts an anti-correlation between the gas fraction (or X-ray luminosity) of cool cores and the central gas entropy, which may be tested by X-ray observations.

In this paper, we report the peculiar HI morphology of the cluster spiral galaxy NGC 6145, which has a 150 kpc HI filament on one side that is nearly parallel to its major axis. This filament is made up of several HI clouds and the diffuse HI gas between them, with no optical counterparts. We compare its HI distribution with other one-sided HI distributions in the literature, and find that the overall HI distribution is very different from the typical tidal and ram-pressure stripped HI shape, and its morphology is inconsistent with being a pure accretion event. Only about 30% of the total HI gas is anchored on the stellar disk, while most of HI gas forms the filament in the west. At a projected distance of 122 kpc, we find a massive elliptical companion (NGC 6146) with extended radio emission, whose axis points to an HI gap in NGC 6145. The velocity of the HI filament shows an overall light-of- sight motion of 80 to 180 km/s with respect to NGC 6145. Using the long-slit spectra of NGC 6145 along its major stellar axis, we find that some outer regions show enhanced star formation, while in contrast, almost no star formation activities are found in its center (less than 2 kpc). Pure accretion, tidal or ram-pressure stripping is not likely to produce the observed HI filament. An alternative explanation is the jet-stripping from NGC 6146, although direct evidence for a jet-cold gas interaction has not been found.

Kink instabilities are likely to occur in the current-carrying magnetized plasma jets. Recent observations of the blazar radiation and polarization signatures suggest that the blazar emission region may be considerably magnetized. While the kink instability has been studied with first-principle magnetohydrodynamic (MHD) simulations, the corresponding time-dependent radiation and polarization signatures have not been investigated. In this paper, we perform comprehensive polarization-dependent radiation modeling of the kink instability in the blazar emission region based on relativistic MHD (RMHD) simulations. We find that the kink instability may give rise to strong flares with polarization angle (PA) swings or weak flares with polarization fluctuations, depending on the initial magnetic topology and magnetization. These findings are consistent with observations. Compared with the shock model, the kink model generates polarization signatures that are in better agreement with the general polarization observations. Therefore, we suggest that kink instabilities may widely exist in the jet environment, and provide an efficient way to convert the magnetic energy and produce multiwavelength flares and polarization variations.

Simulations suggest collisionless steady-state magnetic reconnection of Harris-type current sheets proceeds with a rate of order 0.1, independent of dissipation mechanism. We argue this long-standing puzzle is a result of constraints at the magnetohydrodynamic (MHD) scale. We perform a scaling analysis of the reconnection rate as a function of the opening angle made by the upstream magnetic fields, finding a maximum reconnection rate close to 0.2. The predictions compare favorably to particle-in-cell simulations of relativistic electron-positron and non-relativistic electron-proton reconnection. The fact that simulated reconnection rates are close to the predicted maximum suggests reconnection proceeds near the most efficient state allowed at the MHD-scale. The rate near the maximum is relatively insensitive to the opening angle, potentially explaining why reconnection has a similar fast rate in differing models.

The nature and origin of the Fermi bubbles detected in the inner Galaxy remain elusive. In this paper, we briefly discuss some recent theoretical and observational developments, with a focus on the AGN jet model. Analogous to radio lobes observed in massive galaxies, the Fermi bubbles could be naturally produced by a pair of opposing jets emanating nearly along the Galaxy's rotation axis from the Galactic center. Our two-fluid hydrodynamic simulations reproduce quite well the bubble location and shape, and interface instabilities at the bubble surface could be effectively suppressed by shear viscosity. We briefly comment on some potential issues related to our model, which may lead to future progress.

Type III and type-III-like radio bursts are produced by energetic electron beams guided along coronal magnetic fields. As a variant of type III bursts, Type N bursts appear as the letter "N" in the radio dynamic spectrum and reveal a magnetic mirror effect in coronal loops. Here, we report a well-observed N-shaped burst consisting of three successive branches at metric wavelength with both fundamental and harmonic components and a high brightness temperature ($>$10$^9$ K). We verify the burst as a true type N burst generated by the same electron beam from three aspects of the data. First, durations of the three branches at a given frequency increase gradually, may due to the dispersion of the beam along its path. Second, the flare site, as the only possible source of non-thermal electrons, is near the western feet of large-scale closed loops. Third, the first branch and the following two branches are localized at different legs of the loops with opposite sense of polarization. We also find that the sense of polarization of the radio burst is in contradiction to the O-mode and there exists a fairly large time delay ($\sim$3-5 s) between the fundamental and harmonic components. Possible explanations accounting for these observations are presented. Assuming the classical plasma emission mechanism, we can infer coronal parameters such as electron density and magnetic field near the radio source and make diagnostics on the magnetic mirror process.

A model of global magnetic reconnection rate in relativistic collisionless plasmas is developed and validated by the fully kinetic simulation. Through considering the force balance at the upstream and downstream of the diffusion region, we show that the global rate is bounded by a value $\sim 0.3$ even when the local rate goes up to $\sim O(1)$ and the local inflow speed approaches the speed of light in strongly magnetized plasmas. The derived model is general and can be applied to magnetic reconnection under widely different circumstances.

Radio-mode AGN feedback plays a key role in the evolution of galaxy groups and clusters. Its physical origin lies in the kpc-scale interaction of AGN jets with the intracluster medium. Large-scale jet simulations often initiate light internally-supersonic jets with density contrast $0.01<\eta<1$. Here we argue for the first time for the importance of very-light ($\eta<0.01$) internally-subsonic jets. We investigated the shapes of young X-ray cavities produced in a suite of hydrodynamic simulations, and found that bottom-wide cavities are always produced by internally-subsonic jets, while internally-supersonic jets inflate cylindrical, center-wide, or top-wide cavities. We found examples of real cavities with shapes analogous to those inflated in our simulations by internally-subsonic and internally-supersonic jets, suggesting a dichotomy of AGN jets according to their internal Mach numbers. We further studied the long-term cavity evolution, and found that old cavities resulted from light jets spread along the jet direction, while those produced by very light jets are significantly elongated along the perpendicular direction. The northwestern ghost cavity in Perseus is pancake-shaped, providing tentative evidence for the existence of very light jets. Our simulations show that very-light internally-subsonic jets decelerate faster and rise much slower in the ICM than light internally-supersonic jets, possibly depositing a larger fraction of jet energy to cluster cores and alleviating the problem of low coupling efficiencies found previously. The internal Mach number points to the jet's energy content, and internally-subsonic jets are energetically dominated by non-kinetic energy, such as thermal energy, cosmic rays, or magnetic fields.

Plasma energization through magnetic reconnection in the magnetically-dominated regime featured by low plasma beta ($\beta = 8 \pi nkT_0/B^2 \ll 1$) and/or high magnetization ($\sigma = B^2/(4 \pi nmc^2) \gg 1$) is important in a series of astrophysical systems such as solar flares, pulsar wind nebula, and relativistic jets from black holes, etc. In this paper, we review the recent progress on kinetic simulations of this process and further discuss plasma dynamics and particle acceleration in a low-$\beta$ reconnection layer that consists of electron-positron pairs. We also examine the effect of different initial thermal temperatures on the resulting particle energy spectra. While earlier papers have concluded that the spectral index is smaller for higher $\sigma$, our simulations show that the spectral index approaches $p=1$ for sufficiently low plasma $\beta$, even if $\sigma \sim 1$. Since this predicted spectral index in the idealized limit is harder than most observations, it is important to consider effects that can lead to a softer spectrum such as open boundary simulations. We also remark that the effects of 3D reconnection physics and turbulence on reconnection need to be addressed in the future.

With a test-particle simulation, we investigate the effect of large-scale coronal magnetic fields on electron acceleration at an outward-propagating coronal shock with a circular front. The coronal field is approximated by an analytical solution with a streamer-like magnetic field featured by partially open magnetic field and a current sheet at the equator atop the closed region. We show that the large-scale shock-field configuration, especially the relative curvature of the shock and the magnetic field line across which the shock is sweeping, plays an important role in the efficiency of electron acceleration. At low shock altitudes, when the shock curvature is larger than that of magnetic field lines, the electrons are mainly accelerated at the shock flanks; at higher altitudes, when the shock curvature is smaller, the electrons are mainly accelerated at the shock nose around the top of closed field lines. The above process reveals the shift of efficient electron acceleration region along the shock front during its propagation. It is also found that in general the electron acceleration at the shock flank is not so efficient as that at the top of closed field since at the top a collapsing magnetic trap can be formed. In addition, we find that the energy spectra of electrons is power-law like, first hardening then softening with the spectral index varying in a range of -3 to -6. Physical interpretations of the results and implications on the study of solar radio bursts are discussed.

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization $\sigma_0=B^2/(4\pi n(m_i+m_e)c^2)$, the magnetization for each species is $\sigma_i \sim \sigma_0$ and $\sigma_e \sim (m_i/m_e) \sigma_0$, respectively. We have studied the magnetically dominated regime by varying $\sigma_{e} = 10^3 - 10^5$ with initial ion and electron temperatures $T_i = T_e = 5 - 20 m_ec^2$ and mass ratio $m_i/m_e = 1 - 1836$. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several ($2-3$) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1<s<2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be $\gamma_{be} \sim \sigma_e$ and $\gamma_{bi} \sim \sigma_i$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices $s_p$ are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

By means of fully kinetic simulations, we investigate electron acceleration during magnetic reconnection in a nonrelativistic proton--electron plasma with conditions similar to solar corona and flares. We demonstrate that reconnection leads to a nonthermally dominated electron acceleration with a power-law energy distribution in the nonrelativistic low-$\beta$ regime but not in the high-$\beta$ regime, where $\beta$ is the ratio of the plasma thermal pressure and the magnetic pressure. The accelerated electrons contain most of the dissipated magnetic energy in the low-$\beta$ regime. A guiding-center current description is used to reveal the role of electron drift motions during the bulk nonthermal energization. We find that the main acceleration mechanism is a \textit{Fermi}-type acceleration accomplished by the particle curvature drift motion along the electric field induced by the reconnection outflows. Although the acceleration mechanism is similar for different plasma $\beta$, low-$\beta$ reconnection drives fast acceleration on Alfv\'enic timescales and develops power laws out of thermal distribution. The nonthermally dominated acceleration resulting from magnetic reconnection in low-$\beta$ plasma may have strong implications for the highly efficient electron acceleration in solar flares and other astrophysical systems.

Magnetic reconnection is thought to be the driver for many explosive phenomena in the universe. The energy release and particle acceleration during reconnection have been proposed as a mechanism for producing high-energy emissions and cosmic rays. We carry out two- and three-dimensional kinetic simulations to investigate relativistic magnetic reconnection and the associated particle acceleration. The simulations focus on electron-positron plasmas starting with a magnetically dominated, force-free current sheet ($\sigma \equiv B^2/(4\pi n_e m_e c^2) \gg 1$). For this limit, we demonstrate that relativistic reconnection is highly efficient at accelerating particles through a first-order Fermi process accomplished by the curvature drift of particles along the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra $f \propto (\gamma-1)^{-p}$ and approaches $p = 1$ for sufficiently large $\sigma$ and system size. Eventually most of the available magnetic free energy is converted into nonthermal particle kinetic energy. An analytic model is presented to explain the key results and predict a general condition for the formation of power-law distributions. The development of reconnection in these regimes leads to relativistic inflow and outflow speeds and enhanced reconnection rates relative to non-relativistic regimes. In the three-dimensional simulation, the interplay between secondary kink and tearing instabilities leads to strong magnetic turbulence, but does not significantly change the energy conversion, reconnection rate, or particle acceleration. This study suggests that relativistic reconnection sites are strong sources of nonthermal particles, which may have important implications to a variety of high-energy astrophysical problems.

X-ray observations of galaxy clusters have detected numerous X-ray cavities, evolved from the interaction of AGN jets with the intracluster medium (ICM) and providing compelling evidence for the importance of jet-mode AGN feedback. Here we argue for the physical importance of the cavity shape, which we characterize with two geometric parameters: radial elongation $\tau$ and top wideness $b$. We study the cavity shape with 16 hydrodynamic jet simulations in two representative clusters, and find that the shapes of young cavities are mainly determined by various jet properties. Our simulations successfully reproduce two observed types of young cavities elongated along either the jet ($\tau>1$; type-II) or perpendicular ($\tau\leq1$; type-I) direction. Bottom-wide type-I cavities are produced by very light internally-subsonic jets, while top-wide type-II cavities are produced by heavier, internally-supersonic jets, which may also produce center-wide cavities with $\tau\sim 1$ if the jets are only slightly supersonic. Bottom-wide type-II cavities can be produced by very light jets with very long durations and cylindrical cavities are produced by very light internally-supersonic jets. While not appreciably affecting the shapes of young cavities, viscosity significantly affects the long-term cavity evolution, suppressing both interface instabilities and the formation of torus-like morphology. We encourage observers to study the shapes of young and old X-ray cavities separately, the former probing the properties of AGN jets and the latter potentially probing the ICM viscosity level.

The polarization signatures of the blazar emissions are known to be highly variable. In addition to small fluctuations of the polarization angle around a mean value, sometimes large (> 180^o) polarization angle swings are observed. We suggest that such p henomena can be interpreted as arising from light-travel-time effects within an underlying axisymmetric emission region. We present the first simultaneous fitting of the multi-wavelength spectrum, variability and time-dependent polarization features of a correlated optical and gamma-ray flaring event of the prominent blazar 3C279, which was accompanied by a drastic change of its polarization signatures. This unprecedented combination of spectral, variability, and polarization information in a coherent physical model allows us to place stringent constraints on the particle acceleration and magnetic-field topology in the relativistic jet of a blazar, strongly favoring a scenario in which magnetic energy dissipation is the primary driver of the flare event.

We perform a numerical-simulation study of the acceleration of electrons at shocks that propagate through a prespecified, kinematically defined turbulent magnetic field. The turbulence consists of broadband magnetic fluctuations that are embedded in the plasma and cover a range of wavelengths, the smallest of which is larger than the gyroadii of electrons that are initially injected into the system. We find that when the variance of the turbulent component of the upstream magnetic field is sufficiently large -- $\sigma^2 \sim$ 10 $B_0^2$, where $B_0$ is the strength of the background magnetic field -- electrons can be efficiently accelerated at a collisionless shock regardless of the orientation of the mean upstream magnetic field relative to the shock-normal direction. Since the local angle between the incident magnetic-field vector and the shock-normal vector can be quite large, electrons can be accelerated through shock-drift acceleration at the shock front. In the upstream region, electrons are mirrored back to the shock front leading to multiple shock encounters. Eventually the accelerated electrons are energetic enough that their gyroradii are of the same order as the wavelength of waves that are included in our description of the turbulent magnetic field. Our results are consistent with recent \textit{in situ} observations at Saturn's bow shock. The study may help understand the acceleration of electrons at shocks in space and astrophysical systems.

Recent observations suggested that star formation quenching in galaxies is related to galaxy structure. Here we propose a new mechanism to explain the physical origin of this correlation. We assume that while quiescent galaxies are maintained quenched by a feedback mechanism, cooling flows in the hot halo gas can still develop intermittently. We study cooling flows in a large suite of around 90 hydrodynamic simulations of an isolated galaxy group, and find that the flow development depends significantly on the gravitational potential well in the central galaxy. If the galaxy's gravity is not strong enough, cooling flows result in a central cooling catastrophe, supplying cold gas and feeding star formation to galactic bulges. When the bulge grows prominent enough, compressional heating starts to offset radiative cooling and maintains cooling flows in a long-term hot mode without producing cooling catastrophe. Our model thus describes a self-limited growth channel for galaxy bulges, and naturally explains the connection between quenching and bulge prominence. In particular, we explicitly demonstrate that $M_{*}/R_{\rm eff}^{1.5}$ is a good structural predictor of quenching. We further find that the gravity from the central supermassive black hole also affects the bimodal fate of cooling flows, and predict a more general quenching predictor to be $M_{\rm bh}^{1.6}M_{*}/R_{\rm eff}^{1.5}$, which may be tested in future observational studies.

Magnetic reconnection is believed to be the dominant energy release mechanism in solar flares. The standard flare model predicts both downward and upward outflow plasmas with speeds close to the coronal Alfv\'{e}n speed. Yet, spectroscopic observations of such outflows, especially the downflows, are extremely rare. With observations of the newly launched Interface Region Imaging Spectrograph (IRIS), we report the detection of greatly redshifted ($\sim$125 km s$^{-1}$ along line of sight) Fe {\sc{xxi}} 1354.08\AA{} emission line with a $\sim$100 km s$^{-1}$ nonthermal width at the reconnection site of a flare. The redshifted Fe {\sc{xxi}} feature coincides spatially with the loop-top X-Ray source observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). We interpret this large redshift as the signature of downward-moving reconnection outflow/hot retracting loops. Imaging observations from both IRIS and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) also reveal the eruption and reconnection processes. Fast downward-propagating blobs along these loops are also found from cool emission lines (e.g., Si {\sc{iv}}, O {\sc{iv}}, C {\sc{ii}}, Mg {\sc{ii}}) and images of AIA and IRIS. Furthermore, the entire Fe {\sc{xxi}} line is blueshifted by $\sim$260 km s$^{-1}$ at the loop footpoints, where the cool lines mentioned above all exhibit obvious redshift, a result that is consistent with the scenario of chromospheric evaporation induced by downward-propagating nonthermal electrons from the reconnection site.

Two solar type II radio bursts, separated by ~24 hours in time, are examined together. Both events are associated with coronal mass ejections (CMEs) erupting from the same active region (NOAA 11176) beneath a well-observed helmet streamer. We find that the type II emissions in both events ended once the CME/shock fronts passed the white-light streamer tip, which is presumably the magnetic cusp of the streamer. This leads us to conjecture that the closed magnetic arcades of the streamer may play a role in electron acceleration and type II excitation at coronal shocks. To examine such a conjecture, we conduct a test-particle simulation for electron dynamics within a large-scale partially-closed streamer magnetic configuration swept by a coronal shock. We find that the closed field lines play the role of an electron trap, via which the electrons are sent back to the shock front for multiple times, and therefore accelerated to high energies by the shock. Electrons with an initial energy of 300eV can be accelerated to tens of keV concentrating at the loop apex close to the shock front with a counter-streaming distribution at most locations. These electrons are energetic enough to excite Langmuir waves and radio bursts. Considering the fact that most solar eruptions originate from closed field regions, we suggest that the scenario may be important to the generation of more metric type IIs. This study also provides an explanation to the general ending frequencies of metric type IIs at or above 20-30 MHz and the disconnection issue between metric and interplanetary type IIs.

Using fully kinetic simulations, we demonstrate that magnetic reconnection in relativistic plasmas is highly efficient at accelerating particles through a first-order Fermi process resulting from the curvature drift of particles in the direction of the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra in parameter regimes where the energy density in the reconnecting field exceeds the rest mass energy density $\sigma \equiv B^2/(4 \pi n m_ec^2) > 1$ and when the system size is sufficiently large. In the limit $\sigma \gg 1$, the spectral index approaches $p=1$ and most of the available energy is converted into non-thermal particles. A simple analytic model is proposed which explains these key features and predicts a general condition under which hard power-law spectra will be generated from magnetic reconnection.

Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter $\sigma$ and approaches the light speed when $\sigma > O(100)$, leading to an enhanced reconnection rate. In all regimes, the divergence of pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region remains $\sim 0.1$ in both the non-relativistic and relativistic limits.

Simultaneous radio and extreme ultraviolet (EUV)/white-light imaging data are examined for a solar type II radio burst occurring on 2010 March 18 to deduce its source location. Using a bow-shock model, we reconstruct the 3-dimensional EUV wave front (presumably the type-II emitting shock) based on the imaging data of the two STEREO spacecraft. It is then combined with the Nan\c{c}ay radio imaging data to infer the 3-dimensional position of the type II source. It is found that the type II source coincides with the interface between the CME EUV wave front and a nearby coronal ray structure, providing evidence that the type II emission is physically related to the CME-ray interaction. This result, consistent with those of previous studies, is based on simultaneous radio and EUV imaging data for the first time.

We report observations of the acceleration and trapping of energetic ions and electrons between a pair of corotating interaction regions (CIRs). The event occurred in Carrington Rotation 2060. Observed at spacecraft STEREO-B, the two CIRs were separated by less than 5 days. In contrast to other CIR events, the fluxes of energetic ions and electrons in this event reached their maxima between the trailing-edge of the first CIR and the leading edge of the second CIR. The radial magnetic field (Br) reversed its sense and the anisotropy of the flux also changed from sunward to anti-sunward between the two CIRs. Furthermore, there was an extended period of counter-streaming suprathermal electrons between the two CIRs. Similar observations for this event were also obtained for ACE and STEREO-A. We conjecture that these observations were due to a "U-shape" large scale magnetic field topology connecting the reverse shock of the first CIR and the forward shock of the second CIR. Such a disconnected U-shaped magnetic field topology may have formed due to magnetic reconnection in the upper corona.

Feedback heating from active galactic nuclei (AGNs) has been commonly invoked to suppress cooling flows predicted in hot gas in elliptical galaxies, galaxy groups and clusters. Previous studies have focused on if and how AGN feedback heats the gas, but little paid attention to its triggering mechanism. Using spherically symmetric simulations, we investigate how large-scale cooling flows are accreted by central supermassive black holes (SMBHs) in eight well-observed systems and find an interesting dichotomy. In massive clusters, the gas develops a central cooling catastrophe within about the cooling time (typically ~ 100 - 300 Myr), resulting in a cold-mode accretion onto SMBHs. However, in our four simulated systems on group and galaxy scales at a low metallicity Z=0.3Z_{sun}, the gas quickly settles into a long-term state which has a cuspy central temperature profile extending to several tens to about 100 pc. At the more realistic solar metallicity, two groups (with R_{e} ~ 4 kpc) still host the long-term hot-mode accretion. Both accretion modes naturally appear in our idealized calculations where only cooling, gas inflow, and compressional heating are considered. The long-term hot-mode accretion is maintained by the quickly-established closeness between the timescales of these processes, preferably in systems with low gas densities, low gas metallicities, and importantly, compact central galaxies, which result in strong gravitational acceleration and compressional heating at the intermediate radii. Our calculations predict that central cuspy temperature profiles appear more often in smaller systems than galaxy clusters, which instead often host significant cold gas and star formation.

Using numerical simulations of charged-particles propagating in the heliospheric magnetic field, we study small-scale gradients, or "dropouts", in the intensity of solar energetic particles seen at 1 AU. We use two turbulence models, the foot-point random motion model (Jokipii & Parker 1969; Giacalone et al. 2006) and two-component model (Matthaeus et al. 1990), to generate fluctuating magnetic fields similar to spacecraft observations at 1 AU. The turbulence models include a Kolmogorov-like magnetic field power spectrum containing a broad range of spatial scales from those that lead to large-scale field-line random walk to small scales leading to resonant pitch-angle scattering of energetic particles. We release energetic protons (20 keV - 10 MeV) from a spatially compact and instantaneous source. The trajectories of energetic charged particles in turbulent magnetic fields are numerically integrated. Spacecraft observations are mimicked by collecting particles in small windows when they pass the windows at a distance of 1 AU. We show that small-scale gradients in the intensity of energetic particles and velocity dispersions observed by spacecraft can be reproduced using the foot-point random motion model. However, no dropouts is seen in simulations using the two-component magnetic turbulence model. We also show that particle scattering in the solar wind magnetic field needs to be infrequent for intensity dropouts to form.

We present three-dimensional hybrid simulations of collisionless shocks that propagate parallel to the background magnetic field to study the acceleration of protons that forms a high-energy tail on the distribution. We focus on the initial acceleration of thermal protons and compare it with results from one-dimensional simulations. We find that for both one- and three-dimensional simulations, particles that end up in the high-energy tail of the distribution later in the simulation gained their initial energy right at the shock. This confirms previous results but is the first to demonstrate this using fully three-dimensional fields. The result is not consistent with the "thermal leakage" model. We also show that the gyrocenters of protons in the three-dimensional simulation can drift away from the magnetic field lines on which they started due to the removal of ignorable coordinates that exist in one- and two-dimensional simulations. Our study clarifies the injection problem for diffusive shock acceleration.

The observed bimodality in radio luminosity in galaxy clusters is puzzling. We investigate the possibility that cosmic-ray (CR) streaming in the intra-cluster medium can 'switch off' hadronically induced radio and gamma-ray emission. For self-confined CRs, this depends on the source of MHD wave damping: if only non-linear Landau damping operates, then CRs stream on the slow Alfvenic timescale, but if turbulent wave damping operates, super-Alfvenic streaming is possible. As turbulence increases, it promotes outward streaming more than it enables inward turbulent advection. Curiously, the CR flux is independent of $\nabla f$ (as long as it is non-zero) and depends only on plasma parameters; this enables radio halos with flat inferred CR profiles to turn off. We perform 1D time-dependent calculations of a radio mini-halo (Perseus) and giant radio halo (Coma) and find that both diminish in radio luminosity by an order of magnitude in several hundred Myr, given plausible estimates for the magnetic field in the outskirts of the cluster. Due to the energy dependence of CR streaming, spectral curvature develops, and radio halos turn off more slowly at low frequencies -- properties consistent with observations. Similarly, CR streaming rapidly turns off gamma-ray emission at the high-energies probed by Cherenkov telescopes, but not at the low energies probed by Fermi. CR mediated wave-heating of the ICM is unaffected, as it is dominated by ~GeV CRs which stream Alfvenically.

In recent studies (Feng et al., 2012; Kong et al., 2012), we proposed that source properties of type II radio bursts can be inferred through a causal relationship between the special shape of the type II dynamic spectrum (e.g., bump or break) and simultaneous extreme ultraviolet (EUV)/white light imaging observations (e.g., CME-shock crossing streamer structures). As a further extension of these studies, in this paper we examine the CME event dated on December 31 2007 associated with a multiple type II radio burst. We identify the presence of two spectral bump features on the observed dynamic spectrum. By combining observational analyses of the radio spectral observations and the EUV-white light imaging data, we conclude that the two spectral bumps are resulted from CME-shock propagating across dense streamers on the southern and northern sides of the CME, respectively. It is inferred that the corresponding two type II emissions originate separately from the two CME-shock flanks where the shock geometries are likely quasi-perpendicular or oblique. Since the emission lanes are bumped as a whole within a relatively short time, it suggests that the type IIs with bumps of the study are emitted from spatially confined sources (with a projected lateral dimension smaller than 0.05-0.1 R$_\odot$ at a fundamental frequency level of 20-30 MHz).

After introduction we focus on: the transport of charged particles, the acceleration of ions at shocks, and the acceleration of electrons at shocks. Chapter 2 studies the propagation of solar energetic particles(SEPs) in turbulent magnetic fields. Particle trajectories in turbulent magnetic fields are numerically integrated. The turbulence includes a Kolmogorov-like power spectrum containing a broad range of scales. Small-scale variations in particle intensities(dropouts) and velocity dispersions can be reproduced. The result gives a constraint on the error of onset analysis for inferring SEP informations. We find that dropouts are rarely produced using the two-component model(Matthaeus et al., 1990). The result questions the turbulence model. Chapter 3 studies the acceleration of ions. We use 3-D hybrid simulations to study the acceleration of low-energy particles at parallel shocks. We find that particles gain energy by reflection at the shock. The protons can move off field lines in 3-D electric and magnetic fields. We also use a stochastic integration method to study diffusive shock acceleration including large-scale magnetic variations. The results can explain the observations of anomalous cosmic rays by Voyager 1. Chapter 4 studies electron acceleration at a shock in a turbulent magnetic field by combining hybrid simulations and test-particle simulations. The acceleration is enhanced by including large-scale turbulence. Since electrons mainly follow field lines, the field-line braiding allows electrons interacting with shock many times. Ripples also contribute to acceleration by mirroring electrons. The process favors perpendicular shocks. We discuss the implication to SEPs by comparing the acceleration of electrons with that of protons. The intensity correlation of electrons and ions in SEPs implies perpendicular shocks play important roles in accelerating particles.

The Fermi Gamma-ray Space Telescope reveals two large gamma-ray bubbles in the Galaxy, which extend about 50 degrees (~ 10 kpc) above and below the Galactic center (GC) and are symmetric about the Galactic plane. Using axisymmetric hydrodynamic simulations with a self-consistent treatment of the dynamical cosmic ray (CR) - gas interaction, we show that the bubbles can be created with a recent active galactic nucleus (AGN) jet activity about 1 - 3 Myr ago, which was active for a duration of ~ 0.1 - 0.5 Myr. The bipolar jets were ejected into the Galactic halo along the rotation axis of the Galaxy. Near the GC, the jets must be moderately light with a typical density contrast 0.001 <~ \eta <~ 0.1 relative to the ambient hot gas. The jets are energetically dominated by kinetic energy, and over-pressured with either CR or thermal pressure which induces lateral jet expansion, creating fat CR bubbles as observed. The sharp edges of the bubbles imply that CR diffusion across the bubble surface is strongly suppressed. The jet activity induces a strong shock, which heats and compresses the ambient gas in the Galactic halo, potentially explaining the ROSAT X-ray shell features surrounding the bubbles. The Fermi bubbles provide plausible evidence for a recent powerful AGN jet activity in our Galaxy, shedding new insights into the origin of the halo CR population and the channel through which massive black holes in disk galaxies release feedback energy during their growth.

The origin of the Fermi bubbles recently detected by the Fermi Gamma-ray Space Telescope in the inner Galaxy is mysterious. In the companion paper Guo & Mathews (Paper I), we use hydrodynamic simulations to show that they could be produced by a recent powerful AGN jet event. Here we further explore this scenario to study the potential roles of shear viscosity and cosmic ray (CR) diffusion on the morphology and CR distribution of the bubbles. We show that even a relatively low level of viscosity (\mu_{visc} >~ 3 g cm^{-1} s^{-1}, or ~0.1% - 1% of Braginskii viscosity in this context) could effectively suppress the development of Kelvin-Helmholtz instabilities at the bubble surface, resulting in smooth bubble edges as observed. Furthermore, viscosity reduces circulating motions within the bubbles, which would otherwise mix the CR-carrying jet backflow near bubble edges with the bubble interior. Thus viscosity naturally produces an edge-favored CR distribution, an important ingredient to produce the observed flat gamma-ray surface brightness distribution. Generically, such a CR distribution often produces a limb-brightened gamma-ray intensity distribution. However, we show that by incorporating CR diffusion which is strongly suppressed across the bubble surface (as inferred from sharp bubble edges) but is close to canonical values in the bubble interior, we obtain a reasonably flat gamma-ray intensity profile. The similarity of the resulting CR bubble with the observed Fermi bubbles strengthens our previous result in Paper I that the Fermi bubbles were produced by a recent AGN jet event. Studies of the nearby Fermi bubbles may provide a unique opportunity to study the potential roles of plasma viscosity and CR diffusion on the evolution of AGN jets and bubbles.

We present a study of the acceleration of electrons at a perpendicular shock that propagates through a turbulent magnetic field. The energization process of electrons is investigated by utilizing a combination of hybrid (kinetic ions and fluid electron) simulations and test-particle electron simulations. In this method, the motions of the test-particle electrons are numerically integrated in the time-dependent electric and magnetic fields generated by two-dimensional hybrid simulations. We show that large-scale magnetic fluctuations effect electrons in a number of ways and lead to efficient and rapid energization at the shock front. Since the electrons mainly follow along magnetic lines of force, the large-scale braiding of field lines in space allows the fast-moving electrons to interact with the shock front and get accelerated multiple times. Ripples in the shock front occurring at various scales will also contribute to the acceleration by mirroring the electrons. Our calculation shows that this process favors electron acceleration at perpendicular shocks. The acceleration efficiency is critically dependent on the turbulence amplitude and coherence length. We also discuss the implication of this study for solar energetic particles (SEPs) by comparing the acceleration of electrons with that of protons. Their correlation indicates that perpendicular shocks play an important role in SEP events.

We describe approximate axisymmetric computations of the dynamical evolution of material inside radio lobes and X-ray cluster gas cavities in Fanaroff-Riley II sources such as Cygnus A. All energy is delivered by a jet to the lobe/cavity via a moving hotspot where jet energy dissipates in a reverse shock. Our calculations describe the evolution of hot plasma, cosmic rays (CRs) and toroidal magnetic fields flowing from the hotspot into the cavity. Many observed features are explained. Gas, CRs and field flow back along the cavity surface in a "boundary backflow" consistent with detailed FRII observations. Computed ages of backflowing CRs are consistent with observed radio-synchrotron age variations only if shear instabilities in the boundary backflow are damped and we assume this is done with viscosity of unknown origin. Magnetic fields estimated from synchrotron self-Compton (SSC) X-radiation observed near the hotspot evolve into radio lobe fields. Computed profiles of radio synchrotron lobe emission perpendicular to the jet are dramatically limb-brightened in excellent agreement with FRII observations although computed lobe fields exceed those observed. Strong winds flowing from hotspots naturally create kpc-sized spatial offsets between hotspot inverse Compton (IC-CMB) X-ray emission and radio synchrotron emission that peaks 1-2 kpc ahead where the field increases due to wind compression. In our computed version of Cygnus A, nonthermal X-ray emission increases from the hotspot (some IC-CMB, mostly SSC) toward the offset radio synchrotron peak (mostly SSC). A faint thermal jet along the symmetry axis may be responsible for redirecting the Cygnus A non-thermal jet.

Steady accretion of a radiating gas onto a central mass point is described and compared to classic Bondi accretion. Radiation losses are essential for accretion flows to be observed. Unlike Bondi flows, radiating Bondi flows pass through a sonic point at a finite radius and become supersonic near the center. The morphology of all radiating Bondi flows is described by a single dimensionless parameter. In radiating Bondi flows the mass accretion rate varies approximately as the first power of the central mass -- this differs significantly from the quadratic dependence on the central mass in classical Bondi flows. Mass accretion rates onto galaxy or cluster-centered black holes estimated from traditional and radiating Bondi flows are significantly different. In radiating Bondi flows the gas temperature increases at large radii, as in the cores of many galaxy groups and clusters, allowing radiating Bondi flows to merge naturally with gas arriving from their cluster environments. Some radiating flows cool completely before reaching the center of the flow, and this also occurs in cooling site flows in which there is no central gravitating mass.

We investigate the acceleration of charged particles (both electrons and protons) at collisionless shocks predicted to exist in the vicinity of solar flares. The existence of standing termination shocks has been examined by flare models and numerical simulations e.g., Shibata,Forbes. We study electron energization by numerically integrating the equations of motion of a large number of test-particle electrons in the time-dependent two-dimensional electric and magnetic fields generated from hybrid simulations (kinetic ions and fluid electron) using parameters typical of the solar flare plasma environment. The shock is produced by injecting plasma flow toward a rigid piston. Large-scale magnetic fluctuations -- known to exist in plasmas and known to have important effects on the nonthermal electron acceleration at shocks -- are also included in our simulations. For the parameters characteristic of the flaring region, our calculations suggest that the termination shock formed in the reconnection outflow region (above post-flare loops) could accelerate electrons to a kinetic energy of a few MeV within 100 ion cyclotron periods, which is of the order of a millisecond. Given a sufficient turbulence amplitude level ($\delta B^2/B_0^2 \sim 0.3$), about 10% of thermal test-particle electrons are accelerated to more than 15 keV. We find that protons are also accelerated, but not to as high energy in the available time and the energy spectra are considerably steeper than that of the electrons for the parameters used in our simulations. Our results are qualitatively consistent with the observed hard X-ray emissions in solar flares.

We discuss an intriguing type II radio burst that occurred on 2011 March 27. The dynamic spectrum was featured by a sudden break at about 43 MHz on the well-observed harmonic branch. Before the break, the spectrum drifted gradually with a mean rate of about -0.05 MHz/s. Following the break, the spectrum jumped to lower frequencies. The post-break emission lasted for about three minutes. It consisted of an overall slow drift which appeared to have a few fast drift sub-bands. Simultaneous observations from the Solar TErrestrial RElations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) were also available and are examined for this event. We suggest that the slow-drift period before the break was generated inside a streamer by a coronal eruption driven shock, and the spectral break as well as the relatively wide spectrum after the break is a consequence of the shock crossing the streamer boundary where density drops abruptly. It is suggested that this type of radio bursts can be taken as a unique diagnostic tool for inferring the coronal density structure, as well as the radio emitting source region.

We have performed extensive two-dimensional magnetohydrodynamic simulations to study the amplification of magnetic fields when a supernova blast wave propagates into a turbulent interstellar plasma. The blast wave is driven by injecting high pressure in the simulation domain. The interstellar magnetic field can be amplified by two different processes, occurring in different regions. One is facilitated by the fluid vorticity generated by the ``rippled" shock front interacting with the background turbulence. The resulting turbulent flow keeps amplifying the magnetic field, consistent with earlier work \citep{Giacalone2007}. The other process is facilitated by the growth of the Rayleigh-Taylor instability at the contact discontinuity between the ejecta and the shocked medium. This can efficiently amplify the magnetic field and tends to produce the highest magnetic field. We investigate the dependence of the amplification on numerical parameters such as grid-cell size and on various physical parameters. We show the magnetic field has a characteristic radial profile that the downstream magnetic field gets progressively stronger away from the shock. This is because the downstream magnetic field needs a finite time to reach the efficient amplification, and will get further amplified in the Rayleigh-Taylor region. In our simulation we do not observe a systematic strong magnetic field within a small distance to the shock. This indicates that if the magnetic-field amplification in supernova remnants indeed occurs near the shock front, other processes such as three-dimensional instabilities, plasma kinetics and/or cosmic ray effect may need to be considered to explain the strong magnetic field in supernova remnants.

The total feedback energy injected into hot gas in galaxy clusters by central black holes can be estimated by comparing the potential energy of observed cluster gas profiles with the potential energy of non-radiating, feedback-free hot gas atmospheres resulting from gravitational collapse in clusters of the same total mass. Feedback energy from cluster-centered black holes expands the cluster gas, lowering the gas-to-dark matter mass ratio below the cosmic value. Feedback energy is unnecessarily delivered by radio-emitting jets to distant gas far beyond the cooling radius where the cooling time equals the cluster lifetime. For clusters of mass 4-11 X 10^14 Msun estimates of the total feedback energy, 1-3 X 10^63 ergs, far exceed feedback energies estimated from observations of X-ray cavities and shocks in the cluster gas, energies gained from supernovae, and energies lost from cluster gas by radiation. The time-averaged mean feedback luminosity is comparable to those of powerful quasars, implying that some significant fraction of this energy may arise from the spin of the black hole. The universal entropy profile in feedback-free gaseous atmospheres in NFW cluster halos can be recovered by multiplying the observed gas entropy profile of any relaxed cluster by a factor involving the gas fraction profile. While the feedback energy and associated mass outflow in the clusters we consider far exceed that necessary to stop cooling inflow, the time-averaged mass outflow at the cooling radius almost exactly balances the mass that cools within this radius, an essential condition to shut down cluster cooling flows.

The pair of large radio lobes in the Virgo cluster, each about 23 kpc in radius, have curiously sharp outer edges where the radio-synchrotron continuum flux declines abruptly. However, just adjacent to this sharp transition, the radio flux increases. This radio limb-brightening is observed over at least half of the perimeter of both lobes. We describe slowly propagating steady state diffusion fronts that explain these counterintuitive features. Because of the natural buoyancy of radio lobes, the magnetic field is largely tangent to the lobe boundary, an alignment that polarizes the radio emission and dramatically reduces the diffusion coefficient of relativistic electrons. As cosmic ray electrons diffuse slowly into the cluster gas, the local magnetic field and gas density are reduced as gas flows back toward the radio lobe. Radio emission peaks can occur because the synchrotron emissivity increases with magnetic field and then decreases with the density of non-thermal electrons. A detailed comparison of steady diffusion fronts with quantitative radio observations may reveal information about the spatial variation of magnetic fields and the diffusion coefficient of relativistic electrons. On larger scales, some reduction of the gas density inside the Virgo lobes due to cosmic ray pressure must occur and may be measurable. Such X-ray observations could reveal important information about the presence of otherwise unobservable non-thermal components such as relativistic electrons of low energy or proton cosmic rays.

It is widely accepted that feedback from active galactic nuclei (AGN) plays a key role in the evolution of gas in groups and clusters of galaxies. Unequivocal evidence comes from quasi-spherical X-ray cavities observed near cluster centers having sizes ranging from a few to tens of kpc, some containing radio emission. Cavities apparently evolve from the interaction of AGN jets with the intracluster medium (ICM). However, in numerical simulations it has been difficult to create such fat cavities from narrow jets. Ultra-hot thermal jets dominated by kinetic energy typically penetrate deep into the ICM, forming radially elongated cavities at large radii unlike those observed. Here, we study very light jets dominated energetically by relativistic cosmic rays (CRs) with axisymmetric hydrodynamic simulations, investigating the jet evolution both when they are active and when they are later turned off. We find that, when the thermal gas density in a CR-dominated jet is sufficiently low, the jet has a correspondingly low inertia, and thus decelerates quickly in the ICM. Furthermore, CR pressure causes the jet to expand laterally, encounter and displace more decelerating ICM gas, naturally producing fat cavities near cluster centers similar to those observed. Our calculations of cavity formation imply that AGN jets responsible for creating fat X-ray cavities (radio bubbles) are very light, and dominated by CRs. This scenario is consistent with radio observations of Fanaroff-Riley I jets that appear to decelerate rapidly, produce strong synchrotron emission and expand typically at distances of a few kpc from the central AGN.

In Cygnus A and other classical FR II double radio sources, powerful opposing jets from the cores of halo-centered galaxies drive out into the surrounding cluster gas, forming hotspots of shocked and compressed cluster gas at the jet extremities. The moving hotspots are sandwiched between two shocks. An inner-facing shock receives momentum and cosmic rays from the jet and creates additional cosmic rays that form a radio lobe elongated along the jet axis. An outer-facing bow shock moves directly into the undisturbed group or cluster gas, creating a cocoon of shocked gas enclosing the radio lobe. We describe computations that follow the self-consistent dynamical evolution of the shocked cluster gas and the relativistic synchrotron-emitting gas inside the lobes. Relativistic and non-relativistic components exchange momentum by interacting with small magnetic fields having dynamically negligible energy densities. The evolution of Cygnus A is governed almost entirely by cosmic ray energy flowing from the hotspots. Mass flowing into hotspots from the jets is assumed to be small, greatly reducing the mass of gas flowing back along the jet, common in previous calculations, that would disrupt the spatial segregation of synchrotron-loss ages observed inside FR II radio lobes. We compute the evolution of the cocoon when the velocity and cosmic ray luminosity of the hotspots are constant and when they vary with time. If cosmic rays mix with cluster gas in hotspots before flowing into the radio lobe, the thermal gas is heated to mildly relativistic temperatures, producing an unobserved pressure inside the lobe.

We study the physics of electron acceleration at collisionless shocks that move through a plasma containing large-scale magnetic fluctuations. We numerically integrate the trajectories of a large number of electrons, which are treated as test particles moving in the time dependent electric and magnetic fields determined from 2-D hybrid simulations (kinetic ions, fluid electron). The large-scale magnetic fluctuations effect the electrons in a number of ways and lead to efficient and rapid energization at the shock front. Since the electrons mainly follow along magnetic lines of force, the large-scale braiding of field lines in space allows the fast-moving electrons to cross the shock front several times, leading to efficient acceleration. Ripples in the shock front occuring at various scales will also contribute to the acceleration by mirroring the electrons. Our calculation shows that this process favors electron acceleration at perpendicular shocks. The current study is also helpful in understanding the injection problem for electron acceleration by collisionless shocks. It is also shown that the spatial distribution of energetic electrons is similar to in-situ observations (e.g., Bale et al. 1999; Simnett et al. 2005). The process may be important to our understanding of energetic electrons in planetary bow shocks and interplanetary shocks, and explaining herringbone structures seen in some type II solar radio bursts.

Diffusive shock acceleration at collisionless shocks is thought to be the source of many of the energetic particles observed in space. Large-scale spatial variations of the magnetic field has been shown to be important in understanding observations. The effects are complex, so here we consider a simple, illustrative model. Here, we solve numerically the Parker transport equation for a shock in the presence of large-scale sinusoidal magnetic-field variations. We demonstrate that the familiar planar-shock results can be significantly altered as a consequence of large-scale, meandering magnetic lines of force. Because perpendicular diffusion coefficient $\kappa_\perp$ is generally much smaller than parallel diffusion coefficient $\kappa_\parallel$, the energetic charged particles are trapped and preferentially accelerated along the shock front in the regions where the connection points of magnetic field lines intersecting the shock surface converge, and thus create the "hot spots" of the accelerated particles. For the regions where the connection points separate from each other, the acceleration to high energies will be suppressed. Further, the particles diffuse away from the "hot spot" regions and modify the spectra of downstream particle distribution. These features are qualitatively similar to the recent Voyager's observation in the Heliosheath. These results are potentially important for particle acceleration at shocks propagating in turbulent magnetized plasmas as well as those which contain large-scale nonplanar structures. Examples include anomalous cosmic rays accelerated by the solar wind termination shock, energetic particles observed in propagating heliospheric shocks, and galactic cosmic rays accelerated by supernova blast waves, etc.

Recent X-ray observations of galaxy clusters suggest that cluster populations are bimodally distributed according to central gas entropy and are separated into two distinct classes: cool core (CC) and non-cool core (NCC) clusters. While it is widely accepted that AGN feedback plays a key role in offsetting radiative losses and maintaining many clusters in the CC state, the origin of NCC clusters is much less clear. At the same time, a handful of extremely powerful AGN outbursts have recently been detected in clusters, with a total energy ~10^{61}-10^{62} erg. Using two dimensional hydrodynamic simulations, we show that if a large fraction of this energy is deposited near the centers of CC clusters, which is likely common due to dense cores, these AGN outbursts can completely remove CCs, transforming them to NCC clusters. Our model also has interesting implications for cluster abundance profiles, which usually show a central peak in CC systems. Our calculations indicate that during the CC to NCC transformation, AGN outbursts efficiently mix metals in cluster central regions, and may even remove central abundance peaks if they are not broad enough. For CC clusters with broad central abundance peaks, AGN outbursts decrease peak abundances, but can not effectively destroy the peaks. Our model may simultaneously explain the contradictory (possibly bimodal) results of abundance profiles in NCC clusters, some of which are nearly flat, while others have strong central peaks similar to those in CC clusters. A statistical analysis of the sizes of central abundance peaks and their redshift evolution may shed interesting insights on the origin of both types of NCC clusters and the evolution history of thermodynamics and AGN activity in clusters.

Recent X-ray observations of hot gas in the galaxy cluster MS 0735.6+7421 reveal huge radio-bright, quasi-bipolar X-ray cavities having a total energy ~10^{62} ergs, the most energetic AGN outburst currently known. We investigate the evolution of this outburst with two-dimensional axisymmetric gasdynamical calculations in which the cavities are inflated by relativistic cosmic rays. Many key observational features of the cavities and associated shocks are successfully reproduced. The radial elongation of the cavities indicates that cosmic rays were injected into the cluster gas by a (jet) source moving out from the central AGN. AGN jets of this magnitude must be almost perfectly identically bipolar. The relativistic momentum of a single jet would cause a central AGN black hole of mass 10^9 M_{sun} to recoil at ~6000 km s^{-1}, exceeding kick velocities during black hole mergers, and be ejected from the cluster-center galaxy. When the cavity inflation is complete, 4PV underestimates the total energy received by the cluster gas. Deviations of the cluster gas from hydrostatic equilibrium are most pronounced during the early cavity evolution when the integrated cluster mass found from the observed gas pressure gradient can have systematic errors near the cavities of ~10-30%. The creation of the cavity with cosmic rays generates a long-lasting global cluster expansion that reduces the total gas thermal energy below that received from the cavity shock. One Gyr after this single outburst, a gas mass of ~ 6 \times 10^{11} M_{sun} is transported out beyond a cluster radius of 500 kpc. Such post-cavity outflows can naturally produce the discrepancy observed between the cluster gas mass fraction and the universal baryon fraction inferred from WMAP observations. (Abridged)

We investigate a series of steady-state models of galaxy clusters, in which the hot intracluster gas is efficiently heated by active galactic nucleus (AGN) feedback and thermal conduction, and in which the mass accretion rates are highly reduced compared to those predicted by the standard cooling flow models. We perform a global Lagrangian stability analysis. We show for the first time that the global radial instability in cool core clusters can be suppressed by the AGN feedback mechanism, provided that the feedback efficiency exceeds a critical lower limit. Furthermore, our analysis naturally shows that the clusters can exist in two distinct forms. Globally stable clusters are expected to have either: 1) cool cores stabilized by both AGN feedback and conduction, or 2) non-cool cores stabilized primarily by conduction. Intermediate central temperatures typically lead to globally unstable solutions. This bimodality is consistent with the recently observed anticorrelation between the flatness of the temperature profiles and the AGN activity (Dunn & Fabian 2008) and the observation by Rafferty et al. (2008) that the shorter central cooling times tend to correspond to significantly younger AGN X-ray cavities.

The origin of the bimodality in cluster core entropy is still unknown. At the same time, recent work has shown that thermal conduction in clusters is likely a time-variable phenomenon. We consider if time-variable conduction and AGN outbursts could be responsible for the cool-core (CC), non cool-core (NCC) dichotomy. We show that strong AGN heating can bring a CC cluster to a NCC state, which can be stably maintained by conductive heating from the cluster outskirts. On the other hand, if conduction is shut off by the heat-flux driven buoyancy instability, then the cluster will cool to the CC state again, where it is stabilized by low-level AGN heating. Thus, the cluster cycles between CC and NCC states. In contrast with massive clusters, we predict the CC/NCC bimodality should vanish in groups, due to the lesser role of conductive heating there. We find tentative support from the distribution of central entropy in groups, though firm conclusions require a larger sample carefully controlled for selection effects.

In this paper we address the consumption of energy by astronomers while performing their professional duties. Although we find that astronomy uses a negligible fraction of the US energy budget, the rate at which energy is consumed by an average astronomer is similar to that of a typical high-flying businessperson. We review some of the ways in which astronomers are already acting to reduce their energy consumption. In the coming decades, all citizens will have to reduce their energy consumption to conserve fossil fuel reserves and to help avert a potentially catastrophic change in the Earth's climate. The challenges are the same for astronomers as they are for everyone: decreasing the distances we travel and investing in energy-efficient infrastructure. The high profile of astronomy in the media, and the great public interest in our field, can play a role in promoting energy-awareness to the wider population. Our specific recommendations are therefore to 1) reduce travel when possible, through efficient meeting organization, and by investing in high-bandwidth video conference facilities and virtual-world software, 2) create energy-efficient observatories, computing centers and workplaces, powered by sustainable energy resources, and 3) actively publicize these pursuits.

We investigate a series of steady-state models of galaxy clusters, in which the hot intracluster gas is efficiently heated by active galactic nucleus (AGN) feedback and thermal conduction, and in which the mass accretion rates are highly reduced compared to those predicted by the standard cooling flow models. We perform a global Lagrangian stability analysis. We show for the first time that the global radial instability in cool core clusters can be suppressed by the AGN feedback mechanism, provided that the feedback efficiency exceeds a critical lower limit. Furthermore, our analysis naturally shows that the clusters can exist in two distinct forms. Globally stable clusters are expected to have either: 1) cool cores stabilized by both AGN feedback and conduction, or 2) non-cool cores stabilized primarily by conduction. Intermediate central temperatures typically lead to globally unstable solutions. This bimodality is consistent with the recently observed anticorrelation between the flatness of the temperature profiles and the AGN activity (Dunn & Fabian 2008) and the observation by Rafferty et al. (2008) that the shorter central cooling times tend to correspond to significantly younger AGN X-ray cavities.

Recent observations show that the cooling flows in the central regions of galaxy clusters are highly suppressed. Observed AGN-induced cavities/bubbles are a leading candidate for suppressing cooling, usually via some form of mechanical heating. At the same time, observed X-ray cavities and synchrotron emission point toward a significant non-thermal particle population. Previous studies have focused on the dynamical effects of cosmic-ray pressure support, but none have built successful models in which cosmic-ray heating is significant. Here we investigate a new model of AGN heating, in which the intracluster medium is efficiently heated by cosmic-rays, which are injected into the ICM through diffusion or the shredding of the bubbles by Rayleigh-Taylor or Kelvin-Helmholtz instabilities. We include thermal conduction as well. Using numerical simulations, we show that the cooling catastrophe is efficiently suppressed. The cluster quickly relaxes to a quasi-equilibrium state with a highly reduced accretion rate and temperature and density profiles which match observations. Unlike the conduction-only case, no fine-tuning of the Spitzer conduction suppression factor f is needed. The cosmic ray pressure, P_c/P_g <~ 0.1 and dP_c/dr <~ 0.1 \rho g, is well within observational bounds. Cosmic ray heating is a very attractive alternative to mechanical heating, and may become particularly compelling if GLAST detects the gamma-ray signature of cosmic-rays in clusters.