Theoretical Study of the Regioselectivity of [2 + 2] Photocycloaddition Reactions of Acrolein with Olefins

Minimal Active Space for Diradicals Using Multistate Density Functional Theory

This work explores the electronic structure as well as the reactivity of singlet diradicals, making use of multistate density functional theory (MSDFT). In particular, we show that a minimal active space of two electrons in two orbitals is adequate to treat the relative energies of the singlet and triplet adiabatic ground state as well as the first singlet excited state in many cases. This is plausible because dynamic correlation is included in the first place in the optimization of orbitals in each determinant state via block-localized Kohn-Sham density functional theory. In addition, molecular fragment, i.e., block-localized Kohn-Sham orbitals, are optimized separately for each determinant, providing a variational diabatic representation of valence bond-like states, which are subsequently used in nonorthogonal state interactions (NOSIs). The computational procedure and its performance are illustrated on some prototypical diradical species. It is shown that NOSI calculations in MSDFT can be used to model bond dissociation and hydrogen-atom transfer reactions, employing a minimal number of configuration state functions as the basis states. For p- and s-types of diradicals, the closed-shell diradicals are found to be more reactive than the open-shell ones due to a larger diabatic coupling with the final product state. Such a diabatic representation may be useful to define reaction coordinates for electron transfer, proton transfer and coupled electron and proton transfer reactions in condensed-phase simulations.

Unifying Conceptual Density Functional and Valence Bond Theory: The Hardness-Softness Conundrum Associated with Protonation Reactions and Uncovering Complementary Reactivity Modes

In this study, we address the long-standing issue-arising prominently from conceptual density functional theory (CDFT)-of the relative importance of electrostatic, i.e., "hard-hard", versus spin-pairing, i.e., "soft-soft", interactions in determining regiochemical preferences. We do so from a valence bond (VB) perspective and demonstrate that VB theory readily enables a clear-cut resolution of both of these contributions to the bond formation/breaking process. Our calculations indicate that appropriate local reactivity descriptors can be used to gauge the magnitude of both interactions individually, e.g., Fukui functions or HOMO/LUMO orbitals for the spin-pairing/(frontier) orbital interactions and molecular electrostatic potentials (and/or partial charges) for the electrostatic interactions. In contrast to previous reports, we find that protonation reactions cannot generally be classified as either charge- or frontier orbital-controlled; instead, our results indicate that these two bonding contributions generally interplay in more subtle patterns, only giving the impression of a clear-cut dichotomy. Finally, we demonstrate that important covalent, i.e., spin pairing, reactivity modes can be missed when only a single spin-pairing/orbital interaction descriptor is considered. This study constitutes an important step in the unification of CDFT and VB theory.

How Do Local Reactivity Descriptors Shape the Potential Energy Surface Associated with Chemical Reactions? The Valence Bond Delocalization Perspective

How do local reactivity descriptors, such as the Fukui function and the local spin density distribution, shape the potential energy surface (PES) associated with chemical reactions and thus govern reactivity trends and regioselective preferences? This is the question that is addressed here through a qualitative valence bond (VB) analysis. We demonstrate that common density functional theory (DFT)-based local reactivity descriptors can essentially be regarded-in one way or another-as indirect measures of delocalization, i.e., resonance stabilization, of the reactants within VB theory. The inherent connection between (spatial) delocalization and (energetic) resonance stabilization embedded in VB theory provides a natural and elegant framework for analyzing and comprehending the impact of individual local reactivity descriptors on the global PES. Our analysis provides new insights into the role played by local reactivity descriptors and illustrates under which conditions they can sometimes fail to predict reactivity trends and regioselective preferences, e.g., in the case of ambident reactivity. This treatment constitutes a first step toward a unification of VB theory and conceptual DFT.

Real-Space Approach to the Reaction Force: Understanding the Origin of Synchronicity/Nonsynchronicity in Multibond Chemical Reactions

