Papers for Monday, Dec 18 2023

Analysing next-generation cosmological data requires balancing accurate modeling of non-linear gravitational structure formation and computational demands. We propose a solution by introducing a machine learning-based field-level emulator, within the Hamiltonian Monte Carlo-based Bayesian Origin Reconstruction from Galaxies (BORG) inference algorithm. Built on a V-net neural network architecture, the emulator enhances the predictions by first-order Lagrangian perturbation theory to be accurately aligned with full $N$-body simulations while significantly reducing evaluation time. We test its incorporation in BORG for sampling cosmic initial conditions using mock data based on non-linear large-scale structures from $N$-body simulations and Gaussian noise. The method efficiently and accurately explores the high-dimensional parameter space of initial conditions, fully extracting the cross-correlation information of the data field binned at a resolution of $1.95h^{-1}$ Mpc. Percent-level agreement with the ground truth in the power spectrum and bispectrum is achieved up to the Nyquist frequency $k_\mathrm{N} \approx 2.79h \; \mathrm{Mpc}^{-1}$. Posterior resimulations - using the inferred initial conditions for $N$-body simulations - show that the recovery of information in the initial conditions is sufficient to accurately reproduce halo properties. In particular, we show highly accurate $M_{200\mathrm{c}}$ halo mass function and stacked density profiles of haloes in different mass bins $[0.853,16]\times 10^{14}M_{\odot}h^{-1}$. As all available cross-correlation information is extracted, we acknowledge that limitations in recovering the initial conditions stem from the noise level and data grid resolution. This is promising as it underscores the significance of accurate non-linear modeling, indicating the potential for extracting additional information at smaller scales.

We present the new Guided Moments ($\texttt{GM}$) formalism for neutrino modeling in astrophysical scenarios like core-collapse supernovae and neutron star mergers. The truncated moments approximation ($\texttt{M1}$) and Monte-Carlo ($\texttt{MC}$) schemes have been proven to be robust and accurate in solving the Boltzmann's equation for neutrino transport. However, it is well-known that each method exhibits specific strengths and weaknesses in various physical scenarios. The $\texttt{GM}$ formalism effectively solves these problems, providing a comprehensive scheme capable of accurately capturing the optically thick limit through the exact $\texttt{M1}$ closure and the optically thin limit through a $\texttt{MC}$ based approach. In addition, the $\texttt{GM}$ method also approximates the neutrino distribution function with a reasonable computational cost, which is crucial for the correct estimation of the different neutrino-fluid interactions. Our work provides a comprehensive discussion of the formulation and application of the $\texttt{GM}$ method, concluding with a thorough comparison across several test problems involving the three schemes ($\texttt{M1}$, $\texttt{MC}$, $\texttt{GM}$) under consideration.

We have determined the proper motions (PMs) of 12 dwarf galaxies in the Local Group (LG), ranging from the outer Milky Way (MW) halo to the edge of the LG. We used HST as the first and Gaia as the second epoch using the GaiaHub software. For Leo A and Sag DIG we also used multi-epoch HST measurements relative to background galaxies. Orbital histories derived using these PMs show that two-thirds of the galaxies in our sample are on first infall with $>$90\% certainty. The observed star formation histories (SFHs) of these first-infall dwarfs are generally consistent with infalling dwarfs in simulations. The remaining four galaxies have crossed the virial radius of either the MW or M31. When we compare their star formation (SF) and orbital histories we find tentative agreement between the inferred pattern of SF with the timing of dynamical events in the orbital histories. For Leo~I, SF activity rises as the dwarf crosses the MW's virial radius, culminating in a burst of SF shortly before pericenter ($\approx1.7$~Gyr ago). The SF then declines after pericenter, but with some smaller bursts before its recent quenching ($\approx0.3$~Gyr ago). This shows that even small dwarfs like Leo~I can hold on to gas reservoirs and avoid quenching for several Gyrs after falling into their host, which is longer than generally found in simulations. Leo~II, NGC~6822, and IC~10 are also qualitatively consistent with this SF pattern in relation to their orbit, but more tentatively due to larger uncertainties.

In this work, we constrain the star-forming properties of all possible sites of incipient high-mass star formation in the Milky Way's Galactic Center. We identify dense structures using the CMZoom 1.3mm dust continuum catalog of objects with typical radii of $\sim$0.1pc, and measure their association with tracers of high-mass star formation. We incorporate compact emission at 8, 21, 24, 25, and 70um from MSX, Spitzer, Herschel, and SOFIA, catalogued young stellar objects, and water and methanol masers to characterize each source. We find an incipient star formation rate (SFR) for the CMZ of ~0.08 Msun yr^{-1} over the next few 10^5 yr. We calculate upper and lower limits on the CMZ's incipient SFR of ~0.45 Msun yr^{-1} and ~0.05 Msun yr^{-1} respectively, spanning between roughly equal to and several times greater than other estimates of CMZ's recent SFR. Despite substantial uncertainties, our results suggest the incipient SFR in the CMZ may be higher than previously estimated. We find that the prevalence of star formation tracers does not correlate with source volume density, but instead ~75% of high-mass star formation is found in regions above a column density ratio (N_{SMA}/N_{Herschel}) of ~1.5. Finally, we highlight the detection of ``atoll sources'', a reoccurring morphology of cold dust encircling evolved infrared sources, possibly representing HII regions in the process of destroying their envelopes.

Recent measurements of the $4$-point correlation functions (4PCF) from spectroscopic surveys provide evidence for parity-violations in the large-scale structure of the Universe. If physical in origin, this could point to exotic physics during the epoch of inflation. However, searching for parity-violations in the 4PCF signal relies on a large suite of simulations to perform a rank test, or an accurate model of the 4PCF covariance to claim a detection, and this approach is incapable of extracting parity information from the higher-order $N$-point functions. In this work we present an unsupervised method which overcomes these issues, before demonstrating the approach is capable of detecting parity-violations in a few toy models using convolutional neural networks. This technique is complementary to the 4-point method and could be used to discover parity-violations in several upcoming surveys including DESI, Euclid and Roman.

Gravitational lensing may render individual high-mass stars detectable out to cosmological distances, and several extremely magnified stars have in recent years been detected out to redshifts $z\approx 6$. Here, we present Muspelheim, a model for the evolving spectral energy distributions of both metal-enriched and metal-free stars at high redshifts. Using this model, we argue that lensed stars should form a highly biased sample of the intrinsic distribution of stars across the Hertzsprung-Russell diagram, and that this bias will typically tend to favour the detection of lensed stars in evolved stages characterized by low effective temperatures, even though stars only spend a minor fraction of their lifetimes in such states. We also explore the prospects of detecting individual, lensed metal-free (Population III) stars at high redshifts using the James Webb Space Telescope (JWST). We find that very massive ($\gtrsim 100\ M_\odot$) Population III stars at $z\gtrsim 6$ may potentially be detected by JWST in surveys covering large numbers of strong lensing clusters, provided that the Population III stellar initial mass function is sufficiently top-heavy, that these stars evolve to effective temperatures $\leq 15000$ K, and that the cosmic star formation rate density of Pop III stars reaches $\gtrsim 10^{-4}\ M_\odot$ cMpc$^{-3}$ yr$^{-1}$ at $z\approx$ 6-10. Various ways to distinguish metal-free lensed stars from metal-enriched ones are also discussed.

We introduce the Wide-field Retrieval of Astrodata Program (WRAP), a tool created to aid astronomers in gathering photometric and astrometric data for point sources that may confuse simple cross-matching algorithms because of their faintness or motion. WRAP allows astronomers to correctly cross-identify objects with proper motion across multiple surveys by wedding the catalog data with its underlying images, thus providing visual confirmation of cross-associations in real time. Developed within the Backyard Worlds: Planet 9 citizen science project, WRAP aims to aid in the characterization of faint, high motion sources by this collaboration (and others).

We study the evolution of isolated self-interacting dark matter (SIDM) halos that undergo gravothermal collapse and are driven deep into the short-mean-free-path regime. We assume spherical Navarro-Frenk-White (NFW) halos as initial conditions and allow for elastic dark matter self-interactions. We discuss the structure of the halo core deep in the core-collapsed regime and how it depends on the particle physics properties of dark matter, in particular, the velocity dependence of the self-interaction cross section. We find an approximate universality deep in this regime that allows us to connect the evolution in the short- and long-mean-free-path regimes, and approximately map the velocity-dependent self-interaction cross sections to constant ones for the full gravothermal evolution. We provide a semi-analytic prescription based on our numerical results for halo evolution deep in the core-collapsed regime. Our results are essential for estimating the masses of the black holes that are likely to be left in the core of SIDM halos.

