Papers for Monday, Feb 19 2024

The ALMA interferometer has played a key role in revealing a new component of the Sun-like star forming process: the molecular streamers, i.e. structures up to thousands of au long funneling material non-axisymmetrically to disks. In the context of the FAUST ALMA LP, the archetypical VLA1623-2417 protostellar cluster has been imaged at 1.3 mm in the SO(5$_6$--4$_5$), SO(6$_6$--5$_5$), and SiO(5--4) line emission at the spatial resolution of 50 au. We detect extended SO emission, peaking towards the A and B protostars. Emission blue-shifted down to 6.6 km s$^{-1}$ reveals for the first time a long ($\sim$ 2000 au) accelerating streamer plausibly feeding the VLA1623 B protostar. Using SO, we derive for the first time an estimate of the excitation temperature of an accreting streamer: 33$\pm$9 K. The SO column density is $\sim$ 10$^{14}$ cm$^{-2}$, and the SO/H$_2$ abundance ratio is $\sim$ 10$^{-8}$. The total mass of the streamer is 3 $\times$ 10$^{-3}$ $Msun$, while its accretion rate is 3--5 $\times$ 10$^{-7}$ Msun yr$^{-1}$. This is close to the mass accretion rate of VLA1623 B, in the 0.6--3 $\times$ 10$^{-7}$ Msun yr$^{-1}$ range, showing the importance of the streamer in contributing to the mass of protostellar disks. The highest blue- and red-shifted SO velocities behave as the SiO(5--4) emission, the latter species detected for the first time in VLA1623-2417: the emission is compact (100-200 au), and associated only with the B protostar. The SO excitation temperature is $\sim$ 100 K, supporting the occurrence of shocks associated with the jet, traced by SiO.

One of the most frequently used indicators to characterize the magnetic field's influence on star formation is the relation between magnetic field strength and gas density (B-$\rho$ relation), usually expressed as $B \propto \rho^{\kappa}$. The value of $\kappa$ is an indication of the dynamical importance of the magnetic field during gas compression. Investigating the global magnetic field's impact on this relation and its evolution, we conduct MHD simulations of Milky-Way-like galaxies including gravity, star formation, and supernova feedback along with non-equilibrium chemistry up to $H_2$ formation fueling star formation. Two initial magnetic field morphologies are studied: one completely ordered (toroidal) and the other completely random. In these models, we study the dynamical importance of the magnetic field through the plasma $\beta$ and the B-$\rho$ relation. For both magnetic morphologies, low-density regions are thermally supported, while high-density regions are magnetically dominated. Equipartition is reached earlier and at lower densities in the toroidal model. However, the B-$\rho$ relation is not unique even within the same galaxy, as it consistently includes two different branches for a given density, with $\kappa$ ranging from about 0.2 to 0.8. The mean value of $\kappa$ for each model also displays significant variations over time, which supersede the differences between the two models. While our findings suggest that the magnetic field morphology does influence the galactic B-$\rho$ relation, its impact is transient, since time-averaged differences between the models fall within the large temporal scatter. The context and time-dependent nature of the B-$\rho$ relation underscore the need for comprehensive research and observations to understand the intricate role of magnetic fields in star formation processes across diverse galactic environments.

We analyze the far-ultraviolet spectroscopy of 20 young and massive star clusters (YSCs) in 11 nearby star-forming galaxies. We probe the interstellar gas intervening along the line of sight, detecting several metal absorption lines of a wide range of ionization potentials, from 6.0 eV to 77.5 eV. Multiple-component Voigt fits to the absorption lines are used to study the kinematics of the gas. We find that nearly all targets in the sample feature gas outflowing from 30 up to 190 km per second, often both in the neutral and ionized phase. The outflow velocities correlate with the underlying stellar population properties directly linked to the feedback: the mass of the YSCs, the photon production rate and the instantaneous mechanical luminosity produced by stellar winds and SNe. We detect a neutral inflow in 4 targets, which we interpret as likely not associated with the star cluster but tracing larger scale gas kinematics. A comparison between the outflows energy and that produced by the associated young stellar populations suggests an average coupling efficiency of 10 per cent with a broad scatter. Our results extend the relation found in previous works between galactic outflows and the host galaxy star-formation rate to smaller scales, pointing towards the key role that clustered star formation and feedback play in regulating galaxy growth.

We study a sample of 936 early-type galaxies located in 48 low-z regular galaxy clusters with $M_{200}\geq 10^{14}~ M_\odot$ at $z< 0.1$. We examine variations in the Kormendy relation (KR) according to their location in the projected phase space (PPS) of the clusters. We have used a combination of Bayesian statistical methods to identify possible differences between the fitted relations. Our results indicate that the overall KR is better fitted when we take into account the information about PPS regions. We also find that objects with time since infall $\geq 6.5$ Gyr have a significant statistical difference of the KR coefficients relative to objects that are more recent in the cluster environment. We show that giant central ellipticals are responsible for tilting the KR relation towards smaller slopes. These galaxies present a late growth probably due to cumulative preprocessing during infall, plus cannibalism and accretion of smaller stripped objects near the center of the clusters.

The intermediate neutron capture process (i-process) can develop during proton ingestion events (PIE), potentially during the early stages of low-mass low-metallicity asymptotic giant branch (AGB) stars. We examine the impact of overshoot mixing on the triggering and development of i-process nucleosynthesis in AGB stars of various initial masses and metallicities. We computed AGB stellar models, with initial masses of 1, 2, 3, and 4 M$_{\odot}$ and metallicities in the $-2.5 \le $ [Fe/H] $\le 0$ range, using the stellar evolution code STAREVOL with a network of 1160 nuclei coupled to the transport equations. We considered different overshooting profiles below and above the thermal pulses, and below the convective envelope. The occurrence of PIEs is found to be primarily governed by the amount of overshooting at the top of pulse ($f_{\rm top}$) and to increase with rising $f_{\rm top}$. For $f_{\rm top} =$ 0, 0.02, 0.04, and 0.1, we find that 0 %, 6 %, 24 %, and 86 % of our 21 AGB models with $-2<$ [Fe/H] $<0$ experience a PIE, respectively. We also find that PIEs leave a $^{13}$C-pocket at the bottom of the pulse that can give rise to an additional radiative s-process nucleosynthesis, and ultimately produce a noticeable mixed i+s chemical signature at the surface. Finally, the chemical abundance patterns of 22 observed r/s-stars candidates with $-2<$ [Fe/H] $<-1$ are found to be in reasonable agreement with our AGB model predictions. The binary status of the dwarfs/giants being unclear, we suggest that these stars have acquired their chemical pattern either from the mass transfer of a now-extinct AGB companion or from an early generation AGB star that polluted the natal cloud. Stricter constraints from multi-dimensional hydrodynamical models on overshoot coefficients could deliver new insights into the contribution of AGB stars to heavy elements in the Universe.

We present the results from the NICER observation data of MAXI J1803-298 across the entire 2021 outburst. In the intermediate and soft state, we detect significant absorption lines at $\sim 7.0$ keV and $\sim 6.7$ keV, arising from the X-ray disk wind outflowing with a velocity of hundreds of km per second along our line of sight. The fitting results from photoionized model suggest that the wind is driven by thermal pressure and the mass-loss rate is low. We find a clear transition for iron from predominantly H-like to predominantly He-like during the intermediate-to-soft state transition. Our results indicate this transition for iron is caused by the evolution of the illuminating spectrum and the slow change of the geometric properties of the disk wind together. The coexistence of disk wind and QPOs features in intermediate state is also reported. Our study makes MAXI J1803-298 the first source in which a transition from optical wind to X-ray wind is detected, offering new insights into the evolution of disk winds across an entire outburst and long-term coupling of accretion disks and mass outflows around accreting black holes.

