Abstract
The James Webb Space Telescope (JWST) opened a new era for the identification of molecular systems in the interstellar medium (ISM) by vibrational features. One group of molecules of increasing interest is cyano-derivatives of aromatic organic molecules, which have already been identified in the ISM on the basis of the analysis of rotational signatures, and so, are plausible candidates for the detection by the JWST. Benzonitrile considered in this work represents a suitable example for the validation of a computational strategy, which can be further applied for different, larger, and not-yet observed molecules. For this purpose, anharmonic simulations of infrared (IR) spectra have been compared with recent FTIR experimental studies. The anharmonic computations using the generalized second-order vibrational perturbation theory (GVPT2) in conjunction with a hybrid force field combining the harmonic part of revDSD-PBEP86-D3/jun-cc-pVTZ with anharmonic corrections from B3LYP-D3/SNSD show very good agreement with those in the experiment, with a mean error of for all fundamental transitions overall and only for the stretching fundamental at 4.49 . The inclusion of overtones up to three-quanta transitions also allowed the prediction of spectra in the near-infrared region, which shows distinct features due to overtones at the 2.26 and 1.52 . The remarkable accuracy of the GVPT2 results opens a pathway for the reliable prediction of spectra for a broader range of cyano-astroCOMs.
1 Introduction
Despite being initially considered too harsh and diluted for molecule creation and survival, it is now recognized that several chemical processes can happen in astronomic environments leading to more than 310 molecules already detected in the interstellar medium (ISM) or circumstellar shells (Endres et al., 2016). Several of these molecules have been detected in the last decade, in particular, a number of complex organic species (COMs), defined as astrophysically relevant organic molecules consisting of six or more atoms. Some COMs that contain C, H, N, O, and possibly also S atoms can be claimed as prebiotic (Fulvio et al., 2021). To date, more than 30 prebiotic molecules have been detected in Taurus molecular cloud 1 (TMC-1), a dust-enshrouded gaseous cloud located 400 light-years from the Sun in the Taurus constellation. Such progress and fast increase in new detections has become possible due to advances in instrumentation, in particular the sub-millimeter and radio domains, allowing the analysis of the lowest rotational lines where radiation can pass through dust-enshrouded clouds (Guélin and Cernicharo, 2022).
However, rotational transitions are not always suitable and accessible for a study, for instance, for exoplanet atmospheres or dense dark regions, for which vibrational spectroscopies are often the methods of choice. The increasing importance of vibrational transitions in astrochemical context is clearly represented by the James Webb Space Telescope (JWST), which, with the infrared (IR) (Frost et al., 2022) observations obtained by the mid-infrared (MIR) and near-IR (NIR) instruments integrated within the Integrated Science Instrument Module (ISIM, Near-Infrared Camera [NIRCam], Near-Infrared Spectrograph [NIRSpec], Mid-Infrared Instrument [MIRI], and Jet Propulsion Laboratory [JPL]), already provides superlative sensitivity, spectral resolution, and wavelength coverage compared to previous space telescopes, such as Herschel, that observe in the visible and ultraviolet spectra (see Figure 1). These emerging experimental data already provide new information about the major ices in molecular cloud cores just prior to their collapse to form protostars (McClure et al., 2023). Currently, the Mars 2020 Perseverance rover also searches for signs of organic matter, in the contexts of the emergence of life, as well as its consideration as a habitable planet, performing, among others, NIR studies using SuperCam (Perez et al., 2017; Eigenbrode et al., 2018; Lasne, 2021; Sharma et al., 2023; McIntosh et al., 2024). These missions yield a huge amount of new significant results, which need to be analyzed, also highlighting the urgent need for accurate reference data in the MIR-to-NIR (0.6–28 ) range (Öberg, 2016; Nazari et al., 2021; Zapata Trujillo et al., 2023; Fortenberry, 2024b).
FIGURE 1

