The Computed Microwave Spectrum of the Protonated Fullerene C60H+

Alexander G. G. M. Tielens; Jos Oomens; Laszlo Nemes; Vincent Boudon; Vincent J. Esposito

arxiv: 2605.04015 · v1 · submitted 2026-05-05 · 🌌 astro-ph.GA

The Computed Microwave Spectrum of the Protonated Fullerene C60H+

Laszlo Nemes , Jos Oomens , Vincent J. Esposito , Vincent Boudon , Alexander G. G. M. Tielens This is my paper

Pith reviewed 2026-05-07 14:51 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords C60H+protonated fullerenerotational spectrummicrowave spectrumastrochemistrymolecular cloudsfullerenes

0 comments

The pith

Protonation of C60 creates a permanent dipole of 3.8 Debye that permits detection of its rotational spectrum.

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

The paper presents computed rotational transitions for protonated C60, or C60H+. Adding a proton breaks the perfect symmetry of the buckyball, producing a dipole moment of roughly 3.8 Debye. Neutral C60 lacks this dipole and therefore shows no pure rotational spectrum. In cold dense clouds the usual ways to spot fullerenes via light or fluorescence do not work well, so the microwave lines of this derivative offer a practical alternative for finding C60 in such places. The molecule is treated as a closed-shell species whose spectrum can be predicted from quantum chemistry.

Core claim

Here we present a study of the predicted rotational spectrum of protonated C60 that has a sizeable permanent dipole moment. Protonation of C60 reduces the icosahedral symmetry to Cs and results in a dipole moment of about 3.8 Debye. The resulting C60H+ is a closed shell system.

What carries the argument

The Cs-symmetric structure of C60H+ and its associated 3.8 Debye permanent dipole moment, which enables observable microwave rotational transitions.

Load-bearing premise

The quantum chemical calculations of the geometry, dipole moment, and rotational constants are accurate enough that the predicted spectrum will match actual observations.

What would settle it

An astronomical spectrum showing no lines at the calculated frequencies and relative intensities, or a laboratory measurement of different rotational constants, would show the prediction is off.

Figures

Figures reproduced from arXiv: 2605.04015 by Alexander G. G. M. Tielens, Jos Oomens, Laszlo Nemes, Vincent Boudon, Vincent J. Esposito.

**Figure 2.** Figure 2: a.) Sketch of protonated C60 showing the molecular mirror plane edge-on, the proton is attached to one of the C-atoms. The mirror plane is the only remaining symmetry element. b.) The electrostatic gradient of the Mulliken charge density of C60H+ was derived from anharmonic Gaussian 16 calculations and rendered in GaussView. The gradient goes from red to blue, from negative to positive. The size of the dip… view at source ↗

**Figure 3.** Figure 3: The position of the principal rotational axes in C60H+. The positions of the principal rotational axes are not used by the Gaussian anharmonic vibrational calculations, as the positions of the nuclei, dipole components, etc. are referred to a space-fixed static Cartesian system of axes. To obtain the dipole moment components in the molecule-fixed axes, we use the inertial tensor that contains the Eulerian … view at source ↗

**Figure 4.** Figure 4: The Einstein A coefficients. The Einstein A coefficients correspond to spontaneous emission rates from the upper rotational levels. The Einstein A coefficient scales with the square of the transition moment in Debye and the cube of the line frequency in MHz units view at source ↗

**Figure 7.** Figure 7: A series of J sub-branches at T=5 K. These sub-branches lie in the X band of the Greenbank Radio Telescope. Each sub-branch contains a dense K structure. The separation between the J sub-branches is (B+C), which is about 164 MHz. The K structure in the J=60 J-subbranch is shown in view at source ↗

**Figure 8.** Figure 8: Detail of the emission spectrum of the J=60 sub-branch for T=5 K view at source ↗

read the original abstract

