Thionitroxyl Radical (H₂NS) Isomers: Structures, Vibrational Spectroscopy, Electronic States and Photochemistry

Abstract

The thionitroxyl radical (H₂NS) isomers are characterized using advanced ab initio methodologies. Computations are done using standard and explicitly correlated coupled cluster, CASSCF and MRCI approaches in conjunction with large basis sets, extrapolated to the complete basis set (CBS) limit. The lowest electronic states of different isomers are mapped along the stretching coordinates, thereby confirming the existence of the four already known ground state structures, namely H₂NS, H₂SN, cis-HNSH and trans-HNSH. Also, it is shown that only the lowest electronic excited states are stable, whereas the upper electronic states may undergo unimolecular decomposition processes forming H + HNS/HSN or the HN + SH or N + H₂S or S + NH₂ fragments. These data allow an assignment of the deep blue glow observed after reactions between “active nitrogen” and H₂S at the beginning of the XX^th century. For stable species, a set of accurate structural and spectroscopic parameters are provided. Since small nitrogen-sulfur molecular species are of astrophysical relevance, this work may help for identifying the thionitroxyl radical isomers in astrophysical media and in the laboratory.

Introduction

To date, several molecular species containing H, N and S atoms are identified in astrophysical media. These include the Galactic center cloud Sgr B2, cold dark clouds, R Andromedae, circumstellar envelopes of evolved stars, regions of massive star formation, planetary atmospheres, brown dwarfs, stellar atmospheres, and comae (Espin, 1890; Neugebauer et al., 1965; Wöhl, 1971; Gottlieb et al., 1975; Meyer and Roth, 1991; van Dishoeck et al., 1993; Wagenblast et al., 1993; McGonagle et al., 1994; McGonagle and Irvine, 1997; Irvine et al., 2000; Yamamura et al., 2000; Ziurys, 2006; Phuong et al., 2018). The list of identified species contains the NS, SH, NH, and NH₂ radicals as well as the H₂S molecule. Other species are also very likely to be present there, such as HSN/HNS and their ions (Ben Yaghlane et al., 2014; Ajili et al., 2016; Trabelsi et al., 2016). In these media, these molecules and the H, N, and S atoms may undergo complex chemistry that may lead to the formation of larger molecular species, already detected or possibly present. All of them are closely related to the molecule investigated here, i.e., the thionitroxyl radical (H₂NS), since they may be involved in its formation or fragmentation.

At the beginning of the XX^th century, Strutt, Fowler and Vaidya (Strutt and Fowler, 1912; Strutt, 1913; Fowler and Vaidya, 1931) investigated the reactions of “active nitrogen” with hydrogen sulfide (i.e., N + H₂S) producing NS + H₂. They suggested that these products are formed by decomposition of an N.H₂S tetratomic intermediate. These findings were confirmed in 1960 by Westbury and Winkler (1960) who studied the kinetics of such processes in the 85–440°C temperature range. A deep blue glow was observed in the reaction flame, attributed a priori to the N.H₂S long-lived intermediate. In the literature, no solid confirmation of this assignment is found.

In the laboratory, two isomers of the thionitroxyl radical (namely HNSH and H₂SN) and their positively- and negatively-charged ions were detected in 1994 by a combination of collisional activation and neutralization-reionization mass spectrometry, complemented by theoretical computations based on Möller-Plesset (MP2 and MP4) approaches (Nguyen et al., 1994). These computations showed that the most stable [H,H,N,S] species corresponds to H₂NS (C_2v symmetry) and not to the already detected HNSH and H₂SN isomers. In 1998, the most stable [H,H,N,S] form was observed by Habara et al. (1998) using source-modulation and Fourier-transform microwave spectrometers. These authors recorded the rotational spectral lines of the X²B₁ electronic ground state of the H₂NS and D₂NS radicals in the 150–390 GHz region. They provided accurate equilibrium (r_e) and zero-point average structures (r_z), and estimated the 6 fundamental vibrational modes (ν₁₋₆). Theoretically, earlier works revealed the existence of four stable isomers with molecular formula H₂NS, namely H₂NS (C_2v symmetry), trans-HNSH, cis-HNSH and H₂SN (C_s symmetry). In 2010, Pereira et al. (2010) used the CCSD/6-311++G(d,p) technique to study the ¹NH + ²SH → cis-²HNSH → trans-²HNSH → ²H₂NS → ²NS → + H₂ reactive channels. The H₂NS species are found to be intermediates in these reactions. Although this work is performed at a relatively high computational level, it does not provide structural parameters for the four H₂NS isomers with the accuracy required for their detection in the laboratory and in astrophysical surveys. In addition, it does not consider the H + HNS/HSN reactive channels, which are likely to be involved in the formation of the NS radical in the interstellar clouds, since its formation cannot be fully explained by the electron capture and following fragmentation of the protonated NS. The deduced abundances of NS differ by several orders of magnitude from the measured ones in interstellar dark clouds (Millar et al., 1991; McGonagle et al., 1994; Irvine et al., 2000). The involvement of large H-, N- and S-containing molecules is not invoked. However, these compounds may play key roles in these processes, as intermediates. In this context, H₂NS species are most likely present, as stable molecules or intermediates, in the astrophysical media mentioned and should be considered in the chemical models for the study of their physicochemical properties. Nevertheless, the main limitation of the identification of such compounds is the lack of spectroscopic data. The present work provides a set of structural, vibrational and electronic parameters for the [H,H,N,S] molecular species. It should help for their identification in astrophysical media and in the laboratory. We also show that these tetratomics undergo spin-orbit-induced predissociation and photodissociation upon absorbing UV-Vis light. These restrain the astrophysical environments where they can be formed and detected.

