A Scanning Tunneling Microscope (STM) is one of the most important scanning probe tools available to study and manipulate matter at the nanoscale. In a STM, a tip is scanned on top of a surface with a separation of a few Å. Often, the tunneling current between the tip and the sample is maintained constant by modifying the distance between the tip apex and the surface through a feedback mechanism acting on a piezoelectric transducer. This produces very detailed images of the electronic properties of the surface. The feedback mechanism is nearly always made using a digital processing circuit separate from the user computer. Here, we discuss another approach using a computer and data acquisition through the universal serial bus port. We find that it allows successful ultralow noise studies of surfaces at cryogenic temperatures. We show results on different compounds including a type II Weyl semimetal (WTe2), a quasi-two-dimensional dichalcogenide superconductor (2H–NbSe2), a magnetic Weyl semimetal (Co3Sn2S2), and an iron pnictide superconductor (FeSe).

Binnig et al. invented the Scanning Tunneling Microscope (STM) back in 1981.1 The technique was soon extended to other probes based on a tip scanning a sample surface, revolutionizing the study and manipulation of materials at the nanoscale. The STM uses the quantum tunneling effect between an atomically sharp tip and a flat conducting sample through a vacuum barrier to sense the surface properties and the piezoelectric effect that provides precise subnanometric positioning of the tip over the sample. Usually, we measure the tunneling current between the tip and the sample and maintain its value constant using a feedback loop that acts on the piezo that controls the z-position of the tip. When scanning the tip in the x–y plane, the feedback signal provides the z-position of the tip. Thus, we can build two-dimensional (2D) maps of the surface at a constant tunneling current.

There are many different hardware designs to control scanning probe microscopes.2–6 All STM designs are based on a kernel that includes a feedback mechanism that keeps a constant tunneling current during the scan. The feedback can be either analog7–13 or digital.14–18 It has been long thought that the digital feedback needs a real-time acquisition and control system using a digital signal processor (DSP) or a field programmable gate array (FPGA) that operates separately from the user computer. This indeed allows using a clock and taking data at fixed time intervals. However, the STM often requires measurement of a parameter as a function of the position, not as a function of time. Here, we show that one can successfully operate a STM using a simple universal serial bus (USB) based data acquisition on a computer running the Windows® operating system. This eliminates the need for defining fixed time intervals for data acquisition. We also show that our setup allows for very low noise measurements. Addressing this problem usually requires careful filtering of the high frequency noise created at the control electronics.5,19–26

Let us remind that Tersoff and Hamann applied Bardeen’s tunneling theory to the STM27–29 and showed that under several simplifying assumptions and at zero temperature, the tunneling current vs bias voltage I(VBias) is given by

(1)

where z is the distance between the tip and the sample; zϕ8mϕ, where ϕ is the average work function, m is the electron effective mass and Ns(EeV) and Nt(E) are the energy dependent densities of states of the tip and the sample, respectively (E being the energy). The tip is atomically sharp, and its Nt(E) does not depend on the position of the tip over the sample. Thus, the 2D maps of the surface at a constant tunneling current provide maps of the local density of states Ns(x, y). The tunneling conductance dIdV as a function of the bias voltage is often proportional to the energy dependent density Ns(E = eV, x, y). dIdV can be recorded as a function of the position when scanning the tip on top of the surface. This leads to a set of 2D images, providing Ns(E = eV, x, y). Controlling the noise level is particularly important when there is a need to resolve sharp features in the tunneling conductance, dIdV, caused by variations of Ns(E) in a small energy range. This is often the case in the cryogenic STM setups used to study materials with large variations in the electronic band structure at low energy scales.

The simplified control electronics described here helps achieving the goal of low noise measurements more easily. We show results obtained in relevant materials and in several different setups that use the device described here. The full schematic of the electronics and the software are available in Refs. 30 and 31, respectively.

We describe the hardware of our system in Fig. 1. The system has two 16 bit Digital-to-Analog Converter (DAC) chips (DAC8734), each with four analog output lines, which provide voltage signals from −10 to +10 V. These signals can be amplified by a factor of 14 using the high precision and ultralow noise amplifiers described in Ref. 26. Acquisition of analog voltage signals is carried out using a 16 bit Analog-to-Digital Converter (ADC) chip (ADC7606) with eight separate input channels.

FIG. 1.

(a) Block diagram of the digital and analog part of the STM electronics. The components are marked by red boxes. Black dashed lines mark the different parts. Digital and analog parts are in a single enclosure. (b) Picture of the device. (c) Front end of the device. The USB digital interface (d) is screwed to the front panel. There are five ADC channels for reading signals (1), three auxiliary DAC signals (2), a digital auxiliary output (3), a USB port entry (4), an input for the amplifier used for coarse motion (5), and an auxiliary input for the Z control (6). All high voltage outputs (7) are on the back side (e).

FIG. 1.

(a) Block diagram of the digital and analog part of the STM electronics. The components are marked by red boxes. Black dashed lines mark the different parts. Digital and analog parts are in a single enclosure. (b) Picture of the device. (c) Front end of the device. The USB digital interface (d) is screwed to the front panel. There are five ADC channels for reading signals (1), three auxiliary DAC signals (2), a digital auxiliary output (3), a USB port entry (4), an input for the amplifier used for coarse motion (5), and an auxiliary input for the Z control (6). All high voltage outputs (7) are on the back side (e).

Close modal

