Ústav fyziky kondenzovaných látek Přírodovědecká fakulta, Masarykova univerzita, Brno Fyzikální praktikum 2 5. Magnetické pole Úkoly k měření Povinná část • Měření horizontální složky intenzity magnetického pole Země Gaussovým magnetometrem. Varianty povinně volitelné části A. Magnetická odezva feromagnetického materiálu (hysterezní smyčka). B. Stínění magnetického pole ve válcové dutině. Povinná část Teorie Znalost průběhu magnetického pole v okolí Země je důležitá pro mnoho oborů jako je například geografie, geologie a podobně. Vlastnosti magnetického pole Země popisuje intenzita magnetického pole, obvykle značená H. V každém bodě můžeme vektor intenzity rozdělit na horizontální a vertikální složku, v dalším se soustředíme jen na měření horizontální složky Hz. Princip metody měření Gaussovým magnetometrem spočívá v porovnání intenzity zemského magnetického pole s intenzitou permanentního magnetu pomocí magnetické střelky jako detektoru směru lokálního magnetického pole. Magnetické pole v okolí magnetického dipólu s dipólovým momentem m je H(r) = 1 4πr3 3(r · m)r r2 − m , (5.1) kde r je polohový vektor vzhledem k poloze magnetického dipólu. V reálném případě se ovšem rozměry permanentního magnetu vzhledem ke vzdálenosti, ve které měříme, nedají zanedbat. Proto je třeba tento vztah integrovat přes celý magnet s danou objemovou hustotou dipólového momentu. Přibližně lze výpočet magnetického pole provést nahrazením tyčového permanentního magnetu dvěma fiktivními magnetickými monopóly o magnetickém množství +p a −p ve vzdálenosti l od sebe, jak je znázorněno na obrázku 5.1. Intenzita magnetického pole se pak spočte pomocí analogie s elektrostatickým dipólem magnetostatickou obdobou Coulombova zákona. Je však třeba zdůraznit, že tyto magnetické monopóly jsou pouze fiktivní a ve skutečnosti jako takové neexistují, slouží pouze jako pomůcka k usnadnění výpočtu. 5. Magnetické pole 42 r l H2 − Hz h +p P2 h+ ϕ2 −p H l +p r z 1 P1 1 ϕ −p H a) 1. Gaussova poloha b) 2. Gaussova poloha magneticke pole Zeme Hz Obrázek 5.1: Schéma experimentálního uspořádání. Magnetické pole v Gaussových polohách (P1 první Gaussova poloha, P2 druhá) v okolí permanentního magnetu a jeho skládání s magnetickým polem Země. Permanentní magnet je vždy orientován kolmo ke směru magnetického pole Země. První Gaussova poloha označuje případ, kdy měříme pole v ose permanentního magnetu. Magnetická intenzita v bodě P1 je dána vztahem H1 = 1 4πµ0 p (r − l/2)2 − p (r + l/2)2 , (5.2) kde r je vzdálenost od středu magnetu a l jeho redukovaná délka. Po úpravě dostaneme vztah H1 = 1 4πµ0 2M r3(1 − λ2)2 , (5.3) kde λ = l 2r a M = pl je magnetický moment magnetu. Magnetické pole v druhé Gaussově poloze P2, v přímce vedoucí středem magnetu a kolmé k jeho ose, sečteme z polí h+ a h−. h+ = h− = 1 4πµ0 p r2 + l2/4 = 1 4πµ0 p r2(1 + λ2) . (5.4) Poměr intenzity H2 k h+ je dán vztahem H2 h+ = l r √ 1 + λ2 . (5.5) Magnetická intenzita v druhé Gaussově poloze se pak spočte jako H2 = 1 4πµ0 M r3(1 + λ2)3/2 . (5.6) Známe tedy intenzitu magnetického pole v bodech P1 a P2. Umístíme magnet tak, aby jeho osa směřovala kolmo ke směru magnetického pole Země. Výchylka magnetky v první Gaussově poloze z jejího původního směru k magnetickému pólu Země je ϕ1, přičemž platí tan ϕ1 = H1 Hz = 1 4πµ0Hz 2M r3(1 − λ2)2 . (5.7) 5. Magnetické pole 43 Obdobně v místě P2 se střelka vychýlí o úhel ϕ2 tan ϕ2 = H2 Hz = 1 4πµ0Hz M r3(1 + λ2)3/2 . (5.8) Z každého z těchto vztahů lze již určit velikost magnetického pole Země, známe-li redukovanou délku magnetu l a velikost magnetického momentu M. Kombinací obou vztahů však můžeme dospět k vyjádření, kde redukovaná délka magnetu přímo nevystupuje. Umocníme-li vztah (5.7) na třetí mocninu a (5.8) na čtvrtou, dostaneme M 4πµ0Hz 3 = 1 8 r9 tan3 ϕ1(1 − λ2 )6 (5.9) M 4πµ0Hz 4 = r12 tan4 ϕ2(1 + λ2 )6 . (5.10) Vzájemným vynásobením těchto vztahů dostaneme M 4πµ0Hz 7 = 1 8 r21 tan3 ϕ1 tan4 ϕ2(1 − λ4 )6 . (5.11) Měříme-li ve vzdálenosti mnohem větší než je délka magnetu, platí r ≫ l a tedy i λ4 ≪ 1, pak platí M Hz = 4πµ0r3 7 tan ϕ1 2 3 tan4 ϕ2. (5.12) Na odmocninu na pravé straně se můžeme dívat jako na geometrický průměr, který můžeme nahradit aritmetickým a dostaneme tak zjednodušený vztah M Hz = 4πµ0r3 7 3 tan ϕ1 2 + 4 tan ϕ2 . (5.13) Tento výraz se od předchozího vztahu (5.12) liší o veličinu řádu λ4, kterou můžeme zanedbat. Magnetický moment magnetu určíme z periody kmitů magnetu v magnetickém poli Země. Je-li osa magnetu stočena vůči magnetickému poli Země o úhel ϕ, pak na něj působí magnetický moment velikosti MHz sin ϕ ≈ MHzϕ. Pohybová rovnice magnetu je pak dána vztahem J d2ϕ dt2 + MHzϕ + Dϕ = 0, (5.14) kde J je moment setrvačnosti magnetu a D je torzní moment závěsu. Používáme vlákno s velmi malým torzním momentem, který můžeme vzhledem k velikosti magnetického silového momentu zanedbat. Magnet potom harmonicky kmitá s kruhovou frekvencí ω danou vztahem ω2 = MHz J . (5.15) Vyjádříme frekvenci pomocí doby kyvu magnetu τ = T/2, kde T je perioda kmitů, a dostaneme MHz = π2J τ2 . (5.16) Moment setrvačnosti válcového magnetu je dán vztahem J = m 4 R2 + l2 3 , (5.17) 5. Magnetické pole 44 −p +p l ϕ H Obrázek 5.2: Kmity permanentního magnetu v magnetickém pole Země. kde m je hmotnost magnetu, R jeho poloměr a l délka. Pro magnet tvaru hranolu je jeho moment setrvačnosti J = m 12 l2 + b2 , (5.18) kde b je šířka magnetu a na výšce nezáleží. Vztahy (5.13) a (5.16) nám udávají veličiny A = M/Hz a B = MHz. Z těchto veličin určíme velikost horizontální složky intenzity magnetického pole Země jako Hz = B A . (5.19) Magnetický moment permanentního magnetu můžeme obdobně určit jako M = √ AB. (5.20) Úkoly 1. Změřte výchylku střelky v obou Gaussových polohách magnetu pro tři různé vzdálenosti r od středu magnetu. Měření provádějte na obě strany od magnetu a také pro magnet otočený o 180◦. 2. Změřte periodu kmitů magnetu v magnetickém poli Země, rozměry a hmotnost magnetu. 3. Určete velikost horizontální složky magnetické pole Země pomocí vztahů (5.13), (5.16) a (5.19). Varianta A: Magnetická odezva feromagnetického materiálu. Teorie Vztah mezi magnetickou intenzitou H a magnetickou indukcí B je dán vztahem B = µ0 (H + M) , (5.21) kde M je vektor magnetizace, který udává objemovou hustotu magnetického momentu. V případě paramagnetických a diamagnetických materiálů v e slabém magnetickém poli můžeme závislost magnetizace na okolním poli předpokládat v lineárním tvaru M = χµ0H, (5.22) kde χ je magnetická susceptibilita, která je kladná pro paramagnetické a záporná pro diamagnetické materiály. Pro většinu materiálů s výjímkou přechodových kovů a jejich sloučenim je 5. Magnetické pole 45 susceptibilita velmi malá okolo 10−6 až 10−9. Zřejmě též platí B = (1+ χ)µ0H = µrµ0H, kde µr je relativní permeabilita. V obecném případě je susceptibilita tenzorem a vektory magnetizace a intenzity nemusejí mít stejný směr. Pro feromagnetické materiály však není závislost magnetické indukce na intenzitě pole lineární a vykazuje hysterezní závislost, jejíž typický průběh ukazuje obrázek 5.3. H B B H C BS R Obrázek 5.3: Typický průběh magnetické hysterezní smyčky. Základní odlišnost feromagnetických materiálů od ostatních je schopnost vykazovat magnetizaci bez vnějšího magnetického pole. Magnetizace každého materiálu může dosahovat pouze jisté maximální hodnoty, kdy jsou všechny přítomné magnetické momenty orientovány stejným směrem. Takováto magnetizace se nazývá nasycená (saturační) Ms a její velikost je dána přibližně součinem koncentrace atomů a magnetického momentu každého atomu. Po odstranění vnějšího magnetického pole zůstává v materiálu remanentní (zbytková) magnetizace MR. Hysterezní křivku dále popisuje veličina zvaná koercitivní pole (koercitivní síla) HC, která udává velikost vnějšího pole, při kterém je celková magnetická indukce v materiálu nulová. Koercitivní pole udává informaci o velikosti pole potřebného ke změně orientace magnetického pole v materiálu. Materiály dělíme podle velikosti koercitivního pole na magneticky měkké (pro HC menší než přibližně 103 A/m) a magneticky tvrdé (pro HC větší než přibližně 104 A/m). Měření budeme provádět na feromagnetickém jádře s dvěma vinutími (transformátoru) buzeném střídavým elektrickým proudem zapojeném podle schématu na obrázku 5.4. Primární vinutí slouží k buzení magnetického pole a na sekundárním snímáme indukované napětí. Intenzitu mag- 1 2 1 2 Obrázek 5.4: Schéma obvodu pro měření magnetického pole ve feromagnetu. netického pole můžeme spočíst podle Ampérova zákona L H · dl = S j · n dS , (5.23) kde integrace na levé straně probíhá podél uzavřené křivky L, na pravé straně přes plochu S jí ohraničenou a j je proudová hustota tekoucí plochou. V případě toroidu je řešení jednoduché, 5. Magnetické pole 46 r rma min H r Obrázek 5.5: Schéma řezu toroidní cívkou. Kroužky uvnitř a vně naznačují průběh proudových vodičů. schématicky je naznačeno na obrázku 5.5. Integraci provedeme podél kružnice s poloměrem r. Z důvodu symetrie má intenzita H podél kružnice všude stejnou velikost a předchozí rovnice pak přejde do tvaru 2πrH = N1I , H = N1I 2πr , (5.24) kde N1 je počet závitů primárního vinutí a I proud tekoucí každým z nich. Magnetická intenzita je tedy přímo úměrná proudu, který měříme jako napětí U1 na rezistoru R1 připojeném do série s proudovou cívkou. Hodnota magnetické intenzity v toroidu je rovna H(t) = N1 2πrR1 U1(t) . (5.25) Pokud je rozdíl vnitřního a vnějšího poloměru dostatečně malý, můžeme považovat hodnotu magnetické intenzity nezávislou na poloze v toroidu a za poloměr r dosadit jeho průměrnou hodnotu r = rmin+rmax 2 . Při buzení střídavým proudem se mění s časem též magnetická indukce. Časová změna magnetické indukce B indukuje v sekundárním vinutí elektromotorické napětí E2 podle Faradayova zákona E2(t) = − dΦ dt = −N2S dB dt , (5.26) kde Φ je celkový magnetický tok sekundární cívkou. Jestliže průřez jádra toroidu je S a počet závitů sekundárního vinutí N2, pak je magnetický tok roven Φ = N2SB. Indukované napětí je úměrné časové změně magnetické indukce. Abychom mohli měřit přímo napětí úměrné magnetické indukci, je v obvodu zařazen integrační RC člen. Průběh napětí na kondenzátoru o kapacitě C získáme z druhého Kirchhoffova zákona E2 = RI2 + UC, UC = Q C , I2 = dQ dt , (5.27) kde I2 je proud tekoucí obvodem a Q je náboj na kondenzátoru. Po úpravě získáme diferenciální rovnici pro náboj Q dQ dt + Q RC = 1 R E2(t). (5.28) Tato rovnice má řešení ve tvaru Q(t) = − 1 R ∞ 0 E2(t − τ)e− τ RC dτ . (5.29) Průběh napětí na kondenzátoru je potom dán vztahem UC(t) = − 1 RC ∞ 0 E2(t − τ)e− τ RC dτ. (5.30) 5. Magnetické pole 47 Je-li časová konstanta integračního obvodu RC mnohem větší než perioda budicího střídavého proudu, lze exponenciální člen v integrálu položit přibližně roven 1. Potom po dosazení z rovnice (5.26) do vztahu (5.30) dostaneme výraz pro napětí UC UC(t) ≈ 1 RC t −∞ N2S dB dt τ dτ , UC(t) ≈ N2S RC B(t). (5.31) Po převedení dostaneme vztah pro magnetickou indukci B(t) = RC N2S UC(t). (5.32) V zapojení podle schématu na obrázku 5.4 nastavíme osciloskop do tzv. X-Y režimu, kdy zobrazujeme vzájemnou závislost napětí na jednotlivých vstupech. Jelikož podle vztahu (5.25) je napětí na prvním vstupu úměrné intenzitě magnetického pole a napětí na druhém vstupu je podle vztahu (5.32) úměrné indukci magnetického pole, zobrazujeme přímo hysterezní smyčku, tedy závislost indukce na intenzitě magnetického pole. Napětí naměřená na osciloskopu pak již převedeme na indukci a intenzitu magnetického pole ve zvolených bodech hysterezní smyčky pomocí výše zmíněných vztahů (5.25) a (5.32). Magnetizaci můžeme snadno spočíst z magnetické indukce s použitím vztahu (5.21) jako M = B µ0 − H. (5.33) Úkoly 1. Zapojte obvod podle schématu. 2. Z osciloskopu odečtěte napětí odpovídající koercitivnímu poli, remanentní a saturační mag- netizaci. 3. Změřte rozměry jádra transformátoru. 4. Určete velikost koercitivního pole, saturační a remanentní magnetizace pro zadaný materiál podle vztahů (5.25) a (5.32). Varianta B: Stínění magnetického pole v dutém válci. Teorie Magnetická permeabilita materiálů µ vyjadřuje vztah mezi magnetickou indukcí a magnetickou intenzitou B = µrµ0H a lze studovat prostřednictvím stínění magnetického pole. Pro neferomagnetické materiály nabývá relativní permeabilita µr hodnot velmi blízkých jedné takže jejich odezva v magnetickém poli se příliš neliší od vakua. Permeabilita feromagnetik souvisí s hysterezní křivkou; permeabilita je úměrná směrnici tečny k hysterezní křivce. Pro feromagnetika může nabývat velmi vysokých hodnot, avšak silně závisí na velkosti magnetického pole. Magneticky měkká feromagnetika (malá koercitivní pole) mají vysokou permeabilitu, zatímco magneticky tvrdé materiály mají permeabilitu nízkou. Magneticky tvrdé materiály mají permeabilitu v řádu desítek až stovek, zatímco magneticky měkké speciální materiály s vysokou permebilitou mohou dosahovat hodnot až 106. Pro malá magnetická pole a magneticky měkké materiály můžeme předpokládat lineární závislost B = µrµ0H s konstantní permeabilitou. 5. Magnetické pole 48 Stínění magnetického pole ve válcové dutině Umístíme-li dutý válec o poloměru R do homogenního magnetického pole velikosti Bo kolmého na osu válce je pole uvnitř válce rovněž homogenní o nižší velikosti Bi. Kompletní výpočet je poněkud zdlouhavý [2, 3], zde se omezíme na uvedení předpokladů: • Magnetická intenzita a indukce splňují Maxwellovy rovnice bez přítomnosti vnějších proudů divB = 0, rotH = 0. (5.34) • Magnetická indukce a intenzita jsou svázány lineárním materiálovým vztahem B = µrµ0H. (5.35) • Tečná složka magnetické intenzity je spojitá na rozhraní dvou prostředí. • Normálová složka magnetické indukce je spojitá na rozhraní dvou prostředí. • Magnetická indukce ve velké vzdálenosti r od osy válce je B(r ≫ R) = Bo. (5.36) Výsledný průběh magnetického pole je znázorněn v obrázku 5.6. Poměr indukce vně Bo a uvnitř trubice Bi vyjadřuje stínící koeficient S, pro nějž platí vztah [2, 3] S = Bo Bi = (µr + 1)2 − b2 a2 (µr − 1)2 4µr , (5.37) kde a je vnější poloměr a b je vnitřní poloměr dutého válce. Pro vysoké hodnoty magnetické permeability µr ≫ 1 a malou tloušťku stěny trubice d vzhledem k jeho poloměru d ≪ R můžeme použít aproximativní vztah S = Bo Bi ≈ 1 + µrd 2R . (5.38) Výše uvedené vztahy platí pro malá magnetická pole. Obzvláště uvnitř materiálů s vysokou permeabilitou může maximální hodnota magnetické indukce (přibližně rovna Bmax ≈ µrBo) snadno překročit saturační magnetizaci materiálu (pro železo asi 2,2 T) a celkový stínící koeficient pak vyjde efektivně nižší. Tato vlastnost se projeví jako závislost stínícího koeficientu na vnějším poli, který s větším vnějším polem klesá. Hodnotu permeability pro nízká pole získáme z hodnot stínícího koeficientu pro nízká pole, kdy u stínící koeficient nezávisí na intenzitě pole. Homogenní magnetické pole v Helmholtzových cívkách Nejjednodušší možnost vytvoření homogeního pole představují tzv. Helmholtzovy cívky. Jsou to dvě cívky o stejném počtu závitů a poloměru R umístěné na společné ose ve vzdálenosti jejich poloměru R od sebe. Magnetické pole jedné cívky můžeme vypočíst pomocí Biotova–Savartova zákona H = I 4πr3 r × dl, (5.39) kde I je proud protékající vodičem, r vzdálenost délkového elementu dl od místa měření pole. Magnetické pole na ose úzké cívky poloměru R o N závitech ve vzdálenosti z od středu cívky se spočte snadno jako H(z) = NIR2 4π(R2 + z2)3/2 2π 0 dϕ = NIR2 2(R2 + z2)3/2 . (5.40) Velikost magnetického ve středu dutiny Helmholtzových cívek získáme jako součet příspěvku obou cívek (vzdálenost středu dutiny od středu každé cívky je z = R/2) H = 2 NIR2 2(R2 + (R/2)2)3/2 = 4 5 3/2 NI R . (5.41) 5. Magnetické pole 49 Z ro Hall trubka Obrázek 5.6: Vlevo: Průběh siločar magnetického pole v okolí a uvnitř dutého válce z feromagnetického materálu s permeabilitou µr = 10. V blízkém okolí válce je homogenita magnetického pole poněkud narušena. Vpravo: Schéma zapojení Helmholtzových cívek. Vyznačena je poloha stínící trubky a Hallovy sondy pro měření magnetického pole. − 2 22 2 2 2 2 poloha cm) Hm) ivka 1 ivka 2 poloha cm) polohacm) .. 1. 1 sa civek ivka 1 ivka 2 − − −2 2 − − −2 2 Obrázek 5.7: Vlevo: Průběh intenzity magnetického pole na ose Helmholtzových cívek s poloměrem R = 5 cm. Vpravo: Rozložení v rovině osy Helmholtzových cívek. Zobrazeny jsou vrstevnice pro hodnoty 0.90, 0.95, 0.99 a 1.01 hodnoty ve středu dutiny. Ve středové oblasti hvězdicovitého tvaru je odchylka velikosti magnetického pole menší než 1 %. 5. Magnetické pole 50 B I F Obrázek 5.8: Princip Hallova jevu. Měření magnetického pole Hallovou sondou K měření velikosti magnetického pole použijema Hallova jevu. Při pohybu nositelů náboje ve vzorku v magnetickém poli (elektrony či díry v polovodiči) na ně působí Lorentzova síla kolmo ke směru jejich pohybu F L = qvd × B, (5.42) kde q je jejich náboj a vd driftová rychlost jejich pohybu. V ustáleném stavu vzniká elektrické pole EH, které eliminuje vliv Lorentzovy síly F H = qEH = −F L. (5.43) Dosadíme za driftovou rychlost vd = j nq = 1 nq I wd, kde j je proudová hustota, n koncentrace nositelů náboje, d tloušťka a w šířka vzorku. Pak porovnáním těchto vztahů dostaneme vztah pro Hallovo napětí UH = EHw = RH d I B, (5.44) kde RH = 1 nq je Hallova konstanta a d je tloušťka vzorku. Znaménko Hallovy konstanty odpovídá znaménku nositelů náboje, umožňuje nám tedy určit typ vodivosti a měřit koncentraci nositelů náboje. Naopak Hallova sonda známých parametrů může sloužit k měření magnetické indukce. V našem případě použijeme komerční Hallovu sondu s integrovaným proudovým zdrojem a zesilovací elektronikou neznámých paramatrů. Postup měření Provedeme měření magnetického pole bez vložené trubky, změříme Hallovo napětí Uo úměrné vnějšímu magnetickému poli Bo. S vloženou trubicí změříme napětí Ui úměrné magnetické indukci uvnitř trubice Bi. Stínící koeficient S je roven podílu napětí S = Uo Ui . (5.45) Pro eliminaci případného špatného nastavení nuly (nulové napětí nemusí odpovídat stavu bez magnetického pole) provádímě měření pro obě komutace proudu (+ a −) a výsledné napětí pak průměrujeme s ohledem na znaménko U = 1 2 (U+ − U−) . (5.46) Úkoly 1. Zapojte Helmholtzovy cívky do obvodu. 2. Změřte stínící koeficient S a rozměry sady poskytnutých válcových trubek. Vnější pole Bo měřte ve středu dutiny bez zasunuté stínící trubky, hodnotu Bi po umístění stínící trubky. Měření proveďte pro několik hodnot proudu procházející cívkami (doporučené hodnoty 0,5 A, 1,5 A a 2,5 A) a zjistěte zda je stínící koeficient nezávislý na intenzitě vnějšího pole. Měřte pro oba směry komutace proudu. 5. Magnetické pole 51 3. Vypočtěte jejich permeabilitu podle vztahu (5.38), případně (5.37). Užití v praxi: Měření magnetického pole má význačné praktické aplikace. Lokální magnetické pole Země je ovlivněno také geologickými poměry a jeho měření se využívá při geofyzikálním průzkumu např. pohybu litosférických desek. Feromagnetické materiály mají také mnoho praktických fyzikálních a elektrotechnických aplikací, kdy je podstatná znalost jejich hysterezní křivky. Magneticky tvrdé materiály se používají jako permanentní magnety, zatímco magneticky měkké materiály se používají při aplikacích vyžadujících snadnou změnu magnetizace jako jsou elektromagnety nebo transformátory. Magneticky měkké materiály se používají rovněž k odstínění vnějšího magnetického pole. Obzvláště důležité je stínění v elektronových mikroskopech, kde by parazitní vnější magnetické pole ovlivňovalo elektronovou optiku mikroskopu. Hallovy sondy měření magnetického pole jsou velmi rozšířeným typem měření a detekce magnetického pole. Hallova jevu se také užívá pro měření koncentrace nositelů naboje např. v polovodičové technologii, detailním studiem tohoto jevu se zabývá úloha 9 předmětu F6390 Praktikum z pevných látek 2(b). Literatura: [1] R.P. Feynman, R.B. Leighton, M. Sands: Feynmanovy přednášky z fyziky s řešenými příklady 2/3, Fragment (2006). [2] J. Perry, Proc. Phys. Soc. London 13, 227 (1894). [3] J. D. Jackson, Classical electrodynamics, John Wiley & Sons, Inc. (1998), kap. 5.