Cyano-derivatives in the ISM and the year of detection by rotational (left panel) or vibrational (right panel) features. Left panel: astro- COMs and Herschel’s view of TMC-1 (Herschel Space Observatory, 2017) (credit: ESA/Herschel/NASA/JPL-Caltech TMC-1). Right panel: recently detected by the JWST and NGC 604 NIRCam Image (Webb Space Telescope, 2024) (credit: NASA/ESA/CSA/STScI).
Spectroscopic techniques are the key for the analysis of astronomical observations and the detection of molecules in the interstellar medium and other astrochemical environments, such as atmospheres or soils of exoplanets or planetary moons (Loru et al., 2022; Lemmens et al., 2022; Mackie et al., 2022; Peeters et al., 2021; K. Lemmens et al., 2023; Petrignani and Candian, 2022; Puzzarini et al., 2023). However, the investigation of the chemical composition of “astrochemical samples” is complicated due to the concomitant presence of many, possibly unknown species. The detection of molecules is based on the comparison of spectra from spectroscopic observations in space with reference experiments from the laboratory (Puzzarini et al., 2023). Unfortunately, the latter might be limited, incomplete, or difficult to obtain under appropriate conditions (Barone et al., 2015c; Fortenberry, 2024b). The optimal strategy is represented by the combination of experiments with the theoretical approaches (Barone et al., 2021; Puzzarini, 2022), supporting and/or complementing laboratory studies (Biczysko et al., 2018a; Barone and Puzzarini, 2023; Fortenberry, 2024a). An increasing role of computational spectroscopy is related to its increasing accuracy (Yang et al., 2021; Barone and Puzzarini, 2023), as well as the possibility to match different, often extreme environmental conditions possible in space (Barone et al., 2015c; Biczysko et al., 2018b; Zapata Trujillo et al., 2023; Fortenberry, 2024b).
Among all the molecules detected to date in the ISM, over 30% bear a nitrogen atom, which usually bonds to carbon, in a large fraction (over 80% of species) by the triple CN bond (Endres et al., 2016). This is in line with the increased dipole moment that facilitates detection by rotational features, as well as the abundance of the CN radical in the ISM (McGuire, 2022) and its extreme reactivity (McGuire et al., 2018; 2021) leading to the abundance of cyano-substituted derivatives (cyano-astroCOMs and CN-astroCOMs) in the ISM. Indeed, benzonitrile (B, cyano-benzene, H5CN), was the first six-membered aromatic compound detected in 2018 in the ISM (McGuire et al., 2018), toward the dark molecular cloud TMC-1, a well-studied region where most molecules were first observed. This detection was possible as the CN group attachment creates a permanent dipole moment, which is null in otherwise “silent” (poly-)aromatic hydrocarbons (PAHs). This breakthrough intensified the search for aromatic CN-astroCOMs, leading to the detection in the TMC-1 of 1- and 2-cyanonaphthalene (H7CN) isomers (McGuire et al., 2021) in 2021 and 2-cyanoindene (H7CN) (Sita et al., 2022) in 2022. Following the observation of the six-membered rings, 1- and 2-cyanocyclopentadiene (H5CN) (McCarthy et al., 2021; Lee et al., 2021), which contains a five-membered cycle, were also detected in the same source (see Figure 1). Further searches extending to other hetero-aromatic and bridged bicyclic compounds (Martin-Drumel et al., 2023) also aim at establishing the abundance ratio between CN-PAHs and their parent PAHs (Barone and Lazzari, 2023). To date, none of these molecules have been identified by their vibrational features, but the first cyano-compound, i.e., the O- anion, which is one of the most well-detected species in astrophysical ices, has been detected in the low-mass star-forming region ChameleonI by the JWST based on the peak centered at (4.62 ) (McClure et al., 2023) (see Figure 1).
In this work, we focus on the first observed cyano-astroCOM, i.e., benzonitrile (McGuire et al., 2018). The experimental microwave spectrum of benzonitrile was first studied simultaneously by Erlandsson (1954) and Lide (1954). Bak et al. conducted comprehensive analyses of benzonitrile utilizing infrared and centimeter-wave techniques, also including feasible mono-substituted isotopologs (Bak and Nielsen, 1960; Bak et al., 1962). This led to the first derivation of a substitution structure, , which was further refined by Casado et al. by incorporating multiple Q-branch transitions for the majority of mono-substituted isotopologs (Casado et al., 1971). Green and Harrison subsequently enhanced the analysis of the experimental IR spectrum, taking into account the lowest wavenumber modes (Green and Harrison, 1976). Recently, these initial spectroscopic investigations have been improved, also in view of the increasing accuracy and spectrum coverage requirements due to the astrochemical importance of B (Yamamoto et al., 2000; Kwon et al., 2003; Burova and Anashkin, 2007; Rajasekhar et al., 2022). The accurate equilibrium structure of benzonitrile has been determined by Rudolph et al. (2013) using two different, complementary techniques, namely, the theoretical estimate and semi-experimental . The most extended study is the 2022 synchrotron investigation of B in a very broad spectrum range, up to (11.1 eV), which also included the new detection of the gas-phase IR spectra at a resolution of (Rajasekhar et al., 2022). This work, along with the high-resolution far-infrared spectra collected in the range using synchrotron radiation at the SOLEIL facility (Zdanovskaia et al., 2022), represents the reference for the anharmonic computations of IR spectra.
This work focuses on the IR spectra in the MIR-NIR region, first comparing the simulated infrared spectra with the available experimental counterparts and then providing the prediction of not yet available spectral data, which can support either laboratory or astrochemical studies.
2 Computational details
In order to simulate the spectroscopic parameters of benzonitrile in its electronic ground state , the geometry optimization and harmonic and anharmonic vibrational computations are performed using GAUSSIAN 16 (Frisch et al., 2016). In geometry optimization, the tight convergence criteria (maximum forces and displacements smaller than Hartree/Bohr and Å, respectively), as required for the anharmonic computations, are used. The equilibrium structure, harmonic force constants, and first-order electric dipole moment derivatives have been computed using the double-hybrid density functional revDSD-PBEP86 (Santra et al., 2019), which has been recommended for spectroscopic studies of medium-sized biomolecules (Barone et al., 2020; Yang et al., 2021; Mehta et al., 2023). These computations have been performed in conjunction with the jun-cc-pVTZ (denoted hereafter as junTZ) basis set (Papajak et al., 2011), which provides the optimal accuracy/cost ratio, as recently discussed by Xu et al. (2024). Moreover, the B3LYP (Becke, 1993)/SNSD (Barone et al., 2014) level has been used in the anharmonic computations. For both density functional theory (DFT) functionals, the dispersion correction proposed by Grimme (2011) has been added using the D3 (Grimme et al., 2010) version with Becke–Johnson (BJ) damping (Grimme et al., 2011; Najibi and Goerigk, 2018). For brevity, hereafter, the revDSD-PBE86-D3/jun-cc-pVTZ and B3LYP-D3/SNSD levels will be denoted as revDSD and B3LYP.
Computations of the third- and fourth-order derivatives of the potential energy surface have been performed at the B3LYP level by the numerical differentiation (Barone, 2005; Bloino, 2015) of analytic second-order derivatives, while the cubic electric dipole moment surfaces have been obtained through numerical differentiations of the dipole moment derivatives. The revDSD equilibrium and harmonic computations have been combined with B3LYP anharmonic computations to create a hybrid model used in spectroscopic simulations. The consistency of these two sets of data has been checked automatically, as implemented in GAUSSIAN 16. In other words, the overlap between two sets of normal modes (two different levels of theory) is defined using the linear transformation, as proposed by Duschinsky (1937):
where and represent the two sets of mass-weighted normal coordinates. The Duschinsky matrix describes the projection of normal coordinate basis vectors on those of the other, allowing the automatic check of the normal mode consistency between the two levels of theory used to define the hybrid method. To ensure that the two sets of normal modes computed at different levels of theory are equivalent, a 90% cut-off for each coordinate was required.
This cost-effective (Xu et al., 2024) hybrid revDSD/B3LYP scheme has been further used to compute spectroscopic parameters using the second-order vibrational perturbation theory (VPT2) (Nielsen, 1951; Mills, 1972). The ground vibrational state rotational constants have been obtained from the revDSD equilibrium structure by adding vibrational corrections computed at the revDSD/B3LYP level, which also provided data allowing the determination of the quartic and sextic centrifugal-distortion constants (Puzzarini et al., 2010; Puzzarini, 2013). For the vibrational spectra, it is also necessary to account for the possible presence of anharmonic resonances (Amos et al., 1991; Martin et al., 1995; Barone, 2005; Vázquez and Stanton, 2007; Rosnik and Polik, 2014; Barone et al., 2014; Bloino et al., 2015; Bloino, 2015; Krasnoshchekov et al., 2015; Franke et al., 2021; Mendolicchio et al., 2021; Franke et al., 2021) by the generalized VPT2 (GVPT2) model (Bloino and Barone, 2012; Bloino et al., 2015), where nearly resonant contributions are removed from the perturbative treatment (leading to the deperturbed model, DVPT2) and treated in a second step variationally. Resonance definition and general recommendations on the applied computational procedures are described in detail in the tutorial review by Bloino et al. (2016). It should be noted that, although improved criteria to define automatic resonances have been proposed recently, they would have negligible impacts on the energies (Yang and Bloino, 2022). Overall, the GVPT2 scheme employed in this work has been successfully applied to medium-sized or larger biomolecules with up to 100 atoms (Fornaro et al., 2015; Fusè et al., 2019; Yatsyna et al., 2019; Green and Improta, 2020; Yang et al., 2021; Sheng et al., 2021), also in the astrochemical context (Zhao et al., 2021; McIntosh et al., 2024; Alberini et al., 2024), so it is a valuable tool to be employed also for significantly larger cyano-astroCOMs.
3 Results and discussion
3.1 Equilibrium structure and rotational parameters
Selected equilibrium structural parameters of benzonitrile calculated at the revDSD/junTZ level are shown in Figure 2, while Table 1 compares all equilibrium parameters with the semi-experimental equilibrium structure (Pulay et al., 1978) derived in the reference (Rudolph et al., 2013) by combining the experimental ground-state rotational constants for a set of isotopologs with rovibrational corrections derived from cubic force fields determined at the B3LYP level (Piccardo et al., 2015). Moreover, Table 1 reports the best estimated theoretical structure obtained using the composite scheme employing all-electron CCSD(T) and MP2 geometry optimizations, with basis sets up to the quintuple-zeta, reaching this way the complete basis set (CBS) limit, as well as computations at the CCSD(T)/ANO1 level (Rajasekhar et al., 2022). Structural parameters computed by all combinations of revDSD and B3LYP functionals with the junTZ and SNSD basis sets are provided in Supplementary Material, along with the Cartesian coordinates by revDSD/junTZ. All DFT structures agree very well with the reference, among which revDSD/junTZ shows the smallest mean absolute errors (MAEs) of approximately 0.0030 for bond lengths and 0.08° for the angles, which, in terms of the largest discrepancies, correspond to approximately 0.005 and 0.15°, respectively. The good quality structure with MAEs of 0.0045 and 0.13°, respectively, is also observed for B3LYP/SNSD, justifying its application in the hybrid scheme. Interestingly, the revDSD/junTZ structure is closer to the reference than the CCSD(T)/ANO1 structure, further proving the reliability of revDSD as a cost-effective computational model, allowing to derive accurate geometrical parameters (Ceselin et al., 2021; Barone and Lazzari, 2023).
FIGURE 2

