REVIEW 3 major objections 3 minor 36 references
Defect physics in $Yb^{3+}$-doped $CaF_2$ from first-principles calculation
T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper argues that fluorine-rich growth conditions can enhance Yb:CaF2 luminescence by suppressing deep trap defects while promoting beneficial Yb clustering, based on first-principles defect energetics.
desk verdict Solid first-principles defect dataset for Yb:CaF2, but the F-rich growth recommendation rests on a Fermi-level proxy for the Yb3+ levels that needs a firmer total-energy basis. 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 key machinery is the supercell defect-formation-energy and thermodynamic-transition-level formalism of Freysoldt et al. (equation 1), combined with a Fermi-energy assignment of the Yb3+ ground and excited manifolds and with pair-defect binding energies defined as $E_b = H_f[A] + H_f[B] - H_f[AB]$. These quantities let the authors compare defect stability across fluorine-rich and fluorine-poor chemical-potential limits, identify which defects have trap levels overlapping the Yb optical levels, and quantify whether isolated defects prefer to aggregate into the experimentally favored monomers and dimers.
What would settle it
A decisive experiment would be to grow Yb:CaF2 under carefully controlled fluorine-rich and fluorine-poor conditions, then measure the relative concentrations of $\mathrm{Yb_{Ca}}$ and $\mathrm{Yb_i{-}F_i}$ (for example by site-selective spectroscopy or paramagnetic resonance) alongside the up/down-conversion quantum efficiency. If fluorine-rich samples do not show reduced quenching or enhanced Yb clustering, or if direct spectroscopy of the 4f manifold places the $^{2}F_{5/2}$ level well away from the predicted trap levels, the central claim would be contradicted.
Extended reading notes
Core claim
The paper's central claim is that the luminescence-quenching defects in Yb-doped CaF2 are the isolated antisite $\mathrm{Yb_{Ca}}$ and the pair $\mathrm{Yb_i{-}F_i}$, both of which create deep trap levels near the optically excited $^{2}F_{5/2}$ level of Yb3+. Using density functional theory with PBEsol and HSE06 functionals, the authors place the Yb3+ ground level $^{2}F_{7/2}$ at Fermi energies of 1.03 eV (PBEsol) and 2.76 eV (HSE06) above the valence band maximum, with $^{2}F_{5/2}$ located 1.26 eV higher, and then compare defect transition levels against this optical window. They find that under fluorine-rich growth the formation of $\mathrm{Yb_i}$ is strongly suppressed, while the aggregation of $\mathrm{Yb_{Ca}}$ with interstitial fluorine into $\mathrm{Yb_{Ca}{-}F_i}$ and its dimer is exothermic, with binding energies of 0.68 eV and 1.93 eV for the monomer and the $\mathrm{Yb_i{-}F_i}$ pair respectively. The conclusion is that fluorine-rich conditions simultaneously remove the quenchers and promote the experimentally observed Yb clustering, providing a growth-chemistry guideline for efficient up/down-conversion.
Load-bearing premise
The placement of the Yb3+ ground and optically excited levels ($^{2}F_{7/2}$ and $^{2}F_{5/2}$) inside the DFT band gap is inferred from the Fermi energy of defective supercells rather than from a many-body treatment of the 4f manifold; if those level positions are wrong, the designation of $\mathrm{Yb_{Ca}}$ and $\mathrm{Yb_i{-}F_i}$ as quenchers and the fluorine-rich growth recommendation would collapse.
Editorial extensions
If this is right
- Fluorine-rich (calcium-poor) synthesis should be adopted for Yb:CaF2 up- and down-conversion phosphors, since it suppresses both identified quencher defects while leaving the beneficial Yb clusters intact.
- The same formation-energy and transition-level machinery can be applied to other trivalent lanthanide dopants in CaF2 to predict which growth conditions minimize trapping and maximize clustering.
- The binding-energy analysis implies that n-type doping should favor $\mathrm{Yb_{Ca}{-}2F_i}$ while p-type doping should favor $\mathrm{Yb_{Ca}{-}F_i}$, offering a Fermi-level handle for cluster control beyond growth chemistry.
- If the quencher assignment is correct, minimizing $\mathrm{Yb_{Ca}}$ and $\mathrm{Yb_i{-}F_i}$ should directly show up as increased luminescence quantum efficiency and reduced non-radiative decay in Yb:CaF2 materials.
- The positive binding energies for dimerization rationalize why Yb clustering persists even at low doping concentrations, matching long-standing experimental observations.
Reading between the lines
- An extension not made in the paper is that the same fluorine-rich logic likely applies to other trivalent lanthanides in CaF2, since the charge-compensating interstitial fluorine mechanism is not unique to ytterbium.
- The Fermi-energy placement of the Yb3+ 4f levels is a rough single-particle proxy; a many-body treatment of the 4f manifold could shift the trap-level assignments and change which defects are classified as quenchers.
- A testable next step would be to compute carrier capture cross-sections or non-radiative recombination rates from the calculated trap levels, converting level positions into quantitative lifetime and efficiency predictions.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a first-principles defect study of Yb-doped CaF2 using PBEsol and hybrid HSE06 (with 50% Hartree-Fock exchange) total-energy calculations. The authors compute thermodynamic transition levels, formation energies, and binding energies for point defects (YbCa, Ybi, Yi, Fi, VCa, VF, Cai) and pair defects (Ybi-Fi, YbCa-Fi, YbCa-2Fi, 2(YbCa-Fi)). They place the Yb3+ 2F7/2 ground level at the Fermi energy of supercells containing nominally charged Yb3+ (1.03 eV above VBM in PBEsol, 2.76 eV in HSE06) and the 2F5/2 excited level 1.26 eV above it. Based on the resulting yellow-quenching window in Fig. 2, they identify YbCa and Ybi-Fi as bulk quenchers with deep trap levels, and recommend fluorine-rich growth to suppress these defects and promote Yb clustering, supported by positive binding energies.
Significance. If the central claim holds, the paper provides a practical and quantitative guideline for controlling defect clustering and luminescence efficiency in Yb:CaF2, a widely used up/down-conversion material. The formation-energy formalism, chemical-potential constraints, and finite-size corrections follow established practice (Freysoldt et al.), and the positive binding energies for YbCa-Fi and 2(YbCa-Fi) are consistent with the experimentally observed clustering tendency. The paper also proposes n-/p-type doping as a clustering-control strategy. However, the central quencher assignment and the F-rich recommendation rest on the placement of the Yb3+ 4f levels, which is inferred from a Kohn-Sham Fermi-energy proxy rather than from total-energy transition levels, and on a nonstandard 50% Hartree-Fock admixture that is not validated for 4f states. These issues make the main conclusion conditional rather than fully established.
major comments (3)
- [Section 2 and Fig. 2] The placement of the Yb3+ ground state (2F7/2) at 1.03 eV (PBEsol) and 2.76 eV (HSE06) above the VBM is set equal to the Fermi energy of defective supercells containing a nominally charged Yb3+ ion. This is a Kohn-Sham eigenvalue proxy, not a total-energy transition level, whereas the defect transition levels in the same figure are computed from Eq. (1) using total-energy differences. For a localized 4f state, these two quantities need not coincide, especially at the nonstandard 50% Hartree-Fock mixing used here. The paper itself concedes in the Summary that the 2F7/2 and 2F5/2 levels were "presumably determined." Because a shift of roughly 0.3 eV in this anchor can move YbCa's ε(1+/0) and Ybi-Fi's ε(2+/1+) in or out of the quenching window, the identification of these defects as quenchers and the resulting F-rich growth recommendation are not robust. I recommend computing the 4f-related level from total-energy differences (e.g., constrained DFT, ΔSCF, or a many-body method) or demonstrating that the quencher assignment is insensitive to the anchor position.
- [Computational methods and Fig. 2] The HSE06 functional is used with 50% nonlocal Hartree-Fock exchange, which is not the standard 25% mixing and is a significant departure from established practice. The paper justifies this choice by the underestimated PBEsol gap and the self-interaction error of Yb 4f orbitals, but no validation is provided for how this admixture affects the absolute positions of the Yb 4f levels or the defect transition energies. The 50% mixing also affects the Fermi-energy anchor (2.76 eV vs 1.03 eV), making the position of the yellow window in Fig. 2 functionally dependent on an unvalidated functional parameter. I request a sensitivity test (e.g., standard HSE06 with 25% exchange, or HSE06 with a Hubbard U correction) to show that the central conclusion is not an artifact of this choice.
- [Section "Fermi energy" and Table S2] The key numerical inputs for the Yb3+ level anchor, namely the Fermi energies 1.03 eV and 2.76 eV, are stated in the text but are documented only in the supplementary Table S2, which is not available in the arXiv version. The construction of the "defective supercells containing nominally charged Yb3+" is not described, and the Fermi energy is not a well-defined output of a finite supercell calculation unless a specific definition (e.g., the highest occupied Kohn-Sham eigenvalue or the electron chemical potential) is given. These values are load-bearing for the central claim, so they should be reported in the main text with a clear definition and the underlying data should be available for independent checking.
minor comments (3)
- [References] Reference [26] (Perdew et al.) is missing the journal name; it should be Phys. Rev. Lett. 100, 136406 (2008).
- [Abstract, final body paragraph] The abstract says "n- or p-type doping" while the body uses "p-/n-type doping" (final paragraph of the binding-energy section); please use consistent ordering throughout.
- [Fig. 2 caption] The yellow regions in Fig. 2 are described as representing the Fermi energy values between the 2F7/2 and 2F5/2 levels; for clarity, state explicitly that this range is the assumed quenching window, not a computed Fermi level of a particular sample.
Circularity Check
No significant circularity: the defect transition levels, formation energies, and binding energies are computed from DFT total-energy differences, and the Yb3+ level placement is an interpretive Fermi-energy proxy rather than a fitted input.
full rationale
The paper's central quantitative claims—transition levels via Eq. (1), formation energies via the standard chemical-potential formalism, and binding energies as Hf[A]+Hf[B]−Hf[AB]—are all total-energy differences obtained from DFT supercells, not quantities fitted to the conclusions they support. The quencher assignment for YbCa and Ybi-Fi is made by comparing these total-energy transition levels with a yellow '2F7/2–2F5/2' window in Fig. 2, where the ground level is approximated by the Fermi energy of a nominally Yb3+-charged supercell (1.03 eV and 2.76 eV above VBM, Table S2) and the excited level is placed 1.26 eV higher. That anchor is a physical approximation and the paper itself labels the multielectronic levels as 'presumably determined'; this is a robustness/correctness concern, not a circular one, because the transition levels are computed independently from total-energy differences rather than chosen to reproduce the quenching verdict. The only self-citation (ref [31], for potential alignment) refers to a standard finite-size correction method and is not the load-bearing prediction; it is corroborated by the Freysoldt review (ref [30]). No equation in the paper reduces a predicted outcome to a fitted input, and no uniqueness theorem or ansatz is imported from the authors' prior work to force the central conclusion. The binding energies also independently reproduce the experimentally known Yb-clustering tendency, consistent with an externally grounded calculation. Overall circularity burden is minimal.
Assumptions & free parameters
free parameters (1)
- Hartree-Fock exchange fraction in HSE06 =
50% (standard HSE06 uses 25%)
assumptions (5)
- domain assumption Kohn-Sham DFT with PBEsol and HSE06 is accurate enough to rank defect formation energies and binding energies quantitatively.
- standard math The defect formation energy and finite-size correction formalism of Freysoldt et al. applies to these supercells.
- domain assumption Chemical potential bounds under F-rich and F-poor conditions are set by elemental bulk and gas references plus the YbF3 secondary-phase constraint.
- ad hoc to paper The ground and optically excited states of Yb3+ can be placed in the DFT band structure using the Fermi energy of a supercell containing nominally charged Yb3+.
- domain assumption Treating Yb 4f and 6s electrons as valence electrons with scalar relativistic PAW pseudopotentials is sufficient.
Cite this review
Pith. "Pith review of Defect physics in $Yb^{3+}$-doped $CaF_2$ from first-principles calculation." pith.science (2026). https://pith.science/paper/GIM32IRK
@misc{pith2026190902397,
author = {Pith},
title = {Pith review of: Defect physics in $Yb^3+$-doped $CaF_2$ from first-principles calculation},
year = {2026},
howpublished = {\url{https://pith.science/paper/GIM32IRK}},
note = {Machine review of arXiv:1909.02397}
}
abstract
Calcium fluoride has been widely used for light up-/down-conversion luminescence by accommodating lanthanide ions as sensitizers or activators. Especially, Yb-doped \ce{CaF2} exhibits unique defect physics, causing various effects on the luminescence. This makes it vital for high efficiency of devices to control the defect-clustering, but theoretically principal guidelines for this are rarely provided. Here we perform the first-principles study on defect physics in Yb-doped \ce{CaF2} to reveal the thermodynamic transition levels and formation energies of possible defects. We suggest that the fluorine rich growth condition can play a key role in enhancing the luminescence efficiency by facilitating the Yb-clustering and suppressing the defect quenchers in bulk. Detailed energetics of defect aggregation not only well explains the experimentally favored Yb-clustering but also presents $n$- or $p$-type doping method for the cluster control.
Figures
Reference graph
Works this paper leans on
- [1]
-
[2]
A. Bagheri, H. Arandiyan, C. Boyer, M. Lim, Lanthanide-doped upcon- version nanoparticles: emerging intelligent light-activated drug delivery systems, Adv. Sci. 3 (2016) 1500437
work page 2016
-
[3]
X. Ai, C. J. H. Ho, J. Aw, A. B. E. Attia, J. Mu, Y . Wang, X. Wang, Y . Wang, X. Liu, H. Chen, M. Gao, X. Chen, E. K. L. Yeow, G. Liu, M. Olivo, B. Xing, In vivo covalent cross-linking of photon-converted rare-earth nanostructures for tumour localization and theranostics, Nat. Commun. 7 (2016) 10432
work page 2016
-
[4]
Y .-T. Wong, S.-Y . Pang, M.-K. Tsang, Y . Liu, H. Huang, S.-F. Yu, J. Hao, Electrochemically assisted flexible lanthanide upconversion lu- minescence sensing of heavy metal contamination with high sensitivity and selectivity, Nanoscale Adv. 1 (2019) 265–272
work page 2019
- [5]
-
[6]
B. Zhou, L. Tao, Y . H. Tsang, W. Jin, Coreshell nanoarchitecture: a strategy to significantly enhance white-light upconversion of lanthanide- doped nanoparticles, J. Mater. Chem. C 1 (2013) 4313–4318
work page 2013
-
[7]
Y . Li, Y . Zhou, X. Li, J. Sun, Z. Ren, W. Wen, X. Yang, G. Hana, A facile approach to upconversion crystalline CaF2:Yb 3+,Tm3+mSiO2 nanospheres for tumor therapy, RSC Adv. 6 (2016) 38365
work page 2016
-
[8]
K. Prorok, A. Bednarkiewicz, B. Cichy, A. Gnach, M. Misiak, M. Sobczykd, W. Strek, The impact of shell host (NaYF4/CaF2) and shell deposition methods on the up-conversion enhancement in Tb 3+, Yb 3+ codoped colloidalα-NaYF4 core-shell nanoparticles, Nanoscale 6 (2014) 1855–1864
work page 2014
Show all 36 references
-
[9]
Corish, C
J. Corish, C. R. A. Catlow, P. W. M. Jacobs, S. H. Ong, Defect aggregation in anion-excess fluorites. Dopant monomers and dimers, Phys. Rev. B 25 (1982) 6425–6438
1982
-
[10]
R. H. Petit, P. Evesque, J. Duran, Dimers and clusters in CaF2:Pr3+. Laser selective excitation and time-resolved spectroscopy, J. Phys. C: Solid State Phys. 14 (1981) 5081–5090
1981
-
[11]
Petit, P
V . Petit, P. Camy, J.-L. Doualan, X. Portier, R. Moncorg ´e, Spectroscopy of Yb3+:CaF2: From isolated centers to clusters, Phys. Rev. B 78 (2008) 085131
2008
-
[12]
Lyberis, A
A. Lyberis, A. J. Stevenson, A. Suganuma, S. Ricaud, F. Druon, F. Herbst, D. Vivien, P. Gredin, M. Mortier, Effect of Yb3+ concentration on optical properties of Yb:CaF2 transparent ceramics, Optical Materials 34 (2012) 965–968
2012
-
[13]
Aidilibike, J
T. Aidilibike, J. Guo, L. Wang, X. Liu, Y . Lia, W. Qin, Ultraviolet upcon- version emission of Pb 2+ ions sensitized by Yb 3+-trimers in CaF 2, RSC Adv. 7 (2017) 2676–2681
2017
-
[14]
Y . Li, J. Guo, X. Liu, T. Aidilibike, W. Qin, White upconversion lumi- nescence in CaF2: Yb3+/Eu3+ powders via incorporating Y3+ ions, Phys. Chem. Chem. Phys. 18 (2016) 16094–16097
2016
-
[15]
Y . Li, J. Guo, X. Liu, T. Aidilibike, W. Qin, Pure red upconversion emis- sion from CaF2:Yb3+/Eu3+, J. Lumin 185 (2017) 247–250
2017
-
[16]
Dong, L.-D
H. Dong, L.-D. Sun, C.-H. Yan, Basic understanding of the lanthanide related upconversion emissions, Nanoscale 5 (2013) 5703–5714
2013
-
[17]
Molina, V
P. Molina, V . Vasyliev, E. G. V´ıllora, K. Shimamura, Tb3+-Yb3+ cooper- ative down and up conversion processes in Tb0.81Ca0.19F2.81:Yb3+ single crystals, J. Appl. Phys. 110 (2011) 123527
2011
-
[18]
J. T. van Wijngaarden, S. Scheidelaar, T. J. H. Vlugt, M. F. Reid, A. Mei- jerink, Energy transfer mechanism for downconversion in the (Pr 3+, Yb3+) couple, Phys. Rev. B 81 (2010) 155112
2010
-
[19]
Ricaud, F
S. Ricaud, F. Druon, D. N. Papadopoulos, A. Pellegrina, F. Balembois, P. Georges, A. Courjaud, P. Camy, J. L. Doualan, R. Moncorg ´e, High- power diode-pumped cryogenically cooled Yb:CaF2 laser with extremely low quantum defects, Opt. Lett. 36 (2011) 1602–1604
2011
-
[20]
Serrano, A
D. Serrano, A. Braud, J. L. Doualan, P. Camy, R. Moncorg ´e, Pr3+ cluster management in CaF 2 by codoping with Lu 3+ or Yb3+ for visible lasers and quantum down-converters, J. Opt. Soc. Am. B 29 (2012) 1854–1862
2012
-
[21]
G. Chen, C. Yang, P. N. Prasad, Nanophotonics and Nanochemistry: Controlling the Excitation Dynamics for Frequency Up- and Down- Conversion in Lanthanide-Doped Nanoparticles, Acc. Chem. Res. 46 (2013) 1474–1486. 4
2013
-
[22]
C. Chen, C. Li, Z. Shi, Current Advances in Lanthanide-Doped Upconver- sion Nanostructures for Detection and Bioapplication, Adv. Sci. 3 (2016) 1600029
2016
-
[23]
Serrano, A
D. Serrano, A. Braud, J.-L. Doualan, P. Camy, R. Moncorg ´e, Highly ef- ficient energy transer in Pr 3+, Yb3+ codoped CaF2 for luminescent solar converters, J. Opt. Soc. Am. B 28 (2011) 1760–1765
2011
-
[24]
P. J. Bendall, C. R. A. Catlow, J. Corish, P. W. M. Jacobs, Defect ag- gregation in anion-excess fluorites II. Clusters containing more than two impurity atoms., J. Solid State Chem. 51 (1984) 159
1984
-
[25]
Giannozzi and S
P. Giannozzi and S. Baroni and N. Bonini and M. Calandra and R. Car, et al., QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials, J. Phys.:Condens. Matter 21 (2009) 395502
2009
-
[26]
J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, K. Burke 100 (2008) 136406
2008
-
[27]
J. Heyd, G. E. Scuseria, M. Ernzerhof, Erratum: Hybrid Functionals Based on a Screened Coulomb Potential [J. Chem. Phys. 118, 8207 (2003)], J. Chem. Phys. 124 (2006) 219906
2003
-
[28]
J. L. F. Da Silva, M. V . Ganduglia-Pirovano, J. Sauer, V . Bayer, G. Kresse, Formation of the cerium orthovanadate CeVO 4: DFT +U study, Phys. Rev. B 75 (2007) 045121
2007
-
[29]
B. G. Janesko, T. M. Henderson, G. E. Scuseria, Screened hybrid func- tionals for solid-state chemistry and physics, Phys. Chem. Chem. Phys. 11 (2009) 443
2009
-
[30]
Freysoldt, B
C. Freysoldt, B. Grabowski, T. Hickel, J. Neugebauer, G. Kresse, A. Jan- otti, C. G. Van de Walle, First-principles calculations for point defects in solids, Rev. Mod. Phys. 86 (2014) 253–305
2014
-
[31]
Kye, C.-J
Y .-H. Kye, C.-J. Yu, U.-G. Jong, K.-C. Ri, J.-S. Kim, S.-H. Choe, S.- N. Hong, S. Li, J. N. Wilson, A. Walsh, Vacancy-Driven Stabilization of the Cubic Perovskite Polymorph of CsPbI3, J. Phys. Chem. C 123 (2019) 9735–9744
2019
-
[32]
D. N. Batchelder, R. O. Simmons, Lattice constants and thermal expan- sivities of silicon and calcium fluoride between 6◦C and 322◦C, J. Chem. Phys. 41 (1964) 2324–2329
1964
-
[33]
G. W. Rublo ff, Far-ultraviolet reflectance spectra and the electronic struc- ture of ionic crystals, Phys. Rev. B. 5 (1972) 662–684
1972
-
[34]
M. L. Agiorgousis, Y .-Y . Sun, H. Zeng, S. Zhang, Strong Covalency- Induced Recombination Centers in Perovskite Solar Cell Material CH3NH3PbI3, J. Am. Chem. Soc. 136 (2014) 14570–14575
2014
-
[35]
A. M. Ibraheem, M. A. Khalafalla, M. H. Eisa, First principle calculation of accurate native defect levels in CaF2, Eur. Phys. J. B 90 (2017) 40
2017
-
[36]
Du, Using DFT Methods to Study Activators in Optical Materials, ECS J
M.-H. Du, Using DFT Methods to Study Activators in Optical Materials, ECS J. Solid State Sc. 5(1) (2016) 30073018. 5
2016
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.