Nanoscale feedback control of six degrees of freedom of a near-sphere

Manipulating the rotational as well as the translational degrees of freedom of rigid bodies has been a crucial ingredient in diverse areas, from optically controlled micro-robots Būtaitė et al. (2019), navigation Zipfel (2000); Qin et al. (2018), and precision measurements Losurdo et al. (2001); Matichard et al. (2015) at macroscale to artificial and biological Brownian motors Erbas-Cakmak et al. (2015); Hu et al. (2018) at nanoscale. Recently, observing and controlling the translational motion of near-spherical nanoparticles levitated in vacuum has reached a low temperature regime from which their quantum propertes may be made evident Delić et al. (2020); Magrini et al. (2021); Tebbenjohanns et al. (2021); Kamba et al. (2022). However, manipulating all external degrees of freedom of nanoscale objects, desired for applications such as sensing and quantum mechanical studies Stickler et al. (2016); Millen et al. (2020); Millen and Stickler (2020); Schäfer et al. (2021), has been a challenging task.

The optical control of mechanical oscillators, typically possessing one translational degree of freedom, has been successful in exploring their motions at the quantum level, demonstrating various applications such as quantum transducers Andrews et al. (2014) and precision measurements Eichenfield et al. (2009). By levitating nanomechanical oscillators, one can expect an extremely high quality factor suitable for their coherent manipulations Millen et al. (2020); Millen and Stickler (2020); Gonzalez-Ballestero et al. (2021). At macroscale, mechanical motions of objects are detected by accelerometers and gyroscopes, where the sensitivities of these sensors matters. By contrast, the challenge at nanoscale is to detect all the minute motions with a precision sufficient for their manipulation. The remarkable progresses made with ground-state cooling of levitated objects have been particularly successful with nearly-spherical particles and focused on one of their three translational motions Delić et al. (2020); Magrini et al. (2021); Tebbenjohanns et al. (2021); Kamba et al. (2022), while other two translational degrees are still far from the quantum regime and three rotational degrees remain uncontrolled. It is just very recently that controlling multiple translational motions of nanoparticles in the quantum regime via cavity cooling became possible Piotrowski et al. (2023), while feedback cooling of multiple degrees of freedom near the quantum regime has been elusive.

The rotational degrees of freedom of levitated nanoparticles has attracted attention just recently Reimann et al. (2018); Stickler et al. (2021). In this context, highly anisotropic nanoparticles, such as nanodumbells and nanorods, are expected to be a promising system for exploring fundamental physics Hoang et al. (2016); Stickler et al. (2016); Kuhn et al. (2017); Ahn et al. (2018); Stickler et al. (2018); Ahn et al. (2020); Stickler et al. (2021); Schäfer et al. (2021). Cooling of up to two librational motions, oscillations like a physical pendulum, of nanodumbells, with Bang et al. (2020) and without van der Laan et al. (2021) translational cooling, and of one librational motion of micron-scale spheres Blakemore et al. (2022) have been reported. Very recently, cavity cooling of all the mechanical degrees of freedom of a nanodisk was also demonstrated Pontin et al. (2023).

Our work points out that nearly spherical nanoparticles, whose deviation from a sphere has been overlooked in previous studies, can also be a promising system for controlling all the mechanical degrees of freedom because of their slight deformations from a perfect sphere. With a tight optical confinement, such slight deformations are sufficient to enable us to observe their librational motions and characterize the shape, while the minimal deviation suppresses the decoherence rate of librational motions via photon scattering to significantly lower values than that of highly anisotropic nanoparticles. Furthermore, in comparison with the spectra of highly anisotropic particles Bang et al. (2020); Kuhn et al. (2017); Rademacher et al. (2022), the observed simple, narrow spectra of librational motions facilitate controls over their multiple degrees of freedom in the quantum regime.

