Terahertz pulse time-domain holography (THz PTDH) is the powerful technique for high-resolution amplitude and phase THz imaging that allows mapping spectroscopic information across the imaged object. In this paper, we consider most sought after applications of phase imaging provided by this technique and experimentally demonstrate the ability of the method to reconstruct smooth and stepped relief features of an object that is transparent in THz region. Unlike the amplitude distribution, which does not contain any significant information in this case, phase distribution not only reveals the object qualitatively, but also allows the reconstruction of the object thicknesses pattern, even in low signal-to-noise registration conditions. Main limitations of the proposed method, such as transverse resolution and low signal-to-noise environment are carefully studied and mitigated.

In this paper, we consider computational super-resolution inverse diffraction phase retrieval. The optical setup is lensless, with a spatial light modulator for aperture phase coding. The paper is focused on experimental tests of the super-resolution sparse phase amplitude retrieval algorithm. We start from simulations and proceed to physical experiments. Both simulation tests and experiments demonstrate good-quality imaging for super-resolution with a factor of 4 and a serious advantage over diffraction-limited resolution as defined by Abbe’s criterion.

Terahertz pulse time-domain holography is the ultimate technique allowing the evaluating a propagation of pulse broadband terahertz wavefronts and analyze their spatial, temporal and spectral evolution. We have numerically analyzed pulsed broadband terahertz Gauss-Bessel beam’s both spatio-temporal and spatio-spectral evolution in the non-paraxial approach. We have characterized two-dimensional spatio-temporal beam behavior and demonstrated all stages of pulse reshaping during the propagation, including X-shape pulse forming. The reshaping is also illustrated by the energy transfer dynamics, where the pulse energy flows from leading edge to trailing edge. This behavior illustrates strong spatio-temporal coupling effect when spatio-temporal distribution of Bessel beam’s wavefront depends on propagation distance. The spatio-temporal and spatio-spectral profiles for different spectral components clearly illustrate the model where the Bessel beam’s wavefront at the exit from the axicon can be divided into radial segments for which the wave vectors intersect. Phase velocity via propagation distance is estimated and compared with existing experimantal results. Results of the phase velocity calculation depend strongly on distance increment value, thus demonstrating superluminal or subluminal behavior.

Terahertz phase retrieval is a promising technique able to assess the complex diffracted wave properties through an iterative processing algorithm. In this Letter, we demonstrate the implementation of this technique in reflection geometry with a continuous wave acquisition system working at 0.287 THz. To ensure a high signal-to-noise ratio in the measured dataset, we proposed a double parallel recording scheme with one detector and two lock-in amplifiers operating with the complimentary sensitivity setting. This provided a higher numerical aperture than conventional raster-scanning focal plane imaging. A specialized digital interferometric postprocessing procedure was applied to obtain a surface height map from the reconstructed phase distribution in the object’s irradiated area.

In this paper, we present a novel numerical approach for increasing the resolution of retrieved images of objects after their diffraction patterns are recorded via terahertz pulse time-domain holography (THz PTDH). THz PTDH allows for spectrally resolved imaging with high spatial resolution and does not require the fine alignment of complex optics in the THz path. The proposed data post-processing method opens up the possibility to reconstruct holograms recorded with spatially restricted THz detectors, and overcome the diffraction limit even for the lower-frequency spectral components. The method involves an iterative procedure of backward-forward wavefront propagation to simulate the field distribution beyond the initially recorded hologram area. We show significant improvement in both the object reconstruction and contrast across the whole spectrum, with qualitative resolution enhancement at lower frequency spectral components.