Atomistic-to-continuum coupling approximation of a one-dimensional toy model for density functional theory

We consider an atomistic model defined through an interaction field satisfying a variational principle and which can therefore be considered a toy model of (orbital-free) density functional theory. We investigate atomistic-to-continuum coupling mechanisms for this atomistic model, paying special attention to the dependence of the atomistic subproblem on the atomistic region boundary and the boundary conditions. We rigorously prove first-order error estimates for two related coupling mechanisms

A Dynamic Atomistic-Continuum Method for the Simulation of Crystalline Materials

We present a coupled atomistic-continuum method for the modeling of defects and interface dynamics of crystalline materials. The method uses atomistic models such as molecular dynamics near defects and interfaces, and continuum models away from defects and interfaces. We propose a new class of matching conditions between the atomistic and continuum regions. These conditions ensure the accurate passage of large scale information between the atomistic and continuum regions and at the same time minimize the reflection of phonons at the atomistic-continuum interface. They can be made adaptive if we choose appropriate weight functions. We present applications to dislocation dynamics, friction between two-dimensional crystal surfaces and fracture dynamics. We compare results of the coupled method and the detailed atomistic model.Comment: 48 pages, 20 figure

Integrated atomistic process and device simulation of decananometre MOSFETs

In this paper we present a methodology for the integrated atomistic process and device simulation of decananometre MOSFETs. The atomistic process simulations were carried out using the kinetic Monte Carlo process simulator DADOS, which is now integrated into the Synopsys 3D process and device simulation suite Taurus. The device simulations were performed using the Glasgow 3D statistical atomistic simulator, which incorporates density gradient quantum corrections. The overall methodology is illustrated in the atomistic process and device simulation of a well behaved 35 nm physical gate length MOSFET reported by Toshiba

A multi-resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand-binding site

In multi-resolution simulations, different system components are simultaneously modelled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of substrate binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modelled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically, without loss of accuracy in the atomistic description. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity via an exploration of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function

Stability of a force-based hybrid method with planar sharp interface

We study a force-based hybrid method that couples atomistic model with Cauchy-Born elasticity model with sharp transition interface. We identify stability conditions that guarantee the convergence of the hybrid scheme to the solution of the atomistic model with second order accuracy, as the ratio between lattice parameter and the characteristic length scale of the deformation tends to zero. Convergence is established for hybrid schemes with planar sharp interface for system without defects, with general finite range atomistic potential and simple lattice structure. The key ingredient of the proof is regularity and stability analysis of elliptic systems of difference equations. We apply the results to atomistic-to-continuum scheme for a 2D triangular lattice with planar interface.Comment: 21 pages with a 5 page supplementary materia

Consistent Energy-based Atomistic/Continuum Coupling for Two-body Potentials in Three Dimensions

Very few works exist to date on development of a consistent energy-based coupling of atomistic and continuum models of materials in more than one dimension. The difficulty in constructing such a coupling consists in defining a coupled energy whose minimizers are free from uncontrollable errors on the atomistic/continuum interface. In this paper a consistent coupling in three dimensions is proposed. The main achievement of this work is to identify and efficiently treat a modified Cauchy-Born continuum model which can be coupled to the exact atomistic model. The convergence and stability of the method is confirmed with numerical tests.Comment: 29 pages, 1 Matlab code. Typos corrected, exposition improve