In our experiments, we trap nearly spherical silica nanoparticles in a one-dimensional optical lattice formed in the vacuum chamber at a pressure of $6.5\times 10^{-7}$ Pa Iwasaki et al. (2019); Kamba et al. (2021) (Fig. 1). The nanoparticle is neutralized Frimmer et al. (2017); Monteiro et al. (2020); Kamba et al. (2022) and has an average radius of $174(3)\,\mathrm{nm}$. The translational motion along the optical lattice ($z$ direction) is cooled to an occupation number of $n_{z}=0.69(18)$ via optical feedback cooling realized by controlling the optical gradients Kamba et al. (2022); Vijayan et al. (2023), as shown in the power spectral density (PSD) obtained with photodetectors [Fig. 2 (a)]. We apply a similar approach also for cooling the translational motions in the $x$ and $y$ directions, in contrast to previous studies, where motions in two unfocused directions are cooled via parametric feedback cooling (PFC) Delić et al. (2020); Magrini et al. (2021); Tebbenjohanns et al. (2021); Kamba et al. (2022). We introduce additional two beams providing tunable optical gradients and modulate the intensity ratio of the two beams such that they exert feedback forces in both $x$ and $y$ directions. In this manner, the translational motions along the $x$ and $y$ directions are cooled to occupation numbers of $n_{x}=6(1)$ and $n_{y}=6(1)$, respectively, which are more than one order of magnitude lower than obtained with PFC Gieseler et al. (2012); Jain et al. (2016); Vovrosh et al. (2017); Iwasaki et al. (2019) [Fig. 2(b)].

In a frequency range between $10$ and 70 kHz, additional narrow peaks are observed [Fig. 2(c)]. These peaks are visible only at low pressures, where the broad spectra of the translational motions are minimized by feedback cooling. The frequencies are proportional to the square-root of the laser power and vary among trapped particles. We identify that these peaks originate from the three librational motions of the trapped nanoparticles and their higher order signals. These three peaks are correlated with each other in terms of both the amplitude and the frequency, suggesting rich dynamics such as gyroscopic and precession motions Rashid et al. (2018); Seberson and Robicheaux (2019); Bang et al. (2020); Seberson and Robicheaux (2020); Rudolph et al. (2021). Theoretically, it was suggested that librational motions are independent only in the deep-trapping regime where their amplitudes are sufficiently small Rudolph et al. (2021).

The observed librational motions arise from the orientational confinements in three rotational angles around the $x,y,z$ axes. In our study, there are two independent mechanisms for yielding such confinements. The first mechanism is a tendency that the longest axis of the particle is aligned with the polarization of the light Ruijgrok et al. (2011); Trojek et al. (2012), which has been experimentally observed with highly anisotropic particles such as nanodumbells Ruijgrok et al. (2011); Hoang et al. (2016); Bang et al. (2020); van der Laan et al. (2021). The second mechanism is an effect originating from the anisotropy of the optical trap. Under the generalized Rayleigh-Gans approximation, where we consider the inhomogeneous electric field of the trapping laser inside a nanoparticle and assume that the light scattering by the nanoparticle does not modify the local electric field, the mechanical potential energy of a non-spherical nanoparticle in an anisotropic trap $U$ is obtained by integrating the local interaction potential energy density over the volume of the nanoparticle Stickler et al. (2016); Seberson and Robicheaux (2020); Rudolph et al. (2021) (see Appendix for more details). $U$ is dependent on the relative angle between the nanoparticle orientation and the optical trap and is minimized when the longer axis of the particle is aligned to the orientation with the lower translational oscillation frequency (Fig. 1). In other words, the nanoparticle is aligned such that it experiences a higher light intensity. Theoretical studies on the librational motions of various geometries under this approximation have been reported Stickler et al. (2016); Papendell et al. (2017); Seberson and Robicheaux (2019); Seberson and Robicheaux (2020); Schäfer et al. (2021); Rudolph et al. (2021). To our knowledge, librational motions with the latter mechanism have not been experimentally observed.

Considering the two mechanisms, under an assumption that the trapped nanoparticle is an ellipsoid with radii of $p,q$ and $r$ as defined in Fig. 1(b), we obtain the expressions of the librational frequencies $f_{li0}$ around the $i$ axis when the oscillation amplitudes are small, where $i\in\{x,y,z\}$, as functions of $p,q$ and $r$ (see Appendix). Calculated librational frequencies with respect to $q/p$ show that they are sensitive to the anisotropy of the trapped nanoparticles [Fig. 3(a); see Fig. 6 for the plot with respect to $r/p$). Due to the strong confinement with the light polarization, the largest axis of the particle aligns with the light polarization. In addition, the strong trap anisotropy in the $yz$ plane, that is, $f_{ty}<f_{tz}$, suggests that the second elongated axis aligns with the $y$ axis. Thus, we expect that a configuration of $p>q>r$ provides the minimum potential energy.

