The R Coronae Borealis (RCB) stars are extremely hydrogen-deficient carbon stars which produce large amounts of dust, causing sudden deep declines in brightness. They are believed to be formed primarily through white dwarf mergers. In this paper, we use MESA to investigate how post-merger objects with a range of He-burning shell temperatures from 2.1 - 5.4 x 10^8 K with solar and subsolar metallicities evolve into RCB stars. The most successful model of these has subsolar metallicity and a temperature near 3 x 10^8 K. We find a strong dependence on He-burning shell temperature for surface abundances of elements involved in the CNO cycle, as well as differences in effective temperature and radius of RCBs. We find around 1 dex diminished surface abundances for the models with 10% subsolar metallicity, which is expected. Models with subsolar metallicities also exhibit longer lifetimes than their solar counterparts. Additionally, we find that convective mixing of the burned material occurs only in the first few years of post-merger evolution, after which the surface abundances are constant during and after the RCB phase, providing evidence for why these stars show a strong enhancement of partial He-burning products.

We explore the possibility that the star alpha Orionis (Betelgeuse) is the outcome of a merger that occurred in a low mass ratio (q = M2/M1 = 0.07 - 0.25) binary system some time in the past hundreds of thousands of years. To that goal, we present a simple analytical model to approximate the perturbed internal structure of a post-merger object following the coalescence of a secondary in the mass range 1-4 Msun into the envelope of a 15-17 Msun primary. We then compute the long-term evolution of post-merger objects for a grid of initial conditions and make predictions about their surface properties for evolutionary stages that are consistent with the observed location of Betelgeuse in the Hertzsprung-Russell diagram. We find that if a merger occurred after the end of the primary's main-sequence phase, while it was expanding toward becoming a red supergiant star and typically with radius ~200 - 300 Rsun, then it's envelope is spun-up to values which remain in a range consistent with the Betelgeuse observations for thousands of years of evolution. We argue that the best scenario that can explain both the fast rotation of Betelgeuse and its observed large space velocity is one where a binary was dynamically ejected by its parent cluster a few million years ago and then subsequently merged. An alternative scenario in which the progenitor of Betelgeuse was spun up by accretion in a binary and released by the supernova explosion of the companion requires a finely tuned set of conditions but cannot be ruled out.

The R Coronae Borealis (RCB) stars are rare hydrogen--deficient, carbon--rich supergiants. They undergo extreme, irregular declines in brightness of many magnitudes due to the formation of thick clouds of carbon dust. It is thought that RCB stars result from the mergers of CO/He white dwarf (WD) binaries. We constructed post--merger spherically asymmetric models computed with the MESA code, and then followed the evolution into the region of the HR diagram where the RCB stars are located. We also investigated nucleosynthesis in the dynamically accreting material of CO/He WD mergers which may provide a suitable environment for significant production of 18O and the very low 16O/18O values observed. We have also discovered that the N abundance depends sensitively on the peak temperature in the He--burning shell. Our MESA modeling consists of engineering the star by adding He--WD material to an initial CO--WD model, and then following the post--merger evolution using a nuclear--reaction network to match the observed RCB abundances as it expands and cools to become an RCB star. These new models are more physical because they include rotation, mixing, mass-loss, and nucleosynthesis within MESA. We follow the later evolution beyond the RCB phase to determine the stars' likely lifetimes. The relative numbers of known RCB and Extreme Helium (EHe) stars correspond well to the lifetimes predicted from the MESA models. In addition, most of computed abundances agree very well with the observed range of abundances for the RCB class.

For the ordinary classical nova QZ Aurigae (which erupted in 1964), we report 1317 magnitudes from 1912--2016, including four eclipses detected on archival photographic plates from long before the eruption. We have accurate and robust measures of the orbital period both pre-eruption and post-eruption, and we find that the orbital period decreased, with a fractional change of -290.71+-0.28 parts-per-million across the eruption, with the orbit necessarily getting smaller. Further, we find that the light curve outside of eclipses and eruption is flat at near B=17.14 from 1912--1981, whereupon the average light curve starts fading down to B=17.49 with large variability. QZ Aur is a robust counter-example against the Hibernation model for the evolution of cataclysmic variables, where the model requires that all novae have their period increase across eruptions. Large period decreases across eruptions can easily arise from mass imbalances in the ejecta, as are commonly seen in asymmetric nova shells.

