ブラックホール回転エネルギーの! 電磁場による因果的引抜き 熊本大学理学部 小出眞路 ラックホールの回転エネルギーを電磁場により引き抜くさまざまな 機構が提案されてきた.ブラックホールの地平面では物質,エネ ルギー,情報がブラックホールの外から内側に移動・伝播するだ けで電磁場といえどもエネルギーはブラックホールの外から内に しか輸送されない.この因果律を満たしてブラックホールの回転 エネルギーを引き抜くには『負の電磁エネルギー』をブラックホー ルの地平面を横切り落とし込むしかない.ここでは,この既知の 観点からブランドフォード"ナエク機構,#$%ペンローズ機構,超 放射などの電磁場を使ったブラックホール回転エネルギーの引き 抜き機構を簡単な式を用いて見直す. 第7階ブラックホール磁気圏勉強会 &'()*+*,(水)-熊本大学 話 題! • 電磁場によるブラックホール回転エネ ルギー引き抜き機構の表式! • ブラックホールまわりの電磁波の伝播 の数値シミュレーション! ./01234526!7/85!49!8:/!;943/<5/! #4<1B/0C! .DE<4@F/G! (HHI!! ! =('! >! X-ray, optical, radio emission ! ∼!+! JD004K12D9! DL!DFM0DN O66/0/<12D9! DL!P015K1! !-rays ! =(''! >! ?/3/<10!04@:8!A/1<5! ?/3/<10!(') ('6 04@:8!A/1<5! AGN Quasar (QSO) Gamma-ray burst (GRB) Micro-quasar! Failed supernova (µQSO) ~ light years Relativistic ~ 1AU jet ~ 1 km Gravitational QD<6/!DL!166/0/<12D9!19E!6D004K12D9R! Forming Spinning collapse V!#1@9/26100A" (S!!#1@9/26!T/0E!LD<6/! Black hole (?) E<43/9!W/8 &S!!.1E412D9!P</55F</! +S!!U15!P</55F</ X016Y!$D0/!#1@9/8D5P:/</! (JD<D91) #1@9/26!Q4/0E!\49/5! Z015K1! Z015K1!%45Y! X016Y! $D0/! X016Y!$D0/!#1@9/8D5P:/</! [<@D5P:/</! U.#$%!]@/9/<10!</01234526!#$%S! ! ^E/10!U.#$%!]</5452348A!G/<DS! ./54523/!U.#$% JD31<4198!LD<K!DL!4E/10!U.#$%!/_F12D95! Unit system c = 1 µ0 = 1 • General relativistic equations of conservation laws: P<DP/<!P1<260/!9FKB/<!E/9548A! !" ( nU " ) = 0 (particle number) !" T µ" )"3/0D648A = 0 (energy and momentum) [9/<@A"KDK/98FK!8/95D< #1`N/00!/_F12D95R! %F10!8/95D<!DL!!µ! Q4/0E!58</9@8:!8/95D< )"6F<</98! E/9548A !" *F µ" = 0 Ohm s law with zero resistivity (ideal MHD condition): Fµ!U ! = 0 Space-time around rotating black hole ←Kerr Metric: \1P5/!LF962D9R ! = h02 + ! ( hi" i ) i (gravitational time delay, gravitational potential) 2 ?:4a!3/68D<R !i = hi" i / # ! = ( !1, ! 2 , !3 ) (velocity of dragged frame) Frame-dragging effect and ergosphere !ˆ cbc ˆ O e/<<!! B016Y!:D0/! cbcd!c cbc=c rˆ Q4EF6410!DB5/<3/<! ]Q^%gS!L<1K/!! $D<4GD9! [<@D5P:/</! Q<1K/"E<1@@49@!/f/68! In the ergosphere, any matter, energy, and information can t move or propagate in the opposite direction of the black hole rotation. +h(!QD<K1045K!DL!^E/10!U.#$%![_F12D9 ∼!!54K401<!8D!9D9</01234526!4E/10!#$%! ]6D95/<3123/!LD<KS! (conservation of particle number) Special relativistic mass density, !"# !D = "# $ '(% D ( vˆ + & ) )* !t general relativistic effect Special relativistic total momentum density ! Pˆ (equation of motion) = "# $ '(% Tˆ + & Pˆ )* " ( D + + ) #% + % fcurv " Pˆ : !t special relativistic effect 2 2 ˆ ˆ " % B E 2 ˆ ˆ ˆ ˆ ˆ P = h! v + E " B T = h! vˆvˆ + $ p + + ' I ( Bˆ Bˆ ( Eˆ Eˆ ( ) # Special relativistic total energy density 2 & 2 !" (equation of energy) = #$ % ()& Pˆ # Dvˆ + "' *+ # Pˆ %$& # Tˆ : , !t ˆ $ E ! Bˆ ˆ '! Bˆ + # & Eˆ ) ! J + " # + = % & = "# $ '(% Eˆ " & $ Bˆ )* e ?4K401<48A!N48:!#$%!:15!BDD58/E!8:/!</6/98! ( * !t $t E/3/0DPK/98!DL!4E/10!U.#$%!54KF012D95j! ˆ (Maxwell equations) !e = "# $ E ! " Bˆ = 0 iD!6DFP049@!䈊 ( ) ( ) ) ( ) N48:!D8:/<![_5*! Eˆ + vˆ ! Bˆ = 0 where ( (ideal MHD condition) ! = h02 + # ( hi" i ) i 2 R!]\1P5/!LF962D9SC! !i = hi" i R!]5:4a!3/68D<S! # ! : shear of " Motivation: Extraction of rotational energy from black hole GRMHD simulation of relativistic jet formation eDK4551<D3!&'', #6e499/A!l!U1KK4/!&'')! #1@9/26!5F<L16/5! tk'! tk&'''c/(GM) L158m0DN!E/9548A! DFM0DN 50DNm:4@:8! E/9548A!DFM0DN qD<F5 U.#$%!54KF012D9!DL!F94LD<K0A! K1@9/2G/E!P015K1!1<DF9E!<D8129@! B016Y!:D0/! ?6:N1<G56:40E!<1E4F5! t = 30! S q:/!L158C!0DN! E/9548A!W/8!45! P<DEF6/E!BA! 8:/!B016Y!:D0/! <D812D9* !!!eDK4551<D3!]&'',S!R!n#$%!Z/9<D5/!P<D6/55!45!8<1954/98*!^9!T910C! 58/1EA!518/C!PF</0A!/0/68<DK1@9/26!K/6:1945K!]X019ELD<E" o91W/Y!K/6:1945KS!/`8<1685!<D812D9!/9/<@A!DL!8:/!B016Y!:D0/*p! ジェットのエネルギー源 高速・低密度アウトフロー もともとは,ブラックホールに落ち込んで! 来た物質の重力エネルギー.! 低速・高密度 アウトフロー • ブラックホールの地平面に落ち込む前にその! 自由エネルギーをジェットのエネルギーに変換.! r降着円盤にあるとき! • ブラックホールの地平面に落ち込み,! ブラックホールの回転エネルギーから! ジェットのエネルギーに変換.! #6e499/A!]&''sS数値計算: ジェットの2重構造! ! 外側 !V5FB</01234526!r降着円盤! ! 内側(軸付近)V</01234526rブラックホール・エルゴ領域 X019ELD<E"o91W/Y!K/6:1945K! X019ELD<EC!o91W/Y!! ](HttS!!#i.O?! ・!QD<6/"L<//!]E!fluid << E!EMS!! ・!?8/1EA"5818/! ・ O`4"5AKK/8D<A! ・ XDF9E1<A!6D9E42D9!18!:D<4GD9R! Bˆ$ h$ = (! F " ! H ) P ˆ # B ・ X019ELD<E"o91W/Y!5D0F2D9R! a << 1 ]ブラックホールの回転パラメータS! 2 2 r r ˆ B rˆ = B0 g , B! = " g ( $ F " $ H ) B0 sin 2 % ! #R ・![9/<@A!DF8PF8R! #1@9/26!T/0E! 049/! K1v/<! /9/<@A! r 49LD<K12D9! 4960FE49@! 8/954D9! LD<6/j! ! $D<4GD9 e/<<!B016Y! J1F51048A H :D0/! Rgeu ]ZF950A!l! !JD<D942C!(HH'S S ! ]o91W/Y!(HttS! !F [3/98! :D<4GD9! !F S r # !F (!H " !F ) 18!:D<4GD9! 0 < "F < "H ! S r > 0 /9/<@A!<1E412D9!L<DK! :D<4GD9!8:<DF@:!T/0ER! Z<DB0/K!DL!61F51048A! ]ZF950A!l!JD<D942!(HH' S! q<195PD<8!DL!/0/68<DK1@9/26!/9/<@A!19E! 19@F01<!KDK/98FK! ;948! 5A58/K c =1 µ0 = 1 • General relativistic equations of conservation laws: ! µT µ" = 0 (energy and momentum) [9/<@A"KDK/98FK!8/95D< equations: Q4/0E!58</9@8:!8/95D< • Maxwell ! µ * F µ" = 0 EF10!8/95D<!DL!Qµ$ ! µ F µ" = #J " ! • Force-free condition: Fµ! J = 0 </5452348A Ohm s law with resistivity: Fµ!U ! = " J µ + (U! J ! )U µ ( Kerr Metric: \1P5/!LF962D9R ! = h02 + ! ( hi" i ) i ]@<134812D910!2K/!E/01AS 2 )"6F<</98! E/9548A ! " # : 316FFK ! " 0 : 4E/10!#$% ( x0 , x1, x2 , x3 ) = (t,r,! ," ) ) ?:4a!3/68D<R !i = hi" i / # ! = ( !1, ! 2 , !3 ) ]3/0D648A!DL!E<1@@/E!L<1K/S Several coordinates around rotating black hole ・!\1BD<18D<A!L<1K/C!?6:N1<G56:40E"$40B/<8!5A58/KC!XDA/<" \49E_F458!6DD<E4918/5!!⇒!6DD<E4918/5!DL!@0DB10!L<1K/! µ ds = gµ! dx dx 2 ! !H Kerr BH c!" ZAMO ・!iD<K10!L<1K/C!gB5/<3/<"18"49T948Am?A96:<D9DF5!6DD<E4918/!5A58/K! frame ]g^?S!]5P12100A!DB04_F/!49!@/9/<10S! ⇒!3/68D<!19E!8/95D<! 2 i j ! ! ! ds = !dt + " gij dx dx 2 r!数値計算に有効 i ・ 「空間規格直交iD<K10!L<1K/」C!Q4`/E!0D610!\D</98G!6DD<E4918/!5A58/KC! Q4EF6410!DB5/<3/<!]Q^%gS!L<1K/C!o/<D"19@F01<"KDK/98FK!DB5/<3/<! ]oO#gS!L<1K/!]uS!]5P12100A!D<8:D@D910S! ⇒!3/68D<!19E!8/95D<! ( ) ds = !dtˆ + " dxˆ 2 2 i i 2 (Similar to that of Minkowski metric) ・!JD"KD349@!L<1K/!!⇒56101<!31<41B0/5! r!直感的理解に有効! ]物理量が直感的にS q<195PD<8!DL!/9/<@A!19E!19@F01<!KDK/98FK !" = (0, 0, 0,1) e40049@!3/68D<!LD<!e/<<!5P16/"2K/!wx:! " = (1, 0, 0, 0) 6D95/<312D9!DL!/9/<@A!19E!19@F01<!KDK/98FK: ( ) e! = "#$% T 0% = # eˆ + & ' lˆ' = e!EM + e!hyd l ! = "#$% T 0% = l !EM + l !hyd S i = !"#$ T i$ = S iEM + S ihyd M i = !"#$ T i$ = M iEM + M ihyd !e" = #$ % S !t (/9/<@A!6D95/<312D9) !l " = #$ % M !t $/</C e!EM ( (/0/68<DK1@9/26! ! # Bˆ 2 Eˆ 2 & ˆ ˆ * * ="% + ( + ) +l EM /9/<@A!E/9548A) l EM = h3 E " B 2 ' $ 2 S EM = ! 2 ( Eˆ " # $ Bˆ ) $ ( Bˆ + # $ Eˆ ) (/0/68<DK1@9/26!/9/<@A!zF`! E/9548A) M EM ) ! (19@F01<! KDK/98FK! 6D95/<312D9!) (/0/68<DK1@9/26! 19@F01<!KDK/98FK! E/9548A) ," Bˆ 2 Eˆ 2 % i 3 / ˆ ) Eˆ E ˆ+ E ˆ *B ˆ +1 = ! h3 .$ + ' ( ) Bˆ 3B 3 3 2 & .-# 2 10 2 y!K/915!_F1922/5! DL!Q4EF6410!DB5/<3/<! L<1K/ ( ) (/0/68<DK1@9/26! 19@F01<!KDK/98FK! zF`!E/9548A) ^ B ^ BP {:/9!N/!619!N<48/Eˆ = ! vˆ F" # Bˆ ˆ vˆ / / B), ˆ ( vˆ F = vˆ F! + vˆ F// vˆ F! ! B, F// ^ vF // ^ vF ^ vF eφ ^ ⊥ ^ vFφ ⊥ ^9!@/9/<10,N/!:13/ /0/68<DK1@9/26!/9/<@A"18"49T948AR ! eEM ( ) # Bˆ 2 Eˆ 2 & 2 0 2 -1 ˆ ="% + ( + ) * lEM = " B / 1+ ( vF+ ) + , * vˆF+ 2 2 ' $ 2 .2 1 /0/68<DK1@9/26!/9/<@A!zF`!E/9548AR ( ) ( SEM = ! 2 Eˆ " # $ Bˆ $ Bˆ + # $ Eˆ ) 2 ! % = ! eEM (vˆF& + # ) + 2 (1" vF&2 ) '( Bˆ 2 (vˆF& + # ) " 2 # i Bˆ Bˆ )* ( ) Bφ ZDA929@!zF`!DL!/0/68<DK1@9/26!N13/ {/!6D954E/<![0/68<DK1@9/26!N13/!49!316FFK ( ˆ Bˆ ! exp ( ik " r # i$ t ) = exp i { k ( r,% ) r + l ( r,% )% + m& } # i$ t E, r {/!155FK/!5:D<8!N13/0/9@8:!1PP<D`4K12D9!" k ! 1 $ hi ^9!1!Q^%g!]0D610S!L<1K/C!OKP/</|5!01N!A4/0E5C! # !g jk % !x i '& ) :10L"K4<<D<u floulessence? wave with negative electromagnetic energy injection wave ! d eEM (i) ! !i"ˆ Eˆ = ikˆ # Bˆ ! eEM < 0 8:/9C!N/!:13/! nˆ = kˆ / !ˆ ! ! eEM (o) %45P/<54D9!</012D9!A4/0E5!! nˆ = kˆ / !ˆ = 1 superradiance " $/</C!N/!T9E!!! ˆ S = ! e n + # wave EM EM ! 2 2 ! ! eEM = " (1+ nˆ # $ ) Bˆ = " + % & h& n& Bˆ ! 2 ˆ ( r ! rH ) ' (H RH n& B ! black hole ! ˆ = 1 k = m , !"ˆ = kˆ = 1 k + " $ k = ! 1 " + " $ m k ! ! 0 ! h! h! # 0 # $ # # ergosphere ˆ ! k! # m # m Q49100AC!N/!DB8149nˆ! = = % ! "ˆ h! " $ m" ! RH " $ m&H ( r ! rH ) ! {:/9!}""} d'C!# d'!19E!9/@123/!/9/<@A!DL!/0/68<DK1@9/26!N13/!45!</104G/E!18!8:/!:D<4GD9j ( ) ( ( ) c c ( ) )( ) [9/<@A!zF`!49!8:/!615/!DL!QD<6/"L<// (X019ELD<E"o91W/Y!K/6:1945K) #1`N/00![_F12D95!19E!LD<6/"L<//! 6D9E42D9!A4/0E !v F = h! (# $ % )e , ˆE = ! vˆ # B, ˆ F ! ! F" vˆF! = " vˆF vˆF = " 2 1+ vˆF2 1+ ( Bˆ / Bˆ P ) O8r ! rH (:D<4GD9)C! O8 " # $H r ! rH P EM S = !e " EM ! " # $H ( vˆ F# BP B ^ vF // ^ vF eφ Bˆ! h# 19E!N/!619!N<48/ " (%F & %H ) P ˆ $ B vˆF ! " vˆF! " 1 ^ ^ + $) ^ vF ⊥ ^ ^ Bφ vFφ ⊥ vˆF! Bˆ P = ˆ vˆF B ]o91W/Y!6D9E42D9S i/1<!8:/!:D<4D9,3Q⊥!10N1A5! E4</685!8DN1<E!8:/!:D<4GD9* K1@9/26! ΩF > ΩH ΩF < ΩH 磁力線 磁力線 T/0E!049/5 ^ vF ^ vF ⊥ ^ vF ^ horizon eφ horizon vF ⊥ eφ J1F510!/`8<162D9!DL!B016Y!:D0/!/9/<@A!49!X019ELD<E" o91W/Y!K/6:1945K:/0/68<DK1@9/26!9/@123/!/9/<@A i/1<!8:/!:D<4GD9: ^ BP ^ B ^ ^ 2 -F &1 v ! 2 % % ) ˆ2 F⊥ // v F ˆ eEM = " ( (1+ vˆF# ) + $ vˆF# + B , " B -F . -H '2 * ^φ ^ R B ^ Bˆ = ( Bˆ ! )2 + ( Bˆ P )2 "| Bˆ ! |= ($F % $H ) Bˆ P eφ vF vFφ⊥ # ZD0D4E10!6DKPD9/98!DL!K1@9/26!T/0E ( ) ! EM e ( ) 2 2 1 R & ) H = " Bˆ ( 1+ ( vF# ) + $ % vˆ F# + , -F (-F . -H )( Bˆ PH )2 '2 * " 2 " vˆF! Bˆ P = ˆ 2ˆ v # Bˆ P & F! vˆFB " vˆF! = % ˆ ( vˆF = 1+ vˆF2 $ B' P " SEM = ! eEM ( vˆ F# + $ ) % RH2&F (&F ' &H )( BˆPH )2 ( vˆ F# + $ ) 0 < !F < !H ^9!8:45!615/C!/0/68<DK1@9/26!/9/<@A!45!<1E418/E*!q:45!45!N/00" Y9DN!6D9E42D9!N:/</!X019ELD<E"o91W/Y!K/6:1945K!ND<Y5. ! QF<8:/<KD</C!N/!LDF9E eEM < 0 4*/*!/9/<@A!45!9/@123/.q:45!K/195!/3/9!49!8:/!X019ELD<E"o91W/Y! K/6:1945KC!]/0/68<DK1@9/26S!9/@123/!/9/<@A!45!F204G/E!LD<!8:/!/9/<@A! /`8<162D9!DL!B016Y!:D0/* $DN!1BDF8!^E/10!#$%!615/u! ^E/10!U.#$%!/_F12D95R ;948! 5A58/K c =1 µ0 = 1 • General relativistic equations of conservation laws: proper particle number density ! µ ( nU " ) = 0 ! µT • Maxwell (particle number) )"3/0D648A µ" =0 (energy and momentum) [9/<@A"KDK/98FK!8/95D< equations: Q4/0E!58</9@8:!8/95D< )"6F<</98! E/9548A ! µ * F µ" = 0 Fµ!U ! = 0 Ideal MHD condition: Kerr Metric: \1P5/!LF962D9R ! = h02 + ! ( hi" i ) i ]@<134812D910!2K/!E/01AS 2 ?:4a!3/68D<R !i = hi" i / # ! = ( !1, ! 2 , !3 ) ]3/0D648A!DL!E<1@@/E!L<1K/S JD9581985!10D9@!8:/!K1@9/26!5F<L16/!49!8:/! 5812D91<AC!1`45"5AKK/8<46!#$%!615/ (t, s,!," ) 6DD<E4918/5R ds 2 = !ht2 dt 2 + hs2 ds 2 + h"2 d" 2 + h#2 d# 2 ! 2h#2$ # dtd# . M! (!) = h" h! #$Uˆ , (P1<260/!9FKB/<! 6D95/<312D9) B s (!) = h! h" Bˆ s (E43!Xk') s # !F (") = h$ ' $ $ ) vˆ + % & ( Bˆ $ s * vˆ , s ˆ B + (Q1<1E1A!0DN+4E/10!#$%!6D9E42D9) s B #1@9/26!5F<L16/ sH black hole v ! ˆ B =0 plasma flow vˆ ! = 0 horizon ! ! = =0 & h ˆ " % ˆ" ) L(!) = h" ( U $ B + !t !" (19@F01<!KDK/98FK!6D95/<312D9) ! # m ' * h # ˆ & 2 # 2 Bˆ & & & ˆ& H (!) = # ($ˆ + % U ) ' (B ) + (vˆ + % & ) " "$ˆ m! M! "Uˆ s m! = s = ! ˆ s (/9/<@A!6D95/<312D9) B B JD9581985!10D9@!8:/!K1@9/26!5F<L16/!49!8:/! 5812D91<AC!1`45"5AKK/8<46!#$%!615/ ds 2 = !ht2 dt 2 + hs2 ds 2 + h"2 d" 2 + h#2 d# 2 ! 2h#2$ # dtd# . M! (!) = h" h! #$Uˆ , (P1<260/!9FKB/<! 6D95/<312D9) B s (!) = h! h" Bˆ s (E43!Xk') s # !F (") = h$ ' $ $ ) vˆ + % & ( Bˆ $ s * vˆ , s ˆ B + (Q1<1E1A!0DN+4E/10!#$%!6D9E42D9) O8, r ! rH s B #1@9/26!5F<L16/ sH black hole ˆ ! h! B , " ($F % $H ) s ˆ # B ˆ = ! vˆ # B, ˆ QF<8:/<KD</C!N/!:13/ E F" h! ! $/</, v F = (#F $ %ˆ ! )e! , " v ! ˆ B =0 plasma flow vˆ ! = 0 horizon ! ! = =0 !t !" M! "Uˆ s m! = s = ! ˆ s B B O8,r ! rH ! EM e vF! " 1 R ! RH Bˆ s ! Bˆ Hs " ! # $H RH2 " $F ($F % $H )( Bˆ Hs )2 , # S EM = !F (!F " !H )RH2 ( Bˆ Hs )2 (v! F# + $ ) ! <0 {:/9 0 < !F < !H C!8:/9 eEM ^E/92610!N48:!LD<6/"L<//! 615/(X019ELD<E"o91W/Y! K/6:1945K)! [0/68<DK1@9/26!/9/<@A!45!<1E418/E. $DN/3/<C!49!8:/!#$%!615/C!N/!:13/!8D! 105D!6D954E/<!:AE<DEA91K46!6D98<4BF2D9! !:AE!LD<!8:/!9/8!/9/<@A!<1E412D9* S = S EM + S hyd . " " ). S hyd = ! (ehyd + ! p)(v! + # ehyd # ( 2 ,% / = ! .' " + p * +ˆ $ p 1 + 2 3 l 3hyd # $1 ) -& 0 K1@9/26! ΩF > ΩH ΩF < ΩH 磁力線 磁力線 T/0E!049/5 ^ ^ vF vF ^ vF ⊥ ^ horizon eφ horizon vF ⊥ eφ ?FKK1<A!DL!K/6:1945K5!DL!/9/<@A! /`8<162D9!L<DK!B016Y!:D0/ Table 1. Classification of various mechanisms of energy extraction from black hole. mechanism form of negative energy torque for redistribution of angular output energy references momentum Penrose process magnetic Penrose process Blandford-Znajek mechanism mechanical energy of electrically charged particles electromagnetic energy electromagnetic energy mechanism MHD Penrose process energy extraction with magnetic reconnection superradiance mechanical energy of plasma mechanical energy of plasmoid force of particle fission mechanical energy of particle Penrose (1969) electromagnetic force mechanical energy Wagh (1989) electromagnetic energy Blandford & Znajek (1977) electromagnetic tension force electromagnetic energy and Takahashi et al. (1990); Koi (MHD) kinetic energy (Alfven wave) (2003); Komissarov (2005) Lorentz force (magnetic electromagnetic energy and tension, MHD) kinetic energy (Alfven wave) electromagnetic tension force (force-free) magnetic tension due to magnetic reconnection electromagnetic energy of electromagnetic wave mechanical energy of plasmoid electromagnetic energy of ? electromagnetic wave – 26 – MHD Blandford-Znajek mechanical energy of particle Takahashi et al. (1990); Hirotani et al. (1992); Koid al. (2002); Koide (2003) Koide (2009) Press & Teukolsky (1972); Teukolsky & Press (1974); Lightman et al. (1975) 回転するブラックホール周辺の! 電磁場の伝播の数値計算 放射 プラズマ アウトフロー 磁力線 回転する ブラック ホール エルゴ 領域 超臨界降着 円盤プラズマ • 放射エネルギーの輸送の数値計算,! 目的: • 5FP/<"<1E4196/の数値計算! 方法: • 電磁波を直接計算する。! ただし,波長やスペクトルは気にしない。 回転するブラックホール周辺の電磁場の方程式 { 電磁場テンソル(電場と磁場を表す4元テンソル) 一般相対論的 Maxwell方程式 f ij;k + f jk ;i + f ki ; j = 0 ij i j i = &$ c' 4元電流密度, f ; j = µ0 j % 電磁波の数値計算 ( ここで座標系は x 0 x1 ) x3 = (t x2 x j j x y y z ) (c = 1) xy平面(z=0)上での電磁場の方程式 z ~ B ブラックホール $ &E 1 & 2 ~3 ~1 ~ ~ = % ( J + ( ( B + (' E % (' E ! 1 2 2 1 & t & y ! E !! &E~ 2 & 2 ~3 ~2 ~ ~ = %(J % ( ( B + ('1 E2 % ('2 E1 # x & t & x 電磁波 ! ~ ! &B 3 & ~ & ~ 1 ~3 2 ~3 = % ( E + ' B + ( E % ' B ! 2 1 &x &y !" &t i ( ":時間の遅れ,!:空間のひきずり速度 ) [ ] [ [ j #! z" ] ] [ ] y 数値計算結果 古海亜衣子 (2014) 回転するブラックホール 電磁波を発生させる交流電流 Bz Bz エルゴ領域 エルゴ領域 y y x x a =1 a = !1 (最速で左回りに回転するブラックホール) (最速で右回りに回転するブラックホール) 数値計算結果 a=0 a =1 t = 5.0rg エルゴ領域 「空間の引きずり効果」により 電磁波が引きずり込まれている 数値計算結果 a=0 a =1 a = !1 t = 20.0rg エルゴ領域 地平面 エルゴ領域 回転によって電磁場が引きずられ、電磁波の伝わる速 度が異なっている。 ?FKK1<A " {/!601<4T/E!8:18!/0/68<DK1@9/26!/`8<162D9!DL!B016Y!:D0/! /9/<@A!45!</104G/E!BA!9/@123/!/0/68<DK1@9/26!/9/<@A*!^8! K1Y/5!61F510!/9/<@A!/`8<162D9!L<DK!B016Y!:D0/*!O8!8:/! :D<4GD9C!N/!LDF9E!8:/!</012D9!B/8N//9!8:/! /0/68<DK1@9/26!/9/<@A!E/9548A!19E!/9/<@A!zF`!E/9548A!45! ! P " SEM = ! eEM vˆ F# + $ ! ˆ N:/</!! Eˆ = ! vˆ F" # B. " {/!618/@D<4G/E!8:/!/9/<@A!/`8<162D9!K/6:1945K5!DL! B016Y!:D0/!<D812D910!/9/<@A!N48:!9/@123/!/9/<@A!DN9/<! 19E!8D<_F/!LD<!19@F01<!KDK/98FK!</E458<4BF2D9*! " {/!8<4/E!8D!P/<LD<K!9FK/<4610!54KF012D9!DL!<1E412D9! 8<195L/<!N48:!E4</68!6106F012D9!DL!/0/68<DK1@9/26!N13/5* ( )
© Copyright 2024 ExpyDoc