Mitigating Freezing Artifacts in Polar Coarse-Grained Water through Rebalancing Electrostatic and van der Waals Interactions: A Case Study of the pSPICA Force Field.
5 affiliations
- 1Department of Chemical Engineering, National Cheng Kung University, Tainan701, Taiwan.
- 2Research Institute for Interdisciplinary Science, Okayama University, Okayama700-8530, Japan.
- 3Center for Resilience and Intelligence on Sustainable Energy Research (RiSER), National Cheng Kung University, Tainan701, Taiwan.
- 4Program on Smart and Sustainable Manufacturing, Academy of Innovative Semiconductor and Sustainable Manufacturing, National Cheng Kung University, Tainan701, Taiwan.
- 5Fire Protection and Safety Research Center, National Cheng Kung University, Tainan711, Taiwan.
Abstract
Coarse-grained (CG) water models with explicit electrostatic interactions are essential for capturing the dielectric response of biomolecular systems. However, many CG water models exhibit pronounced freezing artifacts under ambient conditions. Using the SPICA force-field family as a case study, we examined the influences of charge topology, background dielectric screening (εbg), and Lennard-Jones (LJ) interaction strength on the freezing tendency of CG water. Among the models examined here, which all reproduce the target density and surface tension at 298 K, changing the charge topology from a linear dipole to V-shaped geometries did not substantially lower the melting temperature (Tm). Instead, Tm was more sensitive to the LJ interaction strength. Reducing εbg strengthens the electrostatic contribution, allowing the LJ attraction to be weakened while maintaining the target properties. On this basis, we propose a parametrization framework for introducing a permanent dipole into a nonpolar CG water model while mitigating its freezing tendency. Using this framework, we developed pSPICA2, which reduces Tm from 384 to 302 K while retaining accurate density, surface tension, and dielectric response at 298 K. Although Tm remains above the experimental value of 273 K in this work, the proposed framework offers a practical route for developing polar CG water models with reduced freezing artifacts.
Related
- Freezing-induced biomineralization of calcium associated with amino-acid-like groups in a cold-tolerant cactus
- Accidental Origin of the Mpemba Effect in Freezing: Reversals in the Freezing Order Are Stochastic—Experimentation and Theory
- Reactivity-Guided Detection Strategies for Isocyanide Natural Products Reveal Acid-Induced Hazimycin Artifacts.
- Photocatalytic strategies for mitigating microplastic pollution and greenhouse gas emissions
- Polarization Is Required to Model Electrostatic Preorganization: The Case of Chorismate Mutase
BibTeX
@article{tsai2026mitigating,
title = {{Mitigating Freezing Artifacts in Polar Coarse-Grained Water through Rebalancing Electrostatic and van der Waals Interactions: A Case Study of the pSPICA Force Field.}},
author = {Yi-Chen Tsai and Qun-Yan Cheng and Wataru Shinoda and Chi-Cheng Chiu},
journal = {The journal of physical chemistry. B},
year = {2026},
month = oct,
doi = {10.1021/acs.jpcb.6c06050},
url = {https://europepmc.org/article/MED/42910870},
note = {CC BY},
}
RIS
TY - JOUR
TI - Mitigating Freezing Artifacts in Polar Coarse-Grained Water through Rebalancing Electrostatic and van der Waals Interactions: A Case Study of the pSPICA Force Field.
AU - Yi-Chen Tsai
AU - Qun-Yan Cheng
AU - Wataru Shinoda
AU - Chi-Cheng Chiu
PY - 2026
DA - 2026/10/01
JO - The journal of physical chemistry. B
DO - 10.1021/acs.jpcb.6c06050
UR - https://europepmc.org/article/MED/42910870
AB - Coarse-grained (CG) water models with explicit electrostatic interactions are essential for capturing the dielectric response of biomolecular systems. However, many CG water models exhibit pronounced freezing artifacts under ambient conditions. Using the SPICA force-field family as a case study, we examined the influences of charge topology, background dielectric screening (εbg), and Lennard-Jones (LJ) interaction strength on the freezing tendency of CG water. Among the models examined here, which all reproduce the target density and surface tension at 298 K, changing the charge topology from a linear dipole to V-shaped geometries did not substantially lower the melting temperature (Tm). Instead, Tm was more sensitive to the LJ interaction strength. Reducing εbg strengthens the electrostatic contribution, allowing the LJ attraction to be weakened while maintaining the target properties. On this basis, we propose a parametrization framework for introducing a permanent dipole into a nonpolar CG water model while mitigating its freezing tendency. Using this framework, we developed pSPICA2, which reduces Tm from 384 to 302 K while retaining accurate density, surface tension, and dielectric response at 298 K. Although Tm remains above the experimental value of 273 K in this work, the proposed framework offers a practical route for developing polar CG water models with reduced freezing artifacts.
SN - 1520-6106
SN - 1520-5207
ER -