An issue that is characteristic of chiroptical simulations is the choice of gauge for the representation of the dipole (and quadrupole) operator. The representation in terms of the position operator is called the length gauge (LG). This is the natural choice of standard absorption simulations to calculate the oscillator strength. However, since the electric dipole interacts with the magnetic dipole (and electric quadrupole) in the rotatory strength, the length gauge formulation leads to values that are origin dependent. In OR calculations, this is commonly solved by using basis functions that are explicitly dependent on the magnetic field, known as London orbitals or gauge including atomic orbitals (GIAOs). These can be used for transition properties as well, but they lead to cumbersome expressions. Therefore, most programs report the rotatory strength in the velocity gauge (VG), where the electric dipole is represented by the momentum operator. The two choices of gauge are numerically equivalent only for exact methods or in the complete basis set limit for variational methods like DFT (although we showed that gauge invariance with approximate functionals does not imply accuracy, as different functionals converge to different results ). Recently, we have proposed a strategy to overcome the origin-dependence issue of the LG rotatory strength without using GIAOs; a method we called origin-invariant LG or LG(OI). Both the VG and LG(OI) approaches were implemented for the OR calculations with DFT-PBC and can be easily extended to TDDFT-PBC simulations of ECD and CPL spectra.