# More accurate and efficient bath spectral densities from super-resolution

###### Abstract

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.

## I Introduction

Irreversible processes such as solvation, energy transfer, and chemical binding have received renewed interest in recent years. Because these processes involve large systems with many degrees of freedom, the typical approach to studying these processes is the open quantum systems formalism, in which the degrees of freedom are partitioned into a system of interest and a bath held at thermal equilibrium Breuer:2007wp (); Shenvi:2008tr (). It is commonly assumed that the system only couples weakly to the bath, making the precise nature of the bath a secondary concern in the physical theory. For example, in studying the energy transfer dynamics in a system of chromophores embedded in a protein framework, each chromophore is individually coupled to many thousands of atoms in the protein, but the system-bath formalism dramatically simplifies all of these couplings in order to make the dynamics tractable 2013JChPh.138k4102B (); Berkelbach:2013vw (); Singh:2011fl (). Renewed interest in the strong and intermediate coupling region, relevant for energy transfer in the exciton dynamics of light-harvesting complexes, has lead to various studies Jang:2007cz (); Mohseni:2008gp (); Plenio:2008ff (); Rebentrost:2009hu (); Cao:2009cc (); Ishizaki:2009tt (); Ishizaki:2011cx (); Sarovar:2010hs (); Abramavicius:2010et (); Wu:2010bg (); Moix:2011gi (); Kreisbeck:2011dh (); Skochdopole:2011gh (); Ritschel:2011ht (); Rebentrost:2011hc (); Singh:2011fl (); Pachon:2012fm (); Vlaming:2012hv (); Caruso:2012gf (); Zhu:2011ea (); Roden:2009gr (); Olbrich:2010ce (); Olbrich:2011hc (); Hein:2012vn () on the precise influence of the bath on the higher systems. Higher order phonon processes, non-Markovian effects and structures in the exciton-phonon coupling change the energy transfer doi:10.1021/jp304649c (); Chin:2013uh (); Kreisbeck:2012ui (). Thus, details in the bath are relevant and need to be taken into account in realistic simulations. Accordingly, our goal in this paper is to apply a recent signal-processing technique known as super-resolution to obtain realistic atomistic models of environments containing thousands of atoms at feasible computational expense. With these atomistic bath models in hand, one can begin to evaluate the importance of a realistic bath model in a physical theory.

In the approach to open quantum systems employed in this work, we model the bath by an ensemble of noninteracting harmonic oscillators. The central mathematical object of such a model is the spectral density, , which gives the frequency-dependent strength of system-bath coupling. The spectral density can be understood as the density of bath oscillator states at each frequency. Owing to computational limitations, most studies of open quantum systems assume an extremely simple functional form for the spectral density, such as a single broad peak covering all relevant excitonic transitions of the system. With the goal of providing more physically accurate bath models and dynamics, Valleau et al. has previously obtained atomistic spectral densities for the Fenna-Matthews-Olson (FMO) complex from combined Molecular Dynamics (MD) Valleau:2012ig (); Tuckerman:2010wy () and Time-Dependent Density Functional Theory (TDDFT) Runge:1984us () simulations. However, the difficulty of this more realistic approach is the high computational cost of running expensive TDDFT calculations at every step in an MD simulation. In order to obtain a spectral density of sufficient resolution, the MD-TDDFT simulation must be run for over 40 picoseconds (ps) Mallat:2008vn (), which may become computationally intractable for larger systems.

