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.

Due to its powerful radiation over the entire electromagnetic spectrum and its radio jet activity, the blazar 3C273 offers the opportunity to study the physics of $\gamma$-ray emission from active galactic nuclei. Since a historically strong outburst in 2009, 3C273 showed relatively weak emission in the gamma-ray band over multiple years. However, recent Fermi-Large Area Telescope observations indicate elevated activity during 2015-2019. We aim at constraining the origin of the gamma-ray outbursts towards 3C273 and investigate their connection to the parsec-scale jet. We generate Fermi-LAT gamma-ray light curves with multiple binning intervals and study the spectral properties of the gamma-ray emission. Using a 3-mm ALMA light curve, we study the correlation between radio and gamma-ray emission. Relevant activity in the parsec-scale jet of 3C273 is investigated with 7-mm VLBA observations obtained close in time to notable gamma-ray outbursts. We find two prominent gamma-ray outbursts in 2016 (MJD 57382) and 2017 (MJD 57883) accompanied by mm-wavelength flaring activity. The gamma-ray photon index time series show a weak hump-like feature around the gamma-ray outbursts. The monthly gamma-ray flux-index plot indicates a transition from softer-when-brighter to harder-when-brighter at $1.03\times10^{-7}\rm\,ph\,cm^{-2}\,s^{-1}$. A significant correlation between the gamma-ray and mm-wavelength emission is found, with the radio lagging the gamma-rays by about 105-112 days. The 43-GHz jet images reveal the known stationary features (i.e., the core, S1, and S2) in a region upstream of the jet. We find indication for a propagating disturbance and a polarized knot between the stationary components around the times of both gamma-ray outbursts.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

The flat spectrum radio quasar 4C 38.41 showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional monitoring observations at various bands across the electromagnetic spectrum. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components and the evolution of the physical properties seem to be in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1 mG, do not seem to dramatically change during the flares, and are much smaller, by a factor 10,000, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

PKS 1510--089 is a bright and active $\gamma$-ray source that showed strong and complex $\gamma$-ray flares in mid-2015 during which the Major Atmospheric Gamma Imaging Cherenkov telescopes detected variable very high energy (VHE; photon energies $>$100 GeV) emission. We present long-term multi-frequency radio, optical, and $\gamma$-ray light curves of PKS 1510--089 from 2013 to 2018, and results of an analysis of the jet kinematics and linear polarization using 43 GHz Very Long Baseline Array data observed between late 2015 and mid-2017. We find that a strong radio flare trails the $\gamma$-ray flares in 2015, showing an optically thick spectrum at the beginning and becoming optically thin over time. Two laterally separated knots of emission are observed to emerge from the radio core nearly simultaneously during the $\gamma$-ray flares. We detect an edge-brightened linear polarization near the core in the active jet state in 2016, similar to the quiescent jet state in 2008--2013. These observations indicate that the $\gamma$-ray flares may originate from compression of the knots by a standing shock in the core and the jet might consist of multiple complex layers showing time-dependent behavior, rather than of a simple structure of a fast jet spine and a slow jet sheath.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.

A long standing problem in the study of Active Galactic Nuclei (AGNs) is that the observed VLBI core is in fact a blending of the actual AGN core (classically defined by the $\tau=1$ surface) and the upstream regions of the jet or optically thin emitting region flows. This blending may cause some biases towards the observables of the core, such as its flux density, size or brightness temperature, which may lead to misleading interpretation of the derived quantities and physics. We study the effects of such blending under the view of the Korean VLBI Network (KVN) for a sample of AGNs at 43 GHz by comparing their observed properties with observations with the Very Large Baseline Array (VLBA). Our results suggest that the observed core sizes are a factor $\sim11$ larger than these of VLBA, which is similar to the factor expected by considering the different resolutions of the two facilities. We suggest the use of this factor to consider blending effects in KVN measurements. Other parameters, such as flux density or brightness temperature, seem to possess a more complicated dependence.

We study Faraday rotation in the jet of M87 inside the Bondi radius using eight Very Long Baseline Array data sets, one at 8 GHz, four at 5 GHz, and three at 2 GHz. We obtain Faraday rotation measures (RMs) measured across the bandwidth of each data set. We find that the magnitude of RM systematically decreases with increasing distance from the black hole from 5,000 to 200,000 Schwarzschild radii. The data, showing predominantly negative RM sign without significant difference of the RMs on the northern and southern jet edges, suggest that the spatial extent of the Faraday screen is much larger than the jet. We apply models of hot accretion flows, thought to be prevalent in active galactic nuclei having relatively low luminosity such as M87, and find that the decrease of RM is described well by a gas density profile $\rho \propto r^{-1}$. This behavior matches the theoretically expected signature of substantial winds, nonrelativistic un-collimated gas outflows from hot accretion flows, which is consistent with the results of various numerical simulations. The pressure profile inferred from the density profile is flat enough to collimate the jet, which can result in gradual acceleration of the jet in a magneto-hydrodynamical process. This picture is in good agreement with the observed gradual collimation and acceleration of the M87 jet inside the Bondi radius. The dominance of negative RMs suggests that jet and wind axis are misaligned such that the jet emission exposes only one side of the toroidal magnetic fields permeating the winds.

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.

Recent Fermi-Large Area Telescope (LAT) light curves indicate an active $\gamma$-ray state spanning about five months from 2016 June to 2016 October in the BL Lac object 1749+096 (OT 081). During this period, we find two notable $\gamma$-ray events: an exceptionally strong outburst followed by a significant enhancement (local peak). In this study, we analyze multi-waveband light curves (radio, optical, X-ray, and $\gamma$-ray) plus very-long baseline interferometry (VLBI) data to investigate the nature of the $\gamma$-ray events. The $\gamma$-ray outburst coincides with flux maxima at longer wavelengths. We find a spectral hardening of the $\gamma$-ray photon index during the $\gamma$-ray outburst. The photon index shows a transition from a softer-when-brighter to a harder-when-brighter trend at around 1.8 $\times$ $10^{-7}$ ph cm$^{-2}$ s$^{-1}$. We see indication that both the $\gamma$-ray outburst and the subsequent enhancement precede the propagation of a polarized knot in a region near the VLBI core. The highest polarized intensity, 230\,mJy, and an electric vector position angle rotation, by $\sim$32$^{\circ}$, are detected about 12 days after the $\gamma$-ray outburst. We conclude that both $\gamma$-ray events are caused by the propagation of a disturbance in the mm-wave core.

We study the linear polarization of the radio cores of eight blazars simultaneously at 22, 43, and 86 GHz with observations obtained by the Korean VLBI Network (KVN) in three epochs between late 2016 and early 2017 in the frame of the Plasma-physics of Active Galactic Nuclei (PAGaN) project. We investigate the Faraday rotation measure (RM) of the cores; the RM is expected to increase with observing frequency if core positions depend on frequency due to synchrotron self-absorption. We find a systematic increase of RMs at higher observing frequencies in our targets. The RM--$\nu$ relations follow power-laws with indices distributed around 2, indicating conically expanding outflows serving as Faraday rotating media. Comparing our KVN data with contemporaneous optical polarization data from the Steward Observatory for a few sources, we find indication that the increase of RM with frequency saturates at frequencies of a few hundreds GHz. This suggests that blazar cores are physical structures rather than simple $\tau=1$ surfaces. A single region, e.g. a recollimation shock, might dominate the jet emission downstream of the jet launching region. We detect a sign change in the observed RMs of CTA 102 on a time scale of $\approx$1 month, which might be related to new superluminal components emerging from its core undergoing acceleration/deceleration and/or bending. We see indication for quasars having higher core RMs than BL Lac objects, which could be due to denser inflows/outflows in quasars.

The flat spectrum radio quasar 1633+382 (4C~38.41) showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional radio monitoring observations with the OVRO 40m telescope, the Boston University VLBI program, and the Submillimeter Array. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components. Analysis of the spectral energy distribution indicates a higher turnover frequency after the flaring events. The evolution of the flare in the turnover frequency - turnover flux density plane probes the adiabatic losses in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1mG, do not seem to dramatically change during the flares, and are much smaller, by a factor $10^4$, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

