Next Article in Journal
The Stability of a Mixed-Phase Barium Cerium Iron Oxide under Reducing Conditions in the Presence of Hydrogen
Next Article in Special Issue
Probing the Mechanisms of Inhibitors Binding to Presenilin Homologue Using Molecular Dynamics Simulations
Previous Article in Journal
A Practical Approach to Bicyclic Carbamoyl Pyridones with Application to the Synthesis of HIV-1 Integrase Strand Transfer Inhibitors
Previous Article in Special Issue
Exploration of the Product Specificity of chitosanase CsnMY002 and Mutants Using Molecular Dynamics Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental FTIR-MI and Theoretical Studies of Isocyanic Acid Aggregates

1
Faculty of Chemistry, University of Wrocław, F. Joliot-Curie 14, 50-383 Wrocław, Poland
2
Department of Chemistry, University of Jyväskylä, FI-40014 Jyväskylä, Finland
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(3), 1430; https://doi.org/10.3390/molecules28031430
Submission received: 6 January 2023 / Revised: 29 January 2023 / Accepted: 31 January 2023 / Published: 2 February 2023
(This article belongs to the Special Issue Recent Advances in Molecular Computational Chemistry)

Abstract

:
Homoaggregates of isocyanic acid (HNCO) were studied using FTIR spectroscopy combined with a low-temperature matrix isolation technique and quantum chemical calculations. Computationally, the structures of the HNCO dimers and trimers were optimized at the MP2, B3LYPD3 and B2PLYPD3 levels of theory employing the 6-311++G(3df,3pd) basis set. Topological analysis of the electron density (AIM) was used to identify the type of non-covalent interactions in the studied aggregates. Five stable minima were located on the potential energy surface for (HNCO)2, and nine were located on the potential energy surface for (HNCO)3. The most stable dimer (D1) involves a weak, almost linear N-H⋯N hydrogen bond. Other structures are bound by a N-H⋯O hydrogen bond or by O⋯C or N⋯N van der Waals interactions. Similar types of interactions as in (HNCO)2 were found in the case of HNCO trimers. Among nine stable (HNCO)3 structures, five represent cyclic forms. The most stable T1 trimer structure is characterized by a six-membered ring formed by three N-H⋯N hydrogen bonds and representing high symmetry (C3h). The analysis of the HNCO/Ar spectra after deposition indicates that the N-H⋯O hydrogen-bonded dimers are especially prevalent. Upon annealing, HNCO trimers were observed as well. Identification of the experimentally observed species relied on previous experimental data on HNCO complexes as well as computed data on HNCO homoaggregates’ vibrational spectra.

1. Introduction

Isocyanic acid (HNCO) is an atmospheric pollutant emitted into the air during different processes such as fossil fuel combustion and biomass burning [1,2,3,4,5]. As a trace gas in ambient air, HNCO poses a health risk for humans since, when inhaled, the molecule participates in protein carbamylation reactions that lead to development of serious diseases [6]. In addition to research related to the chemistry of the atmosphere and the impact of HNCO on human health, properties and photochemical reactivity of isocyanic acid and its isomers have been extensively studied both theoretically and experimentally [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]. Four open-chain isomers were previously identified for CHNO species: fulminic acid (HCNO), isofulminic acid (HONC), cyanic acid (HOCN) and the most stable form, isocyanic acid (HNCO). Teles et al. [14] reported infrared spectra of all these four CHNO isomers together with their isotopomers isolated in low-temperature argon matrices. Theoretically, it has been shown by Shapley and Bacskay [23] that in addition to chain forms, several cyclic and branched CHNO structures of much higher energies are stable. Photolysis of HNCO has also been the subject of extensive studies both in argon and xenon matrices leading to isomerization and photodecomposition processes [8,17,18,26,27].
Although the properties of HNCO and its isomers are relatively well understood, much less data is available on molecular complexes of isocyanic acid. Such studies are useful to obtain information on how spectral properties and chemical reactivity of the molecule change upon intermolecular interactions. Raunier et al. [28] investigated the thermal reactivity of HNCO with water ice and the 1:1 HNCO complex with H2O isolated in an argon matrix using infrared spectroscopy. It was concluded, based on both experimental findings and the results of MP2/6-31G(d,p) calculations, that HNCO acts as a proton donor in the complex with H2O to form a N-H⋯O hydrogen bond. In turn, photo-induced decomposition of formohydroxamic acid (HCONHOH) isolated in solid argon led to the formation of two different forms of the 1:1 HNCO complex with water in which HNCO acted either as a proton donor or proton acceptor [29]. Molecular complexes of HNCO with NH2OH or CH3OH have been identified upon photodecomposition of N-hydroxyurea and acetohydroxamic acid, respectively [30,31]. Keresztes et al. [32] performed photolysis using 220 nm light on 1,2,5- and 1,3,4-oxadiazoles, and HCN⋯HNCO and HCN⋯HOCN complexes were produced upon photolysis. Subsequent photoreactions using excitation with a hydrogen lamp led to the formation of three complexes of isocyanic acid: HCNO⋯HNC, HNCO⋯HNC and HNC⋯HOCN. Two papers on HNCO complexes with dinitrogen and sulfur dioxide isolated in argon matrices have been published [33,34]. In both cases, 1:1 complexes exhibited either N-H⋯N or N-H⋯O hydrogen bonds. Along with them, complexes bound only by van der Waals interactions were detected. Very recently, Zhao et al. [35] reported computational studies on various hydrogen-bonded HNCO complexes, including the HNCO dimer, related to atmospheric chemistry.
Non-covalent interactions are considered to have a significant influence on properties of chemical and biological systems. They also have an impact on atmospheric chemistry, affecting existing reactions and contributing to new reaction channels. Molecular complexes formed in dense environments at low temperature can affect chemistry of interstellar media and some planetary atmospheres. Here we report results of isocyanic acid aggregation in low-temperature argon matrices. The formation of dimers and trimers was experimentally followed by FTIR spectroscopy, and it was thereafter identified with the help of quantum chemical calculations. To our knowledge, such aggregates have not been the subject of experimental studies and could be of interest for better understanding of chemical reactivity of tropospheric gases, especially via intermolecular interactions of N-H-containing molecules. In our studies on HNCO aggregates, an interesting question arose as to which atoms in the acid are the best proton acceptor positions for hydrogen bonds involving the N-H group in the HNCO molecule.

