收藏切换
Self-Navigated 3D Acoustic Tweezers in Complex Media Based on Time Reversal
收藏切换
PDF
Ye Yang1, 2, Teng Ma1, 2, Sinan Li3, Qi Zhang1, Jiqing Huang1, Yifei Liu1, 2, Jianwei Zhuang1, Yongchuan Li1, Xuemin Du1, 2, Lili Niu1, 2, Yang Xiao1, 2, Congzhi Wang1, 2, Feiyan Cai1, 2, Hairong Zheng1, 2
Research. Vol 2021 Article ID 9781394
Less
收藏切换
Research. Vol 2021 Article ID 9781394
• Research Article •
Self-Navigated 3D Acoustic Tweezers in Complex Media Based on Time Reversal
Full
Ye Yang1, 2, Teng Ma1, 2, Sinan Li3, Qi Zhang1, Jiqing Huang1, Yifei Liu1, 2, Jianwei Zhuang1, Yongchuan Li1, Xuemin Du1, 2, Lili Niu1, 2, Yang Xiao1, 2, Congzhi Wang1, 2, Feiyan Cai1, 2, Hairong Zheng1, 2
Affiliations
  • 1Paul C. Lauterbur Research Center for Biomedical Imaging, Key Laboratory for Magnetic Resonance and Multimodality Imaging of Guangdong Province, Shenzhen Key Laboratory of Ultrasound Imaging and Therapy, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
  • 2Shenzhen College of Advanced Technology, University of the Chinese Academy of Sciences, Beijing 100049, China
  • 3Verasonics, Inc., WA 98034 USA
About Author:

Authors’ Contributions

T. Ma, F. Cai, and H. Zheng conceived the idea and designed the research methodology. T. Ma, J. Huang, Y. Li, and Q. Zhang designed and fabricated the device. Y. Yang and S. Li performed the experiment and analysed the data under the supervision of T. Ma and H. Zheng. X. Du fabricated the PDMS particles. L. Niu prepared the ex vivo macaque and human skulls. Y. Yang, T. Ma, S. Li, Y. Xiao, and C. Wang developed the time reversal acoustic tweezer algorithm. Y. Yang led the simulation with assistance from J. Zhuang. Y. Yang, T. Ma, and S. Li wrote the manuscript with assistance from Y. Li, and H. Zheng contributed to the revision. Ye Yang, Teng Ma, and Sinan Li contributed equally to this work.

Published: 2021-01-20 doi: 10.34133/2021/9781394
Outline
收藏切换

Acoustic tweezers have great application prospects because they allow noncontact and noninvasive manipulation of microparticles in a wide range of media. However, the nontransparency and heterogeneity of media in practical applications complicate particle trapping and manipulation. In this study, we designed a 1.04 MHz 256-element 2D matrix array for 3D acoustic tweezers to guide and monitor the entire process using real-time 3D ultrasonic images, thereby enabling acoustic manipulation in nontransparent media. Furthermore, we successfully performed dynamic 3D manipulations on multiple microparticles using multifoci and vortex traps. We achieved 3D particle manipulation in heterogeneous media (through resin baffle and ex vivo macaque and human skulls) by introducing a method based on the time reversal principle to correct the phase and amplitude distortions of the acoustic waves. Our results suggest cutting-edge applications of acoustic tweezers such as acoustical drug delivery, controlled micromachine transfer, and precise treatment.

