J Streque

M.-A Soussou

F Bartoli

S Hage-Ali

A Taguett

T Aubert

T Mazingue

M Lomello

O Elmazria

  are viewed as potential candidates for high temperature SAW applications, including langasite (LGS) and AlN. Only a few metals or alloys can withstand such temperatures, among them platinum and iridium based alloys or nanocomposites. These materials cannot be considered when high quality factors are required, due to their higher resistivity and density. The (Cr, V)-Si system was previously studied, but has never been proposed as a solution for SAW electrodes, in spite of its high thermal stability, lower density (about 6 g/cm 3 ) and decent resistivity. This study focuses on the development of Cr 5 Si 3 alloy deposition for its use in LGSbased resonators, aiming at wireless applications. Temperature compatible piezoelectric substrates -Langasite LGS (High stability perovskite, parabolic temperature behavior) -Aluminum nitride (Large bandgap, high resistivity semiconductor, linear TCF) Looking for trade-off for the choice of IDT electrodes -Low resistivity, decently stable with temperature -High thermal stability / Passivation -Low density, low hardness -Inexpensive at 500°C and 700°C Poster Presented at IUS 2018, Kobe, Japan Frequency behavior of LGS-based resonators: Cr 5 Si 3 vs Al electrodes Van der Pauw Resistivity on thin films of Cr 5 Si 3 Room Temperature and High Temperature behaviors of LGS resonators with Cr 5 Si 3 electrodes High temperature, harsh environmentsExploration of Cr 5 Si 3 for SAW resonator electrodes -Approach: reaching a good trade-off between temperature range, low resistivity, low density, thermal stability… -Cr 5 Si 3 shows good thermal stability up to 600°C -Resistivity remains high, but rather steady with temperatureFuture strategies for the improvement of Cr 5 Si 3 electrodes -Initial anneal of electrodes, for better thermal stability -Integration of a passivation layer -Study of Cr 5 Si 3 with substrates less prone to interdiffusion (AlN / Sapphire…) -Extension to Cr-V-Si ternary alloysFuture work on integration of Cr 5 Si 3 electrodes -Optimization of LGS-based and AlN-based resonators equipped with Cr 5 Si 3 electrodes -Further work on Cr 5 Si 3 aging behaviorOther study on AlN / Sapphire optimization for wireless HT SAW sensorsLower density than most metals / semi-metal alloys More resistive than Pt… …But more steady with temperature than Al Might be prone to diffusion… …But expected to have good thermal stability Trade-off to be investigated Cr 5 Si 3 alloy First investigations on bulk alloys -Investigated system: Cr-Si -Alloys prepared in electric arc furnace -Samples followed by Thermogravimetry Investigation of Cr 5 Si 3 for thin films -Prepared by DC sputtering from Cr-Si target -Optimization of parameters: -Power on target -Floating / Set potentials on substrate -Agence Nationale de la Recherche, ANR (SALSA Project) French PIA project "Lorraine Université d'Excellence" (ANR-15-IDEX-04-LUE)


Ying Tang

Seungkyu Ha

Thomas Begou

Julien Lumeau

Paul Urbach

Nynke H Dekker

J L Aurèle

Adam

H Paul Urbach

Aurele J L Adam

H Paul Urbach

Birefringent metamaterial nanocuboids with tailored optical constants for force and torque transduction

Keywords: metamaterials, dielectric multilayers, nanoparticles, optical trapping, force spectroscopy, torque spectroscopy

ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

T he development of optical trapping has led to numerous micro-and nanoscale applications over the past decades. [START_REF] Gao | Optical Manipulation from the Microscale to the Nanoscale: Fundamentals, Advances and Prospects[END_REF][START_REF] Polimeno | Optical Tweezers and Their Applications[END_REF][START_REF] Ding | Photonic Tractor Beams: A Review[END_REF] The ability to perform contact-free optical manipulation and accurate detection of linear and angular motion has enabled technologies such as force and torque spectroscopy of single biomolecules, [START_REF] Capitanio | Interrogating Biology with Force: Single Molecule High-Resolution Measurements with Optical Tweezers[END_REF][START_REF] Forth | Torque Measurement at the Single-Molecule Level[END_REF][START_REF] Lipfert | Torque Spectroscopy for the Study of Rotary Motion in Biological Systems[END_REF][START_REF] Van Oene | Quantifying the Precision of Single-Molecule Torque and Twist Measurements Using Allan Variance[END_REF] quantum optomechanics, [START_REF] Arita | Laser-Induced Rotation and Cooling of a Trapped Microgyroscope in Vacuum[END_REF][START_REF] Ahn | Optically Levitated Nanodumbbell Torsion Balance and GHz Nanomechanical Rotor[END_REF] colloidal self-assembly, [START_REF] Nan | Silver-Nanowire-Based Interferometric Optical Tweezers for Enhanced Optical Trapping and Binding of Nanoparticles[END_REF] and optically driven nanomachines. [START_REF] Williams | Transmission of Torque at the Nanoscale[END_REF] For such applications, dielectric materials, e.g., silica and polystyrene (PS), have been popular choices as force transducers. [START_REF] Capitanio | Interrogating Biology with Force: Single Molecule High-Resolution Measurements with Optical Tweezers[END_REF] For the combined manipulation of force and torque, optically birefringent dielectric probes have proven their unique advantages over metallic probes, [START_REF] Shao | Gold Nanorod Rotary Motors Driven by Resonant Light Scattering[END_REF][START_REF] Shen | Trapping and Rotating of a Metallic Particle Trimer with Optical Vortex[END_REF] enabling reduced scattering, less heating, straightforward torque detection via direct measurement of spin angular momentum transfer, and fully confined rotational degrees of freedom. [START_REF] Lipfert | Torque Spectroscopy for the Study of Rotary Motion in Biological Systems[END_REF] Various uniaxial birefringent crystals (Figure 1a) have been successfully applied as optical force and torque transducers, including quartz SiO 2 , [START_REF] Deufel | Nanofabricated Quartz Cylinders for Angular Trapping: DNA Supercoiling Torque Detection[END_REF][START_REF] Gutieŕrez-Medina | An Optical Apparatus for Rotation and Trapping[END_REF][START_REF] Huang | Electron Beam Fabrication of Birefringent Microcylinders[END_REF][START_REF] Li | Fabrication of Birefringent Nanocylinders for Single-Molecule Force and Torque Measurement[END_REF][START_REF] Santybayeva | Fabrication of Quartz Microcylinders by Laser Interference Lithography for Angular Optical Tweezers[END_REF] calcite CaCO 3 , [START_REF] Friese | Optical Alignment and Spinning of Laser-Trapped Microscopic Particles[END_REF] vaterite CaCO 3 , 8,20 liquid crystal RM257, [START_REF] Ramaiya | Kinesin Rotates Unidirectionally and Generates Torque While Walking on Microtubules[END_REF] and rutile TiO 2 . [START_REF] Van Oene | Quantifying the Precision of Single-Molecule Torque and Twist Measurements Using Allan Variance[END_REF][START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF][START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] Nonetheless, all these dielectrics come with particular limitations. Silica and PS are commercially available as highly uniform and chemically stable microspheres with a wide range of diameters and surface chemistries. However, they are not suitable as spin angular momentum-based torque transducers due to their absence of optical birefringence. Another chemically stable material, quartz SiO 2 , is optically birefringent and thus permits straightforward transfer and detection of torque. Moreover, the ability to produce rod-shaped birefringent quartz SiO 2 particles using top-down fabrication [START_REF] Deufel | Nanofabricated Quartz Cylinders for Angular Trapping: DNA Supercoiling Torque Detection[END_REF][START_REF] Gutieŕrez-Medina | An Optical Apparatus for Rotation and Trapping[END_REF][START_REF] Huang | Electron Beam Fabrication of Birefringent Microcylinders[END_REF][START_REF] Li | Fabrication of Birefringent Nanocylinders for Single-Molecule Force and Torque Measurement[END_REF][START_REF] Santybayeva | Fabrication of Quartz Microcylinders by Laser Interference Lithography for Angular Optical Tweezers[END_REF] is highly beneficial for torsional control, as such particle geometry allows for fully constrained rotational

