One of the easier methods of wastewater treatment is adsorption due to its simplicity in implementation, environmental friendliness, and economic feasibility. Polyvinyl alcohol (PVA) looks promising as an adsorbent due to its biocompatible, non-toxic, water-soluble and eco-friendly nature. The investigation of PVA for its potential in the adsorption of pollutants has been reported in many studies but the mechanistic understanding of the adsorption is poor. The present study used a theoretical approach through density functional theory and Monte Carlo with molecular dynamics simulations to investigate the adsorption mechanism behaviors of model organic molecules (bromothymol blue (BTB), methylene blue (MB), metronidazole (MNZ) and tetracycline (TC)) on PVA surface. The quantum chemical calculations result showed that with the increase in PVA chains (2, 4, 8, 16, and 32 units), the zero-point energy decreases (from −308.79 to −4922.93 kcal/mol) while the dipole moment increases (from 4.37 to 87.52 Debye). Temperature effect on the PVA chain structures showed the same trends for all the chain units and with the increase in temperature (50–600 K), there are no appreciable changes in zero-point energy, enthalpy energy increases while Gibbs free energy decreases. Considering PVA-pollutant complexes, the effects of temperature on the structures showed that there are no appreciable changes in the zero-point energy, Gibbs free and thermal energies increase with an increase in temperature while the kinetic rate of reactions decreases with an increase in temperature. The enthalpy of the reaction showed different trends with antibiotic and dye complexes. In all the thermodynamic properties investigated and the rate of reaction, the order of affinity of the pollutants with PVA followed TC > MNZ > MB > BTB. Monte Carlo simulation was used to investigate the adsorption behavior of the pollutants on the surface of PVA. The negative adsorption energies (−366.56 to −2266.81 kcal/mol) in terms of affinity towards the pollutants on the surface of PVA followed the sequence TC > MNZ > BTB > MB and the molecular dynamic simulation results followed the same order. The obtained results give valuable insights into the mechanism and performance of PVA as an adsorbent. Most of these computational observations are in good agreement with the available experimental results.