Interrogating the ballistic regime in liquids with rotational optical tweezers
Feb. 13 , 2025SJ_Zhang

Abstract

Accessing the ballistic regime of single particles in liquids remains an experimental challenge that shrouds our understanding of particle–liquid interactions on exceedingly short time scales. We demonstrate rotational ballistic measurements of microspheres in liquids by observing the thermalization of the angular velocity. This study uses sensitive high-bandwidth polarization measurements from light scattered by orientation-locked birefringent probes trapped with rotational optical tweezers. The particle–liquid interactions in the ballistic regime are decoupled from the optical potential allowing direct studies of single-particle rotational dynamics in a previously inaccessible parameter space. This enabled us to observe and validate rotational hydrodynamic effects and perform ultrafast calibration-free viscometry using less than 50 ms of data. This methodology establishes a unique tool for studying microscopic rotational dynamics and probing highly dynamic microenvironments, including systems out of equilibrium.

1. INTRODUCTION

Complex systems with rapidly changing spatiotemporal dynamics underpin biological processes [15], fluid dynamics [68], fundamental thermodynamic laws [911], and, more generally, systems out of equilibrium [1216], yet there is a lack of sufficient measurement techniques to elucidate the physical and statistical properties governing them. Ballistic tweezers provide a unique tool to exploit particle–liquid interactions in microscopic systems to gain new insight into such properties [1719]. Sensing both the translational and rotational inertia is crucial to fully capture these interactions and interrogate systems out of equilibrium on exceedingly short time scales. The rotational ballistic motion in microscopic systems has remained a largely theoretical research topic and requires experimental techniques to corroborate results. In this work, we developed a technique that improves upon previous rotational optical tweezers methods and gains the capability to probe the complete transition from diffusive to ballistic rotational dynamics.

Optical tweezers (OTs) are a versatile tool in applied and fundamental physics as they can exert and sense positions, forces, and torques in nanoscale and microscale systems with great precision. The technique allows for direct manipulation of a system and to extract its physical characteristics through the use of optically trapped probes. This has led to notable success in studies of the mechanical properties governing complex biological processes. This includes DNA mechanics [2023], cellular adhesion forces [24,25], single molecule sensing [26,27], and cell viscoelasticity [2833], to cite a few.

The physical properties of the environment, including viscosity, can be investigated by analyzing the trajectory of an optically trapped probe, providing excellent spatially and temporally resolved microrheology measurements. Several well-established optical tweezers methods exist that rely on analyzing the stochastic dynamics of the probe [28,3439]. However, OTs can also access energy scales beyond stochastic dynamics and probe the transition into Newtonian dynamics as the temporal resolution improves [1719]. This opens a new measurement space to study particle–liquid interactions, specifically those on exceedingly short time scales down to milliseconds.

In liquids, micrometer-sized particles move due to thermal motion where kinetic energy is gained from random collisions with fluid molecules. The kinetic energy quickly dissipates giving rise to the apparent instantaneous displacements observed in Brownian motion. This is defined as the diffusive regime. However, the dissipation of the particle’s velocity can be observed with sufficient spatial and temporal resolution allowing access to the ballistic regime where the velocity can be resolved. Achieving this resolution is limited by the signal-to-noise ratio (SNR) of the detection system. However, improved SNR can be achieved using optically based methods.

Optical tweezers routinely interrogate the diffusive regime, and more recently, have successfully resolved the ballistic regime [1719] in translational geometry in liquids. Trapping experiments are performed most commonly in the diffusive regime and measure the position relaxation of the probe. This requires stiff optical traps and measurement times of the order of seconds. Improving the sensitivity to access the ballistic regime and instead measuring the velocity dissipation of the probe allows for unique studies of particle–liquid interactions on a spatiotemporal scale that was previously inaccessible. This is particularly important in the context of mobility and microrheometry, which differ between time scales [3537]. High-bandwidth measurements in liquids have resolved translational ballistic motion, enabling validation of hydrodynamic effects [17,18] and ultrafast viscometry measurements [19]. However, to fully characterize a system and understand the complete transition between continuum and stochastic mechanics, it is also crucial to elucidate the rotational dynamics. We demonstrate measurements in the rotational ballistic regime, where the dissipation of angular velocities is observed.

