Subversion Repositories svnkaklik

Compare Revisions

Ignore whitespace Rev 944 → Rev 945

/dokumenty/skolni/PRA2/GammaSpektra/Gamma_spektra.tex
37,9 → 37,9
\begin{tabular}{|l|l|}
\hline
\multicolumn{ 2}{|c|}{\Large \bfseries FYZIKÁLNÍ PRAKTIKUM FJFI ČVUT V PRAZE \huge\strut} \\ \hline
\textbf{Datum měření:} {18.3.2011} & \textbf{Jméno:} {Jakub Kákona} \\ \hline
\textbf{Pracovní skupina:} {4} & \textbf{Ročník a kroužek:} {Pa 9:30} \\ \hline
\textbf{Spolupracovníci:} {Jana Navrátilová} & \textbf{Hodnocení:} \\ \hline
\textbf{Datum měření:} {12.3.2011} & \textbf{Jméno:} {Jakub Kákona} \\ \hline
\textbf{Pracovní skupina:} {2} & \textbf{Hodina:} {Po 7:30} \\ \hline
\textbf{Spolupracovníci: Viktor Polák} {} & \textbf{Hodnocení:} \\ \hline
\end{tabular}
\end {center}
\end {table}
54,13 → 54,23
\subsection{Zadání}
\begin{enumerate}
\item Pozorujte osciloskopem impulsy přiložených zářičů na výstupu jednokanálového spektrometru. Pokuste se odhadnout tvar spektra.(Osciloskop ukazuje tvary a amplitudy jednotlivých pulsů. Počet pulsů je dán intenzitou barvy a energie výškou impulsu.)
\item Naměřte spektrum impulsů $ ^{137}\rm Cs$ pomocí manuálního měření. Okno volte o šířce 100mV. Zpracujte hodnoty do grafu.
\item Naměřte spektrum impulsů $ ^{137}\rm Cs$ pomocí mnohokanálového analyzátoru. Dobu měření volte alespoň 400s. Přiložte graf a porovnejte s předchozími dvěma metodami.
\item Zkalibrujte X osu 3 bodovou kalibrací pomocí dvojice zářičů $ ^{137}\rm Cs+^{60}\rm Co$.
\item Změřte spektra všech přiložených zářičů multikanálovým analyzátorem. Dobu měření volte 15min. Grafy přiložte do protokolu.
\item Určete neznámý zářič zjištěním polohy hlavního píku a porovnáním s tabulkou.
\item Změřte radiační pozadí v místnosti(zářiče je třeba dát do trezoru). Okomentujte, zda má šum vliv na tvar vámi změřených spekter.
\item Určete rozlišovací schopnost spektrometru pro energii spektrální čáry $ ^{137}\rm Cs$(viz poznámky).
 
\item Naměřte spektrum impulsů $ ^{137}\rm Cs$ pomocí manuálního měření. Okno volte o šířce 100mV. Spektrum graficky zpracujte.
 
\item Naměřte spektrum impulsů $ ^{137}\rm Cs$ jednokanálovým analyzátorem pomocí automatického měření. Okno volte o šířce 100mV. Spektrum graficky zpracujte.
 
\item Mnoho kanálovým analyzátorem naměřte jednotlivá spektra přiložených zářičů $ ^{137}\rm Cs, ^{60}\rm Co, ^{241}\rm Am, ^{133}\rm Ba$. (Spektrum nabírejte 10 minut.)
 
\item Pomocí naměřených spekter najděte kalibrační křivku spektrometru, závislost rozlišení spektrometru na energii záření.
 
\item Z naměřeného spektra $ ^{137}\rm Cs$ určete hodnotu píku zpětného rozptylu, Comptonovy hrany, energii rentgenového píku a energii součtového píku.
 
\item Mnohokanálovým analyzátorem naměřte spektrum neznámého zářiče. Určete tento zářič, pozorujte a zaznamenejte další jevy v jeho spektru. (Spektrum nabírejte 10minut.)
 
\item Mnohokanálovým analyzátorem naměřte spektrum pozadí v místnosti (zářiče uschovejte v trezoru), Najděte v pozadí přirozené zářiče a toto pozadí odečtěte od všech zaznamenaných spekter ještě před jejich vyhodnocením. (Pozadí nabírejte 10minut)
 
\item Graficky určete závislost koeficientu útlumu olova na energii gama záření. (Použijte všechny zářiče současně, jednotlivá spektra nabírejte 5minut)
 
\end{enumerate}
 
\section{Experimentální uspořádání a metody}
67,6 → 77,7
 
