I am currently facing a Sylvester equation

$AX+XB = C$

where $A$, $B$, $C$ are all symmetric and a special constraint here is that $X$ should be orthogonal. The Sylvester equation itself may not hold perfectly and hence the least-square solution is also ok.