• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.331-334
• /
• 2009

In this paper, we present a surface relaxation model in atomistic calculations for thin nanofilms. This surface relaxation model is very simple model which have only two degrees of freedoms to determine the atomic positions of nanofilms. Whereas in conventional molecular statics simulations, the same number of degrees of freedoms at all atom positions are used as unknown variables. In order to prove the reliability of the presented model, we present the results of self-equilibrium strain calculations with the surface parameters obtained from this model.

• Bulletin of the Korean Chemical Society
• /
• v.33 no.6
• /
• pp.1913-1918
• /
• 2012

The nitramine explosives with $-NH_2$ and -F groups were optimized to obtain their molecular geometries and electronic structures at DFT-B3LYP/6-31+G(d) level. The theoretical molecular density (${\rho}$), heat of formation (HOF), detonation velocity ($D$) and detonation pressure ($P$), estimated using Kamlet-Jacobs equations, showed that the detonation properties of these compounds were excellent. Based on the frequencies scaled by 0.96 and the principle of statistic thermodynamics, the thermodynamic properties were evaluated, which were respectively related with the temperature. The simulation results reveal that 1,3,5,7-tetranitro-1,3,5,7-tetrazocan-2-amine (molecule B1) performs similarly to the famous explosive HMX, and 2-fluoro-1,3,5-trinitro-1,3,5-triazinane (molecule C1) and 2-fluoro-1,3,5,7-tetranitro-1,3,5,7-tetrazocane (molecule D1) outperform HMX. According to the quantitative standard of energetics and stability as an HEDC (high energy density compound), molecules C1 and D1 essentially satisfy this requirement. These results provide basic information for molecular design of novel high energetic density compounds.

• Bulletin of the Korean Chemical Society
• /
• v.24 no.11
• /
• pp.1590-1598
• /
• 2003

In this paper we have presented the results for viscosity and self-diffusion constants of model systems for four liquid n-alkanes ($C_{12}, C_{20}, C_{32}, and C_{44}$) in a canonical ensemble at several temperatures using molecular dynamics (MD) simulations. The small chains of these n-alkanes are clearly $<{R_{ee}}^2>/6<{R_g}^2>>1$, which leads to the conclusion that the liquid n-alkanes over the whole temperatures considered are far away from the Rouse regime. Calculated viscosity ${\eta}$ and self-diffusion constants D are comparable with experimental results and the temperature dependence of both ${\eta}$ and D is suitably described by the Arrhenius plot. The behavior of both activation energies, $E_{\eta}$ and $E_D$, with increasing chain length indicates that the activation energies approach asymptotic values as n increases to the higher value, which is experimentally observed. Two calculated monomeric friction constants ${\zeta}$ and ${\zeta}_D$ give a correct qualitative trend: decrease with increasing temperature and increase with increasing chain length n. Comparison of the time auto-correlation functions of the end-to-end vector calculated from the Rouse model for n-dodecane ($C_{12}$) at 273 K and for n-tetratetracontane ($C_{44}$) at 473 K with those extracted directly from our MD simulations confirms that the short chain n-alkanes considered in this study are far away from the Rouse regime.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.17 no.11
• /
• pp.1252-1257
• /
• 2004

This paper describes the information for quantitative simulation of weakly ionized plasma. In previous paper, we calculated the electron transport coefficients by using two-term approximation of Boltzmann equation. But there is difference between the result of the two-term approximation of the Boltzmann equation and experiments in pure CF$_4$ molecular gas and in CF$_4$+Ar gas mixture. Therefore, In this paper, we calculated the electron drift velocity (W) in pure CF$_4$ molecular gas and CF$_4$+Ar gas mixture (1 %, 5 %, 10 %) for range of E/N values from 0.17～300 Td at the temperature was 300 K and pressure was 1 Torr by multi-term approximation of the Boltzmann equation by Robson and Ness. The results of two-term and multi-term approximation of the Boltzmann equation have been compared with each other for a range of E/N.

• Coupled systems mechanics
• /
• v.6 no.1
• /
• pp.41-54
• /
• 2017

Shape-memory alloys (SMA) have interesting behaviors and important mechanical properties due to the solid-solid phase transformation. These phenomena are dominated by the evolution of microstructures. In recent years, the microstructures in SMAs have been studied extensively and modeled using molecular dynamics (MD) simulations. However, it remains difficult to identify the crystal variants in the simulation results, which consist of large numbers of atoms. In the present work, a method is developed to identify the austenite phase and the monoclinic martensite crystal variants in MD results. The transformation matrix of each lattice is calculated to determine the corresponding crystal variant. Evolution of the volume fraction of the crystal variants and the microstructure in Ni-Ti SMAs under thermal and mechanical boundary conditions are examined. The method is validated by comparing MD-simulated interface normals with theoretical solutions. In addition, the results show that, in certain cases, the interatomic potential used in the current study leads to inconsistent monoclinic lattices compared with crystallographic theory. Thus, a specific modification is applied and the applicability of the potential is discussed.

