Search

Shocks driven by coronal mass ejections (CMEs) are primary drivers of gradual solar energetic particle (SEP) events, posing significant risks to space technology and astronauts. Concurrently, particles accelerated at these shocks may also propagate back to the Sun, potentially generating gamma-ray emissions through pion decay. We incorporated advanced modeling and multi-messenger observations to explore the role of CME-driven shocks in gamma-ray emissions and SEPs. Motivated by Fermi-LAT long-duration solar flares, we used the AWSoM MHD model to investigate the connection between the shocks and the properties of observed gamma-ray emissions. By coupling the AWSoM with iPATH model, we evaluate the impact of shock evolution complexity near the Sun on SEP intensity and spectra. Our result points to the importance of accurate background coronal and solar wind modeling, as well as detailed observations of CME source regions, in advancing our understanding of CME-driven shocks and the dynamics of associated energetic particles., Comment: 11 pages, 6 figures, to appear in the Proceedings of IAU Symposium No. 388 - Solar and Stellar Coronal Mass Ejections

Gamma-ray bursts (GRBs) can be probes of the early universe, but currently, only 26% of GRBs observed by the Neil Gehrels Swift Observatory GRBs have known redshifts ($z$) due to observational limitations. To address this, we estimated the GRB redshift (distance) via a supervised statistical learning model that uses optical afterglow observed by Swift and ground-based telescopes. The inferred redshifts are strongly correlated (a Pearson coefficient of 0.93) with the observed redshifts, thus proving the reliability of this method. The inferred and observed redshifts allow us to estimate the number of GRBs occurring at a given redshift (GRB rate) to be 8.47-9 $yr^{-1} Gpc^{-1}$ for $1.9

Gamma-Ray Bursts (GRBs), due to their high luminosities are detected up to redshift 10, and thus have the potential to be vital cosmological probes of early processes in the universe. Fulfilling this potential requires a large sample of GRBs with known redshifts, but due to observational limitations, only 11\% have known redshifts ($z$). There have been numerous attempts to estimate redshifts via correlation studies, most of which have led to inaccurate predictions. To overcome this, we estimated GRB redshift via an ensemble supervised machine learning model that uses X-ray afterglows of long-duration GRBs observed by the Neil Gehrels Swift Observatory. The estimated redshifts are strongly correlated (a Pearson coefficient of 0.93) and have a root mean square error, namely the square root of the average squared error $\langle\Delta z^2\rangle$, of 0.46 with the observed redshifts showing the reliability of this method. The addition of GRB afterglow parameters improves the predictions considerably by 63\% compared to previous results in peer-reviewed literature. Finally, we use our machine learning model to infer the redshifts of 154 GRBs, which increase the known redshifts of long GRBs with plateaus by 94\%, a significant milestone for enhancing GRB population studies that require large samples with redshift., Comment: 12 Figures, 24 pages. Accepted for publication at The Astrophysical Journal Supplement Series

Understanding the nature of energetic particles in the solar atmosphere is one of the most important outstanding problems in heliophysics. Flare-accelerated particles compose a huge fraction of the flare energy budget; they have large influences on how events develop; they are an important source of high-energy particles found in the heliosphere; and they are the single most important corollary to other areas of high-energy astrophysics. Despite the importance of this area of study, this topic has in the past decade received only a small fraction of the resources necessary for a full investigation. For example, NASA has selected no new Explorer-class instrument in the past two decades that is capable of examining this topic. The advances that are currently being made in understanding flare-accelerated electrons are largely undertaken with data from EOVSA (NSF), STIX (ESA), and NuSTAR (NASA Astrophysics). This is despite the inclusion in the previous Heliophysics decadal survey of the FOXSI concept as part of the SEE2020 mission, and also despite NASA's having invested heavily in readying the technology for such an instrument via four flights of the FOXSI sounding rocket experiment. Due to that investment, the instrumentation stands ready to implement a hard X-ray mission to investigate flare-accelerated electrons. This white paper describes the scientific motivation for why this venture should be undertaken soon., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 15 pages, 5 figures

The quiet Sun, i.e. in its non-flaring state or non-flaring regions, emits thermal radiation from radio to ultraviolet. The quiet Sun produces also non-thermal radiation observed in gamma rays due to interactions of Galactic Cosmic Rays (GCR) with the solar gas and photons. We report on a new component: the synchrotron emission by GCR electrons in the solar magnetic field. To the best of our knowledge this is the first time this emission has been theoretically claimed and modeled. We find that the measured GCR electrons with energies from tens of GeV to a few TeV produce synchrotron emission in X-rays, which is a few orders of magnitude lower than current upper limits of the quiet Sun set by RHESSI and FOXSI. For a radially decreasing solar magnetic field we find the expected synchrotron intensity to be almost constant in the solar disk, to peak in the close proximity of the Sun, and to quickly drop away from the Sun. We also estimate the synchrotron emission from radio to gamma rays and we compare it with current observations, especially with LOFAR. While it is negligible from radio to UV compared to the solar thermal radiation, this emission can potentially be observed at high energies with NuSTAR and more promising future FOXSI observations. This could potentially allow for constraining CR densities and magnetic-field intensities at the Sun. This study provides a more complete description and a possible new way for understanding the quite Sun and its environment., Comment: 8 pages, 4 figures

Gamma rays are produced by cosmic ray (CR) protons interacting with the particles at solar photosphere and by cosmic ray electrons and positrons (CRes) via inverse Compton scattering of solar photons. The former come from the solar disk while the latter extend beyond the disk. Evaluation of these emissions requires the flux and spectrum of CRs in the vicinity of the Sun, while most observations provide flux and spectra near the Earth, at around 1 AU from the Sun. Past estimates of the quiet Sun gamma-ray emission use phenomenological modulation procedures to estimate spectra near the Sun (see review by Orlando and Strong 2021 and references therein). We show that CRe transport in the inner heliosphere requires a kinetic approach and use a novel approximation to determine the variation of CRe flux and spectrum from 1 AU to the Sun including effects of (1) the structure of large scale magnetic field, (2) small scale turbulence in the solar wind from several in situ measurements, in particular, those by Parker Solar Probe that extend this information to 0.1 AU, and (3) most importantly, energy losses due to synchrotron and inverse Compton processes. We present results on the flux and spectrum variation of CRes from 1 AU to the Sun for several transport models. In forthcoming papers we will use these results for a more accurate estimate of quiet Sun inverse Compton gamma-ray spectra, and, for the first time, the spectrum of extreme ultraviolet to hard X-ray photons produced by synchrotron emission. These can be compared with the quiet Sun gamma-ray observation by Fermi (see, e.g.~Fermi-LAT Collaboration, 2011) and X-ray upper limits set by RHESSI (Hannah et al., 2010)., Comment: 25 pages, 11 figures, in press Astrophysical Journal

