skip to main content

Title: Intersystem crossing in tunneling regime: T 1 → S 0 relaxation in thiophosgene
The T 1 excited state relaxation in thiophosgene has attracted much attention as a relatively simple model for the intersystem crossing (ISC) transitions in polyatomic molecules. The very short (20–40 ps) T 1 lifetime predicted in several theoretical studies strongly disagrees with the experimental values (∼20 ns) indicating that the kinetics of T 1 → S 0 ISC is not well understood. We use the nonadiabatic transition state theory (NA-TST) with the Zhu–Nakamura transition probability and the multireference perturbation theory (CASPT2) to show that the T 1 → S 0 ISC occurs in the quantum tunneling regime. We also introduce a new zero-point vibrational energy correction scheme that improves the accuracy of the predicted ISC rate constants at low internal energies. The predicted lifetimes of the T 1 vibrational states are between one and two orders of magnitude larger than the experimental values. This overestimation is attributed to the multidimensional nature of quantum tunneling that facilitates ISC transitions along the non-minimum energy path and is not accounted for in the one-dimensional NA-TST.
Award ID(s):
Publication Date:
Journal Name:
Physical Chemistry Chemical Physics
Page Range or eLocation-ID:
5500 to 5508
Sponsoring Org:
National Science Foundation
More Like this
  1. The infrared spectra of jet-cooled methyl anthranilate (MA) and the MA–H 2 O complex are reported in both S 0 and S 1 states, recorded using fluorescence-dip infrared (FDIR) spectroscopy under jet-cooled conditions. Using a combination of local mode CH stretch modeling and scaled harmonic vibrational character, a near-complete assignment of the infrared spectra is possible over the 1400–3700 cm −1 region. While the NH stretch fundamentals are easily observed in the S 0 spectrum, in the S 1 state, the hydrogen bonded NH stretch shift is not readily apparent. Scaled harmonic calculations predict this fundamental at just below 2900 cm −1 with an intensity around 400 km mol −1 . However, the experimental spectrum shows no evidence of this transition. A local mode theory is developed in which the NH stretch vibration is treated adiabatically. Minimizing the energy of the corresponding stretch state with one quantum of excitation leads to a dislocation of the H atom where there is equal sharing between N and O atoms. The sharing occurs as a result of significant molecular arrangement due to strong coupling of this NH stretch to other internal degrees of freedom and in particular to the contiguous HNC bend. Amore »two-dimensional model of the coupling between the NH stretch and this bend highlights important nonlinear effects that are not captured by low order vibrational perturbation theory. In particular, the model predicts a dramatic dilution of the NH stretch oscillator strength over many transitions spread over more than 1000 cm −1 , making it difficult to observe experimentally.« less
  2. The kinetic energy dependent reactions of Re + with SO 2 were studied with guided ion beam tandem mass spectrometry. ReO + , ReO 2 + , and OReS + species were observed as products, all in endothermic reactions. Modeling of the kinetic energy dependent cross sections yields 0 K bond dissociation energies (BDEs, in eV) of 4.78 ± 0.06 (Re + –O), 5.75 ± 0.02 (Re + –O 2 ), and 4.35 ± 0.14 (Re + –SO). The latter two values can be combined with other information to derive the additional values 6.05 ± 0.05 (ORe + –O) and 4.89 ± 0.19 (ORe + –S). BDEs of ReO + and ReO 2 + agree with literature values whereas the values for OReS + are the first measurements. The former result is obtained even though formation of ground state ReO + is spin-forbidden. Quantum mechanical calculations at the B3LYP level of theory with a def2-TZVPPD basis set yield results that agree reasonably well with experimental values. Additional calculations at the BP86 and CCSD(T) levels of theory using def2-QZVPPD and aug-cc-pVxZ (x = T, Q, and 5) basis sets were performed to compare thermochemistry with experiment to determine that ReO 2more »+ rather than the isobaric ReS + is formed. Product ground states are 3 Δ 3 (ReO + ), 3 B 1 (OReO + ), 5 Π −1 (ReS + ), and 3 A′′(OReS + ) after including empirical spin–orbit corrections, which means that formation of ground state products is spin-forbidden for all three product channels. The potential energy surfaces for the ReSO 2 + system were also explored at the B3LYP/def2-TZVPPD level and exhibited no barriers in excess of the endothermicities for all products. BDEs for rhenium oxide and sulfide diatomics and triatomics are compared and discussed. The present result for formation of ReO + is compared to that for formation of ReO + in the reactions of Re + + O 2 and CO, where the former system exhibited interesting dual cross section features. Results are consistent with the hypothesis that the distinction of in-plane and out-of-plane C S symmetry in the triatomic systems might be the explanation for the two endothermic features observed in the Re + + O 2 reaction.« less
  3. Resonant tunneling diodes (RTDs) have come full-circle in the past 10 years after their demonstration in the early 1990s as the fastest room-temperature semiconductor oscillator, displaying experimental results up to 712 GHz and fmax values exceeding 1.0 THz [1]. Now the RTD is once again the preeminent electronic oscillator above 1.0 THz and is being implemented as a coherent source [2] and a self-oscillating mixer [3], amongst other applications. This paper concerns RTD electroluminescence – an effect that has been studied very little in the past 30+ years of RTD development, and not at room temperature. We present experiments and modeling of an n-type In0.53Ga0.47As/AlAs double-barrier RTD operating as a cross-gap light emitter at ~300K. The MBE-growth stack is shown in Fig. 1(a). A 15-μm-diam-mesa device was defined by standard planar processing including a top annular ohmic contact with a 5-μm-diam pinhole in the center to couple out enough of the internal emission for accurate free-space power measurements [4]. The emission spectra have the behavior displayed in Fig. 1(b), parameterized by bias voltage (VB). The long wavelength emission edge is at  = 1684 nm - close to the In0.53Ga0.47As bandgap energy of Ug ≈ 0.75 eV at 300 K.more »The spectral peaks for VB = 2.8 and 3.0 V both occur around  = 1550 nm (h = 0.75 eV), so blue-shifted relative to the peak of the “ideal”, bulk InGaAs emission spectrum shown in Fig. 1(b) [5]. These results are consistent with the model displayed in Fig. 1(c), whereby the broad emission peak is attributed to the radiative recombination between electrons accumulated on the emitter side, and holes generated on the emitter side by interband tunneling with current density Jinter. The blue-shifted main peak is attributed to the quantum-size effect on the emitter side, which creates a radiative recombination rate RN,2 comparable to the band-edge cross-gap rate RN,1. Further support for this model is provided by the shorter wavelength and weaker emission peak shown in Fig. 1(b) around = 1148 nm. Our quantum mechanical calculations attribute this to radiative recombination RR,3 in the RTD quantum well between the electron ground-state level E1,e, and the hole level E1,h. To further test the model and estimate quantum efficiencies, we conducted optical power measurements using a large-area Ge photodiode located ≈3 mm away from the RTD pinhole, and having spectral response between 800 and 1800 nm with a peak responsivity of ≈0.85 A/W at  =1550 nm. Simultaneous I-V and L-V plots were obtained and are plotted in Fig. 2(a) with positive bias on the top contact (emitter on the bottom). The I-V curve displays a pronounced NDR region having a current peak-to-valley current ratio of 10.7 (typical for In0.53Ga0.47As RTDs). The external quantum efficiency (EQE) was calculated from EQE = e∙IP/(∙IE∙h) where IP is the photodiode dc current and IE the RTD current. The plot of EQE is shown in Fig. 2(b) where we see a very rapid rise with VB, but a maximum value (at VB= 3.0 V) of only ≈2×10-5. To extract the internal quantum efficiency (IQE), we use the expression EQE= c ∙i ∙r ≡ c∙IQE where ci, and r are the optical-coupling, electrical-injection, and radiative recombination efficiencies, respectively [6]. Our separate optical calculations yield c≈3.4×10-4 (limited primarily by the small pinhole) from which we obtain the curve of IQE plotted in Fig. 2(b) (right-hand scale). The maximum value of IQE (again at VB = 3.0 V) is 6.0%. From the implicit definition of IQE in terms of i and r given above, and the fact that the recombination efficiency in In0.53Ga0.47As is likely limited by Auger scattering, this result for IQE suggests that i might be significantly high. To estimate i, we have used the experimental total current of Fig. 2(a), the Kane two-band model of interband tunneling [7] computed in conjunction with a solution to Poisson’s equation across the entire structure, and a rate-equation model of Auger recombination on the emitter side [6] assuming a free-electron density of 2×1018 cm3. We focus on the high-bias regime above VB = 2.5 V of Fig. 2(a) where most of the interband tunneling should occur in the depletion region on the collector side [Jinter,2 in Fig. 1(c)]. And because of the high-quality of the InGaAs/AlAs heterostructure (very few traps or deep levels), most of the holes should reach the emitter side by some combination of drift, diffusion, and tunneling through the valence-band double barriers (Type-I offset) between InGaAs and AlAs. The computed interband current density Jinter is shown in Fig. 3(a) along with the total current density Jtot. At the maximum Jinter (at VB=3.0 V) of 7.4×102 A/cm2, we get i = Jinter/Jtot = 0.18, which is surprisingly high considering there is no p-type doping in the device. When combined with the Auger-limited r of 0.41 and c ≈ 3.4×10-4, we find a model value of IQE = 7.4% in good agreement with experiment. This leads to the model values for EQE plotted in Fig. 2(b) - also in good agreement with experiment. Finally, we address the high Jinter and consider a possible universal nature of the light-emission mechanism. Fig. 3(b) shows the tunneling probability T according to the Kane two-band model in the three materials, In0.53Ga0.47As, GaAs, and GaN, following our observation of a similar electroluminescence mechanism in GaN/AlN RTDs (due to strong polarization field of wurtzite structures) [8]. The expression is Tinter = (2/9)∙exp[(-2 ∙Ug 2 ∙me)/(2h∙P∙E)], where Ug is the bandgap energy, P is the valence-to-conduction-band momentum matrix element, and E is the electric field. Values for the highest calculated internal E fields for the InGaAs and GaN are also shown, indicating that Tinter in those structures approaches values of ~10-5. As shown, a GaAs RTD would require an internal field of ~6×105 V/cm, which is rarely realized in standard GaAs RTDs, perhaps explaining why there have been few if any reports of room-temperature electroluminescence in the GaAs devices. [1] E.R. Brown,et al., Appl. Phys. Lett., vol. 58, 2291, 1991. [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [2] M. Feiginov et al., Appl. Phys. Lett., 99, 233506, 2011. [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [3] Y. Nishida et al., Nature Sci. Reports, 9, 18125, 2019. [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [4] P. Fakhimi, et al., 2019 DRC Conference Digest. [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018).« less
  4. Highly reactive benzyl radicals are generated by electron dissociative attachment to benzyl chloride doped into a neon–hydrogen–helium discharge and immediately cooled to T rot = 15 K in a high density, supersonic slit expansion environment. The sub-Doppler spectra are fit to an asymmetric-top rotational Hamiltonian, thereby yielding spectroscopic constants for the ground ( v = 0) and first excited ( v = 1, ν 3 , ν 4 ) vibrational levels of the ground electronic state. The rotational constants obtained for the ground state are in good agreement with previous laser induced fluorescence measurements (LIF), with vibrational band origins ( ν 3 = 3073.2350 ± 0.0006 cm −1 , ν 4 = 3067.0576 ± 0.0006 cm −1 ) in agreement with anharmonically corrected density functional theory calculations. To assist in detection of benzyl radical in the interstellar medium, we have also significantly improved the precision of the ground state rotational constants through combined analysis of the ground state IR and LIF combination differences. Of dynamical interest, there is no evidence in the sub-Doppler spectra for tunneling splittings due to internal rotation of the CH 2 methylene subunit, which implies a significant rotational barrier consistent with partial double bond character inmore »the CC bond. This is further confirmed with high level ab initio calculations at the CCSD(T)-f12b/ccpVdZ-f12 level, which predict a zero-point energy corrected barrier to internal rotation of Δ E tun ≈ 11.45 kcal mol −1 or 4005 cm −1 . In summary, the high-resolution infrared spectra are in excellent agreement with simple physical organic chemistry pictures of a strongly resonance-stabilized benzyl radical with a nearly rigid planar structure due to electron delocalization around the aromatic ring.« less
  5. The kinetic energy dependences of the reactions of Pt + ( 2 D 5/2 ) with SO 2 were studied using a guided ion beam tandem mass spectrometer and theory. The observed cationic products are PtO + and PtSO + , with small amounts of PtS + , all formed in endothermic reactions. Modeling the kinetic energy dependent product cross sections allows determination of the product bond dissociation energies (BDEs): D 0 (Pt + –O) = 3.14 ± 0.11 eV, D 0 (Pt + –S) = 3.68 ± 0.31 eV, and D 0 (Pt + –SO) = 3.03 ± 0.12 eV. The oxide BDE agrees well with more precise literature values, whereas the latter two results are the first such measurements. Quantum mechanical calculations were performed for PtO + , PtS + , PtO 2 + , and PtSO + at the B3LYP and coupled-cluster with single, double, and perturbative triple [CCSD(T)] levels of theory using the def2-XZVPPD (X = T, Q) and aug-cc-pVXZ (X = T, Q, 5) basis sets and complete basis set extrapolations. These theoretical BDEs agree well with the experimental values. After including empirical spin–orbit corrections, the product ground states are determined as PtO + ( 4 Σ 3/2more »), PtS + ( 4 Σ 3/2 ), PtO 2 + ( 2 Σ g + ), and PtSO + ( 2 A′). Potential energy profiles including intermediates and transition states for each reaction were also calculated at the B3LYP/def2-TZVPPD level. Periodic trends in the thermochemistry of the group 9 metal chalcogenide cations are compared, and the formation of PtO + from the Pt + + SO 2 reaction is compared with those from the Pt + + O 2 , CO 2 , CO, and NO reactions.« less