scholarly journals Molecular Dynamics Simulation of the Crystal Orientation and Temperature Influences in the Hardness on Monocrystalline Silicon

2014 ◽  
Vol 2014 ◽  
pp. 1-8 ◽  
Author(s):  
Hongwei Zhao ◽  
Peng Zhang ◽  
Chengli Shi ◽  
Chuang Liu ◽  
Lei Han ◽  
...  

A nanoindentation simulation using molecular dynamic (MD) method was carried out to investigate the hardness behavior of monocrystalline silicon with a spherical diamond indenter. In this study, Tersoff potential was used to model the interaction of silicon atoms in the specimen, and Morse potential was used to model the interaction between silicon atoms in the specimen and carbon atoms in the indenter. Simulation results indicate that the silicon in the indentation zone undergoes phase transformation from diamond cubic structure to body-centred tetragonal and amorphous structure upon loading of the diamond indenter. After the unloading of the indenter, the crystal lattice reconstructs, and the indented surface with a residual dimple forms due to unrecoverable plastic deformation. Comparison of the hardness of three different crystal surfaces of monocrystalline silicon shows that the (0 0 1) surface behaves the hardest, and the (1 1 1) surface behaves the softest. As for the influence of the indentation temperature, simulation results show that the silicon material softens and adhesiveness of silicon increases at higher indentation temperatures.

2013 ◽  
Vol 562-565 ◽  
pp. 802-808 ◽  
Author(s):  
Lin Zhang ◽  
Hong Wei Zhao ◽  
Zhi Chao Ma ◽  
Hu Huang ◽  
Chun Yang Geng ◽  
...  

A series of three-dimensional molecular dynamics (MD) simulations of nanoindentation are conducted to investigate the deformation behavior and phase transformation of monocrystalline silicon with different size hemispherical diamond indenters on (010) crystal plane. The technique of coordination number (CN) is employed to elucidate the detailed mechanism of phase transformation in the monocrystalline silicon. The simulation results show that the phase transformation varies according to the different radii indenters. In the phase transformation region beneath the indenter, the crystalline structures of Si-II, Si-XIII, and amorphous phase structures are observed. In addition, the results indicate that phase transformation with large indenters is not same with the small indenter. The six-coordinated silicon phase, Si-XIII, transformed from Si-I is identified. The phases of Si-II and Si-XIII, which have the same coordinate number, are successfully extracted from the transformation region during nanoindentation and amorphous phase will emerge upon unloading.


2014 ◽  
Vol 513-517 ◽  
pp. 113-116
Author(s):  
Jen Ching Huang ◽  
Fu Jen Cheng ◽  
Chun Song Yang

The Youngs modulus of multilayered nanothin films is an important property. This paper focused to investigate the Youngs Modulus of Multilayered Ni/Cu Multilayered nanoThin Films under different condition by Molecular Dynamics Simulation. The NVT ensemble and COMPASS potential function were employed in the simulation. The multilayered nanothin film contained the Ni and Cu thin films in sequence. From simulation results, it is found that the Youngs modulus of Cu/Ni multilayered nanothin film is different at different lattice orientations, temperatures and strain rate. After experiments, it can be found that the Youngs modulus of multilayered nanothin film in the plane (100) is highest. As thickness of the thin film and system temperature rises, Youngs modulus of multilayered nanothin film is reduced instead. And, the strain rate increases, the Youngs modulus of Cu/Ni multilayered nanothin film will also increase.


2004 ◽  
Vol 15 (05) ◽  
pp. 649-658 ◽  
Author(s):  
SHI-WEI REN

By using the microcanonical molecular dynamics simulation, the melting processes of the clusters bound by Morse potential are investigated. It is found that these clusters show a multi-step melting process as long as the range of the Morse potential is a suitable value. The origins of this multi-step process are analyzed. I find that not only the repulsive core of the potential but also the attractive tail range of the potential influences the melting process. Moreover, the occurrence of the multi-step melting process is more sensitive to the tail region of the Morse potential.