To address the DAC and ADC systems, we use a USB 2.0 FIFO (first in first out) circuit (FTDI2232H), which addresses the serial peripheral interface (SPI). The USB is optocoupled to the board to avoid interference from the main computer. All SPI signals are filtered using one pole RC filters with a cutoff frequency slightly above the frequency of the clock required to address the DAC and ADC systems with the SPI (2 MHz). Outputs are carefully low pass filtered using analog one pole RC filters cutting at 10 kHz [Fig. 1(a)]. The digital part of the electronics [Fig. 1(d)] is mounted on a multilayer board, taking care of separating lines as far as possible, where the SPI bus runs from the signal input and output lines. Furthermore, it has a separate power supply and is mounted very close to the high voltage amplifiers to avoid picking up noise. Grounding is designed to avoid induction from digital signal lines by minimizing the impedance with large ground plates and a full connection to the metallic enclosure. In addition, microtransformers separate communication bus from signals to avoid interference. The noise of the output of the DAC is below 1 μV rms, and the ADC has an analog filter and a digital sampling averaging system, giving a high signal-to-noise ratio. The wiring between the electronics and the microscope is made by carefully mixing cables that hold the signals with ground cables, firmly attached to the electronics and the cryostat.

Five of the amplified output voltage signals are used to drive the scanner piezotube, four for the X and Y motion and one for the Z motion. The X and Y motion is performed using two DACs mounted in cascade, achieving thus μV resolution. This setup provides enough high voltage lines to operate a usual cryogenic microscope, as the ones discussed in Refs. 5 and 19–26. Another amplified voltage output line is used to drive the coarse approach motor. There are three further low voltage lines for the bias voltage and additional needs.

The electronics hardware is handy [Fig. 1(b)], and the whole system can be easily transported and connected to the USB port of a computer. Different systems can be used on the same computer by using separate USB ports.

We have written communication routines in Delphi.31 The feedback is made by reading the current and sending an output signal calculated using a PI (proportional integral) algorithm [Fig. 2(a)]. The latency of the USB port allows us to access the port with a bandwidth reaching a few kHz. This is below the typical resonance frequency of the STM, as required to avoid noise amplification and oscillations.4 We show an example of an obtained image on Co3Sn2S2 using a Au tip in Figs. 2(b)2(e). Further examples of similarly acquired images in other compounds and under different conditions are given below.

FIG. 2.

(a) Block diagram of the controlling system. (b)–(e) Example of a scanning experiment where the feedback loop is functioning to maintain a constant current, with the computer running as usual. In (b), we show the feedback signal as a function of the position converted to the z-position of the tip in Å. We can see terraces and atomic resolution obtained on Co3Sn2S2 scanned with a tip of Au (temperature of 4.2 K, magnetic field of 14 T, bias voltage of 100 mV, and current of 1.6 nA). The height color scale is given on the right. In (c), we show the current taken at each position of the image in (b). In (d), we show the height along one scan line, and in (e), the current read along the same line [yellow line in (b) and (c)]. Black scale bars in (b) and (c) are 5 nm long. (f) Interval between two points when measuring a signal as a function of time during a 30 s period approximately each 3.2 ms.

FIG. 2.

(a) Block diagram of the controlling system. (b)–(e) Example of a scanning experiment where the feedback loop is functioning to maintain a constant current, with the computer running as usual. In (b), we show the feedback signal as a function of the position converted to the z-position of the tip in Å. We can see terraces and atomic resolution obtained on Co3Sn2S2 scanned with a tip of Au (temperature of 4.2 K, magnetic field of 14 T, bias voltage of 100 mV, and current of 1.6 nA). The height color scale is given on the right. In (c), we show the current taken at each position of the image in (b). In (d), we show the height along one scan line, and in (e), the current read along the same line [yellow line in (b) and (c)]. Black scale bars in (b) and (c) are 5 nm long. (f) Interval between two points when measuring a signal as a function of time during a 30 s period approximately each 3.2 ms.

Close modal

To have an optimized time base for execution, we use threaded timer and program its access in such a way that it is regularly called during time consuming operations such as loops for data saving.32 This provides smooth operation without interruptions.33,34 We have followed data acquisition and PI operation for many hours, acquiring billions of points and observing random interruptions of the latency, which were significantly below 1 ms and just a few interruptions requiring up to 10 ms. None of these interruptions lead to a tip instability, as shown, for example, in Figs. 2(b)2(e), which provide raw data (no corrected points). To better quantify the jitter in the latency, we have measured a signal as a function of time with a wanted time interval of 3.2 ms. In Fig. 2(f), we show the interval between two measurements. We see that we have a deviation of 0.2 ms from the wanted interval in less than 1% of the acquired points and deviations larger than 0.4 ms in far less than 0.1% of the acquired points. Thus, although the system is not a priori thought for real-time measurements, it can be used in measurements as a function of time where a constant time interval is not critical in the acquisition.

We show a block diagram of the components of the software (blue background) and how these are linked to the hardware (red and white backgrounds) in Fig. 3. The software interacts with the hardware through a main window (“WinSPM”). The main window includes two additional windows. One is used to make spectroscopy (“spectroscopy IV/IZ”) or any other operation that requires cutting the feedback loop, and the other (“STM topography”) is used to show and control the scanning process. The feedback (“digital feedback”) operates continuously in the background, reading the current value and acting on the Z electrode of the piezotube through the corresponding amplifier, whose value is also sent to the “STM topography” to monitor the Z as a function of the position of the tip during scanning. The main window handles the USB communication through the interface containing the DAC and ADC (“LHA-DigiMod”). One converter is used for setting the bias voltage, and others are amplified to set the x-, y-, and z-positions. The ADC is used, for instance, to read the current from the current to voltage transimpedance amplifier (“I/V converter”). In addition, another high voltage amplifier is used to drive the coarse approach motor. During imaging, we can stop at each point and perform IV or bias voltage vs conductance curves to obtain images of the current or the conductance as a function of the bias voltage (“scanning tunneling spectroscopy”).

