Papers for Thursday, May 18 2023

In the classic Hubble spiral sequence, arm windiness correlates with bulge size; Sa type spiral galaxies with larger bulges also have the most tightly wound spirals. Exceptions to this have long been known, and in recent work using Galaxy Zoo morphologies no strong correlation was seen in a volume limited sample. In this Research Note, we explore the impact of galaxy mass and integrated color upon this correlation in the Galaxy Zoo sample, finding that bluer and lower mass spirals show the ``expected" correlation; however, it becomes slightly negative for redder and/or more massive spiral galaxies.

Both direct and indirect methods of exoplanet detection rely upon detailed knowledge of the potential host stars. Such stellar characterization allows for accurate extraction of planetary properties, as well as contributing to our overall understanding of exoplanetary system architecture. In this analysis, we examine the photometry of 264 known exoplanet host stars (harboring 337 planetary companions) that were observed during the TESS Prime Mission. We identify periodic signatures in the light curves of these stars and make possible connections to stellar pulsations and their rotation periods, and compare the stellar variability to the published planetary orbital periods. From these comparisons, we quantify the effects of stellar variability on exoplanet detection, confirming that exoplanets detection is biased toward lower variability stars, but larger exoplanets dominate the population of exoplanets around variable stars. Exoplanet detection methods represented among these systems are distinct between stellar spectral types across the main sequence, though notable outliers exist. In addition, biases present in both the sourced data from TESS and the host star selection process, which strongly influences the representation of both stellar and planetary characteristics in the final populations. We also determine whether the host star's photometric variability affects or mimics the behavior or properties of the system's planets. These results are discussed in the context of how the behavior of the host star is responsible for how we observe exoplanet characteristics, most notably their radii and atmospheric properties, and how the activity may alter our measurements or impact the evolution of planetary properties.

In high dispersion spectra of rotating bodies such as stars and planets, the rotation contributes significantly to, and sometimes dominates, the line broadening. We present a simple method for rotationally broadening large wavelength ranges of high-dispersion spectra. The broadening is rapid and scales linearly with the length of the spectrum array. For large wavelength ranges, the method is much faster than the popular convolution-based broadening. We provide the code implementation of this method in a publicly accessible repository.

It is believed that the spectral state transitions of the outbursts in X-ray binaries (XRBs) are triggered by the rise of the mass accretion rate due to underlying disc instabilities. Recent observations found that characteristics of disc winds are probably connected with the different spectral states, but the theoretical underpinnings of it are highly ambiguous. To understand the correlation between disc winds and the dynamics of the accretion flow, we have performed General Relativistic Magneto-hydrodynamic (GRMHD) simulations of an axisymmetric thin accretion disc with different accretion rates and magnetic field strengths. Our simulations have shown that the dynamics and the temperature properties depend on both accretion rates and magnetic field strengths. We later found that these properties greatly influence spectral properties. We calculated the average coronal temperature for different simulation models, which is correlated with high-energy Compton emission. Our simulation models reveal that the average coronal temperature is anti-correlated with the accretion rates, which is correlated with the magnetic field strengths. We also found that the structured component of the disc winds (Blandford-Payne disc wind) predominates as the accretion rates and magnetic field strengths increase. In contrast, the turbulent component of the disc winds ($B_{\rm tor}$ disc wind) predominates as the accretion rates and magnetic field strengths decrease. Our results suggest that the disc winds during an outburst in XRBs can only be understood if the magnetic field contribution varies over time (e.g., MAXI J1820+070).

The data produced by the future space-based millihertz gravitational-wave detector LISA will require nontrivial pre-processing, which might affect the science results. It is crucial to demonstrate the feasibility of such processing algorithms and assess their performance and impact on the science. We are building an end-to-end pipeline that includes state-of-the-art simulations and noise reduction algorithms. The simulations must include a detailed model of the full measurement chain, capturing the main features that affect the instrument performance and processing algorithms. In particular, we include in these simulations, for the first time, proper relativistic treatment of reference frames with realistic numerically-optimized orbits; a model for onboard clocks and clock synchronization measurements; proper modeling of total laser frequencies, including laser locking, frequency planning and Doppler shifts; and a better treatment of onboard processing. Using these simulated data, we show that our pipeline is able to reduce the most critical noises and form synchronized observables. By injecting signals from a verification binary, we demonstrate that good parameter estimation can be obtained on this more realistic setup, extending existing results from previous LISA Data Challenges.

We present NOEMA CO (2-1) line and ALMA 870 $\mu$m continuum observations of a main-sequence galaxy at $z=1.38$. The galaxy was initially selected as a "gas-giant", based on the gas mass derived from sub-mm continuum (log$(M_{\rm gas}/M_{\odot})=11.20\pm0.20$), however the gas mass derived from CO (2-1) luminosity brings down the gas mass to a value consistent with typical star-forming galaxies at that redshift (log$(M_{\rm gas}/M_{\odot})=10.84\pm0.03$). Despite that the dust-to-stellar mass ratio remains elevated above the scaling relations by a factor of 5. We explore the potential physical picture and consider an underestimated stellar mass and optically thick dust as possible causes. Based on the updated gas-to-stellar mass ratio we rule out the former, and while the latter can contribute to the dust mass overestimate it is still not sufficient to explain the observed physical picture. Instead, possible explanations include enhanced HI reservoirs, CO-dark H$_2$ gas, an unusually high metallicity, or the presence of an optically dark, dusty contaminant. Using the ALMA data at 870 $\mu$m coupled with $HST$/ACS imaging, we find extended morphology in dust continuum and clumpy star-formation in rest-frame UV in this galaxy, and a tentative $\sim 10$ kpc dusty arm is found bridging the galaxy center and a clump in F814W image. The galaxy shows levels of dust obscuration similar to the so-called $HST$-dark galaxies at higher redshifts, and would fall into the optically faint/dark $JWST$ color-color selection at $z>2$. It is therefore possible that our object could serve as low-$z$ analog of the $HST$-dark populations. This galaxy serves as a caveat to the gas masses based on the continuum alone, with a larger sample required to unveil the full picture.

Upcoming large-scale structure surveys will measure the matter power spectrum to approximately percent level accuracy with the aim of searching for evidence for new physics beyond the standard model of cosmology. In order to avoid biasing our conclusions, the theoretical predictions need to be at least as accurate as the measurements for a given choice of cosmological parameters. However, recent theoretical work has shown that complex physical processes associated with galaxy formation (particularly energetic feedback processes associated with stars and especially supermassive black holes) can alter the predictions by many times larger than the required accuracy. Here we present $\texttt{SP(k)}$, a model for the effects of baryon physics on the non-linear matter power spectrum based on a new large suite of hydrodynamical simulations. Specifically, the ANTILLES suite consists of 400 simulations spanning a very wide range of the "feedback landscape" and show that the effects of baryons on the matter power spectrum can be understood at approaching the percent level in terms of the mean baryon fraction of haloes, at scales of up to $k \lesssim 10 \, h \, $Mpc$^{-1}$ and redshifts up to $z=3$. For the range of scales and redshifts that will be probed by forthcoming cosmic shear measurements, most of the effects are driven by galaxy group-mass haloes ($M \sim 10^{13-14}$ M$_\odot$). We present a simple Python implementation of our model, available at $\href{https://github.com/jemme07/pyspk}{\mathrm{https{:}//github.com/jemme07/pyspk}}$, which can be used to incorporate baryon effects in standard gravity-only predictions, allowing for marginalisation over baryon physics within cosmological pipelines.

We present an empirical analysis of the properties of dust-continuum emission in a sample of 17 galaxies in the early Universe ($4 < z < 8$) with well-sampled far-infrared (FIR) spectral energy distributions (SEDs) compiled from the literature. We place our results into context by self-consistently comparing to samples of nearby star-forming galaxies, luminous infrared galaxies (LIRGs), and quasars. With the exception of two sources, we find no significant evolution in the dust emissivity index across cosmic time, measuring a consistent value of $\beta_\text{IR} = 1.8 \pm 0.3$ at $z > 4$, suggesting the effective dust properties do not change dramatically for most galaxies. Despite having comparable stellar masses, we find the high-redshift galaxies to be similar to, or even more extreme than, LIRGs in the HERUS sample in terms of dust temperature ($T_\text{dust} > 40 \, \mathrm{K}$) and IR luminosity ($L_\text{IR} > 10^{11} \, \mathrm{L_\odot}$). We find the dust temperature evolves mildly towards high redshift, though the LIRGs and quasars exhibit elevated temperatures indicating a more efficient and/or additional heating mechanism. Where available, we compare stellar-mass estimates to our inferred dust masses, whose degeneracy with dust temperature can only be mitigated with a well-constrained SED. In merely half of the cases the dust yield may be explained by supernovae alone, with four sources ($44\%$) significantly exceeding a highly optimistic yield where $M_\text{dust} \approx 0.01 M_*$. We discuss possible explanations for this apparent inconsistency and potential observational biases in the measurements of the dust properties of high-redshift galaxies, including in the current IR-bright sample.

We examine the properties of $\sim50,000$ rotational variables from the ASAS-SN survey using distances, stellar properties, and probes of binarity from $\textit{Gaia}$ DR3 and the SDSS APOGEE survey. They have high amplitudes and span a broader period range than previously studied $\textit{Kepler}$ rotators. We find they divide into three groups of main sequence stars (MS1, MS2s, MS2b) and four of giants (G1/3, G2, G4s, and G4b). MS1 stars are slowly rotating (10-30 days), likely single stars with a limited range of temperatures. MS2s stars are more rapidly rotating (days) single stars spanning the lower main sequence up to the Kraft break. There is a clear period gap (or minimum) between MS1 and MS2s, similar to that seen for lower temperatures in the $\textit{Kepler}$ samples. MS2b stars are tidally locked binaries with periods of days. G1/3 stars are heavily spotted, tidally locked RS CVn with periods of tens of days. G2 stars are less luminous, heavily spotted, tidally locked sub-subgiants with periods of $\sim10$ days. G4s stars have intermediate luminosities to G1/3 and G2, slow rotation periods (approaching 100 days) and are almost certainly all merger remnants. G4b stars have similar rotation periods and luminosities to G4s, but consist of sub-synchronously rotating binaries. We see no difference in indicators for the presence of very wide binary companions between any of these groups and control samples of photometric twin stars built for each group.

X-ray quasi-periodic eruptions (QPEs) are intense repeating soft X-ray bursts from the nuclei of nearby galaxies. Their physical origin is still largely unconstrained, and several theoretical models have been proposed to date. We present here results from a recent XMM-Newton observation of GSN 069, the galactic nucleus where QPEs were first discovered. After about two years of absence, QPEs have reappeared in GSN 069, and we detect two consecutive QPEs separated by a much shorter recurrence time than ever before. Moreover, their intensities and peak temperatures are remarkably different, a novel addition to the QPE phenomenology. We study the QPE spectral properties from all XMM-Newton observations assuming QPEs to either represent an additional emission component superimposed on that from the disc, or the transient evolution of the disc emission itself. In the former scenario, QPEs are consistent with black-body emission from a region that expands by a factor of 2-3 during the individual QPE evolution with radius of the order of that of the Sun at QPE peak. In the alternative non-additive scenario, QPEs originate from a region with an area 6-30 times smaller than the quiescent state X-ray emission, with the smallest regions corresponding to the hottest and most luminous eruptions. The QPE reappearance reveals that eruptions are only present below a quiescent luminosity threshold corresponding to an Eddington ratio of the order of 0.4 for a one million solar masses black hole. The disappearance of QPEs above threshold is most likely driven by the ratio of QPE to quiescence temperatures approaching unity at high quiescent luminosity, making QPE detection challenging, if not impossible, above threshold. We briefly discuss some of the consequences of our results on the proposed models for the QPE physical origin. [abridged]

The first direct measurements of the HI intensity mapping power spectrum were recently made using the MeerKAT telescope. These measurements are on nonlinear scales, at redshifts 0.32 and 0.44. We develop a formalism for modelling small-scale power in redshift space, within the context of the mass-weighted HI halo model framework. This model is consistent with the latest findings from surveys on the HI-halo mass relation. In order to model nonlinear scales, we include the 1-halo, shot-noise and finger-of-god effects. Then we apply the model to the MeerKAT auto-correlation data, finding that the model provides a good fit to the data at redshift 0.32, but the data may indicate some evidence for an adjustment at $z \sim 0.44$. Such an adjustment can be achieved by an increase in the HI halo model bias.

Observations of the 21cm line of neutral hydrogen are poised to revolutionize our knowledge of cosmic reionization and the high-redshift population of galaxies. However, harnessing such information requires robust and comprehensive theoretical modeling. We study the non-linear effects of hydrodynamics and astrophysical feedback processes, including stellar and AGN feedback, on the 21cm signal by post-processing three existing cosmological hydrodynamical simulations of galaxy formation: Illustris, IllustrisTNG, and Eagle. Overall and despite their different underlying galaxy-formation models, the three simulations return similar predictions for the global 21cm rightness temperature and its power spectrum. At fixed redshift, most differences are attributable to differences in the history of reionization, in turn driven by differences in the build-up of stellar sources of radiation. However, the impact of astrophysics is imprinted in the 21cm power spectrum through several unique signatures. First, we find significant small scale ($k \geq 10\, \rm {Mpc}^{-1}$) differences between Illustris and IllustrisTNG, where higher velocity winds generated by supernova feedback soften density peaks and lead to lower 21cm power in TNG. Second, we find more 21cm power at intermediate scales ($k \approx 0.8\, \rm {Mpc}^{-1}$) in Eagle, due to differences in ionization driven by highly effective stellar feedback, leading to lower star formation, older and redder stellar populations, and thus lower ionizing luminosities. Though subtle, these features could allow future observations of the 21cm signal, in conjunction with other reionization observables, to constrain theoretical models for galactic feedback at high redshift.

We report the discovery of four blue lurkers with low and extremely low-mass white dwarf (ELM WDs) companions in the Galactic globular cluster NGC 362 using AstroSat Ultra Violet Imaging Telescope (UVIT). We analyzed the multi-wavelength spectral energy distribution (SED) of FUV-bright MS stars using data from the UVIT, UVOT, GAIA EDR3, and 2.2m ESO/MPI telescopes. Two each of low-mass WDs and ELM WDs are found as companions for the four blue lurkers by the fitting of two-component SED models. The effective temperatures, radii, luminosities, and masses of two low-mass WDs are (35000, 23000) K, (0.04, 0.05) Rsun , (1.45, 0.22) Lsun , and (0.2, 0.2) Msun, while the two ELM WDs are (14750, 14750) K, (0.09, 0.10) Rsun, (0.34, 0.40) Lsun, and (0.18, 0.18) Msun. The position of blue lurkers within the cluster shows that they originated via the Case A/B mass-transfer mechanism in a low-density environment. This is the first detection of blue lurkers with low-mass WDs and ELM WDs as companions in a globular cluster. The companion cooling age is less than 4 Myr, which suggests that they were just recently formed. These binary systems might have originated due to the cluster recent core collapse.

Betelgeuse, the nearest red supergiant, dimmed to an unprecedented level in early 2020. The star emerged from this Great Dimming episode with its typical, roughly 400-day pulsation cycle halved, and a new dominant period of around 200 days. The dimming event has been attributed to a surface mass ejection, in which rising material drove shocks through the stellar atmosphere and expelled some material, partially obscuring the star as it formed molecules and dust. In this paper, we use hydrodynamic simulations to reveal the connections between Betelgeuse's vigorously convective envelope, the surface mass ejection, and the pulsation mode switching that ensued. An anomalously hot convective plume, generated rarely but naturally in the star's turbulent envelope, can rise and break free from the surface, powering an upwelling that becomes the surface mass ejection. The rising plume also breaks the phase coherence of the star's pulsation, causing the surface to keep expanding even as the deeper layers contract. This drives a switch from the 400-day fundamental mode of pulsation, in which the whole star expands and contracts synchronously, to the 200-day first overtone, where a radial node separates the interior and exterior of the envelope moving in opposite phase. We predict that the star's convective motions will damp the overtone oscillation and Betelgeuse will return to its previous, 400-day fundamental mode pulsation in the next 5-10 years. With its resolved surface and unprecedentedly detailed characterization, Betelgeuse opens a window to episodic surface mass ejection in the late-stage evolution of massive stars.

We study the redshift evolution of the baryon budget in a large set of galaxy clusters from the {\it Magneticum} suite of SPH cosmological simulations. At high redshifts, we obtain "closed box" systems independently by the mass of the systems on radii greater than $3R_{500,\mathrm c}$, whereas at lower redshifts, only the most massive halos could be considered as `"closed box". The baryon fraction shows a general decrease with the redshift and, for less massive objects, we observe a much more prominent decrease than for massive halos. The gas depletion parameter $Y_{\rm gas}$ shows a steeper and highly scattered radial distribution in the central regions of less massive halos with respect to massive objects at all redshifts, while on larger radii the gas fraction distributions are independent of the masses or the redshifts. The hot component of the gas traces well the total amount of gas at low redshifts. At higher redshifts, the cold component provides a not negligible contribution to the total amount of baryon in our systems. Moreover, the behaviour of the baryonic, entire gas, and hot gas phase depletion parameters as a function of radius, mass, and redshift are described by some functional forms. The evolution of metallicity and stellar mass in halos suggests that the early enrichment process is dominant. We investigate correlations between the time evolution of AGN feedback and the depletion parameters. We demonstrate that the energy injected by the AGN activity shows a particularly strong positive correlation with $Y_{\rm bar}$, $Y_{\rm cold}$,$Y_{\rm star}$ and a negative one with $Y_{\rm hot}$, $Z_{\rm Tot}$. These trends are consistent with previous works, meaning that our results, combined with findings derived from current and future X-rays observations, represent possible proxies to test the AGN feedback models used in different suites of numerical simulations.

Magnetic fluctuations in the solar wind are often observed to maintain constant magnitude of the magnetic field in a manner consistent with spherically-polarized large-amplitude Alfv\'en waves. We investigate the effect of spherical polarization on the magnetic spectral index through a statistical survey of magnetic fluctuations observed by Parker Solar Probe between 20$R_\odot$ and 200$R_\odot$. We find that deviations from spherical polarization, i.e., changes in $|\mathbf{B}|$ (compressive fluctuations) and one-dimensional discontinuities, have a dramatic effect on the scaling behavior of the turbulent fluctuations. We show that shallow $k^{-3/2}$ spectra are only observed for constant magnetic field strength, three-dimensional structures, which we identify as large amplitude Alfv\'en waves. The presence of compressive fluctuations coincides with a steepening of the spectrum up to $k^{-5/3}$. Steeper power law scalings approaching $k^{-2}$ are observed when the fluctuations are dominated by discontinuities. Near-sun fluctuations are found to be the most spherically polarized, suggesting that this spherical state is fundamental to the generation of the solar wind. With increasing distance from the Sun, fluctuations are found to become less three dimensional and more compressive, which may indicate the breakdown of the Alfv\'enic equilibrium state.

For natural guide start adaptive optics (AO) systems, pyramid wavefront sensors (PWFSs) can provide significant increase in sensitivity over the traditional Shack-Hartmann, but at the cost of a reduced linear range. When using a linear reconstructor, non-linearities result in wavefront estimation errors, which can have a significant impact on the image quality delivered by the AO system. Here we simulate a wavefront passing through a PWFS under varying observing conditions to explore the possibility of using a non-linear machine learning model to estimate wavefront errors better than a linear reconstruction. We find significant improvement even with light-weight models, underscoring the need for further investigation of this approach.

If they exist, networks of superconducting cosmic strings are capable of injecting copious amounts of electromagnetic energy into the background over a broad range of frequencies. We study this injection both analytically, as well as numerically using the thermalization code CosmoTherm. With our refined analytic formalism, we update constraints from CMB spectral distortions by following the injection of entropy, as well as energy, on the amplitude of the $\mu$-distortion, leading to a significant improvement in those limits. Furthermore, we utilize the full shape of the distorted spectrum from CosmoTherm to include constraints from non-$\mu$, non-$y$ type distortions. Additionally, we use the outputs for the ionization history and global 21cm signal to derive and update constraints on string model parameters using measurements from other datasets. Analysis of CMB anisotropies provides the most stringent constraints, though with a slightly modified shape and strength when compared to previous results. Modifications of the reionization history provide new bounds in the high current domain, and we also find that the observations of the low-frequency radio background probe a small region of parameter space not explored by other datasets. We also analyze global $21$-cm constraints, and find that the inclusion of soft photon heating plays a crucial role, essentially removing any constraints in the considered parameter domain. Spectral distortion measurements from COBE/FIRAS are covered by other constraints, but our conservative forecast shows that a PIXIE-type satellite would probe important unexplored regions of parameter space.

Cosmological simulations predict dark matter shapes that deviate from spherical symmetry. The exact shape depends on the prescription of the simulation and the interplay between dark matter and baryons. This signature is most pronounced in the diffuse galactic haloes that can be observationally probed with planetary nebulae and globular clusters (GCs). The kinematic observations of these halo tracers support intrinsic triaxial shape for the mass generating the gravitational potential. With discrete axisymmetric modelling of GCs as the halo tracers of NGC 5128 we investigate the overall mass distribution of this nearby giant elliptical galaxy. Our modelling approach constrains $c_{200}$, $(M/L)_{\star, B}$ and inclination. We derive a preliminary $M_{200}\sim 1 \times 10^{12}$ M$_\odot$ and flattening $q_{\mathrm{DM}}\sim 1.3$ indicative of prolate/triaxial halo for NGC 5128.

Weak galaxy lensing surveys have consistently reported low values of the $S_8$ parameter compared to the $\textit{Planck}\ \Lambda\rm{CDM}$ cosmology. Amon & Efstathiou (2022) used KiDS-1000 cosmic shear measurements to propose that this tension can be reconciled if the matter fluctuation spectrum is suppressed more strongly on non-linear scales than assumed in state-of-the-art hydrodynamical simulations. In this paper, we investigate cosmic shear data from the Dark Energy Survey (DES) Year 3. The non-linear suppression of the matter power spectrum required to resolve the $S_8$ tension between DES and the $\textit{Planck}\ \Lambda\rm{CDM}$ model is not as strong as inferred using KiDS data, but is still more extreme than predictions from recent numerical simulations. An alternative possibility is that non-standard dark matter contributes to the required suppression. We investigate the redshift and scale dependence of the suppression of the matter power spectrum. If our proposed explanation of the $S_8$ tension is correct, the required suppression must extend into the mildly non-linear regime to wavenumbers $k\sim 0.2 h {\rm Mpc}^{-1}$. In addition, all measures of $S_8$ using linear scales should agree with the $\textit{Planck}\ \Lambda\rm{CDM}$ cosmology, an expectation that will be testable to high precision in the near future.

We point out the dominant importance of plasma injection effects for relativistic winds from pulsars and black holes. We demonstrate that outside the light cylinder the magnetically dominated outflows while sliding along the helical magnetic field move in fact nearly radially with very large Lorentz factors $\gamma_0 \gg 1 $, imprinted into the flow during pair production within the gaps. Only at larger distances, $r \geq \gamma_0 (c/\Omega)$, the MHD acceleration $\Gamma \propto r$ takes over. As a result, Blandford-Znajek (BZ) driven outflows would produce spine-brightened images, contrary to observations of the edge-brightened jet in M87. We conclude that M87 jet is not BZ-driven. Other implications include: (i) variability time scale of emission even from non-blazar AGNe like M87 and Cen A can be much shorter than the BH spin period; (ii) the model explains blazar-like phenomena in clearly non-aligned AGNe; (iii) in PIC simulations of pulsars' and black holes' magnetospheres with pair production, the pairs should be injected with large Lorentz factors - this would drastically change the resulting magnetospheric and inner-wind structure, if compared with pair injection at rest; (iv) the model resolves the problem of cyclotron absorption in the Crab pulsar wind.

The thermal Sunyaev-Zel'dovich (tSZ) effect is a spectral distortion of the cosmic microwave background (CMB) resulting from inverse Compton scattering of CMB photons with electrons in the medium of galaxy clusters. The spectrum of the tSZ effect is typically calculated assuming the spectrum of the CMB is a blackbody. However, energy or photon number injection at any epoch after photon creation processes become inefficient will distort the blackbody, potentially leading to a chemical potential or $\mu$-distortion for early injection. These $primordial$ spectral distortions will therefore introduce a change in the tSZ effect, effectively a distortion of a distortion. While this effect is small for an individual cluster's spectrum, upcoming and proposed CMB surveys expect to detect tens of thousands of clusters with the tSZ effect. In this paper, we forecast constraints on the $\mu$-distortion monopole from the distortion of the tSZ spectrum of clusters measured by CMB surveys. We find that planned experiments have the raw sensitivity to place constraints on $\mu$ that are comparable to or better than existing constraints but control over foregrounds and other systematics will be critical.

We identify cosmic voids from galaxy density fields under the theory of void-cluster correspondence. We extend the previous novel void-identification method developed for the matter density field to the galaxy density field for practical applications. From cosmological N-body simulations, we construct galaxy number- and mass-weighted density fields to identify cosmic voids that are counterparts of galaxy clusters of specific mass. The parameters for the cluster-counterpart void identification such as Gaussian smoothing scale, density threshold, and core volume fraction are found for galaxy density fields. We achieve about $60$--$67\%$ of completeness and reliability for identifying the voids of corresponding cluster mass above $3\times10^{14}h^{-1}M_{\odot}$ from a galaxy sample with the mean number density, $\bar{n}=4.4\times10^{-3} (h^{-1}{\rm Mpc})^{-3}$. When the mean density is increased to $\bar{n}=10^{-2} (h^{-1}{\rm Mpc})^{-3}$, the detection rate is enhanced by $\sim2$--$7\%$ depending on the `mass scale' of voids. We find that the detectability is insensitive to the density weighting scheme applied to generate the density field. Our result demonstrates that we can apply this method to the galaxy redshift survey data to identify cosmic voids corresponding statistically to the galaxy clusters in a given mass range.

Recently, it was shown that the formation of a photon-trapping surface might not be sufficient to ensure unimpeded super-Eddington (SE) accretion. In light of this finding, here we derive a condition such that sustained and unimpeded SE accretion could be achieved in optically thick slim accretion disks surrounding neutron stars (NSs) and black holes (BHs). For this, we calculate a semi-analytic approximation of the self-similar global radial velocity expression for an advection-dominated flow. Neglecting the influence of relativistic jets on the accretion flow, we find that for Eddington fraction $\dot{m} \gtrsim 1.5 (\epsilon/0.1)^{3/5}$ ($\epsilon$ being the accretion efficiency) sustained SE accretion might be possible in slim disks around BHs irrespective of their spin. The same condition holds for NSs when $\epsilon > 0.03$. The presence of a surface magnetic field might truncate the disk at the magnetosphere of the NS, resulting in lower efficiencies and consequently changing the condition to $\dot{m} > 0.013 \epsilon^{-19/31}$. Our approach suggests that sustained SE accretion might almost always be possible around NSs and BHs hosting accretion disks.

Recent measurements of Jupiter's gravitational moments by the Juno spacecraft and seismology of Saturn's rings suggest that the primordial composition gradients in the deep interior of these planets have persisted since their formation. One possible explanation is the presence of a double-diffusive staircase below the planet's outer convection zone, which inhibits mixing across the deeper layers. However, hydrodynamic simulations have shown that these staircases are not long-lasting and can be disrupted by overshooting convection. In this paper we suggests that planetary rotation could be another factor for the longevity of primordial composition gradients. Using rotational mixing-length theory and 3D hydrodynamic simulations, we demonstrate that rotation significantly reduces both the convective velocity and the mixing of primordial composition gradients. In particular, for Jovian conditions at $t\sim10^{8}~\mathrm{yrs}$ after formation, rotation reduces the convective velocity by a factor of 6, and in turn, the kinetic energy flux available for mixing gets reduced by a factor of $6^3\sim 200$. This leads to an entrainment timescale that is more than two orders of magnitude longer than without rotation. We encourage future hydrodynamic models of Jupiter and other gas giants to include rapid rotation, because the decrease in the mixing efficiency could explain why Jupiter and Saturn are not fully mixed.

Observation of several very high energy (VHE) flaring events of the BL Lac object VER J0521+211 were reported by the VERITAS and MAGIC collaborations between 2009 and 2014. The redshift of this source is uncertain and several analysis have derived different limits for it. In the framework of the photohadronic model, and using three different extragalactic background light (EBL) models, we analyze seven independent VHE spectra of VER J0521+211 and determine the limiting values on its redshift. It is observed that the photohadronic scenario provides excellent fits to the reported observations. It is further observed that the photohadronic scenario, along with the EBL model of Dominguez et al., puts the most restrictive limits on the redshift z of VER J0521+211: 0.29 <= z <= 0.31 from the confidence level (CL) intervals at 2 sigma, or a more conservative 0.28 <= z <= 0.33, at 3 sigma.

We present 850 $\mu$m polarization and $\rm C^{18}O (3-2)$ molecular line observations toward the X-shaped nebula in the California molecular cloud using the JCMT SCUBA-2/POL-2 and HARP instruments. The 850 $\mu$m emission shows that the observed region includes two elongated filamentary structures (Fil1 and Fil2) having chains of regularly spaced cores. We measured the mass per unit length of the filament and found that Fil1 and Fil2 are thermally super- and subcritical, respectively, but both are subcritical if nonthermal turbulence is considered. The mean projected spacings ($\Delta\bar S$) of cores in Fil1 and Fil2 are 0.13 and 0.16 pc, respectively. $\Delta\bar S$ are smaller than $4\times$filament width expected in the classical cylinder fragmentation model. The large-scale magnetic field orientations shown by Planck are perpendicular to the long axes of Fil1 and Fil2, while those in the filaments obtained from the high-resolution polarization data of JCMT are disturbed, but those in Fil1 tend to have longitudinal orientations. Using the modified Davis-Chandrasekhar-Fermi (DCF) method, we estimated the magnetic field strengths ($B_{\rm pos}$) of filaments which are 110$\pm$80 and 90$\pm$60 $\mu$G. We calculated the gravitational, kinematic, and magnetic energies of the filaments, and found that the fraction of magnetic energy is larger than 60 % in both filaments. We propose that a dominant magnetic energy may lead the filament to be fragmented into aligned cores as suggested by Tang et al., and a shorter core spacing can be due to a projection effect via the inclined geometry of filaments or due to a non-negligible, longitudinal magnetic fields in case of Fil1.

We investigated the inner buried nucleus of a nearby luminous infrared galaxy NGC 4418 using high-resolution spectroscopy of fundamental carbon monoxide (CO) ro-vibrational absorptions around $4.67 \mu$m for the first time. This method allowed us to examine the physical and kinematical properties in the hot inner region of this nucleus. We detected a series of both very deep (partly saturated) $^{12}$CO and moderately deep (optically thin) $^{13}$CO absorption lines and inferred a large column density ($N_\mathrm{H2}=(5\pm3)\times10^{23}$ cm$^{-2}$ in front of the $5 \mu$m photosphere) of warm ($T_\mathrm{ex}\simeq170$ K) molecular gas by assuming an isothermal plane-parallel slab illuminated by a compact background MIR-emitting source. We modeled that the warm CO absorber almost covers the central heating source and that it is an inner layer around the $5 \mu$m photosphere (at $r=$several pc) of a compact shroud of gas and dust ($d\sim100$ pc). The width of the absorption lines ($110$ km s$^{-1}$) and their small deviation from the systemic velocity ($<10$ km s$^{-1}$) are consistent with a warm and turbulent layer with little bulk motion in the radial direction.

Using magnetohydrodynamic simulations of Fluctuation dynamos in turbulent flows with rms Mach numbers $\mathcal{M} \approx 0.2, 1.1$ and $3$, we show that magnetic pressure forces play a crucial role in dynamo saturation in supersonic flows. Firstly, as expected when pressure forces oppose compression, an increase in anti-correlation between density and magnetic field strengths obtains even in subsonic flows with the anti-correlation arising from the intense but rarer magnetic structures. In supersonic flows, due to stronger compressive motions density and magnetic field strength continue to maintain a positive correlation. However, the degree of positive correlation decreases as the dynamo saturates. Secondly, we find that the unit vectors of $\nabla\rho$ and $\nabla B^{2}$ are preferentially anti-parallel to each other in subsonic flows. This is indicative of magnetic pressure opposing compression. This anti-parallel alignment persists in transonic and supersonic flows at dynamo saturation. However, compressive motions also lead to the emergence of a parallel alignment in these flows. Finally, we consider the work done against the components of the Lorentz force and the different sources of magnetic energy growth and dissipation. We show that while in subsonic flows, suppression of field line stretching is dominant in saturating the dynamo, the picture is different in supersonic flows. Both field line stretching and compression amplifies the field initially. But the growing magnetic pressure opposes further compression of magnetic flux which then dominates the saturation of the dynamo.

We construct a numerical model of steady-state, general relativistic (GR) super-Eddington accretion flows in an optically thick, advection-dominated regime, motivated by tidal disruption events wherein super-Eddington accretion assumes a pivotal role. Our model takes into account the loss of angular momentum due to radiation and the scale-height derivative in the basic equations of the GR slim disk. For comparison purposes, we also provide a new analytical solution for a radiation-pressure-dominant GR slim disk, which neglects the angular momentum loss due to radiation and the scale-height derivative. We find that the radiation pressure enhances by incorporating the scale height derivative into the basic equations. As a result, the surface density near the disk's inner edge decreases, whereas the disk temperature and scale height increase, brightening the disk spectrum in the soft X-ray waveband. Notably, an extremely high mass accretion rate significantly enhances the effect of the scale-height derivative, affecting the entire disk. In contrast, the inclusion of the radiation-driven angular momentum loss only slightly influences the disk surface density and temperature compared with the case of the scale-height derivative inclusion. The X-ray luminosity increases significantly due to scale height derivative for $\dot{M}/\dot{M}_{\rm Edd} \gtrsim 2$. In addition, the increment is higher for the non-spinning black hole than the spinning black hole case, resulting in a one-order of magnitude difference for $\dot{M}/\dot{M}_{\rm Edd}\gtrsim100$. We conclude that incorporating the scale-height derivative into a GR slim disk model is crucial as it impacts the disk structure and its resultant spectrum, particularly on a soft-X-ray waveband.

We present the intra-cluster kinematics and dynamics of three open clusters: NGC 1193, NGC 2355, and King 12 by incorporating kinematical and photometric data from Gaia DR3, as well as a ground-based telescope. After selecting cluster members based on proper motion data, clusters' fundamental and structural parameters are investigated. We found the clusters at distances of 4.45, 1.97, and 3.34 kpc from the Sun in the direction of the Galactic anticenter. The luminosity function of the cluster NGC 1193 is flat, whereas it advances towards the fainter ends of the other two clusters. We observed a dip in the luminosity function of King 12. The mass function slopes for all three clusters differ from the solar neighbourhood reported by Salpeter, with NGC 1193 and NGC 2355 being flatter and King 12 having a higher value than the Salpeter value. The intra-cluster kinematics depict that stars in King 12 are moving outwards due to tidal forces from the Galactic disc, which we confirmed by plotting the cluster's orbit in the Galaxy. Stars in NGC 2355 are moving with smaller relative velocities and have zero mean relative motion, which signifies that the cluster is neither contracting nor evaporating. The Galactic orbits of NGC 1193 suggest that it is orbiting farther from the Galactic disc, and so is less impacted by the Galactic tidal forces.

Cosmological models and their parameters are widely debated, especially about whether the current discrepancy between the values of the Hubble constant, $H_{0}$, obtained by type Ia supernovae (SNe Ia), and the Planck data from the Cosmic Microwave Background Radiation could be alleviated when alternative cosmological models are considered. Thus, combining high-redshift probes, such as Gamma-Ray Bursts (GRBs) and Quasars (QSOs), together with Baryon Acoustic Oscillations (BAO) and SNe Ia is important to assess the viability of these alternative models and if they can cast further light on the Hubble tension. In this work, for GRBs, we use a 3-dimensional relation between the peak prompt luminosity, the rest-frame time at the end of the X-ray plateau, and its corresponding luminosity in X-rays: the 3D Dainotti fundamental plane relation. Regarding QSOs, we use the Risaliti-Lusso relation among the UV and X-ray luminosities for a sample of 2421 sources. We correct both the QSO and GRB relations by accounting for selection and evolutionary effects with a reliable statistical method. We here use both the traditional Gaussian likelihoods ($\cal L_G$) and the new best-fit likelihoods ($\cal L_N$) to infer cosmological parameters of a non-flat $\Lambda$CDM and flat $w$CDM models. We obtain for all the parameters reduced uncertainties, up to $35\%$ for $H_{0}$, when applying the new $\cal L_N$ likelihoods in place of the Gaussian ones. Our results remain consistent with a flat $\Lambda$CDM model, although with a shift of the dark energy parameter $w$ toward $w<-1$ and a curvature density parameter toward $\Omega_k<0$.

The Ly$\alpha$ line is a powerful probe of distant galaxies, which contains information about inflowing/outflowing gas through which Ly$\alpha$ photons scatter. To develop our understanding of this probe, we post-process a zoom-in radiation-hydrodynamics simulation of a low-mass ($M_* \sim 10^9 M_\odot$) galaxy to construct 22500 mock spectra in 300 directions from $z = 3$ to 4. Remarkably, we show that one galaxy can reproduce the variety of a large sample of spectroscopically observed Ly$\alpha$ line profiles. While most mock spectra exhibit double-peak profiles with a dominant red peak, their shapes cover a large parameter space in terms of peak velocities, peak separation and flux ratio. This diversity originates from radiative transfer effects at ISM and CGM scales, and depends on galaxy inclination and evolutionary phase. Red-dominated lines preferentially arise in face-on directions during post-starburst outflows and are bright. Conversely, accretion phases usually yield symmetric double peaks in the edge-on direction and are fainter. While resonant scattering effects at $< 0.2\times R_{\rm vir}$ are responsible for the broadening and velocity shift of the red peak, the extended CGM acts as a screen and impacts the observed peak separation. The ability of simulations to reproduce observed Ly$\alpha$ profiles and link their properties with galaxy physical parameters offers new perspectives to use Ly$\alpha$ to constrain the mechanisms that regulate galaxy formation and evolution. Notably, our study implies that deeper Ly$\alpha$ surveys may unveil a new population of blue-dominated lines tracing inflowing gas.

Coronal Mass Ejections (CMEs) correspond to dramatic expulsions of plasma and magnetic field from the solar corona into the heliosphere. CMEs are scientifically relevant because they are involved in the physical mechanisms characterizing the active Sun. However, more recently CMEs have attracted attention for their impact on space weather, as they are correlated to geomagnetic storms and may induce the generation of Solar Energetic Particles streams. In this space weather framework, the present paper introduces a physics-driven artificial intelligence (AI) approach to the prediction of CMEs travel time, in which the deterministic drag-based model is exploited to improve the training phase of a cascade of two neural networks fed with both remote sensing and in-situ data. This study shows that the use of physical information in the AI architecture significantly improves both the accuracy and the robustness of the travel time prediction.

This paper investigates a poorly studied open cluster, NGC 5288, using 2MASS JHKS and the recently released Gaia DR3 astrometric and photometric data. The mean proper motions in Right Ascension and Declination are estimated as (-3.840 +/- 0.230) and (-1.934 +/- 0.162) mas/yr, respectively. We also derive the age and distance of the cluster as 510 +/- 190 Myr and 2.64 +/- 0.11 kpc, using colour-magnitude diagrams (CMDs). We have also obtained distance as 2.77 +/- 0.42 kpc using the parallax method. Interstellar reddening E(B-V) in the direction of the cluster is determined as 0.45 mag using the ((J - H), (J - K)) colour-colour diagram. We have found the mass function slope for main-sequence stars as 1.39 +/- 0.29 within the mass range 1.0 - 2.7 solar mass, which agrees with Salpeter's value within uncertainty. Galactic orbits are derived using the Galactic potential model, indicating that NGC 5288 follows a circular path around the Galactic center.

Simulations of the radio emission from extensive air showers have been key in establishing radio detection as a mature and competitive technique. In particular, microscopic Monte Carlo simulations have proven to very accurately describe the emission physics and are at the heart of practically all analysis approaches. Yet with new applications -- for example very inclined air showers, cross-media showers, extreme antenna densities, and higher-frequency measurements -- come new challenges for accurate and efficient simulations. I will review the state of the art of the existing simulation approaches and discuss where further improvements might be needed and how they can be achieved.

Our knowledge of Giga-electron volt (GeV) radio galaxies has been revolutionized by the Fermi-LAT Telescope, which provides an excellent opportunity to study the physical properties of GeV radio galaxies. According to the radio power and morphology, radio galaxies can be separated into Fanaroff-Riley Type I radio galaxies (FR-Is) and Type II radio galaxies (FR-IIs). In this paper, we consider the unification of FR-Is and BL Lacertae objects (BL Lacs), and assume FR-Is to be a standard candle to discuss the beaming effect for Fermi-LAT-detected FR-Is. Our main conclusions are as follows: (1) The estimated Doppler factors for 30 Fermi-LAT-detected FR-Is are in a range of $\delta_{\rm{I}}=0.88-7.49$. The average Doppler factor ($<\delta_{\rm{I}}>=2.56\pm0.30$) of the 30 FR-Is is smaller than that ($<\delta_{\rm{BL}}>=10.28\pm2.03$) of the 126 Fermi-LAT-detected BL Lacs, supporting the unification model that FR-Is are regarded as the misaligned BL Lacs with smaller Doppler factors; (2) We propose that different regions of FR-Is in the plot of the $\gamma$-ray luminosity against the photon spectral index $(\log L_{\gamma}-\alpha_{\rm{ph}})$ may indicate the different beaming effects; (3) The average Doppler factor of the 6 TeV FR-Is is similar to that of the 24 non-TeV FR-Is, which implies that the difference between the TeV and GeV emissions is not driven by the beaming effect in the Fermi-LAT-detected FR-I samples.

As part of the ongoing AugerPrime upgrade of the Pierre Auger Observatory, we are deploying short aperiodic loaded loop antennas measuring radio signals from extensive air showers in the 30-80 MHz band on each of the 1,660 surface detector stations. This new Radio Detector of the Observatory allows us to measure the energy in the electromagnetic cascade of inclined air showers with zenith angles larger than $\sim 65^\circ$. The water-Cherenkov detectors, in turn, perform a virtually pure measurement of the muon component of inclined air showers. The combination of both thus extends the mass-composition sensitivity of the upgraded Observatory to high zenith angles and therefore enlarges the sky coverage of mass-sensitive measurements at the highest energies while at the same time allowing us to cross-check the performance of the established detectors with an additional measurement technique. In this contribution, we outline the concept and design of the Radio Detector, report on its current status and initial results from the first deployed stations, and illustrate its expected performance with a detailed, end-to-end simulation study.

Molecular contamination is a well-known problem in space flight. Water is the most common contaminant and alters numerous properties of a cryogenic optical system. Too much ice means that Euclid's calibration requirements and science goals cannot be met. Euclid must then be thermally decontaminated, a long and risky process. We need to understand how iced optics affect the data and when a decontamination is required. This is essential to build adequate calibration and survey plans, yet a comprehensive analysis in the context of an astrophysical space survey has not been done before. In this paper we look at other spacecraft with well-documented outgassing records, and we review the formation of thin ice films. A mix of amorphous and crystalline ices is expected for Euclid. Their surface topography depends on the competing energetic needs of the substrate-water and the water-water interfaces, and is hard to predict with current theories. We illustrate that with scanning-tunnelling and atomic-force microscope images. Industrial tools exist to estimate contamination, and we must understand their uncertainties. We find considerable knowledge errors on the diffusion and sublimation coefficients, limiting the accuracy of these tools. We developed a water transport model to compute contamination rates in Euclid, and find general agreement with industry estimates. Tests of the Euclid flight hardware in space simulators did not pick up contamination signals; our in-flight calibrations observations will be much more sensitive. We must understand the link between the amount of ice on the optics and its effect on Euclid's data. Little research is available about this link, possibly because other spacecraft can decontaminate easily, quenching the need for a deeper understanding. In our second paper we quantify the various effects of iced optics on spectrophotometric data.

Coded mask instruments have been used in high-energy astronomy for the last forty years now and designs for future hard X-ray/low gamma-ray telescopes are still based on this technique when they need to reach moderate angular resolutions over large field of views, particularly for observations dedicated to the, now flourishing, field of time domain astrophysics. However these systems are somehow unfamiliar to the general astronomers as they actually are two-step imaging devices where the recorded picture is very different from the imaged object and the data processing takes a crucial part in the reconstruction of the sky image. Here we present the concepts of these optical systems applied to high-energy astronomy, the basic reconstruction methods including some useful formulae and the trend of the expected and observed performances as function of the system designs. We review the historical developments and recall the flown space-borne coded mask instruments along with the description of a few relevant examples of major successful implementations and future projects in space astronomy.

In this paper we extend the Tiered Radio Extragalactic Continuum Simulation (T-RECS) to include HI emission. The HI T-RECS model is based on the most recent HI mass function estimates, combined with prescriptions to convert HI mass to total integrated HI flux. It further models source size, morphology and kinematics, including rotational velocity and HI line width. The continuum T-RECS model is updated to improve the agreement with deeper number counts available at 150\,MHz. The model for star-forming galaxies (SFGs) is also modified according to the most recent indications of a star formation rate (SFR)--radio luminosity relation, which depends primarily on stellar mass rather than redshift. We further introduce prescriptions to associate an HI mass to the T-RECS radio continuum SFG and Active Galactic Nuclei (AGN) populations. This gives us a way to meaningfully associate counterparts between HI and continuum catalogues, thus building HI $\times$ continuum simulated observations. Clustering properties of the sources in both HI and continuum are reproduced by associating the galaxies to dark matter haloes of a cosmological simulation. We deliver a set of mock catalogues, as well as the code to produce them, which can be used for simulating observations and predicting results from radio surveys with existing and forthcoming radio facilities, such as the Square Kilometre Array (SKA)

The analysis of absorption features along the line of sight to distant sources is an invaluable tool for observational cosmology, giving a direct insight into the physical and chemical state of the inter/circumgalactic medium. Such endeavour entails the accessibility of bright QSOs as background beacons, and the availability of software tools to extract the information in a reproducible way. In this article, we will present the latest results we obtained in both directions within the QUBRICS project: we will describe how machine learning techniques were applied to detect hundreds of previously unknown QSOs in the southern hemisphere, and how state-of-the art software like QSFit and Astrocook was integrated in the analysis of the targets, opening up new possibilities for the next era of observations.

Taking advantage of the reduced levels of noise and systematics in the data of the latest Planck release (PR4, also known as NPIPE), we construct a new all-sky Compton-$y$ parameter map (hereafter, $y$-map) of the thermal Sunyaev-Zeldovich (SZ) effect from the Planck PR4 data. A tailored Needlet Internal Linear Combination (NILC) pipeline, first validated on detailed sky simulations, is applied to the nine single-frequency Planck PR4 sky maps, ranging from $30$ to $857$ GHz, to produce the PR4 $y$-map over 98% of the sky. Using map comparisons, angular power spectra and one-point statistics we show that the PR4 NILC $y$-map is of improved quality compared to that of the previous PR2 release. The new $y$-map shows reduced levels of large-scale striations associated with $1/f$ noise in the scan direction. Regions near the Galactic plane also show lower residual contamination by Galactic thermal dust emission. At small angular scales, the residual contamination by thermal noise and cosmic infrared background (CIB) emission is found to be reduced by around 7% and 34%, respectively, in the PR4 $y$-map. The PR4 NILC $y$-map is made publicly available for astrophysical and cosmological analyses of the thermal SZ effect.

We report on the second observation of the radio-quiet active galactic nucleus (AGN) MCG-05-23-16 performed with the Imaging X-ray Polarimetry Explorer (IXPE). The observation started on 2022 November 6 for a net observing time of 640 ks, and was partly simultaneous with NuSTAR (86 ks). After combining these data with those obtained in the first IXPE pointing on May 2022 (simultaneous with XMM-Newton and NuSTAR) we find a 2-8 keV polarization degree $\Pi$ = 1.6 $\pm$ 0.7 (at 68 per cent confidence level), which corresponds to an upper limit $\Pi$ = 3.2 per cent (at 99 per cent confidence level). We then compare the polarization results with Monte Carlo simulations obtained with the MONK code, with which different coronal geometries have been explored (spherical lamppost, conical, slab and wedge). Furthermore, the allowed range of inclination angles is found for each geometry. If the best fit inclination value from a spectroscopic analysis is considered, a cone-shaped corona along the disc axis is disfavoured.

The Wide-field Infrared Survey Explorer (WISE) has detected hundreds of millions of sources over the entire sky. However, classifying them reliably is a great challenge due to degeneracies in WISE multicolor space and low detection levels in its two longest-wavelength bandpasses. In this paper, the deep learning classification network, IICnet (Infrared Image Classification network), is designed to classify sources from WISE images to achieve a more accurate classification goal. IICnet shows good ability on the feature extraction of the WISE sources. Experiments demonstrates that the classification results of IICnet are superior to some other methods; it has obtained 96.2% accuracy for galaxies, 97.9% accuracy for quasars, and 96.4% accuracy for stars, and the Area Under Curve (AUC) of the IICnet classifier can reach more than 99%. In addition, the superiority of IICnet in processing infrared images has been demonstrated in the comparisons with VGG16, GoogleNet, ResNet34, MobileNet, EfficientNetV2, and RepVGG-fewer parameters and faster inference. The above proves that IICnet is an effective method to classify infrared sources.

We present simulated galaxy spectral energy distributions (SEDs) from the far ultraviolet through the far infrared, created using hydrodynamic simulations and radiative transfer calculations, suitable for the validation of SED modeling techniques. SED modeling is an essential tool for inferring star formation histories from nearby galaxy observations, but is fraught with difficulty due to our incomplete understanding of stellar populations, chemical enrichment processes, and the non-linear, geometry dependent effects of dust on our observations. Our simulated SEDs will allow us to assess the accuracy of these inferences against galaxies with known ground truth. To create the SEDs, we use simulated galaxies from the NIHAO suite and the radiative transfer code SKIRT. We explore different sub-grid post-processing recipes, using color distributions and their dependence on axis ratio of galaxies in the nearby universe to tune and validate them. We find that sub-grid post-processing recipes that mitigate limitations in the temporal and spatial resolution of the simulations are required for producing FUV to FIR photometry that statistically reproduce the colors of galaxies in the nearby universe. With this paper we release resolved photometry and spatially integrated spectra for our sample galaxies, each from a range of different viewing angles. Our simulations predict that there is a large variation in attenuation laws among galaxies, and that from any particular viewing angle that energy balance between dust attenuation and reemission can be violated by up to a factor of 3. These features are likely to affect SED modeling accuracy.

Mapping molecular line emission beyond the bright low-J CO transitions is still challenging in extragalactic studies, even with the latest generation of (sub-)mm interferometers, such as ALMA and NOEMA. We summarise and test a spectral stacking method that has been used in the literature to recover low-intensity molecular line emission, such as HCN(1-0), HCO+(1-0), and even fainter lines in external galaxies. The goal is to study the capabilities and limitations of the stacking technique when applied to imaged interferometric observations. The core idea of spectral stacking is to align spectra of the low S/N spectral lines to a known velocity field calculated from a higher S/N line expected to share the kinematics of the fainter line, e.g., CO(1-0) or 21-cm emission. Then these aligned spectra can be coherently averaged to produce potentially high S/N spectral stacks. Here, we use imaged simulated interferometric and total power observations at different signal-to-noise levels, based on real CO observations. For the combined interferometric and total power data, we find that the spectral stacking technique is capable of recovering the integrated intensities even at low S/N levels across most of the region where the high S/N prior is detected. However, when stacking interferometer-only data for low S/N emission, the stacks can miss up to 50% of the emission from the fainter line. A key result of this analysis is that the spectral stacking method is able to recover the true mean line intensities in low S/N cubes and to accurately measure the statistical significance of the recovered lines. To facilitate the application of this technique we provide a public Python package, called PyStacker.

Spectroscopic phase curves of transiting hot Jupiters are spectral measurements at multiple orbital phases, giving a set of disc-averaged spectra that probe multiple hemispheres. By fitting model phase curves to observations, we can constrain the atmospheric properties of hot Jupiters such as molecular abundance, aerosol distribution and thermal structure, which offer insights into their dynamics, chemistry, and formation. In this work, we propose a novel 2D temperature scheme consisting of a dayside and a nightside to retrieve information from near-infrared phase curves, and apply the scheme to phase curves of WASP-43b observed by HST/WFC3 and Spitzer/IRAC. In our scheme, temperature is constant on isobars on the nightside and varies with cos$^n$(longitude/$\epsilon$) on isobars on the dayside, where $n$ and $\epsilon$ are free parameters. We fit all orbital phases simultaneously using the radiative transfer code NEMESISPY coupled to a Bayesian inference code. We first validate the performance of our retrieval scheme with synthetic phase curves generated from a GCM, and find our 2D scheme can accurately retrieve the latitudinally-averaged thermal structure and constrain the abundance of H$_2$O and CH$_4$. We then apply our 2D scheme to the observed phase curves of WASP-43b and find: (1) the dayside temperature-pressure profiles do not vary strongly with longitude and are non-inverted; (2) the retrieved nightside temperatures are extremely low, suggesting significant nightside cloud coverage; (3) the H$_2$O volume mixing ratio is constrained to $5.6\times10^{-5}$--$4.0\times10^{-4}$, and we retrieve an upper bound for CH$_4$ at $\sim$10$^{-6}$.

We present a new model of anisotropic cosmic ray propagation in the Milky Way, where cosmic rays are injected at discrete transient sources in the disc and propagated in the Galactic magnetic field. In the framework of our model, we show that the cosmic ray spectrum is time-dependent and space-dependent around the energy of the knee. It has a major contribution of one or a few nearby recent sources at any given location in the Galaxy, in particular at the position of the Solar system. We find that the distribution of $\sim$ PeV cosmic rays in our Galaxy is significantly clumpy and inhomogeneous, and therefore substantially different from the smoother distribution of GeV cosmic rays. Our findings have important implications for the calculation and future interpretation of the diffuse Galactic gamma-ray and neutrino fluxes at very high energies.

Despite having data for over 10^9 stars from Gaia, only less than 10^4 star clusters and candidates have been discovered. Particularly, distant star clusters are rarely identified, due to the challenges posed by heavy extinction and great distance. However, Gaia data has continued to improve, enabling even fainter cluster members to be distinguished from field stars. In this work, we will introduce a star cluster search method based on the DBSCAN algorithm; we have made improvements to make it better suited for identifying clusters on dimmer and more distant stars. After removing member stars of known Gaia-based clusters, we have identified 2086 objects with |b|<10 deg, of which 1488 are highly reliable open star clusters, along with 569 candidates, 28 globular cluster candidates and 1 irregular galaxy IC 10 at low Galactic latitudes. We found that the proper motion of IC10 is similar yet slightly different from the water maser observations, which is an important result for the comparison with Gaia and VLBA. Besides, when compared with the star clusters appearing in Gaia DR2/EDR3, we have found nearly three times as many new objects above a distance of 5 kpc, including hundreds of them above Av > 5 mag. And it has enabled us to detect a higher number of old clusters, over a billion years old, that are difficult to detect due to observational limitations. Our findings significantly expand the remote cluster sample and enhance our understanding of the limits of Gaia DR3 data in stellar aggregates research. The full figure set for 2085 clusters can be seen in \url{https://nadc.china-vo.org/res/r101258/}

The Weakly Interacting Massive Particles(WIMPs) have long been the favored CDM candidate in the standard $\Lambda$CDM model. However, owing to great improvement in the experimental sensitivity in the past decade, some parameter space of the SUSY-based WIMP model is ruled out. In addition, WIMP as the CDM particle is also at variance with other astrophysical observables at small scales. We consider a model that addresses both these issues. In the model, the WIMP decays into a massive particle and radiation. We study the background evolution and the first order perturbation theory (coupled Einstein-Boltzmann equations) for this model and show that the dynamics can be captured by a single parameter $r=m_L/q$, which is the ratio of the lighter mass and the comoving momentum of the decay particle. We incorporate the relevant equations in the existing Boltzmann code CLASS to compute the matter power spectra and CMB angular power spectra. The decaying WIMP model is akin to a non-thermal Warm Dark Matter(WDM) model and suppresses matter power at small scales, which could alleviate several issues that plague the CDM model. We compare the predictions of the model with CMB, galaxy clustering, and high-z HI data. Both these data sets yield $r\gtrsim 10^6$, which can be translated into the bounds on other parameters. In particular, we obtain the following lower bounds on the self-annihilation cross-section of WIMPs $\sigma$, and the lighter mass $m_L$: $\sigma\gtrsim 5\times 10^{-44}\,\rm cm^2$ and $m_L\gtrsim 2.4\,\rm keV$. The lower limit on $m_L$ is comparable to constraints on the mass of thermally-produced WDM particle. The limit on the self-annihilation cross-section greatly expands the available parameter space as compared to the stable WIMP scenario.

Observations of the vicinity of a variety of galactic gamma-ray sources have indicated a local suppression of diffusivity of cosmic rays by up to three orders of magnitude. However, the impact of these low-diffusion zones on \emph{global} properties of cosmic-ray transport is however only poorly understood. Here, we argue that cosmic-ray nuclear ratios, like the boron-to-carbon ratio and relative abundances of Beryllium isotopes are sensitive to the filling fraction of such low-diffusion zones and hence their measurements can be used to constrain the typical sizes and ages of such regions. We have performed a careful parameter study of a cosmic-ray transport model that allows for different diffusion coefficients $\kappa_{\mathrm{disk}}$ and $\kappa_{\mathrm{halo}}$ in the galactic disk and halo, respectively. Making use of preliminary data from the AMS-02 experiment on the ratio of Beryllium isotopes, we find a $3.5 \sigma$ preference for a suppression of the diffusion coefficient in the disk with a best-fit value of $\kappa_{\mathrm{disk}}/\kappa_{\mathrm{halo}} = 0.20^{+0.10}_{-0.06}$. We forecast that with upcoming data from the HELIX balloon experiment, the significance could increase to $6.8 \sigma$. Adopting a coarse-graining approach, we find that such a strong suppression could be realised if the filling fraction of low-diffusion zones in the disk was $\sim 66 \, \%$. We conclude that the impact of regions of suppressed diffusion might be larger than usually assumed and ought to be taken into account in models of Galactic cosmic ray transport.

Ultralight dark matter (ULDM) is an interesting alternative to the cold dark matter (CDM) paradigm. Due to the extremely low mass of the constituent particle ($\sim 10^{-22}$ eV), ULDM can exhibit quantum effects up to kiloparsec scales. In particular, runaway collapse in the centres of ULDM halos is prevented by quantum pressure, providing a possible resolution to the 'core-cusp problem' of CDM. However, the the detailed relationship between the ULDM core mass and that of the overall halo is poorly understood. We simulate the collapse of both spherical and aspherical isolated ULDM overdensities using AxioNyx, finding that the central cores of collapsed halos undergo sustained oscillatory behaviour which affects both their peak density and overall morphology. The variability in core morphology increases with the asphericity of the initial overdensity and remnants of initial asphericity persist long after collapse. Furthermore, the peak central densities are higher in spherical configurations. Consequently, astrophysically realistic halos may exhibit substantial departures from theoretical core-halo profiles and we would expect a significant variance of the properties of halos with the same mass.

While radial velocity surveys have demonstrated that the population of gas giants peaks around $3~\text{au}$, the most recent high-contrast imaging surveys have only been sensitive to planets beyond $\sim~10~\text{au}$. Sensitivity at small angular separations from stars is currently limited by the variability of the point spread function. We demonstrate how moderate-resolution integral field spectrographs can detect planets at smaller separations ($\lesssim~0.3$ arcseconds) by detecting the distinct spectral signature of planets compared to the host star. Using OSIRIS ($R$ $\approx$ 4000) at the W. M. Keck Observatory, we present the results of a planet search via this methodology around 20 young targets in the Ophiuchus and Taurus star-forming regions. We show that OSIRIS can outperform high-contrast coronagraphic instruments equipped with extreme adaptive optics and non-redundant masking in the $0.05-0.3$ arcsecond regime. As a proof of concept, we present the $34\sigma$ detection of a high-contrast M dwarf companion at $\approx0.1$" with a flux ratio of $\approx0.92\%$ around the field F2 star HD 148352. We developed an open-source Python package, breads, for the analysis of moderate-resolution integral field spectroscopy data in which the planet and the host star signal are jointly modeled. The diffracted starlight continuum is forward-modeled using a spline model, which removes the need for prior high-pass filtering or continuum normalization. The code allows for analytic marginalization of linear hyperparameters, simplifying posterior sampling of other parameters (e.g., radial velocity, effective temperature). This technique could prove very powerful when applied to integral field spectrographs like NIRSpec on the JWST and other upcoming first-light instruments on the future Extremely Large Telescopes.

We have collected a catalog of 1095 debris disks with properties and classification (resolved, planet, gas) information. From the catalog, we defined a less biased sample with 612 objects and presented the distributions of their stellar and disk properties to search for correlations between disks and stars. We found debris disks were widely distributed from B to M-type stars while planets were mostly found around solar-type stars, gases were easier to detect around early-type stars and resolved disks were mostly distributed from A to G- type stars. The fractional luminosity dropped off with stellar age and planets were mostly found around old stars while gas-detected disks were much younger. The dust temperature of both one-belt systems and cold components in two-belt systems increased with distance while decreasing with stellar age. In addition, we defined a less biased planet sample with 211 stars with debris disks but no planets and 35 stars with debris disks and planets and found the stars with debris disks and planets had higher metallicities than stars with debris disks but no planets. Among the 35 stars with debris disks and planets, we found the stars with disks and cool Jupiters were widely distributed with age from 10 Myr to 10 Gyr and metallicity from -1.56 to 0.28 while the other three groups tended to be old (> 4Gyr) and metal-rich (> -0.3). Besides, the eccentricities of cool Jupiters are distributed from 0 to 0.932 wider than the other three types of planets (< 0.3).

Recently, the first JWST measurement of thermal emission from a rocky exoplanet was reported. The inferred dayside brightness temperature of TRAPPIST-1 b at 15 $\mu$m is consistent with the planet having no atmosphere and therefore no mechanism by which to circulate heat to its nightside. In this Letter, we compare the measured secondary eclipse depth of TRAPPIST-1 b to predictions from a suite of self-consistent radiative-convective equilibrium models in order to quantify the maximum atmospheric thickness consistent with the observation. We find that plausible atmospheres (i.e., those that contain at least 100 ppm CO$_2$) with surface pressures greater than 0.01 bar (0.1 bar) are ruled out at 1$\sigma$ (3$\sigma$), regardless of the choice of background atmosphere. Thicker atmospheres of up to 10 bar (100 bar) at 1$\sigma$ (3$\sigma$) are only allowed if the atmosphere lacks any strong absorbers across the mid-IR wavelength range, a scenario that we deem unlikely. We additionally model the emission spectra for bare-rock planets of various compositions. We find that a variety of silicate surfaces match the measured eclipse depth to within 1$\sigma$, and the best-fit grey albedo is $0.02 \pm 0.11$. We conclude that planned secondary eclipse observations at 12.8 $\mu$m will serve to validate the high observed brightness temperature of TRAPPIST-1 b, but are unlikely to further distinguish among the consistent atmospheric and bare-rock scenarios.

Millions of quasar spectra will be collected by the Dark Energy Spectroscopic Instrument (DESI), leading to a four-fold increase in the number of known quasars. High accuracy quasar classification is essential to tighten constraints on cosmological parameters measured at the highest redshifts DESI observes ($z>2.0$). We present the spectral templates for identification and redshift estimation of quasars in the DESI Year 1 data release. The quasar templates are comprised of two quasar eigenspectra sets, trained on spectra from the Sloan Digital Sky Survey. The sets are specialized to reconstruct quasar spectral variation observed over separate yet overlapping redshift ranges and, together, are capable of identifying DESI quasars from $0.05 < z <7.0$. The new quasar templates show significant improvement over the previous DESI quasar templates regarding catastrophic failure rates, redshift precision and accuracy, quasar completeness, and the contamination fraction in the final quasar sample.

We present an improved analytical model to predict the Lyman-alpha forest at the field level in redshift space from the dark matter field, expanding upon the widely-used Fluctuating Gunn-Peterson approximation (FGPA). In particular, we introduce the dependence on the cosmic web environment (knots, filaments, sheets, voids) in the model, thereby effectively accounting for non-local bias. Furthermore, we include a detailed treatment of velocity bias in the redshift space distortions modelling, allowing the velocity bias to be cosmic-web dependent. We find evidence for a significant difference of the same model parameters in different environments, suggesting that for the investigated setup the simple standard FGPA is not able to adequately predict the Lyman-alpha forest in the different cosmic web regimes. We reproduce the summary statistics of the reference cosmological hydrodynamic simulation we use for comparison, yielding accurate mean transmitted flux, probability distribution function, 3D power spectrum, and bispectrum. In particular, we achieve maximum deviation and average deviations accuracy in the Lyman-alpha forest 3D power spectrum of $\sim 3\%$ and $\sim 0.1\%$ up to $k\sim 0.4 \, h \, {\rm Mpc}^{-1}$, $\sim 5\%$ and $\sim 1.8\%$ up to $k \sim 1.4 \, h \, {\rm Mpc}^{-1}$. Our new model outperforms previous analytical efforts to predict the Lyman-alpha forest at the field level in all the probed summary statistics, and has the potential to become instrumental in the generation of fast accurate mocks for covariance matrices estimation in the context of current and forthcoming Lyman-alpha forest surveys.

We develop a frequency-domain method for calculating the self-force acting on a scalar charge on a fixed scattering geodesic in Schwarzschild spacetime. Existing frequency-domain methods, which are tailored for bound orbits, are inadequate here for several reasons. One must account for the continuous spectrum in the scattering problem, deal with slowly-convergent radial integrals that are hard to evaluate numerically, and confront the inapplicability of the standard self-force method of "extended homogeneous solutions", which only works for compactly supported sources. We tackle each of these issues in turn, and then present a full numerical implementation, in which we calculate the self-force correction to the scatter angle due to scalar-field back-reaction. We perform a range of internal validation tests, as well as ones based on comparison with existing time-domain results. We discuss the merits and remaining limitations of our method, and outline directions for future work.

The gravitational positivity bound gives quantitative ``swampland'' constraints on low-energy effective theories inside theories of quantum gravity. We give a comprehensive discussion of this bound for those interested in applications to phenomenological model building. We present a practical recipe for deriving the bound, and discuss subtleties relevant for realistic models. As an illustration, we study the positivity bound on the scattering of the massive gauge bosons in the Higgs/St\"{u}ckelberg mechanism. Under certain assumptions on gravitational amplitudes at high energy, we obtain a lower bound $m_{V} \gtrsim \Lambda_\mathrm{UV}^2 /g M_\mathrm{Pl}$ on the gauge boson mass $m_V$, where $g$ is the coupling constant of the gauge field, $M_\mathrm{Pl}$ is the reduced Planck mass and $\Lambda_\mathrm{UV}$ is the ultraviolet cutoff of the effective field theory. This bound can strongly constrain new physics models involving a massive gauge boson.

In recent times, astounding observations of both over- and under-luminous type Ia supernovae have emerged. These peculiar observations hint not only at surpassing the Chandrasekhar limit but may also suggest potential modifications in the physical attributes of their progenitors, such as their cooling rate. This, in turn, can influence their temporal assessments and provide a compelling explanation for these intriguing observations. In this spirit, we investigate here the cooling process of white dwarfs in $f(R,T)$ gravity with the simplest model $f(R,T) = R + \lambda T$, where $\lambda$ is the model parameter. Our modelling suggests that the cooling timescale of white dwarfs exhibits an inverse relationship with the model parameter $\lambda$, which implies that for identical initial conditions, white dwarfs in $f(R,T)$ gravity cool faster. This further unveils that in the realm of $f(R,T)$ gravity, the energy release rate for white dwarfs increases as $\lambda$ increases. Furthermore, we also report that the luminosity of the white dwarfs also depends on $\lambda$ and an upswing in $\lambda$ leads to an amplification in the luminosity, and consequently a larger white dwarf in general relativity can exhibit comparable luminosity to a smaller white dwarf in $f(R,T)$ gravity.

Dark matter is undoubtedly one of the fundamental, albeit unknown, components of the standard cosmological model. The failure to detect WIMPs, the most promising candidate particle for cold dark matter, actually opens the way for the exploration of viable alternatives, of which ultralight bosonic particles with masses $\sim 10^{-21}$ eV represent one of the most encouraging. N-body simulations have shown that such particles form solitonic cores in the innermost parts of virialized galactic halos that are supported by internal quantum pressure on characteristic $\sim$kpc de Broglie scales. In the Galaxy, this halo region can be probed by means of S-stars orbiting the supermassive black hole Sagittarius A* to unveil the presence of such a solitonic core and, ultimately, to bound the boson mass $m_\psi$. Employing a Monte Carlo Markov Chain algorithm, we compare the predicted orbital motion of S2 with publicly available data and set an upper bound $m_\psi \lesssim 3.2\times 10^{-19}$ eV on the boson mass, at 95% confidence level. When combined with other galactic and cosmological probes, our constraints help to reduce the allowed range of the bosonic mass to $(2.0 \lesssim m_\psi \lesssim 32.2)\times 10^{-20}$ eV, at the 95% confidence level, which opens the way to precision measurements of the mass of the ultralight bosonic dark matter.

The metastable hypermassive neutron star produced in the coalescence of two neutron stars can copiously produce axions that radiatively decay into $\mathcal{O}(100)$~MeV photons. These photons can form a fireball with characteristic temperature smaller than $1\rm\, MeV$. By relying on X-ray observations of GW170817/GRB 170817A with CALET CGBM, Konus-Wind, and Insight-HXMT/HE, we present new bounds on the axion-photon coupling for axion masses in the range $1$-$400\,\rm MeV$. We exclude couplings down to $5\times 10^{-11}\,\rm GeV^{-1}$, complementing and surpassing existing constraints. Our approach can be extended to any feebly-interacting particle decaying into photons.