One- and two-photon absorption spectra of organoboron complexes: vibronic and environmental effects.

We synthesized a series of four parent aza-ß-ketoiminate organoboron complexes and performed spectroscopic studies using both experimental and computational techniques. We studied how benzannulation influences the vibronic structure of the UV/Vis absorption bands with a focus on the bright lowest-energy π → π* electronic excitation. Theoretical simulations, accounting for inhomogeneous broadening effects using different embedding schemes, allowed gaining in-depth insights into the observed differences in band shapes induced by structural modifications. We observed huge variations in the distributions of vibronic transitions depending on the position of benzannulation. By and large, the harmonic approximation combined with the adiabatic hessian model delivers qualitatively correct band shapes for the one-photon absorption spectra, except in one case. We also assessed the importance of non-Condon effects (accounted for by the linear term in Herzberg-Teller expansion of the dipole moment) for S0 → S1 band shapes. It turned out that non-Condon contributions have no effect on the band shape in one-photon absorption spectra. In contrast, these effects significantly change the Franck-Condon band shapes of the two-photon absorption spectra. For one of the studied organoboron complexes we also performed a preliminary exploration of mechanical anharmonicity, resulting in an increase of the intensity of the 0-0 transition, which improves the agreement with the experimental data compared to the harmonic model.

Ti(III) Catalysts for CO2/Epoxide Copolymerization at Unusual Ambient Pressure Conditions.

Titanium compounds in low oxidation states are highly reducing species and hence powerful tools for the functionalization of small molecules. However, their potential has not yet been fully realized because harnessing these highly reactive complexes for productive reactivity is generally challenging. Advancing this field, herein we provide a detailed route for the formation of titanium(III) orthophenylendiamido (PDA) species using [LiBHEt3] as a reducing agent. Initially, the corresponding lithium PDA compounds [Li2(ArPDA)(thf)3] (Ar = 2,4,6-trimethylphenyl (MesPDA), 2,6-diisopropylphenyl (iPrPDA)) are combined with [TiCl4(thf)2] to form the heterobimetallic complexes [{TiCl(ArPDA)}(µ-ArPDA){Li(thf)n}] (n = 1, Ar = iPr 3 and n = 2, Ar = Mes 4). Compound 4 evolves to species [Ti(MesPDA)2] (6) via thermal treatment. In contrast, the transformation of 3 into [Ti(iPrPDA)2] (5) only occurs in the presence of [LiNMe2], through a lithium-assisted process, as revealed by density functional theory (DFT). Finally, the Ti(IV) compounds 3-6 react with [LiBHEt3] to give rise to the Ti(III) species [Li(thf)4][Ti(ArPDA)2] (Ar = iPr 8, Mes 9). These low-valent compounds in combination with [PPN]Cl (PPN = bis(triphenylphosphine)iminium) are proved to be highly selective catalysts for the copolymerization of CO2 and cyclohexene epoxide. Reactions occur at 1 bar pressure with activity/selectivity levels similar to Salen-Cr(III) compounds.

A new computational tool for interpreting the infrared spectra of molecular complexes.

The popularity of infrared (IR) spectroscopy is due to its high interpretive power. This study presents a new computational tool for analyzing the IR spectra of molecular complexes in terms of intermolecular interaction energy components. In particular, the proposed scheme enables associating the changes in the IR spectra occurring upon complex formation with individual types of intermolecular interactions (electrostatic, exchange, induction, and dispersion), thus providing a completely new insight into the relations between the spectral features and the nature of interactions in molecular complexes. To demonstrate its interpretive power, we analyze, for selected vibrational modes, which interaction types rule the IR intensity changes upon the formation of two different types of complexes, namely π⋯π stacked (benzene⋯1,3,5-trifluorobenzene) and hydrogen-bonded (HCN⋯HNC) systems. The exemplary applications of the new scheme to these two molecular complexes revealed that the interplay of interaction energy components governing their stability might be very different from that behind the IR intensity changes. For example, in the case of the dispersion-bound π⋯π-type complex, dispersion contributions to the interaction induced IR intensity of the selected modes are notably smaller than their first-order (electrostatic and exchange) counterparts.

