pith. sign in

arxiv: 2407.13224 · v1 · submitted 2024-07-18 · ❄️ cond-mat.mtrl-sci · physics.comp-ph

Vacancy-Induced Boron Nitride Monolayers as Multifunctional Materials for Metal Ion Batteries and Hydrogen Storage Applications

Wadha Alfalasi , Wael Othman , Tanveer Hussain , Nacir Tit This is my paper

Pith reviewed 2026-05-06 18:53 UTC · model claude-opus-4-7

classification ❄️ cond-mat.mtrl-sci physics.comp-ph
keywords storagemetalrespectivelyanalysisapplicationschargeelectronicenergies
0
0 comments X

The pith

Boron-vacancy porous BN monolayers decorated with Li/Na/K are computationally predicted to deliver 490–1822 mAh/g anode capacity and 9.4–10.7 wt% H2 storage at the DFT level.

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

The authors take a recently proposed porous boron-nitride monolayer (a BN analog of holey graphyne, with periodic 6- and 8-vertex rings) and remove a boron atom to create a defect (BN:VB). Using density functional theory with a van der Waals correction, they relax structures, compute band structures, density of states, Bader charges, ab initio molecular dynamics at 400 K, and nudged-elastic-band-style diffusion barriers.

For battery anodes they sequentially add Li, Na, or K atoms until binding energies fall below the bulk cohesive energy of the metal (used as the no-clustering criterion), reaching 44 Li, 28 Na, or 20 K per simulation cell. From these counts they compute specific capacities (1821, 786, 491 mAh/g) and open-circuit voltages (0.15, 0.25, 0.32 V). Diffusion barriers drop with loading, reaching 0.08 eV for Na.

For hydrogen storage they decorate the monolayer with 4 metal atoms and add H2 molecules until average adsorption energy crosses ~0.16 eV, giving 20–24 H2 per cell and 9.4–10.7 wt%. A Langmuir-style grand-canonical partition function with experimental ΔH+TΔS values is then used to estimate effective uptake at adsorption (30 atm, 25 °C) and desorption (3 atm, 100 °C) conditions, yielding 6.3–9.9 wt%.

The work is a fairly standard 2D-materials computational screening exercise; novelty lies in combining the specific BN:VB substrate with both ion-battery and H2-storage analyses.

Core claim

Boron-vacancy porous BN monolayers decorated with Li/Na/K provide theoretical specific capacities of 1821.53, 786.11, and 490.51 mAh/g respectively with low OCVs (0.15–0.32 V) and low diffusion barriers (down to 0.08 eV for Na), and 4-metal-decorated BN:VB stores 9.38–10.72 wt% H2, with effective reversible capacity 6.34–9.89 wt% under practical (30 atm/25 °C ↔ 3 atm/100 °C) conditions — exceeding the DOE 2025 5.5 wt% target.

Load-bearing premise

The maximum metal loading is determined purely by the criterion |E_bind| > |E_cohesive(bulk metal)|, used as a static no-clustering proxy. This neglects (a) kinetic/synthesis pathways to actually deposit 44 Li or 28 Na on a single primitive cell, (b) lateral metal–metal interactions and 2D/3D cluster formation off the surface, (c) electrolyte interactions, and (d) the assumption that the parent porous BN — itself a theoretical structure not yet synthesized — and its B-vacancy variant are realizable. The Langmuir-model H2 effective capacity additionally assumes non-interacting adsorption sites with site energies taken from the DFT static optimization, which is a known overestimate compared to experiment for metal-decorated 2D sorbents.

read the original abstract

This study comprehensively examined the structural, electronic, electrochemical, and energy storage properties of boron-vacancy induced porous boron nitride monolayers (BN:VB) as multifunctional materials, anodes for MIBs and H2 storage applications. Our computational approaches, density functional theory (DFT), ab initio molecular dynamics (AIMD), and thermodynamic analysis, revealed exceptionally high energy and gravimetric densities for MIBs and H2 storage, respectively. We investigated the interactions of Li, Na, and K atoms on BN:VB, which strongly bonded with binding energies stronger than their bulk cohesive energies, which ensured structural stability and the absence of metal clustering. Electronic properties, analyzed through spin-polarized partial density of states (PDOS), band structure, and Bader charge analysis, revealed significant charge transfers from the metal atoms to BN:VB, enhancing the electronic conductivity of the latter. Theoretical specific capacities were calculated as 1821.53, 786.11, and 490.51 mA h/g for Li, Na, and K, respectively, which comfortably exceeded the conventional anodes, such as graphite. Average open-circuit voltages (OCVs) were found as 0.15, 0.25, and 0.32 V, for Li, Na, and K, respectively, indicating strong electrochemical stability. Diffusion studies showed lower barriers of 0.47, 0.08, and 0.60 eV for Li, Na, and K, respectively, with increased metal loadings, suggesting enhanced mobilities and charge/discharge rates. On the other side, the metal-functionalized BN:VB monolayers exhibited remarkably high H2 gravimetric capacities, supported by Langmuir adsorption model-based statistical thermodynamic analysis. Average adsorption energies of H2 on 4Li-, 4Na-, and 4K@BN:VB, were found in perfect range for practical storage applications.

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.

