Celestia representation for the binary star+planet in
Ortiz et al. (arXiv:1608.00963 [astro-ph.EP]). Note this planet is supplied in extrasolar.ssc for a single star representation, so the relevant sections must be commented out of that file first otherwise you'll end up with one version of the planet orbiting the barycentre, the other orbiting the primary star.
The secondary could be a main sequence star or a white dwarf, here I've gone with the main sequence star version. I've used the orbital inclination and ascending node given in the figure 8. I estimated the stellar properties for star B using
isochrones and the radius of the planet is estimated using
forecaster.
Version for Celestia with units support/Celestia.Sci:.stc file
Code:
Replace Barycenter 36616 "HIP 36616 AB:HD 59686 AB"
{
RA<deg> 112.951650
Dec<deg> 17.086047
Distance<pc> 96.90
}
"HIP 36616 A:HD 59686 A"
{
AppMag 5.45
Radius<rS> 13.2
Temperature 4658
SpectralType "K2III"
OrbitBarycenter "HIP 36616 AB"
EllipticalOrbit {
Period<d> 11680
SemiMajorAxis<AU> 3.6917
Eccentricity 0.729
ArgOfPericenter<deg> 284.07
Inclination<deg> 55.56
AscendingNode<deg> 235.09
MeanAnomaly<deg> 259
Epoch 2452621.84
}
}
"HIP 36616 B:HD 59686 B"
{
AbsMag 7.88 # from isochrones
Temperature 4250 # from isochrones
SpectralType "K6V" # closest Celestia type that matches temperature
OrbitBarycenter "HIP 36616 AB"
EllipticalOrbit {
Period<d> 11680
SemiMajorAxis<AU> 10.1473
Eccentricity 0.729
ArgOfPericenter<deg> 104.07
Inclination<deg> 55.56
AscendingNode<deg> 235.09
MeanAnomaly<deg> 259
Epoch 2452621.84
}
}
.ssc file
Code:
"b" "HD 59686 A"
{
Radius<rJ> 1.11 # from Forecaster
Mass 2871
Texture "exo-class3.*"
EllipticalOrbit {
Period<d> 299.36
SemiMajorAxis<AU> 1.0860
Eccentricity 0.05
ArgOfPericenter<deg> 255.67
Inclination<deg> 55.56 # to match stellar binary
AscendingNode<deg> 235.09 # to match stellar binary
MeanAnomaly<deg> 301
Epoch 2452621.84
}
UniformRotation {
Period<h> 10
Inclination<deg> 55.56 # to match orbit
AscendingNode<deg> 235.09 # to match orbit
}
}
Version for Celestia without units support.stc file
Code:
"HIP 36616 A:HD 59686 A"
{
AppMag 5.45
Radius 9187200
SpectralType "K2III"
OrbitBarycenter "HIP 36616 AB"
EllipticalOrbit {
Period 31.978
SemiMajorAxis 3.6917
Eccentricity 0.729
ArgOfPericenter 284.07
Inclination 55.56
AscendingNode 235.09
MeanAnomaly 259
Epoch 2452621.84
}
}
"HIP 36616 B:HD 59686 B"
{
AbsMag 7.88 # from isochrones
SpectralType "K6V" # closest Celestia type that matches isochrones temperature
OrbitBarycenter "HIP 36616 AB"
EllipticalOrbit {
Period 31.978
SemiMajorAxis 10.1473
Eccentricity 0.729
ArgOfPericenter 104.07
Inclination 55.56
AscendingNode 235.09
MeanAnomaly 259
Epoch 2452621.84
}
}
.ssc file
Code:
"b" "HD 59686 A"
{
Radius 79356 # from Forecaster
Mass 2871
Texture "exo-class3.*"
EllipticalOrbit {
Period 0.81960
SemiMajorAxis 1.0860
Eccentricity 0.05
ArgOfPericenter 255.67
Inclination 55.56 # to match stellar binary
AscendingNode 235.09 # to match stellar binary
MeanAnomaly 301
Epoch 2452621.84
}
UniformRotation {
Period 10
Inclination 55.56 # to match orbit
AscendingNode 235.09 # to match orbit
}
}