Krystalloptikk

Krystalloptikk  omhandler optikk i anisotrope materialer. Vanligvis betyr det læren om hvordan lys går gjennom og brytes i krystaller,  men kan også mer generelt beskrive utbredelse av elektromagnetisk stråling i andre medier. Dette kan være flytende krystaller eller materialer som er gjort anisotrope ved pålagte,  mekaniske spenninger eller elektromagnetiske felt. Slike fenomen omtales vanligvis som dobbeltbrytning.
I alminnelighet vil lys i et anisotropt medium ha en brytningsindeks som er avhengig av både dets retning og  dets polarisasjon. En upolarisert lysstråle vil generelt splittes opp i to som hver er polarisert vinkelrett på hverandre. Mens det i en  enakset krystall finnes en spesiell retning som ikke gir noen slik dobbeltbrytning, er det to slike spesielle retninger i en toakset krystall.
Krystalloptikk ble eksperimentelt og teoretisk utforsket i Frankrike på begynnelsen av 1800-tallet. Spesielt var det arbeidene til Augustin Fresnel som la grunnlaget for den beskrivelse av denne delen av optikken som fremdeles gjelder i dag. Dette skjedde på tross av at man på den tiden ikke hadde noen forståelse av lys som et elektromagnetisk fenomen, men bygde på forestillingen at det bestod av svingninger i en mekanisk eter.


== Indeksellipsoiden ==
Det elektriske feltet E og er forbundet med forskyvningsfeltet D gjennom permittiviteten ε ved den lineære sammenhengen D = ε E. I isotrope materialer er denne den samme i alle retninger og har derfor kun en verdi. Men i et anisotropt materiale er vektorene E og D ikke lenger parallelle. Det tilsvarer at permittiviteten har forskjellig verdier i ulike retninger slik at hver komponent av forskyvningsfeltet er en lineærkombinasjon av de elektriske feltkomponentene,

  
    
      
        
          D
          
            i
          
        
        =
        
          ε
          
            i
            x
          
        
        
          E
          
            x
          
        
        +
        
          ε
          
            i
            y
          
        
        
          E
          
            y
          
        
        +
        
          ε
          
            i
            z
          
        
        
          E
          
            z
          
        
      
    
    {\displaystyle D_{i}=\varepsilon _{ix}E_{x}+\varepsilon _{iy}E_{y}+\varepsilon _{iz}E_{z}}
  Da både D og E er vektorfelt, må størrelsene εij  utgjøre komponentene til en tensor. Denne kan representeres ved en 3×3 matrise. Bevarelse av elektromagnetisk energi uttrykt ved Poyntings vektor medfører at denne matrisen er symmetrisk i den forstand at εij  = εji . Det betyr at man kan finne et annet, kartesisk koordinatsystem hvor matrisen tar den diagonale formen

  
    
      
        
          
            [
            
              
                
                  
                    ε
                    
                      x
                      x
                    
                  
                
                
                  
                    ε
                    
                      x
                      y
                    
                  
                
                
                  
                    ε
                    
                      x
                      z
                    
                  
                
              
              
                
                  
                    ε
                    
                      y
                      x
                    
                  
                
                
                  
                    ε
                    
                      y
                      y
                    
                  
                
                
                  
                    ε
                    
                      y
                      z
                    
                  
                
              
              
                
                  
                    ε
                    
                      z
                      x
                    
                  
                
                
                  
                    ε
                    
                      z
                      y
                    
                  
                
                
                  
                    ε
                    
                      z
                      z
                    
                  
                
              
            
            ]
          
        
        ⟹
        
          
            [
            
              
                
                  
                    ε
                    
                      x
                    
                  
                
                
                  0
                
                
                  0
                
              
              
                
                  0
                
                
                  
                    ε
                    
                      y
                    
                  
                
                
                  0
                
              
              
                
                  0
                
                
                  0
                
                
                  
                    ε
                    
                      z
                    
                  
                
              
            
            ]
          
        
      
    
    {\displaystyle {\begin{bmatrix}\varepsilon _{xx}&\varepsilon _{xy}&\varepsilon _{xz}\\\varepsilon _{yx}&\varepsilon _{yy}&\varepsilon _{yz}\\\varepsilon _{zx}&\varepsilon _{zy}&\varepsilon _{zz}\end{bmatrix}}\Longrightarrow {\begin{bmatrix}\varepsilon _{x}&0&0\\0&\varepsilon _{y}&0\\0&0&\varepsilon _{z}\end{bmatrix}}}
  I dette spesielle koordinatsystemet har derfor permittiviteten kun tre verdier tilsvarende de tre hovedretningene til materialet. Komponentene til forskyvningsfeltet blir dermed Dx = εx Ex,  Dy = εy Ey og Dz = εz Ez.

Da permittiviteten i alle medier bestemmer dets brytningsindeks, er det hensiktsmessig å skrive dens komponenter i et anisotropt medium som εk = ε0nk2  hvor ε0 er permittiviteten i vakuum og de dimensjonsløse tallene nk kan betraktes som brytningsindeksene for et anisotropt materiale i de tre hovedretningene.Den elektriske feltenergitettheten som funksjon av dette feltet blir dermed

  
    
      
        
          u
          
            e
            l
          
        
        =
        
          
            1
            2
          
        
        
          E
        
        ⋅
        
          D
        
        =
        
          
            1
            
              2
              
                ε
                
                  0
                
              
            
          
        
        
          
            (
          
        
        
          
            
              D
              
                x
              
              
                2
              
            
            
              n
              
                x
              
              
                2
              
            
          
        
        +
        
          
            
              D
              
                y
              
              
                2
              
            
            
              n
              
                y
              
              
                2
              
            
          
        
        +
        
          
            
              D
              
                z
              
              
                2
              
            
            
              n
              
                z
              
              
                2
              
            
          
        
        
          
            )
          
        
      
    
    {\displaystyle u_{el}={1 \over 2}\mathbf {E} \cdot \mathbf {D} ={1 \over 2\varepsilon _{0}}{\Big (}{D_{x}^{2} \over n_{x}^{2}}+{D_{y}^{2} \over n_{y}^{2}}+{D_{z}^{2} \over n_{z}^{2}}{\Big )}}
  I et kartesisk koordinatsystem med koordinater (x, y, z) proporsjonale med (Dx, Dy, Dz) vil derfor energitettheten være konstant på ellipsoiden

  
    
      
        
          
            
              x
              
                2
              
            
            
              n
              
                x
              
              
                2
              
            
          
        
        +
        
          
            
              y
              
                2
              
            
            
              n
              
                y
              
              
                2
              
            
          
        
        +
        
          
            
              z
              
                2
              
            
            
              n
              
                z
              
              
                2
              
            
          
        
        =
        1
      
    
    {\displaystyle {x^{2} \over n_{x}^{2}}+{y^{2} \over n_{y}^{2}}+{z^{2} \over n_{z}^{2}}=1}
  med halvakser nx, ny og nz. Den kalles for materialets indeksellipsoide og bestemmer dets optiske egenskaper. For noen krystaller er indeksene langs to av aksene like store. Ellipsoiden er da rotasjonssymmetrisk om den tredje aksen som sies å være krystallens optiske akse.  I det generelle tilfellet hvor de tre indeksene alle er forskjellige, vil krystallen ha to optiske akser med spesielle egenskaper for lysets gang.


