Next Article in Journal
Heat Transfer Analysis of Nanostructured Material Flow over an Exponentially Stretching Surface: A Comparative Study
Next Article in Special Issue
Adsorption Free Energy of Cellulose Nanocrystal on Water–Oil Interface
Previous Article in Journal
Graphene Oxide-Coated Metal–Insulator–Metal SERS Substrates for Trace Melamine Detection
Previous Article in Special Issue
Molecular Dynamics Simulation Study of the Self-Assembly of Phenylalanine Peptide Nanotubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Revisiting the Rate-Dependent Mechanical Response of Typical Silicon Structures via Molecular Dynamics

1
Institute of Photovoltaics, Nanchang University, Nanchang 330031, China
2
Institute for Advanced Study, Nanchang University, Nanchang 330031, China
3
School of Automotive Application, Hunan Automotive Engineering Vocational College, Zhuzhou 412001, China
4
School of Material Science and Engineering, Nanchang University, Nanchang 330031, China
5
Department of Science and Technology, Nanchang University College of Science and Technology, Jiujiang 332020, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Nanomaterials 2022, 12(7), 1203; https://doi.org/10.3390/nano12071203
Submission received: 11 February 2022 / Revised: 31 March 2022 / Accepted: 1 April 2022 / Published: 3 April 2022
(This article belongs to the Special Issue Nanomaterials Investigation by Molecular Dynamics Simulation)

Abstract

:
Strain rate is a critical parameter in the mechanical application of nano-devices. A comparative atomistic study on both perfect monocrystalline silicon crystal and silicon nanowire was performed to investigate how the strain rate affects the mechanical response of these silicon structures. Using a rate response model, the strain rate sensitivity and the critical strain rate of two structures were given. The rate-dependent dislocation activities in the fracture process were also discussed, from which the dislocation nucleation and motion were found to play an important role in the low strain rate deformations. Finally, through the comparison of five equivalent stresses, the von Mises stress was verified as a robust yield criterion of the two silicon structures under the strain rate effects.

1. Introduction

Monocrystalline silicon (Mono-Si) has attracted great attention in recent decades due to its material properties and promising applications in the field of semiconductors, photovoltaics and Micro Electro Mechanical Systems [1,2]. Among these attractive properties, mechanical performance is crucial for future applications of silicon-based devices because it is a prerequisite to fulfill other functionalities. Although perfect Mono-Si crystals can be produced through Floating Zone [3] or Czochralski [4] methods, traditional mechanical experiments are still difficult to execute due to the inherent brittleness of silicon. Moreover, the produced Mono-Si wafers have also exhibited edge collapse, hidden cracks and subsurface damage [5,6] in wiresaw cutting, which have slowed its commercial expansion in the photovoltaic industry. Therefore, learning the mechanical properties of silicon is still highly desirable, more and more efforts have been continuously carried out to solve these problems.
Among the relevant studies about the mechanical properties of silicon, a particularly intriguing question is how the applied strain rate affects its mechanical performance, fracture behaviors and deformation mechanisms. The corresponding solutions will not only deepen our understandings of damage tolerance and deformation behaviors [7,8,9], but may also create new opportunities and challenges for silicon nanodevices [10,11] and silicon-based anodes [12,13]. Usually, three widely known methods are used to explore the mechanical response of silicon under strain rate effects. The first comprises theoretical efforts to conclude and predicate the rate-dependent mechanical response, such as constitutive equations [14,15,16], dislocation plasticity models [17,18] and rate response models [19]. However, the material properties of silicon, such as its widely known brittle-to-ductile transition characteristics, often make theoretical study complicated. The second is low strain rate nanomechanical experiments [20,21,22] (mostly <100/s). For example, in situ tensile experiments on amorphous silicon nanowires indicated that brittle-to-ductile transitions were strongly rate-dependent [23]. Smith et al. [24] studied the compression behaviors of monocrystalline silicon, from which the rate- and orientation-dependency were discussed. The rate-dependence of hardness [22] in amorphous silicon and the effects of scribing rate on monocrystalline silicon [25] have been reported using nanoindentation tests. The promising engineering prospects given by these experimental methods have stimulated the efforts to understand the strain rate-dependency of silicon properties. However, for some unavailable experimental conditions, such as defect-free and extremely high strain rates, researchers have often tried to resort the third method: high strain rate molecular dynamics simulations (mostly >108/s). Recent simulations [26,27,28,29,30,31] have shown that the amount of plasticity remains constant for brittle glassy nanowires as the strain rate decreases, but that plasticity decreases for the ductile counterparts [26]. Chen et al. [28] suggested that the tensile strength and the fracture pattern of ideal silicon nanorod were strongly rate-dependent. The rate-dependent tensile response obtained through different potentials also exhibited great consistency, namely, the tensile strength would reduce as applied strain rate decrease [28,29].
However, a huge gap in the timescale between simulations and experiments [32] always results in discrepancies in learning the rate-dependent mechanical response, especially when the deformation mechanisms show rate-dependency. To overcome this shortage of molecular dynamics methods, simulations covering timescales over ten orders of magnitude should be considered [33]. Although there are abundant literatures about the rate-dependent mechanical response of silicon, most of these studies have focused on amorphous silicon structures or surface indentation, resulting in a relatively poor understanding of some typical silicon structures, such as the monocrystalline silicon wafer and silicon nanowire [34,35]. Meanwhile, for the yield criterion of silicon, the effects of strain rate have not been systematically discussed through the aspects of equivalent stresses. The discovery of the yield criteria of typical silicon structures (single crystal, nanowire, nanotube, etc.) under strain rate effects would certainly benefit their engineering applications. Starting from these concerns, we are trying to figure out some unsolved questions regarding the rate-dependent mechanical response of silicon on the basis of present researches, especially concentrated on ignored points, such as surface effects and yield criterion.
With the aid of molecular dynamics, the rate-dependent mechanical response of two typical monocrystalline silicon structures were investigated in the present paper. Two important parameters about the rate-dependent mechanical response were calculated according to a rate response model. Meanwhile, the hidden rate-dependent deformation mechanism and yield criterion were discussed through analyzing dislocation activities and comparing equivalent stresses, respectively.

