Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres

Size: px
Start display at page:

Download "Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres"

Transcription

1 Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres Jiří Kubát Astronomický ústav AV ČR Ondřejov Zářivě (magneto)hydrodynamický seminář Ondřejov p.

2 Outline 1. Basic equations for stellar atmospheres 2. Radiative transfer equation in spherical symmetry 3. Formal solution 4. Equations of statistical equilibrium 5. Λ-iteration and accelerated Λ-iteration 6. Model atmosphere problem 7. Linearization 8. Applications 9. Conclusions Zářivě (magneto)hydrodynamický seminář Ondřejov p.

3 Role of radiation in stellar atmospheres source of information about star and stellar atmosphere influence on matter in stellar atmosphere non-local (long distance) interaction (photon mean free path particle mean free path) change of the population numbers (non-equilibrium values) Zářivě (magneto)hydrodynamický seminář Ondřejov p.

4 Basic equations assume static atmosphere ( v = 0) stationary atmosphere ( / t = 0) 1-dimensional spherically symmetric atmosphere equations to be solved (numerically): radiative transfer equation equations of statistical equilibrium equation of radiative equilibrium equation of hydrostatic equilibrium radial depth equation Zářivě (magneto)hydrodynamický seminář Ondřejov p.

5 Basic equations assume static atmosphere ( v = 0) stationary atmosphere ( / t = 0) 1-dimensional spherically symmetric atmosphere equations to be solved (numerically): radiative transfer equation equations of statistical equilibrium equation of radiative equilibrium equation of hydrostatic equilibrium radial depth equation account for full angle and frequency dependence (I (µ, ) I µ ) Zářivě (magneto)hydrodynamický seminář Ondřejov p.

6 Radiative transfer equation µ I µ(r) r + 1 µ2 r I µ (r) µ = χ (r)i µ (r) + η (r) η (z) emissivity χ (z) opacity θ µ = cos θ r Zářivě (magneto)hydrodynamický seminář Ondřejov p.

7 Radiative transfer equation µ I µ(r) r + 1 µ2 r I µ (r) µ = χ (r)i µ (r) + η (r) = χ (r) [I µ (r) S (r)] η (z) emissivity χ (z) opacity S (z) = η (z) χ (z) source function Zářivě (magneto)hydrodynamický seminář Ondřejov p.

8 Radiative transfer equation with scattering µ I µ(z) r + 1 µ2 r I µ (r) µ = [κ (r) + σ (r)]i µ (r) + η (r) + σ (r)j (r) Zářivě (magneto)hydrodynamický seminář Ondřejov p.

9 Radiative transfer equation with scattering µ I µ(z) r + 1 µ2 r I µ (r) µ = [κ (r) + σ (r)]i µ (r) + η (r) + σ (r) 0 I µ (r)dµ Zářivě (magneto)hydrodynamický seminář Ondřejov p.

10 Moment radiative transfer equations ( ) 1 d r 2 H r 2 dr = χ J + η + σ J. where J = 1 2 we introduce f = K J 1 dq q dr 1 3 f = r dk dr + 3K J = χ H, r 1 1 I µ dµ, H = µi µ dµ, K = 1 2 variable Eddington factor the sphericity function 1 1 µ2 I µ dµ. d (f q J ) dx = r 2 H Zářivě (magneto)hydrodynamický seminář Ondřejov p.

11 Moment radiative transfer equations d 2 (f q J ) dx 2 where dx = q χ /r 2 dr boundary conditions = r4 q ( J η ) + σ J, χ where g = d (f q J ) = r 2 ( g J H ) dx d (f q J ) = r 2 ( H + ) g J dx ( ) 1 0 µj µ dµ / (J ) at the surface at depth Zářivě (magneto)hydrodynamický seminář Ondřejov p.

12 Formal solution of the RTE solution for given χ and η (given S ) relatively simple crucial for the total accuracy of the whole problem solution Zářivě (magneto)hydrodynamický seminář Ondřejov p.

13 Formal solution of the RTE ray method Hummer & Rybicki (1971, MNRAS 152, 1) µ = cos ϑ µ r + 1 µ2 r µ = d ds ± di± µ(r) ds = χ (r)i ± µ(r) + η (r) + σ (r)j (r) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

14 Feautrier solution of the RTE solution along a ray (specific intensities I + and I ) introduce Feautrier variables j µ = 1 2 h µ = 1 2 ( I + µ + I µ) ( I + µ I µ) the transfer equation d 2 j µ dτ 2 s = j µ S where S = η + σ J κ + σ, dτ s = (κ + σ ) ds (µ absorbed to ds) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

15 Feautrier solution of the RTE boundary conditions upper lower dj µ dτ s = j µ I µ dj µ dτ s = I + µ j µ dj µ dτ s = 0 core rays (at depth) tangent rays (in the middle of the ray). Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

16 Formal solution of the RTE with scattering 1. determine j µ by solution along each ray (J given): d 2 j µ dτ 2 s = j µ η + σ J κ + σ 2. f = 1 0 µ2 j µ dµ 1 0 j µ dµ, g = 1 0 µj µ dµ 1 0 j µ dµ, 1 dq q dr = 3 f 1 r 3. determine J by solution: 4. check if δf < ε d 2 (f q J ) dx 2 = r4 q ( J η ) + σ J, χ Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

17 Opacities and emissivities opacity χ = i emissivity l>i [ n i g ] i n l α il () + g l i k ( n i n i e h kt n e n k α kk (, T) ) α ik ()+ ( 1 e h kt ) + n e σ e η = 2h3 c 2 [ i l>i n l g i g l α il () + i n i α ik ()e h kt + k n e n k α kk (, T)e h kt ] Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

18 Opacities and emissivities opacity χ = i emissivity l>i [ n i g ] i n l α il () + g l i k ( n i n i e h kt n e n k α kk (, T) ) α ik ()+ ( 1 e h kt ) + n e σ e η = 2h3 c 2 [ i l>i n l g i g l α il () + i n i α ik ()e h kt + n e n k α kk (, T)e h kt ] k Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

19 Equations of statistical equilibrium in LTE: n i = f (n e, T) Saha-Boltzmann distribution Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

20 Equations of statistical equilibrium outside LTE: n i = f (n e, T, J ) statistical equilibrium for i = 1,...NL {n l [R li + C li ] n i (R il + C il )} = 0 l i Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

21 Equations of statistical equilibrium outside LTE: n i = f (n e, T, J ) statistical equilibrium for i = 1,...NL collisional rates {n l [R li + C li ] n i (R il + C il )} = 0 l i n i C il = n i n e q il (T), Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

