Search

The first very long baseline interferometry (VLBI) detections at 870$\mu$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$\lambda$ corresponding to an angular resolution, or fringe spacing, of 19$\mu$as. The Allan deviation of the visibility phase at 870$\mu$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$\mu$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time., Comment: Corresponding author: S. Doeleman

The active galactic nucleus (AGN) feeding and feedback process in the centers of galaxy clusters and groups is still not well understood. NGC 5044 is the ideal system in which to study AGN feedback. It hosts the largest known reservoir of cold gas in any cool-core galaxy group, and features several past epochs of AGN feedback imprinted as cavities in the X-ray bright intragroup medium (IGrM), as well as parsec scale jets. We present Submillimeter Array (SMA), Karl G. Jansky Very Large Array (VLA), James Clerk Maxwell Telescope (JCMT), and Atacama Large Millimeter/submillimeter Array (ALMA) high frequency observations of NGC 5044 to assess the time variability of the mmwaveband emission from the accretion disk, and quantify the Spectral Energy Distribution (SED) from the radio to sub-millimeter band. The SED is well described by advection dominated accretion flow (ADAF) model and self-absorbed jet emission from an aging plasma with tau = 1kyr. We find a characteristic variability timescale of 150 days, which constrains the ADAF emission region to about 0.1pc, and the magnetic field to 4.7mG in the jets and and 870G in the accretion disk. A longer monitoring/sampling will allow to understand if the underlying process is truly periodic in nature., Comment: 17 pages, 7 figures, accepted for publication in ApJ

We present multiwavelength polarization measurements of the luminous blazar Mrk~501 over a 14-month period. The 2--8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100-ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical-infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at $>3\sigma$ above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions., Comment: Accepted for publication in ApJ

We present multi-frequency (5-345 GHz) and multi-resolution radio observations of 1ES 1927+654, widely considered one of the most unusual and extreme changing-look active galactic nuclei (CL-AGN). The source was first designated a CL-AGN after an optical outburst in late 2017 and has since displayed considerable changes in X-ray emission, including the destruction and rebuilding of the X-ray corona in 2019-2020. Radio observations prior to 2023 show a faint and compact radio source typical of radio-quiet AGN. Starting in February 2023, 1ES 1927+654 began exhibiting a radio flare with a steep exponential rise, reaching a peak 60 times previous flux levels, and has maintained this higher level of radio emission for over a year to date. The 5-23 GHz spectrum is broadly similar to gigahertz-peaked radio sources, which are understood to be young radio jets less than ~1000 years old. Recent high-resolution VLBA observations at 23.5 GHz now show resolved extensions on either side of the core, with a separation of ~0.15 pc, consistent with a new and mildly relativistic bipolar outflow. A steady increase in the soft X-ray band (0.3-2 keV) concurrent with the radio may be consistent with jet-driven shocked gas, though further observations are needed to test alternate scenarios. This source joins a growing number of CL-AGN and tidal disruption events which show late-time radio activity, years after the initial outburst., Comment: Submitted to ApJ Letters 24 June 2024; Accepted 14 October 2024

The Submillimeter Array (SMA) is an array of 8 antennas operating at millimeter and submillimeter wavelengths on Maunakea, Hawaii, operated by the Smithsonian Astrophysical Observatory and Academia Sinica Institute of Astronomy and Astrophysics, Taiwan. Over the past several years, we have been preparing a major upgrade to the SMA that will replace the aging original receiver cryostats and receiver cartridges with all new cryostats and new 230 and 345 GHz receiver designs. This wideband upgrade (wSMA) will also include significantly increased instantaneous bandwidth, improved sensitivity, and greater capabilities for dual frequency observations. In this paper, we will describe the wSMA receiver upgrade and status, as well as the future upgrades that will be enabled by the deployment of the wSMA receivers., Comment: To be published in the proceedings of SPIE Astronomical Telescopes + Instrumentation 2024, paper number 13096-122

We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars., Comment: 17 pages, 10 figures, 4 tables, Accepted for publication in A&A

We present a new SMA survey of 47 Class II sources in the Taurus-Auriga region. Our observations made 12 independent samples of flux densities over the 200-400 GHz frequency range. We tightly constrained the spectral indices of most sources to a narrow range of $2.0\pm0.2$; only a handful of spatially resolved (e.g., diameter $>$250 au) disks present larger spectral indices. The simplest interpretation for this result is that the (sub)millimeter luminosities of all of the observed target sources are dominated by very optically thick (e.g., $\tau\gtrsim$5) dust thermal emission. Some previous works that were based on the optically thin assumption thus might have underestimated optical depths by at least one order of magnitude. Assuming DSHARP dust opacities, this corresponds to underestimates of dust masses by a similar factor. Moreover, some population synthesis models show that to explain the observed, narrowly distributed spectral indices, the disks in our selected sample need to have very similar dust temperatures ($T_{\small{dust}}$). Given a specific assumption of median $T_{\small{dust}}$, the maximum grain sizes ($a_{\small{max}}$) can also be constrained, which is a few times smaller than 0.1 mm for $T_{\small{dust}}\sim$100 K and a few mm for $T_{\small{dust}}\sim$24 K. The results may indicate that dust grain growth outside the water snowline is limited by the bouncing/fragmentation barriers. In the Class II disks, the dust mass budget outside of the water snowline may be largely retained instead of being mostly consumed by planet formation. While Class II disks still possess sufficient dust masses to feed planet formation at a later time, it is unknown whether or not dust coagulation and planet formation can be efficient or natural outside of the water snowline., Comment: accepted by Astrophysical Journal Supplement