In this article, we present a complementary analysis based on the reaction force F(ξ)/reaction force constant κ(ξ) and noncovalent interactions (NCI) index to characterize the energetics (kinetic and thermodynamics) and mechanistic pathways of two sets of multibond chemical reactions, namely, two double-proton transfer and two Diels-Alder cycloaddition reactions. This approach offers a very straightforward and useful way to delve into a deeper understanding of this type of process. While F(ξ) allows the partition of the whole pathway into three regions or phases, κ(ξ) describes how orchestrated are the bond-breaking and bond-formation events. In turn, NCI indicates how the inter- and intramolecular bonds evolve. The most innovative aspect is the inclusion of the formation of the reactant complex along the pathway, which, by means of NCI, unveils the early molecular recognition and the comprehension of its role in determining the degree of the synchronicity/nonsynchronicity of one-step processes. This approach should be a useful and alternative tool to characterize the energetics and the mechanism of general chemical reactions.

Further understanding of the Ru-centered [2+2] cycloreversion/cycloaddition involved into the interconversion of ruthenacyclobutane using the Grubbs catalysts from a reaction force analysis

The chemical reactivity of the first- and second-generation Grubbs catalysts has always been a significant issue in olefin metathesis. In the present work, we study the [2+2] cycloreversion/cycloaddition and the alkylidene rotation involved into the interconversion of the ruthenacyclobutane intermediate, through the reaction force and reaction force constant analysis. It has been found that the structural contribution controls the barrier energy in the interconversion of ruthenacyclobutane via [2+2] cycloreversion/cycloaddition, which is slightly lower in the second generation of Grubbs catalysts while its electronic contribution is slightly higher, which unveils a major rigidity and donor/acceptor properties of the NHC. This finding explains a greater structural contribution in the rate constant. Moreover, on the basis of the reaction force constant, the process can be classified as "two-stage"-concerted reactions, noting a more asynchronous process when the first generation is used as a catalyst.Finally, a similar analysis into the alkylidene rotation was performed. It was determined that [2+2] cycloreversion and alkylidene rotations take place in a sequential manner, the energy barrier is again controlled by structural reorganization, and the pathway is less asynchronous.

Effect of the exchange–correlation functional on the synchronicity/nonsynchronicity in bond formation in Diels–Alder reactions: a reaction force constant analysis

The performance of 24 KS-DFT-based methods (GGA, MGGA, HGGA, HMGGA, and DHGGA) was assessed, finding that M11 and M06-2X (HMGGA) predicting reliable TS geometries, energetics, and (a)synchronicities in Diels–Alder reactions.

Scrutinizing the substituent effect on Mo-based electrocatalysts for molecular hydrogen release through axial–equatorial decomposition: a DFT study

The effect of electron-donating and -withdrawing groups is studied in terms of axial and equatorial locations at the p-pyridine rings in the [(PY₅Me₂)MoO]²⁺ for the H₂ release in aqueous media. It opens a way to modify its kinetics and thermodynamics.

Hetero-intermolecular [2+2] photocycloaddition of 1,4-dihydropyridines: a combined experimental and DFT study

A detailed DFT mechanistic study of hetero-intermolecular [2+2] photocycloaddition of 1,4-DHPs was carried out for investigating the regioselectivity of such transformation.

Solvent effect on the degree of (a)synchronicity in polar Diels-Alder reactions from the perspective of the reaction force constant analysis

