• Bulletin of the Korean Chemical Society
• /
• 제35권12호
• /
• pp.3527-3531
• /
• 2014

In this paper, we report thermodynamic and transport properties (diffusion coefficient, viscosity, and thermal conductivity) of diatomic gases ($H_2$, $N_2$, $O_2$, and $Cl_2$) at 273.15 K and 1.00 atm by performing molecular dynamics simulations using Lennard-Jones intermolecular potential and modified Green-Kubo formulas. The results of self-diffusion coefficients of diatomic gases obtained from velocity auto-correlation functions by Green-Kubo relation are in good agreement with those obtained from mean square displacements by Einstein relation. While the results for viscosities of diatomic gases obtained from stress auto-correlation functions underestimate the experimental results, those for thermal conductivities obtained from heat flux auto-correlation functions overestimate the experimental data except $H_2$.

• Korean Membrane Journal
• /
• 제6권1호
• /
• pp.30-36
• /
• 2004

The gas permeation performance of commercially available polyetherimide (Ultem$\^$/) is simulated by means of molecular dynamics methods. By the observation of trajectory, long distance hopping of gas molecules is needed to transverse from top to bottom of membrane. Two possibilities mechanism of diffusion phenomena through glassy polymers can be issued. Diffusion coefficients were calculated by Einstein relation equation. In solubility simulation, the value of the constants C'$\_$H/ and b for O$_2$ at 300 K were calculated. The diffusion and solubility coefficient of He for PEI were simulated in this simulation work. the permeability coefficient is 9.88 Barrer. This value is closed to experimental value of 9.4 Barrer.

• Bulletin of the Korean Chemical Society
• /
• 제23권11호
• /
• pp.1595-1603
• /
• 2002

In this paper we have presented the results of diffusion behavior of model systems for eight liquid n-alkanes ($C_{12}$-$C_{44}$) in a canonical (NVT) ensemble at several temperatures using molecular dynamics simulations. For these n-alkanes of small chain length n, the chains are clearly <$R_{ee}^2$>/6<$R_g^2$>>1 and non-Gaussian. This result implies that the liquid n-alkanes over the whole temperatures considered are far away from the Rouse regime, though the ratio becomes close to the unity as n increases. Calculated self-diffusion constants $D_{self}$ are comparable with experimental results and the Arrhenius plot of self-diffusion constants versus inverse temperature shows a different temperature dependence of diffusion on the chain length. The global rotational motion of n-alkanes is examined by characterizing the orientation relaxation of the end-to-end vector and it is found that the ratio ${\tau}1/{\tau}2$ is less than 3, the value expected for a isotropically diffusive rotational process. The friction constants ${\xi}$of the whole molecules of n-alkanes are calculated directly from the force auto-correlation (FAC) functions and compared with the monomeric friction constants ${\xi}_D$ extracted from $D_{self}$. Both the friction constants give a correct qualitative trends: decrease with increasing temperature and increase with increasing chain length. The friction constant calculated from the FAC's decreases very slowly with increasing temperature, while the monomeric friction constant varies rapidly with temperature. By considering the orientation relaxation of local vectors and diffusion of each site, it is found that rotational and translational diffusions of the ends are faster than those of the center.

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

• Bulletin of the Korean Chemical Society
• /
• 제4권5호
• /
• pp.223-227
• /
• 1983

A molecular dynamic computer experiment was performed on a system of 108 particles composed of a single polymer chain and solvent molecules. The state considered was in the immediate neighborhood of the triple point of the system. The polymer itself is an analog of a freely jointed chain. The Lennard-Jones potential was used to represent the interactions between all particles except for that between the chain elements forming a bond in the polymer chain, for which the interaction was expressed by a harmonic potential. The self-diffusion coefficient and velocity autocorrelation function (VACF) of a polymer were calculated at various chain lengths $N_p$, and various interaction strengths between solvent molecules and a polymer chain element. For self-diffusion coefficients D, the Einstein relation holds good; as chain length $N_p$ increases the D value decreases, and D also decreases as ${\varepsilon}_{cs}$ (the interaction parameter between the chain element and solvent molecules) increases. The relaxation time of velocity autocorrelation decreases as ${\varepsilon}_{cs}$ increases, and it is constant for various chain lengths. The diffusion coefficients in various conditions reveal that our systems are in a free draining limit as is well known from the behavior of low molecular weight polymers, this also agrees with the Kirkwood-Riesman theory.

• Bulletin of the Korean Chemical Society
• /
• 제32권4호
• /
• pp.1336-1340
• /
• 2011

Molecular dynamics simulations have been performed to investigate the transport properties of self-diffusion coefficients in the penetrable-sphere model system. The resulting simulation data for the product of the packing fraction and the self-diffusion coefficient exhibit a transition from an increasing function of density in lower repulsive systems, where the soft-type collisions are dominant, to a decreasing function in higher repulsive systems, where most particle collisions are the hard-type reflections due to the low-penetrability effects. A modified Enskog-like equation implemented by the effective packing fraction with the mean-field energy correction is also proposed, and this heuristic approximation yields a reasonably good result even in systems of high densities and high repulsive energy barriers.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2014년도 제46회 동계 정기학술대회 초록집
• /
• pp.115.1-115.1
• /
• 2014

Atom-thick 2-dimensional materials such as graphene, h-BN and MoS2 hold substantial potential for applications in future molecular-scale integrated electronics, transparent conducting membranes, nanocomposites, etc. From a fundamental point of view, 2-dim crystal-solid substrates can also serve as a unique system to study various physicochemical phenomena occurring at low dimensions or interfaces. In this talk, I will present our recent Raman spectroscopy studies on the surface science problems of graphene: interfacial charge transfer, molecular diffusion in confined space and structural deformation.

• Bulletin of the Korean Chemical Society
• /
• 제21권1호
• /
• pp.93-96
• /
• 2000

A diffusion-influenced pseudo-first order reversible reaction A + B ⇔C + B is investigated by the molecular dynamics (MD) simulation method. Theoretical finding that the temporal evolution of reactants [conditional probabilities] in the reversible system can be expressed by the irreversible survival probability with an effective rate parameter is confirmed even in the presence of solvent particles. We carry out molecular dynamics simulations for both the irreversible and the reversible cases to evaluate the survival and the conditional probabilities for each cases. When the resultant irreversible survival probability is inserted into the proposed relation, the conditional probabilities given by the simulation are exactly reproduced.

• Bulletin of the Korean Chemical Society
• /
• 제14권1호
• /
• pp.110-112
• /
• 1993

Molecular reorientation of oblate symmetric top molecules in the presence of internal rotation is investigated and an analytic expression for the overall reorientational correlation time is obtained. The overall reorientation of the symmetric top is treated by the anisotropic rotational diffusion and the internal rotation is analyzed by employing a model which describes jumps between several discrete states with different lifetimes. The lifetimes thus obtained can be compared with the internal angular momentum correlation time which appears when the internal rotation is treated by a modified extended rotational diffusion model.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.204-207
• /
• 2008

Stochastic molecular dynamics simulation is a variation of standard molecular dynamics simulation that basically omits water molecules. The omission of water molecules, occupying a majority of space, enables flow simulation at microscale. This study reports our stochastic molecular dynamics simulation of particles diffusing in rectangular microchannels. We interestingly found that diffusion patterns in channels with a very small aspect ratio differ by dimensions. We will also discuss the future direction of our research toward a more realistic simulation of micromixing.