Search

Controlling molecular reactivity by shaped laser pulses is a long-standing goal in chemistry. Here, we suggest a direct optimal control approach that combines external pulse optimization with other control parameters arising in the upcoming field of vibro-polaritonic chemistry for enhanced controllability. The direct optimal control approach is characterized by a simultaneous simulation and optimization paradigm, meaning that the equations of motion are discretized and converted into a set of holonomic constraints for a nonlinear optimization problem given by the control functional. Compared with indirect optimal control, this procedure offers great flexibility, such as final time or Hamiltonian parameter optimization. A simultaneous direct optimal control theory will be applied to a model system describing H-atom transfer in a lossy Fabry–Pérot cavity under vibrational strong coupling conditions. Specifically, optimization of the cavity coupling strength and, thus, of the control landscape will be demonstrated. [ABSTRACT FROM AUTHOR]

The gas phase vibrational spectroscopy of the protonated ammonia dimer N2H7+, a prototypical system for strong hydrogen bonding, is studied in the spectral region from 330 to 1650 cm-1 by combining infrared multiple photon dissociation and multidimensional quantum mechanical simulations. The fundamental transition of the antisymmetric proton stretching vibration is observed at 374 cm-1 and assigned on the basis of a six-dimensional model Hamiltonian, which predicts this transition at 471 cm-1. Photodissociation spectra of the larger protonated ammonia clusters NH4+(NH3)n with n=2–4 are also reported for the range from 1050 to 1575 cm-1. The main absorption features can be assigned within the harmonic approximation, supporting earlier evidence that hydrogen bonding in these clusters is considerably weaker than for n=1. [ABSTRACT FROM AUTHOR]

The gas phase vibrational spectroscopy of the protonated ammonia dimer N2H7+, a prototypical system for strong hydrogen bonding, is studied in the spectral region from 330 to 1650 cm-1 by combining infrared multiple photon dissociation and multidimensional quantum mechanical simulations. The fundamental transition of the antisymmetric proton stretching vibration is observed at 374 cm-1 and assigned on the basis of a six-dimensional model Hamiltonian, which predicts this transition at 471 cm-1. Photodissociation spectra of the larger protonated ammonia clusters NH4+(NH3)n with n=2–4 are also reported for the range from 1050 to 1575 cm-1. The main absorption features can be assigned within the harmonic approximation, supporting earlier evidence that hydrogen bonding in these clusters is considerably weaker than for n=1. [ABSTRACT FROM AUTHOR]

There has been some controversy concerning the assignment of measured tunneling splittings for the formic acid dimer in the vibrational ground state and the asymmetric CO-stretching excited state. The discussion is intimately related to the question whether the fundamental excitation of the CO-vibration promotes or hinders tunneling. Here we will address this issue on the basis of a five-dimensional reaction space Hamiltonian which includes three large amplitude coordinates as well as two harmonic modes whose linear superposition reproduces the asymmetric CO-vibrational mode. Within density functional theory using the B3LYP functional together with a 6-311++G(3df,3pd) basis set we obtain a ground state tunneling splitting which is about 2.4 larger than the one for the CO-stretching excited state. [ABSTRACT FROM AUTHOR]

Femtosecond dynamics of molecular vibrations as well as cage motions in the B←X transition of Cl2 in solid Ar have been investigated. We observed molecular vibrational wave-packet motion in experimental pump-probe spectra and an additional oscillation with a 500 fs period which is assigned to the zone-boundary phonon of the Ar crystal. The cage motion is impulsively driven by the B←X transition due to the expansion of the electronic cloud of the chromophore. To clarify the underlying mechanism, we performed simulations based on the diatomics-in-molecules method which takes into account the different shapes of the Cl2 electronic wave function in the B and X states as well as the anisotropic interaction with the matrix. The simulation results show that Ar atom motion in the (100) plane is initiated by the electronic transition and that only those Ar atoms oscillate coherently with an ∼500 fs period which are essentially decoupled from the molecular vibration. Their phase and time evolution are in good agreement with the experimentally observed oscillation, supporting the assignment as a displacive excitation of coherent phonons. [ABSTRACT FROM AUTHOR]

