Phase separation in binary mixtures containing polymers: A quantitative comparison of single-chain-in-mean-field simulations and computer simulations of the corresponding multichain systems

We discuss a phenomenological, coarse-grained simulation scheme, single-chain-in-mean-field (SCMF) simulation, for investigating the kinetics of phase separation in dense polymer blends and mixtures of polymers and solvents. In the spirit of self-consistent-field calculations, we approximate the interacting multichain problem by that of a single chain in an external field, which, in turn, depends on the local densities of the components. To study the time evolution of the mixture, we perform an explicit Monte Carlo (MC) simulation of an ensemble of independent chains in the external field and periodically calculate the average densities and update the external field. Unlike dynamic self-consistent-field theory, these SCMF simulations do not assume that the chain conformations relax much more quickly than the density and incorporate the single-chain dynamics explicitly rather than via an Onsager coefficient. This allows us to study systems with large spatial inhomogeneities and dynamic asymmetries. To assess the accuracy and limitations of the simulation scheme, we compare the results of SCMF simulations using a discretized Edwards Hamiltonian with computer simulations of the corresponding multichain system for (1) the early stages of spinodal decomposition of a symmetric binary polymer blend in response to a quench from χN = 0.314 to χN = 5 (where χ is the Flory–Huggins parameter and N is the number of segments), for which the growth rate of composition fluctuations is compared with MC simulations of the bond fluctuation model and alternative dynamic self-consistent-field calculations, and (2) the evaporation of a solvent from a low-molecular-weight thin polymer film, for which a comparison is made with molecular dynamics (MD) simulations of a bead-necklace model with a monomeric solvent. In the latter case, the polymer conformations are extracted from MD simulations and modeled in the SCMF simulations by a discretized Edwards Hamiltonian augmented by a chain-bending potential. From the MD simulations of thin polymer films in equilibrium with its vapor, phase coexistence has been determined, and the second- and third-order virial coefficients in the SCMF simulations have been adjusted accordingly. Finally, MD simulations of bulk solutions of a polymer and a solvent over a range of compositions, as well as the pure solvent at various densities, have been performed to determine self-diffusion coefficients that enter the SCMF simulations in the form of density-dependent segmental mobilities. A comparison of the polymer and solvent profiles in a thin film as a function of time and the fraction of the solvent evaporating from a solvent-swollen film, as obtained from MD simulations and parameterized SCMF simulations, shows satisfactory agreement for this simple mapping procedure. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 934–958, 2005     

Monte Carlo simulations of cylinder-forming ABC triblock terpolymer thin films

We systematically study the cylinder-forming ABC triblock terpolymer thin films using canonical ensemble Monte Carlo simulations. The simulated annealing procedure is applied to the self-assembling process. By judicious choice of the system dimensions, we elaborately investigate the effect of film thickness on the orientation of the cylinders. This confined triblock terpolymer system exhibits different phase behavior under the weak and strong surface fields. In addition, we also investigate the ensemble-averaged chain orientations and relative density profiles.     

Molecular simulations of confined liquids: an alternative to the grand canonical Monte Carlo simulations

Commonly, the confinement effects are studied from the grand canonical Monte Carlo (GCMC) simulations from the computation of the density of liquid in the confined phase. The GCMC modeling and chemical potential (μ) calculations are based on the insertion/deletion of the real and ghost particle, respectively. At high density, i.e., at high pressure or low temperature, the insertions fail from the Widom insertions while the performing methods as expanded method or perturbation approach are not efficient to treat the large and complex molecules. To overcome this problem we use a simple and efficient method to compute the liquid's density in the confined medium. This method does not require the precalculation of μ and is an alternative to the GCMC simulations. From the isothermal-isosurface-isobaric statistical ensemble we consider the explicit framework/liquid external interface to model an explicit liquid's reservoir. In this procedure only the liquid molecules undergo the volume changes while the volume of the framework is kept constant. Therefore, this method is described in the Np(n)AV(f)T statistical ensemble, where N is the number of particles, p(n) is the normal pressure, V(f) is the volume of framework, A is the surface of the solid/fluid interface, and T is the temperature. This approach is applied and validated from the computation of the density of the methanol and water confined in the mesoporous cylindrical silica nanopores and the MIL-53(Cr) metal organic framework type, respectively.     

