Papers for Friday, Feb 28 2025

We present a spectroscopic study of the ``Disco Ball'' (SMDG0038365-064207), a rotationally-supported, red-sequence, ultra-diffuse galaxy (UDG) with a nuclear star cluster (NSC), multiple stellar clusters, and active star-forming regions using data obtained with KCWI on the Keck II Telescope. We calculate that the galaxy hosts $34\pm11$ ``globular" clusters. Kinematic measurements confirm rotation with a peak rotational velocity of at least 53km s$^{-1}$ and a dynamical mass within $r_{\rm e}$ of at least $10^{9.3 \pm 0.2}M_{\odot}$. Our dynamical estimates of the halo mass are consistent with that obtained using the number of globular clusters and together suggest $M_{\rm h}=10^{11.1\pm0.2}M_{\odot}$. The NSC may exhibit signatures of weak AGN activity. Our findings challenge two common assumptions: (1) clusters in some UDGs may be younger than generally assumed, and thus more luminous than standard globular clusters (GCs), affecting GC counts and the derived GC luminosity function in these UDGs, and (2) quiescent UDGs can be rotationally supported, making kinematic measurements viewing-angle dependent in such cases. The Disco Ball, while unremarkable in mass, size, projected structural properties, or color, reveals surprising complexity, highlighting the need for detailed studies of more UDGs.

We present ultraviolet to infrared observations of the extraordinary Type IIn supernova 2023zkd (SN 2023zkd). Photometrically, it exhibits persistent and luminous precursor emission spanning $\sim$4 years preceding discovery ($M_r\approx-15$ mag, 1,500~days in the observer frame), followed by a secondary stage of gradual brightening in its final year. Post-discovery, it exhibits two photometric peaks of comparable brightness ($M_r\lesssim-18.7$ mag and $M_r\approx-18.4$ mag, respectively) separated by 240 days. Spectroscopically, SN 2023zkd exhibits highly asymmetric and multi-component Balmer and He I profiles that we attribute to ejecta interaction with fast-moving ($1,\!000-2,\!000\;\mathrm{km}\;\mathrm{s}^{-1}$) He-rich polar material and slow-moving ($\sim$$400\;\mathrm{km}\;\mathrm{s}^{-1}$) equatorially-distributed H-rich material. He II features also appear during the second light curve peak and evolve rapidly. Shock-driven models fit to the multi-band photometry suggest that the event is powered by interaction with $\sim$$5-6\;M_{\odot}$ of CSM, with $2-3\;M_{\odot}$ associated with each light curve peak, expelled during mass-loss episodes $\sim$$3-4$ and $\sim$$1-2$ years prior to explosion. The observed precursor emission, combined with the extreme mass-loss rates required to power each light curve peak, favors either super-Eddington accretion onto a black hole or multiple long-lived eruptions from a massive star to luminosities that have not been previously observed. We consider multiple progenitor scenarios for SN 2023zkd, and find that the brightening optical precursor and inferred explosion properties are most consistent with a massive ($M_{\mathrm{ZAMS}}\geq30\;M_{\odot}$) and partially-stripped He star undergoing an instability-induced merger with a black hole companion.

The study of scaling relations of disc galaxies and their evolution across cosmic time requires accurate estimates of galaxy stellar masses $M_\star$ over broad redshift ranges. While photometric $M_\star$ estimates ($M_{\rm phot}$) based on spectral energy distribution (SED) modelling methods are employed routinely at high-$z$, it is unclear to what extent these are compatible with dynamical $M_\star$ estimates ($M_{\rm dyn}$), available for nearby galaxies. Here we compare newly determined, SED-model based $M_{\rm phot}$ with previously obtained $M_{\rm dyn}$ inferred via rotation curve decomposition techniques in a sample of $\sim100$ nearby galaxies from the SPARC database. We find that the two mass estimates show a systematic agreement at the $\sim12\%$ ($0.05$ dex) level and a $\sim55\%$ ($0.22$ dex) scatter across almost $5$ dex in $M_\star$. Our $M_{\rm phot}$ estimates correspond to mass-to-light ratios in the $3.6\mu$m band that increase gradually with $3.6\mu$m luminosity, as a consequence of the earlier (later) assembly history of high-mass (low-mass) disc galaxies. The choice of using either $M_{\rm dyn}$ or $M_{\rm phot}$ has only a marginal impact on the slope and zero-point of the Tully-Fisher and Fall relations: the observed orthogonal scatter in both relations is virtually the same for the two methods, and indistinguishable from that derived using a constant mass-to-light ratio in the $3.6\mu$m band. $M_\star$ estimates based on the assumption that discs are marginally stable lead to the largest scatter in the scaling relations.

The recent astrometric discovery of the nearby (590 pc) massive ($33 M_\odot$) dormant black hole candidate Gaia BH3 offers the possibility to angularly resolve the black hole from its companion star by using optical interferometry. Our aim is to detect emission in the near-infrared K band from the close-in environment of Gaia BH3 caused by accretion. Gaia BH3 was observed with the GRAVITY instrument using the four 8-meter Unit Telescopes of the VLT Interferometer. We searched for the signature of emission from the black hole in the interferometric data using the CANDID, PMOIRED, and exoGravity tools. With a present separation of 18 mas, the Gaia BH3 system can be well resolved angularly by GRAVITY. We did not detect emission from the black hole at a contrast level of $\Delta m = 6.8$ mag with respect to the companion star, that is, $f_\mathrm{BH}/f_* < 0.2\%$. This corresponds to an upper limit on the continuum flux density of $f_\mathrm{BH} < 1.9 \times 10^{-16}$ W m$^{-2}$ ${\mu}$m$^{-1}$ in the K band. In addition, we did not detect emission from the black hole in the hydrogen Br{\gamma} line. The non-detection of near-infrared emission from the black hole in Gaia BH3 indicates that its accretion of the giant star wind is presently occurring at most at a very low rate. This is consistent with the limit of $f_\mathrm{Edd} < 4.9 \times 10^{-7}$ derived previously on the Eddington ratio for an advection-dominated accretion flow. Deeper observations with GRAVITY may be able to detect the black hole as the companion star approaches periastron around 2030.

We study the response of mono-energetic stellar populations with initially isotropic kinematics to impulsive and adiabatic changes to an underlying dark matter potential. Half-light radii expand and velocity dispersions decrease as enclosed dark matter is removed. The details of this expansion and cooling depend on the time scale on which the underlying potential changes. In the adiabatic regime, the product of half-light radius and average velocity dispersion is conserved. We show that the stellar populations maintain centrally isotropic kinematics throughout their adiabatic evolution, and their densities can be approximated by a family of analytical radial profiles. Metallicity gradients within the galaxy flatten as dark matter is slowly removed. In the case of strong impulsive perturbations, stellar populations develop power-law-like density tails with radially biased kinematics. We show that the distribution of stellar binding energies within the dark matter halo substantially widens after an impulsive perturbation, no matter the sign of the perturbation. This allows initially energetically separated stellar populations to mix, to the extent that previously chemo-dynamically distinct populations may masquerade as a single population with large metallicity and energy spread. Finally, we show that in response to an impulsive perturbation, stellar populations that are deeply embedded in cored dark matter halos undergo a series of damped oscillations before reaching a virialised equilibrium state, driven by inefficient phase mixing in the harmonic potentials of cored halos. This slow return to equilibrium adds substantial systematic uncertainty to dynamical masses estimated from Jeans modeling or the virial theorem.

We use the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA) to detect $\rm CO(1-0)$, $\rm CO(3-2)$, and rest-frame 349-GHz continuum emission from an HI-selected galaxy, DLA1020+2733g, at $z\approx2.3568$ in the field of the $z=2.3553$ damped Lyman-$\alpha$ absorber (DLA) towards QSO J1020+2733. The VLA $\rm CO(1-0)$ detection yields a molecular gas mass of $(2.84\pm0.42)\times10^{11}\times(\alpha_{\rm CO}/4.36)\,\rm{M_\odot}$, the largest ever measured in an HI-selected galaxy. The DLA metallicity is $+0.28 \pm 0.16$, from the ZnII$\, \lambda2026$Å absorption line detected in a Keck Echellette Spectrograph and Imager spectrum. This continues the trend of high-metallicity DLAs being frequently associated with massive galaxies. We obtain a star-formation rate (SFR) of $\lesssim400~$M$_\odot~$yr$^{-1}$ from the rest-frame 349-GHz continuum emission, and a relatively long molecular gas depletion timescale of $\gtrsim0.6~$Gyr. The excitation of the J=3 rotational level is sub-thermal, with $r_{31}\equiv{L'_{\rm{CO(3-2)}}/L'_{\rm{CO(1-0)}}}=0.513\pm0.081$, suggesting that DLA1020+2733g has a low SFR surface density. The large velocity spread of the CO lines, $\approx500~\rm km~s^{-1}$, and the long molecular gas depletion timescale suggest that DLA1020+2733g is likely to be a cold rotating-disk galaxy.

Using data from the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES), we tentatively detect two molecular hydrogen (H$_2$) fluorescent emission features in high-redshift galaxies at 2.3 and 3.1$\sigma$. These features consist of many blended emission lines that result from the de-excitation cascade of H$_2$ molecules that have absorbed Lyman-Werner band photons. Our study targets galaxies at redshifts $z \geq 7$ galaxies of the early Universe, as they host some of the most extreme conditions in terms of star formation, molecular gas content, and the possible presence of outflows driven by starbursts and active galactic nuclei. To enhance the signal-to-noise ratio of H$_2$ emission features in the rest-frame wavelength range of 155-191 nm, we stack JWST/NIRSpec spectra from $z=7.0-13.4$. These stacked spectra also exhibit atomic emission features, such as CIV emission with a P-Cygni profile, as well as CIII], OIII], and HeII. The presence of these features and the slightly blue-shifted fluorescent H$_2$ lines suggest active multiphase atomic and molecular outflows may be common at these redshifts. Our results suggest that it is possible to search for FUV fluorescent H$_2$ lines in high redshift galaxies using JWST, which would enable the characterization of interstellar radiation fields, densities, and temperatures in extreme photodissociation environments.

