Papers for Friday, Jun 21 2024

Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, from pulsar timing, and from eclipsing binary stars in the Solar Neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density, and from this, by accounting for the mass in stars, gas, and dust, we can infer the amount of dark matter. Leveraging the FIRE-2 simulations of Milky Way-mass galaxies, we compare vertical acceleration profiles between cold dark matter (CDM) and self-interacting dark matter (SIDM) with constant cross-section of 1 cm$^2$ g$^{-1}$ across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos consistently exhibit higher local stellar and dark matter densities and steeper vertical acceleration gradients, up to 30% steeper near the Solar Neighborhood. SIDM halos also manifest a more oblate halo shape in the Solar Neighborhood. Furthermore, enhanced precision in acceleration measurements and larger datasets promise to provide better constraints on the local dark matter density, complementing our understanding from kinematic analysis of their distribution within galaxies.

The elusive properties of the first (Pop III) stars can be indirectly unveiled by uncovering their true descendants. To this aim, we exploit our data-calibrated model for the best-studied ultra-faint dwarf (UFD) galaxy, Boötes I, which tracks the chemical evolution (from carbon to zinc) of individual stars from their formation to the present day. We explore the chemical imprint of Pop III supernovae (SNe), with different explosion energies and masses, showing that they leave distinct chemical signatures in their descendants. We find that UFDs are strongly affected by SNe-driven feedback resulting in a very low fraction of metals retained by their gravitational potential well (< 2.5 %). Furthermore, the higher the Pop III SN explosion energy, the lower the fraction of metals retained. Thus, the probability to find descendants of energetic Pair Instability SNe is extremely low in these systems. Conversely, UFDs are ideal cosmic laboratories to identify the fingerprints of less massive and energetic Pop III SNe through their [X/Fe] abundance ratios. Digging into the literature data of Boötes I, we uncover three hidden Pop III descendants: one mono-enriched and two multi-enriched. These stars show the chemical signature of Pop III SNe in the mass range $[20-60]\rm M_{\odot}$, spanning a wide range in explosion energies $[0.3-5] 10^{51}$ erg. In conclusion, Pop III descendants are hidden in ancient UFDs but those mono-enriched by a single Pop III SN are extremely rare. Thus, self-consistent models such as the one presented here are required to uncover these precious fossils and probe the properties of the first Pop III supernovae.

The kinematics of the stellar halo hold important clues to the assembly history and mass distribution of the Galaxy. In this study, we map the kinematics of stars across the Galactic halo with the H3 Survey. We find a complex distribution that breaks both azimuthal symmetry about the $Z$-axis and mirror symmetry about the Galactic plane. This asymmetry manifests as large variations in the radial velocity dispersion $\sigma_r$ from as ``cold'' as 70 $\text{km}\text{ s}^{-1}$ to as ``hot'' as 160 $\text{km}\text{ s}^{-1}$. We use stellar chemistry to distinguish accreted stars from in-situ stars in the halo, and find that the accreted population has higher $\sigma_r$ and radially biased orbits, while the in-situ population has lower $\sigma_r$ and isotropic orbits. As a result, the Galactic halo kinematics are highly heterogeneous and poorly approximated as being spherical or axisymmetric. We measure radial profiles of $\sigma_r$ and the anisotropy parameter $\beta$ over Galactocentric radii $10-80\text{ kpc}$, and find that discrepancies in the literature are due to the nonspherical geometry and heterogeneous nature of the halo. Investigating the effect of strongly asymmetric $\sigma_r$ and $\beta$ on equilibrium models is a path forward to accurately constraining the Galactic gravitational field, including its total mass.

One well known method of generating a large blue spectral index for axionic isocurvature perturbations is through a flat direction not having a quartic potential term for the radial partner of the axion field. In this work, we show how one can obtain a large blue spectral index even with a quartic potential term associated with the Peccei-Quinn symmetry breaking radial partner. We use the fact that a large radial direction with a quartic term can naturally induce a conformal limit which generates an isocurvature spectral index of 3. We point out that this conformal representation is intrinsically different from both the ordinary equilibrium axion scenario or massless fields in Minkowski spacetime. Another way to view this limit is as a scenario where the angular momentum of the initial conditions slows down the radial field or as a superfluid limit. Quantization of the non-static system in which derivative of the radial field and the derivative of the angular field do not commute is treated with great care to compute the vacuum state. The parametric region consistent with axion dark matter and isocurvature cosmology is discussed.

We investigate low-mass central galaxies with Mstar = $10^{9.5}-10^{10}$ Msun/h, located near massive groups and galaxy clusters using the TNG300 and MDPL2-SAG simulations. We set out to study their evolution, aiming to find hints about the large-scale conformity signal they produce. We also use a control sample of low-mass central galaxies located far away from massive structures. For both samples, we find a sub-population of galaxies that were accreted by another halo in the past but are now considered central galaxies; we refer to these objects as former satellites. The fraction of former satellites is higher for quenched central galaxies near massive systems: 45% in TNG300 and 17% in MDPL2-SAG. Our results in TNG300 show that former satellites were typically hosted by massive dark matter halos (M200 $\geq 10^{13}$ Msun/h) at z$\sim$0.3, followed by a drop in halo mass at lower redshifts. In addition, we find a strong drop in the total gas mass at z$\leq$1 for quenched central galaxies near galaxy groups and clusters produced by these former satellites as well. By removing former satellites, the evolution of quenched central galaxies is fairly similar to those of the quenched control galaxies, showing small differences at low-z. For MDPL2-SAG, former satellites were hosted by less massive halos, with a mean halo mass around $10^{11}$ Msun/h, and the evolution remains equal before and after removing former satellites. We also measure the two-halo conformity, i.e., the correlation in the specific SFR between low-mass central galaxies and their neighbors at Mpc scales, and how former satellites contribute to the signal at z=0, 0.3, and 1. The conformity signal decreases from z=0 to z=1 in MDPL2-SAG but it increases in TNG300. However, after removing former satellites in TNG300, the signal is strongly reduced but almost does not change at z$\leq$0.3, and it disappears at z=1 (abridged).