Computational Details

We started our computations by mapping the lowest doublet potential energy surface (PES) of the [H,H,N,S] molecular system to locate the stationary points (both minimal-energy and transition state structures). These geometry optimizations were done without symmetry constraints (i.e., in the C₁ point group), using the partially spin restricted coupled cluster (CC) method with a perturbative treatment of triple excitations [RCCSD(T)] (Hampel et al., 1992; Deegan and Knowles, 1994; Knowles et al., 2000) as well as the explicitly correlated version of this method [RCCSD(T)-F12] (Adler et al., 2007; Werner et al., 2007; Knizia and Werner, 2008; Knizia et al., 2009) and complete active space self-consistent field (CASSCF) (Knowles and Werner, 1985; Werner and Knowles, 1985) approaches as implemented in the MOLPRO suite of programs (Werner et al., 2015). The use of CC approaches is justified by the mono configurational nature of the wavefunctions of H₂NS ground state species close to equilibrium as certified by the T1 and D1 diagnostics close to 0.02 and 0.06 in Coupled Clusters computations and by the dominance of a unique configuration in their CASSCF wavefunctions (weight > 0.95). In CASSCF, the active space is formed by the whole set of valence orbitals, in which all the possible valence electron excitations are considered. For the RCCSD(T) and CASSCF computations, we used the aug-cc-pV(X+d)Z (X = D, T, Q) basis sets for sulfur and the aug-cc-pVXZ (X = D, T, Q) for hydrogen and nitrogen atoms (Kendall and Dunning, 1992; Woon and Dunning, 1993; Dunning et al., 2001). The use of tight d functions for sulfur improves the quality of the results, as established in the benchmark computations presented in Refs (BenYaghlane et al., 2013; Ben Yaghlane et al., 2014; Trabelsi et al., 2018), since it allows for a better description of the electron density of this polarizable atom and can provide accurate descriptions of the electronic structures of its compounds. In the RCCSD(T)-F12 calculations, the cc-pVXZ-F12 (X = D, T, Q) basis sets (Peterson et al., 2008) together with the corresponding auxiliary basis sets and density-fitting functions, as implemented in MOLPRO (Klopper, 2001; Weigend, 2002; Hättig, 2005), were used. These optimizations are followed by frequency computations to establish the nature of the stationary points, either minima or saddle points. Subsequently, we extrapolated the geometrical values and energetics to the CBS limit using the two-parameter equation: A = B + C/X³, in conjunction with a least square procedure, where A is the extrapolated quantity, X is the cardinal number of the basis set and B and C are fitting parameters (Helgaker et al., 1997). For the RCCSD(T)/CBS and RCCSD(T)-F12/CBS geometrical parameters, we evaluated the core-valence correlation corrections at the RCCSD(T)/cc-pCVTZ and RCCSD(T)-F12/cc-pCVTZ-F12 levels, respectively. Within a given level, this was achieved by taking the difference between the optimized geometrical parameters obtained from a calculation where all the electrons are correlated and a calculation where the core electrons are kept frozen.

The PESs of the low-lying electronic states of the H₂NS species were mapped along the stretching and out-of-plane torsion coordinates with the remaining internal coordinates fixed at their respective equilibrium ground state values. These computations were carried out with the CASSCF approach followed by the internally contracted multi-reference configuration interaction (MRCI) method (Werner and Knowles, 1988; Shamasundar et al., 2011), where the H, N, and S atoms were described with the aug-cc-pV(5+d)Z basis set. In the MRCI procedure, all single and double electron excitations from the configurations of the CASSCF wavefunction having a CI coefficient larger than 0.005 were performed. The resulting MRCI wavefunctions are composed of ~1.7 × 10⁷ contracted (~1 × 10⁹ uncontracted) configurations per C_s symmetry. Furthermore, we optimized the equilibrium geometries of the long-lived and stable electronic excited states at the CASSCF/aug-cc-pV(X+d)Z (X= T, Q, 5, CBS) level. The minimal-energy nature of these structures is indicated by all positive frequencies.

On the Electronic States of H₂NS

For the thionitroxyl molecule, as established in the literature, our computations confirm the existence of only four stable forms in the ground state PES. Unlike H₂SN, which is located in the lowest ²A′ PES, H₂NS, cis-HNSH and trans-HNSH are found in the lowest ²A″ PES. These equilibrium structures are displayed in Figure 1. For each of them Figure 2 displays the corresponding singly occupied molecular orbital (SOMO), the next-to-SOMO (SOMO-1) and the lowest unoccupied molecular orbital (LUMO) as computed at the CASSCF/aug-cc-pVTZ level. This figure shows that the SOMOs of H₂NS, cis-HNSH and trans-HNSH are π^* orbitals, while the SOMO-1s are of σ type. All these orbitals are mainly located on the NS bond. Whereas, the SOMO-1 and the SOMO of H₂SN correspond to π orbitals mainly located on the nitrogen atom. For all species, the LUMOs are σ^* SN bonds.

Figure 1

Figure 2

