REFLUX - SIM Binary Simulator


REFLUX is a program that utilizes the Eclipsing Light Curve (ELC) code (Orosz, J.A., & Hauschildt, P.H. 2000, A&A, 364, 265) to model the flux-weighted astrometric reflex motion of a binary system with the specified parameters. Input the parameters of your system below, hit submit, and the program will generate plots of the center of light reflex orbital motion, x and y componenets versus orbital phase, the spectral energy distibution for each object and the system, as well as 3D model animation of your system, (examples on right for SS Cyg). A file that contains descriptions for all the files is also included.

ELC's convention is that star 1 is closest to the observer at phase 0.0. Thus, for sufficiently high inclinations, star 2 would be eclipsed at phase 0.0. Star 2 is typically the higher mass object. If a disk is present, it goes around star 2. Thus in CV systems, star 1 is usually the red dwarf, and star 2 the white dwarf / black hole / compact object.

If you have any questions, run into any problems, or have suggestions for improvement, please contact me at jlcough@nmsu.edu. This seems to work best with FireFox web browser. If you use Safari or IE, it is possible that you will not get a % complete dialog, however, the program will still work, and the results will still pop up eventually - just be patient, or use Firefox. Also note that model atmospheres are only available for main-sequence and giant stars - do not use on disks or compact objects or the code may choke. (In this case, or others such as accidently setting the inner disk radius bigger than the outer, the program will still run, but your resutls will all be values of 0 or nan.)

The default values below are for the SS Cyg system.



MAIN SYSTEM PARAMETERS


Name of System - No spaces! Will be used to name output files and plots
Distance to the system in parsecs
Inclination of the system in degrees
Period of the system in days
Effective Temperature of Star 1 in Kelvin
Effective Temperature of Star 2 in Kelvin (Set this to a negatve number to make star 2 invisible, as in the case of a black hole.)
No   Yes   ⇐ Use a disk around star 2?
Blackbody   Free-free   ⇐ Type of disk? If free-free, temperature power-law exponent is automatically set to 0.0.
No   Yes for all   Yes only for Star 1   ⇐ Use model atmospheres?
Phase resolution in degrees (i.e. 360 divided by this value = number points that will be generated.)


STELLAR MASSES


Mass ratio (M2/M1)
Semi-major axis of the system in Solar Radii
             OR
The mass of Star 1 in solar masses - set both this and mass of Star 2 to greater than zero to use this instead of the mass ratio and semi-major axis.
The mass of Star 2 in solar masses - set both this and mass of Star 1 to greater than zero to use this instead of the mass ratio and semi-major axis.


STELLAR SIZES


Fill factor of Star 1 (Must be ≤ 1.0)
Fill factor of Star 2 (Must be ≤ 1.0)
             OR
The effective radius (Eggleton) of Star 1 in solar radii - Set to greater than zero to use this instead of the fill factor
The effective radius (Eggleton) of Star 2 in solar radii - Set to greater than zero to use this instead of the fill factor


DISK PARAMETERS


Opening angle of the disk in degrees

Inner radius of the disk in units of Star 2's roche lobe
Outer radius of the disk in units of Star 2's roche lobe
OR
The inner disk radius in units of solar radii - Set to greater than zero to use instead of the units of Star 2's roche lobe used above.
The outer disk radius in units of solar radii - Set to greater than zero to use instead of the units of Star 2's roche lobe used above.

Temperature of the innermost part of the disk in Kelvin
Power law exponent of disk temperature where T_disk(R) = T_inner*(R/r_inner)^x, where R is the distance from Star 2, r_inner is the inner disk distance from Star 2, and x is the specfied exponent.
Value to multiply the disk flux by for SED fitting purposes.


SPOT PARAMETERS


STAR 1
Temperature factor of the spot on Star 1 - Set to a negative number to turn off spot.
Latitude of the spot on Star 1 in degrees. (0° = north pole, 90° = equator, 180° = south pole.) - Set to a negative number to turn off spot.
Longitdue of the spot on Star 1 in degrees. (0° = Faces inner Lagrangian point, 90° = Side of star seen at phase 90°, 180° = Back of star facing away from other star. 270° = Side of star seen at phase 270°.) - Set to a negative number to turn off spot.
Angular size of the spot on Star 1. (90° degrees covers exactly half the star.) - Set to a negative number to turn off spot.

STAR 2
Temperature factor of the spot on Star 2 - Set to a negative number to turn off spot.
Latitude of the spot on Star 2 in degrees. (0° = north pole, 90° = equator, 180° = south pole.) - Set to a negative number to turn off spot
Longitdue of the spot on Star 2 in degrees. (0° = Faces inner Lagrangian point, 90° = Side of star seen at phase 270°, 180° = Back of star facing away from other star. 270° = Side of star seen at phase 90°.) - Set to a negative number to turn off spot.
Angular size of the spot on Star 2. (90° degrees covers exactly half the star.) - Set to a negative number to turn off spot.

DISK
Temperature factor of the spot on the Disk - Set to a negative number to turn off spot.
Azimuth of the spot on the Disk in degrees. (0° = seen directly at phase 0°, 90° = seen directly at phase 90°, etc.) - Set to a negative number to turn off spot.
Radial cutoff of the spot on the Disk. (0.9 = The spot fills the outer 10% of the disk sector, ect.) - Set to a negative number to turn off spot.
Angular size of the spot on the Disk. - Set to a negative number to turn off spot.