RADIATIVE TRANSFER IN AXIAL SYMMETRY
|
|
- Hilda Sharp
- 5 years ago
- Views:
Transcription
1 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 a method for solving the radiative transfer equation in axial symmetry both in static case and in moving media. A usage of this method is shown here for stellar atmospheres, stellar winds, and accretion discs. The effects of limb darkening and stellar rotation are discussed, too. Introduction The astrophysical radiative transfer problem is a highly nonlinear task, which has not been solved in its complexity until now. The analytical solution of this problem is impossible, but also the numerical calculations are very time consuming. Due to this reason we have to use simplifications, which are based on physical conditions of a given object. One of the key parts of this problem is the solution of the radiative transfer equation (RTE). For optically thick regions we can use the diffusion approximation method (Kneer & Heasley, 979) for the solution of the radiative transfer equation. In an optically thin region it is sometimes possible only to sum the radiation from various regions neglecting absorption, or we can use the Monte Carlo method (Boissé, 99), which is very fast in these situations. An additional problem is the velocity field. With the exception of supernovae and accretion discs around black holes we can ignore relativistic effects on the radiation field. But even if the problem is nonrelativistic, in most astrophysical situations its solution is not trivial, since the opacity and emissivity can change very much along a ray due to the Doppler shift. If the velocity gradient is sufficiently high, we can use the Sobolev approximation (Sobolev, 947). If the Sobolev condition is not fulfilled, we have to calculate with the additional term on the left hand side of the radiative transfer equation. For the solution of this non-sobolev problem we can use either the observer frame or the comoving frame, the latter being more general. Astronomický ústav, Akademie věd České republiky, CZ Ondřejov, Czech Republic c EDP Sciences 26 DOI: (will be inserted later)
2 2 Title : will be set by the publisher We want to describe the main features of the stellar wind and stellar rotation. This allows us to use the approximation of axial symmetry, which is appropriate also for the disc geometry, as is shown later. The more detailed description of this method is published in Korčáková & Kubát (25a). Here we show some new applications of this method. 2 Description of the method The basic simplifying idea of our method is in the solution of the radiative transfer equation not in the whole star, but separately in planes intersecting the star. In every plane, the radiative transfer equation is solved using a combination of short and long characteristic method. The main advantage of the combination of characteristics is in better description of the radiation field than the short characteristic method and does not consume as much computing time as the long characteristic method. We are not going into details in this paper, since an extensive description as well as the comparison of this method with another independent models is published in Korčáková & Kubát (25a). We section a star into longitudinal planes (see Fig., left panel). The radiative θ radial grid lines rotation axis r φ concentric circles R Fig.. left panel: A stellar sectioning and the coordinate system of longitudinal planes. right panel: Obtaining the whole radiation field by rotating the longitudinal planes (which is possible due to the axial symmetry). transfer equation is solved in every plane independently. The whole radiation field is obtained by rotation of these planes around the rotation axis of the star (see Fig., right panel). In every longitudinal plane we introduce the polar coordinate system. The solution in the given plane must start from the upper (i.e. outer)
3 Korčáková & Kubát: Radiative Transfer in Axial Symmetry 3 boundary. In every grid point we choose several (up to 9) rays per quadrant (see Fig. 2, left panel). Along these rays we solve the radiative transfer equation τ(ab) I (B) = I (A) e τ (AB) + S(t)e [ ( τ(ab) t)] dt. (2.) The quantities used have their usual meaning. The interval AB is a section of the ray within each cell, in which the source function S is assumed to be a linear function of the optical depth. radial grid line α grid circle zone Fig. 2. The schema for the solution of the radiative transfer equation in longitudinal planes. The rays (characteristics) are denoted using the dashed line. left panel: downward (inward) solution, right panel: upward (outward) solution. If the intensity in the central grid circle is known, we must check whether the longitudinal plane intersects the lower (inner) boundary. If not, we must solve the radiative transfer also in the central circle of this plane. In this way we obtain complete downward radiation intensity. The upward solution is very similar to the downward one. In every grid circle we choose the rays under the same angles as before, which is depicted in the Fig. 2 (right panel). We can see from this figure why it is necessary to start from the upper boundary. Since some rays can start and end on the same circle, we need to know the intensity at the beginning of the characteristic (which is the downward intensity) to solve the equation 2. upwards. If velocity is present, we can assume that the velocity vector is constant in the given cell and that it changes only at the cell boundaries. In this assumption we can use the Lorentz invariance of the radiative transfer equation to solve the static equation in given cells. The transformation of frequency is performed only at the cell boundaries (it is possible to neglect the transformation of intensity in case of nonrelativistic velocities).
4 4 Title : will be set by the publisher 3 Selected results In this section we show some calculations based on our method, which were presented at this workshop. Since most of these results have already been published, we show here the main ideas only and refer to corresponding papers for details. 3. limb darkening Limb darkening can be used as one of the most suitable tests of our stellar atmosphere models, since it depends on temperature and density distribution in stellar atmosphere. Its exact knowledge affects the quality of information that can be obtained from interferometric observations, as well as determination of rotation velocities in stars (see next section 3.2). Limb darkening is frequency dependent, which is shown in the Fig. 3. Left limb darkening specific intensity limb darkening in line e e e x [R ] e+4 relative intensity x [R ] specific intensity 4e-5 5e-5 4e-5 continuum line center 3e-5 2e-5 e x [r/r ] e e e e4 specific intensity 3e-5 2e-5 e x [r/r ] Fig. 3. Limb darkening for a thin stellar atmosphere (top) and for an extended stellar atmosphere (bottom). The x-axis is the distance from the center of the star in units of stellar radius (x = for the stellar radius). left panels: 3D plot of limb darkening. right panels: Intensity variation for the central line frequency. panels show the 3D plots of limb darkening. Right panels emphasize in detail an intensity dependence across the stellar disk for the central frequency of a line,
5 Korčáková & Kubát: Radiative Transfer in Axial Symmetry 5 where limb brightening instead of limb darkening is present. This is similar to the effect of flash spectra in solar chromosphere. This effect is more important for extended stellar atmospheres. The limb darkening (brightening) is important for the evaluation of rotation velocities of these stars, since the easiest method, convolution, can produce a huge error. The results in Fig. 3 are calculated for a model atmosphere of a hot B star with effective temperature T eff = 7 3 K, gravitation acceleration log g = 4.2, and radius 3.26R. The stellar radius is artificially distended for the case of the extended atmosphere. 3.2 stellar rotation Stellar rotation is one of the most important features in stars. Unfortunately, until now a lot of unresolved questions remain. In stellar spectra modelling one uses very often rather crude simplifications. The most commonly used one is the convolution of the static profile (H()) with a rotating profile (G()), F = H() G() = H( ) G( ) d. (3.) F c The rotation profile depends on limb darkening. Since limb darkening is strongly frequency dependent, the obtained results can suffer from a huge systematic error. In Fig. 4 we show the comparison of a profile obtained from our axially symmetrical code with the profile obtained using convolution 3.. For the latter case, the limb darkening is taken in the form (see Gray, 976) I(x) = ( ɛ) + ɛ( x 2 ) /2. (3.2) Here x means the distance from the center of star in units of stellar radius (x = for r = R, where R is the stellar radius). The calculations are presented for a cool B-type main sequence star. For more details about the chosen model and results see Korčáková & Kubát (25b). Since our geometrical approximation is axial symmetry, we can naturally include gravity darkening and differential rotation. Obtained line profiles for the case of thin stellar atmosphere (cool main-sequence B type star) and extended atmosphere are plotted in Fig. 5. A detailed description of the stellar atmosphere models is published in Korčáková & Kubát (25b). 3.3 stellar wind There is a longstanding question about the origin of the stellar wind of hot Be stars. To answer this question, we need a model, which is able to include both deep stellar atmosphere layers and the wind region. The presented method is able to solve the radiative transfer equation in a stellar atmosphere where the gradient of global velocity is very small (and the Sobolev approximation is not valid), as well as in the stellar wind region, where the velocity gradient can be large.
6 6 Title : will be set by the publisher.98 relative flux ε = ε = "rot-grav-dif" "rot+grav+dif" 4.56e e e e+4 Fig. 4. The comparison of line profiles of cool B-type main sequence star calculated using convolution (ɛ = and ɛ =, see Eq. 3.2) with two line profiles obtained from our code. One presents the rigid body solution and the other one the case including gravity darkening and differential rotation relative flux rot-grav-dif rot-grav-dif-velocity rot+grav-dif rot-grav+dif rot+grav+dif rot+grav+dif-velocity relative flux rot-grav-dif rot-grav-dif-velocity rot+grav-dif rot-grav+dif rot+grav+dif rot+grav+dif-velocity e e e e e e e e+4 Fig. 5. The Hα line profiles for a thin stellar atmosphere (left panel) and for an extended stellar atmosphere (right panel). The line profiles are calculated with gravity darkening (+grav), differential rotation (+dif ) included, or excluded ( dif, grav). The lines indicated by velocity are obtained by neglecting the velocity field in the solution of the radiative transfer equation. The velocity field is present only in the flux calculation in this case. For example, in the Fig. 6 we plotted the Hα line profiles for the beta wind velocity law (see, e.g., Lamers & Cassinelli, 999) v(r) = v { [ ( vr v ) ] β } β R (3.3) r for various β parameters (β =.5,, 2). The input stellar parameters are the same as for limb darkening (section 3.) with the velocity in the photosphere v R = 2 km s and a terminal velocity of v = 2 km s. Since our method is based on the Local Lorentz Transformation (see Korčáková & Kubát
7 Korčáková & Kubát: Radiative Transfer in Axial Symmetry 7 23) we can calculate the decelerating velocity field, which is also plotted in Fig. 6. More details about this model and solution are in Korčáková & Kubát (25a)..95 relative flux β =.5 β = β = 2 deceleration 4.565e e e+4 Fig. 6. Profile of the Hα line for the case of stellar wind for three values of the parameter β =.5,, 2 (see Eq. 3.3). For comparison, the line profile affected by decelerating velocity field is also plotted (fine dotted line). Though we show results for the stellar wind of hot stars here, this method is applicable to cool stars (Korčáková et al., 24) as well. 3.4 accretion disc Since we have a flexibility in selection of the latitudinal grid (angle θ), it is possible to adapt our method also to accretion discs. The grid in Fig. 7 (left panel) allows to describe an optically thin or optically thick disc, as well as the polar wind and hot corona. The radiation from the underlying star and from the outer boundary region is naturally included here, too. We test this method for a case of the cataclysmic variable HT Cas. The model of this system is taken from Horne et al. (99), Williams (98), and Williams & Shipman (988). The detailed description of the input model, as well as some preliminary results are described in Korčáková et al. (25). Here we show only the intensity map of the accretion disc (see Fig. 7, right panel). 4 Conclusion A method of the solution of the radiative transfer equation in axial symmetry is presented here. Since the detailed description of this method, as well as its tests in a number of stellar atmosphere models appeared in Korčáková & Kubát (25a) already, we merely emphasized the main ideas of this method. The greatest advantage of this method is in the simultaneous solution of the radiative transfer
8 8 Title : will be set by the publisher longitudinal plane Fig. 7. left panel: The latitudinal grid for solving the radiative transfer problem in accretion discs. right panel: The intensity map of the accretion disc for the model corresponding to HT Cas. equation in static media and in media with velocity gradients. The case of large velocity gradients can be handled using a finer space grid. Results from the first calculations are presented here. This method is very useful for an accurate description of limb darkening. It allows to calculate line profiles from rotating stars with gravity darkening and differential rotation. This is very useful for rapidly rotating stars, where the convolution method cannot be used, and flux calculation from non-distorted stars is also inappropriate. For winds we can take advantage of the solution of the wind region together with the stellar atmosphere layers. Since the grid can be finer near the equator, it is possible to calculate the radiative transfer equation in disks. In this case we can include not only the radiation of the disk, but also the central object, the boundary region, hot corona, or fast polar wind, if it is present. We are now working on a generalization of our method, which will consist mainly of the possibility of treating additional opacity sources, and of improved treatment of NLTE calculations. The authors would like to thank Júlia Sokolovičová for her comments to the manuscript. This research was supported by grants 25/4/P224 (GA ČR) and B3635 (GA AV ČR). The Astronomical Institute Ondřejov is supported by a project AVZ35. References Boissé, P. 99, A&A 228, 483 Gray, D. F. 976, Observation and Analysis of Stellar Photospheres, John Wiley & Sons, New York Horne, K., Wood, J. H., & Stiening, R. F. 99, ApJ, 378, 27 Kneer, F., & Heasley, J. N. 979, A&A, 79, 4 Korčáková, D., Kubát, J. 23, A&A, 4, 49
9 Korčáková & Kubát: Radiative Transfer in Axial Symmetry 9 Korčáková, D., Kubát, J. 25a, A&A, 44, 75 Korčáková, D., Kubát, J. 25b, Mem S. A. It., in press Korčáková, D., Kubát, J., Krtička, J., & Šlechta, M. 24, in The A-Star Puzzle, IAU Symp. 224, J. Zverko, J. Žižňovský, S. J. Adelman & W. W. Weiss eds., Cambridge, Univ. Press, p. 533 Korčáková, D., Kubát, J., Kawka, A., 25, in 4th European Workshop on White Dwarfs, D. Koester & S. Moehler eds., ASP Conference Series, 334, 43 Lamers, H. J. G. L. M., & Cassinelli, J. P. 999, Introduction to Stellar Winds, Cambridge Univ. Press, Cambridge Sobolev, V. 946, Dvizhushchiesia obolochki zvedz, Leningr. Gos. Univ., Leningrad Williams, R., E. 98, ApJ, 235, 939 Williams, R., E., & Shipman, H. 988, ApJ, 326, 738
arxiv:astro-ph/ v1 16 Jun 2005
Astronomy & Astrophysics manuscript no. 992 November 5, 208 (DOI: will be inserted by hand later) Radiative transfer in moving media II. Solution of the radiative transfer equation in axial symmetry Daniela
More 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 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 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 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 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 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 informationInfluence of Mass Flows on the Energy Balance and Structure of the Solar Transition Region
**TITLE** ASP Conference Series, Vol. **VOLUME***, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Influence of Mass Flows on the Energy Balance and Structure of the Solar Transition Region E. H. Avrett and
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 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 informationStellar Magnetospheres part deux: Magnetic Hot Stars. Stan Owocki
Stellar Magnetospheres part deux: Magnetic Hot Stars Stan Owocki Key concepts from lec. 1 MagRe# --> inf => ideal => frozen flux breaks down at small scales: reconnection Lorentz force ~ mag. pressure
More informationA study of accretion disk wind emission
Mem. S.A.It. Vol. 83, 525 c SAIt 2012 Memorie della A study of accretion disk wind emission R. E. Puebla 1, M. P. Diaz 1, and D. J. Hillier 2 1 Departamento de Astronomia, Instituto de Astronomia, Geofísica
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 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 informationThree Dimensional Radiative Transfer in Winds of Massive Stars: Wind3D
Three Dimensional Radiative Transfer in Winds of Massive Stars: A. LOBEL 1 and R. BLOMME 1 arxiv:0707.3726v1 [astro-ph] 25 Jul 2007 1 Royal Observatory of Belgium, Ringlaan 3, Brussels, B-1180, Belgium
More informationExercise: A Toy Model for Dust-driven Winds
Astrofysikalisk dynamik, VT 00 Exercise: A Toy Model for Dust-driven Winds Susanne Höfner Department of Physics and Astronomy, Uppsala University Cool luminous giants stars, in particular pulsating AGB
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 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 informationThe Connection between Planets and White Dwarfs. Gilles Fontaine Université de Montréal
The Connection between Planets and White Dwarfs Gilles Fontaine Université de Montréal Astrophysical Context 1) Planets around white dwarfs (WD s) have yet to be found through the more standard methods:
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 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 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 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 informationAnalyzing X-Ray Pulses from Stellar Cores Pencil & Paper Version
Analyzing X-Ray Pulses from Stellar Cores Pencil & Paper Version Purpose: To determine if two end products of stellar evolution GK Per and Cen X-3 could be white dwarfs or neutron stars by calculating
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 informationAdvanced Stellar Astrophysics
v Advanced Stellar Astrophysics William K. Rose University of Maryland College Park CAMBRIDGE UNIVERSITY PRESS Contents Preface xiii Star formation and stellar evolution: an overview 1 1 A short history
More informationB B E D B E D A C A D D C
B B A C E E E E C E D E B B A D B E E A E E B C C A B B E D B E D A C A D D C E D Assigned Reading Read Chapters 8.1 and 8.2 Colonel Cady Coleman, Astronaut: Lessons from Space Lead Straight Back to Earth,
More informationToday The Sun. Events
Today The Sun Events Last class! Homework due now - will count best 5 of 6 Final exam Dec. 20 @ 12:00 noon here Review this Course! www.case.edu/utech/course-evaluations/ The Sun the main show in the solar
More informationAGN Central Engines. Supermassive Black Holes (SMBHs) Masses and Accretion Rates SMBH Mass Determinations Accretion Disks
AGN Central Engines Supermassive Black Holes (SMBHs) Masses and Accretion Rates SMBH Mass Determinations Accretion Disks 1 Supermassive Black Holes Need to generate L > 10 43 ergs/sec inside radius < 10
More informationRadiation-hydrodynamic Models for ULXs and ULX-pulsars
Radiation-hydrodynamic Models for ULXs and ULX-pulsars Tomohisa KAWASHIMA Division of Theoretical Astrophysics, NAOJ in collaboration with Ken OHSUGA, Hiroyuki TAKAHASHI (NAOJ) Shin MINESHIGE, Takumi OGAWA
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 informationConvection When the radial flux of energy is carried by radiation, we derived an expression for the temperature gradient: dt dr = - 3
Convection When the radial flux of energy is carried by radiation, we derived an expression for the temperature gradient: dt dr = - 3 4ac kr L T 3 4pr 2 Large luminosity and / or a large opacity k implies
More informationHigh-Energy Astrophysics Lecture 6: Black holes in galaxies and the fundamentals of accretion. Overview
High-Energy Astrophysics Lecture 6: Black holes in galaxies and the fundamentals of accretion Robert Laing Overview Evidence for black holes in galaxies and techniques for estimating their mass Simple
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 informationBROAD SPECTRAL LINE AND CONTINUUM VARIABILITIES IN QSO SPECTRA INDUCED BY MICROLENSING:METHODS OF COMPUTATION
Proceedings of the IX Bulgarian-Serbian Astronomical Conference: Astroinformatics (IX BSACA) Sofia, Bulgaria, July -,, Editors: M. K. Tsvetkov, M. S. Dimitrijević, O. Kounchev, D. Jevremović andk. Tsvetkova
More informationCentrifugal forces. Equipotential surfaces. Critical rotation velocity ROTATION AND STELLAR STRUCTURE. STELLAR ROTATION and EVOLUTION.
STELLAR ROTATION and EVOLUTION Henny J.G.L.M. Lamers Astronomical Institute, Utrecht University 22/09/09 Lect 1: Rotation and stellar structure 22/09/09 Lect 2: Rotation and stellar winds 24/09/09 Lect
More informationTypes of Stars 1/31/14 O B A F G K M. 8-6 Luminosity. 8-7 Stellar Temperatures
Astronomy 113 Dr. Joseph E. Pesce, Ph.D. The Nature of Stars For nearby stars - measure distances with parallax 1 AU d p 8-2 Parallax A January ³ d = 1/p (arcsec) [pc] ³ 1pc when p=1arcsec; 1pc=206,265AU=3
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 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 informationNUMERICAL METHODS IN ASTROPHYSICS An Introduction
-1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,
More informationCharacteristic temperatures
Characteristic temperatures Effective temperature Most sources are only roughly blackbodies (if that). So we integrate the flux over frequency and define: F = I cosθ dω d = σ T e 4 i.e. a source of effective
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 informationCosine of emission angle: Energy (kev)
EFFECTS OF STRONG GRAVITY ON THE X-RAY SPECTRA OF AGN AND BHCs A. MARTOCCHIA 1, V. KARAS 2 and G. MATT 3 (1) SISSA-ISAS Via Beirut 2/4, I-34014 Trieste (Italy) (2) Astronomical Institute, Charles University
More informationOUTLINE: P. Kotrč (1), P. Heinzel (1) and O. Procházka (2)
On measurements of continuum flux in solar flares. Instrument and first results. P. Kotrč (1), P. Heinzel (1) and O. Procházka (2) (1) - Astronomical Institute, AS CR, v.v.i. Ondřejov, Czech Republic (2)
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 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 informationThe Sun: A Star of Our Own ASTR 2110 Sarazin
The Sun: A Star of Our Own ASTR 2110 Sarazin Sarazin Travel Wednesday, September 19 afternoon Friday, September 21 Will miss class Friday, September 21 TA Molly Finn will be guest lecturer Cancel Office
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 informationPulsars ASTR2110 Sarazin. Crab Pulsar in X-rays
Pulsars ASTR2110 Sarazin Crab Pulsar in X-rays Test #2 Monday, November 13, 11-11:50 am Ruffner G006 (classroom) Bring pencils, paper, calculator You may not consult the text, your notes, or any other
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 informationAccretion disks. AGN-7:HR-2007 p. 1. AGN-7:HR-2007 p. 2
Accretion disks AGN-7:HR-2007 p. 1 AGN-7:HR-2007 p. 2 1 Quantitative overview Gas orbits in nearly circular fashion Each gas element has a small inward motion due to viscous torques, resulting in an outward
More informationProblem set: solar irradiance and solar wind
Problem set: solar irradiance and solar wind Karel Schrijver July 3, 203 Stratification of a static atmosphere within a force-free magnetic field Problem: Write down the general MHD force-balance equation
More informationThe Solar Temperature Minimum and Chromosphere
**TITLE** ASP Conference Series, Vol. **VOLUME***, **YEAR OF PUBLICATION** **NAMES OF EDITORS** The Solar Temperature Minimum and Chromosphere Eugene H. Avrett Smithsonian Astrophysical Observatory, Harvard-Smithsonian
More informationStar Formation and Protostars
Stellar Objects: Star Formation and Protostars 1 Star Formation and Protostars 1 Preliminaries Objects on the way to become stars, but extract energy primarily from gravitational contraction are called
More informationarxiv:astro-ph/ v1 25 Mar 1998
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 8 August 2018 (MN LATEX style file v1.4) A new method of determining the inclination angle in interacting binaries Tariq Shahbaz University of Oxford,
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 informationGravitational light bending prevents gamma-gamma absorption in gravitational lenses
Gravitational light bending prevents gamma-gamma absorption in gravitational lenses NWU, South Africa E-mail: 26559226@nwu.ac.za Markus Böttcher NWU, South Africa E-mail: Markus.Bottcher@nwu.ac.za The
More informationAstrophysical Quantities
Astr 8300 Resources Web page: http://www.astro.gsu.edu/~crenshaw/astr8300.html Electronic papers: http://adsabs.harvard.edu/abstract_service.html (ApJ, AJ, MNRAS, A&A, PASP, ARAA, etc.) General astronomy-type
More informationOutflowing disk formation in B[e] supergiants due to rotation and bi-stability in radiation driven winds
A&A 437, 929 933 (2005) DOI: 10.1051/0004-6361:20052686 c ESO 2005 Astronomy & Astrophysics Outflowing disk formation in B[e] supergiants due to rotation and bi-stability in radiation driven winds M. Curé
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 informationOverview spherical accretion
Spherical accretion - AGN generates energy by accretion, i.e., capture of ambient matter in gravitational potential of black hole -Potential energy can be released as radiation, and (some of) this can
More information9.3 Worked Examples Circular Motion
9.3 Worked Examples Circular Motion Example 9.1 Geosynchronous Orbit A geostationary satellite goes around the earth once every 3 hours 56 minutes and 4 seconds, (a sidereal day, shorter than the noon-to-noon
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 informationRadiative MHD. in Massive Star Formation and Accretion Disks. Rolf Kuiper, Hubert Klahr, Mario Flock, Henrik Beuther, Thomas Henning
Radiative MHD in Massive Star Formation and Accretion Disks, Hubert Klahr, Mario Flock, Henrik Beuther, Thomas Henning, Radiative MHD with Makemake and Pluto : We developed a fast 3D frequency-dependent
More informationThe Sun. October 21, ) H-R diagram 2) Solar Structure 3) Nuclear Fusion 4) Solar Neutrinos 5) Solar Wind/Sunspots
The Sun October 21, 2002 1) H-R diagram 2) Solar Structure 3) Nuclear Fusion 4) Solar Neutrinos 5) Solar Wind/Sunspots Review Blackbody radiation Measuring stars distance luminosity brightness and distance
More informationThe Sun Our Extraordinary Ordinary Star
The Sun Our Extraordinary Ordinary Star 1 Guiding Questions 1. What is the source of the Sun s energy? 2. What is the internal structure of the Sun? 3. How can astronomers measure the properties of the
More informationAn Overview of the Details
The Sun Our Extraordinary Ordinary Star 1 Guiding Questions 1. What is the source of the Sun s energy? 2. What is the internal structure of the Sun? 3. How can astronomers measure the properties of the
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 informationB.V. Gudiksen. 1. Introduction. Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della
Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della À Ø Ò Ø ËÓÐ Ö ÓÖÓÒ B.V. Gudiksen Institute of Theoretical Astrophysics, University of Oslo, Norway e-mail:boris@astro.uio.no Abstract. The heating mechanism
More informationarxiv:astro-ph/ v3 27 May 2005
X-ray Variability of AGN and the Flare Model arxiv:astro-ph/0410079v3 27 May 2005 R.W. Goosmann 1, B. Czerny 2, A.-M. Dumont 1, M. Mouchet 1, and A. Różańska 2 1 LUTH, Observatoire de Paris, Meudon, France
More informationO 5+ at a heliocentric distance of about 2.5 R.
EFFECT OF THE LINE-OF-SIGHT INTEGRATION ON THE PROFILES OF CORONAL LINES N.-E. Raouafi and S. K. Solanki Max-Planck-Institut für Aeronomie, 37191 Katlenburg-Lindau, Germany E-mail: Raouafi@linmpi.mpg.de;
More informationThe Sun. Chapter 12. Properties of the Sun. Properties of the Sun. The Structure of the Sun. Properties of the Sun.
Chapter 12 The Sun, Our Star 1 With a radius 100 and a mass of 300,000 that of Earth, the Sun must expend a large amount of energy to withstand its own gravitational desire to collapse To understand this
More informationarxiv:astro-ph/ v1 4 Oct 2004
X-ray Variability of AGN and the Flare Model arxiv:astro-ph/0410079v1 4 Oct 2004 R.W. Goosmann 1, B. Czerny 2, A.-M. Dumont 1, M. Mouchet 1, and A. Różańska 2 1 LUTH, Observatoire de Paris, Meudon, France
More informationHigh-Energy Astrophysics
Oxford Physics: Part C Major Option Astrophysics High-Energy Astrophysics Garret Cotter garret@astro.ox.ac.uk Office 756 DWB Michaelmas 2011 Lecture 7 Today s lecture: Accretion Discs Part I The Eddington
More informationThe total luminosity of a disk with the viscous dissipation rate D(R) is
Chapter 10 Advanced Accretion Disks The total luminosity of a disk with the viscous dissipation rate D(R) is L disk = 2π D(R)RdR = 1 R 2 GM Ṁ. (10.1) R The disk luminosity is half of the total accretion
More informationExplain how the sun converts matter into energy in its core. Describe the three layers of the sun s atmosphere.
Chapter 29 and 30 Explain how the sun converts matter into energy in its core. Describe the three layers of the sun s atmosphere. Explain how sunspots are related to powerful magnetic fields on the sun.
More informationAn Overview of the Details
Guiding Questions The Sun Our Extraordinary Ordinary Star 1. What is the source of the Sun s energy? 2. What is the internal structure of the Sun? 3. How can astronomers measure the properties of the Sun
More informationF q. Gas at radius R (cylindrical) and height z above the disk midplane. F z. central mass M
Accretion Disks Luminosity of AGN derives from gravitational potential energy of gas spiraling inward through an accretion disk. Derive structure of the disk, and characteristic temperatures of the gas.
More informationRedshifted Broad Absorption Troughs in Quasars
Redshifted Broad Absorption Troughs in Quasars Pat Hall, York University with Niel Brandt and Nur Filiz Ak (PSU), Patrick Petitjean (IAP), and the SDSS-III/BOSS team If due to fallback or rotation, challenge
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 informationBinary Stars (continued) ASTR 2120 Sarazin. γ Caeli - Binary Star System
Binary Stars (continued) ASTR 2120 Sarazin γ Caeli - Binary Star System Visual Binaries: Types of Binary Stars Spectroscopic Binaries: Eclipsing Binaries: Periodic changes in brightness, stars block one
More informationScience Opportunities in Stellar Physics. Douglas R. Gies CHARA, Georgia State Univ., and the Stellar Physics Working Group
Science Opportunities in Stellar Physics Douglas R. Gies CHARA, Georgia State Univ., gies@chara.gsu.edu and the Stellar Physics Working Group General Themes! Fundamental properties! Interior structure
More informationMonte Carlo Simulator to Study High Mass X-ray Binary System
SLAC-PUB-11350 Monte Carlo Simulator to Study High Mass X-ray Binary System S. Watanabe, F. Nagase, T. Takahashi ISAS/JAXA, Sagamihara, Kanagawa 229-8510, Japan M. Sako, S.M. Kahn KIPAC/Stanford, Stanford,
More informationCambridge University Press Advanced Stellar Astrophysics William K. Rose Frontmatter More information
In the last two decades, remarkable progress has been made in understanding stars. This graduate-level textbook provides a systematic, self-contained and lucid introduction to the physical processes and
More informationPhys 100 Astronomy (Dr. Ilias Fernini) Review Questions for Chapter 8
Phys 100 Astronomy (Dr. Ilias Fernini) Review Questions for Chapter 8 MULTIPLE CHOICE 1. Granulation is caused by a. sunspots. * b. rising gas below the photosphere. c. shock waves in the corona. d. the
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 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 informationThe Sun. the main show in the solar system. 99.8% of the mass % of the energy. Homework due next time - will count best 5 of 6
The Sun the main show in the solar system 99.8% of the mass 99.9999...% of the energy 2007 Pearson Education Inc., publishing as Pearson Addison-Wesley Homework due next time - will count best 5 of 6 The
More informationWeek 8: Stellar winds So far, we have been discussing stars as though they have constant masses throughout their lifetimes. On the other hand, toward
Week 8: Stellar winds So far, we have been discussing stars as though they have constant masses throughout their lifetimes. On the other hand, toward the end of the discussion of what happens for post-main
More informationHigh Energy Astrophysics
High Energy Astrophysics Accretion Giampaolo Pisano Jodrell Bank Centre for Astrophysics - University of Manchester giampaolo.pisano@manchester.ac.uk April 01 Accretion - Accretion efficiency - Eddington
More informationEvolution of High Mass stars
Evolution of High Mass stars Neutron Stars A supernova explosion of a M > 8 M Sun star blows away its outer layers. The central core will collapse into a compact object of ~ a few M Sun. Pressure becomes
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 informationPropagation of very high energy γ-rays inside massive binaries LS 5039 and LSI
Mon. Not. R. Astron. Soc. 368, 579 591 (2006) doi:10.1111/j.1365-2966.2006.10121.x Propagation of very high energy γ-rays inside massive binaries LS 5039 and LSI +61 303 W. Bednarek Department of Experimental
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 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 informationFUSE Observations of a Mira Variable Star
Proceedings of 12th Cambridge Workshop on Cool Stars, Stellar Systems, & The Sun, 2003 University of Colorado. FUSE Observations of a Mira Variable Star Donald G. Luttermoser 1, Michael W. Castelaz 2 Abstract.
More informationCataclysmic variables
Cataclysmic variables Sander Bus Kapteyn Astronomical Institute Groningen October 6, 2011 Overview Types of cataclysmic stars How to form a cataclysmic variable X-ray production Variation in outburst lightcurve,
More informationAGN Physics of the Ionized Gas Physical conditions in the NLR Physical conditions in the BLR LINERs Emission-Line Diagnostics High-Energy Effects
AGN Physics of the Ionized Gas Physical conditions in the NLR Physical conditions in the BLR LINERs Emission-Line Diagnostics High-Energy Effects 1 Evidence for Photoionization - continuum and Hβ luminosity
More informationX-ray Spectroscopy of Massive Star Winds: Shocks, Mass-Loss Rates, and Clumping
X-ray Spectroscopy of Massive Star Winds: Shocks, Mass-Loss Rates, and Clumping David Cohen Department of Physics & Astronomy Swarthmore College Maurice Leutenegger (GSFC), Jon Sundqvist (Madrid), Stan
More informationSolar photosphere. Michal Sobotka Astronomical Institute AS CR, Ondřejov, CZ. ISWI Summer School, August 2011, Tatranská Lomnica
Solar photosphere Michal Sobotka Astronomical Institute AS CR, Ondřejov, CZ ISWI Summer School, August 2011, Tatranská Lomnica Contents General characteristics Structure Small-scale magnetic fields Sunspots
More information