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 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 study the kinematics and the stellar populations of the bulge and disc of the spiral galaxy NGC 3521. At each position in the field of view, we separate the contributions of the bulge and the disc from the total observed spectrum and study their kinematics, age, and metallicities independently. Their properties are clearly distinct: the bulge rotates more slowly, has a higher velocity dispersion, and is less luminous than the disc. We identify three main populations of stars in NGC 3521: old ($\geq7$ Gyr), intermediate ($\approx$ 3 Gyr), and young ($\leq$1 Gyr). The mass and light of NGC 3521 are dominated by the intermediate stellar population. The youngest population contributes mostly to the disc component and its contribution increases with radius. We also study the luminosity-weighed properties of the stars in NGC 3521. Along the photometric major axis, we find: i) no age gradient for the stars in the bulge, and a negative age gradient for the stars in the disc; ii) negative metallicity gradients and sub-solar $\alpha$-enhancement for both the bulge and the disc. We propose the following picture for the formation of NGC 3521: initial formation a long time ago ($\geq 7$ Gyr), followed by a second burst of star formation or a merger ($\approx$ 3 Gyrs ago), which contributed predominantly to the mass-build up of the bulge. Recently ($\leq 1$ Gyr), the disc of NGC 3521 experienced an additional episode of star formation that started in the innermost regions.

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.

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 present our recent results about the extended narrow-line region (ENLR) of two nearby Seyfert 2 galaxies (IC 5063 and NGC 7212) obtained by modelling the observed line profiles and spectra with composite models (photoionization+shocks) in the different regions surrounding the AGN. Then, we compare the Seyfert 2 ENLRs with the very extended one recently discovered in the narrow-line Seyfert 1 (NLS1) galaxy Mrk 783. We have found several evidences of interaction between the ISM of the galaxies and their radio jets, such as a) the contribution of shocks in ionizing the high velocity gas, b) the complex kinematics showed by the profile of the emission lines, c) the high fragmentation of matter, etc. The results suggest that the ENLR of IC 5063 have a hollow bi-conical shape, with one edge aligned to the galaxy disk, which may cause some kind of dependence on velocity of the ionization parameter. Regarding the Mrk 783 properties, it is found that the extension of the optical emission is almost twice the size of the radio one and it seems due to the AGN activity, although there is contamination by star formation around 12 arcsec from the nucleus. Diagnostic diagrams excluded the contribution of star formation in IC 5063 and NGC 7212, while the shock contribution was used to explain the spectra emitted by their high velocity gas.

Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the MAGIC collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent~CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of $E_{\mathrm{QG}_1} > 5.5\cdot 10^{17}$GeV ($4.5\cdot 10^{17}$GeV) for a linear, and $E_{\mathrm{QG}_2} > 5.9\cdot 10^{10}$GeV ($5.3\cdot 10^{10}$GeV) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, found to worsen the statistical limits by about 36--42\%. Our constraints would have resulted much more competitive if the intrinsic pulse shape of the pulsar between 200GeV and 400GeV was understood in sufficient detail and allowed inclusion of events well below 400GeV.

MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we develop a dedicated Moon-adapted analysis to characterize the performance of MAGIC under moonlight. We evaluate energy threshold, angular resolution and sensitivity of MAGIC under different background light levels, based on Crab Nebula observations and tuned Monte Carlo simulations. This study includes observations taken under non-standard hardware configurations, such as reducing the camera photomultiplier tubes gain by a factor ~1.7 (Reduced HV settings) with respect to standard settings (Nominal HV) or using UV-pass filters to strongly reduce the amount of moonlight reaching the cameras of the telescopes. The Crab Nebula spectrum is correctly reconstructed in all the studied illumination levels, that reach up to 30 times brighter than under dark conditions. The main effect of moonlight is an increase in the analysis energy threshold and in the systematic uncertainties on the flux normalization. The sensitivity degradation is constrained to be below 10%, within 15-30% and between 60 and 80% for Nominal HV, Reduced HV and UV-pass filter observations, respectively. No worsening of the angular resolution was found. Thanks to observations during moonlight, the maximal duty cycle of MAGIC can be increased from ~18%, under dark nights only, to up to ~40% in total with only moderate performance degradation.

We present new stringent limits on the mass $M_{bh}$ of the central supermassive black hole for a sample of 7 nearby galaxies. Our $M_{bh}$ estimates are based on the dynamical modeling of the central width of the nebular emission lines measured over subarcsecond apertures with the Hubble Space Telescope. The central stellar velocity dispersion $\sigma_c$ of the sample galaxies is derived from new long-slit spectra from ground-based observations and the bulge effective radius is obtained from a two-dimensional photometric decomposition of the i-band images from the Sloan Digital Sky Survey. The derived stringent $M_{bh}$ limits run parallel and above the $M_{bh}-\sigma_c$ relation with no systematic trend depending on the galaxy distance or morphology. This gives further support to previous findings suggesting that the nuclear gravitational potential is remarkably well traced by the width of the nebular lines when the gas is centrally peaked. With our investigation, the number of galaxies with stringent $M_{bh}$ limits obtained from nebular-line width increases to 114 and can be used for studying the scaling relations between $M_{bh}$ and properties of their host galaxies.

B1957+20 is a millisecond pulsar located in a black widow type compact binary system with a low mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high energy {\gamma} rays in the inverse Compton process. We performed extensive observations with the MAGIC telescopes of B1957+20. We interpret results in the framework of a few different models, namely emission from the vicinity of the millisecond pulsar, the interaction of the pulsar and stellar companion wind region, or bow shock nebula. No significant steady very high energy {\gamma}-ray emission was found. We derived a 95% confidence level upper limit of 3.0 x 10 -12 cm -2 s -1 on the average {\gamma}-ray emission from the binary system above 200 GeV. The upper limits obtained with MAGIC constrain, for the first time, different models of the high-energy emission in B1957+20. In particular, in the inner mixed wind nebula model with mono-energetic injection of electrons, the acceleration efficiency of electrons is constrained to be below ~(2-10)% of the pulsar spin down power. For the pulsar emission, the obtained upper limits for each emission peak are well above the exponential cut-off fits to the Fermi-LAT data, extrapolated to energies above 50 GeV. The MAGIC upper limits can rule out a simple power-law tail extension through the sub-TeV energy range for the main peak seen at radio frequencies.

The extragalactic VHE gamma-ray sky is rich in blazars. These are jetted active galactic nuclei viewed at a small angle to the line-of-sight. Only a handful of objects viewed at a larger angle are known so far to emit above 100 GeV. Multi-wavelength studies of such objects up to the highest energies provide new insights into the particle and radiation processes of active galactic nuclei. We report the results from the first multi-wavelength campaign observing the TeV detected nucleus of the active galaxy IC 310, whose jet is observed at a moderate viewing angle of 10 deg - 20 deg. The multi-instrument campaign was conducted between 2012 Nov. and 2013 Jan., and involved observations with MAGIC, Fermi, INTEGRAL, Swift, OVRO, MOJAVE and EVN. These observations were complemented with archival data from the AllWISE and 2MASS catalogs. A one-zone synchrotron self-Compton model was applied to describe the broad-band spectral energy distribution. IC 310 showed an extraordinary TeV flare at the beginning of the campaign, followed by a low, but still detectable TeV flux. Compared to previous measurements, the spectral shape was found to be steeper during the low emission state. Simultaneous observations in the soft X-ray band showed an enhanced energy flux state and a harder-when-brighter spectral shape behaviour. No strong correlated flux variability was found in other frequency regimes. The broad-band spectral energy distribution obtained from these observations supports the hypothesis of a double-hump structure. The harder-when-brighter trend in the X-ray and VHE emission is consistent with the behaviour expected from a synchrotron self-Compton scenario. The contemporaneous broad-band spectral energy distribution is well described with a one-zone synchrotron self-Compton model using parameters that are comparable to those found for other gamma-ray-emitting misaligned blazars.

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

In this work we present data from observations with the MAGIC telescopes of SN 2014J detected in January 21 2014, the closest Type Ia supernova since Imaging Air Cherenkov Telescopes started to operate. We probe the possibility of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early stages of Type Ia supernova explosions. We performed follow-up observations after this supernova explosion for 5 days, between January 27 and February 2 in 2014. We search for gamma-ray signal in the energy range between 100 GeV and several TeV from the location of SN 2014J using data from a total of $\sim5.5$ hours of observations. Prospects for observing gamma-rays of hadronic origin from SN 2014J in the near future are also being addressed. No significant excess was detected from the direction of SN 2014J. Upper limits at 95$\%$ confidence level on the integral flux, assuming a power-law spectrum, d$F/$d$E\propto E^{-\Gamma}$, with a spectral index of $\Gamma=2.6$, for energies higher than 300 GeV and 700 GeV, are established at $1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$, respectively. For the first time, upper limits on the VHE emission of a Type Ia supernova are established. The energy fraction isotropically emitted into TeV gamma rays during the first $\sim10$ days after the supernova explosion for energies greater than 300 GeV is limited to $10^{-6}$ of the total available energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario, the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be detectable in the future by any current or planned generation of Imaging Atmospheric Cherenkov Telescopes.

New photometric and long-slit spectroscopic observations are presented for NGC 7113, PGC 1852, and PGC 67207 which are three bright galaxies residing in low-density environments. The surface-brightness distribution is analysed from the K_S-band images taken with adaptive optics at the Gemini North Telescope and the ugriz-band images from the Sloan Digital Sky Survey while the line-of-sight stellar velocity distribution and line-strength Lick indices inside the effective radius are measured along several position angles. The age, metallicity, and alpha-element abundance of the galaxies are estimated from single stellar-population models. In spite of the available morphological classification, images show that PGC 1852 is a barred spiral which we do not further consider for mass modelling. The structural parameters of the two early-type galaxies NGC 7113 and PGC 67207 are obtained from a two-dimensional photometric decomposition and the mass-to-light ratio of all the (luminous and dark) mass that follows the light is derived from orbit-based axisymmetric dynamical modelling together with the mass density of the dark matter halo. The dynamically derived mass that follows the light is about a factor of 2 larger than the stellar mass derived using stellar-population models with Kroupa initial mass function. Both galaxies have a lower content of halo dark matter with respect to early-type galaxies in high-density environments and in agreement with the predictions of semi-analytical models of galaxy formation.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components, spatially indistinguishable by gamma-ray instruments, but separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energy distribution of QSO B0218+357 can give information on the energetics of z ~ 1 very high energy gamma- ray sources. Moreover the gamma-ray emission can also be used as a probe of the extragalactic background light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during the expected arrival time of the delayed component of the emission. The MAGIC and Fermi-LAT observations were accompanied by quasi-simultaneous optical data from the KVA telescope and X-ray observations by Swift-XRT. We construct a multiwavelength spectral energy distribution of QSO B0218+357 and use it to model the source. The GeV and sub-TeV data, obtained by Fermi-LAT and MAGIC, are used to set constraints on the extragalactic background light. Results. Very high energy gamma-ray emission was detected from the direction of QSO B0218+357 by the MAGIC telescopes during the expected time of arrival of the trailing component of the flare, making it the farthest very high energy gamma-ray sources detected to date. The observed emission spans the energy range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy distribution of QSO B0218+357 is consistent with current extragalactic background light models. The broad band emission can be modeled in the framework of a two zone external Compton scenario, where the GeV emission comes from an emission region in the jet, located outside the broad line region.

We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Mets\"ahovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.

The gamma-ray binary LS I +61$^{\circ}$303 is a well established source from centimeter radio up to very high energy (VHE; E$>$100 GeV). Its broadband emission shows a periodicity of $\sim$26.5 days, coincident with the orbital period. A longer (super-orbital) period of 1667 $\pm$ 8 days was discovered in radio and confirmed in optical and high energy (HE; E>100 MeV) gamma-ray observations. We present a four-year campaign performed by MAGIC together with archival data concentrating on a search for a long timescale signature in the VHE emission. We focus on the search for super-orbital modulation of the VHE peak and on the search for correlations between TeV emission and optical determination of the extension of the circumstellar disk. A four-year campaign has been carried out by MAGIC. The source was observed during the orbital phases when the periodic VHE outbursts have occurred ($\phi$=0.55-0.75). Additionally, we included archival MAGIC observations and data published by the VERITAS collaboration in these studies. For the correlation studies, LS I +61$^{\circ}$303 has also been observed during the orbital phases where sporadic VHE emission had been detected in the past ($\phi$=0.75-1.0). These MAGIC observations were simultaneous with optical spectroscopy from the LIVERPOOL telescope. The TeV flux of the periodical outburst in orbital phases $\phi$=0.5--0.75 was found to show yearly variability consistent with the $\sim$4.5 years long-term modulation found in the radio band. This modulation of the TeV flux can be well described by a sine function with the best fit period of $1610\pm 58$ days. The complete dataset span two super-orbital periods. There is no evidence for a correlation between the TeV emission and the mass-loss rate of the Be star but this may be affected by the strong, short timescale (as short as intra-day) variation displayed by the H$\alpha$ fluxes.

Small kinematically-decoupled stellar discs with scalelengths of a few tens of parsec are known to reside in the centre of galaxies. Different mechanisms have been proposed to explain how they form, including gas dissipation and merging of globular clusters. Using archival Hubble Space Telescope imaging and ground-based integral-field spectroscopy, we investigated the structure and stellar populations of the nuclear stellar disc hosted in the interacting SB0 galaxy NGC 1023. The stars of the nuclear disc are remarkably younger and more metal rich with respect to the host bulge. These findings support a scenario in which the nuclear disc is the end result of star formation in metal enriched gas piled up in the galaxy centre. The gas can be of either internal or external origin, i.e. from either the main disc of NGC 1023 or the nearby satellite galaxy NGC 1023A. The dissipationless formation of the nuclear disc from already formed stars, through the migration and accretion of star clusters into the galactic centre is rejected.

The old and metal rich open cluster NGC 6253 was observed with the FLAMES multi-object spectrograph during an extensive radial velocity campaign monitoring 317 stars with a median of 15 epochs per object. All the targeted stars are located along the upper main sequence of the cluster between 14.8 $<$ V $<$ 16.5. Fifty nine stars are confirmed cluster members both by radial velocities and proper motions and do not show evidence of variability. We detected 45 variable stars among which 25 belong to NGC 6253. We were able to derive an orbital solution for 4 cluster members (and for 2 field stars) yielding minimum masses in between $\sim$90 M$\rm_J$ and $\sim$460 M$\rm_J$ and periods between 3 and 220 days. Simulations demonstrated that this survey was sensitive to objects down to 30 M$\rm_J$ at 10 days orbital periods with a detection efficiency equal to 50%. On the basis of these results we concluded that the observed frequency of binaries down to the hydrogen burning limit and up to 20 days orbital period is around (1.5$\pm$1.3)% in NGC 6253. The overall observed frequency of binaries around the sample of cluster stars is (13$\pm$3)%. The median radial velocity precision achieved by the GIRAFFE spectrograph in this magnitude range was around $\sim$240m$\rm\,s^{-1}$ ($\sim$180 m$\rm\,s^{-1}$ for UVES). Based on a limited follow-up analysis of 7 stars in our sample with the HARPS spectrograph we determined that a precision of 35 m $\rm s^{-1}$ can be reached in this magnitude range, offering the possibility to further extend the variability analysis into the substellar domain. Prospects are even more favourable once considering the upcoming ESPRESSO spectrograph at VLT.

Context. In the last five years the Fermi Large Area Telescope (LAT) instrument detected GeV {\gamma}-ray emission from five novae. The GeV emission can be interpreted in terms of an inverse Compton process of electrons accelerated in a shock. In this case it is expected that protons in the same conditions can be accelerated to much higher energies. Consequently they may produce a second component in the {\gamma}-ray spectrum at TeV energies. Aims. We aim to explore the very-high-energy domain to search for {\gamma}-ray emission above 50 GeV and to shed light on the acceleration process of leptons and hadrons in nova explosions. Methods. We have performed observations with the MAGIC telescopes of the classical nova V339 Del shortly after the 2013 outburst, triggered by optical and subsequent GeV {\gamma}-ray detec- tions. We also briefly report on VHE observations of the symbiotic nova YY Her and the dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the analysis of con- temporaneous Fermi-LAT data of the sources. The TeV and GeV observations are compared in order to evaluate the acceleration parameters for leptons and hadrons. Results. No significant TeV emission was found from the studied sources. We computed upper limits on the spectrum and night-by-night flux. The combined GeV and TeV observations of V339 Del limit the ratio of proton to electron luminosities to Lp<~0.15 Le.

We disentangle two counter-rotating stellar components in NGC 4191 and characterize their physical properties (kinematics, morphology, age, metallicity, and abundance ratio). We performed a spectroscopic decomposition on integral field data to separate the contribution of two stellar components to the observed galaxy spectrum across the field of view. We also performed a photometric decomposition, modelling the galaxy with a S\'ersic bulge and two exponential disks of different scale length, with the aim of associating these structural components with the kinematic components. We measured the equivalent width of the absorption line indices on the best fit that represent the kinematic components and compared our measurements to the predictions of stellar population models. We have evidence that the line-of-sight velocity distributions (LOSVDs) are consistent with the presence of two distinct kinematic components. The combined information of the intensity of the LOSVDs and photometry allows us to associate the S\'ersic bulge and the outer disk with the main kinematic component, and the inner disk with the secondary kinematic component. The two kinematic stellar components counter-rotate with respect to each other. The main component is the most luminous and massive, and it rotates slower than the secondary component, which rotates along the same direction as the ionized gas. We also found that the two kinematic components have the same solar metallicity and sub-solar abundance ratio, without the presence of significant radial gradients. On the other hand, their ages show strong negative gradients and the possible indication that the secondary component is the youngest. We interpret our results in light of recent cosmological simulations and suggest gas accretion along two filaments as the formation mechanism of the stellar counter-rotating components in NGC 4191 (Abridged).

We investigate the stellar population properties of a sample of 24 massive quenched galaxies at $1.25<z_\mathrm{spec}<2.09$ identified in the COSMOS field with our Subaru/MOIRCS near-IR spectroscopic observations. Tracing the stellar population properties as close to their major formation epoch as possible, we try to put constraints on the star formation history, post-quenching evolution, and possible progenitor star-forming populations for such massive quenched galaxies. By using a set of Lick absorption line indices on a rest-frame optical composite spectrum, the average age, metallicity [Z/H], and $\alpha$-to-iron element abundance ratio [$\alpha$/Fe] are derived as $\log(\mathrm{age}/\mathrm{Gyr})=0.04_{-0.08}^{+0.10}$, $\mathrm{[Z/H]}=0.24_{-0.14}^{+0.20}$, and $[\alpha/\mathrm{Fe}]=0.31_{-0.12}^{+0.12}$, respectively. If our sample of quenched galaxies at $\langle z \rangle = 1.6$ is evolved passively to $z=0$, their stellar population properties will align in excellent agreement with local counterparts at similar stellar velocity dispersions, which qualifies them as progenitors of local massive early-type galaxies. Redshift evolution of stellar population ages in quenched galaxies combined with low redshift measurements from the literature suggests a formation redshift of $z_\mathrm{f} \sim 2.3$ around which the bulk of stars in these galaxies have been formed. The measured [$\alpha$/Fe] value indicates a star formation timescale of $\lesssim 1$ Gyr, which can be translated into a specific star formation rate of $\simeq 1\,\mathrm{Gyr}^{-1}$ prior to quenching. Based on these findings, we discuss identifying possible progenitor star-forming galaxies at $z \simeq 2.3$. We identify normal star-forming galaxies, i.e, those on the star-forming main sequence, followed by a rapid quenching event, as likely precursors of the quenched galaxies at $\langle z \rangle = 1.6$ presented here.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

The high-frequency-peaked BL Lac (HBL) 1ES 0806+524 (z = 0.138) was discovered in VHE $\gamma$ rays in 2008. Until now, the broad-band spectrum of 1ES 0806+524 has been only poorly characterized, in particular at high energies. We analysed multiwavelength observations from $\gamma$ rays to radio performed from 2011 January to March, which were triggered by the high activity detected at optical frequencies. These observations constitute the most precise determination of the broad-band emission of 1ES 0806+524 to date. The stereoscopic MAGIC observations yielded a $\gamma$-ray signal above 250 GeV of $(3.7 \pm 0.7)$ per cent of the Crab Nebula flux with a statistical significance of 9.9 $\sigma$. The multiwavelength observations showed significant variability in essentially all energy bands, including a VHE $\gamma$-ray flare that lasted less than one night, which provided unprecedented evidence for short-term variability in 1ES 0806+524. The spectrum of this flare is well described by a power law with a photon index of $2.97 \pm 0.29$ between $\sim$150 GeV and 1 TeV and an integral flux of $(9.3 \pm 1.9)$ per cent of the Crab Nebula flux above 250 GeV. The spectrum during the non-flaring VHE activity is compatible with the only available VHE observation performed in 2008 with VERITAS when the source was in a low optical state. The broad-band spectral energy distribution can be described with a one-zone Synchrotron Self Compton model with parameters typical for HBLs, indicating that 1ES 0806+524 is not substantially different from the HBLs previously detected.