The Distinctive Electronic Structures of Rhenium ... - ACS Publications

21 Dec 2016 - Edward N. Brothers,. ‡ and Michael B. Hall*,†. † ... methods including complete active space self-consistent field. (CASSCF and CA...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/IC

The Distinctive Electronic Structures of Rhenium Tris(thiolate) Complexes, an Unexpected Contrast to the Valence Isoelectronic Ruthenium Tris(thiolate) Complexes Hao Tang,† Edward N. Brothers,‡ and Michael B. Hall*,† †

Department of Chemistry, Texas A&M University, College Station, Texas 77845, United States Science Program, Texas A&M University at Qatar, Doha, Qatar



S Supporting Information *

ABSTRACT: The noninnocent 2-diphenylphosphino-benzenethiolate (DPPBT) ligand containing both phosphorus and sulfur donors delocalizes the electron density in a manner reminiscent of dithiolenes. The electronic structure of the [ReL3]n (L = DPPBT, n = 0, 1+, 2+) complexes was probed with density-functional theory (DFT) and high-level ab initio methods including complete active space self-consistent field (CASSCF and CASPT2) and coupled cluster (CCSD and CCSD(T)). DFT predicts a slight preference for a closed-shell (CS) singlet ground state for the neutral [ReL3]0 and stronger preferences for low-spin ground states for the oxidized [ReL3]+ and [ReL3]2+. High-level ab initio methods confirm a CS singlet with a Re(III) (d4, S = 0) center as the ground state of [ReL3]0. Thus, this neutral Re species has considerably less thiyl radical character than the valence isoelectronic [RuL3]+, which is mainly a Ru(III) (d5, S = 1/2) anti-ferromagnetically (AF) coupled to a thiyl radical (S = 1/2). However, the oxidized derivatives [ReL3]+ and [ReL3]2+ show significant metal-stabilized thiyl radical character like [RuL3]+. Both [ReL3]+ and [ReL3]2+ have major contributions from Re(III) (d4, S = 1) centers AF coupled to thiyl and dithiyl DPPBT ligands. These findings are consistent with the experimental chemistry as [RuL3]+, [ReL3]+, and [ReL3]2+ can add ethylene to form the new C−S bonds, but [ReL3]0 cannot. The thiyl radicals on the S2 position (the S trans to a P donor) serve as the intrinsic electron acceptors in the actual ethylene addition reactions with Ru and Re tris(thiolate) complexes. Chart 1. Ru and Re Tris(thiolate) Complexesa

I. INTRODUCTION Transition-metal complexes with noninnocent ligands are currently an area of intense study.1−4 The chemistry of dithiolene complexes has played a key role in understanding such ligands.1,5−12 Complexes with other noninnocent ligands, especially those with complicated electron-transfer properties, have received wide attention for a variety of applications from catalysis to optoelectronics.2−4 However, the presence of noninnocent ligands challenges the conventional oxidation state assignments and complicates our ability to explain and predict the behavior of their complexes, because of the strong mixing between the ligand and metal orbitals.1−12 The 2-diphenylphosphino-benzenethiolate (DPPBT) ligand containing both phosphorus and sulfur as the donor atoms was reported by Dilworth et al., and it can delocalize the electron density like dithiolenes.13 Recently, the Ru and Re tris(thiolate) complexes [Ru(DPPBT)3]n ([RuL3]n, n = −, 0, +)14 and [Re(DPPBT)3]n ([ReL3]n, n = 0, +, 2+),15 have been synthesized and structurally as well as spectroscopically characterized, wherein the metal ion is chelated meridionally by the DPPBT ligands (Chart 1). These Ru and Re tris(thiolate) complexes can act as electrocatalysts for the uptake and release of olefins and © XXXX American Chemical Society

a

These complexes are nearly octahedral geometries, where the d orbital splits into t2g orbitals (point between ligands) and eg orbitals (point directly at the ligands).

hydrogen,10,14c,e,f,15 like similar reactions in metal dithiolenes.4−10,12 Our knowledge of the noninnocent DPPBT ligand and how it supports this chemistry remains scant. Among the [RuL3]n complexes, it is the cationic species [RuL3]+ that is known to readily undergo addition reactions with a variety of unsaturated organic compounds.14 Although the anionic [RuL3]‑ complex is Received: October 12, 2016

A

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 1. Schematic MO diagram for [RuL3]+ and alternative electronic structures for its ground state.

suite of program.24 The geometric structures of the molecules [ReL3]0/+/2+, where L = DPPBT = 2-diphenylphosphinobenzenethiolate (full ligand), were optimized in the gas phase (GP) with ω-B97XD25 functional in conjunction with the SDD(f)26(Re), 6-311G* (P, S, C, and H atoms), and 6-31G (C and H atoms on phenyl group) basis set (B1). Frequency calculations were performed for all the optimized stationary points to ensure that they are true minima. Solvation was incorporated by using the SMD27 method for all the singlepoint calculations with dichloromethane as the solvent. The SDD(f) (Re) and 6-311++G**(all other atoms) basis set (B2) with the same functional was employed for the single-point energy calculations based on the optimized structures. The above methodology was tested in a variety of ways. The various tests and their results are summarized in the Supporting Information, while the results and discussion in the text refer to free energies (including the B2 electronic energies in solution, GP thermal correction to the Gibbs free energy, and solvation free energy) based on ω-B97XD/B2/SMD//ω-B97XD/B1/ GP for the full-ligand molecules [ReL3]0/+/2+. DFT Calculations for Truncated Models. Because of the high computational expense of the higher-level ab initio methods for the full-ligand structures of [ReL3]0/+/2+, we employed the truncated models (denoted [ReL′3]0/+/2+) as shown in Chart 2. Geometries, frequencies, and free energies at 298.15 K for the truncated models [ReL′3]0/+/2+ were also computed using the ω-B97XD functional in combination with the SDD(f) (Re)/6-311G*(rest) basis set (B4). The singlepoint calculations were performed with the SDD(f) (Re)/6-

clearly a d6 Ru(II) trithiolate, the diamagnetic [RuL3]+ species is ambiguous with respect to the oxidation state of the metal ion as well as that of the DPPBT ligands. Because density functional theory (DFT) calculations suggested nearly degenerate closed-shell (CS) singlet and triplet states for [RuL3]+ (Figure 1a,c), Grapperhaus et al. initially proposed a singlet diradical, with the electronic configuration (S2− S3)2(π*xz)1(π*yz)1 (Figure 1b) as the ground state.14c,h In an alternative assignment that is more obviously consistent with the alkene addition reaction, an S2−S3 diradical with a low-spin Ru(II) was also proposed (Figure 1d).14h Very recently, highlevel ab initio methods showed that the ground state for [RuL3]+ was actually a singlet diradical state from the antiferromagnetic (AF) coupling between Ru(III) (S2 = 1/2) and S pz (S = 1/2; Figure 1e).16 These studies provided information about alternative spin-state and the relationship between the chemistry of the [RuL3]+ complex and this new electronic description as a Ru(III)−thiyl diradical. Although [RuL3]+ readily adds ethylene, the valence isoelectronic [ReL3]0 species does not.14c,e,15a In contrast, it is the one- and two-electron oxidized derivatives [ReL3]+ and [ReL3]2+ that bind to ethylene.14c,15a Knowledge of the differences in the electronic structures of the related rhenium derivatives [ReL3]0/+/2+ that might explain the different reactivity are still missing. In this work, detailed explorations of the electronic structures of [ReL3]0/+/2+ by DFT and high-level ab initio methods including CASSCF,17 CASPT2,18 MP2,19 MP3,20 MP4D, MP4DQ, MP4SDQ,21 CCSD, CCSD(T),22 QCISD, and QCISD(T)20,23 provide insight into the relationship of these complexes with the Ru(III)−thiyl diradical [RuL3]+.16 In addition to clarifying the geometric structure, electronic structure, and spin-state ordering of [ReL3]0/+/2+, the calculations provide explanations for the different reactivity observed for the ethylene addition reactions with [RuL3]+ and [ReL3]0/+/2+.

Chart 2. Truncated Models

II. COMPUTATIONAL METHODS DFT Calculations for Full-Ligand Molecules. All the DFT calculations were performed by using the Gaussian 09 B

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 2. DFT(ω-B97XD)-optimized and experimental (Exp)13a geometric parameters are given with bond lengths (Å) and angles (deg) for the full-ligand molecules [ReL3]0/+/2+ (A−C) and the truncated models [ReL′3]0/+/2+ (D−F). The reported energies (kcal/mol) include ω-B97XD/ SMD electronic energy using dichloromethane as the solvent with GP thermal free energy correction at 298 K.

311++G**(rest) basis set (B2) with the same ω-B97XD functional based on the optimized structures. The continuum solvation model and solvent were same as those for the fullligand molecules. The reported energy in the main text for the truncated model [ReL′3]0 includes the B2 electronic energies in solution, GP thermal correction to the Gibbs free energy, and solvation free energy based on ω-B97XD/B2/SMD//ωB97XD/B4/GP. The other tested TPSS and M06 functionals showed that the trends for the relative energies of the low and high spin states and geometric parameters were identical to those obtained from the ω-B97XD functional (Tables S11 and S12). Ab Initio Calculations. Single-point calculations using the Hartree−Fock (HF), PUHF,28 MP2, PMP2, MP3, PMP3, MP4D, MP4DQ, MP4SDQ, CCSD, CCSD(T), QCISD, and QCISD(T) methods were done in Gaussian 09 suite of program. Multireference ab initio CASSCF and CASPT2 calculations were performed with the MOLPRO suite of program29 using the basis set SDD(f) (Re)/6-311G*(rest) (B4). All the CASSCF calculations were based on the geometries optimized at the ω-B97XD/B4 level. Boys localization was applied to CASSCF active space orbitals, and then the CASPT2 calculations were performed with these localized orbitals. A level shift of 0.2 au was employed in the CASPT2 calculation to avoid the intruder state problem.18 The details of the CASSCF and CASPT2 calculations of [ReL′3]0/+/2+ were relegated to the Supporting Information.

III. RESULTS DFT Results. Like the [RuL 3] − case, wherein the experimental geometric parameters of the full-ligand molecule [RuL3]− are well-reproduced by the DFT calculations,16 comparisons of the geometric parameters of the CS singlet and triplet states for the full-ligand molecule [ReL3]0 against the experimental values (Figure 2A) clearly show that the computed geometric values for the CS singlet state are in closer accord with experimental ones than those of the triplet state. Upon oxidation, the Re−S1 bond distances in the low-spin state for [ReL3]0, [ReL3]+, and [ReL3]2+ become significantly shorter, changing from 2.54 to 2.39 and 2.29 Å, respectively. In contrast, the Re−S2 and Re−S3 bond distances for [ReL3]0, [ReL3]+, and [ReL3]2+ are almost unaffected by the change in the oxidation state. However, the reverse trend was observed in the high-spin state, wherein the Re−S1 bond distances for [ReL3]0, [ReL3]+, and [ReL3]2+ were almost unchanged, but both the Re−S2 and Re−S3 bond distances significantly shortened upon oxidation. Inspection of the free energies between the low- and highspin states (Figure 2) shows that the neutral species [ReL3]0 has an S = 0 ground state with a close-lying S = 1 state (2.5 kcal/mol less stable by ω-B97XD). By comparison, the cationic species [ReL3]+ has an S = 1/2 ground state and an S = 3/2 excited state that is more than 10 kcal/mol higher in energy. For the dicationic species [ReL3]2+, the S = 0 ground state is preferred over the high-spin state by 7.6 kcal/mol. Similar C

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 3. Schematic DFT MO diagrams for (a) CS singlet and (b) triplet states of [ReL3]0, (c) doublet state of [ReL3]+, and (d) CS singlet state of [ReL3]2+.

S1 py orbitals, the Re−S1 bond distance progressively shortens in the low-spin state from [ReL3]0 to [ReL3]2+ (Figure 2). As expected from its anticipated lower oxidation state, the 5d metal complex [ReL3]0 has more metal character in its higherenergy π*xz/π*yz orbitals and correspondingly more ligand S in its πxz/πyz orbitals than the valence isoelectronic 4d metal complex, [RuL3]+ (compare Figure 1 and Figure 3). This difference is reflected by the metal/sulfur ratio in the πxz orbital of 4/1 and 1/2 for [RuL3]+ and [ReL3]0 (Table S10). Likewise, the energy difference between the nonbonding dxy and the nonbonding S2−S3 orbitals decreases from [RuL3]+ (1.13 eV) to [ReL3]0 (0.28 eV), and this difference is smaller than that in [RuL3]+ even for the more oxidized [ReL3]+ (α: 0.55 eV and β: 0.40 eV) and [ReL3]2+ (0.68 eV). High-Level ab Initio Method Results. To reduce the computational cost of the high-level ab initio methods for the full-ligand systems [ReL3]0/+/2+, we used the truncated models [ReL′3]0/+/2+ (Figure 2) that was successfully employed in our previous study of [RuL3]+.16 The optimized structural parameters of the truncated models [ReL′3]0/+/2+ are in reasonable agreement with the full-ligand molecules as presented in Figure 2. The Re−S and Re−P bonds of the truncated models [ReL′3]0/+/2+ calculated at the ω-B97XD/B4 level differ from the full-ligand systems [ReL3]0/+/2+ by ≤0.08 Å, and the P2−Re−P3 angle differs by ≤7° (Table S12 compares these deviations calculated by both TPSS and M06 functionals). The truncated models follow the same trends of decreasing Re−S1 bond lengths and almost unchanged Re−S2/ S3 bond length as the complexes are oxidized from the neutral to the dication, and they reproduce the various trans effects as observed in the full-ligand molecules for different spin states (Figure 2). Moreover, the calculated free-energy difference between the low- and high-spin states for the truncated model [ReL′3]0 (3.1 kcal/mol) is in good agreement with the lowand high-spin state difference in the full-ligand molecule [ReL3]0 (Figure 2). In addition, there are no appreciable differences in the frontier MOs between the full and truncated molecules by visual inspection of the Kohn−Sham orbitals

trends for these energetic gaps were obtained from the other tested density functionals (Table S1). As expected, more HF exact exchange (with the same correlation functional) tends to stabilize the triplet state relative to the CS singlet state (Table S1). Among all the tested DFT functionals (Table S1), MN12SX predicts the most stable CS singlet with the triplet 5.9 kcal/mol higher, a result consistent with the high-level ab initio methods (see below). Because of the large energy separations for the cations, we will only consider the low-spin states for [ReL3]+ and [ReL3]2+ below. The performance of DFT functionals was also assessed by calculating the redox couples [ReL3]+/0 and [ReL3]2+/+ relative to Fc+/0, and the B97D functional performs best among all the tested functionals (Table S9). Like the DFT molecular orbitals (MOs) in [RuL3]+ (Figure 1), the upper valence MOs in [ReL3]0/+/2+ involve strong mixing between the Re 5d and the S1, S2, and S3 3p atomic orbitals (AOs). As shown in Figure 3, one pair of orbitals, πxz/ π*xz, is built from the bonding/antibonding combination of Re dxz and the coplanar, in-phase combination of S2 pz and S3 px. The other pair, πyz/π*yz, is built from the bonding/antibonding combination of Re dyz and S1 py. The Re dx2−y2 and Re dz2 are two metal-centered antibonding orbitals that have σ symmetry with respect to the P and S donor atoms. The Re dxy and the out-of-phase S2−S3 are essentially nonbonding. The πxz/π*xz orbitals in the xz plane are almost orthogonal to the πyz/π*yz orbitals in the yz plane. As mentioned above, the CS singlet and triplet state in [ReL3]0 (Figure 3a,b) are close in energy, an issue that will be examined below. Regardless of the ground state of [ReL3]0, oxidation to [ReL3]+ creates a doublet state with one electron in the π*yz orbital (Figure 3c), while, as mentioned above, the corresponding higher-spin quartet is much higher in energy. Further oxidation to [ReL3]2+ removes the remaining electron from the π*yz, yielding a CS ground state (Figure 3d) with a triplet somewhat higher in energy than that for [ReL3]0. Consistent with the nature of the π*yz, antibonding between the Re dyz and D

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry (Figure S1). Therefore, the truncated models [ReL′3]0/+/2+, which are used in the high-level methods to assess the intrinsic electronic structure of the Re center, should be regarded as good models for the full-ligand molecules [ReL3]0/+/2+. The degree of multireference character in these rhenium complexes is shown by the T1 diagnostic values for these singlet and triplet states of the neutral species, denoted 1,3 [ReL′3]0; doublet state of the cationic species, denoted 2 [ReL′3]+; and singlet state of the dicationic species, denoted 1 [ReL′3]2+ in Table 1. The values, all somewhat larger than 0.02, suggest some multireference character in all of these species.

electronic configurations, the electron occupancy numbers, key natural orbitals, and relative energies for the CASSCF(16,13) calculations for 1,3[ReL′3]0 based on the DFT-optimized geometries. The singlet state of [ReL′3]0 is dominated by two important configurations (Ψ1 and Ψ2) that are related by double (two electrons) excitation from orbital ψ6 to orbital ψ5 (Scheme 1B). Although this appears to be somewhat like the valence isoelectronic [RuL′3]+ (Scheme 1A),16 there are some important differences. First, from an inspection of the key orbitals for the singlet state, as shown in Scheme 1A,B, it is seen that the orbitals denoted as ψ5 and ψ6 in 1[ReL′3]0 are antibonding and bonding MOs between the Re dxz and the inphase combination of the S2 pz and S3 px orbitals, while the corresponding orbitals in 1[RuL′3]+ involve strong double excitation from the bonding MO between the Ru dxz and S2 pz to an antibonding MO.16 Second, both a comparison of the coefficients for the two most leading configurations and of the resulting electron occupancy numbers between 1[RuL′3]+ and 1 [ReL′3]0 (Scheme 1A,B) are consistent with a weaker contribution from this double-excitation for [ReL′3]0 than for [RuL′3]+. Thus, the overall electronic structure for the neutral [ReL′3]0 is best assigned as a CS singlet state with a Re(III) center with much less thiyl radical character than the valence isoelectronic [RuL′3]+. The CASSCF wave function for the low-lying triplet state of [ReL′3]0, as shown in Scheme 1D, has one major configuration, in which the singly occupied orbitals ψ5 and ψ6 are mainly metal Re dxz and Re dyz. This result is in contrast with that of the lowlying triplet state of [RuL′3]+, wherein the singly occupied orbitals ψ3 and ψ4 have mainly metal Re dxz and ligand S1 py characters (Scheme 1C).16 Therefore, the triplet state of the neutral species [ReL′3]0 involves a triplet Re(III) center with three thiolate ligands, while that for [RuL′3]+ involves a doublet Ru(III) center with two thiolate ligands and one thiyl ligand, a result that parallels the overall character of their ground states. As shown in Scheme 1, the CASSCF(16,13) calculations predict the singlet ground state of the neutral [ReL′3]0 to be preferred by 13.5 kcal/mol over the triplet state. The smaller CASSCF(14,12) and CASPT2(14,12) calculations, formed by removing a doubly occupied orbital, predict that the singlet state is favored over the triplet state by 16.4 and 24.5 kcal/mol, respectively (Figure S6 and Table S17). Predictions of this energy difference by other ab initio methods on the ω-B97XDoptimized geometries are collected in Table 2. As expected, HF predicts the triplet as the most stable, while MP2 appears to overestimate the stability of the singlet relative to the triplet state, and then MP3 underestimates this difference. However, like the CASSCF and CASPT2, the other (higher-order) methods: MP4SDQ, CCSD, CCSD(T), QCISD, and QCISD(T) predict somewhat similar singlet stabilities. Electronic Structure of the Cationic Species 2[ReL′3]+. For the cation, [ReL′3]+, Scheme 2 depicts the coefficients of the main electronic configurations, the electron occupation numbers, key natural orbitals, and localized orbitals for the

Table 1. T1 Diagnostic Values from CCSD(T1diag) Calculations for [ReL′3]0/+/2+ 1

T1

[ReL′3]0

3

[ReL′3]0

0.028

0.032

2

[ReL′3]+ 0.037

1

[ReL′3]2+ 0.027

Electronic Structures of the Neutral Species [ReL′3]0. Scheme 1 displays the coefficients of the main

1,3

Scheme 1. Main Electronic Configurations (Ψi), Electron Occupancy Numbers (blue), Key Natural Orbitals (ψi), and Relative Energies for CASSCF(16,13) of the Singlet and Triplet States for [ReL′3]0 a

a

To compare the electronic structures between [RuL’3]+ and [ReL’3]0 at the similar level, the CASSCF(16,11) with natural orbitals for the singlet and triplet states of [RuL′3]+ were recalled from ref 16.

Table 2. Relative Single-Point Energies (kcal/mol) of the Singlet and Triplet States of [ReL′3]0 Calculated by Various Electronic Structure Methodsa singlet triplet a

HF

MP2

MP3

MP4SDQ

CCSD

CCSD(T)

QCISD

QCISD(T)

13.9 0.0

0.0 19.8

0.0 0.3

0.0 13.3

0.0 4.7

0.0 8.5

0.0 4.5

0.0 8.8

Also see Table S28. E

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

complicated multi-configuration character of 2[ReL′3]+ in the localized-orbital representation (CASPT2, Scheme 2D) shows a second important configuration with the coefficient of 0.41 that represents a Re(IV)−thiolate contribution with the doubly occupied ϕ3 orbital and virtual ϕ1 orbital (for further details see Figure S7 and Table S19). The occupancy numbers of the CASPT2 localized orbitals depicted in Scheme 2D show that the net unpaired spin remains in the Re dyz orbital (ϕ2), while the occupation numbers of 0.87 and 1.40 for the respective Re dxz (ϕ1) and S2 pz (ϕ3) also suggest a strong mixture of Re(IV)−thiolate and Re(III)−thiyl character in the 2[ReL′3]+ ground state. Thus, the overall 2[ReL′3]+ has much less thiyl character than 1[RuL′3]+. Interestingly, this prediction is in excellent agreement with the recent mass spectrometric study on the reactivity and radical natures of the cationic species [ReL3]+ and [RuL3]+, in which the radical of [ReL3]+ resides in the same S2 position as [RuL3]+, but the reactivity of [ReL3]+ with dimethyl disulfide was lower in comparison to that of [RuL3]+.14i Electronic Structure of the Dicationic Species 1 [ReL3]2+. The multi-configuration character of the dicationic species [ReL′3]2+ as reflected by the CASSCF(10,7) and CASPT2(10,7) calculations is illustrated in Scheme 3, which

Scheme 2. Main Electronic Configurations (Ψi), Electron Occupancy Numbers (blue), Key Natural Orbitals (ψi and Φi) for CASSCF(11,7), and Localized Orbitals (ϕi) for CASPT2(11,7) of Doublet State for [ReL′3]+ a

Scheme 3. Main Electronic Configurations (Ψi), Electron Occupancy Numbers (blue), Key Natural Orbitals (ψi and Φi) for CASSCF(10,7), and Localized Orbitals (ϕi) for CASPT2(10,7) of Singlet State for [ReL′3]2+

a

To compare the electronic structures between [RuL’3]+ and [ReL’3]+ at the similar level, the CASSCF(10,6) with natural orbitals and CASPT2(10,6) with localized orbitals for [RuL’3]+ were recalled from ref 16.

CASSCF(11,7) and CASPT2(11,7) calculations based on the DFT-optimized geometries. Inspection of the CASSCF results with natural orbitals as shown in Scheme 2B, the doublet state of [ReL′3]+ is well-represented by the two most significant configurations Ψ1 and Ψ2 with the respective coefficients of 0.97 and −0.20, which is the contribution from the doubleexcitation from the ψ3 orbital to ψ1 orbital. As shown in Scheme 2B, the orbitals ψ3 and ψ1 are the bonding and antibonding π interaction between the Re dxz and S2 pz orbitals in the xz plane. The orbital ψ2, the Re dyz orbital with a small antibonding contribution from the S1 py orbital, holds essentially all the excess spin as it is singly occupied in both configurations. Although both [ReL′3]+ and [RuL′3]+ are related by a similar double-excitation from the metal dxz and S2 pz bonding MO to its antibonding counterpart, [ReL′3]+ has less multi-configurational character than [RuL3]+, as indicated by the smaller coefficient for the configuration Ψ2 in [ReL′3]+ relative to that in [RuL′3]+.16 Consistent with this two-states description, the electron occupancy numbers of the antibonding MO ψ1 for [ReL′3]+ is smaller than that for [RuL′3]+ (compare Scheme 2A,B). Additionally, the leading configurations and the important orbitals from the larger CASSCF(15,13) calculation for the doublet state of [ReL′3]+ (Figure S3 and Table S14) are essentially identical to those obtained from the CASSCF(11,7) level. While the CASPT2 results for the singlet [RuL′3]+ show a single dominative configuration that represents the AF coupling of the Ru(III)−thiyl diradical (Scheme 2C),16 the doublet [ReL′3]+ has less thiyl radical character, as reflected by the smaller coefficient for configuration Φ1, a Re(III) (S = 1) AF coupled to the thiyl S (S = 1/2; Scheme 2D). The more

displays the coefficients of the main electronic configurations, the electron occupancy numbers, and key natural and localized orbitals. The CASSCF(10,7) results for the singlet state of [ReL′3]2+ (Scheme 3A) show three important configurations (Ψ1−Ψ3). The leading configuration has two highly delocalized Re−S π bonding orbitals (ψ3 = Re dyz + S1 py and ψ4 = Re dxz + F

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Scheme 4. Electronic Structures of 1,3[RuL′3]+, 1,3[ReL′3]0, 2 [ReL′3]+, and 1[ReL′3]2+ from CASPT2 Calculation

S2 pz) doubly occupied, while the next two most important configurations represent double excitations from these bonding MOs to their antibonding counterparts (ψ 1 and ψ 2 , respectively); see Scheme 3A. These double excitations result in the electron occupancy numbers of 1.87 (ψ3) and 1.89 (ψ4) for the bonding MOs, and 0.13 (ψ1) and 0.13 (ψ2) for the antibonding MOs. The larger CASSCF(16,12) calculation for the singlet state of [ReL′3]2+ (Figure S4 and Table S15) has leading configurations and key orbitals much like those described above. Analysis of the localized orbitals of [ReL′3]2+ from the CASPT2 calculation presents a more prominent multiconfigurational character of the ground state with three major configurations and a large number of other significant contributions (Figure S8 and Table S21). As shown in Scheme 3B, the principal configuration Φ1 with a coefficient of 0.52 is a Re(III)−thiyl double diradical with two singly occupied Re d orbitals (Re dyz (ϕ1) and Re dxz (ϕ2)) AF coupled to two singly occupied S orbitals (S1 py (ϕ3) and S2 pz (ϕ4)). The next two most important configurations, Φ2 and Φ3, with comparable coefficients of −0.35 and −0.34, both represent Re(IV) contributions, in which first one and then the other of the Re(III)−thiyl diradicals becomes a thiolate. Several of the next most important configurations (Figure S8 and Table S21) are other Re(III)−thiyl double diradicals. The CASPT2 calculation produces electron occupancy numbers of 0.84, 1.20, 1.02, and 1.30 for Re dyz (ϕ1), S1 py (ϕ3), Re dxz (ϕ2), and S2 pz (ϕ4), respectively, which leads to the conclusion that the doublet [ReL′3]2+ is best described as a Re(III) center with two singly occupied Re dxz and dyz orbitals AF coupled to an overall DPPBT ligand framework with two singly occupied S orbitals (mainly S2 pz and S1 py). Although the cation [ReL′3]+ and dication [ReL′3]2+ show similar thiyl radical character at S2, the dication is predicted to have significant thiyl-radical character at S1.

one electron of a Re(III) (d4, S = 1) center to produce an overall doublet state). Finally for the [ReL′3]2+, imagine further oxidization of the Re center to Re(IV), which again is so unstable that an electron moves from S1 to Re, creating a Re(III) center (d4, S = 1) with both unpaired Re electrons AF coupled back to thiyl radicals as a “double diradical” (Scheme 4f). In general, all the [ReL′3]0/+/2+ complexes are mainly d4 Re(III), while [RuL′3]+ is mainly a d5 Ru(III). As elaborated above, [RuL′3]+, [ReL′3]+, and [ReL′3]2+ are therefore assigned as the metal-stabilized thiyl radicals, but [ReL′3]0 is not, even though it is valence-isoelectronic with [RuL′3]+. This difference originates in the stability of the d4 Re(III) relative to a d4 Ru(IV). Differential Reactivity of Ethylene Addition to [RuL3]+ and [ReL3]0/+/2+. Experimentally, [RuL3]+, [ReL3]+, and [ReL3]2+ can add ethylene across the S2/S3 atoms to form the new C−S bonds, but [ReL3]0 cannot.14c,e The differential reactivity of these complexes toward olefins and their multireference (thiyl-radical) character are both connected with (i) the different electronegativity of Re versus Ru (group 7 versus 8), (ii) the different row (larger 5d versus 4d orbitals), and (iii) the charge and resulting oxidation state. Although [ReL3]0 and [RuL3]+ are valence-isoelectronic, the former is a CS singlet ground state with little radical character on the ligand, while the latter has significant radical character in ligand S2, where its reactivity originates. However, consistent with the experimental results, the cationic [ReL3]+ and dicationic [ReL3]2+ have sufficient radical character on ligand S2 to support their reaction with ethylene. These results indicate that the pivotal role of the metal-stabilized thiyl radicals in affecting the reactivity of the alkenes addition reaction with the metal− thiolate complexes. Furthermore, the lower electronegativity of

