Monte carlo simulation of self-assembled ordered hybrid materials

Lattice Monte Carlo simulations were performed to study the structure of hybrid organic-inorganic materials. Several cases were modeled where the composition corresponds to high surfactant concentration phases similar to that obtained from the synthesis of hybrid materials resulting from a phase separation. When using hybrid inorganic precursors, comparable to organosilicas, we observe that the organic segment is well mixed with the inorganic precursor and surfactant heads and no preferential location of the organic groups is observed. We show that the behavior of surfactant/hybrid precursor systems is analogous to those where co-surfactants or co-solvents are used, and that the lack of ordering in some cases can be explained by the change in solvent quality when using hybrid precursors. A comparison of structural characterization of the different phases using several tools, such as aggregate size distribution, density profiles, and pair radial distribution function is presented.     

Monte carlo simulations of polydisperse polymers grafted on spherical surfaces

This article presents effects of polydispersity in polymers grafted on spherical surfaces on grafted polymer chain conformations, grafted layer thickness, and free‐end monomer distribution within the grafted layer. At brush‐like grafting densities, as polydispersity index (PDI) increases, the scaling exponent of radius of gyration of grafted chains approaches that of a single chain grafted on the same nanoparticle, because polydispersity alleviates monomer crowding within the brush. At high PDI, the chains shorter than the number average chain length, N_n, have more compressed conformations, and the chains longer than N_n overall stretch less than in the monodisperse case. As seen in polydisperse flat brushes at high grafting densities, the grafted layer thickness on spherical nanoparticle increases with PDI. Polydispersity eliminates the region near the surface devoid of free‐end monomers seen in monodisperse cases, and it reduces the width of free‐end monomer distribution and shifts the free‐end monomer distribution close to the surface. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Monte carlo simulations of gamma images: application to plutonium rich mixed carbide fuel

A novel approach for quantitave characterisation of nuclear fuel for homogeneity based on comparison of autoradiographs with simulated standards is proposed. Autoradiographs which are self generated radiation images are ususally only qualitatively evaluated due to the absence of calibration standards for reference which are not practically possible to synthesize in all possible variations owing to the difficulty in handling plutonium. Monte Carlo simulation for radiation transport has been attempted for the first time to estimate the minimum detectable limit (MDL) of pure PuC agglomerate and generate virtual autoradiographic images of plutonium rich mixed carbide fuel for visual comparison reference purpose.

     

Monte carlo simulations of micellization in model ionic surfactants: application to sodium dodecyl sulfate

A lattice model for ionic surfactants with explicit counterions is proposed for which the micellization behavior can be accurately determined from grand canonical Monte Carlo simulations. The model is characterized by a few parameters that can be adjusted to represent various linear surfactants with ionic headgroups. The model parameters have a clear physical interpretation and can be obtained from experimental data unrelated to micellization, namely, geometric information and solubilities of tail segments. As a specific example, parameter values for sodium dodecyl sulfate were obtained by optimizing for the solubility of hydrocarbons in water and the structural properties of dodecane. The critical micelle concentration (cmc), average aggregation number, degree of counterion binding, and their dependence on temperature were determined from histogram reweighting grand canonical Monte Carlo simulations and were compared to experimental results. The model gives the correct trend and order of magnitude for all quantities but underpredicts the cmc and aggregation number. We suggest ways to modify the model that may improve agreement with experimental values.     

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.     

Monte carlo simulation of three-body ion-ion recombination

The extent to which the use of models in discussing spectroscopic measurements of molecular reorientation in isotropic fluids is necessary, and possibly misleading, is discussed by considering the basic conditions imposed by molecular symmetry on reorientation.     

Monte carlo studies of two-dimensional commensurate-incommensurate transitions

A Monte Carlo study has been made of commensurate-incommensurate behaviour for Lennard-Jones particles on a square commensurate lattice. For a fixed barrier height a pronounced minimum is found in plots of R, the mean distance of the particles from the nearest commensurate site, against ε^*, where ε^* is the depth of the potential well. The initial fall in R with increase in ε^* is shown to be associated with the repulsive part of the inter-particle interaction.     

Monte carlo simulation of many-arm star polymers

The statistics of many-arm star polymers in dilute solution is studied by means of the enrichment algorithm [K. Ohno and K. Binder, J. Stat. Phys. 64 , 781 (1991)] using a Monte Carlo simulation technique. In this study, typically 50,000 samples of star polymers with 12–32 arms each of which has 50 bonds are generated. The results for the configuration-number exponent γ(f) are γ(12) = -3.4 ± 0.3, γ(18) = -8.9 ± 0.2 and γ(32) = -29 ± 2.     

Atomistic monte carlo simulations of polymethylene fluids

Atomistic Monte Carlo simulations of polymethylene fluids in various environments have been of great importance for improving our understanding of the molecular arrangements and conformations in systems of chain molecules. This paper describes models, calculations and principal results of a number of simulations performed in the last two decades for liquid alkanes in bulk and in the proximity of solid surfaces, and for systems of long chain amphiphilic molecules in monolayers, double layers, micelles and droplets.     

Phase diagram of self-assembled rigid rods on two-dimensional lattices: theory and Monte Carlo simulations

Monte Carlo simulations and finite-size scaling analysis have been carried out to study the critical behavior in a two-dimensional system of particles with two bonding sites that, by decreasing temperature or increasing density, polymerize reversibly into chains with discrete orientational degrees of freedom and, at the same time, undergo a continuous isotropic-nematic (IN) transition. A complete phase diagram was obtained as a function of temperature and density. The numerical results were compared with mean field (MF) and real space renormalization group (RSRG) analytical predictions about the IN transformation. While the RSRG approach supports the continuous nature of the transition, the MF solution predicts a first-order transition line and a tricritical point, at variance with the simulation results.     

Monte carlo data on ion-hydration shell characteristics

Ufa Aviation Institute. Translated from Zhurnal Strukturnoi Khimii, Vol. 29, No. 5, pp. 160–162, September–October, 1988.     

Mass spectrometry of surfactant aggregates

In contrast with the enormous amount of literature produced during many decades in the field of surfactant aggregation in liquid, liquid crystalline and solid phases, only a few investigations concerning surfactant self- assembling in the gas phase as charged aggregates have been carried out until now. This lack of interest is disappointing in view of the remarkable theoretical and practical importance of the inherent knowledge. The absence of surfactant-solvent interactions makes it easier to study the role of surfactant-surfactant forces in determining their peculiar self-assembling features as well as the ability of these assemblies to incorporate selected solubilizate molecules. Thus, the study of gas-phase surfactant and surfactant-solubilizate aggregates is a research subject which has exciting potential, including mass and energy transport in the atmosphere, origin of life and simulation of supramolecular aggregation in interstellar space. On the other hand, the structural and dynamic properties of surfactant aggregates in the gas phase could be exploited in a number of interesting applications such as atmospheric cleaning agents, transport and protection of pulmonary drugs or biomolecules and as nanoreactors for specialized chemical reactions in confined space. Spectrometric techniques, together with molecular dynamics simulations, have been the principal investigative tools in this field and appearto be particularly suited to gaining fundamental information on the structure and stability of surfactant-based supramolecular aggregates, charge state effects, entrapment of solubilizate molecules, preferential solubilization sites and chemical reactions localized in a single organized aggregate. The main aim of this review is to present the actual state of the art in this novel and exciting research field underlining the knowledge acquired up to now as well as the aspects needing a more deep understanding. Moreover, intriguing departures of the behavior of surfactant solutions under electrospray ionization conditions from that of ionic, polar and apolar analytes will be discussed.     

Monte Carlo simulations of Lennard-Jones nonionic surfactant adsorption at the liquid/vapor interface

We present Monte Carlo simulations of nonionic surfactant adsorption at the liquid/vapor interface of a monatomic solvent. All molecules in the system, solvent and surfactant, are characterized by the Lennard-Jones (LJ) potential using differing interaction parameters. Surfactant molecules consist of an amphiphilic chain with a solvophilic head and a solvophobic tail. Adjacent atoms along the surfactant chain are connected by finitely extensible harmonic springs. Solvent molecules move via the Metropolis random-walk algorithm, whereas surfactant molecules move according to the continuum configurational bias Monte Carlo (CBMC) method. We generate quantitative thermodynamic adsorption and surface tension isotherms in addition to surfactant radius of gyration, tilt angles, and potentials of mean force. Surface tension simulations compared to those calculated from the simulated adsorbed amounts and the Gibbs adsorption isotherm agree confirming equilibrium in our simulations. We find that the classical Langmuir isotherm is obeyed for our LJ surfactants over the range of head and tail lengths studied. Although simulated surfactant chains in the bulk solution exhibit random orientations, surfactant chains at the interface orient roughly perpendicular and the tails elongate compared to bulk chains even in the submonolayer adsorption regime. At a critical surfactant concentration, designated as the critical aggregation concentration (CAC), we find aggregates in the solution away from the interface. At higher concentrations, simulated surface tensions remain practically constant. Using the simulated potential of mean force in the submonolayer regime and an estimate of the surfactant footprint at the CAC, we predict a priori the Langmuir adsorption constant, KL, and the maximum monolayer adsorption, Gammam. Adsorption is driven not by proclivity of the surfactant for the interface, but by the dislike of the surfactant tails for the solvent, that is by a "solvophobic" effect. Accordingly, we establish that a coarse-grained LJ surfactant system mimics well the expected equilibrium behavior of aqueous nonionic surfactants adsorbing at the air/water interface.     

水-甲醇体系的Monte Carlo分子膜拟

本文应用Monte Carlo分子膜拟方法对水、甲醇水-甲醇的1:1混合物、甲醇无限稀释时的水溶液和水无限稀释时的甲醇溶液等五个体系进行了研究。采用TIP分子位能函数, 得到了上述体系的热力学性质、原子径向分布函数、分子氢键配位数分布。并以Monte Carlo分子模拟获得的结构函数与X射线衍射实验结果进行了比较。     

Monte carlo simulation of carboxylic acid phase equilibria

Configurational-bias Monte Carlo simulations were carried out in the Gibbs ensemble to generate phase equilibrium data for several carboxylic acids. Pure component coexistence densities and saturated vapor pressures were determined for acetic acid, propanoic acid, 2-methylpropanoic acid, and pentanoic acid, and binary vapor-liquid equilibrium (VLE) data for the propanoic acid + pentanoic acid and 2-methylpropanoic acid + pentanoic acid systems. The TraPPE-UA force field was used, as it has recently been extended to include parameters for carboxylic acids. To simulate the branched compound 2-methylpropanoic acid, certain minor assumptions were necessary regarding angle and torsion terms involving the -CH- pseudo-atom, since parameters for these terms do not exist in the TraPPE-UA force field. The pure component data showed good agreement with the available experimental data, particularly with regard to the saturated liquid densities (mean absolute errors were less than 1.1%). On average, the predicted critical temperature and density were within 1% of the experimental values. All of the binary simulations showed good agreement with the experimental x-y data. However, the TraPPE-UA force field predicts saturated vapor pressures of pure components that are larger than the experimental values, and consequently the P-x-y and T-x-y data of the binary systems also deviate from the measured data.     

Monte carlo simulation of macromolecular ionization by nanoelectrospray

Electrospray ionization (ESI) is commonly used in macromolecular mass spectrometry, yet the dynamics of macromolecules in ESI droplets are not well understood. In this study, a Monte Carlo based model was developed, which can predict the efficiency of electrospray ionization for macromolecules, i.e., the number of macromolecular ions produced per macromolecules electrosprayed. The model takes into account ESI droplet evaporation, macromolecular diffusion within the droplet, droplet fissions, and the statistical nature of the ESI process. Two idealized representations of macromolecular analytes were developed, describing cluster prone, droplet surface inactive macromolecules and droplet surface active macromolecules, respectively. It was found that surface active macromolecules are preferentially ionized over surface inactive cluster prone macromolecules when the initial droplet size is large and the analyte concentration in solution is high. Simulations showed that ESI efficiency decreases with increasing initial droplet size and analyte molecular weight, and is influenced by analyte surface activity, the properties of the solvent, and the variance of the droplet size distribution. Model predictions are qualitatively supported by experimental measurements of macromolecular electrospray ionization made previously. Overall, this study demonstrates the potential capabilities of Monte Carlo based ESI models. Future developments in such models will allow for more accurate predictions of macromolecular ESI intensity.     

Programmed construction of discrete self-assembled cyclic aggregates

Examples for programmed self-assembly of alpha,beta-unsaturated ketoximes carrying a terminal pyridine subunit are described. The solid-state structures of a series of pyridinyl alpha,beta-unsaturated ketoximes 1 were investigated by X-ray structural analysis; this demonstrated that molecules 3-5 predictably form discrete cyclic aggregates stabilized by a network of hydrogen bonds and aromatic interactions.