Papers for Wednesday, Jul 31 2024

More than fifty years have been elapsed from the first discovery of a Gamma Ray Burst (GRB) with the American Vela satellites, and more than twenty-five years from the discovery with the BeppoSAX satellite of the first X-ray afterglow of a GRB. Thanks to the afterglow discovery and to the possibility given to the optical and radio astronomers to discover the GRB optical counterparts, the long-time mystery about the origin of these events was solved. Now we know that GRBs are huge explosions, mainly ultra relativistic jets, in galaxies at cosmological distances. Starting from the first GRB detection with the Vela satellites, I will review the story of these discoveries, those obtained with BeppoSAX, the contribution to GRBs by other satellites and ground experiments, among them Venera, Compton Gamma Ray Observatory, HETE2, Swift, Fermi, AGILE, MAGIC, H.E.S.S., which were, and some of them are still, very important for the study of GRB properties. Then I will review the main results obtained thus far and the still open problems and prospects of the GRB astronomy.

The National Science Foundation's Daniel K. Inouye Solar Telescope (DKIST) will provide high-resolution, multi-line spectropolarimetric observations that are poised to revolutionize our understanding of the Sun. Given the massive data volume, novel inference techniques are required to unlock its full potential. Here, we provide an overview of our "SPIn4D" project, which aims to develop deep convolutional neural networks (CNNs) for estimating the physical properties of the solar photosphere from DKIST spectropolarimetric observations. We describe the magnetohydrodynamic (MHD) modeling and the Stokes profile synthesis pipeline that produce the simulated output and input data, respectively. These data will be used to train a set of CNNs that can rapidly infer the four-dimensional MHD state vectors by exploiting the spatiotemporally coherent patterns in the Stokes profile time series. Specifically, our radiative MHD model simulates the small-scale dynamo actions that are prevalent in quiet-Sun and plage regions. Six cases with different mean magnetic fields have been conducted; each case covers six solar-hours, totaling 109 TB in data volume. The simulation domain covers at least $25\times25\times8$ Mm with $16\times16\times12$ km spatial resolution, extending from the upper convection zone up to the temperature minimum region. The outputs are stored at a 40 s cadence. We forward model the Stokes profile of two sets of Fe I lines at 630 and 1565 nm, which will be simultaneously observed by DKIST and can better constrain the parameter variations along the line of sight. The MHD model output and the synthetic Stokes profiles are publicly available.

We present follow-up SCExAO/CHARIS $H$ and $K$-band (R $\sim$ 70) high-contrast integral field spectroscopy and Keck/NIRC2 photometry of directly-imaged brown dwarf companion HD 33632 Ab and new radial-velocity data for the system from the SOPHIE spectrograph, complemented by Hipparcos and Gaia astrometry. These data enable more robust spectral characterization compared to lower-resolution spectra from the discovery paper and more than double the available astrometric and radial-velocity baseline. HD 33632 Ab's spectrum is well reproduced by a field L8.5--L9.5 dwarf. Using the Exo-REM atmosphere models, we derive a best-fit temperature, surface gravity and radius of $T_{\rm eff}$ = 1250 $K$, log(g) = 5, and $R$ = 0.97 $R_{\rm J}$ and a solar C/O ratio. Adding the SOPHIE radial-velocity data enables far tighter constraints on the companion's orbital properties (e.g. $i$=${46.6}_{-5.7}^{+2.9}$$^{o}$) and dynamical mass (${51.7}_{-2.5}^{+2.6}$$M_{\rm J}$) than derived from imaging data and \textit{Gaia} eDR3 astrometry data alone. HD 33632 Ab should be a prime target for multi-band imaging and spectroscopy with the James Webb Space Telescope and the Roman Space Telescope's Coronagraphic Instrument, shedding detailed light on HD 33632 Ab's clouds and chemistry and providing a key reference point for understanding young exoplanet atmospheres.

We present initial results from the Polarisation Sky Survey of the Universe's Magnetism (POSSUM), analysing 22,817 Faraday Rotation Measures (RMs) with median uncertainties of 1.2 rad m^-2 across 1,520 square degrees to study magnetised gas associated with 55 nearby galaxy groups (z less than 0.025) with halo masses between 10^12.5 and 10^14.0 M_sun. We identify two distinct gas phases: the Intragroup Medium (IGrM) within 0-2 splashback radii and the Warm-Hot Intergalactic Medium (WHIM) extending from 2 to 7 splashback radii. These phases enhance the standard deviation of residual (i.e., Galactic foreground RM-subtracted) RMs by 6.9 +/- 1.8 rad m^-2 and 4.2 +/- 1.2 rad m^-2, respectively. Estimated magnetic field strengths are several microGauss within the IGrM and 0.1-1 microGauss in the WHIM. We estimate the plasma beta in both phases and show that magnetic pressure might be more dynamically important than in the ICM of more massive clusters or sparse cosmic web filaments. Our findings indicate that "missing baryons" in the WHIM likely extend beyond the gravitational radii of group-mass halos to Mpc scales, consistent with large-scale, outflow-driven "magnetised bubbles" seen in cosmological simulations. We demonstrate that RM grids are an effective method for detecting magnetised thermal gas at galaxy group interfaces and within the cosmic web. This approach complements X-ray and Sunyaev-Zel'dovich effect methods, and when combined with Fast Radio Burst Dispersion Measures, data from the full POSSUM survey, comprising approximately a million RMs, will allow direct magnetic field measurements to further our understanding of baryon circulation in these environments and the magnetised universe.

We report the serendipitous discovery of a faint (M_UV > -12.2), low-metallicity (Z ~ 0.02 Zsun), ionizing source (dubbed T2c) with a spectroscopic redshift of z=6.146. T2c is part of a larger structure amplified by the Hubble Frontier Field galaxy cluster MACSJ0416, and was observed with JWST/NIRSpec IFU. Stacking the short-wavelength NIRCam data reveals no stellar continuum detection down to a magnitude limit of m_UV ~ 31.0 (3 sigma). However, prominent Hb, [OIII]4959,5007, and Ha emissions are detected, with equivalent widths exceeding 200A, 800A, and 1300A (3 sigma), respectively. The corresponding intrinsic (magnification-corrected x23 +/- 3) ultraviolet and optical rest-frame magnitudes exceed 34.4 and 33.9 (corresponding to M_uv and M_opt fainter than -12.2 and -12.8, at lambda_rest ~ 2000A and ~5000A, respectively), suggesting a stellar mass lower than a few 10^4 Msun under an instantaneous burst scenario. The inferred ionizing photon production efficiency (xi_ion) is high, xi_ion >~ 26.08(25.86) 3(5)sigma, assuming no dust attenuation and no Lyman continuum leakage, indicating the presence of massive stars despite the low mass of the object. The very poor sampling of the initial mass function at such low mass star-forming complex suggests that the formation of very massive stars might be favored in very low metallicity environments. T2c is surrounded by Balmer and weak oxygen emission on a spatial scale of a few hundred parsecs after correcting for lensing effects. This system resembles an HII region potentially powered by currently undetected, extremely efficient, low-metallicity star complexes or clusters. We propose that massive O-type stars populate this low-mass and metallicity high-redshift satellites, likely caught in an early and short formation phase, contributing to the ionization of the surrounding medium.

We present a model to estimate the average primary masses, companion mass ranges, the inclination limit for recognizing a rotational variable, and the primary mass spreads for populations of binary stars. The model fits a population's binary mass function distribution and allows for a probability that some mass functions are incorrectly estimated. Using tests with synthetic data, we assess the model's sensitivity to each parameter, finding that we are most sensitive to the average primary mass and the minimum companion mass, with less sensitivity to the inclination limit and little to no sensitivity to the primary mass spread. We apply the model to five populations of binary spotted rotational variables identified in ASAS-SN, computing their binary mass functions using RV data from APOGEE. Their average primary mass estimates are consistent with our expectations based on their CMD locations ($\sim0.75 M_{\odot}$ for lower main sequence primaries and $\sim 0.9$--$1.2 M_{\odot}$ for RS CVn and sub-subgiants). Their companion mass range estimates allow companion masses down to $M_2/M_1\simeq0.1$, although the main sequence population may have a higher minimum mass fraction ($\sim0.4$). We see weak evidence of an inclination limit $\gtrsim50^{\circ}$ for the main sequence and sub-subgiant groups and no evidence of an inclination limit in the other groups. No groups show strong evidence for a preferred primary mass spread. We conclude by demonstrating that the approach will provide significantly better estimates of the primary mass and the minimum mass ratio and reasonable sensitivity to the inclination limit with 10 times as many systems.