Rotation can be introduced to the system by transferring optical angular momentum from light to the probe in addition to the linear momentum transfer required for trapping, giving birth to rotational optical tweezers (ROTs) that are capable of exerting and precisely measuring torques [37,4045]. Typically, ROTs use a birefringent microsphere as a probe to transfer spin angular momentum from light. There is a measurable change in the polarization state of the trapping beam, which provides a useful and precise method to sense both torques and angular displacements. A rotating probe in ROTs does not move far from the center of the trap. This enables studies of DNA twisting mechanics [46], adhesion [47,48], and microrheometry [32,37,4042,49]. Recently, ROTs were used for intracellular viscometry [30,50,51], highlighting their unique applicability in confined spaces or near surfaces [52,53].

Rotational measurements based on polarization detection have inherent advantages over the typical deflection-based techniques employed in translational measurements. The latter responds to changes in the outgoing light and is most greatly affected by interactions at the beam focus. Unfortunately, deflected photons make up a small portion of the collected light, limiting the SNR in the detection and obscuring the ballistic regime. The signal can be improved by increasing the trapping power, using high-refractive-index particles, or by interferometric methods that filter out the trapping beam [18,19]. However, these methods have increased complexity compared to traditional measurements which limits their use.

The detection of angular position and torque using polarization can achieve a greater SNR as each photon carries information on the beam polarization state. The SNR is limited by the amount of spin angular momentum transferred to the particle (which is determined by the birefringence) and the noise in the detection apparatus. By introducing a balanced photodetector (BPD) to the measurement system—that cuts out common mode noise—we show that ROTs become sensitive enough to detect the transition into rotational ballistic dynamics. Improved noise characteristics in our experiments allow information to be obtained from analysis of the angular velocity power spectrum (AVPS) [19]—which has not been possible in our previous rotation experiments. This enhanced capability enables measurements of particle–liquid interactions on extremely short time scales.

2. ROTATIONAL BALLISTIC REGIME

Consider a spherical birefringent probe that is optically trapped in a purely viscous liquid. The rotational and translational components are decoupled due to the symmetry of the system; hence, we can examine the rotational dynamics alone. Vaterite microspheres are excellent high birefringence optical probes (calcium carbonate, positive uniaxial birefringent crystals with Δ?≈0.1) that can be reliably synthesized with good control over their morphology [54,55]. The interaction of a linearly polarized beam with their birefringence produces an optical torque that aligns the extraordinary axis of the probe with the polarization state of the incident beam. At small displacements, near equilibrium, this can be approximated as a harmonic potential [42]. As vaterite possesses positive uniaxial birefringence, circularly polarized light causes the trapped probe to rotate at a constant rate limited by the drag—with the Stokes drag coefficient, ?0=8???3, where ? is the viscosity of the fluid and ? is the radius of the probe.

Thermal torque impulses lead to stochastic motion that is opposed by the position-dependent optical torque and the velocity-dependent drag torque. The equation of motion that results is given by ??¨(?)=?optical+?thermal+?drag, where ? is the moment of inertia, and ? represents the torques in the system. As the reactionary optical torque aligns the probe to an equilibrium orientation, a characteristic time scale, ??, emerges associated with the angular position relaxation. Similarly, the drag torque from friction with the fluid opposes the angular velocities and a second characteristic time scale, ??, emerges relating to the dissipation of rotational inertia. At low Reynold’s numbers ????? which causes the inertial term to be neglected in most OTs experiments. These two parameters introduce corners into the angular velocity power spectrum at frequencies of ??=2?/?? and ??=2?/??, respectively, shown in Fig. 1. The former is commonly referred to as “the corner frequency.” Here, we have used an alternative naming convention linked to the physics they parameterize to discuss the two corners more easily.