Ye Yang, Teng Ma, Sinan Li, Qi Zhang, Jiqing Huang, Yifei Liu, Jianwei Zhuang, Yongchuan Li, Xuemin Du, Lili Niu, Yang Xiao, Congzhi Wang, Feiyan Cai, Hairong Zheng. Self-Navigated 3D Acoustic Tweezers in Complex Media Based on Time Reversal[J]. Research, 2021 , 2021 (1) : 9781394 - . DOI: 10.34133/2021/9781394
The trapping or transportation of micro- or nanoparticles is critical when attempting to catalyse a desired macroscopic effect via the manipulation of substances in the microscopic world. Optical tweezers [1, 2] are widely used tools for particle manipulation, and this has led to the development of a few cutting-edge applications in the field of biology and physics, such as evaluating the shear modulus of the human erythrocyte membrane [3], measuring the DNA spring constant [4], studying the dynamic characteristics of the kinesin movement [5, 6], measuring the instantaneous velocity of a Brownian particle [7], and trapping cold atoms [8]. Despite these noteworthy applications, optical tweezers suffer from a fundamental constraint in that light cannot traverse nor focus inside optically turbid media such as the human body. Therefore, optical tweezers can only be applied to a limited range of applications.
Acoustic tweezers recently emerged as an alternative tool for manipulating particles in a much broader range of media. Acoustic waves exert acoustic radiation forces on objects because of the momentum transfer that arises from sound scattering or attenuation. Acoustic radiation force has been demonstrated to manipulate several types of microparticles [9–11] in a noninvasive or noncontact manner. Because the velocity of acoustic waves is significantly smaller than that of light, the acoustic radiation force is typically five orders of magnitude greater than optical radiation force for a given input power; this enables the manipulation of more massive particles without increasing the risk of tissue ablation [12, 13]. Moreover, acoustic waves can propagate in transparent as well as optically turbid media irrespective of their optical refractive index, which considerably broadens their range of applications in comparison to that of optical tweezers (for example, in vivo particle manipulation [14, 15]). Currently, the manipulation of microbubbles is a typical application direction of acoustic tweezers, such as acoustic trapping and releasing of microbubbles in complex environments [16], transporting microbubbles using surface acoustic waves [17], and patterning of microbubbles [18].
Acoustic tweezers can be broadly classified into two categories [19], i.e., collective tweezers, which manipulate multiple particles undergoing the same behaviour [20], and selective tweezers, which create various traps that independently manipulate a subgroup of particles undergoing different behaviours [21–23]. Recently, artificial acoustic structures have been applied for complex beam generation and acoustic tweezers, such as waveguiding obtained by a superhydrophobic acoustic metasurface [24], generation of acoustic acceleration beams via flexible active surfaces [25], superoscillation wave packets using a metalens [26], realization of sound deceleration with helical-structured acoustic metamaterials [27], and fine acoustic manipulation via lossy metamaterials [28]. Although the acoustic fields obtained by artificial structures have a high spatial resolution, they are generally fixed and lack flexible dynamic characteristics. However, both collective and selective tweezers require a dynamic update of the acoustic wave fields. A phased array is preferred over single-element transducers because the amplitude and phase of the emitted signal from each element can be adjusted independently, thereby enabling arbitrary acoustic wave fields to be synthesized, updated, and moved in real time; this helps in enriching the types and modes of particle manipulation. For example, a controllable Bessel beam has been used for the precise manipulation of the mixture and stirring liquid microparticles [29]. Further, acoustic vortices have been reported to trap Mie particles in the air [30]. Meanwhile, the concept of holography has been applied to acoustic tweezers to levitate objects with twin, vortex, and bottle traps using single-sided arrays [21] and to dynamically manipulate multiple particles simultaneously and independently in a 3D space using a double-sided arrangement of two opposed arrays in the air [31]. Recently, particle manipulation which combines phased array and artificial acoustic structures has been proposed to integrate the advantages of the two methods [32].
While the 3D acoustic tweezers showed effective performance in the air [21, 23, 31, 33], for several reasons, it was rarely replicated in water, wherein the acoustic impedance is close to that of the human tissue. Currently, acoustic tweezers in water are mainly implemented using standing surface acoustic waves (SSAWs) [34, 35], which are commonly generated by interdigitated transducers. The SSAWs are typically operated at a considerably higher frequency range (tens to hundreds of megahertz [36]) to separate and manipulate tiny particles and cells. However, SSAWs encounter challenges such as the need for additional microfluidic systems, the limitations of the acoustic field patterns, and the small space for manipulating particles, especially in the vertical direction due to sound waves travel along the surface [37]. Acoustic streaming, which is generated by boundary-layer effects, was also used to manipulate microparticles in water environment, such as dynamic particle concentration and translation in an open fluid well [38] and acoustofluidic salivary exosome isolation [39]. To manipulate particles on a large scale in a 3D space, bulk acoustic waves generated by a matrix array need to be applied, which provide a critical basis for developing 3D acoustic tweezers in water. Sound attenuation drastically reduces in water; therefore, the operable frequency of the acoustic tweezers can be increased significantly to improve the precision of manipulation. This results in a considerably finer element pitch required for the matrix array to match the same trapping beam steerability as that of the airborne acoustic tweezers. Currently, the challenge exists in manufacturing an appropriate matrix array and a driving electronic system.
There are several more challenges for 3D acoustic tweezers in vivo. On the one hand, optical imaging cannot be applied to deep tissues, which makes it challenging to track the particle locations inside the human body—a prerequisite for acoustic tweezer beamforming. On the other hand, in vivo structures are complex and therefore require acoustic trapping beams carefully designed considering possible phase abbreviations or distortions caused by heterogeneous tissue structures such as human bones and skulls.
In this study, we aim to tackle these challenges by developing a self-navigated, time-reversal-based 3D acoustic tweezer system. The system consists of a 1.04 MHz, 256-element 2D matrix array; programmable 256-channel pulser-receiver system; and navigation system that dictates acoustic tweezer beamforming (Figure 1). The matrix array first transmitted imaging pulses. A volumetric ultrasound image was reconstructed with the received echo signals to indicate the current locations of the particles; these locations were then fed into the navigation system in real time to formulate the trapping beams to seize the particles (Figure 1, black loops in the middle). The trapping beams were formed based on the time reversal principle. For a complex medium, such as when a skull was present in the acoustic travel path, the navigation system first emulated received signals as if originated from the desired focal spots. It then controlled the ultrasound electronics to transmit the emulated waveforms from each element in the time-reversed order; thus, the trapping beams were focused at the desired locations by assuming the reciprocity of the sound travel paths (Figure 1, red loops). For homogeneous media such as water, beamforming is considerably simplified because the acoustic travel time from each element to the focal spot is known. The navigation system, in this case, simply calculated the appropriate phase delays from each element and passed them to the ultrasound transmission electronics (Figure 1, blue loops). Once the particles had been trapped at the initial location, the navigation system started to dictate the ultrasound electronics to manipulate the particles. The navigation system divided the desired path into a finite number of sections, and within each section, trapping beams—similar to the aforementioned process—were formed based on the time reversal principle. Imaging pulses were interleaved between the trapping beams, and they reported the real-time locations of particles as feedback to the navigation system. Thus, real-time imaging enabled the visualization and monitoring of the entire process of particle manipulation.
In this section, we presented the results of manipulating two different types of particles—polydimethylsiloxane (PDMS) and polystyrene (PS)—in water using the self-navigated 3D acoustic tweezer system (Figures 2(a)–2(i)). Particles made of PDMS and PS are considered “soft” and “rigid,” respectively, as determined by their relative acoustic impedance to water.
As shown in Figure 2(a), a PDMS particle fell from an arbitrary location above the 2D matrix array within its footprint and sank freely in water. The matrix array first transmitted imaging pulses to track the particle location (P0 to P1). Once the particle landed at the desired depth (P1), the navigation system immediately switched the ultrasound electronics to the particle-trapping mode, thereby forming a focused wave field at P1 to trap the particle temporally. Meanwhile, the navigation system was loaded in a user-defined path from the trapping position (P1) to the target position (P4), and it subsequently dictated the ultrasound electronics to manipulate the particles to move along the path. After arriving at the target position (P4), the particle was trapped there to levitate until the next command was given. This result showed that unlike conventional acoustic manipulation, particles need not be artificially placed near the traps initially. The volumetric ultrasound images—acquired using the coherent plane-wave compounding (CPWC) [40]—clearly display the position of the particle in 3D (Figure 2(b)), thereby automating the ultrasonic electronics for acoustic tweezer beamforming. To cross-validate the ultrasonic navigation, we acquired the optical images shown in Figures 2(c) and 2(d)) and Movie S1. The simulated trapping beam profiles at the nodes of the path (P1 to P4) are shown in Figure 2(e). Further, a PS particle with a diameter of 0.5 mm was trapped by the vortex acoustic field and transferred in a circular path by changing the centre of the vortex (Figures 2(f) and 2(g)). The simulated acoustic wave fields of the vortex are shown in Figures 2(h) and 2(i).
In this section, we showed the result of manipulating multiple particles in 3D space using the 2D matrix array. We started by manipulating two PDMS particles along an oblique elliptical path (Figure 3(a)); the particles were 1.5 mm in diameter. Figures 3(b) and 3(d) show snapshots of the particle movement projected in the x-y and y-z planes. Figures 3(c) and 3(e) are the simulated acoustic wave fields corresponding to Figures 3(b) and 3(d), respectively. We showed the collective manipulation of ten particles spinning around a circle (Figures 3(i) and 3(j)) and selective particle manipulation (Figures 3(f)–3(h)), wherein each of the ten particles was moved along the radial direction, transforming the circular pattern into a pattern of a plum blossom.
A large number of PS particles with approximate diameters of 400–800  μm were laid flat on the film to display the acoustic field pattern. The radiation force exerted on the PS particles pushed the particles away from the areas with intense acoustic pressure, thereby displaying acoustic field patterns in blank areas (Figures 3(k)–3(o)). The acoustic manipulation of two types of particles is shown in Movie S2.
Heterogeneity along acoustic travel paths defocuses ultrasound beams, making it challenging to implement acoustic tweezers in complex media. To mitigate this challenge, we applied the concept of time reversal in beamforming and restored beam focusing in heterogeneous media (see Materials and Methods). To validate this method, we made a complex-shaped baffle (Figure 4(b)), through which the focused beams formed assuming that the homogeneity of media (Figure 4(f)) were shown to defocus because of the phase abbreviation and distortion induced by the baffle (Figures 4(g)–4(i)). Shown in Figure 4(j) was the simulated signals originated from a target focal point located behind the baffle, received by our 256-element matrix array located on the other side of the baffle. Compared with Figure 4(e), where the baffle was absent, the waveforms and wavefronts in Figure 4(j) were distorted due to the heterogeneity. When we transmitted the signals back to the media in a time-reversed order (Figure 4(k)), the beams were refocused at the target focal spot (Figures 4(l)–4(n)), showing the effectiveness of the time-reversal-based beamforming in heterogeneous media. As an experimental validation, we showed in Figures 4(c) and 4(d) and Movie S3 that a PDMS particle placed behind the baffle was successfully trapped using the time-reversal-based method. The particle cannot be trapped with beams formed assuming homogeneity of media in the water.
Further, we implemented the time-reversal-based acoustic tweezing through an ex vivo human skull (Figures 5(c) and 5(d)). Shown in Figure 5(a)was the computed tomography (CT) reconstruction of the skull, which was input into our simulation that outputted the transmit waveforms we used in the experiment. In the experiment, we predefined a 3D manipulation path, which the particle was transferred along an isosceles triangular path in the y-z plane and then along a circular path in the x-y plane (Figure 5(a)). We also interleaved imaging pulses between the trapping beams and produced volumetric ultrasound images (Figure 5(b)) as a feedback of the particle manipulation process. Cross-validated by the optical imaging shown in Figures 5(c) and 5(d) and Movie S3, we manipulated a PDMS particle to move along the preset path through the skull; the distance between each movement step was 1 mm, and the time interval between steps was 0.5 second. Using the same approach, we also replicated the experiment through an ex vivo macaque skull (see Figure S10(a–d) and Movie S3).
In previous studies, 3D airborne array-based acoustic tweezers were developed for the manipulation of multiple particles. Three-dimensional underwater acoustic tweezers have a higher practical application value because the acoustic impedance of water is similar to that of human tissue. Further, the unique advantage of ultrasound in water is its imaging capability without the requirement of an optically transparent medium. Ultrasonic waves attenuate rapidly in air, which makes imaging by reflecting or transmitting signals difficult. The unique imaging ability of acoustic tweezers can be used for particle navigation during trapping. Usually, selective acoustic manipulation requires the artificial placement of particles near traps, which is exceedingly difficult, especially for cells, bacteria, or other small particles. With the help of ultrasonic imaging, freely falling PDMS particles can be trapped at any initial location and then transferred to the specified position. This offers a feasible and convenient method for selective acoustic manipulation, which is more suitable for practical applications. Moreover, by combining ultrasound imaging, the volumetric ultrasonic images of particle manipulation along the 3D path through the ex vivo human skull can be obtained (Figure 5(b)), which indicates that the 3D ultrasonic images can record the entire manipulation process in the heterogeneous medium. Despite resolution limitations constrained by the low driving frequency and structure of the transducer, we can also observe the motion state and position of the particle through the skull. These results demonstrated the possibility of using acoustic tweezers in in vivo applications.
In this study, the system’s response time from one frame of ultrasonic images to the next frame mainly includes the time for the transmission and reception of imaging pulses, transmission of manipulation pulses, data transfer, and image reconstruction and processing. This will greatly influence the tracking of particle manipulation. Here, we used the Verasonics Vantage 256 system (Verasonics Inc., WA, USA) as the electronic system, which has a sufficiently high system clock frequency (250 MHz) and pulse repetition rate (up to 10 kHz). First, during the transmission and reception of imaging pulses, as shown in Figure S1, four images can be obtained by transmitting 16 imaging pulses (1 cycle each pulse) with an interval of 327  μs (the total transmitting time was approximately 5.2 ms). This time interval is determined by the maximum imaging depth (132 mm) and the acoustic velocity, and it is independent of the working frequency of the array. Second, during the transmission of manipulation pulses, in one manipulation sequence, ninety trapping pulses (5 cycles each pulse) were transmitted with an interval of 360  μs (the total time was approximately 32.8 ms). This time interval can be adjusted for different moving steps and speed requirements under the stable manipulation condition. Third, the time required for data transfer, as well as image reconstruction and processing, largely depends on the working frequency of the array, system hardware condition, and imaging method and scope. A higher working frequency can generate smaller focal spots, making it possible to manipulate smaller particles; however, this requires a high system sampling rate, which may increase the size of the stored data, thus reducing the speed of data transfer and image reconstruction and processing. The working frequency of the array was 1.04 MHz, and the sampling rate was 4.16 MHz (quadruple sampling). When the working frequency is increased, for example, by five times, the sampling frequency should also be set to increase by five times to ensure the sampling accuracy, which increases the length of the radio frequency data to five times the original length. However, the attenuation of ultrasound will become larger with the increase in the working frequency. A smaller maximum imaging depth is usually set for higher frequencies, which reduces the number of pixels in the imaging space. Therefore, the total time of data transfer and image reconstruction and processing for high working frequencies depends on both the sampling frequency and imaging depth. The hardware condition such as the RAM memory and CPU or GPU speed of the electronic system can also be upgraded to improve the total time. The expansion of RAM space is helpful for decreasing the time required to read and store data in the real-time imaging process. The use of GPU or CPU parallel computing can improve the speed of imaging reconstruction. Further, the imaging algorithm can be optimized by reducing the numbers of plane wave compounding, and the scope of imaging acquisition can also offer some improvements, thus reducing the transmitting and receiving time of imaging, as well as the time for data transmission and reconstruction.
As shown in the results, when manipulating ten PDMS particles in a plane, the distance and time steps between the movements of each particle were approximately 1 mm and 1 s, respectively. The minimum distance between particles was approximately 5 mm. When the distance between two particles was smaller, the focal spots merged and hindered independent manipulation. The Rayleigh limit [41] for this configuration was approximately 2.5 mm (1.22λ F/A, where λ =1.5 mm, F = 60 mm, and A=44.8 mm denote the wavelength, focal distance, and aperture of the array, respectively); however, the ideal distance between traps can only be achieved with a smaller array pitch (≤0.5λ). Moreover, we studied the relationship between the acoustic radiation force exerted on PDMS particles and the number of focal spots and elements, respectively. The simulations indicated that the maximum acoustic radiation force was inversely proportional to the number of focal points (Figure S4)), and it can be increased by reducing the pitch of the array (Figure S5). However, no further improvements were observed when the pitch of the array was smaller than λ/2. Both the axial and transverse radiation forces followed the same trend; however, the axial forces were smaller because of the smaller axial acoustic field gradient.
We designed a low-frequency (1.04 MHz) array with a large pitch (2λ) for the experiments to achieve an array with a large aperture (30λ) that can provide a large manipulation space and trapping force, thereby enabling in vivo applications. The size of the controllable particles is related to the frequency of the array; based on the low-frequency array and size of the traps (the diameters of the focal spot and vortex centre were approximately 3 and 1 mm, respectively), we selected the PDMS and PS particles with a millimetre scale for manipulation. The results also demonstrated that the trapping force and Rayleigh limit can be improved by decreasing the pitch of the array; the pitch should be reduced to near or less than λ/2 if possible. However, when the aperture is invariant, a smaller pitch means more array elements, which imposes new challenges to the design and manufacture of arrays and the support of excitation systems. In this study, under the premise of obtaining sufficient trapping force and manipulating space, we selected a pitch of 2λ to overcome the trade-off between the small pitch and the existing excitation systems and the manufacturing technology of the array. The proposed method of single-array acoustic manipulation was found to be superior compared to array-based acoustic tweezers that use the standing wave methods in terms of setup flexibility and cost reduction. Further, our method prevented the appearance of secondary nodes along the z-axis, which resolved a pressing problem in holographic acoustic tweezers [31]. However, there were limitations to our method. Particle manipulation in the beam propagation direction was less accurate and stable because the axial restraint only depended on the balance of gravity, buoyancy, and acoustic radiation force, while the axial radiation force was exerted upwards instead of converging at the trapping position. The instability of motion was also evident in the heterogeneous medium, where the particles floated upwards and downwards in the axial direction during movement. The distances from the transducer plane to the stably trapped particles are related to the preset focal distance and the energy exerted on the particles. In this study, the preset focal depth from 4 cm to 12 cm (beam width 2.1 mm to 3.3 mm) was selected for trapping to reduce the impact of the grating lobes and to maintain the good focusing ability of the array (Figure S6). The relationships between the driving voltage, preset focal distance, and actual trapping position were explored with particles of different sizes (0.61λ, 0.93λ, and 1.41λ), respectively. As shown in Figure S7, when the focal distance is constant, the trapped particles move upward with an increasing voltage. Once the driving voltage is increased to a certain extent, the particles begin to vibrate and cannot be trapped stably. By adjusting the combination of the preset focal distance and the driving voltage, the particles may be trapped more accurately at a specified height. The results also indicated that the maximum and the minimum distances from the transducer plane to the particles that can be trapped were 17.1 cm and 3.7 cm, respectively. The relationships among the duty factor, pulse repetition frequency (PRF), and position of the manipulated particle had also been studied. The duty factor used in this paper was 1.3%, and the PRF was 2.7 kHz. The influence of the duty factor on manipulation is mainly reflected in the energy that is exerted on the manipulated particles. When the PRF is constant, the particles are driven to move up with the increasing duty factor owing to the increased energy (Figure S8). However, when the duty factor is constant, increasing PRF has little influence on the trapping position (Figure S9).
In this study, a self-designed 1.04 MHz, 256-element 2D matrix array was fabricated for the underwater 3D acoustic tweezers. Its ability of 3D ultrasonic imaging in water was successfully used for guiding and monitoring the process of particle acoustic manipulation in real time. The 3D images tracked a freely falling soft particle from an arbitrary initial position and guided the system to trap and manipulate the particle according to the designed route. Moreover, we dynamically manipulated multiple relatively soft and rigid microparticles in 3D using different traps (multifoci and vortex traps). Furthermore, we introduced the time-reversal-based method to correct the phase and amplitude distortion caused by a heterogeneous medium and successfully manipulated the soft particles through a SIAT baffle and ex vivo macaque and human skulls. These results preliminarily address the existing challenges in the in vivo applications of acoustic tweezers and facilitate their development for practical biomedical applications. Future studies will be performed by exploring new acoustic traps and using arrays with higher frequencies, higher number of elements, smaller pitch, or higher transmitting output power, which will enable the manipulation of considerably smaller particles or cells and improve the capabilities of acoustic tweezers.
As shown in Figure 1, we realized 3D acoustic tweezers in water using a self-fabricated 1.04 MHz, 256-element (16×16) ultrasonic 2D matrix array with a pitch of 2.8 mm and kerf of 0.2 mm. Each element contained matching layers, piezocomposites, Cr/Au electrodes, flexible circuits, and backing layers. The 3D acoustic tweezer system was a programmable 256-channel pulser-receiver system designed by Verasonics Vantage 256 system. The system was used for exciting arrays using trapping and imaging pulse signals to obtain specified acoustic fields and for receiving echo signals to reconstruct images. The phase, amplitude, and waveform of a signal exerted on each element can be adjusted dynamically and independently. The navigation system referred to the high-performance workstation T5810 (Dell Inc., TX, USA) capable of running simulation programs, which dictated acoustic tweezer beamforming by calculating the specified pulse signals and guiding the 3D acoustic tweezer system to excite the array according to the signals.
The matrix array was first driven by imaging pulses. The echo signals from the particle were reconstructed for the 3D ultrasound images to indicate the current location of the particle (Figure 1, P0 position). The navigation system dictated the trapping beams to seize particles by importing the position information in real time (Figure 1, black loops in the middle). The imaging sequences were interleaved between the manipulation sequences (see Figure S1 and Note S3 for details regarding sequence arrangement). The on-time duration of the half cycle and the number of cycles for each pulse signal in these two sequences were different; a higher on-time duration (0.9) and greater number of cycles (five cycles) provide higher energy in the manipulation sequence, while a lower on-time duration (0.8) and smaller number of cycles (one cycle) provide higher resolution in the imaging sequence. The transmitting signals of each element in the manipulation sequence were different, and the generated sound waves converged to form the target multifoci acoustic field. When trapped, the particle can be transferred along a user-defined path by changing the position of the focal spots. The path was divided into a finite number of steps by the navigation system, and the distance of each step of the particle movement was approximately 0.5 mm. The CPWC method was used for the imaging sequence, which improved the image quality of the single plane wave imaging mode [40]. Four plane waves with different inclination angles were transmitted, and each frame of the 3D image was obtained by coherently compounding the images obtained with the waves. The entire 3D imaging space of the self-navigated 3D acoustic tweezers is shown in Figure S2.
In a homogeneous medium, the navigation system calculated the appropriate phase delays using the acoustic travel time, which is equal to the distance between each element and focal spot divided by the wave velocity, and it then passed them to the 3D acoustic tweezer system. When the particle was trapped, the navigation system can divide the path into a finite number of steps by inputting the motion trajectory information (Figure 1, from P0 through P1 to P2, blue dashed line); each step can update the trapping beams to manipulate the particle along the path (Figure 1, blue loops).
Particles with different acoustic contrast factors exist in water [42]. Here, we used the PDMS and PS particles to demonstrate the capability to trap particles with negative and positive contrast factors, respectively. To manipulate the particles, we implemented multifoci (Figure S3(a)) and vortex (Figure S3(d)>) acoustic traps, respectively. The directions of the acoustic radiation force exerted on the PDMS (Figure S3(b–c)) and PS (Figure Figure S3(e–f)) particles were opposite (Note S1). The radiation force exerted on the PDMS particles at multiple focal spots and that exerted on the PS particles at the centre of the vortex converges in the lateral directions and directs upward in the axial direction, thereby enabling the trapping and manipulation of the particles. The transmission signals of each element were determined using the pseudoinverse algorithm (PINV) [43, 44] to generate the multifoci field by initializing the position of the focal spots (Note S2). The vortex was realized by adjusting the transmission phase of each element to generate a screw dislocation in the wavefront [45, 46] (Note S2). When particles were trapped at the focal spots or the vortex centre, the transmitting signals were dynamically changed to move the spots or vortex and hence, the particles. In this study, we manipulated the particles initially resting on a film, thereby enabling the acoustic waves to propagate while suffering little reflection and maintaining travelling-wave character.
We manipulated two PDMS particles in a 3D space simultaneously, whose motion trajectory was preset as a spatial ellipse (the angle between the plane of motion and z-axis was 45° and the diameter of the circular path in the top view was 12 mm). Ten particles were manipulated to transform a circle (diameter=20 mm) to a pentagon (outer circle diameter=20 mm) and rotate the circle in the clockwise and counterclockwise directions. One PS particle (diameter = 500  μm) was trapped to move in a circular path (diameter=10 mm). A large number of PS particles (diameter=400–800  μm) were manipulated to display the patterns of the vortex (topological charge=1), pentagon (40 closely arranged focal spots, outer circle diameter=24 mm), and snowflake (48 closely arranged focal spots, outer circle diameter=24 mm) acoustic fields.
To achieve acoustic manipulation in heterogeneous media, the trapping beams should be carefully designed by considering the possible phase abbreviations or distortions caused by the media. The navigation system can simulate the wave propagation process by importing the medium information (Figure 1, red dashed line) and perform beamforming using signals generated by the time-reversal-based method. This enables the trapping beams to focus at the target position (Figure 1, red loops).
The time-reversal-based method was used for the phase and amplitude correction [47]. The steps in the method are as follows: (1) use the array and heterogeneous medium information (e.g., density, acoustic velocity, and geometry) for acoustic wave propagation modelling based on the navigation system and simulate ultrasound wave signals which propagate through the medium by setting up a virtual sound source at the target trapping position; (2) reverse the simulated signals received by each element along the timeline, and divide them by the amplitude correction coefficient which is equal to the maximum signal amplitude of each element divided by the maximum signal amplitude of all elements; and (3) obtain the transmitting signals by extracting the phase, amplitude, and waveform information of the reversed and modified simulated signals. When the array was driven by the transmitting signals, the acoustic waves propagated through the heterogeneous medium and refocused at the position of the original virtual sound source. Then, the PDMS particles were trapped at the focal spot in the heterogeneous medium.
In this study, we successively added a SIAT baffle, an ex vivo macaque (Figure S10(a)), and an ex vivo human skull between the array and the PDMS particle to demonstrate the heterogeneous medium. The SIAT baffle was 3D-printed with photosensitive resin (Figure 4(b)). The SIAT baffle had the following dimensions: 80 mm×80 mm×11.2 mm. Wave stripes were positioned perpendicular to each other on the upper and lower surfaces (the radius of curvature of each wave stripe was 7.07 mm, and the central angle was 90°), and SIAT grooves were engraved on the wave stripes on the top surface (the bottom of the groove was 5.6 mm away from the top surface) to add more irregularities to the baffle. The sound velocity and density of the baffle were measured to be 2368 m/s and 1123 kg/m3, respectively. The absorption parameter was measured to be 7.3 dB/cm at 1.04 MHz. The designed SIAT baffle can be directly modelled in the simulation by importing information related to its geometry and the acoustic parameters. The acoustic pressure maps obtained through the SIAT baffle with and without applying the time-reversal-based method were simulated and measured to explore the effectiveness of the method.
The macaque and human skulls were obtained from the corpses of a macaque and human, respectively, following Chinese ethical laws. The macaque skull was cut transversely, whereas the human skull was cut longitudinally and opened using a saw. The skulls were prepared by adopting standard anatomical procedures; they were cleaned and disinfected subsequently. The ex vivo macaque skull had dimensions of 65 mm×68 mm×38 mm and an irregular quarter spherical shell shape with an average thickness of approximately 1.5 mm. We used the right half of the ex vivo human skull with dimensions of 19 cm×14 cm×8 cm and average thickness of approximately 12 mm; this was approximately eight times thicker than the ex vivo macaque skull. Before the experiment, the skulls were immersed in degassed water for over 48 h. The experiments were performed in the degassed water at room temperature.
k-Wave, a MATLAB toolbox for the time-domain simulation of acoustic wave fields, was used to perform all simulations in the heterogeneous medium [48]. The core advantage of k-Wave is its advanced numerical model that can account for the arbitrary distribution of heterogeneous medium parameters along with linear and nonlinear wave propagation.
In the case of acoustic manipulation through the skulls, because the geometry of the skulls was not designed by us, we used CT to scan the entire setup of the experiment (including the array, skulls, and bracket); next, we imported the obtained relative position information of the skulls into the simulation. The sound velocity and density information at each simulated grid point of the skulls was obtained by mapping the CT values of the corresponding pixels in the CT images [49]. When excluding the unnecessary parts, we set the minimum CT value corresponding to the sound velocity and density (3100 m/s and 2200 kg/m3, respectively) of the skull and the maximum CT value corresponding to the sound velocity and density (1540 m/s and 1000 kg/m3, respectively) of water. At the simulation points where the CT values varied between the maximum and minimum values, the sound velocity and density information can be obtained by linearly corresponding to the sound velocity (1540–3100 m/s) and density intervals (1000–2200 kg/m3), respectively.
The 256-element (16×16) 2D matrix array was fabricated using the conventional transducer technology. A piece of PZT-5H bulk ceramic was used to prepare the piezoelectric 1-3 composite based on a dice-and-fill technique. In this study, a 170  μm thick dicing blade was selected to achieve a kerf of 200  μm and a pitch of 700  μm. Then, the prepared PZT-5H 1-3 composite was lapped down to the designed thickness based on the PiezoCAD simulation result. Both polished surfaces were sputtered with a chrome/gold (Cr/Au) electrode using a sputtering system with the designed thicknesses of 200 nm/500 nm. The front side of the PZT-5H composite was cast with a matching material made of aluminium oxide powder (23  μm) and Epo-Tek 301 epoxy (Epoxy Technologies, MA, USA). After curing at room temperature for 24 h, the matching layer was polished to an optimized thickness of 750  μm. The prepared acoustic stack with the matching layer was diced into the designed dimensions of 44.8 mm×44.8 mm, which consisted of 256 array elements. Then, the back side of the PZT-5H 1-3 composite was diced with an electrode division kerf of 200  μm, a depth of 200  μm, and an element width of 2.6 mm×2.6 mm for the array assembly process. Four pieces of customized flexible circuits fabricated using polyimide with 64 electrode traces (trace width=254  μm and thickness=118  μm) were soldered to the array elements. Then, the array elements attached to the flexible circuits were fixed inside an acrylic housing, and the gap between the acoustic stack and housing was filled with an insulating epoxy. A total of 256 coaxial cables were sealed with the flexible circuit for the further characterization and application of acoustic tweezers. Finally, a backing material made of Epo-Tek 301 epoxy with 3–5  μm aluminium nitride powder was filled into the housing to reduce the ring-down effect.
The employed particles were 0.9–3 mm diameter PDMS spheres and 0.4–1 mm diameter PS spheres (Huge Biotechnology Co., Ltd., Shanghai, China). The PDMS elastomer kit (Sylgard 184) was purchased from Dow Corning, and Span 80 and hexane were obtained from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). The uniform PDMS microspheres were prepared using a modified mechanical emulsification method [50–52]. The prepolymer base elastomer and the curing agent were mixed with a weight ratio of 10 : 1 under intense stirring. Then, the obtained mixtures were diluted by hexane with a volume ratio of 1 : 9, followed by mechanical stirring and vacuum pumping for 15 min. Finally, the mixture was transferred to a syringe (10  μl). Subsequently, the resulting solution was slowly injected into an aqueous solution (10 ml) that contained the surfactant, Span 80 (5% w/w), without any agitation. Finally, the obtained PDMS microspheres were cured at room temperature for 24 h.
  • Shenzhen Key Laboratory of ultrasound imaging and therapy(ZDSYS20180206180631473)
  • Chinese Academy of Sciences(YJKYYQ20190077)
  • Shenzhen Peacock Plan(KQTD20170810160424889)
  • National Natural Science Foundation of China(81827802; 81827807; 11774370; 61571431; 81527901)
