Test Particle Insertion Method Research Articles

Measuring Oxygen Solubility in Amorphous and Semicrystalline Polyolefins Using Test Particle Insertion: A Comparative Study of Polyethylene and Isotactic Polypropylene.

The test particle insertion method is used to study the solubility of oxygen in two commonly used polymers: polyethylene (PE) and isotactic polypropylene (iPP). Amorphous samples for both polymers were prepared by means of Monte Carlo and molecular dynamics simulations, and the oxygen solubility was measured across different temperatures. The solubility-temperature dependence for iPP proved to be nonmonotonic due to the interplay between binding and reorganizational enthalpy, while for PE, it appeared to be monotonic based on the available data in the studied temperature range. A broad comparison to experiments and simulations is included. Further oxygen insertions were also performed in semicrystalline PE and iPP samples at ambient temperature, and the obtained values were compared to a linear relationship which correlates the solubility in the purely amorphous phase with the solubility in the crystalline phase. The solubility of PE closely follows the linear relationship, while iPP exhibits some divergence. All the semicrystalline samples were previously annealed at elevated temperatures for long periods (a few μs), and a strong effect of annealing was observed on the structure and the solubility of iPP. A well-developed iPP lamellar structure emerged at longer annealing times, while PE develops that structure already in the early crystallization stages. The solubility of semicrystalline iPP samples with lamellar morphology exhibited better agreement with extrapolated solubility values of the amorphous state─the extrapolation was made using a linear relationship connecting solubility in the purely amorphous phase and solubility in mixed phases (amorphous and crystalline). Results on the correlation of the solubility with the local structural ordering are also present.

Molecular Simulations of Hydrogen Sorption in Semicrystalline High-Density Polyethylene: The Impact of the Surface Fraction of Tie-Chains.

The modeling of the barrier properties of semicrystalline polymers has gained interest following the possible application of such materials as protective liners for the safe supply of pressurized hydrogen. The mass transport in such systems is intimately related to the complex intercalation between the crystal and amorphous phases, which was approached in this work through an all-atom representation of high-density polyethylene structures with a tailored fraction of amorphous-crystalline connections (tie-chains). Simulations of the polymer pressure-volume-temperature data and hydrogen sorption were performed by means of molecular dynamics and the Widom test particle insertion method. The discretization of the simulation domains of the semicrystalline structures allowed us to obtain profiles of density, degree of order, and gas solubility. The results indicated that the gas sorption in the crystalline regions is negligible and that the confinement of the amorphous phase between crystals induces a significant increase in density and a drop in the sorption capacity, even in the absence of tie-chains. Adding ties between the crystal and the amorphous phase results in further densification, an increase of the lamella tilt angle, and a decrease in the degree of crystallinity and hydrogen sorption coefficient, in agreement with several literature references.

Calculating Thermodynamic Factors for Diffusion Using the Continuous Fractional Component Monte Carlo Method.

Thermodynamic factors for diffusion connect the Fick and Maxwell-Stefan diffusion coefficients used to quantify mass transfer. Activity coefficient models or equations of state can be fitted to experimental or simulation data, from which thermodynamic factors can be obtained by differentiation. The accuracy of thermodynamic factors determined using indirect routes is dictated by the specific choice of an activity coefficient model or an equation of state. The Permuted Widom's Test Particle Insertion (PWTPI) method developed by Balaji et al. enables direct determination of thermodynamic factors in binary and multicomponent systems. For highly dense systems, for example, typical liquids, it is well known that molecular test insertion methods fail. In this article, we use the Continuous Fractional Component Monte Carlo (CFCMC) method to directly calculate thermodynamic factors by adopting the PWTPI method. The CFCMC method uses fractional molecules whose interactions with their surrounding molecules are modulated by a coupling parameter. Even in highly dense systems, the CFCMC method efficiently handles molecule insertions and removals, overcoming the limitations of the PWTPI method. We show excellent agreement between the results of the PWTPI and CFCMC methods for the calculation of thermodynamic factors in binary systems of Lennard-Jones molecules and ternary systems of Weeks-Chandler-Andersen molecules. The CFCMC method applied to calculate the thermodynamic factors of realistic molecular systems consisting of binary mixtures of carbon dioxide and hydrogen agrees well with the NIST REFPROP database. Our study highlights the effectiveness of the CFCMC method in determining thermodynamic factors for diffusion, even in densely packed systems, using relatively small numbers of molecules.