The Local Volume Complete Cluster Survey (LoVoCCS) is an on-going program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts $0.03<z<0.12$ and X-ray luminosities in the 0.1-2.4 keV band $L_{X500c}>10^{44}$ erg/s) with the Dark Energy Camera (DECam), capturing data in $u,g,r,i,z$ bands with a $5\sigma$ point source depth of approximately 25-26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19 deg within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32 deg within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales.

We carry out three-dimensional computations of the accretion rate onto an object (of size $R_{\rm sink}$ and mass $m$) as it moves through a uniform medium at a subsonic speed $v_{\infty}$. The object is treated as a fully-absorbing boundary (e.g. a black hole). In contrast to early conjectures, we show that when $R_{\rm sink}\ll R_{A}=2Gm/v^2$ the accretion rate is independent of $v_{\infty}$ and only depends on the entropy of the ambient medium, its adiabatic index, and $m$. Our numerical simulations are conducted using two different numerical schemes via the Athena++ and Arepo hydrodynamics solvers, which reach nearly identical steady-state solutions. We find that pressure gradients generated by the isentropic compression of the flow near the accretor are sufficient to suspend much of the surrounding gas in a near-hydrostatic equilibrium, just as predicted from the spherical Bondi-Hoyle calculation. Indeed, the accretion rates for steady flow match the Bondi-Hoyle rate, and are indicative of isentropic flow for subsonic motion where no shocks occur. We also find that the accretion drag may be predicted using the Safronov number, $\Theta=R_{A}/R_{\rm sink}$, and is much less than the dynamical friction for sufficiently small accretors ($R_{\rm sink}\ll R_{A}$).

We describe a novel method for modeling the global, steady solar wind using photospheric magnetic fields as a driving boundary condition. Prior wind models in this class include both rapid heuristic methods that use potential field extrapolation and variants thereof, trading rigor for computation speed, and detailed 3D magnetohydrodynamic (MHD) models that attempt to simulate the entire solar corona with a degree of physical rigor, but require large amounts of computation. The Field Line Universal relaXer (FLUX), an open-source numerical code which implements the 'fluxon' semi-lagrangian approach to MHD modeling, provides an intermediate approach between these two general classes. In particular, the fluxon approach to MHD describes the magnetic field through discrete analogues of magnetic field lines, relaxing these structures to a stationary state of force balance. In this work we introduce a one-dimensional solar wind solution along each fieldline, providing an ensemble of solutions that are interpolated back onto a uniform grid at an outer boundary surface. This provides advantages in physical rigor over heuristic semi-analytic techniques, and in computational efficiency over full 3D MHD techniques. Here we describe the underlying methodology and the FLUXPipe modeling pipeline process.

The Imaging X-ray Polarimetry Explorer (IXPE), launched on December 9, 2021, enables X-ray polarimetric observations with unprecedented sensitivity in the 2-8 keV energy range. X-ray polarization allows us to test accretion disk, corona, and emission models of stellar-mass black holes in X-ray binaries found predominantly in soft and hard states of accretion. Every state is dominated by a combination of thermal disk, coronal, or reflected emission$\unicode{x2014}$each type containing different information about the environment around the black hole that can be understood through polarization. In 2022, IXPE measured the polarization signatures of four stellar-mass black holes: Cygnus X-1, 4U 1630-47, Cygnus X-3, and LMC X-1. We report on the physical consequences of these first IXPE observations and the science driven by this new chapter in X-ray polarimetry.

Previous studies on the fitting of spectral energy distributions (SEDs) often apply the external-Compton process to interpret the high-energy peak of low-synchrotron-peaked (LSP) BL Lac objects (LBLs), despite the lack of strong broad emission lines observed for LBLs. In this work, we collect quasi-simultaneous multi-wavelength data of 15 LBLs from the Fermi fourth LAT AGN catalog (4LAC). We propose an analytical method to assess the necessity of external photon fields in the framework of one-zone scenario. Following derived analytical results, we fit the SEDs of these LBLs with the conventional one-zone leptonic model and study their jet physical properties. Our main results can be summarized as follows. (1)We find that most LBLs cannot be fitted by the one-zone synchrotron self-Compton (SSC) model. This indicates that external photons play a crucial role in the high-energy emission of LBLs, therefore we suggest that LBLs are masquerading BL Lacs. (2) We suggest that the $\gamma$-ray emitting regions of LBLs are located outside the broad-line region and within the dusty torus. (3) By extending the analytical method to all types of LSPs in Fermi-4LAC (using historical data), we find that the high-energy peaks of some flat spectrum radio quasars and blazar candidates of unknown types can be attributed to the SSC emission, implying that the importance of external photons could be minor. We suggest that the variability timescale may help distinguish the origin of the high-energy peak.

Weak-line quasars (WLQs) are a notable group of active galactic nuclei (AGNs) that show unusually weak UV lines even though their optical-UV continuum shapes are similar to those of typical quasars. The physical mechanism for WLQs is an unsolved puzzle in the AGN unified model. We explore the properties of UV emission lines by performing extensive photoionization calculations based on Cloudy simulation with different spectral energy distributions (SEDs) of AGNs. The AGN continua are built from several observational empirical correlations, where the black-body emission from the cold disk, the power-law emission from the hot corona, and a soft X-ray excess component are considered. We find that the equivalent width (EW) of C {\footnotesize IV} from our models is systematically lower than observational values if the component of soft X-ray excess is neglected. The EW will increase several times and is roughly consistent with the observations after considering the soft X-ray excess component as constrained from normal type I AGNs. We find that the UV lines are weak for QSOs with quite large BH mass (e.g., $M_{\rm BH}>10^9M_{\odot}$) and weak soft X-ray emission due to the deficit of ionizing photons. As an example, we present the strength of C {\footnotesize IV} based on the multi-band SEDs for three nearby weak-line AGNs, where the weaker soft X-ray emission normally predicts the weaker lines.

We report spectroscopic surveys of planetary nebulae (PNe) in the Milky Way and Andromeda (M31), using the 10.4-m Gran Telescopio Canarias (GTC). The spectra are of high quality and cover the whole optical range, mostly from 3650 \r{A} to beyond 1 micron, enabling detection of nebular emission lines critical for spectral analysis as well as photoionization modeling. We obtained GTC spectra of 24 compact (angular diameter <5 arcsec) PNe located in the Galactic disk, ~3-20 kpc from the Galactic centre, and can be used to constrain stellar evolution models and derive radial abundance gradients of the Milky Way. We have observed 30 PNe in the outer halo of M31 using the GTC. These halo PNe are uniformly metal-rich and probably all evolved from low-mass stars, consistent with the conjecture that they all formed from the metal-rich gas in M31 disk but displaced to their present locations due to galaxy interactions.

