音響ピンセット

Manipulation of very small objects with sound waves

音響ピンセット(音響ピンセットとも呼ばれる)は、音波を使用して、直径 100ナノメートルから 10ミリメートルの非常に小さな物体の位置と動きを操作するツール セットです。この方法で浮上した物体の最大密度は 5.7 g/cm^3 です。物体を浮上させるのに使用される音は 20 kHz以上の範囲で、通常、ほとんどの消費者向けピンセットと浮上装置では 40 kHzが使用されます。

[1]厳密に言えば、単一ビームベースの構成のみが音響ピンセットと呼ばれます。ただし、音響ピンセットの広い概念には、2つのビーム構成が含まれます。1つは単一の音ビーム、もう1つは定在波を作成するための音の反射器、もう1つは互いに直接向けられた2つの音ビームです。この技術は、音圧の節と腹の位置と距離を制御することで機能します。[2]これにより、物体が音響放射圧のために平均圧力が低い節に引き寄せられますが、物体が波長の10%以下の場合は、ポンデロモーティブ力が音響放射に打ち勝ち、物体は腹に移動します。[3]対象物体は、使用する音の波長よりもかなり小さくなければなりません。ただし、特定の状況が軽視されている場合は除きます。節間の距離と対象の物体に使用する波長を調整することで、使用波長よりもはるかに大きな物体を浮上させることができますが、これには慎重な計算と多くの試行錯誤が必要です。 1次元定在波を用いた微粒子操作は、1982年の研究論文「繊維懸濁液の超音波検査」で初めて報告されました。[4]

音波は生物学的対象に対して安全であることが証明されており、生物医学的用途に最適です。 [5]最近では、フローサイトメトリー、細胞分離、細胞捕捉、単一細胞操作、ナノ材料操作など、ミリメートル未満の物体の操作に音響ピンセットが応用されています。 [6]

方法

音響ピンセットを使用して、これら 10 個の細胞を数百ミクロンのスケールで文字「A」と「T」のパターンに配置しました

定在 音場において、物体は音響放射力を受け、音場の特定の領域に移動します。[2]物体の密度や圧縮率などの特性に応じて、音圧の節(最小圧力領域)または圧力の腹(最大圧力領域)に移動するように誘導できます。[3]その結果、これらの節の位置を制御することで、音波を用いた物体の正確な移動が可能になります。音響ピンセットは、高価な機器や複雑な実験装置を必要としません。[要出典]

基礎理論

音場中の粒子は、音波、流体、粒子間の相互作用によって生じる力によって動かされます。これらの力(音響放射力、粒子間の二次場力、ストークス抗力など)は、音響ピンセット技術の基礎となる 音響泳動現象を生み出します。

音響放射力

粒子が音波場中に浮遊している場合、音波の散乱によって生じる音響放射力が粒子に作用します。この現象は、1934年にルイス・キングによって理想流体中の非圧縮性粒子を対象に初めてモデル化され、解析されました。[7]吉岡と川島は1955年に平面波における圧縮性粒子に作用する音響放射力を計算しました。[8]ゴルコフは先行研究を要約し、粒子の大きさが音波の波長よりもはるかに小さい場合、任意の音響場において粒子に作用する平均的な力を決定するための式を提案しました。[2]最近、ブルースはこの問題を再検討し、音響放射力の詳細な導出を行いました。[9]

粒子に働く力の線画
 1:微粒子に対する音響放射力

図1に示すように、微小粒子に対する音響放射力は、粒子周囲の近接場領域における不均一な運動量フラックス から生じます。この運動量フラックスは、入射する音波と、音波が粒子表面を伝播する際の散乱によって生じます。理想流体中の音波の波長よりもはるかに小さい直径を持つ圧縮性球状粒子の場合、音響放射力は で計算できます。ここでは音響ポテンシャルエネルギーとも呼ばれる所定の量です。[2] [9]音響ポテンシャルエネルギーは次のように表されます。 F r a d = U {\displaystyle \mathbf {\mathit {F^{rad}}} =-\nabla {\mathit {U}}} F r a d = U {\displaystyle \mathbf {\mathit {F^{rad}}} =-\nabla {\mathit {U}}} U {\displaystyle {\mathit {U}}}

U = V 0 ( p i n 2 ¯ 2 ρ f c f 2 f 1 3 ρ f v i n 2 ¯ 4 f 2 ) {\displaystyle U={V_{0}}({{\overline {p_{in}^{2}}} \over {2{\rho _{f}}c_{f}^{2}}}{f_{1}}-{{3{\rho _{f}}{\overline {v_{in}^{2}}}} \over 4}{f_{2}})}

どこ

  • V 0 {\displaystyle {\mathit {V_{\mathit {0}}}}} 粒子の体積、
  • p i n {\displaystyle {\mathit {p_{\mathit {in}}}}} 音圧は
  • v i n {\displaystyle {\mathit {v_{\mathit {in}}}}} 音響粒子の速度であり、
  • ρ f {\displaystyle {\mathit {\rho }}_{\mathit {f}}} は流体の質量密度であり、
  • c f {\displaystyle {\mathit {c_{\mathit {f}}}}} 流体の音速である。
  • < > ¯ {\displaystyle {\overline {<*>}}} は時間平均項であり、 1 T 0 T ( ) d t {\displaystyle {1 \over T}\int \limits _{0}^{T}{(*)}dt}

係数とは次の式で計算できる。 f 1 {\displaystyle {\mathit {f_{\mathit {1}}}}} f 2 {\displaystyle {\mathit {f_{\mathit {2}}}}} f 1 = 1 ρ f c f 2 ρ p c p 2 {\displaystyle {f_{1}}=1-{{{\rho _{f}}c_{f}^{2}} \over {{\rho _{p}}c_{p}^{2}}}} f 2 = 2 ( ρ p ρ f ) 2 ρ p + ρ f {\displaystyle {f_{2}}={{2({\rho _{p}}-{\rho _{f}})} \over {2{\rho _{p}}+{\rho _{f}}}}}

どこ

  • ρ p {\displaystyle {\mathit {\rho }}_{\mathit {p}}} は粒子の質量密度であり、
  • c p {\displaystyle {\mathit {c}}_{\mathit {p}}} 粒子の音速です。

定在波における音響放射力

定在波は安定した音響ポテンシャルエネルギー場を形成するため、多くの音響ピンセット用途において望ましい、安定した音響放射力分布を作り出すことができます。1次元平面定在波の場合、音響場は次のように表されます。[9]

A i n ( x , t ) = P 0 ρ f ω sin ( k x ) sin ( ω t ) {\displaystyle {A_{in}}(x,t)={{-P_{0}} \over {{\rho _{f}}\omega }}\sin(kx)\sin(\omega t)}

p i n ( x , t ) = P 0 cos ( k x ) sin ( ω t ) {\displaystyle {p_{in}}(x,t)={P_{0}}\cos(kx)\sin(\omega t)}

