EOS EOS

Size: px
Start display at page:

Download "EOS EOS"

Transcription

1 EOS EOS

2

3 SN1987A Crab nebula (SN1054) A.K.Mann "Shadow of a Star" From

4 重力崩壊型超新星爆発 中性子星 ブラックホールの誕生 重元素の起源 rプロセス 爆発的元素合成 銀河の進化 宇宙線 ガンマ線バースト 重力波 爆発メカニズムの解明 重力崩壊 コアバウンス 3 10 Explosion Dynamics Trapping, heating 2 10 ν νν ν Pt, Au ν ν SUBAR U Mass number Origin of Elements r-process, ν-process Super-Kamiokande ν ν ν Dense Matter Astronomy Probe inside

5 Fe core ρ c ~10 10 g/cm 3 T c ~1 MeV ν ν- ρ c ~10 12 g/cm 3 T c ~2 MeV ν ν ν ν 1000 km ρ c ~3x10 14 g/cm 3 T c ~5 MeV ν ν ν ρ c ~4x10 14 g/cm 3 T c ~10 MeV ν ν ν Shockwave NS ~ 10 km NS T~0 MeV Y e <0.1

6 M Fe ~M NS ~ 1.4M solar, R Fe ~10 4 km R NS ~10 km ΔE Grav = GM 2 GM 2 ~ erg R Fe R NS E exp ~ erg E ν ~ erg Macrophysics Microphysics

7 多次元流体力学計算 超新星爆発多次元流体計算の例 2D calc with rotation Shen-EOS 超新星爆発残骸 形状, 元素合成, 偏光 中性子星キック 非対称性 球対称:1D (非)軸対称:2D, 3D 回転, 磁場 流体力学的不安定性 対流, 組成混合 Jun Wang 2002 Jan Kotake et al. ApJ 595 (2003) 304

8 / ν ν f(t,x,y,z,p x,p y,p z ) 6+1 ν ν ν ( ) / /

9 / ν ν f(t,x,y,z,p x,p y,p z ) 6+1 ν ν f(e ν f ν (E ν ) ν ( ) 1.0 ρ c =10 14 g/cm 3 & 0.8 bounce ρ c =10 12 / / ρ c = g/cm MeV 100 MeV 1000 E ν [MeV] ν energy

10 9 SciDAC Review Spring 2006

11 ν (ρ, T, Y e ) : ε, p, S, µ n, µ p, X i, m *, ρ ν

12 : σ ~ A 2 ν / : e - + p ν + n e - + A ν + A e + + n ν + p : ν i + N ν i + N ν i + A ν i + A ν i + e ν i + e : e - + e + ν e + ν e γ ν i + ν i N + N N + N + ν i ν i i=e, µ, τ

13 Lattimer-Swesty EOS (1991) Skyrme-Hartree-Fock approach Relativistic EOS (Shen-Oyamatsu-Toki-Sumiyoshi, 1998) Relativistic Mean Field approach Hyperon EOS (Ishizuka-Ohnishi, 2006) Variational Approach (Kanzawa-Takano, 2005~) - 3D Skyrme Hartree-Fock EOS (Newton, Stone Mezzacappa, 2007) A=460 A=1400 (n b =0.08 fm -3, T=2.5 MeV, Y p =0.3) Fuller, Fowler & Newman, Bruenn (1989) Langanke-Pinedo (45<A<112) T. Otsuka & T. Suzuki, Fujita-RCNP Reddy, Horowitz, Yamada-Toki,

14

15 EOS n,p,e ρ ρ0 21pSA T=0 EOS

16 Approaches to obtain the EOS approach starts from ingredients Theory/Model empirical the parametrized EOS nuclear mass, size,... Phenomenological effective NN int. nuclear mass, size,... Liquid-Drop Model Droplet Model Thomas-Fermi Theory... Skyrme HF RMF AMD... many-body bare NN int. (AV18, Bonn, Paris,...) NN scattering,... Variational Calc. DBHF...

17 SN-EOS for supernova simulations Empirical J.M. Lattimer, F.D. Swesty NPA535(1991). Phenomenological (RMF-TM1) H. Shen, H. Toki, K. Oyamatsu, K. Sumiyoshi, NPA637(1998), PTP100(1998). many-body (AV18) under construction (21pSA-11) H. Kanzawa, K. Oyamatsu, K. Sumiyoshi, M. Takano, NPA791(2007) Mukherjee, Pandharipande, Phys. Rev. C 75, (2007).

18 The Nuclear EOS for Supernova Simulations! based on the Realistic Nuclear Potential! H. Kanzawa a, K. Oyamatsu b, K. Sumiyoshi c and M. Takano a! a Waseda Univ., b Aichi Shukutoku Univ., c Numazu CT! 1) The uniform EOS: AV18+UIX,! Vatiational many-body calculations! 2) The Thomas-Fermi calculations for atomic nuclei reproduce their empirical data.! NPA791 (2007) 232! 3) The SN-EOS is under construction!

