[PIPE2D-511] Parameter interpolation of model templates with a single 4D RBF model Created: 29/Jan/20  Updated: 20/Dec/21  Resolved: 20/Dec/21

Status: Done
Project: DRP 2-D Pipeline
Component/s: None
Affects Version/s: None
Fix Version/s: None

Type: Story Priority: Normal
Reporter: Takuji Yamashita Assignee: Takuji Yamashita
Resolution: Done Votes: 0
Labels: flux-calibration, model-templates
Remaining Estimate: Not Specified
Time Spent: Not Specified
Original Estimate: Not Specified

Attachments: PNG File evaluation_single4drbf_p7000_g+4.0_m0.0_t01_z+1.00_a+0.20.png     PNG File RMS_single4DRBF_2020March30.png     PNG File RMS_single4DRBF_2020March30_teff7000_logg4.png    
Issue Links:
Blocks
blocks PIPE2D-365 Interpolate between model templates Done
is blocked by PIPE2D-781 Organize the drp_fstar repository Done
Story Points: 4
Epic Link: flux calibration
Sprint: 2DDRP-2020 A
Reviewers: hassan

 Description   

Attempt to make parameter interpolations of stellar model templates using a single 4D radial basis function (RBF) which is modeled over all the original model templates. 

The 4D RBF interpolations using neighboring model templates reproduce model spectra with high quality (<1%) as reported in PIPE2D-460 but the processes are time-consuming because RBF models are generated at every parameter grid. In this ticket, we try to reduce the processing time by using only a single RBF. 



 Comments   
Comment by Takuji Yamashita [ 14/Feb/20 ]

The interpolation using a single 4D RBF has been tested at an existing parameter point. The parameters are Teff, log(g), Z, and alpha elements. The result is as follows:

The processing time to make an RBF model: ~ 12 hours
The accuracy of a produced spectrum: 0.02 % (for a test sample)

Once the RBF model is generated, we can get an interpolated spectrum immediately. The accuracy is good enough.

Comment by Takuji Yamashita [ 14/Feb/20 ]

I will check the accuracies at several parameter points as well.

Comment by Takuji Yamashita [ 07/Apr/20 ]

The accuracies of this interpolation has been evaluated at 192 existing parameter points. An AMBRE original spectrum is removed when an interpolation is made there. The parameter points for this evaluation were selected to cover the whole grids. To save the processing time, the interpolations were made at 380-500nm, where the accuracies are expected to be worse than redder wavelengths. 

The RMS (%) is r.m.s of residuals between an AMBRE spectrum and a generated interpolation. The spectra are convolved with a quasi LSF (a Gaussian with 0.2 nm FWHM).  

Almost all the interpolations were generated in high accuracies (less than 1%), but 12% show >1% RMS. 1% show >10% RMS. I cannot any clear correlation of RMS to any parameters. When Teff and log(g) are fixed, at a specific metallicity interpolation spectra systematically show worse RMS. The reason on this and treatment for these low-accuracy interpolations are under discussion. 

 

RMS

 

RMS at Teff=7000 and log(g) = 4.0

 

Comment by Takuji Yamashita [ 20/Dec/21 ]

We can close this ticket because we have developed the parameter interpolation and we have evaluated this method in PIPE2D-460. The code implementation should be done in PIPE2D-781.

Generated at Sat Feb 10 15:54:16 JST 2024 using Jira 8.3.4#803005-sha1:1f96e09b3c60279a408a2ae47be3c745f571388b.