Papers for Wednesday, Mar 22 2023

The stellar companion to the weak-line T Tauri star DI Tau A was first discovered by the lunar occultation technique in 1989 and was subsequently confirmed by a speckle imaging observation in 1991. It has not been detected since, despite being targeted by five different studies that used a variety of methods and spanned more than 20 years. Here, we report the serendipitous rediscovery of DI Tau B during our Young Exoplanets Spectroscopic Survey (YESS). Using radial velocity data from YESS spanning 17 years, new adaptive optics observations from Keck II, and a variety of other data from the literature, we derive a preliminary orbital solution for the system that effectively explains the detection and (almost all of the) non-detection history of DI Tau B. We estimate the dynamical masses of both components, finding that the large mass difference (q $\sim$0.17) and long orbital period ($\gtrsim$35 years) make DI Tau system a noteworthy and valuable addition to studies of stellar evolution and pre-main-sequence models. With a long orbital period and a small flux ratio (f2/f1) between DI Tau A and B, additional measurements are needed for a better comparison between these observational results and pre-main-sequence models. Finally, we report an average surface magnetic field strength ($\bar B$) for DI Tau A, of $\sim$0.55 kG, which is unusually low in the context of young active stars.

Using color-magnitude diagrams from deep archival Hubble Space Telescope imaging, we self-consistently measure the star formation history of Eridanus II (Eri II), the lowest-mass galaxy ($M_{\star}(z=0) \sim 10^5 M_{\odot}$) known to host a globular cluster (GC), and the age, mass, and metallicity of its GC. The GC ($\sim13.2\pm0.4$ Gyr, $\langle$[Fe/H]$\rangle = -2.75\pm0.2$ dex) and field (mean age $\sim13.5\pm0.3$ Gyr, $\langle$[Fe/H]$\rangle = -2.6\pm0.15$ dex) have similar ages and metallicities. Both are reionization-era relics that formed before the peak of cosmic star and GC formation ($z\sim2-4$). The ancient star formation properties of Eri II are not extreme and appear similar to $z=0$ dwarf galaxies. We find that the GC was $\lesssim4$ times more massive at birth than today and was $\sim$10% of the galaxy's stellar mass at birth. At formation, we estimate that the progenitor of Eri II and its GC had $M_{\rm UV} \sim -7$ to $-12$, making it one of the most common type of galaxy in the early Universe, though it is fainter than direct detection limits, absent gravitational lensing. Archaeological studies of GCs in nearby low-mass galaxies may be the only way to constrain GC formation in such low-mass systems. We discuss the strengths and limitations in comparing archaeological and high redshift studies of cluster formation, including challenges stemming from the Hubble Tension, which introduces uncertainties into the mapping between age and redshift.

We report the discovery of a low-mass $z=5.200\pm 0.002$ galaxy that is in the process of ceasing its star formation. The galaxy, MACS0417-z5PSB, is multiply imaged with magnification factors $\sim40$ by the galaxy cluster MACS J0417.5-1154, observed as part of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). Using observations of MACS0417-z5PSB with a JWST/NIRSpec Prism spectrum and NIRCam imaging, we investigate the mechanism responsible for the cessation of star formation of the galaxy, and speculate about possibilities for its future. Using spectrophotometric fitting, we find a remarkably low stellar mass of $\rm{M_*}=4.3\pm^{0.9}_{0.8} \times 10^{7} \rm{M_{\odot}}$, less than 1% of the characteristic stellar mass at $z\sim5$. We measure a de-lensed rest-UV half-light radius in the source plane of $30\pm^{7}_{5}$ pc, and measure a star formation rate from H$\alpha$ of $0.14\pm^{0.17}_{0.12}$ $\rm{M_{\odot}/yr}$. We find that under the assumption of a double power law star formation history, MACS0417-z5PSB has seen a recent rise in star formation, peaking $\sim10-30$ Myr ago and declining precipitously since then. Together, these measurements reveal a low-mass, extremely compact galaxy which is in the process of ceasing star formation. We investigate the possibilities of mechanisms that have led to the cessation of star formation in MACS0417-z5PSB, considering stellar and AGN feedback, and environmental processes. We can likely rule out an AGN and most environmental processes, but leave open the possibility that MACS0417-z5PSB could be a star forming galaxy in the lull of a bursty star formation history.

The inner cores of massive neutron stars contain strongly interacting matter at the highest densities reached in our Universe. Under these conditions the cores may undergo a phase transition to deconfined quark matter, which exhibits approximate conformal symmetry. Using a Bayesian inference setup that utilizes all available neutron-star measurements and state-of-the-art theoretical calculations, we demonstrate that in the cores of the most massive stars the equation of state is consistent with the presence of deconfined quark matter. We do this by (i) establishing an effective conformal symmetry restoration with 88% credence at the highest densities probed in these objects, and (ii) demonstrating that the number of active degrees of freedom favors an interpretation of this finding in terms of the presence of deconfined matter. The remaining probability for purely hadronic maximal-mass stars arises from equation-of-state behavior featuring small sound-speed and polytropic-index values, consistent with a first-order phase transition.

Compact Symmetric Objects (CSOs) are jetted Active Galactic Nuclei (AGN) with overall projected size < 1 kpc. The classification was introduced to distinguish these objects from the majority of known compact jetted-AGN, where the observed emission is relativistically boosted towards the observer. The original classification criteria for CSOs were: (i) evidence of emission on both sides of the center of activity, and (ii) overall size < 1 kpc. However some relativistically boosted objects with jet axes close to the line of sight appear symmetric and have been mis-classified as CSOs in the literature, thereby undermining the CSO classification. We introduce two new CSO classification criteria based on (iii) flux density variability, and (iv) the apparent velocity of components moving along the jets. As a first step towards creating a comprehensive catalog of "bona fide" CSOs, in this paper we identify 79 bona fide CSOs that meet our expanded CSO selection criteria. This sample of bona fide CSOs can be used for astrophysical studies of CSOs without fear of contamination by objects incorrectly identified as CSOs. We define three complete sub-samples of the 79 CSOs, which are suitable for statistical tests, and show that the fraction of CSOs in flux density limited samples with $\rm S_{5~GHz}$ > 700 mJy is between 6.4% and 8.2%.

Several studies have shown that a number of stars pulsating in p modes lie between the $\beta$ Cep and $\delta$ Sct instability strips in the Hertzsprung-Russell (HR) Diagram. At present, there is no certain understanding of how p~modes can be excited in this $T_{\rm eff}$ range. The goal of this work is to disprove the conjecture that all stars pulsating in p modes and lying in this $T_{\rm eff}$ range are the result of incorrect measurements of $T_{\rm eff}$, contamination, or the presence of unseen cooler companions lying in the $\delta$ Sct instability strip (given the high binary fraction of stars in this region of the HR Diagram). Using TESS data, we show that the A0Vnne star HD 42477 has a single p mode coupled to several r modes and/or g modes. We rule out a contaminating background star with a pixel-by-pixel examination, and we essentially rule out the possibility of a companion $\delta$ Sct star in a binary. We model the pulsations in HD 42477 and suggest that the g modes are excited by overstable convective core modes. We also conjecture that the single p mode is driven by coupling with the g modes, or that the oblateness of this rapidly-rotating star permits driving by He II ionization in the equatorial region.

Compact Symmetric Objects (CSOs) are a class of compact, jetted Active Galactic Nuclei (AGN) whose jet axes are not aligned close to the line of sight, and whose observed emission is not predominantly relativistically boosted towards us. Using complete samples of CSOs, we present three independent lines of evidence, based on their relative numbers, their redshift distributions, and their size distributions, which show conclusively that most CSOs do not evolve into larger-scale radio sources. Thus CSOs belong to a distinct population of jetted-AGN. This population should be characterized as "short-lived", as opposed to "young". We show that there is a sharp upper cutoff in the CSO size distribution at $\approx$ 500 pc, which cannot result from random episodic fueling events. There is clearly something that limits the fueling to $\lesssim 100 M_\odot$. Possible origins of CSOs, if not related to the fueling, must be related to the accretion disk, or the collimation of the relativistic jets. CSOs may well have a variety of origins, with each of the above mechanisms producing subsets of CSOs. Whatever the physical mechanism(s) might be, the distinct differences between CSOs and other jetted-AGN provide crucial insights into the formation and evolution of relativistic jets in AGN and the supermassive black holes that drive them.

According to the $\Lambda$CDM cosmology, present-day massive galaxies should contain a sizable fraction of dark matter within their stellar body. Models indicate that in massive early-type galaxies (ETGs) with $M_\star\approx1.5\times10^{11}\,{\rm M}_\odot$ dark matter should account for $\sim60\%$ of the dynamical mass within five effective radii ($5\,R_{\rm e}$). Most massive ETGs have been shaped through a two-phase process: the rapid growth of a compact core was followed by the accretion of an extended envelope through mergers. The exceedingly rare galaxies that have avoided the second phase, the so-called relic galaxies, are thought to be the frozen remains of the massive ETG population at $z\gtrsim2$. The best relic galaxy candidate discovered to date is NGC 1277, in the Perseus cluster. We used deep integral field GCMS data to revisit NGC 1277 out to an unprecedented radius of 6 kpc (corresponding to $5\,R_{\rm e}$). By using Jeans anisotropic models we find a negligible dark matter fraction within $5\,R_{\rm e}$ ($f_{\rm DM}(5\,R_{\rm e})<0.07$; two-sigma confidence level), which is in strong tension with the $\Lambda$CDM expectation. Since the lack of an extended envelope would reduce dynamical friction and prevent the accretion of an envelope, we propose that NGC 1277 lost its dark matter very early or that it was dark matter deficient ab initio. We discuss our discovery in the framework of recent proposals suggesting that some relic galaxies may result from dark matter stripping within galaxy clusters. Alternatively, NGC 1277 might have been born in a high-velocity collision of gas-rich proto-galactic fragments, where dark matter left behind a disc of dissipative baryons. We speculate that the relative velocities of $\sim2000\,{\rm km\,s^{-1}}$ required for the latter process to happen were possible in the progenitors of the present-day rich galaxy clusters.

