Far-infrared (FIR) dust polarimetry enables the study of interstellar magnetic fields via tracing of the polarized emission from dust grains that are partially aligned with the direction of the field. The advent of high quality polarimetric data has permitted the use of statistical methods to extract both the direction and magnitude of the magnetic field. This work presents an effort to produce maps (integrated along the line of sight) of the magnetic field strength in both the line-of-sight (LOS) direction and in the plane-of-sky (POS) in the Orion Molecular Cloud (OMC-1) using all four HAWC+ polarimetry bands (53, 89, 154 and 214 $\micron$). Using the statistics of the polarization data (angles and fractions) along with previously-reported Zeeman measurements, information about the three-dimensional magnetic field configuration was inferred over the OMC-1 region. POS magnetic field strengths of up to 2 mG were determined in OMC-1. The strongest fields were observed near the BN/KL object, while the OMC-1 bar shows strengths of up to few hundreds of $\mu$G. The LOS magnetic field has been determined to have comparable strengths to the POS components but it is differently distributed across the cloud. These estimates of the magnetic field components were used to produce more reliable maps of the mass-to-magnetic flux ratio ($M/\Phi$) -- a proxy for probing the conditions for star formation in molecular clouds -- and determine regions of sub- and super-criticality in OMC-1. Such maps can provide invaluable input and comparison to MHD simulations of star formation processes in filamentary structures of molecular clouds.

Exoplanets orbiting very close to their host star are strongly irradiated. This can lead the upper atmospheric layers to expand and evaporate into space. The metastable helium (HeI) triplet at 1083.3nm has recently been shown to be a powerful diagnostic to probe extended and escaping exoplanetary atmosphere. We perform high-resolution transmission spectroscopy of the transiting hot Jupiter HD189733b with the GIARPS (GIANO-B + HARPS-N) observing mode of the Telescopio Nazionale Galileo, taking advantage of the simultaneous optical+near infrared spectral coverage to detect HeI in the planet's extended atmosphere and to gauge the impact of stellar magnetic activity on the planetary absorption signal. Observations were performed during five transit events of HD189733b. By comparison of the in- and out-of-transit GIANO-B observations we compute high-resolution transmission spectra, on which we perform equivalent width measurements and light-curves analyses to gauge the excess in-transit absorption in the HeI triplet. We detect an absorption signal during all five transits. The mean in-transit absorption depth amounts to 0.75+/-0.03%. We detect night-to-night variations in the HeI absorption signal likely due to the transit events occurring in presence of stellar surface inhomogeneities. We evaluate the impact of stellar-activity pseudo-signals on the true planetary absorption using a comparative analysis of the HeI and the H$\alpha$ lines. We interpret the time-series of the HeI absorption lines in the three nights not affected by stellar contamination -exhibiting a mean in-transit absorption depth of 0.77+/-0.04%- using a 3-d atmospheric code. Our simulations suggest that the helium layers only fill part of the Roche lobe. Observations can be explained with a thermosphere heated to $\sim$12000 K, expanding up to $\sim$1.2 planetary radii, and losing $\sim$1 g/s of metastable helium.

The existence of hot Jupiters is still not well understood. Two main channels are thought to be responsible for their current location: a smooth planet migration through the proto-planetary disk or the circularization of an initial high eccentric orbit by tidal dissipation leading to a strong decrease of the semimajor axis. Different formation scenarios result in different observable effects, such as orbital parameters (obliquity/eccentricity), or frequency of planets at different stellar ages. In the context of the GAPS Young-Objects project, we are carrying out a radial velocity survey with the aim to search and characterize young hot-Jupiter planets. Our purpose is to put constraints on evolutionary models and establish statistical properties, such as the frequency of these planets from a homogeneous sample. Since young stars are in general magnetically very active, we performed multi-band (visible and near-infrared) spectroscopy with simultaneous GIANO-B + HARPS-N (GIARPS) observing mode at TNG. This helps to deal with stellar activity and distinguish the nature of radial velocity variations: stellar activity will introduce a wavelength-dependent radial velocity amplitude, whereas a Keplerian signal is achromatic. As a pilot study, we present here the cases of two already claimed hot Jupiters orbiting young stars: HD285507 b and AD Leo b. Our analysis of simultaneous high-precision GIARPS spectroscopic data confirms the Keplerian nature of HD285507's radial velocities variation and refines the orbital parameters of the hot Jupiter, obtaining an eccentricity consistent with a circular orbit. On the other hand, our analysis does not confirm the signal previously attributed to a planet orbiting AD Leo. This demonstrates the power of the multi-band spectroscopic technique when observing active stars.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

We investigate the relation between characteristics of coronal mass ejections and parameterizations of the eruptive capability of solar active regions widely used in solar flare prediction schemes. These parameters, some of which are explored for the first time, are properties related to topological features, namely, magnetic polarity inversion lines (MPILs) that indicate large amounts of stored non-potential (i.e. free) magnetic energy. We utilize the Space Weather Database of Notifications, Knowledge, Information (DONKI) and the Large Angle and Spectrometric Coronograph (LASCO) databases to find flare-associated coronal mass ejections and their kinematic characteristics while properties of MPILs are extracted from Helioseismic and Magnetic Imager (HMI) vector magnetic-field observations of active regions to extract the properties of source-region MPILs. The correlation between all properties and the characteristics of CMEs ranges from moderate to very strong. More significant correlations hold particularly for fast CMEs, which are most important in terms of adverse space-weather manifestations. Non-neutralized currents and the length of the main MPIL exhibit significantly stronger correlations than the rest of the properties. This finding supports a causal relationship between coronal mass ejections and non-neutralized electric currents in highly sheared, conspicuous MPILs. In addition, non-neutralized currents and MPIL length carry distinct, independent information as to the eruptive potential of active regions. The combined total amount of non-neutralized electric currents and the length of the main polarity inversion line, therefore, reflect more efficiently than other parameters the eruptive capacity of solar active regions and the CME kinematic characteristics stemming from these regions.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

(Abridged) We perform a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D$_{\rm 1}$ D$_{\rm 2}$, and He I D$_{\rm 3}$) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. The value of the solar rotation period is found in all the activity indicators, with the only exception being H$\delta$. The derived values vary from 26.29 days (H$\gamma$ line) to 31.23 days (He I). From an analysis of sliding periodograms we find that in most of the activity indicators the spectral power is split into several "bands" of periods around 26 and 30 days, that might be explained by the migration of active regions between the equator and a latitude of $\sim$ 30$^{\circ}$, spot evolution or a combination of both effects. In agreement with previous works a typical lifetime of active regions of $\sim$ ten rotation periods is inferred from the pooled variance diagrams. We find that H$\alpha$, H$\beta$, H$\gamma$, H$\epsilon$, and He I show a significant correlation with the S index. Significant correlations between the contrast, bisector span, and the heliocentric radial velocity with the activity indexes are also found. We show that the full width at half maximum, the bisector, and the disc-integrated magnetic field correlate with the radial velocity variations. The correlation of the S index and H$\alpha$ changes with time, increasing with larger sun spot numbers and solar irradiance. A similar tendency with the S index - radial velocity correlation is also present in the data. Our results are consistent with a scenario in which higher activity favours the correlation between the S index and the H$\alpha$ activity indicators and between the S index and radial velocity variations.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

GIARPS (GIAno \& haRPS) is a project devoted to have on the same focal station of the Telescopio Nazionale Galileo (TNG) both high resolution spectrographs, HARPS-N (VIS) and GIANO-B (NIR), working simultaneously. This could be considered the first and unique worldwide instrument providing cross-dispersed echelle spectroscopy at a resolution of 50,000 in the NIR range and 115,000 in the VIS and over in a wide spectral range ($0.383 - 2.45\ \mu$m) in a single exposure. The science case is very broad, given the versatility of such an instrument and its large wavelength range. A number of outstanding science cases encompassing mainly extra-solar planet science starting from rocky planets search and hot Jupiters to atmosphere characterization can be considered. Furthermore both instruments can measure high precision radial velocities by means the simultaneous thorium technique (HARPS-N) and absorbing cell technique (GIANO-B) in a single exposure. Other science cases are also possible. GIARPS, as a brand new observing mode of the TNG started after the moving of GIANO-A (fiber fed spectrograph) from Nasmyth-A to Nasmyth-B where it was re-born as GIANO-B (no more fiber feed spectrograph). The official Commissioning finished on March 2017 and then it was offered to the community. Despite the work is not finished yet. In this paper we describe the preliminary scientific results obtained with GIANO-B and GIARPS observing mode with data taken during commissioning and first open time observations.

Solar active regions (ARs) that produce major flares typically exhibit strong plasma shear flows around photospheric magnetic polarity inversion lines (MPILs). It is therefore important to quantitatively measure such photospheric shear flows in ARs for a better understanding of their relation to flare occurrence. Photospheric flow fields were determined by applying the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) method to a large data set of 2,548 co-aligned pairs of AR vector magnetograms with 12-min separation over the period 2012-2016. From each AR flow-field map, three shear-flow parameters were derived corresponding to the mean (<S>), maximum (S_max) and integral (S_sum) shear-flow speeds along strong-gradient, strong-field MPIL segments. We calculated flaring rates within 24 hr as a function of each shear-flow parameter, and also investigated the relation between the parameters and the waiting time ({\tau}) until the next major flare (class M1.0 or above) after the parameter observation. In general, it is found that the larger S_sum an AR has, the more likely it is for the AR to produce flares within 24 hr. It is also found that among ARs which produce major flares, if one has a larger value of S_sum then {\tau} generally gets shorter. These results suggest that large ARs with widespread and/or strong shear flows along MPILs tend to not only be more flare productive, but also produce major flares within 24 hr or less.

Efficient prediction of solar flares relies on parameters that quantify the eruptive capability of solar active regions. Several such quantitative predictors have been proposed in the literature, inferred mostly from photospheric magnetograms and/or white-light observations. Two of them are the Ising energy and the sum of the total horizontal magnetic field gradient. The former has been developed from line-of-sight magnetograms, while the latter uses sunspot detections and characteristics, based on continuum images. Aiming to include these parameters in an automated prediction scheme, we test their applicability on regular photospheric magnetic field observations provided by the Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). We test their efficiency as predictors of flaring activity on a representative sample of active regions and investigate possible modifications of these quantities. The Ising energy appears to be an efficient predictor, and the efficiency is even improved if it is modified to describe interacting magnetic partitions or sunspot umbrae. The sum of the horizontal magnetic field gradient appears to be slightly more promising than the three variations of the Ising energy we implement in this article. The new predictors are also compared with two very promising predictors: the effective connected magnetic field strength and the total unsigned non-neutralized current. Our analysis shows that the efficiency of morphological predictors depends on projection effects in a nontrivial way. All four new predictors are found useful for inclusion in an automated flare forecasting facility, such as the Flare Likelihood and Region Eruption Forecasting (FLARECAST), but their utility, among others, will ultimately be determined by the validation effort underway in the framework of the FLARECAST project.

The SPHERE (spectro-photometric exoplanet research) extreme-AO planet hunter saw first light at the VLT observatory on Mount Paranal in May 2014 after ten years of development. Great efforts were put into modelling its performance, particularly in terms of achievable contrast, and to budgeting instrumental features such as wave front errors and optical transmission to each of the instrument's three focal planes, the near infrared dual imaging camera IRDIS, the near infrared integral field spectrograph IFS and the visible polarimetric camera ZIMPOL. In this paper we aim at comparing predicted performance with measured performance. In addition to comparing on-sky contrast curves and calibrated transmission measurements, we also compare the PSD-based wave front error budget with in-situ wave front maps obtained thanks to a Zernike phase mask, ZELDA, implemented in the infrared coronagraph wheel. One of the most critical elements of the SPHERE system is its high-order deformable mirror, a prototype 40x40 actuator piezo stack design developed in parallel with the instrument itself. The development was a success, as witnessed by the instrument performance, in spite of some bad surprises discovered on the way. The devastating effects of operating without taking properly into account the loss of several actuators and the thermally and temporally induced variations in the DM shape will be analysed, and the actions taken to mitigate these defects through the introduction of specially designed Lyot stops and activation of one of the mirrors in the optical train will be described.

Context. Stellar activity is currently challenging the detection of young planets via the radial velocity (RV) technique. Aims. We attempt to definitively discriminate the nature of the RV variations for the young active K5 star BD+20 1790, for which visible (VIS) RV measurements show divergent results on the existence of a substellar companion. Methods. We compare VIS data with high precision RVs in the near infrared (NIR) range by using the GIANO - B and IGRINS spectrographs. In addition, we present for the first time simultaneous VIS-NIR observations obtained with GIARPS (GIANO - B and HARPS - N) at Telescopio Nazionale Galileo (TNG). Orbital RVs are achromatic, so the RV amplitude does not change at different wavelengths, while stellar activity induces wavelength-dependent RV variations, which are significantly reduced in the NIR range with respect to the VIS. Results. The NIR radial velocity measurements from GIANO - B and IGRINS show an average amplitude of about one quarter with respect to previously published VIS data, as expected when the RV jitter is due to stellar activity. Coeval multi-band photometry surprisingly shows larger amplitudes in the NIR range, explainable with a mixture of cool and hot spots in the same active region. Conclusions. In this work, the claimed massive planet around BD+20 1790 is ruled out by our data. We exploited the crucial role of multi- wavelength spectroscopy when observing young active stars: thanks to facilities like GIARPS that provide simultaneous observations, this method can reach its maximum potential.

Coronal mass ejections (CMEs) and other solar eruptive phenomena can be physically linked by combining data from a multitude of ground-based and space-based instruments alongside models, however this can be challenging for automated operational systems. The EU Framework Package 7 HELCATS project provides catalogues of CME observations and properties from the Helio- spheric Imagers onboard the two NASA/STEREO spacecraft in order to track the evolution of CMEs in the inner heliosphere. From the main HICAT catalogue of over 2,000 CME detections, an automated algorithm has been developed to connect the CMEs observed by STEREO to any corresponding solar flares and active region (AR) sources on the solar surface. CME kinematic properties, such as speed and angular width, are compared with AR magnetic field properties, such as magnetic flux, area, and neutral line characteristics. The resulting LOWCAT catalogue is also compared to the extensive AR property database created by the EU Horizon 2020 FLARECAST project, which provides more complex magnetic field parameters derived from vector magnetograms. Initial statistical analysis has been undertaken on the new data to provide insight into the link between flare and CME events, and characteristics of eruptive ARs. Warning thresholds determined from analysis of the evolution of these parameters is shown to be a useful output for operational space weather purposes. Parameters of particular interest for further analysis include total unsigned flux, vertical current, and current helicity. The automated method developed to create the LOWCAT catalogue may also be useful for future efforts to develop operational CME forecasting.

We propose a forecasting approach for solar flares based on data from Solar Cycle 24, taken by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) mission. In particular, we use the Space-weather HMI Active Region Patches (SHARP) product that facilitates cut-out magnetograms of solar active regions (AR) in the Sun in near-real-time (NRT), taken over a five-year interval (2012 - 2016). Our approach utilizes a set of thirteen predictors, which are not included in the SHARP metadata, extracted from line-of-sight and vector photospheric magnetograms. We exploit several Machine Learning (ML) and Conventional Statistics techniques to predict flares of peak magnitude >M1 and >C1, within a 24 h forecast window. The ML methods used are multi-layer perceptrons (MLP), support vector machines (SVM) and random forests (RF). We conclude that random forests could be the prediction technique of choice for our sample, with the second best method being multi-layer perceptrons, subject to an entropy objective function. A Monte Carlo simulation showed that the best performing method gives accuracy ACC=0.93(0.00), true skill statistic TSS=0.74(0.02) and Heidke skill score HSS=0.49(0.01) for >M1 flare prediction with probability threshold 15% and ACC=0.84(0.00), TSS=0.60(0.01) and HSS=0.59(0.01) for >C1 flare prediction with probability threshold 35%.

The effect of using two representations of the normal-to-surface magnetic field to calculate photospheric measures that are related to active region (AR) potential for flaring is presented. Several AR properties were computed using line-of-sight ($B_{\rm los}$) and spherical-radial ($B_{r}$) magnetograms from the Space-weather HMI Active Region Patch (SHARP) products of the \emph{Solar Dynamics Observatory}, characterizing the presence and features of magnetic polarity inversion lines, fractality, and magnetic connectivity of the AR photospheric field. The data analyzed corresponds to $\approx$ 4,000 AR observations, achieved by randomly selecting 25 % of days between September 2012 and May 2016 for analysis at 6-hr cadence. Results from this statistical study include: i) the $B_{r}$ component results in a slight upwards shift of property values in a manner consistent with a field-strength underestimation by the $B_{\rm los}$ component; ii) using the $B_{r}$ component results in significantly lower inter-property correlation in one-third of the cases, implying more independent information about the state of the AR photospheric magnetic field; iii) flaring rates for each property vary between the field components in a manner consistent with the differences in property-value ranges resulting from the components; iv) flaring rates generally increase for higher values of properties, except Fourier spectral power index that has flare rates peaking around a value of $5/3$. These findings indicate that there may be advantages in using $B_{r}$ rather than $B_{\rm los}$ in calculating flare-related AR magnetic properties, especially for regions located far from central meridian.

We explore the association of non-neutralized currents with solar flare occurrence in a sizable sample of observations, aiming to show the potential of such currents in solar flare prediction. We use the regularly produced high quality vector magnetograms by the Helioseismic Magnetic Imager and more specifically the Space weather HMI Active Region Patches (SHARP). Through a newly established method, that incorporates detailed error analysis, we calculate the non-neutralized currents contained in active regions (AR). Two predictors are produced, namely the total and the maximum unsigned non-neutralized current. Both are tested in AR time-series and a representative sample of point-in-time observations during the interval 2012-2016. The average values of non-neutralized currents in flaring active regions are by more than an order of magnitude higher than in non-flaring regions and correlate very well with the corresponding flare index. The temporal evolution of these parameters appears to be connected to physical processes, such as flux emergence and/or magnetic polarity inversion line formation that are associated with increased solar flare activity. Using Bayesian inference of flaring probabilities, it is shown that the total unsigned non-neutralized current outperforms significantly the total unsigned magnetic flux and other well established current-related predictors. Thus, it shows good prospects for inclusion in an operational flare forecasting service. We plan to use the new predictor in the framework of the FLARECAST project along with other highly performing predictors.

Since 2012, thanks to the installation of the high resolution echelle spectrograph in the optical range HARPS-N, the Italian telescope TNG (La Palma) became one of the key facilities for the study of the extrasolar planets. In 2014 TNG also offered GIANO to the scientific community, providing a near-infrared (NIR) cross-dispersed echelle spectroscopy covering 0.97 - 2.45 micron at a resolution of 50,000. GIANO, although designed for direct light-feed from the telescope at the Nasmyth-B focus, was provisionally mounted on the rotating building and connected via fibers to only available interface at the Nasmyth-A focal plane. The synergy between these two instruments is particularly appealing for a wide range of science cases, especially for the search for exoplanets around young and active stars and the characterisation of their atmosphere. Through the funding scheme "WOW" (a Way to Others Worlds), the Italian National Institute for Astrophysics (INAF) proposed to position GIANO at the focal station for which it was originally designed and the simultaneous use of these spectrographs with the aim to achieve high-resolution spectroscopy in a wide wavelength range (0.383-2.45 micron) obtained in a single exposure, giving rise to the project called GIARPS (GIANO-B & HARPS-N). Because of its characteristics GIARPS can be considered the first and unique worldwide instrument providing not only high resolution in a large wavelength band, but also a high precision radial velocity measurement both in the visible and in the NIR arm, since in the next future GIANO-B will be equipped with gas absorption cells.

We re-examine the possibility of ultrahigh energy cosmic rays being accelerated in nearby dormant quasars. We particularize our study to heavy nuclei to accommodate the spectrum and nuclear composition recently reported by the Pierre Auger Collaboration. Particle acceleration is driven by the Blandford-Znajek mechanism, which wires the dormant spinning black holes as Faraday unipolar dynamos. We demonstrate that energy losses are dominated by photonuclear interactions on the ambient photon fields. We argue that the local dark fossils of the past quasar activity can be classified on the basis of how source parameters (mass of the central engine and photon background surrounding the accelerator) impact the photonuclear interaction. In this classification it is possible to distinguish two unequivocal type of sources: those in which nuclei are completely photodisintegrated before escaping the acceleration region and those in which photopion production is the major energy damping mechanism. We further argue that the secondary nucleons from the photodisintegrated nuclei (which have a steep spectral index at injection) can populate the energy region below "the ankle" feature in the cosmic ray spectrum, whereas heavy and medium mass nuclei (with a harder spectral index) populate the energy region beyond "the ankle", all the way to the high energy end of the spectrum. In addition, we show that five potential quasar remnants from our cosmic backyard correlate with the hot-spot observed by the Telescope Array.

GIARPS (GIAno & haRPS) is a project devoted to have on the same focal station of the Telescopio Nazionale Galileo (TNG) both the high resolution spectrographs HARPS-N (VIS) and GIANO (NIR) working simultaneously. This could be considered the first and unique worldwide instrument providing cross-dispersed echelle spectroscopy at a high resolution (R=115,000 in the visual and R=50,000 in the IR) and over in a wide spectral range (0.383 - 2.45 micron) in a single exposure. The science case is very broad, given the versatility of such an instrument and the large wavelength range. A number of outstanding science cases encompassing mainly extra-solar planet science starting from rocky planet search and hot Jupiters, atmosphere characterization can be considered. Furthermore both instrument can measure high precision radial velocity by means the simultaneous thorium technique (HARPS - N) and absorbing cell technique (GIANO) in a single exposure. Other science cases are also possible. Young stars and proto-planetary disks, cool stars and stellar populations, moving minor bodies in the solar system, bursting young stellar objects, cataclysmic variables and X-ray binary transients in our Galaxy, supernovae up to gamma-ray bursts in the very distant and young Universe, can take advantage of the unicity of this facility both in terms of contemporaneous wide wavelength range and high resolution spectroscopy.

We present the results from the first ensemble prediction model for major solar flares (M and X classes). The primary aim of this investigation is to explore the construction of an ensemble for an initial prototyping of this new concept. Using the probabilistic forecasts from three models hosted at the Community Coordinated Modeling Center (NASA-GSFC) and the NOAA forecasts, we developed an ensemble forecast by linearly combining the flaring probabilities from all four methods. Performance-based combination weights were calculated using a Monte-Carlo-type algorithm that applies a decision threshold $P_{th}$ to the combined probabilities and maximizing the Heidke Skill Score (HSS). Using the data for 13 recent solar active regions between years 2012 - 2014, we found that linear combination methods can improve the overall probabilistic prediction and improve the categorical prediction for certain values of decision thresholds. Combination weights vary with the applied threshold and none of the tested individual forecasting models seem to provide more accurate predictions than the others for all values of $P_{th}$. According to the maximum values of HSS, a performance-based weights calculated by averaging over the sample, performed similarly to a equally weighted model. The values $P_{th}$ for which the ensemble forecast performs the best are 25 \% for M-class flares and 15 \% for X-class flares. When the human-adjusted probabilities from NOAA are excluded from the ensemble, the ensemble performance in terms of the Heidke score, is reduced.

We present deep near-infrared (NIR) J, Ks photometry of the old, metal-poor Galactic globular cluster M\,15 obtained with images collected with the LUCI1 and PISCES cameras available at the Large Binocular Telescope (LBT). We show how the use of First Light Adaptive Optics system coupled with the (FLAO) PISCES camera allows us to improve the limiting magnitude by ~2 mag in Ks. By analyzing archival HST data, we demonstrate that the quality of the LBT/PISCES color magnitude diagram is fully comparable with analogous space-based data. The smaller field of view is balanced by the shorter exposure time required to reach a similar photometric limit. We investigated the absolute age of M\,15 by means of two methods: i) by determining the age from the position of the main sequence turn-off; and ii) by the magnitude difference between the MSTO and the well-defined knee detected along the faint portion of the MS. We derive consistent values of the absolute age of M15, that is 12.9+-2.6 Gyr and 13.3+-1.1 Gyr, respectively.

We present an overview of and status report on the WEAVE next-generation spectroscopy facility for the William Herschel Telescope (WHT). WEAVE principally targets optical ground-based follow up of upcoming ground-based (LOFAR) and space-based (Gaia) surveys. WEAVE is a multi-object and multi-IFU facility utilizing a new 2-degree prime focus field of view at the WHT, with a buffered pick-and-place positioner system hosting 1000 multi-object (MOS) fibres, 20 integral field units, or a single large IFU for each observation. The fibres are fed to a single spectrograph, with a pair of 8k(spectral) x 6k (spatial) pixel cameras, located within the WHT GHRIL enclosure on the telescope Nasmyth platform, supporting observations at R~5000 over the full 370-1000nm wavelength range in a single exposure, or a high resolution mode with limited coverage in each arm at R~20000. The project is now in the final design and early procurement phase, with commissioning at the telescope expected in 2017.

We study structure and dynamics of the turbulent photospheric magnetic field in the active region NOAA 11158 by characterizing the spatial and temporal scaling properties of the line-of-sight (LOS) component. Using high-resolution high-cadence LOS magnetograms from SDO/HMI, we measured power-law exponents $\alpha$ and $\beta$ describing wavenumber- ($k$) and frequency-domain ($f$) Fourier power spectra, respectively, and investigated their evolution during the passage of the active region through the field of view of HMI. Flaring active region NOAA 11158 produces an average one-dimensional spatial power spectral density that follows approximately a $k^{-2}$ power law -- a spectrum that suggests parallel MHD fluctuations in an anisotropic turbulent medium. In addition, we found that values of $\alpha$ capture systematic changes in the configuration of LOS photospheric magnetic field during flaring activity in the corona. Position-dependent values of the temporal scaling exponent $\beta$ showed that, on average, the core of the active region scales with $\beta >$ 3 surrounded by a diffusive region with an approximately $f^{-2}$-type spectrum. Our results indicate that only about 1 - 3 \% of the studied LOS photospheric magnetic flux displays $\beta\approx\alpha$, implying that Taylor's hypothesis of frozen-in-flow turbulence is typically invalid for this scalar in presence of turbulent photospheric flows. In consequence, both spatial and temporal variations of the plasma and magnetic field must be included in a complete description of the turbulent evolution of active regions.

We have performed H and Ks band observations of the planetary system around HR 8799 using the new AO system at the Large Binocular Telescope and the PISCES Camera. The excellent instrument performance (Strehl ratios up to 80% in H band) enabled detection the inner planet HR8799e in the H band for the first time. The H and Ks magnitudes of HR8799e are similar to those of planets c and d, with planet e slightly brighter. Therefore, HR8799e is likely slightly more massive than c and d. We also explored possible orbital configurations and their orbital stability. We confirm that the orbits of planets b, c and e are consistent with being circular and coplanar; planet d should have either an orbital eccentricity of about 0.1 or be non-coplanar with respect to b and c. Planet e can not be in circular and coplanar orbit in a 4:2:1 mean motion resonances with c and d, while coplanar and circular orbits are allowed for a 5:2 resonance. The analysis of dynamical stability shows that the system is highly unstable or chaotic when planetary masses of about 5 MJup for b and 7 MJup for the other planets are adopted. Significant regions of dynamical stability for timescales of tens of Myr are found when adopting planetary masses of about 3.5, 5, 5, and 5 Mjup for HR 8799 b, c, d, and e respectively. These masses are below the current estimates based on the stellar age (30 Myr) and theoretical models of substellar objects.

We present diffraction-limited \ks band and \lprime adaptive optics images of the edge-on debris disk around the nearby F2 star HD 15115, obtained with a single 8.4 m primary mirror at the Large Binocular Telescope. At \ks band the disk is detected at signal-to-noise per resolution element (SNRE) \about 3-8 from \about 1-2\fasec 5 (45-113 AU) on the western side, and from \about 1.2-2\fasec 1 (63-90 AU) on the east. At \lprime the disk is detected at SNRE \about 2.5 from \about 1-1\fasec 45 (45-90 AU) on both sides, implying more symmetric disk structure at 3.8 \microns . At both wavelengths the disk has a bow-like shape and is offset from the star to the north by a few AU. A surface brightness asymmetry exists between the two sides of the disk at \ks band, but not at \lprime . The surface brightness at \ks band declines inside 1\asec (\about 45 AU), which may be indicative of a gap in the disk near 1\asec. The \ks - \lprime disk color, after removal of the stellar color, is mostly grey for both sides of the disk. This suggests that scattered light is coming from large dust grains, with 3-10 \microns -sized grains on the east side and 1-10 \microns dust grains on the west. This may suggest that the west side is composed of smaller dust grains than the east side, which would support the interpretation that the disk is being dynamically affected by interactions with the local interstellar medium.

As the only directly imaged multiple planet system, HR 8799 provides a unique opportunity to study the physical properties of several planets in parallel. In this paper, we image all four of the HR 8799 planets at H-band and 3.3 microns with the new LBT adaptive optics system, PISCES, and LBTI/LMIRCam. Our images offer an unprecedented view of the system, allowing us to obtain H and 3.3$ micron photometry of the innermost planet (for the first time) and put strong upper-limits on the presence of a hypothetical fifth companion. We find that all four planets are unexpectedly bright at 3.3 microns compared to the equilibrium chemistry models used for field brown dwarfs, which predict that planets should be faint at 3.3 microns due to CH4 opacity. We attempt to model the planets with thick-cloudy, non-equilibrium chemistry atmospheres, but find that removing CH4 to fit the 3.3 micron photometry increases the predicted L' (3.8 microns) flux enough that it is inconsistent with observations. In an effort to fit the SED of the HR 8799 planets, we construct mixtures of cloudy atmospheres, which are intended to represent planets covered by clouds of varying opacity. In this scenario, regions with low opacity look hot and bright, while regions with high opacity look faint, similar to the patchy cloud structures on Jupiter and L/T transition brown-dwarfs. Our mixed cloud models reproduce all of the available data, but self-consistent models are still necessary to demonstrate their viability.

The Large Binocular Telescope (LBT) is a unique telescope featuring two co-mounted optical trains with 8.4m primary mirrors. The telescope Adaptive Optics (AO) system uses two innovative key components, namely an adaptive secondary mirror with 672 actuators and a high-order pyramid wave-front sensor. During the on-sky commissioning such a system reached performances never achieved before on large ground-based optical telescopes. Images with 40mas resolution and Strehl Ratios higher than 80% have been acquired in H band (1.6 micron). Such images showed a contrast as high as 10e-4. Based on these results, we compare the performances offered by a Natural Guide Star (NGS) system upgraded with the state-of-the-art technology and those delivered by existing Laser Guide Star (LGS) systems. The comparison, in terms of sky coverage and performances, suggests rethinking the current role ascribed to NGS and LGS in the next generation of AO systems for the 8-10 meter class telescopes and Extremely Large Telescopes (ELTs).

The new 8.4m LBT adaptive secondary AO system, with its novel pyramid wavefront sensor, was used to produce very high Strehl (75% at 2.16 microns) near infrared narrowband (Br gamma: 2.16 microns and [FeII]: 1.64 microns) images of 47 young (~1 Myr) Orion Trapezium theta1 Ori cluster members. The inner ~41x53" of the cluster was imaged at spatial resolutions of ~0.050" (at 1.64 microns). A combination of high spatial resolution and high S/N yielded relative binary positions to ~0.5 mas accuracies. Including previous speckle data, we analyse a 15 year baseline of high-resolution observations of this cluster. We are now sensitive to relative proper motions of just ~0.3 mas/yr (0.6 km/s at 450 pc) this is a ~7x improvement in orbital velocity accuracy compared to previous efforts. We now detect clear orbital motions in the theta1 Ori B2/B3 system of 4.9+/-0.3 km/s and 7.2+/-0.8 km/s in the theta1 Ori A1/A2 system (with correlations of PA vs. time at >99% confidence). All five members of the theta1 Ori B system appear likely as a gravitationally bound "mini-cluster". The very lowest mass member of the theta1 Ori B system (B4; mass ~0.2 Msun) has, for the first time, a clearly detected motion (at 4.3+/-2.0 km/s; correlation=99.7%) w.r.t B1. However, B4 is most likely in an long-term unstable (non-hierarchical) orbit and may "soon" be ejected from this "mini-cluster". This "ejection" process could play a major role in the formation of low mass stars and brown dwarfs.