Line-intensity mapping (LIM) offers an approach to obtain three-dimensional maps of the large-scale structure by collecting the aggregate emission from all emitters along the line of sight. The procedure hinges on reconstructing the radial positions of sources by relating the observed frequency to the rest-frame frequency of a target emission line. However, this step is hindered by `interloper-line' emission from different cosmological volumes that redshifts into the same observed frequency. In this work, we propose a model-independent technique to remove the contamination of line interlopers using their statistical correlation with external tracers of the large-scale structure, and identify the weights that minimize the variance of the cleaned field. Furthermore, we derive expressions for the resulting power spectra after applying our cleaning procedure, and validate them against simulations. We find that the cleaning performance improves as the correlation between the line interlopers and the external tracer increases, resulting in a gain in the signal-to-noise ratio of up to a factor 6 (2) for the auto- (cross-)power spectrum in idealized scenarios. This approach has the advantage of being model-independent, and is highly complementary to other techniques, as it removes large-scale clustering modes instead of individually masking the brightest sources of contamination.

We present ALMA Band 7 molecular line observations of the protostars within the VLA 1623 system. We map C$^{17}$O (3 - 2) in the circumbinary disk around VLA 1623A and the outflow cavity walls of the collimated outflow. We further detect red-shifted and blue-shifted velocity gradients in the circumstellar disks around VLA 1623B and VLA 1623W that are consistent with Keplerian rotation. We use the radiative transfer modeling code, pdspy, and simple flared disk models to measure stellar masses of $0.27 \pm 0.03$ M$_\odot$, $1.9^{+0.3}_{-0.2}$ M$_\odot$, and $0.64 \pm 0.06$ M$_\odot$ for the VLA 1623A binary, VLA 1623B, and VLA 1623W, respectively. These results represent the strongest constraints on stellar mass for both VLA 1623B and VLA 1623W, and the first measurement of mass for all stellar components using the same tracer and methodology. We use these masses to discuss the relationship between the young stellar objects (YSOs) in the VLA 1623 system. We find that VLA 1623W is unlikely to be an ejected YSO, as has been previously proposed. While we cannot rule out that VLA 1623W is a unrelated YSO, we propose that it is a true companion star to the VLA 1623A/B system and that the these stars formed in situ through turbulent fragmentation and have had only some dynamical interactions since their inception.

The frequency spectrum of the cosmic microwave background (CMB) is a relatively untapped source of data which can allow us to peer beyond the surface of last scattering. Small deviations away from a perfect blackbody shape will encode valuable information about the state of the primordial Universe which may not be accessible by other means. Here, we briefly review some key science goals of CMB spectral distortions, with an emphasis on how future generations of experiments can be used in tandem with complementary observational probes to perform model discrimination of exotic physics scenarios. We focus here on synergies between spectral distortions, gravitational waves, and 21cm cosmology.

Protoclusters at $z>2$ are gas-rich regions characterized by high star-formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of SMBHs. We aim to test this scenario by probing the AGN content of SPT2349-56, a massive, gas-rich, and highly star-forming protocluster core at $z=4.3$ discovered as an overdensity of DSFGs, via Chandra (200 ks) observations, and comparing the results with the field environment. We detected two protocluster members, corresponding to an AGN fraction among DSFGs of $\approx10\%$. This value is consistent with other protoclusters at $z=2-4$, but higher than the AGN incidence among DSFGs in the field environment. Both AGN are heavily obscured sources and hosted in star-forming galaxies with $\approx3\times10^{10}\,\mathrm{M_\odot}$ stellar masses. We estimate that the ISM in the host galaxies can contribute significantly to the nuclear obscuration. One of the two AGN is highly luminous ($L_X=2\times10^{45}\,\mathrm{erg\,s^{-1}}$) and Compton-thick ($N_H=2\times10^{24}\,\mathrm{cm^{-2}}$), and likely powered by a $M_{BH}>6\times10^8\,\mathrm{M_\odot}$ SMBH. Its high accretion rate suggests that it is in the phase of efficient growth required to explain the presence of extremely massive SMBHs in the centers of local galaxy clusters. Considering SPT2349-56 and DRC, a similar protocuster at $z=4$, we find that gas-rich protocluster cores at $z\approx4$ enhance the triggering of luminous (log$\frac{L_X}{\mathrm{erg\,s^{-1}}}=45-46$) AGN by 3-5 orders of magnitude with respect to the field environment. Our results indicate that gas-rich protoclusters at high redshift boost the growth of SMBHs, which will likely impact the subsequent evolution of the structures, and thus represent key science targets to obtain a complete understanding of the relation between environment and galaxy evolution.

As part of the Young Exoplanets Spectroscopic Survey (YESS), this study explores the spot variability of 13 T Tauri Stars (TTSs) in the near-infrared $H$ band, using spectra from the Immersion GRating INfrared Spectrometer (IGRINS). By analyzing effective temperature ($T_{\rm eff}$) sensitive lines of atomic FeI at ~1.56259 um and ~1.56362 um, and molecular OH at ~1.56310 um and ~1.56317 um, we develop an empirical equivalent width ratio (EWR) relationship for $T_{\rm eff}$ in the range of 3400-5000 K. This relationship allows for precise relative $T_{\rm eff}$ estimates to within tens of Kelvin and demonstrates compatibility with solar metallicity target models. However, discrepancies between observational data and model predictions limit the extension of the $T_{\rm eff}$-EWR relationship to a broader parameter space. Our study reveals that both classical and weak-line TTSs can exhibit $T_{\rm eff}$ variations exceeding 150 K over a span of two years. The detection of a quarter-phase delay between the EWR and radial velocity phase curves in TTSs indicates spot-driven signals. A phase delay of 0.06 $\pm$ 0.13 for CI Tau, however, suggests additional dynamics, potentially caused by planetary interaction, inferred from a posited 1:1 commensurability between the rotation period and orbital period. Moreover, a positive correlation between $T_{\rm eff}$ variation amplitude and stellar inclination angle support the existence of high-latitude spots on TTSs, further enriching our understanding of stellar surface activity in young stars.

The Mrk 463 system is known to host two powerful sources separated by about 4 kpc, both identified as active galactic nuclei (AGN). This makes the Mrk 463 system a unique laboratory to study the geometry and dynamics of galaxy merging and its relation to AGN duty cycles. The eastern nuclei, Mrk 463E, is the brightest of the two and thus a prime target for a polarimetric study. It is classified as a Seyfert 2 galaxy, meaning that one could expect large polarization degrees from scattering off electrons and dust in the polar winds. In the continuity of our series of papers, we reduced archived and previously unpublished polarization observations obtained with the Faint Object Camera (FOC) onboard the Hubble Space Telescope (HST), to obtain a high resolution near ultraviolet (near-UV) polarization map of the Mrk 463E nuclei. We coupled this map to near infrared (NIR) and X-ray observations to get a clear picture of the geometric arrangement of matter around the core of Mrk 463E. We found that the nucleus location is further South from the optical peak flux than previously estimated. The strongly polarized conical wind has a half-opening angle of ~15° and display three main periods of mass ejection. Its polarization allowed us to estimate the AGN inclination towards the observer (~55°) Finally, our maps revealed a gas streamer connecting Mrk 463E and Mrk 463W, with a tentative detection of a large kpc-scale ordered magnetic field connecting both galaxies. This unpublished observation turned out to offer more than the original proposal asked for and allowed to derive tight geometric and dynamical constraints for Mrk 463E. High resolution radio maps and IR polarimetry are now necessary to further study the jet and the newly discovered gas streamer.

Magnetic flux ropes within interplanetary coronal mass ejections are often characterized as simplistic cylindrical or toroidal tubes with field lines that twist around the cylinder or torus axis. Recent multi-point observations suggest that the overall geometry of these large-scale structures may be significantly more complex, so that the contemporary modeling approaches would be, in some cases, insufficient to properly understand the global structure of any interplanetary coronal mass ejection. In an attempt to partially rectify this issue, we have developed a novel magnetic flux rope model that allows for the description of arbitrary distortions of the cross-section or deformation of the magnetic axis. The distorted magnetic flux rope model is a fully analytic flux rope model, that can be used to describe significantly more complex geometries and is numerically efficient enough to be used for large ensemble simulations. To demonstrate the usefulness of our model, we focus on a specific implementation of our model and apply it to an ICME event that was observed \textit{in situ} on 2023 April 23 at the L1 point by the Wind spacecraft and also by the STEREO-A spacecraft that was $10.2^\circ$ further east and $0.9^\circ$ south in heliographic coordinates. We demonstrate that our model can accurately reconstruct each observation individual and also gives a fair reconstruction of both events simultaneously using a multi-point reconstruction approach, which results in a geometry that is not fully constistent with a cylindrical or toroidal approximation.

Depending principally on mass, the compact object remnant left behind after a star's life may be a white dwarf (WD), neutron star (NS), or black hole (BH). While we have large samples of each of these remnants, we lack knowledge of the exact conditions separating these outcomes. The boundary between low-mass BHs and massive NSs is particularly poorly understood, as few objects between 2-5 $M_\odot$ are known. To probe this regime, we search the APOGEE DR17 dataset of 657,000 unique stars for binary systems with one stellar component that exhibit large radial velocity shifts over multiple observations. We identify 4751 likely binary systems, and estimate a minimum mass for each system's "invisible companion" under the assumption of tidal synchronization. Two systems have companion masses $\gtrsim$ 2 $M_\odot$, although we conclude that neither are good candidates for possessing a mass-gap NS or BH companions.

In 2021, the high-mass X-ray binary EXO 2030+375 underwent a giant X-ray outburst, the first since 2006, that reached a peak flux of ${\sim}600\,\mathrm{mCrab}$ (3-50\,keV). The goal of this work is to study the spectral evolution over the course of the outburst, search for possible cyclotron resonance scattering features (CRSFs), and to associate spectral components with the emission pattern of the accretion column. We used broadband spectra taken with the Nuclear Spectroscopic Telescope Array (NuSTAR), the Neutron Star Interior Composition Explorer (NICER), and Chandra near the peak and during the decline phase of the outburst. We describe the data with established empirical continuum models and perform pulse-phase-resolved spectroscopy. We compare the spectral evolution with pulse phase using a proposed geometrical emission model. We find a significant spectral hardening toward lower luminosity, a behavior that is expected for super-critical sources. The continuum shape and evolution cannot be described by a simple power-law model with exponential cutoff; it requires additional absorption or emission components. We can confirm the presence of a narrow absorption feature at ${\sim}10\,\mathrm{keV}$ in both NuSTAR observations. The absence of harmonics puts into question the interpretation of this feature as a CRSF. The empirical spectral components cannot be directly associated with identified emission components from the accretion column.

Following the recommendations to NASA and ESA, the search for life on exoplanets will be a priority in the next decades. Two direct imaging space mission concepts are being developed: the Habitable Worlds Observatory (HWO) and the Large Interferometer for Exoplanets (LIFE). HWO focuses on reflected light spectra in the ultraviolet/visible/near-infrared (UV/VIS/NIR), while LIFE captures the mid-infrared (MIR) emission of temperate exoplanets. We assess the potential of HWO and LIFE in characterizing a cloud-free Earth twin orbiting a Sun-like star at 10 pc, both separately and synergistically, aiming to quantify the increase in information from joint atmospheric retrievals on a habitable planet. We perform Bayesian retrievals on simulated data from an HWO-like and a LIFE-like mission separately, then jointly, considering the baseline spectral resolutions currently assumed for these concepts and using two increasingly complex noise simulations. HWO would constrain H$_2$O, O$_2$, and O$_3$, in the atmosphere, with ~ 100 K uncertainty on the temperature profile. LIFE would constrain CO$_2$, H$_2$O, O$_3$ and provide constraints on the thermal atmospheric structure and surface temperature (~ 10 K uncertainty). Both missions would provide an upper limit on CH$_4$. Joint retrievals on HWO and LIFE data would accurately define the atmospheric thermal profile and planetary parameters, decisively constrain CO$_2$, H$_2$O, O$_2$, and O$_3$, and weakly constrain CO and CH$_4$. The detection significance is greater or equal to single-instrument retrievals. Both missions provide specific information to characterize a terrestrial habitable exoplanet, but the scientific yield is maximized with synergistic UV/VIS/NIR+MIR observations. Using HWO and LIFE together will provide stronger constraints on biosignatures and life indicators, potentially transforming the search for life in the universe.

We present two frameworks to infer some of the properties of neutron stars from their electromagnetic radiation and the emission of continuous gravitational waves due to r-mode oscillations. In the first framework, assuming a distance measurement via electromagnetic observations, we infer three neutron star properties: the moment of inertia, a parameter related to the r-mode saturation amplitude, and the component of magnetic dipole moment perpendicular to the rotation axis. Unlike signals from mountains, r-mode oscillations provide additional information through a parameter (\kappa) that satisfies a universal relation with the star's compactness. In the second framework, we utilize this and the relation between the moment of inertia and compactness, in addition to assuming an equation of state and utilizing pulsar frequency measurements, to directly measure the neutron star's distance along with the aforementioned parameters. We employ a Fisher information matrix-based approach for quantitative error estimation in both frameworks. We find that the error in the distance measurement dominates the errors in the first framework for any reasonable observation time. In contrast, due to the low errors in pulsar frequency measurements, parameters can be inferred accurately via the second framework but work only in a restricted parameter space. We finally address potential ways to overcome critical drawbacks of our analyses and discuss directions for future work.

The galactic center excess is a possible non-gravitational observation of dark matter; however, the canonical dark matter model (thermal freeze-out) is in conflict with other gamma-ray observations, in particular those made of the Milky Way's satellite dwarf galaxies. Here we consider the effects of a two-component dark matter model which results in minimally boosted particles that must remain bound to their host galaxy in order to produce an observational signal. This leads to a signal that is heavily dependent on galactic scale and can help reconcile the differences in the galactic center and dwarf galaxy measurements under the dark matter paradigm.

More than 300 supermassive black holes have been detected at redshifts larger than 6, and they are abundant in the centers of local galaxies. Their formation mechanisms, however, are still rather unconstrained. A possible origin of these supermassive black holes could be through mergers in dense black hole clusters, forming as a result of mass segregation within nuclear star clusters in the center of galaxies. In this study, we present the first systematic investigation of the evolution of such black hole clusters where the effect of an external potential is taken into account. Such a potential could be the result of gas inflows into the central region, for example as a result of galaxy mergers. We show here that the efficiency for the formation of a massive central object is mostly regulated by the ratio of cluster velocity dispersion divided by the speed of light, potentially reaching efficiencies of 0.05-0.08 in realistic systems. Our results show that this scenario is potentially feasible and may provide seeds black hole of at least 10^3 solar masses. We conclude that the formation of seed black holes via this channel should be taken into account in statistical assessments of the black hole population.

A detection of the 21 cm signal can provide a unique window of opportunity for uncovering complex astrophysical phenomena at the epoch of reionization and placing constraints on cosmology at high redshifts, which are usually elusive to large-scale structure surveys. In this work, we provide a theoretical model based on a quadratic bias expansion capable of recovering the 21 cm power spectrum with high accuracy sufficient for upcoming ground-based radio interferometer experiments. In particular, we develop a hybrid effective field theory (HEFT) model in redshift space that leverages the accuracy of $N$-body simulations with the predictive power of analytical bias expansion models, and test it against the Thesan suite of radiative transfer hydrodynamical simulations. We make predictions of the 21 cm brightness temperature field at several distinct redshifts, ranging between $z = 6.5$ and 11, thus probing a large fraction of the reionization history of the Universe ($x_{\rm HI} = 0.3 \sim 0.9$), and compare our model to the `true' 21 cm brightness in terms of the correlation coefficient, power spectrum and modeling error. We find percent-level agreement at large and intermediate scales, $k \lesssim 0.5 h/{\rm Mpc}$, and favorable behavior down to small scales, $k \sim 1 h/{\rm Mpc}$, outperforming pure perturbation-theory-based models. To put our findings into context, we show that even in the absence of any foreground contamination the thermal noise of a futuristic HERA-like experiment is comparable with the theoretical uncertainty in our model in the allowed `wedge' of observations, providing further evidence in support of using HEFT-based models to approximate a range of cosmological observables.

We present deep 1.4-4.8 um JWST-NIRCam imaging of the Serpens Main star-forming region and identify 20 candidate protostellar outflows, most with bipolar structure and identified driving sources. The outflow position angles (PAs) are strongly correlated, and aligned within +/- 24 degrees of the major axis of the Serpens filament. These orientations are further aligned with the angular momentum vectors of the two disk shadows in this region. We estimate that the probability of this number of young stars being co-aligned if sampled from a uniform PA distribution is 10^-4. This in turn suggests that the aligned protostars, which seem to be at similar evolutionary stages based on their outflow dynamics, formed at similar times with a similar spin inherited from a local cloud filament. Further, there is tentative evidence for a systematic change in average position angle between the north-western and south-eastern cluster, as well as increased scatter in the PAs of the south-eastern protostars. SOFIA-HAWC+ archival dust polarization observations of Serpens Main at 154 and 214 um are perpendicular to the dominant jet orientation in NW region in particular. We measure and locate shock knots and edges for all of the outflows and provide an identifying catalog. We suggest that Serpens main is a cluster that formed from an isolated filament, and due to its youth retains its primordial outflow alignment.

Coronal mass ejections (CMEs) and their driven shocks are a major source of large geomagnetic storms due to their large and long-lasting, southward component of magnetic field in the sheath and the flux rope (e.g., magnetic cloud). Predicting the strength and arrival time of southward fields accurately thus plays a key role in space weather predictions. To address this problem, we have developed a new model, which combines the global three-dimensional, time-dependent, magnetohydrodynamic (MHD), data-driven model (G3DMHD) and a self-contained magnetic flux-rope model [1]. As a demonstration and validation, here we simulate the evolution of a Sun-Earth-directed CME that erupted on 2012-July-12. The computational domain spans from 2.5 solar radii (Rs) from the surface of the Sun, where the flux rope is injected, to 245 Rs. We compare the time profiles of the simulated MHD parameters (Density, velocity, temperature, and magnetic field) with in situ solar wind observations acquired at ~1 AU by the Wind spacecraft and the result is encouraging. The model successfully reproduces the shock, sheath, and flux rope similar to those observed by Wind.