Rapid screening of gas solubility in ionic liquids using biased particle insertions with pre-sampled liquid trajectories

ABSTRACTWe present an efficient, general-purpose variant of the Widom test particle insertion method for computing chemical potentials of gaseous solutes in fluids or porous solids. The method is implemented in the Monte Carlo molecular simulation engine Cassandra, but receiving phase configurations are independent of this process and may be pre-sampled by other molecular simulation engines such as molecular dynamics codes. Efficiency enhancements present in this method include configurational biasing and accelerated atomic overlap detection. When applied to the estimation of Henry's law constants of atomistic difluoromethane and pentafluoroethane in ionic liquids, the accelerated overlap detection results in a speedup of more than an order of magnitude compared to conventional methods without sacrificing accuracy. We found good agreement between this method and Hamiltonian replica exchange (HREX) for Henry's law constant and absorption isotherm estimation. This embarrassingly parallel method is especially well suited for screening Henry's law constants of many small gases in the same solvents, since a liquid trajectory can be reused for as many solutes as desired.

Thermodynamics of cavity formation in different solvents: Enthalpy, entropy, and the solvophobic effects

For a long time thermodynamics of cavity formation in non-aqueous solvents remained virtually unexplored. We evaluated the entropic and enthalpic contributions into the Gibbs free energy of cavity formation in 22 structurally different organic solvents and water at 298 K. The calculations were carried out using 50–100 ns long molecular dynamics trajectories of solvents and Widom test particle insertion method with up to 1011 insertions per trajectory. It is shown that protic solvents with three-dimensional hydrogen bond networks, namely water, formamide, glycerol, ethylene glycol, propylene glycol, and diethylene glycol, behave differently from other considered solvents. They are characterized with lower values of the entropy of cavity formation than other considered solvents, which also leads to the higher Gibbs free energies. Water has a particularly low entropy as well as the lowest enthalpy of cavity formation at room temperature. We point out the similar tendencies observed for the thermodynamic functions of the cavities and solvation of real molecules in the considered solvents using entropy vs enthalpy and Gibbs free energy vs enthalpy plots. The cavity term governs the difference in solvation properties of protic hydrogen-bonded and aprotic solvents giving rise to the hydrophobic and more general solvophobic effects.

Molecular Simulations and Mechanistic Analysis of the Effect of CO2 Sorption on Thermodynamics, Structure, and Local Dynamics of Molten Atactic Polystyrene.

A simulation strategy encompassing different scales was applied to the systematic study of the effects of CO2 uptake on the properties of atactic polystyrene (aPS) melts. The analysis accounted for the influence of temperature between 450 and 550 K, polymer molecular weights (Mw) between 2100 and 31000 g/mol, and CO2 pressures up to 20 MPa on the volumetric, swelling, structural, and dynamic properties of the polymer as well as on the CO2 solubility and diffusivity by performing molecular dynamics (MD) simulations of the system in a fully atomistic representation. A hierarchical scheme was used for the generation of the higher Mw polymer systems, which consisted of equilibration at a coarse-grained level of representation through efficient connectivity-altering Monte Carlo simulations, and reverse-mapping back to the atomistic representation, obtaining the configurations used for subsequent MD simulations. Sorption isotherms and associated swelling effects were determined by using an iterative procedure that incorporated a series of MD simulations in the NPT ensemble and the Widom test particle insertion method, while CO2 diffusion coefficients were extracted from long MD runs in the NVE ensemble. Solubility and diffusivity compared favorably with experimental results and with predictions of the Sanchez–Lacombe equation of state, which was reparametrized to capture the Mw dependence of polymer properties with greater accuracy. Structural features of the polymer matrix were correctly reproduced by the simulations, and the effects of gas concentration and Mw on structure and local dynamics were thoroughly investigated. In the presence of CO2, a significant acceleration of the segmental dynamics of the polymer occurred, more pronouncedly at low Mw. The speed-up effect caused by the swelling agent was not limited to the chain ends but affected the whole chain in a similar fashion.

