1
|
Fukuda I, Moritsugu K, Higo J, Fukunishi Y. A cutoff-based method with charge-distribution-data driven pair potentials for efficiently estimating electrostatic interactions in molecular systems. J Chem Phys 2023; 159:234116. [PMID: 38112509 DOI: 10.1063/5.0172270] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 11/27/2023] [Indexed: 12/21/2023] Open
Abstract
We introduce a simple cutoff-based method for precise electrostatic energy calculations in the molecular dynamics (MD) simulations of point-particle systems. Our method employs a theoretically derived smooth pair potential function to define electrostatic energy, offering stability and computational efficiency in MD simulations. Instead of imposing specific physical conditions, such as dielectric environments or charge neutrality, we focus on the relationship represented by a single summation formula of charge-weighted pair potentials. This approach allows an accurate energy approximation for each particle, enabling a straightforward error analysis. The resulting particle-dependent pair potential captures the charge distribution information, making it suitable for heterogeneous systems and ensuring an enhanced accuracy through distant information inclusion. Numerical investigations of the Madelung constants of crystalline systems validate the method's accuracy.
Collapse
Affiliation(s)
- Ikuo Fukuda
- Graduate School of Science, Osaka Metropolitan University, 1-2 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8231, Japan
- Japan Biological Informatics Consortium, 2-4-32 Aomi, Koto-ku, Tokyo, 135-8073, Japan
| | - Kei Moritsugu
- Graduate School of Science, Osaka Metropolitan University, 1-2 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8231, Japan
| | - Junichi Higo
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima Minamimachi, Chuo-ku, Kobe, Hyogo650-0047, Japan
- Research Organization of Science and Technology, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 2-3-26, Aomi, Koto-ku, Tokyo 135-0064, Japan
| |
Collapse
|
2
|
Brumby PE, Kowaguchi A, Nozawa T, Yasuoka K, Wensink HH. Pre-Smectic Ordering and the Unwinding Helix in Monte Carlo Simulations of Cholesteric Liquid-Crystals. J Phys Chem B 2023; 127:7194-7204. [PMID: 37540189 DOI: 10.1021/acs.jpcb.3c02018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/05/2023]
Abstract
In this paper, molecular chirality is studied for liquid-crystal fluids represented by hard rods with the addition of an attractive chiral dispersion term. Chiral forces between molecular pairs are assumed to be long-ranged and are described in terms of the pseudotensor of Goossens [W. J. A. Goossens, Mol. Cryst. Liq. Cryst. 1971, 12, 237-244]. Following Varga and Jackson [S. Varga and G. Jackson, Chem. Phys. Lett. 2003, 377, 6-12], this is combined with a hard-spherocylinder core. We investigate the relationship between molecular chirality and the helical pitch of the system, which occurs in the absence of full three-dimensional periodic boundary conditions. The dependence of the wavenumber of this pitch on the thermodynamic variables, temperature, and density is measured. We also explore the use of a novel surface boundary interaction model. As a result of this approach, we are able to lower the temperature of the system without the occurrence of nematic droplets, which would interfere with the formation of a uniaxial pitch. Regarding the theoretical predictions of Wensink and Jackson [H. H. Wensink and G. Jackson, J. Chem. Phys. 2009, 130, 234911], on the one hand, we have qualitative agreement with the observed non-monotonic density dependence of the wavenumber. Initially increasing with density, the wavenumber reaches a maximum, before falling as the density moves toward the point of phase transition from cholesteric to smectic. However, further analysis for shorter rods, in the presence of novel boundary conditions, reveals some disagreement with the theory, at least in this case; the unwinding of the cholesteric helix in the cholesteric phase occurs simultaneously with subtle increases in smectic ordering. These pre-smectic fluctuations have not been accounted for so far in theories on cholesterics but turn out to play a key role in controlling the pitch of cholesteric phases of rod-shaped mesogens with a small to moderate aspect ratio.
Collapse
Affiliation(s)
- Paul E Brumby
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Akie Kowaguchi
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Takuma Nozawa
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Henricus H Wensink
- Laboratoire de Physique des Solides─UMR 8502, Université Paris-Saclay & CNRS, Orsay 91405, France
| |
Collapse
|
3
|
Renner J, Schmidt M, de Las Heras D. Reduced-variance orientational distribution functions from torque sampling. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:235901. [PMID: 36974000 DOI: 10.1088/1361-648x/acc522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Accepted: 03/08/2023] [Indexed: 06/18/2023]
Abstract
We introduce a method to sample the orientational distribution function in computer simulations. The method is based on the exact torque balance equation for classical many-body systems of interacting anisotropic particles in equilibrium. Instead of the traditional counting of events, we reconstruct the orientational distribution function via an orientational integral of the torque acting on the particles. We test the torque sampling method in two- and three-dimensions, using both Langevin dynamics and overdamped Brownian dynamics, and with two interparticle interaction potentials. In all cases the torque sampling method produces profiles of the orientational distribution function with better accuracy than those obtained with the traditional counting method. The accuracy of the torque sampling method is independent of the bin size, and hence it is possible to resolve the orientational distribution function with arbitrarily small angular resolutions.
Collapse
Affiliation(s)
- Johannes Renner
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
4
|
Eckert T, Schmidt M, de Las Heras D. Effect of sample height and particle elongation in the sedimentation of colloidal rods. SOFT MATTER 2023; 19:2214-2223. [PMID: 36883340 DOI: 10.1039/d3sm00191a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
We study theoretically the effect of a gravitational field on the equilibrium behaviour of a colloidal suspension of rods with different length-to-width aspect ratios. The bulk phases of the system are described with analytical equations of state. The gravitational field is then incorporated via sedimentation path theory, which assumes a local equilibrium condition at each altitude of the sample. The bulk phenomenology is significantly enriched by the presence of the gravitational field. In a suspension of elongated rods with five stable phases in bulk, the gravitational field stabilizes up to fifteen different stacking sequences. The sample height has a non-trivial effect on the stable stacking sequence. New layers of distinct bulk phases appear either at the top, at the bottom, or simultaneously at the top and the bottom when increasing the sample height at constant colloidal concentration. We also study sedimentation in a mass-polydisperse suspension in which all rods have the same shape but different buoyant masses.
Collapse
Affiliation(s)
- Tobias Eckert
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany.
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany.
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany.
| |
Collapse
|
5
|
Fukuda I, Nakamura H. Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference. Biophys Rev 2022; 14:1315-1340. [PMID: 36659982 PMCID: PMC9842848 DOI: 10.1007/s12551-022-01029-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2022] [Accepted: 11/30/2022] [Indexed: 01/13/2023] Open
Abstract
In molecular simulations, it is essential to properly calculate the electrostatic interactions of particles in the physical system of interest. Here we consider a method called the non-Ewald method, which does not rely on the standard Ewald method with periodic boundary conditions, but instead relies on the cutoff-based techniques. We focus on the physicochemical and mathematical conceptual aspects of the method in order to gain a deeper understanding of the simulation methodology. In particular, we take into account the reaction field (RF) method, the isotropic periodic sum (IPS) method, and the zero-multipole summation method (ZMM). These cutoff-based methods are based on different physical ideas and are completely distinguishable in their underlying concepts. The RF and IPS methods are "additive" methods that incorporate information outside the cutoff region, via dielectric medium and isotropic boundary condition, respectively. In contrast, the ZMM is a "subtraction" method that tries to remove the artificial effects, generated near the boundary, from the cutoff sphere. Nonetheless, we find physical and/or mathematical similarities between these methods. In particular, the modified RF method can be derived by the principle of neutralization utilized in the ZMM, and we also found a direct relationship between IPS and ZMM.
Collapse
Affiliation(s)
- Ikuo Fukuda
- Graduate School of Information Science, University of Hyogo, 7-1-28 Minatojima, Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047 Japan
| | - Haruki Nakamura
- Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871 Japan
| |
Collapse
|
6
|
Falcón-González JM, Contreras-Aburto C, Lara-Peña M, Heinen M, Avendaño C, Gil-Villegas A, Castañeda-Priego R. Assessment of the Wolf method using the Stillinger-Lovett sum rules: From strong electrolytes to weakly charged colloidal dispersions. J Chem Phys 2020; 153:234901. [PMID: 33353329 DOI: 10.1063/5.0033561] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The Ewald method has been the cornerstone in molecular simulations for modeling electrostatic interactions of charge-stabilized many-body systems. In the late 1990s, Wolf and collaborators developed an alternative route to describe the long-range nature of electrostatic interactions; from a computational perspective, this method provides a more efficient and straightforward way to implement long-range electrostatic interactions than the Ewald method. Despite these advantages, the validity of the Wolf potential to account for the electrostatic contribution in charged fluids remains controversial. To alleviate this situation, in this contribution, we implement the Wolf summation method to both electrolyte solutions and charged colloids with moderate size and charge asymmetries in order to assess the accuracy and validity of the method. To this end, we verify that the proper selection of parameters within the Wolf method leads to results that are in good agreement with those obtained through the standard Ewald method and the theory of integral equations of simple liquids within the so-called hypernetted chain approximation. Furthermore, we show that the results obtained with the original Wolf method do satisfy the moment conditions described by the Stillinger-Lovett sum rules, which are directly related to the local electroneutrality condition and the electrostatic screening in the Debye-Hückel regime. Hence, the fact that the solution provided by the Wolf method satisfies the first and second moments of Stillinger-Lovett proves, for the first time, the reliability of the method to correctly incorporate the electrostatic contribution in charge-stabilized fluids. This makes the Wolf method a powerful alternative compared to more demanding computational approaches.
Collapse
Affiliation(s)
- José Marcos Falcón-González
- Unidad Profesional Interdisciplinaria de Ingeniería, Campus Guanajuato, Instituto Politécnico Nacional, Av. Mineral de Valenciana No. 200, Col. Fraccionamiento Industrial Puerto Interior, C.P. 36275 Silao de la Victoria, Guanajuato, Mexico
| | - Claudio Contreras-Aburto
- Facultad de Ciencias en Física y Matemáticas, Universidad Autónoma de Chiapas, 29050 Tuxtla Gutiérrez, Mexico
| | - Mayra Lara-Peña
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Mexico
| | - Marco Heinen
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Mexico
| | - Carlos Avendaño
- Department of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester M13 9PL, United Kingdom
| | - Alejandro Gil-Villegas
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Mexico
| | - Ramón Castañeda-Priego
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Mexico
| |
Collapse
|
7
|
Viveros-Méndez PX, Gil-Villegas A, Aranda Espinoza S. Assessment by Monte Carlo computer simulations of the phase behavior of hard spherocylinders confined within cylindrical cavities. J Chem Phys 2017; 147:234902. [DOI: 10.1063/1.5017844] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Perla X. Viveros-Méndez
- Unidad Académica de Física, Universidad Autónoma de Zacatecas, Calzada Solidaridad esq. Paseo, La Bufa s/n, 98060 Zacatecas, Zacatecas, México
| | - Alejandro Gil-Villegas
- Departamento de Ingeniería Física, División de Ciencias e Ingenierías Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Guanajuato, México
| | - Said Aranda Espinoza
- Instituto de Física ’Manuel Sandoval Vallarta’, Universidad Autónoma de San Luis Potosí, Álvaro Obregon 64, 78000 San Luis Potosí, SLP, México
| |
Collapse
|