Computer Lab

Unix Training

• Getting Started With Linux
• UNIX Tutorial

MATLAB Training

Hands-on tutorials by Instructor

• MATLAB: Getting Started
• MATLAB: Sparse Matrices

Hands-on Lab tutorials by MathWorks

• MATLAB Student Center
• Rapid code iteration using cells in the Editor

Reference

• MATLAB Brief List of Commands
• MATLAB Documentation
• Some Common MATLAB Programming Pitfalls and How to Avoid Them

Books and notes

• C. Moler: Numerical Computing with MATLAB (SIAM, Philadelphia, 2004). [PDF]
• UD Crash Course on Matlab

Implementation Tools

• Running MATLAB m-files in the background under Linux
• Multithreaded MATLAB performance on multicores

MATLAB Scripts

Electron density in nanowires using equilibrium density matrix

• electron_density.m
• Discretization of 1D Hamiltonian

DOS of disordered nanowire using eigenvalues and Anderson localization of eigenfunctions

• disordered_nanowire.m

Density of states using equilibrium retarded Green function

• dos_negf_closed.m computes DOS for finite 1D wire
• dos_negf_open.m computes DOS for finite 1D wire attached to a single reservoir
• graphene_dos.m computes DOS for a supercell of graphene with periodic boundary conditions

Subband structure of graphene nanoribbons using tight-binding models

• 8zgnr.m (pedestrian code for 8-ZGNR only)

Quantum transport in 1D nanowires using NEGF

• qt_1d.m (code to compute the conductance and total and local density of states of a 1D nanowire, with possible potential barriers or impurities, attached to two semi-infinite electrodes)

Tunneling magnetoresistance in 1D models

• tmr_1d.m (code to compute conductance of F/I/F junction as a function of angle between magnetizations using 1D tight-binding Hamiltonian)

Spin-transfer torque in 1D models

• stt_1d.m (code to compute in-plan torque component in F/I/F junction using 1D tight-binding Hamiltonian)

Quantum transport in graphene nanostructures using NEGF

• gnr_cond_recursive.m,bstruct.m, blocktosparse.m, sparsetoblock.m, h_zigzag.m, invnn.m, Self.m, (code to compute the conductance of a finite graphene nanoribbon attached to two semi-infinite graphene electrodes)

• M.-H. Liu and K. Richter, Efficient quantum transport simulation for bulk graphene heterojunctions, arXiv:1206.0266.

MATLAB functions

• matrix_exp.m (Exponential, or any other function with small changed in the code, of a Hermitian matrix)
• visual_graphene_H.m (For a given tight-binding Hamiltonian on the honeycomb lattice, function plots position of carbon atoms and draws blue lines to represent hoppings between them; red circles to represent on-site potential between them; and cyan lines to represent the periodic boundary conditions; it can be used to test if the tight-binding Hamiltonian of graphene is set correctly); This function calls another three function which should be placed in the same directory (or in the path): atomCoord.m, atomPosition.m, and constrainView.m
• self_energy.m (Self-energy of the semi-infinite ideal metallic lead modeled on the square tight-binding lattice - the code shows how to convert analytical formulas of the lead surface Green function into a working program)

First-principles calculations using GPAW

• Getting started with GPAW (structure and binding energies of simple molecules)
• Basics of GPAW calculations
• Band structure of Fe
• DOS of Fe
• Graphene nanoribbons

NEGF-DFT using GPAW

Zero-bias transmission function of single-molecule nanojunctions

• Pt-$\mathrm {H} _{2}$-Pt junction

I-V curve of magnetically ordered zigzag GNR

• D. A. Areshkin and B. K. Nikolić, I-V curve signatures of nonequilibrium-driven band gap collapse in magnetically ordered zigzag graphene nanoribbon two-terminal devices, Phys. Rev. B 79, 205430 (2009). [PDF]
• J. Chen, K. S. Thygesen, and K. W. Jacobsen, Phys. Rev. B 85, 155140 (2012). [PDF]