MBCAA Observatory

V1412 Aql: a simulation

Michel Bonnardeau
First version 24 Aug 2009
Updated 2 Sep 2009 (more observations)
Updated and corrected 3 Sep 2009 (more observations; uncertainties on the ToM of the 2 eclipses taken into account at the suggestion of Wolfgang Renz; correction of the number of solutions)
Updated 5 Sep 2009 (more observations; simulation with the period down to 1 d; inclination)
Corrected 10 Sep 2009 (uncertainties on the ToM of the 2 eclipses)
Updated 15 Sep 2009 (more observations)
Updated 20 Sep 2009 (more observations, slightly different ranges for the simulation)
Updated 21 Sep 2009 (inclination taken into account in the Monte Carlo)
Updated 25 Sep 2009 (more observations, smaller mass range and cos(i) at the suggestion of Christian Knigge)
Updated 11 Oct 2009 (better timing for the 1985 eclipse)
Updated 21 Oct 2009 (more observations)
Updated 30 Oct 2009 (more observations)
Updated 28 Jan, 14 Apr, 30 Jun (more observations), 4 Jul 2010 (eclipse duration vs period)
Updated 8, 16 Jul, 2 Aug 2010 (more observations)
Updated 24 Sep 2010 (more observations), 6 Jan 2011 (ephemeris for 2011)


A Monte Carlo simulation which may be useful to search for eclipses is proposed.


The white dwarf G24-9 (or V1412 Aql) was observed to be unexpectedly faint on 2 occasions, in 1985 and 1988. This is interpreted as eclipses due to a dark, substellar companion (Zuckerman & Becklin (1988)).

According to Bergeron et al (2001) G24-9 has a parallax of 0.0447" (then a distance of 22.4 pc), a temperature of 6400K, a mass M=0.65+/-0.04 solar mass, a gravity g=10^8.13cm/s2 (then a radius sqrt(G*M/g)=8018km).

Feb 2009, the AAVSO launched a campaign to observe this object so as to detect its eclipses (AAVSO Special Notice #148).

I propose here a Monte Carlo simulation which may be useful to speed up the discovery of eclipses.


The 2 observed eclipses are:

19851007.11 (Landolt (1985)), t85=2,446,345.61828 HJD (C. Knigge, personal communication). The uncertainty is considered as negligeable et85=0;

19880715.3 (Carilli & Conner (1988)), heliocentric correction 441.74s, t88=2,447,357.805 HJD. The uncertainty is taken as et88=0.05 HJD.

The AAVSO has 2509 negative (i.e. no eclipse) observations on 24 Sep 2010.

Arne Henden obtained 14 time-series showing no eclipses, from 30 Oct to 8 Dec 2008, total duration 34.1 h. And also 7 short time-series, in 2002 and 2008, total duration 2.6 h.

Robert Fridrich obtained a time-series on 1 Sep 2009, duration 0.3 h.

Christian Knigge obtained 478 negative measurements in 2009.

M. Wardak obtained 325 negative measurements in 2009 (6 time-series, total duration 17.8 h).


The orbital period is P=(t88-t85)/n=(1012.187 days)/n where n is an integer (assuming it is constant, i.e. no mass transfer, no third body).

G29-4 is a white dwarf, so it has a small size (about that of Earth), then the eclipse duration tau is given mostly by the diameter D of the occulting body:

with M the mass of the system, i the inclination, G the gravitational constant (taking for the eccentricity e=0).

The computer simulation is a Monte Carlo one where a large number of random sets of n, M, D, i are used to derived ephemeris. The ephemeris that are retained are those that give the 2 observed eclipses and that do not give eclipses for the negative observations (with i large enough to give an eclipse).

The algorithm works the following way:

  • 100,000,000 random sets of n, M, D, i are generated, with n between 1 and 1000, M between 0.57 and 0.80 solar masses (2 sigma uncertainty for the white dwarf, and up to 0.07 for the occulter (a brown dwarf)), D between 0.15 and 5 jovian diameters,cos(i) between cos(85°) and 0;
  • when i is too small for an eclipse to occur, the set is rejected;
  • for each set the period P0=(t88-t85)/n and the eclipse duration tau is computed;
  • the ephemeris is HJD(E)=T+P*E with T an random number between t85-tau and t85+tau, and P a random number between P0-tau/n and P0+tau/n;
  • the ephemeris than do not give the eclipses at t85 and t88 are rejected;
  • the ephemeris that give an eclipse for one of the negative observations or for one of the 14+7+1 time-series are rejected.
  • The spectra of n, M, D and i solutions are:

    The "probability" is actually the number of acceptable solutions from the simulation.

    A close-up view. n>500 (or P<2 days) is very unlikely. P<1 day is ruled out.

    The average for the inclinations is 89.49 .

    The probability for the eclipse duration:

    The average for the eclipse durations is 45.0 mn.

    and the period:

    The average for the periods is 26.4 days.

    The eclipse durations versus the periods:

    The probability for future eclipses may be computed:

    There are 2 larges peaks in 15 September 2010 and 25 June 2013 because, whatever n is, there should be eclipses at these dates.

    However, the uncertainties et85 and et88 introduce an uncertainty on their timing of (where t is the current time):
    dt = et85+(et85+et88)*(t-t85)/(t88-t85)
    dt = 0.5 day

    The 15 September 2010 predicted eclipse

    Whatever n is, there should have been an eclipse on 15 September 2010. This would be the 8th*n eclipse since the one observed in 2008.

    V1412 Aql was observed around that date by AAVSO observers in the USA and Europe, but no eclipse was found. May be the eclipse was missed. There was a lack of coverage from the Asian longitudes and the event duration is probably less than 1 hour. Actually, the probability that there was an eclipse is still high:

    Ephemeris for the next eclipses

    The best opportunities to detect an eclipse in 2011 are (with the uncertainty dt=0.5day):

    THJD2011 date
    587.4 01 25.9
    594.8 02 02.3
    629.8 03 09.3
    631.3 03 10.8
    657.7 04 06.2

    THJD2011 date
    657.7 04 06.2
    664.7 04 13.2
    675.3 04 23.8
    699.5 05 18.0
    719.3 06 06.8
    734.5 06 22.0

    THJD2011 date
    763.3 07 20.8
    769.3 07 26.8
    792.8 08 19.3
    804.2 08 30.7
    807.3 09 02.8
    839.2 10 04.7

    THJD2011 date
    839.2 10 04.7
    851.3 10 16.8
    860.1 10 25.6
    874.1 11 08.6
    895.3 11 29.8
    908.9 12 13.4


    Bergeron P., Leggett S.K., Ruiz M.T. (2001) ApJ Supp 133 413.

    Carilli C., Conner S. (1988) IAU Circ. 4648.

    Landolt A.U. (1985) IAU Circ. 4125.

    Zuckerman B., Becklin E. (1988) IAU Circ. 4652.

    Site map


    Copyright notice

    MBCAA Observatory. Copyright 2004-2011