Next Article in Journal
Polyploidy and Myc Proto-Oncogenes Promote Stress Adaptation via Epigenetic Plasticity and Gene Regulatory Network Rewiring
Previous Article in Journal
Comparative Proteomic Analysis of Liver Tissues and Serum in db/db Mice
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 17
10.3390/ijms23179694
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor

by
Jianzheng Ma
1,
Di Zhao
1,
Chenwei Jiang
1,*,
Zhenggang Lan
2,* and
Fuli Li
1
1
Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi’an Jiaotong University, Xi’an 710049, China
2
Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(17), 9694; https://doi.org/10.3390/ijms23179694
Submission received: 31 July 2022 / Revised: 20 August 2022 / Accepted: 23 August 2022 / Published: 26 August 2022
(This article belongs to the Section Molecular Biophysics)
Browse Figures
Versions Notes

Abstract

:
The working mechanism of conventional light-driven molecular rotary motors, especially Feringa-type motors, contains two photoisomerization steps and two thermal helix inversion steps. Due to the existence of a thermal helix inversion step, both the ability to work at lower temperatures and the rotation speed are limited. In this work, a two-stroke light-driven molecular rotary motor, 2-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDPY), is proposed, which is capable of performing unidirectional and repetitive rotation by only two photoisomerization (EP→ZP and ZP→EP) steps. With trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EP nonadiabatic dynamics of DDPY were systematically studied at different temperatures. Both EP→ZP and ZP→EP photoisomerizations are on an ultrafast timescale (ca. 200–300 fs). The decay mode of EP→ZP photoisomerization is approximately bi-exponential, while that of ZP→EP photoisomerization is found to be periodic. For EP and ZP isomers of DDPY, after the S 0 →S 1 excitation, the dynamical processes of nonadiabatic decay are both followed by twisting about the central C=C double bond and the pyramidalization of the C atom at the stator-axle linkage. The effect of temperature on the nonadiabatic dynamics of EP→ZP and ZP→EP photoisomerizations of DDPY has been systematically investigated. The average lifetimes of the S 1 excited state and quantum yields for both EP→ZP and ZP→EP photoisomerization are almost temperature-independent, while the corresponding unidirectionality of rotation is significantly increased (e.g., 74% for EP→ZP and 72% for ZP→EP at 300 K vs 100% for EP→ZP and 94% for ZP→EP at 50 K) with lowering the temperature.

Graphical Abstract

1. Introduction

Molecular machines [1,2] are a kind of nanomachine that can perform controllable and continuous changes of molecular structure under the stimulation of external energy source, so as to complete specific tasks. Such archetypal machines are molecular motors [3,4,5,6,7,8] which are capable of performing unidirectional and repetitive rotation at the nanoscale upon the stimulation of external energy. Light-driven molecular rotary motors (LDMRMs) [2,5,6], which can utilize the photoinduced E-Z (or cis-trans) isomerization of the double bond (to date, carbon-carbon double bond [9,10,11,12,13,14,15] or carbon-nitrogen double bond [16,17,18]) to complete a full 360 ∘ rotation by the absorption of UV or visible light, have attracted considerable interest in recent decades because of the high efficiency and cleanliness.
The complete 360 ∘ rotation of conventional LDMRMs, especially Feringa-type motors [9,10,11,12,13], is achieved through two steps of photoisomerization combined with two steps of thermal helix inversion (THI). The timescale of the THI steps (from nanoseconds to days) is much longer than that of photoisomerization steps (from femtoseconds to picoseconds) [6]. Thus, reducing the THI steps can not only improve the rotation speed but also enable LDMRM to operate at lower temperatures [6,19,20]. Some developments have been achieved on reducing the THI steps theoretically [21,22,23,24,25] and experimentally [15,26] in recent years.
From a theoretical aspect, by generating a chiral hydrogen bond environment, a two-stroke LDMRM which involves only two steps of photoisomerization was proposed by García-Iriepa et al. [21]. By redesigning a bio-inspired 4-hydroxybenzylidene-1,2-dimethylimidazolinone-based molecular photoswitch [27], Filatov et al. [22] proposed a family of two-stroke photon-only LDMRMs recently. With the nonadiabatic molecular dynamics (NAMD) simulations at the SSR-BH&HLYP/6-31G(d) level, these molecular motors are predicted to have very high quantum yields (about 0.91–0.97) and a sufficiently high degree (0.94–1.00) of unidirectionality. A visible-light responsive Schiff-based LDMRM, which is able to complete a 360 ∘ unidirectional rotation by only two photoisomerization steps, was designed by Wang et al. [23] recently. The quantum yields were predicted to be almost 70% for its individual E to Z and Z to E photoisomerizations using NAMD simulations. A novel molecular motor in which the rotation is induced by the electric coupling of chromophores was suggested by Majumdar et al. [24] recently, which was predicted to achieve unidirectional rotation on a subnanosecond time scale using the power of a single photon. A three-stroke LDMRM, 2-(2,7-dimethyl-2,3-dihydro-1Hinden-1-ylidene)-1,2-dihydro-3H-pyrro- l-3-one (DDIY), was proposed by our group [25] very recently, which is capable of completing a unidirectional rotation by two photoisomerization steps and one thermal helix inversion step at room temperature.
On the experimental side, Gerwien et al. [15] designed a three-stroke photon-only hemithi-oindigo-based molecular motor recently, which interconverts three different isomeric states in a fixed sequence upon visible light irradiation, without thermal ratcheting in the ground state. Three new second-generation molecular motors featuring a phosphorus center in the lower half have been reported by Boursalian et al. [26] recently. Four diastereomeric states of these molecular motors can interconvert solely photochemically. All-photochemical unidirectional rotation of the new molecular motors was confirmed by kinetic analysis and modeling.
Although some excellent developments on photon-only molecular motors [21,22,23] have been achieved from a computational perspective, Feringa et al. [6] pointed out very recently that synthesizing of these molecular motors in experiments is often highly challenging. Designing a light-driven molecular motor with fewer operational steps based on synthesized molecular systems may be an effective way. Based on some easy-to-synthesize oxindole-based molecular motors studied by Roke et al. [12] and Pooler et al. [13] recently, by further reducing the steric hindrance in the fjord region of the LDMRM DDIY proposed in our previous work [25], a new two-stroke LDMRM, 2-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDPY), was designed, in which only two photoisomerization steps (EP→ZP and ZP→EP) are involved to complete a full 360 ∘ rotation. The photoinduced isomerization dynamics of this two-stroke LDMRM were systematically investigated with trajectory surface-hopping molecular dynamics at the semi-empirical OM2/MRCI level.
It is well known that temperature plays a very crucial role in chemical reactions, e.g., affecting the rate and direction of a chemical reaction. The effect of temperature on the thermal helix inversion steps of LDMRMs has been extensively studied [5,6]. However, to the best of our knowledge, much less is known about the effect of temperature on the photoisomerization processes of LDMRMs, especially the unidirectionality [28] of LDMRMs. Very recently, Gerwien et al. [15] found that the monodirectionality [28] and quantum yield increase with lowering the temperature for a three-step photon-only LDMRM. To figure out the effect of temperature on photoisomerization dynamics of LDMRM DDPY, non-adiabatic dynamics simulations at different temperatures were systemically performed. The dynamic results show that, as the temperature decreases, both the average S 1 lifetime and quantum yield of EP→ZP and ZP→EP photoisomerization processes of LDMRM DDPY are almost temperature-independent, while the corresponding unidirectionality increases significantly, e.g., 74% for EP→ZP and 72% for ZP→EP at 300 K vs. 100% for EP→ZP and 94% for ZP→EP at 50 K.

