BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Christoph Ortner

Christoph Ortner contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath-phmath.APmath.MPcond-mat.mtrl-sciphysics.comp-phcond-mat.other

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
8topics
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

21 published item(s)

preprint2024arXiv

Atomic Cluster Expansion without Self-Interaction

The Atomic Cluster Expansion (ACE) (Drautz, Phys. Rev. B 99, 2019) has been widely applied in high energy physics, quantum mechanics and atomistic modeling to construct many-body interaction models respecting physical symmetries. Computational efficiency is achieved by allowing non-physical self-interaction terms in the model. We propose and analyze an efficient method to evaluate and parameterize an orthogonal, or, non-self-interacting cluster expansion model. We present numerical experiments demonstrating improved conditioning and more robust approximation properties than the original expansion in regression tasks both in simplified toy problems and in applications in the machine learning of interatomic potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A framework for a generalisation analysis of machine-learned interatomic potentials

Machine-learned interatomic potentials (MLIPs) and force fields (i.e. interaction laws for atoms and molecules) are typically trained on limited data-sets that cover only a very small section of the full space of possible input structures. MLIPs are nevertheless capable of making accurate predictions of forces and energies in simulations involving (seemingly) much more complex structures. In this article we propose a framework within which this kind of generalisation can be rigorously understood. As a prototypical example, we apply the framework to the case of simulating point defects in a crystalline solid. Here, we demonstrate how the accuracy of the simulation depends explicitly on the size of the training structures, on the kind of observations (e.g., energies, forces, force constants, virials) to which the model has been fitted, and on the fit accuracy. The new theoretical insights we gain partially justify current best practices in the MLIP literature and in addition suggest a new approach to the collection of training data and the design of loss functions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Atomic cluster expansion and wave function representations

The atomic cluster expansion (ACE) has been highly successful for the parameterisation of symmetric (invariant or equivariant) properties of many-particle systems. Here, we generalize its derivation to anti-symmetric functions. We show how numerous well-known linear representations of wave functions naturally arise within this framework and we explore how recent successful nonlinear parameterisations can be further enhanced by employing ACE methodology. From this analysis we propose a wide design space of promising wave function representations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Equivariant analytical mapping of first principles Hamiltonians to accurate and transferable materials models

We propose a scheme to construct predictive models for Hamiltonian matrices in atomic orbital representation from ab initio data as a function of atomic and bond environments. The scheme goes beyond conventional tight binding descriptions as it represents the ab initio model to full order, rather than in two-centre or three-centre approximations. We achieve this by introducing an extension to the Atomic Cluster Expansion (ACE) descriptor that represents Hamiltonian matrix blocks that transform equivariantly with respect to the full rotation group. The approach produces analytical linear models for the Hamiltonian and overlap matrices. Through an application to aluminium, we demonstrate that it is possible to train models from a handful of structures computed with density functional theory, and apply them to produce accurate predictions for the electronic structure. The model generalises well and is able to predict defects accurately from only bulk training data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal Evaluation of Symmetry-Adapted $n$-Correlations Via Recursive Contraction of Sparse Symmetric Tensors

We present a comprehensive analysis of an algorithm for evaluating high-dimensional polynomials that are invariant under permutations and rotations. The key bottleneck is the contraction of a high-dimensional symmetric and sparse tensor with a specific sparsity pattern that is directly related to the symmetries imposed on the polynomial. We propose an explicit construction of a recursive evaluation strategy and show that it is optimal in the limit of infinite polynomial degree.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stability of the Minimum Energy Path

The minimum energy path (MEP) is the most probable transition path that connects two equilibrium states of a potential energy landscape. It has been widely used to study transition mechanisms as well as transition rates in the fields of chemistry, physics, and materials science. In this paper, we derive a novel result establishing the stability of MEPs under perturbations of the energy landscape. The result also represents a crucial step towards studying the convergence of various numerical approximations of MEPs, such as the nudged elastic band and string methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Theoretical Study of Elastic Far-Field Decay from Dislocations in Multilattices

We precisely and rigorously characterise the decay of elastic fields generated by dislocations in crystalline materials, focusing specifically on the role of multilattices. Concretely, we establish that the elastic field generated by a dislocation in a multilattice can be decomposed into a continuum field predicted by a linearised Cauchy-Born elasticity theory, and a discrete and nonlinear core corrector representing the defect core. We demonstrate both analytically and numerically the consequences of this result for cell size effects in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Asymptotic Expansion of the Elastic Far-Field of a Crystalline Defect