2. Methodologies

2.1. Model and Method

All simulations were carried out by LAMMPS software (Large-scale Atomic/Molecular Massive Parallel Simulator, version 5 September 2018). First, a monocrystalline silicon crystal and a cylindrical silicon nanowire were generated, their schematics are shown in Figure 1. The X, Y and Z axes of the two structures corresponded to [100], [010] and [001] orientations of silicon, respectively. Periodic conditions were applied in all directions of systematic boundaries to eliminate additional surface effects.
After the establishment of these two structures, they were relaxed under a 300 K, 1 atm standard NPT (isobaric/isothermal constant number of particles, constant pressure and constant temperature) ensemble. The pressures of the X, Y and Z axes were controlled close to zero by a Berendsen [36] barostat in the relax stage. However, in the deformation process, only the pressures of the Y and Z axes were controlled close to 100 MPa by the Berendsen barostat, in order to ensure the ambient pressure. The system temperature was controlled close to 300 K via a Berendsen thermostat during the deformation process. The selection of interatomic potential is critical for the simulation results. To obtain accurate descriptions about the elastic properties of silicon, we adopted Tersoff’s approach that developed in 1988 [37]. The comparison of interatomic potentials about material properties was given in Table A1 (see Appendix A for details).
Then, for each simulation case, a uniaxial tensile strain rate ε ˙ was loaded along the X axis to commence the deformation process. Ten strain rates were simulated for each structure, which ranged from 10−1/ps to 5 × 10−6/ps (e.g., the first was 10−1/ps, the second was 5 × 10−2/ps, the third was 10−3/ps, etc.) to cover the possible range of molecular dynamics simulations with calculation consumption as low as possible. In addition, we used the logarithm to express the strain rate, which was intended to ensure clarity and readability.

2.2. Relevant Theories

In this study, a modified rate response model [19] based on FCC metals was introduced to examine the rate-dependence of silicon. For applied strain rate ε ˙ , the modified model takes the form below:
σ s = σ q u a s i - s t a t i c 1 + c × ε ˙ ε ˙ c p ,
where σs is the rate-dependent tensile strength and σquasi-static is the quasi-static value (quasi-static strain rate range: 10−17/ps to 10−14/ps [38]) of tensile strength. p represents strain rate sensitivity (for Cu, p ≈ 0.64 and for Ni, p ≈ 0.5 [19]). ε ˙ c is the critical strain rate. The tensile strength changes a little at low strain rate and increases rapidly after the critical strain rate. Under this condition, however, ε ˙ c cannot be treated as the critical strain rate, because the tensile strength is no longer rate-independent. The original form of Equation (1) was first proposed by Cowper and Symonds [39], but Guo et al. [19] modified the original form by adding parameter c, allowing the modified model to be adjusted according to ε ˙ c . According to their article, the modified parameter c is defined by Equation (2):
c = σ c σ q u a s i - s t a t i c 1 ,
where σc is the critical tensile strength under the corresponding ε ˙ c . Therefore, the modified rate response model is introduced, and it will be discussed in the next section.

3. Results

3.1. Strain Rate Sensitivity

The tensile stress–strain curves of Mono-Si and SiNW are plotted in Figure 2a,a1,b,b1, from which we can see a smooth stress increase and then a sudden drop. Both the tensile strength and the maximum tensile strain reduce, while Young’s modulus stays roughly at 166.37 GPa for Mono-Si and 20.79 GPa for SiNW as the applied strain rate decreases from 10−1/ps to 5 × 10−6/ps. We also found that the toughness of Mono-Si is stronger than SiNW by comparing the maximum tensile strain under the same strain rate. In addition, it should be noted in Figure 2 that both Mono-Si and SiNW match the expectations of rate-dependent materials [40] under the test strain rate, according to the stress–strain curves and the variations of tensile strength. So, the theoretical rate response model introduced in Section 2.2 could be applied after obtaining ε ˙ c and σquasi-static.
First, we plot lg ε ˙ versus σs in Figure 3a to compute the quasi-static tensile strength. The exponential relation is used to fit lg ε ˙ versus σs for both structures, and it takes the form below:
σ s = σ 0 + A × e R × lg ε ˙ ,
where A and R are automatically fitting parameters and σ0 represents the tensile strength under the ultra-low strain rate ( ε ˙     10−6/ps) in this fitting. In this study, σ0 is approximately treated as σquasi-static based on the results given in Figure 3a. Then, we can rewrite Equation (1) as:
lg σ s σ q u a s i - s t a t i c 1 = p lg ε ˙ + lg c p lg ε ˙ c ,
where σquasi-static is equal to 17.335 GPa for Mono-Si and 2.43 GPa for SiNW. To figure out the strain rate sensitivity p, lg ε ˙ versus lg ( σ s σ q u a s i - s t a t i c 1 ) is plotted in Figure 3b, and the corresponding fittings show two linear relations. So, the rate sensitivity p is solved by calculating the slopes of linear fittings. The strain rate sensitivities of Mono-Si and SiNW are 0.267 and 0.226, respectively. However, the constant c and critical strain rate ε ˙ c in the rate response model are still unknown. The calculation of ε ˙ c is primarily considered using the structural variations. Figure 3c,c1 show the variations of cubic diamond and non-diamond atoms after fracture, which represent the undamaged structures and fractured structures, respectively. It is obvious that these structural variations follow the same trend with σs and the values of non-diamond percent and cubic diamond percent became stable under 10−4/ps. Combined with the analysis about the variations of σs, ε ˙ c could be roughly regarded as 10−4/ps. In that case, constant c is 0.0274 (Mono-Si) and 0.0493 (SiNW) through Equation (2). After obtaining all required parameters in the rate response model, we can conclude that ε ˙ c is the same, while p is different for the two silicon structures, which means the strain rate sensitivity is weakened by the additional surface effects of SiNW. In addition, the calculations of strain rate sensitivity and critical strain rate are also beneficial for predicting the rate-dependent mechanical response of these silicon structures.

