Photoelectron Spectroscopy Study of Quinonimides - ACS Publications

Jul 21, 2017 - The Department of Chemistry Purdue University West Lafayette, Indiana ... University of Southern California Los Angeles, 90089, United ...
1 downloads 0 Views 975KB Size
Subscriber access provided by UNIV OF NEWCASTLE

Article

Photoelectron Spectroscopy Study of Quinonimides Ekram Hossain, Shihu M. Deng, Samer Gozem, Anna I. Krylov, Xue-Bin Wang, and Paul G Wenthold J. Am. Chem. Soc., Just Accepted Manuscript • DOI: 10.1021/jacs.7b05197 • Publication Date (Web): 21 Jul 2017 Downloaded from http://pubs.acs.org on July 22, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of the American Chemical Society is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Photoelectron Spectroscopy Study of Quinonimides By

Ekram Hossain,† Shihu M. Deng,‡ Samer Gozem,§ Anna I. Krylov,§ Xue-Bin Wang‡* and Paul G. Wenthold†* Contribution from

The Department of Chemistry Purdue University West Lafayette, IN 47906 and Physical Sciences Division, Pacific Northwest National Laboratory P. O. Box 999, MS k8-88 Richland, WA 99352, USA and Department of Chemistry University of Southern California Los Angeles, 90089, USA Submitted to The Journal of the American Chemical Society May 2017 † ‡ §

Purdue University PNNL USC

1 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Structures and energetics of o-, m- and p-quinonimide anions (OC6H4N−) and quinoniminyl radicals have been investigated by using negative ion photoelectron spectroscopy. Modeling of the photoelectron spectrum of the ortho isomer shows that the ground state of the anion is a triplet, while the quinoniminyl radical has a doublet ground state with a doubletquartet splitting of 35.5 kcal/mol. The para radical has doublet ground state, but a band for a quartet state is missing from the photoelectron spectrum indicating that the anion has a singlet ground state, in contrast to previously reported calculations. The theoretical modeling is revisited here, and it is shown that accurate predictions for the electronic structure of the para quinonimide anion require both an accurate account of electron correlation and a sufficiently diffuse basis set. Electron affinities of o- and p-quinoniminyl radicals are measured to be 1.715 ± 0.010 and 1.675 ± 0.010 eV, respectively. The photoelectron spectrum of the m-quinonimide anion shows that the ion undergoes several different rearrangements, including a rearrangement to the energetically favorable para isomer. Such rearrangements preclude a meaningful analysis of the experimental spectrum.

2 ACS Paragon Plus Environment

Page 2 of 38

Page 3 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Introduction Benzoquinones (1o, 1m, 1p) are important chemical structures that are found in a large variety of systems with different applications. As convenient electron acceptors, quinones play

essential roles in critical biological processes including photosynthesis and respiration.1-11 However, quinones are also involved in several harmful processes in biological systems. Orthoand para-quinones can be easily formed by oxidation of the corresponding catechols, and can even be formed from the oxidation of aromatic molecules, such as benzene. Ortho-quinones formed in this manner have been proposed as the sources of benzene toxicity12 and carcinogenicity.13 Similarly, ortho-quinones formed in the oxidation of urushiol cause the contact-dermatitis associated with plants like poison ivy.14 Benzoquinones are also utilized in various chemical applications, serving as dyestuffs,15-17 battery electrodes18-20 and hydrogen storage materials,21 and as oxidizing agents22-24 and dienophiles25 in organic synthesis. The electronic structures of benzoquinones are highly dependent on their topology. The stable quinones encountered in the situations described above are the ortho- and para-isomers, and the stabilities are readily apparent in their Kekulé valence structures, as shown in 1o and 1p. Consistent with this view, both 1o and 1p have closed-shell singlet ground states,26,27 with triplet-states lying 1.68 and 2.32 eV higher in energy, respectively.28 In contrast, the metaisomer, 1m, has a non-Kekulé structure and an open-shell, triplet ground-state28-31 9.0 kcal/mol lower in energy than the lowest-energy singlet state.30 The relative stabilities of the ortho- and para-isomers compared to 1m are reflected in the absolute enthalpies of formation. The 3 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

measured enthalpies of formation of 1o and 1p are -21.0 ± 3.1 and -29.3 ± 0.9 kcal/mol, respectively, whereas ∆Hf,298(1m) = 8.1 ± 3.4 kcal/mol.29,32 The higher enthalpy of formation of 1o compared to 1p has been attributed to electron-pair repulsion in the ortho-isomer that is not present in the para-isomer.29 The benzoquinones 1o – 1p are just one type within the class of molecules with quinoidal structures. As shown in Figure 1, the substituents on the aromatic ring can be varied to create quinone analogs. The most common of these would include the quinomethanes (X = O, Y = CR2, Figure 1)33-39 or the quinodimethanes (“xylylenes,”)40-48 where both X and Y are carboncentered substituents. However, essentially all multivalent atoms can be incorporated, and in almost any combination. Whereas each variation will have unique properties, the common

Figure 1. Some variations of quinoidal molecules. theme is that the ortho- and para-isomers will have stable,33,49 Kekulé structures with closedshell singlet electronic states, whereas the meta-isomers are higher-energy,50 non-Kekulé structures with high-spin electron coupling in the ground state, which is why meta-phenylene derivatives are commonly used as the basis for high-spin materials.51-54 We recently reported a flowing afterglow study of quinoidal imides, including quinonimides (2, X = O, Y = N−) and quinomethanimides (3, X = CH2, N−).55

4 ACS Paragon Plus Environment

Page 4 of 38

Page 5 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Significant differences were observed in the reactions of 2p and 3p with nitric oxide. Whereas 3p reacts with NO merely by addition, 2p reacts by addition and by a novel, N-O exchange reaction to form the semiquinone radical anion, as shown in eq 1.

The observed differences in reactivity were attributed to differences in the electronic structures. The addition reaction that occurs with 3p is consistent with what is expected for the closed-shell, imide anion.56 In contrast, the N-O exchange reaction has been observed previously with nitrene radical anions57 and is interpreted as indicating an open-shell electronic structure.58 Consequently, our previous study indicates that while 3p has a robust ground-state singlet, 2p is either a ground-state triplet or a singlet with a thermally accessible triplet state.55 A low-lying triplet state for 2p is somewhat surprising considering that it is isoelectronic with p-benzoquinone, 1p, which is a singlet with a large (> 2 eV) singlet-triplet splitting.28 However, it can be rationalized by considering that 2p can alternatively be viewed as an oxosubstituted phenylnitrene, as shown in Scheme 1. Considering that the energies of closed-shell singlet states in aromatic nitrenes are approximately 30 kcal/mol higher than the triplets,59,60 the fact that the singlet is close in energy to the triplet indicates that the singlet quinonimide is highly 5 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

stabilized, at least for a closed-shell singlet nitrene. Triplet nitrene 2p has also been proposed as an intermediate in the photolysis of p-azidophenoxide in aprotic conditions,61 affirming that the triplet state is chemically accessible.

The conclusions about the electronic structures of 2p and 3p were consistent with reported electronic structure calculations55 that predict 3p to be a ground-state singlet with a singlet-triplet splitting of about 5 kcal/mol, with 2p a ground state triplet, with a singlet that is 1 – 2 kcal/mol higher in energy. In this work, we describe a negative ion photoelectron spectroscopy (NIPES) study of the quinonimides, 2o – p. Negative ion photoelectron spectroscopy provides a means to investigate the electronic properties of ions as well as those of the neutral molecules obtained upon electron detachment.62 Therefore, in addition to being a study of the quinonimides, this study also provides new spectroscopic information for quinoniminyl radicals, 4.

