Nuclear Spin Statistical Weights Research Articles

Rotational spectroscopy of urea up to 500 GHz: The ground and eight excited vibrational states

Urea, ((NH2)2CO) is an important industrial and biological molecule. The new generation of the submillimeter telescopes, in particular ALMA, allows detection of extraterrestrial molecular spectra at unprecedented spatial and spectral resolutions. Some of the data presented here facilitated the recent definitive identification of urea in the interstellar medium (Belloche et al. Astron. & Astrophys. 628 (2019) A10:1-62). Its detection yields valuable insight into the mechanisms governing formation of complex prebiotic interstellar molecules. This paper reports spectroscopic analysis of the broadband rotational spectra of urea in the 210–500 GHz spectral range for the vibrational ground state and eight lowest excited vibrational states. The lowest excited vibrational state is only 61 cm−1 above the ground state and could be a useful temperature marker for astrophysical urea. The next two lowest pairs of excited states show significant interstate coupling, which has been successfully fitted with a Coriolis coupling model. Higher vibrational states also show evidence of significant coupling. Intensities of some rotational transitions of urea display nuclear spin statistical weights arising from the equivalence of the two ammonia groups, which enabled the assignment of symmetric and antisymmetric vibrational symmetry to the six lowest studied excited states. The vibrational energies of the new states have been determined from extensive relative intensity measurements and, together with vibrational symmetry, give insight into the highly anharmonic vibrational energy level structure.

Multipole-moment effects in ion-molecule reactions at low temperatures: part I - ion-dipole enhancement of the rate coefficients of the He+ + NH3 and He+ + ND3 reactions at collisional energies Ecoll/kB near 0 K.

The energy dependence of the rates of the reactions between He+ and ammonia (NY3, Y = {H,D}), forming NY2+, Y and He as well as NY+, Y2 and He, and the corresponding product branching ratios have been measured at low collision energies Ecoll between 0 and kB·40 K using a recently developed merged-beam technique [Allmendinger et al., ChemPhysChem, 2016, 17, 3596]. To avoid heating of the ions by stray electric fields, the reactions are observed within the large orbit of a highly excited Rydberg electron. A beam of He Rydberg atoms was merged with a supersonic beam of ammonia using a curved surface-electrode Rydberg–Stark deflector, which is also used for adjusting the final velocity of the He Rydberg atoms, and thus the collision energy. A collision-energy resolution of about 200 mK was reached at the lowest Ecoll values. The reaction rate coefficients exhibit a sharp increase at collision energies below ∼kB·5 K and pronounced deviations from Langevin-capture behaviour. The experimental results are interpreted in terms of an adiabatic capture model describing the rotational-state-dependent orientation of the ammonia molecules by the electric field of the He+ atom. The model faithfully describes the experimental observations and enables the identification of three classes of |JKMp〉 rotational states of the ammonia molecules showing different low-energy capture behaviour: (A) high-field-seeking states with |KM| ≥ 1 correlating to the lower component of the umbrella-motion tunnelling doublet at low fields. These states undergo a negative linear Stark shift, which leads to strongly enhanced rate coefficients; (B) high-field-seeking states subject to a quadratic Stark shift at low fields and which exhibit only weak rate enhancements; and (C) low-field-seeking states with |KM| ≥ 1. These states exhibit a positive Stark shift at low fields, which completely suppresses the reactions at low collision energies. Marked differences in the low-energy reactivity of NH3 and ND3—the rate enhancements in ND3 are more pronounced than in NH3—are quantitatively explained by the model. They result from the reduced magnitudes of the tunnelling splitting and rotational intervals in ND3 and the different occupations of the rotational levels in the supersonic beam caused by the different nuclear-spin statistical weights. Thermal capture rate constants are derived from the model for the temperature range between 0 and 10 K relevant for astrochemistry. Comparison of the calculated thermal capture rate coefficients with the absolute reaction rates measured above 27 K by Marquette et al. (Chem. Phys. Lett., 1985, 122, 431) suggests that only 40% of the close collisions are reactive.

Nonrigid water octamer: Computations with the 8-cube.