=== Elektrisk felt som normalvektor ===
På samme måte som man fra ligningen til en ellipse kan man finne dens tangent i hvert punkt, er tangentplanet til indeksellipsoiden i et punkt (x0,y0,z0) gitt som

  
    
      
        
          
            
              x
              
                x
                
                  0
                
              
            
            
              n
              
                x
              
              
                2
              
            
          
        
        +
        
          
            
              y
              
                y
                
                  0
                
              
            
            
              n
              
                y
              
              
                2
              
            
          
        
        +
        
          
            
              z
              
                z
                
                  0
                
              
            
            
              n
              
                z
              
              
                2
              
            
          
        
        =
        1
      
    
    {\displaystyle {xx_{0} \over n_{x}^{2}}+{yy_{0} \over n_{y}^{2}}+{zz_{0} \over n_{z}^{2}}=1}
  Sammenligner man dette med formelen for et generelt plan, ser man at vektoren

  
    
      
        
          N
        
        =
        
          
            (
          
        
        
          
            
              x
              
                0
              
            
            
              n
              
                x
              
              
                2
              
            
          
        
        ,
        
          
            
              y
              
                0
              
            
            
              n
              
                y
              
              
                2
              
            
          
        
        ,
        
          
            
              z
              
                0
              
            
            
              n
              
                z
              
              
                2
              
            
          
        
        
          
            )
          
        
      
    
    {\displaystyle \mathbf {N} ={\Big (}{x_{0} \over n_{x}^{2}},{y_{0} \over n_{y}^{2}},{z_{0} \over n_{z}^{2}}{\Big )}}
  er en normal til ellipsoideflaten. Da koordinatene som benyttes her er proporsjonale med komponentene til forskyvningsfeltet D, betyr det at  normalvektorene N er proporsjonal med det elektriske feltet E. Bortsett fra langs hovedaksene, danner de to vektorene E og D en vinkel med hverandre siden de i allminnelighet ikke er proporsjonale med hverandre i et anisotropt materiale.


=== Fresnel-ellipsoiden ===
Indeksellipsoiden fremkommer fra den elektriske energitettheten i et kartesisk koordinatsystem hvor koordinatene er proporsjonale med komponentene til forskyvningsfeltet D. Man kunne alternativt også uttrykke den samme energitettheten ved komponentene til det elektriske feltet E. Det ville gi opphav til den relaterte ellipsoiden

  
    
      
        
          n
          
            x
          
          
            2
          
        
        
          x
          
            2
          
        
        +
        
          n
          
            y
          
          
            2
          
        
        
          y
          
            2
          
        
        +
        
          n
          
            z
          
          
            2
          
        
        
          z
          
            2
          
        
        =
        1
      
    
    {\displaystyle n_{x}^{2}x^{2}+n_{y}^{2}y^{2}+n_{z}^{2}z^{2}=1}
  med halvakser 1/nx, 1/ny og 1/nz. Den spiller også en viktig rolle for utbredelse av lys i krystaller. Da den  ble først benyttet av Augustin Fresnel, er hans navn knyttet til den for å unngå forveksling med indeksellipsoiden.


== Plane bølger ==

Lys består av elektromagnetiske bølger. De forskjellige feltene som det inneholder, er koblet sammen gjennom de to  Maxwell-ligningene

  
    
      
        
          ∇
        
        ×
        
          E
        
        =
        −
        
          
            
              ∂
              
                B
              
            
            
              ∂
              t
            
          
        
        ,
        
        
        
        
          ∇
        
        ×
        
          H
        
        =
        
          
            
              ∂
              
                D
              
            
            
              ∂
              t
            
          
        
      
    
    {\displaystyle {\boldsymbol {\nabla }}\times \mathbf {E} =-{\frac {\partial \mathbf {B} }{\partial t}},\;\;\;{\boldsymbol {\nabla }}\times \mathbf {H} ={\frac {\partial \mathbf {D} }{\partial t}}}
  Her er de to magnetiske feltene forbundet ved B = μ H via permeabiliteten μ. For at lyset skal kunne fritt utbre seg, må det befinne seg i et dielektrisk materiale som er gjennomsiktig.  Slike stoffer er umagnetiske, og man kan derfor med god nøyaktighet sette deres permeabilitet lik med den i vakuum. Det betyr at μ = μ0 = 4π ⋅10-7 N/A2  som også kalles den «magnetiske konstanten».Det elektriske feltet i en plan bølge med bølgevektor k og  vinkelfrekvens ω  kan skrives som en fasevektor

  
    
      
        
          E
        
        (
        
          r
        
        ,
        t
        )
        =
        
          
            E
          
          
            0
          
        
        
          e
          
            i
            (
            
              k
            
            ⋅
            
              r
            
            −
            ω
            t
            )
          
        
      
    
    {\displaystyle \mathbf {E} (\mathbf {r} ,t)=\mathbf {E} _{0}e^{i(\mathbf {k} \cdot \mathbf {r} -\omega t)}}
  hvis man betrakter det elektriske feltet i bølgen. De andre feltene tilhørende den samme bølgen kan skrives på samme form. Benytter man at en tidsderivasjon ∂/∂t  nå vil tilsvare en multiplikasjon med -iω og nabla-operatoren ∇ gir en tilsvarende faktor i k, reduseres de to Maxwell-ligningene til

  
    
      
        
          k
        
        ×
        
          E
        
        =
        ω
        
          B
        
        ,
        
        
        
        
          k
        
        ×
        
          H
        
        =
        −
        ω
        
          D
        
      
    
    {\displaystyle \mathbf {k} \times \mathbf {E} =\omega \mathbf {B} ,\;\;\;\mathbf {k} \times \mathbf {H} =-\omega \mathbf {D} }
  Med B = μ0H i den første ligningen, finner man det magnetiske feltet H uttrykt ved det elektriske feltet E.  Innsatt i den andre ligningen, får man da resultatet

  
    
      
        (
        
          k
        
        ×
        
          E
        
        )
        ×
        
          k
        
        =
        
          k
          
            2
          
        
        
          E
        
        −
        
          k
        
        (
        
          k
        
        ⋅
        
          E
        
        )
        =
        
          μ
          
            0
          
        
        
          ω
          
            2
          
        
        
          D
        
      
    
    {\displaystyle (\mathbf {k} \times \mathbf {E} )\times \mathbf {k} =k^{2}\mathbf {E} -\mathbf {k} (\mathbf {k} \cdot \mathbf {E} )=\mu _{0}\omega ^{2}\mathbf {D} }
  som kan benyttes til å bestemme de elektriske feltene for bølgen. Den er ekvivalent med den elektromagnetiske bølgeligningen i det anisotrope materialet.


=== Fresnel-ligningen ===
En løsning av bølgeligningen vil gi hvordan vinkelfrekvensen ω varierer med bølgevektoren k til den plane bølgen. Dette gjøres mest direkte i det spesielle koordinatsystemet som har akser langs hovedaksene til krystallen. Komponentene til de to feltene er da forbundet med den enkle relasjonen Dj = ε0nj2Ej  hvor indeksen j tar de tre verdiene x, y og z. Dermed reduseres den sektorielle bølgeligningen til tre vanlige ligninger

  
    
      
        
          k
          
            2
          
        
        
          E
          
            j
          
        
        −
        (
        
          k
        
        ⋅
        
          E
        
        )
        
          k
          
            j
          
        
        =
        
          μ
          
            0
          
        
        
          ω
          
            2
          
        
        
          D
          
            j
          
        
      
    
    {\displaystyle k^{2}E_{j}-(\mathbf {k} \cdot \mathbf {E} )k_{j}=\mu _{0}\omega ^{2}D_{j}}
  De danner et homogent, lineært ligningssystem for det elektriske feltet. Ikke-trivielle løsninger for feltkomponentene vil da bare eksistere når determinanten til koeffisientene i ligningene er null. Det gir en ligning som kalles Fresnel-ligningen. Den kan benyttes til å finne dispersjonsrelasjonen ω =  ω(k).
