Every week several articles are published where ATK was used to perform the calculations - in fact, in 2011 there were over 150 publications based on ATK. You can read the abstracts of all these papers (as well as from earlier years - there are over 450 in total!) in our Publication List.

Although we are often very busy implementing new features in ATK and supporting users, once in a while we at QuantumWise also find time to make some serious studies of interesting atomic-scale transport systems. Most recently, we have looked at the interface between graphene and nickel, and these results have now been published in Physical Review B 85 165442 (2012). You can also find a preview copy in arXiv (arXiv:1202.0171v2 [cond-mat.mes-hall]).

For this article we performed DFT calculations of electron transport across a nickel-graphene interface. Four different geometries (see figure below) were considered, where the contact area, graphene and nickel surface orientations, and the passivation of the terminating graphene edge were varied. We found covalent bond formation between the graphene layer and the nickel surface, in agreement with other theoretical studies. We then calculated the energy-dependent electron transmission for the four systems and find that the systems have very similar edge contact resistance, which is furthermore independent of the contact area between nickel and graphene.

Calculations for graphene bonded with model surfaces show that the results are generic for strongly bonded graphene surfaces, and the minimum attainable edge contact resistance is twice the ideal edge quantum contact resistance of graphene.

The computed contact resistance is in excellent quantitative agreement with recent experimental results by K. Nagashio et al. [Appl. Phys. Lett. 97, 143514 (2010)].

In addition to drawing some important conclusions about the metal-graphene interfaces, the article also showcases a new feature in ATK 12.2: the current density, seen in the plot below, which allows you to visualize how the current propagates in a nanoscale junction.

Top view (B-C plane) and side view (A-C plane) of the four systems investigated in this paper. The transport direction is along the C direction. (a) Zigzag edge graphene on top of a Ni(100) surface with 4 Å binding overlap. (b) Similar to (a) but with 8 Å overlap. (c) Similar to (a) but with a hydrogen-passivated terminal edge. (d) Armchair edge graphene on top of a Ni(110) surface with 4 Å overlap, and hydrogen-passivated terminal edge. The inset (red curve) in each of the figures shows the average electrostatic potential in the vacuum region.

Contour plot of the C-component of the linear response current density (arbitrary units) in the A-C plane at 0.05 V bias, averaged in the B-direction.