Double Hybrid Functionals and the Π-System Bond Length Alternation Challenge: Rivaling Accuracy of Post-HF Methods

Predicting accurate bond length alternations (BLAs) in long conjugated oligomers has been a significant challenge for electronic-structure methods for many decades, made particularly important by the close relationships between BLA and the rich optoelectronic properties of π-delocalized systems. Here, we test the accuracy of recently developed, and increasingly popular, double hybrid (DH) functionals, positioned at the top of Jacobs Ladder of DFT methods of increasing sophistication, computational cost, and accuracy, due to incorporation of MP2 correlation energy. Our test systems comprise oligomeric series of polyacetylene, polymethineimine, and polysilaacetylene up to six units long. MP2 calculations reveal a pronounced shift in BLAs between the 6-31G(d) basis set used in many studies of BLA to date and the larger cc-pVTZ basis set, but only modest shifts between cc-pVTZ and aug-cc-pVQZ results. We hence perform new reference CCSD(T)/cc-pVTZ calculations for all three series of oligomers against which we assess the performance of several families of DH functionals based on BLYP, PBE, and TPSS, along with lower-rung relatives including global- and range-separated hybrids. Our results show that DH functionals systematically improve the accuracy of BLAs relative to single hybrid functionals. xDH-PBE0 (N4 scaling using SOS-MP2) emerges as a DH functional rivaling the BLA accuracy of SCS-MP2 (N5 scaling), which was found to offer the best compromise between computational cost and accuracy the last time the BLA accuracy of DFT- and wave function-based methods was systematically investigated. Interestingly, xDH-PBE0 (XYG3), which differs to other DHs in that its MP2 term uses PBE0 (B3LYP) orbitals that are not self-consistent with the DH functional, is an outlier of trends of decreasing average BLA errors with increasing fractions of MP2 correlation and HF exchange.

Enantioselective Synthesis of Succinimides by Michael Addition of Aldehydes to Maleimides Organocatalyzed by Chiral Primary Amine-Guanidines

The monoguanylation of (1S,2S)- and (1R,2R)-cyclohexane-1,2-diamine affords chiral primary amine-guanidines that are used as chiral organocatalysts in the enantioselective Michael addition of aldehydes, particularly α,α-disubstituted aldehydes, to maleimides. The reaction is carried out in the presence of imidazole, as an additive, in aqueous N,N-dimethylformamide, as the solvent, and affords the corresponding enantioenriched succinimides in high or quantitative yields with enantioselectivities up to 96 % ee. Theoretical calculations (DFT and M06–2X) suggest a different hydrogen-bonding coordination pattern between the maleimide (C=O) and the catalyst (NH groups) is responsible for the enantioinduction switch that is observed when the reaction is carried out using primary amine-guanidines versus primary amine-thioureas as the organocatalysts.

Polarized–unpolarized ground state of small polycyclic aromatic hydrocarbons

Do polyacenes, circumacenes, periacenes, nanographenes, and graphene nanoribbons show a spin polarized ground state? In this work, we present monodeterminantal (Hartree–Fock (HF) and density functional theory (DFT) types), and multideterminantal calculations (Møller–Plesset and Coupled Cluster), for several families of unsaturated organic molecules (n-Acenes, n-Periacenes and n-Circumacenes). All HF calculations and many DFT show a spin-polarized (antiferromagnetic) ground state, in agreement with previous calculations. Nevertheless, the multideterminantal calculations, carried out with perturbative and variational wavefunctions, show that the more stable state is obtained starting from the unpolarized HF wavefunction. The trend of the stabilization of wavefunctions (polarized or unpolarized) with respect to exchange and correlation potentials, and to the number of benzene rings, has been analyzed. A study of the spin (〈Ŝ2〉) and the spin density on the carbon atoms has also been carried out.

Theoretical Rationalization of the Singlet–Triplet Gap in OLEDs Materials: Impact of Charge-Transfer Character

New materials for OLED applications with low singlet–triplet energy splitting have been recently synthesized in order to allow for the conversion of triplet into singlet excitons (emitting light) via a Thermally Activated Delayed Fluorescence (TADF) process, which involves excited-states with a non-negligible amount of Charge-Transfer (CT). The accurate modeling of these states with Time-Dependent Density Functional Theory (TD-DFT), the most used method so far because of the favorable trade-off between accuracy and computational cost, is however particularly challenging. We carefully address this issue here by considering materials with small (high) singlet–triplet gap acting as emitter (host) in OLEDs and by comparing the accuracy of TD-DFT and the corresponding Tamm-Dancoff Approximation (TDA), which is found to greatly reduce error bars with respect to experiments thanks to better estimates for the lowest singlet–triplet transition. Finally, we quantitatively correlate the singlet–triplet splitting values with the extent of CT, using for it a simple metric extracted from calculations with double-hybrid functionals, that might be applied in further molecular engineering studies.

Modelling liquid–liquid equilibria for island type ternary systems

