Papers for Thursday, Dec 05 2024

The conventional picture of supermassive black-hole growth in the standard model had already been seriously challenged by the emergence of $\sim 10^9\;M_\odot$ quasars at $z\sim 7.5$, conflicting with the predicted formation of structure in the early $\Lambda$CDM Universe. But the most recent {\it JWST} discovery of a $\sim 10^8\;M_\odot$ source at $z\sim 10.1$ argues even more strongly against the possibility that these black holes were created in Pop II or III supernovae, followed by Eddington-limited accretion. Attempts at resolving this anomaly have largely focused on the formation of seeds via an exotic, direct collapse of primordial gas to an initial mass $\sim 10^5\;M_\odot$ -- a process that has never been seen anywhere in the cosmos. Our goal in this {\it Letter} is to demonstrate that the emergence of these black holes is instead fully consistent with standard astrophysics in the context of the alternative Friedmann-Lemaître-Robertson-Walker cosmology known as the $R_{\rm h}=ct$ universe. We show that, while the predicted evolution in the standard model is overly compressed, the creation, growth and appearance of such high-$z$ quasars fall comfortably within the evolutionary history in this cosmology, thereby adding considerable observational support to the existing body of evidence favoring it over the standard scenario.

Weak gravitational lensing (WL) convergence peaks contain valuable cosmological information in the regime of non-linear collapse. Using the FLAMINGO suite of cosmological hydrodynamical simulations, we study the physical origin and redshift distributions of the objects generating WL peaks selected from a WL convergence map mimicking a $\textit{Euclid}$ signal. We match peaks to individual haloes and show that the high signal-to-noise ratio (SNR$~>~5$) WL peaks measured by Stage IV WL surveys primarily trace $M_{\mathrm{200c}} > 10^{14}~\mathrm{M_\odot}$ haloes. We find that the WL peak sample can compete with the purity and completeness of state-of-the-art X-ray and Sunyaev-Zel'dovich cluster abundance inferences. By comparing the distributions predicted by simulation variations that have been calibrated to the observed gas fractions of local clusters and the present-day galaxy stellar mass function, or shifted versions of these, we illustrate that the shape of the redshift distribution of SNR$~>~5$ peaks is insensitive to baryonic physics while it does change with cosmology. The difference highlights the potential of using WL peaks to constrain cosmology. As the number density of WL peaks scales differently with cosmology and baryonic feedback, WL peak statistics can simultaneously calibrate baryonic feedback and constrain cosmology.

Mass estimates derived from galactic rotation curves must account for past accretion histories of galaxies. There are several lines of evidence indicating that M31 experienced a major merger 2 to 3 Gyr ago. Here, we have generated a dynamical model of M31 as a merger remnant that reproduces most of its properties, from the central bar to the outskirts. The model accounts for the past major merger, and reproduces the details of M31's rotation curve, including its 14 kpc bump and the observed increase of velocity beyond 25 kpc. Furthermore, we find non-equilibrium and oscillatory motions in the gas of the merger-remnant outskirts caused by material in a tidal tail returning to the merger remnant. A total dynamical M31 mass of 4.5 $\times 10^{11} M_{\odot}$ within 137 kpc has been obtained after scaling it to the observed HI rotation curve. Within this radial distance, 68% of the total dynamical mass is dark.

At the core of the Gaia mission is a multi-epoch survey consisting of astrometric, photometric, spectrophotometric, and spectroscopic measurements. The astrometric time series provides parallaxes and proper motions, along with information on astrometric binary systems. The photometric time series offers a means to investigate the variability of the sources. Due to their whole-sky, multi-epoch nature, multiple instruments, their magnitude range covering 21 magnitudes, and their remarkable photometric precision, these data allow us to describe the variability of celestial phenomena in an unprecedented manner. For the third Gaia Data Release (DR3), the data collection spanned 34 months, with a median number of field-of-view measurements in the G band of about 44, reaching up to 270. At publication time, DR3 delivered the largest collection of variable sources with an associated classification across the entire sky. All these sources have their $G$, $G_{BP}$, $G_{RP}$ epoch data published and accessible in the Gaia ESA archive. In summary, there are 10.5 million variable sources, including 9.5 million variable stars and 1 million QSOs. Additionally, 2.5 million galaxies were identified thanks to spurious variability caused by the non-axisymmetric nature of galaxies and the way Gaia collects data. Moreover, all the epoch data and time series of nearly 1.3 million sources in a pencil beam around the Andromeda galaxy are published, regardless of their status (constant or variable); This dataset is known as the Gaia Andromeda Photometric Survey. We also introduce the citizen science project, GaiaVari to classify variable stars, the Focused Product Release delivered on October 10, 2023. In the future, DR4 will cover 66 months, and we hope DR5 will have accumulated 10.5 years of data.

The high angular resolution and sensitivity of Very Long Baseline Interferometry (VLBI) offer a unique tool to identify and study active galactic nuclei and star-formation activity over cosmic time. However, despite recent technical advances, such as multiple phase centre correlation, VLBI surveys have thus far been limited to either a few well-studied deep-fields or wide-areas to a relatively shallow depth. To enter the era of extensive statistical studies at high angular resolution, a significantly larger area of the sky must be observed to much better sensitivity with VLBI. The Synoptic Wide-field EVN--e-MERLIN Public Survey (SWEEPS) is a proposed commensal observing mode for the EVN and e-MERLIN, where single-target principle investigator-led observations are re-correlated at the position of known radio sources within 12 arcmin of the pointing centre. Here, we demonstrate a proof-of-concept of this methodology by detecting a 5.6 mJy core-jet object at 1.7 GHz that would have otherwise been lost from the parent data set. This is the first object to be recovered as part of the SWEEPS pilot programme, which highlights the potential for increasing sample sizes of VLBI-detected radio sources with commensal observing modes in the near future.

Due to poor observational constraints on the low-mass end of the subhalo mass function, the detection of dark matter (DM) subhalos on sub-galactic scales would provide valuable information about the nature of DM. Stellar wakes, induced by passing DM subhalos, encode information about the mass of the inducing perturber and thus serve as an indirect probe for the DM substructure within the Milky Way (MW). Our aim is to assess the viability and performance of deep learning searches for stellar wakes in the Galactic stellar halo caused by DM subhalos of varying mass. We simulate massive objects (subhalos) moving through a homogeneous medium of DM and star particles, with phase-space parameters tailored to replicate the conditions of the Galaxy at a specific distance from the Galactic center. The simulation data is used to train deep neural networks with the purpose of inferring both the presence and mass of the moving perturber, and assess subhalo detectability in varying conditions of the Galactic stellar and DM halos. We find that our binary classifier is able to infer the presence of subhalos, showing non-trivial performance down to a subhalo mass of $5 \times 10^7 \rm \, M_\odot$. We also find that our binary classifier is generalisable to datasets describing subhalo orbits at different Galactocentric distances. In a multiple-hypothesis case, we are able to discern between samples containing subhalos of different masses. Out of the phase-space observables available to us, we conclude that overdensity and velocity divergence are the most important features for subhalo detection performance.

We present a novel multi-messenger approach for probing nonstandard neutrino properties through the detection of gravitational waves (GWs) from collapsing stellar cores and associated supernova explosions. We show that neutrino flavor conversion inside the proto-neutron star (PNS), motivated by physics Beyond the Standard Model (BSM), can significantly boost PNS convection. This effect leads to large-amplitude GW emission over a wide frequency range during an otherwise relatively quiescent GW phase shortly after core bounce. Such a signal provides a promising new avenue for exploring nonstandard neutrino phenomena and other BSM physics impacting PNS convection.

Distant quasars (QSOs) in galaxy overdensities are considered key actors in the evolution of the early Universe. In this work, we studied the kinematic and physical properties of the BR1202-0725 system at z=4.7, one of the most overdense fields known in the early Universe, consisting of a QSO, a submillimeter galaxy (SMG), and three Lyman-$\alpha$ emitters. We used data from the JWST/NIRSpec Integral Field Unit (IFU) to analyze the rest-frame optical emission of each source in the system. We estimated a bolometric luminosity of log($L_{\rm bol}/$[erg/s]) = 47.2 $\pm$ 0.4 and a black hole mass of log($M_{\rm BH}/M_\odot$) = 10.1 $\pm$ 0.5 for the QSO, which are consistent with previous measurements obtained with ground-based observations. The NIRSpec spectra of the SMG revealed instead unexpected [OIII] and H$\alpha$+[NII] profiles. The overall [OIII] line profile is blue-shifted by more than 700 km/s relative to the systemic velocity of the galaxy. Additionally, both the [OIII] and H$\alpha$+[NII] lines show prominent broad (1300 km/s), blueshifted wings associated with outflowing ionized gas. The analysis of NIRSpec and X-ray observations indicates that the SMG likely hosts an accreting supermassive black hole as supported by the following results: (i) the excitation diagnostic diagram is consistent with ionization from an active galactic nucleus (AGN); (ii) the X-ray luminosity is higher than $10^{44}$ erg/s; and (iii) it hosts a fast outflow ($v_{\rm out}$ = 5000 km/s), comparable to those observed in luminous QSOs. Therefore, the QSO-SMG pair represents one of the highest-redshift double AGN to date, with a projected separation of 24 kpc. Finally, we investigated the environment of this system and found four new galaxies at the same redshift of the QSO and within a projected distance of 5 kpc from it. This overdense system includes at least ten galaxies in only 980 kpc$^2$.

This paper investigates the origin and orbital evolution of S-stars in the Galactic Center using models of binary disruption and relaxation processes. We focus on explaining the recently discovered "zone of avoidance" in S-star orbital parameters, defined as a region where no S-stars are observed with pericenters $\log(r_p / {\rm AU}) \leq 1.57 + 2.6(1 - e)$ pc. We demonstrate that the observed S-star orbital distributions, including this zone of avoidance and their thermal eccentricity distribution, can be largely explained by continuous disruption of binaries near the central supermassive black hole, followed by orbital relaxation. Our models consider binaries originating from large scales (5--100 pc) and incorporate empirical distributions of binary properties. We simulate close encounters between binaries and the black hole, tracking the remnant stars' orbits. The initially highly eccentric orbits of disrupted binary remnants evolve due to non-resonant and resonant relaxation in the Galactic Center potential. While our results provide insights into the formation mechanism of S-stars, there are limitations, such as uncertainties in the initial binary population and mass-function and simplifications in our relaxation models. Despite these caveats, our study demonstrates the power of using S-star distributions to probe the dynamical history and environment of the central parsec of our Galaxy.