Compact Symmetric Objects (CSOs) form a distinct class of jetted active galactic nuclei (jetted-AGN). We examine a carefully selected sample of 54 CSOs, and confirm that there are two unrelated classes: an edge-dimmed, low-luminosity class (CSO 1), and an edge-brightened, high-luminosity class (CSO 2). Using statistically significant blind tests, we show that CSO 2s themselves consist of two morphologically distinct classes: CSO 2.0, having prominent hot-spots at the leading edges of narrow jets and/or narrow lobes; and CSO 2.2, without prominent hot-spots, and with broad jets and/or lobes. An intermediate class, CSO 2.1, exhibits mixed properties. The four classes occupy different, overlapping, portions of the luminosity-size plane, with the sizes of largest CSOs being $\sim 500$pc. We advance the hypothesis that CSO 2.0s are young and evolve through CSO 2.1s into CSO 2.2s, which are old (up to $\sim 5000$ yr). Thus CSOs do not evolve into larger types of jetted-AGN, but spend their whole life cycle as CSOs. The radio emission region energies in the CSO 2s we have studied range from $\sim 10^{-4}\, M_\odot {c}^2$ to $\sim 7 \, M_\odot {c}^2$. We show that the transient nature of CSO 2s, and their birthrate, can be explained through ignition in the tidal disruption events of giant stars. We also consider possibilities for tapping the spin energy of the supermassive black hole, and tapping the energy of the accretion disk, in a manner similar to, but not the same as, that which occurs in dwarf novae. Our results demonstrate conclusively that CSOs constitute a large family of AGN in which we have thus far studied only the brightest. More comprehensive radio studies, with higher sensitivity, resolution, and dynamic range, will revolutionize our understanding of AGN and the central engines that power them.

We study the co-evolution of dark matter halos, galaxies and supermassive black holes using an empirical galaxy evolution model from $z=0$ -- $10$. We demonstrate that by connecting dark matter structure evolution with simple empirical prescriptions for baryonic processes, we are able to faithfully reproduce key observations in the relation between galaxies and their supermassive black holes. By assuming a physically-motivated, direct relationship between the galaxy and supermassive black hole properties to the mass of their host halo, we construct expressions for the galaxy stellar mass function, galaxy UV luminosity function, active black hole mass function and quasar bolometric luminosity function. We calibrate the baryonic prescriptions using a fully Bayesian approach in order to reproduce observed population statistics. The obtained parametrizations are then used to study the relation between galaxy and black hole properties, as well as their evolution with redshift. The galaxy stellar mass -- UV luminosity relation, black hole mass -- stellar mass relation, black hole mass -- AGN luminosity relation, and redshift evolution of these quantities obtained from the model are qualitatively consistent with observations. Based on these results, we present upper limits on the expected number of sources for $z=5$ up to $z=15$ for scheduled JWST and \textit{Euclid} surveys, thus showcasing that empirical models can offer qualitative as well as quantitative prediction in a fast, easy and flexible manner that complements more computationally expensive approaches.

Dynamical stellar-evolution modeling through the AGB phase reveals that radial pulsations with very fast-growing amplitudes develop if the luminosity to mass ratio of stars with tenuous envelopes exceeds a critical limit. An instability going nonlinear already after a few pulsation cycles might qualify as a source of the superwind - postulated to shed a substantial part of a star's envelope over a very short time - of hitherto persistently mysterious nature.

The abundance of the short-lived radioisotopes 26-Al and 60-Fe in the early Solar system is usually explained by the Sun either forming from pre-enriched material, or the Sun's protosolar disc being polluted by a nearby supernova explosion from a massive star. Both hypotheses suffer from significant drawbacks: the former does not account for the dynamical evolution of star-forming regions, while in the latter the time for massive stars to explode as supernovae can be similar to, or even longer than, the lifetime of protoplanetary discs. In this paper, we extend the disc enrichment scenario to include the contribution of 26-Al from the winds of massive stars before they explode as supernovae. We use N-body simulations and a post-processing analysis to calculate the amount of enrichment in each disc, and we vary the stellar density of the star-forming regions. We find that stellar winds contribute to disc enrichment to such an extent that the Solar system's 26-Al/60-Fe ratio is reproduced in up to 50 per cent of discs in dense (rho = 1000Msun pc^-3) star-forming regions. When winds are a significant contributor to the SLR enrichment, we find that Solar system levels of enrichment can occur much earlier (before 2.5 Myr) than when enrichment occurs from supernovae, which start to explode at later ages (>4 Myr). We find that Solar system levels of enrichment all but disappear in low-density star-forming regions (rho < 10Msun pc^-3), implying that the Solar system must have formed in a dense, populous star-forming region if 26-Al and 60-Fe were delivered directly to the protosolar disc from massive-star winds and supernovae.

The efficiency of particle acceleration at shock waves in relativistic, magnetized astrophysical outflows is a debated topic with far-reaching implications. Here, for the first time, we study the impact of turbulence in the pre-shock plasma. Our simulations demonstrate that, for a mildly relativistic, magnetized pair shock (Lorentz factor $\gamma_{\rm sh} \simeq 2.7$, magnetization level $\sigma \simeq 0.01$), strong turbulence can revive particle acceleration in a superluminal configuration that otherwise prohibits it. Depending on the initial plasma temperature and magnetization, stochastic-shock-drift or diffusive-type acceleration governs particle energization, producing powerlaw spectra $\mathrm{d}N/\mathrm{d}\gamma \propto \gamma^{-s}$ with $s \sim 2.5-3.5$. At larger magnetization levels, stochastic acceleration within the pre-shock turbulence becomes competitive and can even take over shock acceleration.

We conduct a statistical analysis of the factors responsible for the variation in the ionization parameter (U) of high-redshift star-forming galaxies based on medium resolution JWST/NIRSpec observations obtained by the Cosmic Evolution Early Release Science (CEERS) survey. The sample consists of 48 galaxies with spectroscopic redshifts z=2.7-6.3 which are largely representative of typical star-forming galaxies at these redshifts. The [SII] 6718, 6733 doublet is used to estimate electron densities (n_e), and dust-corrected Ha luminosities are used to compute total ionizing photon rates (Q). Using composite spectra of galaxies in bins of [OIII]/[OII] (i.e., O32) as a proxy for U, we determine that galaxies with higher O32 have <n_e> ~ 500 cm^-3 that are at least a factor of ~5 larger than that of lower-O32 galaxies. We do not find a significant difference in <Q> between low- and high-O32 galaxies. Photoionization modeling of all available strong rest-frame optical emission lines is used to simultaneously constrain U and oxygen abundance (Z_neb). We find a large spread in log U of ~1.5 dex at a fixed Z_neb. On the other hand, the data indicate a highly significant correlation between U and star-formation-rate surface density (Sigma_SFR) which appears to be redshift invariant at z~1.6-6.3, and possibly up to z~9.5. We consider several avenues through which metallicity and Sigma_SFR (or gas density) may influence U, including variations in n_e and Q that are tied to metallicity and gas density, internal dust extinction of ionizing photons, and the effects of gas density on the volume filling fraction of dense clumps in HII regions and the escape fraction of ionizing photons. Based on these considerations, we conclude that gas density may play a more central role than metallicity in modulating U at these redshifts.

We investigate the link between the bar rotation rate and dark matter content in barred galaxies by concentrating on the cases of the lenticular galaxies NGC4264 and NGC4277. These two gas-poor galaxies have similar morphologies, sizes, and luminosities. But, NGC4264 hosts a fast bar, which extends to nearly the corotation, while the bar embedded in NGC4277 is slow and falls short of corotation. We derive the fraction of dark matter $f_{\rm DM, bar}$ within the bar region from Jeans axisymmetric dynamical models by matching the stellar kinematics obtained with the MUSE integral-field spectrograph and using SDSS images to recover the stellar mass distribution. We build mass-follows-light models as well as mass models with a spherical halo of dark matter, which is not tied to the stars. We find that the inner regions of NGC4277 host a larger fraction of dark matter ($f_{\rm DM, bar} = 0.53 \pm 0.02$) with respect to NGC4264 ($f_{\rm DM, bar} = 0.33 \pm 0.04$) in agreement with the predictions of theoretical works and the findings of numerical simulations, which have found that fast bars live in baryon-dominated discs, whereas slow bars experienced a strong drag from the dynamical friction due to a dense DM halo. This is the first time that the bar rotation rate is coupled to $f_{\rm DM, bar}$ derived from dynamical modelling.

