Accelerating Posterior Inference from Pulsar Light Curves via Learned Latent Representations and Local Simulator-Guided Optimization

📅 2026-02-16
📈 Citations: 0
Influential: 0
📄 PDF
🤖 AI Summary
Traditional MCMC methods for posterior inference of pulsar light curves are computationally prohibitive. This work proposes an efficient alternative that first constructs a latent embedding space of simulated data using a pretrained masked U-Net, enabling rapid initial posterior estimates via nearest-neighbor retrieval. Subsequently, a forward physical simulator guides local hill-climbing optimization within this latent space to refine the likelihood. By uniquely integrating learned latent representations with simulator-informed local optimization, the method achieves a 120-fold acceleration—reducing inference time on real data from PSR J0030+0451 from 24 hours to 12 minutes—while preserving posterior accuracy consistent with full MCMC sampling.

Technology Category

Application Category

📝 Abstract
Posterior inference from pulsar observations in the form of light curves is commonly performed using Markov chain Monte Carlo methods, which are accurate but computationally expensive. We introduce a framework that accelerates posterior inference while maintaining accuracy by combining learned latent representations with local simulator-guided optimization. A masked U-Net is first pretrained to reconstruct complete light curves from partial observations and to produce informative latent embeddings. Given a query light curve, we identify similar simulated light curves from the simulation bank by measuring similarity in the learned embedding space produced by pretrained U-Net encoder, yielding an initial empirical approximation to the posterior over parameters. This initialization is then refined using a local optimization procedure using hill-climbing updates, guided by a forward simulator, progressively shifting the empirical posterior toward higher-likelihood parameter regions. Experiments on the observed light curve of PSR J0030+0451, captured by NASA's Neutron Star Interior Composition Explorer (NICER), show that our method closely matches posterior estimates obtained using traditional MCMC methods while achieving 120 times reduction in inference time (from 24 hours to 12 minutes), demonstrating the effectiveness of learned representations and simulator-guided optimization for accelerated posterior inference.
Problem

Research questions and friction points this paper is trying to address.

posterior inference
pulsar light curves
computational efficiency
Markov chain Monte Carlo
parameter estimation
Innovation

Methods, ideas, or system contributions that make the work stand out.

learned latent representations
simulator-guided optimization
posterior inference acceleration
masked U-Net
pulsar light curves
F
Farhana Taiyebah
Los Alamos National Laboratory, Los Alamos, NM, USA; Department of Scientific Computing, Florida State University, Tallahassee, FL, USA
A
Abu Bucker Siddik
Los Alamos National Laboratory, Los Alamos, NM, USA
I
Indronil Bhattacharjee
Los Alamos National Laboratory, Los Alamos, NM, USA; Department of Computer Science, New Mexico State University, Las Cruces, NM, USA
Diane Oyen
Diane Oyen
Los Alamos National Laboratory
Machine Learning
S
Soumi De
Los Alamos National Laboratory, Los Alamos, NM, USA
G
Greg Olmschenk
NASA Goddard Space Flight Center, Greenbelt, MD, USA
C
Constantinos Kalapotharakos
NASA Goddard Space Flight Center, Greenbelt, MD, USA