Denne ligningen kan utledes enklere ved å løse ut komponentene til forskyvningsfeltet fra de tre ligningene. Det gir

  
    
      
        
          D
          
            j
          
        
        =
        
          
            
              
                ε
                
                  0
                
              
              
                k
                
                  j
                
              
              
                n
                
                  j
                
                
                  2
                
              
            
            
              
                k
                
                  2
                
              
              −
              (
              ω
              
                n
                
                  j
                
              
              
                /
              
              c
              
                )
                
                  2
                
              
            
          
        
        (
        
          k
        
        ⋅
        
          E
        
        )
      
    
    {\displaystyle D_{j}={\varepsilon _{0}k_{j}n_{j}^{2} \over k^{2}-(\omega n_{j}/c)^{2}}(\mathbf {k} \cdot \mathbf {E} )}
  etter å ha innført den kvadrerte lyshastigheten c2 = 1/ε0μ0 i vakuum. Men Gauss' lov ∇⋅D = 0 tar nå formen k⋅D = kxDx + kyDy + kzDz = 0. Det betyr at vektoren D står vinkelrett på bølgevektoren k. Dette gjelder derfor også i et anisotropt materiale, men ikke for den elektriske feltvektoren E. På denne måten har man dermed Fresnel-ligningen

  
    
      
        
          
            
              
                k
                
                  x
                
                
                  2
                
              
              
                n
                
                  x
                
                
                  2
                
              
            
            
              
                k
                
                  2
                
              
              −
              (
              ω
              
                n
                
                  x
                
              
              
                /
              
              c
              
                )
                
                  2
                
              
            
          
        
        +
        
          
            
              
                k
                
                  y
                
                
                  2
                
              
              
                n
                
                  y
                
                
                  2
                
              
            
            
              
                k
                
                  2
                
              
              −
              (
              ω
              
                n
                
                  y
                
              
              
                /
              
              c
              
                )
                
                  2
                
              
            
          
        
        +
        
          
            
              
                k
                
                  z
                
                
                  2
                
              
              
                n
                
                  z
                
                
                  2
                
              
            
            
              
                k
                
                  2
                
              
              −
              (
              ω
              
                n
                
                  z
                
              
              
                /
              
              c
              
                )
                
                  2
                
              
            
          
        
        =
        0
      
    
    {\displaystyle {k_{x}^{2}n_{x}^{2} \over k^{2}-(\omega n_{x}/c)^{2}}+{k_{y}^{2}n_{y}^{2} \over k^{2}-(\omega n_{y}/c)^{2}}+{k_{z}^{2}n_{z}^{2} \over k^{2}-(\omega n_{z}/c)^{2}}=0}
  Dette resultatet kan omskrives på forskjellige former. Benytter man at kx2 + ky2 + kz2 = k2, finner man den ekvivalente ligningen

  
    
      
        
          
            
              k
              
                x
              
              
                2
              
            
            
              
                k
                
                  2
                
              
              −
              (
              ω
              
                n
                
                  x
                
              
              
                /
              
              c
              
                )
                
                  2
                
              
            
          
        
        +
        
          
            
              k
              
                y
              
              
                2
              
            
            
              
                k
                
                  2
                
              
              −
              (
              ω
              
                n
                
                  y
                
              
              
                /
              
              c
              
                )
                
                  2
                
              
            
          
        
        +
        
          
            
              k
              
                z
              
              
                2
              
            
            
              
                k
                
                  2
                
              
              −
              (
              ω
              
                n
                
                  z
                
              
              
                /
              
              c
              
                )
                
                  2
                
              
            
          
        
        =
        1
      
    
    {\displaystyle {k_{x}^{2} \over k^{2}-(\omega n_{x}/c)^{2}}+{k_{y}^{2} \over k^{2}-(\omega n_{y}/c)^{2}}+{k_{z}^{2} \over k^{2}-(\omega n_{z}/c)^{2}}=1}
  Når ligningen skrives om som et polynom, fremkommer en andregradsligning for ω2. For hver verdi av bølgevektoren k, vil det derfor være to løsninger ω1 og ω2. Det at fortegnet til vinkelfrekvensen ikke spiller noen rolle, tilsvarer at for hver slik løsning vil det være en tilsvarende for - k.


== Ortogonale løsninger ==

Hver av de to løsningene ω1(k) og ω2(k) av Fresnel-ligningen kan fremstilles som deler av en sammensatt flate i k-rommet med kartesiske koordinater kx, ky og kz. I dette rommet vil en vektor k skjære flaten i to punkt som tilsvarer de to løsningene. Denne flaten har i litteraturen ulike navn, men kan kalles «bølgevektorflaten», «bølgenormalflaten» eller ganske enkelt bølgeflaten.For hver slik bestemmelse av de to tillatte vinkelfrekvensene for gitt verdi av bølgatallet, kan så de elektriske feltene E1 og  E2 bestemmes fra de opprinnelige ligningene. De tilsvarende forskyvningsfeltene  D1 og  D2 er da ortogonale som betyr at de står vinkelrett på hverandre. Det følger fra

  
    
      
        
          
            
              
                
                  
                    D
                  
                  
                    1
                  
                
                ⋅
                
                  
                    D
                  
                  
                    2
                  
                
              
              
                
                =
                
                  ε
                  
                    0
                  
                  
                    2
                  
                
                
                  ∑
                  
                    j
                    =
                    x
                    ,
                    y
                    ,
                    z
                  
                
                
                  
                    
                      
                        k
                        
                          j
                        
                      
                      
                        n
                        
                          j
                        
                        
                          2
                        
                      
                    
                    
                      
                        k
                        
                          2
                        
                      
                      −
                      (
                      
                        ω
                        
                          1
                        
                      
                      
                        n
                        
                          j
                        
                      
                      
                        /
                      
                      c
                      
                        )
                        
                          2
                        
                      
                    
                  
                
                ⋅
                
                  
                    
                      
                        k
                        
                          j
                        
                      
                      
                        n
                        
                          j
                        
                        
                          2
                        
                      
                    
                    
                      
                        k
                        
                          2
                        
                      
                      −
                      (
                      
                        ω
                        
                          2
                        
                      
                      
                        n
                        
                          j
                        
                      
                      
                        /
                      
                      c
                      
                        )
                        
                          2
                        
                      
                    
                  
                
                (
                
                  k
                
                ⋅
                
                  
                    E
                  
                  
                    1
                  
                
                )
                (
                
                  k
                
                ⋅
                
                  
                    E
                  
                  
                    2
                  
                
                )
              
            
            
              
              
                
                =
                
                  
                    
                      
                        ε
                        
                          0
                        
                        
                          2
                        
                      
                      
                        c
                        
                          2
                        
                      
                    
                    
                      
                        ω
                        
                          1
                        
                      
                      −
                      
                        ω
                        
                          2
                        
                      
                    
                  
                
                
                  ∑
                  
                    j
                    =
                    x
                    ,
                    y
                    ,
                    z
                  
                
                
                  [
                  
                    
                      
                        
                          
                            k
                            
                              j
                            
                            
                              2
                            
                          
                          
                            n
                            
                              j
                            
                            
                              2
                            
                          
                        
                        
                          
                            k
                            
                              2
                            
                          
                          −
                          (
                          
                            ω
                            
                              1
                            
                          
                          
                            n
                            
                              j
                            
                          
                          
                            /
                          
                          c
                          
                            )
                            
                              2
                            
                          
                        
                      
                    
                    −
                    
                      
                        
                          
                            k
                            
                              j
                            
                            
                              2
                            
                          
                          
                            n
                            
                              j
                            
                            
                              2
                            
                          
                        
                        
                          
                            k
                            
                              2
                            
                          
                          −
                          (
                          
                            ω
                            
                              2
                            
                          
                          
                            n
                            
                              j
                            
                          
                          
                            /
                          
                          c
                          
                            )
                            
                              2
                            
                          
                        
                      
                    
                  
                  ]
                
                (
                
                  k
                
                ⋅
                
                  
                    E
                  
                  
                    1
                  
                
                )
                (
                
                  k
                
                ⋅
                
                  
                    E
                  
                  
                    2
                  
                
                )
              
            
          
        
      
    
    {\displaystyle {\begin{aligned}\mathbf {D} _{1}\cdot \mathbf {D} _{2}&=\varepsilon _{0}^{2}\sum _{j=x,y,z}{k_{j}n_{j}^{2} \over k^{2}-(\omega _{1}n_{j}/c)^{2}}\cdot {k_{j}n_{j}^{2} \over k^{2}-(\omega _{2}n_{j}/c)^{2}}(\mathbf {k} \cdot \mathbf {E} _{1})(\mathbf {k} \cdot \mathbf {E} _{2})\\&={\varepsilon _{0}^{2}c^{2} \over \omega _{1}-\omega _{2}}\sum _{j=x,y,z}\left[{k_{j}^{2}n_{j}^{2} \over k^{2}-(\omega _{1}n_{j}/c)^{2}}-{k_{j}^{2}n_{j}^{2} \over k^{2}-(\omega _{2}n_{j}/c)^{2}}\right](\mathbf {k} \cdot \mathbf {E} _{1})(\mathbf {k} \cdot \mathbf {E} _{2})\end{aligned}}}
  Men hver av summene her er null på grunn av Fresnels ligning. Dermed har man at D1⋅D2 = 0. De to løsningene representerer derfor bølger som er polariserte i plan som danner 90° med hverandre.
