• EDISON SW 활용 경진대회 논문집
• /
• 제4회(2015년)
• /
• pp.91-98
• /
• 2015

Formamide의 중성가수분해 mechanism은 QM/MM (quantum mecahnics/molecular mechanics) molecular dynamics simulations 및 CPMD과 같은 방법으로 연구되어왔다. 본 연구에서는. Umbrella sampling을 이용한 QM/MM-MD simulation을 사용하여 4가지 반응의 free energy surface를 도출해냈다. 전체적으로, 가장 선호되는 메커니즘은 two step으로 구성된 water assisted stepwise mechanism이었으며 모든 mechanism은 ab-initio calculation과 QM/MM-MD simulation이 수행되었다. water assisted stepwise mechanism을 살펴보면, 첫 번째 step에서 formamide의 carbonyl group이 hydrate되면서 gem-diol intermediate를 형성한다. 다음 step에서, intermediate의 hydroxyl group으로부터 amino group으로 water-assisted proton transfer이 일어난다. 두 반응 모두에서 물이 proton transfer를 직접적으로 도와주는 것을 관찰할 수 있었다. 특히, ab-initio calculation과는 다르게 QM/MM-MD에서는 gem-diol intermediate가 안정화되는 것으로 solvent effect를 잘 보여준다.

• 한국재료학회지
• /
• 제20권9호
• /
• pp.478-482
• /
• 2010

A series of molecular dynamic (MD), finite element (FE) and ab initio simulations are carried out to establish suitable modeling schemes for the continuum-based analysis of aluminum matrix nanocomposites reinforced with carbon nanotubes (CNTs). From a comparison of the MD with FE models and inferences based on bond structures and electron distributions, we propose that the effective thickness of a CNT wall for its continuum representation should be related to the graphitic inter-planar spacing of 3.4${\AA}$. We also show that shell element representation of a CNT structure in the FE models properly simulated the carbon-carbon covalent bonding and long-range interactions in terms of the load-displacement behaviors. Estimation of the effective interfacial elastic properties by ab initio simulations showed that the in-plane interfacial bond strength is negligibly weaker than the normal counterpart due to the nature of the weak secondary bonding at the CNT-Al interface. Therefore, we suggest that a third-phase solid element representation of the CNT-Al interface in nanocomposites is not physically meaningful and that spring or bar element representation of the weak interfacial bonding would be more appropriate as in the cases of polymer matrix counterparts. The possibility of treating the interface as a simply contacted phase boundary is also discussed.

• Bulletin of the Korean Chemical Society
• /
• 제11권6호
• /
• pp.511-515
• /
• 1990

Ab initio calculations with MP3 and CISD method were performed for the FCO and $FCO^+$. The equilibrium geometry of FCO shows the reasonable agreement with experimental values. $FCO^+$ has a linear geometry with $R_{CF}=1.213\;and\;R_{CO}=1.118{\AA}$. The quadratic force constants of $FCO^+$ are 23.21 md ${\AA}^{-1}$ for CO stretch and 12.38 md ${\AA}^{-1}$ for CF stretch. The cubic force constants and the other molecular constants are also calculated.

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

Peptide bond hydrolysis는 세포 내외의 생화학반응에 있어서 핵심이다. 하지만 amide Hydrolysis Mechanism은 아직 명확하게 규명되지 않았다. pH가 중성인 물에서의 비 촉매 가수분해가 발생하는 몇몇 실험적 증거가 있지만, 해당 반응 매커니즘은 4 가지(non-assisted concerted, non-assisted step-wise, assisted concerted, assisted step-wise)로 여전히 논란이 있다. 이번 연구에서는, Formamide의 가능한 Hydrolysis Mechanism을 자세히 연구해보고자 한다. 먼저, Ab-initio 계산을 통해 4가지 반응 메커니즘의 다시 한번 확인하고, quantum chemical calculations과 quantum mechanical molecular dynamic이 결합된 (QMMD) simulation을 통하여 water solvent에서의 반응 메커니즘의 에너지관계를 규명하였다. 결론적으로 아직 계산이 끝나지 않은 supported concerted mechanism을 제외한 모든 계산에서 non-supported, supported 두 system 모두에서 step-wise가 일어나기 쉬웠고, non-supported 보다 supported mechanism이 선호됨을 보였다. Intermediate인 amino-gem-diol의 수용액 상에서 안정화 또한 나타났다. 이는 Ab-initio 계산만 통해서는 정확하게 산출할 수 없는 엔트로피의 영향을 잘 보여준다.

• 대한화학회지
• /
• 제56권6호
• /
• pp.669-672
• /
• 2012

Classical molecular dynamics (MD) simulations using a scaled OSS2 potential originally derived from ab initio calculations are used to study the radial distribution functions of water. The original OSS2 water potential is shown to represent a glassy or an ice at ambient temperature, but the diffusion coefficient increases on increasing the temperature of the system or decreasing the density. This suggests scaling the OSS2 potential. The O-O, O-H, and H-H radial distribution functions and the corresponding coordination numbers for the scaled OSS2 potential, obtained by MD simulation, are in good agreement with the experiment results and calculations for the SPC/E water potential over a range of temperatures.

• Bulletin of the Korean Chemical Society
• /
• 제20권8호
• /
• pp.948-952
• /
• 1999

Electrophilic substitutions on β-position of 2,5-dimethyl pyrrole have been investigated theoretically. The electron donating methyl groups enrich electron densities on C-3, C-4 positions and π* interactions with methyl groups substituted on C-2 and C-5 positions pushed up the HOMO level of the pyrroles consequently induce rapid substitutions on C-3, C-4 sites. Substitution of phenylsulfonyl group on nitrogen stabilized LUMO levels through weak π bonding interactions. Unexpected deoxidation reaction underwent on the sulfonyl group substituted at C-3 position. The structures were solved by X-ray crystallography. Meanwhile, gas phase HF/6-31G* and density functional method (B3LYP/6-31G*) calculations gave favorable energies for 1-phenylsulfinyl pyrrole (6) over 3-phenylsulfinyl pyrrole (5) by 3.6-4.7 kcal/mol which is contrary to the experimental result. However the methods involve the effects of molecular polarizability and solvent, molecular dynamics (MD) and ab-initio self consistent reaction field (SCRF) calculations showed same trend as experiments. According to MD calculations, compound 5 is more stable than compound 6 by 4.15 kcal/mol and the SCRF, HF/6-31G* calculations gave more stable energy value for structure 5 than 6 by 0.03 kcal/mol.

• 한국세라믹학회지
• /
• 제53권2호
• /
• pp.184-193
• /
• 2016

This review deals with computational treatments of subatomic levels of matter based on density functional theory (DFT), and tries to identify several current trends, which are largely consequences of the ever-increasing power of computers, which has substantially extended the performance of conventional DFT beyond its original scope. This review mainly focuses on the conceptual outline, rather than on lines of equations, highlighting several examples of calculations for each topic. It should be noted that these issues are hardly new to leading groups in the field, but certainly are for materials people in general. It should also be noted that the on-going efforts will continue and lead to a larger system size, a longer time scale, a higher accuracy, and a better efficiency of calculation for years to come.

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

Half-metallic nanostructure is highly applicable in the field of Spintronics and electronic device technology. We examine the electronic properties of a ferrocene-based nanowire as a possible candidate for a half-metallic nanostructure using VASP and SIESTA. Ferrocene-based nanowire shows high stability in both binding energy simulation and molecular dynamics (MD) simulation. The density of states (DOS) and the projected DOS of the ferrocene-based nanowire indicate that one-dimensional clustering of ferrocene molecules can be explained because of p-d orbital hybridization between iron and carbon. Half-metallic property and energy dispersion at the Fermi level due to one-dimensional structure is also observed from the DOS results.

• EDISON SW 활용 경진대회 논문집
• /
• 제5회(2016년)
• /
• pp.120-125
• /
• 2016

용매화 된 이온들간의 상호작용은 화학을 비롯한 여러 분야에서 중요하게 다뤄지고 있다. 계산화학에 있어 다양한 용매에 용매화 된 이온들간의 상호작용은 분자동력학 시뮬레이션(MD)을 이용한 PMF 계산을 통해 많이 연구되어 왔는데, 본 그룹에서는 물을 표현할 수 있는 다양한 모델과 Salt solution 하에 Na+ 이온과 Cl- 이온 사이 상호작용을 PMF 계산을 통해 알아보았다. 첫 번째로 Polarizable water model인 Tip4p-fq를 사용한 PMF 계산 결과는 아직 진행 중에 있으며 결과적으로는 기본적으로 CHARMM water force field에서 제공해주는 Tip3p 모델 및 양자계산 방법 중 하나인 Ab initio를 이용한 결과와 비교, 분석 할 예정이다. 두 번째로 1M/2M/5M Salt solution의 서로 다른 농도를 사용하여 시뮬레이션 한 결과 1M과 2M에서는 물 속에서 시뮬레이션 한 결과와 큰 차이가 없었으나 농도가 짙은 5M 용매 내에선 $Na^+$이온과 $Cl^-$이온 사이 거리가 멀어질수록 주변 Salt에 영향을 받는 것을 확인 할 수 있었다. 결론적으로 본 연구를 진행함으로써 우리는 서로 다른 전하를 갖는 이온 사이 상호작용에 대한 이해를 더 높일 수 있었다.

• Journal of Photoscience
• /
• 제9권2호
• /
• pp.51-54
• /
• 2002

We propose a novel mechanism (Twist Sharing Mechanism) for the cis-trans photoisomerization of rhodopsin, based on the molecular dynamics (MD) simulation study. New things devised in our simulations are (1) the adoption of Mt. Fuji potentials in the excited state for twisting of the three bonds C9=C10, C11=C12 and C13=14 which are modeled using the detailed ab initio quantum chemical calculations and (2) to use the rhodopsin structure which was resolved recently by the X-ray crystallographic study. As a result, we found the followings: Due to the intramolecular steric hindrance between 20-methyl and 10-H in the retinal chromophore, the C12-C13 and C10-C11 bonds are considerably twisted counterclockwise in rhodopsin, allowing only counterclockwise rotation of the C11 =C12 in the excited state. The movement of 19-methyl in rhodopsin is blocked by the surrounding three amino acids, Thr 118, Met 207 and Tyr 268, prohibiting the rotation of C9=C10. As a result only all-trans form of the chromophore is obtainable as a photoproduct. At the 90$^{\circ}$ twisting of C11=C12 in the course of photoisomerization, twisting energies of the other bonds amount to about 20 kcal/mol. If the transition state for the thermal isomerization is assumed to be similar to this structure, the activation energy for the thermal isomerization around C11=C12'in rhodopsin is elevated by about 20 kcal/mol and the thermal isomerization rate is decelerated by 10$\^$-14/ times than that of the retinal chromophore in solution, protecting photosignal from the thermal noise.