Optimal Energy Transfer in Light-Harvesting Systems

Finite temperature dynamics in a polarized sub-Ohmic heat bath: A hierarchical equations of motion-tensor train study

The dynamics of the sub-Ohmic spin-boson model under polarized initial conditions at finite temperatures is investigated by employing both analytical tools and the numerically accurate hierarchical equations of motion-tensor train method. By analyzing the features of nonequilibrium dynamics, we discovered a bifurcation phenomenon, which separates two regimes of the dynamics. It is found that before the bifurcation time, increasing temperature slows down the population dynamics, while the opposite effect occurs after the bifurcation time. The dynamics is highly sensitive to both initial preparation of the bath and thermal effects.

Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra

Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes in molecules of biological and technological significance. The nonlinear femtosecond signals detected for such systems must be interpreted, but it is a challenging task for which theoretical simulations are often indispensable. Accurate simulations of transient absorption or two-dimensional electronic spectra are, however, computationally very expensive, prohibiting the wider adoption of existing first-principles methods. Here, we report an artificial-intelligence-enhanced protocol to drastically reduce the computational cost of simulating nonlinear time-resolved electronic spectra, which makes such simulations affordable for polyatomic molecules of increasing size. The protocol is based on the doorway-window approach for the on-the-fly surface-hopping simulations. We show its applicability for the prototypical molecule of pyrazine for which it produces spectra with high precision with respect to ab initio reference while cutting the computational cost by at least 95% compared to pure first-principles simulations.

Exploiting the photophysical features of DMAN template in ITQ-51 zeotype in the search for FRET energy transfer

The combination between photoactive molecules and inorganic structures is of great interest for the development of advanced materials in the field of optics. Particularly, zeotypes with extra-large pore size are attractive because they allow the encapsulation of bulky dyes. The microporous aluminophoshate Mg-ITQ-51 (IFO-type structure) represents an ideal candidate because of the synergic combination of two crucial features: the IFO framework itself, which is composed of non-interconnected one-dimensional extra-large elliptical channels with a diameter up to 11 Å able to host bulky guest species, and the particular organic structure-directing agent used for the synthesis (1,8-bis(dimethylamino)naphthalene, DMAN), which efficiently fills the IFO pores, and is itself a photoactive molecule with interesting fluorescence properties in the blue range of the visible spectrum, thus providing a densely-incorporated donor species for FRET processes. Besides, occlusion of DMAN dye in the framework triggers a notable improvement of its fluorescence properties by confinement effect. To extend the action of the material and to mimic processes such as photosynthesis in which FRET is essential, two robust laser dyes with bulky size, rhodamine 123 and Nile Blue, have been encapsulated for the first time in a zeolitic framework, together with DMAN, in a straightforward one-pot synthesis. Thus, photoactive systems with emission in the entire visible range have been achieved due to a partial FRET between organic chromophores protected in a rigid aluminophosphate matrix.

Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach

In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.

The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions

This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.

Determining interchromophore effects for energy transport in molecular networks using machine-learning algorithms

Nature uses chromophore networks, with highly optimized structural and energetic characteristics, to perform important chemical functions. Due to its modularity, predictable aggregation characteristics, and established synthetic protocols, structural DNA nanotechnology is a promising medium for arranging chromophore networks with analogous structural and energetic controls. However, this high level of control creates a greater need to know how to optimize the systems precisely. This study uses the system's modularity to produce variations of a coupled 14-Site chromophore network. It uses machine-learning algorithms and spectroscopy measurements to reveal the energy-transport roles of these Sites, paying particular attention to the cooperative and inhibitive effects they impose on each other for transport across the network. The physical significance of these patterns is contextualized, using molecular dynamics simulations and energy-transport modeling. This analysis yields insights about how energy transfers across the Donor-Relay and Relay-Acceptor interfaces, as well as the energy-transport pathways through the homogeneous Relay segment. Overall, this report establishes an approach that uses machine-learning methods to understand, in fine detail, the role that each Site plays in an optoelectronic molecular network.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Calculating non-linear response functions for multi-dimensional electronic spectroscopy using dyadic non-Markovian quantum state diffusion

We present a methodology for simulating multi-dimensional electronic spectra of molecular aggregates with coupling of electronic excitation to a structured environment using the stochastic non-Markovian quantum state diffusion (NMQSD) method in combination with perturbation theory for the response functions. A crucial aspect of our approach is that we propagate the NMQSD equation in a doubled system Hilbert space, but with the same noise. We demonstrate that our approach shows fast convergence with respect to the number of stochastic trajectories, providing a promising technique for numerical calculation of two-dimensional electronic spectra of large molecular aggregates.

Simulation of absorption spectra of molecular aggregates: A hierarchy of stochastic pure state approach

Simulation of spectroscopic observables for molecular aggregates with strong and structured coupling of electronic excitation to vibrational degrees of freedom is an important but challenging task. The Hierarchy of Pure States (HOPS) provides a formally exact solution based on local, stochastic trajectories. Exploiting the localization of HOPS for the simulation of absorption spectra in large aggregates requires a formulation in terms of normalized trajectories. Here, we provide a normalized dyadic equation where the ket- and bra-states are propagated in different electronic Hilbert spaces. This work opens the door to applying adaptive HOPS methods for the simulation of absorption spectra.