3C 84 (NGC 1275) is a well-studied mis-aligned Active Galactic Nucleus (AGN), which has been active in Gamma rays since at least 2008. We have monitored the source at four wavelengths (14 mm, 7 mm, 3 mm and 2 mm) using the Korean VLBI network (KVN) since 2013 as part of the interferometric monitoring of $\gamma$-ray bright AGN (iMOGABA) program. 3C 84 exhibits bright radio emission both near the central supermassive black hole (SMBH) feature known as C1 and from a moving feature located to the south known as C3. Other facilities have also detected these short-term variations above a slowly rising trend at shorter wavelengths, such as in Gamma ray and 1 mm total intensity light-curves. We find that the variations in the $\gamma$ rays and 1 mm total intensity light-curves are correlated, with the $\gamma$ rays leading and lagging the radio emission. Analysis of the 2 mm KVN data shows that both the Gamma rays and 1 mm total intensity short-term variations are better correlated with the SMBH region than C3, likely placing the short-term variations in C1. We interpret the emission as being due to the random alignment of spatially separated emission regions. We place the slowly rising trend in C3, consistent with previous results. Additionally, we report that since mid-2015, a large mm-wave radio flare has been occurring in C3, with a large Gamma ray flare coincident with the onset of this flare at all radio wavelengths.

Atmospheric propagation effects at millimeter wavelengths can significantly alter the phases of radio signals and reduce the coherence time, putting tight constraints on high frequency Very Long Baseline Interferometry (VLBI) observations. In previous works, it has been shown that non-dispersive (e.g. tropospheric) effects can be calibrated with the frequency phase transfer (FPT) technique. The coherence time can thus be significantly extended. Ionospheric effects, which can still be significant, remain however uncalibrated after FPT as well as the instrumental effects. In this work, we implement a further phase transfer between two FPT residuals (i.e. so-called FPT-square) to calibrate the ionospheric effects based on their frequency dependence. We show that after FPT-square, the coherence time at 3 mm can be further extended beyond 8~hours, and the residual phase errors can be sufficiently canceled by applying the calibration of another source, which can have a large angular separation from the target (>20 deg) and significant temporal gaps. Calibrations for all-sky distributed sources with a few calibrators are also possible after FPT-square. One of the strengths and uniqueness of this calibration strategy is the suitability for high-frequency all-sky survey observations including very weak sources. We discuss the introduction of a pulse calibration system in the future to calibrate the remaining instrumental effects and allowing the possibility of imaging the source structure at high frequencies with FPT-square, where all phases are fully calibrated without involving any additional sources.

We present multi-frequency simultaneous VLBI radio observations of the flat spectrum radio quasar 1633+382 (4C~38.41) as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program combined with additional observations in radio, optical, X-rays and $\gamma-$rays carried out between the period 2012 March - 2015 August. The monitoring of this source reveals a significant long-lived increase in its activity since approximately two years in the radio bands, which correlates with a similar increase on all other bands from sub-millimeter to $\gamma-$rays. A significant correlation is also found between radio fluxes and simultaneous spectral indices during this period. The study of the discrete correlation function (DCF) indicates time lags smaller than the $\sim40$ days uncertainties among both radio bands and also high-energy bands, and a time lag of $\sim$70 days, with $\gamma-$rays leading radio. We interpret that the high-energy and radio fluxes are arising from different emitting regions, located at $1\pm12$ and $40\pm13$ pc from the central engine respectively.

We present a study of the inexplicit connection between radio jet activity and gamma-ray emission of BL Lacertae (BL Lac; 2200+420). We analyze the long-term millimeter activity of BL Lac via interferometric observations with the Korean VLBI Network (KVN) obtained at 22, 43, 86, and 129 GHz simultaneously over three years (from January 2013 to March 2016); during this time, two gamma-ray outbursts (in November 2013 and March 2015) can be seen in gamma-ray light curves obtained from Fermi observations. The KVN radio core is optically thick at least up to 86 GHz; there is indication that it might be optically thin at higher frequencies. To first order, the radio light curves decay exponentially over the time span covered by our observations, with decay timescales of 411+/-85 days, 352+/-79 days, 310+/-57 days, and 283+/-55 days at 22, 43, 86, and 129 GHz, respectively. Assuming synchrotron cooling, a cooling time of around one year is consistent with magnetic field strengths B~2microT and electron Lorentz factors gamma~10,000. Taking into account that our formal measurement errors include intrinsic variability and thus over-estimate the statistical uncertainties, we find that the decay timescale tau scales with frequency nu like tau~nu^{-0.2}. This relation is much shallower than the one expected from opacity effects (core shift), but in agreement with the (sub-)mm radio core being a standing recollimation shock. We do not find convincing radio flux counterparts to the gamma$ray outbursts. The spectral evolution is consistent with the `generalized shock model' of Valtaoja et al. (1992). A temporary increase in the core opacity and the emergence of a knot around the time of the second gamma-ray event indicate that this gamma-ray outburst might be an `orphan' flare powered by the `ring of fire' mechanism.

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in North America and Europe and to expand existing global mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b) investigating the dynamics of circumstellar gas in star-forming regions and late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In the first half of this White Paper, we highlight scientific accomplishments of the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are described, and future visions of our VLBI science are briefly presented. Maser science of star formation, stellar evolution, and physics of accretion disks around AGNs are also discussed. A new vision for conducting multi-transition maser studies using mm/submm VLBI together with the Atacama Large Millimeter/Submillimeter Array (ALMA) is described. In the second half of this White Paper, we describe the EA community's vision for using mm/submm VLBI arrays in the framework or extended version of the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the aforementioned goal will maximize the overall scientific outcomes of mm/submm VLBI in the world.

We study the long-term ($\approx$30 years) radio variability of 43 radio bright AGNs by exploiting the data base of the University of Michigan Radio Astronomy Observatory (UMRAO) monitoring program. We model the periodograms (temporal power spectra) of the observed lightcurves as simple power-law noise (red noise, spectral power $P(f)\propto f^{-\beta}$) using Monte Carlo simulations, taking into account windowing effects (red-noise leak, aliasing). The power spectra of 39 (out of 43) sources are in good agreement with the models, yielding a range in power spectral index ($\beta$) from $\approx$1 to $\approx$3. We fit a Gaussian function to each flare in a given lightcurve to obtain the flare duration. We discover a correlation between $\beta$ and the median duration of the flares. We use the derivative of a lightcurve to obtain a characteristic variability timescale which does not depend on the assumed functional form of the flares, incomplete fitting, and so on. We find that, once the effects of relativistic Doppler boosting are corrected for, the variability timescales of our sources are proportional to the accretion rate to the power of $0.25\pm 0.03$ over five orders of magnitude in accretion rate, regardless of source type. We further find that modeling the periodograms of four of our sources requires the assumption of broken powerlaw spectra. From simulating lightcurves as superpositions of exponential flares we conclude that strong overlap of flares leads to featureless simple power-law periodograms of AGNs at radio wavelengths in most cases.

We present results of single-epoch very long baseline interferometry (VLBI) observations of gamma-ray bright active galactic nuclei (AGNs) using the Korean VLBI Network (KVN) at 22, 43, 86, and 129~GHz bands, which are part of a KVN key science program, Interferometric Monitoring of Gamma-ray Bright AGNs (iMOGABA). We selected a total of 34 radio-loud AGNs of which 30 sources are gamma-ray bright AGNs with flux densities of $>6\times10^{-10}$~ph~cm$^{-2}$~s$^{-1}$. Single-epoch multi-frequency VLBI observations of the target sources were conducted during a 24-hr session on 2013 November 19 and 20. All observed sources were detected and imaged at all frequency bands with or without a frequency phase transfer technique which enabled the imaging of 12 faint sources at 129~GHz, except for one source. Many of the target sources are resolved on milliarcsecond scales, yielding a core-jet structure with the VLBI core dominating the synchrotron emission on the milliarcsecond scale. CLEAN flux densities of the target sources are 0.43-28~Jy, 0.32-21~Jy, 0.18-11~Jy, and 0.35-8.0~Jy in the 22, 43, 86, and 129~GHz bands, respectively.

