Extreme Stars Research Articles

Neutron star mergers as sites of r-process nucleosynthesis and short gamma-ray bursts

Neutron star mergers have been long considered as promising sites of heavy [Formula: see text]-process nucleosynthesis. We overview the observational evidence supporting this scenario including: the total amount of [Formula: see text]-process elements in the galaxy, extreme metal-poor stars, geological radioactive elemental abundances, dwarf galaxies and short gamma-ray bursts (sGRBs). Recently, the advanced LIGO and Virgo observatories discovered a gravitational-wave signal of a neutron star merger, GW170817, as well as accompanying multi-wavelength electromagnetic (EM) counterparts. The ultra-violet, optical and near infrared (n/R) observations point to [Formula: see text]-process elements that have been synthesized in the merger ejecta. The rate and ejected mass inferred from GW170817 and the EM counterparts are consistent with other observations. We however, find that, within the simple one zone chemical evolution models (based on merger rates with reasonable delay time distributions as expected from evolutionary models, or from observations of sGRBs), it is difficult to reconcile the current observations of the Eu abundance history of galactic stars for [Fe/H] [Formula: see text]. This implies that to account for the role of mergers in the galactic chemical evolution, we need a galactic model with multiple populations that have different spatial distributions and/or varying formation rates.

Constraints on the Cosmic-Ray Ionization Rate in the z ∼ 2.3 Lensed Galaxies SMM J2135–0102 and SDP 17b from Observations of OH+ and H2O+

Abstract Cosmic rays are predominantly accelerated in shocks associated with star formation such as supernova remnants and stellar wind bubbles, so the cosmic-ray flux and thus cosmic-ray ionization rate, ζ H, should correlate with the star formation rate in a galaxy. Submillimeter bright galaxies (SMGs) are some of the most prolific star-forming galaxies in the universe, and gravitationally lensed SMGs provide bright continuum sources suitable for absorption line studies. Abundances of OH+ and H2O+ are useful for inferring ζ H when combined with chemical models, and have been used for this purpose within the Milky Way. At redshifts z ≳ 2 transitions out of the ground rotational states of OH+ and H2O+ are observable with ALMA, and we present observations of both molecules in absorption toward the lensed SMGs SMM J2135−0102 and SDP 17b. These detections enable an exploration of ζ H in galaxies with extreme star formation and high supernova rates, both of which should significantly enhance cosmic-ray production. The observed OH+ and H2O+ absorption is thought to arise in massive, extended halos of cool, diffuse gas that surround these galaxies. Using a chemical model designed to focus on the reaction network important to both species, we infer cosmic-ray ionization rates of ζ H ∼ 10−16–10−14 s−1 in these extended gaseous halos. Because our estimates come from gas that is far away from the sites of cosmic-ray acceleration, they imply that cosmic-ray ionization rates in the compact regions where star formation occurs in these galaxies are orders of magnitude higher.

Post-common envelope binary stars, radiative levitation, and blue large-amplitude pulsators

Following the discovery of blue large-amplitude pulsators (BLAPs), single star evolu- tion models of post red giant branch stars that have undergone a common envelope (CE) ejection in the form of a high mass loss rate have been constructed and analysed for pulsation stability. The effects of atomic diffusion, particularly radiative levitation, have been examined. Two principal models were considered, being post-CE stars of 0.31 and 0.46 M$_{\odot}$. Such stars are likely, in turn, to become either low-mass helium white dwarfs or core helium-burning extreme horizontal-branch stars. The inclusion of radiative levitation leads to opacity driven pulsations in both types of post-CE object when their effective temperatures are comparable to those of BLAPs, with similar periods. The extent of the instability region for models in these simulations, which are not in thermal balance, is larger than that found for static models, in agreement with previous theory. By comparing to observations, and making some simple evolutionary assumptions, we conclude the 0.31 M$_{\odot}$ star is the more likely candidate for BLAPs. The rate of period change is negative for both cases, so the origin of BLAPs with positive rates of period change remain uncertain.

The End: Witnessing the Death of Extreme Carbon Stars

