{"id":"7d819355-b8ff-4866-8b7c-fbf917f5a904","arxiv_id":"2506.00153","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"TIP4P/epsilon performs as well as TIP4P/2005 on the Vega-Abascal water benchmark and, unlike TIP4P/2005, combines with unscaled Joung-Cheatham NaCl ions to give accurate solution densities.","lead":"Researchers simulated water and salt solutions with three common rigid-charge water models and found that TIP4P/epsilon matches the widely used TIP4P/2005 on a standard water benchmark while also reproducing the dielectric constant. This suggests that targeting the dielectric constant can yield transferable force fields for electrolytes without rescaling ionic charges.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The electrolyte conclusion does not isolate the dielectric constant: Table IV omits the JC water-model variant, so the NaCl density agreement may reflect cross-term provenance rather than ε.","rationale":"The reader's weakest assumption correctly identifies the unspecified JC water-model variant as the key gap in the electrolyte argument. I agree with that choice and sharpen it: the problem is not only a missing citation but a non-isolated causal inference. The paper's core empirical content—long TIP4P/ε simulations of pure water, dielectric constants, densities, transport properties, and ice polymorphs—appears carefully done, and the authors deserve credit for including the 'could be accidental' caveat and for reporting the poor solution viscosities. However, the headline claim that 'the accurate prediction of dielectric constants allows TIP4P/ε to describe densities of NaCl solutions' requires that ε be the variable that changes between the successful and unsuccessful cases. In the present comparison, both the water model and the ion-water cross terms change, so the experiment cannot distinguish the effect of ε from the effect of a different set of fitted cross interactions. Additionally, the assertion that the JC ions are fitted only to crystal and melt properties is inaccurate as stated, because the original JC parametrization also targeted solvation free energies. This does not require rejecting the paper: the empirical density curves are plausible, and the authors are appropriately tentative. But the causal claim should be treated as a hypothesis until the variant is documented and a controlled cross-term comparison is performed. The reader's CONDITIONAL verdict remains appropriate; my read does not move it to ACCEPT or REJECT, but it does reinforce the condition.","tokens_in":23439,"tokens_out":11113,"duration_ms":117634,"concrete_test":"Rerun the NaCl solution density and temperature-of-maximum-density calculations with exactly the same ion-ion and ion-water parameters from Table IV in both TIP4P/ε and TIP4P/2005, and also run the alternative JC parameter set intended for the reported water model once its variant is identified. If the TIP4P/ε advantage disappears or reverses when cross terms are held fixed, the dielectric constant is not the controlling factor; if the advantage persists with identical cross terms and with the JC variant explicitly stated, the claim would be supported. The authors should also state in the text which JC water-model variant Table IV corresponds to.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central causal claim—that reproducing the dielectric constant is what allows full-charge Joung-Cheatham ions to give acceptable NaCl solution densities without rescaling or reparametrization (Abstract; §V.A; Conclusions)—is not actually isolated by the presented simulations. Table IV lists NaCl and NaCl-water parameters labeled only as 'Joung-Cheatham NaCl parameters combined with TIP4P/ε water model,' without stating which water-model variant of the JC force field was used. This matters because JC parameters are host-water specific: the original 2008 parametrization optimized ion parameters for several water models and targeted both solvation free energies and crystal/melt properties. The manuscript's repeated characterization of the ions as 'targeted to their crystal and melt properties only' is therefore incomplete. If the cross terms in Table IV happen to be the JC/TIP4P set, those terms were fitted to a water model with a lower dielectric constant and to solvation data, not to TIP4P/ε. Moreover, Fig. 7(a) compares JC-TIP4P/ε against JC-TIP4P/2005 using a different set of NaCl-water cross interactions (Table II of Ref. 89), so any improvement in density could be due to the cross-term choice rather than to the dielectric constant of the host water. The authors' own parenthetical caveat that the better electrolyte performance 'could be, nevertheless, accidental' concedes this gap. The pure-water benchmark and the direct density data are not invalidated, but the causal inference that the dielectric constant is the operative variable is not supported without a controlled comparison.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript evaluates the TIP4P/epsilon water model against TIP4P/2005 and OPC using the Vega-Abascal benchmark for pure-water properties, and then tests its transferability to NaCl electrolyte solutions using the Joung-Cheatham (JC) ion force field. The authors report that TIP4P/epsilon matches TIP4P/2005 overall accuracy on the benchmark, largely outperforms OPC, and, when combined with full-charge JC ions, gives good solution densities and dielectric constants without charge rescaling or reparameterization, which they interpret as evidence that accurate dielectric constants are a valuable target property in water modeling.","tokens_in":81,"tokens_out":2586,"duration_ms":33966,"significance":"If the results hold, they challenge a common interpretation of the charge-scaling paradigm by showing that a rigid point-charge water model can simultaneously reproduce the dielectric constant and remain competitive with TIP4P/2005 on a wide range of properties, while also providing reasonable electrolyte thermodynamics with unmodified full-charge ions. The paper includes a substantial amount of original simulation data, particularly for dielectric constants, transport properties, melting point, ice polymorph densities, and electrolyte density/TMD, and the direct-coexistence melting point calculation for TIP4P/epsilon is a useful addition. The comparison with OPC and TIP4P/2005 is framed within a well-established benchmark, and the authors are candid about limitations, including the failure of solution viscosities.","major_comments":[{"comment":"The central causal claim--that reproducing the dielectric constant is what allows full-charge JC ions to give acceptable NaCl solution densities--is not isolated by the presented simulations because Table IV does not state which water-model variant of the JC force field was used for the Na+-O and Cl--O cross interactions. JC parameters are host-water specific, and the manuscript's characterization of the ions as 'targeted to their crystal and melt properties only' is incomplete if the cross terms originate from a different JC water parameter set. Since Fig. 7(a) compares JC-TIP4P/epsilon against JC-TIP4P/2005 using different water-ion cross interactions (Table IV vs. Table II of Ref. 89), the observed density improvement could arise from cross-term provenance rather than from the dielectric constant of the host water. The authors must specify the exact JC water-model variant used, justify its choice, and ideally run a control where the cross interactions are kept fixed while only the water model changes, or at least discuss why such a control is not conclusive.","section":"V.A, Table IV"},{"comment":"The conclusion that TIP4P/epsilon and TIP4P/2005 'perform equally well' rests on a comparison in which the TIP4P/2005 and OPC values are taken from Ref. 41 rather than recomputed with the same simulation protocols, and several TIP4P/epsilon values are taken from the original parametrization paper (marked with asterisks). This makes the reported benchmark scores (7.59 vs 7.54) appear more precise than the underlying data support. The authors should state the estimated statistical uncertainties in the VA scores and clarify whether the literature values were obtained with the same treatment of long-range corrections, system sizes, and ensemble conditions; otherwise the claim of equal performance is only as strong as the consistency of the compiled references.","section":"IV.F, Table III"},{"comment":"The paper reports that JC-TIP4P/epsilon yields viscosities that are far too large across the entire electrolyte concentration range (Fig. 9), yet the conclusions still refer to the model as 'a very robust model for preliminary studies of solution properties.' This is not internally contradictory, but the transport failure is a significant limitation that should be weighed more explicitly in the summary, particularly because the abstract highlights the advantage of dielectric constants for solution properties without mentioning that the dynamic properties are not transferable. The authors should temper the 'robust model' claim or explicitly state the domain of applicability (thermodynamic properties only) in the abstract and conclusions.","section":"V.E, Conclusions"}],"minor_comments":[{"comment":"The non-local screening model in Eqs. (4)-(6) is presented as a physical motivation for charge scaling, but it is not used in the simulations and is admittedly a 'mere caricature.' This section could be tightened to avoid giving the impression that the electrolyte results depend on this model; consider moving it to a more clearly labeled conceptual discussion or shortening it.","section":"II, Eq. (6)"},{"comment":"The comparison of dielectric constants for TIP4P/2005 and TIP4P/epsilon would benefit from error bars or shaded uncertainty bands on the simulation data, especially since the models are claimed to differ by about 25% and the reader cannot assess sampling uncertainty from the figure.","section":"IV.A, Fig. 1"},{"comment":"The Yeh-Hummer correction is applied to the diffusion coefficient, but the text does not state which viscosity value is used for the correction in each state point; please specify whether the simulated or experimental viscosity was used.","section":"III.D, Eq. (10)"},{"comment":"There are several typographical and formatting issues, including inconsistent rendering of 'TIP4P/epsilon' as 'TIP4P ε' in places, the use of 'gK' vs. 'g_K', and an apparent duplicate entry for Ref. 77; a careful proofread is needed.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and likely of interest to the water-model and force-field community. The main concern is not the soundness of the pure-water simulations but the strength of the causal inference in the electrolyte portion; the authors have already included a caveat, but the manuscript needs to address the missing JC water-variant information and ideally provide a control simulation to support the title's central question. Given that the issue is identifiable and fixable by further analysis and more careful phrasing, I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a careful, useful benchmark. The central pure-water result—TIP4P/ε scores essentially identically to TIP4P/2005 on the Vega-Abascal test while reproducing ε, with OPC clearly behind—holds up and is worth knowing. The electrolyte section is suggestive but does not isolate ε as the cause, and the abstract overclaims.\n\nWhat’s new: nobody had run TIP4P/ε through the Abascal-Vega battery against TIP4P/2005 and OPC, and the JC/TIP4P/ε density data are new. The methods are standard and clearly described; the melting point calculation with 10k molecules is honest work; the Kirkwood-factor analysis is a nice touch.\n\nSoft spots, in decreasing order:\n\n1. The electrolyte causal claim. Table IV never states which Joung-Cheatham variant the ion-water cross terms come from. JC parameters are host-water specific and were fitted to solvation free energies as well as crystal/melt properties, so the abstract’s “targeted to their crystal and melt properties only” is inaccurate. If these cross terms were optimized for TIP4P (ε≈53), their transfer to TIP4P/ε (ε≈79) could look fine for reasons unrelated to ε. The paper’s own “could be accidental” caveat concedes as much. That doesn’t kill the paper, but the abstract should be toned down.\n\n2. No error bars on most properties, and half the benchmark numbers are imported from Ref. 41 without recomputation. For a benchmark this is forgivable, but it means small differences between models are not significant.\n\n3. The theoretical section on non-local screening is suggestive but unused; fine as framing, but don’t lean on it.\n\nThe pure-water benchmark is the real contribution and it’s solid. The electrolyte part is a pilot result that needs a controlled comparison—e.g., same JC variant in TIP4P/2005 vs TIP4P/ε, or a JC set re-fit to TIP4P/ε—before the ε-causality story is credible.\n\nWho this is for: anyone choosing a rigid water model for large-scale simulations or working on force-field transferability. It deserves serious peer review; the authors should be asked to specify the JC variant, add error estimates, and soften the causal wording. Yes, engage.","headline":"Solid pure-water benchmark showing TIP4P/ε ties TIP4P/2005 while reproducing ε; the electrolyte causality claim is softer than the abstract suggests.","tokens_in":24330,"tokens_out":2906,"would_cite":true,"duration_ms":29719,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A water model tuned to the dielectric constant matches TIP4P/2005 across most properties and, without any charge rescaling, reproduces NaCl solution densities.","keywords":["TIP4P/ε","TIP4P/2005","dielectric constant","charge scaling","Joung–Cheatham NaCl","NaCl solution densities","Kirkwood factor","rigid point-charge water models"],"falsifier":"Re-fit the Na⁺–O and Cl⁻–O cross parameters self-consistently for TIP4P/ε from crystal and melt data and rerun the molality–density curves; if the densities shift by more than about 1%, the reported agreement was inherited from the original cross parameters. Alternatively, measure the ambient Kirkwood factor $g_K$ experimentally: if the true value is near 3.2 rather than 3.85, TIP4P/ε achieves its dielectric constant through incorrect orientational correlations.","tokens_in":23170,"feed_emoji":"💧","tokens_out":10092,"duration_ms":108295,"temperature":0.7,"pith_summary":"The paper asks whether targeting the static dielectric constant in a non-polarizable water force field necessarily sacrifices other properties, as the charge-scaling view of effective point charges would suggest. It reports that the TIP4P/ε model, whose charge geometry was tuned to reproduce ε, scores essentially identically to TIP4P/2005 on the standard benchmark for rigid water models (7.54 vs 7.59), while reproducing ε across wide temperature and pressure ranges where TIP4P/2005 is about 25% too low and clearly outscoring OPC (6.26). The same dielectric accuracy carries over to electrolytes: crystal-and-melt-tuned NaCl ions with full integer charges, inserted with no rescaling or mixing-rule adjustment, give solution densities and temperatures of maximum density close to experiment, comparable to results obtained with a scaled-charge electrolyte force field that was specifically parametrized for water. The payoff is that a simple rigid point-charge water model can be both accurate and transferable without invoking charge scaling, although electrolyte viscosities remain far too high. A sympathetic reading is that the dielectric constant is a legitimate and useful parametrization target rather than a symptom of a spoiled force field.","feed_headline":"Dielectric-tuned water model matches TIP4P/2005","feed_subtitle":"It also reproduces NaCl solution densities at full ion charges, with no parameter retuning.","key_machinery":"The load-bearing mechanism is the Kirkwood-factor decomposition of the static dielectric constant, $\\varepsilon = 1 + \\frac{4\\pi\\rho}{3 k_B T}\\mu^2 g_K$, which splits ε into a single-molecule dipole μ and a collective orientational correlation factor $g_K$. TIP4P/ε achieves its dielectric accuracy by placing the negative dummy site much closer to the oxygen than TIP4P/2005 does ($d_{\\mathrm{OM}} = 0.105$ Å versus 0.1546 Å), which raises the molecular dipole only about 5% (2.43 D versus 2.31 D) yet changes liquid-state orientational correlations enough to raise $g_K$ from 3.2 to 3.85 at ambient conditions. The paper pairs this identity with the standard benchmark used to rate rigid water models and with the unmodified Joung–Cheatham ion model to test transferability. The same identity also exposes the model's weak spot: in ice Ih at 273 K the Kirkwood factor collapses to 2.0, so TIP4P/ε still predicts the wrong trend of ε on freezing, like other point-charge models.","core_discovery":"The paper's central claim is that TIP4P/ε, a rigid four-point water model parametrized to reproduce the experimental dielectric constant, performs as well as TIP4P/2005 for most pure-water and ice properties and, because of its accurate dielectric response, transfers to electrolyte solutions without further fitting. Concretely, the two models earn nearly equal overall scores on the standard benchmark (7.59 for TIP4P/2005 versus 7.54 for TIP4P/ε, with OPC at 6.26); TIP4P/ε reproduces the static dielectric constant of liquid water over a wide temperature range at 1 bar, at saturation, and at 500 bar; and unmodified Joung–Cheatham NaCl at full integer charges yields solution densities and temperatures of maximum density close to experiment. The dielectric accuracy comes not from large charges (TIP4P/ε charges are only about 5% smaller than TIP4P/2005's) but from stronger orientational correlations, reflected in a Kirkwood factor $g_K \\approx 3.85$ versus 3.2. The authors conclude that a robust parametrization can be achieved without charge scaling, while acknowledging that the same point-charge limitations reappear in solution transport, with JC–TIP4P/ε viscosities far too large.","pith_inferences":["A testable extension of the paper's logic: if accurate ε is the causal factor, other rigid water models with accurate dielectric constants should also inherit the unmodified-JC-ion transferability; repeating the NaCl-density calculation with OPC water would probe this directly.","A conservative caveat: the JC ion–water cross parameters used here were fitted for a different host water model, and the paper does not state which variant; if that host had a low dielectric constant, the cross terms may already encode compensation, so the observed agreement would not isolate ε as the cause.","The Kirkwood-factor difference suggests a falsifiable physical picture: TIP4P/ε implies larger orientational correlations in ambient water than TIP4P/2005, and any experiment able to constrain $g_K$ could discriminate which model has the right collective structure.","A practical extension would be to map the dielectric decrement of JC–TIP4P/ε over the full molality range, since the reported solution densities are good at low concentration while the concentrated regime is exactly where the viscosity failure appears."],"forward_implications":["TIP4P/ε can be used as a rigid water model where accurate static dielectric response matters, matching TIP4P/2005 on densities, transport, ice polymorph densities, and surface tension while fixing the roughly 25% dielectric deficit.","Joung–Cheatham NaCl at full integer charges, with no charge rescaling and no Lorentz–Berthelot adjustment, gives acceptable solution densities and temperatures of maximum density in TIP4P/ε, at a level comparable to a specifically parametrized scaled-charge electrolyte force field.","For phase-boundary studies, TIP4P/2005 remains the better choice, since it predicts melting and critical points more accurately; TIP4P/ε melts at about 238 K.","Electrolyte transport remains a shared failure of both point-charge models: JC–TIP4P/ε viscosities are far too high, pointing to missing local charge fluctuations as the next bottleneck.","The benchmark scores imply that OPC, despite its good dielectric constant, is clearly worse overall, so dielectric accuracy alone does not guarantee a good water model."],"supporting_citations":[{"why":"Defines the TIP4P/ε model and its parametrization strategy, the central object under evaluation.","marker":"[39]"},{"why":"Defines TIP4P/2005, the principal comparator and the parent model from which TIP4P/ε is derived.","marker":"[66]"},{"why":"Supplies the benchmark rating scheme used to compare the three water models.","marker":"[16]"},{"why":"Provides the TIP4P/2005 and OPC benchmark values and the earlier finding that OPC underperforms TIP4P/2005.","marker":"[41]"},{"why":"Provides the Joung–Cheatham NaCl parameters used untuned in the electrolyte transferability test.","marker":"[30]"},{"why":"States the charge-scaling rationale that the paper argues against, predicting that dielectric targeting should spoil the force field.","marker":"[17,18]"},{"why":"Defines the OPC model, the third comparator with a dielectric-tuned but less successful parametrization.","marker":"[38]"},{"why":"Documents the scaled-charge electrolyte force-field strategy and the viscosity failure of rigid point-charge electrolytes used as the solution comparison.","marker":"[25]"}],"fun_headline_variants":["Dielectric-focused water model matches TIP4P/2005, plus salt mix","Water model tuned to dielectric constant excels in benchmarks and salt solutions","TIP4P/ε: dielectric accuracy yields salt solution densities without rescaling","Why dielectric constant matters: TIP4P/ε ties TIP4P/2005, wins for NaCl","Full-charge NaCl works with dielectric-matched water model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The electrolyte conclusion rests on the assumption that the Joung–Cheatham ion–water cross-interaction parameters used with TIP4P/ε remain valid for that model even though they were fitted for a different host water model; if those cross terms already compensate for the dielectric deficiency of their original host, the agreement with experiment would not isolate the dielectric constant as the cause.","fun_headline_variants_meta":{"raw":{"variants":["Dielectric-focused water model matches TIP4P/2005, plus salt mix","Water model tuned to dielectric constant excels in benchmarks and salt solutions","TIP4P/ε: dielectric accuracy yields salt solution densities without rescaling","Why dielectric constant matters: TIP4P/ε ties TIP4P/2005, wins for NaCl","Full-charge NaCl works with dielectric-matched water model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001115,"raw_usage":{"total_tokens":4760,"prompt_tokens":1177,"completion_tokens":3583,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":793,"completion_tokens_details":{"reasoning_tokens":3480}},"tokens_in":793,"tokens_out":3583,"duration_ms":30781,"temperature":1.0,"reasoning_tokens":3480,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:11:35.296554+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the Na⁺–O and Cl⁻–O cross parameters self-consistently for TIP4P/ε from crystal and melt data and rerun the molality–density curves; if the densities shift by more than about 1%, the reported agreement was inherited from the original cross parameters. Alternatively, measure the ambient Kirkwood factor $g_K$ experimentally: if the true value is near 3.2 rather than 3.85, TIP4P/ε achieves its dielectric constant through incorrect orientational correlations.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the benchmark rating scheme used to compare the three water models."},{"cited_title":"The motivation behind the choice of the JC model is three- fold","cited_arxiv_id":null,"evidence_quote":"Provides the Joung–Cheatham NaCl parameters used untuned in the electrolyte transferability test."},{"cited_title":"Blazquez et al","cited_arxiv_id":null,"evidence_quote":"Documents the scaled-charge electrolyte force-field strategy and the viscosity failure of rigid point-charge electrolytes used as the solution comparison."}],"review_version":1}