Optomechanical vector sensing of new forces at 6 micron separation

TL;DR

This paper improves the search for new gravity-like forces at very small scales using optically levitated microspheres, achieving much higher sensitivity than before.

Contribution

The study introduces multi-component force vector sensing and achieves a 100-fold sensitivity improvement in detecting new forces at micrometer separations.

Findings

New force strength is constrained to less than 10^7 at a Yukawa range of 5 micrometers.

Sensitivity improved by a factor of ~100 compared to prior measurements.

Results support efforts to study quantum gravity using micrometer-scale objects.

Abstract

The search for new gravity-like interactions at the sub-millimeter scale is a compelling area of research, with important implications for the understanding of classical gravity and its connections with quantum physics. We report improved constraints on Yukawa-type interactions in the \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$10\,\mathrm {\mu m}$$\end{document} regime using optically levitated dielectric microspheres as test masses. The search is performed, for the first time, sensing multiple spatial components of the force vector, and with sensitivity improved by a factor of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sim 100$$\end{document} with respect to previous measurements using the same technique. The resulting upper limit on the strength of a hypothetical new force is \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$10^7$$\end{document} at a Yukawa range \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\lambda \simeq 5\;\mu$$\end{document}m and close to \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$10^6$$\end{document} for \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\lambda \gtrsim 10\;\mu$$\end{document}m. This result also advances our efforts to measure gravitational effects using micrometer-size objects, with important implications for embryonic ideas to investigate the quantum nature of gravity.

Figures

Top: diagram of the experimental setup also defining the coordinate system used and providing the most important physical dimensions. A MS is held at the focus of the trapping beam while the density-patterned attractor oscillates by $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$170\,\mathrm {\mu m _{pp}}$$\end{document}$ at $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsf

Measurements of the *x* and *y* components of the force vector at the three harmonics with the largest observed power. The filled (unfilled) diamonds show the expected signal magnitudes and phases for an attractive (repulsive) Yukawa interaction with $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfo

Force spectra showing the measured backgrounds from three sources (from three different MSs). Top left: vibration backgrounds measured when the attractor is driven in a retracted position; top right: in situ measurement of electromagnetic backgrounds at the sidebands of the bias modulation frequency (note the different centering of the frequency scale). The measured forces have been projected along the phase of the signal template to reflect only the component limiting sensitivity. Bottom: scattered light backgrounds measured by the QPD null stream in the 2021 setup^[23](#CR23)^ and the curren

Constraints on the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha -\lambda$$\end{document}$ parameter space obtained in this work and the prior, 2021 result^[23](#CR23)^. The yellow shaded region indicates the region of parameter space excluded by other experiments^[4](#CR4),[12](#CR12),[14](#CR14),[15](#CR15)^.

Full text

Introduction

At large distances, gravity is the most familiar and best characterized of the four fundamental interactions. The motion of bodies in the Solar System has provided a platform for exquisite tests of Newtonian gravity^1^, while precision tests of General Relativity validate its use to describe phenomena also in the strong gravity regime^2,3^. On short length scales, however, the nature of the gravitational force is far less well explored. Indeed, it has never been directly measured for face-to-face separations between test masses smaller than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$52\,\mathrm {\mu m}$$\end{document} ^4^, and constraints on modifications to the inverse square law (ISL) weaken by orders of magnitude at distances shorter than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10~\mu$$\end{document} m^5^. Yet, the Standard Model (SM) of fundamental interactions implicitly assumes that the ISL applies unchanged all the way down to the Planck scale^6^. Models beyond the SM include a rich landscape of predictions that would manifest as modifications to gravity on short length scales, including the existence of extra dimensions or of new light bosons^7–11^. This broad class of new physics is typically parameterized by adding a Yukawa term to the gravitational potential between two point masses \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_1$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_2$$\end{document} ,

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} V(r) = -\frac{Gm_1m_2}{r}(1+\alpha e^{-r/\lambda }), \end{aligned}$$\end{document}

where G is Newton’s gravitational constant, r is the distance between the point masses, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} are the strength and length scale, respectively, of the Yukawa modification. The most stringent constraints on modifications to the ISL of this type have so far been established by measuring responses of mechanical oscillators, often torsion pendulums^4,12–15^. These techniques search for variations of one component of a temporally modulated force vector at one or more harmonics of the modulation frequency. In contrast, the technique presented here uses the time-dependent behavior of the full 3-dimensional force vector to search for new interactions. By conducting the search for a potential signal across multiple spatial dimensions and harmonics, the unique spectral fingerprint of a potential interaction can be exploited, especially to constrain and confirm a possible discovery.