Application of a two-length-scale field theory to the solvation of neutral and charged molecules

We develop a continuous self-consistent theory of solute-water interactions that allows determination of the hydrophobic layer around a solute molecule of any geometry, with an explicit account of solvent structure described by its correlation function. We compute the mean solvent density profile n(r) surrounding the solute molecule as well as its solvation free energy deltaG. We compare the two-length-scale field theory to the numerical data of Monte Carlo simulations found in the literature for spherical molecules and discuss the possibility of self-consistent adjustment of the free parameters of the theory. In the framework of this approach, we compute the solvation free energies of alkane molecules and the free energy of interaction of two spheres of radius R separated by the distance D. We describe the general setting of the self-consistent account of electrostatic interactions in the framework of our model where the water is considered not as a continuous medium but as a gas of dipoles. We analyze the limiting cases where the proposed theory coincides with the electrostatics of a continuous medium.     

Ordering of Diblock Copolymer Materials on Patterned Substrates: a Single Chain in Mean Field Simulation Study

Summary: The directed assembly of diblock copolymers on patterned substrates is a way to create nanoscopically structured materials. We study the structure and kinetics of diblock copolymers on patterned substrates by simulating a large ensemble of independent chains in an external field. This external field depends on the density created by the ensemble of molecules and it is frequently updated as to mimic the instantaneous interactions of a molecule with its neighbors. This approximate, particle-based field theoretical method allows (i) to incorporate arbitrary chain architecture (ii) to include fluctuations and (iii) the explicit propagation of the chain conformations in time permits us to study the kinetics of structure formation. The factors that control the accuracy of the method are quantitatively discussed and the reconstruction of the soft morphology at substrate patterns that deviate from the periodic morphology of the diblock in the bulk are illustrated.     

Mean field theory-based calculation of FLC polarization

Mean field theory and Monte Carlo sampling are applied to the calculation of the spontaneous polarization density of ferroelectric liquid crystals by the ensemble averaging of single molecules confined in mean field potentials reflecting the SmC environment. Molecules are modelled with atomistic detail, using intramolecular interaction potentials derived from ab initio quantum mechanical calculations. This technique is applied to thirteen members of a family of novel fluoro ether-based compounds. Comparison with experiment shows that the observed variation of polarization density with chemical structure is well reproduced in most cases, but that the observed temperature dependence of polarization density is not captured by our model. The features of molecular organization responsible for the discrepancies between theory and experiment are discussed.     

Mean field theory-based calculation of FLC polarization

Mean field theory and Monte Carlo sampling are applied to the calculation of the spontaneous polarization density of ferroelectric liquid crystals by the ensemble averaging of single molecules confined in mean field potentials reflecting the SmC environment. Molecules are modelled with atomistic detail, using intramolecular interaction potentials derived from ab initio quantum mechanical calculations. This technique is applied to thirteen members of a family of novel fluoro ether-based compounds. Comparison with experiment shows that the observed variation of polarization density with chemical structure is well reproduced in most cases, but that the observed temperature dependence of polarization density is not captured by our model. The features of molecular organization responsible for the discrepancies between theory and experiment are discussed.     

Prediction of adsorption of small molecules in porous materials based on ab initio force field method

Computational prediction of adsorption of small molecules in porous materials has great impact on the basic and applied research in chemical engineering and material sciences. In this work,we report an approach based on grand canonical ensemble Monte Carlo(GCMC) simulations and ab initio force fields. We calculated the adsorption curves of ammonia in ZSM-5 zeolite and hydrogen in MOF-5(a metal-organic-framework material). The predictions agree well with experimental data. Because the predictions are based on the first principle force fields,this approach can be used for the adsorption prediction of new molecules or materials without experimental data as guidance.     

基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

