Theoretically, the running of the cosmological constant in the IR region is not ruled out. On the other hand, from the QFT viewpoint, the energy released due to the variation of the cosmological constant in the late universe cannot go to the matter sector. For this reason, the phenomenological bounds on such a running are not sufficiently restrictive. The situation can be different in the early universe when the gravitational field was sufficiently strong to provide an efficient creation of particles from the vacuum. We develop a framework for systematically exploring this ossibility. It is supposed that the running occurs in the epoch when the Dark Matter already decoupled and is expanding adiabatically, while baryons are approximately massless and can be abundantly created from vacuum due to the decay of vacuum energy. By using the handy model of Reduced Relativistic Gas for describing the Dark Matter, we consider the dynamics of both cosmic background and linear perturbations and evaluate the impact of the vacuum decay on the matter power spectrum and to the first CMB peak. Additionally, using the combined data of CMB+BAO+SNIa we find the best fit values for the free parameters of our model.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

Reduced Relativistic Gas (RRG) is a useful approach to describe the warm dark matter (WDM) or the warmness of baryonic matter in the approximation when the interaction between the particles is irrelevant. The use of Maxwell distribution leads to the complicated equation of state of the J\"{u}ttner model of relativistic ideal gas. The RRG enables one to reproduce the same physical situation but in a much simpler form. For this reason RRG can be a useful tool for the theories with some sort of a "new Physics". On the other hand, even without the qualitatively new physical implementations, the RRG can be useful to describe the general features of WDM in a model-independent way. In this sense one can see, in particular, to which extent the cosmological manifestations of WDM may be dependent on its Particle Physics background. In the present work RRG is used as a complementary approach to derive the main observational exponents for the WDM in a model-independent way. The only assumption concerns a non-negligible velocity $v$ for dark matter particles which is parameterized by the warmness parameter $b$. The relatively high values of $b$ ( $b^2\gtrsim 10^{-6}$) erase the radiation (photons and neutrinos) dominated epoch and cause an early warm matter domination after inflation. Furthermore, RRG approach enables one to quantify the lack of power in linear matter spectrum at small scales and in particular, reproduces the relative transfer function commonly used in context of WDM with accuracy of $\lesssim 1\%$. A warmness with $b^2\lesssim 10^{-6}$ (equivalent to $v\lesssim 300 km/s$) does not alter significantly the CMB power spectrum and is in agreement with the background observational tests.

Quantum effects derived through conformal anomaly lead to an inflationary model that can be either stable or unstable. The unstable version requires a large dimensionless coefficient of about $5\times 10^8$ in front of the $R^2$ term that results in the inflationary regime in the $R+R^2$ ("Starobinsky") model being a generic intermediate attractor. In this case the non-local terms in the effective action are practically irrelevant, and there is a 'graceful exit' to a low curvature matter-like dominated stage driven by high-frequency oscillations of $R$ -- scalarons, which later decay to pairs of all particles and antiparticles, with the amount of primordial scalar (density) perturbations required by observations. The stable version is a genuine generic attractor, so there is no exit from it. We discuss a possible transition from stable to unstable phases of inflation. It is shown that this transition is automatic if the sharp cut-off approximation is assumed for quantum corrections in the period of transition. Furthermore, we describe two different quantum mechanisms that may provide a required large $\,R^2$-term in the transition period.

The Earth-Moon-Sun system has traditionally provided the best laboratory for testing the strong equivalence principle. For a decade, the Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) has been producing the world's best lunar laser ranging data. At present, a single observing session of about an hour yields a distance measurement with uncertainty of about 2~mm, an order of magnitude advance over the best pre-APOLLO lunar laser ranging data. However, these superb data have not yet yielded scientific results commensurate with their accuracy, number, and temporal distribution. There are two reasons for this. First, even in the relatively clean environment of the Earth-Moon system, a large number of effects modify the measured distance importantly and thus need to be included in the analysis model. The second reason is more complicated. The traditional problem with the analysis of solar-system metric data is that the physical model must be truncated to avoid extra parameters that would increase the condition number of the estimator. Even in a typical APOLLO analysis that does not include parameters of gravity physics, the condition number is very high: $8 \times 10^{10}$.

We review the radio very long baseline interferometry (VLBI) observations of the guide star, IM Peg, and three compact extragalactic reference sources, 3C 454.3, B2250+194, and B2252+172, made in support of the NASA/Stanford gyroscope relativity mission, GP-B. The main goal of the observations was the determination of the proper motion of IM Peg relative to the distant universe. VLBI observations made between 1997 and 2005 yield a proper motion of IM Peg of -20.83 $\pm$ 0.09 mas yr$^{-1}$ in RA and -27.27 $\pm$ 0.09 mas yr$^{-1}$ in dec, in a celestial reference frame of extragalactic radio galaxies and quasars virtually identical to the International Celestial Reference Frame 2 (ICRF2). They also yield a parallax for IM Peg of 10.37 $\pm$ 0.07 mas, corresponding to a distance of 96.4 $\pm$ 0.7 pc. The uncertainties are standard errors with statistical and estimated systematic contributions added in quadrature. These results met the pre-launch requirements of the GP-B mission to not discernibly degrade the estimates of the geodetic and frame-dragging effects. The paper also reports on a 1$\sigma$ upper limit on the magnitude of the components of the proper motion of the 'core' of 3C 454.3 relative to the ICRF2 of 46 and 56 $\mu$as yr$^{-1}$ in RA and dec, respectively, and presents densely sampled ellipses of the parallax and the orbit of the giant of the binary system. It further gives a sequence of images of the flickering radio emission relative to the disk of the giant. For a 'movie of a star,' see http://www.yorku.ca/bartel/impeg.mpg .

The Reduced Relativistic Gas (RRG) model was introduced by A. Sakharov in 1965 for deriving the cosmic microwave background (CMB) spectrum. It was recently reinvented by some of us to achieve an interpolation between the radiation and dust epochs in the evolution of the Universe. This model circumvents the complicated structure of the Boltzmann-Einstein system of equations and admits a transparent description of warm-dark-matter effects. It is extended here to include, on a phenomenological basis, an out-of-equilibrium interaction between radiation and baryons which is supposed to account for relevant aspects of pre-recombination physics in a simplified manner. Furthermore, we use the tight-coupling approximation to explore the influence of both this interaction and of the RRG warmness parameter on the anisotropy spectrum of the CMB. The predictions of the model are very similar to those of the {\Lambda}CDM model if both the interaction and the dark-matter warmness parameters are of the order of $10^{-4}$ or smaller. As far as the warmness parameter is concerned, this is in good agreement with previous estimations on the basis of results from structure formation.