We present a detailed JWST/NIRSpec and NIRCam analysis of a gravitationally-lensed galaxy ($\rm \mu=17-21$) at redshift 6.14 magnified by the Hubble Frontier Field galaxy cluster MACS J0416. The target galaxy is overall a typical compact and UV-faint ($\rm M_{UV}=-17.8$) Lyman-$\alpha$ emitter; yet, the large magnification allows the detailed characterisation of structures on sub-galactic (down to few parsec) scales. Prominent optical $\rm H\alpha$, $\rm H\beta$ and [OIII]$\lambda\lambda4959,5007$ lines are spatially resolved with the high spectral resolution grating (G395H, R~2700), with large equivalent widths, EW($\rm H\beta$+[OIII])$\gtrsim1000$ Å, and elevated ionising photon production efficiencies $\rm log(\xi_{ion}/erg^{-1}Hz)=25.2-25.7$. NIRCam deep imaging reveals the presence of compact rest-UV bright regions along with individual star clusters of sizes $\rm R_{eff}=3-8~pc$ and masses $\rm M\sim2\cdot10^5-5\cdot10^{6}~M_\odot$ These clusters are characterised by steep UV slopes, $\rm\beta_{UV}\lesssim-2.5$, in some cases associated with a dearth of line emission, indicating possible leaking of the ionising radiation, as also supported by a Lyman-$\rm \alpha$ emission peaking at $\rm \sim100~km~s^{-1}$ from the systemic redshift. While the entire system is characterised by low-metallicity, $\sim0.1~Z_\odot$, the NIRSpec-IFU map also reveals the presence of a low-luminosity, metal-poor region with $\rm Z\lesssim2\%~Z_\odot$, barely detected in NIRCam imaging; this region is displaced by $\rm >200~pc$ from one of the UV brightest structures of the system, and it would have been too faint to detect if not for the large magnification of the system.

We present {\sc Venice}, an operator splitting algorithm to integrate a numerical model on a hierarchy of timescales. {\sc Venice} allows a wide variety of different physical processes operating a different scales to be coupled on individual and adaptive time-steps. It therewith mediates the development of complex multi-scale and multi-physics simulation environments with a wide variety of independent components. The coupling between various physical models and scales is dynamic, and realized through (Strang) operators splitting using adaptive time steps. We demonstrate the functionality and performance of this algorithm using astrophysical models of a stellar cluster, first coupling gravitational dynamics and stellar evolution, then coupling internal gravitational dynamics with dynamics within a galactic background potential, and finally combining these models while also introducing dwarf galaxy-like perturbers. These tests show numerical convergence for decreasing coupling timescales, demonstrate how {\sc Venice} can improve the performance of a simulation by shortening coupling timescales when appropriate, and provide a case study of how {\sc Venice} can be used to gradually build up and tune a complex multi-physics model. Although the examples couple complete numerical models, {\sc Venice} can also be used to efficiently solve systems of stiff differential equations.

We analyze a long-lived hyper-massive neutron star merger remnant (post-merger lifetime $>250$ ms) that has been obtained via large eddy simulations with a gradient subgrid-scale model. We find a clear helicoidal magnetic field structure that is governed by the toroidal component of the magnetic field. Although no jet emerges during the simulation time, we observe at late times a significant increase of the poloidal component of the magnetic field at all scales. We also compare with the results of several binary neutron star simulations with moderate resolution of $120$~m, that are evolved up to $50$~ms after the merger, which differ in terms of the initial topology and strength of the magnetic field. We find that the best choice is an isotropic small-scale magnetic field distribution that mimics the turbulent state that generically develops during the merger. This initial configuration reaches a closer agreement with our high-resolution simulation results than the purely dipolar large-scale fields that are commonly employed in these type of simulations. This provides a recipe to perform such simulations avoiding the computationally expensive grids required to faithfully capture the amplification of the magnetic field by Kelvin-Helmholtz instabilities.

The Hartmann test is a method used to measure the wavefront error in a focal optical system, wherein a mask with a pattern of small holes is placed at the system's aperture stop. By taking an image at a defocused plane, the differences between the ideal and real positions of the reimaged holes (called the transverse ray aberrations) can be measured, which can then be used to estimate the wavefront error. However, the Hartmann test is usually used with an on-axis field. In this paper, we present a wavefront sensing method which generalizes the classical Hartmann test for off-axis field angles and arbitrary reference wavefronts. Our method involves taking images at two defocused planes, and then using the real reimaged hole positions on both planes to estimate the trajectories of rays from the system's exit pupil, at which the reference wavefront is situated. We then propagate the rays forward from the reference wavefront to one of the two defocused planes, in order to find the ideal reimaged hole positions, from which we can compute transverse ray aberrations. We derive and solve a pair of nonlinear partial differential equations relating transverse ray aberrations to wavefront error, using Zernike decomposition and nonlinear least squares. Our method has been verified on simulated data from the 7-lens f/2.25 red camera system of the GMT-Consortium Large Earth Finder (G-CLEF), a high resolution optical echelle spectrograph which will be a first light instrument for the Giant Magellan Telescope (GMT).

Superluminous supernovae (SLSNe) radiate $\gtrsim 10-100$ times more energy than ordinary stellar explosions, implicating a novel power source behind these enigmatic events. One frequently discussed source, particularly for hydrogen-poor (Type I) SLSNe, is a central engine such as a millisecond magnetar or accreting black hole. Both black hole and magnetar engines are expected to channel a fraction of their luminosity into a collimated relativistic jet. Using 3D relativistic hydrodynamical simulations, we explore the interaction of a relativistic jet, endowed with a luminosity $L_{\rm j}\approx10^{45}\,{\rm erg\,s^{-1}}$ and duration $t_{\rm eng} \approx$ days compatible with those needed to power SLSNe, launched into the envelope of the exploding star. The jet successfully breaks through the expanding ejecta and its shocked cocoon powers ultraviolet/optical emission lasting several days after the explosion and reaching a peak luminosity $\gtrsim 10^{44}$ erg s$^{-1}$, corresponding to a sizable fraction of $L_{\rm j}$. This high radiative efficiency is the result of the modest adiabatic losses the cocoon experiences due to the low optical depths of the enlarged ejecta at these late times, e.g., compared to the more compact stars in gamma-ray bursts. The luminosity and temperature of the cocoon emission match those of the ``bumps'' in SLSNe light curves observed weeks prior to the optical maximum in many SLSNe. Confirmation of jet breakout signatures by future observations (e.g., days to weeks-long internal X-ray emission from the jet for on-axis observers, spectroscopy confirming large photosphere velocities $v/c \gtrsim 0.1$, or detection of a radio afterglow) would offer strong evidence for central engines powering SLSNe.

It is often believed that isolated, compact groups of galaxies (CGs) are special systems. But CGs have rarely been compared to regular groups, which are selected very differently. We study the global properties of a subsample of 80 groups of four members (CG4s), as well as their correlations, within the HMCG Hickson-like sample of compact groups, and compare them with those of the regular groups of Lim et al.. The latter are split into three control samples, complete in redshift and magnitude: one with the four brightest members, one with the four closest members to the brightest group galaxy (BGG) with less than three magnitudes in range (Control4Cs), and one with exactly four members (RG4s). The vast majority of CG4s are located within regular groups, and a large preponderance of the BGGs of these CG4s are the same as those of their host groups. CG4s are smaller than the groups of all other samples and more luminous than RG4s, both results as expected from their selection as high surface brightness systems. However, CG4s have similar luminosities as Control4Cs. CG4s also have higher velocity dispersions, probably because of a too permissive redshift accordance criterion. The BGGs of CG4s are not more dominant in luminosity than those of RG4s, but they are significantly more offset relative to the group size, because the Lim groups are built around their BGGs. In summary, compact groups have similar properties to the regular groups of four galaxies and to the cores of regular groups, once selection criteria of CGs are considered. A large fraction of CGs are the cores of regular groups, isolated on the sky by construction, but rarely isolated in real space (from simulations), indicating that they are often plagued by chance alignments of host groups galaxies along the line of sight.

Intrinsic variability was searched for in arrival times of six gamma-ray bursts (GRBs) at high energies -- between 30 MeV and 2 GeV -- detected by the Fermi satellite's Large Area Telescope (LAT). The GRBs were selected from the Fermi LAT catalog with preference for events with numerous photons, a strong initial pulse, and measured redshifts. Three long GRBs and three short GRBs were selected and tested. Two different variability-detection algorithms were deployed, one counting photons in pairs, and the other multiplying time gaps between photons. In both tests, a real GRB was compared to 1000 Monte-Carlo versions of itself smoothed over a wide range of different timescales. The minimum detected variability timescales for long bursts (GRB 080916C, GRB 090926A, GRB 131108A) was found to be (0.005, 10.0, 10.0) seconds for the photon pair test and (2.0, 20.0, 10.0) seconds for the time-gap multiplication test. Additionally, the minimum detected variability timescales for the short bursts (GRB 090510, GRB 140619B, GRB 160709A) was found to be (0.05, 0.01, 20.0) seconds for the photon pair test and (0.05, 0.01, 20.0) seconds for the gap multiplication test. Statistical uncertainties in these times are about a factor of 2. The durations of these variability timescales may be used to constrain the geometry, dynamics, speed, cosmological dispersion, Lorentz-invariance violations, weak equivalence principle violations, and GRB models.