22 Equations of statistical equilibrium outside LTE: n i = f (n e, T, J ) statistical equilibrium for i = 1,...NL {n l [R li + C li ] n i (R il + C il )} = 0 l i radiative rates αil () n i R il = n i 4π h J d ( g i αij () 2h 3 n l R li = n l 4π g j h c 2 + J ) d Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

23 Equations of statistical equilibrium outside LTE: n i = f (n e, T, J ) statistical equilibrium for i = 1,...NL {n l [R li (J ) + C li ] n i (R il (J ) + C il )} = 0 radiative rates l i αil () n i R il = n i 4π h J d ( g i αij () 2h 3 n l R li = n l 4π g j h c 2 + J ) d Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

24 Final radiative transfer equation d 2 (f q J ) dx 2 f = K J = r4 q K = 1 dq q dr = 3 f 1 r χ = f χ (, n e, T, n i ) Z 1 [ J η ] + σ J, χ 1 «µ 2 I µ dµ η = f η (, n e, T, n i ) X {n l [R li (J ) + C li ] n i (R il (J ) + C il )} = 0 l i Z n i R il = n i 4π Z g i n l R li = n l 4π g j αil () h αij () h J d 2h 3 c 2 + J «d Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

25 Λ iteration radiative transfer equation for J d 2 (f q J ) dx 2 = r4 q [ J η ] + σ J, χ formally may be written as J = Λ S iteration scheme: J ESE S RTE J ESE S l i J (n+1) { [ ] n (n) l R li (J (n) ) + C li = Λ S (n) n (n) i [ R il (J (n) ) + C il ]} = 0 converges extremely slowly for stellar atmospheres (τ 1) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

26 Accelerated lambda iteration simple Λ-iteration convergence problems huge matrices in complete linearization laborious and computationally expensive way out > Accelerated lambda iteration (ALI) formal solution (Cannon 1973, JQSRT 14, 627; ApJ 185, 621) J = Λ S Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

27 Accelerated lambda iteration simple Λ-iteration convergence problems huge matrices in complete linearization laborious and computationally expensive way out > Accelerated lambda iteration (ALI) formal solution (Cannon 1973, JQSRT 14, 627; ApJ 185, 621) J = Λ S + Λ S Λ S Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

28 Accelerated lambda iteration simple Λ-iteration convergence problems huge matrices in complete linearization laborious and computationally expensive way out > Accelerated lambda iteration (ALI) formal solution (Cannon 1973, JQSRT 14, 627; ApJ 185, 621) J = Λ S + (Λ Λ )S Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

29 Accelerated lambda iteration simple Λ-iteration convergence problems huge matrices in complete linearization laborious and computationally expensive way out > Accelerated lambda iteration (ALI) formal solution iteration scheme (Cannon 1973, JQSRT 14, 627; ApJ 185, 621) J = Λ S + (Λ Λ )S J (n+1) = Λ S (n+1) + (Λ Λ )S (n) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