2. Methods and Materials

2.1. Density Functional Methods

The geometrical optimization and frequency calculation of ground state and transition state of molecular motor DDPY by density functional theory (DFT) were performed at the B3LYP/6-31G(d), CAM-B3LYP/6-31G(d), and B3LYP-D3/6-31+G(d) levels. The GAUSSIAN 09 program [29] was utilized to carry out all the DFT calculations.

2.2. Semiempirical Methods

The OM2/MRCI method implemented in the development version of the MNDO program [30] was utilized to perform all the semi-empirical calculations. This method can balance the computational cost and accuracy well, as confirmed by many benchmark calculations [31,32,33], and has been applied to investigate many photoinduced processes [34,35,36,37,38,39,40,41,42,43,44,45] successfully.
For geometry optimizations and dynamics simulations, all required energies, gradients, and nonadiabatic coupling elements were computed analytically. The self-consistent field (SCF) calculations were performed in the restricted open-shell Hartree–Fock (ROHF) formalism, as it provided a better description of the excited-state wave functions. For the multireference configuration interaction (MRCI) treatment, three reference configurations were chosen, which includes the closed-shell ground-state configuration and single and double excitations from the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular orbital (LUMO). The active space in the MRCI calculations included 10 electrons in nine orbitals, which comprises four highest doubly occupied orbitals, two singly occupied orbitals, and three lowest unoccupied orbitals. The Lagrangian-Newton approach [46] was used to locate the S 1 /S 0 minimum-energy conical intersections (CIs) geometries.
The nonadiabatic photoisomerization dynamics of molecular motor DDPY were investigated in the gas phase by the trajectory surface-hopping (TSH) simulations with Tully’s fewest-switches algorithm [47,48,49,50,51]. The initial structures and velocities for the nonadiabatic dynamics simulations were selected randomly from a 5 ps ground state trajectory at specific temperature, and then chosen using the filtering procedure implemented in the MNDO program [30] according to the computed S 0 -S 1 transition probabilities. An empirical decoherence correction (0.1 a.u.) suggested by Granucci et al. [52] was employed. A constant time step of 0.1 fs was chosen to solve the nuclear motion equation, while a 100 times smaller time step was selected for the time-dependent electronic propagation.

3. Results and Discussion

3.1. Equilibrium Structures

Four local minima geometries of DDPY in the ground state were obtained based on the OM2/MRCI, B3LYP/6-31G(d), CAM-B3LYP/6-31G(d), and B3LYP-D3/6-31+G(d) level of theories. According to the conformation and helicity, these four equilibrium structures are called EP, EM, ZP, and ZM, respectively. The approach of helicity definition proposed by Karnik et al. [53] was adopted. Geometry of the most stable isomer ZP-DDPY is presented in Figure 1a, while geometries of the other three isomers are shown in Figure S1 (see ESI†). The corresponding geometrical parameters of the four isomers are listed in Table S1 (see ESI†). As can be seen, the optimized geometries obtained from the different theoretical methods above are consistent with each other.
The transition state of EM-EP-TS (between EM and EP isomers) and ZM-ZP-TS (between ZM and ZP isomers) in the ground state were also optimized at the OM2/MRCI, B3LYP/6-31G(d), CAM-B3LYP/6-31G(d), and B3LYP-D3/6-31+G(d) levels. The optimized transition state geometries obtained with the OM2/MRCI method are presented in Figure S2 (see ESI†), while the corresponding geometrical parameters are summarized in Table S2 (see ESI†). As can be seen, these geometrical parameters obtained from the different methods above are in good agreement with each other. According to the obtained transition state, the energy barriers from EM to EP isomers and from ZM to ZP isomers were calculated based on the OM2/MRCI, B3LYP/6-31G(d), CAM-B3LYP/6-31G(d), and B3LYP-D3/6-31+G(d) methods, as shown in Table 1. As we can see, both energy barriers from EM to EP and from ZM to ZP isomers are very low regardless of the method, e.g., just about 0.31 kcal/mol from EM to EP and 0.47 kcal/mol from ZM to ZP obtained from the B3LYP-D3/6-31+G(d) calculations.
The schematic diagram of S 0 and S 1 potential energy profiles along the reaction coordinate is shown in Figure 1b. As we can see, after the S 0 →S 1 optical excitation of EP or ZP isomer, molecular motor DDPY rotates around the central C=C double bond in a counterclockwise direction, and relaxes to the S 0 state through the S 1 /S 0 conical intersections (CIs), then it arrives at the metastable ZM or EM isomer. Motivated by our previous work [25] and Filatov’s work [22], due to the very low energy barriers from ZM to ZP and from EM to EP isomers, as shown in Table 1, we can expect that the molecular motor DDPY may exceed the barriers in timescales of femtoseconds and arrive at more stable ZP or EP isomers without staying at the ZM or EM metastable isomer at room or even lower temperature. Thus, the molecular motor DDPY could complete a full 360 ∘ rotation by only two photoisomerization steps (EP→ZP and ZP→EP) at room or even lower temperature. The schematic diagram of working cycle for photon-only two-stroke LDMRM DDPY is presented in Figure 1c.

3.2. Nonadiabatic Molecular Dynamics Simulations

To verify whether DDPY can work as a photon-only two-stroke LDMRM, nonadiabatic molecular dynamics simulation was systematically carried out. A total of 297 and 325 trajectories starting from the S 1 excited state of EP and ZP isomers at 300 K were firstly studied on the OM2/MRCI level. In addition, 82 of 297 trajectories experienced EP→ZP photoisomerization, in which 61 trajectories finished the photoisomerization through counterclockwise rotation. Meanwhile, 233 of 325 trajectories underwent ZP→EP photoisomerization, in which 168 trajectories finished the photoisomerization through counterclockwise rotation. Thus, the unidirectionalities of successful EP→ZP and ZP→EP photoisomerization processes were estimated to be about 74% and 72%, respectively.
The low unidirectionality may be oweing to the less steric repulsion between the lower half and the upper half of molecular motor DDPY, which results in nearly planar ground state conformations. Low unidirectionality of DDPY at ambient temperature may reduce its application potential as a light-driven molecular rotary motor. How can the unidirectionality of this LDMRM be improved, especially using the physical method? Does temperature influence the nonadiabatic dynamics of molecular motor DDPY? To answer these questions, the ground state dynamics samplings at different temperatures for EP and ZP isomers were carried out. A total of 328, 341, 291, and 297 geometries for EP isomer and 334, 310, 304, and 325 geometries for ZP isomer at 50 K, 100 K, 200 K and 300 K, respectively, were randomly selected and chosen by a filtering procedure according to the computed S 0 -S 1 transition probabilities. The distributions of C4-C2-C1-C23 and C2-C1-C23-N24 dihedral angles of all sampled geometries at different temperatures are illustrated in Figure 2. As we can see, distributions of the sampled geometries become closer to the stable EP and ZP geometries as the temperature decreases. We conjecture that the unidirectionality of this LDMRM may be improved by reducing the temperature. On the basis of trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EP nonadiabatic photoisomerization dynamics of DDPY at different temperatures were systematically studied in the following.

