Path integral Monte Carlo is an exact simulation method for calculating thermodynamic properties of bosonic systems. Properties such as superfluidity and bose condensation are directly related to multiparticle exchange cycles of individual particle paths. Such calculations of bosonic systems in confined geometries, such as helium and hydrogen on surfaces and in droplets are reviewed.