Parker Solar Probe is a mission designed to explore properties of the solar wind closer than ever before. Detailed particle observations from the Solar Probe Cup (SPC) have primarily focused on examining the proton population in the solar wind. However, several periods throughout the Parker mission have indicated that SPC has observed a pronounced and distinctive population of fully ionized helium, He$^{2+}$. Minor ions are imprinted with properties of the solar wind's source region, as well as mechanisms active during outflow, making them sensitive markers of its origin and formation at the Sun. Through a detailed analysis of the He$^{2+}$ velocity distributions functions, this work examines periods where significant and persistent He$^{2+}$ peaks are observed with SPC. We compute the helium abundance and examine the stream's bulk speed, density, temperature, magnetic field topology, and electron strahl properties to identify distinctive solar wind features that can provide insight to their solar source. We find that nearly all periods exhibit an elevated mean helium composition ($8.34\%$) compared to typical solar wind and a majority ($\sim87\%$) of these periods are connected to coronal mass ejections, with the highest abundance reaching $23.1\%$. The helium abundance and number of events increases as the solar cycle approaches maximum with a weak dependence on speed. Additionally, the events not associated with a CME are clustered near the heliospheric current sheet suggesting they are connected to streamer belt outflows. However, there are currently no theoretical explanations that fully describe the range of depleted and elevated helium abundances observed.

We present the luminosity function and volumetric rate of a sample of Type IIP supernovae (SNe) from the Zwicky Transient Facility Census of the Local Universe survey (CLU). This is the largest sample of Type IIP SNe from a systematic volume-limited survey to-date. The final sample includes 330 Type IIP SNe and 36 low-luminosity Type II (LLIIP) SNe with $M_{\textrm{r,peak}}>-16$ mag, which triples the literature sample of LLIIP SNe. The fraction of LLIIP SNe is $19^{+3}_{-4}\%$ of the total CLU Type IIP SNe population ($8^{+1}_{-2}\%$ of all core-collapse SNe). This implies that while LLIIP SNe likely represent the fate of core-collapse SNe of $8-12$ \Msun\ progenitors, they alone cannot account for the fate of all massive stars in this mass range. To derive an absolute rate, we estimate the ZTF pipeline efficiency as a function of the apparent magnitude and the local surface brightness. We derive a volumetric rate of $(3.9_{-0.4}^{+0.4}) \times 10^{4}\ \textrm{Gpc}^{-3}\ \textrm{yr}^{-1}$ for Type IIP SNe and $(7.3_{-0.6}^{+0.6}) \times 10^{3}\ \textrm{Gpc}^{-3}\ \textrm{yr}^{-1}$ for LLIIP SNe. Now that the rate of LLIIP SNe is robustly derived, the unresolved discrepancy between core-collapse SN rates and star-formation rates cannot be explained by LLIIP SNe alone.

Faint $\gamma$-ray signatures emerge in Fermi-LAT data stacked scaled to the characteristic $\theta_{500}$ angles of MCXC galaxy clusters. After Paper I of this series thus discovered virial shocks, later supported in other bands, this second paper focuses on cluster cores. Stacking $1$-$100$ GeV source-masked data around clusters shows a significant ($4.7\sigma$ for 75 clusters) and extended central excess, inconsistent with central point sources. The resolved signal is best fit ($3.7\sigma$ TS-test) as hadronic emission from a cosmic-ray ion (CRI) distribution that is flat both spectrally ($p\equiv1-d\ln u/d\ln E=2.0\pm0.3$) and spatially (CRI-to-gas index $\sigma=0.1\pm0.4$), carrying an energy density $du(0.1\theta_{500})/d\ln E=10^{-13.6\pm0.5}$ erg cm$^{-3}$ at $E=100$ GeV energy; insufficient resolution would raise $p$ and $\sigma$. Such CRI match the long-predicted distribution needed to power diffuse intracluster radio emission in its various forms (mini-halos, giant halos, standard relics, their transitional forms, and mega-halos), disfavoring models that invoke electron (re)acceleration in weak shocks or turbulence. Stringent upper limits on residual $\gamma$-ray emission, e.g. from dark-matter annihilation, are imposed.

In the coming years, surveys such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) are expected to increase the number of observed Tidal Disruption Events (TDEs) substantially. We employ Monte Carlo integration to calculate the unlensed and lensed TDE rate as a function of limiting magnitude in $u$, $g$, $r$, and $i$-bands. We investigate the impact of multiple luminosity models, black hole mass functions (BHMFs), and flare temperatures on the TDE rate. Notably, this includes a semi-analytical model, which enables the determination of the TDE temperature in terms of black hole (BH) mass. We predict the highest unlensed TDE rate to be in $g$-band. It ranges from $16$ to $5,440\;\mathrm{yr}^{-1}\;(20,000\;\mathrm{deg}^2)^{-1}$ for the Zwicky Transient Facility, being more consistent with the observed rate at the low end. For LSST, we expect a rate in $g$-band between $3,580$ and $82,060\;\mathrm{yr}^{-1}\;(20,000\;\mathrm{deg}^2)^{-1}$. A higher theoretical prediction is understandable, as we do not consider observational effects such as completeness. The unlensed and lensed TDE rates are insensitive to the redshift evolution of the BHMF, even for LSST limiting magnitudes. The best band for detecting lensed TDEs is also $g$-band. Its predicted rates range from $0.43$ to $15\;\mathrm{yr}^{-1}\;(20,000\;\mathrm{deg}^2)^{-1}$ for LSST. The scatter of predicted rates reduces when we consider the fraction of lensed TDEs; that is, a few in ten thousand TDEs will be lensed. Despite the large scatter in the rates of lensed TDEs, our comprehensive considerations of multiple models suggest that lensed TDEs will occur in the $10$-year LSST lifetime, providing an exciting prospect for detecting such events. We expect the median redshift of a lensed TDE to be between $1.5$ and $2$. In this paper, we additionally report on lensed TDE properties, such as the BH mass and time delays.

White dwarf stars are the most common final stage of stellar evolution. Since the serendipitous discovery of the first white dwarf by William Herschel and the first physical models by Subrahmanyan Chandrasekhar and Arthur Eddington, there have been a lot of advances in the field fueled by new observational data. With new astrometric measurements enabling us to identify hundreds of thousands of white dwarf candidates, and spectroscopic surveys revealing a plethora of chemical elements in white dwarf atmospheres pointing at spectral evolution and interaction with planetary bodies, there is a lot we can learn from the characterization of observed white dwarfs. Here we provide an observational overview of white dwarf stars, describing how they are identified and characterized, and the main properties of the observed population.

Massive quiescent galaxies at high redshift are significantly more compact than their present-day counterparts. We investigate the roles, in determining this evolution, of major and minor mergers, and of the accretion of diffuse envelopes of stars and dark matter. We model the evolution in stellar mass ($M_\star$), effective radius ($R_{\rm e}$) and effective stellar velocity dispersion ($\sigma_{\rm e}$) of a representative massive quiescent galaxy from $z=3$ to $z=0$, and compare the model with the observed redshift-dependent $R_{\rm e}$-$M_\star$ and $\sigma_{\rm e}$-$M_\star$ relations. In the model we account for the effects of collisionless (dry) major (satellite-to-main galaxy mass ratio $\xi>1/4$) and minor ($1/10<\xi<1/4$) mergers, using analytic recipes consistent with the results of N-body simulations of binary mergers. For the poorly constrained mini mergers ($\xi<1/10$) we explore both a "standard" model (based on the same assumptions used in the case of higher-$\xi$ mergers), and a heuristic "envelope accretion" model, aimed at describing the case in which diffuse satellites are completely disrupted in the galaxy outskirts. Major and minor dry mergers, at rates estimated observationally from galaxy-pair counts, induce relatively small variations of $R_{\rm e}$ and $\sigma_{\rm e}$, accounting only for $\approx 6\%$ of the size evolution and $\approx 40\%$ of the velocity-dispersion evolution observed from $z=3$ to $z=0$. As addition to major and minor dry mergers, envelope accretion performs better than standard mini mergers at reproducing the redshift-dependence of the $R_{\rm e}$-$M_\star$ and $\sigma_{\rm e}$-$M_\star$ relations, being also consistent with plausible evolutionary scenarios of scaling relations involving the mass of the central supermassive black hole.

We investigate the impact of neutrino emission via Hawking radiation from primordial black holes (PBHs) on the cosmological effective number of neutrino species, $N_{\mathrm{eff}}$, after neutrino decoupling. By comparing this effect with observational limits, we derive bounds on the abundance of light PBHs. Our analysis incorporates two previously unaccounted-for effects: the emission of secondary neutrinos from unstable particles, which increases $N_{\mathrm{eff}}$, and the modification of the neutrino-photon temperature ratio due to particle emission heating the photon plasma, which lowers $N_{\mathrm{eff}}$. Overall, including these effects allows us to impose constraints on PBH masses with initial masses in the range $10^9~{\rm g}\lesssim M_{\rm ini} \lesssim 10^{13}~{\rm g}$. However, our limits remain less stringent than those derived from Big Bang Nucleosynthesis.

We build upon the state-of-the-art semi-analytic model \texttt{FEGA24} (Formation and Evolution of GAlaxies, \citealt{contini2024d}), which integrates the latest prescriptions relevant to galaxy formation and evolution, alongside a comprehensive AGN feedback model. This model incorporates three modes of feedback: negative (preventing excessive cooling), positive (enhancing star formation), and hot gas ejection (expelling gas beyond the virial radius of halos). These modes operate in a coordinated manner: the negative mode regulates the cooling process, the positive mode promotes bursts of star formation, and the hot gas ejection mode expels gas beyond the virial radius when the AGN is sufficiently powerful. Our updated semi-analytic model, \texttt{FEGA25}, retains the qualitative and quantitative consistency of the analyses presented in \cite{contini2024d}, while delivering more robust results. Notably, \texttt{FEGA25} provides a more detailed characterization of the fraction of red galaxies as a function of stellar mass, predicts a main sequence of star-forming galaxies more consistent with observations, and estimates the fraction of hot gas in halos closer to observed values. These findings underscore the importance of a physical mechanism capable of ejecting hot gas beyond the virialized region of dark matter halos without significantly altering the stellar and cold gas components. Such a mechanism is crucial to ensure the proper functioning of other processes, such as cooling and star formation. Since supernova feedback is already modeled at its maximum efficiency, AGN feedback emerges as the natural candidate for this role.

