banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Auto- and cross-power spectral analysis of dual trap optical tweezer experiments using Bayesian inference
Rent this article for


Image of FIG. 1.
FIG. 1.

(a) Schematic of the experimental setup described in Sec. II. (b) Zoom-in into the sample chamber: The notation for the coordinates of the beads (blue) relative to the trap centers (gray) is used throughout the manuscript.

Image of FIG. 2.
FIG. 2.

Cross sections of the flow field around an oscillating no-slip sphere in an incompressible fluid for various driving frequencies: (a) ω = 10−4 a , (b) ω = 10−2 a , (c) ω = 1/τ a . Streamlines indicate the direction of the velocity field and the color coding its magnitude, where F 0 is the driving force amplitude and μ0 = 1/(6πηa) the quasi-static sphere mobility. Spatial coordinates are given in units of the sphere radius a, frequencies in units of , the kinematic viscosity of the fluid being denoted by ν. Amplitudes of the driving force acting on and the velocity of the sphere are shown as red and blue arrows within the spheres. The figures above show the configuration of the flow-field, when the sphere velocity is maximal; animations showing the time evolution of these flow fields are available online.27

Image of FIG. 3.
FIG. 3.

Dependence of real and imaginary parts of the self-mobility of a sphere on the driving frequency ω. In the low-frequency limit, the quasi-stationary self-mobility μ0 ≡ 1/(6πηa) of a sphere of radius a in a fluid of shear viscosity η is recovered. On the lower axis, the angular frequency ω is rescaled by τ a a 2/ν, the upper frequency scale corresponds to the experimental case where τ a ≈ 0.26 μs (beads of radius a = 500 nm in water). The exact result by Stokes (dashed lines, Eq. (10)) and the approximate expressions (solid lines, Eq. (12)) nicely match over the entire range of frequencies resolved in our experiment (f ≲ 100 kHz).

Image of FIG. 4.
FIG. 4.

Frequency dependence of cross-mobilities parallel (γ = ∥) and perpendicular (γ = ⊥) to the inter-bead axis for different ratios of inter-bead separation R to bead radius a: (a) R/a = 20, (b) R/a = 3. Cross-mobilities for spheres of finite radius (solid lines, Eq. (13)) are compared to the Oseen result for point-like particles (dashed lines, Eq. (15)) and to quasi-static HI (horizontal arrows, Eq. (14)). Cross-mobilities are given in units of the quasi-static Stokes self-mobility μ0 ≡ 1/(6πηa) for a sphere of radius a in a fluid of viscosity η; the angular frequency is rescaled by the characteristic time scale τ a for vortex diffusion over the length scale of the bead's radius a, the upper frequency scale corresponds to experimental conditions (τ a ≈ 0.26 μs).

Image of FIG. 5.
FIG. 5.

Visualization of the signals in a typical dual trap optical tweezer experiment as described in Secs. V and VI A. (a) Signal processing as described by Eqs. (17) and (18); averaged periodograms obtained via DFT from the discretely sampled stochastic trajectories yield experimental estimates for the aliased PSDs . (b) Relation between auto- and cross-PSDs of the beads' thermal motion and aliased electrical PSDs according to Eqs. (22) and (23). Unknown parameters in the PSDs can be determined via a global maximum likelihood fit to as explained in Sec. VI B.

Image of FIG. 6.
FIG. 6.

Results of the three calibration measurements described in Sec. VII: (a) no trapped beads, (b) first trap occupied and second one empty, (c) vice versa. Averaged experimental periodograms are denoted by symbols, theoretical aliased PSDs , based on a global maximum likelihood fit to the data as described in the text, are shown as lines, where broken lines correspond to negative values of the PSDs; fit values are given in Table I. Results for the spatial coordinates parallel to the inter-trap axis (γ = ∥) are shown in the upper row, those for the perpendicular direction (γ = ⊥) in the lower one.

Image of FIG. 7.
FIG. 7.

Auto- and cross-PSDs of the fluctuations parallel to the inter-trap axis for various trap separations R: Symbols denote averaged experimental periodograms , lines the theoretical predictions for the aliased PSDs according to Eqs. (22) and (23): black lines include the full frequency dependence of HI (Eq. (13)), light blue lines correspond to instantaneous HI (Eq. (14)). Positive values of the averaged experimental periodograms are denoted by circles, negative ones by squares; similarly, solid and broken lines, respectively, denote positive and negative PSD values.

Image of FIG. 8.
FIG. 8.

Same as Fig. 7 but for the spatial component perpendicular (⊥) to the inter-trap axis.

Image of FIG. 9.
FIG. 9.

Influence of polarization crosstalk: Theoretical PSDs accounting for polarization crosstalk are shown as black lines (same as in Fig. 8), theoretical PSDs neglecting polarization crosstalk are denoted as purple lines; since auto-PSDs are only marginally affected by polarization crosstalk, black and purple lines for S 11 and S 22 overlap. Symbols denote experimental data (same as in Fig. 8).

Image of FIG. 10.
FIG. 10.

Attenuation resulting from various filtering sources as discussed in Appendix B: (a) parasitic filtering of the position sensing devices (Eq. (B1)) using the parameters from Table I for the ⊥ direction, (b) different electric low-pass filters (Eqs. (B3) and (B4)), (c) uniform averaging of signals over a time window of duration (Eq. (B5)), (d) time delay between the two recorded signals (Eq. (B7)). Positive values of are drawn as solid lines, negative ones are dashed.


Generic image for table
Table I.

Calibration results: Best fit parameters obtained from global maximum likelihood fits to the averaged periodograms in Fig. 6 as described in Sec. VII.

Generic image for table
Table II.

Overview of the notation used throughout the manuscript; sub- and super-scripts γ are omitted in the text, if the discussion equally applies to all spatial directions.


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Auto- and cross-power spectral analysis of dual trap optical tweezer experiments using Bayesian inference