• Coupled systems mechanics
• /
• v.3 no.4
• /
• pp.329-344
• /
• 2014

In this work, quantum molecular dynamics simulations (QMD) are preformed to study the hydrogen molecules in three types of carbon nanostructures, $C_{60}$ fullerene, (5,5) and (9,0) carbon nanotubes and graphene layers. Interactions between hydrogen and the nanostructures is of importance to understand hydrogen storage for the development of hydrogen economy. The QMD method overcomes the difficulties with empirical interatomic potentials to model the interaction among hydrogen and carbon atoms in the confined geometry. In QMD, the interatomic forces are calculated by solving the Schrodinger's equation with the density functional theory (DFT) formulation, and the positions of the atomic nucleus are calculated with the Newton's second law in accordance with the Born-Oppenheimer approximation. It is found that the number of hydrogen atoms that is less than 58 can be stored in the $C_{60}$ fullerene. With larger carbon fullerenes, more hydrogen may be stored. For hydrogen molecules passing though the fullerene, a particular orientation is required to obtain least energy barrier. For carbon nanotubes and graphene, adsorption may adhere hydrogen atoms to carbon atoms. In addition, hydrogen molecules can also be stored inside the nanotubes or between the adjacent layers in graphite, multi-layer graphene.

• Current Optics and Photonics
• /
• v.8 no.1
• /
• pp.86-96
• /
• 2024

This research investigation employs a terahertz (THz) time-domain spectroscopy system to study the terahertz spectral characteristics of five different citrates in both solution and solid state. The citrates under examination are lithium citrate, monosodium citrate, disodium citrate, trisodium citrate, and potassium citrate. The results show that the THz absorption coefficients of the first four citrate solutions exhibit a decreasing trend with increasing concentration. However, the potassium citrate solution shows an opposite phenomenon. At the same time, the absorption coefficients of lithium citrate, trisodium citrate, and potassium citrate solutions are compared at the same concentration. The results indicate that the absorption coefficient of citrate solution increases in proportion to the increase of metal cation radius, which is explained from the perspective of the influence of metal cations on hydrogen bonds. In addition, we also study the absorption peaks of solid citrates, and characterize the formation mechanism of the absorption peaks by molecular dynamics simulations. This methodology can be further extended to the study of multitudinous salts, presenting theoretical foundations for the detection in food and medicine industries.

• Interaction and multiscale mechanics
• /
• v.6 no.3
• /
• pp.301-316
• /
• 2013

The influence of indenter geometry on nanoindentation was studied using a static molecular dynamics simulation. Dislocation nucleation, dislocation locks, and dislocation movements during nanoindentation into Al (001) were studied. Spherical, rectangular, and Berkovich indenters were modeled to study the material behaviors and dislocation activities induced by their different shapes. We found that the elastic responses for the three cases agreed well with those predicted from elastic contact theory. Complicated stress fields were generated by the rectangular and Berkovich indenters, leading to a few uncommon nucleation and dislocation processes. The calculated mean critical resolved shear stresses for the Berkovich and rectangular indenters were lower than the theoretical strength. In the Berkovich indenter case, an amorphous region was observed directly below the indenter tip. In the rectangular indenter case, we observed that some dislocation loops nucleated on the plane. Furthermore, a prismatic loop originating from inside the material glided upward to create a mesa on the indenting surface. We observed an unusual softening phenomenon in the rectangular indenter case and proposed that heterogeneously nucleating dislocations are responsible for this.

• Bulletin of the Korean Chemical Society
• /
• v.25 no.4
• /
• pp.533-538
• /
• 2004

Canonical sampling method has been presented to generate the initial conditions for reactive flux studies of organic reactions in water. Velocity Verlet version of Nose-Hoover chain dynamics algorithm has been employed to sample the initial conditions according to canonical distribution. The unstable normal mode of a transition state has been introduced to define a dividing plane separating reactant and product regions in reaction processes. This method has been implemented and tested for the case iels-Alder reaction of methyl vinyl ketone (MVK) and cyclopentadiene (CPD) in water, providing a reliable tool for further reactive flux molecular dynamics studies in condensed media.

• Bulletin of the Korean Chemical Society
• /
• v.35 no.5
• /
• pp.1501-1505
• /
• 2014

Molecular dynamics simulations were performed to understand the structural arrangement of water molecules around highly charged nanoparticles under aqueous conditions. The effect of two highly charged nanoparticles on the solvation charge asymmetry has been examined. We calculated the radial distribution functions of the components of water molecules around nanoparticles which have four charge types at two different salt concentrations. Even though the distributions of water molecules surrounding a sodium ion and a chloride ion are hardly affected by the charges of nanoparticles and the salt concentrations, those around highly charged nanoparticles are strongly influenced by the charges of nanoparticles, but hardly by the charges of nanoparticles and salt concentrations. We find that the distributions of hydrogen atoms in water molecules around one highly charged nanoparticle are dependent on the magnitude of the nanoparticle charge.