Papers for Thursday, Nov 24 2022

We present results from the first public data release of the MaNGA-ARO Survey of CO Targets (MASCOT), focussing our study on galaxies whose star-formation rates and stellar masses place them below the ridge of the star-forming Main Sequence. In optically-selected type 2 AGN/LINERs/Composites, we find an empirical relation between gas-phase metallicity gradients $\nabla Z$ and global molecular gas depletion times $t_\mathrm{dep} = M_{H_2}$/SFR with "more quenched" systems showing flatter/positive gradients. Our results are based on the O3N2 metallicity diagnostic (applied to star-forming regions within a given galaxy) which was recently suggested to also be robust against emission by diffuse ionised gas (DIG) and low-ionisation nuclear emission regions (LINERs). We conduct a systematic investigation into possible drivers of the observed $\nabla Z$ - $t_\mathrm{dep}$ relation (ouflows, gas accretion, in-situ star formation, mergers, and morphology). We find a strong relation between $\nabla Z$ or $t_\mathrm{dep}$ and centralised outflow strength traced by the [OIII] velocity broadening. We also find signatures of suppressed star-formation in the outskirts in AGN-like galaxies with long depletion times and an enhancement of metals in the outer regions. We find no evidence of inflows impacting the metallicity gradients, and none of our results are found to be significantly affected by merger activity or morphology. We thus conclude that the observed $\nabla Z$ - $t_\mathrm{dep}$ relation may stem from a combination of metal redistribution via weak feedback, and a connection to in-situ star formation via a resolved mass-metallicity-SFR relation.

We develop a pipeline to set new constraints on scale-independent modified gravity, from the galaxy power spectrum in redshift space of BOSS DR12. The latter is modelled using the effective field theory of large-scale structure up to 1-loop order in perturbation theory. We test our pipeline on synthetic and simulated data, to assess systematic biases on the inferred cosmological parameters due to marginalization and theoretical errors, and we apply it to the normal branch of the DGP model with a $\Lambda$CDM background. We observe biased posteriors due to the strong degeneracy between the nDGP parameter $\Omega_{\rm rc}$ and the primordial amplitude of fluctuations $A_s$. Fixing the latter to the Planck central value, we obtain $\Omega_{\rm rc}\lesssim 0.2$ at 95$\%$ C.L. We also discuss a procedure to alleviate the prior dependence of this bound.

Models of prestellar cores often assume that the cores are isolated from their environment - material outside the core boundary plays no role in the subsequent evolution. This is unlikely to be the case in reality, where cores are located within hierarchically substructured molecular clouds. We investigate the dynamical and chemical evolution of prestellar cores, modelled as Bonnor-Ebert spheres, and show that the density of the ambient medium has a large impact on the resulting chemical properties of the cores. Models embedded in high-density, low-temperature surroundings have greatly enhanced abundances of several molecules, such as CO and CS, compared to models with more diffuse surroundings, corresponding to relatively isolated cores. The predicted intensities and profile shapes of molecular lines are also affected. The density of the ambient medium has a stronger effect on the chemical evolution than whether the cores are initially in or out of equilibrium. This suggests that the impact of environment cannot be neglected when modelling chemistry in prestellar cores; the results of these models are highly sensitive to the assumptions made about the core surroundings.

Herbig stars can be classified into group I and group II depending on the shape of the far-IR excess from the spectral energy distribution. This separation may be evolutionary and related to the vertical structure of these disks. We aim to determine the emission height of Herbig disks and compare the resulting vertical extent of both groups. ALMA Band 6 observations of 12CO emission lines at sufficient velocity and spatial resolution of eight Herbig disks (four group I and four group II sources) are used to determine the emission heights from the channel maps via geometrical methods developed in other works. We find that all group I disks are vertically extended with a height to radius ratio of at least 0.25, and for three of the disks the gas emission profile can be traced out to 200-500 au. The group II disks are divided between MWC 480 and HD 163296 which have similar emission height profiles as the group I disks, and AK Sco and HD 142666 which are very flat (not exceeding a height of 10 au) and more compact (<200 au in size). The brightness temperatures show no differences between the disks when the luminosity of the host star is accounted for. Our findings agree with previous work suggesting that group I disks are vertically extended and that group II disks are either large and self-shadowed or compact. Both MWC 480 and HD 163296 could be precursors of group I disks, which we see now before a cavity has formed that would allow irradiation of the outer parts of the disk. The very flat disks AK Sco and HD 142666 could be due to significant settling because of the advanced age of these disks (~20 instead of <10 Myr). These large differences in vertical structures are not reflected in the spectral energy distributions of these disks. More and deeper observations at higher spatial and velocity resolution are necessary to further characterize the Herbig sub-groups.

We aim to reveal the structure and kinematics of the Outer-Scutum-Centaurus (OSC) arm located on the far side of the Milky Way through very long baseline interferometry (VLBI) astrometry using KaVA, which is composed of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We report the proper motion of a 22 GHz H$_{2}$O maser source, which is associated with the star-forming region G034.84$-$00.95, to be ($\mu_{\alpha} \rm{cos}\delta$, $\mu_{\delta}$) = ($-$1.61$\pm$0.18, $-$4.29$\pm$0.16) mas yr$^{-1}$ in equatorial coordinates (J2000). We estimate the 2D kinematic distance to the source to be 18.6$\pm$1.0 kpc, which is derived from the variance-weighted average of kinematic distances with LSR velocity and the Galactic-longitude component of the measured proper motion. Our result places the source in the OSC arm and implies that G034.84$-$00.95 is moving away from the Galactic plane with a vertical velocity of $-$38$\pm$16 km s$^{-1}$. Since the H I supershell GS033+06$-$49 is located at a kinematic distance roughly equal to that of G034.84$-$00.95, it is expected that gas circulation occurs between the outer Galactic disk around G034.84$-$00.95 with a Galactocentric distance of 12.8$^{+1.0}_{-0.9}$ kpc and halo. We evaluate possible origins of the fast vertical motion of G034.84$-$00.95, which are (1) supernova explosions and (2) cloud collisions with the Galactic disk. However, neither of the possibilities are matched with the results of VLBI astrometry as well as spatial distributions of H II regions and H I gas.

Electromagnetic radiation plays a crucial role in various physical and chemical processes. Hence, almost all astrophysical simulations require some form of radiative transfer model. Despite many innovations in radiative transfer algorithms and their implementation, realistic radiative transfer models remain very computationally expensive, such that one often has to resort to approximate descriptions. The complexity of these models makes it difficult to assess the validity of any approximation and to quantify uncertainties on the model results. This impedes scientific rigour, in particular, when comparing models to observations, or when using their results as input for other models. We present a probabilistic numerical approach to address these issues by treating radiative transfer as a Bayesian linear regression problem. This allows us to model uncertainties on the input and output of the model with the variances of the associated probability distributions. Furthermore, this approach naturally allows us to create reduced-order radiative transfer models with a quantifiable accuracy. These are approximate solutions to exact radiative transfer models, in contrast to the exact solutions to approximate models that are often used. As a first demonstration, we derive a probabilistic version of the method of characteristics, a commonly-used technique to solve radiative transfer problems.

The demographics of the production and escape of ionizing photons from UV-faint early galaxies is a key unknown in discovering the primary drivers of reionization. With the advent of JWST it is finally possible to observe the rest-frame optical nebular emission from individual sub-L$^*$ z>3 galaxies to measure the production of ionizing photons, $\xi_\mathrm{ion}$. Here we study a sample of 380 z~3-7 galaxies spanning -23 <M$_\mathrm{UV}$ < -15.5 (median M$_\mathrm{UV}\approx$ -18) with deep multi-band HST and JWST/NIRCam photometry covering the rest-UV to optical from the GLASS and UNCOVER JWST surveys. Our sample includes 109 galaxies with Lyman-alpha emission detected in MUSE spectroscopy. We use H-alpha fluxes inferred from NIRCam photometry to estimate the production rate of ionizing photons which do not escape these galaxies $\xi_\mathrm{ion}(1-f_\mathrm{esc})$. We find median $\log_{10}\xi_\mathrm{ion}(1-f_\mathrm{esc})=25.33\pm 0.47$, with a broad intrinsic scatter 0.42 dex, implying a broad range of galaxy properties and ages in our UV-faint sample. Galaxies detected with Lyman-alpha have ~0.1 dex higher $\xi_\mathrm{ion}(1-f_\mathrm{esc})$, which is explained by their higher H-alpha EW distribution, implying younger ages, higher sSFR and thus more O/B stars. We find significant trends of increasing $\xi_\mathrm{ion}(1-f_\mathrm{esc})$ with increasing H-alpha EW, decreasing UV luminosity, and decreasing UV slope, implying the production of ionizing photons is enhanced in young, low metallicity galaxies. We find no significant evidence for sources with very high ionizing escape fraction ($f_\mathrm{esc}$>0.5) in our sample, based on their photometric properties, even amongst the Lyman-alpha selected galaxies. This work demonstrates that considering the full distribution of $\xi_\mathrm{ion}$ across galaxy properties is important for assessing the primary drivers of reionization.

High-energy particles hitting the upper atmosphere of the Earth produce extensive air showers that can be detected from the ground level using imaging atmospheric Cherenkov telescopes. The images recorded by Cherenkov telescopes can be analyzed to separate gamma-ray events from the background hadron events. Many of the methods of analysis require simulation of massive amounts of events and the corresponding images by the Monte Carlo method. However, Monte Carlo simulation is computationally expensive. The data simulated by the Monte Carlo method can be augmented by images generated using faster machine learning methods such as generative adversarial networks or conditional variational autoencoders. We use a conditional variational autoencoder to generate images of gamma events from a Cherenkov telescope of the TAIGA experiment. The variational autoencoder is trained on a set of Monte Carlo events with the image size, or the sum of the amplitudes of the pixels, used as the conditional parameter. We used the trained variational autoencoder to generate new images with the same distribution of the conditional parameter as the size distribution of the Monte Carlo-simulated images of gamma events. The generated images are similar to the Monte Carlo images: a classifier neural network trained on gamma and proton events assigns them the average gamma score 0.984, with less than 3% of the events being assigned the gamma score below 0.999. At the same time, the sizes of the generated images do not match the conditional parameter used in their generation, with the average error 0.33.