figure: Fig. 1.

Fig. 1. The translational (cyan) and rotational (purple) velocity power spectra for a µ2µm radius microsphere normalized by their respective diffusion constant and with ??=10Hz, ??,trans.≈61kHz, and ??,trans.≈202kHz. The dashed lines indicate the models without hydrodynamic modifications [equivalent to ?(?)→?0 in Eq. (2)] for each case. The black line indicates a slope ∝?−3/2.

At short time scales nearing ??, the probe’s motion is complicated by the hydrodynamic interaction with the fluid and cannot be treated as linearly proportional to velocity. Instead, the drag depends on the past motion of the probe and results in a frequency-dependent drag coefficient, ?(?). The rotational drag for a rotating sphere is given by

\gamma (f) = {\gamma _0}\left({1 + \frac{2}{3}\frac{{(1 - i)(f/{f_v}{)^{3/2}} - i(f/{f_v})}}{{1 + 2(f/{f_v}{)^{1/2}} + 2(f/{f_v})}}} \right),

where we define a new characteristic frequency, ??, associated with the dissipation rate of the fluid’s inertia to parametrize the equation. The corresponding time scale, ??=2?/??, is of the same order of magnitude as ?? and can be defined in terms of ??. When ????, ?(?) approaches the Stokes drag constant for a rotating sphere at low Reynold’s numbers, ?0. The inertia of the fluid envelope surrounding the probe modifies the mass term to an effective mass that includes the displaced fluid such that ?=?particle+?fluid/2 [18,56]. Hence, ?=?particle+?fluid/2=(2/5)??2 becomes an effective moment of inertia. Compressibility effects are neglected because they dissipate significantly faster than the measurement rate [18,19].

The transition between diffusive and ballistic behavior is best observed in the angular velocity power spectrum, AVPS=(2??)|?~(?)|2, where ?~(?) is the unilateral Fourier transform of the angular displacement from equilibrium, ?(?). The velocity power spectra (VPS) for the translational and rotational cases are shown in Fig. 1, and can both be described by

{\rm VPS} = \frac{{2D\,{\rm Re}\,\left[{\gamma (f)/{\gamma _0}} \right]}}{{{{\left({{f_p}/f - f/{f_i} + {\rm Im}\left[{\gamma (f)/{\gamma _0}} \right]} \right)}^2} + {{\left({{\rm Re}\left[{\gamma (f)/{\gamma _0}} \right]} \right)}^2}}},

with appropriate substitutions of each parameter for the two cases. These parameters are defined in Table 1 for both translational and rotational motion. The translational ?(?) is not listed in this table but the function is well established [19,34]. The characteristic frequencies of inertia dissipation in the translational and rotational cases are related by ??,rot.=(10/3)??,trans.. The similar magnitudes imply that a similar level of sensitivity is required for resolving ballistic motion in both cases.

Tables Icon

Table 1. Parameters in the Translational and Rotational Velocity Power Spectra [Eqs. (1) and (2)], Defined in Terms of the Viscosity (?), Radius (?), Linear (?) and Angular (?) Trap Stiffnesses, Boltzmann’s Constant (??), and Temperature (?)

 

Figure 1 compares the linear and angular velocity power spectra and highlights the relative hydrodynamic contribution in each. The spectra have been normalized by their respective diffusion constant and the position relaxation frequency has been set to ??=10Hz. There is a negligible difference in the low-frequency limit (diffusive regime) as both cases have harmonic potentials and the drag is linearly proportional to the velocity. In the high-frequency limit (ballistic regime), the spectra decay proportionally to ?−3/2 due to the modified drag term. Vortices are formed in response to the motion of the probe and decay at a rate related to ??. Although ?? is similar in magnitude for both cases, linear displacements generate vortices in the surrounding fluid, whereas angular displacements generate a line of vortices along the rotation axis originating within the probe. The vortices in the fluid propagate further and take longer to decay, leading to greater hydrodynamic effects observed in the translational case between ?? and ??. Ballistic tweezers have experimentally verified the translational drag model for micrometer-sized particles in liquids [18,19], whereas the rotational analog with its significantly smaller deviation has yet to be investigated. The rotational ballistic measurements described in this study provide a viable technique for experimentally verifying these dynamics.

