Generic implementation of semi-analytical CI gradients for NDDO-type methods

Conformational Control of the Photodynamics of a Bilirubin Dipyrrinone Subunit: Femtosecond Spectroscopy Combined with Nonadiabatic Simulations

The photochemistry of bilirubin has been extensively studied due to its importance in the phototherapy of hyperbilirubinemia. In the present work, we investigated the ultrafast photodynamics of a bilirubin dipyrrinone subunit, vinylneoxanthobilirubic acid methyl ester. The photoisomerization and photocyclization reactions of its (E) and (Z) isomers were studied using femtosecond transient absorption spectroscopy and by multireference electronic structure theory, where the nonadiabatic dynamics was modeled with a Landau-Zener surface hopping technique. The following picture has emerged from the combined theoretical and experimental approach. Upon excitation, dipyrrinone undergoes a very fast vibrational relaxation, followed by an internal conversion on a picosecond time scale. The internal conversion leads either to photoisomerization or regeneration of the starting material. Further relaxation dynamics on the order of tens of picoseconds was observed in the ground state. The nonadiabatic simulations revealed a strong conformational control of the photodynamics. The ultrafast formation of a cyclic photochemical product from a less-populated conformer of the studied subunit was predicted by our calculations. We discuss the relevance of the present finding for the photochemistry of native bilirubin. The work has also pointed to the limits of semiclassical nonadiabatic simulations for simulating longer photochemical processes, probably due to the zero-point leakage issue.

metaFALCON: A Program Package for Automatic Sampling of Conical Intersection Seams Using Multistate Metadynamics

The multistate metadynamics for automatic exploration of conical intersection seams and systematic location of minimum energy crossing points in molecular systems and its implementation into the software package metaFALCON is presented. Based on a locally modified energy gap between two Born-Oppenheimer electronic states as a collective variable, multistate metadynamics trajectories are driven toward an intersection point starting from an arbitrary ground state geometry and are subsequently forced to explore the conical intersection seam landscape. For this purpose, an additional collective variable capable of distinguishing structures within the seam needs to be defined and an additional bias is introduced into the off-diagonal elements of an extended (multistate) electronic Hamiltonian. We demonstrate the performance of the algorithm on the examples of the 1,3-butadiene, benzene, and 9H-adenine molecules, where multiple minimum energy crossing points could be systematically located using the Wiener number or Cremer-Pople parameters as collective variables. Finally, with the example of 9H-adenine, we show that the multistate metadynamics potential can be used to obtain a global picture of a conical intersection seam. Our method can be straightforwardly connected with any ab initio or semiempirical electronic structure theory that provides energies and gradients of the respective electronic states and can serve for systematic elucidation of the role of conical intersections in the photophysics and photochemistry of complex molecular systems, thus complementing nonadiabatic dynamics simulations.

Excited state dynamics and time-resolved photoelectron spectroscopy of para-xylylene

We investigated the excited-state dynamics of para-xylylene using a combination of field-induced surface hopping (FISH) simulations and time-resolved ionisation experiments. Our simulations predict an ultrafast decay of the initially excited bright state (S2/S3) to the S1 state on a sub-100 fs time scale, followed by return to the ground state within ∼1 ps. This is accompanied by a transient change of the biradical character of the molecule, as monitored by calculating natural orbital occupation numbers. Specifically, the initially low biradicality is increased by electronic excitation as well as by vibrational activation. Experimentally, para-xylylene was generated by pyrolysis from [2,2]paracyclophane and excited with 266 nm radiation into the S2/S3 bright state. The subsequent dynamics were followed using ionisation as the probe step, with both mass spectra and photoelectron spectra recorded as a function of pump-probe delay. The observed decay of photoelectron and photoion intensities closely matches the theoretical predictions and is consistent with the sequential mechanism found in the simulations. This mechanism exhibits characteristic signatures in both time-resolved mass and photoelectron spectra, in particular in the appearance of fragment ions that are exclusively generated from the S1 state. This allows for a separation of the S2 and S1 dynamics in the photoelectron and mass spectra. An excellent agreement between the observed and the simulated ion signal is observed.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Theory, Implementation, and Parameters

Semiempirical orthogonalization-corrected methods (OM1, OM2, and OM3) go beyond the standard MNDO model by explicitly including additional interactions into the Fock matrix in an approximate manner (Pauli repulsion, penetration effects, and core–valence interactions), which yields systematic improvements both for ground-state and excited-state properties. In this Article, we describe the underlying theoretical formalism of the OMx methods and their implementation in full detail, and we report all relevant OMx parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. For a standard set of mostly organic molecules commonly used in semiempirical method development, the OMx results are found to be superior to those from standard MNDO-type methods. Parametrized Grimme-type dispersion corrections can be added to OM2 and OM3 energies to provide a realistic treatment of noncovalent interaction energies, as demonstrated for the complexes in the S22 and S66×8 test sets.

Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.

UV photodissociation of proline-containing peptide ions: insights from molecular dynamics

UV photodissociation of proline-containing peptide ions leads to unusual product ions. In this paper, we report laser-induced dissociation of a series of proline-containing peptides at 213 nm. We observe specific fragmentation pathways corresponding to the formation of (y-2), (a + 2) and (b + 2) fragment ions. This was not observed at 266 nm or for peptides which do not contain proline residues. In order to obtain insights into the fragmentation dynamics at 213 nm, a small peptide (RPK for arginine-proline-lysine) was studied both theoretically and experimentally. Calculations of absorption spectra and non-adiabatic molecular dynamics (MD) were made. Second and third excited singlet states, S(2) and S(3), lie close to 213 nm. Non-adiabatic MD simulation starting from S(2) and S(3) shows that these transitions are followed by C-C and C-N bond activation close to the proline residue. After this first relaxation step, consecutive rearrangements and proton transfers are required to produce unusual (y-2), (a + 2) and (b + 2) fragment ions. These fragmentation mechanisms were confirmed by H/D exchange experiments.

How shaped light discriminates nearly identical biochromophores

We present a general mechanism for successful discrimination of spectroscopically indistinguishable biochromophores by shaped light. For this purpose we use nonadiabatic dynamics in excited electronic states in the frame of the field-induced surface hopping method driven by the experimentally shaped laser fields. Our findings show that optimal laser fields drive low-frequency vibrational modes localized in the side chains of two biochromophores, thus selecting the parts of their potential energy surfaces characterized by different transition dipole moments leading to different ionization probabilities. The presented mechanism leads to selective fluorescence depletion which serves as a discrimination signal. Our findings offer a promising perspective for using optimally shaped laser pulses in bioanalytical applications by increasing the selectivity beyond the current capability.