3.2. Rate-Dependent Dislocation Activities

The fracture patterns of two structure types are plotted in Figure 4 via Ovito [41], from which the high strain rate deformation was found activating more crystal structures in the fracture process, resulting in the increase in disordered structures after fracture. On the other hand, the fractured structures gradually decrease to a low percent once the applied strain rate is lower than the critical strain rate. If we compare the fracture patterns of Mono-Si and SiNW, it is clear that the slip fracture tends to occur in SiNW while the cleavage fracture is more likely to appear in Mono-Si. This result shows great consistency with the existing literature [42,43].
In order to indicate the formations of these fracture patterns and the low tensile strength under low strain rate deformation, an automated dislocation extraction (DXA) algorithm [44] was used to identify the dislocations in each simulation and plot the results in Figure 5 and Figure 6. As shown in Figure 5a,a1,c,c1, intensive small defect clusters were generated in the high strain rate deformation of Mono-Si, while only a dislocation loop and some relatively larger defect clusters were generated in the low strain rate deformation. Compared with the Mono-Si counterparts, a lot of short dislocation lines were generated on the surface of SiNW in Figure 6a,a1, but these dislocations lacked the ability for expansion or motion. On the other hand, in Figure 5b,c and Figure 6b,c, a lower strain rate allowed the dislocation-driven slip fracture to expand along the {100} cleavage planes. Therefore, for Mono-Si, it is concluded that a low strain rate may benefit dislocation activities such as nucleations and motions, resulting in the low tensile strength of the tested Mono-Si crystals. As for the SiNW, the increase in strain rate would certainly reduce its dislocation motion abilities but increase the dislocation nucleation events in its surface. Besides, low coordination and high interface energy on the SiNW surface could also benefit the dislocation nucleation in the earlier stages of stress–strain curves, indicating that the introduction of additional surfaces causes the difference of tensile strength between Mono-Si and SiNW.

3.3. Comparison of Equivalent Stresses

To explore the hidden mechanisms of the rate-dependent mechanical response, the equivalent stresses were considered according to reference [42]. However, the stress components were required in the calculation of equivalent stresses. Here, a single atom was treated as a microelement so that its stress components could be represented by Virial stress [45]. First, Virial stress tensors [46,47] were computed for each atom. The computation details of these atomic Virial stress tensors are given in Appendix B for reference. Then, five equivalent stresses were calculated using the atomic Virial stress tensors following:
σ h = ( σ x + σ y + σ z ) 3 ,
σ r = σ 1 ,
σ s = σ 1 v ( σ 2 + σ 3 ) ,
σ t = σ 1 σ 3 2 ,
σ m = 1 2 ( σ x σ y ) 2 + ( σ y σ z ) 2 + ( σ z σ x ) 2 + 6 τ x y 2 + τ y z 2 + τ z x 2 ,
where σx, σy, σz, τxy, τyz and τxz are the Virial stress tensors of each atom. σ1, σ2 and σ3 are the three principal stresses computed from the atomic Virial stress tensors. v is the Poisson’s ratio. σh represents the hydrostatic stress and σr represents the maximum principal stress given by Rankine (Rankine stress). σs represents the equivalent stress based on the maximum principal strain given by St. Venant (St. Venant stress). σt represents the Tresca stress while σm represents the von Mises stress. The calculations of three principal stresses and five equivalent stresses follow the methods stated in reference [48]. Thus, five equivalent stresses were computed for each atom. Since fracture is the statistical results of atomic stress conditions, the average values of these equivalent stresses are given in Figure 7, and the maximum values are also given for reference.
As shown in Figure 7, the five equivalent stresses exhibited nearly the same variations for both Mono-Si and SiNW, so they jointly shared the same internal stress conditions in the deformation process. After confirming the consistency of the two structures in the stress conditions, the corresponding yield criteria were verified by comparing the five equivalent stresses in Figure 7. It is clear that the Rankine stress, St. Venant stress and hydrostatic stress were inappropriate for describing the yield criteria of Mono-Si and SiNW under strain rate effects due to their rate-independent variations. If we compare the Tresca stress and the von Mises stress, it is obvious that the von Mises stress showed better relevance in rate-dependency than the Tresca stress, because it gradually decreased as the applied strain rate increased. The von Mises stress not only considers the combination effects of tensile stress and shear stress, but also describes the stress conditions better than Tresca stress under the critical strain rate. Thus, through the comparison of five equivalent stresses, the von Mises stress was demonstrated to be the best in describing the yield criterion under various strain rates.
Another conclusion that can be carried out from Figure 7a,b is that the average von Mises stress was almost linear with lg ε ˙ when the applied strain rate exceeded ε ˙ c . Then, the average von Mises stress gradually started to approach a constant value under ε ˙ c , from which ε ˙ c acted as the division of such variations. The constant von Mises stress value could be considered as the stress threshold of fracture. Irreversible fracture is activated once the average von Mises stress reaches the threshold. In this way, the equivalent stresses and the corresponding yield criteria of the two silicon structures are analyzed and discussed under the strain rate effects. However, it is still unknown how the strain rate affects von Mises stress in the deformation process. To investigate this incomplete point, the distributions of von Mises stress are given in Figure 8 (Mono-Si) and Figure 9 (SiNW) in the form of probability density.
Obviously, the numerical distributions of von Mises stress followed the normal distribution. The stress peaks gradually shifted to a high value as the applied strain varied. Interestingly, the stress peaks exhibited different shapes under different strain rates. It seems that the atomic Virial stress was separated to a wide range of distribution interval by the high applied strain rate. Meanwhile, the width of the distribution peak (e.g., full width at half maximum) was narrowed by the low strain rate deformation. Once the stress distributions of three designated strain rates were compared at a same applied strain in Figure 8b and Figure 9b, it was found that the peak positions (average stress) were different. Accordingly, there were two parameters controlling the average von Mises stress in the deformation process. The first was the applied strain, which affected the peak positions of stress distributions. The second was the applied strain rate, which affected the peak shapes of stress distributions.
Different from the other rate-independent equivalent stresses, the average von Mises stress decreases as the applied strain rate increases. Combined with the colored subfigures about the crystal stress in Figure 8, the low strain rate deformation was found suffering a more serious von Mises stress concentration than the high strain rate deformation. Thus, the rate-dependent distributions of von Mises stress and the variations of average von Mises stress indicated an unexpected fact: compared with the high strain rate deformation, the stress fields in the low strain rate deformation had a higher possibility of exploring the potential stress threshold of fracture, or experiencing more frequent barrier crossing events [12] under the same applied strain, which could also be demonstrated by the variations of maximum von Mises stress shown in Figure 7a1,b1. Therefore, the fracture process will appear earlier in the stress–strain curves of low strain rate deformation because of the strengthened von Mises stress concentration.

