A new approach to the electronic structure of two dimensional disordered alloys

We propose a novel method to calculate the electronic Density of States (DOS) of a two dimensional disordered binary alloy. The method is highly reliable and numerically efficient, and Short Range Order (SRO) correlations can be included with no extra computational cost. The approach devised rests on one dimensional calculations and is applied to very long stripes of finite width, the bulk regime being achieved with a relatively small number of chains in the disordered case. Our approach is exact for the pure case and predicts the correct DOS structure in important limits, such as the segregated, random, and ordered alloy regimes. We also suggest important extensions of the present work. © 1995.

Molecular and Electronic Structure of the Peptide Subunit of Geobacter sulfurreducens Conductive Phi from First Principles

The respiration of metal oxides by the bacterium Geobacter sulfurreducens requires the assembly of a small peptide (the GS pilin) into conductive filaments termed pili. We gained insights into the contribution of the GS pilin to the pilus conductivity by developing a homology model and performing molecular dynamics simulations of the pilin peptide in vacuo and in solution. The results were consistent with a predominantly helical peptide containing the conserved a-helix region required for pilin assembly but carrying a short carboxy-terminal random-coiled segment rather than the large globular head of other bacterial pilins. The electronic structure of the pain was also explored from first principles and revealed a biphasic charge distribution along the pilin and a low electronic HOMO-LUMO gap, even in a wet environment. The low electronic band gap was the result of strong electrostatic fields generated by the alignment of the peptide bond dipoles in the pilin's alpha-helix and by charges from ions in solution and amino acids in the protein. The electronic structure also revealed some level of orbital delocalization in regions of the pilin containing aromatic amino acids and in spatial regions of high resonance where the HOMO and LUMO states are, which could provide an optimal environment for the hopping of electrons under thermal fluctuations. Hence, the structural and electronic features of the pilin revealed in these studies support the notion of a pilin peptide environment optimized for electron conduction.

New insights on reaction pathway selectivity promoted by crown ether phase-transfer catalysis: Model ab initio calculations of nucleophilic fluorination

Crown ethers have the ability of solubilizing inorganic salts in apolar solvents and to promote chemical reactions by phase-transfer catalysis. However, details on how crown ethers catalyze ionic S(N)2 reactions and control selectivity are not well understood. In this work, we have used high level theoretical calculations to shed light on the details of phase-transfer catalysis mechanism of KF reaction with alkyl halides promoted by 18-crown-6. A complete analysis of the of the model reaction between KF(18-crown-6) and ethyl bromide reveals that the calculations can accurately predict the product ratio and the overall kinetics. Our results point out the importance of the K* ion and of the crown ether ring in determining product selectivity. While the K* ion favors the S(N)2 over the E2 anti pathway, the crown ether ring favors the S(N)2 over E2 syn route. The combination effects lead to a predicted 94% for the S(N)2 pathway in excellent agreement with the experimental value of 92%. A detailed analysis of the overall mechanism of the reaction under phase-transfer conditions also reveals that the KBr product generated in the nucleophilic fluorination acts as an inhibitor of the 18-crown-6 catalyst and it is responsible for the observed slow reaction rate. (C) 2012 Elsevier B.V. All rights reserved.

Ab initio studies of pristine and oxidized Cu3Au(100) and (111) surfaces

We review the previous literature and our recent work on first-principles studies of Cu3Au(100) and (111) surfaces, with focus on the segregation of atomic species to the surface at pristine conditions and in the presence of oxygen. In particular, the combined use of experimental and theoretical tools to achieve chemical identification at an atomic level of the surface species is emphasized and discussed.

ESTIMATING REACTION CONSTANTS BY AB INITIO MOLECULAR MODELING: A STUDY ON THE OXIDATION OF PHENOL TO CATECHOL AND HYDROQUINONE IN ADVANCED OXIDATION PROCESSES

Molecular modeling is growing as a research tool in Chemical Engineering studies, as can be seen by a simple research on the latest publications in the field. Molecular investigations retrieve information on properties often accessible only by expensive and time-consuming experimental techniques, such as those involved in the study of radical-based chain reactions. In this work, different quantum chemical techniques were used to study phenol oxidation by hydroxyl radicals in Advanced Oxidation Processes used for wastewater treatment. The results obtained by applying a DFT-based model showed good agreement with experimental values available, as well as qualitative insights into the mechanism of the overall reaction chain. Solvation models were also tried, but were found to be limited for this reaction system within the considered theoretical level without further parameterization.

In-plane mapping of buried InGaAs quantum rings and hybridization effects on the electronic structure

This work reports the investigation on the structural differences between InAs quantum rings and their precursor quantum dots species as well as on the presence of piezoelectric fields and asymmetries in these nanostructures. The experimental results show significant reduction in the ring dimensions when the sizes of capped and uncapped ring and dot samples are compared. The iso-lattice parameter mapped by grazing-incidence x-ray diffraction has revealed the lateral extent of strained regions in the buried rings. A comparison between strain and composition of dot and ring structures allows inferring on how the ring formation and its final configuration may affect optical response parameters. Based on the experimental observations, a discussion has been introduced on the effective potential profile to emulate theoretically the ring-shape confinement. The effects of confinement and strain field modulation on electron and hole band structures are simulated by a multiband k.p calculation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4733964]

The influence of different indium-composition profiles on the electronic structure of lens-shaped InxGa1-xAs quantum dots

We present effective-mass calculations of the bound-state energy levels of electrons confined inside lens-shaped InxGa1-xAs quantum dots (QDs) embedded in a GaAs matrix, taking into account the strain as well as the In gradient inside the QDs due to the strong In segregation and In-Ga intermixing present in the InxGa1-xAs/GaAs system. In order to perform the calculations, we used a continuum model for the strain, and the QDs and wetting layer were divided into their constituting monolayers, each one with a different In concentration, to be able to produce a specific composition profile. Our results clearly show that the introduction of such effects is very important if one desires to correctly reproduce or predict the optoelectronic properties of these nanostructures.

Local order and electronic structure of Pb1-xLaxZr0.40Ti0.60O3 materials and its relation with ferroelectric properties

The knowledge of electronic and local structures is a fundamental step towards understanding the properties of ferroelectric ceramics. X-ray absorption near-edge structure (XANES) of Pb1-xLaxZr0.40Ti0.60O3 ferroelectric samples was measured in order to know how the local order and electronic structure are related to their ferroelectric property, which was tailored by the substitution of lead by lanthanum atoms. The analysis of XANES spectra collected at Ti K- and L-edges XANES showed that the substitution of Pb by La leads to a decrement of local distortion around Ti atoms on the TiO6 octahedron. The analysis of O K-edge XANES spectra showed that the hybridization between O 2p and Pb 6sp states is related to the displacement of Ti atoms in the TiO6 octahedra. Based on these results, it is possible to determine that the degree of ferroelectricity in these samples and the manifestation of relaxor behavior are directly related to the weakening of O 2p and Pb 6sp hybridization. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4720472]

Estimating reaction constants by ab initio molecular modeling: a study on the oxidation of phenol to catechol and hydroquinone in advanced oxidation processes

Molecular modeling is growing as a research tool in Chemical Engineering studies, as can be seen by a simple research on the latest publications in the field. Molecular investigations retrieve information on properties often accessible only by expensive and time-consuming experimental techniques, such as those involved in the study of radical-based chain reactions. In this work, different quantum chemical techniques were used to study phenol oxidation by hydroxyl radicals in Advanced Oxidation Processes used for wastewater treatment. The results obtained by applying a DFT-based model showed good agreement with experimental values available, as well as qualitative insights into the mechanism of the overall reaction chain. Solvation models were also tried, but were found to be limited for this reaction system within the considered theoretical level without further parameterization.

Electronic Structure Calculations in a 2D SixGe1-x Alloy Under an Applied Electric Field.

The recent advances and promises in nanoscience and nanotechnology have been focused on hexagonal materials, mainly on carbon-based nanostructures. Recently, new candidates have been raised, where the greatest efforts are devoted to a new hexagonal and buckled material made of silicon, named Silicene. This new material presents an energy gap due to spin-orbit interaction of approximately 1.5 meV, where the measurement of quantum spin Hall effect(QSHE) can be made experimentally. Some investigations also show that the QSHE in 2D low-buckled hexagonal structures of germanium is present. Since the similarities, and at the same time the differences, between Si and Ge, over the years, have motivated a lot of investigations in these materials. In this work we performed systematic investigations on the electronic structure and band topology in both ordered and disordered SixGe1-x alloys monolayer with 2D honeycomb geometry by first-principles calculations. We show that an applied electric field can tune the gap size for both alloys. However, as a function of electric field, the disordered alloy presents a W-shaped behavior, similarly to the pure Si or Ge, whereas for the ordered alloy a V-shaped behavior is observed.

Ab initio quantum mechanical investigation of structural and chemical-physical properties of selected minerals for minero-petrological, structural ceramic and biomaterial applications

The purpose of this thesis is the atomic-scale simulation of the crystal-chemical and physical (phonon, energetic) properties of some strategically important minerals for structural ceramics, biomedical and petrological applications. These properties affect the thermodynamic stability and rule the mineral-environment interface phenomena, with important economical, (bio)technological, petrological and environmental implications. The minerals of interest belong to the family of phyllosilicates (talc, pyrophyllite and muscovite) and apatite (OHAp), chosen for their importance in industrial and biomedical applications (structural ceramics) and petrophysics. In this thesis work we have applicated quantum mechanics methods, formulas and knowledge to the resolution of mineralogical problems ("Quantum Mineralogy”). The chosen theoretical approach is the Density Functional Theory (DFT), along with periodic boundary conditions to limit the portion of the mineral in analysis to the crystallographic cell and the hybrid functional B3LYP. The crystalline orbitals were simulated by linear combination of Gaussian functions (GTO). The dispersive forces, which are important for the structural determination of phyllosilicates and not properly con-sidered in pure DFT method, have been included by means of a semi-empirical correction. The phonon and the mechanical properties were also calculated. The equation of state, both in athermal conditions and in a wide temperature range, has been obtained by means of variations in the volume of the cell and quasi-harmonic approximation. Some thermo-chemical properties of the minerals (isochoric and isobaric thermal capacity) were calculated, because of their considerable applicative importance. For the first time three-dimensional charts related to these properties at different pressures and temperatures were provided. The hydroxylapatite has been studied from the standpoint of structural and phonon properties for its biotechnological role. In fact, biological apatite represents the inorganic phase of vertebrate hard tissues. Numerous carbonated (hydroxyl)apatite structures were modelled by QM to cover the broadest spectrum of possible biological structural variations to fulfil bioceramics applications.

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.