A New High-Precision Correction Method of Temperature Distribution in Model Stellar Atmospheres

A. Sapar 1 , R. Poolamäe 1 ,  and L. Sapar 1
  • 1 Tartu Observatory, 61602 Tõravere, Estonia


The main features of the temperature correction methods, suggested and used in modeling of plane-parallel stellar atmospheres, are discussed. The main features of the new method are described. Derivation of the formulae for a version of the Unsöld-Lucy method, used by us in the SMART (Stellar Model Atmospheres and Radiative Transport) software for modeling stellar atmospheres, is presented. The method is based on a correction of the model temperature distribution based on minimizing differences of flux from its accepted constant value and on the requirement of the lack of its gradient, meaning that local source and sink terms of radiation must be equal. The final relative flux constancy obtainable by the method with the SMART code turned out to have the precision of the order of 0.5 %. Some of the rapidly converging iteration steps can be useful before starting the high-precision model correction. The corrections of both the flux value and of its gradient, like in Unsöld-Lucy method, are unavoidably needed to obtain high-precision flux constancy. A new temperature correction method to obtain high-precision flux constancy for plane-parallel LTE model stellar atmospheres is proposed and studied. The non-linear optimization is carried out by the least squares, in which the Levenberg-Marquardt correction method and thereafter additional correction by the Broyden iteration loop were applied. Small finite differences of temperature (δT/T = 10−3) are used in the computations. A single Jacobian step appears to be mostly sufficient to get flux constancy of the order 10−2 %. The dual numbers and their generalization – the dual complex numbers (the duplex numbers) – enable automatically to get the derivatives in the nilpotent part of the dual numbers. A version of the SMART software is in the stage of refactorization to dual and duplex numbers, what enables to get rid of the finite differences, as an additional source of lowering precision of the computed results.

If the inline PDF is not rendering correctly, you can download the PDF file here.

  • Aret A., Sapar A. 2002, AN, 323, 21

  • Aret A., Sapar A., Poolamäe R., Sapar L. 2008, in The Art of Modelling Stars in the 21st Century, IAU Symp. 252, 41

  • Broyden C. G. 1965, Math. Comp., 19, 577

  • Castelli F., Kurucz R. L. 2003, in Modelling of Stellar Atmospheres, IAU Symp. 210, Poster A20; arXiv:astro-ph/0405087v1, 2004

  • Castelli F. 2005, MmSAI, 8, 25

  • Gustafsson B., Edvardsson B., Eriksson K. et al. 2008, A&A, 486, 951

  • Hauschildt P. H., Barman T. S., Baron E., Allard F. 2003, in Stellar Atmosphere Modeling, ASPC, 228, 227

  • Kirby E. N. 2011, PASP, 123, 531

  • Koesterke L., Hamann W.-R., Kosmol P. 1992, A&A, 255, 490

  • Kurucz R. L. 1970, ATLAS: A Computer Program for Calculating Model Stellar Atmospheres, SAO Special Report, No. 309

  • Levenberg K. 1944, Quart. Appl. Math., 2, 164

  • Lourakis M. I. A. 2004, http://www.ics.forth.gr/~lourakis/levmar/

  • Lucy L. B. 1964, SAO Special Report, No. 167, 93

  • Marquardt D. 1963, J. Soc. Indust. Appl. Math., 11 (2), 431

  • Mihalas D., Auer L. H. 1969, ApJ, 158, 641

  • Sakhibullin N. 1997, Methods of Modelling in Astrophysics, Fen Press, Kazan (in Russian)

  • Sapar A., Poolamäe R. 2003, in Stellar Atmosphere Modeling, ASPC, 288, 95

  • Sapar A., Aret A., Sapar L., Poolamäe R. 2007a, in Spectroscopic Methods in Modern Astrophysics, INASAN, p. 220 (in Russian)

  • Sapar A., Poolamäe R., Sapar L., Aret A. 2007b, in Spectroscopic Methods in Modern Astrophysics, INASAN, p. 236 (in Russian)

  • Sapar A., Aret A., Sapar L., Poolamäe R. 2008, in Precision Spectroscopy in Astrophysics: ESO Astrophysics Symposia, Lisbon Aveiro, p. 145

  • Sapar A., Aret A., Sapar L., Poolamäe R. 2009, New Astronomy Reviews, 53, 240

  • Sapar A., Poolamäe R. 2012, Baltic Astronomy, 21, 243

  • Unsöld A. 1951, Naturwiss., 38, 525

  • Unsöld A. 1955, Physik der Sternatmosphären, Springer-Verlag

  • Werner K., Dreizler S. J. 1999, Comput. Appl. Math., 109, 65

  • Werner K., Deetjen J. L., Dreizler S. et al. 2003, ASPC, 288, 31


Journal + Issues

Open Astronomy is a peer-reviewed, fully open access online journal, that publishes research, reviews and news spanning all aspects of astronomy and astrophysics. Open Astronomy (formerly Baltic Astronomy) was founded by the Institute of Theoretical Physics and Astronomy (Lithuania) in 1992 and acquired by De Gruyter in 2016.