The quinoniminyl radicals are isoelectronic with the benzoquinone radical cations.63,64 On the basis of isotropic proton coupling in ESR experiments, Chandra et al.63 proposed that ionization of 1p occurs from a σ orbital. By anology, 4p would be a σ-radical, with a closed-shell π system, leading to a doublet ground state, and the ortho-radical (4o) should be similar. High6 ACS Paragon Plus Environment

Page 6 of 38

Page 7 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

spin coupling of the meta-phenylene arrangement in 4m, however, should favor a quartet ground-state.65 The para-quinoniminyl radical, 4p, has previously been suggested as a possible intermediate in the Gibbs reaction between N-haloquinonimines and phenoxides.66 In this study, we use NIPES to examine the ortho-, meta- and para-isomers of the quinonimides. In addition to determining electron affinities of the quinoniminyl radicals, we also use modeling of the spectroscopic features to show that ortho-quinonimide, 2o, has a triplet ground state, consistent with electronic structure calculations.

For the para-isomer, 2p, the

spectroscopy results are best interpreted as indicating a singlet ground state, in contrast to previous electronic structure predictions, which have been revisited in this work. The spectra and calculations presented for 2p, combined with the reactivity studies reported previously,55 suggest that whereas the isomer has a singlet ground state, it also has a low-lying triplet state that can be accessed during reaction.

Experimental Methods All NIPES experiments were carried out at PNNL using an apparatus consisting of an electrospray ionization source (ESI), a cryogenic ion trap, and a magnetic-bottle time-of-flight (TOF) photoelectron spectrometer.67

Azidophenols were synthesized from o-, m- and p-

aminophenols by using procedures similar to those described previously,68,69 and used without purification. A 5×10-4 M acetonitrile solution of each o-, m- and p-azidophenol (N3-C6H4-OH) titrated dropwise with 0.01 M aqueous KOH was prepared in the glove box under an N2 atmosphere. The electrospray ionization conditions were optimized to afford maximum intensity for the dissociation (M –H −N2) product anions, i.e. 2o, 2m, and 2p, m/z 106. The anions generated by ESI were guided by quadrupole ion guides into an ion trap, where they were

7 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

accumulated and cooled for 20-100 ms by collisions with cold buffer gas at 20 K, before being transferred into the extraction zone of a TOF mass spectrometer. The cooling of the anions to 20 K eliminated the possibility of the appearance of hot bands in the NIPES spectra. Each of 2o, 2m, or 2p anions were then mass selected, and maximally decelerated before being photodetached with 266 (4.661 eV), 355 (3.496 eV) and 532 nm (2.331 eV) photons from a Nd:YAG laser. The laser was operated at a 20 Hz repetition rate, with the ion beam off at alternating laser shots, in order to enable shot-to-shot background subtraction to be carried out. Photoelectrons were collected with ca. 100% efficiency with the magnetic bottle and analyzed in a 5.2 m long electron flight tube. The recorded TOF photoelectron spectrum was converted into an electron kinetic energy spectrum by calibration with the known NIPE spectra of I– and OsCl2–. The electron binding energy spectrum was obtained by subtracting the electron kinetic energy spectrum from the energy of the detaching photons. The energy resolution was about 2%, i.e., ∼20 meV for 1 eV kinetic energy electrons.

Computational Methods All geometries and frequencies reported in this work are computed with density functional theory (DFT) with the B5050LYP functional70,71 and 6-311++G(d,p) basis set. B5050LYP is a hybrid functional that uses 50% Hartree-Fock, 8% Slater, and 42% Becke exchange and 19% VWN+81% LYP correlation, and is used because extensive benchmarks indicate that it yields improved geometries, energies, and vibrational frequencies in combination with the spin-flip (SF)-DFT approach.71 We employ SF-DFT for low-spin (singlet or doublet) states in ortho and para isomers (2o, 2p, 4o, 4p), but not for meta (2m, 4m) because in that case SF states were strongly spin-contaminated. For high-spin (triplet and quartet) states and for all 8 ACS Paragon Plus Environment

Page 8 of 38

Page 9 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

meta isomer calculations, standard DFT was used. The DFT or SF-DFT geometries were then used for subsequent energy calculations with coupled-cluster-based methods. Singlet-triplet and doublet-quartet energy splittings were computed with the spin-flip equation-of-motion coupled-cluster (EOM-SF-CC) method.72,73 The SF approach employs a high-spin reference state and uses a spin-flipping operator that allows description of both the high-spin and low-spin target states on an equal footing.72,74-76 As needed for the systems presented in this study, the SF approach is capable of describing multi-configurational character encountered in open-shell singlets and doublet radicals, which may not be described correctly in a single-configuration formalism. Previous studies have shown that the SF variant of EOMCC72,73,77 gives accurate singlet-triplet splittings in aromatic nitrenes. EOM-SF-CC energies are calculated with singles and doubles excitations, denoted EOMSF-CCSD. In addition, the effect of triple excitations was included by perturbative triples correction, denoted EOM-SF-CCSD(dT).77,78 Coupled-cluster calculations were carried out using B3LYP orbitals as the reference orbitals, as described previously,60,72 to mitigate spincontamination.79 Electron attachment energies were computed between the triplet and quartet states by taking the difference in coupled-cluster energies with single and double excitations and perturbative triples correction, CCSD(T). Electron attachment energies for other states are then derived on the basis of triplet-quartet CCSD(T) energy difference and EOM-SF-CCSD(dT) splittings, as described above.80 All reported energies are adiabatic (i.e., accounting for geometry relaxation of different states) and are corrected for zero-point energies (ZPE) computed from DFT or SF-DFT unscaled vibrational frequencies at the corresponding geometry. For all coupled-cluster calculations, we employ Dunning’s correlation consistent basis sets.81

9 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 38

Unless stated otherwise, the values for EOM-SF-CCSD(T)/(aug-)cc-pVTZ calculations were estimated by using the additive approach shown in eq 2.

EPROJECTED=E(CCSD/(aug-)cc-pVTZ)-E(CCSD/(aug-)cc-pVDZ)

(2)

+E(CCSD(T)/(aug-)cc-pVDZ)

Previous studies,60 have found that this approach gives singlet-triplet energy differences within ca. 0.5 kcal/mol of those obtained by a direct EOM-SF-CCSD(T)/cc-pVTZ calculation. In this work, the projected energies will be indicated as p-EOM-SF-CCSD(T)/(aug-)cc-pVTZ. All the calculations are performed using QCHEM ab initio package.82,83 Photoelectron spectra, including Frank-Condon vibrational progressions and detachment cross-sections, have been calculated using procedures described previously.

Frank-Condon

factors are calculated from DFT or SF-DFT normal modes and frequencies (depending on the isomer and state; see above). The Frank-Condon factors for detachment are calculated using ezSpectrum.84 Detachment cross sections were calculated using ezDyson.85,86 Photoionization and photodetachment cross-sections can formally be computed from the photoelectron dipole matrix element, DkIF, eq 3, where ΨN and ΨN-1 are the wavefunctions of the initial (N-electron) and final (N-1-electron) states, respectively, ru is the dipole operator for the ionizing radiation, and Ψel is the wave function of the ejected electron.87,88

 = 〈Ψ | |   ∙ Ψ 〉

(3)

10 ACS Paragon Plus Environment

Page 11 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Owing to electron indistinguishability, DkIF can be expressed in terms of a Dyson orbital, ϕd, such that it is now calculated from one-electron functions, as shown in eq 4.89,90

 = 〈  | |Ψ 〉

