Analiza geometryczna biegu promieni w płytce

Ponieważ - przy opisie zjawiska dwójłomności występują współczynniki załamania światła - jako wprowadzenie przypomnijmy klasyczne zadanie z optyki geometrycznej z płytką płasko-równoległą. Posłuży nam ono jednocześnie do sprawdzenia poprawności działania pierwszego fragmentu programu (listing 1).

 

Promień światła monochromatycznego pada pod kątem α  = 30° na płasko-równoległą płytkę szklaną o grubości d  = 3 cm. Współczynnik załamania materiału płytki n  = 1,5. Ile wynosi przesunięcie promienia świetlnego w płytce? (Ryc. 1).

Odp. 0,58 cm [4].

Ryc.1

obrazek

Listing 12

'DO RYSUNKU 1
'PŁYTKA PŁASKO-RÓWNOLEGŁA

Cls

'WPROWADZANIE DANYCH

Dim As Double f, f1, BD, BC, DE   

Const pi = 3.1415926535
Const n = 1.5
Const d = 3.0

For f = pi / 2 To 0 Step -.01
  f1 = Asin(Sin(f) / n)
  BD = d * Tan(f)
  BC = d * Tan(f1)
  DE = (BD - BC) * Cos(f)
  Print f, DE
  Sleep
Next
End
Ryc. 1

stałe liczbowe i „materiałowe”

zmienne występujące w programie (φ ↔ f}

n – współczynnik załamania światła w ośrodku
d – grubość płytki

 

 

Pętla 1

prawo załamania dla pierwszej krawędzi łamiącej
funkcje trygonometryczne trójkąta prostokątnego

przesunięcie promienia po wyjściu z płytki

kontrola algorytmu

Program

W przypadku zjawiska dwójłomności mamy – jak wspomniano wyżej - do czynienia z dwoma współczynnikami załamania promieni świetlnych. Jako punkt wyjścia przypomnijmy tutaj zjawisko, które geometrycznie przypomina bieg promieni świetlnych w płytce anizotropowej. Chodzi o zależność współczynnika załamania od długości fali (dyspersja światła). Poniższe zadanie przybliży nas wizualnie do zjawiska dwójłomności i jednocześnie pozwoli przetestować następny fragment tworzonego programu (listing 2).

 

Promień światła białego pada pod kątem α  = 30° na płasko-równoległą płytkę szklaną o współczynnikach załamania dla promieni fioletowych i czerwonych równych odpowiednio nf  = 1,53 i ncz  = 1,50. Grubość płytki wynosi d  = 3 cm. Jaka jest szerokość wiązki światła wychodząca płytki? (Ryc. 2).

Odp. 0,2 mm [4]

2  

W niektórych listingach stosujemy polecenie: line (-300,0)-(300,0):line (0,175)-(0,-175): color [foreground],[,background]], które umożliwia ”ręczne” sterowanie kolorami wykresów i tła. Linie układu osi mogą być pomocne dla uchwycenia niektórych właściwości otrzymywanych figur (symetrie, szerokość ciemnych krzyży, itp.) czytelnik zechce sam dokonać wyboru tej opcji.
W prezentowanych w całym artykule figurach ciemne krzyże są przedstawiane w kolorze białym. Zdecydowały względy plastyczne i wizualne. Użytkownik programów nie musi być tego samego zdania.

Ryc.2

obrazek

Listing 2

'DO RYSUNKU 2
'PŁYTKA PŁASKO-RÓWNOLEGŁA

Cls

'WPROWADZANIE DANYCH

Dim As Double f, f1, f2, BC, BC1, EE1   

Const pi = 3.1415926535
Const n1 = 1.53
Const n2 = 1.5
Const d  = 3.0

For f = pi / 2 To 0 Step -.01
  f1  = Asin(Sin(f) / n1)
  f2  = Asin(Sin(f) / n2)
  BC  = d * Tan(f1)
  BC1 = d * Tan(f2)
  EE1 = (BC1 - BC) * Cos(f)
  Print f,EE1
  Sleep
