Papers for Wednesday, Apr 26 2023

An intriguing reference to the existence of a self-portrait by Galileo Galilei is contained in the biography of the scientist by Thomas Salusbury dated ca. 1665, of which only one incomplete and inaccessible copy exists. Galileo grew up in a Renaissance atmosphere, acquiring an artistic touch. He was a musician, a writer and also a painter, as reported by Viviani and documented by his watercolours of the Moon and drawings of solar spots. Recently a new portrait with a remarkable similarity to the portraits of Galileo Galilei by Santi di Tito (1601), Domenico Tintoretto (ca. 1604), and Furini (ca. 1612) has been found and examined using sophisticated face recognition techniques. If the identity could be confirmed, other elements, such as the young age of Galileo or the seam in the canvas revealed by infrared and X-ray analysis, may suggest a possible link with the self-portrait mentioned by Salusbury.

Over 200 quasars have now been discovered at $z >$ 6, including nine at $z >$ 7. They are thought to form from the collapse of supermassive primordial stars to 10$^4$ - 10$^5$ M$_{\odot}$ black holes at $z \sim$ 20 - 25, which then rapidly grow in the low-shear environments of rare, massive halos fed by strong accretion flows. Sensitive new radio telescopes such as the Next-Generation Very Large Array (ngVLA) and the Square Kilometer Array (SKA) could probe the growth of these objects at much earliest stages than is now possible. Here, we estimate radio flux from the first quasars at $z \sim$ 6 - 15 at 1 - 10 GHz. We find that a quasar with properties similar to that of ULAS J1120+0641, a 2.1 $\times$ 10$^9$ M$_{\odot}$ black hole at $z =$ 7.1, could be detected at up to $z \sim$ 16 by the SKA and at $z \sim$ 14 by the ngVLA. The advent of these new observatories, together with the James Webb Space Telescope (JWST), Euclid, and the Roman Space Telescope (RST), will inaugurate the era of $z \lesssim$ 15 quasar astronomy in the coming decade.

An era of kination occurs when the Universe's energy density is dominated by a fast-rolling scalar field. Dark matter that is thermally produced during an era of kination requires larger-than-canonical annihilation cross sections to generate the observed dark matter relic abundance. Furthermore, dark matter density perturbations that enter the horizon during an era of kination grow linearly with the scale factor prior to radiation domination. We show how the resulting enhancement to the small-scale matter power spectrum increases the microhalo abundance and boosts the dark matter annihilation rate. We then use gamma-ray observations to constrain thermal dark matter production during kination. The annihilation boost factor depends on the minimum halo mass, which is determined by the small-scale cutoff in the matter power spectrum. Therefore, observational limits on the dark matter annihilation rate imply a minimum cutoff scale for a given dark matter particle mass and kination scenario. For dark matter that was once in thermal equilibrium with the Standard Model, this constraint establishes a maximum allowed kinetic decoupling temperature for the dark matter. This bound on the decoupling temperature implies that the growth of perturbations during kination cannot appreciably boost the dark matter annihilation rate if dark matter was once in thermal equilibrium with the Standard Model.

The radial structure of debris discs can encode important information about their dynamical and collisional history. In this paper we present a 3-phase analytical model to analyse the collisional evolution of solids in debris discs, focusing on their joint radial and temporal dependence. Consistent with previous models, we find that as the largest planetesimals reach collisional equilibrium in the inner regions, the surface density of dust and solids becomes proportional to $\sim r^{2}$ within a certain critical radius. We present simple equations to estimate the critical radius and surface density of dust as a function of the maximum planetesimal size and initial surface density in solids (and vice versa). We apply this model to ALMA observations of 7 wide debris discs. We use both parametric and non-parametric modelling to test if their inner edges are shallow and consistent with collisional evolution. We find that 4 out of 7 have inner edges consistent with collisional evolution. Three of these would require small maximum planetesimal sizes below 10 km, with HR 8799's disc potentially lacking solids larger than a few centimeters. The remaining systems have inner edges that are much sharper, which requires maximum planetesimal sizes $\gtrsim10$ km. Their sharp inner edges suggest they could have been truncated by planets, which JWST could detect. In the context of our model, we find that the 7 discs require surface densities below a Minimum Mass Solar Nebula, avoiding the so-called disc mass problem. Finally, during the modelling of HD 107146 we discover that its wide gap is split into two narrower ones, which could be due to two low-mass planets formed within the disc.

We study weak gravitational lensing convergence maps produced from the MillenniumTNG (MTNG) simulations by direct projection of the mass distribution on the past backwards lightcone of a fiducial observer. We explore the lensing maps over a large dynamic range in simulation mass and angular resolution, allowing us to establish a clear assessment of numerical convergence. By comparing full physics hydrodynamical simulations with corresponding dark-matter-only runs we quantify the impact of baryonic physics on the most important weak lensing statistics. Likewise, we predict the impact of massive neutrinos reliably far into the non-linear regime. We also demonstrate that the "fixed & paired" variance suppression technique increases the statistical robustness of the simulation predictions on large scales not only for time slices but also for continuously output lightcone data. We find that both baryonic and neutrino effects substantially impact weak lensing shear measurements, with the latter dominating over the former on large angular scales. Thus, both effects must explicitly be included to obtain sufficiently accurate predictions for stage IV lensing surveys. Reassuringly, our results agree accurately with other simulation results where available, supporting the promise of simulation modelling for precision cosmology far into the non-linear regime.

The massive end of the gas-phase mass--metallicity relation (MZR) is a sensitive probe of active galactic nuclei (AGN) feedback that is a crucial but highly uncertain component of galaxy evolution models. In this paper, we extend the $z\sim0.7$ MZR by $\sim$0.5 dex up to log$(M_\star/\textrm{M}_\odot)\sim11.1$. We use extremely deep VLT VIMOS spectra from the Large Early Galaxy Astrophysics Census (LEGA-C) survey to measure metallicities for 145 galaxies. The LEGA-C MZR matches the normalization of the $z\sim0.8$ DEEP2 MZR where they overlap, so we combine the two to create an MZR spanning from 9.3 to 11.1 log$(M_\star/\textrm{M}_\odot)$. The LEGA-C+DEEP2 MZR at $z\sim0.7$ is offset to slightly lower metallicities (0.05-0.13 dex) than the $z\sim0$ MZR, but it otherwise mirrors the established power law rise at low/intermediate stellar masses and asymptotic flattening at high stellar masses. We compare the LEGA-C+DEEP2 MZR to the MZR from two cosmological simulations (IllustrisTNG and SIMBA), which predict qualitatively different metallicity trends for high-mass galaxies. This comparison highlights that our extended MZR provides a crucial observational constraint for galaxy evolution models in a mass regime where the MZR is very sensitive to choices about the implementation of AGN feedback.

The intrinsic alignment (IA) of observed galaxy shapes with the underlying cosmic web is a source of contamination in weak lensing surveys. Sensitive methods to identify the IA signal will therefore need to be included in the upcoming weak lensing analysis pipelines. Hydrodynamical cosmological simulations allow us to directly measure the intrinsic ellipticities of galaxies and thus provide a powerful approach to predict and understand the IA signal. Here we employ the novel, large-volume hydrodynamical simulation MTNG740, a product of the MillenniumTNG (MTNG) project, to study the IA of galaxies. We measure the projected correlation functions between the intrinsic shape/shear of galaxies and various tracers of large-scale structure, $w_{+g},\ w_{+m},\ w_{++}$ over the radial range $r_{\rm p} \in [0.02 , 200]\,h^{-1}{\rm Mpc}$ and at redshifts $z=0.0$, $0.5$ and $1.0$. We detect significant signal-to-noise IA signals with the density field for both elliptical and spiral galaxies. We also find significant intrinsic shear-shear correlations for ellipticals. We further examine correlations of the intrinsic shape of galaxies with the local tidal field. Here we find a significant IA signal for elliptical galaxies assuming a linear model. We also detect a weak IA signal for spiral galaxies under a quadratic tidal torquing model. Lastly, we measure the alignment between central galaxies and their host dark-matter halos, finding small to moderate misalignments between their principal axes that decline with halo mass.

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn.