AbstractA Last affiliation changed 3 to 4 against MS. Please check and confirm if it is fine. small number of the sample of 184 carbon stars in the Magellanic Clouds show signs that they are in the act of evolving off of the asymptotic giant branch. Most carbon stars grow progressively redder in all infrared colors and develop stronger pulsation amplitudes as their circumstellar dust shells become optically thicker. The reddest sources, however, have unexpectedly low pulsation amplitudes, and some even show blue excesses that could point to deviations from spherical symmetry as they eject the last of their envelopes. Previously, all dusty carbon-rich AGB stars have been labeled “extreme,” but that term should be reserved for the truly extreme carbon stars. These objects may well hold the clues needed to disentangle what actually happens when a star ejects the last of its envelope and evolves off of the AGB.

A new method to identify subclasses among AGB stars using Gaia and 2MASS photometry

Aims. We explore the wealth of high-quality photometric data provided by data release 2 (DR2) of the Gaia mission for long-period variables (LPVs) in the Large Magellanic Cloud (LMC). Our goal is to identify stars of various types and masses along the asymptotic giant branch. Methods. For this endeavour, we developed a new multi-band approach combining Wesenheit functions WRP,BP−RP and WKs,J−Ks in the Gaia BP, RP, and 2MASS J, Ks spectral ranges, respectively, and use a new diagram, (WRP,BP−RP − WKs,J−Ks) versus Ks, to distinguish between different kinds of stars in our sample of LPVs. We used stellar population synthesis models to validate our approach. Results. We demonstrate the ability of the new diagram to discriminate between O- and C-rich objects, and to identify low-mass, intermediate-mass, and massive O-rich red giants, as well as extreme C-rich stars. Stellar evolution and population synthesis models guide the interpretation of the results, highlighting the diagnostic power of the new tool to discriminate between stellar initial masses, chemical properties, and evolutionary stages.

Where are the most ancient stars in the Milky Way?

The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disc galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars (z form ≳ 5) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at z = 0 than stars formed at later times as a result of two processes. First, the majority of the oldest stars are not formed in situ but are accreted during hierarchical assembly. These ex situ stars are deposited on dispersion-supported, halo-like orbits but dominate over old stars formed in situ in the solar neighbourhood, and in some simulations, even in the galactic centre. Secondly, old stars formed in situ are driven outwards by bursty star formation and energetic feedback processes that create a time-varying gravitational potential at z ≳ 2, similar to the process that creates dark matter cores and expands stellar orbits in bursty dwarf galaxies. The total fraction of stars that are ancient is more than an order of magnitude higher for sight lines away from the bulge and inner halo than for inward-looking sight lines. Although the task of identifying specific stars as ancient remains challenging, we anticipate that million-star spectral surveys and photometric surveys targeting metal-poor stars already include hundreds of stars formed before z = 5. We predict most of these targets to have higher metallicity (-3 < [Fe/H] < -2) than the most extreme metal-poor stars.

Self-enrichment in globular clusters: the extreme He-rich populationof NGC 2808

Almost several decades after the discovery of the first multiple populations in galactic globular clusters (GC) the debate on their formation is still extremely current and NGC2808 remains one of the best benchmark to test any scenario for their origin and the evolution. In this work we focus on the chemical composition of stars belonging to the extreme He-rich population populated by stars with the most extreme abundance of Mg, Al, Na, O and Si. We checked whether the most recent measures are consistent with the AGB yields of stars of $6.5-8~M_{\odot}$. These stars evolve on time scales of the order of 40-60 Myr and eject matter strongly enriched in helium, owing to a deep penetration of the surface convective zone down to regions touched by CNO nucleosynthesis occurring after the core He-burning phase. Since the big unknown of the AGB phase of massive stars is the mass loss, we propose a new approch that takes into account the effects of the radiation pressure on dust particles. We show that this more realistic description is able to reproduce the observed abundances of Mg, Al, Na and Si in these extreme stars. The large spread in the oxygen abundances is explained by invoking deep mixing during the RGB phase. It will be possible to check this work hypothesis as soon as the oxygen measurements of the main sequence stars of NGC2808 will be available.

Molecular gas in AzTEC/C159: a star-forming disk galaxy 1.3 Gyr after the Big Bang