FIG. 3.

Block diagram of the different operations of the data acquisition software and its communication with the digital electronics. There is a coarse positioning module (“Z′ positioning”), a spectroscopy module to take current vs voltage (IV) or current vs distance (IZ) curves (“Spectroscopy IV/IZ”) or other curves (such as the conductance vs bias voltage), and an image module (“STM topography”). The program controls the electronics (red box) through a USB controlled interface called LHA-DigiMod, which is located on the front of the electronics [Fig. 1(c)]. The STM setup, with the piezotube, approach system, tip, and sample, is indicated by the dark green box marked “STM.” Some additional equipment such as the voltage divider and the current to voltage converter are indicated in white boxes. The digital feedback module, schematically shown in light blue, runs to maintain a constant current during operation, fed by the current measurement and acting on the Z-piezo. The signal sent to the Z-piezo is used to trace the STM topography.

FIG. 3.

Block diagram of the different operations of the data acquisition software and its communication with the digital electronics. There is a coarse positioning module (“Z′ positioning”), a spectroscopy module to take current vs voltage (IV) or current vs distance (IZ) curves (“Spectroscopy IV/IZ”) or other curves (such as the conductance vs bias voltage), and an image module (“STM topography”). The program controls the electronics (red box) through a USB controlled interface called LHA-DigiMod, which is located on the front of the electronics [Fig. 1(c)]. The STM setup, with the piezotube, approach system, tip, and sample, is indicated by the dark green box marked “STM.” Some additional equipment such as the voltage divider and the current to voltage converter are indicated in white boxes. The digital feedback module, schematically shown in light blue, runs to maintain a constant current during operation, fed by the current measurement and acting on the Z-piezo. The signal sent to the Z-piezo is used to trace the STM topography.

Close modal

We provide an example of screenshots during a typical experiment in Fig. 4. The experiment is made at 100 mK with a tip of Au and a sample of FeSe. In Fig. 4(a), we show the main working window, where we have pasted an image showing steps in FeSe. The large square shows an image of an area of several nm square, showing steps in FeSe, and another image taken around the central area, showing atomic resolution. The blue cross in the center of the square is the position of the STM tip. In Fig. 4(d), we show a typical scan (red line in the left panel) and the image being built up on the right panels. We can clearly see the profiles of Se atoms in the line scan and the square atomic Se lattice in the images. In Fig. 4(e), we show a current vs voltage curve (blue and red overlapped lines) obtained in FeSe under a magnetic field of 14.5 T. We see that the curve is slightly nonlinear due to the superconducting gap opening.

FIG. 4.

Screenshots from the data acquisition software. (a) Scanner interface. The blue cross and blue square mark the position of the tip and the scanning area, respectively. We depict terraces and an atomically resolved image on one terrace as colored images (from black to yellow to white). Data are obtained on FeSe at 100 mK and at a magnetic field of 14.5 T with a bias voltage of 5 mV and tunneling current of 4 nA. (b) Digital PI (proportional integral) feedback interface window. (c) Coarse vertical motion (Z′) control window. (d) Interface for the live preview of the scan. The central panel shows the live tip height profile, and the progress of the topographic images made by scanning on the ±Y direction is shown in the right panels (top panel for scanning along one direction in Y and bottom panel for scanning back). (e) Curve acquisition interface, showing the current vs voltage as blue and red lines, respectively.

FIG. 4.

Screenshots from the data acquisition software. (a) Scanner interface. The blue cross and blue square mark the position of the tip and the scanning area, respectively. We depict terraces and an atomically resolved image on one terrace as colored images (from black to yellow to white). Data are obtained on FeSe at 100 mK and at a magnetic field of 14.5 T with a bias voltage of 5 mV and tunneling current of 4 nA. (b) Digital PI (proportional integral) feedback interface window. (c) Coarse vertical motion (Z′) control window. (d) Interface for the live preview of the scan. The central panel shows the live tip height profile, and the progress of the topographic images made by scanning on the ±Y direction is shown in the right panels (top panel for scanning along one direction in Y and bottom panel for scanning back). (e) Curve acquisition interface, showing the current vs voltage as blue and red lines, respectively.

Close modal

To obtain maps of IV or tunneling conductance vs bias voltage curves at each point of an image, we disconnect the feedback loop and perform the measurement at each point of the scan. We store separately the topography file with usual scanning information and the IV curves.

The file containing all IV (or tunneling conductance vs bias voltage) curves is actually a 3D matrix with an IV (or tunneling conductance) curve at each pixel of a 2D space that spans the real space image. We analyze these files using software based on the Matlab® environment.35 Related software based on Matlab can be found in Ref. 36. Our software includes numerous features required to analyze the STM images, such as Fourier transform spectroscopy, rotation and manipulation of images, plot of conductance maps at any bias voltage, and obtaining and analyzing tunneling conductance curves from anywhere in the images.

We show a block diagram of the data treatment software in Fig. 5. This includes a module for single curves, a reading module, and an image processing module.

FIG. 5.