Pitfall in simulations of vibronic TD-DFT spectra: diagnosis and assessment.

In this Communication, we study the effect of spurious oscillations in the profiles of energy derivatives with respect to nuclear coordinates calculated with density functional approximations (DFAs) for formaldehyde, pyridine, and furan in their ground and electronic excited states. These spurious oscillations, which can only be removed using extensive integration grids that increase enormously the CPU cost of DFA calculations, are significant in the case of third- and fourth-order energy derivatives of the ground and excited states computed by M06-2X and ωB97X functionals. The errors in question propagate to anharmonic vibronic spectra computed under the Franck-Condon approximation, i.e., positions and intensities of vibronic transitions are affected to a large extent (shifts as significant as hundreds of cm-1 were observed). On the other hand, the LC-BLYP and CAM-B3LYP functionals show a much less pronounced effect due to spurious oscillations. Based on the results presented herein, we recommend either LC-BLYP or CAM-B3LYP with integration grids (250, 974) (or larger) for numerically stable simulations of vibronic spectra including anharmonic effects.

Decoding the infrared spectra changes upon formation of molecular complexes: the case of halogen bonding in pyridine⋯perfluorohaloarene complexes.

A recently developed computational scheme is employed to interpret changes in the infrared spectra of halogen-bonded systems in terms of intermolecular interaction energy components (electrostatic, exchange, induction, dispersion) taking pyridine⋯perfluorohaloarene complexes as examples. For all complexes, we find a strong linear correlation between the different terms of the interaction-induced changes of the IR band associated with an intermolecular halogen bond stretching mode and the corresponding terms of the interaction energy, which implies that the interaction components play similar roles in both properties. This is not true for other vibrational modes localized in one of the monomers studied here, for which the corresponding interaction-induced changes in IR bands may present a completely different decomposition than the interaction energy.

Effect of external electric fields in the charge transfer rates of donor-acceptor dyads: A straightforward computational evaluation.

We present a straightforward and low-cost computational protocol to estimate the variation of the charge transfer rate constant, kCT, in a molecular donor-acceptor caused by an external electric field. The proposed protocol also allows for determining the strength and direction of the field that maximize the kCT. The application of this external electric field results in up to a >4000-fold increase in the kCT for one of the systems studied. Our method allows the identification of field-induced charge-transfer processes that would not occur without the perturbation caused by an external electric field. In addition, the proposed protocol can be used to predict the effect on the kCT due to the presence of charged functional groups, which may allow for the rational design of more efficient donor-acceptor dyads.

Characterization of a Ferryl Flip in Electronically Tuned Nonheme Complexes. Consequences in Hydrogen Atom Transfer Reactivity.

Two oxoiron(IV) isomers (R 2a and R 2b) of general formula [FeIV (O)(R PyNMe3 )(CH3 CN)]2+ are obtained by reaction of their iron(II) precursor with NBu4 IO4 . The two isomers differ in the position of the oxo ligand, cis and trans to the pyridine donor. The mechanism of isomerization between R 2a and R 2b has been determined by kinetic and computational analyses uncovering an unprecedented path for interconversion of geometrical oxoiron(IV) isomers. The activity of the two oxoiron(IV) isomers in hydrogen atom transfer (HAT) reactions shows that R 2a reacts one order of magnitude faster than R 2b, which is explained by a repulsive noncovalent interaction between the ligand and the substrate in R 2b. Interestingly, the electronic properties of the R substituent in the ligand pyridine ring do not have a significant effect on reaction rates. Overall, the intrinsic structural aspects of each isomer define their relative HAT reactivity, overcoming changes in electronic properties of the ligand.

Carboxylic Acid Directed γ-Lactonization of Unactivated Primary C-H Bonds Catalyzed by Mn Complexes: Application to Stereoselective Natural Product Diversification.

