Papers for Friday, Dec 08 2023

There are more than 200 moons in our Solar System, but their relatively small radii make similarly sized extrasolar moons very hard to detect with current instruments. The best exomoon candidates so far are two nearly Neptune-sized bodies orbiting the Jupiter-sized transiting exoplanets Kepler-1625 b and Kepler-1708 b, but their existence has been contested. Here we reanalyse the Hubble and Kepler data used to identify the two exomoon candidates employing nested sampling and Bayesian inference techniques coupled with a fully automated photodynamical transit model. We find that the evidence for the Kepler-1625 b exomoon candidate comes almost entirely from the shallowness of one transit observed with Hubble. We interpret this as a fitting artifact in which a moon transit is used to compensate for the unconstrained stellar limb darkening. We also find much lower statistical evidence for the exomoon candidate around Kepler-1708 b than previously reported. We suggest that visual evidence of the claimed exomoon transits is corrupted by stellar activity in the Kepler light curve. Our injection-retrieval experiments of simulated transits in the original Kepler data reveal false positive rates of 10.9% and 1.6% for Kepler-1625 b and Kepler-1708 b, respectively. Moreover, genuine transit signals of large exomoons would tend to exhibit much higher Bayesian evidence than these two claims. We conclude that neither Kepler-1625 b nor Kepler-1708 b are likely to be orbited by a large exomoon.

Despite the promise of Neural Posterior Estimation (NPE) methods in astronomy, the adaptation of NPE into the routine inference workflow has been slow. We identify three critical issues: the need for custom featurizer networks tailored to the observed data, the inference inexactness, and the under-specification of physical forward models. To address the first two issues, we introduce a new framework and open-source software nbi (Neural Bayesian Inference), which supports both amortized and sequential NPE. First, nbi provides built-in "featurizer" networks with demonstrated efficacy on sequential data, such as light curve and spectra, thus obviating the need for this customization on the user end. Second, we introduce a modified algorithm SNPE-IS, which facilities asymptotically exact inference by using the surrogate posterior under NPE only as a proposal distribution for importance sampling. These features allow nbi to be applied off-the-shelf to astronomical inference problems involving light curves and spectra. We discuss how nbi may serve as an effective alternative to existing methods such as Nested Sampling. Our package is at https://github.com/kmzzhang/nbi.

The total disk mass sets the formation potential for exoplanets. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing for the possibility of using CO as a more effective total gas mass tracer. Using ALMA archival data and new NOEMA data of 12CO, 13CO, and C18O transitions of 35 Herbig disks within 450 pc, the masses are determined using the thermo-chemical code Dust And LInes (DALI). The majority of Herbig disks for which 13CO and C18O are detected are optically thick in both. Computing the gas mass using a simple optically thin relation between line flux and column density results in an underestimate of the gas mass of at least an order of magnitude compared to the masses obtained with DALI. The inferred gas masses with DALI are consistent with a gas-to-dust ratio of at least 100. These gas-to-dust ratios are two orders of magnitude higher compared to those found for T Tauri disks using similar techniques, even over multiple orders of magnitude in dust mass, illustrating the importance of chemical conversion of CO in colder T Tauri disks. Similar high gas-to-dust ratios are found for Herbig group I and II disks. Since group II disks have dust masses comparable to T Tauri disks, their higher CO gas masses illustrate the determining role of temperature. Compared to debris disks, Herbig disks have four orders of magnitude higher gas masses. At least one Herbig disk, HD 163296, has a detected molecular disk wind, but our investigation has not turned up other detections of the CO disk wind in spite of similar sensitivities.

The flavor evolution of neutrinos in dense astrophysical sources, such as core-collapse supernovae or compact binary mergers, is non-linear due to the coherent forward scattering of neutrinos among themselves. Recent work in this context has been addressed to figure out whether flavor equipartition could be a generic flavor outcome of fast flavor conversion. We investigate the flavor conversion physics injecting random perturbations in the neutrino field in two simulation setups: 1. a spherically symmetric simulation shell without periodic boundaries, with angular distributions evolving dynamically thanks to non-forward scatterings of neutrinos with the background medium, and neutrino advection; 2. a periodic simulation shell, with angular distributions of neutrinos defined a priori and neutrino advection. We find that, independent of the exact initial flavor configuration and type of perturbations, flavor equipartition is generally achieved in the system with periodic boundaries; in this case, perturbations aid the diffusion of flavor structures to smaller and smaller scales. However, flavor equipartition is not a general outcome in the simulation shell without periodic boundaries, where the inhomogeneities induced perturbing the neutrino field affect the flavor evolution, but do not facilitate the diffusion of flavor waves. This work highlights the importance of the choice of the simulation boundary conditions in the exploration of fast flavor conversion physics.

AT 2019azh is a H+He tidal disruption event (TDE) with one of the most extensive ultraviolet and optical datasets available to date. We present our photometric and spectroscopic observations of this event starting several weeks before and out to approximately two years after g-band peak brightness and combine them with public photometric data. This extensive dataset robustly reveals a change in the light-curve slope and a bump in the rising light curve of a TDE for the first time, which may indicate more than one dominant emission mechanism contributing to the pre-peak light curve. We further confirm the relation seen in previous TDEs whereby the redder emission peaks later than the bluer emission. The post-peak bolometric light curve of AT 2019azh is better described by an exponential decline than by the canonical t^{-5/3} (and in fact any) power-law decline. We find a possible mid-infrared excess around peak optical luminosity, but cannot determine its origin. In addition, we provide the earliest measurements of the Halpha emission-line evolution and find no significant time delay between the peak of the V-band light curve and that of the Halpha luminosity. These results can be used to constrain future models of TDE line formation and emission mechanisms in general. More pre-peak 1-2 day cadence observations of TDEs are required to determine whether the characteristics observed here are common among TDEs. More importantly, detailed emission models are needed to fully exploit such observations for understanding the emission physics of TDEs.

This article reports the first detection of a periodic light curve whose modulation is unambiguously due to rotation in a polluted white dwarf. Recent TESS Sector 68 time-series photometry of WD 2138-332, at 16.1 pc distance, exhibits a single periodogram peak with amplitude 0.39 per cent at frequency 3.876 1/d, corresponding to a period of 6.19 h. While this rotation is relatively rapid for isolated white dwarfs, it falls within the range of spin periods commonly found among those with detectable magnetic fields, where WD 2138-332 is notably both metal-rich and weakly magnetic with a 50 kG field. For the other 16 polluted white dwarfs in the local 20 pc volume -- five of which also exhibit magnetism -- multi-sector TESS data reveal no significant periodicities. In contrast, among the 23 magnetic and metal-free white dwarfs within 20 pc, six have light curves consistent with rotation periods between 0.7 and 35 h (three of which are new discoveries). These findings indicate the variable light curve of WD 2138-332 is almost certainly due to magnetism, as opposed to an inhomogeneous distribution of metals. From 13 metallic and magnetic white dwarfs with acceptable TESS data, a single detection of periodicity suggests that polluted white dwarfs are not rotating as rapidly as their magnetic counterparts. Planet ingestion is thus unlikely to be a significant channel for rapid rotation in white dwarfs, but may contribute in exceptional cases such as WD 2138-332.

Utilizing Planck polarized dust emission maps at 353 GHz and new large-area maps of the neutral hydrogen (HI) cold neutral medium (CNM) fraction ($f_\mathrm{CNM}$), we investigate the relationship between dust polarization fraction ($p_{353}$) and $f_\mathrm{CNM}$ in the diffuse high latitude ($|b|>30^{\circ}$) sky. We find that the correlation between $p_{353}$ and $f_\mathrm{CNM}$ is qualitatively distinct from the $p_{353}$-HI column density ($N_\mathrm{H\,I}$) relationship. At low column densities ($N_\mathrm{H\,I}<4\times10^{20}~\mathrm{cm^{-2}}$) where $p_{353}$ and $N_\mathrm{H\,I}$ are uncorrelated, there is a strong positive $p_{353}$-$f_\mathrm{CNM}$ correlation. We fit the $p_{353}$-$f_{\rm CNM}$ correlation with data-driven models to constrain the degree of magnetic field disorder between phases along the line-of-sight. We argue that an increased magnetic field disorder in the warm neutral medium (WNM) relative to the CNM best explains the positive $p_{353}$-$f_\mathrm{CNM}$ correlation in diffuse regions. Modeling the CNM-associated dust column as being maximally polarized, with a polarization fraction $p_{\rm CNM} \sim$ 0.2, we find that the best-fit mean polarization fraction in the WNM-associated dust column is 0.22$p_{\rm CNM}$. The model further suggests that a significant $f_{\rm CNM}$-correlated fraction of the non-CNM column (an additional ~18.4\% of the HI mass on average) is also more magnetically ordered, and we speculate that the additional column is associated with the unstable medium (UNM). Our results constitute the first large-area constraint on the average relative disorder of magnetic fields between the neutral phases of the ISM, and are consistent with the physical picture of a more magnetically aligned CNM column forming out of a disordered WNM.

The oldest stars in the Milky Way (born in the first few billion years) are expected to have a high density in the inner few kpc, spatially overlapping with the Galactic bulge. We use spectroscopic data from the Pristine Inner Galaxy Survey (PIGS) to study the dynamical properties of ancient, metal-poor inner Galaxy stars. We compute distances using StarHorse, and orbital properties in a barred Galactic potential. With this paper, we release the spectroscopic AAT/PIGS catalogue (13 235 stars). We find that most PIGS stars have orbits typical for a pressure-supported population. The fraction of stars confined to the inner Galaxy decreases with decreasing metallicity, but many very metal-poor stars (VMP, [Fe/H] < -2.0) stay confined (~ 60% stay within 5 kpc). The azimuthal velocity v$_\phi$ also decreases between [Fe/H] = -1.0 and -2.0, but is constant for VMP stars (at ~ 40 km/s). The carbon-enhanced metal-poor (CEMP) stars in PIGS appear to have similar orbital properties compared to normal VMP stars. Our results suggest a possible transition between two spheroidal components - a more metal-rich, more concentrated, faster rotating component, and a more metal-poor, more extended and slower/non-rotating component. We propose that the former may be connected to pre-disc in-situ stars (or those born in large building blocks), whereas the latter may be dominated by contributions from smaller galaxies. This is an exciting era where large metal-poor samples, such as in this work (as well as upcoming surveys, e.g., 4MOST), shed light on the earliest evolution of our Galaxy.

