RPA, an Accurate and Fast Method for the Computation of Static Nonlinear Optical Properties.

The accurate computation of static nonlinear optical properties (SNLOPs) in large polymers requires accounting for electronic correlation effects with a reasonable computational cost. The Random Phase Approximation (RPA) used in the adiabatic connection fluctuation theorem is known to be a reliable and cost-effective method to render electronic correlation effects when combined with density-fitting techniques and integration over imaginary frequencies. We explore the ability of the RPA energy expression to predict SNLOPs by evaluating RPA electronic energies in the presence of finite electric fields to obtain (using the finite difference method) static polarizabilities and hyperpolarizabilities. We show that the RPA based on hybrid functional self-consistent field calculations yields accurate SNLOPs as the best-tuned double-hybrid functionals developed today, with the additional advantage that the RPA avoids any system-specific adjustment.

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.

Are Accelerated and Enhanced Wave Function Methods Accurate to Compute Static Linear and Nonlinear Optical Properties?

Key components of organic-based electro-optic devices are challenging to design or optimize because they exhibit nonlinear optical responses, which are difficult to model or rationalize. Computational chemistry furnishes the tools to investigate extensive collections of molecules in the quest for target compounds. Among the electronic structure methods that provide static nonlinear optical properties (SNLOPs), density functional approximations (DFAs) are often preferred because of their low cost/accuracy ratio. However, the accuracy of the SNLOPs critically depends on the amount of exact exchange and electron correlation included in the DFA, precluding the reliable calculation of many molecular systems. In this scenario, wave function methods such as MP2, CCSD, and CCSD(T) constitute a reliable alternative to compute SNLOPs. Unfortunately, the computational cost of these methods significantly restricts the size of molecules to study, a limitation that hampers the identification of molecules with significant nonlinear optical responses. This paper analyzes various flavors and alternatives to MP2, CCSD, and CCSD(T) methods that either drastically reduce the computational cost or improve their performance but were scarcely and unsystematically employed to compute SNLOPs. In particular, we have tested RI-MP2, RIJK-MP2, RIJCOSX-MP2 (with GridX2 and GridX4 setups), LMP2, SCS-MP2, SOS-MP2, DLPNO-MP2, LNO-CCSD, LNO-CCSD(T), DLPNO-CCSD, DLPNO-CCSD(T0), and DLPNO-CCSD(T1). Our results indicate that all these methods can be safely employed to calculate the dipole moment and the polarizability with average relative errors below 5% with respect to CCSD(T). On the other hand, the calculation of higher-order properties represents a challenge for LNO and DLPNO methods, which present severe numerical instabilities in computing the single-point field-dependent energies. RI-MP2, RIJK-MP2, or RIJCOSX-MP2 are cost-effective methods to compute first and second hyperpolarizabilities with a marginal average error with respect to canonical MP2 (up to 5% for ß and up to 11% for γ). More accurate hyperpolarizabilities can be obtained with DLPNO-CCSD(T1); however, this method cannot be employed to obtain reliable second hyperpolarizabilities. These results open the way to obtain accurate nonlinear optical properties at a computational cost that can compete with current DFAs.

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.

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.

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.

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.

Theoretical characterization of sulfur-to-selenium substitution in an emissive RNA alphabet: impact on H-bonding potential and photophysical properties.

We employ density functional theory (DFT) and time-dependent DFT (TDDFT) calculations to investigate the structural, energetic and optical properties of a new computationally designed RNA alphabet, where the nucleobases, tsA, tsG, tsC, and tsU (ts-bases), have been derived by replacing sulfur with selenium in the previously reported tz-bases, based on the isothiazolo[4,3-d]pyrimidine heterocycle core. We find out that the modeled non-natural bases have minimal impact on the geometry and energetics of the classical Watson-Crick base pairs, thus potentially mimicking the natural bases in a RNA duplex in terms of H-bonding. In contrast, our calculations indicate that H-bonded base pairs involving the Hoogsteen edge of purines are destabilized as compared to their natural counterparts. We also focus on the photophysical properties of the non-natural bases and correlate their absorption/emission peaks to the strong impact of the modification on the energy of the lowest unoccupied molecular orbital. It is indeed stabilized by roughly 1.1-1.6 eV as compared to the natural analogues, resulting in a reduction of the gap between the highest occupied and the lowest unoccupied molecular orbital from 5.3-5.5 eV in the natural bases to 3.9-4.2 eV in the modified ones, with a consequent bathochromic shift in the absorption and emission spectra. Overall, our analysis clearly indicates that the newly modelled ts-bases are expected to exhibit better fluorescent properties as compared to the previously reported tz-bases, while retaining similar H-bonding properties. In addition, we show that a new RNA alphabet based on size-extended benzo-homologated ts-bases can also form stable Watson-Crick base pairs with the natural complementary nucleobases.

The preference for dual-gold(i) catalysis in the hydro(alkoxylation vs. phenoxylation) of alkynes.

Dinuclear gold complexes and their use in catalysis have received significant recent attention, but there are few critical comparisons of mono- versus dual gold-catalysed pathways. Herein we study the hydroalkoxylation and hydrophenoxylation of alkynes using density functional theory calculations, and compare two possible mechanisms that have been proposed previously on the basis of theoretical and experimental studies, which unravel different preferences because of both the nature of the alkyne and alcohol and the non-innocent role of the counter-anion of a dual gold based catalyst. Entropy is found to have a significant effect, rendering the nucleophilic attack of the monoaurated intermediate [Au(L)(η2-alkyne)]+ difficult both kinetically and thermodynamically; this mechanism cannot easily form only the trans-alkene product that is observed experimentally. Instead, a reaction via a dual gold catalysed mechanism presents much lower barriers. In addition, for the sake of direct comparison with recent results by Belanzoni and Zuccaccia, oversimplification of the N-heterocyclic carbene (NHC) ligand in the calculations might decrease the enthalpy barrier and lead to results that are not directly applicable to experiments. Moreover, the alkylic or arylic nature of the alkyne and/or alcohol is also tested.