The launches of Parker Solar Probe (Parker) and Solar Orbiter (SolO) are enabling a new era of solar wind studies that track the solar wind from its origin at the photosphere, through the corona, to multiple vantage points in the inner heliosphere. A key ingredient for these models is the input photospheric magnetic field map that provides the boundary condition for the coronal portion of many heliospheric models. In this paper, we perform steady-state, data-driven magnetohydrodynamic (MHD) simulations of the solar wind during Carrington rotation 2258 with the GAMERA model. We use the ADAPT and AFT flux transport models and quantitatively assess how well each model matches in-situ measurements from Parker, SolO, and Earth. We find that both models reproduce the magnetic field components at Parker quantitatively well. At SolO and Earth, the magnetic field is reproduced relatively well, though not as well as at Parker, and the density is reproduced extremely poorly. The velocity is overpredicted at Parker, but not at SolO or Earth, hinting that the Wang-Sheeley-Arge (WSA) relation, fine-tuned for Earth, misses the deceleration of the solar wind near the Sun. We conclude that AFT performs quantitatively similarly to ADAPT in all cases and that both models are comparable to a purely WSA heliospheric treatment with no MHD component. Finally, we trace field lines from SolO back to an active region outflow that was observed by Hinode/EIS, and which shows evidence of elevated charge state ratios.

We present an observational study of wind acceleration based on four low-ionization broad absorption line (LoBAL) quasars (J0136, J1238, J1259, J1344). J0136 and J1344 (group-1) are radio quiet and show large BAL-velocity shifts as opposed to stable line-locking associated absorption lines (AALs). Notably, J1344 displays a linear relation between BAL-velocity shift and time interval over three consecutive epochs, characteristic of compelling evidence for BAL acceleration. J1238 and J1259 (group 2) exhibit small BAL-velocity shifts along with steep-spectrum, weak radio emission at 3.0 and 1.4 GHz. All four quasars have spectral energy distributions (SEDs) with a peak at $\lambda_{\rm rest }\sim10~\mu$m, suggesting a link between the BAL acceleration and hot dust emission. The group-2 quasars are redder than group-1 quasars and have a steeper rise at $1<\lambda_{\rm rest }<3~\mu$m in their SEDs. All but J1238 exhibit a steep rise followed by a plateau-like time evolution in BAL-velocity shift. Our investigations, combined with previous studies of BAL acceleration, indicate that (1) the BAL-ISM coupling process is one of the major avenues for the origin of quasar reddening and patchy obscuration, (2) AAL outflows are ubiquitous and likely signify large-scale remnants of BAL winds coupled to interstellar medium (ISM), and (3) wind deceleration that is closely linked to the BAL-ISM coupling process may produce weak radio emission in otherwise radio-quiet quasars.

Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, atmospheric properties, and orbital decay rate. The planet was modeled as a triaxial ellipsoid parameterized by the second-order fluid Love number, $h_2$, which quantifies its radial deformation and provides insight into the interior structure. We measured the tidal deformation of WASP-12b and estimated a Love number of $h_2=1.55_{-0.49}^{+0.45}$ (at 3.2$\sigma$) from its phase curve. We measured occultation depths of $333\pm24$ppm and $493\pm29$ppm in the CHEOPS and TESS bands, respectively, while the dayside emission spectrum indicates that CHEOPS and TESS probe similar pressure levels in the atmosphere at a temperature of 2900K. We also estimated low geometric albedos of $0.086\pm0.017$ and $0.01\pm0.023$ in the CHEOPS and TESS passbands, respectively, suggesting the absence of reflective clouds in the dayside of the WASP-12b. The CHEOPS occultations do not show strong evidence for variability in the dayside atmosphere of the planet. Finally, we refine the orbital decay rate by 12% to a value of -30.23$\pm$0.82 ms/yr. WASP-12b becomes the second exoplanet, after WASP-103b, for which the Love number has been measured (at 3$sigma$) from the effect of tidal deformation in the light curve. However, constraining the core mass fraction of the planet requires measuring $h_2$ with a higher precision. This can be achieved with high signal-to-noise observations with JWST since the phase curve amplitude, and consequently the induced tidal deformation effect, is higher in the infrared.

Fast blue optical transients (FBOTs) represent a new class of highly energetic sources observed from radio to X-rays. High luminosity, light curves and spectra of the sources can be understood if they are associated with supernova-like or tidal disruption events. Radio observations of the transient sources revealed a mildly relativistic expansion of some of the remnants. The high power and mildly relativistic shock velocities are providing favorable conditions for very high energy particle acceleration. In this paper we present a model of particle acceleration in mildly relativistic magnetohydrodynamic (MHD) outflow of the transient source. To construct the non-thermal radiation and cosmic ray spectra in a broad range of energies we combined the microscopic particle-in-cell (PIC) simulations of electron and proton injection at mildly relativistic shock with Monte Carlo technique for high energy particle transport and acceleration. The kinetic PIC simulations provided the energy partition parameter $\epsilon_{e}$ used to fit the observed non-thermal radio emission using the magnetic field amplification mechanisms modelled with Monte Carlo simulations. The model allowed to describe the radio-spectrum of CSS161010 and it's X-ray luminosity. The high X-ray luminosity of AT2018 and AT2020mrf detected during the first weeks can be connected to the jet interaction with the stellar companion in a binary system. The model predicts that FBOTs can accelerate cosmic rays to energies above 10 PeV with a possible upper limit of maximum energy of 100 PeV. With the expected event rate of FBOTs they can contribute to the very high energy cosmic rays population in galaxies.

The study analyzes five models derived from the Pacif parametrization scheme, a generalized form of the Hubble parameter ($H$), yielding various forms of the deceleration parameter (DP) that vary linearly, quadratically, cubically, quartically, and quintically. We aim to explore the impact of these DP variations on late-time evolution and their potential to alleviate cosmological tensions. Additionally, we enhance model constraints by introducing non-diagonal elements into the covariance matrix, thereby simulating correlations between data points and better capturing the true statistical properties of the dataset. Furthermore, our objective is to test the sensitivity of $H_{0}$ and $r_{d}$ to the Pacif parametrization scheme. To avoid imposing a prior from CMB, we treat the sound horizon $r_{d}$ as a free parameter, allowing the late-time data itself to constrain the value of $r_{d}$ along with other cosmological parameters. This is achieved by incorporating the most recent measurements, such as Baryon Acoustic Oscillations (BAO) from the latest galaxy surveys. The considered redshift range spans from (0.106 < z < 2.33), encompassing Hubble measurements obtained through Cosmic Chronometers Methods, Type Ia Supernovae (SNIa), Gamma-Ray Bursts (GRBs), and Quasars (Q). Additionally, an extra constraint is introduced by integrating the latest Hubble constant measurement from Riess in 2022. Our analysis yields optimal fit values for the Hubble parameter $H_{0}$ and sound horizon $r_{d}$. The results highlight a notable consistency between the values derived from the Pacif parametrization models and those obtained from Planck CMB data. Finally, we employ the Akaike information criterion approach to compare the three models, revealing that none of them can be discounted based on the latest observational measurements.

The strengthening of tensions in the cosmological parameters has led to a reconsideration of fundamental aspects of standard cosmology. The tension in the Hubble constant can also be viewed as a tension between local and early Universe constraints on the absolute magnitude $M_B$ of Type Ia supernova. In this work, we reconsider the possibility of a variation of this parameter in a model-independent way. We employ neural networks to agnostically constrain the value of the absolute magnitude as well as assess the impact and statistical significance of a variation in $M_B$ with redshift from the Pantheon+ compilation, together with a thorough analysis of the neural network architecture. We find an indication for a transition redshift at the $z\approx 1$ region.

