ビリアル応力のソースを表示

{{孤立|date=2021年3月}}
{{翻訳直後|[[:en:Special:Redirect/revision/950409155|Virial stress]]|date=2020年12月}}

'''ビリアル応力'''（ビリアルおうりょく、{{Lang-en-short|virial stress}}）は、均一系に対して用いられる[[原子]]スケールの機械的[[応力]]の尺度の一つ。 局所的なビリアル応力の表式は、分子系の自由エネルギーを[[変形]]テンソルについて汎関数微分したものとして導出できる<ref name="morante06">Morante, S., G. C. Rossi, and M. Testa. "The stress tensor of a molecular system: An exercise in statistical mechanics." The Journal of chemical physics 125.3 (2006): 034101, http://aip.scitation.org/doi/abs/10.1063/1.2214719.
</ref>。

== 体積平均ビリアル応力の定義 ==
体積平均ビリアル応力の瞬時値は次の式により与えられる。

: <math>\tau_{ij} = \frac{1}{\Omega} \sum_{k \in \Omega} \left(-m^{(k)} (u_i^{(k)}- \bar{u}_i) (u_j^{(k)}- \bar{u}_j) + \frac{1}{2} \sum_{\ell \in \Omega} ( x_i^{(\ell)} - x_i^{(k)}) f_j^{(k\ell)}\right)</math>

ここで、以下の[[変数 (数学)|変数]]を用いた。

* ''{{Mvar|k}}'' および {{Mvar|&ell;}} ：領域内の原子を表わす添字
* {{Math|&Omega;}} ：領域の[[体積]]
* {{Math|''m''{{sup|(''k'')}}}}：原子''{{Mvar|k}}''の[[質量]]
* {{Math|''u''{{su|p=(''k'')|b=''i''}}}}：原子''{{Mvar|k}}の''[[速度]]ベクトルの''{{Mvar|i}}''成分
* {{Math|{{overline|''u''}}{{sub|''j''}}}} ：領域内における平均原子速度ベクトルの{{Mvar|j}}成分
* {{Math|''x''{{su|p=(''k'')|b=''i''}}}}：原子''{{Mvar|k}}''の[[位置]]ベクトルの{{Mvar|i}}成分
* {{Math|''f''{{su|p=(''k&ell;'')|b=''i''}}}}：原子{{Mvar|k}}に別の原子{{Mvar|&ell;}}から及ぼされる[[力 (物理学)|力]]の{{Mvar|i}}成分

[[絶対零度]]における表式は、すべての速度がゼロになるため以下の通りに単純化できる。

: <math>\tau_{ij} = \frac{1}{2\Omega} \sum_{k,\ell \in \Omega} ( x_i^{(\ell)} - x_i^{(k)}) f_j^{(k\ell)}</math>

この量は次のように解釈することができる。{{Math|''&tau;''{{sub|11}}}}成分は {{Math|''x''{{sub|1}}}} 方向の力を、この方向に垂直な面の面積で割った値である。このような平面で区切られた2つの隣接する領域について考えると、その境界面における応力の11成分は、面の両側にある原子間に働く相互作用の総計となる。

したがって、体積平均ビリアル応力は、ビリアル応力の瞬時値を体積平均したものの{{仮リンク|アンサンブル平均 (統計力学)|en|Ensemble average (statistical mechanics)|label=アンサンブル平均}}である。

3次元等方系の平衡状態においては、原子スケールの「瞬時」圧力は、負の応力テンソルの対角成分の平均により定義される。

: <math>\mathcal{P}_{at} = -\frac{1}{3} \operatorname{Tr}(\tau)</math>

そして、圧力は瞬時圧力のアンサンブル平均により定義される<ref name="allen91">{{Cite book|last=Allen|first=MP|last2=Tildesley|first2=DJ|year=1991|title=Computer Simulations of Liquids|editor=Clarendon Press|location=Oxford|pages=46–50}}</ref>。

: <math>P_{at} =\langle \mathcal{P}_{at} \rangle</math>

これが領域 {{Math|&Omega;}} 内の平均圧力である。