We investigate the distribution of the lithium abundances, A(Li), of metal-poor dwarf and subgiant stars within the limits 5500 K < Teff < 6700 K, -6.0 < [Fe/H] < -1.5, and logg > ~3.5 (a superset of parameters first adopted by Spite and Spite), using literature data for some 200 stars. We address the problem of the several methods that yield Teff differences up to 350 K, and hence uncertainties of 0.3 dex in [Fe/H] and A(Li), by anchoring Teff to the Infrared Flux Method. We seek to understand the behaviour of A(Li) as a function of [Fe/H] -- small dispersion at highest [Fe/H], ``meltdown'' at intermediate values (i.e. large spread in Li below the Spite Plateau), and extreme variations at lowest [Fe/H]. Decreasing A(Li) is accompanied by increasing dispersion. Insofar as [Fe/H] increases as the universe ages, the behavior of A(Li) reflects chaotic star formation involving destruction of primordial Li, which settles to the classic Spite Plateau, with A(Li) ~2.3, by the time the Galactic halo reaches [Fe/H] ~ -3.0. We consider three phases: (1) first star formation in C-rich environments ([C/Fe] > 2.3), with depleted Li; (2) silicates-dominated star formation and destruction of primordial Li during pre-main-sequence evolution; and (3) materials from these two phases co-existing and coalescing to form C-rich stars with A(Li) below the Spite Plateau, leading to a toy model with the potential to explain the ``meltdown''. We comment on the results of Mucciarelli et al. on the Lower RGB, and the suggestion of Aguado et al. favouring a lower primordial lithium abundance than generally accepted.

We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $\mu$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity gradient unlike what was found in [CII] 158 $\mu$m at a larger scale, suggesting random bulk motion of ionized gas clouds inside the galaxy. In contrast, dust emission appears as two individual clumps apparently separating or bridging the [OIII]/stellar clumps. The cross correlation coefficient between dust and ultraviolet-related emission (i.e., [OIII] and ultraviolet continuum) is unity on a galactic scale, while it drops at < 1 kpc, suggesting well mixed geometry of multi-phase interstellar media on sub-kpc scales. If the cutoff scale characterizes different stages of star formation, the cutoff scale can be explained by gravitational instability of turbulent gas. We also report on a kpc-scale off-center cavity embedded in the dust continuum image. This could be a superbubble producing galactic-scale outflows, since the energy injection from the 4 Myr starburst suggested by a spectral energy distribution analysis is large enough to push the surrounding media creating a kpc-scale cavity.

The faint CO gases in debris disks are easily dissolved into C by UV irradiation, while CO can be reformed via reactions with hydrogen. The abundance ratio of C/CO could thus be a probe of the amount of hydrogen in the debris disks. We conduct radiative transfer calculations with chemical reactions for debris disks. For a typical dust-to-gas mass ratio of debris disks, CO formation proceeds without the involvement of H$_2$ because a small amount of dust grains makes H$_2$ formation inefficient. We find that the CO to C number density ratio depends on a combination of $n_\mathrm{H}Z^{0.4}\chi^{-1.1}$, where $n_\mathrm{H}$ is the hydrogen nucleus number density, $Z$ is the metallicity, and $\chi$ is the FUV flux normalized by the Habing flux. Using an analytic formula for the CO number density, we give constraints on the amount of hydrogen and metallicity for debris disks. CO formation is accelerated by excited H$_2$ either when the dust-to-gas mass ratio is increased or the energy barrier of chemisorption of hydrogen on the dust surface is decreased. This acceleration of CO formation occurs only when the shielding effects of CO are insignificant. In shielded regions, the CO fractions are almost independent of the parameters of dust grains.

The electromagnetic counterparts to gravitational wave (GW) merger events hold immense scientific value, but are difficult to detect due to the typically large localisation errors associated with GW events. The Low-Frequency Array (LOFAR) is an attractive GW follow-up instrument owing to its high sensitivity, large instantaneous field of view, and ability to automatically trigger on events to probe potential prompt emission within minutes. Here, we report on 144-MHz LOFAR observations of three GW merger events containing at least one neutron star that were detected during the third GW observing run (S190426c, S191213g and S200213t). While these events are not particularly significant, we use multi-epoch LOFAR data to devise a sensitive wide-field GW follow-up strategy to be used in future GW observing runs. In particular, we improve on our previously published strategy by implementing direction dependent calibration and mosaicing, resulting in nearly an order of magnitude increase in sensitivity and more uniform coverage. We achieve a uniform $5\sigma$ sensitivity of $870$ $\mu$Jy across a single instantaneous LOFAR pointing's 21 deg$^{2}$ core, and a median sensitivity of 1.1 mJy when including the full 89 deg$^{2}$ hexagonal beam pattern. We also place the deepest transient surface density limits yet on of order month timescales for surveys between 60--340 MHz (0.017 deg$^{-2}$ above $2.0$ mJy and 0.073 deg$^{-2}$ above $1.5$ mJy).

We describe the design, validation, and commissioning of a new correlator termed "MWAX" for the Murchison Widefield Array (MWA) low-frequency radio telescope. MWAX replaces an earlier generation MWA correlator, extending correlation capabilities and providing greater flexibility, scalability, and maintainability. MWAX is designed to exploit current and future Phase II/III upgrades to MWA infrastructure, most notably the simultaneous correlation of all 256 of the MWA's antenna tiles (and potentially more in future). MWAX is a fully software-programmable correlator based around an ethernet multicast architecture. At its core is a cluster of 24 high-performance GPU-enabled commercial-off-the-shelf compute servers that together process in real-time up to 24 coarse channels of 1.28 MHz bandwidth each. The system is highly flexible and scalable in terms of the number of antenna tiles and number of coarse channels to be correlated, and it offers a wide range of frequency / time resolution combinations to users. We conclude with a roadmap of future enhancements and extensions that we anticipate will be progressively rolled out over time.

It has been proposed that accretion-induced collapse (AIC) of massive white dwarfs (WDs) is an indispensable path for the formation of neutron star (NS) binaries. Although there are still no direct evidence for the existence of AIC events, several kinds of NS systems are suggested to originate from the AIC processes. One of the representative evidence is the detection of the strong magnetic field and slow spin NSs with ultra-light companions (<=0.1 Msun) in close orbits. However, previous studies cannot explain the formation of AIC events with such low-mass companions. In the present work, we evolve a series of ONe WD+He WD systems to the formation of AIC events (named as the He WD donor channel), and the NS binaries behave as ultra-compact X-ray binaries when the He WDs refill their Roche-lobes. We found that the ONe WD+He WD systems a possible channel for the formation of the newly formed NS+ultra-light companion systems just after the AIC event. Although there are some other inconsistent properties, the detected companion mass and orbital period of 4U 1626-67 (one of the newly formed NS binaries with ultra-light companions) can be reproduced by the He WD donor channel. In addition, combined with previous asteroseismology results, we speculate that an UCXB source (XTE J1751-305) may originate from the He WD donor channel.

We studied the broadband X-ray timing and spectral behaviors of the newly confirmed accreting millisecond X-ray pulsar MAXI J1816$-$195 during its 2022 outburst. We used the data from Insight-HXMT ME/HE, NICER and NuSTAR which cover the energy range between 0.8--210 keV. Coherent timing analysis of Insight-HXMT ME/HE and NICER data revealed a complex behavior at the early stage of the outburst, between MJD 59737.0-59741.9, and beyond between MJD 59741.9-59760.6 a significant spin-up $\dot{\nu}$ of $(9.0\pm2.1)\times10^{-14}~{\rm Hz~s^{-1}}$. We detected the X-ray pulsations up to $\sim95$ keV in a combination of Insight-HXMT HE observations. The pulse profiles were quite stable over the whole outburst and could be well described by a truncated Fourier series using two harmonics, the fundamental and the first overtone. Both components kept alignment in the range 0.8--64 keV. The joint and time-averaged NICER and Insight-HXMT spectra in the energy range 1--150 keV were well fitted by the absorbed Comptonization model compps plus disk blackbody with two additional Gaussian components. Based on the bolometric flux during the spin-up epoch, we determined a magnetic field strength of $(0.23-1.11)\times10^8$ G in MAXI J1816--195.

Magnetic fields (B-fields) play an important role in molecular cloud fragmentation and star formation, but are very difficult to detect. The temporal correlation between the field strength (B) and gas density (n) of an isolated cloud has been suggested as an indication of the dynamical importance of B-fields relative to self-gravity. This temporal B-n relation is, however, unobservable. What can be observed using Zeeman measurements are the "spatial B-n relations" from the current plane of the sky. Nevertheless, the temporal B-n relation argument has still been widely used to interpret observations. Here we present the first numerical test of the legitimacy of this interpretation. From a simulation that can reproduce the observed Zeeman spatial B~n^2/3 relation, we found that temporal B-n relations of individual cores bear no resemblance to the spatial B-n relations. This result inspired us to discover that the true mechanism behind the 2/3 index is random turbulence compression instead of symmetrical gravitational contraction.

With Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations, we detected highly variable helium emission lines from the optical counterpart of the supersoft ultraluminous X-ray source (ULX) NGC 247 ULX-1. No Balmer lines can be seen in the source spectrum. This is the first evidence for the presence of a helium donor star in ULXs, consistent with a prediction that helium donor stars may be popular in ULXs. The helium lines with a FWHM of about 200 km/s are likely produced on the outer accretion disk. Their strong variation implies that the central X-ray source can be significantly obscured to the outer disk. Also, a ring or a double-ring structure is revealed in the MUSE image. It is unknown whether or not it is related to the progenitor of the ULX binary.

The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a neutral hydrogen survey (HI) that is running on the Australian SKA Pathfinder (ASKAP), a precursor telescope for the Square Kilometre Array (SKA). The goal of WALLABY is to use ASKAP's powerful wide-field phased array feed technology to observe three quarters of the entire sky at the 21 cm neutral hydrogen line with an angular resolution of 30 arcseconds. Post-processing activities at the Australian SKA Regional Centre (AusSRC), Canadian Initiative for Radio Astronomy Data Analysis (CIRADA) and Spanish SKA Regional Centre prototype (SPSRC) will then produce publicly available advanced data products in the form of source catalogues, kinematic models and image cutouts, respectively. These advanced data products will be generated locally at each site and distributed across the network. Over the course of the full survey we expect to replicate data up to 10 MB per source detection, which could imply an ingestion of tens of GB to be consolidated in the other locations near real time. Here, we explore the use of an asymmetric database replication model and strategy, using PostgreSQL as the engine and Bucardo as the asynchronous replication service to enable robust multi-source pools operations with data products from WALLABY. This work would serve to evaluate this type of data distribution solution across globally distributed sites. Furthermore, a set of benchmarks have been developed to confirm that the deployed model is sufficient for future scalability and remote collaboration needs.

We compute the evolution and rotational periods of young stars, using the MESA code, starting from a stellar seed, and take protostellar accretion, stellar winds, and the magnetic star-disk interaction into account. Furthermore, we add a certain fraction of the energy of accreted material into the stellar interior as additional heat and combine the resulting effects on stellar evolution with the stellar spin model. For different combinations of parameters, stellar periods at an age of 1~Myr range between 0.6~days and 12.9~days. Thus, during the relatively short time period of 1~Myr, a significant amount of stellar angular momentum can already be removed by the interaction between the star and its accretion disk. The amount of additional heat added into the stellar interior, the accretion history, and the presence of disk and stellar winds have the strongest impact on the stellar spin evolution during the first million years. The slowest stellar rotations result from a combination of strong magnetic fields, a large amount of additional heat, and effective winds. The fastest rotators combine weak magnetic fields and ineffective winds or result from a small amount of additional heat added to the star. Scenarios that could lead to such configurations are discussed. Different initial rotation periods of the stellar seed, on the other hand, quickly converges and do not affect the stellar period at all. Our model matches up to 90\% of the observed rotation periods in six young ($\lesssim 3$~Myr) clusters. Based on these intriguing results, we motivate to combine our model with a hydrodynamic disk evolution code to self-consistently include several important aspects such as episodic accretion events, magnetic disk winds, internal, and external photo-evaporation. (Shortened)

We apply the full set of most update dynamical and geometrical data in cosmology to the nonextensive Barrow entropic holographic dark energy. We show that the data point towards an extensive Gibbs-like entropic behaviour for the cosmological horizons, which is the extreme case of the Barrow entropy, with the entropy parameter being $\Delta > 0.86$, close to the maximum threshold of $\Delta =1$ where the fractal dimension of the area-horizon becomes almost or just the volume and the intensivity is recovered. Futhermore, we find that the standard Bekenstein area-entropy limit ($\Delta = 0$) is excluded by the set of our data. This contradicts the bounds obtained recently from early universe tests such as the baryon asymmetry, the big-bang nucleosynthesis, and the inflation limiting $\Delta< 0.008$ at the most extreme case.

We use the Random Forest (RF) algorithm to develop a tool for automated activity classification of galaxies into 5 different classes: Star-forming (SF), AGN, LINER, Composite, and Passive. We train the algorithm on a combination of mid-IR (WISE) and optical photometric data while the true labels (activity classes) are based on emission line ratios. Our classifier is built to be redshift-agnostic and it is applicable to objects up to z $\sim$0.1. It reaches a completeness $>$80 % for SF and Passive galaxies, and $\sim$60 % for AGN. Applying it to an all-sky galaxy catalog (HECATE) reveals a large population of low-luminosity AGNs outside the AGN locus in the standard mid-IR diagnostics.

Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. The Mid-Infrared Instrument (MIRI) onboard JWST offers an excellent chance to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN population in the distant universe. We estimate the source properties of MIRI-selected objects by utilizing spectral energy distribution (SED) modelling, and classify them into star-forming galaxies (SF), SF-AGN mixed objects, and AGN. The source numbers of these types are 418, 111, and 31, respectively, from 4 MIRI pointings covering $\sim 9$ arcmin$^2$. The sample spans a redshift range of $\approx 0$--5. We derive the median SEDs for all three source types, respectively, and publicly release them. The median MIRI SED of AGN is similar to the typical SEDs of hot dust-obscured galaxies and Seyfert 2s, for which the mid-IR SEDs are dominantly from AGN-heated hot dust. Based on our SED-fit results, we estimate the black-hole accretion density (BHAD; i.e., total BH growth rate per comoving volume) as a function of redshift. At $z<3$, the resulting BHAD agrees with the X-ray measurements in general. At $z>3$, we identify a total of 28 AGN and SF-AGN mixed objects, leading to that our high-$z$ BHAD is substantially higher than the X-ray results ($\sim 1$ dex at $z \approx 4$--5). This difference indicates MIRI can identify a large population of heavily obscured AGN missed by X-ray surveys at high redshifts.

Low-altitude twisted magnetic fields may be relevant to atmospheric heating in the quiet Sun, but the exact role, topology, and formation of these twisted fields remains to be studied. We investigate the formation and evolution of a preflare flux rope in a stratified, 3D MHD simulation. One puzzle is that this modelled flux rope does not form by the usual mechanisms at work in larger flares such as flux emergence, flux cancellation, or tether-cutting. Using Lagrangian markers to trace representative field lines, we follow the spatiotemporal evolution of the flux rope. We isolate flux bundles associated with reconnecting field line pairs by focusing on thin current sheets within the flux system. We also analyze the time-varying distribution of the force-free parameter as the rope relaxes. Lastly, we compare different seeding methods for magnetic fields and discuss their relevance. We show that the modeled flux rope is gradually built from coalescing, current-carrying flux tubes. This occurs through a series of component reconnections that are driven by flows in the underlying convection zone. These reconnections lead to an inverse cascade of helicity from small to larger scales. We also find that the system attempts to relax toward a linear force-free field, but that the convective drivers and eventual nanoflare prevent full relaxation. Using a self-consistently driven simulation of a nanoflare event, we show for the first time an inverse helicity cascade tending toward a Taylor relaxation in the Sun's corona, resulting in a well-ordered flux rope that later reconnects with surrounding fields. This provides clues toward understanding the buildup of nanoflare events in the quiet Sun through incomplete Taylor relaxations when no flux emergence or cancellation is observed.

We have monitored the Didymos-Dimorphos binary asteroid in spectropolarimetric mode in the optical range before and after the DART impact. The ultimate goal was to obtain constraints on the characteristics of the ejected dust for modelling purposes. Before impact, Didymos exhibited a linear polarization rapidly increasing with phase angle, reaching a level of about 5% in the blue and about 4.5 in the red. The shape of the polarization spectrum was anti-correlated with that of its reflectance spectrum, which appeared typical of an S-class asteroid. After impact, the level of polarization dropped by about 1 percentage point (pp) in the blue band and about 0.5 pp in the red band, then continued to linearly increase with phase angle, with a slope similar to that measured prior to impact. The polarization spectra, once normalised by their values at an arbitrary wavelength, show very little or no change over the course of all observations, before and after impact. The lack of any remarkable change in the shape of the polarization spectrum after impact suggests that the way in which polarization varies with wavelength depends on the composition of the scattering material, rather than on its structure, be this a surface or a debris cloud.

Axion dark matter can be converted into photons in the magnetospheres of neutron stars leading to a spectral line centred on the Compton wavelength of the axion. Due to the rotation of the star and the plasma effects in the magnetosphere the signal is predicted to be periodic with significant time variation - a unique smoking gun for axion dark matter. As a proof of principle and to develop the methodology, we carry out the first time domain search of the signal using data from PSR J2144$-$3933 taken as part of the MeerTIME project on MeerKAT telescope. We search for specific signal templates using a matched filter technique and discuss when a time-domain analysis (as is typically the case in pulsar observations) gives greater sensitivity to the axion-coupling to photons in comparison to a simple time-averaged total flux study. We do not find any candidate signals and, hence, impose an upper limit on the axion-to-photon coupling of $g_{a\gamma\gamma}<4\times 10^{-11}\,{\rm GeV}^{-1}$ over the mass range $m_{\rm a}=3.9-4.7\,\mu{\rm eV}$ using this data. This limit relies on PSR J2144$-$3933 not being an extremely aligned rotator, as strongly supported by simple arguments based on the observed pulse profile width. We discuss the possibilities of improving this limit using future observations with MeerKAT and also SKA1-mid and the possibility of using other objects. Finally, to evade modelling uncertainties in axion radio signals, we also carry out a generic ``any periodic-signal search" in the data, finding no evidence for an axion signal.