An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all sources. The new reported AGNs include 587 blazar candidates and four radio galaxies. We describe the properties of the new sample and outline changes affecting the previously published one. We also introduce two new parameters in this release, namely the peak energy of the SED high-energy component and the corresponding flux. These parameters allow an assessment of the Compton dominance, the ratio of the Inverse-Compton to the synchrotron peak luminosities, without relying on X-ray data., Comment: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal Supplement Series. Fits files are available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/4LACDR3/ or alternatively at https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-h.fits and https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-l.fits

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources., Comment: accepted in ApJS; follow-up paper to 1902.10045

High-energy solar flares have shown to have at least two distinct phases: prompt-impulsive and delayed-gradual. Identifying the mechanism responsible for accelerating the electrons and ions and the site at which it occurs during these two phases is one of the outstanding questions in solar physics. Many advances have been made over the past decade thanks to new observational data and refined simulations that together help to shed light on this topic. For example, the detection by Fermi Large Area Telescope (LAT) of GeV emission from solar flares originating from behind the visible solar limb and >100 MeV emission lasting for more than 20 hours have suggested the need for a spatially extended source of acceleration during the delayed emission phase. In this work we will review some of the major results from Fermi-LAT observations of the 24th solar cycle and how this new observational channel combined with observations from across the electromagnetic spectrum can provide a unique opportunity to diagnose the mechanisms of high-energy emission and particle acceleration in solar flares., Comment: 21 pages, 15 figures, 37th International Cosmic Ray Conference (ICRC2021) proceedings

The Large Area Telescope (LAT) on board the \Fermi Gamma-ray Space Telescope (\Fermi) shows long-lasting high-energy emission in many gamma-ray bursts (GRBs), similar to X-ray afterglows observed by the Neil Gehrels Swift Observatory \citep[\textit{Swift};][]{gehrels2004}. Some LAT light curves (LCs) show a late-time flattening reminiscent of X-ray plateaus. We explore the presence of plateaus in LAT temporally extended emission analyzing GRBs from the second \lat GRB Catalog \citep[2FLGC;][]{Ajello2019apj} from 2008 to May 2016 with known redshifts, and check whether they follow closure relations corresponding to 4 distinct astrophysical environments predicted by the external forward shock (ES) model. We find that three LCs can be fit by the same phenomenological model used to fit X-ray plateaus \citep{Willingale2007} and show tentative evidence for the existence of plateaus in their high-energy extended emission. The most favorable scenario is a slow cooling regime, whereas the preferred density profile for each GRBs varies from a constant density ISM to a $r^{-2}$ wind environment. We also compare the end time of the plateaus in $\gamma$-rays and X-rays using a statistical comparison with 222 \textit{Swift} GRBs with plateaus and known redshifts from January 2005 to August 2019. Within this comparison, the case of GRB 090510 shows an indication of chromaticity at the end time of the plateau. Finally, we update the 3-D fundamental plane relation among the rest frame end time of the plateau, its correspondent luminosity, and the peak prompt luminosity for 222 GRBs observed by \textit{Swift}. We find that these three LAT GRBs follow this relation., Comment: 28 pages, 10 pictures, 3 Tables to appear on ApJ Supplement

Originating from neutron star-neutron star (NS-NS) or neutron star-black hole (NS-BH) mergers, short gamma-ray bursts (SGRBs) are the first electromagnetic emitters associated with gravitational waves. This association makes the determination of SGRB formation rate (FR) a critical issue. We determine the true SGRB FR and its relation to the cosmic star formation rate (SFR). This can help in determining the expected Gravitation Wave (GW) rate involving small mass mergers. We present non-parametric methods for the determination of the evolutions of the luminosity function (LF) and the FR using SGRBs observed by {\it Swift}, without any assumptions. These are powerful tools for small samples, such as our sample of 68 SGRBs. We combine SGRBs with and without extended emission (SEE), assuming that both descend from the same progenitor. To overcome the incompleteness introduced by redshift measurements we use the Kolmogorov-Smirnov (KS) test to find flux thresholds yielding a sample of sources with a redshift drawn from the parent sample including all sources. Using two subsamples of SGRBs with flux limits of $4.57 \times 10^{-7}$ and $2.15 \times 10^{-7}$ erg cm$^{-2}$ s$^{-1}$ with respective KS {\it p=(1, 0.9)}, we find a 3 $\sigma$ evidence for luminosity evolution (LE), a broken power-law LF with significant steepening at $L\sim 10^{50}$ erg s$^{-1}$, and a FR evolution that decreases monotonically with redshift (independent of LE and the thresholds). Thus, SGRBs may have been more luminous in the past with a FR delayed relative to the SFR as expected in the merger scenario., Comment: 11 pages, 13 figures contained in 4 panel, 1 table. The paper is about to appear on ApJL

Multi-wavelength analyses of spectra of active galactic nuclei (AGNs) provide useful information on the physical processes in the accretion disk and jets of black holes. This, however, is limited to bright sources and may not represent the population as a whole. Another approach is through the investigation of the cosmological evolution of the luminosity function (LF) that shows varied evolution (luminosity and density) at different wavelengths. These differences and the correlations between luminosities can shed light on the Jet-accretion disk connection. Most such studies use forward fitting parametric methods that involve several functions and many parameters. We use non-parametric, non-binning methods developed by Efron and Petrosian, and Lynden-Bell, for obtaining unbiased description of the evolution of the LF, from data truncated by observational selection effects. We present analysis of the evolution of gamma-ray LF of blazars with main focus on flat spectrum radio quasars (FSRQs). This requires analysis of both gamma-ray and optical data, essential for redshift measurements, and a description of the joint LF. We use a new approach which divides the sample into two sub-samples, each with its own flux limit. We use the Fermi-LAT and GAIA observations, and present results on the gamma-ray LF and its evolution, and determine the intrinsic correlation between the gamma-ray and optical luminosities corrected for the well known false correlation induced by their similar redshift dependence and evolution of the two luminosities. We also present a direct estimation of the contribution of blazars to the spectrum of the extragalactic gamma-ray background., Comment: 46 pages, 11 figures, 2 tables, accepted for publication in ApJ

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity., Comment: submitted to Astrophysics Decadal Survey as a State of the Profession white paper

