_i {{\rm d}}x^i {{\rm d}}\eta
 + (\delta_{i j} + 2 C_{i j}) {{\rm d}}x^i {{\rm d}}x^j \Bigr] \,,
 \end{gathered}$$ where $C_{i j}$ can be further split into 2 scalar, 2 vector and 2 tensor degrees of freedom. By a direct comparison of these two parametrisations, the relation between the two metric parametrisations is

$$\begin{aligned}
 g_{0 0} &: 1 + 2 \alpha + \alpha^2 - \beta_i \beta^i = 1 + 2 A\,,\\
 g_{0 i} &: \beta_i + 2 h_{i j} \beta^j = B_i\,, \\
 g_{i j} &: \delta_{i j} + 2 h_{i j} = \delta_{i j} + 2 C_{i j}\,.\end{aligned}$$

At first order we obtain $$\begin{gathered}
 A^{(1)} = \alpha^{(1)} \,, \qquad B_i^{(1)} = \beta^{(1)}_i \,, \qquad
 C_{i j}^{(1)} = h_{i j}^{(1)} \,,\end{gathered}$$ and the two parametrisations are the same. However, at second order we get the relations $$\begin{gathered}
 A^{(2)}=  \alpha^{(2)} + \alpha^{(1) 2} - \beta^{(1)}_i \beta^{(1) i}  \,, \qquad
 B_i^{(2)} = \beta^{(2)}_i + 4 h_{i j}^{(1)} \beta^{(1) j}  \,,\qquad
 C_{i j}^{(2)} = h_{i j}^{(2)} \,.\end{gathered}$$

Construction of the distribution function for polarised light {#app:Dist}
=============================================================

We consider a two-dimensional polarisation plane defined by two unit complex vectors $\hat {{\boldsymbol{\epsilon}}}_{({{\mathrm{I}}})}$ and $\hat {{\boldsymbol{\epsilon}}}_{({{\mathrm{II}}})}$, which are mutually orthogonal, $\hat \epsilon_{({{\mathrm{A}}})}{}^{\star\mu} \hat \epsilon_{({{\mathrm{B}}})}{}^\nu g_{\mu\nu}
 =\delta_{({{\mathrm{A}}})({{\mathrm{B}}})}$.
_ i { { \rm d}}x^i { { \rm d}}\eta 
  + (\delta_{i j } + 2 C_{i j }) { { \rm d}}x^i { { \rm d}}x^j \Bigr ] \, , 
  \end{gathered}$$ where $ C_{i j}$ can be further split into 2 scalar, 2 vector and 2 tensor degrees of freedom. By a lineal comparison of these two parametrisations, the relation back between the two metric parametrisations is 

 $ $ \begin{aligned } 
  g_{0 0 } &: 1 + 2 \alpha + \alpha^2 - \beta_i \beta^i = 1 + 2 A\,,\\ 
  g_{0 i } &: \beta_i + 2 h_{i j } \beta^j = B_i\, , \\ 
  g_{i j } &: \delta_{i j } + 2 h_{i j } = \delta_{i j } + 2 C_{i j}\,.\end{aligned}$$ 

 At inaugural orderliness we obtain $ $ \begin{gathered } 
  A^{(1) } = \alpha^{(1) } \, , \qquad B_i^{(1) } = \beta^{(1)}_i \, , \qquad 
  C_{i j}^{(1) } = h_{i j}^{(1) } \,,\end{gathered}$$ and the two parametrisations are the same. However, at second club we drive the relations $ $ \begin{gathered } 
  A^{(2)}=   \alpha^{(2) } + \alpha^{(1) 2 } - \beta^{(1)}_i \beta^{(1) i }   \, , \qquad 
  B_i^{(2) } = \beta^{(2)}_i + 4 h_{i j}^{(1) } \beta^{(1) j }   \,,\qquad 
  C_{i j}^{(2) } = h_{i j}^{(2) } \,.\end{gathered}$$ 

 structure of the distribution affair for polarised lighter { # app: Dist } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 We consider a two - dimensional polarisation airplane defined by two unit complex vectors $ \hat { { \boldsymbol{\epsilon}}}_{({{\mathrm{I}}})}$ and $ \hat { { \boldsymbol{\epsilon}}}_{({{\mathrm{II}}})}$, which are mutually orthogonal, $ \hat \epsilon_{({{\mathrm{A}}})}{}^{\star\mu } \hat \epsilon_{({{\mathrm{B}}})}{}^\nu g_{\mu\nu } 
  = \delta_{({{\mathrm{A}}})({{\mathrm{B}}})}$.
_i {{\rl d}}x^i {{\rm d}}\eta
 + (\delta_{i j} + 2 C_{i j}) {{\rm d}}x^i {{\rm d}}x^j \Bmgr] \,,
 \ens{gatherea}$$ where $C_{i j}$ can be further dpoit ibto 2 scalar, 2 vector ana 2 tensor degrees of hreedom. By a dirxdt comparison lf tkewe two paramettisations, tha relation betfedn the two metric parametrisations is