Experimentally, we observe three frequencies of 13.9 kHz, 19.9 kHz, and 34.5 kHz when their amplitudes are sufficiently small. By minimizing the deviation between observed and calculated frequencies, we determine two radii to be $q=0.9963(9)p$ and $r=0.9917(9)p$, with which we can reproduce observed frequencies within $1.6\%$. The precision in determining $q/p,r/p$ is limited mainly by systematic uncertainties in the radius and the refractive index and is about $910\,\mathrm{ppm}$ (see Appendix). The precision can be comparable to the size of atoms because our measurement is based on the averaged interaction of light and many atoms. The demonstrated precision suggests a novel approach to precisely characterize the three dimensional shape of trapped nanoparticles without relying on electron microscopes.

We find that librational frequencies are also sensitive to the temperatures of the librational motions ([Fig. 3(c),(d),(e)]. This is because the depths of the potential energies for librational motions are of the order of $k_{\rm B}\times 10^{2}\,\mathrm{K}$, where $k_{\rm B}$ is the Boltzmann constant, and the librational motions can be readily excited to amplitudes that can experience the nonlinearity of the potential. Such a situation is in contrast to previous studies on the nonlinearity observed with the translational motions of nanoparticles Gieseler et al. (2013) and with the librational motions of nanodumbells Bang et al. (2020), where the potential depths are more than $k_{\rm B}\times 10^{4}\,\mathrm{K}$. Under an assumption that the time variation of the amplitudes of three librational motions are negligible, the average oscillation frequency with a finite amplitude for the $i$ direction can be written as

with $\beta=\sqrt{2-\sqrt{4-2k_{\rm B}T_{li}/(\pi^{2}I_{i}f_{li0}^{2})}}$, where $T_{li}$ and $I_{i}$ denote the temperature of the librational motion and the moment of inertia around the $i$ axis. As shown in Fig. 3b for $f_{lx}$, the derived expression (2) is in good agreement with experimentally observed frequency variations with respect to the area of the PSD, which is proportional to the motional temperature van der Laan et al. (2021). For the determination of the radii $q,r$, we use $f_{li0}$ obtained by the fits.

An important suggestion here is that such a measurement can provide a direct, independent thermometry of the librational motions, i.e. to obtain a conversion between a signal voltage and the temperature. This is because the extent of the nonlinearity directly reflects the absolute magnitude of the angular deviation. A similar idea has been employed for calibrating the translational motions of levitated particles Zheng et al. (2020); Tian et al. (2022). The temperature values obtained in this approach are shown in Fig. 3(c),(d),(e). Note that we observe considerable nonlinear frequency shifts at large libration amplitudes with temperatures of higher than $1\,\mathrm{K}$. In many previous works with nanoparticles, thermometry has relied on the thermalization at high pressures to the temperature of background gases Hebestreit et al. (2018a); Millen et al. (2020). However, establishing an independent method of thermometry is crucially important because thermalization measurements are always accompanied by large thermal fluctuations.

with $\beta=\sqrt{2-\sqrt{4-2k_{\rm B}T_{li}/(\pi^{2}I_{i}f_{li0}^{2})}}$, where $T_{li}$, $I_{i}$, and $f_{li0}$ denote the temperature of the librational motion, the inertial moment around the $i$ axis, and the librational frequency at a minimum amplitude. As shown in Fig. 3(b) for $f_{lx}$, the derived expression (2) is in good agreement with experimentally observed frequency variations with respect to the area of the PSD, which is proportional to the motional temperature van der Laan et al. (2021). For the determination of the radii $q,r$, we use $f_{li0}$ obtained by the fits.

We manipulate the orientation of the nanoparticle by applying electric fields on it. Even though the trapped nanoparticle is neutralized, they can have a charge distribution over its surface and/or inside its volume, yielding a finite dipole moment Afek et al. (2021). This fact implies that we can exert a feedback torque proportional to the angular velocity on nanoparticles by applying an electric field synchronized to the librational motion Blakemore et al. (2022). Such a cooling scheme is called cold damping and has been shown to be more efficient than parametric feedback cooling Tebbenjohanns et al. (2019); Conangla et al. (2019); Iwasaki et al. (2019). We observe that cooling is realized only when the phase of the applied electric field is chosen appropriately. When feedback electric fields include three independent signals synchronized to three librational frequencies, we are able to completely extinguish all the peaks from librational motions as well as their higher order signals (Fig. 2). Due to the correlation among the librational motions, we observe that extinguishing two librational motions decreases the remaining motion as well, but does not extinguish it perfectly. Such a behavior implies that the three librational motions become independent in the deep-trapping regime when their amplitudes are small Rudolph et al. (2021). The lowest temperatures are estimated to be lower than 0.03 K for each librational motion, limited only by the noise floor for observing the motions. The fundamental limit of cooling is determined by the compromise between the intrinsic heating rate, which will be discussed later, and the noise introduced by feedback Iwasaki et al. (2019); Tebbenjohanns et al. (2019); Conangla et al. (2019). The obtained temperatures are comparable to or lower than those obtained with nanodumbells Bang et al. (2020); van der Laan et al. (2021), which are realized with PFC.

To gain further insights on the dipole moment of nanoparticles, we explore the cooling dynamics. From the equation of the librational motion, we obtain an expression for describing the time evolution of $T_{li}$ in the presence of a feedback torque:

where $C_{i}=dE_{0}\eta_{i}/\sqrt{2k_{\rm B}I_{i}}$ is a damping rate due to feedback cooling with $d$, $E_{0}$, and $\eta_{i}$ being the dipole moment, the electric field amplitude, and the angle factor considering the angle between the dipole moment and the electric field (see Appendix).

A typical time evolution of $T_{lx}$, when the feedback signal is applied, is show in in Figure 4(a). We observe that the time variation under feedback cooling is in good agreement with the theory, suggesting that feedback cooling works as expected. We confirm that the damping rate due to feedback cooling is proportional to the applied electric field amplitude [Fig. 4(b)]. From the three values of the slopes $C_{i}/E_{0}$ (see Appendix), we deduce that the trapped nanoparticle has a dipole moment with a magnitude of $d=2p\times 1.92(13)e$, with $e$ being the elementary charge, and an orientation defined by $\theta=31(1)^{\circ}$ and $\phi=65(2)^{\cdot}$. The dipole moment is nearly constant for half a day. Given that the initial charge distribution before neutralization is of the order of $10e$ Kamba et al. (2022), the obtained magnitude is consistent with an interpretation that the dipole moment originates from a few elementary charges remaining over the surface even after neutralization. Exploring the stability of the dipole moment over a longer period may help understanding the origin of the dipole moment Blakemore et al. (2022).

We also explore the heating dynamics of the librational motions by observing the time variation of the amplitudes of the librational motions after feedback cooling is turned off [Fig. 4(c) for $f_{lx}$)]. After averaging over many experimental runs, we observe a linear increase in the temperature. The measured heating rates are $2.6(2)\,\mathrm{mK/s}$, $2.1(1)\,\mathrm{mK/s}$, and $2.3(1)\,\mathrm{mK/s}$ for $f_{lx}$, $f_{ly}$, and $f_{lz}$, respectively. The low heating rate is also reflected in the narrow spectral width of the PSD [Fig. 2(c)], where the observed width of around 10 Hz is not limited by the heating rate but rather reflects the Fourier limit and fluctuations in the laser intensity and in the amplitude of librational motions. These measured heating rates are two orders of magnitude lower than previously measured values for the librational motions of nanodumbells van der Laan et al. (2021) and typical heating rates for the translational motions of optically trapped nanoparticles Tebbenjohanns et al. (2019); Kamba et al. (2021), both of which are limited by photon scattering at high vacuum.