To make progress, we first observe that a typical vibrational bath is not an arbitrary function but rather a relatively sparse collection of damped harmonic oscillators. Sparsity enables us to apply a novel numerical technique known as super-resolution in order to reconstruct the spectral density from much shorter MD-TDDFT simulations. Super-resolution has been applied to a broad range of scientific problems, including image Freeman:2002va () and video compression Patti:1997vq (), image denoising Elad:1997un (), astronomy 2005A&A…436..373P (), microscopy Mccutchen:1967uf (), and medical imaging 5193030 (). To our knowledge, this paper is the first application of super-resolution to quantum dynamics. Super-resolution provides a provably convergent algorithm for the reconstruction of signals from limited time-domain measurements using a total variation minimization procedure. Super-resolution is related to compressed sensing 5272200 (); Donoho:ci (); 5288845 (); 4770164 (); 4959603 (); MRM:MRM21391 (); 6426647 (); Tuma:2009gb (); Shabani:2009de (); 5419072 (); 4472247 (); Coulter:2010wx (). Compressed sensing is a technique designed to recover sparse signals from randomly-sampled data by minimizing the norm of an underdetermined system of linear equations. Compressed sensing works by finding the sparsest signals consistent with the underdetermined system of equations. This usually involves an optimization problem. Despite its success in many applications, the -norm minimization of compressed sensing can result in spurious signals as it emphasizes the sparsity of the solution only. Super-resolution is a numerical method that shares the spirit of compressed sensing. The difference between superresolution and compressed sensing stems from both the choice of objective function and sampling technique. It was developed to recover sparse signals from nonrandomly undersampled data. By minimizing the -norm of the gradient of the function in addition to the norm of the function itself, super-resolution allows for smoother solutions to the sampling problem CPA:CPA21455 (); BioucasDias:2007tp (); 4378902 ().

Because of the ample experimental and theoretical data to compare against Valleau:2012ig (); Moix:2011gi (); Vulto:1998bu (); Olbrich:2011hc (); Adolphs:2006ey (); Mohseni:2011wg (); Chen:2013wz (); YuenZhou:2012vf (); Ritschel:2011ht (); Shim2012649 (); Pachon:2012fm (); Vlaming:2012hv (); Plenio:2008ff (); Singh:2011fl (); Kolli:2011vy (); Rebentrost:2009hu (); Skochdopole:2011gh (); Kreisbeck:2011dh (); Kreisbeck:2012ui (); Zhu:2011ea (); Cao:2009cc (); Sarovar:2010hs (); Ishizaki:2009tt (); Jang:2011vc (); Ishizaki:2009uh (), we apply super-resolution to the FMO light-harvesting complex of C. tepidium but emphasize that this technique is broadly applicable. While this paper focuses on a vibrational bath which perturbs the energies of molecular electronic states, the techniques we introduce are generic for any model of a bath which is based on time-correlation functions.

## Ii Super-resolution of Spectral Densities

In this section, we briefly review the procedure for simulating the dynamics of open quantum systems and computing spectral densities from combined MD-TDDFT simulations. We then apply the theory of super-resolution to accelerate and improve the accuracy of these computations. Computing spectral densities from atomistic calculations, rather than from semi-empirical functional forms, enables the inclusion of molecular vibrations and other physical effects (such as solvation effects) to produce a more realistic bath model Valleau:2012ig (). Super-resolution, in turn, brings the construction of these atomistic bath models into the realm of computational feasibility.

Armed with our more realistic bath model, we will employ a second-order time-convolutionless master equation (TCL-2) to simulate the dynamics of FMO monomer, allowing us to evaluate the physical impact of different approximations to the spectral density. TCL-2 includes non-Markovian effects up to second order in the system-bath coupling. By comparing TCL-2 with exact methods like the hierarchical equations of motion (HEOM) JPSJ.58.101 () we show that most of the relevant effects of the structured spectral density of the FMO complex are captured by TCL-2. Here we use TCL-2, since it is numerically more treatable than HEOM, in particular for structured spectral densities where HEOM becomes cumbersome and requires a high performance GPU implementation Kreisbeck:2011dh (); Kreisbeck:2012ui (); 16106 (). We employ the equation of motion Breuer:2007wp (); Rebentrost:2011vh (); Rebentrost:2009hu (); Shim2012649 (); Mohseni:2008gp (); Zhu:2011ea (); 2006JChPh.125j4906P (); 2011PhRvB..83k5416T (); Ahn:1994ww (); HeinzPeter:2000ul (); Shabani:2005wa (); Smirne:2010te ():

(II.1) | |||

(II.2) | |||

where is the system Hamiltonian, is the system density matrix, is our bath kernel, the subscript indicates that we are in the interaction picture, the summation runs over all sites, and is the spectral density computed via super-resolution Kolli:2011vy (); Kolli:2011ki (); Jang:2011vc (); Kleinekathofer:2004tx (); Kreisbeck:2011dh (). The bath kernel is heavily dependent on our spectral density, causing it to play a central role in our dynamics. Therefore, a more physical bath picture should provide more physically intuitive dynamics.

