919 research outputs found

    An improved Monte Carlo method for direct calculation of the density of states

    Full text link
    We present an efficient Monte Carlo algorithm for determining the density of states which is based on the statistics of transition probabilities between states. By measuring the infinite temperature transition probabilities--that is, the probabilities associated with move proposal only--we are able to extract excellent estimates of the density of states. When this estimator is used in conjunction with a Wang-Landau sampling scheme [F. Wang and D. P. Landau, Phys. Rev. Lett. 86, 2050 (2001)], we quickly achieve uniform sampling of macrostates (e.g., energies) and systematically refine the calculated density of states. This approach requires only potential energy evaluations, continues to improve the statistical quality of its results as the simulation time is extended, and is applicable to both lattice and continuum systems. We test the algorithm on the Lennard-Jones liquid and demonstrate good statistical convergence properties.Comment: 7 pages, 4 figures. to appear in Journal of Chemical Physic

    Optimized intermolecular potential for nitriles based on Anisotropic United Atoms model

    Get PDF
    An extension of the Anisotropic United Atoms intermolecular potential model is proposed for nitriles. The electrostatic part of the intermolecular potential is calculated using atomic charges obtained by a simple Mulliken population analysis. The repulsion-dispersion interaction parameters for methyl and methylene groups are taken from transferable AUA4 literature parameters [Ungerer et al., J. Chem. Phys., 2000, 112, 5499]. Non-bonding Lennard-Jones intermolecular potential parameters are regressed for the carbon and nitrogen atoms of the nitrile group (–C≡N) from experimental vapor-liquid equilibrium data of acetonitrile. Gibbs Ensemble Monte Carlo simulations and experimental data agreement is very good for acetonitrile, and better than previous molecular potential proposed by Hloucha et al. [J. Chem. Phys., 2000, 113, 5401]. The transferability of the resulting potential is then successfully tested, without any further readjustment, to predict vapor-liquid phase equilibrium of propionitrile and n-butyronitrile

    On the Wang-Landau Method for Off-Lattice Simulations in the "Uniform" Ensemble

    Full text link
    We present a rigorous derivation for off-lattice implementations of the so-called "random-walk" algorithm recently introduced by Wang and Landau [PRL 86, 2050 (2001)]. Originally developed for discrete systems, the algorithm samples configurations according to their inverse density of states using Monte-Carlo moves; the estimate for the density of states is refined at each simulation step and is ultimately used to calculate thermodynamic properties. We present an implementation for atomic systems based on a rigorous separation of kinetic and configurational contributions to the density of states. By constructing a "uniform" ensemble for configurational degrees of freedom--in which all potential energies, volumes, and numbers of particles are equally probable--we establish a framework for the correct implementation of simulation acceptance criteria and calculation of thermodynamic averages in the continuum case. To demonstrate the generality of our approach, we perform sample calculations for the Lennard-Jones fluid using two implementation variants and in both cases find good agreement with established literature values for the vapor-liquid coexistence locus.Comment: 21 pages, 4 figure

    Crowding of Polymer Coils and Demixing in Nanoparticle-Polymer Mixtures

    Full text link
    The Asakura-Oosawa-Vrij (AOV) model of colloid-polymer mixtures idealizes nonadsorbing polymers as effective spheres that are fixed in size and impenetrable to hard particles. Real polymer coils, however, are intrinsically polydisperse in size (radius of gyration) and may be penetrated by smaller particles. Crowding by nanoparticles can affect the size distribution of polymer coils, thereby modifying effective depletion interactions and thermodynamic stability. To analyse the influence of crowding on polymer conformations and demixing phase behaviour, we adapt the AOV model to mixtures of nanoparticles and ideal, penetrable polymer coils that can vary in size. We perform Gibbs ensemble Monte Carlo simulations, including trial nanoparticle-polymer overlaps and variations in radius of gyration. Results are compared with predictions of free-volume theory. Simulation and theory consistently predict that ideal polymers are compressed by nanoparticles and that compressibility and penetrability stabilise nanoparticle-polymer mixtures.Comment: 18 pages, 4 figure

    Thermodynamics of Electrolytes on Anisotropic Lattices

    Full text link
    The phase behavior of ionic fluids on simple cubic and tetragonal (anisotropic) lattices has been studied by grand canonical Monte Carlo simulations. Systems with both the true lattice Coulombic potential and continuous-space 1/r1/r electrostatic interactions have been investigated. At all degrees of anisotropy, only coexistence between a disordered low-density phase and an ordered high-density phase with the structure similar to ionic crystal was found, in contrast to recent theoretical predictions. Tricritical parameters were determined to be monotonously increasing functions of anisotropy parameters which is consistent with theoretical calculations based on the Debye-H\"uckel approach. At large anisotropies a two-dimensional-like behavior is observed, from which we estimated the dimensionless tricritical temperature and density for the two-dimensional square lattice electrolyte to be Ttri=0.14T^*_{tri}=0.14 and ρtri=0.70\rho^*_{tri} = 0.70.Comment: submitted to PR

    Universality class of criticality in the restricted primitive model electrolyte

    Full text link
    The 1:1 equisized hard-sphere electrolyte or restricted primitive model has been simulated via grand-canonical fine-discretization Monte Carlo. Newly devised unbiased finite-size extrapolation methods using temperature-density, (T, rho), loci of inflections, Q = ^2/ maxima, canonical and C_V criticality, yield estimates of (T_c, rho_c) to +- (0.04, 3)%. Extrapolated exponents and Q-ratio are (gamma, nu, Q_c) = [1.24(3), 0.63(3); 0.624(2)] which support Ising (n = 1) behavior with (1.23_9, 0.630_3; 0.623_6), but exclude classical, XY (n = 2), SAW (n = 0), and n = 1 criticality with potentials phi(r)>Phi/r^{4.9} when r \to \infty
    corecore