Source author record

Mitchell Luskin

Mitchell Luskin 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

39works
14topics
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

39 published item(s)

preprint2022arXiv

Non-Abelian topological defects and strain mapping in 2D moiré materials

We present a general method to analyze the topological nature of the domain boundary connectivity that appeared in relaxed moiré superlattice patterns at the interface of 2-dimensional (2D) van der Waals (vdW) materials. At large enough moiré lengths, all moiré systems relax into commensurated 2D domains separated by networks of dislocation lines. The nodes of the 2D dislocation line network can be considered as vortex-like topological defects. We find that a simple analogy to common topological systems with an $S^1$ order parameter, such as a superconductor or planar ferromagnet, cannot correctly capture the topological nature of these defects. For example, in twisted bilayer graphene, the order parameter space for the relaxed moiré system is homotopy equivalent to a punctured torus. Here, the nodes of the 2D dislocation network can be characterized as elements of the fundamental group of the punctured torus, the free group on two generators, endowing these network nodes with non-Abelian properties. Extending this analysis to consider moiré patterns generated from any relative strain, we find that antivortices occur in the presence of anisotropic heterostrain, such as shear or anisotropic expansion, while arrays of vortices appear under twist or isotropic expansion between vdW materials. Experimentally, utilizing the dark field imaging capability of transmission electron microscopy (TEM), we demonstrate the existence of vortex and antivortex pair formation in a moiré system, caused by competition between different types of heterostrains in the vdW interfaces. We also present a methodology for mapping the underlying heterostrain of a moiré structure from experimental TEM data, which provides a quantitative relation between the various components of heterostrain and vortex-antivortex density in moiré systems.

preprint2022arXiv

Seeing moiré: convolutional network learning applied to twistronics

Moiré patterns made of two-dimensional (2D) materials represent highly tunable electronic Hamiltonians, allowing a wide range of quantum phases to emerge in a single material. Current modeling techniques for moiré electrons requires significant technical work specific to each material, impeding large-scale searches for useful moiré materials. In order to address this difficulty, we have developed a material-agnostic machine learning approach and test it here on prototypical one-dimensional (1D) moiré tight-binding models. We utilize the stacking dependence of the local density of states (SD-LDOS) to convert information about electronic bandstructure into physically relevant images. We then train a neural network that successfully predicts moiré electronic structure from the easily computed SD-LDOS of aligned bilayers. This network can satisfactorily predict moiré electronic structures, even for materials that are not included in its training data.

preprint2021arXiv

Existence of the first magic angle for the chiral model of bilayer graphene

We consider the chiral model of twisted bilayer graphene introduced by Tarnopolsky-Kruchkov-Vishwanath (TKV). TKV have proved that for inverse twist angles $α$ such that the effective Fermi velocity at the moiré $K$ point vanishes, the chiral model has a perfectly flat band at zero energy over the whole Brillouin zone. By a formal expansion, TKV found that the Fermi velocity vanishes at $α\approx .586$. In this work, we give a proof that the Fermi velocity vanishes for at least one $α$ between $.57$ and $.61$ by rigorously justifying TKV's formal expansion of the Fermi velocity over a sufficiently large interval of $α$ values. The idea of the proof is to project the TKV Hamiltonian onto a finite dimensional subspace, and then expand the Fermi velocity in terms of explicitly computable linear combinations of modes in the subspace, while controlling the error. The proof relies on two propositions whose proofs are computer-assisted, i.e., numerical computation together with worst-case estimates on the accumulation of round-off error which show that round-off error cannot possibly change the conclusion of the computation. The propositions give a bound below on the spectral gap of the projected Hamiltonian, an Hermitian $80 \times 80$ matrix whose spectrum is symmetric about $0,$ and verify that two real 18th order polynomials, which approximate the numerator of the Fermi velocity, take values with definite sign when evaluated at specific values of $α$. Together with TKV's work our result proves existence of at least one perfectly flat band of the chiral model.

preprint2020arXiv

Duality between atomic configurations and Bloch states in twistronic materials

