Context. We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Aims. A summary of the contents of Gaia EDR3 is presented, accompanied by a discussion on the differences with respect to Gaia DR2 and an overview of the main limitations which are present in the survey. Recommendations are made on the responsible use of Gaia EDR3 results. Methods. The raw data collected with the Gaia instruments during the first 34 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium and turned into this early third data release, which represents a major advance with respect to Gaia DR2 in terms of astrometric and photometric precision, accuracy, and homogeneity. Results. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the (GBP − GRP) colour are also available. The passbands for G, GBP, and GRP are provided as part of the release. For ease of use, the 7 million radial velocities from Gaia DR2 are included in this release, after the removal of a small number of spurious values. New radial velocities will appear as part of Gaia DR3. Finally, Gaia EDR3 represents an updated materialisation of the celestial reference frame (CRF) in the optical, the Gaia-CRF3, which is based solely on extragalactic sources. The creation of the source list for Gaia EDR3 includes enhancements that make it more robust with respect to high proper motion stars, and the disturbing effects of spurious and partially resolved sources. The source list is largely the same as that for Gaia DR2, but it does feature new sources and there are some notable changes. The source list will not change for Gaia DR3. Conclusions. Gaia EDR3 represents a significant advance over Gaia DR2, with parallax precisions increased by 30 per cent, proper motion precisions increased by a factor of 2, and the systematic errors in the astrometry suppressed by 30–40% for the parallaxes and by a factor ~2.5 for the proper motions. The photometry also features increased precision, but above all much better homogeneity across colour, magnitude, and celestial position. A single passband for G, GBP, and GRP is valid over the entire magnitude and colour range, with no systematics above the 1% level

We present multiband photometry of 185 type-Ia supernovae (SNe Ia), with over 11,500 observations. These were acquired between 2001 and 2008 at the F. L. Whipple Observatory of the Harvard-Smithsonian Center for Astrophysics (CfA). This sample contains the largest number of homogeneously observed and reduced nearby SNe Ia (z 0.08) published to date. It more than doubles the nearby sample, bringing SN Ia cosmology to the point where systematic uncertainties dominate. Our natural system photometry has a precision of 0.02 mag in BV RI r i and 0.04 mag in U for points brighter than 17.5 mag. We also estimate a systematic uncertainty of 0.03 mag in our SN Ia standard system BV RI r i photometry and 0.07 mag for U. Comparisons of our standard system photometry with published SN Ia light curves and comparison stars, where available for the same SN, reveal agreement at the level of a few hundredths mag in most cases. We find that 1991bg-like SNe Ia are sufficiently distinct from other SNe Ia in their color and light-curve-shape/ luminosity relation that they should be treated separately in light-curve/distance fitter training samples. The CfA3 sample will contribute to the development of better light-curve/distance fitters, particularly in the few dozen cases where near-infrared photometry has been obtained and, together, can help disentangle host-galaxy reddening from intrinsic supernova color, reducing the systematic uncertainty in SN Ia distances due to dust.

Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray observations, but clear counterparts at longer wavelengths have not been seen. We present a flare observed over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, the Caltech Submillimeter Observatory, and the Submillimeter Array. All telescopes observed strong flare events, but the submillimeter peak is found to occur nearly 100 minutes after the X-ray peak. Submillimeter polarization data show linear polarization in the excess flare emission, increasing from 9% to 17% as the flare passes through its peak, consistent with a transition from optically thick to thin synchrotron emission. The temporal and spectral behavior of the flare require that the energetic electrons responsible for the emission cool faster than expected from their radiative output. This is consistent with adiabatic cooling in an expanding emission region, with X-rays produced through self-Compton scattering, although not consistent with the simplest model of such expansion. We also present a submillimeter flare that followed a bright IR flare on 2005 July 31. Compared to 2006, this event had a larger peak IR flux and similar submillimeter flux, but it lacked measurable X-ray emission. It also showed a shorter delay between the IR and submillimeter peaks. Based on these events we propose a synchrotron and self-Compton model to relate the submillimeter lag and the variable IR/X-ray luminosity ratio.