A 8-cube model of the fully nonrigid water octamer is considered within the 8-dimensional hyperoctahedral wreath product group with 10,321,920 operations and 185 irreducible representations by employing computational and mathematical techniques. For the two lowest-lying isomers of (H2 O)8 with D2d and S4 symmetries of a rigid (H2 O)8 , correlation tables and nuclear spin statistics are constructed for the tunneling splittings of the rotational levels are computed by a computational matrix polynomial generating function technique combined with Möbius inversion, and the relationship to the 8-cube multinomials are pointed out. Multinomial generating functions combined with the induced representation techniques are employed to compute and construct the nuclear spin species, nuclear spin statistical weights and tunneling splittings of rovibronic levels. We have also computed the spin statistical weights and tunneling splittings of the rotational levels for a semirigid water octamer within the wreath product Oh [S2 ] consisting of 12,288 operations.

LIMAO: Cross-platform software for simulating laser-induced alignment and orientation dynamics of linear-, symmetric- and asymmetric tops

A user-friendly and cross-platform software called Laser-Induced Molecular Alignment and Orientation simulator (LIMAO) has been developed. The program can be used to simulate within the rigid rotor approximation the rotational dynamics of gas phase molecules induced by linearly polarized intense laser fields at a given temperature. The software is implemented in the Java and Mathematica programming languages. The primary aim of LIMAO is to aid experimental scientists in predicting and analyzing experimental data representing laser-induced spatial alignment and orientation of molecules. Program summaryProgram title: LIMAOProgram Files doi: http://dx.doi.org/10.17632/x7h287gc9p.1Licensing provisions: GNU General Public License 3 (GPL)Programming languages: JAVA, MathematicaSupplementary material: LIMAO_SM.pdfNature of problem: LIMAO is a user-friendly software, which can be used to compute the temporal evolution of the laser-induced alignment and orientation [1-3] of gas phase molecules at a given temperature within the rigid rotor approximation. The user is asked to provide (i) molecular parameters such as rotational constant(s), a dipole moment, polarizabilities, and nuclear spin-statistical weights for the rotational levels, (ii) external field parameters such as time when the laser pulse becomes maximum, the pulse width, the intensity and carrier frequency of laser pulse(s), and (iii) temperature as input data. As an output, LIMAO generates the temporal evolution of the expectation values of physical quantities of molecules representing the spatial alignment and orientation.Solution method: LIMAO solves numerically the time-dependent Schrödinger equation of a rigid rotor molecule interacting with a time-dependent external field. The time-dependent rotational wave packet is expanded by the field-free stationary rotational eigenstates of the system. The solution is obtained by the following four steps: (i) the rotational eigenstates are derived and their symmetries within the corresponding molecular rotation group are assigned, (ii) the matrix representation of the time-dependent Hamiltonian, containing the interaction terms with an external light field as well as the matrix representations of physical quantities to be observed are constructed using the rotational eigenstates, (iii) the initial populations in the respective rotational eigenstates are determined by the Boltzmann distribution at a given temperature, and (iv) the wave packets are propagated in time and the expectation values of the relevant physical quantities are computed as a function of time.Additional comments including Restrictions and Unusual features: Simulations with the LIMAO software are limited to the rotational degrees of freedom. Light-induced vibrational excitations are neglected and light-induced electronic excitations are only accounted for in a perturbative manner through molecular polarizabilities. Therefore, the results obtained by LIMAO simulations should be treated with caution when dynamical processes such as light-induced deformation of geometrical structure, resonant vibrational excitation, and electronic excitation beyond the perturbative regime need to be considered.

Deuteration of ammonia in the starless core Ophiuchus/H-MM1

Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.

Nuclear spin species, spin statistical weights and rovibronic tunnelling splittings of non-rigid water hexamer