\subsection{Teoretický úvod}
Radioaktivita je charakterizována, jako jev při němž se jádro atomu určitého prvku samovolně přemění na jádro jiného prvku, tento proces je často doprovázen emisí vysokoenergetického záření. Jádra s těmito vlastnostmi se nazývají radionuklidy - radioaktivní zářiče. Rozeznáváme $\alpha, \beta, \gamma $ zářiče. Aktivita je veličina charakterizující počet jader, které se přemění za 1s, jednotkou je 1 becquerel [Bq]. Počet jader ve vzorku se řídí rovnicí
 
\begin{equation} \displaystyle N_{(t)}=N_0 \rm {e}^{-\lambda t}, \end{equation}
 
kde $N_0$ je počet jader na počátku měření a $\lambda$ je střední pravděpodobnost rozpadu jádra. Dále platí rovnost
73,140 → 84,191
 
\begin{equation} \displaystyle \lambda=\frac{\rm {ln}2}{T_{\frac{1}{2}}},\end{equation}
kde $T_{\frac{1}{2}}$ je poločas rozpadu, doba za kterou se rozpadne právě $\frac{1}{2}$ jader.
\\
Průchod záření látkou: \\
Comptonův rozptyl: Jedná se vlastně o pružnou srážku fotonu a elektronu. Foton změní svůj směr díky absorpci části energie volným elektronem. \\
Fotoefekt: Foton všechnu energii předá elektronu, který byl vázán v atomu, a nyní je vyražen a pohybuje se.\\
Tvorba elektron/pozitronový pár: Foton má dostatečnou energii ($ \succ 2m_ec^2$), při průletu silným polem zanikne za vniku elektron-pozitronového páru. Následnou anihilací vzniknou dva fotony o energii 511keV.
 
 
\subsubsection{Průchod záření látkou a jeho detekce}
 
Vzhledem k tomu, že gama záření je elektromagnetické záření, které nelze detekovat přímými metodami. Tak je nutné využít vlastností jeho interakcí s látkou.
 
\begin{description}
\item[Comptonův rozptyl] - Jedná se vlastně o pružnou srážku fotonu a elektronu. Foton změní svůj směr díky absorpci části energie volným elektronem.
 
\item[Fotoefekt] - Foton všechnu energii předá elektronu, který byl vázán v atomu ten je nyní vyražen a pohybuje se s kinetickou energií rovnou energii původního fotonu zmenšenou o vazebnou energii elektronu v atomu.
 
\item[Tvorba elektron/pozitronový párů] - Foton má dostatečnou energii ($ > 2m_ec^2$), při průletu silným polem zanikne za vniku elektron-pozitronového páru. Přitom předá část svojí hybnosti částici, která pole vytvořila. Obvykle tento proces nastává poblíž atomových jader. Následnou anihilací vzniknou dva fotony o energii 511 keV.
\end{description}
 
\subsubsection{Stínění gama záření}
 
Při průletu gama fotonů látkou se nemění jejich energie, ale následkem srážek se postupně zmenšuje proud fotonů. Zeslabení monoenergetického svazku v takovém případě probíhá podle exponenciálního zákona
 
\begin{equation}
I(d) = I_0 e^{\mu d}
\end{equation}
 
kde $I(d)$ je intenzita svazku prošlého materiálem o tlouštce $d$, $I_0$ je počáteční intenzita a $\mu $ se nazývá lineární koeficient útlumu. Tohoto exponenciálního útlumu intenzity se využívá při stínění gama záření.
 
\subsection{Pomůcky}
Scintilační detektor, zdroj vysokého napětí NL2410, jednokanálový analyzátor PHYWE, čítač impulsů NL2301, multikanálový analyzátor PHYWE, osciloskop, osobní počítač, zdroje gama záření, USB link PASCO 2100, program pro datový sběr Data Studio, program MEASURE.
Scintilační detektor, zdroj vysokého napětí NL2410, jednokanálový analyzátor PHYWE, čítač impulsů NL2301, multikanálový analyzátor PHYWE, osciloskop, počítač, zdroje gama záření, USB link PASCO 2100, program pro datový sběr Data Studio, program MEASURE.
 
\section{Výsledky a postup měření}
\subsection{Manuální měření se jednokanálovým analyzátorem}
Nejdříve jsme změřili spektrum zářiče $^{137}\rm Cs$ jednokanálovým analyzátorem. Šířku okna detektoru jsme nastavili na 100mV. V tomto okně jsme pak čítačem měřili četnost impulzů. Naměřené hodnoty jsou vidět v grafu.
 