The relative orientation (twist) of successive layers of stacked two-dimensional (2D) materials creates variations in the interlayer atomic registry. The variations often form a super lattice, called a moiré pattern, which can alter electronic properties. In this work we introduce a classification of the single-particle electronic structures that can occur in twisted stacks of 2D layers by characterizing them as "moiré molecules" or "moiré crystals". The molecules generate localized electronic states and moiré flat bands, while the crystals are sometimes unconventional and produce electronic banding in the configuration basis. The underpinning of this classification is the duality between interlayer configuration and monolayer Bloch momentum in moiré Hamiltonians. We apply this understanding to diagrams of local electron density in untwisted geometries to produce intuitive and quantitative predictions of twistronic properties. We provide a conceptual introduction to this framework through a one-dimensional model, and then apply it to 2D twisted bilayers of the semi-metal graphene, and of MoS$_2$, a representative material of the transition metal dichalcogenide (TMDC) family of semiconductors. This level of thorough understanding of twistronic phenomena is vital in the search for new material platforms for localized moiré electrons.

preprint2020arXiv

Finite-size effects in wave transmission through plasmonic crystals: A tale of two scales

We study optical coefficients that characterize wave propagation through layered structures called plasmonic crystals. These consist of a finite number of stacked metallic sheets embedded in dielectric hosts with a subwavelength spacing. By adjustment of the frequency, spacing, number as well as geometry of the layers, these structures may exhibit appealing transmission properties in a range of frequencies from the terahertz to the mid-infrared regime. Our approach uses a blend of analytical and numerical methods for the distinct geometries with infinite, translation invariant, flat sheets and nanoribbons. We describe the transmission of plane waves through a plasmonic crystal in comparison to an effective dielectric slab of equal total thickness that emerges from homogenization, in the limit of zero interlayer spacing. We demonstrate numerically that the replacement of the discrete plasmonic crystal by its homogenized counterpart can accurately capture a transmission coefficient akin to the extinction spectrum, even for a relatively small number of layers. We point out the role of a geometry-dependent corrector field, which expresses the effect of subwavelength surface plasmons. In particular, by use of the corrector we describe lateral resonances inherent to the nanoribbon geometry.

preprint2020arXiv

Homogenization of hydrodynamic transport in Dirac fluids

Large-scale electrical and thermal currents in ordinary metals are well approximated by effective medium theory: global transport properties are governed by the solution to homogenized coupled diffusion equations. In some metals, including the Dirac fluid of nearly charge neutral graphene, microscopic transport is not governed by diffusion, but by a more complicated set of linearized hydrodynamic equations, which form a system of degenerate elliptic equations coupled with the Stokes equation for fluid velocity. In sufficiently inhomogeneous media, these hydrodynamic equations reduce to homogenized diffusion equations. We re-cast the hydrodynamic transport equations as the infimum of a functional over conserved currents, and present a functional framework to model and compute the homogenized diffusion tensor relating electrical and thermal currents to charge and temperature gradients. We generalize to this system two well-known results in homogenization theory: Tartar's proof of local convergence to the homogenized theory in periodic and highly oscillatory media, and sub-additivity of the above functional in random media with highly oscillatory, stationary and ergodic coefficients.

preprint2020arXiv

Moiré of Moiré: Modeling Mechanical Relaxation in Incommensurate Trilayer van der Waals Heterostructures

The incommensurate stacking of multi-layered two-dimensional materials is a challenging problem from a theoretical perspective and an intriguing avenue for manipulating their physical properties. Here we present a multi-scale model to obtain the mechanical relaxation pattern of twisted trilayer van der Waals (vdW) heterostructures with two independent twist angles, a generally incommensurate system without a supercell description. We adopt the configuration space as a natural description of such incommensurate layered materials, based on the local environment of atomic positions, bypassing the need for commensurate approximations. To obtain the relaxation pattern, we perform energy minimization with respect to the relaxation displacement vectors. We use a continuum model in combination with the Generalized Stacking Fault energy to describe the interlayer coupling, obtained from first-principles calculations based on Density Functional Theory. We show that the relaxation patterns of twisted trilayer graphene and $\mathrm{WSe_2}$ are "moiré of moiré", as a result of the incommensurate coupling two bilayer moiré patterns. We also show that, in contrast to the symmetry-preserving in-plane relaxation in twisted bilayers, trilayer relaxation can break the two-fold rotational symmetry about the xy-plane when the two twist angles are equal.

preprint2019arXiv

Efficient Computation of Kubo Conductivity for Incommensurate 2D Heterostructures