In our atomistic bath model, molecular vibrations in the environment (e.g. a protein framework or solvation effects) create fluctuations in the energy gaps between the ground and excited states of the system (e.g. a set of chromophores). These time-dependent energy gaps are computed from TDDFT calculations run on each of the chromophores at each step of the MD simulation. The key object in the computation of spectral densities is the correlation function of the energy gap time series,

(II.3) |

where is the time-dependent energy gap between the ground and the first excited state of the system (as calculated with TDDFT), is the density matrix of the bath at thermal equilibrium, and is the correlation function obtained after tracing over all the modes of the bath. We discretize this equation by using an unbiased autocorrelation function,

(II.4) |

where is the mean energy gap and and denote discrete time indices. Note that involves comparing energy gaps that are time steps apart ( and ), and is the total number of included comparisons.

The frequency-dependent spectral density, , is typically obtained by computing the Fourier transform of the correlation function Valleau:2012ig (). From the definition of above, it is easy to check that the correlation function is real and symmetric ( = ), which implies that the Fourier transform should be real and symmetric as well. Because quantum mechanical spectral densities must instead be antisymmetric and obey detailed balance, it is necessary to introduce a prefactor that enforces these two properties. Many choices are possible Berens:1981uy (), but Valleau et al. have previously shown that a harmonic prefactor, , produces the most physical temperature dependence Valleau:2012ig (). With this choice, the spectral density becomes the cosine transform

(II.5) |

which characterizes the frequency-dependent coupling strength of the system to all of the nuclear vibrational modes.

The standard approach to performing this integral is the fast Fourier transform. Unfortunately, the fast Fourier transform requires sampling on a uniform grid at the Shannon sampling rate. This means that a relatively long time series, , must be computed in order to obtain good resolution of the spectral density in the frequency domain Candes:eq (); Mallat:2008vn (). Given the computational cost of MD simulations, and the even greater expense of running TDDFT calculations on top of these simulations, any method which can reduce the required length of the time series unplugs the computational bottleneck in deriving physically-accurate atomistic spectral densities. That is our main goal in this paper.

While reducing the amount of time required to reproduce we also choose a basis of functions which has a convenient physical form. When decomposed into a basis of damped cosines,

(II.6) |

the function is smooth and sparse. This allows for the use of the machinery of super-resolution.

To apply the super-resolution method, we discretize in time and cast our task as an inversion problem

(II.7) |

where we seek the basis expansion coefficients and have assumed Einstein summation convention over repeated indices. This can be rewritten as

(II.8) |

where

(II.9) |

is a matrix of damped cosines, and is the set of basis coefficients we seek to recover.

The central idea of super-resolution is that the sparsity of enables its full recovery even when the system is underdetermined, which is to say the number of time samples is significantly smaller than the number of total expansion coefficients we seek to recover. Hence, we can recover the expansion coefficients on a dense grid of frequencies and damping coefficients from fewer time samples . Of the many possible solutions to our underdetermined system, super-resolution simply selects a balance between the smoothest and sparsest (with an emphasis on smoothness) set of basis expansion coefficients. Formally, this is done by finding the vector that minimizes

(II.10) | |||

where the subscript 1 represents the norm (sum of absolute values), represents a sparsity penalty, represents the total variation norm, and represents the solution tolerance. By minimizing , or total variation term, we are enforcing smoothness in the time domain on the reconstructed signal. This throws out the “peaky” solutions that can appear with compressed sensing Andrade28082012 (); Sanders:2012tk (). The total variation norm also provides us with a provably exact technique for recovering peak position at the expense of peak amplitude CPA:CPA21455 (), which solves one of the issues seen previously with compressed sensing Andrade28082012 ().

Recovering the expansion coefficients in this manner by solving an underdetermined matrix inversion problem takes advantage of the natural sparsity of the problem and, as we will see in the next section, enables the construction of a well-resolved spectral density with far less time-domain data. Even more attractive, with the coefficients in hand, it is possible to construct an analytical representation of the spectral density by taking the cosine transform of the basis functions and applying the appropriate prefactors:

(II.11) |