This is an open access article published under a Creative Commons Non-Commercial No Derivative Works (CC-BY-NC-ND) Attribution License, which permits copying and redistribution of the article, and creation of adaptations, all for non-commercial purposes. degrees of freedom in three dimensions by extra lateral torques. [START_REF] Singer | Orientation of Optically Trapped Nonspherical Birefringent Particles[END_REF][START_REF] Bareil | Angular and Position Stability of a Nanorod Trapped in an Optical Tweezers[END_REF] However, its very low birefringence (Δn = 0.009) [START_REF] Ghosh | Dispersion-Equation Coefficients for the Refractive Index and Birefringence of Calcite and Quartz Crystals[END_REF] severely limits the maximum torque that can be transferred to it. Although the torque transfer efficiency can be maximized by enlarging the probe, this not only decreases the trapping stability in single-beam optical traps but also reduces the maximum rotational speed and spatiotemporal resolution due to the increased hydrodynamic drag. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] Although calcite/ vaterite CaCO 3 crystals have higher birefringence (Δn = 0.1-0.16), [START_REF] Ghosh | Dispersion-Equation Coefficients for the Refractive Index and Birefringence of Calcite and Quartz Crystals[END_REF][START_REF] Christy | A Review of the Structures of Vaterite: The Impossible, the Possible, and the Likely[END_REF] they are unstable in water and therefore less suited to biological and microfluidic environments. In addition, the possible particle shapes have been limited to irregular fragments for the case of calcite CaCO 3 [START_REF] Friese | Optical Alignment and Spinning of Laser-Trapped Microscopic Particles[END_REF] and chemically synthesized spheres for vaterite CaCO 3 . [START_REF] Arita | Laser-Induced Rotation and Cooling of a Trapped Microgyroscope in Vacuum[END_REF][START_REF] Arita | Rotational Dynamics and Heating of Trapped Nanovaterite Particles[END_REF] Liquid crystal RM257 is stable in aqueous solutions while having similarly high birefringence (Δn = 0.18). [START_REF] Ramaiya | Kinesin Rotates Unidirectionally and Generates Torque While Walking on Microtubules[END_REF] However, as in the case of vaterite CaCO 3 , the possible particle shape of RM257 is limited to spheres. [START_REF] Ramaiya | Kinesin Rotates Unidirectionally and Generates Torque While Walking on Microtubules[END_REF] Meanwhile, another chemically stable material, rutile TiO 2 , can be fabricated in top-down fashion into rod-shape particles. [START_REF] Van Oene | Quantifying the Precision of Single-Molecule Torque and Twist Measurements Using Allan Variance[END_REF][START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF][START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] Rutile TiO 2 has the highest birefringence among all the known naturally occurring crystalline materials (Δn = 0.26), [START_REF] Devore | Refractive Indices of Rutile and Sphalerite[END_REF] and its average refractive index (n ̅ = 2.6) [START_REF] Devore | Refractive Indices of Rutile and Sphalerite[END_REF] also exceeds those of other common dielectric probes (n ̅ = 1.4-1.6). [START_REF] Malitson | Interspecimen Comparison of the Refractive Index of Fused Silica[END_REF][START_REF] Sultanova | Dispersion Properties of Optical Polymers[END_REF] As we have demonstrated in ref 23, this combination of high optical index and high birefringence makes it possible, even at only moderate beam powers, to achieve high levels of force and torque transfer efficiencies. Despite its superior properties, rutile TiO 2 does not provide a universal solution for all applications. For example, while the high refractive index of rutile TiO 2 makes it possible to generate high linear stiffness and force, it also results in a large scattering force. To compensate this high scattering force with optical gradient force, a tightly focused trapping beam with a high numerical aperture (NA) objective lens is therefore necessary. Even then, the probe size should be limited to guarantee stable threedimensional (3D) trapping. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] These stringent requirements hamper its wide usage, such as in compact and affordable optical trapping systems that often possess low NA optics. [START_REF] Stallinga | Optical Trapping at Low Numerical Aperture[END_REF][START_REF] Hwu | Hacking CD/DVD/Blu-Ray for Biosensing[END_REF] Furthermore, as the 3D-trapping stability of rutile TiO 2 probes is highly size-and shape-dependent, their use in applications regarding larger sizes or exotic shapes, such as gear-shaped microfluidic rotors, [START_REF] Higurashi | Optically Induced Rotation of a Trapped Micro-Object About an Axis Perpendicular to the Laser Beam Axis[END_REF][START_REF] Neale | All-Optical Control of Microfluidic Components Using Form Birefringence[END_REF] would be impractical. Moreover, high trap stiffnesses are not always desirable. Lower trap stiffnesses are intentionally used in certain applications, e.g., photonic force microscopy, where a "soft" probe is beneficial for scanning fragile samples such as cell surfaces. [START_REF] Friedrich | Surface Imaging beyond the Diffraction Limit with Optically Trapped Spheres[END_REF] Therefore, as the optical properties of the currently available birefringent probes do not cover the full application space of optical trapping, there exists a need for a joint force and torque transducer that possesses designed optical constants. Such designed particles could allow one to have any combination of refractive index and birefringence that provides the optimal balance between force and torque transducer properties per (blue) with its long axis aligned to the propagation direction (z-axis) of the trapping beam (red). Given an input beam linearly polarized along the x-axis, the particle's optic axis (n ⊥ ) aligns with the y-axis, while the axes associated with the higher refractive index (n ∥ ) align with the xand z-axes. The resulting optical birefringence in the x-y plane enables active rotary control around the z-axis (indicated by the curved arrow). The cross section in the x-y plane is a square of width W, and the aspect ratio AR sets the height H along the z-axis. (e) Scanning electron microscopy (SEM) image of MM particles with ρ = 0.3, W = 400 nm, and AR = 5. The three differently oriented particles are marked: bottom (1), side (2), and top (3) (as defined in Figure S10). specific application. In addition, they should be chemically stable for demanding applications in aqueous environments. Furthermore, their fabrication should accommodate high uniformity and yield, as well as flexibility in particle geometry. For the previously developed top-down fabrication of quartz SiO 2 and rutile TiO 2 nanoparticles, [START_REF] Deufel | Nanofabricated Quartz Cylinders for Angular Trapping: DNA Supercoiling Torque Detection[END_REF][START_REF] Gutieŕrez-Medina | An Optical Apparatus for Rotation and Trapping[END_REF][START_REF] Huang | Electron Beam Fabrication of Birefringent Microcylinders[END_REF][START_REF] Li | Fabrication of Birefringent Nanocylinders for Single-Molecule Force and Torque Measurement[END_REF][START_REF] Santybayeva | Fabrication of Quartz Microcylinders by Laser Interference Lithography for Angular Optical Tweezers[END_REF][START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF][START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] a crystalline substrate was etched vertically into rod-shaped pillars, which were converted into particles by physically cleaving the base of the pillars. This approach comes with several drawbacks: first, the use of crystalline substrates hampers high volume production, as they are often expensive and only come in small sizes; second, using vertical etching it is challenging to achieve a wide range of aspect ratios for rod-shaped particles; [START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF] and third, physical cleavage may result in unexpected damage to the particles, degrading their uniformity and reducing yield.

We address all the above-mentioned issues by switching to designed metamaterial particles that can be fabricated in a straightforward manner using robust processes. The choice of metamaterials is motivated firstly by the possibility of tailoring their physical properties beyond the range that exists in naturally occurring materials. For example, metamaterials with artificially defined values of the refractive index [START_REF] Padilla | Negative Refractive Index Metamaterials[END_REF][START_REF] Choi | A Terahertz Metamaterial with Unnaturally High Refractive Index[END_REF] and birefringence [START_REF] Min | Dual Metamaterial with Large Birefringence[END_REF] have been demonstrated. Within this context, multilayer metamaterials (MMs) provide a particularly flexible, straightforward means to design and implement desired values of optical constants, [START_REF] Guo | Applications of Hyperbolic Metamaterial Substrates[END_REF][START_REF] Jahani | All-Dielectric Metamaterials[END_REF] as the optical properties of MM particles can be customized by pairing materials from the extensive library of diverse amorphous dielectrics and by the layer thickness ratio. As the latter can be tuned in a continuous manner, this allows us to design particles with a nearly arbitrary combination of refractive index and optical birefringence. The choice of metamaterials is motivated secondly by the simplicity of fabrication of MM particles. As we demonstrate below, the fabrication process is robust, compatible with conventional semiconductor fabrication processes, and capable of largevolume production with high uniformity. A multilayer is achieved by sequentially stacking different material layers on a silicon substrate using various thin-film deposition techniques such as sputtering, [START_REF] Lang | Gold-Silicon Metamaterial with Hyperbolic Transition in near Infrared[END_REF] evaporation, [START_REF] Kruk | Magnetic Hyperbolic Optical Metamaterials[END_REF] and atomic layer deposition (ALD). [START_REF] Kelly | Designing Optical Metamaterial with Hyperbolic Dispersion Based on an Al:ZnO/ ZnO Nano-Layered Structure Using the Atomic Layer Deposition Technique[END_REF] Then, the multilayer undergoes lithography and etching to define the dimensions of nanoparticles, which are then harvested in their free-floating form by dissolving the underlying sacrificial layer. Such sacrificial layer-based harvesting is in principle free from undesired deformations in particle geometry, [START_REF] Higurashi | Optically Induced Rotation of a Trapped Micro-Object About an Axis Perpendicular to the Laser Beam Axis[END_REF][START_REF] Neale | All-Optical Control of Microfluidic Components Using Form Birefringence[END_REF][START_REF] Hernandez | Colloidal Alphabet Soup: Monodisperse Dispersions of Shape-Designed Lithoparticles[END_REF][START_REF] Clerico | Water-Dispersible Three-Dimensional LC-Nanoresonators[END_REF][START_REF] Badman | High Trap Stiffness Microcylinders for Nanophotonic Trapping[END_REF] in contrast to the m e c h a n i c a l c l e a v a g e -b a s e d p r o t o c o l s de s c r i b e d above. [START_REF] Deufel | Nanofabricated Quartz Cylinders for Angular Trapping: DNA Supercoiling Torque Detection[END_REF][START_REF] Gutieŕrez-Medina | An Optical Apparatus for Rotation and Trapping[END_REF][START_REF] Huang | Electron Beam Fabrication of Birefringent Microcylinders[END_REF][START_REF] Li | Fabrication of Birefringent Nanocylinders for Single-Molecule Force and Torque Measurement[END_REF][START_REF] Santybayeva | Fabrication of Quartz Microcylinders by Laser Interference Lithography for Angular Optical Tweezers[END_REF][START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF][START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] As the geometry along the long axis of a rod-shaped particle can be lithographically defined and laterally etched, higher aspect ratios can be readily obtained, without the need for thicker or harder etch masks. [START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF] This also allows for particles with geometries other than rod-shaped, and even includes exotic shapes. [START_REF] Higurashi | Optically Induced Rotation of a Trapped Micro-Object About an Axis Perpendicular to the Laser Beam Axis[END_REF][START_REF] Neale | All-Optical Control of Microfluidic Components Using Form Birefringence[END_REF][START_REF] Hernandez | Colloidal Alphabet Soup: Monodisperse Dispersions of Shape-Designed Lithoparticles[END_REF] It goes without saying that the use of large-area and cost-effective standard silicon substrates for our MM particles is better suited for mass production compared to natural crystalline substrates.

