Modelling stellar atmospheres with full Zeeman treatment
|
|
- Preston Ford
- 6 years ago
- Views:
Transcription
1 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 workshop 12 th of September 2007
2 2 / 16 Contents 1 2 3
3 2 / 16 Contents 1 2 3
4 2 / 16 Contents 1 2 3
5 3 / 16 Outline of the problem Outline of the problem Bagnulo et al., 2001, A&A 369, 889 Properties of magnetic CP stars: upper main sequence stars T eff ranging from to K spectrum and photometric variability peculiar and stratified abundances magnetic fields, B ranging from 100 G to 35 kg, up to 100 kg at magnetic poles
6 3 / 16 Outline of the problem Outline of the problem Bagnulo et al., 2001, A&A 369, 889 Properties of magnetic CP stars: upper main sequence stars T eff ranging from to K spectrum and photometric variability peculiar and stratified abundances magnetic fields, B ranging from 100 G to 35 kg, up to 100 kg at magnetic poles
7 3 / 16 Outline of the problem Outline of the problem Bagnulo et al., 2001, A&A 369, 889 Properties of magnetic CP stars: upper main sequence stars T eff ranging from to K spectrum and photometric variability peculiar and stratified abundances magnetic fields, B ranging from 100 G to 35 kg, up to 100 kg at magnetic poles
8 3 / 16 Outline of the problem Outline of the problem Bagnulo et al., 2001, A&A 369, 889 Properties of magnetic CP stars: upper main sequence stars T eff ranging from to K spectrum and photometric variability peculiar and stratified abundances magnetic fields, B ranging from 100 G to 35 kg, up to 100 kg at magnetic poles
9 4 / 16 Outline of the problem Outline of the problem requirements specification for the atmospheric model code temperature range (of interest): to K peculiar and stratified abundances arbitrary inclination of the field in plane parallel models full Zeeman treatment, polarised Feautrier solver hydrostatic equilibrium with magnetic pressure VALD line data, CoCoS simplifications: plane parallel, Kurucz continuum routines, no dynamic phenomena considered, no microturbulence.
10 5 / 16 Software engineering Polarised radiation transfer Temperature correction Program features: ATLAS12 continua for comparability with standard models consistent with spectral synthesis (COSSAM) and radiative diffusion (CARAT) code written in Ada95 thread parallel modularised, object oriented
11 5 / 16 Software engineering Polarised radiation transfer Temperature correction Program features: ATLAS12 continua for comparability with standard models consistent with spectral synthesis (COSSAM) and radiative diffusion (CARAT) code written in Ada95 thread parallel modularised, object oriented
12 5 / 16 Software engineering Polarised radiation transfer Temperature correction Program features: ATLAS12 continua for comparability with standard models consistent with spectral synthesis (COSSAM) and radiative diffusion (CARAT) code written in Ada95 thread parallel modularised, object oriented
13 5 / 16 Software engineering Polarised radiation transfer Temperature correction Program features: ATLAS12 continua for comparability with standard models consistent with spectral synthesis (COSSAM) and radiative diffusion (CARAT) code written in Ada95 thread parallel modularised, object oriented
14 6 / 16 Software engineering Polarised radiation transfer Temperature correction Polarised radiation transfer equation d dz I = K I + K (S, 0, 0, 0) Stokes vector I = (I, Q, U, V ) absorption matrix K = κ c 1 + κ o Φ line absorption matrix Φ = Show details φ I φ Q φ U φ V φ Q φ I φ V φ U φ U φ V φ I φ Q φ V φ U φ Q φ I z, I B γ χ x y
15 7 / 16 Zeeman Feautrier solver Software engineering Polarised radiation transfer Temperature correction Feautrier equation can be generalised to the magnetic case in the presence of blends (see Alecian and Stift, 2004, A&A 416, 703) d J dτ 5000 = XH dh and dτ 5000 = X( J S) with X := K κ 5000 µ d dτ 5000 (X 1 d J dτ 5000 ) = X( J S) (inner points) dx 1 dτ 5000 d J dτ X 1 d 2 J dτ = X( J S) (boundary condition) N equations, N unknowns B 1 J1 C 1 J2 = L 1 A n Jn 1 + B n Jn C n Jn+1 = L n A N JN 1 + B N JN = L N
16 Temperature correction Software engineering Polarised radiation transfer Temperature correction Dreizler s Lucy Unsöld scheme adapted for polarised radiation transport equations: (see Dreizler, 2003, ASPC 288, 69) momenta of the polarised radiation transport equation differences ( X = X Equilibrium X Model ) of flux, intensity and radiation pressure two flux criteria local balance of emitted versus absorbed energy d(h I ) z dτ 5000 = 0 (constant flux) cannot be used in deep layers where the atmosphere becomes diffusive and S J local Planck function. nonlocal condition of constant flux: 0 Show momenta HI dωdν σ 4π T 4 eff = 0 (desired value of the flux) inefficient in regions with small opacities 8 / 16
17 9 / 16 Combining the flux criteria Software engineering Polarised radiation transfer Temperature correction T = π 4σT 3 R S d [Kν «J 0 ν] I dν 1 R [Kν] 0 I,IS S νdν {z } local energy conservation S R + d [Kν J 0 ν] I dν 2 R J I [Kν] 0 I,IS νdν f 0 H I (0) fg {z } surface flux S R + d [Kν J 0 ν] I dν 3 R J I [Kν] 0 I,IS νdν Z τ 0 1 f R 0 [Kν H ν] I dν H I κ 5000 H I dτ 5000 {z } global energy conservation «
18 10 / 16 Results and comparisons Comparison with Carpenter s results Comparisons with LLMODELS Summary Two effects enhanced line blanketing magnetic pressure Comparison with the results of Carpenter (enhanced line blanketing and magnetic pressure see Carpenter, 1985, ApJ 289, 660, Carpenter, 1983, PhD Thesis ) LLMODELS (enhanced line blanketing, see Kochukhov et al., 2005, A&A 433, 671, Khan and Shulyak, 2006a, A&A 448, 1153, Khan and Shulyak, 2006b, A&A 454, 933)
19 11 / 16 Comparison with Carpenter s results Comparisons with LLMODELS Summary Comparison with Carpenter s results CAMAS results Carpenter, 1985, ApJ 289, 660
20 12 / 16 Comparisons with LLMODELS Comparison with Carpenter s results Comparisons with LLMODELS Summary differences between the field less and the isotropic model CAMAS results Kochukhov et al., 2005, A&A 433, 671
21 12 / 16 Comparisons with LLMODELS Comparison with Carpenter s results Comparisons with LLMODELS Summary differences between the field less and the isotropic model CAMAS results Kochukhov et al., 2005, A&A 433, 671
22 13 / 16 Depth scales Comparison with Carpenter s results Comparisons with LLMODELS Summary τ ross is affected by (magnetic) line blanketing, τ 5000 is not directly affected
23 14 / 16 Depth scales Comparison with Carpenter s results Comparisons with LLMODELS Summary τ ross is affected by (magnetic) line blanketing, τ 5000 is not directly affected
24 15 / 16 Comparison with LLMODELS Comparison with Carpenter s results Comparisons with LLMODELS Summary differences between the anisotropic and the isotropic model CAMAS results
25 15 / 16 Comparison with LLMODELS Comparison with Carpenter s results Comparisons with LLMODELS Summary differences between the anisotropic and the isotropic model Khan and Shulyak, 2006b, A&A 454, 933
26 Comparison with LLMODELS Comparison with Carpenter s results Comparisons with LLMODELS Summary differences between the anisotropic and the isotropic model CAMAS results Khan and Shulyak, 2006b, A&A 454, 933 Zoom 15 / 16
27 16 / 16 Summary Comparison with Carpenter s results Comparisons with LLMODELS Summary The Zeeman effect enhances the opacity. This affects the optical depth and the temperature structure. CAMAS essentially confirms the results of LLMODELS. The isotropic models are largely identical, however the anisotropic models with strong magnetic fields exhibit notable differences that need to be clarified. The atmospheres computed with CAMAS will be used (among others) for the modelling of diffusion processes. Acknowledgements: this work was supported by the Austrian Science Funds (FWF project P16003 of M.J.Stift). KMB wants to thank Bob Kurucz for the useful discussions on the ATLAS code, Christian Stütz and the VALD team for the line data and Holger Pikall for tree weeks of CPU-time.
28 17 / 16 Bibliography I Appendix References Alecian, G. and Stift, M. J.: 2004, A&A 416, 703 Bagnulo, S., Wade, G. A., Donati, J.-F., Landstreet, J. D., Leone, F., Monin, D. N., and Stift, M. J.: 2001, A&A 369, 889 Carpenter, K. G.: 1983, Ph.D. thesis, AA(Ohio State Univ., Columbus.) Carpenter, K. G.: 1985, ApJ 289, 660 Dreizler, S.: 2003, in I. Hubeny, D. Mihalas, and K. Werner (eds.), ASP Conf. Ser. 288: Stellar Atmosphere Modeling, Vol. 288 of ASP Conference Series, p. 69, Review of temperature correction schemes, Damping Khan, S. A. and Shulyak, D. V.: 2006a, A&A 448, 1153 Khan, S. A. and Shulyak, D. V.: 2006b, A&A 454, 933 Kochukhov, O., Khan, S., and Shulyak, D.: 2005, A&A 433, 671
29 Appendix References Comparison with LLMODELS differences between the anisotropic and the isotropic model CAMAS results Khan and Shulyak, 2006b, A&A 454, 933 Return 18 / 16
30 19 / 16 first depth point Appendix References d J dτ = X 1 ( J 1 (1 e τ 1X 1 ) S 1 ) J1 = d J 2 τ J (2,1) dτ 1 ( τ (2,1)) 2 d 2 J 2 dτ ( forward Taylor series) A 1 = 0 B 1 = 1 + τ (2,1) X 1 ( τ (2,1)) 2 dx 1 X 1 1 X 1 + ( τ (2,1)) 2 X dτ C 1 = 1 L1 = ( τ (2,1) ( τ (2,1)) 2 dx 1 X 1 1 )X 1 (1 e τ 1X 1 ) S 1 + ( τ (2,1)) 2 X 2 1S 1 2 dτ
31 20 / 16 inner depth points Appendix References d d J dτ X 1 = dτ n X 1 Jn+1 J n n+ 1 τ X 1 Jn J n 1 (n+1,n) n 1 τ (n,n 1) 2 2 τ (n+1,n) + τ (n,n 1) 2 A n = B n = C n = (X n 1 + X n ) 1 τ (n,n 1) τ (n+1,n 1) (X n 1 X n ) 1 (X 1 n X n) 1 + X n τ (n,n 1) τ (n+1,n 1) τ (n+1,n) τ (n+1,n 1) (X n+1 + X n ) 1 τ (n+1,n) τ (n+1,n 1) Ln = X n Sn
32 Flowchart of CAMAS Appendix References read input, initialise enter radiative equilibrium loop check convection recompute equidistant τ 5000 and interpolate model data EXIT if convergence criteria are met pretabulate continua if necessary reselect lines if necessary enter integration loop (> 90% of computation time) independent tasks use a magnetic Feautrier solver for all wavelength points and add the results calculate temperature corrections calculate and apply pressure corrections search for noise in the flux distribution and check convergence apply temperature corrections if not converged yet apply Ng acceleration (optional) 21 / 16
33 22 / 16 Feautrier solver Appendix References solve radiation transport equation as boundary value problem di ± dτ 5000 = I ± S define flux like and intensity like quantities combine equations for outward and inward rays discretize and add boundary conditions
34 22 / 16 Feautrier solver Appendix References solve radiation transport equation as boundary value problem define flux like and intensity like quantities H n = I+ n I n 2, J n = I+ n + I n 2 combine equations for outward and inward rays discretize and add boundary conditions
35 22 / 16 Feautrier solver Appendix References solve radiation transport equation as boundary value problem di ± dτ 5000 = I ± S define flux like and intensity like quantities combine equations for outward and inward rays dj dh = H, = J S dτ 5000 dτ }{{ 5000 } d 2 J dτ = J S discretize and add boundary conditions
36 22 / 16 Feautrier solver Appendix References solve radiation transport equation as boundary value problem define flux like and intensity like quantities combine equations for outward and inward rays discretize and add boundary conditions I (τ = 0) = 0, no incident radiation on the surface dj dτ 5000 τ=0 = J(τ = 0) H N = I + N J N, diffusive at innermost depth-point dj N = I + N dτ J N with I + N = S N + ds N 5000 dτ 5000
37 Appendix References Angle dependent opacity line absorption terms φ I = φ Q = φ U = φ V = φ p, b, r Faraday terms φ Q = φ U = φ V = φ p, b, r Return 1 4 (2 φ p sin 2 γ + (φ r + φ b )(1 + cos 2 γ)) 1 4 (2 φ p (φ r + φ b )) sin 2 γ cos 2χ 1 4 (2 φ p (φ r + φ b )) sin 2 γ sin 2χ 1 2 (φ r φ b ) cos γ... line absorption profiles 1 4 (2 φ p (φ r + φ b )) sin2 γ cos 2χ 1 4 (2 φ p (φ r + φ b )) sin2 γ sin 2χ 1 2 (φ r φ b ) cos γ... anomalous dispersion profiles z B I β θ x ψ y z, I B γ x χ y 23 / 16
38 24 / 16 Appendix References Combining the flux criteria Return 0th moment: d H I 4σT 3 π dh eq dτ R 5000 [K 0 νxν,eq] I dν X eq S T = 4σT 3 π R 0 R 0 dτ 5000 = [Kν J ν] I dν κ 5000 = 0, 0th moment for dhmod dτ R 5000 [K 0 νxν,mod] I dν X model T = S [Kν J ν] I dν R 0 [Kν] I,IS S + S νdν R d K I 0 dτ 5000 = [Kν H ν] I dν κ 5000 integration variable Eddington factors f τ = Kτ 1st moment: f J I }{{} K I = f0 g H I (0) + τ 0 R 0 R 0 R 0 [Kν] I,IS νdν [Kν] I,I S νdν κ 5000 R 0 J τ and g = H 0 [Kν H ν] I dν κ 5000 dτ 5000 [Kν J ν] I dν J I J I J 0
Observing and modelling stellar magnetic fields. 2. Models
Observing and modelling stellar magnetic fields. 2. Models John D Landstreet Department of Physics & Astronomy University of Western Ontario London, Canada West In the previous episode... We explored the
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 informationAtmospheres of CP stars: magnetic field effects
Mem. S.A.It. Vol. 7, 99 c SAIt 2005 Memorie della Supplementi Atmospheres of CP stars: magnetic field effects D. Shulyak 1,2, G. Valyavin 3,4, O. Kochukhov 5, S. Khan 6, and V. Tsymbal 1,2 1 Tavrian National
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 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 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 informationRadiative 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 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 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 informationTHE OBSERVATION AND ANALYSIS OF STELLAR PHOTOSPHERES
THE OBSERVATION AND ANALYSIS OF STELLAR PHOTOSPHERES DAVID F. GRAY University of Western Ontario, London, Ontario, Canada CAMBRIDGE UNIVERSITY PRESS Contents Preface to the first edition Preface to the
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 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 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 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 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 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 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 informationAbundance structure of the atmospheres of magnetic CP stars
Abundance structure of the atmospheres of magnetic CP stars T.A. Ryabchikova Institute of Astronomy, Russian Academy of Sciences Moscow, Russia Institute for Astronomy, Vienna University Vienna, Austria
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 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 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 3: Specific Intensity, Flux and Optical Depth
Lecture 3: Specific Intensity, Flux and Optical Depth We begin a more detailed look at stellar atmospheres by defining the fundamental variable, which is called the Specific Intensity. It may be specified
More informationLine-by-line opacity stellar model atmospheres
A&A 428, 993 1000 (2004) DOI: 10.1051/0004-6361:20034169 c ESO 2004 Astronomy & Astrophysics Line-by-line opacity stellar model atmospheres D. Shulyak 1,2, V. Tsymbal 1,2, T. Ryabchikova 2,3, Ch. Stütz
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 informationShort/Simple Definitions:
Eric Joseph Bubar Stellar Atmosphere/Interiors Portfolio CHAPTER : CURVES OF GROWTH Short/Simple Definitions: Curve of Growth: Plot of equivalent widths versus number of absorbing atoms that created that
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 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 informationSolution 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 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 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 informationarxiv:astro-ph/ v1 25 Sep 2006
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Are Wolf-Rayet winds driven by radiation? arxiv:astro-ph/0609675v1 25 Sep 2006 Götz Gräfener & Wolf-Rainer
More informationMagnetic Doppler imaging considering atmospheric structure modifications due to local abundances: a luxury or a necessity?
Mon. Not. R. Astron. Soc. 421, 3004 3018 (2012) doi:10.1111/j.1365-2966.2012.20526.x Magnetic Doppler imaging considering atmospheric structure modifications due to local abundances: a luxury or a necessity?
More information1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER
1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER 5.1 Introduction As radiation struggles to make its way upwards through a stellar atmosphere, it may be weakened by absorption
More information1. Runaway Greenhouse
1. Runaway Greenhouse Making the simple complicated is commonplace; making the complicated simple, awesomely simple, that s creativity. - Charlie Mingus (Jazz Musician) Initially we assume that the radiative
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 informationarxiv: v1 [astro-ph] 5 Dec 2007
Contrib. Astron. Obs. Skalnaté Pleso 1, 1 6, (2007) Magnetic, Chemical and Rotational Properties of the Herbig Ae/Be Binary System HD 72106 arxiv:0712.0771v1 [astro-ph] 5 Dec 2007 C.P. Folsom 1,2,3, G.A.Wade
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 informationSome HI is in reasonably well defined clouds. Motions inside the cloud, and motion of the cloud will broaden and shift the observed lines!
Some HI is in reasonably well defined clouds. Motions inside the cloud, and motion of the cloud will broaden and shift the observed lines Idealized 21cm spectra Example observed 21cm spectra HI densities
More informationarxiv: v1 [astro-ph.sr] 16 Nov 2016
Astronomy & Astrophysics manuscript no. 29768 c ESO 2018 July 19, 2018 Doppler imaging of chemical spots on magnetic Ap/Bp stars Numerical tests and assessment of systematic errors O. Kochukhov Department
More informationThe effect of the surface distribution of elements on measuring the magnetic field of chemically peculiar stars. The case of the roap star HD 24712
A&A 425, 271 280 (2004) DOI: 10.1051/0004-6361:20047180 c ESO 2004 Astronomy & Astrophysics The effect of the surface distribution of elements on measuring the magnetic field of chemically peculiar stars
More informationDiffuse Interstellar Medium
Diffuse Interstellar Medium Basics, velocity widths H I 21-cm radiation (emission) Interstellar absorption lines Radiative transfer Resolved Lines, column densities Unresolved lines, curve of growth Abundances,
More informationStellar Atmospheres: A Computational Approach. Senior Project by Kathy Geise
Stellar Atmospheres: A Computational Approach Senior Project by Kathy Geise University of Denver March 2008 Abstract The stellar atmosphere contains the outermost layers of the star and includes the region
More informationt KH = GM2 RL Pressure Supported Core for a Massive Star Consider a dense core supported by pressure. This core must satisfy the equation:
1 The Kelvin-Helmholtz Time The Kelvin-Helmhotz time, or t KH, is simply the cooling time for a pressure supported (i.e. in hydrostatic equilibrium), optically thick object. In other words, a pre-main
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 informationNew Results of Spectral Observations of CP Stars in the Li I 6708 Å Spectral Region with the 6 m BTA Telescope
Magnetic Stars, 2011, pp. 374 381 New Results of Spectral Observations of CP Stars in the Li I 6708 Å Spectral Region with the 6 m BTA Telescope Polosukhina N. 1, Shavrina A. 2, Drake N. 3, Kudryavtsev
More informationSurface magnetic fields across the HR Diagram
Surface magnetic fields across the HR Diagram John Landstreet University of Western Ontario London, Upper Canada & Armagh Observatory Armagh, Northern Ireland Objective of this talk... This meeting combines
More informationRadiative transfer in a floating stratus irradiated by a luminous source I. Basics and a plane flat source
MNRAS 436, 2550 2559 (2013) Advance Access publication 2013 October 22 doi:10.1093/mnras/stt1751 Radiative transfer in a floating stratus irradiated by a luminous source I. Basics a plane flat source J.
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 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 informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
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 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 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 informationThe Sun. Basic Properties. Radius: Mass: Luminosity: Effective Temperature:
The Sun Basic Properties Radius: Mass: 5 R Sun = 6.96 km 9 R M Sun 5 30 = 1.99 kg 3.33 M ρ Sun = 1.41g cm 3 Luminosity: L Sun = 3.86 26 W Effective Temperature: L Sun 2 4 = 4πRSunσTe Te 5770 K The Sun
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 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 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 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 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 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 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 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 informationASTROPHYSICS. K D Abhyankar. Universities Press S T A R S A ND G A L A X I E S
ASTROPHYSICS S T A R S A ND G A L A X I E S K D Abhyankar Universities Press Contents Foreword vii Preface ix 1 Introduction 1 1.1 ' Astronomy and astrophysics 1 1.2 Importance of astronomy 2 1.3 Methods
More informationAstr 323: Extragalactic Astronomy and Cosmology. Spring Quarter 2014, University of Washington, Željko Ivezić. Lecture 1:
Astr 323: Extragalactic Astronomy and Cosmology Spring Quarter 2014, University of Washington, Željko Ivezić Lecture 1: Review of Stellar Astrophysics 1 Understanding Galaxy Properties and Cosmology The
More informationINFRAMET. 2.1 Basic laws
tel: 048 60844873, fax 48 6668780. Basic laws.. Planck law All objects above the temperature of absolute zero emit thermal radiation due to thermal motion of the atoms and the molecules. The hotter they
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 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 informationAtmospheric Sciences 321. Science of Climate. Lecture 6: Radiation Transfer
Atmospheric Sciences 321 Science of Climate Lecture 6: Radiation Transfer Community Business Check the assignments Moving on to Chapter 3 of book HW #2 due next Wednesday Brief quiz at the end of class
More informationThree-dimensional magnetic and abundance mapping of the cool Ap star HD 24712
A&A 73, A123 (21) DOI: 1.11/4-6361/214249 c ESO 21 Astronomy & Astrophysics Three-dimensional magnetic and abundance mapping of the cool Ap star HD 24712 II. Two-dimensional magnetic Doppler imaging in
More informationThe Formation of Spectral Lines. I. Line Absorption Coefficient II. Line Transfer Equation
The Formation of Spectral Lines I. Line Absorption Coefficient II. Line Transfer Equation Line Absorption Coefficient Main processes 1. Natural Atomic Absorption 2. Pressure Broadening 3. Thermal Doppler
More informationDescription of radiation field
Description of radiation field Qualitatively, we know that characterization should involve energy/time frequency all functions of x,t. direction We also now that radiation is not altered by passing through
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 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 informationASTR240: Radio Astronomy
AST24: adio Astronomy HW#1 Due Feb 6, 213 Problem 1 (6 points) (Adapted from Kraus Ch 8) A radio source has flux densities of S 1 12.1 Jy and S 2 8.3 Jy at frequencies of ν 1 6 MHz and ν 2 1415 MHz, respectively.
More informationarxiv:astro-ph/ v1 5 May 2004
Modelling of Stellar Atmospheres IAU Symposium, Vol. xxx, xxxx N. E. Piskunov, W. W. Weiss, D. F. Gray, eds. New Grids of ATLAS9 Model Atmospheres arxiv:astro-ph/0405087v1 5 May 2004 Fiorella Castelli
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 information3D atmospheric structure of the prototypical roap star HD (HR1217)
Contrib. Astron. Obs. Skalnaté Pleso 38, 335 340, (2008) 3D atmospheric structure of the prototypical roap star HD 24712 (HR1217) T. Lüftinger 1, O. Kochukhov 2, T. Ryabchikova 1,3, N. Piskunov 2, W.W.
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 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 informationAstr 323: Extragalactic Astronomy and Cosmology. Spring Quarter 2012, University of Washington, Željko Ivezić. Lecture 1:
Astr 323: Extragalactic Astronomy and Cosmology Spring Quarter 2012, University of Washington, Željko Ivezić Lecture 1: Review of Stellar Astrophysics 1 Understanding Galaxy Properties and Cosmology The
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 informationThe nature of the light variations of chemically peculiar stars
The nature of the light variations of chemically peculiar stars Jiří Krtička, Zdeněk Mikulášek Masaryk University, Brno, Czech Republic Juraj Zverko, Jozef Žižňovský Astronomical Institute SAV, Tatranská
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 informationINTRODUCTION TO MICROWAVE REMOTE SENSING - II. Dr. A. Bhattacharya
1 INTRODUCTION TO MICROWAVE REMOTE SENSING - II Dr. A. Bhattacharya The Radiation Framework The information about features on the Earth s surface using RS depends on measuring energy emanating from the
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature10389 Observations The Kepler satellite gathered 30 days of uninterrupted data for HD 187547 (KIC 7548479) between October 19 and November 18, 2009, corresponding to Kepler s observing
More informationAstronomy. Astrophysics. Magnetic field topology and chemical spot distributions in the extreme Ap star HD 75049
DOI: 10.1051/0004-6361/201425065 c ESO 2015 Astronomy & Astrophysics Magnetic field topology and chemical spot distributions in the extreme Ap star HD 75049 O. Kochukhov 1, N. Rusomarov 1,J.A.Valenti 2,H.C.Stempels
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 Electromagnetic Spectrum
Astr 102: Introduction to Astronomy Fall Quarter 2009, University of Washington, Željko Ivezić Lecture 4: The Electromagnetic Spectrum 1 Understanding Stellar and Galaxy Properties, and Cosmology Four
More informationThe null result of a search for pulsational variations of the surface magnetic field in the roap star γ Equulei
Mon. Not. R. Astron. Soc. 351, L34 L38 (2004) doi:10.1111/j.1365-2966.2004.07946.x The null result of a search for pulsational variations of the surface magnetic field in the roap star γ Equulei O. Kochukhov,
More informationStefan-Boltzmann law for the Earth as a black body (or perfect radiator) gives:
2. Derivation of IPCC expression ΔF = 5.35 ln (C/C 0 ) 2.1 Derivation One The assumptions we will make allow us to represent the real atmosphere. This remarkably reasonable representation of the real atmosphere
More informationRadiation in the atmosphere
Radiation in the atmosphere Flux and intensity Blackbody radiation in a nutshell Solar constant Interaction of radiation with matter Absorption of solar radiation Scattering Radiative transfer Irradiance
More informationLecture 2 Line Radiative Transfer for the ISM
Lecture 2 Line Radiative Transfer for the ISM Absorption lines in the optical & UV Equation of transfer Absorption & emission coefficients Line broadening Equivalent width and curve of growth Observations
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 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 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 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 informationAstronomy. Astrophysics. Isotopic anomaly and stratification of Ca in magnetic Ap stars. T. Ryabchikova 1,2, O. Kochukhov 3, and S.
A&A 480, 811 823 (2008) DOI: 10.1051/0004-6361:20077834 c ESO 2008 Astronomy & Astrophysics Isotopic anomaly and stratification of Ca in magnetic Ap stars T. Ryabchikova 1,2, O. Kochukhov 3, and S. Bagnulo
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 informationarxiv: v1 [astro-ph.sr] 24 Apr 2013
Astronomy & Astrophysics manuscript no. 21467 c ESO 2013 April 26, 2013 Are there tangled magnetic fields on HgMn stars? O. Kochukhov 1, V. Makaganiuk 1, N. Piskunov 1, S. V. Jeffers 2, C. M. Johns-Krull
More informationOpacity. requirement (aim): radiative equilibrium: near surface: Opacity
(Gray) Diffusion approximation to radiative transport: (assumes isotropy valid only in the deep stellar interior) - opacity is a function of frequency (wave length ). - aim: to reduce the (rather complex)
More information