A Potential New Dimension of Astrophysical Aromatic Functionalization: Phosphaethynyl (CP) Derivatives

arXiv:2610.30509 · astro-ph.GA, astro-ph.SR · Submitted 2026-10-21 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: We're looking at the paper: "A Potential New Dimension of Astrophysical Aromatic Functionalization".

Jocelyn: Here's the gist of the paper. The gist The quantum chemical computations outline a thermodynamically viable, low-temperature, gas-phase formation pathway for phosphaethynyl cyclopropenylidene (c-C3HCP), resulting from cyclopropenylidene (c-C3H2) and the phosphaethynyl (CP) radical.

Vera: Formation Pathways "Quantum chemical computations outline a thermodynamically viable, low-temperature, gas-phase formation pathway for phosphaethynyl cyclopropenylidene (c-C3HCP), resulting from cyclopropenylidene (c-C3H2) and the phosphaethynyl (CP) radical." "The reaction of CP and benzene (C6H6) to produce phosphaethynyl benzene (C6H5CP) is also explored." "A marginally emerged, but tentative, transition state is identified for this reaction using certain quantum chemical approaches but not others." "The barrierless formation of c-C3HCP and possibly barrierless synthesis of C6H5CP stand as examples for future investigations of other CP-substituted aromatic molecules, marking the potential beginning of a new frontier of astronomical aromatic functionalization." Computational Methods The study employed multiple quantum chemical approaches to compute reaction pathways, including: 1.

Jocelyn: High-level coupled cluster theory at the singles, doubles, and perturbative triples level within the explicitly correlation formalism

CCSD(T)-F12b: as implemented in the Molpro2026.1 program for geometry optimizations and harmonic frequency computations with the cc-pVTZ-F12 basis set 2.

Vera: Range-separated double-hybrid density functional theory (DFT) ωB97X-2(TQZ) method, including Grimme’s D3 empirical dispersion correction with Becke–Johnson damping (D3(BJ)) and the minimally augmented triple-ζ ma-def2-TZVPP basis set for stationary point corrections 3.

Title and authors: Vera: So we're looking at this paper titled "A Potential New Dimension of Astrophysical Aromatic Functionalization: Phosphaethynyl (CP) Derivatives," and it really focuses on phosphaethynyl cyclopropenylidene, or c-C3HCP.

Jocelyn: Yeah, so what the main point is that they've done these quantum chemical computations to show a thermodynamically viable, low-temperature gas-phase formation pathway for c-C3HCP from cyclopropenylidene and the phosphaethynyl radical.

Subrahmanyan: That's interesting because it suggests c-C3HCP could be a potentially present astro-molecule that astronomers might be able to search for in phosphorus rich environments, given the detections of both c-C3H2 and the CP radical toward circumstellar shells along with this provided reaction pathway.

Vera: Right, but they also looked at other reactions too, specifically looking at the reaction of CP and benzene to produce phosphaethynyl benzene, or C6H5CP.

Jocelyn: And they found a marginally emerged but tentative transition state for that reaction using certain quantum chemical approaches but not others.

Subrahmanyan: That tentatively identified transition state is what makes it interesting because the paper is suggesting that barrierless formation of c-C3HCP and possibly barrierless synthesis of C6H5CP stand as examples for future investigations of other CP-substituted aromatic molecules, marking a potential beginning of a new frontier in astronomical aromatic functionalization.

Vera: They used quite a few computational methods to get these results, including high-level coupled cluster theory at the singles, doubles, and perturbative triples level within the explicitly correlation formalism called CCSD(T)-F12b as implemented in the Molpro2026 point 1 program for geometry optimizations and harmonic frequency computations with the cc-pVTZ-F12 basis set <ref:2610.30509#pg1>.

Jocelyn: And they also used range-separated double-hybrid density functional theory, specifically the ωB97X-two(TQZ) method including Grimme’s D3 empirical dispersion correction with Becke–Johnson damping, and the minimally augmented triple-ζ ma-def2-TZVPP basis set for stationary point corrections <ref:2610.30509#pg1>.

Subrahmanyan: They also explored various DFT approaches like ωB97X with Grimme’s improved D4 empirical dispersion correction, ωB97M-D4, M06-twoX-D3(zero), and B2PLYP-D3(BJ) with the ma-def2-TZVPP basis set for transition state searches <ref:2610.30509#pg2>.