Using these ingredients, we here demonstrate the realization of birefringent MM particles with tailored optical constants. As a proof of concept, we have designed MMs consisting of niobium pentoxide (Nb 2 O 5 ) and silicon dioxide (SiO 2 ) that have a large birefringence at only moderate refractive index. The optical constants of the resulting MM particles present all the advantages of calcite/vaterite CaCO 3 crystals, while overcoming their known disadvantages such as chemical instability in aqueous environments and the lack of fabrication options for rod-shaped particles. We show using optical trapping experiments that one can achieve highly efficient torque transfer with such 3D-trapped MM particles dispersed in water and compare the experimental performance with theoretical calculations. The utilization of designed MM particles is further facilitated by our robust, sacrificial layerbased top-down fabrication route, which allows high yield, uniformity, and flexibility in particle geometry. Such universal birefringent MM particles, which can be designed either to emulate the properties of existing crystals or to introduce entirely novel ones, would further broaden the applicability of force and torque transducers in optical traps.

RESULTS AND DISCUSSION

Design of Birefringent MM Particles. The design of the artificially birefringent MMs is based on effective medium theory (EMT). [START_REF] Cai | Optical Metamaterials[END_REF] EMT predicts that a multilayer of alternating isotropic dielectric materials exhibits an effective birefringence. When the thickness of each layer is much smaller than the wavelength of the incident electromagnetic wave, the effective optical permittivity of the structure can be approximated by a tensor:

ϵ̲ = ϵ ϵ ϵ ⊥ Ä Ç Å Å Å Å Å Å Å Å Å Å Å Å Å Å Å Å Å É Ö Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ Ñ 0 0 0 0 0 0 (1)
where ϵ ∥ and ϵ ⊥ are the permittivity components parallel and perpendicular to the interfacial surfaces of the multilayer stack, respectively. Their values are given by

ρ ρ ϵ = ϵ + -ϵ (1 ) 1 2
(2)

ρ ρ ϵ = ϵ + -ϵ ⊥ - - 1/( (1 ) ) 1 1 2 1 (3)
where ϵ 1 and ϵ 2 are the permittivity values of the higher-and lower-index materials, respectively, and the material filling ratio ρ = d 1 /(d 1 + d 2 ) is defined as the thickness ratio of the higherindex material layer (d 1 ) and the unit layer-pair (d 1 + d 2 ). As a primary method to validate the accuracy of the EMT approximation for the optical properties of MM, we compare the Fresnel coefficients calculated from the approximated homogeneous single-layered structure, whose effective optical constants are calculated using EMT (eqs 1-3), with those from the actual multilayered structure [START_REF] El Gawhary | Extension of the Classical Fabry-Perot Formula to 1D Multilayered Structures[END_REF] (Figure S1). We find that the Fresnel coefficients calculated from the effective single-layer geometry approximate those of the multilayer geometry for both p-and s-polarizations of the incident beam (Figure S1) provided that the number of layer-pairs exceeds 2, confirming that EMT is sufficiently precise to describe our MMs. We additionally validated the EMT prediction by calculating the optical trapping forces and torques using finite element methods (FEM) (Figure S2; Methods).

For the constituent materials of MMs, we focus on dielectric materials that are commonly used in thin-film deposition techniques (e.g., evaporation, sputtering, and ALD) and able to form stable multilayer compositions. For most of such dielectric materials, there are well-developed surface functionalization protocols [START_REF] Sperling | Surface Modification, Functionalization and Bioconjugation of Colloidal Inorganic Nanoparticles[END_REF][START_REF] Pujari | Covalent Surface Modification of Oxide Surfaces[END_REF][START_REF] Reimhult | Design of Surface Modifications for Nanoscale Sensor Applications[END_REF] that could facilitate the use of MMs in diverse areas including biological applications. As a proof of concept, we searched for promising combinations of high-and low-index material pairs for MMs, with the goal of obtaining an

ACS Nano

MM with a large birefringence at only moderate refractive index, a combination of optical properties also found in CaCO 3 crystals (Figure S3). We examined several high-index materials as potential partners for the low-index material SiO 2 and ultimately chose Nb 2 O 5 . The resulting Nb 2 O 5 /SiO 2 multilayer emulates the optical constants of both calcite and vaterite CaCO 3 crystals (Figure S3) and comes with reliable fabrication (as demonstrated in other optical applications [START_REF] Lumeau | Angularly Tunable Bandpass Filter: Design, Fabrication, and Characterization[END_REF] ). Using plasma-assisted reactive magnetron sputtering, we are able to deposit Nb 2 O 5 /SiO 2 multilayers with precisely controlled thickness and refractive indices (n = 1.47 for SiO 2 and n = 2.26 for Nb 2 O 5 at 1064 nm) (Methods).

Using EMT, the optical properties of this multilayer combination can be predicted as a function of material filling ratio ρ (eqs 1-3). One can observe in Figure 1b ), indicating that the Nb 2 O 5 /SiO 2 multilayer emulates a negative uniaxial birefringence (n e < n o ). From these indices, the multilayer's full tunable range for the magnitude of the effective birefringence (Δn = |n en o |) and the corresponding average refractive index (n ̅ = (n e + n o )/2) can be obtained. The value of the birefringence (Δn) can be fixed anywhere between zero and its maximum value (Δn = 0.17 when ρ = 0.55), while the value of index (n ̅ ) varies between 1.47 (that of SiO 2 ) and 2.26 (that of Nb 2 O 5 ). Except at the maximum birefringence, there are two distinct filling ratios (and hence two different average indices) that achieve the same degree of birefringence (Figure 1b). This allows for a trade-off between 3D-trapping stability and maximum achievable force: a smaller index renders the MM particle more easily trapped in 3D due to the reduced scattering force, while a larger index enhances the gradient force due to the increased index contrast with the surrounding medium (i.e., water, in our experiments). Depending on the target application or experimental situation, one can determine the suitable filling ratio to obtain a desired combination of refractive index and birefringence for MM particles.

We have selected two different filling ratios for Nb 2 O 5 /SiO 2 multilayers (Figure 1b) to demonstrate the tunability of the optical constants in our MMs and thus the ability to emulate diverse birefringent crystals. We first select ρ = 0.5 (close to the value of ρ = 0.55 that maximizes birefringence) to demonstrate a birefringence (Δn = 0.16) as large as that of calcite CaCO 3 . We also select ρ = 0.3 with the goal of achieving a lower refractive index (n ̅ = 1.68, compared to n ̅ = 1.82 at ρ = 0.5) while maintaining a high level of birefringence (Δn = 0.13). Further decrease of the material filling ratio to ρ = 0.2 should result in lower values for the index (n ̅ = 1.61) and birefringence (Δn = 0.10) that are comparable to those of 1d. All trapping parameters presented are calculated for the case of an input beam that is linearly polarized along the x-axis and the specific configuration of the employed OTW setup (Methods) and then normalized by the theoretically calculated input laser beam power (Methods). In all maps, the solid blue lines represent contour lines for nanocuboids of constant height that range from 500 nm (leftmost) to 3000 nm (rightmost) at intervals of 500 nm, and the pixel size is ΔAR = 0.1, ΔW = 5 nm. The particle dimensions that correspond to the maximum for each map are indicated by the empty magenta circles. The black pixels in the maps of ρ = 0.5 (panels (b, d, f, h)) indicate the particle dimensions that cannot be trapped in 3D due to excessive scattering forces. The 10 different particle dimensions (A3-5, B3-5, and C3-6) chosen for demonstration are indicated by the crossing-points of the black dash-dotted lines, overlaid with corresponding symbols in each map, and their labels are displayed in panels (a) and (b).

vaterite CaCO 3 . We note that for Nb 2 O 5 /SiO 2 MMs, values of ρ < 0.55 realize the combination of moderate index and high birefringence, similar to calcite/vaterite CaCO 3 crystals (n ̅ ≈ 1.6; Δn = 0.1-0.16); [START_REF] Ghosh | Dispersion-Equation Coefficients for the Refractive Index and Birefringence of Calcite and Quartz Crystals[END_REF][START_REF] Christy | A Review of the Structures of Vaterite: The Impossible, the Possible, and the Likely[END_REF] if higher values of the linear trapping stiffness are desired for the same degree of birefringence, then a higher filling ratio (ρ > 0.55) should be employed (Figure 1b). As can be seen in Figure 1c, numerous other material combinations can also be employed to fabricate MMs with properties either similar to those of existing crystals or going beyond them. For example, the predicted optical constants of MMs composed of aluminum oxide (Al 2 O 3 ) and SiO 2 multilayer pairs can mimic those of quartz SiO 2 . Meanwhile, MMs composed of amorphous silicon (a-Si) and SiO 2 multilayer pairs could either emulate the optical properties of rutile TiO 2 (the most birefringent crystal found in nature, with Δn = 0.26) or access even 3-fold higher values of the birefringence (Δn = 0.86).

