BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stefan Goedecker

Stefan Goedecker appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.comp-phphysics.atm-clusphysics.chem-phcond-mat.mes-hallcond-mat.othercond-mat.stat-mechcond-mat.supr-concond-matquant-ph

Catalog footprint

What is connected

42works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2026arXiv

The transformation mechanisms among cuboctahedra, Ino's decahedra and icosahedra structures of magic-size gold nanoclusters

Gold nanoclusters possess multiple competing structural motifs with small energy differences, enabling structural coexistence and interconversion. Using a high-accuracy machine learned potential trained on some 20'000 density functional theory reference data points, we investigate transformation pathways connecting both high-symmetry and amorphous cuboctahedra, Ino's decahedra and icosahedra for Au55, Au147, Au309 and Au561 nanoclusters. Our saddle point searches reveal that high-symmetry transformations from cuboctahedra and Ino's decahedra to icosahedra proceed through a single barrier and represent soft-mode-driven jitterbug-type and slip-dislocation motions. In addition, we identify lower-barrier asymmetric transformation pathways that drive the system into disordered, Jahn-Teller-stabilized amorphous icosahedra. Minima Hopping sampling further uncovers, in this context, many such low-symmetry minima. Some of the newly identified global minima for Au309 and Au561 have energies that are up to 2.8 eV lower than the previously reported global minima. Hence, both the shapes and the transformation pathways studied in previous investigations are not the physically relevant ones. In contrast to the previously studied pathways, our transformation pathways give reasonable transformation times that are in rough agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Targeting high symmetry in structure predictions by biasing the potential energy surface

Ground state structures found in nature are in many cases of high symmetry. But structure prediction methods typically render only a small fraction of high symmetry structures. Especially for large crystalline unit cells there are many low energy defect structures. For this reason methods have been developed where either preferentially high symmetry structures are used as input or where the whole structural search is done within a certain symmetry group. In both cases it is necessary to specify the correct symmetry group beforehand. However it can in general not be predicted which symmetry group is the correct one leading to the ground state. For this reason we introduce a potential energy biasing scheme that favors symmetry and where it is not necessary to specify any symmetry group beforehand. On this biased potential energy surface, high symmetry structures will be found much faster than on an unbiased surface and independently of the symmetry group to which they belong. For our two test cases, a $C_{60}$ fullerene and bulk silicon carbide, we get a speedups of 25 and 63. In our data we also find a clear correlation between the similarity of the atomic environments and the energy. In low energy structures all the atoms of a species tend to have similar environments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Manifolds of quasi-constant SOAP and ACSF fingerprints and the resulting failure to machine learn four-body interactions

Atomic fingerprints are commonly used for the characterization of local environments of atoms in machine learning and other contexts. In this work, we study the behavior of two widely used fingerprints, namely the smooth overlap of atomic positions (SOAP) and the atom-centered symmetry functions (ACSF), under finite changes of atomic positions and demonstrate the existence of manifolds of quasi-constant fingerprints. These manifolds are found numerically by following eigenvectors of the sensitivity matrix with quasi-zero eigenvalues. The existence of such manifolds in ACSF and SOAP causes a failure to machine learn four-body interactions such as torsional energies that are part of standard force fields. No such manifolds can be found for the Overlap Matrix (OM) fingerprint due to its intrinsic many-body character.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Potential energy surface study of X@Si$_{32}$X$^-_{44}$(X=Cl, Br) clusters to decipher the stabilization process of Si$_{20}$ fullerene

Efforts toward stabilization of the Si$_{20}$ fullerene through different schemes have failed despite several theoretical predictions. However, recently Tillmann {\it et. al.} succeeded to stabilize the Si$_{20}$ fullerene through exohedral decoration with eight Cl substituents and twelve SiCl$_3$ groups on the surface and enclosing Cl$^-$ ion. A deeper understanding on what factors lead to stabilization will open the path for stabilizing other systems of interest. Here, we employ the minima hopping method within density functional theory to understand the potential energy surface. The study shows that the exo-endo halide decoration of the cage alters the glassy nature of the potential energy surface of pure cage to structure seeker. Further analysis of different properties of the global minima, reveal that the extra electron instead of residing on the central encapsulated atom in the cage, it is distributed on the cage and increases the encapsulation energy; thereby stabilizing the system. We also provide estimates of the stability for different kind of exo-endo halide decorations and their feasible realization in experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An assessment of the structural resolution of various fingerprints commonly used in machine learning

