Victoria H. J.
Clark
* and
Sergei N.
Yurchenko
*
Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BT London, UK. E-mail: v.clark.17@ucl.ac.uk; s.yurchenko@ucl.ac.uk
First published on 29th April 2021
This paper sets out a robust methodology for modelling spectra of polyatomic molecules produced in reactive or dissociative environments, with vibrational populations outside local thermal equilibrium (LTE). The methodology is based on accurate, extensive ro-vibrational line lists containing transitions with high vibrational excitations and relies on the detailed ro-vibrational assignments. The developed methodology is applied to model non-LTE IR and visible spectra of silylene (SiH2) produced in a decomposition of disilane (Si2H6), a reaction of technological importance. Two approaches for non-LTE vibrational populations of the product SiH2 are introduced: a simplistic 1D approach based on the Harmonic approximation and a full 3D model incorporating accurate vibrational wavefunctions of SiH2 computed variationally with the TROVE (Theoretical ROVibrational Energy) program. We show how their non-LTE spectral signatures can be used to trace different reaction channels of molecular dissociations.
The modern day study of non-LTE spectroscopy can be traced as far back as the 1930s to the original papers of Milne,14 and the many key papers from the decades following.15–24 The non-LTE spectroscopic effects play important role in high-resolution applications and there exist a number of accurate non-LTE spectroscopic and radiative transfer codes, see van der Tak et al.,25 Funke et al.,26 Pannier and Laux27 and references therein. As such, non-LTE spectra are often vital for the modelling of astrophysical problems, including planetary atmospheric properties,28 stellar atmospheres of solar system and exoplanets29–36 and the ISM.25,37,38
The non-local thermodynamic effects of the spectra of molecules has been of interest to chemists and astronomers alike for many years. A notable example of this is the 2011 work by Ferus et al.2 who studied the isomers of HCN within acetonitrile, formamide, and BrCN discharge. Using the features for both the HCN and HNC molecules from these spectra, Ferus et al.2 were able to calculate the ratios of molecules within the reactions and also the reaction path taken by the HNC molecule during the isomerization. In this work we explore this idea to study non-LTE spectral signatures of silylene (SiH2) produced from disilane (Si2H6).
Reaction properties of silylene, silane and disilane such as the rate constants for the formation, destruction and chemical pathways are important for plasma physics aspects such as silicon deposition.39–41 The ease of hydrogen transfer and high barriers in the saturated silicon system, leading to the ready formation of three-centre interactions and consequently the isomerisation reactions of Si2H6, are just as important to study as the elementary reactions.42 Silane containing reactions are also of importance for astrophysics, with the presence of SiH4 in IRC + 10216 discussed by Goldhaber and Betz,43 Kaiser and Osamura.44
The complexity of the silane containing systems has been discussed elsewhere,42,45 with quantitative calculations proving particularly difficult. There were a number of ab initio studies of the structural properties of Si2H6e.g.,46–53 as well as of the formation and reactions involving this molecule.54–58 Agrawal et al.53 produced a global ab initio potential energy surface of disilane and used it to investigate the dissociation dynamics with classical trajectories. Márquez et al.50 reported a force field for Si2H6. The main isomer of Si2H6 has a staggered, ethane-like, structure (see Fig. 1) with a low barrier (1.2 kcal mol−1) to the eclipsed conformation.59 Si2H6 has been shown to have a local minimum as an inverted stable structure with one of the Si–H3 ‘umbrellas’ pointing to the center as well as a transition state with a similar inverted configuration.46,60 These structures are nominally asymmetric (Cs symmetry) but essentially acquiring the C3v symmetry.
![]() | ||
Fig. 1 Schematic reaction path for the silylene. Molecular SiH2 is the starting material, fragment SiH2 is the product. |
Thermal decomposition of Si2H6 has been extensively studied, both theoretically (mostly using RRKM, Rice–Ramsperger–Kassel–Marcus) and experimentally, with the reaction Si2H6 → SiH2 + SiH4 as the most common53,56,57,61 and important decomposition process of the excited disilane,53,57,62–66 and where Arrhenius parameters and rate constants have been reported (e.g. Bowrey and Purnell,61 Martin et al.,63 Mick et al.,65 Roenigk et al.67). However, it is also possible for the disilane molecule to dissociate homolytically, as was originally thought to be the main pathway owing to disilanes similarities with ethane, and form ,53,56,68 or to undergo dehydrogenation to H2Si–SiH2, H3Si–SiH, or H3Si–H2.53,54,57 Disilane can even undergo double dehydrogenation to form Si2H2, however to our knowledge this has only been reported as the main product when undergoing photolysis at 193 nm.69 Yoshida et al.57 notes that the transition state for Si2H6 → SiH2 + SiH4 is 8.48 kcal mol−1 lower than the transition state for Si2H6 → H3SiSiH + H2, at 43.38 kcal mol−1 compared to 51.86 kcal mol−1.
The spectroscopy of SiH2 has been used to monitor the SiH2 + SiH4 → Si2H6 and Si2H6 → SiH2 + SiH4 reactions (see Fig. 1) and measure the corresponding rate constants and Arrhenius parameters by spectroscopically tracking electronic (Ö) transitions of SiH2.39,46,65,70 In these studies, the reconstructions of the amount of SiH2 relied on the assumption of the Boltzmann thermal distribution when estimating the population of the lower state. No account of the possible non-LTE population of SiH2 molecules after dissociation was made, which could potentially hamper the count of the SiH2 molecules and affect the reaction rates estimated. A similar experimental technique was used in Hertl and Jolly71 to monitor SiH2 in SiH4 plasma.
It is the second, dissociation, part of the reaction shown in Fig. 1 (Si2H6 → SiH2 + SiH4) we study in this work. More specifically, we show that (i) the (vibrational) populations of the molecules produced in reactions can be very different from the Boltzmann distribution and is important to take into account when interpreting spectroscopic measurements. That (ii) spectral shapes of the dissociated SiH2 can bear strong non-LTE character, very different from the LTE spectrum of an LTE SiH2 sample making it possible to distinguish between different reaction stages and even between different dissociation channels the silylene molecules it is produced from. In this work non-LTE spectra of SiH2 under conditions similar to dissociation processes expected in these experiments are modelled.
Recently we have computed an accurate ro-vibrational line list for SiH2, named CATS.72 It covers a large range of rotational and vibrational excitation, capable of modelling very hot spectra of this molecule (up to T = 2000 K) as part of the ExoMol database.73 The CATS line list was produced using the program TROVE,74,75 which solves the nuclear motion Schrödinger equation variationally. The ro-vibrational energies and corresponding wavefunctions were computed using an accurate, empirically refined potential energy surface (PES) of silane and a high-level ab initio dipole moment surface (DMS). The ro-vibrational probabilities (in the form of Einstein A coefficients) were computed using a high level ab initio dipole moment surface.
The study by Clark et al.72 forms the basis for the present work, where we utilize the CATS line list, wavefunctions, purpose-built numerical basis set, and the CATS computational TROVE setup to model non-LTE spectroscopic properties of SiH2 produced from dissociation of Si2H6 through different reaction channels. Using a simplified 1D Harmonic oscillator wavefunctions (see Pastorek et al.6) and more sophisticated 3D vibrational CATS wavefunctions from accurate variational calculations, the non-LTE ro-vibrational populations of SiH2 are generated and used to produce non-LTE spectroscopic spectra of different dissociation channels of disilane. To this end we investigate reaction topology connecting the global minimum of Si2H6 with the closest saddle points and local minima as well as the corresponding structural properties using a high level ab initio theory cc-pVTZ-F12/CCSD(T)-F12b76,77 employing the program MOLPRO2015.78
Theoretically, the non-LTE properties of dissociating molecules were studied by Band and Freed.79 In the present work we use general approach of Berry80 and Band and Freed,79 which assume no significant structural changes between the reactant and product nuclear configuration, along with the slow vibrational relaxation of the product13 to investigate non-LTE spectroscopic signatures of SiH2 produced from dissociation of Si2H6. In this paper we specifically consider situations where the vibrational relaxations are not achieved during the time of the experiment, so that the molecules still hold the memory of the structure during the reaction or dissociation. The rotation relaxation time however is much shorter and the rotational degrees of freedom can be usually assumed to satisfy the Bolzmann equilibrium.13
We also investigate possible non-LTE impact on the electronic Ã(0,2,0)–(0,0,0) spectrum of SiH2. This is a favorite spectroscopic system for the detection of SiH2 due to the large Franck–Condon factor and the availability of suitable laser.39,70,81–91
The non-LTE absorption spectra of SiH2 are simulated using the (non-LTE) ExoCross program,92 where a new feature of non-Boltzmann populations was added. ExoCross has been previously used to model spectra of molecules in environments that can be characterized using two temperatures, vibrational and rotational.93–96
The paper is structured as follows. In Section 2 we describe the calculations of potential energy surfaces for the disilane and silylene structures. The theory used in this paper is described in Section 3. In Section 4.1 we calculate the 1D harmonic wavefunction population and use them to produce non-LTE spectra of SiH2 corresponding to different dissociation routes. In Section 4.2 we calculate the populations and subsequent non-LTE ro-vibrational spectra of SiH2 using the full 3D wavefunctions and describe the new TROVE methodology. A non-LTE electronic Ã(0,2,0)–(0,0,0) spectrum of SiH2 is presented in Section 4.4. Conclusions are offered in Section 5.
We shall refer to different disilane isomers as dGM, dLM and dTS to distinguish them from the SiH2 fragments GM, LM and TS as discussed below.
The global minimum of SiH2 (dGM) has a symmetrical, staggered D3d structure. The closest local minimum (dLM) has an inverted, C3v structure. The lowest transition state, which will be referred to as dTS (TS1 in Becerra et al.46 and TS2 in Tonokura et al.103), has also inverted structure, just a few kJ mol−1 above dLM.46 These structures together with the corresponding optimized parameters for three geometries most relevant to our work are collected in Table 1. Our structural parameters compare well with that from the literature. The structure of the dGM has also been determined spectroscopically, with the equilibrium bond lengths rSi–H = 1.492 Å, rSi–Si = 2.331 Å and bond angles β∠HSiSi = 110.3° and α∠HSiSi = 108.6°.104 The structure of the deuterated isotopologue Si2H5D was reported as rSi–H = 1.4874(17) Å, rSi–Si = 2.3317(15) Å, β∠HSiSi = 110.66(16)°.105
Isomer | LM | TS | GM | Becerra et al.46 | ||
---|---|---|---|---|---|---|
Molecule | LM1 | TS1 | GM | |||
Si–Si (Å) | 2.467 | 2.489 | 2.341 | 2.470 | 2.485 | 2.346 |
SiR–HR (Å) | 1.532 | 1.499 | 1.482 | 1.525 | 1.493 | 1.479 |
SiL–HL (Å) | 1.480 | 1.480 | 1.482 | 1.478 | 1.478 | 1.479 |
∠HRSiSi (°) | 55.4 | 64.9 | 110.2 | — | — | 110.24 |
∠HLSiSi (°) | 110.7 | 110.0 | 110.2 | 110.7 | 110.2 | 110.24 |
∠HRSiHR (°) | 90.9 | 100.6 | 108.6 | 91.4 | 102.6 | — |
∠HLSiHL (°) | 108.3 | 109.7 | 108.6 | 108.2 | 108.8 | — |
ω b (cm−1) | 940 | 940 | 940 | |||
ω s(L) (cm−1) | 2241.69 | 2245.36 | 2240.10 | |||
ω s(R) (cm−1) | 2039.12 | 2170.66 | 2240.10 |
The reaction path connecting the disilane isomers dGM, dLM and dTS is shown in Fig. 2. A zoom of the dLM side is shown as inset. The energies of the global and local minima are 178 kJ mol−1 and 4.1 kJ mol−1 below the transition state, respectively. The energy and geometry information for dLM, dGM, dTS are collected in Table 2. The results calculated compare well with both the results of Becerra et al.46 and Sakai and Nakamura,60 albeit both vary for the dGM structure by 36 kJ mol−1.
Energy, Eh | Rel. energy,d kJ mol−1 | |
---|---|---|
a TS2 in Sakai and Nakamura60 is the same as TS1 in Becerra et al.46 b Complex 2 in Sakai and Nakamura60 is the same as LM1 in Becerra et al.46 c Relative to SiH2 + SiH4, kJ mol−1. d Relative to dTS. | ||
dGM | ||
Calculated | −581.728134 | −178.00 |
Becerra | −236.0 | −214.50 |
Saki | −581.53962 | −214.68715 |
dTS | ||
Calculated | −581.660338 | 0.00 |
Becerra (TS1)a | −48.5 | 0.00 |
Saki (TS2)a | −581.45785 | 0.00 |
dLM | ||
Calculated | −581.661906 | −4.12 |
Becerra (LM1)b | −51.3 | −3.00 |
Saki (compl. 2)b | −581.45837 | −1.3652601 |
The dLM isomer of Si2H6 has a shallow potential with a very low barrier to dTS of 344 cm−1. It can be also recovered using lower levels of theory, for example, using MP2/6-311G(d,p)46 and even with the UFF force fields implemented in Avogadro 1.2.0106via the steepest descent method and 4 steps per update.
dGM, D3d | dTS, C3v | dLM, C3v | |||||||
---|---|---|---|---|---|---|---|---|---|
Calculated | Literature46 | Calculated | Literature46,67 | Calculated | |||||
Mode | Freq. | Mode | Freq. | Mode | Freq. | Mode | Freq. | Mode | Freq. |
a Confirmed by the experimental gas phase Raman spectrum of Durig and Church.59 b Although there are three lines of 925 cm−1, there are no triply degenerate states in the C3v symmetry. Roenigk et al.67 assign one 925 cm−1 to ν15, which in the dGM symmetry is an E state degenerate with ν16. The dTS ν16 can be found at 960 cm−1. The other two 925 cm−1 line were assigned to ν9 and ν10 (E). | |||||||||
E u | 2247.78 | E u | 2179 | E | 2252.06 | A | 2181 | E | 2244.30 |
A 1g | 2241.03 | A 1g | 2163a | A | 2238.66 | E | 2169 | A | 2239.08 |
E g | 2239.28 | E u | 2155 | E | 2176.47 | A | 2105 | A | 2050.89 |
A 1u | 2232.30 | A 1u | 2154 | A | 2161.85 | A | 2087 | E | 2027.35 |
E u | 962.35 | E g | 941a | E | 995.37 | A | 1585 | E | 1089.15 |
E g | 948.37 | E u | 940 | E | 966.95 | A | 960 | E | 978.33 |
A 1g | 931.91 | A 1g | 920a | A | 902.23 | A | 949 | A | 956.33 |
A 1u | 859.02 | A 1u | 844 | E | 570.16 | A | 927 | A | 667.30 |
E g | 636.87 | E g | 628a | A | 353.88 | A & Eb | 925 | E | 618.92 |
A 1g | 438.76 | A 1g | 432a | E | 319.45 | E | 145 | A | 432.13 |
E u | 372.85 | E u | 379 | A | 316.43 | A | 128 | A | 389.31 |
A 1u | 136.56 | A 1u | 128 | A | −551.59 | E | 93 | E | 277.91 |
In our description of the non-LTE vibrational population of the dissociated molecule we follow the polyatomic Franck–Condon type approximation by Band and Freed79 and Berry80 based on the structural differences between reactant and product assuming no significant change in nuclear configuration of the molecule. In order to connect the product (gas phase SiH2) to an initial structure of Si2H6 through the dissociation process, we assume that the dissociation happens instantaneously, i.e. the initial configuration of the product SiH2 corresponds to the structural parameters (bond lengths Si–H and inter-bond angles ∠HSiH) of SiH2 as a fragment of Si2H6, for which the parameters collected in Table 1 are used. For example, for the dissociation from the dGM structure, the initial configuration of the gas phase SiH2 is assumed to be rSi–H = 1.482 Å, α∠HSiH = 108.6°. Naturally, this is a very deformed geometry comparing to the equilibrium structure of the gas phase SiH2, re = 1.5137 ± 0.003 Å and αe = 92.04 ± 0.05°.107 After being dissociated, in relaxing to be a free molecule, the fragment SiH2 has added vibrational energy and is hence in a situation when it's vibrational populations do not match the LTE distribution for corresponding temperature of the surroundings, at least for the vibrational degrees of freedom.
Experience shows that the rotational equillibration time is usually very short and we can thus safely assume the LTE conditions for the rotational degrees of freedom with the rotational temperature the same as the temperature of the surroundings (see also e.g. Dudás et al.13). The vibrational population however is not in the thermal equilibrium and therefore no sensible vibrational temperature could be associated with the corresponding population.
The non-LTE population NJ,k,ν(T) of a ro-vibrational state |J, k, ν〉 is then given by:6
![]() | (1) |
ẼJ,k,ν = Ẽvibν + ẼrotJ,k, | (2) |
![]() | (3) |
SiH2 fragment name | Formation & structure |
---|---|
GM | From dGM, two equivalent SiH3 |
TS-R | From dTS, SiH3 pointing inside |
TS-L | From dTS, SiH3 pointing outside |
LM-R | From dLM, SiH3 pointing inside |
LM-L | From dLM, SiH3 pointing outside |
Let us assume that the Si2H6 is LTE and hence is in its ground vibrational state at the moment of dissociation, while the fragment SiH can end up in any vibrationally excited state ≡ (ν1,ν3) with some transition probability giving rise to the vibrational population Nvibν. On top of that we also assume a full separation of the stretching Si–H and bending ∠HSiH modes inside Si2H6 in its ground vibrational state. Possible consequences of deviation from these approximations are discussed below.
Under the assumptions made we define the vibrational population of the SiH2 fragment as a Franck–Condon factor for a vertical transition from the ground vibrational state of disilane Si2H6 to SiH2 + SiH4 with the gas phase (g.ph.) SiH2 transferred to some vibrational state |(g.ph.)〉 = |ν1,ν2,ν3〉. In the approximation of the full separation of the fragment SiH2 from the rest of Si2H6, the population of SiH2 can be represented as an overlap between the ground state wavefunction |ν′′ = 0(fragment)〉 of a fragment SiH2 and that of the corresponding vibrational state |ν′ = (g.ph.)〉 of the gas phase SiH2 as given by
N = |〈0(fragment)|ν(g.ph.)〉|2. | (4) |
The calculated temperature dependent populations NJ,k,ν(T) in eqn (1) can be then combined with a molecular line list for SiH2 to simulate absorption or emission spectra of this molecule under the non-LTE conditions as defined by NJ,k,ν(T) in eqn (1). Here we use the ExoMol line list CATS by Clark et al.72 as provided by ExoMol (http://www.exomol.com). Technically this is done by incorporating the non-LTE vibrational densities N into the ExoMol States file as described in Section 4 (the ExoMol file formats are discussed extensively elsewhere73). A non-LTE spectrum of SiH2 for given T and P is then calculated using CATS' Einstein A coefficients with the ExoCross program,92 where a new non-LTE option has been implemented as part of this work. The rotational populations are assumed to be in LTE according with eqn (1).
![]() | (5) |
![]() | (6) |
![]() | (7) |
![]() | (8) |
![]() | (9) |
The constants as and ab correspond to inverse masses of the vibrational part of a free three-atomic molecule expressed in terms of the internal coordinates r1, r2, α (see, e.g. Sutcliffe and Tennyson,108 Yurchenko et al.109).
A 1D population for the corresponding mode of the gas phase SiH2 molecule is given by eqn (4) with |ν〉 = Ψv(x). The different disilane fragments have different structural parameters re, αe, ωs and ωb, see Tables 1 and 3, and thus lead to different ground state vibrational 1D wavefunctions |0(fragment)〉 (stretching or bending) and hence result in different vibrational populations Nν of the gas phase SiH2 according with eqn (4).
For each of the three modes (two stretching and one bending), 1D wavefunctions of the gas phase SiH2 for 30 vibrational states from ν0 up to ν29 were calculated. These 1D wavefunctions were then numerically integrated with the corresponding ground state 1D wavefunctions |ν = 0〉 of the fragment in question.
The total vibrational population Nν1,ν2,ν3 in this approximation is then given by a product
![]() | (10) |
where and
are obtained using the stretching harmonic oscillators wavefunctions |ν1〉 and |ν3〉, while
is obtained using the corresponding bending harmonic oscillator wavefunction |ν2〉. The independent treatment of the two stretching populations is partly justified by the local mode character of the vibrational degrees of freedom of SiH2 due to the near 90° bond angle (see, e.g. Jensen110 and Clark et al.72). The asymmetric vibrational modes of SiH2 (B2 in C2v) are non populated in this 1D approximation. This is because for the parallel nature of the Franck–Condon transitions from the ground vibrational state of disilane, which is fully symmetric (A1), only symmetric states of SiH2 give rise to non-zero integrals in eqn (4). For example, the vibrational population of the (2,1,0) state (2 stretching and 1 bending quanta of A1) of the GM fragment is obtained as a product of
,
and
, respectively, resulting in N(ν2,ν1,ν0) = 0.088, while the population of the B2-type (0,0,1) vibratitonal state is assumed to be zero. The populations Nν1,ν2,ν3 are pre-calculated for each vibrational state (ν1,ν2,ν3) of SiH2 and added to the CATS State file to be used in non-LTE simulations (see below, Section 4.3). All the vibrational populations computed and used as part of this work are provided in the ESI.†
![]() | ||
Fig. 3 The 1D bending mode harmonic wavefunctions for the five SiH2 fragments described in Table 1 as a function of bond angle in degrees. Ground state fragment (black) compared with the LTE |ν = 0〉 (blue), |ν = 1〉 (green), |ν = 2〉 (red). |
![]() | ||
Fig. 4 The 1D stretching mode harmonic wavefunctions for the five SiH2 fragments described in Table 1 as a function of bond length in Angstrom. Ground state fragment (black) compared with the LTE |ν = 0〉 (blue), |ν = 1〉 (green), |ν = 2〉 (red). |
The shapes and positions of the curves in Fig. 3 and 4 match with the parameters from Table 1. The larger ω used for the bending modes manifests itself as wider curves, while the black curves are all centred around the equilibrium bond angles and length listed in Table 1.
The similarity between the populations of the LM-L, TS-L and GM fragments is expected owing to their similar structural parameters. It is interesting to see the most populated vibrational levels of LM-R and TS-R are always lower than the vibrational levels of LM-L, TS-L and GM.
Fig. 6 and 7 show the 1000 cm−1 and 2000 cm−1 (10 μm and 5 μm) bands for the SiH2 absorption spectrum, respectively, simulated using the non-LTE densities from Fig. 5 for all five cases considered and compared to the LTE scenario assuming the (rotational) temperature of T = 296 K and using the CATS line list. The strongest bands are indicated using different colours. Fig. 6 focuses on the 1000 cm−1 band. Most of the non-LTE spectra contain bending hot bands (020)–(010) and (020)–(010), which are stronger than the fundamental band (000)–(000). It can be seen that the P and R branches of the non-LTE spectra are shifted to lower wavenumbers in the GM, LM-L and TS-L spectra. In the TS-R and LM-R spectra the bands are not shifted, with the LM-R spectrum having only the fundamental (010)–(000) band visible.
The plots in Fig. 7 show the 2000 cm−1 band in the region of the polyad (100)/(020)/(001) for the five fragments, with the strongest fundamental band (001)–(000). The non-LTE intensities of the hot bands (011)–(010), (200)–(020) are found to be comparable to the intensities of the (001)–(000) band. The Q branch is clearly shifted for the GM, TS-L and LM-L molecules. The band is less shifted for the TS-R and LM-R fragments, but owing to the increased similarity between the fragment and molecular structures with TS-R and LM-R this is to be expected. Only the main polyad system (100)/(020)/(001) is visible for the TS-R spectrum (indicated as (001)–(000) in Fig. 7).
With the equilibrium structures of the TS-R and LM-R fragments being similar to the equilibrium structure of SiH2, their non-LTE spectra are expected to be a similar spectrum to LTE. Indeed, for the 1D harmonic approach their P, Q and R branches maintain the expected LTE intensities for both the 1000 cm−1 and 2000 cm−1 bands.
ϕSi2H6 = |0,0,0(SiH2)〉 |0,0,0,0,0,0,0,0,0,0,0,0,0,0,0(SiH4)〉. |
The vibrational populations are then modelled using the Franck–Condon integrals as follows:
Nν1,ν2,ν3 = |〈0,0,0(fragment)|ν1, ν2, ν3(g.ph.)〉|2, | (11) |
In order to simplify the 3D integration in eqn (11), the variational wavefunction |0,0,0(fragment)〉 is obtained using the same vibrational basis set as the variational solution of the gas phase SiH2. By taking advantage of the compatibility of the orthogonality of the basis sets, the Franck–Condon factors are then given by as a sum of products
![]() | (12) |
TROVE uses optimized non-standard vibrational basis sets, generated numerically by solving 1D Schrödinger equations for realistic 1D potentials.74 This procedure allows producing compact basis functions optimized for a specific problem. In our case, the PESs of the corresponding five fragments and of the gas phase SiH2 are different and therefore the generated basis sets would be different and even not orthogonal. We therefore implemented a feature in TROVE allowing to read and use externally generated basis functions. Of course all relevant calculation setups must be compatible, including the numerical grids used for the stretching and bending modes and their sizes. Using foreign basis sets certainly degrades their quality. However, since we are only interested in fragments' ground state wavefunctions, this degradation can be mitigated by including enough basis functions. Our typical 1D basis sets contain 12–24 functions (see details below), which should be more than enough to obtain a converged ground state solution even with non-optimized basis sets.
VSi2H6 = VSiH2(r1, r2, α) + Vrest, |
VSiH2(r1, r2, α) = fstr(r1) + fstr(r2) + fbnd(α). | (13) |
![]() | (14) |
![]() | (15) |
Parameter | GM | TS-L | TS-R | LM-L | LM- R |
---|---|---|---|---|---|
r e, Å | 1.48218 | 1.47971 | 1.49888 | 1.47982 | 1.51611 |
α e, deg | 108.701 | 108.938 | 103.268 | 108.264 | 92.287 |
a, Å−1 | 1.29065 | 1.29110 | 1.22364 | 1.27911 | 1.27050 |
g 1 | 0 | 0 | 0 | 0 | −8132.4642 |
g 2 | 80964.377 | 90039.635 | 101613.84 | 93608.413 | 168296.56 |
g 3 | −117331.83 | −126539.98 | −159054.10 | −129420.43 | −299843.15 |
g 4 | 139212.98 | 168090.27 | 182459.45 | 305247.60 | 324286.87 |
f 1 | 0 | 0 | 0 | 0 | −793.07685 |
f 2 | 43617.985 | 43796.459 | 45462.423 | 22217.055 | 19382.230 |
f 3 | −6594.4391 | −6937.7586 | −8917.729 | −3824.8354 | −2298.2220 |
f 4 | 2809.6655 | 2714.6696 | 3750.0966 | 1407.7267 | 2233.5205 |
Fig. 8 illustrates the ab initio PESs of different fragments as 1D cuts for the stretching and bending modes compared to the corresponding cuts of the gas phase SiH2 molecule. The bending cuts have especially different shapes with shifts to larger equilibrium angles and much steeper PESs. The differences in the stretching cuts are less pronounced. These features are important for the non-LTE behaviour of the corresponding excited states populations, with the bending degree of freedom to have stronger non-LTE character than stretching.
P= 2(ν1 + ν3) + ν2 ≤ 24, | (16) |
For the gas phase SiH2 calculations we employed the empirically refined PES by Clark et al.72 For the five SiH2 fragments only the vibrational ground state wavefunctions |0,0,0〉 were computed using the same setup and utilizing the basis functions from the g.ph. calculations as described above, but for the fragments' ab initio PESs from eqn (13).
The 3D results generally agree with the 1D case. As an example, Fig. 9 shows the TS-L case: the bending primitive wavefunctions of the gas phase SiH2 are compared to the ν = 0 bending wavefunction generated for the PES of TS-L. Fig. 10 shows the corresponding non-LTE vibrational populations as a function of the corresponding energies for the same TS-L scenario of the 3D vibrational populations. These populations are very different from the Boltzmann distribution, also shown on this figure for the SiH2 vibrational states at T = 1000 K, which exhibits an exponential decay with ν = 0 at its maximum. It nicely demonstrates that it would not be possible to associate a single vibrational temperature for this bell-shaped distribution.
![]() | ||
Fig. 9 3D Bending primitive wavefunctions of SiH2 compared to the ground state 3D bending wavefunctions for Fragment TS-L. Calculated with TROVE. |
![]() | (17) |
In order to simulate absorption spectra of the gas phase SiH2 assuming a non-LTE vibrational populations Nν1,ν2,ν3, we use the line positions and Einstein A coefficients from the ExoMol CATS line list72 employing the ExoCross program.92 The ExoMol line lists are formatted as two files, a States file and a Transition file. It is described extensively elsewhere73 and in this paper we shall only discuss how the ExoMol format has been adapted for use in non-LTE situations.
To adapt the CATS States file72 for non-LTE applications an additional ‘density’ column Nν1,ν2,ν3 was added as a final column, see an extract from the States file in Table 6. This column contains the weightings to the transition probabilities as populations of the vibrational levels occupied by the gas phase SiH2. The ‘density’ column is specific for the calculation of line lists for non-LTE molecules is not routinely included into the ExoMol States files.
This column is read by ExoCross and used to give the population weighting to each line intensity. Evaluation of the non-LTE population NJ,k,ν(T) given for each ro-vibrational state as in eqn (1) is based on the knowledge of the corresponding vibrational state (ν1,ν2,ν3) as well as the rotational energy contribution ẼrotJ,k. Therefore for this approach to work it is mandatory for all ro-vibrational states to be vibrationally assigned in order to be able to subtract the vibrational contribution Ẽvvib from the total energy according with eqn (2). All vibrational quantum numbers (ν1,ν2,ν3) are not required, only a vibrational index indicating the vibrational state in question. In our model, all asymmetric vibrational states (B2) are not populated due to the zero overlap with the ground state of the A1 symmetry in eqn (11), as part of the completely vertical Franck–Condon approximation.
The stark differences of the non-LTE spectra offer an ability for experiment to distinguish between SiH2 molecules produced from fragmenting Si2H6, and even to indicate dissociation channels involved.
In the RENNER calculations, the size of the basis set originally used in Yurchenko et al.118 was reduced in order to be able to increase the rotational excitations. The main purpose of this exercise is to show a qualitative impact of the non-LTE populations on the spectral shape of the electronic band and not so much the quality of the line positions, and therefore a smaller basis set is justified. We used 18 and 12 bending basis functions for the and à electronic states, respectively, for every |k| block (where |k| ≤ J). The
1A1 electronic state basis set included NA = 12 stretching functions of the A1 symmetry and NB = 10 stretching functions of the B2 symmetry. For the Ã1B1 state NA = 10 and NB = 8 stretching functions were used. These stretching functions were constructed from the Morse oscillator functions |n1〉|n3〉 with n1 + n3 ≤ Nstretch = 12.
A rovibronic line list for the Ã1B1–1A1 of SiH2 was generated covering the rotational excitations up to Jmax = 15 with the lower state energies (
) truncated at hc·25
000 cm−1 and the upper state energies truncated at hc·28
000 cm−1.
For the non-LTE simulations we used the 1D vibrational population model with the structural parameters corresponding to TS-L from Table 3. Fig. 12 shows a non-LTE electronic spectrum of SiH2 in the region of the band Ã(0,2,0)←(0,0,0), assuming the rotational temperature Trot = 500 K, compared to an LTE spectrum of T = 500 K. The non-LTE spectrum contains the hot band Ã(0,3,0)←
(0,1,0) which can be used to identify the non-LTE character of the system. The rovibronic line 101←110 belonging to this band used in a number of experimental studies involving SiH2 as a reaction product71,89,90 to estimate reaction rates. It is common for such studies to assume the Boltzmann equilibrium at different stages of the analysis of the measurements. For example, the partition function of SiH2 is required to estimate the number density of SiH2 in its lower, ground electronic state,70,71,119 which is directly affected by the LTE assumption. As our calculations show, the number densities of SiH2 as a reaction product can vary significantly under non-LTE depending on the reaction pathway and impact the experimental rates.
![]() | ||
Fig. 12 An Ã1B1–![]() |
There are two stable isomers of disilane, a local minimum structure and a global minimum structure, with a third transition state structure also known. Non-LTE spectra of SiH2 corresponding to dissociation of disilane from different sides of the three disilane isomer were computed. We have shown that the non-LTE spectra of SiH2 are different in most cases. This is important as the spectrum of SiH2 is used to monitor the quantity present in a reaction as a means to track the progress of SiH2 + SiH4 → Si2H6 and Si2H6 → SiH2 + SiH4 when calculating the corresponding rate constant. If Si2H6 is decomposing at a rate slower than it is being formed, then tracking the quantity of SiH2 can give a rate constant that is not reflective of the speed of reaction, and merely an indication of the equilibrium balance of the two species SiH2 and Si2H6.
In two approaches considered, 1D and the 3D, we assume that the rotational degrees of freedom are equillibrated quickly once the dissociation from disilane occurs, hence we use the Boltzmann distribution for the rotational degrees of freedom. We also assume that the SiH2 fragment during the instantaneous dissociation and is fully decoupled from the rest of the Si2H6 molecule, i.e. can be described by a 3D wavefunction in its lowest, relaxed vibrational configuration and has the same structural parameters as the Si2H6 molecule.
We have shown that the non-LTE spectra of SiH2 can be calculated by the new TROVE methodology and existing ExoMol line list, and it compares well to the simpler 1D harmonic approximation published previously. The method could be applied to the non-LTE spectroscopy of other small molecules including SiH4, which has not been explored here.
We have also shown that despite the many approximations used in the 1D approximation (separability of the modes, Harmonic approximate etc.), the results compare well to the results obtained using the full 3D approach. This lends confidence in using the simplified but robust 1D approach in similar non-LTE studies, as e.g. we have used to model the CO non-LTE spectra,6 recently, which are planning to explore in the future.
The methods described here can be used model the intensity distribution of the reaction products and to ascertain from what molecule the SiH2 dissociated from. The equilibrium structured parameters (bond lengths and angles) can be treated as effective parameters to be adjusted to reproduce the experimental spectra.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1cp00839k |
This journal is © the Owner Societies 2021 |