Based on this multilayer structure, the MM particle geometry was designed. To ensure full three-dimensional rotational confinement of each MM particle, the particle is rodshaped with a high aspect ratio (AR). The high AR generates sufficiently large lateral torques to maintain the particle's orientation with its long axis aligned to the direction of laser beam propagation (z-axis) in an optical trap as shown in Figure 1d, unless their sizes are much smaller than the waist of the focused beam. [START_REF] Singer | Orientation of Optically Trapped Nonspherical Birefringent Particles[END_REF][START_REF] Bareil | Angular and Position Stability of a Nanorod Trapped in an Optical Tweezers[END_REF] As the particle's optic axis (n ⊥ ) is perpendicular to its longer side (Figure 1d), active rotary control of the particle around the z-axis can then be achieved by exploiting its optical birefringence in the x-y plane and rotating the linear polarization of the trapping beam. [START_REF] Lipfert | Torque Spectroscopy for the Study of Rotary Motion in Biological Systems[END_REF] In our top-down fabrication (Figure 1e, Methods), we have confined our particle shape to be square-cuboid to assess only the torque induced by optical birefringence. While the cross section of a square-cuboid in the x-y plane is not fully rotationally symmetric (unlike that of a cylinder [START_REF] Huang | Electron Beam Fabrication of Birefringent Microcylinders[END_REF][START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF] ), compared to a rectangular-cuboid or other possible geometries it is relatively free of unwanted geometrical torques that would affect rotation about the z-axis and hinder precise angular control and detection.

Numerically Calculated Trapping Performance of MM Particles. We have performed FEM calculations of the linear and angular optical trapping landscapes (Methods) to underpin our design of MM particles. Here, we demonstrate the results for the case of Nb 2 O 5 /SiO 2 MM particles, with two chosen material filling ratios (ρ = 0.3 and ρ = 0.5; Figure 1b). We vary the width W of nanocuboids over a range of 250-450 nm, while we examine AR ranging from 1 (cubes) to 7 (square cuboids with increasing heights).

To predict the 3D-trapping stability for MM nanocuboid particles, we have calculated the axial stiffness κ z for ρ = 0.3 and ρ = 0.5 and plotted this in a heat map (Figure 2a,b). For each particle dimension, κ z is defined as the slope of the axial force curve (Figure S2) at the equilibrium trapping position z eq (Figure S4). The nanocuboids with ρ = 0.3 are trappable over the entire range of widths and aspect ratios in the map. In contrast, the nanocuboids with ρ = 0.5 include nontrappable dimensions which appear at W ≳ 375 nm and are displayed as black pixels in Figure 2b; in such cases, the axial stiffnesses cannot be defined, as the dominant scattering forces result in an absence of z eq . Calculations at the corresponding dimensions were not carried through and are thus also displayed as black pixels in Figure 2d,f,h. Additional calculations in which W is extended up to 550 nm (Figure S5) show that the nanocuboids with ρ = 0.3 become nontrappable at W ≳ 525 nm. The much narrower range of 3D-trappable dimensions for MM particles with ρ = 0.5 (also reflected by the occurrence of maximal stiffness values at far smaller widths; compare the positions of the magenta circles in Figure 2a,b) results from the larger scattering forces that derive from their higher average refractive index (n ̅ = 1.82, compared to n ̅ = 1.68 for ρ = 0.3).

We also calculated the lateral stiffness κ y (Figure 2c,d) at z eq , which generally exhibits larger values than κ x when the input beam is linearly polarized along the x-axis and a subwavelength-sized particle is trapped. [START_REF] Madadi | Polarization-Induced Stiffness Asymmetry of Optical Tweezers[END_REF] Overall, the lateral stiffness κ y is almost an order of magnitude larger than the axial stiffness κ z and maximized around W ≈ 300 nm (magenta circles in Figure 2c,d), regardless of the filling ratio. We attribute such a common optimal particle width (∼300 nm) for the maximal lateral forces to an ideal degree of overlap between particle and field gradient profile of the focused beam in the lateral direction.

The calculated maps of angular stiffness κ θ (Figure 2e,f) at z eq represent the maximum torque (τ o = κ θ /2) experienced by each particle rotated 45°with respect to the direction of the input beam linear polarization. [START_REF] Santybayeva | Optical Torque Wrench Design and Calibration[END_REF] Overall, within the range of dimensions shown in the maps, a larger angular stiffness can be achieved with particles with higher aspect ratios for a given width, which can be explained by the dependence of torque transfer efficiency on particle height. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] When we compare the values of the angular stiffness for particles with ρ = 0.5 versus ρ = 0.3 (with otherwise identical dimensions), we find that the maximum angular stiffness in each map (magenta circles in Figure 2e,f) is larger for particles with ρ = 0.5 (877 pN•nm/ rad/mW) versus ρ = 0.3 (653 pN•nm/rad/mW). This difference results from their higher birefringence (by 27%) and ensures a higher torque transfer efficiency for these particles.

Although the angular stiffness of the ρ = 0.5 nanocuboids can be maximized by selecting particle dimensions in the finger-shaped regions [START_REF] Muradoglu | Optical Trapping Map of Dielectric Spheres[END_REF] at the top of the map (e.g., magenta circle in Figure 2f), in practice their utilization in 3D trapping is difficult. Due to variation in particle dimensions within a batch, some fraction of particles will lie outside of the fingershaped region, adopting dimensions that do not support 3D trapping. Also, the nonidealities of the actual trapping system (e.g., aberrations of the trapping beam) will modify the trapping landscape (e.g., shift the finger-shaped regions or threshold width for 3D trapping); [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] hence, even particles fabricated to high uniformity might have reduced stiffness or become nontrappable. In contrast, due to their lower refractive index, nanocuboids with ρ = 0.3 are predicted to be trappable for all dimensions shown in the map, including those with the dimensions that maximize angular stiffness (e.g., magenta circle in Figure 2e).

Using these FEM-calculated maps of angular stiffness κ θ (Figure 2e,f), together with the map for the rotational drag coefficient γ θ (Figure S6), we generate maps for the predicted maximum angular speed f o (= (κ θ /2)/(2πγ θ )) of the nanocuboids (Figure 2g,h). Overall, higher angular speeds are achieved for particles with smaller widths at fixed heights (and thus similar torque transfer efficiencies), which makes sense as the rotation of larger particles is hindered by increased rotational drag in the surrounding medium (water). The angular speeds of nanocuboids with ρ = 0.5 are higher compared to similarly sized nanocuboids with ρ = 0.3

ACS Nano

(compare the magenta circles at 61 Hz/mW and 45 Hz/mW in Figure 2h,g, respectively), reflecting the enhanced angular stiffnesses of particles with ρ = 0.5.

On the basis of these results, we selected our MM particle dimensions for each filling ratio and designed three different multilayer wafers for their fabrication. We deposited two multilayers with ρ = 0.3, with widths W of 300 nm (wafer A) and 400 nm (wafer B), respectively, and one multilayer with ρ = 0.5, with a width W of 300 nm (wafer C) (Table S1). The smaller width W = 300 nm is selected because for both filling ratios it combines stable axial trapping (Figure 2a,b) with maximal lateral trapping (Figure 2c,d). It thus allows for an examination of trapping behavior as a function of filling ratio for identically sized, 3D-trapped particles. For example, ρ = 0.5 particles, having higher birefringence than ρ = 0.3 particles, are expected to achieve both higher torques and angular speeds at W = 300 nm (Figure 2e,f,g,h). To explore the trapping behavior of particles with much larger volumes, we also selected the larger width W = 400 nm for ρ = 0.3 (Figure 2a) (particles of this width not being stably 3D-trappable for ρ = 0.5; Figure 2b). Particles with this width are predicted to exhibit larger torques but slower angular speeds than particles with W = 300 nm (Figure 2e,f,g,h). For all three wafers, we fixed the thickness of the unit layer-pair to 50 nm to enable the use of EMT in predicting the effective optical constants for the designed multilayer particle widths by allowing a sufficiently large number of layer-pairs (6 or 8, Figure S1) and to provide reliable control of the selected filling ratios by defining the thinnest layer to exceed ∼10 nm (Table S1), large compared to the nanometer-scale precision of our multilayer deposition (Methods). Regarding the particle height, we fabricated particles with AR of 3, 4, and 5 for all three wafers, where the range is selected to probe the height-dependence of the torque transfer efficiency while ensuring full three-dimensional rotational confinement. The resulting particle heights span the range of H ≈ 1-2 μm and thus include heights predicted to achieve torque transfer efficiencies close to 100%, i.e., H = 2.0 μm for ρ = 0.3 (W = 400 nm, AR = 5) and H = 1.5 μm for ρ = 0.5 (W = 300 nm, AR = 5) (Figure S7). Lastly, to explore the torsional behavior of particles with heights exceeding those A, B, andC) are indicated by the color shading and the schematics of particle cross section along the x-y plane (defined as in Figure 1d) that illustrate different film thicknesses and material filling ratios (color-coded as in Figure 1d and structured as in Table S1). S2.

predicted to optimize torque transfer efficiency, we also fabricated a particle batch with H = 1.8 μm for ρ = 0.5 (W = 300 nm, AR = 6). The dimensions for 10 designed particle batches (A3-5, B3-5, and C3-6, where the letter and number represent the wafer and the aspect ratio, respectively) are indicated on the maps in Figure 2a,b.

