QuantumCumulants.jl: A Julia framework for generalized mean-field equations in open quantum systems

A full quantum mechanical treatment of open quantum systems via a Master equation is often limited by the size of the underlying Hilbert space. As an alternative, the dynamics can also be formulated in terms of systems of coupled differential equations for operators in the Heisenberg picture. This typically leads to an infinite hierarchy of equations for products of operators. A well-established approach to truncate this infinite set at the level of expectation values is to neglect quantum correlations of high order. This is systematically realized with a so-called cumulant expansion, which decomposes expectation values of operator products into products of a given lower order, leading to a closed set of equations. Here we present an open-source framework that fully automizes this approach: first, the equations of motion of operators up to a desired order are derived symbolically using predefined canonical commutation relations. Next, the resulting equations for the expectation values are expanded employing the cumulant expansion approach, where moments up to a chosen order specified by the user are included. Finally, a numerical solution can be directly obtained from the symbolic equations. After reviewing the theory we present the framework and showcase its usefulness in a few example problems.

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 72 is April 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Rapid Microscopic Fractional Anisotropy Imaging via an Optimized Kurtosis Formulation

Water diffusion anisotropy in the human brain is affected by disease, trauma, and development. Microscopic fractional anisotropy (μFA) is a diffusion MRI (dMRI) metric that can quantify water diffusion anisotropy independent of neuron fiber orientation dispersion. However, there are several different techniques to estimate μFA and few have demonstrated full brain imaging capabilities within clinically viable scan times and resolutions. Here, we present an optimized spherical tensor encoding (STE) technique to acquire μFA directly from the 2nd order cumulant expansion of the dMRI signal (i.e. diffusion kurtosis) which requires fewer powder-averaged signals than other STE fitting techniques and can be rapidly computed. We found that the optimal dMRI parameters for white matter μFA imaging were a maximum b-value of 2000 s/mm2 and a ratio of isotropic to linear tensor encoded acquisitions of 1.7 for our system specifications. We then compared two implementations of the direct approach to the well-established gamma model in 4 healthy volunteers on a 3 Tesla system. One implementation of the direct cumulant approach used mean diffusivity (D) obtained from a 2nd order fit of the cumulant expansion, while the other used a linear estimation of D from the low b-values. Both implementations of the direct approach showed strong linear correlations with the gamma model (ρ=0.97 and ρ=0.90) but mean biases of −0.11 and −0.02 relative to the gamma model were also observed, respectively. All three μFA measurements showed good test-retest reliability (ρ≥0.79 and bias=0). To demonstrate the potential scan time advantage of the direct approach, 2 mm isotropic resolution μFA was demonstrated over a 10 cm slab using a subsampled data set with fewer powder-averaged signals that would correspond to a 3.3-minute scan. Accordingly, our results introduce an optimization procedure that has enabled clinically relevant, nearly full brain μFA in only several minutes.HighlightsDemonstrated method to acquire optimal parameters for regression μFA imagingμFA measured using an optimized linear regression method at 3TFirst μFA comparison between direct regression approach and the gamma modelBoth approaches correlated strongly in white matter in healthy volunteersNearly full brain μFA demonstrated in a 3.3-minute scan at 2 mm isotropic resolution

Absorption and Circular Dichroism Spectra of Molecular Aggregates with the Full Cumulant Expansion.

The exciton Hamiltonian of multichromophoric aggregates can be probed by spectroscopic techniques such as linear absorption and circular dichroism. In order to compare calculated Hamiltonians to experiments, a lineshape theory is needed, which takes into account the coupling of the excitons with inter- and intramolecular vibrations. This coupling is normally introduced in a perturbative way through the cumulant expansion formalism, and further approximated by assuming a Markovian exciton dynamics, for example with the modified Redfield theory.<br><br>Here we present an implementation of the full cumulant expansion (FCE) formalism [Ma and Cao, <i>J. Chem. Phys.</i> <b>2015</b>, 142, 094106 ] to efficiently compute absorption and circular dichroism spectra of molecular aggregates beyond the Markov approximation, without restrictions on the form of the exciton-phonon coupling. By employing the LH2 system of purple bacteria as a challenging test case, we compare the FCE lineshapes with the Markovian lineshapes obtained with the modified Redfield theory, showing that the latter present a much poorer agreement with experiments. The FCE approach instead accurately describes the lineshapes, especially in the vibronic sideband of the B800 peak. We envision that the FCE approach will become a valuable tool for accurately comparing model exciton Hamiltonians with optical spectroscopy experiments.

Absorption and Circular Dichroism Spectra of Molecular Aggregates with the Full Cumulant Expansion.

The exciton Hamiltonian of multichromophoric aggregates can be probed by spectroscopic techniques such as linear absorption and circular dichroism. In order to compare calculated Hamiltonians to experiments, a lineshape theory is needed, which takes into account the coupling of the excitons with inter- and intramolecular vibrations. This coupling is normally introduced in a perturbative way through the cumulant expansion formalism, and further approximated by assuming a Markovian exciton dynamics, for example with the modified Redfield theory.<br><br>Here we present an implementation of the full cumulant expansion (FCE) formalism [Ma and Cao, <i>J. Chem. Phys.</i> <b>2015</b>, 142, 094106 ] to efficiently compute absorption and circular dichroism spectra of molecular aggregates beyond the Markov approximation, without restrictions on the form of the exciton-phonon coupling. By employing the LH2 system of purple bacteria as a challenging test case, we compare the FCE lineshapes with the Markovian lineshapes obtained with the modified Redfield theory, showing that the latter present a much poorer agreement with experiments. The FCE approach instead accurately describes the lineshapes, especially in the vibronic sideband of the B800 peak. We envision that the FCE approach will become a valuable tool for accurately comparing model exciton Hamiltonians with optical spectroscopy experiments.

Nonlinear Spectroscopy in the Condensed Phase: The Role of Duschinsky Rotations and Third Order Cumulant Contributions

First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with a specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and Methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.