arxiv: 2604.26264 · v1 · submitted 2026-04-29 · ❄️ cond-mat.mtrl-sci

Recognition: unknown

Accelerated Prediction of Surface Stability and Particle Morphology in Ionic Crystals via Electrostatic Screening

Sourav Baiju , Payam Kaghazchi

Authors on Pith no claims yet

Pith reviewed 2026-05-07 13:19 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci

keywords ionic crystalssurface stabilityelectrostatic screeningcrystal morphologystoichiometric slabspolar surfaceshigh-throughput screeningWulff construction

0 comments

The pith

Electrostatic interactions alone predict the most stable surfaces and shapes of ionic crystals.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors show that ranking surfaces by their electrostatic energy from simple stoichiometric slab models is sufficient to identify the facets that dominate real crystal morphologies in ionic materials. They build candidate terminations, compute electrostatic energies, flag polar surfaces by dipole moment, and use replica-exchange Monte Carlo to reconstruct them into neutral configurations. This bypasses full density-functional calculations and scales to high-index and large-system cases. Validation against published DFT and experimental data for both bulk 3D ionic compounds and 2D layered oxides confirms that the electrostatic ranking reproduces observed trends. The work therefore supplies a fast filter for morphology prediction and a starting point for more accurate studies of energy materials.

Core claim

Electrostatic screening of stoichiometric slab terminations captures the dominant trends in relative surface stability across ionic materials. Polar surfaces are identified by their dipole moments and stabilized through electrostatics-guided reconstructions performed with replica-exchange Monte Carlo. The resulting surface energies yield equilibrium morphologies whose dominant facets match those reported from density functional theory and from experiment on both simple and complex three-dimensional crystals as well as two-dimensional layered oxides. High-index surfaces are shown to play a measurable role in some equilibrium shapes.

What carries the argument

The electrostatic energy of charge-neutral stoichiometric slab models, augmented by surface dipole calculations to detect polarity and replica-exchange Monte Carlo to find minimal-energy reconstructions.

If this is right

High-throughput screening becomes feasible for thousands of surface orientations and large supercells.
Equilibrium crystal shapes can be estimated directly from electrostatic rankings.
High-index surfaces must sometimes be included to obtain correct particle morphologies.
The electrostatic results supply reliable initial structures for subsequent DFT optimizations.
Two-dimensional layered oxides can be treated on the same footing as three-dimensional ionic compounds.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Materials dominated by covalent rather than ionic bonding may require additional correction terms beyond pure electrostatics.
Rapid electrostatic screening could be used to pre-select compositions and facets for experimental synthesis trials.
The same dipole and reconstruction logic might be adapted to predict how surface energies change under applied electric fields or in electrolyte environments.
Combining the method with machine-learned interatomic potentials would further extend its reach to even larger nanoparticles.

Load-bearing premise

The relative stability of surfaces is governed primarily by their electrostatic energies computed on stoichiometric, dipole-corrected slabs, with other quantum-mechanical contributions being secondary.

What would settle it

Finding an ionic material where the surface predicted to have the lowest electrostatic energy is absent from both experimental particle shapes and from explicit DFT calculations of surface energies would falsify the central claim.

read the original abstract

This work presents a fast and scalable approach for predicting surface stability and equilibrium crystal morphology in ionic materials using electrostatic analysis. The method constructs stoichiometric slab terminations and evaluates their electrostatic energies, enabling high-throughput screening of surface configurations at a fraction of the cost of conventional approaches. Polar surfaces are identified through surface dipole moment calculations and stabilized via electrostatics-based reconstruction using replica-exchange Monte Carlo simulations. The surface dipole moment further emerges as an effective descriptor to distinguish the behavior of different classes of materials. By bypassing expensive Density Functional Theory (DFT) calculations, the approach extends naturally to large systems and high-index surfaces that are typically inaccessible to DFT. Electrostatic interactions are shown to capture the dominant trends in relative surface stability across diverse material systems. The method is validated on simple and complex 3D materials as well as 2D layered oxides, where the predicted dominant facets are consistent with reported density functional theory and experimental observations. Importantly, the framework also reveals cases where high-index surfaces play a non-negligible role in the equilibrium morphology. These results establish electrostatics as a fast and reliable route for high-throughput prediction of surface stability and particle morphology, opening a pathway for accelerated materials discovery and providing a robust starting point for more detailed calculations in complex energy materials.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Electrostatic screening gives quick surface stability estimates for ionic crystals that match DFT trends in the tested cases, but the lack of quantitative comparisons leaves the limits unclear when relaxation or covalency matters.

