Changfeng Fang,
Baswanth Oruganti and
Bo Durbeej*
Division of Computational Physics, IFM, Linköping University, SE-581 83 Linköping, Sweden. E-mail: bodur@ifm.liu.se
First published on 3rd February 2014
In recent years, much progress has been made in the design, synthesis and operation of light-driven rotary molecular motors based on chiral overcrowded alkenes. Through consecutive cis–trans photoisomerization and thermal helix inversion steps, where the latter dictate the overall rate of rotation, these motors achieve a full 360° unidirectional rotation around the carbon–carbon double bond connecting the two (rotator and stator) alkene halves. In this work, we report quantum chemical calculations indicating that a particularly fast-rotating overcrowded alkene-based motor capable of reaching the MHz regime, can be made to rotate even faster by the substitution of a rotator methyl group with a methoxy group. Specifically, using density functional theory methods that reproduce the rate-limiting ∼35 kJ mol−1 thermal free-energy barriers shown by the methyl-bearing motor with errors of ∼5 kJ mol−1 only, it is predicted that this substitution reduces these barriers by a significant 15–20 kJ mol−1. This prediction is preceded by a series of benchmark calculations for assessing how well density functional theory methods account for available experimental data (crystallographic, UV-vis absorption, thermodynamic) on the rotary cycles of overcrowded alkenes, and a detailed examination of the thermal and photochemical reaction mechanisms of the original motor of this type.
While Nature's biological machinery contains a number of complex protein assemblies that convert the energy stored in chemical bonds into directed rotary motion,11 such as ATP synthase,12 the first synthetic unidirectional rotary molecular motors were developed by Kelly13 and Feringa14,15 and their co-workers in the late nineties. Upon uptake of chemical and light energy, respectively, these systems produce motion consisting of rotation around a covalent bond: the former motor around a carbon–carbon single bond in a triptycene derivative13 and the latter around a carbon–carbon double bond in a sterically overcrowded alkene.15 For both motors, chirality is an essential feature for the unidirectional rotary motion.
Especially the light-driven design by Feringa has inspired a wealth of subsequent experimental research.6,16–27 The original biphenanthrylidene motor,15 shown in Scheme 1 and denoted 1, features two identical halves connected by a central carbon–carbon double bond. The configuration at the methyl-bearing stereogenic center is (R) for both halves. P denotes the right-handed helicity in each half of the motor, whereas M (used later) denotes left-handed helicity. The upper motor moiety is taken as the “rotator” that can rotate around the double-bond axle, whereas the lower “stator” moiety can be immobilized on a surface.28–31
![]() | ||
Scheme 1 Chemical structure and atom numbering of Feringa's first-generation (3R,3′R)-(P,P)-trans-1,1′,2,2′,3,3′,4,4′-octahydro-3,3′-dimethyl-4,4′-biphenanthrylidene rotary molecular motor 1. |
Each 360° rotation achieved by 1 involves four discrete steps activated by UV light or a change in temperature of the system.15 The overall process can be described as follows (Fig. 1).15,22 In the first step, irradiation of the (P,P)-trans-1 isomer with UV light (λ ≥ 280 nm) triggers a trans → cis photoisomerization around the central carbon–carbon double bond to produce the (M,M)-cis-2 isomer. This reaction occurs concomitantly with (P,P) → (M,M) inversion of the helicities of the two motor halves, and a change in orientation of the methyl substituents carried by the two stereogenic centers from favorable axial in (P,P)-trans-1 to strained equatorial in (M,M)-cis-2.
In (M,M)-cis-2, further isomerization continuing in the same direction as the initial trans → cis isomerization is prevented by steric hindrance. However, if the temperature is high enough, this hindrance can be overcome in an energetically downhill thermal (M,M) → (P,P) helix inversion process, in which the methyl groups regain their preferred axial orientation. This process constitutes the second step of the rotary cycle, produces the (P,P)-cis-2 isomer, and completes the first 180° of the rotation.
In the third step, irradiation of (P,P)-cis-2 with UV light (λ ≥ 280 nm) triggers a cis → trans photoisomerization that generates the (M,M)-trans-1 isomer, in which the methyl groups once more are forced to adopt a strained equatorial orientation. In further analogy with the first step, this reaction again changes the helicities of the motor halves from (P,P) to (M,M). Because of the preceding thermal step, the cis → trans photoisomerization can only proceed in the same direction as the initial trans → cis photoisomerization. Thus, these two photoisomerizations occur in a unidirectional fashion and produce truly rotary motion.
The fourth step, finally, is analogous to the second step in that it entails a spontaneous thermal (M,M) → (P,P) helix inversion and restores the preferred axial orientation for the methyl substituents. This step returns the system to the initial (P,P)-trans-1 state, thus completing the full 360° rotation and allowing for a new cycle to begin.
As is clear from the foregoing dissection, the rotary cycle of 1 proceeds through sequential photochemical and thermal steps. The unidirectionality of the two photoisomerizations relies on the steric hindrance that the methyl substituents at the stereogenic centers introduce in the so-called fjord regions (see Scheme 1), and on the exergonicity of the thermal helix inversion steps, which effectively block back rotations.15,22 Based on quantum chemical calculations, a similar ratchet-like mechanism was recently found to be in operation in a photosensory protein – Anabaena sensory rhodopsin – that naturally converts light energy into unidirectional rotary motion.32
From the viewpoint of future applications of synthetic rotary molecular motors, which may include rotation of objects much larger than the motors themselves33 and molecular transport,34 two key challenges are reaching high rates of rotation under ambient conditions22,35 and mounting the motors on surfaces.7,28–31 In the first of these respects, however, motor 1 is not an ideal system, because the rate-limiting thermal helix inversion steps are very slow for this molecule.22 A major experimental effort has therefore been invested in the development of second-generation overcrowded alkene-based rotary molecular motors capable of operating at higher rotational frequencies than 1.6,16–27 Indeed, by ingenious design and careful organic synthesis, it has been possible to lower the free-energy barriers of the thermal steps to such an extent that MHz frequencies are now realizable.25,27 One such motor,25 henceforth denoted 2, is shown in Scheme 2.
While it is clear that computational chemistry methods may fruitfully complement experimental endeavors to design faster-rotating molecular motors based on sterically overcrowded alkenes, no computational study in the existing literature seems to have been in-depth devoted to this particular topic, although a few other studies have used quantum chemical methods to put forward suggestions for improved light-driven rotary motors outside of the Feringa design.36–38 The present work is an attempt to help filling this important gap in the literature.
Among those other studies, Frutos and co-workers have recently proposed a motor that provides full 360° unidirectional rotation from a chiral hydrogen-bond environment without the intermediacy of thermal steps.38 Our investigation is also related to a number of previous computational studies of Feringa-type motors that have rather focused on the photochemical steps of the rotary cycle,39–42 or on the ground-state conformational dynamics of these systems.43
Specifically, this work reports quantum chemical calculations using density functional theory (DFT), including time-dependent density functional theory (TD-DFT)44–50 for the treatment of excited states, that illustrate the potential of such calculations to aid the development of faster-rotating overcrowded alkene-based molecular motors. First, we explore the rotary cycle of 1, including both photochemical and thermal steps, and find that a number of density functionals are able to demonstrate the unidirectionality of the two photoisomerizations, and to reproduce experimentally obtained structural, spectroscopic and thermodynamic data of the rotary cycle with reasonable accuracy. Furthermore, these calculations support experimental indications20 that the thermal helix inversion steps occur by a stepwise rather than concerted mechanism. Second, we consider the second-generation motor 2, and compute free-energy barriers for the thermal steps that closely reproduce the corresponding experimental estimates, and are consistent with the observation that this motor is able to reach much higher rotational frequencies than 1.25
Finally, having demonstrated the appropriateness of the computational methodology, we apply it to our own suggestion for a new second-generation motor (henceforth denoted 3 and also shown in Scheme 2) that is very similar to 2, but whose thermal barriers, we reasoned, should be lower by the replacement of the C5 methyl group (in 2) with a methoxy group (in 3). Thereby, we find that the thermal barriers for 3 are indeed a sizable 15–20 kJ mol−1 lower, and are thus able to predict that this molecule may well surpass the MHz rotational frequencies achieved by motor 2.25
To account for the fact that part of the experimental characterization of (P,P)-trans-1 and (P,P)-cis-2 was done in a hexane or methanol solution,14 and to assess the magnitude of solvent effects, the B3LYP and ωB97X-D optimizations were carried out both in the gas phase and using the integral equation formulation56 of the polarizable continuum model (PCM)57 to describe the solvent. For the resulting geometries, frequency calculations were performed at the corresponding levels of theory to ensure that these structures are potential-energy minima, as well as to derive Gibbs free energies (at room temperature).
Starting from the ground-state structures of (P,P)-trans-1 and (P,P)-cis-2, the photoisomerization paths in the lowest excited singlet state (S1) toward the (M,M)-cis-2 and (M,M)-trans-1 isomers were computed in the following way, considering both gas phase and solution environments. First, TD-DFT singlepoint calculations using ωB97X-D/SVP were carried out to obtain the vertically excited Franck–Condon (FC) point of the respective path. These calculations were followed by TD-DFT geometry optimizations, enabled by the availability of analytic TD-DFT gradients,58–61 at the same level of theory to model the subsequent relaxation from the FC points to excited-state minima henceforth denoted (P,P)-trans-1* and (P,P)-cis-2*. This relaxation initiates the photoisomerizations at the C4–C4′ double bond and defines a direction (clockwise or counterclockwise) for the rotary motion of 1 during its photocycle.
Starting from (P,P)-trans-1* and (P,P)-cis-2*, additional points along the photoisomerization paths were then obtained by performing a series of constrained TD-DFT optimizations, where for each optimization all other geometric parameters than the C4a–C4–C4′–C4′a dihedral angle (henceforth denoted α) were allowed to relax. The set of α values considered for these calculations cover, in steps of 10°, the full range (−180° to 180°) of possible dihedral angles. Based on a comparison of the FC relaxation in the gas phase and in solution (discussed in detail below) and for computational expedience, all of the constrained TD-DFT optimizations were carried out in the gas phase.
Based on earlier computational studies of photoisomerizations of organic molecules,62–67 the decay processes from the excited state to the ground state that precede the formation of the photoproducts [i.e., the (M,M)-cis-2 and (M,M)-trans-1 isomers] are expected to be mediated by conical intersections at highly twisted geometries where the two states are degenerate. However, since, despite recent progress,68–72 conical intersections are yet to lend themselves easily to location by means of DFT methods and since the multi-reference ab initio approaches best suited for such calculations (e.g., CASSCF73) remain too expensive for straightforward application to large conjugated systems, no attempt was made to demonstrate the occurrence of conical intersections along the current photoisomerization paths. Instead, the (M,M)-cis-2 and (M,M)-trans-1 isomers were located directly by starting ground-state ωB97X-D/SVP geometry optimizations at highly twisted (−110° ≤ α ≤ −70° and 70° ≤ α ≤ 110°, respectively) structures along the photoisomerization paths.
Importantly, irrespective of which twisted starting point was used from the respective path, these optimizations produced a number of identical structures of both (M,M)-cis-2 and (M,M)-trans-1. Furthermore, starting ground-state geometry optimizations from less twisted excited-state structures invariably returned the system to the parent (P,P)-trans-1 and (P,P)-cis-2 isomers. Thereby, it could be ascertained that the computed photoisomerization paths do indeed connect the parent isomers to the (M,M)-cis-2 and (M,M)-trans-1 isomers, which is important in light of the fact that the photochemical reaction coordinate (α) in our approach is presumed rather than defined by minimum energy path calculations.74
As for modeling the photoisomerization steps of the rotary cycle with a long-range-corrected hybrid functional (ωB97X-D), this choice is largely based on two findings. First, it has been reported that standard functionals like B3LYP may face problems in accounting for partial π-bond breaking during double-bond isomerization reactions.75 Second, it has been found that long-range-corrected hybrids, which in part have been designed to reduce well-known errors49,76 in TD-DFT excitation energies for states with appreciable charge-transfer character (by allowing the fraction of included exact exchange to vary with the interelectronic distance; between 22% at short range and 100% at long range for the case of ωB97X-D),54,77–81 also improve the TD-DFT description of photoisomerizations around carbon–carbon double bonds in conjugated systems.82,83 Specifically, studying the cis–trans photoisomerization of a retinal chromophore, Rostov et al. found that a number of long-range-corrected hybrids, including ωB97X-D, consistently produce more accurate excited-state PESs for this process than B3LYP.82,83 Another advantage of using ωB97X-D for the present molecules is that this functional includes dispersion,54 accounted for in the form of empirical atom–atom dispersion corrections.84,85
In addition to using ωB97X-D, parts (the FC relaxation) of the calculations for the S1 photoisomerization paths were also carried out with a set of complementary methods. This set included BP86, BLYP [pure density functionals based on the generalized gradient approximation (GGA)],51,86,87 B3LYP, PBE0 (global hybrid GGAs with 20 and 25% exact exchange, respectively),51–53,88 M06-2X (a global hybrid meta-GGA with 54% exact exchange),89 ωB97X (the long-range-corrected hybrid GGA from which ωB97X-D was developed),81 and the ab initio configuration interaction singles (CIS) method.90
Having located stationary points on the S0 PES of 1 corresponding to the four (P,P)-trans-1, (M,M)-cis-2, (P,P)-cis-2 and (M,M)-trans-1 isomers, the kinetics of the two thermal helix inversion (M,M)-cis-2 → (P,P)-cis-2 and (M,M)-trans-1 → (P,P)-trans-1 steps were investigated by mapping regions of the S0 PES encompassing both stepwise and concerted mechanisms. Specifically, continuing to use B3LYP and ωB97X-D in combination with the SVP basis set and examining the reactions in both the gas phase and in solution, this investigation involved location of transition structures and frequency and intrinsic reaction coordinate (IRC) calculations91 to verify that located transition structures do indeed mediate helix inversions.
Finally, as for the calculations on Feringa's second-generation rotary molecular motor 2 and our suggestion for a redesigned version 3 thereof, these were done with the same exact protocol as the calculations on motor 1 just described, but with the PCM settings modified to represent the dichloromethane solvent used in the experimental reference study.25
Except where otherwise noted, all calculations were performed with the Gaussian 09 suite of programs.92
Isomer | Method | Medium | C4–C4′ | Dihedral angleb | ΔG | ||||
---|---|---|---|---|---|---|---|---|---|
α | β | β′ | γ | γ′ | |||||
a Bond lengths in Å, dihedral angles in degrees, and energies in kJ mol−1.b With reference to Scheme 1, dihedral angles are defined as follows: α = C4a–C4–C4′–C4′a, β = C4–C4′–C4′a –C4′b, β′ = C4′–C4–C4a–C4b, γ = C4–C4′–C3′–C3′a, and γ′ = C4′–C4–C3–C3a.c Experimental data from ref. 14.d Experimental data not given in ref. 14.e Experimental data not available. | |||||||||
(P,P)-trans-1 | ωB97X-D | Gas phase | 1.349 | 171.0 | 68.1 | 68.1 | −121.3 | −121.3 | 3.9 |
B3LYP | Gas phase | 1.360 | 171.1 | 68.7 | 68.7 | −119.6 | −119.6 | 0.0 | |
ωB97X-D | Hexane | 1.350 | 170.6 | 67.9 | 67.9 | −121.8 | −121.8 | 9.5 | |
B3LYP | Hexane | 1.360 | 170.1 | 68.7 | 68.7 | −120.2 | −120.3 | 0.0 | |
ωB97X-D | Methanol | 1.349 | 170.7 | 68.1 | 68.1 | −121.7 | −121.7 | 16.8 | |
B3LYP | Methanol | 1.360 | 170.5 | 69.0 | 69.0 | −119.6 | −119.6 | 3.6 | |
Crystallographyc | 1.345 | 174.2 | 61.8 | 61.8 | –d | –d | –e | ||
(P,P)-cis-2 | ωB97X-D | Gas phase | 1.357 | −5.5 | 51.9 | 51.9 | −93.7 | −93.7 | 0.0 |
B3LYP | Gas phase | 1.364 | −0.1 | 56.1 | 56.1 | −95.5 | −95.5 | 10.6 | |
ωB97X-D | Hexane | 1.356 | −3.9 | 52.8 | 52.8 | −94.5 | −94.5 | 0.0 | |
B3LYP | Hexane | 1.365 | 0.1 | 55.8 | 55.8 | −95.3 | −95.3 | 3.2 | |
ωB97X-D | Methanol | 1.360 | −5.1 | 52.2 | 52.2 | −93.8 | −93.9 | 0.0 | |
B3LYP | Methanol | 1.365 | 0.3 | 56.2 | 56.2 | −95.3 | −95.3 | 0.0 | |
Crystallographyc | 1.347 | −3.2 | 54.4 | 54.4 | –d | –d | –e |
Starting with the comparison between the two functionals and focusing this analysis on the gas-phase results, which are affected by fewer computational factors than the PCM results, it should first be noted that such a comparison is particularly warranted since long-range-corrected functionals like ωB97X-D aimed especially at yielding accurate TD-DFT excitation energies have sometimes been found to accomplish this goal only at the expense of accuracy in ground-state properties.93 In this light, it is interesting to note that there is no substantial difference between the ωB97X-D and B3LYP geometries of (P,P)-trans-1 and (P,P)-cis-2, as indicated by the corresponding gas-phase parameters in Table 1 [and made clear by a comparison of the full geometries included in part 10 of the ESI†]. Furthermore, it is pleasing that both functionals reproduce the crystallographic data of ref. 14 quite well, with the calculated α and β dihedrals deviating by at most ∼7° from these data.
One result on which the functionals differ, however, is the relative free energies of the two isomers: ωB97X-D predicts (P,P)-cis-2 to be 3.9 kJ mol−1 more stable than (P,P)-trans-1 in the gas phase, whereas B3LYP conversely predicts (P,P)-cis-2 to be 10.6 kJ mol−1 less stable. Although smaller than the computational accuracy needed for the conclusions drawn in this work, this effect warrants further examination of how sensitive calculated free energies of the ground-state stationary points of 1 are to the choice of functional and basis set (beyond SVP). The results of such an investigation are summarized and discussed in the ESI (Tables S1 and S2†). From this, it can be inferred that these energies are not very sensitive to the level of theory, and argued that the fact that ωB97X-D and B3LYP slightly differ in some respects, yet without ambiguity support the same overall conclusions, helps solidifying these conclusions.
As for solvent effects, finally, the selected results in Table 1 are a clear reflection of what can be deduced more thoroughly from part 10 of the ESI:† the gas-phase geometries of the various isomers of 1 are very similar to the ones obtained by performing geometry optimizations in an experimentally relevant14 solvent. This is likely a consequence of 1 being a neutral species. The solvent exerts a somewhat more noticeable effect on the calculated free energies, but not to an extent that its presence appears to be of qualitative importance. For (P,P)-trans-1 and (P,P)-cis-2, the solvent stabilizes the latter more than the former, and increasingly so when going from hexane (5.6–7.4 kJ mol−1) to methanol (12.9–14.2 kJ mol−1).
Isomer | Energies | Geometric parameters (S0 → S1) | |||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
VEE | AEE | λmaxb | C4–C4′ | α | β | β′ | γ | γ′ | |||||||
a Excitation energies in eV, bond lengths in Å, and dihedral angles in degrees.b Experimental absorption maximum (in eV) from ref. 14.c Experimental data not available. | |||||||||||||||
(P,P)-trans-1 | |||||||||||||||
Gas phase | 4.33 | 3.26 | —c | 1.349 | 1.456 | 171.0 | −97.3 | 68.1 | 12.2 | 68.1 | 12.3 | −121.3 | −89.1 | −121.3 | −89.2 |
Hexane | 4.26 | 3.18 | —c | 1.350 | 1.460 | 170.6 | −97.4 | 67.9 | 14.0 | 67.9 | 14.0 | −121.8 | −91.4 | −121.8 | −90.8 |
Methanol | 4.26 | 3.09 | 3.80 | 1.349 | 1.465 | 170.7 | −111.8 | 68.1 | 20.4 | 68.1 | 20.4 | −121.7 | −96.8 | −121.7 | −96.8 |
(P,P)-cis-2 | |||||||||||||||
Gas phase | 4.08 | 3.33 | —c | 1.357 | 1.427 | −5.5 | 11.8 | 51.9 | 27.3 | 51.9 | 34.9 | −93.7 | −108.9 | −93.7 | −84.6 |
Hexane | 4.08 | 3.34 | 4.11 | 1.356 | 1.430 | −3.9 | 12.6 | 52.8 | 26.6 | 52.8 | 35.0 | −94.5 | −107.6 | −94.5 | −84.8 |
Methanol | 4.03 | 3.34 | —c | 1.360 | 1.438 | −5.1 | 16.6 | 52.2 | 24.7 | 52.2 | 34.1 | −93.8 | −105.5 | −93.9 | −86.5 |
As can be seen from Table 2, the ωB97X-D vertical excitation energies are in reasonable agreement with the experimental absorption maxima, and are not appreciably shifted between the different environments. In contrast to the experiments, however, these excitation energies do not yield a blue shift for the absorption of (P,P)-cis-2 relative to (P,P)-trans-1, but a red shift. A similar result has been reported in a previous computational study40 using the BH&HLYP global hybrid functional (that incorporates 50% exact exchange)94 in the framework of TD-DFT, as well as in combination with a state-averaged variant of the spin-restricted ensemble-referenced Kohn–Sham method.95
One possible explanation for this discrepancy is that geometric relaxation effects weaken the (standard) assumption that experimental absorption maxima correspond to vertical transitions (the vertical approximation), because the calculated adiabatic excitation energies in Table 2 do yield a blue shift for (P,P)-cis-2 relative to (P,P)-trans-1. Another possibility is simply that the ωB97X-D vertical excitation energies are erroneous in this particular regard. Testing whether other methods would perform differently, Table 3 lists vertical and adiabatic excitation energies of the two isomers obtained with six alternative density functionals (BP86, BLYP, B3LYP, PBE0, M06-2X and ωB97X) and CIS. Based on the ωB97X-D results, these calculations were carried out in the gas phase. Interestingly, analogously to ωB97X-D, all of these methods yield a red shift for the absorption of (P,P)-cis-2 relative to (P,P)-trans-1 if vertical excitation energies are considered. As far as this test is concerned, then, there seems to be no obvious reason to use another functional than ωB97X-D for the photochemical modeling. Further support for this thesis comes from complementary calculations summarized in Table S3 of the ESI,† whereby also the ab initio approximate coupled-cluster singles and doubles (CC2)96 method red-shifts the vertical excitation energy of (P,P)-cis-2 relative to (P,P)-trans-1, and whereby it is shown that the ωB97X-D results are not changed much by the use of a larger basis set than SVP.
Isomer | Method | Energies | Geometric parameters (S0 → S1) | ||||||
---|---|---|---|---|---|---|---|---|---|
VEE | AEE | C4–C4′ | ΔC4–C4′ | α | Δα | ||||
a Excitation energies in eV, bond lengths in Å, and dihedral angles in degrees.b Experimental absorption maximum (in eV) from ref. 14. | |||||||||
(P,P)-trans-1 | BP86 | 3.03 | 2.84 | 1.372 | 1.406 | 0.034 | 171.0 | 177.5 | 6.5 |
BLYP | 3.02 | 2.84 | 1.374 | 1.404 | 0.030 | 171.4 | 176.7 | 5.3 | |
B3LYP | 3.70 | 3.11 | 1.360 | 1.433 | 0.073 | 171.1 | −162.5 | 26.4 | |
PBE0 | 3.84 | 3.53 | 1.356 | 1.409 | 0.053 | 170.9 | 179.7 | 8.8 | |
M06-2X | 4.30 | 3.15 | 1.350 | 1.451 | 0.101 | 171.2 | −97.9 | 90.9 | |
ωB97X | 4.52 | 3.42 | 1.350 | 1.458 | 0.108 | 170.3 | −99.3 | 90.4 | |
ωB97X-D | 4.33 | 3.26 | 1.349 | 1.456 | 0.107 | 171.0 | −97.3 | 91.7 | |
CIS | 5.06 | 4.36 | 1.339 | 1.446 | 0.107 | 170.6 | −139.0 | 50.4 | |
Exp. λmaxb | 3.80 | — | — | — | — | — | — | ||
(P,P)-cis-2 | BP86 | 2.98 | 2.73 | 1.377 | 1.428 | 0.051 | 0.7 | 2.6 | 1.9 |
BLYP | 2.97 | 2.79 | 1.378 | 1.418 | 0.040 | 0.6 | 8.2 | 7.6 | |
B3LYP | 3.60 | 3.13 | 1.364 | 1.430 | 0.066 | −0.1 | 11.1 | 11.2 | |
PBE0 | 3.71 | 3.16 | 1.361 | 1.423 | 0.062 | −1.1 | 3.7 | 3.16 | |
M06-2X | 4.00 | 3.31 | 1.357 | 1.423 | 0.066 | −5.5 | 5.7 | 11.2 | |
ωB97X | 4.38 | 3.68 | 1.355 | 1.428 | 0.073 | −3.0 | 6.0 | 9.0 | |
ωB97X-D | 4.08 | 3.33 | 1.357 | 1.427 | 0.070 | −5.5 | 11.8 | 17.3 | |
CIS | 4.97 | 4.49 | 1.344 | 1.419 | 0.075 | −3.1 | 6.5 | 9.6 | |
Exp. λmaxb | 4.11 | — | — | — | — | — | — |
Turning to the structural changes during the FC relaxation and considering first the ωB97X-D data in Table 2, a pronounced lengthening of the central C4–C4′ bond is observed, that ranges from 0.07–0.08 Å for (P,P)-cis-2 to 0.11–0.12 Å for (P,P)-trans-1 and shows no particular sensitivity to the presence of a solvent environment. Such bond stretching facilitates photoisomerization and is for both isomers accompanied by substantial torsional motion along the α photoisomerization coordinate. Indeed, during the FC relaxation of (P,P)-trans-1, the upper motor part is rotated around the C4–C4′ bond by more than 90°, which suggests that much of the photochemical motion of this isomer proceeds without an energy barrier to be overcome. In this connection, fluorescence up-conversion measurements predict that the photochemical processes of this type of molecular motors occur within hundreds of fs, i.e., essentially in a barrierless fashion.97 For (P,P)-cis-2, in turn, the upper motor part is rotated by 17–22°. Notably, from the associated changes in the α dihedral angles, it can be deduced that the direction of photoinduced torsional motion is the same – counterclockwise – for both isomers, which means that the two photoisomerizations occur in a unidirectional fashion and produce rotary motion during the reaction cycle of 1.
From Table 3 and the calculations with other density functionals than ωB97X-D, it is found that these support exactly the same conclusions on the FC relaxation as ωB97X-D, albeit with some quantitative differences between them as to the extent of the photoinduced rotation around the C4–C4′ bond. This sensitivity to the choice of method is likely to be a consequence of the flatness of the S1 PES along the α coordinate.
In Fig. 2, all ωB97X-D S1 data points calculated along the photoisomerization coordinate are shown. Based on the foregoing analysis, these calculations were performed in the gas phase. While the full range of possible α dihedral angles are considered [−180° to 0° for (P,P)-trans-1; 0° to 180° for (P,P)-cis-2], it should be pointed out that not all of them are photochemically relevant and that, as discussed in Section 2, the photoisomerizations are likely to involve conical intersections at highly twisted geometries where decay to the ground state takes place. Furthermore, although ωB97X-D has been found to perform well for photoisomerizations around carbon–carbon double bonds in conjugated systems,82,83 it is a single-reference method and as such is less well suited for explicit location of conical intersections than multi-reference methods like CASSCF. Instead, as outlined in Section 2, the (M,M)-cis-2 and (M,M)-trans-1 species whose formation is mediated by the presumed conical intersections were located by starting ground-state geometry optimizations at highly twisted excited-state structures along the photoisomerization paths.
![]() | ||
Fig. 2 ωB97X-D gas-phase photoisomerization paths for (P,P)-trans-1 and (P,P)-cis-2 (the direction of photoisomerization is indicated by arrows). |
The key feature of the calculated paths in Fig. 2 is their flatness, which is indicative of the extraordinarily fast excited-state dynamics that these systems exhibit.97 Furthermore, following population of the vertically excited FC point, there is a considerable driving force to reach the respective excited-state minimum: 103.2 kJ mol−1 for (P,P)-trans-1* and 72.4 kJ mol−1 for (P,P)-cis-2*. Loosely, the somewhat more favorable – both in terms of driving force and kinetics – energetics for the (P,P)-trans-1 reaction is consistent with the experimentally observed photoequilibrium “product-over-reactant” ratios for (M,M)-cis-2 to (P,P)-trans-1 being higher (95:
5) than that for (M,M)-trans-1 to (P,P)-cis-2 (90
:
10).14,15
![]() | ||
Fig. 3 Tentative mechanism for the rotary cycle of motor 1 with the thermal helix inversions proceeding in a stepwise fashion. |
Some mechanistic aspects of the thermal helix inversions have been explored in previous computational studies. For example, in a gas-phase study focusing primarily on modeling the photochemical steps of 1 by means of Car–Parrinello molecular dynamics simulations but investigating also the thermal reactions using static quantum chemical calculations (the thermal reactions are much too slow to be modeled by Carr–Parinello simulations), Grimm et al. considered a stepwise mechanism and explored the relevant parts of the S0 PES at the semiempirical AM1 level of theory.39 While these researchers did not explicitly locate transition structures for the thermal steps, as we have done in the present work, they performed a series of constrained geometry optimizations to obtain a two-dimensional PES with respect to torsional motion along the β and β′ coordinates defined in Table 1. Thereby, they estimated that the potential energy (thus neglecting zero-point vibrational and thermal corrections to the potential energy) barriers for steps 2 and 3 in Fig. 3 amount to 69–88 and 50–59 kJ mol−1, respectively.39 As a comparison, the experimental study by ter Wiel et al. found that the (M,M)-cis-2 → (P,P)-cis-2 conversion (i.e., steps 2 + 3) of 1 in hexane has an overall free-energy barrier of 91 kJ mol−1.20 For the (M,M)-trans-1 → (P,P)-trans-1 reaction, in turn, they reported an overall barrier of 107 kJ mol−1.20 In another computational study, Pérez-Hernández and González carried out an exhaustive Monte Carlo-like conformational search for a second-generation motor and a redesigned version thereof.43
In this work, we first embarked on locating all of the S0 stationary points of the stepwise reaction in Fig. 3 by performing all of the requisite calculations (geometry optimizations and frequency and IRC calculations) using both ωB97X-D and B3LYP, and considering both gas phase and solution environments. Besides the parent (P,P)-trans-1 and (P,P)-cis-2 isomers, the photoisomerized (M,M)-cis-2 and (M,M)-trans-1 isomers, and the (M,P)-cis-2 and (M,P)-trans-1 intermediates, these stationary points also include the relevant transition structures, which are denoted TS2 (step 2), TS3 (step 3), TS5 (step 5) and TS6 (step 6). The results of these calculations are presented in Fig. 4. Since ωB97X-D and B3LYP again were found to produce similar energetics and solvent effects again were found to be small, Fig. 4 shows only the ωB97X-D gas-phase results, together with the corresponding results for the photochemical steps. The results obtained at the other levels of theory are included in Table S4 of the ESI.†
Starting from the photoisomerized (M,M)-cis-2 isomer in the first half of the rotary cycle, the first helix inversion to form the (M,P)-cis-2 intermediate via TS2 is estimated to have a free-energy barrier of ∼101 kJ mol−1, which is higher than the ∼83 kJ mol−1 barrier for the subsequent helix inversion that, via TS3, brings the motor to the (P,P)-cis-2 station and completes the first half of the rotary cycle. Thus, for the (M,M)-cis-2 → (M,P)-cis-2 → (P,P)-cis-2 process, the calculations predict that the first step is the rate-determining one, and ascribe to this step a free-energy barrier (∼101 kJ mol−1) that agrees well with the aforementioned experimental barrier of 91 kJ mol−1 between (M,M)-cis-2 and (P,P)-cis-2.20
As for the calculated free-energy lowering upon the change in orientation of the C3 and C3′ methyl groups from strained equatorial in (M,M)-cis-2 to favorable axial in (P,P)-cis-2, this amounts to ∼55 kJ mol−1. Such a sizable driving force guarantees that this reaction occurs spontaneously, and plays a key role for the unidirectional motion of the motor by rendering the (P,P)-cis-2 → (M,M)-cis-2 back reaction thermodynamically unviable.22 Thereby, the rotary cycle becomes irreversible.
While, as noted above, experimental data for an i-propyl substituted motor has been taken as evidence for a stepwise mechanism for the thermal helix inversions,20 there is to the best of our knowledge no experimental data on how much higher in energy the (M,P)-cis-2 intermediate lies than the (P,P)-cis-2 isomer, that can serve as reference for the calculated value of ∼68 kJ mol−1. However, this value is in very good accord with a previous DFT-based estimate of 63–65 kJ mol−1.39
Turning to the second half of the rotary cycle in Fig. 4 and the thermal transformation of the photoisomerized (M,M)-trans-1 isomer into (P,P)-trans-1, there are both differences and similarities between this free-energy profile and that of the first half cycle just described. As for differences, the initial helix inversion that forms the (M,P)-trans-1 intermediate via TS5 is estimated to have a lower free-energy barrier (∼124 kJ mol−1) than the subsequent helix inversion that completes the full rotary cycle via TS6 (∼136 kJ mol−1). Thus, for this half cycle, the calculations suggest that the second helix inversion is the rate-determining step. The ∼136 kJ mol−1 free-energy barrier associated with this step is in reasonable agreement with the ∼107 kJ mol−1 barrier predicted by experiments for the overall (M,M)-trans-1 → (P,P)-trans-1 process.20
As for similarities between the two half cycles, the calculated ∼46 kJ mol−1 free-energy lowering from the two thermal steps of the second half cycle compares well with the corresponding value of ∼55 kJ mol−1 for the two thermal steps of the first half cycle. Furthermore, the calculated ∼26 kJ mol−1 free-energy difference between the (M,P)-trans-1 intermediate and the (P,P)-trans-1 isomer agrees just as well with what was reported in the DFT-based study referred to above (28–31 kJ mol−1),39 as does the calculated free-energy difference between the (M,P)-cis-2 intermediate and the (P,P)-cis-2 isomer.
Besides yielding rate-determining free-energy barriers of ∼101 and ∼136 kJ mol−1 for the first and second half cycles that agree well and reasonably, respectively, with the experimental values of ∼91 and ∼107 kJ mol−1,20 it is clear that the ωB97X-D calculations also reinforce the experimental prediction that the second half cycle is slower than the first.20 This conclusion is corroborated by the calculation of all ground-state stationary points of motor 1 using a range of alternative density functionals, as detailed in Table S2 of the ESI.† Indeed, similarly to ωB97X-D, all of these methods give higher barriers for the helix inversions of the second half cycle. Moreover, in further agreement with the ωB97X-D results, they uniformly predict that the first helix inversion is rate-determining in the first half cycle, and that the second helix inversion is rate-determining in the second half cycle.
Interestingly, for the i-propyl-substituted motor for which the stepwise mechanism explored by the preceding calculations was originally implicated,20 it was possible to obtain experimental estimates of the free-energy barriers for the individual helix inversions of the second half cycle that, in line with our calculated results for motor 1, indicate that the second inversion (barrier of 131 kJ mol−1) is slower than the first (barrier of 124 kJ mol−1).20 Hence, it is of interest to test whether our computational approach can reproduce also this finding. Such calculations on the i-propyl-substituted motor are summarized in Table S5 of the ESI.† Pleasingly, it is found that this is indeed the case, albeit that the 124 kJ mol−1 barrier for the (M,M)-trans-1 → (M,P)-trans-1 inversion is underestimated by 16–28 kJ mol−1 and the 131 kJ mol−1 barrier for the (M,P)-trans-1 → (P,P)-trans-1 inversion is overestimated by 13–29 kJ mol−1.
Finally, as for the concerted mechanism, our calculations did provide two indications supporting the experimental view that this mechanism is less relevant than the stepwise mechanism,20 albeit not in the form of explicitly located concerted transition structures. The first indication comes from the two-dimensional ωB97X-D and B3LYP S0 PESs for torsional motion along the β and β′ coordinates in Fig. S1 of the ESI,† which clearly suggest that the energy barriers for a concerted process are higher than those for a stepwise process. The second, related, indication is the observation that any attempt to optimize a concerted transition structure starting from a structure “close” to the highest-energy point along a tentative concerted reaction path in Fig. S1,† invariably converged to a transition structure along the stepwise reaction path. Thus, this path appears energetically preferable over the concerted ditto.
In this section, we will further illustrate this point by reporting calculations on the second-generation overcrowded alkene-based rotary motor 2 shown in Scheme 2, whose thermal barriers are much lower than those of 1 and which is able to sustain MHz rotational frequencies at ambient temperatures.25 Moreover, we will also report calculations on our own suggestion 3 for a slightly modified version of motor 2 that, we reasoned, should exhibit even lower thermal barriers. Interestingly, it will be inferred from these calculations that this prediction appears to be correct. As can be seen in Scheme 2, the only difference between motors 2 and 3 is that the C5 methyl group of 2 is replaced by a methoxy group in 3. Although introducing a group of similar size, the idea was that such a substitution should nonetheless reduce the steric hindrance to be overcome during the thermal steps, by virtue of the methoxy group being positioned further away from the stator than the methyl group.
In second-generation rotary motors, a number of different approaches are utilized to improve the kinetics of the thermal helix inversion steps,6,16–27 including some of which are implemented in motor 2. For example, by contracting the six-membered rings fused to the rotation axis into five-membered rings, or by replacing the naphthalene moieties with phenyl groups, the steric demands on the helix inversions can be mitigated. An analogous effect can be achieved by variation of the stereogenic substituents, which may also alter the thermal barriers by stabilizing or destabilizing some isomers of the rotary cycle more than others. It is also possible to influence the thermal barriers through the introduction of different bridging rotator and/or stator atoms (e.g., the sulfur atom in 2) with varying electron-donating/withdrawing capabilities.
The key results from the calculations on motors 2 and 3 are summarized in Table 4. In the interest of brevity, this summary focuses exclusively on the rate-determining thermal barriers in the rotary half cycles, where trans → cis denotes the half cycle initiated by light absorption of the parent trans isomer, and cis → trans denotes the half cycle initiated by light absorption of the parent cis isomer. Results from excited-state calculations indicating that the photoisomerizations of the two parent isomers occur in a unidirectional fashion, and thus produce rotary motion, are included in Fig. S2 of the ESI.†
Motor | Half cycle | ωB97X-D | B3LYP | Exp.b | ||
---|---|---|---|---|---|---|
Gasc | Sol.d | Gasc | Sol.d | Sol.d | ||
a Similar results obtained with larger basis sets than SVP are presented in Table S6 of the ESI.b Experimental data from ref. 25.c Gas-phase environment.d Dichloromethane solvent environment.e Experimental data not available. | ||||||
2 | trans → cis | 41.4 | 37.1 | 34.8 | 30.8 | 35.0 |
cis → trans | 40.0 | 37.1 | 35.2 | 30.0 | 34.2 | |
3 | trans → cis | 23.9 | 17.1 | 19.7 | 12.6 | —e |
cis → trans | 22.9 | 15.2 | 17.4 | 10.4 | —e |
From Table 4, one first notes that the calculations reproduce the experimental observation25 that motor 2 is able to reach much higher rotational frequencies than the original motor 1. Specifically, for both half cycles, all calculated estimates are within ∼5 kJ mol−1 agreement with the observed rate-determining free-energy barrier, which for both half cycles is appreciably lower (∼35 kJ mol−1) than what experiments20 and our foregoing calculations assigned to 1. The fact that the barriers are virtually identical in the two half cycles, whereas for 1 they were quite different, is easily understood through inspection of the optimized reactant and transition structures of motor 2 in Fig. 5. Indeed, from these it is clear that the stator methoxy group is sterically inactive not only in the reactants that precede the thermal helix inversions (shown in red color), but also in the rate-determining transition structures (shown in green color). Hence, while this methoxy group makes the stator asymmetric, the stator is anyhow “symmetric” in terms of how its terminal phenyl moieties interact sterically. As a consequence, the thermal barriers in one of the half cycles are very similar to those in the other. For motor 1, on the other hand, the steric interactions are different in the two half cycles, whereby the barriers are different too.
![]() | ||
Fig. 5 Optimized ωB97X-D gas-phase reactant (in red color) and transition (in green color) structures for the rate-determining thermal helix inversion steps of motors 2 and 3. |
Finally, as for the calculations on motor 3, the results in Table 4 do predict, as discussed above, a lowering of the rate-determining free-energy barriers upon replacement of the C5 methyl group (in 2) with a methoxy group (in 3). In fact, all levels of theory estimate that the barriers for 3 are 15–20 kJ mol−1 lower. Given the small errors with which the calculations on 2 reproduce the corresponding experimental barriers, we believe it is well-founded to conclude that 3 is a promising candidate for an overcrowded alkene-based rotary molecular motor capable of surpassing 2 in rotational frequency.
Exploring the rotary cycle of 1 for benchmark purposes and employing ωB97X-D and B3LYP as “workhorse” methods (but using also a number of alternative density functionals for some of the calculations), it is first found that these methods are able to reproduce crystallographic14 and UV-vis absorption data14 for the parent (P,P)-trans-1 and (P,P)-cis-2 isomers with reasonable accuracy. Furthermore, the distinctly directional light-induced torsional motion of (P,P)-trans-1 and (P,P)-cis-2 with respect to the central C4–C4′ double bond shown by the calculations is consistent with the fact that the photoisomerizations of these isomers afford a full 360° unidirectional rotation.
As for the thermal helix inversion steps, which are clearly what limit the overall rotation rates attainable by overcrowded alkene-based motors97 and hold key to their rate acceleration, the calculations on motor 1 support and complement available experimental data20 on the preference of a stepwise mechanism over a concerted ditto. Moreover, the calculations predict that the first helix inversion is rate-determining in the (P,P)-trans-1 → (P,P)-cis-2 rotary half cycle and the second helix inversion in the (P,P)-cis-2 → (P,P)-trans-1 half cycle, and provide free-energy barriers that are in qualitative accordance with the corresponding experimental estimates.20
Second, investigating the MHz-capable25 second-generation motor 2 to further assess the merits of the computational methodology, it is found that both ωB97X-D and B3LYP can accurately reproduce the much lower rate-determining thermal barriers of around ∼35 kJ mol−1 ascribed to this rotary cycle. Indeed, the calculated barriers agree to within ∼5 kJ mol−1 with the experimental values,25 although partly because of a fortuitous cancellation of errors.
Finally, applying the computational methodology to our own suggestion for a new second-generation motor (motor 3) hypothesized to exhibit even lower thermal barriers than 2 because of the substitution of the rotator methyl group with a methoxy group, this hypothesis does indeed seem correct. Specifically, it is predicted that such a substitution would lower the rate-determining barriers by a further 15–20 kJ mol−1. Hence, if readily synthesizable, 3 appears to be a potent candidate for a fast-rotating molecular motor.
Footnote |
† Electronic supplementary information (ESI) available: Cartesian coordinates for stationary points, Tables S1–S6, and Fig. S1 and S2. See DOI: 10.1039/c3ra46880a |
This journal is © The Royal Society of Chemistry 2014 |