• Interaction and multiscale mechanics
• /
• v.6 no.4
• /
• pp.395-409
• /
• 2013

Molecular dynamics (MD) systems are highly nonlinear and nonlocal, and the conventional model order reduction methods are ineffective for MD systems. The RBF-POD method (Lee and Chen, 2013) employed a radial basis function (RBF) approximated potential energies and inter-atomic forces of MD systems under the framework of the proper orthogonal decomposition (POD) method for the reduced-order modeling of MD systems. In this work, we focus on the numerical procedures of the RBF-POD method and demonstrate how to apply this approach to the modeling of ds-DNA molecules under stretching and bending conditions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.2
• /
• pp.101-106
• /
• 2009

The disjoining pressure is an important physical property in modeling the small-scale transport phenomena on thin film. It is a very useful definition in characterizing the non-continuum effects that are not negligible in heat and mass transport of the film thinner than submicro-scales. We present the calculated values of disjoining pressure of He, Kr and Xe thin films absorbed on graphite substrate using Molecular Dynamics Simulation (MD). The disjoining pressure is accurately calculated in the resolution of a molecular scale of the film thickness. The characteristics of the pressure are discussed regarding the molecular nature of the fluid system such as molecular diameter and intermolecular interaction parameters. The MD results are also compared with those based on the continuum approximation of the slab-like density profile and the results on other novel gases in the previous study. The discrepancies of the continuum model with MD results are shown in all three configurations and discussed in the view point of molecular features.

• Bulletin of the Korean Chemical Society
• /
• v.18 no.5
• /
• pp.501-509
• /
• 1997

In recent papers[Bull. Kor. Chem. Soc. 1996, 17, 735; ibid 1997, 18, 478] we reported results of molecular dynamics (MD) simulations for the thermodynamic, structural, and dynamic properties of liquid normal alkanes, from n-butane to n-heptadecane, using three different models. Two of the three classes of models are collapsed atomic models while the third class is an atomistically detailed model. In the present paper we present results of MD simulations for the corresponding properties of liquid branched-chain alkanes using the same models. The thermodynamic property reflects that the intermolecular interactions become weaker as the shape of the molecule tends to approach that of a sphere and the surface area decreases with branching. Not like observed in the straight-chain alkanes, the structural properties of model Ⅲ from the site-site radial distribution function, the distribution functions of the average end-to-end distance and the root-mean-squared radii of gyration are not much different from those of models Ⅰ and Ⅱ. The branching effect on the self diffusion of liquid alkanes is well predicted from our MD simulation results but not on the viscosity and thermal conductivity.

• Journal of Korean Vacuum Science & Technology
• /
• v.6 no.1
• /
• pp.22-29
• /
• 2002

The process of energetic deposition has been simulated by using molecular dynamics (W) simulation and kinetic Monte Carlo (KMC) simulation. The atomistic mechanisms of film growth in energetic deposition are discussed.

• Tribology and Lubricants
• /
• v.32 no.1
• /
• pp.9-17
• /
• 2016

Contact between solid surfaces is one of the most important factors that influence dynamic behavior in micro/nanoscale. Although numerous theories and experimental results on contact behavior have been proposed, a thorough investigation for nanomaterials is still not available owing to technical difficulties. Therefore, molecular dynamics simulation was performed to investigate the contact behavior of nanomaterials, and the application of conventional contact theories to nanoscale was assessed in this work. Particularly, the contact characteristics of cylindrical nanowires were examined via simulation and contact theories. For theoretical analysis, various contact models were utilized and work of adhesion, Hamaker constant and elastic modulus those are required for calculation of the models were obtained from both indentation simulation and tensile simulation. The contact area of the cylindrical nanowire was assessed directly through molecular dynamics simulation and compared with the results obtained from the theories. Determination of the contact area of the nanowires was carried out via simulation by counting each atom, which is within the equilibrium length. The results of the simulation and theoretical calculations were compared, and it was estimated that the discrepancy in the results calculated between the simulation and the theories was less than 10 except in the case of the smallest nanowires. As the result, it was revealed that contact models can be effectively utilized to assess the contact area of nanomaterials.

• 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 the Korean Magnetic Resonance Society
• /
• v.19 no.1
• /
• pp.1-10
• /
• 2015

Refinements with atomistic molecular dynamics (MD) simulation have contributed to improving the qualities of NMR structures. In most cases, the calculations with atomistic MD simulation for NMR structures employ generalized-Born implicit solvent model (GBIS) to take into accounts solvation effects. Developments in algorithms and computational capacities have ameliorated GBIS to approximate solvation effects that explicit solvents bring about. However, the quantitative comparison of NMR structures in the latest GBIS and explicit solvents is lacking. In this study, we report the direct comparison of NMR structures that atomistic MD simulation coupled with GBIS and water molecules refined. Two model proteins, GB1 and ubiquitin, were recalculated with experimental distance and torsion angle restraints, under a series of simulated annealing time steps. Whereas the root mean square deviations of the resulting structures were apparently similar, AMBER energies, the most favored regions in Ramachandran plot, and MolProbity clash scores witnessed that GBIS-refined structures had the better geometries. The outperformance by GBIS was distinct in the structure calculations with sparse experimental restraints. We show that the superiority stemmed, at least in parts, from the inclusion of all the pairs of non-bonded interactions. The shorter computational times with GBIS than those for explicit solvents makes GBIS a powerful method for improving structural qualities particularly under the conditions that experimental restraints are insufficient. We also propose a method to separate the native-like folds from non-violating diverged structures.

• Journal of the Korea Institute of Information Security & Cryptology
• /
• v.13 no.4
• /
• pp.139-149
• /
• 2003

In this paper a hardware design of area-efficient hash processor which implements MD5 algorithm using hardware sharing and carry-save addition schemes is described. To reduce area, the processor adopts hardware sharing scheme in which 1 step operation is divided into 2 substeps and then each substep is executed using the same hardware. Also to increase clock frequency, three serial additions of substep operation are transformed into two carry-save additions and one carry propagation addition. The MD5 hash processor is designed using 0.25 $\mu\textrm{m}$CMOS technology and consists of about 13,000 gates. From timing simulation results, the designed MD5 hash processor has 465 Mbps hash rates for 512-bit input message data under 120 MHz operating frequency.

• Bulletin of the Korean Chemical Society
• /
• v.30 no.5
• /
• pp.1139-1142
• /
• 2009

Refinement of NMR structures by molecular dynamics (MD) simulations with a solvent model has improved the structural quality. In this study, we applied MD refinement with the generalized Born (GB) implicit solvent model to protein structure determined under membrane-like environments. Despite popularity of the GB model, its applications to the refinement of NMR structures of hydrophobic proteins, in which detergents or organic solvents enclose proteins, are limited, and there is little information on the use of another GB parameter for these cases. We carried out MD refinement of crambin NMR structure in dodecylphosphocholine (DPC) micelles (Ahn et al., J. Am. Chem. Soc. 2006, 128, 4398-4404) with GB/Surface area model and two different surface tension coefficients, one for aquatic and the other for hydrophobic conditions. Our data show that, of two structures by MD refinement with GB model, the one refined with the parameter to consider hydrophobic condition had the better qualities in terms of precision and solvent accessibility.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.137-145
• /
• 2014

We present a design sensitivity analysis(DSA) method for multiscale problems based on bridging scale decomposition. In this paper, we utilize a bridging scale method for the coupled system analysis. Since the analysis of full MD systems requires huge amount of computational costs, a coupled system of MD-level and continuum-level simulation is usually preferred. The information exchange between the MD and continuum levels is taken place at the MD-continuum boundary. In the bridging scale method, a generalized Langevin equation(GLE) is introduced for the reduced MD system and the GLE force using a time history kernel is applied at the boundary atoms in the MD system. Therefore, we can separately analyze the MD and continuum level simulations, which can accelerate the computing process. Once the simulation of coupled problems is successful, the need for the DSA is naturally arising for the optimization of macro-scale design, where the macro scale performance of the system is maximized considering the micro scale effects. The finite difference sensitivity is impractical for the gradient based optimization of large scale problems due to the restriction of computing costs but the analytical sensitivity for the coupled system is always accurate. In this study, we derive the analytical design sensitivity to verify the accuracy and applicability to the design optimization of the coupled system.