Search

The reactivity and mechanistic intricacies of azatrienes in Diels-Alder reactions have been relatively unexplored despite their intriguing potential applications. In this study, we employ Molecular Electron Density Theory to theoretically investigate the hetero-Diels-Alder reaction involving azatrienes with ethyl vinyl ether and allenyl methyl ether. Analysis of Conceptual Density Functional Theory, energetic profiles, and the topological characteristics is conducted to elucidate the reactions. The revealed mechanism manifests as a polar one-step two-stages process under kinetic control. We establish a clear relationship of between the periselectivity, regioselectivity, and stereoselectivity on one hand and the characteristics of the reactions mechanism on the other hand. The influence of weak interactions on reaction activation barriers and bonding evolution are discussed in detail. We demonstrate that substituents enhancing the reverse electron density flux facilitate the feasibility of the reactions. The results lay ground for a meticulous control of the reaction of azatriene in similar synthetic scenarios., Comment: accepted in Journal of Molecular Graphics and Modelling

Among the silicon bearing species discovered in the interstellar medium, SiS and SiO stand out as key tracers due to their distinct chemistry and abundances in interstellar and circumstellar environments. Our objective is to enhance the network of Si- and S-bearing chemical reactions for a gas-grain model in molecular clouds, encompassing both low and high metallicities. We have calculated the energies and rate coefficients for 6 neutral atom-diatom reactions involved in the SiCS triatomic system, with a special focus on the C+SiS and S+SiC collisions. We employ the coupled cluster method with single and double substitutions and a perturbative treatment of triple substitutions (CCSD(T)) refined at the explicitly correlated CCSD(T)-F12 level. With these computational results in conjunction with data from the literature, we construct an extended network of neutral-neutral chemical reactions. We performed time-dependent models employing the Nautilus gas-grain code, setting the gas temperature to 10 K and the density to 2x10$^4$ cm$^{-3}$. The temperature dependence for the reactions involving SiS were modelled using $k(T)=\alpha \left( T/300 \right)^{\beta} \exp{(-\gamma/T)}$. The high-metallicity models significantly boost the SiS production. Higher initial abundances of C, S, and Si, roughly $\sim$2, 190, and 210 times higher, respectively, contribute to this. Around 10$^3$ yr, destruction mechanisms become relevant. The proposed production reaction S + SiC $\rightarrow$ C + SiS, mitigates these effects. By expanding the gas reaction network using a high metallicity model, we derived estimates for the abundances of interstellar molecules. The inclusion of neutral-neutral mechanisms, particularly via Si+HS and S+SiC channels, played a pivotal role in determining SiS abundance. These mechanisms carry a significance on a par with the well-known and fast ion-neutral reactions., Comment: 11 pages, 5 figures, accepted for publication in Astronomy and Astrophysics

In the interstellar cold gas, the chemistry of formaldehyde (H$_2$CO) can be essential to explain the formation of complex organic molecules. On this matter, the massive and energetic protostellar object G331 is still unexplored and, hence, we carried out a comprehensive study of the isotopologues of H$_2$CO and formyl cation (HCO$^+$), and of protonated formaldehyde (H$_2$COH$^+$) through the APEX observations in the spectral window $\sim$159-356~GHz. We employed observational and theoretical methods to derive the physical properties of the molecular gas combining LTE and non-LTE analyses. Formaldehyde was characterized via 35 lines of H$_2$CO, H$_2^{13}$CO, HDCO and H$_2$C$^{18}$O. The formyl cation was detected via 8 lines of HCO$^+$, H$^{13}$CO$^+$, HC$^{18}$O$^+$ and HC$^{17}$O$^+$. Deuterium was clearly detected via HDCO, whereas DCO$^+$ remained undetected. H$_2$COH$^+$ was detected through 3 clean lines. According to the radiative analysis, formaldehyde appears to be embedded in a bulk gas with a wide range of temperatures ($T\sim$20-90 K), while HCO$^+$ and H$_2$COH$^+$ are primarily associated with a colder gas ($T\lesssim$ 30 K). The reaction H$_2$CO+HCO$^+ \rightarrow$ H$_2$COH$^+$ + CO is crucial for the balance of the three species. We used Nautilus gas-grain code to predict the evolution of their molecular abundances relative to H$_2$ which values at time scales $\sim$10$^3$ yr matched with the observations in G331: [H$_2$CO] = (0.2-2) $\times$10$^{-8}$, [HCO$^+$] = (0.5-4) $\times$10$^{-9}$ and [H$_2$COH$^+$] = (0.2-2) $\times$10$^{-10}$. Based on the molecular evolution of H$_2$CO, HCO$^+$ and H$_2$COH$^+$, we hypothesized about the young lifetime of G331, which is consistent with the active gas-grain chemistry of massive protostellar objects., Comment: 24 pages, 8 figures, 8 tables. Accepted for publication in The Astrophysical Journal

