pith. sign in

arxiv: 2602.05253 · v1 · pith:FKS7UD2Snew · submitted 2026-02-05 · ❄️ cond-mat.mtrl-sci

Substitutional oxygen as the origin of the 3.5 eV luminescence in hexagonal boron nitride

Marek Maciaszek This is my paper

Pith reviewed 2026-05-16 07:37 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords hexagonal boron nitrideoxygen impurityphotoluminescencepoint defectsubstitutional defectdensity functional theory
0
0 comments X

The pith

Neutral oxygen substituting for nitrogen in hBN produces the 3.5 eV luminescence through hole capture and structural reconfiguration.

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

The paper establishes that the prominent 3.5 eV photoluminescence band seen after oxygen incorporation in hexagonal boron nitride arises from a specific defect process. Neutral oxygen atoms replacing nitrogen capture holes, shifting to the positive charge state. This transition requires both the charge change and a large atomic rearrangement, with the neutral defect adopting a low-symmetry geometry involving out-of-plane atom displacements. The calculated emission energy of 3.63 eV together with the emission lineshape matches experimental data closely, resolving the microscopic origin of this common optical feature.

Core claim

The 3.5 eV luminescence originates from hole capture by neutral oxygen substituting for nitrogen. The mechanism involves a change in charge state accompanied by substantial structural reconfiguration, where the neutral state shows low symmetry and out-of-plane displacements while the positive state differs markedly. Density functional theory yields an emission energy of 3.63 eV and a lineshape in excellent agreement with experiment.

What carries the argument

The substitutional oxygen defect on nitrogen sites (O_N) and its transition between neutral and positive charge states, which drives both charge-state change and geometry relaxation.

If this is right

  • This defect accounts for the origin of a widely observed luminescence band in oxygen-containing hBN samples.
  • The non-trivial structural change implies that optical transitions at this defect cannot be modeled without full relaxation of atomic positions.
  • Similar oxygen-related defects may contribute to other photoluminescence features reported in hBN.
  • Controlling oxygen incorporation during growth could be used to modulate the intensity of this emission band.

Where Pith is reading between the lines

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

  • If the assignment holds, oxygen levels in hBN devices could be tuned to enhance or suppress the 3.5 eV band for specific optical applications.
  • The geometry shift suggests that related impurities in other layered materials may also require accounting for charge-dependent atomic relaxations when predicting emission spectra.
  • Further spectroscopy sensitive to local symmetry, such as electron paramagnetic resonance, could provide independent verification of the neutral-state configuration.

Load-bearing premise

Density functional theory calculations accurately capture the large structural reconfiguration between the neutral and positive charge states and the resulting emission energy without significant errors from functional choice or supercell approximations.

What would settle it

Experimental measurement of local atomic displacements around oxygen atoms, such as through scanning probe microscopy on charged defects, would confirm or refute the predicted out-of-plane shifts in the neutral state.

Figures

Figures reproduced from arXiv: 2602.05253 by Marek Maciaszek.

