Source author record

S. Alireza Ghasemi

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

11 published item(s)

preprint2021arXiv

An automated approach for developing neural network interatomic potentials with FLAME

The performance of machine learning interatomic potentials relies on the quality of the training dataset. In this work, we present an approach for generating diverse and representative training data points which initiates with \it{ab initio} calculations for bulk structures. The data generation and potential construction further proceed side-by-side in a cyclic process of training the neural network and crystal structure prediction based on the developed interatomic potentials. All steps of the data generation and potential development are performed with minimal human intervention. We show the reliability of our approach by assessing the performance of neural network potentials developed for two inorganic systems.

preprint2019arXiv

FLAME: a library of atomistic modeling environments

FLAME is a software package to perform a wide range of atomistic simulations for exploring the potential energy surfaces (PES) of complex condensed matter systems. The range of methods include molecular dynamics simulations to sample free energy landscapes, saddle point searches to identify transition states, and gradient relaxations to find dynamically stable geometries. In addition to such common tasks, FLAME implements a structure prediction algorithm based on the minima hopping method (MHM) to identify the ground state structure of any system given solely the chemical composition, and a framework to train a neural network potential to reproduce the PES from $\textit{ab initio}$ calculations. The combination of neural network potentials with the MHM in FLAME allows a highly efficient and reliable identification of the ground state as well as metastable structures of molecules and crystals, as well as of nano structures, including surfaces, interfaces, and two-dimensional materials. In this manuscript, we provide detailed descriptions of the methods implemented in the FLAME code and its capabilities, together with several illustrative examples.

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.

preprint2015arXiv

Highly accurate and efficient algorithm for electrostatic interaction of charged particles confined by metallic parallel plates

We present an accurate and efficient algorithm to calculate the electrostatic interaction of charged point particles with partially periodic boundary conditions that are confined along the nonperiodic direction by two metallic parallel plates. The method preserves the original boundary conditions and hence, it does not introduce any kind of artifacts. In addition, it enjoys the quasilinear complexity of $\mathcal{O}(N\ln(N))$, where $N$ being the number of particles in the simulation box. In fact, based on the superposition principle in electrostatics, the problem is split into two electrostatic problems where each one can be calculated by the appropriate Poisson solver. In this paper we apply the method to sodium chloride ultrathin films and investigate its dielectric response with respect to external bias voltage. We show how accurately in this method one can obtain the total charge induced on metallic boundaries to an arbitrary precision.

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.

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.

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.

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.

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.

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.

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.