Aims: To investigate the extension of the very-high-energy spectral tail of the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours of good quality data of Crab with the MAGIC telescope, obtained from February 2007 until April 2014. Results: We report the most energetic pulsed emission ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with the ones measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC telescopes the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least $5\times 10^6$. These results strongly suggest IC scattering off low energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder.

We present results from an optical study of two young Balmer-dominated remnants of SNIa in the Large Magellanic Cloud, 0509$-$67.5 and 0519$-$69.0, in an attempt to search for signatures of efficient cosmic-ray (CR) acceleration. We combine proper motion measurements from HST with corresponding optical spectroscopic measurements of the H$\alpha$ line at multiple rim positions from VLT/FORS2 and SALT/RSS and compare our results to published Balmer shock models. Analysis of the optical spectra result in broad H$\alpha$ widths in the range of 1800-4000 km s$^{-1}$ for twelve separate Balmer-dominated filaments that show no evidence for forbidden line emission, the corresponding shock speeds from proper motion measurements span a range of 1700-8500 km s$^{-1}$. Our measured values of shock speeds and broad H$\alpha$ widths in 0509$-$67.5 and 0519$-$69.0 are fit well with a Balmer shock model that does not include effects of efficient CR acceleration. We determine an upper limit of 7%/$\chi$ (95% confidence) on the CR acceleration efficiency for our ensemble of data points, where $\chi$ is the ionization fraction of the pre-shock gas. The upper limits on the individual remnants are 6%/$\chi$ (0509$-$67.5) and 11%/$\chi$ (0519$-$69.0). These upper limits are below the integrated CR acceleration efficiency in the Tycho supernova remnant, where the shocks predominantly show little H$\alpha$ emission, indicating that Balmer-dominated shocks are not efficient CR accelerators.

We have undertaken a comprehensive search for both Lyman Alpha Emitters (LAEs) and Lyman Break Galaxies (LBGs) in the SHARDS Survey of the GOODS-N field. SHARDS is a deep imaging survey, made with the 10.4 m Gran Telescopio Canarias (GTC), employing 25 medium band filters in the range from 500 to 941 nm. This is the first time that both LAEs and LBGs are surveyed simultaneously in a systematic way in a large field. We draw a sample of 1558 sources; 528 of them are LAEs. Most of the sources (1434) show rest-frame UV continua. A minority of them (124) are pure LAEs with virtually no continuum detected in SHARDS. We study these sources from $z\sim3.35$ up to $z\sim6.8$, well into the epoch of reionization. Note that surveys done with just one or two narrow band filters lack the possibility to spot the rest-frame UV continuum present in most of our LAEs. We derive redshifts, Star Formation Rates (SFRs), Ly$\alpha$ Equivalent Widths (EWs) and Luminosity Functions (LFs). Grouping within our sample is also studied, finding 92 pairs or small groups of galaxies at the same redshift separated by less than 60 comoving kpc. In addition, we relate 87 and 55 UV-selected objects with two known overdensities at $z=4.05$ and $z=5.198$, respectively. Finally, we show that surveys made with broad band filters are prone to introduce many unwanted sources ($\sim20$% interlopers), which means that previous studies may be overestimating the calculated LFs, specially at the faint end.

We present a large set of spectral lines detected in the $40''$ central region of the starburst galaxy NGC 253. Observations were obtained with the three instruments SPIRE, PACS and HIFI on board the Herschel Space Observatory, upGREAT on board of the SOFIA airborne observatory, and the ground based APEX telescope. Combining the spectral and photometry products of SPIRE and PACS we model the dust continuum Spectral Energy Distribution (SED) and the most complete $^{12}$CO Line SED reported so far toward the nuclear region of NGC 253. Properties and excitation of the molecular gas were derived from a three-component non-LTE radiative transfer model, using the SPIRE $^{13}$CO lines and ground based observations of the lower-$J$ $^{13}$CO and HCN lines, to constrain the model parameters. Three dust temperatures were identified from the continuum emission, and three components are needed to fit the full CO LSED. Only the third CO component (fitting mostly the HCN and PACS $^{12}$CO lines) is consistent with a shock/mechanical heating scenario. A hot core chemistry is also argued as a plausible scenario to explain the high-$J$ $^{12}$CO lines detected with PACS. The effect of enhanced cosmic ray ionization rates, however, cannot be ruled out, and is expected to play a significant role in the diffuse and dense gas chemistry. This is supported by the detection of ionic species like OH$^+$ and H$_2$O$^+$, as well as the enhanced fluxes of the OH lines with respect to those of H$_2$O lines detected in both PACS and SPIRE spectrum.

