VOSA. Help and Documentation

Version 6.0, Apr 2018

1. Introduction
2. Input files
2.1. Upload files
2.2. VOSA file format
2.3. Single object
2.4. Manage files
2.5. Archiving
2.6. Filters
3. Objects
3.1. Coordinates
3.2. Distances
3.3. Extinction
4. Build SEDs
4.1. VO photometry
4.2. SED
4.3. Excess
5. Analysis
5.1. Model Fit
5.2. Bayes analysis
5.3. Template Fit
5.4. Templates Bayes
5.5. HR diagram
5.5.1. Interpolation
5.5.2. Flags
5.5.3. Example
5.6. Upper Limits
5.7. Statistics
6. Save results
6.1. Download
6.2. SAMP
6.3. References
6.4. Log file
6.5. Plots
7. VOSA Architecture
8. Phys. Constants
9. FAQ
10. Use Case
11. Quality
11.1. Stellar libraries
11.2. VO photometry
12. Credits
13. Helpdesk
14. About

HR diagram

VOSA offers the option to estimate values for the age and the mass of the objects. In order to do that, the (Teff,log(L)) values obtained from the chi-square fit are used as starting points for interpolating collections of theoretical isochrones and evolutionary tracks obtained from the VO. Then, a HR diagram is displayed showing the data points, isochrones and evolutionary tracks.

For each object, only the theoretical isochrones and evolutionary tracks more adequate to the model that best fits the observed photometry are used in the process. For instance, in the case where this model is "Kurucz" the Siess isochrones are used.

In the case that several collections are used (because we use one for some objects and another one for other objects) a HR plot will be generated for each collection, showing the isochrones, tracks and the points corresponding to the objects analysed using that collection.

You can play with the plots, decide to plot more or less information, locate the objects in it, etc.

Error estimation

In order to make an error estimation, the errors coming from the chi-square fit for Teff and LogL are used to generate a small grid with 9 points.

For each of these 9 points we make the interpolation as explained below.

The final values for (Age,Mass) will be the ones obtained for the point (Teff,LogL). But in some cases, the interpolated value of Age or Mass is different for some of the other 8 points. Thus, in the results table we show the minimum and maximum value obtained for each parameter when using any of the 9 points in this small grid.


(Below, all the explanations are given for the case of obtaining an estimation of the object age interpolating on isochrones. Everything is valid also for the case of obtaining an estimation of the mass interpolating on evolutionary tracks.)

The interpolation between isochrones involves to find the two closer isochrones to the (Teff,log(L)) point (one to each side of the point), calculate the distance from the point to each of the curves and, then evaluate a weighted average between the values of t for each isochrone.

Fig 1   Fig 2   $t=\frac{t_2 D_1+t_1 D_2}{D_1+D_2}$

In order to do this it is necessary to design an algorithm able to estimate the distance from a point to a curve defined by discrete points (note that we do not have an analytical curve but just a series of points that are assumed to define a curve).

Distance from a point to a curve

1.-The main method that we use to estimate the distance form the point to an isochrone is as follows:

  • First, given the point P (Teff,log(L)) find the closest point in the curve (point P1 at a distance D1)
  • Second, find the projection of P to the line defined by P1 and either the next point in the curve or the prior one.
  • If either of the projections lies into the interval defined by the two points (P1 and P2), then calculate the distance from the point P to the projection.
  • This value will be taken as the distance from P to the isochrone.
Fig 2    

2.-In some cases, it is not possible to use the above method because none of the proyections lie inside the interval between the two points that define the line.

Fig 2    

When that is the case, we can estimate the distance to the curve as the distance D1 from P to the closest point in the curve P1

Note that we consider this a worse approximation in general. Actually it is highly probable to be bad when P1 is the fist or last point in the curve.

Fig 2  

That is why we this method will only be used if the first one fails and if the closest point P1 is not the fist or last point in the curve.

Interpolated value for the age

If we have been able to find a curve on each side of the point P and the distance from that point to each curve, we can use the inverse of the distance as weights: $$t=\frac{\frac{1}{D_1}t_1+\frac{1}{D_2}t_2}{\frac{1}{D_1}+\frac{1}{D_2}}=\frac{t_2 D_1+t_1 D_2}{D_1+D_2}$$

In some cases, we are able to determine only the distance to one curve, but we know that there exist an isochrone on each side of the point. If that happens we just show a range of values for the age using the ones corresponding to each isochrone as lower and upper limits.

Finally, if the point lies outside the area covered by the isochrones, we do not even try to estimate a value for the age or the mass of the object.

outside the area

Flags for interpolated values

Whenever we are not able to find a value for the age or the mass of an object or it has been determined using a worse approximation that the one that we consider the best (See above) a flag is shown right to the value.

These are the possible flags and their meanings:

[1]The distance to one of the closer curves has been estimated as the one to the closest point in the curve
[2]The distance to both the closest curves has been estimated as the one to the closest point in each curve
[3]Only a range of values can be estimated
[4]The point lies outside the area covered by the isochrones
[5]No estimation has been posible


We have made a chi-square model fit for a set of objects. For some of the objects the best fit model was Kurucz and for another ones it was BT-Settl-CIFIST. Thus, when we enter the "HR diagram" tab we see that two different collections of isochrones and tracks are going to be used: BHAC15 (for those objects fitted by the BT-Settl-CIFIST model) and Siess (for the ones fitted by Kurucz model).

We could play with the ranges of parameters, restricting the values of the age and mass to be considered in the analysis. But we prefer to keep the full range and click the "Make the HR Diagram" button.

The interpolation process, to obtain the best values (and ranges) for the age and mass of each object, is performed asynchronously so that you don't need to stay in front of the computer waiting for the results. You can close your browser and come back later. If the process is not finished, VOSA will give you some estimation of the status of the operation and the remaining time.

When the process is finished, you can see the list of objects with the interpolation results, and two HR plots, one for each collection of isochrones and tracks.

If you click in any graph, VOSA will locate the object closest to the click point and will show you its properties.

If, instead, you click on one object name in the list, VOSA will locate that object in the corresponding graph.

You also can play with the plots. There are options to zoom to the objects range or to the models range. Other options allow you to define the exact range of each coordinate. And you also can decide what isochrones or tracks you want to display.