Experimental setup

The sensor at the heart of this work is an optically levitated silica microsphere (MS)^16^ trapped in a single-beam, vertically oriented optical tweezer. An early generation of the setup is described in^17^. The experimental geometry is illustrated in Fig. 1. A \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1064\,\textrm{nm}$$\end{document} laser beam is focused with an off-axis parabolic mirror to a waist of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 3.5\,\mathrm {\mu m}$$\end{document} (Numerical Aperture \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{NA}\sim 0.1$$\end{document} , Rayleigh range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z_{\textrm{R}} \sim 30\,\mathrm {\mu m}$$\end{document} ). The MS is trapped near the focus, allowing the closest distance of approach of a source mass (“attractor”) with minimal interaction with the trapping beam. The stabilization of laser intensity and angular jitter by real-time feedback allows the MS to be isolated from extraneous disturbances at high vacuum ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 10^{-7}\,\textrm{hPa}$$\end{document} ). The light scattered off the MS in the forward direction is re-collimated by a second off-axis parabolic mirror and used to reconstruct the x and y positions of the MS^18^ using a quadrant photodiode (QPD). The phase of the retroreflected field from the MS is interferometrically sensed to reconstruct the z coordinate.

The trap is surrounded by six electrodes (three pairs, one for each Cartesian axis, separated by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 8\,\textrm{mm}$$\end{document} ). The electrodes are spray-coated with AquaDAG^19^ to reduce reflections of stray laser light. The monopole electric charge on the MS can be controlled with single-electron resolution^20,21^, allowing for a bias on the electrodes to exert a known force on the MS, thereby enabling the calibration of the optical position readout to force. The entire readout chain is found to be linear up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-13}\,\textrm{N}$$\end{document} ^21^, well beyond the forces of interest here.

A density-patterned attractor^22^ acts as a source for interactions coupling to mass. Gold and silicon are used as the two materials generating a contrast in density, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document} , with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho _{\textrm{Au}} \sim 8\rho _{\textrm{Si}}$$\end{document} . By scanning the position of the attractor along y in close proximity to the MS at a frequency \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_0=3~\textrm{Hz}$$\end{document} , the MS is subject to a characteristic force with a distinctive spectral fingerprint in each Cartesian axis that is unlikely to be mimicked by backgrounds. As electromagnetic forces constitute an important source of background, several design choices were made to mitigate their contribution. To prevent time-varying magnetic fields in the trap, the attractor is fabricated from materials that are not permanently magnetizable. The MS can be made to have zero overall charge, though dipole and higher order moments persist and can couple to time-varying electric field gradients as the attractor is scanned. In order to mitigate this effect, a stationary “shield” (a microfabricated silicon fence structure, conformally coated with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 100\,\textrm{nm}$$\end{document} gold) is positioned between the MS and the scanning attractor. In addition, using the six electrodes surrounding the trap, a rotating electric field can be applied to the dipole moment inducing the rotation of the MS, as explained in Section “Analysis”.