We explore the phenomenology of nontrivial quantum effects on low-energy gravity. These effects come from the running of the gravitational coupling parameter G and the cosmological constant L in the Einstein-Hilbert action, as induced by the Renormalization Group (RG). The Renormalization Group corrected General Relativity (RGGR model) is used to parametrize these quantum effects, and it is assumed that the dominant dark matter-like effects inside galaxies is due to these nontrivial RG effects. Here we present additional details on the RGGR model application, in particular on the Poisson equation extension that defines the effective potential, also we re-analyse the ordinary elliptical galaxy NGC 4494 using a slightly different model for its baryonic contribution, and explicit solutions are presented for the running of G and L. The values of the NGC 4494 parameters as shown here have a better agreement with the general RGGR picture for galaxies, and suggest a larger radial anisotropy than the previously published result.

The paper is about possible effects of infrared quantum contributions to General Relativity on disk and elliptical galaxies. The Renormalization Group corrected General Relativity (RGGR model) is used to parametrize these quantum effects. The new RGGR results presented here concern the elliptical galaxy NGC 4374 and the dwarf disk galaxy DDO 47. Using the effective approach to Quantum Field Theory in curved background, one can argue that the proper RG energy scale, in the weak field limit, should be related to the Newtonian potential. In the context of galaxies, this led to a remarkably small variation of the gravitational coupling G, while also capable of generating galaxy rotation and dispersion curves of similar quality to the the best dark matter profiles (i.e., the profiles that have a core).

When VLBI observations are used to determine the position or motion of a radio source relative to reference sources nearby on the sky, the astrometric information is usually obtained via: (i) phase-referenced maps; or (ii) parametric model fits to measured fringe phases or multiband delays. In this paper we describe a "merged" analysis technique which combines some of the most important advantages of these other two approaches. In particular, our merged technique combines the superior model-correction capabilities of parametric model fits with the ability of phase-referenced maps to yield astrometric measurements of sources that are too weak to be used in parametric model fits. We compare the results from this merged technique with the results from phase-referenced maps and from parametric model fits in the analysis of astrometric VLBI observations of the radio-bright star IM Pegasi (HR 8703) and the radio source B2252+172 nearby on the sky. In these studies we use central-core components of radio sources 3C 454.3 and B2250+194 as our positional references. We obtain astrometric results for IM Peg with our merged technique even when the source is too weak to be used in parametric model fits, and we find that our merged technique yields superior astrometric results to the phase-referenced mapping technique. We used our merged technique to estimate the proper motion and other astrometric parameters of IM Peg in support of the NASA/Stanford Gravity Probe B mission.

We present measurements of the total radio flux density as well as very-long-baseline interferometry (VLBI) images of the star, IM Pegasi, which was used as the guide star for the NASA/Stanford relativity mission Gravity Probe B. We obtained flux densities and images from 35 sessions of observations at 8.4 GHz (wavelength = 3.6 cm) between 1997 January and 2005 July. The observations were accurately phase-referenced to several extragalactic reference sources, and we present the images in a star-centered frame, aligned by the position of the star as derived from our fits to its orbital motion, parallax, and proper motion. Both the flux density and the morphology of IM Peg are variable. For most sessions, the emission region has a single-peaked structure, but 25% of the time, we observed a two-peaked (and on one occasion perhaps a three-peaked) structure. On average, the emission region is elongated by 1.4 +- 0.4 mas (FWHM), with the average direction of elongation being close to that of the sky projection of the orbit normal. The average length of the emission region is approximately equal to the diameter of the primary star. No significant correlation with the orbital phase is found for either the flux density or the direction of elongation, and no preference for any particular longitude on the star is shown by the emission region.

We present a physical interpretation for the locations of the sources of radio emission in IM Pegasi (IM Peg, HR 8703), the guide star for the NASA/Stanford relativity mission Gravity Probe B. This emission is seen in each of our 35 epochs of 8.4-GHz VLBI observations taken from 1997 to 2005. We found that the mean position of the radio emission is at or near the projected center of the primary to within about 27% of its radius, identifying this active star as the radio emitter. The positions of the radio brightness peaks are scattered across the disk of the primary and slightly beyond, preferentially along an axis with position angle, p.a. = (-38 +- 8) deg, which is closely aligned with the sky projections of the orbit normal (p.a. = -49.5 +- 8.6 deg) and the expected spin axis of the primary. Comparison with simulations suggests that brightness peaks are 3.6 (+0.4,-0.7) times more likely to occur (per unit surface area) near the pole regions of the primary (|latitude| >= 70 deg) than near the equator (|latitude| <= 20 deg), and to also occur close to the surface with ~2/3 of them at altitudes not higher than 25% of the radius of the primary.

We used 8.4 GHz VLBI images obtained at up to 35 epochs between 1997 and 2005 to examine the radio structures of the main reference source, 3C 454.3, and two secondary reference sources, B2250+194 and B2252+172, for the guide star for the NASA/Stanford relativity mission Gravity Probe B (GP-B). For one epoch in 2004 May, we also obtained images at 5.0 and 15.4 GHz. The 35 8.4 GHz images for quasar 3C 454.3 confirm a complex, evolving, core-jet structure. We identified at each epoch a component, C1, near the easternmost edge of the core region. Simulations of the core region showed that C1 is located, on average, 0.18 +- 0.06 mas west of the unresolved "core" identified in 43 GHz images. We also identified in 3C 454.3 at 8.4 GHz several additional components which moved away from C1 with proper motions ranging in magnitude between 0.9c and 5c. The detailed motions of the components exhibit two distinct bends in the jet axis located ~3 and ~5.5 mas west of C1. The spectra between 5.0 and 15.4 GHz for the "moving" components are steeper than that for C1. The 8.4 GHz images of B2250+194 and B2252+172, in contrast to those of 3C 454.3, reveal compact structures. The spectrum between 5.0 and 15.4 GHz for B2250+194 is inverted while that for B2252+172 is flat. Based on its position near the easternmost edge of the 8.4 GHz radio structure, close spatial association with the 43 GHz core, and relatively flat spectrum, we believe 3C 454.3 component C1 to be the best choice for the ultimate reference point for the GP-B guide star. The compact structures and inverted to flat spectra of B2250+194 and B2252+172 make these objects valuable secondary reference sources

We present the principal astrometric results of the very-long-baseline interferometry (VLBI) program undertaken in support of the Gravity Probe B (GP-B) relativity mission. VLBI observations of the GP-B guide star, the RS CVn binary IM Pegasi (HR 8703), yielded positions at 35 epochs between 1997 and 2005. We discuss the statistical assumptions behind these results and our methods for estimating the systematic errors. We find the proper motion of IM Peg in an extragalactic reference frame closely related to the International Celestial Reference Frame 2 (ICRF2) to be -20.83 +- 0.03 +- 0.09 mas/yr in right ascension and -27.27 +- 0.03 +- 0.09 mas/yr in declination. For each component the first uncertainty is the statistical standard error and the second is the total standard error (SE) including plausible systematic errors. We also obtain a parallax of 10.37 +- 0.07 mas (distance: 96.4 +- 0.7 pc), for which there is no evidence of any significant contribution of systematic error. Our parameter estimates for the ~25-day-period orbital motion of the stellar radio emission have SEs corresponding to ~0.10 mas on the sky in each coordinate. The total SE of our estimate of IM Peg's proper motion is ~30% smaller than the accuracy goal set by the GP-B project before launch: 0.14 mas/yr for each coordinate of IM Peg's proper motion. Our results ensure that the uncertainty in IM Peg's proper motion makes only a very small contribution to the uncertainty of the GP-B relativity tests.

We made VLBI observations at 8.4 GHz between 1997 and 2005 to estimate the coordinates of the "core" component of the superluminal quasar, 3C 454.3, the ultimate reference point in the distant universe for the NASA/Stanford Gyroscope Relativity Mission, Gravity Probe B. These coordinates are determined relative to those of the brightness peaks of two other compact extragalactic sources, B2250+194 and B2252+172, nearby on the sky, and within a celestial reference frame (CRF), defined by a large suite of compact extragalactic radio sources, and nearly identical to the International Celestial Reference Frame 2 (ICRF2). We find that B2250+194 and B2252+172 are stationary relative to each other, and also in the CRF, to within 1-sigma upper limits of 15 and 30 micro-arcsec/yr in RA and decl., respectively. The core of 3C 454.3 appears to jitter in its position along the jet direction over ~0.2 mas, likely due to activity close to the putative supermassive black hole nearby, but on average is stationary in the CRF within 1-sigma upper limits on its proper motion of 39 micro-arcsec/yr (1.0c) and 30 micro-arcsec/yr (0.8c) in RA and decl., respectively, for the period 2002 - 2005. Our corresponding limit over the longer interval, 1998 - 2005, of more importance to GP-B, is 46 and 56 micro-arcsec/yr in RA and decl., respectively. Some of 3C 454.3's jet components show significantly superluminal motion with speeds of up to ~200 micro-arcsec/yr or 5c in the CRF. The core of 3C 454.3 thus provides for Gravity Probe B a sufficiently stable reference in the distant universe.

We describe the NASA/Stanford gyroscope relativity mission, Gravity Probe B (GP-B), and provide an overview of the following series of six astrometric and astrophysical papers that report on our radio observations and analyses made in support of this mission. The main goal of this 8.5 year program of differential VLBI astrometry was to determine the proper motion of the guide star of the GP-B mission, the RS CVn binary IM Pegasi (IM Peg; HR 8703). This proper motion is determined with respect to compact, extragalactic reference sources. The results are: -20.833 +- 0.090 mas/yr and -27.267 +- 0.095 mas/yr for, respectively, the right ascension and declination, in local Cartesian coordinates, of IM Peg's proper motion, and 10.370 +- 0.074 mas (i.e., 96.43 +- 0.69 pc) for its parallax (and distance). Each quoted uncertainty is meant to represent an ~70% confidence interval that includes the estimated contribution from systematic error. These results are accurate enough not to discernibly degrade the GP-B estimates of its gyroscopes' relativistic precessions: the frame-dragging and geodetic effects.

The quantum contributions to the gravitational action are relatively easy to calculate in the higher derivative sector of the theory. However, the applications to the post-inflationary cosmology and astrophysics require the corrections to the Einstein-Hilbert action and to the cosmological constant, and those we can not derive yet in a consistent and safe way. At the same time, if we assume that these quantum terms are covariant and that they have relevant magnitude, their functional form can be defined up to a single free parameter, which can be defined on the phenomenological basis. It turns out that the quantum correction may lead, in principle, to surprisingly strong and interesting effects in astrophysics and cosmology.

We use the framework of a recently proposed model of reduced relativistic gas (RRG) to obtain the bounds for $\Omega$'s of Dark Matter and Dark Energy (in the present case, a cosmological constant), taking into consideration an arbitrary warmness of Dark Matter. An equivalent equation of state has been used by Sakharov to predict the oscillations in the matter power spectrum. Two kind of tests are accounted for in what follows, namely the ones coming from the dynamics of the conformal factor of the homogeneous and isotropic metric and also the ones based on linear cosmic perturbations. The RRG model demonstrated its high effectiveness, permitting to explore a large volume in the space of mentioned parameters in a rather economic way. Taking together the results of such tests as Supernova type Ia (Union2 sample), $H(z)$, CMB ($R$ factor), BAO and LSS (2dfGRS data), we confirm that $\La$CDM is the most favored model. At the same time, for the 2dfGRS data alone we found that an alternative model with a very small quantity of a Dark Matter is also viable. This output is potentially relevant in view of the fact that the LSS is the only test which can not be affected by the possible quantum contributions to the low-energy gravitational action.

We review our contribution to infrared Renormalization Group (RG) effects to General Relativity in the context of galaxies. Considering the effective action approach to Quantum Field Theory in curved background, we argued that the proper RG energy scale, in the weak field limit, should be related to the Newtonian potential. In the galaxy context, even without dark matter, this led to a remarkably small gravitational coupling G variation (about or less than 10^{-12} of its value per light-year), while also capable of generating galaxy rotation curves about as good as the best phenomenological dark matter profiles (considering both the rotation curve shape and the expected mass-to-light ratios). Here we also comment on related developments, open issues and perspectives.

Recently it was shown that quantum corrections to the Newton potential can explain the rotation curves in spiral galaxies without introducing the Dark Matter halo. The unique phenomenological parameter $\al\nu$ of the theory grows with the mass of the galaxy. In order to better investigate the mass-dependence of $\al\nu$ one needs to check the upper bound for $\al\nu$ at a smaller scale. Here we perform the corresponding calculation by analyzing the dynamics of the Laplace-Runge-Lenz vector. The resulting limitation on quantum corrections is quite severe, suggesting a strong mass-dependence of $\al\nu$.

