Adsorption of Ternary Mixtures in the Presence of Multisite Occupancy: Theory and Monte Carlo Simulations

Adsorption of multicomponent mixtures on solid substrates is essential to numerous technological processes and provides key insights into surface phenomena. Despite advancements in theoretical modeling, many approaches still assume that each adsorbate occupies a single site, thereby neglecting important effects arising from molecules that span multiple adsorption sites. In this work, we broaden the theoretical description of such systems by considering the adsorption of j distinct polyatomic species on triangular lattices. Our approach is based on (i) exact thermodynamic results for polyatomic gases on one-dimensional lattices, extended here to account for substrates with higher coordination numbers, and (ii) the "0D cavity" functional theory originally developed by Lafuente and Cuesta, which reduces to the well-known Guggenheim-DiMarzio model in the limit of rigid rods. As a case study, we explore the behavior of a three-component system consisting of dimers, linear trimers, and triangular trimers adsorbing onto a triangular lattice. This model captures the interplay between structural simplicity, multisite occupancy, configurational diversity, and competition for space, key factors in many practical scenarios involving size-asymmetric molecules. We characterize the system using total and partial isotherms, energy of adsorption, and configurational entropy of the adsorbed phase. To ensure the reliability of our theoretical predictions, we perform Monte Carlo simulations, which show excellent agreement with the analytical approaches. Our findings demonstrate that even complex adsorption systems can be efficiently described using this generalized framework, offering new insights into multicomponent surface adsorption.