Abstract
We introduce a mean-field theoretical framework for generalizing isotropic pair potentials to anisotropic shapes. This method is suitable for generating pair potentials that can be used in both Monte Carlo and molecular dynamics simulations. We demonstrate the application of this theory by deriving a Lennard-Jones (LJ)-like potential for arbitrary geometries along with a Weeks-Chandler-Anderson-like repulsive variant, showing that the resulting potentials behave very similarly to standard LJ potentials while also providing a nearly conformal mapping of the underlying shape. We then describe an implementation of this potential in the simulation engine HOOMD-blue and discuss the challenges that must be overcome to achieve a sufficiently robust and performant implementation. The resulting potential can be applied to smooth geometries like ellipsoids and to convex polytopes. We contextualize these applications with reference to the existing methods for simulating such particles. The pair potential is validated using standard criteria, and its performance is compared to existing methods for comparable simulations. Finally, we show the results of self-assembly simulations, demonstrating that this method can be used to study the assembly of anisotropic particles into crystal structures.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.