Synthesis, Characterization, and Computational Studies on Gallium(III) and Iron(III) Complexes with a Pentadentate Macrocyclic bis-Phosphinate Chelator and Their Investigation As Molecular Scaffolds for 18F Binding

With the aim of obtaining improved molecular scaffolds for 18F binding to use in PET imaging, gallium(III) and iron(III) complexes with a macrocyclic bis-phosphinate chelator have been synthesized and their properties, including their fluoride binding ability, investigated. Reaction of Bn-tacn (1-benzyl-1,4,7-triazacyclononane) with paraformaldehyde and PhP(OR)2 (R = Me or Et) in refluxing THF, followed by acid hydrolysis, yields the macrocyclic bis(phosphinic acid) derivative, H2(Bn-NODP) (1-benzyl-4,7-phenylphosphinic acid-1,4,7-triazacyclononane), which is isolated as its protonated form, H2(Bn-NODP)·2HCl·4H2O, at low pH (HClaq), its disodium salt, Na2(Bn-NODP)·5H2O at pH 12 (NaOHaq), or the neutral H2(Bn-NODP) under mildly basic conditions (Et3N). A crystal structure of H2(Bn-NODP)·2HCl·H2O confirmed the ligand's identity. The mononuclear [GaCl(Bn-NODP)] complex was prepared by treatment of either the HCl or sodium salt with Ga(NO3)3·9H2O or GaCl3, while treatment of H2(Bn-NODP)·2HCl·4H2O with FeCl3 in aqueous HCl gives [FeCl(Bn-NODP)]. The addition of 1 mol. equiv of aqueous KF to these chloro complexes readily forms the [MF(Bn-NODP)] analogues. Spectroscopic analysis on these complexes confirms pentadentate coordination of the doubly deprotonated (bis-phosphinate) macrocycle via its N3O2 donor set, with the halide ligand completing a distorted octahedral geometry; this is further confirmed through a crystal structure analysis on [GaF(Bn-NODP)]·4H2O. The complex adopts the geometric isomer in which the phosphinate arms are coordinated unsymmetrically (isomer 1) and with the stereochemistry of the three N atoms of the tacn ring in the RRS configuration, denoted (N)RRS, and the phosphinate groups in the RR stereochemistry, denoted (P)RR, (isomer 1/RR), together with its (N)SSR (P)SS enantiomer. The greater thermodynamic stability of isomer 1/RR over the other possible isomers is also indicated by density functional theory (DFT) calculations. Radiofluorination experiments on the [MCl(Bn-NODP)] complexes in partially aqueous MeCN/NaOAcaq (Ga) or EtOH (Ga or Fe; i.e. without buffer) with 18F- target water at 80 °C/10 min lead to high radiochemical incorporation (radiochemical yields 60-80% at 1 mg/mL, or ∼1.5 μM, concentration of the precursor). While the [Fe18F(n-NODP)] is unstable (loss of 18F-) in both H2O/EtOH and PBS/EtOH (PBS = phosphate buffered saline), the [Ga18F(Bn-NODP)] radioproduct shows excellent stability, RCP = 99% at t = 4 h (RCP = radiochemical purity) when formulated in 90%:10% H2O/EtOH and ca. 95% RCP over 4 h when formulated in 90%:10% PBS/EtOH. This indicates that the new "GaIII(Bn-NODP)" moiety is a considerably superior fluoride binding scaffold than the previously reported [Ga18F(Bn-NODA)] (Bn-NODA = 1-benzyl-4,7-dicarboxylate-1,4,7-triazacyclononane), which undergoes rapid and complete hydrolysis in PBS/EtOH (refer to Chem. Eur. J. 2015, 21, 4688-4694).

A study of the thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) with atomic chlorine (Cl) in the absence and presence of water, using electronic structure methods

The thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) + atomic chlorine (Cl) were investigated in the absence and presence of a single water molecule, using electronic structure methods. Stationary points on each reaction surface were located using density functional theory (DFT) with the M06-2X functional with aug-cc-pVDZ (aVDZ) and aug-cc-pVTZ (aVTZ) basis sets. Then fixed point calculations were carried out using the UM06-2X/aVTZ optimised stationary point geometries, with aug-cc-pVnZ basis sets (n = T and Q), using the coupled cluster method [CCSD(T)], as well as the domain-based local pair natural orbitals coupled cluster [DLPNO-UCCSD(T)] approach. Four reaction channels are possible, formation of (A) CH₃SCH₂ + HCl, (B) CH₃S + CH₃Cl, (C) CH₃SCl + CH₃, and (C') CH₃S(Cl)CH₃. The results show that, in the absence of water, channels A and C' are the dominant channels. In the presence of water, the calculations show that the reaction mechanisms for A and C formation change significantly. Channel A occurs via submerged TSs and is expected to be rapid. Channel B occurs via TSs which present significant energy barriers indicating that this channel is not significant in the presence of water relative to CH₃SCH₂ + HCl and DMS·Cl adduct formation, as is the case in the absence of water. Channel C was not considered as it is endothermic in the absence of water. In the presence of water, pathways which proceed via (a) DMS·H₂O + Cl, (b) Cl·H₂O + DMS and (c) DMS·Cl + H₂O were considered. It was found that under tropospheric conditions, reactions via pathway (b) are of minor importance relative to those that proceed via pathways (a) and (c). This study has shown that water changes the mechanisms of the DMS + Cl reactions significantly but the presence of water is not expected to affect the overall reaction rate coefficient under atmospheric conditions as the DMS + Cl reaction has a rate coefficient at room temperature close to the collisional limit.

Neutral and Cationic Complexes of Silicon(IV) Halides with Phosphine Ligands

The reaction of SiI₄ and PMe₃ in n-hexane produced the yellow salt, [SiI₃(PMe₃)₂]I, confirmed from its X-ray structure, containing a trigonal bipyramidal cation with trans-phosphines. This contrasts with the six-coordination found in (the known) trans-[SiX₄(PMe₃)₂] (X = Cl, Br) complexes. The diphosphines o-C₆H₄(PMe₂)₂ and Et₂P(CH₂)₂PEt₂ form six-coordinate cis-[SiI₄(diphosphine)], which were also characterized by X-ray crystallography, multinuclear NMR, and IR spectroscopy. Reaction of trans-[SiX₄(PMe₃)₂] (X = Cl, Br) with Na[BAr^F] (BAr^F = [B{3,5-(CF₃)₂C₆H₃}₄]) produced five-coordinate [SiX₃(PMe₃)₂][BAr^F], but while Me₃SiO₃SCF₃ also abstracted chloride from trans-[SiCl₄(PMe₃)₂], the reaction products were six-coordinate complexes [SiCl₃(PMe₃)₂(OTf)] and [SiCl₂(PMe₃)₂(OTf)₂] with the triflate coordinated. X-ray crystal structures were obtained for [SiCl₃(PMe₃)₂][BAr^F] and [SiCl₂(PMe₃)₂(OTf)₂]. The charge distribution across the silicon species was also examined by natural bond orbital (NBO) analyses of the computed density functional theory (DFT) wavefunctions. For the [SiX₄(PMe₃)₂] and [SiX₃(PMe₃)₂]⁺ complexes, the positive charge on Si decreases and the negative charge on X decreases going from X = F to X = I. Upon going from [SiX₄(PMe₃)₂] to [SiX₃(PMe₃)₂]⁺, i.e., removal of X^–, there is an increase in positive charge on Si and a decrease in negative charge on the X centers (except for the case X = F). The positive charge on P shows a slight decrease.

Unusual neutral and monocationic silicon(IV) iodide complexes with soft phosphine coordination are formed in hydrocarbon solvents, while treatment of [SiCl₄(PMe₃)₂] with NaBAr^F forms the [SiCl₃(PMe₃)₂]⁺ cation; the electronic structures of [SiX₄(PMe₃)₂] and [SiX₃(PMe₃)₂]⁺ (X = F, Cl, Br, I) are probed using DFT calculations.

An experimental and theoretical characterization of the electronic structure of doubly ionised disulfur

Using time-of-flight multiple electron and ion coincidence techniques in combination with a helium gas discharge lamp and synchrotron radiation, the double ionisation spectrum of disulfur (S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2) and the subsequent fragmentation dynamics of its dication are investigated. The S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 sample was produced by heating mercury sulfide (HgS), whose vapour at a suitably chosen temperature consists primarily of two constituents: S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 and atomic Hg. A multi-particle-coincidence technique is thus particularly useful for retrieving spectra of S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 from ionisation of the mixed vapour. The results obtained are compared with detailed calculations of the electronic structure and potential energy curves of S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2^{2+}$$\end{document}22+ which are also presented. These computations are carried out using configuration interaction methodology. The experimental results are interpreted with and strongly supported by the computational results.

Single photon double and triple ionization of allene

Double and triple ionization of allene are investigated using electron-electron, ion-ion, electron-electron-ion and electron-electron-ion-ion (ee, ii, eei, eeii) coincidence spectroscopies at selected photon energies. The results provide supporting evidence for a previously proposed roaming mechanism in H₃⁺ formation by double ionization. The lowest vertical double ionization energy is found to be 27.9 eV, while adiabatic double ionization is not accessed by vertical ionization at the neutral geometry. The triple ionization energy is found to be close to 50 eV in agreement with theoretical predictions. The doubly charged parent ion is stable up to about 2 eV above the threshold, after which dissociations by charge separation and by double charge retention occur with comparable intensities. Fragmentation to H⁺ + C₃H₃⁺ starts immediately above the threshold as a slow (metastable) decay with 130.5 ± 9.9 ns mean lifetime.

Pyramidal Dicationic Ge(II) Complexes with Homoleptic Neutral Pnictine Coordination: A Combined Experimental and Density Functional Theory Study

An unusual series of Ge(II) dicationic species with homoleptic phosphine and arsine coordination, [Ge(L)][OTf]₂, L = 3 × PMe₃, triphos (MeC(CH₂PPh₂)₃), triars (MeC(CH₂AsMe₂)₃), or κ³-tetraphos (P(CH₂CH₂PPh₂)₃) (OTf^- = O₃SCF₃^-) have been prepared by reaction of [GeCl₂(dioxane)] with L and 2 mol equiv of Me₃SiOTf in anhydrous CH₂Cl₂ (or MeCN for L = triars, triphos). X-ray crystal structures are reported for [Ge(PMe₃)₃][OTf]₂, [Ge(triars)][OTf]₂, and [Ge(κ³-tetraphos)][OTf]₂, confirming homoleptic P₃- or As₃-coordination at Ge(II) in each case and with the discrete OTf^- anions providing a charge balance. The Ge-P/As bond lengths are significantly shorter than those in neutral germanium(II) dihalide complexes with diphosphine or diarsine coordination. Solution NMR spectroscopic data indicate that the complexes are labile in solution. Using excess AsMe₃ and [GeCl₂(dioxane)] gives only the neutral product, [Ge(AsMe₂)₂(OTf)₂], the crystal structure of which shows four coordination at Ge(II), via two As donor atoms and an O atom from two κ¹-OTf^- ligands; further weak, long-range intermolecular interactions give a chain polymer. The electronic structure of the [Ge(PMe₃)₃]²⁺ dication has been investigated using density functional theory (DFT) calculations. The computed geometrical parameters for this dication are in good agreement with the experimental X-ray crystallographic values in [Ge(PMe₃)₃][OTf]₂. The results also indicate that the pyramidal arrangement of the [Ge(PMe₃)₃]²⁺ (computed P-Ge-P angle 96.8° at the B3LYP-D3 level) arises from a balance between electronic energy (E_elec) contributions, which favor a lower P-Ge-P angle, and nuclear-nuclear contributions (E_nn), which favor a higher P-Ge-P angle, to the total energy (E_TOT). An Atoms in Molecules (AIM) analysis reveals that one reason why E_elec decreases as the P-Ge-P angle decreases is because of C···H and H···H interactions between atoms on different CH₃ groups. The stability of the [Ge(PMe₃)₃]²⁺ dication is enhanced by the distribution of a significant part of the positive charge on Ge²⁺ to the atomic centers of the PMe₃ ligands. Similar results were obtained for [Ge(AsMe₃)₃][OTf]₂, showing the tris-AsMe₃ complex to be less stable compared to the PMe₃ analogue. Related calculations were also performed for the neutral [Ge(PMe₃)₂(OTf)₂] and [Ge(AsMe₃)₂(OTf)₂] complexes.

Comment on "Impact of water on the BrO + HO2 gas-phase reaction: mechanism, kinetics and products" by N. T. Tsona, S. Tang and L. Du, Phys. Chem. Chem. Phys., 2019, , 20296

The reaction, BrO + HO2 → HOBr + O2, is exothermic and can produce O2 in both its ground state (X[combining tilde]3∑g-) and its first excited state (ã1Δg). As a result, this reaction can proceed on both a singlet and a triplet potential energy surface. Recently, Tsona, Tang and Du published a paper entitled "Impact of water on the BrO + HO2 gas-phase reaction: mechanism, kinetics and products (Phys. Chem. Chem. Phys. 2019, 21, 20296-203072). The results of this work showed significant differences from those published earlier on this reaction by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569). Further calculations performed in this present work, combined with higher level calculations published by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569), demonstrate that the work of Tsona et al. is flawed because the integration grid size used in their lowest singlet and triplet calculations is too small, and a closed-shell wavefunction, rather than an open-shell wavefunction, has been used for the singlet surface. The major conclusion in the work of Tsona et al. that the lowest singlet and triplet channels are barrierless is shown to be incorrect. Also, the computed rate coefficients of Tsona et al. showed a positive temperature dependence, which is inconsistent with the experimentally observed negative temperature dependence, whereas the singlet rate coefficients computed by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569) showed a negative temperature dependence consistent with experiment.

Spectroscopy on the wing: Investigating possible differences in protein secondary structures in feather shafts of birds using Raman spectroscopy

The central shaft of a bird's flight feather bears most of the aerodynamic load during flight and exhibits some remarkable mechanical properties. The shaft comprises two parts, the calamus and the rachis. The calamus is at the base of the shaft, while the rachis is the longer upper part which supports the vanes. The shaft is composed of a fibrous outer cortex, and an inner foam-like core. Recent nanoindentation experiments have indicated that reduced modulus values, E_r, for the inner and outer regions of the cortex can vary, with the E_r values of the inner region slightly greater than those of the outer region. In this work, Raman spectroscopy is used to investigate the protein secondary structures in the inner and outer regions of the feather cortex. Analysis of the Amide I region of Raman spectra taken from four birds (Swan, Gull, Mallard and Kestrel) shows that the β-sheet structural component decreases between the inner and outer region, relative to the protein side-chain components. This finding is consistent with the proposal that E_r values are greater in the inner region than the outer region. This work has shown that Raman spectroscopy can be used effectively to study the change in protein secondary structure between the inner and outer regions of a feather shaft.

Tertiary Phosphine and Arsine Complexes of Phosphorus Pentafluoride: Synthesis, Properties, and Electronic Structures

The reaction of PMe₃ or PPh₃ with PF₅ in anhydrous CH₂Cl₂ or hexane forms the white, moisture-sensitive complexes [PF₅(PR₃)] (R = Me, Ph). Similar reactions involving the diphosphines o-C₆H₄(PR₂)₂ afford the complexes [PF₄{o-C₆H₄(PR₂)₂}][PF₆]. The X-ray structures of [PF₅(PR₃)] and [PF₄{o-C₆H₄(PMe₂)₂}][PF₆] show pseudo-octahedral fluorophosphorus centers. Multinuclear NMR spectra (¹H, ¹⁹F{¹H}, ³¹P{¹H}) show that in solution in CH₂Cl₂/CD₂Cl₂ the structures determined crystallographically are the only species present for [PF₅(PMe₃)] and [PF₄{o-C₆H₄(PMe₂)₂}][PF₆] but that [PF₅(PPh₃)] and [PF₄{o-C₆H₄(PPh₂)₂}][PF₆] exhibit reversible dissociation of the phosphine at ambient temperatures, although exchange slows at low temperatures. The complex ¹⁹F{¹H} and ³¹P{¹H} NMR spectra have been analyzed, including those of the cation [PF₄{o-C₆H₄(PMe₂)₂}]⁺, which is a second-order AA'XX'B₂M spin system. The unstable [PF₅(AsMe₃)], which decomposes in a few hours at ambient temperatures, has also been isolated and spectroscopically characterized; neither AsPh₃ nor SbEt₃ forms similar complexes. The electronic structures of the PF₅ complexes have been explored by DFT calculations. The DFT optimized geometries for [PF₅(PMe₃)], [PF₅(PPh₃)], and [PF₄{o-C₆H₄(PMe₂)₂}]⁺ are in good agreement with their respective crystal structure geometries. DFT calculations on the PF₅-L complexes reveal the P-L bond strength falls with L in the order PMe₃ > PPh₃ > AsMe₃, consistent with the experimentally observed stabilities, and in the PF₅-L complexes, electron transfer from L to PF₅ on forming these complexes also follows the order PMe₃ > PPh₃ ≈ AsMe₃.

Photoionization studies of reactive intermediates using synchrotron radiation

Photoionization with synchrotron radiation enables sensitive and selective monitoring of reactive intermediates in environments such as flames and plasmas.

A study of the Group 1 metal tetra-aza macrocyclic complexes [M(Me4cyclen)(L)]+ using electronic structure calculations

Metal-cyclen complexes have a number of important applications. However, the coordination chemistry between metal ions and cyclen-based macrocycles is much less well studied compared to their metal ion-crown ether analogues. This work, which makes a contribution to address this imbalance by studying complex ions of the type [M(Me₄cyclen)(L)]⁺, was initiated by results of an experimental study which prepared some Group 1 metal cyclen complexes, namely [Li(Me₄cyclen)(H₂O)][BAr^F] and [Na(Me₄cyclen)(THF)][BAr^F] and obtained their X-ray crystal structures [J. M. Dyke, W. Levason, M. E. Light, D. Pugh, G. Reid, H. Bhakhoa, P. Ramasami, and L. Rhyman, Dalton Trans., 2015, 44, 13853]. The lowest [M(Me₄cyclen)(L)]⁺ minimum energy structures (M = Li, Na, K, and L = H₂O, THF, DEE, MeOH, DCM) are studied using density functional theory (DFT) calculations. The geometry of each [M(Me₄cyclen)(L)]⁺ structure and, in particular, the conformation of L are found to be mainly governed by steric hindrance which decreases as the size of the ionic radius increases from Li⁺ → Na⁺ → K⁺. Good agreement of computed geometrical parameters of [Li(Me₄cyclen)(H₂O)]⁺ and [Na(Me₄cyclen)(THF)]⁺ with the corresponding geometrical parameters derived from the crystal structures [Li(Me₄cyclen)(H₂O)]⁺[BAr^F]^- and [Na(Me₄cyclen)(THF)]⁺[BAr^F]^- is obtained. Bonding analysis indicates that the stability of the [M(Me₄cyclen)(L)]⁺ structures originates mainly from ionic interaction between the Me₄cyclen/L ligands and the M⁺ centres. The experimental observation that [M(Me₄cyclen)(L)]⁺[BAr^F]^- complexes could be prepared in crystalline form for M⁺ = Li⁺ and Na⁺, but that experiments aimed at synthesising the corresponding K⁺, Rb⁺, and Cs⁺ complexes failed resulting in formation of [Me₄cyclenH][BAr^F] is investigated using DFT and explicitly correlated calculations, and explained by considering production of [Me₄cyclenH]⁺ by a hydrolysis reaction, involving traces of water, which competes with [M(Me₄cyclen)(L)]⁺ formation. [Me₄cyclenH]⁺ formation dominates for M⁺ = K⁺, Rb⁺, and Cs⁺ whereas formation of [M(Me₄cyclen)(L)]⁺ is energetically favoured for M⁺ = Li⁺ and Na⁺. The results indicate that the number and type of ligands, play a key role in stabilising the [M(Me₄cyclen)]⁺ complexes and it is hoped that this work will encourage experimentalists to prepare and characterise other [M(Me₄cyclen)(L)]⁺ complexes.

Can Cyclen Bind Alkali Metal Azides? A DFT Study as a Precursor to Synthesis

Can cyclen (1,4,7,10-tetraazacyclododecane) bind alkali metal azides? This question is addressed by studying the geometric and electronic structures of the alkali metal azide-cyclen [M(cyclen)N3] complexes using density functional theory (DFT). The effects of adding a second cyclen ring to form the sandwich alkali metal azide-cyclen [M(cyclen)2N3] complexes are also investigated. N3(-) is found to bind to a M(+) (cyclen) template to give both end-on and side-on structures. In the end-on structures, the terminal nitrogen atom of the azide group (N1) bonds to the metal as well as to a hydrogen atom of the cyclen ring through a hydrogen bond in an end-on configuration to the cyclen ring. In the side-on structures, the N3 unit is bonded (in a side-on configuration to the cyclen ring) to the metal through the terminal nitrogen atom of the azide group (N1), and through the other terminal nitrogen atom (N3) of the azide group by a hydrogen bond to a hydrogen atom of the cyclen ring. For all the alkali metals, the N3-side-on structure is lowest in energy. Addition of a second cyclen unit to [M(cyclen)N3] to form the sandwich compounds [M(cyclen)2N3] causes the bond strength between the metal and the N3 unit to decrease. It is hoped that this computational study will be a precursor to the synthesis and experimental study of these new macrocyclic compounds; structural parameters and infrared spectra were computed, which will assist future experimental work.

A theoretical study of the atmospherically important radical–radical reaction BrO + HO2; the product channel O2(a1Δg) + HOBr is formed with the highest rate

A theoretical study has been made of the BrO + HO₂ reaction, a radical–radical reaction which contributes to ozone depletion in the atmosphere via production of HOBr.

A theoretical study of the mechanism of the atmospherically relevant reaction of chlorine atoms with methyl nitrate, and calculation of the reaction rate coefficients at temperatures relevant to the troposphere

Computed rate coefficients of the atmospherically important Cl + CH₃ONO₂ → HCl + CH₂ONO₂ reaction reported for the first time.

Simulated photodetachment spectra of AlH2(-)

We have carried out high-level ab initio calculations on AlH2 and its anion, as well as Franck-Condon factor calculations, which include anharmonicity and Duschinsky rotation, to simulate the photodetachment spectrum of AlH2(-), with the aim of assigning the very recently reported photodetachment spectrum of AlH2(-) [X. Zhang, H. Wang, E. Collins, A. Lim, G. Ganteför, B. Kiran, H. Schnöckel, B. Eichhorn, and K. Bowen, J. Chem. Phys. 138, 124303 (2013)]. However, our simulated spectra do not support the assignment of the reported experimental spectrum to AlH2(-).

Determination of the photolysis rate coefficient of monochlorodimethyl sulfide (MClDMS) in the atmosphere and its implications for the enhancement of SO2 production from the DMS + Cl2 reaction

In this work, the photolysis rate coefficient of CH3SCH2Cl (MClDMS) in the lower atmosphere has been determined and has been used in a marine boundary layer (MBL) box model to determine the enhancement of SO2 production arising from the reaction DMS + Cl2. Absorption cross sections measured in the 28000-34000 cm(-1) region have been used to determine photolysis rate coefficients of MClDMS in the troposphere at 10 solar zenith angles (SZAs). These have been used to determine the lifetimes of MClDMS in the troposphere. At 0° SZA, a photolysis lifetime of 3-4 h has been obtained. The results show that the photolysis lifetime of MClDMS is significantly smaller than the lifetimes with respect to reaction with OH (≈ 4.6 days) and with Cl atoms (≈ 1.2 days). It has also been shown, using experimentally derived dissociation energies with supporting quantum-chemical calculations, that the dominant photodissocation route of MClDMS is dissociation of the C-S bond to give CH3S and CH2Cl. MBL box modeling calculations show that buildup of MClDMS at night from the Cl2 + DMS reaction leads to enhanced SO2 production during the day. The extra SO2 arises from photolysis of MClDMS to give CH3S and CH2Cl, followed by subsequent oxidation of CH3S.

Franck-Condon simulation, including anharmonicity, of the photodetachment spectrum of P2H(-): restricted-spin coupled-cluster single-double plus perturbative triple and unrestricted-spin coupled-cluster single-double plus perturbative triple -F12x potential energy functions of P2H and P2H(-)

Geometry optimization and harmonic vibrational frequency calculations have been carried out on the X̃(2)A(') state of P(2)H and the X̃(1)A(') state of P(2)H(-) using the restricted-spin coupled-cluster single-double plus perturbative triple excitation [RCCSD(T)] and explicitly correlated unrestricted-spin coupled-cluster single-double plus perturbative triple excitation [UCCSD(T)-F12x] methods. For RCCSD(T) calculations, basis sets of up to the augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality were employed, and contributions from extrapolation to the complete basis set limit and from core correlation of the P 2s(2)2p(6) electrons were also included. For UCCSD(T)-F12x calculations, different atomic orbital basis sets of triple-zeta quality with different associated complementary auxiliary basis sets and different geminal Slater exponents were used. When the P 2s(2)2p(6) core electrons were correlated in these F12x calculations, appropriate core-valence basis sets were employed. In addition, potential energy functions (PEFs) of the X̃(2)A(') state of P(2)H and the X̃(1)A(') state of P(2)H(-) were computed at different RCCSD(T) and UCCSD(T)-F12x levels, and were used in variational calculations of anharmonic vibrational wavefunctions, which were then utilized to calculate Franck-Condon factors (FCFs) between these two states, employing a method which includes allowance for anharmonicity and Duschinsky rotation. The photodetachment spectrum of P(2)H(-) was then simulated using the computed FCFs. Simulated spectra obtained using the RCCSD(T)/aug-cc-pV5Z and UCCSD(T)-F12x(x = a or b)/aug-cc-pCVTZ PEFs are compared and found to be essentially identical. Based on the computed FCFs, a more detailed assignment of the observed vibrational structure than previously reported, which includes "hot bands," has been proposed. Comparison between simulated and available experimental spectra has been made, and the currently most reliable sets of equilibrium geometrical parameters for P(2)H and its anion have been derived. The photodetachment spectrum of P(2)D, yet to be recorded, has also been simulated.

The enthalpies of formation of AsX(n) molecules, where X=H, F or Cl, and n=1, 2 or 3, by RCCSD(T) and UCCSD(T)-F12x calculations

RCCSD(T) and UCCSD(T)-F12x calculations were performed on AsX(n) molecules, where X = H, F or Cl, and n = 1, 2 or 3, and related species, in order to evaluate their enthalpies of formation (ΔH(f)(Ø)). The recommended ΔH(f)(Ø) values obtained from the present investigation are AsH, 57.7(2); AsF, -7.9(3); AsCl, 27.2(4); AsH(2), 39.8(4); AsF(2), -96.6(9); AsCl(2), -17.8(10); AsH(3), 17.1(4); AsF(3)-196.0(5) and AsCl(3), -59.1(27) kcal mole(-1). These values are anchored only on one thermodynamic quantity, namely, ΔH(f)(Ø)(As) (= 70.3 kcal mole(-1)). In the calculations, the fully-relativistic small-core effective core potential (ECP10MDF) was used for As. Contributions from outer core correlation of As 3d(10) electrons were computed explicitly in both RCCSD(T) and UCCSD(T)-F12 calculations with additional tight basis functions designed for As 3d(10) electrons. Basis sets of up to augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality were used in RCCSD(T) calculations and computed relative electronic energies were extrapolated to the complete basis set (CBS) limit. For the simplified, explicitly correlated UCCSD(T)-F12x calculations, basis sets of up to quadruple-zeta (QZ) quality were employed. Based on the RCCSD(T)/CBS benchmark values, the reliability of available theoretical and experimental values have been assessed.

Franck-Condon simulation of the photoelectron spectrum of AsCl₂ and the photodetachment spectrum of AsCl ₂⁻ employing UCCSD(T)-F12a potential energy functions: IE and EA of AsCl₂

The currently most reliable theoretical estimates of the adiabatic ionization energies (AIE(0)) from the X̃(2)B(1) state of AsCl(2) to the X̃(1)A(1) and ã(3)B(1) states of AsCl 2+, and the electron affinity (EA(0)) of AsCl(2) , including ΔZPE corrections, are calculated as 8.687(11), 11.320(23), and 1.845(12) eV, respectively (estimated uncertainties based on basis-set effects at the RCCSD(T) level). State-of-the-art ab initio calculations, which include RCCSD(T), CASSCF/MRCI, and explicitly correlated RHF/UCCSD(T)-F12x (x = a or b) calculations with basis sets of up to quintuple-zeta quality, have been carried out on the X̃(2)B(1) state of AsCl(2) , the X̃(1)A(1) , ã(3)B(1) , and Ã(1)B(1) states of AsCl 2+, and the X̃(1)A(1) state of AsCl 2-. Relativistic, core correlation and complete basis-set (CBS) effects have been considered. In addition, computed UCCSD(T)-F12a potential energy functions of relevant electronic states of AsCl(2) , AsCl (2)(+), and AsCl( 2)(-) were used to calculate Franck-Condon factors, which were then used to simulate the valence photoelectron spectrum of AsCl(2) and the photodetachment spectrum of AsCl (2)(-), both yet to be recorded. Lastly, we have also computed the AIE and EA values for NCl(2) , PCl(2) , and AsCl(2) at the G4 level and for SbCl(2) at the RCCSD(T)/CBS level. The trends in the AIE and EA values of the group V pnictogen dichlorides, PnCl(2) , where Pn = N, P, As, and Sb, were examined. The AIE and EA of PCl(2) were found to be smaller than those of AsCl(2) , contrary to the order expected from the IE values of P and As.

Photoionization of iodine atoms: angular distributions and relative partial photoionization cross-sections in the energy region 11.0-23.0 eV

Relative partial photoionization cross-sections and angular distribution parameters, beta, have been measured for the first, I(+)((3)P(2))<--I((2)P(3/2)), and fourth, I(+)((1)D(2))<--I((2)P(3/2)), (5p)(-1) photoelectron (PE) bands of atomic iodine, by performing angle-resolved constant-ionic-state (CIS) measurements on these PE bands in the photon energy range 11.0-23.0 eV. Three Rydberg series, two ns and one nd series, which converge to the I(+) (3)P(1) limit at 11.33 eV and four Rydberg series, two ns and two nd series, which converge to the I(+) (1)D(2) limit at 12.15 eV were observed in the first PE band CIS spectra. The fourth band CIS spectrum showed structure in the 12.9-14.1 eV photon energy range, which is also seen in the first band CIS spectra. This structure arises from excitation to ns and nd Rydberg states that are parts of series converging to the I(+) (1)S(0) limit we reported on earlier, as well as 5s-->5p excitations in the photon energy range 17.5-22.5 eV. These atomic iodine CIS spectra show reasonably good agreement with the equivalent spectra obtained for atomic bromine. The beta-plots for the first PE band recorded up to the I(+) (3)P(1) and I(+) (1)D(2) limits only show resonances corresponding to some of the 5p-->nd excitations observed in the first band CIS spectra scanned to the I(+) (1)D(2) limit (12.15 eV). These plots are interpreted in terms of an angular momentum transfer model with the positive values of beta obtained on resonances corresponding to parity allowed j(t)=1 and 3 channels and the off-resonance negative beta values corresponding to parity unfavored channels, where j(t) is the quantum number for angular momentum transfer between the molecule, and the ion and photoelectron. The beta-plots recorded for iodine are significantly different from those obtained for atomic bromine. Comparison of the experimental CIS spectra and beta-plots with available theoretical results highlights the need for higher level calculations which include factors such as configuration interaction in the initial and final states, relativistic effects including spin-orbit interaction, and autoionization via resonant Rydberg states.

High-level ab initio calculations on HGeCl and the equilibrium geometry of the A1A'' state derived from Franck-Condon analysis of the single-vibronic-level emission spectra of HGeCl and DGeCl

CCSD(T) and/or CASSCF/MRCI calculations have been carried out on the X(1)A' and A(1)A'' states of HGeCl. The fully relativistic effective core potential, ECP10MDF, and associated standard valence basis sets of up to the aug-cc-pV5Z quality were employed for Ge. Contributions from core correlation and extrapolation to the complete basis set limit were included in determining the computed equilibrium geometrical parameters and relative electronic energy of these two states of HGeCl. Based on the currently, most systematic CCSD(T) calculations performed in this study, the best theoretical geometrical parameters of the X(1)A' state are r(e)(HGe) = 1.580 +/- 0.001 A, theta(e) = 93.88 +/- 0.01 degrees and r(e)(GeCl) = 2.170 +/- 0.001 A. In addition, Franck-Condon factors including allowance for anharmonicity and Duschinsky rotation between these two states of HGeCl and DGeCl were calculated employing CCSD(T) and CASSCF/MRCI potential energy functions, and were used to simulate A(1)A'' --> X(1)A' SVL emission spectra of HGeCl and DGeCl. The iterative Franck-Condon analysis (IFCA) procedure was carried out to determine the equilibrium geometrical parameters of the A(1)A'' state of HGeCl by matching the simulated, and available experimental SVL emission spectra of HGeCl and DGeCl of Tackett et al., J Chem Phys 2006, 124, 124320, using the available, estimated experimental equilibrium (r(e)(z)) structure for the X(1)A' state, while varying the equilibrium geometrical parameters of the A(1)A'' state systematically. Employing the derived IFCA geometry of r(e)(HGe) = 1.590 A, r(e)(GeCl) = 2.155 A and theta(e)(HGeCl) = 112.7 degrees for the A(1)A'' state of HGeCl in the spectral simulation, the simulated absorption and SVL emission spectra of HGeCl and DGeCl agree very well with the available experimental LIF and SVL emission spectra, respectively.