Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres
|
|
- Gerard McBride
- 5 years ago
- Views:
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 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 informationPHAS3135 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 informationOpacity 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 informationLecture 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 informationEnergy 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 informationLecture 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 informationStars 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 informationRADIATIVE 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 informationRadiative 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 informationTRANSFER 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 information2. 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 informationarxiv: 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 information7. 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 informationStellar 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 information2. 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 informationNLTE 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 information7. 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]
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 informationIn 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 informationRadiative 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 informationModel 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 informationThe 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 informationSIMPLE 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 informationModelling 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 informationAtomic 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 informationAstro 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 informationEquations 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 information2. 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 informationOn 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 informationStellar 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 information9.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 information2. 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 informationLecture 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 informationAssignment 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 informationFundamental 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 informationTopics 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 informationElectromagnetic 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 information11.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 informationRadiative 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 informationarxiv: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 informationXXIX 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 informationFundamental 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 informationVII. 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 informationASTRONOMY 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 informationLimb 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 information2. 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 informationAstronomy 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 informationOverview 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 informationAstronomy. 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 informationModel 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 informationIntroduction 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 informationStellar 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 informationSPA7023P/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 information2. 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 informationStellar 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 informationAtomic 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 informationBeer-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 informationRadiation 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 informationBlack 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 informationRadiative 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 informationSpectroscopy 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 informationStellar 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 informationLinear 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 informationLecture 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 informationNearly 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 informationSection 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 informationObservational 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 informationEnvironment 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 informationSISD 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 informationNon-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 informationRadiative 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 informationStellar 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 informationReview 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 informationRadiative 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 informationPHYS 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 informationSubstellar 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 informationThe 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 informationStellar 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 informationSummary 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 informationI ν. 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 informationAy 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 information6. 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 information23 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.
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 informationASTM109 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 informationCOX & 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 informationPhysics 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 informationASTRONOMY 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 informationThe 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 informationarxiv: 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 informationBremsstrahlung. 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 informationarxiv: 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 informationStellar 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 informationPHYS-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 informationSolar 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 informationP 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 informationComputing 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 informationIf 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 information2 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 informationIntroduction 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