The time evolution of the dependence of the mass accretion rate with the stellar mass and the disk mass represents a fundamental way to understand the evolution of protoplanetary disks and the formation of planets. In this work, we present observations with X-Shooter of 26 Class II very low-mass stars and brown dwarfs in the Ophiuchus, Cha-I, and Upper Scorpius star-forming regions (SFRs). These new observations extend down to SpT M9 ($\sim$0.02 $M_\odot$) the measurement of the mass accretion rate in Ophiuchus and Cha-I and add 11 very-low-mass stars to the sample of objects studied with broadband spectroscopy in Upper Scorpius. We obtained their SpT, extinction and physical parameters, and we used the intensity of various emission lines to derive their accretion luminosity and mass accretion rates. Combining these new observations with data from the literature, we compare relations between accretion and stellar and disk properties of four different SFRs with different ages: Ophiuchus (1 Myr), Lupus (2 Myr), Cha-I (3 Myr), and Upper Scorpius (5-12 Myr). We find the slopes of the $L_*-L\mathrm{_{acc}}$ and $M_*-\dot{M}\mathrm{_{acc}}$ relationships to steepen between Ophiuchus, Lupus, and Cha-I and that both relationships may be better described with a single power law. We also find the relationship between the disk mass and the mass accretion rate of the stellar population to steepen with time down to the age of Upper Scorpius. Overall, we observe hints of a faster evolution into low accretion rates of low-mass stars and brown dwarfs. We also find that brown dwarfs present higher $M\mathrm{_{disk}}/\dot{M}\mathrm{_{acc}}$ ratios (i.e., longer accretion depletion timescales) than stars in Ophiuchus, Lupus, and Cha-I. This apparently contradictory result may imply that the evolution of protoplanetary disks around brown dwarfs is different from what is seen in the stellar regime.

Similar to the cases of anemone jets, two-sided loop solar jets could also be produced by either flux emergence from the solar interior or small scale filament eruptions. Using the high-quality data from the Solar Dynamic Observatory (SDO), we analyzed a two-sided loop solar jet triggered by the eruption of a small filament in this paper. The jet was occurred in a pre-existing big filament channel. The detailed processes involved in the small filament eruption, the interaction between the erupted filament and the big filament channel, and the launch of the two-sided loop jet are presented. The observations further revealed notable asymmetry between the two branches of the jet spire, with the northeastern branch is narrow and short, while the southern branch is wide and long and accompanied by discernible untwisting motions. We explored the unique appearance of the jet by employing the local potential field extrapolation to calculate the coronal magnetic field configuration around the jet. The photospheric magnetic flux below the small filament underwent cancellation for approximately 7 hours before the filament eruption, and the negative flux near the southern foot-point of the filament decreased by about 56 percent during this interval. Therefore, we proposed that the primary photospheric driver of the filament eruption and the associated two-sided loop jet in this event is flux cancellation rather than flux emergence.

Out of the more than 5,000 detected exoplanets a considerable number belongs to a category called 'mini-Neptunes'. Interior models of these planets suggest that they have some primordial, H-He dominated atmosphere. As this type of planet does not occur in the solar system, understanding their formation is a key challenge in planet formation theory. Unfortunately, quantifying the H-He, based on their observed mass and radius, is impossible due to the degeneracy of interior models. We explore the effects that different assumptions on planet formation have on the nebular gas accretion rate, particularly by exploring the way in which solid material interacts with the envelope. This allows us to estimate the range of possible post-formation primordial envelopes. Thereby we demonstrate the importance of envelope enrichment on the initial primordial envelope which can be used in evolution models. We apply formation models that include different solid accretion rate prescriptions. Our assumption is that mini-Neptunes form beyond the ice-line and migrate inward after formation, thus we form planets in-situ at 3 and 5 au. We consider that the envelope can be enriched by the accreted solids in the form of water. We study how different assumptions and parameters influence the ratio between the planet's total mass and the fraction of primordial gas. The primordial envelope fractions for small- and intermediate-mass planets (total mass below 15 M$_{\oplus}$) can range from 0.1% to 50%. Envelope enrichment can lead to higher primordial mass fractions. We find that the solid accretion rate timescale has the largest influence on the primordial envelope size. Primordial gas accretion rates can span many orders of magnitude. Planet formation models need to use a self-consistent gas accretion prescription.

Aims: This work aims to identify the mechanism driving pulsations in hard X-ray (HXR) and microwave emission during solar flares. Here, by using combined HXR and microwave observations from Solar Orbiter/STIX and EOVSA we investigate an X1.3 GOES class flare, 2022-03-30T17:21:00, which displays pulsations on timescales evolving from ~ 7 s in the impulsive phase to ~ 35 s later in the flare. Methods: The temporal, spatial and spectral evolution of the HXR and microwave pulsations during the impulsive phase of the flare are analysed. Images are reconstructed for individual peaks in the impulsive phase and spectral fitting is performed at high cadence throughout the first phase of pulsations. Results: Imaging analysis demonstrates that the HXR and microwave emission originates from multiple sites along the flare ribbons. The brightest sources and the location of the emission changes in time. Through HXR spectral analysis, the electron spectral index is found to be anti-correlated with the HXR flux showing a "soft-hard-soft" spectral index evolution for each pulsation. The timing of the associated filament eruption coincides with the early impulsive phase. Conclusions: Our results indicate that periodic acceleration and/or injection of electrons from multiple sites along the flare arcade is responsible for the pulsations observed in HXR and microwave. The evolution of pulsation timescales is likely a result of changes in the 3D magnetic field configuration in time related to the associated filament eruption.

Quasar winds can shock and sweep up ambient interstellar medium (ISM) gas, contributing to galactic quenching. We combine and extend past models of energy-conserving shock bubbles around quasars, investigate model implications from an observational standpoint, and test model predictions using new high-resolution spectroscopic observations of the broad absorption line quasar SDSS J030000.56+004828.0 (J0300). Even with constant energy input from the wind, a bubble's expansion decelerates over time as more ISM gas is swept up. Our new observations enable a direct search for this deceleration. We obtain the tightest reported 3-sigma limit on the average rest-frame deceleration (or acceleration) of a quasar outflow: |a|$<$0.1 km s$^{-1}$ yr$^{-1}$ ($<3 \times 10^{-4}$ cm s$^{-2}$) in the relatively low-velocity Ca II outflow of J0300 over 9.65 rest-frame years. We can satisfy these limits with certain parameter choices in our model, but the large velocity range of the Ca II absorption in J0300 rules out the hypothesis that such gas shares the velocity of the swept-up ISM gas in a self-similar shock bubble. We investigate the possibility of ram-pressure acceleration of preexisting ISM clouds and conclude that the velocity range seen in Ca II in J0300 is potentially consistent with such an explanation. The Ca II-absorbing gas clouds in J0300 have been inferred to have high densities by Choi et al., in which case they can only have been accelerated to their current speeds if they were originally at least an order of magnitude less dense than they are today.

Recently, a lack of supernova-associated long-duration gamma-ray burst (GRB 230307A) at such a low redshift $z=0.065$, but associated with a possible kilonova emission, has attracted great attention. Its heavy element nucleosynthesis and the characteristic of soft X-ray emission suggests that the central engine of GRB 230307A is magnetar which is originated from a binary compact star merger. The calculated lower value of $\varepsilon \sim 0.05$ suggests that the GRB 230307A seems to be with ambiguous progenitor. The lower value of $f_{\rm eff}=1.23$ implies that the GRB 230307A is not likely to be from the effect of "tip of iceberg". We adopt the magnetar central engine model to fit the observed soft X-ray emission with a varying efficiency and find that the parameters constraints of magnetar fall into a reasonable range, i.e., $B<9.4\times10^{15}$ G and $P<2.5$ ms for $\Gamma_{\rm sat} = 10^3$, and $B<3.6\times10^{15}$ G and $P<1.05$ ms for $\Gamma_{\rm sat} = 10^4$. Whether the progenitor of GBR 230307A is from the mergers of neutron star - white dwarf (NS - WD) or neutron star - neutron star (NS - NS) remains unknown.