In this work, we computationally evaluated the influence of six different molecular solvents, described as a polarizable continuum model at the M06-2X/6-31+G(d,p) level, on the activation barrier/reaction rate, overall energy change, TS geometry, and degree of (a)synchronicity of two concerted Diels-Alder cycloadditions of acrolein (R1) and its complex with Lewis acid acrolein···BH₃ (R2) to cyclopentadiene. In gas-phase, we found that both exothermicity and activation barrier are only reduced by about 2.0 kcal mol^-1, and the asynchronicity character of the mechanism is accentuated when BH₃ is included. An increment in the solvent's polarity lowers the activation energy of R1 by 1.3 kcal mol^-1, while for R2 the reaction rate is enhanced by more than 2000 times at room temperature (i.e., the activation energy decreases by 4.5 kcal mol^-1) if the highest polar media is employed. Therefore, a synergistic effect is achieved when both external agents, i.e., Lewis acid catalyst and polar solvent, are included together. This effect was ascribed to the ability of the solvent to favor the encounter between cyclopentadiene and acrolein···BH₃. This was validated by the asymmetry of the TS which becomes highly pronounced when either both or just BH₃ is considered or the solvent's polarity is increased. Finally, the reaction force constant κ(ξ) reveals that an increment in the solvent's polarity is able to turn a moderate asynchronous mechanism of the formation of the new C-C σ-bonds into a highly asynchronous one. Graphical abstract A synergistic effect is achieved when both external agents, i.e., Lewis acid catalyst and polar solvent, are included together: lowered energy barriers and increased asynchronicities.

Mechanistic Investigation of Visible-Light-Induced Intermolecular [2 + 2] Photocycloaddition Catalyzed with Chiral Thioxanthone

The recent thioxanthone-sensitizer-catalyzed intermolecular [2 + 2] cycloaddition induced by visible-light irradiation set the stage for the future development of feasible photocycloadditions. Nonetheless, the mechanism of this reaction still remains under debate, especially on the activation mode of the thioxanthone photosensitizer (energy transfer, bielectron exchange, and hydrogen transfer are all possible mechanisms). To settle this issue, systematic density functional theory calculations have been carried out. The results indicate that the energy-transfer pathway is more favorable than the bielectron-exchange and the hydrogen-transfer pathways. Meanwhile, the overall transformations involve the complexation and excitation of photosensitizer, the first C-C bond formation, the dissociation of the sensitizer, the triplet-to-singlet electronic state crossing, and the second C-C bond formation. The first C-C bond formation is the rate- and selectivity-determining step, and synergistic energy and electron transfer from photosensitizer to substrate moieties takes place along this process. On this basis, the effect of olefin substrates (ethyl vinyl ketone vs vinyl acetate) on the stereoselectivity was finally analyzed.

Cavity Ringdown Spectrum of 2-Cyclohexen-1-one in the CO/Alkenyl CC Stretch Region of the S1(n, π*)-S0 Vibronic Band System

The 2-cyclohexen-1-one (2CHO) molecule serves as a prototype for understanding the photochemical properties of conjugated enones. We have recorded the cavity ringdown (CRD) absorption spectrum of 2CHO vapor at room temperature over the 360-380 nm range. This portion of the spectrum encompasses the S₁(n,π*) ← S₀ vibronic band system in the region of the C═C and C═O stretch fundamentals. We have assigned about 40 vibronically resolved features in the spectrum, affording fundamental frequencies for 7 different vibrational modes in the S₁(n,π*) state, including the C═C (1554 cm^-1) and OC-CH (1449 cm^-1) stretch modes. The C═O stretch character is spread over at least four different vibrational modes in the S₁(n,π*) state, with fundamentals spanning the 1340-1430 cm^-1 interval. This finding stems from a significant reduction in C═O bond order upon excitation, which leads to near-coincidence of the C═O stretch and several CH₂ wag frequencies. Such complexities make 2CHO an ideal candidate for testing excited-state computational methods. We have used the present spectroscopic results to test EOM-EE-CCSD harmonic-frequency predictions for the S₁(n,π*) state. We have also benchmarked the performance of less costly computational methods, including CIS(D) and TDDFT. For certain density functionals (e.g., B3LYP and PBE0), we find that the accuracy of TDDFT frequency predictions can approach but not meet that of EOM-EE-CCSD.

Lowest triplet (n,π*) state of 2-cyclohexen-1-one: characterization by cavity ringdown spectroscopy and quantum-chemical calculations

