Search

Following the observation of coherent oscillations in non-linear spectra of photosynthetic pigment protein complexes, particularly phycobilliprotein such as PC645, coherent vibronic transport has been suggested as a design principle for novel light harvesting materials operating at room temperature. Vibronic transport between energetically remote pigments is coherent when the presence of a resonant vibration supports transient delocalization between the pair of electronic excited states. Here, we establish the mechanism of vibronic transport for a model heterodimer across a wide range of molecular parameter values. The resulting mechanistic map demonstrates that the molecular parameters of phycobiliproteins in fact support incoherent vibronic transport. This result points to an important design principle: incoherent vibronic transport is more efficient than a coherent mechanism when energetic disorder exceeds the coupling between the donor and vibrationally excited acceptor states. Finally, our results suggest that the role of coherent vibronic transport in pigment protein complexes should be reevaluated.

In this work we introduce PHOENICS, a probabilistic global optimization algorithm combining ideas from Bayesian optimization with concepts from Bayesian kernel density estimation. We propose an inexpensive acquisition function balancing the explorative and exploitative behavior of the algorithm. This acquisition function enables intuitive sampling strategies for an efficient parallel search of global minima. The performance of PHOENICS is assessed via an exhaustive benchmark study on a set of 15 discrete, quasi-discrete and continuous multidimensional functions. Unlike optimization methods based on Gaussian processes (GP) and random forests (RF), we show that PHOENICS is less sensitive to the nature of the co-domain, and outperforms GP and RF optimizations. We illustrate the performance of PHOENICS on the Oregonator, a difficult case-study describing a complex chemical reaction network. We demonstrate that only PHOENICS was able to reproduce qualitatively and quantitatively the target dynamic behavior of this nonlinear reaction dynamics. We recommend PHOENICS for rapid optimization of scalar, possibly non-convex, black-box unknown objective functions.

Automatic Differentiation (AD) is a powerful tool that allows calculating derivatives of implemented algorithms with respect to all of their parameters up to machine precision, without the need to explicitly add any additional functions. Thus, AD has great potential in quantum chemistry, where gradients are omnipresent but also difficult to obtain, and researchers typically spend a considerable amount of time finding suitable analytical forms when implementing derivatives. Here, we demonstrate that automatic differentiation can be used to compute gradients with respect to any parameter throughout a complete quantum chemistry method. We implement DiffiQult, a fully autodifferentiable Hartree-Fock (HF) algorithm, which serves as a proof-of-concept that illustrates the capabilities of AD for quantum chemistry. We leverage the obtained gradients to optimize the parameters of one-particle basis sets in the context of the floating Gaussian framework.

Self-assembled aggregates of pigment molecules are potential building blocks for excitonic circuits that find their application in energy conversion and optical signal processing. Recent experimental studies of one-dimensional heterotriangulene supramolecular aggregates suggested that singlet excitons in these structures can propagate on several micron distances. We explore this possibility theoretically by combining electronic structure calculations with microscopic models for exciton transport. Detailed characterization of the structural disorder and exciton decoherence is provided. We argue that conventional exciton transport models give about an order of magnitude shorter estimates for the exciton propagation length which suggest that there are other possible explanations of the experimental results., Comment: 22 pages, 4 figures

Understanding the relationship between the structure of light-harvesting systems and their excitation energy transfer properties is of fundamental importance in many applications including the development of next generation photovoltaics. Natural light harvesting in photosynthesis shows remarkable excitation energy transfer properties, which suggests that pigment-protein complexes could serve as blueprints for the design of nature inspired devices. Mechanistic insights into energy transport dynamics can be gained by leveraging numerically involved propagation schemes such as the hierarchical equations of motion (HEOM). Solving these equations, however, is computationally costly due to the adverse scaling with the number of pigments. Therefore virtual high-throughput screening, which has become a powerful tool in material discovery, is less readily applicable for the search of novel excitonic devices. We propose the use of artificial neural networks to bypass the computational limitations of established techniques for exploring the structure-dynamics relation in excitonic systems. Once trained, our neural networks reduce computational costs by several orders of magnitudes. Our predicted transfer times and transfer efficiencies exhibit similar or even higher accuracies than frequently used approximate methods such as secular Redfield theory

We present a study on the evolution of the Fenna-Matthews-Olson bacterial photosynthetic pigment-protein complex. This protein complex functions as an antenna. It transports absorbed photons - excitons - to a reaction center where photosynthetic reactions initiate. The efficiency of exciton transport is therefore fundamental for the photosynthetic bacterium's survival. We have reconstructed an ancestor of the complex to establish whether coherence in the exciton transport was selected for or optimized over time. We have also investigated on the role of optimizing free energy variation upon folding in evolution. We studied whether mutations which connect the ancestor to current day species were stabilizing or destabilizing from a thermodynamic view point. From this study, we established that most of these mutations were thermodynamically neutral. Furthermore, we did not see a large change in exciton transport efficiency or coherence and thus our results predict that exciton coherence was not specifically selected for.

Mechanisms controlling excitation energy transport (EET) in light-harvesting complexes remain controversial. Following the observation of long-lived beats in two-dimensional electronic spectroscopy of PC645, vibronic coherence, the delocalization of excited states between pigments supported by a resonant vibration, has been proposed to enable direct down-conversion from the highest-energy states to the lowest-energy pigments. Here, we instead show that for phycobiliprotein PC645 an incoherent vibronic transport mechanism is at play. We quantify the solvation dynamics of individual pigments using ab initio QM/MM nuclear dynamics. Our atomistic spectral densities reproduce experimental observations ranging from absorption and fluorescence spectra to the timescales and selectivity of down-conversion observed in transient absorption measurements. We demonstrate that bilin solvation controls EET pathways and that direct down-conversion proceeds incoherently, enhanced by large reorganization energies and a broad collection of high-frequency vibrations. We thus suggest that engineering local solvation dynamics represents a potential design principle for nanoscale control of EET.

We have performed time-dependent wave packet simulations of realistic Aharonov-Bohm (AB) devices with a quantum dot embedded in one of the arms of the interferometer. The AB ring can function as a measurement device for the intrinsic transmission phase through the quantum dot, however, care has to be taken in analyzing the influence of scattering processes in the junctions of the interferometer arms. We consider a harmonic quantum dot and show how the Darwin-Fock spectrum emerges as a unique pattern in the interference fringes of the AB oscillations., Comment: 17 pages, 5 figures. Published version

The investigation of energy transfer properties in photosynthetic multi-protein networks gives insight into their underlying design principles.Here, we discuss excitonic energy transfer mechanisms of the photosystem II (PS-II) C$_2$S$_2$M$_2$ supercomplex, which is the largest isolated functional unit of the photosynthetic apparatus of higher plants.Despite the lack of a decisive energy gradient in C$_2$S$_2$M$_2$, we show that the energy transfer is directed by relaxation to low energy states. C$_2$S$_2$M$_2$ is not organized to form pathways with strict energetic downhill transfer, which has direct consequences on the transfer efficiency, transfer pathways and transfer limiting steps. The exciton dynamics is sensitive to small structural changes, which, for instance, are induced by the reorganization of vibrational coordinates. In order to incorporate the reorganization process in our numerical simulations, we go beyond rate equations and use the hierarchically coupled equation of motion approach (HEOM). While transfer from the peripherical antenna to the proteins in proximity to the reaction center occurs on a faster time scale, the final step of the energy transfer to the RC core is rather slow, and thus the limiting step in the transfer chain. Our findings suggest that the structure of the PS-II supercomplex guarantees photoprotection rather than optimized efficiency., Comment: 23 pages, 6 figures

With quantum computers being out of reach for now, quantum simulators are the alternative devices for efficient and more exact simulation of problems that are challenging on conventional computers. Quantum simulators are classified into analog and digital, with the possibility of constructing "hybrid" simulators by combining both techniques. In this paper, we focus on analog quantum simulators of open quantum systems and address the limit that they can beat classical computers. In particular, as an example, we discuss simulation of the chlorosome light-harvesting antenna from green sulfur bacteria with over 250 phonon modes coupled to each electronic state. Furthermore, we propose physical setups that can be used to reproduce the quantum dynamics of a standard and multiple-mode Holstein model. The proposed scheme is based on currently available technology of superconducting circuits consist of flux qubits and quantum oscillators.

The theoretical and experimental study of energy transfer in photosynthesis has revealed an interesting transport regime, which lies at the borderline between classical transport dynamics and quantum-mechanical interference effects. Dissipation is caused by the coupling of electronic degrees of freedom to vibrational modes and leads to a directional energy transfer from the antenna complex to the target reaction-center. The dissipative driving is robust and does not rely on fine-tuning of specific vibrational modes. For the parameter regime encountered in the biological systems new theoretical tools are required to directly compare theoretical results with experimental spectroscopy data. The calculations require to utilize massively parallel graphics processor units (GPUs) for efficient and exact computations., Comment: 20 pages, submitted to the AIP conference proceedings of the Latin American School of Physics Marcos Moshinsky (ELAF 2013)

A numerically exact Monte Carlo scheme for calculation of open quantum system dynamics is proposed and implemented. The method consists of a Monte-Carlo summation of a perturbation expansion in terms of trajectories in Liouville phase-space with respect to the coupling between the excited states of the molecule. The trajectories are weighted by a complex decoherence factor based on the second-order cumulant expansion of the environmental evolution. The method can be used with an arbitrary environment characterized by a general correlation function and arbitrary coupling strength. It is formally exact for harmonic environments, and it can be used with arbitrary temperature. Time evolution of an optically excited Frenkel exciton dimer representing a molecular exciton interacting with a charge transfer state is calculated by the proposed method. We calculate the evolution of the optical coherence elements of the density matrix and linear absorption spectrum, and compare them with the predictions of standard simulation methods., Comment: 11 pages, 6 figures

Quantum transport and other phenomena are typically modeled by coupling the system of interest to an environment, or bath, held at thermal equilibrium. Realistic bath models are at least as challenging to construct as models for the quantum systems themselves, since they must incorporate many degrees of freedom that interact with the system on a wide range of timescales. Owing to computational limitations, the environment is often modeled with simple functional forms, with a few parameters fit to experiment to yield semi-quantitative results. Growing computational resources have enabled the construction of more realistic bath models from molecular dynamics (MD) simulations. In this paper, we develop a numerical technique to construct these atomistic bath models with better accuracy and decreased cost. We apply a novel signal processing technique, known as super-resolution, combined with a dictionary of physically-motivated bath modes to derive spectral densities from MD simulations. Our approach reduces the required simulation time and provides a more accurate spectral density than can be obtained via standard Fourier transform methods. Moreover, the spectral density is provided as a convenient closed-form expression which yields an analytic time-dependent bath kernel. Exciton dynamics of the Fenna-Matthews-Olsen light-harvesting complex are simulated with a second order time-convolutionless master equation, and spectral densities constructed via super-resolution are shown to reproduce the dynamics using only a quarter of the amount of MD data., Comment: 8 pages, 3 figures

The prevalence of long-lasting oscillatory signals in the 2d echo-spectroscopy of light-harvesting complexes has led to a search for possible mechanisms. We investigate how two causes of oscillatory signals are intertwined: (i) electronic coherences supporting delocalized wave-like motion, and (ii) narrow bands in the vibronic spectral density. To disentangle the vibronic and electronic contributions we introduce a time-windowed Fourier transform of the signal amplitude. We find that 2d spectra can be dominated by excitations of pathways which are absent in excitonic energy transport. This leads to an underestimation of the life-time of electronic coherences by 2d spectra., Comment: 6 pages, 6 figures

The observed prevalence of oscillatory signals in the spectroscopy of biological light-harvesting complexes at ambient temperatures has led to a search for mechanisms supporting coherent transport through larger molecules in noisy environments. We demonstrate a generic mechanism supporting long-lasting electronic coherence up to 0.3 ps at a temperature of 277 K. The mechanism relies on two properties of the spectral density: (i) a large dissipative coupling to a continuum of higher-frequency vibrations required for efficient transport and (ii) a small slope of the spectral density at zero frequency., Comment: 7 pages

Recent experimental observations of time-dependent beatings in the two-dimensional echo-spectra of light-harvesting complexes at ambient temperatures have opened up the question whether coherence and wave-like behaviour plays a significant role in photosynthesis. We perform a numerical study of the absorption and echo-spectra of the Fenna-Matthews-Olson (FMO) complex in chlorobium tepidum and analyse the requirements in the theoretical model needed to reproduce beatings in the calculated spectra. The energy transfer in the FMO pigment-protein complex is theoretically described by an exciton Hamiltonian coupled to a phonon bath which account for the pigments electronic and vibrational excitations respectively. We use the hierarchical equations of motions method to treat the strong couplings in a non-perturbative way. We show that the oscillations in the two-dimensional echo-spectra persist in the presence of thermal noise and static disorder., Comment: updated text and references, corrected figures 4-7, 20 pages, 7 figures

Excitonic models of light-harvesting complexes, where the vibrational degrees of freedom are treated as a bath, are commonly used to describe the motion of the electronic excitation through a molecule. Recent experiments point toward the possibility of memory effects in this process and require to consider time non-local propagation techniques. The hierarchical equations of motion (HEOM) were proposed by Ishizaki and Fleming to describe the site-dependent reorganization dynamics of protein environments (J. Chem. Phys., 130, p. 234111, 2009), which plays a significant role in photosynthetic electronic energy transfer. HEOM are often used as a reference for other approximate methods, but have been implemented only for small systems due to their adverse computational scaling with the system size. Here, we show that HEOM are also solvable for larger systems, since the underlying algorithm is ideally suited for the usage of graphics processing units (GPU). The tremendous reduction in computational time due to the GPU allows us to perform a systematic study of the energy-transfer efficiency in the Fenna-Matthews-Olson (FMO) light-harvesting complex at physiological temperature under full consideration of memory-effects. We find that approximative methods differ qualitatively and quantitatively from the HEOM results and discuss the importance of finite temperature to achieve high energy-transfer efficiencies., Comment: 14 pages; Journal of Chemical Theory and Computation (2011)

Quantum coherent properties of electrons can be studied in Aharonov-Bohm (AB) interferometers. We investigate both experimentally and theoretically the transmission phase evolution in a four-terminal quasi-one-dimensional AlGaAs/GaAs-based waveguide AB ring. As main control parameter besides the magnetic field, we tune the Fermi wave number along the pathways using a top-gate. Our experimental results and theoretical calculations demonstrate the strong influence of the measurement configuration upon the AB-resistance-oscillation phase in a four-terminal device. While the non-local setup displays continuous phase shifts of the AB oscillations, the phase remains rigid in the local voltage-probe setup. Abrupt phase jumps are found in all measurement configurations. We analyze the phase shifts as functions of the magnetic field and the Fermi energy and provide a detailed theoretical model of the device. Scattering and reflections in the arms of the ring are the source of abrupt phase jumps by $\pi$., Comment: 8 pages, 5 figures

Wave packets provide a well established and versatile tool for studying time-dependent effects in molecular physics. Here, we demonstrate the application of wave packets to mesoscopic nanodevices at low temperatures. The electronic transport in the devices is expressed in terms of scattering and transmission coefficients, which are efficiently obtained by solving an initial value problem (IVP) using the time-dependent Schroedinger equation. The formulation as an IVP makes non-trivial device topologies accessible and by tuning the wave packet parameters one can extract the scattering properties for a large range of energies., Comment: 12 pages, 4 figures

We study the quantum Hall effect (QHE) in graphene based on the current injection model. In our model, the presence of disorder, the edge-state picture, extended states and localized states, which are believed to be indispensable ingredients in describing the QHE, do not play an important role. Instead the boundary conditions during the injection into the graphene sheet, which are enforced by the presence of the Ohmic contacts, determine the current-voltage characteristics., Comment: 4 pages, 3 figures, rewritten, role of contacts for boundary conditions in small devices