On the Nature of Interplay among Major Flexibility Channels in Molecular Rotors

As a part of our interest in the excited-state dynamics of flexible materials, we have undertaken a theoretical investigation to the photo-induced reactions of 2-[4-(dimethylamino)benzylidene]malononitrile (BMN) by a combination of the density functional theory, its extended time-dependent (TD-DFT) single reference, and ab initio molecular dynamic (MD) simulations. )e results showed that double-bond twisting and the neighbor single-bond twisting togetherness in the excited singlet state is the most important nonradiative deactivation channel to the ground state. Doubleand single-bond twisting insert clear intersections among the potential energy surfaces of the singlet states (especially S1/S0) leading to fluorescence quenching. Furthermore, effects of molecular dynamic simulations on molecular properties in the femtosecond to picosecond time domain are studied to validate the results. In agreement with the experimental results, the findings conclude the existence of a flexible geometry-dependent single emission band. Such a studymay give information on how themolecule could be externally modified/fixed to yield a desired effect, i.e., more fluorescence or more nonradiative decay.

e emission properties are potentially environment-dependent, which makes TICT-based fluorophores ideal sensors for medium or chemical species nature 34].Several TICT-based materials have been reported to become fluorescent upon aggregation or confinement in solid polymer films due to inhibition of twisting [30].Furthermore, various recent studies in organic optoelectronics, nonlinear optics, and solar energy conversions utilized the concept of twisting to modulate the electronic-state mixing and coupling on charge transfer states.Most fluorescence molecular motors are based on a twisted intramolecular charge transfer (TICT) mechanism [34] in which the nonradiative decay of an excited state can be altered by the surrounding viscosity and have been developed for their rapid responses and high special resolution [35].
Some models of nonradiative decay were proposed [31][32][33] on the basis of experimentally obtained photophysical data of a series of intramolecular charge-transfer complexes of the varied structure (2-[4-(dimethylamino) benzylidene]malononitriles and coumarin dyes).e models specifically describe the relative contributions of the possible internal relaxation channels to nonradiative energy dissipation of the S 1,ct states via the free rotor mechanism.
e results were supported by simple theoretical predictions based on the semiempirically calculated electron density distributions and the bond-order alternation in the S 1,ct state.e photo-induced nonadiabatic decay dynamics of BMN in the gas phase was investigated by iele et al. using surface hopping simulations [33].
e semiempirical orthogonalization model 2 (OM2) combined with multireference configuration interaction (MRCI) was employed [31][32][33].In the nonadiabatic dynamics of DMN, the twisting around the double bond acts as the driving coordinate toward the lowest-energy S 0 /S 1 conical intersection (CI), which mediates the internal conversion to the electronic ground state.is conical intersection is also characterized by a pronounced pyramidalization at the C10 atom that accompanies the double-bond twisting.
e work of iele et al. provides the motivation to check the theoretical insights and the proposed dynamics scenario for the gas phase BMN by high-level ab initio calculations.
Despite all of the reported experimental and theoretical results and the fact that numerous examples of flexible ICT over a wide range of applications exist, it is reassuring to know that only a relatively small number of design principles are involved.us, the future for flexible ICT is promising, and there is still a room for further new research.In a related system, the interplay between the different modes was recently investigated experimentally [36].Generally, it is concluded that the photophysics of these potentially useful probes and functionally similar malononitriles are still partly comprehended and actively debated.
In this investigation, we report on the DFT theoretically obtained results of donor-acceptor flexible ICT 2-[4-(dimethylamino)benzylidene]malononitrile (1) molecule in which three main internal relaxation channels marked in Figure 1 (C4-C9 single bond (S1), C1-N1 single bond (S2), and C9-C10 double bond (D)) are most probable.Simulations are carried out in BMN ground (S 0 ) and few excited singlet states.e molecular rotor BMN could be appraised as an α, β-enenitrile derivative with an increase in the complexity and the size.e objectives are to simulate, in a simple and sufficient accurate way, internal dynamics including twisting modes, reflecting the role played by most important deactivation channel(s) and to evaluate the influence of three modes on the molecular and spectroscopic properties.In other words, twisting around different bonds should be echoed in molecular and spectroscopic characteristics, which will be shown in this contribution.Moreover, the energetic requirement for twisting about different possible bonds in BMN will be calculated using different models.Furthermore, to validate the results, we carried out ab initio molecular dynamics (MD) [37,38], which involve nuclear dynamics to determine the time evolution of the molecular geometry in combination with electronic structure methods.It is capable of computing electronic ground and excited state potential energy surfaces (PESs).A large number of trajectories are usually produced from the MD simulations.Such simulations directly include all nuclear degrees of freedom, which provides a rather rich picture of the time-evolved microscopic processes at a reasonable computing cost.It is well-known that simulations using multireference alternative methods such as the complete active space self-consistent field (CASSCF) method and the various multireference dynamic correlation calculations are very time-consuming [39].To be concise, however, applying TD-DFT single reference ab initio MD to our system will provide a simplified picture about proximity and intersections that might be occurred among PESs and could help explanation of the nonradiative nature of molecular rotors in a relatively less costly and simplified way.We will discuss the photo-induced dynamics in the nonadiabatic channel of the charge-transfer feature material BMN of the molecular rotor type.e obtained results should not only clarify the general understanding of how photo-induced properties change but also to help evolving useful controllable functional materials photochemically by controlling mode(s) of relaxation.