(4)

The Dyson orbitals are correlated extensions of Koopmans’ theorem. Here Dyson orbitals are calculated by using the ionization potential EOM-CCSD (EOM-IP-CCSD) approach,85,90 which accounts for orbital relaxation and electron correlation. The wave function for reference triplet state of the anion was computed at the CCSD/aug-cc-pVDZ level of theory, and the neutral state wave functions were calculated by using the EOM-IP-CCSD/aug-cc-pVDZ method. To analyze diradical character, we also compute the Head-Gordon’s index, using the EOM-SF-CCSD/aug-cc-pVDZ wave functions. The Head-Gordon index is computed from state density matrices (or, more precisely, populations of natural orbitals) and therefore does not depend on the choice of molecular orbitals. In this work, we use the non-linear variant (eq 5),91,92

!

 , =   2 − 



(5)

"

where M is the total number of orbitals, and ni is the occupation number for each natural orbital. The value of Nu,nl is a measure of the effective number of unpaired electrons, such that for a closed shell system it is 0, a radical gives 1, and a perfect biradical gives 2, etc.

11 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Computational Results Before reporting the results of the photoelectron measurements, we provide a discussion of the electronic structures of the ions and corresponding neutrals. The electronic structures of the anions, 2o, 2m and 2p, can be understood from the framework of substituted phenyl nitrenes. The non-bonding molecular orbitals in nitrenes usually include the benzylic-like π orbital and the σ (p) orbital localized on nitrogen (Scheme 2).

An O– substituent adds an additional benzylic-like π orbital that can interact with the π orbital of the nitrene. For the O– in the ortho and para positions, the two π orbitals have significantly different energies, resulting from quinoidal bonding and anti-bonding combinations as shown at the top of Figure 2. Consequently, the bonding π orbitals are expected to be significantly stabilized with respect to the σ orbitals, while the anti-bonding orbitals are destabilized. For the meta isomer, the combinations of the two benzylic orbitals would be nonbonding molecular orbitals, as shown at the bottom of Figure 2, resembling the orbitals in other meta-phenylene-type systems,

12 ACS Paragon Plus Environment

Page 12 of 38

Page 13 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure 2. Symmetry-adapted π molecular orbitals in quinonimides

such as 1,3-quinodimethane (m-xylylene). However, because they lack the bonding interactions present in the orbitals for 2o and 2p, there is no benefit to the combinations for 2m, such that the π orbitals in that isomer are generally the localized benzylic orbitals, as in Scheme 2. Anions 2o, 2m and 2p are created by adding 4 electrons to the three orbitals (σ and two π orbitals). For 2o and 2p, the two lowest energy electrons occupy the bonding π orbital, and the other two electrons can occupy either the σ or π* orbitals. In our previous study, 2p was calculated to be a ground state triplet, with the σ2 singlet state only 1 – 2 kcal/mol higher in energy.55 Ion 2o is calculated to be a more robust ground-state triplet, with a significantly larger singlet-triplet splitting of 8.7 kcal/mol at the p-EOM-SF-CCSD(T)/aug-cc-pVTZ level of theory. The increased singlet-triplet energy gap is likely the result of destabilization of the σ2 state in the ortho isomer resulting presumably from repulsion between σ electron pairs on the oxygen and nitrogen and, possibly, from 1,4-interaction in the π orbital. In fact, the destabilization is so

13 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

significant that it causes a non-planar geometry for the singlet state when using a spin-restrictedsingle determinant wave function, such as RHF, B3LYP or BLYP. With unrestricted and spinflip methods, however, the state is predicted to be planar. This phenomenon is a manifestation of the pseudo Jahn-Teller effect and is often encountered in open-shell species.93-95 The predicted electronic structure of the meta anion isomer is very different from those of the ortho and para isomers. Because the orbitals are all very similar in energy, there are multiple options for orbital occupancy. However, the electronic state with the lowest energy is predicted to have two electrons in a phenoxide benzylic orbital, with the other two electrons occupying the σ orbital and the π benzylic orbital on nitrogen. Therefore, in 2m the π orbitals localize into their benzylic forms (linear combinations of π1 and π2), presumably because it puts the negative charge predominantly on the oxygen. Consequently, 2m is best described as a phenoxide anion orthogonal to a phenyl nitrene. Finally, because it is a phenyl nitrene distinct from the phenoxide, it is not surprising that it has a triplet ground state. Similarly, the lowest energy singlet state of 2m is likely the singlet nitrene. The singlet-triplet energy gap computed with EOM-SF-CCSD/aug-cc-pVTZ is 17.3 kcal/mol. The difference in the electronic structure of 2m compared to those of 2o and 2p is readily observed in Head-Gordon’s indices for the lowest-lying singlet state of each isomer. While Nu,nl for both 2o and 2p is 0.04, indicating closed shell character consistent with the π2σ2 assignment discussed above, the same index is 2.04 for 2m, indicating that even the lowest energy singlet state has significant diradical character. The neutral quinoniminyl radicals, 4, formed upon detachment of 2 are nominally triradicals, in that they have 3 electrons in three nearly degenerate orbitals. However, because of the strong bonding interactions that stabilize the π orbitals in 4o and 4p, these isomers have

14 ACS Paragon Plus Environment

Page 14 of 38

Page 15 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

doublet ground states, with high-energy quartet states (Table 1). In contrast, the ground state of 4m is predicted to be a quartet, with multiple low-lying doublet states. The relevant computed state energies are shown in Table 1, including computed electron affinities and doublet-quartet splittings.

Experimental Results In this section, we report the measured negative ion photoelectron spectra for the [M – N2 – H]− ions obtained upon electrospray ionization of basic solutions of o-, m- and p-azidophenol. Survey spectra, measured at 266 nm, are shown in Figure 3. Detailed descriptions of each spectrum are provided in the sections below.

ortho (2o) The 266 nm spectrum for the ortho isomer, Figure 3a, shows three identifiable features. The lowest energy feature is a narrow band with an origin at 1.715 ±0.010 eV. The highest energy feature is dominated by the origin peak at 3.256 ± 0.005 eV, indicating a nearly vertical transition. The center band begins at 2.54 eV. However, the relative intensity of this feature is variable and affected by source conditions. Consequently, the center band results from detachment of a different ion from that giving the other features in the spectrum. The 532 nm photoelectron spectrum of 2o (vide infra) only contains the lowest energy feature. However, because the photodetached electrons have lower energies, they are better resolved than in the 266 nm spectrum. The identifiable peaks for the ortho isomer are listed at the top of Table 1.

15 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. 266 nm photoelectron spectra of the m/z 106 ions derived from a) ortho-, b) meta- and c) para-azidophenol. The * in (a) and (b) indicate bands whose relative intensity depends on source conditions.

para (2p) The 266 nm spectrum of the para-substituted ion, shown in Figure 3c, has only a single intense feature, with an origin at 1.675 ± 0.010 eV. As with the ortho isomer, the 532 nm spectrum (vide infra) shows a well-resolved band. The observed peaks are shown at the bottom of Table 1.

meta (2m) The 266 nm spectrum of 2m, shown in Figure 3b, has three features. However, the lowest energy band is identical to the lowest energy feature for the para isomer, in Figure 3c,

16 ACS Paragon Plus Environment

Page 16 of 38

Page 17 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Table 1. Observed peaks in the photoelectron spectra of 2o, 2m and 2p Ion and neutral Vibrational electronic state Peak eBE, eV Relative Energy, cm-1 Assignment 2o doublet