Figures 3–6 present one-dimensional cuts of the 6D-PESs of the doublet and quartet electronic states of H₂NS, cis-HNSH, trans-HNSH and H₂SN along the stretching coordinates, with the remaining internal coordinates kept fixed at their ground state values. We give also the evolution along the out-of-plane / torsion angle. These potentials are computed at the CASSCF/MRCI/aug-cc-pV(5+d)Z level of theory. In these cuts, energies are given with respect to the H₂NS ground state energy at equilibrium. These potentials correlate adiabatically to the HNS/HSN + H or NH + HS or NH₂ + S or H₂S + N dissociation limits at internuclear separations. Table 1 lists the electronic states of H₂NS species located in the 0–6.5 eV internal energy range and gives also their CASSCF/MRCI/aug-cc-pV(5+d)Z vertical electronic excitation energies and dominant electronic configurations, as quoted at the equilibrium geometry of the corresponding electronic state.

Figure 3

Figure 4

Figure 5

Figure 6

For each [H,H,N,S] isomer, Figures 3–6 reveal that the two lowest electronic states have deep potential wells along the stretching coordinates. Besides, the ground electronic states are well-separated from the upper ones in the molecular region and free from interactions that could lead to mixing of their wavefunctions. Therefore, these electronic states are monoconfigurational in nature and can be accurately described by the monoconfigurational post-Hartree-Fock methods used in the present investigation (cf. infra).

Table 1 shows that the first excited states are obtained upon SOMO-1 → SOMO electron promotion. In the Franck-Condon region accessed from the corresponding ground states, the first excited states correspond to 1²B₂, 1²A′, 1²A′, and 1²A″ space symmetry species for H₂NS, cis-HNSH, trans-HNSH, and H₂SN, respectively. They are located at 1.15 – 1.5 eV for H₂NS(1²B₂) and H₂SN(1²A″) and at 3–3.4 eV for cis-HNSH(1²A′) and trans-HNSH(1²A′) with respect to their corresponding electronic ground state. At equilibrium, they possess non-planar configurations with C_s symmetry for H₂NS and H₂SN and C₁ symmetry for HNSH. For cis-HNSH and trans-HNSH a unique minimum is found in the 2²A state PES, whereas both cis and trans forms exist in the ground state potential. Moreover, the lowest electronic excited states interact with the other states for non-planar configurations. For instance, Figure 4 shows that the local minimum in the HNSH excited PES is due to an avoided crossing between the lowest doublets at τ ~100°. In this region the wavefunctions of both states are mixed and acquire multiconfigurational character. Such states should be investigated using multi configurational approaches (cf. infra).

For the low-lying doublet excited states of Table 1, we calculated the oscillator strengths from the vertical excitation energies and transition moments evaluated at the CASSCF/MRCI/aug-cc-pV(5+d)Z level. The calculated values show that rather intense transitions are expected from the X²B₁ to the 2²B₁ state of H₂NS and from the ground to the 2²A″ states of trans-HNSH, cis- HNSH and H₂SN. The other transitions range from very weak to moderately weak.

In the Franck-Condon region accessed from the corresponding electronic ground states, the upper electronic states of H₂NS, trans-HNSH, cis- HNSH and H₂SN are located above 4 eV. They correspond to electron excitations from the occupied inner valence MOs to the LUMO and LUMO+1 unoccupied MOs. The shapes of the resulting potentials are far from the standard Morse-like ones, since most of them are repulsive at large distances. Local minima can also exist at short distances and for non-planar configurations. Moreover, there is a high density of electronic states, which favors the mixing of their wavefunctions by vibronic and spin-orbit couplings. For instance, the quartets are crossing the doublets at intermediate R distances where spin-orbit conversion may occur. At these distances, the quartets are mostly repulsive and lead to fragmentations. For instance, the lowest quartet is crossing the electronic excited doublets at different nuclear geometries, where spin-orbit-induced predissociations may happen. At large internuclear separations, the reactions evolve along the repulsive part of this quartet state, which correlates adiabatically to lower dissociation limits. This is the case indeed for H₂NS, where a crossing between the H₂NS(2²A′) and the H₂NS(1⁴A″) states occurs at NH ~1.6 Å (Figure 3). For this doublet, we note also an avoided crossing with the 1²A′ state at NH ~1.5 Å (Figure 3), where vibronic coupling may occur. Interestingly this avoided crossing is responsible for the double well-shape of H₂NS(2²A′) and cis-HNSH(2²A′) along the SN distance. Also, the cuts of the PESs show that the second electronic excited state presents a flat potential along the torsion/out-of-plane coordinate resulting in a planar structure for H₂NS, whereas two minima can be seen for HNSH (cis and trans) and for H₂SN (C_2v and C_s). In all cuts displayed in Figures 3–6, the potential of such state has generally an odd shape, which may favor intramolecular isomerizations converting H₂NS into HNSH and H₂SN, and vice versa.

Structure and Spectroscopy of Lowest Electronic States of the H₂NS Species

For the investigation of the structure and vibrational spectroscopy of [H,H,N,S] in their electronic ground states, we used monoconfigurational approaches since their electronic wavefunctions, as discussed above, can be dominantly described by a unique configuration close to equilibrium. Thus, we characterized the stable forms of this molecular system in the ground state PES at the RCCSD(T)/aug-cc-pV(X+d)Z (X=D, T, Q) and RCCSD(T)-F12/cc-pVXZ-F12 (X=D, T, Q) levels of theory. Supplementary Tables 1–4 report the corresponding optimized geometrical parameters, harmonic frequencies and total energies for the minimal structures and Supplementary Table 5 those of the transition states. These tables show that the computed bond distances present slight variations by increasing the size of the basis set, where a systematic shortening of the bond distances is noticed because of the refined description of the electron correlation with the larger basis sets. Besides, variations in the bond angles are almost negligible upon changing the basis set.

