Researcher profile

Oliver Kühn

Oliver Kühn contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2023arXiv

Formally Exact Simulations of Mesoscale Exciton Diffusion in a Photosynthetic Aggregate

The photosynthetic apparatus of plants and bacteria combine atomically precise pigment-protein complexes with dynamic membrane architectures to control energy transfer on the 10-100 nm length scales. Recently, synthetic materials have integrated photosynthetic antennae proteins to enhance exciton transport, though the influence of artificial packing on the excited-state dynamics in these biohybrid materials remains unclear. Here, we use the adaptive Hierarchy of Pure States (adHOPS) to perform a formally exact simulation of excitation energy transfer within artificial aggregates of light harvesting complex 2 (LH2) with a range of packing densities. We find that LH2 aggregates support a remarkable exciton diffusion length ranging from 100 nm at a biological packing density to 300 nm at the densest packing previously suggested in an artificial aggregate. The unprecedented scale of these calculations also underscores the efficiency with which adHOPS simulates excited-state processes in realistic molecular materials.

preprint2020arXiv

Exciton Transfer Using Rates Extracted From the "Hierarchical Equations of Motion''

Frenkel exciton population dynamics of an excitonic dimer is studied by comparing results from a quantum master equation (QME) involving rates from second-order perturbative treatment with respect to the excitonic coupling with non-perturbative results from ``Hierarchical Equations of Motion'' (HEOM). By formulating generic Liouville-space expressions for the rates, we can choose to evaluate them either via HEOM propagations or by applying cumulant expansion. The coupling of electronic transitions to bath modes is modeled either as overdamped oscillators for description of thermal bath components or as underdamped oscillators to account for intramolecular vibrations. Cases of initial nonequilibrium and equilibrium vibrations are discussed. In case of HEOM initial equilibration enters via a polaron transformation. Pointing out the differences between the nonequilibrium and equilibrium approach in the context of the projection operator formalism, we identify a further description, where the transfer dynamics is driven only by fluctuations without involvement of dissipation. Despite this approximation, also this approach can yield meaningful results in certain parameter regimes. While for the chosen model HEOM has no technical advantage for evaluation of the rate expressions compared to cumulant expansion, there are situations where only evaluation with HEOM is applicable. For instance, a separation of reference and interaction Hamiltonian via a polaron transformation to account for the interplay between Coulomb coupling and vibrational oscillations of the bath at the level of a second-order treatment can be adjusted for a treatment with HEOM.

preprint2016arXiv

A Network Approach to Unravel Correlated Ion Pair Dynamics in Protic Ionic Liquids. The Case of Triethylammonium Nitrate

The intermolecular interactions in the title compound are investigated using self-consistent charge density functional based tight binding molecular dynamics. Emphasis is put on the analysis of correlated motions of ion pairs using ideas of network theory. At equilibrium such correlations are not very pronounced on average. However, there exist sizeable local correlations for cases where two cations share the same anion via two NHO-hydrogen bonds. The effect of an external perturbation, which artificially introduces a sudden local heating of an NH-bond, is investigated using nonequilibrium molecular dynamics. Here, it is found that the average N-H bond vibrational relaxation time is about 5.3 ~ps. This energy redistribution is rather nonspecific with respect to the ion pairs and does not lead to long-range correlations spreading from the initially excited ion pair.

preprint2016arXiv

Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) Approach to the Correlated Exciton-Vibrational Dynamics in the FMO Complex

The coupled quantum dynamics of excitonic and vibrational degrees of freedom is investigated for high-dimensional models of the Fenna-Matthews-Olson (FMO) complex. This includes a seven and an eight-site model with 518 and 592 harmonic vibrational modes, respectively. The coupling between local electronic transitions and vibrations is described within the Huang-Rhys model using parameters that are obtained by discretization of an experimental spectral density. Different pathways of excitation energy flow are analyzed in terms of the reduced one-exciton density matrix, focussing on the role of vibrational and vibronic excitation. Distinct features due to both competing time scales of vibrational and exciton motion and vibronically-assisted transfer are observed. The question of the effect of initial state preparation is addressed by comparing the case of an instantaneous Franck-Condon excitation at a single site with that of a laser field excitation.