Block diagram of the data analysis software. The data are stored in a binary file with the required headers to return image size, position, etc. There is a curve reading module (“single curves”), which allows rendering the data and performing basic operations to understand tunneling current curves, such as derivation, temperature deconvolution, and plots including many curves. The “scanning tunneling spectroscopy” module is separated into a “reading” module, which creates current and conductance matrices and, from these, conductance and Fourier transform maps at given bias voltage values. The “image processing” module allows rendering the images in a user-friendly manner and making different operations required to understand the real and Fourier space images. The latter is easily extended using Matlab functions.

FIG. 5.

Block diagram of the data analysis software. The data are stored in a binary file with the required headers to return image size, position, etc. There is a curve reading module (“single curves”), which allows rendering the data and performing basic operations to understand tunneling current curves, such as derivation, temperature deconvolution, and plots including many curves. The “scanning tunneling spectroscopy” module is separated into a “reading” module, which creates current and conductance matrices and, from these, conductance and Fourier transform maps at given bias voltage values. The “image processing” module allows rendering the images in a user-friendly manner and making different operations required to understand the real and Fourier space images. The latter is easily extended using Matlab functions.

Close modal

The reading module needs two main inputs: the n by m dimensions of the topographic image and the corresponding binary file containing the associated spectroscopic curves and the starting point in the binary file. We then average data using adjacent points around an interval δ(V). Normalization with respect to the conductance value at a certain bias voltage range can be used to eliminate setpoint effects, as discussed in Ref. 37. Data can be saved into a structure object that stores all the variables needed for a further image processing analysis.

We show an example and screenshots in Fig. 6. In Fig. 6(a), we show the acquisition module, which we use to open files and set a few relevant parameters. In Fig. 6(b), we show the IV curve reader with an IV curve taken in β–Bi2Pd at 100 mK and 0.1 T, shown as a blue line. We show the IV curve (left panel) and the tunneling conductance vs bias voltage (right panel). The red and blue vertical lines in the right panel show the voltage range that we use to normalize curves in this particular case. The resulting maps of the tunneling conductance are shown in Figs. 6(c) and 6(d), where we plot the zero bias conductance as a function of the position [Fig. 6(c)] and its Fourier transform [Fig. 6(d)]. The images show the superconducting vortex lattice of β–Bi2Pd at 0.1 T and 100 mK.38–40 We can visualize results at different bias voltages, identify vortices, or triangulate vortex lattices.41–43 We can also use the same modules to visualize oscillations in the tunneling conductance due to impurity scattering, called quasiparticle interference. Quasiparticle interference measures the electronic band structure.44–53 We can perform different operations of quasiparticle interference, such as rotations, symmetrization, filtering, or extraction of profiles. The module for single curves is shown in Fig. 6(d), with an example of an IV curve and the corresponding conductance obtained by using Al as the tip and the sample.

FIG. 6.

Screenshots from the data analysis software. (a) Main window of the program with the different modules separated in menu items. (b) Conductance curve calculation window. We can see the superconducting gap at a magnetic field of 0.1 T obtained in β–Bi2Pd as blue lines. (c) Window for viewing conductance maps. We show a picture of the vortex lattice in β–Bi2Pd obtained from the zero bias conductance map. (d) Window showing the Fourier transform of the image in (c). A number of operations can be performed on the Fourier transform, such as rotating, filtering, and symmetrizing. (e) Main window to plot and make calculations on tunneling current and tunneling conductance vs bias voltage curves. (f) Tunneling current (left panel) and tunneling conductance (right panel) obtained using the Al tip and the sample below 100 mK.

FIG. 6.

Screenshots from the data analysis software. (a) Main window of the program with the different modules separated in menu items. (b) Conductance curve calculation window. We can see the superconducting gap at a magnetic field of 0.1 T obtained in β–Bi2Pd as blue lines. (c) Window for viewing conductance maps. We show a picture of the vortex lattice in β–Bi2Pd obtained from the zero bias conductance map. (d) Window showing the Fourier transform of the image in (c). A number of operations can be performed on the Fourier transform, such as rotating, filtering, and symmetrizing. (e) Main window to plot and make calculations on tunneling current and tunneling conductance vs bias voltage curves. (f) Tunneling current (left panel) and tunneling conductance (right panel) obtained using the Al tip and the sample below 100 mK.

Close modal

As an example of data treatment, we can discuss tunneling spectroscopy vs temperature measurements in Pb using a tip of Au, as shown in Fig. 7. The experimental curves obtained are shown as colored lines in Fig. 7(a). To understand these curves, we construct the density of states as a function of the energy N(E) [shown in Fig. 7(b), bottom light blue line]. N(E) closely follows the BCS expression [N(E)=EE2Δ2, where Δ is the superconducting gap], with the value of the superconducting gap expected for Pb (Δ = 1.37 meV), except that the divergence at eV = Δ is substituted by non-infinite quasiparticle peaks. This can be due to a combination of inelastic nonequilibrium quasiparticles, the anisotropy of the superconducting gap of Pb, and the finite energy resolution of our experiment.5,19–26,54,55 Using this N(E), we can then calculate the tunneling conductance for all temperatures by convoluting with the derivative of the Fermi function. The results are shown as black lines in Fig. 7(a). To obtain the tunneling conductance, we modify N(E) at each temperature until we reproduce the experiment. From the position of the quasiparticle peaks in N(E), we obtain Δ(T). We see that it follows the BCS superconducting gap dependence, Δ(T), as shown in Fig. 7(c).

FIG. 7.

