Jason Eastman
@exofast.bsky.social
Follow me for updates to EXOFASTv2, an IDL code to fit exoplanet data.
http://arxiv.org/abs/1206.5798
http://arxiv.org/abs/1907.09480
http://arxiv.org/abs/1206.5798
http://arxiv.org/abs/1907.09480
A few things I think are obsolete and were probably never used have been removed. Most significantly, I removed REFINESTAR, which was not included in my testing, difficult to maintain, and I think is better done by the user when it's occasionally useful anyway.
June 21, 2024 at 5:39 PM
A few things I think are obsolete and were probably never used have been removed. Most significantly, I removed REFINESTAR, which was not included in my testing, difficult to maintain, and I think is better done by the user when it's occasionally useful anyway.
The inputs have been grouped and ordered more logically. Tutorials and example templates are probably still the easiest way to get started.
June 21, 2024 at 5:39 PM
The inputs have been grouped and ordered more logically. Tutorials and example templates are probably still the easiest way to get started.
The documentation in $EXOFAST_PATH/exofastv2.pro was very outdated, and occasionally in conflict with $EXOFAST_PATH/mkss.pro. Overlapping documentation in mkss.pro has been moved to exofastv2.pro, which has been updated.
June 21, 2024 at 5:39 PM
The documentation in $EXOFAST_PATH/exofastv2.pro was very outdated, and occasionally in conflict with $EXOFAST_PATH/mkss.pro. Overlapping documentation in mkss.pro has been moved to exofastv2.pro, which has been updated.
$EXOFAST_PATH/getdata.py (python!) is a new utility function that uses lightkurve to retrieve TESS, Kepler, or K2 data and format it for EXOFASTv2. It can optionally undeblend the TESS LCs so you can re-deblend them based on multi-star SED models.
June 21, 2024 at 5:38 PM
$EXOFAST_PATH/getdata.py (python!) is a new utility function that uses lightkurve to retrieve TESS, Kepler, or K2 data and format it for EXOFASTv2. It can optionally undeblend the TESS LCs so you can re-deblend them based on multi-star SED models.
For high SNR LCs with a known planetary eccentricity (where sigma_rho_LC ~< 5%), rstar is better constrained by the density (see Mahajan et al, 2024) rather than MANNRAD.
June 21, 2024 at 5:38 PM
For high SNR LCs with a known planetary eccentricity (where sigma_rho_LC ~< 5%), rstar is better constrained by the density (see Mahajan et al, 2024) rather than MANNRAD.
@amann.bsky.social 's MANNRAD and/or MANNMASS uses the Mann+ 2015 and Mann+ 2019 relations to constrain the radius and/or mass of M dwarfs, respectively, based on the apparent K-band magnitude (and derived absolute K-band magnitude).
June 21, 2024 at 5:38 PM
@amann.bsky.social 's MANNRAD and/or MANNMASS uses the Mann+ 2015 and Mann+ 2019 relations to constrain the radius and/or mass of M dwarfs, respectively, based on the apparent K-band magnitude (and derived absolute K-band magnitude).
We can also now model ramps (intended for JWST data). See /FITRAMP.
June 21, 2024 at 5:38 PM
We can also now model ramps (intended for JWST data). See /FITRAMP.
If the bandname can be interpreted as a number with "_" instead of "." (e..g, 8_0), it will be assumed to be the wavelength in microns, and will be plotted at the end. This, along with TDELTAVS, enables transmission spectroscopy modeling.
June 21, 2024 at 5:38 PM
If the bandname can be interpreted as a number with "_" instead of "." (e..g, 8_0), it will be assumed to be the wavelength in microns, and will be plotted at the end. This, along with TDELTAVS, enables transmission spectroscopy modeling.
NOCLARET is now an NTRANSITS array, and the transit bandname no longer needs to be among the allowed names (if CLARET is disabled for the corresponding transit).
June 21, 2024 at 5:38 PM
NOCLARET is now an NTRANSITS array, and the transit bandname no longer needs to be among the allowed names (if CLARET is disabled for the corresponding transit).
Note to the note: if correlations between linked parameters are not perfectly linear, extremely tight bounds can dramatically slow convergence times.
June 21, 2024 at 5:37 PM
Note to the note: if correlations between linked parameters are not perfectly linear, extremely tight bounds can dramatically slow convergence times.
Note: only fitted parameters can be fixed and you cannot do math on them or specify prior widths or bounds based on other parameters (yet). If you want to link derived parameters, specify a Gaussian width that is small relative to the error.
June 21, 2024 at 5:37 PM
Note: only fitted parameters can be fixed and you cannot do math on them or specify prior widths or bounds based on other parameters (yet). If you want to link derived parameters, specify a Gaussian width that is small relative to the error.
You could also model EBs by linking the "planet" mass and radius to the stellar mass and radius.
Arielle Frommer has a paper coming out soon with more details about that.
Arielle Frommer has a paper coming out soon with more details about that.
June 21, 2024 at 5:37 PM
You could also model EBs by linking the "planet" mass and radius to the stellar mass and radius.
Arielle Frommer has a paper coming out soon with more details about that.
Arielle Frommer has a paper coming out soon with more details about that.
This could lead to significantly better age constraints on all stars (and maybe even detailed cluster modeling). An isochrone is created for each subset of stars with the same age and initial metallicity, with that subset of stars plotted on top.
June 21, 2024 at 5:37 PM
This could lead to significantly better age constraints on all stars (and maybe even detailed cluster modeling). An isochrone is created for each subset of stars with the same age and initial metallicity, with that subset of stars plotted on top.
You can link parameters together in the prior file. If you're modeling a bound binary, you'll want to link their initial metalicities, distances, extinctions, and ages:
initfeh_1 initfeh_0 0
distance_1 distance_0 0
av_1 av_0 0
age_1 age_0 0
initfeh_1 initfeh_0 0
distance_1 distance_0 0
av_1 av_0 0
age_1 age_0 0
June 21, 2024 at 5:37 PM
You can link parameters together in the prior file. If you're modeling a bound binary, you'll want to link their initial metalicities, distances, extinctions, and ages:
initfeh_1 initfeh_0 0
distance_1 distance_0 0
av_1 av_0 0
age_1 age_0 0
initfeh_1 initfeh_0 0
distance_1 distance_0 0
av_1 av_0 0
age_1 age_0 0
The MISTSEDFILE and SEDFILE models now plot the atmosphere and residuals (inspired by Sam Yee). Note: the plotted MISTSEDFILE atmospheres come from the NEXTGEN models, not the MIST BC's C3K models. This inconsistency is purely aesthetic.
June 21, 2024 at 5:36 PM
The MISTSEDFILE and SEDFILE models now plot the atmosphere and residuals (inspired by Sam Yee). Note: the plotted MISTSEDFILE atmospheres come from the NEXTGEN models, not the MIST BC's C3K models. This inconsistency is purely aesthetic.
Note: the inferior "FLUXFILE=" SED model which uses a tophat for each filter response does not support multiple SEDs or differential photometry. I expect the FLUXFILE SED model development to be frozen, as it's not much slower and far less accurate than the new one.
June 21, 2024 at 5:36 PM
Note: the inferior "FLUXFILE=" SED model which uses a tophat for each filter response does not support multiple SEDs or differential photometry. I expect the FLUXFILE SED model development to be frozen, as it's not much slower and far less accurate than the new one.
Thanks to Sam Yee for EXOFASTv2's first substantive pull request that backports the SEDFILE model based differential photometry to the "MISTSEDFILE=" SED model!
June 21, 2024 at 5:36 PM
Thanks to Sam Yee for EXOFASTv2's first substantive pull request that backports the SEDFILE model based differential photometry to the "MISTSEDFILE=" SED model!
A separate error scaling term is fit for each spectrophotometry data set, and every fit I've seen has large errors. It's unclear to me whether this is due to errors in the model atmospheres or in the absolute flux calibration of the spectrum. Probably both.
June 21, 2024 at 5:36 PM
A separate error scaling term is fit for each spectrophotometry data set, and every fit I've seen has large errors. It's unclear to me whether this is due to errors in the model atmospheres or in the absolute flux calibration of the spectrum. Probably both.
This new model supports an arbitrary number of spectrophotometry files (SPECPHOTPATH=). MKTICSED automatically grabs Gaia spectrophotometry, but you can also supply your own (e.g., from JWST).
June 21, 2024 at 5:36 PM
This new model supports an arbitrary number of spectrophotometry files (SPECPHOTPATH=). MKTICSED automatically grabs Gaia spectrophotometry, but you can also supply your own (e.g., from JWST).
There is a new (third) SED model that uses NEXTGEN atmospheres and detailed filter transmission curves (which can be supplied by the user). Specify SEDFILE= (instead of MISTSEDFILE= or FLUXFILE=) to use this model. This is great for AO bands not supported by MIST BC tables.
June 21, 2024 at 5:36 PM
There is a new (third) SED model that uses NEXTGEN atmospheres and detailed filter transmission curves (which can be supplied by the user). Specify SEDFILE= (instead of MISTSEDFILE= or FLUXFILE=) to use this model. This is great for AO bands not supported by MIST BC tables.
The SED can be used to automatically deblend a transit lightcurve (see SEDDEBLEND).
June 21, 2024 at 5:35 PM
The SED can be used to automatically deblend a transit lightcurve (see SEDDEBLEND).