Gamma-ray bursts (GRBs) are spectacularly energetic events, with the potential to inform on the early universe and its evolution, once their redshifts are known. Unfortunately, determining redshifts is a painstaking procedure requiring detailed follow-up multi-wavelength observations often involving various astronomical facilities, which have to be rapidly pointed at these serendipitous events. Here we use Machine Learning algorithms to infer redshifts from a collection of observed temporal and spectral features of GRBs. We obtained a very high correlation coefficient ($0.96$) between the inferred and the observed redshifts, and a small dispersion (with a mean square error of $0.003$) in the test set. The addition of plateau afterglow parameters improves the predictions by $61.4\%$ compared to previous results. The GRB luminosity function and cumulative density rate evolutions, obtained from predicted and observed redshift are in excellent agreement indicating that GRBs are effective distance indicators and a reliable step for the cosmic distance ladder., Comment: 22 pages, 5 figures to be submitted

A knowledge of the particle escape time from the acceleration regions of many space and astrophysical sources is of critical importance in the analysis of emission signatures produced by these particles and in the determination of the acceleration and transport mechanisms at work. This paper addresses this general problem, in particular in solar flares, where in addition to scattering by turbulence, the magnetic field convergence from the acceleration region towards its boundaries also influences the particle escape. We test an (approximate) analytic relation between escape and scattering times, and the field convergence rate, based on the work of Malyshkin and Kulsrud (2001), valid for both strong and weak diffusion limits and isotropic pitch angle distribution of the injected particles, with a numerical model of particle transport. To this end, a kinetic Fokker-Planck transport model of particles is solved with a stochastic differential equation scheme assuming different initial pitch angle distributions. This approach enables further insights into the phase-space dynamics of the transport process, which would otherwise not be accessible. We find that in general the numerical results agree well with the analytic equation for the isotropic case, however, there are significant differences weak diffusion regime for non-isotopic cases, especially for distributions beamed along the magnetic field lines. The results are important in the interpretation of observations of energetic particles in solar flares, and other similar space and astrophysical acceleration sites and for the determination of acceleration-transport coefficients, commonly used in Fokker-Planck type kinetic equations., Comment: 8 pages, 2 figures, submitted

The Fermi gamma-ray Space Telescope (Fermi) has detected hard X-ray (HXR) and gamma-ray photons from three flares, which according to \stereo occurred in active regions behind the limb of the Sun as delineated by near Earth instruments. For two of these flares \r has provided HXR images with sources located just above the limb, presumably from the loop top (LT) region of a relatively large loop. Fermi-Gamma-ray Burst Monitor has detected HXRs and gamma-rays, and RSTN has detected microwaves emissions with similar light curves. This paper presents a quantitative analysis of these multi-wavelength observations assuming that HXRs and microwaves are produced by electrons accelerated at the LT source, with emphasize on the importance of the proper treatment of escape of the particles from the acceleration-source region and the trans-relativistic nature of the analysis. The observed spectra are used to determine the magnetic field and relativistic electron spectra. It is found that a simple power-law in momentum (with cut off above a few 100 MeV) agrees with all observations, but in energy space a broken power law spectrum (steepening at rest mass energy) may be required. It is also shown that the production of the $>100$ MeV photons detected by The Fermi-Large Area Telescope at the LT source would require more energy compared to photospheric emission. These energies are smaller than that required for electrons, so that the possibility that all the emissions originate in the LT cannot be ruled out on energetic grounds. However, the differences in the light curves and emission centroids of HXRs and $>100$ MeV gamma-rays favour a different source for the latter., Comment: 14 pages with 10 figures, In press ApJ

The discovery of the Fermi bubbles---a huge bilobular structure seen in GeV gamma-rays above and below the Galactic center---implies the presence of a large reservoir of high energy particles at $\sim 10 \, \text{kpc}$ from the disk. The absence of evidence for a strong shock coinciding with the edge of the bubbles, and constraints from multi-wavelength observations point towards stochastic acceleration by turbulence as a likely mechanism of acceleration. We have investigated the time-dependent acceleration of electrons in a large-scale outflow from the Galactic centre. For the first time, we present a detailed numerical solution of the particle kinetic equation that includes the acceleration, transport and relevant energy loss processes. We also take into account the addition of shock acceleration of electrons at the bubble's blast wave. Fitting to the observed spectrum and surface brightness distribution of the bubbles allows determining the transport coefficients, thereby shedding light on the origin of the Fermi bubbles., Comment: 14 pages, 11 figures; submitted to A&A

Recent detections of high-energy $\gamma$-rays from behind-the-limb (BTL) solar flares by the \emph{Fermi $\gamma$-ray Space Telescope} pose a puzzle and challenge on the particle acceleration and transport mechanisms. In such events, the $\gamma$-ray emission region is located away from the BTL flare site by up to tens of degrees in heliogrpahic longitude. It is thus hypothesized that particles are accelerated at the shock driven by the coronal mass ejection (CME) and then travel from the shock downstream back to the front side of the Sun to produce the observed $\gamma$-rays. To test this scenario, we performed data-driven, global magnetohydrodynamics simulations of the CME associated with a well-observed BTL flare on 2014 September 1. We found that part of the CME-driven shock develops magnetic connectivity with the $\gamma$-ray emission region, facilitating transport of particles back to the Sun. Moreover, the observed increase in $\gamma$-ray flux is temporally correlated with (1) the increase of the shock compression ratio and (2) the presence of a quasi-perpendicular shock over the area that is magnetically connected to the $\gamma$-ray emitting region, both conditions favoring the diffusive shock acceleration (DSA) of particles. These results support the above hypothesis and can help resolve another puzzle, i.e., long-duration (up to 20 hours) $\gamma$-rays flares. We suggest that, in addition to DSA, stochastic acceleration by plasma turbulence may also play a role, especially in the shock downstream region and during the early stage when the shock Alfv\'{e}n Mach number is small., Comment: 29 pages, 8 figures, accepted by ApJ

Flares close to the solar limb, where the footpoints are occulted, can reveal the spectrum and structure of the coronal loop-top source in X-rays. We aim at studying the properties of the corresponding energetic electrons near their acceleration site, without footpoint contamination. To this end, a statistical study of partially occulted flares observed with RHESSI is presented here, covering a large part of solar cycles 23 and 24. We perform a detailed spectra, imaging and light curve analysis for 116 flares and include contextual observations from SDO and STEREO when available, providing further insights into flare emission that was previously not accessible. We find that most spectra are fitted well with a thermal component plus a broken power-law, non-thermal component. A thin-target kappa distribution model gives satisfactory fits after the addition of a thermal component. X-rays imaging reveals small spatial separation between the thermal and non-thermal components, except for a few flares with a richer coronal source structure. A comprehensive light curve analysis shows a very good correlation between the derivative of the soft X-ray flux (from GOES) and the hard X-rays for a substantial number of flares, indicative of the Neupert effect. The results confirm that non-thermal particles are accelerated in the corona and estimated timescales support the validity of a thin-target scenario with similar magnitudes of thermal and non-thermal energy fluxes., Comment: 17 pages, 11 figures, accepted for publication in ApJ