We present a significantly updated CO2 altitude profile for Venus (64.2-0.9 km) and provide support for a potential deep lower atmospheric haze of particles (17 km and lower). We extracted this information by developing a new analytical model for mass spectra obtained by the Pioneer Venus Large Probe (PVLP) Neutral Mass Spectrometer (LNMS). Our model accounts for changes in LNMS configuration and output during descent and enables the disentanglement of isobaric species via a data fitting routine that adjusts for mass-dependent changes in peak shape. The model yields CO2 in units of density (kg m-3), isotope ratios for 13C/12C and 18O/16O, and 14 measures of CO2 density across 55.4-0.9 km, which represents the most complete altitude profile for CO2 at 60 km towards the surface to date. The CO2 density profile is also consistent with the pressure, temperature, and volumetric gas measurements from the PVLP and VeNeRa spacecraft. Nominal and low-noise operations for the LNMS mass analyzer are supported by the behaviors (e.g., ionization yields, fragmentation yields, and peak shapes) of several internal standards (e.g., CH3+, CH4+, 40Ar+, 136Xe2+, and 136Xe+), which were tracked across the descent. Lastly, our review of the CO2 profile and LNMS spectra reveals hitherto unreported partial and rapidly clearing clogs of the inlet in the lower atmosphere, along with several ensuing data spikes at multiple masses. Together, these observations suggest that atmospheric intake was impacted by particles at 17 km (and lower) and that rapid particle degradation at the inlet yielded a temporary influx of mass signals into the LNMS.

In the warm dark matter scenario, the Press-Schechter formalism is valid only for galaxy masses greater than the "velocity dispersion cut-off". In this work we extend the predictions to masses below the velocity dispersion cut-off, and thereby address the "Missing Satellites Problem", and the rest-frame ultra-violet luminosity cut-off required to not exceed the measured reionization optical depth. We find agreement between predictions and observations of these two phenomena. As a by-product, we obtain the empirical Tully-Fisher relation from first principles.

We find that the chemical abundances and dynamics of APOGEE and GALAH stars in the local stellar halo are inconsistent with a scenario in which the inner halo is primarily composed of debris from a single, massive, ancient merger event, as has been proposed to explain the Gaia-Enceladus/Gaia Sausage (GSE) structure. The data contains trends of chemical composition with energy which are opposite to expectations for a single massive, ancient merger event, and multiple chemical evolution paths with distinct dynamics are present. We use a Bayesian Gaussian mixture model regression algorithm to characterize the local stellar halo, and find that the data is best fit by a model with four components. We interpret these components as the VRM, Cronus, Nereus, and Thamnos; however, Nereus and Thamnos likely represent more than one accretion event because the chemical abundance distributions of their member stars contain many peaks. Although the Cronus and Thamnos components have different dynamics, their chemical abundances suggest they may be related. We show that the distinct low and high alpha halo populations from Nissen & Schuster (2010) are explained by VRM and Cronus stars, as well as some in-situ stars. Because the local stellar halo contains multiple substructures, different popular methods of selecting GSE stars will actually select different mixtures of these substructures, which may change the apparent chemodynamic properties of the selected stars. We also find that the Splash stars in the solar region are shifted to higher v_phi and slightly lower [Fe/H] than previously reported.

We present a two-epoch Hubble Space Telescope (HST) near-infrared (NIR) study of NGC 2071 IR highlighting HOPS 361-C, a protostar producing an arced 0.2 parsec-scale jet. Proper motions for the brightest knots decrease from 350 to 100 km/s with increasing distance from the source. The [Fe II] and Pa$\beta$ emission line intensity ratio gives a velocity jump through each knot of 40-50 km/s. We show a new [O I] 63 $\rm \mu$m spectrum taken with the German REciever for Astronomy at Terahertz frequencies (GREAT) instrument aboard Stratospheric Observatory for Infrared Astronomy (SOFIA), which give a low jet inclination. Proper motions and jump velocities then estimate total flow speed throughout the jet. We model knot positions and speeds with a precessing jet that decelerates within the host molecular cloud. The measurements are matched with a precession period of a few thousand years and half opening angle of 15$\rm\deg$. The [Fe II] 1.26 $\rm \mu$m to 1.64 $\rm \mu$m line intensity ratio gives the extinction to each knot ranging from 5-30 mag. Relative to $\sim$14 mag of extinction through the cloud from C$^{18}$O emission maps, the jet is well embedded at a fractional depth from 1/5 to 4/5, and can interact with the cloud. Our model suggests the jet is locally dissipated over 0.2 pc. This may be because knots sweep through a wide angle, giving the cloud time to fill in cavities opened by the jet. This contrasts with nearly unidirectional protostellar jets that puncture host clouds and can propagate significantly further than a quarter pc.

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range of $z=6-7$ is poorly constrained. We present new constraints on $\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$ quasars at $5.73<z<7.09$. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of "dark" pixels identified in the Lyman $\alpha$ (Ly$\alpha$) and Lyman $\beta$ (Ly$\beta$) forests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction at $z\sim6.2-6.8$ using quasar absorption measurements. Our results give upper limits of $\bar{x}_{\rm HI}(z=6.3) < 0.79\pm0.04$ (1$\sigma$), $\bar{x}_{\rm HI} (z=6.5) < 0.87\pm0.03$ (1$\sigma$), and $\bar{x}_{\rm HI} (z=6.7) < 0.94^{+0.06}_{-0.09}$ (1$\sigma$). The dark pixel fractions at $z>6.1$ are consistent with the redshift evolution of the neutral fraction of the IGM derived from the Planck 2018.

The ALMA observatory is now putting more focus on high-frequency observations (frequencies from 275-950 GHz). However, high-frequency observations often suffer from rapid variations in atmospheric opacity that directly affect the system temperature $T_{sys}$. Current observations perform discrete atmospheric calibrations (Atm-cals) every few minutes, with typically 10-20 occurring per hour for high frequency observation and each taking 30-40 seconds. In order to obtain more accurate flux measurements and reduce the number of atmospheric calibrations (Atm-cals), a new method to monitor $T_{sys}$ continuously is proposed using existing data in the measurement set. In this work, we demonstrate the viability of using water vapor radiometer (WVR) data to track the $T_{sys}$ continuously. We find a tight linear correlation between $T_{sys}$ measured using the traditional method and $T_{sys}$ extrapolated based on WVR data with scatter of 0.5-3%. Although the exact form of the linear relation varies among different data sets and spectral windows, we can use a small number of discrete $T_{sys}$ measurements to fit the linear relation and use this heuristic relationship to derive $T_{sys}$ every 10 seconds. Furthermore, we successfully reproduce the observed correlation using atmospheric transmission at microwave (ATM) modeling and demonstrate the viability of a more general method to directly derive the $T_{sys}$ from the modeling. We apply the semi-continuous $T_{sys}$ from heuristic fitting on a few data sets from Band 7 to Band 10 and compare the flux measured using these methods. We find the discrete and continuous $T_{sys}$ methods give us consistent flux measurements with differences up to 5%. Furthermore, this method has significantly reduced the flux uncertainty due to $T_{sys}$ variability for one dataset, which has large precipitable water vapor (PWV) fluctuation, from 10% to 0.7%.

We present 65 extremely dust-obscured galaxies from the UltraVISTA DR3 survey of the COSMOS field at $1<z<4$. In contrast to other studies of dusty galaxies, we select our sample based on dust attenuation measured by UV-MIR spectral energy distribution (SED) modeling that allows for extreme attenuation levels. We construct our sample by making cuts at $1 \le z \le 4$, A$_V \ge 3$, and log(M$_*$/M$_\odot$)$ \ge 10.5$. This method reliably selects galaxies exhibiting independent indicators of significant dust content, including FIR detection rates. We perform panchromatic SED modeling with matched $Herschel$ photometry and find stellar and dust properties that differ from typical sub-millimeter galaxy (SMG) samples as well as $Herschel$ sources matched in redshift and stellar mass. Our sources have lower star formation rates and higher A$_V$ than SMGs, but comparable total IR luminosities. Most of our sample falls on or near the star-forming main sequence for this redshift range. Finally, we perform a morphological analysis with GALFIT using the $K_S$-band images and $Hubble$ $F814W$ and $F160W$ imaging when available. Typical axis ratios of $\sim 0.4$ suggest disk-like morphology for the majority of our sources, and we note only three apparent merging systems. Our sample generally agrees with the size-mass relation for star-forming galaxies, with a tail extending to smaller sizes. We conclude that the most heavily obscured galaxies in this redshift range share many characteristics with typical star-forming galaxies, forming a population of dusty galaxies that overlaps, but is not encompassed by, those selected through dust emission.

Accurate modeling of the Sun's coronal magnetic field and solar wind structures require inputs of the solar global magnetic field, including both the near and far sides, but the Sun's far-side magnetic field cannot be directly observed. However, the Sun's far-side active regions are routinely monitored by helioseismic imaging methods, which only require continuous near-side observations. It is therefore both feasible and useful to estimate the far-side magnetic-flux maps using the far-side helioseismic images despite their relatively low spatial resolution and large uncertainties. In this work, we train two machine-learning models to achieve this goal. The first machine-learning training pairs simultaneous SDO/HMI-observed magnetic-flux maps and SDO/AIA-observed EUV 304$\r{A}$ images, and the resulting model can convert 304$\r{A}$ images into magnetic-flux maps. This model is then applied on the STEREO/EUVI-observed far-side 304$\r{A}$ images, available for about 4.3 years, for the far-side magnetic-flux maps. These EUV-converted magnetic-flux maps are then paired with simultaneous far-side helioseismic images for a second machine-learning training, and the resulting model can convert far-side helioseismic images into magnetic-flux maps. These helioseismically derived far-side magnetic-flux maps, despite their limitations in spatial resolution and accuracy, can be routinely available on a daily basis, providing useful magnetic information on the Sun's far side using only the near-side observations.