Gravitational coupling between a protoplanetary disc and an embedded eccentric planet is an important, long-standing problem, which has been not yet been conclusively explored. Here we study the torque and associated orbital evolution of an eccentric planet in a two-dimensional disc via the semi-analytical, fully global linear approach. Our methodology has the advantage that the spatial structure of the density waves launched by the planet is solved for fully. This allows us to account for the possibility of torque excitation over an extended radial interval for each Fourier harmonic of the perturbation, as opposed to earlier approximate treatments localized around Lindblad and corotation resonances. We systematically explore the torque behaviour across the space of disc properties (assuming power law profiles for the disc surface density and temperature), including the aspect ratio. Crucially, we examine the torque variation as the orbital eccentricity becomes supersonic relative to the gas motion (when planetary eccentricity is of order the disc aspect ratio), finding that the torque robustly reverses its sign near this transition. We also find that for shallow surface density gradients planetary migration may become outwards beyond this transition, although the rapid eccentricity damping (which is typically $\sim 10^2$ times faster than the orbital migration rate) would quickly restore inwards migration as the planet circularizes. Our self-consistently computed torques are in qualitative agreement with past numerical studies of eccentric planet-disc coupling. We provide our torque data for different disc parameters to the community for future testing and implementation in population synthesis studies.

We examine the local density environments around 67 quasars at $z\sim3$, by combining the imaging data of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and Canada-France-Hawaii Telescope Large Area U-band Survey (CLAUDS) over about 20 deg$^{2}$. Our measurements exploit $U$-dropout galaxies in the vicinities of quasars taken from Sloan Digital Sky Survey (SDSS). We find that the quasars have the indistinguishable surrounding density distribution from the $U$-dropout galaxies, and that three quasars are associated with protocluster candidates within a projected separation of 3 arcmin. According to a halo evolutionary model, our results suggest that quasars at this epoch occupy haloes with a typical mass of $1.3^{+1.4}_{-0.9} \times 10^{13} h^{-1} \mathrm{M_{\odot}}$. We also investigate the dependence of the local galaxy overdensity on ultraviolet (UV) luminosities, black hole (BH) masses, and proximity zone sizes of the quasars, but no statistically-significant correlation was found. Finally, we find that the local density of faint $U$-dropout galaxies are lower than that of bright $U$-dropout galaxies within a projected distance of $0.51\pm0.05$ physical Mpc, where the quasar UV radiation is 30 times intenser than backgroung UV radiation. We argue that photoevaporation may suppress galaxy formation at short distances where the quasar UV intensity is strong, even in massive haloes.

Star-forming regions, characterized by dense environments, experience frequent encounters that significantly influence binary systems, leading to their hardening, softening, or ionization. We extend the Hut \& Bahcall formalism to derive an analytical expression for the ionization cross-section in extreme mass ratio binary systems, allowing us to investigate the orbital evolution and population dynamics of binary planets and binary brown dwarfs in star clusters, while considering ongoing binary system formation. Our findings reveal that for low-mass soft binaries, the semi-major axis distribution asymptotes to a universal power law between $\propto a^{-8/3}$ and $\propto a^{-5/3}$ over the derived ionization timescale. We also discuss the implications of our results for the candidate Jupiter-mass binary objects putatively reported in the Trapezium cluster. We demonstrate that if their existence is verified, they likely form continuously with a spectrum proportional to $a^{1}$, aligning better with the ejection mechanism than with the in-situ formation mechanism, which predicts a distribution roughly proportional to $a^{-1}$. However, this implies an impractically high ejection formation rate. Alternatively, if these objects are binary brown dwarfs, continuous in-situ formation ($\propto a^{-1}$) with an initial minimal semi-major axis around 20 AU and a formation rate of 100 Myr$^{-1}$ plausibly matches the observed number of single objects, binary number, binary fraction, and semi-major axis distribution.

Recent observations made using the James Webb Space Telescope have identified a number of high-redshift galaxies that are unexpectedly luminous. In light of this, it is clear that a more detailed understanding of the high redshift, pre-reionization universe is required for us to obtain the complete story of galaxy formation. This study is the first in a series that seeks to tell the story of galaxy formation at $z \gtrsim 10$ using a suite of large-scale adaptive mesh refinement cosmological simulations. Our machine-learning-accelerated surrogate model for Population III star formation and feedback, StarNet, gives us an unprecedented ability to obtain physically accurate, inhomogeneous chemical initial conditions for a statistically significant number of galaxies. We find that of the 12,423 halos in the mass range of $10^6\,\,M_\odot < M_\mathrm{vir} < 10^9\,\, M_\odot$ that form in our fiducial simulation, $16\%$ are chemically enriched by Population III supernovae by $z\sim12$. We then profile and compare various cooling processes at the centers of halos, and find a complete absence of atomic cooling halos. All of our halos with central cooling gas are dominated by H$_2$ cooling, metal cooling, or a mixture of the two, even in the presence of a strong H$_2$-photodissociating Lyman-Werner background. We also find that gas accretion through the virial radius is not driven by cooling. We find that gas virialization in halos with $M_\mathrm{vir}\gtrsim10^7\,\,M_\odot$ is supported by bulk turbulent flows, and that thermal energy accounts for only a small fraction of the total kinetic energy. Because of this, the mean gas temperature is well below the virial temperature for these halos. We then compute the mass of gas that is available for Population II star formation, and infer star formation rates for each potential star-forming halo.

The James Webb Space Telescope (JWST) has opened up a new window to study highly reddened explosive transients. We present results from late-time (1421 days post-explosion) JWST follow-up spectroscopic observations with NIRSpec and MIRI LRS of the intermediate luminosity red transient (ILRT) AT 2019abn located in the nearby Messier 51 galaxy (8.6 Mpc). ILRTs represent a mysterious class of transients which exhibit peak luminosities between those of classical novae and supernovae and which are known to be highly dust obscured. Similar to the prototypical examples of this class of objects, NGC 300 2008-OT and SN 2008S, AT 2019abn has an extremely red and dusty progenitor detected only in pre-explosion Spitzer/IRAC imaging at 3.6 and 4.5 micron and not in deep optical or near-infrared HST images. We find that late time observations of AT 2019abn from NEOWISE and JWST are consistent with the late time evolution of SN 2008S. In part because they are so obscured by dust, it is unknown what produces an ILRT with hypotheses ranging from high mass stellar merger events, non-terminal stellar outbursts, or terminal supernovae explosions through electron-capture in super-AGB stars. Our JWST observations show strong mid-IR Class C PAH features at 6.3 and 8.25 micron typical of carbon-rich post-AGB sources. These features suggest the dust around AT 2019abn, either pre-existing or newly formed in the ejecta, is composed of carbonaceous grains which are not typically observed around red supergiants. However, depending on the strength and temperature of hot bottom burning, SAGBs may be expected to exhibit a carbon-rich chemistry. Thus our JWST observations are consistent with AT 2019abn having an SAGB progenitor.

Non-Gaussian statistics of the projected weak lensing field are powerful estimators that can outperform the constraining power of the two-point functions in inferring cosmological parameters. This is because these estimators extract the non-Gaussian information contained in the small scales. However, fully leveraging the statistical precision of such estimators is hampered by theoretical uncertainties, such as those arising from baryonic physics. Moreover, as non-Gaussian estimators mix different scales, there exists no natural cut-off scale below which baryonic feedback can be completely removed. We therefore present a Bayesian solution for accounting for baryonic feedback uncertainty in weak lensing non-Gaussianity inference. Our solution implements Bayesian model averaging (BMA), a statistical framework that accounts for model uncertainty and combines the strengths of different models to produce more robust and reliable parameter inferences. We demonstrate the effectiveness of this approach in a Stage IV convergence peak counts analysis, including three baryonic feedback models. We find that the resulting BMA posterior distribution safeguards parameter inference against biases due to baryonic feedback, and therefore provides a robust framework for obtaining accurate cosmological constraints at Stage IV precision under model uncertainty scenarios.

Stellar tidal streams are a key tracer of galaxy evolution and have the potential to provide an indirect means for tracing dark matter. For the Local Group, many diffuse substructures have been identified and their link to galaxy evolution has been traced. However, an analysis of a larger sample is required to better probe the frequency and characteristics of these streams to verify the predictions of the Lambda-CDM model and its implementation in cosmological simulations. For that purpose, we are carrying out the first systematic survey of faint stellar debris from tidally disrupted dwarf satellites around nearby galaxies up to a distance of 100 Mpc. In this paper, we present a catalogue with the results of the first harvest of stellar tidal streams found by visual inspection in deep images of ~ 700 galaxies from the Dark Energy Survey (DES). We include a photometric characterisation of the streams obtained by measuring their surface brightnesses and colours. We found a total of 63 streams in our sample at distances between 40 and 100 Mpc, including 59 which were not previously reported. We measured their average surface brightness for the grz bands, to be 28.35+/-0.20, 27.81+/-0.13 and 27.62+/-0.09 mag arcsec^-2, respectively. By applying a statistical analysis to our findings, we obtained a stream detection frequency of 9.1% +/- 1.1% , in agreement with previous studies. We identified stream progenitors in 5-14% of our stream sample, depending on the confidence level. The first catalogue of streams in the Local Universe presented here will be complemented by future stream surveys within the Stellar Stream Legacy Survey. In this work we have learnt that the faintest measured stream surface brightness can be significantly brighter than the surface brightness limit of an image measured at pixel level, mainly due to correlated noise present in the images.