The high-multiplicity exoplanet systems are generally more tightly packed when compared to the solar system. Such compact multi-planet systems are often susceptible to dynamical instability. We investigate the impact of dynamical instability on the final orbital architectures of multi-planet systems using N-body simulations. Our models initially consist of eight planets placed randomly according to a power-law distribution of mutual Hill separations. We find that almost all of our model planetary systems go through at least one phase of dynamical instability, losing at least one planet. The orbital architecture, including the distributions of mutual Hill separations, planetary masses, orbital periods, and period ratios, of the transit-detectable model planetary systems closely resemble those for the multi-planet systems detected by Kepler. We find that without any formation-dependent input, a dynamically active past can naturally reproduce important observed trends including multiplicity-dependent eccentricity distribution, smaller eccentricities for larger planets, and intra-system uniformity. These findings indicate that dynamical instabilities may have played a vital role in the final assembly of sub-Jovian planets.

RR Lyrae (RRL) are old, low-mass radially pulsating variable stars in their core helium burning phase. They are popular stellar tracers and primary distance indicators, since they obey to well defined period-luminosity relations in the near-infrared regime. Their photometric identification is not trivial, indeed, RRL samples can be contaminated by eclipsing binaries, especially in large datasets produced by fully automatic pipelines. Interpretable machine-learning approaches for separating eclipsing binaries from RRL are thus needed. Ideally, they should be able to achieve high precision in identifying RRL while generalizing to new data from different instruments. In this paper, we train a simple logistic regression classifier on Catalina Sky Survey (CSS) light curves. It achieves a precision of 87% at 78% recall for the RRL class on unseen CSS light curves. It generalizes on out-of-sample data (ASAS/ASAS-SN light curves) with a precision of 85% at 96% recall. We also considered a L1-regularized version of our classifier, which reaches 90% sparsity in the light-curve features with a limited trade-off in accuracy on our CSS validation set and -- remarkably -- also on the ASAS/ASAS-SN light curve test set. Logistic regression is natively interpretable, and regularization allows us to point out the parts of the light curves that matter the most in classification. We thus achieved both good generalization and full interpretability.

Filamentary structures are often identified in column density maps of molecular clouds, and appear to be important for both low- and high-mass star formation. Theoretically, these structures are expected to form in regions where the supersonic cloud-scale turbulent velocity field converges. While this model of filament formation successfully reproduces several of their properties derived from column densities, it is unclear whether it can also reproduce their kinematic features. We use a combination of hydrodynamical, chemical and radiative transfer modelling to predict the emission properties of these dynamically-forming filaments in the $^{13}$CO, HCN and N$_2$H$^+$ $J=1-0$ rotational lines. The results are largely in agreement with observations; in particular, line widths are typically subsonic to transonic, even for filaments which have formed from highly supersonic inflows. If the observed filaments are formed dynamically, as our results suggest, no equilibrium analysis is possible, and simulations which presuppose the existence of a filament are likely to produce unrealistic results.

Identification of Gamma Ray Burst (GRB) progenitors based on the duration of their prompt emission ($T_{90}$) has faced several roadblocks recently. Long-duration GRBs (with $T_{90} > 2s$) have traditionally been thought to be originating from the collapse of massive stars, and the short-duration ones (with $T_{90} < 2s$) from compact binary mergers. However, recent observations of a long GRB associated with a kilonova (KN) and a short GRB with supernova (SN) association demand a more detailed classification of the GRB population. In this {\it Letter}, we focus on GRBs associated with KNe, believed to be originating from mergers of binaries involving neutron stars (NS). We make use of the GRB prompt emission light curves of {\it Swift}-BAT 2022 GRB catalog and employ machine learning algorithms to study the classification of GRB progenitors. Our analysis reveals that there are five distinct clusters of GRBs, of which the KN-associated GRBs are located in two separate clusters indicating they may have been produced by different progenitors. We argue that these clusters may be due to subclasses of binary neutron star (BNS) and/or neutron star--black hole (NS-BH) mergers. We also discuss the implications of these findings for future gravitational-wave (GW) observations and how those observations may help in understanding these clusters better.

Early observations of transient explosions can provide vital clues to their progenitor origins. In this paper we present the nearby Type Iax (02cx-like) supernova (SN), SN 2020udy that was discovered within hours ($\sim$7 hr) of estimated first light. An extensive dataset of ultra-violet, optical, and near-infrared observations was obtained, covering out to $\sim$150 d after explosion. SN 2020udy peaked at -17.86$\pm$0.43 mag in the r band and evolved similarly to other 'luminous' SNe Iax, such as SNe 2005hk and 2012Z. Its well-sampled early light curve allows strict limits on companion interaction to be placed. Main-sequence companion stars with masses of 2 and 6 M$_\odot$ are ruled out at all viewing angles, while a helium-star companion is allowed from a narrow range of angles (140-180$^\circ$ away from the companion). The spectra and light curves of SN2020udy are in good agreement with those of the 'N5def' deflagration model of a near Chandrasekhar-mass carbon-oxygen white dwarf. However, as has been seen in previous studies of similar luminosity events, SN 2020udy evolves slower than the model. Broad-band linear polarisation measurements taken at and after peak are consistent with no polarisation, in agreement with the predictions of the companion-star configuration from the early light curve measurements. The host galaxy environment is low metallicity and is consistent with a young stellar population. Overall, we find the most plausible explosion scenario to be the incomplete disruption of a CO white dwarf near the Chandrasekhar-mass limit, with a helium-star companion.

Disks of low-mass bodies scattered by giant planets to large semi-major axis and constant periapsis orbits are vulnerable to a buckling instability. This instability exponentially grows orbital inclinations, raises periapsis distances, and coherently tilts orbits resulting in clustering of arguments of periapsis. The dynamically hot system is then susceptible to the formation of a lopsided mode. Here we show that the timescale of the buckling instability decreases as the radial surface density of the population becomes more centrally dense, i.e., steeper scattered disks buckle faster. Accounting for differential apsidal precession driven by giant planets, we find that $\sim\!10\,M_\oplus$ is sufficient for a primordial scattered disk in the trans-Neptunian region to have been unstable if $dN \propto a^{-2.5} da$.

Context. Current research has established magnetised disc winds as a promising way of driving accretion in protoplanetary discs. Aims. We investigate the evolution of large protoplanetary disc populations under the influence of magnetically driven disc winds as well as internal and external photoevaporation. We aim to constrain magnetic disc wind models through comparisons with observations. Methods. We ran 1D vertically integrated evolutionary simulations for low-viscosity discs, including magnetic braking and various outflows. The initial conditions were varied and chosen to produce populations that are representative of actual disc populations inferred from observations. We then compared the observables from the simulations (e.g. stellar accretion rate, disc mass evolution, disc lifetime, etc.) with observational data. Results. Our simulations show that to reach stellar accretion rates comparable to those found by observations $\sim 10^{-8}\mathrm{M}_\odot / \mathrm{yr}$, it is necessary to have access not only to strong magnetic torques, but weak magnetic winds as well. The presence of a strong magnetic disc wind, in combination with internal photoevaporation, leads to the rapid opening of an inner cavity early on, allowing the stellar accretion rate to drop while the disc is still massive. Furthermore, our model supports the notion that external photoevaporation via the ambient far-ultraviolet radiation of surrounding stars is a driving force in disc evolution and could potentially exert a strong influence on planetary formation. Conclusions. Our disc population syntheses show that for a subset of magnetohydrodynamic wind models (weak disc wind, strong torque), it is possible to reproduce important statistical observational constraints. The magnetic disc wind paradigm thus represents a novel and appealing alternative to the classical $\alpha$-viscosity scenario.

