PhysicsarXiv

Ranked by formulaVerdict: Notable

Microcaustics: Fast GPU-Accelerated Dynamical Microlensing of Time-Dependent Extended Sources

Gravitational microlensing probes unresolved quasar accretion disc structure and compact objects in lens galaxies and is an important systematic in strong-lens time-delay cosmography.

By Fagin, Stone, Chan +7

Score████░░░░░░4.3

Key numbers

  • 780 times higher than that of

Caveats

  • Preprint; not yet peer reviewed.

VerdictWorth a reader's time today.

Read the originalPDFCode

Abstract

Gravitational microlensing probes unresolved quasar accretion disc structure and compact objects in lens galaxies and is an important systematic in strong-lens time-delay cosmography. Most analyses use static microlensing map simulations because evolving individual microlenses and caustics is computationally prohibitive. We present Microcaustics, a GPU-accelerated Python package built on PyTorch for dynamical microlensing of time-dependent extended sources in strong gravitational lenses. It combines a local-exact/far-field deflection solver, source-directed lens-plane scouting with temporal reuse, and biquadratic inverse polygon mapping (IPM). GPU batching and fused Triton rasterization accelerate the calculation, while binary caustic labels track source-centred crossings. Across three simulated quadruply lensed quasar systems, our package generates a 10-year, 147-frame multiband light curve in an amortized median of 0.68 s (0.88 s with labels) on a consumer GPU with a median light-curve RMSE of 0.35 mmag, while recovering every reference source-centred caustic crossing. At comparable map accuracy, the batched throughput is 780 times higher than that of a recent CUDA IPM implementation and extrapolates to a speedup of millions relative to classic GPU-D. We demonstrate that our IPM converges as N_ray⁻0.84, approaching the idealized N_ray⁻¹ convergence rate of piecewise-affine IPM and improving upon the N_ray⁻1/2 inverse ray-shooting rate. We also introduce a GPU-resident general relativistic and reverberating quasar source model and demonstrate our method with an evolving supernova source. These advances make fully dynamical stellar fields practical for Bayesian and simulation-based inference.

Joshua Fagin, Connor Stone, James Hung-Hsu Chan, Sophia Miskiewicz, Henry Best, Paras Sharma, Matthew O'Dowd, Laurence Perreault-Levasseur, Yashar Hezaveh, Simon Birrer

The editor's rubric

Heuristic review

DimensionLevelWeightWhat that level means
Leverage███░░ 318%A method or resource many groups across the field will adopt within a year.
Magnitude███░░ 320%Large gain: roughly 2x, or a clear new state of the art on a hard, unsaturated problem.
Evidence███░░ 322%Solid: multiple benchmarks or cohorts, ablations, fair baselines, released code or data.
Novelty██░░░ 222%A new combination of known ideas.
Trajectory██░░░ 210%Some room to improve with obvious engineering.
Stakes██░░░ 28%Benefits a professional community (practitioners, clinicians, engineers).

Editor’s rationale

Heuristic triage from title and abstract text only, not a reading of the paper. Cues found: method (we propose, we report); gains (x-fold); verification (error bars, code released); scale (efficient); stakes (global scale).

How the score was computed

rank-2026-10-07

Score████░░░░░░4.3

Score = 10 × (75% × adjusted merit / 10 + 15% × attention + 10% × freshness)

Merit
5.2 / 10
Weighted rubric, evidence-gated.
Adjusted merit
4.5 / 10
Shrunk toward the desk prior by editor confidence (42%).
Attention
0%
Citations, upvotes, points, mentions.
Freshness
91%
Half-life decay since publication.

No attention signals recorded yet.

The record

  • Reviewed by heuristic-v6 on Oct 15, 2026, 07:29 UTC. Paper type: method.
  • Categories: astro-ph.GA, astro-ph.IM
  • BRIEF, No.8 in the Physics edition of October 15, 2026.