This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We present the first detections of CH3SH, C3H+, C3N, HCOOH, CH2CHCN, and H2CN in an extragalactic source. Namely the spiral arm of a galaxy located at z=0.89 on the line of sight to the radio-loud quasar PKS 1830-211. OCS, SO2, and NH2CN were also detected, raising the total number of molecular species identified in that early time galaxy to 54, not counting isotopologues. The detections were made in absorption against the SW quasar image, at 2 kpc from the galaxy centre, over the course of a Q band spectral line survey made with the Yebes 40 m telescope (rest-frame frequencies: 58.7-93.5 GHz). We derived the rotational temperatures and column densities of those species, which are found to be subthermally excited. The molecular abundances, and in particular the large abundances of C3H+ and of several previously reported cations, are characteristic of diffuse or translucent clouds with enhanced UV radiation or strong shocks.

The largest Key Science Program of the RadioAstron space VLBI mission is a survey of active galactic nuclei (AGN). The main goal of the survey is to measure and study the brightness of AGN cores in order to better understand the physics of their emission while taking interstellar scattering into consideration. In this paper we present detection statistics for observations on ground-space baselines of a complete sample of radio-strong AGN at the wavelengths of 18, 6, and 1.3 cm. Two-thirds of them are indeed detected by RadioAstron and are found to contain extremely compact, tens to hundreds of $\mu$as structures within their cores.

OH231.8+4.2 is an archetypal pre-planetary nebulae (pPN). It is a binary system surrounded by bipolar nebula. Some years ago the authors extensively studied it and performed several VLBI observations from which they obtained mas-resolution maps of the SiO (7 mm) and H2O (1.3 cm) maser emissions. H2O masers were found to be distributed in two areas along the symmetry axis of the nebulae oriented nearly north-south delineating a bipolar outflow and their astrometric positions were accurately measured. SiO masers, indicating the position of the Mira component of the binary system, form a structure perpendicular to the axis of the nebulae. The general picture of the source looked satisfactory, except for the relative position of the two masers. Surprisingly, SiO masers, were tentatively placed 250 mas away (370 AU) from the apparent center of the outflow. Using the ALMA we observed the SiO maser emission at 86 GHz and accurately derived the position of the Mira component. Combining our previous VLBA data and our new ALMA observations we are now able to give a more complete and detailed description of the inner part of this amazing pPN.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

We report results from a deep polarization imaging of the nearby radio galaxy 3C$\,$84 (NGC$\,$1275). The source was observed with the Global Millimeter VLBI Array (GMVA) at 86$\,$GHz at an ultra-high angular resolution of $50\mu$as (corresponding to 250$R_{s}$). We also add complementary multi-wavelength data from the Very Long Baseline Array (VLBA; 15 & 43$\,$GHz) and from the Atacama Large Millimeter/submillimeter Array (ALMA; 97.5, 233.0, and 343.5$\,$GHz). At 86$\,$GHz, we measure a fractional linear polarization of $\sim2$% in the VLBI core region. The polarization morphology suggests that the emission is associated with an underlying limb-brightened jet. The fractional linear polarization is lower at 43 and 15$\,$GHz ($\sim0.3-0.7$% and $<0.1$%, respectively). This suggests an increasing linear polarization degree towards shorter wavelengths on VLBI scales. We also obtain a large rotation measure (RM) of $\sim10^{5-6}~{\rm rad/m^{2}}$ in the core at $\gtrsim$43$\,$GHz. Moreover, the VLBA 43$\,$GHz observations show a variable RM in the VLBI core region during a small flare in 2015. Faraday depolarization and Faraday conversion in an inhomogeneous and mildly relativistic plasma could explain the observed linear polarization characteristics and the previously measured frequency dependence of the circular polarization. Our Faraday depolarization modeling suggests that the RM most likely originates from an external screen with a highly uniform RM distribution. To explain the large RM value, the uniform RM distribution, and the RM variability, we suggest that the Faraday rotation is caused by a boundary layer in a transversely stratified jet. Based on the RM and the synchrotron spectrum of the core, we provide an estimate for the magnetic field strength and the electron density of the jet plasma.

