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.
Speeding addition to fullerenes: Similar to the Diels–Alder (DA) cycloaddition, the Bingel–Hirsch (BH) addition of bromomalonate to C60 is accelerated by the presence of monocations (Li+, Na+, K+) encapsulated in the fullerene. However, at variance with the DA cycloaddition, the performance of the BH deteriorates when changing the monocations to dications (Mg2+, Ca2+) owing to the increase of the energy barriers in the last step of the BH addition.
In the oriented external electric-field-driven catalysis, the reaction rates and the selectivity of chemical reactions can be tuned at will. The activation barriers of chemical reactions within ...external electric fields of several strengths and directions can be computationally modeled. However, the calculation of all of the required field-dependent transition states and reactants is computationally demanding, especially for large systems. Herein, we present a method based on the Taylor expansion of the field-dependent energy of the reactants and transition states in terms of their field-free dipole moments and electrical (hyper)polarizabilities. This approach, called field-dependent energy barrier (FDBβ), allows systematic one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) representations of the activation energy barriers for any strength and direction of the external electric field. The calculation of the field-dependent FDBβ energy barriers has a computational cost several orders of magnitude lower than the explicit electric field optimizations, and the errors of the FDBβ barriers are within the range of only 1–2 kcal·mol–1. The achieved accuracy is sufficient for a fast-screening tool to study and predict potential electric-field-induced catalysis, regioselectivity, and stereoselectivity. As illustrative examples, four cycloadditions (1,3-dipolar and Diels–Alder) are studied.
C–F bonds are one of the most inert functionalities. Nevertheless, some Cu2O22+ 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 Cu2O22+ 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 CuII 2(η2:η2-O2) or CuIII 2(μ-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)22+ 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)22+ catalyst for the C–F activation over C–Cl and C–H activations.
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.
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.
Among the nine functionals benchmarked, the most accurate
γ
are obtained by Tα-LC-BLYP, reducing about half the errors of LC-BLYP.
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,
A,
G,
C, and
U (ts-bases), have been derived by replacing sulfur with selenium in the previously reported tz-bases, based on the isothiazolo4,3-dpyrimidine 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 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.
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 (
k
CT
) is needed to speed this development process up. In this context, the performance of time-dependent density functional theory for the calculation of
k
CT
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 (
V
DA
), and Gibbs energies (Δ
G
0
CT
), 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.
Charge-transfer rate constants (
k
CT
) are required for the design of fullerene-based solar cells. We have assessed the performance of six long-range corrected hybrid functionals and a meta-hybrid-GGA functional for the calculation of
k
CT
with the Marcus equation.
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@C
) 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@C
(M=Ø, Li
, Na
, K
, Mg
, Ca
, and Cl
; Ø@C
stands for C
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
@C
is the one that leads to the fastest BH reaction, whereas the slowest reaction is given by Cl
@C
.