The properties of galaxies in low-density regions of the universe suggest an interplay between galaxy formation and environment. However, the specific reason why this particular large-scale environment influences the evolution of galaxies remains unclear. This paper examines the properties and evolutionary paths of galaxies within cosmic voids using the Illustris TNG300 simulation. The population of void galaxies at z = 0 has a higher star formation rate, a smaller stellar-to-halo-mass ratio, higher gas metallicity, and lower stellar metallicity in comparison with non-void galaxies at fixed stellar mass. Our analysis shows that these differences are mainly due to the characteristics of galaxies classified as satellites, for which the largest differences between void and non-void samples are found. Although the mean number of mergers is similar between void and non-void samples at a fixed stellar mass, void galaxies tend to experience mergers at later times, resulting in a more recent accumulation of accreted stellar mass. While the mean net accreted mass is comparable for high mass galaxies, low mass void galaxies tend to exhibit higher fractions of accreted stars than non-void galaxies. This finding challenges the common notion that void galaxies predominantly experience growth with infrequent mergers or interactions.

The goal of this work is to obtain a Hubble constant estimate through the study of the quadruply lensed, variable QSO SDSSJ1433+6007. To achieve this we combine multi-filter, archival $\textit{HST}$ data for lens modelling and a dedicated time delay monitoring campaign with the 2.1m Fraunhofer telescope at the $\textit{Wendelstein Observatory}$. The lens modelling is carried out with the public $\texttt{lenstronomy}$ Python package for each of the filters individually. Through this approach, we find that the data in one of the $\textit{HST}$ filters (F160W) contain a light contaminant, that would, if remained undetected, have severely biased the lensing potentials and thus our cosmological inference. After rejecting these data we obtain a combined posterior for the Fermat potential differences from the lens modelling in the remaining filters (F475X, F814W, F105W and F140W) with a precision of $\sim6\%$. The analysis of the $\textit{g'}$-band Wendelstein light curve data is carried out with a free-knot spline fitting method implemented in the public Python $\texttt{PyCS3}$ tools. The precision of the time delays between the QSO images has a range between 7.5 and 9.8$\%$ depending on the brightness of the images and their time delay. We then combine the posteriors for the Fermat potential differences and time delays. Assuming a flat $\Lambda$CDM cosmology, we infer a Hubble parameter of $H_0=76.6^{+7.7}_{-7.0}\frac{\mathrm{km}}{\mathrm{Mpc\;s}}$, reaching $9.6\%$ uncertainty for a single system.

The Davis-Chandrasekhar-Fermi (DCF) method is widely used to evaluate magnetic fields in star-forming regions. Yet it remains unclear how well DCF equations estimate the mean plane-of-the-sky field strength in a map region. To address this question, five DCF equations are applied to an idealized cloud map. Its polarization angles have a normal distribution with dispersion ${\sigma}_{\theta}$,and its density and velocity dispersion have negligible variation. Each DCF equation specifies a global field strength $B_{DCF}$ and a distribution of local DCF estimates. The "most-likely" DCF field strength $B_{ml}$ is the distribution mode (Chen et al. 2022), for which a correction factor ${\beta}_{ml}$ = $B_{ml}$/$B_{DCF}$ is calculated analytically. For each equation ${\beta}_{ml}$ < 1, indicating that $B_{DCF}$ is a biased estimator of $B_{ml}$. The values of ${\beta}_{ml}$ are ${\beta}_{ml}\approx$ 0.7 when $B_{DCF} \propto {{\sigma}_{\theta}}^{-1}$ due to turbulent excitation of Afv\'enic MHD waves, and ${\beta}_{ml}\approx$ 0.9 when $B_{DCF} \propto {{\sigma}_{\theta}}^{-1/2}$ due to non-Alfv\'enic MHD waves. These statistical correction factors ${\beta}_{ml}$ have partial agreement with correction factors ${\beta}_{sim}$ obtained from MHD simulations. The relative importance of the statistical correction is estimated by assuming that each simulation correction has both a statistical and a physical component. Then the standard, structure function, and original DCF equations appear most accurate because they require the least physical correction. Their relative physical correction factors are 0.1, 0.3, and 0.4 on a scale from 0 to 1. In contrast the large-angle and parallel-${\delta}B$ equations have physical correction factors 0.6 and 0.7. These results may be useful in selecting DCF equations, within model limitations.

Catalina Sky Survey (CSS) is a major survey of Near-Earth Objects (NEOs). In a recent work, we used CSS observations from 2005-2012 to develop a new population model of NEOs (NEOMOD). CSS's G96 telescope was upgraded in 2016 and detected over 10,000 unique NEOs since then. Here we characterize the NEO detection efficiency of G96 and use G96's NEO detections from 2013-2022 to update NEOMOD. This resolves previous model inconsistencies related to the population of large NEOs. We estimate there are 936+/-29 NEOs with absolute magnitude H<17.75 (diameter D>1 km for the reference albedo p_V=0.14). The slope of the NEO size distribution for H=25-28 is found to be relatively shallow (cumulative index 2.6) and the number of H<28 NEOs (D>9 m) is determined to be (1.20+/-0.04)x10^7. Small NEOs have a different orbital distribution and higher impact probabilities than large NEOs. We estimate 0.034+/-0.002 impacts of H<28 NEOs on the Earth per year, which is near the low end of the impact flux range inferred from atmospheric bolide observations. Relative to a model where all NEOs are delivered directly from the main belt, the population of small NEOs detected by G96 shows an excess of low-eccentricity orbits with a=1--1.6 au that appears to increase with H. We suggest that the population of very small NEOs is boosted by tidal disruption of large NEOs during close encounters to the terrestrial planets. When the effect of tidal disruption is (approximately) accounted for in the model, we estimate 0.06+/-0.01 impacts of H<28 NEOs on the Earth per year, which is more in line with the bolide data.

In this paper, we calibrate the luminosity relation of gamma-ray bursts (GRBs) with the machine learning (ML) methods for reconstructing distance-redshift relation from the Pantheon+ sample of type Ia supernovae (SNe Ia) in a cosmology-independent way. The A219 GRB data set at low redshift are used to calibrate the Amati relation ($E_{\rm p}$-$E_{\rm iso}$) relation by the ML methods selected with the best performance %and the calibrated results are extrapolated to the high redshift data to construct the Hubble diagram at high redshift. We constrain cosmological models via the Markov Chain Monte Carlo numerical method with the GRBs at high redshift and the latest observational Hubble data (OHD). By the K-Nearest Neighbors (KNN) methods, we obtain $\Omega_{\rm m}$ = $0.29^{+0.09}_{-0.06}$, $h$ = $0.66^{+0.04}_{-0.07}$ , $w_0$ = $-0.83^{+0.66}_{-0.31}$, $w_a$ = $-0.91^{+0.87}_{-0.46}$ at 1-$\sigma$ confidence level for the Chevallier-Polarski-Linder (CPL) model in a flat space, which favor the dark energy with a possible evolution ($w_a\neq0$) at 1-$\sigma$. These results are consistent with those obtained from GRBs calibrated via the Gaussian Process.

Recent attempts at measuring the variation of $c$ using an assortment of standard candles and the redshift-dependent Hubble expansion rate inferred from the currently available catalog of cosmic chronometers have tended to show that the speed of light appears to be constant, at least up to $z\sim 2$. A notable exception is the use of high-redshift UV $+$ X-ray quasars, whose Hubble diagram seems to indicate a $\sim 2.7\sigma$ deviation of c from its value $c_0$ ($\equiv 2.99792458 \times 10^{10}$ cm s$^{-1}$) on Earth. We show in this paper, however, that this anomaly is due to an error in the derived relation between the luminosity distance, $D_L$, and $H(z)$ when $c$ is allowed to vary with redshift, and an imprecise calibration of the quasar catalog. When these deficiences are addressed correctly, one finds that $c/c_0=0.95 \pm 0.14$ in the redshift range $0\lesssim z\lesssim 2$, fully consistent with zero variation within the measurement errors.

We study the effects of general relativity (GR) on the evolution and alignment of circumbinary disks around binaries on all scales. We implement relativistic apsidal precession of the binary into the hydrodynamics code {\sc phantom}. We find that the effects of GR can suppress the stable polar alignment of a circumbinary disk, depending on how the relativistic binary apsidal precession timescale compares to the disk nodal precession timescale. Studies of circumbinary disk evolution typically ignore the effects of GR which is an appropriate simplification for low mass or widely separated binary systems. In this case, polar alignment occurs providing that the disks initial misalignment is sufficiently large. However, systems with a very short relativistic precession timescale cannot polar align and instead move toward coplanar alignment. In the intermediate regime where the timescales are similar, the outcome depends upon the properties of the disk. Polar alignment is more likely in the wavelike disk regime (where the disk viscosity parameter is less than the aspect ratio, $\alpha<H/r$) since the disk is in good radial communication. In the viscous disk regime disk breaking is more likely. Multiple rings can destructively interact with one another resulting in short disk lifetimes, and the disk moving towards coplanar alignment. Around main-sequence star or stellar mass black hole binaries, polar alignment may be suppressed far from the binary but in general the inner parts of the disk can align to polar. Polar alignment may be completely suppressed for disks around supermassive black holes for close binary separations.

We present the discovery of a new H I structure in the NGC 7194 group from the observations using the Karl G. Jansky Very Large Array. NGC 7194 group is a nearby (z ~ 0.027) small galaxy group with five quiescent members. The observations reveal a 200 kpc-long H I plume that spans the entire group with a total mass of M$_{HI}$ = 3.4 x 10$^{10}$ M$_{\odot}$. The line-of-sight velocity of the H I gas gradually increases from south (7200 km s$^{-1}$) to north (8200 km $^{-1}$), and the local velocity dispersion is up to 70 km s$^{-1}$. The structure is not spatially coincident with any member galaxies but it shows close associations with a number of blue star-forming knots. Intragroup H I gas is not rare, but this particular structure is still one of the unusual cases in the sense that it does not show any clear connection with sizable galaxies in the group. We discuss the potential origins of this large-scale H I gas in the NGC 7194 group and its relation with the intergalactic star-forming knots. We propose that this HI feature could have originated from tidal interactions among group members or the infall of a late-type galaxy into the group. Alternatively, it might be leftover gas from flyby intruders.

Continuum reverberation mapping with high-cadence, long-term UV/optical monitoring of Active Galactic Nuclei (AGNs) enables us to resolve the AGN central engine sizes on different timescales. The frequency-resolved time lags of NGC 5548 (the target for the AGN STORM I campaign) are inconsistent with the X-ray reprocessing of the classical Shakura $\&$ Sunyaev disk model. Here we show that the frequency-resolved time lags in NGC 5548 can be well produced by the Corona-Heated Accretion-disk Reprocessing (CHAR) model. Moreover, we make the CHAR model predictions of the frequency-resolved time lags for Mrk 817, the source of the AGN STORM II campaign. We also obtain the frequency-resolved time lags as a function of the black-hole mass and Eddington ratio, which is valid for black-hole masses from $10^{6.5}$ to $10^9\ M_{\odot}$, and Eddington ratios from 0.01 to 1. Moreover, we demonstrate that, with the time spans of current continuum reverberation-mapping campaigns, the lag-luminosity relation of the CHAR model can be $\tau_{\mathrm{gz}}\propto L_{\mathrm{5100}}^{0.55\pm0.04}$, which is consistent with observations. Future observations can test our results and shed new light on resolving the AGN central engine.

We present the results of the spectroscopic study of a chemically peculiar star HD 152564. Using medium-resolution ($R=37000$) observations obtained with the HRS spectrograph mounted on the South African Large Telescope we determined atmospheric parameters T$_{eff}$=11950$\pm200$ K and log g=3.6$\pm0.2$ dex. Abundance analysis revealed mild deficiency of the light elements and an overabundance of up to $\sim$2 dex of metals with greatest excess for the silicon. With these characteristics HD 152564 is a typical member of the silicon subgroup of Ap stars. Rotational modulation of the light curve and line profiles of HD 152564 are typical for inhomogeneous surface distribution of elements in its atmosphere. We performed multi-element Doppler imaging of the HD 152564 surface. Abundance maps constructed for He, O, Mg, Si, and Fe revealed concentration of these elements in a sequence of equatorial spots as well as in the circumpolar rings. The photometric maximum of the light curve coincided with the visibility of two most overabundant silicon spots. Abundances determined from the different ionisation stages of Fe and Si show clear evidence for vertical stratification of these elements in HD 152564 atmosphere. Meanwhile, the horizontal distribution of silicon reconstructed from the lines of different ionisation stages and excitation energies appeared to be identical with increasing average abundance deeper in atmosphere.

The black hole X-ray binary source 4U 1543--47 experienced a super-Eddington outburst in 2021, reaching a peak flux of up to $\sim1.96\times10^{-7}\rm erg\ \rm cm^{-2}\ \rm s^{-1}$ ($\sim 8.2$ Crab) in the 2--10\,keV band. Soon after the outburst began, it rapidly transitioned into the soft state. Our goal is to understand how the accretion disk structure deviates from a standard thin disk when the accretion rate is near Eddington. To do so, we analyzed spectra obtained from quasi-simultaneous observations conducted by the Hard X-ray Modulation Telescope (Insight-HXMT), the Nuclear Spectroscopic Telescope Array (NuSTAR), and the Neil Gehrels Swift Observatory (Swift). These spectra are well-fitted by a model comprising a disk, a weak corona, and a reflection component. We suggest that the reflection component is caused by disk self-irradiation, that is by photons emitted from the inner disk which return to the accretion disk surface, as their trajectories are bent by the strong gravity field. In this scenario, the best-fitting parameters imply that the reflected flux represents more than half of the total flux. Using general relativistic ray-tracing simulations, we show that this scenario is viable when the disk becomes geometrically thick, with a funnel-like shape, as the accretion rate is near or above the Eddington limit. In the specific case of 4U 1543--47, an angle $\gtrsim$ 45 deg between the disk surface and the equatorial plane can explain the required amount of self-irradiation.

Radio-loud active galactic nuclei (RLAGNs) are a unique AGN population and were thought to be preferentially associated with supermassive black holes (SMBHs) at low accretion rates. They could impact the host galaxy evolution by expelling cold gas through the jet-mode feedback. In this work, we studied CO(6-5) line emission in a high-redshift radio galaxy, MRC 0152-209, at z=1.92 using ALMA up to a $0.024''$-resolution (corresponding to ~200 pc). This system is a starburst major merger constituted of two galaxies: the northwest (NW) one hosting the RLAGN with jet kinetic power $L_{\rm jet}\gtrsim2\times10^{46}$ erg/s and the southeast (SE) one. Based on the SED fitting for the entire system (NW+SE galaxies), we found AGN bolometric luminosity $L_{\rm AGN,bol}\sim0.9-3\times10^{46}$ erg/s for the RLAGN. We estimated BH mass through $M_{\rm BH}-M_\star$ scaling relations and found an Eddington ratio of $\sim0.7-4$ conservatively. These results suggest that the RLAGN is radiatively efficient and the powerful jets could be launched from a super-Eddington accretion disc. ALMA reveals a massive ($M_{\rm H_2}\sim2\times10^9$ Msun), compact ($\sim500$ pc), and lopsided molecular outflow perpendicular to the jet axis. The mass outflow rate (~1200-2600 Msun/yr) is comparable with the star formation rate of ~2000-3000 Msun/yr. The outflow kinetic power/$L_{\rm AGN,bol}$ ratio of ~0.008-0.02 and momentum boost factor ~3-24 agree with the radiative-mode AGN feedback. On the other hand, the jets can also drive the molecular outflow within its lifetime of $\sim2\times10^5$ yr without additional energy supply from AGN radiation. The jets then could remove all cold gas from the host galaxy through long-term, episodic launching. Our study reveals a unique object where starburst, powerful jets, and rapid BH growth co-exist, which may represent a fundamental stage of AGN-host galaxy co-evolution.

Variations in chemical abundances with evolutionary phase have been identified among stars in globular and open clusters with a wide range of metallicities. In the metal-poor clusters, these variations compare well with predictions from stellar structure and evolution models considering the internal diffusive motions of atoms and ions, collectively known as atomic diffusion, when moderated by an additional mixing process with a fine-tuned efficiency. We present here an investigation of these effects in the Galactic globular cluster NGC 6121 (M4) ([Fe/H] = -1.13) through a detailed chemical abundance analysis of 86 stars using high-resolution ESO/VLT FLAMES spectroscopy. The stars range from the main-sequence turnoff point (TOP) to the red giant branch (RGB) just above the bump. We identify C-N-O and Mg-Al-Si abundance anti-correlations, and confirm the presence of a bimodal population differing by 1 dex in nitrogen abundance. The composition of the second-generation stars imply pollution from both massive (20-40 Msol) and asymptotic giant branch stars. We find evolutionary variations in chemical abundances between the TOP and RGB, which are robust to uncertainties in stellar parameters and modelling assumptions. The variations are weak, but match predictions well when employing efficient additional mixing. Without correcting for Galactic production of lithium, we derive an initial lithium abundance 2.63+-0.10, which is marginally lower than the predicted primordial BBN value.

Nearly continuous, densely sampled, space-based photometry allows us to recover the finest details in the light variations of stars. The number of such light curves have been increasing rapidly in the last few years thanks to the extended mission of the Kepler space telescope and the launch of the TESS mission. This new era brings us new perspectives in RR Lyrae and Cepheid studies, where low amplitude phenomena can be studied in a wide range of individual stars and on a statistical basis. In this proceedings I review the recent investigations of the Kepler, K2 and TESS fields, as well as the challenges in accurately reducing high-quality photometry.

We present 24 new dense lightcurves of the near-Earth asteroids (3103) Eger, (161989) Cacus, (2100) Ra-Shalom and (12711) Tukmit, obtained with the Instituto Astrof\'isico Canarias 80 and Telescopio Abierto Remoto 2 telescopes at the Teide Observatory (Tenerife, Spain) during 2021 and 2022, in the framework of projects visible NEAs observations survey and NEO Rapid Observation, Characterization and Key Simulations. The shape models and rotation state parameters ($P$, $\lambda$, $\beta$) were computed by applying the lightcurve inversion method to the new data altogether with the archival data. For (3013) Eger and (161989) Cacus, our shape models and rotation state parameters agree with previous works, though they have smaller uncertainties. For (2100) Ra-Shalom, our results also agree with previous studies. Still, we find that a Yarkovsky - O'Keefe - Radzievskii - Paddack acceleration of $\upsilon = (0.223\pm0.237)\times10^{-8}$ rad d$^{-2}$ slightly improves the fit of the lightcurves, suggesting that (2100) Ra-Shalom could be affected by this acceleration. We also present for the first time a shape model for (12711) Tukmit, along with its rotation state parameters ($P=3.484900 \pm 0.000031$ hr, $\lambda = 27^{\circ}\pm 8^{\circ}$, $\beta = 9^{\circ} \pm 15^{\circ}$).

The possibility of watching the Universe expand in real time and in a model-independent way, first envisaged by Allan Sandage more than 60 years ago and known as the redshift drift, is within reach of forthcoming astrophysical facilities, particularly the Extremely Large Telescope (ELT) and the Square Kilometre Array Observatory (SKAO). The latter, probing lower redshifts, enables us to watch the Universe's acceleration era in real time, while the former does the same for the matter era. We use Fisher Matrix Analysis techniques, which we show to give comparable results to those of a Markov Chain Monte Carlo approach, to discuss forecasts for SKAO measurements of the redshift drift and their cosmological impact. We consider specific fiducial cosmological models but mainly rely on a more agnostic cosmographic series (which includes the deceleration and jerk parameters), and we also discuss prospects for measurements of the drift of the drift. Overall, our analysis shows that SKAO measurements, with a reasonable amount of observing time, can provide a competitive probe of the low-redshift accelerating Universe.

We investigate non-gravitational signals of dark energy within the framework of gauge symmetry in the dark energy sector. Traditionally, dark energy has been primarily studied through gravitational effects within general relativity or its extensions. On the other hand, the gauge principles have played a central role in the standard model sector and dark matter sector. If the dark energy field operates under a gauge symmetry, it introduces the possibility of studying all major components of the present universe under the same gauge principle. This approach marks a significant shift from conventional methodologies, offering a new avenue to explore dark energy.

Context. The determination of accurate photometric redshifts (photo-zs) in large imaging galaxy surveys is key for cosmological studies. One of the most common approaches are machine learning techniques. These methods require a spectroscopic or reference sample to train the algorithms. Attention has to be paid to the quality and properties of these samples since they are key factors in the estimation of reliable photo-zs. Aims. The goal of this work is to calculate the photo-zs for the Y3 DES Deep Fields catalogue using the DNF machine learning algorithm. Moreover, we want to develop techniques to assess the incompleteness of the training sample and metrics to study how incompleteness affects the quality of photometric redshifts. Finally, we are interested in comparing the performance obtained with respect to the EAzY template fitting approach on Y3 DES Deep Fields catalogue. Methods. We have emulated -- at brighter magnitude -- the training incompleteness with a spectroscopic sample whose redshifts are known to have a measurable view of the problem. We have used a principal component analysis to graphically assess incompleteness and to relate it with the performance parameters provided by DNF. Finally, we have applied the results about the incompleteness to the photo-z computation on Y3 DES Deep Fields with DNF and estimated its performance. Results. The photo-zs for the galaxies on DES Deep Fields have been computed with the DNF algorithm and added to the Y3 DES Deep Fields catalogue. They are available at https://des.ncsa.illinois.edu/releases/y3a2/Y3deepfields. Some techniques have been developed to evaluate the performance in the absence of "true" redshift and to assess completeness. We have studied... (Partial abstract)

We present new high-sensitivity e-MERLIN and VLA radio images of the prototypical Seyfert 2 galaxy NGC 1068 at 5, 10 and 21 GHz. We image the radio jet, from the compact components NE, C, S1 and S2 to the faint double-lobed jet structure of the NE and SW jet lobes. Furthermore, we map the jet between by combining e-MERLIN and VLA data for the first time. Components NE, C and S2 have steep spectra indicative of optically-thin non-thermal emission domination between 5 and 21 GHz. Component S1, which is where the AGN resides, has a flat radio spectrum. We report a new component, S2a, a part of the southern jet. We compare these new data with the MERLIN and VLA data observed in 1983, 1992 and 1995 and report a flux decrease by a factor of 2 in component C, suggesting variability of this jet component. With the high angular resolution e-MERLIN maps, we detect the bow shocks in the NE jet lobe that coincide with the molecular gas outflows observed with ALMA. The NE jet lobe has enough radio power considered to be responsible for driving out the dense molecular gas observed with ALMA around the same region.

A semi-empirical IR/Vis line list, SOLIS, for the sulphur monoxide molecule $^{32}$S$^{16}$O is presented. SOLIS includes accurate empirical rovibrational energy levels, uncertainties, lifetimes, quantum number assignments, and transition probabilities in the form of Einstein $A$ coefficients covering the $X\,{}^{3}\Sigma^{-}$, $a\,{}^{1}\Delta^{ }$, $b\,{}^{1}\Sigma^{+}$, $A\,{}^{3}\Pi$, $B\,{}^{3}\Sigma^{-}$, $X\,{}^{\prime\prime3}\Sigma^{+}$, $A\,{}^{\prime 3}\Delta$ and $e\,{}^{1}\Pi$ systems and wavenumber range up to 43303.5 cm$^{-1}$ ($\geq 230.93$ nm) with $J\le 69$. SOLIS has been computed by solving the rovibronic Schr\"{o}dinger equation for diatomics using the general purpose variational code Duo and starting from a published ab initio spectroscopic model of SO (including potential energy curves, coupling curves, (transition) dipole moment curves) which is refined to experimental data. To this end, a database of 50106 experimental transitions, 48972 being non-redundant, has been compiled through the analysis of 29 experimental sources, and a self-consistent network of 8558 rovibronic energy levels for the $X$, $a$, $b$, $A$, $B$, and $C$ electronic states has been generated with the MARVEL algorithm covering rotational and vibrational quantum numbers $J \leq 69$ and $v \leq 30$ and energies up to 52350.40 cm$^{-1}$. No observed transitions connect to the $ B\,{}^{3}\Sigma^{-} (v = 0)$ state which is required to model perturbations correctly, so we leave fitting the $B\,{}^3\Sigma^-$ and $C\,{}^3\Pi$ state UV model to a future project. The SO line list is available at ExoMol from www.exomol.com.

The structure of the Small Magellanic Cloud (SMC) outside of its main body is characterised by tidal branches resulting from its interactions mainly with the Large Magellanic Cloud (LMC). Characterising the stellar populations in these tidal components helps to understand the dynamical history of this galaxy and of the Magellanic system in general. We provide full phase-space vector information for Southern Bridge clusters. We performed a photometric and spectroscopic analysis of twelve SMC clusters, doubling the number of SMC clusters with full phase-space vector information known to date. We reclassify the sample considering 3D distances and 3D velocities. We found that some of the clusters classified as Southern Bridge objects according to the projected 2D classification actually belong to the Main Body and Counter-Bridge in the background. The comparison of the kinematics of the genuine foreground Bridge clusters with those previously analysed in the same way reveals that Southern Bridge clusters are moving towards the LMC and share the kinematics of the Northern Bridge. Adding to our sample clusters from the literature with CaT metallicity determinations we compare the age-metallicity relation of the Southern Bridge with the one of the Northern Bridge. We reinforce the idea that both regions do not seem to have experienced the same chemical enrichment history and that there is a clear absence of clusters in the Northern Bridge older than 3Gyr and more metal-poor than -1.1, which would not seem to be due to a selection effect.

Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems.

Context. Cassiopeia A occupies an important place among supernova remnants (SNRs) in low-frequency radio astronomy. The analysis of its continuum spectrum from low frequency observations reveals the evolution of the SNR absorption properties over time and suggests a method for probing unshocked ejecta and the SNR interaction with the circumstellar medium (CSM). Aims. In this paper we present low-frequency measurements of the integrated spectrum of Cassiopeia A to find the typical values of free-free absorption parameters towards this SNR in the middle of 2023. We also add new results to track its slowly evolving and decreasing integrated flux density. Methods. We used the New Extension in Nan\c{c}ay Upgrading LOFAR (NenuFAR) and the Ukrainian Radio Interferometer of NASU (URAN-2, Poltava) for measuring the continuum spectrum of Cassiopeia A within the frequency range of 8-66 MHz. The radio flux density of Cassiopeia A has been obtained on June-July, 2023 with two sub-arrays for each radio telescope, used as a two-element correlation interferometer. Results. We measured magnitudes of emission measure, electron temperature and an average number of charges of the ions for both internal and external absorbing ionized gas towards Cassiopeia A from its integrated spectrum. Generally, their values are comparable to those presented by Stanislavsky et al. (2023), but their slight changes show the evolution of free-free absorption parameters in this SNR. Based on high accuracy of the measurements, we have detected the SNR-CSM interaction. Conclusions. The integrated flux-density spectrum of Cassiopeia A obtained with the NenuFAR and URAN-2 interferometric observations opens up new possibilities for continuous monitoring the ionized gas properties in and around Cassiopeia A to observe theevolution of unshocked ejecta and the SNR-CSM interaction in future studies.

XTE J1946+274 is a Be/X-ray binary with a 15.8s spin period and 172 d orbital period. Using $\textit{RXTE/PCA}$ data of the 1998 outburst, a cyclotron line around 37 keV was reported. The presence of this line, its dependence on the pulse phase, and its variation with luminosity have been of some debate since. In this work, we present the reanalysis of two $\textit{AstroSat}$ observations: one made during the rising phase of the 2018 outburst and the other during the declining phase of the 2021 outburst. We also present a new analysis of the $\textit{Insight}$-HXMT observations of the source at the peak of the 2018 outburst. We find the source to be spinning up over the course of the outburst and spinning down between the two outbursts. We report the presence of a higher cyclotron line energy using the 2018 $\textit{AstroSat}$ observation ($\sim 45$ keV) and 2018 $\textit{Insight}$-HXMT observation ($\sim$ 50 keV) and a line at $\sim$ 40 keV during the declining phase of the 2021 outburst using data from $\textit{AstroSat}$. We also investigate the pulse phase dependence of the cyclotron line parameters and find that the line is significantly detected in all the phases of both $\textit{AstroSat}$ observations, along with showing variation with the pulse phase. This differs from the previous results reported using $\textit{BeppoSAX}$ and $\textit{NuSTAR}$. We explain this behaviour of the cyclotron line to be due to photon spawning and different accretion column radii at the two poles of this neutron star.

We present the first results from "Surveying the Whirlpool at Arcseconds with NOEMA" (SWAN), an IRAM Northern Extended Millimetre Array (NOEMA)+30m large program that maps emission from several molecular lines at 90 and 110 GHz in the iconic nearby grand-design spiral galaxy M~51 at cloud-scale resolution ($\sim$3\arcsec=125\,pc). As part of this work, we have obtained the first sensitive cloud-scale map of N$_2$H$^+$(1-0) of the inner $\sim5\,\times 7\,$kpc of a normal star-forming galaxy, which we compare to HCN(1-0) and CO(1-0) emission to test their ability in tracing dense, star-forming gas. The average N$_2$H$^+$-to-HCN line ratio of our total FoV is $0.20\pm0.09$, with strong regional variations of a factor of $\gtrsim 2$ throughout the disk, including the south-western spiral arm and the center. The central $\sim1\,$kpc exhibits elevated HCN emission compared to N$_2$H$^+$, probably caused by AGN-driven excitation effects. We find that HCN and N$_2$H$^+$ are strongly super-linearily correlated in intensity ($\rho_\mathrm{Sp}\sim 0.8$), with an average scatter of $\sim0.14\,$dex over a span of $\gtrsim 1.5\,$dex in intensity. When excluding the central region, the data is best described by a power-law of exponent $1.2$, indicating that there is more N$_2$H$^+$ per unit HCN in brighter regions. Our observations demonstrate that the HCN-to-CO line ratio is a sensitive tracer of gas density in agreement with findings of recent Galactic studies which utilize N$_2$H$^+$. The peculiar line ratios present near the AGN and the scatter of the power-law fit in the disk suggest that in addition to a first-order correlation with gas density, second-order physics (such as optical depth, gas temperature) or chemistry (abundance variations) are encoded in the N$_2$H$^+$/CO, HCN/CO and N$_2$H$^+$/HCN ratios.

Machine learning has rose to become an important research tool in the past decade, its application has been expanded to almost if not all disciplines known to mankind. Particularly, the use of machine learning in astrophysics research had a humble beginning in the early 1980s, it has rose and become widely used in many sub-fields today, driven by the vast availability of free astronomical data online. In this short review, we narrow our discussion to a single topic in astrophysics - the estimation of photometric redshifts of galaxies and quasars, where we discuss its background, significance, and how machine learning has been used to improve its estimation methods in the past 20 years. We also show examples of some recent machine learning photometric redshift work done in Malaysia, affirming that machine learning is a viable and easy way a developing nation can contribute towards general research in astronomy and astrophysics.

We present an updated analysis of eleven cosmological models that may help reduce the Hubble tension, which now reaches the $6\sigma$ level when considering the latest SH0ES measurement versus recent CMB and BAO data, assuming $\Lambda$CDM. Specifically, we look at five classical extensions of $\Lambda$CDM (with massive neutrinos, spatial curvature, free-streaming or self-interacting relativistic relics, or dynamical dark energy) and six elaborate models featuring either a time-varying electron mass, early dark energy or some non-trivial interactions in the neutrino sector triggered by a light Majoron. We improve over previous works in several ways. We include the latest data from the South Pole Telescope as well as the most recent measurement of the Hubble rate by the SH0ES collaboration. We treat the summed neutrino mass as a free parameter in most of our models, which reveals interesting degeneracies and constraints. We define additional metrics to assess the potential of a model to reduce or even solve the Hubble tension. We validate an emulator that uses active learning to train itself during each parameter inference run for any arbitrary model. We find that the time-varying electron mass and the Majoron models are now ruled out at more than $3\sigma$. Models with a time-varying electron mass and spatial curvature or with early dark energy reduce the tension to 1.0-2.9$\sigma$. Nevertheless, none of the models considered in this work is favored with enough statistical significance to become the next concordance model of Cosmology.

In X-ray observations of hard state black hole X-ray binaries, rapid variations in accretion disc and coronal power-law emission are correlated and show Fourier-frequency-dependent time lags. On short (~0.1 s) time-scales, these lags are thought to be due to reverberation and therefore may depend strongly on the geometry of the corona. Low-frequency quasi-periodic oscillations (QPOs) are variations in X-ray flux that have been suggested to arise because of geometric changes in the corona, possibly due to General Relativistic Lense-Thirring precession. Therefore one might expect the short-term time lags to vary on the QPO time-scale. We performed novel spectral-timing analyses on NICER observations of the black hole X-ray binary MAXI J1820+070 during the hard state of its outburst in 2018 to investigate how the short-term time lags between a disc-dominated and a coronal power-law-dominated energy band vary on different time-scales. Our method can distinguish between variability due to the QPO and broadband noise, and we find a linear correlation between the power-law flux and lag amplitude that is strongest at the QPO frequency. We also introduce a new method to resolve the QPO signal and determine the QPO-phase-dependence of the flux and lag variations, finding that both are very similar. Our results are consistent with a geometric origin of QPOs, but also provide evidence for a dynamic corona with a geometry varying in a similar way over a broad range of time-scales, not just the QPO time-scale.

Aims. We aim to investigate the stellar population properties, ages, and metal content of the globular clusters (GCs) in NGC 3311, the central galaxy of the Hydra I cluster, to better constrain its evolution history. Methods. We used integral-field spectroscopic data from the Multi-Unit Spectroscopic Explorer (MUSE) to identify 680 sources in the central region of NGC 3311 and extract their 1D spectra. An analysis of these sources in terms of morphologies, radial velocities, and emission lines allowed us to narrow down our selection to 49 bona fide GC candidates. We split these candidates into two groups depending on their projected distance to the galaxy center (R), namely inner (R<20arcsec) and outer (R>20arcsec) GCs. We stacked the extracted 1D spectra of the inner and outer GC populations to increase the signal-to-noise ratios (S/Ns) of the resulting spectra and hence allow full-spectrum fitting. In addition, we also created a stacked spectrum of all GCs in NGC 3311 and one of the two most central GC candidates. Using the code pPXF, we performed a stellar population analysis on the four stacked 1D spectra, obtaining mass-weighted integrated ages, metallicities, and [$\alpha$/Fe] abundances. Results. All GCs are old, with ages >13.5 Gyr, and they have super-solar metallicities. Looking at the color distribution, we find that the inner ones tend to be redder and more metal rich than the outer ones. This is consistent with the two-phase formation scenario. Looking at the full-spectral fitting results, at face value the outer GCs have a larger [$\alpha$/Fe] ratio, which is in line with what is found for the stars that dominate the surface brightness profile at the same radii. However, the values for outer and inner GCs are consistent within the uncertainties.

The disagreement between low- and high-redshift measurements of the Hubble parameter is emerging as a serious tension that challenges the standard model of cosmology. Here we focus not on the tension itself, but on a model-independent view of the Hubble parameter in a general spacetime. We define a covariant Hubble parameter that is in principle observable on the past lightcone. Then we present a careful analysis of the transformation of the Hubble parameter under a change of observer. To leading order in relative velocity, a moving observer will measure the same monopole and quadrupole as a matter-frame observer, but will also detect a dipole, generated by Doppler and aberration effects. The aberrated contribution to the dipole is neglected in a standard perturbative analysis. In order to connect the Hubble parameter to observables, we develop a covariant analysis of cosmic distances and their behaviour under boosts. This reveals that the standard relation between the Hubble parameter and luminosity distance is not invariant under boosts and holds only in the matter frame. A moving observer should correct the luminosity distance by a redshift factor -- otherwise an incorrect dipole and a spurious octupole are detected. In a standard perturbative analysis, this error leads to a false prediction of {\em no} dipole.

Hot subdwarf B (sdB) stars are helium core burning stars that have lost almost their entire hydrogen envelope due to binary interaction. Their assumed canonical mass of $\rm M_{\mathrm{sdB}}\sim0.47 M_{\odot}$ has recently been debated given a broad range found both from observations as well as from the simulations. Here, we revise and refine the mass range for sdBs derived two decades ago with the Eggleton code, using the stellar evolution code MESA, and discuss the effects of metallicity and the inclusion of core overshooting during the main sequence. We find an excellent agreement for low-mass progenitors, up to $\sim2.0 \rm M_{\odot}$. For stars more massive than $\sim2.5 \rm M_{\odot}$ we obtain a wider range of sdB masses compared to the simulations from the literature. Our MESA models for the lower metallicity predict, on average, slightly more massive sdBs. Finally, we show the results for the sdB lifetime as a function of sdB mass and discuss the effect this might have in the comparison between simulations and observational samples. This study paves the way for reproducing the observed Galactic mass distribution of sdB binaries.

The accurate reconstruction of Cosmic Microwave Background (CMB) maps and the measurement of its power spectrum are crucial for studying the early universe. In this paper, we implement a convolutional neural network to apply the Wiener Filter to CMB temperature maps, and use it intensively to compute an optimal quadratic estimation of the power spectrum. Our neural network has a UNet architecture as that implemented in WienerNet, but with novel aspects such as being written in python 3 and TensorFlow 2. It also includes an extra channel for the noise variance map, to account for inhomogeneous noise, and a channel for the mask. The network is very efficient, overcoming the bottleneck that is typically found in standard methods to compute the Wiener Filter, such as those that apply the conjugate gradient. It scales efficiently with the size of the map, making it a useful tool to include in CMB data analysis. The accuracy of the Wiener Filter reconstruction is satisfactory, as compared with the standard method. We heavily use this approach to efficiently estimate the power spectrum, by performing a simulation-based analysis of the optimal quadratic estimator. We further evaluate the quality of the reconstructed maps in terms of the power spectrum and find that we can properly recover the statistical properties of the signal. We find that the proposed architecture can account for inhomogeneous noise efficiently. Furthermore, increasing the complexity of the variance map presents a more significant challenge for the convergence of the network than the noise level does.

We determine necessary and sufficient conditions under which a large class of relativistic generalizations of Braginskii's magnetohydrodynamics, described using Israel-Stewart theory, are causal and strongly hyperbolic in the fully nonlinear regime in curved spacetime. Our new nonlinear analysis provides stricter constraints on the dynamical variables that cannot be obtained via a standard linear expansion around equilibrium. Causality severely constrains the size of shear-viscous corrections, placing a bound on the far-from-equilibrium dynamics of magnetized weakly collisional relativistic plasmas, which rules out the onset of the firehose instability in such systems.

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (\NGDEEP) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). To date, the James Webb Space Telescope (JWST) has observed 50$\%$ of the allocated orbits (Epoch 1) of this program (\NGDEEPA). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, \NGDEEPA\ has achieved a 3$\sigma$ sensitivity limit of 2 $\times$ 10$^{-18}$ erg/s/cm$^2$. We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the dust content, star-formation rates, and metallicity ([OIII]/H$\beta$) of galaxies at $1<z<3.5$. Further, we identify the presence of active galactic nuclei (AGN) and infer the mass of their supermassive black holes (SMBHs) using broadened restframe MgII and H$\beta$ emission lines. The spectroscopic results are then compared with the physical properties of galaxies extrapolated from fitting spectral energy distribution (SED) models to photometry alone. The results clearly demonstrate the unique power and efficiency of WFSS at near-infrared wavelengths over other methods to determine the properties of galaxies across a broad range of redshifts.