Accretion-induced collapse (AIC) or merger-induced collapse (MIC) of white dwarfs (WDs) in binary systems is an interesting path to neutron star (NS) and magnetar formation, alternative to stellar core collapse and NS mergers. Such events could add a population of compact remnants in globular clusters, they are expected to produce yet unidentified electromagnetic transients including gamma-ray and radio bursts, and to act as sources of trans-iron elements, neutrinos, and gravitational waves. Here we present the first long-term (>5 s post bounce) hydrodynamical simulations in axi-symmetry (2D), using energy- and velocity-dependent three-flavor neutrino transport based on a two-moment scheme. Our set of six models includes initial WD configurations for different masses, central densities, rotation rates, and angular momentum profiles. Our simulations demonstrate that rotation plays a crucial role for the proto-neutron star (PNS) evolution and ejecta properties. We find early neutron-rich ejecta and an increasingly proton-rich neutrino-driven wind at later times in a non-rotating model, in agreement with electron-capture supernova models. In contrast to that and different from previous results, our rotating models eject proton-rich material initially and increasingly more neutron-rich matter as time advances, because an extended accretion torus forms around the PNS and feeds neutrino-driven bipolar outflows for many seconds. AIC and MIC events are thus potential sites of r-process element production, which may imply constraints on their occurrence rates. Finally, our simulations neglect the effects of triaxial deformation and magnetic fields, serving as a temporary benchmark for more comprehensive future studies.

Gaia has revealed a clear signal of disequilibrium in the solar neighborhood in the form of a spiral (or snail) feature in the vertical phase-space distribution. We investigate the possibility that this structure emerges from ongoing perturbations by dark $\left(10^{6} M_{\odot} - 10^8 M_{\odot}\right)$ Galactic subhalos. We develop a probabilistic model for generating subhalo orbits based on a semi-analytic model of structure formation, and combine this framework with an approximate prescription for calculating the response of the disk to external perturbations. We also develop a phenomenological treatment for the diffusion of phase-space spirals caused by gravitational scattering between stars and giant molecular clouds, a process that erases the kinematic signatures of old ($t \gtrsim 0.6$ Gyr) events. Perturbations caused by dark subhalos are, on average, orders of magnitude weaker than those caused by luminous satellite galaxies, but the ubiquity of dark halos predicted by cold dark matter makes them a more probable source of strong perturbation to the dynamics of the solar neighborhood. Dark subhalos alone do not cause enough disturbance to explain the Gaia snail, but they excite fluctuations of $\sim 0.1-0.5 \ \rm{km} \ \rm{s^{-1}}$ in the mean vertical velocity of stars near the Galactic midplane that should persist to the present day. Subhalos also produce correlations between vertical frequency and orbital angle that could be mistaken as originating from a single past disturbance. Our results motivate investigation of the Milky Way's dark satellites by characterizing their kinematic signatures in phase-space spirals across the Galaxy.

The cross-correlation of cosmic voids with the lensing convergence ($\kappa$) map of the Cosmic Microwave Background (CMB) fluctuations provides a powerful tool to refine our understanding of the cosmological model. However, several studies have reported a moderate tension between the lensing imprint of cosmic voids on the observed CMB and the simulated $\mathrm{\Lambda}$CDM signal. To address this "lensing-is-low" tension and to obtain new, precise measurements, we exploit the large DESI Legacy Survey Luminous Red Galaxy (LRG) dataset, covering approximately 19,500 $°^2$ of the sky and including about 10 million LRGs at $z < 1.05$. Our $\mathrm{\Lambda}$CDM template was created using the Buzzard mocks, which we specifically calibrated to match the clustering properties of the observed galaxy sample by exploiting more than one million DESI spectra. We identified our catalogs of 3D voids in the range $0.35 < z < 0.95$, dividing the sample into bins according to the redshift and $\lambda_\mathrm{v}$ values of the voids. We report a 14$\sigma$ detection of the lensing signal, with $A_\kappa = 1.016 \pm 0.054$, which increases to 17$\sigma$ when considering the void-in-void ($A_\kappa = 0.944 \pm 0.064$) and the void-in-cloud ($A_\kappa = 0.975 \pm 0.060$) populations individually, the highest detection significance for studies of this kind. We observe a full agreement between the observations and $\mathrm{\Lambda}$CDM predictions across all redshift bins, sky regions, and void populations considered. In addition to these findings, our analysis highlights the importance of matching sparseness and redshift error distributions between mocks and observations, as well as the role of $\lambda_\mathrm{v}$ in enhancing the signal-to-noise ratio.

Recent James Webb Space Telescope observations have unveiled that the first supermassive black holes (SMBHs) were in place at z $\geq$ 10, a few hundred Myrs after the Big Bang. These discoveries are providing strong constraints on the seeding of BHs and the nature of the first objects in the Universe. Here, we study the impact of the freeze-out electron fractions ($f_e$) at the end of the epoch of cosmic recombination on the formation of the first structures in the Universe. At $f_e$ below the current fiducial cosmic values of $\rm \sim 10^{-4}$, the baryonic collapse is delayed due to the lack of molecular hydrogen cooling until the host halo masses are increased by one to two orders of magnitude compared to the standard case and reach the atomic cooling limit. This results in an enhanced enclosed gas mass by more than an order of magnitude and higher inflow rates of up to $0.1~M_{\odot}/{yr}$. Such conditions are conducive to the formation of massive seed BHs with $\sim 10^{4}$ M$_{\odot}$. Our results reveal a new pathway for the formation of massive BH seeds which may naturally arise from free

We present results from a Keck/DEIMOS survey to study satellite quenching in group environments at $z \sim 0.8$ within the Extended Groth Strip (EGS). We target $11$ groups in the EGS with extended X-ray emission. We obtain high-quality spectroscopic redshifts for group member candidates, extending to depths over an order of magnitude fainter than existing DEEP2/DEEP3 spectroscopy. This depth enables the first spectroscopic measurement of the satellite quiescent fraction down to stellar masses of $\sim 10^{9.5}~{\rm M}_{\odot}$ at this redshift. By combining an infall-based environmental quenching model, constrained by the observed quiescent fractions, with infall histories of simulated groups from the IllustrisTNG100-1-Dark simulation, we estimate environmental quenching timescales ($\tau_{\mathrm{quench}}$) for the observed group population. At high stellar masses (${M}_{\star}=10^{10.5}~{\rm M}_{\odot}$) we find that $\tau_{\mathrm{quench}} = 2.4\substack{+0.2 \\ -0.2}$ Gyr, which is consistent with previous estimates at this epoch. At lower stellar masses (${M}_{\star}=10^{9.5}~{\rm M}_{\odot}$), we find that $\tau_{\mathrm{quench}}=3.1\substack{+0.5 \\ -0.4}$ Gyr, which is shorter than prior estimates from photometry-based investigations. These timescales are consistent with satellite quenching via starvation, provided the hot gas envelope of infalling satellites is not stripped away. We find that the evolution in the quenching timescale between $0 \lt z \lt 1$ aligns with the evolution in the dynamical time of the host halo and the total cold gas depletion time. This suggests that the doubling of the quenching timescale in groups since $z\sim1$ could be related to the dynamical evolution of groups or a decrease in quenching efficiency via starvation with decreasing redshift.

The origins of Uranus and Neptune are not fully understood. Their inclined rotation axes -- obliquities -- suggest that they experienced giant impacts during their formation histories. Simulations modeling their accretion from giant impacts among ~5 Earth masses planetary embryos -- with roughly unity impactors' mass ratios -- have been able to broadly match their current masses, final mass ratio, and obliquity. However, due to angular momentum conservation, planets produced in these impacts tend to rotate too fast, compared to Uranus and Neptune. One potential solution for this problem consists of invoking instead collisions of objects with large mass ratios (e.g. a proto-Uranus with 13 Mearth and an embryo of 1 Mearth). Smooth-particle hydrodynamics simulations show that in this scenario final planets tend to have rotation periods more consistent with those of Uranus and Neptune. Here we performed a large suite of N-body numerical simulations modelling the formation of Uranus and Neptune to compare these different dynamical views. Our simulations start with a population of protoplanets and account for the effects of type-I migration, inclination and eccentricity tidal damping. Our results show that although scenarios allowing for large impactors' mass ratio favour slower rotating planets, the probability of occurring collisions in these specific simulations is significantly low. This is because gas tidal damping is relatively less efficient for low-mass embryos (<~1 Merath) and, consequently, such objects are mostly scattered by more massive objects (~13 Mearth) instead of colliding with them. Altogether, our results show that the probability of broadly matching the masses, mass ratio, and rotation periods of Uranus and Neptune in these two competing formation scenarios is broadly similar, within a factor of ~2, with overall probabilities of the order of ~0.1-1%.

Ultraviolet absorption line spectroscopy is a sensitive diagnostic for the properties of interstellar and circumgalactic gas. Down-the-barrel observations, where the absorption is measured against the galaxy itself, are commonly used to study feedback from galactic outflows and to make predictions about the leakage of HI ionizing photons into the intergalactic medium. Nonetheless, the interpretation of these observations is challenging and observational compromises are often made in terms of signal-to-noise, spectral resolution, or the use of stacking analyses. In this paper, we present a novel quantitative assessment of UV absorption line measurement techniques by using mock observations of a hydrodynamical simulation. We use a simulated galaxy to create 22,500 spectra in the commonly used SiII lines while also modeling the signal-to-noise and spectral resolution of recent rest-frame UV galaxy surveys at both high and low redshifts. We show that the residual flux of absorption features is easily overestimated for single line measurements and for stacked spectra. Additionally, we explore the robustness of the partial covering model for estimating column densities from spectra and find under-predictions on average of 1.25 dex. We show that the under-prediction is likely caused by high-column-density sight-lines that are optically-thick to dust making them invisible in UV spectra.

Galaxy mergers have been shown to trigger AGN in the nearby universe, but the timescale over which this process happens remains unconstrained. The Multi-Model Merger Identifier (MUMMI) machine vision pipeline has been demonstrated to provide reliable predictions of time post-merger (T_PM) for galaxies selected from the Ultraviolet Near Infrared and Optical Northern Survey (UNIONS) up to T_PM=1.76 Gyr after coalescence. By combining the post-mergers identified in UNIONS with pre-coalescence galaxy pairs, we can study the triggering of AGN throughout the merger sequence. AGN are identified using a range of complementary metrics: mid-IR colours, narrow emission lines and broad emission lines, which can be combined to provide insight into the demographics of dust and luminosity of the AGN population. Our main results are: 1) Regardless of the metric used, we find that the peak AGN excess (compared with a matched control sample) occurs immediately after coalescence, at 0 < T_PM < 0.16 Gyr. 2) The excess of AGN is observed until long after coalescence; both the mid-IR selected AGN and broad line AGN are more common than in the control sample even in the longest time bin of our sample (0.96 < T_PM < 1.76 Gyr). 3) The AGN excess is larger for more luminous and bolometrically dominant AGN, and we find that AGN in post-mergers are generally more luminous than secularly triggered events. 4) A deficit of broad line AGN in the pre-merger phase, that evolves into an excess in post-mergers is consistent with evolution of the covering fraction of nuclear obscuring material. Before coalescence, tidally triggered inflows increase the covering fraction of nuclear dust; in the post-merger regime feedback from the AGN clears (at least some of) this material. 5) The statistical peak in the triggering of starbursts occurs contemporaneously with AGN, within 0.16 Gyr of coalescence.

The sensitivity and field of view of the MeerKAT radio telescope provides excellent opportunities for commensal transient searches. We carry out a commensal transient search in supernova and short gamma-ray burst fields using methodologies established in Chastain et al. (2023). We search for transients in MeerKAT L-band images with 30 minute integration times, finding 13 variable sources. We compare these sources to the VLASS and RACS survey data, and examine possible explanation for the variability. We find that 12 of these sources are consistent with variability due to interstellar scintillation. The remaining source could possibly have some intrinsic variability. We also split the MeerKAT L-band into an upper and lower half, and search for transients in images with an 8 second integration time. We find a source with a duration of 8 to 16 seconds that is highly polarized at the lowest frequencies. We conclude that this source may be consistent with a stellar flare. Finally, we calculate accurate upper and lower limits on the transient rate using transient simulations.