Experimental and Computational Methods
2.1.Electronic and FTIR Spectra of BMN.Measurements of absorption and fluorescence spectra were carried out in acetonitrile solutions (2.0 × 10 −5 mole•L −1 of BMN).e FTIR was measured in the solid KBr tablet.e material was supplied by Dr. Rafik O. Loutfy, Toronto, Canada.
Single point energy computations on optimized structures obtained at different constraints were performed using the cost-effective B3LYP/def2-SV(P) model for unexpansive computing of thermodynamic parameters using Spartan'16 software [41,42].e built-in graphics user interface (GUI) of Spartan'16 is used for visualizing the results of Spartan computations.
To simulate bond twist, we performed relaxed scans on the S 1 state geometry by applying the CIS/6-31G * method.
ORCA version 4.0.1.2[39] allows the combination of the scan feature with CIS or TD-DFT. is can be used to map out the excited state potential energy surfaces as a function of twisting dihedral angle.
e output of the trajectory run automatically contains the excited state energies in addition to the ground state energy.Our calculation utilizes the range-separated hybrid DFT basis: def2-TZVPD [43,44] and the auxiliary basis: def2/JK.e used functional is BP86 with D3 correction and RI approximation [45].
is method speeds up the computation and reasonably estimates the ICT state.Overlap-fitted RIJCOSX approximation was used as a speed-up option that leads to very large speedups [46,47] at virtually no loss of accuracy [48].
Moreover, Orca enabled us to carry out potential energy surface scans along normal coordinates of IR modes obtained from the S 0 computations, which should validate scanning dihedrals due to its direct relationship to molecular dynamics.
e ground and excited state potential energy surfaces can be mapped as a function of normal coordinates.First, we ran a frequency job using BP86/def2-SV(P) and auxiliary basis def2/J to generate the Hessian matrix.e Hessian matrix was used to construct normal mode displacements [39] (known as motions of the normal mode coordinate).is is followed by using TD-DFT single point scan calculations for a range of the displaced geometries (q).e trajectories were constructed so that corresponding normal coordinates are varied in the chosen displacement range.e DFT-MD simulations were done at the TD-DFT Cam-B3LYP/def2-SVP with RIJCOSX approximation and BP86-D3/def2-SVP levels of theory with RI approximation.Both methods return almost the same results.Moreover, in the trajectory simulation, the temperature is kept constant at 350 K, and the run was submitted for a number of runs aiming at exploring a conformational time domain of 50 fs, 100, 500, and 1000 of 0.5, 1, and 2 fs step, respectively.Because the main emphasis is on the flexible skeleton of the molecule studied, the chosen 2 fs step used for 500 runs is mainly to reduce the computational cost.Orca results are visualized by ChemCraft visualization software.
A Broadberry workstation (40 cores) (UK) and a Mac Pro (12 core) workstation were used.

Spectroscopic Characteristics of BMN.
Experimentally determined electronic spectra in acetonitrile and FTIR vibrational spectrum are depicted in Figures 2 and 3, respectively.e noticed differences between experimental IR data in the solid phase (KBr disk) and theoretically computed ones in the gaseous phase is expected in such a flexible molecule suggesting that the geometry in solid is slightly disagreed with that computed in the gas phase.However, it is wellknown that the vibrational analysis by Gaussian can only give some rough ideas on normal modes.Vibrational modes (mainly stretching modes) involving the bonds responsible for molecular flexibility are identified and used to carry out PES scans along normal coordinates (normal mode displacements), which will be discussed in Section 3.4.Moreover, experimentally obtained fluorescence spectrum in acetonitrile is in excellent agreement with theoretical simulation discussed in Section 3.3.