Context. Accreting ultracompact binaries contain a white dwarf that is accreting from a degenerate object and have orbital periods shorter than 65 minutes. Aims. The aims of this letter are to report the discovery and the orbital period of four new eclipsing accreting ultracompact binaries found using the Zwicky Transient Facility, and to discuss their photometric properties. Methods. We searched through a list of 4171 dwarf novae compiled using the Zwicky Transient Facility and used the Box Least Square method to search for periodic signals in the data. Results. We found four new eclipsing accreting ultracompact binaries with orbital periods between 25.9-56 minutes, one of which is previously published as an AM CVn, while the other three systems are new discoveries. The other two shorter period systems are likely also AM CVn systems, while the longest period system with a period of 56 minutes shows multiple super-outbursts observed in two years which is more consistent with it being a Helium-CV.

The lack of power anomaly is an unexpected feature observed at large angular scales in the CMB maps by the COBE, WMAP and Planck satellites. This signature, which consists in a missing of power with respect to that predicted by the $\Lambda$CDM model, might hint at a new cosmological phase before the standard inflationary era. The main point of this paper is taking the latest Planck polarisation data into account to investigate how CMB polarisation improves the understanding of this feature. With this aim, we apply to the last Planck data, both PR3 (2018) and PR4 (2020) releases, a new class of estimators able to evaluate this anomaly considering temperature and polarisation data both separately and in a jointly way. This is the first time that the PR4 dataset is used to study this anomaly. In order to critically evaluate this feature, taking into account the residuals of known systematic effects present in the Planck datasets, we analyse the cleaned CMB maps using different combinations of sky masks, harmonic range and binning on the CMB multipoles. Our analysis shows that the estimator based only on temperature data confirms the presence of a lack of power with a lower-tail-probability (LTP), depending on the component separation method, $\leq 0.33\%$ and $\leq 1.76\%$, for PR3 and PR4 respectively. To our knowledge the $LTP \leq 0.33\%$ for the PR3 dataset is the lowest one present in the literature obtained from Planck 2018 data considering the Planck confidence mask. We find significant differences between these two datasets when polarisation is taken into account. However, we also show that for the PR3 dataset the inclusion of the subdominant polarisation information provides estimates which are less likely accepted in a $\Lambda$CDM cosmological model than the only-temperature analysis on the whole harmonic-range considered.