Atomic environment fingerprints are widely used in computational materials science, from machine learning potentials to the quantification of similarities between atomic configurations. Many approaches to the construction of such fingerprints, also called structural descriptors, have been proposed. In this work, we compare the performance of fingerprints based on the Overlap Matrix(OM), the Smooth Overlap of Atomic Positions (SOAP), Behler-Parrinello atom-centered symmetry functions (ACSF), modified Behler-Parrinello symmetry functions (MBSF) used in the ANI-1ccx potential and the Faber-Christensen-Huang-Lilienfeld (FCHL) fingerprint under various aspects. We study their ability to resolve differences in local environments and in particular examine whether there are certain atomic movements that leave the fingerprints exactly or nearly invariant. For this purpose, we introduce a sensitivity matrix whose eigenvalues quantify the effect of atomic displacement modes on the fingerprint. Further, we check whether these displacements correlate with the variation of localized physical quantities such as forces. Finally, we extend our examination to the correlation between molecular fingerprints obtained from the atomic fingerprints and global quantities of entire molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detecting non-local effects in the electronic structure of a simple covalent system with machine learning methods

Using methods borrowed from machine learning we detect in a fully algorithmic way long range effects on local physical properties in a simple covalent system of carbon atoms. The fact that these long range effects exist for many configurations implies that atomistic simulation methods, such as force fields or modern machine learning schemes, that are based on locality assumptions, are limited in accuracy. We show that the basic driving mechanism for the long range effects is charge transfer. If the charge transfer is known, locality can be recovered for certain quantities such as the band structure energy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Funnel Hopping Monte Carlo: An efficient method to overcome broken ergodicity

Monte Carlo simulations are a powerful tool to investigate the thermodynamic properties of atomic systems. In practice however, sampling of the complete configuration space is often hindered by high energy barriers between different regions of configuration space which can make ergodic sampling completely infeasible within accessible simulation times. Although several extensions to the conventional Monte Carlo scheme have been developed, that enable the treatment of such systems, these extensions often entail substantial computational cost or rely on the harmonic approximation. In this work we propose an exact method called Funnel Hopping Monte Carlo (FHMC) that is inspired by the the ideas of smart darting but is more efficient. Gaussian mixtures are used to approximate the Boltzmann distribution around local energy minima which are then used to propose high quality Monte Carlo moves that enable the Monte Carlo simulation to directly jump between different funnels. To fit the Gaussian mixtures we developed an extended version of the expectation-maximization algorithm that is able to take advantage of the high symmetry present in many low energy configurations. We demonstrate the methods performance on the example of the 38 as well as the 75 atom Lennard-Jones clusters which are well known for their double funnel energy landscapes that prevent ergodic sampling with conventional Monte Carlo simulations. By integrating FHMC into the parallel tempering scheme we were able to reduce the number of steps required until convergence of the simulation significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximum volume simplex method for automatic selection and classification of atomic environments and environment descriptor compression

Fingerprint distances, which measure the similarity of atomic environments, are commonly calculated from atomic environment fingerprint vectors. In this work we present the simplex method which can perform the inverse operation, i.e. calculating fingerprint vectors from fingerprint distances. The fingerprint vectors found in this way point to the corners of a simplex. For a large data set of fingerprints, we can find a particular largest volume simplex, whose dimension gives the effective dimension of the fingerprint vector space. We show that the corners of this simplex correspond to landmark environments that can by used in a fully automatic way to analyse structures. In this way we can for instance detect atoms in grain boundaries or on edges of carbon flakes without any human input about the expected environment. By projecting fingerprints on the largest volume simplex we can also obtain fingerprint vectors that are considerably shorter than the original ones but whose information content is not significantly reduced.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computationally Efficient Characterization of Potential Energy Surfaces Based on Fingerprint Distances