The previously developed instanton theory [G. V. Mil’nikov and H. Nakamura, J. Chem. Phys. 122, 124311 (2005)] is applied to the calculation of vibrationally assisted tunneling splitting of the deuterated formic acid dimer (DCOOH)2 with all the degrees of freedom taken into account. The ground-state tunnel splitting is determined by the density-functional theory combined with coupled cluster level of quantum chemistry to be 0.0038 cm-1 which is comparable to the experimental value of 0.0029 cm-1. Further, the tunnel splittings of fundamental excitations are estimated for frequencies below 300 cm-1. In this energy range it is found that the excitation modes may either enhance or suppress tunneling as compared to the ground state. For the higher-frequency modes a rapid growth of the tunnel splitting is observed. At frequencies above 1000 cm-1 the semiclassical solution becomes unstable and no reliable tunneling splittings can be obtained. This is in vast contrast to the adiabatic approximation to the instanton theory in which the tunnel splittings can be retrieved up to 3000 cm-1. We discuss this disparity from the viewpoint of the multidimensional character of tunneling in hydrogen bonds and the adiabatic approximation is concluded to be inaccurate. [ABSTRACT FROM AUTHOR]

A dissipative bistable system presents the simplest model to describe condensed phase reaction dynamics. Using a quantum master equation approach to calculate multitime dipole correlation functions we demonstrate how the dissipative dynamics can be characterized by time-resolved third-order infrared spectroscopy. Thereby we incorporate bilinear and linear-quadratic systembath interaction into the Redfield relaxation tensor. Investigating equilibrium and nonequilibrium initial conditions for a symmetric system it is shown that bath-induced coherence transfer can have a dramatic influence on the two-dimensional signals. This occurs when the inverse of the ground state tunneling splitting is of the order of the coherence transfer time. [ABSTRACT FROM AUTHOR]

The linear and third-order polarizabilities of the light-harvesting antennae of photosynthetic bacteria and green plants are calculated using an equation of motion approach which maps the system onto a coupled set of anharmonic excitonic oscillators. The oscillator representation is shown to have several advantages over the traditional picture based on properties of individual global eigenstates of the aggregate; besides a considerable reduction of computational effort, the dynamics of excitations in the two-exciton band is conveniently analyzed in terms of single-exciton Green’s functions and the two-exciton scattering matrix. © 1996 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The density matrix theory is utilized for the description of ultra fast optical properties and related vibrational wave packet dynamics of molecular systems in condensed media. As an example, optically induced vibrational wave packets in the so-called curve–crossing system are considered. Such a system goes beyond the standard treatment of optical phenomena since the vibrational wave packet moves in a double well potential and is subject to environmental influences like wave function dephasing and relaxation. The complete theoretical description has been carried out in a representation of the vibrational wave functions of the diabatic states which refer to the two coupled vibrational surfaces. Solving the corresponding density matrix equations by numerical methods allows us to incorporate the static coupling between the crossed surfaces in a nonperturbative manner. Standard projection operator technique is used to treat environmental contributions up to the second order. For the case of a bilinear coupling between the molecular system and the environment we determine the time development of an initially prepared Gaussian wave packet. Corresponding time-resolved spectra of a pump–probe configuration are also derived. The developed formalism is finally applied to the system of the coupled electronic states c’ 1Σ+u and b’ 1Σ+u of the N2 molecule. © 1994 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The dynamics of cyclic systems with four equivalent potential minima is studied here from two different points of view. The solution of the time-dependent Schrödinger equation provides insight into the coherent wave packet motion. The resulting reaction mechanism involves relocalization between opposite, not neighboring potential minima. The inclusion of an environment within a density matrix description leads to dissipation and therefore to a transition from coherent to incoherent dynamics. The theoretical considerations are applied to a simple model of the cyclic motion of a proton in a molecular framework. [ABSTRACT FROM AUTHOR]