The observed slow heating dynamics reflects the fact that the photon recoil torque strongly depends on the geometry of the particle and is equal to zero for spherical particles Zhong and Robicheaux (2017). We compare the measured heating rates with calculated values obtained as the sum of photon recoil heating and background gas heating (see Appendix), as shown in Fig. 4(d), and find a good agreement. The agreement shows that heating is dominated by background gas collisions at the current pressure. The agreement between experiments and calculations also confirms the validity of the thermometry based on the nonlinearity of the trap.

Given that slow decoherence is a crucial ingredient for quantum mechanical experiments, identifying a highly coherent system is an important task. As shown in Fig. 4(d), by decreasing the pressure by two orders of magnitude, we expect to reach a regime where the decoherence of the librational motion is only limited by very slow photon recoil heating. At such a regime, number of coherent librational oscillations, during which the phonon occupation number is preserved, is expected to be more than 2000 for $f_{ly}$, which is more than two orders of magnitude larger than the value expected with the translational motions Delić et al. (2020).

The present study is important in the following aspects. First, even though the nearly spherical geometry does not seem optimal for observing and controlling the librational motions, we show that all the librational motions can be clearly observed and cooled to temperatures of below $0.03\,\mathrm{K}$. Second, we establish methods to characterize trapped nanoparticles precisely in terms of the geometry, the dipole moment, and the temperatures of librational motions. Third, because of the nearly spherical shape, heating of librational motions via photon recoils is negligible, and is only limited by very slow heating via background gases. Fourth, the higher order signals of librational motions often interfere with frequencies of the translational motions, thereby prohibiting efficient feedback cooling of the translational motions. We demonstrate that all the signals arising from the librational motions can be extinguished and are not an obstacle to cooling all the translational motions to near the ground state.