Sunspots and active regions observed on the solar surface are widely believed to be manifestations of compact predominantly-toroidal magnetic field structures (``flux tubes") that emerge by magnetic buoyancy from the deeper interior of the Sun. Much work has examined the evolution of such magnetic structures, typically considering them as idealized isolated magnetic entities and not as more realistic magnetic concentrations in a volume-filling background magnetic field. Here, we report results that explore the buoyant rise dynamics of magnetic concentrations in a volume-filling field in the full three dimensions. Earlier 2.5D work in this series (arXiv:1805.08806, arXiv:2101.03472, arXiv:2204.13078) established the remarkable fact that the twist orientation of a flux concentration relative to the background field affected it's likelihood to rise and emerge, regardless of whether the buoyant rise took place in the absence or presence of convection. The contrasting dynamics between structures with differing orientations leads to a selection mechanism that reproduces characteristics of the ``solar hemispheric helicity rule(s)" (SHHR) observations strikingly well. Here, we show that this two-dimensional selection mechanism persists in the face of the added complexity of three-dimensional dynamics. Arching of the magnetic structure in the third dimension, as might be expected in the solar application, is introduced. The role of tension force leading to this selection mechanism is elucidated and subtle differences that arise due to the three-dimensional geometry are discussed.

The unstable mass transfer situation in binary systems will asymptotically cause the adiabatic expansion of the donor star and finally lead to the common envelope phase. This process could happen in helium binary systems once the helium donor star fills its Roche-lobe. We have calculated the adiabatic mass loss model of naked helium stars with a mass range of 0.35\,$M_{\odot}$ to 10\,$M_{\odot}$, and every mass sequence evolved from the He-ZAMS to the cooling track of white dwarf or carbon ignition. In consideration of the influence of stellar wind, massive helium stars are not considered in this paper. Comparing stellar radius with the evolution of the Roche-lobe under the assumption of conservative mass transfer, we give the critical mass ratio $q_{\textrm{crit}}=M_{\textrm{He}}/M_{\textrm{accretor}}$ as the binary stability criteria of low and intermediate-mass helium binary stars. On He-MS, the result shows $1.0<q_{\textrm{crit}}<2.6$, which is more unstable than the classical result of polytropic model $q_{\textrm{crit}}=3$. After early He-HG, the $q_{\textrm{crit}}$ quickly increases even larger than 10 (more stable compared with widely used result $q_{\textrm{crit}}=4$), which is dominated by the expansion of radiative envelope. Our result could be useful for these quick mass transfer binary systems such as AM CVns, UCXBs, and helium novae, and it could guide the binary population synthesis for the formation of special objects such as SNe Ia and GW sources.

The SLAC Microresonator Radio Frequency (SMuRF) electronics is being deployed as the readout for the Cosmic Microwave Background (CMB) telescopes of the Simons Observatory (SO). A Radio Frequency System-on-Chip (RFSoC) based readout of microwave frequency resonator based cryogenic sensors is under development at SLAC as an upgrade path for SMuRF with simplified RF hardware, a more compact footprint, and lower total power consumption. The high-speed integrated data converters and digital data path in RFSoC enable direct RF sampling without analog up and down conversion for RF frequencies up to 6 GHz. A comprehensive optimization and characterization study has been performed for direct RF sampling for microwave SQUID multiplexers, which covers noise level, RF dynamic range, and linearity using a prototype implementation. The SMuRF firmware, including the implementation of closed-loop tone tracking, has been ported to the RFSoC platform and interfaced with the quadrature mixers for digital up and down conversion in the data converter data path to realize a full microwave SQUID multiplexer readout. In this paper, a selection of the performance characterization results of direct RF sampling for microwave SQUID multiplexer readout will be summarized and compared with science-driven requirements. Preliminary results demonstrating the read out of cryogenic sensors using the prototype system will also be presented here. We anticipate our new RFSoC-based SMuRF system will be an enabling readout for on-going and future experiments in astronomy and cosmology, which rely on large arrays of cryogenic sensors to achieve their science goals.

Quasars are generally divided into jetted radio-loud and non-jetted radio-quiet ones, but why only 10% quasars are radio loud has been puzzling for decades. Other than jet-induced-phenomena, black hole mass, or Eddington ratio, prominent difference between jetted and non-jetted quasars has scarcely been detected. Here we show a unique distinction between them and the mystery of jet launching could be disclosed by a prominent excess of radio emission in extremely stable quasars (ESQs, i.e., type 1 quasars with extremely weak variability in UV/optical over 10 years). Specifically, we find that $>$ 25% of the ESQs are detected by the FIRST/VLASS radio survey, while only $\sim$ 6-8% of the control sample, matched in redshift, luminosity, and Eddington ratio, are radio-detected. The excess of radio detection in ESQs has a significance of 4.4 $\sigma$ (99.9995%), and dominantly occurs at intermediate radio loudness with R $\sim$ 10 - 60. The radio detection fraction of ESQs also tends to increase in the ESQ samples selected with more stringent thresholds. Our results are in contrast to the common view that RL quasars are likely more variable in UV/optical due to jet contribution. New clues/challenge posed by our findings highlight the importance of extensive follow-up observations to probe the nature of jets in ESQs, and theoretical studies on the link between jet launching and ESQs. Moreover, our results makes ESQs, an essential population which has never been explored, unique targets in the burgeoning era of time domain astronomy, like their opposite counterparts of quasars exhibiting extreme variability or changing-look features.

The Atacama Large Millimeter/submillimeter Array remains the largest mm radio interferometer observatory world-wide. It is now conducting its 11th observing cycle. In our previous paper presented at this conference series in 2020, we outlined a number of possible improvements to the ALMA end-to-end observing and data processing procedures which could further optimize the uv coverage and thus the image quality while at the same time improving the observing efficiency. Here we report an update of our results refining our proposed adjustments to the scheduling and quality assurance processes. In particular we present new results on ways to assess the uv coverage of a given observation efficiently, methods to define and measure the maximum recoverable angular scale, and on the robustness of the deconvolution in the final interferometric imaging process w.r.t. defects in the uv coverage. Finally we present the outline of a design for integrating uv coverage assessment into the control and processing loop of observation scheduling. The results are applicable to all radio interferometers with more than approx. 10 antennas.

Jupiter's Grat Red Spot (GRS) is the largest and longest-lived vortex of all solar system planets but its lifetime is debated and its formation mechanism remains hidden. G. D. Cassini discovered in 1665 the presence of a dark oval at the GRS latitude, known as the "Permanent Spot" (PS) that was observed until 1713. We show from historical observations of its size evolution and motions that PS is unlikely to correspond to the current GRS, that was first observed in 1831. New numerical simulations rule out that the GRS formed by the merging of vortices or by a superstorm, but most likely formed from a flow disturbance between the two opposed Jovian zonal jets north and south of it. If so, the aearly GRS should have had a low tangential velocity so that its rotation velocity has increased over time as it shrunk.

The star formation history (SFH) is a key issue in the evolution of galaxies. In this work, we developed a model based on a Gaussian and gamma function mixture to fit SFHs with varying numbers of components. Our primary objective was to use this model to reveal the shape of SFHs and the corresponding physical driving factors. Specifically, we applied this model to fit SFHs from the TNG100-1 simulation. Our study led to the following findings: 1) Our model fits with TNG star formation histories well, especially for high-mass and red galaxies; 2) A clear relationship exists between the number and shape of fitted components and the mass and color of galaxies, with notable differences observed between central/isolated and satellite galaxies. 3) Our model allowed us to extract different episodes of star formation within star formation histories with ease and analyze the duration and timing of each star formation episode. Our findings indicated a strong relationship between the timing of each star formation episode and galaxy mass and color.

Flux spectrum, event rate, and experimental sensitivity are investigated for the diffuse supernova neutrino background (DSNB), which is originated from past stellar collapses and also known as supernova relic neutrino background. For this purpose, the contribution of collapses that lead to successful supernova (SN) explosion and black hole (BH) formation simultaneously, which are suggested to be a non-negligible population from the perspective of Galactic chemical evolution, is taken into account. If the BH-forming SNe involve the matter fallback onto the protoneutron star for the long term, their total emitted neutrino energy becomes much larger than that of ordinary SNe and failed SNe (BH formation without explosion). The expected event rate according to the current DSNB model is enhanced by up to a factor of two due to the BH-forming SNe, depending on their fraction and the neutrino mass hierarchy. In any case, the operation time required to detect the DSNB at Hyper-Kamiokande would be reduced by such contribution.

For the precise spectral analysis of hot stars, advanced stellar-atmosphere models that consider deviations from the local thermodynamic equilibrium are mandatory. This requires accurate atomic data to calculate all transition rates and occupation numbers for atomic levels in the considered model atoms, not only for a few prominent lines exhibited in an observation. The critical evaluation of atomic data is a challenge because it requires precise laboratory measurements. Ultraviolet spectroscopy of hot stars with high resolving power provide such "laboratory" spectra. We compare observed, isolated lines of the iron group (here calcium to nickel) with our synthetic line profiles to judge the accuracy of the respective oscillator strengths. This will verify them or yield individual correction values to improve the spectral analysis, i.e., the determination of, e.g., effective temperature and abundances. To minimize the error propagation from uncertainties in effective temperature, surface gravity (g), and abundance determination, we start with a precise reanalysis of three hot subdwarf stars, namely EC 11481-2303, Feige 110, and PG 0909+276. Then, we measure the abundances of the iron-group elements individually. Based on identified, isolated lines of these elements, we compare observation and models to measure their deviation in strength (equivalent width). For EC 11481-2303 and Feige 110, we confirmed the previously determined effective temperatures and log g values within their error limits. For all three stars, we fine-tuned all metal abundances to achieve the best reproduction of the observation. For more than 450 isolated absorption lines of the iron group, we compared modeled and observed line strengths. We selected strong, reliable isolated absorption lines, which we recommend to use as reference lines for abundance determinations in related objects.

The IceCube Collaboration has recently reported compelling evidence of high-energy neutrino emission from NGC~1068, and also mild excesses for NGC 4151 and CGCG420-015, local Seyfert galaxies. This has increased the interest along neutrino emission from hot-corona surrounding the super massive black holes of Seyfert Galaxies. In this paper, we revisit phenomenological constraints on the neutrino emission from hot-coronae of seyfert galaxies, using an assumption of equi-ripartition between cosmic-rays and magnetic energy densities. We show that not only these sources are consistent with such an assumption but also that the data point towards low beta plasma parameters inside Seyfert Galaxies. We exploit this finding to constrain the Seyfert diffuse neutrino flux and we obtain that, in order not to overproduce neutrinos, not all the sources can be in an equi-ripartition state. We conclude (along with previous findings) that seyfert galaxies cannot explain the diffuse neutrino spectrum above $\sim 100\, \rm TeV$, allowing space for other astrophysical sources.

Recent sensitive wide-field radio surveys, such as the LOFAR Two Meter Sky Survey (LoTSS), the LOFAR LBA Sky Survey (LoLSS), and the Very Large Array Sky Survey (VLASS), enable the selection of statistically large samples of peaked-spectrum (PS) sources. PS sources are radio sources that have a peak in their radio continuum spectrum and are observed to be compact. They are often considered to be the precursors to large radio galaxies. We present a sample of 8,032 gigahertz-peaked spectrum (GPS) sources with spectral turnovers near 1400 MHz, and a sample of 506 megahertz-peaked spectrum (MPS) sources with turnovers near 144 MHz. Our GPS sample is over five times larger than any previously known sample of PS sources. These large sample sizes allow us to make a robust comparison between GPS sources and MPS sources, such that we can investigate the differences between these types of sources, and study their lifetimes. The shape of the source counts of both samples match that of the general radio-loud active galactic nuclei (AGN) samples, scaled down by a factor 44 $\pm$ 2 for the MPS sample, and a factor 28 $\pm$ 1 for the GPS sample. Assuming no cosmological evolution, these offsets imply that both MPS and GPS sources have shorter duration than general radio-loud AGN, with MPS sources having an $\approx$1.6 times shorter lifespan than GPS sources. The shorter duration of MPS sources relative to GPS sources can be explained by the transition between GPS and MPS sources coinciding with the jet breakout phase of PS sources, such that GPS sources traverse through the surrounding medium at a lower speed than MPS sources. Such evolution has been observed in simulations of PS source evolution.

We present a new perturbative full-shape analysis of BOSS galaxy clustering data, including the full combination of the galaxy power spectrum and bispectrum multipoles, baryon acoustic oscillations, and cross-correlations with the gravitational lensing of cosmic microwave background measured from \textit{Planck}. Assuming the $\Lambda$CDM model, we constrain the matter density fraction $\Omega_m = 0.3154\pm 0.0089$, the Hubble constant $H_0=68.34\pm 0.77\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, and the mass fluctuation amplitude $\sigma_8=0.686\pm 0.027$ (equivalent to $S_8 = 0.704\pm 0.031$). Cosmic structure at low redshifts appears suppressed with respect to the Planck $\Lambda$CDM concordance model at $4.5\sigma$. We explore whether this tension can be explained by the recent DESI preference for dynamical dark energy (DDE): the BOSS data combine with DESI BAO and PantheonPlus supernovae competitively compared to the CMB, yielding no preference for DDE, but the same $\sim 10\%$ suppression of structure, with dark energy being consistent with a cosmological constant at 68\% CL. Our results suggest that either the data contains residual systematics, or more model-building efforts may be required to restore cosmological concordance.

The High-Altitude Water Cherenkov Telescope (HAWC) has detected TeV halos associated with two nearby pulsars/pulsar wind nebulae (PWN) -- Geminga and B0656+14. These TeV halos extend up to tens of pc from the central accelerators, indicating that the diffusion of ultrarelativistic electrons and positrons in the interstellar medium has been suppressed by two orders of magnitude. Although Geminga and B0656+14 are at similar distances and in the same field of view, they have distinct histories. Notably, B0656+14 probably still resides within its parent supernova remnant, the Monogem Ring, which can be observed in X-rays. In this work, we perform high-resolution simulations of the propagation and emission of relativistic lepton pairs around B0656+14 using a two-zone diffusion model using the GALPROP numerical code. We compared the predicted inverse-Compton spectrum to the observations made by HAWC and Fermi-LAT and found physically plausible model parameters that resulted in a good fit to the data. Additionally, we estimated the contribution of this TeV-halo to the positron flux observed on Earth and found it to be smaller than 10\% of the measured flux. We conclude that future observations of the TeV halo and its synchrotron emission counterpart in radio and X-ray frequencies will be crucial to distinguish between various possible models.

The diatomic metal monoxides whose optical spectra define the classification of AGB stars along the sequence M-MS-S-SC to carbon stars, that is, TiO, ZrO, LaO and YO, have the unusual property that their ionization energy is below their dissociation limit. The cations of these metal monoxides can be efficiently produced via associative ionization of their constituent ground state atoms and are long-lived. We present a simple model that can explain the observed relative abundance of these metal oxides as a function of the C/O ratio.

Carbon fusion is important to understand the late stages in the evolution of a massive star. Astronomically interesting energy ranges for the 12C+12C reactions have been, however, poorly constrained by experiments. Theoretical studies on stellar evolution have relied on reaction rates that are extrapolated from those measured in higher energies. In this work, we update the carbon fusion reaction rates by fitting the astrophysical S-factor data obtained from direct measurements based on the Fowler, Caughlan, & Zimmerman (1975) formula. We examine the evolution of a 20 M_sun star with the updated 12C+12C reaction rates performing simulations with the MESA (Modules for Experiments for Stellar Astrophysics) code. Between 0.5 and 1 GK, the updated reaction rates are 0.35 to 0.5 times less than the rates suggested by Caughlan and Fowler (1988). The updated rates result in the increase of core temperature by about 7% and of the neutrino cooling by about a factor of three. Moreover, the carbon-burning lifetime is reduced by a factor of 2.7. The updated carbon fusion reaction rates lead to some changes in the details of the stellar evolution model, their impact seems relatively minor compared to other uncertain physical factors like convection, overshooting, rotation, and mass-loss history. The astrophysical S-factor measurements in lower energies have large errors below the Coulomb barrier. More precise measurements in lower energies for the carbon burning would be useful to improve our study and to understand the evolution of a massive star.