1.715 ± 0.010 1.780 1.845 1.865 1.905 1.985

impurity

2.545 2.65

quartet

3.256 ± 0.005 3.310 3.430

2m X

A 2p doublet

Calculated Energya

EA 525 1050 1210 1530 2160

origin ν8 2ν8 ν20 ν23 ν8+ν23

1.61 eVb 532

+35.5 kcal/mol 440 1410

origin ν5 ν24 and ν25

+36.0 kcal/molc 429 1386, 1428

EA 445 845 1250 1655 2100

origin ν5 ν11, 2ν5 3ν5, ν5+ν11 2ν5+ν11 ν5+2ν11

1.47 eVb 442 721, 883

1186 1508

2.11 ± 0.02 2.22 2.35 2.48 3.15 (vertical)

1.675 ± 0.010 1.730 1.780 1.830 1.880 1.935

+50.2 kcal/molc

quartet a

Vibrational energies, in cm-1, calculated at the (SF)-B50H50/6-311++G(d,p) level of theory, scaled by 0.92; electron affinities and doublet-quartet splittings are adiabatic energy differences, and include ZPE correction b Calculated from p-CCSD(T)/aug-cc-pVTZ triplet-quartet energy differences and p-EOM-SFCCSD(dT)/aug-cc-pVTZ splittings. c p-EOM-SF-CCSD(dT)/aug-cc-pVTZ 17 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

indicating that the meta isomer undergoes rearrangement prior to photodetachment. The most intense feature in the spectrum of 2m has an origin at 2.11 eV, with many closely-spaced peaks. The identifiable peaks are listed in Table 1. A third, weaker feature is also observed in the spectrum. Although it is unresolved, it has a maximum at a vertical detachment energy of about 3.15 eV.

Discussion In this section, we present an analysis of the measured photoelectron spectra of the ions derived from ortho-, meta- and para-azidophenoxide, with a focus on determining the electronic states of the anions and the radicals.

ortho The spectrum of the ortho isomer contains features attributed to two different ions. The high and low energy features result from a common ion, and are readily assigned to formation of the doublet and quartet states of 4o, respectively, formed upon detachment of 2o. The measured electron affinity, 1.715 eV, agrees reasonably well with the calculated EA for formation of the triplet anion, 1.61 eV (Table 1). More importantly, the energy difference between the bands, 35.5 kcal/mol, is in very good agreement with the predicted doublet-quartet splitting of 36.0 kcal/mol, obtained at the p-EOM-SF-CCSD(T)/aug-cc-pVTZ level of theory (Table 1). The photoelectron spectrum of 2o shown in Figure 3a also clearly demonstrates that the ion has a triplet electronic state. The most important evidence for the triplet state of 2o is the fact that detachment to both the doublet and quartet states are observed. Given that the selection rule

18 ACS Paragon Plus Environment

Page 18 of 38

Page 19 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

for photodetachment is ∆S = ± 1, the only way to form both the doublet and quartet states of 4o is to detach from the triplet state of 2o. To confirm that the experimental peaks indeed are assigned to doublet and quartet states formed from the triplet anion, we refer to quantum chemical calculations. In this way, we have unequivocal evidence that the ground state is a triplet. In particular, structural differences between the triplet and singlet states are manifested in the shapes of the features in the spectrum. The 532 nm photoelectron spectrum of the ion 2o, shown in Figure 4, has a well-resolved vibrational progression for the doublet state. For comparison, we have

Figure 4. Comparison of the low-energy feature in the 532 nm photoelectron spectrum of 2o with the band predicted from (SF)-DFT B5050LYP/6-311++G(d,p) geometries and frequencies of the (a) triplet state and (b) singlet state of the anion and the doublet state of the radical, 4o. Calculated harmonic frequencies were scaled by 0.92.

calculated the expected Frank-Condon profile for formation of the doublet from the triplet (Figure 4a) and singlet (Figure 4b) states of the ion, using geometries and frequencies calculated at the DFT or SF-DFT B5050LYP levels of theory.84 The calculated spectrum obtained when

19 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

using the triplet anion is in excellent agreement with the measured spectrum, in both the positions of the peaks and in their intensities. In contrast, the calculated band for detachment from the singlet state has poor resemblance to the measured spectrum. Therefore, the FrankCondon calculation further confirms that the ion has a triplet ground state. Using this approach, we have also calculated the shape of the quartet feature, shown in Figure 5. As with the doublet feature, there is outstanding agreement between the predicted relative peak positions and intensities for detachment from the triplet anion, and those measured experimentally. The band for the doublet state in Figure 5 is the same as that shown in Figure 4a, but with significantly broader peaks (450 cm-1 vs 150 cm-1 full-width at half-maximum) due to the lower resolution of the band in the 266 nm spectrum. Overall, these results create a consistent picture for the photoelectron spectrum, indicating detachment of the triplet state of the anion forming the doublet and quartet states of the neutral.

Figure 5. Comparison of measured (266 nm) and calculated spectra for 2o. The calculated spectra correspond to detachment of the triplet state of 2o. Calculated harmonic frequencies were scaled by 0.92.

The middle band observed in the spectrum has an identifiable origin peak at 2.54 eV, but the band has a relative intensity that is dependent on source conditions, and therefore results 20 ACS Paragon Plus Environment

Page 20 of 38

Page 21 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

from detachment of an isobaric ion. Given the azidophenoxide precursor, the most likely ionic formula is either C6H4NO– (same as 2o) or C5H4N3–. Either would result from more complex dissociation of the precursor. We have not been able to determine the actual structure of the ion leading to the observed signal, but the electron binding energy suggests an oxygen-based anion, such as those shown in Figure 6. Cyanocyclopentadiene has been previously found to be formed by rearrangement of phenylnitrene in the gas-phase.96 Simple B3LYP calculations predict cyclopentadienyl structures such as those in Figure 6 to be up to 2 eV more lower in energy than 2o.

Figure 6. Possible structures of the second ion in the spectrum of 2o

meta The photoelectron spectrum of the ion obtained from the meta isomer, shown in Figure 3b, contains a feature that can be identified as coming from detachment of the para anion, 2p, which indicates that the ion is capable of rearrangement. At the p-CCSD(T)/aug-cc-pVTZ level of theory, 2p is calculated to be 9.4 kcal/mol lower in energy than the meta-quinonimide, which is the driving force for the rearrangement.

Unfortunately, the fact that the ion rearranges makes

it difficult to analyze the rest of the spectrum. The binding energy of 2.03 eV for the most intense band is significantly lower than the predicted value of 2.52 eV for the quartet state of 4m. Moreover, the transition from 2m to the quartet state of 4m is predicted to be close to vertical, but the observed band is much broader than expected for that transition. Therefore, it does not appear that the signal can be attributed to detachment of the triplet quinonimide, 2m. 21 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Interestingly, the computed detachment to the doublet state is at 3.16 eV, which is in very good agreement with the weak feature appearing at around 3.15 eV in the experimental spectrum. This may be an indication that some features from 2m appear in the spectrum, but are overshadowed by peaks from rearranged products. As with the second band in the spectrum of 2o, the electron binding energy of the main feature in the spectrum of the meta-isomer suggests an alkoxide, and there are rational possibilities. For example, unimolecular rearrangement of aromatic nitrenes is known to result in formation of cyclic ketenimines (eq 6a).97 A similar rearrangement (eq 6b) of 2m would lead to structures such as 5a and 5b. The computed electron binding energies of 5a and 5b are 2.06 and 1.71 eV, respectively, at the B3LYP/6-31+G* level of theory, and so are in the range of the measured band. However, neither structure alone accounts for the observed feature, as the bands for 5a or 5b are predicted to be too narrow to account for the observed spectrum. Therefore, the low-energy feature in the spectrum of the meta isomer is either a mixture of ions such as 5a and 5b, or due to some other undetermined rearrangement product. Although the ketenimine structures 5a and 5b are calculated to be slightly higher in energy than triplet 2m, they lower in energy than the singlet state.