2. Results and Discussion

2.1. Structure and Energetics of HNCO Dimers

Five minima (D1-D5) were found on the potential energy surface for the HNCO dimer at the MP2 and B3LYPD3 levels of theory. At the B3LYP and B2PLYPD3 levels, one of the structures denoted D4 converged to D2. All MP2-computed energy minima structures are shown in Figure 1 together with the atom numbering used in this work. Table 1 shows the values of two topological AIM parameters: the electron density ρ(r) and its Laplacian ∇2ρ(r) at the bond critical points. Table 1 also contains the MP2 calculated values of intermolecular distances and angles. Table 2 presents energetic parameters obtained for HNCO dimers using the MP2, B3LYP, B3LYPD3 and B2PLYPD3 methods. Based on the computed structures shown in Figure 1 and the values of the AIM parameters in Table 1, it was found that three structures (D1, D2 and D4) comprised hydrogen bonds of the N-H⋯N or N-H⋯O type, and all three possessed C1 symmetry. Two other structures (D3 and D5) were connected by van der Waals interactions and had C2h symmetry. The most stable structure with the largest interaction energy (Eint) appeared to be an open dimer D1 with two subunits bound by an almost linear N-H⋯N hydrogen bond, with the N⋯N distance of ca. 3.1 Å and the N-H⋯N angle of ca. 174°.
The two molecular moieties in D1 did not lie in the same plane but were twisted with respect to each other. The MP2 method gave a C3N2N6C7 dihedral angle of 62.0°, whereas B3LYPD3 and B2PLYPD3 predicted the D1 structure to be deviated from planarity by 78.3° and 70.7°, respectively. The D2 and D3 dimers were of almost equal stability compared with each other. The D2 form was an elongated structure in which the N-H group of one molecular moiety interacts with the terminal oxygen atom of the second subunit. The N-H⋯O hydrogen bond formed was close to linearity with the N-H⋯O angle of ca. 170°. The NCO backbones of the two moieties in D2 were twisted by ca. 25° from one another. A cyclic planar D3 structure with C2h symmetry had, according to the AIM results, HNCO subunits connected by two C⋯O van der Waals interactions, and their backbones were anti-parallel to each other. The D4 dimer was also planar and contained, like D2, the N-H⋯O type hydrogen bond. However, in this case, the hydrogen bond appeared strongly bent with the N-H⋯O angle of ca. 131° and the backbones of the two moieties almost perpendicular to each other (the N2C7N6 angle of ca. 115°). From the energy point of view, the least stable D5 dimer of C2h symmetry was characterized by anti-parallel NCO backbones, and hydrogen atoms were oriented away from the center of symmetry of the dimer. The two HNCO moieties were bound by one identified N⋯N van der Waals interaction.
Based on computational data on energetics for the HNCO dimers presented in Table 2, some observations can be made. At all levels of calculations employed, the D1 dimer was the most stable one and was characterized by the highest interaction energy. The three dimers D2, D3 and D4 were of similar energetic stability with respect to each other, and they were only slightly less stable than the global minimum D1. The D5 form was the least energetically favored, as this structure was bound only by van der Waals interactions. Very recently, four HNCO dimers of similar structures optimized at B3LYPD3/aug-cc-pVTZ level were reported [35]. These structures resemble the planar structures D1, D2 and D3 in this study. The two structures D4 and D5 were identified for the first time in this work and were not considered in the previous study [35].
It is interesting to compare the geometry of (HNCO)2 with the results obtained earlier for dimers of the sulfur analogue of isocyanic acid, HNCS [36]. Computational studies revealed three structures to be stable for (HNCS)2 (see Figure S1 in Supplementary Materials). Two of them exhibited N-H⋯N and N-H⋯S hydrogen bonds corresponding to the binding patterns found for D1 and D4 isocyanic acid dimers. Although the mutual arrangement of the two molecular moieties in these dimers is slightly different, the type of interaction is analogous. The other (HNCO)2 structures considered here do not have their counterparts in the case of isothiocyanic acid. There was a cyclic structure similar to structure D3 here with anti-parallel HNCS subunits, but the interaction was through N-H⋯S hydrogen bonds opposite to the C⋯O van der Waals interactions appearing in the HNCO dimer. Moreover, for the HNCS dimers [36] the contribution of dispersion to stabilize the system is crucial, and computational methods taking it into account are the right choice, while the results obtained at the more traditional B3LYP level are less reliable. To check whether this is also the case for HNCO dimers, i.e., if the dispersion effects have a significant contribution to the total energy of the interaction, we performed calculations using the B3LYP functional with the 6-311++G(3df,3pd) basis set. These data, shown in Table 2, indicated the highest contribution of dispersion one should expect for the cyclic D3 and D5 dimers characterized by two O⋯C or one N⋯N van der Waals interactions. The numbers in parentheses in Table 2 indicate the percentage contribution of dispersion to the total interaction energy.

2.2. Structure and Energetics of HNCO Trimers

Nine structures were found to be stable at all three levels of theory for HNCO trimers. Their geometries are shown in Figure 2, and their energetic parameters are presented in Table 3. The values of the topological AIM parameters and values of intermolecular distances and angles calculated at the MP2/6-311++G(3df,3pd) level for the trimers are gathered in Table S1 in the Supplementary Materials.
The geometry of (HNCO)3 is, as in the case of dimers, determined by the presence of hydrogen bonds and/or van der Waals interactions. Of the nine structures, the most stable form was a cyclic trimer T1 characterized by three N-H⋯N hydrogen bonds and C3h symmetry. Trimers T2, T3 and T4 exhibited cyclic structures with a Cs, Cs and C3 symmetry, respectively. They had very similar relative energies and were slightly less stable than T1. These trimers included N-H⋯O or both N-H⋯N and N-H⋯O hydrogen bonds. Another cyclic structure T7 characterized by C1 symmetry was the least stable of all hydrogen-bonded forms.
There were also four trimers that represent open, elongated structures with three moieties connected by hydrogen bonds (T5 and T6, both with C1 symmetry) or van der Waals contacts (T9, Cs) or both (T8, C1). Optimized trimer structures indicated the existence of energetic co-operativity in the species. This was also illustrated by the calculated interaction energies. For example, for the most stable trimer (T1), the total interaction energy was about three times the interaction energy associated with a single hydrogen bond found for the D1 dimer. Other trimers with lower symmetry (T2 and T3), representing cyclic structures and containing both N-H⋯N and N-H⋯O hydrogen bonds, had similar energies as T1 and T4 structures. When the cyclic structures were open, the cooperativity effect mediated by the linked hydrogen bonds was lost, and the trimer structures appeared less stable and more weakly bound.

