Chaoyuan Zhu

Non-Born-Oppenheimer trajectories with self-consistent decay of mixing

A semiclassical trajectory method, called the self-consistent decay of mixing (SCDM) method, is presented for the treatment of electronically nonadiabatic dynamics. The SCDM method is a modification of the semiclassical Ehrenfest (SE) method (also called the semiclassical time-dependent self-consistent-field method) that solves the problem of unphysical mixed final states by including decay-of-mixing terms in the equations for the evolution of the electronic state populations. These terms generate a force, called the decoherent force (or dephasing force), that drives the electronic component of each trajectory toward a pure state. Results for several mixed quantum-classical methods, in particular the SCDM, SE, and natural-decay-of-mixing methods and several trajectory surface hopping methods, are compared to the results of accurate quantum mechanical calculations for 12 cases involving five different fully dimensional triatomic model systems. The SCDM method is found to be the most accurate of the methods tested. The method should be useful for the simulation of photochemical reactions.

Intersystem Crossing Pathway in Quinoline–Pyrazole Isomerism: A Time-Dependent Density Functional Theory Study on Excited-State Intramolecular Proton Transfer

The dynamics of the excited-state intramolecular proton-transfer (ESIPT) reaction of quinoline-pyrazole (QP) isomers, designated as QP-I and QP-II, has been investigated by means of time-dependent density functional theory (TDDFT). A lower barrier has been found in the potential energy curve for the lowest singlet excited state (S1) along the proton-transfer coordinate of QP-II compared with that of QP-I; however, this is at variance with a recent experimental report [J. Phys. Chem. A 2010, 114, 7886-7891], in which the authors proposed that the ESIPT reaction would only proceed in QP-I due to the absence of a PT emission for QP-II. Therefore, several deactivating pathways have been investigated to determine whether fluorescence quenching occurs in the PT form of QP-II (PT-II). The S1 state of PT-II has nπ* character, which is a well-known dark state. Moreover, the energy gap between the S1 and T2 states is only 0.29 eV, implying that an intersystem crossing (ISC) process would occur rapidly following the ESIPT reaction. Therefore, it is demonstrated that the ESIPT could successfully proceed in QP-II and that the PT emission would be quenched by the ISC process.

Algorithmic decoherence time for decay-of-mixing non–Born–Oppenheimer dynamics

The performance of an analytical expression for algorithmic decoherence time is investigated for non-Born-Oppenheimer molecular dynamics. There are two terms in the function that represents the dependence of the decoherence time on the system parameters; one represents decoherence due to the quantum time-energy uncertainty principle and the other represents a back reaction from the decoherent force on the classical trajectory. We particularly examine the question of whether the first term should dominate. Five one-dimensional two-state model systems that represent limits of multidimensional nonadiabatic dynamics are designed for testing mixed quantum-classical methods and for comparing semiclassical calculations with exact quantum calculations. Simulations are carried out with the semiclassical Ehrenfest method (SE), Tullys fewest switch version (TFS) of the trajectory surface hopping method, and the decay-of-mixing method with natural switching, coherent switching (CSDM), and coherent switching with reinitiation (CSDM-D). The CSDM method is demonstrated to be the most accurate method, and it has several desirable features: (i) It behaves like the representation-independent SE method in the strong nonadiabatic coupling regions; (ii) it behaves physically like the TFS method in noninteractive region; and (iii) the trajectories are continuous with continuous momenta. The CSDM method is also demonstrated to balance coherence well with decoherence, and the results are nearly independent of whether one uses the adiabatic or diabatic representation. The present results provide new insight into the formulation of a physically correct decoherence time to be used with the CSDM method for non-Born-Oppenheimer molecular dynamic simulations.

C/B codoping effect on band gap narrowing and optical performance of TiO2 photocatalyst: a spin-polarized DFT study

