The electrophilic aromatic bromination of benzenes: mechanistic and regioselective insights from density functional theory.

The HBr-assisted electrophilic aromatic bromination of benzene, anisole and nitrobenzene was investigated using static DFT calculations in gas phase and implicit apolar (CCl4) and polar (acetonitrile) solvent models at the ωB97X-D/cc-pVTZ level of theory. The reaction profiles corresponding to either a direct substitution reaction or an addition-elimination process were constructed and insight into the preferred regioselectivity was provided using a combination of conceptual DFT reactivity indices, aromaticity indices, Wiberg bond indices and the non-covalent interaction index. Our results show that under the considered reaction conditions the bromination reaction preferentially occurs through an addition-elimination mechanism and without formation of a stable charged Wheland intermediate. The ortho/para directing effect of the electron-donating methoxy-group in anisole was ascribed to a synergy between strong electron delocalisation and attractive interactions. In contrast, the preferred meta-addition on nitrobenzene could not be traced back to any of these effects, nor to the intrinsic reactivity property of the reactant. In this case, an electrostatic clash between the ipso-carbon of the ring and the nitrogen atom resulting from the later nature of the rate-determining step, with respect to anisole, appeared to play a crucial role.

Understanding the Reactivity of Supported Late Transition Metals on a Bare Anatase (101) Surface: A Periodic Conceptual DFT Investigation.

The rapidly growing interest for new heterogeneous catalytic systems providing high atomic efficiency along with high stability and reactivity triggered an impressive progress in the field of single-atom catalysis. Nevertheless, unravelling the factors governing the interaction strength between the support and the adsorbed metal atoms remains a major challenge. Based on periodic density functional theory (DFT) calculations, this paper provides insight into the adsorption of single late transition metals on a defect-free anatase surface. The obtained adsorption energies fluctuate, with the exception of Pd, between -3.11 and -3.80 eV and are indicative of a strong interaction. Depending on the considered transition metal, we could attribute the strength of this interaction with the support to i) an electron transfer towards anatase (Ru, Rh, Ni), ii) s-d orbital hybridisation effects (Pt), or iii) a synergistic effect between both factors (Fe, Co, Os, Ir). The driving forces behind the adsorption were also found to be strongly related to Klechkowsky's rule for orbital filling. In contrast, the deviating behaviour of Pd is most likely associated with the lower dissociation enthalpy of the Pd-O bond. Additionally, the reactivity of these systems was evaluated using the Fermi weighted density of states approach. The resulting softness values can be clearly related to the electron configuration of the catalytic systems as well as with the net charge on the transition metal. Finally, these indices were used to construct a model that predicts the adsorption strength of CO on these anatase-supported d-metal atoms. The values obtained from this regression model show, within a 95 % probability interval, a correlation of 84 % with the explicitly calculated CO adsorption energies.

Synergistic Effects in the Activity of Nano-Transition-Metal Clusters Pt12 M (M=Ir, Ru or Rh) for NO Dissociation.

The dissociation of environmentally hazardous NO through dissociative adsorption on metallic clusters supported by oxides, is receiving growing attention. Building on previous research on monometallic M13 clusters [The Journal of Physical Chemistry C 2019, 123 (33), 20314-20318], this work considers bimetallic Pt12 M (M=Rh, Ru or Ir) clusters. The adsorption energy and activation energy of NO dissociation on the clusters have been calculated in vacuum using Kohn-Sham DFT, while their trends were rationalized using reactivity indices such as molecular electrostatic potential and global Fermi softness. The results show that doping of the Pt clusters lowered the adsorption energy as well as the activation energy for NO dissociation. Furthermore, reactivity indices were calculated as a first estimate of the performance of the clusters in realistic amorphous silica pores (MCM-41) through ab initio molecular dynamics simulations.

Reactivity of Single Transition Metal Atoms on a Hydroxylated Amorphous Silica Surface: A Periodic Conceptual DFT Investigation.

