Generating wavelet using amplitude spectra Equation (4) of a recent paper (An improved peak-frequency-shift method for Q estimation) shows how to generate a wavelet using a zero-phase Ricker wavelet.

$|B_1(f)| = |B_0(f)|\exp(-\pi t f / Q)$

In the above, $|B_1(f)|$ is the amplitude spectra of the wavelet that I want to calculate, $|B_0(f)|$ is the amplitude spectra of the reference (Ricker wavelet), $\pi$ is the ubiquitous constant, $t$ is the time separation between wavelets, $f$ is the frequency, and $Q$ is a number (the Quality factor).

I know how to calculate the zero-phase Ricker wavelet in the time domain as $B_0(t)$, and I know how to calculate $|B_0(f)|$ as the square root of the sum of the real and imaginary parts squared of the frequency-domain reference Ricker wavelet.  The real and imaginary parts are found using the FFT.

However, now that I've calculated $|B_1(f)|$, how do I transform $|B_1(f)|$ back into the 
time domain?  I know $|B_1(f)|$, but can I find the real and imaginary parts of $|B_1(f)|$ so that I can use the IFFT to compute $B_1(t)$ in the time domain, or is some other information missing and required?

Is the "without phase changes" words in the recent paper (http://www.seg.org/documents/10161/997244/1218.pdf) a clue that might help with the numerical implementation?

So how did the authors generate the wavelets in Figure 1 of the paper?  Is something other than Equation (4) being used here? 

Perhaps there is nothing more happening here than an issue with notation. Looking at another paper (Estimation of quality factors from CMP records), the frequency spectrum of a Ricker wavelet is computed using Equation (1):

$B(f) = \frac{2}{{\sqrt \pi  }}\frac{{{f^2}}}{{f_m^2}}{e^{\frac{{ - {f^2}}}{{f_m^2}}}}$

In the above, $f_m$ is the dominant frequency.  However, there doesn't seem to be any complex part of this expression.  Does this imply that the complex part is comprised of all zeros?  To numerically evaluate this equation, how do I compute the $f$ variable?  For example, using Matlab, $f$ can be computed using the colon operator:

f = fmin:deltaf:fmax


But what is fmin and fmax?  Now $B(f)$ can be computed and taken back into the time domain using the IDFT (or the IFFT), but what if I want a real time-domain discrete signal for the $B(t)$ function?

Equation (5) of the same paper gives an expression that I think can be used to calculate the Q factor wavelet (note that this looks the same as the equation at the very top of this post):

$B(f,t) = B(f){e^{\frac{{ - \pi ft}}{Q}}}$

However, what is meant here by the "amplitude spectrum"?  How do I get the real and complex parts of this spectrum?

A time domain version equation of the Ricker wavelet can be found here: Wikipedia link, but then what is the time domain expression used to calculate the Q factor wavelet?  As an aside, it would have been easier to have the time domain version of:

$B(f,t) = B(f){e^{\frac{{ - \pi ft}}{Q}}}$

But what is the "amplitude spectrum" in this context, and how can I compute the wavelet in the time domain?
The complex result of a DFT of a signal codifies amplitude and phase information. When you take a magnitude of the resulting complex vector, you automatically throw away any phase information.

In so many words, the DFT result is informing you of what amplitudes and phases of all its frequency components are required to reconstruct your signal. If you strip away phase information you are lining up all your complex sinusoids and then adding them up coherently. If you lose phase how can you reconstruct your signal? 

That being said, what you can get in this case is the autocorrelation function of $b_1(t)$, by taking the IDFT of $|B_1(f)|^2$, but you need phase information to completely reconstruct the original signal.