
==== Front
Front Chem
Front Chem
Front. Chem.
Frontiers in Chemistry
2296-2646
Frontiers Media S.A.

1439194
10.3389/fchem.2024.1439194
Chemistry
Original Research
Toward the identification of cyano-astroCOMs via vibrational features: benzonitrile as a test case
Xu and Biczysko
10.3389/fchem.2024.1439194
Xu Yanting

Biczysko Malgorzata *

International Centre for Quantum and Molecular Structures, Department of Physics, College of Sciences, Shanghai University, Shanghai, China
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

*Correspondence: Malgorzata Biczysko, malgo.biczysko@gmail.com
03 9 2024
2024
12 143919427 5 2024
12 8 2024
Copyright © 2024 Xu and Biczysko.
2024
Xu and Biczysko
https://creativecommons.org/licenses/by/4.0/ 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.
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 11cm−1 for all fundamental transitions overall and only 2cm−1 for the C≡N stretching fundamental at 4.49 μm . 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 C≡N overtones at the 2.26 μm and 1.52 μm . The remarkable accuracy of the GVPT2 results opens a pathway for the reliable prediction of spectra for a broader range of cyano-astroCOMs.

vibrational spectra
density functional theory
second-order vibrational perturbation theory
near-infrared
James Webb Space Telescope
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).section-at-acceptanceChemical Physics and Physical Chemistry
==== Body
pmc1 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 μm ) 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: OCN− 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 C ≡ N 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 C ≡ N-astroCOMs) in the ISM. Indeed, benzonitrile (B C≡N , cyano-benzene, C6 H5C ≡ N), 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 C ≡ N group attachment creates a permanent dipole moment, which is null in otherwise “silent” (poly-)aromatic hydrocarbons (PAHs). This breakthrough intensified the search for aromatic C ≡ N-astroCOMs, leading to the detection in the TMC-1 of 1- and 2-cyanonaphthalene ( C10 H7C ≡ N) isomers (McGuire et al., 2021) in 2021 and 2-cyanoindene ( C9 H7C ≡ N) (Sita et al., 2022) in 2022. Following the observation of the six-membered rings, 1- and 2-cyanocyclopentadiene ( C5 H5C ≡ N) (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 C ≡ N-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 C≡N - 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 2165.9cm−1 (4.62 μm ) (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, rs , 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 C≡N (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 rBO and semi-experimental reSE . The most extended study is the 2022 synchrotron investigation of B C≡N in a very broad spectrum range, up to 90′000cm−1 (11.1 eV), which also included the new detection of the gas-phase IR spectra at a resolution of 0.5cm−1 (Rajasekhar et al., 2022). This work, along with the high-resolution far-infrared spectra collected in the 65−−695cm−1 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 X~1A1 , 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 1.5×10−5 Hartree/Bohr and 6.0×10−5 Å, 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): Q=JQ′+K,

where Q and Q′ represent the two sets of mass-weighted normal coordinates. The Duschinsky matrix J 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 reSE (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 rBO 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 reSE 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 reSE 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 Equilibrium parameters of benzonitrile (bond lengths in Å; angles in degrees).

	reSE a	reBO b	CCSD(T)/ANO1 c	rDSD/junTZ	
C1 - C2	1.3968	1.3962	1.4012	1.3993	
C2 - C3	1.3884	1.3882	1.3934	1.3905	
C3 - C4	1.3917	1.3917	1.3967	1.3941	
C1 - C12	1.4347	1.4359	1.4393	1.4360	
C2 -H8	1.0780	1.0803	1.0823	1.0826	
C3 -H7	1.0799	1.0803	1.0824	1.0830	
C4 -H9	1.0800	1.0806	1.0828	1.0832	
C12 -N13	1.1582	1.1583	1.1646	1.1629	
| MAX |d	.	0.0023	0.0064	0.0047	
MAEd	.	0.0007	0.0043	0.0030	
C1 C2C3	119.42	119.52	119.54	119.55	
C2 C3C4	120.27	120.13	120.13	120.17	
C3 C4C6	120.07	120.22	120.19	120.16	
C2 C1C5	120.55	120.50	120.46	120.41	
C1 C2H8	119.77	119.61	119.61	119.64	
C4 C3H7	120.10	120.16	120.16	120.14	
| MAX |d	.	0.16	0.16	0.15	
MAEd	.	0.09	0.09	0.08	
aSemi-experimental reSE 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.

b Theoretical best estimated from reBO 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.

c Theoretical structure obtained at the CCSD(T)/ANO1 level from Zdanovskaia et al. (2022).

d Largest absolute ( | MAX | ) and mean absolute errors (MAEs) of the bond length and angles compared to the reSE 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, Ir 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 ΦJ 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 Rotational spectroscopic constants for the ground vibrational state of benzonitrile.

	Experimenta	CCSD(T)/ANO1b	rDSD/junTZ//B3LYP/SNSD	
A0 (MHz)	5,655.265,428	5,616.	5,638.	
B0 (MHz)	1,546.8757715	1535.	1,541.	
C0 (MHz)	1,214.4040832	1,205.	1,210.	
ΔJ (kHz)	0.0452858	0.0437	0.0433	
ΔJK (kHz)	0.937983	0.923	0.922	
ΔK (kHz)	0.24411	0.241	0.226	
δJ (kHz)	0.01101116	0.0106	0.0106	
δK (kHz)	0.609187	0.593	0.592	
ΦJ (Hz)	0.000002486	0.00000230	0.00000226	
ΦJK (Hz)	0.0015586	0.00149	0.00150	
ΦKJ (Hz)	-0.007863	-0.00761	-0.00769	
ΦK (Hz)	[0.0066915]	0.0066915	0.0067529	
ϕJ (Hz)	0.000001159	0.00000110	0.00000105	
ϕJK (Hz)	0.0007398	0.000755	0.000757	
ϕK (Hz)	0.007480	0.00712	0.00714	
a 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 cm−1 and largest differences of approximately 18cm−1 , 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 ν18 and ν19 . 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 Harmonic wavenumbers (ω,cm−1) and IR intensities (km/mol) compared with reference computed data.

Sym	CCSD(T)/ANO1 a	rDSD/junTZ	
ω	IR int.	ω	IR int.	
A1	
ν1	3,219	3.53	3,221	3.85	
ν2	3,206	5.43	3,208	6.29	
ν3	3,189	0.01	3,189	0.00	
ν4	2,277	6.14	2,263	11.32	
ν5	1,645	0.12	1,653	0.12	
ν6	1,520	9.82	1,530	8.75	
ν7	1,213	0.16	1,221	0.30	
ν8	1,195	0.55	1,202	0.72	
ν9	1,043	2.67	1,051	3.34	
ν10	1,010	0.12	1,018	0.23	
ν11	763	1.38	769	1.49	
ν12	459	0.00	462	0.00	
A2	
ν13	987	0.00	996	0.00	
ν14	860	0.00	866	0.00	
ν15	402	0.00	405	0.00	
B1	
ν16	1,006	0.01	1,009	0.00	
ν17	939	2.81	945	2.83	
ν18	769	48.33	770	56.01	
ν19	696	26.77	679	19.74	
ν20	550	15.42	556	15.05	
ν21	378	0.69	382	0.64	
ν22	143	1.70	144	1.79	
B2	
ν23	3,214	6.27	3,217	6.86	
ν24	3,198	1.94	3,200	2.52	
ν25	1,620	1.33	1,628	1.24	
ν26	1,471	6.71	1,481	6.18	
ν27	1,348	0.96	1,358	1.15	
ν28	1,314	1.97	1,329	1.46	
ν29	1,175	0.26	1,182	0.24	
ν30	1,096	2.96	1,104	3.60	
ν31	629	0.12	633	0.12	
ν32	546	0.14	552	0.24	
ν33	162	4.44	162	4.65	
MAX	-		15	7.7	
MIN	-		-18	-7.0	
MAE	-		6	0.8	
a Ref. (Zdanovskaia et al., 2022).

b 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 cm−1 –4,000 cm−1 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 0.5cm−1 . 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 11cm−1 and the largest positive and negative errors of approximately ± 44cm−1 . 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 12cm−1 with the most intense peak assigned as ν2 , 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 8cm−1 and maximum discrepancies within 31cm−1 . The most important result is extremely good agreement, within 2cm−1 , for the νC≡N , the fingerprint vibration of benzonitrile, which is predicted at 2227cm−1 (4.49 μm ). 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 νC≡N of 2298cm−1 . Indeed, GVPT2 B3LYP-D3/N07D computations in GAUSSIAN 16 lead to a similar result, with νC≡N of 2305cm−1 . Such a huge discrepancy of approximately 70cm−1 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 νC≡N vibrations normally occurring in the broader region 2,200–2,400 cm−1 (Császár and Fogarasi, 1989).

TABLE 4 Fundamental wavenumbers (υ,cm−1) and IR intensities (km/mol) computed at the GVPT2//revDSD/junTZ//B3LYP/SNSD level for benzonitrile compared with reference experimental and computed data.

Sym	Mode description a
Wilson notation
	Ref. a	Ref. b	Current work	Assignments (PED)c	
Expt	Theory d	Expt	Theory e	IR Int. scaled	υ	IR int	
A1	
ν1	20a	ν CH	3,080	3,207	3,106	3,210	2.40	3,148	2.19	ν CH	
ν2	2	ν CH	3,071	3,196	3,066	3,208	6.54	3,078	3.90	ν CH	
ν3	7a	ν CH	3,042	3,178	3,043	3,198	8.46	2,999	1.20	ν CH	
ν4	ν CN	ν CN	2,232	2,332	2,229	2,323	34.02	2,227	7.44	ν CN, ν CC	
ν5	8a	ν CC	1,599	1,641	1,599	1,643	0.62	1,610	0.01	ν CC, β CCH	
ν6	19a	ν CC	1,492	1,582	1,491	1,509	8.39	1,496	6.98	β HCC, β R sym, ν CC	
ν7	13	X-sens	1,191	1,220	1,193	1,214	0.33	1,202	0.08	ν CC, β CCH, β R sym	
ν8	9a	Ring	1,178	1,203	1,178	1,192	0.93	1,185	0.84	β CCH, ν CC	
ν9	18a	β CH	1,027	1,050	1,027	1,045	3.79	1,039	1.90	ν CC, β R sym, β CCH	
ν10	12	β CH	1,001	1,019	1,001	1,010	0.00	1,008	0.14	ω CCCH, τ R asy	
ν11	1	X-sens	769	774	758	769	2.05	764	1.11	β R tri, ν CC	
ν12	6a	X-sens	461	467	459	463	0.00	458	0.00	β R tri, ν CC	
A2	
ν13	17a	γ CH	978	1,002	975	984	0.00	968	0.00	ω CCCH, τ R asy	
ν14	10a	γ CH	848	863	844	849	0.00	841	0.00	ω CCCH	
ν15	16a	ϕ CC	401	410	398	408	0.00	395	0.00	τ R asy, ω CCCH	
B1	
ν16	5	γ CH	987	1,021	1,001	1,008	0.37	977	0.00	β R asy, ν CC	
ν17	17b	γ CH	925	954	926	935	2.50	917	2.86	ω CCCH	
ν18	11	γ CH	758	781	758	773	34.79	746	50.72	ω CCCH, τ R tri, ω CCCC	
ν19	4	ϕ CC	686	706	687	704	38.87	656	23.67	τ R asy, ω CCCH	
ν20	16b	X-sens	548	573	547	566	17.79	531	12.09	β CCN, ω CCCH, τ R asy	
ν21	γ CN	X-sens	381	392	372	388	0.61	371	0.72	τ R asy, β CCN	
ν22	10b	γ CN	172	147	141	147	1.93	138	1.90	ω CCCC, β CCN, τ R tri	
B2	
ν23	20b	ν CH	3,039	3,188	3,093	3,208	6.54	3,093	7.49	ν CH	
ν24	7b	ν CH	3,027	3,204	3,027	3,190	3.95	3,068	0.86	ν CH	
ν25	8b	ν CC	1,584	1,615	1,583	1,614	0.80	1,591	0.89	ν CC	
ν26	19b	ν CC	1,448	1,481	1,448	1,462	6.69	1,454	5.73	β CCH, ν CC	
ν27	14	ν CC	1,337	1,361	1,335	1,351	1.79	1,337	0.25	ν CC, β CCH	
ν28	3	β CH	1,298	1,319	1,288	1,322	0.46	1,298	1.58	β CCH, ν CC	
ν29	9b	β CH	1,163	1,188	1,163	1,176	0.16	1,172	0.13	β CCH, ν CC	
ν30	18b	β CH	1,071	1,105	1,071	1,099	4.35	1,093	2.52	β CCH, ν CC	
ν31	6b	α CCC	629	641		633	0.12	628	0.11	β R sym	
ν32	β CN	β CN	551	570	547	559	0.27	547	0.09	β CCC, ω CCCN	
ν33	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								
a Ref (Kwon et al., 2003).

b Ref (Rajasekhar et al., 2022).

c 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.

d Scaled harmonic computations at the B3LYP/6-311++G(2df,2pd) level of theory.

e Scaled harmonic computations at the B3LYP/aug-cc-pVDZ level of theory.

f Largest positive (MAX), negative (MIN), and mean absolute errors (MAEs) of the benzonitrile fundamental wavenumbers compared with the experiment by Rajasekhar et al. (2022).

g 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 Non-fundamental band wavenumbers (υ,cm−1) and IR intensities (km/mol) computed at the GVPT2//revDSD/junTZ//B3LYP/SNSD level for benzonitrile compared with reference experimental data.

Experiment	Current work	
Assign.	υ	Assign.	υ	IR int.	
2 ν22	282 a	2 ν22	275	0.007	
ν22 + ν33	303a	ν22 + ν33	296	0	
2 ν33	323a	2 ν33	315	0.01	
2 ν32 +2 ν22	1,393b	ν18 + ν19	1,397	0.41	
2 ν18 + ν22	1,688b	ν17 + ν18	1,658	1.37	
ν26 + ν22	1,769b	ν14 + ν17	1,754	1.30	
ν26 +3 ν22	1,816b	ν13 + ν14	1,803	1.69	
ν26 + ν12	1,900b	ν16 + ν17 , ν13 + ν17	1,883	2.83	
ν26 + ν32	1,970b	ν13 + ν16,2ν15	1,942	3.23	
ν32 + ν18	2,178b	2 ν30	2,178	0.02	
a High-resolution IR spectra in the gas phase obtained by Zdanovskaia et al. (2022).

b Gas-phase IR spectra obtained by Rajasekhar et al. (2022).

FIGURE 3 IR spectrum of benzonitrile in the range 500–4,000 cm−1 . Computed stick spectra were broadened by Lorentzian functions with half-width at half-maximum (HWHM) of 2 cm−1 . 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 1650−1980cm−1 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 ν18 band at the 758cm−1 should be noted (computed as 746cm−1 ), while the experimental spectra show two similar intensity peaks in this range, the second one being ν19 at 686cm−1 . However, this discrepancy needs to be linked to the harmonic values, which already predict the intensity of ν18 as twice that of ν19 . 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 4000−6500cm−1 and 6500−10000cm−1 , respectively (the whole 100−9000cm−1 spectrum is also reported in Supplementary Material). The most pronounced bands in 4040−4300cm−1 and 6000−6180cm−1 are related to the combinations of νCH with in-plane ring deformations and 2νCH overtones, respectively. Similarly, at the higher energies, there are combinations of νCH two quanta transitions with in-plane ring deformations ( 7070−7350cm−1 ) and 3 νCH second overtones ( 9000−9200cm−1 ). 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 1cm−1 has been obtained for fundamentals and first overtones of ν6 and ν2 , while a lower accuracy of ν4 of approximately 10cm−1 transfers to errors of 20cm−1 and 26cm−1 for 2 ν4 and 3 ν4 , respectively. This allows us to provide a reliable prediction of νC≡N overtones, with distinct 2 νC≡N transition at 4,426 ±2cm−1 (2.26 μm ) and a significantly weaker second overtone 3 νC≡N at 6,598 ±5cm−1 (1.52 μm ).

FIGURE 4 NIR spectra simulated at the GVPT2//rDSD/junTZ//B3LYP/SNSD level. The spectra were broadened by Lorentzian functions with HWHM of 2 cm−1 .

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 C≡N 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 35′000−90′000cm−1 (4.3–11.1 eV, 0.111–0.286 μm ), 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 μm 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 ( CH2 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 C≡N -astroCOM species, are deferred to subsequent works, within the framework of development of the COSY-ASTRO dataset (COSY COST Action CA21101, 2024).

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.

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.

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
==== Refs
References

Alberini A. Fornaro T. García-Florentino C. Biczysko M. Poblacion I. Aramendia J. (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 38987581
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
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 Saija F. Cassone G. (Cambridge, UK: The Royal Society of Chemistry), 9–36. 10.1039/9781839164354-00009
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
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
Barone V. (2005). Anharmonic vibrational properties by a fully automated second-order perturbative approach. J. Chem. Phys. 122 , 014108. 10.1063/1.1824881
Barone V. Alessandrini S. Biczysko M. Cheeseman J. R. Clary D. C. McCoy A. B. (2021). Computational molecular spectroscopy. Nat. Rev. Methods Prim. 1 , 38–27. 10.1038/s43586-021-00034-1
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 24346191
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 26575928
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 26575928
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 23247893
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 25894724
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 33195078
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. A 127 , 10517–10527. 10.1021/acs.jpca.3c06649 38033327
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 36413540
Becke A. D. (1993). Density-functional thermochemistry. iii. the role of exact exchange. J. Chem. Phys. 98 , 5648–5652. 10.1063/1.464913
Biczysko M. Bloino J. Puzzarini C. (2018a). Computational challenges in astrochemistry. WIREs Comput. Mol. Sci. 8 , e1349. 10.1002/wcms.1349
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 30246841
Bloino J. (2015). A vpt2 route to near-infrared spectroscopy: the role of mechanical and electrical anharmonicity. J. Phys. Chem. A 119 , 5269–5287. 10.1021/jp509985u 25535769
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
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 22462836
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. A 119 , 11862–11874. 10.1021/acs.jpca.5b10067 26580121
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
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
Cernicharo J. Fuentetaja R. Cabezas C. Agúndez M. Marcelino N. Tercero B. (2022). Discovery of five cyano derivatives of propene with the QUIJOTE line survey. Astron. Astrophys. 663 , L5. 10.1051/0004-6361/202244255
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 34666488
Charmet A. P. Stoppa P. Tasinato N. Giorgianni S. Barone V. Biczysko M. (2013). An integrated experimental and quantum-chemical investigation on the vibrational spectra of chlorofluoromethane. J. Chem. Phys. 139 , 164302. 10.1063/1.4825380 24182024
COSY COST Action CA21101 (2024). COSY Confined molecular systems: from a new generation of materials to the stars CA21101 COST Action.
Császár A. G. Fogarasi G. (1989). Scaled quantum mechanical (SQM) force field and theoretical vibrational spectrum for benzonitrile. Spectrochim. Acta A 45 , 845–854. 10.1016/0584-8539(89)80222-3
[Dataset] Herschel Space Observatory (2017). Herschel’s view of the Taurus molecular cloud annotated. Available at: https://www.herschel.caltech.edu/image/nhsc2017-003a.
Duschinsky F. (1937). The importance of the electron spectrum in multi atomic molecules. concerning the franck-condon principle. Acta Physicochim. URSS 7 , 551.
Eigenbrode J. L. Summons R. E. Steele A. Freissinet C. Millan M. Navarro-González R. (2018). Organic matter preserved in 3-billion-year-old mudstones at Gale crater, Mars. Science 360 , 1096–1101. 10.1126/science.aas9185 29880683
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
Erlandsson G. (1954). Microwave spectrum of benzonitrile. J. Chem. Phys. 22 , 1152. 10.1063/1.1740313
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
Fornaro T. Brucato J. R. Poggiali G. Corazzi M. A. Biczysko M. Jaber M. (2020). UV irradiation and near infrared characterization of laboratory mars soil analog samples. Front. Astron. Space Sci. 7 . 10.3389/fspas.2020.539289
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. A 119 , 5313–5326. 10.1021/jp510101y 25474755
Fortenberry R. C. (2024a). Quantum chemistry and astrochemistry: a match made in the heavens. J. Phys. Chem. A 128 , 1555–1565. 10.1021/acs.jpca.3c07601 38381079
Fortenberry R. C. (2024b). A vision for the future of astrochemistry in the interstellar medium by 2050. ACS Phys. Chem. Au 4 , 31–39. 10.1021/acsphyschemau.3c00043 38283789
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. A 125 , 1301–1324. 10.1021/acs.jpca.0c09526 33506678
Frisch M. J. Trucks G. W. Schlegel H. B. Scuseria G. E. Robb M. A. Cheeseman J. R. (2016). Gaussian 16 rev. C.01. Wallingford CT: Gaussian, Inc.
Frost A. J. Lau R. M. Burtscher L. Packham C. Tasker E. J. Rees G. A. (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
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. Life 11 , 568. 10.3390/life11060568 34204233
Fusè M. Mazzeo G. Longhi G. Abbate S. Masi M. Evidente A. (2019). Unbiased determination of absolute configurations by vis-à-vis comparison of experimental and simulated spectra: the challenging case of diplopyrone. J. Phys. Chem. B 123 , 9230–9237. PMID: 31580674. 10.1021/acs.jpcb.9b08375 31580674
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 32104818
Green J. H. S. Harrison D. J. (1976). Vibrational spectra of benzene derivatives—XVII. Benzonitrile and substituted benzonitriles. Spectrochim. Acta A 32 , 1279–1286. 10.1016/0584-8539(76)80166-3
Grimme S. (2011). Density functional theory with london dispersion corrections. WIREs Comput. Mol. Sci. 1 , 211–228. 10.1002/wcms.30
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 20423165
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 21370243
Guélin M. Cernicharo J. (2022). Organic molecules in interstellar space: latest advances. Front. Astron. Space Sci. 9 . 10.3389/fspas.2022.787567
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. A 119 , 6723–6737. 10.1021/acs.jpca.5b03241 26020099
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. A 107 , 10969–10975. 10.1021/jp036391w
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
Lee K. L. K. Changala P. B. Loomis R. A. Burkhardt A. M. Xue C. Cordiner M. A. (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
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 35695165
Lemmens K. Mackie C. Candian A. Lee J. Tielens M. Anouk M. R. (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 37294543
Lide D. R. (1954). Microwave spectrum and structure of benzonitrile. J. Chem. Phys. 22 , 1577–1578. 10.1063/1.1740461
Loru D. Steber A. L. Thunnissen J. M. M. Rap D. B. Lemmens A. K. Rijs A. M. (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
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. A 126 , 3198–3209. 10.1021/acs.jpca.2c01849 35544706
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
Martin-Drumel M.-A. Spaniol J.-T. Hölzel H. Agúndez M. Cernicharo J. Moth-Poulsen K. (2023). Searches for bridged bicyclic molecules in space—norbornadiene and its cyano derivatives. Faraday Discuss. 245 , 284–297. 10.1039/D3FD00016H 37305958
McCarthy M. C. Lee K. L. K. Loomis R. A. Burkhardt A. M. Shingledecker C. N. Charnley S. B. (2021). Interstellar detection of the highly polar five-membered ring cyanocyclopentadiene. Nat. Astron. 5 , 176–180. 10.1038/s41550-020-01213-y
McClure M. K. Rocha W. R. M. Pontoppidan K. M. Crouzet N. Chu L. E. U. Dartois E. (2023). An Ice Age JWST inventory of dense molecular cloud ices. Nat. Astron. 7 , 431–443. 10.1038/s41550-022-01875-w
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
McGuire B. A. Burkhardt A. M. Kalenskii S. Shingledecker C. N. Remijan A. J. Herbst E. (2018). Detection of the aromatic molecule benzonitrile (c-C6H5CN) in the interstellar medium. Science 359 , 202–205. 10.1126/science.aao4890 29326270
McGuire B. A. Loomis R. A. Burkhardt A. M. Lee K. L. K. Shingledecker C. N. Charnley S. B. (2021). Detection of two interstellar polycyclic aromatic hydrocarbons via spectral matched filtering. Science 371 , 1265–1269. 10.1126/science.abb7535 33737489
McIntosh O. García-Florentino C. Fornaro T. Marabello D. Alberini A. Siljeström S. (2024). Undecanoic acid and L-phenylalanine in vermiculite: detection, characterization, and UV degradation studies for biosignature identification on mars. Astrobiology 24 , 518–537. 10.1089/ast.2023.0088 38669050
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
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 34085530
Mills I. M. (1972). Molecular spectroscopy: modern research (academic press, New York), chap. Vibration-rotation structure in asymmetric- and symmetric-top molecules, 115–140.
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 30299953
Nazari P. van Gelder M. L. van Dishoeck E. F. Tabone B. van ’t Hoff M. L. R. Ligterink N. F. W. (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
Nielsen H. H. (1951). The vibration-rotation energies of molecules. Rev. Mod. Phys. 23 , 90–136. 10.1103/RevModPhys.23.90
Öberg K. I. (2016). Photochemistry and astrochemistry: photochemical pathways to interstellar complex organic molecules. Chem. Rev. 116 , 9631–9663. 10.1021/acs.chemrev.5b00694 27099922
Palmer M. H. Ridley T. Hoffmann S. V. Jones N. C. Coreno M. de Simone M. (2015a). Interpretation of the vacuum ultraviolet photoabsorption spectrum of iodobenzene by ab initio computations. J. Chem. Phys. 142 , 134302. 10.1063/1.4916121 25854238
Palmer M. H. Ridley T. Hoffmann S. V. Jones N. C. Coreno M. de Simone M. (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 25854237
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 26598144
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 33780617
Perez R. Parès L. P. Newell R. Robinson S. Bernardi P. Réess J.-M. (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
Petrignani A. Candian A. (2022). “Chapter 3 - astrochemistry: ingredients of life in space,” in New Frontiers in astrobiology. Editors Thombre R. Vaishampayan P. (Elsevier), 49–66. 10.1016/B978-0-12-824162-2.00007-5
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. A 119 , 2058–2082. 10.1021/jp511432m 25648634
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. Molecules 27 , 748. 10.3390/molecules27030748 35164013
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
Puzzarini C. (2013). Rotational spectroscopy meets theory. Phys. Chem. Chem. Phys. 15 , 6595–6607. 10.1039/C3CP44301A 23487179
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
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
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
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
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
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. A 117 , 12969–12982. 10.1021/jp408208s 24160689
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. A 123 , 5129–5143. 10.1021/acs.jpca.9b03157 31136709
Sharma S. Roppel R. D. Murphy A. E. Beegle L. W. Bhartia R. Steele A. (2023). Diverse organic-mineral associations in Jezero crater, Mars. Nature 619 , 724–732. 10.1038/s41586-023-06143-z 37438522
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. A 125 , 9099–9114. 10.1021/acs.jpca.1c06993 34623165
Sita M. L. Changala P. B. Xue C. Burkhardt A. M. Shingledecker C. N. Lee K. L. K. (2022). Discovery of interstellar 2-cyanoindene (2-C9H7CN) in GOTHAM observations of TMC-1. Astrophys. J. Lett. 938 , L12. 10.3847/2041-8213/ac92f4
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. A 126 , 5328–5342. 10.1021/acs.jpca.2c04680 35930010
Vago J. L. Westall F. Coates A. J. Jaumann R. Korablev O. Ciarletti V. (2017). Habitability on early mars and the search for biosignatures with the ExoMars rover. Astrobiology 17 , 471–510. 10.1089/ast.2016.1533 31067287
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
Webb Space Telescope (2024). NGC 604 (NIRCam image). Available at: https://webbtelescope.org/contents/media/images/2024/110/01HQNV4GP6PR6E7ZSJXRRBQQDS.
Williford K. H. Farley K. A. Stack K. M. Allwood A. C. Beaty D. Beegle L. W. (2018). “Chapter 11 - the NASA mars 2020 rover mission and the search for extraterrestrial life,” in From habitability to life on Mars. Editors Cabrol N. A. Grin E. A. (Elsevier), 275–308. 10.1016/B978-0-12-809935-3.00010-4
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 38682874
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
Yang Q. Bloino J. (2022). An effective and automated processing of resonances in vibrational perturbation theory applied to spectroscopy. J. Phys. Chem. A 126 , 9276–9302. 10.1021/acs.jpca.2c06460 36450042
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
Yatsyna V. Mallat R. Gorn T. Schmitt M. Feifel R. Rijs A. M. (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 30869702
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
Zdanovskaia M. A. Martin-Drumel M.-A. Kisiel Z. Pirali O. Esselman B. J. Woods R. C. (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
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
