Papers for Tuesday, Jun 13 2023

The progenitor of SN 2023ixf was a $\sim$10$^{4.8}$ to $10^{5.0}\rm ~L_\odot$ star ($\sim$9 to $14 \rm~M_\odot$ at birth) obscured by a dusty $\dot{M} \simeq 10^{-5} \rm~M_\odot~yr^{-1}$ wind with a visual optical depth of $\tau_V \simeq 13$. This is required by the progenitor SED, the post-SN X-ray and H$\alpha$ luminosities, and the X-ray column density estimates. In Large Binocular Telescope (LBT) data spanning 5600 to 400 d before the SN, there is no evidence for optical variability at the level of $\sim$10$^3\rm~L_\odot$ in $R$ band, roughly 3 times the predicted luminosity of the obscured progenitor. This constrains direct observation of any pre-SN optical outbursts where there are LBT observations. However, models of the effects of any pre-SN outburst on the dusty wind show that an outburst of essentially any duration exceeding $\sim$5 times the luminosity of the progenitor would have detectable effects on the dust optical depth for decades. While the dust obscuration here is high, all red supergiants have dusty winds, and the destruction (or formation) of dust by even short-lived transients will always have long term effects on the observed brightness of the star because changes in the dust optical depths after a luminous transient occur very slowly.

We investigate the ultra-diffuse galaxy (UDG) UGC 9050-Dw1, which was selected because of its disturbed morphology as part of a larger sample of UDGs that display evidence for significant interactions. We use the Hubble Space Telescope's Advanced Camera for Surveys to identify globular clusters (GCs) associated with UGC 9050-Dw1, and the Jansky Very Large Array to measure its $\mathrm{H}\mathrm{I}$ content. UGC 9050-Dw1, a neighbor to the low surface brightness spiral UGC 9050, exhibits a unique UV bright central ``clump'' with clearly associated $\mathrm{H}{\mathrm{I}}$ gas and an extended stellar tidal plume to the north. We identify $52^{+4}_{-6}$ GCs, implying a specific frequency $S_\mathrm{N} = 122_{-24}^{+30}$, one of the highest reported for a UDG of this luminosity. Additionally, $\sim 20\%$ of the total light of the galaxy is contributed by GCs. Nearly uniform GC colors suggest they were formed during a single intense episode of star formation. We propose that UGC 9050-Dw1 formed via a rare dwarf merger event where induced, clumpy star formation led to its current observed properties.

The enduring tension between local and distant measurements of $H_0$ remains unresolved. It was recently pointed out that cosmic microwave background (CMB) and large-scale structure (LSS) observables are invariant under a uniform rescaling of the gravitational free-fall rates of all species present and the Thomson scattering rate between photons and electrons. We show that a unique variation of the fine-structure constant $\alpha$ and the electron mass $m_{\rm e}$ can leverage this scaling transformation to reconcile the CMB and LSS data with a broad spectrum of Hubble constant values, encompassing those inferred from local measurements. Importantly, this study demonstrates that the constraints on the variation of fundamental constants imposed by the specific recombination history are not as stringent as previously assumed. Our work highlights the critical role of the Thomson scattering rate in the existing Hubble tension and offers a distinct avenue of exploration for particle model builders.

Planets, embedded in their natal discs, harbour hot envelopes. When pebbles are accreted by these planets, the contained volatile components may sublimate, enriching the envelope and potentially changing its thermodynamical properties. However, the envelopes of embedded planets actively exchange material with the disc, which would limit the buildup of a vapour-rich atmosphere. To properly investigate these processes, we have developed a new phase change module to treat the sublimation process with hydrodynamical simultions. Combined with the recently developed multi-dust fluid approach, we conduct 2D self-consistent hydrodynamic simulations to study how pebble sublimation influences the water content of super-Earths and sub-Neptunes. We find the extent and the amount of vapour that a planet is able to hold on to is determined by the relative size of the sublimation front and the atmosphere. When the sublimation front lies far inside the atmosphere, vapour tends to be locked deep in the atmosphere and keeps accumulating through a positive feedback mechanism. On the other hand, when the sublimation front exceeds the (bound) atmosphere, the ice component of incoming pebbles can be fully recycled and the vapour content reaches a low, steady value. Low disc temperature, small planet mass and high pebble flux (omitting accretion heating by pebbles) render the planet atmosphere vapour-rich while the reverse changes render it vapour-poor. The phase change module introduced here can in future studies also be employed to model the chemical composition of the gas in the vicinity of accreting planets and around snowlines.

BL Lac objects detected at TeV energies preferentially belong to the subclass called 'high-frequency-peaked' BL Lacs (HBLs). Parsec-scale radio jets in these TeV-HBLs often show dominant, slow moving radio knots that are at most mildly superluminal. We report the first systematic campaign to characterise the Intra-Night Optical Variability (INOV) of TeV-HBLs using a representative sample of 6 such sources, all showing a fairly high degree of optical polarization. Our campaign consists of high-sensitivity monitoring of this sample in 24 sessions of more than 3 hour duration each. For these TeV-HBLs, we find a striking lack of INOV and based on this, we discuss the importance of superluminal motion of the radio knots vis-a-vis the optical polarization, as the key diagnostic for INOV detection.

The evolution of the quasar luminosity function (QLF) is fundamental to understanding the cosmic evolution of black holes (BHs) through their accretion phases. In the era of the James Webb Space Telescope (JWST), Euclid, and Nancy Grace Roman Space Telescope, their unprecedented detection sensitivity and wide survey area can unveil the low-luminosity quasar and low-mass BH population, and provide new insights into quasar host galaxies. We present a theoretical model describing BH growth from initial seeding at $z>20$ to $z\sim 4$,incorporating the duration of accretion episodes, the distribution of Eddington ratios, and the mass dependency of BH accretion rates. By constraining the model parameters with the observed QLFs at $4\leq z\leq6$ across a wide UV luminosity range ($-29<M_{\rm 1450}<-24$), we find that the high-redshift BH population grows rapidly at $z\gtrsim6$, and decelerates the pace in subsequent epochs. Toward lower redshifts ($z<6$), mass-dependent accretion inhibits the growth of high-mass BHs with $M_{\bullet}>10^8~M_\odot$, leading to mass saturation at $M_\bullet\gtrsim 10^{10}~M_\odot$. We predict the BH mass function down to $M_{\bullet}\sim 10^6~M_\odot$ for both unobscured and obscured quasar populations at $4\leq z \leq 11$, offering a benchmark for future observational tests. Our model accounts for the presence of both bright and faint quasars at $z>4$, including those discovered by JWST. Furthermore, our findings suggest two distinct pathways for the early assembly of the BH-galaxy mass correlation: the population with a BH-to-stellar mass ratio near the local value of $M_\bullet/M_{\star}\simeq5\times10^{-3}$ maintains a proximity to the relation through its evolution via moderate growth, while the population that begins to grow above the local relation accretes mass rapidly and becomes as overmassive as $M_\bullet/M_\star \sim 0.01-0.1$ by $z\sim 6$.

SU Aurigae is a widely studied T Tauri star and here we present original state-of-the-art interferometric observations with better uv and baseline coverage than previous studies. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The MIRC-X instrument at CHARA is a 6 telescope optical beam combiner offering baselines up to 331 m. We undertook image reconstruction for model-independent analysis, and fitted geometric models such as Gaussian and ring distributions. Additionally, the fitting of radiative transfer models constrains the physical parameters of the disk. Image reconstruction reveals a highly inclined disk with a slight asymmetry consistent with inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission/scattering of the far-side. Geometric models find that the underlying brightness distribution is best modelled as a Gaussian with a FWHM of $1.53\pm0.01 \mathrm{mas}$ at an inclination of $56.9\pm0.4^\circ$ and minor axis position angle of $55.9\pm0.5^\circ$. Radiative transfer modelling shows a flared disk with an inner radius at 0.16 au which implies a grain size of $0.14 \mathrm{\mu m}$ assuming astronomical silicates and a scale height of 9.0 au at 100 au. In agreement with literature, only the dusty disk wind successfully accounts for the NIR excess by introducing dust above the mid-plane. Our results confirm and provide better constraints than previous inner disk studies of SU Aurigae. We confirm the presence of a dusty disk wind in the cicumstellar environment, the strength of which is enhanced by a late infall event which also causes very strong misalignments between the inner and outer disks.

We present a robust sample of very high-redshift galaxy candidates from the first epoch of {\it JWST}/NIRCam imaging from the Next Generation Extragalactic Exploratory Deep (NGDEEP) Survey. The NGDEEP NIRCam imaging in the Hubble Ultra Deep Field Parallel Field 2 (HUDF-Par2) reaches $m=30.4$ (5$\sigma$, point-source) in F277W, making it the deepest public {\it JWST} GO imaging dataset to date. We describe our detailed data reduction process of the six-filter broad-band {\it JWST}/NIRCam imaging, incorporating custom corrections for systematic effects to produce high-quality calibrated images. Using robust photometric redshift selection criteria, we identify a sample of 38 $z \gtrsim 9$ galaxy candidates. These objects span a redshift range of $z=8.5-15.8$, and apparent magnitudes of $m_\mathrm{F277W} = 27-30.5$ AB mag, reaching $\sim 1.5$ mag deeper than previous public {\it JWST} imaging surveys. We calculate the rest-frame ultraviolet (UV) luminosity function at $z \sim$ 9 and 11, and present a new measurement of the luminosity function faint-end slope at $z \sim 11$. There is no significant evolution in the faint-end slope and number density from $z=9$ to 11. Comparing our results with theoretical predictions, we find that some models produce better agreement at the faint end than the bright end. These results will help to constrain how stellar feedback impacts star formation at these early epochs.

We present uniformly measured stellar metallicities for 463 stars in 13 Milky Way (MW) ultra-faint dwarf galaxies (UFDs; $M_V = -7.1$ to $-0.8$) using narrowband CaHK (F395N) imaging taken with the Hubble Space Telescope (HST). This represents the largest homogeneous set of stellar metallicities in UFDs, increasing the number of metallicities in these 13 galaxies by a factor of 5 and doubling the number of metallicities in all known MW UFDs. We provide the first well-populated MDFs for all galaxies in this sample, with [Fe/H] ranging from -3.0 dex to -2.0 dex, and $\sigma_{\mbox{[Fe/H]}}$ ranging from 0.3 dex to 0.7 dex. We find a nearly constant [Fe/H] ~ -2.6 over 3 decades in luminosity (~$10^2 - 10^5 L_{\odot}$), suggesting that the mass-metallicity relationship does not hold for such faint systems. We find a larger fraction (24%) of EMP ([Fe/H]<-3.0) stars across our sample compared to the literature (14%), but note that uncertainties in our most metal-poor measurements make this an upper limit. We find 19% of stars in our UFD sample to be metal-rich ([Fe/H]>-2.0), consistent with the sum of literature spectroscopic studies. MW UFDs are known to be predominantly >13 Gyr old, meaning that all stars in our sample are truly ancient, unlike metal-poor stars in the MW, which have a range of possible ages. Our UFD metallicities are not well-matched to known streams in the MW, providing further evidence that known MW substructures are not related to UFDs. We include a catalog of our stars to encourage community follow-up studies, including priority targets for ELT-era observations.