Accretion disks surrounding stellar mass black holes (BHs) have been suggested as potential locations for the nucleosynthesis of light elements, which are our primary observational discriminant of multiple stellar populations within globular clusters. The population of enriched stars in globular clusters are enhanced in N14, Na23, and sometimes in Al27 and/or in K39. In this study, our aim is to investigate the feasibility of initiating nucleosynthesis for these four elements in BH accretion disks, considering various internal parameters such as the temperature of the gas and timescale of the accretion. To achieve this, we employed a 132-species reaction network. We used the slim disk model, suitable for the Super-Eddington mass accretion rate and for geometrically and optically thick disks. We explored the conditions related to the mass, mass accretion rate, viscosity, and radius of the BH-accretion disk system that would allow for the creation of N14, Na23, Al27, and K39 before the gas is accreted onto the central object. Our findings reveal that there is no region in the parameter space where the formation of Na23 can occur and only a very limited region where the formation of N14, Al27, and K39 is plausible. Specifically, this occurs for BHs with masses lower than 10 solar masses, with a preference toward even lower mass values and extremely low viscosity parameters ($\alpha <10^{-3}$). Such values are highly unlikely based on current observations of stellar mass BHs. However, such low mass BHs could actually exist in the early universe, as so-called primordial BHs. In conclusion, our study suggests that the nucleosynthesis within BH accretion disks of four elements of interest for the multiple stellar populations is improbable, but not impossible, using the slim disk model.

Solar flare ribbons are intense brightenings of principally chromospheric material that are responsible for a large fraction of the chromospheric emission in solar and stellar flares. We present an on-disc observation of flare ribbon substructures in an X9.3-class flare observed by the Swedish 1-m Solar Telescope. We identify categories of ribbon substructures seen in the Ca II 8542 \AA, H$\alpha$, and Ca II K lines, focusing on their spatial locations and their (spectro-)polarimetric properties. Color Collapsed Plotting (COCOPLOT) software is used to assist in identifying areas of interest. We present five categories of spectral profiles within the general body of the flare ribbon: (1) Extremely broadened spectral line profiles, where the standard Fabry-Perot interferometer wavelength windows ($\approx 70$ km s$^{-1}$) are insufficiently wide to allow for a complete analysis of the dynamics and atmospheric conditions. The mechanisms causing this degree of this broadening are not yet clearly understood. (2) Long-lived, dense kernels that manifest as more saturated chromospheric line profiles with lower signal in both Stokes parameters. They are interpreted as footpoints of bunched magnetic field loops, whose chromospheric lines form at greater heights than the nearby areas. (3) Doppler-shifted leading edges of the flare ribbon in regions that transiently display lower Stokes signals due to the emission dominating at greater heights in the atmosphere. (4) Condensed coronal rain overlapping the flare ribbons in the line of sight, producing exceptionally high Doppler shifts near the footpoints...

Very recently, a particularly long gamma-ray burst (GRB) 230307A was reported and proposed to originate from a compact binary merger based on its host galaxy property, kilonova, and heavy elements. More intriguingly, a very early plateau followed by a rapid decline in soft X-ray band was detected in its light curve by the Lobster Eye Imager for Astronomy, indicating strong evidence of the existence of a magnetar as the merger product. This work explores that the Magnetar Wind Internal Gradual MAgnetic Dissipation (MIGMAD) model, in which the radiative efficiency evolves over time, successfully fits it to the observed data. Our results reinforce the notion that the X-ray plateau serves as a powerful indicator of a magnetar and imply that an evolving efficiency is likely to be a common feature in X-ray plateaus of GRB afterglows. In addition, we also discuss the explanations for the prompt emission, GRB afterglows, as well as kilonova, and predict possible kilonova afterglows in a magnetar central engine.

Here I present a somewhat personal review of the results of my research carried out in the last couple of decades, mostly focused on the fields of stellar astrophysics, the interstellar medium, and the high energy sources. They have been obtained mainly from dedicated interferometric observations in the radio band at low frequencies of the electromagnetic spectrum.

The classical theory of differentiation states that due to the heat generated by the decay of radioactive elements, some asteroids form an iron core, an olivine-rich mantle, and a crust. The collisional breakup of these differentiated bodies is expected to lead to exposed mantle fragments, creating families of newly-formed asteroids. Among these new objects, some are expected to show an olivine-rich composition in spectroscopic observations. However, several years of spectrophotometric surveys have led to the conclusion that olivine-rich asteroids are rare in the asteroid main belt, and no significant concentration of olivine-rich bodies in any asteroid family has been detected to date. Using ESA's Gaia DR3 reflectance spectra, we show that the family (36256) 1999 XT17 presents a prominence of objects that are likely to present an olivine-rich composition (A-type spectroscopic class). If S-complex asteroids as the second most prominent spectroscopic class in the family are real family members, then arguably the 1999 XT17 family has originated from the break-up of a partially differentiated parent body. Alternatively, if the S-complex asteroids are interlopers, then the 1999 XT17 family could have originated from the breakup of an olivine-rich body. This body could have been part of the mantle of a differentiated planetesimal, which may have broken up in a different region of the Solar System, and one of its fragments (i.e. the parent body of the 1999 XT17 family) could have been dynamically implanted in the main belt.

We simulate atmospheric fractionation in escaping planetary atmospheres using IsoFATE, a new open-source numerical model. We expand the parameter space studied previously to planets with tenuous atmospheres that exhibit the greatest helium and deuterium enhancement. We simulate the effects of EUV-driven photoevaporation and core-powered mass loss on deuterium-hydrogen and helium-hydrogen fractionation of sub-Neptune atmospheres around G, K, and M stars. Our simulations predict prominent populations of deuterium- and helium-enhanced planets along the upper edge of the radius valley with mean equilibrium temperatures of 370 K and as low as 150 K across stellar types. We find that fractionation is mechanism-dependent, so constraining He/H and D/H abundances in sub-Neptune atmospheres offers a unique strategy to investigate the origin of the radius valley around low-mass stars. Fractionation is also strongly dependent on retained atmospheric mass, offering a proxy for planetary surface pressure as well as a way to distinguish between desiccated enveloped terrestrials and water worlds. Deuterium-enhanced planets tend to be helium-dominated and CH4-depleted, providing a promising strategy to observe HDO in the 3.7 um window. We present a list of promising targets for observational follow-up.

The coalescence of the most massive black hole (MBH) binaries releases gravitational waves (GWs) within the detectable frequency range of Pulsar Timing Arrays (PTAs) $(10^{-9} - 10^{-6})$ Hz. The incoherent superposition of GWs from MBH mergers, the stochastic Gravitational Wave Background (GWB), can provide unique information on MBH parameters and the large-scale structure of the Universe. The recent evidence for a GWB reported by the PTAs opens an exciting new window onto MBHs and their host galaxies. However, the astrophysical interpretation of the GWB requires accurate estimations of MBH merger timescales for a statistically representative sample of galaxy mergers. This is numerically challenging; a high numerical resolution is required to avoid spurious relaxation and stochastic effects whilst a large number of simulations is needed to sample a cosmologically representative volume. Here, we present a new multi-mass modelling method to increase the central resolution of a galaxy model at a fixed particle number. We follow mergers of galaxies hosting central MBHs with the Fast Multiple Method code Griffin at two reference resolutions and with two refinement schemes. We show that both refinement schemes are effective at increasing central resolution, reducing spurious relaxation and stochastic effects. A particle number of $N\geq 10^{6}$ within a radius of 5 times the sphere of influence of the MBHs is required to reduce numerical scatter in the binary eccentricity and the coalescence timescale to <30$\%$; a resolution that can only be reached at present with the mass refinement scheme.