Molecules and particles make up $\sim 40 - 70\%$ of carbon in the interstellar medium, yet the exact chemical structure of these constituents remains unknown. We present carbon K-shell absorption spectroscopy of the Galactic Interstellar Medium obtained with the Low Energy Transmission Grating Spectrometer on the {\it Chandra} Observatory, that directly addresses this question. We probe several lines of sight, using bright AGN as backlighters. We make our measurements differentially with respect to the bright source Mrk 421, in order to take the significant carbon K absorption in the instrument into account. In the spectrum of the blazar 1ES 1553+113 we find evidence for a novel feature: strong extinction on the low-energy side of the neutral C $1s-2p$ resonance, which is indicative of scattering by graphite particles. We find evidence for characteristic particle radii of order $0.1-0.15$ $\mu$m. If this explanation for the feature is correct, limits on the mass of the available carbon along the line of sight may imply that the grains are partially aligned, and the X-rays from the source may have intrinsic polarization.

The incidence of young but fading radio sources provides important information on the life cycle of radio emission in radio-loud active galactic nuclei. Despite its importance for constraining the models of radio source evolution, there are no systematic studies of remnants in complete samples of young radio sources. We report results on the study of 18 compact steep-spectrum (CSS) radio sources, selected from the statistically complete B3-VLA CSS sample, characterized by a steep optically-thin spectrum (alpha > 1.0) and no core detection in earlier studies. Our deep multi-frequency Very Large Array (VLA), pc-scale Very Long Baseline Array (VLBA), and eMERLIN observations allowed us to locate the core component in 10 objects. In 3 CSS sources there is no clear evidence of present-time active regions, suggesting they are likely in a remnant phase. Among sources with core detection, we find 3 objects that have no clear active regions (hotspots) at the edges of the radio structure, suggesting that the radio emission may have just restarted. Our results support a power-law distribution of the source ages, although the poor statistics prevents us from setting solid constraints on the percentage of remnants and restarted sources in sub-populations of radio sources.

Despite its crucial role in galaxy evolution, the complex circumgalactic medium (CGM) remains underexplored. Although it is known to be multi-phase, the importance of the molecular gas phase to the total CGM mass budget is, to date, unconstrained. We present the first constraints on the molecular gas covering fraction in the CGM of low-redshift galaxies, using measurements of CO column densities along sightlines towards mm-bright background quasars with intervening galaxies. We do not detect molecular absorption against the background quasars. For the individual, low-redshift, 'normal' galaxy haloes probed here, we can therefore rule out the presence of an extremely molecular gas-rich CGM, as recently reported in high-redshift protoclusters and around luminous active galactic nuclei. We also set statistical limits on the volume filling factor of molecular material in the CGM as a whole, and as a function of radius. ISM-like molecular clouds of ~30 pc in radius with column densities of N(CO) >~ 10^16 cm^-2 have volume filling factors of less than 0.2 per cent. Large-scale smooth gas reservoirs are ruled out much more stringently. The development of this technique in the future will allow deeper constraining limits to be set on the importance (or unimportance) of molecular gas in the CGM.

Goal 1 of the National Academies of Science, Engineering and Mathematics Exoplanet Science Strategy is "to understand the formation and evolution of planetary systems as products of the process of star formation, and characterize and explain the diversity of planetary system architectures, planetary compositions, and planetary environments produced by these processes", with the finding that "Current knowledge of the demographics and characteristics of planets and their systems is substantially incomplete." One significant roadblock to our ongoing efforts to improve our demographics analyses is the lack of comprehensive meta-data accompanying published exoplanet surveys. The Exoplanet Program Analysis Group (ExoPAG) Science Interest Group 2: Exoplanet Demographics has prepared this document to provide guidance to survey architects, authors, referees and funding agencies as to the most valuable such data products for five different exoplanet detection techniques - transit, radial velocity, direct imaging, microlensing and astrometry. We find that making these additional data easily available would greatly enhance the community's ability to perform robust, reproducible demographics analyses, and make progress on achieving the most important goals identified by the exoplanet and wider astronomical community.

Molecular oxygen is a strong indicator of life on Earth, and may indicate biological processes on exoplanets too. Recent studies proposed that Earth-like O$_\mathrm{2}$ levels might be detectable on nearby exoplanets using high-resolution spectrographs on future extremely large telescopes (ELTs). However, these studies did not consider constraints like relative velocities, planet occurrence rates, and target observability. We expanded on past studies by creating a homogeneous catalog of 286,391 main-sequence stars within 120 pc using Gaia DR3, and used the Bioverse framework to simulate the likelihood of finding nearby transiting Earth analogs. We also simulated a survey of M dwarfs within 20 pc accounting for $\eta_{\oplus}$ estimates, transit probabilities, relative velocities, and target observability to determine how long ELTs and theoretical 50-100 meter ground-based telescopes need to observe to probe for Earth-like O$_\mathrm{2}$ levels with an $R=100,000$ spectrograph. This would only be possible within 50 years for up to $\sim$21% of nearby M dwarf systems if a suitable transiting habitable zone Earth-analog was discovered, assuming signals from every observable partial transit from each ELT can be combined. If so, Earth-like O$_\mathrm{2}$ levels could be detectable on TRAPPIST-1 d-g within 16 to 55 years, respectively, and about half that time with an $R=500,000$ spectrograph. These results have important implications for whether ELTs can survey nearby habitable zone Earth analogs for O$_\mathrm{2}$ via transmission spectroscopy. Our work provides the most comprehensive assessment to date of the ground-based capabilities to search for life beyond the solar system.

We present model spectral energy distribution (SED) fits to ultraviolet/optical/infrared observations for the 258 nearby galaxies in the Local Volume Legacy survey, a sample dominated by lower-luminosity dwarf irregular systems. The data for each galaxy include up to 26 spatially-integrated broadband and narrowband fluxes from the Galaxy Evolution Explorer, Spitzer Space Telescope, and Infrared Astronomical Satellite space-based platforms and from the Sloan Digital Sky Survey, Two Micron All Sky Survey, and other ground-based efforts. The CIGALE SED fitting package is employed using a delayed star formation history with an optional late burst or quenching episode to constrain 11 different free parameters that characterize the properties of each galaxy's stellar and dust emission, with the overriding constraint that the ultraviolet/optical emission absorbed by interstellar dust grains is emitted in equal energy portions at infrared wavelengths. The main results are: i) 94% of the SED fits yield reduced chi^2 values less than 3; ii) the modeled stellar masses agree with those derived from 3.6um-based measures with a scatter of 0.07 dex; iii) for a typical galaxy in the sample the SED-derived star formation rate averaged over the past 100 Myr is about 88% of the value derived from standard hybrid indicators on similar timescales; and iv) there is a statistically significant inverse relation between the stellar mass fraction appearing in the late burst and the total stellar mass. These results build upon prior SED modeling efforts in the local volume and lay the groundwork for future studies of more distant low-metallicity galaxies with JWST.

Close encounters between stars in star forming regions are important as they can perturb or destroy protoplanetary discs, young planetary systems, and stellar multiple systems. We simulate simple, viralised, equal-mass $N$-body star clusters and find that both the rate and total number of encounters between stars varies by factors of several in statistically identical clusters due to the stochastic/chaotic details of orbits and stellar dynamics. Encounters tend to rapidly `saturate' in the core of a cluster, with stars there each having many encounters, while more distant stars have none. However, we find that the fraction of stars that have had at least one encounter within a particular distance grows in the same way (scaling with crossing time and half-mass radius) in all clusters, and we present a new (empirical) way of estimating the fraction of stars that have had at least one encounter at a particular distance.

This work advances the (galaxy morphology)-dependent (black hole mass, $M_{\rm bh}$)-(spheroid/galaxy stellar mass, $M_*$) scaling relations by introducing `dust bins' for lenticular (S0) galaxies. Doing so has led to the discovery of $M_{\rm bh}$-$M_{\rm *,sph}$ and $M_{\rm bh}$-$M_{\rm *,gal}$ relations for dusty S0 galaxies - built by major wet mergers and comprising half the S0 sample - offset from the distribution of dust-poor S0 galaxies. The situation is reminiscent of how major dry mergers of massive S0 galaxies have created an offset population of ellicular and elliptical galaxies. For a given $M_{\rm bh}$, the dust-rich S0 galaxies have 3 to 4 times higher $M_{\rm *,sph}$ than the dust-poor S0 galaxies, and the steep distributions of both populations in the $M_{\rm bh}$-$M_{\rm *,sph}$ diagram bracket the $M_{\rm bh} \propto M_{\rm *,sph}^{2.27+/-0.48}$ relation defined by the spiral galaxies, themselves renovated through minor mergers. The new relations offer refined means to estimate $M_{\rm bh}$ in other galaxies and should aid with: (i) constructing (galaxy morphology)-dependent black hole mass functions; (ii) estimating the masses of black holes associated with tidal disruption events; (iii) better quantifying evolution in the scaling relations via improved comparisons with high-$z$ data by alleviating the pickle of apples versus oranges; (iv) mergers and long-wavelength gravitational wave science; (v) simulations of galaxy/black hole coevolution and semi-analytic works involving galaxy speciation; plus (vi) facilitating improved extrapolations into the intermediate-mass black hole landscape. The role of the galaxy's environment is also discussed, and many potential projects that can further explore the morphological divisions are mentioned.

