Density Profiles of Collapsed Rotating Massive Stars Favor Long Gamma-Ray Bursts by Goni Halevi et al. on Wednesday 23 November

Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated

outflows -- jets -- that drill their way out of collapsing massive stars.

Accurately modeling this process requires realistic stellar profiles for the

jets to propagate through and break out of. Most previous studies have used

simple power laws or pre-collapse models for massive stars. However, the

relevant stellar profile for lGRB models is in fact that of a star after its

core has collapsed to form a compact object. To self-consistently compute such

a stellar profile, we use the open-source code GR1D to simulate the

core-collapse process for a suite of low-metallicity, rotating, massive stellar

progenitors that have undergone chemically homogeneous evolution. Our models

span a range of zero-age main sequence (ZAMS) masses: $M_\mathrm{ZAMS} = 13,

18, 21, 25, 35, 40$, and $45 M_\odot$. All of these models, at the onset of

core-collapse, feature steep density profiles, $\rho \propto r^{-\alpha}$ with

$\alpha\approx 2.5$, which would result in jets that are inconsistent with lGRB

observables. We follow the collapse of four out of our seven models until they

form BHs and the other three proto-neutron stars (PNSs). We find, across all

models, that the density profile outside of the newly-formed BH or PNS is

well-represented by a flatter power law with $\alpha \approx 1.35{-}1.55$. Such

flat density profiles are conducive to successful formation and breakout of

BH-powered jets and, in fact, required to reproduce observable properties of

lGRBs. Future models of lGRBs should be initialized with shallower

\textit{post-collapse} stellar profiles like those presented here instead of

the much steeper pre-collapse profiles that are typically used.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11781v1

The dark side of galaxy stellar populations II: The dependence of star formation histories on halo mass and on the scatter of the main sequence by Laura Scholz-Diaz et al. on Wednesday 23 November

Nearby galaxies are the end result of their cosmological evolution, which is

predicted to be influenced by the growth of their host dark matter halos. This

co-evolution potentially leaves signatures in present-day observed galaxy

properties, which might be essential to further understand how the growth and

properties of galaxies are connected to those of their host halos. In this

work, we study the evolutionary histories of nearby galaxies both in terms of

their host halos and the scatter of the star-forming main sequence by

investigating their time-resolved stellar populations using absorption optical

spectra drawn from the Sloan Digital Sky Survey. We find that galaxy star

formation histories depend on the masses of their host halos, and hence they

shape the evolution of the star-forming main sequence over cosmic time.

Additionally, we also find that the scatter around the z=0 star-forming main

sequence is not (entirely) stochastic, as galaxies with currently different

star formation rates have experienced, on average, different star formation

histories. Our findings suggest that dark matter halos might play a key role in

modulating the evolution of star formation in galaxies, and thus of the main

sequence, and further demonstrate that galaxies at different evolutionary

stages contribute to the observed scatter of this relation.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11779v1

A dichotomy in group II Herbig disks: ALMA gas disk height measurements show both shadowed large vertically extended disks and compact flat disks by L. M. Stapper et al. on Wednesday 23 November

Herbig stars can be classified into group I and group II depending on the

shape of the far-IR excess from the spectral energy distribution. This

separation may be evolutionary and related to the vertical structure of these

disks. We aim to determine the emission height of Herbig disks and compare the

resulting vertical extent of both groups. ALMA Band 6 observations of 12CO

emission lines at sufficient velocity and spatial resolution of eight Herbig

disks (four group I and four group II sources) are used to determine the

emission heights from the channel maps via geometrical methods developed in

other works. We find that all group I disks are vertically extended with a

height to radius ratio of at least 0.25, and for three of the disks the gas

emission profile can be traced out to 200-500 au. The group II disks are

divided between MWC 480 and HD 163296 which have similar emission height

profiles as the group I disks, and AK Sco and HD 142666 which are very flat

(not exceeding a height of 10 au) and more compact (<200 au in size). The

brightness temperatures show no differences between the disks when the

luminosity of the host star is accounted for. Our findings agree with previous

work suggesting that group I disks are vertically extended and that group II

disks are either large and self-shadowed or compact. Both MWC 480 and HD 163296

could be precursors of group I disks, which we see now before a cavity has

formed that would allow irradiation of the outer parts of the disk. The very

flat disks AK Sco and HD 142666 could be due to significant settling because of

the advanced age of these disks (~20 instead of <10 Myr). These large

differences in vertical structures are not reflected in the spectral energy

distributions of these disks. More and deeper observations at higher spatial

and velocity resolution are necessary to further characterize the Herbig

sub-groups.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12532v1