Hypervelocity stars (HVSs) unbound to the Galaxy can be formed with extreme stellar interactions. Observational evidence comes from measurements of radial velocities (RVs) of objects crossing the Galactic halo and of tangential velocities based on high proper motions (HPMs) and distances of relatively nearby stars. I searched for new HVS candidates and reviewed known objects using their Gaia astrometric measurements. Candidates were selected with significant Gaia parallaxes of >0.1mas, proper motions of >20mas/yr, and computed Galactocentric tangential velocities $vtan\_g$>500km/s. The DR2 and DR3 samples of several thousand HVS candidates were studied with respect to their proper motions, sky distribution, number of observations, location in crowded fields, colour-magnitude diagrams, selection effects with magnitude, and RVs in DR3. The most extreme ($vtan\_g$>700km/s) and nearest (within 4kpc) 72 DR3 HVS candidates were investigated with respect to detected close neighbours, flags and astrometric quality parameters of objects of similar magnitudes in DR3. The quality checks involved HPM objects in a global comparison and all objects in the vicinity of each target. Spurious HPMs in the Galactic centre region led to false HVS interpretations in Gaia DR2 and are still present in DR3, although to a lesser extent. Otherwise there is good agreement between the HPMs of HVS candidates in DR2 and DR3. However, HVS candidates selected from DR2 tend to have larger parallaxes hence lower tangential velocities in DR3. Most DR3 RVs are much lower than the tangential velocities, indicating that the DR3 HVS candidates are still affected by underestimated parallaxes. None of the 72 extreme nearby DR3 HVS candidates, including three D$^6$ stars, passed all quality checks. Their tangential velocities may turn out to be smaller, but at least some of them still appear unbound to the Galaxy. (abbreviated)