Vinkelfrekvensen ω bestemmer  brytningsindeksen n til bølgen via fasehastigheten u = ω/k = c /n.  Da løsningene ω1,2 avhenger av bølgevektoren k, vil brytningsindeksen vanligvis i et anisotropt materiale derfor være avhengig av retningen til bølgene. De to løsningene sies å være «ekstraordinære» for å skille dem fra ordinære bølger i isotrope medier hvor brytningsindeksen er den samme i alle retninger.


=== Geometrisk konstruksjon ===
De to planpolariserte løsningene av Fresnel-ligningen kan finnes ved en geometrisk konstruksjon basert på indeksellipsoiden. Et plan gjennom sentrum som har en  flatenormal i samme retning som bølgevektoren k, er beskrevet ved ligningen k⋅x = 0.  Da vil skjæringsflaten være en ellipse hvor de to halvaksene peker i de ortogonale retningene tilsvarende vektorene D1 og D2 som angir polarisjonen til bølgene . Samtidig vil lengden av halvaksene være de to tilsvarende brytningsindeksene  n1 og n2.
For en generell krystall med hovedbrytningsindekser nx ≠ ny ≠ nz vil planet i to spesielle retninger føre til at snittellipsen blir en sirkel. Disse retningene definerer krystallens  to «optiske akser». De to mulige bølgene langs disse spesielle retningene vil da kunne forplante seg med samme brytningsindeks. Krystallen sies i dette tilfellet å være toakset.


== Enaksete krystaller ==
I mange krystaller er to av hovedindeksene like store. Det er da vanlig å benytte et koordinatsystem slik at disse tilsvarer x- og y-retningene. En slik krystall vil da bare ha en optisk akse som ved dette koordinatvalget tilsvarer z-aksen. I tillegg kalles nx = ny for den «ordinære» brytningsindeksen  no, mens  nz  er den «ekstraordinære» indeksen  ne. Et kvantitativt uttrykk for størrelsen til dobbeltbrytningen en slik krystall kan frembringe er gitt ved differensen

  
    
      
        Δ
        n
        =
        
          n
          
            e
          
        
        −
        
          n
          
            o
          
        
      
    
    {\displaystyle \Delta n=n_{e}-n_{o}}
  Er denne større enn null, sies krystallen å være postiv. I motsatt fall er den negativ.

Bølgeflaten for en enakset krystall fremkommer fra Fresnel-ligningen som i dette tilfellet kan faktoriseres til å gi

  
    
      
        
          (
          
            
              
                
                  ω
                  
                    2
                  
                
                
                  c
                  
                    2
                  
                
              
            
            −
            
              
                
                  k
                  
                    x
                  
                  
                    2
                  
                
                
                  n
                  
                    o
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  k
                  
                    y
                  
                  
                    2
                  
                
                
                  n
                  
                    o
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  k
                  
                    z
                  
                  
                    2
                  
                
                
                  n
                  
                    o
                  
                  
                    2
                  
                
              
            
          
          )
        
        
          (
          
            
              
                
                  ω
                  
                    2
                  
                
                
                  c
                  
                    2
                  
                
              
            
            −
            
              
                
                  k
                  
                    x
                  
                  
                    2
                  
                
                
                  n
                  
                    e
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  k
                  
                    y
                  
                  
                    2
                  
                
                
                  n
                  
                    e
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  k
                  
                    z
                  
                  
                    2
                  
                
                
                  n
                  
                    o
                  
                  
                    2
                  
                
              
            
          
          )
        
        =
        0
      
    
    {\displaystyle \left({\frac {\omega ^{2}}{c^{2}}}-{\frac {k_{x}^{2}}{n_{o}^{2}}}-{\frac {k_{y}^{2}}{n_{o}^{2}}}-{\frac {k_{z}^{2}}{n_{o}^{2}}}\right)\left({\frac {\omega ^{2}}{c^{2}}}-{\frac {k_{x}^{2}}{n_{e}^{2}}}-{\frac {k_{y}^{2}}{n_{e}^{2}}}-{\frac {k_{z}^{2}}{n_{o}^{2}}}\right)=0}
  Den første faktoren her gir en løsning med brytningsindeks no som er uavhengig av rettning til vektoren k gitt ved komponentene (kx, ky, kz). Denne delen av bølgeflaten er derfor en kule som tilsvarer den ordinære bølgen.
Men i tillegg har bølgeflaten et annet lag som kommer fra den andre faktoren på venstresiden av ligningen. Settes den lik null, har man

  
    
      
        
          
            
              ω
              
                2
              
            
            
              c
              
                2
              
            
          
        
        =
        
          
            
              
                k
                
                  x
                
                
                  2
                
              
              +
              
                k
                
                  y
                
                
                  2
                
              
            
            
              n
              
                e
              
              
                2
              
            
          
        
        +
        
          
            
              k
              
                z
              
              
                2
              
            
            
              n
              
                o
              
              
                2
              
            
          
        
      
    
    {\displaystyle {\omega ^{2} \over c^{2}}={k_{x}^{2}+k_{y}^{2} \over n_{e}^{2}}+{k_{z}^{2} \over n_{o}^{2}}}
  som fremstiller en ny, rotasjonssymmetrisk ellipsoide i k-rommet. For en postiv krystall som har ne > no, er den sammentrykket ved polene (oblat). Derimot er ellipsoiden for en negativ krystall sammentrukket ved ekvator (prolat). Den fulle bølgeflatene er sammensatt av disse to flatene med den ene inni den andre slik at de berører hverandre ved polene. For positive krystaller ligger kuleflaten inni ellipsoiden og omvendt for negative krystaller.