read the letter

This paper gives a fast electrostatic screening method for surface stability and particle morphology in ionic crystals. It constructs stoichiometric slabs, calculates their electrostatic energies, uses surface dipole to spot polar surfaces, and applies replica-exchange Monte Carlo to find stable reconstructions. The dipole also serves as a descriptor for different material classes. For the materials checked, the predicted dominant facets line up with DFT and experiments, and it shows high-index surfaces can be important sometimes. What is new is the end-to-end workflow that skips DFT for the initial screening and handles large systems and high-index faces. It does well at being computationally cheap while still producing results that are consistent with more expensive methods on the examples given. The use of Monte Carlo for electrostatic stabilization is a sensible choice for polar cases. The soft spots are around the assumption that electrostatic terms dominate the relative stabilities. The paper reports consistency but does not show how close the energies are or provide cases where the method disagrees with DFT due to relaxation or bonding effects. Without those, it is hard to know the boundaries of when this works. The concern that covalent or relaxation effects could invert the rankings is reasonable and not directly addressed. This is aimed at researchers doing high-throughput materials discovery for energy applications who need quick filters for surfaces. Someone working on ionic crystals or oxides would get practical value from the method as a starting point, even if they follow up with DFT. It deserves peer review. The idea is grounded and the validation, though limited, is on diverse materials. Referees can ask for more detailed comparisons and quantitative metrics to strengthen it.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a fast electrostatic screening method to predict surface stability and equilibrium morphologies for ionic crystals. Stoichiometric slab terminations are constructed and their electrostatic energies evaluated to rank surfaces; polar surfaces are identified via dipole-moment calculations and stabilized through electrostatics-based reconstructions performed with replica-exchange Monte Carlo. The approach is applied to simple and complex 3D ionic materials as well as 2D layered oxides, with the claim that predicted dominant facets match those obtained from prior DFT calculations and experiments. The work additionally highlights the non-negligible contribution of certain high-index surfaces to the equilibrium shape.

Significance. If the electrostatic ranking proves robust, the method would enable high-throughput screening of surface properties for systems and facet indices that remain inaccessible to routine DFT, supplying inexpensive initial guesses for subsequent electronic-structure refinement in energy materials. The absence of fitted parameters and the direct use of Coulombic slab energies constitute a clear methodological advantage when the approximation holds.

major comments (2)

[Results] The central claim that electrostatic interactions capture the dominant trends in relative surface stability rests on the validation against DFT and experiment, yet the manuscript reports only qualitative consistency of dominant facets without quantitative surface-energy tables, correlation coefficients, or error metrics (Results section). This omission prevents assessment of whether the electrostatic term alone reproduces DFT rankings or whether post-hoc adjustments were required.
[Methods] The method assumes that electrostatic energies computed on stoichiometric slabs (with dipole-based reconstructions for polar cases) are sufficient to determine relative stabilities. The manuscript does not present explicit tests or counter-examples in which covalent bonding, surface relaxation, or charge transfer invert the electrostatic ordering relative to full DFT (Methods and Discussion sections). Such tests are load-bearing for the claim of broad applicability across diverse ionic chemistries.

minor comments (2)

The workflow diagram (if present) or a concise pseudocode block would improve clarity of the slab-construction and Monte Carlo reconstruction steps.
[Methods] Ensure that all acronyms (DFT, MC, etc.) are defined on first use and that the precise definition of the electrostatic energy functional is stated explicitly, including any truncation or periodicity conventions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed feedback. We address each major comment below and have revised the manuscript to improve the quantitative validation and clarify the scope of applicability.

read point-by-point responses

Referee: [Results] The central claim that electrostatic interactions capture the dominant trends in relative surface stability rests on the validation against DFT and experiment, yet the manuscript reports only qualitative consistency of dominant facets without quantitative surface-energy tables, correlation coefficients, or error metrics (Results section). This omission prevents assessment of whether the electrostatic term alone reproduces DFT rankings or whether post-hoc adjustments were required.

Authors: We agree that quantitative comparisons would strengthen the validation. In the revised manuscript, we have added a table in the Results section listing electrostatic surface energies alongside available DFT values for the studied materials. We also report Pearson correlation coefficients between the electrostatic and DFT rankings for multi-facet systems, confirming that the electrostatic term reproduces the DFT ordering without post-hoc adjustments. revision: yes
Referee: [Methods] The method assumes that electrostatic energies computed on stoichiometric slabs (with dipole-based reconstructions for polar cases) are sufficient to determine relative stabilities. The manuscript does not present explicit tests or counter-examples in which covalent bonding, surface relaxation, or charge transfer invert the electrostatic ordering relative to full DFT (Methods and Discussion sections). Such tests are load-bearing for the claim of broad applicability across diverse ionic chemistries.