An analysis of the network defined by the potential energy minima of multi-atomic systems and their connectivity via reaction pathways that go through transition states allows to understand important characteristics like thermodynamic, dynamic and structural properties. Unfortunately computing the transition states and reaction pathways in addition to the significant energetically low-lying local minima is a computationally demanding task. We here introduce a computationally efficient method that is based on a combination of the minima hopping global optimization method and the insight that uphill barriers tend to increase with increasing structural distances of the educt and product states. This method allows to replace the exact connectivity information and transition state energies with alternative and approximate concepts. Without adding any significant additional cost to the minima hopping global optimization approach, this method allows to generate an approximate network of the minima, their connectivity and a rough measure for the energy needed for their interconversion. This can be used to obtain a first qualitative idea on important physical and chemical properties by means of a disconnectivity graph analysis. Besides the physical insight obtained by such an analysis, the gained knowledge can be used to make a decision if it is worthwhile or not to invest computational resources for an exact computation of the transition states and the reaction pathways. Furthermore it is demonstrated that the here presented method can be used for finding physically reasonable interconversion pathways that are promising input pathways for methods like transition path sampling or discrete path sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Emergence of superconductivity in doped H$_2$O ice at high pressure

We investigate the possibility of achieving high-temperature superconductivity in hydrides under pressure by inducing metallization of otherwise insulating phases through doping, a path previously used to render standard semiconductors superconducting at ambient pressure. Following this idea, we study H$_2$O, one of the most abundant and well-studied substances, we identify nitrogen as the most likely and promising substitution/dopant. We show that for realistic levels of doping of a few percent, the phase X of ice becomes superconducting with a critical temperature of about 60 K at 150GPa. In view of the vast number of hydrides that are strongly covalent bonded, but that remain insulating until rather large pressures, our results open a series of new possibilities in the quest for novel high-temperature superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ultralow Thermal Conductivity in Full-Heusler Semiconductors

Semiconducting half- and, to a lesser extent, full-Heusler compounds are promising thermoelectric materials due to their compelling electronic properties with large power factors. However, intrinsically high thermal conductivity resulting in a limited thermoelectric efficiency has so far impeded their widespread use in practical applications. Here, we report the computational discovery of a class of hitherto unknown stable semiconducting full-Heusler compounds with ten valence electrons ($X_2YZ$, $X$=Ca, Sr, and Ba; $Y$= Au and Hg; $Z$=Sn, Pb, As, Sb, and Bi) through high-throughput $ab-initio$ screening. These new compounds exhibit ultralow lattice thermal conductivity $κ_{\text{L}}$ close to the theoretical minimum due to strong anharmonic rattling of the heavy noble metals, while preserving high power factors, thus resulting in excellent phonon-glass electron-crystal materials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A fingerprint based metric for measuring similarities of crystalline structures

Measuring similarities/dissimilarities between atomic structures is important for the exploration of potential energy landscapes. However, the cell vectors together with the coordinates of the atoms, which are generally used to describe periodic systems, are quantities not suitable as fingerprints to distinguish structures. Based on a characterization of the local environment of all atoms in a cell we introduce crystal fingerprints that can be calculated easily and allow to define configurational distances between crystalline structures that satisfy the mathematical properties of a metric. This distance between two configurations is a measure of their similarity/dissimilarity and it allows in particular to distinguish structures. The new method is an useful tool within various energy landscape exploration schemes, such as minima hopping, random search, swarm intelligence algorithms and high-throughput screenings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate and efficient linear scaling DFT calculations with universal applicability

Density Functional Theory calculations traditionally suffer from an inherent cubic scaling with respect to the size of the system, making big calculations extremely expensive. This cubic scaling can be avoided by the use of so-called linear scaling algorithms, which have been developed during the last few decades. In this way it becomes possible to perform ab-initio calculations for several tens of thousands of atoms or even more within a reasonable time frame. However, even though the use of linear scaling algorithms is physically well justified, their implementation often introduces some small errors. Consequently most implementations offering such a linear complexity either yield only a limited accuracy or, if one wants to go beyond this restriction, require a tedious fine tuning of many parameters. In our linear scaling approach within the BigDFT package, we were able to overcome this restriction. Using an ansatz based on localized support functions expressed in an underlying Daubechies wavelet basis -- which offers ideal properties for accurate linear scaling calculations -- we obtain an amazingly high accuracy and a universal applicability while still keeping the possibility of simulating large systems with only a moderate demand of computing resources.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Identification of novel Cu, Ag, and Au ternary oxides from global structural prediction