19 Hix (Proc. Nuclei in the Cosmos X to be published) Protoneutron star size 20% 200ms Shock location 20% lepton, entropy gradient g/cm 3 BH neutrino 21pSA9,10

20 (flux, spectra) : (40M solar ) Sumiyoshi et al., PRL (2006) ν ν Shen EOS ν 50 ν cf. Hyperon WW95: 40M solar ν 40 ν e ν e ν µ < E ν > [MeV] ν e ν e 10 Hyperon- 0.7 Shen- 1.3 ν time after bounce [sec] 1.5 Sumiyoshi, Yamada & Suzuki ApJ (2007) Sumiyoshi et al. (2008) in preparation

21 (Prakash) 0.5 ρ0 ρ0 GDR ρ0 2-3 ρ0 Heavy Ion Collision Bao et al., Phys. Rep.(2008) in press.(arxiv: ) ρ > 2-3 ρ0 Hyperon, Meson Condensate, Qurak Matter (

22 T=0 w n, x w 0 + K 0 18n 0 2 n n x 2 S 0 + L 3n 0 n n 0 ds( n) L = 3n 0 dn n=n0 S 0 = S( n ) 0

23 L, K0 EOS K. Oyamatsu and K. Iida, PTP 109, , 2003.

24 Efforts being done Estimate L value from global behavior of nuclear mass and size in nuclear chart with Prof. Iida (Kouchi U.), Drs. Kohama(RIKEN), Koura(JAEA), Abu-Ibrahim(Cairo U.) For nuclear size, we need to directly compare calculations and cross section measurements. Kurotama (Black sphere) model(21asb7) Glauber calculation of cross sections Preliminary result from nuclear mass L value is relatively large. (closer to EOS C than to EOS G)

25 (Shen EOS) Wigner-Seitz cell free n, p, α Thomas-Fermi

26

27 arxiv: v1 (astro-ph)

28 RMF calculation by Maruyama, Chiba, Tatsumi 12$&2($&*30&3*"$(%#(425"*!!"#$%&'()*+,-"$(%#(./(0"--! Pasta structures appear.! Y p =0.5! Y p =0.1!!!

29 QMD calculations by Sonoda, Watanabe, Sato, Takiwaki, Yasuoka, Ebisuzaki Pasta phases at zero temperature Cooling of Hot nuclear matter(~10 MeV) down to 0.1 MeV Snapshots calculated for model 2 Sphere Rod Slab 0.100! ! ! 0 Red!Proton Blue!Neutron! 0 =0.168 fm -3 "Nuclear density#! Rod-like Bubbles 0.490! 0 Spherical Bubbles 0.575! 0 14

30 Shen EOS table size log10 (T(MeV) ) 31 grids log10(yp) grids log10(ρb (g/cm3)) grids ,904 data points

31 EOS L / GDR K0 Heavy ion collision Nuclear Statistical Equilibrium Hartree-Fock RMF

32 hyperons, meson condensate, quark matter,...

33 EOS T, ρ, Yp Mega Grids nucastrodata.org, bigbangonline.org

1401_HPCI-lecture4.SNEOS.pptx

1401_HPCI-lecture4.SNEOS.pptx 極限物質の性質を決めるには? 超新星 : 状態方程式データテーブル 中性子星と超新星の状態方程式 中性子星 密度のみの関数 ほぼ中性子物質 ゼロ温度 冷えた中性子星 多くの状態方程式 原子核実験 中性子星質量 半径 超新星 密度だけでなく 電子の割合が変わる 有限温度 超新星爆発時 少ないデータテーブル 数値シミュレーション 中性子星合体にも 極限状態での物質の性質 状態方程式 (Equation

More information

PowerPoint Presentation

PowerPoint Presentation 2010 KEK (Japan) (Japan) (Japan) Cheoun, Myun -ki Soongsil (Korea) Ryu,, Chung-Yoe Soongsil (Korea) 1. S.Reddy, M.Prakash and J.M. Lattimer, P.R.D58 #013009 (1998) Magnetar : ~ 10 15 G ~ 10 17 19 G (?)

More information

@ τ weak 1 σ weak n target v relative sec ( T 10MeV ) 2 ( σ weak 4G2 F h2 c 2 π T 2, n target ρ m u, v relative c ( τ dyn 1 0.4msec

@ τ weak 1 σ weak n target v relative sec ( T 10MeV ) 2 ( σ weak 4G2 F h2 c 2 π T 2, n target ρ m u, v relative c ( τ dyn 1 0.4msec @ 2011.12.26 τ weak 1 σ weak n target v relative 8 10 8 sec ( T 10MeV ) 2 ( σ weak 4G2 F h2 c 2 π T 2, n target ρ m u, v relative c ( τ dyn 1 0.4msec Gρ ) 1/2 ρ 10 14 g/cm 3 ) 1 ρ 10 14 g/cm 3 (T 10MeV,

More information

EOS and Collision Dynamics Energy of nuclear matter E(ρ, δ)/a = E(ρ, )/A + E sym (ρ)δ 2 δ = (ρ n ρ p )/ρ 1 6 E(ρ, ) (Symmetric matter ρ n = ρ p ) E sy

EOS and Collision Dynamics Energy of nuclear matter E(ρ, δ)/a = E(ρ, )/A + E sym (ρ)δ 2 δ = (ρ n ρ p )/ρ 1 6 E(ρ, ) (Symmetric matter ρ n = ρ p ) E sy Nuclear collision dynamics and the equation of state We want to measure EOS. Measure T, P and ρ of matter... Prepare matter in the state we want to measure HI collisions What are taking place in collisions?

More information

ニュートリノ駆動型 超新星爆発シミュレーション 3Dと2Dの比較

ニュートリノ駆動型 超新星爆発シミュレーション 3Dと2Dの比較 2016/01/05 宇宙の歴史をひもとく地下素粒子原子核研究 第二回超新星ニュートリノ研究会 ニュートリノ観測による超新星爆発メカニズムの探求 滝脇知也 ( 理化学研究所 ) 1 2 Three phases of supernovae Presupernova phase Burst phase Accretion phase Cooling phase Kato+2015 吉田敬さん石徹白さん

More information

W 1983 W ± Z cm 10 cm 50 MeV TAC - ADC ADC [ (µs)] = [] (2.08 ± 0.36) 10 6 s 3 χ µ + µ 8 = (1.20 ± 0.1) 10 5 (Ge

W 1983 W ± Z cm 10 cm 50 MeV TAC - ADC ADC [ (µs)] = [] (2.08 ± 0.36) 10 6 s 3 χ µ + µ 8 = (1.20 ± 0.1) 10 5 (Ge 22 2 24 W 1983 W ± Z 0 3 10 cm 10 cm 50 MeV TAC - ADC 65000 18 ADC [ (µs)] = 0.0207[] 0.0151 (2.08 ± 0.36) 10 6 s 3 χ 2 2 1 20 µ + µ 8 = (1.20 ± 0.1) 10 5 (GeV) 2 G µ ( hc) 3 1 1 7 1.1.............................

More information

BH BH BH BH Typeset by FoilTEX 2

BH BH BH BH Typeset by FoilTEX 2 GR BH BH 2015.10.10 BH at 2015.09.07 NICT 2015.05.26 Typeset by FoilTEX 1 BH BH BH BH Typeset by FoilTEX 2 1. BH 1.1 1 Typeset by FoilTEX 3 1.2 2 A B A B t = 0 A: m a [kg] B: m b [kg] t = t f star free

More information

Microsoft PowerPoint - kotake

Microsoft PowerPoint - kotake 超新星爆発メカニズムの現状 固武慶 ( 国立天文台 ) 超新星爆発とニュートリノ原子核反応 阪大 RCNP 研究会 2007 年 2 月 目次及び時刻表 第 1 章 Introduction (~5 分 ) 第 2 章超新星の物理 (~5 分 ) (standard supernova scenario) 第 3 章超新星爆発メカニズム最前線 (~30 分 ) Asymmetry と爆発メカニズム

More information

( )

( ) 1. 2. 3. 4. 5. ( ) () http://www-astro.physics.ox.ac.uk/~wjs/apm_grey.gif http://antwrp.gsfc.nasa.gov/apod/ap950917.html ( ) SDSS : d 2 r i dt 2 = Gm jr ij j i rij 3 = Newton 3 0.1% 19 20 20 2 ( ) 3 3

More information

positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) MeV : thermalization m psec 100

positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) MeV : thermalization m psec 100 positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) 0.5 1.5MeV : thermalization 10 100 m psec 100psec nsec E total = 2mc 2 + E e + + E e Ee+ Ee-c mc

More information

1/2 ( ) 1 * 1 2/3 *2 up charm top -1/3 down strange bottom 6 (ν e, ν µ, ν τ ) -1 (e) (µ) (τ) 6 ( 2 ) 6 6 I II III u d ν e e c s ν µ µ t b ν τ τ (2a) (

1/2 ( ) 1 * 1 2/3 *2 up charm top -1/3 down strange bottom 6 (ν e, ν µ, ν τ ) -1 (e) (µ) (τ) 6 ( 2 ) 6 6 I II III u d ν e e c s ν µ µ t b ν τ τ (2a) ( August 26, 2005 1 1 1.1...................................... 1 1.2......................... 4 1.3....................... 5 1.4.............. 7 1.5.................... 8 1.6 GIM..........................

More information

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H 199 1 1 199 1 1. Vx) m e V cos x π x π Vx) = x < π, x > π V i) x = Vx) V 1 x /)) n n d f dξ ξ d f dξ + n f = H n ξ) ii) H n ξ) = 1) n expξ ) dn dξ n exp ξ )) H n ξ)h m ξ) exp ξ )dξ = π n n!δ n,m x = Vx)

More information

42 3 u = (37) MeV/c 2 (3.4) [1] u amu m p m n [1] m H [2] m p = (4) MeV/c 2 = (13) u m n = (4) MeV/c 2 =

42 3 u = (37) MeV/c 2 (3.4) [1] u amu m p m n [1] m H [2] m p = (4) MeV/c 2 = (13) u m n = (4) MeV/c 2 = 3 3.1 3.1.1 kg m s J = kg m 2 s 2 MeV MeV [1] 1MeV=1 6 ev = 1.62 176 462 (63) 1 13 J (3.1) [1] 1MeV/c 2 =1.782 661 731 (7) 1 3 kg (3.2) c =1 MeV (atomic mass unit) 12 C u = 1 12 M(12 C) (3.3) 41 42 3 u

More information

スライド 1

スライド 1 Matsuura Laboratory SiC SiC 13 2004 10 21 22 H-SiC ( C-SiC HOY Matsuura Laboratory n E C E D ( E F E T Matsuura Laboratory Matsuura Laboratory DLTS Osaka Electro-Communication University Unoped n 3C-SiC

More information

1 2 1 a(=,incident particle A(target nucleus) b (projectile B( product nucleus, residual nucleus, ) ; a + A B + b a A B b 1: A(a,b)B A=B,a=b 2 1. ( 10

1 2 1 a(=,incident particle A(target nucleus) b (projectile B( product nucleus, residual nucleus, ) ; a + A B + b a A B b 1: A(a,b)B A=B,a=b 2 1. ( 10 1 2 1 a(=,incident particle A(target nucleus) b (projectile B( product nucleus, residual nucleus, ) ; a + A B + b a A B b 1: A(a,b)B A=B,a=b 2 1. ( 10 14 m) ( 10 10 m) 2., 3 1 =reaction-text20181101b.tex

More information

TQFT_yokota

TQFT_yokota , TY, Naito, Phys. Rev. B 99, 115106 (2019),, 2019 9 2 1 (DFT) (DFT)? HΨ(x 1,, x N ) = EΨ(x 1,, x N ) N DFT! Hohenberg, Kohn, PR (1964) Kohn, Sham, PRA (1965) (EDF) E[ρ] = F[ρ] + dxv(x)ρ(x) δe[ρ] δρ(x)

More information

g µν g µν G µν = 8πG c 4 T µν (1) G µν T µν G c µ ν 0 3 (1) T µν T µν (1) G µν g µν 2 (1) g µν 1 1 描

g µν g µν G µν = 8πG c 4 T µν (1) G µν T µν G c µ ν 0 3 (1) T µν T µν (1) G µν g µν 2 (1) g µν 1 1 描 419 特集 宇宙における新しい流体力学 - ブラックホールと SASI- SASI Study of SASI in Black Hole Accretion Flows by Employing General Relativistic Compressive Hydrodynamics Hiroki NAGAKURA, Yukawa Institute for Theoretical Physics,

More information

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) * * 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *1 2004 1 1 ( ) ( ) 1.1 140 MeV 1.2 ( ) ( ) 1.3 2.6 10 8 s 7.6 10 17 s? Λ 2.5 10 10 s 6 10 24 s 1.4 ( m

More information

Mott散乱によるParity対称性の破れを検証

Mott散乱によるParity対称性の破れを検証 Mott Parity P2 Mott target Mott Parity Parity Γ = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 t P P ),,, ( 3 2 1 0 1 γ γ γ γ γ γ ν ν µ µ = = Γ 1 : : : Γ P P P P x x P ν ν µ µ vector axial vector ν ν µ µ γ γ Γ ν γ

More information

1 a b cc b * 1 Helioseismology * * r/r r/r a 1.3 FTD 9 11 Ω B ϕ α B p FTD 2 b Ω * 1 r, θ, ϕ ϕ * 2 *

1 a b cc b * 1 Helioseismology * * r/r r/r a 1.3 FTD 9 11 Ω B ϕ α B p FTD 2 b Ω * 1 r, θ, ϕ ϕ * 2 * 448 8542 1 e-mail: [email protected] 1. 400 400 1.1 10 1 1 5 1 11 2 3 4 656 2015 10 1 a b cc b 22 5 1.2 * 1 Helioseismology * 2 6 8 * 3 1 0.7 r/r 1.0 2 r/r 0.7 3 4 2a 1.3 FTD 9 11 Ω B ϕ α B

More information

C el = 3 2 Nk B (2.14) c el = 3k B C el = 3 2 Nk B

C el = 3 2 Nk B (2.14) c el = 3k B C el = 3 2 Nk B I [email protected] 217 11 14 4 4.1 2 2.4 C el = 3 2 Nk B (2.14) c el = 3k B 2 3 3.15 C el = 3 2 Nk B 3.15 39 2 1925 (Wolfgang Pauli) (Pauli exclusion principle) T E = p2 2m p T N 4 Pauli Sommerfeld

More information

2 X-ray 6 gamma-ray 7 1 17.1 0:38m 0:77m nm 17.2 Hz Hz 1 E p E E = h = ch= (17.2) p = E=c = h=c = h= (17.3) continuum continuous spectrum line spectru

2 X-ray 6 gamma-ray 7 1 17.1 0:38m 0:77m nm 17.2 Hz Hz 1 E p E E = h = ch= (17.2) p = E=c = h=c = h= (17.3) continuum continuous spectrum line spectru 1 17 object 1 observation 17.1 X electromagnetic wave photon 1 = c (17.1) c =3 10 8 ms ;1 m mm = 10 ;3 m m =10 ;6 m nm = 10 ;9 m 1 Hz 17.1 spectrum radio 2 infrared 3 visual light optical light 4 ultraviolet

More information

(Blackbody Radiation) (Stefan-Boltzmann s Law) (Wien s Displacement Law)

(Blackbody Radiation) (Stefan-Boltzmann s Law) (Wien s Displacement Law) ( ) ( ) 2002.11 1 1 1.1 (Blackbody Radiation).............................. 1 1.2 (Stefan-Boltzmann s Law)................ 1 1.3 (Wien s Displacement Law)....................... 2 1.4 (Kirchhoff s Law)...........................

More information

0406_total.pdf

0406_total.pdf 59 7 7.1 σ-ω σ-ω σ ω σ = σ(r), ω µ = δ µ,0 ω(r) (6-4) (iγ µ µ m U(r) γ 0 V (r))ψ(x) = 0 (7-1) U(r) = g σ σ(r), V (r) = g ω ω(r) σ(r) ω(r) (6-3) ( 2 + m 2 σ)σ(r) = g σ ψψ (7-2) ( 2 + m 2 ω)ω(r) = g ω ψγ

More information

64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () m/s : : a) b) kg/m kg/m k

64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () m/s : : a) b) kg/m kg/m k 63 3 Section 3.1 g 3.1 3.1: : 64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () 3 9.8 m/s 2 3.2 3.2: : a) b) 5 15 4 1 1. 1 3 14. 1 3 kg/m 3 2 3.3 1 3 5.8 1 3 kg/m 3 3 2.65 1 3 kg/m 3 4 6 m 3.1. 65 5

More information

The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Sh

The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Sh The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Shigeru Kitamura, Member Masaaki Sakuma Genya Aoki, Member

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc.com/ 1 19 3 19.1................... 3 19.............................. 4 19.3............................... 6 19.4.............................. 8 19.5.............................

More information

イメージング分光によるMeVガンマ線天文学の展望

イメージング分光によるMeVガンマ線天文学の展望 髙田淳史 ( 京大理 ) 元素合成 SNR : 放射性同位体銀河面 : 26 Al 電子陽電子対消滅線粒子加速ジェット (AGN) : シンクロトロン + 逆コンプトン強い重力場 Black hole : 降着円盤, π 0 Etc. ガンマ線パルサー, 太陽フレア 1-30 MeV MeV sky map CGRO/COMPTEL Bad Sensitivity Good erg / (cm 2

More information

Microsoft PowerPoint - hoshino_part1.pptx

Microsoft PowerPoint - hoshino_part1.pptx 天体プラズマにおける粒子加速機構 無衝突系粒子シミュレーション 星野真弘理学系研究科 地惑惑星科学専攻 内容 ** PIC (Particle In Cell) シミュレーション ** プラズマ輸送係数の研究 (MHD では現象論的に扱う粘性 電気抵抗 熱伝導に関わる物理 ) 非熱的プラズマ ( 局所的熱平衡ではない物理 ) 電子とイオンのエネルギー分配 宇宙での高エネルギー粒子の観測 数値チェレンコフ問題

More information

(e ) (µ ) (τ ) ( (ν e,e ) e- (ν µ,µ ) µ- (ν τ,τ ) τ- ) ( ) ( ) ( ) (SU(2) ) (W +,Z 0,W ) * 1) [ ] [ ] [ ] ν e ν µ ν τ e µ τ, e R,µ R,τ R (2.1a

(e ) (µ ) (τ ) ( (ν e,e ) e- (ν µ,µ ) µ- (ν τ,τ ) τ- ) ( ) ( ) ( ) (SU(2) ) (W +,Z 0,W ) * 1) [ ] [ ] [ ] ν e ν µ ν τ e µ τ, e R,µ R,τ R (2.1a 1 2 2.1 (e ) (µ ) (τ ) ( (ν e,e ) e- (ν µ,µ ) µ- (ν τ,τ ) τ- ) ( ) ( ) ( ) (SU(2) ) (W +,Z 0,W ) * 1) [ ] [ ] [ ] ν e ν µ ν τ e µ τ, e R,µ R,τ R (2.1a) L ( ) ) * 2) W Z 1/2 ( - ) d u + e + ν e 1 1 0 0

More information

MUFFIN3

MUFFIN3 MUFFIN - MUltiFarious FIeld simulator for Non-equilibrium system - ( ) MUFFIN WG3 - - JCII, - ( ) - ( ) - ( ) - (JSR) - - MUFFIN sec -3 msec -6 sec GOURMET SUSHI MUFFIN -9 nsec PASTA -1 psec -15 fsec COGNAC

More information

19 σ = P/A o σ B Maximum tensile strength σ % 0.2% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional

19 σ = P/A o σ B Maximum tensile strength σ % 0.2% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional 19 σ = P/A o σ B Maximum tensile strength σ 0. 0.% 0.% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional limit ε p = 0.% ε e = σ 0. /E plastic strain ε = ε e

More information

news

news ETL NEWS 1999.9 ETL NEWS 1999.11 Establishment of an Evaluation Technique for Laser Pulse Timing Fluctuations Optoelectronics Division Hidemi Tsuchida e-mail:[email protected] A new technique has been

More information

Kaluza-Klein(KK) SO(11) KK 1 2 1

Kaluza-Klein(KK) SO(11) KK 1 2 1 Maskawa Institute, Kyoto Sangyo University Naoki Yamatsu 2016 4 12 ( ) @ Kaluza-Klein(KK) SO(11) KK 1 2 1 1. 2. 3. 4. 2 1. 標準理論 物質場 ( フェルミオン ) スカラー ゲージ場 クォーク ヒッグス u d s b ν c レプトン ν t ν e μ τ e μ τ e h

More information

Nosé Hoover 1.2 ( 1) (a) (b) 1:

Nosé Hoover 1.2 ( 1) (a) (b) 1: 1 [email protected] 1 1.1 Nosé Hoover 1. ( 1) (a) (b) 1: T ( f(p x, p y, p z ) exp p x + p y + p ) z (1) mk B T p x p y p = = z = 1 m m m k BT () k B T = 1.3 0.04 0.03 0.0 0.01 0-5 -4-3 - -1 0

More information

36 th IChO : - 3 ( ) , G O O D L U C K final 1

36 th IChO : - 3 ( ) , G O O D L U C K final 1 36 th ICh - - 5 - - : - 3 ( ) - 169 - -, - - - - - - - G D L U C K final 1 1 1.01 2 e 4.00 3 Li 6.94 4 Be 9.01 5 B 10.81 6 C 12.01 7 N 14.01 8 16.00 9 F 19.00 10 Ne 20.18 11 Na 22.99 12 Mg 24.31 Periodic

More information

数値相対論シミュレーション

数値相対論シミュレーション 数値相対論シミュレーション 関口雄一郎 ( 京都大学基礎物理学研究所 ) 内容 数値相対論の概要 最近の研究成果 - アニメーションを中心にー 数値相対論とは? John von Neumann 非線形問題の解法として 高速電子計算機を用いる時代が必ず到来する 数値相対論のターゲット 強重力場の動的変動を伴う現象 ブラックホール形成の瞬間 大質量星の重力崩壊 コンパクト天体連星の合体, etc 重力波放出現象

More information

スケーリング理論とはなにか? - --尺度を変えて見えること--

スケーリング理論とはなにか?  - --尺度を変えて見えること-- ? URL: http://maildbs.c.u-tokyo.ac.jp/ fukushima mailto:[email protected] DEX-SMI @ 2006 12 17 ( ) What is scaling theory? DEX-SMI 1 / 40 Outline Outline 1 2 3 4 ( ) What is scaling theory?

More information

φ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1)

φ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1) φ 4 Minimal subtraction scheme 2-loop ε 28 University of Tokyo Atsuo Kuniba version 2/Apr/28 Formulas Γ n + ɛ = n n! ɛ + ψn + + Oɛ n =,, 2, ψn + = + 2 + + γ, 2 n ψ = γ =.5772... Euler const, log + ax x

