Context: It has been suggested that some supernovae (SNe) may be powered by a magnetar formed at the moment of the explosion. While this scenario has mostly been applied to hydrogen-free events, it may be possible also for hydrogen-rich objects. Aims: We explore the effect of including a magnetar on the light curves of supernovae with H-rich progenitors. Methods: We have applied a version of our one-dimensional LTE radiation hydrodynamics code that takes into account the relativistic motion of the ejecta caused by the extra energy provided by the magnetar. For a fixed red supergiant (RSG) progenitor, we have obtained a set of light curves that corresponds to different values of the magnetar initial rotation energy and the spin-down timescale. The model is applied to SN~2004em and OGLE-2014-SN-073, two peculiar Type II SNe with long-rising SN1987A-like light curves, although with much larger luminosities. Results: The presence of a plateau phase in either normal or superluminous supernovae is one possible outcome, even if a magnetar is continuously injecting energy into the ejecta.In other cases, the light curve shows a peak but not a plateau. Also, there are intermediate events with a first peak followed by a slow decline and a late break of the declining slope. Our models show that bright and long rising morphologies are possible even assuming RSG structures. Conclusions: A large number of supernova discoveries per year reveal unexpected new types of explosions. According to our results, SLSNe II-P are to be expected, as well as a variety of light curve morphologies that can all be possible signs of a newly born magnetar.

It is difficult to establish the properties of massive stars that explode as supernovae. The electromagnetic emission during the first minutes to hours after the emergence of the shock from the stellar surface conveys important information about the final evolution and structure of the exploding star. However, the unpredictable nature of supernova events hinders the detection of this brief initial phase. Here we report the serendipitous discovery of a newly born, normal type IIb supernova (SN 2016gkg), which reveals a rapid brightening at optical wavelengths of about 40 magnitudes per day. The very frequent sampling of the observations allowed us to study in detail the outermost structure of the progenitor of the supernova and the physics of the emergence of the shock. We develop hydrodynamical models of the explosion that naturally account for the complete evolution of the supernova over distinct phases regulated by different physical processes. This result suggests that it is appropriate to decouple the treatment of the shock propagation from the unknown mechanism that triggers the explosion.

We analyse configurations of compact stars in the so-called R-squared gravity in the Palatini formalism. Using a realistic equation of state we show that the mass-radius configurations are lighter than their counterparts in General Relativity. We also obtain the internal profiles, which run in strong correlation with the derivatives of the equation of state, leading to regions where the mass parameter decreases with the radial coordinate in a counter-intuitive way. In order to analyse such correlation, we introduce a parametrisation of the equation of state given by multiple polytropes, which allows us to explicitly control its derivatives. We show that, even in a limiting case where hard phase transitions in matter are allowed, the internal profile of the mass parameter still presents strange features and the calculated M-R configurations also yield NSs lighter than those obtained in General Relativity.

Two recently discovered very luminous supernovae (SNe) present stimulating cases to explore the extents of the available theoretical models. SN 2011kl represents the first detection of a supernova explosion associated with an ultra-long duration gamma ray burst. ASASSN-15lh was even claimed as the most luminous SN ever discovered, challenging the scenarios so far proposed for stellar explosions. Here we use our radiation hydrodynamics code in order to simulate magnetar powered SNe. To avoid explicitly assuming neutron star properties we adopt the magnetar luminosity and spin-down timescale as free parameters of the model. We find that the light curve (LC) of SN 2011kl is consistent with a magnetar power source, as previously proposed, but we note that some amount of 56^Ni (> 0.08 M_sun) is necessary to explain the low contrast between the LC peak and tail. For the case of ASASSN-15lh we find physically plausible magnetar parameters that reproduce the overall shape of the LC provided the progenitor mass is relatively large (a mass of the ejecta approx 6 M_sun). The ejecta hydrodynamics of this event is dominated by the magnetar input, while the effect is more moderate for SN 2011kl. We conclude that a magnetar model may be used for the interpretation of these events and that the hydrodynamic modeling is necessary to derive the properties of powerful magnetars and their progenitors.

We investigate the plausibility of detecting X-ray emission from a stellar jet that impacts against a dense molecular cloud. This scenario may be usual for classical T Tauri stars with jets in dense star-forming complexes. We first model the impact of a jet against a dense cloud by 2D axisymmetric hydrodynamic simulations, exploring different configurations of the ambient environment. Then, we compare our results with XMM-Newton observations of the Herbig-Haro object HH 248, where extended X-ray emission aligned with the optical knots is detected at the edge of the nearby IC 434 cloud. Our simulations show that a jet can produce plasma with temperatures up to 10 MK, consistent with production of X-ray emission, after impacting a dense cloud. We find that jets denser than the ambient medium but less dense than the cloud produce detectable X-ray emission only at the impact onto the cloud. From the exploration of the model parameter space, we constrain the physical conditions (jet density and velocity, cloud density) that reproduce well the intrinsic luminosity and emission measure of the X-ray source possibly associated with HH 248. Thus, we suggest that the extended X-ray source close to HH 248 corresponds to the jet impacting on a dense cloud.