3. MATERIALS AND METHODS

A. Rotational Optical Tweezers

The optical trapping component of the setup has been described previously [30,37] and a simplified schematic is depicted in Fig. 2. Briefly, trapping is achieved using a 1064 nm NdYAG laser with powers at the trap between 20 and 200 mW. A half-wave plate controls the equilibrium orientation of the vaterite when the trap is linearly polarized. A removable quarter-wave plate enables switching between linearly and circularly polarized traps. The degree of circular polarization of the forward scattered light is measured to calibrate the system and determine the optical torque [30].

figure: Fig. 2.

Fig. 2. Simplified schematic of the rotational optical tweezers setup, adapted from [30]. Photodetectors D1 and D2 (PDA36A-EC, Thorlabs) measure the initial power of each beam. The beams are expanded and a series of wave plates controls the incident polarization before entering the high NA objective to form the optical trap in the sample. The outgoing light is collected by a condenser where D2 and D3 measure the circular polarization states of the IR beam to determine the optical torque and the balanced photodetector (BPD) measures the linear polarization state of the He–Ne beam to determine angular positions. The BPD outputs the corresponding voltage (+ and - output channels) and an amplified voltage signal of the difference between them (diff.). The right illustrates the force and torque a vaterite experiences in a linearly polarized trap where the beam axis is coming out of the page.

A low-powered (<5mW) and weakly focused helium–neon laser (HL210LB, Thorlabs) measures angular displacements of the vaterite. The He–Ne beam provides an independent measurement basis with negligible effects to the optical trap and has a superior noise spectrum compared to our trapping beam such that it is capable of resolving a signal at higher frequencies. The outgoing linear polarization states are measured using a balanced photodetector (HBPR-100M-60K-SI, Femto). The signal from the photodetector is captured using an oscilloscope (DHO4240, Rigol) at 1 GHz to minimize the noise contribution from the device and ensure a sufficient acquisition rate to avoid aliasing in the measurement bandwidth.

Introducing the balanced photodetector enables access to higher frequency measurements previously unobtainable in ROTs experiments. As opposed to typical photodetectors that output a voltage proportional to the intensity of the light, the BPD measures the difference in intensity between two sources. This is ideal for measuring the signal from the orientation-locked vaterite which causes small changes to the polarization as it fluctuates under rotational thermal motion. In our setup, a polarizing beam splitter separates the light into orthogonal linear polarization states that are each incident on the BPD. Measuring the difference between intensities suppresses the common-mode noise up to 40 dB and effectively gains an order of magnitude in temporal resolution compared to previous measurements. The maximum sensitivity occurs when the two states have equal powers, and this is achieved by orientating the vaterite such that the polarization states are equal. Shot noise will vary slightly between the two detectors, which places an upper limit on the effectiveness. However, the SNR in these experiments was limited by laser noise and vaterite morphology.

B. Sample Preparation

Vaterite microspheres were synthesized according to the procedure outlined by [54]. Vaterite is a positive uniaxial birefringent crystal with a density of 2540kg/m3 which has an approximately hyperbolic internal structure when synthesized as a sphere [57] with refractive indices of ??=1.55 and ??=1.65. Briefly, 1 mol of aqueous solutions of CaCl2 (5 mL), MgSO4 (2 mL), and K2CO3 (1.5 mL) were mixed for 3 min followed by the addition of µ50µL AGFA Agepon to stabilize the reaction and prevent further crystal growth. The vaterites were washed three times with ethanol. There are a variety of methods for controllable synthesis of vaterites with size ranges varying from µ0.5µm to µ25µm in diameter [55,5861]; however, much of this research does not consider the birefringence of vaterites, but rather focuses on their capability for drug loading and release. Our synthesis method enabled us to produce vaterite microspheres with a highly reproducible size range and birefringence, which are the key parameters for ROTs experiments. The vaterite microspheres produced for this work had radii ranging from µ0.7µm to µ2.5µm. The vaterites were slightly elliptical, where the ratio of the short axis to the long axis was 0.95±0.3.