Pulsar timing array experiments have recently uncovered evidence for a nanohertz gravitational wave background by precisely timing an ensemble of millisecond pulsars. The next significant milestones for these experiments include characterizing the detected background with greater precision, identifying its source(s), and detecting continuous gravitational waves from individual supermassive black hole binaries. To achieve these objectives, generating accurate and precise times of arrival of pulses from pulsar observations is crucial. Incorrect polarization calibration of the observed pulsar profiles may introduce errors in the measured times of arrival. Further, previous studies (e.g., van Straten 2013; Manchester et al. 2013) have demonstrated that robust polarization calibration of pulsar profiles can reduce noise in the pulsar timing data and improve timing solutions. In this paper, we investigate and compare the impact of different polarization calibration methods on pulsar timing precision using three distinct calibration techniques: the Ideal Feed Assumption (IFA), Measurement Equation Modeling (MEM), and Measurement Equation Template Matching (METM). Three NANOGrav pulsars-PSRs J1643$-$1224, J1744$-$1134, and J1909$-$3744-observed with the 800 MHz and 1.5 GHz receivers at the Green Bank Telescope (GBT) are utilized for our analysis. Our findings reveal that all three calibration methods enhance timing precision compared to scenarios where no polarization calibration is performed. Additionally, among the three calibration methods, the IFA approach generally provides the best results for timing analysis of pulsars observed with the GBT receiver system. We attribute the comparatively poorer performance of the MEM and METM methods to potential instabilities in the reference noise diode coupled to the receiver and temporal variations in the profile of the reference pulsar, respectively.

Trajectories of photons of cosmic microwave background (CMB) from the surface of last scattering to us could be deflected by extremely low frequency primordial gravitational wave (PGW). With large scale structure (LSS) producing a smoothing of the acoustic peaks in the power spectrum of the CMB anisotropies through weak lensing, the presence of extremely low frequency PGW could enhance the effect of weak lensing on CMB due to the coupling of extremely low frequency PGW and LSS, thus, give rise to much more smoothing of the spectrum. This may be an natural explanation for the lensing amplitude anomaly observed by Planck, meaning that lensing amplitude anomaly may be the evidence of extremely low frequency PGW.

We consider the three-dimensional rotating motions of neutron stars blown by the "axion wind". Neutron star precession and spin can change from the magnetic moment coupling to the oscillating axion background field, in analogy to the gyroscope motions with a driving force and the laboratory Nuclear Magnetic Resonance(NMR) detections of the axion. This effect modulates the pulse arrival time of the pulsar timing arrays. It shows up as a signal on the timing residual and two-point correlation function on the recent data of Nanograv and PPTA. The current measurement of PTAs can thus cast constraints on the axion-nucleon coupling as g_{ann} ~ 10^{-12}{GeV}^{-1}.

We examine deep optical images of edge-on galaxies selected from the Sloan Digital Sky Survey (SDSS) Stripe\,82. The entire sample consists of over 800 genuine edge-on galaxies with spectroscopic redshifts out to $z\sim0.2$. To discern the faintest details around the galaxies, we use three different data sources with a photometric depth of down to 30 mag\,arcsec$^{-2}$ in the $r$ band: SDSS Stripe\,82, Hyper Suprime-Cam Strategic Program, and DESI Legacy Imaging Surveys. Our analysis of the deep images reveals a variety of low surface brightness features. 49 galaxies exhibit prominent tidal structures, including tidal tails, stellar streams, bridges, and diffuse shells. Additionally, 56 galaxies demonstrate peculiar structural features such as lopsided discs, faint warps, and dim polar rings. Overall, we detect low surface brightness structures in 94 galaxies out of 838, accounting for 11\% of the sample. Notably, the fraction of tidal structures is only 5.8\%, which is significantly lower than that obtained in modern cosmological simulations and observations. Previous studies have shown that strongly interacting galaxies have stellar discs about 1.5--2 times thicker than those without apparent interactions. In an analysis where tidal features are carefully masked for precise disc axis ratio measurements, we show that discs of galaxies with tidal features are 1.33 times thicker, on average, than control galaxies that do not have visible tidal features. Furthermore, we find that edge-on galaxies with tidal structures tend to have a higher fraction of oval and boxy discs than galaxies without tidal features.

Since a few years, we have finally entered the era of discoveries of multiply-imaged gravitationally lensed supernovae. To date, all cluster lensed supernovae have been found from space, while those deflected by individual galaxies were identified with wide-field ground-based surveys through the magnification of "standard candles" method, i.e., without the need of spatially resolving the individual images. We review the challenges in identifying these extremely rare events, as well as the unique opportunities they offer for time-delay cosmography and the study of the properties of the deflecting bodies acting as lenses.

We investigate the correlation between stellar mass (M*) and star formation rate (SFR) across the stellar mass range log10(M*/Msun)~6-11. We consider almost 50,000 star-forming galaxies at z~3-7, leveraging data from COSMOS/SMUVS, JADES/GOODS-SOUTH, and MIDIS/XDF. This is the first study spanning such a wide stellar mass range without relying on gravitational lensing effects. We locate our galaxies on the SFR-M* plane to assess how the location of galaxies in the star-formation main sequence (MS) and starburst (SB) region evolves with stellar mass and redshift. We find that the two star-forming modes tend to converge at log10(M*/Msun) < 7, with all galaxies found in the SB mode. By dissecting our galaxy sample in stellar mass and redshift, we show that the emergence of the star-formation MS is stellar-mass dependent: while in galaxies with log10(M*/Msun) > 9 the MS is already well in place at z = 5-7, for galaxies with log10(M*/Msun)~7-8 it only becomes significant at z<4. Overall, our results are in line with previous findings that the SB mode dominates amongst low stellar-mass galaxies. The earlier emergence of the MS for massive galaxies is consistent with galaxy downsizing.

The search for the polarized imprint of primordial gravitational waves in the cosmic microwave background (CMB) as direct evidence of cosmic inflation requires exquisite sensitivity and control over systematics. The next-generation CMB-S4 project intends to improve upon current-generation experiments by deploying a significantly greater number of highly-sensitive detectors, combined with refined instrument components based on designs from field-proven instruments. The Precursor Small Aperture Telescope (PreSAT) is envisioned as an early step to this next generation, which will test prototype CMB-S4 components and technologies within an existing BICEP Array receiver, with the aim of enabling full-stack laboratory testing and early risk retirement, along with direct correlation of laboratory component-level performance measurements with deployed system performance. The instrument will utilize new 95/155GHz dichroic dual-linear-polarization prototype detectors developed for CMB-S4, cooled to 100mK via the installation of an adiabatic demagnetization refrigerator, along with a prototype readout chain and prototype optics manufactured with wide-band anti-reflection coatings. The experience gained by integrating, deploying, and calibrating PreSAT will also help inform planning for CMB-S4 small aperture telescope commissioning, calibration, and operations well in advance of the fabrication of CMB-S4 production hardware.

We perform a large-scale hydrodynamic simulation of a massive star cluster whose stellar population mimics that of the Cygnus OB2 association. The main-sequence stars are first simulated during 1.6 Myr, until a quasi-stationary state is reached. At this time the three Wolf-Rayet stars observed in Cygnus OB2 are added to the simulation, which continues to 2 Myr. Using a high-resolution grid in the centre of the domain, we can resolve the most massive stars individually, which allows us to probe the kinetic structures at small (parsec) scales. We find that, although the cluster excavates a spherical "superbubble" cavity, the stellar population is too loosely distributed to blow a large-scale cluster wind termination shock, and that collective effects from wind-wind interactions are much less efficient than usually assumed. This challenges our understanding of the ultra-high energy emission observed from the region.

It was known that an ideal spherically symmetric stellar system with isotropic velocities and an inner core cannot reside in a Navarro, Frenk, and White (NFW) gravitational potential. The incompatibility can be pinned down to the radial gradient of the NFW potential in the very center of the system, which differs from zero. The gradient is identically zero in an Einasto potential, also an alternative representation of the dark matter (DM) halos created by the kind of cold DM (CDM) defining the current cosmological model. Here we show that, despite the inner gradient being zero, stellar cores are also inconsistent with Einasto potentials. This result may have implications to constrain the nature of DM through interpreting the stellar cores often observed in dwarf galaxies.

[Abridged] ESO 141-G55 is a nearby X-ray bright BLS1, which has been classified as a bare AGN due to the lack of warm absorption along its line-of-sight, providing an unhampered view into its disc-corona system. We aim to probe its disc-corona system thanks to the first simultaneous XMM-Newton and NuSTAR observation obtained on October 1-2, 2022. We carry out the X-ray broadband spectral analysis to determine the dominant process(es) at work, as well as the SED analysis to determine the disc-corona properties. The simultaneous broadband X-ray spectrum of ESO 141-G55 is characterised by the presence of a prominent smooth soft X-ray excess, a broad Fe K emission line and a significant Compton hump. The RGS spectra confirmed the lack of intrinsic warm-absorbing gas along our line of sight in the AGN rest frame, confirming that it is still in a bare state. However, soft X-ray emission lines are observed indicating substantial warm gas out of our line of sight. The intermediate inclination of the disc-corona system, ~43°, may offer us a favourable configuration to observe UFOs from the disc, but none is found in this 2022 observation, contrary to a previous 2007 XMM-Newton one. Relativistic reflection alone on a standard disc is ruled out from the X-ray broadband analysis, while a combination of soft and hard Comptonisation by a warm and hot corona (relagn), plus relativistic reflection (reflkerrd) reproduces its SED quite well. The hot corona temperature is very hot, ~140 keV, much higher than about 80% of the AGNs, whereas the warm corona temperature, ~0.3 keV, is similar to the values found in other sub-Eddington AGNs. ESO 141-G55 is accreting at a moderate Eddington accretion rate (~10--20%). Our analysis points to a significant contribution of an optically-thick warm corona to both the soft X-ray and UV emission in ESO 141-G55.

Despite being pivotal to the habitability of our planet, the process by which Earth gained its present-day hydrogen budget is unclear. Due to their isotopic similarity to terrestrial rocks across a range of elements, enstatite chondrites (ECs) are thought to be the meteorites that best represent Earth's building blocks. Because of ECs' nominally anhydrous mineralogy, these building blocks have long been presumed to have supplied negligible hydrogen to the proto-Earth. Instead, hydrogen has been proposed to have been delivered to our planet after its main stage of formation by impacts from hydrated asteroids. In this case, our planet's habitability would have its origins in a stochastic process. However, ECs have recently been found to unexpectedly contain enough hydrogen to readily explain Earth's present-day water budget. Although this result would transform the processes we believe are required for rocky planets to be suitable to life, the mineralogical source of ~80% of hydrogen in these meteorites was previously unknown. As such, the reason ECs are seemingly rich in hydrogen was unclear. Here, we apply sulfur X-ray absorption near edge structure (S-XANES) spectroscopy to ECs, finding that most (~70%) of their hydrogen is bonded to sulfur. Moreover, the concentration of the S-H bond is intimately linked to the abundance of micrometre-scale pyrrhotite (Fe1-xS, 0<x<0.125), suggesting most hydrogen in these meteorites is carried in this phase. These findings elucidate the presence of hydrogen in Earth's building blocks, providing the key evidence that unlocks a systematic, rather than stochastic, origin of Earth's hydrogen.

Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and K-mouflage. In particular, we compare the predictions from full $N$-body simulations, two $N$-body codes with approximate time integration schemes, a parametrised modified $N$-body implementation and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in $2\%$ agreement for $z\leq1$ and $k\leq 1~h/{\rm Mpc}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.

The Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2) flew on May 13$^{\text{th}}$ and 14$^{\text{th}}$ of 2023. Consisting of two novel optical telescopes, the payload utilized next-generation instrumentation for the observations of extensive air showers from near space. One instrument, the fluorescence telescope (FT) searched for Ultra-High Energy Cosmic Rays (UHECRs) by recording the atmosphere below the balloon in the near-UV with a 1~$\mu$s time resolution using 108 multi-anode photomultiplier tubes with a total of 6,912 channels. Validated by pre-flight measurements during a field campaign, the energy threshold was estimated around 2~EeV with an expected event rate of approximately 1 event per 10 hours of observation. Based on the limited time afloat, the expected number of UHECR observations throughout the flight is between 0 and 2. Consistent with this expectation, no UHECR candidate events have been found. The majority of events appear to be detector artifacts that were not rejected properly due to a shortened commissioning phase. Despite the earlier-than-expected termination of the flight, data were recorded which provide insights into the detectors stability in the near-space environment as well as the diffuse ultraviolet emissivity of the atmosphere, both of which are impactful to future experiments.

In multi-spectral images made by Earth observation satellites that use push-broom scanning, such as those operated by Planet Labs Corp., moving objects can be identified by the appearance of the object at a different locations in each spectral band. The apparent velocity can be measured if the relative acquisition time between images in different spectral bands is known to millisecond accuracy. The images in the Planet Labs archive are mosaics of individual exposures acquired at different times. Thus there is not a unique acquisition time for each spectral band. In an earlier paper, we proposed a method to determine the relative acquisition times from the information in the images themselves. High altitude balloons provide excellent targets to test our proposed method because of their high apparent velocity due to the orbital velocity of the satellite and geometric parallax in images aligned to the level of the ground. We use images of the Chinese balloon that crossed the US in February, 2024 as well as images of an identical balloon over Colombia to test our method. Our proposed method appears to be successful and allows the measurement of the apparent velocity of moving objects from the information available in the archive.

Fast radio bursts (FRBs) are brilliant short-duration flashes of radio emission originating at cosmological distances. The vast diversity in the properties of currently known FRBs, and the fleeting nature of these events make it difficult to understand their progenitors and emission mechanism(s). Here we report high time resolution polarization properties of FRB 20210912A, a highly energetic event detected by the Australian Square Kilometre Array Pathfinder (ASKAP) in the Commensal Real-time ASKAP Fast Transients (CRAFT) survey, which show intra-burst PA variation similar to Galactic pulsars and unusual variation of Faraday Rotation Measure (RM) across its two sub-bursts. The observed intra-burst PA variation and apparent RM variation pattern in FRB 20210912A may be explained by a rapidly-spinning neutron star origin, with rest-frame spin periods of ~1.1 ms. This rotation timescale is comparable to the shortest known rotation period of a pulsar, and close to the shortest possible rotation period of a neutron star. Curiously, FRB 20210912A exhibits a remarkable resemblance with the previously reported FRB 20181112A, including similar rest-frame emission timescales and polarization profiles. These observations suggest that these two FRBs may have similar origins.