We present the discovery of Andromeda XXXV, the faintest Andromeda satellite galaxy discovered to date, identified as an overdensity of stars in the Pan-Andromeda Archaeological Survey and confirmed via Hubble Space Telescope imaging. Located at a heliocentric distance of $927^{+76}_{-63}$ kpc and $158^{+57}_{-45}$ kpc from Andromeda, Andromeda XXXV is an extended ($r_h = 53\,^{+13}_{-11}$ pc), elliptical ($\epsilon = 0.4\, \pm 0.2$), metal-poor ($[\text{Fe}/\text{H}]\sim-1.9$) system, and the least luminous ($M_V=-5.2 \pm 0.3$) of Andromeda's dwarf satellites discovered so far. Andromeda XXXV's properties are consistent with the known population of dwarf galaxies around the Local Group, bearing close structural resemblance to the Canes Venatici II and Hydra II Milky Way (MW) dwarf satellite galaxies. Its stellar population, characterized by a red horizontal branch or a red clump feature, mirrors that of other Andromeda satellite galaxies in showing evidence for a spread in age and metallicity, with no signs of younger stellar generations. This age-metallicity spread is not observed in MW satellites of comparable stellar mass, and highlights the persistent differences between the satellite systems of Andromeda and the MW, extending even into the ultrafaint regime.

We present an analysis of the cluster X-ray morphology and active galactic nucleus (AGN) activity in nine $z\sim1$ galaxy clusters from the Massive and Distant Clusters of $WISE$ Survey (MaDCoWS) observed with $Chandra$. Using photon asymmetry ($A_{\text{phot}}$) to quantify X-ray morphologies, we find evidence that the four most dynamically disturbed clusters are likely to be mergers. Employing a luminosity cut of $7.6\times10^{42}$ erg/s to identify AGN in the 0.7-7.0 keV, we show that the majority of these clusters host excess AGN compared to the local field. We use the cumulative number-count ($\log N-\log S$) model to predict AGN incidence in cluster isophotes under this luminosity cut. Our analysis finds evidence (at $> 2\sigma$) of a positive correlation between AGN surface densities and photon asymmetry, suggesting that a disturbed cluster environment plays a pivotal role in regulating AGN triggering. Studying AGN incidence in cluster X-ray isophotes equivalent in area to $1.0r_{500}$, we find that the AGN space density inversely scales with cluster mass as $\sim M^{-0.5^{+0.18}_{-0.18}}$ at the 3.18$\sigma$ level. Finally, when we separately explore the cluster mass dependence of excess AGN surface density in disturbed and relaxed clusters, we see tentative evidence that the two morphologically distinct sub-populations exhibit diverging trends, especially near the outskirts, likely due to cluster merger-driven AGN triggering/suppression.

Current and future surveys rely on machine learning classification to obtain large and complete samples of transients. Many of these algorithms are restricted by training samples that contain a limited number of spectroscopically confirmed events. Here, we present the first real-time application of Active Learning to optimise spectroscopic follow-up with the goal of improving training sets of early type Ia supernovae (SNe Ia) classifiers. Using a photometric classifier for early SN Ia, we apply an Active Learning strategy for follow-up optimisation using the real-time FINK broker processing of the ZTF public stream. We perform follow-up observations at the ANU 2.3m telescope in Australia and obtain 92 spectroscopic classified events that are incorporated in our training set. We show that our follow-up strategy yields a training set that, with 25% less spectra, improves classification metrics when compared to publicly reported spectra. Our strategy selects in average fainter events and, not only supernovae types, but also microlensing events and flaring stars which are usually not incorporated on training sets. Our results confirm the effectiveness of active learning strategies to construct optimal training samples for astronomical classifiers. With the Rubin Observatory LSST soon online, we propose improvements to obtain earlier candidates and optimise follow-up. This work paves the way to the deployment of real-time AL follow-up strategies in the era of large surveys.

Since the recent discovery of the directly imaged super-Jovian planet AF Lep b, several studies have been conducted to characterize its atmosphere and constrain its orbital parameters. AF Lep b has a measured dynamical mass of $3.68 \pm 0.48$ MJup, a radius of $1.3 \pm 0.15$ RJup, a nearly circular orbit in spin-orbit alignment with the host star, a relatively high metallicity, and a near-solar to super-solar C/O ratio. However, key parameters such as the rotational velocity and radial velocity could not be estimated as they require high-resolution spectroscopic data that is impossible to obtain with classical spectrographs. AF Lep b was recently observed with the new HiRISE visitor instrument at the VLT, with the goal of obtaining high-resolution (R~140,000) spectroscopic observations to better constrain the orbital and atmospheric parameters of the young giant exoplanet. We compare the extracted spectrum of AF Lep b to self-consistent atmospheric models using ForMoSA. We then use our measurements of the radial velocity of the planet to provide new constraints on the orbit of the planet. From the forward modeling, we find a C/O ratio that aligns with previous low-resolution analyses, and we confirm the super-solar metallicity. We also confirm unambiguously the presence of methane in the atmosphere of the companion. Based on all available relative astrometry and radial velocity measurements of the host star, we show that two distinct orbital populations are possible for the companion. We derive the radial velocity of AF Lep b to be $10.51 \pm 1.03$ km/s, and show that this value agrees well with one of the two orbital solutions, allowing us to rule out an entire family of orbits. Additionally, assuming that the rotation and orbit are coplanar, the derived planet's rotation rate is consistent with the observed trend of increasing spin velocity with higher planet mass.

Stellar contamination has long been recognized as a major bottleneck in transmission spectroscopy, limiting our ability to accurately characterize exoplanet atmospheres, particularly for terrestrial worlds. In response, significant observational efforts have shifted toward emission spectroscopy as a potentially more robust alternative, exemplified by initiatives such as the 50 hour JWST "Rocky Worlds" Director's Discretionary Time (DDT) program. However, the extent to which emission spectroscopy may be affected by stellar effects remain mostly unexplored, in stark contrast with exploration and mitigation work for transmission spectroscopy. In this study, we assess the impact of imperfect knowledge of stellar spectra on exoplanet atmospheric retrievals from emission spectroscopy. For TRAPPIST-1, a 29 ppm precision on the 15 um eclipse depth is required to differentiate surface albedo with precision of 0.1 at 3 sigma confidence. Yet, we find that current 15 um eclipse depth estimations using different stellar models introduce a 40 to 80 ppm uncertainty. Moreover, such model precision should not be mistaken for a model accuracy, as we also find that JWST MIRI photometric observations of the star in-eclipse at 12.8 and 15 um, acquired with a 90 ppm 1 sigma precision, fall outside the range defined by existing models by 7 to 21 %. This lack of precision and accuracy leads to a degeneracy in retrievals of planetary albedo and impacts the search for atmospheres, effects that could be solved by acquiring the true mid-IR spectrum of the star. We therefore recommend a revised observation design for JWST secondary eclipse measurements in which mid-IR stellar spectra should be systematically acquired alongside eclipse data.

The rotational distribution of asteroids as a function of their size is used {as a diagnostic of} their physical properties and evolution. Recent photometric surveys from the Gaia mission, allowing observation of asteroids with long spin periods (for example $\geq 24$h), found an excessive group of slow rotators and a gap separating them from faster rotators, which is unexplained by current theories. Here we developed an asteroid rotational evolution model capable of reproducing the observed distribution. {We suggest that this distribution is regulated by the competition between collisions and internal friction dampening of "tumblers" -asteroids with unstable rotation vectors, and that the slow rotator group is mainly populated by tumblers.} {We constrain the product of the rigidity and quality factor, which relates to the body's viscosity, to $\mu Q \sim 4 \times 10^9~$Pa. This number, two orders of magnitude smaller than the one assumed for monolithic boulders,} implies that {rubble pile} asteroids could have a porous structure or a thick regolith layer, and undergo stronger tidal effects.

The origin of the slow solar wind is not well understood, unlike the fast solar wind which originates from coronal holes. In-situ elemental abundances of the slow solar wind suggest that it originates from initially closed field lines that become open. Coronal hole boundary regions are a potential source of slow solar wind as there open field lines interact with the closed loops through interchange reconnection. Our primary aim is to quantify the role of interchange reconnection at the boundaries of coronal holes. To this end, we have measured the relative abundances of different elements at these boundaries. Reconnection is expected to modulate the relative abundances through the first ionization potential (FIP) effect. For our analysis we used spectroscopic data from the extreme ultraviolet imaging spectrometer (EIS) on board Hinode. To account for the temperature structure of the observed region we computed the differential emission measure (DEM). Using the DEM we were able to infer the ratio between coronal and photospheric abundances, known as the FIP bias. By examining the variation of the FIP bias moving from the coronal hole to the quiet Sun, we have been able to constrain models of interchange reconnection. The FIP bias variation in the boundary region around the coronal hole has an approximate width of 30-50 Mm, comparable to the size of supergranules. This boundary region is also a source of open flux into interplanetary space. We find that there is an additional ~30% open flux that originates from this boundary region.

During thermonuclear bursts, it is suspected that {\bf the cooling of the corona by the burst emission} may be the cause of hard X-ray {\bf deficits}. Although this {\bf deficit} has been observed in nine sources, it has not been observed {\bf from} 4U~1608--52, a nearby prolific burster. Therefore, the authenticity and universality of the hard X-ray {\bf deficit} may be in question. To investigate this suspicion, Insight-HXMT performed cadence observations during the low/hard state of 4U~1608--52 in September 2022 and detected 10 thermonuclear X-ray bursts. Two of these bursts show a double-peaked structure in the soft X-ray band, which could be caused by the high temperature of the burst emission and a marginal photospheric radius expansion (PRE) around the burst peak time. This is indicated by their peak fluxes being up to the Eddington limit and having a large color factor at the peak of the bursts. The hard X-ray deficit is significantly observed during bursts at $>$ 30 keV. Furthermore, the fraction of this deficit shows saturation at 50\% for the first 8 bursts. This saturation may indicate that the corona is layered and only a part of the corona is cooled by the bursts. For example, the part close to the NS surface is cooled while the rest remains intact during bursts. This result provides a clue to the geometry of the corona, e.g., a possible scenario is that the corona has two forms: a quasi-spheric corona between the NS and the disk, and a disk-corona on both surfaces of the disk.

IceCube has observed a diffuse astrophysical neutrino flux over the energy region from a few TeV to a few PeV. At PeV energies, the spectral shape is not yet well measured due to the low statistics of the data. This analysis probes the gap between 1 PeV and 10 PeV by using high-energy downgoing muon neutrinos. To reject the large atmospheric muon background, two complementary techniques are combined. The first technique selects events with high stochasticity to reject atmospheric muon bundles whose stochastic energy losses are smoothed due to high muon multiplicity. The second technique vetoes atmospheric muons with the IceTop surface array. Using 9 years of data, we found two neutrino candidate events in the signal region, consistent with expectation from background, each with relatively high signal probabilities. A joint maximum likelihood estimation is performed using this sample and an independent 9.5-year sample of tracks to measure the neutrino spectrum. A likelihood ratio test is done to compare the single power-law (SPL) vs. SPL+cutoff hypothesis; the SPL+cutoff model is not significantly better than the SPL. High-energy astrophysical objects from four source catalogs are also checked around the direction of the two events. No significant coincidence was found.

Solar flares are regularly observed in extreme ultraviolet (EUV), soft X-rays (SXR), and hard X-rays (HXR). However, those in near and mid-UV are sparse. The Solar Ultraviolet Imaging Telescope (SUIT) onboard the Aditya-L1, launched on 2nd September, 2023 provides regular observations in the 200-400 nm wavelength range through eleven filters. Here, we report the observation of the X6.3 flare on Feb 22, 2024 using eight narrow band (NB) filters of SUIT. We have also used co-spatiotemporal observations from SDO/AIA, Solar Orbiter/STIX, GONG H$\alpha$, Aditya-L1/SoLEXS and GOES. We obtained light curves over the flaring region from AIA 1600, 1700 Å and GONG H$\alpha$ and compared them with the disk-integrated lightcurve obtained from GOES and SoLEXS SXR and STIX HXR. We find that the flare peaks in SUIT NB01, NB03, NB04, and NB08 filters simultaneously with HXR, 1600, and 1700 Å along with the peak temperature obtained from SoLEXS. In contrast, in NB02 and NB05, the flare peaks $\sim$ 2 minutes later than the HXR peak, while in NB06 and NB07, the flare peaks $\sim$ 3 minutes after the GOES soft X-ray peak. To the best of our knowledge, this is the first observation of a flare in these wavelengths (except in NB03, NB04 and NB05). Moreover, for the first time, we show the presence of a bright kernel in NB02. These results demonstrate the capabilities of SUIT observations in flare studies.

The sensitivity of axion dark matter searches depends on the signal window that results from the velocity dispersion of axion dark matter. Since the ratio of signal windows is about 6500 between the standard halo and the big flow axion dark matter, each axion dark matter search usually uses a separate data acquisition (DAQ) channel with a different frequency resolution bandwidth (RBW). In this work, we demonstrate axion dark matter searches covering the standard halo, the tidal stream, and the big flow employing a DAQ channel starting with a single high resolution RBW, without sacrificing the DAQ efficiency, where the DAQ process includes online fast Fourier transforms and writing the outputs to disk. Assuming the total amount of data is sensitive to Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) axion dark matter that follows the standard halo model and makes up 100\% of the local dark matter density, the same data can also be used for the tidal stream and the big flow axion dark matter searches that would be sensitive to DFSZ axion dark matter that constitute 19.2\% and 12.4\% of the local dark matter densities, respectively, at a 90\% confidence level. We also report that the filtering of the individual power spectra acquired with a relatively high resolution RBW e.g., for the big flow search can prevent a possible significant degradation in the signal to noise ratio from the searches in the lower resolution RBW's, i.e., the standard halo and tidal stream searches.