30 Accelerated lambda iteration simple Λ-iteration convergence problems huge matrices in complete linearization laborious and computationally expensive way out > Accelerated lambda iteration (ALI) formal solution (Cannon 1973, JQSRT 14, 627; ApJ 185, 621) J = Λ S + (Λ Λ )S iteration scheme J (n+1) = Λ S (n+1) + (Λ Λ )S (n) }{{} J (n) correction term Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

31 Accelerated lambda iteration simple Λ-iteration convergence problems huge matrices in complete linearization laborious and computationally expensive way out > Accelerated lambda iteration (ALI) formal solution (Cannon 1973, JQSRT 14, 627; ApJ 185, 621) J = Λ S + (Λ Λ )S iteration scheme J (n+1) = Λ S (n+1) + (Λ Λ )S (n) }{{} J (n) correction term compare: J (n+1) = Λ S (n) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

32 Accelerated lambda iteration iteration scheme J (n+1) = Λ S (n+1) + (Λ Λ ) S (n) = Λ S (n+1) + J (n) J (n) correction term S RTE J, J ESE+ALI S solution of the radiative transfer equation is transferred to the solution of the statistical equilibrium equations Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

33 Accelerated lambda iteration iteration scheme J (n+1) = Λ S (n+1) + (Λ Λ ) S (n) = Λ S (n+1) + J (n) J (n+1) radiative rates αil () n i R il = n i 4π h g i αij () n l R li = n l 4π g j h [ Λ S (n+1) ( 2h 3 c 2 + ] + J (n) d [ Λ S (n+1) ] ) + J (n) d Zářivě (magneto)hydrodynamický seminář Ondřejov p. 1

34 Accelerated lambda iteration equations of statistical equilibrium {n l [R li (n i, n l ) + C li ] n i (R il (n i, n l ) + C il )} = 0 l i for i = 1,...NL ALI eliminated radiation field from the explicit solution equations of statistical equilibrium are nonlinear in n i, n l savings important, for a typical problem NF Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

35 Accelerated lambda iteration equations of statistical equilibrium {n l [R li (n i, n l ) + C li ] n i (R il (n i, n l ) + C il )} = 0 l i for i = 1,...NL ALI eliminated radiation field from the explicit solution equations of statistical equilibrium are nonlinear in n i, n l savings important, for a typical problem NF solution by Newton-Raphson method (linearization) preconditioning Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

36 Model atmosphere problem additional structural equations hydrostatic equilibrium radiative equilibrium radial depth preconditioning does not help for them, linearization necessary Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

37 Radial optical depth equation modified column mass depth ( ) R 2 dm = ρ r 2 dr m independent variable in model atmosphere calculation Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

38 Hydrostatic equilibrium equation upper boundary condition: dp g dm = GM R 2 4π c 0 1 d [q f J ] q dm d p 1 = GM m 1 R 2 4π c ( r1 ) 2 R 0 χ 1 ρ 1 [ g1 J 1 H ] d g 1 = 1 0 µj 1µ dµ/ 1 0 j 1µ dµ H is the incident flux at the upper boundary Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

39 Hydrostatic equilibrium equation dp g dm = GM R 2 4π c upper boundary condition: 0 1 d [q f (Λ S + J )] q dm d p 1 = GM m 1 R 2 4π c ( r1 ) 2 R 0 χ 1 ρ 1 [ g1 (Λ S 1 + J 1 ) H ] d g 1 = 1 0 µj 1µ dµ/ 1 0 j 1µ dµ H is the incident flux at the upper boundary J was eliminated using ALO: J = Λ S + J Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

40 Equation of radiative equilibrium integral form 0 [κ J η ] d = 0 κ [J S ] d = 0 differential form (from F = 0) L (4πR) 2 = 0 ρ d [q f J ] q χ dm d Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

41 Equation of radiative equilibrium integral form 0 [κ (Λ S + J ) η ] d = 0 κ [(Λ 1)S + J ] d = 0 differential form (from F = 0) L (4πR) 2 = 0 ρ d [q f (Λ S + J )] q χ dm d J was eliminated using ALO: J = Λ S + J Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

42 Complete linearization to stellar atmospheres introduced by Auer & Mihalas (1969, ApJ 158, 641) for the case of NLTE model atmospheres solution of equations: radiative transfer (J ) hydrostatic equilibrium (ρ) radiative equilibrium (T ) statistical equilibrium (n i ) radial depth (r) vector of solution ψ = (n e, T, r, n 1,...,n NL, J 1,...,J NF ) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

43 Complete linearization to stellar atmospheres introduced by Auer & Mihalas (1969, ApJ 158, 641) for the case of NLTE model atmospheres solution of equations: radiative transfer (J ) included with a help of ALO hydrostatic equilibrium (ρ) radiative equilibrium (T ) statistical equilibrium (n i ) radial depth (r) vector of solution ψ = (n e, T, r, n 1,...,n NL, J 1,...,J NF ) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

44 Summary of equations equations of statistical equilibrium X {n l [R li (n i, n l ) + C li ] n i (R il (n i, n l ) + C il )} = 0 i = 1,... NL l i n i R il = n i 4π R α il h [Λ g S + J ]d, n l R li = n i l 4π R α ij g j h 2h 3 c 2 + [Λ S + J ] d equation of hydrostatic equilibrium dp g dm = GM R 2 4π c equation of radiative equilibrium Z 0 1 d[q f (Λ S + J )] q dm d Z 0 κ [(Λ 1) S + J ] d = 0. or H 0 = Z 0 ρ d[q f (Λ S + J )] q χ dm d. radial depth equation R 2 «dm = ρ r 2 dr Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

45 Solution by linearization vector of solution ψ = (n e, T, r, n 1,...,n NL ), dimension 3 + NL formally F ( ) ψ = 0 current estimate ψ 0 correct solution ψ = ψ 0 + δψ corrections δ ψ = [ F ψ ( ψ0 ) ] 1 F ( ψ0 ) matrix NL + 3 for each depth point d Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

46 Accelerated linearization if we express the linearized b-factors as δb i = b i n e δn e + b i T δt + b i r δr from equations of statistical equilibrium (cf. Auer 1973, ApJ 180, 469; Anderson 1987, NRT, 163), [ b i x = A ij + A im b j b m B i b j ] 1 ( Ajm x b m B j x ) (x stands for n e, T, r) further reduction of the number of explicitly linearized variables Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

47 Accelerated linearization solution of equations: hydrostatic equilibrium (ρ) radiative equilibrium (T ) radial depth (r) statistical equilibrium (n i ) vector of solution ψ = (n e, T, r, n 1,...,n NL ) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

48 Accelerated linearization solution of equations: hydrostatic equilibrium (ρ) radiative equilibrium (T ) radial depth (r) statistical equilibrium (n i ) included with a help of implicit linearization of b-factors vector of solution ψ = (n e, T, r, n 1,...,n NL ) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 2

49 Solution by linearization vector of solution ψ = (n e, T, r), dimension: 3 significant savings of computing time Zářivě (magneto)hydrodynamický seminář Ondřejov p. 3

50 Application of static spherically symmetric model atmospheres generalization of plane-parallel atmospheres in extended atmospheres large radiaton pressure O and B stars strong stellar wind static approximation unusable Zářivě (magneto)hydrodynamický seminář Ondřejov p. 3

51 Application of static spherically symmetric model atmospheres generalization of plane-parallel atmospheres in extended atmospheres large radiaton pressure O and B stars strong stellar wind static approximation unusable Do static extended atmospheres exist? Zářivě (magneto)hydrodynamický seminář Ondřejov p. 3

52 Application of static spherically symmetric model atmospheres generalization of plane-parallel atmospheres in extended atmospheres large radiaton pressure O and B stars strong stellar wind static approximation unusable Do static extended atmospheres exist? they existed! Zářivě (magneto)hydrodynamický seminář Ondřejov p. 3

53 Application of static spherically symmetric model atmospheres generalization of plane-parallel atmospheres in extended atmospheres large radiaton pressure O and B stars strong stellar wind static approximation unusable Do static extended atmospheres exist? they existed! first stars in the Universe consisted only of primordial hydrogen and helium Krtička & Kubát (2006, A&A 446, 1039) showed that such stars can not have stellar winds Zářivě (magneto)hydrodynamický seminář Ondřejov p. 3

54 Application to moving atmospheres static radiative transfer equation may be used for continuum radiative transfer in moving stellar atmospheres, as it was done by Krtička & Kubát (2004, A&A 417, 1003) static spherically symmetric model atmospheres can be used as a lower boundary condition in calculations of stellar winds (core-halo approximation) Zářivě (magneto)hydrodynamický seminář Ondřejov p. 3

55 Final remarks 1. static spherically symmetric model atmospheres are the first step towards more complicated structures 2. for stellar wind calculations they provide the lower boundary condition 3. for axially symmetric problems they give the zeroth order approximation 4. for the first stars in the Universe they provide full description Zářivě (magneto)hydrodynamický seminář Ondřejov p. 3

Solution of the radiative transfer equation in NLTE stellar atmospheres

Solution of the radiative transfer equation in NLTE stellar atmospheres Solution of the radiative transfer equation in NLTE stellar atmospheres Jiří Kubát kubat@sunstel.asu.cas.cz Astronomický ústav AV ČR Ondřejov Non-LTE Line Formation for Trace Elements in Stellar Atmospheres,

More information

PHAS3135 The Physics of Stars

PHAS3135 The Physics of Stars PHAS3135 The Physics of Stars Exam 2013 (Zane/Howarth) Answer ALL SIX questions from Section A, and ANY TWO questions from Section B The numbers in square brackets in the right-hand margin indicate the

More information

Opacity and Optical Depth

Opacity and Optical Depth Opacity and Optical Depth Absorption dominated intensity change can be written as di λ = κ λ ρ I λ ds with κ λ the absorption coefficient, or opacity The initial intensity I λ 0 of a light beam will be

More information

Lecture 3 Numerical Solutions to the Transport Equation

Lecture 3 Numerical Solutions to the Transport Equation Lecture 3 Numerical Solutions to the Transport Equation Introduction I There are numerous methods for solving the transport problem numerically. First we must recognize that we need to solve two problems:

More information

Energy transport: convection

Energy transport: convection Outline Introduction: Modern astronomy and the power of quantitative spectroscopy Basic assumptions for classic stellar atmospheres: geometry, hydrostatic equilibrium, conservation of momentum-mass-energy,

More information

Lecture 2 Solutions to the Transport Equation

Lecture 2 Solutions to the Transport Equation Lecture 2 Solutions to the Transport Equation Equation along a ray I In general we can solve the static transfer equation along a ray in some particular direction. Since photons move in straight lines

More information

Stars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11)