v in ( x , t ) = P 0 ρ f c f sin ( k x ) cos ( ω t ) ) e x {\displaystyle {\textbf {v}}_{\textbf {in}}(x,t)={\frac {-P_{0}}{\rho _{f}c_{f}}}\sin(kx)\cos(\omega t)){\textbf {e}}_{\textbf {x}}}

どこ

  • A i n {\displaystyle {\mathit {A_{\mathit {in}}}}} 音響粒子の変位であり、
  • P 0 {\displaystyle {\mathit {P_{\mathit {0}}}}} 音圧振幅、
  • ω {\displaystyle {\mathit {\omega }}} は角速度であり
  • k {\displaystyle {\mathit {k}}} 波数です。

これらのフィールドを用いて、時間平均項を得ることができます。これらは次のとおりです。

p i n 2 ¯ = 1 2 P 0 2 cos 2 ( k x ) {\displaystyle {\overline {p_{in}^{2}}}={1 \over 2}P_{0}^{2}{\cos ^{2}}(kx)}

v i n 2 ¯ = P 0 2 2 ρ f 2 c f 2 sin 2 ( k x ) {\displaystyle {\overline {v_{in}^{2}}}={{P_{0}^{2}} \over {2\rho _{f}^{2}c_{f}^{2}}}{\sin ^{2}}(kx)}

したがって、音響の位置エネルギーは次のようになります。

U = V 0 P 0 2 4 ρ f c f 2 [ cos 2 ( k x ) f 1 3 2 sin 2 ( k x ) f 2 ] {\displaystyle U={V_{0}}{{P_{0}^{2}} \over {4{\rho _{f}}c_{f}^{2}}}[{\cos ^{2}}(kx){f_{1}}-{3 \over 2}{\sin ^{2}}(kx){f_{2}}]}

次に、音響放射力を微分によって求めます。

F x r a d = x U = V 0 k E a c sin ( 2 k x ) Φ {\displaystyle {F_{x}}^{rad}=-{\partial _{x}}U={V_{0}}k{E_{ac}}\sin(2kx)\Phi }

E a c = P 0 2 4 ρ f c f 2 {\displaystyle {E_{ac}}={{P_{0}^{2}} \over {4{\rho _{f}}c_{f}^{2}}}} 、、 Φ = f 1 + 3 2 f 2 = 5 ρ p 2 ρ f 2 ρ p + ρ f ρ f c f 2 ρ p c p 2 {\displaystyle \Phi ={f_{1}}+{3 \over 2}{f_{2}}={{5{\rho _{p}}-2{\rho _{f}}} \over {2{\rho _{p}}+{\rho _{f}}}}-{{{\rho _{f}}c_{f}^{2}} \over {{\rho _{p}}c_{p}^{2}}}}

音響放射力の線図
図2: 音圧波形に沿った圧力節(PN)と腹(AN)の位置

どこ

項は、放射力の周期が圧力周期の半分であることを示しています。また、コントラスト係数は、粒子と流体の特性に応じて正または負の値をとります。 の値が正の場合、図2に示すように、放射力は圧力腹から圧力節へと向かい、粒子は圧力節へと押し出されます。 sin ( 2 k x ) {\displaystyle \sin({\mathit {2kx}})} Φ {\displaystyle \Phi }

二次音響力

懸濁液中の複数の粒子が定在波場にさらされると、音響放射力だけでなく、他の粒子による散乱波によって生じる二次的な音響力も作用します。粒子間力はビャクネス力と呼ばれることもあります。同一粒子間の粒子間力を表す簡略化された式は以下のとおりです。[10] [11]

F B ( x ) = 4 π a 6 [ ( ρ p ρ f ) 2 ( 3 cos 2 θ 1 ) 6 ρ f d 4 v i n 2 ( x ) ω 2 ρ f 9 d 2 ( 1 ρ p c p 2 1 ρ f c f 2 ) 2 p i n 2 ( x ) ] {\displaystyle {F_{B}}(x)=4\pi {a^{6}}[{{{{({\rho _{p}}-{\rho _{f}})}^{2}}(3{{\cos }^{2}}\theta -1)} \over {6{\rho _{f}}{d^{4}}}}{v_{in}}^{2}(x)-{{\omega ^{2}{\rho _{f}}} \over {9{d^{2}}}}{({1 \over {{\rho _{p}}c_{p}^{2}}}-{1 \over {{\rho _{f}}c_{f}^{2}}})^{2}}p_{in}^{2}(x)]}

どこ

  • a {\displaystyle {\mathit {a}}} 粒子の半径、
  • d {\displaystyle {\mathit {d}}} 粒子間の距離であり、
  • θ {\displaystyle {\mathit {\theta }}} 粒子の中心線と入射音波の伝播方向との間の角度です。

力の符号は力の方向を表します。負の符号は引力、正の符号は斥力です。方程式の左側は音響粒子の速度振幅に依存し、右側は音響圧力振幅に依存します。速度依存項は、粒子が波の伝播と揃っている場合 (Θ=0°) 斥力となり、波の伝播と垂直な場合 (Θ=90°) は負になります。圧力依存項は粒子の向きに影響されず、常に引力となります。正のコントラスト係数の場合、気泡や脂質小胞の場合のように、粒子が速度の節 (圧力の腹) に追いやられるにつれて、速度依存項は減少します。同様に、圧力依存項は、水溶液中のほとんどの固体粒子のように、粒子が圧力の節 (速度の腹) に追いやられるにつれて減少します。 v i n ( x ) {\displaystyle {\mathit {v}}_{\mathit {in}}({\mathit {x}})} p i n ( x ) {\displaystyle {\mathit {p}}_{\mathit {in}}({\mathit {x}})}

散乱に関連する二次音響力に加えて、各粒子の音響境界層によって生成されるさまざまな音響ストリーミング場の相互作用から生じる流れ場(マイクロストリーミングと呼ばれることもある)は、各粒子表面に追加の粘性せん断力を誘発する可能性があります。その結果、完全に粘性のある定式化において二次音響力に追加の寄与が生じます。[12]二次音響力に対する粘性効果は、上で例示した完全流体定式化と比較すると重要になり、特定の限界ケースでは支配的になることさえあり、非粘性理論によって予測されるものよりも量的にも質的にも異なる結果をもたらします。[13]粘性寄与の関連性は、調査されている特定のケースによって大きく異なるため、特定のシナリオに対して適切な二次音響力モデルを選択する際には細心の注意を払う必要があります。

二次力の影響は通常非常に弱く、粒子間の距離が非常に小さい場合にのみ効果を発揮します。凝集および沈降の応用において、粒子は音響放射力によって最初にノードに集まります。粒子間の距離が小さくなるにつれて、二次力はさらなる凝集を促し、クラスターが沈降を開始するのに十分な重さになるまで続きます。

音響ストリーミング

音響流は、音響場における振動の非線形成分によって生成される定常流である。[詳細な説明が必要]音響流は、そのメカニズムに応じて、エッカート流とレイリー流の2つの一般的なタイプに分類できる。[14] [15]エッカート流は、高振幅の音波が流体中を伝播・減衰する際に生じる時間平均運動量流束によって駆動される。レイリー流は「境界駆動流」とも呼ばれ、粘性境界層におけるレイノルズ応力によって強制される。どちらの駆動メカニズムも、時間平均非線形効果に由来する。