The cavity ringdown (CRD) absorption spectrum of 2-cyclohexen-1-one (2CHO) was recorded over the range 401.5-410.5 nm in a room-temperature gas cell. The very weak band system (ε ≤ 0.1 M(-1) cm(-1)) in this spectral region is due to the T1(n, π*) ← S0 electronic transition. The 0(0)(0) origin band was assigned to the feature observed at 24,558.8 ± 0.3 cm(-1). We have assigned 46 vibronic transitions in a region extending from -200 to +350 cm(-1) relative to the origin band. For the majority of these transitions, we have made corresponding assignments in the spectrum of the deuterated derivative 2CHO-2,6,6-d3. From the assignments, we determined fundamental frequencies for several vibrational modes in the T1(n, π*) excited state of 2CHO, including the lowest ring-twisting (99.6 cm(-1)) and ring-bending (262.2 cm(-1)) modes. These values compare to fundamentals of 122.2 cm(-1) and 251.9 cm(-1), respectively, determined previously for the isoconfigurational S1(n, π*) excited state of 2CHO and 99 cm(-1) and 248 cm(-1), respectively, for the S0 ground state. With the aid of quantum-mechanical calculations, we have also ascertained descriptions for these two modes, thereby resolving ambiguities appearing in the previous literature. The ring-twisting mode (ν39) contains a significant contribution from O=C-C=C torsion, whereas the ring-bending mode (ν38 in the ground state) involves mainly the motion of C-5 with respect to the plane containing the other heavy atoms. The CRD spectroscopic data for the T1(n, π*) state have allowed us to benchmark several computational methods for treating excited states, including time-dependent density functional theory and an equation-of-motion coupled cluster method. In turn, the computational results provide an explanation for observed differences in the T1(n, π*) vs. S1(n, π*) ring frequencies.

Lowest triplet (n, π*) electronic state of acrolein: determination of structural parameters by cavity ringdown spectroscopy and quantum-chemical methods

The cavity ringdown absorption spectrum of acrolein (propenal, CH(2)=CH-CH=O) was recorded near 412 nm, under bulk-gas conditions at room temperature and in a free-jet expansion. The measured spectral region includes the 0(0)(0) band of the T(1)(n, π*) ← S(0) system. We analyzed the 0(0)(0) rotational contour by using the STROTA computer program [R. H. Judge et al., J. Chem. Phys. 103, 5343 (1995)], which incorporates an asymmetric rotor Hamiltonian for simulating and fitting singlet-triplet spectra. We used the program to fit T(1)(n, π*) inertial constants to the room-temperature contour. The determined values (cm(-1)), with 2σ confidence intervals, are A = 1.662 ± 0.003, B = 0.1485 ± 0.0006, C = 0.1363 ± 0.0004. Linewidth analysis of the jet-cooled spectrum yielded a value of 14 ± 2 ps for the lifetime of isolated acrolein molecules in the T(1)(n, π*), v = 0 state. We discuss the observed lifetime in the context of previous computational work on acrolein photochemistry. The spectroscopically derived inertial constants for the T(1)(n, π*) state were used to benchmark a variety of computational methods. One focus was on complete active space methods, such as complete active space self-consistent field (CASSCF) and second-order perturbation theory with a CASSCF reference function (CASPT2), which are applicable to excited states. We also examined the equation-of-motion coupled-cluster and time-dependent density function theory excited-state methods, and finally unrestricted ground-state techniques, including unrestricted density functional theory and unrestricted coupled-cluster theory with single and double and perturbative triple excitations. For each of the above methods, we or others [O. S. Bokareva et al., Int. J. Quantum Chem. 108, 2719 (2008)] used a triple zeta-quality basis set to optimize the T(1)(n, π*) geometry of acrolein. We find that the multiconfigurational methods provide the best agreement with fitted inertial constants, while the economical unrestricted Perdew-Burke-Ernzerhof exchange-correlation hybrid functional (UPBE0) technique performs nearly as well.