search

nespinoza / juliet

A versatile modelling tool for transiting and non-transiting (single and multiple) exoplanetary systems
MIT License
53 stars 31 forks source link

canot fit secondary eclipses #44

Open MohammadEftekhar opened 3 years ago

MohammadEftekhar commented 3 years ago

For fitting secondary eclipse, How could I set parameter of 'r1_p1','r2_p1" ?

nespinoza commented 3 years ago

Hi @MohammadEftekhar,

Currently juliet does not directly fit for secondary eclipses, but you can of course deactivate the limb-darkening by setting q1 equal to zero to mimic a secondary eclipse. There, the time-of-transit would really be the time-of-secondary eclipse, of course.

Let me know if this helps. If you need further help, please provide me with a minimal working example so you can help me help you!

Néstor

MohammadEftekhar commented 3 years ago

Hi Néstor Thank You for response. Yes as You say I did that. What about q2? I setted q1 and q2 equal to zero. Thank You Mohammad

On Tue, Nov 10, 2020 at 8:15 AM Néstor Espinoza notifications@github.com wrote:

Hi @MohammadEftekhar https://github.com/MohammadEftekhar,

Currently juliet does not directly fit for secondary eclipses, but you can of course deactivate the limb-darkening by setting q1 equal to zero to mimic a secondary eclipse. There, the time-of-transit would really be the time-of-secondary eclipse, of course.

Let me know if this helps. If you need further help, please provide me with a minimal working example so you can help me help you!

Néstor

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/nespinoza/juliet/issues/44#issuecomment-724447331, or unsubscribe https://github.com/notifications/unsubscribe-auth/AP77OM2ROABX56B5NBJJ5YLSPDAO5ANCNFSM4TNDYYDA .

nespinoza commented 3 years ago

Hi Mohammad,

That would work as well. Remember that the limb-darkening law is not defined by q1 and q2, but by u1 and u2 such that, e.g., for the quadratic law, u1 = 2.*np.sqrt(q1)*q2 and u2 = np.sqrt(q1)*(1.-2.*q2). So if you set both q1 and q2 to zero, you get u1 = 0 and u2 = 0. Given the limb-darkening law is I(mu) = 1-u1*(1-mu)-u2*(1-mu)**2, this gives I(mu) = 1, i.e., a uniform intenstiy profile, which is what we wanted to emulate here.

Let me know if this works! Néstor

MohammadEftekhar commented 3 years ago

Thank You Néstor