This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

High sensitivity, levitated microsphere apparatus for short-distance force measurements

Akio Kawasaki [email protected] Department of Physics, Stanford University, Stanford, California 94305, USA W. W. Hansen Experimental Physics Laboratory, Stanford University, Stanford, California 94305, USA    Alexander Fieguth Department of Physics, Stanford University, Stanford, California 94305, USA    Nadav Priel Department of Physics, Stanford University, Stanford, California 94305, USA    Charles P. Blakemore Department of Physics, Stanford University, Stanford, California 94305, USA    Denzal Martin Department of Physics, Stanford University, Stanford, California 94305, USA    Giorgio Gratta Department of Physics, Stanford University, Stanford, California 94305, USA W. W. Hansen Experimental Physics Laboratory, Stanford University, Stanford, California 94305, USA
Abstract

A high sensitivity force sensor based on dielectric microspheres in vacuum, optically trapped by a single, upward-propagating laser beam, is described. Off-axis parabolic mirrors are used both to focus the 1064 nm trapping beam and to recollimate it to provide information on the horizontal position of the microsphere. The vertical degree of freedom is readout by forming an interferometer between the light retroreflected by the microsphere and a reference beam, hence eliminating the need for auxiliary beams. The focus of the trapping beam has a 1/e2 radius of 3.2 μ\upmum and small non-Gaussian tails, suitable for bringing devices close to the trapped microsphere without disturbing the optical field. Electrodes surrounding the trapping region provide excellent control of the electric field, which can be used to drive the translational degrees of freedom of a charged microsphere and the rotational degrees of freedom of a neutral microsphere, coupling to its electric dipole moment. With this control, the charge state can be determined with single electron precision, the mass of individual microspheres can be measured, and empirical calibrations of the force sensitivity can be made for each microsphere. A force noise of <1×1017<1\times 10^{-17} N/Hz\sqrt{\rm Hz}, which is comparable to previous reports, is measured on all three degrees of freedom for 4.7 μ\upmum diameter, 84 pg silica microspheres. Various devices have been brought within 1.6μ1.6~{}\upmum of the surface of a trapped microsphere. Metrology in the trapping region is provided by two custom-designed microscopes providing views in the horizontal and one of the vertical planes. The apparatus opens the way to performing high sensitivity three-dimensional force measurements at short distance.

I Introduction

Force and acceleration measurements find numerous applications in science and technology. In recent times, traditional sensors based on mechanical springs have been complemented by more elaborate systems based on non-contact electromagnetic effects. These include atomic force microscopes Binnig et al. (1986), optical Abbott et al. (2009); Aasi et al. (2015); Acernese et al. (2014); Akutsu et al. (2019) and atom interferometers Peters et al. (1999, 2001); Biedermann et al. (2015); Hu et al. (2013), and falling corner-cube-gravimeters Niebauer et al. (1995), each emphasizing a different aspect of the measurement and generally achieving excellent sensitivity and accuracy. Since the pioneering work of Ashkin Ashkin (1970); Ashkin and Dziedzic (1971, 1977); Ashkin et al. (1986), optical tweezers have also been used as force sensors, with initial applications to particles in liquid suspensions, primarily in biology Fazal and Block (2011); Gutiérrez-Medina et al. (2010); Neuman and Block (2004) and polymer science Perkins et al. (1994). In the last decade, the path to trapping dielectric microspheres (MSs) in high vacuum was established Li et al. (2011), and various groups have used the technique to obtain highly sensitive force sensors that are well isolated from the environment Ranjit et al. (2016); Hempston et al. (2017); Blakemore et al. (2019a); Gieseler et al. (2012); Yin et al. (2013); Gieseler et al. (2013); Hempston et al. (2017); Monteiro et al. (2020); Ahn et al. (2020); Delić et al. (2020).

Force sensors using optically trapped MSs have the ability to carry out measurements at distances of sub-millimeter scale. This can be achieved by inducing a force between the MSs and specially designed sources of the desired interaction, placed in close proximity to the trapped MS. This has been demonstrated in a few cases Blakemore et al. (2019a); Rider et al. (2016); Winstone et al. (2018); Diehl et al. (2018), where separations between MSs and attractors have reached the scale of several micrometers.

In this paper, a multipurpose optical tweezer system, evolved from the apparatus described in Ref. Rider et al. (2018) and optimized to search for new fundamental interactions at the micrometer scale, is described. The system is currently used to trap silica MSs of diameters 4.7 and 7.6 μ\upmum and has achieved separations of 1.6 μ\upmum between the surfaces of a MS and a nearby device with nominal noise conditions. The trap uses a single beam of 1064 nm wavelength with interferometric readout on all three degrees of freedom, as demonstrated in Ref. Rider et al. (2018). A number of features, most notably a final focusing and recollimation employing off-axis parabolic mirrors, have been introduced to minimize beam halos at the focus, enabling closer access to the trapping region with minimal distortion of the optical field. A large vacuum chamber allows for the introduction of several motorized actuators, important for the manipulation of devices near the trapped MSs under vacuum. Three-dimensional metrology of these devices around the trap region is provided by two orthogonal microscopes. The readout of the polarization state of the trapping light after its interaction with the MS is used to measure the rotation of trapped MSs, owing to their residual birefringence. Rotation of the MS can be induced by producing a rotating electric field that couples with the permanent electric dipole moment generally present in the silica MSs used Rider et al. (2019). Trap stabilization against long term drifts of the interferometric platform affords a noise spectrum that is flat down to 1{\sim}1 Hz.

Along with its primary motivation to search for new interactions at the micrometer scale Geraci et al. (2010), the system described may be used for the investigation of Casimir forces Casimir (1948); Lamoreaux (1997); Mohideen and Roy (1998); Ederth (2000); Mohideen and Roy (1998); Bressi et al. (2002); Decca et al. (2003); Ether et al. (2015) and other applications Arvanitaki and Geraci (2013); Kawasaki (2020, 2019) requiring extreme force sensitivity.

II Optics Setup

