meeting 2024 01 24 n315 - JacobPilawa/TriaxSchwarzschild_wiki_5 GitHub Wiki

Context

Final Set of Diagnostics

I've made the master spectra plot showing the fits, spectra, residuals for our 249 bins:
Plot
And here's the point-symmetric and non-point symmetric kinematics/maps. I've included Irina's published kinematic plots for reference. Note that her results use the symmetric binning scheme and are not point-symmetric.
Irina's Kinematics (MASSIVE XIII)
non point-symmetric point-symmetric

Test Model

  • I've worked through Emily's new modeling workflow, and I think I got things working! I ran a single model with the following important parameters:

    • nDither = 1: just to make it run quickly/for testing purposes

    • BH = 2e9

    • M/L = 2.0

    • gnfw profile with:

      • rho0 = 2.4e9 (randomly chosen, same as N57)
      • rs = 15 kpc (randomly chosen, same as N57)

    • angles randomly chosen:

      • theta = 86.3
      • phi = 84.3
      • psi = -89.9
      • This translates to (T, Tmaj, Tmin) = (0.416, 0.01, 0.005) [I realize now I should have picked something better, but this is fine for this test model.]
      • (u,p,q) = (0.9988, 0.8754, 0.6622)

    • nI2 = 18 for loop library, nI2 = 9 for box

    • nI3 = 9 for both libraries

  • The resulting kinematic fits actually look quite nice given the random parameter choice and the fact we're using nDith =1.

    • Note that this was simply a test model, and there are a few other things I need to do. Specifically, I still:
      • need to include dummy momemnts for h9-h12
      • decide what to do about the few overlapping Mitchell bins
Kinematics
Beta
  • I also plotted the mass in each radial shell as a function of radius, and noticed that this might need some revision. Specifically, the fit mass and mge mass track each other perfectly, but I think we need to increase the distance of the outermost radial shell since the mass has not seemed to level off yet. Additionally, summing the masses over the 9 radial shells gives ~0.86 instead of ~1.0.
  • I don't have a sense of what increase we do need to do this, but maybe a factor of 2 in the outermost radial shell distance?
Mass vs. Shell
I've also plotted a heatmap of the chi2 so we can get a sense of if there are any systematic issues/weird bins. Note that the h6 Mitchell data seems exceptionally poorly fit.

⚠️ **GitHub.com Fallback** ⚠️