Hello,
The spatial join tool in ArcMap is rather limited. It only looks for intersecting polygons or 'fall completely inside of' polygons:

I want to spatially join a polygon to another polygon based on whether the centroid falls within it. Is there another way to do this? Python script maybe?