µ50µL of vaterite suspension was added to 1 mL of the medium of interest, which for the measurements included in this work, was water and water–glycerol solutions. The solution was sonicated to separate any vaterite microspheres that clumped together. Care must be taken to ensure the vaterite concentration of the system is not too dense to avoid trapping multiple particles.

4. RESULTS AND DISCUSSION

A. Accessing the Ballistic Regime

The AVPS in Fig. 3A demonstrates the full transition between the rotational diffusive and ballistic regime. It is experimentally observed as evidenced by the two corners in the curve becoming resolvable at the characteristic frequencies of ?? and ??. The data is fitted to the AVPS model in Eq. (2) with fitting parameters for ?, ??, and ??, as well as an additional noise term characterized empirically from the signal measured from an empty trap. The magnitude of the noise floor varies due to acquisition settings and the birefringence of the vaterite. For this wideband measurement that captures the full transition, the signal becomes dominated by noise at frequencies above ??. Our typical high-frequency measurements can resolve the power spectra beyond 1 MHz as they only determine ?? and not ??, as discussed in the following section.

figure: Fig. 3.

Fig. 3. A, wideband AVPS measurement for a µ2.3µm sized probe measured at 10 MHz for 1 s showing the transition from the diffusive to the ballistic regime. The data was fitted to Eq. (2) (black line) with an empirically determined noise floor and used to determine ?? and ?? (blue and red circles on the curve). B, a 2 ms snapshot of the measured angular position, ?(?), and angular velocity, ?˙(?), time traces for µ1.4µm radius probe measured over 50 ms with the data low-pass filtered at 10 MHz. Further windowing of ?˙(?) with µ5µs is shown in the black line. C, a histogram of ?˙ captured over the full 50 ms and fitted with a normal distribution (black line).

Importantly, the spectrum plateaus between ?? and ?? as predicted by the hydrodynamically correct model shown in Fig. 1. This indicates that the signal from our rotational measurement is decoupled from the translational dynamics as translational hydrodynamic contributions would cause a deviation at lower frequencies.

The high-frequency measurements aimed at determining ?? and were measured at 100 MHz for 50 ms. The short measurement time prevents the diffusive regime and ?? from being resolved but provides sufficient time to capture the rotational ballistic dynamics and measure angular velocities. Figure 3B shows a 2 ms time trace of the observed angular position and velocity of the probe. The increased sensitivity of our detection systems means ?(?) no longer appears as a stochastic signal and instead becomes differentiable enabling the angular velocity to be determined, ?˙(?)≈Δ?/Δ?. The measured signal is low-pass filtered at 10 MHz to remove the noise-dominated spectrum, which reveals that the slope ?(?) is correlated with ?˙(?). The angular velocity thermalizes to reveal a Maxwell–Boltzmann distribution demonstrated in Fig. 3C. The measured rms?rms=(2.744±0.004)×103rad2/s2 reasonably agrees with the expected modified-equipartition theorem with rms?rms≈2.88×103rad2/s2. Observing the correlation in the angular velocity trace and its distribution is an important step in verifying a successful ballistic measurement [17].

Vaterites are useful and convenient optical probes frequently used in ROTs experiments as they can be easily synthesized, are close to spherical, and have strong birefringence. Birefringence is important for generating optical torque and achieving high SNR because it causes the spin angular momentum transfer from the light and thus affects the sensitivity of the polarization-dependent measurements of angular position. We observed a variation in the birefringence between otherwise equal vaterite probes, which changes the resolving power of the angular displacements. Vaterite has low absorbance at 633 nm and 1064 nm, which is crucial to minimize localized heating around the probe which would otherwise change the local fluid properties.