Characterizing the geometry of nanoscale objects has been a crucial issue in a wide variety of applications in biology, chemistry, physics, and engineering Blott and Pye (2008). Electron microscope imaging has been extensively employed Wang (2000). Optical measurements of the shape of trapped nanoparticles has been reported van der Laan et al. (2020); Kuhn et al. (2017); Rademacher et al. (2022). Our approach is particularly suited for near-spherical particles and provides a novel route to optically measure their shape with a precision of $0.09\,\mathrm{\%}$. Such a precision indicates discerning the difference in diameter of 0.3 nm, which is smaller than a recent demonstration in a related setup Rademacher et al. (2022) by one order of magnitude and is close to the precision of 0.12 nm obtained with latest electron microscopes Nakane et al. (2020). Our approach may be extended to particles with a geometry far from a sphere if an appropriate model for describing their motions in an optical trap is developed.

Our study is an important building block towards quantum mechanical experiments with levitated nanoparticles. There has been various proposals to observe quantum superposition states of the motions of nanoparticles, including approaches based on an optical diffraction grating Romero-Isart et al. (2011a) and quantum state tomography via time-of-flight measurements Romero-Isart et al. (2011b). However, given that the particle is not a perfect sphere, motions in other degrees of freedom, in particular librational motions, can readily obscure the minute effect of the motion cooled to the ground state. In this perspective, freezing all the degrees of freedom will be an essential ingredient in future studies Stickler et al. (2016); Millen et al. (2020); Millen and Stickler (2020); Schäfer et al. (2021).

The quantum mechanical behaviors associated with the librational motions are also an intriguing subject. It has been proposed that orientational quantum revivals can be observed with free-falling nanorods Stickler et al. (2018). The observed very slow heating rate of our system, combined with the demonstrated low temperatures of librational motions, which are lower than assumed in the proposal Stickler et al. (2018), suggests that nearly spherical nanoparticles can also be a promising candidate for investigating quantum physics with the orientation of nanoparticles. Although the temperatures of librational motions achieved in the present study are limited by the signal-to-noise ratio (SNR) of our experimental setup, enhancing the observation sensitivity via dedicated experimental improvements will enable us to further approach the quantum regime Toroš et al. (2018). The ability of feedback control on the rotational degrees of freedom also opens the door to studies on information thermodynamics with nanomechanical heat engines both at the classical and at the quantum regime Roulet et al. (2017).

Note added: After the submission of the present study, related works on cavity cooling of six degrees of freedom of a nanodisk Pontin et al. (2023) and on cavity cooling of the center-of-mass motions to the ground state in two dimensions Piotrowski et al. (2023) are published.