The drive to develop maximal atom-efficient catalysts coupled to the continuous striving for more sustainable reactions has led to an ever-increasing interest in single-atom catalysis. Based on a periodic conceptual density functional theory (cDFT) approach, fundamental insights into the reactivity and adsorption of single late transition metal atoms supported on a fully hydroxylated amorphous silica surface have been acquired. In particular, this investigation revealed that the influence of van der Waals dispersion forces is especially significant for a silver (98 %) or gold (78 %) atom, whereas the oxophilicity of the Group 8-10 transition metals plays a major role in the interaction strength of these atoms on the irreducible SiO2 support. The adsorption energies for the less-electronegative row 4 elements (Fe, Co, Ni) ranged from -1.40 to -1.92 eV, whereas for the heavier row 5 and 6 metals, with the exception of Pd, these values are between -2.20 and -2.92 eV. The deviating behavior of Pd can be attributed to a fully filled d-shell and, hence, the absence of the hybridization effects. Through a systematic analysis of cDFT descriptors determined by using three different theoretical schemes, the Fermi weighted density of states approach was identified as the most suitable for describing the reactivity of the studied systems. The main advantage of this scheme is the fact that it is not influenced by fictitious Coulomb interactions between successive, charged reciprocal cells. Moreover, the contribution of the energy levels to the reactivity is simultaneously scaled based on their position relative to the Fermi level. Finally, the obtained Fermi weighted density of states reactivity trends show a good agreement with the chemical characteristics of the investigated metal atoms as well as the experimental data.

Synthesis and Reactivity of Spirocarbocycles as Scaffolds for Nucleoside Analogues.

A novel class of substituted spiro[3.4]octanes can be accessed via a [2 + 2]-cycloaddition of dichloroketene on a readily prepared exo-methylene cyclopentane building block. This reaction sequence was found to be robust on a multigram scale and afforded a central spirocyclobutanone scaffold for carbocyclic nucleosides. The reactivity of this constrained building block was evaluated and compared to the corresponding 4'-spirocyclic furanose analogues. Density functional theory calculations were performed to support the observed selectivity in the carbonyl reduction of spirocyclobutanone building blocks. Starting from novel spirocyclic intermediates, we exemplified the preparation of an undescribed class of carbocyclic nucleoside analogues and provided a proof of concept for application as inhibitors for the protein methyltransferase target PRMT5.

Stereoselective Reductions of 3-Substituted Cyclobutanones: A Comparison between Experiment and Theory.

The stereoselective reduction of carbonyls is of key importance in the total synthesis of natural products and in medicinal chemistry. Nevertheless, models for rationalizing the stereoselectivity of the hydride reductions of cyclobutanones toward cyclobutanols are largely lacking, unlike cyclohexanone reductions. In order to elucidate the factors that control the stereoselectivity of these reductions, we have investigated the effect of the reaction temperature, solvent, substituent, and type of reducing agent using a synergistic experimental-computational approach. On the experimental side, the hydride reduction of 3-substituted cyclobutanones was proven to be highly selective for the formation of cis alcohol (>90%), irrespective of the size of the hydride reagent. The pronounced selectivity can be further enhanced by lowering the reaction temperature or decreasing the solvent polarity. On the computational side, density functional theory and noncovalent interaction analysis reveal that torsional strain plays a major role in the preference for the antifacial hydride approach, consistent with the Felkin-Anh model. In the presence of the benzyloxy substituent, the high selectivity for the cis isomer is also driven by repulsive electrostatic interactions in the case of a syn-facial hydride attack. The computed cis/trans ratios are in good agreement with the experimental ones and thus show the potential of computational chemistry for predicting and rationalizing the stereoselectivity of hydride reductions of cyclobutanones.

Organogermanium(II) Hydrides as a Source of Highly Soluble LiH.

The reactions of monomeric C,N-chelated organogermanium(II) hydride L(H)Ge⋅BH3 with organolithium salts RLi yielded lithium hydrogermanatoborates (Li(THF)2 {BH3 [L(H)GeR]})2 . Compound (Li(THF)2 {BH3 [L(H)GePh]})2 was used as a source of LiH for the reduction of organic C=O or C=N bonds in nonpolar solvents accompanied by the elimination of a neutral complex L(Ph)Ge⋅BH3 . The interaction of (Li(THF)2 {BH3 [L(H)GePh]})2 with the polar C=O bond was further investigated by computational studies revealing a plausible geometry of a pre-reactive intermediate. The experimental and theoretical studies suggest that, although the Li atom of (Li(THF)2 {BH3 [L(H)GePh]})2 coordinates the C=O bond, the GeH fragment is the active species in the reduction reaction. Finally, benzaldehyde was reduced by a mixture of L(H)Ge⋅BH3 with PhLi in nonpolar solvents.