Vera: When it comes to predicting spectroscopic properties, they computed vibrational frequencies and spectroscopic constants using quartic force fields, which are fourthorder Taylor series approximations of the internuclear Watson Hamiltonian coupled to rotational and vibrational perturbation theory at second order, or VPT2.

Title and authors: Jocelyn: The first DFT quartic force fields were computed in two different ways; one used B3LYP with the N07D basis set, which augments the six-31G basis set with additional dispersion and polarization functions designed to provide higher accuracy for anharmonic vibrational frequency computations of large aromatic molecules.

Subrahmanyan: The analysis of c-C3HCP showed that none of its fundamental vibrational frequencies display any exceptional intensity, implying that unique observation in the infrared would be challenging, although the CP stretch fundamental is computed to be one thousand four hundred fifty-three point one cm−one.

Vera: That means if you were looking for it using infrared astronomy, finding a uniquely intense line would be pretty difficult because none of its fundamental frequencies show any exceptional intensity.

Jocelyn: But they also noted that the rotational spectrum of c-C3HCP should be much more observable with a dipole moment of three point two three D <ref:2610.30509#pg1>.

Subrahmanyan: That dipole moment is important because it suggests potential for radio astronomical searches, which could be a way to detect this molecule, even if the infrared detection is tricky.

Vera: Moving on to the reaction with benzene, they found that the energy of that tentative transition state for C6H5CP formation was computed to be zero point five kcal mol−one emerged with respect to the reactants.

Jocelyn: That zero point five kcal mol−one barrier is considered low enough because it suggests a barrier that could be submerged by minus zero point eight kcal mol−one which points toward the possibility of a reaction happening at very low temperatures.

Subrahmanyan: The paper also states that the addition of CP to c-C3H2 and the subsequent hydrogen atom elimination are each barrierless regardless of the formation pathway, which means those two specific steps are easy to achieve without needing a high energy input.

Vera: That barrierless addition of CP to c-C3H2 is what's enabling this whole process for forming c-C3HCP in the first place.

Jocelyn: And they also explored how these DFT reactions compare to the pure and corrected coupled cluster methods, using ωB97X-two(TQZ)/ma-def2-TZVPP to predict polarizabilities of c-C3H2 and CP for use in collision limit rate coefficient calculations on these reactions <ref:2610.30509#pg1>.

Subrahmanyan: The paper states that while more rigorous kinetic calculations could be performed, previous work conducted for the barrierless additions of CN to pyrene, coronene and ethynylbenzene, as well as C2H to benzonitrile all converge to their respective collision limits at temperatures relevant to cold astrophysical environments.

Vera: So the calculated collision limit rate coefficients provided in this work are likely to sufficiently approximate what would be achieved in more rigorous kinetic studies.

Title and authors: Jocelyn: That gives us confidence that the computational methods used here are reliable enough for predicting how these reactions might actually happen in a cold environment.

Subrahmanyan: Considering the context of other work, like the study on Stardust Galaxies at z>nine this finding connects to how we understand dust-origin transitions behind those UV-bright galaxies.

Vera: It connects by suggesting that phosphorus containing molecules are observable within the ISM and circumstellar shells, as demonstrated by lines for seven out of twenty candidates in that same paper.

Jocelyn: That's a good point because it shows we aren't just looking at carbon and oxygen chemistry anymore; phosphorus is also a player in these environments.

Subrahmanyan: And the fact that they are detecting things like silicon phosphide or tentatively detected cyanophosphaethyne suggests there’s a whole class of molecules out there that we need to keep an eye on as our models get more complex.

Vera: So, to wrap up this paper "A Potential New Dimension of Astrophysical Aromatic Functionalization: Phosphaethynyl (CP) Derivatives," the main implication is that c-C3HCP detection would bring about a new dimension of astrophysical aromatic functionalization in addition to the CN- and CCH-aromatic derivatives we've already seen.

Jocelyn: And it points toward the fact that the CP functionalization of c-C3H2 could be an efficient process in low-temperature astrophysical environments.

Subrahmanyan: I just want to say that while c-C3HCP detection might be difficult in infrared, its radio astronomical features make it an interesting target for future investigations because of that large dipole moment.