At hard X-ray energies, the bright footpoint emission from solar flare loops often prevents a detailed analysis of the weaker loop-top source morphology due to the limited dynamic range available for X-ray imaging. Here, we study the X1.3 April 25, 2014 flare with the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). This partially occulted limb flare allows the analysis of the loop-top emission in isolation. We present results on the flare light curve at different energies, the source morphology from X-ray imaging and a detailed spectral analysis of the different source components by imaging spectroscopy. The loop-top source, a likely site of particle acceleration, shows a clear composition of different emission components. The results indicate the opportunities that detailed imaging of hard X-rays can provide to learn about particle acceleration, transport and heating processes in solar flares., Comment: 5 pages, 4 figures, submitted as conference proceedings for the 15th Annual International Astrophysics Conference

Observations relating the characteristics of electrons seen near Earth (SEPs) and those producing flare radiation show that in certain (prompt) events the origin of both population appears to be the flare site, which show strong correlation between the number and spectral index of SEP and hard X-ray radiating electrons, but in others(delayed), which are associated with fast CMEs, this relation is complex and SEPs tend to be harder. Prompt event spectral relation disagrees with that expected in thick or thin target models. We show that using a a more accurate treatment of the transport of the accelerated electrons to the footpoints and to the Earth can account for this discrepancy. Our results are consistent with those found by Chen and Petrosian (2013) for two flares using non-parametric inversion methods, according to which we have weak diffusion conditions, and trapping mediated by magnetic field convergence. The weaker correlations and harder spectra of delayed events can come about by re-acceleration of electrons in the CME shock environment. We describe under what conditions such a hardening can be achieved. Using this (acceleration at the flare and re-acceleration in the CME) scenario we show that we can describe the similar dichotomy that exists between the so called impulsive, highly enriched ($^3$He and heavy ions) and softer SEP events, and stronger more gradual SEP events with near normal ionic abundances and harder spectra. These methods can be used to distinguish the acceleration mechanisms and to constrain their characteristics., Comment: 26 pages, 10 figures, ApJ in press

We investigate heating and acceleration of protons from a thermal gas with a generic diffusion and acceleration model, and subject to Coulomb scattering and energy loss, as was carried out in Petrosian & East (2008) for electrons. As protons gain energy their loss to electrons becomes important. Thus, we need to solve the coupled proton-electron kinetic equation. We numerically solve the coupled Fokker-Plank equations and computes the time evolution of the spectra of both particles. We show that this can lead to a quasi-thermal component plus a high energy nonthermal tail. We determine the evolution of nonthermal tail and the quasi-thermal component. The results may be used to explore the possibility of inverse bremsstrahlung radiation as a source of hard X-ray emissions from hot sources such as solar flares, accretion disk coronas and the intracluster medium of galaxy clusters. We find that emergence of nonthermal protons is accompanied by excessive heating of the entire plasma, unless the turbulence needed for scattering and acceleration is steeper than Kolmogorov and the acceleration parameters, the duration of the acceleration, and/or the initial distributions are significantly fine-tuned. These results severely constraint the feasibility of nonthermal inverse bremsstrahlung process producing hard X-ray emissions. However the nonthermal tail may be the seed particles for further re-acceleration to relativistic energies, say by a shock. In the Appendix we present some tests of the integrity of the algorithm used and present a new formula for the energy loss rate due to inelastic proton-proton interactions., Comment: 21 pages, 9 figures, accepted for publication in ApJ

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in MgII h&k, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DST/IBIS) in H$\alpha$ 6563 \AA\ and CaII 8542 \AA, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multi-threaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated MgII line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of $\approx$800 km below the transition region can turn the simulated MgII line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares., Comment: 14 pages, 18 figures, accepted in ApJ

Fermi-LAT >30 MeV observations have increased the number of detected solar flares by almost a factor of 10 with respect to previous space observations. These sample both the impulsive and long duration phases of GOES M and X class flares. Of particular interest is the recent detections of three solar flares whose position behind the limb was confirmed by the STEREO-B spacecraft. While gamma-ray emission up to tens of MeV resulting from proton interactions has been detected before from occulted solar flares, the significance of these particular events lies in the fact that these are the first detections of >100 MeV gamma-ray emission from footpoint-occulted flares. We will present the Fermi-LAT, RHESSI and STEREO observations of these flares and discuss the various emission scenarios for these sources and implications for the particle acceleration mechanisms., Comment: 7 pages, 5 figures, in Proceedings of the 34th International Cosmic Ray Conference (ICRC 2015), The Hague (The Netherlands)

We report the first detection of >100 MeV gamma rays associated with a behind-the-limb solar flare, which presents a unique opportunity to probe the underlying physics of high-energy flare emission and particle acceleration. On 2013 October 11 a GOES M1.5 class solar flare occurred ~ 9.9 degrees behind the solar limb as observed by STEREO-B. RHESSI observed hard X-ray emission above the limb, most likely from the flare loop-top, as the footpoints were occulted. Surprisingly, the Fermi Large Area Telescope (LAT) detected >100 MeV gamma-rays for ~30 minutes with energies up to GeV. The LAT emission centroid is consistent with the RHESSI hard X-ray source, but its uncertainty does not constrain the source to be located there. The gamma-ray spectra can be adequately described by bremsstrahlung radiation from relativistic electrons having a relatively hard power-law spectrum with a high-energy exponential cutoff, or by the decay of pions produced by accelerated protons and ions with an isotropic pitch-angle distribution and a power-law spectrum with a number index of ~3.8. We show that high optical depths rule out the gamma rays originating from the flare site and a high-corona trap model requires very unusual conditions, so a scenario in which some of the particles accelerated by the CME shock travel to the visible side of the Sun to produce the observed gamma rays may be at work., Comment: 6 pages, 3 figures, accepted for publication on ApJ Letters