The Si+SO$_2$ reaction is investigated to verify its impact on the abundances of molecules with astrochemical interest, such as SiS, SiO, SO and others. According to our results Si($^3$P) and SO$_2$ react barrierlessly yielding only the monoxides SO and SiO as products. No favourable pathway has been found leading to other products, and this reaction should not contribute to SiS abundance. Furthermore, it is predicted that SiS is stable in collisions with O$_2$, and that S($^3$P)+SiO$_2$ and O($^3$P)+OSiS will also produce SO+SiO. Using these results and gathering further experimental and computational data from the literature, we provide an extended network of neutral-neutral reactions involving Si- and S-bearing molecules. The effects of these reactions were examined in a protostellar shock model, using the Nautilus gas-grain code. This consisted in simulating the physicochemical conditions of a shocked gas evolving from $i.$ primeval cold core, $ii.$ the shock region itself, $iii.$ and finally the gas bulk conditions after the passage of the shock. Emphasising on the cloud ages and including systematically these chemical reactions, we found that [SiS/H$_2$] can be of the order of $\sim$ 10$^{-8}$ in shocks that evolves from clouds of $t=1\times 10^6$ yr, whose values are mostly affected by the SiS+O $\longrightarrow$SiO+S reaction. Perspectives on further models along with observations are discussed in the context of sources harbouring molecular outflows., Comment: 10 pages, 7 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

In most cases, excited state quantum phase transitions can be associated with the existence of critical points (local extrema or saddle points) in a system's classical limit energy functional. However, an excited-state quantum phase transition might also stem from the lowering of the asymptotic energy of the corresponding energy functional. One such example occurs in the 2D limit of the vibron model, once an anharmonic term in the form of a quadratic bosonic number operator is added to the Hamiltonian. The study of this case in the broken-symmetry phase was presented in Phys. Rev. A. 81 050101 (2010). In the present work, we delve further into the nature of this excited-state quantum phase transition and we characterize it in the, previously overlooked, symmetric phase of the model making use of quantities such as the effective frequency, the expected value of the quantum number operator, the participation ratio, the density of states, and the quantum fidelity susceptibility. In addition to this, we extend the usage of the quasilinearity parameter, introduced in molecular physics, to characterize the phases in the spectrum of the anharmonic 2D limit of the vibron model and a down-to-earth analysis has been included with the characterization of the critical energies for the linear isomers HCN/HNC.

Isocyanic acid (HNCO) is a simple molecule with a potential to form prebiotic and complex organic species. Using a spectral survey collected with the Atacama Pathfinder EXperiment (APEX), in this work we report the detection of 42 transitions of HNCO in the hot molecular core/outflow G331.512-0.103 (hereafter G331). The spectral lines were observed in the frequency interval $\sim$ 160 - 355 GHz. By means of Local Thermodynamic Equilibrium (LTE) analyses, applying the rotational diagram method, we studied the excitation conditions of HNCO. The excitation temperature and column density are estimated to be $T_{ex}$ = 58.8 $\pm$ 2.7 K and $N$ = (3.7 $\pm$ 0.5) $\times$ 10$^{15}$ cm$^{-2}$, considering beam dilution effects. The derived relative abundance is between (3.8 $\pm$ 0.5) $\times $10$^{-9}$ and (1.4 $\pm$ 0.2) $\times $10$^{-8}$. In comparison with other hot molecular cores, our column densities and abundances are in agreement. An update of the internal partition functions of the four CHNO isomers: HNCO; cyanic acid, HOCN; fulminic acid, HCNO; and isofulminic acid, HONC is provided. We also used the astrochemical code Nautilus to model and discuss HNCO abundances. The simulations could reproduce the abundances with a simple zero-dimensional model at a temperature of 60 K and for a chemical age of $\sim$ 10$^5$ years, which is larger than the estimated dynamical age for G331. This result could suggest the need for a more robust model and even the revision of chemical reactions associated with HNCO., Comment: 17 pages, 8 figures , Accepted in MNRAS