Referee Report

5 major / 9 minor

Summary. The authors use spin-polarized DFT (VASP, PBE-GGA, DFT-D3), AIMD, and a Langmuir-model thermodynamic analysis to evaluate a B-vacancy porous boron nitride monolayer (BN:VB) as (i) an anode for Li/Na/K-ion batteries and (ii) a substrate (4MA@BN:VB, MA=Li,Na,K) for H2 storage. Key claimed numbers: theoretical specific capacities 1821.53/786.11/490.51 mAh/g (Li/Na/K); average OCVs 0.15/0.25/0.32 V; diffusion barriers reduced to 0.47/0.08/0.60 eV at higher loadings; H2 gravimetric capacities 10.64/10.72/9.38 wt% with effective (Langmuir) capacities 9.22/9.89/6.34 wt%. The work concludes BN:VB is a multifunctional energy-storage material exceeding graphite-anode capacities and the DOE 5.5 wt% H2 target.

Significance. If the central numbers are reproducible and the OCV definition is corrected, the manuscript provides a useful comparative DFT screening of vacancy-engineered porous BN against a fairly broad benchmark set (Table 5 is a genuinely useful compilation). Strengths include: (i) inclusion of single-atom, partial- (4MA), and saturation-loading regimes with AIMD checks at 400 K; (ii) explicit Bader and ELF analyses; (iii) the (P,T) Langmuir analysis to convert 0-K capacities to operationally relevant adsorption/desorption windows, which many comparable papers omit; (iv) coverage of both intercalation electrochemistry and physisorptive H2 storage on the same substrate, allowing internal consistency checks. The novelty over Mahamiya et al. (Pr-BN, ref. [49]) is the V_B variant and the metal-loading/H2 study, which is incremental but reasonable.

major comments (5)
  1. [Eq. (3), §3.3.1, Table 2, Fig. 7(b), Abstract] Internal inconsistency in the OCV calculation that is load-bearing for the central anode claim. Eq. (3) defines E_MA as the energy of an 'isolated MA'; with that reference, OCV = |E_bind|/e, which from §3.2.1/Fig. 7(a) gives 1.80/1.36/1.26 V at maximum loading (and 4.45/3.79/4.27 V at single-atom loading) — neither matches the 0.15/0.25/0.32 V quoted in the abstract and Fig. 7(b), nor the 2.22/2.82/2.04 V listed in Table 2. The abstract values are reproduced almost exactly by |E_bind|−|E_cohesive| (i.e., a bulk-metal reference, the standard MIB convention), so Eq. (3) is mis-stated. Please correct the equation, reconcile Table 2 with Fig. 7(b), and clarify which loading the reported OCV averages over. The viability of BN:VB as an anode (OCV ≲ 0.5 V) hinges entirely on this.
  2. [§3.2.1 and Fig. 4(a)] The criterion 'no clustering' is asserted from |E_bind| > |E_cohesive(bulk)| evaluated at sequential layer addition (n up to 4). This is a necessary but not sufficient condition: it does not preclude lateral 2D islanding on the surface, nor does it test stability of the very high saturation loadings (44 Li, 28 Na, 20 K per primitive cell). AIMD validation is shown only for 4MA@BN:VB at 400 K for 8 ps; please add AIMD (and ideally nudged-elastic-band cluster-formation tests) at the saturation loadings used to compute the headline 1821.53 mAh/g, since that capacity is the figure of merit being advertised.
  3. [§3.3.2, Fig. 8] The interpretation that 'increasing the loading of MA in the pore significantly lowers the diffusion barriers' (down to 0.08 eV for Na) needs more care. NEB barriers between increasingly destabilized sites can drop simply because the initial state becomes higher in energy due to MA–MA repulsion, not because transport actually accelerates. Please report (a) the absolute energies of initial/final states relative to the dilute limit, (b) backward barriers, and (c) whether the relevant rate-limiting step at operating concentration is intra-pore hopping or pore-to-pore transfer. Otherwise the 'enhanced charge/discharge rate' claim is not supported.
  4. [§3.4 and Eq. (4)] Average H2 adsorption energies of −0.19 to −0.22 eV with PBE+D3 on metal-cation-decorated 2D substrates are known to be sensitive to the dispersion treatment and to whether the metal–H2 interaction is dominated by Kubas-like polarization (often poorly captured by PBE). Given that the effective capacity in Table 4 is computed by inserting these per-H2 DFT energies directly into Eq. (6), small systematic errors of ~50 meV translate into large changes in N_a−N_d at 25–100 °C. Please test sensitivity to functional (e.g., a vdW-DF or single-point HSE06 on a few configurations) and report how C_E shifts, since the 'exceeds DOE 5.5 wt%' claim is what the title advertises.
  5. [§3.1 and Introduction] The parent Pr-BN structure (ref. [49]) is itself a theoretically predicted, not experimentally synthesized, monolayer, and BN:VB adds a further structural assumption (ordered B-vacancy in every primitive cell). The manuscript should state this caveat plainly in the abstract/conclusion and discuss the realistic vacancy concentration range from refs. [45], [59], rather than treating one ordered V_B per 24-atom cell as given. As written, the comparison to graphite (an experimentally established anode) is somewhat asymmetric.