In this study, we present a method for detecting and analyzing the velocities of moving objects in Earth observation satellite images, specifically using data from Planet Labs' push broom scanning satellites. By exploiting the sequential acquisition of multi-spectral images, we estimate the relative differences in acquisition times between spectral bands. This allows us to determine the velocities of moving objects, such as aircraft, even without precise timestamp information from the image archive. We validate our method by comparing the velocities of aircraft observed in satellite images with those reported by onboard ADS-B transponders. The results demonstrate the potential, despite challenges posed by proprietary data limitations, of a new, useful application of commercial satellite data originally intended as an ongoing, once-daily survey of single images covering the entire land-area of the Earth. Our approach extends the applicability of satellite survey imagery for dynamic object tracking and contributes to the broader use of commercial satellite data in scientific research.

I compare the dark matter content within stellar half-mass radius expected in a $\Lambda$CDM-based galaxy formation model with existing observational estimates for the observed dwarf satellites of the Milky Way and ultra-diffuse galaxies (UDGs). The model reproduces the main properties and scaling relations of dwarf galaxies, in particular their stellar mass-size relation. I show that the model also reproduces the relation between the dark matter mass within the half-mass radius, $M_{\rm dm}(<r_{1/2})$, and stellar mass exhibited by the observed dwarf galaxies. The scatter in the $M_{\rm dm}(<r_{1/2})-M_\star$ relation is driven primarily by the broad range of sizes of galaxies of a given stellar mass. I also show the $M_{\rm dm}(<r_{1/2})$ of UDGs are within the range expected in the model for their stellar mass, but they tend to lie above the median relation due to their large sizes. The upper limits on $M_{\rm dm}(<r_{1/2})$ for the dark matter deficient UDGs are also consistent with the range of dark matter masses expected in the model. The most dark matter-deficient galaxies of a given size correspond to halos with the smallest concentrations and the largest ratios of $M_\star/M_{\rm 200c}$. Conversely, the most dark matter-dominated galaxies are hosted by the highest concentration halos with the smallest $M_\star/M_{\rm 200c}$ ratios. The model indicates that the scatter between $M_{\rm dm}(<r_{1/2})$ and $M_{\rm 200c}$ is large, which renders inference of the virial mass from $M_{\rm dm}(<r_{1/2})$ uncertain and dependent on specific assumptions about the halo mass profile. Results presented in this paper indicate that dark matter-deficient UDGs may represent a tail of the expected dark matter profiles, especially if the effect of feedback on these profiles is taken into account.

We identify voids as maximal non-overlapping spheres within the haloes of the Uchuu simulation and three smaller halo simulation boxes with smaller volume and different $\sigma_{8}$ values, and galaxies with redshift in the range $0.02<z<0.132$ and absolute magnitude in the $r-$band $M_{r}<-20.5$ of 32 Uchuu-SDSS simulated lightcones the seventh release of \textit{The Sloan Digital Sky Survey} (SDSS DR7) survey. We compute the Void Probability Function and the abundance of voids larger than $r$ predicted by the theoretical framework used in this work and we check that it predicts successfully both void functions for the halo simulation boxes. Next, we asses the potential of this theoretical framework to constrain cosmological parameters using Uchuu-SDSS void statistics, and we calculate the confidence levels using Monte Carlo Markov Chain techniques to infer the values of $\sigma_{8}$, $\Omega_{\rm m}$ and H$_{0}$ from the SDSS sample used. The constraints we obtain from the SDSS survey sample used. The results are: $\sigma_{8}=1.028^{+0.273}_{-0.305}$, $\Omega_{\rm m}=0.296^{+0.110}_{-0.102}$, H$_{0}=83.43\pm^{+29.27}_{-27.70}$, $\Gamma=0.1947^{+0.0578}_{-0.0516}$ and S$_{8}$=1.017$^{+0.363}_{-0.359}$. If we combine these constraints with KiDS-1000+DESY3, we get $\sigma_{8}=0.858^{+0.040}_{-0.040}$, $\Omega_{\rm m}=0.257\pm^{+0.023}_{-0.020}$, H$_{0}=74.17^{+4.66}_{-4.66}$ and S$_{8}$=0.794$^{+0.016}_{-0.016}$. The combined uncertainties are approximately a factor 2-3 smaller than only-Weak-Lensing uncertainties. This is a consequence of the orientation of the confidence level contours of SDSS voids and Weak Lensing in the plane $\sigma_{8}-\Omega_{\rm m}$, which are almost orthogonal (abridged).

