Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA numerical method to efficiently calculate the transport properties of large systems: an algorithm optimized for sparse linear solvers
194
0
0.0
(
0
)
Added by
Tatiane Pereira Dos Santos
Publication date
2018
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We present a self-contained description of the wave-function matching (WFM) method to calculate electronic quantum transport properties of nanostructures using the Landauer-Buttiker approach. The method is based on a partition of the system between a central region (conductor) containing $N_S$ sites and an asymptotic region (leads) characterized by $N_P$ open channels. The two subsystems are linearly coupled and solved simultaneously using an efficient sparse linear solver. Invoking the sparsity of the Hamiltonian matrix representation of the central region, we show that the number of operations required by the WFM method in conductance calculations scales with $sim N_Stimes N_P$ for large $N_S$.
rate research
Read More
As the characteristic lengths of advanced electronic devices are approaching the atomic scale, ab initio simulation method, with fully consideration of quantum mechanical effects, becomes essential to study the quantum transport phenomenon in them. However, current widely used non-equilibrium Greens function (NEGF) approach is based on atomic basis set, which usually can only study small system with less than 1000 atoms in practice. Here we present a large-scale quantum transport simulation method using plane waves basis, based on the previously developed plane wave approach (Phys. Rev. B 72, 045417). By applying several high-efficiency parallel algorithms, such as linear-scale ground-state density function theory (DFT) algorithm, folded spectrum method, and filtering technique, we demonstrate that our new method can simulate the system with several thousands of atoms. We also use this method to study several nanowires with about 4000 copper atoms, and show how the shape and point defect affect the transport properties of them. Such quantum simulation method will be useful to investigate and design nanoscale devices, especially the on-die interconnects.
Vibrational spectra can be computed without storing full-dimensional vectors by using low-rank sum-of-products (SOP) basis functions. We introduce symmetry constraints in the SOP basis functions to make it possible to separately calculate states in different symmetry subgroups. This is done using a power method to compute eigenvalues and an alternating least squares method to optimize basis functions. Owing to the fact that the power method favours the convergence of the lowest states, one must be careful not to exclude basis functions of some symmetries. Exploiting symmetry facilitates making assignments and improves the accuracy. The method is applied to the acetonitrile molecule.
We present an approach to calculate total energies of nanoclusters based on first principles estimates. For very large clusters the total energy can be separated into surface, edge and corner energies, in addition to bulk contributions. Using this separation and estimating these with direct, first principles calculations, together with the relevant chemical potentials, we have calculated the total energies of Cu and CdSe tetrahedrons containing a large number of atoms. In our work we consider polyhedral clusters so that in addition our work provides direct information on relaxation. For Cu the effects are very small and the clusters vary uniformly from very small to very large sizes. For CdSe there are important variations in surface and edge structures for specific sizes; nevertheless, the approach can be used to extrapolate to large non-stoichiometric clusters with polar surfaces.
Density Functional Theory (DFT) has become the quasi-standard for ab-initio simulations for a wide range of applications. While the intrinsic cubic scaling of DFT was for a long time limiting the accessible system size to some hundred atoms, the recent progress with respect to linear scaling DFT methods has allowed to tackle problems that are larger by many orders of magnitudes. However, as these linear scaling methods were developed for insulators, they cannot, in general, be straightforwardly applied to metals, as a finite temperature is needed to ensure locality of the density matrix. In this paper we show that, once finite electronic temperature is employed, the linear scaling version of the BigDFT code is able to exploit this locality to provide a computational treatment that scales linearly with respect to the number of atoms of a metallic system. We provide prototype examples based on bulk Tungsten, which plays a key role in finding safe and long-lasting materials for Fusion Reactors. We believe that such an approach might help in opening the path towards novel approaches for investigating the electronic structure of such materials, in particular when large supercells are required.
We present an implementation of the steady state Keldysh approach in a Greens function multiple scattering scheme to calculate the non-equilibrium spin density. This density is used to obtain the spin transfer torque in junctions showing the magnetoresistance effect. We use our implementation to study the spin transfer torque in metallic Co/Cu/Co junctions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