We present an extension of the Hamiltonian of the two dimensional limit of the vibron model encompassing all possible interactions up to four-body operators. We apply this Hamiltonian to the modeling of the experimental bending spectrum of fourteen molecules. The bending degrees of freedom of the selected molecular species include all possible situations: linear, bent, and nonrigid equilibrium structures; demonstrating the flexibility of the algebraic approach, that allows for the consideration of utterly different physical cases with a general formalism and a single Hamiltonian. For each case, we compute predicted term values used to depict the quantum monodromy diagram, the Birge-Sponer plot, the participation ratio. We also show the bending energy functional obtained using the coherent --or intrinsic-- state formalism., Comment: 67 pages, 18 tables and 15 figures

Recent works have shown that the spectroscopic access to highly-excited states provides enough information to characterize transition states in isomerization reactions. Here, we show that the transition state of the bond breaking HCN-HNC isomerization reaction can also be achieved with the two-dimensional limit of the algebraic vibron model. We describe the system's bending vibration with the algebraic Hamiltonian and use its classical limit to characterize the transition state. Using either the coherent state formalism or a recently proposed approach by Baraban et al. [ Science 2015 , 350 , 1338], we obtain an accurate description of the isomerization transition state. In addition, we show that the energy level dynamics and the transition state wave function structure indicate that the spectrum in the vicinity of the isomerization saddle point can be understood in terms of the formalism for excited state quantum phase transitions., Comment: 16 pages, 1 figure, and two tables. Supplementary Material included. Accepted for publication in The Journal of Physical Chemistry A

We emphasize that the completeness of the partition function, that is, the use of a converged partition function at the typical temperature range of the survey, is very important to decrease the uncertainty on this quantity and thus to derive reliable interstellar molecular densities. In that context, we show how the use of different approximations for the rovibrational partition function together with some interpolation and/or extrapolation procedures may affect the estimate of the interstellar molecular column density. For that purpose, we apply the partition function calculations to astronomical observations performed with the IRAM-30m telescope towards the NGC7538-IRS1 source of two N-bearing molecules: isocyanic acid (HNCO, a quasilinear molecule) and methyl cyanide (CH$_3$CN, a symmetric top molecule). The case of methyl formate (HCOOCH$_3$), which is an asymmetric top O-bearing molecule containing an internal rotor is also discussed. Our analysis shows that the use of different partition function approximations leads to relative differences in the resulting column densities in the range 9 to 43\%. Thus, we expect this work to be relevant for surveys of sources with temperatures higher than 300~K and to observations in the infrared., Comment: Accepted for publication in A&A

The torsional Raman spectra of two astrophysically detected isotopologues of dimethyl-ether, ($^{12}$CH$_3$O$^{12}$CH$_3$ and $^{13}$CH$_3$O$^{12}$CH$_3$), have been recorded at room temperature and cooled in supersonic jet, and interpreted with the help of highly correlated ab initio calculations. Dimethyl-ether displays excited torsional and vibrational levels at low energy that can be populated at the temperatures of the star forming regions, obliging to extend the analysis of the rotational spectrum over the ground state. Its spectrum in the THz region is rather complex due to the coupling of the torsional overtones $2\nu_{11}$ and $2\nu_{15}$ with the COC bending mode, and the presence of many hot bands. The torsional overtones are set here at $2\nu_{11}=385.2$~cm$^{-1}$ and $2\nu_{15}=482.0$~cm$^{-1}$ for $^{12}$CH$_3$O$^{12}$CH$_3$, and $2\nu_{11}=385.0$~cm$^{-1}$ and $2\nu_{15}=481.1$~cm$^{-1}$ for $^{13}$CH$_3$O$^{12}$CH$_3$. The new assignment of $2\nu_{11}$ is downshifted around $\sim 10$~cm$^{-1}$ with respect to the literature. All the other (hot) bands have been re-assigned consistently. In addition, the infrared-forbidden torsional fundamental band $\nu_{11}$ is observed here at 197.8~cm$^{-1}$. The new spectral characterization in the THz region reported here provides improved values of the Hamiltonian parameters, to be used in the analysis of the rotational spectra of DME isotopologues for further astrophysical detections., Comment: Accepted for publication in The Astrophysical Journal Supplement Series

