Linearscaling techniques for first principles calculations of stationary and dynamic systems
Access Status
Authors
Date
2010Supervisor
Type
Education Level
Metadata
Show full item recordAbstract
First principles calculations can be a computationally intensive task when studying large systems. Linearscaling methods must be employed to find the electronic structure of systems consisting of thousands of atoms and greater. The goal of this thesis is to combine the linearscaling divideandconquer (D&C) method with the linearscaling capabilities of the SIESTA (Spanish Initiative for Electronic Simulations with Thousands of Atoms) density functional theory (DFT) methodology and present this union as a viable approach to largescale first principles calculations. In particular, the density matrix version of the D&C method is implemented into the SIESTA package. This implementation can accommodate high quality calculations consisting of atom numbers in the tens of thousands using moderate computing resources. Low quality calculations have been tested up to half million atoms using reasonably sized computing resources.The D&C method is extended to better handle atomic dynamics simulations. First, by alleviating issues caused by discontinuities in the potential energy surface, with the application of a switching function on the Hamiltonian and overlap matrices. This allows for a smooth potential energy surface to be generated. The switching function has the additional benefit of accelerating the selfconsistent field (SCF) process. Secondly, the D&C frozen density matrix (FDM) is modified to allow for improved charge transfer between the active and constrained regions of the system. This modification is found to reduce both the number of SCF iterations required for selfconsistency and the number of relaxation steps in a local geometry optimisation. The D&C paradigm is applied to the realtime approach of timedependent density functional theory (TDDFT). The method is tested on a linear alkane molecule with varying levels of success.Divergences in the induced dipole moment occur when the external excitation field is aligned parallel to the axis of the molecule. The method succeeds in producing accurate dipole moments when the external field is aligned perpendicular to the molecule. Various techniques are tested to improve the proposed method. Finally, the performance and effectiveness of the current D&C implementation is evaluated by studying three current systems. The first two systems consist of two different DNA sequences and the last system is the large ZIF100 zeolitic imidazolate framework (ZIF).
Department
Collections
Related items
Showing items related by title, author, creator and subject.

Cankurtaran, Burak O.; Gale, Julian; Ford, M. (2008)The densitymatrix divideandconquer technique for the solution of KohnSham density functional theory has been implemented within the framework of the SIESTA methodology. Implementation details are provided where the ...

Dovesi, R.; Orlando, R.; Erba, A.; ZicovichWilson, C.; Civalleri, B.; Casassa, S.; Maschio, L.; Ferrabone, M.; De La Pierre, M; D'Arco, P.; Noël, Y.; Causà, M.; Rérat, M.; Kirtman, B. (2014)The capabilities of the CRYSTAL14 program are presented, and the improvements made with respect to the previous CRYSTAL09 version discussed. CRYSTAL14 is an ab initio code that uses a Gaussiantype basis set: both ...

Blanchard, M.; Wright, Kate; Gale, Julian; Catlow, R. (2007)We have performed firstprinciples calculations based on density functional theory (DFT), to model interactions of arsenic with FeS2 pyrite. An understanding of pyrite reactivity on the atomic scale is needed to predict ...