Stars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11) Stars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11) Kenneth Wood, Room 316 kw25@st-andrews.ac.uk http://www-star.st-and.ac.uk/~kw25 What is a Stellar Atmosphere? Transition from dense

More information

RADIATIVE TRANSFER IN AXIAL SYMMETRY

RADIATIVE TRANSFER IN AXIAL SYMMETRY Title : will be set by the publisher Editors : will be set by the publisher EAS Publications Series, Vol.?, 26 RADIATIVE TRANSFER IN AXIAL SYMMETRY Daniela Korčáková and Jiří Kubát Abstract. We present

More information

Radiative Transfer Plane-Parallel Frequency-Dependent

Radiative Transfer Plane-Parallel Frequency-Dependent 4 Radiative Transfer Plane-Parallel Frequency-Dependent variables I ν J ν H ν K ν in astronomy 4. Basic Equations After Fukue, J. 2, PASJ, 63, in press We here assume the followings: i) The disk is steady

More information

TRANSFER OF RADIATION

TRANSFER OF RADIATION TRANSFER OF RADIATION Under LTE Local Thermodynamic Equilibrium) condition radiation has a Planck black body) distribution. Radiation energy density is given as U r,ν = 8πh c 3 ν 3, LTE), tr.1) e hν/kt

More information

2. Basic Assumptions for Stellar Atmospheres

2. Basic Assumptions for Stellar Atmospheres 2. Basic Assumptions for Stellar Atmospheres 1. geometry, stationarity 2. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres!

More information

arxiv: v1 [astro-ph.sr] 14 May 2010

arxiv: v1 [astro-ph.sr] 14 May 2010 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres Editors : will be set by the publisher EAS Publications Series, Vol.?, 2010 arxiv:1005.2458v1 [astro-ph.sr] 14 May 2010 STATISTICAL EQUILIBRIUM

More information

7. Non-LTE basic concepts

7. Non-LTE basic concepts 7. Non-LTE basic concepts LTE vs NLTE occupation numbers rate equation transition probabilities: collisional and radiative examples: hot stars, A supergiants 10/13/2003 Spring 2016 LTE LTE vs NLTE each

More information

Stellar atmospheres: an overview

Stellar atmospheres: an overview Stellar atmospheres: an overview Core M = 2x10 33 g R = 7x10 10 cm 50 M o 20 R o L = 4x10 33 erg/s 10 6 L o 10 4 (PN) 10 6 (HII) 10 12 (QSO) L o Photosphere Envelope Chromosphere/Corona R = 200 km ~ 3x10

More information

2. Basic assumptions for stellar atmospheres

2. Basic assumptions for stellar atmospheres . Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres

More information

NLTE solar flare models with stationary velocity fields

NLTE solar flare models with stationary velocity fields ASTRONOMY & ASTROPHYSICS FEBRUARY I 1998, PAGE 607 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 127, 607-618 (1998) NLTE solar flare models with stationary velocity fields T. Nejezchleba Astronomical

More information

7. Non-LTE basic concepts

7. Non-LTE basic concepts 7. Non-LTE basic concepts LTE vs NLTE occupation numbers rate equation transition probabilities: collisional and radiative examples: hot stars, A supergiants 1 Equilibrium: LTE vs NLTE LTE each volume

More information

(c) Sketch the ratio of electron to gas pressure for main sequence stars versus effective temperature. [1.5]

(c) Sketch the ratio of electron to gas pressure for main sequence stars versus effective temperature. [1.5] 1. (a) The Saha equation may be written in the form N + n e N = C u+ u T 3/2 exp ( ) χ kt where C = 4.83 1 21 m 3. Discuss its importance in the study of stellar atmospheres. Carefully explain the meaning

More information

In this method, one defines

In this method, one defines 1 Feautrier s Method for Radiative Transfer This project is to use Feautrier s method to solve the Gray Atmosphere problem for a plane-parallel atmosphere Although we have developed alternate techniques

More information

Radiative Transfer in Axial Symmetry

Radiative Transfer in Axial Symmetry Radiative Transfer in Axial Symmetry Daniela Korčáková kor@sunstel.asu.cas.cz Astronomical Institute of the Academy of Sciences, CZ collaborators: T. Nagel, K. Werner, V. Suleymanov, V. Votruba, J. Kubát

More information

Model Photospheres with Accelerated Lambda Iteration

Model Photospheres with Accelerated Lambda Iteration Stellar Atmosphere Modeling ASP Conference Series, Vol. xxx, 23 I. Hubeny, D. Mihalas, K. Werner Model Photospheres with Accelerated Lambda Iteration Klaus Werner, Jochen L. Deetjen, Stefan Dreizler, Thorsten

More information

The Stellar Opacity. F ν = D U = 1 3 vl n = 1 3. and that, when integrated over all energies,

The Stellar Opacity. F ν = D U = 1 3 vl n = 1 3. and that, when integrated over all energies, The Stellar Opacity The mean absorption coefficient, κ, is not a constant; it is dependent on frequency, and is therefore frequently written as κ ν. Inside a star, several different sources of opacity

More information

SIMPLE RADIATIVE TRANSFER

SIMPLE RADIATIVE TRANSFER ASTR 511/O Connell Lec 4 1 SIMPLE RADIATIVE TRANSFER The theory of radiative transfer provides the means for determining the emergent EM spectrum of a cosmic source and also for describing the effects

More information

Modelling stellar atmospheres with full Zeeman treatment

Modelling stellar atmospheres with full Zeeman treatment 1 / 16 Modelling stellar atmospheres with full Zeeman treatment Katharina M. Bischof, Martin J. Stift M. J. Stift s Supercomputing Group FWF project P16003 Institute f. Astronomy Vienna, Austria CP#AP

More information

Atomic Physics 3 ASTR 2110 Sarazin

Atomic Physics 3 ASTR 2110 Sarazin Atomic Physics 3 ASTR 2110 Sarazin Homework #5 Due Wednesday, October 4 due to fall break Test #1 Monday, October 9, 11-11:50 am Ruffner G006 (classroom) You may not consult the text, your notes, or any

More information

Astro 305 Lecture Notes Wayne Hu

Astro 305 Lecture Notes Wayne Hu Astro 305 Lecture Notes Wayne Hu Set 1: Radiative Transfer Radiation Observables From an empiricist s point of view there are 4 observables for radiation Energy Flux Direction Color Polarization Energy

More information

Equations of Stellar Structure

Equations of Stellar Structure Equations of Stellar Structure Stellar structure and evolution can be calculated via a series of differential equations involving mass, pressure, temperature, and density. For simplicity, we will assume

More information

2. Basic assumptions for stellar atmospheres

2. Basic assumptions for stellar atmospheres . Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres

More information

On the NLTE plane-parallel and spherically symmetric model atmospheres of helium rich central stars of planetary nebulae

On the NLTE plane-parallel and spherically symmetric model atmospheres of helium rich central stars of planetary nebulae Astron. Astrophys. 323, 524 528 (1997) ASTRONOMY AND ASTROPHYSICS Research Note On the NLTE plane-parallel and spherically symmetric model atmospheres of helium rich central stars of planetary nebulae

More information

Stellar Atmospheres: Basic Processes and Equations

Stellar Atmospheres: Basic Processes and Equations Stellar Atmospheres: Basic Processes and Equations Giovanni Catanzaro Abstract The content of this chapter is a very quick summary of key concepts that concern the interaction between photons created in

More information

9.1 Introduction. 9.2 Static Models STELLAR MODELS

9.1 Introduction. 9.2 Static Models STELLAR MODELS M. Pettini: Structure and Evolution of Stars Lecture 9 STELLAR MODELS 9.1 Introduction Stars are complex physical systems, but not too complex to be modelled numerically and, with some simplifying assumptions,

More information

2. Basic assumptions for stellar atmospheres

2. Basic assumptions for stellar atmospheres . Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres

More information

Lecture Notes on Radiation Transport for Spectroscopy

Lecture Notes on Radiation Transport for Spectroscopy Lecture Notes on Radiation Transport for Spectroscopy ICTP-IAEA School on Atomic Processes in Plasmas 27 February 3 March 217 Trieste, Italy LLNL-PRES-668311 This work was performed under the auspices

More information

Assignment 4 Solutions [Revision : 1.4]

Assignment 4 Solutions [Revision : 1.4] Assignment 4 Solutions [Revision : 1.4] Q9.7 We typically see a optical distance τ 2/3 through an opaque medium. Using τ = κρs, for constant κ = 0.03 m 2 kg 1 and ρ = 1.2 kgm 3, gives a physical distance

More information

Fundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres

Fundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial

More information

Topics ASTR 3730: Fall 2003

Topics ASTR 3730: Fall 2003 Topics Qualitative questions: might cover any of the main topics (since 2nd midterm: star formation, extrasolar planets, supernovae / neutron stars, black holes). Quantitative questions: worthwhile to

More information

Electromagnetic Spectra. AST443, Lecture 13 Stanimir Metchev

Electromagnetic Spectra. AST443, Lecture 13 Stanimir Metchev Electromagnetic Spectra AST443, Lecture 13 Stanimir Metchev Administrative Homework 2: problem 5.4 extension: until Mon, Nov 2 Reading: Bradt, chapter 11 Howell, chapter 6 Tenagra data: see bottom of Assignments

More information

11.1 Local Thermodynamic Equilibrium. 1. the electron and ion velocity distributions are Maxwellian,

11.1 Local Thermodynamic Equilibrium. 1. the electron and ion velocity distributions are Maxwellian, Section 11 LTE Basic treatments of stellar atmospheres adopt, as a starting point, the assumptions of Local Thermodynamic Equilibrium (LTE), and hydrostatic equilibrium. The former deals with the microscopic

More information

Radiative Transfer and Stellar Atmospheres

Radiative Transfer and Stellar Atmospheres Radiative Transfer and Stellar Atmospheres 4 lectures within the first IMPRS advanced course Joachim Puls Institute for Astronomy & Astrophysics, Munich Contents quantitative spectroscopy: the astrophysical

More information

arxiv:astro-ph/ v1 16 Jun 2005

arxiv:astro-ph/ v1 16 Jun 2005 Astronomy & Astrophysics manuscript no. 992 November 5, 208 (DOI: will be inserted by hand later) Radiative transfer in moving media II. Solution of the radiative transfer equation in axial symmetry Daniela

More information

XXIX IAC Winter School of Astrophysics

XXIX IAC Winter School of Astrophysics XXIX IAC Winter School of Astrophysics APPLICATIONS OF RADIATIVE TRANSFER TO STELLAR AND PLANETARY ATMOSPHERES Fundamental physical aspects of radiatie transfer: II.- Statistical Equilibrium Equations

More information

Fundamental Stellar Parameters

Fundamental Stellar Parameters Fundamental Stellar Parameters Radiative Transfer Specific Intensity, Radiative Flux and Stellar Luminosity Observed Flux, Emission and Absorption of Radiation Radiative Transfer Equation, Solution and

More information

VII. Hydrodynamic theory of stellar winds

VII. Hydrodynamic theory of stellar winds VII. Hydrodynamic theory of stellar winds observations winds exist everywhere in the HRD hydrodynamic theory needed to describe stellar atmospheres with winds Unified Model Atmospheres: - based on the

More information

ASTRONOMY QUALIFYING EXAM August Possibly Useful Quantities

ASTRONOMY QUALIFYING EXAM August Possibly Useful Quantities L = 3.9 x 10 33 erg s 1 M = 2 x 10 33 g M bol = 4.74 R = 7 x 10 10 cm 1 A.U. = 1.5 x 10 13 cm 1 pc = 3.26 l.y. = 3.1 x 10 18 cm a = 7.56 x 10 15 erg cm 3 K 4 c= 3.0 x 10 10 cm s 1 σ = ac/4 = 5.7 x 10 5

More information

Limb Darkening. Limb Darkening. Limb Darkening. Limb Darkening. Empirical Limb Darkening. Betelgeuse. At centre see hotter gas than at edges

Limb Darkening. Limb Darkening. Limb Darkening. Limb Darkening. Empirical Limb Darkening. Betelgeuse. At centre see hotter gas than at edges Limb Darkening Sun Betelgeuse Limb Darkening Stars are both redder and dimmer at the edges Sun Limb Darkening Betelgeuse Limb Darkening Can also be understood in terms of temperature within the solar photosphere.

More information

2. Stellar atmospheres: Structure

2. Stellar atmospheres: Structure 2. Stellar atmospheres: Structure 2.1. Assumptions Plane-parallel geometry Hydrostatic equilibrium, i.e. o no large-scale accelerations comparable to surface gravity o no dynamically significant mass loss

More information

Astronomy 421. Lecture 14: Stellar Atmospheres III

Astronomy 421. Lecture 14: Stellar Atmospheres III Astronomy 421 Lecture 14: Stellar Atmospheres III 1 Lecture 14 - Key concepts: Spectral line widths and shapes Curve of growth 2 There exists a stronger jump, the Lyman limit, occurring at the wavelength

More information