Closed miscibility gaps in ternary liquid mixtures, at constant temperature and pressure, are obtained if phase separations occur only in the ternary region, whilst all binary mixtures involved in the system are completely miscible. This type of behaviour, although not very frequent, has been observed for a certain number of systems. Nevertheless, we have found no information about the applicability of the common activity coefficient models, as NRTL and UNIQUAC, for these types of ternary systems. Moreover, any of the island type systems published in the most common liquid–liquid equilibrium data collections, are correlated with any model. In this paper, the applicability of the NRTL equation to model the LLE of island type systems is assessed using topological concepts related to the Gibbs stability test. A first attempt to correlate experimental LLE data for two island type ternary systems is also presented.

Phase equilibria of the water + 1-butanol + toluene ternary system at 101.3 kPa

Isobaric vapour–liquid and vapour–liquid–liquid equilibrium data for the water + 1-butanol + toluene ternary system were measured at 101.3 kPa with a modified VLE 602 Fischer apparatus. In addition, the liquid–liquid equilibrium data at 313.15 K were measured and compared with data from other authors at different temperatures. The system exhibits a ternary heterogeneous azeotrope whose temperature and composition have been determined by interpolation. The thermodynamic consistency of the experimental vapour–liquid and vapour–liquid–liquid data was checked by means of the Wisniak’s Li/Wi consistency test. Moreover, the vapour–liquid and the liquid–liquid equilibrium correlation for the ternary system with NRTL and UNIQUAC models, together with the prediction made with the UNIFAC model, were studied and discussed.

Electronic properties of transition metal atoms on Cu2N/Cu(100)

We study the nature of spin excitations of individual transition metal atoms (Ti, V, Cr, Mn, Fe, Co, and Ni) deposited on a Cu2N/Cu(100) surface using both spin-polarized density functional theory (DFT) and exact diagonalization of an Anderson model derived from DFT. We use DFT to compare the structural, electronic, and magnetic properties of different transition metal adatoms on the surface. We find that the average occupation of the transition metal d shell, main contributor to the magnetic moment, is not quantized, in contrast with the quantized spin in the model Hamiltonians that successfully describe spin excitations in this system. In order to reconcile these two pictures, we build a zero bandwidth multi-orbital Anderson Hamiltonian for the d shell of the transition metal hybridized with the p orbitals of the adjacent nitrogen atoms, by means of maximally localized Wannier function representation of the DFT Hamiltonian. The exact solutions of this model have quantized total spin, without quantized charge at the d shell. We propose that the quantized spin of the models actually belongs to many-body states with two different charge configurations in the d shell, hybridized with the p orbital of the adjacent nitrogen atoms. This scenario implies that the measured spin excitations are not fully localized at the transition metal.

Synthetic scope and DFT analysis of the chiral binap-gold (I) complex-catalyzed 1,3-dipolar cycloaddition of azlactones with alkenes

The 1,3-dipolar cycloaddition between glycine-derived azlactones with maleimides is efficiently catalyzed by the dimeric chiral complex [(Sa)-Binap·AuTFA]2. The alanine-derived oxazolone only reacts with tert-butyl acrylate giving anomalous regiochemistry, which is explained and supported by Natural Resonance Theory and Nucleus Independent Chemical Shifts calculations. The origin of the high enantiodiscrimination observed with maleimides and tert-butyl acrylate is analyzed using DFT computed at M06/Lanl2dz//ONIOM(b3lyp/Lanl2dz:UFF) level. Several applications of these cycloadducts in the synthesis of new proline derivatives with a 2,5-trans-arrangement and in the preparation of complex fused polycyclic molecules are described.

Improvements in DFT Calculations of Spin–Spin Coupling Constants

Different types of spin–spin coupling constants (SSCCs) for several representative small molecules are evaluated and analyzed using a combination of 10 exchange functionals with 12 correlation functionals. For comparison, calculations performed using MCSCF, SOPPA, other common DFT methods, and also experimental data are considered. A detailed study of the percentage of Hartree–Fock exchange energy in SSCCs and in its four contributions is carried out. From the above analysis, a combined functional formed with local Slater (34%), Hartree–Fock exchange (66%), and P86 correlation functional (S66P86) is proposed in this paper. The accuracy of the values obtained with this hybrid functional (mean absolute deviation of 4.5 Hz) is similar to that of the SOPPA method (mean absolute deviation of 4.6 Hz).

Cleavage of the C–C Bond in the Ethanol Oxidation Reaction on Platinum. Insight from Experiments and Calculations

Using a combination of experimental and computational methods, mainly FTIR and DFT calculations, new insights are provided here in order to better understand the cleavage of the C–C bond taking place during the complete oxidation of ethanol on platinum stepped surfaces. First, new experimental results pointing out that platinum stepped surfaces having (111) terraces promote the C–C bond breaking are presented. Second, it is computationally shown that the special adsorption properties of the atoms in the step are able to promote the C–C scission, provided that no other adsorbed species are present on the step, which is in agreement with the experimental results. In comparison with the (111) terrace, the cleavage of the C–C bond on the step has a significantly lower activation energy, which would provide an explanation for the observed experimental results. Finally, reactivity differences under acidic and alkaline conditions are discussed using the new experimental and theoretical evidence.