Solar flares involve complex processes that are coupled and span a wide range of temporal, spatial, and energy scales. Modeling such processes self-consistently has been a challenge in the past. Here we present results from simulations that couple particle kinetics with hydrodynamics of the atmospheric plasma. We combine the Stanford unified Fokker-Planck code that models particle acceleration and transport with the RADYN hydrodynamic code that models the atmospheric response to collisional heating by accelerated electrons through detailed radiative transfer calculations. We perform simulations using two different electron spectra, one an {\it ad hoc} power law and the other predicted by the model of stochastic acceleration by turbulence or plasma waves. Surprisingly, the later model, even with energy flux $\ll 10^{10}$ erg s$^{-1}$ cm$^{-2}$, can cause "explosive" chromospheric evaporation and drive stronger up- and downflows (and hydrodynamic shocks). This is partly because our acceleration model, like many others, produces a spectrum consisting of a quasi-thermal component plus a power-law tail. We synthesize emission line profiles covering different heights in the lower atmosphere, including H$\alpha$ 6563 \AA, HeII 304 \AA, CaII K 3934 \AA\ and SiIV 1393 \AA. One interesting result is the unusual high temperature (up to a few $10^5$ K) of the formation site of HeII 304 \AA, which is expected due to photoionization-recombination under flare conditions, compared to those in the quiet Sun dominated by collisional excitation. When compared with observations, our results can constrain the properties of non-thermal electrons and thus the poorly understood particle acceleration mechanism., Comment: 17 pages, 18 figures. Accepted in ApJ

Gamma-ray bursts (GRBs) by virtue of their high luminosities can be detected up to very high redshifts and therefore can be excellent probes of the early universe. This task is hampered by the fact that most of their characteristics have a broad range so that we first need to obtain an accurate description of the distribution of these characteristics, and specially, their cosmological evolution. We use a sample of about 200 \swift long GRBs with known redshift to determine the luminosity and formation rate evolutions and the general shape of the luminosity function. In contrast to most other forward fitting methods of treating this problem we use the Efron Petrosian methods which allow a non-parametric determination of above quantities. We find a relatively strong luminosity evolution, a luminosity function that can be fitted to a broken power law, and an unusually high rate of formation rate at low redshifts, a rate more than one order of magnitude higher than the star formation rate (SFR). On the other hand, our results seem to agree with the almost constant SFR in redshifts 1 to 3 and the decline above this redshift., Comment: 12 pages, six figures, accepted for publication in ApJ

Active galactic nuclei jets are detected via their radio and/or gamma-ray emissions while the accretion disks are detected by their optical and UV radiation. Observations of the radio and optical luminosities show a strong correlation between the two luminosities. However, part of this correlation is due to the redshift or distances of the sources that enter in calculating the luminosities from the observed fluxes and part of it could be due to the differences in the cosmological evolution of luminosities. Thus, the determination of the intrinsic correlations between the luminosities is not straightforward. It is affected by the observational selection effects and other factors that truncate the data, sometimes in a complex manner (e.g. Antonucci (2011) and Pavildou et al. (2010)). In this paper we describe methods that allow us to determine the evolution of the radio and optical luminosities, and determine the true intrinsic correlation between the two luminosities. We find a much weaker correlation than observed and sub-linear relations between the luminosities. This has a significant implication for jet and accretion disk physics., Comment: 7 pages, 4 figures, published in Proc IAU #314, F. Massaro, C.C. Cheung, E. Lopez, A. Siemiginowska, eds

Solar flares involve impulsive energy release, which results in enhanced radiation over a broad spectral range and a wide range of heights. In particular, line emission from the chromosphere can provide critical diagnostics of plasma heating processes. Thus, a direct comparison between high-resolution spectroscopic observations and advanced numerical modeling results could be extremely valuable, but has not yet been attempted. In this paper, we present in this paper such a self-consistent investigation of an M3.0 flare observed by the Dunn Solar Telescope's Interferometric Bi-dimensional Spectrometer (IBIS) on 2011 September 24 which we have modeled with the radiative hydrodynamic code RADYN. We obtained images and spectra of the flaring region with IBIS in H$\alpha$ 6563 \AA\ and Ca II 8542 \AA, and with the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in X-rays. The latter observations were used to infer the non-thermal electron population, which was passed to RADYN to simulate the atmospheric response to electron collisional heating. We then synthesized spectral lines and compared their shapes and intensities to those observed by IBIS and found a general agreement. In particular, the synthetic Ca II 8542 \AA\ profile fits well to the observed profile, while the synthetic H$\alpha$ profile is fainter in the core than for the observation. This indicates that H$\alpha$ emission is more responsive to the non-thermal electron flux than the Ca II 8542 \AA\ emission. We suggest that it is necessary to refine the energy input and other processes is necessary to resolve this discrepancy., Comment: 16 pages, 13 figures. Accepted in ApJ

We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called $\kappa$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $\kappa$-distribution has a power-law tail, instead of an exponential cutoff, well above the thermal peak momentum. So there are a larger number of potential seed particles with momentum, above that required for participation in the DSA process. As a result, the injection fraction for the $\kappa$-distribution depends on the shock Mach number much less severely compared to that for the Maxwellian distribution. Thus, the existence of $\kappa$-like suprathermal tails at shocks would ease the problem of extremely low injection fractions, especially for electrons and especially at weak shocks such as those found in the intracluster medium. We suggest that the injection fraction for protons ranges $10^{-4}-10^{-3}$ for a $\kappa$-distribution with $10 < \kappa_p < 30$ at quasi-parallel shocks, while the injection fraction for electrons becomes $10^{-6}-10^{-5}$ for a $\kappa$-distribution with $\kappa_e < 2$ at quasi-perpendicular shocks. For such $\kappa$ values the ratio of cosmic ray electrons to protons naturally becomes $K_{e/p}\sim 10^{-3}-10^{-2}$, which is required to explain the observed ratio for Galactic cosmic rays., Comment: 10 pages, 4 figures, submitted to Astrophysical Journal

We have developed an inversion method for determination of the characteristics of the acceleration mechanism directly and non-parametrically from observations, in contrast to the usual forward fitting of parametric model variables to observations. In two recent papers (Petrosian & Chen 2010, Chen & Petrosian 2013), we demonstrate the efficacy of this inversion method by its application to acceleration of electrons in solar flares based on stochastic acceleration by turbulence. Here we explore its application for determining the characteristics of shock acceleration in supernova remnants (SNRs) based on the electron spectra deduced from the observed nonthermal radiation from SNRs and the spectrum of the cosmic ray electrons observed near the Earth. These spectra are related by the process of escape of the electrons from SNRs and energy loss during their transport in the galaxy. Thus, these observations allow us to determine spectral characteristics of the momentum and pitch angle diffusion coefficients, which play crucial roles in both direct acceleration by turbulence and in high Mach number shocks. Assuming that the average electron spectrum deduced from a few well known SNRs is representative of those in the solar neighborhood we find interesting discrepancies between our deduced forms for these coefficients and those expected from well known wave-particle interactions. This may indicate that the standard assumptions made in treatment of shock acceleration need revision. In particular, the escape of particles from SNRs may be more complex than generally assumed., Comment: 16 pages, 4 figures, accepted to PRD

