Polarons in reduced cesium tungsten bronzes studied using the DFT + U method