Global non-hydrostatic atmospheric models are becoming increasingly important for studying the climates of planets and exoplanets. However, such models suffer from computational difficulties due to the large aspect ratio between the horizontal and vertical directions. To overcome this problem, we developed a global model using a vertically-implicit-correction (VIC) scheme in which the integration time step is no longer limited by the propagation of acoustic waves in the vertical. We proved that our model, based on the $\rm Athena^{++}$ framework and its extension for planetary atmospheres - SNAP (Simulating Non-hydrostatic Atmosphere on Planets), rigorously conserves mass and energy in finite volume simulations. We found that traditional numerical stabilizers such as hyper-viscosity and divergence damping are not needed when using the VIC scheme, which greatly simplifies the numerical implementation and improves stability. We present simulation results ranging from 1D linear waves to 3D global circulations with and without the VIC scheme. These tests demonstrate that our formulation correctly tracks local turbulent motions, produces Kelvin-Helmholtz instability, and generates a super-rotating jet on hot Jupiters. Employing this VIC scheme improves the computational efficiency of global simulations by more than two orders of magnitude compared to an explicit model and facilitates the capability of simulating a wide range of planetary atmospheres both regionally and globally.

Binary evolution is indispensable in stellar evolution to understand the formation and evolution of most peculiar and energetic objects, such as binary compact objects, Type Ia supernovae, X-ray binaries, cataclysmic variables, blue stragglers, hot subdwarfs, and central binaries in planetary nebulae. Mass transfer in binary stars can change the evolutionary path and fate of the corresponding objects relative to what is expected from single stellar evolution. What is the critical mass ratio at which unstable mass transfer occurs is an unsolved fundamental problem in binary evolution. To resolve this issue, we construct the thermal equilibrium mass loss model and derive critical mass ratios for both thermal timescale mass transfer and unstable mass transfer, the latter of which occurs when the outer Lagrangian point, L2, is overfilled. Using several 3.2 Msun stellar models as examples, we study the stellar response to thermal equilibrium mass loss and present the thresholds for thermal timescale mass transfer. We study the possible mass transfer channels of binary systems containing a 3.2 Msun donor star, taking into account thermal timescale mass transfer, unstable mass transfer through L2, and dynamical timescale mass transfer. We repeat this simulation for a grid of donor stars with different masses (from 0.1 to 100 Msun with Z = 0.02) and at different evolutionary stages, and present our results. The results show that unstable mass transfer due to the overfilling of the outer Lagrangian point may also play an essential role in the formation of common envelopes for late red giant branch and asymptotic giant branch donors.

Rotational modulations are observed on brown dwarfs and directly imaged exoplanets, but the underlying mechanism is not well understood. Here, we analyze Jupiter's rotational light curves at 12 wavelengths from the ultraviolet (UV) to the mid-infrared (mid-IR). Peak-to-peak amplitudes of Jupiter's light curves range from sub percent to 4% at most wavelengths, but the amplitude exceeds 20% at 5 $\rm \mu m$, a wavelength sensing Jupiter's deep troposphere. Jupiter's rotational modulations are primarily caused by discrete patterns in the cloudless belts instead of the cloudy zones. The light-curve amplitude is dominated by the sizes and brightness contrasts of the Great Red Spot (GRS), expansions of the North Equatorial Belt (NEB), patchy clouds in the North Temperate Belt (NTB) and a train of hot spots in the NEB. In reflection, the contrast is controlled by upper tropospheric and stratospheric hazes, clouds, and chromophores in the clouds. In thermal emission, the small rotational variability is caused by the spatial distribution of temperature and opacities of gas and aerosols; the large variation is caused by the $\rm NH_{3}$ cloud holes and thin-thick clouds. The methane-band light curves exhibit opposite-shape behavior compared with the UV and visible wavelengths, caused by wavelength-dependent brightness change of the GRS. Light-curve evolution is induced by periodic events in the belts and longitudinal drifting of the GRS and patchy clouds in the NTB. This study suggests several interesting mechanisms related to distributions of temperature, gas, hazes, and clouds for understanding the observed rotational modulations on brown dwarfs and exoplanets.

Context. The companions of the exploding carbon-oxygen white dwarfs (CO WDs) for producing type Ia supernovae (SNe Ia) are still not conclusively confirmed. A red-giant (RG) star has been suggested to be the mass donor of the exploding WD, named as the symbiotic channel. However, previous studies on the this channel gave a relatively low rate of SNe Ia. Aims. We aim to systematically investigate the parameter space, Galactic rates and delay time distributions of SNe Ia from the symbiotic channel by employing a revised mass-transfer prescription. Methods. We adopted an integrated mass-transfer prescription to calculate the mass-transfer process from a RG star onto the WD. In this prescription, the mass-transfer rate varies with the local material states. Results. We evolved a large number of WD+RG systems, and found that the parameter space of WD+RG systems for producing SNe Ia is significantly enlarged. This channel could produce SNe Ia with intermediate and old ages, contributing to at most 5% of all SNe Ia in the Galaxy. Our model increases the SN Ia rate from this channel by a factor of 5. We suggest that the symbiotic systems RS Oph and T CrB are strong candidates for the progenitors of SNe Ia.

