• Macromolecular Research
• /
• 제12권2호
• /
• pp.178-188
• /
• 2004

We performed molecular dynamics (MD) simulations at 300 K on a series of poly(benzyl ether) (PBE) dendrimers having a different core functionalities. We used the rotational isomeric state Metropolis Monte Carlo (RMMC) method to construct the initial configuration in a periodic boundary cell (PBC) before the MD simulations were undertaken. To elucidate the effects that the structural features have on the chain dimension, the overall internal structure, and the morphology, we monitored the radii of gyration, R$\sub$g/ and the conformational changes during the simulations. The PBE dendrimers in a glassy state adopted less-extended structures when compared with the conformations obtained from the RMMC calculations. We found that R$\sub$g/ of the PBE dendrimer depends on the molecular weight, M, according to the relation, R$\sub$g/∼M$\^$0.22/. The radial distributions of the dendrimers were developed identically in the PBC, irrespective of the core functionality. A gradual decrease in radial density resulted from the fact that the terminal branch ends are distributed all over the molecule, except for the core region.

• 셀메드
• /
• 제6권4호
• /
• pp.24.1-24.4
• /
• 2016

The emergence of influenza virus and antigenic drift are potential cause of world-wide pandemic. There are some commercially available drugs in the market to treat influenza. During past decade, however, critical resistances have been raised for biological targets. Because of structural complexity and flexibility of target proteins, applying a computational modeling tool is very beneficial for developing alternative anti-influenza drugs. In this review, we introduced molecular dynamics (MD) simulations approach to reflect full conformational flexibility of proteins during molecular modeling works. Case studies of MD works were summarized for the drug discovery and drug resistance mechanism of anti-influenza pharmaceuticals.

• Bulletin of the Korean Chemical Society
• /
• 제12권3호
• /
• pp.309-315
• /
• 1991

The thermal transport coefficients-the self-diffusion coefficient, shear viscosity, and thermal conductivity-of liquid argon at 94.4 K and 1 atm are calculated by non-equilibrium molecular dynamics (NEMD) simulations of a Lennard-Jones potential and compared with those obtained from Green-Kubo relations using equilibrium molecular dynamics (EMD) simulations and with experimental data. The time-correlation functions-the velocity, pressure, and heat flux auto-correlation functions-of liquid argon obtained from the EMD simulations show well-behaved smooth curves which are not oscillating and decaying fast around 1.5 ps. The calculated self-diffusion coefficient from our NEMD simulation is found to be approximately 40% higher than the experimental result. The Lagrange extrapolated shear viscosity is in good agreement with the experimental result and the asymptotic formula of the calculated shear viscosities seems to be an exponential form rather than the square-root form predicted by other NEMD studies of shear viscosity. The agreement for thermal conductivity between the simulation results (NEMD and EMD) and the experimental result is within statistical error. In conclusion, through our NEMD and EMD simulations, the overall agreement is quite good, which means that the Green-Kubo relations and the NEMD algorithms of thermal transport coefficients for simple liquids are valid.

• Bulletin of the Korean Chemical Society
• /
• 제26권8호
• /
• pp.1203-1208
• /
• 2005

A computational study based on molecular dynamics (MD) simulations was performed in order to explain the difference in aqueous solubilities of two flavonoid/$\beta$-cyclodextrin ($\beta$-CD) complexes, hesperetin/$\beta$-CD and naringenin/$\beta$-CD. The aqueous solubility of each flavonoid/$\beta$-CD complex could be characterized by complexwater interaction not by flavonoid-CD interaction. The radial distribution of water around each inclusion complex elucidated the difference of an experimentally observed solubility of each flavonoid/$\beta$-CD complex. The analyzed results suggested that a bulky hydrophobic moiety (-$OCH_3$) of B-ring of hesperetin nearby primary rim of $\beta$-CD was responsible for lower aqueous solubility of the hesperetin/$\beta$-CD complex.

• Bulletin of the Korean Chemical Society
• /
• 제15권10호
• /
• pp.852-856
• /
• 1994

A series of aliphatic hydrocarbons, methane to hexane in the liquid state, are modeled with the molecular mechanical potential parameters treating all hydrogen degrees of freedom explicitly. Thermodynamic properties (heat capacities and heats of vaporization) are calculated from relatively short (20ps) molecular dynamics trajectories. The liquid state structures are also examined through various radial distribution functions. Molecular dynamics simulations reproduce experimentally measured properties within a few percent errors, thus indicate that the present set of all-hydrogen parameters is suitable for simulating macromolecular systems in bulk.

• 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.

• Bulletin of the Korean Chemical Society
• /
• 제23권3호
• /
• pp.441-446
• /
• 2002

We have presented the results of thermodynamic, structural and dynamic properties of liquid benzene, toluene, and p-xylene in canonical (NVT) ensemble at 293.15 K by molecular dynamics (MD) simulations. The molecular model adopted for these molecules is a combination of the rigid body treatment for the benzene ring and an atomistically detailed model for the methyl hydrogen atoms. The calculated pressures are too low in the NVT ensemble MD simulations. The various thermodynamic properties reflect that the intermolecular interactions become stronger as the number of methyl group attached into the benzene ring increases. The pronounced nearest neighbor peak in the center of mass g(r) of liquid benzene at 293.15 K, provides the interpretation that nearest neighbors tend to be perpendicular. Two self-diffusion coefficients of liquid benzene at 293.15 K calculated from MSD and VAC function are in excellent agreement with the experimental measures. The self-diffusion coefficients of liquid toluene also agree well with the experimental ones for toluene in benzene and for toluene in cyclohexane.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.127-136
• /
• 2013

In this study, atomistic simulations are performed to study the effect of Al solute on the behaviour of edge dislocation in Mg alloys. After the dissociation of an Mg basal edge dislocation into two Shockley partials using molecular mechanics, the interaction between the dislocation and Al solute at different temperatures is studied using molecular dynamics. It appears from the simulations that the critical shear stress increases with the Al solute concentration. Comparing with the solute effect at T = 0 K, however, the critical shear stress at a finite temperature is lower since the kinetic energy of the atoms can help the dislocation conquer the energy barriers created by the Al atoms. The velocity of the edge dislocation decreases as the Al concentration increases when the external shear stress is relatively small regardless of temperature. The Al concentration effect on the dislocation velocity is not significant at very high shear stress level when the solute concentration is below 4.0 at%. Drag coefficient B increases with the Al concentration when the stress to temperature ratio is below 0.3 MPa/K, although the effect is more significant at low temperatures.

• 한국분무공학회지
• /
• 제13권1호
• /
• pp.34-41
• /
• 2008

Molecular dynamic simulations have been carried out to study the effect of the nano-structure substrate and its temperature on cluster laminating. The interaction between substrate molecules and liquid molecules was modeled in the molecular scale and simulated by the molecular dynamics method in order to understand behaviors of the liquid cluster on nano-structure substrate. In the present model, the Lennard-Jones potential is applied to mono-atomic molecules of argon as liquid and platinum as nano-structure substrate to perform simulations of molecular dynamics. The effect of wettability on a substrate was investigated for the various beta of Lennard-Jones potential. The behavior of the liquid cluster and nano-structure substrate depends on interface wettability and function of molecules force, such as attraction and repulsion, in the collision progress. Furthermore, nano-structure substrate temperature and beta of Lennard-Jones potential have effect on the accumulation ratio. These results of simulation will be the foundation of coating application technology for micro fabrication manufacturing.

• Bulletin of the Korean Chemical Society
• /
• 제29권8호
• /
• pp.1505-1509
• /
• 2008

Comparative molecular simulations were performed to establish molecular interaction and inhibitory effect of flavonoid myricetin on formation of amyloid fibris. For computational comparison, the conformational stability of myricetin with amyloid $\beta$ -peptide (A$\beta$ ) and $\beta$ -amyloid fibrils (fA$\beta$) were traced with multiple molecular dynamics simulations (MD) using the CHARMM program from Monte Carlo docked structures. Simulations showed that the inhibition by myricetin involves binding of the flavonoid to fA$\beta$ rather than A$\beta$ . Even in MD simulations over 5 ns at 300 K, myricetin/fA$\beta$ complex remained stable in compact conformation for multiple trajectories. In contrast, myricetin/A$\beta$ complex mostly turned into the dissociated conformation during the MD simulations at 300 K. These multiple MD simulations provide a theoretical basis for the higher inhibitory effect of myricetin on fibrillogenesis of fA$\beta$ relative to A$\beta$ . Significant binding between myricetin and fA$\beta$ observed from the computational simulations clearly reflects the previous experimental results in which only fA$\beta$ had bound to the myricetin molecules.