\subsection{Manuální měření s jednokanálovým analyzátorem}
 
\begin{figure}
\begin{center}
\label{amplituda}
\includegraphics [width=100mm] {Cs137_manualne.png}
\caption{Manuální měření $^{137}\rm Cs$}
\end{center}
\end{figure}
Tento krok jsme na pokyn asistenta přeskočili. Na začátku měření jsme pouze stihli připojit osciloskop k měřící aparatuře a zobrazit na něm průběhy výstupních impulzů ze scintilačního detektoru.
 
\subsection{Automatický mnohakanálový analyzátor}
 
Pomocí mnohakanálového analyzátoru připojeného k počítači jsme změřili gamma spektra zářičů $^{241}\rm Am$, $^{60}\rm Co$ a $^{133}\rm Ba$, které jsme pak graficky zpracovali.
Pomocí mnohakanálového analyzátoru připojeného k počítači jsme změřili gamma spektra zářičů $^{241}\rm Am$, $^{60}\rm Co$ a $^{133}\rm Ba$, které jsme pak graficky zpracovali. Spektrum každého ze zářičů bylo nabíráno $(600 \pm 10)$s.
 
\begin{figure}
\begin{center}
\label{amplituda}
\includegraphics [width=100mm] {Am241.png}
\caption{Změřené spektrum $^{241}\rm Am$}
\includegraphics [width=80mm] {Am241.png}
\includegraphics [width=80mm] {Ba133.png}
\includegraphics [width=80mm] {Co60.png}
\includegraphics [width=80mm] {Cs137.png}
\end{center}
\end{figure}
 
 
V naměřeném spektru cesia jsme pak identifikovali další jevy, jako špičku zpětného rozptylu na 218 keV, součtový pík na 1341,31 keV a Comptonovu hranu na 447,306 keV.
 
\begin{figure}
\begin{center}
\label{amplituda}
\begin{center}
\includegraphics [width=100mm] {Ba133.png}
\includegraphics [width=150mm] {Cs137.png}
\end{center}
\caption{Změřené spektrum $^{133}\rm Ba$}
\end{figure}
 
\subsection{Identifikace neznámého zářiče}
 
Ve spektru neznámého zářiče jsme nalezli dva charakteristické píky první s energií 528,4 keV a druhý s mnohem nižší intenzitou a energií 1275,2 keV.
Neznámý zářič jsme určili jako $^{22}\rm Na$, podle charakteristického píku 1274,537keV.
\begin{figure}
\label{amplituda}
\begin{center}
\includegraphics [width=100mm] {Co60.png}
\includegraphics [width=100mm] {neznamy.png}
\end{center}
\caption{Změřené spektrum $^{60}\rm Co$}
\caption{Změřené spektrum neznámého zářiče identifikovaného, jako $^{22}\rm Na$}
\end{figure}
 
Ve spektru je ještě slabý náznak špičky na 695,6 keV, což ale pravděpodobně je artefakt od nedostatečně stíněného cesia.
 
\subsection{Kalibrace detektoru}
 
Ze známých energií zářičů bylo možné kanálům analyzátoru přiřadit konkrétní hodnotu energie a tím získat kalibrační křivku. $ \sigma = (0,6575 \pm 0,0067) / E - (16,4873 \pm 9,031) $
 
\begin{table}[h]
\centering
\begin{tabular}{ccc}
\hline
Kanál & $E_{tab}$ [keV] & $\sigma$\\ \hline
2040 & 1332,492 & 48,6 \\
1813 & 1173,228 & 46,7 \\
97,7 & 59,5409 & 13,8 \\
584 & 356,0129 & 28,4 \\
1040 & 661,657 & 47,3 \\
\hline
\end{tabular}
\caption{Tabulkové hodnoty energií přiřazené jednotlivým kanálům a jejich směrodatné odchylky.}
\label{tkal}
\end{table}
 
 
\begin{figure}
\label{amplituda}
\begin{center}
\includegraphics [width=100mm] {Cs137.png}
\includegraphics [width=100mm] {kalibrace.png}
\end{center}
\caption{Změřené spektrum $^{137}\rm Cs$}
\caption{Kalibrační křivka scintilačního detektoru}
\end{figure}
 
Dále jsme pro píky zářičů určili jejich rozptyl. A z této hodnoty a jejich energií určili energetickou rozlišovací schopnost detektoru. Která je uvedena v grafu \ref{rozliseni}. Nafitovaná funkce je $ \sigma = 11,6 / E + 0,039$
 
\subsection{Identifikace neznámého zářiče}
 
