Completing Chou and Yang's exact solution of the carbon-60 vibration problem: the $Q_{4\pm}$ and $Q_{5\pm}$ characteristic polynomials, normal-mode eigenvectors, and inversion of the measured spectrum to force constants

In 1997 Chou and Yang reduced the harmonic vibration problem of the $C_{60}$ molecule - a generalized eigenvalue problem of order 180 - to irreducible-sector secular determinants of order at most 8, using an Abelian fivefold-rotation subgroup reduction and projection operators, and obtained exact algebraic characteristic polynomials $Q_{\mu p}(\omega^2)$ for the one- and three-dimensional irreducible representations of the icosahedral group as explicit functions of four force constants (single- and double-bond stretching and two angle-bending constants). Three items were explicitly left for later work: the polynomials $Q_{4\pm}$ and $Q_{5\pm}$ of the four- and five-dimensional irreducible representations (the even five-dimensional sector carries eight of the ten Raman-active modes), which were omitted on account of their length; the explicit normal-mode eigenvectors; and the numerical inversion of the measured infrared and Raman spectrum to the four force constants, with its uniqueness, stability and sensitivity questions. The authors' own 2001 sequel, using seven force constants, again does not present the polynomials, gives no eigenvectors, and only computes frequencies forward with constants that are not disclosed. The open problem is to complete this exact algebraic solution: derive and publish $Q_{4\pm}$ and $Q_{5\pm}$, obtain the explicit symmetry-adapted normal-mode eigenvectors, and carry out the force-constant inversion from the measured optically active spectrum together with a well-posedness analysis.

Background

Buckminsterfullerene, $C_{60}$, is a hollow carbon cage whose 60 atoms sit at the vertices of a truncated icosahedron - the soccer-ball pattern of 12 pentagons and 20 hexagons. Its equilibrium configuration is invariant under the full icosahedral group $I_h$ (the 60 proper rotations of the icosahedron combined with spatial inversion), the largest finite molecular point group. In the harmonic approximation the small oscillations of the atoms decompose into 174 normal modes (180 displacement degrees of freedom minus 6 rigid-body translations and rotations). In the classical force-constant ('ball-and-spring') treatment the potential energy is quadratic in internal coordinates - the changes of the 90 bond lengths and of the bond angles - with one spring constant per inequivalent interaction, and the squared frequencies $\lambda = \omega^2$ are the roots of a secular determinant (the characteristic determinant whose vanishing yields the allowed frequencies) of order $3N = 180$.

The cage has 90 bonds of two kinds - 60 edges shared by a pentagon and a hexagon ('single' bonds) and 30 edges shared by two hexagons ('double' bonds) - and two inequivalent classes of bond angles, with the third angular constant equal to one of them by symmetry. The minimal physically sensible model, used by Chou and Yang, therefore has four force constants: single-bond stretching $k_S$, double-bond stretching $k_D$, and two angular restoring constants $k_\alpha$ and $k_\beta$ (with $k_\gamma = k_\beta$). Direct expansion of the $180\times 180$ determinant for arbitrary constants is infeasible; before 1997 the problem was treated by mode counting and numerical diagonalization for chosen constants (Wu, Jelski and George 1987; Weeks and Harter 1989, who also obtained eigenfrequencies and normal-mode eigenvectors in special single-constant limits).