We use ab initio global structural prediction, and specifically the minima hopping method, to explore the periodic table in search of novel oxide phases. In total, we study 183 different compositions of the form MXO2, where M=(Cu, Ag, Au) and X is an element of the periodic table. This set includes the well-known Cu delafossite compounds that are, up to now, the best p-type transparent conductive oxides known to mankind. Our calculations discover 81 stable compositions, out of which only 36 are included in available databases. Some of these new phases are potentially good candidates for transparent electrodes. These results demonstrate, on one hand, how incomplete is still our knowledge of the phase-space of stable ternary materials. On the other hand, we show that structural prediction combined with high-throughput approaches is a powerful tool to extend that knowledge, paving the way for the experimental discovery of new materials on a large scale.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interatomic potentials for ionic systems with density functional accuracy based on charge densities obtained by a neural network

Based on an analysis of the short range chemical environment of each atom in a system, standard machine learning based approaches to the construction of interatomic potentials aim at determining directly the central quantity which is the total energy. This prevents for instance an accurate description of the energetics of systems where long range charge transfer is important as well as of ionized systems. We propose therefore not to target directly with machine learning methods the total energy but an intermediate physical quantity namely the charge density, which then in turn allows to determine the total energy. By allowing the electronic charge to distribute itself in an optimal way over the system, we can describe not only neutral but also ionized systems with unprecedented accuracy. We demonstrate the power of our approach for both neutral and ionized NaCl clusters where charge redistribution plays a decisive role for the energetics. We are able to obtain chemical accuracy, i.e. errors of less than a milli Hartree per atom compared to the reference density functional results. The introduction of physically motivated quantities which are determined by the short range atomic environment via a neural network leads also to an increased stability of the machine learning process and transferability of the potential.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low-density silicon allotropes for photovoltaic applications

Silicon materials play a key role in many technologically relevant fields, ranging from the electronic to the photovoltaic industry. A systematic search for silicon allotropes was performed by employing a modified ab initio minima hopping crystal structure prediction method. The algorithm was optimized to specifically investigate the hitherto barely explored low-density regime of the silicon phase diagram by imitating the guest-host concept of clathrate compounds. In total 44 metastable phases are presented, of which 11 exhibit direct or quasi-direct band-gaps in the range of $\approx$1.0-1.8 eV, close to the optimal Shockley-Queisser limit of $\approx$1.4 eV, with a stronger overlap of the absorption spectra with the solar spectrum compared to conventional diamond silicon. Due to the structural resemblance to known clathrate compounds it is expected that the predicted phases can be synthesized.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Superconductivity in metastable phases of phosphorus-hydride compounds under high pressure

Hydrogen-rich compounds have been extensively studied both theoretically and experimentally in the quest for novel high-temperature superconductors. Reports on sulfur-hydride attaining metallicity under pressure and exhibiting superconductivity at temperatures as high as 200 K have spurred an intense search for room-temperature superconductors in hydride materials. Recently, compressed phosphine was reported to metallize at pressures above 45 GPa, reaching a superconducting transition temperature (T$_{c}$) of 100 K at 200 GPa. However, neither the exact composition nor the crystal structure of the superconducting phase have been conclusively determined. In this work the phase diagram of PH$_n$ ($n=1,2,3,4,5,6$) was extensively explored by means of {\it ab initio} crystal structure predictions using the Minima Hopping Method (MHM). The results do not support the existence of thermodynamically stable PH$_n$ compounds, which exhibit a tendency for elemental decomposition at high pressure even when vibrational contributions to the free energies are taken into account. Although the lowest energy phases of PH$_{1,2,3}$ display T$_{c}$'s comparable to experiments, it remains uncertain if the measured values of T$_{c}$ can be fully attributed to a phase-pure compound of PH$_n$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Daubechies Wavelets for Linear Scaling Density Functional Theory

We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First-principles predicted low-energy structures of NaSc(BH4)4