Overview of Astronomical Concepts III. Stellar Atmospheres; Spectroscopy. PHY 688, Lecture 5 Stanimir Metchev

Overview of Astronomical Concepts III. Stellar Atmospheres; Spectroscopy. PHY 688, Lecture 5 Stanimir Metchev Overview of Astronomical Concepts III. Stellar Atmospheres; Spectroscopy PHY 688, Lecture 5 Stanimir Metchev Outline Review of previous lecture Stellar atmospheres spectral lines line profiles; broadening

More information

Astronomy. Astrophysics. On the importance of the wind emission to the optical continuum of OB supergiants. M. Kraus 1, J. Kubát 1, and J.

Astronomy. Astrophysics. On the importance of the wind emission to the optical continuum of OB supergiants. M. Kraus 1, J. Kubát 1, and J. A&A 481, 499 506 (2008) DOI: 10.1051/0004-6361:20078991 c ESO 2008 Astronomy & Astrophysics On the importance of the wind emission to the optical continuum of OB supergiants M. Kraus 1, J. Kubát 1, and

More information

Model Atmospheres. Model Atmosphere Assumptions

Model Atmospheres. Model Atmosphere Assumptions Model Atmospheres Problem: Construct a numerical model of the atmosphere to estimate (a) Variation of physical variables (T, P) with depth (b) Emergent spectrum in continuum and lines Compare calculated

More information

Introduction to the School

Introduction to the School Lucio Crivellari Instituto de Astrofísica de Canarias D.pto de Astrofísica, Universidad de La Laguna & INAF Osservatorio Astronomico di Trieste (Italy) Introduction to the School 10/11/17 1 Setting the

More information

Stellar Structure. Observationally, we can determine: Can we explain all these observations?

Stellar Structure. Observationally, we can determine: Can we explain all these observations? Stellar Structure Observationally, we can determine: Flux Mass Distance Luminosity Temperature Radius Spectral Type Composition Can we explain all these observations? Stellar Structure Plan: Use our general

More information

SPA7023P/SPA7023U/ASTM109 Stellar Structure and Evolution Duration: 2.5 hours

SPA7023P/SPA7023U/ASTM109 Stellar Structure and Evolution Duration: 2.5 hours MSc/MSci Examination Day 28th April 2015 18:30 21:00 SPA7023P/SPA7023U/ASTM109 Stellar Structure and Evolution Duration: 2.5 hours YOU ARE NOT PERMITTED TO READ THE CONTENTS OF THIS QUESTION PAPER UNTIL

More information

2. NOTES ON RADIATIVE TRANSFER The specific intensity I ν

2. NOTES ON RADIATIVE TRANSFER The specific intensity I ν 1 2. NOTES ON RADIATIVE TRANSFER 2.1. The specific intensity I ν Let f(x, p) be the photon distribution function in phase space, summed over the two polarization states. Then fdxdp is the number of photons

More information

Stellar Spectra ASTR 2120 Sarazin. Solar Spectrum

Stellar Spectra ASTR 2120 Sarazin. Solar Spectrum Stellar Spectra ASTR 2120 Sarazin Solar Spectrum Solar Prominence Sep. 14, 1999 Solar Activity Due to rotation, convection, and magnetic field (Section 7.2 review) Charged Particles in Magnetic Fields

More information

Atomic Physics ASTR 2110 Sarazin

Atomic Physics ASTR 2110 Sarazin Atomic Physics ASTR 2110 Sarazin Homework #5 Due Wednesday, October 4 due to fall break Test #1 Monday, October 9, 11-11:50 am Ruffner G006 (classroom) You may not consult the text, your notes, or any

More information

Beer-Lambert (cont.)

Beer-Lambert (cont.) The Beer-Lambert Law: Optical Depth Consider the following process: F(x) Absorbed flux df abs F(x + dx) Scattered flux df scat x x + dx The absorption or scattering of radiation by an optically active

More information

Radiation Transfer and Radiation Hydrodynamics

Radiation Transfer and Radiation Hydrodynamics 2018/02/02 Winter School @ SHAO Radiation Transfer and Radiation Hydrodynamics Kengo TOMIDA ( 富田賢吾, Osaka University) Some Slides are provided by Yan-Fei Jiang (UCSB) Textbooks Rybicki & Lightman, Radiative

More information

Black Hole Formation in the Early Universe

Black Hole Formation in the Early Universe 1 Black Hole Formation in the Early Universe Luo Yang (Osaka University) XC30 2 1. Introduction Observations of high-redshift quasars reveal that some supermassive black holes (SMBH) with masses exceeding

More information

Radiative transfer in planetary atmospheres

Radiative transfer in planetary atmospheres Chapter 9 Radiative transfer in planetary atmospheres Warning: This chapter is just a very rough introduction to planetary atmospheres. Since the discovery of planets around other stars than the sun, the

More information

Spectroscopy Lecture 2

Spectroscopy Lecture 2 Spectroscopy Lecture 2 I. Atomic excitation and ionization II. Radiation Terms III. Absorption and emission coefficients IV. Einstein coefficients V. Black Body radiation I. Atomic excitation and ionization

More information

Stellar structure Conservation of mass Conservation of energy Equation of hydrostatic equilibrium Equation of energy transport Equation of state