We compare the core-collapse evolution of a pair of 15.8 $M_\odot$ stars with significantly different internal structures, a consequence of bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79 $M_\odot $ progenitors have core masses of 1.47 and 1.78 $M_\odot$ and compactness parameters $\xi_{1.75}$ of 0.302 and 0.604. The core collapse simulations are carried out in 2D to nearly 3 s post-bounce and show substantial differences in the times of shock revival and explosion energies. The 15.78 $M_\odot$ model explodes promptly at 120 ms post-bounce when a strong density decrement at the Si--Si/O shell interface encounters the stalled shock. The 15.79 $M_\odot$ model, which lacks the density decrement, takes 100 ms longer to explode but ultimately produces a more powerful explosion. Larger mass accretion rate of the 15.79 $M_\odot$ model during the first 0.8 s post-bounce results in larger $\nu_{e}$/$\bar \nu_{e}$ luminosities and rms energies. The $\nu_{e}$/$\bar \nu_{e}$ luminosities and rms energies arising from the inner core are also larger in the 15.79 $M_\odot$ model throughout due to the larger negative temperature gradient of this core due to greater adiabatic compression. Larger luminosities and rms energies in the 15.79 $M_\odot$ model and a flatter and higher density heating region, result in more energy deposition behind the shock and more ejected matter with higher enthalpy. We find the ejected $^{56}$Ni mass of the 15.79 $M_\odot$ model is more than double that of the 15.78 $M_\odot$ model. Most of the ejecta in both models is moderately proton-rich, though counterintuitively the highest electron fraction ($Y_e=0.61$) ejecta in either model is in the less energetic 15.78 $M_\odot$ model while the lowest electron fraction ($Y_e=0.45$) ejecta in either model is in the 15.79 $M_\odot$ model.

We present new isochrone fits to the colour--magnitude diagrams of the Galactic globular clusters NGC\,6362 and NGC\,6723. We utilize 22 and 26 photometric filters for NGC\,6362 and NGC\,6723, respectively, from the ultraviolet to mid-infrared using data sets from {\it HST}, {\it Gaia}, unWISE, and other photometric sources. We use models and isochrones from the Dartmouth Stellar Evolution Database (DSED) and Bag of Stellar Tracks and Isochrones (BaSTI) for $\alpha$--enhanced [$\alpha$/Fe]$=+0.4$ and different helium abundances. The metallicities [Fe/H]$=-1.04\pm0.07$ and $-1.09\pm0.06$ are derived from the red giant branch slopes in our fitting for NGC\,6362 and NGC\,6723, respectively. They agree with spectroscopic estimates from the literature. We find a differential reddening up to $\Delta E(B-V)=0.13$ mag in the NGC\,6723 field due to the adjacent Corona Australis cloud complex. We derive the following for NGC\,6362 and NGC\,6723, respectively: distances $7.75\pm0.03\pm0.15$ (statistic and systematic error) and $8.15\pm0.04\pm0.15$ kpc; ages $12.0\pm0.1\pm0.8$ and $12.4\pm0.1\pm0.8$ Gyr; extinctions $A_\mathrm{V}=0.19\pm0.04\pm0.06$ and $0.24\pm0.03\pm0.06$ mag; reddenings $E(B-V)=0.056\pm0.01\pm0.02$ and $0.068\pm0.01\pm0.02$ mag. DSED provides systematically lower [Fe/H] and higher reddenings than BaSTI. However, the models agree in their relative estimates: NGC\,6723 is $0.44\pm0.04$ kpc further, $0.5\pm0.1$ Gyr older, $\Delta E(B-V)=0.007\pm0.002$ more reddened, and with $0.05\pm0.01$ dex lower [Fe/H] than NGC\,6362. The lower metallicity and greater age of NGC\,6723 with respect to NGC\,6362 explain their horizontal branch morphology difference. This confirms age as the second parameter for these clusters. We provide lists of the cluster members from the {\it Gaia} Data Release 3.

Soliton in the hostile turbulent $\Psi$DM halo of a galaxy agitates with various kinds of excitation, and the soliton even breathes heavily under great stress. A theory of collective excitation for a $\Psi$DM soliton is presented. The collective excitation has different degrees of coupling to negative energy modes, where lower-order excitation generally necessitates more negative energy coupling. A constrained variational principle is developed to assess the frequencies and mode structures of small-amplitude perturbations. The predicted frequencies are in good agreement with those found in simulations. Soliton breathing at amplitudes on the verge of breakup is also a highlight of this work. Even in this extreme nonlinear regime, the wave function perturbation amplitudes are moderate. The simulation data shows a stable oscillation with frequency weakly dependent on the oscillation amplitude, and hints a self-consistent quasi-linear model for the wave function that accounts for modifications in the ground state wave function and the equilibrium density. The mock solution, constructed from the simulation data, can shed lights on the dynamics of the large-amplitude breathing soliton and supports the quasi-linear model, as evidenced by its ability to well predict the nonlinear eigenfrequency shifts and large-amplitude breathing frequency observed in simulations.

Extracting the maximum amount of cosmological and astrophysical information from upcoming large-scale surveys remains a challenge. This includes evaluating the exact likelihood, parameter inference and generating new diverse synthetic examples of the incoming high-dimensional data sets. In this work, we propose the use of normalizing flows as a generative model of the neutral hydrogen (HI) maps from the CAMELS project. Normalizing flows have been very successful at parameter inference and generating new, realistic examples. Our model utilizes the spatial structure of the HI maps in order to faithfully follow the statistics of the data, allowing for high-fidelity sample generation and efficient parameter inference.

The Audible Universe project aims at making dialogue between two scientific domains investigating two distinct research objects, briefly said, Stars and Sound. It has been instantiated within a collaborative workshop that started to mutually acculturate both communities, by sharing and transmitting respective knowledge, skills and practices. One main outcome of this exchange was a global view on the astronomical data sonification paradigm that allowed to observe either the diversity of tools, uses and users (including visually-impaired people), but also the current limitations and potential ways of improvement. From this perspective, the current paper presents basic elements gathered and contextualised by sound experts in their respective fields (sound perception / cognition, sound design, psychoacoustics, experimental psychology), in order to anchor sonification for astronomy in a more well-informed, methodological and creative process.

In this Letter, we investigate jet-launching abilities of Bondi-like accretion flows with zero or low specific angular momentum by performing 3D general relativistic magnetohydrodynamic simulations. In order to check if relativistic jets can be launched magnetically, we thread the accretion flow with large-scale poloidal magnetic field, and choose a rapidly spinning black hole. We demonstrate that the magnitude of the initial gas specific angular momentum primarily controls whether the disk can reach and sustain the magnetically arrested disk (MAD) state that launches very powerful jets, at $\gtrsim 100\%$ energy efficiency. We find that MAD forms in the presence of even a very small amount of specific angular momentum, and episodic jets with an average energy efficiency of $\sim 10\%$ can still form even when the gas has zero initial angular momentum. Our results give plausible explanations to why jets can be produced from various astrophysical systems that lack large gas specific angular momenta, such as Sgr A*, wind-fed X-ray binaries, tidal disruption events, and long-duration gamma-ray bursts.

In this study, we continue our project to search for unresolved binary and multiple systems in open clusters exploiting the photometric diagram (H-W2)-W1 vs W2-(BP-K) firstly introduced in \citet{Malofeeva+2022}. In particular, here we estimate the binary and multiple star ratios and the distribution of the component mass ratio $q$ in the Galactic clusters Alpha Persei, Praesepe, and NGC 1039. We have modified the procedure outlined in our first study \citep{Malofeeva+2022} making star counts automatic and by introducing bootstrapping for error estimation. Basing on this, we re-investigated the Pleiades star cluster in the same mass range as in our previous work and corrected an inaccuracy in the mass ratio $q$ distribution. The binary and multiple star ratio in the four clusters is then found to lie between 0.45$\pm$0.03 and 0.73$\pm$0.03. On the other hand, the ratio of systems with multiplicity more than 2 is between 0.06$\pm$0.01 and 0.09$\pm$0.02. The distribution of the component mass ratio $q$ is well fitted with a Gaussian having the mode between 0.22$\pm$0.04 and 0.52$\pm$0.01 and the dispersion between 0.10$\pm$0.02 and 0.35$\pm$0.07. All clusters show a large number of the very low-mass secondary components in the binary systems with primary components below 0.5 $M_{\odot}$.

We present a pilot study to assess the potential of Hyper Suprime-Cam Public Data Release 2 (HSC-PDR2) images for the analysis of extended faint structures within groups of galaxies. We examine the intra-group light (IGL) of the group 400138 ($M_{\rm{dyn}}= 1.3 \pm 0.5 \times 10^{13} $M$_{\odot}$, $z\sim 0.2$) from the Galaxy And Mass Assembly (GAMA) survey using Hyper-Suprime Cam Subaru Strategic Program Public Data Release 2 (HSC-PDR2) images in $g$, $r$, and $i$ bands. We present the most extended IGL measurement to date, reaching down to $\mu_{g}^{\rm{lim}}=30.76$ mag arcsec$^{-2}$ ($3 \sigma$; $10 \times 10$ arcsec$^{2}$) at a semi-major axis of 275 kpc. The IGL shows mean colour values of $g-i=0.92$, $g-r=0.60$, and $r-i=0.32$ ($\pm$0.01). The IGL stellar populations are younger ($2-2.5$ Gyr) and less metal-rich ([Fe/H] $ \sim -$0.4) than those of the host group galaxies. We find a range of IGL fractions as a function of total group luminosity of $\sim 2-36 \%$ depending on the definition of IGL, with larger fractions the bluer the observation wavelength. The early-type to late-type galaxy ratio suggests that 400138 is a more evolved group, dominated by ETGs, and the IGL fraction agrees with that of other similarly evolved groups. These results are consistent with tidal stripping of the outer parts of Milky Way-like galaxies as the main driver of the IGL build-up. This is supported by the detection of substructure in the IGL towards the galaxy member 1660615 suggesting a recent interaction ($<1$ Gyr ago) of that galaxy with the core of the group.