We performed a detailed broadband spectral and timing analysis of a small flaring event of ~120 ks in a narrow-line Seyfert 1 galaxy NGC 4051 using simultaneous XMM-Newton and NuSTAR observations. The ~300 ks long NuSTAR observation and the overlapping XMM-Newton exposure were segregated into pre-flare, flare, and post-flare segments. During the flare, the NuSTAR count rate peaked at 2.5 times the mean count rate before the flare. Using various physical and phenomenological models, we examined the 0.3-50 keV X-ray spectrum, which consists of a primary continuum, reprocessed emission, warm absorber and ultra-fast outflows in different timescales. The mass of the central black hole is estimated to be >1.32*10^5 solar mass from the spectral analysis. The absence of correlation between flux in 6-7 keV and 10-50 keV bands suggests different origins of the iron emission line and the Compton hump. From the spectral analysis, we found that the reflection fraction drops significantly during the flare, accompanied by an increase in the coronal height above the disc. The spectrum became soft during the flare supporting the "softer when brighter" nature of the source. After the alleviation of the flare, the coronal height drops and the corona heats up. This indicates that there could be inflation of the corona during the flare. We found no significant change in the inner accretion disc or the seed photon temperature. These results suggest that the flaring event occurred due to the change in the coronal properties rather than any notable change in the accretion disc.

Gravitational lensing refers to the deflection of light by the gravity of celestial bodies, often predominantly composed of dark matter. Seen through a gravitational lens, the images of distant galaxies appear distorted. In this paper we discuss simulation of the image distortion by gravitational lensing. The objective is to enhance our understanding of how gravitational lensing works through a simple tool to visualise hypotheses. The simulator can also generate synthetic data for the purpose of machine learning, which will hopefully allow us to invert the distortion function, something which is not analytically possible at present.

We present the discovery of a transiting massive giant planet around TOI-4603, a sub-giant F-type star from NASA's Transiting Exoplanet Survey Satellite (TESS). The newly discovered planet has a radius of $1.042^{+0.038}_{-0.035}$ $R_{J}$, and an orbital period of $7.24599^{+0.00022}_{-0.00021}$ days. Using radial velocity measurements with the PARAS {and TRES} spectrographs, we determined the planet's mass to be $12.89^{+0.58}_{-0.57}$ $M_{J}$, resulting in a bulk density of $14.1^{+1.7}_{-1.6}$ g ${cm^{-3}}$. This makes it one of the few massive giant planets with extreme density and lies in the transition mass region of massive giant planets and low-mass brown dwarfs, an important addition to the population of less than five objects in this mass range. The eccentricity of $0.325\pm0.020$ and an orbital separation of $0.0888\pm0.0010$ AU from its host star suggest that the planet is likely undergoing high eccentricity tidal (HET) migration. We find a fraction of heavy elements of $0.13^{+0.05}_{-0.06}$ and metal enrichment of the planet ($Z_{P}/Z_{star}$) of $4.2^{+1.6}_{-2.0}$. Detection of such systems will offer us to gain valuable insights into the governing mechanisms of massive planets and improve our understanding of their dominant formation and migration mechanisms.

Super-orbital signals and negative superhumps are thought to be related to the reverse precession of the nodal line in a tilted disk, but the evidence is lacking. Our results provide new evidence for the precession of the tilted disk. Based on the TESS and K2 photometry, we investigate the super-orbital signals, negative superhumps, positive superhumps, and eclipse characteristics of the long-period eclipsing cataclysmic variable star SDSS J0812. We find super-orbital signals, negative superhumps, and positive superhumps with periods of 3.0451(5) d, 0.152047(2) d, and 0.174686(7) d, respectively, in the K2 photometry, but all disappear in the TESS photometry, where the positive superhumps are present only in the first half of the same campaign, confirming that none of them is permanently present in SDSS J0812. In addition, we find for the first time a cyclic variation of the O-C of minima, eclipse depth, and negative superhumps amplitudes for 3.045(8) d, 3.040(6) d, and 3.053(8) d in SDSS J0812, respectively, and all reach the maximum at ~ 0.75 precession phases of the tilted disk, which provides new evidence for the precession of the tilted disk. We suggest that the O-C and eclipse depth variations may come from a shift of the brightness center of the precession tilted disk. Our first finding on the periodic variation of negative superhumps amplitude with the super-orbital signals is significant evidence that the origin of negative superhumps is related to the precession of the tilted disk.

Gamma-ray observations have recently shifted the focus to higher and higher energies, with capable ground-based instruments enabling measurements in the TeV to PeV domain. While a clear prevalence of diffuse emission is observed in the GeV sky, energy-dependent cosmic-ray transport suggests a reversal of this hierarchy at higher energies. Measurements, however, are at strife regarding this question. While imaging atmospheric Cherenkov telescopes (IACTs) see a source-dominated Galactic plane, air-shower particle detectors (ASPDs) report a dominance of diffuse emission. Reconciling these claims might require a closer look at the involved instrument limitations: IACTs have a small field of view, resulting in poorer performance for large-scale emission due to the applied background subtraction technique. ASPDs have reduced resolution capabilities, resulting in unresolved sources contributing to the measurable diffuse emission signal. Here we contribute to this controversy by investigating the amount of unresolved sources in current TeV measurements in a population synthesis approach and discuss the unique capabilities for high-resolution diffuse-emission measurements with IACTs and their possibilities for overcoming their background limitations.

Primordial magnetic fields (PMF) can enhance baryon perturbations on scales below the photon mean free path. However, a magnetically driven baryon fluid becomes turbulent near recombination, thereby damping out baryon perturbations below the turbulence scale. In this letter, we show that the growth of baryon perturbations is imprinted in the dark matter perturbations, which are unaffected by turbulence and eventually collapse to form $10^{-11}-10^3\ M_{\odot}$ dark matter minihalos. In the process, we analytically derive the evolution of the PMF power spectrum in the viscous drag regime. If the magnetic fields purportedly detected in the blazar observations are PMFs generated after inflation and have a Batchelor spectrum, then such PMFs should also produce minihalos.

As miniaturized spacecraft (e.g., cubesats and smallsats) and instrumentation become an increasingly indispensable part of space exploration and scientific investigations, it is important to understand their potential susceptibility to space weather impacts resulting from the sometimes volatile space environment. There are multitude of complexities involved in how space environment interacts with different space hardware/electronics. Measurements of such impacts, however, have been lacking. Therefore, we recommend developing and/or procuring low-cost, low-power consumption, and compact sensor suites (mainly for space weather and impact purposes) and flying them on all future NASA (and U.S in general) missions in order to measure and quantify space weather impacts, in addition to the main instrumentation.

Recently multiscale imaging approaches such as DoG-HiT were developed to solve the VLBI imaging problem and showed a promising performance: they are fast, accurate, unbiased and automatic. We extend the multiscalar imaging approach to polarimetric imaging, reconstructions of dynamically evolving sources and finally to dynamic polarimetric reconstructions. These extensions (mr-support imaging) utilize a multiscalar approach. The time-averaged Stokes I image is decomposed by a wavelet transform into single subbands. We use the set of statistically significant wavelet coefficients, the multiresolution support, computed by DoG-HiT as a prior in a constrained minimization manner: we fit the single-frame (polarimetric) observables by only varying the coefficients in the multiresolution support. The EHT is a VLBI array imaging supermassive black holes. We demonstrate on synthetic data that mr-support imaging offers ample regularization and is able to recover simple geometric dynamics at the horizon scale in a typical EHT setup. The approach is relatively lightweight, fast and largely automatic and data driven. The ngEHT is a planned extension of the EHT designed to recover movies at the event horizon scales of a supermassive black hole. We benchmark the performance of mr-support imaging for the denser ngEHT configuration demonstrating the major improvements the additional ngEHT antennas will bring to dynamic, polarimetric reconstructions. Current and upcoming instruments offer the observational possibility to do polarimetric imaging of dynamically evolving structural patterns with highest spatial and temporal resolution. State-of-the-art dynamic reconstruction methods can capture this motion with a range of temporal regularizers and priors. With this work, we add an additional, simpler regularizer to the list: constraining the reconstruction to the multiresolution support.

$\require{mediawiki-texvc}$ We study accretion disk emission from eight Seyfert $1 - 1.5$ active galactic nuclei (AGN) using far ultra-violet ($1300-1800$ ${\AA}$) slit-less grating spectra acquired with AstroSat/UVIT. We correct for the Galactic and intrinsic extinction, contamination from the host galaxies, narrow and broad-line regions, Fe II emission and Balmer continuum, and derive the intrinsic continua. We use HST COS/FOS spectra to account for the emission/absorption lines in the low-resolution UVIT spectra. We find generally redder power-law ($f_\nu \propto \nu^{\alpha}$) slopes ($\alpha \sim -1.1 - 0.3$) in the far UV band than predicted by the standard accretion disk model in the optical/UV band. We fit accretion disk models such as the multi-temperature disk blackbody ($\texttt{DISKBB}$) and relativistic disk ($\texttt{ZKERRBB}$, $\texttt{OPTXAGNF}$) models to the observed intrinsic continuum emission. We measure the inner disk temperatures using the $\texttt{DISKBB}$ model for seven AGN. These temperatures in the range $\sim 3.6-5.8$ eV are lower than the peak temperatures predicted for standard disks around maximally spinning super-massive black holes accreting at Eddington rates. The inner disks in two AGN, NGC 7469 and Mrk 352, appear to be truncated at $\sim 35-125r_{g}$ and $50-135r_{g}$, respectively. While our results show that the intrinsic FUV emission from the AGN are consistent with the standard disks, it is possible that UV continua may be affected by the presence of soft X-ray excess emission, X-ray reprocessing, and thermal Comptonisation in the hot corona. Joint spectral modeling of simultaneously acquired UV/X-ray data may be necessary to further investigate the nature of accretion disks in AGN.