Observations of young star clusters in a variety of galaxies have been used to constrain basic properties related to star-formation, such as the fraction of stars found in clusters (Gam) and the shape of the cluster mass function. However, the results can depend heavily on the reliability of the cluster age-dating process and other assumptions. One of the biggest challenges for successful age-dating lies in breaking the age-reddening degeneracy, where older, dust-free clusters and young, reddened clusters can have similar broad-band colors. While this degeneracy affects cluster populations in all galaxies, it is particularly challenging in dusty, extreme star-forming environments systems. We study the cluster demographics in the luminous infrared galaxy NGC1614 using Hubble imaging taken in 8 optical-NIR passbands. For age-dating, we adopt a spectral energy distribution fitting process that limits the maximum allowed reddening by region, and includes Ha photometry directly. We find that without these assumptions, essentially all clusters in the dust-free UV-bright arm which should have ages 50-250Myr are incorrectly assigned ages younger than 10Myr. We find this method greatly reduces the number of clusters in the youngest (tau<10Myrs) age bin and shows a fairly uniform distribution of massive clusters, the most massive being few10^7M. A maximum likelihood fit shows that the cluster mass function is well fitted by a power-law with an index -1.8, with no statistically significant high-mass cutoff. We calculate the fraction of stars born in clusters to be Gam1-10=22.4+_5.7%. The fraction of stars in clusters decreases quickly over time, with Gam10-100= 4.5+_1.1% and Gam100-400=1.7+_0.4%, suggesting that clusters dissolve rapidly over the first ~0.5Gyr. The decreasing fraction of stars in clusters is consistent with the declining shape observed for the cluster age distribution.

Exomoon discovery is on the horizon. Although several exomoon candidates exist around single stars, there are currently no candidates around circumbinary planets (CBPs). Circumbinary planets are predicted to have migrated to their current location from a farther region of the protoplanetary disc where they formed. An exomoon of a CBP therefore represents a fascinating yet complex and evolving four-body system. Their existence (or absence) would shed light on the robustness of moon formation in more dynamically active planetary systems. In this work, we simulate exomoons around migrating CBPs. We show that for fully migrated CBPs a moon is capable of surviving the migration if it is formed within ~5-10% of the planet's Hill Radius. Roughly 30-40% of the surviving moons are in the habitable zone, giving credence to circumbinary habitability even if the known CBPs all gas giants. The majority of moons fall off of their host planet early in the migration and become long-period CBPs (i.e a multi-planet circumbinary system). A subset of exomoons are ejected from the system entirely. This last class presents a new mechanism for producing free-floating planetary mass objects, like those discovered recently and expected in the Roman microlensing survey.

The presence of dips in the gravity modes period spacing versus period diagram of gamma Doradus stars is now well established by recent asteroseismic studies. Such Lorentzian-shaped inertial dips arise from the interaction of gravito-inertial modes in the radiative envelope of intermediate-mass main sequence stars with pure inertial modes in their convective core. They allow to study stellar internal properties. This window on stellar internal dynamics is extremely valuable in the context of the understanding of angular momentum transport inside stars, since it allows us to probe rotation in their core. We investigate the signature and the detectability of a differential rotation between the convective core and the near-core region inside gamma Doradus stars from the inertial dip properties. We study the coupling between gravito-inertial modes in the radiative zone and pure inertial modes in the convective core in the sub-inertial regime, allowing for a two-zones differential rotation from the two sides of the boundary. Taking a bi-layer rotation profile, we derive the wave equations in the convective core and the radiative envelope. We solve the coupling equation numerically and match the result to an analytical derivation of the Lorentzian dip. We then use typical values of measured near-core rotation and buoyancy travel time to infer ranges of parameters for which differential core to near-core rotation would be detectable in current Kepler data. We show that increasing the convective core rotation with respect to the near-core rotation leads to a shift of the period of the observed dip to lower periods. In addition, the dip gets deeper and thinner as the convective core rotation increases. We demonstrate that such a signature is detectable in Kepler data, given appropriate dip parameter ranges and near-core structural properties.

The binary pulsar J0348+0432 was previously shown to have a mass of approximately 2\,${\rm M_\odot}$, based on the combination of radial-velocity and model-dependent mass parameters derived from high-resolution optical spectroscopy of its white-dwarf companion. We present follow-up timing observations that combine archival observations with data acquired by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) pulsar instrument. We find that the inclusion of CHIME/Pulsar data yields an improved measurement of general-relativistic orbital decay in the system that falls within 1.2 $\sigma$ of the original values published by Antoniadis et al. (2013) while being roughly 6 times more precise due to the extended baseline. When we combine this new orbital evolution rate with the mass ratio determined from optical spectroscopy, we determine a pulsar mass of 1.806(37)\,${\rm M_\odot}$. For the first time for this pulsar, timing alone significantly constrains the pulsar mass. We explain why the new mass for the pulsar is $10\%$ lower and discuss how the mis-modeling of the initial observations of the white dwarf companion likely led to an inaccurate determination of the pulsar mass.

Recent multi-wavelength observations have highlighted magnetars as significant sources of cosmic rays, particularly through their gamma-ray emissions. This study examines three magnetar regions - CXOU J171405.7-31031, Swift J1834-0846, and SGR 1806-20 - known for emitting detectable electromagnetic signals. We assess the detectability of these regions using the upcoming Cherenkov Telescope Array Observatory (CTAO) by conducting an ON/OFF spectral analysis and compare the expected results with existing observations. Our findings indicate that CTAO will detect gamma-ray emissions from these three magnetar regions with significantly reduced emission flux errors compared to current instruments. In special, the study shows that the CXOUJ1714-3810 and SwiftJ1834-0846 magnetar regions can be observed by the full southern and northern CTAO arrays in just five hours of observation, with mean significances above $10 \,\sigma$ and $30 \,\sigma$, respectively. This paper discusses the regions analyzed, presents key results, and concludes with insights drawn from the study.

We report integral field spectroscopy observations with the Near-Infrared Spectrograph on board JWST targeting the 60 kpc environment surrounding the most luminous quasar known at $z=4.6$. We detect ionized gas filaments on 40 kpc scales connecting a network of merging galaxies likely to form a cluster. We find regions of low ionization consistent with large-scale shock excitation surrounding the central dust-obscured quasar, out to distances nearly eight times the effective stellar radius of the quasar host galaxy. In the nuclear region, we find an ionized outflow driven by the quasar with velocities reaching 13,000 km s$^{-1}$, one of the fastest discovered to date with an outflow rate of 2000 M$_\odot$ yr$^{-1}$ and a kinetic luminosity of 6$\times10^{46}$ erg s$^{-1}$ resulting in coupling efficiency between the bolometric luminosity of the quasar and the outflow of 5%. The kinetic luminosity of the outflow is sufficient to power the turbulent motion of the gas on galactic and circumgalactic scales and is likely the primary driver of the radiative shocks on interstellar medium and circumgalactic medium scales. This provides compelling evidence supporting long-standing theoretical predictions that powerful quasar outflows are a main driver in regulating the heating and accretion rate of gas onto massive central cluster galaxies.

Astrophysical evidence suggests that the Sun was born near 5 kpc from the Galactic center, within the corotation radius of the Galactic bar, around 6-7 kpc. This presents challenges for outward migration due to the Jacobi energy constraint, preventing stars from easily overcoming the corotation barrier. In this study, we use test particle simulations to explore two possible migration pathways for the Sun: a "trapped" scenario, where the Sun's orbit was influenced by a slowing Galactic bar, and an "untrapped" scenario driven by dynamic spiral arms. Our results demonstrate that both mechanisms can explain how the Sun migrated from its birth radius (approximately 5 kpc) to its current orbital radius around 8.5-9 kpc. Furthermore, we investigate the environmental changes experienced by the Sun along these migration pathways, focusing on variations in radiation hazards and comet fluxes, which may have impacted planetary habitability. These findings highlight the dynamic nature of galactic habitability, emphasizing that the path a star takes within the Milky Way can significantly affect its surrounding environment and the potential for life. We propose a new concept of "Galactic habitable orbits," which accounts for evolving galactic structures and their effects on stellar and planetary systems. This work contributes to a deeper understanding of the solar system's migration and its implications for habitability within the Milky Way.

We compiled the X-ray and soft gamma-ray observations of the Galactic black hole binary XTE J1859+226 in the 1999--2000 outburst from RXTE, ASCA, BeppoSAX and CGRO. Throughout systematic spectral analysis using a two-component model consisting of a multi-temperature accretion disk plus a fraction of its flux convolved with an empirical Comptonized powerlaw component, we found that the innermost radius ($r_{\rm in}$) and temperature (Tin) of the disk are very variable with time in the rising phase of soft X-ray flux where Type-A/B/C low-frequency quasi-periodic oscillations (QPOs) were found. After this phase, $r_{\rm in}$ remains constant at around 60 km assuming a distance of 8 kpc and an inclination angle of 67$^{\circ}$, and Tin smoothly decays with time. The constant $r_{\rm in}$ suggests a presence of the innermost stable circular orbit (ISCO), with $r_{\rm in}$ repeatedly moving closer and farther away from the ISCO in the rising phase. Both disk parameters are remarkably correlated with independently analyzed timing properties such as QPO frequency and rms variability. Type-A/B QPOs are seen only when $r_{\rm in}$ is close to the ISCO, while Type-C are seen when $r_{\rm in}$ is truncated and the frequency changes with a relation of $r^{-1.0}_{\rm in}$, supporting that Type-C QPOs occur at the inner edge of the truncated disk. Accurate determinations of the frequency--$r_{\rm in}$ relation for various objects should be a powerful tool to discriminate plausible Type-C QPO models. Furthermore, we suggest that jet ejection events may occur when $r_{\rm in}$ rapidly approaches to the ISCO, along with rapid changes of the disk flux, the rms variability and the hardness ratio. A rapid shrinkage of $r_{\rm in}$ down to the ISCO can be a useful index as a precursor of radio flares for triggering Target-of-Opportunity observations and would provide constraints on jet launching mechanisms.

Atomic hydrogen (HI) is a vital player in the star-formation process in galaxies. It is the raw fuel for making molecules, an important shielding agent against interstellar radiation, and a buffer that soaks up the energy and momentum of stellar feedback. While for many years detailed studies of the HI thermal structure have been possible only in the Milky Way, the SKA pathfinders are expanding our view beyond the Solar neighborhood allowing for crucial tests of the HI heating and cooling processes under a wide range of physical conditions. This overview article highlights a few recent results and emphasizes areas where future observations can make large contributions. We show that the cold HI disk of the Milky Way is extended and flared, yet appears spatially coupled to the molecular gas. The cold neutral medium (CNM) in the Milky Way is colder and more abundant at higher optical extinctions due to more intense cooling and shielding. A comparison between the Milky Way, the Small Magellanic Cloud, the Large Magellanic Cloud and NGC 6822 shows a good agreement with predictions from recent numerical simulations on how the CNM fraction depends on metalliciity. The fraction and spatial distribution of the thermally unstable HI remain as open questions and observationally have not been studied beyond the Milky Way. The excitation temperature of the warm neutral medium (WNM) is not well understood. Observations suggest higher WNM temperatures than what is seen in numerical simulations. The SKA, once operational, will take the HI studies into a new era of dense HI absorption grids for the Milky Way and nearby galaxies, and the ability to study the HI thermal structure beyond the Local Group.