Gamma-ray bursts (GRBs), which have been observed up to redshifts $z \approx 9.5$ can be good probes of the early universe and have the potential of testing cosmological models. The analysis by Dainotti of GRB Swift afterglow lightcurves with known redshifts and definite X-ray plateau shows an anti-correlation between the \underline{rest frame} time when the plateau ends (the plateau end time) and the calculated luminosity at that time (or approximately an anti-correlation between plateau duration and luminosity). We present here an update of this correlation with a larger data sample of 101 GRBs with good lightcurves. Since some of this correlation could result from the redshift dependences of these intrinsic parameters, namely their cosmological evolution we use the Efron-Petrosian method to estimate the luminosity and time evolution and to correct for this effects to determine the intrinsic nature of this correlation., Comment: "7th Huntsville Gamma-Ray Burst Symposium, GRB 2013: paper 25 in eConf Proceedings C1304143"

Gamma-ray bursts (GRBs), which have been observed up to redshifts z approx 9.5 can be good probes of the early universe and have the potential of testing cosmological models. The analysis by Dainotti of GRB Swift afterglow lightcurves with known redshifts and definite X-ray plateau shows an anti-correlation between the rest frame time when the plateau ends (the plateau end time) and the calculated luminosity at that time (or approximately an anti-correlation between plateau duration and luminosity). We present here an update of this correlation with a larger data sample of 101 GRBs with good lightcurves. Since some of this correlation could result from the redshift dependences of these intrinsic parameters, namely their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic and describe how we recover it and how this can be used to constrain physical models of the plateau emission, whose origin is still unknown. The present result could help clarifing the debated issue about the nature of the plateau emission., Comment: Astrophysical Journal accepted

Following our recent paper (Petrosian & Chen 2010), we have developed an inversion method to determine the basic characteristics of the particle acceleration mechanism directly and non-parametrically from observations under the leaky box framework. In the above paper, we demonstrated this method for obtaining the energy dependence of the escape time. Here, by converting the Fokker-Planck equation to its integral form, we derive the energy dependences of the energy diffusion coefficient and direct acceleration rate for stochastic acceleration in terms of the accelerated and escaping particle spectra. Combining the regularized inversion method of Piana et al. 2007 and our procedure, we relate the acceleration characteristics in solar flares directly to the count visibility data from RHESSI. We determine the timescales for electron escape, pitch angle scattering, energy diffusion, and direct acceleration at the loop top acceleration region for two intense solar flares based on the regularized electron flux spectral images. The X3.9 class event shows dramatically different energy dependences for the acceleration and scattering timescales, while the M2.1 class event shows a milder difference. The M2.1 class event could be consistent with the stochastic acceleration model with a very steep turbulence spectrum. A likely explanation of the X3.9 class event could be that the escape of electrons from the acceleration region is not governed by a random walk process, but instead is affected by magnetic mirroring, in which the scattering time is proportional to the escape time and has an energy dependence similar to the energy diffusion time., Comment: 11 pages, 5 figures, accepted for publication in ApJ

Where particle acceleration and plasma heating take place in relation to magnetic reconnection is a fundamental question for solar flares. We report analysis of an M7.7 flare on 2012 July 19 observed by SDO/AIA and RHESSI. Bi-directional outflows in forms of plasmoid ejections and contracting cusp-shaped loops originate between an erupting flux rope and underlying flare loops at speeds of typically 200-300 km/s up to 1050 km/s. These outflows are associated with spatially separated double coronal X-ray sources with centroid separation decreasing with energy. The highest temperature is located near the nonthermal X-ray loop-top source well below the original heights of contracting cusps near the inferred reconnection site. These observations suggest that the primary loci of particle acceleration and plasma heating are in the reconnection outflow regions, rather than the reconnection site itself. In addition, there is an initial ascent of the X-ray and EUV loop-top source prior to its recently recognized descent, which we ascribe to the interplay among multiple processes including the upward development of reconnection and the downward contractions of reconnected loops. The impulsive phase onset is delayed by 10 minutes from the start of the descent, but coincides with the rapid speed increases of the upward plasmoids, the individual loop shrinkages, and the overall loop-top descent, suggestive of an intimate relation of the energy release rate and reconnection outflow speed., Comment: Accepted by ApJ, March 04, 2013; 18 pages, 14 figures, online movies at http://lmsal.com/~weiliu/my_papers/2012-08_flare-cusp-SADs/movies.zip

The subject of this paper is stochastic acceleration by plasma turbulence, a process akin to the original model proposed by Fermi. We review the relative merits of different acceleration models, in particular the so called first order Fermi acceleration by shocks and second order Fermi by stochastic processes, and point out that plasma waves or turbulence play an important role in all mechanisms of acceleration. Thus, stochastic acceleration by turbulence is active in most situations. We also show that it is the most efficient mechanism of acceleration of relatively cool non relativistic thermal background plasma particles. In addition, it can preferentially accelerate electrons relative to protons as is needed in many astrophysical radiating sources, where usually there are no indications of presence of shocks. We also point out that a hybrid acceleration mechanism consisting of initial acceleration by turbulence of background particles followed by a second stage acceleration by a shock has many attractive features. It is demonstrated that the above scenarios can account for many signatures of the accelerated electrons, protons and other ions, in particular $^3$He and $^4$He, seen directly as Solar Energetic Particles and through the radiation they produce in solar flares., Comment: 29 pages 7 figures for proceedings of ISSI-Bern workshop on Particle Acceleration 2011