The MM nanocuboid dimensions are not limited to the above-mentioned size ranges and can be enlarged or reduced as necessary. The upper boundary of particle dimensions, of interest in, for example, diverse microfluidic components, [START_REF] Higurashi | Optically Induced Rotation of a Trapped Micro-Object About an Axis Perpendicular to the Laser Beam Axis[END_REF][START_REF] Neale | All-Optical Control of Microfluidic Components Using Form Birefringence[END_REF] is set by the 3D-trappable regime for each chosen filling ratio. The lower boundary, which can be beneficial to achieving improved precision in single-biomolecule measurements, [START_REF] Van Oene | Quantifying the Precision of Single-Molecule Torque and Twist Measurements Using Allan Variance[END_REF] is determined by the specific design constraints for MM particles, namely, (i) the thinnest layer of a multilayer stack should be sufficiently thick to ensure reliable and reproducible multilayer deposition (Methods, Table S1); (ii) the number of layer-pairs should exceed 2 to permit usage of the EMT approximation in particle design (Figure S1); and (iii) the particle length should not be substantially smaller than the waist of the focused beam and should exceed the particle width (i.e., maintain a high aspect ratio) to achieve the proper trapping orientation (Figure 1d). Conforming to these criteria, the precision of multilayer deposition (∼1 nm) and the resolution of patterning (∼10 nm) in our fabrication environment (Methods) are sufficient to achieve even sub-100 nm particle width while a sufficiently long particle height is readily defined by lithography.

Linear and Angular Trapping Behavior of MM Particles. We have fabricated and characterized 10 batches of Nb 2 O 5 /SiO 2 MM particles. The fabrication process (Methods; Figures S8-S10, Table S1) was tuned to produce particles with actual dimensions as close as possible to the design parameters with high uniformity, resulting in <3% relative standard deviations (RSD = SD/mean × 100) in each dimension (Figure S10). To measure the linear and angular trapping properties of these different particles, we then introduce them into an optical torque wrench (OTW) setup [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] (Methods; Figures S11,S12, Table S2).

We first confirmed the tight 3D-trapping of the MM nanocuboids by measuring linear trapping properties along all principal axes (x, y, and z) (Figure S12). Overall, the nanocuboids show the highest stiffness along the y-axis (κ y ), and we use this parameter for further analysis and comparison (Figure 3a). We next characterized the angular operation of each 3D-trapped nanocuboid about the z-axis, which is controlled by the linear polarization direction of the input beam, and computed the resulting values of the angular stiffness (Figure 3b). Within each batch, the particles behave similarly, with average RSD values of ∼6% in the measured linear and angular parameters (RSD ≈ 1-14% for stiffness and RSD ≈ 2-13% for drag). These RSD values are as low as those obtained for highly uniform commercially available PS beads in the same setup (RSD ≈ 7-11% for stiffness and RSD ≈ 7-9% for drag), [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] which confirms that the fabricated nanocuboids are highly uniform in both optical properties and geometry.

We then compared the OTW-measured linear and angular trapping properties of the MM nanocuboids with our FEM calculations. The measured linear and angular stiffnesses show quantitative agreement with the corresponding values calculated using FEM, provided that the latter are scaled by 0.63 ± 0.31 for κ y (Figure 3a) and 0.55 ± 0.04 for κ θ (Figure 3b) (mean ± SD; for the measured particles in all batches). Such scaling factors are expected, as similar factors have been observed for other types of particles probed in the same setup (∼0.58 for PS beads and ∼0.42 for rutile TiO 2 cylinders). [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] We attribute the necessity for scaling to the fact that the FEM calculations do not take into account optical aberrations, which can substantially distort the field gradient of the tightly focused trapping beam. [START_REF] Reihani | Significant Improvement of Optical Traps by Tuning Standard Water Immersion Objectives[END_REF] Furthermore, for the particle batches from wafer B, the increased sensitivity to optical aberrations of linear trapping over angular trapping accounts for the increased variability in the scaling factor for lateral stiffness (Figure 3a) over angular stiffness (Figure 3b). This can be attributed to (i) their larger diameters (400 nm), which make it likely that the particle-medium interface lies at the boundary of the focused beam where the field gradient is more susceptible to optical aberrations and (ii) their likely more negative values for equilibrium trapping positions z eq (Figure S4), which leads them to experience a larger field gradient than those trapped at positive z eq values if optical aberrations render the shape of the focused beam axially asymmetric. [START_REF] Rohrbach | Trapping Forces, Force Constants, and Potential Depths for Dielectric Spheres in the Presence of Spherical Aberrations[END_REF] Scaling factors are also required to achieve quantitative agreement between experimentally measured linear and angular drag coefficients and the predictions made by FEM: the calculated drag coefficients are scaled by 0.84 ± 0.14 for γ y and 0.70 ± 0.08 for γ θ (Figure S12). Here, the requirement for scaling likely results from differences between the actual geometries of nanocuboids (which have rounded edges when fabricated; Figure S8) compared to the idealized geometries with sharp edges in our FEM modeling (Methods). Indeed, rotationally symmetric particle geometries measured in the same setup [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] exhibit near-unity scaling factors for the drag coefficients (e.g., spherical beads, scaling factor ∼1.02; cylindrical rods, scaling factor ∼1.04).

The highly consistent agreement between our experiments and calculations for the linear/angular drag coefficients and angular trapping stiffness in all particle batches (Figure S12) proves that the particle geometries and optical properties conform to design, reinforcing the notion that mismatches found in the experimental measurements of linear trapping stiffnesses for several particle batches result from optical aberrations. This also proves that all of the tested nanocuboids are actually trapped in an optical trap with the orientation as shown in Figure 1d, [START_REF] Singer | Orientation of Optically Trapped Nonspherical Birefringent Particles[END_REF][START_REF] Bareil | Angular and Position Stability of a Nanorod Trapped in an Optical Tweezers[END_REF] as otherwise the experimental results would have been inconsistent with the calculations that assume the idealized situation depicted there.

By examining our experimental results for MM particles with different material filling ratios, we could clearly demonstrate the ability to tune their optical constants. The particles in batches A3-5 and C3-5 have the same geometry (W = 300 nm, AR = 3-5) and hence in principle the same drag coefficients; indeed, the ratio of their measured drag coefficients is 1.00 ± 0.13 ((batches C3-5)/(batches A3-5); Figure S12). However, the particles in batches C3-5 (ρ = 0.5) exhibit measured linear and angular stiffnesses that are ∼29% and ∼40% higher than particles in batches A3-5 (ρ = 0.3), respectively (Figure 3a,b). As these results are in accordance with the predictions for the optical constants from EMT (Figure 1b), we attribute these differences to their different material filling ratios.

We illustrate how the ability to tune optical constants of MM particles can be beneficial by showing how it can be used to maximize achievable torques. For example, increasing values of the AR results in higher torques for particles with ρ = 0.5 (batches C3-5 in Figure 3b); however, the efficiency of torque transfer subsequently saturates, as shown by the modest ACS Nano increase in angular stiffness upon further increase of the AR (compare batches C5 and C6 in Figure 3b). If higher torques are desired, then increasing the particle width W can be considered. This approach is limited for particles with ρ = 0.5, as only a narrow range of widths is compatible with stable 3Dtrapping (Figure 2b). Tuning the optical constants provides an alternative: the lower index of multilayers with ρ = 0.3 expands the range of 3D-trappable particle dimensions while retaining high birefringence, allowing the stable trapping of larger particles and hence the generation of larger torques. Thus, particles with ρ = 0.3, the larger width W = 400 nm, and a height close to the theoretical optimum for 100% torque transfer efficiency (batch B5) allow for a ∼23% larger measured torque compared to similarly optimized particles with ρ = 0.5 and smaller width W = 300 nm (batches C5-6), as shown in Figure 3b.

A further parameter with which we can compare the trapping performance of MM particles is the response time t c , which is defined as the ratio of the measured drag and stiffness (i.e., the slope of the fitted lines in Figure 3c,d). We fit the data sets for MM particles of different material filling ratios and benchmark their response times against those of PS beads and rutile TiO 2 cylinders measured in the same setup. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] For facile comparison, each response time is shown for a reference value of 100 mW beam power in the optical trap, scaled from the measured values obtained at 92 mW. By comparing the linear response times for MM particles with ρ = 0.5 (t c,y ≈ 11 μs) and ρ = 0.3 (t c,y ≈ 14 μs), we confirm that these reflect the differences in their refractive indices (n ̅ = 1.82 and 1.68, respectively). As expected, both linear response times are shorter than that of PS beads with a lower refractive index (t c,y ≈ 26 μs for n = 1.45). Meanwhile, differences in the angular response times of MM particles reflect differences in their birefringence, as shown by the 2-fold faster response of MM particles with ρ = 0.5 (t c,θ ≈ 18 μs; Δn = 0.16) compared to those with ρ = 0.3 (t c,θ ≈ 35 μs; Δn = 0.13). Both angular response times are longer than that of rutile TiO 2 cylinders with a higher birefringence (t c,θ ≈ 15 μs; Δn = 0.26). Importantly, the distinguishable linear and angular response times of the MM particles with ρ = 0.3 and ρ = 0.5 again confirm that the index and the birefringence of multilayer structures can both be tuned by choosing different material filling ratios.