Early-type galaxies (ETGs) are supposed to follow the virial relation $M = k_e \sigma_*^2 R_e / G$, with $M$ being the mass, $\sigma_*$ being the stellar velocity dispersion, $R_e$ being the effective radius, $G$ being Newton's constant, and $k_e$ being the virial factor, a geometry factor of order unity. Applying this relation to (a) the ATLAS3D sample of Cappellari et al. (2013) and (b) the sample of Saglia et al. (2016) gives ensemble-averaged factors $\langle k_e\rangle =5.15\pm0.09$ and $\langle k_e\rangle =4.01\pm0.18$, respectively, with the difference arising from different definitions of effective velocity dispersions. The two datasets reveal a statistically significant tilt of the empirical relation relative to the theoretical virial relation such that $M\propto(\sigma_*^2R_e)^{0.92}$. This tilt disappears when replacing $R_e$ with the semi-major axis of the projected half-light ellipse, $a$. All best-fit scaling relations show zero intrinsic scatter, implying that the mass plane of ETGs is fully determined by the virial relation. Whenever a comparison is possible, my results are consistent with, and confirm, the results by Cappellari et al. (2013). The difference between the relations using either $a$ or $R_e$ arises from a known lack of highly elliptical high-mass galaxies; this leads to a scaling $(1-\epsilon) \propto M^{0.12}$, with $\epsilon$ being the ellipticity and $R_e = a\sqrt{1-\epsilon}$. Accordingly, $a$, not $R_e$, is the correct proxy for the scale radius of ETGs. By geometry, this implies that early-type galaxies are axisymmetric and oblate in general, in agreement with published results from modeling based on kinematics and light distributions.

Understanding the interaction of primordial gravitational waves (GWs) with the Cosmic Microwave Background (CMB) plasma is important for observational cosmology. In this article, we provide an analysis of an effect apparently overlooked as yet. We consider a single free electric charge and suppose that it can be agitated by primordial GWs propagating through the CMB plasma, resulting in periodic, regular motion along particular directions. Light reflected by the charge will be partially polarized, and this will imprint a characteristic pattern on the CMB. We study this effect by considering a simple model in which anisotropic incident electromagnetic (EM) radiation is rescattered by a charge sitting in spacetime perturbed by GWs and becomes polarized. As the charge is driven to move along particular directions, we calculate its dipole moment to determine the leading-order rescattered EM radiation. The Stokes parameters of the rescattered radiation exhibit a net linear polarization. We investigate how this polarization effect can be schematically represented out of the Stokes parameters. We work out the representations of gradient modes (E-modes) and curl modes (B-modes) to produce polarization maps. Although the polarization effect results from GWs, we find that its representations, the E- and B-modes, do not practically reflect the GW properties such as strain amplitude, frequency and polarization states.

Active Galactic Nuclei (AGN) with bright radio jets offer the opportunity to study the structure of and physical conditions in relativistic outflows. For such studies, multi-frequency polarimetric very long baseline interferometric (VLBI) observations are important as they directly probe particle densities, magnetic field geometries, and several other parameters. We present results from first-epoch data obtained by the Korean VLBI Network (KVN) within the frame of the Plasma Physics of Active Galactic Nuclei (PAGaN) project. We observed seven radio-bright nearby AGN at frequencies of 22, 43, 86, and 129 GHz in dual polarization mode. Our observations constrain apparent brightness temperatures of jet components and radio cores in our sample to $>10^{8.01}$ K and $>10^{9.86}$ K, respectively. Degrees of linear polarization $m_{L}$ are relatively low overall: less than 10%. This indicates suppression of polarization by strong turbulence in the jets. We found an exceptionally high degree of polarization in a jet component of BL Lac at 43 GHz, with $m_{L} \sim$ 40%. Assuming a transverse shock front propagating downstream along the jet, the shock front being almost parallel to the line of sight can explain the high degree of polarization.

Active galactic nuclei (AGN) are known for irregular variability on all time scales, down to intra-day variability with relative variations of a few percent within minutes to hours. On such short timescales, unexplored territory, such as the possible existence of a shortest characteristic time scale of activity and the shape of the high frequency end of AGN power spectra, still exists. We present the results of AGN single-dish fast photometry performed with the Korean VLBI Network (KVN). Observations were done in a "anti-correlated" mode using two antennas, with always at least one antenna pointing at the target. This results in an effective time resolution of less than three minutes. We used all four KVN frequencies, 22, 43, 86, and 129 GHz, in order to trace spectral variability, if any. We were able to derive high-quality light curves for 3C 111, 3C 454.3, and BL Lacertae at 22 and 43 GHz, and for 3C 279 at 86 GHz, between May 2012 and April 2013. We performed a detailed statistical analysis in order to assess the levels of variability and the corresponding upper limits. We found upper limits on flux variability ranging from $\sim$1.6% to $\sim$7.6%. The upper limits on the derived brightness temperatures exceed the inverse Compton limit by three to six orders of magnitude. From our results, plus comparison with data obtained by the University of Michigan Radio Astronomy Observatory, we conclude that we have not detected source-intrinsic variability which would have to occur at sub-per cent levels.

We report first results from KVN and VERA Array (KaVA) VLBI observations obtained in the frame of our Plasma-physics of Active Galactic Nuclei (PAGaN) project. We observed eight selected AGN at 22 and 43 GHz in single polarization (LCP) between March 2014 and April 2015. Each source was observed for 6 to 8 hours per observing run to maximize the $uv$ coverage. We obtained a total of 15 deep high-resolution images permitting the identification of individual circular Gaussian jet components and three spectral index maps of BL Lac, 3C 111 and 3C 345 from simultaneous dual-frequency observations. The spectral index maps show trends in agreement with general expectations -- flat core and steep jets -- while the actual value of the spectral index for jets shows indications for a dependence on AGN type. We analyzed the kinematics of jet components of BL Lac and 3C 111, detecting superluminal proper motions with maximum apparent speeds of about $5c$. This constrains the lower limits of the intrinsic component velocities to $\sim0.98c$ and the upper limits of the angle between jet and line of sight to $\sim$20$\deg$. In agreement with global jet expansion, jet components show systematically larger diameters $d$ at larger core distances $r$, following the global relation $d\approx0.2r$, albeit within substantial scatter.

The Interferometric Monitoring of Gamma-ray Bright Active galactic nuclei (iMOGABA) program provides not only simultaneous multifrequency observations of bright gamma-ray detected active galactic nuclei (AGN), but also covers the highest Very Large Baseline Interferometry (VLBI) frequencies ever being systematically monitored, up to 129 GHz. However, observation and imaging of weak sources at the highest observed frequencies is very challenging. In the second paper in this series, we evaluate the viability of the frequency phase transfer technique to iMOGABA in order to obtain larger coherence time at the higher frequencies of this program (86 and 129 GHz) and image additional sources that were not detected using standard techniques. We find that this method is applicable to the iMOGABA program even under non-optimal weather conditions.

We report the first detection of linearly polarized emission at an observing wavelength of 350 $\mu$m from the radio-loud Active Galactic Nucleus 3C 279. We conducted polarization observations for 3C 279 using the SHARP polarimeter in the Caltech Submillimeter Observatory on 2014 March 13 and 14. For the first time, we detected the linear polarization with the degree of polarization of 13.3\%$\pm$3.4\% ($3.9\sigma$) and the Electric Vector Position Angle (EVPA) of 34.7$^\circ\pm5.6^\circ$. We also observed 3C~279 simultaneously at 13, 7, and 3.5 mm in dual polarization with the Korean very long baseline interferometry (VLBI) Network on 2014 March 6 (single dish) and imaged in milliarcsecond (mas) scales at 13, 7, 3.5, and 2.3 mm on March 22 (VLBI). We found that the degree of linear polarization increases from 10\% to 13\% at 13 mm to 350 $\mu$m and the EVPAs at all observing frequencies are parallel within $<10^\circ$ to the direction of the jet at mas scale, implying that the integrated magnetic fields are perpendicular to the jet in the innermost regions. We also found that the Faraday rotation measures RM are in a range of $-6.5\times10^2 \sim -2.7\times10^3$ rad m$^{-2}$ between 13 and 3.5 mm, and are scaled as a function of wavelength: $|{\rm RM}|\propto\lambda^{-2.2}$. These results indicate that the mm and sub-mm polarization emission are generated in the compact jet within 1~mas scale and affected by a Faraday screen in or in the close proximity of the jet.

