1
|
Yonetani Y. Unsolved problem of long-range interactions: dipolar spin-ice study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:175401. [PMID: 38270229 DOI: 10.1088/1361-648x/ad1ca6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/15/2023] [Accepted: 01/09/2024] [Indexed: 01/26/2024]
Abstract
Long-range interactions derive various strange phenomena. As illustrated by cutoff simulations of water, increasing cutoff length does not improve the simulation result necessarily; on the contrary, it makes the result worse. In the extreme situation, the structure of water transforms into a layer structure. In this study, to explore the underlying mechanism of this phenomenon, we performed Monte Carlo simulations on dipolar spins arranged on a pyrochlore spin-ice lattice. Like the water case, the present dipolar spin system also showed cutoff-induced dipole ordering and layer formation. The width of the layers depended on the cutoff length; and longer cutoff length led to a broader layer. These features are certainly consistent with the previous water case. This indicates that layer formation is the general behavior of dipolar systems whose interactions are truncated within a finite distance. The result is important for future exploration of the relationship between long-range interactions and resulting structures. In addition, it emphasizes the necessity of rigorous treatment of long-range interactions because increasing the cutoff length prevents convergence and provides an entirely different result from the rigorous Ewald calculation.
Collapse
Affiliation(s)
- Yoshiteru Yonetani
- Kansai Institute for Photon Science and Institute for Quantum Life Science, National Institutes for Quantum Science and Technology (QST), 8-1-7 Umemidai, Kizugawa, Kyoto, 619-0215, Japan
| |
Collapse
|
2
|
Strange layer structure of dipolar spins formed on the spin-ice lattice. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
|
3
|
Yonetani Y. Dielectric continuum model examination of real-space electrostatic treatments. J Chem Phys 2021; 154:044103. [PMID: 33514106 DOI: 10.1063/5.0033053] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Electrostatic interaction is long ranged; thus, the accurate calculation is not an easy task in molecular dynamics or Monte Carlo simulations. Though the rigorous Ewald method based on the reciprocal space has been established, real-space treatments have recently become an attractive alternative because of the efficient calculation. However, the construction is not yet completed and is now a challenging subject. In an earlier theoretical study, Neumann and Steinhauser employed the Onsager dielectric continuum model to explain how simple real-space cutoff produces artificial dipolar orientation. In the present study, we employ this continuum model to explore the fundamental properties of the recently developed real-space treatments of three shifting schemes. The result of the distance-dependent Kirkwood function GK(R) showed that the simple bare cutoff produces a well-known hole-shaped artifact, whereas the shift treatments do not. Two-dimensional mapping of electric field well explained how these shift treatments remove the hole-shaped artifact. Still, the shift treatments are not sufficient because they do not produce a flat GK(R) profile unlike ideal no-cutoff treatment. To test the continuum model results, we also performed Monte Carlo simulations of dipolar particles. The results found that the continuum model could predict the qualitative tendency as to whether each electrostatic treatment produces the hole-shaped artifact of GK(R) or not. We expect that the present study using the continuum model offers a stringent criterion to judge whether the primitive electrostatic behavior is correctly described or not, which will be useful for future construction of electrostatic treatments.
Collapse
Affiliation(s)
- Yoshiteru Yonetani
- Quantum Beam Science Research Directorate and Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| |
Collapse
|
4
|
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
|
5
|
Nguyen D, Macchi P, Volkov A. Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. III. Application to crystal structures via the Ewald and direct summation methods. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2020; 76:630-651. [PMID: 33125348 DOI: 10.1107/s2053273320009584] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2020] [Accepted: 07/13/2020] [Indexed: 12/28/2022]
Abstract
The previously reported exact potential and multipole moment (EP/MM) method for fast and accurate evaluation of the intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density [Volkov, Koritsanszky & Coppens (2004). Chem. Phys. Lett. 391, 170-175; Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536; Nguyen & Volkov (2019). Acta Cryst. A75, 448-464] is extended to the calculation of electrostatic interaction energies in molecular crystals using two newly developed implementations: (i) the Ewald summation (ES), which includes interactions up to the hexadecapolar level and the EP correction to account for short-range electron-density penetration effects, and (ii) the enhanced EP/MM-based direct summation (DS), which at sufficiently large intermolecular separations replaces the atomic multipole moment approximation to the electrostatic energy with that based on the molecular multipole moments. As in the previous study [Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536], the EP electron repulsion integral is evaluated analytically using the Löwdin α-function approach. The resulting techniques, incorporated in the XDPROP module of the software package XD2016, have been tested on several small-molecule crystal systems (benzene, L-dopa, paracetamol, amino acids etc.) and the crystal structure of a 181-atom decapeptide molecule (Z = 4) using electron densities constructed via the University at Buffalo Aspherical Pseudoatom Databank [Volkov, Li, Koritsanszky & Coppens (2004). J. Phys. Chem. A, 108, 4283-4300]. Using a 2015 2.8 GHz Intel Xeon E3-1505M v5 computer processor, a 64-bit implementation of the Löwdin α-function and one of the higher optimization levels in the GNU Fortran compiler, the ES method evaluates the electrostatic interaction energy with a numerical precision of at least 10-5 kJ mol-1 in under 6 s for any of the tested small-molecule crystal structures, and in 48.5 s for the decapeptide structure. The DS approach is competitive in terms of precision and speed with the ES technique only for crystal structures of small molecules that do not carry a large molecular dipole moment. The electron-density penetration effects, correctly accounted for by the two described methods, contribute 28-64% to the total electrostatic interaction energy in the examined systems, and thus cannot be neglected.
Collapse
Affiliation(s)
- Daniel Nguyen
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Via Mancinelli 7, Milano 20131, Italy
| | - Anatoliy Volkov
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| |
Collapse
|
6
|
Stenqvist B, Aspelin V, Lund M. Generalized Moment Correction for Long-Ranged Electrostatics. J Chem Theory Comput 2020; 16:3737-3745. [PMID: 32315176 PMCID: PMC7588037 DOI: 10.1021/acs.jctc.9b01003] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
![]()
Describing
long-ranged electrostatics using short-ranged pair potentials
is appealing because the computational complexity scales linearly
with the number of particles. The foundation of the approach presented
here is to mimic the long-ranged medium response by cancelling electric
multipoles within a small cutoff sphere. We propose a rigorous and
formally exact new method that cancels up to infinitely many multipole moments and is free of operational damping parameters
often required in existing theories. Using molecular dynamics simulations
of water with and without added salt, we discuss radial distribution
functions, Kirkwood–Buff integrals, dielectrics, diffusion
coefficients, and angular correlations in relation to existing electrostatic
models. We find that the proposed method is an efficient and accurate
alternative for handling long-ranged electrostatics as compared to
Ewald summation schemes. The methodology and proposed parameterization
are applicable also for dipole–dipole interactions.
Collapse
Affiliation(s)
- Björn Stenqvist
- Division of Physical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden.,Division of Theoretical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | - Vidar Aspelin
- Division of Theoretical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | - Mikael Lund
- Division of Theoretical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden.,LINXS-Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, SE-223 70 Lund, Sweden
| |
Collapse
|
7
|
Cuetos A, Morillo N, Martı Nez Haya B. Coadsorption of Counterionic Colloids at Fluid Interfaces: A Coarse-Grained Simulation Study of Gibbs Monolayers. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:2877-2885. [PMID: 32118442 DOI: 10.1021/acs.langmuir.9b03886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Monolayers of oppositely charged colloids form versatile self-organizing substrates, with a recognized potential to tailor functional interfaces. In this study, a coarse-grained Monte Carlo simulation approach is laid out to assess the structural properties of Gibbs monolayers, in which one of the counterionic species is partially soluble. It is shown that the composition of this type of monolayer varies in a nontrivial way with surface coverage, as a result of a subtle competition between steric and attractive forces. In the regime of weak electrostatic interactions, the monolayer is depleted of soluble colloids as the surface coverage is increased. At sufficiently strong interactions, the incorporation of soluble colloids is favored at high surface coverage, leading to a re-entrant-type behavior in the expansion/compression isotherms. Strong electrostatic interactions also favor the clustering of the colloids, leading to a range of aggregated configurations, qualitatively resembling those obtained in previous experimental studies. At sufficiently high surface coverage, the clusters collapse into a gel-like percolated mesoscopic structure and eventually into a square crystal lattice configuration. Such interfacial structures are in good agreement with the ones observed in the few experimental investigations available for these systems, showing that the simple methodology introduced in this study provides a valuable predictive framework to anticipate the landscape of interfacial structures that may be produced with oppositely charged colloids, through the modulation of pair interactions and thermodynamical conditions.
Collapse
Affiliation(s)
- Alejandro Cuetos
- Department of Physical, Chemical and Natural Systems, Universidad Pablo de Olavide, 41013 Seville, Spain
| | - Neftali Morillo
- Department of Physical, Chemical and Natural Systems, Universidad Pablo de Olavide, 41013 Seville, Spain
| | - Bruno Martı Nez Haya
- Department of Physical, Chemical and Natural Systems, Universidad Pablo de Olavide, 41013 Seville, Spain
| |
Collapse
|
8
|
Stenqvist B, Lund M. On short-ranged pair-potentials for long-range electrostatics. Phys Chem Chem Phys 2019; 21:24787-24792. [DOI: 10.1039/c9cp03875b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Fast and accurate summation of long-range electrostatics by using a short-ranged pair-potential that ensures moment cancellation in the cutoff sphere.
Collapse
Affiliation(s)
- Björn Stenqvist
- Division of Physical Chemistry
- Lund University
- SE-22100 Lund
- Sweden
| | - Mikael Lund
- Division of Theoretical Chemistry
- Lund University
- SE-22100 Lund
- Sweden
- LINXS – Lund Institute of Advanced Neutron and X-ray Science
| |
Collapse
|
9
|
Waibel C, Feinler MS, Gross J. A Modified Shifted Force Approach to the Wolf Summation. J Chem Theory Comput 2018; 15:572-583. [PMID: 30418767 DOI: 10.1021/acs.jctc.8b00343] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Wolf method for calculation of electrostatic interactions in molecular simulations is known to describe the energy well, whereas the forces have discontinuities. For a more reliable description of the forces this method can be extended with a shifted force approach. This leads to a good description of the forces and precise molecular dynamics simulation, but the description of the energy becomes poorer. In this study we propose a modification of a shifted force extension to describe the energy as well as the forces in better agreement to reference data as determined from the Ewald summation. We show that vapor-liquid phase equilibria (VLE) calculated with Monte Carlo simulations in the grand canonical ensemble and dynamic properties calculated with molecular dynamics simulations can be calculated reliably using this modification to describe the electrostatic interactions.
Collapse
Affiliation(s)
- Christian Waibel
- Institute of Thermodynamics and Thermal Process Engineering , University of Stuttgart , Pfaffenwaldring 9 , 70569 Stuttgart , Germany
| | - Mathias Simon Feinler
- Institute of Thermodynamics and Thermal Process Engineering , University of Stuttgart , Pfaffenwaldring 9 , 70569 Stuttgart , Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering , University of Stuttgart , Pfaffenwaldring 9 , 70569 Stuttgart , Germany
| |
Collapse
|
10
|
Waibel C, Gross J. Modification of the Wolf Method and Evaluation for Molecular Simulation of Vapor–Liquid Equilibria. J Chem Theory Comput 2018; 14:2198-2206. [DOI: 10.1021/acs.jctc.7b01190] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Christian Waibel
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
11
|
Lamichhane M, Parsons T, Newman KE, Gezelter JD. Real space electrostatics for multipoles. III. Dielectric properties. J Chem Phys 2016; 145:074108. [PMID: 27544088 DOI: 10.1063/1.4960957] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In Papers I and II, we developed new shifted potential, gradient shifted force, and Taylor shifted force real-space methods for multipole interactions in condensed phase simulations. Here, we discuss the dielectric properties of fluids that emerge from simulations using these methods. Most electrostatic methods (including the Ewald sum) require correction to the conducting boundary fluctuation formula for the static dielectric constants, and we discuss the derivation of these corrections for the new real space methods. For quadrupolar fluids, the analogous material property is the quadrupolar susceptibility. As in the dipolar case, the fluctuation formula for the quadrupolar susceptibility has corrections that depend on the electrostatic method being utilized. One of the most important effects measured by both the static dielectric and quadrupolar susceptibility is the ability to screen charges embedded in the fluid. We use potentials of mean force between solvated ions to discuss how geometric factors can lead to distance-dependent screening in both quadrupolar and dipolar fluids.
Collapse
Affiliation(s)
- Madan Lamichhane
- Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Thomas Parsons
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Kathie E Newman
- Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - J Daniel Gezelter
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| |
Collapse
|