Superthin galaxies are observed to have stellar disks with extremely small minor-to-major axis ratios. In this work, we investigate the formation of superthin galaxies in the TNG100 simulation. We trace the merger history and investigate the evolution of galaxy properties of a selected sample of superthin galaxies and a control sample of galaxies that share the same joint probability distribution in the stellar-mass and color diagram. Through making comparisons between the two galaxy samples, we find that present-day superthin galaxies had similar morphologies as the control sample counterparts at higher redshifts, but have developed extended flat `superthin' morphologies since $z \sim 1$. During this latter evolution stage, superthin galaxies undergo overwhelmingly higher frequency of prograde mergers (with orbit-spin angle $\theta_{\rm orb} \leqslant 40^\circ$). Accordingly the spins of their dark matter halos have grown significantly and become noticeably higher than that of their normal disk counterparts. This further results in the buildup of their stellar disks at larger distances much beyond the regimes of normal disk galaxies. We also discuss the formation scenario of those superthin galaxies that live in larger dark matter halos as satellite galaxies therein.

Gravitational waves (GWs) are signals that propagate across large distances in the Universe, and thus, they bring information on the cosmic history. GW sources are at the same time distance indicators and tracers of the matter field. Events generated by binary systems can be divided into bright standard sirens, when followed by electromagnetic transients from which the redshift of the source can be measured, and the more numerous dark standard sirens, when counterparts are not available. In this proceeding, I will discuss some methods for testing the cosmological model using either bright or dark sirens and their combinations with other cosmological probes, focusing on some of my own recent contributions.

Average magnetic field measurements are presented for 62 M-dwarf members of the Pleiades open cluster, derived from Zeeman-enhanced Fe I lines in the H-band. An MCMC methodology was employed to model magnetic filling factors using SDSS-IV APOGEE high-resolution spectra, along with the radiative transfer code SYNMAST, MARCS stellar atmosphere models, and the APOGEE DR17 spectral line list. There is a positive correlation between mean magnetic fields and stellar rotation, with slow-rotator stars (Rossby number, Ro$>$0.13) exhibiting a steeper slope than rapid-rotators (Ro$<$0.13). However, the latter sample still shows a positive trend between Ro and magnetic fields, which is given by $<$B$>$ = 1604 $\times$ Ro$^{-0.20}$. The derived stellar radii, when compared with physical isochrones, show that on average, our sample shows radius inflation, with median enhanced radii ranging from +3.0$\%$ to +7.0$\%$, depending on the model. There is a positive correlation between magnetic field strength and radius inflation, as well as with stellar spot coverage, correlations that together indicate that stellar spot-filling factors generated by strong magnetic fields might be the mechanism that drives radius inflation in these stars. We also compare our derived magnetic fields with chromospheric emission lines (H$\alpha$, H$\beta$ and Ca II K), as well as with X-ray and H$\alpha$ to bolometric luminosity ratios, and find that stars with higher chromospheric and coronal activity tend to be more magnetic.

Blazars constitute the most numerous source class in the known extragalactic population of very high energy (VHE) gamma-ray sources. However, determining their redshifts is often challenging due to weak or non-existent emission lines in their spectra. This study focuses on two BL Lacs, KUV 00311-1938 and S2 0109+22, where previous attempts at redshift determination have faced difficulties. By combining spectroscopic observations with photometric redshift estimates, we tentatively assign a redshift of z = 0.634 to KUV 00311-1938 and a likely redshift of z = 0.49 to S2 0109+22. Establishing redshift estimates for high-redshift blazars is crucial for understanding extragalactic VHE gamma-ray sources and their interactions with the surrounding universe.

We present the first high-resolution, high-frequency radio continuum survey that fully maps an extragalactic deep field: the 10GHz survey of the Great Observatories Origins Deep Survey-North (GOODS-N) field. This is a Large Program of the Karl G. Jansky Very Large Array that allocated 380 hours of observations using the X-band ($8-12$GHz) receivers, leading to a 10GHz mosaic of the GOODS-field with an average rms noise $\sigma_{\rm n}=671\,\rm nJy\,beam^{-1}$ and angular resolution $\theta_{1/2}=0.22$arcsec across 297$\rm arcmin^2$. To maximize the brightness sensitivity we also produce a low-resolution mosaic with $\theta_{1/2}=1.0$arcsec and $\sigma_{\rm n}=968\,\rm nJy\,beam^{-1}$, from which we derive our master catalog containing 256 radio sources detected with peak signal-to-noise ratio $\geq 5$. Radio source size and flux density estimates from the high-resolution mosaic are provided in the master catalog as well. The total fraction of spurious sources in the catalog is 0.75%. Monte Carlo simulations are performed to derive completeness corrections of the catalog. We find that the 10GHz radio source counts in the GOODS-N field agree, in general, with predictions from numerical simulations/models and expectations from 1.4 and 3GHz radio counts.

Luminous Fast Blue Optical Transients (LFBOTs) are a class of extragalactic transients notable for their rapid rise and fade times, blue colour and accompanying luminous X-ray and radio emission. Only a handful have been studied in detail since the prototypical example AT2018cow. Their origins are currently unknown, but ongoing observations of previous and new events are placing ever stronger constraints on their progenitors. We aim to put further constraints on the LFBOT AT2023fhn, and LFBOTs as a class, using information from the multi-wavelength transient light-curve, its host galaxy and local environment. Our primary results are obtained by fitting galaxy models to the spectral energy distribution of AT2023fhn's host and local environment, and by modelling the radio light-curve of AT2023fhn as due to synchrotron self-absorbed emission from an expanding blast-wave in the circumstellar medium. We find that the neither the host galaxy nor circumstellar environment of AT2023fhn are unusual compared with previous LFBOTs, but that AT2023fhn has a much lower X-ray to ultraviolet luminosity ratio than previous events. We argue that the variety in ultraviolet-optical to X-ray luminosity ratios among LFBOTs is likely due to viewing angle differences, and that the diffuse, yet young local environment of AT2023fhn - combined with a similar circumstellar medium to previous events - favours a progenitor system containing a massive star with strong winds. Plausible progenitor models in this interpretation therefore include black hole/Wolf-Rayet mergers or failed supernovae.

CMB-S4, the next-generation CMB observatory, will deploy hundreds of thousands of detectors to enable mapping the millimeter-wavelength sky with unprecedented speed. The large aperture telescopes for CMB-S4 consist of six-meter diameter crossed Dragone designs and a five-meter diameter three-mirror anastigmat. The two-mirror crossed Dragone design requires astigmatism corrections in the refractive optics to achieve diffraction-limited performance. We present biconic lens corrections for the CMB-S4 crossed Dragone camera optics and compare these designs to the camera optics for the three mirror anastigmat, as the optical designs of the cameras for these telescopes are being prototyped.

We calculate the extent to which collisionless dark matter impacts the stability of supermassive stars $(M\gtrsim10^4\,M_\odot)$. We find that, depending on the star's mass, a dark matter content in excess of ${\sim}1\%$ by mass throughout the entire star can raise the critical central density for the onset general relativistic instability, in some cases by orders of magnitude. We consider implications of this effect for the onset of nuclear burning and significant neutrino energy losses.

Volume-phase holographic gratings (VPHGs) are widely used in astronomical spectrographs due to their adaptability and high diffraction efficiency. Most VPHGs in operation use dichromated gelatin as a recording material, whose performance is sensitive to the coating and development process, especially in the near-UV. In this letter, we present the characterization of two UV-blue VPHG prototypes for the BlueMUSE integral field spectrograph on the VLT, based on dichromated gelatin and the Bayfol$\circledR$HX photopolymer film as recording materials. Our measurements show that both prototypes meet the required diffraction efficiency and exhibit similar performance with a wavelength-average exceeding 70% in the 350-580 nm range. Deviations from theoretical models increase towards 350 nm, consistently with previous studies on similar gratings. We also report similar performances in terms spatial uniformity and grating-to-grating consistency. Likewise, no significant differences in wavefront error or scattered light are observed between the prototypes.

BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its specific capabilities. The BlueMUSE consortium includes 9 institutes located in 7 countries and is led by the Centre de Recherche Astrophysique de Lyon (CRAL). The BlueMUSE project development is currently in Phase A, with an expected first light at the VLT in 2031. We introduce here the Top Level Requirements (TLRs) derived from the main science cases, and then present an overview of the BlueMUSE system and its subsystems fulfilling these TLRs. We specifically emphasize the tradeoffs that are made and the key distinctions compared to the MUSE instrument, upon which the system architecture is built.

BlueMUSE is a blue, medium spectral resolution, panoramic integral-field spectrograph under development for the Very Large Telescope (VLT). We demonstrate and discuss an early End-To-End simulation software for final BlueMUSE datacube products. Early access to such simulations is key to a number of aspects already in the development stage of a new major instrument. We outline the software design choices, including lessons learned from the MUSE instrument in operation at the VLT since 2014. The current simulation software package is utilized to evaluate some of the technical specifications of BlueMUSE as well as giving assistance in the assessment of certain trade offs regarding instrument capabilities, e.g., spatial and spectral resolution and sampling. By providing simulations of the end-user product including realistic environmental conditions such as sky contamination and seeing, BlueSi can be used to devise and prepare the science of the instrument by individual research teams.

Full spectrum fitting is a powerful tool for estimating the stellar populations of galaxies, but the fitting results are often significantly influenced by internal dust attenuation. For understanding how the choice of the internal dust correction method affects the detailed stellar populations estimated from the full spectrum fitting, we analyze the Sydney-Australian Astronomical Observatory Multi-object Integral field spectrograph (SAMI) galaxy survey data using the Penalized PiXel-Fitting (PPXF) package. Three choices are compared: (Choice-1) using the PPXF reddening option, (Choice-2) using the multiplicative Legendre polynomial, and (Choice-3) using none of them (no dust correction). In any case, the total mean stellar populations show reasonable mass-age and mass-metallicity relations (MTR and MZR), although the correlations appear to be strongest for Choice-1 (MTR) and Choice-2 (MZR). When we compare the age-divided mean stellar populations, the MZR of young (< 10^9.5 yr ~ 3.2 Gyr) stellar components in Choice-2 is consistent with the gas-phase MZR, whereas those in the other two choices hardly are. On the other hand, the MTR of old (>= 10^9.5 yr) stellar components in Choice-1 seems to be more reasonable than that in Choice-2, because the old stellar components in low-mass galaxies tend to be relatively younger than those in massive galaxies. Based on the results, we provide empirical guidelines for choosing the optimal options for dust correction.

We study the anisotropy of centroid and integrated intensity maps with synthetic observations. We perform post-process radiative transfer including the optically thick regime that was not covered in Hernández-Padilla et al. (2020). We consider the emission in various CO molecular lines, that range from optically thin to optically thick ($\mathrm{^{12}CO}$, $\mathrm{^{13}CO}$, $\mathrm{C^{18}O}$, and $\mathrm{C^{17}O}$). The results for the velocity centroids are similar to those in the optically thin case. For instance, the anisotropy observed can be attributed to the Alfvén mode, which dominates over the slow and fast modes when the line of sight is at a high inclination with respect to the mean magnetic field. A few differences arise in the models with higher opacity, where some dependence on the sonic Mach number becomes evident. In contrast to the optically thin case, maps of integrated intensity become more anisotropic in optically thick lines. In this situation the scales probed are restricted, due to absorption, to smaller scales which are known to be more anisotropic. We discuss how the sonic Mach number can affect the latter results, with highly supersonic cases exhibiting a lower degree of anisotropy.

The source of the Galactic Lithium (Li) has long been a puzzle. With the discovery of Li in novae, extensive research has been conducted. However, there still exists a significant disparity between the observed abundance of lithium in novae and the existing theoretical predictions. Using the Modules for Experiments in Stellar Astrophysics (MESA), we simulate the evolution of nova with element diffusion and appropriately increased the amount of 3^He in the mixtures. Element diffusion enhances the transport efficiency between the nuclear reaction zone and the convective region on the surface of the white dwarf during nova eruptions, which results in more 7^Be to be transmitted to the white dwarf surface and ultimately ejected. Compared to the previous predictions, the abundance of 7^Be in novae simulated in our model significantly increases. And the result is able to explain almost all observed novae. Using the method of population synthesis, we calculate Li yield in the Galaxy. We find that the Galactic occurrence rate of nova is about 130 yr^{-1}, and about 110M Li produced by nova eruption is ejected into the interstellar medium (ISM). About 73\% of Li in the Galactic ISM originates from novae, and approximately 15\%-20\% of the entire Galaxy. It means that novae are the important source of Li in the Galactic.

As the hottest high-energy neutrino spot, NGC 1068 has received much attention in recent years. Here we focus on the central region of the active galactic nuclei (AGN) and propose an outflow-cloud interaction model that could probably explain the observed neutrino data. Considering the accretion process adjacent to the central supermassive black hole (SMBH) of NGC 1068, strong outflows will be generated, which will likely interact with surrounding clouds floating in the corona region. Particles carried by the outflow will be accelerated to very high energy by the shocks forming during the outflow-cloud interactions. For the accelerated high-energy protons, $p\gamma$ interactions with the background photon field of the corona and disk and $pp$ interaction with the surrounding gas will produce considerable high-energy $\gamma$-rays and neutrino. However, because of the extremely dense photon fields in the corona and disk, the newly generated $\gamma$-rays will be significantly attenuated through the $\gamma\gamma$ absorptions. In our scenario, the expected GeV-TeV $\gamma$-ray emission will be suppressed to a much lower level than the neutrino emission, consistent with the observational characteristics of NGC 1068, while the generated 1-30\,TeV neutrino flux can fit the IceCube data very well.

Based on the Gluon Condensation (GC) model, the relationship between the spectra of electrons, $\gamma$ rays, and neutrinos in cosmic rays can be deduced. It has been found that these particles share the same parameter, $\beta_p$, and have an identical GC threshold values. This paper explores the connection between the second excess spectra of electron and the spectra of gamma rays and neutrinos. According to the observed gamma-ray data, it is suggested that the source LHAASO J2108+5157 might contribute to the second excess of electron.

This study proposes an analytical framework for deriving the surface brightness profile and geometry of a geometrically-thin axisymmetric disc from interferometric observation of continuum emission. Such precise modelling facilitates the exploration of faint non-axisymmetric structures, such as spirals and circumplanetary discs. As a demonstration, we simulate interferometric observations of geometrically-thin axisymmetric discs. The proposed method can reasonably recover the injected axisymmetric structures, whereas Gaussian fitting of the same data yielded larger errors in disc orientation estimation. To further test the applicability of the method, it was applied to the mock data for $m=1,2$ spirals and a point source, which are embedded in a bright axisymmetric structure. The injected non-axisymmetric structures were reasonably recovered except for the innermost parts, and the disc geometric parameter estimations were better than Gasussian fitting. The method was then applied to the real data of Elias 20 and AS 209, and it adequately subtracted the axisymmetric component, notably in Elias 20, where substantial residuals remained without our method. We also applied our method to continuum data of PDS 70 to demonstrate the effectiveness of the method. We successfully recovered emission from PDS 70 c consistently with previous studies, and also tentatively discovered new substructures. The current formulation can be applied to any data for disc continuum emission, and aids in the search of spirals and circumplanetary discs, whose detection is still limited.

Recent observations of core-collapse supernovae (SNe) suggest aspherical explosions. Globally aspherical structures in SN explosions are regarded as the key for understanding their explosion mechanism. However, the exact explosion geometries from the inner cores to the outer envelopes are poorly understood. Here, we present photometric, spectroscopic and polarimetric observations of the Type IIP SN 2021yja and discuss its explosion geometry, in comparison to those of other Type IIP SNe that show large-scale aspherical structures in their hydrogen envelopes (SNe 2012aw, 2013ej and 2017gmr). During the plateau phase, SNe 2012aw and 2021yja exhibit high continuum polarization characterized by two components with perpendicular polarization angles. This behavior can be interpreted to be due to a bipolar explosion, composed of a polar (energetic) and an equatorial (bulk) components of the SN ejecta. In such a bipolar explosion, an aspherical axis created by the polar ejecta would be dominating at early phases, while the perpendicular axis along the equatorial ejecta would emerge at late phases after the receding of the photosphere in the polar ejecta. The interpretation of the bipolar explosions in SNe 2012aw and 2021yja is also supported by other observational properties, including the time evolution of the line velocities and the line shapes in the nebular spectra. The polarization of other Type IIP SNe that show large-scale aspherical structures in the hydrogen envelope (SNe 2013ej and 2017gmr) is also consistent with the bipolar-explosion scenario, although this is not conclusive.

High resolution spectroscopy of exoplanet atmospheres provides insights into their composition and dynamics from the resolved line shape and depth of thousands of spectral lines. WASP-127 b is an extremely inflated sub-Saturn (R$_\mathrm{p}$= 1.311 R$_\mathrm{Jup}$, M$_\mathrm{p}$= 0.16 M$_\mathrm{Jup}$) with previously reported detections of H$_2$O, CO$_2$, and Na. However, the seeming absence of the primary carbon reservoir expected at WASP-127 b temperatures (T$_{eq}$ $\sim$ 1400 K) from chemical equilibrium, CO, posed a mystery. In this manuscript, we present the analysis of high resolution observations of WASP-127 b with the Immersion GRating INfrared Spectrometer (IGRINS) on Gemini South. We confirm the presence of H$_2$O (8.67 $\sigma$) and report the detection of CO (4.34 $\sigma$). Additionally, we conduct a suite of Bayesian retrieval analyses covering a hierarchy of model complexity and self-consistency. When freely fitting for the molecular gas volume mixing ratios, we obtain super-solar metal enrichment for H$_2$O abundance of log$_{10}$X$_\mathrm{H_2O}$ = --1.23$^{+0.29}_{-0.49}$ and a lower limit on the CO abundance of log$_{10}$X$_\mathrm{CO}$ $\ge$ --2.20 at 2$\sigma$ confidence. We also report a tentative evidence of photochemistry in WASP-127 b based upon the indicative depletion of H$_2$S. This is also supported by the data preferring models with photochemistry over free-chemistry and thermochemistry. The overall analysis implies a super-solar ($\sim$ 39$\times$ Solar; [M/H] = $1.59^{+0.30}_{-0.30}$) metallicity for the atmosphere of WASP-127 b and an upper limit on its atmospheric C/O ratio as $<$ 0.68.

The Simons Observatory (SO) is a group of modern telescopes dedicated to observing the polarized cosmic microwave background (CMB), transients, and more. The Observatory consists of four telescopes and instruments, with over 60,000 superconducting detectors in total, located at ~5,200 m altitude in the Atacama Desert of Chile. During observations, it is important to ensure the detectors, telescope platforms, calibration and receiver hardware, and site hardware are within operational bounds. To facilitate rapid response when problems arise with any part of the system, it is essential that alerts are generated and distributed to appropriate personnel if components exceed these bounds. Similarly, alerts are generated if the quality of the data has become degraded. In this paper, we describe the SO alarm system we developed within the larger Observatory Control System (OCS) framework, including the data sources, alert architecture, and implementation. We also present results from deploying the alarm system during the commissioning of the SO telescopes and receivers.

In this paper, we introduce a Unet model of deep learning algorithms for reconstructions of the 3D peculiar velocity field, which simplifies the reconstruction process with enhanced precision. We test the adaptability of the Unet model with simulation data under more realistic conditions, including the redshift space distortion (RSD) effect and halo mass threshold. Our results show that the Unet model outperforms the analytical method that runs under ideal conditions, with a 16% improvement in precision, 13% in residuals, 18% in correlation coefficient and 27% in average coherence. The deep learning algorithm exhibits exceptional capacities to capture velocity features in non-linear regions and substantially improve reconstruction precision in boundary regions. We then apply the Unet model trained under SDSS observational conditions to the SDSS DR7 data for observational 3D peculiar velocity reconstructions.

We report our identification of three gigaelectronvolt $\gamma$-ray sources, 4FGL J0502.6+0036, 4FGL J1055.9+6507, and 4FGL J1708.2+5519, as Active Galactic Nuclei (AGNs). They are listed in the latest Fermi-LAT source catalog as unidentified ones. We find that the sources all showed $\gamma$-ray flux variations in recent years. Using different survey catalogs, we are able to find a radio source within the error circle of each source's position. Further analysis of optical sources in the fields allows us to determine the optical counterparts, which showed similar variation patterns to those seen in $\gamma$-rays. The optical counterparts have reported redshifts of 0.6, 1.5, and 2.3, respectively, estimated from photometric measurements. In addition, we also obtain an X-ray spectrum of 4FGL J0502.6+0036 and a flux upper limit on the X-ray emission of 4FGL J1055.9+6507 by analyzing the archival data. The broadband spectral energy distributions of the three sources from radio to $\gamma$-rays are constructed. Comparing mainly the $\gamma$-ray properties of the three sources with those of different sub-classes of AGNs, we tentatively identify them as blazars. Followup optical spectroscopy is highly warranted for obtaining their spectral features and thus verifying the identification.

Hypervelocity stars (HVSs) are stars which have been ejected from the Galactic Centre (GC) at velocities of up to a few thousand km/s. They are tracers of the Galactic potential and can be used to infer properties of the GC, such as the initial-mass function and assembly history. HVSs are rare, however, with only about a dozen promising candidates discovered so far. In this work we use a novel, highly efficient method to identify new HVS candidates in Gaia. This method uses the nearly radial trajectories of HVSs to infer their distances and velocities based on their position and Gaia proper motion alone. Through comparison of inferred distances with Gaia parallaxes and photometry we identified 600 HVS candidates with G<20 including the previously discovered S5-HVS1, out of which we obtained ground-based follow-up observations for 196 stars. As we found no new HVSs based on their radial velocity, we used detailed HVS ejection simulations to significantly improve previous HVS ejection rate constraints. In particular, the ejection rate of HVSs more massive than 1 M$_\odot$ cannot be higher than $10^{-5}$ yr$^{-1}$ at $2\sigma$ significance. Additionally, we predict that there are 5-45 unbound HVSs in the complete Gaia catalogue ($1\sigma$ interval), most of which will be main-sequence stars of a few M$_\odot$ at heliocentric distances of tens to hundreds of kpc. By comparing our results to literature HVS candidates, we find an indication of either a time-dependent ejection rate of HVSs or a non-GC origin of many previously identified HVS candidates.

Making observable predictions for cosmic inflation requires determining when the wavenumbers of astrophysical interest today exited the Hubble radius during the inflationary epoch. These instants are commonly evaluated using the slow-roll approximation and measured in e-folds $\Delta N=N-N_\mathrm{end}$, in reference to the e-fold $N_\mathrm{end}$ at which inflation ended. Slow roll being necessarily violated towards the end of inflation, both the approximated trajectory and $N_\mathrm{end}$ are determined at, typically, one or two e-folds precision. Up to now, such an uncertainty has been innocuous, but this will no longer be the case with the forthcoming cosmological measurements. In this work, we introduce a new and simple analytical method, on top of the usual slow-roll approximation, that reduces uncertainties on $\Delta N$ to less than a tenth of an e-fold.

In this paper, we study the long-term (time scale of several years) orbital evolution of lunar satellites under the sole action of natural forces. In particular, we focus on secular resonances, caused either by the influence of the multipole moments of the lunar potential and/or by the Earth's and Sun's third-body effect on the satellite's long-term orbital evolution. Our study is based on a simplified secular model obtained in `closed form' using the same methodology proposed in the recently published report on the semi-analytical propagator of lunar satellite orbits, SELENA. Contrary to the case of artificial Earth satellites, in which many secular resonances compete in dynamical impact, we give numerical evidence that for lunar satellites only the 2 g resonance affects significantly the orbits at secular timescales. We interpret this as a consequence of the strong effect of lunar mascons. We show that the lifetime of lunar satellites is, in particular, nearly exclusively dictated by the 2 g resonance. By deriving a simple analytic model, we propose a theoretical framework which allows for both qualitative and quantitative interpretation of the structures seen in numerically obtained lifetime maps. This involves explaining the main mechanisms driving eccentricity growth in the orbits of lunar satellites. In fact, we argue that the re-entry process for lunar satellites is not necessarily a chaotic process (as is the case for Earth satellites), but rather due to a sequence of bifurcations leading to a dramatic variation in the structure of the separatrices in the 2 g resonance's phase portrait, as we move from the lowest to the highest limit in inclination (at each altitude) where the 2 g resonance is manifested.

