Mia
Kvåle Løvmo
,
Benedikt
Pressl
,
Gregor
Thalhammer
and
Monika
Ritsch-Marte
*
Institut für Biomedizinische Physik, Medizinische Universität Innsbruck, Müllerstraße 44, 6020 Innsbruck, Austria. E-mail: monika.ritsch-marte@i-med.ac.at; Fax: +43 512 9003 73871; Tel: +43 512 9003 70871
First published on 26th February 2021
In cell biology, recently developed technologies for studying suspended cell clusters, such as organoids or cancer spheroids, hold great promise relative to traditional 2D cell cultures. There is, however, growing awareness that sample confinement, such as fixation on a surface or embedding in a gel, has substantial impact on cell clusters. This creates a need for contact-less tools for 3D manipulation and inspection. This work addresses this demand by presenting a reconfigurable, hybrid sono-optical system for contact-free 3D manipulation and imaging, which is suitable for biological samples up to a few hundreds of micrometers in liquid suspension. In our sono-optical device, three independently addressable MHz transducers, an optically transparent top-transducer for levitation and two side-transducers, provide ultrasound excitation from three orthogonal directions. Steerable holographic optical tweezers give us an additional means of manipulation of the acoustically trapped specimen with high spatial resolution. We demonstrate how to control the reorientation or the spinning of complex samples, for instance for 3D visual inspection or for volumetric reconstruction. Whether continuous rotation or transient reorientation takes place depends on the strength of the acoustic radiation torque, arising from pressure gradients, compared to the acoustic viscous torque, arising from the shear forces at the viscous boundary layer around the particle. Based on numerical simulations and experimental insights, we develop a strategy to achieve a desired alignment or continuous rotation around a chosen axis, by tuning the relative strengths of the transducers and thus adjusting the relative contributions of viscous and radiation torques. The approach is widely applicable, as we discuss in several generic examples, with limitations dictated by size and shape asymmetry of the samples.
Ashkin's optical tweezers,4,5 particularly in the form of holographically shaped and steered tweezers,6–9 provide a versatile tool to accomplish this goal. Dielectric particles in water solution can be held, moved, or rotated in a contact-free way in holographically reconfigurable trapping landscapes. Yet, optical trapping has its limitations: for large and heavy multi-cellular specimens the application of single-beam optical tweezers becomes problematic, because the required optical powers to create sufficiently strong forces may give rise to photodamage. This limit can be extended towards larger objects with optical ‘macro-tweezers’,10 a specially designed optical dual-beam mirror-trap with large field-of-view and trapping volume. With this device it was possible to hold, e.g., swimming micro-organisms of ∼70 μm length and to force them on a trajectory with forces exceeding 100 pN with minor temperature rise (contrary to what is sometimes stated in the literature11). But for significantly larger organisms this approach also reaches its limits.
This is where acoustic trapping in liquids becomes an interesting alternative.12 Coupling acoustic waves into the sample solution causes pressure changes in the fluids, which result in the well-known acoustic radiation forces.13 At frequencies of a few MHz, the wavelength in water is a few 100 μm, which gives strong forces on samples of the same size scale, with very low impact on the viability of living specimens.14
Compared to optics, the spatial resolution of acoustic trapping is reduced due to the larger wavelength. Nevertheless, in acoustics an impressively high level of control of positioning can be reached, e.g., by employing surface acoustic standing waves,15,16 also with dynamic shaping.17,18 Also, phased-array emitters around a liquid-filled chamber19 or diffractive beam shaping,20 are employed, for instance for cell-sorting,21 or patterning in ultrasonic gels.22 Taking holographic optical tweezers (HOT) as a model, holographic acoustic tweezers (HAT) enable control of individual particles levitated in water23 or air.24 For the relatively large samples considered here, it is advantageous to use a simpler approach based on bulk acoustic waves rather than surface acoustic waves, since larger specimen size requires larger chamber volumes.25 Moreover, we are aiming at full 3D flexibility in the manipulation, far away from any surface and with independent control of all directions.
An example of a large living specimen which can be trapped acoustically, but not optically, is the tardigrade (Hypsibius dujardini), known colloquially as water bear or moss piglet, shown in Fig. 1 and ESI† Movie S1. This specimen, which measures about 190 μm by 50 μm in size, clearly would be beyond trapping with optical tweezers.
Fig. 1 Acoustic trapping of a large and actively moving live specimen: a tardigrade or ‘moss piglet’ is levitated and confined in 3D (ESI† Movie S1). |
As the two physical trapping mechanisms complement each others strengths and weaknesses,12 it is obviously highly advantageous to combine them in a single device: holographic optical tweezers operate on a finer scale, while the ultrasound can do the heavy lifting, quite literally. We follow the suggestion of Memoli and Fury,26 and name this combination of acoustic and optical trapping a ‘sono-optical’ device, in order to avoid confusion with acousto-optics or photo-acoustics. We expect a benefit from such devices in the precise handling and inspection of large (>100 μm) living specimens, such as microorganisms, embryos, cancer spheroids, or organoids. One would like to keep complex samples in a natural environment, unhampered by confining structures (e.g., when embedded in gel) or walls (e.g., when stored in micro-well plates or squeezed between glass slides), especially for long term monitoring.
In our previous work we demonstrated such sono-optical devices for the first time, combining acoustic levitation with optical macro-tweezers,68 and later with single beam holographic optical tweezers. In both cases acoustics were mostly used for levitating particles against gravity within a plane, and optical manipulation for precise control within this plane. This approach also opens the path to high-quality optical imaging of levitated trapped specimens and detailed measurements of acoustic forces.27–29 Hybrid sono-optical chips may profit from complementary features of the two modalities,30 which has, e.g., been utilized for improving the separation of sub-populations of white blood cells based on differences in both, mechanical and optical properties.31
When aiming at the manipulation of specimens of increasing size, of 100 μm or more in diameter, one faces new challenges, which we have tackled in this work:
• For increasing specimen size, holographic optical tweezers lose their dominant role over acoustics, and it becomes desirable to have more flexibility in the acoustic forces. Therefore, going beyond the simple 1D standing wave used previously, we have implemented 3D acoustic trapping with three independent transducers for each of the three orthogonal directions.
• As the size (and shape) of the device needs to be scaled to accommodate larger samples, the acoustic modes inside the up-scaled 3D acoustic trapping device are affected and more complex patterns appear, which required investigation in theory and experiment.
• For large samples, retrieving knowledge of the 3D structure is problematic, as standard microscopic imaging is hampered by the specimen itself. When focusing deep into the sample the image quality is degraded due to scattering, absorption, and particle induced aberrations. A simple approach to obtain 3D information is to rotate the particles (and not the imaging system) to obtain images from different viewing directions for optical tomographic imaging, using the extended micro-manipulation possibilities provided by our sono-optical trapping device.
• Acoustic trapping acts on a relatively large spatial scale compared to optical manipulation. Therefore, in our device we have implemented auxiliary holographic optical tweezers. The fact that the ‘heavy weight lifting’ is accomplished by the acoustic fields facilitates their usage, since the optical traps can be switched on and off at ease without ever losing stable trapping conditions for the sample. With this, (re-)positioning can be done on a scale much smaller than the typical acoustic wavelengths or particle sizes we use.
As a major contribution, we report observation of controlled reorientation and controlled continuous rotation of specimens in our devices, either when using only acoustics, or triggered by optical manipulation. We also present practical considerations to enable productive use in research.
The manuscript is organized as follows: sec. 2 explains the design, assembly, and properties of the microfluidic chambers. It is followed by sec. 3, a brief summary of the relevant phenomena giving rise to acoustic torque on a particle. Sec. 4 explains the general strategies to induce rotations and exemplifies our experimental results on the controlled reorientation and the sustained rotation of various specimens.
Fig. 3 Images of assembled sono-optical devices. Top: Assembly based on ceramic plate with side-transducers mounted on metal wedges. Bottom: Monolithic design made of metal. |
Three individually driven orthogonal transducers allow us to independently control the strength of the acoustic confinement in all directions. We preferably select bulk acoustic waves with clear directionality, to obtain good independent control. In the vertical direction we require the strongest confinement to hold particles against gravity. To achieve this, we place an ultrasound transducer made out of lithium niobate (LiNbO3) in direct contact with the fluid on top of the chamber, which provides good coupling for the ultrasound wave. A microscope cover slip on the bottom side seals the chamber and acts as a reflector for the acoustic wave. To excite the ultrasound wave in the horizontal directions we use side-transducers mounted on metal wedges glued to the base device substrate (see upper image in Fig. 3), or, in a monolithic design, attached to tapered metal structures, see lower image in Fig. 3.
The optically transparent LiNbO3 top-transducer provides very good optical access to the chamber,32 supporting optical trapping and imaging. Moreover, the chip is designed to be fully compatible with standard microscopy setups, i.e., it easily fits into an inverted microscope, which we use for standard wide-field imaging. More details about imaging and the details of the optical trapping are described in ESI† S1 and the optical setup is shown in ESI† Fig. S1.
Surface-milled channels at the bottom side (width 1 mm, depth 300 μm) connect to opposing diagonal corners of the chamber and act as inlet and outlet for filling the chamber. The channels end in 1 mm diameter through-holes, which provide a water-tight connection to standard pipette tips. On the bottom side, chamber and channels are sealed with a square 0.17 mm thick coverslip with a side length of 20 mm. The layout is simple enough so that hand-milling with a single 1 mm diameter tool is feasible and efficient.
Since the fabricated chamber has a horizontal cross-section with rounded corners and thus deviates from a perfect cuboid, we have carried out simulations based on finite element methods (FEM) in order to get a general idea of the modes in our chamber. For this we look into the solutions of the wave equation
(1) |
By tuning the individual voltages of the transducers one can control the pattern. Driving a single transducer excites a roughly planar wave, while exciting all three transducers confines particles to well separated trapping spots arranged in a more or less regular grid. We prefer this configuration for trapping of large, living specimens. Each trapping site provides a 3D confinement, where the trap stiffness along each direction can be controlled individually. As we will demonstrate in sec. 4.1, this allows us to control the orientation of trapped particles.
But even when we drive a single transducer, we observe some 3D confinement, e.g., additional weaker forces acting in the horizontal direction, when only the top-transducer is actuated. Furthermore, we observe some deviation of the vertical trapping position with horizontal position from an ideal flat arrangement. As our 3D simulations indicate, we attribute this to a simultaneous excitation of more than a single mode by the top-transducer, including modes with nodes along the horizontal directions, which is facilitated by small asymmetries in our actual resonator. Our calculations predict a quite dense set of modes, with unequal frequency spacing of the resonances with a mean value of 6 kHz and a maximum of 20 kHz in the range of 2 MHz to 2.4 MHz for the 3D case. Due to this dense spacing it is likely to excite several resonances, which are observed to have a width of several kHz. We suppose that this simultaneous excitation of at least two modes by a single frequency is a key ingredient to induce sustained rotation by a mechanism which we discuss in sec. 3.
In cases when we prefer flatter, more uniform horizontal confinement layers from the top-transducer, we sometimes additionally apply frequency dithering to the top-transducer. For this we rapidly modulate the excitation frequency, with frequency deviation of 50 kHz at a rate of 1 kHz. This simultaneously excites several neighbouring acoustic modes within a bandwidth of about 100 kHz with strongly varied phase relations. Individual differences in the trapping landscape are averaged out, which leads to more regular and uniform trapping. Furthermore, the sensitivity to small changes of the resonance frequencies, e.g., due to temperature changes, is reduced.
To overcome this issue, we propose a simple and convenient method by employing actively swimming micro-organisms as test particles.10 They quickly disperse if not confined by acoustic forces. We prefer to use Euglena gracilis, an easy to cultivate elongated protist with a typical length in of about 50 μm (in its elongated form). Each acoustic resonance identified from the electrical impedance measurement can be easily but qualitatively compared with respect to levitation strength and resulting mode pattern. Switching on the top and one of the side-transducers, one can confine the specimens within the horizontal nodal planes into stripes. The frequency where the Euglena can be best restricted in their motion is a good indication where one can find the strongest resonance in the liquid. As the operator receives direct and immediate feedback, one can use this method to first find the strongest resonances and then to visualize the mode pattern at the desired driving frequencies.
We observe experimental patterns that are similar to the simulated ones, as illustrated in Fig. 5 which exemplifies the trapping landscape at a resonance near the 3rd harmonic of the x-transducer at 6.678 MHz (with the top-transducer also active). The results of the simulation and experimental data obtained from the trajectories of the swimming micro-organisms are shown back to back. Due to the imperfections of our machine-milled device, no exact quantitative comparison between experimental and simulation data can be made. Nevertheless, the pattern of this resonance reproduces some of the qualitative behaviour of our experimental findings quite well, such as the wavy horizontal trapping lines and the fact that the trapping lines become less defined towards the top and bottom margins of the chamber, where one is obviously also less successful in confining the swimming protists.
The visualization method we have used here is inspired by a popular approach in medical imaging to visualize structures of low intensity in a specific volume, called minimum intensity projection (MIP).35 Here, a 2D image is generated from a 3D volume where each pixel displays the lowest attenuation value in the column of voxels along a chosen direction. In our case, the ‘third dimension’ is the time in the recorded movie. Analyzing every pixel in the top-view of background-subtracted movie recordings of Euglena swimming in the acoustic landscape and keeping just the darkest (most attenuated) value in the time-line of the given pixel reveals local confinement: this type of ‘projection images’ shows where the specimens spend some time during the recording. Areas of relatively stronger and weaker confinement/forces can be mapped out, although not in a fully quantitative way. The procedure needs to be carried out for each identified resonance, but only once per assembly.
The micro-organism-based approach can be used to map out the strongest resonances in the chamber under 2D and 3D operation. In the result shown in Fig. 5 the transducers along x and z were active, and in Fig. 6a) the other side-transducer along y together with the z-transducer. With all three transducers active, the specimens are confined to a combination of islands and stripes, as shown in Fig. 6b). Due to the limited size of our sensor we cannot image the entire xy-cross-section of the chamber at once and the camera has to be rotated 90° to access the full propagation direction of each side-transducer.
Fig. 6 Visualization of patterns in the resonator by means of micro-organisms: sub-images a) and b) are in the same chamber as shown in Fig. 5, filled with actively swimming Euglena gracilis, in a) the top-transducer (waves travelling along z) and the side-transducer along y are switched on, in b) all transducers are on (x, y, and z). The frequencies of the side-transducers are close to the 3rd harmonic frequency of the piezo. For the visualization every 10th frame of a movie recorded at 50 Hz was taken. The scale bar is shared between images. |
Depending on the location inside the chamber, for a 10 μm polystyrene (PS) tracer bead and 10 V applied voltage, we found maximum forces in the range of several pN to tens of pN with the side-transducers, with the top-transducer we can achieve forces between tens of pN up to 100 pN. Please note that voltage refers to peak-to-peak voltage throughout the paper. Maximum pressure from side-transducers at 10 V ranges from 30 kPa to 200 kPa. For the top-transducer we have measured a range of 100 kPa to 300 kPa at 10 V for the 10 μm bead in different places. The trap stiffness at 10 V is 0.8 pN μm−1 to 5 pN μm−1 for top-transducer, 0.05 pN μm−1 to 0.6 pN μm−1 for side-transducers. Since the acoustic radiation force scales with particle volume, forces of several thousand pN are achievable for larger PS beads, e.g., of 30 μm in diameter.
For long term monitoring of biological samples it is essential to estimate what kind of flow rates are compatible with these forces, for instance to ensure trapping of the particle while changing out the buffer-solution. As calculated from the maximal velocity displacement of a 30 μm bead, even for the weakest excitation direction the forces are at least on the order of 100 pN, and thus it would require only about 10 s to replace the solution, at a fluid flow speed of 400 μm s−1. Typically, the forces we operate with are even higher and a buffer-exchange can be done without affecting trapping.
Applying torque to an object by incident waves generally requires transfer of angular momentum. Optics and acoustics share some common aspects concerning linear and angular momentum flux.36,37 The main differences between light and sound waves in this context are the consequences of acoustic pressure waves in fluids being longitudinal rather than transverse, as well as the direct mechanical influence of the liquid medium, e.g., via acoustic streaming.38–40
This work mainly focuses on acoustic waves to induce torques, but let us, nevertheless, start by stating some facts from optical micro-manipulation as a starting point to understand the commonalities and differences to acoustic wave torque. In optics, rotations can be induced by the absorption or scattering of light, transferring angular momentum encoded on the polarization of the beam (spin) or the phase gradient (orbital angular momentum).41,42 Absorption of photons carrying angular momentum induces a rotation of the object. In some situations, however, absorption may be an impractical way to rotate an object, since it also leads to heating. A non-zero torque can also result if an asymmetry in the scattering is present, e.g., for particles of asymmetric shape or anisotropic material.43,44 And chiral particles may even experience a torque in optical fields which have no angular momentum.45
Non-zero time-averaged torques on an object mediated by acoustic fields have been studied for many years since the 1950s,46–49 theoretically as well as in view of lab-on-a-chip applications. The components of an acoustically generated torque follow the general scaling rule,36,50
T = τ0Q with τ0 = VE0, | (2) |
As an example, a small spheroidal particle orients itself such that the shortest axis will align with the direction of strongest confinement (trap stiffness). An oblate spheroid in a plane standing wave will behave as if squeezed in by the acoustic force field into the nodal plane. For small particles of spheroidal shape (with long axis a and small axis b) one can derive an explicit expression for the acoustic radiation torque.50 Such a particle will experience a maximal torque when oriented at 45° to a plane wave, but no torque when placed at an angle of 0° as well as 90°. When the object shape approaches a sphere the torque is found to vanish. Due to the typically strong viscous damping of the motion of small particles in a fluid, a particle will usually approach a stable equilibrium position and orientation without oscillations. Thus, the acoustic radiation torque provides us with a tool for alignment or reorientation of non-spherical particles.
To describe the viscous torque in such a situation in the framework of eqn (2),
(3) |
In Fig. 7a) the superposition of two orthogonal stripe patterns with eigenfrequencies f1 = 2.237 MHz and f2 = 2.404 MHz is chosen as an example, which comes close to the idealized situation of orthogonal plane waves. For the modal width we have assumed b = 80 kHz (HWHM) which lies in the range of the actual values in our chambers. The small images at the top display the intensity maps of the two eigenmodes and their superposition, and in the phase map below only low-pressure regions are plotted (values above a chosen threshold at 40% of the amplitude are left white) and trapping spots (zero pressure sites) are marked by black dots. The phase map of the pressure field demonstrates the existence of several pressure nodes with phase vortices.
However, realizing such a pure situation in practice requires particular care, e.g., a square chamber with two orthogonal transducers operating at the same frequency.56 Nevertheless, in our simulations we have seen that the occurrence of dark spots with local angular momentum is generic in our type of resonator. As described in sec. 2.3.2, the rounded edges lead to a coupling between the horizontal directions, i.e., any eigenmode contains contributions from plane acoustic waves travelling in orthogonal directions. Similar phase structures, as shown in in Fig. 7b), also occur in the superposition of the first eigenmode (with pronounced stripes) with its adjacent eigenmode, but the vortex structures are less symmetric. Inevitably this nearby mode is also excited to some extent when addressing the first mode due to the finite resonance width and the small spacing of 28 kHz.
Repeatedly changing the orientation of the particles by acoustic radiation torque requires some modulation of the incident field(s). In the literature, different methods have been investigated to rotate elongated particles, short glass fibers, or polymer particle clusters, either in a step-wise way by changing the propagation direction of the ultrasound as a function of time or in a continuous way by sweeping its amplitude, frequency or phase.49,59,60
A particle of roughly elliptical shape will settle into an orientation where its shortest axis aligns with the direction of steepest trap stiffness at the trapping spot. This minimizes the potential energy of the object in the acoustic trapping potential. In our device, we can adjust the orientation of the particle in a controlled way by changing the driving voltages of the three transducers. When the top-transducer is strongest, we find that oblate particles (of relatively uniform acoustic contrast-factor distribution) align flatly within the horizontal trapping nodes, with their shortest axis being vertical, which amounts to maximal cross-section of the particle with the nodal xy-plane.
At low voltages of the side-transducers, the particle's location is confined, but the orientation of the particle within the nodal plane is not strongly locked in. Raising the voltage of one side-transducer reorients the particle's largest axis in a controlled manner along the axis with the weakest confinement, while the vertical position of its shortest axis is still maintained by the (still dominant) top-transducer. When one of the side-transducers becomes the strongest, the particle's shortest axis is now aligned with the corresponding trapping direction, the x direction for the x-transducer, and the y-direction for y-transducer. This is demonstrated in Fig. 8 with a 45 μm long oblate starch particle. The particle is trapped with three active transducers around 6 MHz, and by increasing the voltage of one of the transducers at a time, the orientation of the particle can be controlled. The top-transducer at 6.844 MHz with frequency modulation at 50 kHz is kept at 3 V, while x- and y-transducers at 6.678 MHz and 6.807 MHz, respectively, are set within a range of 1.4 V to 5.2 V and 6.8 V to 17.6 V (low to high voltage).
In practice both, acoustic radiation torque and viscous torque, are present simultaneously. The actual behaviour observed in an experiment is governed by an interplay between these mechanisms and other contributions, leading to a complex and diverse behaviour. We present examples that demonstrate strategies for inducing rotation in a controlled manner by a combined action.
Fig. 9 shows a selection, with examples of increasing complexity. Fig. 9a) shows a rotating single, relatively small pollen grain from Tragopogon pratensis, 9b) a larger one from Sansevieria trifasciata. c) Shows a cluster of starch grains of various sizes, some actively rotating, some not. Finally, in d) we see a cancer spheroid performing a sustained rotation.
Fig. 9 A wide variety of individual and clustered samples (ranging from 35 μm to 170 μm in length) performing a sustained rotation in the acoustic traps around an axis orthogonal to the imaging axis. One pollen grain from a) Tragopogon pratensis and b) Sansevieria trifasciata. In panel c) a potato starch cluster of which several connected grains are rotating. d) Shows a cancer cell spheroid. The scale bar is shared between images (ESI† Movies S2–S5). |
From this we gather that sustained rotational behaviour is commonly observed. However, we find that the behaviour depends on several parameters, not all of which are sufficiently well known in a given experimental setting. The shape of the particle, as well as their inner structure, affect the acoustic forces and torques. Particles with roughly spherical shape are more likely to show sustained rotation. Also the local field at the trapping site plays an important role. We often observe different behaviour of similar objects in neighboring trapping sites. This is due to the intricate force landscape in our resonator, as presented in sec. 2.3 and Fig. 7. Furthermore, the local field can be significantly altered by particles nearby due to secondary acoustic radiation forces (Bjerknes forces). In particular, air bubbles, which strongly scatter the acoustic field, show a strong effect. Let us remark here that planar arrays of localized micro-bubbles which oscillate under acoustic excitation are an alternative way to induce controlled rotation61 of several single cells, but samples with gas-filled cavities and bubbles behave quite differently compared to solid micro-objects62 and will not be treated here.
We always see the combined effects of acoustic radiation torque, acoustic viscous torque, and possibly additional torque due to acoustic bulk streaming around the particle, all acting together. As described in sec. 3.1 the acoustic radiation torque acting on its own ‘locks’ an asymmetric particle into a certain orientation. An approximately spheroidal particle with axes lengths a < b < c is aligned such that the shortest axis a points into the direction with strongest confinement, whereas the longest axis aligns along the weakest confinement. To induce sustained rotation, the largely orientation-independent viscous torque needs to overcome the (maximum) acoustic radiation torque. This can be achieved by tuning the strength of the individual transducers (operating at different frequencies) such that the confinement strength (trap stiffness) is the same in each direction. Close to the trapping site this yields an isotropic trapping potential, where the total acoustic radiation torque vanishes also for an asymmetric particle.
For large particles it is more challenging to reach this condition. One reason is the different scaling properties of the torque contributions, according to eqn (2) and (3). The peak contribution of the acoustic radiation torque (i.e., the peak value for a specific orientation) scales as V ∼ a3, whereas the viscous contribution of the torque as ∼Vδ/a ∼ a2δ. With δ < 1 μm, δ/a is small at driving frequencies of a few MHz, which means that the acoustic radiation torque grows faster with particle size than the viscous torque (given the same acoustophoretic contrast and shape asymmetry). Therefore, for large asymmetric particles the acoustic radiation torque is usually the dominant contribution, but leads only to transient motion, tilting the particle into a certain angular orientation rather than sustained rotation. Furthermore, large particles experience more of the trapping landscape away from the trapping site, where it may become anisotropic.
Movies of spinning particles, such as the spheroidal starch grain of Fig. 10, may be used to investigate the relative portion of ‘accelerating’ and ‘decelerating’ parts in the rotational motion. In the ESI† Movie S6 one can see that the grain spends quite a long time in one orientation and the quickly rotates around to the opposite side, which is a clear indication of the existence of an angle-dependent contribution to the torque.
Fig. 10 Sinogram of an ellipsoidal starch grain rotating with a strongly inhomogeneous rotation velocity. The acoustic radiation torque, which is angle-dependent as it depends on the object's orientation, is accelerating and decelerating the rotation around the long axis induced by the viscous torque. Tracking the small feature attached to the surface (marked by a white arrow) and fitting its trajectory to a curve resulting from a constant and an angle-dependent rotation rate allows one to estimate the ratio of the two torques. In the ESI† Movie S6 one can see that the grain spends quite a long time in one orientation, which is followed by a quick revolution around the opposite side, since the two torques in this example are of similar strength (with the acoustic radiation torque amounting to ∼75% of the viscous torque). |
The sinogram in Fig. 10 was created by applying MIP inside the red rectangle. Tracking the small feature attached to the surface (marked by a white arrow) and fitting the resulting trajectory to a rotation containing a constant and an angle-dependent angular velocity allows one to derive an estimate of the relative strength of acoustic viscous and radiation torque. Assuming that the orientation angle φ(t) for rotation of the particle about its long axis follows the differential equation φ(t) = Ω1 + Ω2sin[2φ(t) + φ0], the projected position of the marked feature as a function of time t is proportional to sinφ(t), found by integrating the differential equation. Here we use the model that the viscous torque is described by an angle-independent contribution Ω1, and the acoustic radiation torque depends on the orientation-angle with a peak value of Ω2 of periodicity π. The best fit of this model, which is shown as a dashed line in Fig. 10, gave us the result that in the present case the two torque contributions are of similar strength, with the acoustic radiation torque amounting to ∼75% of the viscous torque. Rough estimates based on eqn (2) and (3) are consistent with this.
Fig. 11 shows the onset/stopping of rotation and changes of the axis of rotations and the rotation rate 1/T as a function of increasing and then again decreasing driving strength of the x-transducer. Initially, the top-transducer creates the strongest field, the x-transducer is on low voltage, the y-transducer off, and the pollen aligns ‘flatly’ with the horizontal nodal plane of the strong z-transducer, with its shortest axis locked-in vertically by the top-transducer, as expected from above considerations (we remark that in the current example the top-transducer alone also induces some confinement along the y-direction, as a result the initial orientation of the particle is as if under the influence of the y-transducer, although it is not switched on).
Fig. 11 Changing the rotation axis and rotation rate of a Hibiscus pollen grain (trapped at the site indicated in the image on the right) by tuning the strength of one of the transverse transducers. The graph shows that increasing the voltage (diamond-shaped data points in gray) beyond about 7 V there is a threshold for the onset of rotation. Beyond this threshold the alignment of the rotating particle is also different, as shown in the middle row with the blue line indicating the major axis of the grain which is 126 μm long. For a better visualization of the two different orientations of the particle, in the diagram we use green and blue background color-coding. On the way ‘down’, reducing the voltage again (circular data points in black), rotation is seen to persist down to a slightly lower voltage (ESI† Movie S7 shows the particle rotating at 1 Hz and at 2.4 Hz). |
When the voltage of the x-transducer is increased, first no motion occurs, until a certain threshold (just above 7 V) is reached, where the particle spontaneously readjusts, aligning itself according to the potential created by a now stronger x-transducer. This change of orientation is visible as a change in apparent cross-section of the pollen. In Fig. 11 we highlight this by color-coding, with a blue background indicating the alignment under the influence of a dominating top-transducer, and a green one indicating alignment under the influence of the x-transducer (particle alignment corresponding to the graph is shown in the middle row in Fig. 11). The forces generated by the top-transducer in this experiment lie in a range of 10 pN to 60 pN, and comparable forces can only be reached by the x-transducer at the highest investigated powers (reaching 8 pN to 26 pN).
Beyond this threshold, the particle starts to rotate around an axis parallel to the y-axis, as shown in the graph, rotating between alignment of the two largest cross-sections to either top- or x-transducer nodes. Increasing the voltage further (along the gray diamond markers in the curve) in steps, results in faster rotation. Eqn (2) predicts a quadratic dependence of the torque on the pressure p—which we have verified for single transducer operation—but the dependence is more complicated when two transducers are used. Similar to the experiment of Fig. 10, the rotation velocity is not constant, and the particle spends more time with its largest cross-section aligned to the top-transducer nodes. The time in this alignment decreases for larger voltages of the x-transducer.
Decreasing the driving voltage again (along the black line) results in a slightly shifted curve for the rotation rate. On the way down, the pollen stays somewhat longer in the x-transducer alignment (until ∼6.5 V is reached). This delay in settling into the new stable orientation is common in systems with friction: one cannot predict exactly when the flipping into the new equilibrium position is going to take place, it depends on the exact nature of the stochastic micro-environment around the particle. The rotation stops after the reorientation has happened. We often observe small differences in the characteristics of the rotations (onset/stopping, direction of rotation axis) when repeatedly increasing and decreasing the voltages, as some small changes may take place (e.g., a slight translation of the trapped particle or of an air-bubble or particle somewhere), which make it difficult to maintain the exact same operating point (particle location and orientation).
Carrying out the same experiment as shown in Fig. 11 in the other direction (tuning the y-transducer) did not lead to the same kind of readily interpretable results. The influence of the y-transducer seemed to be masked by other effects that were superposed. The overall behaviour was not as stable, and often ended in a rotation around the z-axis for large driving voltage. This asymmetry in behaviour is probably caused by the fact that the frequency spacing between the two transducers (z and y) was smaller than in the first case (z and x), changing from a difference of ∼170 kHz to one of ∼40 kHz, giving rise to crosstalk effects which affected the trapping pattern in 3D. Additionally, since we do not have a xy-symmetric chamber, the situation may be less favourable for inducing a controlled rotation when exciting from the other direction around the given frequency. The bulk streaming excited in the chamber also seemed to affect the particle's motion. The MIP visualization technique introduced in sec. 2.4 provides a useful means to directly demonstrate the presence of streaming in the (wider) periphery of the trapped pollen. In Fig. 12, traces from debris in the liquid show that streaming is taking place. Locally around the particle it largely takes its direction from the direction of propagation of the field emanating from the side-transducers, i.e., x in a) and y in b), but in the region around the examined object the streaming was observed to be more complicated in situation b) than in a).
Let us finally remark our observation that for even larger samples than the Hibiscus pollen, when the diameter exceeds the acoustic wavelength (e.g., 250 μm for 6 MHz) a problem occurs: such samples are practically fixated, because they are sufficiently large to be firmly ‘locked-in’ into two trapping spots along the z-direction, which makes it difficult or even impossible to induce continuous rotation in the described way.
The optical trap in the form of one or more laser spots can be dynamically steered around the acoustically levitated particle in a highly controlled way. If we target a (vertical) laser beam at the edge of an acoustically levitated particle, we expect the vertical component of the optical force to rotate the particle. The horizontal component will lead to a translation (typically towards the beam). To induce a continuous rotation, the radiation pressure must be kept off-centered on the particle. Thus balancing the horizontal optical force component to fix the particle position is advantageous, which can be done both with additional optical traps or with horizontal acoustic confinement. In the case of additional optical traps, for example two traps would be used to hold the particle on a symmetry axis, and a third would be placed off-axis to form a triangular shape of the traps and induce the rotation. Alternatively, the particles horizontal position (in addition to its vertical) is fixed by acoustics and we only need one optical trap to induce a rotation.
Indeed, we observe such a behaviour, as shown in Fig. 13. In this example we used a cancer spheroid. Cancer spheroids are used in Oncology to study the development of various cancer types and to screen the effect of therapeutic measures on the cellular level, in the present case the type used to study the oncogenesis of neuroblastoma at the Medical University of Innsbruck.69,70 Pointing a laser beam (marked by a red circle in Fig. 13) at the edge of the acoustically confined spheroid, induces a continuous rotation, which does not occur with acoustic trapping alone. This allows us to induce rotation for rather large specimens, for which it is more challenging to succeed with acoustics alone. Please note that the laser spot appears larger than its actual size in the images due to residual reflections. Typically the optically triggered rotation immediately stops when the tweezers are no longer overlapping with the sample. This indicates that the rotation is mainly induced by the optical radiation pressure, and the position of the tweezers on the sample determines the particle's orientation and rotation.
Fig. 13 Combined optical/acoustic rotation of spheroid: the spheroid did not rotate under acoustic trapping alone. By placing the HOT (marked by a red circle) on the spheroid, one can change its orientation. With the HOT placed at a specific location on the edge of the spheroid, a continuous rotation could be induced (ESI† Movie S8). |
When performing these particle manipulation experiments, where we want to manipulate a particle confined to a horizontal nodal plane by the top-transducer, we have observed an additional effect: placing laser beams, e.g., at the edge of a spruce (Picea abies) pollen grain induces a rotation and also a streaming of the fluid, as shown in Fig. 14. Here, the laser was split into three traps (marked by red circles) to use two of the spots to hold the particle and a third to rotate (traps in triangular shape). However, the horizontal acoustic force was sufficient to fix the position of the particle, and the three traps are in principle used as one trap in the examples included. To visualize the streaming in a single image, we apply the minimum intensity projection along the time axis, as described in sec. 2.4.2, to show the traces of small test particles. Both rotation and streaming are absent without the laser, see Fig. 14a). The streaming persists independent of whether we point the laser at the particle or at some distance away from the particle. As apparent in Fig. 14c), the streaming within the observation plane points out radially from the laser beam position. Interestingly the rotation also persists.
Fig. 14 Optically induced acoustic streaming: a) an acoustically trapped pollen grain from Picea abies is almost stationary under the given acoustic trapping parameters (ESI† Movie S9, ‘HOT off’). b) Placing the HOT (three spots marked by red circles) at the edge of the sample induced a continuous rotation of the sample (ESI† Movie S9, ‘HOT on’). c) The continuous rotation prevailed when the trap was moved far away from the sample (ESI† Movie S10). An identical treatment (minimum intensity projections of 220 frames) was applied to each of the three movie segments, revealing that outward streaming from the HOT center is present in b) and c), but not in a). |
In this case apparently the induced fluid flow leads to the sustained rotation of the particle via viscous drag, at least indirectly: it could be that the drag torque in combination with the acoustic viscous torque overcomes the acoustic radiation torque acting to lock the particle in a fixed orientation. This effect of inducing a particle rotation by a flow is similar to what happens in ‘thermo-optical tweezers',71 which utilize the convective flow induced by heating the fluid by light absorption at gold-coated cover slips to rotate the single cells around a transverse axis. However, in our case, the flow is not induced by convection due to laser heating, because it is strongly reduced if we keep the laser on but switch off the acoustic field. Nevertheless we suppose that a local temperature rise due to light absorption in the fluid plays a significant role. If we exchange the fluid by heavy water, which exhibits a smaller absorption of the near-infrared laser light, we observe a significant reduction in the optically induce streaming. We believe that the local temperature rise leads to a local modification of the fluid properties, e.g., viscosity, that triggers an enhanced acoustic streaming, similar to what happens at a surface. The exact mechanism is yet unexplored and further quantitative investigations are needed.
Specifically, we have demonstrated the effectiveness of a strategy for aligning or spinning biological samples in a controlled way: by tuning the relative strengths of the transducers (a strong optically transparent top-transducer for levitation and two weaker side-transducers, all operating at MHz frequencies) one can affect the relative strength of the orientation-angle dependent acoustic radiation torque Trad arising from pressure gradients and the acoustic viscous torque Tvisc arising from the shear forces at the viscous boundary layer of the particle. The former locks an asymmetric object into an orientation corresponding to lowest potential energy in the acoustic potential. One can determine the orientation direction the particle will settle into by tailoring the potential in 3D via the resonance modes excited by the active transducers.
For continuous rotation of the over-damped particle in the liquid, the acoustic viscous torque Tvisc has to overcome the acoustic radiation torque Trad. We found that one can facilitate rotation by choosing appropriate relative strengths of top- and side-transducer(s) to create a largely isotropic potential around the trapping location of the particle, since this reduces Trad. The viscous torque Tvisc driving the continuous rotation of the particle is created by the uptake of angular momentum by the viscous liquid around the particle. Our simulations have revealed the optimal conditions for creating angular-momentum carrying trapping nodes when exciting neighboring resonances of overlapping widths. This kind of optimization, however, requires detailed information of the 3D mode patterns in the resonators, which is often not directly available. In order to get at least a qualitative impression of the trapping nodes and of the streaming in the resonance chamber, we have developed some practical techniques utilizing suspensions of swimming micro-organism and visualizing their traces in an approach inspired by the ‘minimum intensity projection’ in medical imaging. Moreover, we have shown that – whenever the condition that the viscous torque exceeds the acoustic radiation torque is not (easily) achievable, for instance for larger and/or more asymmetric particles – auxiliary optical tweezers are suitable for inducing rotation, either by direct optical radiation effects or by optically induced streaming.
We have successfully performed sono-optical manipulation of various types of samples (starch particles, pollen grains, uni-cellular organisms, cell clusters, cancer spheroids and more) of sizes up to 350 μm, using forces per particle volume roughly around 0.1 pN μm−3, and achievable rotation rates of 0.1 Hz to 2 Hz. Controlled reorientation in all 3 directions was demonstrated on an elongated starch granule, and then sustained rotation, where the rotation axis and rate are controlled by the relative driving strength of the transducers, on a ∼125 μm pollen grain. In another example, Tvisc being only slightly larger than Trad led to highly non-uniform rotations rates, with the particle spending a much longer time in one orientation than in others. Tracking a small feature on the particle in a recorded sinogram and fitting the curve to an orbiting motion with an angle-dependent and a uniform contribution provided a way to estimate the ratio of the two torque contributions.
The next step towards long-term biological monitoring will be to adapt the device to the respective requirements (such as temperature control, sterility, options for gas exchange), and afterwards dedicated biocompatibility assays will have to be carried out. But since the pressure amplitudes and forces we are applying do not exceed the commonly used levels, we anticipate no adverse effects.12,14 We believe that our novel sono-optical approach will greatly facilitate the qualitative and quantitative assessments of suspended cell clusters, cancer spheroids and organoids, and possibly make a contribution towards reducing the necessary number of screening experiments on animals.
Footnotes |
† Electronic supplementary information (ESI) available: Including Movies S1–S10. See DOI: 10.1039/d0lc01261k |
‡ Please note that here T refers to a non-zero component of the torque in some specific direction, the full angular dependence of the torque vector in 3D is not addressed at the moment. |
This journal is © The Royal Society of Chemistry 2021 |