REVIEW 3 major objections 7 minor 38 references
BO-graphane and BO-diamane
T0 review · 3 major / 7 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper predicts that fully covering graphene with boron and oxygen atoms on both sides produces two stable two-dimensional materials, BO-graphane and BO-diamane, that are stiffer than existing diamanes and conduct heat nearly as well…
desk verdict A plausible new 2D carbon material prediction, but the thermodynamic stability claim outruns the evidence. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing motif is the alternating adsorption pattern: on each graphene surface, boron atoms bind to every other carbon atom, and oxygen atoms bridge the boron atoms, producing a honeycomb B–O layer with polar covalent B–O bonds (Bader charge transfer of roughly 1.5 e from boron to oxygen and 0.5 e from boron to carbon). This pattern forces the carbon framework into $\text{sp}^3$ hybridization and, in multilayers, locks the interlayer spacing near diamond values. The stability and property claims are carried by a formation-energy comparison against graphite (the paper's Eq. 1), phonon-dispersion calculations for dynamic stability, ab initio molecular dynamics for thermal stability, finite-strain stress-strain curves for elastic moduli, and an iterative Boltzmann transport solution for lattice thermal conductivity.
What would settle it
A global structure search (for example random sampling or cluster expansion) that finds any B–O covering of monolayer or bilayer graphene with a lower formation energy than the reported 1L or 2L-AB phases would overturn the thermodynamic-stability claim; so would an experiment exposing bilayer graphene to boron and oxygen that yields segregated boron-oxide islands rather than a diamond-like carbon layer with an interlayer spacing near 2.06–2.09 Å.
Extended reading notes
Core claim
On its own terms, the central discovery is that the alternating B–O adsorption motif acts as a structural switch: every other carbon of the graphene lattice moves out of plane and bonds to boron, boron bonds to oxygen, and the carbon lattice rehybridizes to $\text{sp}^3$. In a single layer this yields buckled BO-graphane; in AB- or AA-stacked bilayers and in an ABA trilayer, the carbon layers become diamond-like (C–C bond lengths of 1.53–1.57 Å, C–C–C angles near 109.8–109.9°, interlayer spacing 2.065–2.093 Å), forming a two-dimensional diamond film sandwiched between covalent boron-oxygen layers. The paper reports that these compounds have lower formation energies than H-, F-, and OH-diamane; show no imaginary phonon frequencies; remain intact in AIMD at 1000 K; have in-plane elastic moduli of 659–771 GPa for the monolayer/bilayer and up to about 804 GPa (uniaxial) and 907 GPa (biaxial) for the trilayer; have indirect band gaps of 4.22 eV (1L), 3.87–3.91 eV (2L), and 3.55 eV (3L); and have lattice thermal conductivities of 879 and 1260 W/m·K for 1L and 2L-AB. Single-side-covered graphene (ss-1L) breaks down at 300 K, which the authors take as evidence that two-sided B–O coverage is required for stability.
Load-bearing premise
The assumption is that the alternating arrangement of boron and oxygen atoms the authors place by hand on the graphene surface is the lowest-energy arrangement, since no computational search over other B–O patterns or stoichiometries was run.
Editorial extensions
If this is right
- Boron and oxygen functionalization would provide a low-pressure chemical route to stable diamond-like carbon films, in contrast to the high pressures usually required to convert graphite to diamond.
- The predicted elastic moduli (750–771 GPa for bilayer BO-diamane, rising to about 804–907 GPa for the trilayer) exceed those of F-, OH-, and H-diamane of similar thickness, suggesting use in ultra-thin protective coatings and composites.
- Room-temperature lattice thermal conductivities of 879 W/m·K (1L) and 1260 W/m·K (2L-AB) beat BeO, MgO, Al$_2$O$_3$, and F-diamane and approach H-diamane, making these sheets candidates for heat spreading.
- The band gap is indirect and wide (3.55–4.22 eV) through three layers, pointing toward nanoelectronics and photocatalysis applications.
- Because single-side coverage is thermally unstable at 300 K, any successful synthesis would need to passivate both graphene faces.
Reading between the lines
- A natural next step is a global search over other B–O arrangements, because if a different coverage pattern were lower in energy, the predicted phases would not form; the paper does not include such a search.
- The strongly covalent, charge-polarized B–O cap could function as a chemically passivating shell, plausibly giving air or moisture stability similar to that reported for F-diamane; the paper does not test reactivity or environmental stability.
- The combination of a wide electronic gap with very high lattice thermal conductivity is unusual; if confirmed experimentally, BO-diamane could serve simultaneously as a heat spreader and a transparent wide-gap substrate, a use the paper does not discuss.
- A practical synthesis route may be sequential exposure of bilayer graphene to boron and oxygen precursors, and the predicted stability at 1000 K suggests a wide processing window; the paper stops at prediction.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports first-principles DFT predictions for fully boron- and oxygen-functionalized monolayer graphene (BO-graphane) and multilayer graphene (BO-diamane), with a buckled graphane-like single layer and diamond-like carbon layers in the multilayer cases. It claims thermodynamic and dynamic stability based on formation energies, phonon dispersions, and AIMD at 1000 K, and reports high Young's moduli, wide indirect band gaps (3.55-4.22 eV), and high lattice thermal conductivities (879 and 1260 W/m·K for 1L and 2L-AB, respectively). The methods are standard (VASP, PHONOPY, ShengBTE) and the paper includes comparisons to H-, F-, and OH-diamane. The central stability claim, however, rests on formation energies relative only to elemental reservoirs, without considering competing phases such as B2O3 or alternative B/O decoration patterns, and the abstract's headline elastic modulus values do not match the table entries for the named structures.
Significance. If the predicted structures are indeed stable or even kinetically accessible, BO-graphane and BO-diamane would be interesting additions to the family of functionalized graphene/diamane materials, with potentially attractive mechanical and thermal properties for nanoscale applications. The paper's strengths include the systematic use of standard DFT workflows, phonon stability checks, AIMD thermal stability, and full iterative BTE thermal conductivity calculations, as well as explicit comparisons with existing diamanes. The property predictions are plausible as a computational exploration, but the thermodynamic stability claim is under-supported by the current formation-energy analysis, and this gap is load-bearing for the paper's main narrative.
major comments (3)
- [Abstract and Table 2] The abstract states that BO-graphane has a Young's modulus of 750 GPa and BO-diamane 771 GPa, but Table 2 reports for 1L (BO-graphane) values of 659.66 (x) and 670.58 (y) GPa, and for 2L-AB values of 750.63 (x) and 755.58 (y) GPa, while 2L-AA has y = 771.13 GPa. The abstract numbers thus correspond to specific axes of specific multilayer structures, not to the named 1L and generic BO-diamane structures. This internal inconsistency undermines the headline numerical claims and must be corrected.
- [Section 3.1, Eq. (1)] The formation energy is defined relative to graphite, alpha-boron, and O2 only. This does not establish thermodynamic stability because decomposition into more stable compounds, notably B2O3, is not considered. For the C2B2O2 composition, the reaction to graphite plus B2O3 plus elemental boron is likely exothermic given the high stability of B2O3. Also, no global structure search is performed over alternative B/O decoration patterns (e.g., O bonded directly to C, B-O rows, or clustered boron oxides), so a lower-energy functionalization pattern cannot be ruled out. The claim that these structures are 'thermodynamically stable' is therefore overstated; the authors should either add a convex-hull analysis or a comparison to B2O3 and other plausible competing phases, or rephrase the claim to 'metastable' or 'kinetically accessible'.
- [Section 3.2 and Table 2] The lattice thermal conductivity is reported in W/m·K, which requires a well-defined film thickness. The manuscript specifies the thickness convention for elastic moduli (outermost atoms plus van der Waals radii) but does not state whether the same or a different thickness is used for thermal conductivity. Without a clearly stated thickness, the comparisons to bulk BeO, MgO, and Al2O3 are ambiguous. Please specify the thickness used and, if appropriate, also report the 2D thermal conductance in W/(m·K) per layer.
minor comments (7)
- [Abstract] The unit 'GPA' should be 'GPa' in the abstract.
- [References] References [37] and [38] appear to be the same paper (Mortazavi et al., Appl. Surf. Sci. 528, 147035); please deduplicate and renumber.
- [Table 2] The row for OH-AB does not list the DS and TS entries that appear for the other structures, and the meaning of 'yes [37]' for F-AB should be clarified in a footnote.
- [Figure 2 caption] The caption says 'Formation energy ... with respect to graphite,' but Eq. (1) defines formation energy relative to graphite, alpha-B, and O2; the caption and equation should be made consistent.
- [Section 2] For the second-order force constants used in PHONOPY, the supercell size is not specified; please provide this detail.
- [Section 3.5] Band gaps are computed with the PBE functional, which typically underestimates gap sizes; a brief caveat or a hybrid-functional check would strengthen the 'wide band gap' claim.
- [Section 3.4 and Fig. S2(a)] The AIMD simulation for ss-1L is only 20 ps; the statement that it 'does not maintain structural integrity' should acknowledge the finite simulation time.
Circularity Check
No significant circularity: all predicted properties are direct first-principles outputs, with no target-property fitting and no load-bearing self-citations.
full rationale
The paper's load-bearing predictions are direct outputs of DFT geometry optimization, DFPT phonon calculations, AIMD, and iterative BTE thermal-conductivity calculations. Formation energies use Eq. (1), a standard elemental-reference definition; the comparisons among 1L, 2L-AB, 2L-AA, 3L, ss-1L, and known diamanes are obtained by evaluating the same formula for distinct structures, so the ranking is not an input. Elastic moduli are obtained by linear regression of computed stress-strain curves within 2% strain; band gaps are from PBE band structures; thermal conductivity is from ShengBTE with computed second- and third-order force constants. No parameter is fitted to any predicted property, and no uniqueness theorem or disfavored-alternative argument is imported from prior work. The hand-constructed B-O decoration is an ansatz, but it is an assumption about the structure, not a circular derivation of its properties. The manuscript does not perform a global structure search or ternary convex-hull analysis; that is a completeness and correctness limitation, not circularity. Self-citations (e.g., refs. [10] and [18]) are contextual or experimental benchmarks and do not support any derivation by themselves. Therefore no circular step can be exhibited.
Assumptions & free parameters
free parameters (1)
- Thickness definition for 3D elastic modulus =
Not specified; uses outermost atoms plus van der Waals radii
assumptions (2)
- domain assumption PBE exchange-correlation functional and DFT-D3 dispersion correction accurately describe the structure, stability, and properties of these B/O functionalized carbon nanosheets.
- domain assumption The hand-constructed alternating B-O bonding pattern represents the ground-state functionalization arrangement.
invented entities (2)
-
BO-graphane (1L)
-
BO-diamane (2L-AB, 2L-AA, 3L)
Cite this review
Pith. "Pith review of BO-graphane and BO-diamane." pith.science (2026). https://pith.science/paper/6QPBMFNK
@misc{pith2026250605578,
author = {Pith},
title = {Pith review of: BO-graphane and BO-diamane},
year = {2026},
howpublished = {\url{https://pith.science/paper/6QPBMFNK}},
note = {Machine review of arXiv:2506.05578}
}
abstract
The adsorption of boron and oxygen atoms onto mono- and multi-layer graphene leads to the formation of a buckled graphene layer (BO-graphane) and a 2D diamond-like structure (BO-diamane) sandwiched between boron monoxide layers per DFT calculations. BO-graphane has a calculated Young's modulus ($\it{E}$) of 750 GPA and BO-diamane 771 GPa, higher than the calculated $\it{E}$ of -F,-OH, and -H diamanes; this is due to the presence of B-O bonds in the functionalizing layers. Electronic band structure calculations show BO-graphane and BO-diamane are wide band gap semiconductors with an indirect band gap up to a thickness of three layers (3L). Phonon dispersion and $ab-initio$ molecular dynamics (AIMD) simulations confirm dynamic and thermal stability, maintaining structural integrity at 1000 K. The room-temperature lattice thermal conductivity of BO-graphane and BO-diamane is found to be 879 W/m.K and 1260 W/m.K, respectively, surpassing BeO (385 W/m.K), MgO (64 W/m.K), and Al$_2$O$_3$ (36 W/m.K); and F-diamane (377 W/m.K), and comparable to H-diamane (1145-1960 W/m.K), suggesting them as candidates for thermal management in applications.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
R. H. Baughman, H. Eckhardt, M. Kertesz, Structure-property predictions for new planar forms of carbon: Layered phases con- taining sp2 and sp atoms, J. Chem. Phys. 87 (1987) 6687–6699
work page 1987
- [2]
-
[3]
V . H. Crespi, L. X. Benedict, M. L. Cohen, S. G. Louie, Pre- diction of a pure-carbon planar covalent metal, Phys. Rev. B 53 (1996) R13303–R13305
work page 1996
- [4]
-
[5]
and Mao, Ho-kwang and Eng, Peter J
Mao, Wendy L. and Mao, Ho-kwang and Eng, Peter J. and Trainor, Thomas P . and Newville, Matthew and Kao, Chi-chang and Heinz, Dion L. and Shu, Jinfu and Meng, Y ue and Hemley, Russell J., Bonding changes in compressed superhard graphite, Science 302 (2003) 425–427
work page 2003
-
[6]
Q. Song, B. Wang, K. Deng, X. Feng, M. Wagner, J. D. Gale, K. Mullen, L. Zhi, Graphenylene, a unique two-dimensional carbon network with nondelocalized cyclohexatriene units, J. Mater. Chem. C 1 (2013) 38–41
work page 2013
-
[7]
B. R. Sharma, A. Manjanath, A. K. Singh, pentahexoctite: A new two-dimensional allotrope of carbon, Scientific reports 4 (2014)
work page 2014
-
[8]
Z. Wang, X.-F. Zhou, X. Zhang, Q. Zhu, H. Dong, M. Zhao, A. R. Oganov, Phagraphene: A low-energy graphene allotrope composed of 5–6–7 carbon rings with distorted dirac cones, Nano Lett. 15 (2015) 6182–6186
work page 2015
Show all 38 references
-
[9]
Yin, Y .-P
W.-J. Yin, Y .-P . Chen, Y .-E. Xie, L.-M. Liu, S. Zhang, A low- surface energy carbon allotrope: the case for bcc-C6, Phys. Chem. Chem. Phys. 17 (2015) 14083–14087
2015
-
[10]
B. Ram, H. Mizuseki, Tetrahexcarbon: A two-dimensional al- lotrope of carbon, Carbon 137 (2018) 266–273
2018
-
[11]
C 568: A new two-dimensional sp 2-sp3 hybridized allotrope of carbon, Carbon 158 (2020) 827–835
2020
-
[12]
Lavini, M
F. Lavini, M. Rejhon, E. Riedo, Two-dimensional diamonds from sp 2-to-sp3 phase transitions, Nat. Rev. Mater. 7 (2022) 814–832
2022
-
[13]
A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys. 81 (2009) 109–162
2009
-
[14]
J. O. Sofo, A. S. Chaudhari, G. D. Barber, Graphane: A two- dimensional hydrocarbon, Phys. Rev. B 75 (2007) 153401
2007
-
[16]
D. C. Elias, R. R. Nair, T. M. G. Mohiuddin, S. V . Morozov, P . Blake, M. P . Halsall, A. C. Ferrari, D. W. Boukhvalov, M. I. Katsnelson, A. K. Geim, K. S. Novoselov, Control of graphene’s properties by reversible hydrogenation: Evidence for graphane, Science 323 (2009) 610–613
2009
-
[17]
L. A. Chernozatonskii, P . B. Sorokin, A. G. Kvashnin, D. G. Kvashnin, Diamond-like C 2H nanolayer, diamane: Simulation of the structure and properties, JETP Letters 90 (2009) 134–138
2009
-
[18]
P . V . Bakharev, M. Huang, M. Saxena, S. W. Lee, S. H. Joo, S. O. Park, J. Dong, D. C. Camacho-Mojica, S. Jin, Y . Kwon, M. Biswal, F. Ding, S. K. Kwak, Z. Lee, R. S. Ruoff, Chemically induced transformation of chemical vapour deposition grown bi- layer graphene into fluorinate...
2020
-
[19]
A. G. Kvashnin, P . V . Avramov, D. G. Kvashnin, L. A. Cher- nozatonskii, P . B. Sorokin, Features of electronic, mechanical, and electromechanical properties of fluorinated diamond films of nanometer thickness, J. Phys. Chem. C 121 (2017) 28484– 28489
2017
-
[20]
Clark, K.-J
S. Clark, K.-J. Jeon, J.-Y . Chen, C.-S. Y oo, Few-layer graphene under high pressure: Raman and X-ray di ffraction studies, Solid 11 State Commun. 154 (2013) 15–18
2013
-
[21]
R. R. Nair, W. Ren, R. Jalil, I. Riaz, V . G. Kravets, L. Brit- nell, P . Blake, F. Schedin, A. S. Mayorov, S. Y uan, M. I. Kat- snelson, H.-M. Cheng, W. Strupinski, L. G. Bulusheva, A. V . Okotrub, I. V . Grigorieva, A. N. Grigorenko, K. S. Novoselov, A. K. Geim, Fluorographe...
2010
-
[22]
W. Li, J. Carrete, N. A. Katcho, N. Mingo, ShengBTE: A solver of the Boltzmann transport equation for phonons, Com- put. Phys. Commun. 185 (2014) 1747–1758
2014
-
[23]
Kresse, J
G. Kresse, J. Furthm ¨uller, E fficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set, Comput. Mater. Sci. 6 (1996) 15 – 50
1996
-
[24]
Kresse, J
G. Kresse, J. Furthm ¨uller, E fficient iterative schemes for ab ini- tio total-energy calculations using a plane-wave basis set, Phys. Rev. B 54 (1996) 11169–11186
1996
-
[25]
P . E. Bl ¨ochl, Projector augmented-wave method, Phys. Rev. B 50 (1994) 17953–17979
1994
-
[26]
Kresse, D
G. Kresse, D. Joubert, From ultrasoft pseudopotentials to the projector augmented-wave method, Phys. Rev. B 59 (1999) 1758–1775
1999
-
[27]
H. J. Monkhorst, J. D. Pack, Special points for brillouin-zone integrations, Phys. Rev. B 13 (1976) 5188–5192
1976
-
[28]
A. Togo, I. Tanaka, First principles phonon calculations in ma- terials science, Scr. Mater. 108 (2015) 1–5
2015
-
[29]
Grimme, Semiempirical GGA-type density functional con- structed with a long-range dispersion correction, J
S. Grimme, Semiempirical GGA-type density functional con- structed with a long-range dispersion correction, J. Comput. Chem. 27 (2006) 1787–1799
2006
-
[30]
Topsakal, S
M. Topsakal, S. Cahangirov, S. Ciraci, The response of mechan- ical and electronic properties of graphane to the elastic strain, Appl. Phys. Lett. 96 (2010) 091912
2010
-
[31]
Raeisi, B
M. Raeisi, B. Mortazavi, E. V . Podryabinkin, F. Shojaei, X. Zhuang, A. V . Shapeev, High thermal conductivity in semi- conducting janus and non-janus diamanes, Carbon 167 (2020) 51–61
2020
-
[32]
L. Zhu, W. Li, F. Ding, Giant thermal conductivity in diamane and the influence of horizontal reflection symmetry on phonon scattering, Nanoscale 11 (2019) 4248–4257
2019
-
[33]
Mortazavi, F
B. Mortazavi, F. Shojaei, T. Rabczuk, X. Zhuang, High tensile strength and thermal conductivity in BeO monolayer: A first- principles study, FlatChem 28 (2021) 100257
2021
-
[34]
H. Xie, W. Y u, W. C. and, MgO nanofluids: higher thermal conductivity and lower viscosity among ethylene glycol-based nanofluids containing oxide nanoparticles, J. Exp. Nanosci. 5 (2010) 463–472
2010
-
[35]
H. Ying, A. Moore, J. Cui, Y . Liu, D. Li, S. Han, Y . Y ao, Z. Wang, L. Wang, S. Chen, Tailoring the thermal transport properties of monolayer hexagonal boron nitride by grain size engineering, 2d Mater. 7 (2019) 015031
2019
-
[36]
Zhang, D
X. Zhang, D. Sun, Y . Li, G.-H. Lee, X. Cui, D. Chenet, Y . Y ou, T. F. Heinz, J. C. Hone, Measurement of Lateral and Inter- facial Thermal Conductivity of Single- and Bilayer MoS 2 and MoSe2 Using Refined Optothermal Raman Technique, ACS Appl. Mater. Interfaces 7 (2015) 25923–25929
2015
-
[38]
Mortazavi, F
B. Mortazavi, F. Shojaei, B. Javvaji, M. Azizi, H. Zhan, T. Rabczuk, X. Zhuang, First-principles investigation of me- chanical, electronic and optical properties of H-, F- and Cl- diamane, Appl. Surf. Sci. 528 (2020) 147035
2020
-
[39]
P . Zhu, S. Wang, X. Zhang, J. Zhao, W. Y u, H. Zhang, Mechani- cal properties of diamane kirigami under tensile deformation, J. Nanopart. Res. 26 (2024) 91
2024
-
[40]
Piazza, K
F. Piazza, K. Gough, M. Monthioux, P . Puech, I. Gerber, R. Wiens, G. Paredes, C. Ozoria, Low temperature, pressureless sp2 to sp 3 transformation of ultrathin, crystalline carbon films, Carbon 145 (2019) 10–22. 12 Supplementary Information BO-graphane and BO-diamane Babu Ram ∗...
2019
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.