Vera: So we've seen how they used the barrierless formation from c-C3H2 and CP to suggest it could be present in phosphorus-rich regions like IRC +ten thousand two hundred sixteen and now we know the potential spectroscopic signatures.

Jocelyn: It’s definitely a path forward for radio astronomy to look for these kinds of molecules.

Subrahmanyan: I just want to say that the work on "A Potential New Dimension of Astrophysical Aromatic Functionalization: Phosphaethynyl (CP) Derivatives" sets up a framework for how we can predict the existence and search strategy in phosphorus-rich regions, which is important for structure formation models.

Vera: So that's it for this paper. We've looked at the potential to form c-C3HCP, its spectroscopic properties, and where it might be found in space.

The paper's summary: Vera: So, this paper is looking at phosphaethynyl cyclopropenylidene, or c-C3HCP, and they're showing that it can form easily in low-temperature gas and how we might actually find it in space.

Jocelyn: Exactly. The core finding is that c-C3HCP forms barrierlessly from its components, cyclopropenylidene and the phosphaethynyl radical.

Subrahmanyan: And this is huge because it suggests that these molecules could be present in phosphorus-rich regions, like IRC plus one thousand twenty sixteen, which we’ve already seen some evidence for in other work.

Vera: Right, but they also looked at a reaction with benzene to make phosphaethynyl benzene. They found a tentative transition state for that reaction and it seems low enough energy-wise that it might happen in cold environments.

Jocelyn: That barrier is just zero point five kcal per mole above the reactants, which means it's not impossibly high, even if they couldn't find a perfectly stable path yet.

Subrahmanyan: The main thing for structure formation models is that this suggests we need to start thinking about phosphorus chemistry in these environments beyond what we currently model.

Vera: And when they looked at the spectroscopy, it turns out none of the main vibrational frequencies for c-C3HCP are super intense, so seeing it in infrared might be tough.

Jocelyn: But they did say that its rotational spectrum should be pretty observable because of that dipole moment, which is three point two three Debye.

Subrahmanyan: That dipole moment makes it a much better target for radio astronomy searches than some other molecules we've seen, like c-C3HCN.

Vera: So the implication for us on the ground is that even if infrared detection is hard, we should focus our search efforts on radio observations for these phosphorus-containing aromatic molecules.

Jocelyn: It’s a clear path forward for observational astronomy to look for these kinds of molecules in cold interstellar clouds and circumstellar shells.

Subrahmanyan: And this work sets up a framework that connects the chemistry we see in labs to the actual molecular populations we expect to find in those distant, low-temperature regions.

The paper's improvements: Tom: So, this part is about how they can make these computational models better for predicting chemistry in space by using different math tools.

Vera: Right, they’re showing how to compare the results from those high-level coupled cluster calculations against the faster density functional theory methods we used earlier.

Jocelyn: They found that when you use range-separated double-hybrid density functional theory with the ma-def2 basis set, it actually aligns quite well with those more rigorous coupled cluster results for rotational constants.

Subrahmanyan: That comparison is important because it tells us which computational tool is reliable for predicting those rotational constants, like the DFT results being more than ten point zero MHz away from the F12-TcCR+TZ values.

Vera: So that means we have a better idea of how much error to expect when we use the DFT approach versus the coupled cluster approach for these larger molecules.

Jocelyn: And they also showed that these DFT methods are good enough to do rovibrational quantum chemical study for bigger CP-PAHs, which is what we need to see for more complex systems.

Subrahmanyan: It suggests that this computational framework can be used to model the behavior of phosphorus containing aromatic molecules in a way that’s consistent across different theoretical levels.

Vera: So the implication is that we can trust these DFT methods more when we’re looking at the rotational properties of molecules with extra atoms, like phosphorus.

Jocelyn: And they also improved how they predict collision limit rate coefficients by using those polarizabilities from the DFT calculations to run those collision limit rate coefficient calculations on c-C3H2 and CP.

Subrahmanyan: That means we can get a more realistic picture of how these reactions actually happen in cold environments, even if we aren't running the full kinetic studies themselves.

Vera: It’s about getting better predictive power for the physical processes that drive this chemistry in space.

Jocelyn: And they also looked at how to handle those larger aromatic species by analyzing the tentative transition state barrier for C6H5CP formation again, but now with a clearer understanding of the energy landscape.

