Open Access Nano Express

Atomistic aspects of ductile responses of cubic silicon carbide during nanometric cutting

Saurav Goel, Xichun Luo*, Robert L Reuben and Waleed Bin Rashid

Author affiliations

School of Engineering and Physical Sciences, Heriot Watt University, Edinburgh, EH144AS, Scotland, UK

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2011, 6:589  doi:10.1186/1556-276X-6-589


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/6/1/589


Received:2 August 2011
Accepted:11 November 2011
Published:11 November 2011

© 2011 Goel et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Cubic silicon carbide (SiC) is an extremely hard and brittle material having unique blend of material properties which makes it suitable candidate for microelectromechanical systems and nanoelectromechanical systems applications. Although, SiC can be machined in ductile regime at nanoscale through single-point diamond turning process, the root cause of the ductile response of SiC has not been understood yet which impedes significant exploitation of this ceramic material. In this paper, molecular dynamics simulation has been carried out to investigate the atomistic aspects of ductile response of SiC during nanometric cutting process. Simulation results show that cubic SiC undergoes sp3-sp2 order-disorder transition resulting in the formation of SiC-graphene-like substance with a growth rate dependent on the cutting conditions. The disorder transition of SiC causes the ductile response during its nanometric cutting operations. It was further found out that the continuous abrasive action between the diamond tool and SiC causes simultaneous sp3-sp2 order-disorder transition of diamond tool which results in graphitization of diamond and consequent tool wear.

Keywords:
ductile regime nanometric cutting; silicon carbide; diamond tool; tool wear

Introduction

Silicon carbide (SiC) is a promising ceramic material suited for advanced neural interfaces, packaging for long-term implantation, microfabricated neural probe as well as for semiconductor devices used in severe environments, such as in military aircraft, combat vehicles, power generation, and petrochemical industries [1]. SiC is a very hard substance (9 to 9.5 on Mohs scale) having comparable hardness to the hardest material known as diamond (10 on Mohs scale). The unique blend of properties possessed by SiC which makes it suitable for various MEMS, bio-medical, and other applications can be summarised in the form of Table 1[2].

Table 1. Commercial applications of SiC [2]

Moreover, SiC is also capable of meeting the requirements of operation in hostile environments (up to 873 K) where conventional silicon-based electronics (limited to 623 K) cannot function. The National Aeronautics and Space Administration agency, NASA, has recently been making efforts to develop SiC as future material for advanced semiconductor electronic device applications [3].

Single-point diamond turning (SPDT) is now an established ultra-precision machining process to manufacture free-form shapes and mirror-finished machined surfaces [4,5]. SPDT was established by exploiting a so called "brittle-ductile transition" phenomenon which has made various brittle materials, amenable to ultra-precision machining using a diamond cutting tool [6,7]. Investigations on exploring silicon carbide as a diamond-turnable material are thus of scientific and technological interest.

Experimental investigation on the feasibility of ductile regime machining of SiC through SPDT was first reported in 2005 [8]. Common believe on machining mechanism of SiC has been that a nanoscale undeformed chip thickness compounded with slow feed rate helps to achieve high-pressure phase transformations (HPPT) which causes ductile responses from this brittle material [9,10]. However, no such evidence of HPPT during nanometric cutting of SiC has been reported yet and the root cause of ductile response of SiC is still unknown. Molecular dynamics (MD) simulation results have been successful in the past to address number of problems concerning the nanometric cutting processes of brittle materials like silicon [11-14].

In this paper, Tersoff potential energy function [15] was used in the MD simulation to elucidate the atomistic mechanism underlying the ductile responses from the cubic SiC during nanometric cutting. Resorting to the simulation results, a theory has been presented and discussed.

MD simulation

MD simulation model

A schematic diagram of the nanometric cutting simulation model is shown in Figure 1. In the simulation model, the single-crystal diamond cutting tool has been modelled as a deformable body. Accordingly, both the work piece and the cutting tool are divided into three different zones: Newton atoms, thermostatic atoms, and boundary atoms. The boundary atoms are unaffected during the simulation and remain fixed in their initial lattice positions, serving to reduce the boundary effects and maintain the symmetry of the lattice.

