Computational techniques for image reconstruction in terahertz microscopy using Josephson cantilevers

The terahertz microscope is a measurement system designed primarily for imaging of spatial distributions of high-frequency electromagnetic fields in the microwave and terahertz bands. It relies on a Josephson junction as a sensor, which makes it possible to demodulate the signal for low frequency readout while also retaining information about its original frequencies. The junction is located on the tip of a cantilever that can be brought into the near field of the examined sample, allowing the microscope to achieve higher image resolutions then what would be otherwise possible due to the Rayleigh criterion at these wavelengths. This work focuses on the algorithms designed for reconstruction of the irradiating spectra and the images from the sensor’s responses, as well as discussing the electronics required for the measurement itself.

The data is acquired through sampling of the current-voltage characteristics of the junction for each of the target pixels individually. After a brief overview of the terahertz microscopy setup, the different Josephson junction readout schemes used throughout the literature are discussed and the solution designed for this setup is presented. Due to the short measurement times required for imaging, there are tradeoffs to be made in terms of the accuracy and noise of this readout. Therefore, a set of preprocessing steps is employed, including automated offset correction and noise-tolerant reconstruction of the differential resistance of the characteristic, latter required by other processing algorithms.

For small signal powers, spectra can be reconstructed based on closed-form approximations of the small-signal response of the junction. Alternative to the well-known Hilbert spectroscopy, a direct deconvolution method has proven to be well suited for the employed measurement scheme. For larger signal amplitudes, the expected response can be calculated using a number of numerical methods. To facilitate curve fitting for a large number of pixel measurements based on these methods, graphics processing units (GPUs) are employed to accelerate such computations through parallelization.

Finally, a set of exemplary measurements is presented, highlighting the cases where the proposed reconstruction techniques prove superior to the ones used in previous works Das Terahertz-Mikroskop ist ein Messsystem, das in erster Linie für die Abbildung der räumlichen Verteilung hochfrequenter elektromagnetischer Felder im Mikrowellen- und Terahertz-Bereich konzipiert wurde. Dabei kommen Josephson-Kontakte als Sensoren zum Einsatz, die es ermöglichen, das Signal zum niederfrequenten Auslesen zu demodulieren und dabei die Informationen über seine ursprünglichen Frequenzen beizubehalten. Der Kontakt befindet sich auf der Spitze eines Cantilevers, die in das Nahfeld der Probe gebracht werden kann, was eine höhere Bildauflösung, als sonst aufgrund des Rayleigh-Kriteriums möglich wäre, ermöglicht. Diese Arbeit beschäftigt sich mit der Entwicklung von Algorithmen zur Rekonstruktion der Bestrahlungsspektren und der Bilder aus den Messdaten, sowie mit der für deren Aufnahme erforderlichen Elektronik.

Bei den Messungen werden in der Regel die gesamten Strom-Spannungs-Kennlinien des Kontaktes für jeden einzelnen Bildpixel aufgenommen. Nach einem kurzen Überblick über den Aufbau des gesamten Terahertz-Mikroskops werden die üblichen Verfahren zum Auslesen von Josephson-Kontakten diskutiert und die für diesen Aufbau entwickelte Lösung vorgestellt. Aufgrund der kurzen Messzeiten, die für die Bildgebung erforderlich sind, müssen Kompromisse in Bezug auf die Genauigkeit und das Rauschen der Messdaten eingegangen werden. Daher wird eine Reihe von Vorverarbeitungsschritten unternommen, darunter automatische Offsetkorrektur sowie rauschtolerante Rekonstruktion des differentiellen Widerstandes der Kennlinie, welches für weitere Verarbeitungsalgorithmen benötigt wird.

Für kleine Signalstärken können die aufgenommenen Spektren auf Grundlage analytischer Näherungen der Kleinsignalantwort rekonstruiert werden. Als Alternative zur bekannten Hilbert-Spektroskopie hat sich eine direkte Entfaltungsmethode für den verwendeten Messaufbau als vorteilhaft erwiesen. Für größere Amplituden muss die zu erwartende Antwort mit numerischen Verfahren berechnet werden. Es wird gezeigt, wie Grafikprozessoren eingesetzt werden können, um eine schnellere Kurvenanpassung für eine große Anzahl von Pixelmessungen basierend auf solchen Berechnungen zu ermöglichen.

Abschließend wird eine Reihe von Beispielmessungen vorgestellt, die die Vorteile der entwickelten Bildrekonstruktionsverfahren gegenüber den in den früheren Arbeiten verwendeten Heuristiken aufzeigen.