The first constraints on modifications to gravity using a levitated force sensor were reported in^23^. The present iteration of the experiment incorporates several improvements over the first. The amount of stray light scattering off the moving attractor and resulting in a background on the detectors was reduced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gtrsim 100-$$\end{document} fold. This was primarily achieved by adding a spatial filter at the focus of a downstream telescope, and by coating the attractor with a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 3\,\mathrm {\mu m}$$\end{document} -thick layer of Platinum Black^24^, with a reflectivity at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1064\,\textrm{nm}$$\end{document} that is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lesssim 1\%$$\end{document} that of the original gold. In addition, various parts of the vacuum chamber were stiffened to mitigate the effect of mechanical vibrations, and a series of accelerometers and microphones were added to sense the seismic and acoustic environment close to the setup. Optimal (Wiener) filters^25,26^ were then constructed to subtract environmental disturbances coherently sensed by the MS and the environmental sensors. Lastly, the noise in the force measurement was improved by over a factor of 10 along the x and y axes, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 5\times$$\end{document} along the z-axis. This improvement was achieved primarily with better stabilization of the laser intensity delivered to the trapping region (relative intensity noise, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{RIN}\lesssim 5\times 10^{-7}/\sqrt{\textrm{Hz}}$$\end{document} ). With these upgrades, the forward-scatter detection system has a displacement sensitivity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 10^{-11}~\mathrm {m/\sqrt{Hz}}$$\end{document} , which translates to a force sensitivity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 10^{-17}~\mathrm {N/\sqrt{Hz}}$$\end{document} for the typical trap spring constant of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 10^{-6}~\mathrm {N/m}$$\end{document} .Fig. 1. Top: diagram of the experimental setup also defining the coordinate system used and providing the most important physical dimensions. A MS is held at the focus of the trapping beam while the density-patterned attractor oscillates by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$170\,\mathrm {\mu m _{pp}}$$\end{document} at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_0=3\,\textrm{Hz}$$\end{document} behind the stationary shield. The inset at the top shows a side view with the MS diameter and the face-to-face separation between the MS and attractor identified. The Platinum Black coating on the attractor, which hides the underlying density pattern, has been omitted in this diagram. Bottom: time evolution of the force exerted on a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10~\mu \textrm{m}$$\end{document} MS by a hypothetical Yukawa interaction with the attractor ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =10^6$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda =10~\mu \textrm{m}$$\end{document} ). The z component of the force is not shown as it is only present when a vertical offset is introduced between the MS and the attractor. The secondary y axis shows the position of the center of the attractor as it is driven through one full cycle.

Analysis

The results presented in this work were obtained from measurements of three MSs, one with nominal diameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$7.56~\mu \textrm{m}$$\end{document} and the other two \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$9.98~\mu \textrm{m}$$\end{document} . Their masses are measured in situ, following the procedure in^27^. A detailed list of the experimental conditions under which data was collected is available in^28^. The search for new interactions is done by determining the value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} that provides the best fit between signal templates and the measured forces across multiple dimensions and harmonics. The fit includes both the magnitude of the force component and its phase with respect to the oscillation of the attractor. Signal templates are computed for various values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} using a finite-element model of the geometry of the attractor and MS over a spatial grid of the trap volume. An example template in the time domain is shown at the bottom of Fig. 1. The path traversed by the attractor during each measurement is recorded and used to construct the signal template for that particular trajectory. Both positive and negative values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} are considered, enabling a search for attractive and repulsive modifications to the Newtonian gravitational potential. The data presented here were collected with negligible offset between the attractor and MS centers-of-mass in the z direction, maximizing the signal along the x and y directions, at the expense of that in z.Fig. 2. Measurements of the x and y components of the force vector at the three harmonics with the largest observed power. The filled (unfilled) diamonds show the expected signal magnitudes and phases for an attractive (repulsive) Yukawa interaction with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} shown in the legend. The measurements are not consistent with either signal template across all harmonics.

While some power is present above the noise level at some of the harmonics of the 3 Hz scanning frequency, the relative distribution of such power does not match the Yukawa signal template and therefore the experiment found no evidence of a new interaction. To illustrate some features of the background, Fig. 2 shows the magnitude and phase of the effective forces measured for three harmonics with the largest power above the noise. The magnitudes and phases expected for true attractive and repulsive interactions parametrized by Eq. 1 are also shown, allowing for a direct comparison between the measurements and the signal template. Three sources of background are identified: mechanical vibrations, electromagnetic coupling, and scattering of stray light. All three of them substantially improved with respect to the previously reported results.Fig. 3. Force spectra showing the measured backgrounds from three sources (from three different MSs). Top left: vibration backgrounds measured when the attractor is driven in a retracted position; top right: in situ measurement of electromagnetic backgrounds at the sidebands of the bias modulation frequency (note the different centering of the frequency scale). The measured forces have been projected along the phase of the signal template to reflect only the component limiting sensitivity. Bottom: scattered light backgrounds measured by the QPD null stream in the 2021 setup^23^ and the current setup, after having added the spatial filter and coated the attractor with Platinum Black.

