LETTERS PUBLISHED ONLINE: 20 MARCH 2011 | DOI: 10.1038/NPHYS1952 Millikelvin cooling of an optically trapped microsphere in vacuum Tongcang Li, Simon Kheifets and Mark G. Raizen* Cooling of micromechanical resonators towards the quantum mechanical ground state in their centre-of-mass motion has advanced rapidly in recent years1���8. This work is an important step towards the creation of ���Schr��dinger cats���, quantum superpositions of macroscopic observables, and the study of their destruction by decoherence. Here we report optical trapping of glass microspheres in vacuum with high oscillation frequencies, and cooling of the centre-of-mass motion from room temperature to a minimum temperature of about 1.5 mK. This new system eliminates the physical contact inherent to clamped cantilevers, and can allow ground-state cooling from room temperature9���15. More importantly, the optical trap can be switched off, allowing a microsphere to undergo free-fall in vacuum after cooling15. This is ideal for studying the gravitational state reduction16���19, a manifestation of the apparent conflict between general relativity and quantum mechanics16,20. A cooled optically trapped object in vacuum can also be used to search for non-Newtonian gravity forces at small scales21, measure the impact of a single air molecule14 and even produce Schr��dinger cats of living organisms9. Previous experiments demonstrated optical levitation of a 20-��m-diameter sphere in vacuum with a trapping frequency of about 20 Hz, as well as feedback control of a trapped sphere that was used to increase the trapping frequency to several hundred hertz and stabilize its position to within a fraction of one micrometre22,23. However, the resolution of its detection system23 was not sufficient to enable feedback cooling. According to the equipartition theorem, the root-mean-square (r.m.s.) amplitude of a trapped microsphere at thermal equilibrium is xr.m.s. = ��� kBT0/(m��2), where kB is the Boltzmann constant, T0 is the air temperature, m is the mass of the microsphere and �� is the angular trapping frequency. For a 20-��m-diameter sphere trapped at one-hundred hertz, the r.m.s. amplitude is about 0.04 ��m at 300 K, and will be much smaller at lower temperature. It is also important that the trapping frequency be much higher than the frequencies of seismic vibration to achieve significant cooling. We use a dual-beam optical tweezer to trap a 3.0-��m- diameter SiO2 sphere in vacuum with much higher oscillation frequencies (about 10 kHz) to minimize the effects of instrumental vibration. We also demonstrate a detection system to monitor the motion of a trapped microsphere with a sensitivity of about 39 fm Hz-1/2 over a wide frequency range. Using active feedback, we simultaneously cool the three centre-of-mass vibration modes of a microsphere from room temperature to millikelvin temperatures with a minimum mode temperature of 1.5 mK, which corresponds to the reduction of the r.m.s. amplitude of the microsphere from 6.7 nm to 15 pm for that mode. A simplified scheme of our optical trap and feedback set-up is shown in Fig. 1. The optical trap is similar to a trap used to Center for Nonlinear Dynamics and Department of Physics, The University of Texas at Austin, Austin, Texas 78712, USA. *e-mail: raizen@physics.utexas.edu. measure the instantaneous velocity of a Brownian particle in air24. It is created inside a vacuum chamber by two counter-propagating laser beams focused to the same point by two identical aspheric lenses with a focal length of 3.1 mm and a numerical aperture of 0.68. The wavelength of both trapping beams is 1,064 nm. They are orthogonally polarized, and are shifted in frequency to avoid interference. The beams are slightly elliptical and approximately form a harmonic trap with three fundamental vibration modes along the horizontal, vertical and axial directions, denoted X, Y and Z in Fig. 1a. The motion of a trapped bead causes deflection of both trapping beams. We monitor the position of the bead by measuring the deflection of one of the trapping beams with ultrahigh spatial and temporal resolution in all three dimensions (Supplementary Fig. S1). Using the position signal, we can calculate the instantaneous velocity of the bead, and implement feedback cooling by applying a force with a direction opposing the velocity (Fig. 1b). The feedback is generated by scattering forces from three orthogonal 532 nm laser beams along the axes as shown in Fig. 1a. The average intensity of the cooling beams is about 1% of the trapping beams. The optical power of each cooling beam is controlled by an acousto-optic modulator (AOM). Each beam is modulated with a time-varying signal proportional to the instantaneous velocity of the bead, added to an offset. The proportional component generates the required cooling force, and the offset slightly shifts the trap centre. An optically trapped microsphere in non-perfect vacuum will exhibit Brownian motion due to collisions between the microsphere and residual air molecules. When the microsphere is at thermal equilibrium, the power spectrum of centre-of-mass motion along each of the three fundamental mode axes is5,25: Sj ( ) = 2kBT0 m 00 (��j 2 - 2)2 + 2002 (1) where /2�� is the observation frequency, 00 is the viscous damping factor due to the air and ��j (j = 1,2,3) are the resonant frequencies of the optical trap along the X, Y and Z axes. The viscous damping factor due to air can be calculated by kinetic theory. Assuming the reflection of air molecules from the surface of a microsphere is diffusive, and the molecules thermalize with the surface during collisions, we have26: 00 = 6����r m 0.619 0.619+Kn (1+cK ) (2) where �� is the viscosity coefficient of the air, r is the radius of the microsphere and Kn = s/r is the Knudsen number. Here, s is the mean free path of the air molecules. cK = (0.31Kn)/(0.785 + 1.152Kn + Kn2) is a small positive function of NATURE PHYSICS | ADVANCE ONLINE PUBLICATION | www.nature.com/naturephysics 1 �� 2011 Macmillan Publishers Limited. All rights reserved.