Performance of Electronic Structure Methods for the Description of Hückel-Möbius Interconversions in Extended π-Systems.

Expanded porphyrins provide a versatile route to molecular switching devices due to their ability to shift between several π-conjugation topologies encoding distinct properties. DFT remains the workhorse for modeling such extended macrocycles, when taking into account their size and huge conformational flexibility. Nevertheless, the stability of Hückel and Möbius conformers depends on a complex interplay of different factors, such as hydrogen bonding, π···π stacking, steric effects, ring strain, and electron delocalization. As a consequence, the selection of an exchange-correlation functional for describing the energy profile of topological switches is very difficult. For these reasons, we have examined the performance of a variety of wave function methods and density functionals for describing the thermochemistry and kinetics of topology interconversions across a wide range of macrocycles. Especially for hexa- and heptaphyrins, the Möbius structures have a stronger degree of static correlation than the Hückel and twisted-Hückel structures, and as a result the relative energies of singly twisted structures are a challenging test for electronic structure methods. Comparison of limited orbital space full CI calculations with CCSD(T) calculations within the same active spaces shows that post-CCSD(T) correlation contributions to relative energies are very minor. At the same time, relative energies are weakly sensitive to further basis set expansion, as proven by the minor energy differences between the extrapolated MP2/CBS energies estimated from cc-pV{T,Q}Z, diffuse-augmented heavy-aug-cc-pV{T,Q}Z and explicitly correlated MP2-F12/cc-pVDZ-F12 calculations. Hence, our CCSD(T) reference values are reasonably well-converged in both 1-particle and n-particle spaces. While conventional MP2 and MP3 yield very poor results, SCS-MP2 and particularly SOS-MP2 and SCS-MP3 agree to better than 1 kcal mol-1 with the CCSD(T) relative energies. Regarding DFT methods, the range-separated double hybrids, such as ωB97M(2) and B2GP-PLYP, outperform other functionals with RMSDs of 0.6 and 0.8 kcal mol-1, respectively. While the original DSD-PBEP86 double hybrid performs fairly poorly for these extended π-systems, the errors drop down to 1.9 kcal mol-1 for the revised revDOD-PBEP86-NL, which eliminates the same-spin correlation energy. Minnesota meta-GGA functionals with high fractions of exact exchange (M06-2X and M08-HX) also perform reasonably well, outperforming more robust and significantly less empirically parametrized functionals like SCAN0-D3.

A Benchmark of Density Functional Approximations For Thermochemistry and Kinetics of Hydride Reductions of Cyclohexanones.

The performance of density functionals and wavefunction methods for describing the thermodynamics and kinetics of hydride reductions of 2-substituted cyclohexanones has been evaluated for the first time. A variety of exchange correlation functionals ranging from generalized gradient approximations to double hybrids have been tested and their performance to describe the facial selectivity of hydride reductions of cyclohexanones has been carefully assessed relative to the CCSD(T) method. Among the tested methods, an approach in which single-point energy calculations using the double hybrid B2PLYP-D3 functional on ωB97X-D optimized geometries provides the most accurate transition state energies for these kinetically-controlled reactions. Moreover, the role of torsional strain, temperature, solvation, noncovalent interactions on the stereoselectivity of these reductions was elucidated. Our results indicate a prominent role of the substituent on the cis/trans ratios driven by the delicate interplay between torsional strain and dispersion interactions.

Aromatic sulfonation with sulfur trioxide: mechanism and kinetic model.

Electrophilic aromatic sulfonation of benzene with sulfur trioxide is studied with ab initio molecular dynamics simulations in gas phase, and in explicit noncomplexing (CCl3F) and complexing (CH3NO2) solvent models. We investigate different possible reaction pathways, the number of SO3 molecules participating in the reaction, and the influence of the solvent. Our simulations confirm the existence of a low-energy concerted pathway with formation of a cyclic transition state with two SO3 molecules. Based on the simulation results, we propose a sequence of elementary reaction steps and a kinetic model compatible with experimental data. Furthermore, a new alternative reaction pathway is proposed in complexing solvent, involving two SO3 and one CH3NO2.