-
An Efficient Method for Simulating Light Curves of Cosmological Microlensing and Caustic Crossing Events
Authors:
Ashish Kumar Meena,
Ofir Arad,
Adi Zitrin
Abstract:
A new window to observing individual stars and other small sources at cosmological distances was opened recently, with the detection of several caustic-crossing events in galaxy cluster fields. Many more such events are expected soon from dedicated campaigns with the \emph{Hubble Space Telescope} and from the \emph{James Webb Space Telescope}. These events can teach us not only about the lensed so…
▽ More
A new window to observing individual stars and other small sources at cosmological distances was opened recently, with the detection of several caustic-crossing events in galaxy cluster fields. Many more such events are expected soon from dedicated campaigns with the \emph{Hubble Space Telescope} and from the \emph{James Webb Space Telescope}. These events can teach us not only about the lensed sources themselves, such as individual high-redshift stars, star clusters, or accretion disks, but through their light-curves they also hold information about the point-mass function of the lens and thus, potentially, the composition of dark matter. We present here a simple method for simulating light curves of such events, i.e., the change in apparent magnitude of the source as it sweeps over the net of caustics generated by microlenses embedded around the critical region of the lens. The method is recursive and so any reasonably sized small source can be accommodated, down to sub-solar scales, in principle. We compare the method, which we dub \emph{Adaptive Boundary Method}, with other common methods such as simple inverse ray shooting, and demonstrate that it is significantly more efficient and accurate in the small-source and high-magnification regime of interest. A \textsc{python} version of the code is made publicly available in an open-source fashion for simulating future events.
△ Less
Submitted 12 August, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Maximally-hard spectra from diffusive shock-acceleration in pulsar-wind nebulae
Authors:
Ofir Arad,
Assaf Lavi,
Uri Keshet
Abstract:
The processes leading to the exceptionally hard radio spectra of pulsar-wind nebulae (PWNe) are not yet understood. Radio photon spectral indices among $29$ PWNe from the literature show an approximately normal, $α=0.2\pm0.2$ distribution. We present $\sim 3σ$ evidence for a distinct sub-population of PWNe, with a hard spectrum $α=0.01\pm0.06$ near the termination shock and significantly softer el…
▽ More
The processes leading to the exceptionally hard radio spectra of pulsar-wind nebulae (PWNe) are not yet understood. Radio photon spectral indices among $29$ PWNe from the literature show an approximately normal, $α=0.2\pm0.2$ distribution. We present $\sim 3σ$ evidence for a distinct sub-population of PWNe, with a hard spectrum $α=0.01\pm0.06$ near the termination shock and significantly softer elsewhere, possibly due to a recent evacuation of the shock surroundings. Such spectra, especially in the hard sub-population, suggest a Fermi process, such as diffusive shock acceleration, at its extreme, $α=0$ limit. We show that this limit is approached for sufficiently anisotropic small-angle scattering, enhanced on either side of the shock for particles approaching the shock front. In the upstream, the spectral hardening is mostly associated with an enhanced energy gain, possibly driven by the same beamed particles crossing the shock. Downstream, the main effect is a diminished escape probability, but this lowers the acceleration efficiency to $\lesssim25\%$ for $α=0.3$ and $\lesssim1\%$ for $α=0.03$.
△ Less
Submitted 1 December, 2020;
originally announced December 2020.
-
Diffusive shock acceleration in $N$ dimensions
Authors:
Assaf Lavi,
Ofir Arad,
Yotam Nagar,
Uri Keshet
Abstract:
Collisionless shocks are often studied in two spatial dimensions (2D), to gain insights into the 3D case. We analyze diffusive shock acceleration for an arbitrary number $N\in\mathbb{N}$ of dimensions. For a non-relativistic shock of compression ratio $\mathcal{R}$, the spectral index of the accelerated particles is $s_E=1+N/(\mathcal{R}-1)$; this curiously yields, for any $N$, the familiar…
▽ More
Collisionless shocks are often studied in two spatial dimensions (2D), to gain insights into the 3D case. We analyze diffusive shock acceleration for an arbitrary number $N\in\mathbb{N}$ of dimensions. For a non-relativistic shock of compression ratio $\mathcal{R}$, the spectral index of the accelerated particles is $s_E=1+N/(\mathcal{R}-1)$; this curiously yields, for any $N$, the familiar $s_E=2$ (i.e., equal energy per logarithmic particle energy bin) for a strong shock in a mono-atomic gas. A precise relation between $s_E$ and the anisotropy along an arbitrary relativistic shock is derived, and is used to obtain an analytic expression for $s_E$ in the case of isotropic angular diffusion, affirming an analogous result in 3D. In particular, this approach yields $s_E = (1+\sqrt{13})/2 \simeq 2.30$ in the ultra-relativistic shock limit for $N=2$, and $s_E(N\to\infty)=2$ for any strong shock. The angular eigenfunctions of the isotropic-diffusion transport equation reduce in 2D to elliptic cosine functions, providing a rigorous solution to the problem; the first function upstream already yields a remarkably accurate approximation. We show how these and additional results can be used to promote the study of shocks in 3D.
△ Less
Submitted 25 February, 2020;
originally announced February 2020.
-
Diffusive shock-acceleration: breakdown of spatial diffusion and isotropy
Authors:
Uri Keshet,
Ofir Arad,
Yuri Lyubarski
Abstract:
We point out that particles accelerated in a non-relativistic shock of compression ratio $r$ attain the standard, $p=(r+2)/(r-1)$ spectral index only under certain conditions. Previous derivations of the spectrum, based on the approximations of spatial diffusion or negligible anisotropy, are shown to be unjustified for a general scattering function. We explain and demonstrate numerically that in c…
▽ More
We point out that particles accelerated in a non-relativistic shock of compression ratio $r$ attain the standard, $p=(r+2)/(r-1)$ spectral index only under certain conditions. Previous derivations of the spectrum, based on the approximations of spatial diffusion or negligible anisotropy, are shown to be unjustified for a general scattering function. We explain and demonstrate numerically that in contrast to previous claims, $p$ can substantially deviate from the standard result for anisotropic scattering. We prove analytically that the standard approach is nevertheless valid in the limit of an isotropic medium. Additional spectral modifications, for example by motions of scattering modes at intermediate optical depths from the shock, are discussed.
△ Less
Submitted 17 October, 2019;
originally announced October 2019.