Bin-Bin Xie

Photoprotection Mechanism of p-Methoxy Methylcinnamate: A CASPT2 Study

p-Methoxy methylcinnamate (p-MMC) shares the same molecular skeleton with octyl methoxycinnamate sunscreen. It is recently found that adding one water to p-MMC can significantly enhance the photoprotection efficiency. However, the physical origin is elusive. Herein we have employed multireference complete active space self-consistent field (CASSCF) and multistate complete active-space second-order perturbation (MS-CASPT2) methods to scrutinize the photophysical and photochemical mechanism of p-MMC and its one-water complex p-MMC-W. Specifically, we optimize the stationary-point structures on the (1)ππ*, (1)nπ*, and S0 potential energy surfaces to locate the (1)ππ*/S0 and (1)ππ*/(1)nπ* conical intersections and to map (1)ππ* and (1)nπ* excited-state relaxation paths. On the basis of the results, we find that, for the trans p-MMC, the major (1)ππ* deactivation path is decaying to the dark (1)nπ* state via the in-plane (1)ππ*/(1)nπ* crossing point, which only need overcome a small barrier of 2.5 kcal/mol; the minor one is decaying to the S0 state via the (1)ππ*/S0 conical intersection induced by out-of-plane photoisomerization. For the cis p-MMC, these two decay paths are comparable (1)ππ* deactivation paths: one is decaying to the dark (1)nπ* state via the (1)ππ*/(1)nπ* crossing point, and the second is decaying to the ground state via the (1)ππ*/S0 conical intersection. One-water hydration stabilizes the (1)ππ* state and meanwhile destabilizes the (1)nπ* state. As a consequence, the (1)ππ* deactivation path to the dark (1)nπ* state is heavily inhibited. The related barriers are increased to 5.8 and 3.3 kcal/mol for the trans and cis p-MMC-W, respectively. In comparison, the barriers associated with the photoisomerization-induced (1)ππ* decay paths are reduced to 2.5 and 1.3 kcal/mol for the trans and cis p-MMC-W. Therefore, the (1)ππ* decay paths to the S0 state are dominant relaxation channels when adding one water molecule. Finally, the present work contributes a lot of knowledge to understanding the photoprotection mechanism of methylcinnamate derivatives.

Surface-hopping dynamics simulations of malachite green: a triphenylmethane dye.

Malachite green is a typical triphenylmethane dye widely used in fundamental and industrial research; however, its excited-state relaxation dynamics remains elusive. In this work we simulate its photodynamics from the S2 and S1 states using the fewest-switches surface-hopping scheme. In the S2 photodynamics, the system first relaxes to the S2 minimum, which immediately hops to the S1 state via an S2/S1 conical intersection. In the S1 state, 90% trajectories evolve into a structurally symmetric S1 minimum; the remaining ones proceed toward two propeller-like S1 minima. Two kinds of S1 minima then decay to the S0 state via the S1/S0 conical intersections. The S1 photodynamics is overall similar to the S1 excited-state dynamics as a result of the ultrafast S2 → S1 internal conversion in the S2 photodynamics, but the weights of the trajectories that decay to the S0 state via three different S1/S0 conical intersections are variational. Moreover, the S2 relaxation dynamics mainly happens in a concerted synchronous rotation of three phenyl rings. In comparison, in the S1 relaxation dynamics, the rotations of two aminophenyl rings can proceed in the same and opposite directions. In certain trajectories, only the rotation of an aminophenyl ring is active. On the basis of the results, the S2 and S1 excited-state lifetimes of malachite green in vacuo are calculated to be 424 fs and 1.2 ps, respectively. The present work provides important mechanistic insights for similar triphenylmethane dyes.

The Position of the N Atom Plays a Significant Role for Excited-State Decay of Heterocycles

We have employed combined electronic structure calculations and nonadiabatic dynamics simulations to study the S1 radiationless deactivation mechanism of pyrazole. In terms of MS-CASPT2 computed results, we propose that the 1πσ* state-driven nonadiabatic N-N dissociation is a major relaxation path; the ring-puckering deformation path as well as the 1πσ* state-driven N-H dissociation are less favorable. This excited-state decay mechanism is supported by MS-CASPT2 nonadiabatic dynamics simulations. The present study demonstrates that pyrazole has a different excited-state radiationless deactivation mechanism compared with its structural isomer imidazole, in which the 1πσ* state-driven nonadiabatic N-H dissociation plays a more important role. However, such a channel is suppressed in pyrazole; instead, the 1πσ* state-driven nonadiabatic N-N dissociation is dominant.

