2.1 Realization of the chip
The electromigration is mainly dominated by the current density [
9] and by the temperature of the wire [
10]. Both these quantities can be controlled by a proper geometry of the probe [
11] and by the applied voltage waveform [
12]. Using electromigration as technique for creating nanogaps, the wire has to have optimized geometries to facilitate the phenomenon and make it more controllable. From this point, to avoid a too high input current, it is necessary to have a small section of the wire [
11]; moreover, if the section of the wire is too small, the thermal conductance decreases and the temperature of the wire tends to become excessive, leading to the melting of the wire.
For these reasons, it is fundamental to have a model that can simulate the behavior of the physical phenomena during the electromigration. In particular, it is interesting to know and to anticipate the temperature inside the gold wire and for this it is necessary to model the geometry of the probe with a software simulation tool. The objective is therefore to achieve a wire section as small as possible, to trigger the electromigration also at lower current, obtaining a better control and a more regular shape. Figure shows the geometry of our probe described in Comsol Multi-physics, while Figure shows the single wire. To evaluate the temperature behavior during the voltage application, we performed a large number of simulations using different values of wire length, keeping in mind that there is a lower limit for the length of the wire that will allow electromigration to occur (the Blech length [
13-
15]). Figure shows the best parameters for generating a temperature profile quite sharp in a way that should be possible to focus the electromigration phenomenon in the center of the wire. The temperature profile along wire length was modeled with the equations:
Equation (1) is the law of conservation of energy: Q is the power transferred to a point, ρ is the mass volume density, Cp is the heat capacity, T is the temperature, t is the time and K is the thermal conductivity. Equation (2) gives the power dissipation (Joule effect): σ is the electrical conductivity and V is the electrical potential. Finally, Equation (3) is the Ohm's law in local form: J is the current density. All the equations are related to an infinitesimal point of space. Comsol uses the nodes of the mesh geometry to make a spatial sampling and integrates the equations in the volume using the nodes as points of integration: if the sampling is quite dense, the error is negligible. As result we have obtained the plot shown in Figure that represents the temperature variation in the wire length. The geometry created also allows a uniform distribution of current density (Figure ).
For having a useful platform where to produce the nanogaps, a silicon chip was realized, containing eight gold probes as shown in Figure , each of these connectable by bonding. In this way, it is possible to realize on the same chip eight nanogap structures, and each one is independent, so an high number of measurements is individually achievable. The final dimension of the chip is 2.4 × 4.1 mm, giving the possibility to insert it in heads of instruments as a cryostat or FESEM/AFM/STM microscopes, for doing for example measurement in vacuum and at very low temperatures. The chip is also ready to be wire bonded to a PCB (Figure ). It is possible to perform wet analysis too, for molecule characterizations, just spinning on the chip the solution that has to be measured. Obviously molecules in solutions must have some suitable sites for bonding with gold, such as thiol groups, in this way it is possible to obtain the desired M-M-M structure.
The realization of the chip starts from a silicon wafer capped by 200 nm of SiO2; the wafer is, then, inserted in a plasma oxygen machine for increasing the oxygen atoms concentration on the surface. Hydroxylation process is developed with a piranha solution, and so, after rinsing and drying, the surface of the wafer exposes -OH groups, fundamental for the anchoring of the organic compound that is evaporated on. In fact, to promote the adhesion of gold on the SiO2 surface, an insulating layer of MPTMS
(3-Mercaptopropyl)trimethoxysilane is then deposited [
16,
17]. The insulation property is important because has to be avoided any alternative path of current, that has to flow through the gold layer only. After, a gold film of 25 nm is grown on the MPTMS layer using an EVA600 evaporator. The wafer is inserted in a EVG-SUSS MA6 where a photolithographic mask process is performed employing a positive photoresist, afterwards the gold etching is done thanks to Iodine/Potassium solution. A second photolithographic process is performed through a second mask that allows the realization of the chip's pads, built with a thin layer of titanium of 100 nm and an aluminum layer of 700 nm (see Figure ).
The custom hardware
To control the experiments, a custom electronic board, connected to a Linux-embedded system, has been realized. To provide the current density of 10
8A/
cm2, needed for the activation of the electromigration process, the circuit must be able to supply a current of at least 50 mA, because of the dimensions of realized geometries. The front-end must also be able to measure the real time current flowing in the wire, to evaluate resistance variations, from hundreds of mA (when the current is high and the break is not yet created) to some pA (for measuring the tunnel current inside the nanogap). The block diagram of the system is showed in Figure , where it is possible to see that the gold probe is connected to the circuit that receives the signal from the embedded Linux (analog input) and generates the desired current for inducing the electromigration. The current in the probe, and so the resistance, is measured. These values are sent to the embedded Linux in which are implemented the algorithms for controlling the current to inject into the probe. In this way, the fabrication of the nanogap is checked, avoiding thermal runaway that can create too large breaks [
18-
22], not useful for the desired applications.
The custom system consists of a driver board that hosts a digital-to-analog converter that realizes the input voltage waveform; a measurement board composed by a transimpedance amplifier with variable feedback, in order to measure a wide range of currents, and an analog-to-digital converter; a switch board to allow the connection of external instrumentation; a digital board that provides electrical power supplies and the bus connection between all boards. The embedded Linux system is built with a real-time custom kernel, so that the electronic components of the other boards are driven in a deterministic way. We have developed a wireless connection between this board and a host computer for sending the experimental data. This solution allows also the use of the system in chambers where a wire interconnection can create difficulties.
The custom algorithm
Nanogaps are produced using a specific algorithm, customized for an optimal nanogap fabrication and with several goals: current management, feedback controlled breaking, temperature control for avoiding the thermal runaway. A simplified schematic of the process flow is reported in Figure . The software controls the voltage
Vbias applied to the probe and stops when the resistance exceeds the value of 13
kΩ, that means that the nanogap is produced. In fact this resistance value is about the inverse of the quantum conductance 2
e2/
h = 77.6
μS and represents the conductance of a single atom of gold placed between two electrodes [
23]. As it is possible to see always in Figure , there are two feedback mechanisms. The first one performs an absolute control over the initial resistance
R0, and when (
R -
R0)/
R0 >0.02 the
Vbias is set to the 85% of its value; this is done for controlling the temperature of the wire and for preventing melting and surface tension effects, that can be the cause of much larger gaps and gold island formation [
18-
20]. The second mechanism performs a check of the resistance value relative to a circular buffer of n samples, when (
R -
Ri)/
Ri >0.01, with
Ri the average of the last
i samples, the
Vbias is decreased to the 95% of its value. The second feedback mechanism performs a check of the resistance in order to control the electromigration effect: the increase of the resistance is no longer linear as in the Joule Heating, but exponential, and we need to promptly react to stop the phenomenon. We also noted in many experiments performed by us, that, after the activation of the feedback, the resistance value decreases due to the acceleration of the grain growth by Joule heating of the wire [
24]. In this case, the first mechanism becomes useless and then the process is controlled only by the second one. Figure shows the probe resistance as a function of the voltage applied: when the temperature increases, the resistance of the wire tends to increase too.
The higher resistance causes a current reduction, but increasing the voltage in this case creates a mechanism by which the current flow tends to be constant. In fact the first increase in resistance (linear growth) is only due to this heat effect, but, when the temperature reaches high values and the current density is near to 108 A/cm2, the electromigration starts (exponential growth) and the structure begins to change.
The experimental time of electromigration has been estimated to be about 50-60 min. All the experimental data obtained by the electromigration tests are stored in an internal database in the Linux board, but on the host computer too, that collects information through the wireless connection. This makes possible to generate very accurate statistics. To evaluate the outcome of our experiments, we performed the analysis of gap widths with a FESEM microscope. Example of a fabricated nanogap is shown in Figure .