Novel numerical method for calculating the pressure tensor in spherical coordinates for molecular systems

DOI PDF 7 Citations 31 References
  • Takenobu Nakamura
    National Institute of Advanced Industrial Science and Technology (AIST) 1 Nanosystem Research Institute (NRI), , Central-2, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
  • Wataru Shinoda
    National Institute of Advanced Industrial Science and Technology (AIST) 1 Nanosystem Research Institute (NRI), , Central-2, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
  • Tamio Ikeshoji
    Tohoku University 3 Advanced Institute for Materials Research (WPI-AIMR), , 2-1-1 Katahira, Aoba-ku, Sendai 980-8579, Japan

Abstract

<jats:p>We propose a novel method for computing the pressure tensor along the radial axis of a molecular system with spherical symmetry. The proposed method uses the slice averaged pressure to improve the numerical stability and precision significantly. Simplified expressions of the local pressure are derived for a conventional molecular force field including non-bond, bond stretching, angle bending, and torsion interactions; these expressions are advantageous in terms of the computational cost. We also discuss an algorithm to avoid numerical singularity. Finally, the method is successfully applied to three different molecular systems, i.e., a water droplet in oil, a spherical micelle, and a liposome.</jats:p>

Journal

Citations (7)*help

See more

References(31)*help

See more

Related Projects

See more

Details 詳細情報について

Report a problem

Back to top