In the limit of extremely rapid mass transfer, the response of a donor star in an interacting binary becomes asymptotically one of adiabatic expansion. We survey here adiabatic mass loss from Population I stars of mass 0.10 Msun to 100 Msun from the zero age main sequence to the base of the giant branch, or to central hydrogen exhaustion for lower main sequence stars. For intermediate- and high-mass stars, dynamical mass transfer is preceded by an extended phase of thermal time scale mass transfer as the star is stripped of most of its envelope mass. The critical mass ratio qad above which this delayed dynamical instability occurs increases with advancing evolutionary age of the donor star, by ever-increasing factors for more massive donors. Most intermediate- or high-mass binaries with nondegenerate accretors probably evolve into contact before manifesting this instability. As they approach the base of the giant branch, however, and begin developing a convective envelope, qad plummets dramatically among intermediate-mass stars, to values of order unity, and a prompt dynamical instability occurs. Among low-mass stars, the prompt instability prevails throughout main sequence evolution, with q_ad declining with decreasing mass. Our calculated qad agree well with the behavior of timedependent models by Chen & Han (2003) of intermediate-mass stars initiating mass transfer in the Hertzsprung gap. Application of our results to cataclysmic variables, as systems which must be stable against rapid mass transfer, nicely circumscribes the range in qad as a function of orbital period in which they are found. These results are intended to advance the verisimilitude of population synthesis models of close binary evolution.

It is believed that, except for metallicity, some other parameters are needed to explain the horizontal branch (HB) morphology of globular clusters (GCs). Furthermore, these parameters are considered to be correlated with the mass loss of the red giant branch (RGB) stars. In our previous work, we proposed that tidally enhanced stellar wind during binary evolution may affect the HB morphology by enhancing the mass loss of the red giant primary. As a further study, we now investigate the effects of metallicity and age on HB morphology by considering tidally enhanced stellar winds during binary evolution. We incorporated the tidally enhanced-stellar-wind model into Eggleton's stellar evolution code to study the binary evolution. To study the effects of metallicity and age on our final results, we conducted two sets of model calculations: (i) for a fixed age, we used three metallicities, namely Z=0.0001, 0.001, and 0.02. (ii) For a fixed metallicity, Z=0.001, we used five ages in our model calculations: 14, 13, 12, 10, and 7 Gyr. We found that HB morphology of GCs becomes bluer with decreasing metallicity, and old GCs present bluer HB morphology than young ones. These results are consistent with previous work. Although the envelope-mass distributions of zero-age HB stars produced by tidally enhanced stellar wind are similar for different metallicities, the synthetic HB under tidally enhanced stellar wind for Z=0.02 presented a distinct gap between red and blue HB. However, this feature was not seen clearly in the synthetic HB for Z=0.001 and 0.0001. We also found that higher binary fractions may make HB morphology become bluer, and we discussed the results with recent observations.

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

The stability of mass transfer is important in the formation of contact binaries from detached binaries when the primaries of the initially detached binaries fill their Roche lobes. Using Eggleton's stellar evolution code, we investigate the formation and the short-period limit of contact binaries by considering the effect of the instability of mass transfer. It is found that with decreasing initial primary mass from 0.89M$_{\rm \odot}$ to 0.63M$_{\rm \odot}$, the range of the initial mass ratio decreases for detached binaries that experience stable mass transfer and evolve into contact. If the initial primary mass is less than 0.63M$_{\rm \odot}$, detached binaries would experience dynamically unstable mass transfer when the primaries of detached binaries fill their Roche lobes. These systems would evolve into a common envelope situation and probably then to a complete merger of two components on a quite short timescale. This results in a low mass limit at about 0.63M$_{\rm \odot}$ for the primary mass of contact binaries, which might be a main reason why the period distribution of contact binaries has a short limit of about 0.22 days. By comparing the theoretical period distribution of contact binaries with the observational data, it is found that the observed contact binaries are above the low mass limit for the primary mass of contact binaries and no observed contact binaries are below this limit. This suggests that the short-period limit of contact binaries can be explained by the instability of the mass transfer that occurs when the primaries of the initially detached binaries fill their Roche lobes.

The asymptotic response of donor stars in interacting binary systems to very rapid mass loss is characterized by adiabatic expansion throughout their interiors. In this limit, energy generation and heat flow through the stellar interior can be neglected. We model this response by constructing model sequences, beginning with a donor star filling its Roche lobe at an arbitrary point in its evolution, holding its specific entropy and composition profiles fixed as mass is removed from the surface. The stellar interior remains in hydrostatic equilibrium. Luminosity profiles in these adiabatic models of mass-losing stars can be reconstructed from the specific entropy profiles and their gradients. These approximations are validated by comparison with time-dependent binary mass transfer calculations. We describe how adiabatic mass loss sequences can be used to quantify threshold conditions for dynamical time scale mass transfer, and to establish the range of post-common envelope binaries that are allowed energetically.

Roche-lobe overflow and common envelope evolution are very important in binary evolution, which is believed to be the main evolutionary channel to hot subdwarf stars. The details of these processes are difficult to model, but adiabatic expansion provides an excellent approximation to the structure of a donor star undergoing dynamical time scale mass transfer. We can use this model to study the responses of stars of various masses and evolutionary stages as potential donor stars, with the urgent goal of obtaining more accurate stability criteria for dynamical mass transfer in binary population synthesis studies. As examples, we describe here several models with the initial masses equal to 1 Msun and 10 Msun, and identify potential limitations to the use of our results for giant-branch stars.