The masses of supermassive black holes in active galactic nuclei (AGN) can be derived spectroscopically via virial mass estimators based on selected broad optical/ultraviolet emission lines. These estimates commonly use the line width as a proxy for the gas speed and the monochromatic continuum luminosity as a proxy for the radius of the broad line region. However, if the size of the broad line region scales with bolometric rather than monochromatic AGN luminosity, mass estimates based on different emission lines will show a systematic discrepancy which is a function of the color of the AGN continuum. This has actually been observed in mass estimates based on H-alpha / H-beta and C IV lines, indicating that AGN broad line regions indeed scale with bolometric luminosity. Given that this effect seems to have been overlooked as yet, currently used single-epoch mass estimates are likely to be biased.

Massive gravity provides a natural solution for the dark energy problem of cosmology and is also a candidate for resolving the dark matter problem. I demonstrate that, assuming reasonable scaling relations, massive gravity can provide for Milgrom's law of gravity (or "modified Newtonian dynamics") which is known to remove the need for particle dark matter from galactic dynamics. Milgrom's law comes with a characteristic acceleration, Milgrom's constant, which is observationally constrained to $a_0\approx1.1\times10^{-10}$ m/s$^2$. In the derivation presented here, this constant arises naturally from the cosmologically required mass of gravitons like $a_0\propto c\sqrt{\Lambda}\propto c H_0\sqrt{3\Omega_{\Lambda}}$, with $\Lambda$, $H_0$, and $\Omega_{\Lambda}$ being the cosmological constant, the Hubble constant, and the third cosmological parameter, respectively. My derivation suggests that massive gravity could be the mechanism behind both, dark matter and dark energy.

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.

Intensity interferometry, based on the Hanbury Brown-Twiss effect, is a simple and inexpensive method for optical interferometry at microarcsecond angular resolutions; its use in astronomy was abandoned in the 1970s because of low sensitivity. Motivated by recent technical developments, we argue that the sensitivity of large modern intensity interferometers can be improved by factors up to approximately 25,000, corresponding to 11 photometric magnitudes, compared to the pioneering Narrabri Stellar Interferometer. This is made possible by (i) using avalanche photodiodes (APD) as light detectors, (ii) distributing the light received from the source over multiple independent spectral channels, and (iii) use of arrays composed of multiple large light collectors. Our approach permits the construction of large (with baselines ranging from few kilometers to intercontinental distances) optical interferometers at the cost of (very) long-baseline radio interferometers. Realistic intensity interferometer designs are able to achieve limiting R-band magnitudes as good as ~14, sufficient for spatially resolved observations of main-sequence O-type stars in the Magellanic Clouds. Multi-channel intensity interferometers can address a wide variety of science cases: (i) linear radii, effective temperatures, and luminosities of stars; (ii) mass-radius relationships of compact stellar remnants; (iii) stellar rotation; (iv) stellar convection and the interaction of stellar photospheres and magnetic fields; (v) the structure and evolution of multiple stars; (vi) direct measurements of interstellar distances; (vii) the physics of gas accretion onto supermassive black holes; and (viii) calibration of amplitude interferometers by providing a sample of calibrator stars.

We present a GUI-based interactive Python program, VIMAP, which generates radio spectral index maps of active galactic nuclei (AGN) from Very Long Baseline Interferometry (VLBI) maps obtained at different frequencies. VIMAP is a handy tool for the spectral analysis of synchrotron emission from AGN jets, specifically of spectral index distributions, turn-over frequencies, and core-shifts. In general, the required accurate image alignment is difficult to achieve because of a loss of absolute spatial coordinate information during VLBI data reduction (self-calibration) and/or intrinsic variations of source structure as function of frequency. These issues are overcome by VIMAP which in turn is based on the two-dimensional cross-correlation algorithm of Croke and Gabuzda (2008). In this paper, we briefly review the problem of aligning VLBI AGN maps, describe the workflow of VIMAP, and present an analysis of archival VLBI maps of the active nucleus 3C 120.

The jet production efficiency of radio galaxies can be quantified by comparison of their kinetic jet powers P_jet and Bondi accretion powers P_B. These two parameters are known to be related linearly, with the jet power resulting from the Bondi power by multiplication with an efficiency factor of order 1%. Using a recently published (Nemmen + Tchekhovskoy 2014) high-quality sample of 27 radio galaxies, I construct a P_B-P_jet diagram that includes information on optical AGN types as far as available. This diagram indicates that the jet production efficiency is a function of AGN type: Seyfert 2 galaxies seem to be systematically (with a false alarm probability of 0.043%) less efficient, by about one order of magnitude, in powering jets than Seyfert 1 galaxies, LINERs, or the remaining radio galaxies. This suggests an evolutionary sequence from Sy2s to Sy1s and LINERs, controlled by an interplay of jets on the one hand and dust and gas in galactic nuclei on the other hand. When taking this effect into account, the P_B-P_jet relation is probably much tighter intrinsically than currently assumed.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

We present the results of a time series analysis of the long-term radio lightcurves of four blazars: 3C 279, 3C 345, 3C 446, and BL Lacertae. We exploit the data base of the University of Michigan Radio Astronomy Observatory (UMRAO) monitoring program which provides densely sampled lightcurves spanning 32 years in time in three frequency bands located at 4.8, 8, and 14.5 GHz. Our sources show mostly flat or inverted (spectral indices -0.5 < alpha < 0) spectra, in agreement with optically thick emission. All lightcurves show strong variability on all time scales. Analyzing the time lags between the lightcurves from different frequency bands, we find that we can distinguish high-peaking flares and low-peaking flares in accord with the classification of Valtaoja et al. (1992). The periodograms (temporal power spectra) of the observed lightcurves are consistent with random-walk powerlaw noise without any indication of (quasi-)periodic variability. The fact that all four sources studied are in agreement with being random-walk noise emitters at radio wavelengths suggests that such behavior is a general property of blazars.

Since the 1930s, astronomical observations have accumulated evidence that our understanding of the dynamics of galaxies and groups of galaxies is grossly incomplete: assuming the validity of Newton's law of gravity on astronomical scales, the observed mass (stored in stars and interstellar gas) of stellar systems can account only for roughly 10% of the dynamical (gravitating) mass required to explain the high velocities of stars in those systems. The standard approach to this "missing mass problem" has been the postulate of "dark matter", meaning an additional, electromagnetically dark, matter component that provides the missing mass. However, direct observational evidence for dark matter has not been found to date. More importantly, astronomical observations obtained during the last decade indicate that dark matter cannot explain the kinematics of galaxies. Multiple observations show that the discrepancy between observed and dynamical mass is a function of gravitational acceleration (or field strength) but not of other parameters (size, rotation speed, etc.) of a galaxy; the mass discrepancy appears below a characteristic and universal acceleration ("Milgrom's constant"). Consequently, the idea of a modified law of gravity, specifically the ansatz of Modified Newtonian Dynamics (MOND), is becoming increasingly important in astrophysics. MOND has successfully predicted various important empirical relations of galaxy dynamics, including the famous Tully-Fisher and Faber-Jackson relations. MOND is found to be consistent with stellar dynamics from binary stars to clusters of galaxies, thus covering stellar systems spanning eight orders of magnitude in size and 14 orders of magnitude in mass. These developments have the potential to initiate a paradigm shift from dark matter to a modified law of gravity as the physical mechanism behind the missing mass problem.

Polarization is a basic property of light and is fundamentally linked to the internal geometry of a source of radiation. Polarimetry complements photometric, spectroscopic, and imaging analyses of sources of radiation and has made possible multiple astrophysical discoveries. In this article I review (i) the physical basics of polarization: electromagnetic waves, photons, and parameterizations; (ii) astrophysical sources of polarization: scattering, synchrotron radiation, active media, and the Zeeman, Goldreich-Kylafis, and Hanle effects, as well as interactions between polarization and matter (like birefringence, Faraday rotation, or the Chandrasekhar-Fermi effect); (iii) observational methodology: on-sky geometry, influence of atmosphere and instrumental polarization, polarization statistics, and observational techniques for radio, optical, and X/gamma wavelengths; and (iv) science cases for astronomical polarimetry: solar and stellar physics, planetary system bodies, interstellar matter, astrobiology, astronomical masers, pulsars, galactic magnetic fields, gamma-ray bursts, active galactic nuclei, and cosmic microwave background radiation.