According to previous interpretations of experimental data, sodium-scandium double-cation borohydride NaSc(BH$_4$)$_4$ crystallizes in the crystallographic space group $Cmcm$ where each sodium (scandium) atom is surrounded by six scandium (sodium) atoms. A careful investigation of this phase based on \textit{ab initio} calculations indicates that the structure is dynamically unstable and gives rise to an energetically and dynamically more favorable phase with $C222_1$ symmetry and nearly identical x-ray diffraction pattern. By additionally performing extensive structural searches with the minima-hopping method we discover a class of new low-energy structures exhibiting a novel structural motif in which each sodium (scandium) atom is surrounded by four scandium (sodium) atoms arranged at the corners of either a rectangle with nearly equal sides or a tetrahedron. These new phases are all predicted to be insulators with band gaps of $7.9-8.2$ eV. Finally, we estimate the influence of these structures on the hydrogen-storage performance of NaSc(BH$_4$)$_4$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Minima Hopping Guided Path Search: An Efficient Method for Finding Complex Chemical Reaction Pathways

The Minima Hopping global optimization method uses physically realizable molecular dynamics moves in combination with an energy feedback that guarantees the escape from any potential energy funnel. For the purpose of finding reactions pathways, we argue that Minima Hopping is particularly suitable as a guide through the potential energy landscape and as a generator for pairs of minima that can be used as input structures for methods capable of finding transition states between two minima. For Lennard-Jones benchmark systems we compared this Minima Hopping guided path search method to a known approach for the exploration of potential energy landscapes that is based on deterministic mode-following. Although we used a stabilized mode-following technique that reliably allows to follow distinct directions when escaping from a local minimum, we observed that Minima Hopping guided path search is far superior in finding lowest-barrier reaction pathways. We therefore suggest that Minima Hopping guided path search can be used as a simple and efficient way to identify energetically low-lying chemical reaction pathways. Finally we applied the Minima Hopping guided path search approach to 75-atom and 102-atom Lennard Jones systems. For the 75-atom system we found pathways whose highest energies are significantly lower than the highest energy along the previously published lowest-barrier pathway. Furthermore, many of these pathways contain a smaller number of intermediate transition states than the previously publish lowest-barrier pathway. In case of the 102-atom system Minima Hopping guided path search found a previously unknown and energetically low-lying funnel.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stabilized Quasi-Newton Optimization of Noisy Potential Energy Surfaces

Optimizations of atomic positions belong to the most commonly performed tasks in electronic structure calculations. Many simulations like global minimum searches or characterizations of chemical reactions require performing hundreds or thousands of minimizations or saddle computations. To automatize these tasks, optimization algorithms must not only be efficient, but also very reliable. Unfortunately computational noise in forces and energies is inherent to electronic structure codes. This computational noise poses a sever problem to the stability of efficient optimization methods like the limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm. We here present a technique that allows obtaining significant curvature information of noisy potential energy surfaces. We use this technique to construct both, a stabilized quasi-Newton minimization method and a stabilized quasi-Newton saddle finding approach. We demonstrate with the help of benchmarks that both the minimizer and the saddle finding approach are superior to comparable existing methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Customized 3D GPU Poisson Solver for Free Boundary Conditions

A 3-dimensional GPU Poisson solver is developed for all possible combinations of free and periodic boundary conditions (BCs) along the three directions. It is benchmarked for various grid sizes and different BCs and a significant performance gain is observed for problems including one or more free BCs. The GPU Poisson solver is also benchmarked against two different CPU implementations of the same method and a significant amount of acceleration of the computation is observed with the GPU version.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Boron aggregation in the ground states of boron-carbon fullerenes

We present novel structural motifs for boron-carbon nano-cages of the stochiometries B12C48 and B12C50, based on first principle calculations. These configurations are distinct from those proposed so far by the fact that the boron atoms are not isolated and distributed over the entire surface of the cages, but rather aggregate at one location to form a patch. Our putative ground state of B12C48 is 1.8 eV lower in energy than the previously proposed ground state and violates all the suggested empirical rules for constructing low energy fullerenes. The B12C50 configuration is energetically even more favorable than B12C48, showing that structures derived from the C60 Buckminsterfullerene are not necessarily magic sizes for heterofullerenes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conducting boron sheets formed by the reconstruction of the α-boron (111) surface

Systematic ab initio structure prediction was applied for the first time to predict low energy surface reconstructions by employing the minima hopping method on the α-boron (111) surface. Novel reconstruction geometries were identified and carefully characterized in terms of structural and electronic properties. Our calculations predict the formation of a planar, mono-layer sheet at the surface, which is responsible for conductive surface states. Furthermore, the isolated boron sheet is shown to be the ground state 2D-structure in vacuum at a hole density of η=1/5 and is therefore a potential candidate as a precursor for boron nano-structures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Glassy clusters: Relations between their dynamics and characteristic features of their energy landscape

