Emission lines from post-AGB UV radiation

Aperta
#217 0 commenti 0 reazioni 0 assegnatari Vedi su GitHub

Nessuno ha ancora preso questa issue.

Valutazione

Difficoltà
4/5
Tempo stimato
3-5 giorni
Idoneità per principianti
35/100
Tipo di issue
Bug
Chiarezza
Da chiarire
Stato di attività
Ferma
Stack tecnologico
python
Ambito
data

Direzione di ricerca

Iniziate eseguendo lo script Python fornito e analizzando la chiamata a get_spectrum insieme ai parametri pAGB, nebemlineinspec e add_neb_emission. Confrontate gli spettri risultanti per i tre valori di pAGB; l’issue è completa quando è stata identificata la causa dei flussi delle righe di emissione pari a zero e sono stati documentati il comportamento previsto o la configurazione richiesta.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

Hi,

I am running some models to test the UV contribution from the post-AGB stars to the emission lines of a galaxy. To test this, I use SSPs at the age of 3Gyrs (also tried for 0.1Gyrs and 5Gyrs). Models have no attenuation. The SSPs have solar metallicity. I test 3 separate instances which all have the same parameter values, except that each time I change the pAGB fraction to be 0, 0.5 and 1. I can see that I get ionizing radiation by turning on post-AGB stars.

The problem is that I get emission line fluxes that are all 0. Shouldn't this UV give emission lines for pAGB = 0.5 and pAGB = 1 ? Am I missing something obvious here?

From this tests I get the following plot:
3Gyr_SSP_no_att_add_neb_pAGB
0p1Gyr_SSP_no_att_add_neb_pAGB

The srcipt I run is the following:

import fsps
import matplotlib.pyplot as plt
import numpy as np

sp = fsps.StellarPopulation(zcontinuous=1)

sp.params['sfh'] = 0
sp.params['dust1'] = 0.0
sp.params['dust2'] = 0.0
sp.params['add_neb_emission'] = True
sp.params['add_dust_emission'] = True
sp.params['nebemlineinspec'] = True
sp.params['logzsol'] = 0.0
sp.params['gas_logz'] = sp.params['logzsol']

t = 3.0
for p in [0.0,0.5,1.0]:
	sp.params['pagb'] = p
	wave, flux = sp.get_spectrum(tage=t,peraa = True)
	plt.loglog(wave, flux,label='pAGB = ' + str(p))
	plt.xlabel('Wavelength (Angstroms)')
	plt.ylim([1e-7, 1e-2])
	plt.xlim([1e2, 1e5])
	plt.ylabel('Flux')
	plt.title('tage = ' + str(t)+'Gyr')
	plt.legend(fontsize=14)
plt.show()

Best,
Babis

Lingua principale
Python
Stelle
76
Fork
45
Nessuna PR unita negli ultimi 30g

Guida per i contributori

Nessuna guida per i contributori indicizzata per questo repository

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di dfm/python-fsps

Tutte le issue di dfm/python-fsps

Issue simili

Altre issue su Python

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.