SPC/Fw Water Model for Charged Interfaces

• SPC/Fw is a flexible three-site water model used in molecular dynamics to simulate aqueous solid–liquid interfaces and electrolyte solutions with realistic dielectric behavior.
• It employs explicit intramolecular flexibility and Lennard-Jones interactions on oxygen to implicitly account for polarization and achieve consistency with continuum electrostatics.
• The model reliably reproduces salt-dependent dielectric constants up to 0.5 M NaCl, though it shows systematic deviations and ion-pairing effects at higher concentrations.

SPC/Fw is a flexible three-site water model used in molecular dynamics simulations of aqueous solid–liquid interfaces and electrolyte solutions. In the comparative study of a charged silica surface in NaCl solution, SPC/Fw was examined both as a bulk solvent model and as an interfacial model, with particular emphasis on its dielectric response, its representation of Stern-layer charge, and its consistency with continuum electrostatics. Within that study, SPC/Fw emerged as one of the more reliable three-site models for reproducing the salt-dependent dielectric constant up to moderately saline conditions, while showing systematic deviations and ion-pairing effects at higher NaCl concentrations (Tavakol et al., 4 Aug 2025).

1. Model definition and force-field structure

In the reported parametrization, SPC/Fw is represented by a flexible three-site model consisting of one oxygen atom and two hydrogen atoms, with no extra “M” charge site. Lennard-Jones interactions are placed on oxygen only, and no explicit Drude or induced-dipole polarization is used; instead, water polarization is described implicitly through charge redistribution in the force field (Tavakol et al., 4 Aug 2025).

The SPC/Fw parameters reported for the charged-interface comparison are as follows.

Parameter	Value
$\sigma_O$	$3.16549\,\text{\AA}$
$\epsilon_{OO}$	$0.155\,\text{kcal mol}^{-1}$
$q_H$	$+0.41\,e$
$q_O$	$-0.82\,e$
$k_{OH}$	$529.5\,\text{kcal mol}^{-1}\,\text{\AA}^{-2}$
$3.16549\,\text{\AA}$0	$3.16549\,\text{\AA}$1
$3.16549\,\text{\AA}$2	$3.16549\,\text{\AA}$3
$3.16549\,\text{\AA}$4	$3.16549\,\text{\AA}$5

The explicit intramolecular flexibility distinguishes SPC/Fw from rigid three-site models in the same comparison set. In the study, SPC/Fw was evaluated alongside SPC/e, TIPS3p, H2O/DC, TIP3P-Fw, OPC3, TIP3P, TIP3P-FB, TIP3P-ST, FBA/e, and TIPS3p-PPPM. The comparative framing is significant because the paper’s conclusions depend not on isolated bulk-water properties alone, but on how accurately a model transports those properties into charged interfacial simulations (Tavakol et al., 4 Aug 2025).

2. Bulk dielectric response in NaCl solutions

The principal bulk observable used to assess SPC/Fw was the static dielectric constant $3.16549\,\text{\AA}$6 as a function of NaCl concentration. Bulk simulations were performed in a $3.16549\,\text{\AA}$7 box containing 3000 SPC/Fw water molecules at NaCl concentrations of $3.16549\,\text{\AA}$8, $3.16549\,\text{\AA}$9, $\epsilon_{OO}$0, $\epsilon_{OO}$1, $\epsilon_{OO}$2, and $\epsilon_{OO}$3. Long-range electrostatics were treated with the PPPM algorithm, bonds and angles were integrated with a $\epsilon_{OO}$4 timestep, $\epsilon_{OO}$5 of equilibration preceded $\epsilon_{OO}$6 of production, and the dielectric constant was computed from dipole fluctuations (Tavakol et al., 4 Aug 2025).

The reported fluctuation formula is

$\epsilon_{OO}$7

with

$\epsilon_{OO}$8

The SPC/Fw results were compared directly to experimental values from Buchner et al. The reported values are:

NaCl (M)	$\epsilon_{OO}$9	$0.155\,\text{kcal mol}^{-1}$0
0.00	$0.155\,\text{kcal mol}^{-1}$1	78
0.50	$0.155\,\text{kcal mol}^{-1}$2	71
1.00	$0.155\,\text{kcal mol}^{-1}$3	64
1.50	$0.155\,\text{kcal mol}^{-1}$4	57
2.00	$0.155\,\text{kcal mol}^{-1}$5	52
2.50	$0.155\,\text{kcal mol}^{-1}$6	44

These data show that SPC/Fw reproduces $0.155\,\text{kcal mol}^{-1}$7 within approximately $0.155\,\text{kcal mol}^{-1}$8–$0.155\,\text{kcal mol}^{-1}$9 units, or at most $q_H$0, up to $q_H$1. At higher salt concentrations it overestimates $q_H$2 by up to $q_H$3 (Tavakol et al., 4 Aug 2025). For interfacial electrostatics, that pattern is important because the screening length, ion free-energy landscape, and Poisson–Boltzmann comparison all depend on the solvent dielectric response.

3. Charged silica interface simulations

The interfacial benchmark for SPC/Fw was a charged silica slab with a Q$q_H$4 surface at $q_H$5, carrying $q_H$6 groups per $q_H$7. The slab thickness was $q_H$8, the cross-section was $q_H$9, and the fluid chamber length was $+0.41\,e$0. Simulations used periodic boundary conditions in all directions, PPPM for long-range electrostatics, a $+0.41\,e$1 switch cutoff for short-range interactions with CHARMM-style smoothing, and a Nosé–Hoover ensemble at $+0.41\,e$2 and $+0.41\,e$3, followed by $+0.41\,e$4 of production after equilibration (Tavakol et al., 4 Aug 2025).

