In this thesis quantum-classical dynamics is applied to the study of quantum condensed phase processes. This approach is based on the quantum-classical Liouville equation where the dynamics of a small subset of the degrees of freedom are treated quantum mechanically while the remaining degrees of freedom are treated by classical mechanics to a good approximation. We use this approach as it is computationally tractable, and the resulting equation of motion accurately accounts for the quantum and classical dynamics, as well as the coupling between these two components of the system. By recasting the quantum-classical Liouville equation into the form of a generalized master equation we investigate connections to surface-hopping. The link between these approaches is decoherence arising from interaction of the subsystem with the environment. We derive an evolution equation for the subsystem which contains terms accounting for the effects of the environment. One of these terms involves a memory kernel that accounts for the coherent dynamics. If this term decays rapidly, a Markovian approximation can be made. By lifting the resulting subsystem master equation into the full phase space, we obtain a Markovian master equation that prescribes surface-hopping-like dynamics. Our analysis outlines the conditions under which such a description is valid. Next, we consider the calculation of the rate constant for a quantum mechanical barrier crossing process. Starting from the reactive-flux autocorrelation function, we derive a quantum-classical expression for the rate kernel. This expression involves quantum-classical evolution of a species operator averaged over the initial quantum equilibrium structure of the system making it possible to compute the rate constant via computer simulation. Using a simple model for a proton transfer reaction we compare the results of the rate calculation obtained by quantum-classical Liouville dynamics with that of master equation dynamics. The master equation provides a good approximation to the full quantum-classical Liouville calculation for our model and a more stable algorithm results due to the elimination of oscillating phase factors in the simulation. Finally, we make use of the theoretical framework established in this thesis to analyze some aspects of decoherence used in popular surface-hopping techniques.

In this thesis quantum-classical dynamics is applied to the study of quantum condensed phase processes. This approach is based on the quantum-classical Liouville equation where the dynamics of a small subset of the degrees of freedom are treated quantum mechanically while the remaining degrees of freedom are treated by classical mechanics to a good approximation. We use this approach as it is computationally tractable, and the resulting equation of motion accurately accounts for the quantum and classical dynamics, as well as the coupling between these two components of the system. By recasting the quantum-classical Liouville equation into the form of a generalized master equation we investigate connections to surface-hopping. The link between these approaches is decoherence arising from interaction of the subsystem with the environment. We derive an evolution equation for the subsystem which contains terms accounting for the effects of the environment. One of these terms involves a memory kernel that accounts for the coherent dynamics. If this term decays rapidly, a Markovian approximation can be made. By lifting the resulting subsystem master equation into the full phase space, we obtain a Markovian master equation that prescribes surface-hopping-like dynamics. Our analysis outlines the conditions under which such a description is valid. Next, we consider the calculation of the rate constant for a quantum mechanical barrier crossing process. Starting from the reactive-flux autocorrelation function, we derive a quantum-classical expression for the rate kernel. This expression involves quantum-classical evolution of a species operator averaged over the initial quantum equilibrium structure of the system making it possible to compute the rate constant via computer simulation. Using a simple model for a proton transfer reaction we compare the results of the rate calculation obtained by quantum-classical Liouville dynamics with that of master equation dynamics. The master equation provides a good approximation to the full quantum-classical Liouville calculation for our model and a more stable algorithm results due to the elimination of oscillating phase factors in the simulation. Finally, we make use of the theoretical framework established in this thesis to analyze some aspects of decoherence used in popular surface-hopping techniques.

No description available.

Light or photons can excite electrons in a molecule, leading to creation of electronically excited states. Such processes are ubiquitous in nature, such as, vision, photo-protection of DNA/RNA nucleobases, light harvesting, energy and charger transfer etc. This photoexcitation induces nuclear motion on the excited states, leading the excess energy to dissipate either non-radiatively via internal conversion back down to the ground state, isomerization, and dissociation, or radiatively via fluorescence and phosphorescence. In this dissertation, we investigate the non-radiative processes in organic chromophores that ensue in an ultrafast manner, mediated via conical intersections (CoIn). Description of such excited state processes generally require multi-reference treatment because of quasi-degeneracy near CoIns. Hence, most insight about these processes is typically gained by constructing potential energy surface (PES) using multi-reference electronic structure methods along important reaction coordinates. Nonetheless, the aforementioned static treatment fails to provide any dynamical information, such as, excited state lifetime, state populations, branching ratio, quantum yield etc. In this dissertation, we have gone beyond the static treatment by undertaking computationally expensive non-adiabatic excited state molecular dynamics simulations employing trajectory surface hopping (TSH) methodology on PESs created on-the-fly using multi-reference electronic structure methods. This allows us to compare theoretical results to experimental observables, when possible, strengthening the explanations underlying those processes. Our goal is to examine the effect of structure, and of electronic structure methods on the excited state dynamics. We have examined the non-adiabatic excited state dynamics of cis,cis-1,3-cyclooctadiene (cc-COD), a cyclic diene, in an effort to systematically compare and contrast the dynamics of cc-COD to that of other well studied conjugated molecules. Such exploration is very significant, since the majority of the molecules involved in natural photoexcited processes, include an ethylenic double bond or alternating double bonds creating conjugation. Our calculations have revealed ultrafast sub-ps decay for cc-COD, and have illustrated that the internal conversion dynamics is facilitated by CoIns, dominated by twisting of one of the double bonds and pyramidalization of one of the carbons of that double bond, similar to trans-1,3-butadiene and unlike 1,3-cyclohexadiene (CHD). Our high-level electronic structure calculations have also explained the features in the experimental time-resolved photoelectron spectrum of cc-COD. Another molecule of biological importance, uracil, was also investigated using TSH simulations, by systematically increasing dynamical correlation. We have found that the inclusion of dynamical correlation for uracil leads to an almost barrierless PES on S2, leading to a faster decay and no population trap on this state. Uracil also contains a double bond and the simulations have revealed that the ultrafast relaxation is dominated by an ethylenic twist and pyramidalization of a carbon of that bond, increasing importance of such nuclear motion in photoexcited molecular dynamics. A comparison of the molecules studied have illustrated that the rigid molecules, such as uracil, CHD, have a very local CoIn seam space, whereas cc-COD, which is flexible having many low frequency degrees of freedom, has a non-local or extended CoIn seam space. Overall, the work performed in this dissertation, elucidates the significance of structure and conjugation, in the photoinduced coupled electron-nuclear dynamics in organic molecules. / Chemistry

We utilize the exact factorization of the electron-nuclear wave function [Abedi et al., PRL 105 123002 (2010)] to illuminate several aspects of laser-driven molecular dynamics in intense femtosecond laser pulses. Above factorization allows for a splitting of the full molecular wave function and leads to a time-dependent Schrödinger equation for the nuclear subsystem alone which is exact in the sense that the absolute square of the corresponding, purely nuclear, wave function yields the exact nuclear N-body density of the full electron-nuclear system. As one remarkable feature, this factorization provides the exact classical force, the force which contains the highest amount of electron-nuclear correlations that can be retained in the quantum-classical limit of the electron-nuclear system. We re-evaluate the classical limit of the nuclear Schrödinger equation from the perspective of the exact factorization, and address the long-standing question of the validity of the popular quantum-classical surface hopping approach in laserdriven cases. In particular, our access to the exact classical force allows for an elaborate evaluation of the various and completely different potential energy surfaces frequently applied in surface hopping calculations. The highlight of this work consists in a generalization of the exact factorization and its application to the laser-driven molecular wave function in the Floquet picture, where the molecule and the laser form an united quantum system exhibiting its own Hilbert space. This particular factorization enables us to establish an analytic connection between the exact nuclear force and Floquet potential energy surfaces. Complementing above topics, we combine different well-known and proven methods to give a systematic study of molecular dissociation mechanisms for the complicated electric fields provided by modern attosecond laser technology.:Contents Introduction 1 The exact factorization of time-dependent wave functions 1.1 Concern and state of the art 1.2 The exact factorization of the electron-nuclear wave function 1.3 The generalized exact factorization 1.4 The exact factorization for coupled harmonic oscillators 1.5 The exact factorization for a single particle with spin 1.6 The exact factorization of the laser-driven electron-nuclear wave function in the Floquet picture 1.7 Summary and conclusion 2 Quantum-classical molecular dynamics from an exact factorization perspective 2.1 Concern and state of the art 2.2 The exact nuclear TDSE 2.3 The Wigner-Moyal equation for the nuclear TDSE and its classical limit 2.4 The Bohmian formulation of the nuclear TDSE and its classical limit 2.5 Comparative calculations 2.5.1 Scenario 1: stationary states 2.5.2 Scenario 2: laser-driven dynamics 2.6 Summary and conclusion 3 Surface hopping in laser-driven molecular dynamics 3.1 Concern and state of the art 3.2 Surface hopping 3.3 Quantum-classical dynamics on the EPES 3.4 The benchmark model and its potential energy surfaces 3.5 Surface hopping in laser-driven molecular dynamics 3.6 Summary and conclusion 4 Beyond the limit of the Floquet picture: molecular dissociation in few-cycle laser pulses 4.1 Concern and state of the art 4.2 Theoretical few-cycle pulses 4.3 Calculation of dissociation probabilities 4.4 Dissociation in few-cycle pulses 4.4.1 Dissociation in half-cycle pulses 4.4.2 Dissociation in few-cycle pulses 4.5 Dissociation in realistic attosecond pulses 4.6 Summary and conclusion Outlook Appendices A List of abbreviations B Numerical details C Calculating electronic observables within quantum-classical molecular dynamics D Ionization in few-cycle pulses E Modeling an optical attosecond pulse Bibliography

The broad range of applications for photochemical reactions is the result of light-matter interaction at the electronic level. The diverse application of photochemistry in various fields, including photovoltaic materials, molecular switches, and biological systems are due to electronic and structural transformations induced by photoexcitation as well as molecular alteration due to electron and charge transfer. An improved understanding of these photochemical events is dependent on the fundamental theoretical evaluation, to model and analyze the ultrafast processes. The studies discussed in this thesis explore such theoretical implementation in two different frontiers. In the first study, dynamic simulations are performed to model the light-induced bond dissociation of phenyl azide. The surface hopping formalism, implemented under the semiclassical molecular dynamics approach helped in tracing the time evolution of the electronic and structural levels, involved in the photodissociation. In the second study, the time-dependent density functional theory has been applied to generate XA spectra of imidazole solutions. The theoretical assessments support experimental measurements and provide more insight into the core excitations and structural influence on the absorption spectra. / Fotokemiska reaktioner styrs av växelverkan mellan ljus of materia på en elektronisk nivå. I olika fält finns det vitt skilda tillämpningarna av fotokemi. Dessa inkluderar ljusinducerade processer i solceller, molekylära strömbrytare, och biologiska system. Reaktionerna beror av elektroniska och strukturella transformationer som induceras av fotoexcitationen, och kan ge upphov till energi- och laddningsöverföring. För att få utökad förståelse av fotokemiska reaktioner behövs grundläggande teoretiska studier där ultrasnabba processer modelleras och analyseras i jämförelse med experiment. Undersökningar som presenteras i denna avhandling använder sig att teoretiska modeller i två olika områden av fotokemi. I den första studien har vi genomfört dynamiska simuleringar för att modellera ljusinducerad dissociation av fenylazid. Vi have använt en semi-klassisk approximation med hopp mellan olika elektroniska tillstånd vilket gör det möjligt att följa utvecklingen i elektroniska och geometriska frihetsgrader under fotodissociationen. I den andra studien har tidsberoende täthetsfunktionalteori använts för att simulera röntgenabsorptionsspektrum för imidazol i lösning. Genom att utvärdera olika geometriska modeller med teoretiska beräkningar kan vi berika tolkningen av experimentella mätningar, och även få detaljerat insikt i innerskalsexcitationer och hur geometrin påverka röntgenspektrum.

No description available.