Structural parameter of benzonitrile as computed at the rDSD/junTZ level, with bond lengths (black) in Å and angles (blue) in degrees.
TABLE 1
| a | b | CCSD(T)/ANO1c | rDSD/junTZ | |
|---|---|---|---|---|
| - | 1.3968 | 1.3962 | 1.4012 | 1.3993 |
| - | 1.3884 | 1.3882 | 1.3934 | 1.3905 |
| - | 1.3917 | 1.3917 | 1.3967 | 1.3941 |
| - | 1.4347 | 1.4359 | 1.4393 | 1.4360 |
| -H8 | 1.0780 | 1.0803 | 1.0823 | 1.0826 |
| -H7 | 1.0799 | 1.0803 | 1.0824 | 1.0830 |
| -H9 | 1.0800 | 1.0806 | 1.0828 | 1.0832 |
| -N13 | 1.1582 | 1.1583 | 1.1646 | 1.1629 |
| MAX | . | 0.0023 | 0.0064 | 0.0047 |
| MAEd | . | 0.0007 | 0.0043 | 0.0030 |
| C2C3 | 119.42 | 119.52 | 119.54 | 119.55 |
| C3C4 | 120.27 | 120.13 | 120.13 | 120.17 |
| C4C6 | 120.07 | 120.22 | 120.19 | 120.16 |
| C1C5 | 120.55 | 120.50 | 120.46 | 120.41 |
| C2H8 | 119.77 | 119.61 | 119.61 | 119.64 |
| C3H7 | 120.10 | 120.16 | 120.16 | 120.14 |
| MAX | . | 0.16 | 0.16 | 0.15 |
| MAEd | . | 0.09 | 0.09 | 0.08 |
Equilibrium parameters of benzonitrile (bond lengths in Å; angles in degrees).
aSemi-experimental from Rudolph et al. (2013), obtained by combining the experimental ground-state rotational constants for a set of isotopologs with rovibrational corrections derived from cubic force fields determined at the B3LYP level.
Theoretical best estimated from Rudolph et al. (2013), obtained using a composite ab initio approach based on CCSD(T) and MP2 all-electron optimizations with basis sets up to quintuple-zeta quality.
Theoretical structure obtained at the CCSD(T)/ANO1 level from Zdanovskaia et al. (2022).
Largest absolute (MAX) and mean absolute errors (MAEs) of the bond length and angles compared to the from Rudolph et al. (2013).
The final validation is provided by the direct comparison with the experiment, i.e., spectroscopic constants from Watson’s asymmetric rotor Hamiltonian (A-reduction, representation), which are given in Table 2. Interestingly, the vibrational ground-state rotational constants obtained based on the revDSD structures with anharmonic corrections computed at the revDSD/B3LYP level agree with experiment within 0.4%, which is even better than the CCSD(T)/ANO1 obtained by Zdanovskaia et al. (2022). Moreover, good agreement, again similar to the CCSD(T)/ANO1 results, is also obtained for the quartic and sextic centrifugal-distortion constants. Notably, in some cases, such as the sextic constant, very good agreement with the most recent global fit including previous and expanded mm-wave measurements (Zdanovskaia et al., 2022) is obtained, while older experimental data reported values smaller by 50%.
TABLE 2
| Experimenta | CCSD(T)/ANO1b | rDSD/junTZ//B3LYP/SNSD | |
|---|---|---|---|
| (MHz) | 5,655.265,428 | 5,616. | 5,638. |
| (MHz) | 1,546.8757715 | 1535. | 1,541. |
| (MHz) | 1,214.4040832 | 1,205. | 1,210. |
| (kHz) | 0.0452858 | 0.0437 | 0.0433 |
| (kHz) | 0.937983 | 0.923 | 0.922 |
| (kHz) | 0.24411 | 0.241 | 0.226 |
| (kHz) | 0.01101116 | 0.0106 | 0.0106 |
| (kHz) | 0.609187 | 0.593 | 0.592 |
| (Hz) | 0.000002486 | 0.00000230 | 0.00000226 |
| (Hz) | 0.0015586 | 0.00149 | 0.00150 |
| (Hz) | -0.007863 | -0.00761 | -0.00769 |
| (Hz) | [0.0066915] | 0.0066915 | 0.0067529 |
| (Hz) | 0.000001159 | 0.00000110 | 0.00000105 |
| (Hz) | 0.0007398 | 0.000755 | 0.000757 |
| (Hz) | 0.007480 | 0.00712 | 0.00714 |
Rotational spectroscopic constants for the ground vibrational state of benzonitrile.
Spectroscopic constants derived by Zdanovskaia et al. (2022) using the single-state approach based on the ground state
bSpectroscopic constants obtained by Zdanovskaia et al. (2022) using the VPT2 computations at the CCSD(T)/ANO1 level.
3.2 Vibrational properties and IR spectra
Table 3 compares harmonic vibrational wavenumbers and IR intensities with those computed at the CCSD(T)/ANO1 level, showing very good agreement with the average error of approximately 6 and largest differences of approximately , as well as qualitative agreement for IR intensities, with a MAE below 1 km/mol, and largest discrepancies of approximately 7 km/mol observed for the most intense bands and . Thus, the good accuracy of revDSD harmonic wavenumbers is also demonstrated for benzonitrile, in accordance with what has been observed based on the comparison with CCSD(T) results with the CBS extrapolation from MP2 computations (Pietropolli Charmet et al., 2022; Tasinato et al., 2022; Xu et al., 2024). In Table 3 and following, we have adapted mode numbering, mode description, and Wilson notation, as done by Rajasekhar et al. (2022).
TABLE 3
| Sym | CCSD(T)/ANO1a | rDSD/junTZ | ||
|---|---|---|---|---|
| IR int. | IR int. | |||
| A1 | ||||
| 3,219 | 3.53 | 3,221 | 3.85 | |
| 3,206 | 5.43 | 3,208 | 6.29 | |
| 3,189 | 0.01 | 3,189 | 0.00 | |
| 2,277 | 6.14 | 2,263 | 11.32 | |
| 1,645 | 0.12 | 1,653 | 0.12 | |
| 1,520 | 9.82 | 1,530 | 8.75 | |
| 1,213 | 0.16 | 1,221 | 0.30 | |
| 1,195 | 0.55 | 1,202 | 0.72 | |
| 1,043 | 2.67 | 1,051 | 3.34 | |
| 1,010 | 0.12 | 1,018 | 0.23 | |
| 763 | 1.38 | 769 | 1.49 | |
| 459 | 0.00 | 462 | 0.00 | |
| A2 | ||||
| 987 | 0.00 | 996 | 0.00 | |
| 860 | 0.00 | 866 | 0.00 | |
| 402 | 0.00 | 405 | 0.00 | |
| B1 | ||||
| 1,006 | 0.01 | 1,009 | 0.00 | |
| 939 | 2.81 | 945 | 2.83 | |
| 769 | 48.33 | 770 | 56.01 | |
| 696 | 26.77 | 679 | 19.74 | |
| 550 | 15.42 | 556 | 15.05 | |
| 378 | 0.69 | 382 | 0.64 | |
| 143 | 1.70 | 144 | 1.79 | |
| B2 | ||||
| 3,214 | 6.27 | 3,217 | 6.86 | |
| 3,198 | 1.94 | 3,200 | 2.52 | |
| 1,620 | 1.33 | 1,628 | 1.24 | |
| 1,471 | 6.71 | 1,481 | 6.18 | |
| 1,348 | 0.96 | 1,358 | 1.15 | |
| 1,314 | 1.97 | 1,329 | 1.46 | |
| 1,175 | 0.26 | 1,182 | 0.24 | |
| 1,096 | 2.96 | 1,104 | 3.60 | |
| 629 | 0.12 | 633 | 0.12 | |
| 546 | 0.14 | 552 | 0.24 | |
| 162 | 4.44 | 162 | 4.65 | |
| MAX | - | 15 | 7.7 | |
| MIN | - | -18 | -7.0 | |
| MAE | - | 6 | 0.8 | |
Harmonic wavenumbers and IR intensities compared with reference computed data.
Ref. (Zdanovskaia et al., 2022).
Largest positive (MAX), negative (MIN), and mean absolute errors (MAEs) of the harmonic wavenumbers compared with the CCSD(T)/ANO1 reference (Zdanovskaia et al., 2022).
Table 4 lists all fundamental anharmonic wavenumbers and IR intensities of benzonitrile computed at the rDSD/junTZ//B3PLYP/SNSD GVPT2 level, while selected overtones and combination bands are given in Table 5. The accuracy of the simulated IR spectra of benzonitrile in the 500 –4,000 range can be assessed by comparing with experimental results recorded in the gas phase (Kwon et al., 2003; Rajasekhar et al., 2022). Table 4 provides both spectra, with the latter showing a higher resolution of . This increased resolution allows us to identify and assign non-fundamental transitions, as shown in Figure 3. The GVPT2 computations show overall good agreement with the experiment with a MAE of approximately and the largest positive and negative errors of approximately . Moreover, the largest errors are all related to the C-H stretching vibrations, which contribute to the broad band with some additional side peaks, which has not been assigned. Considering that our simulation agrees within with the most intense peak assigned as , it could be expected that further analysis of experimental data, including non-fundamental transitions, would lead to some re-assignments. Excluding all C-H stretching vibrations from the statistics leads to the average errors of and maximum discrepancies within . The most important result is extremely good agreement, within , for the , the fingerprint vibration of benzonitrile, which is predicted at (4.49 ). This result can be compared with the very recent study where the B3LYP/N07D quadratic force field was combined with VPT2 computations with resonances included (see Esposito et al. (2024) for details) but performed with the SPECTRO code, yielding of . Indeed, GVPT2 B3LYP-D3/N07D computations in GAUSSIAN 16 lead to a similar result, with of . Such a huge discrepancy of approximately was not expected based on previous benchmark tests, highlighting the need for a dedicated validation. We hope that the proposed GVPT2 revDSD/B3LYP methodology will allow us to distinguish between the different cyano-astroCOMs observed in the ISM using the vibrations normally occurring in the broader region 2,200–2,400 (Császár and Fogarasi, 1989).
TABLE 4
| Sym | Mode descriptiona Wilson notation |
Ref.a | Ref.b | Current work | Assignments (PED)c | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Expt | Theoryd | Expt | Theorye | IR Int. scaled | IR int | |||||
| A1 | ||||||||||
| 20a | CH | 3,080 | 3,207 | 3,106 | 3,210 | 2.40 | 3,148 | 2.19 | CH | |
| 2 | CH | 3,071 | 3,196 | 3,066 | 3,208 | 6.54 | 3,078 | 3.90 | CH | |
| 7a | CH | 3,042 | 3,178 | 3,043 | 3,198 | 8.46 | 2,999 | 1.20 | CH | |
| CN | CN | 2,232 | 2,332 | 2,229 | 2,323 | 34.02 | 2,227 | 7.44 | CN, CC | |
| 8a | CC | 1,599 | 1,641 | 1,599 | 1,643 | 0.62 | 1,610 | 0.01 | CC, CCH | |
| 19a | CC | 1,492 | 1,582 | 1,491 | 1,509 | 8.39 | 1,496 | 6.98 | HCC, R sym, CC | |
| 13 | X-sens | 1,191 | 1,220 | 1,193 | 1,214 | 0.33 | 1,202 | 0.08 | CC, CCH, R sym | |
| 9a | Ring | 1,178 | 1,203 | 1,178 | 1,192 | 0.93 | 1,185 | 0.84 | CCH, CC | |
| 18a | CH | 1,027 | 1,050 | 1,027 | 1,045 | 3.79 | 1,039 | 1.90 | CC, R sym, CCH | |
| 12 | CH | 1,001 | 1,019 | 1,001 | 1,010 | 0.00 | 1,008 | 0.14 | CCCH, R asy | |
| 1 | X-sens | 769 | 774 | 758 | 769 | 2.05 | 764 | 1.11 | R tri, CC | |
| 6a | X-sens | 461 | 467 | 459 | 463 | 0.00 | 458 | 0.00 | R tri, CC | |
| A2 | ||||||||||
| 17a | CH | 978 | 1,002 | 975 | 984 | 0.00 | 968 | 0.00 | CCCH, R asy | |
| 10a | CH | 848 | 863 | 844 | 849 | 0.00 | 841 | 0.00 | CCCH | |
| 16a | CC | 401 | 410 | 398 | 408 | 0.00 | 395 | 0.00 | R asy, CCCH | |
| B1 | ||||||||||
| 5 | CH | 987 | 1,021 | 1,001 | 1,008 | 0.37 | 977 | 0.00 | R asy, CC | |
| 17b | CH | 925 | 954 | 926 | 935 | 2.50 | 917 | 2.86 | CCCH | |
| 11 | CH | 758 | 781 | 758 | 773 | 34.79 | 746 | 50.72 | CCCH, R tri, CCCC | |
| 4 | CC | 686 | 706 | 687 | 704 | 38.87 | 656 | 23.67 | R asy, CCCH | |
| 16b | X-sens | 548 | 573 | 547 | 566 | 17.79 | 531 | 12.09 | CCN, CCCH, R asy | |
| CN | X-sens | 381 | 392 | 372 | 388 | 0.61 | 371 | 0.72 | R asy, CCN | |
| 10b | CN | 172 | 147 | 141 | 147 | 1.93 | 138 | 1.90 | CCCC, CCN, R tri | |
| B2 | ||||||||||
| 20b | CH | 3,039 | 3,188 | 3,093 | 3,208 | 6.54 | 3,093 | 7.49 | CH | |
| 7b | CH | 3,027 | 3,204 | 3,027 | 3,190 | 3.95 | 3,068 | 0.86 | CH | |
| 8b | CC | 1,584 | 1,615 | 1,583 | 1,614 | 0.80 | 1,591 | 0.89 | CC | |
| 19b | CC | 1,448 | 1,481 | 1,448 | 1,462 | 6.69 | 1,454 | 5.73 | CCH, CC | |
| 14 | CC | 1,337 | 1,361 | 1,335 | 1,351 | 1.79 | 1,337 | 0.25 | CC, CCH | |
| 3 | CH | 1,298 | 1,319 | 1,288 | 1,322 | 0.46 | 1,298 | 1.58 | CCH, CC | |
| 9b | CH | 1,163 | 1,188 | 1,163 | 1,176 | 0.16 | 1,172 | 0.13 | CCH, CC | |
| 18b | CH | 1,071 | 1,105 | 1,071 | 1,099 | 4.35 | 1,093 | 2.52 | CCH, CC | |
| 6b | CCC | 629 | 641 | 633 | 0.12 | 628 | 0.11 | R sym | ||
| CN | CN | 551 | 570 | 547 | 559 | 0.27 | 547 | 0.09 | CCC, CCCN | |
| 15 | X-sens | 162 | 169 | 167 | 4.65 | 158 | 4.69 | R asy, CCCN | ||
| All | MAXf | 42 | ||||||||
| MINf | -44 | |||||||||
| MAEf | 11 | |||||||||
| Exclude all CH | MAXg | 22 | ||||||||
| MINg | -31 | |||||||||
| MAEg | 8 | |||||||||
Fundamental wavenumbers and IR intensities computed at the GVPT2//revDSD/junTZ//B3LYP/SNSD level for benzonitrile compared with reference experimental and computed data.
Ref (Kwon et al., 2003).
Ref (Rajasekhar et al., 2022).
Normal mode assignments, , , , , and tri denote the stretching, in-plane bending, out-of-plane bending, torsion, and trigonal deformation, respectively. “sym” and “asy” stand for symmetrical and asymmetric deformation, respectively.
Scaled harmonic computations at the B3LYP/6-311++G(2df,2pd) level of theory.
Scaled harmonic computations at the B3LYP/aug-cc-pVDZ level of theory.
Largest positive (MAX), negative (MIN), and mean absolute errors (MAEs) of the benzonitrile fundamental wavenumbers compared with the experiment by Rajasekhar et al. (2022).
Largest positive (MAX), negative (MIN), and mean absolute errors (MAEs) in the benzonitrile fundamental wavenumbers, with all C-H stretching excluded, compared with the experiment by Rajasekhar et al. (2022).
TABLE 5
| Experiment | Current work | |||
|---|---|---|---|---|
| Assign. | Assign. | IR int. | ||
| 2 | 282a | 2 | 275 | 0.007 |
| + | 303a | + | 296 | 0 |
| 2 | 323a | 2 | 315 | 0.01 |
| 2+2 | 1,393b | + | 1,397 | 0.41 |
| 2+ | 1,688b | + | 1,658 | 1.37 |
| + | 1,769b | + | 1,754 | 1.30 |
| +3 | 1,816b | + | 1,803 | 1.69 |
| + | 1,900b | +, + | 1,883 | 2.83 |
| + | 1,970b | + | 1,942 | 3.23 |
| + | 2,178b | 2 | 2,178 | 0.02 |
Non-fundamental band wavenumbers and IR intensities computed at the GVPT2//revDSD/junTZ//B3LYP/SNSD level for benzonitrile compared with reference experimental data.
High-resolution IR spectra in the gas phase obtained by Zdanovskaia et al. (2022).
Gas-phase IR spectra obtained by Rajasekhar et al. (2022).
FIGURE 3

