[FEA] Support binary operations between timezone-aware datetime columns

**Is your feature request related to a problem? Please describe.**
Support for timezone-aware datetime columns was largely added in #12813. However, one item remains incomplete: support for binary operations between such columns, or between these columns and timedelta scalars.

**Describe the solution you'd like**
We should work on adding that support.