ABSTRACTPowerful multinomial generating functions together with the character table of the six-dimensional hyperoctahedral wreath product group S6[S2] of the non-rigid water hexamer are employed to obtain the nuclear spin species, nuclear spin multiplets and the total nuclear spin statistical weights of the rovibronic-tunnelling levels for both deuterated and regular forms. The methods are composed of algebraic generating functions involving 531,441 nuclear spin functions for (D2O)6 in a group of 46,080 operations and 65 irreducible representations of the S6[S2] group. Although the deuterated form of the non-rigid water hexamer possesses nuclear spin populations in 58 of the 65 possible symmetries of tunnelling levels, for the regular water hexamer only 23 of the 65 symmetries are populated with nuclear spin functions. The tunnelling splitting correlations of rovibronic levels of the water hexamer have been obtained from the Saykally's semi-rigid (G4) model to the fully non-rigid limit (G46080). The computed nuclear spin statistical weights for protonated form of (H2O)6 call for a reinterpretation of the previous assignment of the observed spectra. The results can also be applied to enumerate the hyperfine patterns using the nuclear spin multiplets and intensities in semi-rigid to fully non-rigid limits.

Reactive intermediates in (4)He nanodroplets: infrared laser Stark spectroscopy of dihydroxycarbene.

Singlet dihydroxycarbene (HOC̈OH) is produced via pyrolytic decomposition of oxalic acid, captured by helium nanodroplets, and probed with infrared laser Stark spectroscopy. Rovibrational bands in the OH stretch region are assigned to either trans,trans- or trans,cis-rotamers on the basis of symmetry type, nuclear spin statistical weights, and comparisons to electronic structure theory calculations. Stark spectroscopy provides the inertial components of the permanent electric dipole moments for these rotamers. The dipole components for trans, trans- and trans, cis-rotamers are (μa, μb) = (0.00, 0.68(6)) and (1.63(3), 1.50(5)), respectively. The infrared spectra lack evidence for the higher energy cis,cis-rotamer, which is consistent with a previously proposed pyrolytic decomposition mechanism of oxalic acid and computations of HOC̈OH torsional interconversion and tautomerization barriers.

The ethyl radical in superfluid helium nanodroplets: Rovibrational spectroscopy and ab initio computations

The ethyl radical has been isolated and spectroscopically characterized in (4)He nanodroplets. The band origins of the five CH stretch fundamentals are shifted by < 2 cm(-1) from those reported for the gas phase species [S. Davis, D. Uy, and D. J. Nesbitt, J. Chem. Phys. 112, 1823 (2000); T. Häber, A. C. Blair, D. J. Nesbitt, and M. D. Schuder, J. Chem. Phys. 124, 054316 (2006)]. The symmetric CH2 stretching band (v1) is rotationally resolved, revealing nuclear spin statistical weights predicted by G12 permutation-inversion group theory. A permanent electric dipole moment of 0.28 (2) D is obtained via the Stark spectrum of the v1 band. The four other CH stretch fundamental bands are significantly broadened in He droplets and lack rotational fine structure. This broadening is attributed to symmetry dependent vibration-to-vibration relaxation facilitated by the He droplet environment. In addition to the five fundamentals, three a1' overtone∕combination bands are observed, and each of these have resolved rotational substructure. These are assigned to the 2v12, v4 + v6, and 2v6 bands through comparisons to anharmonic frequency computations at the CCSD(T)∕cc-pVTZ level of theory.

ORTHO-TO-PARA ABUNDANCE RATIO (OPR) OF AMMONIA IN 15 COMETS: OPRs OF AMMONIA VERSUS14N/15N RATIOS IN CN