22 ACS Paragon Plus Environment

Page 22 of 38

Page 23 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

para The spectrum of the ion obtained from p-azidophenoxide is the simplest of the three, containing only one identifiable band, with an origin peak at 1.67 ± 0.01 eV. The measured electron affinity is in reasonable agreement with the predicted value for the electron affinity of the doublet state (Table 1). Although there is broad signal at higher binding energy, there is no feature that can be clearly identified as resulting from detachment to form the quartet state, which is predicted to be ~2.2 eV higher in energy than the doublet (Table 1), and therefore within the energy range of the experiment. The fact that there is only a single product formed upon detachment of the para anion is significant. As noted in the discussion of the ortho spectrum, detachment of the triplet state of the anion would form both the doublet and quartet state of the neutral product. Therefore, a possible explanation for the absence of a peak for the quartet state is that the ion is in a singlet state. In order to confirm the singlet state, we have carried out an analysis of the Franck-Condon progression similar to what we did for 2o. Figure 7 shows a comparison between the experimental 532 nm photoelectron spectrum of 2p and those predicted for formation of the

23 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

doublet upon detachment of the triplet (Figure 7a) and singlet (Figure 7b) states of the anion.

Figure 7. Comparison of the low-energy feature in the 532 nm photoelectron spectrum of 2p with the band predicted from (SF)-DFT B5050LYP/6-311++G(d,p) geometries and frequencies of the (a) triplet state and (b) singlet state of the anion and the doublet state of the radical, 4p. Calculated harmonic frequencies were scaled by a factor of 0.92. Although the agreement between the calculated and predicted spectra is clearly better for the singlet state, there is not a sufficient difference in the quality of the agreement between the spectra for the singlet and triplet states for the para isomer to make an unequivocal assignment of the electronic state, and the differences observed for the triplet state in Figure 7a are not sufficient to rule it out as the electronic state of the anion. Because the Frank-Condon calculation cannot sufficiently distinguish between the triplet and singlet states of 2p, the assignment of the electronic state depends on the absence of the quartet band. However, this conclusion requires the assumption that the quartet band should be observable to a sufficient extent. Although the spectrum for 2o, for example, suggests that the quartet band should be strong, it is not conclusive. Therefore, we have calculated the expected

24 ACS Paragon Plus Environment

Page 24 of 38

Page 25 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

relative cross sections for formation of the doublet and quartet states using the Dyson orbitals computed for correlated EOM-CCSD wave functions, as described above.85,86 The challenge in this calculation involves the choice for the wave function for the detached electron. For simple atomic systems, the detached electron can be described as a Coulomb wave in the case of photoionization,98,99 or as a plane wave in the detachment of an atomic anion.88 To a first approximation, one can treat the detachment of 2p as similar to that of an atomic anion, and describe it as a plane wave. However, for the ortho isomer, it was found that the relative intensity of the quartet feature was underestimated by using a plane wave approximation (see Supporting Information). Alternatively, it has been found previously that, for larger molecular systems, using a Coulomb wave with a partial charge between 0 and 1 often gives the best agreement with experiment.99,100 Because we do not have a rigorous procedure to determine the best effective charge, we have computed the cross section at Z values of 0.0, 0.2, 0.4, 0.6, 0.8 and 1.0, and used the average values to obtain relative cross sections. This approach has been used successfully to predict the cross sections for ionization of aldehydes and enols,100 and is the approach used to determine the relative intensities for the 2o shown in Figure 5. For the quinonimides, the relative intensities of the quartet states obtained by using a modified Coulomb wave are consistently larger than those obtained by using a plane wave approximation, regardless of the chosen Z value (see Supporting Information). The average calculated total absolute cross sections for detachment of triplet 2p to form the doublet and quartet states of 4p are 4.7 and 6.2 megabarns, respectively, giving relative cross sections of about 1/1.3 for formation of the doublet and quartet, respectively. Therefore, as expected, detachment of the triplet anion should lead to abundant formation of the quartet state. Scaling the integrated Frank-Condon profiles by the calculated total cross sections results in the

25 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

predicted spectra shown in Figure 8, which confirm that the quartet is not being formed as expected if 2p were a triplet anion. Although the predicted relative cross-section for formation of the quartet state is lower when using a plane wave approximation for the detached electron, the band for the quartet state is nonetheless predicted to be prominent in the spectrum (see

Figure 8. Comparison of measured (266 nm) and predicted spectra for 2p. The calculated spectra correspond to detachment of the triplet state of 2p.

Supporting Information). Consequently, the photoelectron spectrum is best interpreted as resulting from detachment of the singlet state of 2p.

Vibrational assignments The calculation of the Frank-Condon intensities from theoretical structures and frequencies can be used to identify the vibrations active in the photoelectron spectra. For example, as shown in Table 1, the active modes in the doublet feature in the spectrum of 4o include ν8 (525 cm-1), ν20 (1210 cm-1) and ν23 (1530 cm-1). The vibrational motions associated

26 ACS Paragon Plus Environment

Page 26 of 38

Page 27 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

with these modes are shown in Scheme 3. The assigned modes are those predicted to be active

in the Frank-Condon calculations. The mode ν8 is a ring-deformation mode, similar to the ν6a mode in benzene, and is commonly active in negative ion photoelectron spectra of aromatic molecules.101-104 The mode ν20 is predominately stretching of the CO and CN carbon atoms, whereas ν23 is the symmetric combination of the CO and CN stretching motions.

In contrast,

the active modes in the quartet state of 4o, ν24 and ν25, are anti-symmetric combinations of the CO and CN stretching motions. The active ring-deformation mode, ν5, involves the CN stretch. The active modes in the doublet state of 4p consist of the ring-deformation mode, ν5, and a ring breathing mode, ν11.

27 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Comparison to reactivity and theoretical results As noted in the introduction, the reactivity of 2p has been previously interpreted as indicating that 2p is either a ground-state triplet or a singlet with a thermally accessible triplet state.55 The 20K photoelectron spectrum indicates that the ion has a ground state singlet, and so if reactivity is due to the triplet state, it must be accessed during the reaction. Kenttämaa and coworkers have found that singlet diradicals are capable of undergoing reactions involving the triplet state,105,106 and therefore the presence of a singlet ground state for 2p is not inconsistent with the apparent radical reactivity. The singlet ground state found for 2p is, however, inconsistent with predictions from electronic structure theory reported previously.55 As noted above, all three quinonimide anions were predicted to have triplet ground states. The calculations leading to the prediction of the triplet ground state for 2p reported in the previous work55 were carried out at the BLYP/6311++G(3df,2p) and EOM-SF-CCSD(T)/cc-pVDZ levels of theory. The former is a DFT method employing a large basis set with both polarization and diffuse functions, while the latter is method that accounts for static (via SF ansatz) and dynamic (up to triple excitations) electron correlation. Nonetheless, both of these approaches failed to predict the correct ground state for 2p. Here we revisit the calculations on the singlet-triplet splitting of 2p, summarized in Table 2.

28 ACS Paragon Plus Environment

Page 28 of 38