The E-type $\alpha$-attractor models of single-field inflation were generalized further in order to accommodate production of primordial black holes (PBH) via adding a near-inflection point to the inflaton scalar potential at smaller scales, in good agreement with measurements of the cosmic microwave background (CMB) radiation. A minimal number of new parameters was used but their fine-tuning was maximized in order to increase possible masses of PBH formed during an ultra-slow-roll phase leading to a large enhancement of the power spectrum of scalar (curvature) perturbations by 6 or 7 orders of magnitude against the power spectrum of perturbations observed in CMB. It was found that extreme fine-tuning of the parameters in our models can lead to a formation of the Earth-size PBH with the masses of approximately $10^{27}$ g, still in agreement with CMB observations. Quantum corrections are known to lead to the perturbative upper bound on the amplitude of large scalar perturbations responsible for PBH production. The quantum (one-loop) corrections in our models were found to be suppressed by one order of magnitude for PBH with the masses of approximately $10^{19}$ g, which may form the whole dark matter in the Universe.

Diffuse gamma-ray emission (DGE) has been discovered over the Galactic disk in the energy range from sub-GeV to sub-PeV. While it is believed to be dominated by the pionic emission of cosmic ray (CR) hadrons via interactions with interstellar medium, unresolved gamma-ray sources may also be potential contributors. TeV gamma-ray halos around middle-aged pulsars have been proposed as such sources. Their contribution to DGE, however, highly depends on the injection rate of electrons and the injection spectral shape, which are not well determined based on current observations. The measured fluxes of DGE can thus provide constraints on the $e^\pm$ injection of the pulsar halo population in turn. In this paper, we estimate the contribution of pulsar halos to DGE based on the ATNF pulsar samples with taking into account the off-beamed pulsars. The recent measurement on DGE by Tibet AS$\gamma$ and an early measurement by MILAGRO are used to constrain the pair injection parameters of the pulsar halo population. Our result may be used to distinguish different models for pulsar halos.

Recent analysis on the stellar kinematics of ultra-faint dwarf (UFD) galaxies has put a stringent upper limit on the self-scattering cross section of dark matter, i.e., $\sigma/m<{\cal O}(0.1)\,{\rm cm^2/g}$ at the scattering velocity of ${\cal O}(10)\,{\rm km/s}$. Resonant self-interacting dark matter (rSIDM) is one possibility that can be consistent with the UFDs and explain the low central densities of rotation-supported galaxies; the cross section is resonantly enhanced to be $\sigma/m = {\cal O}(1)\,{\rm cm^2/g}$ around the scattering velocity of ${\cal O}(100)\,{\rm km/s}$ while being suppressed at lower velocities. To further assess this possibility, since the inferred dark matter distribution of halos from astrophysical observations is usually compared to that in constant-cross section SIDM (cSIDM), whether the structures of rSIDM halos can be approximated by the cSIDM halo profiles needs to be clarified. In this work, we employ the grovothermal fluid method to investigate the structural evolution of rSIDM halos in a wide mass range. We find that except for halos in a specific mass range, the present structures of rSIDM halos are virtually indistinguishable from those of the cSIDM halos. For halos in the specific mass range, the resonant self-scattering renders a break in their density profile. We demonstrate how such a density-profile break appears in astrophysical observations, e.g., rotation curves and line-of-sight velocity dispersion profiles. We show that for halos above the specific mass range, the density-profile break thermalizes to disappear before the present. We demonstrate that such distinctive thermalization dynamics can leave imprints on the orbital classes of stars with similar ages and metallicities.

A tidal disruption event (TDE) occurs when a star is destroyed by a supermassive black hole. Broadband radio spectral observations of TDEs trace the emission from any outflows or jets that are ejected from the vicinity of the supermassive black hole. However, radio detections of TDEs are rare, with less than 20 published to date, and only 11 with multi-epoch broadband coverage. Here we present the radio detection of the TDE AT2020vwl and our subsequent radio monitoring campaign of the outflow that was produced, spanning 1.5 years post-optical flare. We tracked the outflow evolution as it expanded between $10^{16}$ cm to $10^{17}$ cm from the supermassive black hole, deducing it was non-relativistic and launched quasi-simultaneously with the initial optical detection through modelling the evolving synchrotron spectra of the event. We deduce that the outflow is likely to have been launched by material ejected from stream-stream collisions (more likely), the unbound debris stream, or an accretion-induced wind or jet from the supermassive black hole (less likely). AT2020vwl joins a growing number of TDEs with well-characterised prompt radio emission, with future timely radio observations of TDEs required to fully understand the mechanism that produces this type of radio emission in TDEs.

The IceCube Neutrino Observatory sends realtime neutrino alerts with high probability of being astrophysical in origin. We present a new method to correlate these events and possible candidate sources using $2,089$ blazars from the Fermi-LAT 4LAC-DR2 catalog and with $3,413$ AGNs from the Radio Fundamental Catalog. No statistically significant neutrino emission was found in any of the catalog searches. The result is compatible with a small fraction, $<1$%, of AGNs being neutrino emitters and prior evidence for neutrino emission presented by IceCube and other authors from sources such as TXS 0506+056 and PKS 1502+06. We also present cross-checks to other analyses that claim a significant correlation using similar data samples, and we find that adding more information on the neutrino events and more data overall makes the result compatible with background.

Powerful winds at accretion disk scales have been observed in the past 20 years in many AGN, the so called Ultra-Fast Outflows (UFOs). Outflows are intimately related to mass accretion due to the conservation of angular momentum, and therefore are a key ingredient of most accretion disk models around BHs. At the same time, nuclear winds and outflows can provide the feedback which regulates the joint BH and galaxy growth. We reconsider UFO observations in the framework of the Magneto-Hydrodynamic Disk Wind (MHDW) scenario and study their statistical properties. We derive the typical wind-activity history in our sources by assuming that it can be statistically described by population functions. We study the statistical properties of UFOs from the literature and derive the distribution functions of the ratio $\bar \omega$ between the mass outflow and inflow rates, and the ratio $\lambda_w$ between the mass outflow and the Eddington accretion rates. We study the links between $\bar \omega$ and $\lambda_w$ and the Eddington ratio $\lambda={L_{bol}}/{L_{Edd}}$. We find that the distribution functions of $\bar \omega$ and $\lambda_w$ can be described as power laws above some threshold, suggesting that there may be many wind sub-events for each major wind event in each AGN activity cycle, which is a fractal behaviour in agreement with current MHDW and Chaotic Cold Accretion theories. We then introduce a simple cellular automaton to investigate how the dynamical properties of an idealized disk-wind system changes following the introduction of simple feedback rules. We find that without feedback the system is over-critical. Conversely, if feedback is present, the system can be driven toward self organized criticality. Our results corroborate the hypothesis that AGN feedback is a necessary key ingredient in disk-wind systems, and thus, in shaping the co-evolution of galaxies and supermassive BHs.