thumbnailFigure 1. Schematic of MD simulation model.

The motions of the atoms in the Newton and thermostatic zones are assumed to follow the classical Newton's second law of motion which can be computed from the interatomic forces produced by the interaction as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M1">View MathML</a>

(1)

where aix represents the ith atom's acceleration in the x direction, mi is the mass of the ith atom, Fix is the interaction force acting on the ith atom by the jth atom in the x direction, xi indicates the ith atom's x-coordinate and V is the potential energy.

During realistic machining operation, the energy from plastic deformation in the primary shear zone and friction at the tool-chip interface transforms to cutting heat, which is carried away by chips and lubricants. The motion of the thermostatic atoms is therefore scaled using the velocity scaling factor (shown in equation 2) to incorporate the dissipation of heat in the simulation.

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M2">View MathML</a>

(2)

The temperature of atoms during the machining simulation can be calculated using the conversion between the kinetic energy (K.E.) and temperature as shown in Eq. 3.

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M3">View MathML</a>

(3)

where N is the number of atoms, vi represents the velocity of ith atom, kb is the Boltzmann constant which is equal to 1.3806503 × 10−23 J/K and T represents the temperature on atoms. However, the instantaneous fluctuations in K.E. of atoms could be very high so K.E. should be averaged (time and/or spatial) over few timesteps to be converted into equivalent temperature. It shall be noted here that the movement of the tool will also be a contributor to the kinetic energy so the tool displacement was accordingly subtracted before equivalent temperature conversion.

Selection of potential energy function

The interaction potential function governs the accuracy of a molecular dynamic simulation which in turn defines the reliability of simulation results. In this paper, Tersoff potential energy function [15] was used in the simulation for accurate description of the SiC mixture. Being a three-body potential function, Tersoff function is suitable to describe covalent interactions of Si and C atoms. Therefore, it was used to describe Si-Si, C-C, and Si-C interactions for interactions in and among the tool and workpiece as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M4">View MathML</a>

(4)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M5">View MathML</a>

(5)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M6">View MathML</a>

(6)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M7">View MathML</a>

(7)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M8">View MathML</a>

(8)

where Ei is the site energy-the sub-function, Vij describes the energy between two atoms (i and j), (i, j, and k) label the atoms of the system, fR represents a repulsive pair potential, fA represents an attractive pair potential, fC represents a smooth cut-off function to limit the range of the potential, rij is the length of the i-j bond, bij is the bond order term, ζij counts the number of other bonds to atom i besides the i-j bond and θijk is the bond angle between the bonds i-j and i-k. Here the indices ij represents the atom species. The mixing parameters between the two atomic species can be obtained from the following mixing rules:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M9">View MathML</a>

(9)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M10">View MathML</a>

(10)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M11">View MathML</a>

(11)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M12">View MathML</a>

(12)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M13">View MathML</a>

(13)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M14">View MathML</a>

(14)

Parameter χij determines the attractive interactions between two atoms. The potential function parameters used in the study have been listed in Table 2.

Table 2. Tersoff potential parameters [15]

Large-scale atomic/molecular massively parallel simulator software [16] was used to perform the simulation.

Calculation of equilibrium lattice parameter

Using inappropriate lattice parameter will affect the total energy content of the system which may lead to lots of thermal vibrations during equilibration process. Thus, the resulting fluctuations will alter the machining parameters like undeformed chip thickness, nose radius, etc. to a large extent during energy minimization which will produce erroneous simulation results. Goel et al. [17,18] have recently suggested to use the equilibrium lattice parameters to represent realistic MD simulation.

Accordingly, in the current work, equilibrium lattice parameters as shown in Table 3 were used as an input to the MD simulation in order to obtain meaningful simulation results [19,20].

Table 3. Comparison of lattice parameters obtained through experiment and simulation