Reactions that enable selective functionalization of strong aliphatic C-H bonds open new synthetic paths to rapidly increase molecular complexity and expand chemical space. Particularly valuable are reactions where site-selectivity can be directed toward a specific C-H bond by catalyst control. Herein we describe the catalytic site- and stereoselective γ-lactonization of unactivated primary C-H bonds in carboxylic acid substrates. The system relies on a chiral Mn catalyst that activates aqueous hydrogen peroxide to promote intramolecular lactonization under mild conditions, via carboxylate binding to the metal center. The system exhibits high site-selectivity and enables the oxidation of unactivated primary γ-C-H bonds even in the presence of intrinsically weaker and a priori more reactive secondary and tertiary ones at α- and ß-carbons. With substrates bearing nonequivalent γ-C-H bonds, the factors governing site-selectivity have been uncovered. Most remarkably, by manipulating the absolute chirality of the catalyst, γ-lactonization at methyl groups in gem-dimethyl structural units of rigid cyclic and bicyclic carboxylic acids can be achieved with unprecedented levels of diastereoselectivity. Such control has been successfully exploited in the late-stage lactonization of natural products such as camphoric, camphanic, ketopinic, and isoketopinic acids. DFT analysis points toward a rebound type mechanism initiated by intramolecular 1,7-HAT from a primary γ-C-H bond of the bound substrate to a highly reactive MnIV-oxyl intermediate, to deliver a carbon radical that rapidly lactonizes through carboxylate transfer. Intramolecular kinetic deuterium isotope effect and 18O labeling experiments provide strong support to this mechanistic picture.

Unveiling Halogen-Bonding Interactions between a Pyridine-Functionalized Fluoroborate Dye and Perfluorohaloarenes with Fluorescence Spectroscopy.

We have studied the halogen-bonding interactions of a pyridine-functionalized fluoroborate dye with perfluorohaloarenes (C6F6, C6F5Cl, C6F5Br, and C6F5I) in the two-component-only liquid phase using fluorescence spectroscopy. Based on the results of spectroscopic measurements and electronic-structure calculations, we have confirmed the stability only for the complex between C6F5I and the emissive dye, and it has been demonstrated that halogen-bonding interactions are accompanied by significant Stokes shifts for the ππ* band. We also provide experimental evidence that for this complex, the emission is quenched due to a simultaneous decrease of radiative and increase of nonradiative decay rate constants upon halogen-bonding interactions.

Csp2-H Amination Reactions Mediated by Metastable Pseudo-Oh Masked Aryl-CoIII-nitrene Species.

Cobalt-catalyzed C-H amination via M-nitrenoid species is spiking the interest of the research community. Understanding this process at a molecular level is a challenging task, and here we report a well-defined macrocyclic system featuring a pseudo-Oh aryl-CoIII species that reacts with aliphatic azides to effect intramolecular Csp2-N bond formation. Strikingly, a putative aryl-Co═NR nitrenoid intermediate species is formed and is rapidly trapped by a carboxylate ligand to form a carboxylate masked-nitrene, which functions as a shortcut to stabilize and guide the reaction to productive intramolecular Csp2-N bond formation. On one hand, several intermediate species featuring the Csp2-N bond formed have been isolated and structurally characterized, and the essential role of the carboxylate ligand has been proven. Complementarily, a thorough density functional theory study of the Csp2-N bond formation mechanism explains at the molecular level the key role of the carboxylate-masked nitrene species, which is essential to tame the metastability of the putative aryl-CoIII═NR nitrene species to effectively yield the Csp2-N products. The solid molecular mechanistic scheme determined for the Csp2-N bond forming reaction is fully supported by both experimental and computation complementary studies.

Evaluation of charge-transfer rates in fullerene-based donor-acceptor dyads with different density functional approximations.

The shift towards renewable energy is one of the main challenges of this generation. Dye-sensitized solar cells (DSSCs), based on donor-acceptor architectures, can help in this transition as they present excellent photovoltaic efficiencies yet cheap and simple manufacturing. For molecular heterojunction DSSCs, donor-acceptor pairs are linked in a covalent manner, which facilitates their tailoring and rational design. Nevertheless, reliable computational characterization of charge transfer rate constants (kCT) is needed to speed this development process up. In this context, the performance of time-dependent density functional theory for the calculation of kCT values in donor-acceptor fullerene-based dyads has not been benchmarked yet. Herein, we present a detailed analysis on the performance of seven well-known density functional approximations (DFAs) for this type of system, focusing on several parameters such as the reorganization energies (λ), electronic couplings (VDA), and Gibbs energies (ΔG0CT), as well as the final rate constants. The amount of exact exchange at short range (SR) and long range (LR) electron-electron distances (and the transition from the SR to LR) turned out to be key for the success of the prediction. The tuning of these parameters improves significantly the performance of current DFAs.