We model the evaporation histories of the three planets around K2-136, a K-dwarf in the Hyades open cluster with an age of 700 Myr. The star hosts three transiting planets, with radii of 1.0, 3.0 and 1.5 Earth radii, where the middle planet lies above the radius-period valley and the inner and outer planets are below. We use an XMM-Newton observation to measure the XUV radiation environment of the planets, finding that the X-ray activity of K2-136 is lower than predicted by models but typical of similar Hyades members. We estimate the internal structure of each planet, and model their evaporation histories using a range of structure and atmospheric escape formulations. While the precise X-ray irradiation history of the system may be uncertain, we exploit the fact that the three planets must have shared the same history. We find that the Earth-sized K2-136b is most likely rocky, with any primordial gaseous envelope being lost within a few Myr. The sub-Neptune, K2-136c, has an envelope contributing 1-1.7% of its mass that is stable against evaporation thanks to the high mass of its rocky core, whilst the super-Earth, K2-136d, must have a mass at the upper end of the allowed range in order to retain any of its envelope. Our results are consistent with all three planets beginning as sub-Neptunes that have since been sculpted by atmospheric evaporation to their current states, stripping the envelope from planet b and removing most from planet d whilst preserving planet c above the radius-period valley.

The field of radio astronomy is witnessing a boom in the amount of data produced per day due to newly commissioned radio telescopes. One of the most crucial problems in this field is the automatic classification of extragalactic radio sources based on their morphologies. Most recent contributions in the field of morphological classification of extragalactic radio sources have proposed classifiers based on convolutional neural networks. Alternatively, this work proposes gradient boosting machine learning methods accompanied by principal component analysis as data-efficient alternatives to convolutional neural networks. Recent findings have shown the efficacy of gradient boosting methods in outperforming deep learning methods for classification problems with tabular data. The gradient boosting methods considered in this work are based on the XGBoost, LightGBM, and CatBoost implementations. This work also studies the effect of dataset size on classifier performance. A three-class classification problem is considered in this work based on the three main Fanaroff-Riley classes: class 0, class I, and class II, using radio sources from the Best-Heckman sample. All three proposed gradient boosting methods outperformed a state-of-the-art convolutional neural networks-based classifier using less than a quarter of the number of images, with CatBoost having the highest accuracy. This was mainly due to the superior accuracy of gradient boosting methods in classifying Fanaroff-Riley class II sources, with 3--4\% higher recall.

NGC 1313\,X-1 is a mysterious Ultra-luminous X-ray (ULX) source whose X-ray powering mechanism and a bubble-like structure surrounding the source are topics of intense study. Here, we perform the X-ray spectroscopic study of the source using a joint {\it XMM-Newton} and {\it NuSTAR} observations taken during 2012 $-$ 2017. The combined spectra cover the energy band 0.3 $-$ 20 keV. We use the accretion-ejection-based JeTCAF model for spectral analysis. The model fitted disc mass accretion rate varies from 4.6 to 9.6 $\dot M_{\rm Edd}$ and the halo mass accretion rate varies from 4.0 to 6.1 $\dot M_{\rm Edd}$ with a dynamic Comptonizing corona of average size of $\sim 15$ $r_g$. The data fitting is carried out for different black hole (BH) mass values. The goodness of the fit and uncertainties in model parameters improve while using higher BH mass with most probable mass of the compact object to be $133\pm33$ M$_\odot$. We have estimated the mass outflow rate, its velocity and power, and the age of the inflated bubble surrounding the source. Our estimated bubble morphology is in accord with the observed optical bubble and winds found through high-resolution X-ray spectroscopy, suggesting that the bubble expanded by the outflows originating from the central source. Finally, we conclude that the super-Eddington accretion onto a nearly intermediate mass BH may power a ULX when the accretion efficiency is low, though their efficiency increases when jet/outflow is taken into account, in agreement with numerical simulations in the literature.

This work presents the charge-coupled device (CCD) photometric survey of the old open cluster NGC 188. Time-series V-band photometric observations were conducted for ten nights in January 2017 using the Nanshan One-meter Wide-field Telescope (NOWT) to search for variable stars in the field of the cluster field. A total of 25 variable stars, including one new variable star, were detected in the target field. Among the detected variables, 16 are cluster member stars, and the others are identified as field stars. The periods, radial velocities, effective temperatures, and classifications of the detected variables are discussed in this work. Most of the stars' effective temperatures are between 4200 K and 6600 K, indicating their spectral types are G or K. The newly discovered variable is probably a W UMa system. In this study, a known cluster variable star (V21 = V0769 Cep) is classified as an EA-type variable star based on the presence of an 0.5 magnitude eclipse in its light curve.

Large, energy-dependent X-ray polarisation is observed in 4U 1630-47, a black hole in an X-ray binary, in the high-soft emission state. In this state, X-ray emission is believed to be dominated by a thermal, geometrically thin, optically thick accretion disc. However, the observations with the Imaging X-ray Polarimetry Explorer (IXPE) reveal an unexpectedly high polarisation degree, rising from 6% at 2 keV to 10% at 8 keV, which cannot be reconciled with standard models of thin accretion discs. We argue that an accretion disc with an only partially ionised atmosphere flowing away from the disc at mildly relativistic velocities can explain the observations.

Observations point to a correlation between outer giants and inner sub-Neptunes, unexplained by simulations so far. We utilize N-body simulations including pebble and gas accretion as well as planetary migration to investigate how the gas accretion rates influence the formation of systems of inner sub-Neptunes and outer gas giants as well as the eccentricity distribution of the outer giant planets. Less efficient envelope contraction rates allow a more efficient formation of systems with inner sub-Neptunes and outer giants. This is caused by the fact that the cores formed in the inner disc are too small to accrete large envelopes and only cores growing in the outer disc can become giants. As a result, instabilities between the outer giant planets do not necessarily destroy the inner systems of sub-Neptunes unlike simulations where giant planets can form closer in. Our simulations show that up to 50% of the systems of cold Jupiters could have inner sub-Neptunes, in agreement with observations. Our simulations show a good agreement with the eccentricity distribution of giants, even though we find a slight mismatch to the mass and semi-major axes distributions. Synthetic transit observations of the inner systems (r<0.7 AU) reveal an excellent match to the Kepler observations, where our simulations match the period ratios of adjacent planet pairs. Thus, the breaking the chains model for super-Earth and sub-Neptune formation remains consistent with observations even when outer giant planets are present. However, simulations with outer giant planets produce more systems with mostly only one inner planet and with larger eccentricities, in contrast to simulations without outer giants. We thus predict that systems with truly single close-in planets are more likely to host outer gas giants and we consequently suggest RV follow-up observations of these systems to constrain the formation pathway.

The characterisation of giant exoplanets is crucial to constrain giant planet formation and evolution theory and for putting the solar-system's giant planets in perspective. Typically, mass-radius (M-R) measurements of moderately irradiated warm Jupiters are used to estimate the planetary bulk composition, which is an essential quantity for constraining giant planet formation, evolution and structure models. The successful launch of the James Webb Space Telescope (JWST) and the upcoming ARIEL mission open a new era in giant exoplanet characterisation as atmospheric measurements provide key information on the composition and internal structure of giant exoplanets. In this review, we discuss how giant planet evolution models are used to infer the planetary bulk composition, and the connection between the compositions of the interior and atmosphere. We identify the important theoretical uncertainties in evolution models including the equations of state, atmospheric models, chemical composition, interior structure and main energy transport processes. Nevertheless, we show that that atmospheric measurements by JWST and ARIEL and the accurate determination of stellar ages by PLATO can significantly reduce the degeneracy in the inferred bulk composition. Furthermore, we discuss the importance of evolution models for the characterisation of direct-imaged planets. We conclude that giant planet theory has a critical role in the interpretation of observation and emphasise the importance of advancing giant planet theory.