Next
End
Ryc. 2

stałe „materiałowe”

zmienne występujące w programie (φ ↔ f}

n1, n2 – współczynniki załamania światła w ośrodku
dla wybranych długości fali

d – grubość płytki

 

Pętla 1

prawo załamania dla pierwszej krawędzi łamiącej

funkcje trygonometryczne trójkąta prostokątnego

szerkość wiązki po wyjściu z płytki

kontrola algorytmu

Program

Niech teraz n1, n2 - oznaczają współczynniki załamania promienia nadzwyczajnego i zwyczajnego. Geometria zjawiska jest taka sama jak na Ryc. 2, ale będzie nas interesowało przesunięcie fazowe wywołane różnicą przebytych dróg optycznych.

Ryc. 2 wykorzystajmy:

z Δ CDE  → DE  = CDsinφ, ale CD = BD - BC, zatem:

CE = (BD - BC)sinφ.

Przesunięcie fazowe:

δ =  2π (ACn1 - AC1n2 - CE).
λ

Uzupełniamy poprzedni fragment programu (listing 2) i od razu sprawdzamy poprawność jego działania na przykładzie konkretnej sytuacji fizycznej (listing 3).

Listing 3

'RÓŻNICA DRÓG OPTYCZNYCH
'PRZESUNIĘCIE FAZOWE

Cls

'WPROWADZANIE DANYCH

Dim As Double r
Dim As Double f, f1, f2, BC, BC1      
Dim As Double CE, AC, AC1, delta    

Const pi = 3.1415926535
Const n1 = 1.54426
Const n2 = 1.55337
Const d  = 0.04

For f = pi / 2 To 0 Step -.01
  f1  = Asin(Sin(f) / n1)
  f2  = Asin(Sin(f) / n2)
  BC  = d * Tan(f1)
  BC1 = d * Tan(f2)
  CE  = (BC1 - BC) * Sin(f)
  AC  = d / Cos(f1)
  AC1 = d / Cos(f2)
  r   = AC * n1 - AC1 * n2 - CE

  'delta = 2 * pi / l * r

  Print f, r
  Sleep
Next
End
Stałe liczbowe i „materiałowe”

zmienne występujące w programie (φ ↔ f}

n1, n2 – współczynniki załamania światła w ośrodku
d – grubość płytki

 

 

Pętla 1

prawo załamania dla pierwszej krawędzi łamiącej

funkcje trygonometryczne trójkąta prostokątnego

przesunięcie promieniu po wyjściu z płytki

przesunięcie fazowe

r – różnica dróg optycznych

kontrola algorytmu

Program

 

Wiązka światła pada prostopadle na płytkę kwarcu wyciętą równolegle do osi optycznej. Wyznaczyć różnicę dróg optycznych promienia zwyczajnego i nadzwyczajnego po przejściu przez płytkę, jeżeli grubość płytki wynosi 0.04 mm. Współczynniki załamania dla promienia zwyczajnego i nadzwyczajnego są odpowiednio równe 1,54426 oraz 1,55337.

Odp. Δ = d(n1 - n2) = 0,04(1,55337 - 1,54426) = 0,36 µm [5].

 

Niech na kryształ pada zbieżna wiązka promieni świetlnych φ  ∈ <0,π/2>. Możemy wówczas zaobserwować skomplikowane figury interferencyjne powstające wskutek tego, że każdemu kątowi padania wiązki przechodzącej przez kryształ odpowiada inna różnica faz drgań promienia nadzwyczajnego i zwyczajnego. W tej sytuacji dla każdej wartości kąta φ  mamy określone wartości różnic dróg optycznych. Zbadajmy teraz zależność pomiędzy różnicą faz δ, a kątem padania promieni φ  na płytkę; inaczej będziemy szukać takich wartości δ, dla których spełnione są warunki:

δ  = 2  - wzmocnienie wiązki (jasne pierścienie),
δ  = (2k  + 1)π  - wygaszenie światła (ciemne pierścienie).

Przyjmijmy, że natężenie światła fali spolaryzowanej liniowo padającej na płytkę będzie I0 i niech z  i n  oznaczają prostopadłe do siebie kierunki drgań promienia zwyczajnego i nadzwyczajnego w materiale dwójłomnym. Przez α  oznaczymy kąt pomiędzy płaszczyzną polaryzacji drgania padającego na płytkę, a płaszczyzną polaryzacji promienia zwyczajnego z, a przez β  kąt pomiędzy płaszczyzną polaryzacji analizatora z płaszczyzną polaryzacji drgania z (Ryc. 3a).

Ryc. 3

a)

obrazek

b)

obrazek

Amplitudy drgań wytworzone w płytce wynoszą:

 I1 = √I0cosα,  I2 = √I0sinα (1)

Każdemu z punktów płytki odpowiada inna różnica faz między interferującymi ze sobą promieniami świetlnymi. W miarę oddalania się od prostej AB  różnica faz wzrasta. W punktach jednakowo oddalonych od prostej AB  różnica faz δ ma stałą wartość. Punkty te leżą na okręgach o środku A  i mają promień BC  (ryc. 2).

Niech teraz promienie wychodzące z płytki podają na analizator tak ustawiony, aby jego kierunek drgań tworzył z kierunkiem drgań promienia zwyczajnego kąt β (Ry3. b). Wówczas:

I'1 = I1cosβ,    I'2 = I2sinβ.

Uwzględniając (1) otrzymujemy:

 I'1 = √I0cosα cosβ,  I'2 = √I0sinα sinβ (2)

Promienie wychodzące z płytki padają na analizator, który wycina składowe wzdłuż kierunku swojej osi. Drganie te odbywają się w tym samym kierunku, mają ten sam okres i są ze sobą spójne (koherentne), a więc po wyjściu z analizatora mogą z sobą interferować.

W tym momencie pojawia się problem, jak znaleźć wypadkową dwóch drgań o amplitudach I1 oraz I2 odbywających się w jednym kierunku o różnicy faz δ ?. Postępujemy tak jak w kursie fizyki szkoły średniej, rzutując dodawane drgania na odpowiednie osie (Ryc. 4) lub stosując twierdzenie cosinusów. Niezwykle efektowne jest wykonanie tych obliczeń przy pomocy liczb zespolonych, ale nie zawsze jest to możliwe w szkole średniej.

Ryc.4

obrazek

I  = (I12 + I22cosδ)2 + (I22sinδ)2 (3)

Listing 4

'POLARYZATOR+ANALIZATOR

Cls

'WPROWADZANIE DANYCH

Screen 9
Window(-320, 175)-(319, -174)

Dim As Double f, f1, f2, BC, BC1
Dim As Double CE, AC, AC1, delta
Dim As Double a, b, I, I1, I2    

Const pi = 3.14
Const n1 = 1.498
Const n2 = 1.683
Const d  = 1.0
Const l  = 0.0004

Line (-300,0)-(300,0): Line (0,175)-(0,-175)
Color 3

For f = pi / 2 To 0 Step -.01
  f1    = Asin(Sin(f) / n1)
  f2    = Asin(Sin(f) / n2)
  BC    = d * Tan(f1)
  BC1   = d * Tan(f2)
  CE    = (BC1 - BC) * Sin(f)
  AC    = d / Cos(f1)
  AC1   = d / Cos(f2)
  delta = 2 * pi / l * (AC1 * n1 - AC * n2 - CE)
  If delta Mod(2 * pi) = 0 Then 