We assume that dark matter is comprised of axion-like particles (ALPs) generated by the realignment mechanism in the post-inflationary scenario. This leads to isocurvature fluctuations with an amplitude of order one for scales comparable to the horizon at the time when the ALP field starts oscillating. The power spectrum of these fluctuations is flat for small wave numbers, extending to scales relevant for cosmological observables. Denoting the relative isocurvature amplitude at $k_*$ = 0.05 Mpc$^{-1}$ by $f_{\rm iso}$, Planck observations of the cosmic microwave background (CMB) yield $f_{\rm iso}$ < 0.31 at the 2$\sigma$-level. This excludes the hypothesis of post-inflationary ALP dark matter with masses $m_{a}$ < 10$^{-20}-$10$^{-16}$ eV, where the range is due to details of the ALP mass-temperature dependence. Future CMB stage IV and 21-cm intensity mapping experiments may improve these limits by 1$-$2 orders of magnitude in $m_{a}$.

We present the first self-consistent, three dimensional study of hydrodynamic simulations of mass transfer in binary systems with bipolytropic (composite polytropic) components. In certain systems, such as contact binaries or during the common envelope phase, the core-envelope structure of the stars plays an important role in binary interactions. In this paper, we compare mass transfer simulations of bipolytropic binary systems in order to test the suitability of our numerical tools for investigating the dynamical behaviour of such systems. The initial, equilibrium binary models possess a core-envelope structure and are obtained using the bipolytropic self-consistent field technique. We conduct mass transfer simulations using two independent, fully three-dimensional, Eulerian codes - Flow-ER and Octo-tiger. These hydrodynamic codes are compared across binary systems undergoing unstable as well as stable mass transfer, and the former at two resolutions. The initial conditions for each simulation and for each code are chosen to match closely so that the simulations can be used as benchmarks. Although there are some key differences, the detailed comparison of the simulations suggests that there is remarkable agreement between the results obtained using the two codes. This study puts our numerical tools on a secure footing, and enables us to reliably simulate specific mass transfer scenarios of binary systems involving components with a core-envelope structure.

We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two white dwarfs is $\lesssim1.0~{\rm M_\odot}$. Mergers involving a CO and a He white dwarf have been suggested as a possible formation channel for R Coronae Borealis type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of $^{18}{\mathrm O}$, the merger must avoid dredging up large amounts of $^{16}{\mathrm O}$, or else it will be difficult to produce sufficient $^{18}{\mathrm O}$ to explain the oxygen ratio observed to be of order unity. We performed a total of 9 simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes, and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, $>10^{-2}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO white dwarf with a $\sim 0.1~{\rm M_\odot}$ layer of helium on top, we find that no $^{16}{\mathrm O}$ is being dredged up, while in the $q=0.8$ simulation $<10^{-4}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We demonstrate that rapidly rotating bipolytropic (composite polytropic) stars and toroidal disks can be obtained using Hachisu's self consistent field technique. The core and the envelope in such a structure can have different polytropic indices and also different average molecular weights. The models converge for high $T/|W|$ cases, where T is the kinetic energy and W is the gravitational energy of the system. The agreement between our numerical solutions with known analytical as well as previously calculated numerical results is excellent. We show that the uniform rotation lowers the maximum core mass fraction or the Sch$\ddot{\rm{o}}$nberg-Chandrasekhar limit for a bipolytropic sequence. We also discuss the applications of this method to magnetic braking in low mass stars with convective envelopes.

Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms --- a finite-volume "grid" code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations an atmosphere forms around the accretor which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.

We show that in realistic cases of accretion in active galactic nuclei or stellar-mass X-ray binaries, the Lense-Thirring effect breaks the central regions of tilted accretion discs around spinning black holes into a set of distinct planes with only tenuous flows connecting them. If the original misalignment of the outer disc to the spin axis of the hole is $45^{\circ} \lesssim \theta \lesssim 135^{\circ}$, as in $\sim 70$% of randomly oriented accretion events, the continued precession of these discs sets up partially counter-rotating gas flows. This drives rapid infall as angular momentum is cancelled and gas attempts to circularize at smaller radii. Disc breaking close to the black hole leads to direct dynamical accretion, while breaking further out can drive gas down to scales where it can accrete rapidly. For smaller tilt angles breaking can still occur, and may lead to other observable phenomena such as QPOs. For such effects not to appear, the black hole spin must in practice be negligibly small, or be almost precisely aligned with the disc. Qualitatively similar results hold for any accretion disc subject to a forced differential precession, such as an external disc around a misaligned black hole binary.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.