EAVN Astrometry toward the Extreme Outer Galaxy: Kinematic distance with the proper motion of G034 84-00 95 by Nobuyuki Sakai et al. on Wednesday 23 November

We aim to reveal the structure and kinematics of the Outer-Scutum-Centaurus

(OSC) arm located on the far side of the Milky Way through very long baseline

interferometry (VLBI) astrometry using KaVA, which is composed of KVN (Korean

VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We report the

proper motion of a 22 GHz H$_{2}$O maser source, which is associated with the

star-forming region G034.84$-$00.95, to be ($\mu_{\alpha} \rm{cos}\delta$,

$\mu_{\delta}$) = ($-$1.61$\pm$0.18, $-$4.29$\pm$0.16) mas yr$^{-1}$ in

equatorial coordinates (J2000). We estimate the 2D kinematic distance to the

source to be 18.6$\pm$1.0 kpc, which is derived from the variance-weighted

average of kinematic distances with LSR velocity and the Galactic-longitude

component of the measured proper motion. Our result places the source in the

OSC arm and implies that G034.84$-$00.95 is moving away from the Galactic plane

with a vertical velocity of $-$38$\pm$16 km s$^{-1}$. Since the H I supershell

GS033+06$-$49 is located at a kinematic distance roughly equal to that of

G034.84$-$00.95, it is expected that gas circulation occurs between the outer

Galactic disk around G034.84$-$00.95 with a Galactocentric distance of

12.8$^{+1.0}_{-0.9}$ kpc and halo. We evaluate possible origins of the fast

vertical motion of G034.84$-$00.95, which are (1) supernova explosions and (2)

cloud collisions with the Galactic disk. However, neither of the possibilities

are matched with the results of VLBI astrometry as well as spatial

distributions of H II regions and H I gas.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12534v1

Radiative Transfer as a Bayesian Linear Regression problem by Frederik De Ceuster et al. on Wednesday 23 November

Electromagnetic radiation plays a crucial role in various physical and

chemical processes. Hence, almost all astrophysical simulations require some

form of radiative transfer model. Despite many innovations in radiative

transfer algorithms and their implementation, realistic radiative transfer

models remain very computationally expensive, such that one often has to resort

to approximate descriptions. The complexity of these models makes it difficult

to assess the validity of any approximation and to quantify uncertainties on

the model results. This impedes scientific rigour, in particular, when

comparing models to observations, or when using their results as input for

other models. We present a probabilistic numerical approach to address these

issues by treating radiative transfer as a Bayesian linear regression problem.

This allows us to model uncertainties on the input and output of the model with

the variances of the associated probability distributions. Furthermore, this

approach naturally allows us to create reduced-order radiative transfer models

with a quantifiable accuracy. These are approximate solutions to exact

radiative transfer models, in contrast to the exact solutions to approximate

models that are often used. As a first demonstration, we derive a probabilistic

version of the method of characteristics, a commonly-used technique to solve

radiative transfer problems.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12547v1

Quasar Factor Analysis -- An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis by Zechang Sun et al. on Wednesday 23 November

Since their first discovery, quasars have been essential probes of the

distant Universe. However, due to our limited knowledge of its nature,

predicting the intrinsic quasar continua has bottlenecked their usage. Existing

methods of quasar continuum recovery often rely on a limited number of

high-quality quasar spectra, which might not capture the full diversity of the

quasar population. In this study, we propose an unsupervised probabilistic

model, \textit{Quasar Factor Analysis} (QFA), which combines factor analysis

(FA) with physical priors of the intergalactic medium (IGM) to overcome these

limitations. QFA captures the posterior distribution of quasar continua through

generatively modeling quasar spectra. We demonstrate that QFA can achieve the

state-of-the-art performance, $\sim 2\%$ relative error, for continuum

prediction in the Ly$\alpha$ forest region compared to previous methods. We

further fit 90,678 $2<\mathrm{z}<3.5$, SNR$>2$ quasar spectra from Sloan

Digital Sky Survey Data Release 16 and found that for $\sim 30\%$ quasar

spectra where the continua were ill-determined with previous methods, QFA

yields visually more plausible continua. QFA also attains $\lesssim 1\%$ error

in the 1D Ly$\alpha$ power spectrum measurements at $\mathrm{z}\sim 3$ and

$\sim 4\%$ in $\mathrm{z}\sim 2.4$. In addition, QFA determines latent factors

representing more physically motivated than PCA. We investigate the evolution

of the latent factors and report no significant redshift or luminosity

dependency except for the Baldwin effect. The generative nature of QFA also

enables outlier detection robustly; we showed that QFA is effective in

selecting outlying quasar spectra, including damped Ly$\alpha$ systems and

potential Type II quasar spectra.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11784v1

The production of ionizing photons in UV-faint z~3-7 galaxies by Gonzalo Prieto-Lyon et al. on Wednesday 23 November

The demographics of the production and escape of ionizing photons from

UV-faint early galaxies is a key unknown in discovering the primary drivers of

reionization. With the advent of JWST it is finally possible to observe the

rest-frame optical nebular emission from individual sub-L$^*$ z>3 galaxies to

measure the production of ionizing photons, $\xi_\mathrm{ion}$. Here we study a

sample of 380 z~3-7 galaxies spanning -23 0.5) in our sample, based on their photometric properties,

