meeting 2024 06 11 n410 - JacobPilawa/TriaxSchwarzschild_wiki_5 GitHub Wiki

Context

  • Summarizing the kinematic extraction, MGE fitting, and other input choices before running the first set of models.

High Level Information

  • Stats:

    • D = 75.5, Jensen+21 SBF
    • Brightest group member, 29 group members
    • _e = 247 km/s
    • Mk = −25.90, log10(Mstar) = 11.86
    • eps = 0.25, lambda = 0.03 (Slow rotator)
    • log10(Mhalo) = 14.4
    • MASSIVE XIV: Note that MASSIVE XIV says this galaxy is aligned, whereas MASSIVE X says it is misaligned
      • PA GMOS = 311 +/- 22
      • PA Mitchell = 161.0 ± 18.5
      • PA Phot = 34.9 (MASSIEV XIV: 35.8 ± 0.9)
      • Psi GMOS = 4.8 ± 9.3 (MASSIVE XIV)
      • Psi Mitchell = 53.9 (MASSIVE X)
      • Kinematic Twist, with note about low S/N
    • Photometry:
      • has disky isophotes between ∼10–30 arcsec. Outside of this range the isophotes go back to B4≈0. The ellipticity also peaks at ∼30 arcsec, beyond which it is approximately constant

  • GMOS:

    • All files seem to be on BOX.
    • Has both symmetric and non-symmetric binning schemes.
    • Other details:
      • Files on BOX say:
        • +X = 148 deg. E of N. in orbit modeling files
        • +X = +X = -32 deg. E of N. in paperversion figures of kinematic moments
      • PSF listed as: avg. FWHM = 0.39" (Gaussian)
      • 128 bins (symmetric)
      • 131 bins (non-symmetric)

  • Mitchell

    • All Mitchell files for all galaxies are present! There are folded and unfolded versions for everything. I'm currently sticking with our usual unfolded data.
    • 64 bins (unfolded)

Kinematics

Reproducing Irina's Results

Plot

  • Need to do this still for comparison to Irina.

  • Irina lists the settings on BOX as:

Nominal settings for pPXF fits to NGC 410 GMOS data
---------------------------------------------------

Wavelength range (rest):  8420-8730 A

Template list:  barthlist.txt
(15 stars from Barth et al. 2002 -- spectra are from MILES Calcium Triplet library)

Spectrum mask:  /GMOS/galaxy_files/N410/spectra/N410_spectra_mask1_nochipgap.fits
(bright sky lines as assessed from inspection of GMOS spectra)
(GMOS chip gap lies outside the pPXF wavelength range)

Initial velocity shift: 5000 km/s

Multiplicative polynomial: degree 3

Additive polynomial: degree 0  (constant offset)

Bias parameter: 0
And the spectra and resulting fits look beautiful!
Plot

Extending to h8 with symmetric binning

  • Need to do this still for comparison to Irina.

Extending to h8 with non-symmetric binning

  • And lastly, here's what we'll likely use for the modeling. Note that this uses the non-symmetric binning scheme (so it's hard to compare directly to Irina's which is symmetrically binned).
Bin Centers Radial Moments GMOS Maps Mitchell Maps
Spectral fits out to h8 for the 167 asymmetric bins
Plot

Preliminary Final Set of Inputs

  • Note these are the point-symmetrized versions of the kinematics.
Sym vs. Non-Sym Comparison Radial Moments GMOS Maps Mitchell Maps

MGE Fitting

  • Started doing some preliminary MGE fitting and background information gathering for N410 while we try to track down the spectra.
  • I've done our usual MGE fitting (unregularized) for a few different bounds on sigma (0,1,2,3,4) pixels from the center:
    • I likely need to go back and regularize these given the distribution of some of the q' values (reaching sometimes ~0.4 or so, severely limit deprojection).
    • The MGEs look alright, but I think there is some work that can be done. There seems to be a lot of contaminating light in the original image (see below), but we already are applying a mask, so unsure what to do beyond that.