Context. The blazar AO 0235+164, located at redshift $z=0.94$, has displayed interesting and repeating flaring activity in the past, the latest episodes occurring in 2008 and 2015. In 2020, the source brightened again, starting a new flaring episode that peaked in 2021. Aims. We study the origin and properties of the 2021 flare in relation to previous studies and the historical behavior of the source, in particular to the 2008 and 2015 flaring episodes. Methods. We analyze the multi-wavelength photo-polarimetric evolution of the source. From Very Long Baseline Array images, we derive the kinematic parameters of new components associated with the 2021 flare. We use this information to constrain a model for the spectral energy distribution of the emission during the flaring period. We propose an analytical geometric model to test whether the observed wobbling of the jet is consistent with precession. Results. We report the appearance of two new components that are ejected in a different direction than previously, confirming the wobbling of the jet. We find that the direction of ejection is consistent with that of a precessing jet.The derived period independently agrees with the values commonly found in the literature. Modeling of the spectral energy distribution further confirm that the differences between flares can be attributed to geometrical effects., Comment: Accepted 15 May 2024

Titan's atmospheric composition and dynamical state have previously been studied over numerous epochs by both ground- and space-based facilities. However, stratospheric measurements remain sparse during Titan's northern summer and fall. The lack of seasonal symmetry in observations of Titan's temperature field and chemical abundances raises questions about the nature of the middle atmosphere's meridional circulation and evolution over Titan's 29-yr seasonal cycle that can only be answered through long-term monitoring campaigns. Here, we present maps of Titan's stratospheric temperature, acetonitrile (or methyl cyanide; CH$_3$CN), and monodeuterated methane (CH$_3$D) abundances following Titan's northern summer solstice obtained with Band 9 ($\sim0.43$ mm) ALMA observations. We find that increasing temperatures towards high-southern latitudes, currently in winter, resemble those observed during Titan's northern winter by the Cassini mission. Acetonitrile abundances have changed significantly since previous (sub)millimeter observations, and we find that the species is now highly concentrated at high-southern latitudes. The stratospheric CH$_3$D content is found to range between 4-8 ppm in these observations, and we infer the CH$_4$ abundance to vary between $\sim0.9-1.6\%$ through conversion with previously measured D/H values. A global value of CH$_4=1.15\%$ was retrieved, lending further evidence to the temporal and spatial variability of Titan's stratospheric methane when compared with previous measurements. Additional observations are required to determine the cause and magnitude of stratospheric enhancements in methane during these poorly understood seasons on Titan., Comment: Accepted for publication in the Planetary Science Journal. 9 Figures, 1 table

Saturn's moon Titan possesses stratospheric zonal winds that places it among a sparse class of planetary bodies known to have superrotation in their atmospheres. Few measurements have been made of these speeds in the upper stratosphere, leaving their seasonal variations still not well understood. We examined observations made with the extended Submillimeter Array (eSMA) in 2009 March (L\textsubscript{s}=355\textdegree{}) and 2010 February (L\textsubscript{s}=5\textdegree{}), shortly before and after Titan's northern spring equinox. Cassini observations and atmospheric models find equinoctial periods to be especially dynamic. Zonal wind calculations, derived from the Doppler frequency shift of CH\textsubscript{3}CN near 349.4 GHz, yielded speeds of 128 $\pm$ 27 m s\textsuperscript{-1} in 2009 and 209 $\pm$ 48 m s\textsuperscript{-1} in 2010. We estimated the measured emission to originate from vertical altitudes of $336^{+112}_{-88}$ kilometers, equivalent to pressures of $3.8^{+19.2}_{-3.4}$ Pa, commensurate with Titan's upper stratosphere/lower mesosphere. This suggests a possible increase in zonal speeds during this period. The results are then compared to those from previous Cassini-inferred and direct-interferometric observations of winds, as well as general circulation model simulations, to form a more complete picture of the seasonal cycle of stratospheric zonal winds., Comment: Accepted to the Planetary Science Journal; will be published 2024 April 12

The study of planets and small bodies within our Solar System is fundamental for understanding the formation and evolution the Earth and other planets. Compositional and meteorological studies of the giant planets provide a foundation for understanding the nature of the most commonly observed exoplanets, while spectroscopic observations of the atmospheres of terrestrial planets, moons, and comets provide insights into the past and present-day habitability of planetary environments, and the availability of the chemical ingredients for life. While prior and existing (sub)millimeter observations have led to major advances in these areas, progress is hindered by limitations in the dynamic range, spatial and temporal coverage, as well as sensitivity of existing telescopes and interferometers. Here, we summarize some of the key planetary science use cases that factor into the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST), a proposed 50-m class single dish facility: (1) to more fully characterize planetary wind fields and atmospheric thermal structures, (2) to measure the compositions of icy moon atmospheres and plumes, (3) to obtain detections of new, astrobiologically relevant gases and perform isotopic surveys of comets, and (4) to perform synergistic, temporally-resolved measurements in support of dedicated interplanetary space missions. The improved spatial coverage (several arcminutes), resolution ($\sim1.2''-12''$), bandwidth (several tens of GHz), dynamic range ($\sim10^5$) and sensitivity ($\sim1$ mK km s$^{-1}$) required by these science cases would enable new insights into the chemistry and physics of planetary environments, the origins of prebiotic molecules and the habitability of planetary systems in general., Comment: Submitted to Open Research Europe (AtLAST collection). 19 pages