Over the last several decades, the light-harvesting protein complexes of purple bacteria have been among the most popular model systems for energy transport in excitonic systems in the weak and intermediate intermolecular coupling regime. Despite this extensive body of scientific work, significant questions regarding the excitonic states and the photo-induced dynamics remain. Here, we address the low-temperature electronic structure and excitation dynamics in the light-harvesting complex 2 of Rhodopseudomonas acidophila by two-dimensional electronic spectroscopy. We find that, although at cryogenic temperature energy relaxation is very rapid, exciton mobility is limited over a significant range of excitation energies. This points to the presence of a sub-200 fs, spatially local energy-relaxation mechanism and suggests that local trapping might contribute substantially more in cryogenic experiments than under physiological conditions where the thermal energy is comparable to or larger than the static disorder.

Frenkel exciton population dynamics of an excitonic dimer is studied by comparing the results from a quantum master equation involving rates from second-order perturbative treatment with respect to the excitonic coupling with the 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 the cumulant expansion. The coupling of electronic transitions to bath modes is modeled either as overdamped oscillators for the description of thermal bath components or as underdamped oscillators to account for intramolecular vibrations. Cases of initial nonequilibrium and equilibrium vibrations are discussed. In the 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, this approach can also 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.

The Caldeira-Leggett model provides a compact characterization of a thermal environment in terms of a spectral density function, which has led to a variety of numerically exact quantum methods for reduced density matrix propagation. Since spectral densities are often computed from classical molecular dynamics simulations, we investigate in this paper whether quantum effects should be accounted for in the calculations. Therefore, we reformulate the recently developed Fourier method for spectral density calculations from semiclassical simulations which approximately allow for quantum effects. We propose two possible protocols based on either correlation functions or expectation values. These protocols are tested on a generic Calderra-Leggett model for the linearized semiclassical initial-value representation (LSC-IVR), the thawed Gaussian wave packet dynamics (TGWD), and hybrid schemes combining the two with the more accurate Herman-Kluk formula. Surprisingly, spectral densities from the LSC-IVR method, which treats the dynamics completely classically, are found to be extremely accurate, even in the quantum regime, where this method does not give a correct description of the correlation functions and expectation values. In contrast, the TGWD method turns out as too inaccurate for spectral density calculations, and the hybrid schemes perform well only if the system is close to the classical regime. This implies that, if the bath has a Caldeira-Leggett form, spectral densities are insensitive to quantum effects and any effort to approximately account for them rather leads to errors. Hence, in this case, spectral densities can be computed from classical simulations and used in a reduced quantum simulation as well.

Measuring vibronic spectra probes dynamical processes in molecular systems. When interpreted via suitable theoretical tools, the experimental data provides comprehensive information about the system in question. For complex many-body problems, such an approach usually requires the formulation of proper classical-like approximations, which is particularly challenging if multiple electronic states are involved. In this work, we express the imaginary-time shifted time correlation function and, thus, the vibronic spectrum in terms of the so-called Matsubara dynamics, which combines quantum statistics and classical-like dynamics. By applying the Matsubara approximation in the adiabatic limit, we derive a formal generalization of the existing Matsubara dynamics formalism to multiple potential energy surfaces (PESs), which, however, does not feature all the defining properties of its single-PES counterpart though suffering equally from the sign problem. The mathematical analysis for two shifted harmonic oscillators suggests a new modified method to practically simulate the standard correlation function via Matsubara-like dynamics. Importantly, this modified method samples the thermal Wigner function without suffering from the sign problem and yields an accurate approximation to the vibronic absorption spectrum, not only for the harmonic system but also for the anharmonic one.