MD simulation setup

MD simulation model of single-crystal SiC and diamond was built using the periodic boundary conditions along the z direction. This was followed by an energy minimization to avoid overlaps in the positions of the atoms. The simulation model was equilibrated to 300 K under the microcanonical (NVE) ensemble and the initial velocities of the atoms were assigned in accordance with a Maxwell-Boltzmann distribution. During the equilibration process, the trajectories must not be used to compute any properties as the potential energy continues to convert to kinetic energy or vice-versa. This causes the temperature to fluctuate until it becomes stationary. Once sufficient time has been given for equilibration, velocity scaling is removed and the system then follows NVE dynamics. Different variables which were used in the current simulation model have been listed in Table 4.

Table 4. Variables used in the MD simulation model

A model output from the MD simulation after relaxation has been shown in Figure 2 where, ice blue coloura and yellow colour correspond to silicon atoms and carbon atoms in the cubic SiC, respectively. The ochre colour represents carbon atoms in the diamond tool. Visual Molecular Dynamics (VMD) [21]b along with Ovito [22]c were used for the enhanced visualisation of the atomistic data.

thumbnailFigure 2. Orthographic view of model after relaxation.

As shown in Figure 3, SiC has a tetrahedral molecular geometry which has a central carbon atom surrounded at equidistance by four silicon(Si) atoms or vice-versa and thus bond angles are cos−1(−1/3) ≈ 109.5° between them. A tetrahedral geometry can however, be distorted by increasing the bond angles. This distortion of lattice structure and consequent increase in the bond angle is accompanied by the changes in the bond length. This phenomenon is called planarization which causes flattening of tetrahedron structure due to increment in bond angle. During MD simulation, a similar phenomenon was observed which has been discussed further in the subsequent sections.

thumbnailFigure 3. Lattice structure of cubic SiC.

Results and discussions

Observations from MD simulation results

Figure 4 shows the snapshot from the MD simulation after 10 nm of tool advance. The work piece shown in Figure 4 has Si-terminated surfaces having silicon atoms at the edge of the workpiece. The nature of SiC is extremely brittle and hence surface roughness on the machined surface of SiC in Figure 4 appears to be relatively higher as has been observed during the experimental work [23].

thumbnailFigure 4. Chip morphology of β-SiC (cubic) during the nanometric cutting process.

Figure 4 also represents the chip morphology of β-silicon carbide (cubic) during the nanometric cutting process against a deformable diamond tool. It can be seen from Figure 4 that the cutting chips are curly shaped, which suggests that material removal is occurring in ductile regime by deformation rather than fracture. During the machining process, a few carbon atoms from the diamond tool deformed and are separated from the tool which can be seen on the machined surface of SiC.

Temperature during the machining process

Figure 5 represents the temperature distribution on the workpiece and the cutting tool during nanometric cutting after the cutting tool has advanced by 4 nm. It can be seen that the maximum temperature in the primary shear zone of the workpiece approaches to a value of 1,700 K. The high temperature in the primary shear zone on SiC during its nanometric cutting is closer in magnitude to the enabling temperature(1,400°C) for the formation of SiC-graphene on either Si- or C-terminated SiC surfaces as reported [24]. Accordingly, it is plausible to state that high temperature ignites the formation of SiC-graphene-like substance during the nanometric cutting operation of SiC. SiC-graphene being much weaker and slippery compared to cubic SiC causes ductile responses in the cutting zone. The formation of SiC graphene from cubic SiC is essentially an outcome of sp3-sp2 order disorder transition. In order to ensure the presence of a substance like SiC-graphene during cutting, angular distribution functions, and radial distribution functions of SiC before and during the cutting were plotted which have been elaborated in subsequent sections.

thumbnailFigure 5. Temperature distribution on the workpiece and the cutting tool after 4 nm of cut.

Angular distribution function of SiC during cutting