IV. DISCUSSION Nature of the Metal−Thiolate Complexes [ReL3]0/+/2+ and [RuL3]+. A summary analysis of the electronic structures of the [RuL′3]+ and [ReL′3]0/+/2+ species from the CASPT2 calculations with the localized orbitals is shown in Scheme 4. For the singlet state of [RuL′3]+, imagine two successive oxidations on the Ru center beginning from [RuL′3]− (Ru(II) trithiolate), which would produce a Ru(IV) trithiolate; however, Ru(IV) is so unstable that an electron moves from S2 to Ru but remains AF coupled to the S2 electron left behind (Scheme 4a). As shown in Scheme 4b, the instability of Ru(IV) also holds true for its triplet state, where now one electron from S1 moves to Ru to create ferromagnetically (triplet) coupled Ru(III)-S1 electrons. In contrast to this description for [RuL′3]+, the [ReL′3]0 analogue clearly favors the Re(III) d4 electron configuration, as there is no significant ligand radical character for the singlet state (Scheme 4c). This difference may be attributed to the stability of Re(III) oxidation state and the lower electronegativity of Re compared to Ru, so that [ReL′3]0 is a trithiolate. Likewise, the triplet state of [ReL′3]0 simply involves a promotion of an electron on Re, creating a Re(III) (d4, S = 1) center with two singly occupied dxz and dyz orbitals (Scheme 4d). Now imagine forming [ReL′3]+ by oxidization of [ReL′3]0 on the Re center, which would produce a Re(IV) trithiolate, but Re(IV) is so unstable that an electron moves from S2 to Re but remains AF coupled to the S2 electron left behind as shown in Scheme 4e (a thiyl radical AF coupled to G

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry Re and lower charge of [ReL3]0 result in Re−S π-bonds more polarized toward sulfur in comparison with the [RuL3]+ case, so that the Re−S π-antibonding orbitals are higher in energy, more polarized toward Re, and thus, not easily attacked by the ethylene. Upon oxidation of the [ReL3]0 complex, these polarization effects begin to shift the electron distribution progressively toward the [RuL3]+ case. Overall, the better overlap between the S 3p and the Re 5d over that with the Ru 4d results in stronger π-bonding between Re and S compared to those for Ru and in somewhat less diradical character in all the Re complexes.



V. CONCLUSIONS In conclusion, the geometries, electronic structures, and spinstate energetics of the Re tris(thiolate) complexes, [ReL3]0/+/2+, were studied by DFT and high-level ab initio quantum calculations and compared with those of [RuL3]+, which is valence-isoelectronic with [ReL3]0. DFT results reveal that in the neutral [ReL3]0 case, the singlet and triplet states are close in energy with a slight preference for low spin. In the oxidized [ReL3]+ and [ReL3]2+ cases, the low-spin ground states are preferred more strongly. Consistent with the DFT, CCSD(T) and QCISD(T) calculations show that the neutral species [ReL3]0 has a singlet ground state. High-level ab initio methods predict that the neutral species [ReL3]0 is well-described as a CS singlet ground state involving a Re(III) (d4, S = 0) center and bound to three thiolate-DPPBT ligands. Thus, [ReL3]0 has obviously less thiyl-radical character than the valence isoelectronic [RuL3]+, which has its major component arising from a Ru(III) (d5, S = 1/2) AF coupled to one thiyl- with two thiolate-DPPBT ligands.16 In contrast, the oxidized derivatives [ReL3]+ and [ReL3]2+ have significant metal-stabilized thiyl-radical character, arising from Re(III) (d4, S = 1) centers AF coupled to one thiyl- with two thiolateDPPBT ligands for [ReL3]+ and to two thiyl- with one thiolateDPPBT ligand for [ReL3]2+. Thus, the [ReL3]0/+/2+ complexes favor the Re(III) d4 electron configuration, while [RuL3]+ disfavors the Ru(IV) d4 electron configuration in favor of Ru(III), low-spin d5. Our findings clarify experiments that show ethylene reacting with [RuL3]+, [ReL3]+, and [ReL3]2+ to form the new C−S bonds but not reacting with [ReL3]0. Besides the different electronegativity and charge of Re versus Ru, the thiyl-radical character is considered as an important factor contributing to the differential reactivity, wherein [RuL3]+ and [ReL3]+/2+ have significant radical characters on the ligand S2, whereas the neutral [ReL3]0 has little thiyl radical character. In fact, the thiyl radical on the S2 position (the S trans to a P donor) serve as the intrinsic electron acceptor in the actual alkene addition reactions. The detailed mechanisms of the ethylene addition to [RuL3]+ and [ReL3]0/+/2+ are under study. The detailed results of the Re tris(DPPBT) species presented in this work should lead to a better appreciation of the noninnocent nature of the DPPBT ligand and its behavior in promoting various reactions in transition-metal complexes.



Computational methods and references; benchmarking of different DFT functionals for relative free energies, geometric structures optimized in GP and solution for the full-ligand molecules, and redox potentials; relative free energies, geometric structures, and DFT frontier MOs of truncated models; orbital compositions; CASSCF and CASPT2 results at the DFT- and CASSCF-optimized geometries for the active spaces, all the configurations, coefficients, and electron occupancy numbers; relative energies of other single-reference methods; Cartesian coordinates and absolute energies (PDF)

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Author Contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We acknowledge financial support from the Qatar National Research Fund under NPRP Grant No. 05-318-1-063 and the Welch Foundation under Grant No. A-0648. Computer time was provided by the TAMU High Performance Research Computing Facility.



REFERENCES

(1) For recent reviews on dithiolenes, see: (a) Eisenberg, R.; Gray, H. B. Noninnocence in Metal Complexes: A Dithiolene Dawn. Inorg. Chem. 2011, 50, 9741−9751. (b) Eisenberg, R. Trigonal Prismatic Coordination in Tris(dithiolene) Complexes: Guilty or Just Noninnocent? Coord. Chem. Rev. 2011, 255, 825−836. (c) Ray, K.; Petrenko, T.; Wieghardt, K.; Neese, F. Joint Spectroscopic and Theoretical Investigations of Transition Metal Complexes Involving Non-innocent Ligands. Dalton Trans. 2007, 1552−1556. (d) Sproules, S.; Wieghardt, K. Dithiolene radicals: Sulfur K-edge X-ray Absorption Spectroscopy and Harry’s Intuition. Coord. Chem. Rev. 2011, 255, 837−860. (2) For recent reviews on noninnocence, see: (a) Lyaskovskyy, V.; de Bruin, B. Redox Non-Innocent Ligands: Versatile New Tools to Control Catalytic Reactions. ACS Catal. 2012, 2, 270−279. (b) Suarez, A. I. O.; Lyaskovskyy, V.; Reek, J. N.; van der Vlugt, J. I.; de Bruin, B. Complexes with Nitrogen-Centered Radical Ligands: Classification, Spectroscopic Features, Reactivity, and Catalytic Applications. Angew. Chem., Int. Ed. 2013, 52, 12510−12529. (c) van der Vlugt, J. I. Cooperative Catalysis with First-Row Late Transition Metals. Eur. J. Inorg. Chem. 2012, 2012, 363−375. (d) Gunanathan, C.; Milstein, D. Metal−Ligand Cooperation by Aromatization−Dearomatization: A New Paradigm in Bond Activation and “Green” Catalysis. Acc. Chem. Res. 2011, 44, 588−602. (e) Gunanathan, C.; Milstein, D. Applications of Acceptorless Dehydrogenation and Related Transformations in Chemical Synthesis. Science 2013, 341, 1229712. (f) Gunanathan, C.; Milstein, D. Bond Activation and Catalysis by Ruthenium Pincer Complexes. Chem. Rev. 2014, 114, 12024−12087. (g) Annibale, V. T.; Song, D. T. Multidentate Actor Ligands as Versatile Platforms for Small Molecule Activation and Catalysis. RSC Adv. 2013, 3, 11432− 11449. (h) Luca, O. R.; Crabtree, R. H. Redox-active Ligands in Catalysis. Chem. Soc. Rev. 2013, 42, 1440−1459. (i) van der Vlugt, J. I.; Reek, J. N. H. Neutral Tridentate PNP Ligands and Their Hybrid Analogues: Versatile Non-Innocent Scaffolds for Homogeneous Catalysis. Angew. Chem., Int. Ed. 2009, 48, 8832−8846. (j) Tezgerevska,

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.inorgchem.6b02434. H

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry T.; Alley, K. G.; Boskovic, C. Valence Tautomerism in Metal Complexes: Stimulated and Reversible Intramolecular Electron Transfer Between Metal Centers and Organic Ligands. Coord. Chem. Rev. 2014, 268, 23−40. (k) Li, H.; Hall, M. B. Computational Mechanistic Studies on Reactions of Transition Metal Complexes with Noninnocent Pincer Ligands: Aromatization−Dearomatization or Not. ACS Catal. 2015, 5, 1895−1931. (l) Boyer, J. L.; Rochford, J.; Tsai, M.-K.; Muckerman, J. T.; Fujita, E. Ruthenium Complexes with Non-innocent Ligands: Electron Distribution and Implications for Catalysis. Coord. Chem. Rev. 2010, 254, 309−330. (m) Blanchard, S.; Derat, E.; Desage-El Murr, M.; Fensterbank, L.; Malacria, M.; MourièsMansuy, V. Non-Innocent Ligands: New Opportunities in Iron Catalysis. Eur. J. Inorg. Chem. 2012, 2012, 376−389. (n) Ward, M. D.; McCleverty, J. A. Non-innocent Behaviour in Mononuclear and Polynuclear Complexes: Consequences for Redox and Electronic Spectroscopic Properties. J. Chem. Soc., Dalton Trans. 2002, 275−288. (o) Zell, T.; Milstein, D. Hydrogenation and Dehydrogenation Iron Pincer Catalysts Capable of Metal−Ligand Cooperation by Aromatization/Dearomatization. Acc. Chem. Res. 2015, 48, 1979−1994. (3) Forum issue on redox-active ligand, see: Chirik, P. J.; et al. Inorg. Chem. 2011, 50, 9737−9914. (4) Forum issue on redox-active ligand, see: Kaim, W.; et al. Eur. J. Inorg. Chem. 2012, 3, 340−580. (5) (a) Schrauzer, G. N.; Mayweg, V. P. Preparation, Reactions, and Structure of Bisdithio-α-diketone Complexes of Nickel, Palladium, and Platinum. J. Am. Chem. Soc. 1965, 87, 1483−1489. (b) Olson, D. C.; Mayweg, V. P.; Schrauzer, G. N. Polarographic Study of Coordination Compounds with Delocalized Ground States. Substituent Effects in Bis- and Trisdithiodiketone Complexes of Transition Metals. J. Am. Chem. Soc. 1966, 88, 4876−4882. (c) Schrauzer, G. N.; Rabinowitz, H. N. Charge Distribution and Nucleophilic Reactivity in Sulfur Ligand Chelates. Dialkyl Derivatives of Nickel(II), Palladium(II), and Platinum(II) Bis(cis ethylenedithiolates). J. Am. Chem. Soc. 1968, 90, 4297−4302. (d) Schrauzer, G. N.; Ho, R. K. Y.; Murillo, R. P. Structure, Alkylation, and Macrocyclic Derivatives of Bicyclo[2.2.1]hepta-2,5-diene Adducts of Metal Dithienes. J. Am. Chem. Soc. 1970, 92, 3508−3509. (6) (a) Schmitt, R. D.; Wing, R. M.; Maki, A. H. Donor-acceptor Complexes of the Inorganic.pi. Acceptor, Bis-cis-(1,2-perfluoromethylethene-1,2-dithiolato)nickel. J. Am. Chem. Soc. 1969, 91, 4394−4401. (b) Wing, R. M.; Tustin, G. C.; Okamura, W. H. Oxidative Cycloaddition of Metal Dithiolenes to Olefins. Synthesis and Characterization of Norbornadiene-bis-cis-(1,2-perfluoromethylethene-1,2-dithiolato)nickel. J. Am. Chem. Soc. 1970, 92, 1935−1939. (c) Baker, J. R.; Hermann, A.; Wing, R. M. Mechanism of Oxidative Cycloaddition of Olefins to Metal Dithiolenes. J. Am. Chem. Soc. 1971, 93, 6486−6489. (7) For recent examples of bis(dithiolene), see: (a) Kogut, E.; Tang, J. A.; Lough, A. J.; Widdifield, C. M.; Schurko, R. W.; Fekl, U. Neutral High-Potential Nickel Triad Bisdithiolenes: Structure and Solid-State NMR Properties of Pt[S2C2(CF3)2]2. Inorg. Chem. 2006, 45, 8850− 8852. (b) Harrison, D. J.; Nguyen, N.; Lough, A. J.; Fekl, U. New Insight into Reactions of Ni[S2C2(CF3)2]2 with Simple Alkenes: Alkene Adduct versus Dihydrodithiin Product Selectivity Is Controlled by [Ni[S2C2(CF3)2]2]− Anion. J. Am. Chem. Soc. 2006, 128, 11026− 11027. (c) Kerr, M. J.; Harrison, D. J.; Lough, A. J.; Fekl, U. LigandBased Reactivity of a Platinum Bisdithiolene: Double Diene Addition Yields a New C2-Chiral Chelate Ligand. Inorg. Chem. 2009, 48, 9043− 9045. (d) Mogesa, B.; Perera, E.; Rhoda, H. M.; Gibson, J. K.; Oomens, J.; Berden, G.; van Stipdonk, M. J.; Nemykin, V. N.; Basu, P. Solution, Solid, and Gas Phase Studies on a Nickel Dithiolene System: Spectator Metal and Reactor Ligand. Inorg. Chem. 2015, 54, 7703− 7716. (e) Mebrouk, K.; Kaddour, W.; Auban-Senzier, P.; Pasquier, C.; Jeannin, O.; Camerel, F.; Fourmigué, M. Molecular Alloys of Neutral Gold/Nickel Dithiolene Complexes in Single-Component Semiconductors. Inorg. Chem. 2015, 54, 7454−7460. (f) Filatre-Furcate, A.; Auban-Senzier, P.; Fourmigué, M.; Roisnel, T.; Dorcet, V.; Lorcy, D. Gold Dithiolene Complexes: Easy Access to 2-alkylthiothiazoledithiolate Complexes. Dalton Trans. 2015, 44, 15683−15689.

(g) Mebrouk, K.; Chotard, F.; Goff-Gaillard, C. L.; Arlot-Bonnemains, Y.; Fourmigué, M.; Camerel, F. Water-soluble Nickel-Bis(dithiolene) Complexes as Photothermal Agents. Chem. Commun. 2015, 51, 5268− 5270. (h) Chandrasekaran, P.; Greene, A. F.; Lillich, K.; Capone, S.; Mague, J. T.; DeBeer, S.; Donahue, J. P. A Structural and Spectroscopic Investigation of Octahedral Platinum Bis(dithiolene)phosphine Complexes: Platinum Dithiolene Internal Redox Chemistry Induced by Phosphine Association. Inorg. Chem. 2014, 53, 9192−9205. (i) Filatre-Furcate, A.; Bellec, N.; Jeannin, O.; Auban-Senzier, P.; Fourmigué, M.; Vacher, A.; Lorcy, D. Radical or Not Radical: Compared Structures of Metal (M = Ni, Au) Bis-Dithiolene Complexes with a Thiazole Backbone. Inorg. Chem. 2014, 53, 8681− 8690. (j) Ray, K.; Weyhermüller, T.; Neese, F.; Wieghardt, K. Electronic Structure of Square Planar Bis(benzene-1,2-dithiolato)metal Complexes [M(L)2]z (z = 2−, 1−, 0; M = Ni, Pd, Pt, Cu, Au): An Experimental, Density Functional, and Correlated ab Initio Study. Inorg. Chem. 2005, 44, 5345−5360. (k) Ray, S.; DeBeer George, S. D.; Solomon, E.; Wieghardt, K.; Neese, F. Description of the GroundState Covalencies of the Bis(dithiolato) Transition-Metal Complexes from X-ray Absorption Spectroscopy and Time-Dependent DensityFunctional Calculations. Chem. - Eur. J. 2007, 13, 2783−2797. (l) Ha, Y.; Tenderholt, A. L.; Holm, R. H.; Hedman, B.; Hodgson, K. O.; Solomon, E. I. Sulfur K-Edge X-ray Absorption Spectroscopy and Density Functional Theory Calculations on Monooxo MoIV and Bisoxo MoVI Bis-dithiolenes: Insights into the Mechanism of Oxo Transfer in Sulfite Oxidase and Its Relation to the Mechanism of DMSO Reductase. J. Am. Chem. Soc. 2014, 136, 9094−9105. (m) Liu, X.; Hou, G.-L.; Wang, X.; Wang, X.-B. Negative Ion Photoelectron Spectroscopy Reveals Remarkable Noninnocence of Ligands in Nickel Bis(dithiolene) Complexes [Ni(dddt)2]− and [Ni(edo)2]−. J. Phys. Chem. A 2016, 120, 2854−2862. (n) Kennedy, S. R.; Goyal, P.; Kozar, M. N.; Yennawar, H. P.; Hammes-Schiffer, S.; Lear, B. J. Effect of Protonation upon Electronic Coupling in the Mixed Valence and Mixed Protonated Complex, [Ni(2,3-pyrazinedithiol)2]. Inorg. Chem. 2016, 55, 1433−1445. (8) For recent computational examples of dithiolenes, see: (a) Fan, Y.; Hall, M. B. How Electron Flow Controls the Thermochemistry of the Addition of Olefins to Nickel Dithiolenes: Predictions by Density Functional Theory. J. Am. Chem. Soc. 2002, 124, 12076−12077. (b) Szilagyi, R. K.; Lim, B. S.; Glaser, T.; Holm, R. H.; Hedman, B.; Hodgson, K. O.; Solomon, E. I. Description of the Ground State Wave Functions of Ni Dithiolenes Using Sulfur K-edge X-ray Absorption Spectroscopy. J. Am. Chem. Soc. 2003, 125, 9158−9169. (c) Dang, L.; Yang, X.; Brothers, E. N.; Hall, M. B.; Zhou, J. Computational Studies on Ethylene Addition to Nickel Bis(dithiolene). J. Phys. Chem. A 2012, 116, 476−482. (d) Dang, L.; Shibl, M. F.; Yang, X.; Alak, A.; Harrison, D. J.; Fekl, U.; Brothers, E. N.; Hall, M. B. The Mechanism of Alkene Addition to a Nickel Bis(dithiolene) Complex: The Role of the Reduced Metal Complex. J. Am. Chem. Soc. 2012, 134, 4481−4484. (e) Dang, L.; Shibl, M. F.; Yang, X.; Harrison, D. J.; Alak, A.; Lough, A. J.; Fekl, U.; Brothers, E. N.; Hall, M. B. Apparent Anti-Woodward− Hoffmann Addition to a Nickel Bis(dithiolene) Complex: The Reaction Mechanism Involves Reduced, Dimetallic Intermediates. Inorg. Chem. 2013, 52, 3711−3723. (f) Dang, L.; Ni, S. F.; Hall, M. B.; Brothers, E. N. Uptake of One and Two Molecules of 1,3-Butadiene by Platinum Bis(dithiolene): A Theoretical Study. Inorg. Chem. 2014, 53, 9692−9702. (g) Li, H.; Brothers, E. N.; Hall, M. B. Computational Exploration of Alternative Catalysts for Olefin Purification: Cobalt and Copper Analogues Inspired by Nickel Bis(dithiolene) Electrocatalysis. Inorg. Chem. 2014, 53, 9679−9691. (h) Bushnell, E. A. C.; Boyd, R. J. Assessment of Several DFT Functionals in Calculation of the Reduction Potentials for Ni−, Pd−, and Pt−Bis-ethylene-1,2dithiolene and -Diselenolene Complexes. J. Phys. Chem. A 2015, 119, 911−918. (i) Yang, X.; Hall, M. B. Mechanism of Water Splitting and Oxygen-Oxygen Bond Formation by a Mononuclear Ruthenium Complex. J. Am. Chem. Soc. 2010, 132, 120−130. (j) Li, H.; Hall, M. B. Role of the Chemically Non-Innocent Ligand in the Catalytic Formation of Hydrogen and Carbon Dioxide from Methanol and I

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry Water with the Metal as the Spectator. J. Am. Chem. Soc. 2015, 137, 12330−12342. (9) For recent examples of tri(dithiolene), see: (a) Harrison, D. J.; Lough, A. J.; Nguyen, N.; Fekl, U. Push−Pull Molybdenum Trisdithiolenes Allow Rapid Nonconventional Binding of Ethylene at Ligand Sulfur Atoms. Angew. Chem., Int. Ed. 2007, 46, 7644−7647. (b) Tenderholt, A. L.; Szilagyi, R. K.; Holm, R. M.; Hodgson, K. O.; Hedman, B.; Solomon, E. I. Electronic Control of the “Bailar Twist” in Formally d0-d2 Molybdenum Tris(dithiolene) Complexes: A Sulfur Kedge X-ray Absorption Spectroscopy and Density Functional Theory Study. Inorg. Chem. 2008, 47, 6382−6392. (c) Harrison, D. J.; Fekl, U. Catalytic Production of Sulfur Heterocycles (dihydrobenzodithiins): A New Application of Ligand-based Alkene Reactivity. Chem. Commun. 2009, 7572−7574. (d) Morsing, T. J.; MacMillan, S. N.; Uebler, J. W. H.; Brock-Nannestad, T.; Bendix, J.; Lancaster, K. M. Stabilizing Coordinated Radicals via Metal−Ligand Covalency: A Structural, Spectroscopic, and Theoretical Investigation of Group 9 Tris(dithiolene) Complexes. Inorg. Chem. 2015, 54, 3660−3669. (e) Nguyen, N.; Harrison, D. J.; Lough, A. J.; De Crisci, A. G. D.; Fekl, U. Molybdenum Dithiolene Complexes as Structural Models for the Active Sites of Molybdenum(IV) Sulfide Hydrodesulfurization Catalysts. Eur. J. Inorg. Chem. 2010, 2010, 3577−3585. (f) Fekl, U.; Sarkar, B.; Kaim, W.; Zimmer-De Iuliis, M.; Nguyen, N. Tuning of the Spin Distribution between Ligand- and Metal-Based Spin: Electron Paramagnetic Resonance of Mixed-Ligand Molybdenum Tris(dithiolene) Complex Anions. Inorg. Chem. 2011, 50, 8685−8687. (g) Nguyen, N.; Lough, A. J.; Fekl, U. Rapid, Covalent Addition of Phosphine to Dithiolene in a Molybdenum Tris(dithiolene). A New Structural Model for Dimethyl Sulfoxide Reductase. Inorg. Chem. 2012, 51, 6446−6448. (h) Sproules, S.; Weyhermuller, T.; Goddard, R.; Wieghardt, K. The Rhenium Tris(dithiolene) Electron Transfer Series: Calibrating Covalency. Inorg. Chem. 2011, 50, 12623−12631. (i) Sproules, S.; Banerjee, P.; Weyhermuller, T.; Yan, Y.; Donahue, J. P.; Wieghardt, K. Monoanionic Molybdenum and Tungsten Tris(dithiolene) Complexes: A Multifrequency EPR Study. Inorg. Chem. 2011, 50, 7106−7122. (j) Yang, T.-L.; Ni, S.-F.; Zhang, P.; Dang, L. Ligand Effect on the Reactivity Difference of Mo Tris(dithiolene) Complexes Towards Ethylene: A Computational Study. J. Organomet. Chem. 2016, 806, 60−67. (10) Wang, K.; Stiefel, E. I. Toward Separation and Purification of Olefins Using Dithiolene Complexes: An Electrochemical Approach. Science 2001, 291, 106−109. (11) (a) Das, A.; Han, Z.; Brennessel, W.; Holland, P.; Eisenberg, R. Nickel Complexes for Robust Light-Driven and Electrocatalytic Hydrogen Production from Water. ACS Catal. 2015, 5, 1397−1406. (b) Zhang, J.; Du, P.; Schneider, J.; Jarosz, P.; Eisenberg, R. Photogeneration of Hydrogen from Water Using an Integrated System Based on TiO2 and Platinum(II) Diimine Dithiolate Sensitizers. J. Am. Chem. Soc. 2007, 129, 7726−7727. (c) Han, Z.; Eisenberg, R. Fuel from Water: The Photochemical Generation of Hydrogen from Water. Acc. Chem. Res. 2014, 47, 2537−2544. (d) Virca, C. N.; McCormick, T. M. DFT Analysis into the Intermediates of Nickel Pyridinethiolate Catalysed Proton Reduction. Dalton Trans. 2015, 44, 14333−14340. (e) McNamara, W. R.; Han, Z.; Yin, C.-J.; Brennessel, W. W.; Holland, P. L.; Eisenberg, R. Cobaltdithiolene Complexes for the Photocatalytic and Electrocatalytic Reduction of Protons in Aqueous Solutions. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 15594−15599. (f) McNamara, W. R.; Han, Z.; Alperin, P. J.; Brennessel, W. W.; Holland, P. L.; Eisenberg, R. A Cobalt− Dithiolene Complex for the Photocatalytic and Electrocatalytic Reduction of Protons. J. Am. Chem. Soc. 2011, 133, 15368−15371. (g) Das, A.; Han, Z.; Haghighi, M. G.; Eisenberg, R. Photogeneration of Hydrogen from Water Using CdSe Nanocrystals Demonstrating the Importance of Surface Exchange. Proc. Natl. Acad. Sci. U. S. A. 2013, 110, 16716−16723. (h) Downes, C. A.; Marinescu, S. C. Efficient Electrochemical and Photoelectrochemical H2 Production from Water by a Cobalt Dithiolene One-Dimensional Metal−Organic Surface. J. Am. Chem. Soc. 2015, 137, 13740−13743. (i) Solis, B. H.; HammesSchiffer, S. Computational Study of Anomalous Reduction Potentials

for Hydrogen Evolution Catalyzed by Cobalt Dithiolene Complexes. J. Am. Chem. Soc. 2012, 134, 15253−15256. (j) Letko, C. S.; Panetier, J. A.; Head-Gordon, M.; Tilley, T. D. Mechanism of the Electrocatalytic Reduction of Protons with Diaryldithiolene Cobalt Complexes. J. Am. Chem. Soc. 2014, 136, 9364−9376. (k) Panetier, J. A.; Letko, C. S.; Tilley, T. D.; Head-Gordon, M. Computational Characterization of Redox Non-Innocence in Cobalt-Bis(Diaryldithiolene)-Catalyzed Proton Reduction. J. Chem. Theory Comput. 2016, 12, 223−230. (12) (a) Zarkadoulas, A.; Field, M. J.; Papatriantafyllopoulou, C.; Fize, J.; Artero, V.; Mitsopoulou, C. A. Experimental and Theoretical Insight into Electrocatalytic Hydrogen Evolution with Nickel Bis(aryldithiolene) Complexes as Catalysts. Inorg. Chem. 2016, 55, 432− 444. (b) Gan, L.; Groy, T. L.; Tarakeshwar, P.; Mazinani, S. K. S.; Shearer, J.; Mujica, V.; Jones, A. K. A Nickel Phosphine Complex as a Fast and Efficient Hydrogen Production Catalyst. J. Am. Chem. Soc. 2015, 137, 1109−1115. (13) (a) Dilworth, J. R.; Hutson, A. J.; Morton, S.; Harman, M.; Hursthouse, M. B.; Zubieta, J.; Archer, C. M.; Kelly, J. D. The Preparation and Electrochemistry of Technetium and Rhenium Complexes of 2-(diphenylphosphino)benzenethiol. The Crystal and Molecular Structures of [Re(2-Ph 2 PC 6 H 4 S) 3 ] and [Tc(2Ph2PC6H4S)3]. Polyhedron 1992, 11, 2151−2155. (b) Dilworth, J. R.; Zheng, Y.; Lu, S.; Wu, Q. Preparation and Characterization of A Novel Asymmetrically Oxidized Complex of 2-(diphenylphosphino)benzenethiol with Ruthenium. The Crystal and Molecular Structure of [Ru(2-Ph 2 PC 6 H 4 S)·(2-Ph 2 PC 6 H 4 S−OH)(2-Ph 2 PC 6 H 4 SO 2 )]·1/ 2H2O. Transition Met. Chem. 1992, 17, 364−368. (14) (a) Grapperhaus, C. A.; Poturovic, S.; Mashuta, M. S. Dichloromethane Alkylates a Trithiolato-Ruthenium Complex to Yield a Methylene-Bridged Thioether Core. Synthesis and Structural Comparison to the Thiolato-Ruthenium Precursor. Inorg. Chem. 2002, 41, 4309−4311. (b) Grapperhaus, C. A.; Poturovic, S. Electrochemical Investigations of the [Tris(2-(diphenylphosphino)thiaphenolato)ruthenate(II)] Monoanion Reveal Metal- and Ligand-Centered Events: Radical, Reactivity, and Rate. Inorg. Chem. 2004, 43, 3292− 3298. (c) Ouch, K.; Mashuta, M. S.; Grapperhaus, C. A. MetalStabilized Thiyl Radicals as Scaffolds for Reversible Alkene Addition via C−S Bond Formation/Cleavage. Inorg. Chem. 2011, 50, 9904− 9914. (d) Poturovic, S.; Grapperhaus, C. A.; Mashuta, M. S. Carbon− Sulfur Bond Formation between a Ruthenium-Coordinated Thiyl Radical and Methyl Ketones. Angew. Chem., Int. Ed. 2005, 44, 1883− 1887. (e) Grapperhaus, C. A.; Venna, K. B.; Mashuta, M. S. Carbon− Sulfur Bond Formation via Alkene Addition to an Oxidized Ruthenium Thiolate. Inorg. Chem. 2007, 46, 8044−8050. (f) Ouch, K.; Mashuta, M. S.; Grapperhaus, C. A. Alkyne Addition to a MetalStabilized Thiyl Radical: Carbon−Sulfur Bond Formation between 1Octyne and [Ru(SP)3]+. Eur. J. Inorg. Chem. 2012, 2012, 475−478. (g) Sampson, K. O.; Kumar, D.; Mashuta, M. S.; Grapperhaus, C. A. Addition of Polysubstituted Alkenes, Aromatic alkynes, and Dienes to a Metal-stabilized Thiyl Radical via Carbon−Sulfur Bond Formation: Electrochemical, Chemical, and Computational Investigations. Inorg. Chim. Acta 2013, 408, 1−8. (h) Grapperhaus, C. A.; Kozlowski, P. M.; Kumar, D.; Frye, H. N.; Venna, K. B.; Poturovic, S. Singlet Diradical Character of an Oxidized Ruthenium Trithiolate: Electronic Structure and Reactivity. Angew. Chem., Int. Ed. 2007, 46, 4085−4088. (i) Lu, M.; Campbell, J. L.; Chauhan, R.; Grapperhaus, C. A.; Chen, H. Probing the Reactivity and Radical Nature of Oxidized Transition Metal-Thiolate Complexes by Mass Spectrometry. J. Am. Soc. Mass Spectrom. 2013, 24, 502−512. (15) (a) Grapperhaus, C. A.; Ouch, K.; Mashuta, M. S. RedoxRegulated Ethylene Binding to a Rhenium-Thiolate Complex. J. Am. Chem. Soc. 2009, 131, 64−65. (b) Chauhan, R.; Moreno, M.; Banda, D. M.; Zamborini, F. P.; Grapperhaus, C. A. Chemiresistive Metalstabilized Thiyl Radical Films as Highly Selective Ethylene Sensors. RSC Adv. 2014, 4, 46787−46790. (c) Haddad, A. Z.; Kumar, D.; Ouch Sampson, K.; Matzner, A. M.; Mashuta, M. S.; Grapperhaus, C. A. Proposed Ligand-Centered Electrocatalytic Hydrogen Evolution and Hydrogen Oxidation at a Noninnocent Mononuclear Metal−Thiolate. J. Am. Chem. Soc. 2015, 137, 9238−9241. J

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Cioslowski, J.; Fox, D. J. Gaussian 09, revision D. 01; Gaussian, Inc.: Wallingford, CT, 2013. (25) Chai, J.-D.; Head-Gordon, M. Long-range Corrected Hybrid Density Functionals with Damped Atom−atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615−6620. (26) (a) Andrae, D.; Haeussermann, U.; Dolg, M.; Stoll, H.; Preuss, H. Energy-adjustedab initio Pseudopotentials for the Second and Third Row Transition Elements. Theor. Chim. Acta 1990, 77, 123−141. (b) Martin, J. M. L.; Sundermann, A. Correlation Consistent Valence Basis Sets for Use with the Stuttgart−Dresden−Bonn Relativistic Effective Core Potentials: The atoms Ga−Kr and In−Xe. J. Chem. Phys. 2001, 114, 3408−3420. (27) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378−6396. (28) Roothaan, C. C. J. New Developments in Molecular Orbital Theory. Rev. Mod. Phys. 1951, 23, 69−89. (29) Werner, H. J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G.; Shamasundar, K. R.; Adler, T. B.; Amos, R. D.; Bernhardsson, A.; Berning, A.; Cooper, D. L.; Deegan, M. J. O.; Dobbyn, A. J.; Eckert, F.; Goll, E.; Hampel, C.; Hesselmann, A.; Hetzer, G.; Hrenar, T.; Jansen, G.; Köppl, C.; Liu, Y.; Lloyd, A. W.; Mata, R. A.; May, A. J.; McNicholas, S. J.; Meyer, W.; Mura, M. E.; Nicklass, A.; O’Neill, D. P.; Palmieri, P.; Pflüger, K.; Pitzer, R.; Reiher, M.; Shiozaki, T.; Stoll, H.; Stone, A. J.; Tarroni, R.; Thorsteinsson, T.; Wang, M.; Wolf, A. MOLPRO, version 2012, a package of ab initio programs; see http:// www.molpro.net.

(16) Tang, H.; Guan, J.; Hall, M. B. Understanding the Radical Nature of an Oxidized Ruthenium Tris(thiolate) Complex and Its Role in the Chemistry. J. Am. Chem. Soc. 2015, 137, 15616−15619. (17) Roos, B. O. The Complete Active Space Self-Consistent Field Method and its Applications in Electronic Structure Calculations. Adv. Chem. Phys. 1987, 69, 399−445. (18) (a) Andersson, K.; Malmqvist, P.-A.; Roos, B. O. Second-order Perturbation Theory with a Complete Active Space Self-consistent Field Reference Function. J. Chem. Phys. 1992, 96, 1218−1226. (b) Finley, J.; Malmqvist, P.-A.; Roos, B. O.; Serrano-Andres, L. The Multi-state CASPT2Method. Chem. Phys. Lett. 1998, 288, 299−306. (19) (a) Head-Gordon, M.; Pople, J. A.; Frisch, M. J. MP2 Energy Evaluation by Direct Methods. Chem. Phys. Lett. 1988, 153, 503−506. (b) Frisch, M. J.; Head-Gordon, M.; Pople, J. A. A Direct MP2 Gradient Method. Chem. Phys. Lett. 1990, 166, 275−280. (c) Frisch, M. J.; Head-Gordon, M.; Pople, J. A. Semi-direct Algorithms for the MP2 Energy and Gradient. Chem. Phys. Lett. 1990, 166, 281−289. (d) Almlöf, J.; Faegri, K., Jr.; Korsell, K. Principles for a Direct SCF Approach to LICAO−MO ab-initio calculations. J. Comput. Chem. 1982, 3, 385−399. (e) Head-Gordon, M.; Head-Gordon, T. Analytic MP2 Frequencies without Fifth-order Storage. Theory and Application to Bifurcated Hydrogen Bonds in the Water Hexamer. Chem. Phys. Lett. 1994, 220, 122−128. (20) Pople, J. A.; Seeger, R.; Krishnan, R. Variational Configuration Interaction Methods and Comparison with Perturbation Theory. Int. J. Quantum Chem. 1977, 12, 149−163. (21) (a) Krishnan, R.; Pople, J. A. Approximate Fourth-order Perturbation Theory of the Electron Correlation Energy. Int. J. Quantum Chem. 1978, 14, 91−100. (b) Trucks, G. W.; Salter, E. A.; Sosa, C.; Bartlett, R. J. Theory and Implementation of the MBPT Density Matrix. An Application to One-electron Properties. Chem. Phys. Lett. 1988, 147, 359−366. (c) Trucks, G. W.; Watts, J. D.; Salter, E. A.; Bartlett, R. J. Analytical MBPT(4) Gradients. Chem. Phys. Lett. 1988, 153, 490−495. (22) (a) Cizek, J. On the Use of the Cluster Expansion and the Technique of Diagrams in Calculations of Correlation Effects in Atoms and Molecules. Adv. Chem. Phys. 1969, 14, 35−89. (b) Purvis, G. D.; Bartlett, R. J. A Full Coupled-Cluster Singles and Doubles Model: The Inclusion of Disconnected Triples. J. Chem. Phys. 1982, 76, 1910−1918. (c) Scuseria, G. E.; Janssen, C. L.; Schaefer, H. F., III. An Efficient Reformulation of the Closed-shell Coupled Cluster Single and Double Excitation (CCSD) Equations. J. Chem. Phys. 1988, 89, 7382−7387. (d) Scuseria, G. E.; Schaefer, H. F., III. Is Coupled Cluster Singles and Doubles (CCSD) More Computationally Intensive Than Quadratic Configuration Interaction (QCISD)? J. Chem. Phys. 1989, 90, 3700−3703. (e) Pople, J. A.; Head-Gordon, M.; Raghavachari, K. Quadratic Configuration Interaction. A General Technique for Determining Electron Correlation Energies. J. Chem. Phys. 1987, 87, 5968−5975. (23) (a) Krishnan, R.; Schlegel, H. B.; Pople, J. A. Derivative Studies in Configuration Interaction Theory. J. Chem. Phys. 1980, 72, 4654− 4655. (b) Raghavachari, K.; Pople, J. A. Calculation of One-electron Properties Using Limited Configuration Interaction Techniques. Int. J. Quantum Chem. 1981, 20, 1067−1071. (24) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, O.; Foresman, J. B.; Ortiz, J. V.; K

DOI: 10.1021/acs.inorgchem.6b02434 Inorg. Chem. XXXX, XXX, XXX−XXX