Computer simulations of ligand–protein binding with ensembles of protein conformations: A Monte Carlo study of HIV‐1 protease binding energy landscapes

We present the results of molecular docking simulations with HIV‐1 protease for the sb203386 and skf107457 inhibitors by Monte Carlo simulated annealing. A simplified piecewise linear energy function, the standard AMBER force field, and the AMBER force field with solvation and a soft‐core smoothing component are employed in simulations with a single‐protein conformation to determine the relationship between docking simulations with a simple energy function and more realistic force fields. The temperature‐dependent binding free energy profiles of the inhibitors interacting with a single protein conformation provide a detailed picture of relative thermodynamic stability and a distribution of ligand binding modes in agreement with experimental crystallographic data. Using the simplified piecewise linear energy function, we also performed Monte Carlo docking simulations with an ensemble of protein conformations employing preferential biased sampling of low‐energy protein conformations, and the results are analyzed in connection with the free energy profiles. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 73–84, 1999     

Reactive Monte Carlo and grand-canonical Monte Carlo simulations of the propene metathesis reaction system

The influence of silicalite-1 pores on the reaction equilibria and the selectivity of the propene metathesis reaction system in the temperature range between 300 and 600 K and the pressure range from 0.5 to 7 bars has been investigated with molecular simulations. The reactive Monte Carlo (RxMC) technique was applied for bulk-phase simulations in the isobaric-isothermal ensemble and for two phase systems in the Gibbs ensemble. Additionally, Monte Carlo simulations in the grand-canonical ensemble (GCMC) have been carried out with and without using the RxMC technique. The various simulation procedures were combined with the configurational-bias Monte Carlo approach. It was found that the GCMC simulations are superior to the Gibbs ensemble simulations for reactions where the bulk-phase equilibrium can be calculated in advance and does not have to be simulated simultaneously with the molecules inside the pore. The confined environment can increase the conversion significantly. A large change in selectivity between the bulk phase and the pore phase is observed. Pressure and temperature have strong influences on both conversion and selectivity. At low pressure and temperature both conversion and selectivity have the highest values. The effect of confinement decreases as the temperature increases.     

Isobaric-isothermal monte carlo simulations from first principles: application to liquid water at ambient conditions.

A series of first-principles Monte Carlo simulations in the isobaric-isothermal ensemble were carried out for liquid water at ambient conditions (T=298 K and p=1 atm). The Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and norm-conserving Goedecker-Teter-Hutter (GTH) pseudopotentials were employed with the CP2 K simulation package to examine systems consisting of 64 water molecules. The fluctuations in the system volume encountered in simulations in the isobaric-isothermal ensemble require a reconsideration of the suitability of the typical charge-density cutoff and the regular grid-generation method previously used for the computation of the electrostatic energy in first-principles simulations in the microcanonical or canonical ensembles. In particular, it is noted that a much higher cutoff is needed and that the most computationally efficient method of creating grids can result in poor simulations. Analysis of the simulation trajectories using a very large charge-density cutoff at 1200 Ry and four different grid-generation methods point to a significantly underestimated liquid density of about 0.8 g cm-3 resulting in a somewhat understructured liquid (with a value of about 2.7 for the height of the first peak in the oxygen-oxygen radial distribution function) for BLYP-GTH water at ambient conditions. In addition, a simulation using a charge-density cutoff at 280 Ry yields a higher density of 0.9 g cm-3, showing the sensitivity of the simulation outcome to this parameter.     

Computer simulations of two-dimensional trifunctional bimodal networks

Transformation of mean chain dimensions in a bimodal Gaussian network is analyzed by generating a two-dimensional trifunctional network using Monte Carlo technique and deforming it, subject to the balance of mean forces at each junction. Instantaneous fluctuations in chain dimensions are not considered in the simulations. The problem is, therefore, equivalent to the analysis of a two-dimensional deterministic system of two sets, short and long, of linear springs, joined at trifunctional junctions. The molecular deformation of network chains, the distribution of orientation of chain vectors and the reduced stresses are studied as a function of extension ratio and network topology.     