(abridged) Very long baseline interferometry (VLBI) observations at 86$\,$GHz (wavelength, $\lambda = 3\,$mm) reach a resolution of about 50 $\mu$as, probing the collimation and acceleration regions of relativistic outflows in active galactic nuclei. To extend the statistical studies of compact extragalactic jets, a large global 86 GHz VLBI survey of 162 radio sources was conducted in 2010-2011 using the Global Millimeter VLBI Array. The survey data attained a typical baseline sensitivity of 0.1 Jy and a typical image sensitivity of 5 mJy/beam, providing successful detections and images for all of the survey targets. For 138 objects, the survey provides the first ever VLBI images made at 86 GHz. Gaussian model fitting of the visibility data was applied to represent the structure of the sources. The Gaussian model-fit-based estimates of brightness temperature ($T_\mathrm{b}$) at the jet base (core) and in moving regions (jet components) downstream from the core were compared to the estimates of $T_\mathrm{b}$ limits made directly from the visibility data, demonstrating a good agreement between the two methods. The apparent brightness temperature estimates for the jet cores in our sample range from $2.5 \times 10^{9}\,$K to $ 1.3\times 10^{12}\,$K. A population model with a single intrinsic value of brightness temperature, $T_\mathrm{0}$, is applied to reproduce the observed $T_\mathrm{b}$ distribution. It yields $T_\mathrm{0} = (3.77^{+0.10}_{-0.14}) \times 10^{11}\,$K for the jet cores, implying that the inverse Compton losses dominate the emission. In the jet components, $T_\mathrm{0} =(1.42^{+0.16}_{-0.19})\times 10^{11}\,$K is found, slightly higher than the equipartition limit of $\sim5\times 10^{10}\,$K expected for these jet regions. For objects with sufficient structural detail detected, the adiabatic energy losses dominate the observed changes of $T_\mathrm{b}$ along the jet.

Context. Radio-loud active galactic nuclei in the early Universe are rare. The quasars J0906+6930 at redshift z=5.47 and J2102+6015 at z=4.57 stand out from the known sample with their compact emission on milliarcsecond (mas) angular scale with high (0.1-Jy level) flux densities measured at GHz radio frequencies. This makes them ideal targets for very long baseline interferometry (VLBI) observations. Aims. By means of VLBI imaging we can reveal the inner radio structure of quasars and model their brightness distribution to better understand the geometry of the jet and the physics of the sources. Methods. We present sensitive high-resolution VLBI images of J0906+6930 and J2102+6015 at two observing frequencies, 2.3 and 8.6 GHz. The data were taken in an astrometric observing programme involving a global five-element radio telescope array. We combined the data from five different epochs from 2017 February to August. Results. For one of the highest redshift blazars known, J0906+6930, we present the first-ever VLBI image obtained at a frequency below 8 GHz. Based on our images at 2.3 and 8.6 GHz, we confirm that this source has a sharply bent helical inner jet structure within ~3 mas from the core. The quasar J2102+6015 shows an elongated radio structure in the east-west direction within the innermost ~2 mas that can be described with a symmetric three-component brightness distribution model at 8.6 GHz. Because of their non-pointlike mas-scale structure, these sources are not ideal as astrometric reference objects. Our results demonstrate that VLBI observing programmes conducted primarily with astrometric or geodetic goals can be utilized for astrophysical purposes as well.

M87 is one of the nearest radio galaxies with a prominent jet extending from sub-pc to kpc-scales. Because of its proximity and large mass of the central black hole, it is one of the best radio sources to study jet formation. We aim at studying the physical conditions near the jet base at projected separations from the BH of $\sim7-100$ Schwarzschild radii ($R_{\rm sch}$). Global mm-VLBI Array (GMVA) observations at 86 GHz ($\lambda=3.5\,$mm) provide an angular resolution of $\sim50\mu$as, which corresponds to a spatial resolution of only $7~R_{\rm sch}$ and reach the small spatial scale. We use five GMVA data sets of M87 obtained during 2004--2015 and present new high angular resolution VLBI maps at 86GHz. In particular, we focus on the analysis of the brightness temperature, the jet ridge lines, and the jet to counter-jet ratio. The imaging reveals a parabolically expanding limb-brightened jet which emanates from a resolved VLBI core of $\sim(8-13) R_{\rm sch}$ size. The observed brightness temperature of the core at any epoch is $\sim(1-3)\times10^{10}\,$K, which is below the equipartition brightness temperature and suggests magnetic energy dominance at the jet base. We estimate the diameter of the jet at its base to be $\sim5 R_{\rm sch}$ assuming a self-similar jet structure. This suggests that the sheath of the jet may be anchored in the very inner portion of the accretion disk. The image stacking reveals faint emission at the center of the edge-brightened jet on sub-pc scales. We discuss its physical implication within the context of the spine-sheath structure of the jet.

