1. I discovered a mistake in the AGB star model. It is relevant to the optical depth. I will fix it this week.

2. Due to the mistake, the mass loss rate was not correct. I think the mass loss rate is crucial in the formation of circumbinary disk. 10^{-7}M_{\odot}/yr is somewhat small. If there could be a 10^{-5}M_{\odot}/yr AGB wind, the formation of the circumbinary disk will be much more easier. ​Liljegren et al. (2016) did a detailed study on AGB wind recently. I will try their model.

3. I will work with Jonathan to add the self-consistent mass loss calculation in AstroBEAR.