IR spectrum of benzonitrile in the range 500–4,000 . Computed stick spectra were broadened by Lorentzian functions with half-width at half-maximum (HWHM) of 2 . The experimental IR gas-phase spectrum of benzonitrile (Rajasekhar et al., 2022) is shown for comparison.
A direct comparison between the spectra given in Figure 3 highlights that the GVPT2 computation not only correctly predicts fundamental bands but also a pattern of five distinct non-fundamental bands in the range, allowing to correct their assignment with respect to the tentative one reported by Rajasekhar et al. (2022). These non-fundamental bands are reported in Table 5, together with those observed in the far-infrared spectra in a high-resolution FTIR experiment (Zdanovskaia et al., 2022). The overestimated intensity of the band at the should be noted (computed as ), while the experimental spectra show two similar intensity peaks in this range, the second one being at . However, this discrepancy needs to be linked to the harmonic values, which already predict the intensity of as twice that of . Notably, CCSD(T)/ANO-1 harmonic IR intensities yield the same pattern of these two bands as revDSD/junT. In order to provide more information about this discrepancy, a dedicated benchmark analysis, which would require appropriate numerical data on integrated intensities Charmet et al. (2013), not available at present, would be required.
Overall, the good accuracy of our simulations, for both fundamental and non-fundamental transitions, allows us to predict the spectra in the NIR region, which is shown in Figures 4A and B for and , respectively (the whole spectrum is also reported in Supplementary Material). The most pronounced bands in and are related to the combinations of with in-plane ring deformations and overtones, respectively. Similarly, at the higher energies, there are combinations of two quanta transitions with in-plane ring deformations () and 3 second overtones (). Although it is expected that the accuracy of GVPT2 results decreases for the higher-quanta transitions, the error bars for the first and second overtones can be estimated based on the fundamental bands (Barone et al., 2015a). For instance, in the case of formaldehyde (Biczysko et al., 2018a), good accuracy within has been obtained for fundamentals and first overtones of and , while a lower accuracy of of approximately transfers to errors of and for 2 and 3, respectively. This allows us to provide a reliable prediction of overtones, with distinct 2 transition at 4,426 (2.26 ) and a significantly weaker second overtone 3 at 6,598 (1.52 ).
FIGURE 4

