title: Novel numerical method for calculating the pressure tensor in spherical coordinates for molecular systems.
abstract: 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.
authors.affiliation: Nanosystem Research Institute (NRI), National Institute of Advanced Industrial Science and Technology (AIST), Central-2, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan. takenobu.nakamura@aist.go.jp
authors.forename: Takenobu
authors.initials: T
authors.lastname: Nakamura
authors.forename: Wataru
authors.initials: W
authors.lastname: Shinoda
authors.forename: Tamio
authors.initials: T
authors.lastname: Ikeshoji
mesh_terms: D000465:Algorithms; D003198:Computer Simulation; D006586:Hexanes / Q000737:chemistry; D008081:Liposomes / Q000737:chemistry; D008956:Models, Chemical; D008958:Models, Molecular; D009821:Oils / Q000737:chemistry; D011312:Pressure; D014867:Water / Q000737:chemistry
chemical_list: D006586:Hexanes; D008081:Liposomes; D009821:Oils; D014867:Water