We present a first statistical sample of faint type-1 AGNs at $z>4$ identified by JWST/NIRSpec deep spectroscopy. Among the 185 galaxies at $z_\mathrm{spec}=3.8-8.9$ confirmed with NIRSpec, our systematic search for broad-line emission reveals 10 type-1 AGNs at $z=4.015-6.936$ whose broad component is only seen in the permitted H$\alpha$ line and not in the forbidden [OIII]$\lambda$5007 line that is detected with greater significance than H$\alpha$. The broad H$\alpha$ line widths of $\mathrm{FWHM}\simeq1000-6000\ \mathrm{km\ s^{-1}}$ suggest that the AGNs have low-mass black holes with $M_\mathrm{BH}\sim10^6-10^7\ M_\odot$, remarkably lower than those of low-luminosity quasars previously identified at $z>4$ with ground-based telescopes. JWST and HST high-resolution images reveal that the majority of them show extended morphologies indicating significant contribution to the total lights from their host galaxies, except for three compact objects two of which show red SEDs, probably in a transition phase from faint AGNs to low luminosity quasars. Careful AGN-host decomposition analyses show that their host's stellar masses are systematically lower than the local relation between the black hole mass and the stellar mass, implying a fast black hole growth consistent with predictions from theoretical simulations. A high fraction of the broad-line AGNs ($\sim5\%$), higher than $z\sim0$, indicates that a number density of such faint AGNs is higher than an extrapolation of the quasar luminosity function, implying a large population of AGNs including type 1 and type 2 in the early universe. Such faint AGNs contribute to cosmic reionization, while the total contribution is not large, up to $\sim50\%$ at $z\sim6$, because of their faint nature.

Super-Mercuries, rocky exoplanets with bulk iron mass fraction of more than 60 per cent, appear to be preferentially hosted by stars with higher iron mass fraction than the Earth. It is unclear whether these iron-rich planets can form in the disc, or if giant impacts are necessary. Here we investigate the formation of super-Mercuries in their natal protoplanetary discs by taking into account their host stars' abundances (Fe, Mg, Si, S). We employ a disc evolution model which includes the growth, drift, evaporation and recondensation of pebbles to compute the pebble iron mass fraction. The recondensation of outward-drifting iron vapour near the iron evaporation front is the key mechanism that facilitates an increase in the pebble iron mass fraction. We also simulate the growth of planetary seeds around the iron evaporation front using a planet formation model which includes pebble accretion and planet migration, and compute the final composition of the planets. Our simulations are able to reproduce the observed iron compositions of the super-Mercuries provided that all the iron in the disc are locked in pure Fe grains and that the disc viscosity is low. The combined effects of slow orbital migration of planets and long retention time of iron vapour in low-viscosity discs makes it easier to form iron-rich planets. Furthermore, we find that decreasing the stellar Mg/Si ratio results in an increase in the iron mass fraction of the planet due to a reduction in the abundance of Mg2SiO4, which has a very similar condensation temperature as iron, in the disc. Our results thus imply that super-Mercuries are more likely to form around stars with low Mg/Si, in agreement with observational data.

This work discusses a few theories including the interaction of dark matter, cosmic strings, and locally coupled dark energy. The paper also examines mathematical models used to describe the pressure and density within a star, including the polytropic relationship and the Lane-Emden equation. Simulation results from the IllustrisTNG datasets are also presented, providing insights into the interacting dark matter solutions. With the derived solutions this paper, it explores the possible causes for the sudden disappearance of the star PHL293B-LBV

We present ALMA Band 9 continuum observation of the ultraluminous quasi-stellar object (QSO) SDSS J0100+2802, providing a $\sim 10\sigma$ detection at $\sim 670$ GHz. SDSS J0100+2802 is the brightest QSO with the most massive super massive black hole (SMBH) known at $z>6$, and we study its dust spectral energy distribution in order to determine the dust properties and the star formation rate (SFR) of its host-galaxy. We obtain the most accurate estimate so far of the temperature, mass and emissivity index of the dust, having $T_{\rm dust}=48.4\pm2.3$ K, $M_{\rm dust}=(2.29\pm0.83)\times 10^7$ M$_\odot$, $\beta=2.63\pm 0.23$. This allows us to measure the SFR with the smallest statistical error for this QSO, SFR$=265\pm 32\ \rm M_\odot yr^{-1}$. Our results enable us to evaluate the relative growth of the SMBH and host galaxy of J0100+2802, finding that the SMBH is dominating the process of BH-galaxy growth in this QSO at $z=6.327$, when the Universe was $865$ Myr old. Such unprecedented constraints on the host galaxy SFR and dust temperature can only be obtained through high frequency observations, and highlight the importance of ALMA Band 9 to obtain a robust overview of the build-up of the first quasars' host galaxies at $z>6$.

The repeating FRB 20121102A was localized to a star-forming region in a dwarf galaxy and found to be co-located with a persistent radio source (PRS). FRB 20190520B is only the second known source sharing phenomenology akin to FRB 20121102A's, with similar burst activity, host galaxy properties, as well as being associated with a PRS. PRS emission is potentially a calorimeter, allowing us to estimate the energy output of the central FRB engine. Independently of FRB studies, PRSs have been found in dwarf galaxies and interpreted as intermediate mass black holes. To improve our understanding of such sources associated with dwarf galaxies, it is imperative to increase the known sample size. Here, we present a search for compact radio sources coincident with dwarf galaxies, discuss source candidates and planned strategies to differentiate them between candidate FRB hosts and intermediate mass black holes.

Using Physics-Informed Neural Networks (PINNs) to solve a specific boundary value problem is becoming more popular as an alternative to traditional methods. However, depending on the specific problem, they could be computationally expensive and potentially less accurate. The functionality of PINNs for real-world physical problems can significantly improve if they become more flexible and adaptable. To address this, our work explores the idea of training a PINN for general boundary conditions and source terms expressed through a limited number of coefficients, introduced as additional inputs in the network. Although this process increases the dimensionality and is computationally costly, using the trained network to evaluate new general solutions is much faster. Our results indicate that PINN solutions are relatively accurate, reliable, and well-behaved. We applied this idea to the astrophysical scenario of the magnetic field evolution in the interior of a neutron star connected to a force-free magnetosphere. Solving this problem through a global simulation in the entire domain is expensive due to the elliptic solver's needs for the exterior solution. The computational cost with a PINN was more than an order of magnitude lower than the similar case solved with classical methods. These results pave the way for the future extension to 3D of this (or a similar) problem, where generalised boundary conditions are very costly to implement.

We present an analysis of two seasons of archival, multi-frequency VLA monitoring of the quad lens system JVAS B1422+231, the 15-GHz data of which have previously been published. The 8.4- and 15-GHz variability curves show significant variability, especially in polarization, but lack features on short time-scales that would be necessary for an accurate measurement of the very short predicted time delays ($\le$1 d) between the three bright images. Time delays can only realistically be measured to the very faint image D and for the first time we detect its long-term variability and determine its polarization properties. However, image-dependent (extrinsic) variability (including variations on time-scales of hours) is present in multiple images and the magnitude of this is largest in image D at 15 GHz ($\pm$10 per cent). As the variations appear to increase in amplitude with frequency, we suggest that the most likely cause is microlensing by compact objects in the lensing galaxy. Combining the monitoring data allows us to detect a faint arc of emission lying between images B and C and the jets responsible for this are imaged using archival VLBA data. Finally, we have also measured the rotation measure of the three bright images and detected the polarization properties of image D.

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

We carried out an observational search for the recently discovered molecule H2NC, and its more stable isomer H2CN, toward eight cold dense clouds (L1544, L134N, TMC-2, Lupus-1A, L1489, TMC-1 NH3, L1498, and L1641N) and two diffuse clouds (B0415+379 and B0355+508) in an attempt to constrain its abundance in different types of interstellar regions and shed light on its formation mechanism. We detected H2NC in most of the cold dense clouds targeted, 7 out of 8, while H2CN was only detected in 5 out of 8 clouds. The column densities derived for both H2NC and H2CN are in the range 1e11-1e12 cm-2 and the abundance ratio H2NC/H2CN varies between 0.51 and >2.7. The metastable isomer H2NC is therefore widespread in cold dense clouds where it is present with an abundance similar to that of H2CN. We did not detect either H2NC or H2CN in any of the two diffuse clouds targeted, which does not allow to shed light on how the chemistry of H2NC and H2CN varies between dense and diffuse clouds. We found that the column density of H2NC is correlated with that of NH3, which strongly suggests that these two molecules are chemically linked, most likely ammonia being a precursor of H2NC through the C + NH3 reaction. We performed electronic structure and statistical calculations which show that both H2CN and H2NC can be formed in the C + NH3 reaction through two different channels involving two different transition states which lie very close in energy. The predicted product branching ratio H2NC/H2CN is very method dependent but values between 0.5 and 0.8 are the most likely ones. Therefore, both the astronomical observations and the theoretical calculations support that the reaction C + NH3 is the main source of H2NC in interstellar clouds.

