Researcher profile

Stefan Goedecker

Stefan Goedecker contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 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.

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.

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.

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.

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.

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.

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.

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.