Proton conduction along a chain of water molecules. Development of a linear model and quantum dynamical investigations using the multiconfiguration time-dependent Hartree method

Resource Optimization for Quantum Dynamics with Tensor Networks: Quantum and Classical Algorithms

The exponential scaling of the quantum degrees of freedom with the size of the system is one of the biggest challenges in computational chemistry and particularly in quantum dynamics. We present a tensor network approach for the time-evolution of the nuclear degrees of freedom of multiconfigurational chemical systems at a reduced storage and computational complexity. We also present quantum algorithms for the resultant dynamics. To preserve the compression advantage achieved via tensor network decompositions, we present an adaptive algorithm for the regularization of nonphysical bond dimensions, preventing the potentially exponential growth of these with time. While applicable to any quantum dynamical problem, our method is particularly valuable for dynamical simulations of nuclear chemical systems. Our algorithm is demonstrated using ab initio potentials obtained for a symmetric hydrogen-bonded system, namely, the protonated 2,2'-bipyridine, and compared to exact diagonalization numerical results.

Synthesis of Hidden Subgroup Quantum Algorithms and Quantum Chemical Dynamics

We describe a general formalism for quantum dynamics and show how this formalism subsumes several quantum algorithms, including the Deutsch, Deutsch-Jozsa, Bernstein-Vazirani, Simon, and Shor algorithms as well as the conventional approach to quantum dynamics based on tensor networks. The common framework exposes similarities among quantum algorithms and natural quantum phenomena: we illustrate this connection by showing how the correlated behavior of protons in water wire systems that are common in many biological and materials systems parallels the structure of Shor's algorithm.

Hydroxide ion can move faster than an excess proton through one-dimensional water chains in hydrophobic narrow pores

Carbon nanotubes (CNT) are known to provide a hydrophobic, confined environment for water where its structure and dynamics can be very different from those of bulk water. In particular, narrow CNTs of the type (6,6) allow only a single one-dimensional (1D) chain of water molecules inside them, thus providing an idealized scenario to study motion in 1D along water chains. In the present study, we have investigated structural and dynamic behavior of water and also of an excess proton and hydroxide ion in water-filled narrow CNTs by means of ab initio molecular dynamics and combined quantum-classical simulations. The main focus of the present work is on the molecular mechanism and kinetics of hydronium and hydroxide ion migration along 1D water chains of different lengths in confinement. It is found that the hydrogen-bonded structures of water and the excess proton and hydroxide ion in CNTs can be very different from those in bulk, and these altered solvation structures play critical roles in determining the proton-transfer (PT) rates along water chains. For the present 1D chain systems, the hydroxide ion is found to migrate at a slightly faster rate than the excess proton, unlike their relative mobilities in bulk water. This faster migration of the hydroxide ion is found not only in CNTs with periodicity along the tube axis but also in isolated CNTs where the excess proton and the hydroxide ion are allowed to move under the influence of an electric field of an oppositely charged ion. The roles of rotational jumps and hydrogen-bond fluctuations in the PT events are discussed. In addition, the significance of hydrogen-bonding defects on the dynamics of an excess proton and hydroxide ion is also discussed for varying chain lengths.

The cytochrome ba3 oxygen reductase from Thermus thermophilus uses a single input channel for proton delivery to the active site and for proton pumping

The heme-copper oxygen reductases are redox-driven proton pumps that generate a proton motive force in both prokaryotes and mitochondria. These enzymes have been divided into 3 evolutionarily related groups: the A-, B- and C-families. Most experimental work on proton-pumping mechanisms has been performed with members of the A-family. These enzymes require 2 proton input pathways (D- and K-channels) to transfer protons used for oxygen reduction chemistry and for proton pumping, with the D-channel transporting all pumped protons. In this work we use site-directed mutagenesis to demonstrate that the ba(3) oxygen reductase from Thermus thermophilus, a representative of the B-family, does not contain a D-channel. Rather, it utilizes only 1 proton input channel, analogous to that of the A-family K-channel, and it delivers protons to the active site for both O2 chemistry and proton pumping. Comparison of available subunit I sequences reveals that the only structural elements conserved within the oxygen reductase families that could perform these functions are active-site components, namely the covalently linked histidine-tyrosine, the Cu(B) and its ligands, and the active-site heme and its ligands. Therefore, our data suggest that all oxygen reductases perform the same chemical reactions for oxygen reduction and comprise the essential elements of the proton-pumping mechanism (e.g., the proton-loading and kinetic-gating sites). These sites, however, cannot be located within the D-channel. These results along with structural considerations point to the A-propionate region of the active-site heme and surrounding water molecules as the proton-loading site.

Operation of the proton wire in green fluorescent protein. A quantum dynamics simulation

A nuclear quantum dynamical simulation of the proton shuttle operating in the green fluorescent protein has been carried out on a high-quality, high-dimensionality potential energy surface describing the photoactive pipi* excited state, and including motion of both the three protons and of the donor and acceptor atoms of the hydrogen bonds in a closed proton wire. The results of the simulations show that proton transfer along the wire is essentially concerted, synchronous, and very fast, with a substantial amount of the green fluorescent species forming within several tens of femtoseconds. In this regard, analysis of the population of the fluorescent species indicates that at least two dynamical regimes are present for its formation. Within the first hundreds of femtoseconds, dynamics is very fast and impulsive. Later on, a slower pace of formation appears. It is discussed that the two largest decay times for the protonated chromophore reported experimentally (Chattoraj, M.; King, B. A.; Bublitz, G. U.; Boxer, S. G. Proc. Natl. Acad. Sci. U.S.A. 1996, 93, 8362-8367) might correspond to some irreversible process occurring after formation of the fluorescent species, rather than to cleavage of the chromophore's phenolic O-H bond.