'ANALIZATOR

    For a = 0 To 2 * pi Step 0.001
      b  = a + pi / 2
      I  = 100
      I1 = Sqr(I) * Cos(a) * Cos(b)
      I2 = Sqr(I) * Sin(a) * Sin(b)

'DRGANIA WYPADKOWE

      I = (I1 + I2 * Cos(delta))^2 + (I2 * Cos(delta))^2
      If I > 3 Then Pset(200 * BC1 * Sin(a), 200 * BC1 * Cos(a))
    Next

  Endif
Next
Sleep
End
Stałe liczbowe i „materiałowe”

zmienne występujące w programie (φ ↔ f}

n1, n2 – współczynniki załamania światła w ośrodku
d –  grubość płytki,
l – długość fali

 

 

 

 

 

 

 

Pętla 1

prawo załamania dla pierwszej krawędzi łamiącej

funkcje trygonometryczne trójkąta prostokątnego

przesunięcie fazowe po wyjściu promieni z płytki

instrukcja warunkowa do tworzenia obrazów

 

 

Pętla 2

amplituda wiązki wejściowej

składowe wiązki wyjściowej (Ryc. 3 i 4)

amplituda wypadkowa wiązki wyjściowej

instrukcja warunkowa do tworzenia końcowych figur

Program

Przystępujemy do skomentowania przedstawionego wyżej listingu. Rozpatrzymy niejako ”po drodze” kilka szczególnych sytuacji fizycznych, a odpowiadające im modyfikacje programu, mogą posłużyć do sprawdzenia poprawności działania przepisanego programu lub dostosowania algorytmu do innego języka programowania.

W pętli 1 wprowadzamy niezbędne stałe: długość fali świetlnej λ  (l), grubości płytki krystalicznej d, współczynniki załamania dla promienia zwyczajnego n1 (n1) i nadzwyczajnego n2 (n2). Na razie przyjmujemy, że współczynnik n1 ma stałą wartość. W zależności od relacji pomiędzy wspomnianymi współczynnikami rozróżniamy kryształy optycznie ujemne (n1 < n2) oraz optycznie dodatnie (n1 > n2). W omawianym programie ko-lejność wprowadzania tych współczynników nie ma znaczenia z uwagi na symetrię zjawiska.

Dalej, w pętli 1, liczymy kąty załamania φ1 (f1) promienia zwyczajnego i φ2 (f2) promienia nadzwyczajnego; wartości te uzyskujemy z prawa załamania światła.

Przerwijmy na chwilę opisywanie programu i sprawdźmy poprawność działania tej części programu dla następujących sytuacji fizycznych:

Ryc. 5, 6, 7 pokazują figury interferencyjne dla różnych materiałów i różnych parametrów fizycznogeometrycznych. Dane dotyczące rysunków zamieszczamy w Tabeli I, co pozwala w prosty sposób - przez po-równanie - uchwycić związki i zależności pomiędzy stałymi charakteryzującymi substancje dwójłomne [6].

Ryc. 5

a)

obrazek

b)

obrazek

c)

obrazek

Ryc. 6

a)

obrazek

b)

obrazek

c)

obrazek

Ryc. 7

a)

obrazek

b)

obrazek

c)

obrazek

d)

obrazek

Tabela 1

L.p. Grubość
płytki d
[mm]
Współczynnik
załamania
promienia
zwyczajnego n1
Współczynnik
załamania
promienia
nadzwyczajnego n2
Długość
fali
świetlnej λ
[µm]
Barwa
światła
Ryc.5 SZPAT ISLANDZKI
a) 1 1,484 1.653 0.687 Czerwona
b) 1 1,486 1.658 0.527 Zielona
c) 1 1,498 1.683 0.400 Fioletowa
Ryc.6 KWARC
a) 1 1.550 1.541 0.687 Czerwona
b) 1 1.556 1.547 0.527 Zielona
c) 1 1.568 1.558 0,400 Fioletowa
Ryc.7 KWARC
a) 4 1,551 1,542 0,656 Czerwono-
pomarańczowa
b) 2 1,551 1,542 0,656
c) 1 1,551 1,542 0,656
d) 0,5 1,551 1,542 0,656

