Influence of monovalent ion size on colloidal forces probed by Monte Carlo simulations

The present work studies the role of ionic size in the interactions between the electrical double layers of colloids immersed into electrolyte solutions of monovalent ions. Such interactions are studied by means of Monte Carlo (MC) simulations and the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. Despite the omission of the steric effects and some other features of real electrolyte solutions, DLVO theory is known to work qualitatively well for 1:1 electrolyte solutions. However, this affirmation is based on previous tests where an ionic diameter around 0.4 nm was taken for all ionic species. In contrast, some experimental studies suggest that larger hydrated ions should be considered and even specified for each type of ion. In this work, the importance of ionic size is analyzed by applying the primitive model of electrolyte to the intermediate region between a pair of equally charged infinite planar surfaces. The double layer interactions were calculated from the ionic densities at the distance of closest approach to the charged surfaces, this method constitutes an alternative to the traditional calculations at the midplane. Our MC simulations predict the existence of negative net pressures for monovalent electrolytes in the case of zero charge density. In addition, MC simulations reveal some disagreements with theoretical predictions for ionic diameters larger than 0.4 nm. These discrepancies can become significant if surface charge density is large enough due to the restructuration of the double layer. The physical mechanisms for these deviations are also discussed.     

Extensional flow of bulk polymers studied by dynamic Monte Carlo simulations

Dynamic Monte Carlo simulations of bulk lattice polymers driven through planar geometries with sequentially converging,parallel and diverging spaces between two neutrally repulsive solid plates are reported.The spatial profiles of polymer velocity and deformation along the course of such a laminar extensional flow have been carefully analyzed.The results appear consistent with experimental observations in literature.In the entrance and exit regions,a linear dependence of chain extension upon the excess velocity has been observed.Moreover,an annexed shear flow and a molecular-dispersion effect are found.The results demonstrate a useful strategy of this approach to study polymer flows and bring new insights into the non-Newtonian-fluid behaviors of bulk polymers in capillary rheometers and micro-fluidic devices.     

Melting behaviors of icosahedral metal clusters studied by Monte Carlo simulations

We present an atom‐resolved analysis method that traces physical quantities such as the root‐mean‐square bond length fluctuation and coordination number for individual atoms as functions of temperature or time. This method is applied to explain the temperature‐dependent behaviors of three types of Ni_N (N=12,13,14) clusters. The detailed studies for the three types of clusters reveal characteristics as follows: (a) as the temperature increases, all three types of clusters undergo two‐stage melting, irrespective of the existence of vacancy or adatom on the icosahedral surfaces, (b) the melting of icosahedral clusters with vacancy starts with vacancy hopping, which has not been observed for any type of small clusters (N<34), (c) the melting of the icosahedral clusters with adatom (N=14) is initiated by adatom hopping, followed by the site exchange between the adatom and surface atoms. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 380–387, 2000     

Confined crystallization of cylindrical diblock copolymers studied by dynamic Monte Carlo simulations

We report dynamic Monte Carlo simulations of polymer crystallization confined in the cylindrical microdomains of diblock copolymers. The microdomains were prepared via spontaneous microphase separation from homogeneous melt, and the major component was then frozen in a vitreous amorphous state to make a hard confinement to the crystallization of the minor component. We found that during the isothermal crystallization at high temperatures, crystal orientations are dominantly perpendicular to the cylinder axis at the early stage of crystal nucleation and remain to the final state; while if the block junctions are broken before crystallization, crystal orientations are dominantly parallel at the early stage of crystal nucleation, and eventually other orientations take the place of parallel preferences. Analysis of bond orientations in the heterogeneous melts demonstrates the microscopic origin of oriented crystal nucleation.     

Dynamics of polyethylene melts studied by Monte Carlo simulations on a high coordination lattice