We consider the application of quantum corrections computed using renormalization group arguments in the astrophysical domain and show that, for the most natural interpretation of the renormalization group scale parameter, a gravitational coupling parameter $G$ varying $10^{-7}$ of its value across a galaxy (which is roughly a variation of $10^{-12}$ per light-year) is sufficient to generate galaxy rotation curves in agreement with the observations. The quality of the resulting fit is similar to the Isothermal profile quality once both the shape of the rotation curve and the mass-to-light ratios are considered for evaluation. In order to perform the analysis, we use recent high quality data from nine regular disk galaxies. For the sake of comparison, the same set of data is modeled also for the Modified Newtonian Dynamics (MOND) and for the recently proposed Scalar Tensor Vector Gravity (STVG). At face value, the model based on quantum corrections clearly leads to better fits than these two alternative theories.

Cosmologies with running cosmological term (Lambda) and gravitational Newton's coupling (G) may naturally be expected if the evolution of the universe can ultimately be derived from the first principles of Quantum Field Theory or String Theory. In this paper, we derive the general cosmological perturbation equations for models with variable G and Lambda in which the fluctuations in both variables are explicitly included. We demonstrate that, if matter is covariantly conserved, the late growth of matter density perturbations is independent of the wavenumber. Furthermore, if Lambda is negligible at high redshifts and G varies slowly, we find that these cosmologies produce a matter power spectrum with the same shape as that of the concordance LCDM model, thus predicting the same basic features on structure formation. Despite this shape indistinguishability, the free parameters of the variable G and Lambda models can still be effectively constrained from the observational bounds on the spectrum amplitude.

Despite the many outstanding cosmological observations leading to a strong evidence for a nonvanishing cosmological constant (CC) term in the gravitational field equations, the theoretical status of this quantity seems to be lagging well behind the observational successes. It thus seems timely to revisit some fundamental aspects of the CC term in Quantum Field Theory (QFT). We emphasize that, in curved space-time, nothing a priori prevents this term from potentially having a mild running behavior associated to quantum effects. Remarkably, this could be the very origin of the dynamical nature of the Dark Energy, in contrast to many other popular options considered in the literature. In discussing this possibility, we also address some recent criticisms concerning the possibility of such running. Our conclusion is that, while there is no comprehensive proof of the CC running, there is no proof of the non-running either. The problem can be solved only through a deeper understanding of the vacuum contributions of massive quantum fields on a curved spacetime background. We suggest that such investigations are at the heart of one of the most important endeavors of fundamental theoretical cosmology in the years to come.

We have studied the growth of the shell-like radio structure of supernova SN1993J in M81 from September 1993 through October 2003 with very-long-baseline interferometry (VLBI) observations at the wavelengths of 3.6, 6, and 18cm. For this purpose, we have developed a method to accurately determine the outer radius (R) of any circularly symmetric compact radio structure like SN1993J. The source structure of SN1993J remains circularly symmetric (with deviations from circularity under 2%) over almost 4000 days. We characterize the decelerated expansion of SN 1993J through approximately day 1500 after explosion with an expansion parameter $m= 0.845\pm0.005$ ($R \propto t^{m}$). However, from that day onwards the expansion is different when observed at 6 and 18cm. Indeed, at 18cm, the expansion can be well characterized by the same $m$ as before day 1500, while at 6cm the expansion appears more decelerated, and is characterized by another expansion parameter, $m_{6}= 0.788\pm0.015$. Therefore, since about day 1500 on, the radio source size has been progressively smaller at 6cm than at 18cm. These findings are in stark contrast to previous reports by other authors on the details of the expansion. In our interpretation the supernova expands with a single expansion parameter, $m= 0.845\pm0.005$, and the 6cm results beyond day 1500 are due to physical effects, perhaps also coupled to instrumental limitations. Two physical effects may be involved: (a) a changing opacity of the ejecta to the 6cm radiation, and (b) a radial decrease of the magnetic field in the emitting region. (Long abstract cut. Please, read full abstract in manuscript).

A long standing goal in astrophysics is to directly observe the immediate environment of a black hole with angular resolution comparable to the event horizon. Realizing this goal would open a new window on the study of General Relativity in the strong field regime, accretion and outflow processes at the edge of a black hole, the existence of an event horizon, and fundamental black hole physics (e.g., spin). Steady long-term progress on improving the capability of Very Long Baseline Interferometry (VLBI) at short wavelengths has now made it extremely likely that this goal will be achieved within the next decade. The most compelling evidence for this is the recent observation by 1.3mm VLBI of Schwarzschild radius scale structure in SgrA*, the compact source of radio, submm, NIR and xrays at the center of the Milky Way. SgrA* is thought to mark the position of a ~4 million solar mass black hole, and because of its proximity and estimated mass presents the largest apparent event horizon size of any black hole candidate in the Universe. Over the next decade, existing and planned mm/submm facilities will be combined into a high sensitivity, high angular resolution "Event Horizon Telescope" that will bring us as close to the edge of black hole as we will come for decades. This white paper describes the science case for mm/submm VLBI observations of both SgrA* and M87 (a radio loud AGN of a much more luminous class that SgrA*). We emphasize that while there is development and procurement involved, the technical path forward is clear, and the recent successful observations have removed much of the risk that would normally be associated with such an ambitious project.

We report on phase-referenced 23 GHz Very-Long-Baseline-Interferometry (VLBI) observations of the type IIb supernova SN 2008ax, made with the Very Long Baseline Array (VLBA) on 2 April 2008 (33 days after explosion). These observations resulted in a marginal detection of the supernova. The total flux density recovered from our VLBI image is 0.8$\pm$0.3 mJy (one standard deviation). As it appears, the structure may be interpreted as either a core-jet or a double source. However, the supernova structure could be somewhat confused with a possible close by noise peak. In such a case, the recovered flux density would decrease to 0.48$\pm$0.12 mJy, compatible with the flux densities measured with the VLA at epochs close in time to our VLBI observations. The lowest average expansion velocities derived from our observations are $(1.90 \pm 0.30) \times 10^5$ km s$^{-1}$ (case of a double source) and $(5.2 \pm 1.3) \times 10^4$ km s$^{-1}$ (taking the weaker source component as a spurious, close by, noise peak, which is the more likely interpretation). These velocities are 7.3 and 2 times higher, respectively, than the maximum ejecta velocity inferred from optical-line observations.