2.3. Results of Matrix Isolation FTIR Studies

Figure 3 presents the selected regions of the infrared spectra of the HNCO/Ar matrices obtained upon deposition at 15 K (10 K for measurement) at different sample concentrations. Trace (e) in this figure shows the corresponding spectral ranges of the difference spectrum obtained by subtracting the HNCO/Ar = 1/12,000 spectrum from the HNCO/Ar = 1/2000 spectrum. When the HNCO/Ar mixtures were deposited on a cold window, the HNCO monomer bands appeared in the spectra, as previously described [12,14]. Note, that in the νNH stretching region of HNCO, a doublet is present due to the Fermi resonance [14]. In addition, even at high dilutions of the gaseous mixtures, several other absorptions were observed. The intensity of these bands increased relative to the monomer absorptions as the concentration of isocyanic acid in the matrix increased.
Analyzing the spectra after the subsequent annealing processes of the matrices at 26, 28, 30, 33 and 35 K (the highest temperature possible for argon matrices annealing) revealed interesting changes in the bands observed. The temperature increase in the range of 26–30 K did not cause significant changes in the spectra, apart from a slight general increase in the observed intensities of the bands. More pronounced changes were observed upon annealing the matrix at 33 and 35 K. Figure 4 shows selected regions of the spectra of the HNCO/Ar = 1/6000 matrices obtained after deposition at 15 K and subsequent 10 min annealing at 33 K. The trace (c) in Figure 4 shows the corresponding ranges of the difference spectrum obtained by subtracting the spectrum after deposition from the spectrum after annealing at 33 K.
As can be seen in Figure 4 (trace (b)) there are bands which were not present or were relatively weak in the spectrum recorded after deposition (see trace (a)). Moreover, these additional bands grew in intensity when the matrix was annealed up to 33 K. Therefore, they must be related to dimers or higher aggregates of HNCO. Additionally, these bands can be grouped into two sets since they behave differently when the matrix temperature increases from 33 K to 35 K (see Figure S2 in Supplementary Materials).
The group of bands denoted as A was characterized by a slight increase in intensity in the temperature range 26–30 K, with a strong intensity increase upon annealing at 33 K and a decrease in intensity upon annealing at 35 K. On the other hand, the intensities of the B-type bands increased over the entire temperature range, but the increase was much stronger between 33 and 35 K. This different behavior of the A and B bands allowed us to assign them to the HNCO dimer and HNCO trimer structures, respectively.
In Table 4, the selected wavenumber shifts calculated for the five most plausible dimer structures using B3LYPD3 and MP2 methods are compared to the corresponding experimental values. Those obtained at the B2PLYPD3/6-311++G(3df,3pd) level are listed in Table S2 in the Supplementary Materials. The computed intensities of the bands are shown in parentheses in these Tables. Analogous data for HNCO trimers, for the νNH spectral region, are presented in Table S3 in the Supplementary Materials. These computational data were used to support the assignment of the HNCO dimer and trimer species in the following.
The spectral features of group A included several bands located 3519.5, 3465.5, 3462.5, 3452.0 and 3448.5 cm−1 (νNH); 2275.0, 2261.5 and 2257.5 cm−1asNCO); 788.0, 764.5 and 729.0 cm−1 (δHNC) and 599.0 and 560.5 cm−1 (δNCO). These band positions and vibrational shifts in relation to the bands of the HNCO monomer, shown in the right panel of Table 4, were in relatively good agreement with the theoretically predicted values for the D2 and D4 dimers. As shown in Figure 1, the D2 and D4 dimers had very similar geometries, both characterized by an elongated structure involving an N-H⋯O hydrogen bond. However, based on obtained experimental and computational results, it is difficult to unequivocally distinguish between these two (HNCO)2 forms.
Another set of bands in the spectra of the annealed matrices (group B) included those observed at 3430.0, 3418.0 and 3353.0 cm−1 (νNH), 720.0 (δHNC) and 555.5 cm−1 (δNCO). Their vibrational shifts relative to the positions of the corresponding HNCO monomer bands were equal to −81.5, −93.5, −158.5, −50.0 and −18.0 cm−1, respectively. Formation of trimers, which is expected to increase when annealing the matrix sample at elevated temperatures, most probably takes place by attaching another HNCO monomer to an existing dimer in the matrix based on increasing mobility of molecules in the less restrictive lattice of solid argon. The same process can also result from mobilising isolated HNCO monomers in three-body collisions, but this is less probable compared to the monomer–dimer interaction.
From the structures shown in Figure 2, while many trimer structures can be formed, the two elongated aggregates T5 and T6 appear to be the most advantageous. Even though these were not the most stable trimer structures computed, the T5 and T6 trimers seem to be the most plausible candidates based on the local or global mobilisation of HNCO molecules in the matrix. Moreover, since the most stable HNCO dimers identified here were linear structures, they will most likely form another hydrogen bond by adding an additional HNCO subunit. Of the two trimer structures, the theoretical shifts in the νNH range calculated for T5 reproduced the experimental results well, and this complex could be the plausible candidate to explain the observed spectral changes upon high temperature annealing cycles.

3. Experimental and Computational Details

3.1. Matrix Isolation FTIR Studies