Different waveform models can yield notably different conclusions about the properties of individual gravitational wave events. For instance, previous analyses using the SEOBNRv4PHM, IMRPhenomXPHM models, and NRSur7dq4 have led to varying results regarding event properties. This variability complicates the interpretation of the data and understanding of the astrophysical phenomena involved. There is an ongoing need to reassess candidate events with the best available interpretations and models. Current approaches lack efficiency or consistency, making it challenging to perform large-scale reanalyses with updated models or improved techniques. It is imperative that investigations into waveform systematics be reproducible. Frameworks like asimov can facilitate large-scale reanalyses with consistent settings and high-quality results, and can reliably show how different waveform models affect the interpretation of gravitational wave events. This, in combination with other provided tools, allow for reanalysis of several events from the GWTC-3 catalog. We include access to full analysis settings that facilitate public use of GWOSC data on the Open Science Grid, particularly those conducted with the IMRPhenomPv2, SEOBNRv4PHM, SEOBNRv5PHM, and NRSur7dq4 waveform models. Our parameter inference results find similar conclusions to previously published work: for several events, all models largely agree, but for a few exceptional events these models disagree substantially on the nature of the merging binary.

Works on the energy extraction from a rotating black hole via magnetic reconnection attract more attentions in recent years. Discussions on this topic, however, are often based on many simplifications, such as assuming a circularly flowing bulk plasma and a fixed orientation angle. A significant gap remains between theoretical models and the magnetic reconnection occurring in real astrophysical scenarios. In our previous work, we investigated the influence of orientation angle on energy extraction and figured out the differences between the plunging and circularly flowing bulk plasma. We introduced the concept of covering factor to quantify the capability of an accretion system in extracting energy via magnetic reconnection from a rotating black hole. In this work, as an improvement, we extend our discussions by treating the parameters in reconnection models as free parameters, bringing the theoretical model closer to the situations in real astrophysical systems. We employ two reconnection models, in which the geometric index and the guide field fraction are respectively induced. We separately present the dependences of energy extraction on the geometric index and the guide field fraction. More importantly, we propose to define the averaged covering factor weighted by reconnection rate to quantify the overall capability of an accretion system in extracting energy via magnetic reconnection from a rotating black hole, under the assumption that the magnetic reconnection process occurs stochastically and randomly within the ergosphere.

NinjaSat is an X-ray CubeSat designed for agile, long-term continuous observations of bright X-ray sources, with the size of 6U ($100\times200\times300$ mm$^3$) and a mass of 8 kg. NinjaSat is capable of pointing at X-ray sources with an accuracy of less than $0^{\circ}\hspace{-1.0mm}.1$ (2$\sigma$ confidence level) with 3-axis attitude control. The satellite bus is a commercially available NanoAvionics M6P, equipped with two non-imaging gas X-ray detectors covering an energy range of 2-50 keV. A total effective area of 32 cm$^2$ at 6 keV is capable of observing X-ray sources with a flux of approximately 10$^{-10}$ erg cm$^{-2}$ s$^{-1}$. The arrival time of each photon can be tagged with a time resolution of 61 $\mu$s. The two radiation belt monitors continuously measure the fluxes of protons above 5 MeV and electrons above 200 keV trapped in the geomagnetic field, alerting the X-ray detectors when the flux exceeds a threshold. The NinjaSat project started in 2020. Fabrication of the scientific payloads was completed in August 2022, and satellite integration and tests were completed in July 2023. NinjaSat was launched into a Sun-synchronous polar orbit at an altitude of about 530 km on 2023 November 11 by the SpaceX Transporter-9 mission. After about three months of satellite commissioning and payload verification, we observed the Crab Nebula on February 9, 2024, and successfully detected the 33.8262 ms pulsation from the neutron star. With this observation, NinjaSat met the minimum success criterion and stepped forward to scientific observations as initially planned. By the end of November 2024, we successfully observed 21 X-ray sources using NinjaSat. This achievement demonstrates that, with careful target selection, we can conduct scientific observations effectively using CubeSats, contributing to time-domain astronomy.

Density estimation is a fundamental problem that arises in many areas of astronomy, with applications ranging from selecting quasars using color distributions to characterizing stellar abundances. Astronomical observations are inevitably noisy; however, the density of a noise-free feature is often the desired outcome. The extreme-deconvolution (XD) method can be used to deconvolve the noise and obtain noise-free density estimates by fitting a mixture of Gaussians to data where each sample has non-identical (heteroscedastic) Gaussian noise. However, XD does not generalize to cases where some feature dimensions have highly non-Gaussian distribution, and no established method exists to overcome this limitation. We introduce a possible solution using neural networks to perform Gaussian mixture modeling of the Gaussian-like dimensions conditioned on those non-Gaussian features. The result is the CondXD algorithm, a generalization of XD that performs noise-free conditional density estimation. We apply CondXD to a toy model and find that it is more accurate than other approaches. We further test our method on a real-world high redshift quasar versus contaminant classification problem. Specifically, we estimate noise-free densities in flux-ratio (i.e., color) space for contaminants, conditioned on their magnitude. Our results are comparable to the existing method, which divides the samples into magnitude bins and applies XD separately in each bin, and our method is approximately ten times faster. Overall, our method has the potential to significantly improve estimating conditional densities and enable new discoveries in astronomy.

Finding pulsars in binaries are important for measurements of the masses of neutron stars, for tests of gravity theories, and for studies of star evolution. We are carrying out the Galactic Plane Pulsar Snapshot survey (GPPS) by using the the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Here we present the Keplerian parameters for 112 newly discovered pulsars in the FAST GPPS survey, and obtain timing solutions for 27 pulsars. Companions of these pulsars are He white dwarfs, CO/ONe white dwarfs, neutron stars, main sequence stars and ultra light objects or even planets. Our observations uncover eclipses of 8 binary systems. The optical counterpart for the companion of PSR J1908+1036 is identified. The Post-Keplerian parameter \Dot{\omega} for the double neutron star systems PSR J0528+3529 and J1844-0128 have been determined, with which the total masses of the binary systems are determined.

Binary millisecond pulsars with a massive white dwarf (WD) companion are intermediate-mass binary pulsars (IMBPs). They are formed via the Case BB Roche-lobe overflow (RLO) evolution channel if they are in compact orbits with an orbital period of less than 1 day. They are fairly rare in the known pulsar population, only five such IMBPs have been discovered before, and one of them is in a globular cluster. Here we report six IMBPs in a compact orbit, PSRs J0416+5201, J0520+3722, J1919+1341, J1943+2210, J1947+2304 and J2023+2853, discovered during the Galactic Plane Pulsar Snapshot (GPPS) survey by using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), doubling the number of such IMBPs due to the high survey sensitivity in the short survey time of 5 minutes. Follow-up timing observations show that they all have either a CO WD or an ONeMg WD companion with a mass greater than about 0.8~M$_\odot$ in a very circular orbit with an eccentricity in the order of $\lesssim10^{-5}$. PSR J0416+5201 should be an ONeMg WD companion with a remarkable minimum mass of 1.28 M$_\odot$. These massive white dwarf companions lead to a detectable Shapiro delay for PSRs J0416+5201, J0520+3722, J1943+2210, and J2023+2853, indicating that their orbits are highly inclined. From the measurement of the Shapiro delay, the pulsar mass of J1943+2210 was constrained to be 1.84$^{\,+0.11}_{-0.09}$~M$_\odot$, and that of PSR J2023+2853 to be 1.28$^{\,+0.06}_{-0.05}$~M$_\odot$.

The ability to differentiate between different models of inflation through the imprint of primordial non-Gaussianity (PNG) requires a tight constraint on the local PNG parameter $f_{\text{NL}}^{\text{loc}}$. Future large scale structure surveys like \textit{Euclid}, Vera C. Rubin Observatory, and the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer (SPHEREx) will play a crucial role in advancing our understanding of the inflationary epoch. In this light, we present forecasts on PNG with tomographic angular power spectrum from simulations of SPHEREx. We put forward the effects of redshift bin mismatch of galaxies as a source of systematic error when estimating $f_{\text{NL}}^{\text{loc}}$ and galaxy linear halo bias. We %use \texttt{GLASS} to simulate $500$ SPHEREx-like galaxy density fields, and divide the galaxies into $13$ redshift bins assuming Gaussian photometric redshift errors. We show that the misclassification of galaxies in redshift bins can result in strong apparent tensions on $f_{\text{NL}}^{\text{loc}}$ up to $\sim 3\sigma$ and up to $\sim 9\sigma$ on galaxy bias. We propose a scattering matrix formalism to mitigate bin mismatch of galaxies and to obtain unbiased estimates of cosmological parameters from tomographic angular clustering measurements.

The dynamics of the merger of a dwarf disc galaxy with a massive spiral galaxy of the Milky Way type have been studied in detail. The remnant of such interaction after numerous crossings of the satellite through the disc of the main galaxy is a compact stellar core, the characteristics of which are close to small compact elliptical galaxies (cEs) or large ultra-compact dwarfs (UCDs). Such transitional cE/UCD objects with an effective radius of 100-200 pc arise as a result of stripping the outer layers of the stellar core during the destruction of a disc dwarf galaxy. Numerical models of the satellite before interaction include baryonic matter (stars and gas) and dark mass. We use N-body to describe the dynamics of stars and dark matter and Smoothed-Particle Hydrodynamics to model the gas components of both galaxies. The direct method of calculating the gravitational force between all particles provides a qualitative resolution of spatial structures up to 10 pc. The simulated cE/UCD galaxies contain very little gas and dark matter at the end of their evolution.

Nuclear Star Clusters (NSCs) are commonly found in galaxy centers, but their dominant formation mechanisms remain elusive. We perform a consistent analysis of stellar populations of 97 nearby NSCs, based on VLT spectroscopic data. The sample covers a galaxy stellar mass range of 10$^{7}$ to 10$^{11}$ M$_{\odot}$ and is more than 3$\times$ larger than any previous studies. We identify three galaxy stellar mass regimes with distinct NSC properties. In the low-mass regime of $\log M_{\rm host}$ $\lesssim$ 8.5, nearly all NSCs have metallicities lower than circum-NSC host but similar to typical red globular clusters (GCs), supporting the GC inspiral-merger scenario of NSC formation. In the high-mass regime of $\log M_{\rm host}$ $\gtrsim$ 9.5, nearly all NSCs have higher metallicities than circum-NSC host and red GCs, suggesting significant contributions from in-situ star formation (SF). In the intermediate-mass regime, a comparable fraction of NSCs have higher or lower metallicities than circum-NSC host and red GCs, with no clear dependence on NSC mass, suggesting intermittent in-situ SF. The majority of NSCs with higher metallicities than their host exhibit a negative age$-$metallicity correlation, providing clear evidence of long-term chemical enrichment. The average NSC$-$host metallicity difference peaks broadly around $\log M_{\rm host} \sim 9.8$ and declines towards both higher and lower galaxy mass. We find that the efficiency of dynamical friction-driven inspiral of GCs observed in present-day galaxies can explain the NSC mass at $\log M_{\rm host} \lesssim 9.5$ but falls short of observed ones at higher galaxy mass, reinforcing our conclusions based on stellar population analysis.

