Vapour-Liquid Phase Equilibria of n-alkanes by Direct Monte Carlo Simulations

Abstract

We report results of direct Monte Carlo simulations of n-pentane and n-decane at the liquidvapour interface for a number of temperatures. The intermolecular interactions are modeled using the last version of the anisotropic united atom model (AUA4). We have used the local long range correction energy and an algorithm allowing to select randomly with equal probability two different displacements. The liquid and vapour densities are in excellent agreement with experimental data and with those previously calculated using the GEMC method.     

Basic statistics and variational concepts behind the reverse Monte Carlo technique

We revise the statistical foundations of the reverse Monte Carlo (RMC) technique by constructing the associated functional of a variational principle which incorporates, without any ad hoc assumptions, the inherent errors accompanying the simulation and the experimental data. We propose a Bayesian criteria for acceptance/rejection of configurations, in terms of the variations of the functional. The loss function and variational functional minimization approaches are compared.     

Spherical Boundary Conditions: A Finite and System Size Independent Geometry for Simulations of Electrolytic Liquids

The system size dependence of the thermodynamic properties of electrolytic systems in "Spherical Boundary Conditions" (SBC) have been examined in this work. Coulombic systems were simulated with different system sizes (N ) for a wide range of concentrations, different ionic charges and solvent permittivities. The effects of system size upon the mean internal energy ?U? and mean ionic activity coefficient ?γ _±? values were determined. Our results indicated that there was no dependence of the thermodynamic properties upon system size in the non-Euclidean geometry. Different methods of extrapolating the thermodynamic properties to infinite numbers were studied in detail. In contrast to the Euclidean geometry, the method of extrapolating ?U? or ?γ _±? versus N ^-1, N ^-2/3 and N ^-1/3 was not statistically justifiable nor advantageous. Therefore, SBC is well suited to simulating systems involving long-range electrostatic interactions because the results do not dependent upon system size.     

Properties of Confined Square-well Fluids using the Gibbs Ensemble Simulation Technique

In this work we have used the extension of the Gibbs ensemble simulation technique to inhomogeneous fluids [Panagiotopoulos, A.Z. (1987) "Adsorption and capillary condensation of fluid in cylindrical pores by Monte Carlo simulation in the Gibbs ensemble", Mol. Phys. , 62 (3), 701-719], which has been applied to adsorption phenomena of confined fluids. Fluid molecules are described by spherical particles interacting via a square-well potential. The fluid is confined in two types of walls: symmetrical (two hard walls) and non-symmetrical (one square-well wall and one hard wall). In order to analyze the behavior of the confined fluid by varying the potential parameters, we evaluated the bulk and confined densities, the internal energies and the density profiles for different supercritical temperatures. A variety of adsorption profiles can be obtained by using this model. The simulation data reported here complements the available simulation data for this system and can be useful in the development of inhomogeneous fluid theories. Since the square-well parameters can be related to real molecules this system can also be used to understand real adsorption systems.     

A Modified Valley Restrained Monte Carlo Method to Efficiently Search the Low Energy Structures of Peptides

Abstract

A new Monte Carlo sampling scheme, namely the Modified Valley Restrained Monte Carlo procedure, is used to obtain the global energy minimum conformations for polypeptides, such as Met-enkephalin and Melittin. For each peptide, we found close agreement with previous results from both theoretical and experimental studies. The simple idea for controlling the step size according to the Valley Function, provides useful suggestions in searching the global energy minimum structures, and furthermore helps solve the multiple minima problem.     

A simple model of the hydrophobic effect for molecular simulation of interfacial phenomena

A coarse-grained model for simulation of interfacial phenomena in aqueous systems has been developed. The model captures the hydrophobic effect by only considering the structure and cohesiveness of water. Monte Carlo (MC) simulations of water-oil mixtures show that low concentrations of oil are solvated with little perturbation of the hydrogen bonding network structure of the water, while high concentrations of oil are excluded altogether. Analysis of the water structure in the simulations indicates that the water molecules maintain close to four coordination in the presence of solutes and the distribution of bond angles is not markedly affected by the presence of solutes. MC simulations of an alkane oligomer in water and a poly(ethylene oxide) (PEO) oligomer in water indicate that the chains are quite flexible and also do not perturb the network structure of the water phase.     

“Nanostructures in Thin-film Epitaxy: Exploring and Exploiting Substrate-mediated Interactions”

We review recent advances made in understanding the ramifications of substrate-mediated interactions for thin-film growth. Experimental studies and first-principles calculations with density-functional theory (DFT) indicate that substrate-mediated interactions can significantly influence thin-film growth. We review the findings from our kinetic Monte Carlo simulations used to model the growth of thin films, both with and without substrate-mediated interactions. For Ag heteroepitaxy on Pt(1?1?1), the pair interaction energies and adsorbate diffusion barriers were obtained from DFT calculations. Island densities for this system show significant deviations from what is predicted by classical nucleation theory. The electronic interactions created by the adsorbed atoms lead to the formation of repulsive barriers surrounding small islands and, as a result, sharp island-size distributions are produced. The island-size distributions can be manipulated by changing the growth conditions to yield desirable island sizes and shapes.     

Mathematical modelling and computer simulation of Brownian motion and hybridisation of nanoparticle–bioprobe–polymer complexes in the low concentration limit

A wide variety of nano-biotechnological applications are being developed for nanoparticles based on in vitro diagnostic and imaging systems. Some of these systems allow highly sensitive detection of molecular biomarkers. Frequently, the very low concentration of the biomarkers makes very difficult the mathematical simulation of the motion of nanoparticles based on classical, partial differential equations. We address the issue of incubation times for low concentration systems using Monte Carlo simulations. We describe a mathematical model and computer simulation of Brownian motion of nanoparticle–bioprobe–polymer contrast agent complexes and their hybridisation to immobilised targets. We present results for the dependence of incubation times on the number of particles available for detection, and the geometric layout of the DNA-detection assay on the chip.     

Monte Carlo study of sI hydrogen hydrates

In this study, we perform grand canonical Monte Carlo simulations to evaluate the hydrogen storage capacity of structure I (sI) hydrogen hydrates at pressures up to 500 MPa. Initially, we calculate the upper limit of H₂ content of sI hydrates by studying the hypothetical sI hydrate, where H₂ is the single guest component. It is found that the storage capacity of the hypothetical pure H₂ sI hydrate could reach 3.5 wt% at 500 MPa and 274 K. Depending on pressure, the large cavities of the pure H₂ hydrate can accommodate up to three H₂ molecules while the small ones are singly occupied at most, even at pressures as high as 500 MPa, without any double occupancy being observed. Subsequently, the binary H₂–ethylene oxide (EO) hydrate is examined. In this case, the large cavities are occupied by a single EO molecule while the small cavities can accommodate at most a single H₂ molecule. Such configuration results in a maximum H₂ content of only 0.37 wt%. The hydrogen storage capacity does not improve significantly even in case when EO is replaced by a component with smaller molecular weight.