This work proposes a scalable probabilistic latent variable model based on Gaussian processes (Lawrence, 2004) in the context of multiple observation spaces. We focus on an application in astrophysics where data sets typically contain both observed spectral features and scientific properties of astrophysical objects such as galaxies or exoplanets. In our application, we study the spectra of very luminous galaxies known as quasars, along with their properties, such as the mass of their central supermassive black hole, accretion rate, and luminosity-resulting in multiple observation spaces. A single data point is then characterized by different classes of observations, each with different likelihoods. Our proposed model extends the baseline stochastic variational Gaussian process latent variable model (GPLVM) introduced by Lalchand et al. (2022) to this setting, proposing a seamless generative model where the quasar spectra and scientific labels can be generated simultaneously using a shared latent space as input to different sets of Gaussian process decoders, one for each observation space. Additionally, this framework enables training in a missing data setting where a large number of dimensions per data point may be unknown or unobserved. We demonstrate high-fidelity reconstructions of the spectra and scientific labels during test-time inference and briefly discuss the scientific interpretations of the results, along with the significance of such a generative model.

I will demonstrate the effectiveness of Physics-Informed Neural Networks (PINNs) in solving partial differential equations (PDEs) when training data are scarce or noisy. The training data can be located either at the boundaries or within the domain. Additionally, PINNs can be used as an inverse method to determine unknown coefficients in the equations. This study will highlight the application of PINNs in modeling magnetohydrodynamic processes relevant to strongly magnetized plasmas, such as those found in the solar corona.

The dominant fraction of the extragalactic $\gamma$-ray sources are blazars, active galactic nuclei with jets inclined at a small angle to the line of sight. Apart from blazars, a few dozen narrow-line Seyfert 1 galaxies (NLS1) and a number of radio galaxies are associated with $\gamma$-ray sources. The identification of $\gamma$-ray sources requires multiwavelength follow-up observations since several candidates could reside within the relatively large $\gamma$-ray localisation area. The $\gamma$-ray source 4FGL 0959.6+4606 was originally associated with a radio galaxy. However, follow-up multiwavelength work suggested a nearby NLS1 as the more probable origin of the $\gamma$-ray emission. We performed high-resolution very long baseline interferometry (VLBI) observation at 5 GHz of both proposed counterparts of 4FGL 0959.6+4606. We clearly detected the NLS1 source SDSS J095909.51+460014.3 with relativistically boosted jet emission. On the other hand, we did not detect milliarcsecond-scale compact emission in the radio galaxy 2MASX J09591976+4603515. Our VLBI imaging results suggest that the NLS1 object is the origin of the $\gamma$-ray emission in 4FGL 0959.6+4606.

Star-forming galaxies (SFGs) are considered to be an important component of the diffuse extragalactic gamma-ray background (EGB) radiation observed in 0.1 -- 820 GeV, but their quantitative contribution has not yet been precisely determined. In this study, we aim to provide the currently most reliable estimate of the contribution of SFGs based on careful calibration with gamma-ray luminosities of nearby galaxies and physical quantities (star formation rate, stellar mass, and size) of galaxies observed by high-redshift galaxy surveys. Our calculations are based on the latest database of particle collision cross-sections and energy spectra of secondary particles, and take into account not only hadronic but also leptonic processes with various radiation fields in a galaxy. We find that SFGs are not the dominant component of the unresolved EGB measured by Fermi; the largest contribution is around 50% -- 60% in the 1 -- 10 GeV region, and the contribution falls rapidly in lower and higher energy ranges. This result appears to contradict a previous study, which claimed that SFGs are the dominant component of the unresolved EGB, and the origin of the discrepancy is examined. In calculations of cosmic-ray production, propagation, and interaction in a galaxy, we try models developed by two independent groups and find that they have little impact on EGB.

Context. When leaving the main sequence (MS) for the red-giant branch (RGB), subgiant stars undergo fast structural changes. Consequently, their observed oscillation spectra mirror these changes, constituting key tracers of stellar structure and evolution. However, the complexity of their spectra makes their modelling an arduous task, which few authors have undertaken. Gemma (KIC11026764) is a young subgiant with $45$ precise oscillation modes observed with Kepler, making it the ideal benchmark for seismic modelling. Aims. This study is aimed at modelling the subgiant Gemma, taking advantage of most of the precise seismic information available. This approach enables us to pave the way for the seismic modelling of evolved solar-like stars and provide the relevant insights into their structural evolution. Methods. Using our Levenberg-Marquardt stellar modelling tool, we built a family of models representative of Gemma's measured seismic indicators obtained via our seismic tool, EGGMiMoSA. We studied the structural information these indicators hold by carefully varying stellar parameters. We also complemented the characterisation with information held by \who indicators and non-seismic data. Results. From the extensive set of models we built and using most of the seismic information at hand, including two $\ell=1$ and one $\ell=2$ mixed modes, we were able to probe the chemical transition at the hydrogen-burning shell. Indeed, we have demonstrated that among our models, only the ones with the sharpest chemical gradient are able to reproduce all the seismic information considered. One possibility to account for such a gradient is the inclusion of a significant amount of overshooting, namely $\alpha_{\textrm{ov}}=0.17$, which is unexpected for low-mass stars such as Gemma (expected mass of about $1.15~M_{\odot}$).

Machine learning is a useful tool for identifying radio pulses from cosmic-ray air showers and for cleaning such pulses from radio background. This can lower the detection threshold and increase the accuracy for the pulse time and amplitude. We have trained Convolutional Neural Networks (CNNs) using CoREAS simulations and background recorded by a prototype station at the IceTop surface array at the South Pole and have applied them to air-shower measurements by this station. The station consists of 3 SKALA antennas and 8 scintillators, which are used to trigger the readout of the antennas upon a sixfold coincidence. Afterwards, the radio signal is filtered to the band of 70-350 MHz. By applying neural networks to search for radio signals in about four months of data, we find about five times more events than by a traditional method based on a signal-to-noise ratio cut after filtering for radio frequency interferences. Despite the lower threshold, the purity of the selected events seems to improve, and the angular resolution of the radio measurements does not deteriorate, which we have confirmed by a comparison of the reconstructed shower direction with IceTop. This analysis thus provides experimental confirmation that neural networks can indeed be used to clean air-shower radio signals from background and to lower the radio detection threshold of hybrid arrays combining particle and radio detectors.

With a mass, radius, and mean density similar to Earth's, the rocky planet GJ 1132 b is the first truly small planet for which an atmosphere detection was proposed. If confirmed, ultra-reduced magma outgassing is the only mechanism capable of producing HCN and H$_2$O in large enough quantities to match the HST observations. The proposed atmosphere detection, however was challenged by reanalysis of the same HST data by different teams. Recent JWST observations returned ambiguous results due to the unaccounted for variability seen between two different visits. Here we report the analysis of three CRIRES+ transit observations of GJ 1132 b in order to determine the presence or absence of He I, HCN, CH$_4$, and H$_2$O in its atmosphere. We are unable to detect the presence of any of these species in the atmosphere of GJ 1132 b assuming a clear, H$_2$-dominated atmosphere, although we can place upper limits for the volume mixing ratios of CH$_4$, HCN, and H$_2$O using injections tests and atmospheric retrievals. These retrieved upper limits show the capability of CRIRES+ to detecting chemical species in rocky exoplanets, if the atmosphere is H$_2$ dominated. The detection of the atmospheres of small planets with high mean molecular weight, and the capability to distinguish between the variability introduced by stellar activity and/or the planetary atmosphere will require high-resolution spectrographs in the upcoming extremely large telescopes.

