• Bulletin of the Korean Chemical Society
• /
• v.23 no.3
• /
• pp.441-446
• /
• 2002

We have presented the results of thermodynamic, structural and dynamic properties of liquid benzene, toluene, and p-xylene in canonical (NVT) ensemble at 293.15 K by molecular dynamics (MD) simulations. The molecular model adopted for these molecules is a combination of the rigid body treatment for the benzene ring and an atomistically detailed model for the methyl hydrogen atoms. The calculated pressures are too low in the NVT ensemble MD simulations. The various thermodynamic properties reflect that the intermolecular interactions become stronger as the number of methyl group attached into the benzene ring increases. The pronounced nearest neighbor peak in the center of mass g(r) of liquid benzene at 293.15 K, provides the interpretation that nearest neighbors tend to be perpendicular. Two self-diffusion coefficients of liquid benzene at 293.15 K calculated from MSD and VAC function are in excellent agreement with the experimental measures. The self-diffusion coefficients of liquid toluene also agree well with the experimental ones for toluene in benzene and for toluene in cyclohexane.

• Carbon letters
• /
• v.11 no.4
• /
• pp.316-324
• /
• 2010

The effects of geometrical parameters on mechanical properties of graphite-vinylester nanocomposites and their constituents (matrix, reinforcement and interface) are studied using molecular dynamics (MD) simulations. Young's modulii of 1.3 TPa and 1.16 TPa are obtained for graphene layer and for graphite layers respectively. Interfacial shear strength resulting from the molecular dynamic (MD) simulations for graphene-vinylester is found to be 256 MPa compared to 126 MPa for graphitevinylester. MD simulations prove that exfoliation improves mechanical properties of graphite nanoplatelet vinylester nanocomposites. Also, the effects of bromination on the mechanical properties of vinylester and interfacial strength of the graphene.brominated vinylester nanocomposites are investigated. MD simulation revealed that, although there is minimal effect of bromination on mechanical properties of pure vinylester, bromination tends to enhance interfacial shear strength between graphite-brominated vinylester/graphene-brominated vinylester in a considerable magnitude.

• International Journal of Automotive Technology
• /
• v.5 no.3
• /
• pp.215-219
• /
• 2004

A vehicle structure needs to be more precisely analyzed because of complexities and varieties. Structural fatigue which is generated by fluctuations of stresses during the service life of a mechanical system is the primary concern in the structural design for safety. A fatigue life is difficult to obtain in structural components during the service life of mechanical systems since the fluctuating stress contributes to fatigue. This study introduces new procedures to measure the lethargy coefficient and to predict the fatigue life of a mechanical structure by using molecular dynamic simulation. A lethargy coefficient is the total defect-estimating coefficient, which was obtained by using the results of a simple tensile test in this study. With this lethargy coefficient, fatigue life was estimated. The proposed method will be useful in predicting the fatigue life of a structurally-modified vehicle design. The effectiveness of the proposed method using lethargy coefficient measurement to predict the fatigue life of a structure was examined by applying this method to predict the fatigue life of SS41 steel, used extensively as material of vehicle structures. Two types of specimen such as pre-cracked plate and simple plate is discussed. equation of fatigue life using the lethargy coefficient and failure time, both obtained from a simple tensile test, will be useful in engineering. This measurement and prediction technology will be extended for use in analysis of any geometric shapes of modified automotive structures.

• Korean Journal of Metals and Materials
• /
• v.47 no.1
• /
• pp.32-37
• /
• 2009

Molecular dynamic simulations on the structural evolution of the Si(001) vicinal surfaces, which are tilted with respect to [100] and [110] directions were performed by using the empirical Tersoff potential. Tersoff potential was implemented at LAMMPS code and confirmed to describe the properties of Si. When the steps are generated along [100] direction, symmetric dimer rows formed with respect to the step edges. On the other hand, when the steps are generated along [110] direction, alternating dimer rows form with respect to the step edges. The configurational differences between the two vicinal surfaces were discussed in terms of the surface diffusion and the possibility of preventing step bunching for the (001) vicinal surface tilted along [100] direction was suggested.

• Journal of the Korean Vacuum Society
• /
• v.15 no.4
• /
• pp.410-420
• /
• 2006

The dynamic effects, such as the steering and the screening effects during deposition on an epitaxial growth is studied by kinetic Monte Carlo simulation. In the simulation, we incorporates molecular dynamic simulation to rigorously take the interaction of the deposited atom with the substrate atoms into account, We find three characteristic features of the surface morphology developed by grazing angle deposition: (1) enhanced surface roughness, (2) asymmetric mound, and (3) asymmetric slopes of mound sides, Regarding their dependence on both deposition angle and substrate temperature, a reasonable agreement of the simulated results with the previous experimental ones is found. The characteristic growth features by grazing angle deposition are mainly caused by the inhomogeneous deposition flux due to the steering and screening effects, where the steering effects play the major role rather than the screening effects. Newly observed in the present simulation is that the side of mound in each direction is composed of various facets instead of all being in one selected mound angle even if the slope selection is attained, and that the slope selection does not necessarily mean the facet selection.

• Polymer(Korea)
• /
• v.35 no.6
• /
• pp.513-519
• /
• 2011

Molecular simulations via the molecular dynamics method have been carried out to investigate the dynamic collision properties of penetrable-sphere model fluids. The collision frequencies, the mean free paths, the angle distributions of the hard-type reflection and the soft-type penetration, and the effective packing fractions are computed over a wide range of the packing fraction ${\phi}$ and the repulsive energy ${\varepsilon}^*$. The soft-type collisions are dominated for lower repulsive energy systems, while the hardtype collisions for higher repulsive energy systems. Very interestingly, the ratio of the soft-type (or, the hard-type) collision frequency to the total collision frequency is directly related with the Boltzmann factor of acceptance (or rejection) probabilities in the canonical ensemble Monte Carlo calculations. Such dynamic collision properties are shown to be restricted for highly repulsive and dense systems of ${\varepsilon}^*{\geqq}3.0 $and ${\phi}{\geqq}0.7$, indicating the cluster forming structures in the penetrable-sphere model.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2003.10a
• /
• pp.339-346
• /
• 2003

Molecular dynamic simulations of nano indentation on single-crystal silicon (100) surface were performed using diamond indentor. Silicon substrate and diamond indentor were modeled diamond structure with Tersoff potential model. Phase transformation of silicon, incipient plastic deformation, change of incident temperature distribution are investigated through the change of potential energy distribution, displacement-load diagram, the change of kinetic energy distribution and displacements of silicon atoms. Phase transformation is highly localized and consists of a high-density region surrounding the tip. Axial load linearly increased according to the indenting depth. Number of atoms with high kinetic energy increased at the interface between substrate and indentor tip.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.19 no.12
• /
• pp.883-888
• /
• 2007

As a starting point of investigating what molecular dynamic simulations can reveal about the nature of atomic level of heating and cooling process, argon described by the LJ potential is considered. Stepwise heating and cooling of constant rates are simulated in the NPT (constant number, pressure and temperature) ensemble. Hysteresis is found due to the superheating and supercooling. Drastic change of volume and energy is involved with phase change, but the melting point can not be obtained by simply observing the changes of these quantities. Since liquid and solid phases can co-exist at the same temperature, Gibbs free energy should be calculated to find the temperature where the Gibbs free energy of liquid is equal to that of the solid since the equilibrium state is the state of minimum Gibbs free energy. The obtained melting temperature, $T^*=0.685$, is close to that of the experiment with only 2% error.

• 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 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.