preprint2016arXiv

Towards an ab initio theory for metal L-edge soft X-ray spectroscopy of molecular aggregates

The Frenkel exciton model was adapted to describe X-ray absorption and resonant inelastic scattering spectra of polynuclear transition metal complexes by means of restricted active space self-consistent field method. The proposed approach allows to substantially decrease the requirements to computational resources if compared to a full supermolecular quantum chemical treatment. This holds true in particular in cases where the dipole approximation to the electronic transition charge density can be applied. The computational protocol was applied to the calculation of X-ray spectra of the hemin complex, which forms dimers in aqueous solution. The aggregation effects were found to be comparable to the spectral alterations due to the replacement of the axial ligand by solvent molecules.

preprint2016arXiv

Vibrational spectroscopy via the Caldeira-Leggett model with anharmonic system potentials

The Caldeira-Leggett (CL) model, which describes a system bi-linearly coupled to a harmonic bath, has enjoyed popularity in condensed phase spectroscopy owing to its utmost simplicity. However, the applicability of the model to cases with anharmonic system potentials, as it is required for the description of realistic systems in solution, is questionable due to the presence of the invertibility problem [J. Phys. Chem. Lett., \textbf{6}, 2722 (2015)] unless the system itself resembles the CL model form. This might well be the case at surfaces or in the solid regime, which we here confirm for a particular example of an iodine molecule in the atomic argon environment under high pressure. For this purpose we extend the recently proposed Fourier method for parameterizing linear generalized Langevin dynamics[J. Chem. Phys., \textbf{142}, 244110 (2015)] to the non-linear case based on the CL model and perform an extensive error analysis. In order to judge on the applicability of this model in advance, we give handy empirical criteria and discuss the effect of the potential renormalization term. The obtained results provide evidence that the CL model can be used for describing a potentially broad class of systems.

preprint2015arXiv

Explicit Correlated Exciton-Vibrational Dynamics of the FMO Complex

The coupled exciton-vibrational dynamics of a 3-site FMO model is investigated using the numerically exact multilayer multiconfiguration time-dependent Hartree approach. Thereby the vibrational mode specific coupling to local electronic transitions is adapted from a discretized experimental spectral density. The solution of the resulting time-dependent Schrödinger equation including three electronic and 450 vibrational degrees of freedom is analyzed in terms of excitonic populations and coherences. Emphasis is put onto the role of specific ranges of vibrational frequencies. It is observed that modes between 160 and 300 cm$^{-1}$ are responsible for the subpicosecond population and coherence decay. Further, it is found that a mean-field approach with respect to the vibrational degrees of freedom is not applicable.

preprint2015arXiv

Multi-reference approach to the calculation of photoelectron spectra including spin-orbit coupling

X-ray photoelectron spectra provide a wealth of information on the electronic structure. The extraction of molecular details requires adequate theoretical methods, which in case of transition metal complexes has to account for effects due to the multi-configurational and spin-mixed nature of the many-electron wave function. Here, the Restricted Active Space Self-Consistent Field method including spin-orbit coupling is used to cope with this challenge and to calculate valence and core photoelectron spectra. The intensities are estimated within the frameworks of the Dyson orbital formalism and the sudden approximation. Thereby, we utilize an efficient computational algorithm that is based on a biorthonormal basis transformation. The approach is applied to the valence photoionization of the gas phase water molecule and to the core ionization spectrum of the $\text{[Fe(H}_2\text{O)}_6\text{]}^{2+}$ complex. The results show good agreement with the experimental data obtained in this work, whereas the sudden approximation demonstrates distinct deviations from experiments.

preprint2015arXiv