The 1064 nm trapping light is produced using a distributed Bragg reflector laser (Innolume LD-1064-DBR-150) to seed a ytterbium-doped fiber amplifier (Thorlabs YDFA100P), resulting in a maximum power of 100 mW. The production of the trapping beam and the reference beams for the heterodyne detection system makes use of fiber optic components based on single-mode PM980 fiber or equivalent, as shown in Fig. 1. Light from the fiber amplifier first goes through a 50:50 fiber-coupled polarization maintaining (PM) beam splitter. The two output channels of this splitter are independently frequency shifted by 150{\sim}150 MHz to a final frequency difference of 125 kHz, using two fiber-coupled acousto-optic modulators (AOMs, Gooch and Housego T-M150-0.4C2G-3-F2P). One channel, the trapping beam, is then launched to free space for further manipulation, while the second channel is subsequently split into two halves to produce reference beams for the detection of the vertical (zz) and horizontal (xyx-y) MS positions. In order to passively stabilize its temperature, the entire fiber optics system is heat sunk to the 4,000 kg granite table on which the trap is located, and embedded in foam to decouple it from the ambient air, which help to reduce long term drift of the interferometric readout. Continuous temperature measurements show that the laboratory air conditioning system maintains the air temperature at 23.023.0~{}^{\circ}C±0.5\pm 0.5~{}^{\circ}C.

Refer to caption
Figure 1: Schematic depiction of the laser system: all of the components shown are fiber-coupled with PM fibers, and arrows indicate the fiber outcouplers to free space, shown in Fig. 2. DBR: distributed Bragg reflector laser and AOM: fiber-coupled acousto-optic modulator.
Refer to caption
Figure 2: Simplified free-space optics system: only the essential components are drawn, and some mirrors are omitted. Each reference beam has a pair of mirrors (not shown) between the fiber outcouplers and the NBSs to align their wavefronts to that of the trapping beam. The parabolic mirrors and MS are shown from a side view, while the remainder of the components are shown from a top view. The two auxiliary microscopes described in Section VI operate independently from the trap and are not shown here. BS: beam sampler, PBS: polarizing beam splitter, NBS: nonpolarizing beam splitter (50:50), HWP: half waveplate, QPD: quadrant photodiode, PD1-3: photodiodes, and XYM: xyx-y microscope.

Both the input and output free-space optical systems are each mounted on a 60×3060\times 30 cm2 breadboard. The breadboards themselves are the actuated elements of six-axis stages so that, once internal alignment between components on the breadboards is achieved, they can be collectively adjusted relative to the trap, which is directly mounted on the vacuum chamber on the granite table (see Section IV).

The trap and the surrounding optics are schematically illustrated in Fig. 2. The trapping beam is first launched to free space and collimated by a single aspheric lens (Thorlabs PAF2A-A15C), resulting in a 1/e2 beam radius of w=1.35w=1.35 mm. A part of the beam (9%{\sim}9\%) is sent to a photodiode (PD, Thorlabs DET100A2) by a beam sampler (CVI W1-IF-1012-UV-1064-0) in order to monitor the power of the trapping beam. The rest of the beam passes through an optical isolator consisting of two polarizing beam splitters (PBSs, CVI PBS-1064-100), a Faraday rotator (Electro-Optics Technology 110-10299-0001-ROT), and a half waveplate (Newport 10RP02-34). The light is horizontally polarized after the optical isolator.