Ve spektru neznámého zářiče jsme nalezli dva charakteristické píky první s energií 550 keV a druhý s mnohem nižší intenzitou a energií 1270 keV.
Neznámý zářič jsme určili jako Na22, podle charakteristického píku 1274.5keV.
\begin{figure}
\label{amplituda}
\begin{center}
\includegraphics [width=100mm] {neznamy.png}
\includegraphics [width=100mm] {kalibrace_sigma.png}
\end{center}
\caption{Změřené spektrum neznámého zářiče}
\caption{Energetická rozlišovací schopnost}
\label{rozliseni}
\end{figure}
 
\subsection{Kalibrace detektoru}
 
Ze známých energií zářičů $^{60}\rm Co$ a $^{137}\rm Cs$ bylo možné kanálům analyzátoru přiřadit konkrétní hodnotu energie a tím získat kalibrační křivku.
\subsection{Útlum olova}
 
Pro měření jsme použili několik vrstev olověných plátů. A změřili nejdříve spektrum pozadí za vrstvou olova 2mm. Potom jsme přidali zářiče a znovu změřili spektrum.
Dále byla přidána ještě jedna vrstva olova 1,5mm a znovu změřeno pozadí i utlumené spektrum se zářiči. Výsledné grafické zpracování je rozděleno do dvou grafů, kvůli odlišným požadavkům na rozsahy os pro jednotlivé zářiče.
 
\begin{figure}
\label{amplituda}
\label{stineni}
\begin{center}
\includegraphics [width=100mm] {kalibrace.png}
\includegraphics [width=150mm] {stineni.png}
\end{center}
\caption{Kalibrační křivka scintilačního detektoru}
\caption{Útlum vrstvy olova pro Co}
\end{figure}
 
Dále jsme pro píky zářičů $^{60}\rm Co$ a $^{137}\rm Cs$ určili jejich šířku v polovině maxima. A z této hodnoty a jejich energií určili energetickou rozlišovací schopnost detektoru.
 
\begin{figure}
\label{amplituda}
\label{stineni2}
\begin{center}
\includegraphics [width=100mm] {Ekalibrace.png}
\includegraphics [width=150mm] {stineni2.png}
\end{center}
\caption{Energetická rozlišovací schopnost}
\caption{Útlum vrstvy olova pro Cs}
\end{figure}
 
\subsection{Útlum olova}
 
Pro měření jsme použili několik vrstev olověných plátů. A změřili nejdříve spektrum za vrstvou olova 4,7mm. Pak jsme otočili detektor a přidali další vrstvu olova. Ovšem právě tato změna konfigurace měření znemožňuje přesně určit útlum v materiálu, neboť se tím změnila i konfigurace zářičů. Nicméně je z naměřených hodnot patrné, že vyšší energie jsou utlumeny méně, než energie nízké.
\begin{table}[h]
\centering
\begin{tabular}{ccc}
\hline
FWHM & $E_{tab}$ [keV] & $\mu $\\ \hline
114,4 & 1332,492 & 0,244\\
109,9 & 1173,228 & 0,20 \\
111,3 & 661,657 & 0,038 \\
\hline
\end{tabular}
\caption{ Pološířky maxim jednotlivých charakteristických energií a zjištěné koeficienty útlumu $ \mu $ pro tyto energie}
\label{tkal}
\end{table}
 
Z grafů pro útlumy je vidět, že při měření byly získány hodnoty pro tři velikosti $d=0$, $d=2$, $d=3,5$ mm. Ašak pro fit jsou využity pouze dva. Je to z důvodu, že
třetí bod nemá pro měření význam a je hrubou chybou. Důsledkem toho ale je, že fit má dva parametry a je fitován přes dva body, což znemožňuje určit jeho nejistotu.
 
Důvod, proč u poslední nejtlustší vrstvy olova vyšel záporný útlum může být například v konfiguraci scintilačního krystalu před fotonásobičem nebot krystal je pravděpodobně umístěn v detektoru rovnoběžně se stenou detektoru. A my jsme stínění na detektor pokládali také rovnoběžně s touto stěnou, to znamená, že boční stínění bylo minimální. následně vzhledem k tomu, že měření bylo prováděno poměrově vůči "stínému pozadí v radiačnímu místnosti" a dářič byl během této doby pokládán nedaleko v rovině scintilačního krystalu, tak nastala situace, že naměřené radiační pozadí v místnosti je větší, než následné měření stínění ozářeného ze zářičů Co a Cs.
 