On the Applicability of the Caldeira-Leggett Model to Condensed Phase Vibrational Spectroscopy

Formulating a rigorous system-bath partitioning approach remains an open issue. In this context the famous Caldeira-Leggett model that enables quantum and classical treatment of Brownian motion on equal footing has enjoyed popularity. Although this model is by any means a useful theoretical tool, its ability to describe anharmonic dynamics of real systems is often taken for granted. In this Letter we show that the mapping between a molecular system under study and the model cannot be established in a self-consistent way, unless the system part of the potential is taken effectively harmonic.Mathematically, this implies that the mapping is not invertible. This `invertibility problem' is not dependent on the peculiarities of particular molecular systems under study and is rooted in the anharmonicity of the system part of the potential.

preprint2015arXiv

Parametrizing linear generalized Langevin dynamics from explicit molecular dynamics simulations

Fundamental understanding of complex dynamics in many-particle systems on the atomistic level is of utmost importance. Often the systems of interest are of macroscopic size but can be partitioned into few important degrees of freedom which are treated most accurately and others which constitute a thermal bath. Particular attention in this respect attracts the linear generalized Langevin equation (GLE), which can be rigorously derived by means of a linear projection (LP) technique. Within this framework a complicated interaction with the bath can be reduced to a single memory kernel. This memory kernel in turn is parametrized for a particular system studied, usually by means of time-domain methods based on explicit molecular dynamics data. Here we discuss that this task is most naturally achieved in frequency domain and develop a Fourier-based parametrization method that outperforms its time-domain analogues. Very surprisingly, the widely used rigid bond method turns out to be inappropriate in general. Importantly, we show that the rigid bond approach leads to a systematic underestimation of relaxation times, unless the system under study consists of a harmonic bath bi-linearly coupled to the relevant degrees of freedom.

preprint2015arXiv

Quantum chemical study of the electronic properties of an Iridium-based photosensitizer bound to medium-sized silver clusters

The equilibrium structures and electronic excitation spectra of the Ir(III) photosensitizer Ir(ppy)$_2$(bpy)]$^+$ bound to medium-sized silver clusters Ag$_n$ ($n=$19, 20) are investigated using time-dependent density functional theory. The long-range corrected LC-BLYP approach is used with a system-specific range separation parameter. The weak physisorption of the hybrid complexes yield only small changes in the broadened absorption spectra of the hybrid system as compared with its constituents. However, the density of states as well as the fine structure of the spectra is strongly modified upon complexation. It is shown that the standard range separation parameter (0.47 bohr$^{-1}$) cannot predict these properties correctly and the optimized value of 0.16 bohr$^{-1}$ should be used instead.

preprint2015arXiv

Tuning Range-Separated Density Functional Theory for Photocatalytic Water Splitting Systems

We discuss the system-specific optimization of long-range separated density functional theory (DFT) for the prediction of electronic properties relevant for a photocatalytic cycle based on an Ir(III) photosensitizer (IrPS). Special attention is paid to the charge-transfer properties, which are of key importance for the photoexcitation dynamics, but and cannot be correctly described by means of conventional DFT. The optimization of the range-separation parameter using the $Δ$SCF method is discussed for IrPS including its derivatives and complexes with electron donors and acceptors used in photocatalytic hydrogen production. Particular attention is paid to the problems arising for a description of medium effects by means of a polarizable continuum model.

preprint2014arXiv

DFT-D Investigation of the Interaction Between Ir(III) Based Photosensitizers and Small Silver Clusters Ag$_n$ ($n$=2-20, 92)

A dispersion-corrected density functional theory study of the photosensitizer [Ir(ppy)$_2$(bpy)]$^+$ and its derivatives bound to silver clusters Ag$_n$ ($n$=2-20, 92) is performed. The goal is to provide a new system-specific set of $C_{\rm 6}^{}$ interaction parameters for Ag and Ir atoms. To this end a QM:QM scheme is employed using the PBE functional and RPA as well as MP2 calculation as a reference. The obtained $C_{\rm 6}^{}$ coefficients were applied to calculate dissociation curves of selected IrPS-Ag$_n$ complexes and binding energies of derivatives containing oxygen and sulphur as heteroatoms in the ligands. Comparing different $C_{\rm 6}^{}$ parameters it is concluded that RPA-based dispersion correction produces binding energies close to standard D2 and D3 models, whereas MP2-derived parameters overestimate these energies.

preprint2014arXiv

Exciton coupling induces vibronic hyperchromism in light-harvesting complexes

The recently suggested possibility that weak vibronic transitions can be excitonically enhanced in light-harvesting complexes is studied in detail. A vibronic exciton dimer model which includes ground state vibrations is investigated using multi-configuration time-dependent Hartree method with a parameter set typical to photosynthetic light-harvesting complexes. Absorption spectra are discussed in dependence on the Coulomb coupling, the detuning of site energies, and the number of vibrational mode. Calculations of the fluorescence spectra show that the spectral densities obtained from the low temperature fluorescence line narrowing measurements of light-harvesting systems need to be corrected for the exciton effects. For the J-aggregate configuration, as in most of the light-harvesting complexes, the true spectral density has larger amplitude than what is obtained from the measurement.

preprint2014arXiv

Interaction of Polar and Nonpolar Organic Pollutants with Soil Organic Matter: Sorption Experiments and Molecular Dynamics Simulation

The fate of organic pollutants in the environment is influenced by several factors including the type and strength of their interactions with soil components especially SOM. However, a molecular level answer to the question How organic pollutants interact with SOM? is lacking. In order to explore mechanisms of this interaction, we have developed a new SOM model followed by carrying out molecular dynamics (MD) simulations in parallel with sorption experiments. The new SOM model comprises free SOM functional groups (carboxylic acid and naphthalene) as well as SOM cavities (with two different sizes), representing the soil voids, containing the same SOM functional groups. To examine the effect of the hydrophobicity on the interaction, the organic pollutants hexachlorobenzene (HCB, non-polar) and sulfanilamide (SAA, polar) were considered. The experimental and the theoretical outcomes explored four major points regarding sorption of SAA and HCB on soil. 1. The interaction depends on the SOM chemical composition more than the SOM content. 2. The interaction causes a site-specific adsorption on the soil surfaces. 3. Sorption hysteresis occurs, which can be explained by inclusion of these pollutants inside soil voids. 4. It was observed that the hydrophobic HCB is adsorbed on soil stronger than the hydrophilic SAA. Moreover, the theoretical results showed that HCB forms stable complexes with all SOM models in the aqueous solution while most of SAA complexes are accompanied by dissociation. The SOM-cavity modeling showed a significant effect on binding of organic pollutants to SOM.

preprint2014arXiv

Laser-Driven Localization of Collective CO Vibrations in Metal-Carbonyl Complexes

Using the example of a cobalt dicarbonyl complex it is shown that two perpendicularly polarized IR laser pulses can be used to trigger an excitation of the delocalized CO stretching modes, which corresponds to an alternating localization of the vibration within one CO bond. The switching time for localization in either of the two bonds is determined by the energy gap between the symmetric and asymmetric fundamental transition frequencies. The phase of the oscillation between the two local bond excitations can be tuned by the relative phase of the two pulses. The extend of control of bond localization is limited by the anharmonicity of the potential energy surfaces leading to wave packet dispersion. This prevents such a simple pulse scheme from being used for laser-driven bond breaking in the considered example.

preprint2013arXiv

A New Efficient Method for Calculation of Frenkel Exciton Parameters in Molecular Aggregates

The Frenkel exciton Hamiltonian is at the heart of many simulations of excitation energy transfer in molecular aggregates. It separates the aggregate into Coulomb-coupled monomers. Here it is shown that the respective parameters, i.e. monomeric excitation energies and Coulomb couplings between transition densities, can be efficiently calculated using time-dependent tight-binding-based density functional theory (TD-DFTB). Specifically, Coulomb couplings are expressed in terms of self-consistently determined Mulliken transition charges. The determination of the sign of the coupling requires an additional super-molecule calculation. The approach is applied to two dimer systems. First, formaldehyde oxime for which a detailed comparison with standard DFT using the B3LYP and the PBE functionals is provided. Second, the Coulomb coupling is explored in dependence on the intermolecular coordinates for a perylene bisimide dimer. This provides structural evidence for the previously observed biphasic aggregation behavior of this dye.

preprint2013arXiv

Chemical Bonding in Aqueous Ferrocyanide: Experimental and Theoretical X-ray Spectroscopic Study

Resonant inelastic X-ray scattering (RIXS) and X-ray absorption (XA) experiments at the iron L- and nitrogen K-edge are combined with high-level first principles restricted active space self-consistent field (RASSCF) calculations for a systematic investigation of the nature of the chemical bond in potassium ferrocyanide in aqueous solution. The atom- and site-specific RIXS excitations allow for direct observation of ligand-to-metal (Fe L-edge) and metal-to-ligand (N K-edge) charge transfer bands and thereby evidence for strong σ-donation and π-back-donation. The effects are identified by comparing experimental and simulated spectra related to both the unoccupied and occupied molecular orbitals in solution.

preprint2013arXiv

How Soil Organic Matter Composition Controls Hexachlorobenzene-Soil-Interactions: Adsorption Isotherms and Quantum Chemical Modelling

Hazardous persistent organic pollutants (POPs) interact in soil with the soil organic matter (SOM) but this interaction is insufficiently understood at the molecular level. We investigated the adsorption of hexachlorobenzene (HCB) on soil samples with systematically modified SOM. These samples included the original soil, the soil modified by adding a hot water extract (HWE) fraction (soil+3 HWE and soil+6 HWE), and the pyrolyzed soil. The SOM contents increased in the order pyrolyzed soil < original soil < soil+3 HWE < soil+6 HWE. For the latter three samples this order was also valid for the HCB adsorption. The pyrolyzed soil adsorbed more HCB than the other samples at low initial concentrations, but at higher concentrations the HCB adsorption became weaker than in the samples with HWE addition. This adsorption behaviour combined with the differences in the chemical composition between the soil samples suggested that alkylated aromatic, phenol, and lignin monomer compounds contributed most to the HCB adsorption. To obtain a molecular level understanding, a test set has been developed on the basis of elemental analysis which comprises 32 representative soil constituents. The calculated binding free energy for HCB with each representative system shows that HCB binds to SOM stronger than to soil minerals. For SOM, HCB binds to alkylated aromatic, phenols, lignin monomers, and hydrophobic aliphatic compounds stronger than to polar aliphatic compounds confirming the above adsorption isotherms. Moreover, quantitative structure-activity relationship (QSAR) of the binding free energy with independent physical properties of the test set systems for the first time indicated that the polarizability, the partial charge on the carbon atoms, and the molar volume are the most important properties controlling HCB-SOM interactions.

preprint2013arXiv

State-Dependent Electron Delocalization Dynamics at the Solute-Solvent Interface: Soft X-ray Absorption Spectroscopy and Ab Initio Calculations

Non-radiative decay channels in the L-edge fluorescence spectra from transition metal-aqueous solutions give rise to spectral dips in X-ray transmission spectra. Their origin is unraveled here using partial and inverse partial fluorescence yields on the micro-jet combined with multi-reference ab initio electronic structure calculations. Comparing Fe2+, Fe3+, and Co2+ systems we demonstrate unequivocally that spectral dips are due to a state-dependent electron delocalization within the manifold of d-orbitals.

preprint2012arXiv

Laser Control of Dissipative Two-Exciton Dynamics in Molecular Aggregates

There are two types of two-photon transitions in molecular aggregates, that is, non-local excitations of two monomers and local double excitations to some higher excited intra-monomer electronic state. As a consequence of the inter-monomer Coulomb interaction these different excitation states are coupled to each other. Higher excited intra-monomer states are rather short-lived due to efficient internal conversion of electronic into vibrational energy. Combining both processes leads to the annihilation of an electronic excitation state, which is a major loss channel for establishing high excitation densities in molecular aggregates. Applying theoretical pulse optimization techniques to a Frenkel exciton model it is shown that the dynamics of two-exciton states in linear aggregates (dimer to tetramer) can be influenced by ultrafast shaped laser pulses. In particular, it is studied to what extent the decay of the two-exciton population by inter-band transitions can be transiently suppressed. Intra-band dynamics is described by a dissipative hierarchy equation approach, which takes into account strong exciton-vibrational coupling in the non-Markovian regime.

preprint2012arXiv

Quantum Dynamics and Spectroscopy of Excitons in Molecular Aggregates

The theoretical description and the properties of Frenkel excitons in non-covalently bonded molecular aggregates are reviewed from a multi-exciton perspective of dissipative quantum dynamics. First, the photophysical and quantum chemical characterization of the monomeric dye building blocks is discussed, including the important aspect of electron-vibrational coupling within the Huang-Rhys model. Supplementing the model by the Coulombic interactions between monomers, the description of aggregates in terms of excitonic or vibrational-excitonic bands follows. Besides of giving rise to complex absorption and emission line shapes, exciton-vibrational interaction is responsible for energy and phase relaxation and thereby limits the size of coherent excitations in larger aggregates. Throughout, emphasis is put on the electronic three-level model as a minimum requirement to describe nonlinear spectroscopies including effects of two-exciton states such as excited state absorption and exciton-exciton annihilation. The experimentally observed characteristics of stationary absorption and fluorescence spectra of aggregates as well as their temperature dependence are discussed. Examples for ultrafast spectroscopic experiments including pump-probe studies, photon echo and two-dimensional spectroscopy are presented and results on the size of coherence domains and on intra- and interband relaxation are given. Finally, experimental signatures for exciton-exciton annihilation and their analysis with respect to the mobility of excitons are described.

preprint2012arXiv

Size dependent exciton dynamics in one-dimensional perylene bisimide aggregates

The size dependent exciton dynamics of one-dimensional aggregates of substituted perylene bisimides are studied by ultrafast transient absorption spectroscopy and kinetic Monte-Carlo simulations in dependence on the temperature and the excitation density. For low temperatures the aggregates can be treated as infinite chains and the dynamics is dominated by diffusion driven exciton-exciton annihilation. With increasing temperature the aggregates decompose into small fragments consisting of very few monomers. This scenario is also supported by the time dependent anisotropy deduced from polarization dependent experiments.

preprint2011arXiv

Carbonyl Vibrational Wave Packet Circulation in Mn$_2$(CO)$_{10}$ Driven by Ultrashort Polarized Laser Pulses

The excitation of the degenerate $E_1$ carbonyl stretching vibrations in dimanganese decacarbonyl is shown to trigger wave packet circulation in the subspace of these two modes. On the time scale of about 5 picoseconds intramolecular anharmonic couplings do not cause appreciable disturbance, even under conditions where the two $E_1$ modes are excited by up to about two vibrational quanta each. The compactness of the circulating wave packet is shown to depend strongly on the excitation conditions such as pulse duration and field strength. Numerical results for the solution of the seven-dimensional vibrational Schrödinger equation are obtained for a density functional theory based potential energy surface and using the multi-configuration time-dependent Hartree method.

preprint2011arXiv

Energy transfer among distant quantum systems in spatially shaped laser fields: Two H atoms with the internuclear separation of 5.29 nm (100 a.u.)

The quantum dynamics of two distant H atoms excited by ultrashort and spatially shaped laser pulses is studied by the numerical solution of the non-Born-Oppenheimer time-dependent Schrödinger equation within a three-dimensional (3D) model, including the internuclear distance R and the two z coordinates of the electrons, z1 and z2. The two 1D hydrogen atoms, A and B, are assumed to be initially in their ground states with a large (but otherwise arbitrary) internuclear separation of R = 100 a.u. (5.29 nm). Two types of a spatial envelope of a laser field linearly po- larized along the z-axis are considered: (i) a broad Gaussian envelope, such that atom A is excited by the laser field predominantly, and (ii) a narrow envelope, such that practically only atom A is excited by the laser field. With the laser carrier frequency ω = 1.0 a.u. and the pulse duration tp = 5 fs, in both cases an efficient energy transfer from atom A to atom B has been found. The ionization of atom B achieved mostly after the end of the laser pulse is close to or even higher than that of atom A. It is shown that with a narrow spatial envelope of the laser field, the underlying mechanisms of the energy transfer from A to B and the ionization of B are the Coulomb attraction of the laser driven electron by the proton of atom B and a short-range Coulomb repulsion of the two electrons when their wave functions significantly overlap in the domain of atom B. In the case of a broad Gaussian spatial envelope of the laser field, the opposite process also occurs, but with smaller probability: the energy is transferred from the weakly excited atom B to atom A, and the ionization of atom A is also induced by the electron-electron repulsion in the domain of atom A due to a strong overlap of the electronic wave functions.

preprint2011arXiv

Path Integral Approach to the Calculation of Reaction Rates for a Reaction Coordinate Coupled to a Dual Harmonic Bath

We present a new method for the numerical calculation of canonical reaction rate constants in complex molecular systems, which is based on a path integral formulation of the flux-flux correlation function. Central is the partitioning of the total system into a relevant part coupled to a dual bath. The latter consists of two parts: First, a set of strongly coupled harmonic modes, describing, for example, intramolecular degrees of freedom. They are treated on the basis of a reaction surface Hamiltonian approach. Second, a set of bath modes mimicking an unspecific environment modeled by means of a continuous spectral density. After deriving a set of general equations expressing the canonical rate constant in terms of appropriate influence functionals, several approximations are introduced to provide an efficient numerical implementation. Results for an initial application to the H-transfer in 6-Aminofulvene-1-aldimine are discussed.

preprint2009arXiv

Infrared Laser Driven Double Proton Transfer. An Optimal Control Theory Study

Laser control of ultrafast double proton transfer is investigated for a two-dimensional model system describing stepwise and concerted transfer pathways. The pulse design has been done by employing optimal control theory in combination with the multiconfiguration time-dependent Hartree wave packet propagation. The obtained laser fields correspond to multiple pump-dump pulse sequences. Special emphasis is paid to the relative importance of stepwise and concerted transfer pathways for the driven wave packet and its dependence on the parameters of the model Hamiltonian as well as on the propagation time. While stepwise transfer is dominating in all cases considered, for high barrier systems concerted transfer proceeding via tunneling can make a contribution.

preprint2009arXiv

IR Spectrum of the O-H$...$O Hydrogen Bond of Phthalic Acid Monomethylester in Gas Phase and in CCl$_4$ Solution

The absorption spectrum of the title compound in the spectral range of the Hydrogen-bonded OH-stretching vibration has been investigated using a five-dimensional gas phase model as well as a QM/MM classical molecular dynamics simulation in solution. The gas phase model predicts a Fermi-resonance between the OH-stretching fundamental and the first OH-bending overtone transition with considerable oscillator strength redistribution. The anharmonic coupling to a low-frequency vibration of the Hydrogen bond leading to a vibrational progression is studied within a diabatic potential energy curve model. The condensed phase simulation of the dipole-dipole correlation function results in a broad band in the 3000 \cm region in good agreement with experimental data. Further, weaker absorption features around 2600 \cm have been identified as being due to motion of the Hydrogen within the Hydrogen bond.