We argue that resonant friction has a dramatic effect on a disc whose rotation direction is misaligned with that of its host nuclear star cluster. The disc's gravity causes gravitational perturbation of the cluster that in turn exerts a strong torque back onto the disc. We argue that this torque may be responsible for the observed disruption of the clockwise disc of young stars in the Galactic Center, and show in numerical experiments that it produces the observed features in the distribution of the stars' angular momenta. More generally, we speculate that the rotation of nuclear star clusters has a stabilizing effect on the orientation of transient massive accretion discs around the supermassive black holes residing in their centers, and thus on the directions and magnitudes of the black-hole spins.

The Sun is the source of different types of radio bursts that are associated with solar flares, for example. Among the most frequently observed phenomena are type III solar bursts. Their radio images at low frequencies (below 100 MHz) are relatively poorly studied due to the limitations of legacy radio telescopes. We study the general characteristics of types IIIb and U with stria structure solar radio bursts in the frequency range of 20 - 80 MHz, in particular the source size and evolution in different altitudes, as well as the velocity and energy of electron beams responsible for their generation. In this work types IIIb and U with stria structure radio bursts are analyzed using data from the LOFAR telescope including dynamic spectra and imaging observations, as well as data taken in the X-ray range (GOES and RHESSI satellites) and in the extreme ultraviolet (SDO satellite). In this study we determined the source size limited by the actual shape of the contour at particular frequencies of type IIIb and U solar bursts in a relatively wide frequency band from 20 to 80 MHz. Two of the bursts seem to appear at roughly the same place in the studied active region and their source sizes are similar. It is different in the case of another burst, which seems to be related to another part of the magnetic field structure in this active region. The velocities of the electron beams responsible for the generation of the three bursts studied here were also found to be different.

We present a consistent one-loop calculation for the inflationary tensor power spectrum in the presence of an excited spectator scalar field using the in-in formalism. We find that the super-horizon primordial power spectrum of the tensor mode can be scale-invariantly enhanced or reduced by the loop effects of a subhorizon scalar field. Our calculation also includes the scalar-induced gravitational wave spectrum classically computed in the previous literature, which is significant only near the scales where the scalar field is amplified. The super-horizon enhancement is a higher-order effect of the interaction Hamiltonian, which can be understood as a Bogoliubov transformation introduced by nonlinear interactions. On the other hand, the scale-invariant reduction of the tensor power spectrum may occur due to the fourth-order scalar-scalar-tensor-tensor coupling. This phenomenon can be understood as the evolution of an anisotropic Bianchi type-I background in the separate universe approach. Our result suggests that large-scale measurements may indirectly test the dramatic effects of small-scale cosmological perturbations through loop corrections. This possibility opens a new ground in probing the small-scale physics of the primordial Universe through gravitational wave detectors of cosmological scales.

We study the circular polarization of the Mn I resonance lines at 279.56, 279.91, and 280.19 nm (hereafter, UV multiplet) by means of radiative transfer modeling. In 2019, the CLASP2 mission obtained unprecedented spectropolarimetric data in a region of the solar ultraviolet including the Mg II h and k resonance lines and two lines of a subordinate triplet, as well as two Mn I resonance lines. The first analysis of such data, in particular those corresponding to a plage region, allowed the inference of the longitudinal magnetic field from the photosphere to the upper chromosphere just below the transition region. This was achieved by applying the weak field approximation to the circular polarization profiles of the Mg II and Mn I lines. While the applicability of this approximation to the Mg II lines was already demonstrated in previous works, this is not the case for the Mn I UV multiplet. These lines are observed as absorptions between the Mg II h and k lines, a region whose intensity is shaped by their partial frequency redistribution effects. Moreover, the only Mn I stable isotope has nuclear spin $I=5/2$ and thus hyperfine structure must be, a priori, taken into account. Here we study the generation and transfer of the intensity and circular polarization of the Mn I resonance lines accounting for these physical ingredients. We analyze their sensitivity to the magnetic field by means of their response function, and we demonstrate the applicability of the weak field approximation to determine the longitudinal component of the magnetic field.

Optical properties are required for the correct understanding and modelling of protoplanetary and debris discs. By assuming that comets are the most pristine bodies in the solar system, our goal is to derive optical constants of real protoplanetary material. We determine the complex index of refraction of the near-surface material of comet 67P/Churyumov-Gerasimenko by fitting the sub-millimetre/millimetre observations of the thermal emission of the comet's sub-surface made by the Microwave Instrument for the Rosetta Orbiter (MIRO) with synthetic temperatures derived from a thermophysical model and radiative-transfer models. According to the two major formation scenarios of comets, we model the sub-surface layers to consist of pebbles as well as of homogeneously packed dust grains. In the case of a homogeneous dusty surface material, we find a solution for the length-absorption coefficient of $\alpha \approx 0.22~\mathrm{cm^{-1}}$ for a wavelength of 1.594 mm and $\alpha \geq 3.84~\mathrm{cm^{-1}}$ for a wavelength of 0.533 mm and a constant thermal conductivity of $0.006~\mathrm{Wm^{-1}K^{-1}}$. For the pebble scenario, we find for the pebbles and a wavelength of 1.594 mm a complex refractive index of $n = (1.074 - 1.256) + \mathrm{i} \, (2.580 - 7.431)\cdot 10^{-3}$ for pebble radii between 1 mm and 6 mm. Taking into account other constraints, our results point towards a pebble makeup of the cometary sub-surface with pebble radii between 3 mm and 6 mm. The derived real part of the refractive index is used to constrain the composition of the pebbles and their volume filling factor. The optical and physical properties are discussed in the context of protoplanetary and debris disc observations.

The 21 cm spectral line emission of atomic neutral hydrogen (HI) is one of the primary wavelengths observed in radio astronomy. However, the signal is intrinsically faint and the HI content of galaxies depends on the cosmic environment, requiring large survey volumes and survey depth to investigate the HI Universe. As the amount of data coming from these surveys continues to increase with technological improvements, so does the need for automatic techniques for identifying and characterising HI sources while considering the tradeoff between completeness and purity. This study aimed to find the optimal pipeline for finding and masking the most sources with the best mask quality and the fewest artefacts in 3D neutral hydrogen cubes. Various existing methods were explored in an attempt to create a pipeline to optimally identify and mask the sources in 3D neutral hydrogen 21 cm spectral line data cubes. Two traditional source-finding methods were tested, SoFiA and MTObjects, as well as a new supervised deep learning approach, in which a 3D convolutional neural network architecture, known as V-Net was used. These three source-finding methods were further improved by adding a classical machine learning classifier as a post-processing step to remove false positive detections. The pipelines were tested on HI data cubes from the Westerbork Synthesis Radio Telescope with additional inserted mock galaxies. SoFiA combined with a random forest classifier provided the best results, with the V-Net-random forest combination a close second. We suspect this is due to the fact that there are many more mock sources in the training set than real sources. There is, therefore, room to improve the quality of the V-Net network with better-labelled data such that it can potentially outperform SoFiA.

Stealthy Coronal Mass Ejections (CMEs), lacking low coronal signatures, may result in significant geomagnetic storms. However, the mechanism of stealthy CMEs is still highly debated. In this work, we investigate whether there are differences between the stealthy and standard CMEs in terms of their dynamic behaviors. Seven stealthy and eight standard CMEs with slow speeds are selected. We calculate two-dimensional speed distributions of CMEs based on the cross-correlation method, rather than the unidimensional speed, and further obtain more accurate distributions and evolution of CME mechanical energies. Then we derive the CME driving powers and correlate them with CME parameters (total mass, average speed, and acceleration) for standard and stealthy CMEs. Besides, we study the forces that drive CMEs, namely, the Lorentz force, gravitational force, and drag force due to the ambient solar wind near the Sun. The results reveal that both the standard and stealthy CMEs are propelled by the combined action of those forces in the inner corona. The drag force and gravitational force are comparable with the Lorentz force. However, the impact of the drag and Lorentz forces on the global evolution of the stealthy CMEs is significantly weaker than that of the standard CMEs.

Several studies have shown the influence of the relative streaming velocity (SV) between baryons and dark matter on the formation of structures. For the first time, we constrain the local value of the SV in which the Milky Way was formed. We use the semi-analytical model A-SLOTH to simulate the formation of Milky Way-like galaxies. The high resolution in mass and time of the dark matter merger trees from the Caterpillar simulation enables to accurately model star formation in the smallest progenitor halos at high redshift. The efficient semi-analytical nature of A-SLOTH allows us to run many simulations with various values of the local SV. Our investigation on the influence of the SV shows that it delays star formation at high redshift. However, at redshift z=0, the SV has no effect on the total stellar mass in the Milky Way nor its Satellites. We find that extremely metal-poor and ultra metal-poor stars are affected by the SV, and can hence be used to constrain its local value. The local optimal value of the SV is $v_\mathrm{SV} =1.75^{+0.13}_{-0.28}\,\sigma_\mathrm{SV}$, which is based on four independent observables. We further find that the SV decreases the number of luminous Milky Way satellites, but this decrease is not enough to solve the missing satellite problem.

A primordial group of open clusters containing NGC 6871 is confirmed and described through Gaia DR3 data and the previous literature. It is a star-forming complex containing at least six young OCs, including Teutsch 8, FSR 198 and Biurakan 2. Two nearby OCs (Casado 82 and Casado-Hendy 1) are newly identified and studied in detail and found to be also members of the cited group. The parameters of the components are sufficiently similar to postulate the case of at least six clusters born from a single GMC. None of the cluster pairs of the group seems to be an authentic binary cluster, with the possible exception of the candidate pair Teutsch 8/FSR 198. Instead, NGC 6871 seems to be disintegrating, and the primordial group members appear to be dispersing out rapidly. Searching for new open clusters in the vicinity of young or grouped OCs using Gaia data is an efficient strategy to find new associated OCs forming primordial groups.