We complete our investigation of the angular trapping behaviors of Nb 2 O 5 /SiO 2 MM particles by plotting their full torque-speed curves (Figure 4; measured at 92 mW beam power in water) and benchmarking their torque-speed performance relative to single-crystal rutile TiO 2 cylinders (magenta stars in Figure 4; data from ref 23). Here, the optical torque is plotted as a function of the linear polarization rotation frequency (PRF) of the input beam. One initially observes a linear regime in which torque rises as PRF is increased (and the particle rotates more rapidly), followed by a nonlinear regime in which the average torque decreases as PRF is increased (and the particle rotation becomes more asynchronous with PRF). [START_REF] Santybayeva | Optical Torque Wrench Design and Calibration[END_REF] The measured optical torques excellently fit to the analytical equations in both the linear and nonlinear regimes [START_REF] Santybayeva | Optical Torque Wrench Design and Calibration[END_REF] (overlaid black lines in Figure 4). The slope of the linear part of each curve represents the particle's angular drag coefficient (γ θ = τ o /(2πf o )). As evidenced by their nearly identical slopes, particles in batches A3-5 (green squares) and C3-5 (cyan to violet circles) have similar angular drag coefficients, as expected from their identical design geometries. The higher birefringence of wafer C than wafer A (ρ = 0.5 versus ρ = 0.3, respectively), however, allows one to access higher maximal torques and angular speeds for particles from batches C3-5. Meanwhile, as can be deduced by their discernibly steeper slopes, particles from batches B3-5 (yellow to red triangles) have increased angular drag coefficients compared to particles from batches A3-5 and C3-6, as expected from their larger volumes. The highest torques were achieved with particles from wafer B (ρ = 0.3), facilitated by the larger particle sizes compatible with stable 3D-trapping at this filling ratio (but resulting in their concomitantly lower values for the maximal rotation speeds). Thus, we observe that the generation of high torques is facilitated by either high birefringence or large (trappable) particle size, but in the latter case comes at the expense of a reduced temporal response.

When compared to single-crystal rutile TiO 2 cylinders measured under the same conditions, we observe that our Nb 2 O 5 /SiO 2 MM particles can achieve up to ∼2 times higher maximum torque (compare 16.3 nN•nm for MM particles with ρ = 0.3 and 13.5 nN•nm for MM particles with ρ = 0.5, to 8.0 nN•nm for rutile TiO 2 ). While the MM particles have a [START_REF] Santybayeva | Optical Torque Wrench Design and Calibration[END_REF] are overlaid as black solid lines. The maximum torque-speed data points are marked by the corresponding symbols (see legend; color coding of symbols, color shading, and MM nanocuboid cross sections are the same as in Figure 3a,b). Until the maximum torque is reached, the angular speed ( f) of the trapped particle is identical to the PRF of the laser beam, and it is proportional to the optical torque (τ), with the angular drag (γ θ ) as the proportionality constant (τ = 2πγ θ f). When the PRF exceeds the maximum angular speed of a particle, the particle rotation becomes more asynchronous to the PRF, resulting in a lower average torque. [START_REF] Santybayeva | Optical Torque Wrench Design and Calibration[END_REF] For comparison, the maximum torque-speed values of eight differently sized single rutile TiO 2 (RT) nanocylinders are overlaid (magenta stars; measured in ref 23 under identical conditions). The small spikes near zero PRF and the ripples at high PRF are attributed to artifacts of the measurement method. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] lower value of the birefringence (Δn = 0.13-0.16 versus Δn = 0.26 for rutile), this can be compensated through the use of larger volumes. Such larger volumes, which result in a lower maximum angular speed for Nb 2 O 5 /SiO 2 MM particles (3.5 kHz for ρ = 0.3; 5.0 kHz for ρ = 0.5) compared to rutile TiO 2 cylinders (7.6 kHz), are permitted since the threshold width for stable 3D trapping of our Nb 2 O 5 /SiO 2 MM particles (∼525 nm for ρ = 0.3, Figure S5; ∼375 nm for ρ = 0.5, Figure 2b) exceeds the threshold diameter of ∼255 nm for rutile TiO 2 cylindrical particles. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF]

CONCLUSIONS

In this work, we have demonstrated a robust approach to the design and fabrication of multilayer metamaterial nanoparticles for force and torque transduction in optical trapping applications. We have validated the accuracy of the EMT approximation for predicting the optical constants of multilayer structures using analytical and numerical calculations. From among the many potential materials, we selected Nb 2 O 5 and SiO 2 as the multilayer material pair for experimental demonstration of our approach. We modeled the linear and angular trapping properties of square-cuboid MM nanoparticles as a function of particle dimension and material filling ratio using FEM calculations and fabricated selected nanocuboid designs by etching into a multilayer stack composed of alternating Nb 2 O 5 and SiO 2 layers. The nanocuboids are released from their substrate by dissolution of an underlying sacrificial layer to produce freely floating particles that remain stable in aqueous environments. By capturing such particles in an optical torque wrench and measuring their linear and angular trapping properties, we could experimentally deduce their high uniformity in size and composition from the low variability in their trapping properties. Importantly, we could demonstrate the tunability of optical constants. As such, we could realize MM nanoparticles that emulate the optical properties of calcite/ vaterite CaCO 3 crystals, while avoiding the physicochemical limitations of these crystals that have prevented them for achieving their full potential. We could also show that the trapping properties of differently designed MM particles could be tuned and understood: MM particles with higher material filling ratios exhibited higher linear and angular stiffnesses, and shorter response times, compared to particles of identical dimensions but lower material filling ratios; similarly, such particles could be rotated fastest, at speeds up to ∼5 kHz in water. The highest torques (∼16 nN•nm) were obtained for particles with lower material filling ratios, as their lower index enables stable 3D-trapping of larger dimensions. These maximal torques are 2-fold larger compared to that achieved with rutile TiO 2 nanocylinders, which have higher birefringence but require smaller dimensions to remain trappable in 3D, and moreover are more challenging to fabricate.

Our MM particles have many advantages over existing materials for optical trapping applications. First, the abovementioned tunability allows one to circumvent the limited set of optical constants accessible with naturally occurring materials and design new types of nanoparticles to, for example, access broader regimes of force and/or torque transduction. For example, large, highly birefringent, yet 3Dtrappable MM particles, such as those derived from our Nb 2 O 5 /SiO 2 multilayer with ρ = 0.3, are promising as powerful rotating microfluidic components. [START_REF] Neale | All-Optical Control of Microfluidic Components Using Form Birefringence[END_REF] Moreover, there exist specific material combinations beyond the Nb 2 O 5 /SiO 2 pair presented here to further expand the accessible range of optical constants: the a-Si/SiO 2 pair, for example, is expected to achieve exceptionally high birefringence. Second, our MMs allow for chemical stability (by pairing materials stable in aqueous environments), particle harvesting in water, and various options for surface functionalization, [START_REF] Clerico | Water-Dispersible Three-Dimensional LC-Nanoresonators[END_REF][START_REF] Sperling | Surface Modification, Functionalization and Bioconjugation of Colloidal Inorganic Nanoparticles[END_REF][START_REF] Pujari | Covalent Surface Modification of Oxide Surfaces[END_REF][START_REF] Reimhult | Design of Surface Modifications for Nanoscale Sensor Applications[END_REF] properties that are beneficial for microfluidics and biological applications such as studies of rotary motors [START_REF] Nord | Catch Bond Drives Stator Mechanosensitivity in the Bacterial Flagellar Motor[END_REF] and torque-response of biomolecules. [START_REF] Forth | Torque Measurement at the Single-Molecule Level[END_REF][START_REF] Lipfert | Torque Spectroscopy for the Study of Rotary Motion in Biological Systems[END_REF][START_REF] Van Oene | Quantifying the Precision of Single-Molecule Torque and Twist Measurements Using Allan Variance[END_REF] Third, due to the orientation of the optic axis normal to the multilayer substrate surface, the particle geometry along the long axis is lithographically defined, presenting no practical limitation on particle height or shape. [START_REF] Hernandez | Colloidal Alphabet Soup: Monodisperse Dispersions of Shape-Designed Lithoparticles[END_REF] In contrast, the fabrication of particles from birefringent substrates with an in-plane optic axis requires vertical etching, making it a challenge to realize particles with heights exceeding 1-2 μm or with unconventional sidewall shapes, particularly for hard-to-etch crystals. [START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF] The ease of realizing a longer particle is a notable advantage for torque transducers, as the optical torque transfer efficiency is largely dependent on particle height. For example, the maximal achievable torque and rotation frequency are 2-fold larger for a circularly polarized rather than a linearly polarized input beam, provided that the particle height can be made 2-fold longer than the optimal height for a linearly polarized input beam. [START_REF] Lipfert | Torque Spectroscopy for the Study of Rotary Motion in Biological Systems[END_REF] Fourth, unlike single-crystal substrate-based fabrication (e.g., quartz SiO 2 or rutile TiO 2 ), [START_REF] Huang | Electron Beam Fabrication of Birefringent Microcylinders[END_REF][START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF] our sacrificial layer-based fabrication utilizes standard silicon substrates and does not require mechanical cleavage of particles. As such, it guarantees uniform particle geometry and can be extended to mass production. This high yield and uniformity are beneficial for many torquerelated biological [START_REF] Van Oene | Quantifying the Precision of Single-Molecule Torque and Twist Measurements Using Allan Variance[END_REF] and physical [START_REF] Arita | Laser-Induced Rotation and Cooling of a Trapped Microgyroscope in Vacuum[END_REF] experiments, in which a large number of identical particles is desired to obtain statistically sound and reproducible results. Together, the designable optical properties and versatile fabrication of universal birefringent MM particles target them as highly promising candidates for future development of optimized, high-performance optomechanical transducers for use in science and engineering.

METHODS

Finite Element Method Modeling of Optical Trapping

Behavior and Hydrodynamic Drag. The finite element method (COMSOL Multiphysics 5.2a) is used to calculate the optical response of MM particles. A square cuboid-shaped particle is surrounded by a spherical region of uniform medium (water, n = 1.33), which is terminated by a perfectly matched layer shell. The scattering of the particle is calculated under background illumination by a tightly focused Gaussian beam that is linearly polarized along the x-axis and is described by the diffraction integral of Richards and Wolf. [START_REF] Wolf | Electromagnetic Diffraction in Optical Systems -I. An Integral Representation of the Image Field[END_REF][START_REF] Richards | Electromagnetic Diffraction in Optical Systems, II. Structure of the Image Field in an Aplanatic System[END_REF] Parameters used for the description of the focused beam, e.g., the vacuum wavelength of 1064 nm, the NA of 1.2, and the aperture filling ratio of 1.7, are as in our experimental setup [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] and are also used to calculate the theoretical beam power. The optical forces and torques are retrieved by integrating the Maxwell stress tensor on the surface of a virtual sphere enclosing the MM nanocuboid. [START_REF] Chen | Optical Pulling Force[END_REF] These are used to extract the linear and angular trapping stiffness values, which are then normalized by the theoretically calculated input beam power.