We use the Five-hundred-meter Aperture Spherical radio Telescope (FAST) to observe the bright millisecond pulsar (MSP) PSR B1937+21 (J1939+2134) and record the data in the band from 1.0 GHz to 1.5 GHz. We measure the neutral hydrogen (HI) emission and absorption lines near 1420 MHz ($\lambda \simeq 21$ cm). We derive the kinematic distance of the pulsar with the HI observation, and update the upper bound of kinematic distance from the previous $14.8\pm 0.9$ npc in the Outer Arm to the nearer $9.4\pm 0.5$ kpc in the Perseus Arm. By comparing with the archival absorption spectra observed decades ago, we notice possible variations in the absorption spectra towards this pulsar, which corresponds to a possible tiny-scale atomic structure (TSAS) of a few AU in size. We also verify the apparent faster-than-light anomalous dispersion at the HI absorption line of this pulsar previously reported.

IceCube-Gen2, the next generation extension of the IceCube Neutrino Observatory at the South Pole, offers a unique scientific potential for cosmic-ray physics at PeV to EeV energies complementing the main science case of neutrino astronomy. The cosmic-ray science case will be enabled by a surface array on top of an extended optical array deep in the polar ice. The optical array measures TeV muons of air showers, and the surface array primarily measures the electromagnetic shower component and low-energy muons. The design of the surface array foresees scintillation panels providing a full-efficiency threshold for near-vertical proton showers of 0.5 PeV and radio antennas increasing the measurement accuracy for the electromagnetic shower component in the energy range of the Galactic-to-extragalactic transition. Compared to IceCube, the aperture for air showers measured in coincidence with the surface and optical arrays will increase by a factor of 30, due to the larger area and angular acceptance in zenith angle. The science potential includes both, the particle physics of air showers, such as prompt muons, and the astrophysics of the highest energy Galactic cosmic-rays, enabled by the higher sensitivity for the mass composition and anisotropy of cosmic rays, and by the search for PeV photons. This proceeding summarizes the science case and design of the surface array as presented in the recently released IceCube-Gen2 Technical Design Report: this https URL

The Enhanced Resolution Imager and Spectrograph (ERIS) is the new Adaptive-Optics (AO) assisted Infrared instrument at the Very Large Telescope (VLT). Its refurbished Integral Field Spectrograph (IFS) SPIFFIER leverages a new AO module, enabling high-contrast imaging applications and giving access to the orbital and atmospheric characterisation of super-Jovian exoplanets. We test the detection limits of ERIS and demonstrate its scientific potential by exploring the atmospheric composition of the young super-Jovian AF Lep b and improving its orbital solution by measuring its radial velocity relative to its host star. We present new spectroscopic observations of AF Lep b in $K$-band at $R\sim 11000$ obtained with ERIS/SPIFFIER at the VLT. We reduce the data using the standard pipeline together with a custom wavelength calibration routine, and remove the stellar PSF using principal component analysis along the spectral axis. We compute molecular maps by cross-correlating the residuals with molecular spectral templates and measure the radial velocity of the planet relative to the star. Furthermore, we compute contrast grids for molecular mapping by injecting fake planets. We detect a strong signal from H$_{2}$O and CO but not from CH$_{4}$ or CO$_{2}$. This result corroborates the hypothesis of chemical disequilibrium in the atmosphere of AF Lep b. Our measurement of the RV of the planet yields $\Delta v_{\mathrm{R,P\star}} = 7.8 \pm 1.7$ km s$^{-1}$. This enables us to disentangle the degeneracy of the orbital solution, namely the correct longitude of the ascending node is $\Omega=248^{+0.4}_{-0.7}$ deg and the argument of periapsis is $\omega=109^{+13}_{-21}$ deg. Our results demonstrate the competitiveness of the new ERIS/SPIFFIER instrument for the orbital and atmospheric characterisation of exoplanets at high contrast and small angular separation.

Radio antennas have become a standard tool for the detection of cosmic-ray air showers in the energy range above $10^{16}\,$eV. The radio signal of these air showers is generated mostly due to the deflection of electrons and positrons in the geomagnetic field, and contains information about the energy and the depth of the maximum of the air showers. Unlike the traditional air-Cherenkov and air-fluorescence techniques for the electromagnetic shower component, radio detection is not restricted to clear nights, and recent experiments have demonstrated that the measurement accuracy can compete with these traditional techniques. Numerous particle detector arrays for air showers have thus been or will be complemented by radio antennas. In particular when combined with muon detectors, the complementary information provided by the radio antennas can enhance the total accuracy for the arrival direction, energy and mass of the primary cosmic rays. Digitization and computational techniques have been crucial for this recent progress, and radio detection will play an important role in next-generation experiments for ultra-high-energy cosmic rays. Moreover, stand-alone radio experiments are under development and will search for ultra-high-energy photons and neutrinos in addition to cosmic rays. This article provides a brief introduction to the physics of the radio emission of air showers, an overview of air-shower observatories using radio antennas, and highlights some of their recent results.

We present $I$-band time-series photometric variability studies of three known nearby ($\sim$ 140 pc) and young ( $\sim$ 1 Myr) brown dwarfs (BD) in the Taurus star-forming region in the Perseus Molecular Cloud. From 10 nights of observations over a time span of 10 years, with a typical run of 3 to 6 hours each night, we estimated that the BDs show unstable short-scale periodicity from 1.5 to 4.8 hours. Using the long-term photometry from the Transiting Exoplanet Survey Satellite (TESS), we have conducted a time-resolved variability analysis of CFHT-BD-Tau 3 and CFHT-BD-Tau 4, revealing orbital periods of $\sim$ 0.96 days and $\sim$ 3 days respectively, consistent with earlier studies. We also found two superflares in TESS sector 43 data for CFHT-BD-Tau 4 and estimated the flare energies as $7.09\times10^{35}$ erg and $3.75\times10^{36}$ erg. A magnetic field of $\sim3.39 ~kG$ is required to generate such flare energies on this BD. We performed spot modelling analysis on CFHT-BD-Tau 3 and CFHT-BD-Tau 4 to address the variability detected in the data using the package BASSMAN. Spectral energy distribution and infrared colours of the sources suggest that they have a sufficient amount of circumstellar material around them.

A glitch is a rare and sudden increase in the otherwise steadily decreasing rotation rate of a pulsar. Its cause is widely attributed to the transfer of angular momentum to the crust of the star from the array of superfluid vortices enclosed within. The magnitude of such an increase defines the size of the glitch. The distribution of glitch sizes in individual pulsars, the power-law being the most argued for, is shrouded in uncertainty due to the small sample size. From a Bayesian perspective, we revisit the data for PSR J0537-6910, the pulsar with the most glitches, and find a bimodality in the distribution, reminiscent of the Bactrian camel. To understand this bimodality, we use a superfluid vortex simulator and study three independent neutron star paradigms: (i) Annular variation in pinning strength to account for the predicted differences between the crust and the core; (ii) Sectorial triggers to mimic local disturbances; and (iii) Stress-waves to model global disturbances. We find that annular variation in pinning introduces a bimodality in the glitch-size distribution and that sectorial triggers do so weakly. Stress-waves do not lead to any such features for the range of parameters tested. This provides us with new insights into the effects of various perturbations on the vortex dynamics and the glitch statistics of neutron stars.

The characteristic timescale ($\tau$) of continuum variability of the accretion disk in active galactic nuclei (AGNs) is known to be related to the thermal timescale, which is predicted to scale with AGN luminosity ($L$) and restframe wavelength ($\lambda_{\rm RF}$) as $t_{\rm th} \propto L^{0.5} \lambda_{\rm RF}^2$ in the standard disk model. Using multi-epoch spectroscopic data from the Sloan Digital Sky Survey Reverberation Mapping project, we construct ultraviolet ensemble structure functions of luminous AGNs as a function of their luminosity and wavelength. Assuming that AGNs exhibit a single universal structure function when $\Delta t$ is normalized by $\tau$, wherein $\tau \propto L^{a} \lambda_{\rm RF}^{b}$, we find $a=0.50\pm0.03$ and $b=1.42\pm0.09$. While the value of $a$ aligns with the prediction from the standard disk model, $b$ is significantly smaller than expected, suggesting that the radial temperature (color) profile of the accretion disk is significantly steeper (shallower) than the standard disk model. Notably, this discrepancy with theory has been observed in previous studies based on spectroscopic reverberation mapping and gravitational microlensing. Although no current model of accretion disks fully matches our results, our findings provide valuable constraints for testing future physical models.

We present the timing and spectral analyses of the NICER, NuSTAR, and IXPE observations of the magnetar 1E 1841-045 covering 82 days following its August 2024 bursting activity as well as radio observations utilizing MeerKAT and Effelsberg. We supplement our study with a historical NuSTAR and all 2024 pre-outburst NICER observations. The outburst is marked by an X-ray flux enhancement of a factor 1.6 compared to the historical level, predominantly driven by a newly-formed non-thermal emitting component with a photon index $\Gamma=1.5$. This flux showed a 20% decay at the end of our monitoring campaign. The radio monitoring did not reveal any pulsed radio emission with an upper-limit of 20 mJy and 50 mJy ms on the mean flux density and single pulse fluence, respectively. We detect a spin-up glitch at outburst onset with a $\Delta\nu=6.1\times10^{-8}$ Hz and a $\Delta\dot{\nu}=-1.4\times10^{-14}$ Hz s$^{-1}$, consistent with the near-universality of this behavior among the continuously-monitored magnetars. Most intriguingly, the 1E 1841-045 2-10 keV pulse profile is markedly different compared to pre-outburst; it shows a new, narrow (0.1 cycles) peak that appears to shift towards merging with the main, persistently-present, pulse. This is the second case of pulse-peak migration observed in magnetars after SGR 1830$-$0645, and the two sources exhibit a similar rate of phase shift. This implies that this phenomenon is not unique and might present itself in the broader population. The newly-formed peak for 1E 1841-045 is non-thermal, with emission extending to $\gtrsim20$ keV, in contrast to the case of SGR 1830$-$0645. Our results are consistent with an untwisting magnetic field bundle with migration towards the magnetic pole, perhaps accompanied by plastic motion of the crust.

