The map. Define π:GLn→PD, π(A)=AAT. It is continuous, and for U∈O(n),
π(AU)=AUUTAT=AAT=π(A),
so π is constant on left cosets A⋅O(n) and induces a continuous map πˉ:GLn/O(n)→PD. (Note: with A↦ATA one would get the space of right cosets O(n)\GLn.)
Bijective. *Surjective:* for P∈PD, A=P1/2 satisfies AAT=P. *Injective:* if AAT=BBT, put U=B−1A; then UUT=B−1AATB−T=B−1BBTB−T=I, so U∈O(n) and A=BU∈B⋅O(n).
Continuous inverse. The inverse is P↦[P1/2], the composition of the continuous map P↦P1/2=exp(21logP) and the quotient map GLn→GLn/O(n). So πˉ is a homeomorphism.
The case of SLn. Restricting, SLn/SO(n)→{P∈PD:detP=1} is a continuous bijection: if A∈SLn then detAAT=1; conversely, for detP=1 the square root P1/2 has determinant 1 and lies in SLn. Also AAT=BBT with A,B∈SLn gives U=B−1A∈O(n) with detU=1, so U∈SO(n).
Identification with a vector space. By deteS=etrS and the bijection exp:Symn→PD, deteS=1 iff trS=0. So exp restricts to a homeomorphism
{S∈Symn:trS=0} → {P∈PD:detP=1},
and the space of traceless symmetric matrices is a vector space of dimension 2n(n+1)−1.
Check. n=2: SL2/SO(2)≅R2, the hyperbolic plane (as a topological space an open disc).
Remark. Another consequence: GLn deformation retracts onto O(n) and SLn onto SO(n), so topologically SLn≃SO(n).