Authors: The presented validation covers a range of ionic crystals where electrostatics dominate, as shown by agreement with DFT and experiment. We have expanded the Discussion to address potential limits, noting that in more covalent systems other effects could dominate and that the approach serves as a rapid screening tool. While a systematic set of counter-examples would require new DFT calculations outside the current scope, the tested cases and literature context support applicability within the ionic regime. revision: partial

Circularity Check

0 steps flagged

No circularity: direct electrostatic computation independent of target outcomes

full rationale

The paper computes surface electrostatic energies on stoichiometric slabs (with dipole-based reconstructions for polar cases) and uses the resulting dipole moment as a derived descriptor. These quantities are obtained from first-principles electrostatics without fitting parameters to the reported surface stabilities or morphologies. Validation consists of external consistency checks against independent DFT calculations and experiments rather than internal self-consistency or self-citation chains. No step renames a fitted input as a prediction, imports uniqueness from prior author work, or smuggles an ansatz via citation. The derivation chain therefore remains self-contained and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Information extracted from abstract only; no explicit free parameters, invented entities, or detailed axioms visible in the summary.

axioms (1)

domain assumption Electrostatic interactions capture the dominant trends in relative surface stability for ionic crystals
This premise allows bypassing DFT as stated in the abstract.

pith-pipeline@v0.9.0 · 5526 in / 1247 out tokens · 76465 ms · 2026-05-07T13:19:51.824409+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

25 extracted references · 1 canonical work pages

[1]

Physical review82(1), 87 (1951)

Herring, C.: Some theorems on the free energies of crystal surfaces. Physical review82(1), 87 (1951)

1951
[2]

Dantzig, J.: Equilibrium shapes for crystalline solids
[3]

Surface science reports19(3-6), 191–229 (1993)

Van Hove, M.A., Moritz, W., Over, H., Rous, P.J., Wander, A., Barbieri, A., Materer, N., Starke, U., Somorjai, G.A.: Automated determination of complex surface structures by leed. Surface science reports19(3-6), 191–229 (1993)

1993
[4]

Nature communications10(1), 269 (2019)

Luo, Y., Tang, L., Khan, U., Yu, Q., Cheng, H.-M., Zou, X., Liu, B.: Morphology and surface chemistry engineering toward ph-universal catalysts for hydrogen evolution at high current density. Nature communications10(1), 269 (2019)

2019
[5]

Chemical Society Reviews (2025)

Hossain, M.N., Zhang, L., Neagu, R., Sun, S.: Exploring the properties, types, and performance of atomic site catalysts in electrochemical hydrogen evolution reactions. Chemical Society Reviews (2025)

2025
[6]

Physical Review B—Condensed Matter and Materials Physics76(16), 165435 (2007)

Wang, L., Zhou, F., Meng, Y., Ceder, G.: First-principles study of surface prop- erties of li fe po 4: surface energy, structure, wulff shape, and surface redox potential. Physical Review B—Condensed Matter and Materials Physics76(16), 165435 (2007)

2007
[7]

Electrochemical and Solid-State Letters11(6), 94 (2008)

Wang, L., Zhou, F., Ceder, G.: Ab initio study of the surface properties and nanoscale effects of limnpo4. Electrochemical and Solid-State Letters11(6), 94 (2008)

2008
[8]

Physical Chemistry Chemical Physics16(39), 21788–21794 (2014)

Whiteside, A., Fisher, C.A., Parker, S.C., Islam, M.S.: Particle shapes and sur- face structures of olivine nafepo 4 in comparison to lifepo 4. Physical Chemistry Chemical Physics16(39), 21788–21794 (2014)

2014
[9]

Journal of solid state electrochemistry18(5), 1401–1413 (2014)

H¨ ormann, N., Groß, A.: Stability, composition and properties of li2fesio4 surfaces studied by dft. Journal of solid state electrochemistry18(5), 1401–1413 (2014)

2014
[10]

Journal of Physics: Energy3(4), 044007 (2021)

Qu´ erel, E., Seymour, I.D., Cavallaro, A., Ma, Q., Tietz, F., Aguadero, A.: The role of nasicon surface chemistry in stabilizing fast-charging na metal solid-state batteries. Journal of Physics: Energy3(4), 044007 (2021)

2021
[11]

Computational materials science11(1), 16–26 (1998)