For modeling the inner structure of an MM nanocuboid, we use an equivalent single-layer structure with effective optical constants predicted by EMT to calculate the optical trapping maps (Figure 2) rather than the actual multilayer structure, in order to reduce the calculation time. This approach has been validated by calculating the trapping forces and torques from both structures and showing that the resulting force and torque curves have only negligible differences (Figure S2). This comparison also confirms that EMT is a sufficiently precise description of our subwavelength multilayer structures at the operating wavelength of 1064 nm.

The viscous drag coefficients are also calculated by FEM (computational fluid dynamics module of COMSOL). The surrounding medium (water) is set to flow translationally (rotationally), inducing viscous drag force (torque) on the MM nanocuboid located at the center of the calculation domain. Solving the Navier-Stokes equations yields force and torque as a function of the speed of the medium flow, from which the linear (angular) drag coefficients are extracted.

These FEM models were initially validated in our previous report 23 to model rutile TiO 2 cylinders and only altered for the current work to use a square-cuboid instead of a cylinder as the particle geometry.

Sacrificial Layer-Based Top-Down Fabrication of MM Particles. The MM particles are fabricated in a top-down process (Figure S8a) that is adapted from our previous work. [START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF][START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] First, a fourinch silicon (Si) wafer is cleaned with fuming (99.5%) HNO 3 (10 min), followed by thoroughly washing with DI water and spin-drying (Figure S8a, step 1). Next, a 100 nm thick sacrificial chromium (Cr) layer is deposited on the Si wafer by electron-beam evaporation (FC-2000, Temescal) with a deposition rate of 0.5 Å/s at a chamber pressure of ∼3 × 10 -7 Torr (Figure S8a, step 2). This slow deposition ensures a high-quality Cr coating with low surface roughness. Then, the Nb 2 O 5 /SiO 2 multilayer is deposited using plasma-assisted reactive magnetron sputtering (HELIOS, Buḧler) [START_REF] Lumeau | Angularly Tunable Bandpass Filter: Design, Fabrication, and Characterization[END_REF] (Figure S8a, step 3), yielding highly reproducible refractive indices (Figure S9), while the layer thicknesses are controlled to ∼1 nm precision with the aid of in situ optical monitoring [START_REF] Macleod | Thin Film Optical Filters[END_REF] (Figure S9). The multilayer stack is designed to be symmetric by starting and finishing with SiO 2 layers of the same thickness. The detailed design of multilayer composition is described in Table S1.

The square-cuboid MM particles are shaped by lithography and etching with a Cr hard mask. First, the substrate is machine-diced to 1 cm × 1 cm chips, which are then cleaned with HNO 3 and DI water, as in the Si wafer preparation step. The samples are further cleaned in an ultrasonic bath with acetone and isopropyl alcohol (IPA), 5 min each, then spin-dried. A layer of ∼250 nm thick positive-tone electron-beam resist (AR-P 6200.9, Allresist) is spin-coated (Figure S8a, step 4) and patterned by electron-beam lithography (EBPG 5000+, Vistec; patterning resolution of ∼10 nm) with an exposure dose of 240 μC/cm 2 (Figure S8a, step 5). The particles are patterned in a hexagonal array (Figure S8b), with a gap size of 1 μm between adjacent particles. The development of the patterned resist layer is performed in an ultrasonic bath by sequential dipping in pentyl acetate (2 min), a 1:1 mixture of methyl isobutyl ketone and IPA (1 min), and IPA (1 min), then spin-dried. A hard etch mask layer (∼70 nm thick Cr) is deposited by argon (Ar)-plasma sputtering (AC450, Alliance Concept), with a radio frequency (RF) power of 100 W, Ar supply of 20 sccm, and chamber pressure of 100 μbar (Figure S8a, step 6). These sputtering conditions are optimized for conformal deposition of Cr, allowing the enhanced etch mask shape, which results in more vertical sidewalls of particles after etching. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] For liftoff, adhesive tape (Kapton) is initially used to remove most of the top Cr layer, and then the remaining resist layer is removed by dipping the sample in PRS-3000 (J.T. Baker) resist stripper at 80 °C for 30 min (Figure S8a, step 7). After rinsing the sample thoroughly with DI water and spin-drying, a reactive ion etcher (RIE; Fluor Z401S, Leybold Heraeus) is used to etch the multilayer vertically (Figure S8a, step 8). This dry etching uses a mixture of CHF 3 (50 sccm) and O 2 (3 sccm), [START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF] at a chamber pressure and RF power of 50 μbar and 200 W, respectively, resulting in an average multilayer etch rate of ∼50 nm/ min. To ensure the complete removal of the multilayer film within the unmasked region, ∼15 s of additional etching is done after observing the endpoint by laser interferometry (LEM, HORIBA Scientific). In the same RIE machine, Ar-plasma sputter cleaning (20 sccm Ar flow, 10 μbar chamber pressure, and 100 W RF power) is performed for 5 min to remove an oxidized Cr layer from the surface of the hard mask and the sacrificial layer, which was probably induced during the previous dry etching step. This additional step enhances the result of wet etching for removal of remaining Cr layers in a later step, as oxidized Cr is less soluble in Cr etchant solution.

To retrieve the particles, the chips are manually cleaved into quarters, as the resulting smaller chips (approximately 5 mm × 5 mm each) can be inserted into a 2 mL volume plastic tube in a later step. Each chip is soaked in a Cr etchant solution (bright yellow; TechniEtch Cr01, MicroChemicals) for 10 min, where the top Cr mask and bottom sacrificial Cr layer are both dissolved (Figure S8a, step 9). Then, the chip is gently immersed in ample DI water for 30 s for the initial washing of Cr etchant droplets (the color of DI water turns light yellow after washing) (Figure S8a, step 10), followed by a second immersion in another beaker of DI water for 1 min without any agitation, which completely removes all etchant from its surface (after the second washing, the DI water remains as transparent as fresh water) (Figure S8a, step 11). Until this second immersion in DI water, most MM particles remain on the surface. However, they are visibly released from the surface when the chip is very slowly taken out of the second beaker of DI water and passed through the waterair interface while keeping its surface facing upward. The chip and the water droplet on its surface, which contains many released particles, are then jointly transferred into a 2 mL volume plastic tube containing 200 μL of fresh DI water, which is vortexed for 30 s (Figure S8a, step 12). The chip, whose surface is only clean Si after vortexing, is removed from the tube. The remaining particle solution is expected to have a concentration of ∼2 × 10 4 particles per μL, estimated from the solution volume, chip size, and pitch of the particle array (Figure S8b). As the particles are directly dispersible into water, the fraction of particles actually collected is much higher than in a previous top-down protocol for single-crystal substrates of quartz SiO 2 and rutile TiO 2 , where a sacrificial layer cannot be included and particles could be lost or destroyed during their collection with a sharp blade and a pipette tip. [START_REF] Huang | Electron Beam Fabrication of Birefringent Microcylinders[END_REF][START_REF] Ha | Tunable Top-Down Fabrication and Functional Surface Coating of Single-Crystal Titanium Dioxide Nanostructures and Nanoparticles[END_REF] By SEM inspection of the randomly dispersed MM particles on a Si substrate, the complete removal of both the top Cr mask layer and the bottom Cr sacrificial layer has been observed (Figure S8c). This observation is further confirmed by the stable 3D optical trapping and highly homogeneous trapping parameters of MM particles, as 3D-trapping would become impossible or trapping performance would be inhomogeneous due to the high scattering force if any Cr layer remained on the MM particle surface. Although the MM particle fabrication is demonstrated with rather small samples here, our method is not limited to the shown chip sizes and can be scaled up.

Measurements of Linear and Angular Optical Trapping Properties. The optical trapping experiments are conducted in our custom-built optical torque wrench setup. The setup schematic and configuration are shown in Figure S11, and more details can be found in ref 23. For the sample chamber, we use a custom-made flow cell assembled with two borosilicate glass coverslips (No. 1.5H, Marienfeld) separated by a single-layer Parafilm spacer of ∼100 μm thickness. The flow cell channel is completely filled with the MM particle solution without any air bubbles, and both the input and the output of the channel are sealed with vacuum grease (18405, Sigma-Aldrich). For all measurements, the MM particle solution is diluted until the concentration is low enough to have small probability of collision with other particles during measurements. It is notable that the laser beam power delivered at the focal plane inside the sample chamber is precisely calibrated [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] and used to normalize the measured trapping stiffness values.

We measured 10 different MM particle batches (A3-5, B3-5, and C3-6), and 3-7 particles were recorded from each batch (Figure S12, Table S2). For measurement and analysis of linear and angular trapping properties, we used our previously developed methods. [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] Notably, the frequency sweep method [START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF] is used in the torque-speed curve measurements (Figure 4).

ASSOCIATED CONTENT

* sı Supporting Information

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsnano.0c04233.

