Axel Koslowski

Implementation of a general multireference configuration interaction procedure with analytic gradients in a semiempirical context using the graphical unitary group approach

The graphical unitary group approach has been applied in an efficient implementation of a general multireference configuration interaction (MRCI) method for use with small active molecular orbital spaces in a semiempirical framework. Gradients can be computed analytically for molecular orbitals from a closed‐shell or a half‐electron open‐shell Hartree‐Fock calculation. CPU times for single point energy and gradient calculations are reported. The code allows MRCI geometry optimizations of large molecules, as illustrated for the singlet ground state and the four lowest triplet states of fullerene C76.

Cooperating Dinitrogen and Phenyl Rotations in trans-Azobenzene Photoisomerization

Semiempirical OM2/MRCI surface-hopping simulations have been performed to study the trans-to-cis photoisomerization of azobenzene upon excitation to the S1 state. The decay dynamics to the ground state shows an oscillatory pattern that can be attributed to an out-of-plane rotation of the N2 moiety. The reaction is thus initially driven by N2 rotation which triggers phenyl rotations around the C-N bonds. The cis isomer is produced most effectively when the phenyl rings rotate in phase. Mode-specific excitations cause variations in the computed decay times and product yields.

Semiempirical Quantum Chemical Calculations Accelerated on a Hybrid Multicore CPU-GPU Computing Platform.

In this work, we demonstrate that semiempirical quantum chemical calculations can be accelerated significantly by leveraging the graphics processing unit (GPU) as a coprocessor on a hybrid multicore CPU-GPU computing platform. Semiempirical calculations using the MNDO, AM1, PM3, OM1, OM2, and OM3 model Hamiltonians were systematically profiled for three types of test systems (fullerenes, water clusters, and solvated crambin) to identify the most time-consuming sections of the code. The corresponding routines were ported to the GPU and optimized employing both existing library functions and a GPU kernel that carries out a sequence of noniterative Jacobi transformations during pseudodiagonalization. The overall computation times for single-point energy calculations and geometry optimizations of large molecules were reduced by one order of magnitude for all methods, as compared to runs on a single CPU core.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks for Ground-State Properties

The semiempirical orthogonalization-corrected OMx methods (OM1, OM2, and OM3) go beyond the standard MNDO model by including additional interactions in the electronic structure calculation. When augmented with empirical dispersion corrections, the resulting OMx-Dn approaches offer a fast and robust treatment of noncovalent interactions. Here we evaluate the performance of the OMx and OMx-Dn methods for a variety of ground-state properties using a large and diverse collection of benchmark sets from the literature, with a total of 13035 original and derived reference data. Extensive comparisons are made with the results from established semiempirical methods (MNDO, AM1, PM3, PM6, and PM7) that also use the NDDO (neglect of diatomic differential overlap) integral approximation. Statistical evaluations show that the OMx and OMx-Dn methods outperform the other methods for most of the benchmark sets.

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.

Nonequilibrium H/D isotope effects from trajectory-based nonadiabatic dynamics.

Ground-state equilibrium kinetic isotope effects can be treated well in the framework of transition state theory, whereas excited-state nonequilibrium isotope effects are theoretically less explored. In this article we show for the first time that trajectory-based nonadiabatic dynamics simulations are able to reproduce experimental values for nonequilibrium H/D isotope effects in excited-state processes. We use high-level electronic structure calculations (MS-CASPT2, DFT/MRCI, and TDDFT) and full-dimensional OM2/MRCI-based nonadiabatic dynamics simulations to study the ultrafast intramolecular excited-state proton transfer (ESIPT) and the subsequent deactivation of 7-(2-pyridyl)indole (7PyIn) and its deuterated analogue (7PyIn-D). We evaluate a total of 1367 surface-hopping trajectories to establish the differences in the dynamical behavior of 7PyIn and 7PyIn-D. The computed H/D isotope effects for ESIPT and excited-state decay are consistent with recent experimental results from femtosecond pump-probe resonance-enhanced multiphoton ionization spectroscopy. We also analyze the influence of temperature fluctuations in the initially prepared sample on the photodynamics of 7PyIn and 7PyIn-D.