How Many Electrons Does a Molecular Electride Hold?

Electrides are very peculiar ionic compounds where electrons occupy the anionic positions. In a crystal lattice, these isolated electrons often form channels or surfaces, furnishing electrides with many traits with promising technological applications. Despite their huge potential, thus far, only a few stable electrides have been produced because of the intricate synthesis they entail. Due to the difficulty in assessing the presence of isolated electrons, the characterization of electrides also poses some serious challenges. In fact, their properties are expected to depend on the arrangement of these electrons in the molecule. Among the criteria that we can use to characterize electrides, the presence of a non-nuclear attractor (NNA) of the electron density is both the rarest and the most salient feature. Therefore, a correct description of the NNA is crucial to determine the properties of electrides. In this paper, we analyze the NNA and the surrounding region of nine molecular electrides to determine the number of isolated electrons held in the electride. We have seen that the correct description of a molecular electride hinges on the electronic structure method employed for the analyses. In particular, one should employ a basis set with sufficient flexibility to describe the region close to the NNA and a density functional approximation that does not suffer from large delocalization errors. Finally, we have classified these nine molecular electrides according to the most likely number of electrons that we can find in the NNA. We believe this classification highlights the strength of the electride character and will prove useful in designing new electrides.

A Unified Electro- and Photocatalytic CO2 to CO Reduction Mechanism with Aminopyridine Cobalt Complexes.

A mechanistic understanding of electro- and photocatalytic CO2 reduction is crucial to develop strategies to overcome catalytic bottlenecks. In this regard, for a new CO2-to-CO reduction cobalt aminopyridine catalyst, a detailed experimental and theoretical mechanistic study is herein presented toward the identification of bottlenecks and potential strategies to alleviate them. The combination of electrochemistry and in situ spectroelectrochemistry together with spectroscopic techniques led us to identify elusive key electrocatalytic intermediates derived from complex [LN4Co(OTf)2] (1) (LN4 = 1-[2-pyridylmethyl]-4,7-dimethyl-1,4,7-triazacyclononane) such as a highly reactive cobalt(I) (1(I)) and a cobalt(I) carbonyl (1(I)-CO) species. The combination of spectroelectrochemical studies under CO2, 13CO2, and CO with DFT disclosed that 1(I) reacts with CO2 to form the pivotal 1(I)-CO intermediate at the 1(II/I) redox potential. However, at this reduction potential, the formation of 1(I)-CO restricts the electrocatalysis due to the endergonicity of the CO release step. In agreement with the experimentally observed CO2-to-CO electrocatalysis at the CoI/0 redox potential, computational studies suggested that the electrocatalytic cycle involves striking metal carbonyls. In contrast, under photochemical conditions, the catalysis smoothly proceeds at the 1(II/I) redox potential. Under the latter conditions, it is proposed that the electron transfer to form 1(I)-CO from 1(II)-CO is under diffusion control. Then, the CO release from 1(II)-CO is kinetically favored, facilitating the catalysis. Finally, we have found that visible-light irradiation has a positive impact under electrocatalytic conditions. We envision that light irradiation can serve as an effective strategy to circumvent the CO poisoning and improve the performance of CO2 reduction molecular catalysts.

Bingel-Hirsch Addition of Diethyl Bromomalonate to Ion-Encapsulated Fullerenes M@C60 (M=Ø, Li+, Na+, K+, Mg2+, Ca2+, and Cl-).

