Atomistically-informed continuum modeling and isogeometric analysis of 2D materials over holey substrates

Abstract

This work develops, discretizes, and validates a continuum model of a molybdenum disulfide (MoS${}_2$) monolayer interacting with a periodic holey silicon nitride substrate via van der Waals (vdW) forces. The MoS${}_2$ layer is modeled as a geometrically nonlinear Kirchhoff-Love shell, and vdW forces are modeled by a Lennard-Jones potential, simplified using approximations for a smooth substrate topography. The material parameters of the shell model are calibrated by comparing small-strain tensile and bending tests with atomistic simulations. This model is efficiently discretized using isogeometric analysis (IGA) for the shell structure and a pseudo-time continuation method for energy minimization. The IGA shell model is validated against fully-atomistic calculations for several benchmark problems with different substrate geometries. The continuum simulations reproduce deflections, strains and curvatures predicted by atomistic simulations, which are known to strongly affect the electronic properties of MoS${}_2$, with deviations well below the modeling errors suggested by differences between the widely-used reactive empirical bond order and Stillinger-Weber interatomic potentials. Agreement with atomistic results depends on geometric nonlinearity in some cases, but a simple isotropic St. Venant-Kirchhoff model is found to be sufficient to represent material behavior. We find that the IGA discretization of the continuum model has a much lower computational cost than atomistic simulations, and expect that it will enable efficient design space exploration in strain engineering applications. This is demonstrated by studying the dependence of strain and curvature in MoS${}_2$ over a holey substrate as a function of the hole spacing on scales inaccessible to atomistic calculations. The results show an unexpected qualitative change in the deformation pattern below a critical hole separation.

References in corpus (14)

• Electric Field Effect in Atomically Thin Carbon Films
• Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds
• Local strain engineering in atomically thin MoS2
• DOLFIN: Automated Finite Element Computing
• Parametrization of Stillinger-Weber Potential Based on Valence Force Field Model: Application to Single-Layer MoS2 and Black Phosphorus
• A compiler for variational forms
• Assessment and optimization of the fast inertial relaxation engine (FIRE) for energy minimization in atomistic simulations and its implementation in LAMMPS
• From Layers to Nanotubes: Transition Metal Disulfides TMS2
• Mechanical response of all-MoS2 single-layer hetrostructures: A ReaxFF investigation
• TSFC: a structure-preserving form compiler
• A new shell formulation for graphene structures based on existing ab-initio data
• Modal analysis of graphene-based structures for large deformations, contact and material nonlinearities
• A new efficient hyperelastic finite element model for graphene and its application to carbon nanotubes and nanocones
• Comparing quantum, molecular and continuum models for graphene at large deformations