Coronal Mass Ejections (CMEs) are one of the main drivers of disturbances in the interplanetary space. Strong CMEs, when directed towards the Earth, cause geo-magnetic storms upon interacting with the magnetic field of the Earthand can cause significant damage to our planet and affect everyday life. As such, efficient space weather prediction tools are necessary to forecast the arrival and impact of CME eruptions. Recently, a new heliospheric model Icarus was developed based on MPI-AMRVAC, which is a 3D ideal MHD model for the solar wind and CME propagation, and it introduces advanced numerical techniques to make the simulations more efficient. A cone model is used to study the evolution of the CME through the background solar wind and its arrival and impact at Earth. Grid stretching and AMR are combined in the simulations by using multiple refinement criteria. We compare simulation results to the EUFHORIA model. As a result, the simulations were sped up by a factor of 17 for the most optimal configuration in Icarus. For the cone CME model, we found that limiting the AMR to the region around the CME-driven shock yields the best results. The results modelled by the simulations with radial grid stretching and AMR level 4 are similar to the results provided by the original EUHFORIA and Icarus simulations with the 'standard' resolution and equidistant grids. The simulations with 5 AMR levels yielded better results than the simulations with an equidistant grid and standard resolution. Solution AMR is flexible and provides the user the freedom to modify and locally increase the grid resolution according to the purpose of the simulation. The advanced techniques implemented in Icarus can be further used to improve the forecasting procedures, since the reduced simulation time is essential to make physics-based forecasts less computationally expensive.

The dynamical evolution of globular clusters is theoretically described by a series of well known events typical of N-body systems. Still, the identification of observational signatures able to empirically describe the stage of dynamical evolution of a stellar system of the density typical of a globular cluster, represents a challenge. In this paper we study the dynamical age of the globular clusters Rup 106 and IC 4499. To this aim, we study the radial distribution of the Blue Straggler Stars via the A+ parameter and of the slope of the Main Sequence Mass Function. Both tracers show that Rup 106 and IC 4499 are dynamically young clusters where dynamical friction has just started to segregate massive stars towards the clusters' centre. In fact, we observe that the Blue Straggler stars are more centrally concentrated in both clusters than the reference population. On the same line, we find that in both cases the slope of the mass function significantly decreases as a function of the cluster-centric distances. This result provides additional support for the use of the the radial distribution of the blue stragglers as a powerful observationally convenient indicator of the cluster dynamical age.

Context. It is known that Alfv\'en and magnetoacoustic waves both contribute to the heating of the solar chromosphere and drive plasma outflows. In both cases, the thermalization of the wave energy occurs due to ion-neutral collisions, but the obtained rates of plasma heating cannot explain the observational data. The same is true for the magnitudes of the outflows. Aims. The aim of the present paper is to reexamine two-fluid modeling of Alfv\'en and magnetoacoustic waves in the partially ionized solar chromosphere. We attempt to detect variations in the ion temperature, and vertical plasma flows for different wave combinations. Methods. We performed numerical simulations of the generation and evolution of coupled Alfv\'en and magnetoacoustic waves using the JOANNA code, which solves the two-fluid equations for ions (protons)+electrons and neutrals (hydrogen atoms), coupled by collision terms. Results. We confirm that the damping of impulsively generated small-amplitude waves negligibly affects the chromosphere temperature and generates only slow plasma flows. In contrast, waves generated by large-amplitude pulses significantly increase the chromospheric temperature and result in faster plasma outflows. The maximum heating occurs when the pulse is launched from the center of the photosphere, and the magnitude of the related plasma flows increases with the amplitude of the pulse. Conclusions. Large-amplitude coupled two-fluid Alfv\'en and magnetoacoustic waves can significantly contribute to the heating of the solar chromosphere and to the generation of plasma outflows.

One of the most promising probes to constrain the reionization history of the universe is the power spectrum of neutral hydrogen 21 cm emission fluctuations. The corresponding analyses require computationally efficient modelling of reionization, usually achieved through semi-numerical simulations. We investigate the capability of one such semi-numerical code, SCRIPT, to constrain the reionization parameters. Our study involves creating a mock data set corresponding to the upcoming SKA-Low, followed by a Bayesian inference method to constrain the model parameters. In particular, we explore in detail whether the inferred parameters are unbiased with respect to the inputs used for the mock, and also if the inferences are insensitive to the resolution of the simulation. We find that the model is successful on both fronts. We also develop a simple template model of reionization which can mimic the complex physical processes like inhomogeneous recombinations and radiative feedback and show that it can recover the global reionization history reliably with moderate computational cost. Our results are relevant for constraining reionization using high-quality data expected in future telescopes.

In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compared to typical narrow-band methods are, firstly, generation of two-dimensional templates accounting for pulse shape evolution with frequency and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and dispersion measure (DM). This is the first time that wideband timing has been applied to a uniform data set collected with a single large-fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions including initial noise analysis. Precision of our ToA and DM measurements is in the range of 0.005 $-$ 2.08 $\mu$s and (0.043$-$14.24)$\times10^{-4}$ cm$^{-3}$ pc, respectively, with 94% of the pulsars achieving a median ToA uncertainty of less than 1 $\mu$s.

Some observations and numerical simulations of disc-magnetosphere interaction show that accretion can proceed in the propeller regime. When the Alfv\'en radius is beyond the corotation radius, matter climbs up to the high latitudes where the Alfv\'en surface is inside the equilibrium surface and can accrete. We calculate the fraction of the mass flux in the disc that can accrete onto the neutron star depending on the fastness parameter and the inclination angle between rotation and magnetic axis. We find that, for a narrow range of the fastness parameter, the Alfv\'en and the equilibrium surfaces intersect at two different critical latitudes. While the system is transiting from the propeller to the accretion regime (the initial rise of an outburst), accretion proceeds from the region above the higher critical latitude. In transitions from the accretion to the propeller regime (decay of the outburst), accretion of matter can proceed from the disc midplane. As a result, the transition from the propeller to the accretion regime occurs at a luminosity higher than the transition from the accretion to the propeller regime. We discuss the implications of our results for spectral transitions exhibited by low-mass X-ray binaries.

In this paper we investigate a Yukawa gravity modification of the Newtonian gravitational potential in a weak field approximation. For that purpose we derived the corresponding equations of motion and used them to perform two-body simulations of the stellar orbits. In 2020 the GRAVITY Collaboration detected the orbital precession of the S2 star around the supermassive black hole (SMBH) at the Galactic Center (GC) and showed that it is close to the general relativity (GR) prediction. Using this observational fact, we evaluated parameters of the Yukawa gravity (the range of Yukawa interaction $\Lambda$ and universal constant $\delta$) with the Schwarzschild precession of the S-stars assuming that the observed Schwarzschild precession will be equal to their GR estimates. GR provides the most natural way to fit observational data for S-star orbits, however, their precessions can be fitted by Yukawa gravity. Our main goal was to study the possible influence of the strength of Yukawa interaction, i.e. the universal constant $\delta$, on the precessions of S-star orbits. We analyze S-star orbits assuming different strength of Yukawa interaction $\delta$ and find that this parameter has strong influence on range of Yukawa interaction $\Lambda$. Using MCMC simulations we obtain the best-fit values and uncertanties of Yukawa gravity parameters for S-stars. Also, we introduce a new criterion which can be used for classification of gravitational systems in this type of gravity, according to their scales. We demonstrated that performed analysis of the observed S-stars orbits around the GC in the frame of the Yukawa gravity represent a tool for constraining the Yukawa gravity parameters and probing the predictions of gravity theories.

Aims. We validate a semi-analytical model for the covariance of real-space 2-point correlation function of galaxy clusters. Methods. Using 1000 PINOCCHIO light cones mimicking the expected Euclid sample of galaxy clusters, we calibrate a simple model to accurately describe the clustering covariance. Then, we use such a model to quantify the likelihood analysis response to variations of the covariance, and investigate the impact of a cosmology-dependent matrix at the level of statistics expected for the Euclid survey of galaxy clusters. Results. We find that a Gaussian model with Poissonian shot-noise does not correctly predict the covariance of the 2-point correlation function of galaxy clusters. By introducing few additional parameters fitted from simulations, the proposed model reproduces the numerical covariance with 10 per cent accuracy, with differences of about 5 per cent on the figure of merit of the cosmological parameters $\Omega_{\rm m}$ and $\sigma_8$. Also, we find that the cosmology-dependence of the covariance adds valuable information that is not contained in the mean value, significantly improving the constraining power of cluster clustering. Finally, we find that the cosmological figure of merit can be further improved by taking mass binning into account. Our results have significant implications for the derivation of cosmological constraints from the 2-point clustering statistics of the Euclid survey of galaxy clusters.

One of the leading explanations for the origin of Fermi Bubbles is a past jet activity in the Galactic center supermassive black hole Sgr A$^*$. The claimed jets are often assumed to be perpendicular to the Galactic plane. Motivated by the orientation of pc-scale nuclear stellar disk and gas streams, and a low inclination of the accretion disk around Sgr A$^*$ inferred by the Event Horizon Telescope, we perform hydrodynamical simulations of nuclear jets significantly tilted relative to the Galactic rotation axis. The observed axisymmetry and hemisymmetry (north-south symmetry) of Fermi/eROSITA bubbles (FEBs) due to quasi-steady jets in Sgr A$^*$ can be produced if the jet had a super-Eddington power ($\gtrsim 5\times 10^{44}$ erg s$^{-1}$) for a short time (jet active period $\lesssim 6$ kyr) for a reasonable jet opening angle ($\lesssim 10^\circ$). Such powerful explosions are, however, incompatible with the observed O VIII/O VII line ratio towards the bubbles, even after considering electron-proton temperature non-equilibrium. We argue that the only remaining options for producing FEBs are i) a low-luminosity ($\approx 10^{40.5-41}$ erg s$^{-1}$)) magnetically dominated jet or accretion wind from the Sgr A$^*$, and ii) a SNe or TDE driven wind of a similar luminosity from the Galactic center.

We study the formation of the first galaxies in overdense regions modelled by the FOREVER22 simulation project. Our simulations successfully reproduce the star formation rates and the $M_{\rm UV}-M_{\rm star}$ relations of candidate galaxies at $z \sim 10-14$ observed by James Webb Space Telescope (JWST). We suggest that the observed galaxies are hosted by dark-matter haloes with $M_{\rm h} \sim 10^{11}~\rm M_{\odot}$ and are in short-period starburst phases. On the other hand, even simulated massive galaxies in overdense regions cannot reproduce the intense star formation rates and the large stellar masses of observed candidates at $z \sim 17$. Also, we show that the contribution of population III stars to the UV flux decreases as the stellar mass increases and it is a few percent for galaxies with $M_{\rm star} \sim 10^{8}~\rm M_{\odot}$. Therefore, a part of the observed flux by JWST could be the light from population III stars. Our simulations suggest that the UV flux can be dominated by population III stars and the UV-slope shows $\beta \lesssim -3$ if future observations would reach galaxies with $M_{\rm stars} \sim 10^{5}~\rm M_{\odot}$ at $z \sim 20$ of which the mass fraction of population III stars can be greater than 10 percent.

Determining the atomic and molecular content of the interstellar medium (ISM) as a function of environmental parameters is of fundamental importance to understand the star-formation process across the epochs. Although there exist various three-dimensional hydro-chemical codes modelling the ISM at different scales and redshifts, they are computationally expensive and inefficient for studies over a large parameter space. Building on our earlier approach, we present PDFchem, a novel algorithm that models the cold ISM at moderate and large scales using functions connecting the quantities of the local ($A_{\rm V,eff}$) and the observed ($A_{\rm V,obs}$) visual extinctions, and the local number density, $n_{\rm H}$, with probability density functions (PDF) of $A_{\rm V,obs}$ on cloud scales typically tens-to-hundreds of pc as an input. For any given $A_{\rm V,obs}$-PDF, sampled with thousands of clouds, the algorithm instantly computes the average abundances of the most important species (HI, H$_2$, CII, CI, CO, OH, OH$^+$, H$_2$O$^+$, CH, HCO$^+$) and performs radiative transfer calculations to estimate the average emission of the most commonly observed lines ([CII]~$158\mu$m, both [CI] fine-structure lines and the first five rotational transitions of $^{12}$CO). We examine two $A_{\rm V,obs}$-PDFs corresponding to a non star-forming and a star-forming ISM region, under a variety of environmental parameters combinations. These cover FUV intensities in the range of $\chi/\chi_0=10^{-1}-10^3$, cosmic-ray ionization rates in the range of $\zeta_{\rm CR}=10^{-17}-10^{-13}\,{\rm s}^{-1}$ and metallicities in the range of $Z=0.1-2\,{\rm Z}_{\odot}$. PDFchem is fast, easy to use, reproduces the PDR quantities of the time-consuming hydrodynamical models and can be used directly with observed data to understand the evolution of the cold ISM chemistry.

PSR J1023+0038 is the first millisecond pulsar that was ever observed as an optical and UV pulsar. So far, it is the only optical transitional millisecond pulsar. The rotation- and accretion-powered emission mechanisms hardly individually explain the observed characteristics of optical pulsations. A synergistic model, combining these standard emission processes, was proposed to explain the origin of the X-ray/UV/optical pulsations. We study the phase lag between the pulses in the optical and X-ray bands to gain insight into the physical mechanisms that cause it. We performed a detailed timing analysis of simultaneous or quasi-simultaneous observations in the X-ray band, acquired with the XMM-Newton and NICER satellites, and in the optical band, with the fast photometers SiFAP2 (mounted at the 3.6 m Telescopio Nazionale Galileo) and Aqueye+ (mounted at the 1.8 m Copernicus Telescope). We estimated the time lag of the optical pulsation with respect to that in the X-rays by modeling the folded pulse profiles with two harmonic components. Optical pulses lag the X-ray pulses by $\sim$ 150 $\mu$s in observations acquired with instruments (NICER and Aqueye+) whose absolute timing uncertainty is much smaller than the measured lag. We also show that the phase lag between optical and X-ray pulsations lies in a limited range of values, $\delta \phi \in$ (0 $-$ 0.15), which is maintained over timescales of about five years. This indicates that both pulsations originate from the same region, and it supports the hypothesis of a common emission mechanism. Our results are interpreted in the shock-driven mini pulsar nebula scenario. This scenario suggests that optical and X-ray pulses are produced by synchrotron emission from the shock that formed within a few light cylinder radii away ($\sim$ 100 km) from the pulsar, where its striped wind encounters the accretion disk inflow.

The determination of the chemical composition of interplanetary coronal mass ejection (ICME) plasma is an open issue. More specifically, it is not yet fully understood how remote sensing observations of the solar corona plasma during solar disturbances evolve into plasma properties measured in situ away from the Sun. The ambient conditions of the background interplanetary plasma are important for space weather because they influence the evolutions, arrival times, and geo-effectiveness of the disturbances. The Reverse In situ and MHD APproach (RIMAP) is a technique to reconstruct the heliosphere on the ecliptic plane (including the magnetic Parker spiral) directly from in situ measurements acquired at 1 AU. It combines analytical and numerical approaches, preserving the small-scale longitudinal variability of the wind flow lines. In this work, we use RIMAP to test the interaction of an ICME with the interplanetary medium. We model the propagation of a homogeneous non-magnetised (i.e. with no internal flux rope) cloud starting at 800 km s-1 at 0.1 AU out to 1.1 AU. Our 3D magnetohydrodynamics (MHD) simulation made with the PLUTO MHD code shows the formation of a compression front ahead of the ICME, continuously driven by the cloud expansion. Using a passive tracer, we find that the initial ICME material does not fragment behind the front during its propagation, and we quantify the mixing of the propagating plasma cloud with the ambient solar wind plasma, which can be detected at 1 AU.

As a key feature, NASA's Parker Solar Probe (PSP) and ESA-NASA's Solar Orbiter (SO) missions cooperate to trace solar wind and transients from their sources on the Sun to the inner interplanetary space. The goal of this work is to accurately reconstruct the interplanetary Parker spiral and the connection between coronal features observed remotely by the Metis coronagraph on-board SO and those detected in situ by PSP at the time of the first PSP-SO quadrature of January 2021. We use the Reverse In-situ and MHD Approach (RIMAP), a hybrid analytical-numerical method performing data-driven reconstructions of the Parker spiral. RIMAP solves the MHD equations on the equatorial plane with the PLUTO code, using the measurements collected by PSP between 0.1 and 0.2 AU as boundary conditions. Our reconstruction connects density and wind speed measurements provided by Metis (3-6 solar radii) to those acquired by PSP (21.5 solar radii) along a single streamline. The capability of our MHD model to connect the inner corona observed by Metis and the super Alfv\'enic wind measured by PSP, not only confirms the research pathways provided by multi-spacecraft observations, but also the validity and accuracy of RIMAP reconstructions as a possible test bench to verify models of transient phenomena propagating across the heliosphere, such as coronal mass ejections, solar energetic particles and solar wind switchbacks.

The study of post asymptotic giant branch (post-AGB) stars is a valuable tool to study still poorly known aspects of the evolution of the stars through the AGB. This is due to the accurate determination of their surface chemical composition and to the peculiar shape of the SED: the emission from the central star can be easily disentangled from the contribution from the dusty shell, which can then be characterized. The goal of the present study is to reconstruct the dust formation process and more generally the late phases of the evolution of O-rich stars across the AGB phase. This is performed by studying O-rich, post-AGB stars, which are analyzed in terms of their luminosity, Teff and infrared excess. We study sources classified as single, O-rich, post-AGB stars in the Galaxy, which exhibit a double-peaked (shell-type) SED. We use results from stellar evolution modelling combined with dust formation and radiative transfer modelling to reconstruct the late AGB phases and the initial contraction to the post-AGB phase. We also determine the mass-loss and dust formation rates for stars of different mass and chemical composition. The analysis of the IR excess of the post-AGB, O-rich stars examined in this study outlines an interesting complexity, in terms of the correlation between the dust in the surroundings of the stars, the evolutionary status and the progenitor's mass. The sources descending from massive AGBs (>3Msun depending on metallicity) are generally characterized by higher IR excess than the lower mass counterparts, owing to the more intense dust formation taking place during the final AGB phases. From the determination of the location of the dusty regions we deduce that the expanding velocities of the outflow change significantly from star to star. The possibility that radiation pressure is unable to accelerate the wind in the faintest objects is also commented.

The far-IR/sub-mm wavelength range contains a wealth of diagnostic information that is important for understanding the role of radio AGN in galaxy evolution. Here we present the results of Herschel PACS and SPIRE observations of a complete sample of 46 powerful 2Jy radio AGN at intermediate redshifts (0.05 < z < 0.7), which represent the deepest pointed observations of a major sample of radio AGN undertaken by Herschel. In order to assess the importance of non-thermal synchrotron emission at far-IR wavelengths, we also present new APEX sub-mm and ALMA mm data. We find that the overall incidence of non-thermal contamination in the PACS bands ($<$200$\mu$m) is in the range 28 -- 43%; however, this rises to 30 -- 72% for wavelengths ($> $200$\mu$m) sampled by the SPIRE instrument. Non-thermal contamination is strongest in objects with compact CSS/GPS or extended FRI radio morphologies, and in those with type 1 optical spectra. Considering thermal dust emission, we find strong correlations between the 100 and 160$\mu$m monochromatic luminosities and AGN power indicators, providing further evidence that radiation from the AGN may be an important heating source for the far-IR emitting dust. Clearly, AGN contamination -- whether by the direct emission from synchrotron-emitting lobes and cores, or via radiative heating of the cool dust -- needs to be carefully considered when using the far-IR continuum to measure the star formation rates in the host galaxies of radio AGN.

