Source author record

Mario A. Riquelme

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

ResearcherUnclaimed source record

astro-ph.HE astro-ph physics.plasm-ph

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2012arXiv

Local 2D Particle-in-cell simulations of the collisionless MRI

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the rate of Coulomb collisions between particles is very small, making the disk evolve essentially as a collisionless plasma. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell (PIC) plasma simulations. In this initial study we focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field until the Alfvén velocity, v_A, is comparable to the speed of light, c (independent of the initial value of v_A/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the perpendicular pressure being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission -- from the radio to the gamma-rays -- of systems such as Sgr A*.

preprint2011arXiv

Electron Injection by Whistler Waves in Non-relativistic Shocks

Electron acceleration to non-thermal, ultra-relativistic energies (~ 10-100 TeV) is revealed by radio and X-ray observations of shocks in young supernova remnants (SNRs). The diffusive shock acceleration (DSA) mechanism is usually invoked to explain this acceleration, but the way in which electrons are initially energized or 'injected' into this acceleration process starting from thermal energies is an unresolved problem. In this paper we study the initial acceleration of electrons in non-relativistic shocks from first principles, using two- and three-dimensional particle-in-cell (PIC) plasma simulations. We systematically explore the space of shock parameters (the Alfvénic Mach number, M_A, the shock velocity, v_{sh}, the angle between the upstream magnetic field and the shock normal, theta_{Bn}, and the ion to electron mass ratio, m_i/m_e). We find that significant non-thermal acceleration occurs due to the growth of oblique whistler waves in the foot of quasi-perpendicular shocks. The obtained electron energy distributions show power law tails with spectral indices up to alpha ~ 3-4. The maximum energies of the accelerated particles are consistent with the electron Larmor radii being comparable to that of the ions, indicating potential injection into the subsequent DSA process. This injection mechanism, however, requires the shock waves to have fairly low Alfénic Mach numbers, M_A <~ 20, which is consistent with the theoretical conditions for the growth of whistler waves in the shock foot (M_A <~ (m_i/m_e)^{1/2}). Thus, if the whistler mechanism is the only robust electron injection process at work in SNR shocks, then SNRs that display non-thermal emission must have significantly amplified upstream magnetic fields. Such field amplification is likely achieved by the escaping cosmic rays, so electron and proton acceleration in SNR shocks must be interconnected.

preprint2009arXiv

Magnetic Amplification by Magnetized Cosmic Rays in SNR Shocks

(Abridged) X-ray observations of synchrotron rims in supernova remnant (SNR) shocks show evidence of strong magnetic field amplification (a factor of ~100 between the upstream and downstream medium). This amplification may be due to plasma instabilities driven by shock-accelerated cosmic rays (CRs). One candidate is the cosmic ray current-driven (CRCD) instability (Bell 2004), caused by the electric current of large Larmor radii CRs propagating parallel to the upstream magnetic field. Particle-in-cell (PIC) simulations have shown that the back-reaction of the amplified field on CRs would limit the amplification factor of this instability to less than ~10 in galactic SNRs. In this paper, we study the possibility of further amplification driven near shocks by "magnetized" CRs, whose Larmor radii are smaller than the length scale of the field that was previously amplified by the CRCD instability. We find that additional amplification can occur due to a new instability, driven by the CR current perpendicular to the field, which we term the "perpendicular current-driven instability" (PCDI). We derive the growth rate of this instability, and, using PIC simulations, study its non-linear evolution and saturation. We find that PCDI increases the amplification of the field (amplification factor up to ~45, not including the shock compression) and discuss its observational signatures. Our results strengthen the idea of CRs driving a significant part of the magnetic field amplification observed in SNR shocks.

preprint2008arXiv

Non-linear Study of Bell's Cosmic Ray Current-driven Instability

The cosmic ray current-driven (CRCD) instability, predicted by Bell (2004), consists of non-resonant, growing plasma waves driven by the electric current of cosmic rays (CRs) that stream along the magnetic field ahead of both relativistic and non-relativistic shocks. Combining an analytic, kinetic model with one-, two-, and three-dimensional particle-in-cell simulations, we confirm the existence of this instability in the kinetic regime and determine its saturation mechanisms. In the linear regime, we show that, if the background plasma is well magnetized, the CRCD waves grow exponentially at the rates and wavelengths predicted by the analytic dispersion relation. The magnetization condition implies that the growth rate of the instability is much smaller than the ion cyclotron frequency. As the instability becomes non-linear, significant turbulence forms in the plasma. This turbulence reduces the growth rate of the field and damps the shortest wavelength modes, making the dominant wavelength, λ_d, grow proportional to the square of the field. At constant CR current, we find that plasma acceleration along the motion of CRs saturates the instability at the magnetic field level such that v_A ~ v_{d,cr}, where v_A is the Alfven velocity in the amplified field, and v_{d,cr} is the drift velocity of CRs. The instability can also saturate earlier if CRs get strongly deflected by the amplified field, which happens when their Larmor radii get close to λ_d. We apply these results to the case of CRs in the upstream medium of supernova remnants. Considering only the most energetic CRs that escape from the shock, we obtain that the field amplification factor of ~10 can be reached. This confirms the CRCD instability as a potentially important component of magnetic amplification process in astrophysical shocks.

preprint2007arXiv

Separating the Weak Lensing and Kinetic SZ Effects from CMB Temperature Maps

A new generation of CMB experiments will soon make sensitive high resolution maps of the microwave sky. At angular scales less than $\sim$10 arcminutes, most CMB anisotropies are generated at z $< 1000$, rather than at the surface of last scattering. Therefore, these maps potentially contain an enormous amount of information about the evolution of structure. Whereas spectral information can distinguish the thermal Sunyaev-Zeldovich (tSZ) effect from other anisotropies, the spectral form of anisotropies generated by the gravitational lensing and the kinetic Sunyaev-Zeldovich (kSZ) effects are identical. While spectrally identical, the statistical properties of these effects are different. We introduce a new real-space statistic, $<θ(\hat{n})^3 θ(\hat{m})>_c$, and show that it is identically zero for weakly lensed primary anisotropies and, therefore, allows a direct measurement of the kSZ effect. Measuring this statistic can offer a new tool for studing the reionization epoch. Models with the same optical depth, but different reionization histories, can differ by more than a factor of 3 in the amplitude of the kSZ-generated non-Gaussian signal.