=== 等価な別の定義 ===
一部の論文や教科書<ref name="allen91">{{Cite book|last=Allen|first=MP|last2=Tildesley|first2=DJ|year=1991|title=Computer Simulations of Liquids|editor=Clarendon Press|location=Oxford|pages=46–50}}</ref>ではわずかに異なるが同等の定義式を使用していることに注意が必要である。

: <math>\tau_{ij} = \frac{1}{\Omega} \sum_{k \in \Omega} \left(-m^{(k)} (u_i^{(k)}- \bar{u}_i) (u_j^{(k)}- \bar{u}_j) - \frac{1}{2} \sum_{\ell \in \Omega} x_i^{(k\ell)} f_j^{(k\ell)}\right)</math>

ここで、{{Math|''x''{{su|p=(''k&ell;'')|b=''i''}}}} は{{Mvar|&ell;}}番目の原子からみた{{Mvar|k}}番目の原子の位置ベクトルの{{Mvar|i}}成分である。

: <math> x_i^{k\ell} = x_i^{(k)} - x_i^{(\ell)} </math>

これら二つの方程式は厳密に同等であるが、ベクトルの定義が混乱を招く可能性がある。

=== 導出 ===
[[ビリアル定理]]を用い、粒子間に働く力と容器との間の働く力を分割することによりビリアル応力は導出できる<ref>{{Cite journal|last=Navet|first=M.|last2=Jamin|first2=E.|last3=Feix|first3=M. R.|date=1980-02-01|title=" Virial " pressure of the classical one-component plasma|url=https://hal.archives-ouvertes.fr/jpa-00231724|journal=Journal de Physique Lettres|volume=41|issue=3|pages=69–73|language=en|DOI=10.1051/jphyslet:0198000410306900|ISSN=0302-072X}}</ref>。あるいは、圧力の定義式 <math>P=-\frac{\partial F(N,V,T)}{\partial V}</math> を直接用い、スケーリングされた座標を用いて計算することによっても導出できる。

== 不均一系 ==
系がある所与の領域において均一ではない場合、上記の体積平均された圧力は適切な尺度ではなくなる。不均一系における圧力は位置と方向に依存して変化する。したがって、不均一系においては、局所的な圧力を定義する必要がある<ref>Numerical Simulations of a Smectic Lamellar Phase of Amphiphilic Molecules, p. 40, https://books.google.de/books?id=rPpegGthzO4C&lpg=PA40&dq=local%20pressure%20tensor&hl=de&pg=PA40#v=onepage&q=local%20pressure%20tensor&f=false</ref>。圧力が不均一な系の一般的な例として、地球の大気中の圧力が高度に依存することが上げられる。

== 局所ビリアル応力の瞬時値 ==
局所ビリアル応力の瞬時値は次の式により与えられる<ref name="morante06">Morante, S., G. C. Rossi, and M. Testa. "The stress tensor of a molecular system: An exercise in statistical mechanics." The Journal of chemical physics 125.3 (2006): 034101, http://aip.scitation.org/doi/abs/10.1063/1.2214719.
</ref>。

<math>\tau_{ab}(\vec{r})=- \sum_{i=1}^N \delta(\vec{r}-\vec{r}^{(i)}) \left(m^{(i)} u^{(i)}_a u^{(i)}_b + \frac{1}{2} \sum_{j=1, j \neq i}^{N} (\vec{r}^{(i)}-\vec{r}^{(j)})_a \vec{f}^{(ij)}_b \right)</math>

== 分子シミュレーションにおけるビリアル圧力の測定 ==
ビリアル圧力は、上記の式を使用するか、体積をスケーリングして再計算することにより算出される<ref>{{Cite journal|last=Miguel|first=Enrique de|last2=Jackson|first2=George|date=2006-10-30|title=The nature of the calculation of the pressure in molecular simulations of continuous models from volume perturbations|journal=The Journal of Chemical Physics|volume=125|issue=16|pages=164109|language=en|bibcode=2006JChPh.125p4109D|DOI=10.1063/1.2363381|ISSN=0021-9606|PMID=17092065}}</ref>。

== 関連項目 ==

* [[ビリアル定理]]

== 参考文献 ==
{{Reflist}}
[[Category:連続体力学]]
{{DEFAULTSORT:ひりあるおうりよく}}