A New Spin On Atomic Logic
Lead Research Organisation:
King's College London
Department Name: Physics
Abstract
Abstracts are not currently available in GtR for all funded research. This is normally because the abstract was not required at the time of proposal submission, but may be because it included sensitive information such as personal details.
| Description | As a theoretical group, we have been working on a development of a new efficient method of calculating the Scanning Tunnelling Microscopy (STM) images. STM is the main tool used by the Nottingham experimental group in this project and it is essential that a proper theoretical support is provided to adequately interpret their images and workflow. In STM an atomically sharp tip is approached to the sample surface (with the distances between the two being between 2-10 Angstrems); when the bias is applied between the tip and surface, a current flows that is measured. By placing the tip at different lateral positions around the sample surface (its vertical position depends on the experimental mode used), a 2D image of the surface is made that in many cases has molecular and even atomic resolution. Currently, practically one method is being widely used for theoretical interpretation of experimental images, that is due to Tersoff and Hamann (TH). This method is cheap (which explains why it is so popular), as the whole STM image can actually be obtained from a single electronic structure calculation (made usually using a Density Functional Theory based code) of the surface under investigation, without using any tip model. Hence, the method provides a kind of "average" image over a range of possible tips. However, this method is highly approximate: it is not supposed to work at small tip-sample distances (if used in experiment) or at high voltages; basically, it is only expected to work in the limit of zero bias. As a consequence, the TH method gives often quantitatively, but also in many cases, qualitatively wrong results. The "exact" method of calculating the current roots in the Landauer formula derived using the Non-Equilibrium Green's Function (NEGF) method. This approach however has two main problems: (i) it requires a precise atomistic model of the tip; (ii) it is extremely computationally expensive as very costly current calculation should be done on a grid above the surface, basically hundreds or even thousands of such calculations are required. The latter point makes this method impractical. In this project we have been developing a new theoretical tool for calculating the STM images. The expectation is that it will be much more efficient computationally compared to the full NEGF method; however, it will bring much more precision in terms of the range of the tip-sample distances and the bias compared to the TH method. The tool is based on the exact Landauer formula for the tunnelling current, but makes a number of simplifications concerning the tip. As a result, each calculation required for the given tip position above the sample surface is greatly simplified, and hence the second difficulty mentioned above of the NEGF method is essentially relaxed, if not eliminated. With this method, we should be able to provide more realistic STM images at a wider range of tip-sample distance, voltages and even at different temperatures (the TH method can only be applied at zero temperature). Currently: (i) the theoretical part of the method has been finished; (ii) it has been implemented in a code; (iii) a workflow has been created that enables one to run a CP2K code (relation to CP2K-UK award EP/K038222/1) in a script and extract the obtained wavefunctions and calculate the current; (iv) a workflow is being created for high throughput automatic simulations of STM images. We have also run a number of tests on the Au(111) surface with encouraging results. The method is versatile enough to include the spin into the consideration, i.e. the tip model can be, if desired, be spin-sensitive leading to different spin-up and spin-down current channels. The method requires a few parameters associated with the tip model, these can be obtained by fitting the current vs. tip height data obtained from experiment. |
| Exploitation Route | The code will become freely available (under the GNU General Public License (GPL)) for use by a wider scientific community. It will be placed on GitHub with full description, examples and a tutorial. Of course, when we are ready, we expect a number of scientific publications as well, both on the theory and applications to various systems including quantum corals. |
| Sectors | Aerospace, Defence and Marine,Digital/Communication/Information Technologies (including Software),Education,Electronics,Energy,Manufacturing, including Industrial Biotechology |
| Description | Nottingham |
| Organisation | University of Nottingham |
| Country | United Kingdom |
| Sector | Academic/University |
| PI Contribution | our theory provides an input to experiments performed by this partner |
| Collaborator Contribution | experimental results from the partner serve as the basis for our theory |
| Impact | Two papers are ready for submission. The partner uses chemistry synthesis methods, we use physics methods. |
| Start Year | 2014 |
| Title | CP2K |
| Description | CP2K is a program to perform atomistic and molecular simulations of solid state, liquid, molecular, and biological systems. It provides a general framework for different methods such as e.g., density functional theory (DFT) using a mixed Gaussian and plane waves approach (GPW) and classical pair and many-body potentials. |
| Type Of Technology | Software |
| Year Produced | 2014 |
| Open Source License? | Yes |
| Impact | New functionality added and the website improved. |
| URL | http://cp2k.org |