Hvis k-vektoren danner vinkelen θ med z-aksen, er kz = k cosθ, mens kx2 +  ky2 = k2 sin2θ. Ellipsoiden gir derfor opphav til den retningsavhengig  brytningsindeksen

  
    
      
        
          
            1
            
              
                n
                
                  2
                
              
              (
              θ
              )
            
          
        
        =
        
          
            
              
                sin
                
                  2
                
              
              ⁡
              θ
            
            
              n
              
                e
              
              
                2
              
            
          
        
        +
        
          
            
              
                cos
                
                  2
                
              
              ⁡
              θ
            
            
              n
              
                o
              
              
                2
              
            
          
        
      
    
    {\displaystyle {1 \over n^{2}(\theta )}={\sin ^{2}\theta  \over n_{e}^{2}}+{\cos ^{2}\theta  \over n_{o}^{2}}}
  for den ekstraordinære bølgen. Langs den optiske aksen der θ = 0, blir den like stor som for den ordinære bølgen, n = no. Derimot for en retning normalt på den optiske aksen, blir denne ekstraordinære indeksen lik med  ne, mens den ordinære løsningen fortsatt har brytningsindeks no.


=== Polarisasjon ===
Fra den geometriske konstruksjonen av løsningene følger direkte polarisasjonen til de to løsningene. Ordinære bølger vil alltid være polariserte normalt på den optiske aksen. Ligger den langs z-aksen, vil da begge vektorene E og D være i xy-planet. De vil også i dette tilfellet ha samme retning da nx = ny. Samtidig befinner magnetfeltene seg i det ortogonale planet som den optiske aksen danner sammen med bølgevektoren.
Derimot har ekstraordinære bølger sine elektriske vektorer E og D i det planet definert av den optiske aksen og bølgevektoren k. Men i dette tilfellet har E og D en liten vinkel seg i mellom. Det magnetiske feltet er nå i xy-planet slik at vektoren H står vinkelrett både på E og k.


== Poyntings vektor ==
Det elektriske forskyvningsfeltet D er alltid vinkelrett på bølgevektoren k, også i   i en anisotrop krystall. Men da det elektriske feltet E ikke alltid har samme retning som D, vil heller ikke Poyntings vektor   S = E × H generelt peke i samme retning som k, men danne en vinkel med denne som er like stor som vinkelen mellom E og D.
Fra Maxwells to ligninger k × E = ω B  og k × H = -ω D  for en plane bølge med bølgevektoren k, vinkelfrekvensen ω uttrykkes ved Poyntings vektor. Multipliseres den første med H og den andre med E, finner man

  
    
      
        ω
        =
        
          
            
              (
              
                k
              
              ×
              
                E
              
              )
              ⋅
              
                H
              
            
            
              
                B
              
              ⋅
              
                H
              
            
          
        
        =
        −
        
          
            
              (
              
                k
              
              ×
              
                H
              
              )
              ⋅
              
                E
              
            
            
              
                D
              
              ⋅
              
                E
              
            
          
        
      
    
    {\displaystyle \omega ={(\mathbf {k} \times \mathbf {E} )\cdot \mathbf {H}  \over \mathbf {B} \cdot \mathbf {H} }=-{(\mathbf {k} \times \mathbf {H} )\cdot \mathbf {E}  \over \mathbf {D} \cdot \mathbf {E} }}
  Her kan de to trippelproduktene i tellerne skrives om til k⋅(E × H) = k⋅S. Begge nevnerne er derfor like store slik at summen av den elektriske og den magnetiske feltenergitettheten  er ganske enkelt gitt som U = D⋅E. Dermed har man den fundamentale sammenhengen

  
    
      
        ω
        =
        
          
            
              
                k
              
              ⋅
              
                S
              
            
            U
          
        
      
    
    {\displaystyle \omega ={\mathbf {k} \cdot \mathbf {S}  \over U}}
  som selvsagt også gjelder for en elektromagnetisk bølge i et isotropt medium der vektorene k og S er parallelle.
Dette resultatet kan betraktes som en dispersjonsrelasjonen. Den tilhørende gruppehastigheten er dermed

  
    
      
        
          v
        
        =
        
          
            ∇
          
          
            
              k
            
          
        
        ω
        =
        
          
            
              S
            
            U
          
        
      
    
    {\displaystyle \mathbf {v} ={\boldsymbol {\nabla }}_{\mathbf {k} }\omega ={\mathbf {S}  \over U}}
  og er rettet langs Poynting-vektoren som transporterer energien til bølgen. Det er denne retningen som definerer den tilsvarende lysstrålen som kan påvises ved observasjon eller måling. I krystalloptikken blir derfor v vanligvis kalt for strålehastigheten.


=== Spredningsvinkelen ===
Da bølgeflaten ω = konstant kan betraktes som en ekvipotensialflate i k-rommet, betyr det at vektorene v og S står vinkelrett på denne flaten. Da den også står vinkelrett på det elektriske feltet E, er dette en tangentvektor til bølgeflaten. I tillegg står bølgevektoren k vinkelrett på forskyvningsfeltet D slik at  vinkelen mellom k og S den samme som mellom D og E.
Vinkelen omtales noen ganger for «spredningsvinkelen»  da den er et uttrykk for hvor mye den ekstraordinære strålen spredes fra den ordinære ved dobbeltbrytning. Hvis den kalles for α, vil da D⋅E = DE cosα. Størrelsen til Poyntings vektor er S = EH der magnetfeltet er bestemt av Maxwell-ligningen k × H = -ω D  som gir  H = uD  der u = ω/k  er fasehastigheten til bølgen. Dermed er v = EH/DE cosα  slik at

  
    
      
        cos
        ⁡
        α
        =
        
          
            u
            v
          
        
      
    
    {\displaystyle \cos \alpha ={u \over v}}
  Det kan gis et enkelt, geometrisk argument for dette resultatet. Da en plan bølgefront med retning k beveger seg fremover med fasehastigheten u, må lysstrålen til et punkt på denne som ligger en vinkel α til siden, bevege seg med en hastighet v slik at u = v cosα for at fronten skal forbli plan.


=== Enaksete krystaller ===