More information

研究室ガイダンス(H29)福山研v2.pdf

研究室ガイダンス(H29)福山研v2.pdf J.M. Kosterlitz and D.J. Thouless, Phys. 5, L124 (1972); ibid. 6, 1181 (1973) David J. Thouless J. Michael Kosterlitz + ρ T s KT D.J. Bishop and J.D. Reppy, PRL 40, 1727 (1978) ( ) = 2k B m2 T KT π 2 T

More information

EGunGPU

EGunGPU Super Computing in Accelerator simulations - Electron Gun simulation using GPGPU - K. Ohmi, KEK-Accel Accelerator Physics seminar 2009.11.19 Super computers in KEK HITACHI SR11000 POWER5 16 24GB 16 134GFlops,

More information

68 A mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1

68 A mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1 67 A Section A.1 0 1 0 1 Balmer 7 9 1 0.1 0.01 1 9 3 10:09 6 A.1: A.1 1 10 9 68 A 10 9 10 9 1 10 9 10 1 mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1 A.1. 69 5 1 10 15 3 40 0 0 ¾ ¾ É f Á ½ j 30 A.3: A.4: 1/10

More information

1 9 v.0.1 c (2016/10/07) Minoru Suzuki T µ 1 (7.108) f(e ) = 1 e β(e µ) 1 E 1 f(e ) (Bose-Einstein distribution function) *1 (8.1) (9.1)

1 9 v.0.1 c (2016/10/07) Minoru Suzuki T µ 1 (7.108) f(e ) = 1 e β(e µ) 1 E 1 f(e ) (Bose-Einstein distribution function) *1 (8.1) (9.1) 1 9 v..1 c (216/1/7) Minoru Suzuki 1 1 9.1 9.1.1 T µ 1 (7.18) f(e ) = 1 e β(e µ) 1 E 1 f(e ) (Bose-Einstein distribution function) *1 (8.1) (9.1) E E µ = E f(e ) E µ (9.1) µ (9.2) µ 1 e β(e µ) 1 f(e )

More information

From Evans Application Notes

From Evans Application Notes 3 From Evans Application Notes http://www.eaglabs.com From Evans Application Notes http://www.eaglabs.com XPS AES ISS SSIMS ATR-IR 1-10keV µ 1 V() r = kx 2 = 2π µν x mm 1 2 µ= m + m 1 2 1 k ν = OSC 2

More information