We study the utility of the [OI] 6300$\mathring{\mathrm A}$ forbidden line for identifying and interpreting externally driven photoevaporative winds in different environments and at a range of distances. Thermally excited [OI] 6300$\mathring{\mathrm A}$ is a well known tracer of inner disc winds, so any external contribution needs to be distinguishable. In external winds, the line is not thermally excited and instead results from the dissociation of OH and we study how the line luminosity resulting from that process scales with the disc/environmental parameters. We find that the line luminosity increases dramatically with FUV radiation field strength above around 5000 G$_0$. The predicted luminosities from our models are consistent with measurements of the line luminosity of proplyds in the Orion Nebula Cluster. The high luminosity in strong UV environments alone may act as a diagnostic, but a rise in the [OI]-to-accretion luminosity ratio is predicted to better separate the two contributions. This could provide a means of identifying external photoevaporation in distant clusters where the proplyd morphology of evaporating discs cannot be spatially resolved.

With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We present the analysis of K-band spectro-interferometric observations toward the Massive Young Stellar Object IRAS 13481-6124, which is known to host an accreting dusty disc. Using GRAVITY on the VLTI, we trace the crucial au-scales of the warm inner interface between the star and the accretion dusty disc. We detect and spatially resolve the Na I doublet and He I transitions towards an object of this class for the first time. The new observations in combination with our geometric models allowed us to probe the smallest -au- scales of accretion/ejection around an MYSO. We find that Na I originates in the disc at smaller radii than the dust disc and is more compact than any of the other spatially resolved diagnostics (Br$\gamma$, He I, and CO). Our findings suggest that Na I can be a new powerful diagnostic line in tracing the warm star/disc accreting interface of forming (massive) stars, while the similarities between He I and Br$\gamma$ point towards an accretion/ejection origin of He I

To probe the current theory on galaxy formation and evolution, an increased synergy between observations and simulations is necessary. For this reason, in our previous paper of this series, we presented a method to mock SDSS-IV/MaNGA integral-field spectroscopic galaxy observations from cosmological simulations of galaxy formation. Here we present the resulting mock galaxy catalogue. This catalogue consists of 1,000 unique galaxies in IllustrisTNG-50 falling into the SDSS-IV/MaNGA-primary target footprint, defined in the redshift and i-band absolute magnitude space. In this paper, we describe the general characteristics of the catalogue, in terms of morphology, kinematics, and stellar population properties. We also investigate our ability to recover the galaxy characteristics, as given by the simulations, analysing the synthetic spectra. We show that the `intrinsic' and recovered kinematics, age and metallicity are consistent within 1${\sigma}$, with residuals over all tassels ($\sim 8$ million) consistent with $0$ at the $68\%$ confidence level. We also compare `intrinsic' and recovered star formation histories, finding a close resemblance. Therefore, our mocking and spectral fitting processes do not distort intrinsic galaxy properties, hence we can use these results for scientific analysis. In the next papers of this series, we will present a comprehensive comparison and scientific analysis of TNG50 simulations with MaNGA observational results.

We report discoveries of 165 new quasar Ca II absorbers from the Sloan Digital Sky Survey (SDSS) Data Release 7 and 12. Our Ca II rest frame equivalent width distribution supports the weak and strong subpopulations, split at ${W}^{\lambda3934}_{0}=0.7${\AA}. Comparison of both populations' dust depletion shows clear consistency for weak absorber association with halo-type gas in the Milky Way (MW) while strong absorbers have environments consistent with halo and disc-type gas. We probed our high redshift Ca II absorbers for 2175{\AA} dust bumps, discovering 12 2175{\AA} dust absorbers (2DAs). This clearly shows that some Ca II absorbers follow the Large Magellanic Cloud (LMC) extinction law rather than the Small Magellanic Cloud extinction law. About 33% of our strong Ca II absorbers exhibit the 2175{\AA} dust bump while only 6% of weak Ca II absorbers show this bump. 2DA detection further supports the theory that strong Ca II absorbers are associated with disk components and are dustier than the weak population. Comparing average Ca II absorber dust depletion patterns to that of Damped Ly{\alpha} Absorbers (DLAs), Mg II absorbers, and 2DAs shows that Ca II absorbers generally have environments with more dust than DLAs and Mg II absorbers, but less dust than 2DAs. Comparing 2175{\AA} dust bump strengths from different samples and also the MW and LMC, the bump strength appears to grow stronger as the redshift decreases, indicating dust growth and the global chemical enrichment of galaxies in the universe over time.

Identification of Cherenkov light generated by muons has been suggested as a promising way to dramatically improve the background rejection power of Imaging Atmospheric Cherenkov Telescope (IACT) arrays at high energies. However, muon identification remains a challenging task, for which efficient algorithms are still being developed. We present an approach in which, rather than identifying Cherenkov light from muons, we simply consider the presence of Cherenkov light other than the main shower image in IACTs with large mirror area. We show that in the case of the H.E.S.S. array of five telescopes this approach results in background rejection improvements at all energies above 1 TeV. In particular, the rejection power can be improved by a factor $\sim3-4$ at energies above 20 TeV while keeping $\sim90\%$ of the original gamma-ray efficiency.

Local low metallicity dwarf galaxies are relics of the early universe and hold clues into the origins of supermassive black holes (SMBHs). In recent work, coronal lines have been used to unveil a population of candidate accreting black holes in dwarf galaxies with gas phase metallicities and stellar masses well below the host galaxies of any previously known AGNs. Using MUSE/VLT observations, we report the detection of [Fe X] $\lambda$6374 coronal line emission and a broad H$\alpha$ line in the nucleus of SDSS J094401.87$-$003832.1, a nearby ($z=0.0049$) metal poor dwarf galaxy at least fifty times less massive than the LMC. The [Fe X] $\lambda$6374 emission is compact and centered on the brightest nuclear source, with a spatial extent of $\approx$100 pc. The [Fe X] luminosity is $\approx 10^{37}$ erg s$^{-1}$, within the range seen in previously identified AGNs in the dwarf galaxy population. This line has never been observed in gas ionized by hot stars. While it can be produced in supernova ejecta, the [Fe X] flux from SDSS J094401.87$-$003832.1 has persisted over the ~19 year time period between the SDSS and MUSE observations, ruling out supernovae as the origin for the emission. The FWHM of the broad component of the H$\alpha$ line is $446 \pm 17$ km s$^{-1}$ and its luminosity is $\approx 1.5\times10^{38}$ erg s$^{-1}$, lower than the broad line luminosities of previously identified low mass broad line AGNs. These observations, together with previously reported multi-wavelength observations, can most plausibly be explained by the presence of an accreting intermediate mass black hole in a primordial galaxy analog. However, we cannot rule out the possibility that current stellar population models of metal poor stars significantly under-predict the stellar ionizing photon flux, and that metal poor stars can produce an extreme ionizing spectrum similar to that produced by AGNs.

Symbiotic binaries sometimes hide their symbiotic nature for significant periods of time. There is mounting observational evidence that in those symbiotics that are powered solely by accretion of red-giant's wind material onto a white dwarf, without any quasi-steady shell burning on the surface of the white dwarf, the characteristic emission lines in the optical spectrum can vanish, leaving the semblance of an isolated red giant spectrum. Here we present compelling evidence that this disappearance of optical emission lines from the spectrum of RT Cru during 2019 was due to a decrease in the accretion rate, which we derive by modeling the X-ray spectrum. This drop in accretion rate leads to a lower flux of ionizing photons and thus to faint/absent photoionization emission lines in the optical spectrum. We observed the white dwarf symbiotic RT Cru with XMM-Newton and Swift in X-rays and UV and collected ground-based optical spectra and photometry over the last 33 years. This long-term coverage shows that during most of the year 2019, the accretion rate onto the white dwarf was so low, $\dot{M}= (3.2\pm 0.06)\, \times$10$^{-11}$ $M_{\odot}$ yr$^{-1}$ (d/2.52 kpc)$^2$, that the historically detected hard X-ray emission almost vanished, the UV flux faded by roughly 5 magnitudes, the $U$, $B$ and $V$ flickering amplitude decreased, and the Balmer lines virtually disappeared from January through March 2019. Long-lasting low-accretion episodes as the one reported here may hamper the chances of RT Cru experiencing nova-type outburst despite the high-mass of the accreting white dwarf.

The Observatorio Astrof\'isico de Javalambre (OAJ{\dag}1) in Spain is a young astronomical facility, conceived and developed from the beginning as a fully automated observatory with the main goal of optimizing the processes in the scientific and general operation of the Observatory. The OAJ has been particularly conceived for carrying out large sky surveys with two unprecedented telescopes of unusually large fields of view (FoV): the JST/T250, a 2.55m telescope of 3deg field of view, and the JAST/T80, an 83cm telescope of 2deg field of view. The most immediate objective of the two telescopes for the next years is carrying out two unique photometric surveys of several thousands square degrees, J-PAS{\dag}2 and J-PLUS{\dag}3, each of them with a wide range of scientific applications, like e.g. large structure cosmology and Dark Energy, galaxy evolution, supernovae, Milky Way structure, exoplanets, among many others. To do that, JST and JAST are equipped with panoramic cameras under development within the J-PAS collaboration, JPCam and T80Cam respectively, which make use of large format (~ 10k x 10k) CCDs covering the entire focal plane. This paper describes in detail, from operations point of view, a comparison between the detailed cost of the global automation of the Observatory and the standard automation cost for astronomical facilities, in reference to the total investment and highlighting all benefits obtained from this approach and difficulties encountered. The paper also describes the engineering development of the overall facilities and infrastructures for the fully automated observatory and a global overview of current status, pinpointing lessons learned in order to boost observatory operations performance, achieving scientific targets, maintaining quality requirements, but also minimizing operation cost and human resources.