Narrow-line Seyfert 1 (NLS1) galaxies are a peculiar subclass of active galactic nuclei (AGN). Among them, TXS 1206+549 belongs to a small group of radio-loud and gamma-ray-emitting NLS1 galaxies. We focus on the radio properties of this galaxy by analysing archival, high-resolution, very long baseline interferometry (VLBI) imaging observations taken at 8 GHz frequency in six epochs between 1994 and 2018. Using the milliarcsecond-scale radio structure, we can resolve a core and a jet component whose angular separation increases by (0.055 +/- 0.006) mas/yr. This corresponds to an apparent superluminal jet component motion of (3.5 +/- 0.4) c. From the core brightness temperature and the jet component proper motion, we determine the characteristic Doppler-boosting factor, the bulk Lorentz factor, and the jet viewing angle. We find no compelling evidence for a very closely aligned blazar-type jet. The parameters for TXS 1206+549 resemble those of radio-loud quasar jets with a moderate Lorentz factor (Gamma~4) and theta~24 deg inclination to the line of sight.

Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionation, this provides the possible opportunity to probe the ability of isotopes ratios to trace the ages of high star forming regions, such as SSCs. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC253 are correlated with their different ages as derived from stellar population modelling. We have measured abundance ratios of CO, HCN and HCO$^+$ isotopologues in six regions containing SSCs within NGC253 using high spatial resolution (1.6",$\sim 28$pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We have then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star formation tracers (such as Br$\gamma$) to starburst stellar population evolution models. We do not find any significant trend with age for the CO and HCN isotopologue ratios on the timescales for the ages of the SSC* regions observed. The driving factors of these ratios within SSCs could be the Initial Mass Function as well as possibly fractionation effects. To further probe these effects in SSCs over time a larger sample of SSCs must be observed spanning a larger age range.

The prominent radio quasar PKS 2215+020 (J2217+0220) was once labelled as a new laboratory for core--jet physics at redshift z=3.572 because of its exceptionally extended jet structure traceable with very long baseline interferometric (VLBI) observations up to a ~600 pc projected distance from the compact core and a hint of an arcsec-scale radio and an X-ray jet. While the presence of an X-ray jet could not be confirmed later, this active galactic nucleus is still unique at high redshift with its long VLBI jet. Here, we analyse archival multi-epoch VLBI imaging data at five frequency bands from 1.7 to 15.4 GHz covering a period of more than 25 years from 1995 to 2020. We constrain apparent proper motions of jet components in PKS 2215+020 for the first time. Brightness distribution modeling at 8 GHz reveals a nearly 0.02 mas/yr proper motion (moderately superluminal with apparently two times the speed of light), and provides delta=11.5 for the Doppler-boosting factor in the inner relativistic jet that is inclined within 2 deg to the line of sight and has a Gamma=6 bulk Lorentz factor. These values qualify PKS 2215+020 as a blazar, with rather typical jet properties in a small sample of only about 20 objects at z>3.5 that have similar measurements to date. According to the 2-GHz VLBI data, the diffuse and extended outer emission feature at ~60 mas from the core, probably a place where the jet interacts with and decelerated by the ambient galactic medium, is consistent with being stationary, albeit slow motion cannot be excluded based on the presently available data.

Submillimeter single-dish telescopes offer two key advantages compared to interferometers: they can efficiently map larger portions of the sky and recover larger spatial scales. Nonetheless, fluctuations in the atmosphere limit the accurate retrieval of signals from astronomical sources. Therefore, we introduce a user-friendly simulator named ${\tt maria}$ to optimize scanning strategies and instrument designs to efficiently reduce atmospheric noise and filtering effects. We further use this tool to produce synthetic time streams and maps from hydrodynamical simulations, enabling a fair comparison between theory and reality. ${\tt maria}$ has implemented a suite of telescope and instrument designs intended to mimic current and future facilities. To generate synthetic time-ordered data, each mock observatory scans through the atmosphere in a configurable pattern over the celestial object. We generate evolving and location-and-time-specific weather for each of the fiducial sites using a combination of satellite and ground-based measurements. While ${\tt maria}$ is a generic virtual telescope, this study specifically focuses on mimicking broadband bolometers observing at 100 GHz. To validate our virtual telescope, we compare the mock time streams with real MUSTANG-2 observations and find that they are quantitatively similar by conducting a k-sample Anderson-Darling test resulting in p<0.001. Subsequently, we image the time-ordered data to create noise maps and mock observations of clusters of galaxies for both MUSTANG-2 and an instrument concept for the 50m Atacama Large Aperture Submillimeter Telescope (AtLAST). Furthermore, using ${\tt maria}$, we find that a 50m dish provides the highest levels of correlation of atmospheric signals across adjacent detectors compared to smaller apertures (e.g., 42-cm and 6-m survey experiments), facilitating removal of atmospheric signal on large scales.

This letter presents, for the first time, direct constraints on the black-hole-to-halo-mass relation using weak gravitational lensing measurements. We construct type I and type II Active Galactic Nuclei (AGNs) samples from the Sloan Digital Sky Survey (SDSS), with a mean redshift of 0.4 0.1 for type I (type II) AGNs. This sample is cross-correlated with weak lensing shear from the Ultraviolet Near Infrared Northern Survey (UNIONS). We compute the excess surface mass density of the halos associated with $36,181$ AGNs from $94,308,561$ lensed galaxies and fit the halo mass in bins of black-hole mass. We find that more massive AGNs reside in more massive halos. We see no evidence of dependence on AGN type or redshift in the black-hole-to-halo-mass relationship when systematic errors in the measured black-hole masses are included. Our results are consistent with previous measurements for non-AGN galaxies. At a fixed black-hole mass, our weak-lensing halo masses are consistent with galaxy rotation curves, but significantly lower than galaxy clustering measurements. Finally, our results are broadly consistent with state-of-the-art hydro-dynamical cosmological simulations, providing a new constraint for black-hole masses in simulations.

Supernova remnants (SNRs) are profoundly affected by their ambient medium. We present carbon monoxide (CO) observations around two mixed morphology SNRs, VRO 42.05.01 and G 350.0-2.0, that look remarkably similar in continuum radio emission, showing what we refer to as a shell and wing shape. It has been proposed that the shell and wing shape is the result of environmental effects, in the form of a sharp density gradient or discontinuity. Our motivation for studying these two sources jointly is that if the dense molecular environment causes the development of these sources' shell and wing shape, then these two sources' environments must be similar. This is contrary to what we observe. In the case of VRO 42.05.01, we have found direct evidence of an interaction with its molecular environment, in the form of broadened $^{12}$CO line profiles, high $(J=2-1)$ to $(J=1-0)$ line ratios, and arc features in position-velocity space. We interpret some of these features to be associated with the SNR shock, and some of them to be due to the presence of a pre-supernova stellar wind. We have found no such features in the abundant molecular gas surrounding G 350.0-2.0. We have also made a spectral index map of G 350.0-2.0, and we see that the radio spectrum of G 350.0-2.0 steepens significantly at frequencies $<200$~MHz, much like that of VRO 42.05.01. In spite of their spectral and morphological similarities, these two sources look substantially different in their optical and infrared emission. The lack of large-scale correspondence between the radio continuum and the molecular material, as well as the differences in the excitation and morphological properties of the molecular gas surrounding both sources, lead us to conclude that the shell and wing morphology of these two sources is not due to interactions with a similar ambient molecular ISM.

We present a deep Chandra observation of the shell supernova remnant G32.4+0.1, whose featureless X-ray spectrum has led to its classification as an X-ray synchrotron-dominated supernova remnant (SNR). We find a partial shell morphology whose outline is quite circular, with a radius of about 11 pc at an assumed distance of 11 kpc. Thermal and power-law spectral models for three relatively bright regions provided equally good fits, but the absence of spectral lines required ionization timescales from thermal fits that are inconsistent with mean densities derived from emission measures. We thus confirm the nonthermal, i.e., synchrotron, origin of X-rays from G32.4+0.1. Shock velocities needed to accelerate electrons to the required TeV energies are >~1000 km/s, giving remnant ages <~5,000 -- 9,000 yr. There is no obvious X-ray counterpart to the radio pulsar PSR J1850--0026, but its position adjoins a region of X-ray emission whose spectrum is somewhat harder than that of other regions of the shell, and which may be a pulsar-wind nebula (PWN), though its spectrum is steeper than almost all known X-ray PWNe. The distance of the pulsar from the center of symmetry of the shell disfavors a birth in a supernova event at that location only a few thousand years before: either the pulsar (and putative PWN) are not associated with the shell SNR, requiring a coincidence of both position and (roughly) absorbing column density, or the SNR is much older, making the origin of nonthermal emission problematic.

We present the detailed spectral energy distribution (SED) modeling of 14 local (ultra)luminous infrared galaxies (U/LIRGs) with outstanding photometric data from the literature covering ultraviolet-far-infrared and radio bands ($\sim$50\,MHz to $\sim$30\,GHz). We employ the CIGALE SED fitting code for ultraviolet--far-infrared--radio SED modeling. For radio-only SED modeling, we use the UltraNest package, leveraging its nested sampling algorithm. Combining results from our previous study on 11 LIRGs, we discuss the global astrophysical properties of a sample of 25 starburst galaxies ($z<$0.5). Their radio spectra are frequently characterized by bends and turnovers, with no indication of ULIRGs exhibiting more complicated SEDs than LIRGs despite showing more signs of interactions. Including radio measurements in CIGALE modeling constrained the dust luminosity and star formation rate (SFR) estimates by more than one order of magnitude better than previously reported for starburst galaxies. We show that total and nonthermal radio luminosity at 1.4\, and 4.8\,GHz frequencies can be good estimators of the recent SFR for all LIRGs and those ULIRGS with insignificant AGN influence. A weaker but still significant correlation is observed between radio SFR at 1.4\,GHz and old (averaged over 100\,Myr) SFR based on SED modeling, indicative of multiple episodes of starburst activity during their lifetime. The thermal radio luminosity at 4.8\,GHz is a better tracer of recent star formation than at 1.4\,GHz. Our modeled nonthermal radio spectral indices are not statistically significantly correlated with redshift, stellar mass, SFR, specific SFR, and dust mass.

We measure the kinematic signature arising from the Milky Way (MW) disc moving with respect to the outer stellar halo, which is observed as a dipole signal in the kinematics of stellar halo tracers. We quantify how the reflex motion varies as a function of Galactocentric distance, finding that (i) the amplitude of the dipole signal increases as a function of radius, and (ii) the direction moves across the sky. We compare the reflex motion signal against a compilation of published models that follow the MW--LMC interaction. These models show a similar trend of increasing amplitude of the reflex motion as a function of distance, but they do not reproduce the direction of the disc motion with respect to the stellar halo well. We also report mean motions for the stellar halo as a function of distance, finding radial compression in the outer halo and nonzero prograde rotation at all radii. The observed compression signal is also present in MW--LMC models, but the rotation is not, which suggests that the latter is not induced by the LMC. We extensively validate our technique to measure reflex motion against idealised tests. We discuss prospects for directly constraining the mass and orbital history of the LMC through the impact on the motion of the MW stellar disc, and how the modelling of the reflex motion can be improved as more and better data become available.

Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, characterising their ionisation rate is still challenging from an observational point of view. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionization rate $\zeta_2$ from molecular line observations. These have been derived from the chemical kinetics of the involved species, but they have not been validated using synthetic data processed with a standard observative pipeline. We aim to bridge this gap. We perform the radiative transfer on a set of three-dimensional magneto-hydrodynamical simulations of prestellar cores, exploring different initial $\zeta_2$, evolutionary stages, types of radiative transfer (e.g. assuming local-thermodynamic-equilibrium conditions), and telescope responses. We then compute the column densities of the involved tracers to determine $\zeta_2$, using, in particular, the equation proposed by Bovino et. al (2020) and by Caselli et al. (1998) both used nowadays. Our results confirm that the method of Bovino et al. (2020) accurately retrieves the actual $\zeta_2$ within a factor of $2-3$, in the physical conditions explored in our tests. Since we also explore a non-local thermodynamic equilibrium radiative transfer, this work indirectly offers insights into the excitation temperatures of common transitions at moderate volume densities ($n\approx 10^5 \, \rm cm^{-3}$). We have also performed a few tests using the formula proposed by Caselli et al. (1998), which overestimates the actual $\zeta_2$ by at least two orders of magnitudes. We also consider a new derivation of this method, which, however, still leads to large overestimates.

Certain relations among neutron-star observables that are insensitive to the equation of state are known to exist. Such universal relations have been shown to be valid for cold and stationary neutron stars. Here, we study these relations in more dynamic scenarios: protoneutron stars and hypermassive neutron stars. First, we study protoneutron stars. We use an effective equation of state, extracted from three-dimensional core-collapse supernova simulations, to obtain the structure of spherically symmetric protoneutron stars. We then consider nonradial oscillations to compute their $f$-mode frequency ($f$), as well as slow rotation and small tidal deformation, to compute their moment of inertia ($I$), spin-induced quadrupole moment ($Q$), and Love number. We find that well-established universal relations for cold neutron stars involving these observables ($I$-Love-$Q$ and $f$-Love relations) are approximately valid for protoneutron stars, with a deviation below $\approx$ 10$\%$ for a postbounce time above $\approx$ 0.5 s, considering eight different supernova progenitors and the SFHo equation of state. Next, we study hypermassive neutron stars. We obtain a new universal relation between the $f$-mode frequency and the compactness of cold and nonrotating neutron stars, using bulk quantities. We show that this relation has an equation-of-state-variation of $\approx$ $3\%$, considering a set of ten equations of state. Using results from binary neutron star merger simulations, we study the evolution of hypermassive neutron stars on the $f$-$C$ plane, considering two different mass ratios and the SFHo equation of state. We find that the relation between the peak frequency of the gravitational-wave signal and the compactness from these hypermassive neutron stars deviates from the universal $f$-$C$ relation by 70 $-$ 80$\%$, when the peak frequency is taken directly as a proxy for the $f$-mode.

This study analyzes secular dynamics using averaged equations that detail tidal effects on the motion of two extended bodies in Keplerian orbits. It introduces formulas for energy dissipation within each body of a binary system. The equations, particularly in contexts like the Sun-Mercury system, can be delineated into a fast-slow system. A significant contribution of this work is the demonstration of the crucial role complex rheological models play in the capture by spin-orbit resonances. This is particularly evident in the notable enlargement of the basin of attraction for Mercury's current state when transitioning from a single characteristic time rheology to a dual characteristic time model, under the constraint that both models comply with the same estimate of the complex Love number at orbital frequency. The study also underscores the importance of Mercury's elastic rigidity on secular timescales.

In this study, we formulate a set of differential equations for a binary system to describe the secular-tidal evolution of orbital elements, rotational dynamics, and deformation (flattening), under the assumption that one body remains spherical while the other is slightly aspherical throughout the analysis. By applying singular perturbation theory, we analyze the dynamics of both the original and secular equations. Our findings indicate that the secular equations serve as a robust approximation for the entire system, often representing a slow-fast dynamical system. Additionally, we explore the geometric aspects of spin-orbit resonance capture, interpreting it as a manifestation of relaxation oscillations within singularly perturbed systems.

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, $N$-body simulations, such as the public \textsc{AbacusPNG} set presented in this study, consisting of 9 boxes, each of size $L_{\rm box} = 2~{\rm Gpc}/h$ and particle mass of $1.01 \times 10^{10} \ M_\odot/h$, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory and the `separate universe' method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, $b_\phi$. As a science application, we investigate the link between halo assembly bias and $b_\phi$ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between $b_\phi$ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the $b_\phi$ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasi-stellar object (QSO) samples and comment on the effect on $f_{\rm NL}$ constraints for the allowed galaxy-halo models (note that $\sigma [f_{\rm NL}] \propto \frac{\sigma [b_\phi]}{b_\phi}$). We find that the error on $f_{\rm NL}$ is 21, 6, 22 for the LRGs at $z = 0.5$ and $z = 0.8$ and QSOs at $z = 1.4$, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

We study lepto-axiogenesis in theories where the right-handed neutrino is light enough that its dynamics affect the determination of the baryon asymmetry. When compared with theories of high-scale lepto-axiogenesis where the Majorana neutrino mass may be treated as an effective dimension-five operator, we find that the predicted saxion mass is lower. Two distinct scenarios emerge. In the first, processes that generate the baryon asymmetry are in equilibrium down to the mass of the right-handed neutrino. In the second, the relevant processes never reach equilibrium; the baryon number freezes in. We comment on implications for supersymmetric spectra and discuss constraints on late decays of supersymmetric relics and from dark radiation. In contrast to high-scale lepto-axiogenesis, which predicts superpartners with masses of 10-100 TeV or more, we find this scenario is consistent with a wider range of superpartner masses, all the way down to current direct search bounds.

The multi-staged XENON program at INFN Laboratori Nazionali del Gran Sasso aims to detect dark matter with two-phase liquid xenon time projection chambers of increasing size and sensitivity. The XENONnT experiment is the latest detector in the program, planned to be an upgrade of its predecessor XENON1T. It features an active target of 5.9 tonnes of cryogenic liquid xenon (8.5 tonnes total mass in cryostat). The experiment is expected to extend the sensitivity to WIMP dark matter by more than an order of magnitude compared to XENON1T, thanks to the larger active mass and the significantly reduced background, improved by novel systems such as a radon removal plant and a neutron veto. This article describes the XENONnT experiment and its sub-systems in detail and reports on the detector performance during the first science run.

We study the post-Newtonian limit of higher-order scalar-tensor theories that are degenerate in the unitary gauge. They can be conveniently described by the effective field theory (EFT) of dark energy. We determine all the parametrized post-Newtonian (PPN) parameters in terms of the EFT of dark energy parameters. Experimental bounds on the PPN parameters are then translated to constraints on the EFT parameters. We present a Lagrangian of a unitary degenerate higher-order scalar-tensor theory characterized by a single function of the kinetic term of the scalar field whose PPN parameters have the same values as in general relativity.

We developed a mid-infrared optical frequency comb-based Fourier-transform spectrometer and performed a line-shape study of the fundamental vibrational band of CO perturbed by N${_2}$, which is crucial for atmospheric science and astronomical observations. The comb-based FTS enabled us to measure the whole vibrational band with high resolution and precision at several pressures between 10 and 400 Torr. Observed absorption profiles were fitted with the speed-dependent Voigt profile. Collisional broadening, speed-dependent collisional width and shift coefficients are derived. The reliability of our results is established from considerations of systematic errors and comparison with previous studies.

We examine the main properties of gravitational waves (GWs) emitted by transient hyperbolic encounters of black holes. We begin by building the set of basic variables most relevant to setting our problem. After exposing the ranges of masses and eccentricities accessible at a given GW frequency, we analyze the dependence of the gravitational strain on those parameters and determine the trajectories resulting in the most sizeable strains. Some non-trivial behaviors are unveiled, showing that highly eccentric events can be more easily detectable than parabolic ones. In particular, we underline the correct way to extend formulas from hyperbolic to parabolic orbits. Our reasonings are as general as possible, and we make a point of explaining our considerations pedagogically.

"If one could ever prove the existence of gravitational waves, the processes responsible for their generation would probably be much more curious and interesting than even the waves themselves." (Gustav Mie, 1868 - 1957) The discovery of gravitational waves has opened new windows on astrophysics, cosmology and physics beyond the Standard Model (BSM). Measurements by the LIGO, Virgo and KAGRA Collaborations of stellar-mass binaries and neutron star mergers have shown that gravitational waves travel at close to the velocity of light, and also constrain BSM possibilities such as a graviton mass and Lorentz violation in gravitational wave propagation. Follow-up measurements of neutron star mergers have provided evidence for the production of heavy elements, possibly including some essential for human life. The gravitational waves in the nanoHz range observed by Pulsar Timing Arrays (PTAs) may have been emitted by supermassive black hole binaries, but might also have originated from BSM cosmological scenarios such as cosmic strings, or phase transitions in the early Universe. The answer to the question in the title may be provided by gravitational-wave detectors at higher frequencies, such as LISA and atom interferometers.