Observations of linear polarization in the 2-8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating non-thermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy particles within the jet. IXPE observations of the X-ray selected BL Lac-type blazar 1ES 1959+650 in 2022 May 3-4 showed a significant linear polarization degree of $\Pi_\mathrm{x} = 8.0\% \pm 2.3\%$ at an electric-vector position angle $\psi_\mathrm{x} = 123^\circ \pm 8^\circ$. However, in 2022 June 9-12, only an upper limit of $\Pi_\mathrm{x} \leq 5.1\%$ could be derived (at the 99% confidence level). The degree of optical polarization at that time $\Pi_\mathrm{O} \sim 5\%$ is comparable to the X-ray measurement. We investigate possible scenarios for these findings, including temporal and geometrical depolarization effects. Unlike some other X-ray selected BL Lac objects, there is no significant chromatic dependence of the measured polarization in 1ES 1959+650, and its low X-ray polarization may be attributed to turbulence in the jet flow with dynamical timescales shorter than 1 day., Comment: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

X-ray continuum emission of active galactic nuclei (AGNs) may be reflected by circumnuclear dusty tori, producing prominent fluorescence iron lines at X-ray frequencies. Here we discuss the broad-band emission of three radio-loud AGN belonging to the class of compact symmetric objects (CSOs), with detected narrow Fe\,K$\alpha$ lines. CSOs have newly-born radio jets, forming compact radio lobes with projected linear sizes of the order of a few to hundreds of parsecs. We model the radio--to--$\gamma$-ray spectra of compact lobes in {J1407+2827}, J1511+0518, and {J2022+6137}, which are among the nearest and the youngest CSOs known to date, and are characterized by an intrinsic X-ray absorbing column density of $N_{\rm H} > 10^{23}$\,cm$^{-2}$. In addition to the archival data, we analyze the newly acquired \chandra\ X-ray Observatory and Sub-Millimeter Array (SMA) observations, and also refine the $\gamma$-ray upper limits from the \fermi\ Large Area Telescope (LAT) monitoring. The new \chandra\ data exclude the presence of the extended X-ray emission components on scales larger than $1.5^{\prime \prime}$. The SMA data unveil a correlation of the spectral index of the electron distribution in the lobes and $N_{\rm H}$, which can explain the $\gamma$-ray quietness of heavily obscured CSOs. Based on our modeling, we argue that the inverse-Compton emission of compact radio lobes may account for the intrinsic X-ray continuum in all these sources. Furthermore, we propose that the observed iron lines may be produced by a reflection of the lobes' continuum from the surrounding cold dust., Comment: The Astrophysical Journal, in press

Hyper-luminous infrared galaxies (HyLIRGs) are the most extreme star-forming systems observed in the early Universe, and their properties still elude comprehensive understanding. We have undertaken a large XMM-Newton observing program to probe the total accreting black hole population in three HyLIRGs at z = 2.12, 3.25, and 3.55, gravitationally lensed by foreground galaxies. Selected from the Planck All-Sky Survey to Analyze Gravitationally-lensed Extreme Starbursts (PASSAGES), these HyLIRGs have apparent infrared luminosities > E14 Lsun. Our observations revealed X-ray emission in each of them. PJ1336+49 appears to be dominated by high-mass X-ray binaries (HMXBs). Remarkably, the luminosity of this non-AGN X-ray emission exceeds by a factor of about three the value obtained by calibration with local galaxies with much lower star formation rates. This enhanced X-ray emission most likely highlights the efficacy of dynamical HMXB production within compact clusters, which is an important mode of star formation in HyLIRGs. The remaining two (PJ0116-24 and PJ1053+60) morphologically and spectrally exhibit a compact X-ray component in addition to the extended non-AGN X-ray emission, indicating the presence of Active Galactic Nuclei (AGNs). The AGN appears to be centrally located in the reconstructed source plane images of PJ0116-24, which manifests its star-forming activity predominantly within an extended galactic disk. In contrast, the AGN in the field of PJ1053+60 is projected 60 kpc away from the extreme star-forming galaxy and could be ejected from it. These results underline the synergistic potential of deep X-ray observations with strong lensing for the study of high-energy astrophysical phenomena in HyLIRGs., Comment: Just accepted for publication in MNRAS; comments are welcome

We present the multiwavelength flaring activity of the blazar AO 0235+164 during its recent active period from 2013 to 2019. From a discrete correlation function (DCF) analysis, we find a significant (>95%) correlation between radio and $\gamma$-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz VLBA data that was ejected from the radio core on MJD $57246^{+26}_{-30}$ (2015 August 12), during the peak of the 2015 radio flare. From the analysis of the jet component, we derived a Doppler factor of $\delta_{\rm var}=28.5\pm8.4$, a bulk Lorentz factor of $\Gamma=16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $\theta_{\rm v}=1.42^{+1.07}_{-0.52}\textrm{ degrees}$. Investigation of the quasi-simultaneous radio data revealed a partially absorbed spectrum with the turnover frequency varying in the range of $10-70$ GHz and the peak flux density varying in the range of $0.7-4$ Jy. We find the synchrotron self-absorption magnetic field strength to be $B_{\rm SSA}=15.3^{+12.6}_{-14.0}$ mG at the peak of the 2015 radio flare, which is comparable to the equipartition magnetic field strength of $B_{\rm EQ}=43.6^{+10.6}_{-10.4}$ mG calculated for the same epoch. Additional analysis of the radio emission region in the relativistic jet of AO 0235+164 suggests that it did not significantly deviate from equipartition during its recent flaring activity., Comment: 13 pages, 11 figures, 4 tables; Accepted for publication in MNRAS