Page 29 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Table 2. Calculated Values for the Singlet-Triplet Energy Splitting in 2p Computational Method ∆EST, kcal/mola,b B3LYP/6-31+G* B3LYP/6-311++G(3df,2p) BLYP/6-31+G* BLYP/6-311+G(3df,2p) M06-2X/6-31+G* M06-2X/6-311++G(3df,2p) SF-B5050LYP/6-311++G(d,p)

3.2 3.3 1.3 1.5 2.4 2.5 -1.6

EOM-SF-CCSD/cc-pVDZ EOM-SF-CCSD(dT)/cc-pVDZ EOM-SF-CCSD/aug-cc-pVDZ EOM-SF-CCSD(dT)/aug-cc-pVDZ EOM-SF-CCSD/cc-pVTZ EOM-SF-CCSD(dT)/cc-pVTZ EOM-SF-CCSD/aug-cc-pVTZ p-EOM-SF-CCSD(dT)/aug-cc-pVTZ

7.8 2.9 4.8 -0.2 5.7 0.5 4.2 -0.8c -1.0d -1.0e

EOM-SF-CCSD(dT)/aug-cc-pVTZ a

Singlet-state energies, relative to triplet state Positive values indicate a triplet ground state. b DFT energies are for geometries optimized at the respective level of theory; spin-flip energies are for singlet geometries optimized with SF-B5050LYP/6-311++G(d,p) and triplet geometries optimized with B5050LYP/6-311++G(d,p) c Projected value calculated using the EOM-SF-CCSD/aug-cc-pVTZ energy with the dT component calculated with the aug-cc-pVDZ basis set d Projected value calculated using the EOM-SF-CCSD/aug-cc-pVTZ energy with the dT component calculated with the cc-pVTZ basis set e Direct (non-projcted) calculation

We find that standard DFT methods give triplet ground states, regardless of the basis set. As for EOM-SF-CC calculations, we find that increasing the size of the basis set leads to preferential stabilization of the singlet state.

In particular, augmentation of the basis set with diffuse

functions significantly stabilizes the singlet state more than it does the triplet, such that the predicted ground state when using singly-augmented functions is the singlet state. At the best level of theory used in this work, p-EOM-SF-CCSD(dT)/aug-cc-pVTZ, 2p is indeed predicted to have a singlet ground state, with a singlet-triplet energy splitting of 0.8 – 1.0 kcal/mol.

29 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Therefore, accurate prediction of the correct ground state in 2p requires both an accurate account of electron correlation but also a very large basis set that includes diffuse functions. Interestingly, SF-DFT with B5050LYP/6-311++G(d,p) also predicts a singlet ground state, by 1.6 kcal/mol. The predicted energy difference between the singlet and triplet states is sufficiently small that it is possible that there could be some population of triplet state present at thermal temperatures. The spectra shown above were all measured at a temperature of 20K, and therefore, where the population of the excited state should not be significant. We have also measured the 355 nm photoelectron spectrum of 2p at a temperature of 299K. Although there are some differences in the relative peak intensities of the doublet feature, they are not sufficient to be attributed to triplet anion contribution. It is especially difficult to draw any conclusions due to the similarities of the doublet bands obtained from the singlet and triplet states of the anion (see the calculated spectra in Figure 7).

Conclusions Negative ion photoelectron spectroscopy studies have been used to provide ground state assignments for quinonimide anions. Whereas previous experimental and computational studies55 suggested that para-quinonimide, 2p, either has triplet ground state or a low-lying triplet, the photoelectron spectrum presented in this work is best interpreted as showing a ground state singlet. It is likely the case that the triplet state is relatively low in energy, as very high levels of theory, involving an accurate account of electron correlation and triple-zeta basis sets with diffuse functions, are required to predict correctly the singlet ground state, while standard density functional methods are insufficient. However, SF-DFT is able to produce singlet-triplet

30 ACS Paragon Plus Environment

Page 30 of 38

Page 31 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

energy splittings in reasonably good agreement with EOM-SF-CCSD(dT) calculations, and correctly predict a singlet ground state for para. In contrast, the ortho isomer, 2o, is found spectroscopically to have a triplet ground state, which is consistent with destabilization of the singlet state due to electron-pair repulsion, as proposed for ortho-quinone.29 The effect in 2o is likely more pronounced because of the diffuse nature of the lone pair of electrons on the anionic nitrogen. Removing an electron to form the radical, 4o, reduces the electron repulsion. Any effects of repulsion of the in-plane unpaired electron in 4o appear to be minimal, as the relative electron affinities of 4o and 4p are similar to what has been found previously for o- and pbenzoquinone (1o and 1p), 1.90 and 1.85 eV, respectively.28 Similarly, the measured doubletquartet splitting in 4o, 35.5 kcal/mol, is similar to that found for 1o, 36.0 kcal/mol. Consequently, the o- and p-quinoniminyl radicals resemble the corresponding quinones. The meta-quinonimide anion, 4m, is found to rearrange to the para- and other unidentified isomers. The instability of 4m is not surprising. Attempts to generate this anion under flowing afterglow conditions55 were unsuccessful, and reactivity studies of the mquinomethanimide anologue, 6, indicated that it also likely undergoes rearrangement.55

The m-quinonimide anion, 2m, and the radical, 4m, are both meta-phenylene structures and therefore would be expected to have high-spin ground states, which is supported by calculations.

31 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ASSOCIATED CONTENT Supporting Information Calculated spectra for 2o and 2p obtained using a plane wave approximation; comparison of 355 nm spectra of 2p taken at 20K and 299K; calculated structures, energies and frequencies used in this work.

AUTHOR INFORMATION Corresponding Authors [email protected] [email protected] ORCID Paul G. Wenthold: 0000-0002-8257-3907 Xue-Bin Wang: 0000-0001-8326-1780 Ekram Hossain: 0000-0002-0616-3420 Samer Gozem: 0000-0002-6429-2853. Notes The authors declare no competing financial interest.

32 ACS Paragon Plus Environment

Page 32 of 38

Page 33 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

ACKNOWLEDGEMENTS This work was supported by the National Science Foundation (PGW: CHE11-11777 and CHE15-65755). The photoelectron spectroscopy work done at PNNL was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, the Division of Chemical Sciences, Geosciences, and Biosciences (XBW and SHD), and performed using EMSL, a national scientific user facility sponsored by DOE's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. A.I.K. is supported by the Department of Energy through the DE-FG02-05ER15685 grant. A.I.K. is also a grateful recipient of the Bessel Research Award from the Alexander von Humboldt Foundation.

