A How do I determine the surface magnitude of a galaxy from a SED?

  • A
  • Thread starter Thread starter AdrianD
  • Start date Start date
  • Tags Tags
    Galaxy
AI Thread Summary
To determine the surface magnitude of a galaxy from its spectral energy distribution (SED), one can use the BlackBodyNorm1D model in Python's synphot library to generate flux data. By calculating the area of a circular aperture in square arcseconds, the surface brightness can be derived by dividing the flux by this area. The resulting surface brightness is expressed in FLAM per square arcsecond. Adjusting the aperture diameter allows for accurate measurements based on the source's angular size. This method effectively provides the surface brightness needed for further analysis of the galaxy's characteristics.
AdrianD
Messages
7
Reaction score
2
Let's say I have a galaxy SED like this:

[CODE lang="python" title="SEDgen"]from synphot import SourceSpectrum
from synphot.models import BlackBodyNorm1D
sp = SourceSpectrum(BlackBodyNorm1D, temperature=6170.4796)

wave = sp.waveset
flux = sp(wave) #Photlam
flux_flam = flux.to(u.erg / u.s / u.angstrom/u.cm**2,equivalencies=u.spectral_density(wave))


plt.plot(wave, flux_flam)
plt.xlim(1000, 30000)
plt.xlabel('Wavelength (Angstrom)')
plt.ylabel('Flux [FLAM]')[/CODE]

How do I get surface brightness from this?
 
Astronomy news on Phys.org
Python:
import numpy as np

# Assuming a circular aperture with a diameter in arcseconds
aperture_diameter_arcsec = 1.0  # Example: 1 arcsecond diameter aperture

# Calculate the area of the circular aperture in square arcseconds
aperture_area_arcsec2 = np.pi * (aperture_diameter_arcsec / 2)**2

# Calculate the surface brightness in FLAM per square arcsecond
surface_brightness_flam_arcsec2 = flux_flame / aperture_area_arcsec2

# Plot the surface brightness
plt.plot(wave, surface_brightness_flam_arcsec2)
plt.xlim(1000, 30000)
plt.xlabel('Wavelength (Angstrom)')
plt.ylabel('Surface Brightness [FLAM/arcsec^2]')
Replace the aperture_diameter_arcsec variable with the appropriate angular size of your source. This will give you the surface brightness inFLAM per square arcsecond.
 
Last edited by a moderator:
Today at about 4:30 am I saw the conjunction of Venus and Jupiter, where they were about the width of the full moon, or one half degree apart. Did anyone else see it? Edit: The moon is 2,200 miles in diameter and at a distance of 240,000 miles. Thereby it subtends an angle in radians of 2,200/240,000=.01 (approximately). With pi radians being 180 degrees, one radian is 57.3 degrees, so that .01 radians is about .50 degrees (angle subtended by the moon). (.57 to be more exact, but with...
This thread is dedicated to the beauty and awesomeness of our Universe. If you feel like it, please share video clips and photos (or nice animations) of space and objects in space in this thread. Your posts, clips and photos may by all means include scientific information; that does not make it less beautiful to me (n.b. the posts must of course comply with the PF guidelines, i.e. regarding science, only mainstream science is allowed, fringe/pseudoscience is not allowed). n.b. I start this...
Back
Top