4. Conclusions

In the present paper, two silicon structures were tested using high tensile strain rate molecular dynamics simulations. The results are concluded below:
Through a rate response model, the strain rate sensitivity and the critical strain rate of Mono-Si and SiNW were calculated, from which the extra surface of SiNW was found decrease strain rate sensitivity. However, the critical strain rate was the same for these structures, indicating that it is an inherent property which cannot be affected by surface effects. Then, the dislocation activities in the fracture process were examined. It was found that the rate-dependent dislocation nucleations and motions caused the earlier yield behaviors in the stress–strain curves of low strain rate deformation. Finally, five equivalent stresses and their descriptions on the rate-dependent mechanical response were verified; the von Mises stress was proved better than the Tresca stress in describing the yield criteria of the two silicon structures under the strain rate effects. The applicability of von Mises stress in the rate-dependent mechanical response indicates that a high strain rate not only affects internal stress fields, but also influences the elastic strain energy.

Author Contributions

Conceptualization, W.W.; methodology, Y.L., W.W. and C.T.; software, W.W., Z.X.; validation, Y.L., W.W., Q.L. and C.T.; formal analysis, W.W.; investigation, Y.L. and W.W.; resources, C.T.; data curation, Y.L. and W.W.; writing—original draft preparation, W.W.; writing—review and editing, Y.L., W.W. and C.T.; visualization, W.W.; supervision, C.T. and L.Z.; project administration, L.Z.; funding acquisition, C.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 11802112.

Data Availability Statement

Data available on request due to restrictions like data capacity. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to data capacity (According to an uncompleted statistic, the related simulation data in the present paper is about at least 50 GB) which is too big to storage for the website. Please inform the corresponding author for further usage of the related simulation data.

Acknowledgments

Thanks Changxin Tang and Lang Zhou for their supports in the preparation of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

To demonstrate the reason for using Tersoff T2 potential in this study, the descriptions of various methods about the material properties were compared. These listed material properties, such as elastic constants, vacancy formation energy and bulk modulus were calculated under equilibrium state and nearly zero strain. Only the elastic properties were considered in this article due to the research topic. In Table A1, the Tersoff T2 potential was found in well accordance with the experimental results of the elastic properties. The comparisons about other properties, such as thermal properties and crystal defects, were given in the reference [49].
Table A1. Properties of silicon in its equilibrium structure [49] from experiment [50], quantum-mechanical (QM) methods such as density functional theory (DFT) [51], as well as tight binding (TB) [52] calculations and for various analytical potentials (Tersoff T2: Tersoff J. in 1988 [37], Tersoff T3: Tersoff J. in 1989 [53], SW: Stillinger and Weber [54], EDIP: Bazant et al. [55], DS: Dyson and Smith [56]). a0: lattice constant; EC: cohesive energy; B: bulk modulus; cij: elastic constants; S: static shear modulus; EV: vacancy formation energy; ζ: Kleinman parameter.
Table A1. Properties of silicon in its equilibrium structure [49] from experiment [50], quantum-mechanical (QM) methods such as density functional theory (DFT) [51], as well as tight binding (TB) [52] calculations and for various analytical potentials (Tersoff T2: Tersoff J. in 1988 [37], Tersoff T3: Tersoff J. in 1989 [53], SW: Stillinger and Weber [54], EDIP: Bazant et al. [55], DS: Dyson and Smith [56]). a0: lattice constant; EC: cohesive energy; B: bulk modulus; cij: elastic constants; S: static shear modulus; EV: vacancy formation energy; ζ: Kleinman parameter.
PropertiesExpt.QM MethodsAnalytical Potentials
DFTTBTersoff T2Tersoff T3SWEDIPABOPDS
a0 (Å)5.4295.4005.4295.4325.4325.4315.4305.4295.432
EC (eV)−4.63−4.62−4.63−4.62−4.63−4.65−4.63−4.63
C11 (GPa)168159167166143162175167109
C12 (GPa)656167657582626593
C44 (GPa)808575776960717238
B (GPa)99931009898108999998
S (GPa)111119119117112111114
EV (eV)3.173.683.723.702.823.223.20
ζ0.540.530.670.670.630.520.91