3.2.1. The Nonadiabatic Dynamics of EP→ ZP Photoisomerization

A total of 328, 341, 291, and 297 trajectories starting from the S 1 excited state of EP-DDPY were performed at the OM2/MRCI level for 1000 fs at 50 K, 100 K, 200 K, and 300 K, respectively. The excited state of S 1 corresponds to the single-electron excitation from the HOMO (bonding π orbital) to the LUMO (antibonding π * orbital), with the excitation wavelength at about 377 nm based on the OM2/MRCI level. All trajectories reached the S 0 ground state within 1000 fs. In addition, 104, 105, 100, and 82 trajectories underwent EP→ZP photoisomerization at 50 K, 100 K, 200 K, and 300 K, respectively, which means the quantum yields of EP→ZP photoisomerization at 50 K, 100 K, 200 K, and 300 K are estimated to be about 32%, 31%, 34%, and 28%, respectively. This indicates that a decrease in temperature has little influence on the quantum yields of EP→ZP photoisomerization process of DDPY.
The average occupation of electronic states S 0 and S 1 varying with simulation time at different temperatures are shown in Figure 3. The S 1 time-dependent fractional occupation at different temperatures can all be fitted by a bi-exponential function, as shown in Figure S3 (see ESI†), which indicates that the decay modes of molecular motor DDPY are nearly not affected by lowering the temperature. Based on the S 1 excited state lifetimes of all 328, 341, 291, and 297 trajectories at 50 K, 100 K, 200 K, and 300 K, average lifetimes of the S 1 excited state of the EP-DDPY on above temperatures are estimated to be about 192 fs, 210 fs, 206 fs and 191 fs, respectively. The results show that lowering the temperature does not have a significant impact on the decay mode and average S 1 lifetime of EP→ZP photoisomerization process of molecular motor DDPY.
Based on all geometries at the S 1 /S 0 hopping events, four optimized S 1 /S 0 conical intersections (CIs) were obtained at the OM2/MRCI level, as shown in Figure 4, while the corresponding geometrical parameters are summarized in Table S3 (see SI†). According to the characteristic dihedral angle C4-C2-C1-C23 (108.3 ∘ , 55.7 ∘ , −56.5 ∘ and −109.3 ∘ , the atomic labels can be seen in Figure 1a), and the four CIs are called ECI(1), ECI(2), ZCI(1), and ZCI(2), respectively. As can be seen in Table S3, all CIs involve obvious pyramidalization at the C2 atom site. Similar pyramidalization of the carbon atom at the stator-axle linkage was also observed in other molecular rotary motors [13,25,34,54].
It is helpful for us to understand the decay mechanism through distribution of geometrical parameters at hopping events. For the trajectories experienced EP→ZP photoisomerization at 50 K, 100 K, 200 K, and 300 K, the distributions of C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the S 1 → S 0 hopping events are illustrated in Figure 5. For all trajectories at 50 K, 100 K, 200 K, and 300 K, the distributions of C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the S 1 →S 0 hopping events are also illustrated in Figure S5 (see SI†). The corresponding points of the ground state EP-isomer, conical intersections ECI(1), ECI(2), ZCI(1), and ZCI(2) are also presented in Figure 5 and Figure S5.
Taking 300 K as an example, as shown in Figure 5d, most of the trajectories rotate counterclockwise and some trajectories rotate clockwise. The trajectories of counterclockwise rotation experienced EP→ZP photoisomerization were accessed through hops close to the ECI(1). The trajectories of clockwise rotation experienced EP→ZP photoisomerization were accessed through hops close to the ZCI(1) and ZCI(2). Some hops close to the initial EP structure were also observed, as shown in Figure S5d, but all corresponding trajectories returned to the reactant.
As shown in Figure 5, as the temperature decreases, the proportion of trajectories that rotates clockwise became smaller and smaller. For example, at 50 K, all trajectories that experienced EP→ZP photoisomerization underwent counterclockwise rotation. The statistical unidirectionalities of the successful EP→ZP photoisomerization process are 74%, 77%, 95% and 100% at 300K, 200 K, 100 K, and 50 K, respectively. Thus, lowering the temperature can significantly increase the unidirectionality of the EP→ZP photoisomerization process of molecular motor DDPY.
In order to understand the EP→ZP photoisomerization mechanism of DDPY in detail, time-dependent evolutions of central bond length C1-C2, central dihedral angle C4-C2-C1-C23, side dihedral angle C2-C1-C23-N24, and pyramid dihedral angle C2-N18-C4-C1 in five typical trajectories (called trajectories 1–5, respectively) at 50 K are presented in Figure 6 and Figures S7–S10 (see ESI†). The corresponding geometrical parameters of reaction product ZP isomer and S 1 →S 0 hopping time are also shown in the figures.
Take trajectory 1 as an example, as shown in Figure 6, after the excitation from S 0 to S 1 , the central C1=C2 double bond is weakened, increasing from its optimized ground state value of 1.38 Å to about 1.46 Å, varying around 1.44 Å until the nonadiabatic decay at 242 fs, then returning to about 1.38 Å. That is, the excitation from the bonding π orbital of the central C=C bond to the antibonding π * orbital reduces its double bond character obviously. The dihedral angle C4-C2-C1-C23 increased gradually from 9.1 ∘ to about 95.1 ∘ around 242 fs, after the de-excitation, it increased continually to its optimized ground state value of 183.6 ∘ in the ZP structure at about 400 fs. The dihedral angle C2-N18-C4-C1, characterizing the pyramidalization at the C2 atom, increased to 24.5 ∘ when nonadiabatic decay occurred at 242 fs, then decreasing dramatically to 2.3 ∘ at about 310 fs, and varying around 2.3 ∘ until the end of simulation. Both optimized geometries of conical intersection presented in Figure 4 and the time dependence of geometrical parameters shown in Figure 6 verify that, after the S 0 →S 1 excitation, the dynamical process of nonadiabatic decay is followed by twisting about the central C=C double bond and the pyramidalization of the C atom at the stator–axle linkage.
Side dihedral angle C2-C1-C23-N24 is the key geometrical parameter to distinguish ZP and ZM isomers of DDPY, as can be seen in Table S1 (see ESI†). As shown in Figure 6c, the dihedral angle C2-C1-C23-N24 vibrated around 14.9 ∘ (optimized value in EP geometry) until the nonadiabatic decay at 250 fs, then decreasing dramatically to –1.8 ∘ (optimized value in ZM geometry) at about 300 fs, i.e., molecular motor arrived at the ZM geometry. After staying around the ZM geometry for less than 200 fs, the dihedral angle C2-C1-C23-N24 increased continually to 16.6 ∘ (optimized value in ZP geometry) at about 460 fs, and vibrated around this value until the end of simulation. Time dependence of dihedral angle C2-C1-C23-N24 shown in Figure 6c verifies that, after the S 0 →S 1 excitation of EP isomer, molecular motor DDPY arrives at the ZM isomer firstly, then reaching the ZP isomer in a very short time. For example, the EP→ZP photoisomerization process of DDPY can be realized at 50 K or even lower temperatures, which confirms our expectation in the beginning.