We give in Table 2 the data after CBS extrapolation, together with their comparison to previous works. Computations confirm the existence of the already predicted H₂NS, trans-HNSH, cis-HNSH, and H₂SN forms and the transition states connecting them. The nature of these stationary points was carefully checked by harmonic frequency analysis, as given in Supplementary Tables 1–5. No other isomers were located. The relative energies confirm the previously established order: H₂NS followed by trans-HNSH, cis-HNSH and H₂SN. Table 2 shows also that all these compounds have relatively large dipole moments (1.1 < μ_e <3.7 debye). Therefore, these molecular species are susceptible to be observed by high-resolution rotational spectroscopy techniques in the laboratory and by radioastronomy. To-date, only H₂NS (X²B₁) was identified by Habara et al. (1998) using a source-modulation microwave spectrometer and a Fourier-transform microwave spectrometer. These authors deduced an accurate equilibrium molecular structure: R_NS = R₁ = 1.6398 Å, R_NH = R₃ = R₂ = 1.0055 Å and in-plane angle θ₁ = θ₂ = 120.5°. Our RCCSD(T)/CBS and RCCSD(T)-F12/CBS values compare well with this experimental geometry, with deviations of less than 0.6%. As expected, the less accurate previous calculations lead to larger deviations. For the other not yet detected isomers, we expect a similar accuracy. From the RCCSD(T)/CBS and RCCSD(T)-F12/CBS optimized geometries, we deduced the equilibrium rotational constants (A_e, B_e, C_e), to which we added the respective vibrational and centrifugal distortion corrections computed at the RMP2/aug-cc-pV(Q+d)Z level, within the VPT2 approach. Hence, at the RCCSD(T)/CBS level, the corrected (A₀, B₀, C₀) rotational constants of Table 2 are straightforwardly obtained from the following equation:

B₀(RCCSD(T)/CBS) = B_e(RCCSD(T)/CBS) + B₀(RMP2/aug-cc-pV(Q+d)Z) – B_e(RMP2/aug-cc-pV(Q+d)Z)

and similarly for A₀(RCCSD(T)/CBS), C₀(RCCSD(T)/CBS), A₀(RCCSD(T)-F12/CBS), B₀(RCCSD(T)-F12/CBS), and C₀(RCCSD(T)-F12/CBS). Table 2 shows that we have good agreement with the experimental data by Habara et al. (1998).

Table 2 lists the anharmonic frequencies of the fundamental modes (ν_i), which have been evaluated using the RCCSD(T)/aug-cc-pV(Q+d)Z//RMP2/aug-cc-pV(Q+d)Z composite scheme described in Hochlaf et al., 2015). Here, ν_i is equal to ω_i + Δν_i, where ω_i is the RCCSD(T)/aug-cc-pV(Q+d)Z harmonic frequency of the i-th mode, and Δν_i is the corresponding RMP2/aug-cc-pV(Q+d)Z anharmonic correction (Supplementary Table 6) as determined using second-order vibrational perturbation theory (VPT2) (Mills, 1972; Aliev and Watson, 1985), as implemented in GAUSSIAN. The validity of such procedure was widely discussed in ( Carbonnière et al., 2005; McKean et al., 2008; Temelso and Shields, 2011; Zen et al., 2012). Briefly, it allows to predict accurate and fully ab initio anharmonic vibrational frequencies at a low computational cost.

For H₂NS(X²B₁), we compute ν₁ = 3400.5, ν₂ = 1566.2, ν₃ = 938.9, ν₄ = 3503.5, ν₅ = 1015.4 and ν₆ = 412.9 (all values in cm⁻¹). These fundamentals are of a₁ symmetry for ν₁-ν₃, b₂ symmetry for ν₄-ν₅ and b₁ symmetry for ν₆ (out-of-plane mode). They compare favorably with the values by Habara et al. (1998) for the following modes: ν₁ (= 3430 cm⁻¹), ν₃ (= 918 cm⁻¹), ν₄ (= 3553 cm⁻¹) and ν₅ (= 989 cm⁻¹), while the comparison is less satisfactory for ν₂ (= 1328 cm⁻¹) and ν₆ (= 325 cm⁻¹). Also, we found that the ν₆ mode is strongly anharmonic, with inverted anharmonicity. This is related to the flatness of the PES along this normal coordinate. VPT2 treatment is not enough for deducing the ν₆ frequency. Instead, a full variational treatment in a 6D-PES is needed, which is beyond the scope of the present work.

For cis-HNSH and trans-HNSH, ν₁ through ν₅ are of a′ symmetry, and ν₆ is of a″ symmetry, while for H₂SN ν₁ through ν₄ are of a′ symmetry, and ν₅ and ν₆ are of a″ symmetry. Table 2 lists their fundamentals. Moreover, Table 2 lists the respective IR intensities, which show that all the possible modes are potentially detectable in the IR spectra.

