BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christoph Ortner

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

ResearcherUnclaimed source record
math.NAmath.APcond-mat.mtrl-scimath-phmath.MPNumerical Analysisphysics.comp-phcond-mat.othermath.OCmath.PR

Catalog footprint

What is connected

41works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

41 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
preprint2016arXiv

A universal preconditioner for simulating condensed phase materials

We introduce a universal sparse preconditioner that accelerates geometry optimisation and saddle point search tasks that are common in the atomic scale simulation of materials. Our preconditioner is based on the neighbourhood structure and we demonstrate the gain in computational efficiency in a wide range of materials that include metals, insulators and molecular solids. The simple structure of the preconditioner means that the gains can be realised in practice not only when using expensive electronic structure models but also for fast empirical potentials. Even for relatively small systems of a few hundred atoms, we observe speedups of a factor of two or more, and the gain grows with system size. An open source Python implementation within the Atomic Simulation Environment is available, offering interfaces to a wide range of atomistic codes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Analysis of a Predictor-Corrector Method for Computationally Efficient Modeling of Surface Effects in 1D

The regular Cauchy--Born method is a useful and efficient tool for analyzing bulk properties of materials in the absence of defects. However, the method normally fails to capture surface effects, which are essential to determining material properties at small length scales. In this paper, we present a corrector method that improves upon the prediction for material behavior from the Cauchy--Born method over a small boundary layer at the surface of a 1D material by capturing the missed surface effects. We justify the separation of the problem into a bulk response and a localized surface correction by establishing an error estimate, which vanishes in the long wavelength limit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence and Cycling in Walker-type Saddle Search Algorithms

Algorithms for computing local minima of smooth objective functions enjoy a mature theory as well as robust and efficient implementations. By comparison, the theory and practice of saddle search is destitute. In this paper we present results for idealized versions of the dimer and gentlest ascent (GAD) saddle search algorithms that show-case the limitations of what is theoretically achievable within the current class of saddle search algorithms: (1) we present an improved estimate on the region of attraction of saddles; and (2) we construct quasi-periodic solutions which indicate that it is impossible to obtain globally convergent variants of dimer and GAD type algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electronic Density of States for Incommensurate Layers

We prove that the electronic density of states (DOS) for 2D incommensurate layered structures, where Bloch theory does not apply, is well-defined as the thermodynamic limit of finite clusters. In addition, we obtain an explicit representation formula for the DOS as an integral over local configurations. Next, based on this representation formula, we propose a novel algorithm for computing electronic structure properties in incommensurate heterostructures, which overcomes limitations of the common approach to artificially strain a large supercell and then apply Bloch theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On assessing the accuracy of defect free energy computations

We develop a rigorous error analysis for coarse-graining of defect-formation free energy. For a one-dimensional constrained atomistic system, we establish the thermodynamic limit of the defect-formation free energy and obtain explicitly the rate of convergence. We then construct a sequence of coarse-grained energies with the same rate but significantly reduced computational cost. We illustrate our analytical results through explicit computations for the case of harmonic potentials and through numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Regularity and Locality of Point Defects in Multilattices

We formulate a model for a point defect embedded in a homogeneous multilattice crystal with an empirical interatomic potential interaction. Under a natural, phonon stability assumption we quantify the decay of the long-range elastic fields with increasing distance from the defect. These decay estimates are an essential ingredient in quantifying approximation errors in coarse-grained models and in the construction of optimal numerical methods for approximating crystalline defects.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

QM/MM methods for crystalline defects. Part 1: Locality of the tight binding model

The tight binding model is a minimal electronic structure model for molecular modelling and simulation. We show that the total energy in this model can be decomposed into site energies, that is, into contributions from each atomic site whose influence on their environment decays exponentially. This result lays the foundation for a rigorous analysis of QM/MM coupling schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

QM/MM methods for crystalline defects. Part 2: Consistent energy and force-mixing