AGB stars are major contributors to the chemical enrichment of the ISM through nucleosynthesis and extensive mass loss. Most of our current knowledge of AGB atmospheric and circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC+10216. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB CSEs by studying a sample of three carbon stars, IRAS15194-5115, IRAS15082-4808, and IRAS07454-7112, and test the archetypal status often attributed to IRC+10216. We performed spatially resolved, unbiased spectral surveys in ALMA Band 3. We identify a total of 132 rotational transitions from 49 molecular species. There are two main morphologies of the brightness distributions: centrally-peaked (e.g. HCN) and shell-like (e.g. C$_2$H). We estimated the sizes of the molecular emitting regions using azimuthally-averaged radial profiles of the line brightness distributions, and derived abundance estimates. Of the shell distributions, the cyanopolyynes peak at slightly smaller radii than the hydrocarbons, and CN and HNC show the most extended emission. The emitting regions for each species are the smallest for IRAS07454-7112. We find that, within the uncertainties of the analysis, the three stars present similar abundances for most species, also compared to IRC+10216. We find that SiO is more abundant in our three stars compared to IRC+10216. Our estimated isotopic ratios match well the literature values for the sources. The observed circumstellar chemistry appears very similar across our sample and compared to that of IRC+10216, both in terms of the relative location of the emitting regions and molecular abundances. This implies that, to a first approximation, the chemical models tailored to IRC+10216 are able to reproduce the observed chemistry in C-rich envelopes across roughly an order of magnitude in wind density.