Sorption and permeation of gases in hyper-cross-linked hybrid poly(POSS-imide) networks: An in silico study

Networked hybrid materials containing various types of organic imides covalently bonded with polyhedral oligomeric silsesquioxanes (POSS) have recently been developed. One of their possible applications is for gas-separation under extreme conditions of pressure and/or temperature. Experimental results indicate that their permeation properties depend on the specific organic dianhydride precursor used in the synthesis. In this work, atomistic models of 6FDA-based and PMDA-based poly(POSS-imide) crosslinked networks have been used to investigate the sorption, desorption, relaxation and permeation properties of such materials at the molecular level and their dependence on the nature of the imide linker. Results have been compared to available experimental data.A highly-efficient parallelised version of the excluded volume map sampling test particle insertion (EVMS-TPI) method was used to obtain the infinite dilution solubilities of a number of small penetrants (N2, O2, CH4 and CO2) in these model networks. Analyses of the Boltzmann-weighted probability densities for the distribution of the insertion energies showed that very reliable estimates of solubilities could be obtained.Extensive molecular dynamics simulations were then performed to obtain the complete model sorption isotherms for CO2 and CH4 at 35 °C up to pressures of over 60 bar using an iterative technique which matches the chemical potential of the penetrant in the poly(POSS-imide) phase to that in the corresponding gas phase. The sorption-induced changes in volume, the significant space available to penetrants, the mean insertion energies and their distributions were characterized. Comparisons were not only made between both types of poly(POSS-imide) to characterize the effect of the linker, but also with (hypothetical) uncrosslinked mixtures of the POSS and dianhydride molecules in order to elucidate the effects of crosslinking. Given the known conditioning properties of CO2, desorption isotherms were also obtained for this penetrant. These showed significant hysteresis. Prolonged exposure to high concentrations of CO2 was also found to lead to slow relaxations that have significant effects on the solubility.Permeabilities were obtained for methane and carbon dioxide at 100 °C, 200 °C and 300 °C at a pressure of ~ 2 bar. Differences in permeabilities were largely diffusion dependent.

Chemical potentials of water, methanol, carbon dioxide and hydrogen sulphide at low temperatures using continuous fractional component Gibbs ensemble Monte Carlo

Chemical potentials of coexisting gas and liquid phases for water, methanol, hydrogen sulphide and carbon dioxide for the temperature range K to K are computed using two different methodologies: (1) Widom’s test particle insertion (WTPI) method in the conventional Gibbs Ensemble (GE), and (2) the Continuous Fractional Component Gibbs Ensemble Monte Carlo (CFCGE MC) method. It is shown that the WTPI method fails to accurately compute the chemical potentials of water and methanol in the liquid phase at low temperatures, while accurate chemical potentials in the liquid phase are computed using CFCGE MC method. For the CFCGE MC method, the statistical uncertainty for computed chemical potentials of water and methanol in the liquid phase are considerably smaller compared to the WTPI method. For the water models considered in this study (SPC, TIP3P-EW, TIP4P-EW, TIP5P-EW), computed excess chemical potentials based on three-site models are in better agreement with the chemical potentials computed from an empirical equation of state from the NIST database. For water, orientational biasing is applied during test particle insertion to check whether certain orientations of test particle are energetically unfavourable. A two-dimensional Overlapping Distribution Method (ODM) in the NVT ensemble is derived for this purpose. It is shown that failure of the WTPI method for systems with a strong hydrogen bonding network does not depend on orientation of the test molecule in that system. For all systems in this study, the WTPI method breaks down when the void fraction of the system drops below approximately 0.50.