Cyanuric acid powder (Acros Organics, Geel, Belgium, 98%) was heated in a quartz vessel under vacuum to T = 450 °C, and the resulting gaseous isocyanic acid was collected in a liquid nitrogen trap. Next, it was passed several times through P2O5 to remove water traces and was finally stored in a 250 mL glass bulb. The gaseous mixtures were prepared by mixing HNCO with argon (Linde, 5.0, Dublin, Ireland) in a container in a stainless-steel vacuum system. Pressure of the gas mixture and the deposition rate were controlled by piezo transducers (model 902B, MKS Instruments, Andover, MA, USA, accuracy ±1%) installed in the deposition line. Typically, a rate of 2 mbar/min was used. The amount of the deposited mixture varied depending on the concentration of HNCO/Ar from 20 to 120 mbar. Low temperature was obtained using a closed cycle helium cryostat (APD-Cryogenics, Macungie, PA, USA) and measured by a silicon diode sensor coupled with the digital controller (Scientific Instruments, model 9650-1, West Palm Beach, FL, USA). The matrices were deposited onto a CsI window maintained at 15 K. Infrared spectra were collected after cooling the sample to 10 K. Spectra were recorded (averaged over 128 scans) in transmission mode by a Bruker IFS 66 Fourier Transform spectrometer (Bruker Optik, Ettlingen, Germany) equipped with a liquid nitrogen cooled MCT detector with 0.5 cm−1 resolution.

3.2. Computational Methods

Computational studies for HNCO dimers were performed using the Gaussian16 program package [37]. Structures of the HNCO monomer, dimers and trimers were optimized at the MP2 [38,39,40,41], B3LYP [42,43,44,45], B3LYPD3 [46,47] and B2PLYPD3 [48,49,50] levels of theory using the 6-311++G(3df,3pd) [51,52] basis set. The above methods were chosen because the dispersion interaction is expected to play an important role in the stabilization of the studied complexes. The Cartesian coordinates of optimized HNCO aggregates are provided in Table S4 in Supplementary Materials. Energy optimization of the dimer and trimer structures was conducted with the Boys–Bernardi full counterpoise method by Dannenberg [53,54]. Energies of the optimized structures were further refined by undertaken single-point calculations using the CCSD(T) [55] approach. The interaction energies Eint were estimated by subtracting the energies of the isolated monomers with the frozen geometry from the energy of the complexes. The topological analysis of electron density (Atoms in Molecules, AIM) [56] was performed at the MP2/6-311++G(3df,3pd) level using the AIMAll program (Version 19.10.12, Professional) [57] on the HNCO dimers and trimers to recognize the type of interaction present in the studied aggregates.
Vibrational wavenumbers and intensities were computed at the same levels using the harmonic approximation to confirm that the optimized geometries corresponded to the minima on the potential energy surfaces and to support the analysis of experimental spectra. Spectral shifts upon complexation were obtained as the difference between the aggregate and monomer wavenumbers of the respective modes. To verify the quality of the computational approach for the vibrational data adopted here, a check with higher level theory was performed on the HNCO monomer. The monomer structure was optimized at the CCSD(T)/6-311++G(3df,3pd) level of theory, and the harmonic frequencies obtained thereafter were amended either with B3LYPD3/6-311++G(3df,3pd) or MP2/6-311++G(3df,3pd) anharmonic vibrational corrections. The vibrational frequencies obtained this way were in satisfactory agreement with existing matrix isolation spectrum [12,14] (see Table S5). This gives a plausible confidence in the computed overall vibrational spectra and vibrational shifts of complexed HNCO species considered in the following to help in assigning the vibrational features appearing in the experimental data. We could estimate the role of the anharmonicity on the fundamental vibrational modes primarily for the HNCO monomer species. Based on this, we could approximate that the anharmonicity effects were rather similar for the HNCO subunits in the complexes. This gave us directions for using harmonic calculations and how to assess them with respect to the experimental data obtained.

4. Conclusions

Results of computational and experimental FTIR matrix isolation studies of structure and infrared spectra of dimers and trimers of HNCO were presented and discussed. The computational methods used (MP2, B3LYPD3 and B2PLYPD3) revealed five or four stable dimer and nine trimer structures. Based on the results of topological electron density analysis (AIM), it was possible to determine the types of non-covalent interactions in the aggregates under study. The HNCO subunits upon complexation were bound either by N-H⋯N or/and N-H⋯O hydrogen bonds or by van der Waals interactions. Experimentally, HNCO dimers exhibiting the N-H⋯O hydrogen bond were observed and assigned in the spectra of the freshly deposited argon matrices as well as upon low temperature annealing cycles (below 30 K). Annealing at higher temperatures (at 33 and 35 K) led to the formation of new species: HNCO trimers. Both experimental and computational studies indicated that HNCO molecules engage in specific intermolecular interactions, leading to notable changes in infrared spectra. Our results indicate that open, elongated dimers and trimers with N-H⋯O hydrogen bonds are preferred under the conditions of matrix isolation. Therefore, the oxygen atom in the HNCO molecule seems to be a better proton acceptor in solid argon. The observation presented in this work is interesting because such interactions can become essential at low temperatures and appear to be important in molecular aggregation and nucleation processes related to both atmospheric chemistry and astrochemistry.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28031430/s1, Figure S1: The optimized structures of the HNCS dimer [36]; Figure S2: Example plots of integrated intensity versus annealing temperature for bands A and B; Table S1: Interatomic distances (Å), angles (degree) and electron density parameters of the intermolecular bond critical points BCP (Atomic units) and ring critical points RCP (Atomic units) of the HNCO trimers computed at the MP2/6-311++G(3df,3pd) level; Table S2: Selected unscaled wavenumber shifts (cm−1) calculated for the HNCO dimers using the B2PLYPD3 method with the 6-311++G(3df,3pd) basis set. The calculated infrared intensities (km mol−1) of the bands are given in parentheses; Table S3: Selected unscaled wavenumber shifts (cm−1) for the νNH region calculated for the HNCO trimers using the MP2, B3LYP and B2PLYPD3 methods with the 6-311++G(3df,3pd) basis set. The calculated infrared intensities (km mol−1) of the bands are given in parentheses; Table S4: The MP2/6-311++G(3df,3pd) computed Cartesian coordinates (Å) of all optimized HNCO aggregates; Table S5: Harmonic unscaled and anharmonic fundamental vibrational wavenumbers of HNCO, in cm−1.

Author Contributions