Here we introduce a numerical method for computing conductivity via the Kubo Formula for incommensurate 2D bilayer heterostructures using a tight-binding framework. We begin with deriving the momentum space formulation and Kubo Formula from the real space tight-binding model using the appropriate Bloch transformation operator. We further discuss the resulting algorithm along with its convergence rate and computation cost in terms of parameters such as relaxation time and temperature. In particular, we show that for low frequencies, low temperature, and long relaxation times conductivity can be computed very efficiently using momentum space for a wide class of materials. We then demonstrate our method by computing conductivity for twisted bilayer graphene (tBLG) for small twist angles.

preprint2019arXiv

Nonretarded edge plasmon-polaritons in anisotropic two-dimensional materials

By an integral equation approach to the time-harmonic classical Maxwell equations, we describe the dispersion in the nonretarded frequency regime of the edge plasmon-polariton (EPP) on a semi-infinite flat sheet. The sheet has an arbitrary, physically admissible, tensor valued and spatially homogeneous conductivity, and serves as a model for a family of two-dimensional conducting materials. We formulate a system of integral equations for the electric field tangential to the sheet in a homogeneous and isotropic ambient medium. We show how this system is simplified via a length scale separation. This view entails the quasi-electrostatic approximation, by which the tangential electric field is replaced by the gradient of a scalar potential, $φ$. By the Wiener-Hopf method, we solve an integral equation for $φ$ in some generality. The EPP dispersion relation comes from the elimination of a divergent limiting Fourier integral for $φ$ at the edge. We connect the existence, or lack thereof, of the EPP dispersion relation to the index for Wiener-Hopf integral equations, an integer of topological character. We indicate that the values of this index may express an asymmetry due to the material anisotropy in the number of wave modes propagating on the sheet away from the edge with respect to the EPP direction of propagation. We discuss extensions such as the setting of two semi-infinite, coplanar sheets. Our theory forms a generalization of the treatment by Volkov and Mikhailov (1988 Sov. Phys. JETP 67 1639).

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.

preprint2016arXiv

Analysis of rippling in incommensurate one-dimensional coupled chains

Graphene and other recently developed 2D materials exhibit exceptionally strong in-plane stiffness. Relaxation of few-layer structures, either free-standing or on slightly mismatched substrates occurs mostly through out-of-plane bending and the creation of large-scale ripples. In this work, we present a novel double chain model, where we allow relaxation to occur by bending of the incommensurate coupled system of chains. As we will see, this model can be seen as a new application of the well-known Frenkel-Kontorova model for a one-dimensional atomic chain lying in a periodic potential. We focus in particular on modeling and analyzing ripples occurring in ground state configurations, as well as their numerical simulation.

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.

preprint2015arXiv

A Theoretical Examination of Diffusive Molecular Dynamics

Diffusive molecular dynamics is a novel model for materials with atomistic resolution that can reach diffusive time scales. The main ideas of diffusive molecular dynamics are to first minimize an approximate variational Gaussian free energy of the system with respect to the mean atomic coordinates (averaging over many vibrational periods), and to then to perform a diffusive step where atoms and vacancies (or two species in a binary alloy) flow on a diffusive time scale via a master equation. We present a mathematical framework for studying this algorithm based upon relative entropy, or Kullback-Leibler divergence. This adds flexibility in how the algorithm is implemented and interpreted. We then compare our formulation, relying on relative entropy and absolute continuity of measures, to existing formulations. The main difference amongst the equations appears in a model for vacancy diffusion, where additional entropic terms appear in our development.

preprint2015arXiv

On solutions of Maxwell's equations with dipole sources over a thin conducting film

We derive and interpret solutions of time-harmonic Maxwell's equations with a vertical and a horizontal electric dipole near a planar, thin conducting film, e.g. graphene sheet, lying between two unbounded isotropic and non-magnetic media. Exact expressions for all field components are extracted in terms of rapidly convergent series of known transcendental functions when the ambient media have equal permittivities and both the dipole and observation point lie on the plane of the film. These solutions are simplified for all distances from the source when the film surface resistivity is large in magnitude compared to the intrinsic impedance of the ambient space. The formulas reveal the analytical structure of two types of waves that can possibly be excited by the dipoles and propagate on the film. One of these waves is intimately related to the surface plasmon-polariton of transverse-magnetic (TM) polarization of plane waves.

preprint2014arXiv

Accuracy of computation of crystalline defects at finite temperature

The present paper aims at developing a theory of computation of crystalline defects at finite temperature. In a one-dimensional setting we introduce Gibbs distributions corresponding to such defects and rigorously establish their asymptotic expansion. We then give an example of using such asymptotic expansion to compare the accuracy of computations using the free boundary conditions and using an atomistic-to-continuum coupling method.

