GPUQT: An efficient linear-scaling quantum transport code fully implemented on graphics processing units

被引:8
作者
Fan, Zheyong [1 ,2 ]
Vierimaa, Ville [2 ]
Harju, Ari [2 ]
机构
[1] Bohai Univ, Sch Math & Phys, Jinzhou, Peoples R China
[2] Aalto Univ, QTF Ctr Excellence, Dept Appl Phys, FI-00076 Aalto, Finland
基金
芬兰科学院; 中国国家自然科学基金;
关键词
Quantum transport; Linear-scaling; GPU acceleration; CONDUCTIVITY;
D O I
10.1016/j.cpc.2018.04.013
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present GPUQT, a quantum transport code fully implemented on graphics processing units. Using this code, one can obtain intrinsic electronic transport properties of large systems described by a real-space tight-binding Hamiltonian together with one or more types of disorder. The DC Kubo conductivity is represented as a time integral of the velocity auto-correlation or a time derivative of the mean square displacement. Linear scaling (with respect to the total number of orbitals in the system) computation time and memory usage are achieved by using various numerical techniques, including sparse matrix-vector multiplication, random phase approximation of trace, Chebyshev expansion of quantum evolution operator, and kernel polynomial method for quantum resolution operator. We describe the inputs and outputs of GPUQT and give a few examples to demonstrate its usage, paying attention to the interpretations of the results. Program summary Program Title: GPUQT Program Files doi: http://dx.doi.org/10.17632/xbf5kbkzx7.1 Licensing provisions: GPLv3 Programming language: CUDA Nature of problem: Obtain intrinsic electronic transport properties of large systems described by real-space tight-binding Hamiltonians. Solution method: The DC conductivity is represented as a time integral of the velocity auto-correlation (VAC) or a time derivative of the mean square displacement (MSD). The calculations achieve linear scaling (with respect to the number of orbitals in the system) computation time and memory usage by using various numerical techniques, including sparse matrix-vector multiplication, random phase approximation of trace, Chebyshev expansion of quantum evolution operator, and kernel polynomial method for quantum resolution operator. Restrictions: The number of orbitals is restricted to about 20 million due to the limited amount of device memory in current GPUs. (C) 2018 Elsevier B.V. All rights reserved.
引用
收藏
页码:113 / 120
页数:8
相关论文
共 21 条
[21]   Electrical transport in carbon nanotubes: Role of disorder and helical symmetries [J].
Triozon, F ;
Roche, S ;
Rubio, A ;
Mayou, D .
PHYSICAL REVIEW B, 2004, 69 (12)