Distinguishing the roles of energy funnelling and delocalization in photosynthetic light harvesting.

Photosynthetic complexes improve the transfer of excitation energy from peripheral antennas to reaction centers in several ways. In particular, a downward energy funnel can direct excitons in the right direction, while coherent excitonic delocalization can enhance transfer rates through the cooperative phenomenon of supertransfer. However, isolating the role of purely coherent effects is difficult because any change to the delocalization also changes the energy landscape. Here, we show that the relative importance of the two processes can be determined by comparing the natural light-harvesting apparatus with counterfactual models in which the delocalization and the energy landscape are altered. Applied to the example of purple bacteria, our approach shows that although supertransfer does enhance the rates somewhat, the energetic funnelling plays the decisive role. Because delocalization has a minor role (and is sometimes detrimental), it is most likely not adaptive, being a side-effect of the dense chlorophyll packing that evolved to increase light absorption per reaction center.

Spectral dependence of the internal quantum efficiency of organic solar cells: effect of charge generation pathways.

The conventional picture of photocurrent generation in organic solar cells involves photoexcitation of the electron donor, followed by electron transfer to the acceptor via an interfacial charge-transfer state (Channel I). It has been shown that the mirror-image process of acceptor photoexcitation leading to hole transfer to the donor is also an efficient means to generate photocurrent (Channel II). The donor and acceptor components may have overlapping or distinct absorption characteristics. Hence, different excitation wavelengths may preferentially activate one channel or the other, or indeed both. As such, the internal quantum efficiency (IQE) of the solar cell may likewise depend on the excitation wavelength. We show that several model high-efficiency organic solar cell blends, notably PCDTBT:PC70BM and PCPDTBT:PC60/70BM, exhibit flat IQEs across the visible spectrum, suggesting that charge generation is occurring either via a dominant single channel or via both channels but with comparable efficiencies. In contrast, blends of the narrow optical gap copolymer DPP-DTT with PC70BM show two distinct spectrally flat regions in their IQEs, consistent with the two channels operating at different efficiencies. The observed energy dependence of the IQE can be successfully modeled as two parallel photodiodes, each with its own energetics and exciton dynamics but both having the same extraction efficiency. Hence, an excitation-energy dependence of the IQE in this case can be explained as the interplay between two photocurrent-generating channels, without recourse to hot excitons or other exotic processes.

Delocalisation enables efficient charge generation in organic photovoltaics, even with little to no energetic offset.

Organic photovoltaics (OPVs) are promising candidates for solar-energy conversion, with device efficiencies continuing to increase. However, the precise mechanism of how charges separate in OPVs is not well understood because low dielectric constants produce a strong attraction between the charges, which they must overcome to separate. Separation has been thought to require energetic offsets at donor-acceptor interfaces, but recent materials have enabled efficient charge generation with small offsets, or with none at all in neat materials. Here, we extend delocalised kinetic Monte Carlo (dKMC) to develop a three-dimensional model of charge generation that includes disorder, delocalisation, and polaron formation in every step from photoexcitation to charge separation. Our simulations show that delocalisation dramatically increases charge-generation efficiency, partly by enabling excitons to dissociate in the bulk. Therefore, charge generation can be efficient even in devices with little to no energetic offset, including neat materials. Our findings demonstrate that the underlying quantum-mechanical effect that improves the charge-separation kinetics is faster and longer-distance hops between delocalised states, mediated by hybridised states of exciton and charge-transfer character.

Light-Harvesting Efficiency Cannot Depend on Optical Coherence in the Absence of Orientational Order.

The coherence of light has been proposed as a quantum-mechanical control for enhancing light-harvesting efficiency. In particular, optical coherence can be manipulated by changing either the polarization state or the spectral phase of the light. Here, we show that, in weak light, light-harvesting efficiency cannot be controlled using any form of optical coherence in molecular light-harvesting systems and, more broadly, those comprising orientationally disordered subunits and operating on longer-than-ultrafast time scales. Under those conditions, optical coherence does not affect the light-harvesting efficiency, meaning that it cannot be used for control. Specifically, polarization-state control is lost in disordered samples or when the molecules reorient on the time scales of light harvesting, and spectral-phase control is lost when the efficiency is time-averaged over a period longer than the optical coherence time. In practice, efficiency is always averaged over long times, meaning that coherent optical control is only possible through polarization and only in systems with orientational order.

Mechanism of Delocalization-Enhanced Exciton Transport in Disordered Organic Semiconductors.