One of important questions concerning particles which compose the Dark Matter (DM) is their average speed. We consider the model of relativistic weakly interacting massive particles and try to impose an upper bound on their actual and past warmness through the analysis of density perturbations and comparison with the LSS data. It is assumed that the DM can be described by the recently invented model of reduced relativistic gas (RRG). The equation of state of the RRG model is closely reproducing the one of the Maxwell distribution, while being much simpler. This advantage of the RRG model makes our analysis very efficient. As a result we arrive at the rigid and model-independent bound for the DM warmness without using the standard (much more sophisticated) approach based on the Einstein-Boltzmann system of equations.

Using standard quantum field theory, we discuss several theoretical aspects of the possible running of the cosmological constant (CC) term in Einstein's equations. The basic motivation for the present work is to emphasize that this possibility should also be taken into account when considering dynamical models for the dark energy (DE), which are nowadays mainly focused on identifying the DE with the energy density associated to one or more ad hoc scalar fields. At the same time, we address some recent criticisms that have been published (or privately communicated to us) attempting to cast doubts on the fundamental possibility of such running. In this work, we argue that while there is no comprehensive proof of the CC running, there is no rigorous proof of the non-running either. In particular, some purported "non-running theorem" recently adduced in the literature is, in our opinion, completely insubstantial and formally incorrect. The way to the CC running is, therefore, still open and we take here the opportunity to present a pedagogical review of the present state of the art in this field, including a a brief historical account.

We present brief, to great extent pedagogical review on renormalization in curved space-time and of some recent results on the derivation and better understanding of quantum corrections to the action of gravity. The paper is mainly devoted to the semiclassical approach, but we also discuss its importance for quantum gravity and string theory.