preprint2014arXiv

Analysis of an optimization-based atomistic-to-continuum coupling method for point defects

We formulate and analyze an optimization-based Atomistic-to-Continuum (AtC) coupling method for problems with point defects. Near the defect core the method employs a potential-based atomistic model, which enables accurate simulation of the defect. Away from the core, where site energies become nearly independent of the lattice position, the method switches to a more efficient continuum model. The two models are merged by minimizing the mismatch of their states on an overlap region, subject to the atomistic and continuum force balance equations acting independently in their domains. We prove that the optimization problem is well-posed and establish error estimates.

preprint2014arXiv

Analysis of Transition State Theory Rates upon Spatial Coarse-Graining

Spatial multiscale methods have established themselves as useful tools for extending the length scales accessible by conventional statics (i.e., zero temperature molecular dynamics). Recently, extensions of these methods, such as the finite-temperature quasicontinuum (hot-QC) or Coarse-Grained Molecular Dynamics (CGMD) methods, have allowed for multiscale molecular dynamics simulations at finite temperature. Here, we assess the quality of the long-time dynamics these methods generate by considering canonical transition rates. Specifically, we analyze the transition state theory (TST) rates in CGMD and compare them to the corresponding TST rate of the fully atomistic system. The ability of such an approach to reliably reproduce the TST rate is verified through a relative error analysis, which is then used to highlight the major contributions to the error and guide the choice of degrees of freedom. Finally, our analytical results are compared with numerical simulations for the case of a 1-D chain.

preprint2013arXiv

An Optimization-Based Atomistic-to-Continuum Coupling Method

We present a new optimization-based method for atomistic-to-continuum (AtC) coupling. The main idea is to cast the coupling of the atomistic and continuum models as a constrained optimization problem with virtual Dirichlet controls on the interfaces between the atomistic and continuum subdomains. The optimization objective is to minimize the error between the atomistic and continuum solutions on the overlap between the two subdomains, while the atomistic and continuum force balance equations provide the constraints. Splitting of the atomistic and continuum problems instead of blending them and their subsequent use as constraints in the optimization problem distinguishes our approach from the existing AtC formulations. We present and analyze the method in the context of a one-dimensional chain of atoms modeled using a linearized two-body next-nearest neighbor interactions.

preprint2013arXiv

Development of an Optimization-Based Atomistic-to-Continuum Coupling Method

Atomistic-to-Continuum (AtC) coupling methods are a novel means of computing the properties of a discrete crystal structure, such as those containing defects, that combine the accuracy of an atomistic (fully discrete) model with the efficiency of a continuum model. In this note we extend the optimization-based AtC, formulated in arXiv:1304.4976 for linear, one-dimensional problems to multi-dimensional settings and arbitrary interatomic potentials. We conjecture optimal error estimates for the multidimensional AtC, outline an implementation procedure, and provide numerical results to corroborate the conjecture for a 1D Lennard-Jones system with next-nearest neighbor interactions.

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.

preprint2012arXiv

Numerical Analysis of Parallel Replica Dynamics

Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of their time about the minima of the underlying potential, occasionally transitioning into different basins of attraction. The essential idea of parallel replica dynamics is that the exit time distribution from a given well for a single process can be approximated by the minimum of the exit time distributions of $N$ independent identical processes, each run for only 1/N-th the amount of time. While promising, this leads to a series of numerical analysis questions about the accuracy of the exit distributions. Building upon the recent work in Le Bris et al., we prove a unified error estimate on the exit distributions of the algorithm against an unaccelerated process. Furthermore, we study a dephasing mechanism, and prove that it will successfully complete.

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.

preprint2011arXiv

Analysis of Energy-Based Blended Quasicontinuum Approximations