Figure 1 .

 1 Figure 1. Design principles for MM particles. (a) Optical constants of the most common materials used in optical trapping. For each material, the average refractive index n ̅ (circle) and optical birefringence Δn (triangle) are shown. (b) EMT-predicted effective indices and birefringence of the Nb 2 O 5 /SiO 2 multilayer structure as a function of the material filling ratio, ρ. Shown are calculated values for n ∥ (dashed line) and n ⊥ (dash-dotted line), as well as for n ̅ (= (n ∥ + n ⊥ )/2; dotted line) and Δn (= |n ∥ -n ⊥ |; solid line). The filling ratios chosen for demonstration and the corresponding optical constants are indicated by the vertical dotted lines and the overlaid symbols (as in panel (a)), respectively. (c) For each MM of different material combination, the EMT-predicted ranges of average refractive index (left-side bar) and birefringence (right-side bar) are shown, together with the optical constants of a natural crystal (overlaid symbols as in legend) that can be emulated by the corresponding MM. (d) Three-dimensional schematic of an MM square-cuboid particle made of Nb 2 O 5 (green) and SiO 2(blue) with its long axis aligned to the propagation direction (z-axis) of the trapping beam (red). Given an input beam linearly polarized along the x-axis, the particle's optic axis (n ⊥ ) aligns with the y-axis, while the axes associated with the higher refractive index (n ∥ ) align with the xand z-axes. The resulting optical birefringence in the x-y plane enables active rotary control around the z-axis (indicated by the curved arrow). The cross section in the x-y plane is a square of width W, and the aspect ratio AR sets the height H along the z-axis. (e) Scanning electron microscopy (SEM) image of MM particles with ρ = 0.3, W = 400 nm, and AR = 5. The three differently oriented particles are marked: bottom (1), side (2), and top (3) (as defined in FigureS10).

  that the refractive index along the ordinary axes (

Figure 2 .

 2 Figure 2. Calculated maps of the trapping stiffnesses and maximum angular speeds for MM nanocuboids. FEM-calculated maps for (a, b) the axial stiffness κ z , (c, d) the lateral stiffness κ y , (e, f) the angular stiffness κ θ , and (g, h) the maximum angular speed f o are shown as a function of aspect ratio AR and width W of the Nb 2 O 5 /SiO 2 MM nanocuboids with filling ratios of ρ = 0.3 (panels (a, c, e, g)) and ρ = 0.5 (panels (b, d, f, h)). The directions of linear (panels (a-d)) and angular (panels (e-h)) trapping are indicated by the coordinate system and the curved arrow, respectively, as shown in Figure1d. All trapping parameters presented are calculated for the case of an input beam that is linearly polarized along the x-axis and the specific configuration of the employed OTW setup (Methods) and then normalized by the theoretically calculated input laser beam power (Methods). In all maps, the solid blue lines represent contour lines for nanocuboids of constant height that range from 500 nm (leftmost) to 3000 nm (rightmost) at intervals of 500 nm, and the pixel size is ΔAR = 0.1, ΔW = 5 nm. The particle dimensions that correspond to the maximum for each map are indicated by the empty magenta circles. The black pixels in the maps of ρ = 0.5 (panels (b, d, f, h)) indicate the particle dimensions that cannot be trapped in 3D due to excessive scattering forces. The 10 different particle dimensions (A3-5, B3-5, and C3-6) chosen for demonstration are indicated by the crossing-points of the black dash-dotted lines, overlaid with corresponding symbols in each map, and their labels are displayed in panels (a) and (b).

Figure 3 .

 3 Figure 3. Measured linear and angular trapping stiffnesses and drag coefficients for MM nanocuboids in an OTW. The OTW-measured (a) linear stiffness κ y and (b) angular stiffness κ θ are shown for all 10 batches of MM nanocuboids (A3 to C6). The measured stiffness values are normalized by the precisely calibrated actual input laser beam power (Methods). The empty symbols and the error bars denote the mean and the standard deviation for each batch, while the corresponding FEM-predicted values (taken from the maps of Figure 2c,d,e,f) are shown as filled symbols. The three different groups of MM nanocuboid batches (A, B, and C) are indicated by the color shading and the schematics of particle cross section along the x-y plane (defined as in Figure1d) that illustrate different film thicknesses and material filling ratios (color-coded as in Figure1dand structured as in TableS1). (c, d) The OTW-measured stiffnesses (κ) shown in panels (a) and (b) are plotted against the corresponding OTW-measured drag coefficients (γ) (empty symbols and error bars, the same color coding as in panels (a) and (b)). The characteristic response time (t c = γ/κ) for each filling ratio ρ is deduced from the slope of the linear fit to the corresponding data points (forced through the origin; dash-dotted lines for ρ = 0.3; dashed lines for ρ = 0.5). For comparison, the fit to the linear stiffness versus drag data for PS beads is shown as a solid line in panel (c), and the fit for angular data of rutile TiO 2 cylinders is shown as a dotted line in panel (d) (measured in ref 23 under identical conditions). The data in panels (a)-(d) are summarized in TableS2.

  Figure 3. Measured linear and angular trapping stiffnesses and drag coefficients for MM nanocuboids in an OTW. The OTW-measured (a) linear stiffness κ y and (b) angular stiffness κ θ are shown for all 10 batches of MM nanocuboids (A3 to C6). The measured stiffness values are normalized by the precisely calibrated actual input laser beam power (Methods). The empty symbols and the error bars denote the mean and the standard deviation for each batch, while the corresponding FEM-predicted values (taken from the maps of Figure 2c,d,e,f) are shown as filled symbols. The three different groups of MM nanocuboid batches (A, B, and C) are indicated by the color shading and the schematics of particle cross section along the x-y plane (defined as in Figure1d) that illustrate different film thicknesses and material filling ratios (color-coded as in Figure1dand structured as in TableS1). (c, d) The OTW-measured stiffnesses (κ) shown in panels (a) and (b) are plotted against the corresponding OTW-measured drag coefficients (γ) (empty symbols and error bars, the same color coding as in panels (a) and (b)). The characteristic response time (t c = γ/κ) for each filling ratio ρ is deduced from the slope of the linear fit to the corresponding data points (forced through the origin; dash-dotted lines for ρ = 0.3; dashed lines for ρ = 0.5). For comparison, the fit to the linear stiffness versus drag data for PS beads is shown as a solid line in panel (c), and the fit for angular data of rutile TiO 2 cylinders is shown as a dotted line in panel (d) (measured in ref 23 under identical conditions). The data in panels (a)-(d) are summarized in TableS2.

Figure 4 .

 4 Figure 4. Measured optical torque transduced by MM nanocuboids in an OTW. The optical torques transduced by single optically trapped MM nanocuboids in water (one particle per batch) versus input beam polarization rotation frequency (PRF) as measured in an OTW. The laser power at the optical trap is ∼92 mW. The moving-average torque signals are displayed as colored solid lines, and the corresponding fits to the theoretically predicted response[START_REF] Santybayeva | Optical Torque Wrench Design and Calibration[END_REF] are overlaid as black solid lines. The maximum torque-speed data points are marked by the corresponding symbols (see legend; color coding of symbols, color shading, and MM nanocuboid cross sections are the same as in Figure3a,b). Until the maximum torque is reached, the angular speed ( f) of the trapped particle is identical to the PRF of the laser beam, and it is proportional to the optical torque (τ), with the angular drag (γ θ ) as the proportionality constant (τ = 2πγ θ f). When the PRF exceeds the maximum angular speed of a particle, the particle rotation becomes more asynchronous to the PRF, resulting in a lower average torque.[START_REF] Santybayeva | Optical Torque Wrench Design and Calibration[END_REF] For comparison, the maximum torque-speed values of eight differently sized single rutile TiO 2 (RT) nanocylinders are overlaid (magenta stars; measured in ref 23 under identical conditions). The small spikes near zero PRF and the ripples at high PRF are attributed to artifacts of the measurement method.[START_REF] Ha | Single-Crystal Rutile TiO 2 Nanocylinders Are Highly Effective Transducers of Optical Force and Torque[END_REF]

Downloaded via 194.167.230.102 on November 30, 2020 at 12:37:18 (UTC).See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

https://dx.doi.org/10.1021/acsnano.0c04233ACS Nano 2020, 14, 14895-14906

ACKNOWLEDGMENTS

We thank Joseph Braat for advice on optical modeling; Daniel Tam for advice on developing numerical models for drag coefficient calculation; Marc Zuiddam, Charles de Boer, Marco van der Krogt, Hozanna Miro, Arnold van Run, Anja van Langen-Suurling, Wiel Evers, Majid Ahmadi, and Iman Esmaeil Zadeh for technical advice on nanofabrication; Margreet Docter, Roland Dries, and Jacob Kerssemakers for technical advice on the OTW setup; Dimitri de Roos for the production of custom-made mechanical parts for the OTW setup; Theo van Laar and Sacha Khaiboulov for help with laboratory infrastructure; and Kaley McCluskey, Belen Solano, and Jacob Hoogenboom for reviewing the manuscript. This work is supported by the Faculty of Applied Sciences at Delft University of Technology (scienceDATE grant 2017; to Y.T. and S.H.); by funding provided by the China Scholarship Council under grant number 201306160089 (to Y.T.); by NanoNextNL, a micro-and nanotechnology consortium of the Government of The Netherlands and 130 partners (to N.H.D.); by the Department of Bionanoscience at Delft University of Technology (to N.H.D.); by the Foundation for Fundamental Research on Matter (FOM), part of The Netherlands Organisation for Scientific Research (NWO) (to N.H.D); and by the European Research Council (ERC) Consolidator Grant under grant number 312221 (DynGenome; to N.H.D.).

Notes

The authors declare no competing financial interest.

NOTE ADDED AFTER ASAP PUBLICATION

This article was published ASAP on November 10, 2020, with incorrect Supporting Information. The corrected version was reposted on November 11, 2020.

ACS Nano