3.2.2. The Nonadiabatic Dynamics of ZP→EP Photoisomerization

With the same method as the above EP→ZP photoisomerization nonadiabatic dynamics simulation, nonadiabatic dynamics of ZP→EP photoisomerization was systemically investigated. At 50 K, 100 K, 200 K, and 300 K, molecular dynamics simulations of 334, 310, 304, and 325 trajectories starting from the S 1 excited state (with the excitation wavelength at about 353 nm) were performed at the OM2/MRCI level for 1000 fs, respectively. All trajectories decayed to the ground state before the end of simulation. A total of 288, 240, 232, and 233 trajectories underwent ZP→EP photoisomerization at 50 K, 100 K, 200 K, and 300 K, respectively, which means the quantum yields of ZP→EP photoisomerization at the above corresponding temperatures are estimated to be about 86%, 77%, 76%, and 72%, respectively. This indicates that the effect of temperature on the quantum yields of ZP→EP photoisomerization process of molecular motor DDPY is not significant.
The average occupation of electronic states S 0 and S 1 varying with simulation time at different temperatures are shown in Figure 7. As we can see, the S 1 population decay at different temperatures are obviously not exponential. Taking a numerical derivative on the occupation of S 0 state over time at different temperatures, as shown in Figure S4 (see SI†), the decay mode of the S 1 excited state are all found to be periodic. Taking 300 K as an example, as is shown in Figure 7d and Figure S4d, four major hopping event maxima arose around 50 fs, 270 fs, 490 fs and 690 fs, respectively. This indicates that the motion of the molecular motor on the S 1 excited state towards the conical intersection is regulated by a periodic structural change. The periodic intervals of hopping event maxima in Figure S4d are roughly in 200–220 fs range, close to a ground state normal mode of ZP-DDPY (148 cm − 1 , the fifth normal mode, corresponding vibrational duration is 225 fs) involving a swing of phenmethyl ring around the central C=C double bond. Similar periodic decay modes have also been observed in Z-E photoisomerization of some azobenzene-based molecules [40,42,43]. Based on the S 1 excited state lifetimes of all 334, 310, 303, and 324 trajectories at 50 K, 100 K, 200 K, and 300 K in our calculations, average lifetime of the S 1 excited state of the ZP-DDPY is estimated to be about 316 fs, 310 fs, 300 fs, and 322 fs, respectively. The results show that lowering the temperature does not have a significant impact on the decay mode and average S 1 lifetime of ZP→EP photoisomerization dynamics process of molecular motor DDPY.
Based on all geometries at the S 1 /S 0 hopping events, four optimized S 1 /S 0 conical intersections were obtained at the OM2/MRCI level, which are the same as those obtained in the above EP→ZP photoisomerization process. For the trajectories experienced ZP→EP photoisomerization at 50 K, 100 K, 200 K, and 300 K, the distributions of C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the S 1 →S 0 hopping events are illustrated in Figure 8. The distributions of C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the S 1 →S 0 hopping events for all trajectories at 50 K, 100 K, 200 K, and 300 K are presented in Figure S6 (see SI†). The corresponding points of ground state ZP-isomer, conical intersections ECI(1), ECI(2), ZCI(1), and ZCI(2) are also displayed in Figure 8 and Figure S6.
Taking 300 K as an example, as shown in Figure 8d, most of the trajectories rotated counterclockwise and some of the trajectories rotated clockwise. The trajectories of counterclockwise rotation experienced ZP→EP photoisomerization were accessed through hops close to the ZCI(1) and ZCI(2). The trajectories of clockwise rotation experienced ZP→EP photoisomerization were accessed through hops close to the ECI(1) and ECI(2). Although a few hops close to the initial ZP structure were also observed, as shown in Figure S6d, all the corresponding trajectories returned to the reactant ZP isomer and did not experience the ZP→EP photoisomerization. As shown in Figure 8, as the temperature decreases, the proportion of clockwise rotation trajectories became smaller and smaller. The statistical unidirectionalities of the trajectories which experienced ZP→EP photoisomerization are 72%, 76%, 88% and 94% at 300 K, 200 K, 100 K, and 50 K, respectively. The same as ZP→EP photoisomerization of molecular motor DDPY, the unidirectionality of ZP→EP photoisomerization process can also be significantly improved through lowering the temperature.
In order to understand the ZP→EP photoisomerization mechanism of DDPY in detail, time-dependent evolutions of central bond length C1=C2, central dihedral angle C4-C2-C1-C23, side dihedral angle C2-C1-C23-N24, and pyramid dihedral angle C2-N18-C4-C1 in five typical trajectories (named as trajectory 1–5, respectively) at 50 K are presented in Figure 9 and Figures S11–S14 (see ESI†). The corresponding geometrical parameters of reaction product EP isomer and S 1 →S 0 hopping time are also shown in the figures.
Take trajectory 1 as an example, as shown in Figure 9, after the excitation from S 0 to S 1 , the central C1=C2 double bond is weakened, increasing from its optimized ground state value of 1.37 Å to about 1.48 Å, varying around 1.44 Å until the nonadiabatic decay at 317 fs, then returning to about 1.37 Å. The dihedral angle C4-C2-C1-C23 increased gradually from −179.1 ∘ to about −46.7 ∘ around 317 fs; after the de-excitation, it increased continually to its optimized ground state value of 5.8 ∘ in the EP structure at about 450 fs. The dihedral angle C2-N18-C4-C1, characterizing the pyramidalization at the C2 atom, increased to 29.7 ∘ when nonadiabatic decay occurred at 317 fs, then decreasing dramatically to 2.4 ∘ around 350 fs, and varying around 2.4 ∘ until the end of simulation. The time dependence of geometrical parameters shown in Figure 9, together with optimized geometries of conical intersections presented in Figure 4, verify that the dynamical process of nonadiabatic decay is followed by twisting about the central C=C double bond and the pyramidalization of the C atom at the stator–axle linkage.
Side dihedral angle C2-C1-C23-N24 is the key geometrical parameter to distinguish EP and EM isomers of DDPY, as can be seen in Table S1 (see ESI†). As shown in Figure 9, the dihedral angle C2-C1-C23-N24 decreased dramatically from 16.6 ∘ (optimized value in ZP geometry) to −35 ∘ within 75 fs and then increased gradually to 20 ∘ at 317 fs. After the nonadiabatic decay at about 317 fs, the dihedral angle C2-C1-C23-N24 decreased dramatically to −6.1 ∘ (optimized value in EM geometry) at about 340 fs, i.e., molecular motor arrived at the EM geometry. After staying around the EM geometry for about 200 fs, the dihedral angle C2-C1-C23-N24 increased to 14.9 ∘ (optimized value in EP geometry) at about 560 fs and then vibrated around this value until the end of simulation. It suggests that the ZP→EP photoisomerization process of DDPY can also be realized at 50 K or even lower temperatures, which also confirms our expectation in the beginning.
Before the conclusions, average lifetime of the S 1 excited state, quantum yield, and unidirectionality at different temperatures for the EP→ZP and ZP→EP photoisomerization of molecular motor DDPY are summarized in Table 2 for comparison. As we can see, as the temperature decreases, average lifetimes of the S 1 excited state and quantum yield of both EP→ZP and ZP→EP photoisomerization are almost unaffected, while the unidirectionalities are significantly increased.

