Multi-Paradigm Modeling of Fracture of a Silicon Single Crystal under Mode II Shear Loading

Abstract

We report a novel multi-paradigm multi-scale approach based on a combination of the first principles ReaxFF force field with an empirical Tersoff potential. Our hybrid multi-scale simulation model is computationally efficient and capable of treating thousands of atoms with QM accuracy, extending our ability to simulate the dynamical behavior of a wider range of chemically complex materials such as silicon, silica and metal-organic compounds. It is implemented in the Python based Computational Materials Design Facility (CMDF). We exemplify our method in a study focused on a systematic comparison of the fracture dynamics in silicon under mode II shear versus mode I tensile loading. We find that the mode II crack tends to branch at an angle of approximately 45 degrees once the crack speed approaches 38% of the Rayleigh-wave speed. In contrast, the mode I crack continuously propagates in the direction of the initial crack, and only makes a slight change of direction towards 10 degrees once fracture instabilities occur. Our results reveal fundamental differences of fracture dynamics under mode I versus mode II loading.

Keywords

Mechanics chemistry elasticity fracture silicon mode I mode II

References

Freund

L.B.

, Dynamic Fracture Mechanics. 1990: Cambridge University Press, ISBN 0-521-30330-3.

Buehler

M.J.

and Gao

, Dynamical fracture instabilities due to local hyperelasticity at crack tips. Nature, 2006. 439: P. 307–310.

Buehler

M.J.

Abraham

F.F.

, and Gao

, Hyperelasticity governs dynamic fracture at a critical length scale. Nature, 2003. 426: P. 141–146.

Tadmor

E.B.

, and Kaxiras

, From electrons to finite elements: A concurrent multiscale approach for metals. Physical Review B, 2006. 73(2).

Choly

, Multiscale simulations in simple metals: A density-functional-based methodology. Physical Review B, 2005. 71(9).

Becke

A.D.

, Density-function thermochemistry. 3. The role of exact exchange. J. Chem. Phys., 1993. 98(7): P. 5648–5652.

Duin

A.C.T.V.

, ReaxFF: A Reactive Force Field for Hydrocarbons. J. Phys. Chem. A, 2001. 105: P. 9396–9409.

Duin

A.C.T.V.

, ReaxFF SiO: Reactive Force Field for Silicon and Silicon Oxide Systems. J. Phys. Chem. A, 2003. 107: P. 3803–3811.

Buehler

M.J.

Duin

A.C.T.V.

, and Goddard

W.A.

, Multi-paradigm modeling of dynamical crack propagation in silicon using the ReaxFF reactive force field. Phys. Rev. Lett., 2006. 96(9): P. 095505

10.

Deegan

R.D.

, Wavy and rough cracks in silicon. Phys. Rev. E, 2003. 67(6): P. 066209.

11.

Cramer

Wanner

, and Gumbsch

, Energy dissipation and path instabilities in dynamic fracture of silicon single crystals. Phys. Rev. Lett., 2000. 85: P. 788–791.

12.

Cramer

Wanner

, and Gumbsch

, Crack Velocities during Dynamic Fracture of Glass and Single Crystalline Silicon. Phys. Status Solidi A, 1997. 164: P. R5.

13.

Hauch

J.A.

, Dynamic fracture in Single Crystal Silicon. Phys. Rev. Lett., 1999. 82: P. 3823–2826.

14.

Holland

and Marder

, Ideal brittle fracture of silicon studied with molecular dynamics. Phys. Rev. Lett., 1998. 80(4): P. 746.

15.

Abraham

F.F.

, Spanning the length scales in dynamic simulation. Computers in Physics, 1998. 12(6): P. 538–546.

16.

Bailey

N.P.

and Sethna

J.P.

, Macroscopic measure of the cohesive length scale: Fracture of notched single-crystal silicon. Phys. Rev. B, 2003. 68(20): P. 205204.

17.

Bazant

M.Z.

Kaxiras

, and Justo

J.F.

, Environment-Dependent Interatomic Potential for bulk silicon. Physical Review B-Condensed Matter, 1997. 56: P. 8542.

18.

Swadener

J.G.

Baskes

M.I.

, and Nastasi

, Molecular Dynamics Simulation of Brittle Fracture in Silicon. Phys. Rev. Lett., 2002. 89(8): P. 085503.

19.

Tersoff

, Empirical interatomic potentials for carbon, with applications to amorphous carbon. Phys. Rev. Lett., 1988. 61(25): P. 2879–2883.

20.

Stillinger

and Weber

T.A.

, Computer-simulation of local order in condensed phases of silicon. Phys. Rev. B, 1985. 31(8): P. 5262–5271.

21.

Bernstein

and Hess

D.W.

, Lattice trapping barriers to brittle fracture. Physical Review Letters, 2003. 91(2): P. 025501.

22.

Broughton

J.Q.

, Concurrent coupling of length scales: Methodology and application. Physical Review B, 1999. 60(4): P. 2391–2403.

23.

Xiao

S.P.

and Belytschko

, A bridging domain method for coupling continua with molecular dynamics. Computer Methods in Applied Mechanics and Engineering, 2004. 193(17–20): P. 1645–1669.

24.

Wagner

G.J.

and Liu

W.K.

, Coupling of atomistic and continuum simulations using a bridging scale decomposition. Journal Of Computational Physics, 2003. 190(1): P. 249–274.

25.

Park

H.S.

, The bridging scale for two-dimensional atomistic/continumn coupling. Philosophical Magazine, 2005. 85(1): P. 79–113.

26.

Tadmor

E.B.

Ortiz

, and Phillips

, Quasicontinuum analysis of defects in solids. Phil. Mag. A, 1996. 73: P. 1529.

27.

Knap

and Ortiz

, An analysis of the quasicontinuum method. Journal Of The Mechanics And Physics Of Solids, 2001. 49(9): P. 1899–1923.

28.

Dupuy

L.M.

, Finite-temperature quasicontinuum: Molecular dynamics without all the atoms. Physical Review Letters, 2005. 95(6).

29.

Buehler

M.J.

, Threshold Crack Speed Controls Dynamical Fracture of Silicon Single Crystals. Phys. Rev. Lett., 99: P.165502.

30.

Buehler

M.J.

, The Computational Materials Design Facility (CMDF): A powerful framework for multiparadigm multi-scale simulations. Mat. Res. Soc. Proceedings, 2006. 894: P. LL3.8.

31.

Allen

M.P.

and Tildesley

D.J.

, Computer Simulation of Liquids. 1989: Oxford University Press.

32.

Gao

Huang

, and Abraham

F.F.

, Continuum and atomistic studies of intersonic crack propagation. J. Mech. Phys. Solids, 2001. 49: P. 2113–2132.