We use numerical simulations to model Gaia DR3 data with the aim of constraining the Milky Way bar and spiral structure parameters. We show that both the morphology and the velocity field in Milky Way-like galactic disc models are strong functions of time, changing dramatically over a few tens of Myr. This suggests that by finding a good match to the observed radial velocity field, v_R(x,y), we can constrain the bar-spiral orientation. Incorporating uncertainties into our models is necessary to match the data; most importantly, a heliocentric distance uncertainty above 10-15% distorts the bar's shape and $v_R$ quadrupole pattern morphology, and decreases its apparent angle with respect to the Sun-Galactocentric line. An excellent match to the \Gaia DR3 $v_R(x,y)$ field is found for a simulation with a bar length R_b~3.6 kpc. We argue that the data are consistent with a MW bar as short as ~3 kpc, for moderate strength inner disc spiral structure (A_2/A_0~0.25) or, alternatively, with a bar length up to ~5.2 kpc, provided that spiral arms are quite weak (A_2/A_0~0.1), and is most likely in the process of disconnecting from a spiral arm. We demonstrate that the bar angle and distance uncertainty can similarly affect the match between our models and the data - a smaller bar angle (20 deg instead of 30 deg) requires smaller distance uncertainty (20% instead of 30%) to explain the observations. Fourier components of the face-on density distribution of our models suggest that the MW does not have strong m=1 and/or m=3 spirals inside the solar radius.

The detection of radio emission from an exoplanet would constitute the best way to determine its magnetic field. Indeed, the presence of a planetary magnetic field is a necessary condition for radio emission via the Cyclotron Maser Instability. The presence of a magnetic field is, however, not sufficient. At the emission site, the local cyclotron frequency has to be sufficiently high compared to the local plasma frequency. As strong stellar insolation on a low-mass planet can lead to an extended planetary atmosphere, the magnetospheric plasma frequency depends on the planetary mass, its orbital distance, and its host star. We show that an extended planetary atmosphere can quench the radio emission. This seems to be true, in particular, for an important fraction of the planets less massive than approximately two Jupiter masses and with orbital distances below $\sim$0.2 AU. Most of the best candidates suggested by radio scaling laws lie in this parameter space. Taking this effect quenching into account will have important implications for the target selection of observation campaigns. At the same time, this effect will have consequences for the interpretation of observational data.

Cosmological parameters can be measured by comparing peculiar velocities with those predicted from a galaxy density field. Previous work has tested the accuracy of this approach with N-body simulations, but generally on idealised mock galaxy surveys. However, systematic biases may arise solely due to survey selection effects such as flux-limited samples, edge-effects, and complications due to the obscuration of the Galactic plane. In this work, we explore the impact of each of these effects independently and simultaneously, using the semi-analytic models from numerical simulations to generate mock catalogues that mimic the 2M++ density field. We find the reconstruction and analysis methods used for our 2M++ mocks produce a value of fsigma_8 that is biased high by a factor 1.04 \pm 0.01 compared to the true value. Moreover, a cosmic volume matching that of 2M++ has a cosmic variance uncertainty in fsigma_8 of ~5%. The systematic bias is a function of distance: it is unbiased close to the origin but is biased slightly high for distances in the range 100-180 Mpc/h. Correcting for this small bias, we find the linear fsigma_8 = 0.362 \pm 0.023. The predicted peculiar velocities from 2M++ have an error of 170 km/s that slowly increases with distance, exceeding 200 km/s only at distances of 180-200 Mpc/h. Finally, the residual bulk flow speeds found in previous work are shown to be not in conflict with those expected in the LCDM model.

``Extremely red quasars'' (ERQs) are a non-radio-selected, intrinsically luminous population of quasars at cosmic noon selected by their extremely red colour from rest-frame UV to mid-IR. ERQs are uniquely associated with exceptionally broad and blueshifted [\ion{O}{III}]\,$\lambda$5007 emission reaching speeds $>$6000\,km\,s$^{-1}$. We obtained laser-guided adaptive optics integral-field spectroscopic observations using Keck/OSIRIS and Gemini/NIFS of a sample of 10 ERQs with bolometric luminosities (10$^{47.0}$\textendash10$^{47.9}$)\,erg\,s$^{-1}$ at $z\sim$~(2.3\textendash3.0). The goal is to measure the sizes and spatially-resolved kinematics of the [\ion{O}{III}]-emitting regions. We study the surface brightness maps and aperture-extracted spectra and model the point-spread functions. We identify signs of merger activities in the continuum emissions. We identify physically distinct [\ion{O}{III}] kinematic components that are bimodal and respectively trace ERQ-driven outflows of velocity dispersion $\gtrsim$250\,km\,s$^{-1}$ and dynamically quiescent interstellar media. We find that the ERQ-driven ionized outflows are typically at $\sim$1\,kpc scales whereas the quiescent ionized gas extend to a few kpc. Compared to other luminous quasars the extremely fast ERQ-driven [\ion{O}{III}] outflows tend to be more compact, supporting the notion that ERQs are in a young stage of quasar/galaxy evolution and represent systems with unique physical conditions beyond orientation differences with other quasar populations. The kinematically quiescent [\ion{O}{III}] emissions in ERQs tend to be spatially-resolved but less extended than in other luminous quasars, which can be explained by global and patchy dust obscuration. The hint of ionization cones suggests some of the obscuration can be partially explained by a patchy torus.

Convincing evidence of a stochastic gravitational wave background has been found by the NANOGrav collaboration in the 15-Year data set. From this signal, we can evaluate the possibility of its source being from the early universe through the tensor perturbations induced by a massive spin-2 graviton field. We consider a time dependent model of the minimal theory of massive gravity, and find values of the graviton mass, mass cutoff time, and Hubble rate of inflation that amplify the energy spectra of primordial gravitational waves sufficiently to reproduce the signal from the NANOGrav data within 1-3 standard deviation. However, a suppression mechanism for high frequency modes must be introduced to conservatively obey the big bang nucleosynthesis (BBN) bound. While there are regions of the parameter space that reproduces the signal, it remains a challenge to simultaneously respect the BBN and cosmic microwave background (CMB) bounds without making the graviton mass cutoff time too deep into the matter dominated era.

Massive infrared dark clouds (IRDCs) are considered to host the earliest stages of high-mass star formation. In particular, 70 $\mu$m dark IRDCs are the colder and more quiescent clouds. At a scale of about 5000 au using formaldehyde (H2CO) emission, we investigate the kinetic temperature of dense cores in 12 IRDCs obtained from the pilot ALMA Survey of 70 $\mu$m dark High-mass clumps in Early Stages (ASHES). Compared to 1.3 mm dust continuum and other molecular lines, such as C18O and deuterated species, we find that H2CO is mainly sensitive to low-velocity outflow components rather than to quiescent gas expected in the early phases of star formation. The kinetic temperatures of these components range from 26 to 300 K. The Mach number reaches about 15 with an average value of about 4, suggesting that the velocity distribution of gas traced by H2CO is significantly influenced by a supersonic non-thermal component. In addition, we detect warm line emission from HC3N and OCS in 14 protostellar cores, which requires high excitation temperatures (Eu/k ~ 100 K). These results show that some of the embedded cores in the ASHES fields are in an advanced evolutionary stage, previously unexpected for 70 $\mu$m dark IRDCs.

Coronal mass ejections (CMEs) are large-scale eruptions with a typical radial size at 1 au of 0.21 au but their angular width in interplanetary space is still mostly unknown, especially for the magnetic ejecta (ME) part of the CME. We take advantage of STEREO-A angular separation of 20$^\circ$-60$^\circ$ from the Sun-Earth line from October 2020 to August 2022, and perform a two-part study to constrain the angular width of MEs in the ecliptic plane: a) we study all CMEs that are observed remotely to propagate between the Sun-STEREO-A and the Sun-Earth lines and determine how many impact one or both spacecraft in situ, and b) we investigate all in situ measurements at STEREO-A or at L1 of CMEs during the same time period to quantify how many are measured by the two spacecraft. A key finding is that, out of 21 CMEs propagating within 30$^\circ$ of either spacecraft, only four impacted both spacecraft and none provided clean magnetic cloud-like signatures at both spacecraft. Combining the two approaches, we conclude that the typical angular width of a ME at 1 au is $\sim$ 20$^\circ$-30$^\circ$, or 2-3 times less than often assumed and consistent with a 2:1 elliptical cross-section of an ellipsoidal ME. We discuss the consequences of this finding for future multi-spacecraft mission designs and for the coherence of CMEs.

The Breakthrough Listen search for intelligent life is, to date, the most extensive technosignature search of nearby celestial objects. We present a radio technosignature search of the centers of 97 nearby galaxies, observed by Breakthrough Listen at the Robert C. Byrd Green Bank Telescope. We performed a narrowband Doppler drift search using the turboSETI pipeline with a minimum signal-to-noise parameter threshold of 10, across a drift rate range of $\pm$ 4 Hz\ $s^{-1}$, with a spectral resolution of 3 Hz and a time resolution of $\sim$ 18.25 s. We removed radio frequency interference by using an on-source/off-source cadence pattern of six observations and discarding signals with Doppler drift rates of 0. We assess factors affecting the sensitivity of the Breakthrough Listen data reduction and search pipeline using signal injection and recovery techniques and apply new methods for the investigation of the RFI environment. We present results in four frequency bands covering 1 -- 11 GHz, and place constraints on the presence of transmitters with equivalent isotropic radiated power on the order of $10^{26}$ W, corresponding to the theoretical power consumption of Kardashev Type II civilizations.