In 1981, a faint radio source (G') was detected near the center of the lensing galaxy of the famous "twin quasar" Q0957+561. It is still unknown whether this central radio source is a third quasar image or an active nucleus of the lensing galaxy, or a combination of both. In an attempt to resolve this ambiguity, we observed Q0957+561 at radio wavelengths of 13cm, 18cm, and 21cm, using the Very Long Baseline Array in combination with the phased Very Large Array and the Green Bank Telescope. We measured the spectrum of G' for the first time and found it to be significantly different from the spectra of the two bright quasar images. This finding suggests that the central component is primarily or entirely emission from the foreground lens galaxy, but the spectrum is also consistent with the hypothesis of a central quasar image suffering free-free absorption. In addition, we confirm the previously-reported VLBI position of G' just north of the optical center of the lens galaxy. The position slightly favors the hypothesis that G' originates in the lens, but is not conclusive. We discuss the prospects for further clarification of this issue.

We report on 8.4GHz Very Long Baseline Interferometry (VLBI) observations of the type II-P supernova SN2004et in the spiral galaxy NGC 6946, made on 20 February 2005 (151 days after explosion). The Very Large Array (VLA) flux density was 1.23$\pm$0.07 mJy, corresponding to an isotropic luminosity at 8.4GHz of (4.45$\pm$0.3)$\times10^{25}$ erg s$^{-1}$ Hz$^{-1}$ and a brightness temperature of (1.3$\pm$0.3)$\times10^{8}$ K. We also provide an improved source position, accurate to about 0.5 mas in each coordinate. The VLBI image shows a clear asymmetry. From model fitting of the size of the radio emission, we estimate a minimum expansion velocity of 15,700$\pm$2,000 km s$^{-1}$. This velocity is more than twice the expected mean expansion velocity estimated from a synchrotron self-absorbed emission model, thus suggesting that synchrotron self-absorption is not relevant for this supernova. With the benefit of an optical spectrum obtained 12 days after explosion, we favor an emission model which consists of two hot spots on an underlying expanding shell of width comparable to that of SN 1993J.

The dynamics of density and metric perturbations is investigated for the previously developed model where the decay of the vacuum energy into matter (or vice versa) is due to the renormalization group (RG) running of the cosmological constant (CC) term. The evolution of the CC depends on the single parameter \nu, which characterizes the running of the CC produced by the quantum effects of matter fields of the unknown high energy theory below the Planck scale. The sign of \nu indicates whether bosons or fermions dominate in the running. The spectrum of perturbations is computed assuming an adiabatic regime and an isotropic stress tensor. Moreover, the perturbations of the CC term are generated from the simplest covariant form suggested by the RG model under consideration. The corresponding numerical analysis shows that for \nu>0 there is a depletion of the matter power spectrum at low scales (large wave numbers) as compared to the standard LCDM model, whereas for \nu<0 there is an excess of power at low scales. We find that the LSS data rule out the range |\nu|> 10^{-4} while the values |\nu|< 10^{-6} look perfectly acceptable. For \nu<0 the excess of power at low scales grows rapidly and the bound is more severe. From the particle physics viewpoint, the values |\nu|\sim 10^{-6} correspond to the ``desert'' in the mass spectrum above the GUT scale M_X\sim 10^{16} GeV. Our results are consistent with those obtained in other dynamical models admitting an interaction between dark matter and dark energy. We find that the matter power spectrum analysis is a highly efficient method to discover a possible scale dependence of the vacuum energy.

We report on 8.4 GHz VLBI observations of SN2001gd in the spiral galaxy NGC5033 made on 26 June 2002 and 8 April 2003. Our data nominally suggests a relatively strong deceleration for the expansion of SN2001gd, but we cannot dismiss the possibility of a free supernova expansion. From our VLBI observations on 8 April 2003, we inferred a minimum total energy in relativistic particles and magnetic fields in the supernova shell of E_min =(0.3-14) 10^{47} ergs, and a corresponding equipartition average magnetic field of B_min = (50--350) mG. We also present multiwavelength VLA measurements of SN2001gd, which are well fit by an optically thin, synchrotron spectrum, partially absorbed by thermal plasma. We obtain a supernova flux density of (1.02 +/- 0.05) mJy at the observing frequency of 8.4 GHz for the second epoch, which results in an isotropic radio luminosity of (6.0 +/- 0.3) * 10^{36} ergs between 1.4 and 43.3 GHz, at an adopted distance of 13.1 Mpc. Finally, we report on an XMM-Newton X-ray detection of SN2001gd on 18 December 2002. The supernova X-ray spectrum is consistent with optically thin emission from a soft component (associated with emission from the reverse shock) at a temperature around 1 keV. The observed flux corresponds to an isotropic X-ray luminosity of L_X = 1.4 +/- 0.4 * 10^{39} ergs/s in the (0.3-5) keV band. We suggest that both radio and X-ray observations of SN2001gd indicate that a circumstellar interaction similar to that displayed by SN1993J in M81 is taking place.

The renormalization group (RG) approach to cosmology is an efficient method to study the possible evolution of the cosmological parameters from the point of view of quantum field theory in curved space-time. In this work we continue our previous investigations of the RG method based on potential low-energy effects induced from physics at very high energy scales M_X near M_P. In the present instance we assume that both the Newton constant, G, and the cosmological term, \Lambda, can be functions of a scale parameter \mu. It turns out that G(\mu) evolves according to a logarithmic law which may lead to asymptotic freedom of gravity, similar to the gauge coupling in QCD. At the same time \Lambda(\mu) evolves quadratically with \mu. We study the consistency and cosmological consequences of these laws when \mu=H. Furthermore, we propose to extend this method to the astrophysical domain after identifying the local RG scale at the galactic level. It turns out that Kepler's third law of celestial mechanics receives quantum corrections that may help to explain the flat rotation curves of the galaxies without introducing the dark matter hypothesis. The origin of these effects (cosmological and astrophysical) could be linked, in our framework, to physics at M_X= 10^{16-17} GeV.

The high resolution obtained through the use of VLBI gives an unique opportunity to directly observe the interaction of an expanding radio supernova with its surrounding medium. We present here results from our VLBI observations of the young supernovae SN 1979C, SN 1986J, and SN 2001gd.

Within the Quantum Field Theory context the idea of a "cosmological constant" (CC) evolving with time looks quite natural as it just reflects the change of the vacuum energy with the typical energy of the universe. In the particular frame of Ref.[30], a "running CC" at low energies may arise from generic quantum effects near the Planck scale, M_P, provided there is a smooth decoupling of all massive particles below M_P. In this work we further develop the cosmological consequences of a "running CC" by addressing the accelerated evolution of the universe within that model. The rate of change of the CC stays slow, without fine-tuning, and is comparable to H^2 M_P^2. It can be described by a single parameter, \nu, that can be determined from already planned experiments using SNe Ia at high z. The range of allowed values for \nu follow mainly from nucleosynthesis restrictions. Present samples of SNe Ia can not yet distinguish between a "constant" CC or a "running" one. The numerical simulations presented in this work show that SNAP can probe the predicted variation of the CC either ruling out this idea or confirming the evolution hereafter expected.

The idea of a variable dark energy has been entertained many times in the literature and from many different points of view. Quintessence is just a popular way to implement this idea in recent times, but so far with little success. Another possibility is to think of the cosmological term, Lambda, as a ``running quantity'' much in the same way as the electromagnetic coupling constant. However, the fact that Lambda is a dimension-four parameter implies that it may obey a peculiar renormalization group equation, which at low energies could be dominated by ``soft decoupling'' contributions of the form Lambda ~ H^2 M_P^2 stemming from physics near the Planck scale. This value lies in the ballpark of the measurements from CMB and high-z supernovae. A ``renormalized'' FLRW cosmology of this kind may reveal itself as a sound, and testable, proposal for a variable Lambda model within quantum field theory in curved space time.

We construct a semiclassical FLRW cosmological model assuming a running cosmological constant (CC). It turns out that the CC becomes variable at arbitrarily low energies due to the remnant quantum effects of the heaviest particles, e.g. the Planck scale physics. These effects are universal in the sense that they lead to a low-energy structure common to a large class of high-energy theories. Remarkably, the uncertainty concerning the unknown high-energy dynamics is accumulated into a single parameter \nu, such that the model has an essential predictive power. Future Type Ia supernovae experiments (like SNAP) can verify whether this framework is correct. For the flat FLRW case and a moderate value \nu ~0.01, we predict an increase of 10-20% in the value of Omega_{Lambda} at redshifts z=1-1.5 perfectly reachable by SNAP.

Starting from generic quantum effects at the Planck scale M_P, we find that the renormalization group running of the cosmological constant (CC) at low energies is possible if there is a smooth decoupling of all massive particles from M_P to the mass of the lightest neutrino. We discuss the theoretical implications of this running for the ``old'' and ``new'' cosmological constant problems. Interestingly enough, the CC running implies a strong relationship between quantum field theory and cosmology, which should be observable in the near future in experiments such as SNAP through the measurement of a cubic redshift dependence of the CC.

We present VLBI images of the RS CVn binary star HR 5110 (=BH CVn; HD 118216), obtained from observations made at 8.4 GHz on 1994 May 29/30 in support of the NASA/Stanford Gravity Probe B project. Our images show an emission region with a core-halo morphology. The core was 0.39 +/- 0.09 mas (FWHM) in size, or 66% +/- 20% of the 0.6 +/- 0.1 mas diameter of the chromospherically active K subgiant star in the binary system. The halo was 1.95 +/- 0.22 mas (FWHM) in size, or 1.8 +/- 0.2 times the 1.1 +/- 0.1 mas separation of the centers of the K and F stars. The core increased significantly in brightness over the course of the observations and seems to have been the site of flare activity that generated an increase in the total flux density of ~200% in 12 hours. The fractional circular polarization simultaneously decreased from ~10% to 2.5%.

We present a model of inflation based on the anomaly-induced effective action of gravity in the presence of a conformally invariant Hilbert-Einstein term. Our approach is based on the conformal representation of the fields action and on the integration of the corresponding conformal anomaly. In contradistinction to the original Starobinsky's model, inflation can be stable at the beginning and unstable at the end. The instability is caused by a slowing down of inflation due to quantum effects associated to the massive fermion fields. In supersymmetric theories this mechanism can be linked to the breaking of SUSY and suggests a natural way to achieve graceful exit from the inflationary to the FLRW phase.

SN1993J is to date the radio supernova whose evolution has been monitored in greatest detail and the one which holds best promise for a comprehensive theoretical-observational analysis. The shell-like radio structure of SN1993J has expanded in general accord with models of shock excited emission, showing almost circular symmetry for over 8 years, except for a bright feature at the south-eastern region of the shell that has been observed at every epoch. The spectrum of SN1993J has flattened from alpha =-1 to alpha =-0.67 (S_(\nu) propto nu**(alpha)). The decelerated expansion can be modeled well with a single slope but apparently better with two slopes. There are also intriguing hints of structure in the expansion curve. The results by the two VLBI groups carrying out this research show general agreement, but also some differences. A comparison of the optical and VLBI results about the details of the deceleration show some discrepancies.

In quantum field theory the parameters of the vacuum action are subject to renormalization group running. In particular, the ``cosmological constant'' is not a constant in a quantum field theory context, still less should be zero. In this paper we continue with previous work, and derive the particle contributions to the running of the cosmological and gravitational constants in the framework of the Standard Model in curved space-time. At higher energies the calculation is performed in a sharp cut off approximation. We assess, in two different frameworks, whether the scaling dependences of the cosmological and gravitational constants spoil primordial nucleosynthesis. Finally, the cosmological implications of the running of the cosmological constant are discussed.

A method of calculating the vacuum effective action for massive quantum fields in curved space-time is outlined. Our approach is based on the conformal representation of the fields action and on the integration of the corresponding conformal anomaly. As a relevant cosmological application, we find that if taking the masses of the fields into account, then the anomaly-induced inflation automatically slows down. The only relevant massive fields for this purpose turn out to be the fermion fields. So in supersymmetric theories this mechanism can be specially efficient, for it may naturally provide the graceful exit from the inflationary to the FLRW phase. Taking the SUSY breaking into account, the anomaly-induced inflation could be free of the well-known difficulties with the initial data and also with the amplitude of the gravitational waves.

The stable version of the anomaly-induced inflation does not need a fine tuning and leads to sufficient expansion of the Universe. The non-stable version (Starobinsky model) provides the graceful exit to the FRW phase. We indicate the possibility of the inflation which is stable at the beginning and unstable at the end. The effect is due to the soft supersymmetry breaking and the decoupling of the massive sparticles at low energy.

We report the detection of the lens cluster of the gravitational lens (GL) system Q0957+561 from a deep observation with the Advanced CCD Imaging Spectrometer on-board the Chandra X-ray Observatory. Intracluster X-ray emission is found to be centered 4.3 +/- 1.3 arcsec east and 3.5(-0.6,+1.3) arcsec north of image B, nearer than previous estimates. Its spectrum can be modeled well with a thermal plasma model consistent with the emission originating from a cluster at a redshift of 0.36. Our best-fit estimates of the cluster temperature of T_e = 2.09(-0.54,+0.83) keV (90 percent confidence) and mass distribution of the cluster are used to derive the convergence parameter kappa, the ratio of the cluster surface mass density to the critical density required for lensing. We estimate the convergence parameter at the location of the lensed images A and B to be kappa_A = 0.22(+0.14,-0.07) and kappa_B = 0.21(+0.12,-0.07), respectively (90 percent confidence levels). The observed cluster center, mass distribution and convergence parameter kappa provide additional constraints to lens models of this system. Our new results break a mass-sheet degeneracy in GL models of this system and provide better constraints of ~ 29 percent (90 percent confidence levels) on the Hubble constant. We also present results from the detection of the most distant X-ray jet (z = 1.41) detected to date. The jet extends approximately 8 arcsec NE of image A and three knots are resolved along the X-ray jet with flux densities decreasing with distance from the core. The observed radio and optical flux densities of the knots are fitted well with a synchrotron model and the X-ray emission is modeled well with inverse Compton scattering of Cosmic Microwave Background photons by synchrotron-emitting electrons in the jet.

On 18 November 1993, we observed the QSO 1150+812 and the BL Lac object 1803+784, nearly 15\deg apart on the sky, with a VLBI array, simultaneously recording at 8.4 and 2.3 GHz. Using difference astrometry techniques, we determined the coordinates of 1803+784 relative to those of 1150+812 to be \Delta\alpha = 6h 7m 33.18469 arcsec +- 0.00020 arcsec \Delta\delta = -2deg 30' 25''.13557 +- 0''.00075 These standard errors contain estimated contributions from the propagation medium, the effects of source opacity, and the possible misidentification of the reference point in some of the maps. We obtained comparably accurate and consistent relative positions of the radio sources with GPS-based ionosphere corrections, thus demonstrating that dual-frequency observations are not required for state-of-the-art accuracy in VLBI astrometry. These results also demonstrate the feasibility of using single-frequency measurements for sources separated by 15$^\circ$ on the sky, and open the avenue for the application of this technique on a full-sky scale.

Global Positioning System (GPS) satellites broadcast at frequencies of 1,575.42 MHz (L1) and of 1,227.60 MHz (L2). The dispersive property of the ionosphere can be used to combine independent measurements at the two frequencies to estimate the total electron content (TEC) between a GPS receiver site and a broadcasting satellite. Such measurements, made at sites near to Very Long Baseline Interferometry (VLBI) sites, can be used to estimate the ionospheric contribution to VLBI observables. For our 1991.9 astrometric VLBI experiment in which we obtained group-delay observations in the 8.4 and 2.3 GHz bands simultaneously, we found that the GPS and VLBI determinations of the ionosphere delays agreed with root-mean-square differences below 0.15 ns for intercontinental baselines and 0.10 ns for continental ones. We also successfully applied the GPS-based procedure to reduce the ionospheric effect in phase delays used for high precision differenced astrometry at 8.4 GHz for this same experiment.

A large value of the cosmological constant (CC) is induced in the Standard Model (SM) of Elementary Particle Physics because of Spontaneous Symmetry Breaking. To provide a small value of the observable CC one has to introduce the vacuum term which cancels the induced one at some point in the very far infrared cosmic scale. Starting from this point we investigate whether the cancellation is preserved at different energy scales. We find that the running of the Higgs mass, couplings and the vacuum term inevitably result in a scaling dependence of the observable CC value. As a consequence one meets a nonzero CC at an energy scale comparable to the typical electron neutrino mass suggested by some experiments, and the order of magnitude of this constant is roughly the one derived from recent supernovae observations. However the sign of it is negative -- opposite to what is suggested by these observations. This discrepancy may be a hint of the existence of an extra very light scalar, perhaps a Cosmon-like dilaton, which should essentially decouple from the SM Lagrangian, but that it nevertheless could mediate new macroscopic forces in the submillimeter range.

We determined the separations of the radio sources in the triangle formed by the BL Lac objects 1803+784 and 2007+777, and the QSO 1928+738 from intercontinental interferometric observations carried out in November 1991 at the frequencies of 2.3 and 8.4 GHz simultaneously. We successfully connected differenced phase delays over 7^o on the sky at 8.4 GHz at an epoch of maximum solar activity. The effects of the ionosphere on these VLBI data were mostly removed by estimates of the total electron content from observations of GPS satellites. The comparison of the estimated separation of QSO 1928+738 and BL 2007+777 with previous such estimates obtained from data at different epochs leads us to a particular alignment of the maps of QSO 1928+738 at different epochs relative to those of BL 2007+777, although with significant uncertainty. For this alignment, the jet components of QSO 1928+738 show a mean proper motion of 0.32+/-0.10 mas/yr and also suggest an identification for the position of the core of this radio source.

We examine models of the mass distribution for the first known case of gravitational lensing. Several new sets of constraints are used, based on recent observations. We remodel the VLBI observations of the radio jets in the two images of Q0957+561, showing that the previously derived positions and uncertainties were incorrect. We use as additional constraints the candidate lensed pairs of galaxies identified recently with the Hubble Space Telescope. We explore a wider range of lens models than before, and find that the Hubble constant is not tightly constrained once elliptical lens models are considered. We also discuss the systematic uncertainties caused by the cluster containing the lens galaxy. We conclude that additional observations of the candidate lensed galaxies as well as direct measurements of the cluster mass profile are needed before an accurate value of the Hubble constant can be derived from this lens.

Analysis of several recent ROSAT HRI observations of the gravitationally lensed system Q0957+561 has led to the detection at the 3sigma level of the cluster lens containing the primary galaxy G1. The total mass was estimated by applying the equation of hydrostatic equilibrium to the detected hot intracluster gas for a range of cluster core radii, cluster sizes and for different values of the Hubble constant. X-ray estimates of the lensing cluster mass provide a means to determine the cluster contribution to the deflection of rays originating from the quasar Q0957+561. The present mass estimates were used to evaluate the convergence parameter kappa, the ratio of the local surface mass density of the cluster to the critical surface mass density for lensing. The convergence parameter, kappa, calculated in the vicinity of the lensed images, was found to range between 0.07 and 0.21, depending on the assumed cluster core radius and cluster extent. This range of uncertainty in kappa does not include possible systematic errors arising from the estimation of the cluster temperature through the use of the cluster luminosity-temperature relation and the assumption of spherical symmetry of the cluster gas. Applying this range of values of kappa to the lensing model of Grogin & Narayan (1996) for Q0957+561 but not accounting for uncertainties in that model yields a range of values for the Hubble constant:67<H_0<82 km s^-1 Mpc^-1, for a time delay of 1.1 years.

A rarity among supernova, SN 1993J in M81 can be studied with high spatial resolution. Its radio power and distance permit VLBI observations to monitor the expansion of its angular structure. This radio structure was previously revealed to be shell-like and to be undergoing a self-similar expansion at a constant rate. From VLBI observations at the wavelengths of 3.6 and 6 cm in the period 6 to 42 months after explosion, we have discovered that the expansion is decelerating. Our measurement of this deceleration yields estimates of the density profiles of the supernova ejecta and circumstellar material in standard supernova explosion models.

We present long-slit LRIS/Keck spectroscopic observations of the gravitational lens system 0957+561. Averaged over all of our data, the rest-frame velocity dispersion sigma_v of the central lens galaxy G1 is sigma_v = 279 +- 12 km/s. However, there appears to be a significant decrease in sigma_v as a function of distance from the center of G1 that is not typical of brightest cluster galaxies. Within 0.2" of the center of G1, we find the average sigma_v = 316+- 14 km/s, whereas for positions > 0.2" from the center of G1, we find the average sigma_v= 266+-12 km/s. A plausible explanation is that G1 contains a central massive dark object of mass M_MDO ~ 4x10^9 h_100^{-1} M_sun (h_100 = H_0/100km/s/Mpc), which contributes to the central velocity dispersion, and that the outer value of sigma_v is the appropriate measure of the depth of the potential well of G1. The determination of a luminosity-weighted estimate of sigma_v is essential for a determination of H_0 from 0957+561; our accurate measurements remove one of the chief uncertainties in the assumed form of the mass distribution of the lens. Thus, with the recent apparent reduction in the uncertainty in the measurement of the time delay for the images A and B of 0957+561, Delta tau_BA = 417+-3 d (Kundic et al. 1996), we obtain an estimate for the Hubble constant: H_0 = 62+-7 km/s/Mpc. If for some reason the trend of sigma_v with slit position is spurious and we should use the dispersion averaged along the slit, then the estimate of H_0 increases to 67+-8 km/s/Mpc. These standard errors do not, however, include any contribution from any errors in the assumed form of the mass distribution of the lens. In particular, we used the mass model described by Falco, Gorenstein & Shapiro (1991), as updated by Grogin & Narayan (1996a, b). The reduced chi^2

Quadruple gravitational lens systems offer the possibility of measuring time delays for image pairs, microlensing effects, and extinction in distant galaxies. Observations of these systems may be used to obtain estimates of H_o and to study the various mass components of lens galaxies at high redshifts. With HST, we have observed the reddest known gravitational lens system, MG 0414+0534. We used WFPC2/PC1 to obtain deep, high-resolution images with two filters, F675W and F814W. We present a detailed analysis of all of the components, as well as macrolens models. Our main results are: (1) confirmation that MG 0414+0534 is inescapably a gravitational lens system; (2) discovery of a blue arc connecting the 3 brightest images of the QSO central core; (3) accurate positions and apparent brightnesses for all 4 known images of the QSO central core and for the lens galaxy G; (4) a good representation of the brightness distribution of G by elliptical isophotes with a De Vaucouleurs profile, characteristic of an elliptical galaxy; (5) models that consist of simple elliptical potentials and account qualitatively, not quantitatively, for the HST image positions, arc morphology and radio flux ratios for the images of the QSO central core; (6) a possible new test to distinguish between reddening in the host galaxy of the QSO and in the lens galaxy, based on future accurate measurements of spatial variations in the color of the arc; and (7) the suggestion that microlensing is a plausible cause for the differences between the radio and optical flux ratios for the brightest images, A1 and A2. Further observations and measurements such as of the redshift of the lens galaxy, can be used fruitfully to study microlensing for this system.

We present hybrid maps of the A and B images of 0957+561 from each of four sessions of 6 cm VLBI observations that span the six-year interval 1987-1993. The inner- and outer-jets are clearly detected, and confirm the structures reported previously. There is no evidence of change in the separation between the core and inner-jet components, so the prospect of measuring the time delay using differential proper motions is not promising. The flux density in the core of each image peaked between 1989 and 1992. From the variation in these flux densities, we obtain a time-delay estimate of $\sim$1 yr.