minor comments (9)
  1. [Abstract and §3.3.1] Specific capacities are reported to two decimal places (1821.53 mAh/g) — this precision is not justified by the underlying DFT or by the ambiguity in the saturation criterion. Round to 4 significant figures at most.
  2. [Table 1] Footnote '*non-polarized band gap energy' is unclear; please define what 'non-polarized band gap' means here (presumably the gap if spin polarization is suppressed) and explain why two different gap values are tabulated for the same row.
  3. [Eq. (5) and Table 4] Define N_T as 'molecules per simulation cell' explicitly, and clarify the conversion to wt% accounts for the full BN:VB primitive cell (12 B + 12 N − 1 B vacancy + 4 MA). A worked example for one system in the SI would help reproducibility.
  4. [§2] k-mesh of 5×5×1 for total-energy convergence on a 24-atom 10.96 Å cell is on the low side once V_B and metal loading break symmetry. Please report a convergence test, or use the 8×8×1 mesh consistently.
  5. [Table 2] Rows are labeled 'Li@BN', 'Na@BN', 'K@BN' — these should presumably be 'Li@BN:VB', etc., to match §3.3.1.
  6. [Fig. 11 and §3.4.2] The Langmuir model assumes equivalent, non-interacting sites. Given that the per-H2 E_ads varies across the 20–24 H2 adsorbed (Fig. 10(a) shows a distribution), please state whether E_b^i in Eq. (6) uses a single average or the full spectrum of site energies; this affects the entropy term.
  7. [References] Refs. [56] and [62] appear to be duplicates (Keshavarz et al., Condens. Matter 7, 8). Ref. [3] and [64] are also nearly duplicate Kaskhedikar/Maier entries.
  8. [Figures 5, 7, 9, 10] Several axis labels and legends are difficult to read at the printed size; ensure tick labels and energy units are legible after typesetting.
  9. [§3.2.2] The statement 'transition of BN:VB from a semiconductor to metal ... driven by the substantial charge transfer' should be qualified: PBE-GGA is well known to underestimate gaps and to spuriously metallize defective insulators. A single HSE06 check on 4Li@BN:VB would be reassuring.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The paper rests on (i) standard Kohn-Sham DFT with PBE-GGA and Grimme-D3 dispersion as adequate for B/N/alkali-metal/H2 energetics in the 0.1–0.6 eV regime, (ii) the assumption that a B-vacancy variant of the theoretical Pr-BN porous BN structure is synthesizable, (iii) the no-clustering criterion |E_bind|>|E_cohesive(bulk)| as a sufficient stability test, (iv) Langmuir non-interacting site model with experimental ΔH+TΔS for translating static DFT adsorption energies into finite-T,P uptake. No new particles, forces, or fitted free parameters are introduced; the work is a parameter-free DFT calculation in the sense that no constants are fit to the target observables. Reproducibility-relevant choices (k-mesh, cutoff, supercell, AIMD time) are stated.

pith-pipeline@v0.9.0 · 9911 in / 6801 out tokens · 99927 ms · 2026-05-06T18:53:37.773026+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.