The HH 111 protostellar disk has recently been found to host a pair of spiral arms. Here we report the dust polarization results in the disk as well as the inner envelope around it, obtained with the Atacama Large Millimeter/submillimeter Array in continuum at lambda ~ 870 micron and ~ 0. 05" resolution. In the inner envelope, polarization is detected with a polarization degree of ~ 6% and an orientation almost everywhere parallel to the minor axis of the disk, and thus likely to be due to the dust grains magnetically aligned mainly by toroidal fields. In the disk, the polarization orientation is roughly azimuthal on the far side and becomes parallel to the minor axis on the near side, with a polarization gap in between on the far side near the central protostar. The disk polarization degree is ~ 2%. The polarized intensity is higher on the near side than the far side, showing a near-far side asymmetry. More importantly, the polarized intensity and thus polarization degree are lower in the spiral arms, but higher in between the arms, showing an anticorrelation of the polarized intensity with the spiral arms. Our modeling results indicate that this anticorrelation is useful for constraining the polarization mechanism and is consistent with the dust self-scattering by the grains that have grown to a size of ~ 150 micron. The interarms are sandwiched and illuminated by two brighter spiral arms and thus have higher polarized intensity. Our dust self-scattering model can also reproduce the observed polarization orientation parallel to the minor axis on the near side and the observed azimuthal polarization orientation at the two disk edges in the major axis.

Recent advancements have shown tensions between observations and our current understanding of the Universe. Such observations may include the Ho tension, as well as massive and mature galaxies at high redshifts that according to the existing models are older than the Universe itself. Since these observations are based on the redshift as the primary distance indicator, a bias in the redshift may explain these tensions. While the redshift follows an established model, when applied to astronomy it is based on the assumption that the rotational velocity of the Milky Way relative to the rotational velocity of the observed galaxies has a negligible effect on the redshift. But given the mysterious nature of the physics of galaxy rotation, that assumption should be tested. The test is done by comparing the redshift of galaxies that rotate in the same direction relative to the Milky Way to the redshift of galaxies that rotate in the opposite direction relative to the Milky Way. The results show that the mean redshift of galaxies that rotate in the same direction relative to the Milky Way is significantly higher than the mean redshift of galaxies that rotate in the opposite direction. Additionally, the redshift difference becomes larger as the redshift gets higher. The consistency of the analysis was verified by comparing data collected by three different telescopes, annotated using four different methods, released by three different research teams, and cover both the Northern and Southern ends of the galactic pole. All datasets are in excellent agreement with each other, showing consistency in the observed redshift bias. Given the "reproducibility crisis" in science, all datasets used in this study are publicly available, and the results can be easily reproduced. The observation could be a first direct observation for the Zwicky's "tired-light" model.

At metallicities lower than that of the Small Magellanic Cloud, it remains essentially unexplored how fossil magnetic fields, forming large-scale magnetospheres, could affect the evolution of massive stars, thereby impacting the fundamental building blocks of the early Universe. We extend our stellar evolution model grid with representative calculations of main-sequence, single-star models with initial masses of 20 and 60 M$_\odot$, including appropriate changes for low-metallicity environments ($Z = 10^{-3}-10^{-6}$). We scrutinise the magnetic, rotational, and chemical properties of the models. When lowering the metallicity, the rotational velocities can become higher and the tendency towards quasi-chemically homogeneous evolution increases. While magnetic fields aim to prevent the development of this evolutionary channel, the weakening stellar winds lead to less efficient magnetic braking in our models. Since the stellar radius is almost constant during a blueward evolution caused by efficient chemical mixing, the surface magnetic field strength remains unchanged in some models. We find core masses at the terminal-age main sequence between 22 and 52 M$_\odot$ for initially 60 M$_\odot$ models. This large difference is due to the vastly different chemical and rotational evolution. We conclude that in order to explain chemical species and, in particular, high nitrogen abundances in the early Universe, the adopted stellar models need to be under scrutiny. The assumptions regarding wind physics, chemical mixing, and magnetic fields will strongly impact the model predictions.

The stellar-mass black hole surrounded by the neutrino-dominated accretion flow (NDAF) is proposed to be the central engine of the gamma-ray burst (GRB). In this work, we investigate the neutrino/anti-neutrino luminosity and annihilation luminosity from the NDAF and pair-annihilation model, taking into account the neutrino oscillation above the accretion disk. The disk hydrodynamical properties are modelled using the empirical solution previously derived with boundary conditions including the effect of electron degeneracy and neutrino trapping. Our key parameters are the mass accretion rate and the black hole spin given in the range of $\dot{M} = 0.1$ - $10$ $M_\odot$ s$^{-1}$ and $ 0 \le a < 1$, respectively. Without neutrino oscillation, the obtained neutrino/anti-neutrino luminosity is $\sim 10^{51}$ - $10^{53}$ $\textrm{erg s}^{-1}$, while the neutrino annihilation luminosity is found to be $\sim 10^{46}$ - $10^{51}$ $\textrm{erg s}^{-1}$. In the presence of neutrino oscillation in the vacuum limit, the electron neutrino annihilation luminosity decreases by $\lesssim 22\%$ through the flavor transformation, while the muon- and tau-neutrino luminosity can increase up to $\sim 45\%$ and $ 60\%$, respectively. As a result, the total annihilation luminosity can be reduced up to $\sim 19\%$ due to the oscillation process above the disk. Finally, we also investigate the case whereby the CP-violating phase is changed from $\delta \textrm{CP} =$ $0^{\circ}$ to $245^{\circ}$. However, our results reveal that the CP-violating phase has minimal impact on neutrino annihilation luminosity.

Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $\sigma$ and 3\ $\sigma$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav.

We report the spectroscopic confirmation and fundamental properties of TOI-757 b, a mini-Neptune on a 17.5-day orbit transiting a bright star ($V = 9.7$ mag) discovered by the TESS mission. We acquired high-precision radial velocity measurements with the HARPS, ESPRESSO, and PFS spectrographs to confirm the planet detection and determine its mass. We also acquired space-borne transit photometry with the CHEOPS space telescope to place stronger constraints on the planet radius, supported with ground-based LCOGT photometry. WASP and KELT photometry were used to help constrain the stellar rotation period. We also determined the fundamental parameters of the host star. We find that TOI-757 b has a radius of $R_{\mathrm{p}} = 2.5 \pm 0.1 R_{\oplus}$ and a mass of $M_{\mathrm{p}} = 10.5^{+2.2}_{-2.1} M_{\oplus}$, implying a bulk density of $\rho_{\text{p}} = 3.6 \pm 0.8$ g cm$^{-3}$. Our internal composition modeling was unable to constrain the composition of TOI-757 b, highlighting the importance of atmospheric observations for the system. We also find the planet to be highly eccentric with $e$ = 0.39$^{+0.08}_{-0.07}$, making it one of the very few highly eccentric planets among precisely characterized mini-Neptunes. Based on comparisons to other similar eccentric systems, we find a likely scenario for TOI-757 b's formation to be high eccentricity migration due to a distant outer companion. We additionally propose the possibility of a more intrinsic explanation for the high eccentricity due to star-star interactions during the earlier epoch of the Galactic disk formation, given the low metallicity and older age of TOI-757.

We investigate the properties of AGB and post-AGB (PAGB) stars, planetary nebulae, and young stellar objects (YSOs) in our Galaxy through an analysis of observational data covering visual and infrared (IR) wavelengths. Utilizing datasets from IRAS, 2MASS, AllWISE, Gaia DR3, and the SIMBAD object database, we perform an in-depth comparison between observational data and theoretical models. For this comparison, we present various color-magnitude diagrams (CMDs) in visual and IR bands, as well as IR two-color diagrams (2CDs). Our results demonstrate that the CMDs, incorporating the latest distance and extinction data from Gaia DR3 for a majority of sample stars, are effective in distinguishing different classes of stars. To improve the precision of our analysis, we apply theoretical radiative transfer models for dust shells around AGB and PAGB stars. A thorough comparison of the theoretical models with observations across various IR 2CDs and CMDs shows a significant agreement. We find that AGB and PAGB stars are among the brightest classes in visual and IR bands. Furthermore, most YSOs are clearly distinguishable from AGB stars on various IR CMDs, exhibiting fainter absolute magnitudes in IR bands.