(a) Tunneling conductance vs bias voltage taken at 0.15, 0.25, 1, 3, 5, 6, 7, and 7.5 K in Pb using a tip of Au. The black lines are the convolution of the density of states [shown in (b)] with the Fermi function. The position of the quasiparticle peaks provides the value of the superconducting gap, shown in the main panel of (c). The black line is a fit to the BCS temperature dependence of the superconducting gap Δ(T). In the inset, we show the density of states at zero energy as a function of temperature. We see that it remains zero up to the point when superconductivity vanishes.

FIG. 7.

(a) Tunneling conductance vs bias voltage taken at 0.15, 0.25, 1, 3, 5, 6, 7, and 7.5 K in Pb using a tip of Au. The black lines are the convolution of the density of states [shown in (b)] with the Fermi function. The position of the quasiparticle peaks provides the value of the superconducting gap, shown in the main panel of (c). The black line is a fit to the BCS temperature dependence of the superconducting gap Δ(T). In the inset, we show the density of states at zero energy as a function of temperature. We see that it remains zero up to the point when superconductivity vanishes.

Close modal

We now operate in our laboratory several units of the system described here. In Fig. 8, we show results of topographic images at constant current in WTe2, 2H–NbSe2, Co3Sn2S2, and FeSe. We see the one-dimensional Te chains that usually appear on the surface of WTe2 [Fig. 8(a)]. We also see the typical atomic lattice and charge density wave that are characteristic of 2H–NbSe2 [Fig. 8(b)]. The charge density wave is hexagonal and produces a charge modulation, which is nearly commensurate with a period three times the lattice constant. There are two vacancies visible in the surface lattice, which consists of a Se hexagonal lattice. In Co3Sn2S2 [Fig. 8(c)], we observe an atomically resolved surface, consisting of two different atomic species. The continuous hexagonal layer is composed of Sn atoms. Notice their hexagonal shape. The upper incomplete layer is made of S atoms. There is a tendency to form rows, and the S atoms induce shadow-like features in the neighboring Sn atoms. In FeSe [Fig. 8(d)], we observe the square Se lattice. There are defects with an elongated in-plane shape characteristic of the nematic features of the electronic properties of this compound.

FIG. 8.

(a) Atomic resolution image obtained in WTe2 with a bias voltage Vbias= 100 mV and tunneling current Itun = 4 nA. (b) STM image obtained in 2H–NbSe2 with Vbias = 100 mV and Itun = 3.2 nA. (c) STM image obtained in Co3Sn2S2 with Vbias = 100 mV and Itun = 1.6 nA. (d) STM image obtained in FeSe, with Vbias= 10 mV and Itun = 4 nA. All images have been taken at a temperature of 4.2 K (a)–(c) and 100 mK (d). Scale bars in black are 2 nm long. Vertical bars provide the z-changes as a color scale. Images are raw, except for the removal of a plane.

FIG. 8.

(a) Atomic resolution image obtained in WTe2 with a bias voltage Vbias= 100 mV and tunneling current Itun = 4 nA. (b) STM image obtained in 2H–NbSe2 with Vbias = 100 mV and Itun = 3.2 nA. (c) STM image obtained in Co3Sn2S2 with Vbias = 100 mV and Itun = 1.6 nA. (d) STM image obtained in FeSe, with Vbias= 10 mV and Itun = 4 nA. All images have been taken at a temperature of 4.2 K (a)–(c) and 100 mK (d). Scale bars in black are 2 nm long. Vertical bars provide the z-changes as a color scale. Images are raw, except for the removal of a plane.

Close modal

The images are provided raw, implying that, during their acquisition, which needed between half an hour and five hours, the data acquisition was fast enough so that there were no significant interruptions that influenced the stability of the PI algorithm controlling the tip z-position.

In summary, we have described a system where we successfully control a STM microscope using a simple USB based data acquisition. Details of hardware and software are available in Refs. 30 and 31, respectively.

Our approach builds on developments of the so-called soft real-time data acquisition, as opposed to data acquisition on an exactly periodic time basis.33,34 Our results show that for the STM and probably for many other applications too, computers allow for an access latency for data acquisition in usual input and output ports that is enough for successful operation.

Care should be taken to use the method proposed here in combination with a PI algorithm. The system stability should be tested so as to make sure that random interruptions do not influence the PI stability. At least for the systems imaged here and probably for many other systems that require high resolution and thus slow data acquisition, the latency for accessing the port of a usual desktop computer is largely sufficient to maintain the stability of the PI. Furthermore, measurements as a function of time require obviously a stable time basis. This can be circumvented to some extent by measuring the time difference between acquisitions but leads to curves with random changes in the time base. Our results show that these changes are mostly below a ms and only rarely extend over more than a few ms.

There are commercially available choices for data acquisition using the USB port that include data processors.56,57 Using appropriate filtering, a similar noise level as we obtain here could be achieved.15 However, we believe that its implementation is more difficult than using the simple solution proposed here. Furthermore, calculations required to acquire the actual data are made within the data acquisition systems. We believe instead that allowing the user to perform operations on the main controlling computer provides a significant improvement. Multiple core computers with different components running in parallel are becoming largely available. This will produce a shift toward integrating data acquisition within the running operational system environments. The results obtained here could then be applied to all kinds of scanning probe microscopes. The user will be able to considerably increase the possibilities for making complex operations during scanning.