Figure 6 represents angular distribution function of SiC workpiece obtained from the simulation before and during the cutting. It is evident from Figure 6 that before cutting, the peak of angular distribution function of SiC is visible at a bond angle of 109.5° which confirm the perfect sp3 tetrahedral bonding structure of SiC. During cutting, the bond angle of around 5% of atoms reduces at the bond angle of 109.5° with a corresponding increase and a small peak at a value of 120°. A change in value of bond angle from 109.5° to 120° is a strong indication of sp3-sp2 hybridisation in the manner shown in Figure 7.

thumbnailFigure 6. Angular distribution function of SiC before and during cutting.

thumbnailFigure 7. sp3-sp2 hybridisation.

As shown in Figure 7, the sp3 hybridisation has a tetrahedral geometry with a bond angle of 109.5°. However, when this tetrahedral geometry is distorted the bond angle changes to 120° and results in planarization and consequent sp2 hybridisation. This sp3-sp2 transition occurs due to the abrasive action between diamond cutting tool and SiC workpiece. Thus, a change in bond angle from 109.5° to 120° obtained through angular distribution function is an indication of sp3-sp2 order-disorder transition and transformation of cubic silicon carbide to SiC-graphene-like substance.

Radial distribution function and relative tool wear

The radial distribution functions (RDF), g(r), also called pair distribution functions or pair correlation functions, are the primary linkage between macroscopic thermodynamic properties and intermolecular interactions.

As illustrated in Figure 8[25], if the atoms are distributed homogeneously in space, then the RDF, g(r), gives the probability of finding an atom in a shell dr at a distance r from another atom chosen as a reference point. The number of atoms dn(r) at a distance between r and r + dr from a given atom is expressed as follows:

thumbnailFigure 8. Space discretization for the evaluation of the radial distribution function [25] .

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/589/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/589/mathml/M15">View MathML</a>

(15)

where N represents the total number of atoms, V is the model volume and g(r) is the radial distribution function. Thus, g(r) can be employed to observe the variation in magnitude of bond length of atoms during nanometric cutting. Accordingly, radial distribution functions of SiC workpiece and diamond tool were plotted before and during the cutting process which has been shown in Figure 9.

thumbnailFigure 9. Radial distribution function of SiC workpiece.

It is evident from Figure 9 that before cutting, the first peak of RDF in SiC workpiece was visible at a bond length of 1.85 Å which is the equilibrium bond length of cubic SiC. A second small peak at a bond length of 1.9 Å was also evident which represents Si-terminated dangling bonds on the surface. During nanometric cutting, the number of atoms at a bond length of 1.85 Å decreases with corresponding major increase at a bond length of 1.9 Å while few others at a bond length of 1.75 Å. These two bond lengths of 1.9 and 1.75 Å are in excellent correlation with the reported bond length 1.87 Å [26] and 1.77 Å [27] of Si-terminated graphene and Si = C double bond length respectively. Thus, RDF reconfirms the presence of SiC-graphene-like substance during the nanometric cutting process. The high magnitude of abrasion between SiC and diamond tool also results in the sp3-sp2 order-disorder transition and consequent graphitization of diamond tool but at relatively slower rate. The RDF of the diamond tool before and after the cutting is shown in Figure 10.

thumbnailFigure 10. Radial distribution function of diamond tool.

It is widely known that the thermal stability of diamond gets adversely effected in the environment of severe temperature [28]. This causes graphitization/sp3-sp2 order-disorder transition of the diamond tool during nanometric cutting of SiC. It can be seen from Figure 10 that at timestep 0, the RDF of diamond tool shows its first peak at 1.54 Å which is the known bond length of diamond [29] while few bonds(dangling bonds) on the surface shows a small peak at 1.42 Å. During cutting, the small peak continued to grow at a bond length of 1.42 Å with corresponding decrease in the number of atoms at the bond length of 1.54 Å. The bond length of 1.42 Å is the known bond length of another stable allotrope of carbon known as graphite [29] which is much weaker than diamond due to the layered structure. Thus, g(r) confirms the graphitization of the diamond tool during SPDT operation of cubic SiC and consequent wear as earlier observed from Figure 4. The numbers of the SiC atoms in the cutting chips are more than those deformed carbon atoms from the diamond tool and this proves that the rate of sp3-sp2 transition of diamond tool is relatively slower than SiC.