Subrahmanyan: So they’ve refined their method to assess whether a reaction is barrierless or has a slight emerging barrier, giving us more nuanced predictions about CP-substituted benzoid aromatics.

Conclusion: Vera: So we’re wrapping up this look at "A Potential New Dimension of Astrophysical Aromatic Functionalization: Phosphaethynyl (CP) Derivatives." The main point is that c-C3HCP might be a real molecule lurking in phosphorus-rich regions of the cosmos.

Jocelyn: Yeah, because the barrierless formation pathway from cyclopropenylidene and the CP radical seems to make its existence very plausible at low temperatures.

Subrahmanyan: It’s interesting because this connects to how we model structure formation; if these molecules form easily, they could be part of the building blocks we’re trying to understand in those distant galaxies.

Vera: And even though the infrared detection is tough because the vibrational frequencies aren't super intense, that large dipole moment suggests a strong radio astronomical signature.

Jocelyn: That makes it a clear target for future radio searches, which is something my team and I can definitely look into with our pulsar data.

Subrahmanyan: If we can confirm its presence through radio lines, it would give us new constraints on the chemical pathways available in those early universe environments.

Vera: So the takeaway is that even if it’s hard to see directly with a telescope, the predicted radio signal makes this molecule worth chasing with our observational instruments.

Jocelyn: We should keep an eye out for any spectral features matching those predicted rotational constants as we get more data from time-domain surveys.

Subrahmanyan: And for structure formation models, this paper gives us a new chemical ingredient to consider when building those simulations of early cosmic chemistry.

Reace H. J. Willis, Vincent J. Esposito, Thomas H. Speak, Ilsa R. Cooke, Ryan C. Fortenberry

Department of Chemistry, University of British Columbia · Schmid College of Science and Technology, Chapman University Department of Chemistry & Biochemistry, University of Mississippi

astro-ph.GA, astro-ph.SR

Submitted: 2026-10-21

Updated: 2026-10-21

License: http://creativecommons.org/licenses/by/4.0/

The gist: The gist The quantum chemical computations outline a thermodynamically viable, low-temperature, gas-phase formation pathway for phosphaethynyl cyclopropenylidene (c-C3HCP), resulting from

Key concepts

c-C3HCP: Phosphaethynyl cyclopropenylidene is a molecule formed by combining cyclopropenylidene and the phosphaethynyl radical. It is a key focus because its formation pathway suggests it could exist in phosphorus-rich environments like IRC +10216.
Barrierless Formation: This means a chemical reaction can occur without needing to overcome an energy barrier. The paper found that the formation of c-C3HCP from c-C3H2 and CP is barrierless, implying this process is highly efficient at low temperatures.
Spectroscopic Predictions: The study calculated vibrational frequencies and rotational constants for molecules like c-C3HCP. These predictions help researchers understand how these molecules would look in infrared (IR) or radio spectra, guiding observational searches even if the molecule itself is hard to detect in IR.

Summary

The gist The quantum chemical computations outline a thermodynamically viable, low-temperature, gas-phase formation pathway for phosphaethynyl cyclopropenylidene (c-C3HCP), resulting from cyclopropenylidene (c-C3H2) and the phosphaethynyl (CP) radical.

Formation Pathways

"Quantum chemical computations outline a thermodynamically viable, low-temperature, gas-phase formation pathway for phosphaethynyl cyclopropenylidene (c-C3HCP), resulting from cyclopropenylidene (c-C3H2) and the phosphaethynyl (CP) radical."

The reaction of CP and benzene (C6H6) to produce phosphaethynyl benzene (C6H5CP) is also explored.

A marginally emerged, but tentative, transition state is identified for this reaction using certain quantum chemical approaches but not others.

"The barrierless formation of c-C3HCP and possibly barrierless synthesis of C6H5CP stand as examples for future investigations of other CP-substituted aromatic molecules, marking the potential beginning of a new frontier of astronomical aromatic functionalization."

Computational Methods

The study employed multiple quantum chemical approaches to compute reaction pathways, including:

High-level coupled cluster theory at the singles, doubles, and perturbative triples level within the explicitly correlation formalism [CCSD(T)-F12b] as implemented in the Molpro2026.1 program for geometry optimizations and harmonic frequency computations with the cc-pVTZ-F12 basis set
Range-separated double-hybrid density functional theory (DFT) ωB97X-2(TQZ) method, including Grimme’s D3 empirical dispersion correction with Becke–Johnson damping (D3(BJ)) and the minimally augmented triple-ζ ma-def2-TZVPP basis set for stationary point corrections
Various DFT approaches like ωB97X with Grimme’s improved D4 empirical dispersion correction, ωB97M-D4, M06-2X-D3(0), and B2PLYP-D3(BJ) with the ma-def2-TZVPP basis set for transition state searches

Spectroscopic Predictions

The researchers computed vibrational frequencies and spectroscopic constants using quartic force fields (QFFs), fourthorder Taylor series approximations of the internuclear Watson Hamiltonian coupled to rotational and vibrational perturbation theory at second-order (VPT2) The QFF approach involves computing quadratic force constants with CCSD(T)-F12b/cc-pCVTZ-F12 and adding a relativistic correction

The DFT QFFs are computed in two different ways.

"The first uses B3LYP with the N07D basis set, which augments the 6-31G basis set with additional dispersion and polarization functions designed to provide higher accuracy for anharmonic vibrational frequency computations of large aromatic molecules."

The analysis of c-C3HCP showed that none of its fundamental vibrational frequencies display any exceptional intensity, implying that unique observation in the IR would be challenging The CP stretch fundamental is computed to be 1453.1 cm−1

The frequency of this fundamental is computed to be 1453.1 cm−1.

The rotational spectrum of c-C3HCP should be much more observable with a dipole moment of 3.23 D

Aromatic Functionalization Potential

The study explored the reaction of CP and C6H6 to produce phosphaethynyl benzene (C6H5CP) The tentative transition state for this reaction is predicted to be 0.5 kcal mol−1 emerged with respect to the reactants, suggesting a barrier that could be submerged by-0.8 kcal mol−1

The energy of this transition state is computed to be 0.5 kcal mol−1 emerged with respect to the reactants.

The reaction of CP and c-C3H2 is exothermic by 18.4 kcal mol−1, and the addition of CP to c-C3H2 and the subsequent hydrogen atom elimination are each barrierless regardless of the formation pathway

The addition of CP to c-C3H2 and the subsequent hydrogen atom elimination are each barrierless regardless of the formation pathway.

The potential for C6H5CP formation is uncertain because whether an emerged transition state will inhibit the formation of phosphaethynyl substituted benzoid aromatics via CP + aromatic reactions in lowtemperature astrophysical environments is unknown

Conclusion and Outlook

The detection of c-C3HCP would bring about a new dimension of astrophysical aromatic functionalization in addition to the already detected suite of CN- and CCH-aromatic derivatives The rotational spectrum of c-C3HCN has been experimentally measured, it possesses a slightly larger dipole moment than c-C3HCP, and a CN:CP abundance ratio of ∼100 exists in IRC +10216 The work suggests that the CP functionalization of c-C3H2 could be an efficient process in low-temperature astrophysical environments

How it works

The barrierless formation of c-C3HCP occurs through the addition of the CP radical to c-C3H2.

The reaction of CP and C6H6 is likely to act thermodynamically similar, but because of their larger sizes would be much more challenging to acquire trustworthy results for.

The computed rotational constants for c-C3HCP are given in Table 2 and highlight that the extra mass of the P atom render this molecule near-prolate with a κ value of-0.993, slightly higher than that of c-C3HCN at-0.979

"The rotational constants (in MHz: A0 = 34893.2; B0 = 1993.8; C0 = 1884.8) for this molecule are given in Table 2 and highlight that the extra mass of the P atom render this molecule near-prolate with a κ value of-0.993, slightly higher than that of c-C3HCN at-0.979."

The DFT methods utilized here are shown to align well with and potentially surpass the coupled cluster approach, enabling rovibrational quantum chemical study of larger CP-PAHs

Summary