Lattice defects in crystalline materials create long-range elastic fields which can be modelled on the atomistic scale using an infinite system of discrete nonlinear force balance equations. Starting with these equations, this work rigorously derives a novel far-field expansion of these fields: The expansion is computable and is expressed as a sum of continuum correctors and discrete multipole terms which decay with increasing algebraic rate as the order of the expansion increases. Truncating the expansion leaves a remainder describing the defect core structure, which is localised in the sense that it decays with an algebraic rate corresponding to the order at which the truncation occurred.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Rigorous body-order approximations of an electronic structure potential energy landscape

We show that the local density of states (LDOS) of a wide class of tight-binding models has a weak body-order expansion. Specifically, we prove that the resulting body-order expansion for analytic observables such as the electron density or the energy has an exponential rate of convergence both at finite Fermi-temperature as well as for insulators at zero Fermi-temperature. We discuss potential consequences of this observation for modelling the potential energy landscape, as well as for solving the electronic structure problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Posteriori Error Estimates for Adaptive QM/MM Coupling Methods

Hybrid quantum/molecular mechanics models (QM/MM methods) are widely used in material and molecular simulations when MM models do not provide sufficient accuracy but pure QM models are computationally prohibitive. Adaptive QM/MM coupling methods feature on-the-fly classification of atoms during the simulation, allowing the QM and MM subsystems to be updated as needed. In this work, we propose such an adaptive QM/MM method for material defect simulations based on a new residual based it a posteriori error estimator, which provides both lower and upper bounds for the true error. We validate the analysis and illustrate the effectiveness of the new scheme on numerical simulations for material defects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of cell size effects in atomistic crack propagation

We consider crack propagation in a crystalline material in terms of bifurcation analysis. We provide evidence that the stress intensity factor is a natural bifurcation parameter, and that the resulting bifurcation diagram is a periodic "snaking curve". We then prove qualitative properties of the equilibria and convergence rates of finite-cell approximations to the "exact" bifurcation diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Point Defects in Tight Binding Models for Insulators

We consider atomistic geometry relaxation in the context of linear tight binding models for point defects. A limiting model as Fermi-temperature is sent to zero is formulated, and an exponential rate of convergence for the nuclei configuration is established. We also formulate the thermodynamic limit model at zero Fermi-temperature, extending the results of [H. Chen, J. Lu, C. Ortner. Arch. Ration. Mech. Anal., 2018]. We discuss the non-trivial relationship between taking zero temperature and thermodynamic limits in the finite Fermi-temperature models.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Locality of Interatomic Forces in Tight Binding Models for Insulators

The tight binding model is a minimalistic electronic structure model for predicting properties of materials and molecules. For insulators at zero Fermi-temperature we show that the potential energy surface of this model can be decomposed into exponentially localised site energy contributions, thus providing qualitatively sharp estimates on the interatomic interaction range which justifies a range of multi-scale models. For insulators at finite Fermi-temperature we obtain locality estimates that are uniform in the zero-temperature limit. A particular feature of all our results is that they depend only weakly on the point spectrum. Numerical tests confirm our analytical results. This work extends and strengthens (Chen, Ortner 2016) and (Chen, Lu, Ortner 2018) for finite temperature models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of an Energy-based Atomistic/Continuum Coupling Approximation of a Vacancy in the 2D Triangular Lattice

We present a comprehensive a priori error analysis of a practical energy based atomistic/continuum coupling method (Shapeev, arXiv:1010.0512) in two dimensions, for finite-range pair-potential interactions, in the presence of vacancy defects. The majority of the work is devoted to the analysis of consistency and stability of the method. These yield a priori error estimates in the H1-norm and the energy, which depend on the mesh size and the "smoothness" of the atomistic solution in the continuum region. Based on these error estimates, we present heuristics for an optimal choice of the atomistic region and the finite element mesh, which yields convergence rates in terms of the number of degrees of freedom. The analytical predictions are supported by extensive numerical tests.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Linear Stationary Iterative Methods for the Force-based Quasicontinuum Approximation