where the Einstein summation convention has again been assumed. This is an analytical representation of the spectral density in Drude-Lorentz form, and it explicitly provides the oscillation frequencies which characterize the system-bath coupling. We note that the Drude-Lorentz basis naturally provides us with a width parameter, , that can be understood as the lifetime of oscillations in the bath. This is seen by examining the time dependent formula, Eq. (II.7), where this parameter determines the strength of damping. It is important to note that in the limit as , we recover the cosine basis in the time domain and a Dirac delta distribution in the frequency domain. By using this super-resolution technique in concert with the Drude-Lorentz basis, we see that we can recover a small set of peaks with physically-relevant information. Additionally, the parameters that characterize the Drude-Lorentz spectral densities can be input directly into both TCL-2 and HEOM without any additional parameter fitting or numerical integration.

## Iii Numerical Methods

We employ the proposed Drude-Lorentz super-resolution method described above and apply it to a monomer of the Fenna-Matthews-Olsen (FMO) photosynthetic energy transfer complex of the green-sulfur bacterium C. tepidium. The FMO monomer is a system of seven chlorophyll molecules which are excitonically coupled to each other, as well as to the vibrations of the atoms in the protein framework. It functions as a molecular excitonic wire, passing excitons from the light harvesting antenna complex to the reaction center, where a biochemical cascade is initiated.

To create spectral densities for the FMO complex, we use the MD-TDDFT results of Shim et al. Shim2012649 (). The calculations were done in an isothermal-isobaric ensemble at 77 K using the AMBER force field Cornell:1995td (); Ceccarelli:2003ts (). These calculations began with a 2 ns equilibration before performing the production computations. The production steps ran for a total of 40 picoseconds with a 2 femtosecond timestep, and the optical gap was calculated for each fragment every 4 femtoseconds using TDDFT with the BLYP Becke:1988tx (); Miehlich:1989vq (); Lee:1988ub () functional in the 3-21G basis set in Q-Chem Shao:2006wf ().

To perform super-resolution numerically, we require an algorithm which minimizes the total variation norm to solve the minimization problem described by eq. (II.10). In our implementation, we use the two step iterative shrinkage thresholding (TwIST) algorithm BioucasDias:2007tp (); 4378902 (), which combines computational efficiency with strong convergence. To construct the measurement matrix described in eq. (II.9), we must select a grid of possible frequencies () and linewidths (). In our implementation, we use a grid of frequencies ranging from 0 to 2000 cm in 2 cm intervals, and a grid of linewidths ranging from 0 to 160 cm in 6 cm intervals. We assume that our calculations are converged when (in eq. (II.10)), or the solution vector remains constant for 100 iterations. Finally, we perform an minimization of while freezing the recovered nonzero basis functions, allowing us to further minimize the error. We refer to this procedure as debiasing because it partly removes the bias towards sparsity and smoothness introduced by the minimization. This debiasing procedure reduces our solution tolerance to , allowing convergence to a better solution. It is important to note that, in general, the super-resolution technique is robust to an over-complete basis.

## Iv Results

Fig. 1 shows the results of employing the Drude-Lorentz super-resolution method to recover the spectral density for site 1 of FMO. The figure compares Drude-Lorentz super-resolution with 10 ps of MD to a standard fast Fourier transform approach with both 10 and 40 ps of MD. We take the fast Fourier transform with 40 ps of MD as our standard for comparison. By comparing the two methods with 10 ps of MD, it is clear that super-resolution resolves more features of the spectral density than the standard fast Fourier transform from the same amount of time-domain data. Moreover, super-resolution captures most of the features of the fast Fourier transform with the full 40 ps of MD: we see the expected CO stretch at 1600 cm, which we attribute to the amides in the protein scaffold, as well as all of the other major peaks in the spectral density. We attribute a significant amount of the error in our spectral density reconstruction to the fact that the truncated MD series does not explore the phase space as thoroughly in only 10 ps.

The Drude-Lorentz basis also provides significant sparsity gains in comparison to the cosine basis: we require only 56 Drude-Lorentz peaks to create the spectral density given in Fig. 1 whereas. This sparsity provides a significant computational advantage for excitonic propagation in both hierarchical equations of motion (HEOM) JPSJ.58.101 () and second order time-convolutionless master equation Breuer:2007wp ()(TCL-2) approaches because the propagations scale factorially and linearly, respectively, as a function of the number of peaks included. In the excitonically accessible regime of 0-540cm, we recover only 20 Drude-Lorentz peaks, and six of them have amplitudes that are two orders of magnitude smaller than the rest. These Drude-Lorentz peaks can be entered directly into master equation simulations, including HEOM codes, without the need to perform any intermediate fitting 16106 (). In summary, super-resolution yields a well-resolved spectral density using less time-domain data than is required by the standard fast Fourier transform approach and precludes the need for additional fitting.

