Raphael M.
Jay
*ab,
Sebastian
Eckert
ab,
Mattis
Fondell
b,
Piter S.
Miedema†
b,
Jesper
Norell
c,
Annette
Pietzsch
b,
Wilson
Quevedo
b,
Johannes
Niskanen‡
b,
Kristjan
Kunnus
d and
Alexander
Föhlisch
ab
aUniversität Potsdam, Institut für Physik und Astronomie, 14476 Potsdam, Germany. E-mail: rajay@uni-potsdam.de
bHelmholtz-Zentrum Berlin für Materialien und Energie GmbH, Institute for Methods and Instrumentation for Synchrotron Radiation Research, 12489 Berlin, Germany
cStockholm University, Department of Physics, Albanova University Center, 10691 Stockholm, Sweden
dPULSE Institute, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
First published on 13th September 2018
Understanding and controlling properties of transition metal complexes is a crucial step towards tailoring materials for sustainable energy applications. In a systematic approach, we use resonant inelastic X-ray scattering to study the influence of ligand substitution on the valence electronic structure around an aqueous iron(II) center. Exchanging cyanide with 2-2′-bipyridine ligands reshapes frontier orbitals in a way that reduces metal 3d charge delocalization onto the ligands. This net decrease of metal–ligand covalency results in lower metal-centered excited state energies in agreement with previously reported excited state dynamics. Furthermore, traces of solvent-effects were found indicating a varying interaction strength of the solvent with ligands of different character. Our results demonstrate how ligand exchange can be exploited to shape frontier orbitals of transition metal complexes in solution-phase chemistry; insights upon which future efforts can built when tailoring the functionality of photoactive systems for light-harvesting applications.
In a recent series of studies, Gaffney and coworkers demonstrated how exchanging and combining different kinds of ligands allows to manipulate rates and pathways of MLCT state decay in Fe(II) complexes.8–11 By varying the ligands coordinated towards the metal center within a suitable solvent environment, the relative energies between MLCT and metal-centered (MC) states could be manipulated in a way that intersystem crossing channels were quenched. Specifically, it was found that the MLCT lifetime can be increased by a factor of 100 compared to [Fe(bpy)3]2+ by substituting two 2-2′-bipyridine (bpy) ligands with four cyanide (CN) ligands and dissolving the compound in dimethyl sulfoxide.9 By combining time-resolved Kβ X-ray emission and transient absorption spectroscopy, MLCT and MC states could be distinguished, providing a robust description of how molecular structure varies the dynamics following a charge-transfer excitation.
In this study, we demonstrate how these differences can be traced back to the varying nature of frontier orbitals induced by ligand substitution. Along the series [Fe(CN)6]4−, [Fe(bpy)(CN)4]2− and [Fe(bpy)3]2+, we systematically study changes in the ground state valence electronic structure and deduce trends in the valence excited state landscape. We use resonant inelastic X-ray scattering12 (RIXS) in the soft X-ray regime, that, due to its local and element-specific access to valence excited final states, is a probe particularly suited for such investigations. Through the Fe 2p → 3d core-excitation, the method is sensitive to occupied as well as unoccupied orbitals around the central metal,13–18 while the nitrogen (N) 1s → 2p core-excitation probes the analogue ligand orbitals.19–21 Thus, our approach allows to experimentally access frontier orbitals that contribute to the metal–ligand bond as well as participate in solute–solvent interaction. By independently investigating local charge densities and orbital character of Fe 3d as well as CN and bpy derived orbitals, trends in ligand field strength, metal–ligand bonding and solvation are observed that are directly linked to the varying excited-state landscape of charge-transfer Fe complexes with different ligand cages.
The quantum chemical simulations were for the most part performed within the same density functional theory (DFT) framework as by Van Kuiken et al.24 For the discussion of orbitals and charges, the molecular structure of each sample was optimized with the ORCA software25 at the B3LYP level26,27 using the all-electron def2-TZVP(-f) basis set28 for all atoms. All calculations have been additionally performed with the reparametrized B3LYP* functional;29 results for this functional showed a sightly worse agreement with the experimental data, but were in overall good agreement with ones presented here for B3LYP. To simulate the solvent-effects of water the conductor-like polarizable continuum model COSMO30 was employed for all calculations. Additionally, the RIJCOSX method31 was used with the def2-TZV/J as auxiliary basis set32 to speed up calculations. Based on the optimized structures, X-ray absorption and emission spectra at the N K-edge were obtained using time-dependent DFT. Core-excitations were simulated individually for N sites of different chemical environment by using the Pipek–Mezey orbital localization scheme.33 The X-ray absorption spectra were generated by convolution of each transition with a Gaussian function of 0.5 eV FWHM for bpy ligands and 0.9 eV FWHM for CN ligands. A shift of 12.3 eV was applied to the X-ray absorption spectrum of [Fe(bpy)3]2+ to match the experimental partial fluorescence yield. All other absorption spectra were shifted equally. For the X-ray emission spectra, all transitions were convolved with a Gaussian function of 1.5 eV FWHM. The spectra were shifted individually to match the energy transfer of the experimental data.
Fig. 1 (a) Molecular structure and (b) Fe L3/2-edge partial fluorescence yield absorption spectra of [Fe(CN)6]4−, [Fe(bpy)(CN)4]2− and [Fe(bpy)3]2+ in aqueous solution. |
Fig. 1b shows the Fe L3/2-edge partial fluorescence yield (PFY) of the three samples under investigation. The spectra are acquired by integrating over all detected L-edge X-ray emission channels and are normalized to the maximum of the first resonance. By comparison of the detected transitions, the implications of ligand exchange on the local valence electronic structure can be observed, exemplified by two major trends visible in the spectra. By gradually exchanging CN with bpy ligands, a drastic decrease in relative intensity of the second absorption feature at ∼710.5 eV incident energy can be observed. In [Fe(CN)6]4−, this feature has been identified to originate from a mixing of metal t2g MOs with the anti-bonding CN 2π* MOs of equal symmetry. Thereby, metal electron density is delocalized onto the ligands. As the Fe 2p → 2π* core-excitation probes the Fe content of the CN 2π* MO, the resulting feature has been postulated to be a “direct probe of back-bonding”.34 Consequently, the gradually lowered intensity of the feature in [Fe(bpy)(CN)4]2− and [Fe(bpy)3]2+ is reflective of the reduced ability of bpy ligands to accept metal t2g electron density. The remaining absorption in [Fe(bpy)3]2+ at 710 eV incident energy can therefore be rationalized as a combination of remaining Fe 2p → 2π* excitation and multiplet effects that previously have been observed in complexes with similar ligand cages.35–37 Additional multiplet-effects caused by the simultaneous presence of CN and bpy ligands and the resulting deviation from octahedrality could therefore explain the broadened second resonance in [Fe(bpy)(CN)4]2− compared to [Fe(CN)6]4−.
Besides the decrease of Fe 2p → 2π* intensity upon ligand-exchange, the onset of the lowest lying Fe 2p → 3d excitation at the L3-edge (corresponding to a Fe 2p → eg excitation) shifts to lower incident energy. Shifts of the L-edge absorption onset have traditionally been interpreted to reflect the local charge at the metal center, where an increase in electron density enhances the ability to screen the Fe 2p core-hole in the core-excited state.38 This has been proven to be a useful tool to determine the oxidation state of the central metal.39 As previously stated, exclusively Fe(II) complexes are compared within this study. The shift of the absorption onset is thus directly reflective of a change in covalent metal–ligand interaction. The shift to lower incident energy can therefore be interpreted as a reduced delocalization of metal electron density onto the ligands (i.e. backbonding) and subsequent increase of electron density at the metal center, in agreement with the decrease of Fe 2p → 2π* intensity.
For a discussion of the RIXS final states of the three complexes we turn to Fig. 2. The displayed range of incident energy is restricted to the L3-edge as the interpretation of the L2-edge is complicated by the presence of Coster–Kronig features.40 Similarly to the discussion of the L-edge PFY, the RIXS plane of [Fe(CN)6]4− displayed in (a) can serve as a reference, as a detailed assignment of final states has been performed by Kunnus et al.16 The two most intense features can be identified as t52ge1g (709.1 eV incident energy, 3.5 eV energy loss) and t52g2π*1 final states (710.8 eV incident energy, 5 eV energy loss). The decrease in intensity of the latter final state in [Fe(bpy)(CN)4]2− and almost complete disappearance in [Fe(bpy)3]2+ (compare (b) and (c)) is the cause for the trend of decreasing integrated Fe 2p → 2π* intensity observed in the PFYs.
The t52ge1g final state allows for a more quantitative analysis. In octahedral symmetry, the position of the feature in terms of energy loss effectively probes the difference between the t2g and eg levels and can therefore measure trends in 10Dq.41Fig. 2d shows intensities integrated over the respective marked incident energy ranges in the RIXS maps. The final state energy loss therefore decreases from 3.5 eV in [Fe(CN)6]4− to 3.1 eV in [Fe(bpy)(CN)4]2− and 2.2 eV in [Fe(bpy)3]2+. This trend is consistent with the field strength of the two different ligands predicted by the spectrochemical series.42 Furthermore, it agrees with the mechanistic explanation of the extended MLCT lifetime of [Fe(bpy)(CN)4]2− dissolved in dimethyl sulfoxide. There, spin crossover decay channels were successfully quenched by destabilizing MC states through an exchange of bpy with CN ligands, while simultaneously stabilizing MLCT energies by solvents of low Lewis-acidity.9,11
To substantiate the interpretations of the observed spectral trends we performed DFT and time-dependent (TD) DFT calculations. TD-DFT provides a natural link between experimental X-ray spectra and the ground state electronic structure of the complexes, as both the electronic ground state and the X-ray excitations are here described within the same basis of (Kohn–Sham) molecular orbitals. This approach further provides a framework in which the different complexes under investigation can be similarly described, besides their different ligand structures; in contrast to restricted active space approaches,16,43–45 where comparability is only guaranteed if the orbitals within the respective active spaces are the same for all complexes.
As previously stated, the Fe L-edge PFY measurements suggest reduced delocalization of electron density, i.e. a reduced π-backbonding, onto the ligands upon exchanging CN with bpy ligands. This trend is reproduced by our calculations, as can be seen in Fig. 3, where the highest occupied molecular orbitals (HOMOs) for [Fe(CN)6]4−, [Fe(bpy)(CN)4]2− and [Fe(bpy)3]2+ are displayed. In [Fe(CN)6]4− (a), there is significant amplitude visible at the N sites in plane with the Fe 3d MOs. This mixing of metal and ligand orbitals is reduced in [Fe(bpy)(CN)4]2− (b), with smaller contribution on the CN ligands and no contribution visible on the bpy ligand. Consequently, no amplitude can be observed on any of the bpy ligands in [Fe(bpy)3]2+ (c). This degree of mixing is quantified in Table 1 by the average Fe 3d character of the HOMO, HOMO(−1), HOMO(−2) orbitals (representing the Fe t2g orbitals in octahedral symmetry). The Fe 3d character increases as a result of exchanging CN with bpy ligands. This is in agreement with the gradually vanishing Fe 2p → 2π* resonance in the Fe L-edge PFY (compare Fig. 1b), as it measures the mixing of the Fe 3d and ligand 2π* orbitals.
HOMO Fe 3d character | Mulliken | Hirshfeld | Loewdin | |
---|---|---|---|---|
[Fe(CN)6]4− | 75.6 | −0.01 | −1.01 | −0.83 |
[Fe(bpy)(CN)4]2− | 78.0 | −0.15 | −1.12 | −0.62 |
[Fe(bpy)3]2+ | 83.7 | −0.23 | −1.14 | −0.26 |
The trend of reduced delocalization of metal electron density can also be expected to be reflected by a varying local charge of the Fe centers in different ligand cages. Table 1 therefore additionally shows the Mulliken, Hirshfeld and Loewdin charge analysis of the respective Fe sites in the three complexes. Consistent with our interpretation, the Mulliken and Hirshfeld charges are indeed turning more negative in the order [Fe(CN)6]4−, [Fe(bpy)(CN)4]2−, [Fe(bpy)3]2+ as a result of incremental concentration of electronic density at the metal site. This potentially causes the experimentally observed shift of the L-edge PFY absorption onset to lower incident energy (compare Fig. 1b) as a result of increased electronic screening of the Fe 2p core-hole. However, a recent study by Kubin et al.46 suggests that such a picture might be an oversimplification. The Loewdin charges also presented in Table 1, in fact, show the opposite trend. This exemplifies the known difficulty of unambiguously and accurately defining atomic charges from electronic structure calculations. Future in-depth theoretical studies, possibly based on higher level quantum chemical methods, may therefore be necessary to derive a mechanistic explanation also for the observed onset shifts resulting from changes in metal–ligand covalency, similarly as for varying oxidation states.
Fig. 4b shows RIXS spectra of the three compounds gathered by integration over their respective main resonance in the PFY (compare Fig. 4a). To facilitate some degree of comparability between spectra generated by integration over different energy ranges, they are normalized to their elastic peak maximum at 0 eV energy loss. For [Fe(CN)6]4− and [Fe(bpy)3]2+, non-resonant X-ray emission (XES) calculations are additionally shown. Intensity from elastic scattering (0 eV energy loss) is not accounted for by non-resonant XES calculations, where the core-excited electron is ignored. The calculations are normalized and energetically shifted to the experimental emission maximum. The main features of [Fe(CN)6]4− can now be assigned to dominantly result from inelastic scattering from the Fe t2g (4.6 eV energy loss) and CN 1π MOs (8.3 eV energy loss) in agreement with the final state assignment by Kunnus et al.16 The biggest difference between experiment and theory here is the underestimation of the emission signal from the Fe t2g orbitals at 4.6 eV energy loss. Similar to the broadening of the absorption resonance observed in the N K-edge PFY, this could be due to an insufficient treatment of the solvent environment. The earlier mentioned withdrawal of electron density from the CN ligand in an aqueous environment has been reported to subsequently increase π-backdonation in Fe cyanides.11,47,48 Due to its sensitivity to the overlap of the N 1s core-hole with the Fe t2g orbitals, this emission line is a direct probe thereof. Within a simple model similar to Kjær et al.,11 we simulated the solute–solvent interaction by explicitly adding one water molecule to each N site of the six CN ligands. This calculation is additionally shown and indeed exhibits an increase of emission from the t2g orbitals, however to a very small degree. While strict solvent-treatment can therefore be expected to increase accuracy, the remaining disagreement illustrates the limits of DFT to simultaneously and with high accuracy account for metal-, ligand- and solvent structure within the same framework of basis sets and exchange correlation potentials.
When comparing the RIXS of [Fe(CN)6]4− with [Fe(bpy)3]2+, the latter one exhibits significantly more features, directly reflective of the more complicated bpy ring structure. This is confirmed by our calculations, where emission from a variety of orbitals with predominant π character can be assigned to the dominant spectral features. Within the applied normalization the individual emission lines have less intensity than the main feature in [Fe(CN)6]4−. This can be explained by the fact that the bpy π* orbitals are distributed across the pyridine rings and therefore have less overlap with the N core-hole than in a CN ligand.
Having established the characteristic RIXS signatures of the N sites in a CN and a bpy ligand, the signature of [Fe(bpy)(CN)4]2− can be understood as well. Analogously to the N 1s PFYs, the RIXS spectrum of [Fe(bpy)(CN)4]2− is a superposition of the different N sites. By adding the spectra of [Fe(CN)6]4− and [Fe(bpy)3]2+ weighted by the ratio of N sites in CN and bpy ligands (2:1), the spectrum of [Fe(bpy)(CN)4]2− can be reproduced with high accuracy, as can be seen by the gray spectrum superimposed onto the measured spectrum. Different ligands in [Fe(bpy)(CN)4]2− are thus seen as independent moieties from a N 1s spectroscopic perspective. Furthermore, the excellent agreement also justifies the applied normalization on the elastic line.
Footnotes |
† Present address: Deutsches Elektronen-Synchrotron DESY, Photon Science, 22607 Hamburg, Germany. |
‡ Present address: University of Turku, Department of Physics and Astronomy, FI-20014, Finland. |
This journal is © the Owner Societies 2018 |