In the standard model of core-collapse supernova (CCSN), all neutrinos are assumed to be in pure flavor eigenstates in CCSN cores, but the assumption becomes invalid if neutrino distributions are unstable to flavor conversions. In this paper, we present a study of the occurrences of two representative neutrino-flavor instabilities, fast- and collisional flavor instabilities, in the cooling phase of proto-neutron star (PNS) from 1- to 50 seconds. We follow the long-term evolution of a PNS under spherically symmetric and quasi-static approximations, in which the matter profile is determined by solving the Tolman-Oppenheimer-Volkoff equation with neutrino feedback under the treatment of multi-group flux limited diffusion. For the stability analysis of neutrino flavor conversions, we recompute neutrino distributions using Monte Carlo transport in order to obtain the full angular distribution needed to compute the dispersion relations. We find no signs of flavor conversions in our models; the physical reason is thoroughly investigated. We also argue that the negative conclusion in flavor conversions could be changed qualitatively if multi-dimensional effects are included, as similar to cases in the earlier phase of CCSN.

We present the detection and characterization of fluctuations in linearly polarized emission from the atmosphere above the South Pole. These measurements make use of Austral winter survey data from the SPT-3G receiver on the South Pole Telescope in three frequency bands centered at 95, 150, and 220 GHz. We use the cross-correlation between detectors to produce an unbiased estimate of the power in Stokes I, Q, and U parameters on large angular scales. Our results are consistent with the polarized signal being produced by the combination of Rayleigh scattering of thermal radiation from the ground and thermal emission from a population of horizontally aligned ice crystals with an anisotropic distribution described by Kolmogorov turbulence. The signal is most significant at large angular scales, high observing frequency, and low elevation angle. Polarized atmospheric emission has the potential to significantly impact observations on the large angular scales being targeted by searches for inflationary B-mode CMB polarization. We present the distribution of measured angular power spectrum amplitudes in Stokes Q and I for 4 years of winter observations, which can be used to simulate the impact of atmospheric polarization and intensity fluctuations at the South Pole on a specified experiment and observation strategy. For the SPT-3G data, downweighting the small fraction of significantly contaminated observations is an effective mitigation strategy. In addition, we present a strategy for further improving sensitivity on large angular scales where maps made in the 220 GHz band are used to measure and subtract the polarized atmosphere signal from the 150 GHz band maps. In observations with the SPT-3G instrument at the South Pole, the polarized atmospheric signal is a well-understood and sub-dominant contribution to the measured noise after implementing the mitigation strategies described here.

As a byproduct of our search for Galactic stellar systems with gamma-ray emission, we have identified an unrelated cool and evolved star (IRC-10412) that attracted our attention due to its strong radio emission level with a spectral index matching, almost perfectly, the canonical +0.6 value expected from an ionized stellar wind. A follow-up observational analysis was undertaken given that these two properties are hard to reconcile as originating in the same stellar object. As a result, IRC-10412 has been classified as a new semiregular variable of SRb type in the asymptotic giant branch, and different but consistent estimates of its mass-loss parameter are reported. We propose that its unusually high radio emission arises from a ~ 1E-5 Msun/yr stellar wind exposed to an external source of ionizing photons, possibly coming from nearby OB associations.

In this manuscript, strong clues are reported to support the normal broad line AGN SDSS J1617+0638 as the host galaxy harboring a central tidal disruption event (TDE). Through the optical flare in the CSS 8.5years-long light curve and the none-variability in the up-to-date ASAS-SN light curves, the theoretical TDE model described by the MOSFIT code can be applied in SDSS J1617+0638. Meanwhile, considering the assumed central TDE expected continuum emissions not strong enough to describe the continuum emissions in the SDSS spectrum of SDSS J1617+0638, an additional power law component from pre-existing AGN activity should be necessary in SDSS J1617+0638. Furthermore, considering the short time duration to the observed date for the SDSS spectrum from the starting time of the assumed central TDE in SDSS J1617+0638, TDE model expected accreting mass only about 0.03$M_\odot$ can lead to few effects of TDEs debris on the observed broad emission lines in the SDSS spectrum of SDSS J1617+0638, indicating the TDE model determined BH mass simply consistent with the virial BH mass by broad emission lines, as determined results in SDSS J1617+0638. Therefore, through both the photometric variability and the spectroscopic results, a central TDE can be preferred in the normal broad line AGN SDSS J1617+0638 with pre-existing central AGN activity and pre-existing broad emission line regions.

In this work the random forest algorithm GALPRO is implemented to generate photometric redshift posteriors, and its performance when trained and then applied to data from another survey is investigated. The algorithm is initially calibrated using a truth dataset compiled from the DESI Legacy survey. We find that the testing and training datasets must have very similar redshift distributions, with the range of their photometric data overlapping by at least 90% in the appropriate photometric bands in order for the training data to be applicable to the testing data. Then GALPRO is again trained using the DESI dataset and then applied to a sample drawn from the PanSTARRS survey, to explore whether GALPRO can be first trained using a trusted dataset and then applied to an entirely new survey, albeit one that uses a different magnitude system for its photometric bands, thus requiring careful conversion of the measured magnitudes for the new survey before GALPRO can be applied. The results of this further test indicate that GALPRO does not produce accurate photometric redshift posteriors for the new survey, even where the distribution of redshifts for the two datasets overlaps by over 90%. Hence, we conclude that the photometric redshifts generated by GALPRO are not suitable for generating estimates of photometric redshifts and their posterior distribution functions when applied to an entirely new survey, particularly one that uses a different magnitude system. However, our results demonstrate that GALPRO is a useful tool for inferring photometric redshift estimates in the case where a spectroscopic galaxy survey is nearly complete, but is missing some spectroscopic redshift values.

Additive manufacture (AM) involves creating a part layer by layer and is a rapidly evolving manufacturing process. It has multiple strengths that apply to space-based optics, such as the ability to consolidate multiple parts into one, reducing the number of interfaces. The process also allows for greater mass reduction, making parts more cost-effective to launch, achieved by optimising the shape for intended use or creating intricate geometries like lattices. However, previous studies have highlighted issues associated with the AM process. For example, when trying to achieve high-precision optical surfaces on AM parts, the latticing on the underside of mirrors can provide insufficient support during machining, resulting in the quilting effect. This paper builds on previous work and explores such challenges further. This will be implemented by investigating ways to apply AM to a deployable mirror from a CubeSat project called A-DOT. The reflective surface has a spherical radius of curvature of 682 mm and approximate external dimensions of 106 mm x 83 mm. The aim is to produce two mirrors that will take full advantage of AM design benefits and account for the challenges in printing and machining a near-net shape. The designs will have reduced mass by using selected internal lattice designs and topology-optimised connection points, resulting in two mirrors with mass reduction targets of 50% and 70%. Once printed in aluminium using laser powder bed fusion, the reflective surface will be created using single point diamond turning. Finally, an evaluation of the dimensional accuracy will be conducted, using interferometry, to quantify the performance of the reflective surface.

The Southern Photometric Local Universe Survey (S-PLUS) is a project to map $\sim9300$ sq deg of the sky using twelve bands (seven narrow and five broadbands). Observations are performed with the T80-South telescope, a robotic telescope located at the Cerro Tololo Observatory in Chile. The survey footprint consists of several large contiguous areas, including fields at high and low galactic latitudes, and towards the Magellanic Clouds. S-PLUS uses fixed exposure times to reach point source depths of about $21$ mag in the $griz$ and $20$ mag in the $u$ and the narrow filters. This paper describes the S-PLUS Data Release 4 (DR4), which includes calibrated images and derived catalogues for over 3000 sq deg, covering the aforementioned area. The catalogues provide multi-band photometry performed with the tools \texttt{DoPHOT} and \texttt{SExtractor} -- point spread function (\PSF) and aperture photometry, respectively. In addition to the characterization, we also present the scientific potential of the data. We use statistical tools to present and compare the photometry obtained through different methods. Overall we find good agreement between the different methods, with a slight systematic offset of 0.05\,mag between our \PSF and aperture photometry. We show that the astrometry accuracy is equivalent to that obtained in previous S-PLUS data releases, even in very crowded fields where photometric extraction is challenging. The depths of main survey (MS) photometry for a minimum signal-to-noise ratio $S/N = 3$ reach from $\sim19.5$ for the bluer bands to $\sim21.5$ mag on the red. The range of magnitudes over which accurate \PSF photometry is obtained is shallower, reaching $\sim19$ to $\sim20.5$ mag depending on the filter. Based on these photometric data, we provide star-galaxy-quasar classification and photometric redshift for millions of objects.

