Optimizing sparse sampling for 2D electronic spectroscopy

Fisher information for smart sampling in time-domain spectroscopy

Time-domain spectroscopy encompasses a wide range of techniques, such as Fourier-transform infrared, pump-probe, Fourier-transform Raman, and two-dimensional electronic spectroscopies. These methods enable various applications, such as molecule characterization, excited state dynamics studies, or spectral classification. Typically, these techniques rarely use sampling schemes that exploit the prior knowledge scientists typically have before the actual experiment. Indeed, not all sampling coordinates carry the same amount of information, and a careful selection of the sampling points may notably affect the resulting performance. In this work, we rationalize, with examples, the various advantages of using an optimal sampling scheme tailored to the specific experimental characteristics and/or expected results. We show that using a sampling scheme optimizing the Fisher information minimizes the variance of the desired parameters. This can greatly improve, for example, spectral classifications and multidimensional spectroscopy. We demonstrate how smart sampling may reduce the acquisition time of an experiment by one to two orders of magnitude, while still providing a similar level of information.

A Quantum Algorithm from Response Theory: Digital Quantum Simulation of Two-Dimensional Electronic Spectroscopy

Multidimensional optical spectroscopies are powerful techniques to investigate energy transfer pathways in natural and artificial systems. Because of the high information content of the spectra, numerical simulations of the optical response are of primary importance to assist the interpretation of spectral features. However, the increasing complexity of the investigated systems and their quantum dynamics call for the development of novel simulation strategies. In this work, we consider using digital quantum computers. By combining quantum dynamical simulation and nonlinear response theory, we present a quantum algorithm for computing the optical response of molecular systems. The quantum advantage stems from the efficient quantum simulation of the dynamics governed by the molecular Hamiltonian, and it is demonstrated by explicitly considering exciton-vibrational coupling. The protocol is tested on a near-term quantum device, providing the digital quantum simulation of the linear and nonlinear response of simple molecular models.

Two-dimensional spectroscopy beyond the perturbative limit: The influence of finite pulses and detection modes

Ultra-fast and multi-dimensional spectroscopy gives a powerful looking glass into the dynamics of molecular systems. In particular, two-dimensional electronic spectroscopy (2DES) provides a probe of coherence and the flow of energy within quantum systems, which is not possible with more conventional techniques. While heterodyne-detected (HD) 2DES is increasingly common, more recently fluorescence-detected (FD) 2DES offers new opportunities, including single-molecule experiments. However, in both techniques, it can be difficult to unambiguously identify the pathways that dominate the signal. Therefore, the use of numerically modeling of 2DES is vitally important, which, in turn, requires approximating the pulsing scheme to some degree. Here, we employ non-perturbative time evolution to investigate the effects of finite pulse width and amplitude on 2DES signals. In doing so, we identify key differences in the response of HD and FD detection schemes, as well as the regions of parameter space where the signal is obscured by unwanted artifacts in either technique. Mapping out parameter space in this way provides a guide to choosing experimental conditions and also shows in which limits the usual theoretical approximations work well and in which limits more sophisticated approaches are required.

Interplay between coherence-time undersampling and scattered light in two-dimensional electronic spectroscopy

Scanning pulse delays in multi-pulse non-linear optical spectroscopy experiments is a major contributor to lengthy data acquisition. Using large steps for the scan can significantly speed up the experiment. However, an improper choice of step length can cause distortions to the resulting spectra, especially if the light scattered on the sample is mixed into the signal. In this work, we identify potential sources of such distortions and suggest appropriate countermeasures to avoid them while maintaining a faster data collection.

Compressed Sensing for Reconstructing Coherent Multidimensional Spectra

We apply two sparse reconstruction techniques, the least absolute shrinkage and selection operator (LASSO) and the sparse exponential mode analysis (SEMA), to two-dimensional (2D) spectroscopy. The algorithms are first tested on model data, showing that both are able to reconstruct the spectra using only a fraction of the data required by the traditional Fourier-based estimator. Through the analysis of the sparsely sampled experimental fluorescence-detected 2D spectra of LH2 complexes, we conclude that both SEMA and LASSO can be used to significantly reduce the required data, still allowing one to reconstruct the multidimensional spectra. Of the two techniques, it is shown that SEMA offers preferable performance, providing more accurate estimation of the spectral line widths and their positions. Furthermore, SEMA allows for off-grid components, enabling the use of a much smaller dictionary than that of the LASSO, thereby improving both the performance and the lowering of the computational complexity for reconstructing coherent multidimensional spectra.

Optimized reconstructions of compressively sampled two-dimensional infrared spectra

Compressive sampling has the potential to dramatically accelerate the pace of data collection in two-dimensional infrared (2D IR) spectroscopy. We have previously introduced the Generic Iteratively Reweighted Annihilating Filter (GIRAF) reconstruction algorithm to solve the reconstruction in 2D IR compressive sampling. Here, we report a thorough assessment of this method and comparison to our earlier efforts using the Total Variation (TV) algorithm. We show that the GIRAF algorithm has some distinct advantages over TV. Although it is no better or worse in terms of ameliorating the impacts of compressive sampling on the measured 2D IR line shape, we find that the nature of those effects is different for GIRAF than they were for TV. In addition to assessing the impacts on the line shape of a single oscillator, we also test the ability of the algorithm to reconstruct spectra that have transitions from more than one oscillator, such as the coupled carbonyl oscillators in rhodium dicarbonyl. Finally, and perhaps most importantly, we show that the GIRAF algorithm has a distinct denoising effect on the signal-to-noise ratio (SNR) of the 2D IR spectra that can increase the SNR by as much as 4× without any additional signal averaging and collecting fewer data points, which should further enhance the acceleration of data collection that can be achieved using compressive sampling and enable even more challenging experimental measurements.

Optimising magnetic resonance sampling patterns for parametric characterisation

Sampling strategies are often central to experimental design. Choosing efficiently which data to acquire can improve the estimation of parameters and reduce the acquisition time. This work is focused on designing optimal sampling patterns for Nuclear Magnetic Resonance (NMR) applications, illustrated with respect to the best estimate of the parameters characterising a lognormal distribution. Lognormal distributions are commonly used as fitting models for distributions of spin-lattice relaxation time constants, spin-spin relaxation time constants and diffusion coefficients. A method for optimising the choice of points to be sampled is presented which is based on the Cramér-Rao Lower Bound (CRLB) theory. The method's capabilities are demonstrated experimentally by applying it to the problem of estimating the emulsion droplet size distribution from a pulsed field gradient (PFG) NMR diffusion experiment. A difference of <5% is observed between the predictions of CRLB theory and the PFG NMR experimental results. It is shown that CLRB theory is stable down to signal-to-noise ratios of ∼10. A sensitivity analysis for the CRLB theory is also performed. The method of optimizing sampling patterns is easily adapted to distributions other than lognormal and to other aspects of experimental design; case studies of optimising the sampling scheme for a fixed acquisition time and determining the potential for reduction in acquisition time for a fixed parameter estimation accuracy are presented. The experimental acquisition time is typically reduced by a factor of 3 using the proposed method compared to a constant gradient increment approach that would usually be used.