Cation Dominated but Negatively Charged Na2SO4,aq-Graphene Interfaces
Pith reviewed 2026-05-16 15:24 UTC · model grok-4.3
The pith
Na+ ions accumulate at the Na2SO4-graphene interface while the near-surface region carries a net negative charge from ion ratios
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Our results show that Na+ ions accumulate between the outermost and second water layers whereas SO4^{2-} ions accumulate within the second interfacial water layer indicating cation dominated interfaces. We find that the interfacial region (within ~10 Å of the graphene sheet) is negatively charged due to sub-stoichiometric Na+/SO4^{2-} ratio at the interface. Our simulated SFG spectra show enhancement and a red-shift of the spectra in the hydrogen bonded region as a function of Na2SO4 concentration similar to measurements due to SO4^{2-}-induced changes in the orientational order of water molecules in the second interfacial layer.
What carries the argument
Machine learning interatomic potential simulations paired with computed sum frequency generation spectra that track ion positions and water orientations layer by layer
Load-bearing premise
The machine learning interatomic potential reproduces the correct ion layering and water orientations at the interface even though it was not fitted to experimental data for this particular electrolyte-graphene system.
What would settle it
High-resolution X-ray reflectivity or direct surface force measurements that quantify the net charge density within the first 10 Å of the graphene sheet would confirm or contradict the reported sub-stoichiometric Na+/SO4^{2-} ratio.
read the original abstract
The distribution of ions and their impact on the structure of electrolyte interfaces plays an important role in many applications. Interestingly, recent experimental studies have suggested the preferential accumulation of $SO_4^{2-}$ ions at the $Na_2SO_{4,aq}$-graphene interface in disagreement with the generally known tendency of cations to accumulate at graphene-electrolyte interfaces. Herein, we resolve the atomistic structure of the $Na_2SO_{4,aq}$-graphene interfaces in the 0.1-2.0 M concentration range using machine learning interatomic potential-based simulations and simulated sum frequency generation (SFG) spectra to reveal the molecular origins of the conundrum. Our results show that Na+ ions accumulate between the outermost and second water layers whereas $SO_4^{2-}$ ions accumulate within the second interfacial water layer indicating cation dominated interfaces. We find that the interfacial region (within ~10 ${\AA}$ of the graphene sheet) is negatively charged due to sub-stoichiometric $Na^+$/$SO_4^{2-}$ ratio at the interface. Our simulated SFG spectra show enhancement and a red-shift of the spectra in the hydrogen bonded region as a function of $Na_2SO_4$ concentration similar to measurements due to $SO_4^{2-}$-induced changes in the orientational order of water molecules in the second interfacial layer. Our study demonstrates that ion stratification and ion-induced water reorganization are key elements of understanding the electrolyte-graphene interface.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses machine learning interatomic potential molecular dynamics simulations combined with calculated sum frequency generation (SFG) spectra to examine the Na₂SO₄(aq)-graphene interface in the 0.1-2.0 M range. The central finding is that Na⁺ ions accumulate between the outermost and second water layers, SO₄²⁻ ions within the second layer, resulting in cation-dominated interfaces that are net negatively charged within ~10 Å of the graphene due to a sub-stoichiometric Na⁺/SO₄²⁻ ratio. The simulated SFG spectra reproduce the experimental concentration-dependent red-shift and intensity increase in the hydrogen-bonded region, attributed to SO₄²⁻-induced reorientation of water molecules.
Significance. If the ion distributions are accurate, the work provides a molecular-level explanation for the structure of electrolyte-graphene interfaces, reconciling SFG experiments with ion stratification that differs from simple expectations. The use of ML potentials allows access to larger scales while the SFG calculation links simulation to experiment. This is significant for understanding double-layer structure in applications like energy storage and sensing.
major comments (2)
- [Methods (ML potential training and validation)] The fidelity of the machine learning interatomic potential for Na⁺ and SO₄²⁻ distributions at the graphene-water interface is not benchmarked against ab initio calculations or direct experimental probes such as X-ray reflectivity. The SFG spectral agreement validates water orientational order but does not independently confirm the absolute ion number densities or the sub-stoichiometric ratio central to the negative charge claim (see Results on density profiles).
- [Results (ion density profiles and charge distribution)] The claim of negative charge in the ~10 Å interfacial region relies on the integrated ion densities showing sub-stoichiometry. However, no quantitative error analysis or convergence tests with respect to simulation length or system size are provided to establish the statistical significance of the reported Na⁺/SO₄²⁻ ratio deviation.
minor comments (2)
- [Abstract] The phrasing 'cation dominated interfaces' in the abstract could be clarified to avoid confusion with the title's 'Cation Dominated but Negatively Charged'.
- [Figures] Ensure that all concentration-dependent plots include clear legends and error bars where applicable for ion densities.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable comments. Below, we address each major comment in detail and indicate the revisions made to the manuscript.
read point-by-point responses
-
Referee: [Methods (ML potential training and validation)] The fidelity of the machine learning interatomic potential for Na⁺ and SO₄²⁻ distributions at the graphene-water interface is not benchmarked against ab initio calculations or direct experimental probes such as X-ray reflectivity. The SFG spectral agreement validates water orientational order but does not independently confirm the absolute ion number densities or the sub-stoichiometric ratio central to the negative charge claim (see Results on density profiles).
Authors: We agree that additional benchmarking of the MLIP for ion distributions at the specific interface would be beneficial. The potential was trained on ab initio data from bulk solutions and various interfaces, and we have validated it against experimental properties such as radial distribution functions and diffusion coefficients in bulk. However, we acknowledge the lack of direct comparison to AIMD for the graphene interface or X-ray data. In the revised manuscript, we will include a more detailed description of the training set and validation for ion-related properties. The SFG agreement supports the water structure influenced by ions, and the ion distributions are consistent with known behaviors in similar systems. We will add a note on this limitation and suggest future work for direct benchmarking. revision: partial
-
Referee: [Results (ion density profiles and charge distribution)] The claim of negative charge in the ~10 Å interfacial region relies on the integrated ion densities showing sub-stoichiometry. However, no quantitative error analysis or convergence tests with respect to simulation length or system size are provided to establish the statistical significance of the reported Na⁺/SO₄²⁻ ratio deviation.
Authors: We thank the referee for pointing this out. In the original manuscript, we focused on the average profiles but did not include detailed error analysis. We have now performed additional analysis using block averaging over the simulation trajectory to estimate uncertainties in the ion densities and the integrated Na⁺/SO₄²⁻ ratio. Convergence tests with respect to simulation length (extending to 50 ns) and system size (doubling the lateral dimensions) confirm that the sub-stoichiometric ratio is robust and statistically significant, with the deviation exceeding the estimated error bars. We will include these results, error bars on the density profiles, and a discussion of the convergence in the revised manuscript. revision: yes
Circularity Check
No significant circularity; claims are direct simulation outputs
full rationale
The paper derives ion stratification, sub-stoichiometric Na+/SO4^{2-} ratio, and net-negative interfacial charge from MD trajectories generated by an ML interatomic potential, followed by independent SFG spectrum computation. These quantities are simulation predictions, not inputs redefined or fitted parameters relabeled as outputs. No self-definitional equations, load-bearing self-citations, uniqueness theorems, or ansatz smuggling appear in the abstract or described workflow. The SFG-experiment match serves as post-hoc validation rather than a circular constraint on ion densities. The derivation chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The ML interatomic potential is transferable to the graphene-electrolyte interface.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.