Cosmological models with a dynamical dark energy field typically lead to a modified propagation of gravitational waves via an effectively time-varying gravitational coupling $G(t)$. The local variation of this coupling between the time of emission and detection can be probed with standard sirens. Here we discuss the role that Lunar Laser Ranging (LLR) and binary pulsar constraints play in the prospects of constraining $G(t)$ with standard sirens. In particular, we argue that LLR constrains the matter-matter gravitational coupling $G_N(t)$, whereas binary pulsars and standard sirens constrain the quadratic kinetic gravity self-interaction $G_{gw}(t)$. Generically, these two couplings could be different in alternative cosmological models, in which case LLR constraints are irrelevant for standard sirens. We use the Hulse-Taylor pulsar data and show that observations are highly insensitive to time variations of $G_{gw}(t)$ yet highly sensitive to $G_N(t)$. We thus conclude that future gravitational waves data will become the best probe to test $G_{gw}(t)$, and will hence provide novel constraints on dynamical dark energy models.

We consider cosmological models with a dynamical dark energy field, and study the presence of three types of commonly found instabilities, namely ghost (when fields have negative kinetic energy), gradient (negative momentum squared) and tachyon (negative mass squared). In particular, we study the linear scalar perturbations of theories with two interacting scalar fields as a proxy for a dark energy and matter fields, and explicitly show how canonical transformations relate these three types of instabilities with each other. We generically show that low-energy ghosts are equivalent to tachyonic instabilities, and that high-energy ghosts are equivalent to gradient instabilities. Via examples we make evident the fact that whenever one of these fields exhibits an instability then the entire physical system becomes unstable, with an unbounded Hamiltonian. Finally, we discuss the role of interactions between the two fields, and show that whereas most of the time interactions will not determine whether an instability is present or not, they may affect the timescale of the instability. We also find exceptional cases in which the two fields are ghosts and hence the physical system is seemingly unstable, but the presence of interactions actually lead to stable solutions. These results are very important for assessing the viability of dark energy models that may exhibit ghost, gradient or tachyonic modes.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction.

After an optical peak, a classical or recurrent nova settles into a brief (days to years) period of quasi-stable thermonuclear burning in a compact configuration nearly at the white dwarf (WD) radius. During this time, the underlying WD becomes visible as a strong emitter of supersoft X-rays. Observations during this phase have revealed oscillations in the X-ray emission with periods on the order of tens of seconds. A proposed explanation for the source of these oscillations are internal gravity waves excited by nuclear reactions at the base of the hydrogen-burning layer. In this work, we present the first models exhibiting unstable surface $g$-modes with periods similar to oscillation periods found in galactic novae. However, when comparing mode periods of our models to the observed oscillations of several novae, we find that the modes which are excited have periods shorter than that observed.

We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models from the pre main-sequence to core O-depletion at solar and subsolar metallicities for a total of 2000 Monte Carlo stellar models. For each stellar model, we independently and simultaneously sample 665 thermonuclear reaction rates and use them in a $\texttt{MESA}$ in situ reaction network that follows 127 isotopes from $^{1}$H to $^{64}$Zn. With this framework we survey the core mass, burning lifetime, composition, and structural properties at five different evolutionary epochs. At each epoch we measure the probability distribution function of the variations of each property and calculate Spearman Rank-Order Correlation coefficients for each sampled reaction rate to identify which reaction rate has the largest impact on the variations on each property. We find that uncertainties in $^{14}$N$(p,\gamma)^{15}$O, triple-$\alpha$, $^{12}$C$(\alpha,\gamma)^{16}$O, $^{12}$C($^{12}$C,$p$)$^{23}$Na, $^{12}$C($^{16}$O,$p$)$^{27}$Al, $^{16}$O($^{16}$O,$n$)$^{31}$S, $^{16}$O($^{16}$O,$p$)$^{31}$P, and $^{16}$O($^{16}$O,$\alpha$)$^{28}$Si reaction rates dominate the variations of the properties surveyed. We find that variations induced by uncertainties in nuclear reaction rates grow with each passing phase of evolution, and at core H-, He-depletion are of comparable magnitude to the variations induced by choices of mass resolution and network resolution. However, at core C-, Ne-, and O-depletion, the reaction rate uncertainties can dominate the variation causing uncertainty in various properties of the stellar model in the evolution towards iron core-collapse.

We present a measurement of the volumetric rate of superluminous supernovae (SLSNe) at z~1, measured using archival data from the first four years of the Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS). We develop a method for the photometric classification of SLSNe to construct our sample. Our sample includes two previously spectroscopically-identified objects, and a further new candidate selected using our classification technique. We use the point-source recovery efficiencies from Perrett et.al. (2010) and a Monte Carlo approach to calculate the rate based on our SLSN sample. We find that the three identified SLSNe from SNLS give a rate of 91 (+76/-36) SNe/Yr/Gpc^3 at a volume-weighted redshift of z=1.13. This is equivalent to 2.2 (+1.8/-0.9) x10^-4 of the volumetric core collapse supernova rate at the same redshift. When combined with other rate measurements from the literature, we show that the rate of SLSNe increases with redshift in a manner consistent with that of the cosmic star formation history. We also estimate the rate of ultra-long gamma ray bursts (ULGRBs) based on the events discovered by the Swift satellite, and show that it is comparable to the rate of SLSNe, providing further evidence of a possible connection between these two classes of events. We also examine the host galaxies of the SLSNe discovered in SNLS, and find them to be consistent with the stellar-mass distribution of other published samples of SLSNe.

We present observations of four rapidly rising (t_{rise}~10d) transients with peak luminosities between those of supernovae (SNe) and superluminous SNe (M_{peak}~-20) - one discovered and followed by the Palomar Transient Factory (PTF) and three by the Supernova Legacy Survey (SNLS). The light curves resemble those of SN 2011kl, recently shown to be associated with an ultra-long-duration gamma ray burst (GRB), though no GRB was seen to accompany our SNe. The rapid rise to a luminous peak places these events in a unique part of SN phase space, challenging standard SN emission mechanisms. Spectra of the PTF event formally classify it as a Type II SN due to broad Halpha emission, but an unusual absorption feature, which can be interpreted as either high velocity Halpha (though deeper than in previously known cases) or Si II (as seen in Type Ia SNe), is also observed. We find that existing models of white dwarf detonations, CSM interaction, shock breakout in a wind (or steeper CSM) and magnetar spindown can not readily explain the observations. We consider the possibility that a "Type 1.5 SN" scenario could be the origin of our events. More detailed models for these kinds of transients and more constraining observations of future such events should help better determine their nature.

Novae undergo a supersoft X-ray phase of varying duration after the optical outburst. Such transient post-nova supersoft X-ray sources (SSSs) are the majority of the observed SSSs in M31. In this paper, we use the post-nova evolutionary models of Wolf et al. to compute the expected population of post-nova SSSs in M31. We predict that depending on the assumptions about the WD mass distribution in novae, at any instant there are about 250-600 post-nova SSSs in M31 with (unabsorbed) 0.2-1.0 keV luminosity L_x>10^36 erg/s. Their combined unabsorbed luminosity is of the order of ~10^39 erg/s. Their luminosity distribution shows significant steepening around log(L_x)~37.7-38 and becomes zero at L_x~2x10^38 erg/s, the maximum L_x achieved in the post-nova evolutionary tracks. Their effective temperature distribution has a roughly power law shape with differential slope of ~4-6 up to the maximum temperature of T_eff~1.5x10^6 K. We compare our predictions with the results of the XMM-Newton monitoring of the central field of M31 between 2006 and 2009. The predicted number of post-nova SSSs exceed the observed number by a factor of ~2-5, depending on the assumed WD mass distribution in novae. This is good agreement, considering the number and magnitude of uncertainties involved in calculations of the post-nova evolutionary models and their X-ray output. Furthermore, only a moderate circumstellar absorption, with hydrogen column density of the order of ~10^21 cm^-2, will remove the discrepancy.

The expanding ejecta from a classical nova remains hot enough ($\sim10^{4}\, {\rm K}$) to be detected in thermal radio emission for up to years after the cessation of mass loss triggered by a thermonuclear instability on the underlying white dwarf (WD). Nebular spectroscopy of nova remnants confirms the hot temperatures observed in radio observations. During this same period, the unstable thermonuclear burning transitions to a prolonged period of stable burning of the remnant hydrogen-rich envelope, causing the WD to become, temporarily, a super-soft X-ray source. We show that photoionization heating of the expanding ejecta by the hot WD maintains the observed nearly constant temperature of $(1-4)\times10^4\mathrm{~K}$ for up to a year before an eventual decline in temperature due to either the cessation of the supersoft phase or the onset of a predominantly adiabatic expansion. We simulate the expanding ejecta using a one-zone model as well as the Cloudy spectral synthesis code, both incorporating the time-dependent WD effective temperatures for a range of masses from $0.60\ M_{\odot}$ to $1.10\ M_{\odot}$. We show that the duration of the nearly isothermal phase depends most strongly on the velocity and mass of the ejecta and that the ejecta temperature depends on the WD's effective temperature, and hence its mass.

We report the optical identification of the companion to the {\it Fermi} black widow millisecond pulsar PSR J1544+4937. We find a highly variable source on Keck LRIS images at the nominal pulsar position, with 2 magnitude variations over orbital period in the B, g, R, and I bands. The nearly achromatic light curves are difficult to explain with a simply irradiated hemisphere model, and suggest that the optical emission is dominated by a nearly isothermal hot patch on the surface of the companion facing the pulsar. We roughly constrain the distance to PSR J1544+4937 to be between 2 and 5 kpc. A more reliable distance measurement is needed in order to constrain the composition of the companion.

The iPTF detection of the most recent outburst of the recurrent nova system RX J0045.4+4154 in the Andromeda Galaxy has enabled the unprecedented study of a massive ($M>1.3\ M_\odot$) accreting white dwarf (WD). We detected this nova as part of the near daily iPTF monitoring of M31 to a depth of $R\approx 21$\,mag and triggered optical photometry, spectroscopy and soft X-ray monitoring of the outburst. Peaking at an absolute magnitude of $M_R=-6.6$ mag, and with a decay time of 1 mag per day, it is a faint and very fast nova. It shows optical emission lines of He/N and expansion velocities of 1900 to 2600 km s$^{-1}$ 1--4 days after the optical peak. The {\it Swift} monitoring of the X-ray evolution revealed a supersoft source (SSS) with $kT_{\rm eff}\approx 90-110\ {\rm eV}$ that appeared within 5 days after the optical peak, and lasted only 12 days. Most remarkably, this is not the first event from this system, rather it is a recurrent nova with a time between outbursts of approximately 1 year, the shortest known. Recurrent X-ray emission from this binary was detected by ROSAT in 1992 and 1993, and the source was well characterized as a $M>1.3\ M_\odot$ WD SSS. Based on the observed recurrence time between different outbursts, the duration and effective temperature of the SS phase, MESA models of accreting WDs allow us to constrain the accretion rate to $\dot M>1.7\times10^{-7}\ {M_{\odot}\ {\rm yr}}^{-1}$ and WD mass $>1.30\ M_{\odot}$. If the WD keeps $30\%$ of the accreted material, it will take less than a Myr to reach core densities high enough for carbon ignition (if made of C/O) or electron capture (if made of O/Ne) to end the binary evolution.

We examine the properties of white dwarfs (WDs) accreting hydrogen-rich matter in and near the stable burning regime of accretion rates as modeled by time-dependent calculations done with Modules for Experiments in Stellar Astrophysics (MESA). We report the stability boundary for WDs of masses between 0.51 solar masses and 1.34 solar masses as found via time-dependent calculations. We also examine recurrent novae that are accreting at rates close to, but below, the stable burning limit and report their recurrence times and ignition masses. Our dense grid in accretion rates finds the expected minimum possible recurrence times as a function of the WD mass. This enables inferences to be made about the minimum WD mass possible to reach a specific recurrence time. We compare our computational models of post-outburst novae to the stably burning WDs and explicitly calculate the duration and effective temperature (Teff) of the post-novae WD in the supersoft phase. We agree with the measured turnoff time - Teff relation in M31 by Henze and collaborators, infer WD masses in the 1.0-1.3 solar masses range, and predict ejection masses consistent with those observed. We close by commenting on the importance of the hot helium layer generated by stable or unstable hydrogen burning for the short- and long-term evolution of accreting white dwarfs.