Roosen, A.R., McCormack, R.P., Carter, W.C.: Wulffman: A tool for the calcula- tion and display of crystal shapes. Computational materials science11(1), 16–26 (1998)

1998
[12]

The Journal of Physical Chemistry C117(41), 21274–21280 (2013)

Zhang, W.-B., Chen, C., Zhang, S.-Y.: Equilibrium crystal shape of ni from first principles. The Journal of Physical Chemistry C117(41), 21274–21280 (2013)

2013
[13]

15 Physical review letters90(4), 046103 (2003)

Reuter, K., Scheffler, M.: First-principles atomistic thermodynamics for oxidation catalysis:¡? format?¿ surface phase diagrams and catalytically interesting regions. 15 Physical review letters90(4), 046103 (2003)

2003
[14]

npj Computational Materials7(1), 58 (2021)

Yoo, S.-H., Todorova, M., Wickramaratne, D., Weston, L., Walle, C.G.V.d., Neugebauer, J.: Finite-size correction for slab supercell calculations of materials with spontaneous polarization. npj Computational Materials7(1), 58 (2021)

2021
[15]

Journal of Physics C: Solid State Physics12(22), 4977 (1979)

Tasker, P.: The stability of ionic crystal surfaces. Journal of Physics C: Solid State Physics12(22), 4977 (1979)

1979
[16]

https://doi

Jain, A., Ong, S.P., Hautier, G., Chen, W., Richards, W.D., Dacek, S., Cholia, S., Gunter, D., Skinner, D., Ceder, G., Persson, K.a.: The Materials Project: A materials genome approach to accelerating materials innovation. APL Materials 1(1), 011002 (2013) https://doi.org/10.1063/1.4812323

work page doi:10.1063/1.4812323 2013
[17]

Journal of Physics: Condensed Matter29(27), 273002 (2017)

Larsen, A.H., Mortensen, J.J., Blomqvist, J., Castelli, I.E., Christensen, R., Du lak, M., Friis, J., Groves, M.N., Hammer, B., Hargus, C.,et al.: The atomic simulation environment—a python library for working with atoms. Journal of Physics: Condensed Matter29(27), 273002 (2017)

2017
[18]

npj Computational Materials11(1), 202 (2025)

K¨ oster, K., Binninger, T., Kaghazchi, P.: Optimization of coulomb energies in gigantic configurational spaces of multi-element ionic crystals. npj Computational Materials11(1), 202 (2025)

2025
[19]

Beilstein journal of nanotechnology6(1), 361–368 (2015)

Barmparis, G.D., Lodziana, Z., Lopez, N., Remediakis, I.N.: Nanoparticle shapes by using wulff constructions and first-principles calculations. Beilstein journal of nanotechnology6(1), 361–368 (2015)

2015
[20]

Nanoscale6(22), 13391–13395 (2014)

Chen, Y.-X., Kaghazchi, P.: Metalization of li 2 s particle surfaces in li–s batteries. Nanoscale6(22), 13391–13395 (2014)

2014
[21]

Nature Energy2(7), 1–10 (2017)

Tan, G., Xu, R., Xing, Z., Yuan, Y., Lu, J., Wen, J., Liu, C., Ma, L., Zhan, C., Liu, Q.,et al.: Burning lithium in cs2 for high-performing compact li2s–graphene nanocapsules for li–s batteries. Nature Energy2(7), 1–10 (2017)

2017
[22]

Journal of Nanomaterials2016(1), 8123652 (2016)

Liu, Y., Yang, Y.: Recent progress of tio2-based anodes for li ion batteries. Journal of Nanomaterials2016(1), 8123652 (2016)

2016
[23]

Applied Surface Science436, 989–994 (2018)

Jiang, F., Yang, L., Zhou, D., He, G., Zhou, J., Wang, F., Chen, Z.-G.: First- principles atomistic wulff constructions for an equilibrium rutile tio2 shape modeling. Applied Surface Science436, 989–994 (2018)

2018
[24]

Chemistry of Materials21(16), 3799–3809 (2009)

Kramer, D., Ceder, G.: Tailoring the morphology of licoo2: a first principles study. Chemistry of Materials21(16), 3799–3809 (2009)

2009
[25]

ACS Applied Materials & Interfaces4(6), 3062–3068 (2012) 16

Jiang, J., Liu, W., Chen, J., Hou, Y.: Lifepo4 nanocrystals: liquid-phase reduc- tion synthesis and their electrochemical performance. ACS Applied Materials & Interfaces4(6), 3062–3068 (2012) 16

2012