RESUMO
Extensive molecular dynamics simulations in the equilibrium isobaric-isothermal (NPT) ensemble were developed to determine the coexistence temperatures of the water hydrogen mixture using the direct coexistence method. The water molecules were modeled using the four-site TIP4P/Ice analytical potential, and the hydrogen molecules were described using a three-site potential. The simulations were performed at different pressures (p = 900, 1500, 3000, and 4000 bars). At each pressure, a series of simulations were developed at different temperatures (from 230 to 270 K). Our results followed a line parallel to the experimental coexistence temperatures and underestimated these temperatures by approximately 25 K in the investigated range. The final configurations could or could not contain a fluid phase depending on the pressure, in accordance with the phase diagram. In addition, we explored the dynamics of the H2 molecules through clathrate hydrate cages and observed different behaviors of the H2 molecules in the small cages and the large cages of the sII structure.
RESUMO
Molecular simulations in the canonical and isothermal-isobaric ensembles were performed to study the effect of varying the shape of the intermolecular potential on the phase diagram, critical, and interfacial properties of model fluids. The molecular interactions were modeled by the Approximate Non-Conformal (ANC) theory potentials. Unlike the Lennard-Jones or Morse potentials, the ANC interactions incorporate parameters (called softnesses) that modulate the steepness of the potential in their repulsive and attractive parts independently. This feature allowed us to separate unambiguously the role of each region of the potential on setting the thermophysical properties. In particular, we found positive linear correlation between all critical coordinates and the attractive and repulsive softness, except for the critical density and the attractive softness which are negatively correlated. Moreover, we found that the physical properties related to phase coexistence (such as span of the liquid phase between the critical and triple points, variations in the P-T vaporization curve, interface width, and surface tension) are more sensitive to changes in the attractive softness than to the repulsive one. Understanding the different roles of attractive and repulsive forces on phase coexistence may contribute to developing more accurate models of liquids and their mixtures.
RESUMO
Na-montmorillonite hydrates in presence of ethane molecules are studied by means of hybrid Monte Carlo simulations in the NP(zz)T and muP(zz)T ensembles. The NP(zz)T ensemble allows us to study the interlaminar distance as a function of water and ethane content. These data show clear plateaus for lower ethane contents and mainly for water contents consistent with the formation of a single water layer. In addition, from this ensemble the structure for some of these interlaminar compositions were analyzed. For systems containing few ethane molecules and water enough to complete a single layer, it was observed that ethane mainly situates close to the interlayer midplane and adopts a nearly parallel arrangement to the clay surface. On the other hand, the muP(zz)T ensemble allows us to determine the interlaminar distance and water-ethane content for any specific reservoir. Here, some important findings are the following: the partial exchange of water by ethane molecules that enhances for decreasing the water vapor pressure; the obtention of a practically constant interlaminar space distance as a function of the water vapor pressure; the conservation of ion solvation shells; the enhancement of the water-ethane exchange for burial conditions; and finally, the incapability for a dehydrated clay mineral to swell in a dry and rich ethane atmosphere.