The high brightness temperatures, $T_\mathrm{b}\gtrsim 10^{13}$ K, detected in several active galactic nuclei by RadioAstron space VLBI observations challenge theoretical limits. Refractive scattering by the interstellar medium may affect such measurements. We quantify the scattering properties and the sub-mas scale source parameters for the quasar B0529+483. Using RadioAstron correlated flux density measurements at 1.7, 4.8, and 22 GHz on projected baselines up to 240,000 km we find two characteristic angular scales in the quasar core, about 100 $\mu$as and 10 $\mu$as. Some indications of scattering substructure are found. Very high brightness temperatures, $T_\mathrm{b}\geq 10^{13}$ K, are estimated at 4.8 GHz and 22 GHz even taking into account the refractive scattering. Our findings suggest a clear dominance of the particle energy density over the magnetic field energy density in the core of this quasar.

We present a proof of concept on the coupling of radio astronomical receivers and spectrometers with chemical reactorsand the performances of the resulting setup for spectroscopy and chemical simulations in laboratory astrophysics. Several experiments including cold plasma generation and UV photochemistry were performed in a 40\,cm long gas cell placed in the beam path of the Aries 40\,m radio telescope receivers operating in the 41-49 GHz frequency range interfaced with fast Fourier transform spectrometers providing 2 GHz bandwidth and 38 kHz resolution. The impedance matching of the cell windows has been studied using different materials. The choice of the material and its thickness was critical to obtain a sensitivity identical to that of standard radio astronomical observations. Spectroscopic signals arising from very low partial pressures of CH3OH, CH3CH2OH, HCOOH, OCS,CS, SO2 (<1E-03 mbar) were detected in a few seconds. Fast data acquisition was achieved allowing for kinetic measurements in fragmentation experiments using electron impact or UV irradiation. Time evolution of chemical reactions involving OCS, O2 and CS2 was also observed demonstrating that reactive species, such as CS, can be maintained with high abundance in the gas phase during these experiments.

We present an approach to testing the gravitational redshift effect using the RadioAstron satellite. The experiment is based on a modification of the Gravity Probe A scheme of nonrelativistic Doppler compensation and benefits from the highly eccentric orbit and ultra-stable atomic hydrogen maser frequency standard of the RadioAstron satellite. Using the presented techniques we expect to reach an accuracy of the gravitational redshift test of order $10^{-5}$, a magnitude better than that of Gravity Probe A. Data processing is ongoing, our preliminary results agree with the validity of the Einstein Equivalence Principle.

In September 2016, the microquasar Cygnus X-3 underwent a giant radio flare, which was monitored for 6 days with the Medicina Radio Astronomical Station and the Sardinia Radio Telescope. Long observations were performed in order to follow the evolution of the flare on a hourly scale, covering six frequency ranges from 1.5 GHz to 25.6 GHz. The radio emission reached a maximum of 13.2 +/- 0.7 Jy at 7.2 GHz and 10 +/- 1 Jy at 18.6 GHz. Rapid flux variations were observed at high radio frequencies at the peak of the flare, together with rapid evolution of the spectral index: alpha steepened from 0.3 to 0.6 within 5 hours. This is the first time that such fast variations are observed, giving support to the evolution from optically thick to optically thin plasmons in expansion moving outward from the core. Based on the Italian network (Noto, Medicina and SRT) and extended to the European antennas (Torun, Yebes, Onsala), VLBI observations were triggered at 22 GHz on five different occasions, four times prior to the giant flare, and once during its decay phase. Flux variations of 2-hour duration were recorded during the first session. They correspond to a mini-flare that occurred close to the core ten days before the onset of the giant flare. From the latest VLBI observation we infer that four days after the flare peak the jet emission was extended over 30 mas.

The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler measurements of planetary spacecraft using the Mars Express spacecraft. The technique will be used in the framework of Planetary Radio Interferometry and Doppler Experiments in various planetary missions, in particular in fly-by mode. We advanced a novel approach to spacecraft data processing using the techniques of Doppler and phase-referenced very long baseline interferometry spacecraft tracking. We achieved, on average, mHz precision (30 {\mu}m/s at a 10 seconds integration time) for radial three-way Doppler estimates and sub-nanoradian precision for lateral position measurements, which in a linear measure (at a distance of 1.4 AU) corresponds to ~50 m.

Observations of 28SiO v=0 J=1-0 line emission (7-mm wavelength) from AGB stars show in some cases peculiar profiles, composed of a central intense component plus a wider plateau. Very similar profiles have been observed in CO lines from some AGB stars and most post-AGB nebulae and, in these cases, they are clearly associated with the presence of conspicuous axial symmetry and bipolar dynamics. We present systematic observations of 28SiO v=0 J=1-0 emission in 28 evolved stars, performed with the 40~m radio telescope of the IGN in Yebes, Spain. We find that the composite core plus plateau profiles are almost always present in O-rich Miras, OH/IR stars, and red supergiants. They are also found in one S-type Mira ($\chi$ Cyg), as well as in two semiregular variables (X Her and RS Cnc) that are known to show axial symmetry. In the other objects, the profiles are simpler and similar to those of other molecular lines. The composite structure appears in the objects in which SiO emission is thought to come from the very inner circumstellar layers, prior to dust formation. The central spectral feature is found to be systematically composed of a number of narrow spikes, except for X Her and RS Cnc, in which it shows a smooth shape that is very similar to that observed in CO emission. These spikes show a significant (and mostly chaotic) time variation, while in all cases the smooth components remain constant within the uncertainties. The profile shape could come from the superposition of standard wide profiles and a group of weak maser spikes. Alternatively, we speculate that the very similar profiles detected in objects that are axisymmetric may be indicative of the systematic presence of a significant axial symmetry in the very inner circumstellar shells around AGB stars; the presence of such symmetry would be independent of the probable weak maser effects in the central spikes.

Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH+) has been identified through the J=5-4 and J=10-9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 30m telescope. We derive beam averaged column densities for NCCNH+ of (8.6+/-4.4)e10 cm-2 in TMC-1 and (3.9+/-1.8)e10 cm-2 in L483, which translate to fairly low fractional abundances relative to H2, in the range (1-10)e-12. The chemistry of protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio NCCNH+/NCCN of 1e-4, which implies that the abundance of cyanogen in dark clouds could be as high as (1-10)e-8 relative to H2, i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC3N.

The Global mm-VLBI Array (GMVA) is a network of 14 3\,mm and 7\,mm capable telescopes spanning Europe and the United States, with planned extensions to Asia. The array is capable of sensitive maps with angular resolution often exceeding 50\,$\mu$as. Using the GMVA, a large sample of prominent $\gamma$-ray blazars have been observed approximately 6 monthly from later 2008 until now. Combining 3\,mm maps from the GMVA with near-in-time 7\,mm maps from the VLBA-BU-BLAZAR program and 2\,cm maps from the MOJAVE program, we determine the sub-pc morphology and high frequency spectral structure of $\gamma$-ray blazars. The magnetic field strength can be estimated at different locations along the jet under the assumption of equipartition between magnetic field and relativistic particle energies. Making assumptions on the jet magnetic field configuration (e.g. poloidal or toroidal), we can estimate the separation of the mm-wave "core" and the jet base, and estimate the strength of the magnetic field there. The results of this analysis show that on average, the magnetic field strength decreases with a power-law $B \propto r^{-n}$, $n=0.3 \pm 0.2$. This suggests that on average, the mm-wave "core" is $\sim 1-3$\,pc downstream of the de-projected jet apex and that the magnetic field strength is of the order $B_{\rm{apex}} \sim 5-20$\,kG, broadly consistent with the predictions of magnetic jet launching (e.g. via magnetically arrested disks (MAD)).

The gas cloud G2 falling toward Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, is supposed to provide valuable information on the physics of accretion flows and the environment of the black hole. We observed Sgr A* with four European stations of the Global Millimeter Very Long Baseline Interferometry Array (GMVA) at 86 GHz on 1 October 2013 when parts of G2 had already passed the pericenter. We searched for possible transient asymmetric structure -- such as jets or winds from hot accretion flows -- around Sgr A* caused by accretion of material from G2. The interferometric closure phases remained zero within errors during the observation time. We thus conclude that Sgr A* did not show significant asymmetric (in the observer frame) outflows in late 2013. Using simulations, we constrain the size of the outflows that we could have missed to ~2.5 mas along the major axis, ~0.4 mas along the minor axis of the beam, corresponding to approximately 232 and 35 Schwarzschild radii, respectively; we thus probe spatial scales on which the jets of radio galaxies are suspected to convert magnetic into kinetic energy. As probably less than 0.2 Jy of the flux from Sgr A* can be attributed to accretion from G2, one finds an effective accretion rate eta*Mdot < 1.5*10^9 kg/s ~ 7.7*10^-9 Mearth/yr for material from G2. Exploiting the kinetic jet power--accretion power relation of radio galaxies, one finds that the rate of accretion of matter that ends up in jets is limited to Mdot < 10^17 kg/s ~ 0.5 Mearth/yr, less than about 20% of the mass of G2. Accordingly, G2 appears to be largely stable against loss of angular momentum and subsequent (partial) accretion at least on time scales < 1 year.

The Global mm-VLBI Array (GMVA) is the highest angular resolution imaging interferometer currently available as a common user facility. It is capable of angular resolutions on the order of 40 microarcseconds. Currently 14 stations in the United States and Europe participate in global 3 mm VLBI observations. The GMVA is used for continuum and spectroscopic imaging, probing the central regions of active galaxies and the origin of jets as these regions are typically not observable at longer wavelengths due to synchrotron self-absorption. In early 2012, fringes were detected to the three stations of the Korean VLBI Network (KVN), opening the possibility of extending the baseline coverage of the VLBI array to the East. In these proceedings, we will present recent images from a monitoring program of gamma-ray blazars using the GMVA, including the sources 3C454.3 and 0235+164, and an update of its current status and abilities.

Aims. The phase scintillation of the European Space Agency's (ESA) Venus Express (VEX) spacecraft telemetry signal was observed at X-band (\lambda = 3.6 cm) with a number of radio telescopes of the European VLBI Network (EVN) in the period 2009-2013. Methods. We found a phase fluctuation spectrum along the Venus orbit with a nearly constant spectral index of -2.42 +/-0.25 over the full range of solar elongation angles from 0{\deg} to 45{\deg}, which is consistent with Kolmogorov turbulence. Radio astronomical observations of spacecraft signals within the solar system give a unique opportunity to study the temporal behaviour of the signal's phase fluctuations caused by its propagation through the interplanetary plasma and the Earth's ionosphere. This gives complementary data to the classical interplanetary scintillation (IPS) study based on observations of the flux variability of distant natural radio sources. Results. We present here our technique and the results on IPS. We compare these with the total electron content (TEC) for the line of sight through the solar wind. Finally, we evaluate the applicability of the presented technique to phase-referencing Very Long Baseline Interferometry (VLBI) and Doppler observations of currently operational and prospective space missions.

We study the origin of large abundances of complex organic molecules in the Galactic center (GC). We carried out a systematic study of the complex organic molecules CH3OH, C2H5OH, (CH3)2O, HCOOCH3, HCOOH, CH3COOH, H2CO, and CS toward 40 GC molecular clouds. Using the LTE approximation, we derived the physical properties of GC molecular clouds and the abundances of the complex molecules.The CH3OH abundance between clouds varies by nearly two orders of magnitude from 2.4x10^{-8} to 1.1x10^{-6}. The abundance of the other complex organic molecules relative to that of CH3OH is basically independent of the CH3OH abundance, with variations of only a factor 4-8. The abundances of complex organic molecules in the GC are compared with those measured in hot cores and hot corinos, in which these complex molecules are also abundant. We find that both the abundance and the abundance ratios of the complex molecules relative to CH3OH in hot cores are similar to those found in the GC clouds. However, hot corinos show different abundance ratios than observed in hot cores and in GC clouds. The rather constant abundance of all the complex molecules relative to CH3OH suggests that all complex molecules are ejected from grain mantles by shocks. Frequent (similar 10^{5}years) shocks with velocities >6km/s are required to explain the high abundances in gas phase of complex organic molecules in the GC molecular clouds. The rather uniform abundance ratios in the GC clouds and in Galactic hot cores indicate a similar average composition of grain mantles in both kinds of regions. The Sickle and the Thermal Radio Arches, affected by UV radiation, show different relative abundances in the complex organic molecules due to the differentially photodissociation of these molecules.

The nucleus of M82 has been mapped in several 3mm and 1mm lines of CN, HCN, C2H, c-C3H2, CH3C2H, HC3N and HOC+ using the IRAM 30m telescope. These species have been purposely selected as good tracers of photon-dominated chemistry. We have measured [CN]/[HCN] ~ 5 in the inner 650 pc galaxy disk. Furthermore, we have detected the HOC+ 1--0 line with an intensity similar to that of the H13CO+ 1--0 line. This implies a [HCO+]/[HOC+] ratio of ~40. These results corroborate the existence of a giant photo-dissociation region (PDR) in the nucleus of M82. In fact, the low [HCO+]/[HOC+] ratio can only be explained if the nucleus of M82 is formed by small (r<0.02-0.2 pc) and dense (n ~ a few 10^4--10^5 cm^{-3}) clouds immersed in an intense UV field (G_0 ~ 10^4 in units of the Habing field). The detection of the hydrocarbons c-C3H2 and CH3C2H in the nucleus of M82 suggests that a complex carbon chemistry is developing in this giant PDR.

We present images of the J=10-9 rotational lines of HC3N in the vibrationally excited levels 1v7, 1v6 and 1v5 of the hot core (HC) in Orion KL. The images show that the spatial distribution and the size emission from the 1v7 and 1v5 levels are different. While the J=10-9 1v7 line has a size of 4''x 6'' and peaks 1.1'' NE of the 3 mm continuum peak, the J=10--9 1v5 line emission is unresolved (<3'') and peaks 1.3'' south of the 3 mm peak. This is a clear indication that the HC is composed of condensations with very different temperatures (170 K for the 1v7 peak and $>230$ K for the 1v5 peak). The temperature derived from the 1v7 and 1v5 lines increases with the projected distance to the suspected main heating source I. Projection effects along the line of sight could explain the temperature gradient as produced by source I. However, the large luminosity required for source I, >5 10^5 Lsolar, to explain the 1v5 line suggests that external heating by this source may not dominate the heating of the HC. Simple model calculations of the vibrationally excited emission indicate that the HC can be internally heated by a source with a luminosity of 10^5 Lsolar, located 1.2'' SW of the 1v5 line peak (1.8'' south of source I). We also report the first detection of high-velocity gas from vibrationally excited HC3N emission. Based on excitation arguments we conclude that the main heating source is also driving the molecular outflow. We speculate that all the data presented in this letter and the IR images are consistent with a young massive protostar embedded in an edge-on disk.

We present Infrared Space Observatory (ISO) observations of fine structure lines toward the Radio Arc region, close to Sgr A, in the Galactic center region (GCR). The ionization at large scale is dominated by the clusters of massive stars in this region: the Quintuplet and the Arches cluster. Any other possible mechanism as interaction with the magnetic fields revealed by the Non-Thermal Filaments (NTFs) must play a minor role. The radiation in this area is relatively hard ($T_\mathrm{eff}\sim 33-35$ kK) but diluted. It reaches large distances due to the inhomogeneity of the medium.

We present some recent studies on the physical conditions and the chemistry of the interstellar medium (ISM) in the Galactic center region (GCR). Large scale shocks in the context of a bar potential can play a role in explaining the observed properties, at least in some particular regions. Alternatively, we propose that the observed properties (inhomogeneity, high temperatures of the molecular gas, rich chemistry,...) of the GCR ISM can be the consequence of the recent past ($10^6-10^7$ years ago) of massive star formation in the GCR.

(Short version) We present an analysis of selected fine structure lines observed ISO toward the Radio Arc in the Galactic center region (GCR). We study the density of the ionized gas and the large scale ionization structure in a region of ~30x30 pc^2 by means of the NeIII 15.6 mum/NeII 12.8 mum and NIII 57 mum/NII 122 mum line ratios. The electron densities (n_e) derived from the OIII lines ratio indicate the presence of diffuse ionized material with n_e of \~10^{1.8-2.6} cm^{-3}. The NIII/NII ratio show two clear gradients, one pointing toward the Quintuplet cluster and the other pointing toward the Arches cluster. The model used to explain the ionization structure shows that the large scale ionization of the whole region can be accounted for by the UV radiation produced by the Quintuplet and the Arches cluster. Any other ionization mechanism should play a minor role. We also investigate the influence of the clusters on the bubble of warm dust seen in the MSX infrared images. The warm dust is well correlated with the ionized gas indicating that the dust is also heated by the radiation from both clusters. Furthermore, the elliptical rather than circular symmetry of some structures like the Thermal Filaments can also be explained when one considers the combined effects of both the Arches and the Quintuplet clusters. The bubble is filled with continuum emission of hard X rays and with emision from the 6.4 keV line of neutral or low ionized Fe.

We present observations of C2H5OH toward molecular clouds in Sgr A, Sgr B2 and associated with thermal and non-thermal features in the Galactic center. C2H5OH emission in Sgr A and Sgr B2 is widespread, but not uniform. C2H5OH emission is much weaker or it is not detected in some molecular clouds in both complexes, in particular those with radial velocities between 70 and 120 km/s. While most of the clouds associated with the thermal features do not show C2H5OH emission, that associated with the Non-Thermal Radio Arc shows emission. The fractional abundance of C2H5OH in most of the clouds with radial velocities between 0 and 70 km/s in Sgr A and Sgr B2 is relatively high, of few 1e-8. The C2H5OH abundance decreases by more than one order of magnitude (aprox. 1e-9) in the clouds associated with the thermal features. The large abundance of C2H5OH in the gas-phase indicates that C2H5OH has formed in grains and released to gas-phase by shocks in the last aprox. 1e5 years.

We present ISO observations of several H2 pure-rotational lines (from S(0) to S(5)) towards a sample of 16 molecular clouds distributed along the central ~ 500 pc of the Galaxy. We also present C18O and 13CO J=1->0 and J=2->1 observations of these sources made with the IRAM-30m telescope. With the CO data we derive H2 densities of 10e(3.5-4.0) cm-3 and H2 column densities of a few 10e22 cm-2. We have corrected the H2 data for ~ 30 magnitudes of visual extinction using a self-consistent method. In every source, we find that the H2 emission exhibits a large temperature gradient. The S(0) and S(1) lines trace temperatures (T) of ~150 K while the S(4) and S(5) lines indicate temperatures of ~ 600K. The warm H2 column density is typically ~1-2 x 10e22 cm-2, and is predominantly gas with T=150 K. This is the first direct estimate of the total column density of the warm molecular gas in the Galactic center region. These warm H2 column densities represent a fraction of ~ 30 % of the gas traced by the CO isotopes emission. The cooling by H2 in the warm component is comparable to that by CO. Comparing our H2 and CO data with available ammonia NH3 observations from literature one obtains relatively high NH3 abundances of a few 10e(-7) in both the warm and the cold gas. A single shock or Photo-Dissociation Region (PDR) cannot explain all the observed H2 lines. Alternatives for the heating mechanisms are discussed.

We present single dish and interferometric maps of several rotational transitions of HC3N vibrationally excited levels towards Sgr B2. The single dish maps show, for the first time, new four bright clumps (Sgr B2R1 to B2R4) in the ridge connecting the main cores Sgr B2M and Sgr B2N, and a low brightness extended region to the west of the ridge (Sgr B2W). The physical properties and the kinematics of all hot cores have been derived from the HC3N* lines. Our high angular resolution images show that the Sgr B2N hot core breaks in two different hot cores, Sgr B2N1 and N2, with different radial velocities and separated by 2 arcsec in declination. We find that the excitation of the HC3N* emission in all hot cores can be represented by a single temperature and that the linewidth of the HC3N* rotational lines arising from different vibrational levels systematically decreases as the energy of the vibrational level increases. The systematic trend in the linewidth is likely related to the increase of the velocity as the distance to the exciting source increases. The different HC3N abundances in the hot cores reflect different stages of evolution due to time dependent chemistry and/or photo-dissociation by UV radiation from nearby HII regions. We propose that Sgr B2M and B2N contain a cluster of 20-30 hot cores and that the Sgr B2R1-B2R4 hot cores represent isolated formation of massive stars.

In this letter we present a map of the GC in the 1-0 line of SiO covering the same region than the ASCA observations of the Fe 6.4 keV line. The SiO emission is restricted to molecular clouds with radial velocity between 10 and 60 km/s in the Sgr A and Sgr B complexes. We find a correlation between the SiO morphology and the spatial distribution of the Fe 6.4 keV line, on the large scale and also within Sgr A and Sgr B. The SiO abundance increases by a factor of 20 in the regions with strong Fe 6.4 keV line. This indicates that the Fe 6.4 keV line mainly arises from molecular clouds with large gas phase abundance of refractory elements. We discuss the implications of the correlation on the origin of the hard X-rays, and the heating and the chemistry of the molecular clouds in the GC

We present ISO observations of the S(0), S(1), S(2), and S(3) rotational lines of molecular hydrogen from two molecular clouds near the Galactic Center (GC). We have also measured continuum dust emission at infrared wavelengths with ISO and the rotational radio lines J=1-0 of 13CO and C18O and J=2-1 of C18O with the IRAM-30m telescope. Using the dust continuum spectra and the CO lines we derive a total visual extinction of 15-20 magnitudes toward these GC clouds. After correcting the H2 data for extinction, the gas temperatures are \~250 K and the column densities of warm gas are ~2e21 cm-2. This is the first direct measure of the H2 column densities of the warm component; with this, we estimate an NH3 abundance in the warm gas of ~ 2e-7. The column density of warm gas is, at least, a factor of 100 larger than the corresponding column densities derived from the warm dust. The observed ortho-to-para ratio (OTPR) is ~ 1, clearly below the local thermodynamical equilibrium (LTE) OTPR for gas at 250 K of ~3. Low velocity shocks (~ 10 km/s) are the most likely explanation for the column densities of warm gas and dust and the non-LTE H2 OTPR.

We have observed the S(0), S(1), S(2), S(3), S(4) and S(5) rotational lines of molecular hydrogen (H2) towards the peak of the photodissociation region (PDR) associated with the reflection nebula NGC 7023. The observed H2 line ratios show that they arise in warm gas with kinetic temperatures ~300 - 700 K. However, the data cannot be fitted by an ortho- to para- (OTP) ratio of 3. An OTP ratio in the range ~1.5 - 2 is necessary to explain our observations. This is the first detection of a non-equilibrium OTP ratio measured from the H2 pure-rotational lines in a PDR. The existence of a dynamical PDR is discussed as the most likely explanation for this low OTP ratio.

We have mapped the J=1-0 line of SiO in a 1 degree x 12 arcmin. (lxb) region around the Galactic center (GC) with an angular resolution of 2 arcmin. (approx. 4 pc). In contrast to the spatial distribution of other high dipole moment molecules like CS, whose emission is nearly uniform, the SiO emission is very fragmented and it is only associated with some molecular clouds. In particular, it is remarkable that the SiO emission closely follows the non-thermal radio arc in the GC. The SiO clouds are more extended than the beam with typical sizes between 4 and 20 pc. High angular resolution (26 arcsec.) mapping in the J=1-0 line of SiO toward the molecular clouds in Sgr B2 and Sgr A shows that the SiO emission is relatively smooth with structures of typically 2 pc. From the line intensities of the J=1-0, J=3-2 and J=5-4 transitions of SiO we derive H2 densities for these clouds of a few 1e4 cm-3. The SiO fractional abundances are approx. 1e-9 for the SiO clouds and 1e-10 for the other molecular clouds in the GC. The characteristics (size and H2 densities) of the SiO emission in the GC are completely different from those observed in the Galactic disk, where the SiO emission arises from much smaller regions with larger H2 densities. We briefly discuss the implications of the SiO emission in the molecular clouds of the GC. We conclude that the particular chemistry in these clouds is probably related to large scale fast shocks occurring in the Galactic center region.