Molecular Mobility in Cyclic Hydrocarbons: A Simulation Study

Application of Reverse Nonequilibrium Molecular Dynamics to the Calculation of the Mutual Diffusion Coefficient of Alkane Mixtures

In a recent publication, a reverse nonequilibrium molecular dynamics (RNEMD) method was presented for computing the mutual diffusion coefficient of liquid mixtures. A concentration gradient and a subsequent mass flux are induced in the system by suitably exchanging molecules in different regions. The algorithm has been successfully tested on Lennard-Jones mixtures and molecular fluid mixtures with molecules having the same number of particles. In this work, a modification is made to the RNEMD method to determine the mutual diffusion coefficient of binary liquid mixtures with molecules having different sizes and masses. To migrate molecules of a different type, the splitting method has been used in this work. Investigation of the resulting steady-state mass fraction profile allows the evaluation of the mutual diffusion coefficient. For validation, the mutual diffusion coefficients of ethane-propane and ethane-pentane liquid mixtures at different compositions and temperatures have been obtained using this method. The mutual diffusion coefficients obtained from the RNEMD simulations are within the error bars of values obtained by equilibrium molecular dynamics for the identical model and conditions. The excess energy released due to the exchange of molecules is efficiently removed by strongly coupling a local thermostat in the region around the insertion point. There is no heating of the analysis region.

Atomistic picture of isothermal volume relaxation behavior of atactic polystyrene glass provided by a molecular dynamics simulation

Using an atomistic molecular dynamics (MD) simulation, we study the volume relaxation behavior of atactic polystyrene and attempt to correlate this macroscopic behavior with certain microscopic aspects. To this end, the gyration radius, the dimensionless relative shape anisotropy, the mean-squared displacement, and the non-Gaussian parameter are examined simultaneously. Our result shows that the structures characterized at different length scales change in a self-similar way and these changes are intimately correlated to the translational mobility of atoms. The initial incubation of structural changes at the different scales originates from the restriction of mobility due to the cage effect. The applicability of the MD simulation to the investigation of the bulk properties is discussed.

The Soret effect in dilute polymer solutions: influence of chain length, chain stiffness, and solvent quality

Thermal diffusion in dilute polymer solutions is studied by reverse nonequilibrium molecular dynamics. The polymers are represented by a generic bead-spring model. The influence of the solvent quality on the Soret coefficient is investigated. At constant temperature and monomer fraction, a better solvent quality causes a higher affinity for the polymer to the cold region. This may even go to thermal-diffusion-induced phase separation. The sign of the Soret coefficient changes in a symmetric nonideal binary Lennard-Jones solution when the solvent quality switches from good to poor. The known independence of the thermal diffusion coefficients of the molecular weight is reproduced for three groups of polymers with different chain stiffnesses. The thermal diffusion coefficients reach constant values at chain lengths of around two to three times the persistence length. Moreover, rigid polymers have higher Soret coefficients and thermal diffusion coefficients than more flexible polymers.

Reverse nonequilibrium molecular-dynamics calculation of the Soret coefficient in liquid benzene/cyclohexane mixtures

Reverse nonequilibrium molecular dynamics is the method applied here for the investigation of thermal diffusion in realistic molecular fluids. The Soret coefficients of benzene/cyclohexane mixtures are calculated using an all-atom model. The autocorrelation functions indicate that the mole fraction gradient converges much slower than the temperature gradient. Compared to experimental data, the results show the same tendency of the Soret coefficient variation versus the mole fraction. Although a systematic error exists for the magnitude of the Soret coefficient, a meanwhile systematic error for both the mutual diffusion and thermal diffusion coefficients provides some explanation of it; and the calculation with different force field parameters indicates a possibility to annihilate the systematic error. The influences of algorithm variables such as cutoff lengths and perturbation intensities are tested. Furthermore the temperature dependence of the Soret effect is observed, yielding the same trend as previous studies.

Molecular Dynamics of a Polymer in Mixed Solvent: Atactic Polystyrene in a Mixture of Cyclohexane and N,N-Dimethylformamide

Local structural and dynamic properties of atactic polystyrene in a mixed solvent of cyclohexane (CH) and N,N-dimethylformamide (DMF) have been investigated using molecular dynamics simulations. We measure local conformations in the polymer and classify them by distance and angle distribution histograms. End-to-end distances and structure factors are employed to describe the static structure of polystyrene chains. Polystyrene concentration, including 1.6%, 4.8%, and 14% (by weight), and solution temperatures of 300, 330, or 360 K are used to elucidate the concentration and temperature dependencies of the solvation by the two solvents. Both solvent molecules preferentially approach the phenyl rings. At lower temperatures, polystyrene dissolves more favorably in cyclohexane. With rising temperature DMF molecules approach more closely with the oxygen oriented toward the phenyl rings. Additionally, the global and segmental relaxation times of the chains decrease and the system becomes more homogeneous. The segmental and global dynamics of polystyrene show different concentration behaviors: the reorientation times of solvent molecules and segments of polystyrene increase with concentration while the global dynamics of polystyrene decelerates as the concentration is changed from 1.6% to 4.8% but accelerates when the concentration rises to 14%. We conclude that the change of concentration from 4.8% to 14% qualitatively marks the change from a dilute to a semidilute solution. The diffusion constants of the small molecules and corresponding activation energies have also been measured. Our simulation data are compared with available experimental results and we find a satisfactory agreement.

Thermal Conductivities of Molecular Liquids by Reverse Nonequilibrium Molecular Dynamics

The reverse nonequilibrium molecular dynamics method for thermal conductivities is adapted to the investigation of molecular fluids. The method generates a heat flux through the system by suitably exchanging velocities of particles located in different regions. From the resulting temperature gradient, the thermal conductivity is then calculated. Different variants of the algorithm and their combinations with other system parameters are tested: exchange of atomic velocities versus exchange of molecular center-of-mass velocities, different exchange frequencies, molecular models with bond constraints versus models with flexible bonds, united-atom versus all-atom models, and presence versus absence of a thermostat. To help establish the range of applicability, the algorithm is tested on different models of benzene, cyclohexane, water, and n-hexane. We find that the algorithm is robust and that the calculated thermal conductivities are insensitive to variations in its control parameters. The force field, in contrast, has a major influence on the value of the thermal conductivity. While calculated and experimental thermal conductivities fall into the same order of magnitude, in most cases the calculated values are systematically larger. United-atom force fields seem to do better than all-atom force fields, possibly because they remove high-frequency degrees of freedom from the simulation, which, in nature, are quantum-mechanical oscillators in their ground state and do not contribute to heat conduction.