We report the spectroscopic confirmation of the background source of the most distant Einstein ring known to date, the COSMOS-Web ring. This system consists of a complete Einstein ring at $z=5.1$, lensed by a massive early-type galaxy at $z\sim2$. The redshift $z=5.1043\pm0.0004$ is unambiguously identified with our NOEMA and Keck/MOSFIRE spectroscopy, where the NOEMA observations reveal the CO(4-3) and CO(5-4) lines at $>8\,\sigma$, and the MOSFIRE data detect [O\textsc{ii}] at $\sim 6\,\sigma$. Using multi-wavelength photometry spanning near-infrared to radio bands, we find that the lensed galaxy is a dust-obscured starburst ($M_{\star} \sim 1.8\times10^{10}\,{\rm M_{\odot}}$, ${\rm SFR_{IR}\sim 60\,{\rm M_{\odot}} ~yr^{-1}}$) with high star-formation efficiency (gas depletion time $\tau_{\rm dep}<100~$Myr) as indicated by the [C\textsc{i}](1-0) non-detection. The redshift confirmation revalidates that the total lens mass budget within the Einstein radius is fully accounted for by the stellar and dark matter components, without the need of modifying the initial mass function or dark matter distribution profile. This work paves the way for detailed studies and future follow-ups of this unique lensing system, providing an ideal laboratory for studying mass distribution at $z\sim2$ and physical conditions of star formation at $z\sim5$.

We show that the merger tree of dark matter halos is approximately self-similar by investigating the universality of the subhalo peak mass function (PMF) describing the mass distribution of progenitor halos. Using a set of cosmological simulations and identifying subhalos of different merger levels with HBT+, we verify that the level-1 subhalo PMF is close to universal across halo mass, redshift, and cosmology. This approximate self-similarity allows us to analytically derive the subhalo PMF for subhalos accreted at any level (i.e., for sub-sub...halos) through self-convolutions of the level-1 PMF, and the resulting model shows good agreement with simulation measurements. We further derive a number of analytical properties on the hierarchical origin of subhalos, including the level distribution, accretion rate at each level, initial merger ratio distribution, and accretion redshift distribution. We find that higher-level subhalos dominate at progressively lower peak mass in the PMF and are more likely to originate from major mergers than lower-level ones. At a given mass ratio at accretion time, the subhalo accretion rates at each level track the growth rate of the host halo. At a fixed final mass ratio, however, the accretion redshift distribution of subhalos depends on the subhalo level, peak mass, and host mass. Higher-level and higher-mass-ratio subhalos tend to be accreted more recently, and more massive halos also accrete their subhalos more recently. Our model provides a concise summary of simulation results and can serve as a basis for further theoretical understanding of the hierarchical structure formation.

In this paper we study the behavior of test particles on top of a galactic-type of Fuzzy Dark Matter (FDM) structure, characterized by the core-halo density profile found in simulations. Our workhorse structure is an anisotropic, time-dependent, virialized core-tail FDM clump resulting from a multicore merger. For our analysis we allow this structure to keep evolving, which implies that the core oscillates and accretes matter from the halo, while the halo dynamics is dominated by its characteristic high kinetic energy. On top of this time-dependent structure that in turn has a time-dependent gravitational potential, we solve the motion equations of test particles with initial conditions associated to circular orbits at different radii. Our results indicate that: 1) no trajectory remains circular, 2) the trajectories are sensitive to initial conditions and 3) the departure of initially near trajectories has always a positive Lyapunov exponent. A qualitative result is that the motion of test particles is more erratic with a bigger Lyapunov exponent within and near the core than in the halo region, which can be understood in terms of the random motion of the core within the core-halo structure. We expect these results warn on the importance of the anisotropic and time-dependent nature of FDM clumps when studying the motion of test particles.

We report the discovery of a metal-rich red giant star, WINERED-HVS1, which is a candidate for a hyper-velocity star (HVS). Its past trajectory suggests that this star may have been ejected by the Galactic supermassive black hole (SMBH; Sgr A*) about 14 Myr ago, with a modest ejection velocity of at least 500 km/s. Since WINERED-HVS1 is gravitationally bound to the Milky Way and is not necessarily young, it is not an unambiguous HVS candidate from a kinematic perspective, unlike the previously confirmed A-type main-sequence HVS known as S5-HVS1. Therefore, the chemical characterization of this star is essential for understanding its origin. Through a spectroscopic follow-up observation with Magellan/WINERED, we determined its metallicity [Fe/H]=$-0.136^{+0.040}_{-0.038}$, as well as the abundances for Na, Mg, Si, S, K, Ca, Ti, Cr, Ni, and Sr. The high metallicity value suggests that this star was ejected from the Galactic center and is unlikely to be a halo star with a radial orbit. The abundance pattern of this star -- including the pattern of [$\alpha$/Fe], [Mg/Fe], [Si/Mg], and [Ca/Mg] -- is consistent with that of the nuclear star cluster surrounding the SMBH, further supporting our view. This discovery opens a new window to look through the Galactic center environment without the hindrance of dust extinction, by using the HVSs moving in the halo region where dust extinction is minimal.

Accurately inferring black hole spin is crucial for understanding black hole dynamics and their astrophysical environments. In this work, we outline a geometric method for spin estimation by using the interferometric shape of the first photon ring ($n=1$) as an approximation to the critical curve, which, given an assumed value of the black hole inclination, is then mapped to a spin value. While future space-based missions will capture a wealth of data on the first photon ring--including the full angle-dependent diameter, angular brightness profile, and astrometric offset from the $n=0$ ring--our analysis is restricted to using only two angle-dependent diameters to compute its shape asymmetry and infer spin. Focusing on low inclinations and moderate-to-high spins, we test the method across various emission models, baselines, and noise sources, including a mock space-based observation. Although the size of the $n=1$ ring depends on the emission model, its interferometric shape remains a robust spin probe at low inclinations. We find that the inferred asymmetry of the $n=1$ image may be biased by the critical curve morphology, and it can be heavily skewed by the presence of noise, whether astrophysical or instrumental. In low-noise limits at low viewing inclination, significant contributions from the n=0 image at short baselines may lead to a downward bias in asymmetry estimates. While our method can estimate high spins in noise-free time-averaged images, increasing the noise and astrophysical variability degrades the resulting constraints, providing only lower bounds on the spin when applied to synthetic observed data. Remarkably, even using only the ring's asymmetry, we can establish lower bounds on the spin, underscoring the promise of photon ring-based spin inference in future space-based very long baseline interferometry (VLBI) missions, such as the proposed Black Hole Explorer (BHEX).

We identify observational signatures suggesting a history of dynamical instability in 26 out of 34 M-dwarf multi-planet systems containing no large planets. These systems may have primarily formed in a gas-rich environment, potentially hosted more planets and were more compact. We extend previous simulations of the formation of the TRAPPIST-1 system to 100 Myr to test the stability of these systems without gas. We find the absence of a strong mean motion resonance in the innermost planet pair and the absence of three body resonances throughout the system are likely to result in the merging and ejection of planets after the gas disk disperses. The runs that experience such an instability tend to produce final systems with lower multiplicities, period ratios larger than two, increased orbital spacings, higher planetary angular momentum deficits, and slightly smaller mass ratios between adjacent planets. Remarkably, we find these same trends in the observations of M-dwarf multi-planet systems containing no large planets. Our work allows us to identify specific systems that may have experienced an instability and suggests that only ~25% of these systems formed in their current observed state while most systems were likely more compact and multiplicitous earlier in time. Previous research indicates that systems that have experienced a late stage giant impact may host planets potentially more habitable than the systems that did not. With this in mind, we suggest systems around M-dwarfs that contain period ratios larger than two be given priority in the search for habitable worlds.

In this paper, we conduct a detailed study on the effect of Radiative Torque Disruption (RATD) mechanism on the fragmentation of micrometer-sized dust grains into nanoparticles within the heliosphere. We start by estimating the disruption timescales for dust grains under various centrifugal stresses. Our numerical calculations demonstrate that RATD is a highly effective mechanism for breaking down micrometer-sized grains, producing nanoparticles more efficiently than other fragmentation processes. RATD also prevents micrometer-sized grains from being expelled by radiation pressure. Our findings indicate that the location of the present water snow line depends not only on temperature but also on the size of dust grains. For smaller grains, the snow line can shift outward beyond the position defined by thermal sublimation. Furthermore, we model the size distribution of dust grains modified by the RATD mechanism using a simplified model, showing that rotational disruption significantly decreases the number density of micrometer-sized grains while substantially increasing the number density of sub-micrometer-sized grains. However, the fraction of dust grains aligned at high-$J$ attractors by radiative torques less than 80\% can considerably weaken the effect of RATD on the grain size distribution. Finally, we suggest several experiments that could potentially test the RATD mechanism and discuss the uncertainties of our model in more realistic applications to heliospheric dust.

We investigate linear matter density perturbations in the $\Lambda_{\rm s}$CDM, in which the $\Lambda$ is replaced by late-time ($z\sim2$) mirror AdS-dS transition, resulting in distinct growth dynamics. We use two complementary approaches: (i) determining the initial density contrast and its evolution rate for a given collapse scale factor, (ii) computing the collapse scale factor for a specified initial density contrast and evolution rate. We derive analytical solutions for the growth rate $f=\Omega_{\rm m}^\gamma$ and growth index $\gamma$ in both models. Prior to the transition, the AdS-like $\Lambda$ reduces cosmic friction, causing linear matter density perturbations to grow more rapidly than in $\Lambda$CDM; this effect is most pronounced just before the transition, with a growth rate around $15\%$ higher than $\Lambda$CDM around $z\sim2$. After the transition, $\Lambda_{\rm s}$CDM behaves similarly to $\Lambda$CDM but features a larger cosmological constant, leading to higher $H(z)$ and greater cosmic friction that more effectively suppresses growth. Before the transition, $\gamma$ remains below both the $\Lambda$CDM and EdS values ($\gamma\approx6/11$); during the transition, it increases rapidly and then grows gradually, paralleling $\Lambda$CDM while remaining slightly higher in the post-transition era-though overall, it stays near $\gamma\sim0.55$. Using the Planck best-fit values, $\Omega_{\rm m0}=0.28$ for $\Lambda_{\mathrm{s}}$CDM and $\Omega_{\rm m0}=0.32$ for $\Lambda$CDM, we find $f=0.49$ and $f=0.53$, respectively. $\Lambda_{\rm s}$CDM predicts a value of $f=0.48$, recently obtained from LSS data when $\gamma$ is treated as a free parameter in $\Lambda$CDM. This suggests that $\Lambda_{\rm s}$CDM may resolve the structure growth anomaly, without deviating from $\gamma \sim 0.55$.