In the core accretion scenario of planet formation, rocky cores grow by first accreting solids until they are massive enough to accrete gas. For giant planet formation this means that a massive core must form within the lifetime of the gas disk. The accretion of roughly km-sized planetesimals and the accretion of mm-cm sized pebbles are typically discussed separately as the main solid accretion mechanisms. We investigate the interplay between the two accretion processes in a disk containing both pebbles and planetesimals for planet formation in general and in the context of giant planet formation specifically. The goal is to disentangle and understand the fundamental interactions that arise in such hybrid pebble-planetesimal models. We combine a simple model of pebble formation and accretion with a global model of planet formation which considers the accretion of planetesimals. We compare synthetic populations of planets formed in disks composed of different amounts of pebbles and 600 meter sized planetesimals. On a system-level, we study the formation pathway of giant planets in these disks. We find that, in hybrid disks containing both pebbles and planetesimals, the formation of giant planets is strongly suppressed whereas in a pebbles-only or planetesimals-only scenario, giant planets can form. We identify the heating associated with the accretion of up to 100 km sized planetesimals after the pebble accretion period to delay the runaway gas accretion of massive cores. Coupled with strong inward type-I migration acting on these planets, this results in close-in icy sub-Neptunes originating from the outer disk. We conclude that, in hybrid pebble-planetesimal scenarios, the late accretion of planetesimals is a critical factor in the giant planet formation process and that inward migration is more efficient for planets in increasingly pebble dominated disks.

A first hydrostatic core (FHC) is proposed to form after the initial collapse of a prestellar core, as a seed of a Class 0 protostar. FHCs are difficult to observe because they are small, compact, embedded, and short lived. In this work, we explored the physical properties of two well-known FHC candidates, B1-bN and B1-bS, by comparing interferometric data from Submillimeter Array (SMA) 1.1 and 1.3 mm and Atacama Large Millimeter/submillimeter Array (ALMA) 870 $\mu$m observations with simulated synthesis images of the two sources. The simulated images are based on a simple model containing a single, hot compact first-core-like component at the center surrounded by a large-scale, cold and dusty envelope described by a broken power-law density distribution with an index, $\alpha$. Our results show that the hot compact components of B1-bN and B1-bS can be described by temperatures of \sim 500 K with a size of \sim 4 au, which are in agreement with theoretical predictions of an FHC. If the $\alpha$ inside the broken radii is fixed to -1.5, we find $\alpha$ \sim -2.9 and \sim -3.3 outside the broken radii for B1-bN and B1-bS, respectively, consistent with theoretical calculations of a collapsing, bounded envelope and previous observations. Comparing the density and temperature profiles of the two sources with radiation-hydrodynamic simulations of an FHC, we find both sources lie close to, but before, the second collapse stage. We suggest that B1-bS may have started the collapsing process earlier compared to B1-bN, since a larger discontinuity point is found in its density profile.

We present a comprehensive timing and spectral analysis of the HMXB 4U 1538-522 by using the Nuclear Spectroscopic Telescope Array (NuSTAR) observatory data. Using three archived observations made between 2019 and 2021, we have detected $\sim $ 526 s coherent pulsations up to 60 keV. We have found an instantaneous spin-down rate of $\dot{P} = 6.6_{-6.0}^{+2.4} \times 10^{-6}$ s s$^{-1}$ during the first observation. The pulse profiles had a double peaked structure consisting of a broad primary peak and an energy dependent, weak secondary peak. We have also analysed the long-term spin-period evolution of 4U 1538-522 from data spanning more than four decades, including the data from Fermi/GBM. Based on the recent spin trends, we have found that the third torque reversal in 4U 1538-522 happened around MJD 58800. The source is currently spinning up with $\dot{P} = -1.9(1) \times 10^{-9}$ s s$^{-1}$. We also report a periodic fluctuation in the spin-period of 4U 1538-522. The broad-band persistent spectra can be described with a blackbody component and either powerlaw or Comptonization component along with a Fe K$_{\alpha}$ line at 6.4 keV and a cyclotron absorption feature around 22 keV. We have also found a relatively weak absorption feature around 27 keV in the persistent spectra of 4U 1538-522 in all three observations. We have estimated a magnetic field strength of $1.84_{-0.06}^{+0.04} (1+z) \times 10^{12}$ and $2.33_{-0.24}^{+0.15} (1+z) \times 10^{12}$ G for the two features, respectively.