There is increasing evidence that massive stars may exhibit an enhanced mass loss shortly before their termination explosion. Some of them also indicate the enhancement of their circumstellar matter (CSM) is not spherically symmetric. Supernova (SN) interacting with aspherical CSM could induce special polarization signals from multiple radiation components that deviate from spherical symmetry. We investigate the time-evolution of the continuum polarization induced by the SN ejecta interacting with a disk/torus-like CSM. Our calculation suggests that the interaction between the SN ejecta and an immediate disk-like CSM with a thin, homogenous density structure would produce a high continuum polarization, which may reach a peak level of $\sim$15\%. The interplay between the evolving geometry of the emitting regions and the time-variant flux ratio between the polar ejecta and the equatorial CSM interaction may produce a double-peaked feature in the polarization time sequence. A similar trend of the time evolution of the polarization is also found for a radially extended CSM disk that exhibits a wind-like density structure, with an overall relatively lower level of continuum polarization ($<2.5\%$) during the interaction process. We also identify a non-uniform temperature distribution along the radial direction of the CSM disk, which yields a strong wavelength dependence of the continuum polarization. These signatures provide a unique geometric diagnostic to explore the interaction process and the associated extreme mass loss of the progenitors of interacting transients.

We present a new supergiant-rich stellar cluster hidden by extinction and christen it as Barbá 2, in honor of its discoverer Rodolfo Barbá. The cluster is at a distance of $7.39^{+0.65}_{-0.55}$ kpc and contains several supergiants, of which we provide spectral classifications for one blue, one yellow, and five red ones. The cluster extinction indicates an above-average grain-size ($R_{5495} \sim 3.7$), its age has a minimum value of 10 Ma, and its core radius is $0.84\pm0.19$ pc.

The \textit{Vera C. Rubin} Observatory will start operations in 2025. During the first two years, too few visits per target per band will be available, meaning that mean magnitude measurements of variable stars will not be precise and thus, standard candles like RR Lyrae (RRL) will not be usable. Light curve templates (LCTs) can be adopted to estimate the mean magnitude of a variable star with few magnitude measurements, provided that their period (plus amplitude and reference epoch, depending on how the LCT is applied) is known. LSST will provide precise RRL periods within the first six months, allowing to exploit RRLs if LCTs were available. We aim to build LCTs in the LSST bands to enhance the early science with LSST. Using them will provide a 1-2 years advantage with respect to a classical approach, concerning distance measurements. We collected $gri$-band data from the ZTF survey and $z$-band data from DECam to build the LCTs of RRLs. We also adopted synthetic $griz$-band data in the LSST system from pulsation models, plus SDSS, \gaia and OGLE photometry, inspecting the light amplitude ratios in different photometric systems to provide useful conversions to apply the LCTs. We have built LCTs of RRLs in the $griz$ bands of the LSST photometric system; for the $z$ band, we could build only fundamental-mode RRL LCTs. We quantitatively demonstrated that LCTs built with ZTF and DECam data can be adopted on the LSST photometric system. LCTs will decrease by a factor of at least two the uncertainty on distance estimates of RRLs, with respect to a simple average of the available measurements. Finally, within our tests, we have found a brand new behavior of amplitude ratios in the Large Magellanic Cloud.

The 23 Myr young star $\beta$ Pictoris is a laboratory for planet formation studies given its observed debris disk, its two directly imaged super-Jovian planets, and the evidence of transiting extrasolar comets. The most recent evidence of exocometary transits around $\beta$ Pic came from the TESS space mission. Previous analyses of these transits constrained the orbital distribution of the underlying exocomet population to range between about 0.03 and 1.3 AU assuming a fixed transit impact parameter. Here we examine the distribution of the observed transit durations (${\Delta}t$) to infer the orbital surface density distribution ($\delta$) of the underlying exocomet sample. We show that a narrow belt of exocomets around $\beta$ Pic, in which the transit impact parameters are randomized but the orbital semimajor axes are equal, results in a pile-up of long transit durations. This is contrary to observations, which reveal a pile-up of short transit durations (${\Delta}t \approx 0.1$ d) and a tail of only a few transits with ${\Delta}t > 0.4$ d. A flat density distribution of exocomets between about 0.03 and 2.5 AU results in a better match between the resulting ${\Delta}t$ distribution and the observations but the slope of the predicted ${\Delta}t$ histogram is not sufficiently steep. An even better match to the observations can be produced with a $\delta \propto a^{\beta}$ power law. Our modeling reveals a best fit between the observed and predicted ${\Delta}t$ distribution for $\beta = -0.15_{-0.10}^{+0.05}$. A more reasonable scenario in which the exocometary trajectories are modeled as hyperbolic orbits can also reproduce the observed ${\Delta}t$ distribution to some extent. Our results imply that cometary material exists on highly eccentric orbits with a more extended range of semimajor axes than suggested by previous spectroscopic observations.

The formation of galaxies is significantly influenced by galactic winds, possibly driven by cosmic rays due to their long cooling times and better coupling to plasma compared to radiation. In this study, we compare the radio observations of the edge-on galaxy NGC 4217 from the CHANG-ES collaboration catalog with a mock observation of an isolated galaxy based on the arepo simulation that adopts the state-of-the-art two-moment cosmic ray transport treatment and multiphase interstellar medium model. We find significant agreement between the simulated and observed images and spectroscopic data. Specifically, we find that (i) while the radio flux is sensitive to the magnitude of the magnetic field, the shapes of intensity profiles are only weakly dependent on the normalization of the magnetic field and are robustly predicted by our galaxy simulations; (ii) the multi-wavelength spectrum above 0.1~GHz is in agreement with the radio observations and its slope is also only weakly sensitive to the magnetic field strength; (iii) the magnetic field direction exhibits an X-shaped morphology, often seen in edge-on galaxies, which is consistent with the observations and indicates the presence of a galactic-scale outflow. Since the normalization of the intensity profiles and spectrum are degenerate with respect to the magnetic field strength and the normalization of the electron spectrum, this qualitative match is possible provided that the simulated magnetic field is only mildly boosted, which could be attributed to the missing galactic disk interaction with the CGM in our idealized simulation. Our results highlight the importance of incorporating advanced cosmic ray transport models in simulations and provide a deeper understanding of galactic wind dynamics and its impact on galaxy evolution.

Precise physical properties of the known transiting exoplanets are essential for their precise atmospheric characterization using modern and upcoming instruments. Leveraging the large volume of high SNR photometric follow-up data from TESS, highly precise physical properties can be estimated for these systems, especially for those discovered using ground-based instruments prior to the TESS mission. In this work, I have used the publicly available TESS follow-up data for 28 transiting systems with 10 $<$ V$_{mag}$ $<$ 10.5, with an aim to update their known physical properties. The observed lightcurves have been analysed by implementing a state-of-the-art critical noise treatment algorithm to effectively reduce both time-correlated and un-correlated noise components, using sophisticated techniques like wavelet denoising and Gaussian-process regression. Compared with the previous studies, the estimated transit parameters are found to be more precise for most of the targets, including a few cases where a larger space-based instrument like Spitzer, Kepler or CHEOPS has been used in the previous study. The large volume of transit observations used for each target has also resulted in a more accurate estimation of the physical properties, as this overcomes any error in parameter estimations from bias present in a smaller volume of data. Thus, comparing with the literature values, statistically significant improvements in the known physical properties of several targeted systems have been reported from this work. The large volume of transit timing information from the analyses was also used to search for Transit Timing Variation trends in these targets, which has resulted in no significant detection.

Understanding star formation in galaxies requires resolving the physical scale on which star formation often occurs: the scale of star clusters. We present a multiwavelength, eight-parsec resolution study of star formation in the circumnuclear star cluster and molecular gas rings of the early-type spiral NGC 1386. The cluster ring formed simultaneously ~ 4 Myr ago. The clusters have similar properties in terms of mass and star formation rate, resembling those of H II regions in the Milky Way disc. The molecular CO gas resolves into long filaments, which define a secondary ring detached from the cluster ring. Most clusters are in CO voids. Their separation with respect the CO filaments is reminiscent of that seen in galaxy spiral arms. By analogy, we propose that a density wave through the disc of this galaxy may have produced this gap in the central kpc. The CO filaments fragment into strings of dense, unresolved clouds with no evidence of a stellar counterpart. These clouds may be the sites of a future population of clusters in the ring. The free-fall time of these clouds, ~ 10 Myr, is close to the orbital time of the CO ring. This coincidence could lead to a synchronous bursting ring, as is the case for the current ring. The inward spiralling morphology of the CO filaments and co-spatiality with equivalent kpc-scale dust filaments are suggestive of their role as matter carriers from the galaxy outskirts to feed the molecular ring and a moderate active nucleus.