The mechanical coupling of vibrations from the oscillating attractor stage may cause relative motion between the MS, the beam, and the rest of the hardware. This contribution is measured by retracting the attractor along x by a large distance ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 5\,\textrm{mm}$$\end{document} ) from the nominal measurement position. The left panel of Fig. 3 shows the measured spectrum under these conditions, demonstrating that vibration backgrounds are present predominantly at the fundamental frequency. Excluding the fundamental frequency from the analysis renders this background negligible with minimal loss in signal power.

Interactions between the electric dipole moment (EDM) of the MS and electric field gradients^29^ from the attractor were already identified as another possible source of backgrounds. While the shield attenuates direct coupling, substantially better alignment and better control of the MS spinning have reduced this background to a sub-dominant strength. The effects of a contact potential difference between the attractor and shield could still produce field gradients at the location of the MS which would change as the relative position of the attractor and shield does. Measurements of the EDM of MSs using the technique described in^30^ give a typical EDM of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(100\, e\cdot \mu m)$$\end{document} .

A rotating ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_\textrm{rot}\sim 200\,\textrm{kHz}$$\end{document} ) electric field, applied using the electrodes surrounding the trap, is used to confine the EDM of the MS to a plane that minimizes backgrounds. Empirically, drifts in the out-of-plane component of the EDM are typically \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lesssim 10\%$$\end{document} over the course of a measurement. An initial investigation of the remaining electromagnetic background was carried out in situ by applying an amplitude modulation to the bias between the attractor and shield at a frequency \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_\textrm{bias}\gg f_0$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_\textrm{bias}\ll f_{\textrm{rot}}$$\end{document} ). Forces coupling to the EDM, which normally appear as backgrounds at the harmonics of the attractor’s mechanical drive \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n f_0$$\end{document} , would manifest as sidebands at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_\textrm{bias}\pm n f_0$$\end{document} . A Yukawa interaction is not expected to couple to the EDM, ensuring that these sidebands remain clean background witness channels, although potential effects of an induced dipole in the MS requires more investigation. The right panel of Fig. 3 shows the typical fingerprint obtained by this process (with a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$150\,\mathrm {mV_{pp}}$$\end{document} bias modulation and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_\textrm{bias}=40.5~\textrm{Hz}$$\end{document} ), indicating that EDM backgrounds appear to be sub-dominant with the current force sensitivity.

The dominant background currently originates from stray light from the trapping beam scattering off the scanning attractor and being imaged on the QPD. While great effort is taken to ensure good spatial quality of the trapping beam around the focus, ignoring the shield, it is estimated that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(100\,\textrm{ppm})$$\end{document} of light from the Gaussian ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {TEM_{00}}$$\end{document} ) trapping beam exists at the location of the attractor. In addition, it is plausible that a small fraction of light scatters off the AquaDAG-coated electrodes surrounding the trap volume and into the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x-y$$\end{document} (z) QPD (photodiode) detection systems, modulated by the reciprocating attractor. Investigations of this background have qualitatively confirmed several expected properties: (i) it is present in some form when the attractor scans with its nominal cycle but without a MS in the trap, (ii) it decreases as the attractor is retracted away from the MS, (iii) it increases as the attractor moves vertically away from the focus, cutting into more of the beam, and (iv) it was drastically reduced by the combination of using an attractor coated by absorbent Platinum Black and adding a spatial filter ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$50~\mu \textrm{m}$$\end{document} diameter Acktar^31^ coated aperture) to the downstream optics to select the mode containing information about MS motion. The first point allows for some aspects of this background to be studied in the absence of a MS. However, the Mie-scattering distribution from the MS significantly alters the light scattering pattern, preventing this technique from being used for quantitative background modeling. In addition, the background varies in time, making its effective subtraction more challenging. The QPD used to sense forward-scattered light from the MS consists of 4 separate sectors arranged in a square, from which four linearly-independent combinations can be constructed. Three of these are constructed to report the x and y positions of the MS, and the sum of light incident on the QPD. The fourth remaining linear combination is interpreted as a “null stream,” which is sensitive to random spatially fluctuating light fields, but empirically verified to be insensitive to a true force transduced by the MS. This construction of data streams from the QPD quadrants (described in detail in Sec 8.4.1 of^32^) was used only for estimating the contribution of stray light backgrounds. In the search for Yukawa interactions, the standard position data streams, given by Eq 8.1 of^32^ were used. As shown in the bottom panel of Fig. 3, the null stream shows substantial background reduction when a Platinum Black coated attractor is used in combination with the spatial filter, mitigating the effects of stray light.