Thermodynamics and simulation of hard-sphere fluid and solid: Kinetic Monte Carlo method versus standard Metropolis scheme.

The paper aims at a comparison of techniques based on the kinetic Monte Carlo (kMC) and the conventional Metropolis Monte Carlo (MC) methods as applied to the hard-sphere (HS) fluid and solid. In the case of the kMC, an alternative representation of the chemical potential is explored [E. A. Ustinov and D. D. Do, J. Colloid Interface Sci. 366, 216 (2012)], which does not require any external procedure like the Widom test particle insertion method. A direct evaluation of the chemical potential of the fluid and solid without thermodynamic integration is achieved by molecular simulation in an elongated box with an external potential imposed on the system in order to reduce the particle density in the vicinity of the box ends. The existence of rarefied zones allows one to determine the chemical potential of the crystalline phase and substantially increases its accuracy for the disordered dense phase in the central zone of the simulation box. This method is applicable to both the Metropolis MC and the kMC, but in the latter case, the chemical potential is determined with higher accuracy at the same conditions and the number of MC steps. Thermodynamic functions of the disordered fluid and crystalline face-centered cubic (FCC) phase for the hard-sphere system have been evaluated with the kinetic MC and the standard MC coupled with the Widom procedure over a wide range of density. The melting transition parameters have been determined by the point of intersection of the pressure-chemical potential curves for the disordered HS fluid and FCC crystal using the Gibbs-Duhem equation as a constraint. A detailed thermodynamic analysis of the hard-sphere fluid has provided a rigorous verification of the approach, which can be extended to more complex systems.

Liquid-crystal phase equilibria of Lennard-Jones chains

ABSTRACTLiquid-crystal phase equilibria of Lennard-Jones chain fluids and the solubility of a Lennard-Jones gas in the coexisting phases are calculated from Monte Carlo simulations. Direct phase equilibria calculations are performed using an expanded formulation of the Gibbs ensemble. Monomer densities, order parameters, and equilibrium pressures are reported for the coexisting isotropic and nematic phases of: (1) linear Lennard-Jones chains, (2) a partially-flexible Lennard-Jones chain, and (3) a binary mixture of linear Lennard-Jones chains. The effect of chain length is determined by calculating the isotropic-nematic coexistence of linear Lennard-Jones chain fluids made of 8, 10, and 12 segments (8-, 10-, 12-mer). The effect of molecular flexibility on the isotropic-nematic equilibrium is studied for a Lennard-Jones 10-mer chain fluid with one freely-jointed segment at the end of the chain. An isotropic-nematic phase split and fractionation are reported for a binary mixture of linear 7-mer and 12-mer chains. Simulation results are compared with theoretical results as obtained from a recently developed analytical equation of state based on perturbation theory. Excellent agreement between theory and simulations is observed. The solubility of a monomer Lennard-Jones gas in the coexisting isotropic and nematic phases is estimated using the Widom test-particle insertion method. A linear relationship between solubility difference and density difference at isotropic-nematic coexistence is observed. It is shown that gas solubility is independent of the nematic ordering of the fluid, at constant temperature and density conditions.

Molecular Simulation Studies of the Diffusion of Methane, Ethane, Propane, and Propylene in ZIF-8

ZIF-8 is a strong candidate for propane/propylene separation, which is regarded as one of the most industrially demanding. Molecular simulation of this separation must account for the flexibility of the structure, which enables the adsorption and diffusion of molecules with kinetic diameter larger than the apertures of the pores. Moreover, this simulation requires modeling subtle changes since the strong sieving effect upon the mixture depends on the very small differences between propane and propylene molecular sizes (∼0.2 A). In this work, a new force-field for the ZIF-8 structure has been developed from DFT calculations in simplified structures. The new parameter set reproduces structural properties in very good agreement with the experimental measurements reported in literature. Molecular dynamics simulations and the Widom test particle insertion method were then employed for the calculation of diffusivities, activation energies and adsorption properties of propane and propylene. The results are in ag...