The ortho-to-para abundance ratio (OPR) of cometary molecules is considered to be one of the primordial characteristics of cometary ices. We present OPRs of ammonia (NH3) in 15 comets based on optical high-dispersion spectroscopic observations of NH2, which is a photodissociation product of ammonia in the gaseous coma. The observations were mainly carried out with the VLT/UVES. The OPR of ammonia is estimated from the OPR of NH2 based on the observations of the NH2 (0, 9, 0) vibronic band. The absorption lines by the telluric atmosphere are corrected and the cometary C2 emission lines blended with NH2 lines are removed in our analysis. The ammonia OPRs show a cluster between 1.1 and 1.2 (this corresponds to a nuclear spin temperature of ~30 K) for all comets in our sample except for 73P/Schwassmann-Wachmann 3 (73P/SW3). Comet 73P/SW3 (both B- and C-fragments) shows the OPR of ammonia consistent with nuclear spin statistical weight ratio (1.0) that indicates a high-temperature limit as nuclear spin temperature. We compared the ammonia OPRs with other properties (14N/15N ratios in CN, D/H ratios of water, and mixing ratios of volatiles). Comet 73P/SW3 is clearly different from the other comets in the plot of ammonia OPRs versus 14N/15N ratios in CN. The ammonia OPRs of 1.0 and lower 15N-fractionation of CN in comet 73P/SW3 imply that icy materials in this comet formed under warmer conditions than other comets. Comets may be classified into two groups in the plot of ammonia OPRs against 14N/15N ratios in CN.

Vibration–rotation-tunneling states of the benzene dimer: an ab initio study

An improved intermolecular potential surface for the benzene dimer is constructed from interaction energies computed by symmetry-adapted perturbation theory, SAPT(DFT), with the inclusion of third-order contributions. Twelve characteristic points on the surface have been investigated also using the coupled-cluster method with single, double, and perturbative triple excitations, CCSD(T), and triple-zeta quality basis sets with midbond functions. The SAPT and CCSD(T) results are in close agreement and provide the best representation of these points to date. The potential was used in calculations of vibration-rotation-tunneling (VRT) levels of the dimer by a method appropriate for large amplitude intermolecular motions and tunneling between multiple equivalent minima in the potential. The resulting VRT levels were analyzed with the use of the permutation-inversion full cluster tunneling (FCT) group G(576) and a chain of subgroups that starts from the molecular symmetry group C(s)(M) of the rigid dimer at its equilibrium C(s) geometry and leads to G(576) if all possible intermolecular tunneling mechanisms are feasible. Further information was extracted from the calculated wave functions. It was found, in agreement with the experimental data, that for all of the 54 G(576) symmetry species (with different nuclear spin statistical weights) the lower VRT states have a tilted T-shape (TT) structure; states with the parallel-displaced structure are higher in energy than the ground state of A symmetry by at least 30 cm(-1). The dissociation energy D(0) equals 870 cm(-1), while the depth D(e) of the TT minimum in the potential is 975 cm(-1). Hindered rotation of the cap in the TT structure and tilt tunneling lead to level splittings on the order of 1 cm(-1). Also intermolecular vibrations with excitation energies starting at a few cm(-1) were identified. A further small, but probably significant, level splitting was assigned to cap turnover, although in scans of the potential surface we could not find a plausible 'reaction path' for this process. Rotational constants were extracted from energy levels calculated for total angular momentum J = 0 and 1, and from expectation values of the inertia tensor. Although the end-over-end rotational constant B + C agrees well with the measured microwave spectra, there is disagreement with the measurements concerning the (a)symmetric rotor character of the benzene dimer. It is concluded from calculations for the 54 nuclear spin species that the microwave spectrum should show overlapping contributions from many different species. Another interesting conclusion regards the role of the quantum number K, for a prolate near-symmetric rotor the projection of the total angular momentum on the prolate axis. For the benzene dimer, K has a substantial effect on the energy levels associated with the intermolecular motions of the complex.

10.1007/s11449-008-3004-7

The rotational contours of the bands corresponding to electronic and electronic-vibrational transitions of the fluorescence excitation spectrum of jet-cooled carbazole complexes with one, two, and three water molecules have been studied. For the carbazole-(H2O)1 complex, two bands with a spectral shift of 0.57 cm−1 were recorded under exposure to radiation with a spectral width of 0.08 cm−1 at the frequency of the purely electronic transition and some other electronic-vibrational transitions. This is caused by the tunnel effect. The intensities of the shifted low-frequency bands are threefold weaker than those of high-frequency bands due to different values of nuclear spin-statistical weights. In the carbazole-(H2O)2 and carbazole-(H2O)3 complexes, water molecules are combined into a chain by the hydrogen bond, and the two ends of the chain are hydrogen-bonded to the carbazole molecule. The principal axes I A and I B of the moments of inertia in carbazole-(H2O)3 have different orientation compared to the other complexes considered, and this leads to an increase in the intensity of the Q-branch.