Force-based multiphysics coupling methods have become popular since they provide a simple and efficient coupling mechanism, avoiding the difficulties in formulating and implementing a consistent coupling energy. They are also the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. However, the development of efficient and reliable iterative solution methods for the force-based approximation presents a challenge due to the non-symmetric and indefinite structure of the linearized force-based quasicontinuum approximation, as well as to its unusual stability properties. In this paper, we present rigorous numerical analysis and computational experiments to systematically study the stability and convergence rate for a variety of linear stationary iterative methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The role of the patch test in 2D atomistic-to-continuum coupling methods

For a general class of atomistic-to-continuum coupling methods, coupling multi-body interatomic potentials with a P1-finite element discretisation of Cauchy--Born nonlinear elasticity, this paper adresses the question whether patch test consistency (or, absence of ghost forces) implies a first-order error estimate. In two dimensions it is shown that this is indeed true under the following additional technical assumptions: (i) an energy consistency condition, (ii) locality of the interface correction, (iii) volumetric scaling of the interface correction, and (iv) connectedness of the atomistic region. The extent to which these assumptions are necessary is discussed in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Iterative Methods for the Force-based Quasicontinuum Approximation

Force-based atomistic-continuum hybrid methods are the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. For this reason, and due to their algorithmic simplicity, force-based coupling methods have become a popular class of atomistic-continuum hybrid models as well as other types of multiphysics models. However, the recently discovered unusual stability properties of the linearized force-based quasicontinuum (QCF) approximation, especially its indefiniteness, present a challenge to the development of efficient and reliable iterative methods. We present analytic and computational results for the generalized minimal residual (GMRES) solution of the linearized QCF equilibrium equations. We show that the GMRES method accurately reproduces the stability of the force-based approximation and conclude that an appropriately preconditioned GMRES method results in a reliable and efficient solution method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Spectrum of the Force-Based Quasicontinuum Operator for a Homogeneous Periodic Chain

We show under general conditions that the linearized force-based quasicontinuum (QCF) operator has a positive spectrum, which is identical to the spectrum of the quasinonlocal quasicontinuum (QNL) operator in the case of second-neighbour interactions. Moreover, we establish a bound on the condition number of a matrix of eigenvectors that is uniform in the number of atoms and the size of the atomistic region. These results establish the validity of and improve upon recent conjectures ([arXiv:0907.3861, Conjecture 2] and [arXiv:0910.2013, Conjecture 8]) which were based on numerical experiments. As immediate consequences of our results we obtain rigorous estimates for convergence rates of (preconditioned) GMRES algorithms, as well as a new stability estimate for the QCF method.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Sharp Stability Estimates for the Force-based Quasicontinuum Method

A sharp stability analysis of atomistic-to-continuum coupling methods is essential for evaluating their capabilities for predicting the formation and motion of lattice defects. We formulate a simple one-dimensional model problem and give a detailed analysis of the stability of the force-based quasicontinuum (QCF) method. The focus of the analysis is the question whether the QCF method is able to predict a critical load at which fracture occurs. Numerical experiments show that the spectrum of a linearized QCF operator is identical to the spectrum of a linearized energy-based quasi-nonlocal quasicontinuum operator (QNL), which we know from our previous analyses to be positive below the critical load. However, the QCF operator is non-normal and it turns out that it is not generally positive definite, even when all of its eigenvalues are positive. Using a combination of rigorous analysis and numerical experiments, we investigate in detail for which choices of "function spaces" the QCF operator is stable, uniformly in the size of the atomistic system. Force-based multi-physics coupling methods are popular techniques to circumvent the difficulties faced in formulating consistent energy-based coupling pproaches. Even though the QCF method is possibly the simplest coupling method of this kind, we anticipate that many of our observations apply more generally.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

An Analysis of Node-Based Cluster Summation Rules in the Quasicontinuum Method

We investigate two examples of node-based cluster summation rules that have been proposed for the quasicontinuum method: a force-based approach (Knap & Ortiz, J. Mech. Phys. Solids 49, 2001), and an energy-based approach which is a generalization of the non-local quasicontinuum method (Eidel & Stukowski, J. Mech. Phys. Solids, to appear). We show that, even for the case of nearest neighbour interaction in a one-dimensional periodic chain, both of these approaches create large errors when used with graded and, more generally, non-smooth meshes. These errors cannot be removed by increasing the cluster size. We offer some suggestions how the accuracy of (cluster) summation rules may be improved.

0 citations0 reviewsNo code linkedOpen access