QM/MM hybrid methods employ accurate quantum (QM) models only in regions of interest (defects) and switch to computationally cheaper interatomic potential (MM) models to describe the crystalline bulk. We develop two QM/MM hybrid methods for crystalline defect simulations, an energy-based and a force-based formulation, employing a tight binding QM model. Both methods build on two principles: (i) locality of the QM model; and (ii) constructing the MM model as an explicit and controllable approximation of the QM model. This approach enables us to establish explicit convergence rates in terms of the size of QM region.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of Blended Atomistic/Continuum Hybrid Methods

We present a comprehensive error analysis of two prototypical atomistic-to-continuum coupling methods of blending type: the energy-based and the force-based quasicontinuum methods. Our results are valid in two and three dimensions, for finite range many-body interactions (e.g., EAM type), and in the presence of lattice defects (we consider point defects and dislocations). The two key ingredients in the analysis are (i) new force and energy consistency error estimates; and (ii) a new technique for proving energy norm stability of a/c couplings that requires only the assumption that the exact atomistic solution is a stable equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of stable screw dislocation configurations in an anti--plane lattice model

We consider a variational anti-plane lattice model and demonstrate that at zero temperature, there exist locally stable states containing screw dislocations, given conditions on the distance between the dislocations, and on the distance between dislocations and the boundary of the crystal. In proving our results, we introduce approximate solutions which are taken from the theory of dislocations in linear elasticity, and use the inverse function theorem to show that local minimisers lie near them. This gives credence to the commonly held intuition that linear elasticity is essentially correct up to a few spacings from the dislocation core.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

(In-)Stability and Stabilisation of QNL-Type Atomistic-to-Continuum Coupling Methods

We study the stability of ghost force-free energy-based atomistic-to-continuum coupling methods. In 1D we essentially complete the theory by introducing a universally stable a/c coupling as well as a stabilisation mechanism for unstable coupling schemes. We then present a comprehensive study of a two-dimensional scalar planar interface setting, as a step towards a general 2D/3D vectorial analysis. Our results point out various new challenges. For example, we find that none of the ghost force-free methods known to us are universally stable (i.e., stable under general interaction and general loads). We then explore to what extent our 1D stabilisation mechanism can be extended.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Existence and stability of a screw dislocation under anti-plane deformation

We formulate a variational model for a geometrically necessary screw dislocation in an anti-plane lattice model at zero temperature. Invariance of the energy functional under lattice symmetries renders the problem non-coercive. Nevertheless, by establishing coercivity with respect to the elastic strain and a concentration compactness principle, we prove existence of a global energy minimiser and thus demonstrate that dislocations are globally stable equilibria within our model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theory-based Benchmarking of the Blended Force-Based Quasicontinuum Method

We formulate an atomistic-to-continuum coupling method based on blending atomistic and continuum forces. Our precise choice of blending mechanism is informed by theoretical predictions. We present a range of numerical experiments studying the accuracy of the scheme, focusing in particular on its stability. These experiments confirm and extend the theoretical predictions, and demonstrate a superior accuracy of B-QCF over energy-based blending schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Note on Linear Elliptic Systems on $\R^d$

We are concerned with the well-posedness of linear elliptic systems posed on $\mathbb{R}^d$. The concrete problem of interest, for which we require this theory, arises from the linearization of the equations of anisotropic finite elasticity, however, our results are more generally applicable to translation-invariant problem posed on $\mathbb{R}^d$. We describe a variant of homogeneous Sobolev spaces, which are convenient for treating problems of this kind.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Symmetries of 2-lattices and second order accuracy of the Cauchy--Born Model

We show that the Cauchy--Born model of a single-species 2-lattice is second order if the atomistic and continuum kinematics are connected in a novel way. Our proof uses a generalization to 2-lattices of the point symmetry of Bravais lattices. Moreover, by identifying similar symmetries in multi-species pair interaction models, we construct a new stored energy density, using shift-gradients but not strain gradients, that is also second order accurate. These results can be used to develop highly accurate continuum models and atomistic/continuum coupling methods for materials such as graphene, hcp metals, and shape memory alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Computational and Theoretical Investigation of the Accuracy of Quasicontinuum Methods