We studied the molecular gas properties of AzTEC/C159, a star-forming disk galaxy at z = 4.567, in order to better constrain the nature of the high-redshift end of the submillimeter-selected galaxy (SMG) population. We secured 12CO molecular line detections for the J = 2 →1 and J = 5 →4 transitions using the Karl G. Jansky Very Large Array (VLA) and the NOrthern Extended Millimeter Array (NOEMA) interferometer. The broad (FWHM ~ 750 km s−1) and tentative double-peaked profiles of the two 12CO lines are consistent with an extended molecular gas reservoir, which is distributed in a rotating disk, as previously revealed from [CII] 158 μm line observations. Based on the 12CO(2 →1) emission line, we derived L′CO=(3.4±0.6)×1010 K km s−1 pc2, which yields a molecular gas mass of MH2(αCO/4.3)=(1.5±0.3)×1011 M⊙ and unveils a gas-rich system with μgas(αCO/4.3)≡MH2/M⋆=3.3±0.7. The extreme star formation efficiency of AzTEC/C159, parametrized by the ratio LIR/L′CO=(216±80) L⊙ (K km s−1 pc2)−1, is comparable to merger-driven starbursts such as local ultra-luminous infrared galaxies and SMGs. Likewise, the 12CO(5 →4)/CO(2 →1) line brightness temperature ratio of r52 = 0.55 ± 0.15 is consistent with high-excitation conditions as observed in SMGs. Based on mass budget considerations, we constrained the value for the L′CO – H2 mass conversion factor in AzTEC/C159, that is, αCO=3.9−1.3+2.7 M⊙ K−1 km−1 s pc−2, which is consistent with a self-gravitating molecular gas distribution as observed in local star-forming disk galaxies. Cold gas streams from cosmological filaments might be fueling a gravitationally unstable gas-rich disk in AzTEC/C159, which breaks into giant clumps and forms stars as efficiently as in merger-driven systems and generates high gas excitation. These results support the evolutionary connection between AzTEC/C159-like systems and massive quiescent disk galaxies at z ~ 2.

A universal route for the formation of massive star clusters in giant molecular clouds

Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show that GC-sized clusters are naturally produced in radiation-hydrodynamic simulations of isolated 10$^7$ M$_{\odot}$ Giant Molecular Clouds (GMCs) with properties typical of the local universe, even under the influence of radiative feedback. In all cases, these massive clusters grow to GC-level masses within 5 Myr via a roughly equal combination of filamentary gas accretion and mergers with several less massive clusters. Lowering the heavy-element abundance of the GMC by a factor of 10 reduces the opacity of the gas to radiation and better represents the high-redshift formation conditions of GCs. This results in higher gas accretion leading to a mass increase of the largest cluster by a factor of ~4. When combined with simulations of less massive GMCs (10$^{4-6}$ M$_{\odot}$), a clear relation emerges between the maximum YMC mass and the mass of the host GMC. Our results demonstrate that YMCs, and potentially GCs, are a simple extension of local cluster formation to more massive clouds and do not require suggested exotic formation scenarios.

The Growth of Stellar Mass Black Hole Binaries Trapped in the Accretion Disks of Active Galactic Nuclei

Abstract Among the four black hole (BH) binary merger events detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO), six progenitor BHs have masses greater than 20 M ⊙. The existence of such massive BHs suggests that extreme metal-poor stars are the progenitors. An alternative possibility, that a pair of stellar mass BHs each with mass ∼7 M ⊙ increases to &gt;20 M ⊙ via accretion from a disk surrounding a supermassive BH (SMBH) in an active galactic nucleus (AGN), is considered. The growth of mass of the binary and the transfer of orbital angular momentum to the disk accelerates the merger. Based on the recent numerical work of Tang et al., it is found that, in the disk of a low-mass AGN with mass ∼106 M ⊙ and Eddington ratio &gt;0.01, the mass of an individual BH in the binary can grow to &gt;20 M ⊙ before coalescence, provided that accretion takes place at a rate more than 10 times the Eddington value. This mechanism predicts a new class of gravitational wave (GW) sources involving the merger of two extreme Kerr black holes associated with AGNs and a possible electromagnetic wave counterpart.

Evidence for a maximum mass cut-off in the neutron star mass distribution and constraints on the equation of state