We present the results of radio monitoring observations of Sgr A* at 7 mm (i.e. 43 GHz) with VLBI Exploration of Radio Astrometry (VERA), which is a VLBI array in Japan. VERA provides angular resolutions on millisecond scales, resolving structure within ~100 Schwarzschild radii of Sgr A* similar to Very Large Baseline Array (VLBA). We performed multi-epoch observations of Sgr A* in 2005 - 2008, and started monitoring it again with VERA from January 2013 for tracing the current G2 encounter event. Our preliminary results in 2013 show that Sgr A* on mas scales has been in ordinary state as of August 2013, although some fraction of the G2 cloud already passed pericenter of Sgr A* in April 2013. We will continue on monitoring Sgr A* with VERA and newly developed KaVA (KVN and VERA Array).

The empirical mass discrepancy-acceleration (MDA) relation of disk galaxies provides a key test for models of galactic dynamics. In terms of modified laws of gravity and/or inertia, the MDA relation quantifies the transition from Newtonian to modified dynamics at low centripetal accelerations a_c < 10^-10 m/s^2. As yet, neither dynamical models based on dark matter nor proposed modifications of the laws of gravity/inertia have predicted the functional form of the MDA relation. In this work, I revisit the MDA data and compare them to four different theoretical scaling laws. Three of these scaling laws are entirely empirical, the fourth one - the "simple mu" function of Modified Newtonian Dynamics - derives from a toy model of gravity based on massive gravitons (the "graviton picture"). All theoretical MDA relations comprise one free parameter of the dimension of an acceleration, Milgrom's constant a_M. I find that the "simple mu" function provides a good fit to the data free of notable systematic residuals and provides the best fit among the four scaling laws tested. The best-fit value of Milgrom's constant is a_M = (1.06 +/- 0.05)*10^-10 m/s^2. Given the successful prediction of the functional form of the MDA relation, plus an overall agreement with the observed kinematics of stellar systems spanning eight orders of magnitude in size and 14 orders of magnitude in mass, I conclude that the "graviton picture" is sufficient (albeit probably not a necessary or unique approach) to describe galactic dynamics on all scales well beyond the scale of the solar system. This suggests that, at least on galactic scales, gravity behaves as if it was mediated by massive particles.

Modified Newtonian Dynamics (MOND) is a possible solution for the missing mass problem in galactic dynamics; its predictions are in good agreement with observations in the limit of weak accelerations. However, MOND does not derive from a physical mechanism and does not make predictions on the transitional regime from Newtonian to modified dynamics; rather, empirical transition functions have to be constructed from the boundary conditions and comparison to observations. I compare the formalism of classical MOND to the scaling law derived from a toy model of gravity based on virtual massive gravitons (the "graviton picture") which I proposed recently. I conclude that MOND naturally derives from the "graviton picture" at least for the case of non-relativistic, highly symmetric dynamical systems. This suggests that - to first order - the "graviton picture" indeed provides a valid candidate for the physical mechanism behind MOND and gravity on galactic scales in general.

We probe the feasibility of high-frequency radio observations of very rapid flux variations in compact active galactic nuclei (AGN). Our study assumes observations at 230 GHz with a small 6-meter class observatory, using the SNU Radio Astronomical Observatory (SRAO) as example. We find that 33 radio-bright sources are observable with signal-to-noise ratios larger than ten. We derive statistical detection limits via exhaustive Monte Carlo simulations assuming (a) periodic, and (b) episodic flaring flux variations on time-scales as small as tens of minutes. We conclude that a wide range of flux variations is observable. This makes high-frequency radio observations - even with small observatories - a powerful probe of AGN intra-day variability; especially, those observations complement observations at lower radio frequencies with larger observatories like the Korean VLBI Network (KVN).

I present a simple scheme for the treatment of gravitational interactions on galactic scales. In analogy to known mechanisms of quantum field theory, I assume ad hoc that gravitation is mediated by virtual exchange particles - gravitons - with very small but non-zero masses. The resulting density and mass profiles are proportional to the mass of the gravitating body. The mass profile scales with the centripetal acceleration experienced by a test particle orbiting the central mass; this comes at the cost of postulating a universal characteristic acceleration a0 = 4.3*10^{-12} m/s^2 (or 8*pi*a0 = 1.1*10^{-10} m/s^2). The scheme predicts the asymptotic flattening of galactic rotation curves, the Tully-Fisher/Faber-Jackson relations, the mass discrepancy-acceleration relation of galaxies, the surface brightness-acceleration relation of galaxies, the kinematics of galaxy clusters, and "Renzo's rule" correctly; additional (dark) mass components are not required. Given that it is based on various ad-hoc assumptions, and given further limitations, the scheme I present is not yet a consistent theory of gravitation; rather, it is a "toy model" providing a convenient scaling law that simplifies the description of gravity on galactic scales.

We present a test of the emission statistics of active galactic nuclei (AGN), probing the connection between red-noise temporal power spectra and multi-modal flux distributions known from observations. We simulate AGN lightcurves under the assumption of uniform stochastic emission processes for different power-law indices of their power spectra. For sufficiently shallow slopes (power-law indices beta < 1), the flux distributions (histograms) of the resulting lightcurves are approximately Gaussian. For indices corresponding to steeper slopes (beta > 1), the flux distributions become multi-modal. This finding disagrees systematically with results of recent mm/radio observations. Accordingly, we conclude that the emission from AGN does not necessarily originate from uniform stochastic processes even if their power spectra suggest this. Possible mechanisms are transitions between different activity states and/or the presence of multiple, spatially disconnected, emission regions.

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 report a search for linear polarization in the active galactic nucleus (AGN) 3C 84 (NGC 1275) at observed frequencies of 239 GHz and 348 GHz, corresponding to rest-frame frequencies of 243 GHz and 354 GHz. We collected polarization data with the IRAM Plateau de Bure Interferometer via Earth rotation polarimetry. We do not detect linear polarization. Our analysis finds 3-sigma upper limits on the degree of polarization of 0.5% and 1.9% at 239 GHz and 348 GHz, respectively. We regard the influence of Faraday conversion as marginal, leading to expected circular polarizations <0.3%. Assuming depolarization by a local Faraday screen, we constrain the rotation measure, as well as the fluctuations therein, to be 10^6 rad/m^2. From this we estimate line-of-sight magnetic field strengths of >100 microG. Given the physical dimensions of 3C 84 and its observed structure, the Faraday screen appears to show prominent small-scale structure, with \DeltaRM > 10^6 rad/m^2 on projected spatial scales <1 pc.

We present an analysis of the linear polarization of six active galactic nuclei - 0415+379 (3C~111), 0507+179, 0528+134 (OG+134), 0954+658, 1418+546 (OQ+530), and 1637+574 (OS+562). Our targets were monitored from 2007 to 2011 in the observatory-frame frequency range 80-253 GHz, corresponding to a rest-frame frequency range 88-705 GHz. We find average degrees of polarization m_L ~ 2-7%; this indicates that the polarization signals are effectively averaged out by the emitter geometries. We see indication for fairly strong shocks and/or complex, variable emission region geometries in our sources, with compression factors <0.9 and/or changes in viewing angles by >10 deg. An analysis of correlations between source fluxes and polarization parameter points out special cases: the presence of (at least) two distinct emission regions with different levels of polarization (for 0415+379) as well as emission from a single, predominant component (for 0507+179 and 1418+546). Regarding the evolution of flux and polarization, we find good agreement between observations and the signal predicted by "oblique shock in jet" scenarios in one source (1418+546). We attempt to derive rotation measures for all sources, leading to actual measurements for two AGN and upper limits for three sources. We derive values of RM = -39,000 +/- 1,000 (stat) +/- 13,000 (sys) rad/m^2 and RM = 420,000 +/- 10,000 (stat) +/- 110,000 (sys) rad/m^2 for 1418+546 and 1637+574, respectively; these are the highest values reported to date for AGN. These values indicate magnetic field strengths of the order ~0.0001 G. For 0415+379, 0507+179, and 0954+658 we derive upper limits |RM| < 17,000 rad/m^2. From the relation |RM| ~ nu^a we find a = 1.9 +/- 0.3 for 1418+546, in good agreement with a = 2 as expected for a spherical or conical outflow.

We analyze the long-term evolution of the fluxes of six active galactic nuclei (AGN) - 0923+392, 3C 111, 3C 273, 3C 345, 3C 454.3, and 3C 84 - in the frequency range 80-267 GHz using archival calibration data of the IRAM Plateau de Bure Interferometer. Our dataset spans a long timeline of ~14 years with 974-3027 flux measurements per source. We find strong (factors ~2-8) flux variability on timescales of years for all sources. The flux density distributions of five out of six sources show clear signatures of bi- or even multimodality. Our sources show mostly steep (alpha~0.5-1), variable spectral indices that indicate outflow dominated emission; the variability is most probably due to optical depth variations. The power spectra globally correspond to red-noise spectra with five sources being located between the cases of white and flicker noise and one source (3C 111) being closer to the case of random walk noise. For three sources the low-frequency ends of their power spectra appear to be upscaled in spectral power by factors ~2-3 with respect to the overall powerlaws. In two sources, 3C 454.3 and 3C 84, the 1.3-mm emission preceeds the 3-mm emission by ~55 and ~300 days, respectively, probably due to (combinations of) optical depth and emission region geometry effects. We conclude that the source emission cannot be described by uniform stochastic emission processes; instead, a distinction of "quiescent" and (maybe multiple) "flare" states of the source emission appears to be necessary.

In this paper we examine properties of the variable source Sgr A* in the near-infrared (NIR) using a very extensive Ks-band data set from NACO/VLT observations taken 2004 to 2009. We investigate the variability of Sgr A* with two different photometric methods and analyze its flux distribution. We find Sgr A* is continuously emitting and continuously variable in the near-infrared, with some variability occurring on timescales as long as weeks. The flux distribution can be described by a lognormal distribution at low intrinsic fluxes (<~5 mJy, dereddened with A_{Ks}=2.5). The lognormal distribution has a median flux of approximately 1.1 mJy, but above 5 mJy the flux distribution is significantly flatter (high flux events are more common) than expected for the extrapolation of the lognormal distribution to high fluxes. We make a general identification of the low level emission above 5 mJy as flaring emission and of the low level emission as the quiescent state. We also report here the brightest Ks-band flare ever observed (from August 5th, 2008) which reached an intrinsic Ks-band flux of 27.5 mJy (m_{Ks}=13.5). This flare was a factor 27 increase over the median flux of Sgr A*, close to double the brightness of the star S2, and 40% brighter than the next brightest flare ever observed from Sgr~A*.

We present a detailed analysis of high resolution near-infrared imaging and spectroscopy of the potential star cluster IRS13E very close to the massive black hole in the Galactic Center. We detect 19 objects in IRS13E from Ks-band images, 15 of which are also detected reliably in H-band. We derive consistent proper motions for these objects from the two bands. Most objects share a similar westward proper motion. We characterize the objects using spectroscopy (1.45 to 2.45 micrometer) and (narrow-band) imaging from H- (1.66 mircrometer) to L'-band (3.80 micrometer). Nine of the objects detected in both Ks- and H-band are very red, and we find that they are all consistent with being warm dust clumps. The dust emission may be caused by the colliding winds of the two Wolf-Rayet stars in the cluster. Three of the six detected stars do not share the motion or spectral properties of the three bright stars. This leaves only the three bright, early-type stars as potential cluster members. It is unlikely that these stars are a chance configuration. Assuming the presence of an IMBH, a mass of about 14000 solar masses follows from the velocities and positions of these three stars. However, our acceleration limits make such an IMBH nearly as unlikely as a chance occurrence of such a star association. Furthermore, there is no variable X-ray source in IRS13E despite the high density of dust and gas. Therefore, we conclude that is unlikely that IRS13E hosts a black hole massive enough to bind the three stars.

We have studied the linear polarization of 86 active galactic nuclei (AGN) in the observed frequency range 80-267 GHz (3.7-1.1mm in wavelength), corresponding to rest-frame frequencies 82-738 GHz, with the IRAM Plateau de Bure Interferometer (PdBI). The large number of measurements, 441, makes our analysis the largest polarimetric AGN survey in this frequency range to date. We extracted polarization parameters via earth rotation polarimetry with unprecedented median precisions of ~0.1% in polarization fractions and ~1.2 degrees in polarization angles. For 73 of 86 sources we detect polarization at least once. The degrees of polarization are as high as ~19%, with the median over all sources being ~4%. Source fluxes and polarizations are typically highly variable, with fractional variabilities up to ~60%. We find that BLLac sources have on average the highest level of polarization. There appears to be no correlation between degree of polarization and redshift, indicating that there has been no substantial change of polarization properties since z~2.4. Our polarization and spectral index distributions are in good agreement with results found from various samples observed at cm/radio wavelengths; thus our frequency range is likely tracing the signature of synchrotron radiation without noticeable contributions from other emission mechanisms. The "millimeter-break" located at frequencies >1 THz appears to be not detectable in the frequency range covered by our survey.

We summarize recent observations and modeling of the brightest Sgr A* flare to be observed simultaneously in (near)-infrared and X-rays to date. Trying to explain the spectral characteristics of this flare through inverse Compton mechanisms implies physical parameters that are unrealistic for Sgr A*. Instead, a "cooling break" synchrotron model provides a more feasible explanation for the X-ray emission. In a magnetic field of about 5-30 Gauss the X-ray emitting electrons cool very quickly on the typical dynamical timescale while the NIR-emitting electrons cool more slowly. This produces a spectral break in the model between NIR and X-ray wavelengths that can explain the differences in the observed spectral indices.

We present new observations of the nuclear star cluster in the central parsec of the Galaxy with the adaptive optics assisted, integral field spectrograph SINFONI on the ESO/VLT. Our work allows the spectroscopic detection of early and late type stars to m_K >= 16, more than 2 magnitudes deeper than our previous data sets. Our observations result in a total sample of 177 bona fide early-type stars. We find that most of these Wolf Rayet (WR), O- and B- stars reside in two strongly warped disks between 0.8" and 12" from SgrA*, as well as a central compact concentration (the S-star cluster) centered on SgrA*. The later type B stars (m_K>15) in the radial interval between 0.8" and 12" seem to be in a more isotropic distribution outside the disks. The observed dearth of late type stars in the central few arcseconds is puzzling, even when allowing for stellar collisions. The stellar mass function of the disk stars is extremely top heavy with a best fit power law of dN/dm ~ m^(-0.45+/-0.3). Since at least the WR/O-stars were formed in situ in a single star formation event ~6 Myrs ago, this mass function probably reflects the initial mass function (IMF). The mass functions of the S-stars inside 0.8" and of the early-type stars at distances beyond 12" are compatible with a standard Salpeter/Kroupa IMF (best fit power law of dN/dm ~ m^(-2.15+/-0.3)).

In this article we identify and discuss various statistical and systematic effects influencing the astrometric accuracy achievable with MICADO, the near-infrared imaging camera proposed for the 42-metre European Extremely Large Telescope (E-ELT). These effects are instrumental (e.g. geometric distortion), atmospheric (e.g. chromatic differential refraction), and astronomical (reference source selection). We find that there are several phenomena having impact on ~100 micro-arcsec scales, meaning they can be substantially larger than the theoretical statistical astrometric accuracy of an optical/NIR 42m-telescope. Depending on type, these effects need to be controlled via dedicated instrumental design properties or via dedicated calibration procedures. We conclude that if this is done properly, astrometric accuracies of 40 micro-arcsec or better - with 40 micro-arcsec/year in proper motions corresponding to ~20 km/s at 100 kpc distance - can be achieved in one epoch of actual observations

We systematically investigate the error sources for high-precision astrometry from adaptive optics based near-infrared imaging data. We focus on the application in the crowded stellar field in the Galactic Center. We show that at the level of <=100 micro-arcseconds a number of effects are limiting the accuracy. Most important are the imperfectly subtracted seeing halos of neighboring stars, residual image distortions and unrecognized confusion of the target source with fainter sources in the background. Further contributors to the error budget are the uncertainty in estimating the point spread function, the signal-to-noise ratio induced statistical uncertainty, coordinate transformation errors, the chromaticity of refraction in Earth's atmosphere, the post adaptive optics differential tilt jitter and anisoplanatism. For stars as bright as mK=14, residual image distortions limit the astrometry, for fainter stars the limitation is set by the seeing halos of the surrounding stars. In order to improve the astrometry substantially at the current generation of telescopes, an adaptive optics system with high performance and weak seeing halos over a relatively small field (r<=3") is suited best. Furthermore, techniques to estimate or reconstruct the seeing halo could be promising.

This paper reports measurements of Sgr A* made with NACO in L' -band (3.80 um), Ks-band (2.12 um) and H-band (1.66 um) and with VISIR in N-band (11.88 um) at the ESO VLT, as well as with XMM-Newton at X-ray (2-10 keV) wavelengths. On 4 April, 2007, a very bright flare was observed from Sgr A* simultaneously at L'-band and X-ray wavelengths. No emission was detected using VISIR. The resulting SED has a blue slope (beta > 0 for nuL_nu ~ nu^beta, consistent with nuL_nu ~ nu^0.4) between 12 micron and 3.8 micron. For the first time our high quality data allow a detailed comparison of infrared and X-ray light curves with a resolution of a few minutes. The IR and X-ray flares are simultaneous to within 3 minutes. However the IR flare lasts significantly longer than the X-ray flare (both before and after the X-ray peak) and prominent substructures in the 3.8 micron light curve are clearly not seen in the X-ray data. From the shortest timescale variations in the L'-band lightcurve we find that the flaring region must be no more than 1.2 R_S in size. The high X-ray to infrared flux ratio, blue nuL_nu slope MIR to L' -band, and the soft nuL_nu spectral index of the X-ray flare together place strong constraints on possible flare emission mechanisms. We find that it is quantitatively difficult to explain this bright X-ray flare with inverse Compton processes. A synchrotron emission scenario from an electron distribution with a cooling break is a more viable scenario.

The central parsec around the super-massive black hole in the Galactic Center hosts more than 100 young and massive stars. Outside the central cusp (R~1") the majority of these O and Wolf-Rayet (WR) stars reside in a main clockwise system, plus a second, less prominent disk or streamer system at large angles with respect to the main system. Here we present the results from new observations of the Galactic Center with the AO-assisted near-infrared imager NACO and the integral field spectrograph SINFONI on the ESO/VLT. These include the detection of 27 new reliably measured WR/O stars in the central 12" and improved measurements of 63 previously detected stars, with proper motion uncertainties reduced by a factor of four compared to our earlier work. We develop a detailed statistical analysis of their orbital properties and orientations. Half of the WR/O stars are compatible with being members of a clockwise rotating system. The rotation axis of this system shows a strong transition as a function of the projected distance from SgrA*. The main clockwise system either is either a strongly warped single disk with a thickness of about 10 degrees, or consists of a series of streamers with significant radial variation in their orbital planes. 11 out of 61 clockwise moving stars have an angular separation of more than 30 degrees from the clockwise system. The mean eccentricity of the clockwise system is 0.36+/-0.06. The distribution of the counter-clockwise WR/O star is not isotropic at the 98% confidence level. It is compatible with a coherent structure such as stellar filaments, streams, small clusters or possibly a disk in a dissolving state. The observed disk warp and the steep surface density distribution favor in situ star formation in gaseous accretion disks as the origin of the young stars.

Sagittarius A*, the supermassive compact object at the center of the Galaxy, exhibits outbursts in the near infrared and X-ray domains. These flares are likely due to energetic events very close to the central object, on a scale of a few Schwarzschild radii. Optical interferometry will soon be able to provide astrometry with an accuracy of this order (~10 muas). In this article we use recent photometric near infrared data observed with the adaptive optics system NACO at the Very Large Telescope combined with simulations in order to deploy a method to test the nature of the flares and to predict the possible outcome of observations with the Very Large Telescope Interferometer. To accomplish this we implement a hot spot model and investigate its appearance for a remote observer in terms of light curves and centroid tracks, based on general relativistic ray tracing methods. First, we use a simplified model of a small steady source in order to investigate the relativistic effects qualitatively. A more realistic scenario is then being developed by fitting our model to existing flare data. While indications for the spin of the black hole and multiple images due to lensing effects are marginal in the light curves, astrometric measurements offer the possibility to reveal these high-order general relativistic effects. This study makes predictions on these astrometric measurements and leads us to the conclusion that future infrared interferometers will be able to detect proper motion of hot spots in the vicinity of Sagittarius A*.

We present the results of 16 years of monitoring stellar orbits around the massive black hole in center of the Milky Way using high resolution NIR techniques. This work refines our previous analysis mainly by greatly improving the definition of the coordinate system, which reaches a long-term astrometric accuracy of 300 microarcsecond, and by investigating in detail the individual systematic error contributions. The combination of a long time baseline and the excellent astrometric accuracy of adaptive optics data allow us to determine orbits of 28 stars, including the star S2, which has completed a full revolution since our monitoring began. Our main results are: all stellar orbits are fit extremely well by a single point mass potential to within the astrometric uncertainties, which are now 6 times better than in previous studies. The central object mass is (4.31 +- 0.06|stat +- 0.36|R0) * 10^6 M_sun where the fractional statistical error of 1.5 percent is nearly independent from R0 and the main uncertainty is due to the uncertainty in R0. Our current best estimate for the distance to the Galactic Center is R0 = 8.33 +- 0.35 kpc. The dominant errors in this value is systematic. The mass scales with distance as (3.95 +- 0.06) * 10^6 M_sun * (R0/8kpc)^2.19. The orientations of orbital angular momenta for stars in the central arcsecond are random. We identify six of the stars with orbital solutions as late type stars, and six early-type stars as members of the clockwise rotating disk system, as was previously proposed. We constrain the extended dark mass enclosed between the pericenter and apocenter of S2 at less than 0.066, at the 99% confidence level, of the mass of Sgr A*. This is two orders of magnitudes larger than what one would expect from other theoretical and observational estimates.

We aim at a detailed description of the kinematic properties of the old, (several Gyrs) late-type CO-absorption star population among the Galactic centre (GC) cluster stars. We applied AO-assisted, near-infrared imaging and integral-field spectroscopy using the instruments NAOS/CONICA and SINFONI at the VLT. We obtained proper motions for 5445 stars, 3D velocities for 664 stars, and acceleration limits (in the sky plane) for 750 stars. We detect for the first time significant cluster rotation in the sense of the general Galactic rotation in proper motions. Out of the 3D velocity dispersion, we derive an improved statistical parallax for the GC of R0 = 8.07 +/- 0.32 (stat) +/- 0.13 (sys) kpc. The distribution of 3D stellar speeds can be approximated by local Maxwellian distributions. Kinematic modelling provides deprojected 3D kinematic parameters, including the mass profile of the cluster. We find an upper limit of 4% for the amplitude of fluctuations in the phase-space distribution of the cluster stars compared to a uniform, spherical model cluster. Using upper limits on accelerations, we constrain the minimum line-of-sight distances from the plane of SgrA* of five stars located within the innermost few (projected) arcsec. The stars within 0.7'' radius from the star group IRS13E do not co-move with this group, making it unlikely that IRS13E is the core of a substantial star cluster. Overall, the GC late-type cluster is described well as a uniform, isotropic, rotating, dynamically relaxed, phase-mixed system.

We analyze the properties of the star S2 orbiting the supermassive black hole at the center of the Galaxy. A high quality SINFONI H and K band spectrum obtained from coadding 23.5 hours of observation between 2004 and 2007 reveals that S2 is an early B dwarf (B0-2.5V). Using model atmospheres, we constrain its stellar and wind properties. We show that S2 is a genuine massive star, and not the core of a stripped giant star as sometimes speculated to resolve the problem of star formation so close to the supermassive black hole. We give an upper limit on its mass loss rate, and show that it is He enriched, possibly because of the presence of a magnetic field.

We classify 329 late-type giants within 1 parsec of Sgr A*, using the adaptive optics integral field spectrometer SINFONI on the VLT. These observations represent the deepest spectroscopic data set so far obtained for the Galactic Center, reaching a 50% completeness threshold at the approximate magnitude of the helium-burning red clump (Ks ~ 15.5 mag.). Combining our spectroscopic results with NaCo H and Ks photometry, we construct an observed Hertzsprung-Russell diagram, which we quantitatively compare to theoretical distributions of various star formation histories of the inner Galaxy, using a chi-squared analysis. Our best-fit model corresponds to continuous star formation over the last 12 Gyr with a top-heavy initial mass function (IMF). The similarity of this IMF to the IMF observed for the most recent epoch of star formation is intriguing and perhaps suggests a connection between recent star formation and the stars formed throughout the history of the Galactic Center.

In the first two papers of this series, we determined the position of Sgr A* on infrared images, by aligning the positions of red giant stars with positions measured at radio wavelengths for their circumstellar SiO masers. In this paper, we report detections of 5 new stellar SiO masers within 50" (2 pc) of Sgr A* and new and/or improved positions and proper motions of 15 stellar SiO masers. The current accuracies are ~1 mas in position and ~0.3 mas/y in proper motion. We find that the proper motion of the central stellar cluster with respect to Sgr A* is less than 45 km/s. One star, IRS 9, has a three-dimensional speed of ~370 km/s at a projected distance of 0.33 pc from Sgr A*. If IRS 9 is bound to the inner parsec, this requires an enclosed mass that exceeds current estimates of the sum of the mass of Sgr A* and luminous stars in the stellar cusp by ~0.8 x 10^6 Msun. Possible explanations include i) that IRS 9 is not bound to the central parsec and has "fallen" from a radius greater than 9 pc, ii) that a cluster of dark stellar remnants accounts for some of the excess mass, and/or iii) that Ro is considerably greater than 8 kpc.

We study the stellar and wind properties of massive stars in the central cluster of the Galaxy. We use non-LTE atmosphere models including winds and line-blanketing to fit their H and K band spectra obtained with the 3D spectrograph SINFONI on the VLT. We derive the main stellar (Teff, L, abundances, ionizing flux) and wind (mass loss rate, terminal velocity) properties. They are found to be similar to other galactic massive stars. We show that a direct evolutionary link between Ofpe/WN9, WN8 and WN/C stars exists. Using individual SEDs for each massive star, we construct the total spectral energy distribution of the cluster and use it to compute photoionization models. We show that the nebular properties of the central HII region are well reproduced. We conclude that, contrary to previous claims, standard stellar evolution and atmosphere models are well suited to explain the properties of the central cluster. Our results indicate that massive stars in the central cluster do not have a peculiar evolution as could be expected from their proximity to the supermassive black hole SgrA*.

In this article we summarise and discuss the infrared, radio, and X-ray emission from the supermassive black hole in the Galactic Centre, SgrA*. We include new results from near-infrared polarimetric imaging observations obtained on May 31st, 2006. In that night, a strong flare in Ks band (2.08 microns) reaching top fluxes of ~16 mJy could be observed. This flare was highly polarised (up to ~40%) and showed clear sub-structure on a time scale of 15 minutes, including a swing in the polarisation angle of about 70 degrees. For the first time we were able to observe both polarised flux and short-time variability, with high significance in the same flare event. This result adds decisive information to the puzzle of the SgrA* activity. The observed polarisation angle during the flare peak is the same as observed in two events in 2004 and 2005. Our observations strongly support the dynamical emission model of a decaying plasma hotspot orbiting SgrA* on a relativistic orbit. The observed polarisation parameters and their variability with time might allow to constrain the orientation of accretion disc and spin axis with respect to the Galaxy.

Context: We report new polarization measurements of the variable near-infrared emission of the SgrA* counterpart associated with the massive 3--4 10^6 solar masses Black Hole at the Galactic Center. Aims: We investigate the physical processes responsible for the variable emission from SgrA*. Methods: The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope Results: We find that the variable NIR emission of SgrA* is highly polarized and consists of a contribution of a non- or weakly polarized main flare with highly polarized sub-flares. The flare activity shows a possible quasi-periodicity of 20+-3 minutes consistent with previous observations. Conclusions: The highly variable and polarized emission supports that the NIR emission is non-thermal. The observations can be interpreted in a jet or temporary disk model. In the disk model the quasi-periodic flux density variations can be explained due to spots on relativistic orbits around the central MBH. Alternative explanations for the high central mass concentration involving boson or fermion balls are increasingly unlikely.

We report on the spectroscopic monitoring of GCIRS16SW, an Ofpe/WN9 star and LBV candidate in the central parsec of the Galaxy. SINFONI observations show strong daily spectroscopic changes in the K band. Radial velocities are derived from the HeI 2.112 um line complex and vary regularly with a period of 19.45 days, indicating that the star is most likely an eclipsing binary. Under various assumptions, we are able to derive a mass of ~ 50 Msun for each component.

We have observed a bright flare of Sgr A* in the near infrared with the adaptive optics assisted integral field spectrometer SINFON. Within the uncertainties, the observed spectrum is featureless and can be described by a power law. Our data suggest that the spectral index is correlated with the instantaneous flux and that both quantities experience significant changes within less than one hour. We argue that the near infrared flares from Sgr A* are due to synchrotron emission of transiently heated electrons, the emission being affected by orbital dynamics and synchrotron cooling, both acting on timescales of ~20 minutes.

We report the definite spectroscopic identification of 41 OB supergiants, giants and main sequence stars in the central parsec of the Galaxy. Detection of their absorption lines have become possible with the high spatial and spectral resolution and sensitivity of the adaptive optics integral field spectrometer SPIFFI/SINFONI on the ESO VLT. Several of these OB stars appear to be helium and nitrogen rich. Almost all of the ~80 massive stars now known in the central parsec (central arcsecond excluded) reside in one of two somewhat thick (<|h|/R>~0.14) rotating disks. These stellar disks have fairly sharp inner edges (R~1") and surface density profiles that scale as R^{-2}. We do not detect any OB stars outside the central 0.5 pc. The majority of the stars in the clockwise system appear to be on almost circular orbits, whereas most of those in the `counter-clockwise' disk appear to be on eccentric orbits. Based on its stellar surface density distribution and dynamics we propose that IRS 13E is an extremely dense cluster (core density > 3x10^8 sunmass/pc^3), which has formed in the counter-clockwise disk. The stellar contents of both systems are remarkably similar, indicating a common age of ~6+/-2 Myr. The K-band luminosity function of the massive stars suggests a top-heavy mass function and limits the total stellar mass contained in both disks to ~1.5x10^4 sunmass. Our data strongly favor in situ star formation from dense gas accretion disks for the two stellar disks. This conclusion is very clear for the clockwise disk and highly plausible for the counter-clockwise system.

We report new simultaneous near-infrared/sub-millimeter/X-ray observations of the SgrA* counterpart associated with the massive 3-4x10**6 solar mass black hole at the Galactic Center. The main aim is to investigate the physical processes responsible for the variable emission from SgrA*. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope and the ACIS-I instrument aboard the Chandra X-ray Observatory as well as the Submillimeter Array SMA on Mauna Kea, Hawaii, and the Very Large Array in New Mexico. We detected one moderately bright flare event in the X-ray domain and 5 events at infrared wavelengths.

We report the results of time-resolved photometric and spectroscopic near-infrared observations of the Ofpe/WN9 star and LBV candidate GCIRS34W in the Galactic Centre star cluster. Diffraction limited resolution photometric observations obtained in H and K bands show a strong, non-periodic variability on time scales from months to years in both bands accompanied by variations of the stellar colour. Three K band spectra obtained in 1996, 2003 and 2004 with integral field spectrometers are identical within their accuracies and exclude significant spectroscopic variability. The most probable explanation of the stellar photometric variability is obscuration by circumstellar material ejected by the star. The approximated position of GCIRS34W in a HR diagram is located between O supergiants and LBVs, suggesting that this star is a transitional object between these two phases of stellar evolution.

We report 75 milli-arcsec resolution, near-IR imaging spectroscopy within the central 30 light days of the Galactic Center [...]. To a limiting magnitude of K~16, 9 of 10 stars in the central 0.4 arcsec, and 13 of 17 stars out to 0.7 arcsec from the central black hole have spectral properties of B0-B9, main sequence stars. [...] all brighter early type stars have normal rotation velocities, similar to solar neighborhood stars. We [...] derive improved 3d stellar orbits for six of these S-stars in the central 0.5 arcsec. Their orientations in space appear random. Their orbital planes are not co-aligned with those of the two disks of massive young stars 1-10 arcsec from SgrA*. We can thus exclude [...] that the S-stars as a group inhabit the inner regions of these disks. They also cannot have been located/formed in these disks [...]. [...] we conclude that the S-stars were most likely brought into the central light month by strong individual scattering events. The updated estimate of distance to the Galactic center from the S2 orbit fit is Ro = 7.62 +/- 0.32 kpc, resulting in a central mass value of 3.61 +/- 0.32 x 10^6 Msun. We happened to catch two smaller flaring events from SgrA* [...]. The 1.7-2.45 mum spectral energy distributions of these flares are fit by a featureless, red power law [...]. The observed spectral slope is in good agreement with synchrotron models in which the infrared emission comes from [...] radiative inefficient accretion flow in the central R~10 Rs region.