[Abridged]We present analysis of a pair of unusually energetic coronal hard X-ray (HXR) sources detected by RHESSI during the impulsive phase of an X3.9 class solar flare on 2003 November 3, which simultaneously shows two intense footpoint (FP) sources. A distinct loop top (LT) coronal source is detected up to ~150 keV and a second (upper) coronal source up to ~80 keV. These photon energies are much higher than commonly observed in coronal sources and pose grave modeling challenges. The LT source in general appears higher in altitude with increasing energy and exhibits a more limited motion compared to the expansion of the thermal loop. The high energy LT source shows an impulsive time profile and its nonthermal power law spectrum exhibits soft-hard-soft evolution during the impulsive phase, similar to the FP sources. The upper coronal source exhibits an opposite spatial gradient and a similar spectral slope compared to the LT source. These properties are consistent with the model of stochastic acceleration of electrons by plasma waves or turbulence. However, the LT and FP spectral index difference (varying from ~0-1) is much smaller than commonly measured and than that expected from a simple stochastic acceleration model. Additional confinement or trapping mechanisms of high energy electrons in the corona are required. Comprehensive modeling including both kinetic effects and the macroscopic flare structure may shed light on this behavior. These results highlight the importance of imaging spectroscopic observations of the LT and FP sources up to high energies in understanding electron acceleration in solar flares. Finally, we show that the electrons producing the upper coronal HXR source may very likely be responsible for the type III radio bursts at the decimetric/metric wavelength observed during the impulsive phase of this flare., Comment: 13 emulateapj pages, 10 figures; accepted for publication in ApJ

We examine the effects of time dilation on the temporal profiles of gamma-ray burst (GRB) pulses. By using prescriptions for the shape and evolution of prompt gamma-ray spectra, we can generate a simulated population of single pulsed GRBs at a variety of redshifts and observe how their light curves would appear to a gamma-ray detector here on Earth. We find that the observer frame duration of individual pulses does not increase as a function of redshift as one would expect from the cosmological expansion of a Friedman-Lemaitre-Robertson-Walker Universe. In fact, the duration of individual pulses is seen to decrease as their signal-to-noise decreases with increasing redshift, as only the brightest portion of a high redshift GRB's light curve is accessible to the detector. The results of our simulation are consistent with the fact that a systematic broadening of GRB durations as a function of redshift has not materialized in either the Swift or Fermi detected GRBs with known redshift. We show that this fundamental duration bias implies that the measured durations and associated Eiso estimates for GRBs detected near an instrument's detection threshold should be considered lower limits to their true values. We conclude by predicting that the average peak-to-peak time for a large number of multi-pulsed GRBs as a function of redshift may eventually provide the evidence for time dilation that has so far eluded detection., Comment: 6 pages, 8 figures, submitted to ApJ

We review basic theoretical concepts in particle acceleration, with particular emphasis on processes likely to occur in regions of magnetic reconnection. Several new developments are discussed, including detailed studies of reconnection in three-dimensional magnetic field configurations (e.g., current sheets, collapsing traps, separatrix regions) and stochastic acceleration in a turbulent environment. Fluid, test-particle, and particle-in-cell approaches are used and results compared. While these studies show considerable promise in accounting for the various observational manifestations of solar flares, they are limited by a number of factors, mostly relating to available computational power. Not the least of these issues is the need to explicitly incorporate the electrodynamic feedback of the accelerated particles themselves on the environment in which they are accelerated. A brief prognosis for future advancement is offered., Comment: This is a chapter in a monograph on the physics of solar flares, inspired by RHESSI observations. The individual articles are to appear in Space Science Reviews (2011)

The model of stochastic acceleration of particles by turbulence has been successful in explaining many observed features of solar flares. Here we demonstrate a new method to obtain the accelerated electron spectrum and important acceleration model parameters from the high resolution hard X-ray observations provided by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). In our model, electrons accelerated at or very near the loop top produce thin target bremsstrahlung emission there and then escape downward producing thick target emission at the loop footpoints. Based on the electron flux spectral images obtained by the regularized inversion of the RHESSI count visibilities, we derive several important parameters for the acceleration model. We apply this procedure to the 2003 November 03 solar flare, which shows a loop top source up to 100--150 keV in hard X-ray with a relatively flat spectrum in addition to two footpoint sources. The results imply presence of strong scattering and a high density of turbulence energy with a steep spectrum in the acceleration region., Comment: 4 pages, 2 figures; revised manuscript submitted to ApJL

In recent years there has been considerable activity in using gamma-ray bursts as cosmological probes for determining global cosmological parameters complementing results from type Ia supernovae and other methods. This requires a characteristics of the source to be a standard candle. We show that contrary to earlier indications the accumulated data speak against this possibility. Another method would be to use correlation between a distance dependent and a distance independent variable to measure distance and determine cosmological parameters as is done using Cepheid variables and to some extent Type Ia supernovae. Many papers have dealt with the use of so called Amati relation, first predicted by Lloyd, Petrosian and Mallozzi, or the Ghirlanda relation for this purpose. We have argued that these procedure involve many unjustified assumptions which if not true could invalidate the results. In particular, we point out that many evolutionary effects can affect the final outcome. In particular, we demonstrate that the existing data from Swift and other earlier satellites show that the gamma-ray burst may have undergone luminosity evolution. Similar evolution may be present for other variables such as the peak photon energy of the total radiated energy. Another out come of our analysis is determination of the luminosity function and the comoving rate evolution of gamma-ray bursts which does not seem to agree with the cosmic star formation rate. We caution however, that the above result are preliminary and includes primarily the effect of detection threshold. Other selection effects, perhaps less important than this, are also known to be present and must be accounted for. We intend to address these issues in future publications., Comment: 19 pages, 7 figures

Recent observations of cosmic ray electrons from several instruments have revealed various degrees of deviation in the measured electron energy distribution from a simple power-law, in a form of an excess around TeV energies. An even more prominent deviation has been observed in the fraction of cosmic ray positrons around 100 GeV energies. In this paper we show that the observed excesses in the electron spectrum may be easily re-produced without invoking any unusual sources other than the general diffuse Galactic components of cosmic rays. The primary physical effect involved is the Klein-Nishina suppression of the electron cooling rate around TeV energies. With a very reasonable choice of the model parameters characterizing the local interstellar medium, we can reproduce the most recent observations by Fermi and HESS experiments. We also find that high positron fraction increasing with energy, as claimed by the PAMELA experiment, cannot be explained in our model with the conservative set of the model parameters. We are able, however, to reproduce the PAMELA results assuming high values of the starlight and interstellar gas densities, which would be more appropriate for vicinities of supernova remnants. A possible solution to this problem may be that cosmic rays undergo most of their interactions near their sources due to the efficient trapping in the far upstream of supernova shocks by self-generated, cosmic ray-driven turbulence., Comment: 31 pages, accepted for publication in ApJ (abstract abridged for arXiv)

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares. We present here self-consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the NRL flux tube code. We calculated the collisional heating rate from the particle transport code, which is more accurate than those based on approximate analytical solutions. We used a realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a nonthermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (vs. chromosphere), a larger downward conductive flux, and thus a stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced X-ray emission., Comment: Accepted by ApJ (2009 June 12), 15 page, 12 figures