Pulsar Wind Nebulae (PWNe) shine at multi-wavelengths and are expected to constitute the largest class of gamma-ray sources in our Galaxy. They are known to be very efficient particle accelerators: the Crab nebula, the PWNe class prototype, is the unique firmly identified leptonic PeVatron of the Galaxy to date, and most of the PeVatrons recently detected by LHAASO appear to be compatible with a pulsar origin. PWNe have been proved to be associated with the formation of misaligned X-ray tails and TeV halos, as sign of an efficient escape of energetic particles from the PWN into the surrounding medium. With the advent of the Cherenkov Telescope Array we expect that ~200 new PWNe will be detected. Being able to correctly model their multi-wavelength spectral properties, spatial and spectral morphology at gamma-rays is then topical today. This in particular means we should be able to account for their different evolutionary phases, and to correctly determine the influence they have on the spectral properties of the source. This indeed reflects directly on the expectation of how many PWNe will be detected at gamma-rays. Finally, the identification of PWNe in future gamma-ray data, not only is relevant for their scientific importance, but also to allow for the identification of less prominent sources that might be hidden by the background of non-identified PWNe.

Measuring the size distribution of small (km-scale) KBOs can help constrain models of Solar System formation and planetary migration. Such small, distant bodies are hard to detect with current or planned telescopes, but can be identified as sub-second occultations of background stars. We present the analysis of data from the Weizmann Fast Astronomical Survey Telescope (W-FAST), consisting of fast photometry of ~10^6 star hours at a frame rate of 10-25 Hz. Our pipeline utilizes a matched-filter approach with a large template bank, including red-noise treatment, and injection of simulated events for estimating the detection efficiency. The KBO radius at which our survey is 10% (50%) efficient is 1.1 (2.0) km. The data from 2020-2021 observing seasons were analyzed and no occultations were identified. We discuss a sample of sub-second false-positive events, both occultation-like and flare-like, which are still not fully understood but could be instructive for future surveys looking for short-duration events. We use our null-detection result to set limits on the km-scale KBO number density. Our individual radius bin limits are consistent with most previous works, with N(r>1km) <=10^6 deg^-2 (95% confidence limit). Our integrated (all size) limits, assuming a power law normalized to large (~45 km) KBOs gives a power law index q<3.93 (95% confidence limit). Finally, our results are in tension with a recently reported KBO detection from the ground, at the p=4x10^-4 level.

We revisit the one-loop correction in curvature perturbation power spectrum in models of single field inflation which undergo a phase of ultra slow-roll (USR) inflation. We include the contributions from both the cubic and quartic interaction Hamiltonians and calculate the one-loop corrections on the spectrum of the CMB scale modes from the small scale modes which leave the horizon during the USR phase. It is shown that the amplitude of one-loop corrections depends on the sharpness of the transition from the USR phase to the final slow-roll phase. For an arbitrarily sharp transition, the one-loop correction becomes arbitrarily large, invalidating the perturbative treatment of the analysis. We speculate that for a mild transition, the large one-loop corrections are washed out during the subsequent evolution after the USR phase. The implications for primordial black holes formation are briefly reviewed.

Over the past decade, hundreds of nights have been spent on the worlds largest telescopes to search for and directly detect new exoplanets using high-contrast imaging (HCI). Thereby, two scientific goals are of central interest: First, to study the characteristics of the underlying planet population and distinguish between different planet formation and evolution theories. Second, to find and characterize planets in our immediate Solar neighborhood. Both goals heavily rely on the metric used to quantify planet detections and non-detections. Current standards often rely on several explicit or implicit assumptions about the noise. For example, it is often assumed that the residual noise after data post-processing is Gaussian. While being an inseparable part of the metric, these assumptions are rarely verified. This is problematic as any violation of these assumptions can lead to systematic biases. This makes it hard, if not impossible, to compare results across datasets or instruments with different noise characteristics. We revisit the fundamental question of how to quantify detection limits in HCI. We focus our analysis on the error budget resulting from violated assumptions. To this end, we propose a new metric based on bootstrapping that generalizes current standards to non-Gaussian noise. We apply our method to archival HCI data from the NACO-VLT instrument and derive detection limits for different types of noise. Our analysis shows that current standards tend to give detection limit that are about one magnitude too optimistic in the speckle-dominated regime. That is, HCI surveys may have excluded planets that can still exist.

Radio activity in AGN (Active Galactic Nuclei) produce feedback on the host galaxy via the impact of the relativistic jets on the circumnuclear gas. Although radio jets can reach up to several times the optical radius of the host galaxy, in this review we focus on the observation of the feedback deposited locally in the central region of the host galaxies, in the form of outflows due to the jet-gas interaction. We begin by discussing how galaxy mergers and interactions are the most favored scenario for triggering radio AGN after gas accretion to the nuclear supermassive black hole and star formation enhancement in the nuclear region, observed in particular in the most luminous sources. We then discuss observational signatures of the process of jet-gas coupling, in particular the resulting outflows and their effects on the host galaxy. These include the presence of shock signatures and the detection of outflows not only along the radio jet but perpendicular to it in many sources. Although most of the studies are done via the observation of ionized gas, molecular gas is also being increasingly observed in outflow, contributing to the bulk of the mass outflow rate. Even though most radio sources present outflow kinetic powers that do not reach $1\%\,L_{bol}$, and thus do not seem to provide an immediate impact on the host galaxy, they act to heat the ISM gas, preventing star formation, slowing the galaxy mass build-up process and limiting the stellar mass growth, in a ``maintenance mode" feedback.

We report on an X-ray polarimetric observation of the high-mass X-ray binary LMC X-1 in the high/soft state, obtained by the Imaging X-ray Polarimetry Explorer (IXPE) in October 2022. The measured polarization is below the minimum detectable polarization of 1.1 per cent (at the 99 per cent confidence level). Simultaneously, the source was observed with the NICER, NuSTAR and SRG/ART-XC instruments, which enabled spectral decomposition into a dominant thermal component and a Comptonized one. The low 2-8 keV polarization of the source did not allow for strong constraints on the black-hole spin and inclination of the accretion disc. However, if the orbital inclination of about 36 degrees is assumed, then the upper limit is consistent with predictions for pure thermal emission from geometrically thin and optically thick discs. Assuming the polarization degree of the Comptonization component to be 0, 4, or 10 per cent, and oriented perpendicular to the polarization of the disc emission (in turn assumed to be perpendicular to the large scale ionization cone detected in the optical and radio bands), an upper limit to the polarization of the disc emission of 0.5, 1.7, or 3.6 per cent, respectively, is found (at the 99 per cent confidence level).

Despite the appearance of two- and three-dimensional models thanks to the rapid growth of computing performance, numerical hydrocodes used to model radial stellar pulsations still apply a one-dimensional stellar envelope model without any realistic atmosphere, in which a significant improvement was the inclusion of turbulent convection. However, turbulent convection is an inherently multi-dimensional physical process in the vicinity of the ionization zones that generate pulsation. The description of these processes in one dimension can only be approximated based on simplified theoretical considerations involving several undetermined dimensionless parameters. In this work, we confront two one-dimensional numerical codes, namely the Budapest-Florida code (BpF) and the MESA Radial Stellar Pulsations module (RSP), with radial velocity observations of several non-modulated RRab stars of the M3 globular cluster and specified the undetermined convective parameters by the measured data for both codes independently. Our determination shows that some of the parameters depend on the effective temperature, which dependence is established for the first time in this work, and we also found some degeneracy between the parameters. This procedure gives as by-product suggestions for parameters of the publicly available RSP code extensively used recently by researchers through the MESA package. This work is part of the preparatory work to establish a theoretical framework required to make progress based on the results of one-dimensional models to couple them with multi-dimensional ones for further detailed analysis of physical processes.

Semi-numerical simulations are the leading candidates for evolving reionization on cosmological scales. These semi-numerical models are efficient in generating large-scale maps of the 21cm signal, but they are too slow to enable inference at the field level. We present different strategies to train a U-Net to accelerate these simulations. We derive the ionization field directly from the initial density field without using the ionizing sources' location, and hence emulating the radiative transfer process. We find that the U-Net achieves higher accuracy in reconstructing the ionization field if the input includes either white noise or a noisy version of the ionization map beside the density field during training. Our model reconstructs the power spectrum over all scales perfectly well. This work represents a step towards generating large-scale ionization maps with a minimal cost and hence enabling rapid parameter inference at the field level.