The trapping beam is then sent to a piezoelectric deflector that steers the beam along two orthogonal axes and is used to provide feedback to the xyx-y translational degrees of freedom. The deflector, placed in a Fourier plane of the trap, employs a 7 mm diameter mirror (Edmund Optics #34-370) glued to a high bandwidth actuator (Thorlabs ASM003). The mechanical resonant frequency of this assembly is 2{\sim}2 kHz and feedback is applied predominantly at frequencies below 1{\sim}1 kHz. Finally, the beam is expanded with a 1:4 telescope using d=5d=5 cm and d=20d=20 cm lenses (Thorlabs AL2550H-B and Newport PAC32AR.16), with dd being the focal length, yielding w=4.84w=4.84 mm, slightly smaller than the expected 5.4 mm.

All vacuum chamber viewports are made from 5 cm diameter windows (CVI W2-IF-2037-UV-1064-0 for the trapping beam and CVI W2-IF-2037-UV-633-1064-0 for auxiliary imaging, see Section VI) custom-mounted with a tilt of 5 and sealed with O-rings. The trapping beam is focused inside the vacuum chamber by a d=50.8d=50.8 mm off-axis parabolic mirror (Edmund Optics #35-507). The mirror sits on a five-axis stage (Newport 9082-V), where rotation around the xx axis is the non-adjustable degree of freedom. Reflective optics are preferred here, owing to the absence of spurious reflections from optical interfaces in refractive optics. The 1/e2 beam radius at the trap focus w0w_{0} is estimated with a knife-edge method, making use of a nanofabricated device mounted on the main nanopositioning stage (see Section V). The measured beam size is w0=3.4w_{0}=3.4 μ\upmum (3.1 μ\upmum) in the xx (yy) direction, as shown in Fig. 3, averaged to 3.2 μ\upmum corresponding to a Rayleigh range of 31 μ\upmum. This is in reasonable agreement with the theoretical prediction of 3.5 μ\upmum from the measured numerical aperture of 0.095. A careful choice of the optical components, the use of 5 cm aperture optics after expanding the beam radius, and a careful alignment are important to achieve close to nominal performance. It is estimated that the residual imperfections and astigmatism (xx and yy foci displaced by 10μ{\sim}10~{}\upmum in zz) are dominated by nonideal alignment.

The diverging beam emerging from the trapping region is recollimated by an identical parabolic mirror, also mounted on the same type of five-axis stage as the first parabolic mirror. The trapping beam exits the vacuum chamber horizontally polarized, while vertically polarized light is separated with a PBS (Edmund Optics #65-606) and projected onto a PD (Thorlabs DET100A2) in order to monitor the MS rotation Rider et al. (2019). After being extracted from the vacuum chamber, the trapping beam passes through a 2:1 telescope composed of a d=20d=20 cm and a d=10d=10 cm lens (Newport PAC32AR.16 and Thorlabs LA1509-C), and is combined in a 50:50 non-polarizing beam splitter (NBS, Newport 10BC17MB.2) with a reference beam. The resulting superposition of light is projected onto a quadrant photodiode (QPD, Hamamatsu S5980). The 2:1 telescope coarsely matches the mode of light transmitted through the MS to that of the reference beam.

The counter propagating beam, retroreflected by the surface of the MS, is extracted from the rejection port of the optical isolator on the input optics side, and combined with the second reference beam using a 50:50 NBS (Newport 10BC17MB.2) and projected onto a PD (Thorlabs DET100A2).

Both reference beams are launched to free space by outcouplers with adjustable beam size and focal length (Thorlabs ZC618APC-C), which allow further optimization of the mode-matching to the transmitted and reflected components of the trapping beam. Prior to the NBSs, each reference beam undergoes two reflections so that the position and angle of incidence on the PD or the QPD can be adjusted.

As mentioned, the output optics are also mounted on a breadboard that allows for the overall alignment with respect to the trap. Both input and output free-space optics are enclosed in lens tubes to reduce air currents and microphonic effects (except for the two motorized mirrors in the reference beam path for the xyx-y position detection). The entirety of each of the six-axis stages holding the breadboards is enclosed in boxes made of acrylic to further suppress noise.

III Data Acquisition and Feedback

On the PD dedicated to the zz position of the MS (PD2 in Fig. 2), and the QPD dedicated to the xyx-y position of the MS, the incident optical power is modulated by the Δf=125\Delta f=125 kHz frequency shift between the trapping and reference beams. The five photocurrents corresponding to PD2 and each quadrant of the QPD are individually amplified and then digitized at a sampling frequency fs=500f_{s}=500 kS/s, exactly four times Δf\Delta f. By phase-locking the radio frequency (RF) synthesizers driving the AOMs to the master clock driving the analog-to-digital converters (ADCs), real-time estimates of the amplitude and phase of the oscillating photocurrent are obtained. Every sample AiA_{i} is spaced by a quarter wavelength, and thus,

Ai=GRtIphoto{sin(ϕ)i=1,5,cos(ϕ)i=2,6,sin(ϕ)i=3,7,cos(ϕ)i=4,8,,\displaystyle A_{i}=GR_{\rm t}I_{\rm photo}\begin{cases}\text{sin}(\phi)&i=1,5,...\\ -\text{cos}(\phi)&i=2,6,...\\ -\text{sin}(\phi)&i=3,7,...\\ \text{cos}(\phi)&i=4,8,...\end{cases}, (1)

where GG is the unitless voltage-gain of the amplifier circuit, RtR_{\rm t} is the transimpedance resistance used to convert the photocurrent IphotoI_{\rm photo} to a voltage, and ϕ\phi is an arbitrary, but fixed, phase offset between the photocurrent and the digitizer’s master clock. The phase and amplitude of this signal can be different for each quadrant and for the zz-position PD. The amplitude GRtIphotoGR_{\rm t}I_{\rm photo} is estimated as GRtIphoto=(Ai2+Ai12)1/2GR_{\rm t}I_{\rm photo}=(A_{i}^{2}+A_{i-1}^{2})^{1/2}, while the phase is estimated as ϕ=tan1[A1/(A2)]\phi=\tan^{-1}\left[A_{1}/(-A_{2})\right], or as ϕ=tan1[A3/(A2)]\phi=\tan^{-1}\left[-A_{3}/(-A_{2})\right] etc., where the appropriate negation and ratio repeats every fourth sample, as seen in Eq. (1). This procedure is often referred to as digital demodulation.

Since displacements of the MS in the horizontal (xyx-y) plane at the trap produce displacements of the transmitted trapping beam in the plane of the QPD, while the reference beam is fixed, the xyx-y degrees of freedom are read out as imbalances between the amplitude of interference photocurrent across the QPD, normalized by the total photocurrent amplitude. Vertical (zz) displacements of the MS change the optical path length of the retroreflected light, which are read out directly from changes in the phase of the interference photocurrent from PD2.

The aforementioned demodulation and construction of the xx, yy, and zz signals takes place within a field programmable gate array (FPGA), embedded with ADCs (National Instruments PXIe-7858R). This FPGA also computes and generates the stabilizing feedback on similarly embedded digital-to-analog converters, which are sent to the two orthogonal axes of the piezoelectric deflector for xx and yy and to the amplitude modulation port of the RF synthesizer driving the trapping beam AOM for zz. The FPGA is connected to a host computer to which it transfers the demodulated amplitudes and phases of PD2 and each QPD-quadrant, together with the position estimates and the generated feedback for offline analysis. While in the xx and yy directions interferometry is used only to suppress stray light from sources displaced from the trap center, it is an essential feature to measure the zz position that allows the trap to operate with a single beam, affording unimpeded access in the horizontal plane Rider et al. (2018); Blakemore et al. (2019a).

Temperature drifts affect the optical path length of all light used for interferometry. In the xx and yy degrees of freedom, the amplitude of interference photocurrent in each quadrant of the QPD is, at first order, independent of the phase ϕ\phi. Fluctuations in the optical power are also suppressed by normalizing the xx and yy estimates by the total photocurrent amplitude. In the zz degree of freedom, optical path length fluctuations propagate directly into the estimate of the zz position. These fluctuations are attributed to both air currents and residual temperature drifts and have, roughly, a 1/f1/f spectrum. They are suppressed passively by heat-sinking, the insulation, and the enclosures mentioned in the Sec. II, and actively by using the image reconstructed by the yzy-z microscope, described in Section VI. This latter technique addresses temperature fluctuations with periods of minutes.

Refer to caption
Figure 3: Beam profile at the trap for the xx and yy directions: the black lines are the data, and the orange curves are fits to Gaussian functions. The 1/e2 beam radii are 3.4 μ\upmum in the xx direction and 3.1 μ\upmum in the yy direction. During the measurement, the nanofabricated device approaches the beam from the negative side. Coordinates of the xx and yy axes are determined by strain gauges in the piezoelectrically driven flexures (See Section V).

IV Vacuum system and trap mechanics

The optical trap is located within a 42.5×42.5×34.2642.5\times 42.5\times 34.26 cm3 cubical aluminum vacuum chamber sealed on all six sides with International Standard Organization (ISO) flanges and Viton gaskets, which simplify access to the trap. High vacuum is achieved primarily with a 250 l/s turbomolecular pump (TMP, Pfeiffer HiPace 300), roughed by a scroll pump (Edward XDS35i) located in a separate room to reduce the acoustic noise. Although the use of Viton gaskets limits the attainable vacuum level to 108{\gtrsim}10^{-8} hPa, the actual base pressure of 2.4×1072.4\times 10^{-7} hPa achieved after a few days of pumping is thought to be limited by the outgassing of the motorized stages, some of which have stated vacuum compatibility of 10610^{-6} hPa. A system bakeout is impractical with the current setup. A high conductance port is available directly at the top of the chamber for the future installation of a large getter pump (SAES CAPACITORR HV1600) to improve the vacuum level. The current vacuum level does not limit the performance of the system, and data presented here are collected when the pressure reaches 106{\lesssim}10^{-6} hPa.

The bottom ISO-320 flange rests on an adapter to the granite table and contains a zero-length reducer from ISO-320 to ISO-100. The pumping system is connected to this flange and extends to the space below the granite table through a hole. A ceramic break (MDC 9632010) and a 10 cm long bellow are located between the pumping system and the main chamber for electrical and vibration isolation, respectively. The TMP can be isolated from the main chamber by a gate valve (MDC GV-4000M-P-01), while a custom tee between the chamber and gate valve connects a low conductance bypass for slow pumping at low vacuum with only the scroll pump. The bypass is throttled by a manual leak valve (Dunway VLVE-1000) and can be turned off by a pneumatically driven on-off valve (VAT 28324-GE41-0002/0068). Leaking N2 gas into the chamber is accomplished with an electrically controlled leak valve (MKS 148JA53CR1M) together with a pneumatically driven on-off valve (US Solid PSV00032).

The vacuum system is monitored by a full range vacuum gauge (Pfeiffer PKR 251), as well as a residual gas analyzer (MKS eVision+). However, the use of these devices is found to affect the charge state of the trapped MSs, and therefore, during charge-sensitive measurements, a capacitance manometer (MKS 627FU2TLE1B), with a minimum measurable pressure of 2×1052\times 10^{-5} hPa, is used to monitor the vacuum level. All gauges are mounted on conflat (CF) flange ports on an adapter nipple at the top of the main chamber. The rotational dynamics of the MS can also be used to measure the vacuum level, as has been demonstrated with the predecessor apparatus Blakemore et al. (2020).

Refer to caption
Figure 4: A photo of the components inside the vacuum chamber: the gold-coated cube in the middle (A) is the exterior of the six pyramidal electrodes surrounding the trapping region. The black diagonally-cut cylinder above the electrodes (B) is the recollimating parabolic mirror, which sends the collimated trapping beam to the PBS on the right (C). A small part of the optical surface of the focusing parabolic mirror is barely visible (D), under the cube. The motorized stage in the foreground of the PBS (E) supports the dropper (see Section VII), not installed in this image.

Components inside the vacuum chamber are mounted on a 25.4 mm lattice of 1/4-20 screw holes on the interior of the bottom ISO-320 flange. A view of the components inside the chamber is shown in Fig. 4. The trapping region is surrounded by six identical electrodes shaped as truncated pyramids. The electrode faces are 4.3 mm away from the trap center, forming a cubical cavity of 8.6 mm side, with narrow gaps between electrodes. Each electrode is hollowed out and, on the trap end, terminates with a 5.3 mm diameter aperture providing optical and mechanical access to the center. A cross section in the xyx-y plane and at the nominal zz position of the trap is shown in Fig. 5. In addition to the shielding against stray electric fields, each electrode is electrically isolated and can be independently biased in order to apply electrical forces and torques to the MSs. Torque can be applied with a constant magnitude rotating electric field, generated by four phased sinusoids on four coplanar electrodes, that couples to the residual electric dipole moment generally present in the silica MSs used here Rider et al. (2019); Blakemore et al. (2020). The electric field within the trapping region is calculated using finite element analysis for any configuration of electrode biases. The entire electrode assembly and mounting structure are composed of 6061 aluminum alloy and gold coated by electrode-poised plating. Surfaces facing the trap region, as well as the conical cavities of the top and bottom electrodes through which the trapping beam propagates, are further coated with colloidal graphite (Electron Microscopy Sciences 12660) to reduce the scattering of stray light.

V Nanopositioning Stage Ports

The system is designed to provide stable and reproducible access to the trapping region for devices with dimensions in the μ\upmum to mm range. Mechanical access is realized through the holes in two of the four electrodes in the horizontal plane, shown at the bottom and right of Fig. 5. In the current configuration, these house the main nanopositioning stage and an auxiliary nanopositioning stage, respectively.

Refer to caption
Figure 5: Components surrounding the trap in the xyx-y plane: the main NP stage carries the nanofabricated device described in Section V mounted on its end. The y-z microscope objective is described in Section VI, and the fiber for UV light is described in Section VII. (Aux) NP stage: (auxiliary) nanopositioning stage.

The main nanopositioning stage consists of a stack of piezoelectrically driven flexures (Aerotech QNP40Z-100 for zz mounted on top of Aerotech QNP60XY-500 for xyx-y), which has a full range of 500μ500~{}\upmum in xx and yy, and 100μ100~{}\upmum in zz, with a resolution better than 1 nm and reproducibility of 2{\sim}2 nm. The actual positions along each axis are measured by strain gauges within the actuator and subsequently recorded together with the MS position. This actuator is mounted on a custom six-axis stage for coarse alignment to the trap center, with a long travel along the xx axis, required for insertion of the primary device into the trapping region. Four degrees of freedom are manually operated, and one degree of freedom is actuated with a piezoelectric motor (Newport 8301-UHV), while the translation along the xx axis is accomplished with a DC servo motor with 12 mm range (Thorlabs Z812V) for high repeatability. On top of the stack of piezoelectrically driven flexures, a gold-coated aluminum conical cantilever is mounted. To reduce mechanical load and optimize the bandwidth of the flexure’s motion, the conical cantilever is hollowed out with a final wall thickness of 0.66 mm resulting in a mass of 3.13 g. The entire assembly is shown in Fig. 6. With the cantilever installed, the piezoelectrically driven flexure has a measured bandwidth of 80(100){\sim}80~{}(100) Hz in the xx (yy) direction, limited by mechanical resonances, and an estimated bandwidth of 500{\sim}500 Hz in the zz direction.

The auxiliary nanopositioning stage consists of a three-axis piezoelectrically driven flexure (Newport NPXYZ100SGV6) with 100μ100~{}\upmum range in each of the three orthogonal directions. For insertion into the electrode structure, a motorized stage with 12 mm range (Newport AG-LS25V6) is mounted on top of the three-axis actuator. A second, identical conical cantilever is mounted on this motorized stage for a secondary device to be inserted into the the trap orthogonally to the main device. Further coarse alignment and angular adjustment are provided by another manually adjusted platform (Newport 9071-V) onto which this assembly is mounted. The secondary system is intended for static use.

Refer to caption
Figure 6: CAD model of the conical cantilever (A) onto which the primary silicon device is mounted, the piezoelectrically driven flexures (B), and the custom six-axis stage (C) for coarse alignment and device insertion.

Nanofabricated devices, e.g. Ref. Wang et al. (2017) are mounted at the ends of the two cantilevers with a conductive epoxy (Epo-Tek H21D). Typical devices have a high aspect ratio, being wide and long in the two horizontal dimensions (1{\sim}1~{}mm) and thin (10μ~{}\upmum to 25μ~{}\upmum) in the vertical dimension, in order to minimize their effect on the converging and diverging trapping beam. Devices are typically gold-coated to optimally define their electrical potentials (although charge patches are still present Blakemore et al. (2019a)). Both devices can be independently biased, allowing the production of electric fields with large gradients in the immediate vicinity of a trapped MS Blakemore et al. (2019a); Rider et al. (2016). Each can also be used as a knife edge to scan across the beam and characterize its radius, center position, and focal point, as noted in Section II (see Fig. 3) and demonstrated with a previous system Rider et al. (2018). Typically with this technique, the center position of the trapping beam relative to the nanofabricated device is determined with a precision of 0.1μ{\sim}0.1~{}\upmum, where the device’s position is measured with strain gauges within the piezoelectrically driven flexure on which it is mounted. The coordinate system of the device is then registered to the position at which the edge of the device crosses the beam focus (See Fig. 3). Additionally, the size of these particular MSs is known with a precision of 0.04 μ\upmum [49]. Thus, the distance between the surface of the sphere and the surface of the nanofabricated device is determined by the relative positions of the trap and the nearby device, and the size of the MS, assuming the MS is trapped at the beam focus. The closest separation achieved between the surface of a MS and a nearby device is 1.6 μ\upmum, at which position the device can be translated in front of a stably trapped MS. The orientation and position of the devices on the two nanopositioning stages relative to the trap are determined by a combination of their interactions with the trapping beam and the auxiliary yzy-z and xyx-y microscopes, as discussed in Sec. VI. The relative position between the two devices can be registered by bringing them into contact, which is clearly visible by their elastic deformation. This has been found to be reproducible to within 1μ{\lesssim}1~{}\upmum.

VI Auxiliary Imaging Systems

The setup is equipped with two separate systems for auxiliary imaging and metrology near the trapping region. The primary system images the yzy-z plane of the trap through one of the electrodes in the xyx-y plane, as shown in Fig. 5. An infinity-corrected microscope objective (Nikon N10X-PF) is mounted, in vacuum, on a motorized linear stage with 12 mm range (Newport AG-LS25V6) to adjust the focal plane. The image from the objective is extracted through a viewport and focused by a d=20d=20 cm lens onto a CMOS camera (Allied Vision MAKO U-030B). A trapped MS is visible in Fig. 7 due to scattered light from the 1064 nm trapping beam. In the same image, the nanofabricated devices near the trap are illuminated with an 870 nm wavelength light-emitting diode (LED) (Thorlabs LED870E), injected into the imaging system through a 50:50 NBS (Newport 10BC17MB.2) located in the path between the microscope objective within the vacuum chamber and the d=20d=20 cm lens. Additionally, filters can be inserted in the optical path to selectively attenuate different wavelengths and optimize the visibility of various components. The combination of the illumination and a filter of 1064 nm light provides an additional way of imaging the MS, which casts a shadow on the device surfaces behind it.

While one pixel of the camera spans 0.5μ0.5~{}\upmum in the image plane, the resolution of the imaging system is close to the diffraction limit (1.5μ1.5~{}\upmum for 870 nm light), verified using a USAF1951 resolution test chart. The entire imaging system can be cross-calibrated into the coordinate system of the piezoelectrically driven flexures by imaging objects of known size (e.g. the nanofabricated devices) within the trapping region.

This yzy-z microscope is also used in a slow loop of the feedback maintaining a constant vertical position of the MS. The stability of this slow feedback is better than 1μ1~{}\upmum with a single image with an exposure time of 44 μ\upmus taken every 10 s.

The second imaging system provides a view of the xyx-y plane at the zz position of the trap through the path for vertically polarized light, as shown in Fig. 2. Half of the vertically polarized light from the output parabolic mirror is focused through a d=7.5d=7.5 cm lens into an infinity-corrected objective, mounted on a manual linear translational stage to adjust the focal plane. The system is designed such that the objective can be switched to any model with a diameter of less than 3 cm, though a ×10\times 10 magnification objective (Reichert Neoplan 1754) with a tube length of 16 cm is currently in use. Therefore, a broad range of magnifications are achievable. The image from the objective is focused by a second d=3d=3 cm lens onto another CMOS camera (Allied Vision MAKO U-029B). The field of view is illuminated by 870 nm light produced by an LED (Thorlabs LED870E) and injected in the system through a NBS. The resolution in the center of the image is 1.5μ\sim 1.5~{}\upmum, near the diffraction limit. However, substantial aberrations are present away from the center of the image, owing to the lack of correction for the parabolic element and possible imperfections in the system alignment. For this reason, the xyx-y imaging is used only for qualitative assessment and rough alignment. An improved microscope, using custom correction optics, may be installed at a later time.

Refer to caption
Figure 7: Composite image of two nanofabricated devices in the vicinity of the trap, together with a trapped MS, captured by the yzy-z microscope. Device 1 is a silicon and gold beam Wang et al. (2017), mounted on the main nanopositioning stage. Device 2, mounted on the auxiliary nanopositioning stage, is a 1000μ1000~{}\upmum long all-silicon device, which has an L-shaped cross-section when viewed in the xzx-z plane but not visible here, which usually houses device 1, shielding the trapped MS from background forces associated with device 1. In this image, device 1 is translated vertically from its nominal position and is out of focus, the latter reducing the apparent vertical extent. The main frame of the image shows the scattering of the trapping laser by the MS (A), which leads to saturation of several pixels of the camera. The inset is taken with a notch filter that has an optical depth of 6 for 1064 nm light (Thorlabs NF1064-44), in order to demonstrate the shadow of the MS (B) blocking the 870 nm illumination of device 2. The bright and blurry region (C) behind the device is caused by a reflection of the illumination from another part of device 1.

VII Trap Operation

Dielectric MSs are prepared by rubbing a fused silica cover slip (dropper) on a powder of MSs laying on a sheet of glass. The dropper is glued to a piezoelectric transducer (Thorlabs PA4DG) and inserted between the recollimating parabolic mirror and the electrodes, with the face coated with MSs oriented down. The dropper is then vibrated by driving the piezoelectric transducer with an oscillating voltage, chirped between 150 kHz and 400 kHz. The MSs, held on the surface of the dropper by van der Waals forces, are released by the vibration, which is expected to generate kilo-gg or larger accelerations Li (2011). One dropper can be used to refill the trap many times, provided that the drive amplitude is gradually increased, suggesting that MSs are bound to the dropper with varying force. Between the initial preparation of the dropper and its depletion, the RF power driving the transducer typically has to increase from 1 mW to 1 W.

Refer to caption
Figure 8: Typical discharging process: each data point corresponds to a 10-second measurement of the charge state (i.e. the amplitude of a MS response to a driving field), between which a certain number of flashes of the UV lamp occur. Quantization of the charge state in units of ee is observed. When near the neutral charge state, the resolution of the measurement can be increased by increasing the driving voltage.

The loading process has a small but finite efficiency as most MSs fall without being captured in the trap. Typically, 6{\sim}6 hPa pressure of N2 is used to slow the falling MSs and make trapping possible, together with an increased laser power to increase the depth of the optical trap. For example, 4.7μ4.7~{}\upmum diameter silica MSs Bangs Laboratories , with mass 84{\simeq}84 pg Blakemore et al. (2019b) are caught with 9{\sim}9 mW of power in the trapping beam, while they are maintained at the focal point of the trap with 1.7{\sim}1.7 mW. Additionally, 7.6 μ\upmum diameter, 420 pg mass silica MSs microParticles GmbH are stably trapped with 15{\sim}15 mW of power and 18{\sim}18 hPa of N2 pressure. It is expected that even larger MSs could be trapped with increased laser power. In order to minimize accumulation of MSs on the focusing parabolic mirror, a second, larger silica plate (catcher) is inserted between the electrode and the focusing parabolic mirror during loading. The dropper and catcher are inserted and removed independently using motorized stages (Newport AG-LS25-27V6). The trapped MS remains stable in the trap during the removal of the catcher, which is positioned 27 mm below the trap’s focus, without any special effort.

Once a MS is trapped, the low conductance bypass system is used to slowly pump down the vacuum chamber to 0.5 hPa (typically over 25 min), the pressure at which the feedback is initialized, typically with a reduced power of 1.8 mW (for 4.7μ4.7~{}\upmum diameter MSs). Feedback from the position measurements of the three degrees of freedom of the MS is then applied to the laser power (zz) and the pieozelectric deflector (xyx-y). The feedback for the vertical (zz) degree of freedom includes proportional, integral, and derivative terms, and is applied in addition to an independent laser power stabilization which is primarily an integral term together with a proportional term. The slow, zz feedback mentioned in Section VI is only applied during long term measurements that typically last hours or days. The feedback for the horizontal (xyx-y) degrees of freedom is applied after the zz direction is stable and includes only derivative terms in order to damp the MS’s resonant motion in the horizontal plane.

Refer to caption
Figure 9: Transfer functions of a typical MS. Shown are the magnitude (left) and phase (right) of the complex, frequency-dependent matrix 𝐇ij(f)\mathbf{H}_{ij}(f) discussed in Section VIII. The orange solid curves for (i,j)=(x,x),(y,y)(i,j)=(x,x),(y,y) are fits of the measured response to that of a damped harmonic oscillator, whereas the dashed line for (i,j)=(z,z)(i,j)=(z,z) is a quadratic spline interpolation. Off-diagonal interpolations are not shown. Resonant frequencies fiif_{ii} and damping coefficients Γii\Gamma_{ii} obtained from the fit are fxx=301f_{xx}=301 Hz±1\pm 1 Hz, Γxx=50\Gamma_{xx}=50 Hz±2\pm 2 Hz, fyy=292f_{yy}=292 Hz±1\pm 1 Hz, and Γyy=27\Gamma_{yy}=27 Hz±2\pm 2 Hz.

The zz direction stabilization also allows the zz position of the MS to be varied with respect to the trap focus and other mechanical devices in its proximity. At the focus, where the trapping beam intensity is greatest, the optical spring constant confining the MS in the xx and the yy degrees of freedom has a maximum. The stochastic force on the MS from the 0.5 hPa of residual gas impacting the MS has a white frequency spectrum that drives these degrees of freedom, resulting in MS motion with a frequency spectrum well modeled with a driven, damped harmonic oscillator with a clearly observable resonant frequency. For a given zz position, the slight astigmatism produces two distinct values of fxf_{x} (fyf_{y}), the resonant frequencies for the xx (yy) directions, and the final setpoint is determined empirically to minimize the difference between fxf_{x} and fyf_{y}. For the 4.7μ4.7~{}\upmum diameter silica MSs, typical values are fxfy300f_{x}\simeq f_{y}\simeq 300 Hz and |fxfy|10|f_{x}-f_{y}|\simeq 10 Hz. At this location, the harmonic trap for zz is generated by the feedback, with fz300f_{z}\simeq 300 Hz.

With xx, yy, and zz feedbacks on, the vacuum chamber is further pumped down, first with the low conductance bypass and then by opening the gate valve and starting the TMP at <0.1<0.1 hPa. The system reaches a pressure of 106{\sim}10^{-6} hPa in a few minutes, at which point the noise floor of the MS position is no longer dominated by the Brownian motion of the MS. MSs trapped in this fashion have been observed to be stable in the trap indefinitely in the absence of external disturbances.

The charge state of the trapped MS is measured by applying a sinusoidally oscillating voltage difference to a pair of opposing electrodes, and measuring the amplitude of the MS response to this driving field. For sufficiently long integration times, charge quantization is observed with a signal to noise ratio of 20{\sim}20. By increasing the amplitude of the driving field, the neutral state can be characterized exceedingly well as demonstrated in a simpler version of the trap in Ref. Moore et al. (2014).

As they are loaded in the trap, MSs can have either positive or negative overall charge. The charge state of MSs can be changed by flashing ultraviolet (UV) light from a xenon flash lamp (Hamamatsu L9455-13) into the trapping region. The light is brought into the trap with a multimode solarization-resistant UV fiber (Thorlabs UM22-600) and coupled to free space with a 4 mm fused silica ball lens (Edmund Optics #67-385), as shown in the left of Fig. 5. The fiber is brought into the vacuum chamber using the feedthrough discussed in Ref. Abraham and Cornell (1998). This system is capable of both increasing and decreasing the MS charge state. Pulsing the UV lamp with nothing in close proximity to a MS tends to eject electrons from the MS, yielding a more positive charge state. If the sputtered gold surface of a nanofabricated device is placed behind the MS, more electrons appear to be ejected from the gold surface, changing the charge of the MS in the opposite direction. MS charge states upwards of ±500e\pm 500~{}e, with ee the fundamental charge can be obtained in this way or, if desired, net neutrality can be achieved. A typical discharging cycle is shown in Fig. 8. Since conductive structures close to the MS distort the electric field, absolute charge calibration is only performed as the charge increases by removing electrons from the MS, when the nanofabricated device is fully retracted.

VIII Force Calibration

With charge quantization, it is possible to empirically calibrate the response of the MS without assumptions. This is achieved by applying a sinusoidally oscillating electric field E=E0sin(2πft)E=E_{0}\sin(2\pi ft) oriented along the degree of freedom to be calibrated, with E0E_{0} being the amplitude of the electric field and ff being the frequency of the oscillation. The electric field amplitude E0E_{0} at the MS location is calculated by finite element analysis, starting from the electrode geometry and the applied voltage. The force applied is then F=qMSEF=q_{MS}E, where qMS=keq_{MS}=ke with ee being the fundamental charge is deduced by counting the number of quanta kk. Simultaneously, the response of the MS, RR, as determined by the imaging and demodulation procedure described previously, is also known. Thus, for a given frequency of an oscillating electric field, the ratio R/FR/F between the amplitude of MS response in arbitrary units and the amplitude of the applied force in physical units is derived. The same procedure also extracts any phase shift of the response relative to the drive. This can be done independently for each degree of freedom xx, yy, and zz.

To measure the frequency dependence of this calibration, the electric field is applied in the form of a frequency comb

E(t)=E0n=1Nsin(2πnft+ϕn)E(t)=E_{0}\sum^{N}_{n=1}\sin\left(2\pi nft+\phi_{n}\right) (2)

with ϕn=2πn2/N\phi_{n}=2\pi n^{2}/N being a phase shift, implemented in order to avoid large spikes of the electric field (consider the Fourier series of a delta function). Typically, f=7f=7~{}Hz, N=100N=100, and E0=100E_{0}=100 V/m are used for the initial characterization of the system. This electric field drive and the corresponding MS response can be continued indefinitely and averaged until a desired precision is achieved.

This method also provides a means to measure any cross talk inherent to the imaging, as well as the trap itself. With the symmetric electrode geometry employed here, forces are induced along specific coordinate axes, while the MS response along all three degrees of freedom is measured. Writing the applied force as Fj=qMSEjF_{j}=q_{\rm MS}E_{j} where j{x,y,z}j\in\{x,y,z\}, and EjE_{j} of the form of Eq. (2), and writing the MS response as RiR_{i} for i{x,y,z}i\in\{x,y,z\}, each has a Fourier transform F~j(f)\widetilde{F}_{j}(f) and R~i(f)\widetilde{R}_{i}(f), respectively. The frequency dependent transfer function 𝐇ij(f)\mathbf{H}_{ij}(f) is calculated as 𝐇ij(f)=R~i(f)/F~j(f)\mathbf{H}_{ij}(f)=\widetilde{R}_{i}(f)/\widetilde{F}_{j}(f). The complex-valued 𝐇ij(f)\mathbf{H}_{ij}(f) naturally includes amplitude ratios as well as phase shifts for each drive frequency.

For a given ijij, the amplitude and phase of 𝐇ij(f)\mathbf{H}_{ij}(f) are smoothly interpolated (and extrapolated under certain assumptions) for frequencies not explicitly included in Eq. (2). Then, for any given measurement, the force applied to the MS is first constructed in the frequency domain from the measured MS response as F~j(f)=𝐇ij1(f)R~i(f)\widetilde{F}_{j}(f)=\mathbf{H}^{-1}_{ij}(f)\widetilde{R}_{i}(f). With sufficient optical alignment, the off-diagonal components of 𝐇ij\mathbf{H}_{ij} are typically five times smaller than their diagonal components. An estimate of the physical displacement of the MS at low frequencies is then obtained as x~j(f)=(1/mMSωj2)F~j(f)\widetilde{x}_{j}(f)=(1/m_{\rm MS}\omega_{j}^{2})\widetilde{F}_{j}(f), with ωj\omega_{j} being the angular resonant frequency for a particular coordinate axis, expressed in rad/s.

An example of 𝐇ij\mathbf{H}_{ij} is shown in Fig. 9. In both the xx and the yy (horizontal) directions, the optical spring constant kj=mMSωj2k_{j}=m_{\rm MS}\omega_{j}^{2} is 3×107\sim 3\times 10^{-7}~{}N/m, while the translational damping constant Γ\Gamma is approximately 50\sim 50 Hz and can be adjusted over a wide range by the active feedback. The resonant frequency and optical spring constant are in agreement with estimates from the Optical Tweezers Computational Toolbox Nieminen et al. (2007), which assumes perfectly spherical MSs with uniform density, as well as an ideal Gaussian beam and diffraction-limited focus. In the zz (vertical) direction, relatively large proportional and integral gain terms are required to stabilize the MS’s position, and a traditional damped harmonic oscillator model is insufficient to describe the frequency spectrum of the MS motion.

Between successive 10-s integrations, the response of the MS is stable to within 5%, and reproducible on the timescale of days. The response of a charged MS to an oscillating electric field is linear up to 2×1013\sim 2\times 10^{-13}~{}N of applied force, over a wide range of charge states |qMS|500e|q_{\rm MS}|\leq 500~{}e. The linearity of the system described here matches or exceeds that of the previous version of the apparatus Blakemore et al. (2019a).

IX Sensitivity

Refer to caption
Figure 10: Force (left axis) and acceleration (right axis) sensitivity of the system using 4.7 μ\upmum diameter MSs in the xx (upper panel), yy (middle panel), and zz (lower panel) directions, obtained by the procedure described in Section IX. The power spectral density SFjS_{Fj} is estimated from the complex square of the Fourier transformation: F~jF~j\widetilde{F}_{j}^{\ast}\widetilde{F}_{j}.

Utilizing the transfer function measurement described above, the noise spectrum of a neutral, 4.7 μ\upmum diameter MS without any external driving field is converted into physical units (N/Hz{\rm N}/\sqrt{\rm Hz}), as shown in Fig. 10. This illustrates the typical force sensitivity of the system, reaching a level of 1×1017N/Hz1\times 10^{-17}~{}{\rm N}/\sqrt{\rm Hz}. This force sensitivity is translated to an acceleration sensitivity of 12μg/Hz12~{}\upmu g/\sqrt{\rm Hz}, with g=9.8m/s2g=9.8~{}{\rm m/s}^{2}. These values are comparable to the sensitivity previously reported Rider et al. (2018). The observed noise floor is of additive, white Gaussian nature, and thus the sensitivity for any signal would improve with t\sqrt{t}, with tt being the integration time, in the absence of background forces.

The expected Brownian motion of the MS at the vacuum level achieved is estimated to be 2×1019N/Hz2\times 10^{-19}~{}{\rm N}/\sqrt{\rm Hz}. In addition, it is demonstrated that the noise floor is the same at different pressures, below 105\sim 10^{-5} hPa. The shot noise is subdominant as well and estimated to be of the order of 3×1019N/Hz3\times 10^{-19}~{}{\rm N}/\sqrt{\rm Hz} Rider et al. (2018). Hence the observed noise floor is likely technical in nature, dominated by pointing fluctuations of the trapping and reference beams or due to electronics noise of the detection system. The improvement in low frequency 1/f1/f noise as compared to Ref. Rider et al. (2018) is mainly due to the enclosure of the free-space optics, resulting in reduced air currents and microphonics.

X Conclusion

A system using a single optical tweezer to trap silica microspheres with diameters (masses) of 4.7 and 7.6 μ\upmum (84 and 420 pg), intended to measure the force between the microsphere and another device in close proximity, is described. The small radial extent of the beam at the trap focus and the lack of significant non-Gaussian tails opens the possibility of bringing mechanical devices, such as thin silicon cantilevers, as close as 1.6 μ\upmum from the surface of the microsphere. Electrodes surrounding the trap provide shielding, allow control over the translational degrees of freedom of charged microspheres and the rotational degrees of freedom of neutral microspheres, coupling to their electric dipole moment. The same electrode system allows for the calibration of the system’s response to forces applied to the microsphere. The observed force sensitivity of 1×10171\times 10^{-17} N/Hz{\rm N}/\sqrt{\rm Hz} is limited by non-fundamental noise sources, so that further improvements are expected in the future. Auxiliary imaging of the optical trap and devices in its immediate vicinity from both the yzy-z plane and the xyx-y plane, allows precise alignment of devices relative to the trapped microsphere, making use of a number of nanopositioning stages. The system can be utilized for force measurements in the range of 1μ1~{}\upmum-10μ10~{}\upmum, and for other applications requiring high precision, easily adjusted to trap larger or smaller microspheres.

Acknowledgements.
This work was supported by NSF grant PHY1802952, ONR grant N00014-18-1-2409, and the Heising-Simons Foundation. A.K. acknowledges the partial support of a William M. and Jane D. Fairbank Postdoctoral Fellowship of Stanford University. N.P. acknowledges the partial support of the Koret Foundation. We acknowledge regular discussions on the physics of trapped microspheres with the group of Prof. D.C. Moore at Yale. We also thank A.O. Ames, A.D. Rider, B. Sandoval, J. Singh, and T. Yu, who contributed to early developments of the experimental apparatus and the Physics Machine Shop at Stanford for their technical support. The data that support the findings of this study are available from the corresponding author upon reasonable request.

References