4. Conclusions

Based on electronic structure calculation at the B3LYP/6-31G(d), CAM-B3LYP/6-31G(d), B3LYP-D3/6-31+G(d), and OM2/MRCI level, together with nonadiabatic molecular dynamics simulation at the OM2/MRCI level, a two-stroke light-driven molecular rotary motor, 2-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDPY) is proposed, which is capable of completing a unidirectional 360 ∘ rotation by only two photoisomerization (EP→ZP and ZP→EP) steps. The nonadiabatic dynamics of EP→ZP and ZP→EP photoisomerization of DDPY are investigated by trajectory surface-hopping molecular dynamics at the OM2/MRCI level. Both photoisomerization processes are on an ultrafast timescale (ca. 200–300 fs). The decay mode of EP→ZP photoisomerization is approximately bi-exponential, while the decay mode of ZP→EP photoisomerization was found to be periodic. For EP and ZP isomer of DDPY, after the S 0 →S 1 excitation, the dynamical processes of nonadiabatic decay are both followed by twisting about the central C=C double bond and the pyramidalization of the C atom at the stator–axle linkage.
The effect of temperature on the nonadiabatic dynamics of EP→ZP and ZP→EP photoisomerization of DDPY has been systematically investigated. Based on a large number of trajectories starting from the S 1 excited state of EP and ZP isomer at 50 K, 100 K, 200 K, and 300 K, we found that average lifetimes of the S 1 excited state and quantum yields for both EP→ZP and ZP→EP photoisomerization are almost temperature-independent, while the corresponding unidirectionality of rotation is significantly increased as the temperature decreases. Our present computational results not only proposed a new class of two-stroke photon-only light-driven molecular rotary motor, but also supplied a physical way to increase the unidirectionality of molecular motor, which may stimulate further research for the development of more efficient light-driven molecular rotary motors.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23179694/s1.

Author Contributions

