Vertical Excitation Energies and Lifetimes of the Two Lowest Singlet Excited States of Cytosine, 5-Aza-cytosine, and the Triazine Family: Quantum Mechanics-Molecular Mechanics Studies.

A swarm of semi-classical quantum mechanics/molecular mechanics molecular-dynamics simulations where OM2/MNDO is combined with the Gromacs program for consideration of explicit water is performed, solving the time-dependent Schrödinger equation in each step of the trajectories together with the Tully's fewest switches algorithm. Within this stochastic treatment, time dependent probabilities of the three lowest electronic states are determined. The fact that nucleobases are quickly deactivated is confirmed in the cytosine case where our best lifetime estimation is τ1=0.82 ps for the model with 100 water molecules with the SPCE force field and a time step of 0.1 fs. Lifetimes of the remaining molecules are visibly longer: 5-azacytosine, 2,4-diamino-1,3,5-triazine (DT), and 2,4,6-triamino-1,3,5-triazine (TT) molecules have an S1 → S0 de-excitation time of slightly above 10 ps. The lifetimes of the triazine family increases with the increasing number of exocyclic amino groups, that is, s-triazine < 2-amino-1,3,5-triazine < DT < TT. This can be explained by a higher mobility of the carbon-bonded hydrogen atoms in comparison with heavier amino groups since their movement is slowed down due to a substantially higher mass than hydrogen atoms, which can easier reach the out-of-plane positions required in the conical intersection structures. Moreover, bulkier NH2 ligands suffer due to greater friction caused by the surrounding water environment. These mechanical aspects caused a change in the explored lifetime dependences in comparison with our previous gas-phase study.

Newton-X Platform: New Software Developments for Surface Hopping and Nuclear Ensembles.

Newton-X is an open-source computational platform to perform nonadiabatic molecular dynamics based on surface hopping and spectrum simulations using the nuclear ensemble approach. Both are among the most common methodologies in computational chemistry for photophysical and photochemical investigations. This paper describes the main features of these methods and how they are implemented in Newton-X. It emphasizes the newest developments, including zero-point-energy leakage correction, dynamics on complex-valued potential energy surfaces, dynamics induced by incoherent light, dynamics based on machine-learning potentials, exciton dynamics of multiple chromophores, and supervised and unsupervised machine learning techniques. Newton-X is interfaced with several third-party quantum-chemistry programs, spanning a broad spectrum of electronic structure methods.

Theoretical Investigation of the Effect of Alkylation and Bromination on Intersystem Crossing in BODIPY-Based Photosensitizers.

Halogenated and alkylated BODIPY derivatives are reported as suitable candidates for their use as photosensitizers in photodynamic therapy due to their efficient intersystem crossing (ISC) between states of different spin multiplicities. Spin-orbit couplings (SOCs) are evaluated using an effective one-electron spin-orbit Hamiltonian for brominated and alkylated BODIPY derivatives to investigate the quantitative effect of alkyl and bromine substituents on ISC. BODIPY derivatives containing bromine atoms have been found to have significantly stronger SOCs than alkylated BODIPY derivatives outside the Frank-Condon region while they are nearly the same at local minima. Based on calculated time-dependent density functional theory (TD-DFT) vertical excitation energies and SOCs, excited-state dynamics of three BODIPY derivatives were further explored with TD-DFT surface hopping molecular dynamics employing a simple accelerated approach. Derivatives containing bromine atoms have been found to have very similar lifetimes, which are much shorter than those of the derivatives possessing just the alkyl moieties. However, both bromine atoms and alkyl moieties reduce the HOMO/LUMO gap, thus assisting the derivatives to behave as efficient photosensitizers.

Photophysics of BODIPY-Based Photosensitizer for Photodynamic Therapy: Surface Hopping and Classical Molecular Dynamics.

Halogenated BODIPY derivatives are emerging as important candidates for photodynamic therapy of cancer cells due to their high triplet quantum yield. We probed fundamental photophysical properties and interactions with biological environments of such photosensitizers. To this end, we employed static TD-DFT quantum chemical calculations as well as TD-DFT surface hopping molecular dynamics on potential energy surfaces resulting from the eigenstates of the total electronic Hamiltonian including the spin-orbit (SO) coupling. Matrix elements of an effective one-electron spin-orbit Hamiltonian between singlet and triplet configuration interaction singles (CIS) auxiliary wave functions are calculated using a new code capable of dealing with singlets and both restricted and unrestricted triplets built up from up to three different and independent sets of (singlet, alpha, and beta) molecular orbitals. The interaction with a biological environment was addressed by using classical molecular dynamics (MD) in a scheme that implicitly accounts for electronically excited states. For the surface hopping trajectories, an accelerated MD approach was used, in which the SO couplings are scaled up, to make the calculations computationally feasible, and the lifetimes are extrapolated back to unscaled SO couplings. The lifetime of the first excited singlet state estimated by semiclassical surface hopping simulations is 139 ± 75 ps. Classical MD demonstrates that halogenated BODIPY in the ground state, in contrast to the unsubstituted one, is stable in the headgroup region of minimalistic cell membrane models, and while in the triplet state, the molecule relocates to the membrane interior ready for further steps of photodynamic therapy.