LETTERS NATURE PHYSICS DOI: 10.1038/NPHYS1952 a s PBS s p AOM Vacuum RF d dt DAQ Y cooling beam s-polarized trap beam Z cooling beam X cooling beam p-polarized trap beam Position detector b Figure 1 | Schematic of feedback cooling of an optically trapped microsphere. a, Simplified schematic showing a glass microsphere trapped at the focus of a counter-propagating dual-beam optical tweezer, with three laser beams along the axes for cooling. The wavelengths of the trapping beams and the cooling beams are 1064 nm and 532 nm, respectively. b, Diagram of the feedback mechanism for the X axis. The position of a trapped microsphere is monitored by a home-built detection system. The position signal is sent through a band-pass filter (typically 100 Hz���300 kHz) and a derivative circuit (d/dt) to provide a signal proportional to velocity. This velocity signal is used to control the output power of a radiofrequency (RF) AOM driver that modulates the power of the X cooling beam. The data are digitized and stored on a computer by a data acquisition card (DAQ). PBS, polarizing beam splitter. Kn. At low pressures where Kn 1, the viscous damping factor is proportional to the pressure. The behaviour of the system with three-dimensional feedback cooling is straightforward to understand if we assume that there is no coupling between feedback forces and velocities in different directions. In this case, the feedback force in each direction adds an effective cold damping factor 0j fb, and the total damping becomes 0j tot = 00 +0j fb. The power spectrum of the motion of a trapped microsphere with feedback cooling can be described by equation (1) with T0 and 00 replaced by Tjfb and 0j tot, where Tjfb = T000/0j tot is the temperature of the motion with feedback cooling5. The motion can be cooled significantly by applying feedback damping 0j fb 00. The lowest temperature will be limited by the noise in the detection system and feedback circuits, as well as coupling between different directions. Figure 2 shows the linewidths, 00/2��, of the oscillation of a trapped 3-��m microsphere at different pressures without feedback cooling. The powers of the two trapping beams are 120 mW and 100 mW, respectively. The linewidths are obtained by fitting the measured power spectra with equation (1). The measured 10��1 10��1 100 100 101 101 10��2 10��5 101 102 102 103 103 104 105 Pressure (Pa) Frequency (kHz) 9.2 9.6 10 �� 0 /2 �� (Hz) 0.13 Pa S X (nm 2 Hz ��1 ) Figure 2 | Measured linewidths of the oscillation of an optically trapped 3-��m-diameter microsphere at different pressures. The blue curve is the prediction of a kinetic theory (equation (2)). The inset is the measured power spectrum at 0.13 Pa. By fitting the spectrum with equation (1) (red curve), we obtain ��1 = 2��(9756.4��0.3) Hz and 00 = 2��(0.46��0.06) Hz for this example. The same method is used to obtain linewidths for other pressures. Error bars indicate the standard errors in the fitting of the power spectra. linewidths agree very well with the prediction of kinetic theory (equation (2)) from 105 Pa down to 1 Pa. At pressures below 1 Pa, the measured linewidths are larger than the theoretical prediction. This linewidth broadening is due to power fluctuations of the trapping laser. The inset of Fig. 2 shows a power spectrum at 0.13 Pa. The trapping frequency ��1/2�� is 9,756.4��0.3 Hz, and the linewidth is 0.46 ��0.06 Hz, giving a quality factor (Qj = ��j/00) of 2.1��104. This implies the power fluctuation of the trapping laser is smaller than 0.01% during the measurement. An optically trapped microsphere provides a method to directly convert laser power to a frequency signal, which can be measured precisely. Stabilization of laser power to a trapped bead can find applications in laser physics, and can enable a more precise measurement of the Q for a second trapped bead in vacuum. Figure 3 shows experimental results of feedback cooling. Before feedback is turned on, the resonant frequencies ��j/2�� are 8,066��5 Hz, 9,095 �� 4 Hz and 2,072 �� 6 Hz for the fundamental modes at 637 Pa along the X, Y and Z axes, respectively. At this pressure, the peaks in the power spectrum due to the three fundamental modes are distinguishable, and heating effects due to the laser are negligible. We can therefore use the measured power spectra at 637 Pa to calibrate the position detectors for the fundamental modes at room temperature. After we turn on feedback cooling, the temperature of the Y mode changes from 297 to 24 K at 637 Pa. Then we reduce the air pressure while keeping the feedback gain almost constant thus, the heating rate due to collisions from air molecules decreases, while the cooling rate remains constant. As a result, the temperature of the motion drops. At 5.2 mPa, the mode temperatures are 150 �� 8 mK, 1.5 �� 0.2 mK and 68 �� 5 mK for the X, Y and Z modes. The mean thermal occupancy hni= kBTjfb/(����j h ) of the Y mode is reduced from about 6.8��108 at 297 K to about 3,400 at 1.5 mK. Here, �� h is the reduced Planck constant. Figure 3d shows the temperature of the three fundamental modes as a function of pressure. At low pressure and whenthe feedbackgainis constant,themode temperatureshould be proportional to the pressure, which is shown as a straight line with slope 1 in the figure. The temperature of the Y mode agrees with this prediction very well at pressures above 1 Pa. At our lowest temperatures, the power spectra are still much larger than the noise level, and the minimum temperature is achieved at pressures above the minimum pressure we can 2 NATURE PHYSICS | ADVANCE ONLINE PUBLICATION | www.nature.com/naturephysics �� 2011 Macmillan Publishers Limited. All rights reserved.