Impact flashes on the moon are caused by high-speed collisions of celestial bodies with the lunar surface. The study of the impacts is critical for exploring the evolutionary history and formation of the Moon, and for quantifying the risk posed by the impacts to future human activity. Although the impacts have been monitored from the Earth by a few projects in past 20 years, the events occurring on the lunar far side have not been explored systematically so far. We here present an end-to-end image simulator dedicated to detecting and monitoring the impacts from space, which is useful for future mission design. The simulator is designed for modularity and developed in the Python environment, which is mainly composed of four components: the flash temporal radiation, the background emission, the telescope and the detector used to collect and measure the radiation. Briefly speaking, with a set of input parameters, the simulator calculates the flash radiation in the context of the spherical droplet model and the background emission from the lunar surface. The resulting images are then generated by the simulator after considering a series observational effects, including the stray light, transmission of the instrument, point spread function and multiple kinds of noise caused by a CCD/CMOS detector. The simulator is validated by comparing the calculation with the observations taken on the ground. The modular design enables the simulator to be improved and enhanced by including more complex physical models in the future, and to be flexible for other future space missions.

The broad-band radio spectral energy distribution (SED) of star-forming galaxies (SFGs) contains a wealth of complex physics. We aim to determine the physical emission and loss processes causing radio SED curvature and steepening to see which observed global astrophysical properties are correlated with radio SED complexity. We have acquired radio continuum data between 70 MHz and 17 GHz for a sample of 19 southern local (z < 0.04) SFGs. Of this sample 11 are selected to contain low-frequency (< 300 MHz) turnovers (LFTOs) in their SEDs and eight are control galaxies with similar global properties. We model the radio SEDs for our sample using a Bayesian framework whereby radio emission (synchrotron and free-free) and absorption or loss processes are included modularly. We find that without the inclusion of higher frequency data, single synchrotron power-law based models are always preferred for our sample; however, additional processes including free-free absorption (FFA) and synchrotron losses are often required to accurately model radio SED complexity in SFGs. The fitted synchrotron spectral indices range from -0.45 to -1.07 and are strongly anticorrelated with stellar mass suggesting that synchrotron losses are the dominant mechanism acting to steepen the spectral index in larger nearby SFGs. We find that LFTOs in the radio SED are independent from the inclination. The merging systems in our SFG sample have elevated specific star formation rates and flatter fitted spectral indices with unconstrained LFTOs. Lastly, we find no significant separation in global properties between SFGs with or without modelled LFTOs. Overall LFTOs are likely caused by a combination of FFA and ionisation losses in individual recent starburst regions with specific orientations and interstellar medium properties that, when averaged over the entire galaxy, do not correlate with global astrophysical properties.

We revisit the binary and stellar properties of the double-degenerate system NLTT 16249. An analysis of new echelle spectra, supported by a joint study of a DQZ velocity template NLTT 44303, confirms the orbital period and constrains the mass ratio revealing a carbon-polluted DQ white dwarf that is up to ~6 percent more massive than its hydrogen-rich DA companion. Our new model atmosphere analysis of the DA and DQ components, constrained by an accurate Gaia parallax measurement that places the binary at a distance of 57.8 pc, reveals lower mass and temperature than previously estimated for both components, but with higher carbon and nitrogen abundances in the DQ atmosphere. The two components are nearly coeval and could have been generated following a single common envelope event.

Low-frequency radio emission from the Large Magellanic Cloud~(LMC) is assumed to be dominated by nonthermal synchrotron radiation from energy loss of energetic $e^+/e^-$ in magnetic field. Two different kinds of sources of $e^+/e^-$, dark matter~(DM) annihilation and cosmic rays~(CR) related to massive stars, are taken into account in this paper. We fit the multiple low-frequency radio observations, from 19.7 MHz to 1.4 GHz, with a double power-law model $S_{nth} =S_{DM}(\frac {\nu}{\nu_{\star}})^{-\alpha_{DM}}+S_{CR}( \frac {\nu}{\nu_{\star}})^{-\alpha_{CR}} $. $\nu_{\star}$ is set to be $1.4$ GHz and $S_{CR}$ could be determined from the 24 $\mu m$ luminosity based on the global radio-infrared correlation. Our best fit with a fixed $\alpha_{CR}$ changing from $0.80$ to $0.55$ yields $\alpha_{DM}$ ranging from $0.21$ to $0.66$. Given a fixed value of $\alpha_{CR}$, we derive the upper limits of synchrotron emission induced by dark matter annihilation at different radio frequencies. Larger value of $\alpha_{CR}$ represents for a harder $e^+/e^-$ spectrum from cosmic rays, which leads to a smaller value of $\alpha_{DM}$ and allow less synchrotron emission resulted from dark matter annihilation in lower frequency. Under the same assumption on the magnetic field, we find that the lower the frequency, the stronger the restriction on DM parameter space. Meanwhile, as the peak frequency of synchrotron radiation decrease with the energy of $e^+/e^-$, constraints on DM properties obtained from lower frequency are more severe in the case of DM with lower mass. Future low-frequency radio survey should be considered a promising and powerful way to constrain DM.

We reformulate the gain correction problem of the radio interferometry as an optimization problem with regularization, which is solved efficiently with an iterative algorithm. Combining this new method with our previously proposed imaging method, PRIISM, the whole process of the self-calibration of radio interferometry is redefined as a single optimization problem with regularization. As a result, the gains are corrected, and an image is estimated. We tested the new approach with ALMA observation data and found it provides promising results.

Quiet-Sun Ellerman bombs (QSEBs) are small-scale magnetic reconnection events in the lower atmosphere of the quiet Sun. Recent work has shown that a small percentage of them can occur co-spatially and co-temporally to ultraviolet (UV) brightenings in the transition region. We aim to understand how the magnetic topologies associated with closely occurring QSEBs and UV brightenings can facilitate energy transport and connect these events. We used high-resolution H-beta observations from the Swedish 1-m Solar Telescope (SST) and detected QSEBs using k-means clustering. We obtained the magnetic field topology from potential field extrapolations using spectro-polarimetric data in the photospheric Fe I 6173 A line. To detect UV brightenings, we used coordinated and co-aligned data from the Interface Region Imaging Spectrograph (IRIS) and imposed a threshold of 5 sigma above the median background on the (IRIS) 1400 A slit-jaw image channel. We identify four distinct magnetic configurations that associate QSEBs with UV brightenings, including a simple dipole configuration and more complex fan-spine topologies with a three-dimensional (3D) magnetic null point. In the fan-spine topology, the UV brightenings occur near the 3D null point, while QSEBs can be found close to the footpoints of the outer spine, the inner spine, and the fan surface. We find that the height of the 3D null varies between 0.2 Mm to 2.6 Mm, depending on the magnetic field strength in the region. We note that some QSEBs and UV brightenings, though occurring close to each other, are not topologically connected with the same reconnection process. We find that the energy released during QSEBs falls in the range of 10^23 to 10^24 ergs. This study shows that magnetic connectivity and topological features, like 3D null points, are crucial in linking QSEBs in the lower atmosphere with UV brightenings in the transition region.

We use a simple dynamical scheme to simulate the ejecta of type Ia supernova (SN Ia) scenarios with two exploding white dwarfs (WDs) and find that the velocity distribution of the ejecta has difficulties accounting for bimodal emission line profiles with a large separation between the two emission peaks. The essence of the dynamical code is in including the fact that the ejecta does not leave the system instantaneously. We find that the final separation velocity between the centers of masses of the two WDs' ejecta is ~80% of the pre-explosion WDs' orbital velocity, i.e., we find separation velocities of 4200-5400 km/s for two WDs of masses M1=M2=0.94 Mo. The lower separation velocities we find challenge scenarios with two exploding WDs to explain bimodal emission line profiles with observed velocity separations of up to ~7000 km/s. Only the mass in the ejecta of one WD with an explosion velocity lower than the separation velocity contributes to one peak of the bimodal profile; this is the inner ejecta. We find the inner ejecta to be only <15% of the ejecta mass in energetic explosions. Less energetic explosions yield higher inner mass but lower separation velocities. We encourage searching for alternative explanations of bimodal line profiles.

We present the first multi-band centimeter detection of POX 52, a nearby dwarf galaxy believed to habor a robust intermediate mass black hole (IMBH). We conducted the deep observations using the Australia Telescope Compact Array (ATCA), spanning frequencies from 4.5 to 10 GHz, as well as the sensitive observations from the Karl G. Jansky Very Large Array (VLA) operating in its most extended A-configuration at S band (2--4 GHz) and C band (4--8 GHz). In the ATCA observations, the source shows a compact morphology, with only one direction marginally resolved. The higher resolution of the VLA allowed us to slightly resolve the source, fitting it well with a two-dimensional Gaussian model. The detected radio emission confirms the presence of Active Galactic Nucleus (AGN) activity, indicating either a low-power jet or AGN-driven winds/outflows. Our dual-epoch observations with ATCA and VLA, together with previous non-detection flux density upper limits, reveal radio emission variability spanning two decades. In addition, we find that POX 52 aligns well with the low-mass extension of the fundamental plane for high-accretion, radio-quiet massive AGNs.

Surveillance of objects in the cislunar domain is challenging due primarily to the large distances (10x the Geosynchronous orbit radius) and total volume of space to be covered. Ground-based electro-optical observations are further hindered by high background levels due to scattered moonlight. In this paper, we report on ground-based demonstrations of space surveillance for targets in the cislunar domain exploiting the remarkable performance of 36cm, F2.2 class telescopes equipped with current generation, back-side illuminated, full-frame CMOS imager. The demonstrations leverage advantageous viewing conditions for the Artemis Orion vehicle during its return to earth, and the total lunar eclipse of November 8th 2022 for viewing the CAPSTONE vehicle. Estimated g-band magnitudes for vehicles were 19.57 at a range of 4.4e5 km and 15.53 at a range of 3.2e5 km for CAPSTONE and Artemis Orion, respectively. In addition to observations, we present reflectance signature modeling implemented in the LunaTK space-sensing simulation framework and compare calculated apparent magnitudes to observed. Design of RApid Telescopes for Optical Response (RAPTOR) instruments, observing campaigns of Artemis Orion and CAPSTONE missions, and initial comparisons to electro-optical modeling are reviewed.

