Variational Monte Carlo study of the hydrogen electronic structure at ultra-high pressure

Die causa finalis der vorliegenden Arbeit ist das Verständnis des Phasendiagramms von Wasserstoff bei ultrahohen Drücken, welche von nichtleitendem H2 bis hin zu metallischem H reichen. Da die Voraussetzungen für ultrahohen Druck im Labor schwer zu schaffen sind, bilden Computersimulationen ein wichtiges alternatives Untersuchungsinstrument. Allerdings sind solche Berechnungen eine große Herausforderung. Eines der größten Probleme ist die genaue Auswertung des Born-Oppenheimer Potentials, welches sowohl für die nichtleitende als auch für die metallische Phase geeignet sein muss. Außerdem muss es die starken Korrelationen berücksichtigen, die durch die kovalenten H2 Bindungen und die eventuellen Phasenübergänge hervorgerufen werden. Auf dieses Problem haben unsere Anstrengungen abgezielt. Im Kontext von Variationellem Monte Carlo (VMC) ist die Shadow Wave Function (SWF) eine sehr vielversprechende Option. Aufgrund ihrer Flexibilität sowohl lokalisierte als auch delokalisierte Systeme zu beschreiben sowie ihrer Fähigkeit Korrelationen hoher Ordnung zu berücksichtigen, ist sie ein idealer Kandidat für unsere Zwecke. Unglücklicherweise bringt ihre Formulierung ein Vorzeichenproblem mit sich, was die Anwendbarkeit limitiert. Nichtsdestotrotz ist es möglich diese Schwierigkeit zu umgehen indem man die Knotenstruktur a priori festlegt. Durch diesen Formalismus waren wir in der Lage die Beschreibung der Elektronenstruktur von Wasserstoff signifikant zu verbessern, was eine sehr vielversprechende Perspektive bietet. Während dieser Forschung haben wir also die Natur des Vorzeichenproblems untersucht, das sich auf die SWF auswirkt, und dabei ein tieferes Verständnis seines Ursprungs erlangt. Die vorliegende Arbeit ist in vier Kapitel unterteilt. Das erste Kapitel führt VMC und die SWF mit besonderer Ausrichtung auf fermionische Systeme ein. Kapitel 2 skizziert die Literatur über das Phasendiagramm von Wasserstoff bei ultrahohem Druck. Das dritte Kapitel präsentiert die Implementierungen unseres VMC Programms und die erhaltenen Ergebnisse. Zum Abschluss fasst Kapitel 4 unsere Bestrebungen zur Lösung des zur SWF zugehörigen Vorzeichenproblems zusammen.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

A New Crystalline LiPON Electrolyte: Synthesis, Properties, and Electronic Structure

The new crystalline compound, Li2PO2N, was synthesized using high temperature solid state methods starting with a stoichiometric mixture of Li2O, P2O5, and P3N5. Its crystal structure was determined ab initio from powder X-ray diffraction. The compound crystallizes in the orthorhombic space group Cmc2(1) (# 36) with lattice constants a = 9.0692(4) angstrom, b = 53999(2) angstrom, and c = 4.6856(2) angstrom. The crystal structure of SD-Li2PO2N consists of parallel arrangements of anionic chains formed of corner sharing (PO2N2) tetrahedra. The chains are held together by Li+ cations. The structure of the synthesized material is similar to that predicted by Du and Holzwarth on the basis of first principles calculations (Phys. Rev. B 81,184106 (2010)). The compound is chemically and structurally stable in air up to 600 degrees C and in vacuum up to 1050 degrees C. The Arrhenius activation energy of SD-Li2PO2N in pressed pellet form was determined from electrochemical impedance spectroscopy measurements to be 0.6 eV, comparable to that of the glassy electrolyte LiPON developed at Oak Ridge National Laboratory. The minimum activation energies for Li ion vacancy and interstitial migrations are computed to be 0.4 eV and 0.8 eV, respectively. First principles calculations estimate the band gap of SD-Li2PO2N to be larger than 6 eV. (C) 2013 Elsevier B.V. All rights reserved.

Electronic structure of copper nitrides as a function of nitrogen content

he nitrogen content dependence of the electronic properties for copper nitride thin films with an atomic percentage of nitrogen ranging from 26 ± 2 to 33 ± 2 have been studied by means of optical (spectroscopic ellipsometry), thermoelectric (Seebeck), and electrical resistivity measurements. The optical spectra are consistent with direct optical transitions corresponding to the stoichiometric semiconductor Cu3N plus a free-carrier contribution, essentially independent of temperature, which can be tuned in accordance with the N-excess. Deviation of the N content from stoichiometry drives to significant decreases from − 5 to − 50 μV/K in the Seebeck coefficient and to large enhancements, from 10− 3 up to 10 Ω cm, in the electrical resistivity. Band structure and density of states calculations have been carried out on the basis of the density functional theory to account for the experimental results.

Effect of Sb incorporation on the electronic structure of InAs quantum dots

On the basis of optical characterization experiments and an eight band kp model, we have studied the effect of Sb incorporation on the electronic structure of InAs quantum dots (QDs). We have found that Sb incorporation in InAs QDs shifts the hole wave function to the center of the QD from the edges of the QD where it is otherwise pinned down by the effects of shear stress. The observed changes in the ground-state energy cannot merely be explained by a composition change upon Sb exposure but can be accounted for when the change in lateral size is taken into consideration. The Sb distribution inside the QDs produces distinctive changes in the density of states, particularly, in the separation between excitation shells. We find a 50% increase in the thermal escape activation energy compared with reference InAs quantum dots as well as an increment of the fundamental transition decay time with Sb incorporation. Furthermore, we find that Sb incorporation into quantum dots is strongly nonlinear with coverage, saturating at low doses. This suggests the existence of a solubility limit of the Sb incorporation into the quantum dots during growth.

Photochemically induced nuclear spin polarization in reaction centers of photosystem II observed by 13C-solid-state NMR reveals a strongly asymmetric electronic structure of the P680.+ primary donor chlorophyll

We report 13C magic angle spinning NMR observation of photochemically induced dynamic nuclear spin polarization (photo- CIDNP) in the reaction center (RC) of photosystem II (PS2). The light-enhanced NMR signals of the natural abundance 13C provide information on the electronic structure of the primary electron donor P680 (chlorophyll a molecules absorbing around 680 nm) and on the pz spin density pattern in its oxidized form, P680⨥. Most centerband signals can be attributed to a single chlorophyll a (Chl a) cofactor that has little interaction with other pigments. The chemical shift anisotropy of the most intense signals is characteristic for aromatic carbon atoms. The data reveal a pronounced asymmetry of the electronic spin density distribution within the P680⨥. PS2 shows only a single broad and intense emissive signal, which is assigned to both the C-10 and C-15 methine carbon atoms. The spin density appears shifted toward ring III. This shift is remarkable, because, for monomeric Chl a radical cations in solution, the region of highest spin density is around ring II. It leads to a first hypothesis as to how the planet can provide itself with the chemical potential to split water and generate an oxygen atmosphere using the Chl a macroaromatic cycle. A local electrostatic field close to ring III can polarize the electronic charge and associated spin density and increase the redox potential of P680 by stabilizing the highest occupied molecular orbital, without a major change of color. This field could be produced, e.g., by protonation of the keto group of ring V. Finally, the radical cation electronic structure in PS2 is different from that in the bacterial RC, which shows at least four emissive centerbands, indicating a symmetric spin density distribution over the entire bacteriochlorophyll macrocycle.

The CuA center of cytochrome-c oxidase: electronic structure and spectra of models compared to the properties of CuA domains.

The electronic structure and spectrum of several models of the binuclear metal site in soluble CuA domains of cytochrome-c oxidase have been calculated by the use of an extended version of the complete neglect of differential overlap/spectroscopic method. The experimental spectra have two strong transitions of nearly equal intensity around 500 nm and a near-IR transition close to 800 nm. The model that best reproduces these features consists of a dimer of two blue (type 1) copper centers, in which each Cu atom replaces the missing imidazole on the other Cu atom. Thus, both Cu atoms have one cysteine sulfur atom and one imidazole nitrogen atom as ligands, and there are no bridging ligands but a direct Cu-Cu bond. According to the calculations, the two strong bands in the visible region originate from exciton coupling of the dipoles of the two copper monomers, and the near-IR band is a charge-transfer transition between the two Cu atoms. The known amino acid sequence has been used to construct a molecular model of the CuA site by the use of a template and energy minimization. In this model, the two ligand cysteine residues are in one turn of an alpha-helix, whereas one ligand histidine is in a loop following this helix and the other one is in a beta-strand.

Nanocontacts: probing electronic structure under extreme uniaxial strains

Nanometer-sized metallic necks have the unique ability to sustain extreme uniaxial loads (about 20 times greater than the bulk material). We present an experimental and theoretical study of the electronic transport properties under such extreme conditions. Conductance measurements on gold and aluminum necks show a strikingly different behavior: While gold shows the expected conductance decrease with increasing elastic elongation of the neck, aluminum necks behave in the opposite way. We have performed first-principles electronic-structure calculations which reproduce this behavior, showing that it is an intrinsic property of the bulk band structure under high uniaxial strain.

Emergence of half-metallicity in suspended NiO chains: Ab initio electronic structure and quantum transport calculations

Contrary to the antiferromagnetic and insulating character of bulk NiO, one-dimensional chains of this material can become half metallic due to the lower coordination of their atoms. Here we present ab initio electronic structure and quantum transport calculations of ideal infinitely long NiO chains and of more realistic short ones suspended between Ni electrodes. While infinite chains are insulating, short suspended chains are half-metallic minority-spin conductors that displays very large magnetoresistance and a spin-valve behavior controlled by a single atom.

Electronic structure of gated graphene and graphene ribbons

We study the electronic structure of gated graphene sheets. We consider both infinite graphene and finite width ribbons. The effect of Coulomb interactions between the electrically injected carriers and the coupling to the external gate are computed self-consistently in the Hartree approximation. We compute the average density of extra carriers n2D, the number of occupied subbands, and the density profiles as a function of the gate potential Vg. We discuss quantum corrections to the classical capacitance and we calculate the threshold Vg above which semiconducting armchair ribbons conduct. We find that the ideal conductance of perfectly transmitting wide ribbons is proportional to the square root of the gate voltage.