Related papers: Fast Multipole Method for Gravitational Lensing. A…
We propose a novel means of directly measuring cosmological distances using scintillated microlensing of fast radio bursts (FRBs). In standard strong lensing measurements of cosmic expansion, the main source of systematic uncertainty lies…
Vortex element methods are often used to efficiently simulate incompressible flows using Lagrangian techniques. Use of the FMM (Fast Multipole Method) allows considerable speed up of both velocity evaluation and vorticity evolution terms in…
Magnetic resonance fingerprinting (MRF) is a new technique for simultaneously quantifying multiple MR parameters using one temporally resolved MR scan. But its brute-force dictionary generating and searching (DGS) process causes a huge disk…
Recently, there have been two landmark discoveries of gravitationally lensed supernovae: the first multiply-imaged SN, "Refsdal", and the first Type Ia SN resolved into multiple images, SN iPTF16geu. Fitting the multiple light curves of…
It is believed that rays bent inside and outside photon spheres could affect partially the black hole shadow images by the Event Horizon Telescope and the rays near photon spheres would be detected by near-future space observations. The…
Weak gravitational lensing measurements based on photometry are limited by shape noise, the variance in the unknown unlensed orientations of the source galaxies. If the source is a disk galaxy with a well-ordered velocity field, however,…
In this paper, we propose a fast multipole method (FMM) for 3-D linearized Poisson-Boltzmann (PB) equation in layered media. The main framework of the algorithm is analogous to the FMM for Helmholtz and Laplace equation in layered media…
Gravitational lensing is one of the most powerful probes of dark matter, yet creating high-fidelity lensed images at scale remains a bottleneck. Existing tools rely on ray-tracing or forward-modeling pipelines that, while precise, are…
In this work, the fast-convolving reproducing kernel particle method (FC-RKPM) is introduced. This method is hundreds to millions of times faster than the traditional RKPM for 3D meshfree simulations. In this approach, the meshfree…
This paper uses the Schwarzschild metric to derive an effective refractive index and acceleration vector that account for relativistic deflection of light rays, in an otherwise classical kinematic framework. The new refractive index and the…
To take full advantage of the unprecedented power of upcoming weak lensing surveys, understanding the noise, such as cosmic variance and geometry/mask effects, is as important as understanding the signal itself. Accurately quantifying the…
We propose a numerical interferometry method for identification of optical multiply-scattering systems when only intensity can be measured. Our method simplifies the calibration of optical transmission matrices from a quadratic to a linear…
Context: The determination of the stellar initial mass function (IMF) of massive galaxies is one of the open problems in cosmology. Strong gravitational lensing is one of the few methods that allow us to constrain the IMF outside of the…
We present an implementation of the fast multipole method for computing coulombic electrostatic and polarization forces from polarizable force-fields based on induced point dipole moments. We demonstrate the expected $O(N)$ scaling of that…
In strong gravitational lens systems, the light bending is usually dominated by one main galaxy, but may be affected by other mass along the line of sight (LOS). Shear and convergence can be used to approximate the contributions from less…
We propose a new mass-mapping algorithm, specifically designed to recover small-scale information from a combination of gravitational shear and flexion. Including flexion allows us to supplement the shear on small scales in order to…
Purpose: Previously-developed Voxel Spread Function (VSF) method (Yablonskiy, et al, MRM, 2013;70:1283) provides means to correct artifacts induced by macroscopic magnetic field inhomogeneities in the images obtained by…
This paper presents an accelerated quadrature scheme for the evaluation of layer potentials in three dimensions. Our scheme combines a generic, high order quadrature method for singular kernels called Quadrature by Expansion (QBX) with a…
The Initial Mass Function (IMF) for massive galaxies can be constrained by combining stellar dynamics with strong gravitational lensing. However, this method is limited by degeneracies between the density profile of dark matter and the…
In this paper, a fast multipole method (FMM) is proposed for 3-D Laplace equation in layered media. The potential due to charges embedded in layered media is decomposed into a free space component and four types of reaction field…