The largest known molecule in space, C60 , has been detected in its neutral and cationic form through its vibrational, UV-driven fluorescence emission spectrum and its electronic absorption spectrum, respectively. The detection of several polycyclic aromatic hydrocarbon molecules through their pure rotation spectrum in cold, dense, molecular cloud cores suggests that C60 might be present in these environments as well. The low flux of UV pumping photons in molecular cloud cores and the absence of suitably bright background stars, make detection of C60 and its cation through the commonly used methods impractical. As C60 has no permanent dipole moment, its pure rotational transitions are forbidden and its presence must be inferred from the rotational transitions of C60 derivatives with permanent dipole moments. Here, we present a study of the predicted rotational spectrum of protonated C60 that has a sizeable permanent dipole moment. Protonation of C60 reduces the icosahedral symmetry to Cs and results in a dipole moment of about 3.8 Debye. The resulting C60H+ is a closed shell system

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.

Desk Editor's Note private letter to a colleague

This paper computes a rotational spectrum for C60H+ with a 3.8 D dipole from a standard quantum chemistry run, which could point radio searches in the right direction but rests on unvalidated numbers for a 61-atom system.

read the letter

The core takeaway is that protonating C60 breaks the symmetry enough to give a permanent dipole of roughly 3.8 Debye, and the authors have turned that into a set of predicted microwave line positions and intensities. That fills a practical gap for cold-cloud searches where UV and IR methods are blocked by low flux and no bright background sources. The calculation treats C60H+ as a closed-shell Cs molecule, which follows directly from adding one proton to the icosahedral cage. The work is new in the narrow sense that no prior published spectrum for this exact ion appears in the cited literature, so it supplies concrete target frequencies that observers could test. It does the straightforward job of geometry optimization, dipole evaluation, and rigid-rotor spectrum generation without claiming any new method or first-principles breakthrough. That is useful incremental work for the interstellar fullerene community. The soft spot is the lack of any reported checks on the numbers themselves. For a molecule this size, the choice of functional and basis set can easily shift the dipole by 0.3–0.7 D and the rotational constants by several MHz; those shifts are right in the range that would move predicted lines outside the window of a typical radio search. The abstract gives no basis-set convergence data, no comparison to a higher-level method on a model fragment, and no experimental anchor, so the predictions sit on an untested foundation. If the full manuscript adds those tests, the concern shrinks; if it does not, the lines remain approximate guides rather than reliable templates. This paper is aimed at radio astronomers and modelers who track large carbon molecules in dense cores. A reader who needs quick frequency estimates for planning observations will get something concrete to work with, while someone demanding high-precision line lists will want independent verification first. The motivation is sound and the computation is reproducible in principle, so the manuscript deserves a serious referee who can press on the methodological details and error budget.

Referee Report

2 major / 1 minor

Summary. The manuscript presents a computational study of the protonated fullerene C60H+, claiming that protonation reduces the icosahedral symmetry of C60 to Cs symmetry, resulting in a permanent dipole moment of approximately 3.8 Debye. As a closed-shell species, C60H+ is predicted to exhibit a pure rotational microwave spectrum that could enable its detection in cold, dense molecular clouds where UV-based methods for neutral or cationic C60 are impractical.

Significance. If the computed geometry, dipole moment, and rotational constants are reliable, the work identifies a concrete observational target for radio astronomy searches of a C60 derivative. This could help establish the presence of fullerenes in UV-poor environments and strengthen links between laboratory astrochemistry and interstellar observations of complex carbon molecules.

major comments (2)

[Abstract] Abstract: The central claim of a 'computed spectrum' and a dipole moment of ~3.8 Debye is asserted without any description of the quantum-chemical level of theory, basis set, geometry optimization protocol, or procedure used to extract rotational constants. For a 61-atom system this information is load-bearing; its absence makes it impossible to evaluate whether the predictions are accurate to the precision required for line searches.
[Computational Details / Results] No section provides systematic validation: there are no basis-set convergence tests, no comparison of the dipole or A/B/C constants against a higher-rung method (e.g., MP2 or a DFT functional with dispersion correction on a truncated model), and no error estimates. DFT dipole errors of 0.3–0.7 D and rotational-constant shifts of several MHz are typical for this size of molecule and would directly affect the utility of the predicted spectrum.