Probes below µ0.5µm radius can become less spherical and slightly more difficult to trap, which sets a reasonable upper limit for the required sensitivity of our technique (??=1.74MHz in water at 20°C). Vaterites in ROTs experiments are more commonly sized 1–3 µm in radius with corresponding inertia dissipation rates of 0.44–0.04 MHz. The spherical geometry simplifies the equations of motion and helps to decouple the translational and rotational degrees of motion. It should be noted that spherical geometry is less important for measurements aimed at sensing optical torque as the rotational dynamics, including the hydrodynamic effects, can be neglected. Variation in shape could cause coupling of translational and rotational dynamics. However, the aspect ratios of particles used in our experiments were close enough to unity so that coupling was small. The flat plateau in the AVPS curve in Fig. 3A implies negligible translational–rotational coupling. This is strong evidence that our detection technique is sensitive only to rotation.

There are few alternatives to vaterite as a probe in our system; however, the possibility of engineering probes with controllable birefringence and size remains a tantalizing opportunity to improve the sensitivity and precision of rotational ballistic measurements. Vaterite microspheres can also be coated and functionalized [54]. This improves their stability in more environments and can be very beneficial when applied to a large variety of biological systems.

B. Ultrafast Measurements

Ballistic measurements were performed by recording the angular trajectory with 50 ms measurement times captured at 100 MHz. The AVPS data was windowed by taking the mean of consecutive 1 kHz bins (also referred to as “blocking” [34]) before being truncated at high frequencies once the signal-to-noise ratio (SNR)≈1 and information could no longer be retrieved. A frequency-dependent function, ?(?), describes the noise floor and was determined from measurements without a particle in the optical trap and scaled to the total optical power on the detector. Non-linear least-squares fitting of the data to the VPS model [Eq. (2)] with the noise floor was performed. Under these conditions the fitting model becomes

\begin{split}{{\rm AVPS}_{{\rm fit}}}(f) & = \frac{{2D\,{\rm Re}\left[{\gamma (f)/{\gamma _0}} \right]}}{{{{\left({{f_p}/f - f/{f_i} + {\rm Im}\left[{\gamma (f)/{\gamma _0}} \right]} \right)}^2} + {{\left({{\rm Re}\left[{\gamma (f)/{\gamma _0}} \right]} \right)}^2}}} \\ & \quad + N(f),\end{split}

where only ? and ?? were used as free parameters because ?(?)/?0 can be described in terms of ?? and the total measurement time was too short to resolve ??, which was instead determined from longer-time calibration measurements.

Figure 4A shows a typical experimental AVPS that has been fitted to Eq. (3) (blue) and compared with the theoretical prediction calculated from Eq. (2) (red). The deviation from the expected model in red was caused by the signal approaching the noise floor (∝?2), which became dominant beyond 1 MHz. The shape of the noise spectrum, ?(?), is easily characterized; however, its magnitude and hence its prevalence depends on the SNR and the optical power on the detector. The SNR changes slightly between probes due to variations in morphology and birefringence which alters the amount of spin angular momentum that can be transferred from the beam.

figure: Fig. 4.

Fig. 4. A, a typical experimental AVPS (blue dots) of a µ1.3µm radius vaterite in water overlaid with the theoretical model (red line). The data has been windowed into 1 kHz bins and the dashed line indicates the determined angular velocity relaxation rate, ??. B, the estimated value of ?? of a µ1.1µm radius vaterite determined from varying measurement times over multiple 100 ms samples with varying trappings power.

