• 한국광물학회지
• /
• 제26권2호
• /
• pp.119-127
• /
• 2013

핵폐기물을 고화시키는 재료로 사용하는 붕규산염(borosilicate) 유리의 용해는 지층 처분장에 처리된 고준위 방사성 폐기물의 생태계 유출을 결정할 수 있는 중요한 화학반응이다. 습식 실험에서 유리의 용해속도(dissolution rate)는 유리 화학조성에 의해 크게 좌우되는 것이 관찰된다. 유리의 bulk 구조를 규명한 분광분석 실험에 의하면 유리의 화학조성과 분자수준(molecular-level) 구조(예: $SiO_4$ 사면체의 연결구조와 B 원소의 배위구조) 사이의 상관관계가 존재한다. 따라서 화학조성에 따른 유리 용해도의 차이는 조성에 따른 bulk 내부구조의 변화로 이해되어 왔다. 그런데 유리 표면은 수용액과 계면을 이루면서 용해 과정에서 가장 직접적으로 반응하는 부분이기 때문에, 화학조성에 따른 표면구조 변화에 대한 지식 또한 필요하다. 본 논문에서는 분자 동역학(molecular dynamics, MD) 시뮬레이션을 사용하여 4가지의 다른 화학조성을 가지는 소듐붕규산염 유리($xNa_2O{\cdot}B_2O_3{\cdot}ySiO_2$ 화학조성)에 대하여 bulk 구조와 실험으로 얻기 어려운 표면(surface) 구조를 연구하였다. MD 시뮬레이션은 유리 표면의 화학조성과 분자수준 구조가 bulk의 것과 매우 상이한 결과를 보여준다. 본 연구의 MD 시뮬레이션 결과는 화학조성에 따른 유리 용해도(특히 초기 용해과정)는 bulk 구조의 변화보다 유리 표면구조의 변화에 의해 크게 좌우될 수 있다는 표면구조에 대한 이해의 중요성을 역설한다.

• Coupled systems mechanics
• /
• 제6권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.

• Bulletin of the Korean Chemical Society
• /
• 제12권5호
• /
• pp.490-494
• /
• 1991

The ${\beta} to {\alpha}$ phase transition in silver iodide is studied with the (N, V, E) and (N, P, T) molecular dynamics (MD) method. In experiments, the phase transition temperature is 420 K. Upon heating of ${\beta}$ form, the iodine ions undergo hcp to bcc transformation and silver ions become mobile. MD simulations for the ${\beta}$ and ${\alpha}$ phases are carried out at several temperatures and the radial distribution functions (rdf) are obtained at those temperatures in the (N, V, E) ensemble. But the phase transition is not found in our calculation. Next the phase transition is studied with the (N, P, T) MD and we find some evidences of phase transition. At 3 Kbars and 2 Kbars the phase transition temperatu re is about 300 K. For 3.55 Kbars, the phase transition is higher (420 K) than the low pressure case. The phase transition temperature is somewhat dependent on the pressure in our calculations.

• Nuclear Engineering and Technology
• /
• 제41권1호
• /
• pp.11-20
• /
• 2009

Radiation effects on materials are inherently multiscale phenomena in view of the fact that various processes spanning a broad range of time and length scales are involved. A multiscale modeling approach to embrittlement of pressure vessel steels is presented here. The approach includes an investigation of the mechanisms of defect accumulation, microstructure evolution and the corresponding effects on mechanical properties. An understanding of these phenomena is required to predict the behavior of structural materials under irradiation. We used molecular dynamics (MD) simulations at an atomic scale to study the evolution of high-energy displacement cascade reactions. The MD simulations yield quantitative information on primary damage. Using a database of displacement cascades generated by the MD simulations, we can estimate the accumulation of defects over diffusional length and time scales by applying kinetic Monte Carlo simulations. The evolution of the local microstructure under irradiation is responsible for changes in the physical and mechanical properties of materials. Mechanical property changes in irradiated materials are modeled by dislocation dynamics simulations, which simulate a collective motion of dislocations that interact with the defects. In this paper, we present a multi scale modeling methodology that describes reactor pressure vessel embrittlement in a light water reactor environment.

• EDISON SW 활용 경진대회 논문집
• /
• 제3회(2014년)
• /
• pp.169-182
• /
• 2014

이 연구에서는 Lennard-Jones (LJ) particle을 이용하여 상분리 현상을 이해하기 위한 컴퓨터 시뮬레이션 연구를 수행하였다. 초기에 균일하게 분포되어 있는 LJ 입자들을 시뮬레이션 하면 상대적으로 dense phase와 dilute phase로 상분리 현상이 일어나게 된다. 상분리 현상의 첫 번째 단계를 핵 생성 (nucleation) 이라고 한다. 본 연구에서는 Brownian Dynamics (BD) Simulation과 Molecular Dynamics (MD) Simulation을 이용하여 상평형 그림을 구하고 초기에 일어나는 LJ 입자들의 nucleation rates를 구하였다.

• 한국분말재료학회지
• /
• 제22권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.

• 한국진공학회지
• /
• 제4권4호
• /
• pp.367-372
• /
• 1995

We investigate the Au(001) surface reconstruction, numerically, by Molecular Dynamics (MD) simulation. We find that the top-most layer of Au(001) surface is reconstructed to a contracted hexagonal face, and relaxed about 0.05$\AA$ upward at room temperature. The contraction ratio with respect to a unreconstructed Au(111) surface is about 3.5%. The hexagonal layer is slightly distorted and buckled. The surface corrugation is found to be about 0.28$\AA$ on average. In our earlier work we have predicted the in-plane orientation of the reconsturcted layer to be either $0^{\circ}$ or $0.7^{\circ}$ depending on the size of the cluster. However, we find only $0.0^{\circ}$ in this simulation because the size of the cluster correspoding to the $0.7^{\circ}$ orientation is larger than the current limitation of MD simulation. These findings are in good agreement with experimental results.

• Bulletin of the Korean Chemical Society
• /
• 제18권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.

• Bulletin of the Korean Chemical Society
• /
• 제19권10호
• /
• pp.1068-1072
• /
• 1998

We present results of molecular dynamic (MD) simulations for the segmental motion of liquid n-butane as the base case for a consistent study for conformational transition from one rotational isomeric state to another in long chains of liquid n-alkanes. The behavior of the hazard plots for n-butane obtained from our MD simulations are compared with that for n-butane of Brownian dynamics study. The MD results for the conformational transition of n-butane by a Poisson process form the total first passage times are different from those from the separate t-g and g-t first passage times. This poor agreement is probably due to the failure of the detailed balance between the fractions of trans and gauche. The enhancement of the transitions t-g and g-t at short time regions are also discussed.

• 통합자연과학논문집
• /
• 제13권4호
• /
• pp.170-180
• /
• 2020

Abnormal regulation, hyperphosphorylation, and aggregation of the tau protein are the hallmark of several types of dementia, including Alzheimer's Disease. Increased activity of Glycogen Synthase Kinase-3β (GSK-3β) in the Central Nervous System (CNS), increased the tau hyperphosphorylation and caused the neurofibrillary tangles (NFTs) formation in the brain cells. Over the last two decades, numerous adenosine triphosphate (ATP) competitive inhibitors have been discovered that show inhibitory activity against GSK-3β. But these compounds exhibited off-target effects which motivated researchers to find new GSK-3β inhibitors. In the present study, we have collected the dataset of 31 C-Glycosylflavones derivatives that showed inhibitory activity against GSK-3β. Among the dataset, the most active compound was docked with the GSK-3β and molecular dynamics (MD) simulation was performed for 50 ns. Based on the 50 ns MD pose of the most active compound, the other dataset compounds were sketched, minimized, and aligned. The 3D-QSAR based Comparative Molecular Field Analysis (CoMFA) model was developed, which showed a reasonable value of q2=0.664 and r2=0.920. The contour maps generated based on the CoMFA model elaborated on the favorable substitutions at the R2 position. This study could assist in the future development of new GSK-3β inhibitors.