The lines of HOC$^+$, HCO and CO$^+$ are considered good tracers of photon-dominated regions (PDRs) and X-ray dominated regions. We study these tracers towards regions of the Sgr B2 cloud selected to be affected by different heating mechanisms. We find the lowest values of the column density ratios of HCO$^+$ versus HOC$^+$, HCO and CO$^+$ in dense HII gas, where UV photons dominate the heating and chemistry of gas. HOC$^+$, HCO and CO$^+$ abundances and the above ratios are compared with those of chemical modeling, finding that high temperature chemistry, a cosmic-ray ionization rate of 10$^{-16}$ s$^{-1}$ and timescales $>$10$^{5.0}$ years explain well the HOC$^+$ abundances in quiescent Sgr B2 regions, while shocks are also needed to explain the highest HCO abundances derived for these regions. CO$^+$ is mainly formed in PDRs since the highest CO$^+$ abundances of $\sim$(6-10)$\times$10$^{-10}$ are found in HII regions with electron densities $>$540 cm$^{-3}$ and that CO$^+$ emission is undetected in quiescent gas. Between the ratios, the HCO$^+$/HCO ratio is sensitive to the electron density as it shows different values in dense and diffuse HII regions. We compare SiO J=2-1 emission maps of Sgr B2 with X-ray maps from 2004 and 2012. One known spot shown on the 2012 X-ray map is likely associated with molecular gas at velocities of 15-25 km s$^{-1}$. We also derive the X-ray ionization rate of $\sim$10$^{-19}$ s$^{-1}$ for Sgr B2 regions pervaded by X-rays in 2004, which is quite low to affect the chemistry of the molecular gas.

We report an astrochemical study on the evolution of interstellar molecular clouds and consequent star formation in the center of the barred spiral galaxy M83. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to image molecular species indicative of shocks (SiO, CH$_3$OH), dense cores (N$_2$H$^+$), and photodissociation regions (CN and CCH), as well as a radio recombination line (H41$\alpha$) tracing active star-forming regions. M83 has a circumnuclear gas ring that is joined at two areas by gas streams from the leading-edge gas lanes on the bar. We found elevated abundances of the shock and dense-core tracers in one of the orbit-intersecting areas, and found peaks of CN and H41$\alpha$ downstream. At the other orbit-intersection area, we found similar enhancement of the shock tracers, but less variation of other tracers, and no sign of active star formation in the stream. We propose that the observed chemical variation or lack of it is due to the presence or absence of collision-induced evolution of molecular clouds and induced star formation. This work presents the most clear case of the chemical evolution in the circumnuclear rings of barred galaxies, thanks to the ALMA resolution and sensitivity.

We have mapped the G287.84-0.82 cometary globule (with the Treasure Chest cluster embedded in it) in the South Pillars region of Carina (i) in [CII], 63micron [OI], and CO(11-10) using upGREAT on SOFIA and (ii) in J=2-1 transitions of CO, 13CO, C18O and J=3-2 transitions of H2CO using the APEX telescope in Chile. We probe the morphology, kinematics, and physical conditions of the molecular gas and the photon dominated regions (PDRs) in G287.84-0.82. The [CII] and [OI] emission suggest that the overall structure of the pillar (with red-shifted photo evaporating tails) is consistent with the effect of FUV radiation and winds from eta-Car and O stars in Trumpler 16. The gas in the head of the pillar is strongly influenced by the embedded cluster, whose brightest member is an O9.5V star, CPD-59 2661. The emission of the [CII] and [OI] lines peak at a position close to the embedded star, while all other tracers peak at another position lying to the north-east consistent with gas being compressed by the expanding PDR created by the embedded cluster. The molecular gas inside the globule is probed with the J=2-1 transitions of CO and isotopologues as well as H2CO, and analyzed using a non-LTE model (escape-probability approach), while we use PDR models to derive the physical conditions of the PDR. We identify at least two PDR gas components; the diffuse part (~10^4 cm^-3) is traced by [CII], while the dense (n~ 2-8x10^5 cm^-3) part is traced by [CII], [OI], CO(11-10). Using the F=2-1 transition of [13CII] detected at 50 positions in the region, we derive optical depths (0.9-5), excitation temperatures of [CII] (80-255 K), and N(C+) of 0.3-1x10^19 cm^-2. The total mass of the globule is ~1000 Msun, about half of which is traced by [CII]. The dense PDR gas has a thermal pressure of 10^7-10^8 K cm^-3, which is similar to the values observed in other regions.

We present the performance of the upGREAT heterodyne array receivers on the SOFIA telescope after several years of operations. This instrument is a multi-pixel high resolution (R > 10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receivers use 7-pixel subarrays configured in a hexagonal layout around a central pixel. The low frequency array receiver (LFA) has 2x7 pixels (dual polarization), and presently covers the 1.83-2.06 THz frequency range, which allows to observe the [CII] and [OI] lines at 158 um and 145 um wavelengths. The high frequency array (HFA) covers the [OI] line at 63 um and is equipped with one polarization at the moment (7 pixels, which can be upgraded in the near future with a second polarization array). The 4.7 THz array has successfully flown using two separate quantum-cascade laser local oscillators from two different groups. NASA completed the development, integration and testing of a dual-channel closed-cycle cryocooler system, with two independently operable He compressors, aboard SOFIA in early 2017 and since then, both arrays can be operated in parallel using a frequency separating dichroic mirror. This configuration is now the prime GREAT configuration and has been added to SOFIA's instrument suite since observing cycle 6.

We report the first detection in the interstellar medium of the Z-isomer of cyanomethanimine (HNCHCN), an HCN dimer proposed as precursor of adenine. We identified six transitions of Z-cyanomethanimine, along with five transitions of E-cyanomethanimine, using IRAM 30m observations towards the Galactic Center quiescent molecular cloud G+0.693. The Z-isomer has a column density of (2.0$\pm$0.6)$\times$10$^{14}$ cm$^{-2}$ and an abundance of 1.5$\times$10$^{-9}$. The relative abundance ratio between the isomers is [Z/E]$\sim$6. This value cannot be explained by the two chemical formation routes previously proposed (gas-phase and grain surface), which predicts abundances ratios between 0.9 and 1.5. The observed [Z/E] ratio is in good agreement with thermodynamic equilibrium at the gas kinetic temperature (130$-$210 K). Since isomerization is not possible in the ISM, the two species may be formed at high temperature. New chemical models, including surface chemistry on dust grains and gas-phase reactions, should be explored to explain our findings. Whatever the formation mechanism, the high abundance of Z-HNCHCN shows that precursors of adenine are efficiently formed in the ISM.

We present an unbiased spectral line survey toward the Galactic Centre (GC) quiescent giant molecular cloud (QGMC), G+0.693 using the GBT and IRAM 30$\,$ telescopes. Our study highlights an extremely rich organic inventory of abundant amounts of nitrogen (N)-bearing species in a source without signatures of star formation. We report the detection of 17 N-bearing species in this source, of which 8 are complex organic molecules (COMs). A comparison of the derived abundances relative to H$_2$ is made across various galactic and extragalactic environments. We conclude that the unique chemistry in this source is likely to be dominated by low-velocity shocks with X-rays/cosmic rays also playing an important role in the chemistry. Like previous findings obtained for O-bearing molecules, our results for N-bearing species suggest a more efficient hydrogenation of these species on dust grains in G+0.693 than in hot cores in the Galactic disk, as a consequence of the low dust temperatures coupled with energetic processing by X-ray/cosmic ray radiation in the GC.

We present SOFIA-upGREAT observations of [CII] emission of Infrared Dark Cloud (IRDC) G035.39-00.33, designed to trace its atomic gas envelope and thus test models of the origins of such clouds. Several velocity components of [CII] emission are detected, tracing structures that are at a wide range of distances in the Galactic plane. We find a main component that is likely associated with the IRDC and its immediate surroundings. This strongest emission component has a velocity similar to that of the $^{13}$CO(2-1) emission of the IRDC, but offset by $\sim3\:{\rm km\:s}^{-1}$ and with a larger velocity width of $\sim9\:{\rm km\:s}^{-1}$. The spatial distribution of the [CII] emission of this component is also offset predominantly to one side of the dense filamentary structure of the IRDC. The CII column density is estimated to be of the order of $\sim10^{17}-10^{18}\,{\rm cm}^{-2}$. We compare these results to the [CII] emission from numerical simulations of magnetized, dense gas filaments formed from giant molecular cloud (GMC) collisions, finding similar spatial and kinematic offsets. These observations and modeling of [CII] add further to the evidence that IRDC G035.39-00.33 has been formed by a process of GMC-GMC collision, which may thus be an important mechanism for initiating star cluster formation.

Phosphorus (P) is one of the essential elements for life due to its central role in biochemical processes. Recent searches have shown that P-bearing molecules (in particular PN and PO) are present in star-forming regions, although their formation routes remain poorly understood. In this Letter, we report observations of PN and PO towards seven molecular clouds located in the Galactic Center, which are characterized by different types of chemistry. PN is detected in five out of seven sources, whose chemistry is thought to be shock-dominated. The two sources with PN non-detections correspond to clouds exposed to intense UV/X-rays/cosmic-ray radiation. PO is detected only towards the cloud G+0.693$-$0.03, with a PO/PN abundance ratio of $\sim$1.5. We conclude that P-bearing molecules likely form in shocked gas as a result of dust grain sputtering, while are destroyed by intense UV/X-ray/cosmic ray radiation.

Aims: To reveal the morphology, chemical composition, kinematics and to establish the main processes prevalent in the gas at the foot points of the giant molecular loops (GMLs) in the Galactic center region Methods: Using the 22-m Mopra telescope, we mapped the M$-3.8+0.9$ molecular cloud, placed at the foot points of a giant molecular loop, in 3-mm range molecular lines. To derive the molecular hydrogen column density, we also observed the $^{13}$CO $(2-1)$ line at 1 mm using the 12-m APEX telescope. From the 3 mm observations 12 molecular species were detected, namely HCO$^+$, HCN, H$^{13}$CN, HNC, SiO, CS, CH$_3$OH, N$_2$H$^+$, SO, HNCO, OCS, and HC$_3$N. Results: Maps revealing the morphology and kinematics of the M$-3.8+0.9$ molecular cloud in different molecules are presented. We identified six main molecular complexes. We derive fractional abundances in 11 selected positions of the different molecules assuming local thermodynamical equilibrium. Conclusions: Most of the fractional abundances derived for the M$-3.8+0.9$ molecular cloud are very similar over the whole cloud. However, the fractional abundances of some molecules show significant difference with respect to those measured in the central molecular zone (CMZ). The abundances of the shock tracer SiO are very similar between the GMLs and the CMZ. The methanol emission is the most abundant specie in the GMLs. This indicates that the gas is likely affected by moderate $\sim $ 30 km s$^{-1}$ or even high velocity (50 km s$^{-1}$) shocks, consistent with the line profile observed toward one of the studied position. The origin of the shocks is likely related to the flow of the gas throughout the GMLs towards the foot points.

Context: The bulk of the Molecular gas in the Central Molecular Zone (CMZ) of the Galactic center region shows warm kinetic temperatures, ranging from $>20$ K in the coldest and densest regions (n$\sim 10^{4-5}$ cm$^{-3}$) up to more than 100 K for densities of about n$\sim 10^3$ cm$^{-3}$. Recently, a more diffuse, hotter ($n \sim 100$ cm$^{-3}$, $T\sim 250$ K) gas component was discovered through absorption observations of H$_3^+$. This component may be widespread in the Galactic center, and low density gas detectable in absorption may be present even outside the CMZ along sightlines crossing the extended bulge of the Galaxy. Aims: We aim to observe and characterize diffuse and low density gas using observations of 3-mm molecular transitions seen in absorption. Methods: Using the Atacama Large (sub)Millimeter Array (ALMA) we observed the absorption against the quasar J1744-312, which is located towards the Galactic bulge region at $(l,b)=(-2{^\circ}.13,-1{^\circ}.0)$, but outside the main molecular complexes. Results: ALMA observations in absorption against the J1744-312 quasar reveal a rich and complex chemistry in low density molecular and presumably diffuse clouds. We detected three velocity components, at $\sim $ 0, $-153$, and $-192$ km s$^{-1}$. The component at $\sim 0$ km s$^{-1}$ could represent gas in the Galactic disk while the velocity components at $-153$, and $-192$ km s$^{-1}$ likely originate from the Galactic bulge. We detected 12 molecules in the survey, but only 7 in the Galactic bulge gas.

Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,\mu\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations.

We present a new multi-pixel high resolution (R >10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receiver uses 2 x 7-pixel subarrays in orthogonal polarization, each in an hexagonal array around a central pixel. We present the first results for this new instrument after commissioning campaigns in May and December 2015 and after science observations performed in May 2016 . The receiver is designed to ultimately cover the full 1.8-2.5 THz frequency range but in its first implementation, the observing range was limited to observations of the [CII] line at 1.9 THz in 2015 and extended to 1.83-2.07 THz in 2016. The instrument sensitivities are state-of-the-art and the first scientific observations performed shortly after the commissioning confirm that the time efficiency for large scale imaging is improved by more than an order of magnitude as compared to single pixel receivers. An example of large scale mapping around the Horsehead Nebula is presented here illustrating this improvement. The array has been added to SOFIA's instrument suite already for ongoing observing cycle 4.

We study the effects of an asymmetric radiation field on the properties of a molecular cloud envelope. We employ observations of carbon monoxide (12CO and 13CO), atomic carbon, ionized carbon, and atomic hydrogen to analyze the chemical and physical properties of the core and envelope of L1599B, a molecular cloud forming a portion of the ring at approximately 27 pc from the star Lambda Ori. The O III star provides an asymmetric radiation field that produces a moderate enhancement of the external radiation field. Observations of the [CII] fine structure line with the GREAT instrument on SOFIA indicate a significant enhanced emission on the side of the cloud facing the star, while the [Ci], 12CO and 13CO J = 1-0 and 2-1, and 12CO J = 3-2 data from the PMO and APEX telescopes suggest a relatively typical cloud interior. The atomic, ionic, and molecular line centroid velocities track each other very closely, and indicate that the cloud may be undergoing differential radial motion. The HI data from the Arecibo GALFA survey and the SOFIA/GREAT [CII] data do not suggest any systematic motion of the halo gas, relative to the dense central portion of the cloud traced by 12CO and 13CO.

The circumnuclear disk (CND) of the Galactic Center is exposed to many energetic phenomena coming from the supermassive black hole Sgr A* and stellar activities. These energetic activities can affect the chemical composition in the CND by the interaction with UV-photons, cosmic-rays, X-rays, and shock waves. We aim to constrain the physical conditions present in the CND by chemical modeling of observed molecular species detected towards it. We analyzed a selected set of molecular line data taken toward a position in the southwest lobe of the CND with the IRAM 30m and APEX 12-meter telescopes and derived the column density of each molecule using a large velocity gradient (LVG) analysis. The determined chemical composition is compared with a time-dependent gas-grain chemical model based on the UCL\_CHEM code that includes the effects of shock waves with varying physical parameters. Molecules such as CO, HCN, HCO$^+$, HNC, CS, SO, SiO, NO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$ and H$_3$O$^+$ are detected and their column densities are obtained. Total hydrogen densities obtained from LVG analysis range between $2 \times 10^4$ and $1 \times 10^6\,$cm$^{-3}$ and most species indicate values around several $\times 10^5\,$cm$^{-3}$, which are lower than values corresponding to the Roche limit, which shows that the CND is tidally unstable. The chemical models show good agreement with the observations in cases where the density is $\sim10^4\,$cm$^{-3}$, the cosmic-ray ionization rate is high, $>10^{-15} \,$s$^{-1}$, or shocks with velocities $> 40\,$km s$^{-1}$ have occurred. Comparison of models and observations favors a scenario where the cosmic-ray ionization rate in the CND is high, but precise effects of other factors such as shocks, density structures, UV-photons and X-rays from the Sgr A* must be examined with higher spatial resolution data.

The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate ${\zeta}_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.

We have mapped the NGC 2023 reflection nebula in [CII] and CO(11--10) with the heterodyne receiver GREAT on SOFIA and obtained slightly smaller maps in 13CO(3--2), CO(3--2), CO(4--3), CO(6--5), and CO(7--6) with APEX in Chile. We use these data to probe the morphology, kinematics, and physical conditions of the C II region, which is ionized by FUV radiation from the B2 star HD37903. The [CII] emission traces an ellipsoidal shell-like region at a position angle of ~ -50 deg, and is surrounded by a hot molecular shell. In the southeast, where the C II region expands into a dense, clumpy molecular cloud ridge, we see narrow and strong line emission from high-J CO lines, which comes from a thin, hot molecular shell surrounding the [CII] emission. The [CII] lines are broader and show photo evaporating gas flowing into the C II region. Based on the strength of the [13CII] F=2--1 line, the [CII] line appears to be somewhat optically thick over most of the nebula with an optical depth of a few. We model the physical conditions of the surrounding molecular cloud and the PDR emission using both RADEX and simple PDR models. The temperature of the CO emitting PDR shell is ~ 90 -- 120 K, with densities of 10^5 -- 10^6 cm^-3, as deduced from RADEX modeling. Our PDR modeling indicates that the PDR layer where [CII] emission dominates has somewhat lower densities, 10^4 to a few times 10^5 cm^-3

We used the IRAM 30m telescope to observe the frequency range [86-116]GHz towards the central regions of the starburst galaxies M83, M82, and NGC253, the AGNs M51, NGC1068, and NGC7469, and the ULIRGs Arp220 and Mrk231. Assuming LTE conditions, we calculated the column densities of 27 molecules and 10 isotopologues. Among others, we report the first tentative detections of CH3CHO, HNCO, and NS in M82 and, for the first time in the extragalactic medium, HC5N in NGC253. Halpha recombination lines were only found in M82 and NGC253. Vibrationally excited lines of HC3N were only detected in Arp220. CH3CCH emission is only seen in the starburst-dominated galaxies. By comparison of the fractional abundances among the galaxies, we looked for the molecules that are best suited to characterise the chemistry of starbursts, AGNs and ULIRGs, as well as the differences among galaxies within the same group.

During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow regulates star formation by e.g. inducing turbulence and heating the surrounding gas. Understanding the associated shocks is therefore essential to the study of star formation. We present comparisons of shock models with CO, H2, and SiO observations in a 'pure' shock position in the BHR71 bipolar outflow. These comparisons provide an insight into the shock and pre-shock characteristics, and allow us to understand the energetic and chemical feedback of star formation on Galactic scales. New CO (Jup = 16, 11, 7, 6, 4, 3) observations from the shocked regions with the SOFIA and APEX telescopes are presented and combined with earlier H2 and SiO data (from the Spitzer and APEX telescopes). The integrated intensities are compared to a grid of models that were obtained from a magneto-hydrodynamical shock code which calculates the dynamical and chemical structure of these regions combined with a radiative transfer module based on the 'large velocity gradient' approximation. The CO emission leads us to update the conclusions of our previous shock analysis: pre-shock densities of 1e4 cm-3 and shock velocities around 20-25 km s-1 are still constrained, but older ages are inferred ( 4000 years). We evaluate the contribution of shocks to the excitation of CO around forming stars. The SiO observations are compatible with a scenario where less than 4% of the pre-shock SiO belongs to the grain mantles. We infer outflow parameters: a mass of 1.8x1e-2 Msun was measured in our beam, in which a momentum of 0.4 Msun km s-1 is dissipated, for an energy of 4.2x1e43erg. We analyse the energetics of the outflow species by species. Comparing our results with previous studies highlights their dependence on the method: H2 observations only are not sufficient to evaluate the mass of outflows.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

Within a few parsecs around the central Black Hole Sgr A*, chemistry in the dense molecular cloud material of the circumnuclear disk (CND) can be affected by many energetic phenomena such as high UV-flux from the massive central star cluster, X-rays from Sgr A*, shock waves, and an enhanced cosmic-ray flux. Recently, spectroscopic surveys with the IRAM 30 meter and the APEX 12 meter telescopes of substantial parts of the 80--500 GHz frequency range were made toward selected positions in and near the CND. These datasets contain lines from the molecules HCN, HCO$^+$, HNC, CS, SO, SiO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$, H$_3$O$^+$ and others. We conduct Large Velocity Gradient analyses to obtain column densities and total hydrogen densities, $n$, for each species in molecular clouds located in the southwest lobe of CND. The data for the above mentioned molecules indicate 10$^5\,$cm$^{-3} \lesssim n <10^6\,$cm$^{-3}$, which shows that the CND is tidally unstable. The derived chemical composition is compared with a chemical model calculated using the UCL\_CHEM code that includes gas and grain reactions, and the effects of shock waves. Models are run for varying shock velocities, cosmic-ray ionization rates, and number densities. The resulting chemical composition is fitted best to an extremely high value of cosmic-ray ionization rate $\zeta \sim 10^{-14}\,$s$^{-1}$, 3 orders of magnitude higher than the value in regular Galactic molecular clouds, if the pre-shock density is $n=10^5\,$cm$^{-3}$.

We aimed to study the molecular composition of the interstellar medium (ISM) surrounding an Active Galactic Nucleus (AGN), by making an inventory of molecular species and their abundances, as well as to establish a chemical differentiation between starburst galaxies and AGN. We used the IRAM-30 m telescope to observe the central 1.5-2 kpc region of NGC1068, covering the frequencies between 86.2 GHz and 115.6 GHz. Using Boltzmann diagrams, we calculated the column densities of the detected molecules. We used a chemical model to reproduce the abundances found in the AGN, to determine the origin of each detected species, and to test the influence of UV fields, cosmic rays, and shocks on the ISM. We identified 24 different molecular species and isotopologues, among which HC3N, SO, N2H+, CH3CN, NS, 13CN, and HN13C are detected for the first time in NGC1068. We obtained the upper limits to the isotopic ratios 12C/13C=49, 16O/18O=177 and 32S/34S=5. Our chemical models suggest that the chemistry in the nucleus of NGC1068 is strongly influenced by cosmic rays, although high values of both cosmic rays and far ultraviolet (FUV) radiation fields also explain well the observations. The gas in the nucleus of NGC1068 has a different chemical composition as compared to starburst galaxies. The distinct physical processes dominating galaxy nuclei (e.g. C-shocks, UV fields, X-rays, cosmic rays) leave clear imprints in the chemistry of the gas, which allow to characterise the nucleus activity by its molecular abundances.

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.

Context: It is well known that the kinetic temperatures, Tkin, of the molecular clouds in the Galactic center region are higher than in typical disk clouds. However, the Tkin of the molecular complexes found at higher latitudes towards the giant molecular loops in the central region of the Galaxy is so far unknown. The gas of these high latitude molecular clouds (hereafter referred to as halo clouds) is located in a region where the gas in the disk may interact with the gas in the halo in the Galactic center region. Aims: To derive Tkin in the molecular clouds at high latitude and understand the physical process responsible for the heating of the molecular gas both in the Central Molecular Zone (the concentration of molecular gas in the inner 500 pc) and in the giant molecular loops. Methods: We measured the metastable inversion transitions of NH3 from (1,1) to (6,6) toward six positions selected throughout the Galactic central disk and halo. We used rotational diagrams and large velocity gradient modeling to estimate the kinetic temperatures toward all the sources. We also observed other molecules like SiO, HNCO, CS, C34S, C18O, and 13CO, to derive the densities and to trace different physical processes (shocks, photodissociation, dense gas) expected to dominate the heating of the molecular gas. Results: We derive for the first time Tkin of the high latitude clouds interacting with the disk in the Galactic center region. We find high rotational temperatures in all the observed positions. We derive two kinetic temperature components (150 K and 40 K) for the positions in the Central Molecular Zone, and only the warm kinetic temperature component for the clouds toward the giant molecular loops. The fractional abundances derived from the different molecules suggest that shocks provide the main heating mechanism throughout the Galactic center, also at high latitudes

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We surveyed the circumnuclear disk of the Seyfert galaxy NGC1068 between the frequencies 86.2 GHz and 115.6 GHz, and identified 17 different molecules. Using a time and depth dependent chemical model we reproduced the observational results, and show that the column densities of most of the species are better reproduced if the molecular gas is heavily pervaded by a high cosmic ray ionization rate of about 1000 times that of the Milky Way. We discuss how molecules in the NGC1068 nucleus may be influenced by this external radiation, as well as by UV radiation fields.

The F-GAMMA program is among the most comprehensive programs that aim at understanding the physics in active galactic nuclei through the multi-frequency monitoring of Fermi blazars. Here we discuss monthly sampled broad-band radio spectra (2.6 - 142 GHz). Two different studies are presented. (a) We discuss that the variability patterns traced can be classified into two classes: (1) to those showing intense spectral-evolution and (2) those showing a self-similar quasi-achromatic behaviour. We show that a simple two-component model can very well reproduce the observed phenomenologies. (b) We present the cm-to-mm behaviour of three gamma-ray bright Narrow Line Seyfert 1 galaxies over time spans varying between ~1.5 and 3 years and compare their variability characteristics with typical blazars.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of Fermi blazars. In the current study we show and discuss the evolution of broad-band radio spectra, which are measured at ten frequencies between 2.64 and 142 GHz using the Effelsberg 100-m and the IRAM 30-m telescopes. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. It is argued that these can be attributed to only two classes of variability mechanisms. The first four types are dominated by spectral evolution and can be described by a simple two-component system composed of: (a) a steep quiescent spectral component from a large scale jet and (b) a time evolving flare component following the "Shock-in-Jet" evolutionary path. The fifth type is characterised by an achromatic change of the broad band spectrum, which could be attributed to a different mechanism, likely involving differential Doppler boosting caused by geometrical effects. Here we present the classification, the assumed physical scenario and the results of calculations that have been performed for the spectral evolution of flares.

We aimed to study the chemistry of the circumnuclear molecular gas of NGC1068, and to compare it with those of the starburst galaxies M82 and NGC253. Using the IRAM-30m telescope, we observed the inner 2 kpc of NGC1068 between 86.2 GHz and 115.6 GHz. We identified 35 spectral features, corresponding to 24 different molecular species. Among them, HC3N, SO, N2H+, CH3CN, NS, 13CN, and HN13C are detected for the first time in NGC1068. Assuming local thermodynamic equilibrium (LTE), we calculated the column densities of the detected molecules, as well as the upper limits to the column densities of some undetected species. The comparison among the chemistries of NGC1068, M82, and NGC253, suggests that, apart from X-rays, shocks also determine the chemistry of NGC1068. We propose the column density ratio between CH3CCH and HC3N as a prime indicator of the imprints of starburst and AGN environments in the circumnuclear interstellar medium. This ratio is, at least, 64 times larger in M82 than in NGC1068, and, at least, 4 times larger in NGC253 than in NGC1068. Finally, we used the UCL_CHEM and UCL_PDR chemical codes to constrain the origin of the species, as well as to test the influence of UV radiation fields and cosmic rays on the observed abundances.

We discuss the parsec-scale structural variability of the extragalactic jet 3C 111 related to a major radio flux density outburst in 2007. The data analyzed were taken within the scope of the MOJAVE, UMRAO, and F-GAMMA programs, which monitor a large sample of the radio brightest compact extragalactic jets with the VLBA, the University of Michigan 26 m, the Effelsberg 100 m, and the IRAM 30 m radio telescopes. The analysis of the VLBA data is performed by fitting Gaussian model components in the visibility domain. We associate the ejection of bright features in the radio jet with a major flux-density outburst in 2007. The evolution of these features suggests the formation of a leading component and multiple trailing components.

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of {\em Fermi} blazars. The current study is concerned with the broad-band radio spectra composed of measurement at ten frequencies between 2.64 and 142 GHz. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. The first four types are dominated by spectral evolution and can be reproduced by a simple two-component system made of the quiescent spectrum of a large scale jet populated with a flaring event evolving according to Marscher & Gear (1985). The last type is characterized by an achromatic change of the broad-band spectrum which must be attributed to a completely different mechanism. Here are presented, the classification, the assumed physical system and the results of simulations that have been conducted.

Aims: A large scale survey of the Galactic center region in he 3 mm rotational transitions of SiO, HCO+ and H13CO+ (beamsize ~ 3. 6 arcmin) was conducted to provide an estimate of cloud conditions, heating mechanisms, chemistry and other properties. Methods: Using the NANTEN 4m telescope from Nagoya University, a region between -5.75<l<5.6 (degree) and -0.68<b<1.3 (degree) was mapped in the J=1-0 lines of HCO+ and H13CO+ and in the J=2-1 line of SiO with a spacing of 3.75 arcmin (HCO+) and 1.875 arcmin (SiO and H13CO+). Results: Velocity channel maps, longitude-velocity maps and latitude-velocity maps are presented. We identify 51 molecular clouds; 33 of them belong to the Galactic center, and 18 to disk gas. We derive an average of the luminosity ratio of SiO(J=2-1)/CO(J=1-0) in clouds belonging to the Galactic center of 4.9x10^{-3} and for disk clouds of 3.4x 10^{-3}. The luminosity ratio of HCO^+(J=1-0)/CO(J=1\to0) in the Galactic Center is 3.5x 10^{-2}, and for disk clouds is 1.5x 10^{-2}. We can distinguish clearly between regions where the SiO or HCO+ are dominating.

Ams: We study the 12C/13C isotopic ratio in the disk of the central molecular zone and in the halo to trace gas accretion toward the Galactic center region in the Milky Way. Methods: Using the IRAM 30m telescope, we observe the J=1-0 rotational transition of HCO+, HCN, HNC and their 13C isotopic substitutions in order to measure the 12C/13C isotopic ratio. We observe 9 positions selected throughout the Galactic center region, including clouds at high latitude; locations where the X1 and X2 orbits associated with the barred potential are expected to intersect; and typical Galactic center molecular clouds. Results: We find a systematically higher 12C/13C isotopic ratio (>40) toward the halo and the X1 orbits than for the Galactic center molecular clouds (20-25). Our results point out to molecular gas which has undergone a different degree of nuclear processing than that observed in the gas towards the inner Galactic center region. Conclusions: The high isotopic ratios are consistent with the accretion of the gas from the halo and from the outskirts of the Galactic disk.

We study the disk-halo interaction, in the context of orbits and Giant Molecular loops (GMLs) in the Galactic center (GC) region. We present a large scale survey in the central kpc in SiO(J=2-1), HCO$^+$(J=1-0) and H$^{13}$CO$^+$ (J=1-0), observations in 3mm lines toward a region in two clumps M+5.3-0.3 and M-3.8+0.9 placed in the foot point of two molecular loops, and observations toward selected positions, to trace accretion of the gas from the halo