The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

During the first few hundred days after the explosion, core-collapse supernovae (SNe) emit down-scattered X-rays and gamma-rays originating from radioactive line emissions, primarily from the $^{56}$Ni $\rightarrow$ $^{56}$Co $\rightarrow$ $^{56}$Fe chain. We use SN models based on three-dimensional neutrino-driven explosion simulations of single stars and mergers to compute this emission and compare the predictions with observations of SN 1987A. A number of models are clearly excluded, showing that high-energy emission is a powerful way of discriminating between models. The best models are almost consistent with the observations, but differences that cannot be matched by a suitable choice of viewing angle are evident. Therefore, our self-consistent models suggest that neutrino-driven explosions are able to produce, in principle, sufficient mixing, although remaining discrepancies may require small changes to the progenitor structures. The soft X-ray cutoff is primarily determined by the metallicity of the progenitor envelope. The main effect of asymmetries is to vary the flux level by a factor of ${\sim}$3. For the more asymmetric models, the shapes of the light curves also change. In addition to the models of SN 1987A, we investigate two models of Type II-P SNe and one model of a stripped-envelope Type IIb SN. The Type II-P models have similar observables as the models of SN 1987A, but the stripped-envelope SN model is significantly more luminous and evolves faster. Finally, we make simple predictions for future observations of nearby SNe.

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.

Based on the work of Menon & Heger (2017), we present the bolometric light curvesand spectra of the explosions of blue supergiant progenitors from binary mergers. We study SN 1987A and two other peculiar Type IIP supernovae: SN 1998A and SN 2006V. The progenitor models were produced using the stellar evolution code Keplerand then exploded using the 1D radiation hydrodynamic code Crab. The explosions of binary merger models exhibit an overall better fit to the light curve of SN 1987A than previous single star models, because of their lower helium-core masses, larger envelope masses, and smaller radii. The merger model that best matches the observational constraints of the progenitor of SN 1987A and the light curve is a model with a radius of 37 solar radii, an ejecta mass of 20.6 solar masses, an explosion energy of 1.7 Bethe, a nickel mass of 0.073 solar masses, and a nickel mixing velocity of 3,000 km/s. This model also works for SN 1998A and is comparable with earlier estimates from semi-analytic models. In the case of SN 2006V, however, a model with a radius of 150 solar radii and ejecta mass of 19.1 solar masses matches the light curve. These parameters are significantly higher than predictions from semi-analytic models for the progenitor of this supernova.

We present the results of a detailed, systematic stellar evolution study of binary mergers for blue supergiant (BSG) progenitors of Type II supernovae. In particular, these are the first evolutionary models that can simultaneously reproduce nearly all observational aspects of the progenitor of SN 1987A, $\text{Sk}-69\,^{\circ}202$, such as its position in the HR diagram, the enrichment of helium and nitrogen in the triple-ring nebula, and its lifetime before its explosion. The merger model, based on the one proposed by Podsiadlowski 1992 et al. and Podsiadlowski 2007 et al., consists of a main sequence secondary star that dissolves completely in the common envelope of the primary red supergiant at the end of their merger. We empirically explore a large initial parameter space, such as primary masses ($15\,\text{M}_{\odot}$, $16\,\text{M}_{\odot}$, and $17\,\text{M}_{\odot}$), secondary masses ($2\,\text{M}_{\odot}$, $3\,\text{M}_{\odot}$, ..., $8\,\text{M}_{\odot}$) and different depths up to which the secondary penetrates the He core of the primary during the merger. The evolution of the merged star is continued until just before iron-core collapse and the surface properties of the 84 pre-supernova models ($16\,\text{M}_{\odot}-23\,\mathrm{M}_{\odot}$) computed have been made available in this work. Within the parameter space studied, the majority of the pre-supernova models are compact, hot BSGs with effective temperature $>12\,\text{kK}$ and radii of $30\,\text{R}_{\odot}-70\,\mathrm{R}_{\odot}$ of which six match nearly all the observational properties of $\text{Sk}-69\,^{\circ}202$.

We study the viability of regions of large $\tan \beta$ within the framework of Fat Higgs/$\lambda$-SUSY Models. We compute the one-loop effective potential to find the corrections to the Higgs boson mass due to the heavy non-standard Higgs bosons. As the tree level contribution to the Higgs boson mass is suppressed at large $\tan \beta$, these one-loop corrections are crucial to raising the Higgs boson mass to the measured LHC value. By raising the Higgsino and singlino mass parameters, typical electroweak precision constraints can also be avoided. We illustrate these new regions of Fat Higgs/$\lambda$-SUSY parameter space by finding regions of large $\tan \beta$ that are consistent with all experimental constraints including direct dark matter detection experiments, relic density limits and the invisible decay width of the $Z$ boson. We find that there exist regions around $\lambda = 1.25, \tan \beta = 50$ and a uniform psuedo-scalar $4~{\rm TeV} \lsim M_A \lsim 8$~TeV which are consistent will all present phenomenological constraints. In this region the dark matter relic abundance and direct detection limits are satisfied by a lightest neutralino that is mostly bino or singlino. As an interesting aside we also find a region of low $\tan \beta$ and small singlino mass parameter where a well-tempered neutralino avoids all cosmological and direct detection constraints.

The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars that are most likely the result of He-CO WD mergers. They display extremely low oxygen isotopic ratios, 16O/18O ~ 1 - 10, 12C/13C>=100, and enhancements up to 2.6dex in F and in s-process elements from Zn to La, compared to solar. These abundances provide stringent constraints on the physical processes during and after the double-degenerate merger. As shown before O-isotopic ratios observed in RCB stars cannot result from the dynamic double-degenerate merger phase, and we investigate now the role of the long-term 1D spherical post-merger evolution and nucleosynthesis based on realistic hydrodynamic merger progenitor models. We adopt a model for extra envelope mixing to represent processes driven by rotation originating in the dynamical merger. Comprehensive nucleosynthesis post-processing simulations for these stellar evolution models reproduce, for the first time, the full range of the observed abundances for almost all the elements measured in RCB stars: 16O/18O ratios between 9 and 15, C-isotopic ratios above 100, and ~1.4 - 2.35dex F enhancements, along with enrichments in s-process elements. The nucleosynthesis processes in our models constrain the length and temperature in the dynamic merger shell-of-fire feature as well as the envelope mixing in the post-merger phase. s-process elements originate either in the shell-of-fire merger feature or during the post-merger evolution, but the contribution from the AGB progenitors is negligible. The post-merger envelope mixing must eventually cease ~ 10^6yr after the dynamic merger phase, before the star enters the RCB phase.

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

We consider a new scenario for supersymmetric decaying dark matter without R-parity violation in theories with goldstini, which arise if supersymmetry is broken independently by multiple sequestered sectors. The uneaten goldstino naturally has a long lifetime and decays into three-body final states including the gravitino, which escapes detection, and two visible particles. The goldstini low-energy effective interactions are derived, which can be non-universal and allow the dark matter to be leptophilic, in contrast to the case of a single sector supersymmetry breaking. In addition, the three-body decay with a missing particle gives a softer spectrum. Consequently, it is possible to fit both the positron excess observed by the PAMELA and the total e+ + e- measurements by the Fermi-LAT using universal couplings to all three lepton flavors or 100% branching fraction into electrons/positrons, both of which are disfavored in the conventional scenario of dark matter decays into two or four visible particles without missing energy.

In models with one universal extra dimension (UED), the first Kaluza-Klein excitations of the hypercharge gauge boson, B^(1), and the neutral component of isospin gauge boson, W^3(1), are each viable dark matter candidates. In either case, such particles are predicted to accumulate in the core of the Sun, where they annihilate to generate a potentially observable flux of high energy neutrinos. In this article, we calculate the flux of neutrinos produced in this model and determine the constraints that can be placed on the UED parameter space from current IceCube data. For the case of B^(1) dark matter, we find that the present limits from IceCube are stronger than those from direct dark matter detection experiments such as CDMS and XENON10. For W^3(1) dark matter, the present IceCube data provides a constraint slightly weaker than direct detection experiments. In addition, we also present the projected regions of UED parameter space that can be probed by IceCube/DeepCore in the near future and compare them to the prospects for future direct detection experiments.

In this paper we investigate the phenomenological implications of boundary localized terms (BLTs) in the model of Universal Extra Dimensions (UED). In particular, we study the electroweak Kaluza-Klein mass spectrum resulting from BLTs and their effect on electroweak symmetry breaking via the five dimensional Higgs mechanism. We find that the addition of BLTs to massive five dimensional fields induces a non-trivial extra dimensional profile for the zero and non-zero Kaluza-Klein (KK) modes. Hence BLTs generically lead to a modification of Standard Model parameters and are therefore experimentally constrained, even at tree level. We study Standard Model constraints on three representative non-minimal UED models in detail and find that the constraints on BLTs are weak. On the contrary, non-zero BLTs have a major impact on the spectrum and couplings of non-zero KK modes. For example, there are regions of parameter space where the Lightest Kaluza-Klein particle (LKP) is either the Kaluza-Klein Higgs boson or the first KK mode of the W^3.