• Proceedings of the IEEK Conference
• /
• 1998.06a
• /
• pp.421-424
• /
• 1998

In this paper are presented a newly proposed 3D monte carlo (MC) damage model for the dynamic simulation in order to more accurately and consistently predict the implant-induced point defect distributions of the various ions in crystalline silicon. This model was applied to phosphorus implants for the ULSI CMOS technology developement. In additon, a newly applied 3D-trajectory split method has been implemented into our model to reduce the statistical fluctuations of the implanted impurity and the defect profiles in the relatively large implanted area as compared to 1D or 2D simulations. Also, an empirical electronic energy loss model is proposed for phosphorus and silicon implants. The 3D formations of the amorphous region and the ultra-shallow junction around the implanted region could be predicted by using our model, TRICSI(Transport ions into crystal-silicon).

• Journal of the Korean Ceramic Society
• /
• v.44 no.8
• /
• pp.412-418
• /
• 2007

The atomic scale details of the melting of solid argon were monitored with the aid of molecular dynamics simulations. The potential energy distribution is substantially disturbed by an increase in the interatomic distance and the random of set distance from the lattice points, with increasing temperature. The potential energy barriers between the lattice points decrease in magnitude with the temperature. Eventually, at the melting point, these barriers can be overcome by atoms that are excited with the entropy gain acquired when the atoms obtain rotational freedom in their atomic motion, and the rotational freedom leads to the collapse of the crystal structure. Furthermore, it was found that the surface of crystals plays an important role in the melting process: the surface eliminates the barrier for the nucleation of the liquid phase and facilitates the melting process. Moreover, the atomic structure of the surface varies with increasing temperature, first via surface roughening and then, before the bulk melts, via surface melting.

• Carbon letters
• /
• v.16 no.3
• /
• pp.183-191
• /
• 2015

Gas transport through graphene-derived membranes has gained much interest recently due to its promising potential in filtration and separation applications. In this work, we explore Kr-85 gas radionuclide sequestration from natural air in nanoporous graphene oxide membranes in which different sizes and geometries of pores were modeled on the graphene oxide sheet. This was done using atomistic simulations considering mean-squared displacement, diffusion coefficient, number of crossed species of gases through nanoporous graphene oxide, and flow through interlayer galleries. The results showed that the gas features have the densest adsorbed zone in nanoporous graphene oxide, compared with a graphene membrane, and that graphene oxide was more favorable than graphene for Kr separation. The aim of this paper is to show that for the well-defined pore size called P-7, it is possible to separate Kr-85 from a gas mixture containing Kr-85, O₂ and N₂. The results would benefit the oil industry among others.

• Proceedings of the IEEK Conference
• /
• 2003.07b
• /
• pp.731-734
• /
• 2003

An atomistic process modeling, Kinetic Monte Carlo simulation, has the advantage of being both conceptually simple and extremely powerful. Instead of diffusion equations, it is based on the definitions of the interactions between individual atoms and defects. Those interactions can be derived either directly from molecular dynamics, first principles calculations, or from experiment. In this paper, as a simple illustration of the kinetic Monte Carlo we simulate defects (self-interstitials and vacancies) diffusion after ion implantation in Si crystalline.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1083-1088
• /
• 2003

Molecular dynamics simulations on the deformation behavior of single-walled carbon nanotube are performed. Formation energies of CNT's by interatomic potentials are computed and compared with ab initio results. Bending and axial compression are applied under lattice statics and NVT ensemble conditions. Specifically, we focus on the mechanism of kink formation in bending. The simulation results are comprehensively explained in the framework of atomistic energetics. The effects of temperature and chirality on the deformation of carbon nanotube are also studied.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.03a
• /
• pp.59-62
• /
• 1997

Recently, quasimolecular dynamics has been successfully used to simulate the deformation characteristic of actual size material. In quasimolecular dynamics, which is an attempt to bridge the gab between atomistic and continuum simulations, molecules are aggregated into large units, called quasimolecules, to simulate the large scale material behavior. In this paper, a numerical simulation using quasimolecular dynamics has been performed to investigate the laminar composite material fracture and crack propagation behaviors in bending process of laminar composite material which is made of fictitious materials. The simulation of the bending of laminar composite material has clarified the effects of strength of material at outer surface upon the fracture behviors of the specimen.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2010.08a
• /
• pp.93-93
• /
• 2010

In this talk, I will introduce a reactive force field (ReaxFF) molecular dynamics (MD) simulation. In contrast to common MD simulations with empirical FFs, we can predict chemical reactions (bond breaking and formation) in large scale systems with the ReaxFF simulation where all of the ReaxFF parameters are from quantum mechanical calculations such as density functional theory to provide high accuracy. Accordingly, the ReaxFF simulation provides both accuracy of quantum mechanical calculations and description of large scale systems of atomistic simulations at the same time. Here, I will first discuss a theory in the ReaxFF including the differences from other empirical FFs, and then show several applications for studying chemical reactions of SiHx radicals on Si surfaces, which is an important issue in Si process.

• Journal of Powder Materials
• /
• v.22 no.4
• /
• pp.247-253
• /
• 2015

The sintering mechanisms of nanoscale copper powders have been investigated. A molecular dynamics (MD) simulation with the embedded-atom method (EAM) was employed for these simulations. The dimensional changes for initial-stage sintering such as characteristic lengths, neck growth, and neck angle were calculated to understand the densification behavior of copper nano-powders. Factors affecting sintering such as the temperature, powder size, and crystalline misalignment between adjacent powders have also been studied. These results could provide information of setting the processing cycles and material designs applicable to nano-powders. In addition, it is expected that MD simulation will be a foundation for the multi-scale modeling in sintering process.

• Elastomers and Composites
• /
• v.44 no.2
• /
• pp.112-115
• /
• 2009

In nature, some creatures have the ability to walk vertically or upside down on various natural surfaces. One of them, Tokay Gecko has billions of nano-hair at the epidermis of its digital. These nano hairs makes atomistic attractive force to the surface, therefore, it could walk and run on the vertical or upside surfaces. Recently, many researchers tried to fabricate the dry adhesive structure mimicking the nano-hair structure. in this study, I tried to survey these studies to discuss the direction of future fabrication works of dry adhesive structures.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.127-136
• /
• 2013

In this study, atomistic simulations are performed to study the effect of Al solute on the behaviour of edge dislocation in Mg alloys. After the dissociation of an Mg basal edge dislocation into two Shockley partials using molecular mechanics, the interaction between the dislocation and Al solute at different temperatures is studied using molecular dynamics. It appears from the simulations that the critical shear stress increases with the Al solute concentration. Comparing with the solute effect at T = 0 K, however, the critical shear stress at a finite temperature is lower since the kinetic energy of the atoms can help the dislocation conquer the energy barriers created by the Al atoms. The velocity of the edge dislocation decreases as the Al concentration increases when the external shear stress is relatively small regardless of temperature. The Al concentration effect on the dislocation velocity is not significant at very high shear stress level when the solute concentration is below 4.0 at%. Drag coefficient B increases with the Al concentration when the stress to temperature ratio is below 0.3 MPa/K, although the effect is more significant at low temperatures.