Isotropic-nematic phase transition in the Lebwohl-Lasher model from density of states simulations

Density of states Monte Carlo simulations have been performed to study the isotropic-nematic (IN) transition of the Lebwohl-Lasher model for liquid crystals. The IN transition temperature was calculated as a function of system size using expanded ensemble density of states simulations with histogram reweighting. The IN temperature for infinite system size was obtained by extrapolation of three independent measures. A subsequent analysis of the kinetics in the model showed that the transition occurs via spinodal decomposition through aggregation of clusters of liquid crystal molecules.     

Chain conformation and solvent partitioning in reversed-phase liquid chromatography: Monte Carlo simulations for various water/methanol concentrations

Many structural models for the stationary phase in reversed-phase liquid chromatography (RPLC) systems have been suggested from thermodynamic and spectroscopic measurements and theoretical considerations. To provide a molecular picture of chain conformation and solvent partitioning in a typical RPLC system, a particle-based Monte Carlo simulation study is undertaken for a dimethyl octadecyl (C(18)) bonded stationary phase on a model siliceous substrate in contact with mobile phases having different methanol/water concentrations. Following upon previous simulations for gas-liquid chromatography and liquid-liquid phase equilibria, the simulations are conducted using the configurational-bias Monte Carlo method in the Gibbs ensemble and the transferable potentials for phase equilibria force field. The simulations are performed for a chain surface density of 2.9 micromol/m(2), which is a typical bonded-phase coverage for mono-functional alkyl silanes. The solvent concentrations used here are pure water, approximately 33 and 67% mole fraction of methanol and pure methanol. The simulations show that the chain conformation depends only weakly on the solvent composition. Most chains are conformationally disordered and tilt away from the substrate normal. The interfacial width increases with increasing methanol content and, for mixtures, the solvent shows an enhancement of the methanol concentration in a 10 Angstrom region outside the Gibbs dividing surface. Residual surface silanol groups are found to provide hydrogen bonding sites that lead to the formation of substrate bound water and methanol clusters, including bridging clusters that penetrate from the solvent/chain interfacial region all the way to the silica surface.     

Density functional theory of adsorption in pillared slit-like pores

We propose a density functional theory to describe adsorption of Lennard-Jones fluid in pillared slit like pores. Specifically, the pillars are built of chains that are bonded by their ends to the opposite pore walls. The approach we propose combines theory of quenched-annealed systems and theory of nonuniform fluids involving chain molecules. We compare the results of theoretical predictions with grand canonical ensemble Monte Carlo simulations and compute theoretical capillary condensation phase diagrams for several model systems.     

A new method for generating volume changes in isobaric-isothermal Monte Carlo simulations of flexible molecules

A novel method has been developed to generate volume fluctuations as required in isobaric-isothermal Monte Carlo simulations of flexible (or semi-flexible) chain molecules. While conventional global volume moves entail symmetric changes of the simulation box size, in the new method proposed here, volume changes consist of one-dimensional fluctuations of one of two sections of the system. One of the two sections is a thin slab where all molecular rearrangements are restrained to occur. A deterministric method has been implemented to reposition the molecular segments contained within the slab after a volume change is proposed. The performance of the new method is examined by studying the evolution towards equilibrium of homopolymer hard-core systems upon application of a step increase in the pressure. The new method is shown to be an order of magnitude faster than conventional methods in attaining the equilibrium density of our test systems.     

Monte Carlo simulations of liquid crystals near rough walls

The effect of surface roughness on the structure of liquid crystalline fluids near solid substrates is studied by Monte Carlo simulations. The liquid crystal is modeled as a fluid of soft ellipsoidal molecules and the substrate is modeled as a hard wall that excludes the centers of mass of the fluid molecules. Surface roughness is introduced by embedding a number of molecules with random positions and orientations within the wall. It is found that the density and order near the wall are reduced as the wall becomes rougher, i.e., the number of embedded molecules is increased). Anchoring coefficients are determined from fluctuations in the reciprocal space order tensor. It is found that the anchoring strength decreases with increasing surface roughness.