Conceptualization, C.J. and J.M.; methodology, C.J. and J.M.; software, Z.L.; calculation, J.M. and D.Z.; formal analysis, J.M. and D.Z.; resources, C.J.; writing—original draft preparation, J.M.; writing—review and editing, C.J., Z.L. and F.L.; project administration, C.J.; funding acquisition, C.J., Z.L. and F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Natural Science Foundation of China (Grant Nos. 12074307 and 12075177) and the Natural Science Foundation of Shannxi Province (Grant No. 2019JM-410).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was supported by the HPC Platform in Xi’an Jiaotong University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Erbas-Cakmak, S.; Leigh, D.A.; Mcternan, C.T.; Nussbaumer, A.L. Artificial Molecular Machines. Chem. Rev. 2015, 115, 10081–10206. [Google Scholar] [PubMed]
  2. Groppi, J.; Baroncini, M.; Venturi, M.; Silvi, S.; Credi, A. Design of photo-activated molecular machines: Highlights from the past ten years. Chem. Commun. 2019, 55, 12595–12602. [Google Scholar]
  3. Kassem, S.; Leeuwen, T.V.; Lubbe, A.S.; Wilson, M.R.; Feringa, B.L.; Leigh, D.A. Artificial Molecular Motors. Chem. Soc. Rev. 2017, 46, 2592–2621. [Google Scholar] [PubMed]
  4. Kottas, G.S.; Clarke, L.I.; Horinek, D.; Michl, J. Artificial Molecular Rotors. Chem. Rev. 2005, 105, 1281–1376. [Google Scholar]
  5. Baroncini, M.; Silvi, S.; Credi, A. Photo- and Redox-Driven Artificial Molecular Motors. Chem. Rev. 2019, 120, 200–268. [Google Scholar]
  6. Pooler, D.R.S.; Lubbe, A.S.; Crespi, S.; Feringa, B.L. Designing light-driven rotary molecular motors. Chem. Sci. 2021, 12, 14964–14986. [Google Scholar]
  7. Kelly, T.R.; Cai, X.; Damkaci, F.; Panicker, S.B.; Tu, B.; Bushell, S.M.; Cornella, I.; Piggott, M.J.; Salives, R.; Cavero, M.; et al. Progress toward a Rationally Designed, Chemically Powered Rotary Molecular Motor. J. Am. Chem. Soc. 2007, 129, 376–386. [Google Scholar]
  8. Kudernac, T.; Ruangsupapichat, N.; Parschau, M.; Macia, B.; Katsonis, N.; Harutyunyan, S.R.; Ernst, K.-H.; Feringa, B.L. Electrically Driven Directional Motion of a Four-Wheeled Molecule on a Metal Surface. Nature 2011, 479, 208–211. [Google Scholar]
  9. Feringa, B.L. The Art of Building Small: From Molecular Switches to Molecular Motors. J. Org. Chem. 2007, 72, 6635–6652. [Google Scholar]
  10. Koumura, N.; Zijlstra, R.W.J.; Van Delden, R.A.; Harada, N.; Feringa, B.L. Light-driven Monodirectional Molecular Rotor. Nature 1999, 401, 152–155. [Google Scholar]
  11. Filatov, M.; Paolino, M.; Min, S.K.; Kim, K.S. Fulgides as Light-Driven Molecular Rotary Motors: Computational Design of a Prototype Compound. J. Phys. Chem. Lett. 2018, 9, 4995–5001. [Google Scholar] [CrossRef] [PubMed]
  12. Roke, D.; Sen, M.; Danowski, W.; Wezenberg, S.J.; Feringa, B.L. Visible-Light-Driven Tunable Molecular Motors Based on Oxindole. J. Am. Chem. Soc. 2019, 141, 7622–7627. [Google Scholar] [CrossRef] [PubMed]
  13. Pooler, D.R.S.; Pierron, R.; Crespi, S.; Costil, R.; Pfeifer, L.; Leonard, J.; Olivucci, M.; Feringa, B.L. Effect of charge-transfer enhancement on the efficiency and rotary mechanism of an oxindole-based molecular motor. Chem. Sci. 2021, 12, 7486–7497. [Google Scholar]
  14. Wilcken, R.; Schildhauer, M.; Rott, F.; Huber, L.A.; Guentner, M.; Thumser, S.; Hoffmann, K.; Oesterling, S.; De Vivie-Riedle, R.; Riedle, E.; et al. Complete Mechanism of Hemithioindigo Motor Rotation. J. Am. Chem. Soc. 2018, 140, 5311–5318. [Google Scholar] [CrossRef] [PubMed]
  15. Gerwien, A.; Mayer, P.; Dube, H. Photon-Only Molecular Motor with Reverse Temperature-Dependent Efficiency. J. Am. Chem. Soc. 2018, 140, 16442–16445. [Google Scholar] [CrossRef]
  16. Greb, L.; Lehn, J.M. Light-Driven Molecular Motors: Imines as Four-Step or Two-Step Unidirectional Rotors. J. Am. Chem. Soc. 2014, 136, 13114–13117. [Google Scholar] [CrossRef]
  17. Greb, L.; Eichhöfer, A.; Lehn, J.M. Synthetic Molecular Motors: Thermal N Inversion and Directional Photoinduced C=N Bond Rotation of Camphorquinone Imines. Angew. Chem. Int. Ed. 2015, 54, 14345–14348. [Google Scholar] [CrossRef]
  18. Li, Y.; Wang, W.; Liu, F. Exploring the Mechanism of a Chiral N-Alkyl Imine-Based Light-Driven Molecular Rotary Motor at MS-CASPT2//CASSCF and MS-CASPT2//(TD) DFT Levels. Chem. Eur. J. 2019, 25, 4194–4201. [Google Scholar] [CrossRef]
  19. Koumura, N.; Geertsema, E.M.; Van Gelder, M.B.; Meetsma, A.; Feringa, B.L. Second Generation Light-Driven Molecular Motors. Unidirectional Rotation Controlled by a Single Stereogenic Center with Near-Perfect Photoequilibria and Acceleration of the Speed of Rotation by Structural Modification. J. Am. Chem. Soc. 2002, 124, 5037–5051. [Google Scholar] [CrossRef]
  20. Pollard, M.M.; Meetsma, A.; Feringa, B.L. A redesign of light-driven rotary molecular motors. Org. Biomol. Chem. 2008, 6, 507–512. [Google Scholar] [CrossRef]
  21. Garcia-Iriepa, C.; Marazzi, M.; Zapata, F.; Valentini, A.; Sampedro, D.; Frutos, L.M. Chiral Hydrogen Bond Environment Providing Unidirectional Rotation in Photoactive Molecular Motors. J. Phys. Chem. Lett. 2013, 4, 1389–1396. [Google Scholar] [CrossRef] [PubMed]
  22. Filatov, M.; Paolino, M.; Min, S.K.; Choi, C.H. Design and photoisomerization dynamics of a new family of synthetic 2-stroke light driven molecular rotary motors. Chem. Commun. 2019, 55, 5247–5250. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, J.; Durbeej, B. Molecular motors with high quantum efficiency and visible-light responsiveness: Meeting two challenges in one design. Comput. Theor. Chem. 2019, 1148, 27–32. [Google Scholar] [CrossRef]
  24. Majumdar, A.; Jansen, T.L.C. Quantum-Classical Simulation of Molecular Motors Driven Only by Light. J. Phys. Chem. Lett. 2021, 12, 5512–5518. [Google Scholar] [CrossRef] [PubMed]
  25. Ma, J.; Yang, S.; Zhao, D.; Jiang, C.; Lan, Z.; Li, F. Design and Nonadiabatic Photoisomerization Dynamics Study of a Three-Stroke Light-Driven Molecular Rotary Motor. Int. J. Mol. Sci. 2022, 23, 3908. [Google Scholar] [CrossRef]
  26. Boursalian, G.; Nijboer, E.; Dorel, R.; Pfeifer, L.; Markovitch, O.; Blokhuis, A.; Feringa, B.L. All-Photochemical Rotation of Molecular Motors with a Phosphorus Stereoelement. J. Am. Chem. Soc. 2020, 142, 16868–16876. [Google Scholar] [CrossRef] [PubMed]
  27. Paolino, M.; Gueye, M.; Pieri, E.; Manathunga, M.; Fusi, S.; Cappelli, A.; Latterini, L.; Pannacci, D.; Filatov, M.; Léonard, J.; et al. Design, Synthesis, and Dynamics of a Green Fluorescent Protein Fluorophore Mimic with an Ultrafast Switching Function. J. Am. Chem. Soc. 2016, 138, 9807–9825. [Google Scholar] [CrossRef]
  28. The unidirectionality is defined as a fraction of the trajectories propagating in a specific direction, e.g., the counterclockwise direction, while the monodirectionality is defined as preference for one cycle against all others.
  29. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  30. Thiel, W. MNDO Program, Version 6.1; Max-Planck-Institut fur Kohlenforschung: Mulheim, Germany, 2007. [Google Scholar]
  31. Weber, W.; Thiel, W. Orthogonalization Corrections for Semiempirical Methods. Theor. Chem. Acc. 2000, 103, 495–506. [Google Scholar] [CrossRef]
  32. Otte, N.; Scholten, M.; Thiel, W. Looking at Self-Consistent-Charge Density Functional Tight Binding from a Semiempirical Perspective. J. Phys. Chem. A 2007, 111, 5751–5755. [Google Scholar] [CrossRef]
  33. Zhuang, X.; Wang, J.; Lan, Z. Photoinduced Nonadiabatic Decay and Dissociation Dynamics of Dimethylnitramine. J. Phys. Chem. A 2013, 117, 4785–4793. [Google Scholar] [CrossRef]
  34. Pang, X.; Cui, X.; Hu, D.; Jiang, C.; Zhao, D.; Lan, Z.; Li, F. “Watching” the Dark State in Ultrafast Nonadiabatic Photoisomerization Process of a Light-Driven Molecular Rotary Motor. J. Phys. Chem. A. 2017, 121, 1240–1249. [Google Scholar] [CrossRef] [PubMed]
  35. Nikiforov, A.; Gamez, J.A.; Thiel, W.; Filatov, M. Computational Design of a Family of Light-Driven Rotary Molecular Motors with Improved Quantum Efficiency. J. Phys. Chem. Lett. 2016, 7, 105–110. [Google Scholar] [CrossRef] [PubMed]
  36. Che, M.; Gao, Y.J.; Zhang, Y.; Xia, S.H.; Cui, G. Electronic structure calculations and nonadiabatic dynamics simulations of excited-state relaxation of Pigment Yellow 101. Phys. Chem. Chem. Phys. 2018, 20, 6524–6532. [Google Scholar] [CrossRef]
  37. Zhang, Y.H.; Sun, X.W.; Zhang, T.S.; Liu, X.Y.; Cui, G. Nonadiabatic Dynamics Simulations on Early-Time Photochemistry of Spirobenzopyran. J. Phys. Chem. A 2020, 124, 2547–2559. [Google Scholar] [CrossRef] [PubMed]
  38. Gao, Y.J.; Chang, X.P.; Liu, X.Y.; Li, Q.S.; Cui, G. Excited-State Decay Paths in Tetraphenylethene Derivatives. J. Phys. Chem. A 2017, 121, 2572–2579. [Google Scholar] [CrossRef]
  39. Xia, S.H.; Che, M.; Liu, Y.; Zhang, Y.; Cui, G. Photochemical mechanism of 1,5-benzodiazepin-2-one: Electronic structure calculations and nonadiabatic surface-hopping dynamics simulations. Phys. Chem. Chem. Phys. 2019, 21, 10086–10094. [Google Scholar] [CrossRef]
  40. Weingart, O.; Lan, Z.; Koslowski, A.; Thiel, W. Chiral Pathways and Periodic Decay in cis-Azobenzene Photodynamics. J. Phys. Chem. Lett. 2011, 2, 1506–1509. [Google Scholar] [CrossRef]
  41. Lan, Z.; Lu, Y.; Weingart, O.; Thiel, W. Nonadiabatic Decay Dynamics of a Benzylidene Malononitrile. J. Phys. Chem. A 2012, 116, 1510–1518. [Google Scholar] [CrossRef]
  42. Wang, Y.T.; Liu, X.Y.; Cui, G.; Fang, W.H.; Thiel, W. Photoisomerization of Arylazopyrazole Photoswitches: Stereospecific Excited-State Relaxation. Angew. Chem. Int. Ed. 2016, 55, 14009–14013. [Google Scholar] [CrossRef]
  43. Pang, X.; Jiang, C.; Qi, Y.; Yuan, L.; Hu, D.; Zhang, X.; Zhao, D.; Wang, D.; Lan, Z.; Li, F. Ultrafast unidirectional chiral rotation in the Z-E photoisomerization of two azoheteroarene photoswitches. Phys. Chem. Chem. Phys. 2018, 20, 25910–25917. [Google Scholar] [CrossRef]
  44. Cui, G.; Lan, Z.; Thiel, W. Intramolecular Hydrogen Bonding Plays a Crucial Role in the Photophysics and Photochemistry of the GFP Chromophore. J. Am. Chem. Soc. 2012, 134, 1662–1672. [Google Scholar] [CrossRef] [PubMed]
  45. Keal, T.W.; Wanki, M.; Thiel, W. Assessment of Semiempirical Methods for the Photoisomerisation of a Protonated Schiff Base. Theor. Chem. Acc. 2009, 123, 145–156. [Google Scholar] [CrossRef]
  46. Keal, T.W.; Koslowski, A.; Thiel, W. Comparison of Algorithms for Conical Intersection Optimization Using Semiempirical Methods. Theor. Chem. Acc. 2007, 118, 837–844. [Google Scholar] [CrossRef]
  47. Tully, J.C. Molecular Dynamics with Electronic Transitions. J. Chem. Phys. 1990, 93, 1061–1071. [Google Scholar] [CrossRef]
  48. Schiffer, S.H.; Tully, J.C. Proton Transfer in Solution: Molecular Dynamics with Quantum Transitions. J. Chem. Phys. 1994, 101, 4657–4667. [Google Scholar] [CrossRef]
  49. Barbatti, M.; Granucci, G.; Persico, M.; Ruckenbauer, M.; Vazdar, M.; Eckert-Maksic, M.; Lischka, H.J. The On-the-Fly Surface-hopping Program System NEWTON-X: Application to Ab Initio Simulation of the Nonadiabatic Photodynamics of Benchmark Systems. J. Photochem. Photobiol. A 2007, 190, 228–240. [Google Scholar] [CrossRef]
  50. Fabiano, E.; Groenhof, G.; Thiel, W. Approximate Switching Algorithms for Trajectory Surface Hopping. Chem. Phys. 2008, 351, 111–116. [Google Scholar] [CrossRef]
  51. Fabiano, E.; Keal, T.W.; Thiel, W. Implementation of surface hopping molecular dynamics using semiempirical methods. Chem. Phys. 2008, 349, 334–347. [Google Scholar] [CrossRef]
  52. Granucci, G.; Persico, M.; Zoccante, A. Including quantum decoherence in surface hopping. J. Chem. Phys. 2010, 133, 134111. [Google Scholar] [CrossRef]
  53. Karnik, A.V.; Hasan, M. Stereochemistry: A Three-Dimensional Insight; Elements of Chirality and Chiral Stereoisomerism; Elsevier: Amsterdam, The Netherlands, 2021; Chapter 3. [Google Scholar]
  54. Kazaryan, A.; Lan, Z.; Schafer, L.V.; Thiel, W.; Filatov, M. Surface Hopping Excited-State Dynamics Study of the Photoisomerization of a Light-Driven Fluorene Molecular Rotary Motor. J. Chem. Theory Comput. 2011, 7, 2189–2199. [Google Scholar] [CrossRef] [Green Version]