Supplementary Figure 1 and Supplementary Table 5 give the optimized transition states (TS₁-TS₄) geometries, harmonic frequencies and total energies as calculated at the RCCSD(T)/aug-cc-pV(X+d)Z and RCCSD(T)-F12/cc-pVXZ (X=D,T,Q, CBS) levels. TS₁ connects cis-HNSH and trans-HNSH. The corresponding isomerization barrier is computed to be 0.88 eV at the RCCSD(T)/CBS level. TS₂ connects trans-HNSH and H₂NS. It is located at 1.84 eV above trans-HNSH. TS₃ and TS₄ are transition states along the H₂NS and H₂SN → H₂ + NS production channels, respectively. The energy barrier for H₂NS → TS₃ is relatively large (> 3 eV) and that of H₂SN → H₂ + NS is about 1 eV. The relative energies and the structures of these transition states are close and consistent with the data obtained by Pereira et al. (2010) using the CCSD/6-311++G(d,p) approach. The interconversion processes through these transition states may compete with the photophysical processes discussed in Section Photochemistry.

Table 3 lists the CASSCF/CBS optimized geometrical and spectroscopic parameters of the first electronic excited state for each isomer. The full set of data is given in Supplementary Table 7. Here, only three structures are found since cis- and trans-HNSH have the same excited state for non-planar HNSH arrangement. The comparison of the bond distances in H₂NS(X²B₁) and H₂NS(1²A′) reveals that the SN distance is lengthened upon electronic excitation, since the outermost electron is promoted from a bonding MO (SOMO-1) to an antibonding MO (SOMO), as illustrated in Figure 2. On the contrary, the NH bonds are only slightly affected. This is also the case for H₂SN, whereas no variations are seen for the HNSH species.

For the calculated electronic state of each isomer, we evaluated the radiative lifetime (τ_rad) using the formula

1/τ_rad = Σ_i 1/τ_i,rad, where τ_i,rad (= 6.07706 /(|R_e,i|² (ΔE_a,i)³) in μs) is the radiative lifetime populating the lower doublet states. Here, |R_e,i| (in debye) is the transition moment between the upper and the lower states as computed at the CASSCF/aug-cc-pV(5+d)Z level, and ΔE_a,i (in eV) is the corresponding CASSCF/MRCI/aug-cc-pV(5+d)Z adiabatic excitation energy.

Radiative lifetime computations were carried out at the ground state (τ_rad,g) and excited state (τ_rad,e) equilibrium geometries. The results are summarized in Table 3. This table shows that the calculated lifetimes are at least of several hundred ns, i.e., long enough to allow these electronic states to emit if populated.

Photochemistry

In order, to understand the thermodynamic aspects of the different dissociation processes taking place in [H,H,N,S] species or their formation from simple fragments, we give in Figure 7 the fragmentation energies to form diatomic + diatomic and atomic + triatomic entities, which are computed using:

or

where E(D), E(AB), E(CD), E(ABC) and E(ABCD) are the total energies of the D, AB, CD, ABC and ABCD entities. The ZPEs are their zero-point vibrational energies.

Figure 7

Figure 7 shows that the electronic ground state of each isomer is located below the lowest dissociation limit, i.e., HNS(X²A′)+H(²S). Therefore, all the calculated species are expected to be stable with respect to dissociation. All bimolecular reactions may lead to their formation. This is also the case for the first excited states of H₂NS and HNSH (denoted as H₂NS^* and HNSH^* in Figure 7). Such compounds are also expected to be long-lived. Table 3, which gives their radiative lifetimes, confirm such assumptions. Indeed, we compute lifetimes of several hundreds of ns for H₂NS^* and HNSH^*. On the other hand, the lifetime of the first electronic excited state of H₂SN (labeled H₂SN^*) may be reduced by spin-orbit-induced predissociation to populate the SN bond breaking channels HNS(X²A′)+H(²S), NH₂(X²B₁)+S(³P) or H₂S(X¹A₁)+N(⁴S). Figure 6 shows that the latter process is plausible since the H₂SN^*(1²A″) is crossed by the repulsive 1⁴A″ state which leads directly to these fragments. Figure 7 shows also that the second electronic excited states of all species are located above the lowest dissociation limit. Their lifetimes should be reduced by predissociation phenomena via the couplings with the quartet manifolds. Nevertheless, Table 3 shows that their radiative lifetimes are long enough to fluoresce whenever populated. For the upper electronic states (> 4 eV), they are located in energy above the lowest dissociation limit. They may undergo unimolecular decomposition processes forming H+HNS/HNS, H₂+NS or NH+SH. A multi-step mechanism along the PESs of their doublet states is expected, where the wavepackets on these PESs evolve until reaching the regions of conical intersections or avoided crossings. The electronic de-excitation of the H₂NS, trans-HNSH, cis- HNSH and H₂SN molecules may occur via emission of light or more likely by radiationless internal conversion processes. The lowest states are thus populated. Alternatively, a conversion to the quartet manifold takes place at the doublet-quartet crossings, followed by dissociation.