We introduce the NASA/IPAC Extragalactic Database (NED) Local Volume Sample (NED-LVS), a subset of $\sim$1.9 million objects with distances out to 1000~Mpc. We use UV and IR fluxes available in NED from all-sky surveys to derive physical properties, and estimate the completeness relative to the expected local luminosity density. The completeness relative to NIR luminosities (which traces a galaxy's stellar mass) is roughly 100% at $D<$30~Mpc and remains moderate (70%) out to 300~Mpc. For brighter galaxies ($\gtrsim L^{*}$), NED-LVS is $\sim$100% complete out to $\sim$400~Mpc. When compared to other local Universe samples (GLADE and HECATE), all three are $\sim$100% complete below 30~Mpc. At distances beyond $\sim$80~Mpc, NED-LVS is more complete than both GLADE and HECATE by $\sim$10-20%. NED-LVS is the underlying sample for the NED gravitational wave follow-up (NED-GWF) service, which provides prioritized lists of host candidates for GW events within minutes of alerts issued by the LIGO-Virgo-KAGRA collaboration. We test the prioritization of galaxies in the volume of GW170817 by 3 physical properties, where we find that both stellar mass and inverse specific star formation rate place the correct host galaxy in the top ten. In addition, NED-LVS can be used for a wide variety of other astrophysical studies: galaxy evolution, star formation, large-scale structure, galaxy environments, and more. The data in NED are updated regularly, and NED-LVS will be updated concurrently. Consequently, NED-LVS will continue to provide an increasingly complete sample of galaxies for a multitude of astrophysical research areas for years to come.

The discovery of ultra-high energy gamma-ray sources (detected at energies $\geq$ 100 TeV) thanks to highly sensitive observatories such as the High Altitude Water Cherenkov (HAWC) Observatory, the Tibet AS-gamma Experiment, and the Large High Altitude Air Shower Observatory (LHAASO), marked the beginning of the sub--PeV and PeV era in gamma-ray astrophysics (energies $\sim$ 0.1 to 100 PeV). This new astrophysics is closely related to the emerging topic of PeVatrons, in which HAWC has played a remarkable role and will continue to contribute discoveries and relevant studies thanks to the better data provided by the installed outriggers. In this paper, we present a brief overview of the PeVatrons and the HAWC observatory.

Even when the direct view toward the active nucleus is obscured, nuclear emission propagating along other directions can scatter off surrounding material, become polarized and reach the observer. Spectropolarimetry can thus be an important tool in investigating the circumnuclear geometry and kinematics of quasars on scales that cannot yet be probed via direct observations. Here we discuss an intriguing class of quasars where the polarization position angle swings by large amounts (90 deg) within an emission line. We investigate a kinematic model in which the scattering dust or electrons are in an axisymmetric outflow. We propagate Stokes parameters in a variety of geometries of emitter, scatterer and observer. We use these models to predict polarization fraction, line profiles and polarization position angles and compare them to observations. We demonstrate that the swinging polarization angle can be a result of the geometry of the outflow and the orientation of the observer. Polarization properties of a near-Eddington extremely red quasar SDSS J1652 can be successfully explained by a model in which the quasar is surrounded by a geometrically thick disk, whose `skin' is outflowing at 1000 km/s and acts as the scatterer on scales of a few tens of pc. The line of sight to the observer in this source is within or close to the skin of the torus, in agreement with multi-wavelength data. Spectropolarimetric data and models presented here strongly support the thick-disk geometry of circumnuclear material suggested by recent numerical simulations of high-rate accretion flows onto black holes.

The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z<1.1$, 0.39% over the redshift interval $1.1<z<1.9$, and 0.46% over the redshift interval $1.9<z<3.5$.

The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes good-quality spectral information from 466,447 objects targeted as part of the Milky Way Survey, 428,758 as part of the Bright Galaxy Survey, 227,318 as part of the Luminous Red Galaxy sample, 437,664 as part of the Emission Line Galaxy sample, and 76,079 as part of the Quasar sample. In addition, the release includes spectral information from 137,148 objects that expand the scope beyond the primary samples as part of a series of secondary programs. Here, we describe the spectral data, data quality, data products, Large-Scale Structure science catalogs, access to the data, and references that provide relevant background to using these spectra.

The Dark Energy Spectroscopic Instrument (DESI) is currently measuring the spectra of 40\,million galaxies and quasars, the largest such survey ever made to probe the nature of cosmological dark energy. The 4-meter Mayall telescope at Kitt Peak National Observatory has been adapted for DESI, including the construction of a 3.2-degree diameter prime focus corrector that focuses astronomical light onto a 0.8-meter diameter focal surface with excellent image quality over the DESI bandpass of 360-980nm. The wide-field corrector includes six lenses, as large as 1.1-meters in diameter and as heavy as 237\,kilograms, including two counter-rotating wedged lenses that correct for atmospheric dispersion over Zenith angles from 0 to 60 degrees. The lenses, cells, and barrel assembly all meet precise alignment tolerances on the order of tens of microns. The barrel alignment is maintained throughout a range of observing angles and temperature excursions in the Mayall dome by use of a hexapod, which is itself supported by a new cage, ring, and truss structure. In this paper we describe the design, fabrication, and performance of the new corrector and associated structure, focusing on how they meet DESI requirements. In particular we describe the prescription and specifications of the lenses, design choices and error budgeting of the barrel assembly, stray light mitigations, and integration and test at the Mayall telescope. We conclude with some validation highlights that demonstrate the successful corrector on-sky performance, and list some lessons learned during the multi-year fabrication phase.

We present the one-dimensional Lyman-$\alpha$ forest power spectrum measurement using the first data provided by the Dark Energy Spectroscopic Instrument (DESI). The data sample comprises $26,330$ quasar spectra, at redshift $z > 2.1$, contained in the DESI Early Data Release and the first two months of the main survey. We employ a Fast Fourier Transform (FFT) estimator and compare the resulting power spectrum to an alternative likelihood-based method in a companion paper. We investigate methodological and instrumental contaminants associated to the new DESI instrument, applying techniques similar to previous Sloan Digital Sky Survey (SDSS) measurements. We use synthetic data based on log-normal approximation to validate and correct our measurement. We compare our resulting power spectrum with previous SDSS and high-resolution measurements. With relatively small number statistics, we successfully perform the FFT measurement, which is already competitive in terms of the scale range. At the end of the DESI survey, we expect a five times larger Lyman-$\alpha$ forest sample than SDSS, providing an unprecedented precise one-dimensional power spectrum measurement.

We present and validate the catalog of Lyman-{\alpha} forest fluctuations for 3D analyses using the Early Data Release (EDR) from the Dark Energy Spectroscopic Instrument (DESI) survey. We used 96,317 quasars collected from DESI Survey Validation (SV) data and the first two months of the main survey (M2). We present several improvements to the method used to extract the Lyman-{\alpha} absorption fluctuations performed in previous analyses from the Sloan Digital Sky Survey (SDSS). In particular, we modify the weighting scheme and show that it can improve the precision of the correlation function measurement by more than 20%. This catalog can be downloaded from https://data.desi.lbl.gov/public/edr/vac/edr/lya/fuji/v0.3, and it will be used in the near future for the first DESI measurements of the 3D correlations in the Lyman-{\alpha} forest.

We perform a SubHalo Abundance Matching (SHAM) study with two algorithms: \{$\sigma, V_{\rm ceil}, v_{\rm smear}$\}-SHAM and \{$\sigma, V_{\rm ceil},f_{\rm sat}$\}-SHAM. These are designed to reproduce the clustering on 5--30$\,h^{-1}\,{\rm Mpc}$ of DESI One Percent Survey targets: Luminous Red Galaxies (LRGs), Emission Line Galaxies (ELGs) and Quasi-Stellar Objects (QSOs) at $0.4<z<3.5$. These SHAM methods are implemented on (sub)haloes from the dark-matter-only UNIT simulations. $V_{\rm ceil}$ is the incompleteness of the massive host (sub)haloes and is the key to the generalized SHAM. $v_{\rm smear}$ models the uncertainties in the spectroscopic redshift measurement. We cross-check it with results from redshift differences $\Delta v$ of repeat observations, fitted well by Lorentzian profiles in general. A free satellite fraction $f_{\rm sat}$ is necessary to reproduce the clustering of ELGs. We find ELGs present a more complex galaxy--halo mass relation than LRGs due to their weak constraints on $\sigma$. Each tracer shows a certain range of $V_{\rm ceil}$ values, corresponding to the large-stellar-mass incompleteness of LRGs, the quenched star formation of ELGs and the quenched black hole accretion of QSOs. The impact of the statistical redshift uncertainty on ELG clustering is negligible. For LRGs, a Gaussian profile is a better choice for sub-samples at redshift bins. The best-fitting satellite fraction for DESI ELGs is around 4 per cent, which is lower than previous observations. The mean halo mass log$_{10}(\langle M_{\rm vir}\rangle)$ in $h^{-1}\,M_\odot$ for LRGs, ELGs and QSOs are ${13.16\pm0.01}$, ${11.90\pm0.06}$ and ${12.66\pm0.45}$ respectively. We find the mean peak maximum circular velocity $\langle V_{\rm peak}\rangle$ for ELGs is twice smaller than that of LRGs and QSOs. Our SHAM algorithms can also provide galaxy mocks for cosmological tests.

We present the first comprehensive Halo Occupation Distribution (HOD) analysis of the DESI One-Percent survey Luminous Red Galaxy (LRG) and Quasi-Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 Mpc/h in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of fsat = 11+-1%, a mean halo mass of log10 Mh = 13.40+0.02-0.02, and a linear bias of blin = 1.93+0.06-0.04. For LRGs in 0.6 < z < 0.8, we find fsat = 14+-1%, log10 Mh = 13.24+0.02-0.02, and blin = 2.08+0.03-0.03. For QSOs, we infer fsat = 3+8-2%, log10 Mh = 12.65+0.09-0.04, and blin = 2.63+0.37-0.26 in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high-fidelity galaxy mocks. We also study the redshift-evolution of the DESI LRG sample from z = 0.4 up to z = 1.1, revealing significant and interesting trends in mean halo mass, linear bias, and satellite fraction.

The one-dimensional power spectrum $P_{\mathrm{1D}}$ of the Ly$\alpha$ forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of $P_{\mathrm{1D}}$ with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of $54~600$ quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with two-dimensional image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly$\alpha$ emission line. In a companion paper, we present a similar analysis based on the Fast Fourier Transform estimate of the power spectrum. We conclude with a comparison of these two approaches and implications for the upcoming DESI Year 1 analysis.

In current Dark Energy Spectroscopic Instrument (DESI) survey, emission line galaxies (ELGs) and luminous red galaxies (LRGs) are essential for mapping dark matter distribution at $z \sim 1$. We measure the auto and cross correlation functions of ELGs and LRGs at $0.8<z\leq 1.0$ from the DESI One-Percent survey. Following Gao et al. (2022), we construct the galaxy-halo connections for ELGs and LRGs simultaneously. With the stellar-halo mass relation (SHMR) for the whole galaxy population (i.e. normal galaxies), LRGs can be directly selected according to the stellar mass, while ELGs can also be selected randomly based on the observed number density of each stellar mass once the probability $P_{\mathrm{sat}}$ of a satellite galaxy becoming an ELG is determined. We demonstrate that the observed small scale clustering prefers a halo mass-dependent $P_{\mathrm{sat}}$ model rather than a constant. With this model, we can well reproduce the auto correlations of LRGs and the cross correlations between LRGs and ELGs at $r_{\mathrm{p}}>0.1$ $\mathrm{Mpc}\,h^{-1}$. We can also reproduce the auto correlations of ELGs at $r_{\mathrm{p}}>0.3$ $\mathrm{Mpc}\,h^{-1}$ ($s>1$ $\mathrm{Mpc}\,h^{-1}$) in real (redshift) space. Although our model has only seven parameters, we show that it can be extended to higher redshifts and reproduce the observed auto correlations of ELGs in the whole range of $0.8<z<1.6$, which enables us to generate a lightcone ELG mock for DESI. With the above model, we further derive halo occupation distributions (HODs) for ELGs which can be used to produce ELG mocks in coarse simulations without resolving subhalos.

The observations from the Dark Energy Spectroscopic Instrument (DESI) will significantly increase the numbers of known extremely metal-poor stars by a factor of ~ 10, improving the sample statistics to study the early chemical evolution of the Milky Way and the nature of the first stars. In this paper we report high signal-to-noise follow-up observations of 9 metal-poor stars identified during the DESI commissioning with the Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy (OSIRIS) instrument on the 10.4m Gran Telescopio Canarias (GTC). The analysis of the data using a well-vetted methodology confirms the quality of the DESI spectra and the performance of the pipelines developed for the data reduction and analysis of DESI data.

We present the first period-luminosity (PL) and period-Wesenheit (PW) relations in the Sloan-Pans-STARRS gP1rP1iP1 bands for classical fundamental mode Cepheids in the Milky Way. We used a relatively modest number of 76 stars for the PL and 84-85 stars for the PW relations calibration. The data for the project were collected with the network of 40-cm telescopes of Las Cumbres Observatory, and Gaia Data Release 3 parallaxes were used for the calculations. These gri-band PL and PW relations calibrations will be a useful tool for distance determinations in the era of large sky surveys using the Sloan photometric system, especially with the near-future start of the Large Synoptic Survey of Space and Time (LSST).

Stellar ages are critical for understanding the temporal evolution of a galaxy. We calculate the ages of over 6000 red giant branch stars in the Large Magellanic Cloud (LMC) observed with SDSS-IV / APOGEE-S. Ages are derived using multi-band photometry, spectroscopic parameters (T$_\text{eff}$, $\log{g}$, [Fe/H], and [$\alpha$/Fe]) and stellar isochrones and the assumption that the stars lie in a thin inclined plane to get accurate distances. The isochrone age and extinction are varied until a best match is found for the observed photometry. We perform validation using the APOKASC sample, which has asteroseismic masses and accurate ages, and find that our uncertainties are $\sim$20% and range from $\sim$1$-$3 Gyr for the calculated age values. Here we present the LMC age map as well as the age-radius relation and an accurate age-metallicity relation (AMR). The age map and age-radius relation reveal that recent star formation in the galaxy was more centrally located and that there is a slight dichotomy between the north and south with the northern fields being slightly younger. The northern fields that cover a known spiral arm have median ages of $\gtrsim$ 2 Gyr, which is the time when an interaction with the SMC is suggested to have happened. The AMR is mostly flat especially for older ages although recently (about 2.0-2.5 Gyr ago) there is an increase in the median [Fe/H]. Based on the time frame, this might also be attributed to the close interaction between the LMC and SMC.

N$_2$ plays a profound role in supporting processes on the surface and in the atmosphere of Pluto, yet the origin of Pluto's N$_2$ remains a mystery. However, this may begin to change as the $^{14}$N/$^{15}$N ratio of N$_2$ was recently estimated based on a non-detection of HC15N in Pluto's atmosphere, while accounting for $^{14}$N/$^{15}$N fractionation between HCN and N$_2$. Here, I show that, if this latter step of translating isotope ratios is adequately understood, then the derived $^{14}$N/$^{15}$N ratio represents the first distinguishing constraint on the origin of N$_2$. One finding of the present study is that isotopic fractionation between atmospheric N$_2$ and N$_2$-rich ices on Pluto does not appear to be significant. I infer a lower limit of ~197 for the $^{14}$N/$^{15}$N ratio of the dominant reservoir of N$_2$ on Pluto; i.e., mostly contained in Sputnik Planitia. From this lower limit, an endmember ammonia source of Pluto's N$_2$ can be ruled out. I perform N isotope mixing calculations that enable quantitative understanding of the relationships between contributions by primordial N$_2$, NH$_3$, and nitrogen originally sourced in organic materials (Norg) to Pluto's observed N$_2$ inventory. These calculations address how uncertainties in the isotopic composition of Norg and the history of atmospheric escape affect the allowed ranges of primordial N$_2$, NH$_3$, and Norg contributions. While present uncertainties are substantial, I find that a contribution by primordial N$_2$, Norg, or both is implied, and the sum of their contributions should be at least ~45%. Hence, Pluto likely formed from building blocks that were cold enough to trap N$_2$, or Pluto has a thermally processed and dynamic interior that supports generation of N$_2$ from Norg and N$_2$ transport to the surface. The lower limit on $^{14}$N/$^{15}$N suggests that NH$_3$ has been a less significant contributor to the origin of N$_2$ on Pluto than on Titan.

Some gamma-ray bursts (GRBs) have an afterglow in the tera-electronvolt (TeV) band, but the early onset of this afterglow has not been observed. We report observations with the Large High Altitude Air Shower Observatory of the bright GRB 221009A, which serendipitously occurred within the instrument field of view. More than 64,000 photons (above 0.2~TeV) were detected within the first 3000 seconds. The TeV photon flux began several minutes after the GRB trigger, then rose to a flux peak about 10 seconds later. This was followed by a decay phase, which became more rapid at $\sim 650\,{\rm s}$ after the peak. The emission can be explained with a relativistic jet model with half-opening angle $\sim 0.8^\circ$, consistent with the core of a structured jet. This interpretation could explain the high isotropic energy of this GRB.

As a considerable investment of time from various telescope facilities were dedicated toward studying the Spiderweb protocluster at $z=2.2$, it so far remains one of the most extensively studied protocluster. We report here the latest results in this field, adding a new dimension to previous research on cluster formation at high redshift. Previous studies have reported a significant overdensity ($\delta\sim10$) of massive H$\alpha$ (+ [Nii]) -emitting galaxies in 3700 comoving Mpc$^3$. Many of these were previously considered to be dusty, active star-forming galaxies, given their rest-frame optical and infrared features. However, this study argues that a third of them are more likely to be "passively-evolving" galaxies with active galactic nuclei (AGNs) rather than star-forming galaxies, given the multi-wavelength spectral energy distribution (SED) fitting including an AGN component. For their SED-based star formation rates to be valid, bulk of their H$\alpha$ + [Nii] emission should come from the central AGNs. This difference in interpretation between this work and past studies, including ours, is particularly supported by the recent deep Chandra X-ray observation. Furthermore, we have spectroscopically confirmed a quiescent nature for one of these AGNs, with its multiple stellar absorption lines but also low ionisation emission lines. This important update provides new insights into the role of AGNs in forming the cluster red sequence observed in the present-day universe.

The Imaging Atmospheric Cherenkov technique allows to detect very high energy gamma rays from few tens of GeV to hundreds of TeV using ground-based instrumentation. At these energies a gamma ray generates a shower of secondary particles when it enters the Earth's atmosphere. These particles emit Cherenkov light in the visible and near UV ranges. The Cherenkov light produced by the shower reaches the ground as a short pulse of a few nanosecond duration over a large circle of around 100 m radius (a light pool). This pulse of light can be imaged with telescopes provided with fast photodetectors and electronics. Combining the images of several telescopes distributed over this light pool allows to estimate the gamma-ray energy and incident direction, and to reject gamma rays from the strong background of charged cosmic rays. The collection area of an array of a few telescopes is of the order of the area of the light pool, i.e. $>$10$^5$m$^2$. Such an array reaches a sensitivity of a few millicrabs at 100 GeV energies in 50 hours of observations, an angular resolution of $\sim$5 arcmin and a spectral resolution of $\sim$10%. This chapter describes the technical implementation of Imaging Atmospheric Cherenkov telescopes and describes how the data are analyzed to reconstruct the physical parameters of the primary gamma rays.

We examine the reasons for discrepancies between two alternative approaches to modeling small-amplitude tides in binary systems. The 'direct solution' (DS) approach solves the governing differential equations and boundary conditions directly, while the 'modal decomposition' (MD) approach relies on a normal-mode expansion. Applied to a model for the primary star in the heartbeat system KOI-54, the two approaches predict quite different behavior of the secular tidal torque. The MD approach exhibits the pseudosynchronization phenomenon, where the torque due to the equilibrium tide changes sign at a single, well-defined and theoretically predicted stellar rotation rate. The DS approach instead shows 'blurred' pseudosynchronization, where positive and negative torques intermingle over a range of rotation rates. We trace a major source of these differences to an incorrect damping coefficient in the profile functions describing the frequency dependence of the MD expansion coefficients. With this error corrected some differences between the approaches remain; however, both are in agreement that pseudosynchronization is blurred in the KOI-54 system. Our findings generalize to any type of star for which the tidal damping depends explicitly or implicitly on the forcing frequency.

Very metal-poor stars that have $[\text{Fe}/\text{H}]<-2$ and that are enhanced in C relative to Fe ($[\text{C}/\text{Fe}]>0.7$) but have no enhancement of heavy elements ($[\text{Ba}/\text{Fe}]<0$) are known as carbon-enhanced metal-poor normal (CEMP-no) stars. These stars are thought to be produced from a gas that was polluted by the supernova (SN) ejecta of the very first generation (Pop III) massive stars. Although theoretical models of SN explosions from massive Pop III stars can explain the relative abundance pattern, the very high enrichment of C ($A(\text{C})\gtrsim 6$) observed in many of the CEMP-no stars is difficult to explain when a reasonable dilution of the supernova ejecta, that is consistent with detailed simulation of metal mixing in minihaloes, is adopted. We explore rapidly rotating Pop III stars that undergo efficient mixing and reach a quasi-chemically homogeneous (QCH) state. We find that rapidly rotating models that reach the QCH state can eject large amounts of C in the wind and that the resulting dilution of the wind ejecta in the interstellar medium can lead to a C enrichment of $A(\text{C})\lesssim7.75$. This can naturally explain the high C enrichment observed in CEMP-no stars. The core of QCH stars can produce up to an order of magnitude of more C than non-rotating progenitors of similar mass and the resulting SN can lead to a C enrichment of $A(\text{C})\lesssim7$. We find that the abundance pattern from our models that use dilution masses that are consistent with simulations can provide an excellent match to observed abundance patterns in most of the known CEMP-no stars. Our rapidly rotating massive Pop III stars are a promising site for explaining the high C enhancement in the early Galaxy that is the birthplace of the CEMP-no stars. Our work indicates that a substantial fraction of Pop III stars were likely rapid rotators.

We report the results of optical \'echelle spectroscopy with the Southern African Large Telescope (SALT) of the mass donor star BSDL\,923 in the neutron star (NS) high-mass X-ray binary XMMU\,J051342.6$-$672412 associated with the LMC supernova remnant (SNR) MCSNR\,J0513$-$6724. We found that BSDL\,923 is a B0.7\,III star with double peaked emission lines originating in a circumbinary disk-like structure. This classification and the presence of double-peaked emission lines imply that BSDL\,923 is a Be star. Modelling with the stellar atmosphere code {\sc fastwind} was used to derive the effective temperature $T_{\rm eff}=27\pm1$\,kK, surface gravity $\log g=3.22\pm0.10$, projected rotational velocity $v\sin i\approx100\pm45 \, \kms$, colour excess $E(B-V)=0.53\pm0.05$\,mag, and luminosity $\log(L_*/\lsun)=5.46\pm0.10$ of BSDL\,923, as well as to show that the surface of this star is polluted with $\alpha$-elements (O, Mg and Si) from the supernova ejecta. We found also that the NS is orbiting BSDL\,923 in an eccentric ($e=0.158\pm0.061$) orbit with the orbital period of $1.280\pm0.006$\,d and the semi-major axis of $17\pm3 \, \rsun$, and the radius of BSDL\,923 is $25\pm5 \,\rsun$. We speculate that the NS is embedded in the atmosphere of BSDL\,923 either because it was kicked at birth towards this star or because of inflation of BSDL\,923 caused by the energy input from the supernova blast wave. Using long-slit spectroscopy with SALT, we searched for possible signs of the SNR shell in the 2D spectrum, but did not find them. This lack of detection is consistent with the young age ($\approx4^{+2} _{-1}$\,kyr) of MCSNR\,J0513$-$6724, implying that it is still in the adiabatic (non-radiative) phase.

The high energy emission of rotation powered pulsars is supposed to be produced in "gaps" in the pulsar magnetosphere where charges are accelerated and currents are produced. The rest of the magnetosphere is supposed to be mostly a "force-free" plasma without any currents. Two important currents are the main current that flows away from the pulsar, that produces the observed radiation, and the current that returns to the pulsar to maintain charge neutrality. This work attempts to study the return current in the Crab pulsar using the soft X-ray data from the {\it{NICER}} observatory. It is assumed that the two currents vary as a function of time. This would modulate the electric fields in the "gaps", which would affect the observed X-ray flux. These flux variations will show up only in the on-pulse phases, while those caused by the Crab Nebula, instrumental effects, etc. will be present in the off-pulse phases also. This work obtains the correlation coefficient of the flux variations in the two peaks of the Crab pulsar, after removing the off-pulse flux variations. No correlation was observed; its error of $0.000012$ sets an upper limit of $0.036\%$ on the rms variation of correlated X-ray flux in the Crab pulsar. Reasons exist for the return current variations to be correlated, while the main current variations are probably uncorrelated. So the above number is considered an upper limit on correlated return current variations, which may be an important constraint for pulsar magnetospheric structure.

Although low surface brightness galaxies (LSBs) contribute a large fraction to the number density of galaxies, their properties are still poorly known. LSBs are often considered dust poor, based only on a few studies. We use, for the first time, a large sample of LSBs and high surface brightness galaxies (HSBs) with deep observational data to study their dust properties as a function of surface brightness. Our sample consists of 1631 optically selected galaxies at $z < 0.1$ from the North Ecliptic Pole (NEP) wide field. We use the large set of data available in this field, from UV to FIR. We measured the optical size and the surface brightness of the targets, and analyzed their spectral energy distribution using the CIGALE fitting code. We found that the specific star formation rate and specific infrared luminosity (total infrared luminosity per stellar mass) remain mostly flat as a function of surface brightness for both LSBs and HSBs that are star-forming but decline steeply for the quiescent galaxies. The majority of LSBs in our sample have negligible dust attenuation (A$_{V} < 0.1$ mag), except for about 4% of them that show significant attenuation with a mean A$_{V}$ of 0.8 mag. We found that these LSBs also have a high $\textit{r}$-band mass-to-light ratio ($M/L_r>3$ M$_{\odot}$/L$_{\odot}$), and show similarity to the extreme giant LSBs from the literature, indicating a possibly higher dust attenuation in giant LSBs as well. This work provides a large catalog of LSBs and HSBs with detailed measurements of their several optical and infrared physical properties. Our results suggest that the dust content of LSBs is more varied than previously thought, with some of them having significant attenuation making them fainter than their intrinsic value. This will have serious implications for the observation and analysis of LSBs with current/upcoming surveys like JWST and LSST.

With the launch of JWST and other scheduled missions aimed at probing the distant Universe, we are entering a new promising era for high-$z$ astronomy. One of our main goals is the detection of the first population of stars (Population III or Pop III stars), and models suggest that Pop III star formation is allowed well into the Epoch of Reionization (EoR), rendering this an attainable achievement. In this paper, we focus on our chance of detecting massive Pop IIIs at the moment of their death as Pair-Instability Supernovae (PISNe). We estimate the probability of discovering PISNe during the EoR in galaxies with different stellar masses ($7.5 \leq \mathrm{Log}(M_\star/\mathrm{M_\odot}) \leq 10.5$) from six dustyGadget simulations of $50h^{-1}$ cMpc per side. We further assess the expected number of PISNe in surveys with JWST/NIRCam and Roman/WFI. On average, less than one PISN is expected in all examined JWST fields, while $\simeq 1.5 \, \eta_\mathrm{III}$ PISNe may be found in a $\sim 1$ deg$^2$ Roman field, with potential for increased discoveries when considering more optimistic choices for the Pop III initial mass function and formation efficiency $\eta_\mathrm{III}$, and when including the contribution of coarsely-resolved environments. JWST/NIRCam and Roman/WFI allow the detection of massive-progenitor ($\sim 250$ $\mathrm{M_\odot}$) PISNe throughout all the optimal F200W-F356W, F277W-F444W, and F158-F213 colors. PISNe are also predominantly located at the outskirts of their hosting haloes, facilitating the disentangling of underlying stellar emission thanks to the spatial-resolution capabilities of the instruments.

We have completed our observational program to search for wide binary systems with non-coeval components in the southern sky and report our results here. The final set of four systems was spectroscopically investigated in this paper. No binary systems with components of different ages were found among them. Taking into account our previous studies, we estimate the fraction of such binaries (i.e., binaries formed, presumably, by capture) not higher than 0.06~\%. The study will be continued on the northern sky.

The emergence of active regions on the Sun is an integral feature of the solar dynamo mechanism. However, details about the generation of active-region-scale magnetism and the journey of this magnetic flux to the photosphere are still in question. Shifting paradigms are now developing for the source depth of the Sun's large-scale magnetism, the organization of this magnetism into fibril flux tubes, and the role of convection in shaping active-region observables. Here we review the landscape of flux emergence theories and simulations, highlight the role flux emergence plays in the global dynamo process, and make connections between flux emergence on the Sun and other cool stars. As longer-term and higher fidelity observations of both solar active regions and their associated flows are amassed, it is now possible to place new constraints on models of emerging flux. We discuss the outcomes of statistical studies which provide observational evidence that flux emergence may be a more passive process (at least in the upper convection zone); dominated to a greater extent by the influence of convection and to a lesser extent by buoyancy and the Coriolis force acting on rising magnetic flux tubes than previously thought. We also discuss how the relationship between stellar rotation, fractional convection zone depth, and magnetic activity on other stars can help us better understand flux emergence processes. Looking forward, we identify open questions regarding magnetic flux emergence that we anticipate can be addressed in the next decade with further observations and simulations.

Most stars form in multiple star systems. For a better understanding of their formation processes, it is important to resolve the individual protostellar components and the surrounding envelope and disk material at the earliest possible formation epoch because the formation history can be lost in a few orbital timescales. Here we present the ALMA observational results of a young multiple protostellar system, IRAS 04239+2436, where three well-developed large spiral arms were detected in the shocked SO emission. Along the most conspicuous arm, the accretion streamer was also detected in the SO$_2$ emission. The observational results are complemented by numerical magneto-hydrodynamic simulations, where those large arms only appear in magnetically weakened clouds. The numerical simulations also suggest that the large triple spiral arms are the result of gravitational interactions between compact triple protostars and the turbulent infalling envelope.

The galactic component in clusters is commonly thought to be generally nonrotating and in a dynamical state different from that of a collisionally relaxed system. In practice, a test of such a picture is often not available. We consider the member galaxies of two clusters, Abell S1063 and MACS J1206.2$-$0847, and study the possible presence of mean rotation and some properties of their distribution in phase space. We look for empirical evidence of factors normally found in collisionally relaxed systems and others characteristic of violently-relaxed collisionless systems. Starting from the CLASH-VLT data, we obtain positions, stellar masses, and individual line-of-sight velocities for a large number of galaxies (N_{AS1063}=1200 and N_{M1206}=650) extending out to 1.6 (Abell) and 2.5 (MACS) times the radius r_{200}. We study the spatial distribution of the galaxy velocities and the properties of the available galaxy sets when divided in stellar mass bins. To test the presence of velocity dispersion anisotropy we compare the results based on the Jeans equations with those obtained by assuming a specific form of the galaxy distribution function incorporating the picture of violent relaxation, where the total gravitational potential is imposed as set by the available gravitational lensing observations. We find evidence of systematic rotation in both clusters, with significant rotation in each core (within 0.5' from the center) and no signatures of rotation at large radii. While no signs are found of energy equipartition, there is a clear indication of (stellar) mass segregation. Velocity dispersion anisotropy is present and qualitatively similar to that found in violently relaxed collisionless systems; this last conclusion is strengthened by the overall success in matching the observations with the predictions of the physically justified distribution function.

It is known that the weak state of the heliosphere due to diminished solar activity in cycle 24 back-reacted on coronal mass ejections (CMEs) to make them appear wider for a given speed. One of the consequences of the weak state of the heliosphere is that more CMEs appear as halo CMEs (HCMEs), and halos are formed at shorter heliocentric distances. Current predictions for the strength of solar cycle (SC) 25 range from half to twice the strength of SC 24. We compare the HCME occurrence rate and other properties during the rise phase of cycles 23, 24, and 25 to weigh in on the strength of SC 25. We find that HCME and solar wind properties in SC 25 are intermediate between SCs 23 and 24, but closer to SC 24. The HCME occurrence rate, normalized to the sunspot number, is higher in SCs 24 and 25 than in SC 23. The solar wind total pressure in SC 25 is ~35% smaller than that in SC 23. Furthermore, the occurrence rates of high-energy solar energetic particle events and intense geomagnetic storms are well below the corresponding values in SC 23, but similar to those in SC 24. We conclude that cycle 25 is likely to be similar to or slightly stronger than cycle 24, in agreement with polar-field precursor methods for cycle 25 prediction

Starlink Generation 2 Mini satellites are fainter than Gen 1 spacecraft despite their larger size. The mean of apparent magnitudes for satellites in brightness mitigation mode is 7.06 +/- 0.10. When these magnitudes are adjusted to a uniform distance of 1,000 km that mean is 7.87 +/- 0.09. The brightness mitigation mode reduces distance-adjusted satellite luminosity by a factor of 12 relative to spacecraft that are not mitigated.

Context. It is almost banal to say that the interstellar medium (ISM) is structurally and thermodynamically complex. But the variety of the governing processes, including stellar feedback, renders the investigation challenging. High latitude molecular clouds (HLMCs) with no evidence of internal star formation, such as MBM 40, are excellent sites for studying the chemistry and dynamic evolution of the cold neutral ISM. Aims. We used this high latitude cloud as an exemplar for the dynamical and chemical processes in the diffuse interstellar medium. Methods. We analyzed new and archival $^{12}$CO, $^{13}$CO, CH, HCO$^+$, CS, H$_2$CO, HCN data from Five College Radio Observatory (FCRAO), Onsala Space Observatory (OSO), Arizona Radio Observatory (ARO) and W. Gordon telescope (Arecibo) combined with the Galactic Arecibo L-band Feed Array HI (GALFA-HI) HI 21 cm data set, to study the chemistry, thermal state, and dynamics of MBM 40. A new dynamical analytical approach was adopted by considering each line profile as a line of sight Probability Distribution Function (PDF) of the turbulence weighted by gas emissivity. Results. The atomic and molecular gas are smoothly distributed in space and velocity. No steep transition is seen between circumcloud atomic and cloud molecular gas in either radial velocity or structure. We proposed a topology of the cloud from the molecular tracers, a contorted filamentary structure that is shaped by a broad embedding shear flow in the neutral atomic gas. Comparative examination of different molecular tracers shows that $^{13}$CO, H$_2$CO and CS arise from only denser molecular cores, where $^{12}$CO, CH and HCO$^+$ traces diffuse gas with broader range of dynamics.

Modern radio telescopes strongly rely on accurate computational electromagnetic tools for "beam" models. Especially for densely-packed aperture array radio telescopes, the only feasible way to produce accurate models of the individual embedded element patterns is by using electromagnetic codes. In this paper, the accuracy of two models computed by different commercial codes is evaluated for one station of the SKA-Low radio telescope. Except for a couple of critical frequencies, the amplitude and phase errors are low enough to allow a beamformer efficiency higher than 99%.

Atmospheric turbidity is one of the key factors influencing the propagation of artificial light into the environment during cloudless nights. High aerosol loading can reduce the visibility of astronomical objects, and thus information on atmospheric pollution is critical for the prediction of the night sky brightness (NSB) distribution. In particular, the aerosol optical depth (AOD) and asymmetry parameter (g) are among the most important aerosol properties influencing the NSB amplitudes. However, these two parameters are rarely available at astronomical sites. Here, we develop a method for AOD and g retrievals from clear-sky radiometry carried out around sunset or sunrise, shortly before or after night-time observation is intended. The method allows for reducing the number of unknowns needed in the processing and interpretation of night sky radiances, and thus provides an efficient tool for gathering input data to present skyglow simulators. The practice of collecting information about aerosols in this way could become a routine part of astronomical observations, much like observing standard stars to obtain extinction coefficients. If the procedure were conducted around sunset and the data were quickly reduced, it could offer an on-the-spot estimate of the NSB for the night ahead. The error analysis is performed using the theoretical model, while taking into account experimental errors of radiance readings. The capability of the method is demonstrated in a field experiment conducted under cloudless conditions.

As one of the most detectable types of terrestrial planets, lava worlds are highly prioritized targets for exoplanet atmosphere characterization since their atmospheres may reveal what they are made of and how. Our work examines the possibility of true polar wander (TPW) occurring on these ultra-hot tidally-locked planets, powered by mass redistribution from atmospheric flow between the hot permanent day-side and the cold permanent night-side. We find that lava planets within a certain mass and temperature range may undergo TPW, and this likelihood increases with star mass. As a result of TPW, the magma ocean and atmospheric compositions may be less evolved (refractory-enriched) than previously thought and may be validated by exoplanet demographic surveys.

Using the Multi Unit Spectroscopic Explorer (MUSE) instrument on the Very Large Telescope, we identified three bubble nebulae (denoted as A, B, and C) around an ultraluminous X-ray source (ULX) in NGC 55. Bubble A shows a regular elliptical shape surrounding the ULX, with a morphology similar to the canonical ULX bubble around NGC 1313 X-2. It is most likely inflated by the ULX disk wind with a mechanical power close to $10^{39}$ erg s$^{-1}$. Bubble B lies 11\arcsec\ away from the ULX on the sky plane and is not contiguous to Bubble A. It displays a bow shock like morphology, and is likely driven by a collimated dark jet from the ULX with a mechanical power of about $3 \times 10^{38}$ erg s$^{-1}$. If this scenario is correct, we predict that Bubble B should present radio emission with a flux of about $1 - 10^2$ $\mu$Jy at 5 GHz. Bubble C appears within Bubble A, with a velocity and velocity dispersion distinct from the rest of Bubble A. Its nature is unclear and could be part of Bubble A as a result of low local density. The optical counterpart of ULX-1 exhibits a broad H$\alpha$, consistent with emission from a hot disk wind.

The International VLBI Service for Geodesy & Astrometry (IVS) regularly provides high-quality data to produce Earth Orientation Parameters (EOP), and for the maintenance and realization of the International Terrestrial and Celestial Reference Frames, ITRF and ICRF. The first iteration of the celestial reference frame (CRF) at radio wavelengths, the ICRF1, was adopted by the International Astronomical Union (IAU) in 1997 to replace the FK5 optical frame. Soon after, the IVS began official operations and in 2009 there was a significant increase in data sufficient to warrant a second iteration of the CRF, ICRF2. The most recent ICRF3, was adopted by the IAU in 2018. However, due to the geographic distribution of observing stations being concentrated in the Northern hemisphere, CRFs are generally weaker in the South due to there being fewer Southern Hemisphere observations. To increase the Southern Hemisphere observations, and the density, precision of the sources, a series of deep South observing sessions was initiated in 1995. This initiative in 2004 became the IVS Celestial Reference Frame Deep South (IVS-CRDS) observing program. This paper covers the evolution of the CRDS observing program for the period 1995 to 2021, details the data products and results, and concludes with a summary of upcoming improvements to this ongoing project.

We have devised a predominantly Naive Bayes method to classify the optical/IR matches to X-ray sources detected by Chandra in the Cygnus OB2 association into foreground, member, and background objects. We employ a variety of X-ray, optical, and infrared characteristics to construct likelihoods using training sets defined by well-measured sources. Combinations of optical photometry from SDSS (riz) and IPHAS (riHa), IR magnitudes from UKIDSS and 2MASS (JHK), X-ray quantiles and hardness ratios, and estimates of extinction Av are used to compute the relative probabilities that a given source belongs to one of the classes. We use Principal Component Analysis of photometric magnitude combinations to isolate the best axes for classification. We incorporate measurement errors into the classification. We evaluate the accuracy of the classification by inspection and reclassify a number of sources based on IR magnitudes, presence of disks, and X-ray spectral hardness. We also consider systematic errors due to extinction. We find that about 6100 objects are association members, 1400 are background, and 500 are foreground objects. The overall classification accuracy is 95%.

In this paper we present the first set of 3D magnetohydrodynamic (MHD) simulations performed with the riemann geomesh code. We study the dynamics of the magnetically channeled winds of magnetic massive stars in full three dimensions using a code that is uniquely suited to spherical problems. Specifically, we perform isothermal simulations of a smooth wind on a rotating star with a tilted, initially dipolar field. We compare the mass-loss, angular momentum loss, and magnetospheric dynamics of a template star (with the properties that are reminiscent of the O4 supergiant {\zeta} Pup) over a range of rotation rates, magnetic field strengths, and magnetic tilt angles. The simulations are run up to a quasi-steady state and the results are observed to be consistent with the existing literature, showing the episodic centrifugal breakout events of the mass outflow, confined by the magnetic field loops that form the closed magnetosphere of the star. The catalogued results provide perspective on how angular-momentum loss varies for different configurations of rotation rate, magnetic field strength, and large magnetic tilt angles. In agreement with previous 2D MHD studies, we find that high magnetic confinement reduces the overall mass-loss rate, and higher rotation increases the mass-loss rate. This and future studies will be used to estimate the angular-momentum evolution, spin-down time, and mass-loss evolution of magnetic massive stars as a function of magnetic field strength, rotation rate, and dipole tilt.

The kinematic behavior of superluminal components observed at 43 GHz in blazar 3C454.3 were model-fitted with their light curves interpreted in terms of their Doppler-boosting effect. The relation between the flux evolution of superluminal components and their accelerated/decelerated motion or the increase/decrease in their Lorentz/Doppler factor was investigated. The precessing jet-nozzle scenario previously proposed by Qian et al. (1991, 2018a, 2021) and Qian (2018b, 2022a, 2022b) was applied to consistently model-fit the kinematic behavior and light curves for two superluminal components (B4 and B6) measured by Jorstad et al. (2005). For both B4 and B6 which were ascribed respectively to the jet-A and jet-B of the double-jet structure assumed for 3C454.3, their kinematic features were well model-fitted with their bulk Lorentz factor and Doppler factor (as function of time) convincingly derived. It is shown that the light curves of the radio bursts associated with knot B4 and knot B6 can be well explained in terms of their Doppler boosting effect. Similarly, for the knot R3 observed at 15GHz (Qian et al. 2014, Britzen et al. 2013) the interpretation of its kinematic behavior and light curve is presented in the Appendix. We emphisize that the interpretation of the flux evolution of superluminal components combined with the model-fit of their kinematics is important and fruitful. This kind of combined investigation not only can greatly improve the model-simulation of their kinematics with properly selecting the model parameters (especially the bulk Lorentz factor and Doppler factor as function of time), but also their light curves can be well interpreted in terms of their Doppler-boosting effect. Therefore, we can almost completely (or perfectly) understand the physical nature of these components: their kinematic/dynamic characteristics and emission properties.

We present a multi-wavelength (from far ultraviolet to HI emission) study of star formation feedback on the kinematics of the interstellar medium in the Sculptor Galaxy, NGC253. Its three well-known features (a disrupted stellar disc, a previously reported delining rotation curve, and anomalous HI gas) are studied in a common context of disc asymmetries. About 170 h of on-source ATCA observations are collected and reduced into two versions of \HI\ data cubes of different angular resolution (30'' / 2') and HI column density sensitivity (7.4 $\times$ $10^{19}$cm$^{-2}$ / 4$\times$ $10^{18}$cm$^{-2}$). We separate the anomalous gas from the disc using a custom-made line profile fitting toolkit called FMG. Two star formation tracers (H$\alpha$, FUV emission) are carefully processed and studied. We find that at $R > 7.5~\mathrm{kpc}$ the star formation activity is strongly lopsided (SFR$_{NE}$ >SFR$_{SW}$), and investigate several other properties (H$\alpha$/FUV, dust temperature, stellar age, and disc stability parameters). We also find that the declining nature of the rotation curve perceived by previous studies is not intrinsic but a combined effect of kinematical asymmetries at $R = 7.5$--$16~\mathrm{kpc}$. This is likely the consequence of star formation triggered outflow. The mass distribution and the timescale of the anomalous gas also imply that it originates from gas outflow, which is perhaps caused by galaxy-galaxy interaction considering the crowded environment of NGC253.

Thousands of exoplanet detections have been made over the last twenty-five years using Doppler observations, transit photometry, direct imaging, and astrometry. Each of these methods is sensitive to different ranges of orbital separations and planetary radii (or masses). This makes it difficult to fully characterize exoplanet architectures and to place our solar system in context with the wealth of discoveries that have been made. Here, we use the EXtreme PREcision Spectrograph (EXPRES) to reveal planets in previously undetectable regions of the mass-period parameter space for the star $\rho$ Coronae Borealis. We add two new planets to the previously known system with one hot Jupiter in a 39-day orbit and a warm super-Neptune in a 102-day orbit. The new detections include a temperate Neptune planet ($M{\sin{i}} \sim 20$ M$_\oplus$) in a 281.4-day orbit and a hot super-Earth ($M{\sin{i}} = 3.7$ M$_\oplus$) in a 12.95-day orbit. This result shows that details of planetary system architectures have been hiding just below our previous detection limits; this signals an exciting era for the next generation of extreme precision spectrographs.

Variability studies have proven to be a powerful diagnostic tool for understanding the physics and properties of of Active Galactic Nuclei (AGNs). They provide insights into the spatial and temporal distribution of the emitting regions, the structure and dynamics of the accretion disk, and the properties of the central black hole. Here, we have analysed the K-band spectral variability of the Seyfert 1.9/2 galaxy NGC4388 spanning five epochs over a period of ten years. We have performed spectral synthesis of the nuclear region and found that the contribution of warm dust (T~800K) declined by 88% during these 10 years. In the same period, the [CaVIII] coronal line decreased 61%, whereas BrG emission declined 35%. For the HeI and H2, we did not detect any significant variation beyond their uncertainties. Based on the time span of these changes, we estimate that the region where the warm dust is produced is smaller than 0.6pc, which suggests that this spectral feature comes from the innermost part of the region sampled, directly from the AGN torus. On the other hand, the bulk of [CaVIII] is produced in the inner ~2pc and the nuclear BrG region is more extended, spanning a region larger than 3pc. Lastly, HeI and H2 are even more external, with most of the emission probably being produced in the host galaxy rather than in the AGN. This is the first spectroscopic variability study in the NIR for an AGN where the central source is not directly visible.

Given the widespread availability of grids of models for stellar atmospheres, it is necessary to recover intermediate atmospheric models by means of accurate techniques that go beyond simple linear interpolation and capture the intricacies of the data. Our goal is to establish a reliable, precise, lightweight, and fast method for recovering stellar model atmospheres, that is to say the stratification of mass column, temperature, gas pressure, and electronic density with optical depth given any combination of the defining atmospheric specific parameters: metallicity, effective temperature, and surface gravity, as well as the abundances of other key chemical elements. We employed a fully connected deep neural network which in turn uses a 1D convolutional auto-encoder to extract the nonlinearities of a grid using the ATLAS9 and MARCS model atmospheres. This new method we call iNNterpol effectively takes into account the nonlinearities in the relationships of the data as opposed to traditional machine-learning methods, such as the light gradient boosting method (LightGBM), that are repeatedly used for their speed in well-known competitions with reduced datasets. We show a higher precision with a convolutional auto-encoder than using principal component analysis as a feature extractor.We believe it constitutes a useful tool for generating fast and precise stellar model atmospheres, mitigating convergence issues, as well as a framework for future developments. The code and data for both training and direct interpolation are available online at https://github.com/cwestend/iNNterpol for full reproducibility and to serve as a practical starting point for other continuous 1D data in the field and elsewhere.

We study kilonova emission from binary neutron star (BNS) mergers for the case that a remnant massive neutron star (MNS) forms and collapses to a black hole within $20$ ms after the onset of the merger (which we refer to as "a short-lived case") by consistently employing numerical-relativity and nucleosynthesis results. We find that such kilonovae are fainter and last shorter than those for BNSs resulting in the formation of long-lived ($\gg 1\,{\rm s}$) MNSs, in particular in the optical band. The resulting light curves are too faint and last for a too short duration to explain the kilonova observation for the BNS associated with GW170817, indicating that the merger remnant formed in GW170817 is unlikely to have collapsed to a black hole within a short period of time ($\sim 20$ ms) after the onset of the merger. Our present result implies that early observation is necessary to detect kilonovae associated with BNSs leading to short-lived MNS formation in particular for the optical blue band as well as that kilonovae could be hidden by the gamma-ray burst afterglow for nearly face-on observation. We provide a possible approximate scaling law for near-infrared light curves with the given reference time and magnitude when the decline power of the ${\it z}$-band magnitude, $d M_{\it z}/d{\rm log}_{10}t$, reaches $2.5$. This scaling law suggests that the ${\it HK}$-band follow-up observation should be at least $1$ mag deeper than that for the ${\it z}$-band reference magnitude and earlier than 4 times the reference time.

Escaping exoplanet atmospheres have been observed as deep transit signatures in a few specific spectral lines. Detections have been made in the hydrogen Ly-$\alpha$ line, the metastable helium line at 10830 {\AA} and some UV lines of metallic species. Observational challenges, unexpected non-detections and model degeneracies have generally made it difficult to draw definitive conclusions about the escape process for individual planets. Expanding on the suite of spectral tracers used may help to mitigate these challenges. We present a new framework for modeling the transmission spectrum of hydrodynamically escaping atmospheres. We predict FUV to NIR spectra for systems with different planet and stellar types and identify new lines that can potentially be used to study their upper atmospheres. Measuring the radius in the atmosphere at which the strongest lines form puts them into context within the upper atmospheric structure. Targeting a set of complementary spectral lines for the same planet will help us to better constrain the outflow properties.

We report the discovery of a velocity coherent, funnel shaped ^13CO emission feature in the Galactic centre (GC) using data from the SEDIGISM survey. The molecular cloud appears as a low velocity structure (V_LSR=[-3.5, +3.5] km/s) with an angular extent of 0.95{\deg} x 1{\deg}, extending toward positive Galactic latitudes. The structure is offset from Sgr A* toward negative Galactic longitudes and spatially and morphologically correlates well with the northern lobe of the 430 pc radio bubble, believed to be the radio counterpart of the multiwavelength GC chimney. Spectral line observations in the frequency range of 85-116 GHz have been carried out using the IRAM 30 metre telescope toward 12 positions along the funnel-shaped emission. We examine the ^12C/^13C isotopic ratios using various molecules and their isotopologues. The mean ^12C/^13C isotope ratio (30.6+-2.9) is consistent with the structure located within inner 3 kpc of the Galaxy and possibly in the GC. The velocity of the molecular funnel is consistent with previous radio recombination line measurements of the northern lobe of radio bubble. Our multiwavelength analysis suggests that the funnel shaped structure extending over 100 pc above the Galactic plane is the molecular counterpart of the northern GC chimney.

Optical starlight can be partially polarized while propagating through the dusty, magnetized interstellar medium. The polarization efficiency describes the polarization intensity fraction per reddening unit, P$_V$/E($B-V$), related to the interstellar dust grains and magnetic field properties. The maximum value observed, [P$_V$/E$(B-V)]_{max}$, is thus achieved under optimal polarizing conditions of the interstellar medium. Therefore, the analysis of polarization efficiency observations across the Galaxy contributes to the study of magnetic field topology, small-scale magnetic fluctuations, grain-alignment efficiency, and composition. Infrared observations from $Planck$ satellite have set [P$_V$/E$(B-V)]_{max}$ to 13$\%$ mag$^{-1}$. However, recent optical polarization observations in $Planck$'s highly polarized regions showed polarization efficiency values between 13.6$\%$ mag$^{-1}$ and 18.2$\%$ mag$^{-1}$ (depending on the extinction map used), indicating that [P$_V$/E$(B-V)]_{max}$ is not well constrained yet. We used $V$-band polarimetry of the Interstellar Polarization Survey (consisting of $\sim$10500 high-quality observations distributed in 34 fields of $0.3^{\circ}\times0.3^{\circ}$) to accurately estimate the polarization efficiency in the interstellar medium. We estimated the upper limit of P$_V$/E($B-V$) with the weighted $99th$ percentile of the field. In five regions, the polarization efficiency upper limit is above 13$\%$ mag$^{-1}$. Furthermore, we found [P$_V$/E$(B-V)]_{max} = 15.8^{+1.3}_{-0.9}\%$ mag$^{-1}$ using diffuse intermediate latitude ($|b|>7.5^{\circ}$) regions with apparently strong regular Galactic magnetic field in the plane-of-sky. We studied the variations of P$_V$/E($B-V$) across the sky and tested toy models of polarization efficiency with Galactic longitude that showed some correspondence with a uniform spiral magnetic field.

We investigated the origin of the Planetary Nebula (PN) M 1-16 using narrow band optical imaging, and high and low resolution optical spectra to perform a detailed morphokinematic and chemical studies. M 1-16 is revealed to be a multipolar PN that predominantly emits in [O III] in the inner part of the nebula and [N II] in the lobes. A novel spectral unsharp masking technique was applied to the position-velocity maps (PVs) to reveal a set of multiple structures at the centre of M 1-16 spanning radial velocities from $-40km\,s{-1}$ to $20km\,s{-1}$, with respect to the systemic velocity . The morphokinematic model indicates that the deprojected velocity of the lobe outflows are $\geq100km\,s{-1}$, and particularly the larger lobes and knots have a deprojected velocity of $\simeq350km\,s{-1}$; the inner ellipsoidal component has a deprojected velocity of $\simeq29km\,s{-1}$. A kinematical age of $\sim$8700yr has been obtained from the model assuming an homologous velocity expansion law and a distance of 6.2$\pm$1.9kpc. The chemical analysis indicates that M 1-16 is a Type I PN with a central star of PN (CSPN) mass in the range of $\simeq0.618-0.713$M$\odot$ and an initial mass for the progenitor star between 2.0 and 3.0M$\odot$ (depending on metalicity). An $T_\mathrm{eff}\simeq140\,000$K and log($L/{\rm L}_{\odot})$=2.3 was estimated using the 3MdB photoionization models to reproduce the ionisation stage of the PN. All these results lead us to suggest that M 1-16 is an evolved PN, contrary to the scenario of proto-PN suggested in previous studies. We propose that the mechanism responsible for the morphology of M 1-16 is related to the binary (or multiple star) evolution scenario.

Magnetic fields are now widely recognised as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Models that impose galactic magnetic fields cannot reliably capture their dynamical structure nor interactions with phases within the interstellar medium (ISM). Dynamos must be modelled to create such magnetic fields. ISM models exist in which a small-scale dynamo (SSD) drives a turbulent magnetic field. Others model the large-scale dynamo (LSD) with magnetic field organised at the scale of the disc or spiral arms. Separately, neither can fully describe the dynamics of galactic magnetic fields nor represent their topology. We model the LSD and SSD together at high enough resolution to use explicit Lagrangian resistivity and viscosity. The galactic SSD saturates in less than 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr vs. 1--2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due to ${\alpha}$-quenching near the midplane as the growing mean field produces a magnetic ${\alpha}$ that opposes the kinetic ${\alpha}$. The magnetic energy in our models of the LSD shows slightly sublinear response to increasing resolution, indicating that we are converging towards the physical solution at 1 pc resolution. Including clustered supernovae from OB associations increases the rate of growth of both the SSD and the LSD compared to a horizontally uniform distribution.

We study the subhalo and satellite populations in haloes similar to the Milky Way (MW)-Andromeda paired configuration in the Millennium II and P-Millennium simulations. We find subhaloes are $5\%-15\%$ more abundant in paired haloes than their isolated counterparts that have the same halo mass and large-scale environmental density. Paired haloes tend to reside in a more isotropic environment than isolated haloes, the shear tensor of their large-scale tidal field is possibly responsible for this difference. We also study the thickness of the spatial distribution of the top 11 most massive satellite galaxies obtained in the semi-analytic galaxy sample constructed from the Millennium II simulation. Moreover, satellites that have lost their host subhaloes due to the resolution limit of the simulation have been taken into account. As a result, we find that the difference in the distribution of the satellite thickness between isolated and paired haloes is indistinguishable, which suggests that the paired configuration is not responsible for the observed plane of satellites in the Milky Way. The results in this study indicate the paired configuration could bring some nonnegligible effect on the subhalo abundance in the investigation of the Milky Way's satellite problems.

The mechanism of dust emission from a cometary nucleus is still an open question and thermophysical models have problems reproducing outgassing and dust productions rates simultaneously. In this study, we investigate the capabilities of a rather simple thermophysical model to match observations from Rosetta instruments at comet 67P/Churyumov-Gerasimenko and the influence of model variations. We assume a macro-porous surface structure composed of pebbles and investigate the influence of different model assumptions. Besides the scenario in which dust layers are ejected when the vapour pressure overcomes the tensile strength, we use artificial ejection mechanisms, depending on ice-depletion of layers. We find that dust activity following the pressure criterion is only possible for reduced tensile strength values or reduced gas diffusivity and is inconsistent with observed outgassing rates, because activity is driven by CO$_2$. Only when we assume that dust activity is triggered when the layer is completely depleted in H$_2$O, the ratio of CO$_2$ to H$_2$O outgassing rates is in the expected order of magnitude. However, the dust-to-H$_2$O ratio is never reproduced. Only with decreased gas diffusivity, the slope of the H$_2$O outgassing rate is matched, however absolute values are too low. To investigate maximum reachable pressures, we adapted our model equivalent to a gas-impermeable dust structure. Here, pressures exceeding the tensile strength by orders of magnitude are possible. Maximum activity distances of $3.1 \,\mathrm{au}$, $8.2 \,\mathrm{au}$, and $74 \,\mathrm{au}$ were estimated for H$_2$O-, CO$_2$-, and CO-driven activity of $1 \,\mathrm{cm}$-sized dust, respectively. In conclusion, the mechanism behind dust emission remains unclear.

In this work, we investigated Bayesian methodologies for constraining in the Solar System a Yukawa suppression of the Newtonian potential -- which we interpret as the effect of a non-null graviton mass -- by considering its impact on planetary orbits. Complementary to the previous results obtained with INPOP planetary ephemerides, we consider here a Markov Chain Monte Carlo approach associated with a Gaussian Process Regression for improving the resolution of the constraints driven by planetary ephemerides on the graviton mass in the Solar System. At the end of the procedure, a posterior for the mass of the graviton is presented, providing an upper bound at $1.01 \times 10^{-24} \; eV c^{-2}$ (resp. $\lambda_g \geq 122.48 \times 10^{13} \; km$) with a $99.7\%$ confidence level. The threshold value represents an improvement of 1 order of magnitude relative to the previous estimations. This updated determination of the upper bound is mainly due to the Bayesian methodology, although the use of new planetary ephemerides (INPOP21a used here versus INPOP19a used previously) already induces a gain of a factor 3 with respect to the previous limit. The INPOP21a ephemerides is characterized by the addition of new Juno and Mars orbiter data, but also by a better Solar System modeling, with notably a more realistic model of the Kuiper belt. Finally, by testing the sensitivity of our results to the choice of the $\textit{a priori}$ distribution of the graviton mass, it turns out that the selection of a prior more favorable to zero-mass graviton (that is, here, General Relativity) seems to be more supported by the observations than non-zero mass graviton, leading to a possible conclusion that planetary ephemerides are more likely to favor General Relativity.

In this era of exoplanet characterisation with JWST, the need for a fast implementation of classical forward models to understand the chemical and physical processes in exoplanet atmospheres is more important than ever. Notably, the time-dependent ordinary differential equations to be solved by chemical kinetics codes are very time-consuming to compute. In this study, we focus on the implementation of neural networks to replace mathematical frameworks in one-dimensional chemical kinetics codes. Using the gravity profile, temperature-pressure profiles, initial mixing ratios, and stellar flux of a sample of hot-Jupiters atmospheres as free parameters, the neural network is built to predict the mixing ratio outputs in steady state. The architecture of the network is composed of individual autoencoders for each input variable to reduce the input dimensionality, which is then used as the input training data for an LSTM-like neural network. Results show that the autoencoders for the mixing ratios, stellar spectra, and pressure profiles are exceedingly successful in encoding and decoding the data. Our results show that in 90% of the cases, the fully trained model is able to predict the evolved mixing ratios of the species in the hot-Jupiter atmosphere simulations. The fully trained model is ~1000 times faster than the simulations done with the forward, chemical kinetics model while making accurate predictions.

The binary-driven hypernova (BdHN) model address long gamma-ray bursts (GRBs) associated with type Ic supernovae (SNe) through a series of physical episodes that occur in a binary composed of a carbon-oxygen (CO) star (of mass about 10 solar mass) and a neutron star (NS) companion (of mass about 2 solar mass) in a compact orbit. The SN explosion of the CO star triggers sequence of seven events. The BdHN model has followed the traditional picture of the SN from the CO iron's core collapse. However, the lack of a solution to the problem of producing successful SNe leaves room for alternative scenarios. We here show that tidal synchronization of the CO-NS binary can lead the CO star to critical conditions for fission, hence splitting into two stellar remnants, e.g., about 8.5 solar mass + 1.5 solar mass. We give specific examples of the properties of the products for various orbital periods relevant to BdHNe. The astrophysical consequences of this scenario are outlined.

Since the discovery of the first double neutron star (DNS) system, the number of these exotic binaries has reached fifteen. Here we investigate a channel of DNS formation in binary systems with components above the mass limit of type II supernova explosion (SN II), i.e. 8 MSun. We apply a spherically symmetric homologous envelope expansion model to account for mass loss, and follow the dynamical evolution of the system numerically with a high-precision integrator. The first SN occurs in a binary system whose orbital parameters are pre-defined, then, the homologous expansion model is applied again in the newly formed system. Analysing 1 658 880 models we find that DNS formation via subsequent SN II explosions requires a fine-tuning of the initial parameters. Our model can explain DNS systems with a separation greater than 2.95 au. The eccentricity of the DNS systems spans a wide range thanks to the orbital circularisation effect due to the second SN II explosion. The eccentricity of the DNS is sensitive to the initial eccentricity of the binary progenitor and the orbital position of the system preceding the second explosion. In agreement with the majority of the observations of DNS systems, we find the system centre-of mass velocities to be less than 60 km/s. Neutron stars that become unbound in either explosion gain a peculiar velocity in the range of 0.02 - 240 km/s. In our model, the formation of tight DNS systems requires a post-explosion orbit-shrinking mechanism, possibly driven by the ejected envelopes.

Recent JWST observations of the Type Ia supernova (SN Ia) 2021aefx in the nebular phase have paved the way for late-time studies covering the full optical to mid-infrared (MIR) wavelength range, and with it the hope to better constrain SN Ia explosion mechanisms. We investigate whether public SN Ia models covering a broad range of progenitor scenarios and explosion mechanisms can reproduce the full optical-MIR spectrum of SN 2021aefx at $\sim$270 days post explosion. We perform 1D steady-state non-LTE simulations with the radiative-transfer code CMFGEN, and compare the predicted spectra to SN 2021aefx. The models can explain the main features of SN 2021aefx over the full wavelength range. However, no single model, or mechanism, emerges as a preferred match. We discuss possible causes for the mismatch of the models, including ejecta asymmetries and ionisation effects. Our new calculations of the collisional strengths for Ni III have a major impact on the two prominent lines at 7.35 and 11.00 $\mu$m, and highlight the need for more accurate collisional data for forbidden transitions. Using updated atomic data, we identify a strong feature due to [Ca IV] 3.21 $\mu$m, attributed to [Ni I] in previous studies. We also provide a tentative identification of a forbidden line due to [Ne II] 12.81 $\mu$m, whose peaked profile suggests that neon is mixed inwards during the explosion, as predicted for instance in violent merger models. Contrary to previous claims, we show that the [Ar III] 8.99 $\mu$m line can be broader in sub-$M_\mathrm{Ch}$ models compared to near-$M_\mathrm{Ch}$ models. Our models suggest that key physical ingredients are missing from either the explosion models, or the radiative-transfer post-processing, or both. Nonetheless, they also show the potential of the near- and mid-infrared to uncover new spectroscopic diagnostics of SN Ia explosion mechanisms. [Abridged]

The automated identification of extragalactic objects in large surveys provides reliable and reproducible samples of galaxies in less time than procedures involving human interaction. However, regions near the Galactic disc are more challenging due to the dust extinction. We present the methodology for the automatic classification of galaxies and non-galaxies at low Galactic latitude regions using both images and, photometric and morphological near-IR data from the VVVX survey. Using the VVV-NIRGC, we analyse by statistical methods the most relevant features for galaxy identification. This catalogue was used to train a CNN with image data and an XGBoost model with both photometric and morphological data and then to generate a dataset of extragalactic candidates. This allows us to derive probability catalogues used to analyse the completeness and purity as a function of the configuration parameters and to explore the best combinations of the models. As a test case, we apply this methodology to the Northern disc region of the VVVX survey, obtaining 172,396 extragalatic candidates with probabilities of being galaxies. We analyse the performance of our methodology in the VVV disc, reaching an F1-score of 0.67, a 65 per cent purity and a 69 per cent completeness. We present the VVV-NIR Galaxy Catalogue: Northern part of the Galactic disc comprising 1,003 new galaxies, with probabilities greater than 0.6 for either model, with visual inspection and with only 2 previously identified galaxies. In the future, we intend to apply this methodology to other areas of the VVVX survey.

Based on magnetohydrodynamic turbulence simulations, we generate synthetic synchrotron observations to explore the scaling slope of the underlying MHD turbulence. We propose the new $Q$-$U$ cross intensity $X$ and cross-correlation intensity $Y$ to measure the spectral properties of magnetic turbulence, together with statistics of the traditional synchrotron $I$ and polarization $PI$ intensities. By exploring the statistical behavior of these diagnostics, we find that the new statistics $X$ and $Y$ can extend the inertial range of turbulence to improve measurement reliability. When focusing on different Alfv{\'e}nic and sonic turbulence regimes, our results show that the diagnostics proposed in this paper not only reveal the spectral properties of the magnetic turbulence but also gain insight into the individual plasma modes of compressible MHD turbulence. The synergy of multiple statistical methods can extract more reliable turbulence information from the huge amount of observation data from the Low-Frequency Array for Radio astronomy and the Square Kilometer Array.

Gravitational lensing studies of the Bullet Cluster suggested convincingly in favor of the existence of dark matter. However, it was performed without the knowledge of the original orientation of each galaxy before gravitational lensing. A potential improvement to this issue lies in the measurement of the original orientation from the polarization direction of radio waves emitted from each galaxy. In this context, Francfort et al. derived a formula that can utilize the information about the original orientation of each galaxy to obtain what is called ${\it shear}$. However, we demonstrate that shear in their formula should be replaced by ${\it reduced~shear}$ when the change in sizes of images of galaxies is taken into account. As the previous gravitational lensing analysis of the Bullet Cluster used reduced shear, we suggest applying our improved formula directly for the reanalysis once we obtain the polarization direction of radio waves. In particular, we show that our new formula can yield a more accurate analysis than the previous one, if the polarization direction can be measured more precisely than $10^\circ$.

As one of the nearest and most dormant supermassive black holes (SMBHs), M31* provides a rare but promising opportunity for studying the physics of black hole accretion and feedback at the quiescent state. Previous Karl G. Jansky Very Large Array (VLA) observations with an arcsec resolution have detected M31* as a compact radio source over centimeter wavelengths, but the steep radio spectrum suggests optically-thin synchrotron radiation from an outflow driven by a hot accretion flow onto the SMBH. Aiming to probe the putative radio outflow, we have conducted milli-arcsec-resolution very long baseline interferometric (VLBI) observations of M31* in 2016, primarily at 5 GHz and combining the Very Long Baseline Array, Tianma-65m and Shanghai-25m Radio Telescopes. Despite the unprecedented simultaneous resolution and sensitivity achieved, no significant ($\gtrsim 3\sigma$) signal is detected at the putative position of M31* given an RMS level of $\rm 5.9~\mu Jy\ beam^{-1}$, thus ruling out a point-like source with a peak flux density comparable to that ($\sim30~\mu Jy\ beam^{-1}$) measured by the VLA observations taken in 2012. We disfavor the possibility that M31* has substantially faded since 2012, in view that a 2017 VLA observation successfully detected M31* at a historically-high peak flux density ($\sim75~\mu Jy\ beam^{-1}$ at 6 GHz). Instead, the non-detection of the VLBI observations is best interpreted as the arcsec-scale core being resolved out at the milli-arcsec-scale, suggesting an intrinsic size of M31* at 5 GHz larger than $\sim300$ times the Schwarzschild radius. Such extended radio emission may originate from a hot wind driven by the weakly accreting SMBH.

The Australian SKA Pathfinder (ASKAP) is being used to undertake a campaign to rapidly survey the sky in three frequency bands across its operational spectral range. The first pass of the Rapid ASKAP Continuum Survey (RACS) at 887.5 MHz in the low band has already been completed, with images, visibility datasets, and catalogues made available to the wider astronomical community through the CSIRO ASKAP Science Data Archive (CASDA). This work presents details of the second observing pass in the mid band at 1367.5 MHz, RACS-mid, and associated data release comprising images and visibility datasets covering the whole sky south of declination $+$49$^\circ$. This data release incorporates selective peeling to reduce artefacts around bright sources, as well as accurately modelled primary beam responses. The Stokes I images reach a median noise of 198 $\mu$Jy PSF$^{-1}$ with a declination-dependent angular resolution of 8.1 to 47.5 arcsec that fills a niche in the existing ecosystem of large-area astronomical surveys. We also supply Stokes V images after application of a widefield leakage correction, with a median noise of 165 $\mu$Jy PSF$^{-1}$. We find the residual leakage of Stokes I into V to be $\lesssim$ 0.9 to 2.4 % over the survey. This initial RACS-mid data release will be complemented by a future release comprising catalogues of the survey region. As with other RACS data releases, data products from this release will be made available through CASDA.

Accretion occurs across a large range of scales and physical regimes. Despite this diversity in the physics, the observed properties show remarkably similarity. The theory of propagating fluctuations, in which broad-band variability within an accretion disc travel inwards and combine, has long been used to explain these phenomena. Recent numerical work has expanded on the extensive analytical literature but has been restricted to using the 1D diffusion equation for modelling the disc behaviour. In this work we present a novel numerical approach for 2D (vertically integrated), stochastically driven {\alpha}-disc simulations, generalising existing 1D models. We find that the theory of propagating fluctuations translates well to 2D. However, the presence of epicyclic motion in 2D (which cannot be captured within the diffusion equation) is shown to have an important impact on local disc dynamics. Additionally, there are suggestions that for sufficiently thin discs the log-normality of the light-curves changes. As in previous work, we find that the break frequency in the luminosity power spectrum is strongly dependent on the driving timescale of the stochastic perturbations within the disc, providing a possible observational signature for probing the magnetorotational instability (MRI) dynamo. We also find that thinner discs are significantly less variable than thicker ones, providing a compelling explanation for the greater variability seen in the hard state vs the soft state of X-ray binaries. Finally, we consider the wide-ranging applications of our numerical model for use in other simulations.

We present the expected X-ray polarization signal resulting from distant reprocessing material around black holes. Using a central isotropic power-law emission at the center of the simulated model, we add distant equatorial and axially symmetric media that are covering the central accreting sources. We include partial ionization and partial transparency effects, and the impact of various polarization and steepness of the primary radiation spectrum. The results are obtained with the Monte Carlo code STOKES that considers both line and continuum processes and computes the effects of scattering and absorption inside static homogenous wedge-shaped and elliptical toroidal structures, varying in relative size, composition and distance to the source. We provide first order estimates for parsec-scale reprocessing in Compton-thin and Compton-thick active galactic nuclei, as well as winds around accreting stellar-mass compact objects. The resulting polarization can reach tens of % with either parallel or perpendicular orientation with respect to the axis of symmetry, depending on subtle details of the geometry, density and ionization structure. We also show how principal parameters can be constrained from X-ray spectroscopy or polarimetry in other wavelengths to lift the shown degeneracies in the X-ray band. We provide an application example of the broad modelling discussion by revisiting the recent IXPE 2-8 keV X-ray polarimetric observation of the accreting stellar-mass black hole in Cygnus X-3 from the perspective of partial transparency and ionization of the obscuring outflows.

We investigate the impact of a non-minimal coupling of the scalar field with gravity in inflationary models, where a small coupling is allowed. As a concrete example, we consider the Witten-O'Raifeartaigh model, where, in line with other models, the presence of a coupling strength $\xi$ can recover concordance of the inflationary parameters with cosmic microwave background (CMB) constraints, provided by the Planck collaboration. We go beyond the slow-roll regime and investigate the impact in the description of CMB anisotropies by performing a statistical analysis of the model with the most recent Planck + Baryon Acoustic Oscillations (BAO) data to seek for any indication of a non-zero coupling by data within the model. We find that not only the presence of a non-minimal coupling is seen, but the model has a slight statistical preference when compared with the standard $\Lambda$CDM one. We also discuss the results on the minimally-coupled model, which in general, favours the simple setting where the associated mass scale is equal to the reduced Planck mass $M_p$ while being, in general, disfavored concerning the standard model.

We investigate the use of harmonic analysis techniques to perform measurements of the angular power spectrum on mock pulsar timing data for an isotropic stochastic gravitational-wave background (SGWB) with a dimensionless strain amplitude $A_{\text{gw}}=2 \times 10^{-15}$ and spectral index $\gamma_{\text{gw}}=13/3$. We examine the sensitivity of our harmonic analysis to the number of pulsars (50, 100, and 150) and length of pulsar observation time (10, 20, and 30 years) for an isotropic distribution of pulsars. We account for intrinsic pulsar red noise and use an average value of white noise of ~100 ns. We are able to detect the quadrupole for all our mock harmonic analyses, and for the analysis with 150 pulsars observed for 30 years, we are able to detect up to the $\ell = 5$ multipole. We provide scaling laws for the SGWB amplitude, the quadrupole, and $\ell = 3$ as a function of pulsar observation time and as a function of number of pulsars. We estimate the sensitivity of our harmonic approach to deviations of General Relativity that produce subluminal gravitational wave propagation speeds.

We develop precise analytic description of oscillons - long-lived quasiperiodic field lumps - in scalar field theories with nearly quadratic potentials, e.g. the monodromy potential. Such oscillons are essentially nonperturbative due to large amplitudes, and they achieve extreme longevities. Our method is based on a consistent expansion in the anharmonicity of the potential at strong fields, which is made accurate by introducing a field-dependent "running mass." At every order, we compute effective action for the oscillon profile and other parameters. Comparison with explicit numerical simulations in (3+1)-dimensional monodromy model shows that our method is significantly more precise than other analytic approaches.

A direct detection of black hole formation in neutron star mergers would provide invaluable information about matter in neutron star cores and finite temperature effects on the nuclear equation of state. We study black hole formation in neutron star mergers using a set of 196 numerical relativity simulations consisting of long-lived and black hole-forming remnants. The postmerger gravitational-wave spectrum of a long-lived remnant has greatly reduced power at a frequency $f$ greater than $f_{\rm peak}$, for $f \gtrsim 4\,\rm kHz$, with $f_{\rm peak} \in [2.5, 4]\,\rm kHz$. On the other hand, black-hole-forming remnants exhibit excess power in the same large $f$ region and manifest exponential damping in the time domain characteristic of a quasi-normal mode. We demonstrate that the gravitational-wave signal from a collapsed remnant is indeed a quasi-normal ringing. We report on the opportunity for direct detections of black hole formation with next-generation gravitational-wave detectors such as Cosmic Explorer and Einstein Telescope and set forth the tantalizing prospect of such observations up to a distance of 100 Mpc.

We analyze the isospin equilibration properties of neutrino-transparent nuclear ($npe$) matter in the temperature and density range that is relevant to neutron star mergers. Our analysis incorporates neutrino-transparency corrections to the isospin (``beta'') equilibrium condition which become noticeable at $T\gtrsim 1\,$MeV. We find that the isospin relaxation rate rises rapidly as temperature rises, and at $T\approx 5\,$MeV it is comparable to the timescale of the density oscillations that occur immediately after the merger. This produces a resonant peak in the bulk viscosity at $T\sim 5\,$MeV, which causes density oscillations to be damped on the timescale of the merger. We conclude that there is good reason to include isospin relaxation dynamics in merger simulations.

With the aim of exploring the evidence for or against phase transitions in cold and dense baryonic matter, the inference of the sound speed and equation-of-state for dense matter in neutron stars is extended in view of recent new observational data. The impact of the heavy (2.35 $M_\odot$) black widow pulsar PSR J0952-0607 and of the unusually light supernova remnant HESS J1731-347 is inspected. In addition a detailed re-analysis is performed of the low-density constraint based on chiral effective field theory and of the perturbative QCD constraint at asymptotically high densities, in order to clarify the influence of these constraints on the inference procedure. The trace anomaly measure, $\Delta = 1/3 - P/\varepsilon$, is also computed and discussed. A systematic Bayes factor assessment quantifies the evidence (or non-evidence) of a phase transition within the range of densities realised in the core of neutron stars. One of the consequences of including PSR J0952-0607 in the data base is a further stiffening of the equation-of-state, resulting for a typical 2.1 solar-mass neutron star in a reduced central density of less than five times the equilibrium density of normal nuclear matter. The evidence against the occurrence of a first-order phase transition in neutron star cores is further strengthened.

The new (improved) model of inflation and primordial black hole (PBH) formation is proposed by combining the Starobinsky model of inflation, the Appleby-Battye-Starobinsky (ABS) model of dark energy and a quantum correction in the framework of modified $F(R)$ gravity. The energy scale parameter in the ABS model is taken to be close to the inflationary scale, in order to describe double inflation instead of dark energy. The quantum correction is given by the term quartic in the spacetime scalar curvature $R$ with a negative coefficient $\delta$ in the $F(R)$ function. It is demonstrated that perfect (within $1\sigma$) agreement with current measurements of the cosmic microwave background (CMB) radiation can be achieved by choosing the proper value of $\delta$, thus solving the problem of low values of the tilt of CMB scalar perturbations in the earlier proposed model in arXiv:2205.00603. A large enhancement (large peak) in the power spectrum of scalar perturbations is achieved by fine-tuning the parameters of the model. It is found by numerical analysis that it leads to a formation of asteroid-size PBH with the masses up to $10^{20}$ g, which may form dark matter in the current universe.

Stimulated emission is shown to be robust in stars. Through Bose enhancement this produces quantum states of aligned, monochromatic photons similar to a laser. The probability of creating such states is computed. We show that from the solar corona such quantum states would propagate outside of the solar region and through the Solar System without decoherence. For a $1 {\rm m}^2$ detector at the distance of the Earth from the Sun we estimate rates of such quantum states in the few per second thus potentially detectable. The same process should lead to such quantum states also arriving from stars at interstellar distances.

In this work, we systematically derive the Einstein field equations in general relativity and $f(\mathcal{R})$ gravity, the Tolman-Oppenheimer-Volkoff (TOV) equation, and the expressions for axial and polar Tidal Love Numbers (TLNs) for neutron stars. The derivations are sourced from existing literature and elaborated for the ease of comprehension.

Reducing noises and enhancing signal-to-noise ratios (SNRs) have become critical for designing third-generation gravitational-wave (GW) detectors with a GW strain of less than $10^{-23}$/$\rm \sqrt{Hz}$. In this paper, we propose a potential third-generation GW detector based on autocorrelative weak-value amplification (AWVA) for GW detection with a strain of $h_g =$ $4 \times 10^{-25}$/$\rm \sqrt{Hz}$. In our scheme, a GW event induces a phase difference $\Delta \phi$ by passing through an 11-bounce delay line, 10-km arm-length, zero-area Sagnac interferometer illuminated with a 1064-nm laser. Subsequently, $\Delta \phi$ is amplified as the parameter of post-selection by choosing the appropriate pre-selected state and coupling strength in AWVA. In particular, we theoretically investigate the AWVA measurements for GW detection within the frequency band of 200 Hz $\leq$ $f_g$ $\leq$ 800 Hz, considering Gaussian noises with negative-decibel SNRs. The peak response of the AWVA sensitivity $\kappa(f_g)$ occurs at frequency $f_{g, max}$ = 500 Hz, which falls within the frequency band of interest of the current third-generation GW detectors. Our simulation results indicate that AWVA can demonstrate a measurable sensitivity of $\Theta(f_g)$ within the frequency band of interest. Moreover, the robustness of WVA shows promising potential in mitigating the effects of Gaussian noises.

We exploit the great potential offered by Bayesian Neural Networks (BNNs) to directly decipher the internal composition of neutron stars (NSs) based on their macroscopic properties. By analyzing a set of simulated observations, namely NS radius and tidal deformability, we leverage BNNs as effective tools for inferring the proton fraction and sound speed within NS interiors. To achieve this, several BNNs models were developed upon a dataset of $\sim$ 25K nuclear EoS within a relativistic mean-field framework, obtained through Bayesian inference that adheres to minimal low-density constraints. Unlike conventional neural networks, BNNs possess an exceptional quality: they provide a prediction uncertainty measure. To simulate the inherent imperfections present in real-world observations, we have generated four distinct training and testing datasets that replicate specific observational uncertainties. Our initial results demonstrate that BNNs successfully recover the composition with reasonable levels of uncertainty. Furthermore, using mock data prepared with the DD2, a different class of relativistic mean-field model utilized during training, the BNN model effectively retrieves the proton fraction and speed of sound for neutron star matter.

We derive new constraints on the dilaton parameter appearing in the spherically-symmetric black hole solution of Einstein-Maxwell-dilaton-axion gravity, by studying the geodesic motion of the S2 star in the Galactic Center. Einstein-Maxwell-dilaton-axion black holes represent a compelling alternative to the standard black hole paradigm in General Relativity. This theory emerges from the low energy effective action of the heterotic string theory and has been proven to predict peculiar observational features from the direct imaging of black hole shadows. At a fundamental level, Einstein-Maxwell-dilaton-axion includes additional electromagnetic, dilatonic and axionic fields coupled to the space-time metric. When considering charged non-rotating black hole solutions, the additional fields endow the metric with one extra parameter $b$, called dilaton parameter, that is theoretically bound to $0<b<M$. Using publicly available astrometric data for S2 we derive an upper bound on $b\lesssim 12M$ at 95% confidence level and we demonstrate that only including the measurement of the relativistic orbital precession for S2 is sufficient to reduce this bound to $b\lesssim 1.4M$ at the same confidence level. Additionally, using a mock data mimicking future observations of S2 with the GRAVITY interferometer, we show that improved astrometric precision can help further narrow down the allowed dilaton parameter range to $b\lesssim0.033M$ after monitoring the S2 orbit for one and a half period.