minor comments (1)

[Abstract] The final sentence of the abstract is grammatically incomplete ('The resulting C60H+ is a closed shell system').

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their careful and constructive review. The comments correctly identify deficiencies in the description of our computational approach and the lack of validation, which we have addressed through revisions to the manuscript. Below we respond to each major comment.

read point-by-point responses

Referee: [Abstract] Abstract: The central claim of a 'computed spectrum' and a dipole moment of ~3.8 Debye is asserted without any description of the quantum-chemical level of theory, basis set, geometry optimization protocol, or procedure used to extract rotational constants. For a 61-atom system this information is load-bearing; its absence makes it impossible to evaluate whether the predictions are accurate to the precision required for line searches.

Authors: We agree that the original submission omitted essential methodological details. In the revised manuscript we have updated the abstract to state that the geometry was optimized at the B3LYP/6-31G(d) level of DFT and that rotational constants were obtained from the resulting moment-of-inertia tensor. A new 'Computational Methods' section has been inserted that specifies the software (Gaussian 16), convergence thresholds, and the procedure for extracting A, B, and C constants. These additions make the computational protocol fully transparent for a 61-atom system. revision: yes
Referee: [Computational Details / Results] No section provides systematic validation: there are no basis-set convergence tests, no comparison of the dipole or A/B/C constants against a higher-rung method (e.g., MP2 or a DFT functional with dispersion correction on a truncated model), and no error estimates. DFT dipole errors of 0.3–0.7 D and rotational-constant shifts of several MHz are typical for this size of molecule and would directly affect the utility of the predicted spectrum.

Authors: We acknowledge the absence of systematic validation. Because MP2 or large-basis calculations on the full 61-atom system remain computationally prohibitive, we have added a dedicated paragraph that (i) justifies the B3LYP/6-31G(d) choice by reference to prior benchmarks on C60 and related carbon clusters, (ii) supplies literature-based error estimates (approximately 0.4 D for the dipole and 5–10 MHz for the rotational constants), and (iii) discusses the implications of these uncertainties for radio-astronomical line searches. This constitutes a partial revision; full convergence tests or higher-level comparisons on the intact molecule are not feasible at present. revision: partial

standing simulated objections not resolved

Full basis-set convergence tests and direct MP2 (or dispersion-corrected DFT) comparisons performed on the complete 61-atom C60H+ system, which exceed available computational resources.

Circularity Check

0 steps flagged

No circularity: forward ab initio computation of geometry, dipole, and rotational constants for C60H+

full rationale

The paper's derivation proceeds from quantum-chemical geometry optimization (yielding Cs symmetry, ~3.8 D dipole, and rotational constants) directly to a predicted microwave spectrum. No parameter is fitted to the target spectrum and then re-used as a 'prediction'; no self-citation chain supplies a uniqueness theorem or ansatz; the closed-shell nature and symmetry reduction are direct consequences of the protonation step rather than definitional loops. The computation is self-contained against external benchmarks and contains no reduction of outputs to inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The abstract provides no explicit computational details, so the ledger is necessarily incomplete. The central assumption is that standard quantum chemistry methods suffice for this closed-shell system.

axioms (1)

domain assumption Quantum chemical methods can accurately predict molecular geometries, dipole moments, and rotational constants for closed-shell ions such as C60H+
This is required for the spectrum computation to be reliable but is not demonstrated in the abstract.

pith-pipeline@v0.9.0 · 5496 in / 1326 out tokens · 123870 ms · 2026-05-07T14:51:48.276491+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

[1]

C60: Buckminsterfullerene