Order and correlation contributions to the entropy of hydrophobic solvation.

The entropy of hydrophobic solvation has been explained as the result of ordered solvation structures, of hydrogen bonds, of the small size of the water molecule, of dispersion forces, and of solvent density fluctuations. We report a new approach to the calculation of the entropy of hydrophobic solvation, along with tests of and comparisons to several other methods. The methods are assessed in the light of the available thermodynamic and spectroscopic information on the effects of temperature on hydrophobic solvation. Five model hydrophobes in SPC/E water give benchmark solvation entropies via Widom's test-particle insertion method, and other methods and models are tested against these particle-insertion results. Entropies associated with distributions of tetrahedral order, of electric field, and of solvent dipole orientations are examined. We find these contributions are small compared to the benchmark particle-insertion entropy. Competitive with or better than other theories in accuracy, but with no free parameters, is the new estimate of the entropy contributed by correlations between dipole moments. Dipole correlations account for most of the hydrophobic solvation entropy for all models studied and capture the distinctive temperature dependence seen in thermodynamic and spectroscopic experiments. Entropies based on pair and many-body correlations in number density approach the correct magnitudes but fail to describe temperature and size dependences, respectively. Hydrogen-bond definitions and free energies that best reproduce entropies from simulations are reported, but it is difficult to choose one hydrogen bond model that fits a variety of experiments. The use of information theory, scaled-particle theory, and related methods is discussed briefly. Our results provide a test of the Frank-Evans hypothesis that the negative solvation entropy is due to structured water near the solute, complement the spectroscopic detection of that solvation structure by identifying the structural feature responsible for the entropy change, and point to a possible explanation for the observed dependence on length scale. Our key results are that the hydrophobic effect, i.e. the signature, temperature-dependent, solvation entropy of nonpolar molecules in water, is largely due to a dispersion force arising from correlations between rotating permanent dipole moments, that the strength of this force depends on the Kirkwood g-factor, and that the strength of this force may be obtained exactly without simulation.

Isotropic-nematic phase equilibria of hard-sphere chain fluids-Pure components and binary mixtures.

The isotropic-nematic phase equilibria of linear hard-sphere chains and binary mixtures of them are obtained from Monte Carlo simulations. In addition, the infinite dilution solubility of hard spheres in the coexisting isotropic and nematic phases is determined. Phase equilibria calculations are performed in an expanded formulation of the Gibbs ensemble. This method allows us to carry out an extensive simulation study on the phase equilibria of pure linear chains with a length of 7 to 20 beads (7-mer to 20-mer), and binary mixtures of an 8-mer with a 14-, a 16-, and a 19-mer. The effect of molecular flexibility on the isotropic-nematic phase equilibria is assessed on the 8-mer+19-mer mixture by allowing one and two fully flexible beads at the end of the longest molecule. Results for binary mixtures are compared with the theoretical predictions of van Westen et al. [J. Chem. Phys. 140, 034504 (2014)]. Excellent agreement between theory and simulations is observed. The infinite dilution solubility of hard spheres in the hard-sphere fluids is obtained by the Widom test-particle insertion method. As in our previous work, on pure linear hard-sphere chains [B. Oyarzún, T. van Westen, and T. J. H. Vlugt, J. Chem. Phys. 138, 204905 (2013)], a linear relationship between relative infinite dilution solubility (relative to that of hard spheres in a hard-sphere fluid) and packing fraction is found. It is observed that binary mixtures greatly increase the solubility difference between coexisting isotropic and nematic phases compared to pure components.

Temperature dependence of local solubility of hydrophobic molecules in the liquid-vapor interface of water.

One important aspect of the hydrophobic effect is that solubility of small, nonpolar molecules in liquid water decreases with increasing temperature. We investigate here how the characteristic temperature dependence in liquid water persists or changes in the vicinity of the liquid-vapor interface. From the molecular dynamics simulation and the test-particle insertion method, the local solubility Σ of methane in the liquid-vapor interface of water as well as Σ of nonpolar solutes in the interface of simple liquids are calculated as a function of the distance z from the interface. We then examine the temperature dependence of Σ under two conditions: variation of Σ at fixed position z and that at fixed local solvent density around the solute molecule. It is found that the temperature dependence of Σ at fixed z depends on the position z and the system, whereas Σ at fixed local density decreases with increasing temperature for all the model solutions at any fixed density between vapor and liquid phases. The monotonic decrease of Σ under the fixed-density condition in the liquid-vapor interface is in accord with what we know for the solubility of nonpolar molecules in bulk liquid water under the fixed-volume condition but it is much robust since the solvent density to be fixed can be anything between the coexisting vapor and liquid phases. A unique feature found in the water interface is that there is a minimum in the local solubility profile Σ(z) on the liquid side of the interface. We find that with decreasing temperature the minimum of Σ grows and at the same time the first peak in the oscillatory density profile of water develops. It is likely that the minimum of Σ is due to the layering structure of the free interface of water.

Contribution of the acetate anion to CO2 solubility in ionic liquids: theoretical method development and experimental study.

A new theoretical method was developed to compute the Henry's law constant for gas absorption in a solvent through strong nonphysical interactions. The new method was created by expanding the test particle insertion method typically applied to physisorbing systems to account for the strong intermolecular interactions present in chemisorbing systems. By using an ab initio (AI)-based Boltzmann-averaged potential to model the interaction between CO2 and the tetra-n-butylphosphonium acetate ([P4444][CH3COO]) ionic liquid, the total Henrys's law constant at 298 K was computed to be 0.011 to 0.039 bar, reasonably comparable to the experimental value of 0.18 bar measured in this work. Three different AI potentials were used to verify the applicability of this approach. In contrast, when a classical force field (FF) was used to describe the interaction between CO2 and [P4444][CH3COO], the Henry's law constant was computed to be 27 bar, significantly larger than the experimental value. The classical FF underestimates the CO2-[P4444][CH3COO] interaction compared with the AI calculations, which in turn leads to the smaller simulated CO2 solubility. Simulations further indicate that the CO2 interaction with the [CH3COO](-) anion is much stronger than with the [P4444](+) cation. This result strongly suggests that the large CO2 solubility in [P4444][CH3COO] is due to the strong CO2-[CH3COO](-) interaction.

The phase behavior of linear and partially flexible hard-sphere chain fluids and the solubility of hard spheres in hard-sphere chain fluids

The liquid crystal phase behavior of linear and partially flexible hard-sphere chain fluids and the solubility of hard spheres in hard-sphere chain fluids are studied by constant pressure Monte Carlo simulations. An extensive study on the phase behavior of linear fluids with a length of 7, 8, 9, 10, 11, 12, 13, 14, 15, and 20 beads is carried out. The phase behavior of partially flexible fluids with a total length of 8, 10, 14, and 15 beads and with different lengths for the linear part is also determined. A precise description of the reduced pressure and of the packing fraction change at the isotropic-nematic coexistence was achieved by performing long simulation runs. For linear fluids, a maximum in the isotropic to nematic packing fraction change is observed for a chain length of 15 beads. The infinite dilution solubility of hard spheres in linear and partially flexible hard-sphere chain fluids is calculated by the Widom test-particle insertion method. To identify the effect of chain connectivity and molecular anisotropy on free volume, solubility is expressed relative to that of hard spheres in a hard sphere fluid at same packing fraction as relative Henry's law constants. A linear relationship between relative Henry's law constants and packing fraction is observed for all linear fluids. Furthermore, this linearity is independent of liquid crystal ordering and seems to be independent of chain length for linear chains of 10 beads and longer. The same linear relationship was observed for the solubility of hard spheres in nematic forming partially flexible fluids for packing fractions up to a value slightly higher than the nematic packing fraction at the isotropic-nematic coexistence. At higher packing fractions, the small flexibility of these fluids seems to improve solubility in comparison with the linear fluids.

Calculating thermodynamic factors of ternary and multicomponent mixtures using the Permuted Widom test particle insertion method

We present an approach to calculate thermodynamic factors for multicomponent liquid mixtures from molecular simulations by the Permuted Widom test particle insertion method. The thermodynamic factors relate experimentally obtained Fick diffusion coefficients and the Maxwell–Stefan diffusion coefficients computed from molecular simulations/theory. In our earlier work, we presented the Permuted Widom test particle insertion method (Balaji et al. in Mol Phys 2012) for calculating thermodynamic factors for binary systems from a single simulation. In this paper, we have extended this approach to ternary systems consisting of Lennard–Jones particles. The thermodynamic factors computed from the simulations and from numerically differentiating the activity coefficients obtained from the conventional Widom single test particle insertion method were compared. An excellent agreement was observed.

A direct method for calculating thermodynamic factors for liquid mixtures using the Permuted Widom test particle insertion method

Understanding mass transport in liquids by mutual diffusion is an important topic for many applications in chemical engineering. The reason for this is that diffusion is often the rate limiting step in chemical reactors and separators. In multicomponent liquid mixtures, transport diffusion can be described by both generalized Fick's law and the Maxwell–Stefan theory. The Maxwell–Stefan and Fick approaches in an n-component system are related by the so-called thermodynamic factor [R. Taylor and H.A. Kooijman, Chem. Eng. Commun, 102, 87 (1991)]. As Fick diffusivities can be measured in experiments and Maxwell–Stefan diffusivities can be obtained from molecular simulations/theory, the thermodynamic factors bridge the gap between experiments and molecular simulations/theory. It is therefore desirable to be able to compute thermodynamic factors from molecular simulations. Unfortunately, presently used simulation techniques for computing thermodynamic factors are inefficient and often require numerical differentiation of simulation results. In this work, we propose a modified version of the Widom test-particle method to compute thermodynamic factors from a single simulation. This method is found to be more efficient than the conventional Widom test particle insertion method combined with numerical differentiation of simulation results. The approach is tested for binary systems consisting of Lennard-Jones particles. The thermodynamic factors computed from the simulation and from numerically differentiating the activity coefficients obtained from the conventional Widom test particle insertion method are in excellent agreement.

A Molecular Dynamics Technique to Extract Forces in Soft Matter Systems Under Compression With Constant Solvent Chemical Potential.

Molecular dynamics simulations of opposing polymer brushes at varying surface separation distances were performed to develop a method for conducting a static compression of soft matter. As all separation distances were represented by independent simulations, the proper solvent density for every level of compression needed to be determined to acquire realistic data. This was accomplished by maintaining a constant solvent chemical potential for each separation distance. In doing so, each independent simulation is equilibrated with all others, reproducing conditions encountered experimentally in force spectroscopy measurements. Chemical potential was determined using the Widom test particle insertion method. Force information was extracted from pressure profiles, such that unphysical forces occurring within the surface layers were not accounted for in the calculation. Each individual simulation was a canonical ensemble molecular dynamics simulation, but taken together they approximate a grand canonical ensemble for the solvent particles by holding their chemical potential constant.

Predicting Water Sorption and Volume Swelling in Dense Polymer Systems via Computer Simulation

Atomistic model structures of amorphous polyamide 6 (PA-6) and of an adhesive system consisting of the diglycidyl ether of bisphenol A (DGEBA) as epoxy resin and isophorone diamine (IPD) as a curing agent are generated. For the adhesive, we use a new approach for the generation of the cross-linked polymer networks. It takes into account the chemical reaction kinetics of the curing reaction and, therefore, results in more realistic network structures. On the basis of the corresponding model structures, the equilibrium water content and the swelling ratio of amorphous PA-6 and of the DGEBA+IPD networks are calculated via computer simulation for different thermodynamic conditions. A hybrid method is used combining the molecular dynamics technique with an accelerated test particle insertion method. The results are in reasonable agreement with experiments and, in the case of the PA-6 system, with results obtained via other computer simulation methods.