Numerical solution of the Hartree-Fock equation by multilevel tensor-structured methods [Elektronische Ressource] / vorgelegt von Venera Khoromskaia

Sujets

• Mathematics

Informations

• Publié par: technische_universitat_berlin
• Publié le: 01 janvier 2010
• Nombre de lectures: 8
• Langue: English
• Poids de l'ouvrage: 5 Mo

Extrait

Numerical solution of the Hartree-Fock equation
by multilevel tensor-structured methods
vorgelegt von Diplom-Physikerin
Venera Khoromskaia
Stadt Kazan, Russland
Von der Fakult¨at II - Mathematik und Naturwissenschaften
der Technischen Universit¨at Berlin
zur Erlangung des akademisches Grades
Doktor der Naturwissenschaften
Dr.rer.nat.
genehmigte Dissertation
Promotionausschuss:
Vorsitzender: Prof. Dr. J. Blath
Berichter/Gutachter: Prof. Dr. Reinhold Schneider
Berichter/Gutachter: Prof. Dr. Dr. h.c. Wolfgang Hackbusch
zus¨atzlicher Gutachter: Prof. Dr. Eugene Tyrtyshnikov
Tag der mu¨ndlichen Pru¨fung: 10 December 2010
Berlin 2011
D 832Acknowledgements
I would like to express my gratitude to Prof. Dr. Reinhold Schneider for su-
pervising my PhD project, forfruitful discussions and his friendly encouragement
during the work on this project. His interest and expertise in the research topics
related tothe Hartree-Fockequation and the density functional theorymotivated
muchoftherecentprogressinthealgebraictensormethodsinelectronicstructure
calculations.
I would like to thank Prof. Dr. Dr. h.c. Wolfgang Hackbusch for valuable
discussions and excellent conditions for performing research at the Max-Planck-
Institute for Mathematics in the Sciences.
This work is done due to fruitful collaboration with Dr. Heinz-Ju¨rgen Flad. I
appreciateverymuchhiskindencouragement andsupportofmyworkandalways
beneficial and stimulating discussions owing to his thorough expertise in modern
quantum chemistry.
ProfessionalassistanceofPDDrSci. BorisKhoromskijintheresearchontensor
numerical methods helped me to make an active start in a new field. I acknowl-
edge him for the statement of some research problems and for proofreading the
manuscript.
I kindly acknowledge Prof. Eugene Tyrtyshnikov, Prof. Dr. Christian Lu-
bich and Prof. Dr. Lars Grasedyck, for their interest to my work and for the
opportunity to give talks at recent conferences and seminars on tensor methods.
IamverymuchappreciativetoProf. Dr. IvanGavrilyukforhisencouragement
and interest to my work.
I would like to thank my colleagues at the Max-Planck-Institute in Leipzig,
Dr. Ronald Kriemann, Dr. Jan Schneider, Dr. Kishore Kumar Naraparaju,
Dr. Sambasiva Rao Chinnamsetty, Dr. Hongjun Luo, Dr. Thomas Blesgen, Dr.
Lehel Banjai, Dipl.-Math. Konrad Kaltenbach, Dipl.-Math. Stephan Schwinger,
Dipl.-Math. Florian Drechsler and colleagues at the TU Berlin, Prof. Dr. Harry
Yserentant, Dipl.-Math. Fritz Kru¨ger and Dipl.-Math. Andr´e Ushmaev for inter-
esting discussions. I would like to acknowledge the colleagues from the Institute
of Numerical Mathematics of the Russian Academy of Science in Moscow, Dr.
Ivan Oseledets and Dr. Dmitrij Savostyanov for stimulating discussions.
I would like to thank Prof. Vikram Gavini (University of Michigan) for pro-
ductive collaboration and for the data on electron density of large Aluminium
clusters.
Kind assistance of the librarians Mrs. Britta Schneemann and Mrs. Katarzyna
Baier was very helpful during my work. I would like to thank cordially the
3secretaries at the Max-Planck-Institute and TU Berlin, Mrs. Valeria Hu¨nniger
and Mrs. Susan Kosub for their helpful technical support.
4Numerische L¨osung der Hartree-Fock-Gleichung
mit mehrstufigen Tensor-strukturierten Verfahren
Venera Khoromskaia
Abstract der PhD Dissertation
DiegenaueL¨osungderHartree-Fock-Gleichung(HFG),dieeinnichtlinearesEigen-
3wertprobleminR darstellt, istinfolgedernichtlokalen Intergraltransformationen
und der scharfen Peaks in der Elektronendichte und den Moleku¨lorbitalen eine
herausfordernde numerische Aufgabe. Aufgrund der nichtlinearen Abh¨angigkeit
der Hamilton-Matrix von den Eigenvektoren, ist das Problem nur iterativ l¨osbar.
Die traditionelle L¨osung der HFG basiert auf einer analytischen Berechnung der
3auftretenden Faltungsintegrale im R mit Hilfe von dem Problem angepassten
Basen (so genannte Zwei-Elektron-Integrale). Die inh¨arenten Grenzen dieses
Konzepts werden wegen der starken Abh¨angigkeit der numerischen Effizienz von
der Gr¨oße und den Eigenschaften der analytisch separablen Basis sichtbar.
IndieserDissertationwurdenneuegitter-basiertemehrstufigeTensor-strukturierte
Verfahren entwickelt und anhand der numerischen L¨osung der HFG getestet.
Diese Methoden beinhalten effiziente Algorithmen zur Darstellung diskretisierter
3Funktionen und Operatoren in R durch strukturierte Tensoren in den kanonis-
chen,Tucker-undkombiniertenTensorformatenmiteinerkontrollierbarenGenauigkeit
sowie schnelle entsprechenden Operationen fu¨r multilineare Tensoren. Insbeson-
dere wird die beschleunigte Mehrgitter-Rang-Reduktion des Tensors vorgestellt,
die auf der reduzierten Singul¨arwertzerlegung h¨oherer Ordnung basiert.
Der Kern der Anwendung dieser Verfahren fu¨rdie L¨osung der HFG ist die Ver-
wendung strukturierter Tensoren zur genauen Berechnung der Elektronendichte
3und der nichtlinearen Hartree- und (nichtlokalen) Austauschoperatoren in R ,
die in jedem Iterationsschritt auf einer Reihenfolge von n×n×n kartesischen
Gittern darstellt wurden. Somit wurden die entsprechenden sechs-dimensionalen
Integrationen durch multilineare algebraische Operationen wie das Skalar- und
Hadamardprodukt, die dreidimensionale Faltungstransformation und die Rang-
Reduktionfu¨rTensorendritterOrdnungersetzt,dieann¨aherndmitO(n)-Komplexit¨at
implementiert wurden, wobei n die eindimensionale Gittergr¨oße ist. Daher ist
der wesentliche Vorteil unserer Tensor-strukturierten Verfahren, dass die gitter-
dbasierte Berechnung von Integraloperatoren inR , d≥3, lineare Komplexit¨at in
n hat. Man beachte, dass im Sinne der u¨bliche Absch¨atzung mittels des Gitter-
3volumens N =n die Operationen mit strukturierten Tensoren eine sublinearevol
1/3
Komplexit¨at haben, O(N ).vol
Das vorgestellte ”‘grey-box”’-Schema zur L¨osung der HFG erfordert keine an-
alytischen Vorberechnungen der Zwei-Elektron-Integrale. Weiterhin ist dieses
5Schema sehr flexibel hinsichtlich der Wahl der gitter-orientierten Basisfunktio-
nen.
Numerische Berechnungen am Beispiel des “all electron” Falls fu¨r H O, CH2 4
undC H unddesPseudopotentialfalls fu¨rCH OHandC H OHMoleku¨le zeigen2 6 3 2 5
die geforderte hohe Genauigkeit.
Die Tensor-strukturierten Verfahren k¨onnen auch zur L¨osung der Kohn-Sham-
Gleichung angewandt werden, indem anstelle einer problem-unabh¨angigen Basis,
3wie die der ebenen Wellen oder einer großen Anzahl finiter Elemente imR , eine
geringe Anzahl problem-orientierter rang-strukturierter algebraischer Basisfunk-
tionen verwendet werden, die auf einem Tensorgitter dargestellt sind.
6Numerical solution of the Hartree-Fock equation
by multilevel tensor-structured methods
Venera Khoromskaia
Abstract of PhD dissertation
An accurate solution of the Hartree-Fock equation, a nonlinear eigenvalue prob-
3lem inR , is a challenging numerical task due to the presence of nonlocal integral
transforms and strong cusps in the electron density and molecular orbitals. In
view of the nonlinear dependence of the Hamiltonian matrix on the eigenvec-
tors, thisproblem can onlybe solved iteratively, byself-consistent field iterations.
Traditionally, the solution of the Hartree-Fock equation is based on rigorous an-
3alytical precomputation of the arising convolution type integrals in R in the
naturally separable basis (so-called two-electron integrals). Inherent limitations
of this concept are evident because of the strong dependence of the numerical
efficiency on the size and approximation quality of the problem adapted basis
sets.
In this dissertation, novel grid-based multilevel tensor-structured methods are
developedandtestedbyanumericalsolutionoftheHartree-Fockequation. These
methodsinclude efficient algorithmsforthe low-rankrepresentation ofdiscretized
3functions and operators inR , in the canonical, Tucker and mixed tensor formats
with a controllable accuracy, and fast procedures for the corresponding multilin-
ear tensor operations. In particular, a novel multigrid accelerated tensor rank
reduction method is introduced, based on the reduced higher order singular value
decomposition.
The core of our approach to the solution of the Hartree-Fock equation is the
accurate tensor-structured computation of the electron density and the nonlinear
3Hartree and the (nonlocal) exchange operators in R , discretized on a sequence
of n×n×n Cartesian grids, at every step of nonlinear iterations. Hence, the
corresponding six-dimensional integrations are replaced by multilinear algebra
operations such as the scalar and Hadamard products, the 3D convolution trans-
form, and the rank truncation for 3rd order tensors, which are implemented with
an almost O(n)-complexity, where n is the univariate grid size. In this way, the
basic advantage of our tensor-structured methods is the grid-based evaluation of
dintegral operators inR , d≥ 3, with linear complexity in n. Note that in terms
3of usual estimation by volume size N = n , the tensor-structured operationsvol
1/3
are of sublinear complexity, O(N ).
vol
The proposed “grey-box” scheme for the