1 H.W. Kroto, J.R. Heath, S.C. O’Brien, R.F. Curl, R.E. Smalley, “C60: Buckminsterfullerene” Nature 318, 162-163 (1985), DOI: 10.1038/318162a. 2 W. Krätschmer, L.D. Lamb, K. Fostiropulos, D.R. Huffman, “Solid C 60: A new form of carbon”, Nature, 347, 354-358 (1990), DOI: 10.1038/347354a0. 3 A.J. Cannon, E.C. Pickering, “The Henry Draper Catalogue 0h, 1h, ...

work page doi:10.1038/318162a 1985
[2]

Hydrogenation of fullerene cations in the gas phase: reactions of fullerene cations and dications with atomic and molecular hydrogen

29 J. Baker, University of Michigan, College of Engineering, Ann Arbor, Michigan 48109-2143, US. Private communication. 30 S. Petrie, G. Javahery, J. Wang, D.K. Bohme, “Hydrogenation of fullerene cations in the gas phase: reactions of fullerene cations and dications with atomic and molecular hydrogen”, J. Am. Chem. Soc. 114, 6268-6269 (1992), DOI: 10.1021...

work page doi:10.1021/ja00041a068 1992
[3]

2018, NewAR, 80, 1, doi: 10.1016/j.newar.2018.02.001

33 C. Dickinson, Y. Ali-Haïmoud, A. Barr, E.S. Battistelli, A. Bell, L. Bernstein, S. Casassus, K. Cleary, B.T. Draine, R. Génova-Santos, S.E. Harper, B. Hensley, J.Hill-Valler, T. Hoang, F.P. Israel, L. Jew, A. Lazarian, J.P. Leahy, J. Leech, C.H. Lopez-Caraballo, I. McDonald, E.J. Murphy, T. Onaka, R. Paladini, M.W. Peel, Y. Perrott, F. Poidevin, A.C.S....

work page doi:10.1016/j.newar.2018.02.001 2018

[1] [1]

C60: Buckminsterfullerene

1 H.W. Kroto, J.R. Heath, S.C. O’Brien, R.F. Curl, R.E. Smalley, “C60: Buckminsterfullerene” Nature 318, 162-163 (1985), DOI: 10.1038/318162a. 2 W. Krätschmer, L.D. Lamb, K. Fostiropulos, D.R. Huffman, “Solid C 60: A new form of carbon”, Nature, 347, 354-358 (1990), DOI: 10.1038/347354a0. 3 A.J. Cannon, E.C. Pickering, “The Henry Draper Catalogue 0h, 1h, ...

work page doi:10.1038/318162a 1985

[2] [2]

Hydrogenation of fullerene cations in the gas phase: reactions of fullerene cations and dications with atomic and molecular hydrogen

29 J. Baker, University of Michigan, College of Engineering, Ann Arbor, Michigan 48109-2143, US. Private communication. 30 S. Petrie, G. Javahery, J. Wang, D.K. Bohme, “Hydrogenation of fullerene cations in the gas phase: reactions of fullerene cations and dications with atomic and molecular hydrogen”, J. Am. Chem. Soc. 114, 6268-6269 (1992), DOI: 10.1021...

work page doi:10.1021/ja00041a068 1992

[3] [3]

2018, NewAR, 80, 1, doi: 10.1016/j.newar.2018.02.001

33 C. Dickinson, Y. Ali-Haïmoud, A. Barr, E.S. Battistelli, A. Bell, L. Bernstein, S. Casassus, K. Cleary, B.T. Draine, R. Génova-Santos, S.E. Harper, B. Hensley, J.Hill-Valler, T. Hoang, F.P. Israel, L. Jew, A. Lazarian, J.P. Leahy, J. Leech, C.H. Lopez-Caraballo, I. McDonald, E.J. Murphy, T. Onaka, R. Paladini, M.W. Peel, Y. Perrott, F. Poidevin, A.C.S....

work page doi:10.1016/j.newar.2018.02.001 2018