We give computational results to study the accuracy of several quasicontinuum methods for two benchmark problems - the stability of a Lomer dislocation pair under shear and the stability of a lattice to plastic slip under tensile loading. We find that our theoretical analysis of the accuracy near instabilities for one-dimensional model problems can successfully explain most of the computational results for these multi-dimensional benchmark problems. However, we also observe some clear discrepancies, which suggest the need for additional theoretical analysis and benchmark problems to more thoroughly understand the accuracy of quasicontinuum methods.

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

Construction and sharp consistency estimates for atomistic/continuum coupling methods with general interfaces: a 2D model problem

We present a new variant of the geometry reconstruction approach for the formulation of atomistic/continuum coupling methods (a/c methods). For multi-body nearest-neighbour interactions on the 2D triangular lattice, we show that patch test consistent a/c methods can be constructed for arbitrary interface geometries. Moreover, we prove that all methods within this class are first-order consistent at the atomistic/continuum interface and second-order consistent in the interior of the continuum region.

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

Positive-Definiteness of the Blended Force-Based Quasicontinuum Method

The development of consistent and stable quasicontinuum models for multi-dimensional crystalline solids remains a challenge. For example, proving stability of the force-based quasicontinuum (QCF) model remains an open problem. In 1D and 2D, we show that by blending atomistic and Cauchy--Born continuum forces (instead of a sharp transition as in the QCF method) one obtains positive-definite blended force-based quasicontinuum (B-QCF) models. We establish sharp conditions on the required blending width.

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

Accuracy of Quasicontinuum Approximations Near Instabilities

The formation and motion of lattice defects such as cracks, dislocations, or grain boundaries, occurs when the lattice configuration loses stability, that is, when an eigenvalue of the Hessian of the lattice energy functional becomes negative. When the atomistic energy is approximated by a hybrid energy that couples atomistic and continuum models, the accuracy of the approximation can only be guaranteed near deformations where both the atomistic energy as well as the hybrid energy are stable. We propose, therefore, that it is essential for the evaluation of the predictive capability of atomistic-to-continuum coupling methods near instabilities that a theoretical analysis be performed, at least for some representative model problems, that determines whether the hybrid energies remain stable {\em up to the onset of instability of the atomistic energy}. We formulate a one-dimensional model problem with nearest and next-nearest neighbor interactions and use rigorous analysis, asymptotic methods, and numerical experiments to obtain such sharp stability estimates for the basic conservative quasicontinuum (QC) approximations. Our results show that the consistent quasi-nonlocal QC approximation correctly reproduces the stability of the atomistic system, whereas the inconsistent energy-based QC approximation incorrectly predicts instability at a significantly reduced applied load that we describe by an analytic criterion in terms of the derivatives of the atomistic potential.

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
preprint2009arXiv

Stability, Instability, and Error of the Force-based Quasicontinuum Approximation

Due to their algorithmic simplicity and high accuracy, force-based model coupling techniques are an exciting development in computational physics. For example, the force-based quasicontinuum approximation is the only known pointwise consistent quasicontinuum approximation for coupling a general atomistic model with a finite element continuum model. In this paper, we analyze the stability of the force-based quasicontinuum approximation. We then use our stability result to obtain an optimal order error analysis of this coupling method that provides theoretical justification for the high accuracy of the force-based quasicontinuum approximation -- the computational efficiency of continuum modeling can be utilized without the loss of significant accuracy if defects are captured in the atomistic region. The main challenge we need to overcome is the fact (which we prove) that the linearized quasicontinuum operator is typically not positive definite. Moreover, we prove that no uniform inf-sup stability condition holds for discrete versions of the $W^{1,p}$-$W^{1,q}$ "duality pairing" with $1/p+1/q=1$, if $1 \leq p < \infty$. We must therefore derive an inf-sup stability condition for a discrete version of the $W^{1,\infty}$-$W^{1,1}$ "duality pairing" which then leads to optimal order error estimates in a discrete $W^{1,\infty}$-norm.

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