We report the first detection and high angular resolution (1.8" $\times$ 1.1") imaging of acetic acid (CH$_3$COOH) and gGg$^{\prime}$--ethylene glycol (gGg$^{\prime}$(CH$_2$OH)$_2$) towards the Orion Kleinmann--Low nebula. The observations were carried out at $\sim$1.3mm with ALMA during the Cycle~2. A notable result is that the spatial distribution of the acetic acid and ethylene glycol emission differs from that of the other O-bearing molecules within Orion-KL. Indeed, while the typical emission of O-bearing species harbors a morphology associated with a "V-shape" linking the Hot Core region to the Compact Ridge (with an extension towards the BN object), that of acetic acid and ethylene glycol mainly peaks at about 2" southwest from the hot core region (near sources I and n). We find that the measured CH$_3$COOH:aGg$^{\prime}$(CH$_2$OH)$_2$ and CH$_3$COOH:gGg$^{\prime}$(CH$_2$OH)$_2$ ratios differ from the ones measured towards the low-mass protostar IRAS 16293--2422 by more than one order of magnitude. Our best hypothesis to explain these findings is that CH$_3$COOH, aGg$^{\prime}$(CH$_2$OH)$_2$ and gGg$^{\prime}$(CH$_2$OH)$_2$ are formed on the icy-surface of grains and then released into the gas-phase, via co-desorption with water, due to a bullet of matter ejected during the explosive event that occurred in the heart of the Nebula about 500-700 years ago., Comment: Accepted for publication in Astronomy and Astrophysics, section 1. Letters to the Editor

We have surveyed a sample of massive star-forming regions located over a range of distances from the Galactic centre for methyl formate, HCOOCH$_{3}$, and its isotopologues H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$. The observations were carried out with the APEX telescope in the frequency range 283.4--287.4~GHz. Based on the APEX observations, we report tentative detections of the $^{13}$C-methyl formate isotopologue HCOO$^{13}$CH$_{3}$ towards the following four massive star-forming regions: Sgr~B2(N-LMH), NGC~6334~IRS~1, W51 e2 and G19.61-0.23. In addition, we have used the 1~mm ALMA science verification observations of Orion-KL and confirm the detection of the $^{13}$C-methyl formate species in Orion-KL and image its spatial distribution. Our analysis shows that the $^{12}$C/$^{13}$C isotope ratio in methyl formate toward Orion-KL Compact Ridge and Hot Core-SW components (68.4$\pm$10.1 and 71.4$\pm$7.8, respectively) are, for both the $^{13}$C-methyl formate isotopologues, commensurate with the average $^{12}$C/$^{13}$C ratio of CO derived toward Orion-KL. Likewise, regarding the other sources, our results are consistent with the $^{12}$C/$^{13}$C in CO. We also report the spectroscopic characterization, which includes a complete partition function, of the complex H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$ species. New spectroscopic data for both isotopomers H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$, presented in this study, has made it possible to measure this fundamentally important isotope ratio in a large organic molecule for the first time., Comment: Accepted for publication in The Astrophysical Journal Supplement Series

Methyl formate in its first torsionally excited state (vt=1 at 131 cm-1) is detected for the first time toward W51 e2. All transitions from excited methyl formate within the observed spectral range are actually detected (82 transitions) and no strong lines are missing. The column density of the excited state is comparable to that of the ground state. For a source size of 7'' we find that Trot = 104 +/- 14 K and N = 9.4 +4.0/-2.8 x 10^16 cm-2 for the excited state and Trot = 176 +/- 24 K and N = 1.7 +.2/-.2 x 10^17 cm-2 for the ground state. Lines from ethyl cyanide in its two first excited states (vt=1, torsion mode at 212 cm-1) and (vb=1, CCN in-plane bending mode at 206 cm-1) are also present in the observed spectrum. However blending problems prevent a precise estimate of its abundance. With regard to the number of lines of excited methyl formate and ethyl cyanide detected in W51 e2, it appears that excited states of large molecules certainly account for a large number of unidentified lines in spectral survey of molecular clouds.