HST Image
Comparing Deprojected 3D
Plot 0p0 1p0 2p0 3p0 4p0
Large Scale Residuals
Zoomed in Residuals
Large Contours
Zoomed Contours
Major Axis vs. Points
Deprojected
Enclosed Light
And here are all of the components listed for each of these MGEs:
  • 0p0:
  I [Lsun/pc^2]        Sig [arcsec]         q            PA 
        9681.66      0.0324025              1              0
        9050.49       0.347318              1              0
        10309.2       0.633257              1              0
        7491.84        1.11551       0.962657              0
        377.607        1.62733       0.504445              0
        3267.61        2.08795        0.89042              0
         1648.3        3.80717       0.897147              0
        278.569        4.47822       0.379619              0
         322.17        7.58985              1              0
        352.533        7.92285       0.479959              0
        418.602         10.289       0.713776              0
        158.103        19.8302       0.542994              0
         162.11        24.8701       0.777531              0
        53.4346        55.4603       0.718483              0
  • 1p0:
  I [Lsun/pc^2]        Sig [arcsec]         q            PA 
        3011.96            0.1              1              0
        9465.52       0.367336              1              0
        9797.26       0.655894              1              0
         7346.5         1.1368       0.948555              0
        3218.46        2.07872       0.899578              0
        197.405        2.10612        0.37637              0
        1596.32        3.84006       0.910268              0
        360.677        4.20741       0.442922              0
        337.767        8.11675       0.460665              0
        664.404        8.73971       0.855008              0
        133.828         15.954       0.529093              0
        111.016        21.5333       0.589591              0
        141.304        25.2786       0.789614              0
        52.1756        55.9714       0.715807              0
  • 2p0:
  I [Lsun/pc^2]        Sig [arcsec]         q            PA 
        7085.04        0.26999              1              0
        12453.4       0.560416              1              0
        8966.47        1.09309       0.947945              0
        3579.17        2.19181        0.84354              0
        1639.66        4.01745       0.887016              0
        322.078        7.88833       0.386095              0
        729.004         9.0452       0.814569              0
        171.004        19.0639       0.530437              0
         172.22        24.6405       0.770349              0
         54.449        55.0486       0.720493              0
  • 3p0:
  I [Lsun/pc^2]        Sig [arcsec]         q            PA 
        8270.01            0.3              1              0
        11004.6       0.582731              1              0
        7724.13        1.03596       0.979062              0
         807.51        1.44056       0.649154              0
        3436.59        1.92622       0.913712              0
         847.36        3.34952       0.707696              0
        1330.33        3.86266       0.937046              0
        103.423        5.23177       0.232302              0
         369.15        7.93331       0.472814              0
        651.932         8.7723       0.855676              0
        136.689        15.9643       0.527817              0
        110.333        21.6012       0.591412              0
        140.276        25.3063       0.790162              0
         52.079        56.0097       0.715665              0
  • 4p0:
  I [Lsun/pc^2]        Sig [arcsec]         q            PA 
          14032            0.4              1              0
        11813.7       0.887422       0.967784              0
        4718.84        1.81146       0.897153              0
        1939.57        3.63963       0.891352              0
         280.32        4.36675       0.354459              0
        356.382        7.48161              1              0
        375.049        7.86277       0.489978              0
          400.4        10.4518       0.705815              0
        156.468        19.9151       0.544937              0
         161.12        24.8729       0.778398              0
         53.435        55.4645       0.718339              0

Input Choices

  • Here are all the input files used in the modeling and the various choices I've made:
aperture_gmos/aperture_mitchell.dat
  • The only real number that I control in these files are the kin*gridtheta values, which represent the difference between the MGE +x axis and the +x axis of the kinematics.
  • These numbers are computed as PA_MGE - PA_kinematics:
    • The MGE PA comes from the fits abiove, which uses 35.75 E of N
    • The kinematic PAs come from the files on BOX.
      • For GMOS, these values are listed in /MASSIVE/GMOS_kinematics/N0410/N0410_PA_note.txt, which lists the +x axis of the orbit modeling files as +35 deg E of N.
      • For Mitchell, the numbers are located in /MASSIVE/Mitchell_reduced_data/NGC0547/kinematics_paperversion/unfolded/more_files/NGC0777-s2-unfolded-bininfo.txt, which lists the gal pa as +34.9 deg E of N.
      • This means that the input files should have kin1gridtheta and kin2gridtheta as 0.75 and 0.85, respectively.
bins_gmos/bins_mitchell.dat
  • I don't make any choices for these files. These are generated when preparing the orbit model files and associated the apertures with specific kinematics.
kindata.dat
  • This is our normal kindata.dat file, which comes from the Gh8Mh6 minimization described above.
  • This file contains 12 moments total, 8 of which are measured for the GMOS data and are from my own spectral fits, and 6 moments are fit for Mitchell (which I obtained from the files on BOX).
    • I've currently left in all bins, but we can decide to mask the overlapping Mitchell bins if we'd like.
    • I've added dummy moments for GMOS h9 through h12 and Mitchell h7 through h12 by assigning all the values to 0.0, and uncertainties identical to the h7/h8 errors for GMOS and h5/h6 errors for Mitchell.
mge.dat
  • Currently a placeholder MGE so that we can discuss the fits above.
orbitstart_loop/box.dat
  • These are just auxillary files which point to other files in our models. No changes or choices.
orblib_loop/box.in

  • The only specific changes I had to make here relate to the PSFs we assumed for each dataset:

    • I obtained the PSF for GMOS from the file /MASSIVE/GMOS_kinematics/N0410/N410_GMOS_PSF.txt, which says that the average weighted PSF = 0.48". The code accepts the equivalent "sigma" for this PSF, which I obtain as sigma = 0.48/(2sqrt(2ln(2))) = 0.204.
    • I am using the same PSF for Mitchell as we have used for the other galaxies, which has a sigma = 0.5, which is equivalent to a 1.177" PSF.

  • Other notable parameters that I kept constant from past orbit modeling:

    • innermost mass bin starts at 0.1"
    • outermost mass bin starts at 100" and extends to 316"
    • 200 dynamical time integration for the box library; 2000 dynamical times for the loop library
parameters_loop/box.in
  • Many of these parameters are changed on the fly (mass components and angles).
  • I'm using:
    • nEner = 40 for both
    • nI2 = 18, 9 for the loop and box library respectively
    • nI3 = 9, 9 for both
    • logrmin = 10^-2 arcsec (minimum IC radius in arcsec)
    • logrmin = 10^2.5 arcsec (maximum IC radius in arcsec)
    • interplogmin = 10^-8
    • interplogmax = 10^2.6
    • gnfw halo with gamma = 0
triaxmass.in
  • Auxillary file which points to other files.
triaxmassbin.in
  • Auxillary file which points to other files.
  • Also contains the same PSF information present in orblib_loop/box.dat.
triaxnnls.in
  • Mostly points to other files.
  • I've ensured that this file correctly assumes we have 12 total moments.
  • I've also checked to make sure that this file reads in both the loop orbit library and box orbit library.
⚠️ **GitHub.com Fallback** ⚠️