Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The compact, non-thermal emission in DoAr21 has been studied with the VLBA to investigate the possibility that the residuals of the astrometry fitting are due to the reflex motion induced by a possible companion. We find that the fitting of VLBA astrometric observations of DoAr21 improves significantly by adding the orbital motions of three companions. We obtain an improved distance to the source of $134.6\pm1.0$ pc, and estimate that the central star, DoAr21, has a mass of about $2.04\pm0.70~M_{\odot}$. We suggest that DoAr21 represents a unique case where two sub-stellar companions, DoAr21b and DoAr21c ($m_{b}\sim 35.6\pm27.2~M_{jup}$ and $m_{c}\sim44.0\pm13.6~M_{jup}$, respectively), have been found associated to a relatively low mass, pre-main sequence star. In addition, we find that this WTTau star is an astrometric double system, having a low-mass star companion, DoAr21$B$ ($m_{B}\sim0.35\pm0.12~M_{\odot}$), in a relatively eccentric orbit. The orbit of this low-mass stellar companion is compact, while the Brown Dwarfs are located in external orbits. DoAr21$c$ has the strongest astrometric signature in the periodogram, while DoAr21$B$ has a weak but significant signature. On the other hand, the astrometric signature of DoAr21$b$ does not appear in the periodogram, however, this Brown Dwarf was directly detected in some of the VLBA observations. The estimated orbital periods of DoAr21$B$, DoAr21$b$ and DoAr21$c$ are $P_{B}\sim92.92\pm0.02$, $P_{b}\sim450.9\pm3.8$ and $P_{c}\sim1013.5\pm25.3$ days, respectively. Since the estimated age of this young star is about 0.4$-$0.8 Myrs, the detected Brown Dwarf companion is among the youngest companions observed to date.

Aims:We take advantage of the second data release of the Gaia space mission and the state-of-the-art astrometry delivered from very long baseline interferometry observations to revisit the structure and kinematics of the nearby Taurus star-forming region. Methods: We apply a hierarchical clustering algorithm for partitioning the stars in our sample into groups (i.e., clusters) that are associated with the various molecular clouds of the complex, and derive the distance and spatial velocity of individual stars and their corresponding molecular clouds. Results: We show that the molecular clouds are located at different distances and confirm the existence of important depth effects in this region reported in previous studies. For example, we find that the L 1495 molecular cloud is located at $d=129.9^{+0.4}_{-0.3}$ pc, while the filamentary structure connected to it (in the plane of the sky) is at $d=160.0^{+1.2}_{-1.2}$ pc. We report B 215 and L 1558 as the closest ($d=128.5^{+1.6}_{-1.6}$ pc) and most remote ($d=198.1^{+2.5}_{-2.5}$ pc) substructures of the complex, respectively. The median inter-cloud distance is 25 pc and the relative motion of the subgroups is on the order of a few km/s. We find no clear evidence for expansion (or contraction) of the Taurus complex, but signs of the potential effects of a global rotation. Finally, we compare the radial velocity of the stars with the velocity of the underlying $^{13}$CO molecular gas and report a mean difference of $0.04\pm0.12$ km/s (with r.m.s. of 0.63 km/s) confirming that the stars and the gas are tightly coupled.

The GLOSTAR survey will study the star formation in the Galactic plane between $-2^{\circ}< \textit{l}< 85^{\circ}$ and $|b| < 1^{\circ}$ with unprecedented sensitivity in both, flux density ($\sim$ 40 $\mu Jy$ beam$^{-1}$) and range of angular scales ($\sim$ 1."5 to the largest radio structures in the Galaxy). In this paper we present the first results obtained from a radio continuum map of a 16 square degree sized region of the Galactic plane centered on $\textit{l} = 32^{\circ}$ and $b = 0^{\circ}$ ($28^{\circ} < \textit{l} < 36^{\circ}$ and $|b| < 1^{\circ}$). This map has a resolution of 18" and sensitivity of $\sim$ 60-150 $\mu Jy$ beam$^{-1}$. We present data acquired with the VLA in D-configuration. Two 1 GHz wide sub-bands were observed simultaneously and centred at 4.7 and 6.9 GHz. These data were calibrated and imaged using the $\textit{Obit}$ software package. The source extraction has been performed using the BLOBCAT software package and verified through a combination of visual inspection and cross-matching with other radio and mid-infrared surveys. The final catalog consists of 1575 discrete radio sources and 27 large scale structures (including W43 and W44). By cross-matching with other catalogs and calculating the spectral indices ($S(\nu) \propto \nu^\alpha$), we have classified 231 continuum sources as HII regions, 37 as ionization fronts, and 46 as planetary nebulae. The longitude and latitude distribution and negative spectral indices are all consistent with the vast majority of the unclassified sources being extragalactic background sources. We present a catalog of 1575 radio continuum sources and discuss their physical properties, emission nature and relation with previously reported. These first GLOSTAR results have increased the number of reliable HII regions in this part of the Galaxy by a factor of four.

We present Gaia-DR2 astrometry of a sample of YSO candidates in Ophiuchus, Serpens Main and Serpens South/W40 in the Aquila Rift, which had been mainly identified by their infrared excess with Spitzer. We compare the Gaia-DR2 parallaxes against published and new parallaxes obtained from our Very Long Baseline Array (VLBA) program GOBELINS. We obtain consistent results between Gaia and the VLBA for the mean parallaxes in each of the regions analyzed here. We see small offsets, when comparing mean values, of a few tens of micro-arcseconds in the parallaxes, which are either introduced by the Gaia zero-point error or due to a selection effect by Gaia toward the brightest, less obscured stars. Gaia-DR2 data alone conclusively places Serpens Main and Serpens South at the same distance, as we first inferred from VLBA data alone in a previous publication. Thus, Serpens Main, Serpens South and W40 are all part of the same complex of molecular clouds, located at a mean distance of 436+/-9 pc. In Ophiuchus, both Gaia and VLBA suggest a small parallax gradient across the cloud, and the distance changes from 144.2+/-1.3 pc to 138.4+/-2.6 pc when going from L1689 to L1688.

The compact radio source Sagittarius~A$^*$ (Sgr~A$^*$)in the Galactic Center is the primary supermassive black hole candidate. General relativistic magnetohydrodynamical (GRMHD) simulations of the accretion flow around Sgr\,A$^*$ predict the presence of sub-structure at observing wavelengths of $\sim 3$\,mm and below (frequencies of 86\,GHz and above). For very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at this frequency the blurring effect of interstellar scattering becomes subdominant, and arrays such as the High Sensitivity Array (HSA) and the global mm-VLBI Array (GMVA) are now capable of resolving potential sub-structure in the source. Such investigations improve our understanding of the emission geometry of the mm-wave emission of Sgr\,A$^*$, which is crucial for constraining theoretical models and for providing a background to interpret 1\,mm VLBI data from the Event Horizon Telescope (EHT). We performed high-sensitivity very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at 3\,mm using the Very Long Baseline Array (VLBA) and the Large Millimeter Telescope (LMT) in Mexico on two consecutive days in May 2015, with the second epoch including the Green Bank Telescope (GBT). We find an overall source geometry that matches previous findings very closely, showing a deviation in fitted model parameters less than 3\% over a time scale of weeks and suggesting a highly stable global source geometry over time. The reported sub-structure in the 3\,mm emission of Sgr\,A$^*$ is consistent with theoretical expectations of refractive noise on long baselines. However, comparing our findings with recent results from 1\,mm and 7\,mm VLBI observations, which also show evidence for east-west asymmetry, an intrinsic origin cannot be excluded. Confirmation of persistent intrinsic substructure will require further VLBI observations spread out over multiple epochs.

We present an analysis of the astrometric results from Gaia second data release (DR2) to Young Stellar Objects (YSOs) in star-forming regions related to the Gould Belt. These regions are Barnard 59, Lupus 1 to 4, Chamaeleon I and II, $\epsilon$-Chamaeleontis, the Cepheus flare, IC 5146 and Corona Australis. The mean distance to the YSOs in each region are consistent with earlier estimations, though a significant improvement to the final errors was obtained. The mean distances to the star-forming regions were used to fit an ellipsoid of size $(358\pm7)\times(316\pm13)\times(70\pm4)$ pc, and centered at $(X_0,Y_0,Z_0)=(-82\pm15, 39\pm7, -25\pm4)$ pc, consistent with recently determined parameter of the Gould Belt. The mean proper motions were combined with radial velocities from the literature to obtain the three dimensional motion of the star-forming regions, which are consistent with a general expansion of the Gould Belt. We estimate that this expansion is occurring at a velocity of $2.5\pm0.1$ km s$^{-1}$. This is the first time that YSOs motions are used to investigate the kinematic of the Gould Belt. As an interesting side result, we also identified stars with large peculiar velocities.

We derive the distance and structure of the Perseus molecular cloud by combining trigonometric parallaxes from Very Long Baseline Array (VLBA) observations, taken as part of the GOBELINS survey, and Gaia Data Release 2. Based on our VLBA astrometry, we obtain a distance of 321+/-10 pc for IC 348. This is fully consistent with the mean distance of 320+/-26 measured by Gaia. The VLBA observations toward NGC 1333 are insufficient to claim a successful distance measurement to this cluster. Gaia parallaxes, on the other hand, yield a mean distance of 293+/-22 pc. Hence, the distance along the line of sight between the eastern and western edges of the cloud is ~30 pc, which is significantly smaller than previously inferred. We use Gaia proper motions and published radial velocities to derive the spatial velocities of a selected sample of stars. The average velocity vectors with respect to the LSR are (u,v,w) = (-6.1+/-1.6, 6.8+/-1.1, -0.9+/-1.2) and (-6.4+/-1.0, 2.1+/-1.4, -2.4+/-1.0) km/s for IC 348 and NGC 1333, respectively. Finally, our analysis of the kinematics of the stars has shown that there is no clear evidence of expansion, contraction, or rotational motions within the clusters.

We present observations of VLBA 20, a radio source found towards the edge of the Orion Nebula Cluster (ONC). Nonthermal emission dominates the spectral energy distribution of this object from radio to mid-infrared regime, suggesting that VLBA 20 is extragalactic. This source is heavily scattered in the radio regime. Very Long Baseline Array observations resolve it to ~34x19 mas at 5 GHz, and the wavelength dependence of the scattering disk is consistent with $\nu^{-2}$ at other frequencies. The origin of the scattering is most likely the ionized X-ray emitting gas from the winds of the most massive stars of the ONC. The scattering is highly anisotropic, with the axis ratio of 2:1, higher than what is typically observed towards other sources.

IRAS 16293-2422 is a very well studied young stellar system seen in projection towards the L1689N cloud in the Ophiuchus complex. However, its distance is still uncertain with a range of values from 120 pc to 180 pc. Our goal is to measure the trigonometric parallax of this young star by means of H$_2$O maser emission. We use archival data from 15 epochs of VLBA observations of the 22.2 GHz water maser line. By modeling the displacement on the sky of the H$_2$O maser spots, we derived a trigonometric parallax of $7.1\pm1.3$ mas, corresponding to a distance of $141_{-21}^{+30}$ pc. This new distance is in good agreement with recent values obtained for other magnetically active young stars in the L1689 cloud. We relate the kinematics of these masers with the outflows and the recent ejections powered by source A in the system.

The LkH$\alpha$ 101 cluster takes its name from its more massive member, the LkH$\alpha$ 101star, which is an $\sim15$ M$_\odot$ star whose true nature is still unknown. The distance to the LkH$\alpha$ 101 cluster has been controversial for the last few decades, with estimated values ranging from 160 to 800 pc. We have observed members and candidate members of the LkH$\alpha$ 101 cluster with signs of magnetic activity, using the Very Long Baseline Array, in order to measure their trigonometric parallax and, thus, obtain a direct measurement of their distances. A young star member, LkH$\alpha$ 101 VLA J043001.15+351724.6, was detected at four epochs as a single radio source. The best fit to its displacement on the plane of the sky yields a distance of 535$\pm$29 pc. We argue that this is the distance to the LkH$\alpha$ 101 cluster.

We present the first results of the Gould's Belt Distances Survey (GOBELINS), a project aimed at measuring the proper motion and trigonometric parallax of a large sample of young stars in nearby regions using multi-epoch Very Long Baseline Array (VLBA) radio observations. Enough VLBA detections have now been obtained for 16 stellar systems in Ophiuchus to derive their parallax and proper motion. This leads to distance determinations for individual stars with an accuracy of 0.3 to a few percent. In addition, the orbits of 6 multiple systems were modelled by combining absolute positions with VLBA (and in some cases, near infrared) angular separations. Twelve stellar systems are located in the dark cloud Lynds 1688, the individual distances for this sample are highly consistent with one another, and yield a mean parallax for Lynds 1688 of $\varpi=7.28\pm0.06$ mas, corresponding to a distance $d=137.3\pm1.2$ pc. This represents an accuracy better than 1%. Three systems for which astrometric elements could be measured are located in the eastern streamer (Lynds 1689) and yield an estimate of $\varpi=6.79\pm0.16$ mas, corresponding to a distance $d=147.3\pm3.4$ pc. This suggests that the eastern streamer is located about 10 pc farther than the core, but this conclusion needs to be confirmed by observations (currently being collected) of additional sources in the eastern streamer. From the measured proper motions, we estimate the one-dimensional velocity dispersion in Lynds 1688 to be 2.8$\pm$1.8 and 3.0$\pm$2.0 ${\rm km~s}^{-1}$, in R.A. and DEC., respectively, these are larger than, but still consistent within $1\sigma$, with those found in other studies.

We report on new distances and proper motions to seven stars across the Serpens/Aquila complex. The observations were obtained as part of the Gould's Belt Distances Survey (GOBELINS) project between September 2013 and April 2016 with the Very Long Baseline Array (VLBA). One of our targets is the proto-Herbig AeBe object EC 95, which is a binary system embedded in the Serpens Core. For this system, we combined the GOBELINS observations with previous VLBA data to cover a total period of ~8 years, and derive the orbital elements and an updated source distance. The individual distances to sources in the complex are fully consistent with each other, and the mean value corresponds to a distance of $436.0\pm9.2$~pc for the Serpens/W40 complex. Given this new evidence, we argue that Serpens Main, W40 and Serpens South are physically associated and form a single cloud structure.

We present the results of the Gould's Belt Distances Survey (GOBELINS) of young star forming regions towards the Orion Molecular Cloud Complex. We detected 36 YSOs with the Very Large Baseline Array (VLBA), 27 of which have been observed in at least 3 epochs over the course of 2 years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions towards these stars to study the structure and kinematics of the Complex. We measured a distance of 388$\pm$5 pc towards the Orion Nebula Cluster, 428$\pm$10 pc towards the southern portion L1641, 388$\pm$10 pc towards NGC 2068, and roughly $\sim$420 pc towards NGC 2024. Finally, we observed a strong degree of plasma radio scattering towards $\lambda$ Ori.

We present the results of a closure phase analysis of 3 mm very long baseline interferometry (VLBI) measurements performed on Sagittarius A* (Sgr A*). We have analyzed observations made in May 2015 using the Very Long Baseline Array, the Robert C. Byrd Green Bank Telescope and the Large Millimeter Telescope Alfonso Serrano and obtained non-zero closure phase measurements on several station triangles - indicative of a non-point-symmetric source structure. The data are fitted with an asymmetric source structure model in Sgr A*, represented by a simple two-component model, which favours a fainter component due East of the main source. This result is discussed in light of a scattering screen with substructure or an intrinsically asymmetric source.

We report the discovery of variable radio emission associated with the T Tauri star, V830 Tau, which was recently shown to host a hot Jupiter companion. Very Large Array observations at a frequency of 6 GHz reveal a detection on 01 May 2011 with a flux density $919 \pm 26\ \mu$Jy, along with non-detections in two other epochs at $<66$ and $<150\ \mu$Jy. Additionally, Very Long Baseline Array observations include one detection and one non-detection at comparable sensitivity, demonstrating that the emission is nonthermal in origin. The emission is consistent with the gyro-synchrotron or synchrotron mechanism from a region with a magnetic field ${\stackrel{\scriptstyle >}{\scriptstyle \sim}} 30$ G, and is likely driven by an energetic event such as magnetic reconnection that accelerated electrons. With the limited data we have, we are not able to place any constraint on the relationship between the radio emission and the rotational or orbital properties of V830 Tau. This is the first detection of radio emission from a non-degenerate star known to host an exoplanet.

We present a series of sixteen Very Long Baseline Array (VLBA) high angular resolution observations of a cluster of suspected low-mass young stars in the Monoceros R2 region. Four compact and highly variable radio sources are detected; three of them in only one epoch, the fourth one a total of seven times. This latter source is seen in the direction to the previously known \UCHII\ region VLA~1, and has radio properties that resemble those of magnetically active stars; we shall call it VLA~1$^\star$. We model its displacement on the celestial sphere as a combination of proper motion and trigonometric parallax. The fit obtained using a uniform proper motion yields a parallax $\varpi$ = 1.10 $\pm$ 0.18 mas, but with a fairly high post-fit dispersion. If acceleration terms (probably due to an undetected companion) are included, the quality of the fit improves dramatically, and the best estimate of the parallax becomes $\varpi$ = 1.12 $\pm$ 0.05 mas. The magnitude of the fitted acceleration suggest an orbital period of order a decade. The measured parallax corresponds to a distance $d$ = 893$^{+44}_{-40}$ pc, in very good agreement with previous, indirect, determinations.

The radio emission from Sgr A$^\ast$ is thought to be powered by accretion onto a supermassive black hole of $\sim\! 4\times10^6~ \rm{M}_\odot$ at the Galactic Center. At millimeter wavelengths, Very Long Baseline Interferometry (VLBI) observations can directly resolve the bright innermost accretion region of Sgr A$^\ast$. Motivated by the addition of many sensitive, long baselines in the north-south direction, we developed a full VLBI capability at the Large Millimeter Telescope Alfonso Serrano (LMT). We successfully detected Sgr A$^\ast$ at 3.5~mm with an array consisting of 6 Very Long Baseline Array telescopes and the LMT. We model the source as an elliptical Gaussian brightness distribution and estimate the scattered size and orientation of the source from closure amplitude and self-calibration analysis, obtaining consistent results between methods and epochs. We then use the known scattering kernel to determine the intrinsic two dimensional source size at 3.5 mm: $(147\pm7~\mu\rm{as}) \times (120\pm12~\mu\rm{as})$, at position angle $88^\circ\pm7^\circ$ east of north. Finally, we detect non-zero closure phases on some baseline triangles, but we show that these are consistent with being introduced by refractive scattering in the interstellar medium and do not require intrinsic source asymmetry to explain.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

We present multi-epoch, large-scale ($\sim$ 2000 arcmin${}^2$), fairly deep ($\sim$ 16 $\mu$Jy), high-resolution ($\sim$ 1") radio observations of the Perseus star-forming complex obtained with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. These observations were mainly focused on the clouds NGC 1333 and IC 348, although we also observed several fields in other parts of the Perseus complex. We detect a total of 206 sources, 42 of which are associated with young stellar objects (YSOs). The radio properties of about 60% of the YSOs are compatible with a non-thermal radio emission origin. Based on our sample, we find a fairly clear relation between the prevalence of non-thermal radio emission and evolutionary status of the YSOs. By comparing our results with previously reported X-ray observations, we show that YSOs in Perseus follow a G\"udel-Benz relation with $\kappa$ = 0.03 consistent with other regions of star formation. We argue that most of the sources detected in our observations but not associated with known YSOs are extragalactic, but provide a list of 20 unidentified radio sources whose radio properties are consistent with being YSO candidates. Finally we also detect 5 sources with extended emission features which can clearly be associated with radio galaxies.

We investigate the internal and relative motions of the Taurus and Ophiuchus star-forming regions using a sample of young stars with accurately measured radial velocities and proper motions. We find no evidence for expansion or contraction of the Taurus complex, but a clear indication for a global rotation, resulting in velocity gradients, this suggests a common origin, possibly related to that of Gould's Belt.

We present deep ($\sim 17~\mu$Jy) radio continuum observations of the Serpens molecular cloud, the Serpens south cluster, and the W40 region obtained using the Very Large Array in its A configuration. We detect a total of 146 sources, 29 of which are young stellar objects (YSOs), 2 are BV stars and 5 more are associated with phenomena related to YSOs. Based on their radio variability and spectral index, we propose that about 16 of the remaining 110 unclassified sources are also YSOs. For approximately 65% of the known YSOs detected here as radio sources, the emission is most likely non-thermal, and related to stellar coronal activity. As also recently observed in Ophiuchus, our sample of YSOs with X-ray counterparts lies below the fiducial G\"udel & Benz relation. Finally, we analyze the proper motions of 9 sources in the W40 region. This allows us to better constrain the membership of the radio sources in the region.

We present a multi-epoch radio study of the Taurus-Auriga star-forming complex made with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. We detect a total of 610 sources, 59 of which are related to young stellar objects and 18 to field stars. The properties of 56\% of the young stars are compatible with non-thermal radio emission. We also show that the radio emission of more evolved young stellar objects tends to be more non-thermal in origin and, in general, that their radio properties are compatible with those found in other star forming regions. By comparing our results with previously reported X-ray observations, we notice that young stellar objects in Taurus-Auriga follow a G\"{u}del-Benz relation with $\kappa$=0.03, as we previously suggested for other regions of star formation. In general, young stellar objects in Taurus-Auriga and in all the previous studied regions seem to follow this relation with a dispersion of $\sim1$ dex. Finally, we propose that most of the remaining sources are related with extragalactic objects but provide a list of 46 unidentified radio sources whose radio properties are compatible with a YSO nature.

We present results from a high-sensitivity (60 $\mu$Jy), large-scale (2.26 square degree) survey obtained with the Karl G. Jansky Very Large Array as part of the Gould's Belt Survey program. We detected 374 and 354 sources at 4.5 and 7.5 GHz, respectively. Of these, 148 are associated with previously known Young Stellar Objects (YSOs). Another 86 sources previously unclassified at either optical or infrared wavelengths exhibit radio properties that are consistent with those of young stars. The overall properties of our sources at radio wavelengths such as their variability and radio to X-ray luminosity relation are consistent with previous results from the Gould's Belt Survey. Our detections provide target lists for followup VLBA radio observations to determine their distances as YSOs are located in regions of high nebulosity and extinction, making it difficult to measure optical parallaxes.

We present large-scale ($\sim$ 2000 square arcminutes), deep ($\sim$ 20 $\mu$Jy), high-resolution ($\sim$ 1$''$) radio observations of the Ophiuchus star-forming complex obtained with the Karl G.\ Jansky Very Large Array at $\lambda$ = 4 and 6 cm. In total, 189 sources were detected, 56 of them associated with known young stellar sources, and 4 with known extragalactic objects; the other 129 remain unclassified, but most of them are most probably background quasars. The vast majority of the young stars detected at radio wavelengths have spectral types K or M, although we also detect 4 objects of A/F/B types and 2 brown dwarf candidates. At least half of these young stars are non-thermal (gyrosynchrotron) sources, with active coronas characterized by high levels of variability, negative spectral indices, and (in some cases) significant circular polarization. As expected, there is a clear tendency for the fraction of non-thermal sources to increase from the younger (Class 0/I or flat spectrum) to the more evolved (Class III or weak line T Tauri) stars. The young stars detected both in X-rays and at radio wavelengths broadly follow a G\"udel-Benz relation, but with a different normalization than the most radio-active types of stars. Finally, we detect a $\sim$ 70 mJy compact extragalactic source near the center of the Ophiuchus core, which should be used as gain calibrator for any future radio observations of this region.