Ijms 23 09694 g001 550
Figure 1. (a) Optimized geometry of ZP isomer of DDPY. Some atoms around the central C=C double bond are labeled; (b) the schematic diagram of ground and the first excited state potential energy profiles along the reaction coordinate of LDMRM DDPY; (c) the schematic diagram of a working cycle of the photon-only two-stroke LDMRM DDPY.
Figure 1. (a) Optimized geometry of ZP isomer of DDPY. Some atoms around the central C=C double bond are labeled; (b) the schematic diagram of ground and the first excited state potential energy profiles along the reaction coordinate of LDMRM DDPY; (c) the schematic diagram of a working cycle of the photon-only two-stroke LDMRM DDPY.
Ijms 23 09694 g001
Ijms 23 09694 g002 550
Figure 2. The distribution of C4-C2-C1-C23 and C2-C1-C23-N24 dihedral angles of all sampled initial geometries in Franck–Condon region for (a) ZP→EP photoisomerization and (b) EP→ZP photoisomerization at 300 K, 200 K, 100 K and 50 K. The corresponding points of ground state ZP and EP isomers are also presented in this figure.
Figure 2. The distribution of C4-C2-C1-C23 and C2-C1-C23-N24 dihedral angles of all sampled initial geometries in Franck–Condon region for (a) ZP→EP photoisomerization and (b) EP→ZP photoisomerization at 300 K, 200 K, 100 K and 50 K. The corresponding points of ground state ZP and EP isomers are also presented in this figure.
Ijms 23 09694 g002
Ijms 23 09694 g003 550
Figure 3. Average occupation of the electronic states S 0 and S 1 as a function of simulation time in EP→ZP photoisomerization of DDPY at (a) 50 K, (b) 100 K, (c) 200 K, and (d) 300 K, respectively.
Figure 3. Average occupation of the electronic states S 0 and S 1 as a function of simulation time in EP→ZP photoisomerization of DDPY at (a) 50 K, (b) 100 K, (c) 200 K, and (d) 300 K, respectively.
Ijms 23 09694 g003
Ijms 23 09694 g004 550
Figure 4. Optimized geometries of four S 1 /S 0 conical intersections (a) ECI(1), (b) ECI(2), (c) ZCI(1), and (d) ZCI(2) in the EP→ZP and ZP→EP photoisomerization processes calculated with the OM2/MRCI method implemented in the MNDO99 program.
Figure 4. Optimized geometries of four S 1 /S 0 conical intersections (a) ECI(1), (b) ECI(2), (c) ZCI(1), and (d) ZCI(2) in the EP→ZP and ZP→EP photoisomerization processes calculated with the OM2/MRCI method implemented in the MNDO99 program.
Ijms 23 09694 g004
Ijms 23 09694 g005 550
Figure 5. Distribution of the C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the hopping events of successful EP→ZP photoisomerization trajectories starting from the EP structure of DDPY at (a) 50 K, (b) 100 K, (c) 200 K, and (d) 300 K, respectively. The corresponding points of ground state EP isomer, S 1 /S 0 conical intersections ECI(1), ECI(2), ZCI(1), and ZCI(2) are also presented in the figure.
Figure 5. Distribution of the C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the hopping events of successful EP→ZP photoisomerization trajectories starting from the EP structure of DDPY at (a) 50 K, (b) 100 K, (c) 200 K, and (d) 300 K, respectively. The corresponding points of ground state EP isomer, S 1 /S 0 conical intersections ECI(1), ECI(2), ZCI(1), and ZCI(2) are also presented in the figure.
Ijms 23 09694 g005
Ijms 23 09694 g006 550
Figure 6. Time dependence of (a) central bond length C1-C2; (b) central dihedral angle C4-C2-C1-C23; (c) side dihedral angle C2-C1-C23-N24; and (d) pyramid dihedral angle C2-N18-C4-C1 in a representative trajectory (called trajectory 1) of EP→ZP photoisomerization process. The S 1 →S 0 hopping time (blue line) and corresponding geometrical parameters of reaction product ZP isomer (red lines) are also shown in the figure.
Figure 6. Time dependence of (a) central bond length C1-C2; (b) central dihedral angle C4-C2-C1-C23; (c) side dihedral angle C2-C1-C23-N24; and (d) pyramid dihedral angle C2-N18-C4-C1 in a representative trajectory (called trajectory 1) of EP→ZP photoisomerization process. The S 1 →S 0 hopping time (blue line) and corresponding geometrical parameters of reaction product ZP isomer (red lines) are also shown in the figure.
Ijms 23 09694 g006
Ijms 23 09694 g007 550
Figure 7. Average occupation of the electronic states S 0 and S 1 as a function of simulation time in ZP→EP photoisomerization process of DDPY at (a) 50 K, (b) 100 K, (c) 200 K, and (d) 300 K, respectively.
Figure 7. Average occupation of the electronic states S 0 and S 1 as a function of simulation time in ZP→EP photoisomerization process of DDPY at (a) 50 K, (b) 100 K, (c) 200 K, and (d) 300 K, respectively.
Ijms 23 09694 g007
Ijms 23 09694 g008 550
Figure 8. Distribution of the C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the hopping events of successful ZP→EP photoisomerization trajectories starting from the ZP structure of DDPY at (a) 50 K; (b) 100 K; (c) 200 K; and (d) 300 K, respectively. The ground state ZP isomer, ECI(1), ECI(2), ZCI(1), and ZCI(2) are also presented in this figure at different temperatures.
Figure 8. Distribution of the C4-C2-C1-C23 and C2-N18-C4-C1 dihedral angles at the hopping events of successful ZP→EP photoisomerization trajectories starting from the ZP structure of DDPY at (a) 50 K; (b) 100 K; (c) 200 K; and (d) 300 K, respectively. The ground state ZP isomer, ECI(1), ECI(2), ZCI(1), and ZCI(2) are also presented in this figure at different temperatures.
Ijms 23 09694 g008
Ijms 23 09694 g009 550
Figure 9. Time dependence of (a) central bond length C1-C2; (b) central dihedral angle C4-C2-C1-C23; (c) side dihedral angle C2-C1-C23-N24; and (d) pyramid dihedral angle C2-N18-C4-C1 in a representative trajectory (named as trajectory 1) of ZP→EP photoisomerization process. The S 1 →S 0 hopping time (blue line) and correponding geometrical parameters of reaction product EP isomer (red lines) are also shown in the figure.
Figure 9. Time dependence of (a) central bond length C1-C2; (b) central dihedral angle C4-C2-C1-C23; (c) side dihedral angle C2-C1-C23-N24; and (d) pyramid dihedral angle C2-N18-C4-C1 in a representative trajectory (named as trajectory 1) of ZP→EP photoisomerization process. The S 1 →S 0 hopping time (blue line) and correponding geometrical parameters of reaction product EP isomer (red lines) are also shown in the figure.
Ijms 23 09694 g009
Table
Table 1. Energy barriers from EM to EP isomers and from ZM to ZP isomers in the ground state, obtained from the OM2/MRCI, B3LYP/6-31G(d), CAM-B3LYP/6-31G(d), and B3LYP-D3/6-31+G(d) methods. The energy unit is kcal/mol.
Table 1. Energy barriers from EM to EP isomers and from ZM to ZP isomers in the ground state, obtained from the OM2/MRCI, B3LYP/6-31G(d), CAM-B3LYP/6-31G(d), and B3LYP-D3/6-31+G(d) methods. The energy unit is kcal/mol.
OM2/MRCIB3LYP/6-31G(d)CAM-B3LYP/6-31G(d)B3LYP-D3/6-31+G(d)
EM→EP0.170.200.050.31
ZM→ZP0.090.600.560.47
Table
Table 2. The average S 1 lifetime, quantum yield, and unidirectionality for EP→ZP and ZP→EP photoisomerization nonadiabatic dynamics simulation at 300 K, 200 K, 100 K, and 50 K, respectively.
Table 2. The average S 1 lifetime, quantum yield, and unidirectionality for EP→ZP and ZP→EP photoisomerization nonadiabatic dynamics simulation at 300 K, 200 K, 100 K, and 50 K, respectively.
TemperatureAverage S 1 LifetimeQuantum YieldUnidirectionality
EP→ZP300 K191 fs28%74%
200 K206 fs34%77%
100 K210 fs31%95%
50 K192 fs32%100%
ZP→EP300 K322 fs72%72%
200 K300 fs76%76%
100 K310 fs77%88%
50 K316 fs86%94%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, J.; Zhao, D.; Jiang, C.; Lan, Z.; Li, F. Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor. Int. J. Mol. Sci. 2022, 23, 9694. https://doi.org/10.3390/ijms23179694

AMA Style

Ma J, Zhao D, Jiang C, Lan Z, Li F. Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor. International Journal of Molecular Sciences. 2022; 23(17):9694. https://doi.org/10.3390/ijms23179694

Chicago/Turabian Style

Ma, Jianzheng, Di Zhao, Chenwei Jiang, Zhenggang Lan, and Fuli Li. 2022. "Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor" International Journal of Molecular Sciences 23, no. 17: 9694. https://doi.org/10.3390/ijms23179694

APA Style

Ma, J., Zhao, D., Jiang, C., Lan, Z., & Li, F. (2022). Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor. International Journal of Molecular Sciences, 23(17), 9694. https://doi.org/10.3390/ijms23179694

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop