Details

To calculate the real part of the one-sided Fourier transform of a function , we first sample from to , where is the time step and is half the number of steps (this is to ensure an array of even length). The resulting sample vector of the one-sided (red dots) is then changed to a two-sided sample vector such that (red line). In doing so, we drop the repeating and elements. Passing this vector to the built-in Mathematica function Fourier yields the desired one-sided Fourier transform multiplied by two.

Details of the theory of the linear absorption and emission spectra of a single excitation coupled (with coupling strength given by the reorganization energy) to a thermal environment (with bath vibrational frequency and temperature) be can be found in [1], where the well-known Drude–Lorentz, Gaussian, and Lorentzian lineshapes are considered.

Reference

[1] S. Mukamel, Principles of Nonlinear Optical Spectroscopy, New York: Oxford University Press, 1995.