For en enakset krystall er bølgeflaten for den ordinære strålen kuleformet slik at S peker i samme retning som k. Dermed har den også en strålehastighet som er lik med fasehastigheten, v = u. Den ekstraordinære strålen derimot har en ellipsoideformet bølgeflate som er rotasjonssymmetrisk om den optiske aksen som vanligvis legges langs kz -aksen.  Da vil Poyntings vektor S ikke lenger være parallell med k som danner vinkelen θ med den optiske aksen. Betegner man den tilsvarende vinkelen  for S med χ, vil da differansen α = θ -  χ  være spredningsvinkelen.I ky -planet danner bølgeflaten for den ekstraordinære strålen en ellipse som har en halvakse med lengde ne langs kx -aksen,  mens den andre halvaksen  langs kz -aksen har lengde no. Ved å ta den deriverte av ligningen for ellipsen, følger at  no2kx dkx + ne2kz dkz = 0. Nå er dkz/dkx  lik med stigningstallet for tangenten til ellipsen i punktet (kx, kz) der kx/kz = tanθ. Men vektoren S står vinkelrett på bølgeflaten, det vil si på denne tangenten. Bortsett fra fortegnet er da tangens til vinkelen χ som den danner med  kz -aksen,  derfor like stor som dette stigningstallet,

  
    
      
        tan
        ⁡
        χ
        =
        
          
            
              n
              
                o
              
              
                2
              
            
            
              n
              
                e
              
              
                2
              
            
          
        
        tan
        ⁡
        θ
      
    
    {\displaystyle \tan \chi ={n_{o}^{2} \over n_{e}^{2}}\tan \theta }
  Herav kan man finne spredningsvinkelen α uttrykt ved vinkelen θ for bølgevektoren k ved å benytte den trigonometriske identiteten for tangens til en differans av to vinkler. Den gir

  
    
      
        tan
        ⁡
        α
        =
        
          
            
              tan
              ⁡
              θ
              −
              tan
              ⁡
              χ
            
            
              1
              +
              tan
              ⁡
              θ
              tan
              ⁡
              χ
            
          
        
        =
        
          
            
              (
              
                n
                
                  e
                
                
                  2
                
              
              −
              
                n
                
                  o
                
                
                  2
                
              
              )
              tan
              ⁡
              θ
            
            
              
                n
                
                  e
                
                
                  2
                
              
              +
              
                n
                
                  o
                
                
                  2
                
              
              
                tan
                
                  2
                
              
              ⁡
              θ
            
          
        
      
    
    {\displaystyle \tan \alpha ={\tan \theta -\tan \chi  \over 1+\tan \theta \tan \chi }={(n_{e}^{2}-n_{o}^{2})\tan \theta  \over n_{e}^{2}+n_{o}^{2}\tan ^{2}\theta }}
  og er positiv eller negativ avhengig av differansen Δn = ne - no  som også bestemmer om krystallen er positiv eller negativ. Det er dette som avgjør om vektoren k eller S ligger nærmest den optiske aksen.


== Stråleflaten ==
Spredningsvinkelen bestemmer strålehastighten v når fasehastigheten u er kjent. Men denne vinkelen kan i allminnelighet ikke uten videre beregnes fra hovedbrytningsindeksene nx, ny og nz. Men på samme måte som Fresnel-ligningen med den tilhørende bølgeflaten gir de tillatte løsningene for brytningsindeksene og dermed for fasehastighetene, kan det utledes en analog ligning som bestemmer de tillatte strålehastighetene. Denne nye ligningen har løsninger som kan fremstilles som punkter på en stråleflate som på et vis kan sies å være den «resiproke» bølgeflaten.Denne nye flaten kan igjen utledes ved å bruke de samme to Maxwell-ligningene for en plan bølge med bølgenormal k. Istedenfor å uttrykke forskyvningsfeltet D ved det elektriske feltet E, skrives disse nå om slik at man kan finne E fra D. Det kan gjøres ved å benytte at vektorene S, E og D ligger i samme plan. Generelt kan man da skrive vektoren D som lineærkombinasjonen

  
    
      
        
          D
        
        =
        p
        
          E
        
        +
        q
        
          S
        
      
    
    {\displaystyle \mathbf {D} =p\mathbf {E} +q\mathbf {S} }
  hvor koeffisientene p og q må bestemmes. Ved å benytte at S er vinkelrett på E, har man med en gang at S⋅E = 0 slik at q = S⋅D/S2. Dermed har man den lineære sammenhengen S2 D - (S⋅D)S = pS2 E. Bortsett fra den ennå ukjente faktoren p, har den samme struktur som ligningen

  
    
      
        
          k
          
            2
          
        
        
          E
        
        −
        (
        
          k
        
        ⋅
        
          E
        
        )
        
          k
        
        =
        
          μ
          
            0
          
        
        
          ω
          
            2
          
        
        
          D
        
      
    
    {\displaystyle k^{2}\mathbf {E} -(\mathbf {k} \cdot \mathbf {E} )\mathbf {k} =\mu _{0}\omega ^{2}\mathbf {D} }
  som ble brukt til å utlede Fresnel-ligningen. Ved å multiplisere denne skalart med S og den forrige med k, vil en sammenligning av de to resulterende uttrykkene gi at

  
    
      
        p
        
          μ
          
            0
          
        
        
          ω
          
            2
          
        
        
          S
          
            2
          
        
        =
        (
        
          k
        
        ⋅
        
          S
        
        
          )
          
            2
          
        
        =
        
          k
          
            2
          
        
        
          S
          
            2
          
        
        
          cos
          
            2
          
        
        ⁡
        α
      
    
    {\displaystyle p\mu _{0}\omega ^{2}S^{2}=(\mathbf {k} \cdot \mathbf {S} )^{2}=k^{2}S^{2}\cos ^{2}\alpha }
  Dermed er koeffisienten p bestemt. Resultatet kan forenkles ved å skrive cosinus til spredningsvinkelen som cosα = u/v = ω/kv. Dermed er p = 1/μ0v2  slik at den ønskede, lineære sammenhengen blir

  
    
      
        
          S
          
            2
          
        
        
          D
        
        −
        (
        
          S
        
        ⋅
        
          D
        
        )
        
          S
        
        =
        
          
            
              U
              
                2
              
            
            
              μ
              
                0
              
            
          
        
        
          E
        
      
    
    {\displaystyle S^{2}\mathbf {D} -(\mathbf {S} \cdot \mathbf {D} )\mathbf {S} ={U^{2} \over \mu _{0}}\mathbf {E} }
  hvor  energitettheten U = S/v  opptrer. Igjen har man et lineært ligningssystem for de tre komponentene av det elektriske feltet. Det har bestemte løsninger bare for visse verdier av denne tettheten hvis man antar at strålevektoren S er gitt. På denne måten kan også strålehastigheten v finnes.

Når man benytter koordinatsystemet gitt ved hovedaksene til permittivitetstensoren der den har de diagonale komponentene εk = ε0nk2, kan de tre ligningene skrives på formen

  
    
      
        
          E
          
            k
          
        
        =
        
          
            
              
                μ
                
                  0
                
              
              
                S
                
                  k
                
              
              
                c
                
                  k
                
                
                  2
                
              
            
            
              
                S
                
                  2
                
              
              −
              (
              U
              
                c
                
                  k
                
              
              
                )
                
                  2
                
              
            
          
        
        (
        
          S
        
        ⋅
        
          D
        
        )
      
    
    {\displaystyle E_{k}={\mu _{0}S_{k}c_{k}^{2} \over S^{2}-(Uc_{k})^{2}}(\mathbf {S} \cdot \mathbf {D} )}
  etter å ha innført lyshastighetene ck = c /nk langs de tre aksene. Da S⋅E = 0, har man dermed den ene ligningen

  
    
      
        
          
            
              
                S
                
                  x
                
                
                  2
                
              
              
                c
                
                  x
                
                
                  2
                
              
            
            
              
                S
                
                  2
                
              
              −
              (
              U
              
                c
                
                  x
                
              
              
                )
                
                  2
                
              
            
          
        
        +
        
          
            
              
                S
                
                  y
                
                
                  2
                
              
              
                c
                
                  y
                
                
                  2
                
              
            
            
              
                S
                
                  2
                
              
              −
              (
              U
              
                c
                
                  y
                
              
              
                )
                
                  2
                
              
            
          
        
        +
        
          
            
              
                S
                
                  z
                
                
                  2
                
              
              
                c
                
                  z
                
                
                  2
                
              
            
            
              
                S
                
                  2
                
              
              −
              (
              U
              
                c
                
                  z
                
              
              
                )
                
                  2
                
              
            
          
        
        =
        0
      
    
    {\displaystyle {S_{x}^{2}c_{x}^{2} \over S^{2}-(Uc_{x})^{2}}+{S_{y}^{2}c_{y}^{2} \over S^{2}-(Uc_{y})^{2}}+{S_{z}^{2}c_{z}^{2} \over S^{2}-(Uc_{z})^{2}}=0}
  til bestemmelse av strålehastigheten i retningen gitt ved S. Den har helt samme struktur som Fresnel-ligningen i k-rommet og kan sies å være den resiproke (eller duale) til denne. For en gitt vektor S er den en andregradsligning for energitettheten U og har dermed også i allminnelighet to løsninger for denne og den tilhørende strålehastigheten. Løsningene ligger på en sammensatt flate gitt ved U = konstant som definerer stråleflaten i S-rommet.