We carry out 2D viscous hydrodynamical simulations of circumbinary accretion using the moving-mesh code AREPO. We self-consistently compute the accretion flow over a wide range of spatial scales, from the circumbinary disk (CBD) far from the central binary, through accretion streamers, to the disks around individual binary components, resolving the flow down to 2% of the binary separation. We focus on equal-mass binaries with arbitrary eccentricities. We evolve the flow over long (viscous) timescales until a quasi-steady state is reached, in which the mass supply rate at large distancesṀ 0 (assumed constant) equals the time-averaged mass transfer rate across the disk and the total mass accretion rate onto the binary components. This quasi-steady state allows us to compute the secular angular momentum transfer rate onto the binary, J b , and the resulting orbital evolution. Through direct computation of the gravitational and accretional torques on the binary, we find that J b is consistently positive (i.e., the binary gains angular momentum), with l 0 ≡ J b /Ṁ 0 in the range of (0.4 − 0.8)a 2 b Ω b , depending on the binary eccentricity (where a b , Ω b are the binary semi-major axis and angular frequency); we also find that this J b is equal to the net angular momentum current across the CBD, indicating that global angular momentum balance is achieved in our simulations. In addition, we compute the time-averaged rate of change of the binary orbital energy for eccentric binaries, and thus obtain the secular rates ȧ b and ė b . In all cases, ȧ b is positive, i.e., the binary expands while accreting. We discuss the implications of our results for the merger of supermassive binary black holes and for the formation of close stellar binaries.

We carry out numerical simulations of circumbinary discs, solving the viscous hydrodynamics equations on a polar grid covering an extended disc outside the binary co-orbital region. We use carefully controlled outer boundary conditions and longterm integrations to ensure that the disc reaches a quasi-steady state, in which the time-averaged mass accretion rate onto the binary, Ṁ , matches the mass supply rate at the outer disc. We focus on binaries with comparable masses and a wide range of eccentricities (e B ). For e B 0.05, the mass accretion rate of the binary is modulated at about 5 times the binary period; otherwise it is modulated at the binary period. The inner part of the circumbinary disc (r 6a B ) generally becomes coherently eccentric. For low and high e B , the disc line of apsides precesses around the binary, but for intermediate e B (0.2 − 0.4), it instead becomes locked with that of the binary. By considering the balance of angular momentum transport through the disc by advection, viscous stress, and gravitational torque, we determine the time-averaged net angular momentum transfer rate to the binary, J . The specific angular momentum, l 0 = J / Ṁ , depends non-monotonically on e B . Contrary to previous claims, we find that l 0 is positive for most e B , implying that the binary receives net angular momentum, which may cause its separation to grow with time. The minimum l 0 occurs at intermediate e B (0.2 − 0.4), corresponding to the regime where the inner eccentric disc is apsidally aligned with the binary.

Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code Arepo and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems Arepo only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the L 1 norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.

The Lidov-Kozai (LK) mechanism plays an important role in the secular evolution of many hierarchical triple systems. The standard LK mechanism consists of large-amplitude oscillations in eccentricity and inclination of a binary subject to the quadrupole potential from an outer perturber. Recent work has shown that when the octupole terms are included in the potential, the inner binary can reach more extreme eccentricities as well as undergo orientation flips. It is known that pericenter precessions due to short-range effects, such as General Relativity and tidal and rotational distortions, can limit the growth of eccentricity and even suppress standard (quadrupolar) LK oscillations, but their effect on the octupole-level LK mechanism has not been fully explored. In this paper, we systematically study how these short-range forces affect the extreme orbital behaviour found in octupole LK cycles. In general, the influence of the octupole potential is confined to a range of initial mutual inclinations i tot centered around 90• (when the inner binary mass ratio is 1), with this range expanding with increasing octupole strength. We find that, while the short-range forces do not change the width and location of this "window of influence", they impose a strict upper limit on the maximum achievable eccentricity. This limiting eccentricity can be calculated analytically, and its value holds even for strong octupole potential and for the general case of three comparable masses. Short-range forces also affect orbital flips, progressively reducing the range of i tot within which flips are possible as the intensity of these forces increases.

We present numerical simulations of circumbinary accretion onto eccentric and circular binaries using the moving-mesh code AREPO. This is the first set of simulations to tackle the problem of binary accretion using a finite-volume scheme on a freely moving mesh, which allows for accurate measurements of accretion onto individual stars for arbitrary binary eccentricity. While accretion onto a circular binary shows bursts with period of ∼5 times the binary period P b , accretion onto an eccentric binary is predominantly modulated at the period ∼1P b . For an equal-mass circular binary, the accretion rates onto individual stars are quite similar to each other, following the same variable pattern in time. By contrast, for eccentric binaries, one of the stars can accrete at a rate 10-20 times larger than its companion. This "symmetry breaking" between the stars, however, alternates over timescales of order 200P b , and can be attributed to a slowly precessing, eccentric circumbinary disk. Over longer timescales, the net accretion rates onto individual stars are the same, reaching a quasi-steady state with the circumbinary disk. These results have important implications for the accretion behavior of binary T-Tauri stars and supermassive binary black holes.