We infer the mass distribution of neutron stars in binary systems using a flexible Gaussian mixture model and use Bayesian model selection to explore evidence for multi-modality and a sharp cut-off in the mass distribution. We find overwhelming evidence for a bimodal distribution, in agreement with previous literature, and report for the first time positive evidence for a sharp cut-off at a maximum neutron star mass. We measure the maximum mass to be $2.0M_\odot < m_\mathrm{max} < 2.2M_\odot$ (68\%), $2.0M_\odot < m_\mathrm{max}< 2.6M_\odot$ (90\%), and evidence for a cut-off is robust against the choice of model for the mass distribution and to removing the most extreme (highest mass) neutron stars from the dataset. If this sharp cut-off is interpreted as the maximum stable neutron star mass allowed by the equation of state of dense matter, our measurement puts constraints on the equation of state. For a set of realistic equations of state that support $>2M_\odot$ neutron stars, our inference of $m_\mathrm{max}$ is able to distinguish between models at odds ratios of up to $12:1$, whilst under a flexible piecewise polytropic equation of state model our maximum mass measurement improves constraints on the pressure at $3-7\times$ the nuclear saturation density by $\sim 30-50\%$ compared to simply requiring $m_\mathrm{max}> 2M_\odot$. We obtain a lower bound on the maximum sound speed attained inside the neutron star of $c_s^\mathrm{max} > 0.63c$ (99.8\%), ruling out $c_s^\mathrm{max} < c/\sqrt{3}$ at high significance. Our constraints on the maximum neutron star mass strengthen the case for neutron star-neutron star mergers as the primary source of short gamma-ray bursts.

HE 0430–2457: a post-merger extremely low-mass pre-white dwarf in a wide binary posing as an extreme horizontal branch star

Abstract We report the discovery of HE 0430–2457, the first extremely low-mass pre-white dwarf (ELM pre-WD) in a long period binary (P = 771 ± 3 d). The spectroscopic parameters of the primary are determined to be Teff = 26 200 ± 1500 K and log g = 5.40 ± 0.35, placing it in the region occupied by core He-burning hot subdwarf B stars. By comparing the spectroscopic parameters of the K-type companion to stellar models, and using the mass ratio, the mass of the hot primary is determined to be 0.23 M⊙. Given that this is too low for core He-burning, the primary in HE 0430–2457 is not an extreme horizontal branch (EHB) star but a pre-WD of the ELM type. As the lifetime of ELM pre-WDs in this region of the Hertzsprung Russel diagram populated by EHBs is thought to be very short, they are not considered to be part of the observed EHBs. However, the discovery of this system indicates that the percentage of ELM pre-WDs in the observed EHB population might be higher than previously thought. Binary evolution models indicate that HE 0430–2457 is likely formed by a merger of the inner binary in a hierarchical triple system.

Extreme star formation in the Milky Way: luminosity distributions of young stellar objects in W49A and W51

We have compared the star-formation properties of the W49A and W51 regions by using far-infrared data from the Herschel infrared Galactic Plane Survey (Hi-GAL) and 850-um observations from the James Clerk Maxwell Telescope (JCMT) to obtain luminosities and masses, respectively, of associated compact sources. The former are infrared luminosities from the catalogue of Elia et al. (2017), while the latter are from the JCMT Plane survey source catalogue as well as measurements from new data. The clump-mass distributions of the two regions are found to be consistent with each other, as are the clump-formation efficiency and star-formation efficiency analogues. However, the frequency distributions of the luminosities of the young stellar objects are significantly different. While the luminosity distribution in W51 is consistent with Galaxy-wide samples, that of W49A is top-heavy. The differences are not dramatic, and are concentrated in the central regions of W49A. However, they suggest that physical conditions there, which are comparable in part to those in extragalactic starbursts, are significantly affecting the star-formation properties or evolution of the dense clumps in the region.

Noether's stars in [formula omitted] gravity

The Noether Symmetry Approach can be used to construct spherically symmetric solutions in f(R) gravity. Specifically, the Noether conserved quantity is related to the gravitational mass and a gravitational radius that reduces to the Schwarzschild radius in the limit f(R)→R. We show that it is possible to construct the M–R relation for neutron stars depending on the Noether conserved quantity and the associated gravitational radius. This approach enables the recovery of extreme massive stars that could not be stable in the standard Tolman–Oppenheimer–Volkoff based on General Relativity. Examples are given for some power law f(R) gravity models.

The effects of diffusion in hot subdwarf progenitors from the common envelope channel