2021 ◽  
Author(s):  
Zhi Meng Zhang ◽  
Hua Yang ◽  
Jun Xia Shi ◽  
Jia Jun Wang ◽  
Zheng Guo Huang ◽  
...  

Abstract The orientation of hydrocarbon chains plays a key role in the applications of organic materials. And chain folding in the process of molecular orientation is also of great significance for the design of organic molecular thin films. The effect of chain length and simulation temperature on the isothermal orientation of n-alkanes on graphene surface is studied by molecular dynamics simulation in this paper. And the chain folding is also described. The n-alkanes can form perpendicular ordered structure, parallel ordered structure or perpendicular orientation at relative low temperature and parallel orientation at relative high temperature on graphene surface. The chain fold happens when long n-alkanes form perpendicular ordered structure on graphene surface. And the simulation results show the interactions of n-alkane−graphene and n-alkane−n-alkane affect chain fold.


2020 ◽  
Vol 35 ◽  
pp. 18-28
Author(s):  
Muhammad Rubayat Bin Shahadat ◽  
A.K.M.M. Morshed

Non-equilibrium molecular dynamics simulations have been employed to study the explosive boiling phenomena of water over a hot copper plate. The molecular system was comprised of three sections: solid copper wall, liquid water, and water vapor. A few layers of the liquid water were placed on the solid Cu surface. The rest of the simulation box was filled with water vapor. Initially, the water molecules were equilibrated by using Berendsen thermostat at 298 K. Then heat was given to the copper plate at different temperatures so that explosive boiling occurs. After achieving the equilibrium by performing the previous two steps, the liquid water at 298 K is suddenly dropped on the hot plate. NVE ensemble was used in the simulation and the temperature of the copper plate was controlled to different temperatures with phantom atom thermostat. Four temperatures (400K, 500K, 650 K and 1000K) were taken to study the explosive boiling. The simulation results show that, the explosive boiling temperature of water on Cu plate is 500 K temperature. At this point, the energy flux was found 1.79x108 J/m3 which is very promising with the experimental results. Moreover, if the temperature of the surface was increased the explosive boiling occurred at a faster rate. The simulation results also show that explosive boiling occurs earlier for the hydrophilic surface than hydrophobic surface as for the hydrophilic surface the water attracted the Cu plate more than the hydrophobic surface and so the amount of energy transfer is more for the hydrophilic surface.


2012 ◽  
Vol 26 (20) ◽  
pp. 1250117 ◽  
Author(s):  
L. T. VINH ◽  
N. V. HUY ◽  
P. K. HUNG

Molecular dynamics simulation is carried out for liquid SiO 2 at pressure ranged from zero to 30 GPa and by using BKS, Born–Mayer type and Morse–Stretch potentials. The constructed models reproduce well the experimental data in terms of mean coordination number, bond angle and pair radial distribution function. Furthermore, the density of all samples can be expressed by a linear function of fractions SiO x. It is found that the topology of units SiO x and linkages OSi y is unchanged upon compression although the liquid undergoes substantial change in its network structure. Consequently, the partial bond angle distribution for SiO x and OSi y is identical for all samples constructed by the same potential. This result allows to establishing a simple expression between total bond angle distribution (BAD) and fraction of SiO x and OSi y. The simulation shows a good agreement between the calculation and simulation results for both total O–Si–O and Si–O–Si BADs. This supports a technique to estimate amount of units SiO x and linkages OSi y on base of total Si–O–Si and O–Si–O BADs measured experimentally.


Author(s):  
G. T. Gao ◽  
J. D. Schall ◽  
K. Van Workum ◽  
P. T. Mikulski ◽  
J. A. Harrison

A constant tension and constant temperature molecular dynamics simulation method was used in the calculations of the elastic constants of the nanocomposite systems. The nanocomposite systems contain a core of sp3 diamond structure surrounded by an amorphous carbon network. The simulation results show that the elastic properties of nanocomposites of diamond-like carbons are closely related to the size of the sp3 diamond core; the bigger the core, the larger the elastic constants, and the system becomes more anisotropic.


Sign in / Sign up

Export Citation Format

Share Document