ISM Scattering

Ref

Steps

  1. Set up scattering model
  2. Initialize parameters
  3. Sgr A* Parameters
  4. Create Images
  5. Define Flux Cutoffs
  6. Remove any (u,v) points with noise above thermal noise floor
  7. Plot :)

Step 1: Set up scattering model

sm=so.ScatteringModel(scatt_alpha=1.38) #change alpha to the value that you want; our paper gives some values and motivations
scatt_kernel=sm.Ensemble_Average_Kernel_Visibility #forms the scattering kernel in the Ensemble Average regime

Step 2: Initialize Parameters

f_Hz=320e9 #observing frequency
c=3e10 # speed of light in cm/s
wavelength_cm=c/f_Hz

u=np.linspace(35e9,-35e9,200)
v=np.linspace(35e9,-35e9,200)
uu,vv=np.meshgrid(u,v) #array of (u,v) coordinates

Step 3: Sgr A* Parameters

geometric_model = eh.model.Model()
geometric_model = geometric_model.add_thick_mring(F0 = 3.5, d = 50.*eh.RADPERUAS, alpha = 10.*eh.RADPERUAS, x0 = 0.0*eh.RADPERUAS, y0 = 0.0*eh.RADPERUAS)
geometric_model = geometric_model.add_thick_mring(F0 = 0.5, d = 50.*eh.RADPERUAS, alpha = 1.*eh.RADPERUAS, x0 = 0.0*eh.RADPERUAS, y0 = 0.0*eh.RADPERUAS)
geometric_model.rf = f_Hz

ISM Scattering

By Ref Bari

ISM Scattering

  • 5

More from Ref Bari