Understanding the Adsorption Mechanism of Phenol and Para-Chlorophenol onto Sepiolite Clay: A Combined DFT Calculations, Molecular Dynamics Simulations, and Isotherm Analysis

This study integrates molecular dynamics (MD) simulations and density functional theory (DFT) computations to elucidate the unique adsorption characteristics of phenol and para-chlorophenol onto sepiolite by examining structural deformation, electronic properties, and adsorption energetics. The hydroxyl group (-OH) of phenol mainly determines its adsorption process since it has a quite negative Mulliken charge (−0.428) and significant electrophilic reactivity (fi+ = 0.090), therefore enabling strong hydrogen bonding with the silanol (-SiOH) groups of sepiolite. By π-π interactions with the electron-rich siloxane (-Si-O-Si-) surfaces, the aromatic carbons in phenol improve stability. The close molecular structure allows minimum deformation energy (Edef = 94.18 kcal/mol), hence optimizing alignment with the sepiolite surface. The much negative adsorption energy (Eads = −349.26 kcal/mol) of phenol supports its further thermodynamic stability. Conversely, because of its copious chlorine (-Cl) component, para-chlorophenol runs against steric and electrical obstacles. The virtually neutral Mulliken charge (−0.020) limits electrostatic interactions even if the chlorine atom shows great electrophilicity (fi+ = 0.278). Chlorine’s electron-withdrawing action lowers the hydroxyl group’s (fi+ = 0.077) reactivity, hence lowering hydrogen bonding. Moreover, para-chlorophenol shows strong deformation energy (Edef = 102.33 kcal/mol), which causes poor alignment and less access to high-affinity sites. With less negative than phenol, the adsorption energy for para-chlorophenol (Eads = −317.53 kcal/mol) indicates its reduced thermodynamic affinity. Although more evident in para-chlorophenol because of the polarizable chlorine atom, van der Waals interactions do not balance its steric hindrance and reduced electrostatic interactions. With a maximum Qmax = 0.78 mmol/g, isotherm models confirm the remarkable adsorption capability of phenol in contrast to Qmax = 0.66 mmol/g for para-chlorophenol. By hydrogen bonding and π-cation interactions, phenol builds a dense and structured adsorption layer, and para-chlorophenol shows a chaotic organization with reduced site use. Supported by computational approaches and experimental validation, the results provide a comprehensive knowledge of adsorption mechanisms and provide a basis for the design of adsorbents catered for particular organic pollutants.