In the past, high-z AGNs were given a minor role as possible drivers of reionization, despite initial evidences in favor of their large space densities at low luminosities by Chandra and HST. Recent observations from JWST are finding relatively large numbers of faint AGNs at z>4, convincingly confirming these early results. We present a sample of z~5 AGNs (both from wide, shallow ground-based surveys and from deep, pencil-beam observations from JWST), allowing to estimate their space densities with unprecedented accuracy. The bright end (M1450<-26) of the z~5 AGN luminosity function is well constrained, with a rather steep slope. The faint end (M1450>-22) indicates a high space density, the scatter is significant and the knee (M1450~-24) is mostly undetermined. Comparisons with state-of-the-art models find reasonable agreement with the observed AGN luminosity function at z=5, while the predicted space density evolution at higher redshifts appears to be too fast with respect to observational constraints. Given the large variance at the faint end, we consider different options in fitting the luminosity functions and deriving the ionizing emissivity. Even in the most conservative scenario, the photo-ionization rate produced by z~5 AGNs is consistent with the UV background measurements. A slow evolution of the space density of faint AGNs is observed, indicating that active SMBHs are probably producing large amounts of ionizing photons at z>6, well into the epoch of reionization. This is an important indication that high-z AGNs could be the major contributors to the reionization of the Universe.

In our Solar system, spin-orbit coupling is a common phenomenon in binary asteroid systems, where the mutual orbits are no longer invariant due to exchange of angular momentum between translation and rotation. In this work, dynamical structures in phase space are explored for the problem of spin-orbit coupling by taking advantage of analytical and numerical methods. In particular, the technique of Poincaré sections is adopted to reveal numerical structures, which are dependent on the total angular momentum, the Hamiltonian, mass ratio and asphericity parameter. Analytical study based on perturbative treatments shows that high-order and/or secondary spin-orbit resonances are responsible for numerical structures arising in Poincaré sections. Analytical solutions are applied to (65803) Didymos, (80218) ${\rm VO}_{123}$ and (4383) Suruga to reveal their phase-space structures, showing that there is a high possibility for them to locate inside secondary 1:1 spin-orbit resonance.

We explore the science prospects of a 14,000 deg$^2$ Kinematic Lensing (KL) survey with the Dark Energy Spectroscopic Instrument (DESI) and overlapping imaging surveys. KL infers the cosmic shear signal by jointly forward modeling the observed photometric image and velocity field of a disk galaxy. The latter can be constrained by placing multiple DESI fibers along the galaxy's major and minor axis, a concept similar to the DESI Peculiar Velocity Survey. We study multiple subset galaxy samples of the DESI Legacy Survey Data Release 9 catalog and quantify the residual shape noise, $\sigma_\epsilon$, of each sample as a function of cuts in $r$-band magnitude using mock observations. We conduct simulated likelihood analyses for these galaxy samples and find that a DESI-KL program can place highly interesting constraints on structure formation at very low redshifts, i.e. $\sigma_8(z<0.15)$. We conclude that if the $S_8$ tension consolidates as a phenomenon, a DESI-KL survey can provide unique insights into this phenomenon in the very late-time Universe. Given the different footprints of DESI and Rubin Observatory's Legacy Survey of Space and Time (LSST), lensing results from both surveys are highly complementary and can be combined into a joint lensing survey. We further note that DESI-KL benefits multiple additional science cases, e.g. studies of modified gravity models when combined with peculiar velocity surveys, and dark matter studies that are based on galaxy-galaxy lensing of dwarf galaxies.

Prime-Cam is a first-generation science instrument for the CCAT Observatory's six-meter aperture Fred Young Submillimeter Telescope (FYST). FYST's crossed-Dragone design provides high optical throughput to take advantage of its unique site at 5600 m on Cerro Chajnantor in Chile's Atacama Desert to reach mapping speeds over ten times greater than current and near-term submillimeter experiments. Housing up to seven independent instrument modules in its 1.8-meter diameter cryostat, Prime-Cam will combine broadband polarization-sensitive modules and spectrometer modules designed for observations in several frequency windows between 210 GHz and 850 GHz to study a wide range of astrophysical questions from Big Bang cosmology to the formation of stars and galaxies in the Epoch of Reionization and beyond. In order to cover this range of frequencies and observation modes, each of the modules contains a set of cold reimaging optics that is optimized for the science goals of that module. These optical setups include several filters, three or four anti-reflection-coated silicon lenses, and a Lyot stop to control the field of view and illumination of the primary mirror, satisfy a series of mechanical constraints, and maximize optical performance within each passband. We summarize the design considerations and trade-offs for the optics in these modules and provide a status update on the fabrication of the Prime-Cam receiver and the design of its 1 K and 100 mK thermal BUSs.

We report observations of deuterated acetylene (C$_2$HD) at 19.3 $\mu$m (519 cm$^{-1}$) with the Texas Echelon Cross Echelle Spectrograph on the NASA Infrared Telescope Facility in July 2017. Six individual lines from the Q-branch of the $\nu_4$ band were clearly detected with a S/N ratio up to 10. Spectral intervals around 8.0 $\mu$m (745 cm$^{-1}$) and 13.4 $\mu$m (1247 cm$^{-1}$) containing acetylene (C$_2$H$_2$) and methane (CH$_4$) lines respectively, were observed during the same run to constrain the disk-averaged C$_2$H$_2$ abundance profile and temperature profile. Cassini observations with the Composite Infrared Spectrometer (CIRS) were used to improve the flux calibration and help to constrain the atmospheric model. The measured D/H ratio in acetylene, derived from the C$_2$HD/C$_2$H$_2$ abundance ratio, is (1.22$^{+0.27}_{-0.21})$ $\times$ 10$^{-4}$, consistent with that in methane obtained in previous studies. Possible sources of fractionation at different steps of the acetylene photochemistry are investigated.

We present the results of the ALMA Large Program Early Planet Formation in Embedded disks observations of the Class 0 protostar GSS30 IRS3. Our observations included 1.3 mm continuum with a resolution of 0.''05 (7.8 au) and several molecular species including $^{12}$CO, $^{13}$CO, C$^{18}$O, H$_{2}$CO and c-C$_{3}$H$_{2}$. The dust continuum analysis unveiled a disk-shaped structure with a major axis size of $\sim$200 au. We observed an asymmetry in the minor axis of the continuum emission suggesting that the emission is optically thick and the disk is flared. On the other hand, we identified two prominent bumps along the major axis located at distances of 26 and 50 au from the central protostar. The origin of the bumps remains uncertain and might be due to an embedded substructure within the disk or the result of the temperature distribution instead of surface density due to optically thick continuum emission. The $^{12}$CO emission reveals a molecular outflow consisting of three distinct components: a collimated one, an intermediate velocity component exhibiting an hourglass shape, and a wider angle low-velocity component. We associate these components with the coexistence of a jet and a disk-wind. The C$^{18}$O emission traces both a Keplerian rotating circumstellar disk and the infall of the rotating envelope. We measured a stellar dynamical mass of 0.35$\pm$0.09 M$_{\odot}$.

(abridged) We present an account of known cataclysmic variables (CVs) that were detected as X-ray sources in eROSITA X-ray surveys and have Gaia DR3 counterparts. We address standard CVs with main sequence donors and white dwarfs accreting via Roche-lobe overflow (RLOF) and related objects, the double degenerates (DDs), and the symbiotic stars (SySts). We discern between nonmagnetic (dwarf novae and nova-like objects) and magnetic CVs (polars and intermediate polars (IPs)). In the publically available eROSITA catalog from the recent DR1, typically 65\% of known cataloged and classified CVs are detected. This fraction rises to over 90\% if the stack of all \ero X-ray surveys (called S45 in this paper) is considered and the search volume is restricted to a radius of 500\,pc. We examine the various classes of CVs in various diagnostic diagrams relating X-ray and optical properties (luminosity, absolute magnitude, color, X-ray spectral hardness, and optical variability) and establish their average class properties. We derive spectral properties for the 22 brightest polars and confirm an increase in the ratio of soft to hard X-rays with increasing magnetic field in the accretion region. We report three new soft IPs and present a spectral analysis of all soft IPs. Their blackbody temperatures agree well with published values. The DDs represent the bluest and faintest subcategory but reach the same identification fraction as the standard CVs. The SySts are the most distant systems; only 20 (13\%) were detected as X-ray sources in S45, and 7 of those are first-time detections. We investigate their mean properties using an upper limit on the flux of the nondetected CVs. Their X-ray nondetection is indeed a distance effect.

The 21 cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation (EoR). This has led to the construction of low-frequency radio interferometric arrays, such as the Hydrogen Epoch of Reionization Array (HERA), aimed at systematically mapping this emission for the first time. Precision calibration, however, is a requirement in 21 cm radio observations. Due to the spatial compactness of HERA, the array is prone to the effects of mutual coupling, which inevitably lead to non-smooth calibration errors that contaminate the data. When unsmooth gains are used in calibration, intrinsically spectrally-smooth foreground emission begins to contaminate the data in a way that can prohibit a clean detection of the cosmological EoR signal. In this paper, we show that the effects of mutual coupling on calibration quality can be reduced by applying custom time-domain filters to the data prior to calibration. We find that more robust calibration solutions are derived when filtering in this way, which reduces the observed foreground power leakage. Specifically, we find a reduction of foreground power leakage by 2 orders of magnitude at k=0.5.