Photochemical dynamics of E-iPr-furylfulgide

As an important theoretical step towards unraveling the mechanistic details of the photochemical switching processes in molecules of the fulgide type, we carried out a large-scale, full-dimensional computational study of the ring closure reaction of E-iPr-furylfulgide. Simulated static UV spectra and femtosecond transient spectra are in good agreement with their experimental counterparts. Using surface-hopping photodynamics simulations, we identify three major de-excitation pathways and their interplay. The dominant photochemical pathway (70% of the trajectories) allows for ring closure, while the two minor pathways involve E-Z double bond isomerization rather than cyclization. The relative abundance of the pathways is rationalized by arguments linking structure with dynamics. It should be emphasized, however, that the distinction into three pathways is only a simplified interpretational model, since the actual dynamical trajectories do not strictly follow these idealized pathways but often show mixed behaviour, evolving along two or three of them during the course of the simulation.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks of Electronically Excited States

The semiempirical orthogonalization-corrected OMx methods have recently been shown to perform well in extensive ground-state benchmarks. They can also be applied to the computation of electronically excited states when combined with a suitable multireference configuration interaction (MRCI) treatment. We report on a comprehensive evaluation of the performance of the OMx/MRCI methods for electronically excited states. The present benchmarks cover vertical excitation energies, excited-state equilibrium geometries (including an analysis of significant changes between ground- and excited-state geometries), minimum-energy conical intersections, ground- and excited-state zero-point vibrational energies, and 0-0 transition energies for a total of 520 molecular structures and 412 excited states. For comparison, we evaluate the TDDFT/B3LYP method for all benchmark sets, and the CC2, MRCISD, and CASPT2 methods for some of them. We find that the current OMx/MRCI methods perform reasonably well for many of the excited-state properties. However, in comparison to the first-principles methods, there are also a number of shortcomings that should be addressed in future developments.

Fluorophores as Optical Sensors for Local Forces

The main aim of this study is to investigate correlations between the impact of an external mechanical force on the molecular framework of fluorophores and the resultant changes in their fluorescence properties. Taking into account previous theoretical studies, we designed a suitable custom-tailored oligoparaphenylenevinylene derivative (OPV5) with a twisted molecular backbone. Thin foils made of PVC doped with 100 nM OPV were prepared. By applying uniaxial force, the foils were stretched and three major optical effects were observed simultaneously. First, the fluorescence anisotropy increased, which indicates a reorientation of the fluorophores within the matrix. Second, the fluorescence lifetime decreased by approximately 2.5% (25 ps). Finally, we observed an increase in the emission energy of about 0.2% (corresponding to a blue-shift of 1.2 nm). In addition, analogous measurements with Rhodamine 123 as an inert reference dye showed only minor effects, which can be attributed to matrix effects due to refractive index changes. To relate the observed spectroscopic changes to the underlying changes in molecular properties, quantum-chemical calculations were also performed. Semiempirical methods had to be used because of the size of the OPV5 chromophore. Two conformers of OPV5 (C(2) and C(i) symmetry) were considered and both gave very similar results. Both the observed blue-shift of fluorescence and the reduced lifetime of OPV5 under tensile stress are consistent with the results of the semiempirical calculations. Our study proves the feasibility of fluorescence-based local force probes for polymers under tension. Improved optical sensors of this type should in principle be able to monitor local mechanical stress in transparent samples down to the single-molecule level, which harbors promising applications in polymer science and nanotechnology.

Enhanced E → Z photoisomerisation in 2-aminoazobenzene

Theoretical non-adiabatic dynamics simulations at the semiempirical OM2/MRCI level reveal an enhanced E → Z photoisomerisation for 2-aminoazobenzene. The E isomer contains an intramolecular hydrogen bond which weakens the azo N–N bond and thus facilitates isomerisation.