Large exciton diffusion lengths generally improve the performance of organic semiconductor devices, because they enable energy to be transported farther during the exciton lifetime. However, the physics of exciton motion in disordered organic materials is not fully understood, and modeling the transport of quantum-mechanically delocalized excitons in disordered organic semiconductors is a computational challenge. Here, we describe delocalized kinetic Monte Carlo (dKMC), the first model of three-dimensional exciton transport in organic semiconductors that includes delocalization, disorder, and polaron formation. We find that delocalization can dramatically increase exciton transport; for example, delocalization across less than two molecules in each direction can increase the exciton diffusion coefficient by over an order of magnitude. The mechanism for the enhancement is 2-fold: delocalization enables excitons to hop both more frequently and further in each hop. We also quantify the effect of transient delocalization (short-lived periods where excitons become highly delocalized) and show that it depends strongly upon the disorder and transition dipole moments.

Jumping Kinetic Monte Carlo: Fast and Accurate Simulations of Partially Delocalized Charge Transport in Organic Semiconductors.

Developing devices using disordered organic semiconductors requires accurate and practical models of charge transport. In these materials, charge transport occurs through partially delocalized states in an intermediate regime between localized hopping and delocalized band conduction. Partial delocalization can increase mobilities by orders of magnitude compared to those with conventional hopping, making it important for the design of materials and devices. Although delocalization, disorder, and polaron formation can be described using delocalized kinetic Monte Carlo (dKMC), it is a computationally expensive method. Here, we develop jumping kinetic Monte Carlo (jKMC), a model that approaches the accuracy of dKMC for modest amounts of delocalization (such as those found in disordered organic semiconductors), with a computational cost comparable to that of conventional hopping. jKMC achieves its computational performance by modeling conduction using identical spherical polarons, yielding a simple delocalization correction to the Marcus hopping rate that allows polarons to jump over their nearest neighbors. jKMC can be used in regimes of partial delocalization inaccessible to dKMC to show that modest delocalization can increase mobilities by as much as 2 orders of magnitude.

Predicting molecular vibronic spectra using time-domain analog quantum simulation.

Spectroscopy is one of the most accurate probes of the molecular world. However, predicting molecular spectra accurately is computationally difficult because of the presence of entanglement between electronic and nuclear degrees of freedom. Although quantum computers promise to reduce this computational cost, existing quantum approaches rely on combining signals from individual eigenstates, an approach whose cost grows exponentially with molecule size. Here, we introduce a method for scalable analog quantum simulation of molecular spectroscopy: by performing simulations in the time domain, the number of required measurements depends on the desired spectral range and resolution, not molecular size. Our approach can treat more complicated molecular models than previous ones, requires fewer approximations, and can be extended to open quantum systems with minimal overhead. We present a direct mapping of the underlying problem of time-domain simulation of molecular spectra to the degrees of freedom and control fields available in a trapped-ion quantum simulator. We experimentally demonstrate our algorithm on a trapped-ion device, exploiting both intrinsic electronic and motional degrees of freedom, showing excellent quantitative agreement for a single-mode vibronic photoelectron spectrum of SO2.

Simulating chemistry using quantum computers.

The difficulty of simulating quantum systems, well known to quantum chemists, prompted the idea of quantum computation. One can avoid the steep scaling associated with the exact simulation of increasingly large quantum systems on conventional computers, by mapping the quantum system to another, more controllable one. In this review, we discuss to what extent the ideas in quantum computation, now a well-established field, have been applied to chemical problems. We describe algorithms that achieve significant advantages for the electronic-structure problem, the simulation of chemical dynamics, protein folding, and other tasks. Although theory is still ahead of experiment, we outline recent advances that have led to the first chemical calculations on small quantum information processors.

Even a little delocalization produces large kinetic enhancements of charge-separation efficiency in organic photovoltaics.

In organic photovoltaics, charges can separate efficiently even if their Coulomb attraction is an order of magnitude greater than the available thermal energy. Delocalization has been suggested to explain this fact, because it could increase the initial separation of charges in the charge-transfer (CT) state, reducing their attraction. However, understanding the mechanism requires a kinetic model of delocalized charge separation, which has proven difficult because it involves tracking the correlated quantum-mechanical motion of the electron and the hole in large simulation boxes required for disordered materials. Here, we report the first three-dimensional simulations of charge-separation dynamics in the presence of disorder, delocalization, and polaron formation, finding that even slight delocalization, across less than two molecules, can substantially enhance the charge-separation efficiency, even starting with thermalized CT states. Delocalization does not enhance efficiency by reducing the Coulomb attraction; instead, the enhancement is a kinetic effect produced by the increased overlap of electronic states.

Polynomial-time quantum algorithm for the simulation of chemical dynamics.