Generating axion dark matter through the kinetic misalignment mechanism implies the generation of large asymmetries for Standard Model fermions in the early universe. Even if these asymmetries are washed out at later times, they can trigger a chiral plasma instability in the early universe. Similarly, a direct coupling of the axion with the hypercharge gauge field can trigger a tachyonic instability. These instabilities produce helical magnetic fields, which are preserved until the electroweak phase transition. At the electroweak phase transition, these become a source of baryon asymmetry, which can be much more efficient than the original axiogenesis proposal. We discuss constraints on axion dark matter production from the overproduction of the baryon asymmetry as well as a minimal, albeit fine-tuned setup, where both the correct dark matter abundance and baryon asymmetry can be achieved. For a given axion decay constant, this leads to a sharp prediction for the mass of the radial direction of the Peccei Quinn field, which is a soft mass scale in supersymmetric theories.

We present a systematic investigation of the possible locations for the special point (SP), a unique feature of hybrid neutron stars in the mass-radius diagram. The study is performed within the two-phase approach where the high-density (quark matter) phase is described by the covariant nonlocal Nambu--Jona-Lasinio (nlNJL) model equation of state (EOS) which is shown to be equivalent to a constant-sound-speed (CSS) EOS. For the nuclear matter phase around saturation density different relativistic density functional EOSs are used: DD2p00, its excluded-volume modification DD2p40 and the hypernuclear EOS DD2Y-T. In the present contribution we apply the Maxwell construction scheme for the deconfinement transition and demonstrate that a simultaneous variation of the vector and diquark coupling constants results in the occurrence of SP "trains" which are invariant against changing the nuclear matter EOS. We propose that the SP train corresponding to a variation of the diquark coupling at constant vector coupling is special since it serves as a lower bound for the line of maximum masses and accessible radii of massive hybrid stars.

The hypertriton is predicted to have a small binding energy (a weighted average of about 170 keV), consistent with a large matter radius (~ 10 fm), large than the historical11Li halo discovered more than 35 years ago. But the reported experimental values of the binding energy of the hypertriton range from 50 to 500 keV. In this work I discuss the electromagnetic response and interaction radius of the hypertriton and how high energy heavy ion collisions (~ 1 - 2 GeV/nucleon) can help achieving a higher accuracy for the determination of its size and binding energy.

We present a relativistic density functional approach to color superconducting quark matter that mimics quark confinement by a fast growth of the quasiparticle selfenergy in the confining region. The approach is shown to be equivalent to a chiral model of quark matter with medium dependent couplings. While the (pseudo)scalar sector of the model is fitted to the vacuum phenomenology of quantum chromodynamics, the strength of interaction in the vector and diquark channels is varied in order to provide the best agreement with the observational constraints on the mass-radius relation and tidal deformability of neutron stars modelled with our approach. In order to recover the conformal behavior of quark matter at asymptotically high densities we introduce a medium dependence of the vector and diquark couplings motivated by the nonperturbative gluon exchange. Our analysis signals that the onset of deconfinement to color superconducting quark matter is likely to occur in neutron stars with masses below 1.0 $M_\odot$.

Shockwaves in plasma are usually dealt with using Magnetohydrodynamics (MHD). Yet, MHD entails the assumption of a short mean free path, which is not fulfilled in a collisionless plasma. Recently, for pair plasmas, we devised a model allowing to account for kinetic effects within an MHD-like formalism. Its relies on an estimate of the anisotropy generated when crossing the front, with a subsequent assessment of the stability of this anisotropy in the downstream. We solved our model for parallel, perpendicular and switch-on shocks. Here we bridge between all these cases by treating the problem of an arbitrarily, but coplanar, oriented magnetic field. Even though the formalism presented is valid for anisotropic upstream temperatures, only the case of a cold upstream is solved. We find extra solutions which are not part of the MHD catalog, and a density jump that is notably less in the quasi parallel, highly magnetized, regime. Given the complexity of the calculations, this work is mainly devoted to the presentation of the mathematical aspect of our model. A forthcoming article will be devoted to the physics of the shocks here defined.

The rate of energy loss and orbital period decay of quasi-stable compact binary systems are derived in $f(R)$ theory of gravity using the method of a single vertex graviton emission process from a classical source. After linearising the $f(R)$ action written in an equivalent scalar-tensor format in the Einstein frame, we identify the appropriate interaction terms between the massless spin-2 tensor mode, massive scalar mode, and the energy momentum tensor. The definition of the scalar field is related to the $f(R)$ models. Then using the interaction vertex we compute the rate of energy loss due to spin-2 quadrupole radiation, which comes out to be the same as the Peter-Mathews formula with a multiplication factor, and also the energy loss due to the scalar dipole radiation. The total energy loss is the sum of these two contributions. Our derivation is most general as it is applicable for both arbitrary eccentricity of the binary orbits and arbitrary mass of the scalar field. Using the derived theoretical formula for the period decay of the binary systems, we compare the predictions of $f(R)$ gravity and general relativity for the observations of three binary systems, i.e. Hulse-Taylor Binary, PSR J1141-6545 and PSR J1738+0333. Thus we put bound on three well-known $f(R)$ dark energy models, namely the Hu-Sawicki, the Straobinsky, and the Tsujukawa model. We get the best constraint on $f'(R_0)-1$ (where $R_0$ is the scalar curvature of the Universe at the present epoch) from the Tsujikawa model, i.e $\vert f'(R_0)-1\vert < 3.44\times 10^{-4}$. This bound is stronger than those from most of the astrophysical observations and even some cosmological observations.

In a companion article, we discussed the radiometric sensitivity and resolution of a new passive optical sensing technique, Space-Time Projection Optical Tomography (SPOT), to detect and track sub-cm and larger space debris for Space Situational Awareness. SPOT is based on the principle that long synthetic exposure can be achieved if the phase-space trajectory of a hypothetical point-source is precisely predictable within a very wide telescope field-of-view, which is the case for orbiting debris. This article discusses the computational search space for debris mining as well as a recursive measure-and-fit algorithm based on a generalized Hough transform for orbit determination.

We propose an alternative scenario of axion misalignment mechanism based on nontrivial interplay between axion and a light dilaton in the early universe. Dark matter abundance is still sourced by the initial misalignment of axion field, whose motion along the potential kicks the dilaton field away from its minimum, and dilaton starts to oscillate later with a delayed onset time for oscillation and a relatively large misalignment value due to kick, eventually the dilaton dominates over axion in their energy densities and dilaton is identified as dark matter. The kick effect due to axion motion is the most significant if the initial field value of dilaton is near its minimum, therefore we call this scenario axion free-kick misalignment mechanism, where axion plays the role similar to a football player. Dark matter abundance can be obtained with a lower axion decay constant compared to the conventional misalignment mechanism.

The current LIGO-Virgo observing run has been pushing the sensitivity limit to touch the stochastic gravitational-wave backgrounds (SGWBs). However, no significant detection has been reported to date for any single dominated source of SGWBs with a single broken-power-law (BPL) spectrum. Nevertheless, it could equally well escape from existing Bayesian searches from, for example, two comparable dominated sources with two separate BPL spectra (double-peak case) or a single source with a doubly-BPL (DBPL) spectrum (doubly-broken case). In this paper, we put constraints on these two cases from Advanced LIGO-Virgo's first three observing runs. We found strong negative evidence for the double-peak case and hence place 95\% confidence-level (CL) upper limits $\Omega_\mathrm{BPL,1}<2.5\times10^{-7}$ and $\Omega_\mathrm{BPL,2}<9.4\times10^{-8}$ on the two BPL spectra amplitudes with respect to the unresolved compact-binary-coalescence (CBC) amplitude $\Omega_\mathrm{CBC}<5.6\times10^{-9}$. We further found weak negative evidence for the doubly-broken case and hence place 95\% CL upper limit $\Omega_\mathrm{DBPL}<1.7\times10^{-7}$ on the overall amplitude of the DBPL spectrum with respect to $\Omega_\mathrm{CBC}<6.0\times10^{-9}$. The implications of cosmological first-order phase transitions are also discussed.

We study the possibility of generating light Dirac neutrino mass from a radiative seesaw mechanism with dark sector particles going inside the loop, known as the scotogenic framework. The loop suppression and additional free parameters allow large ($\sim\mathcal{O}(1))$ coupling of light Dirac neutrinos with the dark sector particles. Such large Yukawa coupling not only dictates the relic abundance of heavy fermion singlet dark matter but also can lead to the thermalisation of the right chiral part of Dirac neutrinos, generating additional relativistic degrees of freedom ${\rm \Delta{N_{eff}}}$. We find that the parameter space consistent with dark matter phenomenology and neutrino mass bounds can also be probed at future cosmic microwave background experiments like CMB-S4 via precision measurements of ${\rm \Delta{N_{eff}}}$. The same parameter space, while leading to loop-suppressed direct detection cross-section outside future sensitivities, can also have other interesting and complementary observational prospects at colliders, charged lepton flavour violation.

This Multi-Petawatt Physics Prioritization (MP3) Workshop Report captures the outcomes from a community-initiated workshop held April 20-22, 2022 at Sorbonne University in Paris, France. The MP3 workshop aimed at developing science questions to guide research and future experiments in four areas identified by corresponding MP3 working groups: high-field physics and quantum electrodynamics (HFP/QED), laboratory astrophysics and planetary physics (LAPP), laser-driven nuclear physics (LDNP), and particle acceleration and advanced light sources (PAALS).

Asymmetric reheating is a generic requirement for models of dark sectors with light species, but its implementation is usually in tension with unique phenomenologies otherwise possible in compelling theories containing dark copies of the Standard Model. We present a simple module to implement asymmetric reheating during a $\mathbb{Z}_2$-breaking phase above some critical temperature. This reinvigorates the possibility of an exactly degenerate mirror sector and the striking phenomenology of composite particles oscillating into their mirror counterparts.