Appendix B

The LAMMPS software supports multiple kinds of pressure tensors computation. To avoid the pressure tensors computation that cannot completely fill the periodic boundary. We adopted the atomic Virial tensors to compute the system pressure. The computation of atomic Virial tensors stated in reference [46,47] following:
σ i A t o m = 1 V A t o m m v i v i + 1 2 j ( i ) r i j f i j ,
where σ i A t o m represents the Virial tensors of atom i. vi is the velocity of this atom, rij = rirj is the relative displacement between atom i and j, fij is the interparticle force and represents the tensor product. In the molecular dynamics simulation, the widely used Virial stress is considered equivalent to Cauchy stress. If one requires the pressure tensor σSystem of a specific system, then he must average the atomic Virial tensors in the region/volume (2D/3D system) following:
σ S y s t e m = i N σ i a t o m N ,
where σ i A t o m is the atomic Virial tensors that previously calculated using Equation (A1). N is the total atom count contained in this region/volume.

References

  1. Popel, O.S.; Tarasenko, A.B. Modern development trends in photovoltaics. Therm. Eng. 2021, 68, 807–825. [Google Scholar] [CrossRef]
  2. Punturat, J.; Tangwarodomnukun, V.; Dumkum, C. Surface characteristics and damage of monocrystalline silicon induced by wire-EDM. Appl. Surf. Sci. 2014, 320, 83–92. [Google Scholar] [CrossRef]
  3. Ciszek, T.F.; Wang, T.H. Silicon defect and impurity studies using float-zone crystal growth as a tool. J. Cryst. Growth 2002, 237–239, 1685–1691. [Google Scholar] [CrossRef]
  4. Moore, R.M. Czochralski silicon solar cell modules: Present cost and future prospects. Sol. Cells 1982, 5, 313–329. [Google Scholar] [CrossRef]
  5. Arif, M.; Rahman, M.; San, W.Y. A state-of-the-art review of ductile cutting of silicon wafers for semiconductor and microelectronics industries. Int. J. Adv. Manuf. Tech. 2012, 63, 481–504. [Google Scholar] [CrossRef]
  6. Shen, J.Y.; Zhu, X.; Chen, J.B.; Tao, P.; Wu, X. Investigation on the Edge Chipping in Ultrasonic Assisted Sawing of Monocrystalline Silicon. Micromachines 2019, 10, 616. [Google Scholar] [CrossRef] [Green Version]
  7. Ramachandramoorthy, R.; Schwiedrzik, J.; Petho, L.; Guerra-Nuñez, C.; Frey, D.; Breguet, J.M.; Michler, J. Dynamic plasticity and failure of microscale glass: Rate-dependent ductile-brittle-ductile transition. Nano Lett. 2019, 19, 2350–2359. [Google Scholar] [CrossRef]
  8. Ramachandramoorthy, R.; Gao, W.; Bernal, R.; Espinosa, H. High strain rate tensile testing of silver nanowires: Rate-dependent brittle-to-ductile transition. Nano Lett. 2016, 16, 255–263. [Google Scholar] [CrossRef]
  9. Hu, Z.Q.; Shao, J.L.; Xie, Y.F.; Mei, Y. Atomistic study on tensile fracture of densified silica glass and its dependence on strain rate. Chin. Phys. B 2020, 29, 580–587. [Google Scholar] [CrossRef]
  10. Sansa, M. Frequency fluctuations in silicon nanoresonators. Nat. Nanotechnol. 2016, 11, 552. [Google Scholar] [CrossRef] [Green Version]
  11. Xu, T.; Wu, W.K.; Luo, X.; Guo, L.J. Plasmonic nanoresonators for high-resolution colour filtering and spectral imaging. Nat. Commun. 2010, 1, 59. [Google Scholar] [CrossRef] [PubMed]
  12. Darbaniyan, F.; Yan, X.; Sharma, P. An atomistic perspective on the effect of strain rate and lithium fraction on the mechanical behavior of silicon electrodes. J. Appl. Mech. 2020, 87, 031011. [Google Scholar] [CrossRef] [Green Version]
  13. Pharr, M.; Suo, Z.G.; Vlassak, J.J. Variation of stress with charging rate due to strain-rate sensitivity of silicon electrodes of Li-ion batteries. J. Power Sources 2014, 270, 569–575. [Google Scholar] [CrossRef]
  14. Steinberg, D.J.; Cochran, S.G.; Guinan, M.W. A constitutive model for metals applicable at high-strain rate. J. Appl. Phys. 1980, 51, 1498–1504. [Google Scholar] [CrossRef] [Green Version]
  15. Johnson, G.R. A constitutive model and data for materials subjected to large strains, high strain rates, and high temperatures. In Proceedings of the 7th International Symposium on Ballistics, The Hague, the Netherlands, 19–21 April 1983; pp. 541–547. [Google Scholar]
  16. Bodner, S.R.; Partom, Y. Constitutive equations for elastic-viscoplastic strain-hardening materials. J. Appl. Mech. 1975, 42, 385–389. [Google Scholar] [CrossRef]
  17. Zerilli, F.J.; Armstrong, R.W. Dislocation-mechanics-based constitutive relations for material dynamics calculations. J. Appl. Phys. 1987, 61, 1816–1825. [Google Scholar] [CrossRef] [Green Version]
  18. Hansen, B.L.; Beyerlein, I.J.; Bronkhorst, C.A.; Cerreta, E.K.; Dennis-Koller, D. A dislocation-based multi-rate single crystal plasticity model. Int. J. Plasticity 2013, 44, 129–146. [Google Scholar] [CrossRef]
  19. Guo, Y.; Zhuang, Z.; Li, X.Y.; Chen, Z. An investigation of the combined size and rate effects on the mechanical responses of FCC metals. Int. J. Solids Struct. 2007, 44, 1180–1195. [Google Scholar] [CrossRef] [Green Version]
  20. Kulshreshtha, P.K.; Youssef, K.M.; Rozgonyi, G. Nano-indentation: A tool to investigate crack propagation related phase transitions in PV silicon. Sol. Energy Mater. Sol. Cells 2012, 96, 166–172. [Google Scholar] [CrossRef]
  21. Wang, F.; Huang, P.; Xu, K. Strain rate sensitivity of nanoindentation creep in polycrystalline Al film on Silicon substrate. Surf. Coat. Tech. 2006, 201, 5216–5218. [Google Scholar] [CrossRef]
  22. Jarzbek, D.M.; Milczarek, M.; Nosewicz, S.; Bazarnik, P.; Schift, H. Size effects of hardness and strain rate sensitivity in amorphous silicon measured by nanoindentation. Metall. Mater. Trans. A 2020, 51, 1625–1633. [Google Scholar] [CrossRef] [Green Version]
  23. Yue, Y.; Zheng, K. Strong strain rate effect on the plasticity of amorphous silica nanowires. Appl. Phys. Lett. 2014, 104, 231906. [Google Scholar] [CrossRef]
  24. Smith, R.F.; Minich, R.W.; Rudd, R.E.; Eggert, J.H.; Collins, G.W. Orientation and rate dependence in high strain-rate compression of single-crystal silicon. Phys. Rev. B 2012, 86, 6536–6541. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, B.; Melkote, S.N.; Wang, P.Z.; Saraogi, S. Effect of speed on material removal behavior in scribing of monocrystalline silicon. Precis. Eng. 2020, 66, 315–323. [Google Scholar] [CrossRef]
  26. Zhang, Y.M.; Cao, P.H.; Deng, B.H.; Huang, L.P.; Shi, Y.F. Strain rate-dependent tensile response of glassy silicon nanowires studied by accelerated atomistic simulations. J. Appl. Phys. 2021, 130, 085105. [Google Scholar] [CrossRef]
  27. Xu, W.T.; Kim, W.K. Molecular dynamics simulation of the uniaxial tensile test of silicon nanowires using the MEAM potential. Mech. Mater. 2019, 137, 103140. [Google Scholar] [CrossRef]
  28. Chen, T.; Xu, R.; Li, Q.B. Effect of strain rate on tensile strength of defective silicon nanorods. Acta Mech. Solida Sin. 2015, 28, 133–144. [Google Scholar] [CrossRef]
  29. Jeng, Y.R.; Tsai, P.C.; Fang, T.H. Effects of temperature, strain rate, and vacancies on tensile and fatigue behaviors of silicon-based nanotubes. Phys. Rev. B 2005, 71, 085411. [Google Scholar] [CrossRef]
  30. Jing, Y.H.; Meng, Q.Y.; Zhao, W. Molecular dynamics simulations of the tensile and melting behaviours of silicon nanowires. Physica E 2009, 41, 685–689. [Google Scholar] [CrossRef]
  31. Kang, K.W.; Cai, W. Size and temperature effects on the fracture mechanisms of silicon nanowires: Molecular dynamics simulations. Int. J. Plasticity 2010, 26, 1387–1401. [Google Scholar] [CrossRef]
  32. Voter, A.F.; Montalenti, F.; Germann, T.C. Extending the time scale in atomistic simulation of materials. Annu. Rev. Mater. Res. 2012, 32, 321–346. [Google Scholar] [CrossRef] [Green Version]
  33. Yan, X.; Cao, P.; Tao, W.; Sharma, P.; Park, H.S. Atomistic modeling at experimental strain rates and timescales. J. Phys. Appl. Phys. 2016, 49, 493002. [Google Scholar] [CrossRef] [Green Version]
  34. Kerrache, A.; Mousseau, N.; Lewis, L.J. Crystallization of amorphous silicon induced by mechanical shear deformations. Phys. Rev. B 2011, 84, 014110. [Google Scholar] [CrossRef] [Green Version]
  35. Li, P.H.; Guo, X.G.; Yuan, S.; Li, M.; Kang, R.K.; Guo, D.M. Effects of grinding speeds on the subsurface damage of single crystal silicon based on molecular dynamics simulations. Appl. Surf. Sci. 2021, 554, 149668. [Google Scholar] [CrossRef]
  36. Berendsen, H.J.C.; Postma, J.P.M.; VanGunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684. [Google Scholar] [CrossRef] [Green Version]
  37. Tersoff, J. Empirical interatomic potential for silicon with improved elastic properties. Phys. Rev. B 1988, 38, 9902. [Google Scholar] [CrossRef]
  38. Guillonneau, G.; Mieszala, M.; Wehrs, J.; Schwiedrzik, J.; Grop, S.; Frey, D.; Philippe, L.; Breguet, J.M.; Michler, J.; Wheeler, J.M. Nanomechanical testing at high strain rates: New instrumentation for nanoindentation and microcompression. Mater. Design 2018, 148, 39–48. [Google Scholar] [CrossRef]
  39. Cowper, G.R.; Symonds, P.S. Strain-hardening and strain-rate effects in the impact loading of cantilever beams. Small Bus. Econ. 1957, 31, 235–263. [Google Scholar]
  40. Song, Y.Q.; Cheng, Y.C.; Wang, X.W. Mechanical meaning of strain-rate sensitivity index and criterion of it’s measurement under tension condition. Chin. J. Mech. Eng. 2000, 36, 6. [Google Scholar] [CrossRef]
  41. Stukowski, A. Visualization and analysis of atomistic simulation data with ovito–the open visualization tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012. [Google Scholar] [CrossRef]
  42. Wan, W.; Tang, C.X.; Qiu, A.; Xiang, Y.K. The size effects of point defect on the mechanical properties of monocrystalline silicon: A molecular dynamics study. Materials 2021, 14, 3011. [Google Scholar] [CrossRef] [PubMed]
  43. Wan, W.; Tang, C.X.; Zhang, J.J.; Zhou, L. General molecular dynamics approach to understand the mechanical anisotropy of monocrystalline silicon under the nanoscale effects of point defect. Nanomaterials 2021, 11, 1965. [Google Scholar] [CrossRef] [PubMed]
  44. Stukowski, A.; Bulatov, V.V.; Arsenlis, A. Automated identification and indexing of dislocations in crystal interfaces. Model. Simul. Mater. Sci. Eng. 2012, 20, 85007. [Google Scholar] [CrossRef]
  45. Schiøtz, J.; Vegge, T.; Tolla, F.D.D.; Jacobsen, K.W. Atomic-scale simulations of the mechanical deformation of nanocrystalline metals. Phys. Rev. B 1999, 60, 11971–11983. [Google Scholar] [CrossRef] [Green Version]
  46. Thompson, A.P.; Plimpton, S.J.; Mattson, W. General formulation of pressure and stress tensor for arbitrary many-body interaction potentials under periodic boundary conditions. J. Chem. Phys. 2009, 131, 154107. [Google Scholar] [CrossRef] [Green Version]
  47. Heyes, D.M. Pressure tensor of partial-charge and point-dipole lattices with bulk and surface geometries. Phys. Rev. B 1994, 49, 755. [Google Scholar] [CrossRef]
  48. Goel, S.; Luo, X.C.; Agrawal, A.; Reuben, R.L. Diamond machining of silicon: A review of advances in molecular dynamics simulation. Int. J. Mach. Tool. Manu. 2015, 88, 131–164. [Google Scholar] [CrossRef] [Green Version]
  49. Erhart, P.; Albe, K. Analytical potential for atomistic simulations of silicon, carbon, and silicon carbide. Phys. Rev. B 2005, 71, 035211. [Google Scholar] [CrossRef] [Green Version]
  50. Balamane, H.; Halicioglu, T.; Tiller, W.A. Comparative study of silicon empirical interatomic potentials. Phys. Rev. B 1992, 46, 2250. [Google Scholar] [CrossRef]
  51. Nielsen, O.H.; Martin, R.M. Stresses in semiconductors: Ab initio calculations on Si, Ge, and GaAs. Phys. Rev. B 1985, 32, 3792. [Google Scholar] [CrossRef]
  52. Lenosky, T.J.; Kress, J.D.; Kwon, I.; Voter, A.F.; Edwards, B.; Yang, D.F.; Yang, S.; Adams, J.B. Highly optimized tight-binding model of silicon. Phys. Rev. B 1997, 55, 1528. [Google Scholar] [CrossRef]
  53. Tersoff, J. Modeling solid-state chemistry: Interatomic potentials for multicomponent systems. Phys. Rev. B 1989, 39, 5566. [Google Scholar] [CrossRef] [PubMed]
  54. Stillinger, F.H.; Weber, T.A. Computer simulation of local order in condensed phases of silicon. Phys. Rev. B 1985, 31, 5262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Justo, J.F.; Bazant, M.Z.; Kaxiras, E.; Bulatov, V.V.; Yip, S. Interatomic potential for silicon defects and disordered phases. Phys. Rev. B 1998, 58, 2539. [Google Scholar] [CrossRef] [Green Version]
  56. Dyson, A.J.; Smith, P.V. Extension of the Brenner Empirical Interatomic Potential to C–Si–H Systems. Surf. Sci. 1996, 355, 140–150. [Google Scholar] [CrossRef]