Based on a recently introduced metric for measuring distances between configurations, we in- troduce distance-energy (DE) plots to characterize the potential energy surface (PES) of clusters. Producing such plots is computationally feasible on the density functional (DFT) level since it re- quires only a set of a few hundred stable low energy configurations including the global minimum. By comparison with standard criteria based on disconnectivity graphs and on the dynamics of Lennard- Jones clusters we show that the DE plots convey the necessary information about the character of the potential energy surface and allow to distinguish between glassy and non-glassy systems. We then apply this analysis to real systems on the DFT level and show that both glassy and non-glassy clusters can be found in simulations. It however turns out that among our investigated clusters only those can be synthesized experimentally which exhibit a non-glassy landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Low-Energy Polymeric Phases of Alanates

Low-energy structures of alanates are currently known to be described by patterns of isolated, nearly ideal tetrahedral [AlH$_4$] anions and metal cations. We discover that the novel polymeric motif recently proposed for LiAlH$_4$ plays a dominant role in a series of alanates, including LiAlH$_4$, NaAlH$_4$, KAlH$_4$, Mg(AlH$_4$)$_2$, Ca(AlH$_4$)$_2$ and Sr(AlH$_4$)$_2$. In particular, most of the low-energy structures discovered for the whole series are characterized by networks of corner-sharing [AlH$_6$] octahedra, forming wires and/or planes throughout the materials. Finally, for Mg(AlH$_4$)$_2$ and Sr(AlH$_4$)$_2$, we identify two polymeric phases to be lowest in energy at low temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Metrics for measuring distances in configuration spaces

In order to characterize molecular structures we introduce configurational fingerprint vectors which are counterparts of quantities used experimentally to identify structures. The Euclidean distance between the configurational fingerprint vectors satisfies the properties of a metric and can therefore safely be used to measure dissimilarities between configurations in the high dimensional configuration space. We show that these metrics correlate well with the RMSD between two configurations if this RMSD is obtained from a global minimization over all translations, rotations and permutations of atomic indices. We introduce a Monte Carlo approach to obtain this global minimum of the RMSD between configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Norm-conserving pseudopotentials with chemical accuracy compared to all-electron calculations

By adding a non-linear core correction to the well established Dual Space Gaussian type pseudopotentials for the chemical elements up to the third period, we construct improved pseudopotentials for the Perdew Burke Ernzerhof (PBE) functional and demonstrate that they exhibit excellent accuracy. Our benchmarks for the G2-1 test set show average atomization energy errors of only half a kcal/mol. The pseudopotentials also remain highly reliable for high pressure phases of crystalline solids. When supplemented by empirical dispersion corrections the average error in the interaction energy between molecules is also about half a kcal/mol. The accuracy that can be obtained by these pseudopotentials in combination with a systematic basis set is well superior to the accuracy that can be obtained by commonly used medium size Gaussian basis sets in all-electron calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermodynamic stability of alkali metal/zinc double-cation borohydrides at low temperatures

We study the thermodynamic stability at low temperatures of a series of alkali metal/zinc double-cation borohydrides, including LiZn(BH$_4$)$_3$, LiZn$_2$(BH$_4$)$_5$, NaZn(BH$_4$)$_3$, NaZn$_2$(BH$_4$)$_5$, KZn(BH$_4$)$_3$, and KZn$_2$(BH$_4$)$_5$. While LiZn$_2$(BH$_4$)$_5$, NaZn(BH$_4$)$_3$, NaZn$_2$(BH$_4$)$_5$ and KZn(BH$_4$)$_3$ were recently synthesized, LiZn(BH$_4$)$_3$ and KZn$_2$(BH$_4$)$_5$ are hypothetical compounds. Using the minima-hopping method, we discover two new lowest-energy structures for NaZn(BH$_4$)$_3$ and KZn$_2$(BH$_4$)$_5$ which belong to the $C2/c$ and $P2$ space groups, respectively. These structures are predicted to be both thermodynamically stable and dynamically stable, implying that their existence may be possible. On the other hand, the lowest-energy $P1$ structure of LiZn(BH$_4$)$_3$ is predicted to be unstable, suggesting a possible reason elucidating why this compound has not been experimentally identified. In exploring the low-energy structures of these compounds, we find that their energetic ordering is sensitive to the inclusion of the van der Waals interactions. We also find that a proper treatment of these interactions, e.g., as given by a non-local density functional such as vdW-DF2, is necessary to address the stability of the low-energy structures of these compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Multi-scale Approach for Simulations of Kelvin Probe Force Microscopy with Atomic Resolution