This paper explores the dark energy phenomenon within the context of $f(R,L_m)$ gravity theory. Two specific non-linear $f (R, L_m)$ models are considered: $f(R,L_m)=\frac{R}{2}+L_m^\alpha$ and $f(R,L_m)=\frac{R}{2}+(1+\alpha R)L_m$, where the parameter $\alpha$ is free. Here, we adopt a parametrization form for the Hubble parameter in terms of redshift $z$ as $H(z)=H_0 \left[A(1+z)^3+B+\epsilon \log(1+z)\right]^\frac{1}{2}$, which allows for deviations from the standard $\Lambda$CDM model at both low and high redshifts. We then incorporate the Hubble parameter solution into the Friedmann equations for both models. We employ Bayesian analysis to estimate the constraints on the free parameters $H_0$, $A$, $B$, and $\epsilon$ using the Hubble measurements and the Pantheon dataset. Further, we investigate the evolution of key cosmological quantities, such as the deceleration parameter, energy density, pressure, EoS parameter, and energy conditions. The evolution of the deceleration parameter reveals a significant transition from a decelerating phase to an accelerating phase in the Universe. The EoS parameter exhibits quintessence-like behavior for both non-linear $f (R, L_m)$ models.

Microbursts are short duration intensifications in precipitating electron flux that are believed to be a significant contributor to electron losses in the magnetosphere. Microbursts have been observed in the form of bouncing electron packets, which offer a unique opportunity to study the properties of microbursts and their importance as a loss process. We present a collection of bouncing microbursts observed by the HILT instrument on SAMPEX from 1994-2004.We analyze the locations of the bouncing microbursts in L and MLT and find they align well with the properties of relativistic microbursts as a whole. We find that that the majority of bouncing microbursts observed by SAMPEX have scale sizes of 30km at the point of observation, or about 1500km when mapped to the equator.The time separation between the peaks of these bouncing microbursts is usually either half a bounce period or a whole bounce period.

Gravitational wave observations of binary black holes have revealed unexpected structure in the black hole mass distribution. Previous studies of the mass distribution employ physically-motivated phenomenological models and infer the parameters that directly control the features of the mass distribution that are allowed in their model, associating the constraints on those parameters with their physical motivations. In this work, we take an alternative approach in which we introduce a model parameterizing the underlying stellar and core-collapse physics and obtaining the remnant black hole distribution as a derived byproduct. In doing so, we directly constrain the stellar physics necessary to explain the astrophysical distribution of black hole properties under a given model. We apply this approach to modeling the mapping between stellar core mass and remnant black hole mass, including the effects of mass loss due to the pulsational pair instability supernova (PPISN) process, which has been proposed as an explanation for the observed excess of black holes at $\sim 35 M_\odot$. Placing constraints on the nuclear reaction rates necessary to explain the PPISN parameters, we conclude that the peak observed at $\sim 35 M_\odot$ is highly unlikely to be a signature from the PPISN process. This procedure can be applied to modeling any physical process that underlies the astrophysical mass distribution. Allowing the parameters of the core-remnant mass relationship to evolve with redshift permits correlated and physically reasonable changes in the location, shape, and amplitude of features in the mass function. We find that the current data are consistent with no redshift evolution in the core-remnant mass relationship, but ultimately place only weak constraints on the change of these parameters.

We present a $^{12}$CO($J$=2-1) survey of 60 local galaxies using data from the Atacama Large Millimeter/submillimeter Compact Array as part of the Extragalactic Database for Galaxy Evolution: the ACA EDGE survey. These galaxies all have integral field spectroscopy from the CALIFA survey. Compared to other local galaxy surveys, ACA EDGE is designed to mitigate selection effects based on CO brightness and morphological type. Of the 60 galaxies in ACA EDGE, 36 are on the star-formation main sequence, 13 are on the red sequence, and 11 lie in the ``green valley" transition between these sequences. We test how star formation quenching processes affect the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}=$SFR/$M_{\rm mol}$, and related quantities in galaxies with stellar masses $10\leq$log[$M_\star/$M$_\odot$]$\leq11.5$ covering the full range of morphological types. We observe a systematic decrease of the molecular-to-stellar mass fraction ($R^{\rm mol}_{\star}$) with decreasing level of star formation activity, with green valley galaxies having also lower SFE$_{\rm mol}$ than galaxies on the main sequence. On average, we find that the spatially resolved SFE$_{\rm mol}$ within the bulge region of green valley galaxies is lower than in the bulges of main sequence galaxies if we adopt a constant CO-to-H$_2$ conversion factor, $\alpha_{\rm CO}$. While efficiencies in main sequence galaxies remain almost constant with galactocentric radius, in green valley galaxies we note a systematic increase of SFE$_{\rm mol}$, $R^{\rm mol}_{\star}$, and specific star formation rate, sSFR, with increasing radius. Our results suggest that although gas depletion (or removal) seems to be the most important driver of the star-formation quenching in galaxies transiting through the green valley, a reduction in star formation efficiency is also required during this stage.