Using the PMAS Integral Field Unit on the Calar Alto 3.5m telescope we observed the southern component (Markarian 59) of the `cometary' starburst galaxy NGC 4861. Mrk 59 is centred on a giant nebula and concentration of stars 1 kpc in diameter. Strong $\rm H\alpha$ emission points to a star-formation rate (SFR) at least 0.47 $\rm M_{\odot}yr^{-1}$. Mrk 59 has a very high [OIII]$\rm\lambda5007/H\beta$ ratio, reaching 7.35 in the central nebula, with a second peak at a star-forming hotspot further north. Fast outflows are not detected but nebular motion and galaxy rotation produce relative velocities up to 40 km $\rm s^{-1}$. Spectral analysis of different regions with `Fitting Analysis using Differential evolution Optimisation' (FADO) finds that the stars in the central and `spur' nebulae are very young, $\rm \leq125~Myr$ with a large $\rm <10~Myr$ contribution. Older stars ($\rm \sim 1~Gyr$), make up the northern disk component, while the other regions show mixtures of 1 Gyr age with very young stars. This and the high specific SFR $\rm\sim 3.5~Gyr^{-1}$ imply a bimodal star formation history, with Mrk 59 formed in ongoing starbursts fuelled by a huge gas inflow, turning the galaxy into an asymmetric `green pea' or blue compact dwarf. We map the HeII$\lambda4686$ emission, and identify a broad component from the central nebula, consistent with the emission of $\sim 300$ Wolf-Rayet stars. About a third of the HeII$\lambda$4686 flux is a narrow line emitted from a more extended area covering the central and spur nebulae, and may have a different origin.

We report detailed X-ray observations of the unique binary system $\epsilon$ Lupi, the only known short-period binary consisting of two magnetic early-type stars. The components have comparably strong, but anti-aligned magnetic fields. The orbital and magnetic properties of the system imply that the magnetospheres overlap at all orbital phases, suggesting the possibility of variable inter-star magnetospheric interaction due to the non-negligible eccentricity of the orbit. To investigate this effect, we observed the X-ray emission from $\epsilon$ Lupi both near and away from periastron passage, using the Neutron Star Interior Composition Explorer mission (NICER) X-ray Telescope. We find that the system produces excess X-ray emission at the periastron phase, suggesting the presence of variable inter-star magnetospheric interaction. We also discover that the enhancement at periastron is confined to a very narrow orbital phase range ($\approx 5\%$ of the orbital period), but the X-ray properties close to periastron phase are similar to those observed away from periastron. From these observations, we infer that the underlying cause is magnetic reconnection heating the stellar wind plasma, rather than shocks produced by wind-wind collision. Finally, by comparing the behavior of $\epsilon$ Lupi with that observed for cooler magnetic binary systems, we propose that elevated X-ray flux at periastron phase is likely a general characteristic of interacting magnetospheres irrespective of the spectral types of the constituent stars.

V404 Cyg is a Low Mass X-Ray Binary (LMXB) system that has undergone outbursts in 1938, 1989, and 2015. During these events, it has been possible to determine relevant data of the system; such as the masses of the compact object (a black hole, BH) and its companion, the orbital period, the companion spectral type, and luminosity class, among others. Remarkably, the companion star has a metallicity appreciably higher than solar. All these data allow us to construct theoretical models to account for its structure, looking for its initial configuration and predicting its final fate. Assuming that the BH is already formed when the primary star reaches the Zero Age Main Sequence, we used our binary evolution code for such a purpose. We obtained that the present characteristics of the system are nicely accounted for by a model with initial masses of 9 solar masses for the BH, 1.5 solar masses for the companion star, an initial orbital period of 1.5 d and considering that at most 30% of the mass transferred by the donor is accreted by the BH. The metallicity of the donor for our best fit was Z = 0.028 (twice solar metallicity). We also studied the evolution of the BH spin parameter assuming that initially, it is not rotating. Remarkably, the spin of the BHs in our models is far from reaching the available observational determination. This may indicate that the BH in V404 Cyg is initially spinning, a result that may be relevant for understanding the formation BHs in the context of LMXB systems.

Context. The origin of the large star-to-star variation of the [Eu/Fe] ratios observed in the extremely metal-poor (at [Fe/H]$\leq-3$) stars of the Galactic halo is still a matter of debate.\\ Aims. In this paper, we explore this problem by putting our stochastic chemical evolution model in the hierarchical clustering framework, with the aim of explaining the observed spread in the halo.\\ Methods. We compute the chemical enrichment of Eu occurring in the building blocks that have possibly formed the Galactic halo. In this framework, the enrichment from neutron star mergers can be influenced by the dynamics of the binary systems in the gravitational potential of the original host galaxy. In the least massive systems, the neutron stars can merge outside the host galaxy and so only a small fraction of newly produced Eu can be retained by the parent galaxy itself.\\ Results. In the framework of this new scenario, the accreted merging neutron stars are able to explain the presence of stars with sub-solar [Eu/Fe] ratios at [Fe/H]$\leq-3$, but only if we assume a delay time distribution for merging of the neutron stars $\propto t^{-1.5}$. We confirm the correlation between the dispersion of [Eu/Fe] at a given metallicity and the fraction of massive stars which give origin to neutron star mergers. The mixed scenario, where both neutron star mergers and magneto-rotational supernovae do produce Eu, can explain the observed spread in the Eu abundance also for a delay time distribution for mergers going either as $\propto t^{-1}$ or $\propto t^{-1.5}$.

We investigate the angular momentum of mono-abundance populations (MAPs) of the Milky Way thick disk by using a sample of 26,076 giant stars taken from APOGEE DR17 and Gaia EDR3. The vertical and perpendicular angular momentum components, $L_Z$ and $L_P$, of MAPs in narrow bins have significant variations across the [$\alpha$/M]-[M/H] plane. $L_Z$ and $L_P$ systematically change with [M/H] and [$\alpha$/M] and can be alternatively quantified by the chemical gradients: $d[{\rm M/H}]/dL_Z = 1.2\times 10^{-3} $\,dex\,kpc$^{-1}$\,km$^{-1}$\,s, $d{\rm [M/H]}/dL_P = -5.0\times 10^{-3}$\,dec\,kpc$^{-1}$\,km$^{-1}$\,s, and $d[\alpha/{\rm M}]/dL_Z = -3.0\times 10^{-4} $\,dex\,kpc$^{-1}$\,km$^{-1}$\,s, $d[\alpha/{\rm M}]/dL_P = 1.2\times 10^{-3}$\,dec\,kpc$^{-1}$\,km$^{-1}$\,s. These correlations can also be explained as the chemical-dependence of the spatial distribution shape of MAPs. We also exhibit the corresponding age dependence of angular momentum components. Under the assumption that the guiding radius ($R_g$) is proportional to $L_Z$, it provides direct observational evidence of the inside-out structure formation scenario of the thick disk, with $dR_g/dAge = -1.9$\,kpc\,Gyr$^{-1}$. The progressive changes in the disk thickness can be explained by the upside-down formation or/and the consequent kinematical heating.

Despite the universe containing primordial thorium (Th) of sufficient abundance to appear in stellar spectra, detection of Th has to date been tentative and based on just a few weak and blended lines. Here, we present convincing evidence not only for the first Th detection in a magnetic chemically peculiar Ap star but also for the first detection of Th III in a stellar spectrum. CPD-62 2717 was initially recognized as a highly-magnetized Ap star thanks to resolved magnetically split lines captured in $H$-band spectra from the SDSS/APOGEE survey. The star was subsequently pinpointed as extraordinarily peculiar when careful inspection of the $H$-band line content revealed the presence of five lines of Th III, none of which are detected in the other $\sim1500$ APOGEE-observed Ap stars. Follow-up with the VLT+UVES confirmed a similarly peculiar optical spectrum featuring dozens of Th III lines, among other peculiarities. Unlike past claims of Th detection, and owing to high-resolution observations of the strong ($\sim$8$-$12$\,$kG) magnetic field of CPD-62 2717, the detection of Th III can in this case be supported by matches between the observed and theoretical magnetic splitting patterns. Comparison of CPD-62 2717 to stars for which Th overabundances have been previously reported (e.g., Przybylski's Star) indicate that only for CPD-62 2717 is the Th detection certain. Along with the focus on Th III, we use time series measurements of the magnetic field modulus to constrain the rotation period of CPD-62 2717 to $\sim$4.8 years, thus establishing it as a new example of a super-slowly-rotating Ap star.

Asteroseismology has become a powerful tool to study the internal rotation of stars, and its study allows to constrain the internal AM transport processes and better understand their physical nature. In this context, we compared the rotation rates predicted by asteroseismology and by starspots measurements for four main-sequence stars from the Kepler LEGACY sample, considering different AM transport prescriptions, and investigated if some of these prescriptions could be ruled out. We decoupled the modelling of the structure and of the rotational profile, respectively obtained by an asteroseismic characterization and by using rotating models including a detailed treatment of the AM transport. We then compared the mean asteroseismic rotation rate with the surface rotation rate from starspots measurements for each of the AM transport prescriptions. In the hotter part of the HRD (M > ~ 1.2Msun), combining asteroseismic constraints from splittings of pressure modes and surface rotation rates does not allow to conclude on the need for an efficient AM transport in addition to the sole transport by meridional circulation and shear instability. Both prescriptions are indeed consistent with the quasi-solid rotation measured by Benomar et al. (2015) and Nielsen et al. (2017). In the colder part of the HRD, the situation is different due to the efficient braking of the stellar surface by magnetised winds. We find a clear disagreement between the rotational properties of models including only hydrodynamic processes and asteroseismic constraints, while models with magnetic fields correctly reproduce the observations, similarly to the solar case. This shows the existence of a mass regime corresponding to main-sequence stars around ~ 6000 - 6200 K for which it is difficult to constrain the AM transport processes, unlike for hotter, Gamma Dor stars or colder, less massive solar analogs.

$Context.$ The YSO Th 28 possesses a highly collimated jet, which clearly exhibits an asymmetric brightness of its jet lobes at optical and NIR wavelengths. There may be asymmetry in the jet plasma parameters in opposite jet lobes (e.g. electron density, temperature, and outflow velocity). $Aims.$ We examined the Th 28 jet in a 3"x3" where the jet material is collimated and accelerated. Our goal is to map the morphology and determine its physical parameters to determine the physical origin of such asymmetries. $Methods.$ We present $JHK$-spectra of Th 28 obtained with the SINFONI on the (VLT, ESO) in June-July 2015. $Results.$ The [Fe II] emission originates in collimated jet lobes. Two new axial knots are detected at 1" in the blue lobe and 1".2 in the red lobe. The H$_2$ radiation is emitted from an extended region with a radius of $\gtrsim270$ au, which is perpendicular to the jet. The PV diagrams of the bright H$_2$ lines reveal faint H$_2$ emission along both jet lobes as well. The compact and faint H I emission (Pa$\beta$ and Br$\gamma$) comes from two regions, namely from a spherical region around the star and from the jet lobes. The size of the jet launching region is derived as 0".015 ($\sim$3 au at 185 pc), and the initial opening angle of the Th 28 jet is $\sim28^0$, which makes this jet substantially less collimated than most jets from other CTTs. $Conclusions.$ The emission in [Fe II], H$_2$, and H I lines suggests a morphology in which the ionised gas in the disc appears to be disrupted by the jet. The resolved disc-like H$_2$ emission most likely arises in the disc atmosphere from shocks caused by a radial uncollimated wind. The asymmetry of the [Fe II] photocentre shifts with respect to the jet source arises in the immediate vicinity of the driving source of Th28 and suggests that the observed brightness asymmetry is intrinsic as well.

This paper describes the extended data release of the Calar Alto Legacy Integral Field Area (CALIFA) survey (eDR). It comprises science-grade quality data for 895 galaxies obtained with the PMAS/PPak instrument at the 3.5 m telescope at the Calar Alto Observatory along the last 12 years, using the V500 setup (3700-7500{\AA}, 6{\AA}/FWHM) and the CALIFA observing strategy. It includes galaxies of any morphological type, star-formation stage, a wide range of stellar masses ($\sim$10$^7$ 10$^{12}$ Msun ), at an average redshift of $\sim$0.015 (90\% within 0.005$<$z$<$0.05). Primarily selected based on the projected size and apparent magnitude, we demonstrate that it can be volume corrected resulting in a statistically limited but representative sample of the population of galaxies in the nearby Universe. All the data were homogeneous re-reduced, introducing a set of modifications to the previous reduction. The most relevant is the development and implementation of a new cube-reconstruction algorithm that provides with an (almost) seeing-limited spatial resolution (FWHM PSF $\sim$1.0").To illustrate the usability and quality of the data, we extracted two aperture spectra for each galaxy (central 1.5" and fully integrated), and analyze them using pyFIT3D. We obtain a set of observational and physical properties of both the stellar populations and the ionized gas, that have been compared for the two apertures, exploring their distributions as a function of the stellar masses and morphologies of the galaxies, comparing with recent results in the literature. DATA RELEASE: this http URL unam.mx/CALIFA_WEB/public_html/

We present first results from a novel numerical relativity code based on a tetrad formulation of the Einstein-scalar field equations combined with recently introduced gauge/frame invariant diagnostics indicating that inflation does not solve the homogeneity and isotropy problem beginning from generic initial conditions following a big bang.

R-mode oscillations of rotating neutron stars(NS) are promising candidates for continuous gravitational wave (GW) observations. In our recent work(Ghosh et al. 2023), we derived universal relations of the NS parameters, compactness and dimensionless tidal deformability with the r-mode frequency. In this work, we investigate how these universal relations can be used to infer various NS intrinsic parameters following a successful detection of the r-modes. In particular, we show that for targeted r-mode searches, these universal relations along with the "I-Love-Q" relation can be used to estimate both the moment of inertia and the distance of the NS thus breaking the degeneracy of distance measurement for continuous gravitational wave(CGW) observations. We also discuss that with a prior knowledge of the distance of the NS from electromagnetic observations, these universal relations can also be used to constrain the dense matter equation of state (EOS) inside NS. We quantify the accuracy to which such measurements can be done using the Fisher information matrix for a broad range of possible, unknown parameters, for both the a-LIGO and Einstein Telescope (ET) sensitivities.

We use the method of the superpotential to derive exact solutions describing inflationary cosmologies in multi-field models. An example that describes a solution that interpolates between two de Sitter universes is described in detail.

The scalar field sector in low--energy effective field theories motivated by string theory often contains several scalar fields, some of which possess non--standard kinetic terms. In this paper, we study theories with two scalar fields, in which one of the fields has a non--canonical kinetic term. The kinetic coupling is allowed to depend on both fields, going beyond the work in the literature, which usually considers the case of the coupling to depend on the other field only. Our aim is to study adiabatic and isocurvature perturbations in these extended theories. Our results show that the evolution equation for the curvature perturbation does not change when allowing the coupling to depend on both fields, while the effective mass of the entropy perturbation changes. We find expressions for the spectral index and its running at horizon crossing and at the end of inflation. We apply the formalism and study three phenomenological models, with different kinetic couplings.

Running vacuum models and viscous dark matter scenarios beyond perfect fluid idealization are two appealing theoretical strategies that have been separately studied as alternatives to solve some problems rooted in the $\Lambda$CDM cosmological model. In this paper, we combine these two notions in a single cosmological setting and investigate their cosmological implications, paying particular attention in the interplay between these two constituents in different cosmological periods. Specifically, we consider a well-studied running vacuum model inspired by renormalization group, and a recently proposed general parameterization for the bulk viscosity $\xi$. By employing dynamical system analysis, we explore the physical aspects of the new phase space that emerges from the combined models and derive stability conditions that ensure complete cosmological dynamics. We identify four distinct classes of models and find that the critical points of the phase space are non-trivially renewed compared to the single scenarios. We then proceed, in a joint and complementary way to the dynamical system analysis, with a detailed numerical exploration to quantify the impact of both the running parameter and the bulk viscosity coefficient on the cosmological evolution. Thus, for some values of the model parameters, numerical solutions show qualitative differences from the $\Lambda$CDM model, which is phenomenologically appealing in light of cosmological observations.

We discuss in detail the possibility that the "type-II majoron" -- that is, the pseudo Nambu-Goldstone boson that arises in the context of the type-II seesaw mechanism if the lepton number is spontaneously broken by an additional singlet scalar -- account for the dark matter (DM) observed in the universe. We study the requirements the model's parameters have to fulfill in order to reproduce the measured DM relic abundance through two possible production mechanisms in the early universe, freeze-in and misalignment, both during a standard radiation-dominated era and early matter domination. We then study possible signals of type-II majoron DM and the present and expected constraints on the parameter space that can be obtained from cosmological observations, direct detection experiments, and present and future searches for decaying DM at neutrino telescopes and cosmic-ray experiments. We find that -- depending on the majoron mass, the production mechanism, and the value of the vacuum expectation value of the type-II triplet -- all of the three decay modes (photons, electrons, neutrinos) of majoron DM particles can yield observable signals at future indirect searches for DM. Furthermore, in a corner of the parameter space, detection of majoron DM is possible through electron recoil at running and future direct detection experiments.

Since electroweak symmetry is generally broken during inflation, the Standard Model Higgs field can become supermassive even after the end of inflation. In this paper, we study the non-thermal phase space distribution of the Higgs field during reheating, focusing in particular on two different contributions: primordial condensate and stochastic fluctuations. We obtain their analytic formulae, which agree with the previous numerical result. As a possible consequence of the non-thermal Higgs spectrum, we discuss perturbative Higgs decay during reheating for the case it is kinematically allowed. We find that the soft-relativistic and hard spectra are dominant in the decay rate of the stochastic fluctuation and that the primordial condensate and stochastic fluctuations decay almost at the same time.

In this paper, we revise the constraints on the $f(Q)=Q/(8\pi G) - \alpha \ln(Q/Q_0)$, symmetric teleparallel model using local measurements and gravitational waves mock standard sirens. Using observational local SNIa and BAO data and energy conditions, the logarithmic $f(Q)$ model is capable of explaining the cosmic late-time acceleration by geometrical means. This result suggests that the logarithmic symmetric teleparallel model could be a candidate to solve the cosmological constant problem. In the case of the simulated standard siren data, by using the performance of the future ET and LISA detectors, we expect to be able to measure the current Hubble constant $H_0$, and the matter content $\Omega_m$, with a precision better than 1% and 6%, respectively. Furthermore, we explore the predicted $f(Q)$ logarithmic model deviation from the standard GR using ET and LISA mock standard sirens. The ratio $d_L^{\text{gw}}(z)/d_L^{\text{em}}(z)$, which quantifies the deviation from GR gives us a significant deviation higher than 13% at $z=1$, and it continues growing to reach a deviation higher than 18% in its median value. Future standard siren data will be able to quantify the strength of the deviation from GR and hence whether a cosmology like the one implied by this $f(Q)$ model is feasible.

Lorentz invariance is one of the foundations of modern physics; however, Lorentz violation may happen from the perspective of quantum gravity, and plenty of studies on Lorentz violation have arisen in recent years. As a good tool to explore Lorentz violation, Finsler geometry is a natural and fundamental generalization of Riemann geometry. The Finsler structure depends on both coordinates and velocities. Here, we simply introduce the mathematics of Finsler geometry. We review the connection between modified dispersion relations and Finsler geometries and discuss the physical influence from Finsler geometry. We review the connection between Finsler geometries and theories of Lorentz violation, such as the doubly special relativity, the standard-model extension, and the very special relativity.

Light weakly interacting particles could be copiously produced in the Sun which, as a well-understood star, could provide severe constraints on such new physics. In this work, we calculate the solar production rates of light gauge bosons (e.g. dark photon) arising from various $U(1)$ extensions of the standard model. It is known that the dark photon production rate is suppressed by the dark photon mass if it is well below the plasmon mass of the medium. We show that for more general $U(1)$ gauge bosons, this suppression is absent if the couplings are not in alignment with those of the photon. We investigate a few frequently discussed $U(1)$ models including $B-L$, $L_{\mu}-L_{\tau}$, and $L_{e}-L_{\mu(\tau)}$, and derive the stellar cooling bounds for these models.