REVIEW 3 major objections 4 minor 43 references
C$_{60}$ building blocks with tuneable structures for tailored functionalities
T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Four crystal phases built from C60 molecules have tuneable band gaps from 2.03 to 2.46 eV and bright-exciton binding energies above 0.2 eV, making fullerenes a versatile semiconductor platform.
desk verdict Useful new numbers on C60 polymorphs, but two of the four phases are unstable as computed, so the band-gap tuning story needs to be scaled back. 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 central object is the C60 molecule treated as an 'ultimate building block,' whose icosahedral symmetry and orientational degrees of freedom generate distinct crystal symmetries when placed on different lattices. The key computational machinery is the combination of PBE0 hybrid-functional band-structure calculations with time-dependent Hartree-Fock (TDHF) exciton calculations, which the authors argue reproduces GW+BSE-quality results at lower computational cost. The mechanism carrying the argument is the variation in intermolecular bonding (covalent [2+2] cycloaddition in the layered phases versus van der Waals interactions in the cubic phases) and in molecular orientation, which together shift band gaps by roughly 0.4 eV and change the character of the gap from indirect to direct.
What would settle it
A decisive test is an anharmonic phonon or finite-temperature molecular dynamics calculation on the Fm-3 and Pa-3 phases: if the imaginary rotational modes persist at all temperatures up to the 255 K phase transition, the predicted 2.293 eV and 2.461 eV direct gaps describe unrealizable structures. A complementary experiment would be phase-pure optical absorption measurements on each polymorph to look for the predicted sub-gap exciton peaks near 1.9 eV and direct-gap onsets near 2.3–2.5 eV.
Extended reading notes
Core claim
The paper's central discovery is that orientational and lattice degrees of freedom of C60 molecules provide a structural knob for engineering electronic and optical properties. In the orthorhombic Immm phase, C60 cages are linked by [2+2] cycloaddition bonds into a quasi-square 2D lattice with an indirect PBE0 band gap of 2.075 eV; in the trigonal R-3m phase, the same bonding motif forms a quasi-triangular lattice with an indirect gap of 2.027 eV. In contrast, the face-centred cubic Fm-3 phase, in which all molecules share one orientation, has a direct gap of 2.293 eV at X, and the Pa-3 phase with four inequivalent molecular orientations has a direct gap of 2.461 eV at R. The calculated optical absorption shows that excitonic effects are strong, with bright-exciton binding energies between 0.206 and 0.295 eV, and that the symmetry of the lattice controls the polarization dependence of absorption: R-3m absorbs strongly in-plane across the visible spectrum, while Immm has stronger out-of-plane exciton absorption below the gap.
Load-bearing premise
The load-bearing premise is that the four idealized static crystal structures used in the calculations are physically realizable phases; the paper's own phonon spectra show imaginary rotational modes in both cubic phases, so if those structures cannot be stabilized, the predicted direct gaps and optical spectra are properties of artificial crystals.
Editorial extensions
If this is right
- The four fullerene polymorphs provide a single-molecule materials platform with band gaps spanning 2.027–2.461 eV, covering both indirect (layered) and direct (cubic) optical transitions.
- The layered Immm and R-3m phases exhibit polarization-dependent excitonic absorption, with R-3m absorbing strongly in-plane across the visible spectrum, enabling polarisation filters and directional photodetectors.
- Bright-exciton binding energies of 0.206–0.295 eV are large enough that excitons should persist near room temperature, relevant for excitonic devices and solar-energy conversion.
- Adjusting the interfullerene distance and molecular orientation tunes effective masses and band-edge properties, offering a route to optimize carrier mobility in fullerene-based semiconductors.
- The two layered phases are dynamically stable at the harmonic level, supporting their synthesis under high-pressure conditions and their use as robust 2D building blocks.
Reading between the lines
- My inference: the paper's design principle—using molecular orientation as a structural degree of freedom—could generalize to other nearly spherical molecular cages (e.g., substituted or endohedral fullerenes), potentially expanding the accessible band-gap range beyond the roughly 0.4 eV span shown here.
- My inference: because the Fm-3 and Pa-3 phases are dynamically unstable at zero temperature, the experimentally observed high-temperature plastic phase of solid C60, which has orientational disorder, would likely have a band gap closer to an average of the static configurations; comparing the predicted direct gaps against measurements on the disordered cubic phase would clarify whether the static
- My inference: a practical testable extension is to compute the electron and hole effective masses at the band edges for each phase, since the paper identifies carrier mobility as a tunable target but does not report effective-mass values.
- My inference: the roughly 1.9 eV exciton peak in the layered phases sits near the solar-spectrum maximum shown in the paper, so these phases may be worth exploring as photovoltaic absorbers, a conclusion the paper does not explicitly draw.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports first-principles calculations for four C60 polymorphs: two layered structures (orthorhombic Immm and trigonal R-3m) and two cubic van der Waals structures (Fm-3 and Pa-3). Using PBE+DFT-D3 for relaxations, harmonic phonons from finite differences, PBE0 for electronic band structures, and PBE0+TDHF for excitonic effects, it claims that these phases form a tunable platform with band gaps between 2.027 and 2.461 eV, bright-exciton binding energies of 0.206-0.295 eV, and distinct optical absorption spectra. The paper also presents phonon spectra, band structures, and absorption curves for all four phases.
Significance. If the computed properties correspond to physically realizable phases, the reported variation of band gap, effective masses, and excitonic response with space group and molecular orientation could be useful for selecting fullerene-based materials for optoelectronics and photovoltaics. The methodological choices (PBE0 for gaps, TDHF for excitons) are reasonable for molecular crystals, and the comparison of the Pa-3 gap with earlier GW calculations (2.461 eV vs. ~2.4 eV) gives some confidence. However, the central claim of a four-phase tunable platform is weakened by the paper's own finding that both cubic phases are dynamically unstable in the idealized static structures used for the electronic-structure calculations; the Fm-3 model is explicitly not the experimentally observed orientationally disordered Fm-3m plastic phase, and the Pa-3 stabilization is only speculated. The significance therefore rests mainly on the two layered phases unless the cubic phases can be shown to be stabilized as modeled.
major comments (3)
- [Section III.B, Fig. 2(c), Fig. 2(d)] The phonon spectra show imaginary rotational modes throughout the Brillouin zone for both Fm-3 and Pa-3, and the text acknowledges that the actual high-temperature phase is the orientationally disordered Fm-3m structure, not the all-aligned Fm-3 model. For Pa-3, the harmonic instability is particularly concerning because the experimentally reported Pa-3 phase is the low-temperature ordered phase (below 255 K); the statement that anharmonic effects at "increased temperatures" will remove imaginary phonons is directionally inconsistent with the known phase transition and is not supported by any calculation. Consequently, the PBE0 band gaps of 2.293 eV (Fm-3) and 2.461 eV (Pa-3) and their optical absorption spectra are properties of saddle-point configurations, not of demonstrated physical phases, so two of the four data points in the proposed band-gap tuning window (2.027-2.461 eV) are not tied to realizable structures. The authors should either demonstrate dynamical stabilization (e.g., via anharmonic free-energy calculations or molecular-dynamics sampling) or clearly reframe the claim to the two layered phases and discuss the cubic structures as hypothetical idealized models.
- [Section III.A and Table I] The text and Table I contradict each other on cohesive-energy ordering. The text states that R-3m "has the highest cohesive energy," but Table I lists R-3m with -7.65026 eV/atom, which is the least negative value and hence the weakest binding under the standard convention used elsewhere ("lower cohesive energies than the two layered phases" for the fcc phases). The fcc phases are said to have lower cohesive energies than the layered phases, yet Pa-3 (-7.66431) and Fm-3 (-7.66258) are more negative than Immm (-7.66107). Also, the text says "The Fm-3 phase has a slightly larger lattice constant and a higher cohesive energy than the Pa-3 phase," while Table I gives Pa-3 a more negative Ec. The definition of cohesive energy (sign convention) must be stated explicitly, and the text and table must be made consistent, since this discussion underlies the relative stability of the layered versus van der Waals phases.
- [Section II and Tables II-III] No convergence tests are reported for any of the central quantitative results. The plane-wave cutoff (800 eV) and k-point grid (3x3x3) are asserted to be "well-converged," but no evidence is given; more importantly, the PBE0 and TDHF calculations depend on the active-space size (16 highest valence bands and 16 lowest conduction bands) and the k-mesh (8x8x8 or 4x4x4), and no convergence with respect to these parameters is shown. Since the paper reports band gaps to 0.001 eV (e.g., 2.075 eV, 2.027 eV) and exciton binding energies to 0.001 eV (e.g., 0.293 eV, 0.295 eV), the reader cannot assess the numerical uncertainty of the claimed tuning window without convergence data.
minor comments (4)
- [Section III.D] The phrase "independent-parcticle approximation" contains a typo; it should be "independent-particle approximation."
- [Section II] The validation of the PBE0+TDHF method cites Refs. [14,17], which are the authors' own papers; including an independent benchmark beyond the already-cited GW+BSE work would strengthen the methodology discussion.
- [Section III.A] The volume per C60 is listed in Table I but not discussed; stating whether these volumes correspond to the experimental lattice constants (especially for the known Pa-3 phase) would help the reader judge the quality of the relaxed geometries.
- [Introduction] The reference list in the Introduction appears to have a formatting issue ("[6 ? ? –12]") that should be corrected.
Circularity Check
No significant circularity: the quantitative results are parameter-free ab initio calculations, and the only self-citations are non-load-bearing method validations corroborated by independent GW+BSE work.
full rationale
The central results (PBE0 gaps of 2.075, 2.027, 2.293 and 2.461 eV, and the PBE0+TDHF exciton binding energies) are direct outputs of first-principles calculations with no fitted input and no target quantity defined in terms of another target quantity. The self-citations relevant to the derivation appear in the Methods section, where PBE0 and PBE0+TDHF are said to agree with many-body GW+BSE for C60 systems by reference to the authors' earlier papers [14,17]; however, the same sentences cite the independent GW+BSE calculation [30] as the benchmark, so the agreement is an external validation rather than an input to the derivation. The paper's own phonon results showing imaginary rotational modes for Fm-3 and Pa-3 (Sec. III.B) and the speculation that anharmonic effects will stabilize Pa-3 are physical-relevance limitations, not circular reductions: they concern whether the static structures correspond to real phases, while the calculated gaps and spectra remain self-contained predictions for those structures and are compared with independent earlier GW and photoemission data. No fitted parameter is relabeled as a prediction, and no load-bearing uniqueness theorem or ansatz is imported solely from the authors' prior work.
Assumptions & free parameters
assumptions (4)
- domain assumption The PBE0 hybrid functional predicts reliable band gaps for C60 molecular crystals when compared with GW+BSE.
- domain assumption The TDHF approach with 16 valence and 16 conduction bands captures the relevant excitonic effects.
- domain assumption DFT-D3 captures the van der Waals interactions between C60 molecules well enough for structural and energetic comparisons.
- domain assumption The harmonic approximation and finite-displacement supercells reliably indicate dynamic stability for the layered phases.
Cite this review
Pith. "Pith review of C$_{60}$ building blocks with tuneable structures for tailored functionalities." pith.science (2026). https://pith.science/paper/6QK6QW5L
@misc{pith2026250101494,
author = {Pith},
title = {Pith review of: C$_60$ building blocks with tuneable structures for tailored functionalities},
year = {2026},
howpublished = {\url{https://pith.science/paper/6QK6QW5L}},
note = {Machine review of arXiv:2501.01494}
}
abstract
We show that C$_{60}$ fullerene molecules can serve as promising building blocks in the construction of versatile crystal structures with unique symmetries using first-principles calculations. These phases include quasi-2D layered structures and 3D van der Waals crystals where the molecules adopt varied orientations. The interplay of molecular arrangement and lattice symmetry results in a variety of tuneable crystal structures with distinct properties. Specifically, the electronic structures of these phases vary significantly, offering potential for fine-tuning the band gap for electronics and optoelectronics. Additionally, the optical properties of these materials are strongly influenced by their crystalline symmetry and molecular alignment, providing avenues for tailoring optical responses for photonics. Our findings highlight the potential of fullerene-based building blocks in the rational design of functional materials.
Figures
Reference graph
Works this paper leans on
-
[1]
F. Jona, X. Z. Ji, and P. M. Marcus, Phys. Rev. B 68, 172101 (2003)
work page 2003
-
[2]
G. Grimvall, B. Magyari-K¨ ope, V. Ozoli¸ nˇ s, and K. A. Persson, Rev. Mod. Phys. 84, 945 (2012)
work page 2012
-
[3]
E. A. Doud, A. Voevodin, T. J. Hochuli, A. M. Champ- saur, C. Nuckolls, and X. Roy, Nature Reviews Materials 5, 371 (2020)
work page 2020
-
[4]
T. Hueckel, G. M. Hocky, and S. Sacanna, Nature Re- views Materials 6, 1053 (2021)
work page 2021
-
[5]
H. W. Kroto, J. R. Heath, S. C. O’Brien, R. F. Curl, and R. E. Smalley, Nature 318, 162 (1985)
work page 1985
-
[6]
P. W. Stephens, L. Mihaly, P. L. Lee, R. L. Whetten, S.-M. Huang, R. Kaner, F. Deiderich, and K. Holczer, Nature 351, 632 (1991)
work page 1991
-
[7]
P. A. Heiney, Journal of Physics and Chemistry of Solids 53, 1333 (1992)
work page 1992
-
[8]
M. N´ u˜ nez-Regueiro, L. Marques, J. L. Hodeau, O. B´ ethoux, and M. Perroux, Phys. Rev. Lett. 74, 278 (1995)
work page 1995
Show all 43 references
-
[9]
C. H. Xu and G. E. Scuseria, Phys. Rev. Lett. 74, 274 (1995)
1995
-
[10]
Springborg, Phys
M. Springborg, Phys. Rev. B 52, 2935 (1995)
1995
-
[11]
Marques, J
L. Marques, J. L. Hodeau, M. N´ u˜ nez Regueiro, and M. Perroux, Phys. Rev. B 54, R12633 (1996)
1996
-
[12]
´Alvarez Murga and J
M. ´Alvarez Murga and J. Hodeau, Carbon82, 381 (2015). 6
2015
-
[13]
L. Hou, X. Cui, B. Guan, S. Wang, R. Li, Y. Liu, D. Zhu, and J. Zheng, Nature 606, 507 (2022)
2022
-
[14]
B. Peng, J. Am. Chem. Soc. 144, 19921 (2022)
2022
-
[15]
Peng, Nano Lett
B. Peng, Nano Lett. 23, 652 (2023)
2023
-
[16]
Meirzadeh, A
E. Meirzadeh, A. M. Evans, M. Rezaee, M. Milich, C. J. Dionne, T. P. Darlington, S. T. Bao, A. K. Bartholomew, T. Handa, D. J. Rizzo, R. A. Wiscons, M. Reza, A. Zan- giabadi, N. Fardian-Melamed, A. C. Crowther, P. J. Schuck, D. N. Basov, X. Zhu, A. Giri, P. E. Hopkins, P. Kim,...
2023
-
[17]
Jones and B
C. Jones and B. Peng, J. Phys. Chem. Lett. 14, 11768 (2023)
2023
-
[18]
Wu and B
J. Wu and B. Peng, J. Am. Chem. Soc. , (2024)
2024
-
[19]
Momma and F
K. Momma and F. Izumi, Journal of Applied Crystallog- raphy 44, 1272 (2011)
2011
-
[20]
Kresse and J
G. Kresse and J. Furthm¨ uller, Phys. Rev. B 54, 11169 (1996)
1996
-
[21]
Kresse and J
G. Kresse and J. Furthm¨ uller, Computational Materials Science 6, 15 (1996)
1996
-
[22]
P. E. Bl¨ ochl, Phys. Rev. B50, 17953 (1994)
1994
-
[23]
Kresse and D
G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999)
1999
-
[24]
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)
1996
-
[25]
Grimme, J
S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, J. Chem. Phys. 132, 154104 (2010)
2010
-
[26]
A. Togo, F. Oba, and I. Tanaka, Phys. Rev. B78, 134106 (2008)
2008
-
[27]
Togo and I
A. Togo and I. Tanaka, Scripta Materialia 108, 1 (2015)
2015
-
[28]
Le Page and P
Y. Le Page and P. Saxe, Phys. Rev. B 65, 104104 (2002)
2002
-
[29]
Adamo and V
C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 (1999)
1999
-
[30]
Champagne, M
A. Champagne, M. Camarasa-G´ omez, F. Ricci, L. Kro- nik, and J. B. Neaton, Nano Lett. 24, 7033 (2024)
2024
-
[31]
Sander and G
T. Sander and G. Kresse, J. Chem. Phys. 146, 064110 (2017)
2017
-
[32]
Iwasa, T
Y. Iwasa, T. Arima, R. M. Fleming, T. Siegrist, O. Zhou, R. C. Haddon, L. J. Rothberg, K. B. Lyons, H. L. Carter, A. F. Hebard, R. Tycko, G. Dabbagh, J. J. Krajewski, G. A. Thomas, and T. Yagi, Science 264, 1570 (1994)
1994
-
[33]
B. Peng, H. Zhang, H. Shao, Z. Ning, Y. Xu, G. Ni, H. Lu, D. W. Zhang, and H. Zhu, Materials Research Letters 5, 399 (2017)
2017
-
[34]
O. I. Malyi, K. V. Sopiha, and C. Persson, ACS Appl. Mater. Interfaces 11, 24876 (2019)
2019
-
[35]
D. Luo, K. Yin, and R. Dronskowski, J. Am. Chem. Soc. 144, 5155 (2022)
2022
-
[36]
Pallikara, P
I. Pallikara, P. Kayastha, J. M. Skelton, and L. D. Whal- ley, Electron. Struct. 4, 033002 (2022)
2022
-
[37]
Mitrano, A
M. Mitrano, A. Cantaluppi, D. Nicoletti, S. Kaiser, A. Perucchi, S. Lupi, P. Di Pietro, D. Pontiroli, M. Ricc` o, S. R. Clark, D. Jaksch, and A. Cavalleri, Nature 530, 461 (2016)
2016
-
[38]
Gueymard, D
C. Gueymard, D. Myers, and K. Emery, Solar Energy 73, 443 (2002)
2002
-
[39]
E. L. Shirley and S. G. Louie, Phys. Rev. Lett. 71, 133 (1993)
1993
-
[40]
Schwedhelm, L
R. Schwedhelm, L. Kipp, A. Dallmeyer, and M. Ski- bowski, Phys. Rev. B 58, 13176 (1998)
1998
-
[41]
T. L. Makarova, Semiconductors 35, 243 (2001)
2001
-
[42]
E. L. Shirley, L. X. Benedict, and S. G. Louie, Phys. Rev. B 54, 10970 (1996)
1996
-
[43]
J. H. Weaver, J. L. Martins, T. Komeda, Y. Chen, T. R. Ohno, G. H. Kroll, N. Troullier, R. E. Haufler, and R. E. Smalley, Phys. Rev. Lett. 66, 1741 (1991)
1991
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.