In their pioneering work on “active nitrogen,” Strutt, Fowler and Vaidya (Strutt and Fowler, 1912; Strutt, 1913; Fowler and Vaidya, 1931) and later on Westbury and Winkler (Westbury and Winkler, 1960) showed that the reactions of “active nitrogen” with hydrogen sulfide occurring in a flame are associated with a deep blue glow that appeared to fill the entire upper half of the reaction vessel. This light emission was interpreted as emanating from the N.H₂S tetratomic intermediate. Computations show that the bimolecular reactive collisions, in the presence of a third party to absorb the excess energy, should lead to the formation of [H,H,N,S] tetratomics in their electronic ground states as established by Pereira et al. (2010) for the ¹NH + ²SH → cis-²HNSH → trans-²HNSH → ²H₂NS → ²NS + H₂ reactive channels. Also the [H,H,N,S] tetratomics can be formed in their first electronic excited states, which are long-lived enough to emit light in the IR domain (hν <1.2 eV, Table 3). Nevertheless, such light was not observed in the experiment most likely because it was not possible to measure it at the beginning of the XX^th century. Table 3 shows also that the second electronic excited state of [H,H,N,S] species can be a suitable candidate to explain the deep blue glow. Indeed the H₂SN(X²A′) ← H₂SN(2²A′) transition has an energy of ~3.8 eV. As mentioned above the PES of this electronic state is flat, so intramolecular isomerization to form the other isomers is possible. At least, we may obtain the cis-HNSH (2²A′) species, which is long lived enough to emit light with hν ~ 2.8 eV i.e., close to the ~2.7 eV experimentally observed deep blue glow. Thus, our work establishes that the less stable forms HSNH and H₂SN are responsible for the flame color. Figure 6 shows that the upper state can be populated after collisions between N(²D) or N(²P) with H₂S, which means that the N atom is not in its electronic ground state (⁴S). Accordingly, the so-called “active nitrogen” corresponds to electronically excited nitrogen atoms that might be formed in the flame. The reactions evolve on the PESs of the electronic excited states of doublet spin multiplicity.

Conclusions

Using advanced ab initio methodologies, we mapped the PESs of the lowest electronic states of the [H,H,N,S] molecular species. Thence, we characterized the stable forms in the ground and two lowest electronic states. For the ground states, we confirm the existence of the four already known isomers. Since the corresponding experimental data are estimates, we suggest that the IR spectra of these compounds should be measured again for a definitive conclusion. On the excited potentials, the H₂NS, HNSH and H₂SN forms were identified for the first time. Also, we used these potentials to examine the unimolecular decomposition processes undergone by the H₂NS, HNSH and H₂SN species and the bimolecular reactions processes between H and HNS/HSN, NH + SH and H₂S + N. The findings relative to the latter process were used to confirm the observation of H₂SN fluorescent as an intermediate during the reactions between “active nitrogen” and H₂S, known since the beginning of the XX^th century. Our data should help in identifying [H,H,N,S] molecular compounds in the laboratory as well as in astrophysical and environmental media.

References

• 1

  AdlerT. B.KniziaG.WernerH.-J. (2007). A simple and efficient CCSD(T)-F12 approximation. J. Chem. Phys. 127:221106. 10.1063/1.2817618

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AjiliY.Ben AbdallahD.Mogren Al-MogrenM.FranciscoJ. S.HochlafM. (2016). Rotational (de-)excitation of HNS by He: Three-dimensional potential energy surface and collision rate coefficients. Mon. Not. R. Astron. Soc.458, 1581–1589. 10.1093/mnras/stw371

  • CrossRef
  • Google Scholar

• 3

  AlievM. R.WatsonJ. K. G. (1985). Molecular Spectroscopy: Modern Research, Vol. III, edited by RaoK. N.. New York, NY: Academic Press.

  • Google Scholar

• 4

  Ben YaghlaneS.JaidaneN.-E.CottonC. E.FranciscoJ. S.Mogren Al MogrenM.LinguerriR.HochlafM. (2014). Theoretical spectroscopic investigations of HNS^q and HSN^q (q = 0,+1,-1) in the gas phase. J. Chem. Phys. 140:244309. 10.1063/1.4883915

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BenYaghlaneS.CottonC. E.FranciscoJ. S.LinguerriR.HochlafM. (2013). Ab initio structural and spectroscopic study of HPS^x and HSP^x (x = 0,+1,−1) in the gas phase. J. Chem. Phys. 139:174313. 10.1063/1.4827520

  • CrossRef
  • Google Scholar

• 6

  CarbonnièreP.LuccaT.PouchanC.RegaN.BaroneV. (2005). Vibrational computations beyond the harmonic approximation: performances of the B3LYP density functional for semirigid molecules. J. Comput. Chem. 26, 384–388. 10.1002/jcc.20170

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  DeeganM. J. O.KnowlesP. J. (1994). Perturbative corrections to account for triple excitations in closed and open shell coupled cluster theories. Chem. Phys. Lett. 227, 321–326. 10.1016/0009-2614(94)00815-9

  • CrossRef
  • Google Scholar

• 8

  DunningT. H.PetersonK. A.WilsonA. K. (2001). Gaussian basis sets for use in correlated molecular calculations. X. The atoms aluminum through argon revisited. J. Chem. Phys. 114, 9244–9253. 10.1063/1.1367373

  • CrossRef
  • Google Scholar

• 9

  EspinT. E. (1890). Über helle linien im spectrum von R andromedae. Astronomische Nachrichten123:31. 10.1002/asna.18901230106

  • CrossRef
  • Google Scholar

• 10

  FowlerA.VaidyaW. H. (1931). The spectrum of the flame of carbon disulphide. Proc. Roy. Soc. A132, 310–330. 10.1098/rspa.1931.0103

  • CrossRef
  • Google Scholar

• 11

  GottliebC. A.BallJ. A.GottliebE. W.LadaC. J.PenfieldH. (1975). Detection of interstellar nitrogen sulfide. Astrophys. J. 200, L147–L149. 10.1086/181918

  • CrossRef
  • Google Scholar

• 12

  HabaraH. S.Yamamoto SaitoS. (1998). Microwave spectrum and molecular structure of the H₂NS radical. J. Chem. Phys. 109, 2700–2707. 10.1063/1.476876

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  HampelC.PetersonK.WernerH.-J. (1992). A comparison of the efficiency and accuracy of the quadratic configuration interaction (QCISD), coupled cluster (CCSD), and Brueckner coupled cluster (BCCD) methods. Chem. Phys. Lett. 190, 1–12. 10.1016/0009-2614(92)86093-W

  • CrossRef
  • Google Scholar