A detailed comparison is made between the experiment, prior simulations by other groups, and our simulation based on a newly designed dynamic Monte Carlo algorithm, on the dynamics of polyethylene (PE) melts. The new algorithm, namely, noncross random two-bead move has been developed on a high coordination lattice (the 2nnd lattice) for studying the dynamics of realistic polymers. The chain length (molecular weight) in our simulation ranges from C40 (562 Da) to C324 (4538 Da). The effects of finite chain length have been confirmed and significant non-Gaussian statistics evidently results in nonstandard static and dynamic properties of short PE chains. The diffusion coefficients scale with molecular weight (M) to the −1.7 power for short chains and −2.2 for longer chains, which coincides very well with experimental results. No pure Rouse scaling in diffusion has been observed. The transitional molecular weight to the entanglement regime is around 1500 Da. The detailed mean square displacements of middle bead (g1) are presented for several chain lengths. The reptation-like slowdown can be clearly observed only above M ∼ 2400 Da. The slope 0.25 predicted by the theory for the intermediate regime is missing; instead a slope close to 0.4 appears, indicating that additional relaxation mechanism exists in this transitional region. The relaxation times extracted by fitting the autocorrelation function of end-to-end vectors with reptation model scale with M to 2.5 for long chains, which seemingly conflicts with the scaling of diffusion. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2556–2571, 2006     

Macroion solutions in the cell model studied by field theory and Monte Carlo simulations

Aqueous solutions of charged spherical macroions with variable dielectric permittivity and their associated counterions are examined within the cell model using a field theory and Monte Carlo simulations. The field theory is based on separation of fields into short- and long-wavelength terms, which are subjected to different statistical-mechanical treatments. The simulations were performed by using a new, accurate, and fast algorithm for numerical evaluation of the electrostatic polarization interaction. The field theory provides counterion distributions outside a macroion in good agreement with the simulation results over the full range from weak to strong electrostatic coupling. A low-dielectric macroion leads to a displacement of the counterions away from the macroion.     

Adsorption of binary mixtures on two-dimensional surfaces: theory and Monte Carlo simulations

The adsorption of binary mixtures on square lattices is studied by combining theoretical modeling and Monte Carlo (MC) simulations in grand canonical ensemble. The adsorption thermodynamics is analyzed through the total and partial isotherms. Two theoretical models have been used in the present study: (i) the first, which we called cluster approximation (CA), is based on exact calculations of configurations on finite cells. An efficient algorithm allows us to calculate the detailed structure of the configuration space for \(m = l \times l \) cells; and (ii) the second is a generalization of the classical quasi-chemical approximation (QCA) in which the adsorbate is a binary mixture of species \(a\) and \(b\) . Adsorbate–adsorbate lateral interactions are incorporated in the context of the two mentioned approximations. Results from CA and QCA are compared with MC simulations. Close agreement between simulated and theoretical data supports the validity of the theoretical models to describe the adsorption of mixed gases on two-dimensional surfaces.     

Divalent and mixed divalent/monovalent conduction in β″-alumina: A Monte Carlo study

A Monte Carlo method is developed for simulation of mixed ionic conductivity in β″-alumina-type materials. The conduction plane of these materials is represented by a lattice gas model in which monovalent and divalent cation carriers diffuse via a vacancy mechanism and interact through a nearest-neighbor coulombic repulsion. By comparing experimental data for pure Na⁺ and pure Ba²⁺ β″-aluminas with simulation results, it is possible to estimate the near-neighbor interaction energies ε_i and jump barriers U_i for both kinds of ions. On the basis of these estimations the total ionic conductivity of Na⁺Ba²⁺ β″-alumina is calculated as a function of temperature and concentration of carriers. As Ba²⁺ replaces Na⁺, the conductivity initially increases as more vacancies become available. For very high temperatures, this increase continues until exchange is complete; but at lower temperatures, the conductivity reaches a peak for some optimal ${}^{\mathrm{<mtext>Ba2+</mtext>}}{}^{\mathrm{<mtext>Na+</mtext>}}$ composition, and then drops off as the number of Ba²⁺, and hence the strength of ionic correlation, goes up. The presence of ordering in the fully exchanged (all Ba²⁺) case manifests itself in substantial curvature of the Arrhenius plots for conductivity. The activation energy for conductivity as a function of Ba²⁺ mole fraction (X_Ba²⁺) shows a pronounced rise near an X value of $\text{2}\text{3}$, in agreement with recent experimental observations.     

Adsorption of branched-linear polyethyleneimine-ethylene oxide conjugate on hydrophilic silica investigated by ellipsometry and Monte Carlo simulations

The adsorption of and conformation adopted by a branched-linear polymer conjugate to the hydrophilic silica-aqueous solution interface have been studied by in situ null ellipsometry and Monte Carlo simulations. The conjugate is a highly branched polyethyleneimine structure with ethyleneoxide chains grafted to its primary and secondary amino groups. In situ null ellipsometry demonstrated that the polymer conjugate adsorbs to the silica surface from water and aqueous solution of 1 mM asymmetric divalent salt (calcium and magnesium chloride to emulate hard water) over a large pH range. The adsorbed amount is hardly affected by pH and large charge reversal on the negatively charged silica surface occurred at pH = 4.0, due to the adsorption of the cationic polyelectrolyte. The Monte Carlo simulations using an appropriate coarse-grained model of the polymer in solution predicted a core-shell structure with no sharp boundary between the ethyleneimine and ethyleneoxide moieties. The structure at the interface is similar to that in solution when the polymer degree of protonation is low or moderate while at high degree of protonation the strong electrostatic attraction between the ethyleneimine core and oppositely charged silica surface distorts the ethyleneoxide shell so that an "anemone"-like configuration is adopted. The adsorption of alkyl benzene sulfonic acid (LAS) to a preadsorbed polymer layer was also investigated by null ellipsometry. The adsorption data brought additional support for the existence of a strong polymer adsorption and showed the presence of a binding which was further enhanced by the decreased solvency of the surfactant in the salt solution and confirmed the surface charge reversal by the polymer adsorption at pH = 4.0.     

Morphology and mechanical properties of surfactant aggregates at water-silica interfaces: molecular dynamics simulations

Dilute and concentrated surfactant systems at the solid-liquid interface are examined using classical molecular dynamics simulations. Particular emphasis is placed on understanding how surfactants aggregate and form the micellar structure, how micelles change shape at high concentrations in aqueous media and in the presence of hydrophilic surfaces, and at what force this micellar structure breaks apart during indentation of micelle-covered surfaces with a proximal probe microscope tip. The specific system of interest is C12TAB (n-dodecyltrimethylammonium bromide) surfactant in an aqueous medium that is modeled with empirical potentials. The simulations predict that the micelle structure in water is compact and either spherical or elliptical in shape. In the presence of a hydrophilic surface of silica, the structure evolves into a flat elliptical shape, in agreement with experimental findings. The simulated indentation of the micelle/silica system causes the micelle to break apart at an indentation force of about 1 nN and form a surfactant monolayer. The predicted force curve is in excellent agreement with experimental measurements.     

Adsorption of octyl cyanide at the free water surface as studied by Monte Carlo simulation

Monte Carlo simulations of the adsorption layer of octyl cyanide have been performed on the canonical (N, V, T) ensemble at 300 K. The systems simulated cover the range of octyl cyanide surface densities from 0.27 to 7.83 mumol/m2. The surface density value at which the saturation of the adsorption layer occurs is estimated to be 1.7 mumol/m2. At low surface densities, the main driving force of the adsorption is found to be the formation of hydrogen bonds between the water and octyl cyanide molecules, whereas at higher surface concentrations, the dipole-dipole attraction between the neighboring adsorbed octyl cyanide molecules becomes more important. At low surface concentrations, the water-octyl cyanide hydrogen bonds prefer tilted alignments relative to the interface; however, in the case of the saturated adsorption layer, the number of such hydrogen bonds is maximized, leading to the preference of these bonds for the orientation perpendicular to the interface. Contrary to nonionic surfactants of multiple hydrogen bonding abilities (e.g., 1-octanol, C8E3), the increasing surface concentration of octyl cyanide was not found to lead to considerable competition of the molecules for positions of optimal arrangement. As a consequence, the energy and geometry of the water-octyl cyanide hydrogen bonds are found to be insensitive to the octyl cyanide surface concentration.     

Adsorption of 1-octanol at the free water surface as studied by Monte Carlo simulation

The adsorption of 1-octanol at the free water surface has been investigated by Monte Carlo computer simulation. Six different systems, built up by an aqueous and a vapor phase, the latter also containing various number of octanol molecules, have been simulated. The number of the octanol molecules has been chosen in such a way that the octanol surface density varies in a broad range, between 0.27 and 7.83 micromol/m(2) in the six systems simulated. For reference, the interfacial system containing bulk liquid octanol in the apolar phase has also been simulated. The results have shown that the formation of hydrogen bonds between the interfacial water and adsorbed octanol molecules is of key importance in determining the properties of the adsorbed layer. At low octanol surface concentration values all the octanol molecules are strongly (i.e., by hydrogen bonds) bound to the aqueous phase, whereas their interaction with each other is negligibly small. Hence, they are preferentially oriented in such a way that their own binding energy (and thus their own free energy) is minimized. In this preferred orientation the O-H bond of the octanol molecule points flatly toward the aqueous phase, declining by about 30 degrees from the interfacial plane, irrespectively from whether the octanol molecule is the H-donor or the H-acceptor partner in the hydrogen bond. Hence, in its preferred orientation the octanol molecule can form at least two low energy hydrogen bonds with water: one as a H-donor and another one as a H-acceptor. Moreover, the preferred orientation of the hydrogen bonded water partners is close to one of the two preferred interfacial water alignments, in which the plane of the water molecule is parallel with the interface. When increasing the octanol surface density, the water surface gets saturated with hydrogen bonded octanols, and hence any further octanol molecule can just simply condense to the layer of the adsorbed octanols. The surface density value at which this saturation occurs is estimated to be about 1.7 micromol/m(2). Above this surface density value the hydrogen bonded octanols and their water partners are oriented in such a way that the number of the water-octanol hydrogen bonds is maximized. Hence, the preferred alignment of the O...O axes of these hydrogen bonds is perpendicular to the interface. This orientation is far from the optimal alignment of the individual octanol molecules, which is also reflected in the observed fact that, unlike in the case of many other adsorbents, the average molecular binding energy of the adsorbed octanol molecules increases (i.e., becomes less negative) with increasing octanol surface density.     

Dynamics of polypropylene chains in their binary blends of different stereochemical sequences studied by Monte Carlo simulations

The conformational and dynamic properties of polypropylene （PP） for both pure melts and blends with different chain tacticity were investigated by Monte Carlo simulation of isotactic （iPP）, atactic （aPP） and syndiotactic （sPP） polypropylenes. The simulation of coarse-grained PP models was performed on a high coordination lattice incorporating short- and long-range intramolecular interactions from the rotational isomeric state （RIS） model and Lennard-Jones （LJ） potential function of propane pairs, respectively. The dynamics of chains in binary PP/PP mixture were investigated with the composition of C150H302 with different chain taciticity. The diffusion rates of PP with different stereochemistry are generally in the order as： iPP 〉 aPP 〉〉 sPP. For PP/PP blends with 50：50 wt% binary mixtures, immiscibility was observed when sPP was introduced into the mixtures. The diffusion rate of iPP and aPP became slower after mixing, while sPP diffuses significantly faster in the binary mixtures. The mobility of PP chains depends on both intramolecular （molecular size and chain stiffness） and intermolecular （chain packing） interactions. The effect of intramolecular contribution is greater than that of intermolecular contribution for iPP and aPP chains in binary mixtures. For sPP chain, intermolecular interaction has greater influence on the dynamics than intramolecular contribution.     

Surface tension of associating fluids by Monte Carlo simulations

Canonical Monte Carlo (NVT-MC) simulations were performed to obtain surface tension and coexistence densities at the liquid-vapor interface of one-site associating Lennard-Jones and hard-core Yukawa fluids, as functions of association strength and temperature. The method to obtain the components of the pressure tensor from NVT-MC simulations was validated by comparing the equation of state of the associative hard sphere system with that coming from isothermal-isobaric Monte Carlo simulations. Surface tension of the associative Lennard-Jones fluid determined from NVT-MC is compared with previously reported results obtained by molecular dynamics simulations of a pseudomixture model of monomers and dimers. A good agreement was found between both methods. Values of surface tension of associative hard-core Yukawa fluids are presented here for the first time.     

Multilayer polyelectrolyte films as nanofiltration membranes for separating monovalent and divalent cations

Alternating adsorption of polyanions and polycations on porous supports provides a convenient way to prepare ion-selective nanofiltration membranes. This work examines optimization of ultrathin, multilayer polyelectrolyte films for monovalent/divalent cation separations relevant to water softening. Membranes composed of five bilayers of poly(styrene sulfonate)/poly(allylamine hydrochloride) (PSS/PAH) on porous alumina supports allow a solution flux of 0.85 m³/(m² day) at 4.8 bar, and exhibit 95% rejection of MgCl₂ along with a Na⁺/Mg²⁺ selectivity of 22. Similar results were obtained in Na⁺/Ca²⁺ separations. PSS/poly(diallyl-dimethylammonium chloride) (PDADMAC) films permit higher fluxes than PSS/PAH systems due to the higher swelling of films containing PDADMAC, but the Mg²⁺ rejection by PSS/PDADMAC membranes is less than 45%. However, capping PSS/PDADMAC films with a bilayer of PSS/PAH yields Mg²⁺ rejections and Na⁺/Mg²⁺ selectivities that are typical of pure PSS/PAH membranes. Separation performance can be optimized through control over deposition conditions (pH and supporting electrolyte concentration) and the charge of the outer layer since Donnan exclusion is a major factor in monovalent/divalent cation selectivity. Streaming potential measurements demonstrate that the magnitude of positive surface charge increases with increasing concentrations of Mg²⁺ in solution or when the outer polycation layer is deposited from a solution of high ionic strength.     

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.     

Adsorption of adenine derivatives on a gold electrode studied by specular reflectivity measurement

The adsorption of adenine, deoxyadenosine, deoxyadenosine-5′-monophosphate,-diphosphate and-triphosphate on a gold electrode has been studied by specular reflectivity measurement in 0.1 M NaClO₄ solution. In the presence of these compounds, a marked decrease in reflectivity was found on reflectivity-potential curves in the potential region more positive than ?0.8 V vs. Ag/AgCl, the decrease being ascribed to the adsorption of them. The magnitude of change in reflectivity was dependent on both the concentration and the electrode potential. The reflectivity change observed in the negative potential region was analyzed quantitatively according to the procedure previously described. The results were elucidated on the basis of the same isotherm as used by Green and Dahms in their adsorption study of aromatic hydrocarbons, and the number of solvent molecules being replaced through the adsorption of one organic molecule and the free energy change of adsorption were obtained. The former is suggestive of a flat orientation of the adsorbed molecule in contact with its adenine moiety on the electrode surface. It is also suggested from the latter that the presence of phosphate groups leads to a decrease in ΔG_ad⁰ resulting from their hydrophilic properties and a repulsive interaction between these groups and the negative charges on the surface.     

Effect of monovalent and divalent cations on the photoinactivation of bacteria with meso-substituted cationic porphyrins

It is well established that for successful photoinactivation (PI) of gram-negative bacteria a cationic photosensitizer is required. This requirement suggests a charge-dependent interaction between the photosensitizer and the gram-negative bacterium, which may be influenced by the presence of ions in the suspending medium. The aim of the present study was to investigate the effect of cations Na+ and Ca2+ on the efficacy of the PI of the gram-negative Pseudomonas aeruginosa and the gram-positive Staphylococcus aureus. The bacteria were suspended in buffer containing either meso-tetra(N-methyl-4-pyridyl)-porphyrin or meso-mono-phenyl-tri(N-methyl-4-pyridyl)-porphyrin as photosensitizer and various concentrations of Na+ or Ca2+. The cell suspensions were exposed to a broadband light dose of 9 J/cm2. In buffer without added cations, P. aeruginosa and S. aureus were equally sensitive to PI. Addition of cations strongly decreased the sensitivity of both bacteria to PI, with the PI of P. aeruginosa being much more decreased than that of S. aureus, and Ca2+ being more effective than Na+. The decreased sensitivity was accompanied by a reduced binding of the photosensitizers to the bacteria.     

Monte Carlo simulations of low background detectors

An implementation of the Electron Gamma Shower 4 code (EGS4) has been developed to allow convenient simulation of typical gamma ray measurement systems. Coincidence gamma rays, beta spectra, and angular correlations have been added to adequately simulate a complete nuclear decay and provide corrections to experimentally determined detector efficiencies. This code has been used to strip certain low-background spectra for the purpose of extremely low-level assay. Monte Carlo calculations of this sort can be extremely successful since low background detectors are usually free of significant contributions from poorly localized radiation sources, such as cosmic muons, secondary cosmic neutrons, and radioactive construction or shielding materials. Previously, validation of this code has been obtained from a series of comparisons between measurements and blind calculations. An example of the application of this code to an exceedingly low background spectrum stripping will be presented.Pacific Northwest Laboratory is operated by Battelle Memorial Institute for the US Department of Energy under Contract DE-AC06-76RLO 1830.