SYNSPEC SYNSPEC p. 1 SYNSPEC calculation of synthetic spectra of - PowerPoint PPT Presentation

SYNSPEC SYNSPEC – p. 1

SYNSPEC calculation of synthetic spectra of stellar atmospheres and disks input model atmospheres TLUSTY or ATLAS Hubeny, I., Štefl, S., Harmanec, P ., 1985, Bull. Astron. Inst. Czechosl. 36, 214 http://www.physics.muni.cz/˜krticka/synspec49.tgz SYNSPEC – p. 2

Equation to solve radiative transfer equation µ ∂ ∂zI ( z,µ,ν ) = η ( z,ν ) − χ ( z,ν ) I ( z,µ,ν ) fixed structure of atmosphere SYNSPEC – p. 3

Files with the source code http://nova.astro.umd.edu/Synspec49/synspec.html synspec45.f code PARAMS.FOR some basic settings LINDAT.FOR, MODELP .FOR rotin3.f code for the rotational convolution compilation: ifort [-O3 ] synspec45.f SYNSPEC – p. 4

Input files http://nova.astro.umd.edu/ file 1: the same as in TLUSTY file 5: the same as in TLUSTY files of the model ions (as in TLUSTY) nonstandard settings (as in TLUSTY) file 55: basic parameters of the calculation file 56: change of the chemical composition file 19: line list file 8: input model atosphere SYNSPEC – p. 5

File 55 (fort.55) 0 30 1 !IMODE, IDSTD, IPRIN 0 0 0 0 !INPUT MODEL, INTRPL, ICHANGE, ICHEMC 0 !IOPHLI 0 0 0 0 1 !IFREQ, INLTE, ICONTL, INORIG, IFHE2 0 0 0 0 0 0 0 !ILYCS,IBVCS,IHE1,IHE4471,IHE2UV,IH 4000 7000 10 10 1d-6 0.001 ! ALAM0,ALAST,CUTOF0, CUTOFFS,RELOP,SPACE 2.0 ! turbulent velocity (in km/s) SYNSPEC – p. 6

Input parameters IMODE : mode of calculation IDSTD : index of the standard depth (unity optical depth) IPRIN : amount of output INPUT MODEL : (0=ATLAS, 1=TLUSTY) IOPHLI : treatment of hydrogen lines INLTE : treatment of NLTE ALAM0,ALAST : wavelength limits RELOP : minimum relative contribution of lines VTB : turbulent velocity in km/s SYNSPEC – p. 7

Example tun /home/krticka/prakaf/synspec/synspec file syn : DOMA=/home/krticka JOBID=$1 MODEL=$2 LINES=$3 ln -s $MODEL.7 fort.8 ln -s $JOBID.55 fort.55 ln -s $DOMA/GF/$LINES fort.19 $DOMA/prakaf/synspec/synspec/synspec45 <$MODEL.5 >! $JOBID.6 cp fort.7 $JOBID.7 cp fort.17 $JOBID.17 cp fort.12 $JOBID.12 cp fort.18 $JOBID.18 cp fort.13 $JOBID.13 cp fort.10 $JOBID.10 SYNSPEC – p. 8 rm fort.*

Example run /home/krticka/prakaf/synspec/synspec file syn : run syn slun.s slun gf99.dat input files slun.5, slun.7, slun.s.55, gf99.dat output files slun.s.6, slun.s.7, slun.s.17 ... SYNSPEC – p. 9

Output files file 6: basic model description file 7: emergent flux in lines file 17: emergent flux in continuum file 12: line identification (except H and He II) CGS! SYNSPEC – p. 10

File 6 (fort.6) general model output tables with input data error messages SYNSPEC – p. 11

File 7 (fort.7) wavelength-dependent emergent flux ( λ – Å, H λ – erg cm − 2 s − 1 Å − 1 ) SYNSPEC – p. 12

File 17 (fort.17) wavelength-dependent emergent flux in continuum ( λ – Å, H λ – erg cm − 2 s − 1 Å − 1 ) SYNSPEC – p. 13

File 12 (fort.12) line identification table, estimate of the equivalent width, line strength ( * ) SYNSPEC – p. 14

Rotational convolution code rotin3.f input files file 5: parameters of the calculation file 7, file 17 (SYNSPEC output) output file file 11: convolved spectrum (rotation+instrumental) SYNSPEC – p. 15

Input file 5 ’fort.7’ ’fort.17’ ’fort.11’ ! filenames for the 70.0 0.00 0.001 ! vrot, chard, step 0.03001 0.001 ! fwhm (instrumental 4100 6900 1 ! lam0, lam1, irel ( names of input and output files v rot sin i in km/s, 0., step FWHM of instrumental profile wavelength interval, absolute/relative spectrum SYNSPEC – p. 16

Example run file rot: JOBID=$1 SPEK=$2 ln -s $SPEK.7 fort.7 ln -s $SPEK.17 fort.17 rotin3 <$JOBID.5 rm fort.7 fort.17 mv fort.11 $SPEK.11 execute: rot rvs b.s SYNSPEC – p. 17

Example of visual spectrum 1.00 0.95 relative intensity 0.90 Ti II Ca I Fe II Fe II Ti II Fe I 0.85 Ti II Ti II Sc II Fe II Ti II 0.80 0.75 fit observed Fe I Mg II 0.70 4400 4420 4440 4460 4480 4500 wavelength λ [Å] SYNSPEC – p. 18

Example of UV spectrum 2.5 × 10 -10 prediction φ =0.005 flux F λ [erg cm -2 s -1 A -1 ] observation (IUE) 2.0 × 10 -10 1.5 × 10 -10 1.0 × 10 -10 5.0 × 10 -11 C II Si II Si II Si II Al II Al II Fe II Fe II Fe II Fe II Mg II 0 1200 1400 1600 1800 2000 2200 2400 2600 2800 3000 wavelength λ [ ° A] SYNSPEC – p. 19

Abundance analysis minimalization − I syn � 2 I obs � χ 2 = � i i σ 2 i λ i SYNSPEC – p. 20