even amongst the Lyman-alpha selected galaxies. This work demonstrates that

considering the full distribution of $\xi_\mathrm{ion}$ across galaxy

properties is important for assessing the primary drivers of reionization.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12548v1

Reconstructing and Classifying SDSS DR16 Galaxy Spectra with Machine-Learning and Dimensionality Reduction Algorithms by Felix Pat et al. on Wednesday 23 November

Optical spectra of galaxies and quasars from large cosmological surveys are

used to measure redshifts and infer distances. They are also rich with

information on the intrinsic properties of these astronomical objects. However,

their physical interpretation can be challenging due to the substantial number

of degrees of freedom, various sources of noise, and degeneracies between

physical parameters that cause similar spectral characteristics. To gain deeper

insights into these degeneracies, we apply two unsupervised machine learning

frameworks to a sample from the Sloan Digital Sky Survey data release 16 (SDSS

DR16). The first framework is a Probabilistic Auto-Encoder (PAE), a two-stage

deep learning framework consisting of a data compression stage from 1000

elements to 10 parameters and a density estimation stage. The second framework

is a Uniform Manifold Approximation and Projection (UMAP), which we apply to

both the uncompressed and compressed data. Exploring across regions on the

compressed data UMAP, we construct sequences of stacked spectra which show a

gradual transition from star-forming galaxies with narrow emission lines and

blue spectra to passive galaxies with absorption lines and red spectra.

Focusing on galaxies with broad emission lines produced by quasars, we find a

sequence with varying levels of obscuration caused by cosmic dust. The

experiments we present here inform future applications of neural networks and

dimensionality reduction algorithms for large astronomical spectroscopic

surveys.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11783v1

Wideband timing of the Parkes Pulsar Timing Array UWL data by Małgorzata Curyło et al. on Wednesday 23 November

In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on

the 64-m Parkes Radio Telescope enabling observations with an instantaneous

frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a

three-year data set of 35 millisecond pulsars observed with the UWL by the

Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key

differences compared to typical narrow-band methods are, firstly, generation of

two-dimensional templates accounting for pulse shape evolution with frequency

and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and

dispersion measure (DM). This is the first time that wideband timing has been

applied to a uniform data set collected with a single large-fractional

bandwidth receiver, for which such techniques were originally developed. As a

result of our study, we present a set of profile evolution models and new

timing solutions including initial noise analysis. Precision of our ToA and DM

measurements is in the range of 0.005 $-$ 2.08 $\mu$s and

(0.043$-$14.24)$\times10^{-4}$ cm$^{-3}$ pc, respectively, with 94% of the

pulsars achieving a median ToA uncertainty of less than 1 $\mu$s.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12924v1

MeerKAT follow-up of enigmatic GLEAM 4-Jy G4Jy sources by Precious K. Sejake et al. on Wednesday 23 November

We present the results from studying 140 radio sources in the GLEAM (GaLactic

and Extragalactic All-sky MWA [Murchison Widefield Array]) 4-Jy (G4Jy) Sample.

These sources were followed-up with MeerKAT to assess their radio morphology

and enable host-galaxy identification, as existing radio images of 25 to

45-arcsec resolution do not provide sufficient information. We refer to these

sources as the MeerKAT-2019 subset. The aim is to identify the host galaxy of

these sources by visually inspecting the overlays comprising radio data from

four surveys (at 150, 200, 843/1400, and 1300 MHz). Our morphological

classification and host-galaxy identification relies upon the ~7-arcsec

resolution images from MeerKAT (1300 MHz). Through the visual inspection of the

overlays, 14 radio sources in the MeerKAT-2019 subset have wide-angle tail

(WAT) morphology, 10 are head-tail, and 5 have X-, S-/Z-shaped morphology. Most

of the remaining sources have the radio morphology of typical symmetric lobes.

Of 140 sources, we find host galaxies for 98 sources, leaving 42 with no

identified host galaxy. These 42 sources still have ambiguous identification

even with higher resolution images from MeerKAT.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12242v1