• Bulletin of the Korean Chemical Society
• /
• v.34 no.12
• /
• pp.3800-3804
• /
• 2013

In this study, molecular dynamics simulations of SPC/E (extended simple point charge) model have been carried out in the canonical NVT ensemble over the range of temperatures 300 to 550 K with and without Ewald summation. The quaternion method was used for the rotational motion of the rigid water molecule. Radial distribution functions $g_{OO}(r)$, $g_{OH}(r)$, and $g_{HH}(r)$ and self-diffusion coefficients D for SPC/E water were determined at 300-550 K and compared to experimental data. The temperature dependence on the structural and diffusion properties of SPC/E water was discussed.

• Transactions of the Korean hydrogen and new energy society
• /
• v.11 no.3
• /
• pp.127-136
• /
• 2000

Molecular dynamics simulations have been carried out to investigate the scattering and penetration properties of hydrogen ions with the normal incident angle to Ni (100) surface. The initial kinetic energies of hydrogen ions range from 100 to 1,600 eV. The simulation results are used to assess the applicabilities of theoretical predictions based on the binary collision approximation, and, in the high kinetic regime, theoretical results for scattering energies were shown to he a good agreement with molecular simulations. The angle dependencies on both scattering and penetration distributions were found in the longitudinal direction, but not in the azimuthal direction except for the high kinetic energy of 1,600 eV.

• Proceedings of the IEEK Conference
• /
• 2003.07b
• /
• pp.735-738
• /
• 2003

In this paper, we report a molecular dynamics (MD) simulation of the ion implantation for nano-scale devices with ultra-shallow junctions. In order to model the profile of ion distribution in nanometer scale, the molecular dynamics with a damage model has been employed. As an exemplary case, we calculate the dopant profile during the ion implantation of B, As, and Ge.

• Bulletin of the Korean Chemical Society
• /
• v.27 no.8
• /
• pp.1154-1158
• /
• 2006

Anomalously high ionic mobilities of H+ and $OH^-$ are owing to the transfer of $H^+$ by the Grotthus chain mechanism. Molecular dynamics simulations for the system of 215 water including $OH^-$ ion at 298.15 K using the OSS2 model [J. Chem. Phys. 109, 5547 (1998)] as a dissociable water model with the use of Ewald summation were carried out in order to study the dynamics of $OH^-$ in water. The calculated ionic mobility of $OH^-$ is in good agreement with the experimental result and the Grotthus chain mechanism is fully understood.

• Bulletin of the Korean Chemical Society
• /
• v.35 no.12
• /
• pp.3527-3531
• /
• 2014

In this paper, we report thermodynamic and transport properties (diffusion coefficient, viscosity, and thermal conductivity) of diatomic gases ($H_2$, $N_2$, $O_2$, and $Cl_2$) at 273.15 K and 1.00 atm by performing molecular dynamics simulations using Lennard-Jones intermolecular potential and modified Green-Kubo formulas. The results of self-diffusion coefficients of diatomic gases obtained from velocity auto-correlation functions by Green-Kubo relation are in good agreement with those obtained from mean square displacements by Einstein relation. While the results for viscosities of diatomic gases obtained from stress auto-correlation functions underestimate the experimental results, those for thermal conductivities obtained from heat flux auto-correlation functions overestimate the experimental data except $H_2$.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.8
• /
• pp.826-832
• /
• 2007

Fully flexible cell preserves Hamiltonian in structure so that the symplectic time integrator is applicable to the equations of motion. In the direct formulation of fully flexible cell N-Sigma-T ensemble, a generalized leapfrog time integration (GLF) is applicable for fully flexible cell simulation, but the equations of motion by GLF has structure of implicit algorithm. In this paper, the time integration formula is derived for the fully flexible cell molecular dynamics simulation by using the splitting time integration. It separates flexible cell Hamiltonian into terms corresponding to each of Hamiltonian term. Thus the simple and completely explicit recursion formula was obtained. We compare the performance and the result of present splitting time integration with those of the implicit generalized leapfrog time integration.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.32 no.12
• /
• pp.1047-1054
• /
• 2008

In this paper, the molecular dynamics (MD) simulations are performed with single-crystal copper blocks under simple shear and simple tension to investigate the effect of size and crystal orientation. There are many variances to give influences such as deformation path, temperature, specimen size and crystal orientation. Among them, the crystal orientation has a primary influence on the volume averaged stress. The numerical results show that the volume averaged shear stress decreases as the specimen size increases and as the crystal orientation changes from single to octal. Furthermore, the Schmid factor and yield stress for crystal orientation are evaluated by using the MD simulation on the standard triangle of stereographic projection.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.228-232
• /
• 2003

It has recently been shown that the instability of thin film of a nanoscale can be used in the processes of building nano-size structures, which have potential practical importance in nanotechnology. Molecular dynamics simulation is conducted to probe the thin fluid film of a nano-size and its dynamic behavior during destabilization and structure formation. Non-continuum characteristics are shown in the properties like pressure tensor, viscosity, and thermal conductivity. The thermocapillary force induces a slow growth of long waves in the scale considered. A long-range interaction with the solid wall induces vertical structures, whose formation time and space between neighbors are proportional to the strength of the interaction.

• Transactions on Electrical and Electronic Materials
• /
• v.4 no.1
• /
• pp.1-6
• /
• 2003

We have studied the size confinement effect on the properties of melting-like transition of small icosahedral copper clusters using a classical molecular dynamics simulation based on a well fitted empirical potential. We investigated the caloric curves of icosahedron nanoclusters and the significant depression in the melting temperatures of the copper nanoclusters was compared with that of the bulk copper. A structural transitions from decahedral to icosahedral shapes were shown. As the cluster size increased, the melting temperature increased, and the latent heat increased but seem to be saturated. However, the specific heat was unrelated to the cluster size.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.16 no.6
• /
• pp.454-459
• /
• 2003

We investigate the structure and properties of SiC (Silicon Carbide) nanotubes using molecular dynamics simulation based on the Tersoff bond-order potential. For small diameter tubes, the Si-C bond distance of SiC nanotubes decreases as the nanotube diameter is decreased, due to curvature of the nanotube surface. We find that Young's modulus of SiC nanotubes is somewhat smaller than that of the other nanotubes considered so far. However, Young's modulus for SiC nanotubes is larger than that of ${\beta}$-SiC and almost equal to the experimental value for SiC nanorod and SiC whisker. The strain energy of the SiC nanotubes is also lower than that of the other nanotubes. The lower strain energy of SiC nanotubes raises the possibility of synthesis of SiC nanotubes.