非線形音響流現象を解析するために摂動法が用いられる。[16]この問題の支配方程式は質量保存則とナビエ・ストークス方程式である。

t ρ = ( ρ v ) {\displaystyle \partial _{\mathit {t}}{\mathit {\rho }}=-\nabla \cdot ({\mathit {\rho }}{\textbf {v}})}

ρ [ t v + ( v ) v ] = p + μ 2 v + β μ ( v ) ) {\displaystyle {\mathit {\rho }}[\partial _{\mathit {t}}{\textbf {v}}+({\textbf {v}}\cdot \nabla ){\textbf {v}}]=-\nabla {\mathit {p}}+{\mathit {\mu }}\nabla ^{2}{\textbf {v}}+{\mathit {\beta }}{\mathit {\mu }}\nabla (\nabla \cdot {\textbf {v}}))}

どこ

  • ρ {\displaystyle {\mathit {\rho }}} 流体の密度である。
  • v {\displaystyle {\textbf {v}}} 流体粒子の速度であり、
  • p {\displaystyle {\mathit {p}}} プレッシャーは、
  • μ {\displaystyle {\mathit {\mu }}} 流体の動粘性であり、
  • β {\displaystyle {\mathit {\beta }}} 粘度比です。

摂動級数は、、、と表すことができます。これらは高次の項が低次の項よりもはるかに小さい、減少級数です。 p = p 0 + p 1 + p 2 {\displaystyle p={p_{0}}+{p_{1}}+{p_{2}}} v = 0 + v 1 + v 2 {\displaystyle {\textbf {v}}={\textbf {0}}+{\textbf {v}}_{\textbf {1}}+{\textbf {v}}_{\textbf {2}}} ρ = ρ 0 + ρ 1 + ρ 2 {\displaystyle {\mathit {\rho }}={\mathit {\rho }}_{0}+{\mathit {\rho }}_{1}+{\mathit {\rho }}_{2}}

液体は零次状態では静止し均質である。摂動級数を質量保存則とナビエ・ストークス方程式に代入し、の関係を用いると、一次方程式は一次の項をまとめることで得られる。 p 1 = c f 2 ρ 1 {\displaystyle {p_{1}}=c_{f}^{2}{\rho _{1}}}


  
    
      
        
          
          
            
              t
            
          
        
        
          
            
              p
            
          
          
            1
          
        
        =
        
        
          
            
              ρ
            
          
          
            0
          
        
        
          
            
              c
            
          
          
            
              f
            
          
          
            2
          
        
        
        
        
          
            
              v
            
          
          
            
              1
            
          
        
      
    
    {\displaystyle \partial _{\mathit {t}}{\mathit {p}}_{1}=-{\mathit {\rho }}_{0}{\mathit {c}}_{\mathit {f}}^{2}\nabla \cdot {\textbf {v}}_{\textbf {1}}}
  

  
    
      
        
          
            
              ρ
            
          
          
            0
          
        
        
          
          
            
              t
            
          
        
        
          
            
              v
            
          
          
            
              1
            
          
        
        =
        
        
        
          
            
              p
            
          
          
            1
          
        
        +
        
          
            μ
          
        
        
          
          
            2
          
        
        
          
            
              v
            
          
          
            
              1
            
          
        
        +
        
          
            β
          
        
        
          
            μ
          
        
        
        (
        
        
        
          
            
              v
            
          
          
            
              1
            
          
        
        )
      
    
    {\displaystyle {\mathit {\rho }}_{0}\partial _{\mathit {t}}{\textbf {v}}_{\textbf {1}}=-\nabla {\mathit {p}}_{1}+{\mathit {\mu }}\nabla ^{2}{\textbf {v}}_{\textbf {1}}+{\mathit {\beta }}{\mathit {\mu }}\nabla (\nabla \cdot {\textbf {v}}_{\textbf {1}})}
  

同様に、2次方程式も求められます。


  
    
      
        
          
          
            
              t
            
          
        
        
          
            
              ρ
            
          
          
            2
          
        
        =
        
        
          
            
              ρ
            
          
          
            0
          
        
        
        
        
          
            
              v
            
          
          
            
              2
            
          
        
        
        
        
        (
        
          
            
              ρ
            
          
          
            1
          
        
        
          
            
              v
            
          
          
            
              1
            
          
        
        )
      
    
    {\displaystyle \partial _{\mathit {t}}{\mathit {\rho }}_{2}=-{\mathit {\rho }}_{0}\nabla \cdot {\textbf {v}}_{\textbf {2}}-\nabla \cdot ({\mathit {\rho }}_{1}{\textbf {v}}_{\textbf {1}})}
  

  
    
      
        
          
            
              ρ
            
          
          
            0
          
        
        
          
          
            
              t
            
          
        
        
          
            
              v
            
          
          
            
              2
            
          
        
        =
        
        
        
          
            
              p
            
          
          
            2
          
        
        +
        
          
            μ
          
        
        
          
          
            2
          
        
        
          
            
              v
            
          
          
            
              2
            
          
        
        +
        
          
            β
          
        
        
          
            μ
          
        
        
        (
        
        
        
          
            
              v
            
          
          
            
              2
            
          
        
        )
        
        
          
            
              ρ
            
          
          
            1
          
        
        
          
          
            t
          
        
        
          
            
              v
            
          
          
            
              1
            
          
        
        
        
          
            
              ρ
            
          
          
            0
          
        
        (
        
          
            
              v
            
          
          
            
              1
            
          
        
        
        
        )
        
          
            
              v
            
          
          
            
              1
            
          
        
      
    
    {\displaystyle {\mathit {\rho }}_{0}\partial _{\mathit {t}}{\textbf {v}}_{\textbf {2}}=-\nabla {\mathit {p}}_{2}+{\mathit {\mu }}\nabla ^{2}{\textbf {v}}_{\textbf {2}}+{\mathit {\beta }}{\mathit {\mu }}\nabla (\nabla \cdot {\textbf {v}}_{\textbf {2}})-{\mathit {\rho }}_{1}\partial _{t}{\textbf {v}}_{\textbf {1}}-{\mathit {\rho }}_{0}({\textbf {v}}_{\textbf {1}}\cdot \nabla ){\textbf {v}}_{\textbf {1}}}
  

一次方程式については、ナビエ・ストークス方程式の時間微分を取り、質量保存則を挿入すると、次の複合方程式が得られます。

1 c f 2 t 2 p 1 = [ 1 ( 1 β ) μ ρ 0 c f 2 t ] 2 p 1 {\displaystyle {1 \over {c_{f}^{2}}}\partial _{t}^{2}{p_{1}}=[1-{{(1-\beta )\mu } \over {{\rho _{0}}c_{f}^{2}}}{\partial _{t}}]{\nabla ^{2}}{p_{1}}}

これは粘性減衰を伴う音波方程式です。物理的には、圧と音響粒子の速度として解釈できます。 p 1 {\displaystyle {\mathit {p}}_{1}} v 1 {\displaystyle {\textbf {v}}_{\textbf {1}}}

2次方程式は、質量源と力源を持つ流体の運動を記述する支配方程式と考えることができる。一般的に、音響流は定常平均流であり、応答時間スケールは音響振動のそれよりもはるかに小さい。音響流を表現するために、通常、時間平均項が用いられる。 を用いることで、時間平均2次方程式が得られる。 [ ( ρ 1 v 1 ) ) ] {\displaystyle [-\nabla \cdot ({\mathit {\rho }}_{1}{\textbf {v}}_{\textbf {1}}))]} [ ρ 1 t v 1 ρ 0 ( v 1 ) v 1 ] {\displaystyle [-\rho _{1}\partial _{t}{\textbf {v}}_{\textbf {1}}-\rho _{0}({\textbf {v}}_{\textbf {1}}\cdot \nabla ){\textbf {v}}_{\textbf {1}}]} v 2 ¯ {\displaystyle {\overline {{\textbf {v}}_{\textbf {2}}}}} t ρ 2 ¯ = 0 {\displaystyle {\overline {\partial _{\mathit {t}}{\mathit {\rho }}_{2}}}=0}


  
    
      
        
          
            
              ρ
            
          
          
            0
          
        
        
        
        
          
            
              
                
                  v
                
              
              
                
                  2
                
              
            
            ¯
          
        
        =
        
        
        
        
          
            
              (
              
                
                  
                    ρ
                  
                
                
                  1
                
              
              
                
                  
                    v
                  
                
                
                  
                    1
                  
                
              
              )
            
            ¯
          
        
      
    
    {\displaystyle {\mathit {\rho }}_{0}\nabla \cdot {\overline {{\textbf {v}}_{\textbf {2}}}}=-\nabla \cdot {\overline {({\mathit {\rho }}_{1}{\textbf {v}}_{\textbf {1}})}}}
  

  
    
      
        
          
            μ
          
        
        
          
          
            2
          
        
        
          
            
              
                
                  v
                
              
              
                
                  2
                
              
            
            ¯
          
        
        +
        
          
            β
          
        
        
          
            μ
          
        
        
        (
        
        
        
          
            
              
                
                  v
                
              
              
                
                  2
                
              
            
            ¯
          
        
        )
        
        
          
            
              
              
                
                  
                    p
                  
                
                
                  2
                
              
            
            ¯
          
        
        =
        
          
            
              
                
                  
                    ρ
                  
                
                
                  1
                
              
              
                
                
                  
                    t
                  
                
              
              
                
                  
                    v
                  
                
                
                  
                    1
                  
                
              
            
            ¯
          
        
        
        
          
            
              ρ
            
          
          
            0
          
        
        
          
            
              (
              
                
                  
                    v
                  
                
                
                  
                    1
                  
                
              
              
              
              )
              
                
                  
                    v
                  
                
                
                  
                    1
                  
                
              
            
            ¯
          
        
      
    
    {\displaystyle {\mathit {\mu }}\nabla ^{2}{\overline {{\textbf {v}}_{\textbf {2}}}}+{\mathit {\beta }}{\mathit {\mu }}\nabla (\nabla \cdot {\overline {{\textbf {v}}_{\textbf {2}}}})-{\overline {\nabla {\mathit {p}}_{2}}}={\overline {{\mathit {\rho }}_{1}\partial _{\mathit {t}}{\textbf {v}}_{\textbf {1}}}}-{\mathit {\rho }}_{0}{\overline {({\textbf {v}}_{\textbf {1}}\cdot \nabla ){\textbf {v}}_{\textbf {1}}}}}
  
青い波打つ模様に囲まれた白い円
図3: 円筒柱の周りの音響流の断面

純粋な1次項は定義上調和波であるため、時間平均すると打ち消されることに注意することが重要です。つまり、純粋な1次項は純粋な正弦波であり、平均が0であるため、それらを含む項はすべて打ち消されます。ただし、2次項は調和ではないため、時間平均では打ち消されません。[17]これは、音響流を理解する上で最も重要です。単純な振動運動に関連する1次項は、2次項よりもはるかに大きな値を持ち、したがって振動の時間スケールで支配的です。ただし、これらの1次項は純粋な正弦波であり、準定常状態であるため、各振動サイクルの後に繰り返され、正味の流体の流れは生成されません。一方、2次項は調和ではないため、累積効果を持つ可能性があり、小さいにもかかわらず、多くの振動サイクルにわたって追加され、音響流として識別される正味の定常流れの発生につながります。

したがって、音響流を決定する上で、2次方程式が最も重要です。ナビエ・ストークス方程式は単純なケースでのみ解析的に解くことができるため、通常は数値解析手法が用いられ、その中でも有限要素法(FEM)が最も一般的な手法です。FEMは音響流現象のシミュレーションに利用できます。図3は、FEMによって計算された円柱周りの音響流の一例です。

前述のように、音響流は音響減衰に起因する質量源と力源によって駆動されます。しかし、これらが音響流の駆動力の唯一の要因ではありません。境界振動も、特に「境界駆動流」に寄与することがあります。このような場合、境界条件も摂動法を用いて処理し、2次方程式に適切に適用する必要があります。

粒子の運動

音響場において、重力と浮力とが釣り合っている懸濁粒子の運動は、音響放射力ストークス抗力という2つの力によって決まります。ニュートンの法則を適用すると、この運動は次のように記述できます。


  
    
      
        
          
            m
          
        
        
          
            
              d
              
                
                  u
                
              
            
            
              d
              t
            
          
        
        =
        
          
            
              F
            
          
          
            
              rad
            
          
        
        +
        
          
            
              F
            
          
          
            
              drag
            
          
        
      
    
    {\displaystyle {\mathit {m}}{\frac {d{\textit {u}}}{\mathit {dt}}}={\textit {F}}^{\textit {rad}}+{\textit {F}}^{\textit {drag}}}
  

  
    
      
        
          
            
              F
            
          
          
            
              drag
            
          
        
        =
        6
        π
        
          
            a
          
        
        
          
            μ
          
        
        (
        
          
            v
          
        
        
        
          
            u
          
        
        )
      
    
    {\displaystyle {\textit {F}}^{\textit {drag}}=6\pi {\mathit {a}}{\mathit {\mu }}({\textit {v}}-{\textit {u}})}
  

どこ

  • v {\displaystyle {\textit {v}}} 流体の速度、
  • u {\displaystyle {\textit {u}}} 粒子の速度です。

静的流れのアプリケーションでは、流体速度は音響流によって生じます。音響流の強度は、入力の電力と周波数、および流体媒体の特性に依存します。一般的な音響ベースのマイクロデバイスでは、動作周波数はkHzからMHzの範囲です。振動振幅は0.1 nmから1 μmの範囲です。流体として水を使用すると、音響流の強度は1 μm/sから1 mm/sの範囲と推定されます。したがって、ほとんどの連続流アプリケーションでは、音響流は主流よりも小さくなります。これらのアプリケーションでは、抗力は主に主流によって発生します。

アプリケーション

細胞分離

密度や圧縮強度の異なる細胞は、理論的には音響力で分離できます。音響ピンセットを用いて脂質粒子を赤血球から分離できる可能性が示唆されています[18]これは人工心肺装置を用いた心臓手術における問題であり、現在の技術では解決が困難です。提案によれば、チャネルを通過する血漿に音響力を加えると、赤血球は中央の圧力節に集まり、脂質粒子は両側の腹に集まります(図4参照)。チャネルの端では、分離された細胞と粒子は別々の出口から排出されます。

音響法は、異なるサイズの粒子を分離するためにも使用できる可能性があります。一次音響放射力の式によれば、大きな粒子は小さな粒子よりも大きな力を受けます。Shiらは、インターデジタルトランスデューサ(IDT)を用いて、マイクロ流体チャネルの中央に圧力ノードを有する定在表面音響波(SSAW)場を発生させ、異なる直径の微粒子を分離したことを報告しました。[19]チャネルの端から異なるサイズの粒子の混合物を導入すると、大きな粒子は中央に向かってより速く移動し、中央の出口に集められます。小さな粒子は、側面の出口から集められる前に中央の出口に移動することはできません。この実験装置は、血液成分、細菌、ハイドロゲル粒子の分離にも使用されています。[20] [21] [22]

3D細胞フォーカス

蛍光活性化セルソーター(FACS)は、細胞を含む流体の流れを集束させ、個々の細胞からの蛍光を検出し、目的の細胞を他の細胞から分離することで細胞を選別します。スループットは高いものの、購入と維持に費用がかかり、構成が複雑で大型です。また、高いせん断圧、衝撃力、電磁力によって細胞生理機能に影響を与え、細胞や遺伝子に損傷を与える可能性があります。音響力は細胞にとって危険ではなく[要出典]、音響ピンセットと光学/電気モジュールを統合することで、より小型で安価な装置に細胞分析と選別を同時に行う技術開発が進められています。

マイクロ流体デバイスにおける細胞/粒子の3次元集束を実現するために、音響ピンセットが開発されている[23] 。圧電基板上に一対のインターデジタルトランスデューサ(IDT)が堆積され、マイクロ流体チャネルが基板に接合され、2つのIDTの間に配置される。微粒子溶液は圧力駆動流によってマイクロ流体チャネルに注入される。両方のIDTにRF信号が印加されると、2系列の表面弾性波(SAW)がマイクロチャネル内の微粒子懸濁液に向かって反対方向に伝播する。2つのSAWの建設的干渉により、SSAWが形成される。チャネル内に縦モードの漏れ波が発生し、粒子に対して横方向に作用する圧力変動を引き起こす。その結果、チャネル内の懸濁粒子は、粒子と媒体の密度と圧縮率に応じて、圧力節点または腹点に向かって押し出される。チャネル幅が1つの圧力節点(または腹点)のみを覆う場合、粒子はその節点に集束される。

細胞/粒子は水平方向だけでなく、垂直方向にも集束させることができます。[24] SSAWを作動させると、ランダムに分布した粒子は垂直方向に一列に並ぶ流れ(図10c)に集束されます。定在表面音響波(SSAW)をベースとした3D粒子/細胞集束が可能なマイクロデバイスとレーザー誘起蛍光(LIF)検出システムを統合することで、音響ピンセットは高スループット単一細胞解析のためのマイクロフローサイトメーターとして開発されます。

チャープ型インターデジタルトランスデューサ[ 25] [26]の持つ調整可能性により、細胞を1ステップで複数(例えば5つ)の出口チャネルに正確に分別することが可能になります。これは、通常2つの出口チャネルにしか細胞を分別しない既存の多くの分別方法に比べて大きな利点です。

非侵襲的な細胞捕捉とパターン形成

エッチングされた流体チャネルを備えたガラス反射板が、トランスデューサーを保持するPCBに固定されています。チップに注入された細胞は、チャネル内に形成される超音波定在波に捕捉されます。音響力によって細胞はチャネル中央にクラスター状に集束し、挿入図に示されています。捕捉はトランスデューサー表面近くで行われるため、実際の捕捉位置は3D画像に示すように近傍場圧力分布によって決まります。細胞は局所的な圧力最小値付近にクラスター状に捕捉され、捕捉された細胞数に応じて異なるパターンを形成します。グラフのピークは圧力最小値に対応しています。

単一細胞、粒子、または生物の操作

ブロック体を形成する点線
図4:単一細胞操作の痕跡

単一細胞の操作は、細胞微小環境の制御や特定の対象細胞の単離など、多くの生物学研究にとって重要です。音響ピンセットは、マイクロメートルレベルの解像度で個々の細胞を操作できることが実証されています。細胞の直径は通常10~20μmです。単一細胞を操作するための解像度要件を満たすには、短波長の音波を使用する必要があります。この場合、より短い波長の音波(通常200μm未満)を使用できる表面弾性波(SAW)の方がバルク弾性波(BAW)よりも好まれます。[27] Dingらは、規定の経路で単一細胞を操作できるSSAWマイクロデバイスを報告しました。[28]図6は、単一細胞の動きが音響ピンセットで細かく制御できることを実証した記​​録です。デバイスの動作原理は、SSAWフィールド内の圧力ノードの制御された動きにあります。Dingらは、 SSAWフィールド内の圧力ノードの制御された動きを制御入力周波数を変化させることで圧力ノードの位置を調整できるSSAWを生成できるチャープ型インターデジタルトランスデューサ(IDT)を採用した。また、ミリメートルサイズの微生物であるC. eleganを同様の方法で操作できることを示した。さらに、音響処理後の細胞代謝と増殖を調べたところ、対照群と比較して有意差は見られず、音響ベースの操作が非侵襲的であることが示された。チャープ型IDTの使用に加えて、位相シフトに基づく単一粒子/細胞操作も報告されている。[29] [30] [31]

単一生体分子の操作

Sittersらは、音響を用いてDNAやタンパク質などの単一生体分子を操作できることを示した。発明者らが音響力分光法と呼ぶこの方法は、単一分子の力の応答を測定することを可能にする。これは、分子の片側に小さな微小球を付着させ、もう片側にそれを表面に付着させることで実現される。定在音波によって微小球を表面から押し出すことで、分子は効果的に引き伸ばされる。[32]

有機ナノ材料の操作

ポリマー分散液晶(PDLC)ディスプレイは、音響ピンセットを用いて不透明から透明に切り替えることができます。硬化したPDLCフィルムと一対のインターデジタルトランスデューサ(IDT)を圧電基板上に集積することで、SAW駆動のPDLC光シャッターが実証されています。[33]

無機ナノ材料の操作

音響ピンセットは、調整可能なナノワイヤパターニングのための簡便なアプローチを提供します。このアプローチでは、インターデジタルトランスデューサによってSSAWが生成され、圧電基板上に周期的な交流(AC)電界が誘起され、その結果、懸濁液中の金属ナノワイヤがパターン化されます。パターンは、液体が蒸発した後に基板上に堆積されます。SSAW電界の分布を制御することで、金属ナノワイヤは平行配列や垂直配列など、様々なパターンに組み立てられます。ナノワイヤアレイの間隔は、表面弾性波の周波数を制御することで調整できます。[34]

選択的操作

ほとんどの音響ピンセットは多数の物体をまとめて操作できますが[27]、補完的な機能として、隣接する物体を動かさずにクラスター内の単一の粒子を操作できるというものがあります。この目的を達成するには、音響トラップを空間的に局所化する必要があります。最初のアプローチは、高度に集束された音響ビームを用いることです[35] 。多くの対象粒子は音響場の節点に引き寄せられ、焦点から押し出されるため、この種の粒子をトラップするには、強い集束と焦点における圧力振幅の最小化(トラップを形成する強度のリングで囲まれた部分)を両立する特定の波構造が必要です。これらの特定の条件は、位相秩序がゼロより大きいベッセルビーム、いわゆる「音響渦」によって満たされます。この種の波構造を用いることで、プログラム可能な電子機器によって駆動されるトランスデューサーアレイを用いて、2次元[36]および3次元[37] [38]における粒子の選択的操作が実証されています。

あるいは、音響エネルギーを局在化するための別のアプローチは、局所的な音響定在波を生成するためにナノ秒スケールのパルス場を使用することに依存しています。[39]

高周波ピンセットとホログラフィックインターデジタルトランスデューサ(IDT)

微小物体を個別に選択的に操作するには、必要な空間分解能(通常、選択的に操作される物体のサイズと波長が同等でなければならない)に達するのに十分な高周波数で、音響渦(前節参照)などの複雑な音響場を合成する必要がある。複雑な波動場を合成するためのホログラフィック手法は数多く開発されており、トランスデューサーアレイ[40] [41] [36] [42 ] [38] 3Dプリントホログラム[43]メタマテリアル[44]回折格子[45] [46]な​​どがある。しかしながら、これらの手法はすべて比較的低い周波数に限定されており、マイクロメートル単位の粒子、細胞、微生物を個別に処理するには分解能が不十分である。一方、インターディジタルトランスデューサー(IDT)は、GHz周波数までの音響波動場を合成する信頼性の高い技術として知られていた。[47]

参照

参考文献

  1. ^ Ozcelik, A., Rufo, J., Guo, F. 他. 生命科学のための音響ピンセット. Nat Methods 15, 1021–1028 (2018). https://doi.org/10.1038/s41592-018-0222-9
  2. ^ abcd ゴルコフ、LP;ソビエトの物理学 - Doklady、1962、6(9)、773-775。
  3. ^ ab Nilsson, Andreas; Petersson, Filip; Jönsson, Henrik; Laurell, Thomas (2004). 「マイクロ流体チップにおける懸濁粒子の音響制御」. Lab Chip . 4 (2). 英国王立化学協会 (RSC): 131– 135. doi :10.1039/b313493h. ISSN  1473-0197. PMID  15052353.
  4. ^ ディオン、JL;マルッタ、A.シエロ、P. (1982)。 「ファイバーサスペンションの超音波検査」。アメリカ音響学会誌72 (5): 1524–1526書誌コード:1982ASAJ...72.1524D。土井:10.1121/1.388688。
  5. ^ Lin, Sz-Chin Steven; Mao, Xiaole; Huang, Tony Jun (2012). 「表面音響波(SAW)音響泳動:現在と未来」. Lab on a Chip . 12 (16). Royal Society of Chemistry (RSC): 2766– 2770. doi :10.1039/c2lc90076a. ISSN  1473-0197. PMC 3992433. PMID 22781941  . 
  6. ^ Ding, Xiaoyun; Li, Peng; Lin, Sz-Chin Steven; Stratton, Zackary S.; Nama, Nitesh; Guo, Feng; Slotcavage, Daniel; Mao, Xiaole; Shi, Jinjie; Costanzo, Francesco; Huang, Tony Jun (2013). 「表面音響波マイクロ流体工学」. Lab on a Chip . 13 (18). Royal Society of Chemistry (RSC): 3626– 3649. doi :10.1039/c3lc50361e. ISSN  1473-0197. PMC 3992948. PMID 23900527  . 
  7. ^ キング, ルイス・V. (1934年11月15日). 「球面上の音響放射圧について」.ロンドン王立協会紀要. シリーズA, 数学・物理科学. 147 (861). 王立協会: 212– 240. Bibcode :1934RSPSA.147..212K. doi : 10.1098/rspa.1934.0215 . ISSN  2053-9169.
  8. ^ 吉岡和也、川島裕子;アクスティカ、1955、5(3)​​、167-173。
  9. ^ abc Bruus, Henrik (2012). 「音響流体工学 7:微小粒子に対する音響放射力」. Lab on a Chip . 12 (6). 英国王立化学協会 (RSC): 1014– 1021. doi :10.1039/c2lc21068a. ISSN  1473-0197. PMID  22349937.
  10. ^ ワイザー、MAH; RE、Apfel、EA の Neppiras。アクスティカ、1984、56(2)、114-119。
  11. ^ ローレル, トーマス; ピーターソン, フィリップ; ニルソン, アンドレアス (2007). 「細胞と粒子の音響分離および操作のためのチップ統合戦略」.化学協会レビュー. 36 (3). 英国王立化学協会 (RSC): 492– 506. doi :10.1039/b601326k. ISSN  0306-0012. PMID  17325788.
  12. ^ Doinikov, Alexander A. (1994). 「粘性流体中の圧縮性球にかかる音響放射圧」 . Journal of Fluid Mechanics . 267 : 1– 22. Bibcode :1994JFM...267....1D. doi :10.1017/S0022112094001096. S2CID  123245307. 2023年4月10日閲覧
  13. ^ Baasch, Thierry; Pavlic, Alen; Dual, Jürg (2019年12月). 「定在波中の重い粒子に作用する音響放射力は、音響マイクロストリーミングによって支配される可能性がある」 . Phys. Rev. E. 100 ( 6) 061102. Bibcode :2019PhRvE.100f1102B. doi :10.1103/PhysRevE.100.061102. PMID  31962519. S2CID  210870075. 2023年4月10日閲覧
  14. ^ ライトヒル卿 (1978). 「音響ストリーミング」. Journal of Sound and Vibration . 61 (3). Elsevier BV: 391– 418. Bibcode :1978JSV....61..391L. doi :10.1016/0022-460x(78)90388-7. ISSN  0022-460X.
  15. ^ Boluriann, S.およびMorris, PJ; Aeroacoustics, 2003, 2(3), 255-292.
  16. ^ Bruus, Henrik (2012). 「音響流体工学 2:摂動理論と超音波共鳴モード」. Lab Chip . 12 (1). 英国王立化学協会 (RSC): 20–28 . doi :10.1039/c1lc20770a. ISSN  1473-0197. PMID  22105715.
  17. ^ Bruus, Henrik (2014年12月15日). 「摂動理論と超音波共鳴」.マイクロスケール音響流体工学. 英国王立化学協会. pp.  29– 45. doi :10.1039/9781849737067-00029. ISBN 978-1-84973-671-8
  18. ^ Petersson, Filip; Nilsson, Andreas; Holm, Cecilia; Jönsson, Henrik; Laurell, Thomas (2005). 「層流と音響定在波力による赤血球からの脂質粒子の連続分離」 . Lab Chip . 5 (1): 20– 22. doi :10.1039/B405748C. ISSN  1473-0197. PMID  15616735.
  19. ^ Shi, Jinjie; Huang, Hua; Stratton, Zak; Huang, Yiping; Huang, Tony Jun (2009). 「定在表面音響波(SSAW)によるマイクロ流体チャネルにおける連続粒子分離」. Lab on a Chip . 9 (23). Royal Society of Chemistry (RSC): 3354– 3359. doi :10.1039/b915113c. ISSN  1473-0197. PMID  19904400.
  20. ^ Nam, Jeonghun; Lim, Hyunjung; Kim, Dookon; Shin, Sehyun (2011). 「マイクロチャネル内の定在表面音響波を用いた全血からの血小板分離」. Lab on a Chip . 11 (19). Royal Society of Chemistry (RSC): 3361– 3364. doi :10.1039/c1lc20346k. ISSN  1473-0197. PMID  21842070.
  21. ^ Ai, Ye; Sanders, Claire K.; Marrone, Babetta L. (2013-09-09). 「定在表面音響波を用いた末梢血単核細胞からの大腸菌(Escherichia coli)の分離」.分析化学. 85 (19). アメリカ化学会 (ACS): 9126– 9134. doi : 10.1021/ac4017715 . ISSN  0003-2700. PMC 3789253. PMID  23968497 . 
  22. ^ Nam, Jeonghun; Lim, Hyunjung; Kim, Choong; Yoon Kang, Ji; Shin, Sehyun (2012). 「定在表面音響波を用いたマイクロ流体チャネル内のカプセル化細胞の密度依存分離」. Biomicrofluidics . 6 (2). AIP Publishing: 24120– 2412010. doi : 10.1063/1.4718719 . ISSN  1932-1058. PMC 3365908. PMID 22670167  . 
  23. ^ Shi, Jinjie; Mao, Xiaole; Ahmed, Daniel; Colletti, Ashley; Huang, Tony Jun (2008). 「定在表面音響波(SSAW)を用いたマイクロ流体チャネル内の微粒子の集束」. Lab Chip . 8 (2). Royal Society of Chemistry (RSC): 221– 223. doi :10.1039/b716321e. ISSN  1473-0197. PMID  18231658.
  24. ^ Shi, Jinjie; Yazdi, Shahrzad; Steven Lin, Sz-Chin; Ding, Xiaoyun; Chiang, I-Kao; Sharp, Kendra; Huang, Tony Jun (2011). 「定在表面音響波(SSAW)によるマイクロ流体チャネルにおける3次元連続粒子集束」. Lab on a Chip . 11 (14). Royal Society of Chemistry (RSC): 2319–24 . doi :10.1039/c1lc20042a. ISSN  1473-0197. PMC 3997299. PMID 21709881  . 
  25. ^ Li, Sixing; Ding, Xiaoyun; Guo, Feng; Chen, Yuchao; Lapsley, Michael Ian; et al. (2013-05-23). 「定在表面音響波を用いたオンチップ・マルチチャンネル液滴ソーター」.分析化学. 85 (11). アメリカ化学会 (ACS): 5468– 5474. doi :10.1021/ac400548d. ISSN  0003-2700. PMC 3988909. PMID  23647057 . 
  26. ^ Ding, Xiaoyun; Lin, Sz-Chin Steven; Lapsley, Michael Ian; Li, Sixing; Guo, Xiang; Chan, Chung Yu; Chiang, I-Kao; Wang, Lin; McCoy, J. Philip; Huang, Tony Jun (2012). 「定在表面弾性波(SSAW)に基づくマルチチャネル細胞選別」. Lab on a Chip . 12 (21). Royal Society of Chemistry (RSC): 4228–31 . doi :10.1039/c2lc40751e. ISSN  1473-0197. PMC 3956451. PMID  22992833 . 
  27. ^ ab Gedge, Michael; Hill, Martyn (2012). 「Acoustofluidics 17: 粒子操作のための表面弾性波デバイスの理論と応用」(PDF) . Lab on a Chip . 12 (17). 英国王立化学協会 (RSC): 2998– 3007. doi :10.1039/c2lc40565b. ISSN  1473-0197. PMID  22842855.
  28. ^ Ding, X.; Lin, S.-CS; Kiraly, B.; Yue, H.; Li, S.; Chiang, I.-K.; Shi, J.; Benkovic, SJ; Huang, TJ (2012-06-25). 「表面音響波を用いた単一微粒子、細胞、および生物のオンチップ操作」. Proceedings of the National Academy of Sciences . 109 (28): 11105– 11109. Bibcode :2012PNAS..10911105D. doi : 10.1073/pnas.1209288109 . ISSN  0027-8424. PMC 3396524. PMID 22733731  . 
  29. ^ Courtney, Charles RP; Demore, Christine EM; Wu, Hongxiao; Grinenko, Alon; Wilcox, Paul D.; Cochran, Sandy; Drinkwater, Bruce W. (2014-04-14). 「器用な音響ピンセットを用いた微粒子の独立した捕捉と操作」. Applied Physics Letters . 104 (15). AIP Publishing: 154103. Bibcode :2014ApPhL.104o4103C. doi : 10.1063/1.4870489 . ISSN  0003-6951.
  30. ^ Meng, Long; Cai, Feiyan; Chen, Juanjuan; Niu, Lili; Li, Yanming; Wu, Junru; Zheng, Hairong (2012-04-23). 「2次元定在表面音響波によるマイクロバブルの精密かつプログラム可能な操作」. Applied Physics Letters . 100 (17). AIP Publishing: 173701. Bibcode :2012ApPhL.100q3701M. doi :10.1063/1.4704922. ISSN  0003-6951.
  31. ^ Wood, CD; Cunningham, JE; O'Rorke, R.; Wälti, C.; Linfield, EH; Davies, AG; Evans, SD (2009-02-02). 「表面音響波によるマイクロ流体システムにおけるミクロンスケール粒子の2次元アレイの形成と操作」. Applied Physics Letters . 94 (5). AIP Publishing: 054101. Bibcode :2009ApPhL..94e4101W. doi :10.1063/1.3076127. ISSN  0003-6951.
  32. ^ シッターズ、ゲリット;カムズマ、ダウエ。タルハマー、グレゴール。リッチ・マルテ、モニカ。ピーターマン、アーウィン JG。 Wuite、Gijs JL (2014-11-24)。 「音響力分光法」。ネイチャーメソッド12 (1)。 Springer Science and Business Media LLC: 47–50 . doi :10.1038/nmeth.3183。ISSN  1548-7091。PMID  25419961。S2CID 12886472  。
  33. ^ Liu, Yan Jun; Ding, Xiaoyun; Lin, Sz-Chin Steven; Shi, Jinjie; Chiang, I-Kao; Huang, Tony Jun (2011-03-14). 「ポリマー分散液晶を用いた表面弾性波駆動型光シャッター」. Advanced Materials . 23 (14). Wiley: 1656– 1659. Bibcode :2011AdM....23.1656L. doi :10.1002/adma.201003708. ISSN  0935-9648. PMID  21438028. S2CID  205238981.
  34. ^ Chen, Yuchao; Ding, Xiaoyun; Steven Lin, Sz-Chin; Yang, Shikuan; Huang, Po-Hsun; et al. (2013-04-09). 「定在表面弾性波を用いたチューナブルナノワイヤパターニング」. ACS Nano . 7 (4). American Chemical Society (ACS): 3306– 3314. doi :10.1021/nn4000034. ISSN  1936-0851. PMC 3989880. PMID  23540330 . 
  35. ^ Lee, Jungwoo; Teh, Shia-Yen; Lee, Abraham; Kim, Hyung Ham; Lee, Changyang; Shung, K. Kirk (2009-08-17). 「シングルビーム音響トラッピング」. Applied Physics Letters . 95 (7). AIP Publishing: 073701. Bibcode :2009ApPhL..95g3701L. doi :10.1063/1.3206910. ISSN  0003-6951. PMC 2755305. PMID 19798424  . 
  36. ^ ab Courtney, Charles RP; Demore, Christine EM; Wu, Hongxiao; Grinenko, Alon; Wilcox, Paul D.; Cochran, Sandy; Drinkwater, Bruce W. (2014-04-14). 「器用な音響ピンセットを用いた微粒子の独立した捕捉と操作」. Applied Physics Letters . 104 (15): 154103. Bibcode :2014ApPhL.104o4103C. doi : 10.1063/1.4870489 . ISSN  0003-6951.
  37. ^ Baresch, Diego; Thomas, Jean-Louis; Marchiano, Régis (2016-01-11). 「弾性粒子のための単一ビーム勾配力音響トラップの観測:音響ピンセット」(PDF) . Physical Review Letters . 116 (2) 024301. American Physical Society (APS). Bibcode :2016PhRvL.116b4301B. doi :10.1103/physrevlett.116.024301. ISSN  0031-9007. PMID  26824541. S2CID  206267537.
  38. ^ ab Marzo, Asier; Seah, Sue Ann; Drinkwater, Bruce W.; Sahoo, Deepak Ranjan; Long, Benjamin; Subramanian, Sriram (2015-10-27). 「浮遊物体の操作のためのホログラフィック音響素子」. Nature Communications . 6 (1): 8661. Bibcode :2015NatCo...6.8661M. doi :10.1038/ncomms9661. ISSN  2041-1723. PMC 4627579. PMID 26505138  . 
  39. ^ Collins, David J.; Devendran, Citsabehsan; Ma, Zhichao; Ng, Jia Wei; Neild, Adrian; Ai, Ye (2016). 「サブタイムオブフライト領域表面音響波による音響ピンセット」. Science Advances . 2 (7) e1600089. アメリカ科学振興協会 (AAAS). Bibcode :2016SciA....2E0089C. doi : 10.1126/sciadv.1600089 . ISSN  2375-2548. PMC 4956186. PMID 27453940  . 
  40. ^ Hefner, Brian T.; Marston, Philip L. (1999-11-23). 「超音波および水中システムのアライメントに応用可能な音響ヘリコイド波トランスデューサー」アメリカ音響学会誌. 106 (6): 3313– 3316. Bibcode :1999ASAJ..106.3313H. doi : 10.1121/1.428184 . ISSN  0001-4966.
  41. ^ Thomas, Jean-Louis; Marchiano, Régis (2003-12-11). 「非線形音響渦における擬似角運動量とトポロジカル電荷保存則」. Physical Review Letters . 91 (24) 244302. Bibcode :2003PhRvL..91x4302T. doi :10.1103/PhysRevLett.91.244302. PMID :  14683126.
  42. ^ Volke-Sepúlveda, Karen; Santillán, Arturo O.; Boullosa, Ricardo R. (2008-01-16). 「自由空間における音響渦から物質への角運動量の伝達」. Physical Review Letters . 100 (2) 024302. Bibcode :2008PhRvL.100b4302V. doi :10.1103/PhysRevLett.100.024302. PMID  18232874.
  43. ^ メルデ, カイ; マーク, アンドリュー G.; チウ, ティアン; フィッシャー, ピア (2016年9月). 「音響の​​ためのホログラム」 . Nature . 537 (7621): 518– 522. Bibcode :2016Natur.537..518M. doi :10.1038/nature19755. ISSN  1476-4687. PMID  27652563. S2CID  4403584.
  44. ^ Jiang, Xue; Li, Yong; Liang, Bin; Cheng, Jian-chun; Zhang, Likun (2016-07-12). 「音響共鳴から軌道角運動量への変換」. Physical Review Letters . 117 (3) 034301. Bibcode :2016PhRvL.117c4301J. doi : 10.1103/PhysRevLett.117.034301 . PMID  27472113.
  45. ^ Jiménez, Noé; Picó, R.; Sánchez-Morcillo, V.; Romero-García, V.; García-Raffi, LM; Staliunas, K. (2016-11-23). 「スパイラル回折格子による高次音響ベッセルビームの形成」. Physical Review E. 94 ( 5) 053004. arXiv : 1604.08353 . Bibcode :2016PhRvE..94e3004J. doi :10.1103/PhysRevE.94.053004. hdl : 2117/101154 . PMID:  27967159. S2CID  : 27190492.
  46. ^ ヒメネス、ノエ;ロメロ=ガルシア、ヴィセント。ガルシア=ラフィ、ルイス・M.カマレナ、フランシスコ。スタリウナス、ケストゥティス (2018-05-14)。 「フレネルスパイラルゾーンプレートによるシャープな音響渦集束」。応用物理学の手紙112 (20): 204101。書誌コード:2018ApPhL.112t4101J。土井:10.1063/1.5029424。hdl : 2117/118785ISSN  0003-6951。
  47. ^ Yeo, Leslie Y.; Friend, James R. (2014-01-03). 「表面音響波マイクロ流体工学」. Annual Review of Fluid Mechanics . 46 (1): 379– 406. Bibcode :2014AnRFM..46..379Y. doi :10.1146/annurev-fluid-010313-141418. ISSN  0066-4189.
  • 高速音響ピンセット -音響ピンセットの仕組みを説明するYouTube動画
Retrieved from "https://en.wikipedia.org/w/index.php?title=Acoustic_tweezers&oldid=1323819119"