It thus becomes plausible to state that the abrasion between SiC and the diamond tool causes the temperature rise in the cutting zone and consequent sp3-sp2 order-disorder transition of both these materials. The disorder taking place inside the diamond tool has been pictorially showed in Figure 11 where silver colour atoms represent sp3; green colour represents sp2 while yellow colours represent sp arrangement. The regular arrangement of sp2 and sp bonds at the tool cutting edges represents termination of bonds (dangling bonds) on the surface and sides of the cutting tool. However, the sp3-sp2-sp bond transition around the tool nose appears as an initiation point of graphitization of diamond tool.

thumbnailFigure 11. Wear of diamond tool. (a) sp3-sp2 transition during SPDT of SiC, (b) sp3-sp2 transition during polishing of diamond [30]

The phenomenon of sp3-sp2 order-disorder transition during SPDT of SiC also appears to be similar in nature to what occurs during polishing of diamond, which has been explained in details by Pastweka et al. [30] using MD simulation studies.

Conclusions

The MD simulation has been used to gain extensive insights into the atomistic aspects of ductile responses of SiC during nanometric cutting operations. The following conclusions can be drawn accordingly:

During nanometric cutting, the tetrahedral bonding structure of β-SiC work material gets distorted accompanying the change of bond angle from 109.5° to 120° which represents sp3-sp2 order-disorder transition of SiC.

This sp3-sp2 disorder causes the formation of SiC-graphene-like substance which causes ductile response from cubic SiC.

The formation of SiC-graphene-like substance is attributed to the high temperature during the nanometric cutting which is consequent due to the abrasive action between these two ultra-hard materials.

Abrasive action also causes simultaneous sp3-sp2 order-disorder transition of the diamond tool but at relatively slower rate which results in tool wear.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SG did the literature review and carried out the MD simulation. WBR provided some basic inputs to the MD simulation. XL suggested the basic framework of the manuscript and also did the proof reading of the manuscript. RLR helped with the critical corrections wherever necessary. All authors read and approved the final manuscript.

Endnotes

a Readers are requested to refer to the web-based version of this article for correct interpretation of the colour legends.

b Developed at the University of Illinois, USA

c Developed at LLNL, USA