The Mg II k resonance line is commonly used for diagnosing the solar chromosphere. Here we theoretically investigate its intensity and polarization in solar prominences, taking into account the effects of 3D radiative transfer and the Hanle and Zeeman effects. We use an optically thick 3D model representative of a solar prominence and apply several inversion methods to the synthetic Stokes profiles, clarifying their pros and cons for inferring prominence magnetic fields. We conclude that the self-consistent 3D inversion with radiative transfer is necessary to determine the magnetic field vector, although its geometry cannot be inferred with full fidelity. We also demonstrate that more traditional methods, such as those based on the weak field approximation or the constant-property slab assumption, can offer useful information under certain conditions.

Over the last decade, studies of large samples of binary systems identified chemical anomalies, and showed that they might be attributed to planet formation or planet engulfment. However, both scenarios have primarily been tested in pairs without known exoplanets. In this work, we explore these scenarios in the newly detected planet-hosting wide binary TOI-1173 A/B (projected separation $\sim 11,400$ AU) using high-resolution MAROON-X and ARCES spectra. We determined photospheric stellar parameters both by fitting stellar models and via the spectroscopic equilibrium approach. Both analyses agree and suggest that they are cool main sequence stars located in the thin disc. A line-by-line differential analysis between the components (B$-$A) displays an abundance pattern in the condensation temperature plane where the planet-hosting star TOI-1173 A is enhanced in refractory elements such as iron by more than 0.05 dex. This suggests the engulfment of $\sim$18 M$_{\oplus}$ of rocky material in star A. Our hypothesis is supported by the dynamics of the system detailed in our companion paper Yana Galarza et al. 2024, which suggest that the Super-Neptune TOI-1173 A $b$ might have been delivered to its current short period ($\sim7$ days) through circulatization and von Zeipel-Lidov-Kozai mechanisms, thereby triggering the engulfment of inner rocky exoplanets.

Young, self-luminous super-Jovian companions discovered by direct imaging provide a challenging test of planet formation and evolution theories. By spectroscopically characterizing the atmospheric compositions of these super-Jupiters, we can constrain their formation histories. Here we present studies of the recently discovered HIP 99770 b, a 16 MJup high-contrast companion on a 17 au orbit, using the fiber-fed high-resolution spectrograph KPIC (R~35,000) on the Keck II telescope. Our K-band observations led to detections of H2O and CO in the atmosphere of HIP 99770 b. We carried out free retrieval analyses using petitRADTRANS to measure its chemical abundances, including the metallicity and C/O ratio, projected rotation velocity (vsini), and radial velocity (RV). We found that the companion's atmosphere has C/O=0.55(-0.04/+0.06) and [M/H]=0.26(-0.23/+0.24) (1{\sigma} confidence intervals), values consistent with those of the Sun and with a companion formation via gravitational instability or core accretion. The projected rotation velocity < 7.8 km/s is small relative to other directly imaged companions with similar masses and ages. This may imply a near pole-on orientation or effective magnetic braking by a circumplanetary disk. In addition, we added the companion-to-primary relative RV measurement to the orbital fitting and obtained updated constraints on orbital parameters. Detailed characterization of super-Jovian companions within 20 au like HIP 99770 b is critical for understanding the formation histories of this population.

Recent neutron star surface observations corroborate a long-standing theory that neutron stars may be accelerated over extended periods after their birth. We analyze how these prolonged rocket-like accelerations, combined with rapid birth kicks, impact binary orbits. We find that even a small contribution of rocket kicks combined with instantaneous natal kicks can allow binaries to reach period--eccentricity combinations unattainable in standard binary evolution models. We propose these kick + rocket combinations as a new channel to form wide low-eccentricity neutron star binaries such as Gaia NS1, as well as inducing stellar mergers months to years after a supernova to cause peculiar high-energy transients.

We present our current update on the gas-phase chemical reactions and spectral lines of TiO in the spectral synthesis code Cloudy. For this purpose, we have added 229 Ti-related reactions in the chemical network. In addition, we consider 230 fine-structure energy levels, the corresponding 223 radiative transitions, and 444 collisional transitions with ortho and para H$_2$ and predict 66 TiO lines. We perform spectroscopic simulations of TiO emission from the circumstellar region of the oxygen-rich red supergiant VY Canis Majoris to validate our update. Our model reproduces the observed TiO column density. This update is helpful in modeling dust-free astrophysical environments where Ti is in the gas phase and TiO can form.

Using data from the first flight of SPIDER and from Planck HFI, we probe the properties of polarized emission from interstellar dust in the SPIDER observing region. Component separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses spanning the full SPIDER region demonstrate that i) the spectral energy distribution of diffuse Galactic dust emission is broadly consistent with a modified-blackbody (MBB) model with a spectral index of $\beta_\mathrm{d}=1.45\pm0.05$ $(1.47\pm0.06)$ for $E$ ($B$)-mode polarization, slightly lower than that reported by Planck for the full sky; ii) its angular power spectrum is broadly consistent with a power law; and iii) there is no significant detection of line-of-sight decorrelation of the astrophysical polarization. The size of the SPIDER region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature $T_\mathrm{d}=19.6$ K, an analysis of two independent sub-regions of that field results in inferred values of $\beta_\mathrm{d}=1.52\pm0.06$ and $\beta_\mathrm{d}=1.09\pm0.09$, which are inconsistent at the $3.9\,\sigma$ level. Furthermore, a joint analysis of SPIDER and Planck 217 and 353 GHz data within a subset of the SPIDER region is inconsistent with a simple MBB at more than $3\,\sigma$, assuming a common morphology of polarized dust emission over the full range of frequencies. These modeling uncertainties have a small--but non-negligible--impact on limits on the cosmological tensor-to-scalar ratio derived from the \spider dataset. The fidelity of the component separation approaches of future CMB polarization experiments may thus have a significant impact on their constraining power.

The lack of benchmark data on the real minerals, native ammonium carriers in Solar System gives rise to controversial opinions on extraterrestrial ammonium reservoirs. We herein report on discovery of the first mineral carrier of meteoritic ammonium and show its relevance to the compositional and spectral characteristics of cometary and asteroidal bodies. Chemically distant from previously inferred volatile organics or ammoniated phyllosilicates, it is an aqueous metal-ammonium sulfate related to a family of so-called Tutton salts. Nickeloan boussingaultite, (NH4)2(Mg,Ni)(SO4)2 6H2O, occurs in Orgueil, a primitive carbonaceous chondrite closely related to (162173) Ryugu and (101955) Bennu, the C-type asteroids. The available spectroscopic, chemical and mineralogical data signify that natural Tutton salts perfectly fit into the role of ammonium reservoir under conditions of cometary nuclei and carbonaceous asteroids.

We present a formalism to extract the angular power spectrum of fields sampled at a finite number of points with arbitrary positions -- a common situation for several catalog-based astrophysical probes -- through a simple extension of the standard pseudo-$C_\ell$ algorithm. A key complication in this case is the need to handle the shot noise component of the associated discrete angular mask which, for sparse catalogs, can lead to strong coupling between very different angular scales. We show that this problem can be solved easily by estimating this contribution analytically and subtracting it. The resulting estimator is immune to small-scale pixelization effects and aliasing, and, more interestingly, unbiased against the contribution from measurement noise uncorrelated between different sources. We demonstrate the validity of the method in the context of cosmic shear datasets, and showcase its usage in the case of other spin-0 and spin-1 astrophysical fields of interest. We incorporate the method in the public $\texttt{NaMaster}$ code (this https URL).

The 3C 220.3 system is a rare case of a foreground narrow-line radio galaxy ("galaxy A," $z_A = 0.6850$) lensing a background submillimeter galaxy ($z_{\rm SMG1} = 2.221$). New spectra from MMT/Binospec confirm that the companion galaxy ("galaxy B") is part of the lensing system with $z_B = 0.6835$. New three-color HST data reveal a full Einstein ring and allow a more precise lens model. The new HST images also reveal extended emission around galaxy A, and the spectra show extended [OII] emission with irregular morphology and complex velocity structure. All indications are that the two lensing galaxies are a gravitationally interacting pair. Strong [OII] emission from both galaxies A and B suggests current star formation, which could be a consequence of the interaction. This would indicate a younger stellar population than previously assumed and imply smaller stellar masses for the same luminosity. The improved lens model and expanded spectral energy distributions have enabled better stellar-mass estimates for the foreground galaxies. The resulting dark matter fractions are ~0.8, higher than previously calculated. Deeper Chandra imaging shows extended X-ray emission but no evidence for a point X-ray source associated with either galaxy. The detection of X-rays from the radio lobes of 3C 220.3 allows an estimate of ~3 nT for the magnetic fields in the lobes, a factor of ~3 below the equipartition fields, as typical for radio galaxies.