• Bulletin of the Korean Chemical Society
• /
• v.31 no.8
• /
• pp.2402-2405
• /
• 2010

• Bulletin of the Korean Chemical Society
• /
• v.31 no.4
• /
• pp.959-964
• /
• 2010

In this work, equilibrium molecular dynamics (MD) simulations in a microcanonical ensemble are performed to evaluate the friction coefficient of a Brownian particle (BP) in a Lennard-Jones (LJ) solvent. The friction coefficients are determined from the time dependent friction coefficients and the momentum autocorrelation functions of the BP with its infinite mass at various ratios of LJ size parameters of the BP and solvent, ${\sigma}_B/{\sigma}_s$. The determination of the friction coefficients from the decay rates of the momentum autocorrelation functions and from the slopes of the time dependent friction coefficients is difficult due to the fast decay rates of the correlation functions in the momentum-conserved MD simulation and due to the scaling of the slope as 1/N (N: the number of the solvent particle), respectively. On the other hand, the friction coefficient can be determined correctly from the time dependent friction coefficient by measuring the extrapolation of its long time decay to t=0 and also from the decay rate of the momentum autocorrelation function, which is obtained by time integration of the time dependent friction coefficient. It is found that while the friction coefficient increases quadratically with the ratio of ${\sigma}_B/{\sigma}_s$ for all ${\sigma}_B$, for a given ${\sigma}_s$ the friction coefficient increases linearly with ${\sigma}_B$.

• Bulletin of the Korean Chemical Society
• /
• v.24 no.6
• /
• pp.771-776
• /
• 2003

The structural and dynamic properties of small n-alkane clusters embedded in a mesoscopic solvent are investigated. The solvent interactions are taken into account through a multi-particle collision operator that conserves mass, momentum and energy and the solvent dynamics is updated at discrete time intervals. The cluster molecules interact among themselves and with the solvent molecules through intermolecular forces. The properties of n-heptane and n-decane clusters interacting with the mesoscopic solvent molecules through repulsive Lennard-Jones interactions are studied as a function of the number of the mesoscopic solvent molecules. Modifications of both the cluster and solvent structure as a result of cluster-solvent interactions are considered. The cluster-solvent interactions also affect the dynamics of the small n-alkane clusters.

• Bulletin of the Korean Chemical Society
• /
• v.4 no.2
• /
• pp.87-91
• /
• 1983

Molecular dynamic method has been applied to a single polymer chain immersed in a solvent. The interactions for the pairs, of two solvent molecules (SS), of a chain element and a solvent molecules (CS), and of two non-neighbor chain elements (CC) are given by the Lennard-Jones potential, and the interaction between two bonded chain elements is given by a harmonic potential. We changed the CS interaction parameter ${\varepsilon}_{CS}$ to 0.5, 1.0 and 2.0 times of the SS interaction ${\varepsilon}_{SS}$. We calculated the pair correlation functions for the SS, CS, and CC pairs, end-to-end distance and radius of gyration with the varying ${\varepsilon}_{CS}$ parameters. The results showed that a phase separation occurs between the polymer and solvent in the 0.5 system where ${\varepsilon}_{CS}$ = 0.5 {\varepsilon}_{SS}$. The autocorrelation functions for end-to-end distance and radius of gyration were also calculated.

• Bulletin of the Korean Chemical Society
• /
• v.3 no.2
• /
• pp.44-49
• /
• 1982

Dynamic and equilibrium structures of a polymer chain immersed in solvent molecules have been investigated by a molecular dynamic method. The calculation employs the Lennard-Jones potential function to represent the interactions between two solvent molecules (SS) and between a constituent particle (monomer unit) of the polymer chain and a solvent molecule (CS) as well as between two non-nearest neighbor constituent particles of the polymer chain (CC), while the chemical bond for nearest neighbor constituent particles was chosen to follow a harmonic oscillator potential law. The correlation function for the SS, CS and CC pairs, the end-to-end distance square and the radius of gyration square were calculated by varying the chain length (= 5, 10, 15, 20). The computed end-to-end distance square and the radius of gyration square were found to be in a fairly good agreement with the corresponding results from the random-flight model. Unlike earlier works, the present simulation rsesult shows that the autocorrelation function of radius of gyration square decays slower than that of the end-to-end distance square.

• Bulletin of the Korean Chemical Society
• /
• v.4 no.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.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.18 no.5
• /
• pp.60-65
• /
• 2019

Since nanofluids (NFs), which are a mixture of a small amount of nanoparticles and a bulk liquid solvent, were first proposed by Stephen Choi at the Argonne National Lab in 1995, they have been considered for use in many technical studies of power cooling systems and their practical application due to their high thermal conductivity and heat transfer coefficients compared to conventional coolants. Although nanofluids are a well-known form of engineering fluid that show great promise for use in future cooling systems, their underlying physics as demonstrated in experiments remain unclear. One proven method of determining the heat transfer performance of nanofluids is measuring the concentration of nanoparticles in a mixture. However, it is experimentally inefficient to build testbeds to systematically observe particle distributions on a nanoscale. In this paper, we demonstrate the distribution of nanoparticles under a temperature gradient in a solution using molecular dynamics simulations. First, temperature profiles based on substrate temperature are introduced. Following this, the radial pair distribution functions of pairs of nanoparticles, solvents, and substrates are calculated. Finally, the distribution of nanoparticles in different heating regions is determined.

• Bulletin of the Korean Chemical Society
• /
• v.29 no.8
• /
• pp.1554-1560
• /
• 2008

The probe diffusion and friction constants of methyl yellow (MY) in liquid n-alkanes of increasing chain length were calculated by equilibrium molecular dynamics (MD) simulations at temperatures of 318, 418, 518 and 618 K. Lennard-Jones particles with masses of 225 and 114 g/mol are modeled for MY. We observed that the diffusion constant of the probe molecule follows a power law dependence on the molecular weight of nalkanes, DMY${\sim}M^{-\gamma}$ well. As the molecular weight of n-alkanes increases, the exponent $\gamma$ shows sharp transitions near n-dotriacontane ($C_{32}$) for the large probe molecule (MY2) at low temperatures of 318 and 418 K. For the small probe molecule (MY1) $D_{MY1}$ in $C_{12}$ to C80 at all the temperatures are always larger than Dself of n-alkanes and longer chain n-alkanes offer a reduced friction relative to the shorter chain n-alkanes, but this reduction in the microscopic friction for MY1 is not large enough to cause a transition in the power law exponent in the log-log plot of DMY1 vs M of n-alkane. For the large probe molecule (MY2) at high temperatures, the situation is very similar to that for MY1. At low temperatures and at low molecular weights of n-alkanes, $D_{MY2}$ are smaller than $D_{self}$ of n-alkanes due to the relatively large molecular size of MY2, and MY2 experiences the full shear viscosity of the medium. As the molecular weight of n-alkane increases, $D_{self}$ of n-alkanes decreases much faster than $D_{MY2}$ and at the higher molecular weights of n-alkane, MY2 diffuses faster than the solvent fluctuations. Therefore there is a large reduction of friction in longer chains compared to the shorter chains, which enhances the diffusion of MY2. The calculated friction constants of MY1 and MY2 in liquid n-alkanes supported these observations. We deem that this is the origin of the so-called“solventoligomer”transition.