Diffusion of elements in the atmosphere and envelope of a star can drastically alter its surface composition, leading to extreme chemical peculiarities. We consider the case of hot subdwarfs, where surface helium abundances range from practically zero to almost 100 percent. Since hot subdwarfs can form via a number of different evolution channels, a key question concerns how the formation mechanism is connected to the present surface chemistry. A sequence of extreme horizontal branch star models was generated by producing post-common envelope stars from red giants. Evolution was computed with MESA from envelope ejection up to core-helium ignition. Surface abundances were calculated at the zero-age horizontal branch for models with and without diffusion. A number of simulations also included radiative levitation. The goal was to study surface chemistry during evolution from cool giant to hot subdwarf and determine when the characteristic subdwarf surface is established. Only stars leaving the giant branch close to core-helium ignition become hydrogen-rich subdwarfs at the zero-age horizontal branch. Diffusion, including radiative levitation, depletes the initial surface helium in all cases. All subdwarf models rapidly become more depleted than observations allow. Surface abundances of other elements follow observed trends in general, but not in detail. Additional physics is required.

Diffusion in hot subdwarf progenitors from the common envelope channel

Abstract Diffusion of elements in the atmosphere of a star can drastically affect its surface composition, sometimes leading to unusual mixtures. These chemically peculiar stars can be identified fromthe presence of unusual lines in their spectra. Some hot subdwarf stars show extraordinary abundances of elements such as lead, zirconium and strontium, while the abundance of helium ranges from practically zero to almost 100 percent across the hot subdwarf population. A sequence of extreme horizontal branch star models was generated by producing a number of post-common envelope objects from red giants. The evolution of these subdwarf progenitors was computed with the MESA stellar evolution code from immediately after envelope ejection right up to the ignition of helium in the core. Envelope abundances were calculated at the zero age horizontal branch for models both with and without the presence of diffusion. A small number of simulations also looked at the effects on radiative levitation of these abundances, to test how well diffusion physics is able to reproduce observational data.

The star-forming complex LMC-N79 as a future rival to 30 Doradus

Within the early Universe, ‘extreme’ star formation may have been the norm rather than the exception 1,2 . Super star clusters (with masses greater than 105 solar masses) are thought to be the modern-day analogues of globular clusters, relics of a cosmic time (redshift z ≳ 2) when the Universe was filled with vigorously star-forming systems 3 . The giant H ii region 30 Doradus in the Large Magellanic Cloud is often regarded as a benchmark for studies of extreme star formation 4 . Here, we report the discovery of a massive embedded star-forming complex spanning about 500 pc in the unexplored southwest region of the Large Magellanic Cloud, which manifests itself as a younger, embedded twin of 30 Doradus. Previously known as N79, this region has a star-formation efficiency greater than that of 30 Doradus, by a factor of about 2, as measured over the past ~0.5 Myr. Moreover, at the heart of N79 lies the most luminous infrared compact source discovered with large-scale infrared surveys of the Large Magellanic Cloud and Milky Way, possibly a precursor to the central super star cluster of 30 Doradus, R136. The discovery of a nearby candidate super star cluster may provide invaluable information to understand how extreme star formation proceeds in the current and high-redshift Universe. The relatively unexplored southwestern region of the Large Magellanic Cloud is host to a massive, embedded star-forming complex that rivals the star-forming efficiency of 30 Doradus. Its most luminous object could be a super star cluster in formation.

Spectroscopic twin to the hypervelocity sdO star US 708 and three fast sdB stars from the Hyper-MUCHFUSS project