It has been suggested that the bulge-to-total stellar mass ratio or feedback from black holes (BHs), traced by the BH-to-(total stellar) mass ratio, might establish a galaxy's specific star formation rate (sSFR). We reveal that a galaxy's morphology -- reflecting its formation history, particularly accretions and mergers -- is a far better determinant of the sSFR. Consequently, we suggest that galaxy formation models which regulate the sSFR primarily through BH feedback prescriptions or bulge-regulated disc fragmentation consider acquisitions and mergers which establish the galaxy morphology. We additionally make several new observations regarding current ($z\sim0$) star-formation rates. (i) Galaxies with little to no star formation have bulges with an extensive range of stellar masses; bulge mass does not dictate presence/absence on the `star-forming main sequence'. (ii) The (wet merger)-built, dust-rich S0 galaxies are the `green valley' bridging population between elliptical galaxies on the `red sequence' and spiral galaxies on the blue star-forming main sequence. (iii) The dust-poor S0 galaxies are not on the star-forming main sequence nor in the `green valley'. Instead, they wait in the field for gas accretion and/or minor mergers to transform them into spiral galaxies. Mid-infrared sample selection can miss these (primordial) low dust-content and low stellar-luminosity S0 galaxies. Finally, the appearance of the quasi-triangular-shaped galaxy-assembly sequence, previously dubbed the Triangal, which tracks the morphological evolution of galaxies, is revealed in the sSFR-(stellar mass) diagram.

We study the dynamics of the quintom dark energy model using state-of-the-art cosmological observations. The set of equations has been converted into an autonomous system using suitable transformations of the variables. We have discussed the fixed points of the model and the general phase-space behavior, in particular, in finding the existence of the tracker solutions for this model. The observations suggest that at late times the phantom field should dominate the dark energy sector with an approximately 15% share to the quintessence counterpart, and with both fields tracking the background at early times. A Bayesian model comparison with LambdaCDM has also been done by computing the Bayes factor and a positive preference has been obtained for the quintom model. Although not fully resolved, the Hubble tension can be reduced to 2.6{\sigma} when compared with the value of H0 reported in [1] and to 1.6{\sigma} when compared with that of [2].

The Hubble tension persists as a challenge in cosmology. Even early dark energy (EDE) models, initially considered the most promising for alleviating the Hubble tension, fall short of addressing the issue without exacerbating other tensions, such as the $S_8$ tension. Considering that a negative dark matter (DM) equation of state (EoS) parameter is conducive to reduce the value of $\sigma_8$ parameter, in this paper, we extend the axion-like EDE model in this paper by replacing the cold dark matter (CDM) with DM characterized by a constant EoS $w_{\rm dm}$ (referred as WDM hereafter). We then impose constraints on this axion-like EDE extension model, along with three other models: the axion-like EDE model, $\Lambda$WDM, and $\Lambda$CDM. These constraints are derived from a comprehensive analysis incorporating data from the Planck 2018 cosmic microwave background (CMB), baryon acoustic oscillations (BAO), the Pantheon compilation, as well as a prior on $H_0$ (i.e., $H_0=73.04\pm1.04$, based on the latest local measurement by Riess et al.) and a Gaussianized prior on $S_8$ (i.e., $S_8=0.766\pm0.017$, determined through the joint analysis of KID1000+BOSS+2dLenS). We find that although the new model maintains the ability to alleviate the Hubble tension to $\sim$ 1.4$\sigma$, it still exacerbate the $S_8$ tension to a level similar to that of the axion-like EDE model.

The observational data of high redshift galaxies become increasingly abundant, especially since the operation of the James Webb Space Telescope (JWST), which allows us to verify and optimize the galaxy formation model at high redshifts. In this work, we investigate the merging history of massive galaxies at $3 < z < 6$ using a well-developed semi-analytic galaxy formation catalogue. We find that the major merger rate increases with redshift up to 3 and then flattens. The fraction of wet mergers, during which the sum of the cold gas mass is higher than the sum of the stellar mass in two merging galaxies, also increases from $\sim$ 34\% at $z = 0$ to 96\% at $z = 3$. Interestingly, almost all major mergers are wet at $z > 3$ . This can be attributed to the high fraction ($> 50\%$) of cold gas at $z > 3$. In addition, we study some special systems of massive merging galaxies at $3 < z < 6$, including the massive gas-rich major merging systems and extreme dense proto-clusters, and investigate the supermassive black hole-dark matter halo mass relation and dual AGNs. We find that the galaxy formation model reproduces the incidence of those observed massive galaxies, but fails to reproduce the relation between the supermassive black hole mass and the dark matter halo mass at $z \sim 6$. The latter requires more careful estimates of the supermassive black hole masses observationally. Otherwise, it could suggest modifications of the modeling of the supermassive black hole growth at high redshifts.

This book was created as part of the SIRIUS B VERGE program to orient students to astrophysics as a broad field. The 2023-2024 VERGE program and the printing of this book is funded by the Women and Girls in Astronomy Program via the International Astronomical Union's North American Regional Office of Astronomy for Development and the Heising-Simons Foundation; as a result, this document is written by women in astronomy for girls who are looking to pursue the field. However, given its universal nature, the material covered in this guide is useful for anyone interested in pursuing astrophysics professionally.

We present a systematic study of the putative hot gas corona around late-type galaxies (LTGs) residing in the Virgo cluster, based on archival Chandra observations. Our sample consists of 21 nearly edge-on galaxies representing a star formation rate (SFR) range of ($0.2-3\rm~M_\odot~yr^{-1}$) a stellar mass ($M_*$) range of $(0.2-10) \times 10^{10}\rm~M_{\odot}$, the majority of which have not been explored with high-sensitivity X-ray observations so far. Significant extraplanar diffuse X-ray (0.5-2 keV) emission is detected in only three LTGs, which are also the three galaxies with the highest SFR. A stacking analysis is performed for the remaining galaxies without individual detection, dividing the whole sample into two subsets based on SFR, stellar mass, or specific SFR. Only the high-SFR bin yields a significant detection, which has a value of $L\rm_X \sim3\times10^{38}\rm~erg~s^{-1}$ per galaxy. The stacked extraplanar X-ray signals of the Virgo LTGs are consistent with the empirical $L\rm_X - SFR$ and $L\rm_X - M_*$ relations found among highly inclined disk galaxies in the field, but appear to be systematically lower than that of a comparison sample of simulated cluster star-formation galaxies identified from the Illustris-TNG100 simulation. The apparent paucity of hot gas coronae in the sampled Virgo LTGs might be understood as the net outcome of the long-lasting effect of ram pressure stripping exerted by the hot intra-cluster medium and in-disk star-forming activity acting on shorter timescales. A better understanding of the roles of environmental effects in regulating the hot gas content of cluster galaxies invites sensitive X-ray observations for a large galaxy sample.

We present a study of blue straggler stars (BSSs) of open cluster NGC 7142 using AstroSat/UVIT data and other archival data. Using a machine learning-based algorithm, ML-MOC, on Gaia DR3 data, we find 546 sources as cluster members. Based on the location on the Gaia color-magnitude diagram, we identify ten BSS candidates, also detected in UVIT/F148W filter. We study the variable nature of BSSs by constructing their light curves using the TESS data. Two BSSs reported as eclipsing binaries in Gaia DR3 are confirmed to be eclipsing binaries based on our analysis and also show the presence of hot companions as per the multi-wavelength spectral energy distributions (SEDs). The physical parameters of the hot companions of these two BSSs derived by fitting binary models to their light curves and those derived from the SEDs are found to be in good agreement. Additionally, five more BSSs in the cluster shows UV excess, four of which are likely to have a hot companion based on SEDs. The hot companions with the estimated temperatures $\sim$14000 $-$ 28000 K, radii $\sim$0.01 $-$ 0.05 R$_{\odot}$, and luminosities $\sim$0.03 $-$ 0.1 L$_{\odot}$, are inferred to be extremely low mass ($<$ 0.2 M$_{\odot}$), low-mass (0.2 $-$ 0.4 M$_{\odot}$), normal-mass (0.4 $-$ 0.6 M$_{\odot}$), and high-mass ($>$ 0.6 M$_{\odot}$) white dwarfs (WD). For the first time in an open cluster, we find the entire range of masses in WDs found as hot companions of BSSs. These masses imply that the Case-A/Case-B mass transfer as well as merger are responsible for the formation of at least 60$\%$ of the BSSs of this cluster.

Recent gravitational wave (GW) observations of binary black hole (BH) mergers and the stochastic GW background have triggered renewed interest in primordial black holes (PBHs) in the stellar-mass ($\sim 10 - 100\ \rm M_\odot$) and supermassive regimes ($\sim 10^7 - 10^{11}\ \rm M_\odot$). Although only a small fraction ($\lesssim 1\%$) of dark matter (DM) in the form of PBHs is required to explain such observations, these PBHs may play important roles in early structure/star/galaxy formation. In this chapter, we combine semi-analytical analysis and cosmological simulations to explore the possible impact of PBHs on the formation of the first stars and galaxies, taking into account two (competing) effects of PBHs: acceleration of structure formation and gas heating by BH accretion feedback. We find that the impact of stellar-mass PBHs (allowed by existing observational constraints) on primordial star formation is likely minor, although they do alter the properties of the first star-forming halos/clouds and can potentially trigger the formation of massive BHs, while supermassive PBHs serve as seeds of massive structures that can explain the apparent overabundance of massive galaxies in recent observations. Our tentative models and results call for future studies with improved modeling of the interactions between PBHs, particle DM, and baryons to better understand the impact of PBHs on early star/galaxy/structure formation and their imprints in high-redshift observations.

We report a faint non-axisymmetric structure in NGC\,1087 through the use of JWST Near Infrared Camera { (NIRCam)}, with an associated kinematic counterpart observed as an oval distortion in the stellar velocity map, \ha~and CO~$J=2\rightarrow1$ velocity fields. This structure is not evident in the MUSE optical continuum images but only revealed in the near-IR with the F200W and F300M band filters at $2\mu$m and $3\mu$m respectively. Due to its elongation, this structure resembles a stellar bar although with remarkable differences with respect to conventional stellar bars. Most of the near-IR emission is concentrated within $6\arcsec~\sim500$~pc with a maximum extension up to 1.2~kpc. The spatial extension of the large-scale non-circular motions is coincident with the bar, which undoubtedly confirms the presence of a non-axisymmetric perturbation in the potential of NGC\,1087. The oval distortion is enhanced in CO due to its dynamically cold nature rather than in \ha. We found that the kinematics in all phases including stellar, ionized and molecular, can be described simultaneously by a model containing a bisymmetric perturbation; however, we find that an inflow model of gas along the bar major axis is also likely. Furthermore the molecular mass inflow rate associated can explain the observed star formation rate in the bar. This reinforces the idea that bars are mechanisms for transporting gas and triggering star formation. This work contributes to our understanding of non-axisymmetry in galaxies using the most sophisticated data so far.

It is well known that small-scale magnetism dominates the surface magnetic field topologies of active late-type stars. However, little information is available on the spatial distribution of this key magnetic field component. Here, we take advantage of the recently developed magnetic field diagnostic procedure relying on the magnetic intensification of iron atomic lines in the optical. We extend this methodology from measuring a single average field strength value to simultaneous Doppler imaging reconstruction of the two-dimensional maps of temperature and magnetic field strength. We applied this novel surface mapping approach to two spectroscopic data sets of the young active Sun-like star LQ Hya. For both epochs, we found a fairly uniform field strength distribution, apart from a latitudinal trend of the field strength increasing from 1.5-2.0 kG at low latitudes to 3.0-3.5 kG, close to the rotational poles. This distribution of the small-scale field does not display a clear correlation with the locations of temperature spots or the global magnetic field structure reconstructed for the same epochs.

After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a mass-loading factor for every galaxy, which is a first for a SAM. Building on the model's resolved angular-momentum structure of galaxies, we now consider the heating of stellar discs, delivering predictions for disc structure both radially and vertically, another SAM first. We add a further dimension to stellar discs by tracking the distribution of stellar ages in each annulus. Each annulus--age bin has its own velocity dispersion and metallicity evolved in the model. This allows Dark Sage to make structural predictions for galaxies that previously only hydrodynamic simulations could. We present the model as run on the merger trees of the highest-resolution gravity-only simulation of the MillenniumTNG suite. Despite its additional complexity relative to other SAMs, Dark Sage only has three free parameters, the least of any SAM, which we calibrate exclusively against the cosmic star formation history and the $z=0$ stellar and HI mass functions using a particle-swarm optimisation method. The Dark Sage codebase, written in C and Python, is publicly available at https://github.com/arhstevens/DarkSage

General-relativistic magneto-hydrodynamical (GRMHD) simulations of accreting black holes suggest that the accretion flows form toroidal structures embedded in a large scale component of magnetic field, which becomes organized on length-scales exceeding the gravitational radius of the central black hole. Magnetic field grows gradually until a Magnetically Arrested Disc (MAD) develops that diminishes or inhibits further accretion. We study an outflow that develops in the MAD state in 3D GRMHD simulations. We show that the outflow can be accelerated to relativistic velocities and persist over the course of our simulation. We compare the properties of the outflow from MAD discs with those launched by orbiting secondary at close orbit. The main difference is that the orbiting body launches a more coherent, quasiperiodic ultrafast outflow at lower velocities ($v<0.5c$) while the outflow launched in the MAD state (without the body) has a stochastic behaviour and has an approximately flat velocity distribution between lower and higher outflow velocities, $0.2c<v<0.3c$ and $v>0.5c$.

The detections of Lyman-$\alpha$ ($\rm Ly\alpha$) emission in galaxies with redshifts above 5 are of utmost importance for constraining the cosmic reionization timeline, yet such detections are usually based on slit spectroscopy. Here we investigate the significant bias induced by slit placement on the estimate of $\rm Ly\alpha$ escape fraction ( $f_{\rm esc}^{\mathrm{Ly\alpha}}$), by presenting a galaxy (dubbed A2744-z6Lya) at $z=5.66$ where its deep JWST NIRSpec prism spectroscopy completely misses the strong $\rm Ly\alpha$ emission detected in the MUSE data. A2744-z6Lya exhibits a pronounced UV continuum with an extremely steep spectral slope of $\beta=-2.574_{-0.008}^{+0.008}$, and it has a stellar mass of $\mathrm{\sim10^{8.82}~M_\odot}$, a star-formation rate of $\mathrm{\sim8.35~M_\odot yr^{-1}}$ and gas-phase metallicity of $\mathrm{12+log\,(O/H)\sim7.88}$. The observed flux and rest-frame equivalent width of its Ly$\alpha$ from MUSE spectroscopy are $1.2\times \rm 10^{-16} erg~s^{-1}cm^{-2}$ and 75\r{A}, equivalent to $f_{\rm esc}^{\mathrm{Ly\alpha}}=78\pm4 \%$. However, its Ly$\alpha$ non-detection from JWST NIRSpec gives a 5-$\sigma$ upper limit of $<13 \%$, in stark contrast to that derived from MUSE. To explore the reasons for this bias, we perform spatially resolved stellar population analysis of A2744-z6Lya using the JWST NIRCam imaging data to construct 2-dimensional maps of SFR, dust extinction and neutral hydrogen column density. We find that the absence of Ly$\alpha$ in the slit regions probably stems from both the resonance scattering effect of neutral hydrogen and dust extinction. Through analyzing an extreme case in detail, this work highlights the important caveat of inferring $f_{\rm esc}^{\mathrm{Ly\alpha}}$ from slit spectroscopy, particularly when using the JWST multiplexed NIRSpec microshutter assembly.

Magnetic active regions on the Sun are harbingers of space weather events. Understanding the physics of how they form and evolve will improve space weather forecasting. Our aim is to characterise the surface magnetic field and flows for a sample of active regions with persistent magnetic bipoles prior to emergence. We identified 42 emerging active regions (EARs), in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey (Schunker et al. 2016), associated with small magnetic bipoles at least one day before the time of emergence. We then identified a contrasting sample of 42 EARs that emerge more abruptly without bipoles before emergence. We computed the supergranulation scale surface flows using helioseismic holography. We averaged the flow maps and magnetic field maps over all active regions in each sample at each time interval from 2 days before emergence to 1 day after. We found that EARs associated with a persistent pre-emergence bipole evolve to be, on average, lower flux active regions than EARs that emerge more abruptly. Further, we found that the EARs that emerge more abruptly do so with a diverging flow of $(3\pm 0.6) \times 10^{-6}$ s$^{-1}$ on the order of 50-100 ms$^{-1}$. Our results suggest that there is a statistical dependence of the surface flow signature throughout the emergence process on the maximum magnetic flux of the active region.

We report the discovery of a radio relic in the northeastern periphery of the cluster Abell 2108 (A2108). A2108 is a part of the uGMRT LOw-MAss Galaxy Cluster Survey (GLOMACS), where our main aim is to search for diffuse radio emission signatures in very sparsely explored low-mass galaxy clusters using uGMRT band-3 (central frequency 400 MHz). We used our uGMRT band-3 data along with the existing archival band-3 uGMRT data to improve image sensitivity. Along with the previously reported southwestern relic, the discovery of the new relic makes A2108 one of the few low-mass clusters hosting double relics. The new relic spans over a region of 610 kpc $\times$ 310 kpc and, interestingly, differs considerably in size and morphology from the other relic. With XMM-Newton science archive data, we also report the tentative detection of a mildly supersonic shock of Mach number $\mathcal{M}_\mathrm{SB}=1.42$ and $\mathcal{M}_\mathrm{T} = 1.43$ from the surface brightness and temperature discontinuities, respectively near this newly found relic. Both the relics in A2108 are found to be significantly under-luminous compared to other double relic systems in the mass-luminosity plane. Though mild supersonic shocks resulting from an off-axis merger could have influenced their origin, we hypothesize that further local environments have played a crucial role in shaping their morphology.

The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and memory-saving Runge-Kutta (RK) scheme. It handles non-ideal effects through super-time stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a non-reflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up till several thousands of CPUs.

We report on recent results from our successful and pioneering observational program with ALMA to study emerging ultracom pact HII regions of pre-Planetary Nebulae (pPNe) using mm-wavelength recombination lines (mRRLs) as new optimal tracers. We focus on our study of two poster-child pPNe, namely, M2-9 and CRL618. We reveal the structure and kinematics of the en igmatic inner nebular regions of these objects with an unprecedented angular resolution down to 20-30mas (~15-30AU). For both targets, the ionized central regions are elongated along the main symmetry axis of the large-scale nebulae, consiste nt with bipolar winds, and show notable axial velocity gradients with expansion velocities of up to ~100km/s. The H30a pr ofiles exhibit time variability, reflecting changes in the physical properties and kinematics on scales of a few years. O ur ongoing analysis employs 3D, non-LTE radiative transfer modeling, providing a detailed description of the innermost la yers of these well known pPNe with exceptional clarity.

We present the results of a systematic search of the Transiting Exoplanet Survey Satellite (TESS) 2-min cadence data for new rapidly oscillating Ap (roAp) stars observed during the Cycle 2 phase of its mission. We find seven new roAp stars previously unreported as such and present the analysis of a further 25 roAp stars that are already known. Three of the new stars show multiperiodic pulsations, while all new members are rotationally variable stars, leading to almost 70 per cent (22) of the roAp stars presented being $\alpha^2$ CVn-type variable stars. We show that targeted observations of known chemically peculiar stars are likely to overlook many new roAp stars, and demonstrate that multi-epoch observations are necessary to see pulsational behaviour changes. We find a lack of roAp stars close to the blue edge of the theoretical roAp instability strip, and reaffirm that mode instability is observed more frequently with precise, space-based observations. In addition to the Cycle 2 observations, we analyse TESS data for all known roAp stars. This amounts to 18 further roAp stars observed by TESS. Finally, we list six known roAp stars that TESS is yet to observe. We deduce that the incidence of roAp stars amongst the Ap star population is just 5.5 per cent, raising fundamental questions about the conditions required to excite pulsations in Ap stars. This work, coupled with our previous work on roAp stars in Cycle 1 observations, presents the most comprehensive, homogeneous study of the roAp stars in the TESS nominal mission, with a collection of 112 confirmed roAp stars in total.

We analyze photometric observations of the symbiotic star MWC 560 in B and V bands obtained during the period 1990-2023. We estimate the luminosity and the mass accretion rate of the hot component. We find that the luminosity varies in the range from 200 $L_\odot$ to 3000 $L_\odot$, corresponding to a mass accretion rate in the range $1 \times 10^{-7} - 2 \times 10^{-6} \; M_\odot \; yr^{-1}$ (for a 0.9 $M_\odot$ white dwarf and distance 2217 pc). The optical flickering disappears at mass accretion rate of about $1\times 10^{-6}\; M_\odot \; yr^{-1}$, which sets an upper limit for the short-term variability from accreting white dwarf.

The mergers of two neutron stars are exceptional multi-messenger events that enable us to probe fundamental physics in one of the most extreme environments in the Universe. Multi-wavelength follow-up observations are essential in order to probe the physics of the outflows from and remnants of these neutron star mergers, both when detected as short gamma-ray bursts (GRBs) and as gravitational wave events. Rapid follow-up can provide localisations for targeted deep follow-up observations and, ideally, a distance measurement, which constrains for instance the energetics of the merger. A key outstanding question is the remnant's nature: with its expected mass and rapid spin, it could either be a black hole or a supramassive, likely highly magnetised neutron star (a magnetar). Both can power a GRB, but rapidly spinning magnetars are additionally predicted to emit coherent radio bursts following their formation and may constitute a small fraction of the progenitors of fast radio bursts. Black holes, by contrast, are not expected to emit coherent radio bursts in the time following the GRB itself. Here we present rapid follow-up observations of the short GRB 201006A using LOFAR. We have detected a 5.6$\sigma$, short, coherent radio flash at 144 MHz at 76.6 mins post-burst. This radio flash is 27 arcsec offset from the GRB location, which has a probability of occurring by chance of $\sim$0.5% (2.6$\sigma$) when accounting for measurement uncertainties. Despite the offset, we show that the probability of finding an unrelated transient within 40 arcsec of the GRB location is $<10^{-6}$ and conclude that this is likely to be the radio counterpart to GRB 201006A. The radio flash is tentatively (2.5$\sigma$) shown to be highly dispersed, allowing a distance estimate, corresponding to a redshift of $0.58\pm0.06$, that is in the range of typical short GRB distances. Using the estimated distance, the...

This work is part of the BinaMIcS project, the aim of which is to understand the interaction between binarity and magnetism in close binary systems. All the studied spectroscopic binaries targeted by the BinaMIcS project encompass hot massive and intermediate-mass stars on the main sequence, as well as cool stars over a wide range of evolutionary stages. The present paper focuses on the binary system FK Aqr, which is composed of two early M dwarfs. Both stars are already known to be magnetically active based on their light curves and detected flare activity. In addition, the two components have large convective envelopes with masses just above the fully convective limit, making the system an ideal target for studying effect of binarity on stellar dynamos. We use spectropolarimetric observations obtained with ESPaDOnS at CFHT in September 2014. Mean Stokes I and V line profiles are extracted using the least-squares deconvolution (LSD) method. The radial velocities of the two components are measured from the LSD Stokes I profiles and are combined with interferometric measurements in order to constrain the orbital parameters of the system. The longitudinal magnetic fields Bl and chromospheric activity indicators are measured from the LSD mean line profiles. The rotational modulation of the Stokes V profiles is used to reconstruct the surface magnetic field structures of both stars via the Zeeman Doppler imaging (ZDI) inversion technique. Maps of the surface magnetic field structures of both components of FK Aqr are presented for the first time. Our study shows that both components host similar large-scale magnetic fields of moderate intensity (Bmean ~ 0.25 kG); both are predominantly poloidal and feature a strong axisymmetric dipolar component. (abridged)

Cosmological studies have now entered Stage IV according to the Dark Energy Task Force prescription, thanks to new missions (Euclid, Rubin Observatory, SRG/eROSITA) that are expected to provide the required ultimate accuracy in the dark energy (DE) equation of state (EoS). However, none of these projects have the power to systematically unveil the galaxy cluster population at $z>1$. There therefore remains the need for an ATHENA-like mission to run independent cosmological investigations and scrutinise the consistency between the results from the $0<z<1$ and $1<z<2$ epochs. We study the constraints on the DE EoS and on primordial non-Gaussanities for typical X-ray cluster surveys executed by ATHENA. We consider two survey designs: 50 deg$^2$ at 80ks (survey A) and 200 deg$^2$ at 20ks (survey B). We analytically derive cluster counts in a space of observable properties, and predict the cosmological potential of the corresponding samples with a Fisher analysis. The achieved depth allows us to unveil the halo mass function down to the group scale out to $z=2$. We predict the detection of thousands of clusters down to a few 10$^{13} h^{-1} M_{\odot}$, in particular 940 and 1400 clusters for surveys A and B, respectively, at $z>1$. Such samples will allow a detailed modelling of the evolution of cluster physics along with a standalone cosmological analysis. Our results suggest that survey B has the optimal design as it provides greater statistics. Remarkably, high-$z$ clusters, despite representing 15% or less of the full samples, allow a significant reduction of the uncertainty on the cosmological parameters: $\Delta w_a$ is reduced by a factor of 2.3 and $\Delta f_{NL}^{loc}$ by a factor of 3. Inventorying the high-$z$ X-ray cluster population can play a crucial role in ensuring overall cosmological consistency. This will be the major aim of future new-generation ATHENA-like missions.

Gravitational waves, detected a century after they were first theorized, are spacetime distortions caused by some of the most cataclysmic events in the universe, including black hole mergers and supernovae. The successful detection of these waves has been made possible by ingenious detectors designed by human experts. Beyond these successful designs, the vast space of experimental configurations remains largely unexplored, offering an exciting territory potentially rich in innovative and unconventional detection strategies. Here, we demonstrate the application of artificial intelligence (AI) to systematically explore this enormous space, revealing novel topologies for gravitational wave (GW) detectors that outperform current next-generation designs under realistic experimental constraints. Our results span a broad range of astrophysical targets, such as black hole and neutron star mergers, supernovae, and primordial GW sources. Moreover, we are able to conceptualize the initially unorthodox discovered designs, emphasizing the potential of using AI algorithms not only in discovering but also in understanding these novel topologies. We've assembled more than 50 superior solutions in a publicly available Gravitational Wave Detector Zoo which could lead to many new surprising techniques. At a bigger picture, our approach is not limited to gravitational wave detectors and can be extended to AI-driven design of experiments across diverse domains of fundamental physics.

During the Epoch of Reionization (EoR), the ultraviolet radiation from the first stars and galaxies ionised the neutral hydrogen of the intergalactic medium, which can emit radiation through its 21 cm hyperfine transition. This 21 cm signal is a direct probe of the first light sources in the early Universe. Measuring the 21 cm power spectrum is a key science goal for the future Square Kilometre Array (SKA), however, observing and interpreting it is a challenging task. Another high-potential probe of the EoR is the patchy kinetic Sunyaev-Zel'dovich effect (pkSZ), observed as a foreground to the primary cosmic microwave background temperature anisotropies on small scales. Despite recent promising measurements by ground-based telescopes, placing constraints on reionization from pkSZ observations is a non-trivial task, subject to strong model dependence. In this work, we propose to alleviate the difficulties in observing and interpreting the 21 cm and pkSZ power spectra by combining them. With a simple yet effective parametric model that establishes a formal connection between them, we are able to jointly fit mock 21 cm and pkSZ data points. We confirm that these two observables provide complementary information on reionization, leading to significantly improved constraints when combined. We demonstrate that with as few as two measurements of the 21 cm power spectrum with 100 hours of observations with the SKA, as well as a single $\ell=3000$ pkSZ data point, we can reconstruct the reionization history of the Universe and its morphology. We find that the reionization global history (morphology) is better constrained with two 21 cm measurements at different redshifts (scales). Therefore, a combined analysis of the two probes will give access to tighter constraints on cosmic reionization even in the early stages of 21 cm detections.

Calibration of pixel values is a fundamental step for accurate measurements in astronomical imaging. In astronomical jargon this is known as estimating zero point magnitude. Here, we introduce a newly added script in GNU Astronomy Utilities (Gnuastro) version 0.20 for the zero point magnitude estimation, named: astscript-zeropoint. The script offers numerous features, such as the flexibility to use either image(s) or a catalog as the reference dataset. Additionally, steps are parallelized to enhance efficiency for big data. Thanks to Gnuastro's minimal dependencies, the script is both flexible and portable. The figures of this research note are reproducible with Maneage, on the Git commit c89275e.

T Pyx is one of the most enigmatic recurrent novae, and it has been proposed as a potential Galactic type-Ia supernova progenitor. Using spatially-resolved data obtained with MUSE, we characterized the geometrical distribution of the material expelled in previous outbursts surrounding the white dwarf progenitor. We used a 3D model for the ejecta to determine the geometric distribution of the extended remnant. We have also calculated the nebular parallax distance ($d = 3.55 \pm 0.77$ kpc) based on the measured velocity and spatial shift of the 2011 bipolar ejecta. These findings confirm previous results, including data from the GAIA mission. The remnant of T Pyx can be described by a two-component model, consisting of a tilted ring at $i = 63.7$ deg, relative to its normal vector and by fast bipolar ejecta perpendicular to the plane of the equatorial ring. We find an upper limit for the bipolar outflow ejected mass in 2011 of the bipolar outflow of $M_{ej,b} < (3.0 \pm 1.0) \times 10^{-6}$ M$_{\odot}$, which is lower than previous estimates given in the literature. However, only a detailed physical study of the equatorial component could provide an accurate estimate of the total ejecta of the last outburst, a fundamental step to understand if T Pyx will end its life as a type-Ia supernova.

Neutrinos are believed to be the most abundant fermions in the Universe, but their masses are unknown, except for being non-zero but much smaller than other fermions. Cosmological relic neutrinos could also have non-zero chemical potentials (or asymmetries). We develop a fast neutrino-involved N-body simulation code to include both neutrino mass and asymmetry effects. Using it, we investigate the neutrino effects on the matter pairwise velocity, which itself is an interesting probe of cosmology. We find that for light-halo ($[10^{11},10^{13}]\ M_\odot$) mean pairwise velocity, in the transition range ($[4,15]\ \mathrm{Mpc}$), the effects of neutrino masses overwhelm the effects of neutrino asymmetries, while in the two-halo-group range ($[25,50]\ \mathrm{Mpc}$), for both light and heavy haloes ($[10^{13},10^{15}]\ M_\odot$), the effects of neutrino asymmetries dominate, making it possible to disentangle the two effects. We provide fitting formulae to quantify the effects of neutrino mass and asymmetry on halo-halo pairwise velocities.

We present a study on the star formation histories (SFHs) of galaxies covering the range $10^{4}$ < M$_{\star}$/M$_{\odot}$ < $10^{12}$, leveraging full spectral fitting algorithms. Our sample consists of 31 dwarf galaxies from the SAMI-Fornax Survey with stellar masses between $10^{7}$-$10^{9.5} M_{\odot}$, early-type galaxies from the ATLAS$^{3D}$ project with stellar masses between $10^{10}$-$10^{12} M_{\odot}$, and dwarf galaxies that are satellites of Andromeda and the Milky Way, with $10^{4}$ < M$_{\star}$/M$_{\odot}$ < $10^{8}$. We find that galaxies from $10^{7}$-$10^{8} M_{\odot}$ exhibit the smallest star formation rates (SFRs), while the SFR increase as we move down or up in mass. In this sense, we find that some $10^{5} M_{\odot}$ galaxies have cumulative SFHs that are comparable to those of $10^{12} M_{\odot}$ galaxies. Our study shows that the evolution of giant galaxies is primarily governed by their internal properties, with timescales that do not depend on their environmental location. In contrast, dwarf galaxies below $10^{8} M_{\odot}$ can be significantly affected in dense environments, such as the inner regions of a cluster, that severely quench the galaxies before the assembly of their 50% present-day mass. We find that, only dwarfs with stellar masses between $10^{7}$-$10^{9} M_{\odot}$ actively form stars nowadays, while less massive galaxies seem to remain unaffected by the environment due to the expulsion of most of their gas at an early stage in their evolution. Our study highlights and corroborates a critical threshold around $10^{8}-10^{9} M_{\odot}$ in galaxy evolution from previous studies, separating more massive galaxies minimally impacted by the environment from those less massive galaxies quenched by it.

An explanation is proposed for the appearance of slowly drift shadow bursts in the dynamic spectrum of Jupiter against the background of decameter radio emission with a quasi-harmonic structure. Background radio emission is caused by hot ions with a loss cone type distribution function, which generate ion cyclotron waves due to the effect of double plasma resonance. A flow of hot ions with a distribution function of the Maxwell type is injected into the source region, fills the loss cone of generating ions and interrupts the generation of ion cyclotron waves due to the filling of the loss cone. The condition under which instability breaks down is obtained, and the optimal values of the parameters of the injected ions necessary for the occurrence of bursts in absorption are determined.

Extremely large telescopes (ELTs) are considered worldwide to be one of the highest priorities in ground-based astronomy. The European Southern Observatory (ESO) is developing an ELT that will have a 39 m main mirror and will be the largest visible and infrared light telescope in the world. The ELT will be equipped with a lineup of cutting-edge instruments, designed to cover a wide range of scientific possibilities. The leap forwards with the ELT can lead to a paradigm shift in our perception of the Universe, much as Galileo's telescope did 400 years ago. We illustrate here the various components of the ELT, including the dome and main structure, the five mirrors, and the telescope systems. We then describe the ELT instrumentation and some of the astronomical topics it will address. We then conclude by examining the synergies with other astronomical facilities.

PSR J1048$-$5832 (B1046$-$58) is a Vela-like pulsar that has exhibited multiple glitch events. In this study, we analyze the timing data spanning nearly 16 years, acquired from both the Fermi Gamma-ray Space Telescope and the Parkes 64 m radio telescope. As a result, a total of five glitches are detected within this dataset. Among them, a previously unknown small glitch is newly found at MJD 56985(9) (November 24, 2014), which is also the smallest glitch recorded from this source so far. The increments of the spin frequency and its first derivative are $\Delta \nu \approx 2.2(3) \times 10^ {-8} $ Hz, and $\Delta \dot{\nu} \approx 3(2) \times 10^ {-15} $ s$^{-2}$, respectively. Significant changes in the integrated mean pulse profile are observed following two of the five glitch events, notably in the radio band. We discuss the influence of glitches on the pulsar's emission properties due to superfluid dynamics accompanied by crustquake, as well as the constraints on the equation of state.

Context: T CrB (NOVA CrB 1946) is a famous recurrent nova with a recurrence timescale of 80 years. Aims: We aim to estimate the colours, luminosity, and mass-accretion rate for T CrB (NOVA CrB 1946) during and after the superactive state. Methods and Results: We performed and analysed $UBV$ photometry of the recurrent nova T~CrB. For the hot component of T~CrB, we find average dereddened colours of $(U-B)_0 = -0.70 \pm 0.08$ and $(B-V)_0 = 0.23 \pm 0.06$, which correspond to an effective temperature of $9400 \pm 500$~K and an optical luminosity of $40-110~L_\odot$ during the superactive state (2016-2022). After the end of the superactive state, the hot component became significantly redder, $(U-B)_0 \approx -0.3$ and $(B-V)_0 \approx 0.6$ in August 2023, and its luminosity decreased markedly to $20-25$~$L_\odot$ in April-May 2023, and to $8-9~L_\odot$ in August 2023. The total mass accreted during the superactive state from 2014 to 2023 is $\sim 2 \times 10^{-7}$~M$_\odot$. This is a significant fraction of the mass required to cause a thermonuclear runaway (TNR). Overall our results support a model in which a large accretion disc acts as a reservoir with increased accretion rate onto the central white dwarf during disc high states, ultimately leading to a TNR explosion, which now seems to be imminent.

We present the discovery of the most distant OH megamaser to be observed in the main lines, using data from the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. At a newly measured redshift of $z = 0.7092$, the system has strong emission in both the 1665MHz ($L \approx 2500$ L$_{\odot}$) and 1667 MHz ($L \approx 4.5\times10^4$ L$_{\odot}$) transitions, with both narrow and broad components. We interpret the broad line as a high-velocity-dispersion component of the 1667 MHz transition, with velocity $v \sim 330$km s$^{-1}$ with respect to the systemic velocity. The host galaxy has a stellar mass of $M_{\star} = 2.95 \times 10^{10}$ M$_{\odot}$ and a star-formation rate of SFR = 371 M$_{\odot}$yr$^{-1}$, placing it $\sim 1.5$dex above the main sequence for star-forming galaxies at this redshift, and can be classified as an ultra-luminous infrared galaxy. Alongside the optical imaging data, which exhibits evidence for a tidal tail, this suggests that the OH megamaser arises from a system that is currently undergoing a merger, which is stimulating star formation and providing the necessary conditions for pumping the OH molecule to saturation. The OHM is likely to be lensed, with a magnification factor of $\sim 2.5$, and perhaps more if the maser emitting region is compact and suitably offset relative to the centroid of its host galaxy's optical light. This discovery demonstrates that spectral line mapping with the new generation of radio interferometers may provide important information on the cosmic merger history of galaxies.

We present an observational approach for the independent detection of the prompt optical emission of long gamma-ray bursts (GRBs). For this purpose, we explore the potential of the Large Array Survey Telescope (LAST). This array of small optical telescopes can be used to scan a wide region of the sky, and to focus on a smaller field of view with increased sensitivity, as needed. The modularity of the array facilitates dynamic scanning of multiple fields, by shifting telescope pointing directions with high cadence. This can significantly increase the effective sky-coverage of a blind survey on short time scales. For events associated with gamma-ray counterparts, the valuable early-time data can supplement high-energy observations. Regardless of gamma-ray association, detections can potentially be used to explore various phenomena associated with GRBs, such as orphan afterglows; dirty fireballs; and choked jets. We simulate a sample of GRBs and their respective optical signals at early times. After accounting for dynamic cadence, the light curves are given as input to a machine learning classifier, used to identify astrophysical transients. We find that by dedicating half of a LAST array to a blind search, one would expect to discover 7-11 GRBs per year, corresponding to an approximate intrinsic event rate of 0.12 per square degree per year.

We report the discovery of X-ray polarization from the X-ray-bright thread/filament G0.13-0.11 in the Galactic Center region. This filament features a bright hard X-ray source, most plausibly a Pulsar Wind Nebula (PWN), and an extended and structured diffuse component. Combining the polarization signal from IXPE with the imaging/spectroscopic data from Chandra, we find that X-ray emission of G0.13-0.11 is highly polarized PD=$57(\pm18)$\% in the 3-6 keV band, while the polarization angle is PA=$21^\circ(\pm9^\circ)$. This high degree of polarization proves the synchrotron origin of the X-ray emission from G0.13-0.11. In turn, the measured polarization angle implies that the X-ray emission is polarized approximately perpendicular to a sequence of non-thermal radio filaments that may be part of the Galactic Center Radio Arc. The magnetic field of the order of $100\,{\rm\mu G}$ appears to be preferentially ordered along the filaments. The above field strength is the fiducial value that makes our model self-consistent, while the other conclusions are largely model-independent.

The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$, $w$, $r$, $i$, $z$ and $y$ filters from Pan-STARRS as part of the Young Supernova Experiment, spanning $\sim$5,000 days. We find no significant detections in the Pan-STARRS pre-explosion light curve. We train a multilayer perceptron neural network to classify pre-supernova outbursts. We find no evidence of eruptive pre-supernova activity to a limiting absolute magnitude of $-7$. The limiting magnitudes from the full set of $gwrizy$ (average absolute magnitude $\approx$-8) data are consistent with previous pre-explosion studies. We use deep photometry from the literature to constrain the progenitor of SN2023ixf, finding that these data are consistent with a dusty red supergiant (RSG) progenitor with luminosity $\log\left(L/L_\odot\right)$$\approx$5.12 and temperature $\approx$3950K, corresponding to a mass of 14-20 M$_\odot$

In turbulent magnetized plasmas, charged particles can be accelerated to high energies through their interactions with the turbulent motions. As they do so, they draw energy from the turbulence, possibly up to the point where they start modifying the turbulent cascade. Stochastic acceleration then enters a nonlinear regime because turbulence damping back-reacts in turn on the acceleration process. This article develops a phenomenological model to explore this phenomenon and its consequences on the particle and turbulent energy spectra. We determine a criterion that specifies the threshold of nonthermal particle energy density and the characteristic momentum beyond which back-reaction becomes effective. Once the back-reaction sets in, the turbulence cascade becomes damped below a length scale that keeps increasing in time. The accelerated particle momentum distribution develops a near power-law of the form ${\rm d}n/{\rm d}p\propto p^{-s}$ with $s\sim2$ beyond the momentum at which back-reaction first sets in. At very high energies, where the gyroradius of accelerated particles becomes comparable to the outer scale of the turbulence, the energy spectrum can display an even harder spectrum with $s\sim 1.3-1.5$ over a short segment. The low-energy part of the spectrum, below the critical momentum, is expected to be hard ($s\sim 1$ or harder), and shaped by any residual acceleration process in the damped region of the turbulence cascade. This characteristic broken power-law shape with $s\sim 2$ at high energies may find phenomenological applications in various high-energy astrophysical contexts.

We present results from the Very Long Baseline Array (VLBA) multi-frequency (1.6, 4.4, 8.6, 22 GHz), high-sensitivity (~25 microJy beam^-1), sub-parsec scale (<1 pc) observations and Spectral Energy Distributions (SEDs) for a sample of 12 local active galactic nuclei (AGNs), a subset from our previous volume-complete sample with hard X-ray (14-195 keV) luminosities above 10^42 erg s^-1, out to a distance of 40 Mpc. All 12 of the sources presented here were detected in the C (4.4 GHz) and X (8.6 GHz) bands, 75% in the L band(1.6 GHz), and 50% in the K band (22 GHz). Most sources showed compact, resolved/slightly resolved, central sub-parsec scale radio morphology, except a few with extended outflow-like features. A couple of sources have an additional component that may indicate the presence of a dual-core, single or double-sided jet or a more intricate feature, such as radio emission resulting from interaction with nearby ISM. The spectral slopes are mostly GHz-peaked or curved, with a few showing steep, flat, or inverted spectra. We found that in the sub-parsec scale, the GHz-peaked spectra belong to the low-accreting, radio-loud AGNs with a tendency to produce strong outflows, possibly small-scale jet, and/or have a coronal origin. In contrast, flat/inverted spectra suggest compact radio emission from highly-accreting AGNs' central region, possibly associated with radio-quiet AGNs producing winds/shocks or nuclear star formation in the vicinity of black holes.

Interplanetary coronal mass ejections (ICMEs) are defined as ''coherent'' if they are capable of responding to external perturbations in a collective manner. This implies that information must be able to propagate across ICME structures, and if this is not the case, single-point in-situ measurements cannot be considered as indicative of global ICME properties. Here, we investigate the role of Alfvenic fluctuations (AFs) as mediators of ICME coherence. We consider multi-point magnetic field and plasma measurements of 10 ICMEs observed by the ACE and Wind spacecraft at 1 au at longitudinal separations of 0.5{\deg}-0.7{\deg}. For each event, we analyze the Alfvenicity in terms of the residual energy and cross helicity of fluctuations, and the coherence in terms of the magnetic correlation between Wind and ACE. We find that ~65% and 90% of ICME sheaths and magnetic ejecta (MEs), respectively, present extended AFs covering at least 20% of the structure. Cross helicity suggests AFs of solar and interplanetary origin may co-exist in the ICME population at 1 au. AFs are mainly concentrated downstream of shocks and in the back of MEs. The magnetic field is poorly correlated within sheaths, while the correlation decreases from the front to the back of the MEs for most magnetic field components. AFs are also associated with lower magnetic field correlations. This suggests either that ICME coherence is not mediated by Alfven waves, implying that the coherence scale may be smaller than previously predicted, or that the magnetic field correlation is not a measure of coherence.

We estimate the redshift-dependent, anisotropic clustering signal in DESI's Year 1 Survey created by tidal alignments of Luminous Red Galaxies (LRGs) and a selection-induced galaxy orientation bias. To this end, we measured the correlation between LRG shapes and the tidal field with DESI's Year 1 redshifts, as traced by LRGs and Emission-Line Galaxies (ELGs). We also estimate the galaxy orientation bias of LRGs caused by DESI's aperture-based selection, and find it to increase by a factor of seven between redshifts 0.4 - 1.1 due to redder, fainter galaxies falling closer to DESI's imaging selection cuts. These effects combine to dampen measurements of the quadrupole of the correlation function caused by structure growth on scales of 10 - 80 Mpc/h by about 0.15% for low redshifts (0.4<z<0.6) and 0.8% for high (0.8<z<1.1). We provide estimates of the quadrupole signal created by intrinsic alignments that can be used to correct this effect, which is necessary to meet DESI's forecasted precision on measuring the growth rate of structure. While imaging quality varies across DESI's footprint, we find no significant difference in this effect between imaging regions in the Legacy Imaging Survey.

The redshift evolution of bars is an important signpost of the dynamic maturity of disk galaxies. To characterize the intrinsic evolution safe from band-shifting effects, it is necessary to gauge how bar properties vary locally as a function of wavelength. We investigate bar properties in 16 nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G) at ultraviolet, optical and mid-infrared wavebands. Based on the ellipticity and position angle profiles from fitting elliptical isophotes to the two-dimensional light distribution, we find that both bar length and ellipticity - the latter often used as a proxy for bar strength - increase at bluer wavebands. Bars are 9% longer in the B-band than at 3.6 um. Their ellipticity increases typically by 8% in the B-band, with a significant fraction (>40%) displaying an increase up to 35%. We attribute the increase in bar length to the presence of star forming knots at the end of bars: these regions are brighter in bluer bands, stretching the bar signature further out. The increase in bar ellipticity could be driven by the apparent bulge size: the bulge is less prominent at bluer bands, allowing for thinner ellipses within the bar region. Alternatively, it could be due to younger stellar populations associated to the bar. The resulting effect is that bars appear longer and thinner at bluer wavebands. This indicates that band-shifting effects are significant and need to be corrected for high-redshift studies to reliably gauge any intrinsic evolution of the bar properties with redshift.

Inflation, a period of exponential expansion in the early Universe, is considered an important part of the standard $\Lambda$CDM cosmological model, and plays a crucial role in explaining a wide range of current observations. The standard inflationary model predicts a primordial spectrum of fluctuations that is nearly scale-independent, fitting remarkably well the latest observational data. Nevertheless, there is an ongoing discussion surrounding the transition from an initial homogeneous and isotropic quantum state, characterizing the matter fields during inflation, to a classical inhomogeneous/anisotropic one, which gives rise to large-scale structure in the Universe. To tackle this issue, in the present work we explore an inflationary scenario where quantum ``collapse'' (or reduction) occurs naturally during the evolution of the system; this model is inspired in the so called Continuous Spontaneous Localization (CSL) model. Our present work builds upon previous results by considering the primordial power spectrum up to the second order in the Hubble Flow Functions, where we perform an estimation of the model free parameters. By validating the predictions of the model against observational data, we investigate whether this second-order calculation can explain the slight departure from the power law observed in the scalar spectral running index. We hope this research contributes to the understanding of the quantum-to-classical transition and its implications for cosmology.

For the first time we analyse gravitational-wave strain data using waveforms constructed from strong gravity simulations of cosmic string loops collapsing to Schwarzschild black holes; a previously unconsidered source. Since the expected signal is dominated by a black-hole ringdown, it can mimic the observed gravitational waves from high-mass binary black hole mergers. To illustrate this, we consider GW190521, a short duration gravitational-wave event observed in the third LIGO--Virgo--KAGRA observing run. We show that describing this event as a collapsing cosmic string loop is favoured over previous cosmic string analyses by an approximate log Bayes factor of $22$. The binary black hole hypothesis is still preferred, mostly because the cosmic string remnant is non-spinning. It remains an open question whether a spinning remnant could form from loops with angular momentum, but if possible, it would likely bring into contention the binary black hole preference. Finally, we suggest that searches for ringdown-only waveforms would be a viable approach for identifying collapsing cosmic string events. This work opens up an important new direction for the cosmic-string and gravitational-wave communities.

We study the effects of the first nontrivial hexaquark, $d^*$(2380), on the equation of state of dense neutron star matter and investigate the consequences of its existence for neutron stars. The matter in the core regions of neutron stars is described using density-dependent relativistic mean-field theory. Our results show that within the parameter spaces examined in our paper, (i) the critical density at which the $d^*$ condensate emerges lies between 4 and 5 times the nuclear saturation density, (ii) $d^*$ hexaquarks are found to exist only in rather massive neutron stars, (iii) only relatively small fractions of the matter in the core of a massive neutron star may contain hexaquarks.

Dark matter may consist of axion-like particles (ALPs). When polarized electromagnetic radiation passes through the dark-matter media, interaction with background ALPs affects the polarization of photons. The condensate of axionic dark matter experiences periodic oscillations, and the period of the oscillations is of order of years for ultra-light dark matter. This would result in observable periodic changes in the polarization plane, determined by the phases of the ALP field at the Earth and at the source. In this paper, we use recent polarimetric observations of the HL Tauri protoplanetary disk performed in different years to demonstrate the lack of changes of polarization angles, and hence to constrain masses and photon couplings of the hypothetical axion-like ultralight dark matter.

The discovery of gravitational waves (GWs) opens a new window for exploring the physics of the early universe. Identifying the source of GWs and their spectra today turn out to be the important tasks so as to assist the experimental detection of stochastic GWs. In this paper, we investigate the oscillations of the ultralight dark photon (ULDP) into GWs in the dark halo. Assuming dark matter is composed of the ULDP and there are primordial dark magnetic fields (PDMFs) arising from the axion inflation and/or the dark phase transition, then the ULDP can oscillate into the GW when it passes through an environment of PDMFs. We derive the local energy density of GWs in the galaxy cluster induced by the instaneous oscillation of ULDP in the PDMFs. These stochastic local GWs exhibit a pulse-like spectrum, with frequency depending on the mass of the ULDP, and can be detected in Pulsar Timing Arrays (PTAs) or future space-based interferometers. We also find that the low-frequency GW signal observed by the NANOGrav collaboration and other PTA experiments can be explained by the oscillation of the ULDP in the PDMFs in the early universe.

Model-based approaches to imaging, like specialized image enhancements in astronomy, favour physics-based models which facilitate explanations of relationships between observed inputs and computed outputs. While this paper features a tutorial example, inspired by exoplanet imaging, that reveals embedded 2D fast Fourier transforms in an image enhancement model, the work is actually about the tensor algebra and software, or tensor frameworks, available for model-based imaging. The paper proposes a Ricci-notation tensor (RT) framework, comprising an extended Ricci notation, which aligns well with the symbolic dual-index algebra of non-Euclidean geometry, and codesigned object-oriented software, called the RTToolbox for MATLAB. Extensions offer novel representations for entrywise, pagewise, and broadcasting operations popular in extended matrix-vector (EMV) frameworks for imaging. Complementing the EMV algebra computable with MATLAB, the RTToolbox demonstrates programmatic and computational efficiency thanks to careful design of tensor and dual-index classes. Compared to a numeric tensor predecessor, the RT framework enables superior ways to model imaging problems and, thereby, to develop solutions.

The experiments were performed at bremsstrahlung end-point energies of 10-23 MeV with the beam from the MT-25 microtron with the use of the {\gamma}-activation technique. The experimental values of relative yields were compared with theoretical results obtained on the basis of TALYS with the standard parameters and the combined model of photonucleon reactions. Including isospin splitting in the combined model of photonucleon reactions allows to describe experimental data on reactions with proton escape in energies range from 10 to 23 MeV. Therefore, taking into account isospin splitting is necessary for a correct description of the decay of the GDR.

We develop the Brans-Dicke theory of gravity in the context of varying constants of Nature. Using the unimodular formalism of General Relativity, we create a platform to provide physical relational times giving the evolution of physical constants. We therefore review the ideas and experiments behind varying constants, mostly focusing on the speed of light and the gravitational constant. Then, we apply this idea to the energy conservation in cosmology, illustrating the arising patterns. Motivated by a varying gravitational constant resulting from Mach's principle, we develop the unimodular formalism of varying constants in the Brans-Dicke theory. Doing so, we obtain several original results, some of which can be compared with phenomenological observation. Finally, we suggest how a varying Brans-Dicke parameter could be linked to the Cosmological Constant problem.

In this study, we delve into the observational implications of rotating Loop Quantum Black Holes (LQBHs) within an astrophysical framework. We employ semi-analytical General Relativistic Radiative Transfer (GRRT) computations to study the emission from the accretion flow around LQBHs. Our findings indicate that the increase of Loop Quantum Gravity (LQG) effects results in an enlargement of the rings from LQBHs, thereby causing a more circular polarization pattern in the shadow images. We make comparisons with the Event Horizon Telescope (EHT) observations of Sgr\,A$^*$ and M\,87$^*$, which enable us to determine an upper limit for the polymetric function $P$ in LQG. The upper limit for Sgr\,A$^*$ is $0.2$, while for M\,87$^*$ it is $0.07$. Both black holes exhibit a preference for a relatively high spin ($a\gtrsim0.5$ for Sgr\,A$^*$ and $0.5\lesssim a \lesssim 0.7$ for M\,87$^*$). The constraints for Sgr\,A$^*$ are based on black hole spin and ring diameter, whereas for M\,87$^*$, the constraints are further tightened by the polarimetric pattern. In essence, our simulations provide observational constraints on the effect of LQG in supermassive black holes (SMBH), providing the most consistent comparison with observation.

Observations, mainly of outbursts in dwarf novae, imply that the anomalous viscosity in highly ionized accretion discs is magnetic in origin, and requires that the plasma $\beta \sim 1$. Until now most simulations of the magnetic dynamo in accretion discs have used a local approximation (known as the shearing box). While these simulations demonstrate the possibility of a self-sustaining dynamo, the magnetic activity generated in these models saturates at $\beta \gg 1$. This long-standing discrepancy has previously been attributed to the local approximation itself. There have been recent attempts at simulating magnetic activity in global accretion discs with parameters relevant to the dwarf novae. These too find values of $\beta \gg 1$. We speculate that the tension between these models and the observations may be caused by numerical magnetic diffusivity. As a pedagogical example, we present exact time-dependent solutions for the evolution of weak magnetic fields in an incompressible fluid subject to linear shear and magnetic diffusivity. We find that the maximum factor by which the initial magnetic energy can be increased depends on the magnetic Reynolds number as ${\mathcal R}_{\rm m}^{2/3}$. We estimate that current global numerical simulations of dwarf nova discs have numerical magnetic Reynolds numbers around 6 orders of magnitude less than the physical value found in dwarf nova discs of ${\mathcal R}_{\rm m} \sim 10^{10}$. We suggest that, given the current limitations on computing power, expecting to be able to compute realistic dynamo action in observable accretion discs using numerical MHD is, for the time being, a step too far.

We develop a method for accurately calculating vacuum decay rates beyond the thin-wall regime in a pure scalar field theory at the one-loop level of the effective action. It accounts for radiative effects resulting from quantum corrections to the classical bounce, including gradient effects stemming from the inhomogeneity of the bounce background. To achieve this, it is necessary to compute not only the functional determinant of the fluctuation operator in the background of the classical bounce but also its functional derivative evaluated at the classical bounce. The former is efficiently calculated using the Gel'fand-Yaglom method. We illustrate how the latter can also be calculated with the same method, combined with a computation of various Green's functions.