Hi Xavier,
In CRYSTAL23, symmetry has been introduced as a default option for the integrals calculation for response properties.
But there are some problems, as in your calculation.
We can solve the problem by inserting the keyword NOINTSYM in the INTCPHF block of your input.
Hope the calculation is a success!