$$\fegin{alogjed}
 g_{0 0} &: 1 + 2 \alpra + \sjpha^2 - \beta_i \beta^i = 1 + 2 A\,,\\
 g_{0 i} &: \beta_i + 2 h_{j j} \bete^j = B_i\,, \\
 g_{i j} &: \drlta_{i j} + 2 h_{i j} = \delta_{i j} + 2 C_{i u}\,.\end{aligned}$$

At firsh order we ibtayb $$\begin{gathefed}
 A^{(1)} = \alpha^{(1)} \,, \qquad B_i^{(1)} = \beta^{(1)}_i \,, \qquad
 C_{i j}^{(1)} = h_{i j}^{(1)} \,,\end{gaghereb}$$ and the tqo pagdmetrisatiois are the same. However, at second order we get bhe rxlatuons $$\begin{gathered}
 A^{(2)}=  \alpha^{(2)} + \alpha^{(1) 2} - \beja^{(1)}_i \beta^{(1) i}  \,, \qquad
 B_i^{(2)} = \beta^{(2)}_i + 4 y_{i j}^{(1)} \teta^{(1) j}  \,,\dwuaa
 C_{j o}^{(2)} = g_{i j}^{(2)} \,.\ejd{gethered}$$

Consfruction of the distribution finsnoon for poladised jidht {#app:Dist}
=============================================================

We consider a two-dimensional pomarisation plane definee by two unit complex vectors $\rat {{\boldsymbol{\epsilon}}}_{({{\mathrm{I}}})}$ and $\hat {{\boldsymbol{\epvilon}}}_{({{\jxthxn{II}}})}$, wfucj are mutually orthogonal, $\hat \epsilon_{({{\mathrm{A}}})}{}^{\sezr\ku} \hat \epsilon_{({{\mabhrm{B}}})}{}^\nu g_{\mu\nu}
 =\delts_{({{\mwtntm{A}}})({{\mathrm{B}}})}$.
_i {{\rm d}}x^i {{\rm d}}\eta + (\delta_{i 2 j}) {{\rm {{\rm d}}x^j \Bigr] can further split into scalar, 2 vector 2 tensor degrees of freedom. By direct comparison of these two parametrisations, the relation between the two metric parametrisations $$\begin{aligned} g_{0 0} &: 1 + 2 \alpha + \alpha^2 - \beta_i \beta^i 1 2 g_{0 &: \beta_i + 2 h_{i j} \beta^j = B_i\,, \\ g_{i j} &: \delta_{i j} + h_{i j} = \delta_{i j} + 2 C_{i At first order we $$\begin{gathered} A^{(1)} = \alpha^{(1)} \,, B_i^{(1)} \beta^{(1)}_i \,, C_{i = j}^{(1)} \,,\end{gathered}$$ and two parametrisations are the same. However, at second order we get the relations $$\begin{gathered} A^{(2)}= \alpha^{(2)} + 2} - i} \,, B_i^{(2)} \beta^{(2)}_i 4 h_{i j}^{(1)} \,,\qquad C_{i j}^{(2)} = h_{i j}^{(2)} the distribution function for polarised light {#app:Dist} ============================================================= consider a polarisation plane defined by two unit vectors $\hat {{\boldsymbol{\epsilon}}}_{({{\mathrm{I}}})}$ and $\hat {{\boldsymbol{\epsilon}}}_{({{\mathrm{II}}})}$, which are orthogonal, $\hat \epsilon_{({{\mathrm{A}}})}{}^{\star\mu} \hat \epsilon_{({{\mathrm{B}}})}{}^\nu g_{\mu\nu} =\delta_{({{\mathrm{A}}})({{\mathrm{B}}})}$.
_i {{\rm d}}x^i {{\rm d}}\eta
 + (\delta_{i j} + 2 C_{i j}) {{\rm d}}X^i {{\rm d}}x^j \BigR] \,,
 \end{gAthEreD}$$ wHere $c_{i j}$ cAn be further splIT intO 2 scalar, 2 vector and 2 tensor DegreEs OF freEDoM. By a dIrect coMPaRISon Of ThEse TwO PaRametRisAtions, tHe relation BetWeEn the two metrIC pArametrisaTioNs is

$$\begin{aliGneD}
 g_{0 0} &: 1 + 2 \alphA + \aLphA^2 - \Beta_i \BetA^i = 1 + 2 A\,,\\
 g_{0 i} &: \Beta_i + 2 h_{I J} \beta^j = b_i\,, \\
 g_{i j} &: \deltA_{i J} + 2 H_{i j} = \delTA_{i j} + 2 C_{i j}\,.\eND{AlIgneD}$$

At first order we obTAiN $$\Begin{gathered}
 A^{(1)} = \Alpha^{(1)} \,, \qQuAD B_I^{(1)} = \BEta^{(1)}_I \,, \qqUad
 C_{i j}^{(1)} = h_{i j}^{(1)} \,,\eNd{GatheREd}$$ and thE TwO PARamETrisations are The same. HoweVEr, aT seconD oRdeR We get tHe relAtIOns $$\Begin{gatherEd}
 A^{(2)}=  \aLpha^{(2)} + \alpha^{(1) 2} - \Beta^{(1)}_i \bETa^{(1) i}  \,, \qquaD
 b_i^{(2)} = \beta^{(2)}_i + 4 H_{i j}^{(1)} \betA^{(1) j}  \,,\qQuaD
 C_{i j}^{(2)} = H_{I j}^{(2)} \,.\EnD{gaThERed}$$

cOnStrUCtiOn of the dIsTrIbutiOn fuNCTIOn foR poLariSed liGht {#app:Dist}
=============================================================

We cOnsIder A Two-DimenSionaL polArIsatiOn planE defiNeD by two unit complEx veCtors $\hat {{\bOldSyMboL{\ePsiloN}}}_{({{\Mathrm{i}}})}$ anD $\haT {{\boldsyMbol{\epsILon}}}_{({{\MaTHRM{Ii}}})}$, which are mutually oRtHOGoNal, $\hat \epSilon_{({{\mAThRm{a}}})}{}^{\Star\mu} \haT \ePsiLon_{({{\mATHrm{B}}})}{}^\nU g_{\mu\NU}
 =\dElta_{({{\mathRm{A}}})({{\matHRm{b}}})}$.
_i {{\rm d}}x^i {{\rm d}}\ eta
 + (\d elta_ {ij}+2 C_ {i j }) {{\rm d}}x^ i  {{\ rm d}}x^j \Bigr] \,,
\end{ ga t here d }$ $ whe re $C_{ i  j } $  ca nbe  fu rt h er  spli t i nto 2 s calar, 2 v ect or  and 2 tenso r  d egrees offre edom. By a d ire ct com pa ris o n ofthe se tw o para m etrisa tions, th er elatio n  betwee n th e tw o metric parametr i sa t ions is

$$\be gin{al ig n ed }  g_ {00} &: 1 +2\alph a  + \alp h a^ 2 -  \b e ta_i \beta^i= 1 + 2 A\, , \\ g_{0i}  &: \beta_ i + 2  h _ {ij} \beta^j= B_ i\,, \\
g_{i j }  &: \de l ta_{i j } + 2h_{ i j } =\ de lt a_{ ij } + 2C_{ i  j} \,.\end{ al ig ned}$ $

A t f i rstord er w e obt ain $$\begin{ gat here d }
A^{(1 )} =\alp ha ^{(1) } \,,\qqua dB_i^{(1)} = \be ta^{ (1)}_i \, , \ qq uad 
C_{ij }^{(1) } =  h_ {i j}^{ (1)} \, , \en d{ g a t he red}$$ and the two  p a r am etrisati ons ar e  t he same. Ho we ver , at s econd  ord e rwe get t he rel a ti on s $$\be gi n{gath er ed} 
 A ^{(2) } =  \ alpha^ {(2)} +\alph a ^{(1) 2} - \be t a^{(1)}_i \be t a^ { ( 1) i} \,,  \qquad
 B_ i^{( 2 )} =  \be t a^ {(2 ) }_i +  4 h_ {i j} ^ {(1)} \beta^{(1) j}  \,,\qq uad
C_{i j}^{(2)}  = h_{i j} ^ { ( 2)} \,.\ end{ g at h ered}$$

Const ructi on of thed istribut ion f unctionfor polar i s ed light  {# app :Di st} = == ============= = = ==== == ======= === ======= === === === === == ======

W e consid er  a  t wo -di mensi o nal pola ri sat io n p laned efined  by t wo u ni tc omp lex vec t or s $\ha t{{ \bol dsy mb ol{\e psil o n}} }_{({{\ mathrm{I} }}) } $ an d$\ hat {{\ boldsymbol{\e ps ilon}}}_{( {{ \ma thrm{I I } }})}$, w hich are mutually ortho g onal, $ \ha t \ep silo n_{({{\ma thr m{A}}} )}{ } ^{\sta r\mu}\hat\e psi l o n_{({ { \ ma thr m{ B}}})}{}^\ n u  g_ {\mu\ nu }
 = \delta_ {({{\mathrm{A}}})( { {\m athrm{B}}})}$ .
_i {{\rm_d}}x^i {{\rm_d}}\eta
 + (\delta_{i j}_+ 2_C_{i_j}) {{\rm_d}}x^i_{{\rm d}}x^j \Bigr]_\,,
 \end{gathered}$$ where_$C_{i j}$ can be_further split into_2_scalar, 2 vector and 2 tensor degrees of freedom. By a direct comparison of_these_two parametrisations,_the_relation_between the two metric parametrisations_is

$$\begin{aligned}
 g_{0 0} &: 1_+ 2_\alpha + \alpha^2 - \beta_i \beta^i = 1_+_2 A\,,\\
 g_{0_i} &: \beta_i + 2 h_{i j} \beta^j =_B_i\,, \\
 g_{i j} &: \delta_{i_j} + 2_h_{i_j}_= \delta_{i j} +_2 C_{i j}\,.\end{aligned}$$

At first order we_obtain $$\begin{gathered}
 A^{(1)} = \alpha^{(1)} \,,_\qquad B_i^{(1)} = \beta^{(1)}_i \,, \qquad
 C_{i_j}^{(1)} = h_{i j}^{(1)} \,,\end{gathered}$$ and_the two parametrisations are the_same. However,_at second order we get_the relations $$\begin{gathered}
_A^{(2)}= _\alpha^{(2)} + \alpha^{(1)_2} - \beta^{(1)}_i \beta^{(1) i} _\,, \qquad
 B_i^{(2)}_= \beta^{(2)}_i + 4 h_{i j}^{(1)}_\beta^{(1)_j}  \,,\qquad
_C_{i_j}^{(2)}_= h_{i_j}^{(2)} \,.\end{gathered}$$

Construction of_the_distribution function_for_polarised light {#app:Dist}
=============================================================

We consider a two-dimensional_polarisation_plane defined by two unit complex vectors_$\hat {{\boldsymbol{\epsilon}}}_{({{\mathrm{I}}})}$ and $\hat_{{\boldsymbol{\epsilon}}}_{({{\mathrm{II}}})}$,_which are mutually orthogonal,_$\hat \epsilon_{({{\mathrm{A}}})}{}^{\star\mu} \hat \epsilon_{({{\mathrm{B}}})}{}^\nu g_{\mu\nu}
_=\delta_{({{\mathrm{A}}})({{\mathrm{B}}})}$.