The electronic and optical properties of several possible C/B-codoped models of anatase and rutile TiO2 have been investigated systematically using spin-polarized density functional theory calculations. Our calculated results show that the synergistic effect of C/B codoping can induce some hybridized energy states appearing in the forbidden gap and the band gap has a narrowing in anatase and rutile TiO2, which can lead to a decrease of the photon excitation energy and an obvious red-shift of the optical absorption edge. These results lead to an excellent photocatalytic activity in C/B-codoped TiO2. Moreover, with the increase of C and B impurities’ concentration in anatase and rutile TiO2, we find that the intensity of impurity states has a strengthening in the band gap, the position of impurity states changes, and the visible-light absorption performance improves gradually.

Franck-Condon simulation of vibrationally resolved optical spectra for zinc complexes of phthalocyanine and tetrabenzoporphyrin including the Duschinsky and Herzberg-Teller effects

High resolved absorption and fluorescence spectra of zinc complexes of phthalocyanine (ZnPc) and tetrabenzoporphyrin (ZnTBP) in the region of Q states were reported. Few theoretical investigations were performed to simulate the well-resolved spectra and assigned the vibrational bands of the large molecules, especially for high symmetrical characteristic molecules, on account of the difficulties to optimize the excited states and analyze a large number of final vibrational-normal modes. In the present work, the S(0) ↔ S(1) absorption and fluorescence spectra (that is, the Q band) of ZnPc and ZnTBP were simulated using time-dependent density functional theory with the inclusions of Duschinsky and Herzberg-Teller contributions to the electronic transition dipole moments. The theoretical results provide a good description of the optical spectra and are proved to be in excellent agreement with experimental spectra in inert-gas matrices or in supersonic expansion. This study focused attentions on the optical spectral similarities and contrasts between ZnPc and ZnTBP, in particular the noticeable Duschinsky and Herzberg-Teller effects on the high-resolved absorption and fluorescence spectra were considered. Substitution of meso-tetraaza on the porphyrin macrocycle framework could affect the ground state geometry and alter the electron density distributions, the orbital energies that accessible in the Q band region of the spectrum. The results were used to help interpret both the nature of the electronic transitions in Q band region, and the spectral discrepancies between phthalocyanine and porphyrin systems.

Orientation hydrogen-bonding effect on vibronic spectra of isoquinoline in water solvent: Franck-Condon simulation and interpretation

The excited-state orientation hydrogen-bonding dynamics, and vibronic spectra of isoquinoline (IQ) and its cationic form IQc in water have been investigated at the time-dependent density functional theory quantum chemistry level plus Franck-Condon simulation and interpretation. The excited-state orientation hydrogen bond strengthening has been found in IQ:H2O complex due to the charge redistribution upon excitation; this is interpreted by simulated 1:1 mixed absorption spectra of free IQ and IQ:H2O complex having best agreement with experimental results. Conversely, the orientation hydrogen bond in IQc:H2O complex would be strongly weakening in the S1 state and this is interpreted by simulated absorption spectra of free IQc having best agreement with experimental results. By performing Franck-Condon simulation, it reveals that several important vibrational normal modes with frequencies about 1250 cm-1 involving the wagging motion of the hydrogen atoms are very sensitive to the formation of the orientation hydrogen bond for the IQ/IQc:H2O complex and this is confirmed by damped Franck-Condon simulation with free IQ/IQc in water. However, the emission spectra of the IQ and IQc in water have been found differently. Upon the excitation, the simulated fluorescence of IQ in water is dominated by the IQ:H2O complex; thus hydrogen bond between IQ and H2O is much easier to form in the S1 state. While the weakened hydrogen bond in IQc:H2O complex is probably cleaved upon the laser pulse because the simulated emission spectrum of the free IQc is in better agreement with the experimental results.

A TDDFT study on the excited-state double proton transfer reaction of 8-hydroxyquinoline along a hydrogen-bonded bridge