References

  1. Yuan X, Hobbs LW: Influence of interatomic potentials in MD investigation of ordering in a-SiC in mat.

    Mater Res Soc Symp Proc 2001, 650:R3.18.1-R3.18.6. OpenURL

  2. Zorman CA: Silicon carbide as a material for biomedial microsystems. [http://www.youtube.com/watch?v=XodCl3qiiLg] webcite

  3. Silicon Carbide electronics [http://www.grc.nasa.gov/WWW/SiC/] webcite

  4. Wilks J: Performance of diamonds as cutting tools for precision machining.

    Precis Eng 1980, 2(2):57-70. Publisher Full Text OpenURL

  5. Shore P: Machining of Optical Surfaces in Brittle Materials using an Ultra Precision Machine Tool. In PhD Thesis. Cranfield University; 1995. OpenURL

  6. Lawn BR, Marshall DB: Hardness, toughness, and brittleness: an indentation analysis.

    J Am Ceram Soc 1979, 62(7-8):347-350. Publisher Full Text OpenURL

  7. Scattergood RO, Blake N: Ductile-regime machining of germanium and silicon.

    J Am Ceram Soc 1990, 73(4):949-957. Publisher Full Text OpenURL

  8. Patten J, Gao W, Yasuto K: Ductile regime nanomachining of single-crystal silicon carbide.

    J Manuf Sci Eng 2005, 127(3):522-532. Publisher Full Text OpenURL

  9. Patten J, Jacob J: Comparison between numerical simulations and experiments for single-point diamond turning of single-crystal silicon carbide.

    J Manuf Process 2008, 10:28-33. Publisher Full Text OpenURL

  10. Patten JA: High pressure phase transformation analysis and molecular dynamics simulations of single point diamond turning of germanium. In PhD Dissertation. North Carolina State University; 1996. OpenURL

  11. Belak JF, Stowers IF, Boercker DB: Simulation of diamond turning of silicon surfaces.

    Proceedings of 7th American Society Precision Engineering Annual conference 1992, 76-79. OpenURL

  12. Komanduri R, Chandrasekaran N, Raff LM: Effect of tool geometry in nanometric cutting: a molecular dynamics simulation approach.

    Wear 1998, 219(1):84-97. Publisher Full Text OpenURL

  13. Cheng K, Luo X, Holt R, Ward R: Modeling and simulation of the tool wear in nanometric cutting.

    Wear 2003, 255:1427-1432. Publisher Full Text OpenURL

  14. Goel S, Luo X, Reuben RL, Pen H: Wear mechanism of diamond tools against single crystal silicon during the single point diamond turning process.

    Diam Relat Mater 2011, in press. OpenURL

  15. Tersoff J: Chemical order in amorphous silicon carbide.

    Physical Review B 1994, 49(23):16349. Publisher Full Text OpenURL

  16. Plimpton S: Fast parallel algorithms for short-range molecular dynamics.

    J Comput Phys 1995, 117:1-19. Publisher Full Text OpenURL

  17. Goel S, Luo X, Reuben RL: Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide.

    Comput Mater Sci 2012, 51(1):402-408. Publisher Full Text OpenURL

  18. Goel S: Wear Mechanism of Diamond Tools during Ultra Precision Machining: MD Simulation Study, Volume 1.

    Germany: LAP LAMBERT Academic Publishing 2011, 100. OpenURL

  19. Lattice of SiC [http://www.ioffe.rssi.ru/SVA/NSM/Semicond/SiC/basic.html] webcite

  20. Lattice constants [http://www.siliconfareast.com/lattice_constants.htm] webcite

  21. Humphrey W, Dalke A, Schulten K: VMD - visual molecular dynamics.

    J Molec Graph 1996, 14:33-38. Publisher Full Text OpenURL

  22. Stukowski A: Visualization and analysis of atomistic simulation data with OVITO-the Open Visualization Tool.

    Model Simulat Mater Sci Eng 2010, 18(1):015012. Publisher Full Text OpenURL

  23. Yan J, Zhang Z, Kuriyagawa T: Mechanism for material removal in diamond turning of reaction-bonded silicon carbide.

    Int J Mach Tool Manufac 2009, 49(5):366-374. Publisher Full Text OpenURL

  24. Mattausch A, Pankratov O: Ab Initio Study of Graphene on SiC.

    Phys Rev Lett 2007, 99(7):076802. PubMed Abstract | Publisher Full Text OpenURL

  25. g(r) RDF explanation [http://isaacs.sourceforge.net/phys/rdfs.html] webcite

  26. Jahn C: Density functional study of carbonmono- and bilayers on SiC. [http://www14.informatik.tu-muenchen.de] webcite

  27. Huda MN, Yan Y, Al-Jassim MM: On the existence of Si-C double bonded graphene-like layers.

    Chem Phys Lett 2009, 479(4-6):255-258. Publisher Full Text OpenURL

  28. Cai MB, Li XP, Rahman M: Study of the temperature and stress in nanoscale ductile mode cutting of silicon using molecular dynamics simulation.

    J Mater Process Tech 2007, 192-193:607-612. OpenURL

  29. Bond length [http://hypertextbook.com/facts/2001/AliceWarrenGregory.shtml] webcite

  30. Pastewka L, Moser S, Gumbsch PMoseler M: Anisotropic mechanical amorphization drives wear in diamond.

    Nat Mat 2011, 10(1):34-38. Publisher Full Text OpenURL