Conceptualization, J.K. and M.W.; methodology, J.K. and M.W.; software, J.K. and J.L.; validation, J.K., M.W. and J.L.; formal analysis, J.K. and M.W.; investigation, J.K.; resources, M.W. and J.L.; data curation, M.W.; writing—original draft preparation, J.K. and M.W.; writing—review and editing, M.W., J.K. and J.L.; visualization, J.K.; supervision, M.W.; project administration, M.W.; funding acquisition, M.W. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Academy of Finland, grant number 332023, for the work performed at the University of Jyväskylä.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Allocation of computer time from the Wrocław Centre for Networking and Supercomputing (Wrocław, Poland) is gratefully acknowledged. The support of the COST Action CA21101 (COSY) is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Jathar, S.H.; Heppding, C.; Link, M.F.; Farmer, D.K.; Akherati, A.; Kleeman, M.J.; de Gouw, J.A.; Veres, P.R.; Roberts, J.M. Investigating diesel engines as an atmospheric source of isocyanic acid in urban areas. Atmos. Chem. Phys. 2017, 17, 8959–8970. [Google Scholar] [CrossRef]
  2. Perry, R.A.; Siebers, D.L. Rapid reduction of nitrogen oxides in exhaust gas streams. Nature 1986, 324, 657–658. [Google Scholar] [CrossRef]
  3. Wicke, B.G.; Grady, K.A.; Ratcliffe, J.W. Limitations on the rapid reduction of nitrogen oxides in exhaust gas streams. Nature 1989, 338, 492–493. [Google Scholar] [CrossRef]
  4. Lyon, R.K.; Cole, J.A. A reexamination of the RapreNOx process. Combust. Flame 1990, 82, 435–443. [Google Scholar] [CrossRef]
  5. Miller, J.A.; Bowman, C.T. Kinetic modeling of the reduction of nitric oxide in combustion products by isocyanic acid. Int. J. Chem. Kinet. 1991, 23, 289–313. [Google Scholar] [CrossRef]
  6. Leslie, M.D.; Ridoli, M.; Murphy, J.G.; Borduas-Dedekind, N. Isocyanic acid (HNCO) and its fate in the atmosphere: A review. Environ. Sci. Process. Impacts 2019, 21, 793–808. [Google Scholar] [CrossRef] [PubMed]
  7. Kewley, R.; Sastry, K.V.L.N.; Winnewisser, M. The millimeter wave spectra of isocyanic and isothiocyanic acids. J. Mol. Spectrosc. 1963, 10, 418–441. [Google Scholar] [CrossRef]
  8. Jacox, M.E.; Milligan, D.E. Low-Temperature Infrared Study of Intermediates in the Photolysis of HNCO and DNCO. J. Chem. Phys. 1964, 40, 2457–2460. [Google Scholar] [CrossRef]
  9. Steiner, D.A.; Wishah, K.A.; Polo, S.R.; McCubbin, T.K., Jr. Infrared spectrum of isocyanic acid between 465 and 1100 cm−1. J. Mol. Spectrosc. 1979, 76, 341–373. [Google Scholar] [CrossRef]
  10. Drozdoski, W.S.; Baronavski, A.P.; McDonald, J.R. Photodissociation of HNCO at 193 nm. NH(a1∆) internal energy distribution and reaction rate with O2. Chem. Phys. Lett. 1979, 64, 421–425. [Google Scholar] [CrossRef]
  11. Yamada, K. Molecular structure and centrifugal distortion constants of isocyanic acid from the microwave, millimeter wave, and far-infrared spectra. J. Mol. Spectrosc. 1980, 79, 323–344. [Google Scholar] [CrossRef]
  12. Bondybey, V.E.; English, J.H.; Mathews, C.W.; Contolini, R.J. Infrared spectra and isomerization of CHNO species in rare gas matrices. J. Mol. Spectrosc. 1982, 92, 431–442. [Google Scholar] [CrossRef]
  13. Fujimoto, G.T.; Umstead, M.E.; Lin, M.C. Dynamics of CO formation in the photodissociation of HNCO and CH2CO at 193nm. Chem. Phys. 1982, 65, 197–203. [Google Scholar] [CrossRef]
  14. Teles, J.H.; Maier, G.; Hess, B.A., Jr.; Schaad, L.J.; Winnewisser, M.; Winnewisser, B.P. The CHNO Isomers. Chem. Ber. 1989, 122, 753–766. [Google Scholar] [CrossRef]
  15. Ruscic, B.; Berkowitz, J. The H–NCO bond energy and ∆H0f(NCO) from photoionization mass spectrometric studies of HNCO and NCO. J. Chem. Phys. 1994, 100, 4498. [Google Scholar] [CrossRef]
  16. Droz-Georget, T.; Zyrianov, M.; Sanov, A.; Reisler, H. Photodissociation of HNCO: Three competing pathways. Ber. Bunsenges. Phys. Chem. 1997, 101, 469–477. [Google Scholar] [CrossRef]
  17. Pettersson, M.; Khriachtchev, L.; Jolkkonen, S.; Räsänen, M. Photochemistry of HNCO in Solid Xe: Channels of UV Photolysis and Creation of H2NCO Radicals. J. Phys. Chem. A 1999, 103, 9154–9162. [Google Scholar] [CrossRef]
  18. Pettersson, M.; Khriachtchev, L.; Lundell, J.; Jolkkonen, A.S.; Räsänen, M. Photochemistry of HNCO in Solid Xenon: Photoinduced and Thermally Activated Formation of HXeNCO. J. Phys. Chem. A 2000, 104, 3579–3583. [Google Scholar] [CrossRef]
  19. Yu, S.; Su, S.; Dai, D.; Yuan, K.; Yang, X. Vacuum Ultraviolet Photodissociation Dynamics of Isocyanic Acid: The Hydrogen Elimination Channel. J. Phys. Chem. A 2013, 117, 13564–13571. [Google Scholar] [CrossRef]
  20. Yu, S.; Su, S.; Dorenkamp, Y.; Wodtke, A.M.; Dai, D.; Yuan, K.; Yang, X. Competition between Direct and Indirect Dissociation Pathways in Ultraviolet Photodissociation of HNCO. J. Phys. Chem. A 2013, 117, 11673–11678. [Google Scholar] [CrossRef]
  21. Zhang, Z.; Chen, Z.; Huang, C.; Chen, Y.; Dai, D.; Parker, D.H.; Yang, X. Imaging the Pair-Correlated HNCO Photodissociation: The NH(a1∆) + CO(X1Σ+) Channel. J. Phys. Chem. A 2014, 118, 2413–2418. [Google Scholar] [CrossRef] [PubMed]
  22. Holzmeier, F.; Wolf, T.J.A.; Gienger, C.; Wagner, I.; Bozek, J.; Nandi, S.; Nicolas, C.; Fischer, I.; Gühr, M.; Fink, R.F. Normal and resonant Auger spectroscopy of isocyanic acid, HNCO. J. Chem. Phys. 2018, 149, 034308. [Google Scholar] [CrossRef] [PubMed]
  23. Shapley, W.A.; Bacskay, G.B. A Gaussian-2 Quantum Chemical Study of CHNO: Isomerization and Molecular Dissociation Reactions. J. Phys. Chem. A 1999, 103, 6624–6631. [Google Scholar] [CrossRef]
  24. Mladenović, M.; Elhiyani, M.; Lewerenz, M. Electric and magnetic properties of the four most stable CHNO isomers from ab initio CCSD(T) studies. J. Chem. Phys. 2009, 131, 034302. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, C.; Mattila, J.M.; Farmer, D.K.; Arata, C.; Goldstein, A.H.; Abbatt, J.P.D. Behavior of Isocyanic Acid and Other Nitrogen-Containing Volatile Organic Compounds in The Indoor Environment. Environ. Sci. Technol. 2022, 56, 7598–7607. [Google Scholar] [CrossRef]
  26. Milligan, D.E.; Jacox, M.E. Matrix-Isolation Study of the Infrared and Ultraviolet Spectra of the Free Radical NCO. J. Chem. Phys. 1967, 47, 5157. [Google Scholar] [CrossRef]
  27. Maas, M.; Gross, C.; Schurath, U. Vibronic spectra of the matrix-isolated van der Waals complexes NH·N2 and NH·CO. Chem. Phys. 1994, 189, 217–234. [Google Scholar] [CrossRef]
  28. Raunier, S.; Chiavassa, T.; Allouche, A.; Marinelli, F.; Aycard, J.-P. Thermal reactivity of HNCO with water ice: An infrared and theoretical study. Chem. Phys. 2003, 288, 197–210. [Google Scholar] [CrossRef]
  29. Sałdyka, M.; Mielke, Z. Photodecomposition of formohydroxamic acid. Matrix isolation FTIR and DFT studies. Phys. Chem. Chem. Phys. 2003, 5, 4790–4797. [Google Scholar] [CrossRef]
  30. Sałdyka, M. Photodecomposition of N-hydroxyurea in argon matrices. FTIR and theoretical studies. RSC Adv. 2013, 3, 1922–1932. [Google Scholar] [CrossRef]
  31. Sałdyka, M.; Mielke, Z. Photochemistry of acetohydroxamic acid in solid argon. FTIR and theoretical studies. J. Phys. Chem. A 2018, 122, 60–71. [Google Scholar] [CrossRef] [PubMed]
  32. Keresztes, B.; Csonka, I.P.; Lajgút, G.G.; Bazsó, G.; Tarczay, G. Matrix isolation study on the 1:1 hydrogen-bonded complexes [H, C, N]· · · [H, C, N, O]. J. Mol. Struct. 2020, 1219, 128535. [Google Scholar] [CrossRef]
  33. Krupa, J.; Wierzejewska, M.; Lundell, J. Structure and IR spectroscopic properties of HNCO complexes with SO2 isolated in solid argon. Molecules 2021, 26, 6441. [Google Scholar] [CrossRef]
  34. Krupa, J.; Wierzejewska, M.; Lundell, J. Matrix isolation FTIR and theoretical study of weakly bound complexes of isocyanic acid with nitrogen. Molecules 2022, 27, 495. [Google Scholar] [CrossRef] [PubMed]
  35. Zhao, H.; Han, D.; Sun, Y.; Song, X.; Zhang, Y.; Shi, F.; Sheng, X. Interactions between isocyanic acid and atmospheric acidic, neutral and basic species. Comput. Theor. Chem. 2021, 1204, 113384. [Google Scholar] [CrossRef]
  36. Krupa, J.; Wierzejewska, M. Structure of isothiocyanic acid dimers. Theoretical and FTIR matrix isolation studies. Chem. Phys. Lett. 2016, 652, 46–49. [Google Scholar] [CrossRef]
  37. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision C.01; Gaussian, Inc.: Wallingford, CT, USA, 2019. [Google Scholar]
  38. Head-Gordon, M.; Pople, J.A.; Frisch, M.J. MP2 energy evaluation by direct methods. Chem. Phys. Lett. 1988, 153, 503–506. [Google Scholar] [CrossRef]
  39. Head-Gordon, M.; Head-Gordon, T. Analytic MP2 frequencies without fifth-order storage. Theory and application to bifurcated hydrogen bonds in the water hexamer. Chem. Phys. Lett. 1994, 220, 122–128. [Google Scholar] [CrossRef]
  40. Frisch, M.J.; Head-Gordon, M.; Pople, J.A. Semi-direct algorithms for the MP2 energy and gradient. Chem. Phys. Lett. 1990, 166, 281–289. [Google Scholar] [CrossRef]
  41. Sæbø, S.; Almlöf, J. Avoiding the integral storage bottleneck in LCAO calculations of electron correlation. Chem. Phys. Lett. 1989, 154, 83–89. [Google Scholar] [CrossRef]
  42. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  43. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  44. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  45. Vosko, S.H.; Wilk, L.; Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: A critical analysis. Can. J. Phys. 1980, 58, 1200–1211. [Google Scholar] [CrossRef]
  46. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed]
  47. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef] [PubMed]
  48. Grimme, S. Semiempirical hybrid density functional with perturbative second-order correlation. J. Chem. Phys. 2006, 124, 034108. [Google Scholar] [CrossRef]
  49. Schwabe, T.; Grimme, S. Double-hybrid density functionals with long-range dispersion corrections: Higher accuracy and extended applicability. Phys. Chem. Chem. Phys. 2007, 9, 3397–3406. [Google Scholar] [CrossRef]
  50. Biczysko, M.; Panek, P.; Scalmani, G.; Bloino, J.; Barone, V. Harmonic and Anharmonic Vibrational Frequency Calculations with the Double-Hybrid B2PLYP Method: Analytic Second Derivatives and Benchmark Studies. J. Chem. Theory Comput. 2010, 6, 2115–2125. [Google Scholar] [CrossRef]
  51. McLean, A.D.; Chandler, G.S. Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z = 11–18. J. Chem. Phys. 1980, 72, 5639. [Google Scholar] [CrossRef]
  52. Krishnan, R.; Binkley, J.S.; Seeger, R.; Pople, J.A. Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. J. Chem. Phys. 1980, 72, 650. [Google Scholar] [CrossRef]
  53. Boys, S.F.; Bernardi, F. Calculation of small molecular interactions by differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  54. Simon, S.; Duran, M.; Dannenberg, J.J. How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers? J. Chem. Phys. 1996, 105, 11024. [Google Scholar] [CrossRef]
  55. Raghavachari, K.; Trucks, G.W.; Pople, J.A.; Head-Gordon, M. A fifth-order perturbation comparison of electron correlation theories. Chem. Phys. Lett. 1989, 157, 479–483. [Google Scholar] [CrossRef]
  56. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  57. Keith, T.A. AIMAll, Version 19.10.12; TK Gristmill Software, Overland Park, KS, USA. 2019. Available online: aim.tkgristmill.com (accessed on 6 December 2022).