We acknowledge discussions with Nicolás Agraït and Sebastián Vieira. This work was supported by the Spanish Research State Agency (Grant Nos. FIS2017-84 330-R, pPID2020-114071RB-100, CEX2018-000805-M, RYC-2014-15093, MAT2017-87134-C2-2-R, and MAT2017-88693-R), the Comunidad de Madrid through program NANOFRONTMAG-CM (Grant No. S2018/NMT-4321), and the European Research Council PNICTEYES through Grant Agreement No. 679080. We acknowledge collaborations through EU program Cost CA16218 (Nanocohybri). Work at Ames Laboratory was supported by the U.S. Department of Energy, Office of Basic Energy Science, Division of Materials Sciences and Engineering. Ames Laboratory is operated for the U.S. Department of Energy by Iowa State University under Contract No. DE-AC02-07CH11358. N.H.J. and J.S. were also supported by the Gordon and Betty Moore Foundation’s EPiQS Initiative through Grant No. GBMF4411.

The authors have no conflicts to disclose.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
G.
Binnig
,
H.
Rohrer
,
C.
Gerber
, and
E.
Weibel
, “
Tunneling through a controllable vacuum gap
,”
Appl. Phys. Lett.
40
,
178
180
(
1982
).
2.
Y.
Kuk
and
P. J.
Silverman
, “
Scanning tunneling microscope instrumentation
,”
Rev. Sci. Instrum.
60
,
165
180
(
1989
).
3.
R.
Wiesendanger
,
Scanning Probe Microscopy and Spectroscopy: Methods and Applications
(
Cambridge University Press
,
1994
).
4.
B.
Voigtlaender
,
Scanning Probe Microscopy
(
Springer-Verlag Berlin Heidelberg
,
2015
).
5.
Y. J.
Song
,
A. F.
Otte
,
V.
Shvarts
,
Z.
Zhao
,
Y.
Kuk
,
S. R.
Blankenship
,
A.
Band
,
F. M.
Hess
, and
J. A.
Stroscio
, “
Invited review article: A 10 mK scanning probe microscopy facility
,”
Rev. Sci. Instrum.
81
,
121101
(
2010
).
6.
F. J.
Giessibl
, “
Advances in atomic force microscopy
,”
Rev. Mod. Phys.
75
,
949
983
(
2003
).
7.
S.
Grafstrom
,
J.
Kowalski
, and
R.
Neumann
, “
Design and detailed analysis of a scanning tunnelling microscope
,”
Meas. Sci. Technol.
1
,
139
146
(
1990
).
8.
H.
Olin
, “
Design of a scanning probe microscope
,”
Meas. Sci. Technol.
5
,
976
984
(
1994
).
9.
E.
Flaxer
, “
Compact programmable controller for a linear piezo-stepper motor
,”
Mechatronics
16
,
303
308
(
2006
).
10.
E.
Anguiano
,
A. I.
Oliva
,
M.
Aguilar
, and
J. L.
Peña
, “
Analysis of scanning tunneling microscopy feedback system: Experimental determination of parameters
,”
Rev. Sci. Instrum.
67
,
2947
2952
(
1996
).
11.
C. Y.
Nakakura
,
V. M.
Phanse
,
G.
Zheng
,
G.
Bannon
,
E. I.
Altman
, and
K. P.
Lee
, “
A high-speed variable-temperature ultrahigh vacuum scanning tunneling microscope
,”
Rev. Sci. Instrum.
69
,
3251
3258
(
1998
).
12.
L.
Kuipers
,
R. W. M.
Loos
,
H.
Neerings
,
J.
ter Horst
,
G. J.
Ruwiel
,
A. P.
de Jongh
, and
J. W. M.
Frenken
, “
Design and performance of a high temperature, high speed scanning tunneling microscope
,”
Rev. Sci. Instrum.
66
,
4557
4565
(
1995
).
13.
O.
Haase
,
M.
Borbonus
,
P.
Muralt
,
R.
Koch
, and
K. H.
Rieder
, “
A novel ultrahigh vacuum scanning tunneling microscope for surface science studies
,”
Rev. Sci. Instrum.
61
,
1480
1483
(
1990
).
14.
R.
Piner
and
R.
Reifenberger
, “
Computer control of the tunnel barrier width for the scanning tunneling microscope
,”
Rev. Sci. Instrum.
60
,
3123
3127
(
1989
).
15.
S. J.
Yu
,
E.
Fajeau
,
L. Q.
Liu
,
D. J.
Jones
, and
K. W.
Madison
, “
The performance and limitations of FPGA based digital servos for atomic, molecular, and optical physics experiments
,”
Rev. Sci. Instrum.
89
,
025107
(
2018
).
16.
P.
Zahl
,
T.
Wagner
,
R.
Möller
, and
A.
Klust
, “
Open source scanning probe microscopy control software package GXSM
,”
J. Vac. Sci. Technol. B
28
,
C4E39
C4E47
(
2010
).
17.
D. R.
Baselt
,
S. M.
Clark
,
M. G.
Youngquist
,
C. F.
Spence
, and
J. D.
Baldeschwieler
, “
Digital signal processor control of scanned probe microscopes
,”
Rev. Sci. Instrum.
64
,
1874
1882
(
1993
).
18.
I.
Horcas
,
R.
Fernández
,
J. M.
Gómez-Rodríguez
,
J.
Colchero
,
J.
Gómez-Herrero
, and
A. M.
Baro
, “
WSXM: A software for scanning probe microscopy and a tool for nanotechnology
,”
Rev. Sci. Instrum.
78
,
013705
(
2007
).
19.
H.
Suderow
,
I.
Guillamon
, and
S.
Vieira
, “
Compact very low temperature scanning tunneling microscope with mechanically driven horizontal linear positioning stage
,”
Rev. Sci. Instrum.
82
,
033711
(
2011
).
20.
M.
Assig
,
M.
Etzkorn
,
A.
Enders
,
W.
Stiepany
,
C. R.
Ast
, and
K.
Kern
, “
A 10 mK scanning tunneling microscope operating in ultra high vacuum and high magnetic fields
,”
Rev. Sci. Instrum.
84
,
033903
(
2013
).
21.
I.
Battisti
,
G.
Verdoes
,
K.
van Oosten
,
K. M.
Bastiaans
, and
M. P.
Allan
, “
Definition of design guidelines, construction, and performance of an ultra-stable scanning tunneling microscope for spectroscopic imaging
,”
Rev. Sci. Instrum.
89
,
123705
(
2018
).
22.
S.
Misra
,
B. B.
Zhou
,
I. K.
Drozdov
,
J.
Seo
,
L.
Urban
,
A.
Gyenis
,
S. C. J.
Kingsley
,
H.
Jones
, and
A.
Yazdani
, “
Design and performance of an ultra-high vacuum scanning tunneling microscope operating at dilution refrigerator temperatures and high magnetic fields
,”
Rev. Sci. Instrum.
84
,
103903
(
2013
).
23.
T.
Machida
,
Y.
Kohsaka
, and
T.
Hanaguri
, “
A scanning tunneling microscope for spectroscopic imaging below 90 mK in magnetic fields up to 17.5 T
,”
Rev. Sci. Instrum.
89
,
093707
(
2018
).
24.
M.
Marz
,
G.
Goll
, and
H. V.
Löhneysen
, “
A scanning tunneling microscope for a dilution refrigerator
,”
Rev. Sci. Instrum.
81
,
045102
(
2010
).
25.
Q.
Li
,
Q.
Wang
,
Y.
Hou
, and
Q.
Lu
, “
18/20 T high magnetic field scanning tunneling microscope with fully low voltage operability, high current resolution, and large scale searching ability
,”
Rev. Sci. Instrum.
83
,
043706
(
2012
).
26.
J. A.
Galvis
,
E.
Herrera
,
I.
Guillamón
,
J.
Azpeitia
,
R. F.
Luccas
,
C.
Munuera
,
M.
Cuenca
,
J. A.
Higuera
,
N.
Díaz
,
M.
Pazos
,
M.
García-Hernandez
,
A.
Buendía
,
S.
Vieira
, and
H.
Suderow
, “
Three axis vector magnet set-up for cryogenic scanning probe microscopy
,”
Rev. Sci. Instrum.
86
,
013706
(
2015
).
27.
J.
Bardeen
, “
Tunnelling from a many-particle point of view
,”
Phys. Rev. Lett.
6
,
57
59
(
1961
).
28.
J.
Tersoff
and
D. R.
Hamann
, “
Theory and application for the scanning tunneling microscope
,”
Phys. Rev. Lett.
50
,
1998
2001
(
1983
).
29.
J.
Tersoff
and
D. R.
Hamann
, “
Theory of the scanning tunneling microscope
,”
Phys. Rev. B
31
,
805
813
(
1985
).
30.
See https://osf.io/by4aj/ for schematics of the electronics.
31.
See https://github.com/LowTemperaturesUAM/MyScanner for the code of the software.
32.
We use Threaded Timer v.1.24, https://torry.net/pages.php?id=294.
33.
N.
Korver
, “
Adequacy of the universal serial bus for real-time systems
,”
2003
, https://core.ac.uk/download/pdf/11461022.pdf.
34.
G.
Buttazzo
,
L.
Lipari
, and
M.
Abeni
,
Soft Real-Time Systems
, Series in Computer Science (
Springer
,
2005
).
35.
MATLAB R2019b,
The Mathworks, Inc.
,
Natick, MA
,
2019
.
37.
M. J.
Lawler
,
K.
Fujita
,
J.
Lee
,
A. R.
Schmidt
,
Y.
Kohsaka
,
C. K.
Kim
,
H.
Eisaki
,
S.
Uchida
,
J. C.
Davis
,
J. P.
Sethna
, and
E.-A.
Kim
, “
Intra-unit-cell electronic nematicity of the high-Tc copper-oxide pseudogap states
,”
Nature
466
,
347
351
(
2010
).
38.
Ø.
Fischer
,
M.
Kugler
,
I.
Maggio-Aprile
,
C.
Berthod
, and
C.
Renner
, “
Scanning tunneling spectroscopy of high-temperature superconductors
,”
Rev. Mod. Phys.
79
,
353
419
(
2007
).
39.
H.
Suderow
,
I.
Guillamón
,
J. G.
Rodrigo
, and
S.
Vieira
, “
Imaging superconducting vortex cores and lattices with a scanning tunneling microscope
,”
Supercond. Sci. Technol.
27
,
063001
(
2014
).
40.
E.
Herrera
,
I.
Guillamón
,
J. A.
Galvis
,
A.
Correa
,
A.
Fente
,
R. F.
Luccas
,
F. J.
Mompean
,
M.
García-Hernández
,
S.
Vieira
,
J. P.
Brison
, and
H.
Suderow
, “
Magnetic field dependence of the density of states in the multiband superconductor β–Bi2Pd
,”
Phys. Rev. B
92
,
054507
(
2015
).
41.
I.
Guillamón
,
H.
Suderow
,
A.
Fernández-Pacheco
,
J.
Sesé
,
R.
Córdoba
,
J. M.
De Teresa
,
M. R.
Ibarra
, and
S.
Vieira
, “
Direct observation of melting in a two-dimensional superconducting vortex lattice
,”
Nat. Phys.
5
,
651
655
(
2009
).
42.
I.
Guillamón
,
R.
Córdoba
,
J.
Sesé
,
J. M.
De Teresa
,
M. R.
Ibarra
,
S.
Vieira
, and
H.
Suderow
, “
Enhancement of long-range correlations in a 2D vortex lattice by an incommensurate 1D disorder potential
,”
Nat. Phys.
10
,
851
856
(
2014
).
43.
J. B.
Llorens
,
I.
Guillamón
,
I. G.
Serrano
,
R.
Córdoba
,
J.
Sesé
,
J. M.
De Teresa
,
M. R.
Ibarra
,
S.
Vieira
,
M.
Ortuño
, and
H.
Suderow
, “
Disordered hyperuniformity in superconducting vortex lattices
,”
Phys. Rev. Res.
2
,
033133
(
2020
).
44.
J. E.
Hoffman
,
K.
McElroy
,
D.-H.
Lee
,
K. M.
Lang
,
H.
Eisaki
,
S.
Uchida
, and
J. C.
Davis
, “
Imaging quasiparticle interference in Bi2Sr2CaCu2O8+δ
,”
Science
297
,
1148
1151
(
2002
).
45.
P. T.
Sprunger
,
L.
Petersen
,
E. W.
Plummer
,
E.
Lægsgaard
, and
F.
Besenbacher
, “
Giant Friedel oscillations on the beryllium(0001) surface
,”
Science
275
,
1764
1767
(
1997
).
46.
K.
McElroy
,
R. W.
Simmonds
,
J. E.
Hoffman
,
D.-H.
Lee
,
J.
Orenstein
,
H.
Eisaki
,
S.
Uchida
, and
J. C.
Davis
, “
Relating atomic-scale electronic phenomena to wave-like quasiparticle states in superconducting Bi2Sr2CaCu2O8+δ
,”
Nature
422
,
592
596
(
2003
).
47.
A.
Stróżecka
,
A.
Eiguren
, and
J. I.
Pascual
, “
Quasiparticle interference around a magnetic impurity on a surface with strong spin-orbit coupling
,”
Phys. Rev. Lett.
107
,
186805
(
2011
).
48.
J. E.
Hoffman
, “
Spectroscopic scanning tunneling microscopy insights into Fe-based superconductors
,”
Rep. Prog. Phys.
74
,
124513
(
2011
).
49.
J.
Lee
,
K.
Fujita
,
A. R.
Schmidt
,
C. K.
Kim
,
H.
Eisaki
,
S.
Uchida
, and
J. C.
Davis
, “
Spectroscopic fingerprint of phase-incoherent superconductivity in the underdoped Bi2Sr2CaCu2O8+Δ
,”
Science
325
,
1099
1103
(
2009
).
50.
K.
Iwaya
,
Y.
Kohsaka
,
K.
Okawa
,
T.
Machida
,
M. S.
Bahramy
,
T.
Hanaguri
, and
T.
Sasagawa
, “
Full-gap superconductivity in spin-polarised surface states of topological semimetal β-PdBi2
,”
Nat. Commun.
8
,
976
(
2017
).
51.
S.-Y.
Xu
,
C.
Liu
,
S. K.
Kushwaha
,
R.
Sankar
,
J. W.
Krizan
,
I.
Belopolski
,
M.
Neupane
,
G.
Bian
,
N.
Alidoust
,
T.-R.
Chang
,
H.-T.
Jeng
,
C.-Y.
Huang
,
W.-F.
Tsai
,
H.
Lin
,
P. P.
Shibayev
,
F.-C.
Chou
,
R. J.
Cava
, and
M. Z.
Hasan
, “
Observation of Fermi arc surface states in a topological metal
,”
Science
347
,
294
298
(
2015
).
52.
H.
Inoue
,
A.
Gyenis
,
Z.
Wang
,
J.
Li
,
S. W.
Oh
,
S.
Jiang
,
N.
Ni
,
B. A.
Bernevig
, and
A.
Yazdani
, “
Quasiparticle interference of the Fermi arcs and surface-bulk connectivity of a Weyl semimetal
,”
Science
351
,
1184
1187
(
2016
).
53.
J. B.
Llorens
,
E.
Herrera
,
V.
Barrena
,
B.
Wu
,
N.
Heinsdorf
,
V.
Borisov
,
R.
Valentí
,
W. R.
Meier
,
S.
Bud’ko
,
P. C.
Canfield
,
I.
Guillamón
, and
H.
Suderow
, “
Anisotropic superconductivity in the spin-vortex antiferromagnetic superconductor CaK(Fe0.95Ni0.05)4As4
,”
Phys. Rev. B
103
,
L060506
(
2021
).
54.
M.
Ruby
,
B. W.
Heinrich
,
J. I.
Pascual
, and
K. J.
Franke
, “
Experimental demonstration of a two-band superconducting state for lead using scanning tunneling spectroscopy
,”
Phys. Rev. Lett.
114
,
157001
(
2015
).
55.
J. G.
Rodrigo
,
H.
Suderow
, and
S.
Vieira
, “
On the use of STM superconducting tips at very low temperatures
,”
Eur. Phys. J. B
40
,
483
488
(
2004
).
57.
See https://www.redpitaya.com/ for another example of USB data acquisition software.