As mentioned above, the TCL-2 propagation of the exciton dynamics of the FMO complex, with the Hamiltonian coming from Adolphs:2006ey (), was carried out using the Drude-Lorenz spectral densities obtained from super-resolution. We propagated 1 ps of dynamics and obtained the populations of sites 1-3, as well as the coherence between sites 1 and 3.

Fig. 2 shows the coherence between excitonic eigenstates 1 and 3 as a function of time. Compared to the 40 ps fast Fourier transform, we see that the 10 ps Drude-Lorentz super-resolution more faithfully reproduces the coherence dynamics than the 10 ps fast Fourier transform, both in terms of the oscillation frequency and the overall damping. The fast Fourier transform with 10 ps of MD data introduces serious overdamping as well as a significant shift in oscillation frequency. In contrast, the Drude-Lorentz expansion with 10 ps of MD data introduces only a small shift in oscillation frequency, resulting in more accurate coherence dynamics overall. We attribute most of the discrepancies to slight relative differences in the reorganization of each site between spectral densities constructed with 10 and 40 ps of MD data. It appears that while the oscillations are extremely sensitive to the relative reorganization energies between the sites, the damping is more dependent on the fine structure of the spectral densities. The Drude-Lorentz super-resolution (10 ps MD-TDDFT data) reproduces the coherence life-times obtained by fast Fourier transform recovered using all 40 ps of MD-TDDFT data – representing a factor of four improvement.

The contrast between the two approximation techniques becomes even more significant when we simulate dynamics beginning with an exciton fully localized on site 1. In Fig. 3, we have plotted the populations of the first three sites as a function of time. The Drude-Lorentz expansion with 10 ps of MD yields good qualitative agreement with our standard of comparison whereas the fast Fourier transform. The fast Fourier transform on 10 ps overestimates population transfer to site 3 at short times and grows much more quickly from there, whereas the Drude-Lorentz expansion slightly under predicts the population transfer at long times. We attribute these errors in the asymptotic behavior to slight differences in the reorganization energies for the spectral densities of each of the sites, since each site is embedded in a different enviornment, the reorganization process of the individual pigments is different. This sensitivity affects overall dissipation and even small changes in the spectral density of the Drude-Lorentz expansion (10 ps) when compared to the standard of comparison affects energy relaxation. Beyond that, the Drude-Lorentz expansion is capable of reproducing the oscillations at 0.2 and 0.4 ps in the data for sites 1 and 2 whereas the fast Fourier transform reproduces them less faithfully. In summary, the Drude-Lorentz super-resolution technique provides us with much more physical behavior.

## V Conclusions

We have shown that the Drude-Lorentz super-resolution method provides significant computational advantages for the construction of atomistic bath models. In particular, the super-resolution calculations require only 10 ps of MD-TDDFT simulations to obtain reasonable atomistic spectral densities and system dynamics; this is one quarter the amount of data needed in standard fast Fourier transform-based calculations. Ultimately, this will permit the use of more physically accurate calculations or larger systems. Given the computational expense of running TDDFT calculations at every MD simulation step, we believe that the super-resolution method will enable the treatment of larger systems than previously possible.

One of the most significant advantages of our super-resolution method is the decomposition of these atomistic spectral densities into a naturally-sparse basis of Drude-Lorentz oscillators. This makes it easy to perform fast master equation simulations within either the TCL-2 or HEOM formalisms by exploiting analytic integrals of the spectral density. Beyond this, we also directly extract physically-important parameters such as the coherence lifetimes of all the oscillators in the bath. In the future, it is easy to imagine turning this technique on its head to create new spectral densities in a constructive fashion from a set of Drude-Lorentz oscillators.

## Vi Acknowledgements

We acknowledge S. Valleau for useful discussions and computer code. We acknowledge the financial support of Defense Advanced Research Projects Agency grant N66001-10-1-4063 and the Defense Threat Reduction Agency under contract no. HDTRA1-10-1-0046. T.M. acknowledges support from the National Science Foundation (NSF) through the Graduate Research Fellowship Program (GRFP). S.B. acknowledges support from the Department of Energy (DoE) through the Computational Sciences Graduate Fellowship (CSGF). J.N.S. acknowledges support from the Department of Defense (DoD) through the National Defense Science & Engineering Graduate Fellowship (NDSEG) Program. A.A.G. thanks the Corning Foundation.

## References

- (1) H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press, USA, 2007.
- (2) N. Shenvi, J. R. Schmidt, S. T. Edwards, and J. C. Tully, Phys. Rev. A 78, 022502 (2008).
- (3) T. C. Berkelbach, M. S. Hybertsen, and D. R. Reichman, Journal of Chemical Physics 138, 4102 (2013).
- (4) T. C. Berkelbach, M. S. Hybertsen, and D. R. Reichman, J. Chem. Phys. 138, 114103 (2013).
- (5) N. Singh and P. Brumer, Faraday Discuss. 153, 41 (2011).
- (6) S. Jang, M. D. Newton, and R. J. Silbey, J. Phys. Chem. B 111, 6807 (2007).
- (7) M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru-Guzik, J. Chem. Phys. 129, 174106 (2008).
- (8) M. B. Plenio and S. F. Huelga, New J. Phys. 10, 113019 (2008).
- (9) P. Rebentrost, M. Mohseni, I. Kassal, S. Lloyd, and A. Aspuru-Guzik, New J. Phys. 11, 033003 (2009).
- (10) J. Cao and R. J. Silbey, J. Phys. Chem. A 113, 13825 (2009).
- (11) A. Ishizaki and G. R. Fleming, Proc. Natl. Acad. Sci. 106, 17255 (2009).
- (12) A. Ishizaki and G. R. Fleming, J. Phys. Chem. B 115, 6227 (2011).
- (13) M. Sarovar, A. Ishizaki, G. R. Fleming, and K. B. Whaley, Nature Phys. 6, 462 (2010).
- (14) D. Abramavicius and S. Mukamel, J. Chem. Phys. 133, 064510 (2010).
- (15) J. Wu, F. Liu, Y. Shen, J. Cao, and R. J. Silbey, New J. Phys. 12, 105012 (2010).
- (16) J. Moix, J. Wu, P. Huo, D. Coker, and J. Cao, J. Phys. Chem. Lett. 2, 3045 (2011).
- (17) C. Kreisbeck, T. Kramer, M. Rodríguez, and B. Hein, J. Chem. Theory Comput. 7, 2166 (2011).
- (18) N. Skochdopole and D. A. Mazziotti, J. Phys. Chem. Lett. 2, 2989 (2011).
- (19) G. Ritschel, J. Roden, W. T. Strunz, A. Aspuru-Guzik, and A. Eisfeld, J. Phys. Chem. Lett. 2, 2912 (2011).
- (20) P. Rebentrost and A. Aspuru-Guzik, J. Chem. Phys. 134, 101103 (2011).
- (21) L. A. Pachón and P. Brumer, Phys. Chem. Chem. Phys. 14, 10094 (2012).
- (22) S. M. Vlaming and R. J. Silbey, J. Chem. Phys. 136, 055102 (2012).
- (23) F. Caruso, S. K. Saikin, E. Solano, S. F. Huelga, A. Aspuru-Guzik, and M. B. Plenio, Phys. Rev. B 85, 125424 (2012).
- (24) J. Zhu, S. Kais, P. Rebentrost, and A. Aspuru-Guzik, J. Phys. Chem. B 115, 1531 (2011).
- (25) J. Roden, A. Eisfeld, W. Wolff, and W. Strunz, Phys. Rev. Lett. 103, 058301 (2009).
- (26) C. Olbrich and U. Kleinekathöfer, J. Phys. Chem. B 114, 12427 (2010).
- (27) C. Olbrich, T. L. C. Jansen, J. Liebers, M. Aghtar, J. Strümpfer, K. Schulten, J. Knoester, and U. Kleinekathöfer, J. Phys. Chem. B 115, 8609 (2011).
- (28) B. Hein, C. Kreisbeck, T. Kramer, and M. Rodríguez, New J. Phys. 14, 023018 (2012).
- (29) N. Christensson, H. F. Kauffmann, T. Pullerits, and T. Mančal, J. Phys. Chem. B 116, 7449 (2012).
- (30) A. W. Chin, J. Prior, R. Rosenbach, F. Caycedo-Soler, S. F. Huelga, and M. B. Plenio, Nature Phys. 9, 113 (2013).
- (31) C. Kreisbeck and T. Kramer, J. Phys. Chem. Lett. 3, 2828 (2012).
- (32) S. Valleau, A. Eisfeld, and A. Aspuru-Guzik, J. Chem. Phys. 137, 224103 (2012).
- (33) M. Tuckerman, Statistical Mechanics: Theory and Molecular Simulation, OUP Oxford, 2010.
- (34) E. Runge and E. K. Gross, Phys. Rev. Lett. 52, 997 (1984).
- (35) S. Mallat, A Wavelet Tour of Signal Processing, Third Edition: The Sparse Way, Academic Press, 3 edition, 2008.
- (36) W. T. Freeman, T. R. Jones, and E. C. Pasztor, IEEE Comput. Graph. Appl. 22, 56 (2002).
- (37) A. J. Patti, M. I. Sezan, and A. Murat Tekalp, IEEE Trans. Image Process. 6, 1064 (1997).
- (38) M. Elad and A. Feuer, IEEE Trans. Image Process. 6, 1646 (1997).
- (39) K. G. Puschmann and F. Kneer, Astron. Astrophysic. 436, 373 (2005).
- (40) C. W. Mccutchen, J. Opt. Soc. Am. 57, 1190 (1967).
- (41) D. Kouame and M. Ploquin, Super-resolution in medical imaging: An illustrative approach through ultrasound, in Biomedical Imaging: From Nano to Macro, 2009. ISBI ’09. IEEE International Symposium on, pp. 249–252, 2009.
- (42) J. Ma, Instrumentation and Measurement, IEEE Transactions on 59, 1600 (2010).
- (43) D. L. Donoho, IEEE Trans. Inform. Theory 52, 1289 (2006).
- (44) A. Oka and L. Lampe, A compressed sensing receiver for bursty communication with UWB Impulse Radio, in Ultra-Wideband, 2009. ICUWB 2009. IEEE International Conference on, pp. 279–284, 2009.
- (45) M. A. Herman and T. Strohmer, Signal Processing, IEEE Transactions on 57, 2275 (2009).
- (46) C. Qiu, W. Lu, and N. Vaswani, Real-time dynamic MR image reconstruction using Kalman Filtered Compressed Sensing, in Acoustics, Speech and Signal Processing, 2009. ICASSP 2009. IEEE International Conference on, pp. 393–396, 2009.
- (47) M. Lustig, D. Donoho, and J. M. Pauly, Magnetic Resonance in Medicine 58, 1182 (2007).
- (48) M. Nagahara, D. E. Quevedo, and J. Ostergaard, Sparse representations for packetized predictive networked control, in Decision and Control (CDC), 2012 IEEE 51st Annual Conference on, pp. 1362–1367, 2012.
- (49) T. Tuma, S. Rooney, and P. Hurley, On the Applicability of Compressive Sampling in Fine Grained Processor Performance Monitoring, in 2009 14th IEEE International Conference on Engineering of Complex Computer Systems, pp. 210–219, IEEE, 2009.
- (50) A. Shabani, R. L. Kosut, M. Mohseni, H. Rabitz, M. A. Broome, M. P. Almeida, A. Fedrizzi, and A. G. White, arXiv (2009).
- (51) M. Mishali and Y. C. Eldar, Selected Topics in Signal Processing, IEEE Journal of 4, 375 (2010).
- (52) M. F. Duarte, M. A. Davenport, D. Takhar, J. N. Laska, T. Sun, K. F. Kelly, and R. G. Baraniuk, Signal Processing Magazine, IEEE 25, 83 (2008).
- (53) W. K. Coulter, C. J. Hillar, G. Isley, and F. T. Sommer, p. 5494 (2010).
- (54) E. J. Candès and C. Fernandez-Granda, Comm. Pure Appl. Math. (2013).
- (55) J. M. Bioucas-Dias and M. A. Figueiredo, IEEE Trans. Image Process. 16, 2992 (2007).
- (56) J. M. Bioucas-Dias and M. A. T. Figueiredo, Two-step algorithms for linear inverse problems with non-quadratic regularization, in Image Processing, 2007. ICIP 2007. IEEE International Conference on, pp. I–105–I – 108, 2007.
- (57) S. I. E. Vulto, M. A. de Baat, R. J. W. Louwe, H. P. Permentier, T. Neef, M. Miller, H. van Amerongen, and T. J. Aartsma, J. Phys. Chem. B 102, 9577 (1998).
- (58) J. Adolphs and T. Renger, Biophys. J. 91, 2778 (2006).
- (59) M. Mohseni, A. Shabani, S. Lloyd, and H. Rabitz, arXiv (2011).
- (60) X. Chen, J. Cao, and R. J. Silbey, J. Chem. Phys. (2013).
- (61) J. Yuen-Zhou, J. J. Krich, and A. Aspuru-Guzik, J. Chem. Phys. 136, 234501 (2012).
- (62) S. Shim, P. Rebentrost, S. p. Valleau, and A. n Aspuru-Guzik, Biophys. J. 102, 649 (2012).
- (63) A. Kolli, A. Nazir, and A. Olaya-Castro, J. Chem. Phys. 135, 154112 (2011).
- (64) S. Jang, J. Chem. Phys. 135, 034105 (2011).
- (65) A. Ishizaki and G. R. Fleming, J. Chem. Phys. 130, 234111 (2009).
- (66) Y. Tanimura and R. Kubo, Journal of the Physical Society of Japan 58, 101 (1989).
- (67) C. Kreisbeck and T. Kramer, Exciton Dynamics Lab for Light-Harvesting Complexes (GPU-HEOM), 2013.
- (68) P. Rebentrost and A. Aspuru-Guzik, J. Chem. Phys. 134, 101103 (2011).
- (69) A. Pereverzev and E. R. Bittner, Journal of Chemical Physics 125, 4906 (2006).
- (70) C. Timm, Phys. Rev. B 83, 115416 (2011).
- (71) D. Ahn, Phys. Rev. B 50, 8310 (1994).
- (72) B. Heinz-Peter, B. Kappler, and F. Petruccione, Decoherence: Theoretical, Experimental, and Conceptual Problems , 233 (2000).
- (73) A. Shabani and D. A. Lidar, Phys. Rev. A 71, 020101 (2005).
- (74) A. Smirne and B. Vacchini, Phys. Rev. A 82, 022110 (2010).
- (75) A. Kolli, A. Nazir, and A. Olaya-Castro, J. Chem. Phys. 135, 154112 (2011).
- (76) U. Kleinekathöfer, J. Chem. Phys. 121, 2505 (2004).
- (77) P. H. Berens, S. R. White, and K. R. Wilson, J. Chem. Phys. 75, 515 (1981).
- (78) E. J. Candes, J. Romberg, and T. Tao, IEEE Trans. Inform. Theory 52, 489 (2006).
- (79) X. Andrade, J. N. Sanders, and A. Aspuru-Guzik, Proc. Natl. Acad. Sci. 109, 13928 (2012).
- (80) J. N. Sanders, S. K. Saikin, S. Mostame, X. Andrade, J. R. Widom, A. H. Marcus, and A. Aspuru-Guzik, J. Phys. Chem. Lett. 3, 2697 (2012).
- (81) W. D. Cornell, P. Cieplak, C. I. Bayly, I. R. Gould, K. M. Merz, D. M. Ferguson, D. C. Spellmeyer, T. Fox, J. W. Caldwell, and P. A. Kollman, J. Am. Chem. Soc. 117, 5179 (1995).
- (82) M. Ceccarelli, P. Procacci, and M. Marchi, J. Comput. Chem. 24, 129 (2003).
- (83) A. D. Becke, Phys. Rev. A 38, 3098 (1988).
- (84) B. Miehlich, A. Savin, H. Stoll, and H. Preuss, Chem. Phys. Lett. 157, 200 (1989).
- (85) C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 (1988).
- (86) Y. Shao, L. F. Molnar, Y. Jung, J. Kussmann, C. Ochsenfeld, S. T. Brown, A. T. Gilbert, L. V. Slipchenko, S. V. Levchenko, and D. P. O’Neill, Phys. Chem. Chem. Phys. 8, 3172 (2006).