Lees-Edwards 境界条件を用いたずり変形の計算によって、粘度を評価することができます。

n-ｄecane とn-ｈexadecane の結果を以下に示します。J-OCTAを用いて作成した全原子モデルを使用し、n-ｄecane は密度0.613[g/cm3 ]で温度480K、n-ｈexadecane は密度0.753[g./cm3] で温度323K の条件で計算を実施しました。系内の分子数はそれぞれ200 と100 です。ずり速度は1e+12 [1/s]から1e+10[1/s]まで変化させながら粘度の評価を実施しました。分子動力学の計算にはVSOP を使用しました (Figure.1)。

Figure.2 に、ずり速度に対して算出された粘度を両対数でプロットした結果を示します。ずり速度の低下と共に粘度が増加しますが、n-decane ではずり速度1e+10[1/s]以下では粘度が一定に近づいていく様子が見えます。

文献[1]でのGreen-Kubo 法での粘度の評価では、同条件でのn-decane が0.178[E-3Pa・s]、
n-hexadecane が1.17[E-3Pa・s]となっています。