Jeśli zmniejszymy grubość płytki, taką samą różnicę faz δ  znajdujemy dla promieni padających pod większym kątem na płytkę. Wobec tego, w miarę zmniejszania grubości płytki, rosną grubości pierścieni ciemnych i odstępy miedzy nimi (Ryc. 7a-d). załamania promienia zwyczajnego i nadzwyczajnego.

Dla płytki umieszczonej między równoległymi polaroidami uzyskujemy krzyż na przekątnych skośnych. W tym celu wystarczy położyć b  = a w pętli 2 (Ryc. 8).

Ryc.8

a)

obrazek

b)

obrazek

c)

obrazek

d)

obrazek

W pętli 2 występuje skok warunkowy, który ma związek z natężeniem wiązki świetlnej opuszczającej układ. Jak pokazują powyższe przykłady wartości natężenia zależą od wielu parametrów, ale przecież nie uwzględniają pochłaniania i rozpraszania energii świetlnej w płytce krystalicznej. Zjawisko to dostatecznie złożone dla substancji izotropowych komplikuje się jeszcze bardziej w ośrodkach optycznie aktywnych. Wspomnimy tutaj tylko o tzw. pleochronizmie, polegającym na różnym pochłanianiu światła w zależności od orientacji wektora świetlnego (anizotropia absorpcji), czy dichronizmie, gdy jedna ze składowych światła spolaryzowanego jest pochłaniana znaczne silniej niż druga. A przecież straty energii zależą jeszcze od długości fali świetlnej, gęstości substancji, rozproszenia wiązki padającej, emisji wymuszonej itp.[6]. Zagadnienia te wypełniają stronice wielu monografii i dla-tego w niniejszej pracy nie uwzględniano tych zjawisk. Przyjęto stały współczynnik pochłaniania k  = 0,03 (Ryc. 5, 6, 7, 8) i tylko dla porównania k  = 0,1 (Ryc. 9). Widać, że pochłanianie ma na przykład wpływ na szerokość ciemnego krzyża.

Program daje możliwości przewidywania i modelowana zjawisk. Wydaje się, że z dydaktycznego punktu widzenia cenne byłoby następujące podejście. Jak wyglądałaby figura interferencyjna, gdybyśmy rzucili na nią monochromatyczną wiązkę światła nadfioletowego (Ryc. 9a - λ  = 0,2µm) czy podczerwonego (Ryc. 9b - λ  = 1µm) ? Jak zmieniać się będzie ilość prążków?

Uprzednio, dla przejrzystości tworzonego programu, założyliśmy, że współczynnik załamania promienia nadzwyczajnego ma stałą wartość. Tymczasem wiadomo, że istotną cechą substancji optycznie aktywnych jest zależność tego współczynnika od kata padania. Możemy pójść dalej; zależność tę ująć w dowolną formułę ma-tematyczną, spełniająca oczywiście warunki fizyczne. Na przykład: niech wartość współczynnika załamani promienia nadzwyczajnego będzie następującą funkcją kąta padania promieni świetlnych:

- n2 = 1 + φ,

- n2 = 1 + √φ.

Ryc.9

a)

obrazek

b)

obrazek


   I Liceum Ogólnokształcące   
im. Kazimierza Brodzińskiego
w Tarnowie

©2024 mgr Jerzy Wałaszek

Dokument ten rozpowszechniany jest zgodnie z zasadami licencji
GNU Free Documentation License.

Pytania proszę przesyłać na adres email: i-lo@eduinf.waw.pl

W artykułach serwisu są używane cookies. Jeśli nie chcesz ich otrzymywać,
zablokuj je w swojej przeglądarce.
Informacje dodatkowe