Figure 1. The MP2-optimized structures of the HNCO dimers. The H, N, C and O atoms are shown in white, blue, gray and red, respectively. The positions of the bond (3,−1) and ring (3,+1) critical points derived from AIM calculations are shown by small green and red dots, respectively.
Figure 1. The MP2-optimized structures of the HNCO dimers. The H, N, C and O atoms are shown in white, blue, gray and red, respectively. The positions of the bond (3,−1) and ring (3,+1) critical points derived from AIM calculations are shown by small green and red dots, respectively.
Molecules 28 01430 g001
Figure 2. The MP2-optimized structures of the HNCO trimers. The H, N, C and O atoms are shown in white, blue, gray and red, respectively. The positions of the bond (3,−1) and ring (3,+1) critical points derived from AIM calculations are shown by small green and red dots, respectively.
Figure 2. The MP2-optimized structures of the HNCO trimers. The H, N, C and O atoms are shown in white, blue, gray and red, respectively. The positions of the bond (3,−1) and ring (3,+1) critical points derived from AIM calculations are shown by small green and red dots, respectively.
Molecules 28 01430 g002
Figure 3. Bottom panels: the νNH, νasNCO and δNCO regions in the spectra of HNCO/Ar matrices with 1/12,000 (a), 1/6000 (b), 1/3000 (c) and 1/2000 (d) concentrations. Top panels: the difference spectrum obtained by subtracting the HNCO/Ar = 1/12,000 spectrum from the HNCO/Ar = 1/2000 spectrum (e). The numbers show the red (minus) and blue (plus) wavenumber shifts compared to the HNCO monomer values. Letter M denotes HNCO monomer bands.
Figure 3. Bottom panels: the νNH, νasNCO and δNCO regions in the spectra of HNCO/Ar matrices with 1/12,000 (a), 1/6000 (b), 1/3000 (c) and 1/2000 (d) concentrations. Top panels: the difference spectrum obtained by subtracting the HNCO/Ar = 1/12,000 spectrum from the HNCO/Ar = 1/2000 spectrum (e). The numbers show the red (minus) and blue (plus) wavenumber shifts compared to the HNCO monomer values. Letter M denotes HNCO monomer bands.
Molecules 28 01430 g003
Figure 4. Bottom panels: the νNH, νasNCO and δNCO vibrational spectral regions for HNCO/Ar = 1/6000 after deposition (a, blue line) and after subsequent 10 min annealing at 33 K (b, grey line). Top panels: the difference spectrum (c, pink line) obtained by subtracting the spectrum (a) from the spectrum (b).
Figure 4. Bottom panels: the νNH, νasNCO and δNCO vibrational spectral regions for HNCO/Ar = 1/6000 after deposition (a, blue line) and after subsequent 10 min annealing at 33 K (b, grey line). Top panels: the difference spectrum (c, pink line) obtained by subtracting the spectrum (a) from the spectrum (b).
Molecules 28 01430 g004
Table 1. Interatomic distances (Å), angles (degree) and electron density parameters of the intermolecular bond critical points BCP (Atomic units) and ring critical points RCP (Atomic units) of the HNCO dimers computed at the MP2/6-311++G(3df,3pd) level.
Table 1. Interatomic distances (Å), angles (degree) and electron density parameters of the intermolecular bond critical points BCP (Atomic units) and ring critical points RCP (Atomic units) of the HNCO dimers computed at the MP2/6-311++G(3df,3pd) level.
ComplexGeometric
Parameters 1
AIM Parameters
Interatomic
Distance
AngleBCPρ(r)2ρ(r)
H...YX⋯YX–H...Y
D12.0473.055173.5H1⋯N60.0211+0.0663
D22.0393.037169.7H1⋯O80.0170+0.0698
D3 3.036 C3⋯O80.0071+0.0291
3.036 C7⋯O40.0071+0.0291
(4 at.) 20.0070+0.0299
D42.303 3.062131.2 H1⋯O80.0110+0.0437
D5 3.196 N2⋯N60.0073+0.0257
1 X: N or C; Y: N or O. 2 RCP.
Table 2. BSSE-corrected interaction energies Eint and relative energies ΔE (kJ mol−1) of the HNCO dimers calculated at MP2, B3LYP, B3LYPD3, B2PLYPD3 and CCSD(T) levels with 6-311++G(3df,3pd) basis set.
Table 2. BSSE-corrected interaction energies Eint and relative energies ΔE (kJ mol−1) of the HNCO dimers calculated at MP2, B3LYP, B3LYPD3, B2PLYPD3 and CCSD(T) levels with 6-311++G(3df,3pd) basis set.
DIMERMP2CCSD(T)/
MP2 2
B3LYPD3CCSD(T)/
B3LYPD3 2
B2PLYPD3 1CCSD(T)/
B2PLYPD3 2
B3LYP 1
ΔEEint ΔEEint ΔEEint ΔEEint ΔEEint ΔEEint ΔEEint
D10.00−15.950.00−15.070.00−16.870.00−15.110.00−16.370.00−15.070.00−12.48 (26%)
D22.42−13.361.77−13.651.56−15.111.74−13.821.91−14.281.77−13.690.84−11.43 (24%)
D32.47−13.191.06−13.820.74−15.910.90−14.031.85−14.280.98−13.984.05−8.04 (49%)
D43.14−12.432.45−12.522.16−14.282.34−12.64
D57.92−7.547.66−6.919.17−7.167.63−6.998.88−6.997.63−6.9511.03−0.92 (87%)
1 D4 structure is not energetically stable at B2PLYP and B3LYP levels. 2 CCSD(T) calculations are single-point calculations at the MP2, B3LYPD3 or B2PLYPD3 optimized geometries, respectively.
Table 3. BSSE-corrected interaction energies Eint and relative energies ΔE (kJ mol−1) of the HNCO trimers calculated at MP2, B3LYPD3, B2PLYPD3 and CCSD(T) levels with 6-311++G(3df,3pd) basis set.
Table 3. BSSE-corrected interaction energies Eint and relative energies ΔE (kJ mol−1) of the HNCO trimers calculated at MP2, B3LYPD3, B2PLYPD3 and CCSD(T) levels with 6-311++G(3df,3pd) basis set.
TRIMERMP2CCSD(T)/
MP2 1
B3LYPD3CCSD(T)/
B3LYPD3 1
B2PLYPD3CCSD(T)/
B2PLYPD3 1
ΔEEintΔEEintΔEEintΔEEintΔEEintΔEEint
T10.00−47.900.00−44.760.00−49.700.00−45.300.00−48.730.00−44.88
T23.74−45.643.09−43.960.48−50.793.03−44.762.51−47.692.99−44.21
T34.63−43.843.32−41.912.45−48.063.31−42.503.84−45.593.27−42.20
T46.50−43.255.26−42.912.80−48.825.05−43.754.80−45.765.09−43.21
T515.50−33.0813.70−32.2814.37−35.9613.69−32.5314.83−34.5013.66−32.28
T618.26−29.9415.76−30.6116.35−33.6615.75−31.0217.03−31.9915.70−30.69
T717.61−29.5615.44−28.5516.70−32.1515.47−28.6017.84−30.0615.43−28.55
T819.11−28.6416.70−27.8018.37−31.0216.79−27.6319.08−29.3916.77−27.72
T926.31−20.8523.28−20.7725.84−23.1523.11−21.1026.63−21.3523.12−21.02
1 CCSD(T) calculations are single-point calculations at the MP2, B3LYPD3 or B2PLYPD3 optimized geometries, respectively.
Table 4. Selected unscaled wavenumber shifts (cm−1) of vibrational modes calculated for the HNCO dimers using the B3LYPD3 and MP2 methods with the 6-311++G(3df,3pd) basis set and compared to the experimental results. The calculated infrared intensities (km mol−1) of the bands are given in parentheses.
Table 4. Selected unscaled wavenumber shifts (cm−1) of vibrational modes calculated for the HNCO dimers using the B3LYPD3 and MP2 methods with the 6-311++G(3df,3pd) basis set and compared to the experimental results. The calculated infrared intensities (km mol−1) of the bands are given in parentheses.
B3LYPD3 1MP2 1ModeExp. 2Assignment
D1D2D3D4D5D1D2D3D4D5νΔν
−3511−140−2−3811−15−1−6νNH3519.5+8.0D2
(158)(185)(324)(188)(298)(174)(197)(354)(207)(4)
−165−78 −48 −143−62 −50−6 3465.5
3462.5
−46.0
−49.0 (site)
D2 and D4
(868)(694) (253) (840)(692) (259)(319) 3452.0
3448.5
−59.5
−63.0
711113115121048νasNCO2275.0
2261.5
+16.0
+2.5
D2
D4
(544)(113)(1365)(530)(1319)(554)(241)(1200)(544)(1141)
−8−2 −1 −70 −1 2257.5−2.0 −1.5D2 and D4
(1131)(1658) (1033) (912)(1361) (858)
8119−83698412−162814δHNC788.0+18.0D2
(212)(247)(525)(137)(396)(213)(249)(528)(136)(408)
30−41 −26 25−40 −30 729.0−41.0D2
(283)(266) (455) (309)(286) (462)
5139−1033−55235−1025−3δNCO599.0+25.5D2
(28)(67)(180)(140)(103)(27)(72)(220)(157)(111)
25−6 −4 13−3 −4 560.5−13.0D2
(85)(113) (93) (96)(134) (125)
1 The theoretical shifts were calculated relative to the unscaled harmonic frequencies of HNCO monomer at 3681, 2330, 798 and 576 cm−1 (B3LYPD3) and 3729, 2337, 785 and 571 cm−1 (MP2). 2 The experimental shifts were calculated relative to the corresponding monomer band positions at 3511.5, 2259.0, 770.0 and 573.5 cm−1, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Krupa, J.; Wierzejewska, M.; Lundell, J. Experimental FTIR-MI and Theoretical Studies of Isocyanic Acid Aggregates. Molecules 2023, 28, 1430. https://doi.org/10.3390/molecules28031430

AMA Style

Krupa J, Wierzejewska M, Lundell J. Experimental FTIR-MI and Theoretical Studies of Isocyanic Acid Aggregates. Molecules. 2023; 28(3):1430. https://doi.org/10.3390/molecules28031430

Chicago/Turabian Style

Krupa, Justyna, Maria Wierzejewska, and Jan Lundell. 2023. "Experimental FTIR-MI and Theoretical Studies of Isocyanic Acid Aggregates" Molecules 28, no. 3: 1430. https://doi.org/10.3390/molecules28031430

APA Style

Krupa, J., Wierzejewska, M., & Lundell, J. (2023). Experimental FTIR-MI and Theoretical Studies of Isocyanic Acid Aggregates. Molecules, 28(3), 1430. https://doi.org/10.3390/molecules28031430

Article Metrics

Back to TopTop