The mechanism of the excited-state double proton transfer (ESDPT) reactions of 8-hydroxyquinoline (8HQ) along three types of hydrogen-bonded bridges—ammonia (NH3), water (H2O) and acetic acid (AcOH) has been investigated by using time-dependent density functional theory. Based on the analysis of hydrogen bond strengths and the excited-state potential energy surfaces (PES) along the proton transfer coordinates, it is concluded that the hydrogen bonds play the key role in the excited-state multiple proton transfer reaction. Moreover, three different concerted mechanisms have been found in 8HQ·NH3, 8HQ·H2O and 8HQ·AcOH complexes. Upon photoexcitation, the deprotonation of a hydroxyl group in 8HQ would occur first in 8HQ·NH3 due to the stronger hydrogen bond –OH⋯NH3. On the contrary, the hydrogen bond –OH⋯O–C–OH is much weaker in the 8HQ·AcOH complex, and this leads to a fast protonation of –N– in the 8HQ moiety. For the 8HQ·H2O complex, the hydrogen bond strengths are almost the same, so that both protons would transfer simultaneously in a symmetrical and concerted fashion.

A phenothiazine-based colorimetric chemodosimeter for the rapid detection of cyanide anions in organic and aqueous media

A novel, phenothiazine-based chemodosimeteric probe PCP 1 has been developed and studied as a cyanide selective indicator in organic and aqueous media. Complete colour bleaching was observed due to nucleophilic addition of cyanide to the tricyanovinyl moiety of PCP 1, which results in the disruption of the extended conjugate and turn-off the intramolecular charge transfer process.

Exploring the role of varied-length spacers in charge transfer: a theoretical investigation on pyrimidine-bridged porphyrin dyes

A series of dyes based on a porphyrin donor and a cyanoacrylic acid anchor/acceptor group for solar cell application are investigated with regards to varied-length π-spacers affecting the photo-to-electric conversion efficiency (PCE). Investigations are firstly performed on three porphyrin sensitizers with 1–3 conjugated phenylethynyl (PE) units, which have experimentally proved that the efficiency of power conversion decreases systematically with increasing spacer length. The distances and amounts of charge transfer after photoexcitation are calculated. In the PE bridged porphyrin dyes, the calculated electron injection driving forces and the regeneration driving forces gradually decrease as the distance of the π-spacer increases. Our theoretical calculations can reproduce well the experimental conclusion, showing that the photo-to-electric efficiency has a strong distance dependence for the electron-rich phenyl spacer. Then we replace the phenyl group with a pyrimidyl (PM) group to uncover how the characteristics of the π-spacer affect the performance of optical absorption, charge separation, and the regeneration process, to further improve the power conversion. We find that the adoption of electron-deficient pyrimidyl can break and even remove the distance dependence of the π-spacer. Some integral factors affecting the dye performance, such as short-circuit photocurrent, open-circuit voltage and charge collection efficiency are analyzed. It would help to interpret what role the electron deficient π-spacers with varied lengths will play and how they are expected to behave in the performance of sensitizers. In this regard, this study presents us with a promising way to design novel functional dyes and to utilize the potential advantages of the lengthy spacer dyes.

Landscapes of four-enantiomer conical intersections for photoisomerization of stilbene: CASSCF calculation.

The photoisomerization of cis- and trans-stilbene through conical intersections (CI) is mainly governed by four dihedral angles around central C═C double bonds. The two of them are C-C═C-C and H-C═C-H dihedral angles that are found to form a mirror rotation coordinate, and the mirror plane appears at the two dihedral angles equal to zeros with which the middle state is defined through partial optimization. There exist the first-type of hula-twist-CI enantiomers, the second-type of hula-twist-CI enantiomers, the first-type of one-bond-flip-CI enantiomers, and the second type of one-bond-flip-CI enantiomers as well as cis-enantiomers and trans-enantiomers with respect to this mirror plane. The complete active space self-consistent field method is employed to calculate minimum potential energy profile along the mirror rotation coordinate for each enantiomers, and it is found that the left-hand manifold and the right-hand manifold of potential energy surfaces can be energetically transferred via photoisomerization. Furthermore, two-dimensional potential energy surfaces in terms of the branching plane g-h coordinates are constructed at vicinity of each conical intersection, and the landscapes of conical intersections show distinct feature, and in excited-state four potential wells separated in different section of g-h plane related to different conical intersections which indicate different photoisomerization pathways.