Funneling dynamics in a phenylacetylene trimer: Coherent excitation of donor excitonic states and their superposition

Funneling dynamics in conjugated dendrimers has raised great interest in the context of artificial light-harvesting processes. Photoinduced relaxation has been explored by time-resolved spectroscopy and simulations, mainly by semiclassical approaches or referring to open quantum systems methods, within the Redfield approximation. Here, we take the benefit of an ab initio investigation of a phenylacetylene trimer, and in the spirit of a divide-and-conquer approach, we focus on the early dynamics of the hierarchy of interactions. We build a simplified but realistic model by retaining only bright electronic states and selecting the vibrational domain expected to play the dominant role for timescales shorter than 500 fs. We specifically analyze the role of the in-plane high-frequency skeletal vibrational modes involving the triple bonds. Open quantum system non-adiabatic dynamics involving conical intersections is conducted by separating the electronic subsystem from the high-frequency tuning and coupling vibrational baths. This partition is implemented within a robust non-perturbative and non-Markovian method, here the hierarchical equations of motion. We will more precisely analyze the coherent preparation of donor states or of their superposition by short laser pulses with different polarizations. In particular, we extend the π-pulse strategy for the creation of a superposition to a V-type system. We study the relaxation induced by the high-frequency vibrational collective modes and the transitory dissymmetry, which results from the creation of a superposition of electronic donor states.

Simulation of Time- and Frequency-Resolved Four-Wave-Mixing Signals at Finite Temperatures: A Thermo-Field Dynamics Approach

We propose a new approach to simulate four-wave-mixing signals of molecular systems at finite temperatures by combining the multiconfigurational Ehrenfest method with the thermo-field dynamics theory. In our approach, the four-time correlation functions at finite temperatures are mapped onto those at zero temperature in an enlarged Hilbert space with twice the vibrational degrees of freedom. As an illustration, we have simulated three multidimensional spectroscopic signals, time- and frequency-resolved fluorescence spectra, transient-absorption pump-probe spectra, and electronic two-dimensional (2D) spectra at finite temperatures, for a conical intersection-mediated singlet fission model of a rubrene crystal. It is shown that a detailed dynamical picture of the singlet fission process can be extracted from the three spectroscopic signals. An increasing temperature leads to lower intensities of the signals and broadened vibrational peaks, which can be attributed to faster singlet-triplet population transfer and stronger bath-induced electronic dephasing at higher temperatures.

Efficient simulation of time- and frequency-resolved four-wave-mixing signals with a multiconfigurational Ehrenfest approach

We have extended the multiconfigurational Ehrenfest approach to the simulation of four-wave-mixing signals of systems involving multiple electronic and vibrational degrees of freedom. As an illustration, we calculate signals of three widely used spectroscopic techniques, time- and frequency-resolved fluorescence spectroscopy, transient absorption spectroscopy, and two-dimensional (2D) electronic spectroscopy, for a two-electronic-state, twenty-four vibrational-mode conical intersection model. It has been shown that all these three spectroscopic signals characterize fast population transfer from the higher excited electronic state to the lower excited electronic state. While the time- and frequency-resolved spectrum maps the wave packet propagation exclusively on the electronically excited states, the transient absorption and 2D electronic spectra reflect the wave packet dynamics on both electronically excited states and the electronic ground state. Combining trajectory-guided Gaussian basis functions and the nonlinear response function formalism, the present approach provides a promising general technique for the applications of various Gaussian basis methods to the calculations of four-wave-mixing spectra of polyatomic molecules.

Geometry Effects on Light-Harvesting Complex's Light Absorption and Energy Transfer in Purple Bacteria

Light-harvesting complexes (LHC) in photosynthetic organisms perform the major function of light absorption and energy transportation. Optical spectrum of LHC provides a detailed understanding of the molecular mechanisms involved in the excitation energy transfer (EET) processes, which has been widely studied. Here, we study how the geometric property of LHC in Rhodospirillum (Rs.) molischianum would affect its spectral characteristics and energy transfer process. By adopting the effective Hamiltonian and the dipole-dipole approximation, we calculate the exciton level structures for the LH2 ring and LH1 ring and the energy transfer time between different LHCs under various structural parameters and different rotational symmetries. Our numerical results show that the LHC's absorption peaks and the energy transfer time between different LHCs can be modified by changing the geometric configurations. Our study may be beneficial to the applications in designing highly efficient photovoltaic cell and other artificial photosynthetic systems.

Simulation of Femtosecond Phase-Locked Double-Pump Signals of Individual Light-Harvesting Complexes LH2

Recent phase-locked femtosecond double-pump experiments on individual light-harvesting complexes LH2 of purple bacteria at ambient temperature revealed undamped oscillatory responses on a time scale of at least 400 fs [ Hildner et al. Science 2013 , 340 , 1448 ]. Using an excitonic Hamiltonian for LH2 available in the literature, we simulate these signals numerically by a method that treats excitonic couplings and exciton-phonon couplings in a nonperturbative manner. The simulations provide novel insights into the origin of coherent dynamics in individual LH2 complexes.