• 14

  HättigC. (2005). Optimization of auxiliary basis sets for RI-MP2 and RI-CC2 calculations: Core–valence and quintuple-ζ basis sets for H to Ar and QZVPP basis sets for Li to Kr. Phys. Chem. Chem. Phys. 7, 59–66. 10.1039/B415208E

  • CrossRef
  • Google Scholar

• 15

  HelgakerT.KlopperW.KochH.NogaJ. (1997). Basis-set convergence of correlated calculations on water. J. Chem. Phys. 106, 9639–9646. 10.1063/1.473863

  • CrossRef
  • Google Scholar

• 16

  HochlafM.PuzzariniC.SenentM. L. (2015). Towards the computations of accurate spectroscopic parameters and vibrational spectra for organic compounds. Mol. Phys. 113, 1661–1673. 10.1080/00268976.2014.1003986

  • CrossRef
  • Google Scholar

• 17

  IrvineW. M.SenayM.LovellA. J.MatthewsH. E.McGonagleD.MeierR. (2000). Detection of nitrogen sulfide in Comet Hale-Bopp. Icarus143, 412–414. 10.1006/icar.1999.6281

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  KendallR. A.DunningT. H.Jr.HarrisonR. J. (1992). Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 96, 6796–6806. 10.1063/1.462569

  • CrossRef
  • Google Scholar

• 19

  KlopperW. (2001). Highly accurate coupled-cluster singlet and triplet pair energies from explicitly correlated calculations in comparison with extrapolation techniques. Mol. Phys. 99, 481–507. 10.1080/00268970010017315

  • CrossRef
  • Google Scholar

• 20

  KniziaG.AdlerT. B.WernerH.-J. (2009). Simplified CCSD(T)-F12 methods: Theory and benchmarks. J. Chem. Phys. 130:054104. 10.1063/1.3054300

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  KniziaG.WernerH.-J. (2008). Explicitly correlated RMP2 for high-spin open-shell reference states. J. Chem. Phys. 128:154103. 10.1063/1.2889388

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  KnowlesP. J.HampelC.WernerH.-J. (2000). Coupled cluster theory for high spin, open shell reference wave functions. J. Chem. Phys. 99, 5219–5227. 10.1063/1.465990

  • CrossRef
  • Google Scholar

• 23

  KnowlesP. J.WernerH.-J. (1985). An efficient second-order MC SCF method for long configuration expansions. Chem. Phys. Lett. 115, 259–267. 10.1016/0009-2614(85)80025-7

  • CrossRef
  • Google Scholar

• 24

  McGonagleD.IrvineW. M. (1997). Nitrogen sulfide in giant molecular clouds. Astrophys. J. 477, 711–721. 10.1086/303749

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  McGonagleD.IrvineW. M.OhishiM. (1994). Nitrogen sulfide in quiescent dark clouds. Astrophys. J. 422, 621–625. 10.1086/173755

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  McKeanD. C.CraigN. C.LawM. M. (2008). Vibrational Anharmonicity and Harmonic Force Fields for Dichloromethane from Quantum-Chemical Calculations. J. Phys. Chem. A112, 10006–10016. 10.1021/jp803881c

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  MeyerD. M.RothK. C. (1991). Discovery of interstellar NH. Astrophys. J. Letters376:L49–L52. 10.1086/186100

  • CrossRef
  • Google Scholar

• 28

  MillarT. J.HerbstE.CharnleyS. B. (1991). The formation of oxygen-containing organic molecules in the Orion compact ridge. Astrophys. J. 369, 147–156. 10.1086/169745

  • CrossRef
  • Google Scholar

• 29

  MillsI. M. (1972). Molecular Spectroscopy: Modern Research, edited by RaoK.N.MathewsC.W.. New York, NY: Academic Press.

  • Google Scholar

• 30

  MiuraY.AsadaH.KlnoshitaM.OhtaK. (1983). Electron spin resonance study of N-Alkyl-N-(a1kylthio)aminyl Radicals. J. Phys. Chem. 87, 3450–4555. 10.1021/j100241a020

  • CrossRef
  • Google Scholar

• 31

  NeugebauerG.MartzD. E.LeightonR. B. (1965). Observations of extremely cool stars. Astrophys. J. 142, 399–401. 10.1086/148300

  • CrossRef
  • Google Scholar

• 32

  NguyenM. T.VanquickenborneL. G.FlammangR. (1994). The thionitroxyl free radical (H₂NS) and its ionic counterparts (H₂NS⁺ and H₂NS⁻): a theoretical and experimental study. J. Chem. Phys. 101, 4885–4892. 10.1063/1.467410

  • CrossRef
  • Google Scholar

• 33

  PereiraP. S. S.MacedoL. G. M.PimentelA. S. (2010). New insight into the formation of nitrogen sulfide: a quantum chemical study. J. Phys. Chem. A114, 509–515. 10.1021/jp907384d

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  PetersonK. A.AdlerT. B.WernerH.-J. (2008). Systematically convergent basis sets for explicitly correlated wavefunctions: the atoms H, He, B–Ne, and Al–Ar. J. Chem. Phys. 128:084102. 10.1063/1.2831537

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  PhuongN. T.ChapillonE.MajumdarL.DutreyA.GuilloteauS.PiétuV.et al. (2018). First detection of H₂S in a protoplanetary disk. The dense GG Tauri A ring. A&A616:6. 10.1051/0004-6361/201833766

  • CrossRef
  • Google Scholar

