[mesa-users] MESA3 paper 7_binary_with_rot

Tue Jun 16 04:03:55 EDT 2015

Cool,
Very useful stuff, Roni.
Let us know how it goes!

-M

Matteo Cantiello  |  http://matteocantiello.com/ 
Associate Specialist in Astrophysics  | Kavli Institute for Theoretical Physics
Chief Scientist and Board Member, Authorea  |  https://www.authorea.com/

On Jun 16, 2015, at 9:38 AM, Roni Waldman <roni181066 at gmail.com> wrote:

> Sure, and I will.
> 
> ‫בתאריך יום ג׳, 16 ביוני 2015 ב-10:28 מאת ‪Pablo Marchant‬‏ <‪pamarca at gmail.com‬‏>:‬
> No, but the use_other_tsync option is there for you to play with it ;)
> 
> On Tue, Jun 16, 2015 at 9:10 AM, Roni Waldman <roni181066 at gmail.com> wrote:
> Thanks Pablo.
> Do you happen to have a subroutine for the Tassoul formula?
> I'm trying to compare to Yoon et al 2010.
> 
> Roni
> 
> ‫בתאריך יום ג׳, 16 ביוני 2015 ב-10:04 מאת ‪Pablo Marchant‬‏ <‪pamarca at gmail.com‬‏>:‬
> My bad, here it is. Will update the ones on mesastar
> 
> On Tue, Jun 16, 2015 at 8:37 AM, Roni Waldman <roni181066 at gmail.com> wrote:
> Hi,
> More on this one.
> The inlist_project has:
>    ! use_other_tsync will use the Zahn timescale as defined in run_binary_extras.f
>    !use_other_tsync = .true.
> 
> I dont't see any other_tsync in the default run_binary_extras.f, and there's no special run_binary_extras.f in the 7_binary_with_rot directory.
> 
> Roni
> 
> 
> ‫בתאריך יום ב׳, 15 ביוני 2015 ב-21:11 מאת ‪Roni Waldman‬‏ <‪roni181066 at gmail.com‬‏>:‬
> Same here.
> 
> ‫בתאריך יום ב׳, 15 ביוני 2015 ב-19:37 מאת ‪Willie Strickland‬‏ <‪cwskas at gmail.com‬‏>:‬
> Thanks for the confirmation.  I thought they were probably for debugging, but didnt see obvious way to turn them off.
> 
> Willie
> 
> 
>> On Jun 15, 2015, at 11:13 AM, Pablo Marchant <pamarca at gmail.com> wrote:
>> 
>> The w > w_crit messages come from the implicit method that increases the wind strength so the star does not evolve avobe critical rotation. Maybe they could be turned off, I mainly used that info for debugging. And your results are the same as mine.
>> 
>> On Mon, Jun 15, 2015 at 5:59 PM, Willie Strickland <cwskas at gmail.com> wrote:
>> Did you get a lot of w > w_crit lines in the rotating one?  I got the same results as you it sounds like, but had all these redo lines.
>> 
>> 8,821 lines	fix w > w_crit: bisect mdots and redo
>>      89 lines	w > w_crit: revise mdot and redo
>> 
>> My run finished at step 6429, age 1.254e7, 2.45 mass1, 18.7 mass2.
>> 
>> accretor photosphere has reached its Roche Lobe
>>                                      accretor (r-rl)/rl        6429    5.0151227593021562D-02
>> 
>> Willie
>> 
>> 
>> 
>> 
>>> On Jun 15, 2015, at 7:49 AM, Roni Waldman <roni181066 at gmail.com> wrote:
>>> 
>>> Hi Pablo,
>>> 
>>> Oops... right, I put the file in the wrong place.
>>> Rerunning, and let's see what happens.
>>> 
>>> I also ran the rotating binary example.
>>> Ended up with donor mass 2.45 Msun and accretor reaching Roche-lobe.
>>> Does that sound right?
>>> 
>>> Thanks,
>>> Roni
>>> 
>>> ‫בתאריך יום ב׳, 15 ביוני 2015 ב-14:11 מאת ‪Pablo Marchant‬‏ <‪pamarca at gmail.com‬‏>:‬
>>> Roni, together with the inlists there is a run_star_extras.f provided (and in the inlists it is said that you need it). Are you using it? From you output it seems not. This implements the switch from a Kudritski wind to a Hamann one as described in Wellstein & Langer 1999.
>>> 
>>> I guess changing this will make it consistent with that of the paper. Your model is doing some crazy things though.
>>> 
>>> On Mon, Jun 15, 2015 at 9:06 AM, Roni Waldman <roni181066 at gmail.com> wrote:
>>> Hi Pablo,
>>> Too big for one email, will send several.
>>> Thanks,
>>> Roni
>>> 
>>> ‫בתאריך יום א׳, 14 ביוני 2015 ב-22:55 מאת ‪Pablo Marchant‬‏ <‪pamarca at gmail.com‬‏>:‬
>>> Roni, can you send me your history.data for each star and binary_history.data to compare?
>>> 
>>> On Sun, Jun 14, 2015 at 8:35 PM, Roni Waldman <roni181066 at gmail.com> wrote:
>>> Hi,
>>> 
>>> I downloaded the inlists of the MESA3 paper, and ran the example 6_massive_binary, without changing anything in the inlists or the code, version 7624.
>>> The run ends with the donor having 2Msun, the accretor filling its Roche-lobe, and the time steps being super small (1e-12).
>>> This not what seems to happen in the paper.
>>> Am I doing anything wrong?
>>> 
>>> Roni
>>> 
>>> 
>>> ------------------------------------------------------------------------------
>>> 
>>> _______________________________________________
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net
>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>> 
>>> 
>>> 
>>> 
>>> -- 
>>> Pablo Marchant Campos
>>> M.Sc on Astrophysics, Universidad Católica de Chile
>>> PhD student, Argelander-Institut für Astronomie
>>> 
>>> 
>>> 
>>> -- 
>>> Pablo Marchant Campos
>>> M.Sc on Astrophysics, Universidad Católica de Chile
>>> PhD student, Argelander-Institut für Astronomie
>>> ------------------------------------------------------------------------------
>>> _______________________________________________
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net
>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>> 
>> 
>> 
>> 
>> -- 
>> Pablo Marchant Campos
>> M.Sc on Astrophysics, Universidad Católica de Chile
>> PhD student, Argelander-Institut für Astronomie
> 
> ------------------------------------------------------------------------------
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
> 
> 
> 
> -- 
> Pablo Marchant Campos
> M.Sc on Astrophysics, Universidad Católica de Chile
> PhD student, Argelander-Institut für Astronomie
> 
> 
> 
> -- 
> Pablo Marchant Campos
> M.Sc on Astrophysics, Universidad Católica de Chile
> PhD student, Argelander-Institut für Astronomie
> ------------------------------------------------------------------------------
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150616/0fa3f2dd/attachment.html>