Orientation of Laurdan in Phospholipid Bilayers Influences Its Fluorescence: Quantum Mechanics and Classical Molecular Dynamics Study.

Fluidity of lipid membranes is known to play an important role in the functioning of living organisms. The fluorescent probe Laurdan embedded in a lipid membrane is typically used to assess the fluidity state of lipid bilayers by utilizing the sensitivity of Laurdan emission to the properties of its lipid environment. In particular, Laurdan fluorescence is sensitive to gel vs liquid⁻crystalline phases of lipids, which is demonstrated in different emission of the dye in these two phases. Still, the exact mechanism of the environment effects on Laurdan emission is not understood. Herein, we utilize dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylcholine (DOPC) lipid bilayers, which at room temperature represent gel and liquid⁻crystalline phases, respectively. We simulate absorption and emission spectra of Laurdan in both DOPC and DPPC bilayers with quantum chemical and classical molecular dynamics methods. We demonstrate that Laurdan is incorporated in heterogeneous fashion in both DOPC and DPPC bilayers, and that its fluorescence depends on the details of this embedding.

A new approach to molecular dynamics with non-adiabatic and spin-orbit effects with applications to QM/MM simulations of thiophene and selenophene.

We present surface hopping dynamics on potential energy surfaces resulting from the spin-orbit splitting, i.e., surfaces corresponding to the eigenstates of the total electronic Hamiltonian including the spin-orbit coupling. In this approach, difficulties arise because of random phases of degenerate eigenvectors and possibility of crossings of the resulting mixed states. Our implementation solves these problems and allows propagation of the coefficients both in the representation of the spin free Hamiltonian and directly in the "diagonal representation" of the mixed states. We also provide a detailed discussion of the state crossing and point out several peculiarities that were not mentioned in the previous literature. We also incorporate the effect of the environment via the quantum mechanics/molecular mechanics approach. As a test case, we apply our methodology to deactivation of thiophene and selenophene in the gas phase, ethanol solution, and bulk liquid phase. First, 100 trajectories without spin-orbit coupling have been calculated for thiophene starting both in S1 and S2 states. A subset of 32 initial conditions starting in the S2 state was then used for gas phase simulations with spin-orbit coupling utilizing the 3-step integrator of SHARC, our implementation of the 3-step propagator in Newton-X and two new "one-step" approaches. Subsequently, we carried out simulations in ethanol solution and bulk liquid phase for both thiophene and selenophene. For both molecules, the deactivation of the S2 state proceeds via the ring opening pathway. The total population of triplet states reaches around 15% and 40% after 80 fs for thiophene and selenophene, respectively. However, it only begins growing after the ring opening is initiated; hence, the triplet states do not directly contribute to the deactivation mechanism. For thiophene, the resulting deactivation lifetime of the S2 state was 68 fs in the gas phase, 76 fs in ethanol solution, and 78 fs in the liquid phase, in a good agreement with the experimental value of 80 fs (liquid phase). For selenophene, the obtained S2 lifetime was 60 fs in the gas phase and 62 fs for both ethanol solution and liquid phase. The higher rate of intersystem crossing to the triplet states in selenophene is likely the reason for the lower fluorescence observed in selenium containing polymer compounds.

Surface Hopping Dynamics with Correlated Single-Reference Methods: 9H-Adenine as a Case Study.

Surface hopping dynamics methods using the coupled cluster to approximated second order (CC2), the algebraic diagrammatic construction scheme to second order (ADC(2)), and the time-dependent density functional theory (TDDFT) were developed and implemented into the program system Newton-X. These procedures are especially well-suited to simulate nonadiabatic processes involving various excited states of the same multiplicity and the dynamics in the first excited state toward an energetic minimum or up to the region where a crossing with the ground state is found. 9H-adenine in the gas phase was selected as the test case. The results showed that dynamics with ADC(2) is very stable, whereas CC2 dynamics fails within 100 fs, because of numerical instabilities present in the case of quasi-degenerate excited states. ADC(2) dynamics correctly predicts the ultrafast character of the deactivation process. It predicts that C2-puckered conical intersections should be the preferential pathway for internal conversion for low-energy excitation. C6-puckered conical intersection also contributes appreciably to internal conversion, becoming as important as C2-puckered for high-energy excitations. In any case, H-elimination plays only a minor role. TDDFT based on a long-range corrected functional fails to predict the ultrafast deactivation. In the comparison with several other methods previously used for dynamics simulations of adenine, ADC(2) has the best performance, providing the most consistent results so far.

Nonadiabatic molecular dynamics study of the cis-trans photoisomerization of azobenzene excited to the S1 state.

Ab initio nonadiabatic dynamics simulations of cis-to-trans isomerization of azobenzene upon S(1) (n-π*) excitation are carried out employing the fewest-switches surface hopping method. Azobenzene photoisomerization occurs purely as a rotational motion of the central CNNC moiety. Two nonequivalent rotational pathways corresponding to clockwise or counterclockwise rotation are available. The course of the rotational motion is strongly dependent on the initial conditions. The internal conversion occurs via an S(0)/S(1) crossing seam located near the midpoint of both of these rotational pathways. Based on statistical analysis, it is shown that the occurrence of one or other pathway can be completely controlled by selecting adequate initial conditions.