NIR spectra simulated at the GVPT2//rDSD/junTZ//B3LYP/SNSD level. The spectra were broadened by Lorentzian functions with HWHM of 2 .
4 Conclusion and astrochemical implications
Good accuracy of simulated MIR spectra, confirmed by a comparison with available experimental results, allowed us to provide predictions regarding the “missing” data on relevance for the astrochemical observations, in particular, concerning the NIR region.
The availability of NIR references is important for the interpretation of data collected during the Mars 2020 (Williford et al., 2018) space mission by instruments such as SuperCam, for incoming ExoMars 2022 (ESZ-Roscosmos) (Vago et al., 2017), as well as for the JWST observations by the NIRSpec (McClure et al., 2023). The advantage of NIR is its lower spectrum congestion than MIR, which is also clearly visible in Figures 3, 4. However, these reference NIR data are often very scarce and were not available even for benzonitrile prior to this work. It has been already highlighted that anharmonic computations provide significant support in the analysis of experimental results and identification of plausible molecules by NIR features (Fornaro et al., 2020; Alberini et al., 2024).
Computational spectroscopy can also support the identification of B in other spectral regions, also considering effects due to the interactions with cosmic rays (Öberg, 2016; Arumainayagam et al., 2021), relevant for different astrochemical environments, from the ISM to planetary atmospheres or soil. These interactions can initiate different processes within molecules, depending on the photon energies, and can lead to electronic excitation within neutral molecules or create ions by ejecting off the valence or even inner-layer electrons. Interactions with photons can lead to the creation of new molecules, or their damage, but are also relevant for extending observable spectral ranges (Öberg, 2016).
Extensive laboratory experimental investigation of the photoabsorption spectra of benzonitrile recorded using synchrotron radiation in (4.3–11.1 eV, 0.111–0.286 ), which encompasses several neutral and ionic excited states, as recently reported by Rajasekhar et al. (2022). From a computational perspective, these processes can be simulated by means of vibronic computations (Bloino et al., 2016; Barone et al., 2021), which have been demonstrated to allow us to decipher a broad range of energies by the combination of two-state electronic transitions for a series of halogenated benzene (Palmer et al., 2015a; b). Moreover, first-principle spectral simulations also allow us to obtain reference data for unstable species difficult to study in the laboratories and to improve the resolution and predict spectra at a broad range of temperatures (Zhao et al., 2021). Computational spectroscopy studies combining anharmonic vibrational and vibronic simulations covering the broad range from MIR at approximately 20 up to even 20 nm in a high-energy photoelectron range can be extended toward other cyano-astroCOMs, supporting their possible detection.
The most relevant are those based on aromatic systems, such as 1- and 2-cyanocyclopentadiene (McCarthy et al., 2021; Lee et al., 2021), and 1- and 2-cyanonaphthalene (McGuire et al., 2021), which have been already discovered in TMC-1. However, it can be expected that similar accuracy can also be obtained for aliphatic systems, for which several conformers can be present (Barone et al., 2013; 2015b). This situation was highlighted by the recent discovery of five cyano-derivatives of propene (CHCH3), based on the QUIJOTE line survey of TMC-1 (Cernicharo et al., 2022). Such computations, including electronic spectra for benzonitrile and MIR to PES spectra for other -astroCOM species, are deferred to subsequent works, within the framework of development of the COSY-ASTRO dataset (COSY COST Action CA21101, 2024).
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
YX: formal analysis, investigation, validation, writing–original draft, writing–review and editing, data curation, and visualization. MB: formal analysis, investigation, validation, writing–original draft, writing–review and editing, conceptualization, funding acquisition, methodology, project administration, resources, and supervision.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was performed within COST Action CA21101 “Confined molecular systems: from a new generation of materials to the stars” (COSY) supported by COST (European Cooperation in Science and Technology). This work was supported by the National Natural Science Foundation of China (Grant No. 31870738).
Acknowledgments
MB acknowledges support from the COST Action “COSY,” in particular, discussions with Action Chair Prof. María Pilar de Lara-Castells, Action Vice-Chair Prof. Cristina Puzzarini, and WG5 “Confined Systems in Astrochemistry” Leader Prof. Lauri Halonen. MB and XY thank Prof. Julien Bloino for fruitful discussions on anharmonic computations.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2024.1439194/full#supplementary-material
References
1
Alberini A. Fornaro T. García-Florentino C. Biczysko M. Poblacion I. Aramendia J. et al (2024). Investigating the stability of aromatic carboxylic acids in hydrated magnesium sulfate under uv irradiation to assist detection of organics on mars. Sci. Rep.14, 15945. 10.1038/s41598-024-66669-8
2
Amos R. D. Handy N. C. Green W. H. Jayatilaka D. Willetts A. Palmieri P. (1991). Anharmonic vibrational properties of ch2f2: a comparison of theory and experiment. J. Chem. Phys.95, 8323–8336. 10.1063/1.461259
3
Arumainayagam C. R. Herbst E. Heays A. N. Mullikin E. Farrah M. Mavros M. G. (2021). “Extraterrestrial photochemistry: principles and applications,” in Prebiotic photochemistry: from urey–miller-like experiments to recent findings. Editors SaijaF.CassoneG. (Cambridge, UK: The Royal Society of Chemistry), 9–36. 10.1039/9781839164354-00009
4
Bak B. Christensen D. Dixon W. B. Hansen-Nygaard L. Rastrup-Andersen J. (1962). Benzene ring distortion by one substituent. Microwave determination of the complete structure of benzonitrile. J. Chem. Phys.37, 2027–2031. 10.1063/1.1733422
5
Bak Bø. Nielsen J. T. (1960). Infrared absorption spectra of benzonitrile and its monodeuterated derivatives: tentative assignment of fundamentals. Z. Elektrochem. Ber. Bunsenges. Phys. Chem.64, 560–562. 10.1002/bbpc.196000003
6
Barone V. (2005). Anharmonic vibrational properties by a fully automated second-order perturbative approach. J. Chem. Phys.122, 014108. 10.1063/1.1824881
7
Barone V. Alessandrini S. Biczysko M. Cheeseman J. R. Clary D. C. McCoy A. B. et al (2021). Computational molecular spectroscopy. Nat. Rev. Methods Prim.1, 38–27. 10.1038/s43586-021-00034-1
8
Barone V. Biczysko M. Bloino J. (2014). Fully anharmonic ir and Raman spectra of medium-size molecular systems: accuracy and interpretation. Phys. Chem. Chem. Phys.16, 1759–1787. 10.1039/C3CP53413H
9
Barone V. Biczysko M. Bloino J. Cimino P. Penocchio E. Puzzarini C. (2015a). CC/DFT route toward accurate structures and spectroscopic features for observed and elusive conformers of flexible molecules: pyruvic acid as a case study. J. Chem. Theory Comput.11, 4342–4363. 10.1021/acs.jctc.5b00580
10
Barone V. Biczysko M. Bloino J. Cimino P. Penocchio E. Puzzarini C. (2015b). Cc/dft route toward accurate structures and spectroscopic features for observed and elusive conformers of flexible molecules: pyruvic acid as a case study. J. Chem. Theory Comput.11, 4342–4363. 10.1021/acs.jctc.5b00580
11
Barone V. Biczysko M. Bloino J. Puzzarini C. (2013). Glycine conformers: a never-ending story?Phys. Chem. Chem. Phys.15, 1358–1363. 10.1039/C2CP43884D
12
Barone V. Biczysko M. Puzzarini C. (2015c). Quantum chemistry meets spectroscopy for astrochemistry: increasing complexity toward prebiotic molecules. Acc. Chem. Res.48, 1413–1422. 10.1021/ar5003285
13
Barone V. Ceselin G. Fusè M. Tasinato N. (2020). Accuracy meets interpretability for computational spectroscopy by means of hybrid and double-hybrid functionals. Front. Chem.8, 584203. 10.3389/fchem.2020.584203
14
Barone V. Lazzari F. (2023). Hunting for complex organic molecules in the interstellar medium: the role of accurate low-cost theoretical geometries and rotational constants. J. Phys. Chem. A127, 10517–10527. 10.1021/acs.jpca.3c06649
15
Barone V. Puzzarini C. (2023). Gas-phase computational spectroscopy: the challenge of the molecular bricks of life. Annu. Rev. Phys. Chem.74, 29–52. 10.1146/annurev-physchem-082720-103845
16
Becke A. D. (1993). Density-functional thermochemistry. iii. the role of exact exchange. J. Chem. Phys.98, 5648–5652. 10.1063/1.464913
17
Biczysko M. Bloino J. Puzzarini C. (2018a). Computational challenges in astrochemistry. WIREs Comput. Mol. Sci.8, e1349. 10.1002/wcms.1349
18
Biczysko M. Krupa J. Wierzejewska M. (2018b). Theoretical studies of atmospheric molecular complexes interacting with NIR to UV light. Faraday Discuss.212, 421–441. 10.1039/C8FD00094H
19
Bloino J. (2015). A vpt2 route to near-infrared spectroscopy: the role of mechanical and electrical anharmonicity. J. Phys. Chem. A119, 5269–5287. 10.1021/jp509985u
20
Bloino J. Baiardi A. Biczysko M. (2016). Aiming at an accurate prediction of vibrational and electronic spectra for medium-to-large molecules: an overview. Int. J. Quantum Chem.116, 1543–1574. 10.1002/qua.25188
21
Bloino J. Barone V. (2012). A second-order perturbation theory route to vibrational averages and transition properties of molecules: general formulation and application to infrared and vibrational circular dichroism spectroscopies. J. Chem. Phys.136, 124108. 10.1063/1.3695210
22
Bloino J. Biczysko M. Barone V. (2015). Anharmonic effects on vibrational spectra intensities: infrared, Raman, vibrational circular dichroism and Raman optical activity. J. Phys. Chem. A119, 11862–11874. 10.1021/acs.jpca.5b10067
23
Burova T. G. Anashkin A. A. (2007). Quantum-mechanical calculation of the intensity distribution in the resonance Raman spectrum of benzonitrile. Opt. Spectrosc.102, 825–828. 10.1134/S0030400X07060045
24
Casado J. Nygaard L. Sørensen G. O. (1971). Microwave spectra of isotopic benzonitriles. Refined molecular structure of benzonitrile. J. Molec. Struct.8, 211–224. 10.1016/0022-2860(71)80056-X
25
Cernicharo J. Fuentetaja R. Cabezas C. Agúndez M. Marcelino N. Tercero B. et al (2022). Discovery of five cyano derivatives of propene with the QUIJOTE line survey. Astron. Astrophys.663, L5. 10.1051/0004-6361/202244255
26
Ceselin G. Barone V. Tasinato N. (2021). Accurate biomolecular structures by the nano-LEGO approach: pick the bricks and build your geometry. J. Chem. Theory Comput.17, 7290–7311. 10.1021/acs.jctc.1c00788
27
Charmet A. P. Stoppa P. Tasinato N. Giorgianni S. Barone V. Biczysko M. et al (2013). An integrated experimental and quantum-chemical investigation on the vibrational spectra of chlorofluoromethane. J. Chem. Phys.139, 164302. 10.1063/1.4825380
28
COSY COST Action CA21101 (2024). COSY Confined molecular systems: from a new generation of materials to the stars CA21101 COST Action.
29
Császár A. G. Fogarasi G. (1989). Scaled quantum mechanical (SQM) force field and theoretical vibrational spectrum for benzonitrile. Spectrochim. Acta A45, 845–854. 10.1016/0584-8539(89)80222-3
30
[Dataset] Herschel Space Observatory (2017). Herschel’s view of the Taurus molecular cloud annotated. Available at: https://www.herschel.caltech.edu/image/nhsc2017-003a.
31
Duschinsky F. (1937). The importance of the electron spectrum in multi atomic molecules. concerning the franck-condon principle. Acta Physicochim. URSS7, 551.
32
Eigenbrode J. L. Summons R. E. Steele A. Freissinet C. Millan M. Navarro-González R. et al (2018). Organic matter preserved in 3-billion-year-old mudstones at Gale crater, Mars. Science360, 1096–1101. 10.1126/science.aas9185
33
Endres C. P. Schlemmer S. Schilke P. Stutzki J. Müller H. S. P. (2016). The cologne database for molecular spectroscopy, CDMS, in the virtual atomic and molecular data centre, VAMDC. J. Mol. Spectrosc.327, 95–104. 10.1016/j.jms.2016.03.005
34
Erlandsson G. (1954). Microwave spectrum of benzonitrile. J. Chem. Phys.22, 1152. 10.1063/1.1740313
35
Esposito V. J. Fortenberry R. C. Boersma C. Maragkoudakis A. Allamandola L. J. (2024). CN stretches around 4.4 microns dominate the IR absorption spectra of cyano-polycyclic aromatic hydrocarbons. Mon. Not. R. Astron. Soc. Lett.531, L87–L93. 10.1093/mnrasl/slae037
36
Fornaro T. Brucato J. R. Poggiali G. Corazzi M. A. Biczysko M. Jaber M. et al (2020). UV irradiation and near infrared characterization of laboratory mars soil analog samples. Front. Astron. Space Sci.7. 10.3389/fspas.2020.539289
37
Fornaro T. Carnimeo I. Biczysko M. (2015). Toward feasible and comprehensive computational protocol for simulation of the spectroscopic properties of large molecular systems: the anharmonic infrared spectrum of uracil in the solid state by the reduced dimensionality/hybrid vpt2 approach. J. Phys. Chem. A119, 5313–5326. 10.1021/jp510101y
38
Fortenberry R. C. (2024a). Quantum chemistry and astrochemistry: a match made in the heavens. J. Phys. Chem. A128, 1555–1565. 10.1021/acs.jpca.3c07601
39
Fortenberry R. C. (2024b). A vision for the future of astrochemistry in the interstellar medium by 2050. ACS Phys. Chem. Au4, 31–39. 10.1021/acsphyschemau.3c00043
40
Franke P. R. Stanton J. F. Douberly G. E. (2021). How to VPT2: accurate and intuitive simulations of CH stretching infrared spectra using VPT2+K with large effective Hamiltonian resonance treatments. J. Phys. Chem. A125, 1301–1324. 10.1021/acs.jpca.0c09526
41
Frisch M. J. Trucks G. W. Schlegel H. B. Scuseria G. E. Robb M. A. Cheeseman J. R. et al (2016). Gaussian 16 rev. C.01. Wallingford CT: Gaussian, Inc.
42
Frost A. J. Lau R. M. Burtscher L. Packham C. Tasker E. J. Rees G. A. et al (2022). IR 2022: an infrared-bright future for ground-based IR observatories in the era of JWST. Nat. Astron.6, 772–773. 10.1038/s41550-022-01733-9
43
Fulvio D. Potapov A. He J. Henning T. (2021). Astrochemical pathways to complex organic and prebiotic molecules: experimental perspectives for in situ solid-state studies. Life11, 568. 10.3390/life11060568
44
Fusè M. Mazzeo G. Longhi G. Abbate S. Masi M. Evidente A. et al (2019). Unbiased determination of absolute configurations by vis-à-vis comparison of experimental and simulated spectra: the challenging case of diplopyrone. J. Phys. Chem. B123, 9230–9237. PMID: 31580674. 10.1021/acs.jpcb.9b08375
45
Green J. A. Improta R. (2020). Vibrations of the guanine–cytosine pair in chloroform: an anharmonic computational study. Phys. Chem. Chem. Phys.22, 5509–5522. 10.1039/C9CP06373K
46
Green J. H. S. Harrison D. J. (1976). Vibrational spectra of benzene derivatives—XVII. Benzonitrile and substituted benzonitriles. Spectrochim. Acta A32, 1279–1286. 10.1016/0584-8539(76)80166-3
47
Grimme S. (2011). Density functional theory with london dispersion corrections. WIREs Comput. Mol. Sci.1, 211–228. 10.1002/wcms.30
48
Grimme S. Antony J. Ehrlich S. Krieg H. (2010). A consistent and accurate ab initio parametrization of density functional dispersion correction (dft-d) for the 94 elements h-pu. J. Chem. Phys.132, 154104. 10.1063/1.3382344
49
Grimme S. Ehrlich S. Goerigk L. (2011). Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem.32, 1456–1465. 10.1002/jcc.21759
50
Guélin M. Cernicharo J. (2022). Organic molecules in interstellar space: latest advances. Front. Astron. Space Sci.9. 10.3389/fspas.2022.787567
51
Krasnoshchekov S. V. Vogt N. Stepanov N. F. (2015). Ab initio anharmonic analysis of vibrational spectra of uracil using the numerical-analytic implementation of operator Van Vleck perturbation theory. J. Phys. Chem. A119, 6723–6737. 10.1021/acs.jpca.5b03241
52
Kwon C. H. Kim H. L. Kim M. S. (2003). Vibrational analysis of vacuum ultraviolet mass-analyzed threshold ionization spectra of phenylacetylene and benzonitrile. J. Phys. Chem. A107, 10969–10975. 10.1021/jp036391w
53
Lasne J. (2021). Heterogeneous physical chemistry in the atmospheres of earth, mars, and venus: perspectives for rocky exoplanets. ACS Earth Space Chem.5, 149–162. 10.1021/acsearthspacechem.0c00126
54
Lee K. L. K. Changala P. B. Loomis R. A. Burkhardt A. M. Xue C. Cordiner M. A. et al (2021). Interstellar detection of 2-cyanocyclopentadiene, C5H5CN, a second five-membered ring toward TMC-1. Astrophys. J. Lett.910, L2. 10.3847/2041-8213/abe764
55
Lemmens A. K. Rap D. B. Brünken S. Buma W. J. Rijs A. M. (2022). Polycyclic aromatic hydrocarbon growth in a benzene discharge explored by IR-UV action spectroscopy. Phys. Chem. Chem. Phys.24, 14816–14824. 10.1039/D2CP01631A
56
Lemmens K. Mackie C. Candian A. Lee J. Tielens M. Anouk M. R. et al (2023). Size distribution of polycyclic aromatic hydrocarbons in space: an old new light on the 11.2/3.3 m intensity ratio. Faraday Discuss.245, 380–390. 10.1039/D2FD00180B
57
Lide D. R. (1954). Microwave spectrum and structure of benzonitrile. J. Chem. Phys.22, 1577–1578. 10.1063/1.1740461
58
Loru D. Steber A. L. Thunnissen J. M. M. Rap D. B. Lemmens A. K. Rijs A. M. et al (2022). New potential candidates for astronomical searches discovered in the electrical discharge of the PAH naphthalene and acetonitrile. J. Mol. Spectrosc.386, 111629. 10.1016/j.jms.2022.111629
59
Mackie C. J. Candian A. Lee T. J. Tielens A. G. G. M. (2022). Anharmonicity and the IR emission spectrum of neutral interstellar PAH molecules. J. Phys. Chem. A126, 3198–3209. 10.1021/acs.jpca.2c01849
60
Martin J. M. L. Lee T. J. Taylor P. M. François J.-P. (1995). The anharmonic force field of ethylene, c2h4, by means of accurate ab initio calculations. J. Chem. Phys.103, 2589–2602. 10.1063/1.469681
61
Martin-Drumel M.-A. Spaniol J.-T. Hölzel H. Agúndez M. Cernicharo J. Moth-Poulsen K. et al (2023). Searches for bridged bicyclic molecules in space—norbornadiene and its cyano derivatives. Faraday Discuss.245, 284–297. 10.1039/D3FD00016H
62
McCarthy M. C. Lee K. L. K. Loomis R. A. Burkhardt A. M. Shingledecker C. N. Charnley S. B. et al (2021). Interstellar detection of the highly polar five-membered ring cyanocyclopentadiene. Nat. Astron.5, 176–180. 10.1038/s41550-020-01213-y
63
McClure M. K. Rocha W. R. M. Pontoppidan K. M. Crouzet N. Chu L. E. U. Dartois E. et al (2023). An Ice Age JWST inventory of dense molecular cloud ices. Nat. Astron.7, 431–443. 10.1038/s41550-022-01875-w
64
McGuire B. A. (2022). 2021 census of interstellar, circumstellar, extragalactic, protoplanetary disk, and exoplanetary molecules. Astrophys. J. Suppl. Ser.259, 30. 10.3847/1538-4365/ac2a48
65
McGuire B. A. Burkhardt A. M. Kalenskii S. Shingledecker C. N. Remijan A. J. Herbst E. et al (2018). Detection of the aromatic molecule benzonitrile (c-C6H5CN) in the interstellar medium. Science359, 202–205. 10.1126/science.aao4890
66
McGuire B. A. Loomis R. A. Burkhardt A. M. Lee K. L. K. Shingledecker C. N. Charnley S. B. et al (2021). Detection of two interstellar polycyclic aromatic hydrocarbons via spectral matched filtering. Science371, 1265–1269. 10.1126/science.abb7535
67
McIntosh O. García-Florentino C. Fornaro T. Marabello D. Alberini A. Siljeström S. et al (2024). Undecanoic acid and L-phenylalanine in vermiculite: detection, characterization, and UV degradation studies for biosignature identification on mars. Astrobiology24, 518–537. 10.1089/ast.2023.0088
68
Mehta N. Santra G. Martin J. M. (2023). Is explicitly correlated double-hybrid density functional theory advantageous for vibrational frequencies?Can. J. Chem.101, 656–663. 10.1139/cjc-2022-0277
69
Mendolicchio M. Bloino J. Barone V. (2021). General perturb-then-diagonalize model for the vibrational frequencies and intensities of molecules belonging to abelian and non-abelian symmetry groups. J. Chem. Theory Comput.17, 4332–4358. 10.1021/acs.jctc.1c00240
70
Mills I. M. (1972). Molecular spectroscopy: modern research (academic press, New York), chap. Vibration-rotation structure in asymmetric- and symmetric-top molecules, 115–140.
71
Najibi A. Goerigk L. (2018). The nonlocal kernel in van der waals density functionals as an additive correction: An extensive analysis with special emphasis on the B97M-v and B97M-V approaches. J. Chem. Theory Comput.14, 5725–5738. 10.1021/acs.jctc.8b00842
72
Nazari P. van Gelder M. L. van Dishoeck E. F. Tabone B. van ’t Hoff M. L. R. Ligterink N. F. W. et al (2021). Complex organic molecules in low-mass protostars on Solar System scales - II. Nitrogen-bearing species. Astron. Astrophys.650, A150. 10.1051/0004-6361/202039996
73
Nielsen H. H. (1951). The vibration-rotation energies of molecules. Rev. Mod. Phys.23, 90–136. 10.1103/RevModPhys.23.90
74
Öberg K. I. (2016). Photochemistry and astrochemistry: photochemical pathways to interstellar complex organic molecules. Chem. Rev.116, 9631–9663. 10.1021/acs.chemrev.5b00694
75
Palmer M. H. Ridley T. Hoffmann S. V. Jones N. C. Coreno M. de Simone M. et al (2015a). Interpretation of the vacuum ultraviolet photoabsorption spectrum of iodobenzene by ab initio computations. J. Chem. Phys.142, 134302. 10.1063/1.4916121
76
Palmer M. H. Ridley T. Hoffmann S. V. Jones N. C. Coreno M. de Simone M. et al (2015b). The ionic states of iodobenzene studied by photoionization and ab initio configuration interaction and DFT computations. J. Chem. Phys.142, 134301. 10.1063/1.4916120
77
Papajak E. Zheng J. Xu X. Leverentz H. R. Truhlar D. G. (2011). Perspectives on basis sets beautiful: seasonal plantings of diffuse basis functions. J. Chem. Theory Comput.7, 3027–3034. 10.1021/ct200106a
78
Peeters E. Mackie C. Candian A. Tielens A. G. G. M. (2021). A spectroscopic view on cosmic PAH emission. Acc. Chem. Res.54, 1921–1933. 10.1021/acs.accounts.0c00747
79
Perez R. Parès L. P. Newell R. Robinson S. Bernardi P. Réess J.-M. et al (2017). The supercam instrument on the NASA Mars 2020 mission: optical design and performance. nternational Conf. Space Opt. — ICSO 2016 (SPIE)10562, 744–752. 10.1117/12.2296230
80
Petrignani A. Candian A. (2022). “Chapter 3 - astrochemistry: ingredients of life in space,” in New Frontiers in astrobiology. Editors ThombreR.VaishampayanP. (Elsevier), 49–66. 10.1016/B978-0-12-824162-2.00007-5
81
Piccardo M. Penocchio E. Puzzarini C. Biczysko M. Barone V. (2015). Semi-experimental equilibrium structure determinations by employing b3lyp/snsd anharmonic force fields: validation and application to semirigid organic molecules. J. Phys. Chem. A119, 2058–2082. 10.1021/jp511432m
82
Pietropolli Charmet A. Ceselin G. Stoppa P. Tasinato N. (2022). The spectroscopic characterization of halogenated pollutants through the interplay between theory and experiment: application to R1122. Molecules27, 748. 10.3390/molecules27030748
83
Pulay P. Meyer W. Boggs J. E. (1978). Cubic force constants and equilibrium geometry of methane from Hartree–Fock and correlated wavefunctions. J. Chem. Phys.68, 5077–5085. 10.1063/1.435626
84
Puzzarini C. (2013). Rotational spectroscopy meets theory. Phys. Chem. Chem. Phys.15, 6595–6607. 10.1039/C3CP44301A
85
Puzzarini C. (2022). Gas-phase chemistry in the interstellar medium: the role of laboratory astrochemistry. Front. Astron. Space Sci.8. 10.3389/fspas.2021.811342
86
Puzzarini C. Alessandrini S. Bizzocchi L. Melosso M. Rivilla V. M. (2023). From the laboratory to the interstellar medium: a strategy to search for exotic molecules in space. Front. Astron. Space Sci.10. 10.3389/fspas.2023.1211784
87
Puzzarini C. Stanton J. F. Gauss J. (2010). Quantum-chemical calculation of spectroscopic parameters for rotational spectroscopy. Int. Rev. Phys. Chem.29, 273–367. 10.1080/01442351003643401
88
Rajasekhar B. N. Dharmarpu V. Das A. K. Shastri A. Veeraiah A. Krishnakumar S. (2022). A spectroscopic study of benzonitrile. J. Quant. Spectrosc. Radiat. Transf.283, 108159. 10.1016/j.jqsrt.2022.108159
89
Rosnik A. M. Polik W. F. (2014). Vpt2+k spectroscopic constants and matrix elements of the transformed vibrational Hamiltonian of a polyatomic molecule with resonances using Van Vleck perturbation theory. Mol. Phys.112, 261–300. 10.1080/00268976.2013.808386
90
Rudolph H. D. Demaison J. Császár A. G. (2013). Accurate determination of the deformation of the benzene ring upon substitution: equilibrium structures of benzonitrile and phenylacetylene. J. Phys. Chem. A117, 12969–12982. 10.1021/jp408208s
91
Santra G. Sylvetsky N. Martin J. M. L. (2019). Minimally empirical double-hybrid functionals trained against the GMTKN55 database: revDSD-PBEP86-D4, revDOD-PBE-D4, and DOD-SCAN-D4. J. Phys. Chem. A123, 5129–5143. 10.1021/acs.jpca.9b03157
92
Sharma S. Roppel R. D. Murphy A. E. Beegle L. W. Bhartia R. Steele A. et al (2023). Diverse organic-mineral associations in Jezero crater, Mars. Nature619, 724–732. 10.1038/s41586-023-06143-z
93
Sheng M. Silvestrini F. Biczysko M. Puzzarini C. (2021). Structural and vibrational properties of amino acids from composite schemes and double-hybrid dft: hydrogen bonding in serine as a test case. J. Phys. Chem. A125, 9099–9114. 10.1021/acs.jpca.1c06993
94
Sita M. L. Changala P. B. Xue C. Burkhardt A. M. Shingledecker C. N. Lee K. L. K. et al (2022). Discovery of interstellar 2-cyanoindene (2-C9H7CN) in GOTHAM observations of TMC-1. Astrophys. J. Lett.938, L12. 10.3847/2041-8213/ac92f4
95
Tasinato N. Pietropolli Charmet A. Ceselin G. Salta Z. Stoppa P. (2022). In vitro and in silico vibrational–rotational spectroscopic characterization of the next-generation refrigerant HFO-1123. J. Phys. Chem. A126, 5328–5342. 10.1021/acs.jpca.2c04680
96
Vago J. L. Westall F. Coates A. J. Jaumann R. Korablev O. Ciarletti V. et al (2017). Habitability on early mars and the search for biosignatures with the ExoMars rover. Astrobiology17, 471–510. 10.1089/ast.2016.1533
97
Vázquez J. Stanton J. F. (2007). Treatment of fermi resonance effects on transition moments in vibrational perturbation theory. Mol. Phys.105, 101–109. 10.1080/00268970601135784
98
Webb Space Telescope (2024). NGC 604 (NIRCam image). Available at: https://webbtelescope.org/contents/media/images/2024/110/01HQNV4GP6PR6E7ZSJXRRBQQDS.
99
Williford K. H. Farley K. A. Stack K. M. Allwood A. C. Beaty D. Beegle L. W. et al (2018). “Chapter 11 - the NASA mars 2020 rover mission and the search for extraterrestrial life,” in From habitability to life on Mars. Editors CabrolN. A.GrinE. A. (Elsevier), 275–308. 10.1016/B978-0-12-809935-3.00010-4
100
Xu R. Jiang Z. Yang Q. Bloino J. Biczysko M. (2024). Harmonic and anharmonic vibrational computations for biomolecular building blocks: benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers. J. Comput. Chem.45, 1846–1869. 10.1002/jcc.27377
101
Yamamoto R. Ishikawa S. Ebata T. Mikami N. (2000). Vibrational spectra and relaxation of benzonitrile and its clusters using time-resolved stimulated Raman–UV double resonance spectroscopy. J. Raman Spectrosc.31, 295–304. 10.1002/(sici)1097-4555(200004)31:4<295:aid-jrs542>3.0.co;2-y
102
Yang Q. Bloino J. (2022). An effective and automated processing of resonances in vibrational perturbation theory applied to spectroscopy. J. Phys. Chem. A126, 9276–9302. 10.1021/acs.jpca.2c06460
103
Yang Q. Mendolicchio M. Barone V. Bloino J. (2021). Accuracy and reliability in the simulation of vibrational spectra: a comprehensive benchmark of energies and intensities issuing from generalized vibrational perturbation theory to second order (GVPT2). Front. Astron. Space Sci.8. 10.3389/fspas.2021.665232
104
Yatsyna V. Mallat R. Gorn T. Schmitt M. Feifel R. Rijs A. M. et al (2019). Competition between folded and extended structures of alanylalanine (ala-ala) in a molecular beam. Phys. Chem. Chem. Phys.21, 14126–14132. 10.1039/C9CP00140A
105
Zapata Trujillo J. C. Pettyjohn M. M. McKemmish L. K. (2023). High-throughput quantum chemistry: empowering the search for molecular candidates behind unknown spectral signatures in exoplanetary atmospheres. Mon. Not. R. Astron. Soc.524, 361–376. 10.1093/mnras/stad1717
106
Zdanovskaia M. A. Martin-Drumel M.-A. Kisiel Z. Pirali O. Esselman B. J. Woods R. C. et al (2022). The eight lowest-energy vibrational states of benzonitrile: analysis of Coriolis and Darling-Dennison couplings by millimeter-wave and far-infrared spectroscopy. J. Mol. Spectrosc.383, 111568. 10.1016/j.jms.2021.111568
107
Zhao Y. Hochlaf M. Biczysko M. (2021). Identification of DNA bases and their cations in astrochemical environments: computational spectroscopy of thymine as a test case. Front. Astron. Space Sci.8. 10.3389/fspas.2021.757007
Summary
Keywords
vibrational spectra, density functional theory, second-order vibrational perturbation theory, near-infrared, James Webb Space Telescope
Citation
Xu Y and Biczysko M (2024) Toward the identification of cyano-astroCOMs via vibrational features: benzonitrile as a test case. Front. Chem. 12:1439194. doi: 10.3389/fchem.2024.1439194
Received
27 May 2024
Accepted
12 August 2024
Published
03 September 2024
Volume
12 - 2024
Edited by
Jean-Claude Guillemin, UMR6226 Institut des Sciences Chimiques de Rennes (ISCR), France
Reviewed by
Barbara Michela Giuliano, Max Planck Institute for Extraterrestrial Physics, Germany
Nicola Tasinato, Normal School of Pisa, Italy
Updates
Copyright
© 2024 Xu and Biczysko.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Malgorzata Biczysko, malgo.biczysko@gmail.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.