Source author record

Matous Mrovec

Matous Mrovec 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-sci physics.comp-ph

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Effects of thermal, elastic, and surface properties on the stability of SiC polytypes

SiC polytypes have been studied for decades, both experimentally and with atomistic simulations, yet no consensus has been reached on the factors that determine their stability and growth. Proposed governing factors are temperature-dependent differences in the bulk energy, biaxial strain induced through point defects, and surface properties. In this work, we investigate the thermodynamic stability of the 3C, 2H, 4H, and 6H polytypes with density functional theory (DFT) calculations. The small differences of the bulk energies between the polytypes can lead to intricate changes in their energetic ordering depending on the computational method. Therefore, we employ and compare various DFT-codes: VASP, CP2K, and FHI-aims; exchange-correlation functionals: LDA, PBE, PBEsol, PW91, HSE06, SCAN, and RTPSS; and nine different van der Waals (vdW) corrections. At $T=0$~K, 4H-SiC is marginally more stable than 3C-SiC, and the stability further increases with temperature by including entropic effects from lattice vibrations. Neither the most advanced vdW corrections nor strain on the lattice have a significant effect on the relative polytype stability. We further investigate the energies of the (0001) polytype surfaces that are commonly exposed during epitaxial growth. For Si-terminated surfaces, we find 3C-SiC to be significantly more stable than 4H-SiC. We conclude that the difference in surface energy is likely the driving force for 3C-nucleation, whereas the difference in the bulk thermodynamic stability slightly favors the 4H and 6H polytypes. In order to describe the polytype stability during crystal growth correctly, it is thus crucial to take into account both of these effects.

preprint2021arXiv

Performant implementation of the atomic cluster expansion (PACE): Application to copper and silicon

The atomic cluster expansion is a general polynomial expansion of the atomic energy in multi-atom basis functions. Here we implement the atomic cluster expansion in the performant C++ code \verb+PACE+ that is suitable for use in large scale atomistic simulations. We briefly review the atomic cluster expansion and give detailed expressions for energies and forces as well as efficient algorithms for their evaluation. We demonstrate that the atomic cluster expansion as implemented in \verb+PACE+ shifts a previously established Pareto front for machine learning interatomic potentials towards faster and more accurate calculations. Moreover, general purpose parameterizations are presented for copper and silicon and evaluated in detail. We show that the new Cu and Si potentials significantly improve on the best available potentials for highly accurate large-scale atomistic simulations.

preprint2019arXiv

Phase transitions in titanium with an analytic bond-order potential

Titanium is the base material for a number of technologically important alloys for energy conversion and structural applications. Atomic-scale studies of Ti-based metals employing first-principles methods, such as density functional theory, are limited to ensembles of a few hundred atoms. To perform large-scale and/or finite temperature simulations, computationally more efficient interatomic potentials are required. In this work, we coarse grain the tight-binding (TB) approximation to the electronic structure and develop an analytic bond-order potential (BOP) for Ti by fitting to the energies and forces of elementary deformations of simple structures. The BOP predicts the structural properties of the stable and defective phases of Ti with a quality comparable to previous TB parametrizations at a much lower computational cost. The predictive power of the model is demonstrated for simulations of martensitic transformations.

preprint2015arXiv

Atomistically enabled nonsingular anisotropic elastic representation of near-core dislocation stress fields in $α$-iron

The stress fields of dislocations predicted by classical elasticity are known to be unrealistically large approaching the dislocation core, due to the singular nature of the theory. While in many cases this is remedied with the approximation of an effective core radius, inside which ad hoc regularizations are implemented, such approximations lead to a compromise in the accuracy of the calculations. In this work, an anisotropic non-singular elastic representation of dislocation fields is developed to accurately represent the near-core stresses of dislocations in $α$-iron. The regularized stress field is enabled through the use of a non-singular Green's tensor function of Helmholtz-type gradient anisotropic elasticity, which requires only a single characteristic length parameter in addition to the material's elastic constants. Using a novel magnetic bond-order potential to model atomic interactions in iron, molecular statics calculations are performed, and an optimization procedure is developed to extract the required length parameter. Results show the method can accurately replicate the magnitude and decay of the near-core dislocation stresses even for atoms belonging to the core itself. Comparisons with the singular isotropic and anisotropic theories show the non-singular anisotropic theory leads to a substantially more accurate representation of the stresses of both screw and edge dislocations near the core, in some cases showing improvements in accuracy of up to an order of magnitude. The spatial extent of the region in which the singular and non-singular stress differ substantially is also discussed. The general procedure we describe may in principle be applied to accurately model the near-core dislocation stresses of any arbitrarily shaped dislocation in anisotropic cubic media.