As the measured forces are not consistent with the signature of a new Yukawa interaction, we report upper limits on the strength parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} for a range of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} values. A likelihood function is defined as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{aligned} \mathcal {L}(\alpha |\lambda ) = \prod _{i,j,k} \frac{1}{2\pi \sigma _{ijk}^2}&\times \exp \left\{ -\frac{\left[ \Re \left( F_{ijk} - \alpha \,\tau _{ijk}(\lambda )\right) \right] ^2}{2\sigma _{ijk}^2}\right\} \\&\times \exp \left\{ -\frac{\left[ \Im \left( F_{ijk} - \alpha \,\tau _{ijk}(\lambda )\right) \right] ^2}{2\sigma _{ijk}^2}\right\} , \end{aligned} \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{ijk}$$\end{document} is the discrete Fourier transform of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i^\textrm{th}$$\end{document} measurement at harmonic j of the force component along axis k, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau _{ijk}(\lambda )$$\end{document} is the corresponding signal template, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{ijk}$$\end{document} is the RMS noise for the relevant measurement, harmonic, and axis combination, measured with the attractor stationary at its nominal equilibrium position. Due to observed backgrounds which do not exhibit a spectral fingerprint consistent with Eq. 1, but which nonetheless vary significantly in their magnitude and phase, the choice of which harmonics to include in the analysis can affect the resulting limit. Including too many harmonics simply weakens the results by including channels with large observed backgrounds and little signal power, worsening the fit to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} . Using too few harmonics results in discarding a large fraction of the available signal power, and increases the sensitivity to the particular selection of harmonics used, and thus single-frequency systematic effects. To strike a balance, for each polarity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} , we use the combination of six harmonics that maximizes the goodness of fit, quantified by the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2$$\end{document} . This criteria conservatively selects the harmonics that most mimic the spectral behavior of a potential signal while ensuring that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$>50\%$$\end{document} of the spectral power of such a signal is captured. With the likelihood construction described above, the profile likelihood ratio reduces to a simple parabolic form,

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} -2\log {\frac{\mathcal {L}(\alpha |\lambda )}{\mathcal {L}(\hat{\alpha }|\lambda )}} = \frac{(\alpha - \hat{\alpha })^2}{\sigma _\textrm{stat}^2}, \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{\alpha }$$\end{document} is the maximum likelihood estimator for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _\textrm{stat}$$\end{document} describes the statistical uncertainty on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{\alpha }$$\end{document} . Constraints can then be placed on the strength parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} using a test statistic for upper limits, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_\alpha$$\end{document} , defined as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} q_{\alpha } = {\left\{ \begin{array}{ll} \frac{(\alpha - \hat{\alpha })^2}{\sigma _\alpha ^2} & |\alpha |\ge |\hat{\alpha }| \\ 0 & |\alpha |<|\hat{\alpha }| \end{array}\right. }, \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _\alpha ^2 = \sigma _\textrm{stat}^2 + \sigma _\textrm{sys}^2$$\end{document} is the total variance including both statistical and systematic uncertainties. The systematic uncertainties considered, along with their effect on the estimate of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} , are summarized in Table 1.Table 1. Systematic uncertainties in measured parameters and their contribution to the uncertainty on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda =10~\mathrm {\mu m}$$\end{document} .Parameter, p \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta p$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \alpha / \alpha$$\end{document} Transfer function amplitude10%0.1MS mass \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$20~\textrm{pg}$$\end{document} 0.04MS x position \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.7~\mu \textrm{m}$$\end{document} 0.1MS y position \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1~\mu \textrm{m}$$\end{document} 0.01MS z position \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.5~\mu \textrm{m}$$\end{document} 0.01Attractor thickness \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1~\mu \textrm{m}$$\end{document} 0.1Platinum Black thickness \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1~\mu \textrm{m}$$\end{document} 0.01Platinum Black density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10~\mathrm {g/cm^3}$$\end{document} 0.01