References (1) Allen, J. P.; Feher, G.; Yeates, T. O.; Rees, D. C.; Deisenhofer, J.; Michel, H.; Huber, R. Proc. Natl. Acad. Sci. U. S. A. 1986, 83, 8589-93. (2) Lancaster, C. R. D.; Michel, H. J. Mol. Biol. 1999, 286, 883-898. (3) Deisenhofer, J.; Epp, O.; Miki, K.; Huber, R.; Michel, H. Nature 1985, 318, 61824. (4) Feher, G.; Allen, J. P.; Okamura, M. Y.; Rees, D. C. Nature 1989, 339, 111-16. (5) Kamiya, N.; Shen, J.-R. Proc. Natl. Acad. Sci. U. S. A. 2003, 100, 98-103. (6) Zouni, A.; Witt, H.-T.; Kern, J.; Fromme, P.; Krauss, N.; Saenger, W.; Orth, P. Nature 2001, 409, 739-743. (7) Ferreira, K. N.; Iverson, T. M.; Maghlaoui, K.; Barber, J.; Iwata, S. Science 2004, 303, 1831-1838. (8) Jordan, P.; Fromme, P.; Witt, H. T.; Klukas, O.; Seanger, W.; Krauss, N. Nature 2001, 411, 909-917. (9) Krauss, N.; Schubert, W. D.; Klukas, O.; Fromme, P.; Witt, H. T.; Saenger, W. Nat Struct Biol 1996, 3, 965-73. (10) Loll, B.; Kern, J.; Saenger, W.; Zouni, A.; Biesiadka, J. Nature 2005, 438, 10401044. (11) Srinivasan, N.; Golbeck, J. H. Biochim. Biophys. Acta, Bioenerg. 2009, 1787, 1057-1088. (12) Barreto, G.; Madureira, D.; Capani, F.; Aon-Bertolino, L.; Saraceno, E.; AlvarezGiraldez, L. D. Environ. Mol. Mutagen. 2009, 50, 771-780. (13) Cavalieri, E. L.; Rogan, E. G.; Chakravarti, D. Cell. Mol. Life Sci. 2002, 59, 665681. (14) Crosby, D. G. Plants Toxic to the Skin; Oxford University Press: North Carolina, USA, 2004. (15) Gordon, P. F.; Gregory, P. Organic Chemistry in Color; Springer-Verlag: Berlin Heidelberg, 1983. 33 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(16) Matsuoka, M. Quinone dyes, in Infrared Absorbing Dyes; Matsuoka, M., Ed.; Plenum: New York, NY, 1990, p 35-43. (17) Venkataraman, K.; Iyer, V. N. Anthraquinonoid vat dyes, in Chem. Syn. Dyes; Academic Press: 1971; Vol. 5, p 131-240. (18) Barres, A.-L.; Geng, J.; Bonnard, G.; Renault, S.; Gottis, S.; Mentre, O.; Frayret, C.; Dolhem, F.; Poizot, P. Chem. - Eur. J. 2012, 18, 8800-8812, S8800/1-S8800/11. (19) Chen, H.; Armand, M.; Courty, M.; Jiang, M.; Grey, C. P.; Dolhem, F.; Tarascon, J.-M.; Poizot, P. J. Am. Chem. Soc. 2009, 131, 8984-8988. (20) Liang, Y.; Tao, Z.; Chen, J. Adv. Energy Mater. 2012, 2, 742-769. (21) Nagao, M.; Kobayashi, K.; Yamamoto, Y.; Hibino, T. J. Electrochem. Soc. 2015, 162, F410-F418. (22) Popp, B. V.; Stahl, S. S. Top. Organomet. Chem. 2007, 22, 149-189. (23) Shen, C.-Y.; Yang, T.-K. 1,4-benzoquinone, in Handbook of Reagents for Organic Synthesis: Catalytic Oxidation Reagents; Fuchs, P., Ed.; John Wiley & Sons Ltd.: Chichester, UK, 2013, p 49-52. (24) Yang, T.-K.; Shen, C.-Y. 1,4-Benzoquinone, in e-EROS Encyclopedia of Reagents for Organic Synthesis; John Wiley & Sons, Ltd.: Chichester, UK, 2001, p No pp. given. (25) Nawrat, C. C.; Moody, C. J. Angew. Chem., Int. Ed. 2014, 53, 2056-2077. (26) Chiang, Y.; Kresge, A. J.; Zhu, Y. J. Am. Chem. Soc. 2001, 123, 8089-8094. (27) Chiang, Y.; Kresge, A. J.; Zhu, Y. J. Am. Chem. Soc. 2002, 124, 6349-6356. (28) Fu, Q.; Yang, J.; Wang, X.-B. J. Phys. Chem. A 2011, 115, 3201-3207. (29) Fattahi, A.; Kass, S. R.; Liebman, J. F.; Matos, M. A. R.; Miranda, M. S.; Morais, V. M. F. J. Am. Chem. Soc. 2005, 127, 6116-6122. (30) Chen, B.; Hrovat, D. A.; Deng, S. H. M.; Zhang, J.; Wang, X.-B.; Borden, W. T. J. Am. Chem. Soc. 2014, 136, 3589-3596. (31) Fort, R. C., Jr.; Getty, S. J.; Hrovat, D. A.; Lahti, P. M.; Borden, W. T. J. Am. Chem. Soc. 1992, 114, 7549-52. (32) Recalculated from the results from reference 29 but with the revised value of EA(1m) taken from reference 30. (33) Chacko, S. A.; Wenthold, P. G. Int. J. Mass Spectrom. 2007, 267, 277-283. (34) Yamamura, S. Polycyclic arenes, quinomethanes and large ring systems, including heterocycles, in Rodd's Chemistry of Carbon Compounds (2nd Edition); Topical Volumes; Organic Electrochemistry ed.; Elsevier: Amsterdam, Neth, 2002; Vol. 5, p 225-305. (35) Wan, P.; Brousmiche, D. W.; Chen, C. Z.; Cole, J.; Lukeman, M.; Xu, M. Pure Appl. Chem. 2001, 73, 529-534. (36) Thompson, D. C.; Thompson, J. A.; Sugumaran, M.; Moldeus, P. Chem.-Biol. Interact. 1993, 86, 129-62. (37) Toteva, M. M.; Richard, J. P. Adv. Phys. Org. Chem. 2011, 45, 39-91. (38) Wang, P.; Song, Y.; Zhang, L.; He, H.; Zhou, X. Curr. Med. Chem. 2005, 12, 2893-2913. (39) Willis, N. J.; Bray, C. D. Chem. - Eur. J. 2012, 18, 9160-9173. (40) Segura, J. L.; Martin, N. Chem. Rev. 1999, 99, 3199-3246. (41) Yoshida, H.; Ohshita, J.; Kunai, A. Bull. Chem. Soc. Jpn. 2010, 83, 199-219. (42) Atlas, S. M. Nuova Chim. 1971, 47, 41-6. (43) Charlton, J. L.; Alauddin, M. M. Tetrahedron 1987, 43, 2873-89. 34 ACS Paragon Plus Environment

Page 34 of 38

Page 35 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