Rotational contours of electronic and electronic-vibrational bands of carbazole complexes with water cooled in a supersonic jet

The rotational contours of the bands corresponding to electronic and electronic-vibrational transitions of the fluorescence excitation spectrum of jet-cooled carbazole complexes with one, two, and three water molecules have been studied. For the carbazole-(H2O)1 complex, two bands with a spectral shift of 0.57 cm−1 were recorded under exposure to radiation with a spectral width of 0.08 cm−1 at the frequency of the purely electronic transition and some other electronic-vibrational transitions. This is caused by the tunnel effect. The intensities of the shifted low-frequency bands are threefold weaker than those of high-frequency bands due to different values of nuclear spin-statistical weights. In the carbazole-(H2O)2 and carbazole-(H2O)3 complexes, water molecules are combined into a chain by the hydrogen bond, and the two ends of the chain are hydrogen-bonded to the carbazole molecule. The principal axes I A and I B of the moments of inertia in carbazole-(H2O)3 have different orientation compared to the other complexes considered, and this leads to an increase in the intensity of the Q-branch.

High resolution infrared study of the 2ν9 and ν4 bands of 10BF2OH and 11BF2OH: evidence of large amplitude effects for the OH-torsion and OH-bending modes in the 92 and 41 excited states

High resolution (2–3 × 10−3 cm−1) Fourier transform infrared spectra of gas phase 10B and 11B enriched and natural samples of BF2OH (difluoroboric acid) were recorded at Wuppertal and Richland. Starting from the results of previous studies, it was possible to perform the first rovibrational analysis of the 2ν9 (first overtone of ν9, the OH torsion) and ν4 (BOH bending) bands located at about 1043.9 and 961.7 cm−1 and 1042.9 and 961.5 cm−1 for the 10BF2OH and 11BF2OH isotopic species, respectively. Numerous ‘classical’ perturbations were observed in the analysis of the 2ν9 and ν4 bands. The energy levels of the 92 bright state are indeed involved in a B-type Coriolis resonance with those of the 6191 dark state. The 41 levels are perturbed by a B-type Coriolis resonance and by an anharmonic resonance with the levels of the 7191 and the 6171 dark states, respectively. In addition large amplitude effects were observed for the 2ν9 and also, more surprisingly, the ν4 bands. This results in splittings of the energy levels of about 0.005 and 0.0035 cm−1 for the 92 and 41 states, respectively, which are easily observable in the P and R branches for both bands. The theoretical model used to reproduce the experimental levels accounts for the classical vibration–rotation resonances. Also the large amplitude torsional (or bending) effects are accounted for within the frame of the IAM (Internal Axis Method)-like approach. The Coriolis resonances between the two torsional (or bending) substates are taken into account by {Jx, Jz} non orthorhombic terms in the v-diagonal blocks. This means that the z-quantification axis deviates from the a inertial axis by an axis switching effect of ∼35° for the {92, 6191} system and of ∼16.6° for the {41, 7191, 6171} system of interacting vibrational states. The calculation of the relative line intensities for the 2ν9 and ν4 bands accounts for these axis switching effects as well as for the intensity alternation which is due to the nuclear spin statistical weights since the OH large amplitude torsion and/or bending motion results indeed in an exchange of the two fluorine nuclei.

Rovibronic symmetry and nuclear spin analysis of octanitrocubane

Octanitrocubane is shown to possess interesting nuclear spin statistics and rovibronic level patterns in its equilibrium chiral ( D 4) symmetry both due to 17O and 14N nuclear spin species. The populations and distributions of rovibronic levels will be influenced by the chiral symmetry and the nuclear spin species and spin statistical weights due to 17O, 14N and 13C species. We have computed the nuclear spin species and spin statistics and correlations of the rovibronic levels and spin species in different symmetries. Our predicted results could aid in the interpretation of 17O dynamic NMR, vibronic and rovibronic spectra of gas-phase octanitrocubane.