Recent observations suggest that the extended stellar halos of low-redshift massive galaxies are tightly connected to the assembly of their dark matter halos. In this paper, we use the Illustris, IllustrisTNG100, and IllustrisTNG300 simulations to compare how different stellar aperture masses trace halo mass. For massive central galaxies ($M_\star\geq 10^{11.2}M_\odot$), we find that a 2D outskirt stellar mass measured between 50 to 100 kpc ($M_{\star,[50,100]}$) consistently outperforms other aperture-based stellar masses. We further show that $M_{\star,[50,100]}$ correlates better with halo mass than the total amount of accreted stars (the ex situ mass), which suggests that not all accreted stars connect to halo assembly equally. While the galaxy formation recipes are different between Illustris and IllustrisTNG100, the two simulations yield consistent ex situ outskirt fractions for massive galaxies (about 70% in $M_{\star,[50,100]}$). These results demonstrate the potential of using the outskirt stellar mass to deepen our understanding of galaxy-halo connection in massive dark matter halos and trace dark matter halos better.

We report JWST NIRSpec/G395H observations of TOI-1685 b, a hot rocky super-Earth orbiting an M2.5V star, during a full orbit. We obtain transmission and emission spectra of the planet and characterize the properties of the phase curve, including its amplitude and offset. The transmission spectrum rules out clear H$_2$-dominated atmospheres, while secondary atmospheres (made of water, methane, or carbon dioxide) cannot be statistically distinguished from a flat line. The emission spectrum is featureless and consistent with a blackbody-like brightness temperature, helping rule out thick atmospheres with high mean molecular weight. Collecting all evidence, the properties of TOI-1685 b are consistent with a blackbody with no heat redistribution and a low albedo, with a dayside brightness temperature 0.98$\pm$0.07 times that of a perfect blackbody in the NIRSpec NRS2 wavelength range (3.823-5.172 um). Our results add to the growing number of seemingly airless M-star rocky planets, thus constraining the location of the "Cosmic Shoreline". Three independent data reductions have been carried out, all showing a high-amplitude correlated noise component in the white and spectroscopic light curves. The correlated noise properties are different between the NRS1 and NRS2 detectors - importantly the timescales of the strongest components (4.5 hours and 2.5 hours, respectively) - suggesting the noise is from instrumental rather than astrophysical origins. We encourage the community to look into the systematics of NIRSpec for long time-series observations.

This paper presents a comparative analysis of the bulk properties (mass and radius) of transiting giant planets ($\gtrsim$ 8$R_{\oplus}$) orbiting FGKM stars. Our findings suggest that the average mass of M-dwarf Jupiters is lower than that of their solar-type counterparts, primarily due to the scarcity of super-Jupiters ( $\gtrsim$ 2 $M_J$) around M-dwarfs. However, when super-Jupiters are excluded from the analysis, we observe a striking similarity in the average masses of M-dwarf and FGK warm-Jupiters. We propose that these trends can be explained by a minimum disk dust mass threshold required for Jovian formation through core accretion, which is likely to be satisfied more often around higher mass stars. This simplistic explanation suggests that the disk mass has more of an influence on giant planet formation than other factors such as the host star mass, formation location, metallicity, radiation environment, etc., and also accounts for the lower occurrence of giant planets around M-dwarf stars. Additionally, we explore the possibility of an abrupt transition in the ratio of super-Jupiters to Jupiters around F-type stars at the Kraft break, which could be a product of $v$sin$i$ related detection biases, but requires additional data from an unbiased sample with published non-detections to confirm. Overall, our results provide valuable insights into the formation and evolution of giant exoplanets across a diverse range of stellar environments.

This paper presents an investigation of the X-ray emission associated with the Wolf-Rayet star, WR 48-6, using observations from the XMM Newton and Chandra X-ray telescopes covering two epochs separated by eleven months. The X-ray spectrum of WR 48-6 is well explained by a two-temperature plasma model, with cool and hot plasma temperatures of $0.8_{-0.2}^{\,+0.1}$ and $2.86_{-0.66}^{\,+1.01}$ keV. No significant X-ray variability is observed during these two epochs of observations. However, an increase in the local hydrogen column density accompanied by a decrease in the intrinsic X-ray flux between two epochs of observations is seen. Additionally, the intrinsic X-ray luminosity is found to be more than $10^{33} \rm\,erg\,s^{-1}$ during both epochs of observations. Based on the analysis presented, WR 48-6 is a promising colliding wind binary candidate with a possible companion of spectral type O5-O6.

This paper presents 'SpyDust', an improved and extended implementation of the spinning dust emission model based on a Fokker-Planck treatment. 'SpyDust' serves not only as a Python successor to 'spdust', but also incorporates some corrections and extensions. Unlike 'spdust', which is focused on specific grain shapes, 'SpyDust' considers a wider range of grain shapes and provides the corresponding grain dynamics, directional radiation field and angular momentum transports. We recognise the unique effects of different grain shapes on emission, in particular the shape-dependent mapping between rotational frequency and spectral frequency. In addition, we update the expressions for effects of electrical dipole radiation back-reaction and plasma drag on angular momentum dissipation. We also discuss the degeneracies in describing the shape of the spectral energy distribution (SED) of spinning dust grains with the interstellar environmental parameters. Using a typical Cold Neutral Medium (CNM) environment as an example, we perform a perturbative analysis of the model parameters, revealing strong positive or negative correlations between them. A principal component analysis (PCA) shows that four dominant modes can linearly capture most of the SED variations, highlighting the degeneracy in the parameter space of the SED shape in the vicinity of the chosen CNM environment. This opens the possibility for future applications of moment expansion methods to reduce the dimensionality of the encountered SED parameter space.