The computational cost of exact methods for quantum simulation using classical computers grows exponentially with system size. As a consequence, these techniques can be applied only to small systems. By contrast, we demonstrate that quantum computers could exactly simulate chemical reactions in polynomial time. Our algorithm uses the split-operator approach and explicitly simulates all electron-nuclear and interelectronic interactions in quadratic time. Surprisingly, this treatment is not only more accurate than the Born-Oppenheimer approximation but faster and more efficient as well, for all reactions with more than about four atoms. This is the case even though the entire electronic wave function is propagated on a grid with appropriately short time steps. Although the preparation and measurement of arbitrary states on a quantum computer is inefficient, here we demonstrate how to prepare states of chemical interest efficiently. We also show how to efficiently obtain chemically relevant observables, such as state-to-state transition probabilities and thermal reaction rates. Quantum computers using these techniques could outperform current classical computers with 100 qubits.

Environmentally Improved Coherent Light Harvesting.

Coherence-enhanced light harvesting has not been directly observed experimentally, despite theoretical evidence that coherence can significantly enhance light-harvesting performance. The main experimental obstacle has been the difficulty in isolating the effect of coherence in the presence of confounding variables. Recent proposals for externally controlling coherence by manipulating the light's degree of polarization showed that coherent efficiency enhancements would be possible, but they were restricted to light-harvesting systems weakly coupled to their environment. Here, we show that increases in system-bath coupling strength can amplify coherent efficiency enhancements, rather than suppress them. This result dramatically broadens the range of systems that could be used to conclusively demonstrate coherence-enhanced light harvesting or to engineer coherent effects into artificial light-harvesting devices.

Analog quantum simulation of chemical dynamics.

Ultrafast chemical reactions are difficult to simulate because they involve entangled, many-body wavefunctions whose computational complexity grows rapidly with molecular size. In photochemistry, the breakdown of the Born-Oppenheimer approximation further complicates the problem by entangling nuclear and electronic degrees of freedom. Here, we show that analog quantum simulators can efficiently simulate molecular dynamics using commonly available bosonic modes to represent molecular vibrations. Our approach can be implemented in any device with a qudit controllably coupled to bosonic oscillators and with quantum hardware resources that scale linearly with molecular size, and offers significant resource savings compared to digital quantum simulation algorithms. Advantages of our approach include a time resolution orders of magnitude better than ultrafast spectroscopy, the ability to simulate large molecules with limited hardware using a Suzuki-Trotter expansion, and the ability to implement realistic system-bath interactions with only one additional interaction per mode. Our approach can be implemented with current technology; e.g., the conical intersection in pyrazine can be simulated using a single trapped ion. Therefore, we expect our method will enable classically intractable chemical dynamics simulations in the near term.

Classification of Coherent Enhancements of Light-Harvesting Processes.

Several kinds of coherence have recently been shown to affect the performance of light-harvesting systems, in some cases significantly improving their efficiency. Here, we classify the possible mechanisms of coherent efficiency enhancements, based on the types of coherence that can characterize a light-harvesting system and the types of processes these coherences can affect. We show that enhancements are possible only when coherences and dissipative effects are best described in different bases of states. Our classification allows us to predict a previously unreported coherent enhancement mechanism, where coherence between delocalized eigenstates can be used to localize excitons away from dissipation, thus reducing the rate of recombination and increasing efficiency.

Delocalised kinetic Monte Carlo for simulating delocalisation-enhanced charge and exciton transport in disordered materials.

Charge transport is well understood in both highly ordered materials (band conduction) or highly disordered ones (hopping conduction). In moderately disordered materials-including many organic semiconductors-the approximations valid in either extreme break down, making it difficult to accurately model the conduction. In particular, describing wavefunction delocalisation requires a quantum treatment, which is difficult in disordered materials that lack periodicity. Here, we present the first three-dimensional model of partially delocalised charge and exciton transport in materials in the intermediate disorder regime. Our approach is based on polaron-transformed Redfield theory, but overcomes several computational roadblocks by mapping the quantum-mechanical techniques onto kinetic Monte Carlo. Our theory, delocalised kinetic Monte Carlo (dKMC), shows that the fundamental physics of transport in moderately disordered materials is that of charges hopping between partially delocalised electronic states. Our results reveal why standard kinetic Monte Carlo can dramatically underestimate mobilities even in disordered organic semiconductors, where even a little delocalisation can substantially enhance mobilities, as well as showing that three-dimensional calculations capture important delocalisation effects neglected in lower-dimensional approximations.

Quantum algorithm for molecular properties and geometry optimization.