A central construct in the analysis is the Stern layer. In the SPC/Fw simulations, the Stern layer was defined operationally as the region from the solid surface to the position $+0.41\,e$5 of the free-energy minimum of the $+0.41\,e$6 ion, approximately $+0.41\,e$7 from the surface. The total ionic charge in this layer was computed by integrating the ionic number-density difference,

$+0.41\,e$8

For SPC/Fw, the resulting Stern-layer ionic charge was $+0.41\,e$9 per face, while the surface charge was $q_O$0 per face (Tavakol et al., 4 Aug 2025). This partial compensation of the surface charge is essential to the continuum comparison, because the analytic electrostatic potential was not parameterized from the nominal silica charge alone; it used the effective surface charge given by the silica plus the MD-determined Stern-layer charge.

4. Free-energy profiles and continuum correspondence

The one-dimensional ion free energy extracted from SPC/Fw molecular dynamics was defined as

$q_O$1

where $q_O$2 is the local ion density at position $q_O$3 and $q_O$4 is the bulk concentration (Tavakol et al., 4 Aug 2025).

For comparison with continuum electrostatics, the study used a Gouy–Chapman formulation with

$q_O$5

and

$q_O$6

Here $q_O$7 is set by the effective surface charge, namely the silica charge together with the Stern-layer ionic charge (Tavakol et al., 4 Aug 2025).

For SPC/Fw, the reported free-energy minima at the counter-ion peak were:

$q_O$8 (M)	$q_O$9
0.042	$-0.82\,e$0
0.21	$-0.82\,e$1
0.42	$-0.82\,e$2
0.84	$-0.82\,e$3

Using the experimental dielectric constant in the analytic expression produced $-0.82\,e$4 values within $-0.82\,e$5 of the SPC/Fw molecular-dynamics minima at all reported concentrations. If instead one used $-0.82\,e$6 in the analytic model, the predicted well became slightly too deep, by up to $-0.82\,e$7 at $-0.82\,e$8 (Tavakol et al., 4 Aug 2025).

These results show that the continuum comparison is not controlled by dielectric constant alone. Agreement depended on including the MD-determined total charge of ions in the Stern layer and an appropriate dielectric constant. A plausible implication is that SPC/Fw can be used consistently in continuum-coupled interfacial analysis provided the Stern-layer charge is measured directly in MD rather than inferred from nominal surface chemistry alone.

5. Regime of validity and high-salinity limitations

The SPC/Fw results support a clear concentration-dependent interpretation. Up to approximately $-0.82\,e$9 NaCl, SPC/Fw reproduced the dielectric constant with relatively small error and therefore behaved as a reliable three-site solvent model for moderately saline charged interfaces. The paper explicitly characterizes SPC/Fw as among the most reliable three-site models for reproducing $k_{OH}$0 up to about $k_{OH}$1, and therefore as suitable for moderately saline interfaces (Tavakol et al., 4 Aug 2025).

At higher concentrations, however, two limitations became prominent. First, SPC/Fw overestimated the dielectric constant by as much as $k_{OH}$2. Second, for salt concentrations above $k_{OH}$3 NaCl, the simulations exhibited random ion–ion pair formation, which limited reproducibility and the applicability of the analytical method. At $k_{OH}$4, the reported number of ion pairs was approximately $k_{OH}$5 in the simulation box (Tavakol et al., 4 Aug 2025). Because these correlations are not captured by Poisson–Boltzmann theory, quantitative agreement deteriorated in the diffuse tail and in the condensation parameter beyond the Debye length.

The study therefore distinguishes two different claims. One is that SPC/Fw remains useful for direct MD at high concentration, because the simulation naturally captures ion pairing and associated correlations. The other is that continuum reduction becomes less reliable in that regime unless explicit corrections are introduced (Tavakol et al., 4 Aug 2025).

6. Position within the broader literature and nomenclature

Within the supplied arXiv literature, “SPC” is an overloaded acronym. It refers to the Spatial Proximity and Connectivity method in segregation measurement (Roberto, 2015), the stationary point concentration technique in FMCW radar (Park et al., 2020), single parity check coding in concatenated FEC (Lentner et al., 2022), smoothed phase-coded FMCW waveforms (Kumbul et al., 2022), and single-photon-counting detector modes for XFEL imaging (Wang, 2015). In molecular simulation, by contrast, SPC/Fw denotes the flexible water model summarized above, and its meaning is specific to atomistic simulations of water and electrolyte interfaces (Tavakol et al., 4 Aug 2025).

In the charged silica–electrolyte study, SPC/Fw occupied a pragmatic middle ground. It offered good dielectric accuracy, computational efficiency, and transferability for electrified or charged solid–liquid interfaces in the $k_{OH}$6–$k_{OH}$7 NaCl range, while showing manageable but systematic deviations at higher salinity (Tavakol et al., 4 Aug 2025). This suggests a specific methodological role for SPC/Fw: it is well suited to molecular-dynamics studies that require realistic dielectric response and explicit Stern-layer structure, especially when continuum comparisons are calibrated by MD-derived interfacial charge rather than by nominal surface parameters alone.