Jacek Jakowski

Hydrogen Tunneling in an Enzyme Active Site: A Quantum Wavepacket Dynamical Perspective

We study the hydrogen tunneling problem in a model system that represents the active site of the biological enzyme, soybean lipoxygenase-1. Toward this, we utilize quantum wavepacket dynamics performed on potential surfaces obtained by using hybrid density functional theory under the influence of a dynamical active site. The kinetic isotope effect is computed by using the transmission amplitude of the wavepacket, and the experimental value is reproduced. By computing the hydrogen nuclear orbitals (eigenstates) along the reaction coordinate, we note that tunneling for both hydrogen and deuterium occurs through the existence of distorted, spherical s-type proton wave functions and p-type polarized proton wave functions for transfer along the donor-acceptor axis. In addition, there is also a significant population transfer through distorted p-type proton wave functions directed perpendicular to the donor-acceptor axis (via intervening pi-type proton eigenstate interactions) which underlines the three-dimensional nature of the tunneling process. The quantum dynamical evolution indicates a significant contribution from tunneling processes both along the donor-acceptor axis and along directions perpendicular to the donor-acceptor axis. Furthermore, the tunneling process is facilitated by the occurrence of curve crossings and avoided crossings along the proton eigenstate adiabats.

The isotopic effects of deuteration on optoelectronic properties of conducting polymers

The attractive optoelectronic properties of conducting polymers depend sensitively upon intra- and inter-polymer chain interactions, and therefore new methods to manipulate these interactions are continually being pursued. Here, we report a study of the isotopic effects of deuterium substitution on the structure, morphology and optoelectronic properties of regioregular poly(3-hexylthiophene)s with an approach that combines the synthesis of deuterated materials, optoelectronic properties measurements, theoretical simulation and neutron scattering. Selective substitutions of deuterium on the backbone or side-chains of poly(3-hexylthiophene)s result in distinct optoelectronic responses in poly(3-hexylthiophene)/[6,6]-phenyl-C61-butyric acid methyl ester (PCBM) photovoltaics. Specifically, the weak non-covalent intermolecular interactions induced by the main-chain deuteration are shown to change the film crystallinity and morphology of the active layer, consequently reducing the short-circuit current. However, side-chain deuteration does not significantly modify the film morphology but causes a decreased electronic coupling, the formation of a charge transfer state, and increased electron-phonon coupling, leading to a remarkable reduction in the open circuit voltage.

Quantum wave packet ab initio molecular dynamics: an approach to study quantum dynamics in large systems.

A methodology to efficiently conduct simultaneous dynamics of electrons and nuclei is presented. The approach involves quantum wave packet dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator, in conjunction with ab initio molecular dynamics treatment of the electronic and classical nuclear degree of freedom. The latter may be achieved either by using atom-centered density-matrix propagation or by using Born-Oppenheimer dynamics. The two components of the methodology, namely, quantum dynamics and ab initio molecular dynamics, are harnessed together using a time-dependent self-consistent field-like coupling procedure. The quantum wave packet dynamics is made computationally robust by using adaptive grids to achieve optimized sampling. One notable feature of the approach is that important quantum dynamical effects including zero-point effects, tunneling, as well as over-barrier reflections are treated accurately. The electronic degrees of freedom are simultaneously handled at accurate levels of density functional theory, including hybrid or gradient corrected approximations. Benchmark calculations are provided for proton transfer systems and the dynamics results are compared with exact calculations to determine the accuracy of the approach.

Liouville–von Neumann molecular dynamics

We present a novel first principles molecular dynamics scheme, called Liouville-von Neumann molecular dynamics, based on Liouville-von Neumann equation for density matrices propagation and Magnus expansion of the time-evolution operator. The scheme combines formally accurate quantum propagation of electrons represented via density matrices and a classical propagation of nuclei. The method requires a few iterations per each time step where the Fock operator is formed and von Neumann equation is integrated. The algorithm (a) is free of constraint and fictitious parameters, (b) avoids diagonalization of the Fock operator, and (c) can be used in the case of fractional occupation as in metallic systems. The algorithm is very stable, and has a very good conservation of energy even in cases when a good quality conventional Born-Oppenheimer molecular dynamics trajectories is difficult to obtain. Test simulations include initial phase of fullerene formation from gaseous C(2) and retinal system.

The divide–expand–consolidate MP2 scheme goes massively parallel

For large molecular systems conventional implementations of second order Møller–Plesset (MP2) theory encounter a scaling wall, both memory- and time-wise. We describe how this scaling wall can be removed. We present a massively parallel algorithm for calculating MP2 energies and densities using the divide–expand–consolidate scheme where a calculation on a large system is divided into many small fragment calculations employing local orbital spaces. The resulting algorithm is linear-scaling with system size, exhibits near perfect parallel scalability, removes memory bottlenecks and does not involve any I/O. The algorithm employs three levels of parallelisation combined via a dynamic job distribution scheme. Results on two molecular systems containing 528 and 1056 atoms (4278 and 8556 basis functions) using 47,120 and 94,240 cores are presented. The results demonstrate the scalability of the algorithm both with respect to the number of cores and with respect to system size. The presented algorithm is thus highly suited for large super computer architectures and allows MP2 calculations on large molecular systems to be carried out within a few hours – for example, the correlated calculation on the molecular system containing 1056 atoms took 2.37 hours using 94240 cores.

Ab initio based study of the ArO− photoelectron spectra: Selectivity of spin–orbit transitions

A combined ab initio atoms-in-molecule approach was implemented to model the photoelectron spectra of the ArO− anion. The lowest adiabatic states of Σ and Π symmetry of ArO and ArO− were investigated using the fourth-order Moller–Plessett perturbation theory including bond functions. The total energies were dissected into electrostatic, exchange, induction, and dispersion components. The complex of Ar with atomic oxygen is only weakly bound, primarily by dispersion interaction. The Π state possesses a deeper minimum (Re=3.4 A,De=380 μEh) than the Σ state (Re=3.8 A,De=220 μEh). In contrast, the anion complex is fairly strongly bound, primarily by ion-induced dipole induction forces, and the Σ state possesses a deeper minimum at shorter interatomic distances (Re=3.02 A,De=3600 μEh) than the Π state (Re=3.35 A,De=2400 μEh). The Σ–Π splittings in both systems are mainly due to differences in the exchange repulsion terms. Atoms-in-molecule models were used to account for the spin–orbit interaction, and to gene...

Computational Improvements to Quantum Wave Packet ab Initio Molecular Dynamics Using a Potential-Adapted, Time-Dependent Deterministic Sampling Technique

In a recent publication, we introduced a computational approach to treat the simultaneous dynamics of electrons and nuclei. The method is based on a synergy between quantum wave packet dynamics and ab initio molecular dynamics. Atom-centered density-matrix propagation or Born-Oppenheimer dynamics can be used to perform ab initio dynamics. In this paper, wave packet dynamics is conducted using a three-dimensional direct product implementation of the distributed approximating functional free-propagator. A fundamental computational difficulty in this approach is that the interaction potential between the two components of the methodology needs to be calculated frequently. Here, we overcome this problem through the use of a time-dependent deterministic sampling measure that predicts, at every step of the dynamics, regions of the potential which are important. The algorithm, when combined with an on-the-fly interpolation scheme, allows us to determine the quantum dynamical interaction potential and gradients at every dynamics step in an extremely efficient manner. Numerical demonstrations of our sampling algorithm are provided through several examples arranged in a cascading level of complexity. Starting from a simple one-dimensional quantum dynamical treatment of the shared proton in [Cl-H-Cl](-) and [CH3-H-Cl](-) along with simultaneous dynamical treatment of the electrons and classical nuclei, through a complete three-dimensional treatment of the shared proton in [Cl-H-Cl](-) as well as treatment of a hydrogen atom undergoing donor-acceptor transitions in the biological enzyme, soybean lipoxygenase-1 (SLO-1), we benchmark the algorithm thoroughly. Apart from computing various error estimates, we also compare vibrational density of states, inclusive of full quantum effects from the shared proton, using a novel unified velocity-velocity, flux-flux autocorrelation function. In all cases, the potential-adapted, time-dependent sampling procedure is seen to improve the computational scheme tremendously (by orders of magnitude) with minimal loss of accuracy.

Quantum Trajectory-Electronic Structure Approach for Exploring Nuclear Effects in the Dynamics of Nanomaterials.

A massively parallel, direct quantum molecular dynamics method is described. The method combines a quantum trajectory (QT) representation of the nuclear wave function discretized into an ensemble of trajectories with an electronic structure (ES) description of electrons, namely using the density functional tight binding (DFTB) theory. Quantum nuclear effects are included into the dynamics of the nuclei via quantum corrections to the classical forces. To reduce computational cost and increase numerical accuracy, the quantum corrections to dynamics resulting from localization of the nuclear wave function are computed approximately and included into selected degrees of freedom representing light particles where the quantum effects are expected to be the most pronounced. A massively parallel implementation, based on the message passing interface allows for efficient simulations of ensembles of thousands of trajectories at once. The QTES-DFTB dynamics approach is employed to study the role of quantum nuclear effects on the interaction of hydrogen with a model graphene sheet, revealing that neglect of nuclear effects can lead to an overestimation of adsorption.

Detection of hydrogen using graphene

Irradiation dynamics of a single graphene sheet bombarded by hydrogen atoms is studied in the incident energy range of 0.1 to 200 eV. Results for reflection, transmission, and adsorption probabilities, as well as effects of a single adsorbed atom to the electronic properties of graphene, are obtained by the quantum-classical Monte Carlo molecular dynamics within a self-consistent-charge-density functional tight binding formalism We compare these results with those, distinctly different, obtained by the classical molecular dynamics.PACS: 61.80.Az, 61.48.Gh, 61.80.Jh, 34.50.Dy.

Modeling Charge Transfer in Fullerene Collisions via Real-Time Electron Dynamics.

An approach for performing real-time dynamics of electron transfer in a prototype redox reaction that occurs in reactive collisions between neutral and ionic fullerenes is discussed. The quantum dynamical simulations show that the electron transfer occurs within 60 fs directly preceding the collision of the fullerenes, followed by structural changes and relaxation of electron charge. The consequences of real-time electron dynamics are fully elucidated for the far from equilibrium processes of collisions between neutral and multiply charged fullerenes.