Quantum computers, if available, could substantially accelerate quantum simulations. We extend this result to show that the computation of molecular properties (energy derivatives) could also be sped up using quantum computers. We provide a quantum algorithm for the numerical evaluation of molecular properties, whose time cost is a constant multiple of the time needed to compute the molecular energy, regardless of the size of the system. Molecular properties computed with the proposed approach could also be used for the optimization of molecular geometries or other properties. For that purpose, we discuss the benefits of quantum techniques for Newton's method and Householder methods. Finally, global minima for the proposed optimizations can be found using the quantum basin hopper algorithm, which offers an additional quadratic reduction in cost over classical multi-start techniques.

Preparation of many-body states for quantum simulation.

While quantum computers are capable of simulating many quantum systems efficiently, the simulation algorithms must begin with the preparation of an appropriate initial state. We present a method for generating physically relevant quantum states on a lattice in real space. In particular, the present algorithm is able to prepare general pure and mixed many-particle states of any number of particles. It relies on a procedure for converting from a second-quantized state to its first-quantized counterpart. The algorithm is efficient in that it operates in time that is polynomial in all the essential descriptors of the system, the number of particles, the resolution of the lattice, and the inverse of the maximum final error. This scaling holds under the assumption that the wave function to be prepared is bounded or its indefinite integral is known and that the Fock operator of the system is efficiently simulatable.

Why are photosynthetic reaction centres dimeric?

All photosynthetic organisms convert solar energy into chemical energy through charge separation in dimeric reaction centres. It is unknown why early reaction centres dimerised and completely displaced their monomeric ancestors. Here, we discuss several proposed explanations for reaction-centre dimerism and conclude-with only weak assumptions about the primordial dimerisation event-that the most probable explanation for the dimerism is that it arose because it enhanced light-harvesting efficiency by deepening the excitonic trap, i.e., by enhancing the rate of exciton transfer from an antenna complex and decreasing the rate of back transfer. This effect would have outweighed the negative effect dimerisation would have had on charge transfer within the reaction centre. Our argument implies that dimerisation likely occurred after the evolution of the first antennas, and it explains why the lower-energy state of the special pair is bright.

Measuring Energetic Disorder in Organic Semiconductors Using the Photogenerated Charge-Separation Efficiency.

Quantifying energetic disorder in organic semiconductors continues to attract attention because of its significant impact on the transport physics of these technologically important materials. Here, we show that the energetic disorder of organic semiconductors can be determined from the relationship between the internal quantum efficiency of charge generation and the frequency of the incident light. Our results for a number of materials suggest that energetic disorder in organic semiconductors could be greater than previously reported, and we advance ideas as to why this may be the case.

Generalised Marcus theory for multi-molecular delocalised charge transfer.

Although Marcus theory is widely used to describe charge transfer in molecular systems, in its usual form it is restricted to transfer from one molecule to another. If a charge is delocalised across multiple donor molecules, this approach requires us to treat the entire donor aggregate as a unified supermolecule, leading to potentially expensive quantum-chemical calculations and making it more difficult to understand how the aggregate components contribute to the overall transfer. Here, we show that it is possible to describe charge transfer between groups of molecules in terms of the properties of the constituent molecules and couplings between them, obviating the need for expensive supermolecular calculations. We use the resulting theory to show that charge delocalisation between molecules in either the donor or acceptor aggregates can enhance the rate of charge transfer through a process we call supertransfer (or suppress it through subtransfer). The rate can also be enhanced above what is possible with a single donor and a single acceptor by judiciously tuning energy levels and reorganisation energies. We also describe bridge-mediated charge transfer between delocalised molecular aggregates. The equations of generalised Marcus theory are in closed form, providing qualitative insight into the impact of delocalisation on charge dynamics in molecular systems.

Anomalous Exciton Quenching in Organic Semiconductors in the Low-Yield Limit.

The dynamics of exciton quenching are critical to the operational performance of organic optoelectronic devices, but their measurement and elucidation remain ongoing challenges. Here, we present a method for quantifying small photoluminescence quenching efficiencies of organic semiconductors under steady-state conditions. Exciton quenching efficiencies of three different organic semiconductors, PC70BM, P3HT, and PCDTBT, are measured at different bulk quencher densities under continuous low-irradiance illumination. By implementing a steady-state bulk-quenching model, we determine exciton diffusion lengths for the studied materials. At low quencher densities we find that a secondary quenching mechanism is in effect, which is responsible for approximately 20% of the total quenched excitons. This quenching mechanism is observed in all three studied materials and exhibits quenching volumes on the order of several thousand cubic nanometers. The exact origin of this quenching process is not clear, but it may be indicative of delocalized excitons being quenched prior to thermalization.