We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements on 2022 December 6-8. We detected X-ray polarization of Mrk 421 with a degree of $\Pi_{\rm X}$=14$\pm$1$\%$ and an electric-vector position angle $\psi_{\rm X}$=107$\pm$3$^{\circ}$ in the 2-8 keV band. From the time variability analysis, we find a significant episodic variation in $\psi_{\rm X}$. During 7 months from the first IXPE pointing of Mrk 421 in 2022 May, $\psi_{\rm X}$ varied across the range of 0$^{\circ}$ to 180$^{\circ}$, while $\Pi_{\rm X}$ maintained similar values within $\sim$10-15$\%$. Furthermore, a swing in $\psi_{\rm X}$ in 2022 June was accompanied by simultaneous spectral variations. The results of the multiwavelength polarimetry show that the X-ray polarization degree was generally $\sim$2-3 times greater than that at longer wavelengths, while the polarization angle fluctuated. Additionally, based on radio, infrared, and optical polarimetry, we find that rotation of $\psi$ occurred in the opposite direction with respect to the rotation of $\psi_{\rm X}$ over longer timescales at similar epochs. The polarization behavior observed across multiple wavelengths is consistent with previous IXPE findings for HSP blazars. This result favors the energy-stratified shock model developed to explain variable emission in relativistic jets. The accompanying spectral variation during the $\psi_{\rm X}$ rotation can be explained by a fluctuation in the physical conditions, e.g., in the energy distribution of relativistic electrons. The opposite rotation direction of $\psi$ between the X-ray and longer-wavelength polarization accentuates the conclusion that the X-ray emitting region is spatially separated from that at longer wavelengths., Comment: 17 pages, 13 figures, 4 tables; Accepted for publication in A&A

We present polarization measurements in the $2-8 \thinspace \mathrm{keV}$ band from blazar 1ES 0229+200, the first extreme high synchrotron peaked source to be observed by the Imaging X-ray Polarimetry Explorer (IXPE). Combining two exposures separated by about two weeks, we find the degree of polarization to be $\Pi_{X} = 17.9 \pm 2.8 \%$ at an electric-vector position angle $\psi_X = 25.0 \pm 4.6^{\circ}$ using a spectropolarimetric fit from joint IXPE and XMM-Newton observations. There is no evidence for the polarization degree or angle varying significantly with energy or time on both short time scales (hours) or longer time scales (days). The contemporaneous polarization degree at optical wavelengths was $>$7$\times$ lower, making 1ES 0229+200 the most strongly chromatic blazar yet observed. This high X-ray polarization compared to the optical provides further support that X-ray emission in high-peaked blazars originates in shock-accelerated, energy-stratified electron populations, but is in tension with many recent modeling efforts attempting to reproduce the spectral energy distribution of 1ES 0229+200 which attribute the extremely high energy synchrotron and Compton peaks to Fermi acceleration in the vicinity of strongly turbulent magnetic fields., Comment: 17 Pages, 6 Figures, Resubmitted to ApJ after addressing referee comments