The millimeter-wave rotational spectrum of fluorobenzene

The room-temperature rotational spectrum of fluorobenzene was studied in the frequency region 167–318 GHz. Rotational transitions were assigned and measured in the ground vibrational state, and all six excited vibrational states with energies below 600 cm −1, i.e., v 11 = 1, v 11 = 2, v 18 b = 1, v 16 a = 1, v 16 b = 1, and v 6 a = 1. Accurate quartic-level spectroscopic constants were determined for all states, allowing spectral predictions well into the submillimeter region. The states v 18 b = 1 and v 16 a = 1 were found to be connected by a strong Coriolis interaction, which allowed precise determination of their energy separation, Δ E = 7.455088(3) cm −1. Unambiguous assignment of vibrational modes was made on the basis of the calculated inertial defect and nuclear spin statistical weights. Rotational constants for the 13C 4 isotopomer have also been redetermined and two new least-squares determinations of the geometry of fluorobenzene, r 0 and r m ( 1 L ) are reported.

Group theory, nuclear spin statistics and tunneling splittings of 1,3,5-triamino-2,4,6-trinitrobenzene

The symmetry group of the non-rigid 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) is constructed for the internal rotations of the nitro groups that exhibit low rotation barriers of 5 kcal/mole. The permutation group has 48 operations characterized by the wreath product S 3[S 2]. The nuclear spin statistical weights for TATB with 17O, 14N and 1H species are computed and the tunneling splitting patterns of the rotational and rovibronic levels of TATB are constructed for the first time. It is shown that the ground rovibronic state of TATB is split into A 1, A 3, T 1 and T 3 levels.

Nonrigid Group Theory, Tunneling Splittings, and Nuclear Spin Statistics of Water Pentamer: (H2O)5

The character table of the fully nonrigid water pentamer, (H2O)5, is derived for the first time. The group of all feasible permutations is the wreath product group S5[S2] and it consists of 3840 operations divided into 36 conjugacy classes and irreducible representations. We have shown that the full character table can be constructed using elegant matrix type generator algebra. The character table has been applied to the water pentamer by obtaining the nuclear spin statistical weights of the rovibronic levels and tunneling splittings of the fully nonrigid pentamer. We have also obtained the statistical weights and tunneling splittings of a semirigid deuterated pentamer that exhibits pseudorotation with an averaged C5h (G10) symmetry used in the assignment of vibration-rotation−tunneling spectra. It is also shown that the previously considered group G320 for water pentamer of feasible permutations is a subgroup of the full group and is the direct product of wreath product C5[S2] and the inversion group. The correlation tables have been constructed for the semirigid (G10) to nonrigid (G3840) groups for the rotational levels and tunneling levels. The nuclear spin statistical weights have also been derived for both the limits and through the use of subduced representations the corresponding information can be obtained for G(320) as well from G(3840).

Nuclear spin statistics of extended aromatic C 48N 12 azafullerene

We have presented the nuclear spin statistics of the novel extended aromatic C 48N 12 azafullerene. The nuclear spin multiplets and statistical weights of 14N spin-1 bosons are provided. In addition, we have also provided the 13C nuclear spin species and spin statistical weights of 13C 48N 12. The spin statistical weights and spin species show that the presence of 14N nuclei in the aromatic fullerene can provide unique experimental opportunity to investigate the nuclear spin species.

Computer-generated character tables and nuclear spin statistical weights: application to benzene dimer and methane dimer

A short GAP computer program is presented for computing nuclear spin statistical weights of rovibronic symmetry species in polyatomic molecules. As examples, benzene dimer and methane dimer are discussed.

Nuclear spin statistics of extended aromatic C48N12 azafullerene

We have presented the nuclear spin statistics of the novel extended aromatic C48N12 azafullerene. The nuclear spin multiplets and statistical weights of 14N spin-1 bosons are provided. In addition, we have also provided the 13C nuclear spin species and spin statistical weights of 13C48N12. The spin statistical weights and spin species show that the presence of 14N nuclei in the aromatic fullerene can provide unique experimental opportunity to investigate the nuclear spin species.