(44) Scaiano, J. C.; Wintgens, V.; Netto-Ferreira, J. C. Pure Appl. Chem. 1990, 62, 1557-64. (45) Yamamoto, K. Konan Joshi Daigaku Kenkyu Kiyo 1981, 18, 135-44. (46) Neuhaus, P.; Grote, D.; Sander, W. J. Am. Chem. Soc. 2008, 130, 2993-3000. (47) Berson, J. A. meta-Quinonoid compounds, in Chemistry of Quinoidal Compounds; Pt. 1 ed.; Patai, S., Rappoport, Z., Eds.; Wiley: Chichester, UK, 1988; Vol. 2, p 455-36. (48) Berson, J. A. Non-kekule molecules as reactive intermediates, in Reactive Intermediate Chemistry; Moss, R. A., Platz, M. S., Jones, M., Jr, Eds.; John Wiley & Sons, Inc.: Hoboken, N. J., 2004, p 165-203. (49) Pollack, S. K.; Raine, B. C.; Hehre, W. J. J. Am. Chem. Soc. 1981, 103, 63086313. (50) Hammad, L. A.; Wenthold, P. G. J. Am. Chem. Soc. 2000, 122, 11203. (51) Ichimura, A. S.; Lahti, P. M. Mol. Cryst. Liq. Cryst. Sci. Technol., Sect. A 1993, 233, 33-40. (52) Nicolaides, A.; Tomioka, H. J. Photosci. 2003, 10, 165-173. (53) Olankitwanit, A.; Pink, M.; Rajca, S.; Rajca, A. J. Am. Chem. Soc. 2014, 136, 14277-14288. (54) Rajca, A.; Olankitwanit, A.; Rajca, S. J. Am. Chem. Soc. 2011, 133, 4750-4753. (55) Rau, N. J.; Welles, E. A.; Wenthold, P. G. J. Am. Chem. Soc. 2013, 135, 683-690. (56) Rau, N. J.; Wenthold, P. G. Int. J. Mass Spectrom. 2015, 377, 496-501. (57) Wijeratne, N. R.; Wenthold, P. G. J. Org. Chem. 2007, 72, 9518-9522. (58) Chacko, S. A.; Wenthold, P. G. Mass Spectrom. Rev. 2006, 25, 112-126. (59) Winkler, M. J. Phys. Chem. A 2008, 112, 8649-8653. (60) Wenthold, P. G. J. Org. Chem 2012, 77, 208-214. (61) Geiger, U.; Haas, Y.; Grinstein, D. J. Photochem. Photobiol. A: Chem. 2014, 277, 53-61. (62) Wenthold, P. G.; Lineberger, W. C. Acc. Chem. Res. 1999, 32, 597-604. (63) Chandra, H.; Portwood, L.; Symons, M. C. R. J. Chem. Soc., Faraday Trans. 1 1989, 85, 1801-7. (64) Liu, Z. L.; Wu, L. M.; Liu, Y. C. Magn. Reson. Chem. 1988, 26, 577-80. (65) Krylov, A. I. J. Phys. Chem. A 2005, 109, 10638-10645. (66) Pallagi, I.; Toro, A.; Horvath, G. J. Org. Chem. 1999, 64, 6530-6540. (67) Wang, X. B.; Wang, L. S. Rev. Sci. Instrument 2008, 79, 073108. (68) Ryu, B.-Y.; Emrick, T. Macromol. 2011, 44, 5693-5700. (69) Tomioka, H.; Matsushita, T.; Murata, S.; Koseki, S. Liebigs Ann. 1996, 19711980. (70) Grimme, S.; Waletzke, M. J. Chem. Phys. 1999, 111, 5645-5655. (71) Shao, Y.; Head-Gordon, M.; Krylov, A. I. J. Chem. Phys. 2003, 118, 4807-4818. (72) Krylov, A. I. Acc. Chem. Res. 2006, 39, 83-91. (73) Krylov, A. I. Annu. Rev. Phys. Chem. 2008, 59, 433-462. (74) Slipchenko, L. V.; Krylov, A. I. J. Chem. Phys. 2002, 117, 4694-4708. (75) Krylov, A. I. Chem. Phys. Lett. 2001, 350, 522-530. (76) Levchenko, S. V.; Krylov, A. I. J. Chem. Phys. 2004, 120, 175-185. (77) Manohar, P. U.; Krylov, A. I. J. Chem. Phys. 2008, 129, 194105 - 194110. (78) Piecuch, P. J. Chem. Phys. 2005, 123, 224105. 35 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(79) Beran, G. J. O.; Gwaltney, S. R.; Head-Gordon, M. Phys. Chem. Chem. Phys. 2003, 5, 2488-2493. (80) Cristian, A.-M. C.; Shao, Y.; Krylov, A. I. J. Phys. Chem. A 2004, 108, 65816588. (81) Dunning, T. H., Jr. J. Chem. Phys. 1989, 90, 1007-23. (82) Shao, Y.; Molnar, L. F.; Jung, Y.; Kussmann, J.; Ochsenfeld, C.; Brown, S. T.; Gilbert, A. T. B.; Slipchenko, L. V.et al Phys Chem Chem Phys 2006, 8, 3172-3191. (83) Shao, Y.; Gan, Z.; Epifanovsky, E.; Gilbert, A. T. B.; Wormit, M.; Kussmann, J.; Lange, A. W.; Behn, A.et al Mol. Phys. 2015, 113, 184-215. (84) V.A. Mozhayskiy and A.I. Krylov, ezSpectrum, http://iopenshell.usc.edu/downloads (85) Oana, C. M.; Krylov, A. I. J. Chem. Phys. 2009, 131, 124114. (86) Gozem, S.; Krylov, A. I. ezDyson, http://iopenshell.usc.edu/downloads/ezdyson, 2017. (87) Bethe, H. A.; Salpeter, E. E. Quantum Mechanics of One and Two Electron Atoms; Plenum: New York, 1977. (88) Reed, K. J.; Zimmerman, A. H.; Anderson, H. C.; Brauman, J. I. J. Chem. Phys. 1976, 64, 1368-1375. (89) Ortiz, J. V. Toward an Exact One–Electron Picture of Chemical Bonding, in Advances in Quantum Chemistry; Per-Olov, L., Ed.; Academic Press: 1999; Vol. Volume 35, p 33-52. (90) Oana, C. M.; Krylov, A. I. J. Chem. Phys. 2007, 127, 234106. (91) Head-Gordon, M. Chem. Phys. Lett. 2003, 372, 508-511. (92) Plasser, F.; Wormit, M.; Dreuw, A. J. Chem. Phys. 2014, 141, 024106. (93) Munsch, T. E.; Slipchenko, L. V.; Krylov, A. I.; Wenthold, P. G. J. Org. Chem. 2004, 69, 5735-5741. (94) Vanovschi, V.; Krylov, A. I.; Wenthold, P. G. Theor. Chem. Acc. 2008, 120, 4558. (95) Krylov, A. I. The Quantum Chemistry of Open-Shell Species, in Reviews in Computational Chemistry; John Wiley & Sons, Inc.: Hoboken, NJ, 2017, p 151-224. (96) Cullin, D. W.; Soundararajan, N.; Platz, M. S.; Miller, T. A. J. Phys. Chem. 1990, 94, 8890-6. (97) Karney, W. L.; Borden, W. T. J. Am. Chem. Soc. 1997, 119, 1378 - 1387. (98) Landau, L. D.; Lifshitz, E. M. Quantum Mechanics: Non-relativistic theory; Pergamon: Oxford, 1977. (99) Gozem, S.; Guinina, A. O.; Ichino, T.; Osborn, D. L.; Stanton, J. F.; Krylov, A. I. J. Phys. Chem. Lett. 2015, 6, 4532-4540. (100) Abplanalp, M. J.; Gozem, S.; Krylov, A. I.; Shingledecker, C. N.; Herbst, E.; Kaiser, R. I. Proc. Natl. Acad. Sci. U. S. A. 2016, 113, 7727-7732. (101) Gunion, R. F.; Gilles, M. K.; Polak, M. L.; Lineberger, W. C. Int. J. Mass Spectrom. 1992, 117, 601 - 620. (102) Kim, J. B.; Wenthold, P. G.; Lineberger, W. C. J. Phys. Chem. A 1999, 103, 10833 - 10841. (103) Leopold, D. G.; Miller, A. E. S.; Lineberger, W. C. J. Am. Chem. Soc. 1986, 108, 1379 - 1384.

36 ACS Paragon Plus Environment

Page 36 of 38

Page 37 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

(104) Wenthold, P. G.; Squires, R. R.; Lineberger, W. C. J. Am. Chem. Soc. 1998, 120, 5279 - 5290. (105) Amegayibor, F. S.; Nash, J. J.; Kenttaemaa, H. I. J. Am. Chem. Soc. 2003, 125, 14256-14257. (106) Williams, P. E.; Jankiewicz, B. J.; Yang, L.; Kenttamaa, H. I. Chem. Rev. 2013, 113, 6949-6985.

37 ACS Paragon Plus Environment

Journal of the American Chemical Society

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

For Table of Contents Only

38 ACS Paragon Plus Environment

Page 38 of 38