The distance dependence and atomic-scale contrast observed in nominal contact potential difference (CPD) signals recorded by KPFM on surfaces of insulating and semiconducting samples, have stimulated theoretical attempts to explain such effects. We attack this problem in two steps. First, the electrostatics of the macroscopic tip-cantilever-sample system is treated by a finite-difference method on an adjustable nonuniform mesh. Then the resulting electric field under the tip apex is inserted into a series of atomistic wavelet-based density functional theory (DFT) calculations. Results are shown for a realistic neutral but reactive silicon nano-scale tip interacting with a NaCl(001) sample. Bias-dependent forces and resulting atomic displacements are computed to within an unprecedented accuracy. Theoretical expressions for amplitude modulation (AM) and frequency modulation (FM) KPFM signals and for the corresponding local contact potential differences (LCPD) are obtained by combining the macroscopic and atomistic contributions to the electrostatic force component generated at the voltage modulation frequency, and evaluated for several tip oscillation amplitudes A up to 10 nm. Being essentially constant over a few Volts, the slope of atomistic force versus bias is the basic quantity which determines variations of the atomic-scale LCPD contrast. Already above A = 0.1 nm, the LCPD contrasts in both modes exhibit almost the same spatial dependence as the slope. In the AM mode, this contrast is approximately proportional to $A^{-1/2}$, but remains much weaker than the contrast in the FM mode, which drops somewhat faster as A is increased. These trends are a consequence of the macroscopic contributions to the KPFM signal, which are stronger in the AM-mode and especially important if the sample is an insulator even at sub-nanometer separations where atomic-scale contrast appears.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Low-energy silicon allotropes with strong absorption in the visible for photovoltaic applications

We present state-of-the-art first-principle calculations of the electronic and optical properties of silicon allotropes with interesting characteristics for applications in thin-film solar cells. These new phases consist of distorted sp$^3$ silicon networks and have a lower formation energy than other experimentally produced silicon phases. Some of these structures turned out to have quasi-direct and dipole-allowed band gaps in the range 0.8--1.5\,eV, and to display absorption coefficients comparable with those of chalcopyrites used in thin-film record solar cells.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Low-energy structures of zinc borohydride Zn(BH$_4$)$_2$

We present a systematic study of the low-energy structures of zinc borohydride, a crystalline material proposed for the hydrogen storage purpose. In addition to the previously proposed structures, many new low-energy structures of zinc borohydride are found by utilizing the minima-hopping method. We identify a new dynamically stable structure which belongs to the $I4_122$ space group as the most stable phase of zinc borohydride at low temperatures. A low transition barrier between $I4_122$ and $P1$, the two lowest-lying phases of zinc borohydride is predicted, implying that a coexistence of low-lying phases of zinc borohydride is possible at ambient conditions. An analysis based on the simulated X-ray diffraction pattern reveals that the $I4_122$ structure exhibits the same major features as the experimentally synthesized zinc borohydride samples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Prediction of a novel monoclinic carbon allotrope

A novel allotrope of carbon with $P2/m$ symmetry was identified during an \emph{ab-initio} minima-hopping structural search which we call $M10$-carbon. This structure is predicted to be more stable than graphite at pressures above 14.4 GPa and consists purely of $sp^3$ bonds. It has a high bulk modulus and is almost as hard as diamond. A comparison of the simulated X-ray diffraction pattern shows a good agreement with experimental results from cold compressed graphite.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The crystal structure of cold compressed graphite