The 95% confidence level limit on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} for each sign is determined using the asymptotic behavior of the test statistic given by Wilks’ theorem^33^. This statistical procedure has been tested using synthetic datasets created by injecting (in post-processing) a signal into data collected with the attractor stationary, ensuring that a true signal can be reconstructed accurately. Unmodeled backgrounds present in both polarities across multiple harmonics can add to or subtract from a real signal in the data, potentially resulting in misreported statistical coverage. To ensure robustness against this, synthetic signals were combined with datasets containing the measured backgrounds. Multiple statistical frameworks were then used to reconstruct the injected signals, and the framework demonstrating the most accurate coverage was selected.

Results

Limits on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} for new attractive and repulsive interactions for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} ranging from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1-100~\mu \textrm{m}$$\end{document} were computed for each of the datasets^28^ considered. The strongest constraints among these are reported in Fig. 4. These new constraints represent an improvement of a factor of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 50$$\end{document} over previous results for length scales of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda>10~\mu \textrm{m}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gtrsim 100$$\end{document} at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda =2~\mu \textrm{m}$$\end{document} . Moreover, they are the first to use reconstruction of multiple spatial components of the force vector as a function of time, hence exploiting the unique signature of a new force and providing a robust pathway for discovery. For example, a z offset in either direction (c.f. the negligible offset used to collect data analyzed in this work) could be used to verify that a signal with the appropriate signature is transferred to this degree of freedom in case of a discovery.Fig. 4. Constraints on the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha -\lambda$$\end{document} parameter space obtained in this work and the prior, 2021 result^23^. The yellow shaded region indicates the region of parameter space excluded by other experiments^4,12,14,15^.

The sensitivity of this search is limited by backgrounds primarily due to stray laser light scattering off the attractor, creating modulations of the optical power on the QPD at harmonics of the attractor drive frequency. While substantial progress has been made in reducing this and other backgrounds, potential avenues for further improvement have been identified. Work is underway to implement a low-noise, high dynamic range and frame rate 100-pixel sensor, to be used for offline discrimination between true MS motion and scattered light effects. An alternative attractor geometry, consisting of a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$30\,\mathrm {\mu m}$$\end{document} thick, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.5\,\textrm{mm}$$\end{document} diameter rotating disk with azimuthal density modulations^34^, has also been prepared for potential use in the next experimental run. Better free space optics upstream of the trap is expected to produce more stable conditions and backgrounds, specifically addressing slow drifts of the trap focus which appear correlated with small variations in the refractive index of air driven by atmospheric pressure fluctuations. Beyond these immediate steps, a completely enclosed attractor should prevent any direct coupling between its motion and the beam, vastly attenuating optical backgrounds. Design of such an attractor system is a technical challenge to which future engineering efforts will be dedicated. Smaller and more constant stray light effects will facilitate the investigation of the remaining background, using some of the techniques described here.

Conclusion

In this paper we have described the use of optically levitated MSs as force sensors capable of probing new Yukawa interactions that couple to mass at micron-scale separations. In addition to complementing existing techniques, this platform uses the fully-reconstructed force vector to search for the complex temporal signature that would arise from such interactions, a feature that may become crucial should future work result in a discovery. Techniques to probe the main sources of background limiting the experiment, as well as mitigation strategies, have been identified. Apart from the direct physics goal of the experiment, this work may be relevant to other uses of levitated optomechanics for fundamental physics^35^, including searches for dark matter interactions^36,37^, sterile neutrino emission in weak decays^38^, or tests of the neutrality of matter^29^. In addition, this work demonstrates the operation of a setup with microscopic objects levitated a few micrometers away from mechanical structures, while maintaining stable conditions, exploring backgrounds, and performing precision metrology. These are among the challenges that will need to be addressed towards experimental tests of the quantum nature of gravity^39–41^.

Supplementary Information

Supplementary Information.