This paper investigates the potential for phosphaethynyl cyclopropenylidene (c-C3HCP) and its derivatives to exist in astrophysical environments by using quantum chemical computations to model their formation pathways and spectroscopic properties It finds that c-C3HCP can form barrierlessly from c-C3H2 and CP, suggesting it could be present in phosphorus-rich regions like IRC +10216 However, its lack of exceptional vibrational intensity makes IR detection difficult, although its large dipole moment suggests potential for radio astronomical searches The study also predicts a tentative barrier for the formation of CP-substituted benzoid aromatics from CP + C6H6, indicating that the viability of this reaction in low-temperature environments is uncertain due to energy coincidence with reactants Ultimately, the paper concludes that while c-C3HCP detection might be difficult in IR, its radio astronomical features make it an interesting target for future investigations

References

Adler, T. B., Knizia, G., & Werner, H.-J. 2007, J. Chem. Phys., 127, 221106

Ag´undez, M., Cernicharo, J., Decin, L., Encrenaz, P., & Teyssier, D. 2014a, Astrophys. J. Lett., 790, L27

Ag´undez, M., Cernicharo, J., & Gu´elin, M. 2007, Astrophys. J., 662, L91

Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1989, Astrophys. J. Suppl. Ser.

Improvements for AI systems

Improved molecular formation prediction for phosphorus-containing aromatic molecules by incorporating barrierless pathways for CP addition to c-C3H2 and C6H6. The system can now reliably predict whether phosphaethynyl cyclopropenylidene (c-C3HCP) or phosphaethynyl benzene (C6H5CP) can form in low-temperature astrophysical environments based on computed transition state energies.
Enhanced spectroscopic prediction for novel aromatic functionalization by providing accurate vibrational frequencies and rotational constants for CP-substituted molecules. The improved system can now generate predicted infrared and radio astronomical signatures, allowing researchers to search for these molecules in observational data despite the low intrinsic intensities of some modes like the c-C3HCP stretch at 1453.1 cm-1.
Improved comparative analysis between different computational methodologies for molecular properties by benchmarking coupled cluster theory against double-hybrid density functional theory. The system can now determine which level of theory is most reliable for predicting rotational constants, such as the DFT results being more than 10.0 MHz from the F12-TcCR+TZ values, guiding future laboratory refinements.
Specific guidance on astrophysical detection strategies by correlating molecular properties with observational constraints. The improved system can now assess whether molecules like c-C3HCP are more likely to be detected via radio astronomy due to their large dipole moments (3.23 D) compared to those like c-C3HCN (3.46 D).
Refined assessment of reaction viability for larger aromatic species by analyzing tentative transition state barriers. The system can now evaluate the barrierless nature or the existence of a tentative and slightly emerged transition state corresponding to the addition of CP to C6H6, providing a more nuanced prediction about whether CP-substituted benzoid aromatics are present in circumstellar media.

Abstract

Quantum chemical computations outline a thermodynamically viable, low-temperature, gas-phase formation pathway for phosphaethynyl cyclopropenylidene (c-C 3 HCP), resulting from cyclopropenylidene (c-C 3 H 2) and the phosphaethynyl (CP) radical. Given the detections of both c-C 3 H 2 and the CP radical toward circumstellar shells along with the provided reaction pathway, c-C 3 HCP is a potentially present astro-molecule which could be searched for in phosphorus rich environments. The reaction of CP and benzene (C 6 H 6) to produce phosphaethynyl benzene (C 6 H 5 CP) is also explored. A marginally emerged, but tentative, transition state is identified for this reaction using certain quantum chemical approaches but not others. Consequently, the formation of CP-substituted benzoid polycyclic aromatic hydrocarbons (PAHs) in astrophysical environments from CP + PAH reactions is possible but currently uncertain. Additionally, explicitly correlated coupled cluster quartic force field and second-order vibrational perturbation theory rovibrational analysis of the vibrational frequencies and rotational constants of c-C 3 HCP are provided, as well as density functional theory predictions of the same for CP-benzene (C 6 H 5 CP) and both isomers of CP-naphthalene (C 1 0 H 7 CP). Our results suggest that IR searches for these molecules would likely be unsuccessful due to the low intensities and coincidence with PAH features, but their large dipole moments could make them feasible to detect with radio astronomy. The barrierless formation of c-C 3 HCP and possibly barrierless synthesis of C 6 H 5 CP stand as examples for future investigations of other CP-substituted aromatic molecules, marking the potential beginning of a new frontier of astronomical aromatic functionalization.

Sources

Aromatic rings in the Central Molecular Zone: Benzonitrile