Finite-temperature time-dependent variation with multiple Davydov states

The Dirac-Frenkel time-dependent variational approach with Davydov Ansätze is a sophisticated, yet efficient technique to obtain an accurate solution to many-body Schrödinger equations for energy and charge transfer dynamics in molecular aggregates and light-harvesting complexes. We extend this variational approach to finite temperature dynamics of the spin-boson model by adopting a Monte Carlo importance sampling method. In order to demonstrate the applicability of this approach, we compare calculated real-time quantum dynamics of the spin-boson model with that from numerically exact iterative quasiadiabatic propagator path integral (QUAPI) technique. The comparison shows that our variational approach with the single Davydov Ansätze is in excellent agreement with the QUAPI method at high temperatures, while the two differ at low temperatures. Accuracy in dynamics calculations employing a multitude of Davydov trial states is found to improve substantially over the single Davydov Ansatz, especially at low temperatures. At a moderate computational cost, our variational approach with the multiple Davydov Ansatz is shown to provide accurate spin-boson dynamics over a wide range of temperatures and bath spectral densities.

Reassessment of the Four-Point Approach to the Electron-Transfer Marcus-Hush Theory

The Marcus-Hush theory has been successfully applied to describe and predict the activation barriers and hence the electron-transfer (ET) rates in several physicochemical and biological systems. This theory assumes that in the ET reaction, the geometry of the free Gibbs energy landscape is parabolic, with equal curvature near the local minimum for both reactants and products. In spite of its achievements, more realistic models have included the assumption of the two parabolas having not the same curvature. This situation is analyzed by the Nelsen's four-point method. As a benchmark to compare the Marcus-Hush approximation to a precise calculation of the excitation energy, we studied the non-ET process of the electronic excitation of the aluminum dimer that has two local minima (³∑_g ^- and ³∏_u electronic states) and allows to obtain analytically the Marcus-Hush nonsymmetric parameters. We appraise the ability of the Marcus-Hush formula to approximate the analytical results by using several averages of the two reorganization energies associated with the forward and backward transitions and analyze the error. It is observed that the geometric average minimizes the relative error and that the analytical case is recovered. The main results of this paper are obtained by the application of the Nelsen's four-point method to compute the reorganization energies of a large set of potential π-conjugated molecules proposed for organic photovoltaic devices using the above-mentioned averages for the Marcus-Hush formula. The activation energies obtained with the geometric average are significantly larger for some donor-acceptor pairs in comparison with the previously employed arithmetic average, their differences being suitable for experimental testing.

Bio-Optics and Bio-Inspired Optical Materials

Through the use of the limited materials palette, optimally designed micro- and nanostructures, and tightly regulated processes, nature demonstrates exquisite control of light-matter interactions at various length scales. In fact, control of light-matter interactions is an important element in the evolutionary arms race and has led to highly engineered optical materials and systems. In this review, we present a detailed summary of various optical effects found in nature with a particular emphasis on the materials and optical design aspects responsible for their optical functionality. Using several representative examples, we discuss various optical phenomena, including absorption and transparency, diffraction, interference, reflection and antireflection, scattering, light harvesting, wave guiding and lensing, camouflage, and bioluminescence, that are responsible for the unique optical properties of materials and structures found in nature and biology. Great strides in understanding the design principles adapted by nature have led to a tremendous progress in realizing biomimetic and bioinspired optical materials and photonic devices. We discuss the various micro- and nanofabrication techniques that have been employed for realizing advanced biomimetic optical structures.

Fast, Accurate Simulation of Polaron Dynamics and Multidimensional Spectroscopy by Multiple Davydov Trial States

By employing the Dirac-Frenkel time-dependent variational principle, we study the dynamical properties of the Holstein molecular crystal model with diagonal and off-diagonal exciton-phonon coupling. A linear combination of the Davydov D1 (D2) ansatz, referred to as the "multi-D1 ansatz" ("multi-D2 ansatz"), is used as the trial state with enhanced accuracy but without sacrificing efficiency. The time evolution of the exciton probability is found to be in perfect agreement with that of the hierarchy equations of motion, demonstrating the promise the multiple Davydov trial states hold as an efficient, robust description of dynamics of complex quantum systems. In addition to the linear absorption spectra computed for both diagonal and off-diagonal cases, for the first time, 2D spectra have been calculated for systems with off-diagonal exciton-phonon coupling by employing the multiple D2 ansatz to compute the nonlinear response function, testifying to the great potential of the multiple D2 ansatz for fast, accurate implementation of multidimensional spectroscopy. It is found that the signal exhibits a single peak for weak off-diagonal coupling, while a vibronic multipeak structure appears for strong off-diagonal coupling.

Energy transfer cascade in bio-inspired chlorophyll-a/polyacrylamide hydrogel: towards a new class of biomimetic solar cells

Efficient solar energy harvesting in natural photosystem inspired chlorophyll-a/hydrogel based soft, simple system, revealing the effect of coherence-dephasing interpaly.