In the last 30 years, fullerene-based materials have become popular building blocks for devices with a broad range of applications. Among fullerene derivatives, endohedral metallofullerenes (EMFs, M@Cx) have been widely studied owing to their unique properties and reactivity. For real applications, fullerenes and EMFs must be exohedrally functionalized. It has been shown that encapsulated metal cations facilitate the Diels-Alder reaction in fullerenes. Herein, the Bingel-Hirsch (BH) addition of ethyl bromomalonate over a series of ion-encapsulated M@C60 (M=Ø, Li+, Na+, K+, Mg2+, Ca2+, and Cl-; Ø@C60 stands for C60 without any endohedral metal) is quantum mechanically explored to analyze the effect of these ions on the BH addition. The results show that the incarcerated ion has a very important effect on the kinetics and thermodynamics of this reaction. Among the systems studied, K+@C60 is the one that leads to the fastest BH reaction, whereas the slowest reaction is given by Cl-@C60.

Mechanistic Insights into the ortho-Defluorination-Hydroxylation of 2-Halophenolates Promoted by a Bis(µ-oxo)dicopper(III) Complex.

C-F bonds are one of the most inert functionalities. Nevertheless, some [Cu2O2]2+ species are able to defluorinate-hydroxylate ortho-fluorophenolates in a chemoselective manner over other ortho-halophenolates. Albeit it is known that such reactivity is promoted by an electrophilic attack of a [Cu2O2]2+ core over the arene ring, the crucial details of the mechanism that explain the chemo- and regioselectivity of the reaction remain unknown, and it has not being determined either if CuII2(η2:η2-O2) or CuIII2(µ-O)2 species are responsible for the initial attack on the arene. Herein, we present a combined theoretical and experimental mechanistic study to unravel the origin of the chemoselectivity of the ortho-defluorination-hydroxylation of 2-halophenolates by the [Cu2(O)2(DBED)2]2+ complex (DBED = N,N'-di-tert-butylethylenediamine). Our results show that the equilibria between (side-on)peroxo (P) and bis(µ-oxo) (O) isomers plays a key role in the mechanism, with the latter being the reactive species. Furthermore, on the basis of quantum-mechanical calculations, we were able to rationalize the chemoselective preference of the [Cu2(O)2(DBED)2]2+ catalyst for the C-F activation over C-Cl and C-H activations.

A new tuned range-separated density functional for the accurate calculation of second hyperpolarizabilities.

The calculation of nonlinear optical properties (NLOPs) using density functional theory (DFT) remains a challenge in computational chemistry. Although the existing range-separated functionals display the best performance for the calculation of this type of properties, their errors strongly depend on the family of molecules studied. Herein, we have explored a new strategy to empirically tune the range-separated LC-BLYP method to improve the accuracy of the calculation of the second hyperpolarizabilities (γ), which are poorly described by current density functional approximations. First, we benchmarked nine of the most accurate commonly used range-separated hybrid and optimally tuned functionals (i.e. B3LYP, PBE0, BH&HLYP, M06-2X, MN15, ωB97X-D, CAM-B3LYP, LC-BLYP and OT-LC-BLYP) for the calculation of γ using as a reference the CCSD(T) values of a chemically diverse set of 60 molecules. Among these nine functionals, LC-BLYP gives the lowest average errors. We determined the value of the range-separation parameter ω required to reproduce the CCSD(T) second hyperpolarizabilities with the LC-BLYP functional (ωCC) for the set of 60 molecules. Our new tuned range-separated functional, Tα-LC-BLYP, uses a quadratic correlation between ωCC and a molecular descriptor in terms of the linear polarizability and the number of electrons in the molecule. The average error of the γ values obtained with Tα-LC-BLYP is reduced by half or more as compared with the most accurate among the nine density functional approximations benchmarked.

Partitioning of interaction-induced nonlinear optical properties of molecular complexes. II. Halogen-bonded systems.

Following our study on hydrogen-bonded (HB) complexes [Phys. Chem. Chem. Phys., 2018, 20, 19841], the physical nature of interaction-induced (non)linear optical properties of another important class of molecular complexes, namely halogen-bonded (XB) systems, was analyzed in this study. The excess electronic and nuclear relaxation (hyper)polarizabilities of nine representative XB complexes covering a wide range of halogen-bond strengths were computed. The partitioning of the excess properties into individual interaction-energy components (electrostatic, exchange, induction, dispersion) was performed by using the variational-perturbational energy decomposition scheme at the MP2/aug-cc-pVTZ level of theory and further supported by calculations with the SCS-MP2 method. In the case of the electronic interaction-induced properties, the physical composition of Δαel and Δγel was found to be very similar for the two types of bonding, despite the different nature of the binding. For Δßel, the XB complexes exhibit a more systematic interplay of interaction-energy contributions compared to the HB systems studied in the previous work. Our analysis revealed that the patterns of interaction-energy contributions to the interaction-induced nuclear-relaxation contributions to the linear polarizability and the first hyperpolarizability are very similar. For both properties the exchange repulsion term is canceled out by the electrostatic and delocalization terms. The physical composition of these contributions is analogous to those observed for the HB complexes.

Design of Iron Coordination Complexes as Highly Active Homogenous Water Oxidation Catalysts by Deuteration of Oxidation-Sensitive Sites.

The nature of the oxidizing species in water oxidation reactions with chemical oxidants catalyzed by α-[Fe(OTf)2(mcp)] (1α; mcp = N, N'-dimethyl- N, N'-bis(pyridin-2-ylmethyl)cyclohexane-1,2-diamine, OTf = trifluoromethanesulfonate anion) and ß-[Fe(OTf)2(mcp)] (1ß) has been investigated. Mössbauer spectroscopy provides definitive evidence that 1α and 1ß generate oxoiron(IV) species as the resting state. Decomposition paths of the catalysts have been investigated by identifying and quantifying ligand fragments that form upon degradation. This analysis correlates the water oxidation activity of 1α and 1ß with stability against oxidative damage of the ligand via aliphatic C-H oxidation. The site of degradation and the relative stability against oxidative degradation are shown to be dependent on the topology of the catalyst. Furthermore, the mechanisms of catalyst degradation have been rationalized by computational analyses, which also explain why the topology of the catalyst enforces different oxidation-sensitive sites. This information has served in creating catalysts where sensitive C-H bonds have been replaced by C-D bonds. The deuterated analogues D4-α-[Fe(OTf)2(mcp)] (D4-1α), D4-ß-[Fe(OTf)2(mcp)] (D4-1ß), and D6-ß-[Fe(OTf)2(mcp)] (D6-1ß) were prepared, and their catalytic activity has been studied. D4-1α proves to be an extraordinarily active and efficient catalyst (up to 91% of O2 yield); it exhibits initial reaction rates identical with those of its protio analogue, but it is substantially more robust toward oxidative degradation and yields more than 3400 TON ( n(O2)/ n(Fe)). Altogether this evidences that the water oxidation catalytic activity is performed by a well-defined coordination complex and not by iron oxides formed after oxidative degradation of the ligands.

Partition of optical properties into orbital contributions.

Nonlinear optical properties (NLOPs) play a major role in photonics, electro-optics and optoelectronics, and other fields of modern optics. The design of new NLO molecules and materials has benefited from the development of computational tools to analyze the relationship between the electronic structure of molecules and their optical response. In this paper, we present a new means to analyze the response property through the partition of NLOPs in terms of orbital contributions (PNOC). This tool can be used to obtain a real-space representation of the NLOPs, providing a powerful visualization aid to connect the magnitude of the optical property with some parts of the molecule. Unlike other methods to analyze NLOPs, the PNOC decomposes the optical property into orbitals of the unperturbed system, furnishing this method with the ability to assess the performance of single- and multi-determinant electronic structure methods. PNOC can be also used to design small basis sets for an accurate description of large systems, saving a substantial amount of computer time for the calculation of optical properties.

Metal Cluster Electrides: A New Type of Molecular Electride with Delocalised Polyattractor Character.

Electrides are ionic substances containing isolated electrons. These confined electrons are topologically characterised by a quasi-atom, that is, a non-nuclear attractor (NNA) of the electron density. The electronic structure of the octahedral 4 A1g Li6+ and 5 A1g Be6 species shows that these species have a large number of NNAs. These NNAs have highly delocalised electron densities and, as a result, the chemical bonding pattern of these systems is reminiscent of that in solid metals, in which metal cations are surrounded by a "sea" of delocalised valence electrons. We propose the term metal cluster electrides to refer to this new class of compounds. In this study, we establish a computational protocol to identify, characterize, and design metal cluster electrides and we elucidate the intricate bonding patterns of this particular type of species.