Galactic cosmic rays (GCRs) constitute a significant part of the energy budget of our Galaxy, and the study of their accelerators is of high importance in modern astrophysics. Their main sources are likely supernova remnants (SNRs). These objects are capable to convert a part of their mechanical energy into accelerated charged particles. However, even though the mechanical energy reservoir of SNRs is promising, a conversion rate into particle energy of 10 to 20% is necessary to feed the population of GCRs. Such an efficiency is however not guaranteed. Complementary sources deserve thus to be investigated. This communication aims to address the question of the contribution to the acceleration of GCRs by pre-supernova massive stars in binary or higher multiplicity systems

The multiplicity of massive stars is known to be significantly high. Even though the majority of massive stars are located in binary systems, the census of binaries is biased toward shorter periods as longer period systems are more difficult to identify. Alternatively, the search for binary systems with longer periods may proceed differently. As massive binary systems are typically colliding-wind systems, hints for processes occurring in the colliding-wind region could be used as a valuable proxy to identify likely binary systems, and then organize dedicated spectroscopic or interferometric campaigns on a short list of pre-selected targets. In this context, any hint for synchrotron radio emission is seen as a promising indicator of long period binaries, as short period systems undergo severe free-free absorption of the synchrotron emission by the stellar wind material. Usual techniques to identify synchrotron radio emitters constitute thus valid tools to explore that poorly investigated part of the massive binary parameter space. In addition, the identification of a synchrotron emission component in a short period binary can be used as an indicator of the presence of a third companion on a still unrevealed wider orbit in a triple system.

A rather clear problem has remained in black hole physics: localizing black holes. One of the recent theoretical ways proposed to identify black hole mergers' hosts is through multi-messenger gravitational lensing: matching the properties of a lensed galactic host with those of a lensed gravitational wave. This paper reviews the most recent literature and introduces some of the ongoing work on the localization of binary black holes and their host galaxies through lensing of gravitational waves and their electromagnetically-bright hosts.

We propose a novel approach to the detection of point-like sources of high-energy neutrinos. Motivated by evidence for emerging sources in existing data, we focus on the characterisation and interpretation of these sources. The hierarchical Bayesian model is implemented in the Stan platform, enabling computation of the posterior distribution with Hamiltonian Monte Carlo. We simulate a population of weak neutrino sources detected by the IceCube experiment and use the resulting data set to demonstrate and validate our framework. We show that even for the challenging case of sources at the threshold of detection and using limited prior information, it is possible to correctly infer the source properties. Additionally, we demonstrate how modelling flexible connections between similar sources can be used to recover the contribution of sources that would not be detectable individually, going beyond what is possible with existing stacking methods.

Very low-mass stars (those <0.3 solar masses) host orbiting terrestrial planets more frequently than other types of stars, but the compositions of those planets are largely unknown. We use mid-infrared spectroscopy with the James Webb Space Telescope to investigate the chemical composition of the planet-forming disk around ISO-ChaI 147, a 0.11 solar-mass star. The inner disk has a carbon-rich chemistry: we identify emission from 13 carbon-bearing molecules including ethane and benzene. We derive large column densities of hydrocarbons indicating that we probe deep into the disk. The high carbon to oxygen ratio we infer indicates radial transport of material within the disk, which we predict would affect the bulk composition of any planets forming in the disk.