=== Geometrisk konstruksjon ===
De to løsningene for strålehastigheten v kan også finnes ved en geometrisk konstruksjon basert på Fresnel-ellipsoiden på samme måte som fasehastigheten u kunne bestemmes fra indeksellipsoden ved å snitte den med et plan gjenom sentrum. Snittplanet  S⋅x = 0 er vinkelrett på strålevektoren S og går gjennom sentrum av Fresnel-ellipsoiden som det skjærer i en ellipse. Halvaksene til denne ellipsen har lengder som gir de to løsningene for strålehastigheten målt i enheter av den vanlige lyshastigheten c. Samtidig vil retningene til halvaksene angi polarisasjonen for den tilsvarende strålen.
Fresnel-ellipsoiden for en enakset krystall er rotasjonssymmetrisk om den optiske aksen med  nx = ny =  no. Den ene halvaksen til snittellipsen vil da ha en konstant lengde. Denne representerer den ordinære strålen som brer seg utover med den konstante strålehastigheten  c /no. Den andre halvaksen har en lengde som varierer med vinkelen χ som strålevektoren S danner med den optiske aksen,


=== Enaksete krystaller ===

Den duale Fresnel-ligningen  forenkles for enaksete krystaller der cx =  cy =  co = c /no og cz = ce =  c /ne. Den kan da faktoriseres og tar formen

  
    
      
        
          (
          
            
              U
              
                2
              
            
            −
            
              
                
                  S
                  
                    x
                  
                  
                    2
                  
                
                
                  c
                  
                    o
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  S
                  
                    y
                  
                  
                    2
                  
                
                
                  c
                  
                    o
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  S
                  
                    z
                  
                  
                    2
                  
                
                
                  c
                  
                    o
                  
                  
                    2
                  
                
              
            
          
          )
        
        
          (
          
            
              U
              
                2
              
            
            −
            
              
                
                  S
                  
                    x
                  
                  
                    2
                  
                
                
                  c
                  
                    e
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  S
                  
                    y
                  
                  
                    2
                  
                
                
                  c
                  
                    e
                  
                  
                    2
                  
                
              
            
            −
            
              
                
                  S
                  
                    z
                  
                  
                    2
                  
                
                
                  c
                  
                    o
                  
                  
                    2
                  
                
              
            
          
          )
        
        =
        0
      
    
    {\displaystyle \left(U^{2}-{\frac {S_{x}^{2}}{c_{o}^{2}}}-{\frac {S_{y}^{2}}{c_{o}^{2}}}-{\frac {S_{z}^{2}}{c_{o}^{2}}}\right)\left(U^{2}-{\frac {S_{x}^{2}}{c_{e}^{2}}}-{\frac {S_{y}^{2}}{c_{e}^{2}}}-{\frac {S_{z}^{2}}{c_{o}^{2}}}\right)=0}
  Andregradsligningen er derfor et produkt av to faktorer som hver for seg beskriver en del av den fulle stråleflaten. Fra den første faktoren ser man at den ordinære strålen sprer seg ut som en kulebølge med den konstante strålehastigheten  co. Den andre faktoren beskriver en ellipsoide med halvakser ce, ce og co i de tre hovedretningene som er stråleflaten for den ekstraordinære bølgen.
Hvis man innfører vinkelen χ som strålevektoren S danner med den optiske aksen, er dens komponent langs denne aksen Sz = S cosχ . Da vil også Sx2 +  Sy2 = S2 sin2χ slik at strålehastigheten v = S/U i retning χ er gitt som

  
    
      
        
          
            1
            
              
                v
                
                  2
                
              
              (
              χ
              )
            
          
        
        =
        
          
            
              
                sin
                
                  2
                
              
              ⁡
              χ
            
            
              c
              
                e
              
              
                2
              
            
          
        
        +
        
          
            
              
                cos
                
                  2
                
              
              ⁡
              χ
            
            
              c
              
                o
              
              
                2
              
            
          
        
      
    
    {\displaystyle {1 \over v^{2}(\chi )}={\sin ^{2}\chi  \over c_{e}^{2}}+{\cos ^{2}\chi  \over c_{o}^{2}}}
  I hvert plan som inneholder den optiske aksen, brer de ekstraordinære strålene seg derfor utover med en hastighet som beskriver en ellipse. Langs den optiske aksen beveger den ordinære og ekstraordinære strålen seg med samme hastighet. Dette resultatet for strålehastigheten har samme form som uttrykket for den brytningsindeksen n(θ) i retning θ  for den ekstraordinære bølgen.


== Strålegang i enaksete krystaller ==
 
Når en lysstråle går gjennom en krystall, vil den vanligvis splittes i to deler. Hver del er polarisert vinkelrett til den andre og beveger seg med ulike hastigheter bestemt av retningene til de optiske aksene.  Denne strålegangen kan enklest beskrives for enaksete krystaller. Det er også de som er vanligvis blir benyttet i praktiske sammenhenger basert på slik dobbeltbrytning.
I det enkleste tilfellet faller lyset vinkelrett på en sideflate til krystallen med den optiske aksen i innfallsplanet. Er den og loddrett på sideflaten, vil både den ordinære og den ekstraordinære bevege seg langs denne aksen. De vil derfor ha samme hastighet og følge samme, rette linje gjennom krystallen.
Hvis derimot den optiske aksen er parallell med sidekanten, vil de to bølgene bevege seg vinkelrett på den optiske aksen. De to resulterende strålene med motsatt polarisasjon forplanter seg derfor  inn i krystallen i samme retning, men med ulike hastigheter c /no og c /ne. Denne viktige effekten blir gjort nytte av i forskjellige polarisatorer.Når man beskriver utbredelse av bølgen ved bruk av Huygens-Fresnels prinsipp, vil hvert punkt på den innkommende bølgefronten gi opphav til nye elementærbølger som bygger opp en ny bølgefront inni krystallen. De ordinære elementærbølgene beveger seg med samme hastighet c /no i alle retninger og utgjør en kulebølge. Men de ekstraordinære elementærbølgene har en strålehastighet v(χ) som varierer med retningen og vil gi en utbredelse med samme form som den tilsvarende stråleflaten, det vil si en ellipse i innfallsplanet. Likevel vil de i dette tilfellet bygge opp en ny bølgefront som blir parallell med sideflaten. Ved Huygens konstruksjon er den tangenten til alle disse ellipseformede strålefrontene. Den tilsvarende strålevektoren S står derfor også vinkelrett på krystallens sideflate og er parallell med bølgevektoren k.


=== Snells brytningslov ===