It is well-recognized that the presence of magnetic fields will lead to anisotropic energy cascade and dissipation of astrophysical turbulence. With the diffusion approximation and linear dissipation rates, we study the cascade and damping of Alfv\'en-cyclotron fluctuations in solar plasmas numerically. For an isotropic case the steady-state turbulence spectra are nearly isotropic in the inertial range and can be fitted by a single power-law function with a spectral index of -3/2, similar to the Iroshnikov-Kraichnan phenomenology. Beyond the MHD regime the kinetic effects make the spectrum softer at higher wavenumbers. In the dissipation range the turbulence spectrum cuts off at the wavenumber, where the damping rate becomes comparable to the cascade rate, and the cutoff wavenumber changes with the wave propagation direction. The angle averaged turbulence spectrum of the isotropic model resembles a broken power-law. Taking into account the Doppler effects, the model naturally reproduces the broken power-law turbulence spectra observed in the solar wind and predicts that a higher break frequency always comes along with a softer dissipation range spectrum that may be caused by the increase of the turbulence intensity, the reciprocal of the plasma \beta, and/or the angle between the solar wind velocity and the mean magnetic field. These predictions can be tested by detailed comparisons with more accurate observations., Comment: 51 pages, 14 figures, to appear on ApJ

Solar Energetic Particles (SEPs) show a rich variety of spectra and relative abundances of many ionic species and their isotopes. A long standing puzzle has been the extreme enrichments of 3He ions. The most extreme enrichments are observed in low fluence, the so-called impulsive, events which are believed to be produced at the flare site in the solar corona with little scattering and acceleration during transport to the Earth. In two earlier papers (Liu et al. 2004 and 2006) we showed how such extreme enrichments can result in the model developed by Petrosian and Liu (2004), where ions are accelerated stochastically by plasma waves or turbulence. In this paper we address the relative distributions of the fluences of 3He and 4He ions presented by Ho et al. (2005) which show that while the distribution of 4He fluence like many other extensive characteristics of solar flare, is fairly broad, the 3He fluence is limited to a narrow range. Moreover, the ratio of the fluences shows a strong correlation with the 4He fluence. One of the predictions of our model was presence of steep variation of the fluence ratio with the level of turbulence or the rate of acceleration. We show here that this feature of the model can reproduce the observed distribution of the fluences with very few free parameters. The primary reason for the success of the model in both fronts is because fully ionized 3He ion, with its unique charge to mass ratio, can resonantly interact with more plasma modes and accelerate more readily than 4He. Essentially in most flares, all background 3He ions are accelerated to few MeV/nucleon range, while this happens for 4He ions only in very strong events. A much smaller fraction of 4He ions reach such energies in weaker events., Comment: 14 pages, 5 figures, submitted to ApJ

We discuss possible mechanisms of acceleration of particles in solar flares and show that turbulence plays an important role in all the mechanism. It is also argued that stochastic particle acceleration by turbulent plasma waves is the most likely mechanism for production of the high energy electrons and ions responsible for observed radiative signatures of solar flares and for solar energetic particle or SEPs, and that the predictions of this model agrees well with many past and recent high spectral and temporal observations of solar flares. It is shown that, in addition, the model explains many features of SEPs that accompany flares. In particular we show that it can successfully explain the observed extreme enhancement, relative to photospheric values, of $^3$He ions and the relative spectra of $^3$He and $^4$He. It has also the potential of explaining the relative abundances of most ions including the increasing enhancements of heavy ions with ion mass or mass-to-charge ratio., Comment: To appear in proceedings of ISSI Conference on Impulsive Solar Energetic Particle Events

We present a detailed imaging and spectroscopic study of the conjugate hard X-ray (HXR) footpoints (FPs) observed with RHESSI in the 2003 October 29 X10 flare. The double FPs first move toward and then away from each other, mainly parallel and perpendicular to the magnetic neutral line, respectively. The transition of these two phases of FP unshearing motions coincides with the direction reversal of the motion of the loop-top (LT) source, and with the minima of the estimated loop length and LT height. The FPs show temporal correlations between HXR flux, spectral index, and magnetic field strength. The HXR flux exponentially correlates with the magnetic field strength, which also anti-correlates with the spectral index before the second HXR peak's maximum, suggesting that particle acceleration sensitively depends on the magnetic field strength and/or reconnection rate. Asymmetries are observed between the FPs: on average, the eastern FP is 2.2 times brighter in HXR flux and 1.8 times weaker in magnetic field strength, and moves 2.8 times faster away from the neutral line than the western FP; the estimated coronal column density to the eastern FP from the LT source is 1.7 times smaller. The two FPs have marginally different spectral indexes. The eastern-to-western FP HXR flux ratio and magnetic field strength ratio are anti-correlated only before the second HXR peak's maximum. Neither magnetic mirroring nor column density alone can explain the totality of these observations, but their combination, together with other transport effects, might provide a full explanation. We have also developed novel techniques to remove particle contamination from HXR counts and to estimate effects of pulse pileup in imaging spectroscopy, which can be applied to other RHESSI flares in similar circumstances., Comment: 22 pages, 14 figures, 4 tables; ApJ 2009, in press

With the diffusion approximation, we study the cascade and damping of Alfv\'{e}n-cyclotron fluctuations in solar plasmas numerically. Motivated by wave-wave couplings and nonlinear effects, we test several forms of the diffusion tensor. For a general locally anisotropic and inhomogeneous diffusion tensor in the wave vector space, the turbulence spectrum in the inertial range can be fitted with power-laws with the power-law index varying with the wave propagation direction. For several locally isotropic but inhomogeneous diffusion coefficients, the steady-state turbulence spectra are nearly isotropic in the absence of damping and can be fitted by a single power-law function. However, the energy flux is strongly polarized due to the inhomogeneity that leads to an anisotropic cascade. Including the anisotropic thermal damping, the turbulence spectrum cuts off at the wave numbers, where the damping rates become comparable to the cascade rates. The combined anisotropic effects of cascade and damping make this cutoff wave number dependent on the wave propagation direction, and the propagation direction integrated turbulence spectrum resembles a broken power-law, which cuts off at the maximum of the cutoff wave numbers or the $^4$He cyclotron frequency. Taking into account the Doppler effects, the model can naturally reproduce the broken power-law wave spectra observed in the solar wind and predicts that a higher break frequency is aways accompanied with a greater spectral index change that may be caused by the increase of the Alfv\'{e}n Mach number, the reciprocal of the plasma beta, and/or the angle between the solar wind velocity and the mean magnetic field. These predictions can be tested by future observations.