Statistical thermodynamics of adsorbates with nonsymmetrical lateral interactions

Water adsorbate phases on ZnO and impact of vapor pressure on the equilibrium shape of nanoparticles

ZnO nanoparticles are used as catalysts and have potential applications in gas-sensing and solar energy conversion. A fundamental understanding of the exposed crystal facets, their surface chemistry, and stability as a function of environmental conditions is essential for rational design and improvement of synthesis and properties. We study the stability of water adsorbate phases on the non-polar low-index (101¯0) and (112¯0) surfaces from low coverage to multilayers using ab initio thermodynamics. We show that phonon contributions and the entropies due to a 2D lattice gas at low coverage and multiple adsorbate configurations at higher coverage have an important impact on the stability range of water adsorbate phases in the (T,p) phase diagram. Based on this insight, we compute and analyze the possible growth mode of water films for pressures ranging from UHV via ambient conditions to high pressures and the impact of water adsorption on the equilibrium shape of nanoparticles in a humid environment. A 2D variant of the Wulff construction shows that the (101¯0) and (112¯0) surfaces coexist on 12-faceted prismatic ZnO nanoparticles in dry conditions, while in humid environment, the (101¯0) surface is selectively stabilized by water adsorption resulting in hexagonal prisms.

Lattice-gas model of nonadditive interacting particles on nanotube bundles

In the present paper, the adsorption thermodynamics of a lattice-gas model which mimics a nanoporous environment is studied by considering nonadditive interactions between the adsorbed particles. It is assumed that the energy linking a certain atom with any of its nearest neighbors strongly depends on the state of occupancy in the first coordination sphere of such an adatom. By means of Monte Carlo (MC) simulations in the grand canonical ensemble, adsorption isotherms and differential heats of adsorption were calculated. Their striking behaviors were analyzed and discussed in terms of the low temperature phases formed in the system. Finally, the results obtained from MC simulations were compared with the corresponding ones from Bragg-Williams approximation.

Surface phase transitions in one-dimensional channels arranged in a triangular cross-sectional structure: theory and Monte Carlo simulations

Monte Carlo simulations and finite-size scaling analysis have been carried out to study the critical behavior in a submonolayer lattice-gas of interacting monomers adsorbed on one-dimensional channels arranged in a triangular cross-sectional structure. Two kinds of lateral interaction energies have been considered: (1) w(L), interaction energy between nearest-neighbor particles adsorbed along a single channel and (2) w(T), interaction energy between particles adsorbed across nearest-neighbor channels. We focus on the case of repulsive transverse interactions (w(T)>0), where a rich variety of structural orderings are observed in the adlayer, depending on the value of the parameters k(B)Tw(T) (being k(B) the Boltzmann constant) and w(L)w(T). For w(L)w(T)=0, successive planes are uncorrelated, the system is equivalent to the triangular lattice, and the well-known ([square root] 3 x [square root] 3) [([square root] 3 x ([square root] 3)(*)] ordered phase is found at low temperatures and a coverage, theta, of 13. In the more general case (w(L)/w(T) not equal 0), a competition between interactions along a single channel and a transverse coupling between sites in neighboring channels leads to a three-dimensional adsorbed layer. Consequently, the ([square root] 3 x ([square root] 3) and (([square root] 3 x ([square root] 3)(*) structures "propagate" along the channels and new ordered phases appear in the adlayer. Each ordered phase is separated from the disordered state by a continuous order-disorder phase transition occurring at a critical temperature, T(c), which presents an interesting dependence with w(L)/w(T). The Monte Carlo technique was combined with the recently reported free energy minimization criterion approach (FEMCA) [F. Roma et al., Phys. Rev. B 68, 205407 (2003)] to predict the critical temperatures of the order-disorder transformation. The excellent qualitative agreement between simulated data and FEMCA results allows us to interpret the physical meaning of the mechanisms underlying the observed transitions.

Configurational entropy of interacting particles adsorbed on one-dimensional channels arranged in a triangular structure

The configurational entropy of interacting particles adsorbed on one-dimensional channels arranged in a triangular cross-sectional structure is studied by combining Monte Carlo simulation and thermodynamic integration method. Three different energies have been considered in the adsorption process: (1) epsilono, constant interaction energy between a monomer and an adsorption site; (2) wL, interaction energy between nearest-neighbor particles adsorbed along a single channel, and (3) wT, interaction energy between particles adsorbed across nearest-neighbor channels. Special attention is devoted to the case of repulsive transversal interactions (wT>0), for which a rich variety of ordered phases are observed in the adlayer, depending on the value of the parameters kBT/wT (being kB the Boltzmann constant) and wL/wT. The influence of each ordered structure on the configurational entropy of the adlayer has been analyzed and discussed in the context of the lattice-gas model.