I don't believe this is possible. Before you calculate the power spectral density, you must first have a power spectrum. When you convert from a FFT to a power spectrum, your phase information is lost. With magnitude information only, it isn't possible to go back to the time domain from the frequency domain.
One solution would be to use the PXI-5620 frequency domain digitizer with the included Spectral Measurements Toolset. The PXI-5620 is a 14-bit, 64 MSample/sec digitizer which you could use to digitize the raw time domain data. The Spectral Measurements Toolset can be used to turn this data into FFTs, power spectrums, and power spectral density traces, in addition to performing several types of spectral measurements upon them. In addition, since you would
be acquiring the data with a digitizer, you would have access to all the raw time data and could make any number of time domain measurements with built in LabVIEW functionality. With 32 MHz of Nyquist bandwidth on the digitizer (one half of 64 MSample/sec), you could capture all the odd harmonics of the 1 MHz pulse up to 31 MHz which should be enough to capture the pulse adequately.