Type 2 quasars are an important constituent of active galaxies, possibly representing the evolutionary precursors of traditionally studied type 1 quasars. We characterize the black hole mass ($M_{\rm BH}$) and Eddington ratio ($L_{\rm bol}/L_{\rm Edd}$) for 669 type 2 quasars selected from the Sloan Digital Sky Survey, using black hole masses estimated from the $M_{\rm BH}-\sigma_{\ast}$ relation and bolometric corrections scaled from the extinction-corrected $ [\rm O~{\scriptsize III}]~ \lambda 5007$ luminosity. When stellar velocity dispersions cannot be measured directly from the spectra, we estimate them from the core velocity dispersions of the narrow emission lines $ \rm [O~{\scriptsize II}]~ \lambda\lambda 3726, 3729$, $ \rm [S~{\scriptsize II}]~\lambda\lambda 6716, 6731$, and $[ \rm O~{\scriptsize III}]~ \lambda 5007$, which are shown to trace the gravitational potential of the stars. Energy input from the active nucleus still imparts significant perturbations to the gas kinematics, especially to high-velocity, blueshifted wings. Nonvirial motions in the gas become most noticeable in systems with high Eddington ratios. The black hole masses of our sample of type 2 quasars range from $M_{\rm BH} \approx 10^{6.5} $ to $10^{10.4} \, M_\odot$ (median $10^{8.2} \, M_\odot$). Type 2 quasars have characteristically large Eddington ratios ($L_{\rm bol}/L_{\rm Edd}~\approx 10^{-2.9}-10^{1.8}$; median $10^{-0.7}$), slightly higher than in type~1 quasars of similar redshift; the luminosities of $\sim$20% of the sample formally exceed the Eddington limit. The high Eddington ratios may be consistent with the notion that obscured quasars evolve into unobscured quasars.

Solar and stellar dynamos shed small-scale and large-scale magnetic helicity of opposite signs. However, solar wind observations and simulations have shown that some distance above the dynamo, both the small-scale and large-scale magnetic helicities have reversed signs. With realistic simulations of the solar corona above an active region now being available, we have access to the magnetic field and current density along coronal loops. We show that a sign reversal in the horizontal averages of the magnetic helicity occurs when the local maximum of the plasma beta drops below unity and the field becomes nearly fully force free. Hence, this reversal is expected to occur well within the solar corona and would not directly be accessible to in situ measurements with Parker Solar Probe or SolarOrbiter. We also show that the reversal is associated with subtle changes in the relative dominance of structures with positive and negative magnetic helicity.

The emergence of dipolar magnetic features on the solar surface is an idealization. Most of the magnetic flux emergence occurs in complex multipolar regions. Here we show that the surface pattern of magnetic structures alone can reveal the sign of the underlying magnetic helicity in the nearly force-free coronal regions above. The sign of the magnetic helicity can be predicted to good accuracy by considering the three-dimensional position vectors of three spots on the sphere ordered by their relative strengths at the surface and compute from them the skew product. This product, which is a pseudo scalar, is shown to be a good proxy for the sign of the coronal magnetic helicity.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

Lyman-$\alpha$ (Ly$\alpha$) is a powerful astrophysical probe. Not only is it ubiquitous at high redshifts, it is also a resonant line, making Ly$\alpha$ photons scatter. This scattering process depends on the physical conditions of the gas through which Ly$\alpha$ propagates, and these conditions are imprinted on observables such as the Ly$\alpha$ spectrum and its surface brightness profile. In this work, we focus on a less-used observable capable of probing any scattering process: polarization. We implement the density matrix formalism of polarization into the Monte Carlo radiative transfer code tlac. This allows us to treat it as a quantum mechanical process where single photons develop and lose polarization from scatterings in arbitrary gas geometries. We explore static and expanding ellipsoids, biconical outflows, and clumpy multiphase media. We find that photons become increasingly polarized as they scatter and diffuse into the wings of the line profiles, making scattered Ly$\alpha$ polarized in general. The degree and orientation of Ly$\alpha$ polarization depends on the kinematics and distribution of the scattering HI gas. We find that it generally probes spatial or velocity space asymmetries and aligns itself tangentially to the emission source. We show that the mentioned observables, when studied separately, can leave similar signatures for different source models. We conclude by revealing how a joint analysis of the Ly$\alpha$ spectra, surface brightness profiles, and polarization can break these degeneracies and help us extract unique physical information on galaxies and their environments from their strongest, most prominent emission line.

Models and observations suggest that luminous quasar activity is triggered by mergers, so it should preferentially occur in the most massive primordial dark matter haloes, where the frequency of mergers is expected to be the highest. Since the importance of galaxy mergers increases with redshift, we identify the high-redshift Universe as the ideal laboratory for studying dual AGN. Here we present the X-ray properties of two systems of dual quasars at z=3.0-3.3 selected from the SDSS-DR6 at separations of 6-8 arcsec (43-65kpc) and observed by Chandra for 65ks each. Both members of each pair are detected with good photon statistics to allow us to constrain the column density, spectral slope and intrinsic X-ray luminosity. We also include a recently discovered dual quasar at z=5 (separation of 21 arcsec, 136kpc) for which XMM-Newton archival data allow us to detect the two components separately. Using optical spectra we derived bolometric luminosities, BH masses and Eddington ratios that were compared to those of luminous SDSS quasars in the same redshift ranges. We find that the brighter component of both pairs at z=3.0-3.3 has high luminosities compared to the distribution of SDSS quasars at similar redshift, with J1622A having an order magnitude higher luminosity than the median. This source lies at the luminous end of the z~3.3 quasar luminosity function. While we cannot conclusively state that the unusually high luminosities of our sources are related to their having a close companion, for J1622A there is only a 3% probability that it is by chance.

(Abridged) Astronomers recently started discovering exoplanets around binary systems. Therefore, understanding the formation and evolution of circumbinary disks is crucial for a complete scenario of planet formation. The aim of this paper is to present the detection of a circumbinary disk around Oph-IRS67 and analyse its structure. We present high-angular-resolution (0.4", 60 AU) observations of C17O, H13CO+ , C34S, SO2, C2H and c-C3H2 molecular transitions with ALMA at 0.8 mm. The spectrally and spatially resolved maps reveal the kinematics of the circumbinary disk as well as its chemistry. Molecular abundances are estimated using RADEX. The continuum emission reveals the presence of a circumbinary disk around the two sources. This disk has a diameter of ~620 AU and is well traced by C17O and H13CO+ emission. C2H and c-C3H2 trace a higher-density region which is spatially offset from the sources (~430 AU). Finally, SO2 shows compact emission around one of the sources, Oph-IRS67 B. The molecular transitions which trace the circumbinary disk are consistent with a Keplerian profile on disk scales (< 200 AU) and an infalling profile for envelope scales (> 200 AU). The Keplerian fit leads to a mass of 2.2 Msun. Inferred CO abundances w.r.t. H2 are comparable to the canonical ISM value of 2.7e-4. This study proves the first detection of the circumbinary disk associated with Oph-IRS67. The disk is chemically differentiated from the nearby high-density region. The lack of methanol emission suggests the extended disk dominates the mass budget in the inner- most regions of the protostellar envelope, generating a flat density profile where less material is exposed to high temperatures. Thus, complex organic molecules would be associated with lower column densities. Finally, Oph-IRS67 is a promising candidate for the detection of both circumstellar disks with higher-angular-resolution observations.

We aim at studying the causal link between the knotty jet structure in CARMA 7, a young Class 0 protostar in the Serpens South cluster, and episodic accretion in young protostellar disks. We used numerical hydrodynamics simulations to derive the protostellar accretion history in gravitationally unstable disks around solar-mass protostars. We compared the time spacing between luminosity bursts \Delta\tau_mod, caused by dense clumps spiralling on the protostar, with the differences of dynamical timescales between the knots \Delta\tau_obs in CARMA 7. We found that the time spacing between the bursts have a bi-modal distribution caused by isolated and clustered luminosity bursts. The former are characterized by long quiescent periods between the bursts with \Delta\tau_mod = a few * (10^3-10^4) yr, whereas the latter occur in small groups with time spacing between the bursts \Delta\tau_mod= a few * (10-10^2) yr. For the clustered bursts, the distribution of \Delta\tau_mod in our models can be fit reasonably well to the distribution of \Delta\tau_obs in the protostellar jet of CARMA 7, if a certain correction for the (yet unknown) inclination angle with respect to the line of sight is applied. The K-S test on the model and observational data sets suggests the best-fit values for the inclination angles of 55-80 deg., which become narrower (75-80 deg.) if only strong luminosity bursts are considered. The dynamical timescales of the knots in the jet of CARMA 7 are too short for a meaningful comparison with the long time spacings between isolated bursts in our models. The exact sequences of time spacings between the luminosity bursts in our models and knots in the jet of CARMA 7 were found difficult to match. (abridged)

The innermost regions in active galactic nuclei (AGNs) were not being spatially resolved so far but spectropolarimetry can provide us insight about their hidden physics and the geometry. From spectropolarimetric observations in broad emission lines and assuming equatorial scattering as a dominant polarization mechanism, it is possible to estimate the mass of supermassive black holes (SMBHs). We explore the possibilities and limits and to put constraints on the usage of the method for determining SMBH masses using polarization in broad emission lines by providing more in-depth theoretical modeling. Methods. We use the Monte Carlo radiative transfer code STOKES for exploring polarization of Type 1 AGNs. We model equatorial scattering using flared-disk geometry for a set of different SMBH masses assuming Thomson scattering. In addition to the Keplerian motion in the BLR, we also consider cases of additional radial inflows and vertical outflows. We model the profiles of polarization plane position angle, degree of polarization and total unpolarized line for different BLR geometries and different SMBH masses. Our modeling confirms that the method can be widely used for Type-1 AGNs when viewing inclinations are between 25 and 45 degrees. We show that the distance between the BLR and scattering region (SR) has a significant impact on the mass estimates and the best mass estimates are when the SR is situated at the distance 1.5-2.5 times larger than the outer BLR radius. Our models show that if Keplerian motion can be traced through the polarized line profile, then the direct estimation of the mass of the SMBH can be performed. When radial inflows or vertical outflows are present in the BLR, this method can be applied if velocities of the inflow/outflow are less than 500 km/s. We find that models for NGC4051, NGC4151, 3C273 and PG0844+349 are in good agreements with observations.

The very-high-energy (VHE, $\gtrsim 100$ GeV) $\gamma$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $\gamma$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $\gamma$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL).

The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 {\sigma} using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy.

We study the ellipticity of galaxy cluster halos as characterized by the distribution of cluster galaxies and as measured with weak lensing. We use monte-carlo simulations of elliptical cluster density profiles to estimate and correct for Poisson noise bias, edge bias and projection effects. We apply our methodology to 10,428 SDSS clusters identified by the redMaPPer algorithm with richness above 20. We find a mean ellipticity $= 0.271 \pm 0.002$ (stat) $\pm 0.031$ (sys) corresponding to an axis ratio $= 0.573 \pm 0.002$ (stat) $\pm 0.039$ (sys). We compare this ellipticity of the satellites to the halo shape, through a stacked lensing measurement using optimal estimators of the lensing quadrupole based on Clampitt and Jain (2016). We find a best-fit axis ratio of $0.56 \pm 0.09$ (stat) $\pm 0.03$ (sys), consistent with the ellipticity of the satellite distribution. Thus cluster galaxies trace the shape of the dark matter halo to within our estimated uncertainties. Finally, we restack the satellite and lensing ellipticity measurements along the major axis of the cluster central galaxy's light distribution. From the lensing measurements we infer a misalignment angle with an RMS of ${30^\circ \pm 10}^\circ$ when stacking on the central galaxy. We discuss applications of halo shape measurements to test the effects of the baryonic gas and AGN feedback, as well as dark matter and gravity. The major improvements in signal-to-noise expected with the ongoing Dark Energy Survey and future surveys from LSST, Euclid and WFIRST will make halo shapes a useful probe of these effects.

Circumstantial evidence suggests that magnetism and enhanced X-ray emission are likely correlated in early B-type stars: similar fractions of them ($\sim$ 10 %) are strong and hard X-ray sources and possess strong magnetic fields. It is also known that some B-type stars have spots on their surface. Yet up to now no X-ray activity associated with spots on early-type stars was detected. In this Letter we report the detection of a magnetic field on the B2V star $\rho$ Oph A. Previously, we assessed that the X-ray activity of this star is associated with a surface spot, herewith we establish its magnetic origin. We analyzed FORS2 ESO VLT spectra of $\rho$ Oph A taken at two epochs and detected a longitudinal component of the magnetic field of order of $\sim500$ G in one of the datasets. The detection of the magnetic field only at one epoch can be explained by stellar rotation which is also invoked to explain observed periodic X-ray activity. From archival HARPS ESO VLT high resolution spectra we derived the fundamental stellar parameters of $\rho$ Oph A and further constrained its age. We conclude that $\rho$ Oph A provides strong evidence for the presence of active X-ray emitting regions on young magnetized early type stars.

We describe the hCOSMOS redshift survey of the COSMOS field conducted with the Hectospec spectrograph on the MMT. In the central 1~deg$^2$, the hCOS20.6 subset of the survey is $>90\%$ complete to a limiting $r=20.6$. The hCOSMOS survey includes 1701 new redshifts in the COSMOS field. We also use the total of 4362 new and remeasured objects to derive the age sensitive D$_n4000$ index over the entire redshift interval $0.001\lesssim z\lesssim0.6$. For $85\%$ of the quiescent galaxies in hCOS20.6, we measure the central line-of-sight velocity dispersion. To explore potential uses of this survey, we combine previously measured galaxy sizes, profiles and stellar masses with the spectroscopy. The comparison reveals the known relations among structural, kinematic, and stellar population properties. We also compare redshift and D$_n4000$ distributions of hCOS20.6 galaxies with SHELS; a complete spectroscopic survey of 4~deg$^2$ observed to the same depth. The redshift distributions in the two fields are very different but the D$_n4000$ distribution is remarkably similar. The relation between velocity dispersion and stellar mass for massive hCOS20.6 galaxies is consistent with the local relation from SDSS. Using measured velocity dispersions, we test a photometric proxy calibrated to galaxies in the local universe. The systematic differences between the measured and photometric proxy velocity dispersions are correlated with galaxy dynamical and stellar population properties highlighting the importance of direct spectroscopic measurements.

A former member of the Hyades cluster, $\iota$ Horologii ($\iota$ Hor) is a planet-hosting Sun-like star which displays the shortest coronal activity cycle known to date ($P_{\rm cyc}\sim$ 1.6 yr). With an age of $\sim$$625$ Myr, $\iota$ Hor is also the youngest star with a detected activity cycle. The study of its magnetic properties holds the potential to provide fundamental information to understand the origin of cyclic activity and stellar magnetism in late-type stars. In this series of articles, we present the results of a comprehensive project aimed at studying the evolving magnetic field in this star and how this evolution influences its circumstellar environment. This paper summarizes the first stage of this investigation, with results from a long-term observing campaign of $\iota$ Hor using ground-based high-resolution spectropolarimetry. The analysis includes precise measurements of the magnetic activity and radial velocity of the star, and their multiple time-scales of variability. In combination with values reported in the literature, we show that the long-term chromospheric activity evolution of $\iota$ Hor follows a beating pattern, caused by the superposition of two periodic signals of similar amplitude at $P_{1} \simeq 1.97 \pm 0.02$ yr and $P_{2} \simeq 1.41 \pm 0.01$ yr. Additionally, using the most recent parameters for $\iota$ Hor b in combination with our activity and radial velocity measurements, we find that stellar activity dominates the radial velocity residuals, making the detection of additional planets in this system challenging. Finally, we report here the first measurements of the surface longitudinal magnetic field strength of $\iota$ Hor, which displays varying amplitudes within $\pm4$ G and served to estimate the rotation period of the star ($P_{\rm rot} =~7.70^{+0.18}_{-0.67}$ d).

Many young, massive stars are found in close binaries. Using population synthesis simulations we predict the likelihood of a companion star being present when these massive stars end their lives as core-collapse supernovae (SNe). We focus on stripped-envelope SNe, whose progenitors have lost their outer hydrogen and possibly helium layers before explosion. We use these results to interpret new Hubble Space Telescope observations of the site of the broad-lined Type Ic SN 2002ap, 14 years post-explosion. For a subsolar metallicity consistent with SN 2002ap, we expect a main-sequence companion present in about two thirds of all stripped-envelope SNe and a compact companion (likely a stripped helium star or a white dwarf/neutron star/black hole) in about 5% of cases. About a quarter of progenitors are single at explosion (originating from initially single stars, mergers or disrupted systems). All the latter scenarios require a massive progenitor, inconsistent with earlier studies of SN 2002ap. Our new, deeper upper limits exclude the presence of a main-sequence companion star $>8$-$10$ Msun, ruling out about 40% of all stripped-envelope SN channels. The most likely scenario for SN 2002ap includes nonconservative binary interaction of a primary star initially $\lesssim 23$ Msun. Although unlikely ($<$1% of the scenarios), we also discuss the possibility of an exotic reverse merger channel for broad-lined Type Ic events. Finally, we explore how our results depend on the metallicity and the model assumptions and discuss how additional searches for companions can constrain the physics that governs the evolution of SN progenitors.

To better understand how planets form, it is important to study planet occurrence rates as a function of stellar mass. However, estimating masses of field stars is often difficult. Over the past decade, a controversy has arisen about the inferred occurrence rate of gas-giant planets around evolved intermediate-mass stars -- the so-called `retired A-stars'. The high masses of these red-giant planet hosts, derived using spectroscopic information and stellar evolution models, have been called into question. Here we address the controversy by determining the masses of eight evolved planet-hosting stars using asteroseismology. We compare the masses with spectroscopic-based masses from the Exoplanet Orbit Database that were previously adopted to infer properties of the exoplanets and their hosts. We find a significant one-sided offset between the two sets of masses for stars with spectroscopic masses above roughly 1.6Msun, suggestive of an average 15--20% overestimate of the adopted spectroscopic-based masses. The only star in our sample well below this mass limit is also the only one not showing this offset. Finally, we note that the scatter across literature values of spectroscopic-based masses often exceed their formal uncertainties, making it comparable to the offset we report here.

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black-hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV) regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behavior in a multiwavelength context, we compare our results with \textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI}, \textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index $\Gamma=3.2$. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95\% confidence level for energies above 200 GeV at $2.6\times10^{-12}$~photons cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$, respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of the MAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

We investigate the evolution of mass segregation in initially sub-structured young embedded star clusters with two different background potentials mimicking the gas. Our clusters are initially in virial or sub-virial global states and have different initial distributions for the most massive stars: randomly placed, initially mass segregated or even inverse segregation. By means of N-body simulation we follow their evolution for 5 Myr. We measure the mass segregation using the minimum spanning tree method Lambda_MSR and an equivalent restricted method. Despite this variety of different initial conditions, we find that our stellar distributions almost always settle very fast into a mass segregated and more spherical configuration, suggesting that once we see a spherical or nearly spherical embedded star cluster, we can be sure it is mass segregated no matter what the real initial conditions were. We, furthermore, report under which circumstances this process can be more rapid or delayed, respectively.