Getting audio samples from power spectrum
No, the power spectrum only holds half of the information that the DFT'ed signal held, which held all the information that the original audio signal held.

Sadly, that half is important, and not reconstruction is possible. You might come up with an estimator that "guesses" the phase information you're losing by doing the magnitude square, but it will not recover the original signal in the general case. I'm not aware of any such estimators that do something useful.

In short: you can reverse the DFT, it's a reversible operation. You cannot revert calculating the magnitude square of complex numbers – that is a irreversible operation, and hence, you can't recover the original signal from a PSD estimate. 

To illustrate this: 

Imagine you record one hour (let's say from $t=-18000$ to $t=1800) of perfect silence, and in the middle of that is one second of a pure 100 Hz sine tone. 

The DFT of that is pretty simple: In time domain, it's a sine (which is infinite) multiplied with a unit rectangular window. We know the spectrum of that; it's the convolution of the Fourier transform of the sine (which is dirac with "area" $\frac 12j$ at the positive frequency, and $-\frac12j$ at the negative frequency) convolved with the Fourier transform of a rectangular window (which is always some scaled $\text{sinc}$ function). Hence, we know that the DFT of that recording must be the sum of two shifted sincs, one of those being "mirrored" on the frequency axis.

Now, we estimate the PSD by just calculating the magnitude square point-wise, and we end up with something that looks like two positive lumps at the positive and negative frequencies. For obvious reasons, it's symmetric to the magnitude axis. We call the DFT $X(f)$, and that leaves us with the PSD estimate $|X(F)|^2$

We now decide that the 1s of tone shouldn't be at the $t=0$, but let's say at $t=900$; now, instead of recalculating the whole spectrum, we just apply the time-shift property of Fourier transform: a $\Delta t$ shift in time domain is a multiplication with $e^{j2\pi f \Delta t}$ in frequency domain. That way, we just need to multiply with this function in frequency domain.

And lo! Let $X(f)\cdot e^{j2\pi f \Delta t}$ be our spectrum, and $|X(f)\cdot e^{j2\pi f \Delta t}|^2$ be our PSD estimate. Observe that the modulation with the complex $e$xponential doesn't change magnitude square – hence, you can't tell time shifts from your PSD alone.