Galaxy sizes are a key parameter to distinguishing between different galaxy types and morphologies, reflecting their formation and assembly histories. Several methods define galaxy boundaries, often relying on light concentration or isophotal densities. However, these approaches were often constrained by observational limitations and did not necessarily provide a clear physical boundary for galaxy outskirts. With modern deep imaging surveys, a new physically motivated definition has emerged using the radial position of the star formation threshold as the galaxy size, approximated by the stellar mass density contour at 1 Msun pc^-2 (R_1). We test this definition using three state-of-the-art hydrodynamical simulation suites, analyzing stellar surface density profiles across a wide range of stellar masses and redshifts. We measure the galaxy sizes according to this new definition and compare them with the most traditional size metric, the stellar half-mass radius. Our analysis demonstrates that the R_1-M_star relation exhibits consistent behaviour across both low and high-stellar mass galaxies, with remarkably low scatter. This relation is independent of redshift and holds across the three different cosmological hydrodynamical simulation suites, highlighting its robustness to variations in galaxy formation models. Furthermore, we explore the connection between a galaxy's total mass within R1 and its stellar mass, finding very little scatter in this relation. This suggests that R1 could serve as a reliable observational tracer for the galaxy's dynamical mass. The size-stellar mass relation proposed provides a reliable and physically motivated method for defining the outskirts of galaxies. This method remains consistent not only at z=0 but also throughout the evolutionary history of galaxies, offering a robust and meaningful framework for galaxy evolution studies.

A non iterative direct blind deconvolution procedure, previously used successfully to sharpen Hubble Space Telescope imagery, is now found useful in sharpening nanoscale scanning electron microscope (SEM) and helium ion microscope (HIM) images. The method is restricted to images $g(x,y)$, whose Fourier transforms $\hat{g}(\xi,\eta)$ are such that $log~|\hat{g}(\xi,0)|$ is globally monotone decreasing and convex. The method is not applicable to defocus blurs. A point spread function in the form of a Linnik probability density function is postulated, with parameters obtained by least squares fitting the Fourier transform of the preconditioned microscopy image. Deconvolution is implemented in slow motion by marching backward in time, in Fourier space, from $t = 1$ to $t = 0$, in an associated logarithmic diffusion equation. Best results are usually found in a partial deconvolution at time $\bar{t}$, with $0 < \bar{t} < 1$, rather than in total deconvolution at $t=0$. The method requires familarity with microscopy images, as well as interactive search for optimal parameters.

In a series of two papers, we provide a comprehensive agenda of actions the American Astronomical Society (AAS) can take to create a more diverse and inclusive professional system for astronomers, with a focus on women astronomers. This first paper of the series outlines the background and methods, while the recommendations are treated in the second companion paper (Paper II). We take the stance that since the 2020 Decadal Survey (Astro2020) was delivered in 2021, with its first-ever set of recommendations on the State of the Profession, now is the time for the AAS to take decisive action to transform astronomy into a diverse and inclusive profession. In the spring of 2019, the CSWA surveyed the astronomical community to assess the popularity and feasibility of actions that the AAS can take to reduce harassment and advance career development for women in astronomy. Here we present the quantitative results of that survey and a synopsis of the free response sections, which are publicly accessible. By combining the results of our survey, peer-reviewed academic literature, and findings from many of the white papers submitted to Astro2020, the CSWA has developed 26 specific actions that the AAS can take to help end harassment in astronomy, to advance career development for astronomers who are women and who are other members of historically marginalized groups, and intersections of these populations, and to improve the climate and culture of AAS and AAS-sponsored meetings. This paper presents the data we used to make these recommendations, and the recommendations themselves will be presented in Paper II.

This paper, the second in a series of two, provides a set of recommendations that the American Astronomical Society (AAS) can take to create a more diverse and inclusive professional society for astronomers, with a focus on women astronomers. As noted in Paper I, now is the time for the AAS to take decisive action to transform astronomy into a diverse and inclusive profession. By combining the results of our 2019 survey, which is described in Paper I, peer-reviewed academic literature, and findings from many of the white papers submitted to Astro2020, the CSWA has developed 26 specific actions the AAS can take to help end harassment and bullying in astronomy; advance career development for astronomers who are women, members of other underrepresented groups, and intersections of these populations; and improve the climate and culture of AAS meetings. Actions to reduce rates of harassment and bullying include improvements to the AAS's anti-harassment policies and procedures and the development of astronomy-specific anti-harassment training resources. Actions to advance career development include creating a compensation database, improving how jobs are posted in the AAS Job Register, and supporting/enhancing a distance mentorship program. Finally, we call on the AAS to continue improving the accessibility of AAS meetings and to continue to support meeting sessions whose focus is to discuss issues of equity, diversity, and inclusion.

We investigate the potential of IAXO and its intermediate version, BabyIAXO, to detect axions produced in core-collapse supernovae (SNe). Our study demonstrates that these experiments have realistic chances of identifying SN axions, offering crucial insights into both axion physics and SN dynamics. IAXO's sensitivity to SN axions allows for the exploration of regions of the axion parameter space inaccessible through solar observations. In addition, in the event of a nearby SN, $d \sim O(100)$ pc, and sufficiently large axion couplings, $g_{a \gamma} \gtrsim 10^{-11} GeV^{-1}$, IAXO could have a chance to significantly advance our understanding of axion production in nuclear matter and provide valuable information about the physics of SNe, such as pion abundance, the equation of state, and other nuclear processes occurring in extreme environments.

Effective field theories (EFTs) parametrize our ignorance of the underlying UV theory through their Wilson coefficients. However, not all values of these coefficients are consistent with fundamental physical principles. In this paper, we explore the consequences of imposing causal propagation on the comoving curvature perturbation in the EFT of inflation, particularly its impact on the primordial power spectrum and the effective sound speed $c_s^\text{eff}$. We investigate scenarios where $c_s^\text{eff}$ undergoes a transition, remaining consistent with CMB constraints at early times but later experiencing a drastic change, becoming highly subluminal. Such scenarios allow the primordial power spectrum to grow at small scales, potentially leading to the formation of primordial black holes or the generation of scalar-induced gravitational waves. We find the generic feature that in a causal theory, luminal sound speeds imply a free theory, effectively constraining the dynamics. Additionally, we obtain that when considering natural values for the Wilson coefficients, maintaining the validity of the EFT and the weakly coupled regime, and enforcing causal propagation of the EFT modes, the power spectrum cannot increase drastically. This imposes significant constraints on the parameter space of models aiming to produce such features.

Strong first-order phase transitions in a dark sector offer a compelling explanation for the stochastic gravitational wave background in the nano-Hertz range recently detected by pulsar timing arrays (PTAs). We explore the possibility that such a phase transition at the same time gives mass to a stable fermion that accounts for the observed dark matter abundance and leads to testable effects in laboratory experiments. Concretely, we consider a classically conformal dark sector with a hidden $U(1)^\prime$ gauge symmetry that couples to the Standard Model via kinetic mixing. Since the PTA signal requires a phase transition in the MeV temperature range, spontaneous symmetry breaking gives rise to a sub-GeV dark matter candidate that couples to the Standard Model via a dark photon mediator and obtains its relic abundance via annihilations into electrons and dark Higgs bosons. Such a scenario is tightly constrained by laboratory searches for dark photons and cosmological constraints on the decays of dark Higgs bosons after the phase transition. We show that viable parameter regions can be found both for the case that the dark Higgs bosons remain in equilibrium with the Standard Model and that they decouple and only decay much later. In the latter case, the parameter regions preferred by the PTA signal and the dark matter relic abundance can be fully explored by future beam-dump experiments searching for missing energy.

We investigate the role of an axion-like particle (ALP) as a portal between the dark and visible sectors. Unlike conventional studies, which typically assume fermionic dark matter (DM), we explore the phenomenological implications of scalar DM within this ALP portal framework. A key challenge arises from the fact that the interaction between the ALP spacetime derivative and the spin-one current of the scalar DM can be a redundant operator, which may be removed via a field redefinition. However, this interaction reveals a profound connection to the underlying global symmetry that stabilizes the DM particle. We choose a non-Abelian discrete symmetry, ensuring the persistence of the DM-ALP interactions, and in doing so, unveil a rich phenomenology. Working within a general effective field theory approach, we identify the following hallmark features of our scenario: (i) a relic density determined by semi-annihilations, with an abundance independent of ALP couplings to the visible sector; (ii) direct detection rates naturally suppressed; (iii) indirect detection spectra enriched relative to pure annihilation scenarios, with rates also independent of ALP couplings to visible particles. Lastly, we discuss potential microscopic origins for this framework and highlight the broader implications of our results.

The density of the Earth's middle and upper atmosphere is an important question in Earth science and is a critical factor in the design, operation, and orbital determination of low Earth orbit spacecraft. In this study, we employ the Earth Occultation Technique (EOT) combined with Maximum Likelihood Estimation to estimate the neutral atmospheric density by modeling the attenuation of X-ray photons during the occultation process of \textit{Insight}-HXMT observations of Crab Nebula. Based on 83 occultation datasets of the Crab Nebula observed by all three sets of telescopes of \textit{Insight}-HXMT between 2022 and 2024, we derived the atmospheric densities at altitudes ranging from 55\,--130\,km. We find a general agreement between our results and the prediction by the NRLMSIS model within the altitude ranges of 65\,-- 90\,km, 95\,--100\,km and 120\,--130\,km, particularly during periods of enhanced solar activity. However, we also find that the NRLMSIS model overestimates atmospheric density at altitudes 90\,--95\,km and 100\,--120\,km by approximately 20\%. Furthermore, since the atmospheric density measurements at altitudes of 55\,--\,65\,km may be subject to selection bias, we do not report the prediction accuracy of the NRLMSIS model at this altitude.

