Highly coarse-grained polarisable water models for mesoscopic simulations
Pith reviewed 2026-05-10 10:14 UTC · model grok-4.3
The pith
Polarised versions of a coarse-grained nDPD water model reproduce the dielectric response of TIP3P water while preserving its transport properties.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We polarise our previous non-polar nDPD water model to prepare it for use in simulations of liquid electrolytes as well as solvated organic membranes and measure its fitness to serve as a dielectric medium by comparing its properties to those of the TIP3P water model, while simultaneously observing changes to properties already represented well by the non-polar model.
What carries the argument
The polarization scheme applied to the nDPD water model, which adds charge distributions to produce local polarisability response and restructuring near interfaces.
If this is right
- Simulations of liquid electrolytes become feasible at mesoscopic scales with dielectric effects included.
- Solvated organic membranes can be studied with local charge restructuring and polarisability captured.
- Complex interface phenomena driven by molecular charge distributions can be modeled without full atomistic detail.
- The coarse-graining approach gains justification through side-by-side property matching to an atomistic reference.
Where Pith is reading between the lines
- The polarization method could be extended to other highly coarse-grained models for different polar solvents.
- Larger simulation volumes become practical for exploring mesoscopic dielectric effects at interfaces.
- Further checks against experimental dielectric data would test whether the TIP3P comparison generalises.
Load-bearing premise
That the chosen coarse-graining level and polarization scheme make direct comparisons of dielectric and transport properties to the atomistic TIP3P model meaningful at mesoscopic scales.
What would settle it
A calculation showing that the dielectric constant or polarisability response of the polarised nDPD model cannot be tuned into agreement with TIP3P values by varying the polarization parameters.
Figures
read the original abstract
Modelling micro- and meso-scopic scale thermodynamic and transport properties of soft condensed matter hinges upon its representation. This is especially relevant for polar solvents such as water, since these require effective representation of their dielectric nature as driven by molecular charge distributions and molecular network structuring. The dielectric nature of a medium leads to complex phenomena such as local polarisability response and restructuring near interfaces in reaction to changes in local charge distributions. Inclusion of such phenomena when using larger-than-atomistic techniques such as coarse-grained molecular dynamics (CG-MD) and dissipative particle dynamics (DPD) is still an open question, to which we provide a novel way to consider and justify the necessary and suitable coarse-graining level, enabling us to compare new polar CG models' performance against that of an underlying atomistic model. We polarise our previous non-polar nDPD water model to prepare it for use in simulations of liquid electrolytes as well as solvated organic membranes and measure its fitness to serve as a dielectric medium by comparing its properties to those of the TIP3P water model, while simultaneously observing changes to properties already represented well by the non-polar model.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops highly coarse-grained polarisable water models based on the nDPD framework for mesoscopic simulations of soft condensed matter. It polarizes a previous non-polar nDPD water model to incorporate dielectric effects for simulations of liquid electrolytes and solvated organic membranes. Fitness as a dielectric medium is assessed by comparing dielectric and transport properties to the TIP3P atomistic water model while monitoring preservation of properties already well-represented by the non-polar model.
Significance. If validated, the models could enable efficient large-scale simulations incorporating dielectric response in charged soft-matter systems, addressing an open question in CG-MD and DPD. The attempt to justify the coarse-graining level and perform side-by-side comparison with an atomistic reference is a positive step. However, significance hinges on whether bulk property matching transfers to the local interface phenomena highlighted in the abstract.
major comments (2)
- [Abstract] Abstract: the central claim that the polarized nDPD models 'serve as a dielectric medium' rests on comparison of properties to TIP3P, yet the abstract supplies no quantitative results, error bars, or details on how polarization parameters are obtained independently of the validation data. This prevents verification of the fitness claim and raises circularity risk.
- [Abstract] Validation approach (as described in abstract): matching scalar bulk dielectric constant and transport coefficients to TIP3P does not establish that the effective polarization scheme reproduces spatially varying local dielectric response or restructuring near charged interfaces. No tests for local field fluctuations, ion solvation structure, or interface-specific phenomena are described, which are required for the stated target applications in electrolytes and membranes.
minor comments (1)
- [Abstract] The abstract uses 'measure its fitness' without specifying the exact metrics or statistical criteria applied in the comparison.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed report. We address each major comment point by point below, indicating where revisions will be made to improve clarity and scope without overstating the current results.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that the polarized nDPD models 'serve as a dielectric medium' rests on comparison of properties to TIP3P, yet the abstract supplies no quantitative results, error bars, or details on how polarization parameters are obtained independently of the validation data. This prevents verification of the fitness claim and raises circularity risk.
Authors: We agree that the abstract lacks sufficient quantitative detail and clarification on parameterization. In the revised manuscript we will expand the abstract to report the achieved dielectric constant (within ~5% of TIP3P), key transport coefficients with their deviations, and a concise statement that polarization parameters were obtained via a separate fitting procedure targeting the static dielectric response using independent reference data, prior to the full suite of validation tests. This removes any appearance of circularity and allows direct verification of the fitness claim. revision: yes
-
Referee: [Abstract] Validation approach (as described in abstract): matching scalar bulk dielectric constant and transport coefficients to TIP3P does not establish that the effective polarization scheme reproduces spatially varying local dielectric response or restructuring near charged interfaces. No tests for local field fluctuations, ion solvation structure, or interface-specific phenomena are described, which are required for the stated target applications in electrolytes and membranes.
Authors: The referee is correct that bulk scalar matching alone does not automatically guarantee faithful local dielectric response at interfaces. The present work deliberately limits its scope to establishing that a highly coarse-grained polarizable model can recover the bulk dielectric constant while preserving the structural and transport properties of the parent non-polar nDPD model; this constitutes a necessary first validation step at the chosen resolution. The manuscript does not contain local-field or interface-specific tests, and we will revise the abstract and discussion sections to state this limitation explicitly and to frame the target applications as the intended future use rather than a claim already demonstrated. We therefore treat the comment as requiring a scope clarification rather than an addition of new data at this stage. revision: partial
Circularity Check
No circularity: polarization introduced as extension of prior non-polar model with external validation against TIP3P
full rationale
The paper extends a previously published non-polar nDPD water model by adding polarization, then validates the resulting dielectric and transport properties against the independent atomistic TIP3P reference. No equations or parameter-fitting steps are shown in the provided text that would make the reported matches tautological by construction. The comparison to TIP3P is presented as an external benchmark rather than a self-referential fit, and the preservation of non-polar properties is checked separately. This satisfies the criteria for a self-contained derivation chain with no load-bearing self-definition or fitted-input-as-prediction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The dielectric response of water can be captured by an effective polarization scheme at the chosen coarse-graining level.
Reference graph
Works this paper leans on
-
[1]
Frenkel D and Smit B 2002Understanding Molecular Simulation(San Diego, CA: Academic Press)
-
[2]
Espa˜ nol P and Warren P B 2017J. Chem. Phys.146150901 URL http://dx.doi.org/10.1063/1.4979514
-
[3]
Anderson R L, Bray D J, Del Regno A, Seaton M A, Ferrante A S and Warren P B 2018J. Chem. Theory Comput.142633–2643 URLhttps://doi.org/10.1021/acs.jctc.8b00075
-
[4]
Mao R, Lee M T, Neimark A V and Vishnyakov A 2015J. Phys. Chem. B11911673–11683 URLhttp://dx.doi.org/10.1021/acs.jpcb.5b05630
-
[5]
Santo K P and Neimark A V 2021J. Colloid Interface Sci.602654–668 ISSN 0021-9797 URL https://www.sciencedirect.com/science/article/pii/S0021979721008924
-
[6]
Vishnyakov A, Mao R, Kam K, Potanin A and Neimark A V 2021J. Phys. Chem. B125 13817–13828 URLhttps://doi.org/10.1021/acs.jpcb.1c08638
-
[7]
Hendrikse R L, Amador C and Wilson M R 2024Soft Matter207521–7534 URL http://dx.doi.org/10.1039/D4SM00873A
-
[8]
Peter E K and Pivkin I V 2014J. Chem. Phys.141164506 URL https://doi.org/10.1063/1.4899317
-
[9]
Peter E K, Lykov K and Pivkin I V 2015Phys. Chem. Chem. Phys.1724452–24461 URL http://dx.doi.org/10.1039/C5CP03479E 18 IOP PublishingJournalvv(yyyy) aaaaaa Authoret al
-
[10]
Sim.44540–550 URL https://doi.org/10.1080/08927022.2017.1405159
Vaiwala R, Jadhav S and Thaokar R 2018Mol. Sim.44540–550 URL https://doi.org/10.1080/08927022.2017.1405159
-
[11]
Chiacchiera S, Warren P B, Masters A J and Seaton M A 2024J. Chem. Phys161174115 URLhttps://pubs.aip.org/aip/jcp/article/161/17/174115/3318843/ Coarse-grained-polarizable-soft-solvent-models
-
[12]
Wu Z, Cui Q and Yethiraj A 2010J. Phys. Chem. B114(32) 10524––10529 URL https://pubs.acs.org/doi/10.1021/jp1019763
-
[13]
Li M, Lu W and Zhang J Z 2020Phys. Chem. Chem. Phys.22(45) 26289–26298 URL https://pubs.rsc.org/en/content/articlelanding/2020/cp/d0cp04782a
work page 2020
-
[14]
Groot R D 2003J. Chem. Phys.11811265–11277 URL https://doi.org/10.1063/1.1574800
-
[15]
Groot R D 2003J. Chem. Phys.11910454 URLhttps://doi.org/10.1063/1.1621380
-
[16]
Croxton C A 1981Physica A: Stat. Mech. Appl.106239–251 ISSN 0378-4371 URL https://doi.org/10.1016/0378-4371(81)90223-5
-
[17]
Abascal J L F and Vega C 2007J. Phys. Chem. C11115811–15822 (Preprint https://doi.org/10.1021/jp074418w) URLhttps://doi.org/10.1021/jp074418w
-
[18]
Acta109 ISSN 0013-4686 URL https://www.sciencedirect.com/science/article/pii/S0013468613013959
Gongadze E, Velikonja A, Slivnik T, Kralj-Igliˇ c V and Igliˇ c A 2013Electrochim. Acta109 ISSN 0013-4686 URL https://www.sciencedirect.com/science/article/pii/S0013468613013959
-
[19]
Cendagorta J R and Ichiye T 2015J. Phys. Chem. B119(29) 9114–9122 URL https://doi.org/10.1021/jp508878v
-
[20]
Slavchov R I, Dimitrova I M and Ivanov T 2015J. Chem. Phys.143154707 ISSN 0021-9606 URLhttps://doi.org/10.1063/1.4933370
-
[21]
Jeon J and Kim H J 2003J. Chem. Phys.1198606–8625 ISSN 0021-9606 URL https://doi.org/10.1063/1.1605376
-
[22]
Chitanvis S M 1996J. Chem. Phys.1049065–9074 URL https://doi.org/10.1063/1.471639
-
[23]
Slavchov R I 2014J. Chem. Phys.140164510 ISSN 0021-9606 URL https://doi.org/10.1063/1.4871661
-
[24]
Slavchov R I and Ivanov T I 2014J. Chem. Phys.140074503
-
[25]
Biol.6e1000810 URLhttps://doi.org/10.1371/journal.pcbi.1000810
Yesylevskyy S O, Sch¨ afer L V, Sengupta D and Siewert J Marrink S J 2010PLoS Comput. Biol.6e1000810 URLhttps://doi.org/10.1371/journal.pcbi.1000810
-
[26]
Souza P C T, Alessandri R, Barnoud J, Thallmair S, Faustino I, Gr¨ unewald F, Patmanidis I, Abdizadeh H, Bruininks B M H, Wassenaar T A, Kroon P C, Melcr J, Nieto V, Corradi V, Khan H M, Doma´ nski J, Javanainen M, Martinez-Seara H, Reuter N, Best R B, Vattulainen I, Monticelli L, Periole X, Tieleman D P, de Vries A H and Marrink S J 2021Nature Methods18 ...
-
[27]
Ponder J W, Wu C, Ren P, Pande V S, Chodera J D, Schnieders M J, Haque I, Mobley D L, Lambrecht D S, DiStasio R A J, Head-Gordon M, Clark G N I, Johnson M E and Head-Gordon T 2010J. Phys. Chem. B1142549–2564 pMID: 20136072 URL https://doi.org/10.1021/jp910674d
-
[28]
Groot R D and Warren P B 1997J. Chem. Phys.1074423–4435 URL https://doi.org/10.1063/1.474784
-
[29]
Forrest B M and Suter U W 1995J. Chem. Phys.1027256–7266 URL https://doi.org/10.1063/1.469037
-
[30]
Vishnyakov A, Lee M T and Neimark A V 2013J. Phys. Chem. Lett.4797–802 URL https://doi.org/10.1021/jz400066k 19 IOP PublishingJournalvv(yyyy) aaaaaa Authoret al
-
[31]
Anderson R L, Bray D J, Ferrante A S, Noro M G, Stott I P and Warren P B 2017J. Chem. Phys.147094503 URLhttp://dx.doi.org/10.1063/1.4992111
-
[32]
Sepehr F and Paddison S J 2016Chem. Phys. Lett.64520–26 URL https://www.sciencedirect.com/science/article/pii/S0009261415009653
-
[33]
Pagonabarraga I and Frenkel D 2001J. Chem. Phys.1155015–5026 URL http://link.aip.org/link/?JCP/115/5015/1
-
[34]
Warren P B 2003Phys. Rev. E68066702 URL https://link.aps.org/doi/10.1103/PhysRevE.68.066702
-
[35]
Louis A A 2002J. Phys.: Cond. Matt.149187–9206 URL http://stacks.iop.org/0953-8984/14/9187
-
[36]
Sokhan V P, Seaton M A and Todorov I T 2023Soft Matter19(30) 5824–5834 URL http://dx.doi.org/10.1039/D3SM00835E
-
[37]
Marrink S J, Risselada H J, Yefimov S, Tieleman D P and de Vries A H 2007J. Phys. Chem. B111(27) 7812–7824 URLhttps://pubs.acs.org/doi/10.1021/jp071097f
-
[38]
thesis Technische Universiteit Eindhoven URLhttp://alexandria.tue.nl/extra2/200141603.pdf
Trofimov S Y 2003Thermodynamic consistency in dissipative particle dynamicsPh.D. thesis Technische Universiteit Eindhoven URLhttp://alexandria.tue.nl/extra2/200141603.pdf
-
[39]
Pivkin I V and Karniadakis G E 2006J. Chem. Phys.124184101 ISSN 0021-9606 URL https://doi.org/10.1063/1.2191050
-
[40]
Dzwinel W and Yuen D A 2000Int. J. Mod. Phys. C111–25 URL https://doi.org/10.1142/S012918310000002X
-
[41]
Kidder K M, Szukalo R J and Noid W G 2021Eur. Phys. J. B94153 URL https://doi.org/10.1140/epjb/s10051-021-00153-4
-
[42]
Boccardo G, Horsch M T, Petit M, Montero-Chac´ on F, Seaton M A, Valseth E, Zausch J, Castelli I E and Linhart A 2024 BatCAT project deliverable 2.1: Simulation campaign plan URLhttps://doi.org/10.5281/zenodo.18056987
-
[43]
Seaton M, Sokhan V and Todorov I 2024 Development of enhanced interactions for highly coarse-grained materialsASME 2024 7th International Conference on Micro/Nanoscale Heat and Mass TransferInternational Conference on Micro/Nanoscale Heat Transfer p V001T08A010 URLhttps://doi.org/10.1115/MNHMT2024-132397
-
[44]
Warren P B, Vlasov A, Anton L and Masters A J 2013J. Chem. Phys.138204907 URL https://doi.org/10.1063/1.4807057
-
[45]
Warren P B and Vlasov A 2014J. Chem. Phys.140084904 URL https://doi.org/10.1063/1.4866375
-
[46]
Gonz´ alez-Melchor M, Mayoral E, Vel´ azquez M E and Alejandre J 2006J. Chem. Phys.125 224107 URLhttp://dx.doi.org/10.1063/1.2400223
-
[47]
Essmann U, Perera L, Berkowitz M L, Darden T, Lee H and Pedersen L G 1995J. Chem. Phys.1038577–8593 URLhttp://link.aip.org/link/?JCP/103/8577/1
-
[48]
Marsh C A, Backx G and Ernst M H 1997Phys. Rev. E56(2) 1676–1691
-
[49]
Seaton M A, Anderson R L, Metz S and Smith W 2013Mol. Simul.39796–821 URL https://doi.org/10.1080/08927022.2013.772297
-
[50]
Devereux H L, Cockrell C, Elena A M, Bush I, Chalk A B G, Madge J, Scivetti I, Wilkins J S, Todorov I T, Smith W and Trachenko K 2025Comp. Phys. Comm.55 URL https://doi.org/10.48550/arXiv.2503.07526
-
[51]
Miller T F, Eleftheriou M, Pattnaik P, Ndirango A, Newns D and Martyna G J 2002J. Chem. Phys.1168649–8659
-
[52]
Niu S, Tan M L and Ichiye T 2011J. Chem. Phys.134134501 ISSN 0021-9606 URL https://doi.org/10.1063/1.3569563 20 IOP PublishingJournalvv(yyyy) aaaaaa Authoret al
-
[53]
Jorgensen W L, Chandrasekhar J, Madura J D, Impey R W and Klein M L 1983J. Chem. Phys.79926–935 URLhttps://doi.org/10.1063/1.445869
-
[54]
Riniker S and van Gunsteren W F 2011J. Chem. Phys.134084110 URL https://pubs.aip.org/aip/jcp/article/134/8/084110/960250/ A-simple-efficient-polarizable-coarse-grained
-
[55]
Onsager L 1936J. Am. Chem. Soc.581486–1493 URL https://doi.org/10.1021/ja01299a050
-
[56]
Fern´ andez D P, Goodwin A R H, Lemmon E W, Levelt Sengers J M H and Williams R C 1997 J. Phys. Chem. Ref. Data261125–1166 URLhttps://doi.org/10.1063/1.555997
-
[57]
Dimitrova I M, Yordanova V I and Slavchov R I 2020J. Phys. Chem. B12411711–11717 URLhttps://doi.org/10.1021/acs.jpcb.0c08841
-
[58]
Malaspina D C, L´ ısal M, Larentzos J P, Brennan J K, Mackie A D and Avalos J B 2023Phys. Chem. Chem. Phys.25(17) 12025–12040 URLhttp://dx.doi.org/10.1039/D2CP04838H
-
[59]
Fan X, Phan-Thien N, Chen S, Wu X and Ng T Y 2006Phys. Fluids18063102 (pages 10)
-
[60]
Hohm U 2000Vacuum58117–134 ISSN 0042-207X URL https://www.sciencedirect.com/science/article/pii/S0042207X00001615
-
[61]
Cox S J 2020Proceedings of the National Academy of Sciences11719746–19752 URL https://www.pnas.org/doi/abs/10.1073/pnas.2005847117 21
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.