Ab initio implementation of quantum trajectory mean-field approach and dynamical simulation of the N2CO photodissociation

In this work, the recently introduced quantum trajectory mean-field (QTMF) approach is implemented and employed to explore photodissociation dynamics of diazirinone (N2CO), which are based on the high-level ab initio calculation. For comparison, the photodissociation process has been simulated as well with the fewest-switches surface hopping (FSSH) and the ab initio multiple spawning (AIMS) methods. Overall, the dynamical behavior predicted by the three methods is consistent. The N2CO photodissociation at λ > 335 nm is an ultrafast process and the two C-N bonds are broken in a stepwise way, giving birth to CO and N2 as the final products in the ground state. Meanwhile, some noticeable differences were found in the QTMF, FSSH, and AIMS simulated time constants for fission of the C-N bonds, excited-state lifetime, and nonadiabatic transition ratios in different intersection regions. These have been discussed in detail. The present study provides a clear evidence that direct ab initio QTMF approach is one of the reliable tools for simulating nonadiabatic dynamics processes.

Multiple-State Nonadiabatic Dynamics Simulation of Photoisomerization of Acetylacetone with the Direct ab Initio QTMF Approach

In the present work, the quantum trajectory mean-field (QTMF) approach is numerically implemented by ab initio calculation at the level of the complete active space self-consistent field, which is used to simulate photoisomerization of acetylacetone at ∼265 nm. The simulated results shed light on the possible nonadiabatic pathways from the S2 state and mechanism of the photoisomerization. The in-plane proton transfer and the subsequent S2 → S1 transition through the E-S2/S1-1 intersection region is the predominant route to the S1 state. Meanwhile, rotational isomerization occurs in the S2 state, which is followed by internal conversion to the S1 state in the vicinity of the E-S2/S1-2 conical intersection. As a minor pathway, the direct S2 → S1 → S0 transition can take place via the E-S2/S1/S0 three-state intersection region. The rotamerization in the S1 state was determined to be the key step for formation of nonchelated enolic isomers. The final formation yield is predicted to be 0.57 within the simulated period. The time constant for the S2 proton transfer was experimentally inferred to be ∼70.0 fs in the gas phase and ∼50.0 fs in dioxane, acetonitrile, and n-hexane, which is well-reproduced by the present QTMF simulation. The S1 lifetime of 2.11 ps simulated here is in excellent agreement with the experimentally inferred values of 2.12, 2.13, and 2.25 ps in n-hexane, acetonitrile, and dioxane, respectively. The present study provides clear evidence that a direct ab initio QTMF approach is a reliable tool for simulating multiple-state nonadiabatic dynamics processes.

Short-time dynamics and decay mechanism of 2(1H)-pyridinone upon excitation to the light-absorbing S4(21𝝅𝝅*) state

The excited-state structural dynamics and the decay mechanism of 2(1H)-pyridinone (NHP) after excitation to the S4(21ππ*) light-absorbing state were studied using resonance Raman spectroscopy and complete-active space self-consistent field (CASSCF) calculations. The B-band absorption cross-section and the corresponding absolute resonance Raman cross-sections were simulated using a simple model based on time-dependent wave-packet theory. The geometric structures of the singlet electronic excited states and their curve-crossing points were optimized at the CASSCF level of theory. The obtained short-time structural dynamics in easy-to-visualize internal coordinates were then compared with the CASSCF-predicted structural-parameter changes of S4(21ππ*)/S3(21nπ*)-MIN, S4(21ππ*)/S1(11nπ*)-MIN, and S4(21ππ*)-MIN. Our results indicate that the initial population of NHP in the S4 state bifurcates in or near the Franck-Condon region, leading to two predominant (S4S3-MIN and S4S1-MIN) internal conversion pathways. The lowest-lying S2(11ππ*) excited state is finally formed via subsequent internal conversions S3(21nπ*)/S2(11ππ*)-MIN and S1(11nπ*)/S2(11ππ*)-MIN. The enol-keto tautomeric mechanism does not seem to play a role. The decay mechanism in the singlet realm is proposed.