With the advent of modern neutrino and gravitational wave detectors, the promise of multi-messenger detections of the next galactic core-collapse supernova has become very real. Such detections will give insight into the core-collapse supernova mechanism, the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detections, a thorough understanding of the landscape of possible core-collapse supernova events, and their multi-messenger signals, is needed. We present detailed predictions of neutrino and gravitational wave signals from 1D simulations of stellar core collapse, spanning the landscape of core-collapse progenitors from $9-120\,\mathrm{M}_{\odot}$. In order to achieve explosions in 1D, we use the STIR model, which includes the effects of turbulence and convection in 1D supernova simulations to mimic the 3D explosion mechanism. We study the gravitational wave emission from the 1D simulations using an astroseismology analysis of the proto-neutron star. We find that the neutrino and gravitational wave signals are strongly correlated with the structure of the progenitor star and remnant compact object. Using these correlations, future detections of the first few seconds of neutrino and gravitational wave emission from a galactic core-collapse supernova may be able to provide constraints on stellar evolution independent of pre-explosion imaging and the mass of the compact object remnant prior to fallback accretion.

Many Type II supernovae (SNe) show hot early (~30 days) emission, and a diversity in their light curves extending from the Type IIP to the Type IIL, which can be explained by interaction with dense and confined circumstellar material (CSM). We perform hydrodynamical simulations of red supergiants to model the ejection of CSM caused by wave heating during late-stage nuclear burning. Even a small amount of deposited energy ($10^{46}-10^{47}$ erg), which is roughly that expected due to waves excited by convection in the core, is sufficient to change the shapes of SN light curves and bring them into better agreement with observations. As a test case, we consider the specific example of SN 2017eaw, which shows that a nuclear burning episode is able to explain the light curve if it occurs ~150-450 days prior to core-collapse. Due to the long timescale it takes for the low energy shock to traverse the star, this would manifest as a pre-SN outburst ~50-350 days prior to the full-fledged SN. Applying work like this to other SNe will provide a direct connection between the SN and pre-SN outburst properties, which can be tested by future wide field surveys. In addition, we show that our models can qualitatively explain the short lived `flash-ionization' lines seen in the early spectra of many Type II SNe.

We study the gravitational wave signal from eight new 3D core-collapse supernova simulations. We show that the signal is dominated by $f$- and $g$-mode oscillations of the protoneutron star and its frequency evolution encodes the contraction rate of the latter, which, in turn, is known to depend on the star's mass, on the equation of state, and on transport properties in warm nuclear matter. A lower-frequency component of the signal, associated with the standing accretion shock instability, is found in only one of our models. Finally, we show that the energy radiated in gravitational waves is proportional to the amount of turbulent energy accreted by the protoneutron star.

We present extensive optical photometric and spectroscopic observations, from 4 to 482 days after explosion, of the Type II-plateau (II-P) supernova (SN) 2017eaw in NGC 6946. SN 2017eaw is a normal SN II-P intermediate in properties between, for example, SN 1999em and SN 2012aw and the more luminous SN 2004et, also in NGC 6946. We have determined that the extinction to SN 2017eaw is primarily due to the Galactic foreground and that the SN site metallicity is likely subsolar. We have also independently confirmed a tip-of-the-red-giant-branch (TRGB) distance to NGC 6946 of 7.73+/-0.78 Mpc. The distances to the SN that we have also estimated via both the standardized candle method and expanding photosphere method corroborate the TRGB distance. We confirm the SN progenitor identity in pre-explosion archival Hubble Space Telescope (HST) and Spitzer Space Telescope images, via imaging of the SN through our HST Target of Opportunity program. Detailed modeling of the progenitor's spectral energy distribution indicates that the star was a dusty, luminous red supergiant consistent with an initial mass of ~15 Msuns.

Hydrogen-rich, core-collapse supernovae are typically divided into four classes: IIP, IIL, IIn, and IIb. In general, interaction with circumstellar material is only considered for Type IIn supernovae. However, recent hydrodynamic modeling of IIP and IIL supernovae requires circumstellar material to reproduce their early light curves. In this scenario, IIL supernovae experience large amounts of mass loss before exploding. We test this hypothesis on ASASSN-15oz, a Type IIL supernova. With extensive follow-up in the X- ray, UV, optical, IR, and radio we present our search for signs of interaction, and the mass-loss history indicated by their detection. We find evidence of short-lived intense mass-loss just prior to explosion from light curve modeling, amounting in 1.5 M$_{\odot}$ of material within 1800 R$_{\odot}$ of the progenitor. We also detect the supernova in the radio, indicating mass-loss rates of $10^{-6}-10^{-7}$ M$_{\odot}$ yr$^{-1}$ prior to the extreme mass-loss period. Our failure to detect the supernova in the X-ray and the lack of narrow emission lines in the UV, optical, and NIR do not contradict this picture and place an upper limit on the mass-loss rate outside the extreme period of $<10^{-4}$ M$_{\odot}$ yr$^{-1}$. This paper highlights the importance gathering comprehensive data on more Type II supernovae to enable detailed modeling of the progenitor and supernova which can elucidate their mass-loss histories and envelope structures and thus inform stellar evolution models.

The X-ray signal from hydrogen-rich supernovae (SNe II) in the first tens to hundreds of days after the shock breakout encodes important information about the circumstellar material (CSM) surrounding their progenitors before explosion. In this study, we describe a way to generate the SN II X-ray light curves from hydrodynamical simulations performed with the code Athena++, using the X-ray package XSPEC. In addition, we employ a radiation diffusion hydrodynamic code SNEC for generating the optical light curves in different bands. In this numerical setup, we model the X-ray and optical emission from a set of progenitor models, consisting of either two (red supergiant + low density steady wind), or three (red supergiant + dense CSM + low density steady wind) components. We vary the density in the wind and the slope in the CSM to see how these parameters influence the resulting X-ray and optical light curves. Among our models, we identify one that is able to roughly reproduce both optical and X-ray data of the well observed SN 2013ej. In order to achieve this, the slope of the dense CSM in this model should be steeper than the one of a steady wind ($\rho\propto r^{-2}$), and closer to $\rho\propto r^{-5}$. On the other hand, we show that too steep and extended CSM profiles may produce excessive X-ray emission in the first few tens of days, up to a few orders of magnitude larger than observed. We conclude that ability to reproduce the observed X-ray signal from SNe~II together with their optical light curves is crucial in establishing the validity of different CSM models.

Recent modeling of hydrogen-rich Type II supernova (SN II) light curves suggests the presence of dense circumstellar material (CSM) surrounding the exploding progenitor stars. This has important implications for the activity and structure of massive stars near the end of their lives. Since previous work focused on just a few events, here we expand to a larger sample of twenty well-observed SNe II. For each event we are able to constrain the progenitor zero-age main-sequence (ZAMS) mass, explosion energy, and the mass and radial extent of dense CSM. We then study the distribution of each of these properties across the full sample of SNe. The inferred ZAMS masses are found to be largely consistent with a Salpeter distribution with minimum and maximum masses of 10.4 and 22.9 Msun, respectively. We also compare the individual ZAMS masses we measure with specific SNe II that have pre-explosion imaging to check their consistency. Our masses are generally comparable to or larger than the pre-explosion imaging masses, potentially helping ease the red supergiant problem. The explosion energies vary from (0.1-1.3)x10^51 erg, and for ~70% of the SNe we obtain CSM masses in the range between 0.18-0.83 Msun. We see a potential correlation between the CSM mass and explosion energy, which suggests that pre-explosion activity has a strong impact on the structure of the star. This may be important to take into account in future studies of the ability of the neutrino mechanism to explode stars. We also see a possible correlation between the CSM's radial extent and ZAMS mass, which could be related to the time with respect to explosion when the CSM is first generated.

We study gravitational waves (GWs) from a set of two-dimensional multi-group neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Our goal is to systematize the current knowledge about the post-bounce CCSN GW signal and recognize the templatable features that could be used by the ground-based laser interferometers. We demonstrate that starting from ~400ms after core bounce the dominant GW signal represents the fundamental quadrupole (l=2) oscillation mode (f-mode) of the proto-neutron star (PNS), which can be accurately reproduced by a linear perturbation analysis of the angle-averaged PNS profile. Before that, in the time interval between ~200 and ~400ms after bounce, the dominant mode has two radial nodes and represents a g-mode. We associate the high-frequency noise in the GW spectrograms above the main signal with p-modes, while below the dominant frequency there is a region with very little power. The collection of models presented here summarizes the dependence of the CCSN GW signal on the progenitor mass, equation of state, many-body corrections to the neutrino opacity, and rotation. Weak dependence of the dominant GW frequency on the progenitor mass motivates us to provide a simple fit for it as a function of time, which can be used as a prior when looking for CCSN candidates in the LIGO data.

A longstanding problem in the study of supernovae (SNe) has been the relationship between the Type IIP and Type IIL subclasses. Whether they come from distinct progenitors or they are from similar stars with some property that smoothly transitions from one class to another has been the subject of much debate. Here we show using one-dimensional radiation-hydrodynamic SN models that the multi-band light curves of SNe IIL are well fit by ordinary red supergiants surrounded by dense circumstellar material (CSM). The inferred extent of this material, coupled with a typical wind velocity of ~10-100 km/s, suggests enhanced activity by these stars during the last ~months to ~years of their lives, which may be connected with advanced stages of nuclear burning. Furthermore, we find that even for more plateau-like SNe that dense CSM provides a better fit to the first ~20 days of their light curves, indicating that the presence of such material may be more widespread than previously appreciated. Here we choose to model the CSM with a wind-like density profile, but it is unclear whether this just generally represents some other mass distribution, such as a recent mass ejection, thick disk, or even inflated envelope material. Better understanding the exact geometry and density distribution of this material will be an important question for future studies.

During the first several days after explosion, Type Ia supernova light curves probe the outer layers of the exploding star and therefore provide important clues for identifying their progenitors. We investigate how both the shallow $^{56}$Ni distribution and the presence of circumstellar material shape these early light curves. This is performed using a series of numerical experiments with parameterized properties for systematic exploration. Although not all of the considered models may be realized in nature (and indeed there are arguments why some of them should not occur), the spirit of this work is to provide a broader exploration of the diversity of possibilities. We find that shallower $^{56}$Ni leads to steeper, bluer light curves. Differences in the shape of the rise can introduce errors in estimating the explosion time and thus impact efforts to infer upper limits on the progenitor or companion radius from a lack of observed shock cooling emission. Circumstellar material can lead to significant luminosity during the first few days, but its presence can be difficult to identify depending on the degree of nickel mixing. In some cases, the hot emission of circumstellar material may even lead to a signature similar to interaction with a companion, and thus in the future additional diagnostics should be gathered for properly assessing early light curves.

The early rise of Type IIP supernovae (SN IIP) provides important information for constraining the properties of their progenitors. This can in turn be compared to pre-explosion imaging constraints and stellar models to develop a more complete picture of how massive stars evolve and end their lives. Using the SuperNova Explosion Code (SNEC), we model the first 40 days of SNe IIP to better understand what constraints can be derived from their early light curves. We use two sets of red supergiant progenitor models with zero-age main sequence masses in the range between 9 Msol and 20 Msol. We find that the early properties of the light curve depend most sensitively on the radius of the progenitor, and thus provide a relation between the g-band rise time and the radius at the time of explosion. This relation will be useful for deriving constraints on progenitors from future observations, especially in cases where detailed modeling of the entire rise is not practical. When comparing to observed rise times, the radii we find are a factor of a few larger than previous semi-analytic derivations and generally in better agreement with what is found with current stellar evolution calculations.

We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 Msun (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ~20-100 days if >~1.5-2 Msun of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses <~ 20 Msun, hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ~500 Rsun.

We revisit particle acceleration in the polar cap region of a neutron star by taking into account both general relativistic effects and the presence of toroidal oscillations at the star surface. In particular, we address the question of whether toroidal oscillations at the stellar surface can affect the acceleration properties in the polar cap. We solve numerically the relativistic electrodynamics equations in the stationary regime, focusing on the computation of the Lorentz factor of a space-charge-limited electron flow accelerated in the polar cap region of a rotating and oscillating pulsar. To this extent, we adopt the correct expression of the general relativistic Goldreich-Julian charge density in the presence of toroidal oscillations. Depending on the ratio of the actual charge density of the pulsar magnetosphere to the Goldreich-Julian charge density, we distinguish two different regimes of the Lorentz factor of the particle flow, namely an oscillatory regime produced for sub-GJ current density configurations, which does not produce an efficient acceleration, and a true accelerating regime for super-GJ current density configurations. We find that star oscillations may be responsible for a significant asymmetry in the pulse profile that depends on the orientation of the oscillations with respect to the pulsar magnetic field. In particular, significant enhancements of the Lorentz factor are produced by stellar oscillations in the super-GJ current density regime.

Using simple arguments based on photometric light curves and velocity evolution, we propose that some stripped envelope supernovae (SNe) show signs that a significant fraction of their helium is effectively transparent. The main pieces of evidence are the relatively low velocities with little velocity evolution, as are expected deep inside an exploding star, along with temperatures that are too low to ionize helium. This means that the helium should not contribute to the shaping of the main SN light curve, and thus the total helium mass may be difficult to measure from simple light curve modeling. Conversely, such modeling may be more useful for constraining the mass of the carbon/oxygen core of the SN progenitor. Other stripped envelope SNe show higher velocities and larger velocity gradients, which require an additional opacity source (perhaps the mixing of heavier elements or radioactive nickel) to prevent the helium from being transparent. We discuss ways in which similar analysis can provide insights into the differences and similarities between SNe Ib and Ic, which will lead to a better understanding of their respective formation mechanisms.

We try to explain the subpulse drift phenomena adopting the space-charge limited flow (SCLF) model and comparing the plasma drift velocity in the inner region of pulsar magnetospheres with the observed velocity of drifting subpulses. We apply the approach described in a recent paper of van Leeuwen & Timokhin (2012), where it was shown that the standard estimation of the subpulse drift velocity through the total value of the scalar potential drop in the inner gap gives inaccurate results, while the exact expression relating the drift velocity to the gradient of the scalar potential should be used instead. After considering a selected sample of sources taken from the catalog of Weltevrede, Edwards & Stappers (2006) with coherently drifting subpulses and reasonably known observing geometry, we show that their subpulse drift velocities would correspond to the drift of the plasma located very close or above the pair formation front. Moreover, a detailed analysis of PSR B0826-34 and PSR B0818-41 reveals that the variation of the subpulse separation with the pulse longitude can be successfully explained by the dependence of the plasma drift velocity on the angular coordinates.

We derive an analytic phenomenological expression that predicts the final mass of the black-hole remnant resulting from the merger of a generic binary system of black holes on quasi-circular orbits. Besides recovering the correct test-particle limit for extreme mass-ratio binaries, our formula reproduces well the results of all the numerical-relativity simulations published so far, both when applied at separations of a few gravitational radii, and when applied at separations of tens of thousands of gravitational radii. These validations make our formula a useful tool in a variety of contexts ranging from gravitational-wave physics to cosmology. As representative examples, we first illustrate how it can be used to decrease the phase error of the effective-one-body waveforms during the ringdown phase. Second, we show that, when combined with the recently computed self-force correction to the binding energy of nonspinning black-hole binaries, it provides an estimate of the energy emitted during the merger and ringdown. Finally, we use it to calculate the energy radiated in gravitational waves by massive black-hole binaries as a function of redshift, using different models for the seeds of the black-hole population.

We investigate the conditions for radio emission in rotating and oscillating magnetars, by focusing on the main physical processes determining the position of their death-lines in the P-\dot{P} diagram, i.e. of those lines that separate the regions where the neutron star may be radio-loud or radio-quiet. After using the general relativistic expression for the electromagnetic scalar potential in the magnetar magnetosphere, we find that larger compactness parameters of the star as well as larger inclination angles between the rotation axis and the magnetic moment produce death-lines well above the majority of known magnetars. This is consistent with the observational evidence of no regular radio emission from the magnetars in the frequency range typical for the ordinary pulsars. On the contrary, when oscillations of the magnetar are taken into account, the death-lines shift downward and the conditions necessary for the generation of radio emission in the magnetosphere are met. Present observations showing a close connection between the burst activity of magnetars and the generation of the radio emission in the magnetar magnetosphere are naturally accounted for within our interpretation.

We study the structure of electromagnetic field of slowly rotating magnetized star in a Randall-Sundrum II type braneworld. The star is modeled as a sphere consisting of perfect highly magnetized fluid with infinite conductivity and frozen-in dipolar magnetic field. Maxwell's equations for the external magnetic field of the star in the braneworld are analytically solved in approximation of small distance from the surface of the star. We have also found numerical solution for the electric field outside the rotating magnetized neutron star in the braneworld in dependence on brane tension. The influence of brane tension on the electromagnetic energy losses of the rotating magnetized star is underlined. Obtained "brane" corrections are shown to be relevant and have non-negligible values. In comparison with astrophysical observations on pulsars spindown data they may provide an evidence for the brane tension and, thus, serve as a test for the braneworld model of the Universe.

Plasma magnetosphere surrounding rotating magnetized neutron star in the braneworld has been studied. For the simplicity of calculations Goldreich-Julian charge density is analyzed for the aligned neutron star with zero inclination between magnetic field and rotation axis. From the system of Maxwell equations in spacetime of slowly rotating star in braneworld, second-order differential equation for electrostatic potential is derived. Analytical solution of this equation indicates the general relativistic modification of an accelerating electric field and charge density along the open field lines by brane tension. The implication of this effect to the magnetospheric energy loss problem is underlined. It was found that for initially zero potential and field on the surface of a neutron star, the amplitude of the plasma mode created by Goldreich-Julian charge density will increase in the presence of the negative brane charge. Finally we derive the equations of motion of test particles in magnetosphere of slowly rotating star in the braneworld. Then we analyze particle motion in the polar cap and show that brane tension can significantly change conditions for particle acceleration in the polar cap region of the neutron star.

We study the magnetosphere of a slowly rotating magnetized neutron star subject to toroidal oscillations in the relativistic regime. Under the assumption of a zero inclination angle between the magnetic moment and the angular momentum of the star, we analyze the Goldreich-Julian charge density and derive a second-order differential equation for the electrostatic potential. The analytical solution of this equation in the polar cap region of the magnetosphere shows the modification induced by stellar toroidal oscillations on the accelerating electric field and on the charge density. We also find that, after decomposing the oscillation velocity in terms of spherical harmonics, the first few modes with $m=0,1$ are responsible for energy losses that are almost linearly dependent on the amplitude of the oscillation and that, for the mode $(l,m)=(2,1)$, can be a factor $\sim8$ larger than the rotational energy losses, even for a velocity oscillation amplitude at the star surface as small as $\eta=0.05 \ \Omega \ R$. The results obtained in this paper clarify the extent to which stellar oscillations are reflected in the time variation of the physical properties at the surface of the rotating neutron star, mainly by showing the existence of a relation between $P\dot{P}$ and the oscillation amplitude. Finally, we propose a qualitative model for the explanation of the phenomenology of intermittent pulsars in terms of stellar oscillations that are periodically excited by star glitches.

The notion of death line of rotating pulsars is applied to model of oscillating neutron stars. It is shown that the magnetosphere of typical non-rotating oscillating stars may not contain secondary plasma to support the generation of radio emission in the region of open field lines of plasma magnetosphere.

We study magnetospheric structure surrounding rotating magnetized neutron star with nonvanishing NUT (Newman-Tamburino-Unti) parameter. For the simplicity of calculations Goldreich-Julian charge density is analyzed for the aligned neutron star with zero inclination between magnetic field, gravitomagnetic field and rotation axis. From the system of Maxwell equations in spacetime of slowly rotating NUT star, second-order differential equation for electrostatic potential is derived. Analytical solution of this equation indicates the general relativistic modification of an accelerating electric field and charge density along the open field lines by the gravitomagnetic charge. The implication of this effect to the magnetospheric energy loss problem is underlined. In the second part of the paper we derive the equations of motion of test particles in magnetosphere of slowly rotating NUT star. Then we analyze particle motion in the polar cap and show that NUT parameter can significantly change conditions for particle acceleration.

General relativistic quantum interference effects in the slowly rotating NUT space-time as the Sagnac effect and the phase shift effect of interfering particle in neutron interferometer are considered. It was found that in the case of the Sagnac effect the influence of NUT parameter is becoming important due to the fact that the angular velocity of the locally non rotating observer must be larger than one in the Kerr space-time. In the case of neutron interferometry it is found that due to the presence of NUT-parameter an additional term in the phase shift of interfering particle emerges. This term can be, in principle, detected by sensitive interferometer and derived results can be further used in experiments to detect the gravitomagnetic charge. Finally, as an example, we apply the obtained results to the calculation of the UCN (ultra-cold neutrons) energy level modification in the slowly rotating NUT space-time.