Stellar structure Conservation of mass Conservation of energy Equation of hydrostatic equilibrium Equation of energy transport Equation of state Stellar structure For an isolated, static, spherically symmetric star, four basic laws / equations needed to describe structure: Conservation of mass Conservation of energy (at each radius, the change

More information

Linear Theory of Stellar Pulsation

Linear Theory of Stellar Pulsation Linear Theory of Stellar Pulsation Mir Emad Aghili Department of Physics and Astronomy, University of Mississippi, University, MS, 38677 Introduction The first pulsating star to be discovered was omicron

More information

Lecture Notes: Basic Equations

Lecture Notes: Basic Equations Basic Equations Lecture Notes: INTRODUCTION TO NON-LTE RADIATIVE TRANSFER AND ATMOSPHERIC MODELING Eugene H. Avrett Harvard-Smithsonian Center for Astrophysics July 2008 The specific intensity of radiation

More information

Nearly everything that we know about stars comes from the photons they emit into space.

Nearly everything that we know about stars comes from the photons they emit into space. STELLAR ATMOSPHERES Nearly everything that we know about stars comes from the photons they emit into space. In the interior, we saw the photon distribution had a dominant isotropic part I 0, and just a

More information

Section 11.5 and Problem Radiative Transfer. from. Astronomy Methods A Physical Approach to Astronomical Observations Pages , 377

Section 11.5 and Problem Radiative Transfer. from. Astronomy Methods A Physical Approach to Astronomical Observations Pages , 377 Section 11.5 and Problem 11.51 Radiative Transfer from Astronomy Methods A Physical Approach to Astronomical Observations Pages 365-375, 377 Cambridge University Press 24 by Hale Bradt Hale Bradt 24 11.5

More information

Observational Appearance of Black Hole Wind Effect of Electron Scattering

Observational Appearance of Black Hole Wind Effect of Electron Scattering Observational Appearance of Black Hole Wind Effect of Electron Scattering Kazuyuki OGURA Astronomical Institute Osaka Kyoiku Univ. 29 Jun 2013 Meeting of BH Horizon Project @Nagoya Univ. Contents Introduction

More information

Environment of the Radiation Field ...

Environment of the Radiation Field ... Copyright (2003) Geroge W. Collins, II 11 Environment of the Radiation Field... Thus far, we have said little or nothing about the gas through which the radiation is flowing. This constitutes the second

More information

SISD Training Lectures in Spectroscopy

SISD Training Lectures in Spectroscopy SISD Training Lectures in Spectroscopy Anatomy of a Spectrum Visual Spectrum of the Sun Blue Spectrum of the Sun Morphological Features in Spectra λ 2 Line Flux = Fλ dλ λ1 (Units: erg s -1 cm -2 ) Continuum

More information

Non-LTE models for synthetic spectra of Type Ia supernovae

Non-LTE models for synthetic spectra of Type Ia supernovae A&A 569, A62 (24) DOI:.5/4-636/2322265 c ESO 24 Astronomy & Astrophysics Non-LTE models for synthetic spectra of Type Ia supernovae IV. A modified Feautrier scheme for opacity-sampled pseudo-continua at

More information

Radiative Transfer and Molecular Lines Sagan Workshop 2009

Radiative Transfer and Molecular Lines Sagan Workshop 2009 Radiative Transfer and Molecular Lines Sagan Workshop 2009 Sara Seager Trent Schindler Trent Schindler MIT Lecture Contents Overview of Equations for Planetary Atmospheres Radiative Transfer Thermal Inversions

More information

Stellar Winds: Mechanisms and Dynamics

Stellar Winds: Mechanisms and Dynamics Astrofysikalisk dynamik, VT 010 Stellar Winds: Mechanisms and Dynamics Lecture Notes Susanne Höfner Department of Physics and Astronomy Uppsala University 1 Most stars have a stellar wind, i.e. and outflow

More information

Review from last class:

Review from last class: Review from last class: Properties of photons Flux and luminosity, apparent magnitude and absolute magnitude, colors Spectroscopic observations. Doppler s effect and applications Distance measurements

More information

Radiative transfer theory (Pure gas atmosphere)

Radiative transfer theory (Pure gas atmosphere) Radiative transfer theory (Pure gas atmosphere Bruno Carli Radiative transfer equation The radiative transfer equation states that the specific intensity of radiation I( (defined as the energy flux per

More information

PHYS 231 Lecture Notes Week 3

PHYS 231 Lecture Notes Week 3 PHYS 231 Lecture Notes Week 3 Reading from Maoz (2 nd edition): Chapter 2, Sec. 3.1, 3.2 A lot of the material presented in class this week is well covered in Maoz, and we simply reference the book, with

More information

Substellar Atmospheres. PHY 688, Lecture 18 Mar 9, 2009

Substellar Atmospheres. PHY 688, Lecture 18 Mar 9, 2009 Substellar Atmospheres PHY 688, Lecture 18 Mar 9, 2009 Outline Review of previous lecture the Kepler mission launched successfully results P < 1 month planets by September 09 giant planet interiors comparison

More information

The structure and evolution of stars. Introduction and recap

The structure and evolution of stars. Introduction and recap The structure and evolution of stars Lecture 3: The equations of stellar structure 1 Introduction and recap For our stars which are isolated, static, and spherically symmetric there are four basic equations

More information

Stellar Atmospheres. University of Denver, Department of Physics and Astronomy. Physics 2052 Stellar Physics, Winter 2008.

Stellar Atmospheres. University of Denver, Department of Physics and Astronomy. Physics 2052 Stellar Physics, Winter 2008. Stellar Atmospheres University of Denver, Department of Physics and Astronomy Physics 2052 Stellar Physics, Winter 2008 By Kathy Geise Introduction A star does not have a solid surface, so the definition

More information

Summary of Results. 1. Mass Conserva-on. dr dm = 1. dm dr = 4πr 2 ρ. 4πr 2 ρ

Summary of Results. 1. Mass Conserva-on. dr dm = 1. dm dr = 4πr 2 ρ. 4πr 2 ρ Summary of Results 1. Mass Conserva-on dm dr = 4πr 2 ρ eulerian dr dm = 1 4πr 2 ρ lagrangian Summary of Results 2. Hydrosta-c Equilibrium dp dr = ρ Gm eulerian r 2 dp dm = Gm 4πr 4 lagrangian Need equa-on

More information

I ν. di ν. = α ν. = (ndads) σ ν da α ν. = nσ ν = ρκ ν

I ν. di ν. = α ν. = (ndads) σ ν da α ν. = nσ ν = ρκ ν Absorption Consider a beam passing through an absorbing medium. Define the absorption coefficient, α ν, by ie the fractional loss in intensity in travelling a distance ds is α ν ds (convention: positive

More information

Ay Fall 2004 Lecture 6 (given by Tony Travouillon)

Ay Fall 2004 Lecture 6 (given by Tony Travouillon) Ay 122 - Fall 2004 Lecture 6 (given by Tony Travouillon) Stellar atmospheres, classification of stellar spectra (Many slides c/o Phil Armitage) Formation of spectral lines: 1.excitation Two key questions:

More information

6. Stellar spectra. excitation and ionization, Saha s equation stellar spectral classification Balmer jump, H -

6. Stellar spectra. excitation and ionization, Saha s equation stellar spectral classification Balmer jump, H - 6. Stellar spectra excitation and ionization, Saha s equation stellar spectral classification Balmer jump, H - 1 Occupation numbers: LTE case Absorption coefficient: = n i calculation of occupation numbers

More information

23 Astrophysics 23.5 Ionization of the Interstellar Gas near a Star

23 Astrophysics 23.5 Ionization of the Interstellar Gas near a Star 23 Astrophysics 23.5 Ionization of the Interstellar Gas near a Star (8 units) No knowledge of Astrophysics is assumed or required: all relevant equations are defined and explained in the project itself.

More information

= m H 2. The mean mass is 1/2 the mass of the hydrogen atom. Since, the m e << m p, the mean mass of hydrogen is close to 1/2 the mass of a proton.

= m H 2. The mean mass is 1/2 the mass of the hydrogen atom. Since, the m e << m p, the mean mass of hydrogen is close to 1/2 the mass of a proton. From Last Time: Mean mass m m = m e + m p 2 = m H 2 The mean mass is 1/2 the mass of the hydrogen atom. Since, the m e

More information

ASTM109 Stellar Structure and Evolution Duration: 2.5 hours

ASTM109 Stellar Structure and Evolution Duration: 2.5 hours MSc Examination Day 15th May 2014 14:30 17:00 ASTM109 Stellar Structure and Evolution Duration: 2.5 hours YOU ARE NOT PERMITTED TO READ THE CONTENTS OF THIS QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY

More information

COX & GIULI'S PRINCIPLES OF STELLAR STRUCTURE

COX & GIULI'S PRINCIPLES OF STELLAR STRUCTURE COX & GIULI'S PRINCIPLES OF STELLAR STRUCTURE Extended Second Edition A. Weiss, W. Hillebrandt, H.-C. Thomas and H. Ritter Max-Planck-lnstitut fur Astrophysik, Garching, Germany C S P CONTENTS PREFACE

More information

Physics 160: Stellar Astrophysics. Midterm Exam. 27 October 2011 INSTRUCTIONS READ ME!

Physics 160: Stellar Astrophysics. Midterm Exam. 27 October 2011 INSTRUCTIONS READ ME! Physics 160: Stellar Astrophysics 27 October 2011 Name: S O L U T I O N S Student ID #: INSTRUCTIONS READ ME! 1. There are 4 questions on the exam; complete at least 3 of them. 2. You have 80 minutes to

More information

ASTRONOMY AND ASTROPHYSICS Theoretical X-ray spectra of hot H-rich white dwarfs. Impact of new partition functions of iron, Fe V through Fe VII

ASTRONOMY AND ASTROPHYSICS Theoretical X-ray spectra of hot H-rich white dwarfs. Impact of new partition functions of iron, Fe V through Fe VII Astron. Astrophys. 343, 531 535 (1999) ASTRONOMY AND ASTROPHYSICS Theoretical X-ray spectra of hot H-rich white dwarfs. Impact of new partition functions of iron, Fe V through Fe VII J. Madej 1, J. Halenka

More information

The structure of radiative shock waves. V. Hydrogen emission lines

The structure of radiative shock waves. V. Hydrogen emission lines A&A 420, 423 435 (2004) DOI: 10.1051/0004-6361:20040992 c ESO 2004 Astronomy & Astrophysics The structure of radiative shock waves V. Hydrogen emission lines Y. A. Fadeyev 1 and D. Gillet 2 1 Institute

More information

arxiv: v1 [astro-ph.sr] 2 Jan 2009

arxiv: v1 [astro-ph.sr] 2 Jan 2009 Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 22 January 2009 (MN LATEX style file v2.2) NLTE models of line-driven stellar winds III. Influence of X-ray radiation on wind structure of O stars

More information

Bremsstrahlung. Rybicki & Lightman Chapter 5. Free-free Emission Braking Radiation

Bremsstrahlung. Rybicki & Lightman Chapter 5. Free-free Emission Braking Radiation Bremsstrahlung Rybicki & Lightman Chapter 5 Bremsstrahlung Free-free Emission Braking Radiation Radiation due to acceleration of charged particle by the Coulomb field of another charge. Relevant for (i)

More information

arxiv: v1 [astro-ph.sr] 13 Jun 2014

arxiv: v1 [astro-ph.sr] 13 Jun 2014 Current status of NLTE analysis of stellar atmospheres Jiří Kubát arxiv:1406.3559v1 [astro-ph.sr] 13 Jun 2014 Abstract Various available codes for NLTE modeling and analysis of hot star spectra are reviewed.

More information

Stellar Interiors. Hydrostatic Equilibrium. PHY stellar-structures - J. Hedberg

Stellar Interiors. Hydrostatic Equilibrium. PHY stellar-structures - J. Hedberg Stellar Interiors. Hydrostatic Equilibrium 2. Mass continuity 3. Equation of State. The pressure integral 4. Stellar Energy Sources. Where does it come from? 5. Intro to Nuclear Reactions. Fission 2. Fusion

More information

PHYS-633: Introduction to Stellar Astrophysics

PHYS-633: Introduction to Stellar Astrophysics PHYS-633: Introduction to Stellar Astrophysics Version of July 7, 2015 Stan Owocki Contents 1 Course Overview: Stellar Atmospheres vs. Interiors 6 1.1 Differences between Atmospheres and Interiors.................

More information

Solar and Stellar Atmospheric Modeling Using the Pandora Computer Program

Solar and Stellar Atmospheric Modeling Using the Pandora Computer Program **TITLE** ASP Conference Series, Vol. **VOLUME**, **PUBLICATION YEAR** **EDITORS** Solar and Stellar Atmospheric Modeling Using the Pandora Computer Program Eugene H. Avrett and Rudolf Loeser Smithsonian

More information

P M 2 R 4. (3) To determine the luminosity, we now turn to the radiative diffusion equation,

P M 2 R 4. (3) To determine the luminosity, we now turn to the radiative diffusion equation, Astronomy 715 Final Exam Solutions Question 1 (i). The equation of hydrostatic equilibrium is dp dr GM r r 2 ρ. (1) This corresponds to the scaling P M R ρ, (2) R2 where P and rho represent the central

More information

Computing radiative equilibria with Monte Carlo techniques

Computing radiative equilibria with Monte Carlo techniques Astron. Astrophys. 344, 282 288 (999) ASTRONOMY AND ASTROPHYSICS Computing radiative equilibria with Monte Carlo techniques L.B. Lucy Space Telescope European Coordinating Facility, Karl-Schwarzschild-Strasse

More information

If light travels past a system faster than the time scale for which the system evolves then t I ν = 0 and we have then

If light travels past a system faster than the time scale for which the system evolves then t I ν = 0 and we have then 6 LECTURE 2 Equation of Radiative Transfer Condition that I ν is constant along rays means that di ν /dt = 0 = t I ν + ck I ν, (29) where ck = di ν /ds is the ray-path derivative. This is equation is the

More information

2 The Radiative Transfer Equation

2 The Radiative Transfer Equation 9 The Radiative Transfer Equation. Radiative transfer without absorption and scattering Free space or homogeneous space I (r,,) I (r,,) r -r d da da Figure.: Following a pencil of radiation in free space

More information

Introduction to Solar Radiative Transfer II Non-LTE Radiative Transfer

Introduction to Solar Radiative Transfer II Non-LTE Radiative Transfer Introduction to olar Radiative Transfer II Non-LTE Radiative Transfer Han Uitenbroek National olar Observatory/acramento Peak unspot NM Overview I Basic Radiative Transfer Intensity, emission, absorption,

More information