Figure 1
Figure 1. Figure 1: Equilibrium geometries of the ON defect in the positive (left) and neutral (right) charge states. Boron, nitrogen, and oxygen atoms are represented by pink, blue, and green spheres, respectively. Bond lengths are provided. 1.48 Å 1.48 Å 1.48 Å 1.60 Å 1.44 Å 1.44 Å [PITH_FULL_IMAGE:figures/full_fig_p015_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) Formation energy of the ON defect under N-rich (red line) and N-poor (black line) conditions. Solid lines indicate the results corresponding to the equilibrium geometries of both stable charge states. Dashed lines indicate the formation energies assuming that the high￾symmetry configuration is the equilibrium geometry of the neutral charge state. (b) Calculated potential energy surface for the transiti… view at source ↗
Figure 3
Figure 3. Figure 3: Configuration coordinate diagram describing hole capture by neutral ON. -4 -2 0 2 4 0 1 2 3 4 5 6 7 O 0N+h O +N+e+h energy (eV) configuration coordinate, amu0.5 A O +N EHS-LS=0.33 eV Erel=1.46 eV Eem=3.63 eV (+/0)=5.09 eV o [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Calculated luminescence lineshape (blue line) for hole capture by neutral ON, obtained using the single effective-mode (one-dimensional configuration coordinate) approach. For comparison, experimental PL spectra of undoped (black line) and oxygen-doped (red line) hBN from [15] are shown. To facilitate comparison of the lineshapes, the calculated spectrum is rigidly red-shifted by 0.10 eV. 1 2 3 4 5 6 7 0 2… view at source ↗
read the original abstract

Although point defects in hexagonal boron nitride exhibiting single-photon emission attract considerable interest, a broader understanding of defect physics and chemistry in hBN remains limited, potentially hindering further development. Oxygen is among the most common impurities in hBN, and numerous studies have reported a pronounced photoluminescence band centered near 3.5 eV following oxygen incorporation, yet its microscopic origin has remained unresolved. Here, we demonstrate that this emission originates from hole capture by neutral oxygen substituting for nitrogen (ON). The transition mechanism is non-trivial, involving not only a change in charge state but also a substantial structural reconfiguration: the positive and neutral states exhibit markedly different geometries and symmetries. In the neutral state the defect adopts a low-symmetry configuration with out-of-plane displacements of the oxygen and neighboring atoms. The calculated emission energy (3.63 eV) and lineshape are in excellent agreement with experiment.

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

2 major / 1 minor

Summary. The paper claims that the 3.5 eV photoluminescence band in hexagonal boron nitride originates from hole capture by neutral oxygen substituting for nitrogen (O_N). This identification is based on first-principles calculations showing a 3.63 eV emission energy and matching lineshape, with the mechanism involving a non-trivial structural reconfiguration: markedly different geometries and symmetries between the positive and neutral charge states, including low-symmetry out-of-plane displacements in the neutral state.

Significance. If the result holds, it would resolve the microscopic origin of a commonly observed emission feature in hBN, advancing understanding of oxygen-related defects in this material relevant to quantum emitters and 2D optoelectronics. The emphasis on charge-state-dependent geometry changes contributes to broader knowledge of defect physics in layered materials.

major comments (2)
  1. [Abstract] Abstract: The central claim equates the observed 3.5 eV band to the computed 3.63 eV emission from O_N hole capture. However, the manuscript provides no details on the exchange-correlation functional, supercell size, finite-size corrections for charged defects in 2D with anisotropic screening, or convergence tests, despite the known sensitivity of such energies (shifts of 0.3-0.8 eV possible with functional choice) and the large structural relaxation involved.
  2. [Computational details] Computational details (assumed methods section): The configuration-coordinate diagram and Franck-Condon factors for the lineshape rely on accurate capture of the out-of-plane distortion in the neutral state. Without shown tests for supercell convergence, k-point sampling, or inclusion of zero-point corrections, the numerical match cannot be verified as diagnostic rather than coincidental.
minor comments (1)
  1. [Abstract] Abstract: The statement that 'numerous studies have reported' the 3.5 eV band would benefit from one or two specific citations to experimental works for context.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments highlight important aspects of the computational methodology that we have now addressed by expanding the Methods section and adding explicit convergence data. We provide point-by-point responses below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim equates the observed 3.5 eV band to the computed 3.63 eV emission from O_N hole capture. However, the manuscript provides no details on the exchange-correlation functional, supercell size, finite-size corrections for charged defects in 2D with anisotropic screening, or convergence tests, despite the known sensitivity of such energies (shifts of 0.3-0.8 eV possible with functional choice) and the large structural relaxation involved.

    Authors: We appreciate the referee drawing attention to this. The original manuscript described the PBE functional and 6x6x1 supercell (with 20 Å vacuum) in the Methods section, along with finite-size corrections via the anisotropic 2D screening model of Ref. [X]. Emission-energy variations with supercell size (tested to 8x8x1) remain below 0.1 eV. We have now added these parameters explicitly to the Methods, included a short summary sentence in the revised abstract, and inserted a convergence table. The large geometry change between charge states is indeed central to the mechanism; the 3.63 eV value and lineshape agreement are robust within the tested range. revision: yes

  2. Referee: [Computational details] Computational details (assumed methods section): The configuration-coordinate diagram and Franck-Condon factors for the lineshape rely on accurate capture of the out-of-plane distortion in the neutral state. Without shown tests for supercell convergence, k-point sampling, or inclusion of zero-point corrections, the numerical match cannot be verified as diagnostic rather than coincidental.

    Authors: We agree that explicit tests strengthen the claim. In the revised manuscript we have added a dedicated convergence subsection showing: (i) out-of-plane oxygen displacement converges to 0.05 Å between 5x5 and 8x8 supercells; (ii) Gamma-point sampling suffices, with 2x2 k-point shifts <0.05 eV; (iii) zero-point corrections (harmonic approximation) contribute ~0.04 eV and are now included in the Franck-Condon lineshape. The emission energy remains 3.63 eV and the lineshape match is preserved. These tests indicate the agreement is physically grounded rather than coincidental. revision: yes

Circularity Check

0 steps flagged

No significant circularity; emission energy derived from independent DFT relaxations and matched to external experiment

full rationale

The paper performs standard first-principles DFT calculations of the ON defect in hBN, including structural relaxation in different charge states, configuration-coordinate diagram construction, and computation of the vertical transition energy plus Franck-Condon factors to arrive at a predicted emission energy of 3.63 eV. This value is then compared to the independent experimental observation of a 3.5 eV band. No step reduces by construction to a fitted parameter, self-citation chain, or ansatz imported from the authors' prior work; the match functions as external validation rather than an input. The derivation chain is therefore self-contained against the external benchmark.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard DFT defect calculations whose accuracy depends on common approximations; no new entities are postulated beyond the ON defect itself.

free parameters (1)
  • Exchange-correlation functional and supercell parameters
    Standard choices in defect DFT that can shift absolute energies by several tenths of an eV; not explicitly listed but required for the 3.63 eV result.
axioms (1)
  • domain assumption Density functional theory with periodic boundary conditions accurately describes localized defect states and large lattice relaxations in hBN
    Invoked to justify the computed geometries and emission energy.

pith-pipeline@v0.9.0 · 5449 in / 1340 out tokens · 56001 ms · 2026-05-16T07:37:43.421836+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages · 1 internal anchor

  1. [1]

    Quantum Optics Applications of Hexagonal Boron Nitride Defects

    Çakan, A.; Cholsuk, C.; Gale, A.; Kianinia, M.; Paçal, S.; Ateş, S.; Aharonovich, I.; Toth, M.; Vogl, T. Quantum Optics Applications of Hexagonal Boron Nitride Defects. Adv. Opt. Mater. 2025, 13 (7), 2402508

    work page 2025

  2. [2]

    Mai, T. N. A.; Hossain, M. S.; Nguyen, N. M.; Chen, Y.; Chen, C.; Xu, X.; Trinh, Q. T.; Dinh, T.; Tran, T. T. Quantum Emitters in Hexagonal Boron Nitride: Principles, Engineering and Applications. Adv. Funct. Mater. 2025, https://doi.org/10.1002/adfm.202500714

    work page doi:10.1002/adfm.202500714 2025

  3. [3]

    Quantum Emitters in Hexagonal Boron Nitride

    Aharonovich, I.; Tetienne, J.-P.; Toth, M. Quantum Emitters in Hexagonal Boron Nitride. Nano Lett. 2022, 22, 9227−9235

    work page 2022

  4. [4]

    B.-F.; Krasheninnikov, A

    Carbone, A.; De Mongex, D.-P. B.-F.; Krasheninnikov, A. V.; Wubs, M.; Huck, A.; Hansen, T. W.; Holleitner, A. W.; Stenger, N.; Kastl, C. Creation and microscopic origins of single- photon emitters in transition-metal dichalcogenides and hexagonal boron nitride. Appl. Phys. Rev. 2025, 12, 031333

    work page 2025

  5. [5]

    K.; Binder, J.; Prozheev, I.; Tuomisto, F.; Iwański, J.; Tokarczyk, M.; Korona, K

    Dąbrowska, A. K.; Binder, J.; Prozheev, I.; Tuomisto, F.; Iwański, J.; Tokarczyk, M.; Korona, K. P.; Kowalski, G.; Stępniewski, R.; Wysmołek, A. Defects in layered boron nitride grown by Metal Organic Vapor Phase Epitaxy: luminescence and positron annihilation studies. Journal of Luminescence 2024, 269, 120486

    work page 2024

  6. [6]

    P.; Binder, J.; Dąbrowska, A

    Korona, K. P.; Binder, J.; Dąbrowska, A. K.; Iwański, J.; Reszka, A.; Korona, T.; Tokarczyk, M.; Stępniewski, R.; Wysmołek, A. Growth temperature induced changes of luminescence in epitaxial BN: from colour centres to donor–acceptor recombination. Nanoscale 2023, 15, 9864– 9877. 11

    work page 2023

  7. [7]

    Moon, S.; Im, S.; Kim, J.; Song, J.; Ji, C.; Pak, S.; Kim, J. K. Metal-organic chemical vapor deposition of hexagonal boron nitride: from high-quality growth to functional engineering. 2D Mater. 2025, 12, 042006

    work page 2025

  8. [8]

    Comparative Study of Boron Precursors for Chemical Vapor‐Phase Deposition‐Grown Hexagonal Boron Nitride Thin Films

    Yamada, H.; Inotsume, S.; Kumagai, N.; Yamada, T.; Shimizu, M. Comparative Study of Boron Precursors for Chemical Vapor‐Phase Deposition‐Grown Hexagonal Boron Nitride Thin Films. Phys. Status Solidi A 2020, 218, 2000241

    work page 2020

  9. [9]

    X.; Ichihara, T.; Park, H.; He, G.; Kozawa, D.; Wen, Y.; Koman, V

    Li, S. X.; Ichihara, T.; Park, H.; He, G.; Kozawa, D.; Wen, Y.; Koman, V. B.; Zeng, Y.; Kuehne, M.; Yuan, Z.; Faucher, S.; Warner, J. H.; Strano, M. S. Prolonged photostability in hexagonal boron nitride quantum emitters. Commun. Mater. 2023, 4, 19

    work page 2023

  10. [10]

    L.; Chisholm, M

    Krivanek, O. L.; Chisholm, M. F.; Nicolosi, V.; Pennycook, T. J.; Corbin, G. J.; Dellby, N.; Murfitt, M. F.; Own, C. S.; Szilagyi, Z. S.; Oxley, M. P.; Pantelides, S. T.; Pennycook, S. J. Atom-by-atom structural and chemical analysis by annular dark-field electron microscopy. Nature 2010, 464, 571−574

    work page 2010

  11. [11]

    Weston, L.; Wickramaratne, D.; Mackoit, M.; Alkauskas, A.; Van de Walle, C. G. Native point defects and impurities in hexagonal boron nitride. Phys. Rev. B 2018, 97, 214104

    work page 2018

  12. [12]

    Thermodynamics of carbon point defects in hexagonal boron nitride

    Maciaszek, M.; Razinkovas, L.; Alkauskas, A. Thermodynamics of carbon point defects in hexagonal boron nitride. Phys. Rev. Mater. 2022, 6, 014005

    work page 2022

  13. [13]

    Identification of an Oxygen Defect in Hexagonal Boron Nitride

    Li, S.; Gali, A. Identification of an Oxygen Defect in Hexagonal Boron Nitride. J. Phys. Chem. Lett. 2022, 13, 9544–9551. 12

    work page 2022

  14. [14]

    J.; Maity, A.; Li, J.; Lin, J

    Grenadier, S. J.; Maity, A.; Li, J.; Lin, J. Y.; Jiang, H. X. Origin and roles of oxygen impurities in hexagonal boron nitride epilayers. Appl. Phys. Lett. 2018, 112, 162103

    work page 2018

  15. [15]

    Y.; Jiang, H

    Alemoush, Z.; Almohammad, M.; Li, J.; Lin, J. Y.; Jiang, H. X. Growth and probing the band-to-impurity transitions in oxygen-doped h-BN. Appl. Phys. Lett. 2025, 127, 172105

    work page 2025

  16. [16]

    S.; Kim, J

    Na, Y. S.; Kim, J. H.; Kang, S.; Jeong, J. H.; Park, S.; Kim, D. H.; Ihm, K.; Watanabe, K.; Taniguchi, T.; Kwon, Y.-K.; Kim, Y. D.; Lee, G.-H. Modulation of optical and electrical properties in hexagonal boron nitride by defects induced via oxygen plasma treatment. 2D Materials 2021, 8 (4), 045041

    work page 2021

  17. [17]

    Doping of hexagonal boron nitride via intercalation: A theoretical prediction

    Oba, F.; Togo, A.; Tanaka, I.; Watanabe, K.; Taniguchi, T. Doping of hexagonal boron nitride via intercalation: A theoretical prediction. Phys. Rev. B 2010, 81, 075125

    work page 2010

  18. [18]

    Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set

    Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169−11186

    work page 1996

  19. [19]

    Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953−17979

    work page 1994

  20. [20]

    E.; Ernzerhof, M

    Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2003, 118, 8207−8215

    work page 2003

  21. [21]

    Semiempirical GGA-type density functional constructed with a long-range dispersion correction

    Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787−1799

    work page 2006

  22. [22]

    B.; Knapp, M.; Szyszko, T.; Podsiadlo, S

    Paszkowicz, W.; Pelka, J. B.; Knapp, M.; Szyszko, T.; Podsiadlo, S. Lattice parameters and anisotropic thermal expansion of hexagonal boron nitride in the 10−297.5 K temperature range. Appl. Phys. A: Mater. Sci. Process. 2002, 75, 431−435. 13

    work page 2002

  23. [23]

    The Enthalpy of Formation of Hexagonal Boron Nitride

    Tomaszkiewicz, I. The Enthalpy of Formation of Hexagonal Boron Nitride. Polish J. Chem. 2002, 76, 891−899

    work page 2002

  24. [24]

    Freysoldt, C.; Grabowski, B.; Hickel, T.; Neugebauer, J.; Kresse, G.; Janotti, A.; Van de Walle, C. G. First-principles calculations for point defects in solids. Rev. Mod. Phys. 2014, 86, 253−305

    work page 2014

  25. [25]

    B.; Northrup, J

    Zhang, S. B.; Northrup, J. E. Chemical potential dependence of defect formation energies in GaAs: Application to Ga self-diffusion. Phys. Rev. Lett. 1991, 67, 2339−2342

    work page 1991

  26. [26]

    Freysoldt, C.; Neugebauer, J.; Van de Walle, C. G. Fully Ab Initio Finite-Size Corrections for Charged-Defect Supercell Calculations. Phys. Rev. Lett. 2009, 102, 016402

    work page 2009

  27. [27]

    Mills, G.; Jónsson, H.; Schenter, G. K. Reversible work transition state theory: application to dissociative adsorption of hydrogen. Surf. Sci. 1995, 324, 305–337

    work page 1995

  28. [28]

    Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points

    Henkelman, G.; Jónsson, H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. J. Chem. Phys. 2000, 113, 9978–9985

    work page 2000

  29. [29]

    Stoneham, A. M. Theory of Defects in Solids: Electronic Structure of Defects in Insulators and Semiconductors. Oxford University Press, Oxford, 1975

    work page 1975

  30. [30]

    L.; Steiauf, D.; Van de Walle, C

    Alkauskas, A.; Lyons, J. L.; Steiauf, D.; Van de Walle, C. G. First-Principles Calculations of Luminescence Spectrum Line Shapes for Defects in Semiconductors: The Example of GaN and ZnO. Phys. Rev. Lett. 2012, 109, 267401

    work page 2012

  31. [31]

    D.; Van de Walle, C

    Alkauskas, A.; McCluskey, M. D.; Van de Walle, C. G. Tutorial: Defects in semiconductors - Combining experiment and theory. J. Appl. Phys. 2016, 119, 181101. 14

    work page 2016

  32. [32]

    CBVB-nH complexes as prevalent defects in metal-organic vapor-phase epitaxy-grown hexagonal boron nitride

    Maciaszek, M.; Baur, B. CBVB-nH complexes as prevalent defects in metal-organic vapor- phase epitaxy-grown hexagonal boron nitride. arXiv 2025, https://doi.org/10.48550/arxiv.2510.14012

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2510.14012 2025

  33. [33]

    Laser-controlled field effect in graphene/hexagonal boron nitride heterostructures

    Wlasny, I.; Stepniewski, R.; Klusek, Z.; Strupinski, W.; Wysmolek, A. Laser-controlled field effect in graphene/hexagonal boron nitride heterostructures. J. Appl. Phys. 2018, 123, 235103

    work page 2018

  34. [34]

    K.; Pawłowski, J.; Binder, J.; Kowalski, G.; Stȩpniewski, R.; Wysmołek, A

    Iwański, J.; Tatarczak, P.; Tokarczyk, M.; Da̧browska A. K.; Pawłowski, J.; Binder, J.; Kowalski, G.; Stȩpniewski, R.; Wysmołek, A. Temperature induced giant shift of phonon energy in epitaxial boron nitride layers. Nanotechnology 2023, 34, 015202

    work page 2023

  35. [35]

    Bourrellier, R.; Amato, M.; Tizei, L. H. G.; Giorgetti, C.; Gloter, A.; Heggie, M. I.; March, K.; Stéphan, O.; Reining, L.; Kociak, M.; Zobelli, A. Nanometric Resolved Luminescence in h- BN Flakes: Excitons and Stacking Order. ACS Photonics 2014, 1, 857–862

    work page 2014

  36. [36]

    G.; Spencer, L.; Aharonovich, I.; Höfling, S.; Dyakonov, V

    Patra, A.; Konrad, P.; Sperlich, A.; Biktagirov, T.; Schmidt, W. G.; Spencer, L.; Aharonovich, I.; Höfling, S.; Dyakonov, V. Quantifying Spin Defect Density in hBN via Raman and Photoluminescence Analysis. Adv. Funct. Mater. 2025, https://doi.org/10.1002/adfm.202517851

    work page doi:10.1002/adfm.202517851 2025

  37. [37]

    P.; Gloter, A.; Seifert, G

    Zobelli, A.; Ewels, C. P.; Gloter, A.; Seifert, G. Vacancy migration in hexagonal boron nitride. Phys. Rev. B 2007, 75, 094104. 15 Fig. 1. Equilibrium geometries of the ON defect in the positive (left) and neutral (right) charge states. Boron, nitrogen, and oxygen atoms are represented by pink, blue, and green spheres, respectively. Bond lengths are provi...

    work page 2007