Förster theory describes electronic exciton energy migration in molecular assemblies as an incoherent hopping process between donor and acceptor molecules. The rate is expressed in terms of the overlap integral between donor fluorescence and acceptor absorption spectra. Typical time scales for systems like photosynthetic antennae are on the order of a few picoseconds. Prior to transfer, it is assumed that the initially excited donor molecule has equilibrated with respect to the local environment. However, upon excitation and during the equilibration phase, the state of the system needs to be described by the full density matrix, including coherences between donor and acceptor states. While being intuitively clear, addressing this regime experimentally has been a challenge until the recently reported advances in fluorescence detected two-dimensional spectroscopy. Here, we demonstrate using fourth order perturbation theory the conditions for the presence of donor-acceptor coherence induced cross-peaks at zero waiting time between the first and the second pair of pulses. The approach is illustrated for a heterodimer model which facilitates an analytical solution.

The framework to approach quasi-classical dynamics in the electronic ground state is well established and is based on the Kubo-transformed time correlation function (TCF), being the most classical-like quantum TCF. Here we discuss whether the choice of the Kubo-transformed TCF as a starting point for simulating vibronic spectra is as unambiguous as it is for vibrational ones. Employing imaginary-time path integral techniques in combination with the interaction representation allowed us to formulate a method for simulating vibronic spectra in the adiabatic regime that takes nuclear quantum effects and dynamics on multiple potential energy surfaces into account. Further, a generalized quantum TCF is proposed that contains many well-established TCFs, including the Kubo one, as particular cases. Importantly, it also provides a framework to construct new quantum TCFs. Applying the developed methodology to the generalized TCF leads to a plethora of simulation protocols, which are based on the well-known TCFs as well as on new ones. Their performance is investigated on 1D anharmonic model systems at finite temperatures. It is shown that the protocols based on the new TCFs may lead to superior results with respect to those based on the common ones. The strategies to find the optimal approach are discussed.

Modern X-ray spectroscopy has proven itself as a robust tool for probing the electronic structure of atoms in complex environments. Despite working on energy scales that are much larger than those corresponding to nuclear motions, taking nuclear dynamics and the associated nuclear correlations into account may be of importance for X-ray spectroscopy. Recently, we have developed an efficient protocol to account for nuclear dynamics in X-ray absorption and resonant inelastic X-ray scattering spectra [Karsten et al., J. Phys. Chem. Lett. 8, 992 (2017)], based on ground state molecular dynamics accompanied with state-of-the-art calculations of electronic excitation energies and transition dipoles. Here, we present an alternative derivation of the formalism and elaborate on the developed simulation protocol using gas phase and bulk water as examples. The specific spectroscopic features stemming from the nuclear motions are analyzed and traced down to the dynamics of electronic energy gaps and transition dipole correlation functions. The observed tendencies are explained on the basis of a simple harmonic model, and the involved approximations are discussed. The method represents a step forward over the conventional approaches that treat the system in full complexity and provides a reasonable starting point for further improvements.

Effects of thermal fluctuations on the electronic excitation energies and intermonomeric Coulomb couplings are investigated for a perylene-tetracarboxylic-diimide crystal. To this end, time dependent density functional theory based tight binding (TD-DFTB) in the linear response formulation is used in combination with electronic ground state classical molecular dynamics. As a result, a parametrized Frenkel exciton Hamiltonian is obtained, with the effect of exciton-vibrational coupling being described by spectral densities. Employing dynamically defined normal modes, these spectral densities are analyzed in great detail, thus providing insight into the effect of specific intramolecular motions on excitation energies and Coulomb couplings. This distinguishes the present method from approaches using fixed transition densities. The efficiency by which intramolecular contributions to the spectral density can be calculated is a clear advantage of this method as compared with standard TD-DFT.

