Application of the Random Phase Approximation to complex problems in materials science

Die vorliegende Dissertation widmet sich der Bewertung und Anwendung der Random Phase Approximation (RPA) im Rahmen des Adiabatic-Connection Fluctuation-Dissipation Theorems (ACFDT) auf Problemstellungen der Festkörperphysik. Im Theorie- und Methodenteil wird ein Überblick über Dichte Funktional Theorie (DFT) und das ACFD Theorem gegeben. Weiters beinhaltet dieser Teil eine einführende Diskussion des Vielteilchen Problems wie auch eine Beschreiubung der Implementierung der Random Phase Approximation in dem Softwarepaket Vienna Ab-initio Simulation Package (VASP). Im zweiten Teil der Disseration wird die Genauigkeit der RPA uberprüft und mit Resultaten anderer Funktionale verglichen. Experimentelle Ergebnisse werden in Bezug auf Phononen Nullpunkt-Vibrationsenergien, welche ab-initio auf DFT Ebene berechnet wurden, korrigiert. RPA liefert eine sehr genaue Beschreibung aller Bindungsarten und ist daher ein vielversprechender Kandidat für komplexere Problemstellungen der Festkörperphysik. Als erstes Beispiel untersuchen wir die Wechselwirkung von Wassermolekülen mit Kohlenstoffverbindungen an Hand zweier Fälle: Wasser auf Benzen und Wasser auf einer Graphenoberfläche. Unsere Ergebnisse zeigen gute Ubereinstimmung mit einer weiteren korrelierten Methode: Diffusion Monte Carlo(DMC). Wir glauben daher, dass unsere gefundenen Werte als Richtwerte für die Entwicklung weiterer DFT Funktionale zur Beschreibung von Wasser-Kohlenstoff Wechselwirkungen dienen können. Ein Kernstück dieser Dissertation ist sicherlich die erfolgreiche Anwendung der RPA auf das bis dahin ungelöste CO Adsorptions Rätsel. Wir diskutieren die Adsorption eines Kohlenmonoxidmoleküls auf Cu, Ru, Rh, Pd und Pt. RPA ist gegenwärtig die einzige ab-initio Methode, welche sowohl die Adsorptionsenergie, als auch die Oberflächenenergie genau beschreibt und weiters die richtige Position des Moleküls auf der Metalloberfläche für jeden einzelnen Fall korrekt wiedergibt.This thesis is devoted to the assessment and application of the random phase approximation (RPA) in the adiabatic-connection fluctuation-dissipation (ACFD) framework in solid state physics. The first part presents a review of density functional theory (DFT) and the ACFD theorem in the RPA. This includes an introduction to the many-body problem as well as a description of the implementation of the RPA in the Vienna Ab-initio Simulation Package (VASP). In the results part, the quality of the RPA is assessed and its performance compared to three (beyond) DFT functionals. The experimental values are corrected for the effect of phonon zero-point vibrational energies which were calculated at the DFT level from ab-initio. We find that the RPA describes all bonding situations very accurately, making it a promising candidate for more complex problems in solid state physics. In light of these findings, we investigate the carbon-water interaction in two specific cases: the adsorption of water on benzene and the adsorption of water on a graphene layer. We compare our results to a different correlated method: diffusion Monte Carlo (DMC). We find very good agreement and thus believe that our values can serve as a benchmark for the development of other DFT functionals to treat water-carbon interfaces. The highlight of this thesis is the successful application of the RPA to the long-standing and (at DFT level) unsolved CO adsorption puzzle. We show results for CO adsorption on Cu, late 4d metals and Pt. RPA is at present the only ab-initio method that describes adsorption and surface energies accurately at the same time and predicts the correct adsorption site in every single case

Warming Up Density Functional Theory

Density functional theory (DFT) has become the most popular approach to electronic structure across disciplines, especially in material and chemical sciences. Last year, at least 30,000 papers used DFT to make useful predictions or give insight into an enormous diversity of scientific problems, ranging from battery development to solar cell efficiency and far beyond. The success of this field has been driven by usefully accurate approximations based on known exact conditions and careful testing and validation. In the last decade, applications of DFT in a new area, warm dense matter, have exploded. DFT is revolutionizing simulations of warm dense matter including applications in controlled fusion, planetary interiors, and other areas of high energy density physics. Over the past decade or so, molecular dynamics calculations driven by modern density functional theory have played a crucial role in bringing chemical realism to these applications, often (but not always) with excellent agreement with experiment. This chapter summarizes recent work from our group on density functional theory at non-zero temperatures, which we call thermal DFT. We explain the relevance of this work in the context of warm dense matter, and the importance of quantum chemistry to this regime. We illustrate many basic concepts on a simple model system, the asymmetric Hubbard dimer

Communication: Light driven remote control of microgels’ size in the presence of photosensitive surfactant: Complete phase diagram

Here we report on a light triggered remote control of microgel size in the presence of photosensitive surfactant. The hydrophobic tail of the cationic surfactant contains azobenzene group that undergoes a reversible photo-isomerization reaction from a trans- to a cis-state accompanied by a change in the hydrophobicity of the surfactant. We have investigated light assisted behaviour and the complex formation of the microgels with azobenzene containing surfactant over the broad concentrational range starting far below and exceeding several times of the critical micelle concentration (CMC). At small surfactant concentration in solution (far below CMC), the surfactant in the trans-state accommodates within the microgel causing its compaction, while the cis-isomer desorbs out of microgel resulting in its swelling. The process of the microgel size change can be described as swelling on UV irradiation (trans-cis isomerization) and shrinking on irradiation with blue light (cis-trans isomerization). However, at the surfactant concentrations larger than CMC, the opposite behaviour is observed: the microgel swells on blue irradiation and shrinks during exposure to UV light. We explain this behaviour theoretically taking into account isomer dependent micellization of surfactant within the microgels

Random-phase approximation and its applications in computational chemistry and materials science

The random-phase approximation (RPA) as an approach for computing the electronic correlation energy is reviewed. After a brief account of its basic concept and historical development, the paper is devoted to the theoretical formulations of RPA, and its applications to realistic systems. With several illustrating applications, we discuss the implications of RPA for computational chemistry and materials science. The computational cost of RPA is also addressed which is critical for its widespread use in future applications. In addition, current correction schemes going beyond RPA and directions of further development will be discussed.Comment: 25 pages, 11 figures, published online in J. Mater. Sci. (2012

Carbon related defects in irradiated silicon revisited

Electronic structure calculations employing hybrid functionals are used to gain insight into the interaction of carbon (C) atoms, oxygen (O) interstitials, and self-interstitials in silicon (Si). We calculate the formation energies of the C related defects C(i)(Si(I)), C(i)O(i), C(i)C(s), and C(i)O(i)(Si(I)) with respect to the Fermi energy for all possible charge states. The C(i)(Si(I))(2+) state dominates in almost the whole Fermi energy range. The unpaired electron in the C(i)O(i)(+) state is mainly localized on the C interstitial so that spin polarization is able to lower the total energy. The three known atomic configurations of the C(i)C(s) pair are reproduced and it is demonstrated that hybrid functionals yield an improved energetic order for both the A and B-types as compared to previous theoretical studies. Different structures of the C(i)O(i)(Si(I)) cluster result for positive charge states in dramatically distinct electronic states around the Fermi energy and formation energies

CRYSTAL14: A program for the ab initio investigation of crystalline solids

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 Gaussian-type basis set: both pseudopotential and all-electron strategies are permitted; the latter is not much more expensive than the former up to the first-second transition metal rows of the periodic table. A variety of density functionals is available, including as an extreme case Hartree–Fock; hybrids of various nature (global, range-separated, double) can be used. In particular, a very efficient implementation of global hybrids, such as popular B3LYP and PBE0 prescriptions, allows for such calculations to be performed at relatively low computational cost. The program can treat on the same grounds zero-dimensional (molecules), one-dimensional (polymers), two-dimensional (slabs), as well as three-dimensional (3D; crystals) systems. No spurious 3D periodicity is required for low-dimensional systems as happens when plane-waves are used as a basis set. Symmetry is fully exploited at all steps of the calculation; this permits, for example, to investigate nanotubes of increasing radius at a nearly constant cost (better than linear scaling!) or to perform self-consistent-field (SCF) calculations on fullerenes as large as (10,10), with 6000 atoms, 84,000 atomic orbitals, and 20 SCF cycles, on a single core in one day. Three versions of the code exist, serial, parallel, and massive-parallel. In the second one, the most relevant matrices are duplicated, whereas in the third one the matrices in reciprocal space are distributed for diagonalization. All the relevant vectors are now dynamically allocated and deallocated after use, making CRYSTAL14 much more agile than the previous version, in which they were statically allocated.The program now fits more easily in low-memory machines (as many supercomputers nowadays are). CRYSTAL14 can be used on parallel machines up to a high number of cores (benchmarks up to 10,240 cores are documented) with good scalability, the main limitation remaining the diagonalization step. Many tensorial properties can be evaluated in a fully automated way by using a single input keyword: elastic, piezoelectric, photoelastic, dielectric, as well as first and second hyperpolarizabilies, electric field gradients, Born tensors and so forth. Many tools permit a complete analysis of the vibrational properties of crystalline compounds. The infrared and Raman intensities are now computed analytically and related spectra can be generated. Isotopic shifts are easily evaluated, frequencies of only a fragment of a large system computed and nuclear contribution to the dielectric tensor determined. New algorithms have been devised for the investigation of solid solutions and disordered systems. The topological analysis of the electron charge density, according to the Quantum Theory of Atoms in Molecules, is now incorporated in the code via the integrated merge of the TOPOND package. Electron correlation can be evaluated at the Möller–Plesset second-order level (namely MP2) and a set of double-hybrids are presently available via the integrated merge with the CRYSCOR program