The goal is to compare the intricate details of the magnetic and flow fields around two solar pores, where one is part of an active region and the other is an isolated pore, with a secondary goal of demonstrating the scientific capabilities of the GRIS IFU. Two pores were observed with the HiFI and the GRIS IFU at the GREGOR solar telescope on 29 May and 6 June 2019. The GRIS IFU mosaics provide spectropolarimetric data for inversions of the Ca I 1083.9 nm and Si I 1082.7 nm spectral lines, covering the deep and upper photosphere. The t-SNE machine learning algorithm is employed to identify different classes of Si I Stokes-V profiles. The LCT technique derives horizontal proper motions around the pores. Both pores contain a thin light bridge, are stable during the observations, and never develop a penumbra. The isolated pore is three times smaller and significantly darker than the active-region pore, which is not predicted by simulations. The LCT maps show inflows around both pores, with lower velocities for the isolated pore. Both pores are embedded in the photospheric LOS velocity pattern of the granulation but filamentary structures are only visible in the chromospheric LOS maps of the active-region pore. The t-SNE identifies five clusters of Si I Stokes-V profiles, revealing an `onion-peel' magnetic field structure, despite the small size of the pores. The core with strong vertical magnetic fields is surrounded by concentric layers with lower and more inclined magnetic fields. The active-region pore shows some signatures of increased interaction between plasma motions and magnetic fields, which can be considered as early signs of penumbra formation. However, similar physical properties prevail for smaller pores. A statistically meaningful sample of different pore sizes and morphologies is required to distinguish between active-region and isolated pore formation mechanisms.

Multiple magnetic and kinetic solar wind plasma parameters are used to detect coronal mass ejections (CMEs) as they travel through the heliosphere. There are various interplanetary CME (ICME) catalogues, but due to differences between their ICME identification criteria they can significantly vary. In this paper we analyze Richardson and Cane and CCMC CME Scoreboard ICME catalogues and the SRI RAS solar wind types catalogue, and propose an algorithm of merging them. A unified catalogue is constructed for 2010 to 2022. The resulting catalogue is completed with data from the OMNI database. Analysis of the unified catalogue demonstrated high accuracy when merging events present in multiple catalogues and a tendency of events defined in all three initial catalogues to demonstrate greater duration, speed and geoeffectiveness. The catalog is presented on the SINP MSU Space Weather Exchange website: this https URL

The recent Gamma-Ray Burst (GRB) GRB~211211A provides the earliest ($\sim 5$ h) data of a kilonova (KN) event, displaying bright ($\sim10^{42}$ erg s$^{-1}$) and blue early emission. Previously, this KN has been explained using simplistic multi-component fitting methods. Here, in order to understand the physical origin of the KN emission in GRB~211211A, we employ an analytic multi-zone model for r-process powered KN. We find that r-process powered KN models alone cannot explain the fast temporal evolution and the spectral energy distribution (SED) of the observed emission. Specifically, i) r-process models require high ejecta mass to match early luminosity, which overpredicts late-time emission, while ii) red KN models that reproduce late emission underpredict early luminosity. We propose an alternative scenario involving early contributions from the GRB central engine via a late low-power jet, consistent with plateau emission in short GRBs and GeV emission detected by Fermi-LAT at $\sim10^4$ s after GRB 211211A. Such late central engine activity, with an energy budget of $\sim \text{a few }\%$ of that of the prompt jet, combined with a single red-KN ejecta component, can naturally explain the light curve and SED of the observed emission; with the late-jet -- ejecta interaction reproducing the early blue emission and r-process heating reproducing the late red emission. This supports claims that late low-power engine activity after prompt emission may be common. We encourage very early follow-up observations of future nearby GRBs, and compact binary merger events, to reveal more about the central engine of GRBs and r-process events.

In the X-ray spectra of AGNs, a noticeable excess of soft X-rays is typically detected beyond the extrapolation of the power-law trend observed between 2-10 keV. In the scenario of warm Comptonization, observations propose a warm corona temperature ranging from 0.1-1 keV and an optical depth of approximately 10-20. Furthermore, according to radiative constraints derived from spectral analyses employing Comptonization models, it is suggested that the majority of the accretion power is released within the warm corona, while the disk beneath it is largely non-dissipative, emitting mainly the reprocessed radiation from the corona. We test the dissipative warm corona model using the radiative transfer code-TITAN/NOAR on a sample of 82 XMM-Newton observations of AGNs. Through spectral modeling of the X-ray data, we aim to estimate the total amount of internal heating inside the warm corona situated on top of the accretion disk. By modeling the 0.3-10 keV EPIC-pn spectra, we estimate the internal heating and optical depth of the warm corona and check their correlations with global parameters blackhole parameters. From model normalization, we compute the radial extent of warm corona on top of cold accretion disk. Our model infers the presence of dissipative warm corona, with optical depths distributed in the range 6-30 and total internal heating in the range 1-29 x 1e-23 erg/s-cm3. The extent of warm corona is spread across a large range from 7-408 gravitational radii, and we find that warm corona is more extended for larger accretion rates. Soft excess emission is ubiquitous in AGNs across wide mass range and accretion rate. We confirm that warm corona responsible for producing the soft-excess is highly dissipative in nature with larger optical depths being associated with lower internal heating and vice versa. The presence of cold standard accretion disk regulates the extent of warm corona.

Many high-state cataclysmic variables (CVs) exhibit blue-shifted absorption features in their ultraviolet (UV) spectra -- a smoking-gun signature of outflows. However, the impact of these outflows on {\em optical} spectra remains much more uncertain. During its recent outburst, the eclipsing dwarf nova V455 And displayed strong optical emission lines whose cores were narrower than expected from a Keplerian disc. Here, we explore whether disc + wind models developed for matching UV observations of CVs can also account for these optical spectra. Importantly, V455~And was extremely bright at outburst maximum: the accretion rate implied by fitting the optical continuum with a standard disc model is $\dot{M}_{\rm acc} \simeq 10^{-7}~{\rm M}_\odot~{\rm yr^{-1}}$. Allowing for continuum reprocessing in the outflow helps to relax this constraint. A disk wind can also broadly reproduce the optical emission lines, but only if the wind is (i) highly mass-loaded, with a mass-loss rate reaching $\dot{M}_{\rm wind} \simeq 0.4 \dot{M}_{\rm acc}$, and/or (ii) clumpy, with a volume filling factor $f_V \simeq 0.1$. The same models can describe the spectral evolution across the outburst, simply by lowering $\dot{M}_{\rm acc}$ and $\dot{M}_{\rm wind}$. Extending these models to lower inclinations and into the UV produces spectra consistent with those observed in face-on high-state CVs. We also find, for the first time in simulations of this type, P-Cygni-like absorption features in the Balmer series, as have been observed in both CVs and X-ray binaries. Overall, dense disc winds provide a promising framework for explaining multiple observational signatures seen in high-state CVs, but theoretical challenges persist.

Products of stellar mergers are predicted to be common in stellar populations and can potentially explain stars with peculiar properties. When the merger occurs after the initially more massive star has evolved into the Hertzsprung gap (HG), the merger product may remain in the blue part of the Hertzsprung-Russell diagram (HRD) for millions of years. Such objects could, therefore, explain the overabundance of observed blue stars, such as blue supergiants. However, it is currently not straightforward to distinguish merger products from genuine single stars. We make detailed asteroseismic comparisons between models of massive post-main-sequence merger products and genuine single stars to identify which asteroseismic diagnostics can be used to distinguish them. In doing so, we develop tools for the relatively young field of merger seismology. Genuine single stars in the HG are fully radiative, while merger products have a convective He-burning core and convective H-burning shell while occupying similar locations in the HRD. These structural differences are reflected in lower asymptotic period spacing values for merger products and the appearance of deep dips in their period spacing patterns. Our genuine single-star models with masses above roughly 11.4 solar masses develop short-lived intermediate convective zones during their HG evolution. This also leads to deep dips in their period spacing patterns. Because of the lack of a convective core, merger products and genuine single stars can be distinguished based on their asymptotic period spacing value in this mass range. We perform the comparisons with and without the effects of slow rotation included in the pulsation equations and conclude that the two types of stars are seismically distinguishable in both cases. The observability of the distinguishing asteroseismic features of merger products can now be assessed and exploited in practice.

High-contrast imaging (HCI) is a technique designed to observe faint signals near bright sources, such as exoplanets and circumstellar disks. The primary challenge in revealing the faint circumstellar signal near a star is the presence of quasi-static speckles, which can produce patterns on the science images that are as bright, or even brighter, than the signal of interest. Strategies such as angular differential imaging (ADI) or reference-star differential imaging (RDI) aim to provide a means of removing the quasi-static speckles in post-processing. In this paper, we present and discuss the adaptation of state-of-the-art algorithms, initially designed for ADI, to jointly leverage angular and reference-star differential imaging (ARDI) for direct high-contrast imaging of circumstellar disks. Using a collection of high-contrast imaging data sets, we assess the performance of ARDI in comparison to ADI and RDI based on iterative principal component analysis (IPCA). These diverse data sets are acquired under various observing conditions and include the injection of synthetic disk models at various contrast levels. Our results demonstrate that ARDI with IPCA improves the quality of recovered disk images and the sensitivity to planets embedded in disks, compared to ADI or RDI individually. This enhancement is particularly pronounced when dealing with extended sources exhibiting highly ambiguous structures that cannot be accurately retrieved using ADI alone, and when the quality of the reference frames is suboptimal, leading to an underperformance of RDI. We finally apply our method to a sample of real observations of protoplanetary disks taken in star-hopping mode, and propose to revisit the protoplanetary claims associated with these disks.

We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using NICER data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gravitational mass are $12.49_{-0.88}^{+1.28}$ km and $2.073_{-0.069}^{+0.069}$ $M_\odot$ respectively, each reported as the posterior credible interval bounded by the $16\,\%$ and $84\,\%$ quantiles, with an estimated systematic error $\lesssim 0.1$ km. This result was obtained using the best computationally feasible sampler settings providing a strong radius lower limit but a slightly more uncertain radius upper limit. The inferred radius interval is also close to the $R=12.76_{-1.02}^{+1.49}$ km obtained by Dittmann et al. 2024, when they require the radius to be less than $16$ km as we do. The results continue to disfavor very soft equations of state for dense matter, with $R<11.15$ km for this high mass pulsar excluded at the $95\,\%$ probability. The results do not depend significantly on the assumed cross-calibration uncertainty between NICER and XMM-Newton. Using simulated data that resemble the actual observations, we also show that our pipeline is capable of recovering parameters for the inferred models reported in this paper.

PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be $2.08\pm0.07\,\rm M_\odot$. Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. (2021) and Riley et al. (2021) reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 17 April 2020, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 28 December 2021 was presented in Salmi et al. (2022). Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 21 April 2022, totaling $\sim1.1$ Ms additional exposure compared to the data set analyzed in Miller et al. (2021) and Riley et al. (2021), and to data from X-ray Multi-Mirror (XMM-Newton) observations. We find that the equatorial circumferential radius of PSR J0740+6620 is $12.92_{-1.13}^{+2.09}$ km (68% credibility), a fractional uncertainty $\sim83\%$ the width of that reported in Miller et al. (2021), in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al. (2024), then our 68% credible region would become $R=12.76^{+1.49}_{-1.02}$ km, which is close to the headline result of Salmi et al. (2024). Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements.

One of the most promising probes to complement current standard cosmological surveys is the HI intensity map, i.e. the distribution of temperature fluctuations in neutral hydrogen. In this paper we present calculations of the 2-point function between HI (at redshift $z$ < 1) and lensing convergence ($\kappa$). We also construct HI intensity maps from N-body simulations, and measure 2-point functions between HI and lensing convergence. HI intensity mapping requires stringent removal of bright foregrounds, including emission from our galaxy. The removal of large-scale radial modes during this HI foreground removal will reduce the HI-lensing cross-power spectrum signal, as radial modes are integrated to find the convergence; here we wish to characterise this reduction in signal. We find that after a simple model of foreground removal, the cross-correlation signal is reduced by $\sim$50-70\%; we present the angular and redshift dependence of the effect, which is a weak function of these variables. We then calculate S/N of $\kappa$HI detection, including cases with cut sky observations, and noise from radio and lensing measurements. We present Fisher forecasts based on the resulting 2-point functions; these forecasts show that by measuring $\kappa\Delta$$T_\mathrm{HI}$ correlation functions in a sufficient number of redshift bins, constraints on cosmology and HI bias will be possible

We conduct X-ray reverberation mapping and spectral analysis of the radio galaxy Centaurus A to uncover its central structure. We compare the light curve of the hard X-ray continuum from Swift Burst Alert Telescope observations with that of the Fe K$\alpha$ fluorescence line, derived from the Nuclear Spectroscopic Telescope Array (NuSTAR), Suzaku, XMM-Newton, and Swift X-ray Telescope observations. The analysis of the light curves suggests that a top-hat transfer function, commonly employed in reverberation mapping studies, is improbable. Instead, the relation between these light curves can be described by a transfer function featuring two components: one with a lag of $0.19_{- 0.02}^{+ 0.10}~\mathrm{pc}/c$, and another originating at $r > 1.7~\mathrm{pc}$ that produces an almost constant light curve. Further, we analyze the four-epoch NuSTAR and six-epoch Suzaku spectra, considering the time lag of the reflection component relative to the primary continuum. This spectral analysis supports that the reflecting material is Compton-thin, with $N_{\mathrm{H}} = 3.14_{-0.74}^{+0.44} \times 10^{23}~ \mathrm{cm}^{-2}$. These results suggest that the Fe K$\alpha$ emission may originate from Compton-thin circumnuclear material located at sub-parsec scale, likely a dust torus, and materials at a greater distance.

The particle-in-cell approach has proven effective at modeling neutron star and black hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where microphysics operates and the system-size scales due to limited numerical resources. A legitimate concern is whether the scale separation currently achieved is large enough, such that results can be safely extrapolated to realistic scales. In this work, our aim is to explore the effect of scaling physical parameters up, and to check whether salient features uncovered by pure kinetic models at smaller scales are still valid, with a special emphasis on particle acceleration and high-energy radiation emitted beyond the light cylinder. To reach this objective, we develop a new hybrid numerical scheme coupling the ideal force-free and the particle-in-cell methods, to optimize the numerical cost of global models. We propose a domain decomposition of the magnetosphere based on the magnetic field topology using the flux function. The force-free model is enforced along open field lines while the particle-in-cell model is restricted to the reconnecting field line region. As a proof of concept, this new hybrid model is applied to simulate a weak millisecond pulsar magnetosphere with realistic scales using high-resolution axisymmetric simulations. Magnetospheric features reported by previous kinetic models are recovered, and strong synchrotron radiation above 100MeV consistent with the Fermi-LAT gamma-ray pulsar population is successfully reproduced. This work further consolidates the shining reconnecting current sheet scenario as the origin of the gamma-ray emission in pulsars, as well as firmly establishes pulsar magnetospheres as at least TeV particle accelerators.

We present ForSE+, a Python package that produces non-Gaussian diffuse Galactic thermal dust emission maps at arcminute angular scales and that has the capacity to generate random realizations of small scales. This represents an extension of the ForSE (Foreground Scale Extender) package, which was recently proposed to simulate non-Gaussian small scales of thermal dust emission using generative adversarial networks (GANs). With the input of the large-scale polarization maps from observations, ForSE+ has been trained to produce realistic polarized small scales at 3' following the statistical properties, mainly the non-Gaussianity, of observed intensity small scales, which are evaluated through Minkowski functionals. Furthermore, by adding different realizations of random components to the large-scale foregrounds, we show that ForSE+ is able to generate small scales in a stochastic way. In both cases, the output small scales have a similar level of non-Gaussianity compared with real observations and correct amplitude scaling as a power law. These realistic new maps will be useful, in the future, to understand the impact of non-Gaussian foregrounds on the measurements of the cosmic microwave background (CMB) signal, particularly on the lensing reconstruction, de-lensing, and the detection of cosmological gravitational waves in CMB polarization B-modes.

Roman microlensing stands at a crossroads between its originally charted path of cataloging a population of cool planets that has subsequently become well-measured down to super-Earths, and the path of free-floating planets (FFPs), which did not exist when Roman was chosen in 2010, but by now promises revolutionary insights into planet formation and evolution via their possible connection to a spectrum of objects spanning 18 decades in mass. Until now, it was not even realized that the 2 paths are in conflict: Roman strategy was optimized for bound-planet detections, and FFPs were considered only in the context of what could be learned about them given this strategy. We derive a simple equation that mathematically expresses this conflict and explains why the current approach severely depresses detection of 2 of the 5 decades of potential FFP masses, i.e., exactly the two decades, $M_{\rm Pluto}< M <2\,M_{\rm Mars}$, that would tie terrestrial planets to the proto-planetary material out of which they formed. FFPs can be either truly free floating or can be bound in "Wide", "Kuiper", and "Oort" orbits, whose separate identification will allow further insight into planet formation. In the (low-mass) limit that the source radius is much bigger than the Einstein radius, $\theta_*\gg\theta_{\rm E}$, the number of significantly magnified points on the FFP light curve is $N=2\Gamma\theta_*\sqrt{1-z^2}/\mu$ --> 3.0, when normalized to the adopted Roman cadence $\Gamma=4/$hr, and to source radius $\theta_*=0.3\,\mu$as, lens-source proper motion $\mu=6\,$mas/yr, and source impact parameter $z=0.5$, which are all typical values. By contrast $N=6$ are needed for an FFP detection. Thus, unless $\Gamma$ is doubled, FFP detection will be driven into the (large-$\theta_*$, small-$\mu$) corner of parameter space, reducing the detections by a net factor of 2 and cutting off the lowest-mass FFPs.

We present $\texttt{cunusht}$, a general-purpose Python package that wraps a highly efficient CUDA implementation of the nonuniform spin-$0$ spherical harmonic transform. The method is applicable to arbitrary pixelization schemes, including schemes constructed from equally-spaced iso-latitude rings as well as completely nonuniform ones. The algorithm has an asymptotic scaling of $\mathrm{O}{(\ell_{\rm max}^3)}$ for maximum multipole $\ell_{\rm max}$ and achieves machine precision accuracy. While $\texttt{cunusht}$ is developed for applications in cosmology in mind, it is applicable to various other interpolation problems on the sphere. We outperform the fastest available CPU algorithm by a factor of up to 5 for problems with a nonuniform pixelization and $\ell_{\rm max}>4\cdot10^3$ when comparing a single modern GPU to a modern 32-core CPU. This performance is achieved by utilizing the double Fourier sphere method in combination with the nonuniform fast Fourier transform and by avoiding transfers between the host and device. For scenarios without GPU availability, $\texttt{cunusht}$ wraps existing CPU libraries. $\texttt{cunusht}$ is publicly available and includes tests, documentation, and demonstrations.

The recent DESI 2024 Baryon Acoustic Oscillations (BAO) measurements combined with the CMB data from the Planck 18 PR3 dataset and the Planck PR4+ACT DR6 lensing data, with a prior on the sum of the neutrino masses $\sum m_\nu>0$, leads to a strong constraint, $\sum m_\nu<0.072$ eV, which would exclude the inverted neutrino hierarchy and put some tension on even the standard hierarchy. We show that actually this bound gets significantly relaxed when combining the new DESI measurements with the HiLLiPoP+LoLLiPoP likelihoods, based on the Planck 2020 PR4 dataset, and with supernovae datasets. We note that the fact that neutrino masses are pushed towards zero, and even towards negative values, is known to be correlated with the so-called $A_L$ tension, a mismatch between lensing and power spectrum measurements in the Planck PR3 data, which is reduced by HiLLiPoP+LoLLiPoP to less than 1$\sigma$. We find $\sum m_\nu<0.1$ eV and $\sum m_\nu<0.12$ eV, with the supernovae Pantheon+ and DES-SN5YR datasets respectively. The shift caused by these datasets is more compatible with the expectations from neutrino oscillation experiments, and both the normal and inverted hierarchy scenarios remain now viable, even with the $\sum m_\nu>0$ prior. Finally, we analyze neutrino mass bounds in an extension of $\Lambda$CDM that addresses the $H_0$ tension, with extra fluid Dark Radiation, finding that in such models bounds are further relaxed and the posterior probability for $\sum m_\nu$ begins to exhibit a peak at positive values.

JWST gives a unique access to the physical and chemical structure of inner disks ($<10$~au), where the majority of the planets are forming. However, the interpretation of mid-infrared (mid-IR) spectra requires detailed thermo-chemical models able to provide synthetic spectra readily comparable to spectroscopic observations. Our goal is to explore the potential of mid-IR emission of OH to probe H$_2$O photodissociation. We include in the DALI disk model prompt emission of OH following photodissociation of H$_2$O in its $\tilde{B}$ electronic state ($\lambda < 144$~nm). This model allows to compute in a self-consistent manner the thermo-chemical structure of the disk and the resulting mid-IR line intensities of OH and H$_2$O. The OH line intensities in the $9-13~\mu$m range are proportional to the total amount of water photodissociated. As such, these lines are a tracer of the amount of water exposed to the FUV field, which depends on the temperature, density, and strength of the FUV field reaching the upper molecular layers. In particular, the OH line fluxes primarily scale with the FUV field emitted by the star in contrast with H$_2$O lines in the 10-20$~\mu$m range which scale with the bolometric luminosity. OH is therefore a key diagnostic to probe the effect of Ly$\alpha$ and constrain the dust FUV opacity in the upper molecular layers. A strong asymmetry between the A' and A'' components of each rotational quadruplet is also predicted. OH mid-IR emission is a powerful tool to probe H$_2$O photodissociation and infer the physical conditions in disk atmospheres. As such, the inclusion of OH mid-IR lines in the analysis of JWST-MIRI spectra will be key for robustly inferring the composition of planet-forming disks. The interpretation of less excited OH lines requires additional quantum calculations of the formation pumping of OH levels by O+H$_2$ and the collisional rate coefficients.