The isotopes of manganese in the mass range A equal to 53 to 63 are abundant in the core material of high mass stars and are believed to be of prime importance in the progression of the pre collapse phases. During these late evolutionary phases, nuclear processes associated with weak interactions, including decay and electron capture (EC) on these isotopes, significantly alter the Ye (lepton to baryon ratio) of the cores composition. The temporal change of this parameter is one of the basic elements to simulate a successful explosion. Hence, the decay and EC rates of manganese (Mn) nuclides may serve as an important input in the simulation codes of core collapse supernova. In this paper, we focus on the study of the weak decay characteristics of 55 63Mn nuclides. The strength distributions of Gamow Teller transitions in the directions of decay and EC for these isotopes were calculated employing the proton neutron quasi particle random phase approximation (pn QRPA) model. The decay half life values of Mn isotopes under terrestrial conditions were computed and compared with measured data and previous calculations. The pn QRPA estimated half lives are in good agreement with the experimental values. The and EC rates were later calculated covering a large range of stellar temperatures (001 30) 109 K and densities (10 1011) g/cm3.

We reconstruct a holographic dark energy model with a scalar torsion $\phi$ under the assumption of no interaction between dark energy and dark matter. We show that the accelerating expansion can occur by setting the IR cut-off to the Hubble radius as $L=H^{-1}$ with the Hubble parameter $H$. Motivated by two physical choices based on Friedmann equations and the holographic principle, tellingly, spin-induced torsion at non-zero matter density $\rho_m \neq 0$ and $H$-dependence of torsion in empty space $\rho_m =0$, we propose a time-dependent scalar torsion $\phi(t) = k H(t) ({\rho_{m}(t)}/ \rho_c^0)$ with a dimensionless constant $k$ and the current critical density $\rho_c^0$. As a result, we find minima $(\omega_X^{0})_{min}$ of the current equation of state for dark energy for various values of the free parameter $c$, yielding $-1 < (\omega_X^{0})_{min} < -0.778$ as $c$ changes from $1$ to $0.654$. Also, we get the values of the current dimensionless ratio, $(\phi_0/H_{0})_{min}$ corresponding to the minima $(\omega^{0}_{X})_{min}$, which is consistent with weak torsion assumption $|\phi_0| /H_{0} < 1$, and thus it ensuring model viability.

[Background] Bayesian inference frameworks incorporating multi-messenger astrophysical constraints have recently been applied to covariant density functional (CDF) models to constrain their parameters. Among these, frameworks utilizing CDFs with density-dependent meson-nucleon couplings furnishing the equation of state (EoS) of compact star (CS) matter have been explored. [Purpose] The aforementioned inference framework has not yet incorporated astrophysical objects with potentially extreme high masses or ultra-small radii among its constraints, leaving its flexibility and predictive power under such extreme parameters still unknown. [Method] We apply the Bayesian inference framework based on CDFs with density dependent couplings. The astrophysical data is expanded to include not only the latest multi-messenger constraints from NICER and gravitational wave events but also the highest measured mass to date for the ``black widow" pulsar PSR J0952-0607 and the mass-radius estimates for the ultra-compact, low-mass object HESS J1731-347. [Results] Our systematic Bayesian analysis indicates that our CDF models can support higher maximum masses for CSs, reaching up to $2.4$-$2.5\,M_{\odot}$. However, achieving sufficient softening of the EoS in the low-density regime to accommodate the HESS J1731-347 data remains challenging. Nonetheless, we are able to impose tighter constraints on the parameter space of CDF models, ensuring consistency with current nuclear experimental and astrophysical data. [Conclusions] CDF models with density-dependent meson-nucleon couplings encompass a wide range of nuclear and astrophysical phenomena, providing a robust theoretical framework for interpreting compact objects. However, the predicted lower limit for the radii of low-mass stars is approximately 12 km, which stems from the restricted degrees of freedom in the isovector sector.

We study turbulence in self-gravitating superfluids by performing direct numerical simulations of the 3D Gross-Pitaevskii-Poisson (GPP) equation, which is also a model for dark matter haloes around galaxies. In the absence of self-gravity, the spectrally truncated Gross-Pitaevskii (GP) equation shows the emergence of Kolmogorov's $5/3$ scaling in the incompressible kinetic energy spectrum. Introducing self-gravity, we observe the formation of spherically collapsed structures, which introduce a minimum in the kinetic energy spectrum that corresponds to the sizes of these structures. The system shows early convergence towards statistically stationary states, which we show by the onset of thermalisation in the compressible kinetic energy spectrum, where $E_{\rm kin}^c \propto k^2$. We also show that the formation of such large-scale structures suggests that the particles (bosons) move from small to large scales through an inverse cascade, supporting a mechanism for the formation of large-scale structures, such as dark matter haloes, around our galaxy Milky Way.

Scalar-Induced Gravitational Waves (SIGWs), second-order tensor modes sourced by first-order scalar fluctuations in General Relativity (GR), are expected to contribute to the Stochastic Gravitational Wave Background (SGWB) potentially detectable by current and future gravitational wave interferometers. In the framework of GR, this SGWB represents an unavoidable contribution to the gravitational wave spectrum. In this paper, we go beyond GR and we investigate the behavior of SIGWs in f(R) gravity. We explore how the SIGW spectrum is influenced across a broad range of frequencies, from the nano-Hz regime, where the Pulsar Timing Array (PTA) operates, through the milli-Hz band probed by the space-based LISA detector, up to the kilo-Hz frequency range, where the ground-based LIGO/Virgo/KAGRA network is currently operational. Our results indicate that the beyond-GR correction leaves an observational imprint, mainly in the low-frequency part of the spectrum, giving the possibility to use SIGW to constrain GR on scales on which we have limited information.

The exploration of the surrounding world and the universe is an important theme in the legacy of humankind. The detection of gravitational waves is adding a new dimension to this grand effort. What are the fundamental physical laws governing the dynamics of the universe? What is the fundamental composition of the universe? How has the universe evolved in the past and how will it evolve in the future? These are the basic questions that press for answers. The space-based gravitational wave detector TianQin will tune in to gravitational waves in the millihertz frequency range ($10^{-4} \sim 1$ Hz, to be specific), opening a new gravitational wave spectrum window to explore many of the previously hidden sectors of the universe. TianQin will discover many astrophysical systems, populating the universe at different redshifts: some will be of new types that have never been detected before, some will have very high signal-to-noise ratios, and some will have very high parameter estimation precision. The plethora of information collected will bring us to new fronts on which to search for the breaking points of general relativity, the possible violation of established physical laws, the signature of possible new gravitational physics and new fundamental fields, and to improve our knowledge on the expansion history of the universe. In this white paper, we highlight the advances that TianQin can bring to fundamental physics and cosmology.

Convection is the main heat transport mechanism in the Earth's liquid core and is thought to power the dynamo that generates the geomagnetic field. Core convection is strongly constrained by rotation while being turbulent. Given the difficulty in modelling these conditions, some key properties of core convection are still debated, including the dominant energy-carrying lengthscale. Different regimes of rapidly-rotating, unmagnetised, turbulent convection exist depending on the importance of viscous and inertial forces in the dynamics, and hence different theoretical predictions for the dominant flow lengthscale have been proposed. Here we study the transition from viscously-dominated to inertia-dominated regimes using numerical simulations in spherical and planar geometries. We find that the cross-over occurs when the inertial lengthscale approximately equals the viscous lengthscale. This suggests that core convection in the absence of magnetic fields is dominated by the inertial scale, which is hundred times larger than the viscous scale.

We conducted three-dimensional lattice simulations to study the density perturbation and generation of gravitational waves (GWs) through first-order phase transition (FOPTs). Our findings reveal that for the case of $\beta/H = 6$, density perturbations successfully collapse into PBHs. Moreover, smaller values of $\beta/H$ are more conducive to PBH formation, while the strength of the phase transition $\alpha$ has a relatively minor impact on PBH production. We find that for $\alpha > 1$, the dominant contribution to the density perturbation comes from the delay of vacuum decay, whereas for $\alpha < 1$, the bubble walls' moving forward behavior becomes the primary source. Additionally, the power spectrum of density perturbations generated by the phase transition exhibits a slope of $k^3$ at small wavenumbers and $k^{-1.5}$ at large wavenumbers. Furthermore, we calculated the GW power spectra, and the simulation results show a slope of $k^3$ at small wavenumbers and $k^{-2}$ at large wavenumbers. Our numerical simulations confirm that slow PTs can produce PBHs and provide predictions for the GW power spectrum, offering crucial theoretical support for GW detection.

We study the stability properties of multi-state configurations of the Schrödinger-Poisson system without self-interaction, with monopolar and first dipolar components $(1,0,0)$+$(2,1,0)$. We show these configurations studied are stable using numerical simulations, and using criteria of stationarity, unitarity and time dependence consistency. The study covers a range of states with monopolar to dipolar mass ratio between 47 and 0.17. The astrophysical implication of this result is that this type of configurations is at least stable and can be considered physically sound in multi-state ultralight bosonic dark matter.

We have conducted an extensive study using a diverse set of equations of state (EoSs) to uncover strong relationships between neutron star (NS) observables and the underlying EoS parameters using symbolic regression method. These EoS models, derived from a mix of agnostic and physics-based approaches, considered neutron stars composed of nucleons, hyperons, and other exotic degrees of freedom in beta equilibrium. The maximum mass of a NS is found to be strongly correlated with the pressure and baryon density at an energy density of approximately 800 this http URL$^{-3}$. We have also demonstrated that the EoS can be expressed as a function of radius and tidal deformability within the NS mass range 1-2$M_\odot$. These insights offer a promising and efficient framework to decode the dense matter EoS directly from the accurate knowledge of NS observables.

Adaptive mesh refinement efficiently facilitates the computation of gravitational waveforms in numerical relativity. However, determining precisely when, where, and to what extent to refine when solving the Einstein equations poses challenges; several ad hoc refinement criteria have been explored in the literature. This work introduces an optimized resolution baseline derived in situ from the inspiral trajectory (ORBIT). This method uses the binary's orbital frequency as a proxy for anticipated gravitational waves to dynamically refine the grid, satisfying the Nyquist frequency requirements on grid resolution up to a specified spin weighted spherical harmonic order. ORBIT sustains propagation of gravitational waves while avoiding the more costly alternative of maintaining high resolution across an entire simulation, both spatially and temporally. We find that enabling ORBIT decreases waveform noise by an order of magnitude and better resolves high-order wave amplitudes through merger. Combined with WAMR and other improvements, updates to Dendro-GR decrease waveform noise, decrease constraint violations, and boost refinement efficiency each by factors of $\mathcal{O}(100)$, while reducing computational cost by a factor of four. ORBIT and other recent improvements to Dendro-GR begin to prepare us for gravitational wave science with next-generation detectors.