Important tracers for the dark matter halo of the Galaxy are hypervelocity stars (HVSs), which are faster than the local escape velocity of the Galaxy and their slower counterparts, the high-velocity stars in the Galactic halo. Such HVSs are believed to be ejected from the Galactic Centre (GC) through tidal disruption of a binary by the super-massive black hole (Hills mechanism). The Hyper-MUCHFUSS survey aims at finding high-velocity potentially unbound hot subdwarf stars. We present the spectroscopic and kinematical analyses of a He-sdO as well as three candidates among the sdB stars using optical Keck/ESI and VLT (Xshooter, FORS) spectroscopy. Proper motions are determined by combining positions from early-epoch photographic plates with those derived from modern digital sky surveys. The Galactic rest frame velocities range from 203 km\,s$^{-1}$ to 660 km\,s$^{-1}$, indicating that most likely all four stars are gravitationally bound to the Galaxy. With $T_\text{eff}=47000$ K and a surface gravity of $\log g = 5.7$, SDSS J205030.39$-$061957.8 (J2050) is a spectroscopic twin of the hypervelocity He-sdO US 708. As for the latter, the GC is excluded as a place of origin based on the kinematic analysis. Hence, the Hills mechanism can be excluded for J2050. The ejection velocity is much more moderate ($385\pm79$ km\,s$^{-1}$) than that of US 708 ($998\pm68$ km\,s$^{-1}$). The binary thermonuclear supernova scenario suggested for US 708 would explain the observed properties of J2050 very well without pushing the model parameters to their extreme limits, as required for US 708. Accordingly, the star would be the surviving donor of a type Ia supernova. Three sdB stars also showed extreme kinematics; one could be a HVS ejected from the GC, whereas the other two could be ejected from the Galactic disk through the binary supernova mechanism. Alternatively, they might be extreme halo stars.

ZFIRE: using Hα equivalent widths to investigate the in situ initial mass function at z ∼ 2

We use the ZFIRE (http://zfire.swinburne.edu.au) survey to investigate the high-mass slope of the initial mass function (IMF) for a mass-complete (⁠|$\mathrm{log_{10}({\it M}_*/\mathrm{M}_{\odot })\sim 9.3}$|⁠) sample of 102 star-forming galaxies at z ∼ 2 using their Hα equivalent widths (Hα EWs) and rest-frame optical colours. We compare dust-corrected Hα EW distributions with predictions of star formation histories (SFHs) from pegase.2 and starburst|$\scriptstyle{99}$| synthetic stellar population models. We find an excess of high Hα EW galaxies that are up to 0.3–0.5 dex above the model-predicted Salpeter IMF locus and the Hα EW distribution is much broader (10–500 Å) than can easily be explained by a simple monotonic SFH with a standard Salpeter-slope IMF. Though this discrepancy is somewhat alleviated when it is assumed that there is no relative attenuation difference between stars and nebular lines, the result is robust against observational biases, and no single IMF (i.e. non-Salpeter slope) can reproduce the data. We show using both spectral stacking and Monte Carlo simulations that starbursts cannot explain the EW distribution. We investigate other physical mechanisms including models with variations in stellar rotation, binary star evolution, metallicity and the IMF upper-mass cut-off. IMF variations and/or highly rotating extreme metal-poor stars (Z ∼ 0.1 Z⊙) with binary interactions are the most plausible explanations for our data. If the IMF varies, then the highest Hα EWs would require very shallow slopes (Γ > −1.0) with no one slope able to reproduce the data. Thus, the IMF would have to vary stochastically. We conclude that the stellar populations at z ≳ 2 show distinct differences from local populations and there is no simple physical model to explain the large variation in Hα EWs at z ∼ 2.

ALMACAL II: Extreme Star Formation Rate Densities in Dusty Starbursts Revealed by ALMA 20 mas Resolution Imaging

Abstract We present ultrahigh spatial resolution (∼20 mas or 150 pc) ALMA observations of the dust continuum at 920 μm and 1.2 mm in two submillimeter sources at z = 3.442, ALMACAL–1 (A–1: ) and ALMACAL–2 (A–2: ). About half of the star formation in each of these sources is dominated by a single compact clump (FWHM size of ∼350 pc). In A–1, two additional fainter clumps are found. The star formation rate (SFR) surface densities of all these clumps are extremely high, to , the highest rates found in high-redshift galaxies. Given their geometry and identical redshifts, there is a possibility that A–1 and A–2 are the lensed images of a single background source that are gravitationally amplified by the blazar host. If this were the case, the effective radius of the dusty galaxy in the source plane would be and the demagnified SFR surface density would be ∼ 10,000 , comparable with the eastern nucleus of Arp 220. Although we cannot rule out an AGN contribution, our results suggest that a significant percentage of the enormous far-IR luminosity in some dusty starbursts is extremely compact. The high in these sources could only be measured thanks to the ultrahigh-resolution ALMA observations used in this work, demonstrating that long-baseline observations are essential to study and interpret the properties of dusty starbursts in the early Universe.