I begge disse tifellene med lys som faller loddrett inn på en sideflate, er det antatt at bølgeretningen inne i krystallen har samme retning. Dette er en konsekvens av Snells brytningslov når den generaliseres til dette tilfellet hvor den brutte bølgen består av to deler, hver med sin egen brytningsindeks. Den ordinære bølgen som er polarisert vinkelrett på planet som inneholder den optiske aksen, har brytningsindeksen n0 som er en konstant. Men den ekstraordinære bølgen som er polarisert i dette planet, har en brytningsindeks n(θ) som varierer med vinkelen som bølgevektoren danner med den optiske aksen.
For begge bølgene kan den tilsvarende brytningsloven utledes ved bruk av grensebetingelsene for feltene i overgangen fra luft til krystall på samme måte som for vanlig lysbrytning. Setter man brytningsindeksen i luft lik med en, vil da brytningsvinkelen θo2  for den ordinære bølgen være gitt ved

  
    
      
        sin
        ⁡
        
          θ
          
            1
          
        
        =
        
          n
          
            o
          
        
        sin
        ⁡
        
          θ
          
            o
            2
          
        
      
    
    {\displaystyle \sin \theta _{1}=n_{o}\sin \theta _{o2}}
  der θ1 er innfallsvinkelen.  Dette er den vanlige brytningsloven. På samme måte blir brytningsvinkelen θe2  for den ekstraordinære bølgen bestemt ved ligningen

  
    
      
        sin
        ⁡
        
          θ
          
            1
          
        
        =
        n
        (
        θ
        )
        sin
        ⁡
        
          θ
          
            e
            2
          
        
      
    
    {\displaystyle \sin \theta _{1}=n(\theta )\sin \theta _{e2}}
  hvor nå θ vil avhenge av denne brytningsvinkelen. Derfor kan denne ligningen bare løses tilnærmet ved numeriske metoder. Men i det spesielle tilfellet med loddrett innfall, blir begge brytningsvinklene null som allerede antatt.At de to brytningsvinklene er null for loddrett innfall gjelder også når den optiske aksen ligge i samme plan, men danner en vinkel med sideflaten. Denne vinkelen er da også vinkelen θ mellom k-vektor for den ekstraordinære bølgen og den optiske aksen. Den tilsvarende, ekstraordinære strålen vil da danne vinkelen α  med innfallsretningen gitt ved det etablerte resultatet

  
    
      
        tan
        ⁡
        α
        =
        
          
            
              (
              
                n
                
                  e
                
                
                  2
                
              
              −
              
                n
                
                  o
                
                
                  2
                
              
              )
              tan
              ⁡
              θ
            
            
              
                n
                
                  e
                
                
                  2
                
              
              +
              
                n
                
                  o
                
                
                  2
                
              
              
                tan
                
                  2
                
              
              ⁡
              θ
            
          
        
      
    
    {\displaystyle \tan \alpha ={(n_{e}^{2}-n_{o}^{2})\tan \theta  \over n_{e}^{2}+n_{o}^{2}\tan ^{2}\theta }}
  Stråleretningen S ligger derfor nærmere eller lenger fra den optiske aksen enn bølgevektoren k avhengig av om krystallen er positiv eller negativ. I dette spesielle tilfellet med loddrett innfall splittes derfor den innkommende strålen i to ved å trenge inn i krystallen. Den ordinære strålen fortsetter rett frem, mens den ekstraordinære går ut til siden i en retning gitt ved denne vinkelen.


=== Skjevt innfall ===

Mer komplisert er det å beregne retningen til den ekstraordinære strålen når det innfallende lyset danner en vinkel θ1 > 0 med innfallsloddet. Det enkleste tilfellet fremkommer når den optiske aksen står vinkelrett på sideflaten. Bølgevektoren til denne strålen danner da samme vinkel θ2 med innfallsloddet som med den optiske aksen slik at θ = θ2. Ved å benytte uttrykket

  
    
      
        n
        (
        θ
        )
        =
        
          
            
              
                n
                
                  o
                
              
              
                n
                
                  e
                
              
            
            
              
                
                  n
                  
                    o
                  
                  
                    2
                  
                
                
                  sin
                  
                    2
                  
                
                ⁡
                θ
                +
                
                  n
                  
                    e
                  
                  
                    2
                  
                
                
                  cos
                  
                    2
                  
                
                ⁡
                θ
              
            
          
        
        
        ,
      
    
    {\displaystyle n(\theta )={n_{o}n_{e} \over {\sqrt {n_{o}^{2}\sin ^{2}\theta +n_{e}^{2}\cos ^{2}\theta }}}\;,}
  i Snells lov for denne ekstraordinære strålen, har dens bølgevektor inni krystallen en retning θ2 som er implisitt gitt ved ligningen

  
    
      
        sin
        ⁡
        
          θ
          
            1
          
        
        =
        
          
            
              
                n
                
                  o
                
              
              
                n
                
                  e
                
              
              tan
              ⁡
              
                θ
                
                  2
                
              
            
            
              
                
                  n
                  
                    e
                  
                  
                    2
                  
                
                +
                
                  n
                  
                    o
                  
                  
                    2
                  
                
                
                  tan
                  
                    2
                  
                
                ⁡
                
                  θ
                  
                    2
                  
                
              
            
          
        
      
    
    {\displaystyle \sin \theta _{1}={n_{o}n_{e}\tan \theta _{2} \over {\sqrt {n_{e}^{2}+n_{o}^{2}\tan ^{2}\theta _{2}}}}}
  Den kan kvadrereres og dermed lett løses med resultatet

  
    
      
        tan
        ⁡
        
          θ
          
            2
          
        
        =
        
          
            
              
                n
                
                  e
                
              
              sin
              ⁡
              
                θ
                
                  1
                
              
            
            
              
                n
                
                  o
                
              
              
                
                  
                    n
                    
                      e
                    
                    
                      2
                    
                  
                  −
                  
                    sin
                    
                      2
                    
                  
                  ⁡
                  
                    θ
                    
                      1
                    
                  
                
              
            
          
        
      
    
    {\displaystyle \tan \theta _{2}={n_{e}\sin \theta _{1} \over n_{o}{\sqrt {n_{e}^{2}-\sin ^{2}\theta _{1}}}}}
  Når denne vinkelen er kjent, finnes retningen til den ekstraordinære stråleretningen fra den generelle sammenhengen 
  
    
      
        tan
        ⁡
        χ
        =
        (
        
          n
          
            o
          
          
            2
          
        
        
          /
        
        
          n
          
            e
          
          
            2
          
        
        )
        tan
        ⁡
        θ
      
    
    {\displaystyle \tan \chi =(n_{o}^{2}/n_{e}^{2})\tan \theta }
   som gir

  
    
      
        tan
        ⁡
        χ
        =
        
          
            
              
                n
                
                  o
                
              
              sin
              ⁡
              
                θ
                
                  1
                
              
            
            
              
                n
                
                  e
                
              
              
                
                  
                    n
                    
                      e
                    
                    
                      2
                    
                  
                  −
                  
                    sin
                    
                      2
                    
                  
                  ⁡
                  
                    θ
                    
                      1
                    
                  
                
              
            
          
        
      
    
    {\displaystyle \tan \chi ={n_{o}\sin \theta _{1} \over n_{e}{\sqrt {n_{e}^{2}-\sin ^{2}\theta _{1}}}}}
  Samme resultatet kommer man frem til ved å betrakte denne retningen slik som når den bestemmes ved  Huygens konstruksjon.


== Se også ==
Lysbrytning
Dobbeltbrytning


== Referanser ==


== Litteratur ==
E.A. Wood, Crystals and Light: An Introduction to Optical Crystallography, Dover Publications, New York (1977). ISBN 0-486-23431-2.