1
Ashkin A., Dziedzic J. M., Bjorkholm J. E., and Chu S., “Observation of a single-beam gradient force optical trap for dielectric particles,” Optics Letters, vol. 11, no. 5, p. 288, 1986
2
Ashkin A., and Dziedzic J. M., “Optical trapping and manipulation of viruses and bacteria,” Science, vol. 235, no. 4795, pp. 1517–1520, 1987
3
Henon S., Lenormand G., Richert A., and Gallet F., “A new determination of the shear modulus of the human erythrocyte membrane using optical tweezers,” Biophysical Journal, vol. 76, no. 2, pp. 1145–1151, 1999
4
Wang M. D., Yin H., Landick R., Gelles J., and Block S. M., “Stretching DNA with optical tweezers,” Biophysical Journal, vol. 72, no. 3, pp. 1335–1346, 1997
5
Svoboda K., Schmidt C. F., Schnapp B. J., and Block S. M., “Direct observation of kinesin stepping by optical trapping interferometry,” Nature, vol. 365, no. 6448, pp. 721–727, 1993
6
Yildiz A., Tomishige M., Vale R. D., and Selvin P. R., “Kinesin walks hand-over-hand,” Science, vol. 303, no. 5658, pp. 676–678, 2004
7
Li T., Kheifets S., Medellin D., and Raizen M. G., “Measurement of the instantaneous velocity of a Brownian particle,” Science, vol. 328, no. 5986, pp. 1673–1675, 2010
8
Ketterle W., Davis K. B., Joffe M. A., Martin A., and Pritchard D. E., “High densities of cold atoms in a dark spontaneous-force optical trap,” Physical Review Letters, vol. 70, no. 15, pp. 2253–2256, 1993
9
Wu J. R., “Acoustical tweezers,” The Journal of the Acoustical Society of America, vol. 89, no. 5, pp. 2140–2143, 1991
10
Brandt E. H., “Acoustic physics. Suspended by sound,” Nature, vol. 413, no. 6855, pp. 474–475, 2001
11
Bruus H., “Acoustofluidics 7: the acoustic radiation force on small particles,” Lab on a Chip, vol. 12, no. 6, pp. 1014–1021, 2012
12
Lee J., Teh S. Y., Lee A., Kim H. H., Lee C., and Shung K. K., “Single beam acoustic trapping,” Applied Physics Letters, vol. 95, no. 7, article 073701, 2009
13
Lam K. H., Hsu H. S., Li Y., Lee C., Lin A., Zhou Q., Kim E. S., and Shung K. K., “Ultrahigh frequency lensless ultrasonic transducers for acoustic tweezers application,” Biotechnology and Bioengineering, vol. 110, no. 3, pp. 881–886, 2013
14
Shih-Tsung Kang, and Chih-Kuang Yeh, “Potential-well model in acoustic tweezers,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 57, no. 6, pp. 1451–1459, 2010
15
Li Y., Lee C., Chen R., Zhou Q., and Shung K. K., “A feasibility study ofin vivoapplications of single beam acoustic tweezers,” Applied Physics Letters, vol. 105, no. 17, p. 173701, 2014
16
Baresch D., and Garbin V., “Acoustic trapping of microbubbles in complex environments and controlled payload release,” Proceedings of the National Academy of Sciences, vol. 117, no. 27, pp. 15490–15496, 2020
17
O’Rorke R., Wood C. D., Wälti C., Evans S. D., Davies A. G., and Cunningham J. E., “Acousto-microfluidics: transporting microbubble and microparticle arrays in acoustic traps using surface acoustic waves,” Journal of Applied Physics, vol. 111, no. 9, article 094911, 2012
18
Bernassau A., MacPherson P., Beeley J., Drinkwater B., and Cumming D., “Patterning of microspheres and microbubbles in an acoustic tweezers,” Biomedical Microdevices, vol. 15, no. 2, pp. 289–297, 2013
19
Baudoin M., and Thomas J. L., “Acoustic tweezers for particle and fluid micromanipulation,” Annual Review of Fluid Mechanics, vol. 52, no. 1, pp. 205–234, 2020
20
Courtney C. R. P., Ong C. K., Drinkwater B. W., Bernassau A. L., Wilcox P. D., and Cumming D. R. S., “Manipulation of particles in two dimensions using phase controllable ultrasonic standing waves,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 468, no. 2138, pp. 337–360, 2012
21
Marzo A., Seah S. A., Drinkwater B. W., Sahoo D. R., Long B., and Subramanian S., “Holographic acoustic elements for manipulation of levitated objects,” Nature Communications, vol. 6, no. 1, 2015
22
Baresch D., Thomas J. L., and Marchiano R., “Observation of a single-beam gradient force acoustical trap for elastic particles: acoustical tweezers,” Physical Review Letters, vol. 116, no. 2, article 024301, 2016
23
Prisbrey M., and Raeymaekers B., “Ultrasound noncontact particle manipulation of three-dimensional dynamic user-specified patterns of particles in air,” Physical Review Applied, vol. 10, no. 3, article 034066, 2018
24
Tong L., Xiong Z., Shen Y. X., Peng Y. G., Huang X. Y., Ye L., Tang M., Cai F. Y., Zheng H. R., Xu J. B., Cheng G. J., and Zhu X. F., “An acoustic meta-skin insulator,” Advanced Materials, vol. 32, no. 37, article 2002251, 2020
25
Li P. Q., Shen Y. X., Geng Z. G., Cao P. C., Peng Y. G., and Zhu X. F., “Realization of controllable acoustic acceleration beams via flexible active surfaces,” Journal of Physics D-Applied Physics, vol. 53, no. 15, p. 155502, 2020
26
Shen Y. X., Peng Y. G., Cai F., Huang K., Zhao D. G., Qiu C. W., Zheng H., and Zhu X. F., “Ultrasonic super-oscillation wave-packets with an acoustic meta-lens,” Nature Communications, vol. 10, no. 1, p. 3411, 2019
27
Zhu X. F., Li K., Zhang P., Zhu J., Zhang J., Tian C., and Liu S., “Implementation of dispersion-free slow acoustic wave propagation and phase engineering with helical-structured metamaterials,” Nature Communications, vol. 7, no. 1, 2016
28
Zhu Y., Hu J., Fan X., Yang J., Liang B., Zhu X., and Cheng J., “Fine manipulation of sound via lossy metamaterials with independent and arbitrary reflection amplitude and phase,” Nature Communications, vol. 9, no. 1, p. 1632, 2018
29
Hong Z., Zhang J., and Drinkwater B. W., “Observation of orbital angular momentum transfer from Bessel-shaped acoustic vortices to diphasic liquid-microparticle mixtures,” Physical Review Letters, vol. 114, no. 21, p. 214301, 2015
30
Marzo A., Caleap M., and Drinkwater B. W., “Acoustic virtual vortices with tunable orbital angular momentum for trapping of Mie particles,” Physical Review Letters, vol. 120, no. 4, article 044301, 2018
31
Marzo A., and Drinkwater B. W., “Holographic acoustic tweezers,” Proceedings of the National Academy of Sciences, vol. 116, no. 1, pp. 84–89, 2019
32
Cox L., Melde K., Croxford A., Fischer P., and Drinkwater B. W., “Acoustic hologram enhanced phased arrays for ultrasonic particle manipulation,” Physical Review Applied, vol. 12, no. 6, 2019
33
Hoshi T., Ochiai Y., and Rekimoto J., “Three-dimensional noncontact manipulation by opposite ultrasonic phased arrays,” Japanese Journal of Applied Physics, vol. 53, no. 7S, article 07KE07, 2014
34
Guo F., Mao Z., Chen Y., Xie Z., Lata J. P., Li P., Ren L., Liu J., Yang J., Dao M., Suresh S., and Huang T. J., “Three-dimensional manipulation of single cells using surface acoustic waves,” Proceedings of the National Academy of Sciences, vol. 113, no. 6, pp. 1522–1527, 2016
35
Chen K. J., Wu M., Guo F., Li P., Chan C. Y., Mao Z., Li S., Ren L., Zhang R., and Huang T. J., “Rapid formation of size-controllable multicellular spheroids via 3D acoustic tweezers,” Lab on a Chip, vol. 16, no. 14, pp. 2636–2643, 2016
36
Meng L., Cai F., Li F., Zhou W., Niu L., and Zheng H., “Acoustic tweezers,” Journal of Physics D: Applied Physics, vol. 52, no. 27, article 273001, 2019
37
Ozcelik A., Rufo J., Guo F., Gu Y., Li P., Lata J., and Huang T. J., “Acoustic tweezers for the life sciences,” Nature Methods, vol. 15, no. 12, pp. 1021–1028, 2018
38
Bachman H., Gu Y., Rufo J., Yang S., Tian Z., Huang P. H., Yu L., and Huang T. J., “Low-frequency flexural wave based microparticle manipulation,” Lab on a Chip, vol. 20, no. 7, pp. 1281–1289, 2020
39
Wang Z., Li F., Rufo J., Chen C., Yang S., Li L., Zhang J., Cheng J., Kim Y., Wu M., Abemayor E., Tu M., Chia D., Spruce R., Batis N., Mehanna H., Wong D. T. W., and Huang T. J., “Acoustofluidic salivary exosome isolation: a liquid biopsy compatible approach for human papillomavirus–associated oropharyngeal cancer detection,” The Journal of Molecular Diagnostics, vol. 22, no. 1, pp. 50–59, 2020
40
Montaldo G., Tanter M., Bercoff J., Benech N., and Fink M., “Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 56, no. 3, pp. 489–506, 2009
41
Rayleigh L., “XXXI. Investigations in optics, with special reference to the spectroscope,” The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, vol. 8, no. 49, pp. 261–274, 2009
42
Johnson L. M., Gao L., Shields IV C., Smith M., Efimenko K., Cushing K., Genzer J., and López G. P., “Elastomeric microparticles for acoustic mediated bioseparations,” Journal of Nanobiotechnology, vol. 11, no. 1, p. 22, 2013
43
Ebbini E. S., and Cain C. A., “Multiple-focus ultrasound phased-array pattern synthesis: optimal driving-signal distributions for hyperthermia,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 36, no. 5, pp. 540–548, 1989
44
Hertzberg Y., Naor O., Volovick A., and Shoham S., “Towards multifocal ultrasonic neural stimulation: pattern generation algorithms,” Journal of Neural Engineering, vol. 7, no. 5, article 056002, 2010
45
Hefner B. T., and Marston P. L., “An acoustical helicoidal wave transducer with applications for the alignment of ultrasonic and underwater systems,” The Journal of the Acoustical Society of America, vol. 106, no. 6, pp. 3313–3316, 1999
46
Shi C., Dubois M., Wang Y., and Zhang X., “High-speed acoustic communication by multiplexing orbital angular momentum,” Proceedings of the National Academy of Sciences, vol. 114, no. 28, pp. 7250–7253, 2017
47
Fink M., “Time reversal of ultrasonic fields. I. Basic principles,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 39, no. 5, pp. 555–566, 1992
48
Treeby B. E., and Cox B. T., “k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields,” Journal of Biomedical Optics, vol. 15, no. 2, article 021314, 2010
49
Jing Y., Meral F. C., and Clement G. T., “Time-reversal transcranial ultrasound beam focusing using a k-space method,” Physics in Medicine and Biology, vol. 57, no. 4, pp. 901–917, 2012
50
du X., Lei N. Y., Chiu H. M., Ge X., Zhang Z., and Hon-Wah Lam M., “Non-invasive in vivo imaging of the ionic regimes along the gastrointestinal tract of a freshwater vertebrate model organism (Japanese medaka) using responsive photonic crystal beads,” Journal of Materials Chemistry B, vol. 1, no. 11, pp. 1535–1540, 2013
51
du X., Lei N. Y., Hu P., Lei Z., Ong D. H. C., Ge X., Zhang Z., and Lam M. H. W., “In vivo imaging of the morphology and changes in pH along the gastrointestinal tract of Japanese medaka by photonic band-gap hydrogel microspheres,” Analytica Chimica Acta, vol. 787, pp. 193–202, 2013
52
Zhao Q., Cui H., Wang Y., and du X., “Microfluidic platforms toward rational material fabrication for biomedical applications,” Small, vol. 16, no. 9, article 1903798, 2019
Year 2021 volume 2021 Issue 1
PDF
57
32
Cite this Article
BibTeX
Article Info
doi: 10.34133/2021/9781394
  • Receive Date:2020-08-15
  • Online Date:2026-10-27
  • Published:2021-01-20
Article Data
Affiliations
History
  • Received:2020-08-15
  • Accepted:2020-10-08
Funding
Shenzhen Key Laboratory of ultrasound imaging and therapy(ZDSYS20180206180631473)
Chinese Academy of Sciences(YJKYYQ20190077)
Shenzhen Peacock Plan(KQTD20170810160424889)
National Natural Science Foundation of China(81827802; 81827807; 11774370; 61571431; 81527901)
Affiliations
    1Paul C. Lauterbur Research Center for Biomedical Imaging, Key Laboratory for Magnetic Resonance and Multimodality Imaging of Guangdong Province, Shenzhen Key Laboratory of Ultrasound Imaging and Therapy, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
    2Shenzhen College of Advanced Technology, University of the Chinese Academy of Sciences, Beijing 100049, China
    3Verasonics, Inc., WA 98034 USA
References
Share
https://castjournals.cast.org.cn/joweb/research/EN/10.34133/2021/9781394
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)
属
Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT