The uncharted territory of the high-redshift ($z\gtrsim3$) Universe holds the key to understand the evolution of quasars. In an attempt to identify the most extreme members of the quasar population, i.e., blazars, we have carried out a multi-wavelength study of a large sample of radio-loud quasars beyond $z=3$. Our sample consists of 9 $\gamma$-ray detected blazars and 133 candidate blazars selected based on the flatness of their soft X-ray spectra (0.3$-$10 keV photon index $\leq1.75$), including 15 with NuSTAR observations. The application of the likelihood profile stacking technique reveals that the high-redshift blazars are faint $\gamma$-ray emitters with steep spectra. The high-redshift blazars host massive black holes ($\langle \log~M_{\rm BH,~M{\odot}} \rangle>9$) and luminous accretion disks ($\langle L_{\rm disk} \rangle>10^{46}$ erg s$^{-1}$). Their broadband spectral energy distributions are found to be dominated by high-energy radiation indicating their jets to be among the most luminous ones. Focusing on the sources exhibiting resolved X-ray jets (as observed with the Chandra satellite), we find the bulk Lorentz factor to be larger with respect to other $z>3$ blazars, indicating faster moving jets. We conclude that the presented list of the high-redshift blazars may act as a reservoir for follow-up observations, e.g., with NuSTAR, to understand the evolution of relativistic jets at the dawn of the Universe.

BL Lacertae objects (BL Lacs) represent a large fraction (22\%) of gamma-ray sources in the Third Fermi Large Area Telescope catalog (3FGL). Nearly half of the BL Lac population remains without a redshift because of their featureless optical spectra. We aim to increase the number of BL Lacs with a redshift measurement by using the photometric technique. For this work, we have used 6 Swift-UVOT filters and SDSS g,r,i,z optical filters mounted on the 0.65 m SARA-CTIO located in Chile and the 1.0 m SARA-ORM in Canary Islands. A sample of 45 sources was selected from the 3FGL catalog for which photometry was performed in 10 optical and UV filters to obtain redshift measurements. We found 3 sources with z>1.3 while reliable upper limits have been derived for 17 sources. The results presented here bring the total number of high-z (z>1.3) BL Lacs to 29.

High synchrotron peak (HSP; $\nu_{sy}^{pk} > 10^{15}$ Hz) BL Lacs are some of the most extreme accelerators in the Universe. Those found at high redshifts ($z>1$) challenge our understanding of blazar evolution models and are crucial for cosmological measurements of the Extragalactic Background Light. In this paper, we study a high-$z$ BL Lac, 4FGL J2146.5-1344, detected to be at $z$=1.34 using the photometric dropout technique. We collected multi-wavelength data for this source from optical up to $\gamma$-rays, in order to study its spectral energy distribution (SED). In particular, this source was observed for the first time with NuSTAR, which accurately measures the synchrotron emission of this blazar up to 50 keV. Despite being classified as an HSP BL Lac object, the modeling of the SED reveals that this source likely belongs to the "masquerading BL Lac" class, which comprises of FSRQs appearing as disguised BL Lac objects.

With bolometric luminosities exceeding $10^{48}$ erg s$^{-1}$, powerful jets and supermassive black holes at their center, MeV blazars are some of the most extreme sources in the Universe. Recently, the Fermi-Large Area Telescope detected five new $\gamma$-ray emitting MeV blazars beyond redshift $z=3.1$. With the goal of precisely characterizing the jet properties of these extreme sources, we started a multiwavelength campaign to follow them up with joint NuSTAR, Swift and SARA observations. We observe six high-redshift quasars, four of them belonging to the new $\gamma$-ray emitting MeV blazars. Thorough X-ray analysis reveals spectral flattening at soft X-ray for three of these objects. The source NVSS J151002$+$570243 also shows a peculiar re-hardening of the X-ray spectrum at energies $E>6\,\rm keV$. Adopting a one-zone leptonic emission model, this combination of hard X-rays and $\gamma$-rays enables us to determine the location of the Inverse Compton peak and to accurately constrain the jet characteristics. In the context of the jet-accretion disk connection, we find that all six sources have jet powers exceeding accretion disk luminosity, seemingly validating this positive correlation even beyond $z>3$. Our six sources are found to have $10^9 \rm M_{\odot}$ black holes, further raising the space density of supermassive black holes in the redshift bin $z=[3,4]$.

Approximately one-third of the gamma-ray sources in the third Fermi-LAT catalog are unidentified or unassociated with objects at other wavelengths. Observations with Swift-XRT have yielded possible counterparts in $\sim$30% of these source regions. The objective of this work is to identify the nature of these possible counterparts, utilizing their gamma ray properties coupled with the Swift derived X-ray properties. The majority of the known sources in the Fermi catalogs are blazars, which constitute the bulk of the extragalactic gamma-ray source population. The galactic population on the other hand is dominated by pulsars. Blazars and pulsars occupy different parameter space when X-ray fluxes are compared with various gamma-ray properties. In this work, we utilize the X-ray observations performed with the Swift-XRT for the unknown Fermi sources and compare their X-ray and gamma-ray properties to differentiate between the two source classes. We employ two machine learning algorithms, decision tree and random forest classifier, to our high signal-to-noise ratio sample of 217 sources, each of which correspond to Fermi unassociated regions. The accuracy score for both methods were found to be 97% and 99%, respectively. The random forest classifier, which is based on the application of a multitude of decision trees, associated a probability value (P$_{bzr}$) for each source to be a blazar. This yielded 173 blazar candidates with P$_{bzr}$ $\geq$ 90% for each of these sources, and 134 of these possible blazar source associations had P$_{bzr}$ $\geq$ 99%. The results yielded 13 sources with P$_{bzr}$ $\leq$ 10%, which we deemed as reasonable candidates for pulsars, 7 of which result with P$_{bzr}$ $\leq$ 1%. There were 31 sources that exhibited intermediate probabilities and were termed ambiguous due to their unclear characterization as a pulsar or a blazar.

Experience suggests that structural issues in how institutional Astrophysics approaches data-driven science and the development of discovery technology may be hampering the community's ability to respond effectively to a rapidly changing environment in which increasingly complex, heterogeneous datasets are challenging our existing information infrastructure and traditional approaches to analysis. We stand at the confluence of a new epoch of multimessenger science, remote co-location of data and processing power and new observing strategies based on miniaturized spacecraft. Significant effort will be required by the community to adapt to this rapidly evolving range of possible discovery moduses. In the suggested creation of a new Astrophysics element, Advanced Astrophysics Discovery Technology, we offer an affirmative solution that places the visibility of discovery technologies at a level that we suggest is fully commensurate with their importance to the future of the field.

The common observations of multiple radio VLBI stationary knots in high-frequency-peaked BL Lacs (HBLs) can be interpreted as multiple recollimation shocks accelerating particles along jets. This approach can resolve the so-called "bulk Lorentz factor crisis" of sources with high Lorentz factor, deduced from maximum gamma-gamma opacity and fast variability, and apparently inconsistent slow/stationary radio knots. It also suggests that a unique pattern of the non-thermal emission variability should appear after each strong flare. Taking advantage of the 13 years of observation of the HBL Mrk 421 by the X-ray Telescope on the Neil Gehrels Swift Observatory Swift-XRT, we probe for such an intrinsic variability pattern. Its significance is then statistically estimated via comparisons with numerous similar simulated lightcurves. A suggested variability pattern is identified, consistent with a main flare emission zone located in the most upstream 15.3 GHz radio knot at 0.38 mas from the core. Subsequent flux excesses in the lightcurve are consistent with a perturbation crossing all the downstream radio knots with a constant apparent speed of 45 c. The significance of the observed variability pattern not arising from stochastic processes is found above 3 standard deviations, opening a promising path for further investigations in other blazars and with other energy bands. In addition to highlight the role of stationary radio knots as high-energy particle accelerators in jets, the developed method allows estimates of the apparent speed and size of a jet perturbation without the need to directly observe any motion in jets.

We present the results of an identification campaign of unassociated sources from the Fermi Large Area Telescope 3FHL catalog. Out of 200 unidentified sources, we selected 110 sources for which archival Swift-XRT observations were available, 52 of which were found to have exactly one X-ray counterpart within the 3FHL 95\% positional uncertainty. In this work, we report the X-ray, optical, IR, and radio properties of these 52 sources using positional associations with objects in various catalogs. The Wide-field Infrared Survey Explorer color-color plot for sources suggests that most of these belong to the blazar class family. The redshift measurements for these objects range from z=0.277 to z=2.1. Additionally, under the assumption that the majority of these sources are blazars, three machine-learning algorithms are employed to classify the sample into flat spectrum radio quasars or BL Lacertae objects. These suggest that the majority of the previously unassociated sources are BL Lac objects, in agreement with the fact the BL Lac objects represent by far the most numerous population detected above 10 GeV in 3FHL.

We report the first results obtained from our campaign to characterize the intranight-optical variability (INOV) properties of {\it Fermi}~detected blazars, using the observations from the recently commissioned 1.3 m J C Bhattacharya telescope (JCBT). During the first run, we were able to observe 17 blazars in the Bessel $R$ filter for $\sim$137 hrs. Using $C$ and scaled $F$-statistics, we quantify the extent of INOV and derive the duty cycle (DC) which is the fraction of time during which a source exhibits a substantial flux variability. We find a high DC of 40\% for BL Lac objects and the flat spectrum radio quasars are relatively less variable (DC $\sim15$\%). However, when estimated for blazars sub-classes, a high DC of $\sim$59\% is found in low synchrotron peaked (LSP) blazars, whereas, intermediate and high synchrotron peaked objects have a low DC of $\sim$11\% and 13\%, respectively. We find evidences about the association of the high amplitude INOV with the $\gamma$-ray flaring state. We also notice a high polarization during the elevated INOV states (for the sources that have polarimetric data available), thus supporting the jet based origin of the observed variability. We plan to enlarge the sample and utilize the time availability from the small telescopes, such as 1.3 m JCBT, to strengthen/verify the results obtained in this work and those existed in the literature.

We present the results of the optical spectroscopy follow-up of a sample of 28 unclassified blazars from the Third Fermi-LAT Catalog of High-Energy Sources (3FHL). All the spectra were taken with the 4m Mayall telescope at Kitt Peak. With this follow-up program we are able to classify 27 out of 28 objects as BL Lacs, while the remaining one is a flat spectrum radio quasar. We determine a redshift (z) for three of these objects and a lower limit on z for other four sources: the farthest object for which we obtain a redshift has z>0.836. These results are part of a more extended campaign of optical spectroscopy follow-up of 3FHL blazars, aimed to obtain a complete sample of blazars at >10 GeV which will then be used to extend our knowledge on blazar emission mechanisms and on the extragalactic background light.

BL Lacertae (BL Lac) objects are the prominent members of the third {\it Fermi} Large Area Telescope catalog of $\gamma$-ray sources. Half of the BL Lac population ($\sim$ 300) lack redshift measurements, which is due to the absence of lines in their optical spectrum, thereby making it difficult to utilize spectroscopic methods. Our photometric drop-out technique can be used to establish the redshift for a fraction of these sources. This work employed 6 filters mounted on the $Swift$-UVOT and 4 optical filters on two telescopes, the 0.65 m SARA-CTIO in Chile and 1.0 m SARA-ORM in the Canary Islands, Spain. A sample of 15 sources was extracted from the $Swift$ archival data for which 6 filter UVOT observations were conducted. By complementing the {\it Swift} observations with the SARA ones, we were able to discover two high redshift sources: 3FGL J1155.4-3417 and 3FGL J1156.7-2250 at $z=1.83^{+0.10}_{-0.13}$ and $z=1.73^{+0.11}_{-0.19}$ , respectively, resulting from the dropouts in the powerlaw template fits to these data. The discoveries add to the important (26 total) sample of high-redshift BL Lacs. While the sample of high-z BL Lacs is still rather small, these objects do not seem to fit well within known schemes of the blazar population and represent the best probes of the extragalactic background light.

Since its discovery in 2008, the Andromeda galaxy nova M31N 2008-12a has been observed in eruption every single year. This unprecedented frequency indicates an extreme object, with a massive white dwarf and a high accretion rate, which is the most promising candidate for the single-degenerate progenitor of a type-Ia supernova known to date. The previous three eruptions of M31N 2008-12a have displayed remarkably homogeneous multi-wavelength properties: (i) From a faint peak, the optical light curve declined rapidly by two magnitudes in less than two days; (ii) Early spectra showed initial high velocities that slowed down significantly within days and displayed clear He/N lines throughout; (iii) The supersoft X-ray source (SSS) phase of the nova began extremely early, six days after eruption, and only lasted for about two weeks. In contrast, the peculiar 2016 eruption was clearly different. Here we report (i) the considerable delay in the 2016 eruption date, (ii) the significantly shorter SSS phase, and (iii) the brighter optical peak magnitude (with a hitherto unobserved cusp shape). Early theoretical models suggest that these three different effects can be consistently understood as caused by a lower quiescence mass-accretion rate. The corresponding higher ignition mass caused a brighter peak in the free-free emission model. The less-massive accretion disk experienced greater disruption, consequently delaying re-establishment of effective accretion. Without the early refueling, the SSS phase was shortened. Observing the next few eruptions will determine whether the properties of the 2016 outburst make it a genuine outlier in the evolution of M31N 2008-12a.

The detection of significant $\gamma$-ray emission from radio-loud narrow line Seyfert 1 (NLSy1s) galaxies enables us to study jets in environments different than those in blazars. However, due to the small number of known $\gamma$-ray emitting NLSy1 ($\gamma$-NLSy1) galaxies, a comprehensive study could not be performed. Here we report the first detection of significant $\gamma$-ray emission from four active galactic nuclei (AGN), recently classified as NLSy1 from their Sloan Digital Sky Survey (SDSS) optical spectrum. Three flat spectrum radio quasars (FSRQs) present in the third Large Area Telescope AGN catalog (3LAC) are also found as $\gamma$-NLSy1 galaxies. Comparing the $\gamma$-ray properties of these objects with 3LAC blazars reveals their spectral shapes to be similar to FSRQs, however, with low $\gamma$-ray luminosity ($\lesssim10^{46-47}$ erg s$^{-1}$). In the Wide-field Infrared Survey Explorer color-color diagram, these objects occupy a region mainly populated by FSRQs. Using the H$_{\beta}$ emission line parameters, we find that on average $\gamma$-NLSy1 have smaller black hole masses than FSRQs at similar redshifts. In the low-resolution SDSS image of one of the $\gamma$-NLSy1 source, we find the evidence of an extended structure. We conclude by noting that overall many observational properties of $\gamma$-NLSy1 sources are similar to FSRQs and therefore, these objects could be their low black hole mass counterparts, as predicted in the literature.

We obtained an optical spectrum of a star we identify as the optical counterpart of the M31 Chandra source CXO J004318.8+412016, because of prominent emission lines of the Balmer series, of neutral helium, and a He II line at 4686 Angstrom. The continuum energy distribution and the spectral characteristics demonstrate the presence of a red giant of K or earlier spectral type, so we concluded that the binary is likely to be a symbiotic system. CXO J004318.8+412016 has been observed in X-rays as a luminous supersoft source (SSS) since 1979, with effective temperature exceeding 40 eV and variable X-ray luminosity, oscillating between a few times 10(35) erg/s and a few times 10(37) erg/s. The optical, infrared and ultraviolet colors of the optical object are consistent with an an accretion disk around a compact object companion, which may either be a white dwarf, or a black hole, depending on the system parameters. If the origin of the luminous supersoft X-rays is the atmosphere of a white dwarf that is burning hydrogen in shell, it is as hot and luminous as post-thermonuclear flash novae, yet no major optical outburst has ever been observed, suggesting that the white dwarf is very massive (m>1.2 M(sol)) and it is accreting and burning at the high rate (mdot>10(-8)M(sol)/year) expected for type Ia supernovae progenitors. In this case, the X-ray variability may be due to a very short recurrence time of only mildly degenerate thermonuclear flashes.

The most powerful sources among the blazar family are MeV blazars. Often detected at $z>2$, they usually display high X- and \gm-ray luminosities, larger-than-average jet powers and black hole masses $\gtrsim 10^9 M_{\odot}$. In the present work we perform a multiwavelength study of three high redshift blazars: 3FGL J0325.5+2223 ($z=2.06$), 3FGL J0449.0+1121 ($z= 2.15$), and 3FGL J0453.2$-$2808 ($z=2.56$), analysing quasi simultaneous data from GROND, \swift-UVOT and XRT, \nustar, and \fermi-LAT. Our main focus is on the hard X-ray band recently unveiled by \nustar~(3$-$79 keV) where these objects show a hard spectrum which enables us to constrain the inverse Compton peak and the jet power. We found that all three targets resemble the most powerful blazars, with the synchrotron peak located in the sub-millimeter range and the inverse Compton peak in the MeV range, and therefore belong to the MeV blazar class. Using a simple one zone leptonic emission model to reproduce the spectral energy distributions, we conclude that a simple combination of synchrotron and accretion disk emission reproduces the infrared-optical spectra while the X-ray to \gm-ray part is well reproduced by the inverse Compton scattering of low energy photons supplied by the broad line region. The black hole masses for each of the three sources are calculated to be $\gtrsim 4 \times 10^{8} M_{\odot}$. The three studied sources have jet power at the level of, or beyond, the accretion luminosity.

Determining redshifts for BL Lacertae (BL Lac) objects using the traditional spectroscopic method is challenging due to the absence of strong emission lines in their optical spectra. We employ the photometric dropout technique to determine redshifts for this class of blazars using the combined 13 broad-band filters from Swift-UVOT and the multi-channel imager GROND at the MPG 2.2 m telescope at ESO's La Silla Observatory. The wavelength range covered by these 13 filters extends from far ultraviolet to the near-Infrared. We report results on 40 new Fermi detected BL Lacs with the photometric redshifts determinations for 5 sources, with 3FGL J1918.2-4110 being the most distant in our sample at z=2.16. Reliable upper limits are provided for 20 sources in this sample. Using the highest energy photons for these Fermi-LAT sources, we evaluate the consistency with the Gamma-ray horizon due to the extragalactic background light.

The Andromeda Galaxy recurrent nova M31N 2008-12a had been observed in eruption ten times, including yearly eruptions from 2008-2014. With a measured recurrence period of $P_\mathrm{rec}=351\pm13$ days (we believe the true value to be half of this) and a white dwarf very close to the Chandrasekhar limit, M31N 2008-12a has become the leading pre-explosion supernova type Ia progenitor candidate. Following multi-wavelength follow-up observations of the 2013 and 2014 eruptions, we initiated a campaign to ensure early detection of the predicted 2015 eruption, which triggered ambitious ground and space-based follow-up programs. In this paper we present the 2015 detection; visible to near-infrared photometry and visible spectroscopy; and ultraviolet and X-ray observations from the Swift observatory. The LCOGT 2m (Hawaii) discovered the 2015 eruption, estimated to have commenced at Aug. $28.28\pm0.12$ UT. The 2013-2015 eruptions are remarkably similar at all wavelengths. New early spectroscopic observations reveal short-lived emission from material with velocities $\sim13000$ km s$^{-1}$, possibly collimated outflows. Photometric and spectroscopic observations of the eruption provide strong evidence supporting a red giant donor. An apparently stochastic variability during the early super-soft X-ray phase was comparable in amplitude and duration to past eruptions, but the 2013 and 2015 eruptions show evidence of a brief flux dip during this phase. The multi-eruption Swift/XRT spectra show tentative evidence of high-ionization emission lines above a high-temperature continuum. Following Henze et al. (2015a), the updated recurrence period based on all known eruptions is $P_\mathrm{rec}=174\pm10$ d, and we expect the next eruption of M31N 2008-12a to occur around mid-Sep. 2016.

We seek clarification of the nature of X-ray sources detected in M 31. Here we focus on CXOM31 J004253.1+411422, the brightness of which suggests that it belongs to the class of ultraluminous X-ray sources. We determine the X-ray properties of sources detected in the XMM-Newton Chandra monitoring program. We investigate spectral properties and search for periodic or quasi-periodic oscillations. A multi-component model is applied to the spectra obtained from XMM-Newton data to evaluate the relative contributions from thermal and non-thermal emission. The time dependence of this ratio is evaluated over a period of forty days. We simultaneously fit data from XMM-Newton EPIC-pn, MOS1 and MOS2 detectors with (non-thermal) powerlaw and (thermal) multicolored blackbody. The X-ray spectrum is best fit by the combination of a thermal component with kT ~ 1 keV and a powerlaw component with photon index approximately 2.6. From combined analysis of Chandra, Swift and XMM-Newton data, the unabsorbed total luminosity of this source decreases from ~ 3.8 x 10^{39} erg s^{-1} in the first observation to ~ 0.5 x 10^{39} ergs s^{-1} over a period of three months. The decay closely follows an exponential decline with a time constant of 32 days. The source spectrum evolves significantly, exhibiting a faster decline of the thermal component. We do not find evidence of any significant temporal features in the power density spectrum. The presence of a thermal component at kT ~ 1 keV in conjunction with a non-thermal high energy tail, is also consistent with spectral properties of other ULXs in the "high state". Our analysis indicates that the underlying source of this first ULX in M~31 is a black hole of mass, M > 13 M_{sun}, accreting near the Eddington limit, that underwent a transient outburst followed by an exponential decay reminiscent of transients associated with galactic X-ray novae.

The Arecibo L-Band Feed Array Zone of Avoidance (ALFA ZOA) Survey has discovered a nearby galaxy, ALFA ZOA J1952+1428, at a heliocentric velocity of +279 km s-1. The galaxy was discovered at low Galactic latitude by 21-cm emission from neutral hydrogen (Hi). We have obtained follow-up observations with the EVLA and the 0.9-m SARA optical telescope. The Hi distribution overlaps an uncataloged, potential optical counterpart. The Hi linear size is 1.4 kpc at our adopted distance of D = 7 Mpc, but the distance estimate is uncertain as Hubble's law is unreliable at low recessional velocities. The optical counterpart has mB = 16.9 mag and B - R = 0.1 mag. These characteristics, including MHI = 107.0 M\odot and LB = 107.5 L\odot, if at 7 Mpc, indicate that this galaxy is a blue compact dwarf, but this remains uncertain until further follow-up observations are complete. Optical follow-up observations are ongoing and near infrared follow-up observations have been scheduled.