The effects implied for the structure of compact objects by the modification of General Relativity produced by the generalization of the Lagrangian density to the form f(R)=R+\alpha R^2, where R is the Ricci curvature scalar, have been recently explored. It seems likely that this squared-gravity may allow heavier Neutron Stars (NSs) than GR. In addition, these objects can be useful to constrain free parameters of modified-gravity theories. The differences between alternative gravity theories is enhanced in the strong gravitational regime. In this regime, because of the complexity of the field equations, perturbative methods become a good choice to treat the problem. Following previous works in the field, we performed a numerical integration of the structure equations that describe NSs in f(R)-gravity, recovering their mass-radius relations, but focusing on particular features that arise from this approach in the profiles of the NS interior. We show that these profiles run in correlation with the second-order derivative of the analytic approximation to the Equation of State (EoS), which leads to regions where the enclosed mass decreases with the radius in a counter-intuitive way. We reproduce all computations with a simple polytropic EoS to separate zeroth-order modified gravity effects.

We present continuum high resolution Submillimeter Array (SMA) observations of the transition disk object RX J1633.9-2442, which is located in the Ophiuchus molecular cloud and has recently been identified as a likely site of ongoing giant planet formation. The observations were taken at 340 GHz (880 micron) with the SMA in its most extended configuration, resulting in an angular resolution of 0.3" (35 AU at the distance of the target). We find that the disk is highly inclined (i ~50 deg) and has an inner cavity ~25 AU in radius, which is clearly resolved by our observations. We simultaneously model the entire optical to millimeter wavelength spectral energy distribution (SED) and SMA visibilities of RX J1633.9-2442 in order to constrain the structure of its disk. We find that an empty cavity ~25 AU in radius is inconsistent with the excess emission observed at 12, 22, and 24 micron. Instead, the mid-IR excess can be modeled by either a narrow, optically thick ring at ~10 AU or an optically thin region extending from ~7 AU to ~25 AU. The inner disk (r < 5 AU) is mostly depleted of small dust grains as attested by the lack of detectable near-IR excess. We also present deep Keck aperture masking observations in the near-IR, which rule out the presence of a companion up to 500 times fainter than the primary star (in K-band) for projected separations in the 5-20 AU range. We argue that the complex structure of the RX J1633.9-2442 disk is best explained by multiple planets embedded within the disk. We also suggest that the properties and incidence of objects such as RX J1633.9-2442, T Cha, and LkCa 15 (and those of the companions recently identified to these two latter objects) are most consistent with the runaway gas accretion phase of the core accretion model, when giant planets gain their envelopes and suddenly become massive enough to open wide gaps in the disk.

A large set of Spitzer-selected transitional disks in the Ophiuchus molecular cloud was recently examined by Cieza et al (2010), and 4 of the targets were identified as (giant) planet-forming candidates based on the morphology of their Spectral Energy Distributions (SEDs), the apparent lack of stellar companions, and evidence for accretion. Here we characterize the structures of these disks modeling their optical, infrared and (sub)millimeter SEDs. We use the Monte Carlo radiative transfer package RADMC to construct a parametric model of the dust distribution in a flared disk with an inner cavity and calculate the temperature structure consistent with the density profile, in thermal equilibrium with the irradiating star. For each object, we conducted a Bayesian exploration of the parameter space generating Monte Carlo Markov Chains (MCMC) that allow the identification of the best fit parameters and to constrain their range of statistical confidence. Our calculations point to the presence of evacuated cavities with radii ~2-8 AU, consistent with having been carved by embedded giant planets. We found parameter values consistent with those previously given in the literature, indicating a mild degree of grain growth and dust settling, which deserves to be investigated with further modeling and follow up observations. Resolved images with (sub)millimeter interferometers are required to break some of the degeneracies of the models and better constrain the physical properties of these fascinating disks.

The high-mass microquasar Cygnus X-1, the best-established candidate for a stellar-mass black hole in the Galaxy, has been detected in a flaring state at very high energies (VHE), E > 200 GeV, by the Atmospheric Cherenkov Telescope MAGIC. The flare occurred at orbital phase 0.91, where phase 1 is the configuration with the black hole behind the companion high-mass star, when the absorption of gamma-ray photons by photon-photon annihilation with the stellar field is expected to be highest. We aim to set up a model for the high-energy emission and absorption in Cyg X-1 that can explain the nature of the observed gamma-ray flare. We study the gamma-ray opacity due to pair creation along the whole orbit, and for different locations of the emitter. Then we consider a possible mechanism for the production of the VHE emission. We present detailed calculations of the gamma-ray opacity and infer from these calculations the distance from the black hole where the emitting region was located. We suggest that the flare was the result of a jet-clump interaction where the decay products of inelastic proton-proton collisions dominate the VHE outcome. We are able to reproduce the spectrum of Cyg X-1 during the observed flare under reasonable assumptions. The flare may be the first event of jet-cloud interaction ever detected at such high energies.

We have obtained millimeter wavelength photometry, high-resolution optical spectroscopy and adaptive optics near-infrared imaging for a sample of 26 Spitzer-selected transition circumstellar disks. All of our targets are located in the Ophiuchus molecular cloud (d ~125 pc) and have Spectral Energy Distributions (SEDs) suggesting the presence of inner opacity holes. We use these ground-based data to estimate the disk mass, multiplicity, and accretion rate for each object in our sample in order to investigate the mechanisms potentially responsible for their inner holes. We find that transition disks are a heterogeneous group of objects, with disk masses ranging from < 0.6 to 40 M_JUP and accretion rates ranging from <1E-11 1E-7 M_solar/yr, but most tend to have much lower masses and accretion rates than "full disks" (i.e., disks without opacity holes). Eight of our targets have stellar companions: 6 of them are binaries and the other 2 are triple systems. In four cases, the stellar companions are close enough to suspect they are responsible for the inferred inner holes. We find that 9 of our 26 targets have low disk mass (< 2.5 M_JUP) and negligible accretion (< 1E-11) M_solar/yr, and are thus consistent with photoevaporating (or photoevaporated) disks. Four of these 9 non-accreting objects have fractional disk luminosities < 1E-3 and could already be in a debris disk stage. Seventeen of our transition disks are accreting. Thirteen of these accreting objects are consistent with grain growth. The remaining 4 accreting objects have SEDs suggesting the presence of sharp inner holes, and thus are excellent candidates for harboring giant planets.

Centaurus A, the closest active galaxy, has been detected from radio to high-energy gamma-rays. The synchrotron radiation by extremely high energy protons may be a suitable mechanism to explain the MeV to GeV emission detected by the instruments of the Compton Gamma-Ray Observatory, as coming from the inner jets. This scenario requires a relatively large magnetic field of about 10^4 G that could be present only close to the central black hole. We investigate the spectral energy distribution (SED) resulting from a one-zone compact acceleration region, where both leptonic and hadronic relativistic populations arise. We present here results of such a model, where we have considered synchrotron radiation by primary electrons and protons, inverse Compton interactions, and gamma-ray emission originated by the inelastic hadronic interactions between relativistic protons and cold nuclei within the jets themselves. Photo-meson production by relativistic hadrons were also taken into account, as well as the effects of secondary particles injected by all interactions. The internal and external absorption of gamma rays is shown to be of great relevance to shape the observable SED, which was also recently constrained by the results of HESS.

LS I +61 303 is a puzzling object detected from radio up to high-energy gamma-rays. Variability has recently been observed in its high-energy emission. The object is a binary system, with a compact object and a Be star as primary. The nature of the secondary and the origin of the gamma-ray emission are not clearly established at present. Recent VLBA radio data have been used to claim that the system is a Be/neutron star colliding wind binary, instead of a microquasar. We review the main views on the nature of LS I +61 303 and present results of 3D SPH simulations that can shed some light on the nature of the system. Our results support an accretion powered source, compatible with a microquasar interpretation.

Context: It has been proposed that the origin of the very high-energy photons emitted from high-mass X-ray binaries with jet-like features, so-called microquasars (MQs), is related to hadronic interactions between relativistic protons in the jet and cold protons of the stellar wind. Leptonic secondary emission should be calculated in a complete hadronic model that include the effects of pairs from charged pion decays inside the jets and the emission from pairs generated by gamma-ray absorption in the photosphere of the system. Aims: We aim at predicting the broadband spectrum from a general hadronic microquasar model, taking into account the emission from secondaries created by charged pion decay inside the jet. Methods: The particle energy distribution for secondary leptons injected along the jets is consistently derived taking the energy losses into account. We also compute the spectral energy distribution resulting from these leptons is calculated after assuming different values of the magnetic field inside the jets. The spectrum of the gamma-rays produced by neutral pion-decay and processed by electromagnetic cascades under the stellar photon field. Results: We show that the secondary emission can dominate the spectral energy distribution at low energies (~1 MeV). At high energies, the production spectrum can be significantly distorted by the effect of electromagnetic cascades. These effects are phase-dependent, and some variability modulated by the orbital period is predicted.

We present a model to estimate the synchrotron radio emission generated in microquasar (MQ) jets due to secondary pairs created via decay of charged pions produced in proton-proton collisions between stellar wind ions and jet relativistic protons. Signatures of electrons/positrons are obtained from consistent particle energy distributions that take into account energy losses due to synchrotron and inverse Compton (IC) processes, as well as adiabatic expansion. The space parameter for the model is explored and the corresponding spectral energy distributions (SEDs) are presented. We conclude that secondary leptonic emission represents a significant though hardly dominant contribution to the total radio emission in MQs, with observational consequences that can be used to test some still unknown processes occurring in these objects as well as the nature of the matter outflowing in their jets.

LS I +61 303 is a puzzling Be/X-ray binary with variable gamma-ray emission at up TeV energies. The nature of the compact object and the origin of the high-energy emission are unclear. One family of models invokes particle acceleration in shocks from the collision between the B-star wind and a relativistic pulsar wind, while another centers on a relativistic jet powered by accretion. Recent high-resolution radio observations showing a putative "cometary tail" pointing away from the Be star near periastron have been cited as support for the pulsar-wind model. We wish here to carry out a quantitative assessment of these competing models for this extraordinary source. We apply a 3D SPH code for dynamical simulations of both the pulsar-wind-interaction and accretion-jet models. The former yields a description of the shape of the wind-wind interaction surface. The latter provides an estimation of the accretion rate. The results allow critical evaluation of how the two distinct models confront the data in various wavebands under a range of conditions. When one accounts for the 3D dynamical wind interaction under realistic constraints for the relative strength of the B-star and pulsar winds, the resulting form of the interaction front does not match the putative "cometary tail" claimed from radio observations. On the other hand, dynamical simulations of the accretion-jet model indicate that the orbital phase variation of accretion power includes a secondary broad peak well away from periastron, thus providing a plausible way to explain the observed TeV gamma ray emission toward apastron. We conclude that the colliding-wind model is not clearly established for LS I +61 303, while the accretion-jet model can reproduce many key characteristics of the observed TeV gamma-ray emission.

Context: Some high mass X-ray binaries (HMXB) have been recently confirmed as gamma-ray sources by ground based Cherenkov telescopes. In this work, we discuss the gamma-ray emission from X-ray transient sources formed by a Be star and a highly magnetized neutron star. This kind of systems can produce variable hadronic gamma-ray emission through the mechanism proposed by Cheng and Ruderman, where a proton beam accelerated in the pulsar magnetosphere impacts the transient accretion disk. We choose as case of study the best known system of this class: A0535+26. Aims: We aim at making quantitative predictions about the very high-energy radiation generated in Be-X ray binary systems with strongly magnetized neutron stars. Methods: We study the gamma-ray emission generated during a major X-ray outburst of a HMXB adopting for the model the parameters of A0535+26. The emerging photon signal from the disk is determined by the grammage of the disk that modulates the optical depth. The electromagnetic cascades initiated by photons absorbed in the disk are explored, making use of the so-called "Approximation A" to solve the cascade equations. Very high energy photons induce Inverse Compton cascades in the photon field of the massive star. We implemented Monte Carlo simulations of these cascades, in order to estimate the characteristics of the resulting spectrum. Results: TeV emission should be detectable by Cherenkov telescopes during a major X-ray outburst of a binary formed by a Be star and a highly magnetized neutron star. The gamma-ray light curve is found to evolve in anti-correlation with the X-ray signal.

LS I+61 303 has been detected by the Cherenkov telescope MAGIC at very high energies, presenting a variable flux along the orbital motion with a maximum clearly separated from the periastron passage. In the light of the new observational constraints, we revisit the discussion of the production of high-energy gamma rays from particle interactions in the inner jet of this system. The hadronic contribution could represent a major fraction of the TeV emission detected from this source. The spectral energy distribution resulting from p-p interactions is recalculated. Opacity effects introduced by the photon fields of the primary star and the stellar decretion disk are shown to be essential in shaping the high-energy gamma-ray light curve at energies close to 200 GeV. We also present results of Monte Carlo simulations of the electromagnetic cascades developed very close to the periastron passage. We conclude that a hadronic microquasar model for the gamma-ray emission in LS I +61 303 can reproduce the main features of its observed high-energy gamma-ray flux.

We show that high-energy neutrinos can be efficiently produced in X-ray binaries with relativistic jets and high-mass primary stars. We consider a system where the star presents a dense equatorial wind and the jet has a small content of relativistic protons. In this scenario, neutrinos and correlated gamma-rays result from pp interactions and the subsequent pion decays. As a particular example we consider the microquasar LS I +61 303. Above 1 TeV, we obtain a mean-orbital neutrino muon induced luminosity of 5 10^34 erg/s which can be related to an event rate of 4-5 muon-type neutrinos per kilometer-squared per year after considering the signal attenuation due to maximal neutrino oscillations. The maximal neutrino energies here considered will range between 20 and 85 TeV along the orbit. The local infrared photon-field is responsible for opacity effects on the associated gamma radiation at high energies, but below 50 GeV the source could be detected by MAGIC telescope. GLAST observations at photon energies E>100 MeV should also reveal a strong source.

We report the discovery of a new O-type double-lined spectroscopic binary with a short orbital period of 1.4 days. We find the primary component of this binary, FO 15, to have an approximate spectral type O5.5Vz, i.e. a Zero-Age-Main-Sequence star. The secondary appears to be of spectral type O9.5V. We have performed a numerical model fit to the public ASAS photometry, which shows that FO 15 is also an eclipsing binary. We find an orbital inclination of ~ 80 deg. From a simultaneous light-curve and radial velocity solution we find the masses and radii of the two components to be 30 +/- 1 and 16 +/- 1 solar masses and 7.5 +/- 0.5 and 5.3 +/- 0.5 solar radii. These radii, and hence also the luminosities, are smaller than those of normal O-type stars, but similar to recently born ZAMS O-type stars. The absolute magnitudes derived from our analysis locate FO 15 at the same distance as Eta Carinae. From Chandra and XMM X-ray images we also find that there are two close X-ray sources, one coincident with FO 15 and another one without optical counterpart. This latter seems to be a highly variable source, presumably due to a pre-main-sequence stellar neighbour of FO 15.

We present a hadronic model for gamma-ray production in the microquasar LS I $+61 303$. The system is formed by a neutron star that accretes matter from the dense and slow equatorial wind of the Be primary star. We calculate the gamma-ray emission originated in $pp$ interactions between relativistic protons in the jet and cold protons from the wind. After taking into account opacity effects on the gamma-rays introduced by the different photons fields, we present high-energy spectral predictions that can be tested with the new generation Cherenkov telescope MAGIC.

Microquasars are accreting X-ray binary systems with non-thermal radio jets. In some of these systems the jet is expected to be strongly misaligned with the perpendicular to the orbital plane. If the donor star is an early-type star, the jet could collide with the stellar wind producing a standing shock between the compact object and the stellar surface. Relativistic particles injected by the jet can be re-accelerated and isotropized at the colliding region. If the jet has hadronic content, TeV protons will diffuse into the inner, dense wind leading to gamma-ray and neutrino production from interactions with the matter of the wind. In the case of very powerful jets, the wind pressure can be overbalanced and the jet might impact directly onto the stellar surface. We present estimates of the gamma-ray and neutrino luminosities for different sets of parameters in these scenarios and we briefly discuss the effects of this radiation on the donor star and its detectability with current instruments.

Be/X-ray binaries are systems formed by a massive Be star and a magnetized neutron star, usually in an eccentric orbit. The Be star has strong equatorial winds occasionally forming a circumstellar disk. When the neutron star intersects the disk the accretion rate dramatically increases and a transient accretion disk can be formed around the compact object. This disk can last longer than a single orbit in the case of major outbursts. If the disk rotates faster than the neutron star, the Cheng-Ruderman mechanism can produce a current of relativistic protons that would impact onto the disk surface, producing gamma-rays from neutral pion decays and initiating electromagnetic cascades inside the disk. In this paper we present calculations of the evolution of the disk parameters during both major and minor X-ray events, and we discuss the generation of gamma-ray emission at different energies within a variety of models that include both screened and unscreened disks.

A recent analysis of COMPTEL data has revealed the existence of an unidentified and variable gamma-ray source, designated GRO J1411-64, at MeV energies in the galactic plane (Zhang et al. 2002). We discuss the possibility that such a source can be produced in the accreting neutron star system 2S 1417-624 through the bombardment of the accretion disk by relativistic protons accelerated in the magnetosphere by the Cheng-Ruderman mechanism.