\begin{figure}
\label{amplituda}
\label{stineni2}
\begin{center}
\includegraphics [width=100mm] {olovo.png}
\includegraphics [width=150mm] {utlum.png}
\end{center}
\caption{Útlum vrstvy olova}
\caption{Útlum pro jednotlivé vrstvy olova s odečteným pozadím}
\end{figure}
 
 
\section{Diskuse}
\begin{enumerate}
\item Po správném nastavení osciloskopu připojeného na jednokanálový spektrometr bylo vidět průběhy mnoha pulzů překrývajících se přes sebe. Kde ve vyšších amplitudách byla patrná jistá nehomogenita,podle které bylo možno předpokládat existenci charakteristického píku.
 
\item Spektrum $ ^{137}\rm Cs$ jsme naměřili pomocí manuálního měření. Použité okno bylo 100mV. Naměřené spektrum je uvedeno v grafu.
\item Tento bod byl vynechán.
 
\item Mnohokanálovým analyzátorem jsme naměřili podobným postupem i spektra dalších zářičů. $^{241}\rm Am$, $^{60}\rm Co$ a $^{133}\rm Ba$
\item Mnohokanálovým analyzátorem jsme naměřili podobným postupem i spektra zářičů $^{241}\rm Am$, $^{60}\rm Co, ^{137}\rm Cs$ a $^{133}\rm Ba$
 
\item Díky znalosti charakteristických energií $^{137}\rm Cs$ a $^{60}\rm Co$ jsme získali kalibrační křivku detektoru a také jeho rozlišovací schopnost v závislosti na energii záření. Obě tyto charakteristiky jsou vyneseny v grafech.
\item Díky znalosti charakteristických energií jsme získali kalibrační křivku detektoru a také jeho rozlišovací schopnost v závislosti na energii záření. Obě tyto charakteristiky byly vyneseny v grafech.
 
\item K přesnému určení bližších parametrů spektra $ ^{137}\rm Cs$ bohužel nemáme potřebná data, protože o část spektra jsme neplánovaně přišli zřejmě v důsledku chyby v softwaru. Což jsme zjistili až po ukončení měření. Nicméně některé hodnoty jsme přibližně určili z poznámek během měření. Pík zpětného rozptylu byl na 190 keV a Comptonova hrana 440keV.
\item V naměřeném spektru cesia jsme pak identifikovali další jevy, jako špičku zpětného rozptylu na 218 keV, součtový pík na 1341,31 keV a Comptonovu hranu na 447,306 keV. Sekundární rentgenové záření je pravděpodobně utopeno v šumu.
 
\item Neznámý zářič jsme díky charakteristické energii identifikovali jako.
\item Neznámý zářič jsme díky charakteristické energii 1275,2 keV identifikovali jako $^{22}\rm Na$.
 
\item Spektrálním analyzátorem jsem také naměřili přírodní pozadí v místnosti. Avšak naměřené intenzity jsou tak malé, že nemohou nijak výrazně ovlivnit tvar naměřených spekter. Zároveň se také nepodařilo v naměřeném pozadí identifikovat konkrétní zářiče, protože naměřená data neobsahují, žádný identifikovatelný pík.
\item Spektrálním analyzátorem jsem také naměřili přírodní pozadí v místnosti. Ale nepodařilo se v naměřeném pozadí identifikovat konkrétní zářiče, protože naměřená data neobsahují, žádný identifikovatelný pík, kromě špiček od našich zářičů špatně odstíněných během měření.
 
 
\end{enumerate}
 
\section{Závěr}
V měření se podařilo získat spektra zářičů $ ^{137}\rm Cs$ , $^{241}\rm Am$, $^{60}\rm Co$ a $^{133}\rm Ba$ kalibrovat scintilační detektor a zjistit jeho energetické rozlišení. A následně i identifikovat neznámý zářič.
V měření se podařilo získat spektra zářičů $ ^{137}\rm Cs$ , $^{241}\rm Am$, $^{60}\rm Co$ a $^{133}\rm Ba$ kalibrovat scintilační detektor a zjistit jeho energetické rozlišení. A následně i identifikovat neznámý zářič jako $^{22}\rm Na$. Zajímavým výsledkem, je "záporný útlum olova" způsobený pravděpodobně geometrií aparatury a citlivostí detektoru na záření přicházející z boku.
 
 
\begin{thebibliography}{10} %REFERENCE
\bibitem{3} {http://praktika.fjfi.cvut.cz/GammaSpektr/GammaSpektr.pdf}{ -Zadání úlohy}
\bibitem{3} {http://praktikum.fjfi.cvut.cz/mod/resource/view.php?id=196}{ -Zadání úlohy}
\end{thebibliography}
 
\end{document}