As the hydrodynamic effect has a negligible contribution to the AVPS below 10 kHz, 2? can be estimated by averaging those values. There was a maximum error of 6% between the predicted value of ? and the fitted value providing an accurate estimation for the diffusion constant. Similarly, ?? was predicted from AVPS|?=??≈?/1.49 and was within 30% of the fitted value. Measurements where the value of ?? is determined with SNR<1 cannot be resolved adequately and were removed from the analysis. These limitations can be attributed to specific properties of the vaterite probes in those cases, mainly decreased birefringence, which disrupted the effectiveness of the polarization-sensitive measurements (<10% of probes). As ?? can be defined in terms of ?? along with knowledge of the densities, the hydrodynamic contributions were also parameterized with ?? ensuring only two fitting parameters are needed. The agreement of the fit experimentally verifies the rotational hydrodynamic model of a single microsphere in a viscous fluid in a new measurement bandwidth. The technique provides a unique opportunity to validate rotational hydrodynamic models approaching the MHz region for micrometer-scale systems.

figure: Fig. 5.

Fig. 5. A, experimental AVPS captured for the same probe in water with different trapping powers (µ1.8µm radius). B, the estimated values of ?? from fitting with a calibrated value of ?? compared to an uncalibrated model with ??=0. There was a mean difference of 0.07% between the models.

The SNR depends on the sensitivity in resolving angular displacements, which is largely determined by the birefringence of the vaterite probe. An ideal probe for our polarization-based detection is a strongly birefringent microsphere with low absorbance and reflectance. The spherical geometry simplifies the hydrodynamic model, while the latter properties are required to generate an alignment torque and measure the orientation from the polarization changes in the forward scattered light. Since the aim of rotational ballistic measurements is to resolve the angular velocity dissipation rate, lowering ?? away from the noise-dominated high-frequency region by increasing the probes radius or lowering the viscosity achieves this goal.

A key advantage of ballistic measurements is the short total measurement time required to retrieve a measurement of ??. Figure 4B shows the estimated value of ?? for increasing measurement times from the same experiment. Although the data is windowed into 1 kHz bins, measurement times longer than 1 ms are necessary to obtain sufficient time-averaging to adequately resolve the AVPS for fitting. A minimum measurement time of 30 ms was required to constrain the uncertainty of ?? within 10%. This minimum time varies slightly between probes based on their size and birefringence; however, we determined 50 ms of measurement time is sufficient to resolve ?? for all of our appropriate 2–4 µm sized probes.

The angular position trajectory in the ballistic regime is decoupled from the optical potential, which provides a major advantage for ballistic tweezers applications because the position relaxation time, ??, is significantly larger than the total measurement time. Hence, the optical trap in the AVPS has a negligible effect once ????. Figure 5A shows the AVPS from repeated ballistic tweezers measurements in water with trapping powers between 40 and 180 mW to shift the value of ??. Any variation between the spectra beyond 1 kHz is obscured by the noise, demonstrating that there is a negligible influence from changes in the angular trap stiffness. It is important to note that the increased trapping does not cause any observable heating as evidenced by the near-identical spectra. An increase in temperature would reduce the viscosity and shift the velocity thermalization corner at ?? toward lower frequencies.

This result can be exploited to allow for calibration-free measurements where the trap stiffness does not need to be known. This was verified in Fig. 5B which shows that the estimated values of ?? determined from fitting with a calibrated trap (where ?? is known) is equivalent to fitting with an extreme “uncalibrated” case assuming ??=0. The calibration-free fits remain accurate when the sample time is shorter than ??. The SNR is not affected by changes in the trap power because the angular position sensing comes from the independent beam allowing trapping with very low powers without loss of sensitivity. The optical potential is only required to confine the probe and control the equilibrium orientation.

figure: Fig. 6.

Fig. 6. A, experimental AVPS of 0% (blue) and 5% glycerol aqueous (red) solutions and the corresponding ?? values determined from fitting (circles) with µ1.8µm and µ1.4µm radius microspheres respectively. B, ballistic viscometry measurements of glycerol aqueous solutions (red) compared to literature values (blue crosses) [62]. The error bars are propagated from the parameter estimation and the uncertainty in radius.