The development of patch test consistent quasicontinuum energies for multi-dimensional crystalline solids modeled by many-body potentials remains a challenge. The original quasicontinuum energy (QCE) has been implemented for many-body potentials in two and three space dimensions, but it is not patch test consistent. We propose that by blending the atomistic and corresponding Cauchy-Born continuum models of QCE in an interfacial region with thickness of a small number $k$ of blended atoms, a general quasicontinuum energy (BQCE) can be developed with the potential to significantly improve the accuracy of QCE near lattice instabilities such as dislocation formation and motion. In this paper, we give an error analysis of the blended quasicontinuum energy (BQCE) for a periodic one-dimensional chain of atoms with next-nearest neighbor interactions. Our analysis includes the optimization of the blending function for an improved convergence rate. We show that the $\ell^2$ strain error for the non-blended QCE energy (QCE), which has low order $\text{O}(ε^{1/2})$ where $ε$ is the atomistic length scale, can be reduced by a factor of $k^{3/2}$ for an optimized blending function where $k$ is the number of atoms in the blending region. The QCE energy has been further shown to suffer from a O$(1)$ error in the critical strain at which the lattice loses stability. We prove that the error in the critical strain of BQCE can be reduced by a factor of $k^2$ for an optimized blending function, thus demonstrating that the BQCE energy for an optimized blending function has the potential to give an accurate approximation of the deformation near lattice instabilities such as crack growth.

preprint2011arXiv

Analysis of the quasi-nonlocal approximation of linear and circular chains in the plane

We give an analysis of the stability and displacement error for linear and circular atomistic chains in the plane when the atomistic energy is approximated by the Cauchy-Born continuum energy and by the quasi-nonlocal atomistic-to-continuum coupling energy. We consider atomistic energies that include Lennard-Jones type nearest neighbor and next nearest neighbor pair-potential interactions. Previous analyses for linear chains have shown that the Cauchy-Born and quasi-nonlocal approximations reproduce (up to the order of the lattice spacing) the atomistic lattice stability for perturbations that are constrained to the line of the chain. However, we show that the Cauchy-Born and quasi-nonlocal approximations give a finite increase for the lattice stability of a linear or circular chain under compression when general perturbations in the plane are allowed. We also analyze the increase of the lattice stability under compression when pair-potential energies are augmented by bond-angle energies. Our estimates of the largest strain for lattice stability (the critical strain) are sharp (exact up to the order of the lattice scale). We then use these stability estimates and modeling error estimates for the linearized Cauchy-Born and quasi-nonlocal energies to give an optimal order (in the lattice scale) {\em a priori} error analysis for the approximation of the atomistic strain in $\ell^2_ε$ due to an external force.

preprint2011arXiv

Lattice Stability for Atomistic Chains Modeled by Local Approximations of the Embedded Atom Method

The accurate approximation of critical strains for lattice instability is a key criterion for predictive computational modeling of materials. In this paper, we present a comparison of the lattice stability for atomistic chains modeled by the embedded atom method (EAM) with their approximation by local Cauchy-Born models. We find that both the volume-based local model and the reconstruction-based local model can give O(1) errors for the critical strain since the embedding energy density is generally strictly convex. The critical strain predicted by the volume-based model is always larger than that predicted by the atomistic model, but the critical strain for reconstruction-based models can be either larger or smaller than that predicted by the atomistic model.

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.

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.

preprint2010arXiv

A Generalized Quasi-Nonlocal Atomistic-to-Continuum Coupling Method with Finite Range Interaction

The accurate and efficient computation of the deformation of crystalline solids requires the coupling of atomistic models near lattice defects such as cracks and dislocations with coarse-grained models away from the defects. Quasicontinuum methods utilize a strain energy density derived from the Cauchy-Born rule for the coarse-grained model. Several quasicontinuum methods have been proposed to couple the atomistic model with the Cauchy-Born strain energy density. The quasi-nonlocal coupling method is easy to implement and achieves a reasonably accurate coupling for short range interactions. In this paper, we give a new formulation of the quasi-nonlocal method in one space dimension that allows its extension to arbitrary finite range interactions. We also give an analysis of the stability and accuracy of a linearization of our generalized quasi-nonlocal method that holds for strains up to lattice instabilities.

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.

preprint2010arXiv

An Analysis of the Quasi-Nonlocal Quasicontinuum Approximation of the Embedded Atom Model

The quasi-nonlocal quasicontinuum method (QNL) is a consistent hybrid coupling method for atomistic and continuum models. Embedded atom models are empirical many-body potentials that are widely used for FCC metals such as copper and aluminum. In this paper, we consider the QNL method for EAM potentials, and we give a stability and error analysis for a chain with next-nearest neighbor interactions. We identify conditions for the pair potential, electron density function, and embedding function so that the lattice stability of the atomistic and the EAM-QNL models are asymptotically equal.

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.

preprint2009arXiv

An Optimal Order Error Analysis of the One-Dimensional Quasicontinuum Approximation

