We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time.

Due to the high energies and long distances to the sources, astrophysical observations provide a unique opportunity to test possible signatures of Lorentz invariance violation (LIV). Superluminal LIV enables the decay of photons at high energy. The High Altitude Water Cherenkov (HAWC) Observatory is among the most sensitive gamma-ray instruments currently operating above 10 TeV. HAWC finds evidence of 100 TeV photon emission from at least four astrophysical sources. These observations exclude, for the strongest of the limits set, the LIV energy scale to $2.2\times10^{31}$ eV, over 1800 times the Planck energy and an improvement of 1 to 2 orders of magnitude over previous limits.

We present the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy gamma-ray source catalog to date. We report the integral flux of each of these objects. We also report spectra for three highest-energy sources and discuss the possibility that they are PeVatrons.

Twenty-four years after the discoveries of the first exoplanets, the radial-velocity (RV) method is still one of the most productive techniques to detect and confirm exoplanets. But stellar magnetic activity can induce RV variations large enough to make it difficult to disentangle planet signals from the stellar noise. In this context, HD41248 is an interesting planet-host candidate, with RV observations plagued by activity-induced signals. We report on ESPRESSO observations of HD41248 and analyse them together with previous observations from HARPS with the goal of evaluating the presence of orbiting planets. Using different noise models within a general Bayesian framework designed for planet detection in RV data, we test the significance of the various signals present in the HD41248 dataset. We use Gaussian processes as well as a first-order moving average component to try to correct for activity-induced signals. At the same time, we analyse photometry from the TESS mission, searching for transits and rotational modulation in the light curve. The number of significantly detected Keplerian signals depends on the noise model employed, which can range from 0 with the Gaussian process model to 3 with a white noise model. We find that the Gaussian process alone can explain the RV data while allowing for the stellar rotation period and active region evolution timescale to be constrained. The rotation period estimated from the RVs agrees with the value determined from the TESS light curve. Based on the data that is currently available, we conclude that the RV variations of HD41248 can be explained by stellar activity (using the Gaussian process model) in line with the evidence from activity indicators and the TESS photometry.

We present TeV gamma-ray observations of the Crab Nebula, the standard reference source in ground-based gamma-ray astronomy, using data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory. In this analysis we use two independent energy-estimation methods that utilize extensive air shower variables such as the core position, shower angle, and shower lateral energy distribution. In contrast, the previously published HAWC energy spectrum roughly estimated the shower energy with only the number of photomultipliers triggered. This new methodology yields a much improved energy resolution over the previous analysis and extends HAWC's ability to accurately measure gamma-ray energies well beyond 100 TeV. The energy spectrum of the Crab Nebula is well fit to a log parabola shape $\left(\frac{dN}{dE} = \phi_0 \left(E/\textrm{7 TeV}\right)^{-\alpha-\beta\ln\left(E/\textrm{7 TeV}\right)}\right)$ with emission up to at least 100 TeV. For the first estimator, a ground parameter that utilizes fits to the lateral distribution function to measure the charge density 40 meters from the shower axis, the best-fit values are $\phi_o$=(2.35$\pm$0.04$^{+0.20}_{-0.21}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.79$\pm$0.02$^{+0.01}_{-0.03}$, and $\beta$=0.10$\pm$0.01$^{+0.01}_{-0.03}$. For the second estimator, a neural network which uses the charge distribution in annuli around the core and other variables, these values are $\phi_o$=(2.31$\pm$0.02$^{+0.32}_{-0.17}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.73$\pm$0.02$^{+0.03}_{-0.02}$, and $\beta$=0.06$\pm$0.01$\pm$0.02. The first set of uncertainties are statistical; the second set are systematic. Both methods yield compatible results. These measurements are the highest-energy observation of a gamma-ray source to date.

List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

Detached WD+MS PCEBs are perhaps the most suitable objects for testing predictions of close-compact binary-star evolution theories, in particular, CE evolution. The population of WD+MS PCEBs has been simulated by several authors in the past and compared with observations. However, most of those predictions did not take the possible contributions to the envelope ejection from additional sources of energy (mostly recombination energy) into account. Here we update existing binary population models of WD+MS PCEBs by assuming that a fraction of the recombination energy available within the envelope contributes to ejecting the envelope. We performed Monte Carlo simulations of 10^7 MS+MS binaries for 9 different models using standard assumptions for the initial primary mass function, binary separations, and initial-mass-ratio distribution and evolved these systems using the publicly available BSE code. Including a fraction of recombination energy leads to a clear prediction of a large number of long orbital period (>~10 days) systems mostly containing high-mass WDs. The fraction of systems with He-core WD primaries increases with the CE efficiency and the existence of very low-mass He WDs is only predicted for high values of the CE efficiency (>~0.5). All models predict on average longer orbital periods for PCEBs containing C/O-core WDs than for PCEBs containing He WDs. This effect increases with increasing values of both efficiencies. Longer periods after the CE phase are also predicted for systems containing more massive secondary stars. The initial-mass-ratio distribution affects the distribution of orbital periods, especially the distribution of secondary star masses. Our simulations, in combination with a large and homogeneous observational sample, can provide constraints on the values of the CE efficiencies, as well as on the initial-mass-ratio distribution for MS+MS binary stars.

Detached white dwarf + main sequence (WD+MS) systems represent the simplest population of post-common envelope binaries (PCEBs). Since the ensemble properties of this population carries important information about the characteristics of the common-envelope (CE) phase, it deserves close scrutiny. However, most population synthesis studies do not fully take into account the effects of the observational selection biases of the samples used to compare with the theoretical simulations. Here we present the results of a set of detailed Monte Carlo simulations of the population of WD+MS binaries in the Sloan Digital Sky Survey (SDSS) Data Release 7. We used up-to-date stellar evolutionary models, a complete treatment of the Roche lobe overflow episode, and a full implementation of the orbital evolution of the binary systems. Moreover, in our treatment we took into account the selection criteria and all the known observational biases. Our population synthesis study allowed us to make a meaningful comparison with the available observational data. In particular, we examined the CE efficiency, the possible contribution of internal energy, and the initial mass ratio distribution (IMRD) of the binary systems. We found that our simulations correctly reproduce the properties of the observed distribution of WD+MS PCEBs. In particular, we found that once the observational biases are carefully taken into account, the distribution of orbital periods and of masses of the WD and MS stars can be correctly reproduced for several choices of the free parameters and different IMRDs, although models in which a moderate fraction (<=10%) of the internal energy is used to eject the CE and in which a low value of CE efficiency is used (<=0.3) seem to fit better the observational data. We also found that systems with He-core WDs are over-represented in the observed sample, due to selection effects.

We argue that the multiple shells of circumstellar material (CSM) and the supernovae (SN) ejecta interaction with the CSM starting 59 days after the explosion of the Type Ia SN (SN Ia) PTF 11kx, are best described by a violent prompt merger. In this prompt merger scenario the common envelope (CE) phase is terminated by a merger of a WD companion with the hot core of a massive asymptotic giant (AGB) star. In most cases the WD is disrupted and accreted onto the more massive core. However, in the rare cases where the merger takes place when the WD is denser than the core, the core will be disrupted and accreted onto the cooler WD. In such cases the explosion might occur with no appreciable delay, i.e., months to years after the termination of the CE phase. This, we propose, might be the evolutionary route that could lead to the explosion of PTF 11kx. This scenario can account for the very massive CSM within ~1000 AU of the exploding PTF 11kx star, for the presence of hydrogen, and for the presence of shells in the CSM.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the solar neighborhood.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the Solar neighborhood.

NGC 6791 is a well-studied, metal-rich open cluster that is so close to us that can be imaged down to luminosities fainter than that of the termination of its white dwarf cooling sequence, thus allowing for an in-depth study of its white dwarf population. We use a Monte Carlo simulator that employs up-to-date evolutionary cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres, with carbon-oxygen and helium cores. The cooling sequences for carbon-oxygen cores account for the delays introduced by both Ne^22 sedimentation in the liquid phase and by carbon-oxygen phase separation upon crystallization. We do not find evidence for a substantial fraction of helium-core white dwarfs, and hence our results support the suggestion that the origin of the bright peak of the white dwarf luminosity function can only be attributed to a population of unresolved binary white dwarfs. Moreover, our results indicate that the number distribution of secondary masses of the population of unresolved binaries has to increase with increasing mass ratio between the secondary and primary components of the progenitor system. We also find that the observed cooling sequence appears to be able to constrain the presence of progenitor sub-populations with different chemical compositions and the fraction of non-DA white dwarfs. Our simulations place interesting constraints on important characteristics of the stellar populations of NGC 6791. In particular, we find that the fraction of single helium-core white dwarfs must be smaller than 5%, that a sub-population of stars with zero metallicity must be <12%, while if the adopted metallicity of the sub-population is solar the upper limit is ~8%. Finally, we also find that the fraction of non-DA white dwarfs in this particular cluster is surprinsingly small <6%.

The nature of the several microlensing events observed by the MACHO team towards the Large Magellanic Cloud (LMC) is still a subject of debate. Low-mass substellar objects and stars with masses larger than ~M_{sun} have been ruled out as major components of a Massive Astrophysical Compact Halo Object (MACHO) Galactic halo, while stars of half a solar mass seem to be viable candidates. Main sequence stars have been already discarded, and there are tight restrictions on the role played by white dwarfs with hydrogen-dominated atmospheres. In this paper we evaluate the contribution to the dark matter content of the Galaxy of white dwarfs with hydrogen-deficient atmospheres. For this purpose we use a Monte Carlo simulator which incorporates up-to-date evolutionary sequences of white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres. We also take into account detailed descriptions of the thick disk and the halo of our Galaxy as well as of a reliable model of the LMC. We find that the contribution of white dwarfs with hydrogen-deficient atmospheres moderately increases the theoretical estimate of the optical depth with respect to the value obtained when only hydrogen-rich white dwarfs are considered. We also find that the contribuiton of the thick disk population of white dwarfs is comparable to the halo contribution. However, the contributions of both the halo and the thick disk white-dwarf populations are still insufficient to explain the number of events observed by the MACHO team. Finally, we find that the contribution to the halo dark matter of the entire population under study is less than 10% at the 95% conficence level.

The nature of the several microlensing events observed by the MACHO team towards the LMC still remains controversial. Low-mass substellar objects and stars with masses larger than ~1 M_{sun} have been ruled out as major components of a MACHO galactic halo, while stars of half solar masses are the most probable candidates. In this paper we assess jointly the relative contributions of both red dwarfs and white dwarfs to the mass budget of the galactic halo. In doing so we use a Monte Carlo simulator which incorporates up-to-date evolutionary sequences of both red dwarfs and white dwarfs as well as detailed descriptions of our Galaxy and of the LMC. We explore the complete mass range between 0.08 and 1 M_{sun} as possible microlensing candidates and we compare the synthetic populations obtained with our simulator with the results obtained by the MACHO and EROS experiments. Our results indicate that, despite that the contribution of the red dwarf population increases by a factor of 2 the value of the optical depth obtained when taking into account the white dwarf population alone, it is still insufficient to explain the number of events claimed by the MACHO team. Finally, we find that the contribution to the halo dark matter of the whole population under study is smaller than 10% at the 95% conficence level.

The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains controversial. White dwarfs have been proposed to explain these results and, hence, to contribute significantly to the mass budget of our Galaxy. However, several constraints on the role played by regular carbon-oxygen white dwarfs exist. Massivewhite dwarfs are thought to be made of a mixture of oxygen and neon. Correspondingly, their cooling rate is larger than those of typical carbon-oxygen white dwarfs and they fade to invisibility in short timescales. Consequently, they constitute a good candidate for explaining the microlensing results. Here, we examine in detail this hypothesis by using the most recent and up-to-date cooling tracks for massive white dwarfs and a Monte Carlo simulator which takes into account the most relevant Galactic inputs. We find that oxygen-neon white dwarfs cannot account for a substantial fraction of the microlensing depth towards the LMC, independently of the adopted initial mass function, although some microlensing events could be due to oxygen--neon white dwarfs. The white dwarf population contributes at most a 5% to the mass of the Galactic halo.