Bayesian Estimation of Experimental Parameters in Stochastic Inertial Systems: Theory, Simulations, and Experiments with Objects Levitated in Vacuum

TL;DR

This paper presents a Bayesian inference method for characterizing stochastic inertial systems using noisy position data, validated through simulations and experiments with levitated particles, enabling precise force and environment estimation.

Contribution

It introduces a novel Bayesian inference approach that operates solely on position time series, accurately estimating forces and environmental parameters without simplifying assumptions.

Findings

Method accurately characterizes nonlinear and non-conservative forces.

Requires significantly shorter trajectories than existing methods.

Provides guidance on optimal sampling frequency and trajectory length.

Abstract

High-quality nanomechanical oscillators can sensitively probe force, mass, or displacement in experiments bridging the gap between the classical and quantum domain. Dynamics of these stochastic systems is inherently determined by the interplay between acting external forces, viscous dissipation, and random driving by the thermal environment. The importance of inertia then dictates that both position and momentum must, in principle, be known to fully describe the system, which makes its quantitative experimental characterization rather challenging. We introduce a general method of Bayesian inference of the force field and environmental parameters in stochastic inertial systems that operates solely on the time series of recorded noisy positions of the system. The method is first validated on simulated trajectories of model stochastic harmonic and anharmonic oscillators with damping. Subsequently, the method is applied to experimental trajectories of particles levitating in tailored optical fields and used to characterize the dynamics of particle motion in a nonlinear Duffing potential, a static or time-dependent double-well potential, and a non-conservative force field. The presented inference procedure does not make any simplifying assumptions about the nature or symmetry of the acting force field and provides robust results with trajectories two orders of magnitude shorter than those typically required by alternative inference schemes. In addition to being a powerful tool for quantitative data analysis, it can also guide experimentalists in choosing appropriate sampling frequency (at least 20 measured points per single characteristic period) and length of the measured trajectories (at least 10 periods) to estimate the force field and environmental characteristics with a desired accuracy and precision.