Figure 1. Schematic representations of the model settings.
Figure 1. Schematic representations of the model settings.
Nanomaterials 12 01203 g001
Figure 2. Stress–strain curves of the two structures under the simulated strain rates: (a) Mono-Si, strain rate range: 10−1/ps to 10−3/ps; (a1) Mono-Si, strain rate range: 5 × 10−4/ps to 5 × 10−6/ps; (b) SiNW, strain rate range: 10−1/ps to 10−3/ps; (b1) SiNW, strain rate range: 5 × 10−4/ps to 5 × 10−6/ps.
Figure 2. Stress–strain curves of the two structures under the simulated strain rates: (a) Mono-Si, strain rate range: 10−1/ps to 10−3/ps; (a1) Mono-Si, strain rate range: 5 × 10−4/ps to 5 × 10−6/ps; (b) SiNW, strain rate range: 10−1/ps to 10−3/ps; (b1) SiNW, strain rate range: 5 × 10−4/ps to 5 × 10−6/ps.
Nanomaterials 12 01203 g002
Figure 3. (a) lg ε ˙ versus σs; (b) lg ε ˙ versus lg ( σ s σ q u a s i - s t a t i c 1 ) ; (c) variations of non-diamond atoms and cubic diamond atoms of Mono-Si after fracture; (c1) variations of non-diamond atoms and cubic diamond atoms of SiNW after fracture.
Figure 3. (a) lg ε ˙ versus σs; (b) lg ε ˙ versus lg ( σ s σ q u a s i - s t a t i c 1 ) ; (c) variations of non-diamond atoms and cubic diamond atoms of Mono-Si after fracture; (c1) variations of non-diamond atoms and cubic diamond atoms of SiNW after fracture.
Nanomaterials 12 01203 g003
Figure 4. Cleavage fracture of Mono-Si under three strain rates: (a) 10−2/ps; (b) 10−4/ps; (c) 10−6/ps. Slip fracture of SiNW under three strain rates: (a1) 10−2/ps; (b1) 10−4/ps; (c1) 10−6/ps. The SiNWs were sliced to reveal their inherent structures.
Figure 4. Cleavage fracture of Mono-Si under three strain rates: (a) 10−2/ps; (b) 10−4/ps; (c) 10−6/ps. Slip fracture of SiNW under three strain rates: (a1) 10−2/ps; (b1) 10−4/ps; (c1) 10−6/ps. The SiNWs were sliced to reveal their inherent structures.
Nanomaterials 12 01203 g004
Figure 5. Dislocation and defect mesh in the fracture process of Mono-Si under three strain rates: (a,a1) 10−2/ps; (b,b1) 10−4/ps; (c,c1) 10−6/ps. The atoms were removed for clarity.
Figure 5. Dislocation and defect mesh in the fracture process of Mono-Si under three strain rates: (a,a1) 10−2/ps; (b,b1) 10−4/ps; (c,c1) 10−6/ps. The atoms were removed for clarity.
Nanomaterials 12 01203 g005
Figure 6. Dislocation and defect mesh in the fracture process of SiNW under three strain rates: (a,a1) 10−2/ps; (b,b1) 10−4/ps; (c,c1) 10−6/ps. The atoms were removed for clarity.
Figure 6. Dislocation and defect mesh in the fracture process of SiNW under three strain rates: (a,a1) 10−2/ps; (b,b1) 10−4/ps; (c,c1) 10−6/ps. The atoms were removed for clarity.
Nanomaterials 12 01203 g006
Figure 7. lg ε ˙ versus five equivalent stresses: (a) average stress level of Mono-Si; (a1) maximum stress level of Mono-Si; (b) average stress level of SiNW; (b1) maximum stress level of SiNW.
Figure 7. lg ε ˙ versus five equivalent stresses: (a) average stress level of Mono-Si; (a1) maximum stress level of Mono-Si; (b) average stress level of SiNW; (b1) maximum stress level of SiNW.
Nanomaterials 12 01203 g007
Figure 8. Probability density of von Mises stress distributions in Mono-Si: (a) at three different applied strains and equal strain rates; (b) at three different strain rates and equal applied strains.
Figure 8. Probability density of von Mises stress distributions in Mono-Si: (a) at three different applied strains and equal strain rates; (b) at three different strain rates and equal applied strains.
Nanomaterials 12 01203 g008
Figure 9. Probability density of von Mises stress distributions in SiNW: (a) at three different applied strains and equal strain rates; (b) at three different strain rates and equal applied strains.
Figure 9. Probability density of von Mises stress distributions in SiNW: (a) at three different applied strains and equal strain rates; (b) at three different strain rates and equal applied strains.
Nanomaterials 12 01203 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Y.; Wan, W.; Li, Q.; Xiong, Z.; Tang, C.; Zhou, L. Revisiting the Rate-Dependent Mechanical Response of Typical Silicon Structures via Molecular Dynamics. Nanomaterials 2022, 12, 1203. https://doi.org/10.3390/nano12071203

AMA Style

Liu Y, Wan W, Li Q, Xiong Z, Tang C, Zhou L. Revisiting the Rate-Dependent Mechanical Response of Typical Silicon Structures via Molecular Dynamics. Nanomaterials. 2022; 12(7):1203. https://doi.org/10.3390/nano12071203

Chicago/Turabian Style

Liu, Yi, Wei Wan, Quan Li, Zhenkang Xiong, Changxin Tang, and Lang Zhou. 2022. "Revisiting the Rate-Dependent Mechanical Response of Typical Silicon Structures via Molecular Dynamics" Nanomaterials 12, no. 7: 1203. https://doi.org/10.3390/nano12071203

APA Style

Liu, Y., Wan, W., Li, Q., Xiong, Z., Tang, C., & Zhou, L. (2022). Revisiting the Rate-Dependent Mechanical Response of Typical Silicon Structures via Molecular Dynamics. Nanomaterials, 12(7), 1203. https://doi.org/10.3390/nano12071203

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop