Quantum dynamics simulations using Gaussian wavepackets: the vMCG method

Electron Dynamics upon Ionization of Polyatomic Molecules: Coupling to Quantum Nuclear Motion and Decoherence

Knowledge about the electronic motion in molecules is essential for our understanding of chemical reactions and biological processes. The advent of attosecond techniques opens up the possibility to induce electronic motion, observe it in real time, and potentially steer it. A fundamental question remains the factors influencing electronic decoherence and the role played by nuclear motion in this process. Here, we simulate the dynamics upon ionization of the polyatomic molecules paraxylene and modified bismethylene-adamantane, with a quantum mechanical treatment of both electron and nuclear dynamics using the direct dynamics variational multiconfigurational Gaussian method. Our simulations give new important physical insights about the expected decoherence process. We have shown that the decoherence of electron dynamics happens on the time scale of a few femtoseconds, with the interplay of different mechanisms: the dephasing is responsible for the fast decoherence while the nuclear overlap decay may actually help maintain it and is responsible for small revivals.

Accelerating direct quantum dynamics using graphical processing units

The direct dynamics variational multi-configurational Gaussian (DD-vMCG) method is combined with electronic structure calculations accelerated by Graphical Processing Units (GPUs). This is used to identify GPU acceleration will have a significant effect for both ground and excited state simulations.

Ab initio calculation of inelastic scattering

We calculate nonresonant inelastic electron and X-ray scattering cross sections for bound-to-bound transitions in atoms and molecules from ab initio electronic wavefunctions.

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.

Quantum Statistical Ensemble Resilient to Thermalization

The sampling of the wave function within a suitable ensemble is an important tool in the statistical analysis of a molecule interacting with its environment. The uniform statistical distribution of quantum pure states in an active space is often the privileged choice. However, such a distribution with constant average populations of eigenstates is not preserved upon the interaction between quantum systems. This appears as a severe methodological shortcoming, as long as a quantum system can be always considered as the result of interactions among previously isolated subsystems. In the present work we formulate an alternative statistical ensemble of pure states that is robust with respect to interaction, and it is thus preserved when subsystems are merged. It is derived from the condition of invariance of the average populations upon interaction between quantum systems in the same thermal state. These average populations allow a simple identification of the thermodynamic properties of the system. We find that such a statistical distribution is robust with respect to interaction of systems at different temperatures reproducing the thermalization of macroscopic bodies, and for this reason we identify it as the Thermalization Resilient Ensemble.

Challenges facing an understanding of the nature of low-energy excited states in photosynthesis

While the majority of the photochemical states and pathways related to the biological capture of solar energy are now well understood and provide paradigms for artificial device design, additional low-energy states have been discovered in many systems with obscure origins and significance. However, as low-energy states are naively expected to be critical to function, these observations pose important challenges. A review of known properties of low energy states covering eight photochemical systems, and options for their interpretation, are presented. A concerted experimental and theoretical research strategy is suggested and outlined, this being aimed at providing a fully comprehensive understanding.

Non-adiabatic excited state molecular dynamics of phenylene ethynylene dendrimer using a multiconfigurational Ehrenfest approach

Photoinduced dynamics of electronic and vibrational unidirectional energy transfer between meta-linked building blocks in a phenylene ethynylene dendrimer is simulated using a multiconfigurational Ehrenfest in time-dependent diabatic basis (MCE-TDDB) method.