A novel method for investigating the sensitivity of the tip of the red giant branch (TRGB) I band magnitude $M_I$ to stellar input physics is presented. We compute a grid of $\sim$125,000 theoretical stellar models with varying mass, initial helium abundance, and initial metallicity, and train a machine learning emulator to predict $M_I$ as a function of these parameters. First, our emulator can be used to theoretically predict $M_I$ in a given galaxy using Monte Carlo sampling. As an example, we predict $M_I = -3.84^{+0.14}_{-0.12}$ in the Large Magellanic Cloud. Second, our emulator enables a direct comparison of theoretical predictions for $M_I$ with empirical calibrations to constrain stellar modeling parameters using Bayesian Markov Chain Monte Carlo methods. We demonstrate this by using empirical TRGB calibrations to obtain new independent measurements of the metallicity in three galaxies. We find $Z=0.0117^{+0.0083}_{-0.0055}$ in the Large Magellanic Cloud, $Z=0.0077^{+0.0074}_{-0.0038}$ in NGC 4258, and $Z=0.0111^{+0.0083}_{-00.0056}$ in $\omega$-Centauri, consistent with other measurements. Other potential applications of our methodology are discussed.

Massive vector particles are minimal dark matter candidates that motivate a wide range of laboratory searches, primarily exploiting a postulated kinetic mixing with the photon. However, depending on the high energy field content, the dominant vector dark matter (VDM) coupling to visible particles may arise at higher operator dimension, motivating efforts to predict direct detection rates for more general interactions. Here we present the first calculation of VDM absorption through its coupling to electron electric (EDM) or magnetic (MDM) dipole moments, which can be realized in minimal extensions to the Standard Model and yield the observed abundance through a variety of mechanisms across the eV\,-\,MeV mass range. We compute the absorption rate of the MDM and EDM models for a general target, and then derive direct detection constraints from targets currently in use: Si and Ge crystals and Xe and Ar atoms. We find that current experiments are already sensitive to VDM parameter space corresponding to a cosmological freeze-in scenario, and future experiments will be able to completely exclude MDM and EDM freeze-in models with reheat temperatures below the electroweak scale. Additionally, we find that while constraints on the MDM interaction can be related to constraints on axion-like particles, the same is not true for the EDM model, so the latter absorption rate must be computed from first principles. To achieve this, we update the publicly available program EXCEED-DM to perform these new calculations.

Next-generation large-volume detectors, such as GRAND, POEMMA, Trinity, TAROGE-M, and PUEO, have been designed to search for ultra-high-energy cosmic rays (UHECRs) with unprecedented sensitivity. We propose to use these detectors to search for new physics beyond the Standard Model (BSM). By considering the simple case of a right-handed neutrino that mixes exclusively with the active $\tau$ neutrino, we demonstrate that the existence of new physics can increase the probability for UHECRs to propagate through the Earth and produce extensive air showers that will be measurable soon. We compare the fluxes of such showers that would arise from various diffuse and transient sources of high-energy neutrinos, both in the Standard Model and in the presence of a right-handed neutrino. We show that detecting events with emergence angles $\gtrsim 10$ deg is promising to probe the existence of BSM physics, and we study the sensitivity of GRAND and POEMMA to do so. In particular, we show that the hypothesis of a right-handed neutrino with a mass of $\mathcal O(1-16)$ GeV may be probed in the future for mixing angles as small as $|U_{\tau N}|^2 \gtrsim 10^{-7}$, thus competing with existing and projected experimental limits.

Bimetric gravity, is a theory of gravity that posits the existence of two interacting and dynamical metric tensors. The spectrum of bimetric gravity consists of a massless and a massive spin-2 particle. The form of the interactions between the two metrics $g_{\mu\nu}$ and $f_{\mu\nu}$ is constrained by requiring absence of the so called Boulware-Deser ghost. In this work we extend the original bimetric theory to its bimetric-affine counterpart, in which the associated two connections, $\Gamma_{ \mu\,\,\,\nu}^{\,\,\,\rho}(g)$ and $\widetilde{\Gamma}_{ \mu\,\,\,\nu}^{\,\,\,\rho}(f)$, are treated as independent variables. We examine in detail the case of an additional quadratic in the Ricci scalar curvature term $\mathcal{R}^2(g)$ and we find that this theory is free of ghosts for a wide range of the interaction parameters, not excluding the possibility of a Dark Matter interpretation of the massive spin-2 particle.

Well-motivated scenarios of thermally-produced dark matter often result in a population of electrons and positrons within galaxies produced through dark matter annihilation -- often in association with high-energy gamma rays. As they diffuse through galactic magnetic fields, these $e^\pm$ produce synchrotron radio emission. The intensity and morphology of this signal depends on the properties of the interstellar medium through which the $e^\pm$ propagate. Using observations of the Andromeda Galaxy (M31) to construct a model of the gas, magnetic fields, and starlight, we set constraints on dark matter annihilation to $b\bar{b}$ using the morphology of 3.6 cm radio emission. As the emission signal at the center of M31 is very sensitive to the diffusion coefficient and dark matter profile, we base our limits on the differential flux in the region between $0.9-6.9$ kpc from the center. We exclude annihilation cross sections $\gtrsim 3 \times 10^{-25}$ cm$^3$/s in the mass range $10-500$ GeV, with a maximum sensitivity of $7\times 10^{-26}$ cm$^3$/s at $20-40$ GeV. Though these limits are weaker than those found in previous studies of M31, they are robust to variations of the diffusion coefficient.

In dark-matter annihilation channels to hadronic final states, stable particles -- such as positrons, photons, antiprotons, and antineutrinos -- are produced via complex sequences of phenomena including QED/QCD radiation, hadronisation, and hadron decays. These processes are normally modelled by Monte Carlo event generators whose limited accuracy imply intrinsic QCD uncertainties on the predictions for indirect-detection experiments like Fermi-LAT, Pamela, IceCube or AMS-02. In this article, we perform a complete analysis of QCD uncertainties in antimatter spectra from dark-matter annihilation, based on parametric variations of the Pythia 8 event generator. After performing several retunings of light-quark fragmentation functions, we define a set of variations that span a conservative estimate of the QCD uncertainties. We estimate the effects on antimatter spectra for various annihilation channels and final-state particle species, and discuss their impact on fitted values for the dark-matter mass and thermally-averaged annihilation cross section. We find dramatic impacts which can go up to $\mathcal{O}(40)$ GeV for uncertainties on the dark-matter mass and up to $\mathcal{O}(10\%)$ for the annihilation cross section. We provide the spectra in tabulated form including QCD uncertainties and code snippets to perform fast dark-matter fits, in this https://github.com/ajueid/qcd-dm.github.io.git repository.

New feebly interacting particles would emerge from a supernova core with 100-MeV-range energies and produce $\gamma$-rays by subsequent decays. These would contribute to the diffuse cosmic $\gamma$-ray background or would have shown up in the Solar Maximum Mission (SMM) satellite from SN~1987A. However, we show for the example of axion-like particles (ALPs) that, even at distances beyond the progenitor star, the decay photons may not escape, and can instead form a fireball, a plasma shell with $T\lesssim 1$ MeV. Thus, existing arguments do not exclude ALPs with few 10 MeV masses and a two-photon coupling of a few $10^{-10}~{\rm GeV}^{-1}$. However, the energy would have showed up in sub-MeV photons, which were not seen from SN 1987A in the Pioneer Venus Orbiter (PVO), closing again this new window. A careful re-assessment is required for other particles that were constrained in similar ways.

These lecture notes collect the material that I have been using over the years for various short courses on the physics of gravitational waves, first at the Institut d'Astrophysique de Paris (France), and then at SISSA (Italy) and various summer/winter schools. The level should be appropriate for PhD students in physics or for MSc students that have taken a first course in general relativity. The focus is on deriving results from first principles, rather than on astrophysical applications.

We derive the expressions on the observed light-cone for some relevant cosmological gauge invariant variables, such as the Mukhanov-Sasaki variable and $E$- and $B$- modes of the tensor perturbations. Since the structure of the light-cone does not reflect in a direct way the FLRW symmetries, we develop a formalism which is coordinate independent and classifies the perturbations according to their helicities. Even though we work with linear perturbations, our formalism can be readily extended to non-linear theory and put the basis to study the evolution of cosmological perturbations, since the early- until the the late-time Universe, directly along the observed light-cone.

Cosmic gravitons are expected in the MHz-GHz regions that are currently unreachable by the operating wide-band interferometers and where various classes of electromechanical detectors have been proposed through the years. The minimal chirp amplitude detectable by these instruments is often set on the basis of the sensitivities reachable by the detectors currently operating in the audio band. By combining the observations of the pulsar timing arrays, the limits from wide-band detectors and the other phenomenological bounds we show that this requirement is far too generous and even misleading since the actual detection of relic gravitons well above the kHz would demand chirp and spectral amplitudes that are ten or even fifteen orders of magnitude smaller than the ones currently achievable in the audio band, for the same classes of stochastic sources. We then examine more closely the potential high-frequency signals and show that the sensitivity in the chirp and spectral amplitudes must be even smaller than the ones suggested by the direct and indirect constraints on the cosmic gravitons. We finally analyze the high-frequency detectors in the framework of Hanbury-Brown Twiss interferometry and argue that they are actually more essential than the ones operating in the audio band (i.e. between few Hz and few kHz) if we want to investigate the quantumness of the relic gravitons and their associated second-order correlation effects. We suggest, in particular, how the statistical properties of thermal and non-thermal gravitons can be distinguished by studying the corresponding second-order interference effects.