Relative Energy.
Figure 4 shows that twisting around the S1-single bond (the C4-C9 (p-N, N-dimethyl aryl group)) induces largest thermal instability relative to twisting around S2-bond or D-bond.It is obvious that twisting around the C9-C10 double bond induces smallest molecular thermal instability.Figure 5 is a sketch of relative potential energy diagram of BMN in the planar S 0 and twisted S 1,ct states.At 90 °, twist angle of the double bond C9-C10 is largest and smallest destabilization energy is noticed in S 0 and S 1,ct , respectively.Twisting of the anilino moiety (around C4-C9 single bond) induced largest destabilization in the first excited state S 1,ct relative to the nontwist state.In conclusion, based on relative potential energy values (refer to Figures 4  and 5), the favorable tendency for radiation-less deactivation process decreases as follows: D-bond > S2-bond > S1-bond.In other words, double-bond twist is the key factor in nonradiative deactivation of the first excited state, most probably, due to the noticed decrease in the calculated energy gap between states (S 1 -S 0 ) at 90 °twist.

Fluorescence Properties.
Fluorescence deactivation of BMN occurs through twist involving the dicyanomethylene bond (C9-C10 in Figure 1), upon which the ground-and the excited-state potential energy surfaces (PESs) come into sufficiently close proximity, or intersect each other, to allow efficient deactivation (Section 3.4).Figure 6 shows that twistinduced fluorescence wavelength and intensity changes are most remarkable at 90 °twist.Progressive increase in the twist angle results in fluorescence red shift with quenched emission.Similar behavior is noticed in case of simultaneous twists involving different possible combinations of two bonds (see Figure 7).e effect of simultaneous twisting of D-S1 and D-S2 bonds on the relative energy (relative to the energy of the nontwisting geometry) of the molecule is depicted in Figure 8. e 2D drawing reflect the importance of D-S1 twisting in destabilizing the molecule relative to D-S2, which should be reflected in the simulated fluorescence spectra as indicated in Figure 9.
Fluorescence spectra in acetonitrile due to simultaneous two bonds twist angles variation are depicted in Figures 9(a Red shifted emission bands are characterized by low intensity due to larger progressive nonradiative deactivation with increasing bond twisting angles, which is more pronounced in case of D-S1 relative to D-S2. Figure 6 illustrates the effect of bond twisting at least qualitatively.Furthermore, characters of the low-lying electronic states at equilibrium geometries are given in Table 1.All states have mixed character. We previously reported some experimentally determined fluorescence characteristics of BMN in different solvents [31][32][33], and the data are in agreement with the theoretically obtained results.Furthermore, experimentally measured fluorescence spectrum given in Figure 2  e MO involved in the electronic transitions (Figure 10(a)) identifies the nature of transitions in case of the planar nontwisted molecule.It is obvious that the more the HOMO-LUMO involvement in a transition is, the more its CTnature is, as shown in Figure 10(b).(HOMO-LUMO in case of twisting).

PES in Different Electronic
States.In order to identify, at least qualitatively in a simple way, conical intersections (CIs) [49] between electronic potential energy surfaces and the roles played in molecular relaxation processes, we carried out several scans of the dynamics of the BMN material with the help of the molecular mechanical and ab initio trajectories methods.It is well-known that simulations using alternative methods such as the complete active space selfconsistent field (CASSCF) method and the various multireference dynamic correlation calculations, which are proper methods, are very time-consuming.One should also be aware that, with multireference methods, it is very easy to let a large computer run for a long time and still not to produce a meaningful result [37,38,41,42].us, we tried to obtain at least a qualitative picture by carrying out several scans of the dynamics of the BMN material with the help of the molecular mechanical and ab initio trajectories methods.
Potential energy surface scans of different ground and excited states were performed by two techniques: (1) scanning dihedrals and (2) scanning along normal coordinates of IR modes, which should validate scanning dihedrals due to its direct relationship to molecular dynamics.
e ground and excited state potential energy surfaces can be also mapped as a function of normal coordinates.e results PES obtained by scanning dihedrals revealed the presence of apparent intersection points between the singlet state and the ground state, corresponding to C9-C10 double-bond twisting (Figure 11).Apparent intersection in case of single-bond twisting is also noticed (Figure 11(a)).Passage through these apparent intersections on the PES leads to nonradiative energy dissipation due to bond twisting.Recalling, it is well-understood that our single reference computational methodology cannot describe the CI topology.
is is why we clearly see that the S 1 state becomes lower than the S 0 state near 90 °. is feature was discussed before [50][51][52][53].As pointed out before [50], one should not expect good agreement in cases where charge transfer is important, since ionic states are often strongly stabilized by dynamical electron correlation.
However, our simplified approach clearly shows the existence of points of intersections in certain geometries, which    could be helpful in describing photochemical-induced dynamics.More accurate and quantitative examinations of CI would require a CASSCF computation with large active spaces [50][51][52][53].It should be pointed out here that applying the CASSCF method would result in accurate description of the double cone topology of the CI.It is clear that single reference TD-DFT would provide a much more economical alternative, with a possible application to large molecules.

Journal of Chemistry
It could be safely concluded that simulation of molecular dynamics of the complex BMN molecule in case of S1-bond and D-bond twisting reveals that double-bond twisting and the neighbor single-bond twisting togetherness in the excited singlet state is the most important nonradiative deactivation channel to the ground state.
e noticed PES intersections between different excited and ground states strongly suggest that only one radiative transition (one fluorescence band) is most probable.is is consistent with the theoretically and experimentally obtained fluorescence spectra in acetonitrile.
Furthermore, PESs obtained from scanning along normal coordinates of some relevant IR modes (Figure 11) were obtained and analyzed.Although no clear PES intersection is noticed between S 1 and S 0 states, different state trajectories become closer to each other due to vibrational modes mostly related to dynamical motions involving the flexible coordinates, i.e., the S1-, S2-, and D-bonds. is should result in radiation-less relaxation of the excited BMN molecule.

Ab Initio Molecular Dynamics.
In order to further identify intersection points between electronic potential energy surfaces and the roles played in molecular relaxation processes, ab initio molecular dynamics are investigated.MD trajectories were computed at the geometry of singlet/singlet crossing points depicted in Figure 11, i.e., starting from 90 °and 176 °, as shown in Figure 11, for single-and double-bond twisting, respectively.en, the physical quantities (energies, bond lengths, dihedrals, and angles) that are needed in the time evolutions, have been calculated at the CIS/6-31G * level of theory and visualized using the Spartan'16 software packages.Our data graphically visualized in Figures 12 and 13 confirm that the   It is clear from Figure 13(c) that as the time evolves, the molecule becomes more relaxed (relative energy becomes more negative) due to radiation-less deactivation.
Ultrafast internal conversion (IC) (femto-to-picosecond time domain) from the first excited electronic state (S 1 ) to the ground electronic state (S 0 ) via a conical intersection (CI) can play an essential role in the initial steps of the decomposition of energetic materials.Such nonradiative relaxations following electronic excitation quench emission dissipate the excitation energy in the vibrational modes of the ground electronic state.e structure then exhibits irregular periodic energy stabilization-destabilization amounts to 100-25 kJ/mol during the first 50 fs (with a major geometry change due to involvement of bonds twisting, in-plane and out-of-plane bending, as well as H-atoms stretching and bending (rocking and wagging) beside ring deformation and ring C-C stretching).en, in about 200-250 fs further, the molecule enjoys stabilization.In about 800 fs, further stabilization mounts to −40 kJ/mol is achieved.Pyramidalization of C10 (Figure 13(b) is obviously attained in about 12 fs confirming the twisting around the C10-C9 double bond.Figure 13(b) indicates also that double bond and S1-single-bond twisting tends to change almost periodically out of phase with respect to each other, whereas C1-N1 single-bond twisting starts to increase after 30 fs. is suggests that, points of intersections (176 °and 90 °for D and S1-bonds, respectively), both D-bond and S1-bond behave similarly as time evolves and most probably simultaneously play major role in the radiation-less deactivation process of the excited states.

Conclusion
We discussed the photo-induced dynamics in the nonadiabatic channel of the charge-transfer feature material of the molecular rotor type.And we tried to clarify the relaxation channels on the excited state of BMN by using electronic structure calculations on-the-fly dynamics.We report on the DFT, TD-DFT, and ab initio MD theoretically obtained results of stimuli-responsive molecular and spectroscopic properties induced by controlled twisting pathways of the charge transfer molecule, 2-[4-(dimethylamino)benzylidene] malononitrile (BMN) as a model compound of donoracceptor material.Simulations are carried out in BMN ground (S 0 ) and first few excited states in vacuum and in acetonitrile.It is verified that following dynamical displacements of vibrational vectors and intramolecular twisting, the excited states of this flexible molecule return to the ground state either through red-shifted emission or by nonradiative relaxation via conical intersections.e present simulation study shed lights on simultaneous all possible two bond twists in an ICT molecule BMN.Induced changes of molecular properties such as bond length, relative energy, and electronic spectroscopic characteristics (UV-Vis and fluorescence) by simultaneous two bonds twist should deepen our view to the TICT states in flexible molecules of potential applications by identifying, and probably controlling, the roles played by specific bonds twisting.e results are validated by carrying out several scans of the dynamics of the BMN material with the help of the molecular mechanical and ab initio trajectories methods.Simulation of molecular dynamics of the complex BMN molecule showed that double-bond twisting and the neighbor single-bond twisting togetherness in the excited singlet state is the most important nonradiative deactivation channel to the ground state.
Points of intersections between electronic potential energy surfaces are identified by a single reference TDDFT computation and play a crucial role in arguing molecular relaxation pathways.Our data confirms that the propensity for a radiation-less transition is large in the vicinity of nuclear configurations where potential energy surfaces intersect.More sophisticated multireference methods brought about the same conclusions.
Hopefully, this study may provide the foundation for more future investigations on the fluorescence quenching mechanism in D-A materials, which may be reflected in the photochemical molecular design of new efficient materials for a wide range of profitable and new applications.Additionally, it will be interesting to check the theoretical insights and the proposed dynamics scenario for BMN from the current single reference molecular dynamics simulations by high-level ab initio calculations.

Figure 1 :
Figure 1: Optimized planar molecular structure of BMN showing the three main rotational pathways.e molecule is also planar in the first excited state.
validates data of Figure 10(b) at twisting coordinates 126 and 90 °represented by Step 24.

Figure 9 :
Figure 9: Simulated fluorescence spectra in acetonitrile due to simultaneous two bonds twist angles variation: (a) and (b) represent fluorescence due to relaxed scan steps 1-25 for D-S1 (the insets show the dihedral angles of scan grid points); (c) represents the case of D-S2 scans.Computed using TD-DFT (cam-B3LYP/6-31G * model).All spectra are plotted by GaussView 6.

FIGURE 10 :
FIGURE 10: Gaussian'16 results depicts (a) the MO surfaces involved in the electronic transitions from the 1 st and 2 nd singlet excited states to the ground state and shows clearly (b) the CT nature reflected in HOMO-LUMO involved in electronic transition in a twisted state.Both MO's shown in (a) and (b) are obtained by Gaussian'16 package.We used different color codes to illustrate the difference between planar and twisted structures.e isovalue for all surfaces is 0.02.All surfaces are plotted by GaussView 6.(a) Planar configuration.(b) Twisted configuration: scan step 18 shown in the inset Figure 9(b).

Figure 12 :
Figure 12: PESs scanned along some IR normal coordinate modes in the ground state (a) shows PESs of ground and the first two excited singlet states scanned along the displacement coordinate (q) due to different IR modes (from the left to right): mode 63 (at 1597 cm −1 mainly due to D-, S1-, and S2-bond vibration), mode 64 (at 1633 cm −1 due to D-and S2-bonds), and mode 65 (at 1682.8 cm −1 mainly due to D-bond vibration), respectively.(b) An example of most important mode 65 visualizing the displacement vibrational vectors (shown as arrows) involved in the dynamical motions due to the molecular flexibility.

Figure 13 Figure 13 :
Figure 13 summaries the results of the trajectory dynamic simulation in 50 and 1000 fs time domains.Inspection of Figures 13(a) and 13(b) reveals rapid change in dynamics parameters/energies which are noticed even during the first 50 fs time period.Regarding changes in bond length by time evolution, Figure 13(a) indicates that D-bond, S1-bond, and S2-single bond lengths remarkably changed during first 30 fs.S1 single bond length shows small changes while both other two bond lengths as the time evolves change periodically at the points of intersections.Time evolution graphs (Figures13(a) and 13(c)) show a significant potential energy increase, which amounts to over 100 kJ/mol and is attained in around 15 fs.

Table 1 :
e most probable contributions of TD-DFT excited states (singlets).