Through a systematic structural search we found an allotrope of carbon with Cmmm symmetry which we predict to be more stable than graphite for pressures above 10 GPa. This material, which we refer to as Z-carbon, is formed by pure sp3 bonds and is the only carbon allotrope which provides an excellent match to unexplained features in experimental X-ray diffraction and Raman spectra of graphite under pressure. The transition from graphite to Z-carbon can occur through simple sliding and buckling of graphene sheets. Our calculations predict that Z-carbon is a transparent wide band gap semiconductor with a hardness comparable to diamond.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effect of ionization on the global minima of small and medium sized silicon and magnesium clusters

We re-examine the question of whether the geometrical ground state of neutral and ionized clusters are identical. Using a well defined criterion for being "identical" together, the extensive sampling methods on a potential energy surface calculated by density functional theory, we show that the ground states are in general different. This behavior is to be expected whenever there are metastable configurations which are close in energy to the ground state, but it disagrees with previous studies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The energy landscape of fullerene materials: a comparison between boron, boron-nitride and carbon

Using the minima hopping global geometry optimization method on the density functional potential energy surface we study medium size and large boron clusters. Even though for isolated medium size clusters the ground state is a cage like structure they are unstable against external perturbations such as contact with other clusters. The energy landscape of larger boron clusters is glass like and has a large number of structures which are lower in energy than the cages. This is in contrast to carbon and boron nitride systems which can be clearly identified as structure seekers in our minima hopping runs. The differences in the potential energy landscape explain why carbon and boron nitride systems are found in nature whereas pure boron fullerenes have not been found.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystal structure prediction using the Minima Hopping method

A structure prediction method is presented based on the Minima Hopping method. Optimized moves on the configurational enthalpy surface are performed to escape local minima using variable cell shape molecular dynamics by aligning the initial atomic and cell velocities to low curvature directions of the current minimum. The method is applied to both silicon crystals and binary Lennard-Jones mixtures and the results are compared to previous investigations. It is shown that a high success rate is achieved and a reliable prediction of unknown ground state structures is possible.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient moves for global geometry optimization methods and their application to binary systems

We show that molecular dynamics based moves in the Minima Hopping (MH) method are more efficient than saddle point crossing moves which select the lowest possible saddle point. For binary systems we incorporate identity exchange moves in a way that allows to avoid the generation of high energy configurations. Using this modified Minima Hopping method we reexamine the binary Lennard Jones (BLJ) benchmark system with up to 100 atoms and we find a large number of new putative global minima structures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structural Metastability of Endohedral Silicon Fullerenes

Endohedrally doped Si20 fullerenes appear as appealing building blocks for nanoscale materials. We investigate their structural stability with an unbiased and systematic global geometry optimization method within density-functional theory. For a wide range of metal doping atoms, it was sufficient to explore the Born Oppenheimer surface for only a moderate number of local minima to find structures that clearly differ from the initial endohedral cages, but are considerably more favorable in terms of energy. Previously proposed structures are thus all metastable.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Daubechies wavelets as a basis set for density functional pseudopotential calculations

Daubechies wavelets are a powerful systematic basis set for electronic structure calculations because they are orthogonal and localized both in real and Fourier space. We describe in detail how this basis set can be used to obtain a highly efficient and accurate method for density functional electronic structure calculations. An implementation of this method is available in the ABINIT free software package. This code shows high systematic convergence properties, very good performances and an excellent efficiency for parallel calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Optimization and Parallelization of a force field for silicon using OpenMP

The force field by Lenosky and coworkers is the latest force field for silicon which is one of the most studied materials. It has turned out to be highly accurate in a large range of test cases. The optimization and parallelization of this force field using OpenMp and Fortan90 is described here. The optimized program allows us to handle a very large number of silicon atoms in large scale simulations. Since all the parallelization is hidden in a single subroutine that returns the total energies and forces, this subroutine can be called from within a serial program in an user friendly way.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The solution of multi-scale partial differential equations using wavelets

Wavelets are a powerful new mathematical tool which offers the possibility to treat in a natural way quantities characterized by several length scales. In this article we will show how wavelets can be used to solve partial differential equations which exhibit widely varying length scales and which are therefore hardly accessible by other numerical methods. As a benchmark calculation we solve Poisson's equation for a 3-dimensional Uranium dimer. The length scales of the charge distribution vary by 4 orders of magnitude in this case. Using lifted interpolating wavelets the number of iterations is independent of the maximal resolution and the computational effort therefore scales strictly linearly with respect to the size of the system.

0 citations0 reviewsNo code linkedOpen access