We derive a model problem for quasicontinuum approximations that allows a simple, yet insightful, analysis of the optimal-order convergence rate in the continuum limit for both the energy-based quasicontinuum approximation and the quasi-nonlocal quasicontinuum approximation. The optimal-order error estimates for the quasi-nonlocal quasicontinuum approximation are given for all strains up to the continuum limit strain for fracture. The analysis is based on an explicit treatment of the coupling error at the atomistic to continuum interface, combined with an analysis of the error due to atomistic and continuum schemes using the stability of the quasicontinuum approximation.

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.

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.

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.

preprint2008arXiv

An Analysis of the Effect of Ghost Force Oscillation on Quasicontinuum Error

The atomistic to continuum interface for quasicontinuum energies exhibits nonzero forces under uniform strain that have been called ghost forces. In this paper, we prove for a linearization of a one-dimensional quasicontinuum energy around a uniform strain that the effect of the ghost forces on the displacement nearly cancels and has a small effect on the error away from the interface. We give optimal order error estimates that show that the quasicontinuum displacement converges to the atomistic displacement at the optimal rate O($h$) in the discrete $\ell^\infty$ norm and O($h^{1/p}$) in the $w^{1,p}$ norm for $1 \leq p < \infty.$ where $h$ is the interatomic spacing. We also give a proof that the error in the displacement gradient decays away from the interface to O($h$) at distance O($h|\log h|$) in the atomistic region and distance O($h$) in the continuum region. E, Ming, and Yang previously gave a counterexample to convergence in the $w^{1,\infty}$ norm for a harmonic interatomic potential. Our work gives an explicit and simplified form for the decay of the effect of the atomistic to continuum coupling error in terms of a general underlying interatomic potential and gives the estimates described above in the discrete $\ell^\infty$ and $w^{1,p}$ norms.

preprint2008arXiv

Non-ergodicity of Nose-Hoover dynamics

The numerical integration of the Nose-Hoover dynamics gives a deterministic method that is used to sample the canonical Gibbs measure. The Nose-Hoover dynamics extends the physical Hamiltonian dynamics by the addition of a "thermostat" variable, that is coupled nonlinearly with the physical variables. The accuracy of the method depends on the dynamics being ergodic. Numerical experiments have been published earlier that are consistent with non-ergodicity of the dynamics for some model problems. The authors recently proved the non-ergodicity of the Nose-Hoover dynamics for the one-dimensional harmonic oscillator. In this paper, this result is extended to non-harmonic one-dimensional systems. It is also shown for some multidimensional systems that the averaged dynamics for the limit of infinite thermostat "mass" have many invariants, thus giving theoretical support for either non-ergodicity or slow ergodization. Numerical experiments for a two-dimensional central force problem and the one-dimensional pendulum problem give evidence for non-ergodicity.

preprint2007arXiv

Goal-Oriented Adaptive Mesh Refinement for the Quasicontinuum Approximation of a Frenkel-Kontorova Model

The quasicontinuum approximation is a method to reduce the atomistic degrees of freedom of a crystalline solid by piecewise linear interpolation from representative atoms that are nodes for a finite element triangulation. In regions of the crystal with a highly nonuniform deformation such as around defects, every atom must be a representative atom to obtain sufficient accuracy, but the mesh can be coarsened away from such regions to remove atomistic degrees of freedom while retaining sufficient accuracy. We present an error estimator and a related adaptive mesh refinement algorithm for the quasicontinuum approximation of a generalized Frenkel-Kontorova model that enables a quantity of interest to be efficiently computed to a predetermined accuracy.

preprint2006arXiv

Analysis of a Force-Based Quasicontinuum Approximation

We analyze a force-based quasicontinuum approximation to a one-dimensional system of atoms that interact by a classical atomistic potential. This force-based quasicontinuum approximation is derived as the modification of an energy-based quasicontinuum approximation by the addition of nonconservative forces to correct nonphysical ``ghost'' forces that occur in the atomistic to continuum interface. We prove that the force-based quasicontinuum equations have a unique solution under suitable restrictions on the loads. For Lennard-Jones next-nearest-neighbor interactions, we show that unique solutions exist for loads in a symmetric region extending nearly to the tensile limit. We give an analysis of the convergence of the ghost force iteration method to solve the equilibrium equations for the force-based quasicontinuum approximation. We show that the ghost force iteration is a contraction and give an analysis for its convergence rate.