Ultrafast viscometry is an ideal application of ballistic tweezers to take advantage of short measurement times and calibration-free analysis. Viscosity, ?, is usually defined in the long-time limit as a measure of a fluid’s resistance to flow. At low Reynold’s numbers, this is often linked to the diffusion rate or Stoke’s drag of a particle in the fluid. At short times, the hydrodynamic interactions modify the drag on the probe such that it still depends on viscosity but is no longer linearly proportional to velocity. However, it can be determined from the fast thermalization of angular velocity with ?=2??2???/15. It is important to reiterate that this result gives the same viscosity as defined in the long-time limit. As ballistic measurements have negligible influence from the optical trap, we have developed a rotational microviscometry technique decoupled from the optical potential and only requiring short measurement times to resolve.

We performed viscometry of water and aqueous glycerol solutions as proof-of-principle tests for rotational ballistic viscometry. Figure 6 shows two sample curves and compares the measured viscosity values to literature values [62]. The AVPS curves shown in Fig. 6A agree with the expected model. The variations in ? and ?? result from the different viscosities of the two samples and the probe sizes. The slight shift in the magnitude of the noise floor is due to differences in the optical power on the detector. There is excellent agreement for the estimated water viscosity and good agreement for the glycerol solutions, as shown in Fig. 6B. The increased uncertainty was partially caused by the higher ?? values that have lower SNR and the smaller radius of the ballistic probes. The systematic underestimation is likely a result of either a 1°C discrepancy or the hygroscopic behavior of glycerol, which reduces its viscosity. The uncertainty in the radius dominates the measurement uncertainty of ? as Δ?∝2Δ? and ? is seldom determined with better precision than 5% in vaterite trapping experiments [30]. However, this is still an improvement from previous ROTs experiments because those were performed in the diffusive regime where Δ?∝3Δ?. The SNR at ?? determines the ease with which the velocity dissipation rate is resolved and thus also contributes to Δ?.

The technique can be extended to viscoelastic fluids; however, it is unlikely to capture the elastic component because that normally becomes significant on longer time scales. Instead, the improved sensitivity will aid in the accuracy of those measurements and can still provide information about the short-time particle–liquid interactions. Coupling the ballistic measurements with existing active microrheometry would provide an unprecedented measurement bandwidth and sensitivity to perform microrheometry.

5. CONCLUSION

Access to rotational ballistic tweezers improves the speed at which fluid properties are determined, allowing the observation of the transition from diffusive to ballistic rotational motion. We access this regime through improvements in our experimental apparatus that combine measurement basis optimization with reduced detector noise, significantly improving the SNR from previous experiments. In addition to the trapping beam, we incorporate a stable measurement beam which allows for large variations in the trapping power without reducing the SNR. The polarization-sensitive technique is the optimal measurement basis for performing rotational ballistic measurements approaching the standard quantum limit, decoupled from the optical potential when performing ultrafast measurements.

We demonstrated that our methodology is capable of resolving angular velocity thermalization in liquids to access the ballistic regime allowing validation of hydrodynamic models and close-to-real-time measurements of viscosity. Employing this methodology to investigate systems out of equilibrium offers an exciting avenue of research. Applications into complex systems including studies of non-Newtonian fluids, non-equilibrium Brownian dynamics, biological processes, and those favoring rotational geometry, will greatly benefit from the improved sensitivity and fast measurement times of rotational ballistic tweezers measurements and allow for comparison with numerical predictions [30,63,64]. Combining rotational and translational ballistic techniques and extending the analysis to non-spherical particles would provide an opportunity to investigate force–torque coupling. An attractive example of such a measurement would be to investigate the translational dynamics and corresponding force on the surrounding fluid of a micromachine driven by an optical torque. This would provide insight into the dynamics of diverse optically driven nanomachines and micromachines, microswimmers, and other hydrodynamic models [16,53,6568]. Our enhanced measurement sensitivity and reduced measurement time have applications to biophysical measurements of cell–fluid interactions and in dynamic cellular environments where traditional approaches may be too slow and lack the spatial resolution to accurately capture the dynamics.