Sorption of water in poly(e-caprolactone) (PCL), with specific focus on the hydrogen-bonding interactions, has been analyzed by combining ab initio calculations, macroscopic thermodynamics modeling, and relevant features emerging from spectroscopic and gravimetric measurements. Fourier transform infrared (FTIR) data, analyzed by difference spectroscopy, two-dimensional correlation spectroscopy, and least-squares curve-fitting analysis associated with gravimetric determination of water sorption isotherm provided information on the system's behavior and on the molecular interactions established between the polymer and the penetrant. A consistent physical picture emerged pointing to the presence of two spectroscopically discernible water species (first-shell and second-shell layers) that have been quantified. Water molecules are present in the form of dimers within the polymer equilibrated with water vapor up to a relative humidity of 0.65. At higher humidities, clustering of water sorbed molecules starts to take place. The multicomponent nu(OH) band representative of absorbed water has been interpreted with the aid of ab initio calculations performed on suitably chosen model systems. The outcomes of spectroscopic analyses were interpreted at a macroscopic level by modeling the thermodynamics of water sorption in PCL based on a nonrandom compressible lattice theory accounting for hydrogen-bonding (HB) interactions. Starting from the fitting of the gravimetric sorption isotherm, the model provided quantitative estimates for the amount of self- and cross-HBs which compare favorably with the FTIR results.