The applicability of the density functional based tight binding (DFTB) method to the description of hydrogen bond dynamics and infrared (IR) spectroscopy is addressed for the exemplary protic ionic liquid triethylammonium nitrate. Potential energy curves for proton transfer in gas and liquid phases are shown to be comparable to the high level coupled cluster theory in the thermally accessible range of bond lengths. Geometric correlations in the hydrogen bond dynamics are analyzed for a cluster of six ion pairs. Comparing DFTB and DFT data lends further support for the reliability of the DFTB method. Therefore, DFTB bulk simulations are performed to quantify the extent of geometric correlations in terms of Pauling's bond order model. Further, IR absorption spectra are obtained using DFTB and analyzed putting emphasis on the signatures of hydrogen bonding in the NH-stretching and far IR hydrogen bond range.

The coupled quantum dynamics of excitonic and vibrational degrees of freedom is investigated for high-dimensional models of the Fenna-Matthews-Olson 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 excitations. 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.

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 [F. Gottwald et al., J. Phys. Chem. Lett. 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 [F. Gottwald et al., J. Chem. Phys. 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 practical 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.

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-level 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 [Fe(H2O)6](2+) complex. The results show good agreement with the experimental data obtained in this work, whereas the sudden approximation demonstrates distinct deviations from experiments.

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 a 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, which can be rigorously derived by means of a linear projection 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 more 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 overestimation of relaxation times, unless the system under study consists of a harmonic bath bi-linearly coupled to the relevant degrees of freedom.

Using the example of a cobalt dicarbonyl complex it is shown that two perpendicular linearly 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.

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 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 as well as with SCS-CC2 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.

The vertical singlet-singlet and singlet-triplet electronic excitation energies of bis(2-phenylpyridinato-)(2,2'-bipyridine)iridium(III) ([Ir(ppy)(2)(bpy)](+)) are calculated on the basis of a comparative quantum chemical study using wave function methods such as CASSCF∕CASPT2 and density functional theory (TDDFT) with local and range-separated functionals. The TDDFT results show a strong dependence of the charge-transfer transition energies on the amount of the exact exchange in the functional. In general, TDDFT with range-separated functionals provides a good agreement with the experimental spectra. As a result a new assignment of the absorption spectrum of the title compound is proposed.

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 ps, 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., (© 2011 American Institute of Physics)

In this work we present an all-Cartesian reaction surface approach, where the large amplitude coordinates span the so-called reaction plane, that is, the unique plane defined by the two minima and the saddle-point structure of an isomerization reaction. Orthogonal modes are treated within harmonic approximation which gives the total Hamiltonian an almost separable form that is suitable for multidimensional quantum dynamics calculations. The reaction plane Hamiltonian is constructed for the H-atom transfer in tropolone as an example for a system with an intramolecular O...H-O hydrogen bond. We find ground-state tunneling splittings of 3.5 and 0.16 cm(-1) for the normal and deuterated species, respectively. We calculated infrared-absorption spectra for a four-dimensional model focusing on the low-frequency region. Here, we identify a reaction mode which is closely connected to the tautomerization that is reflected in the increase of tunneling splitting to 18 cm(-1) upon excitation.

Based on the Cartesian Reaction Surface framework we construct a four-dimensional potential for the tropolone derivative 3,7-dichlorotropolone, a molecule with an intramolecular O-H...O hydrogen bond. The reduced configuration space involves the in-plane hydrogen atom coordinates, a symmetric O-O vibrational mode, and an antisymmetric mode related to deformations of the seven-membered ring. The system is characterized in terms of quantum mechanical computations of the low-lying eigenstates as well as a classical and semiclassical analysis of spectra obtained via Fourier transforming autocorrelation functions. For the semiclassical analysis we utilize the amplitude-free correlation function method [K. Hotta and K. Takatsuka, J. Phys. A 36, 4785 (2003)]. Our results demonstrate substantial anharmonic couplings leading to highly correlated wave functions even at moderate energies. Furthermore, the importance of dynamical tunneling in tropolone is suggested since many low-lying states--including the ground state--lie above the classical saddle point but nevertheless appear as split pairs.