Solution of the radiative transfer equation in NLTE stellar atmospheres
|
|
- Stella Glenn
- 5 years ago
- Views:
Transcription
1 Solution of the radiative transfer equation in NLTE stellar atmospheres Jiří Kubát Astronomický ústav AV ČR Ondřejov Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 1
2 Outline 1. Basic equations for stellar atmospheres 2. Radiative transfer equation 3. Equations of statistical equilibrium 4. Discretization 5. Formal solution of RTE 6. Lambda iteration 7. Complete linearization 8. Accelerated lambda iteration Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 2
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) radiatively driven stellar wind Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 3
4 Basic equations assume static atmosphere ( v = 0) stationary atmosphere ( / t = 0) 1-dimensional atmosphere given temperature structure T ( r) given density structure ρ( r) equations to be solved (numerically): radiative transfer equation equations of statistical equilibrium Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 4
5 Basic equations assume static atmosphere ( v = 0) stationary atmosphere ( / t = 0) 1-dimensional atmosphere given temperature structure T ( r) given density structure ρ( r) equations to be solved (numerically): radiative transfer equation equations of statistical equilibrium Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 4
6 Basic equations assume static atmosphere ( v = 0) stationary atmosphere ( / t = 0) 1-dimensional atmosphere given temperature structure T ( r) given density structure ρ( r) equations to be solved (numerically): radiative transfer equation equations of statistical equilibrium account for full angle and frequency dependence (I (µ, ) I µ ) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 4
7 Radiative transfer equation µ di µ(z) dz = η (z) χ (z)i µ (z) = χ [S (z) I µ (z)] η (z) emissivity χ (z) opacity S (z) = η (z) χ (z) source function Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 5
8 Radiative transfer equation µ di µ(z) dz = η (z) χ (z)i µ (z) = χ [S (z) I µ (z)] introduce dτ = χ dz optical depth J = 1 1 I 2 1 µ dµ mean intensity K = µ2 I µ dµ f = K /J variable Eddington factor (Auer & Mihalas 1970) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 5
9 Radiative transfer equation µ di µ(z) dz = η (z) χ (z)i µ (z) = χ [S (z) I µ (z)] introduce dτ = χ dz optical depth J = 1 1 I 2 1 µ dµ mean intensity K = µ2 I µ dµ f = K /J variable Eddington factor (Auer & Mihalas 1970) 2nd order equation d 2 [f (z)j (z)] dτ 2 = J (z) S (z) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 5
10 Radiative transfer equation 2nd order equation d 2 [f (z)j (z)] dτ 2 = J (z) S (z) + boundary conditions d [f (z)j (z)] = g (z) H dτ upper d [f (z)j (z)] = H + + g (z) dτ lower g = R 1 0 µ (I + µ I µ) dµ J Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 5
11 Radiative transfer equation 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 ] Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 6
12 Radiative transfer equation 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 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 6
13 Equations of statistical equilibrium in LTE: n i = f (n e, T ) Saha-Boltzmann distribution Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 7
14 Equations of statistical equilibrium outside LTE: n i = f (n e, T, J ) statistical equilibrium for i = 1,... NL l i {n l [R li + C li ] n i (R il + C il )} = 0 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 7
15 Equations of statistical equilibrium outside LTE: n i = f (n e, T, J ) statistical equilibrium for i = 1,... NL l i {n l [R li + C li ] n i (R il + C il )} = 0 collisional rates n i C il = n i n e q il (T ), Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 7
16 Equations of statistical equilibrium outside LTE: n i = f (n e, T, J ) statistical equilibrium for i = 1,... NL l i radiative rates {n l [R li + C li ] n i (R il + C il )} = 0 α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 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 7
17 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 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 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 7
18 Final radiative transfer equation d 2 [f (z)j (z)] dτ 2 = J (z) S (z, J ) S = η χ f = K J K = χ = f χ (, n e, T, n i ) Z 1 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 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 8
19 Two-level atom d 2 [f (z)j (z)] dτ 2 = J (z) S (z, J ) the source function may be written as (Mihalas 1978) S = (1 ε) ϕ J d + εb ε = ε 1 + ε ε = C 21 ( 1 e h/kt ) A 21 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 9
20 Discretization J continuous function of z, depth (z d) equidistant in log m or log τ some 4 5 depth points per decade frequency ( n) not equidistant need to resolve lines continuum edges angle (µ m) 3 directions all quantities expressed as J dn Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 10
21 Discretized radiative transfer equation for a frequency point n and angle m, depth points d = 1,..., ND a d f d 1 J d 1 + (b d f d + 1) J d + c d f d+1 J d+1 = S d (b 1 f 1 + g 1 + 1) J 1 + c 1 f 2 J 2 = H + S 1 a ND f ND 1 J ND 1 + (b ND f ND + g ND + 1) J ND = H + + S ND» 1 1 a d = τ 2 d + τ 12 d+ τ 12 d 12» 1 1 c d = τ 2 d + τ 12 d+ τ 12 d+ 12 b d = a d c d. Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 11
22 Discretized frequency integration for lines, ensure that J() d φ() d = NF n=1 NF n=1 w n J n w n φ n = 1 if not, it is obligatory to renormalize w n for continuum edges place 2 frequency points (on both sides) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 12
23 Discretized angle integration I(µ) dµ NF m=1 w m I m 3 points sufficient if Gaussian quadrature is used also necessary to check normalization Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 13
24 Discretized equations of statistical equilibrium for all d = 1,..., ND, i = 1,..., NL {(n l ) d [(R li ) d + (C li ) d ] (n i ) d [(R il ) d + (C il ) d ]} = 0 l i (n i ) d (R il ) d = (n i ) d 4π NF n=1 (n l ) d (R li ) d = (n l ) d g i g j 4π w n (α il ) n h n NF n=1 J dn w n (α il ) n h n ( 2h 3 n c 2 + J dn ) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 14
25 Formal solution of the RTE solution for given χ and η (given S ) relatively simple d 2 [f (z)j (z)] dτ 2 = J (z) S (z) 2nd order differential equation same numerical scheme as for the Feautrier solution crucial for the total accuracy of the whole problem solution Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 15
26 Λ iteration radiative transfer equation for J d 2 [f (z)j (z)] dτ 2 = J (z) S (z) 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 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 16
27 Complete linearization to stellar atmospheres introduced by Auer & Mihalas (1969) for the case of NLTE model atmospheres solution of equations: radiative transfer (J ) hydrostatic equilibrium (ρ) radiative equilibrium (T ) statistical equilibrium (n i ) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 17
28 Complete linearization to stellar atmospheres introduced by Auer & Mihalas (1969) for the case of NLTE model atmospheres for a restricted problem NLTE line formation (Auer 1973) solution of equations: radiative transfer (J ) statistical equilibrium (n i ) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 17
29 Complete linearization to stellar atmospheres introduced by Auer & Mihalas (1969) for the case of NLTE model atmospheres for a restricted problem NLTE line formation (Auer 1973) solution of equations: radiative transfer (J ) statistical equilibrium (n i ) vector of solution ψ = (J 1,..., J NF, n 1,..., n NL ) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 17
30 Complete linearization vector of solution ψ = (J 1,..., J NF, n 1,..., n NL ), dimension NF + NL formally F ( ) ψ = 0 current estimate ψ 0 correct solution ψ = ψ 0 + δ ψ corrections δ ψ = [ F ψ ( ψ0 ) ] 1 matrix NL + NF for each depth point d F ( ψ0 ) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 18
31 Complete linearization radiative transfer equation (from Mihalas, 1978, Stellar atmospheres) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 19
32 Complete linearization radiative transfer equation J = (J 1,..., J NF ) A d δ J d 1 + B d δ J d + C d δ J d+1 = L d δs nd = NL l=1 S n n l (δn l ) d d Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 19
33 Complete linearization radiative transfer equation J = (J 1,..., J NF ) A d δ J d 1 + B d δ J d + C d δ J d+1 = L d δs nd = NL l=1 S n n l (δn l ) d d A d δj n,d 1 + B d δj n,d + C d δj n,d+1 + D d δn l,d = L n,d Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 19
34 Complete linearization radiative transfer equation A d δj n,d 1 + B d δj n,d + C d δj n,d+1 + D d δn l,d = L n,d statistical equilibrium A n = b, n = (n 1,..., n NL ) [ ] A b n n n + A d δ n d + [ A J n,d n b J n,d ] δj n,d = b d A d n d E d δn l,d + F d δj n,d = K l,d Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 19
35 Complete linearization radiative transfer equation A d δj n,d 1 + B d δj n,d + C d δj n,d+1 + D d δn l,d = L n,d statistical equilibrium E d δn l,d + F d δj n,d = K l,d Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 19
36 Complete linearization radiative transfer equation A d δj n,d 1 + B d δj n,d + C d δj n,d+1 + D d δn l,d = L n,d statistical equilibrium E d δn l,d + F d δj n,d = K l,d n l J n = δn l,d = NL r=1 A 1 lr NF n=1 [ n l J n b r J n δ(j n ) d d NL s=1 A rs J n n s ] Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 19
37 Complete linearization radiative transfer + statistical equilibrium A d δ J d 1 + B d δ J d + C d δ J d+1 = L d J d = (J 1,..., J NF ), d = 1,..., ND A 1 = 0, C ND = 0 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 20
38 Complete linearization reformulation using net radiative bracket (Auer & Heasley 1976) t transition index (δz t ) d = (n i ) d (δr il ) d (n l ) d (δr li ) d δn i = t n i Z t δz t δ Z t + R t n i + S t δ n l = L t ( ) n i n l 1 + R t + S t δz t = Z t Z L t R t t t t n i Z t δz t S t t t n l Z t δz t system is solved for (Z t ) d Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 21
39 Accelerated lambda iteration huge matrices in complete linearization laborious and computationally expensive simple Λ-iteration convergence problems way out > Accelerated lambda iteration (ALI) (Cannon 1973) formal solution J = Λ S Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 22
40 Accelerated lambda iteration huge matrices in complete linearization laborious and computationally expensive simple Λ-iteration convergence problems way out > Accelerated lambda iteration (ALI) (Cannon 1973) formal solution J = Λ S + Λ S Λ S Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 22
41 Accelerated lambda iteration huge matrices in complete linearization laborious and computationally expensive simple Λ-iteration convergence problems way out > Accelerated lambda iteration (ALI) (Cannon 1973) formal solution J = Λ S + (Λ Λ ) S Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 22
42 Accelerated lambda iteration huge matrices in complete linearization laborious and computationally expensive simple Λ-iteration convergence problems way out > Accelerated lambda iteration (ALI) (Cannon 1973) formal solution iteration scheme J = Λ S + (Λ Λ ) S J (n+1) = Λ S (n+1) + (Λ Λ ) S (n) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 22
43 Accelerated lambda iteration huge matrices in complete linearization laborious and computationally expensive simple Λ-iteration convergence problems way out > Accelerated lambda iteration (ALI) (Cannon 1973) formal solution iteration scheme J = Λ S + (Λ Λ ) S J (n+1) = Λ S (n+1) + (Λ Λ ) S (n) }{{} J (n) correction term Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 22
44 Accelerated lambda iteration huge matrices in complete linearization laborious and computationally expensive simple Λ-iteration convergence problems way out > Accelerated lambda iteration (ALI) (Cannon 1973) formal solution iteration scheme J = Λ S + (Λ Λ ) S J (n+1) = Λ S (n+1) + (Λ Λ ) S (n) }{{} J (n) correction term compare: J (n+1) = Λ S (n) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 22
45 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 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 23
46 Accelerated lambda iteration iteration scheme J (n+1) = Λ S (n+1) + (Λ Λ ) S (n) = Λ S (n+1) + J (n) J (n+1) radiative rates n i R il = n i 4π n l R li = n l g i g j 4π αil () h αij () h [ Λ S (n+1) ( 2h 3 c 2 + ] + J (n) d [ Λ S (n+1) ] ) + J (n) d Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 23
47 Accelerated lambda iteration equations of statistical equilibrium l i {n l [R li (n i, n l ) + C li ] n i (R il (n i, n l ) + C il )} = 0 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 Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 24
48 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 solution by Newton-Raphson method (linearization) preconditioning Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 24
49 Summary two basic methods of the solution of the RTE+ESE system 1. complete linearization (Kiel) 2. accelerated lambda iteration (DETAIL) Non-LTE Line Formation for Trace Elements in Stellar Atmospheres, Nice p. 25
Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres
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 20.03.2008 p. Outline 1.
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationCMFGEN. Associate every atomic species with a given identification number
CMFGEN Associate every atomic species with a given identification number For each species, set: atomic mass atomic number solar abundance species abbreviations species identification Read in MODEL_SPEC
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 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 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 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 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 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 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 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 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 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 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 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 informationHale Collage. Spectropolarimetric Diagnostic Techniques!!!!!!!! Rebecca Centeno
Hale Collage Spectropolarimetric Diagnostic Techniques Rebecca Centeno March 1-8, 2016 Today Simple solutions to the RTE Milne-Eddington approximation LTE solution The general inversion problem Spectral
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 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 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 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 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 informationxii ASTROPHYSICAL BACKGROUND
Preface The purpose of this book is to provide an account of the major advances achieved in making quantitative spectroscopic analyses of the observable outer layers of stars. These analyses can now be
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 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 informationLecture 4: Absorption and emission lines
Lecture 4: Absorption and emission lines Senior Astrophysics 2018-03-13 Senior Astrophysics () Lecture 4: Absorption and emission lines 2018-03-13 1 / 35 Outline 1 Absorption and emission line spectra
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 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 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 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
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 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 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 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 informationPhoton Physics. Week 4 26/02/2013
Photon Physics Week 4 6//13 1 Classical atom-field interaction Lorentz oscillator: Classical electron oscillator with frequency ω and damping constant γ Eqn of motion: Final result: Classical atom-field
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 informationTo remain at 0 K heat absorbed by the medium must be removed in the amount of. dq dx = q(l) q(0) dx. Q = [1 2E 3 (τ L )] σ(t T 4 2 ).
294 RADIATIVE HEAT TRANSFER 13.7 Two infinite, isothermal plates at temperatures T 1 and T 2 are separated by a cold, gray medium of optical thickness τ L = L (no scattering). (a) Calculate the radiative
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 informationWINDS OF HOT MASSIVE STARS III Lecture: Quantitative spectroscopy of winds of hot massive stars
WINDS OF HOT MASSIVE STARS III Lecture: Quantitative spectroscopy of winds of hot massive stars 1 Brankica Šurlan 1 Astronomical Institute Ondřejov Selected Topics in Astrophysics Faculty of Mathematics
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 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 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 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 informationRadiation in climate models.
Lecture. Radiation in climate models. Objectives:. A hierarchy of the climate models.. Radiative and radiative-convective equilibrium.. Examples of simple energy balance models.. Radiation in the atmospheric
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 informationTheory of optically thin emission line spectroscopy
Theory of optically thin emission line spectroscopy 1 Important definitions In general the spectrum of a source consists of a continuum and several line components. Processes which give raise to the continuous
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 informationEmitted Spectrum Summary of emission processes Emissivities for emission lines: - Collisionally excited lines - Recombination cascades Emissivities
Emitted Spectrum Summary of emission processes Emissivities for emission lines: - Collisionally excited lines - Recombination cascades Emissivities for continuum processes - recombination - brehmsstrahlung
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 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 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 informationPlasma Spectroscopy Inferences from Line Emission
Plasma Spectroscopy Inferences from Line Emission Ø From line λ, can determine element, ionization state, and energy levels involved Ø From line shape, can determine bulk and thermal velocity and often
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 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 informationNearly everything that we know about the Sun (& other stars) comes from the photons it emits into space.
THE SOLAR ATMOSPHERE Nearly everything that we know about the Sun (& other stars) comes from the photons it emits into space. In the interior, we saw the photon distribution had a dominant isotropic part
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 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 informationLecture 2: Formation of a Stellar Spectrum
Abundances and Kinematics from High- Resolution Spectroscopic Surveys Lecture 2: Formation of a Stellar Spectrum Eline Tolstoy Kapteyn Astronomical Institute, University of Groningen I have a spectrum:
More informationSubstellar Atmospheres II. Dust, Clouds, Meteorology. PHY 688, Lecture 19 Mar 11, 2009
Substellar Atmospheres II. Dust, Clouds, Meteorology PHY 688, Lecture 19 Mar 11, 2009 Outline Review of previous lecture substellar atmospheres: opacity, LTE, chemical species, metallicity Dust, Clouds,
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 information25 years of progress in modeling stellar winds
Universität Kiel September 2001 Gaining Insight into Stellar Atmospheres International Workshop in the honor of the 80 th birthday of Prof. Dr. Kurt Hunger 25 years of progress in modeling stellar winds
More informationQuantum Electronics/Laser Physics Chapter 4 Line Shapes and Line Widths
Quantum Electronics/Laser Physics Chapter 4 Line Shapes and Line Widths 4.1 The Natural Line Shape 4.2 Collisional Broadening 4.3 Doppler Broadening 4.4 Einstein Treatment of Stimulated Processes Width
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 informationSpectroscopy. AST443, Lecture 14 Stanimir Metchev
Spectroscopy AST443, Lecture 14 Stanimir Metchev Administrative Homework 2: problem 5.4 extension: until Mon, Nov 2 Homework 3: problems 8.32, 8.41, 10.31, 11.32 of Bradt due in class Mon, Nov 9 Reading:
More informationAstronomy 421. Lecture 13: Stellar Atmospheres II. Skip Sec 9.4 and radiation pressure gradient part of 9.3
Astronomy 421 Lecture 13: Stellar Atmospheres II Skip Sec 9.4 and radiation pressure gradient part of 9.3 1 Announcements: Homework #4 is due Oct 3 Outline is due October 8 See example on the class web
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 informationRadiative Transfer with Polarization
The Radiative Transfer Equation with Polarization Han Uitenbroek National Solar Observatory/Sacramento Peak Sunspot, USA Hale COLLAGE, Boulder, Feb 16, 2016 Today s Lecture Equation of transfer with polarization
More informationOutline. Today we will learn what is thermal radiation
Thermal Radiation & Outline Today we will learn what is thermal radiation Laws Laws of of themodynamics themodynamics Radiative Radiative Diffusion Diffusion Equation Equation Thermal Thermal Equilibrium
More informationThe art of Stokes inversions. Luis R. Bellot Rubio Instituto de Astrofísica de Andalucía (CSIC), Spain National Astronomical Observatory, Japan
The art of Stokes inversions Instituto de Astrofísica de Andalucía (CSIC), Spain National Astronomical Observatory, Japan Outline What is an inversion technique? How do they work? ME inversions Accounting
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 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 informationPHYS 633: Introduction to Stellar Astrophysics
PHYS 633: Introduction to Stellar Astrophysics Spring Semester 2006 Rich Townsend (rhdt@bartol.udel.edu) Ionization During our derivation of the equation of state for stellar material, we saw that the
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 informationX-ray Emission from O Stars. David Cohen Swarthmore College
X-ray Emission from O Stars David Cohen Swarthmore College Young OB stars are very X-ray X bright L x up to ~10~ 34 ergs s -1 X-ray temperatures: few up to 10+ kev (10s to 100+ million K) K Orion; Chandra
More informationSOLAR SPECTRUM FORMATION
SOLAR SPECTRUM FORMATION Robert J. Rutten [Sterrekundig Instituut Utrecht & Institutt for Teoretisk Astrofysikk Oslo] Texts: R.J. Rutten, Radiative Transfer in Stellar Atmospheres (RTSA), my website D.
More informationOptical Depth & Radiative transfer
University of Naples Federico II, Academic Year 2011-2012 Istituzioni di Astrofisica, read by prof. Massimo Capaccioli Lecture 8 Optical Depth & Radiative transfer Learning outcomes The student will :
More informationStellar Astrophysics: Heat Transfer 1. Heat Transfer
Stellar Astrophysics: Heat Transfer 1 Heat Transfer Update date: September 3, 21 We discuss various energy transfer mechanisms in stars: radiative, conductive, and convective. 1 Radiative Transfer 1.1
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 informationExample: model a star using a two layer model: Radiation starts from the inner layer as blackbody radiation at temperature T in. T out.
Next, consider an optically thick source: Already shown that in the interior, radiation will be described by the Planck function. Radiation escaping from the source will be modified because the temperature
More informationIntroduction. Stellar Objects: Introduction 1. Why should we care about star astrophysics?
Stellar Objects: Introduction 1 Introduction Why should we care about star astrophysics? stars are a major constituent of the visible universe understanding how stars work is probably the earliest major
More informationHydrodynamic model atmospheres for WR stars. Self-consistent modeling of a WC star wind
A&A 432, 633 645 (2005) DOI: 10.1051/0004-6361:20041732 c ESO 2005 Astronomy & Astrophysics Hydrodynamic model atmospheres for WR stars Self-consistent modeling of a WC star wind G. Gräfener 1,2,3 and
More information