The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($\lambda$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region., Comment: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

We present a thermal observation of Callisto's leading hemisphere obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) at 0.87 mm (343 GHz). The angular resolution achieved for this observation was $\sim$$0.16^{\prime\prime}$, which for Callisto at the time of this observation ($D\sim 1.05^{\prime\prime}$) was equivalent to $\sim$6 elements across the surface. Our disk-integrated brightness temperature of 116 $\pm$ 5 K (8.03 $\pm$ 0.40 Jy) is consistent with prior disk-integrated observations. Global surface properties were derived from the observation using a thermophysical model (de Kleer et al. 2021) constrained by spacecraft data. We find that models parameterized by two thermal inertia components more accurately fit the data than single thermal inertia models. Our best-fit global parameters adopt a lower thermal inertia of 15-50 $\text{J}\:\text{m}^{-2}\:\text{K}^{-1}\:\text{s}^{-1/2}$ and a higher thermal inertia component of 1200-2000 $\text{J}\:\text{m}^{-2}\:\text{K}^{-1}\:\text{s}^{-1/2}$, with retrieved millimeter emissivities of 0.89-0.91. We identify several thermally anomalous regions, including spots $\sim$3 K colder than model predictions co-located with the Valhalla impact basin and a complex of craters in the southern hemisphere; this indicates the presence of materials possessing either a higher thermal inertia or a lower emissivity. A warm region confined to the mid-latitudes in these leading hemisphere data may be indicative of regolith property changes due to exogenic sculpting., Comment: 13 pages, 7 figures

The lower energy peak of the spectral energy distribution of blazars has commonly been ascribed to synchrotron radiation from relativistic particles in the jets. Despite the consensus regarding jet emission processes, the particle acceleration mechanism is still debated. Here, we present the first X-ray polarization observations of PG 1553+113, a high-synchrotron-peak blazar observed by the Imaging X-ray Polarimetry Explorer (IXPE). We detect an X-ray polarization degree of $(10\pm2)\%$ along an electric-vector position angle of $\psi_X=86^{\circ}\pm8^{\circ}$. At the same time, the radio and optical polarization degrees are lower by a factor of $\sim$3. During our IXPE pointing, we observed the first orphan optical polarization swing of the IXPE era, as the optical angle of PG 1553+113 underwent a smooth monotonic rotation by about 125$^\circ$, with a rate of $\sim$17 degrees per day. We do not find evidence of a similar rotation in either radio or X-rays, which suggests that the X-ray and optically emitting regions are separate or, at most, partially co-spatial. Our spectro-polarimetric results provide further evidence that the steady-state X-ray emission in blazars originates in a shock-accelerated and energy-stratified electron population., Comment: Accepted for pubblication in The Astrophysical Journal Letters. 3 Figures, 13 pages

We present a multi-wavelength analysis using SMA, JCMT, NOEMA, JWST, HST, and SST of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A1489. These SMA-located sources both lie at z=4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at <2um, even with deep JWST/NIRCam observations. We investigate their stellar, ISM, and dynamical properties, including a pixel-level SED analysis to derive sub-kpc-resolution stellar-mass and Av maps. We find that 850.1 is one of the most massive and highly obscured, Av~5, galaxies known at z>4 with M*~10^11.8 Mo (likely forming at z>6), and 850.2 is one of the least massive and least obscured, Av~1, members of the z>4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z>3-4 based on imaging at <2um, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar mass distribution, Re(mass)~1kpc, but its expected evolution to z~1.5 and then z~0 matches both the properties of massive, quiescent galaxies at z~1.5 and ultra-massive early-type galaxies at z~0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature which we link to the active dynamical environment it resides within., Comment: Submitted to ApJ, comments welcome!

We present 1.3 mm (230 GHz) observations of the recent and nearby Type II supernova, SN2023ixf, obtained with the Submillimeter Array (SMA) at 2.6-18.6 days after explosion. The observations were obtained as part the SMA Large Program POETS (Pursuit of Extragalactic Transients with the SMA). We do not detect any emission at the location of SN2023ixf, with the deepest limits of $L_\nu(230\,{\rm GHz})\lesssim 8.6\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 2.7 and 7.7 days, and $L_\nu(230\,{\rm GHz})\lesssim 3.4\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 18.6 days. These limits are about a factor of 2 times dimmer than the mm emission from SN2011dh (IIb), about an order of magnitude dimmer compared to SN1993J (IIb) and SN2018ivc (IIL), and about 30 times dimmer than the most luminous non-relativistic SNe in the mm-band (Type IIb/Ib/Ic). Using these limits in the context of analytical models that include synchrotron self-absorption and free-free absorption we place constraints on the proximate circumstellar medium around the progenitor star, to a scale of $\sim 2\times 10^{15}$ cm, excluding the range $\dot{M}\sim {\rm few}\times 10^{-6}-10^{-2}$ M$_\odot$ yr$^{-1}$ (for a wind velocity, $v_w=115$ km s$^{-1}$, and ejecta velocity, $v_{\rm eje}\sim (1-2)\times 10^4$ km s$^{-1}$). These results are consistent with an inference of the mass loss rate based on optical spectroscopy ($\sim 2\times 10^{-2}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$), but are in tension with the inference from hard X-rays ($\sim 7\times 10^{-4}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$). This tension may be alleviated by a non-homogeneous and confined CSM, consistent with results from high-resolution optical spectroscopy., Comment: Submitted

We present the results of a radio multi-frequency ($\rm 3-340~GHz$) study of the blazar 3C~454.3. After subtracting the quiescent spectrum corresponding to optically thin emission, we found two individual synchrotron self-absorption (SSA) features in the wide-band spectrum. The one SSA had a relatively low turnover frequency ($\nu_{\rm m}$) in the range of $\rm 3-37~GHz$ (lower $\nu_{\rm m}$ SSA spectrum, LSS), and the other one had a relatively high $\nu_{\rm m}$ of $\rm 55-124~GHz$ (higher $\nu_{\rm m}$ SSA spectrum, HSS). Using the SSA parameters, we estimated magnetic field strengths at the surface where optical depth $\tau=1$. The estimated magnetic field strengths were $\rm >7~mG$ and $\rm >0.2~mG$ for the LSS and HSS, respectively. The LSS emitting region was magnetically dominated before the June 2014 $\gamma$-ray flare. The quasi-stationary component (C), $\sim 0.6~{\rm mas}$ apart from the 43 GHz radio core, became brighter than the core with decreasing observing frequency, and we found that component C was related to the LSS. A decrease in jet width was found near component C. As a moving component, K14 approached component C, and the flux density of the component was enhanced while the angular size decreased. The high intrinsic brightness temperature in the fluid frame was obtained as $T_{\rm B, int} \approx (7.0\pm1.0) \times 10^{11}~{\rm K}$ from the jet component after the 2015 August $\gamma$-ray flare, suggesting that component C is a high-energy emitting region. The observed local minimum of jet width and re-brightening behavior suggest a possible recollimation shock in component C., Comment: 16 pages, 12 figures, 5 tables, Accepted for publication in MNRAS

We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$\gamma$-ray outburst, a 287 ksec observation (2022 November 27-30) was taken using the Imaging X-ray Polarimetry Explorer ({\it IXPE}), together with contemporaneous multiwavelength observations from the Neil Gehrels {\it Swift} observatory and {\it XMM-Newton} in soft X-rays (0.3--10~keV), {\it NuSTAR} in hard X-rays (3--70~keV), and optical polarization from the Calar Alto, and Perkins Telescope observatories. Our contemporaneous X-ray data suggest that the {\it IXPE} energy band is at the crossover between the low- and high-frequency blazar emission humps. The source displays significant variability during the observation, and we measure polarization in three separate time bins. Contemporaneous X-ray spectra allow us to determine the relative contribution from each emission hump. We find $>99\%$ confidence X-ray polarization $\Pi_{2-4{\rm keV}} = 21.7^{+5.6}_{-7.9}\%$ and electric vector polarization angle $\psi_{2-4{\rm keV}} = -28.7 \pm 8.7^{\circ}$ in the time bin with highest estimated synchrotron flux contribution. We discuss possible implications of our observations, including previous {\it IXPE} BL Lacertae pointings, tentatively concluding that synchrotron self-Compton emission dominates over hadronic emission processes during the observed epochs., Comment: Accepted to ApJL. 16 pages, 8 figures

The magnetic field conditions in astrophysical relativistic jets can be probed by multiwavelength polarimetry, which has been recently extended to X-rays. For example, one can track how the magnetic field changes in the flow of the radiating particles by observing rotations of the electric vector position angle $\Psi$. Here we report the discovery of a $\Psi_{\mathrm x}$ rotation in the X-ray band in the blazar Mrk 421 at an average flux state. Across the 5 days of Imaging X-ray Polarimetry Explorer (IXPE) observations of 4-6 and 7-9 June 2022, $\Psi_{\mathrm x}$ rotated in total by $\geq360^\circ$. Over the two respective date ranges, we find constant, within uncertainties, rotation rates ($80 \pm 9$ and $91 \pm 8 ^\circ/\rm day$) and polarization degrees ($\Pi_{\mathrm x}=10\%\pm1\%$). Simulations of a random walk of the polarization vector indicate that it is unlikely that such rotation(s) are produced by a stochastic process. The X-ray emitting site does not completely overlap the radio/infrared/optical emission sites, as no similar rotation of $\Psi$ was observed in quasi-simultaneous data at longer wavelengths. We propose that the observed rotation was caused by a helical magnetic structure in the jet, illuminated in the X-rays by a localized shock propagating along this helix. The optically emitting region likely lies in a sheath surrounding an inner spine where the X-ray radiation is released.

The variability of submillimeter emission provides a useful tool to probe the accretion physics in low-luminosity active galactic nuclei. We accumulate four years of observations using Submillimeter Array for Centaurus A, NGC 4374, NGC 4278, and NGC 5077 and one year of observations for NGC 4552 and NGC 4579. All sources are variable. We measure the characteristic timescale at which the variability is saturated by modeling these sources' light curve as a damped random walk. We detect a timescale for all the sources except NGC 4552. The detected timescales are comparable to the orbital timescale at the event horizon scale for most sources. Combined with previous studies, we show a correlation between the the timescale and the black hole mass over three orders of magnitude. This discovery suggests the sub-mm emission is optically thin with the emission originating from the event horizon. The mass scaling relationship further suggests that a group of radio sources with a broadband spectrum that peaks at submillimeter wavelengths have similar inner accretion physics. Sources that follow this relationship may be good targets for high-resolution imaging with the Event Horizon Telescope., Comment: 18 pages, 7 figures, 4 tables, accepted for publication by ApJ

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties., Comment: Accepted for publication in ApJ

The accretion of matter by compact objects can be inhibited by radiation pressure if the luminosity exceeds the critical value, known as the Eddington limit. Discovery of ultraluminous X-ray sources has shown that accretion can proceed even when the apparent luminosity significantly exceeds this limit. High apparent luminosity might be produced thanks to geometric beaming of the radiation by an outflow. The outflow half-opening angle, which determines the amplification due to beaming, has never been robustly constrained. Using the Imaging X-ray Polarimetry Explorer, we made the measurement of X-ray polarization in the Galactic X-ray binary Cyg X-3. We find high, over 20%, nearly energy-independent linear polarization, orthogonal to the direction of the radio ejections. These properties unambiguously indicate the presence of a collimating outflow in the X-ray binary Cyg~X-3 and constrain its half-opening angle, <15 degrees. Thus, the source can be used as a laboratory for studying the super-critical accretion regime. This finding underscores the importance of X-ray polarimetry in advancing our understanding of accreting sources., Comment: 52 pages, 14 figures, 3 tables. Author version of the article published in Nature Astronomy

We report on the observations of the quasar NRAO 530 with the Event Horizon Telescope (EHT) on 2017 April 5-7, when NRAO 530 was used as a calibrator for the EHT observations of Sagittarius A*. At z=0.902 this is the most distant object imaged by the EHT so far. We reconstruct the first images of the source at 230 GHz, at an unprecedented angular resolution of $\sim$ 20 $\mu$as, both in total intensity and in linear polarization. We do not detect source variability, allowing us to represent the whole data set with static images. The images reveal a bright feature located on the southern end of the jet, which we associate with the core. The feature is linearly polarized, with a fractional polarization of $\sim$5-8% and has a sub-structure consisting of two components. Their observed brightness temperature suggests that the energy density of the jet is dominated by the magnetic field. The jet extends over 60 $\mu$as along a position angle PA$\sim -$28$^\circ$. It includes two features with orthogonal directions of polarization (electric vector position angle, EVPA), parallel and perpendicular to the jet axis, consistent with a helical structure of the magnetic field in the jet. The outermost feature has a particularly high degree of linear polarization, suggestive of a nearly uniform magnetic field. Future EHT observations will probe the variability of the jet structure on ${\mu}$as scales, while simultaneous multi-wavelength monitoring will provide insight into the high energy emission origin.

Blazars are a class of jet-dominated active galactic nuclei with a typical double-humped spectral energy distribution. It is of common consensus the Synchrotron emission to be responsible for the low frequency peak, while the origin of the high frequency hump is still debated. The analysis of X-rays and their polarization can provide a valuable tool to understand the physical mechanisms responsible for the origin of high-energy emission of blazars. We report the first observations of BL Lacertae performed with the Imaging X-ray Polarimetry Explorer ({IXPE}), from which an upper limit to the polarization degree $\Pi_X<$12.6\% was found in the 2-8 keV band. We contemporaneously measured the polarization in radio, infrared, and optical wavelengths. Our multiwavelength polarization analysis disfavors a significant contribution of proton synchrotron radiation to the X-ray emission at these epochs. Instead, it supports a leptonic origin for the X-ray emission in BL Lac., Comment: 17 pages, 5 figures, accepted for publication in ApJL

The blazar J1924-2914 is a primary Event Horizon Telescope (EHT) calibrator for the Galactic Center's black hole Sagittarius A*. Here we present the first total and linearly polarized intensity images of this source obtained with the unprecedented 20 $\mu$as resolution of the EHT. J1924-2914 is a very compact flat-spectrum radio source with strong optical variability and polarization. In April 2017 the source was observed quasi-simultaneously with the EHT (April 5-11), the Global Millimeter VLBI Array (April 3), and the Very Long Baseline Array (April 28), giving a novel view of the source at four observing frequencies, 230, 86, 8.7, and 2.3 GHz. These observations probe jet properties from the subparsec to 100-parsec scales. We combine the multi-frequency images of J1924-2914 to study the source morphology. We find that the jet exhibits a characteristic bending, with a gradual clockwise rotation of the jet projected position angle of about 90 degrees between 2.3 and 230 GHz. Linearly polarized intensity images of J1924-2914 with the extremely fine resolution of the EHT provide evidence for ordered toroidal magnetic fields in the blazar compact core.

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

The nearby Seyfert type galaxy NGC 1275 contains a bright radio nucleus at its center, revealed through high-spatial resolution imaging to be the source of the jets emanating from the galaxy. Coincident with the emergence of a new component C3 in the nucleus since 2005, flux densities from NGC 1275, at least at radio, millimeter (mm), and gamma-ray frequencies, had been increasing up through 2017 and leveled off afterwards. We analyze the long-term light curves of the nucleus that span the rising trend to 2015 July, and find a pair of approximately year-long quasi-periodic oscillations, with periods of $P_l\simeq 345$\,d and $P_h\simeq 386$\,d respectively, in emission at 1.3-mm wavelength. We discuss the case that there would be a long precession period $P_{\rm prec}\simeq 9$\,yr, causing the appearance of $P_h$ that is slightly higher than $P_l$. The accretion disk around the central supermassive black hole (SMBH) would be precessing at $P_{\rm prec}$, induced by either the Lense-Thirring effect or the existence of a companion SMBH. In the two scenarios, $P_l$ would be the jet wobbling timescale or the SMBH binary period respectively. The finding, which could be verified through high-spatial resolution mm imaging, would not only identify the nature of the jet variation but also help reveal the full features of the galaxy., Comment: 12 pages, 8 figures, accepted for publication in ApJ

The origin of gamma-ray flares observed from blazars is one of the major mysteries in jet physics. We have attempted to address this problem following a novel spectral energy distribution (SED) fitting technique that explored the flaring patterns identified in the broadband SEDs of two gamma-ray bright blazars 3C 279 (z=0.54) and 3C 454.3 (z=0.86), using near-simultaneous radio-to-gamma-ray observations. For both sources, the gamma-ray flux strongly correlates with the separation of the SED peaks and the Compton dominance. We propose that spectral hardening of the radiating electron population and/or enhancement of the Doppler factor can naturally explain these observations. In both cases, magnetic reconnection may play a pivotal role in powering the luminous gamma-ray flares., Comment: ApJS, in press

Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below ${\lambda}1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at ${\lambda}1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses., Comment: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

We explore the connection between the $\gamma$-ray and radio emission in the jet of the blazar 0716$+$714 by using 15, 37, and 230 GHz radio and 0.1$-$200 GeV $\gamma$-ray light curves spanning 10.5 years (2008$-$2019). We find significant positive and negative correlations between radio and $\gamma$-ray fluxes in different time ranges. The time delays between radio and $\gamma$-ray emission suggest that the observed $\gamma$-ray flares originated from multiple regions upstream of the radio core, within a few parsecs from the central engine. Using time-resolved 43 GHz VLBA maps we identified 14 jet components moving downstream along the jet. Their apparent speeds range from 6 to 26 $c$, showing notable variations in their position angles upstream the stationary component ($\sim$0.53 mas from the core). The brightness temperature declines as function of distance from the core according to a power-law which becomes shallower at the location of the stationary component. We also find that the periods at which significant correlations between radio and $\gamma$-ray emission occur overlap with the times when the jet was oriented to the north. Our results indicate that the passage of a propagating disturbance (or shock) through the radio core and the orientation of the jet might be responsible for the observed correlation between the radio and $\gamma$-ray variability. We present a scenario that connects the positive correlation and the unusual anti-correlation by combining the production of a flare and a dip at $\gamma$-rays by a strong moving shock at different distances from the jet apex., Comment: 26 pages, 13 figures, 5 tables, Accepted for publication in ApJ

The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas., Comment: Accepted for Publication in ApJ

We report multi-epoch Very Long Baseline Interferometric (VLBI) observations of the compact radio lobe in the radio galaxy 3C84 (NGC1275) during 2016 - 2020. The image sequence of 3C84 reveals that the hotspot in the radio lobe showed the one-year long frustration in 2017 within a compact region of about 0.07 parsec, suggesting a strong collision between the jet and a compact dense cloud with the estimated average density about $(4\-- 6)\times 10^{5}~{\rm cm^{-3}}$. Although the hotspot and the radio lobe began to move south again after its breakout, the radio lobe showed a morphological transition from FR II- to FR I-class radio lobe and its radio flux became fainter. This is the first detection of the dynamical feedback from the cloud to the jet where the cloud located on the jet axis significantly interferes with the jet propagation and evolution at the central one-parsec region in active galactic nucleus., Comment: The Astrophysical Journal Letters, Volume 920, Issue 1, id.L24, 7 pp

Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes., Comment: 15 pages, 4 figures, published in PRD on May 19

Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images., Comment: 29 pages, 11 figures, published in ApJ on May 3

We present thermal observations of Ganymede from the Atacama Large Millimeter Array (ALMA) in 2016-2019 at a spatial resolution of 300-900 km (0.1-0.2'' angular resolution) and frequencies of 97.5, 233, and 343.5 GHz (wavelengths of 3, 1.3, and 0.87 mm); the observations collectively covered all Ganymede longitudes. We determine the global thermophysical properties using a thermal model that considers subsurface emission and depth- and temperature-dependent thermophysical and dielectric properties, in combination with a retrieval algorithm. The data are sensitive to emission from the upper $\sim$0.5 meter of the surface, and we find a millimeter emissivity of 0.75-0.78 and (sub)surface porosities of 10-40%, corresponding to effective thermal inertias of 400-800 J m^{-2} K^{-1} s^{-1/2}. Combined with past infrared results, as well as modeling presented here of a previously-unpublished Galileo PPR nighttime infrared observation, the multi-wavelength constraints are consistent with a compaction profile whereby the porosity drops from ~85% at the surface to 10{+30/-10}% at depth over a compaction length scale of tens of cm. We present maps of temperature residuals from the best-fit global models which indicate localized variations in thermal surface properties at some (but not all) dark terrains and at impact craters, which appear 5-8 K colder than the model. Equatorial regions are warmer than predicted by the model, in particular near the centers of the leading and trailing hemispheres, while the mid-latitudes (~30-60 degrees) are generally colder than predicted; these trends are suggestive of an exogenic origin., Comment: Fixed typo in Eqn 1; article published in Planetary Science Journal, Vol 2, Issue 1, id.5 (Feb 2021)

The detection of phosphine (PH3) in the atmosphere of Venus has been recently reported based on millimeter-wave radio observations (Greaves et al. 2020), and its re-analyses (Greaves et al. 2021a/b). In this Matters Arising we perform an independent reanalysis, identifying several issues in the interpretation of the spectroscopic data. As a result, we determine sensitive upper-limits for PH3 in Venus' atmosphere (>75 km, above the cloud decks) that are discrepant with the findings in G2020 and G2021a/b. The measurements target the fundamental first rotational transition of PH3 (J=1-0) at 266.944513 GHz, which was observed with the James Clerk Maxwell Telescope (JCMT) in June 2017 and with the Atacama Large Millimeter/submillimeter Array (ALMA) in March 2019. This line's center is near the SO2 (J=309,21-318,24) transition at 266.943329 GHz (only 1.3 km/s away from the PH3 line) which represents a potential source of contamination. The JCMT and ALMA data, as presented in G2020, are at spectral resolutions comparable to the frequency separation of the two lines. Moreover, the spectral features identified are several km/s in width, and therefore do not permit distinct spectroscopic separation of the candidate spectral lines of PH3 and SO2. We present the radiative transfer modelling we have performed and then discuss the ALMA and JCMT analyses in turn.

The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources., Comment: Physical Review Letters

We present the results of simultaneous multi-frequency imaging observations at 22, 43, 86, and 129\,GHz of OJ\,287. We used the Korean VLBI Network as part of the Interferometric MOnitoring of GAmma-ray Bright active galactic nuclei (iMOGABA). The iMOGABA observations were performed during 31 epochs from 2013 January 16 to 2016 December 28. We also used 15\,GHz OVRO and 225\,GHz SMA flux density data. We analyzed four flux enhancements in the light curves. The estimated time scales of three flux enhancements were similar with time scales of $\sim$50 days at two frequencies. A fourth flux enhancement had a variability timescale approximately twice as long. We found that 225\,GHz enhancements led the 15\,GHz enhancements by a range of 7 to 30 days in the time delay analysis. We found the fractional variability did not change with frequency between 43 and 86\,GHz. We could reliably measure the turnover frequency, $\nu_{\rm c}$, of the core of the source in three epochs. This was measured to be in a range from 27 to 50\,GHz and a flux density at the turnover frequency, $S_{\rm m}$, ranging from 3-6\,Jy. The derived SSA magnetic fields, $B_{\rm SSA}$, are in a range from $0.157\pm0.104$ to $0.255\pm0.146$ mG. We estimated the equipartition magnetic field strengths to be in a range from $0.95\pm0.15$ to $1.93\pm0.30$ mG. The equipartition magnetic field strengths are up to a factor of 10 higher than the values of $B_{\rm SSA}$. We conclude that the downstream jet may be more particle energy dominated., Comment: 27 pages, 9 figures, Accepted for publication in ApJ

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/., Comment: Accepted to ApJS

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

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