The symmetry route is to decompose the 180-dimensional displacement representation into irreducible representations (irreps) of $I_h = I \times Z_2$ - the inequivalent symmetry species, labeled by the five irreps of the rotation group $I$, of dimensions 1, 3, 3, 4 and 5, and by parity $g/u$ (even/odd under inversion). In the Chou-Yang notation these are $R_{1\pm}$, $R_{3\pm}$, $R_{3'\pm}$, $R_{4\pm}$ (four-dimensional) and $R_{5\pm}$ (five-dimensional). The optically active modes are the infrared-active $R_{3-}$ (four $T_{1u}$ fundamentals) and the Raman-active $R_{1+}$ and $R_{5+}$ (two $A_g$ and eight $H_g$ lines); the remaining modes are silent in first-order spectra.

Chou and Yang, extending the Abelian-subgroup reduction they had used with Deng for the Hückel (tight-binding pi-electron) Hamiltonian of $C_{60}$ (Deng and Yang 1992), organized the 60 atoms as ten images - five rotations about a chosen fivefold axis, each combined or not with inversion - of a basic set of six atoms. Displacements labeled by the fivefold-rotation eigenvalue $\eta^m$ ($\eta = e^{2\pi i/5}$, $m = -2,\ldots,2$) and inversion parity $p = \pm 1$ compress the full problem into six independent $18\times 18$ Hermitian blocks $M_{mp} = V_{mp} - \lambda K$ (sectors with $m$ and $-m$ are complex conjugates). Projection operators (matrices built from group characters that project onto a single irrep) are then applied inside this reduced space, and the secular determinant of each irreducible sector factorizes into a characteristic polynomial $Q_{\mu p}(\lambda)$ in $\lambda = \omega^2$, of degree equal to the multiplicity of that irrep in the 180-dimensional representation - at most 8. The coefficients are exact polynomial expressions in the four force constants and the geometric constant $r = \sqrt{5}$.

The 1997 paper presents $Q_{1\pm}$, $Q_{3\pm}$ and $Q_{3'\pm}$ explicitly, solves the one-dimensional sectors in closed form, verifies the zero-mode structure (translations in $R_{3-}$, rotations in $R_{3+}$) and the zero-frequency counting in special limits, and checks agreement with Wu-Jelski-George and Weeks-Harter in single-constant limits. What it leaves out, explicitly deferred to later work: the polynomials $Q_{4\pm}$ and $Q_{5\pm}$ of the four- and five-dimensional irreps (omitted on account of their length; the $R_{5+}$ sector carries eight of the ten Raman lines), the explicit normal-mode eigenvectors (the symmetry-adapted displacement pattern of each mode), and the numerical determination of the four force constants by fitting the computed roots to the measured infrared and Raman frequencies - including whether that inversion is unique and stable.

Problem Statement

Complete the exact algebraic solution of the vibration problem of the carbon-60 molecule in the Chou-Yang four-force-constant model, at the generality of the 1997 formulation. Concretely, completing this exact solution requires the characteristic polynomials $Q_{4\pm}(\lambda)$ and $Q_{5\pm}(\lambda)$, the symmetry-adapted normal-mode eigenvectors, and the force-constant inversion together with its well-posedness analysis. (1) Derive and present explicitly the characteristic polynomials $Q_{4\pm}(\lambda)$ and $Q_{5\pm}(\lambda)$ - the secular polynomials of the four-dimensional and five-dimensional irreducible representations $R_{4\pm}$ and $R_{5\pm}$ of the icosahedral group $I_h$ - as polynomials in $\lambda = \omega^2$ whose coefficients are exact algebraic expressions in the four force constants $k_S$ (single-bond stretching), $k_D$ (double-bond stretching), $k_\alpha$ and $k_\beta$ (angular restoring constants, $k_\gamma = k_\beta$), with the same conventions ($r = \sqrt{5}$, Abelian $(m,p)$ sector labeling, basic-set coordinates) as the published $Q_{1\pm}$, $Q_{3\pm}$, $Q_{3'\pm}$. (2) Obtain the corresponding explicit normal-mode eigenvectors: for each eigenvalue of each irreducible sector, the symmetry-adapted displacement pattern, expressed in the 18-component basic-set coordinates of the $(m,p)$ reduction and consistent with the phases $\eta^{m n}$ under the fivefold rotations and with inversion parity. (3) Carry out the inversion the paper defers: using the complete polynomial set, determine the four force constants from the experimentally measured infrared ($R_{3-}$) and Raman ($R_{1+}$, $R_{5+}$) frequencies of isolated $C_{60}$ reported in the experimental studies cited below (References 16 and 17), stating the specific measured lines and the measurement uncertainty adopted for each, and analyze the well-posedness of this inverse problem by returning an explicit quantitative answer to each of three questions - the uniqueness or multiplicity of the fitted constants, their stability under those measurement uncertainties, and the sensitivity of the computed spectrum to each constant.

References

1. T. T. Chou and Chen Ning Yang, 'Exact solution of the vibration problem for the carbon-60 molecule', Physics Letters A 235 (1997) 97-104, doi:10.1016/s0375-9601(97)00529-x, https://doi.org/10.1016/S0375-9601(97)00529-X
2. T. T. Chou and Chen Ning Yang, 'Vibration of the Carbon-60 Molecule', in Jingshin Theoretical Physics Symposium in Honor of Professor Ta-You Wu (eds. Jong-Ping Hsu and Leonardo Hsu, World Scientific, 1998) pp. 97-115, doi:10.1142/9789812816566_0009, https://doi.org/10.1142/9789812816566_0009
3. T. T. Chou, 'Symmetries and exact solution of the molecular vibration problem', Nuclear Physics A 684 (2001) 652-654, doi:10.1016/s0375-9474(01)00513-9, https://doi.org/10.1016/S0375-9474(01)00513-9
4. T. T. Chou and Chen Ning Yang, 'Exact solution of the vibration problem for the carbon-60 molecule', in Selected Papers of Chen Ning Yang II (World Scientific, 2013), doi:10.1142/9789814449021_0027, https://doi.org/10.1142/9789814449021_0027 (reprint of the 1997 paper)
5. Z. C. Wu, Daniel A. Jelski and Thomas F. George, 'Vibrational motions of buckminsterfullerene', Chemical Physics Letters 137 (1987) 291-294, doi:10.1016/0009-2614(87)80221-X, https://doi.org/10.1016/0009-2614(87)80221-X
6. David E. Weeks and William G. Harter, 'Rotation-vibration spectra of icosahedral molecules. II. Icosahedral symmetry, vibrational eigenfrequencies, and normal modes of buckminsterfullerene', The Journal of Chemical Physics 90 (1989) 4744-4771, doi:10.1063/1.456571, https://doi.org/10.1063/1.456571
7. Yuefan Deng and Chen Ning Yang, 'Eigenvalues and eigenfunctions of the Hückel Hamiltonian for carbon-60', Physics Letters A 170 (1992) 116-126, doi:10.1016/0375-9601(92)90784-J, https://doi.org/10.1016/0375-9601(92)90784-J
8. Jin-Quan Chen and Jia-Lun Ping, 'Algebraic Expressions of Irreducible Bases for the Molecule C60', Journal of Molecular Spectroscopy 186 (1997) 403-415, doi:10.1006/jmsp.1997.7448, https://doi.org/10.1006/jmsp.1997.7448
9. Jin-Quan Chen, Xiang-Jun Xin, Peng-Dong Fan and Jia-Lun Ping, 'Algebraic expressions of irreducible bases for molecules C20H20, C80, and C240', International Journal of Quantum Chemistry 73 (1999) 283-297, doi:10.1002/(SICI)1097-461X(1999)73:3<283::AID-QUA3>3.0.CO;2-5, https://doi.org/10.1002/(SICI)1097-461X(1999)73:3%3C283::AID-QUA3%3E3.0.CO;2-5
10. Shi-Hai Dong, Xi-Wen Hou and Zhong-Qi Ma, 'Irreducible bases in icosahedral group space', preprint arXiv:math-ph/9811025, https://arxiv.org/abs/math-ph/9811025
11. E. Yépez, A. Calles and J. J. Castro, 'A simple algorithm for the group theoretical classification of quantum states', Applied Mathematics and Computation 133 (2002) 119-130, doi:10.1016/S0096-3003(01)00226-0, https://doi.org/10.1016/S0096-3003(01)00226-0
12. E. Yépez, J. R. Soto, J. J. Castro and A. Calles, 'A simple algorithm for the group theoretical classification of quantum states II. The case of molecular electronic states', Applied Mathematics and Computation 164 (2005) 719-729, doi:10.1016/j.amc.2004.04.026, https://doi.org/10.1016/j.amc.2004.04.026
13. J. R. Soto, A. Calles and J. L. Morán-López, 'Normal modes of vibration of the C60 isotopomers 12C59-13C and 12C58-13C2', The Journal of Chemical Physics 113 (2000) 1055-1060, doi:10.1063/1.481885, https://doi.org/10.1063/1.481885
14. A. Zingoni, 'Group-theoretic exploitations of symmetry in computational solid and structural mechanics', International Journal for Numerical Methods in Engineering 79 (2009) 253-289, doi:10.1002/nme.2576, https://doi.org/10.1002/nme.2576
15. Janette L. Dunn, 'A Pictorial Visualization of Normal Mode Vibrations of the Fullerene (C60) Molecule in Terms of Vibrations of a Hollow Sphere', Journal of Chemical Education 87 (2010) 819-822, doi:10.1021/ed100263f, https://doi.org/10.1021/ed100263f
16. W. Krätschmer, K. Fostiropoulos and Donald R. Huffman, 'The infrared and ultraviolet absorption spectra of laboratory-produced carbon dust: evidence for the presence of the C60 molecule', Chemical Physics Letters 170 (1990) 167-170, doi:10.1016/0009-2614(90)87109-5, https://doi.org/10.1016/0009-2614(90)87109-5
17. Ping Zhou, Kai-An Wang, Ying Wang, P. C. Eklund, M. S. Dresselhaus, G. Dresselhaus and R. A. Jishi, 'Raman scattering in C60 and alkali-metal-saturated C60', Physical Review B 46 (1992) 2595-2605, doi:10.1103/PhysRevB.46.2595, https://doi.org/10.1103/PhysRevB.46.2595
18. A. G. Sukharev, 'Symmetry of Icosahedron and Optical Properties of Fullerene C60', Optics and Spectroscopy 129 (2021) 170-184, doi:10.1134/S0030400X21020168, https://doi.org/10.1134/S0030400X21020168