We present numerical simulations of dynamically unstable mass transfer in a double white dwarf binary with initial mass ratio, q = 0.4. The binary components are approximated as polytropes of index n = 3/2 and the initially synchronously rotating, semi-detached equilibrium binary is evolved hydrodynamically with the gravitational potential being computed through the solution of Poisson's equation. Upon initiating deep contact in our baseline simulation, the mass transfer rate grows by more than an order of magnitude over approximately ten orbits, as would be expected for dynamically unstable mass transfer. However, the mass transfer rate then reaches a peak value, the binary expands and the mass transfer event subsides. The binary must therefore have crossed the critical mass ratio for stability against dynamical mass transfer. Despite the initial loss of orbital angular momentum into the spin of the accreting star, we find that the accretor's spin saturates and angular momentum is returned to the orbit more efficiently than has been previously suspected for binaries in the direct impact accretion mode. To explore this surprising result, we directly measure the critical mass ratio for stability by imposing artificial angular momentum loss at various rates to drive the binary to an equilibrium mass transfer rate. For one of these driven evolutions, we attain equilibrium mass transfer and deduce that effectively q_crit has evolved to approximately 2/3. Despite the absence of a fully developed disk, tidal interactions appear effective in returning excess spin angular momentum to the orbit.

We describe the evolution of double degenerate binary systems, consisting of components obeying the zero temperature mass radius relationship for white dwarf stars, from the onset of mass transfer to one of several possible outcomes including merger, tidal disruption of the donor, or survival as a semi-detached AM CVn system. We use a combination of analytic solutions and numerical integrations of the standard orbit-averaged first-order evolution equations, including direct impact accretion and the evolution of the components due to mass exchange. We include also the effects of mass-loss during super-critical (super-Eddington) mass transfer and the tidal and advective exchanges of angular momentum between the binary components. We cover much the same ground as Marsh et al.(2004) with the additional effects of the advective or consequential angular momentum from the donor and its tidal coupling to the orbit which is expected to be stronger than that of the accretor. With the caveat that our formalism does not include an explicit treatment of common envelope phases, our results suggest that a larger fraction of detached double white dwarfs than what has been hitherto assumed, survive the onset of mass transfer, even if this mass transfer is initially unstable and rises to super-Eddington levels. In addition, as a consequence of the tidal coupling, systems that come into contact near the mass transfer instability boundary undergo a phase of oscillation cycles in their orbital period (and other system parameters). Much of the formalism developed here would also apply to other mass-transferring binaries, and in particular to cataclysmic variables and Algol systems.

Hydrodynamical simulations of semi-detached, polytropic binary stars are presented in an effort to study the onset and stability of dynamical mass transfer events. Initial, synchronously rotating equilibrium models are constructed using a self-consistent-field technique and then evolved with an Eulerian hydrodynamics code in a fully self-consistent manner. We describe code improvements introduced over the past few years that permit us to follow dynamical mass-transfer events through more than 30 orbits. Mass-transfer evolutions are presented for two different initial configurations: A dynamically unstable binary with initial mass ratio (donor/accretor) $q_0 = 1.3$ that leads to a complete merger in $\sim 10$ orbits; and a double-degenerate binary with initial mass ratio $q_0 = 0.5$ that, after some initial unstable growth of mass transfer, tends to separate as the mass-transfer rate levels off.

The masses of supermassive black holes (SBHs) show correlations with bulge properties in disk and elliptical galaxies. We study the formation of galactic structure within flat-core triaxial haloes and show that these correlations can be understood within the framework of a baryonic component modifying the orbital structure in the underlying potential. In particular, we find that terminal properties of bulges and their central SBHs are constrained by the destruction of box orbits in the harmonic cores of dark haloes and the emergence of progressively less eccentric loop orbits there. SBH masses, M_SBH, should exhibit a tighter correlation with bulge velocity dispersions, sigma_B, than with bulge masses, M_B, in accord with observations, if there is a significant scatter in the M_H-sigma_H relation for the halo. In the context of this model the observed M_SBH-sigma_B relation implies that haloes should exhibit a Faber-Jackson type relationship between their masses and velocity dispersions. The most important prediction of our model is that halo properties determine the bulge and SBH parameters. The model also has important implications for galactic morphology and the process of disk formation.

We describe computational tools that have been developed to simulate dynamical mass transfer in semi-detached, polytropic binaries that are initially executing synchronous rotation upon circular orbits. Initial equilibrium models are generated with a self-consistent field algorithm; models are then evolved in time with a parallel, explicit, Eulerian hydrodynamics code with no assumptions made about the symmetry of the system. Poisson's equation is solved along with the equations of ideal fluid mechanics to allow us to treat the nonlinear tidal distortion of the components in a fully self-consistent manner. We present results from several standard numerical experiments that have been conducted to assess the general viability and validity of our tools, and from benchmark simulations that follow the evolution of two detached systems through five full orbits (up to approximately 90 stellar dynamical times). These benchmark runs allow us to gauge the level of quantitative accuracy with which simulations of semi-detached systems can be performed using presently available computing resources. We find that we should be able to resolve mass transfer at levels $\dot{M} / M > few x 10^-5$ per orbit through approximately 20 orbits with each orbit taking about 30 hours of computing time on parallel computing platforms.

Warped H I gas layers in the outer regions of spiral galaxies usually display a noticeably twisted structure. This structure almost certainly arises primarily as a result of differential precession in the H I disk as it settles toward a preferred orientation in an underlying dark halo potential well that is not spherically symmetric. In an attempt to better understand the structure and evolution of these twisted, warped disk structures, we have adopted the "twist-equation" formalism originally developed by Petterson (1977) to study accretion onto compact objects. Utilizing more recent treatments of this formalism, we have generalized the twist-equation to allow for the treatment of non-Keplerian disks and from it have derived a steady-state structure of twisted disks that develops from free precession in a nonspherical, logarithmic halo potential. We have used this steady-state solution to produce H I maps of five galaxies (M83, NGC 300, NGC 2841, NGC 5033, NGC5055), which match the general features of the observed maps of these galaxies quite well. In addition, the model provides an avenue through which the kinematical viscosity of the H I disk and the quadrupole distortion of the dark halo in each galaxy can be quantified. This generalized equation can also be used to examine the time-evolutionary behavior of warped galaxy disks.

The orbital evolution of cataclysmic variables in which the companion is illuminated by a fraction of the accretion luminosity consists of irradiation-driven limit cycles on thermal timescales, superimposed on a secular evolution toward shorter periods due to systemic angular momentum losses. We show that positive orbital period derivatives during bright phases are a natural consequence of the expansion of the companion during high mass transfer phases in the limit cycle. The irradiation instability may be enhanced by consequential angular momentum losses, CAML, accompanying the limit cycle. We investigate the secular evolution of cataclysmic binaries under the combined effects of irradiation and CAML and show that faster than secular transfer fluctuations that occur during these cycles can account for the observed dispersion in disk luminosities or estimated accretion rates at a given orbital period. If indeed irradiation-driven and CAML--assisted mass transfer fluctuations on timescales faster than secular occur, as discussed in this paper, then we may be able to predict the relative abundances of the different types of cataclysmic variable at a given orbital period. For example this mechanism may explain the relative paucity of dwarf novae with respect to nova-like variables between 3 and 4 hours.

We consider the evolutionary state of the black-hole X-ray source GRO J1655-40 in the context of its transient nature. Recent optical observations show that the donor in GRO J1655-40 is an intermediate-mass star (~ 2.3 solar masses) crossing the Hertzsprung gap. Usually in such systems the donor's radius expansion drives a near-Eddington or super-Eddington mass transfer rate which would sustain a persistently bright accretion disk. We show that GRO J1655-40 is close to a narrow parameter range where disk instabilities can occur. This range corresponds to a short-lived evolutionary stage where the secondary's radius expansion stalls (or reverses), with a correspondingly lower mass transfer rate. If GRO J1655-40 belongs to this class of transients the predicted accretion rates imply large populations of luminous persistent and transient sources, which are not seen in X-rays. The transient nature of the related system GRS 1915+105 may reflect spectral variations in a bolometrically persistent source rather than a genuine luminosity increase.

We show that stable disk accretion should be very rare among low-mass X-ray binaries and cataclysmic variables whose evolution is driven by the nuclear expansion of the secondary star on the first giant branch. Stable accretion is confined to neutron-star systems where the secondary is still relatively massive, and some supersoft white dwarf accretors. All other systems, including all black-hole systems, appear as soft X-ray transients or dwarf novae. All long-period neutron-star systems become transient well before most of the envelope mass is transferred, and remain transient until envelope exhaustion. This complicates attempts to compare the numbers of millisecond pulsars in the Galactic disk with their LMXB progenitors, and also means that the pulsar spin rates are fixed in systems which are transient rather than steady, contrary to common assumption. The long-period persistent sources Sco X-2, LMC X-2, Cyg X-2 and V395 Car must have minimum companion masses > 0.75 Msun if they contain neutron stars, and still larger masses if they contain black holes. The companion in the neutron-star transient GRO J1744-2844 must have a mass <0.87 Msun. The existence of any steady sources at all at long periods supports the ideas that (a) the accretion disks in many, if not all, LMXBs are strongly irradiated by the central source, and (b) mass transfer is thermally unstable in long-period supersoft X-ray sources.

We give a global analysis of mass transfer variations in low-mass X-ray binaries and cataclysmic variables whose evolution is driven by the nuclear expansion of the secondary star. We show that limit cycles caused by irradiation of the secondary by the accreting primary are possible in a large class of these binaries. In the high state the companion transfers a large fraction of its envelope mass on a thermal timescale. In most cases this implies super-Eddington transfer rates, and would thus probably lead to common-envelope evolution and the formation of an ultrashort-period binary. Observed systems with (sub)giant secondaries stabilize themselves against this possibility either by being transient, or by shielding the secondary from irradiation in some way.