Source author record

Bastian Schaefer

Bastian Schaefer 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

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

8 published item(s)

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.

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.

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.

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

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.

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.

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.

preprint2010arXiv

Agglomeration and filtration of colloidal suspensions with DVLO interactions in simulation and experiment

Cake filtration is a widely used solid-liquid separation process. However, the high flow resistance of the nanoporous filter cake lowers the efficiency of the process significantly. The structure and thus the permeability of the filter cakes depend on the compressive load acting on the particles, the particles size, and the agglomeration of the particles. The latter is determined by the particle charge and the ionic strength of the suspension, as described by the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. In this paper, we propose a combined stochastic rotation dynamics (SRD) and molecular dynamics (MD) methodology to simulate the cake formation. The simulations give further insight into the dependency of the filter cakes' structure on the agglomeration of the particles, which cannot be accessed experimentally. The permeability, as investigated with lattice Boltzmann (LB) simulations of flow through the discretized cake, depends on the particle size and porosity, and thus on the agglomeration of the particles. Our results agree qualitatively with experimental data obtained from colloidal boehmite suspensions.