We use the Millennium Simulation to study the relation of galaxies and dark matter haloes to the cosmic web. We define the web as the unique, fully connected, percolating object with (unsmoothed) matter density everywhere exceeding 5.25 times the cosmic mean. This object contains 35\% of all cosmic mass but occupies only 0.62\% of all cosmic volume. It contains 26\% of dark matter haloes of mass $10^{11}M_\odot$, rising to 50\% at $10^{12.7}M_\odot$, and to $>90\%$ above $10^{14}M_\odot$. In contrast, it contains 45\% of all galaxies of stellar mass $10^{8.5}M_\odot$, rising to 50\% at $10^{10}M_\odot$, to 60\% at $10^{11}M_\odot$ and to 90\% at $10^{11.5}M_\odot$. This difference arises because a large fraction of all satellite and backsplash galaxies are part of the cosmic web. Indeed, more than 50\% of web galaxies are satellites for stellar masses below that of the Milky Way, rising to about 70\% below $10^{10}M_\odot$, whereas centrals substantially outnumber satellites in the non-web population at all stellar masses. As a result, web galaxies have systematically lower specific star-formation rates (sSFR's) than non-web galaxies. For the latter, the distributions of stellar mass and sSFR are almost independent of web distance. Furthermore, for both central and satellite galaxies, the sSFR distributions at given stellar mass are identical in and outside the web, once differences in backsplash fraction are accounted for. For the galaxy formation model considered here, differences between web and non-web galaxy populations are almost entirely due to the difference in halo mass distribution between the two environments.

We study the rest-frame ultraviolet-optical color gradients of 669 galaxies at $4<z<8$ by characterizing the wavelength dependence of their structural parameters derived from simultaneously fitting the seven-band NIRCam images acquired with the James Webb Space Telescope. Distinct from trends observed at lower redshifts, where most galaxies exhibit negative color gradients whereby galaxy centers are redder than their outskirts, in high-redshift galaxies positive color gradients match or even outnumber negative color gradients. The color gradients principally reflect radial variations in stellar population instead of dust reddening or contamination from active galactic nuclei. The sign and magnitude of the color profile depend systematically on the global properties of the galaxy: positive color gradients, characteristic of centrally concentrated star formation or outside-in growth, preferentially inhabit galaxies of lower stellar mass, smaller size, and bluer spectral energy distribution.

In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high instrumental efficiency ( $>$ 37\%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R $>$ 20, 000 (with a lower-resolution, sky-limited mode of R $\sim$ 7, 000). With the design focusing on maximizing the instrument throughput (ensuring a Signal to Noise Ratio -SNR- $\sim$ 20 per spectral resolution element at 313 nm for U $\sim$ 17.5 mag objects in 1h of observations), it will offer new possibilities in many fields of astrophysics: i) access to key lines of stellar spectra (e.g. lighter elements, in particular Beryllium), extragalactic studies (e.g. circumgalactic medium of distant galaxies, cosmic UV background) and follow-up of explosive transients. We present the CUBES instrument design, currently in Phase-C and approaching the final design review, summarizing the hardware architecture and interfaces between the different subsystems as well as the relevant technical requirements. We describe the optical, mechanical, electrical design of the different subsystems (from the telescope adapter and support structure, through the main opto-mechanical path, including calibration unit, detector devices and cryostat control, main control electronics), detailing peculiar instrument functions like the Active Flexure Compensation (AFC). Furthermore, we outline the AITV concept and the main instrument operations giving an overview of its software ecosystem. Installation at the VLT is planned for 2028-2029 and first science operations in late 2029.

The origin of black hole (BH) spins remains one of the least understood aspects of BHs. Despite many uncertainties, it is commonly assumed that if BHs originated from isolated massive star binaries, their spins should be aligned with the orbital angular momentum of the binary system. This assumption stems from the notion that BHs inherit their spins from their progenitor stars. In this study, we relax this long-held viewpoint and explore various mechanisms that can spin up BHs before or during their formation. In addition to natal spins, we discuss physical processes that can spin BHs isotropically, parallel to natal kicks, and perpendicular to natal kicks. These different mechanisms leave behind distinct imprints on the observable distributions of spin magnitudes, spin-orbit misalignments and the effective inspiral spin of merging binaries. In particular, these mechanisms allow even the binaries originating in the field to exhibit precession and retrograde spin ($\chi_{\rm eff}<0$). This broadens the parameter space allowed for isolated binary evolution into regimes which were previously thought to be exclusive to dynamically assembled binaries.

The occurrence frequency distributions of fluxes (F) and fluences or energies (E) observed in astrophysical observations are found to be consistent with the predictions of the fractal-diffusive self-organized criticality (FD-SOC) model, which predicts power law slopes with universal constants of $\alpha_F=(9/5)=1.80$ for the flux and $\alpha_E=(5/3)\approx 1.67$ for the fluence. The energy integrated over the power law-like (size distribution) energy range is found to be finite for these power law slopes with $\alpha_E < 2$, which refutes earlier claims of a divergent energy integral that has been postulated in the energy budget of solar and stellar nanoflare scenarios. The theoretial FD-SOC model approximates the microscopic cellular automaton models satisfactorily with the macroscopic scaling law of classical diffusion. The universal scaling laws predict the size distributions of numerous astrophysical phenomena, such as solar flares, stellar flares, coronal mass ejections (CME), auroras, blazars, galactic fast radio bursts (FRB), active galactic nuclei (AGN), gamma-ray bursts (GRB), soft gamma-ray repeaters (SGB), and black-hole systems (BH), while coherent solar radio bursts, random radio bursts, solar energetic partices (SEP), cosmic rays, and pulsar glitches require non-standard SOC models.

Building fast and accurate ways to model the distribution of neutral hydrogen during the Epoch of Reionization (EoR) is essential for interpreting upcoming 21 cm observations. A key component of semi-numerical models of reionization is the collapse fraction field $f_{\text{coll}}(\mathbf{x})$, which represents the fraction of mass within dark matter halos at each location. Using high-dynamic range N-body simulations to obtain this is computationally prohibitive and semi-analytical approaches, while being fast, end up compromising on accuracy. In this work, we bridge the gap by developing a machine learning model that can generate $f_{\text{coll}}$ maps by sampling from the full distribution of $f_{\text{coll}}$ conditioned on the dark matter density contrast $\delta$. The conditional distribution functions and the input density field to the model are taken from low-dynamic range N-body simulations that are more efficient to run. We evaluate the performance of our ML model by comparing its predictions to a high-dynamic range N-body simulation. Using these $f_{\text{coll}}$ maps, we compute the HI and HII maps through a semi-numerical code for reionization. We are able to recover the large-scale HI density field power spectra $(k \lesssim 1\ h\,{\rm Mpc}^{-1})$ at the $\lesssim10\%$ level, while the HII density field is reproduced with errors well below 10% across all scales. Compared to existing semi-analytical prescriptions, our approach offers significantly improved accuracy in generating the collapse fraction field, providing a robust and efficient alternative for modeling reionization.

Black hole accretion in active galactic nuclei (AGN) is coupled to the evolution of their host galaxies. Outflowing winds in AGN can play an important role in this evolution through the resulting feedback mechanism. Multi-wavelength spectroscopy is key for probing the intertwined physics of inflows and outflows in AGN. However, with the current spectrometers, crucial properties of the ionized outflows are poorly understood, such as their coupling to the accretion rate, their launching mechanism, and their kinetic power. In this paper we discuss the need for simultaneous X-ray and UV high-resolution spectroscopy for tackling outstanding questions on these outflows in AGN. The instrumental requirements for achieving the scientific objectives are addressed. We demonstrate that these requirements would be facilitated by the proposed Arcus Probe mission concept. The multi-wavelength spectroscopy and timing by Arcus would enable us to establish the kinematics and ionization structure of the entire ionized outflow, extending from the vicinity of the accretion disk to the outskirts of the host galaxy. Arcus would provide key diagnostics on the origin, driving mechanism, and the energetics of the outflows, which are useful benchmarks for testing various theoretical models of outflows and understanding their impact in AGN.

Magnetic fields, together with cosmic rays (CRs), play an important role in the dynamics and evolution of galaxies, but are difficult to estimate. Energy equipartition between magnetic fields and CRs provides a convenient way to approximate magnetic field strength from radio observations. We present a new approach for calculating the equipartition magnetic field strength based on Bayesian methods. In this approach, the magnetic field is a random variable that is distributed according to a posterior distribution conditional on synchrotron emission and the size of the emitting region. It allows the direct application of the general formulas for total and polarized synchrotron radiation without the need to invert these formulas, which has limited the equipartition method to highly simplified cases. We have derived the equipartition condition for the case of different low-energy breaks, slopes, and high-energy cutoffs of power law spectra of the CR proton and electron distributions. The derived formalism was applied in the general case of a magnetic field consisting of both uniform and randomly oriented field components. The applied Bayesian approach naturally provides the uncertainties in the estimated magnetic field strengths resulting from the uncertainties in the observables and the assumed values of the unknown physical parameters. In the examples presented, we used two different Markov Chain Monte Carlo methods to generate the posterior distribution of the magnetic field. We have also developed a web application called BMAG that implements the described approach for different models and observational parameters of real sources.

In this study we are testing whether the power law slopes ($\alpha_F$, $\alpha_E$) of fluxes $(F)$, fluences or energies $(E)$ are universal in their size distributions, $N(F) \propto F^{-\alpha_F}$ and $N(E) \propto E^{-\alpha_E}$, in astrophysical observations of galactic, extragalactic, and black-hole systems. This is a test of fundamental importance for self-organized criticality (SOC) systems. The test decides whether (i) power laws are a natural consequence of the scale-freeness and inherent universality of SOC systems, or (ii) if they depend on more complex physical scaling laws. The former criterion allows quantitative predictions of the power law-like size distributions, while the later criterion requires individual physical modeling for each SOC variable and data set. Our statistical test, carried out with 61 published data sets, yields strong support for the former option, which implies that observed power laws can simply be derived from the scale-freeness and do not require specific physical models to understand their statistical distributions. The observations show a mean and standard deviation of $\alpha_F=1.78\pm0.29$ for SOC fluxes, and $\alpha_E=1.66\pm0.22$ for SOC fluences, and thus are consistent with the prediction of the fractal-diffusive SOC model, with $\alpha_F=1.80$ and $\alpha_E=1.67$.

Many physical models contain nuisance parameters that quantify unknown and irrelevant properties of an experiment. Typically, these cannot be measured except by fitting the models to the data from the experiment, requiring simultaneous measurement of interesting parameters that are our target of inference and nuisance terms that are not directly of interest. A recent example of this is fitting Effective Field Theory (EFT) models to large-scale structure (LSS) data to make cosmological inferences. These models have a large number of nuisance parameters that are typically correlated with cosmological parameters in the posterior, leading to strong dependence on the nuisance parameter priors. We introduce a reparametrization method that leverages Generalized Additive Models (GAMs) to decorrelate nuisance parameters from the parameters of interest in the likelihood, even in the presence of non-linear relationships. This reparametrization forms a natural basis within which to define priors that are independent between nuisance and target parameters: the separation means that the marginal posterior for cosmological parameters does not depend on simple priors placed on nuisance terms. In application to EFT models using LSS data, we demonstrate that the proposed approach leads to robust cosmological inference.

Red supergiants may lose significant mass through steady winds and episodic eruptions in the final 100-1000 years before the core collapses, shaping their circumstellar environment. Interaction between supernova (SN) ejecta and distant circumstellar material (CSM) can generate shocks, which can energize the ejecta and serve as a key power source during the nebular phase of the SN. In the present work, we investigate the nebular spectrum of SN 2023ixf, observed one year post-explosion (at +363 d) with the recently commissioned WEAVE instrument on the 4.2m William Herschel Telescope. This marks the first supernova spectrum captured with WEAVE. In this spectrum, H$\alpha$ exhibits a peculiar evolution, flanked by blueward and redward broad components centred at $\sim\pm 5650\,\mathrm{km\,s^{-1}}$ from the rest velocity of H$\alpha$, which are seen for only a few SNe to date. These features indicate energy deposition from shocks generated by the interaction of ejecta with a CSM expelled nearly 350 $-$ 640 years pre-explosion. Comparisons of the +363 d spectrum with model spectra from the literature, that include varying shock powers, suggest a shock power of at least $\sim 5 \times 10 ^{40}\,\mathrm{erg\,s^{-1}}$ at this epoch. Additionally, analysis of the [O I] doublet, along with other prominent emission lines, provides evidence for clumpiness, dust formation, and asymmetry within the ejecta and/or the surrounding CSM. These emission lines also helped to constrain the oxygen mass ($\approx0.19^{\scriptscriptstyle +0.08}_{\scriptscriptstyle -0.04} M_\odot$), He-core mass ($<3 M_\odot$) and the zero-age main sequence mass ($\lesssim 12 M_\odot$) of the progenitor of SN 2023ixf. The comparison with other Type II SNe highlights SN 2023ixf's unique shock interaction signatures and evidence of dust formation, setting it apart in terms of evolution and dynamics.

Energy dissipation in collisionless shocks is a key mechanism in various astrophysical environments. Its non-linear nature complicates analytical understanding and necessitate Particle-in-Cell (PIC) simulations. This study examines the impact of reducing the ion-to-electron mass ratio ($m_r$), to decrease computational cost, on energy partitioning in 1D3V (one spatial and three velocity-space dimensions) PIC simulations of strong, non-relativistic, parallel electron-ion collisionless shocks using the SHARP code. We compare simulations with a reduced mass ratio ($m_r = 100$) to those with a realistic mass ratio ($m_r = 1836$) for shocks with high ($\mathcal{M}_A = 21.3$) and low ($\mathcal{M}_A = 5.3$) Alfv$\acute{\text{e}}$n Mach numbers. Our findings show that the mass ratio significantly affects particle acceleration and thermal energy dissipation. At high $\mathcal{M}_A$, a reduced mass ratio leads to more efficient electron acceleration and an unrealistically high ion flux at higher momentum. At low $\mathcal{M}_A$, it causes complete suppression of electron acceleration, whereas the realistic mass ratio enables efficient electron acceleration. The reduced mass ratio also results in excessive electron heating and lower heating in downstream ions at both Mach numbers, with slightly more magnetic field amplification at low $\mathcal{M}_A$. Consequently, the electron-to-ion temperature ratio is high at low $\mathcal{M}_A$ due to reduced ion heating and remains high at high $\mathcal{M}_A$ due to increased electron heating. In contrast, simulations with the realistic $m_r$ show that the ion-to-electron temperature ratio is independent of the upstream magnetic field, a result not observed in reduced $m_r$ simulations.

Young and forming planetesimals experience impacts from particles present in a protostellar disk. Using crater scaling laws, we integrate ejecta distributions for oblique impacts. For impacts at 10 to 65 m/s, expected for impacts associated with a disk wind, we estimate the erosion rate and torque exerted on the planetesimal. We find that the mechanism for angular momentum drain proposed by Dobrovolskis and Burns (1984) for asteroids could operate in the low velocity regime of a disk wind. Though spin-down associated with impacts can facilitate planetesimal collapse, we find that the process is inefficient. We find that angular momentum drain via impacts operates in the gravitational focusing regime, though even less efficiently than for lower mass planetesimals. The angular momentum transfer is most effective when the wind speed is low, the projectile density is high compared to the bulk planetesimal density, and the planetesimal is composed of low-strength material. Due to its inefficiency, we find that angular momentum drain due to impacts within a pebble cloud does not by itself facilitate collapse of single planetesimals.

Recent observations by the James Webb Space Telescope (JWST) have revealed the presence of bright and well-formed galaxies at high redshifts, challenging the predictions of the standard Lambda-Cold Dark Matter (LCDM) cosmological model. This paper explores the potential of Modified Gravity (MOG), specifically Scalar-Tensor-Vector Gravity (STVG), to account for the rapid formation of these galaxies in the early universe. By enhancing the gravitational constant through a dimensionless parameter $\alpha$ and incorporating a massive vector field $\phi_\mu$, MOG predicts deeper gravitational wells that can accelerate the collapse of baryonic matter. We present theoretical insights demonstrating how MOG can facilitate the increase in star formation rate and early formation of galaxies, offering a compelling alternative to LCDM. Our findings suggest that MOG provides a viable framework for understanding the rapid growth of galaxies observed by JWST.

For a galaxy, given its observed rotation curve, can one directly infer parameters of the dark matter density profile (such as dark matter particle mass $m$, scaling parameter $s$, core-to-envelope transition radius $r_t$ and NFW scale radius $r_s$), along with Baryonic parameters (such as the stellar mass-to-light ratio $\Upsilon_*$)? In this work, using simulated rotation curves, we train neural networks, which can then be fed observed rotation curves of dark matter dominated dwarf galaxies from the SPARC catalog, to infer parameter values and their uncertainties. Since observed rotation curves have errors, we also explore the very important effect of noise in the training data on the inference. We employ two different methods to quantify uncertainties in the estimated parameters, and compare the results with those obtained using Bayesian methods. We find that the trained neural networks can extract parameters that describe observations well for the galaxies we studied.

In this paper we present the results of a theoretical study of the trajectories of massive particles in the Köttler metric in view of the cosmological constant {\Lambda}. For both negative and positive signs of {\Lambda} a classification of trajectories is proposed, with entries based on different solutions of the trajectory equation, obtained by the expansion of the corresponding algebraic curve in Puiseux series. We also provide some specific types of trajectories which correspond to different values of the cosmological constant. In the case of negative values of the cosmological constant its upper limit is estimated from the galaxy rotation curves

We study Coupling Extended Proca-Nuevo gravity, a non-linear theory extending from dRGT massive gravity with a spin-1 field. This theory is shown to yield reliable, ghost-free cosmological solutions, modeling both the Universe's thermal history and late-time acceleration. By analyzing data from Dark energy spectroscopic instruments (DESI), Cosmic Chronometer (CCh), Gamma Ray Bursts (GRBs), and Type Ia Supernova (SNeIa), we derive parameter constraints with up to 3$\sigma$ confidence, demonstrating good agreement with observations. Our comparison of $BAO$ data from $WiggleZ$ and $DESI$ highlights its constraining power on the Hubble constant. The analysis of the cosmographic parameter, $q$ shows the statistical compatibility with the recent data. Further, this indicates that Universe's current accelerated expansion aligns with quintessential behavior.

The predictions of inflation are usually defined in terms of equal time in-in correlation functions in an accelerating cosmological background. These same observables exist for quantum field theory in other spacetimes, including flat space. In this paper, we will explore how the Wilsonian renormalization group (RG) and effective field theory (EFT) apply to these observables in both flat and de Sitter space. Specifically, we show that matching the short- and long-distance calculations requires additional terms localized at the time of the measurement that are not captured by the effective action of the EFT. These additional terms only correct the local and semi-local terms in the EFT correlators. In flat space, we give an explicit demonstration by matching in-in correlators of light scalars interacting with a heavy field with the EFT result. We then show how these additional terms arise generically via exact RG. We also compare these explicit results in flat space with the corresponding theory in de Sitter and show that the local terms typically redshift away. Our results are closely related to momentum space entanglement that arises from tracing over short-wavelength modes.

Dense neutrino gases can exhibit collective flavor instabilities, triggering large flavor conversions that are driven primarily by neutrino-neutrino refraction. One broadly distinguishes between fast instabilities that exist in the limit of vanishing neutrino masses, and slow ones, that require neutrino mass splittings. In a related series of papers, we have shown that fast instabilities result from the resonant growth of flavor waves, in the same way as turbulent electric fields in an unstable plasma. Here we extend this framework to slow instabilities, focusing on the simplest case of an infinitely homogeneous medium with axisymmetric neutrino distribution. The relevant length and time scales are defined by three parameters: the vacuum oscillation frequency $\omega_E=\delta m^2/2E$, the scale of neutrino-neutrino refraction energy $\mu=\sqrt{2}G_F(n_\nu+n_{\overline\nu})$, and the ratio between lepton and particle number $\epsilon=(n_\nu-n_{\overline\nu})/(n_\nu+n_{\overline\nu})$. We distinguish between two very different regimes: (i) For $\omega_E\ll \mu \epsilon^2$, instabilities occur at small spatial scales of order $(\mu\epsilon)^{-1}$ with a time scale of order $\epsilon \omega_E^{-1}$. This novel branch of slow instability arises from resonant interactions with neutrinos moving along the axis of symmetry. (ii) For $\mu \epsilon^2\ll \omega_E\ll \mu$, the instability is strongly non-resonant, with typical time and length scales of order $1/\sqrt{\omega_E \mu}$. Unstable modes interact with all neutrino directions at once, recovering the characteristic scaling of the traditional studies of slow instabilities. In the inner regions of supernovae and neutron-star mergers, the first regime may be more likely to appear, meaning that slow instabilities in this region may have an entirely different character than usually envisaged.

Some neutron stars known as magnetars possess very strong magnetic fields, with surface fields as large as $10^{15}\,\rm G$ and internal fields that are possibly stronger. Recent observations of the radio pulsar GLEAM-X J1627 suggest it may have a surface field as strong as $10^{16} \,\rm G$. In the presence of a strong magnetic field, the energy levels of electrons and protons are quantized and the Direct Urca process allows neutron stars to cool rapidly, even at low density. For the case of magnetic fields $B \geq 10^{16}\,\rm G$, we find features in the emissivity due to energy quantization that are not captured by the frequently employed quasiclassical approximation where energy levels are treated as nearly continuous. Resonances can result in amplification of the neutrino emissivity at specific densities compared to a calculation that neglects quantization, particularly at low temperature. These effects are not important for the thermal evolution of an entire neutron star, but may be relevant for phenomena that depend on behavior at specific densities. We present a fully relativistic calculation of the Direct Urca rate in a strong magnetic field using the standard V-A weak Lagrangian incorporating mean field nuclear effects and discuss approaches to the numerical challenge the modified wavefunctions present and a new semi-analytic approximation. These tools are also applicable to calculating neutrino opacities in strong magnetic fields in the ejecta of binary neutron star mergers. We calculate the opacities for neutrinos capturing on free nucleons at sub-saturation densities and temperatures exceeding an MeV. We find an enhancement to capture processes of the lowest energy neutrinos by an order of magnitude or more due to suppression of electron Pauli blocking in the case of capture on neutrons, and from the effect of the nucleon magnetic moments in the case of capture on protons.

We investigate the global structures of neutron stars within the framework of general relativity, treating the entire star as a quantum-degenerate system. Rather than relying on the Tolman-Oppenheimer-Volkoff (TOV) equation, we solve the Einstein-Cartan (EC) field equations self-consistently, incorporating the energy-momentum tensor contributions from neutrons. Neutron wave functions are obtained by solving the Dirac equation in a curved spacetime with both torsion and curvature effects. Given that neutron stars contain about 10^57 particles, we adopt a scaled h-bar approach to efficiently describe the quantum state of highly degenerate system.

Cosmography has been extensively utilized to constrain the kinematic state of the Universe using measured distances. In this work, we propose a new method to reconstruct coupling theories using the first kind of Chebyshev polynomial for two variables in which the functional form of the $f(Q,T)$ theory has been obtained. Further, the unknowns that appeared in the series are constrained using the cosmographic parameters. We find the explicit form of the luminosity distance in terms of cosmographic parameters to perform MCMC analysis using the PANTHEON+SH0ES data set. Through the distance modulus function, we observe that the result comes out to be an excellent match to the standard cosmological model and data.

Decomposing theoretical nuclear mass predictions into a liquid-drop parametrization and local shell effects shows that r-process abundances are virtually insensitive to large variations of the masses which originate from nuclear bulk properties of the model, such as the symmetry energy. Therefore, experimental and theoretical studies of masses devoted to r-process applications, such as the nucleosynthesis in the ejecta of neutron star mergers, should focus on the physical origin of local changes in mass trends without necessarily providing highly accurate mass determinations of individual nuclei.

For the central values of the relevant experimental inputs, that is the strong coupling constant and the top quark and Higgs masses, the effective Higgs potential displays two minima, one at the electroweak scale and a deeper one at high energies. We review the phenomenology of the Higgs inflation model, extending the Standard Model to include a non-minimal coupling to gravity; as recently shown, even configurations that would be metastable in the Standard Model, become viable for inflation if the non-minimal coupling is large enough to flatten the Higgs potential at field values below the barrier between the minima.

We investigate the impact of hyperons and phase transition to quark matter on the structural properties of neutron stars within the four-dimensional Einstein-Gauss-Bonnet gravity (EGB). We employ the density-dependent relativistic mean-field model (DDME2) for the hadronic phase and the density-dependent quark mass (DDQM) model for the quark phase to construct hadronic and hybrid equations of state (EoSs) that are consistent with the astrophysical constraints. The presence of hyperons softens the EoS and with a phase transition, the EoS further softens, and the speed of sound squared drops to around 0.2 for the maximum mass configuration which lies in the pure quark phase. Adjusting the Gaussian-Bonnet coupling constant $\alpha$ within its allowed range results in a decrease in the mass-radius relationship for negative $\alpha$, and an increase for positive $\alpha$. In addition, functions are fitted to the maximum mass and its associated radius as a function of constant $\alpha$ to observe its impact on these properties.

Numerous observations on astrophysical and cosmological scales can be interpreted to mean that, in addition to the familiar kind of matter well described by the standard model of elementary particle physics, there exists Dark Matter (DM). The fundamental properties of the elementary particles which make up the DM e.g. particle mass, spin, couplings etc are currently being observationally constrained. In particular, if DM particles have spin zero, there exist recent constraints which suggest a lower limit on its mass which is often a couple of orders of magnitude larger than $10^{-22}$ eV. In this talk, we will (a) argue that these limits are based on the assumption that the self coupling of the spinless DM particles is negligible, and, (b) show how some of these lower limits will get modified in the presence of incredibly feeble self interactions.

Gravitational waves (GWs) from binary neutron stars (BNSs) offer valuable understanding of the nature of compact objects and hadronic matter. However, their analysis requires substantial computational resources due to the challenges in Bayesian stochastic sampling. The third-generation (3G) GW detectors are expected to detect BNS signals with significantly increased signal duration, detection rates, and signal strength, leading to a major computational burden in the 3G era. We demonstrate a machine learning-based workflow capable of producing source parameter estimation and constraints on equations of state (EOSs) for hours-long BNS signals in seconds with minimal hardware costs. We employ efficient compressions on the GW data and EOS using neural networks, based on which we build normalizing flows for inferences. Given that full Bayesian analysis is prohibitively time-intensive, we validate our model against (semi-)analytical predictions. Additionally, we estimate the computational demands of BNS signal analysis in the 3G era, showing that the machine learning methods will be crucial for future catalog-level analysis.

Recent results from DESI combined with cosmic microwave background data give the tightest constraints on the sum of neutrino masses to date. However, these analyses approximate the neutrino mass hierarchy by three degenerate-mass (DM) neutrinos, instead of the normal (NH) and inverted hierarchies (IH) informed by terrestrial neutrino oscillation experiments. Given the stringency of the upper limits from DESI data, we test explicitly whether the inferred neutrino constraints are robust to the choice of neutrino mass ordering using both Bayesian and frequentist methods. For Planck data alone, we find that the DM hierarchy presents a good approximation to the physically motivated hierarchies while showing a strong dependence on the assumed lower bound of the prior, confirming previous studies. For the combined Planck and DESI baryon acoustic oscillation data, we find that assuming NH ($M_\mathrm{tot} < 0.13\,\mathrm{eV}$) or IH ($M_\mathrm{tot} < 0.16\,\mathrm{eV}$) loosens the Bayesian upper limits compared to the DM approximation ($M_\mathrm{tot} < 0.086\,\mathrm{eV}$). The frequentist analysis shows that the different neutrino models fit the data equally well and the loosening of the constraints can thus be attributed to the lower bounds induced by NH and IH. Overall, we find that the DM hierarchy presents a good approximation to the physically motivated hierarchies also for Planck+DESI data as long as the corresponding lower neutrino mass bounds are imposed.