Magnetars form a special class of neutron stars possessing superstrong magnetic fields and demonstrating power flares triggered likely by these fields. Observations of such flares reveal the presence of quasi-periodic oscillations (QPOs) at certain frequencies; they are thought to be excited in the flares. QPOs carry potentially important information on magnetar structure, magnetic field, and mechanisms of magnetar activity. We calculate frequencies of torsional (magneto-elastic) oscillations of the magnetar crust treating the magnetic field effects in the first order of perturbation theory. The theory predicts splitting of non-magnetic oscillation frequencies into Zeeman components. Zeeman splitting of torsional oscillation spectrum of magnetars was suggested, clearly described and estimated by Shaisultanov and Eichler (2009) but their work has not been given considerable attention. To extend it we suggest the technique of calculating oscillation frequencies including Zeeman splitting at not too strong magnetic fields for arbitrary magnetic field configuration. Zeeman splitting enriches the oscillation spectrum and simplifies theoretical interpretation of observations. We calculate several low-frequency oscillations of magnetars with pure dipole magnetic field in the crust. The results qualitatively agree with low-frequency QPOs detected in the hyperflare of SGR 1806--20, and in the giant flare of SGR 1900+14.

The secular Laplace-Lagrange orbital solution, decomposing eccentricities into a set of uniformly precessing eigenmodes is a classical result that is typically solved numerically. However, in the limit where orbits are closely spaced, several simplifications make it possible to make analytical progress. We derive simple expressions for the eccentricity eigenmodes in a co-planar 3-planet system where the middle planet is massless, and show that these approximate the true eigenmodes of more general systems with 3 massive planets in various limits. These results provide intuition for the secular dynamics of real systems, and have applications for understanding the stability boundary for compact multi-planet systems.

We characterize the earliest galaxy population in the JADES Origins Field (JOF), the deepest imaging field observed with JWST. We make use of the ancillary Hubble optical images (5 filters spanning $0.4-0.9\mu\mathrm{m}$) and novel JWST images with 14 filters spanning $0.8-5\mu\mathrm{m}$, including 7 medium-band filters, and reaching total exposure times of up to 46 hours per filter. We combine all the imaging data at $>2\mu\mathrm{m}$ to construct the deepest imaging ever taken at these wavelengths, reaching as deep as $\approx31.4$ AB mag in the stack and 30.1-30.8 AB mag ($5\sigma$, $r=0.1"$ circular aperture) in individual filters. We measure photometric redshifts and use robust selection criteria to identify a sample of eight galaxy candidates at redshifts $z=11.5-15$. These objects show compact half-light radii of $R_{1/2}\sim50-200$pc, stellar masses of $M_\star\sim10^7-10^8M_\odot$, and star-formation rates of $\mathrm{SFR}\sim0.1-1~M_\odot~\mathrm{yr}^{-1}$. Our search finds no candidates at $15<z<20$, placing upper limits at these redshifts. We develop a forward modeling approach to infer the properties of the evolving luminosity function without binning in redshift or luminosity that marginalizes over the photometric redshift uncertainty of our candidate galaxies and incorporates the impact of non-detections. We find a $z=12$ luminosity function in good agreement with prior results, and that the luminosity function normalization and UV luminosity density decline by a factor of $\sim2.5$ from $z=12$ to $z=14$. We discuss the possible implications of our results in the context of theoretical models for evolution of the dark matter halo mass function.

In the string axiverse scenario, primordial black holes (PBHs) can sustain non-negligible spin parameters as they evaporate. We show that tracking both the mass and spin evolution of a PBH in its final hour can yield a purely gravitational probe of new physics beyond the TeV scale, allowing one to determine the number of new scalars, fermions, vector bosons, and spin-3/2 particles. Furthermore, we propose a multi-messenger approach to accurately measure the mass and spin of a PBH from its Hawking photon and neutrino primary emission spectra, which is independent of putative interactions between the new degrees of freedom and the Standard Model particles, as well as from the Earth-PBH distance.

We investigate gravitational-wave backgrounds (GWBs) of primordial origin that would manifest only at ultra-high frequencies, from kilohertz to 100 gigahertz, and leave no signal at either LIGO, Einstein Telescope, Cosmic Explorer, LISA, or pulsar-timing arrays. We focus on GWBs produced by cosmic strings and make predictions for the GW spectra scanning over high-energy scale (beyond $10^{10}$ GeV) particle physics parameters. Signals from local string networks can easily be as large as the Big Bang nucleosynthesis/cosmic microwave background bounds, with a characteristic strain as high as $10^{-26}$ in the 10 kHz band, offering prospects to probe grand unification physics in the $10^{14}-10^{17}$ GeV energy range. In comparison, GWB from axionic strings is suppressed (with maximal characteristic strain $\sim 10^{-31}$) due to the early matter era induced by the associated heavy axions. We estimate the needed reach of hypothetical futuristic GW detectors to probe such GW and, therefore, the corresponding high-energy physics processes. Beyond the information of the symmetry-breaking scale, the high-frequency spectrum encodes the microscopic structure of the strings through the position of the UV cutoffs associated with cusps and/or kinks, as well as potential information about friction forces on the string. The IR slope, on the other hand, reflects the physics responsible for the decay of the string network. We discuss possible strategies for reconstructing the scalar potential, particularly the scalar self-coupling, from the measurement of the UV cutoff of the GW spectrum.

The primordial black holes (PBHs) as all-dark matter (DM) hypothesis has recently been demotivated by the prediction that these objects would source an excessive rate of fast radio bursts (FRBs). However, these predictions were based on several simplifying assumptions to which this rate is highly sensitive. In this article, we improve previous estimates of this rate arising from the capture of PBHs by neutron stars (NSs), aiming to revitalise this theory. We more accurately compute the velocity distribution functions of PBHs and NSs and also consider an enhancement in the NS and DM density profiles at galactic centers due to the presence of a central supermassive black hole. We find that previous estimates of the rate of FRBs sourced by the capture of PBHs by NSs were 3 orders of magnitude too large, concluding that the PBHs as all DM hypothesis remains a viable theory and that the observed FRB rate can only be entirely explained when considering a central, sufficiently spiky PBH density profile.

500,000 to 1 million satellites are expected in the next decades, primarily to build internet constellations called megaconstellations. These megaconstellations are disposable and will constantly re-enter and be replaced, hence creating a layer of conductive particulate. Here it will be shown that the mass of the conductive particles left behind from worldwide distribution of re-entry satellites is already billions of times greater than the mass of the Van Allen Belts. From a preliminary analysis, the Debye length in spaceflight regions is significantly higher than non-spaceflight regions according to CCMC ionosphere data. As the megaconstellations grow, the Debye length of the satellite particulate may exceed that of the cislunar environment and create a conductive layer around the earth worldwide. Thus, satellite reentries may create a global band of plasma dust with a charge higher than the rest of the magnetosphere. Therefore, perturbation of the magnetosphere from conductive satellites and their plasma dust layer should be expected and should be a field of intensive research. Human activity is not only impacting the atmosphere, it is clearly impacting the ionosphere.

We compute first and second-order bulk-viscous transport properties due to weak-interaction processes in $npe$ matter in the neutrino transparent regime. The transport coefficients characterize the out-of-beta-equilibrium pressure corrections, which depend on the weak-interaction rates and the equation of state. We calculate these coefficients for realistic equations of state and show they are sensitive to changes in the nuclear symmetry energy $J$ and its slope $L$.

Massive fields can imprint unique oscillatory features on primordial correlation functions or inflationary correlators, which is dubbed the cosmological collider signal. In this work, we analytically investigate the effects of a time-dependent mass of a scalar field on inflationary correlators, extending previous numerical studies and implementing techniques developed in the cosmological bootstrap program. The time-dependent mass is in general induced by couplings to the slow-roll inflaton background, with particularly significant effects in the case of non-derivative couplings. By linearly approximating the time dependence, the mode function of the massive scalar is computed analytically, on which we derive analytic formulae for two-, three-, and four-point correlators with the tree-level exchange of the massive scalar. The obtained formulae are utilized to discuss the phenomenological impacts on the power spectrum and bispectrum, and it is found that the scaling behavior of the bispectrum in the squeezed configuration, i.e., the cosmological collider signal, is modified from a time-dependent Boltzmann suppression. By investigating the scaling behavior in detail, we are in principle able to determine the non-derivative couplings between the inflaton and the massive particle.

We discuss an extra-dimensional braneworld with a 5th dimension compactified on a circle. As a characteristic feature, the warp factor is hyperbolic and separates the hidden and visible branes by a bulk horizon without a singularity. The two most widely separated scales of 4D physics - the 4D Planck mass and 4D cosmological constant - are determined by two physical scales in the extra dimension, namely: $(i)$ the proper size of the extra dimension, $R$, and, $(ii)$ the distance between the visible brane and the horizon, $R_0$. A realistic scale hierarchy between 4D Planck mass and 4D cosmological constant is obtained for $R/R_0\sim2.34$. The usual fine tuning is not reduced but promoted to a fine tuning of two separate brane energy densities that must approach the fundamental scale of the model with very high precision. Our scenario is based on an exact solution to the 5D Einstein equations with a strictly empty bulk and Friedmann-Lema\^itre-Robertson-Walker metric on the 4D branes. This requires positive 4D brane energy densities and describes an adiabatic runaway solution in agreement with the de Sitter swampland conjecture. The Kaluza-Klein (KK) graviton states are solutions of a modified P\"oschl-Teller potential which permits a discrete graviton spectrum of $\textit{exactly two}$ modes. In addition to the usual massless graviton, our scenario predicts an extra massive spin-2 graviton with a mass gap of $m_1=\sqrt{2}H_0\approx2\times10^{-33}\,\mathrm{eV}$ which might be detectable in the foreseeable future. A KK tower of gravitons, or a possible continuum of massive graviton states, is prohibited by unitarity with respect to the horizon. We discuss hurdles in turning this model into a realistic cosmology at all times, which points us towards 4D brane tensions that that must be raising towards the fundamental scale of the model, while the observable 4D expansion rate is decreasing.

We present dynamos computed using a hybrid QG-3D numerical scheme in a thick spherical shell geometry. Our model is based on a quasi-geostrophic convection code extended with a 3D treatment of heat transport and magnetic induction. We find a collection of self-sustained, multipolar, weak field dynamos with magnetic energy one or two orders of magnitude lower than the kinetic energy. The poloidal magnetic energy is weak and, by construction, there is a lack of equatorially anti-symmetric components in the Buoyancy and Lorentz forces. This leads to configurations where the velocity field is only weakly impacted by the magnetic field, similar to dynamos found in 3D simulations where zonal flows and the $\Omega$-effect dominate. The time-dependence of these dynamos is characterised by quasi-periodic oscillations that we attribute to dynamo waves. The QG-3D dynamos found so far are not Earth-like. The inability of our setup to produce strong, dipole-dominated, magnetic fields likely points to a missing ingredient in our QG flows, and a related lack of helicity and $\alpha$-effect. The models presented here may be more relevant for studying stellar dynamos where zonal flows are known to dominate. This study was carried out at modest control parameters, however moving to lower Ekman numbers, when smaller values of both the magnetic and hydrodynamic Prandtl numbers can be of interest, our approach will be able to gain in efficiency by using relatively coarse grids for the 3D magnetic and temperature fields and a finer grid for the QG velocity field.

The bubble expansion velocity is an important parameter in the prediction of gravitational waves from first order phase transitions. This parameter is difficult to compute, especially in phase transitions in strongly coupled theories. In this work, we present a method to estimate the wall velocity for phase transitions with a large enthalpy jump, valid for weakly and strongly coupled theories. We find that detonations are disfavored in this limit, but wall velocities are not necessarily small. We also investigate the effect of two other features in the equation of state: non-conformal sound speeds and a limited range of temperatures for which the phases exist. We find that the former can increase the wall velocity for a given nucleation temperature, and the latter can restrict the wall velocities to small values. To test our approach, we use holographic phase transitions, which typically display these three features. We find excellent agreement with numerically obtained values of the wall velocity. We also demonstrate that the implications for gravitational waves can be significant.

Does deconfined cold quark matter occur in nature? This is currently one of the fundamental open questions in nuclear astrophysics. In these proceedings, I review the current state-of-the-art techniques to address this question in a model-agnostic manner, by synthesizing inputs from astrophysical observations of neutron stars and their binary mergers, and first-principles calculations within nuclear and particle theory. I highlight recent improvements in perturbative calculations in asymptotically dense cold quark matter, as well as compelling evidence for a conformalizing transition within the cores of massive neutron stars.

We study a group field theory (GFT) for quantum gravity coupled to four massless scalar fields, using these matter fields to define a (relational) coordinate system. We exploit symmetries of the GFT action, in particular under shifts in the values of the scalar fields, to derive a set of classically conserved currents, and show that the same conservation laws hold exactly at the quantum level regardless of the choice of state. We propose a natural interpretation of the conserved currents which implies that the matter fields always satisfy the Klein-Gordon equation in GFT. We then observe that in our matter reference frame, the same conserved currents can be used to extract all components of an effective GFT spacetime metric. Finally, we apply this construction to the simple example of a spatially flat homogeneous and isotropic universe. Our proposal goes substantially beyond the GFT literature in which only specific geometric quantities such as the total volume or volume perturbations could be defined, opening up the possibility to study more general geometries as emerging from GFT.

In this work we introduce PhenomXO4a, the first phenomenological, frequency-domain gravitational waveform model to incorporate multipole asymmetries and precession angles tuned to numerical relativity. We build upon the modeling work that produced the PhenomPNR model and incorporate our additions into the IMRPhenomX framework, retuning the coprecessing frame model and extending the tuned precession angles to higher signal multipoles. We also include, for the first time in frequency-domain models, a recent model for spin-precession-induced multipolar asymmetry in the coprecessing frame to the dominant gravitational-wave multipoles. The accuracy of the full model and its constituent components is assessed through comparison to numerical relativity and numerical relativity surrogate waveforms by computing mismatches and performing parameter estimation studies. We show that, for the dominant signal multipole, we retain the modeling improvements seen in the PhenomPNR model. We find that the relative accuracy of current full IMR models varies depending on location in parameter space and the comparison metric, and on average they are of comparable accuracy. However, we find that variations in the pointwise accuracy do not necessarily translate into large biases in the parameter estimation recoveries.