Optical trap Langevin calibration ($\zeta \dot{x}=-kx+\xi (t),f_c=k/(2\pi \zeta )$)

Add the random thermal force to the deterministic optical trap relaxation time model. In one stationary coordinate of a harmonic optical trap, the overdamped Langevin equation isHere $\xi$ is Gaussian white noise understood through its integrals, $\delta$ is the Dirac delta function, and $\zeta$ is the linear friction coefficient. The covariance coefficient has units ${\mathrm{N}}^2\mathrm{s}$; multiplying by $\delta (t-t^{\prime })$, of units ${\mathrm{s}}^{-1}$, gives squared force. This fluctuation-dissipation relation for a Langevin particle makes the stationary variance ${\sigma }^2=k_{B}T/k$, consistent with the equipartition theorem. The solution is an Ornstein-Uhlenbeck process with autocovariance functionIndeed, at positive lag the noise has zero correlation with the initial position, so $\zeta C_x^{\prime }=-k{C}_x$; its initial value is ${\sigma }^2$. This supplies both a decay time and an amplitude, making a stronger calibration test than a position histogram alone.

Specify the Fourier transform convention before using a power spectral density. For frequency $f$ in hertz, define $S_x^{(2)}(f)={\int }_{-\infty }^{\infty }{C}_x(t)e^{-2\pi ift}dt$ and the one-sided density $S_x^{(1)}(f)=2S_x^{(2)}(f)$ for $f>0$. Integration of the exponential givesThe position power spectral density has units ${\mathrm{m}}^2/Hz$; its integral is a variance, not a standard deviation. At $f=f_c$ it is half its zero-frequency plateau. Angular frequency conventions require a corresponding $2\pi$ conversion in both axis and integration measure. If $x=\beta (V-V_0)$ with position sensitivity $\beta$ in $m/V$, then $S_V^{(1)}=S_x^{(1)}/{\beta }^2$ in ${\mathrm{V}}^2/Hz$ before independent detector noise is added.

For an original invented calibration example, choose $a=0.500\mu \mathrm{m}$, dynamic viscosity $\eta =1.00\times 10^{-3}\mathrm{Pa}\mathrm{s}$, $T=300\mathrm{K}$, $k=0.0200pN/nm=2.00\times 10^{-5}N/m$, and $\beta =100nm/V$. The unbounded no-slip Stokes drag law givesThe corresponding voltage standard deviation is $0.14391\mathrm{V}$, and the position power spectral density plateau is $0.390369\mathrm{n}{\mathrm{m}}^2/Hz$, or $3.90369\times 10^{-5}{\mathrm{V}}^2/Hz$ in the detector. Recovering both $k=k_{B}T/{\sigma }^2$ and $k=2\pi \zeta f_c$ is an internal model check; actual agreement also tests the assumed temperature, drag and detector response.

The generator uses the exact stationary update $x_{j+1}=\rho x_j+\sigma \sqrt{1-{\rho }^2}{Z}_j$, with $\rho =e^{-\Delta t/\tau }$ and independent standard normal random variables $Z_j$. For a demeaned block of $N$ samples and its unnormalized discrete Fourier transform $X_j$, the two-sided periodogram is $\Delta t|X_j{|}^2/N$; doubling positive bins except the zero and Nyquist endpoints gives a one-sided density. Summing those bins times $\Delta f=1/(N\Delta t)$ returns the block mean square by discrete Parseval identity. Independent detector noise, finite exposure and aliasing would alter a measured spectrum and must enter a fit before interpreting $f_c$.