• 36

  ShamasundarK. R.KniziaG.WernerH.-J. (2011). A new internally contracted multi-reference configuration interaction method. J. Chem. Phys. 135:054101. 10.1063/1.3609809

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  StruttR. J. (1913). An active modification of nitrogen, produced by the electric discharge.—V. Proc. Roy. Soc. A88, 539–549. 10.1098/rspa.1913.0049

  • CrossRef
  • Google Scholar

• 38

  StruttR. J.FowlerA. (1912). Spectroscopic investigations in connection with the active modification of nitrogen. II—Spectra of elements and compounds excited by the nitrogen. Proc. Roy. Soc. (London) A86, 105–117. 10.1098/rspa.1912.0005

  • CrossRef
  • Google Scholar

• 39

  TemelsoB.ShieldsG. C. (2011). The role of anharmonicity in hydrogen-bonded systems: the case of water clusters. J. Chem. Theory Comput. 7, 2804–2817. 10.1021/ct2003308

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  TrabelsiT.Ben YaghlaneS.Mogren Al MogrenM.FranciscoJ. S.HochlafM. (2016). HNS⁺ and HSN⁺ cations: electronic states, spin-rovibronic spectroscopy with planetary and biological implications. J. Chem. Phys. 145:084307. 10.1063/1.4961313

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  TrabelsiT.LinguerriR.Ben YaghlaneS.JaidaneN.-E.Mogren Al-MogrenM.FranciscoJ. S.HochlafM. (2015). On the role of HNS and HSN as light-sensitive NO-donors for delivery in biological media. J. Chem. Phys. 143:134301. 10.1063/1.4932084

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  TrabelsiT.Mogren Al-MogrenM.HochlafM.FranciscoJ. S. (2018). Mechanistic study of the photoexcitation, photoconversion, and photodissociation of CS₂. J. Chem. Phys. 149:064304. 10.1063/1.5040141

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  van DishoeckE. F.JansenD. J.SchilkeP.PhillipsT. G. (1993). Detection of the Interstellar NH₂ Radical. Astrophys. J. Letters, 416: L83–L86. 10.1086/187076

  • CrossRef
  • Google Scholar

• 44

  WagenblastR.WilliamsD. A.MillarT. J.NejadL. A. M. (1993). On the origin of NH in diffuse interstellar clouds. Mon. Not. R. Astron. Soc.260, 420–424. 10.1093/mnras/260.2.420

  • CrossRef
  • Google Scholar

• 45

  WeigendF. (2002). A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency. Phys. Chem. Chem. Phys. 4, 4285–4291. 10.1039/b204199p

  • CrossRef
  • Google Scholar

• 46

  WernerH.-J.KnowlesP. J. (1985). A second order multiconfiguration SCF procedure with optimum convergence. J. Chem. Phys. 82, 5053–5063. 10.1063/1.448627

  • CrossRef
  • Google Scholar

• 47

  WernerH.-J.KnowlesP. J. (1988). An efficient internally contracted multiconfiguration–reference configuration interaction method. J. Chem. Phys. 89, 5803–5814. 10.1063/1.455556

  • CrossRef
  • Google Scholar

• 48

  WernerH. J.AdlerT. B.ManbyF. R. (2007). General orbital invariant MP2-F12 theory. J. Chem. Phys. 126:164102. 10.1063/1.2712434

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  WernerH. J.KnowlesJ.et al. (2015). MOLPRO version 2015 is a package of ab initio programs. University of Cardiff Chemistry Consultants (UC3): Cardiff, Wales, UK. Available online at: http://www.molpro.net

  • Google Scholar

• 50

  WestburyR. A.WinklerC. A. (1960). The reactions of active nitrogen with hydrogen sulfide and carbon disulfide. Can. J. Chem. 38, 334–342. 10.1139/v60-049

  • CrossRef
  • Google Scholar

• 51

  WöhlH. (1971). On Molecules in Sunspots. Solar Phys. 16, 362–272. 10.1007/BF00162477

  • CrossRef
  • Google Scholar

• 52

  WoonD. E.DunningT. H.Jr. (1993). Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 98, 1358–1371. 10.1063/1.464303

  • CrossRef
  • Google Scholar

• 53

  XinJ.FanH.IonescuI.AnnesleyC.ReidS. A. (2003). Fluorescence excitation spectroscopy of the A²A₁ - X²B₁ system of jet-cooled NH₂ in the region 2900–4300Å. J. Mol. Spectrosc. 219, 37–44. 10.1016/S0022-2852(03)00008-0

  • CrossRef
  • Google Scholar

• 54

  YamamuraS. T.KawaguchiK.RidgwayS. T. (2000). Identification of SH v=1 Ro-vibrational Lines in R Andromedae. Astrophys. J. 528, L33–L36. 10.1086/312420

  • CrossRef
  • Google Scholar

• 55

  ZenA.ZhelyazovD.GuidoniL. (2012). Optimized structure and vibrational properties by error affected potential energy surfaces. J. Chem. Theory Comput. 8, 4204–4215. 10.1021/ct300576n

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  ZiurysL. M. (2006). The chemistry in circumstellar envelopes of evolved stars: Following the origin of the elements to the origin of life. Proc. Natl. Acad. Sc. U.S.A.103, 12274–12279. 10.1073/pnas.0602277103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar