Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFree energy determination of phase coexistence in model C60: A comprehensive Monte Carlo study
137
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The free energy of the solid and fluid phases of the Girifalco C60 model are determined through extensive Monte Carlo simulations. In this model the molecules interact through a spherical pair potential, characterized by a narrow and attractive well, adjacent to a harshly repulsive core. We have used the Widom test particle method and a mapping from an Einstein crystal, in order to estimate the absolute free energy in the fluid and solid phases, respectively; we have then determined the free energy along several isotherms, and the whole phase diagram, by means of standard thermodynamic integrations. We highlight how the interplay between the liquid-vapor and the liquid-solid coexistence conditions determines the existence of a narrow liquid pocket in the phase diagram, whose stability is assessed and confirmed in agreement with previous studies. In particular, the critical temperature follows closely an extended corresponding-states rule recently outlined by Noro and Frenkel [J. Chem. Phys. 113:2941 (2000)]. We discuss the emerging energetic properties of the system, which drive the phase behavior in systems interacting through short-range forces [A. A. Louis, Phil. Trans. R. Soc. A 359:939 (2001)], in order to explain the discrepancy between the predictions of several structural indicators and the results of full free energy calculations, to locate the fluid phase boundaries. More generally, we aim to provide extended reference data for calculations of the free energy of the C60 fullerite in the low temperature regime, as for the determination of the phase diagram of higher order fullerenes and other fullerene-related materials, whose description is based on the same model adopted in this work.
rate research
Read More
We have investigated the phase diagram of the Girifalco model of C60 fullerene in the framework provided by the MHNC and the SCOZA liquid state theories, and by a Perturbation Theory (PT), for the free energy of the solid phase. We present an extended assessment of such theories as set against a recent Monte Carlo study of the same model [D. Costa et al, J. Chem. Phys. 118:304 (2003)]. We have compared the theoretical predictions with the corresponding simulation results for several thermodynamic properties. Then we have determined the phase diagram of the model, by using either the SCOZA, or the MHNC, or the PT predictions for one of the coexisting phases, and the simulation data for the other phase, in order to separately ascertain the accuracy of each theory. It turns out that the overall appearance of the phase portrait is reproduced fairly well by all theories, with remarkable accuracy as for the melting line and the solid-vapor equilibrium. The MHNC and SCOZA results for the liquid-vapor coexistence, as well as for the corresponding critical points, are quite accurate. All results are discussed in terms of the basic assumptions underlying each theory. We have selected the MHNC for the fluid and the first-order PT for the solid phase, as the most accurate tools to investigate the phase behavior of the model in terms of purely theoretical approaches. The overall results appear as a robust benchmark for further theoretical investigations on higher order C(n>60) fullerenes, as well as on other fullerene-related materials, whose description can be based on a modelization similar to that adopted in this work.
The non-equilibrium dynamics of condensation phenomena in nano-pores is studied via Monte Carlo simulation of a lattice gas model. Hysteretic behavior of the particle density as a function of the density of a reservoir is obtained for various pore geometries in two and three dimensions. The shape of the hysteresis loops depend on the characteristics of the pore geometry. The evaporation of particles from a pore can be fitted to a stretched exponential decay of the particle density $rho_f(t) sim exp [ -(t/tau)^beta]$. Phase separation dynamics inside the pore is effectively described by a random walk of the non-wetting phases. Domain evolution is significantly slowed down in presence of random wall-particle potential and gives rise to a temperature dependent growth exponent. On the other hand roughness of the pore wall only delays the onset of a pure domain growth.
In this paper we present a classical Monte Carlo simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles and periodic boundary conditions. We used a recently developed algorithm which apart from standard Metropolis local moves employs also global moves consisting of displacements of the center of mass of the whole chains in all three spatial directions as well as rotations of the chains around an axis parallel to the crystallographic c-direction. Our simulations are performed in the NpT ensemble, at zero pressure, and extend over the whole range of temperatures in which the orthorhombic phase is experimentally known to be stable (10 - 450 K). In order to investigate the finite-size effects in this extremely anisotropic crystal, we used different system sizes and different chain lengths, ranging from C_12 to C_96 chains, the total number of atoms in the super-cell being between 432 and 3456. We show here the results for structural parameters, such as the orthorhombic cell parameters a,b,c, and the setting angle of the chains, as well as internal parameters of the chains, such as the bond lengths and angles. Among thermodynamic quantities, we present results for thermal expansion coefficients, elastic constants and specific heat. We discuss the temperature dependence of the measured quantities as well as the related finite-size effects. In case of lattice parameters and thermal expansion coefficients, we compare our results to those obtained from other theoretical approaches as well as to some available experimental data. We also suggest some possible ways of extending this study.
A new Monte Carlo approach is proposed to investigate the fluid-solid phase transition of the polydisperse system. By using the extended ensemble, a reversible path was constructed to link the monodisperse and corresponding polydisperse system. Once the fluid-solid coexistence point of the monodisperse system is known, the fluid-solid coexistence point of the polydisperse system can be obtained from the simulation. The validity of the method is checked by the simulation of the fluid-solid phase transition of a size-polydisperse hard sphere colloid. The results are in agreement with the previous studies.
The `dynamic Hubbard Hamiltonian describes interacting fermions on a lattice whose on-site repulsion is modulated by a coupling to a fluctuating bosonic field. We investigate one such model, introduced by Hirsch, using the determinant Quantum Monte Carlo method. Our key result is that the extended s-wave pairing vertex, repulsive in the usual static Hubbard model, becomes attractive as the coupling to the fluctuating Bose field increases. The sign problem prevents us from exploring a low enough temperature to see if a superconducting transition occurs. We also observe a stabilization of antiferromagnetic correlations and the Mott gap near half-filling, and a near linear behavior of the energy as a function of particle density which indicates a tendency toward phase separation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
