Papers for Tuesday, Sep 12 2023

As the classic viscous paradigm for protoplanetary disk accretion is challenged by the observational evidence of low turbulence, the alternative scenario of MHD disk winds is being explored as potentially able to reproduce the same observed features traditionally explained with viscosity. Although the two models lead to different disk properties, none of them has been ruled out by observations - mainly due to instrumental limitations. In this work, we present a viable method to distinguish between the viscous and MHD framework based on the different evolution of the distribution in the disk mass ($M_{\mathrm{d}}$) - accretion rate ($\dot M$) plane of a disk population. With a synergy of analytical calculations and 1D numerical simulations, performed with the population synthesis code \texttt{Diskpop}, we find that both mechanisms predict the spread of the observed ratio $M_{\mathrm{d}}/\dot M$ in a disk population to decrease over time; however, this effect is much less pronounced in MHD-dominated populations as compared to purely viscous populations. Furthermore, we demonstrate that this difference is detectable with the current observational facilities: we show that convolving the intrinsic spread with the observational uncertainties does not affect our result, as the observed spread in the MHD case remains significantly larger than in the viscous scenario. While the most recent data available show a better agreement with the wind model, ongoing and future efforts to obtain direct gas mass measurements with ALMA and ngVLA will cause a reassessment of this comparison in the near future.

Long-term magma ocean phases on rocky exoplanets orbiting closer to their star than the runaway greenhouse threshold - the inner edge of the classical habitable zone - may offer insights into the physical and chemical processes that distinguish potentially habitable worlds from others. Thermal stratification of runaway planets is expected to significantly inflate their atmospheres, potentially providing observational access to the runaway greenhouse transition in the form of a "habitable zone inner edge discontinuity" in radius-density space. Here, we use Bioverse, a statistical framework combining contextual information from the overall planet population with a survey simulator, to assess the ability of ground- and space-based telescopes to test this hypothesis. We find that the demographic imprint of the runaway greenhouse transition is likely detectable with high-precision transit photometry for sample sizes $\gtrsim 100$ planets if at least ~10 % of those orbiting closer than the habitable zone inner edge harbor runaway climates. Our survey simulations suggest that in the near future, ESA's PLATO mission will be the most promising survey to probe the habitable zone inner edge discontinuity. We determine survey strategies that maximize the diagnostic power of the obtained data and identify as key mission design drivers: 1. A follow-up campaign of planetary mass measurements and 2. The fraction of low-mass stars in the target sample. Observational constraints on the runaway greenhouse transition will provide crucial insights into the distribution of atmospheric volatiles among rocky exoplanets, which may help to identify the nearest potentially habitable worlds.

The dominant form of oxygen in cold molecular clouds is gas-phase carbon monoxide (CO) and ice-phase water (H$_2$O). Yet, in planet-forming disks around young stars, gas-phase CO and H$_2$O are less abundant relative to their ISM values, and no other major oxygen-carrying molecules have been detected. Some astrochemical models predict that gas-phase molecular oxygen (O$_2$) should be a major carrier of volatile oxygen in disks. We report a deep search for emission from the isotopologue $^{16}$O$^{18}$O ($N_J=2_1-0_1$ line at 233.946 GHz) in the nearby protoplanetary disk around TW Hya. We used imaging techniques and matched filtering to search for weak emission but do not detect $^{16}$O$^{18}$O. Based on our results, we calculate upper limits on the gas-phase O$_2$ abundance in TW Hya of $(6.4-70)\times10^{-7}$ relative to H, which is $2-3$ orders of magnitude below solar oxygen abundance. We conclude that gas-phase O$_2$ is not a major oxygen-carrier in TW Hya. Two other potential oxygen-carrying molecules, SO and SO$_2$, were covered in our observations, which we also do not detect. Additionally, we report a serendipitous detection of the C$^{15}$N $N_J = 2_{5/2}-1_{3/2}$ hyperfine transitions, $F = 3 - 2$ and $F = 2 - 1$, at 219.9 GHz, which we found via matched filtering and confirm through imaging.

We report the detection of a 13$\sigma$ H$\alpha$ emission line from HDF850.1 at $z=5.188\pm0.001$ using the FRESCO NIRCam F444W grism observations. Detection of H$\alpha$ in HDF850.1 is noteworthy, given its high far-IR luminosity, substantial dust obscuration, and the historical challenges in deriving its redshift. HDF850.1 shows a clear detection in the F444W imaging data, distributed between a northern and southern component, mirroring that seen in [CII] from the Plateau de Bure Interferometer. Modeling the SED of each component separately, we find that the northern component has a higher mass, star formation rate (SFR), and dust extinction than the southern component. The observed H$\alpha$ emission appears to arise entirely from the less-obscured southern component and shows a similar $\Delta$v$\sim$+130 km/s velocity offset to that seen for [CII] relative to the source systemic redshift. Leveraging H$\alpha$-derived redshifts from FRESCO observations, we find that HDF850.1 is forming in one of the richest environments identified to date at $z>5$, with 100 $z=5.17-5.20$ galaxies distributed across 10 structures and a $\sim$(15 cMpc)$^3$ volume. Based on the evolution of analogous structures in cosmological simulations, the $z=5.17-5.20$ structures seem likely to collapse into a single $>$10$^{14}$ $M_{\odot}$ cluster by $z\sim0$. Comparing galaxy properties forming within this overdensity with those outside, we find the masses, SFRs, and $UV$ luminosities inside the overdensity to be clearly higher. The prominence of H$\alpha$ line emission from HDF850.1 and other known highly-obscured $z>5$ galaxies illustrates the potential of NIRCam-grism programs to map both the early build-up of IR-luminous galaxies and overdense structures.

HDF850.1 is the brightest submillimeter galaxy (SMG) in the Hubble Deep Field. It is known as a heavily dust-obscured star-forming galaxy embedded in an overdense environment at $z = 5.18$. With nine-band NIRCam images at 0.8-5.0 $\mu$m obtained through the JWST Advanced Deep Extragalactic Survey (JADES), we detect and resolve the rest-frame UV-optical counterpart of HDF850.1, which splits into two components because of heavy dust obscuration in the center. The southern component leaks UV and H$\alpha$ photons, bringing the galaxy $\sim$100 times above the empirical relation between infrared excess and UV continuum slope (IRX-$\beta_\mathrm{UV}$). The northern component is higher in dust attenuation and thus fainter in UV and H$\alpha$ surface brightness. We construct a spatially resolved dust attenuation map from the NIRCam images, well matched with the dust continuum emission obtained through millimeter interferometry. The whole system hosts a stellar mass of $10^{11.0\pm0.1}\,\mathrm{M}_\odot$ and star-formation rate of $10^{3.0\pm0.2}\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$, placing the galaxy at the massive end of the star-forming main sequence at this epoch. We further confirm that HDF850.1 resides in a complex overdense environment at $z=5.17-5.30$, which hosts another luminous SMG at $z=5.30$ (GN10). The filamentary structures of the overdensity are characterized by 109 H$\alpha$-emitting galaxies confirmed through NIRCam slitless spectroscopy at 3.9-5 $\mu$m, of which only eight were known before the JWST observations. Given the existence of a similar galaxy overdensity in the GOODS-S field, our results suggest that $50\pm20$% of the cosmic star formation at $z=5.1-5.5$ occur in protocluster environments.

Dark energy is frequently modelled as an additional dynamical scalar field component in the Universe, referred to as "quintessence", which drives the late-time acceleration. Furthermore, the quintessence field may be coupled to dark matter and/or baryons, leading to a fifth force. In this paper we explore the consequences for non-linear cosmological structure formation arising from a momentum coupling between the quintessence field and dark matter only. The coupling leads to a modified Euler equation, which we implement in an N-body cosmological simulation. We then analyse the effects of the coupling on the non-linear power spectrum and the properties of the dark matter halos. We find that, for certain quintessence potentials, a positive coupling can lead to significantly reduced structure on small scales and somewhat enhanced structure on large scales, as well as reduced halo density profiles and increased velocity dispersions.

We introduce a model for the large-scale, global 3D structure of molecular clouds. Motivated by the morphological appearance of clouds in surface density maps, we model clouds as cylinders, with the aim of backing out information about the volume density distribution of gas and its relationship to star formation. We test our model by applying it to surface density maps for a sample of nearby clouds and find solutions that fit each of the observed radial surface density profiles remarkably well. Our most salient findings are that clouds with higher central volume densities are more compact and also have lower total mass. These same lower-mass clouds tend to have shorter gas depletion times, regardless of whether we consider their total mass or dense mass. Our analyses lead us to conclude that cylindrical clouds can be characterized by a universal structure that sets the timescale on which they form stars.

This study introduces novel constraints on the free-streaming of thermal relic warm dark matter (WDM) from Lyman-$\alpha$ forest flux power spectra. Our analysis utilises a high-resolution, high-redshift sample of quasar spectra observed using the HIRES and UVES spectrographs ($z=4.2-5.0$). We employ a Bayesian inference framework and a simulation-based likelihood that encompasses various parameters including the free-streaming of dark matter, cosmological parameters, the thermal history of the intergalactic medium, and inhomogeneous reionization, to establish lower limits on the mass of a thermal relic WDM particle of $5.7\;\mathrm{keV}$ (at 95\% C.L.). This result surpasses previous limits from the Lyman-$\alpha$ forest through reduction of the measured uncertainties due to a larger statistical sample and by measuring clustering to smaller scales ($k_{\rm max}=0.2\;\mathrm{km^{-1}\,s}$). The approximately two-fold improvement due to the expanded statistical sample suggests that the effectiveness of Lyman-$\alpha$ forest constraints on WDM models at high redshifts are limited by the availability of high-quality quasar spectra. Restricting the analysis to comparable scales and thermal history priors as in prior studies ($k_{\rm max}<0.1\;\mathrm{km^{-1}\,s}$) lowers the bound on the WDM mass to $4.1\;\mathrm{keV}$. As the precision of the measurements increases, it becomes crucial to examine the instrumental and modelling systematics. On the modelling front, we argue that the impact of the thermal history uncertainty on the WDM particle mass constraint has diminished due to improved independent observations. At the smallest scales, the primary source of modeling systematic arises from the structure in the peculiar velocity of the intergalactic medium and inhomogeneous reionization.

The $\Lambda$CDM paradigm successfully explains the large-scale structure of the Universe, but is less well constrained on sub-galactic scales. Gravitational lens modelling has been used to measure the imprints of dark substructures on lensed arcs, testing the small-scale predictions of $\Lambda$CDM. However, the methods required for these tests are subject to degeneracies among the lens mass model and the source light profile. We present a case study of the unique compound gravitational lens SDSSJ0946+1006, wherein a dark, massive substructure has been detected, whose reported high concentration would be unlikely in a $\Lambda$CDM universe. For the first time, we model the first two background sources in both I- and U-band HST imaging, as well as VLT-MUSE emission line data for the most distant source. We recover a lensing perturber at a $5.9\sigma$ confidence level with mass $\log_{10}(M_\mathrm{sub}/M_{\odot})=9.2^{+0.4}_{-0.1}$ and concentration $\log_{10}c=2.4^{+0.5}_{-0.3}$. The concentration is more consistent with CDM subhalos than previously reported, and the mass is compatible with that of a dwarf satellite galaxy whose flux is undetectable in the data at the location of the perturber. A wandering black hole with mass $\log_{10}(M_\mathrm{BH}/M_{\odot})=8.9^{+0.2}_{-0.1}$ is a viable alternative model. We systematically investigate alternative assumptions about the complexity of the mass distribution and source reconstruction; in all cases the subhalo is detected at around the $\geq5\sigma$ level. However, the detection significance can be altered substantially (up to $11.3\sigma$) by alternative choices for the source regularisation scheme.

We present an observational study of the dark matter fraction in star-forming disk-like galaxies up to redshift $z \sim 2.5$, selected from publicly available integral field spectroscropic surveys, namely KMOS3D}, KGES, and KROSS. We provide novel observational evidence, showing that at a fixed redshift, the dark matter fraction gradually increases with radius, indicating that the outskirts of galaxies are dark matter dominated, similarly to local star-forming disk galaxies. This observed dark matter fraction exhibits a decreasing trend with increasing redshift. However, on average, the fraction within the effective radius (upto outskirts) remains above 50\%, similar to locals. Furthermore, we investigated the relationships between the dark matter, baryon surface density, and circular velocity of galaxies. We observe a decreasing trend in the dark matter fraction as baryon surface densities increase, which is consistent across all stellar masses, redshift ranges, and radii, with a scatter of 0.13 dex. On the other hand, the correlation between the circular velocity at the outermost radius and the dark matter fraction within this radius has a relatively low scatter (0.11 dex), but its slope varies with stellar mass and with redshift, providing observational evidence of the dynamical evolution of the interplay between the baryonic and dark matter distributions with cosmic time. We observe that low stellar mass galaxies ($\log(M_{\star}/\mathrm{M_\odot}) \leq 10.0$) undergo a higher degree of evolution, which may be attributed to the hierarchical merging of galaxies.

Polarization observations through the next-generation large telescopes will be invaluable for exploring the magnetic fields and composition of jets in AGN, multi-messenger transients follow-up, and understanding interstellar dust and magnetic fields. The 25m Giant Magellan Telescope (GMT) is one of the next-generation large telescopes and is expected to have its first light in 2029. The telescope consists of a primary mirror and an adaptive secondary mirror comprising seven circular segments. The telescope supports instruments at both Nasmyth as well as Gregorian focus. However, none of the first or second-generation instruments on GMT has the polarimetric capability. This paper presents a detailed polarimetric modeling of the GMT for both Gregorian and folded ports for astronomical B-K filter bands and a field of view of 5 arc minutes. At 500nm, The instrumental polarization is 0.1% and 3% for the Gregorian and folded port, respectively. The linear to circular crosstalk is 0.1% and 30% for the Gregorian and folded ports, respectively. The Gregorian focus gives the GMT a significant competitive advantage over TMT and ELT for sensitive polarimetry, as these telescopes support instruments only on the Nasmyth platform. We also discuss a list of polarimetric science cases and assess science case requirements vs. the modeling results. Finally, we discuss the possible routes for polarimetry with GMT and show the preliminary optical design of the GMT polarimeter.

A major goal of proposed future space observatories, such as the Habitable World Observatory, is to directly image and characterize Earth-like planets around Sun-like stars to search for habitability signatures requiring the starlight suppression (contrast) of 1e-10. One of the significant aspects affecting this contrast is the polarization aberrations generated from the reflection from mirror surfaces. The polarization aberrations are the phase-dependent amplitude and phase patterns originating from the Fresnel reflections of the mirror surfaces. These aberrations depend on the angle of incidence and coating parameters of the surface. This paper simulates the polarization aberrations for an on-axis and off-axis TMA telescope of a 6.5 m monolithic primary mirror. We analyze the polarization aberrations and their effect on the coronagraphic performance for eight different recipes of mirror coatings for Astronomical filter bands g-I: three single-layer metal coatings and five recipes of protective coatings. First, the Jones pupils are estimated for each coating and filter band using the polarization ray tracing in Zemax. Then, we propagate these Jones pupils through a Vector Vortex Coronagraph and Perfect Coronagraphs using hcipy, a physical optics-based simulation framework. The analysis shows that the two main polarization aberrations generated from the four mirrors are the retardance-defocus and retardance-tilt. The simulations also show that the coating plays a significant role in determining the strength of the aberrations. The bare/oxi-aluminum and Al+18nm LiF coating outperforms all the other coatings by one order of magnitude.

In the development of space-based large telescope systems, having the capability to perform active optics correction allows correcting wavefront aberrations caused by thermal perturbations so as to achieve diffraction-limited performance with relaxed stability requirements. We present a method of active optics correction used for current ground-based telescopes and simulate its effectiveness for a large honeycomb primary mirror in space. We use a finite-element model of the telescope to predict misalignments of the optics and primary mirror surface errors due to thermal gradients. These predicted surface error data are plugged into a Zemax ray trace analysis to produce wavefront error maps at the image plane. For our analysis, we assume that tilt, focus and coma in the wavefront error are corrected by adjusting the pointing of the telescope and moving the secondary mirror. Remaining mid- to high-order errors are corrected through physically bending the primary mirror with actuators. The influences of individual actuators are combined to form bending modes that increase in stiffness from low-order to high-order correction. The number of modes used is a variable that determines the accuracy of correction and magnitude of forces. We explore the degree of correction that can be made within limits on actuator force capacity and stress in the mirror. While remaining within these physical limits, we are able to demonstrate sub-25 nm RMS surface error over 30 hours of simulated data. The results from this simulation will be part of an end-to-end simulation of telescope optical performance that includes dynamic perturbations, wavefront sensing, and active control of alignment and mirror shape with realistic actuator performance.

The Roman Coronagraph is expected to perform its high-order wavefront sensing and control (HOWFSC) with a ground-in-the-loop scheme due to the computational complexity of the Electric-Field-Conjugation (EFC) algorithm. This scheme provides the flexibility to alter the HOWFSC algorithm for given science objectives. A new alternative implicit-EFC algorithm is of particular interest as it requires no optical model to create a dark-hole, making the final contrast independent of the model accuracy. The intended HOWFSC scheme involves running EFC while observing a bright star such as $\zeta$ Puppis to create the initial dark-hole, then slew to the science target while maintaining the contrast with low-order WFSC over the given observation. Given a similar scheme, the efficacy of iEFC is simulated for two coronagraph modes, namely the Hybrid Lyot Coronagraph (HLC) and the wide-field-of-view Shaped-Pupil-Coronagraph (SPC-WFOV). End-to-end physical optics models for each mode serve as the tool for the simulations. Initial monochromatic simulations are presented and compared with monochromatic EFC results obtained with the FALCO software. Various sets of calibration modes are tested to understand the optimal modes to use when generating an iEFC response matrix. Further iEFC simulations are performed using broadband images with the assumption that $\zeta$ Puppis is the stellar object being observed. Shot noise, read noise, and dark current are included in the broadband simulations to determine if iEFC could be a suitable alternative to EFC for the Roman Coronagraph.

In order to reach contrast ratios of $10^{-8}$ and beyond, coronagraph testbeds need source optics that reliably emulate nearly-point-like starlight, with microfabricated pinholes being a compelling solution. To verify, a physical optics model of the Space Coronagraph Optical Bench (SCoOB) source optics, including a finite-difference time-domain (FDTD) pinhole simulation, was created. The results of the FDTD simulation show waveguide-like behavior of pinholes. We designed and fabricated microfabricated pinholes for SCoOB made from an aluminum overcoated silicon nitride film overhanging a silicon wafer substrate, and report characterization of the completed pinholes.

We report {\em JWST}/NIRCam measurements of quasar host galaxy emissions and supermassive black hole (SMBH) masses for six quasars at $5.9<z<7.1$ in the \textit{Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization} (EIGER) project. We obtain deep NIRCam imaging in the F115W, F200W, and F356W bands, as well as F356W grism spectroscopy of the quasars. We use bright unsaturated stars to construct models of the point spread function (PSF) and estimate the errors of these PSFs. We then measure or constrain the fluxes and morphology of the quasar host galaxies by fitting the quasar images as a point source plus an exponential disk. We successfully detect the host galaxy of three quasars, which have host-to-quasar flux ratios of $\sim1\%-5\%$. Spectral Energy Distribution (SED) fitting suggests that these quasar host galaxies have stellar masses of $M_*\gtrsim10^{10}M_\odot$. For quasars with host galaxy non-detections, we estimate the upper limits of their stellar masses. We use the grism spectra to measure the {\hb} line profile and the continuum luminosity, then estimate the SMBH masses for the quasars. Our results indicate that the positive relation between SMBH masses and host galaxy stellar masses already exists at redshift $z\gtrsim6$. The quasars in our sample show a high black hole to stellar mass ratio of $M_\text{BH}/M_*\sim0.15$, which is about $\sim1-2$ dex higher than the local relations. This result suggests that luminous quasars at $z\gtrsim6$ form a biased sample with overmassive black holes, which might have experienced early SMBH growth compared to their host galaxies' star formation.

We present a generalized Non-negative factorization (NMF)-based data reduction pipeline for circumstellar disk and exoplanet detection. By using an adaptable pre-processing routine that applies algorithmic masks and corrections to improper data, we are able to easily offload the computationally-intensive NMF algorithm to a graphics processing unit (GPU), significantly increasing computational efficiency. NMF has been shown to better preserve disk structural features compared to other post-processing approaches and has demonstrated improvements in the analysis of archival data. The adaptive pre-processing routine of this pipeline, which automatically aligns and applies image corrections to the raw data, is shown to significantly improve chromatic halo suppression. Utilizing HST-STIS and JWST-MIRI coronagraphic datasets, we demonstrate a factor of five increase in real-time computational efficiency by using GPUs to perform NMF compared to using CPUs. Additionally, we demonstrate the usefulness of higher numbers of NMF components with SNR and contrast improvements, which necessitates the use of a more computationally efficient approach for data reduction.

A few percent of red giants are enriched in Lithium with $A(\mathrm{Li}) > 1.5$. The evolutionary phase of the Li-rich red giants has remained uncertain because they could be placed both on the red-giant branch (RGB) near the bump luminosity and in the red clump (RC) region. However, thanks to asteroseismology, it has been found that most of them are actually RC stars. Starting at the bump luminosity, RGB progenitors of the RC stars experience extra mixing in the radiative zone separating the H-burning shell from the convective envelope followed by a series of convective He-shell flashes at the RGB tip, known as the He-core flash. Therefore, the He-core flash was proposed to cause fast extra mixing in the stars at the RGB tip that is needed for the Cameron-Fowler mechanism to produce Li. Alternatively, we propose that the RGB stars are getting enriched in Li by the same extra mixing that starts at the bump luminosity and initially leads to a decrease of the surface Li abundance but that is getting enhanced and begins to produce Li when the stars are approaching the RGB tip. We discuss five mechanisms of the RGB extra mixing, namely, the joint operation of rotation-driven meridional circulation and turbulent diffusion, the Azimuthal Magneto-Rotational Instability (AMRI), thermohaline convection, buoyancy of magnetic flux tubes, and internal gravity waves, and, based on results of (magneto-) hydrodynamics simulations, come to the conclusion that it is the mechanism of the AMRI that is most likely to support our hypothesis.

Integrated optical models allow for accurate prediction of the as-built performance of an optical instrument. Optical models are typically composed of a separate ray trace and diffraction model to capture both the geometrical and physical regimes of light. These models are typically separated across both open-source and commercial software that don't interface with each other directly. To bridge the gap between ray trace models and diffraction models, we have built an open-source optical analysis platform in Python called Poke that uses commercial ray tracing APIs and open-source physical optics engines to simultaneously model scalar wavefront error, diffraction, and polarization. Poke operates by storing ray data from a commercial ray tracing engine into a Python object, from which physical optics calculations can be made. We present an introduction to using Poke, and highlight the capabilities of two new propagation physics modules that add to the utility of existing scalar diffraction models. Gaussian Beamlet Decomposition is a ray-based approach to diffraction modeling that allows us to integrate physical optics models with ray trace models to directly capture the influence of ray aberrations in diffraction simulations. Polarization Ray Tracing is a ray-based method of vector field propagation that can diagnose the polarization aberrations in optical systems. Poke has been recently used to study the next generation of astronomical observatories, including the ground-based Extremely Large Telescopes and a 6 meter space telescope early concept for NASA's Habitable Worlds Observatory.

Judge et al. (2021) recently argued that a region of the solar spectrum in the near-UV between about 250 and 290 nm is optimal for studying magnetism in the solar chromosphere due to an abundance of Mg II, Fe II, and Fe I lines that sample various heights in the solar atmosphere. In this paper we derive requirements for spectropolarimetric instruments to observe these lines. We derive a relationship between the desired sensitivity to magnetic field and the signal-to-noise of the measurement from the weak-field approximation of the Zeeman effect. We find that many lines will exhibit observable polarization signals for both longitudinal and transverse magnetic field with reasonable amplitudes.

The role of energetic outflows from galactic nuclei in shaping galaxy formation and evolution is still shrouded in uncertainty. In this study, we shed light on this complex phenomenon by presenting evidence for a large-scale bipolar radio/X-ray-emitting bubble-like structure emanating from the central region of the nearby disk galaxy M106 (NGC 4258). Our findings, based on Low-Frequency Array survey data and Chandra observations, provide a glimpse into the underlying physical processes driving this enigmatic structure. Similar to the eROSITA/Fermi bubbles in our own Galaxy, the M106 bubbles enclose diffuse hot plasma and are partially bounded by prominent radio/X-ray-emitting edges. We constrain the magnetic field and cosmic-ray properties of the structure. The analysis of the X-ray data gives an estimate of the thermal energy of the bubbles as ~8 x 10^56 erg. This energy can be supplied by the jets and perhaps by the wind from the accretion flow of the galaxy's low-luminosity AGN, which most likely has been much more powerful in the recent past, with an average mechanical energy release rate of ~4 x 10^42 erg/s over the last ~ 8 x 10^6 yr -- the estimated age of the structure. We also show evidence for diffuse X-ray emission on larger scales, indicating the presence of a hot galactic corona. Our results provide a clear manifestation of galactic nuclear feedback regulating the gas content and energetics of the circumgalactic medium of disk galaxies similar to our own.

Continuous wavefront sensing on future space telescopes allows relaxation of stability requirements while still allowing on-orbit diffraction-limited optical performance. We consider the suitability of phase retrieval to continuously reconstruct the phase of a wavefront from on-orbit irradiance measurements or point spread function (PSF) images. As phase retrieval algorithms do not require reference optics or complicated calibrations, it is a preferable technique for space observatories, such as the Hubble Space Telescope or the James Webb Space Telescope. To increase the robustness and dynamic range of the phase retrieval algorithm, multiple PSF images with known amount of defocus can be utilized. In this study, we describe a recently constructed testbed including a 97 actuator deformable mirror, changeable entrance pupil stops, and a light source. The aligned system wavefront error is below ~30nm. We applied various methods to generate a known wavefront error, such as defocus and/or other aberrations, and found the accuracy and precision of the root mean squared error of the reconstructed wavefronts to be less than ~10nm and ~2nm, respectively. Further, we discuss the signal-to-noise ratios required for continuous dynamic wavefront sensing. We also simulate the case of spacecraft drifting and verify the performance of the phase retrieval algorithm for continuous wavefront sensing in the presence of realistic disturbances.

When star clusters are formed at low star-formation rates (SFRs), their stellar initial mass function (IMF) can hardly be filled continuously with stars at each mass. This lack holds for massive stars and is observationally verified by the correlation between star-cluster mass and its most massive cluster star. Since galaxy evolution is strongly affected by massive stars, numerical models should account for this lack. Because a filled IMF is mostly applied even when only fractions of massive stars form, here we investigate by 3D chemo-dynamical simulations of isolated dwarf galaxies how deviations from a standard IMF in star clusters affect the evolution. We compare two different IMF recipes, a filled IMF with one truncated at a maximum mass at which a single complete star forms. Attention is given to energetic and chemical feedback by massive stars. Since their energy release is mass dependent but steeper than the negative IMF slope, the energetic feedback retains a positive mass dependence, so that a filled IMF regulates SF stronger than truncated IMFs, though only stellar number fractions exist. The higher SFR of the truncated IMF in the simulation leads to more supernovae II (SNeII), driving galactic winds. Whether this results from the model-inherent larger SFR is questioned and therefore analytically explored. This shows the expected result for Lyman continuum, but that the total SNII energy release is equal for both IMF modes, while the power is smaller for the truncated IMF. Reasonably, the different IMFs leave fingerprints in abundance ratios of massive-to-intermediate-mass star elements.

The Large High Altitude Air Shower Observatory (LHAASO) has detected multiple ultra-high energy (UHE; E$_\gamma \ge$ 100 TeV) gamma-ray sources in the Milky Way Galaxy, which are associated with Galactic ``PeVatrons'' that accelerate particles up to PeV (= 10$^{15}$ eV) energies. Although supernova remnants (SNRs) and pulsar wind nebulae (PWNe), as source classes, are considered the leading candidates, further theoretical and observational efforts are needed to find conclusive proof to confirm the nature of these PeVatrons. This work aims to provide a phenomenological model to account for the emission observed from the direction of LHAASO J0341+5258, an unidentified UHE gamma-ray source observed by LHAASO. 15 years of Fermi-LAT data was analyzed to find the high energy (HE; 100 MeV $\le$ E$_\gamma$ $\le$ 100 GeV) GeV gamma-ray counterpart of LHAASO J0341+5258, in the 4FGL-DR3 catalog. We have explained the spectrum of the closest 4FGL source, 4FGL J0340.4+5302, by a synchro-curvature emission formalism typically used in the case of GeV pulsars. Escape-limited hadronic interaction between protons accelerated in an old, now invisible SNR and cold protons inside associated molecular clouds (MCs) and leptonic emission from a putative TeV halo was explored to explain the multi-wavelength (MWL) spectral energy distribution (SED) observed from the LHAASO source region. We have further discussed possible observational avenues that can be explored in the near future and predicted the outcome of those observational efforts from the model explored in this paper.

Large wind kinetic power of Wolf-Rayet (WR) stars make them ideal targets in low radio frequencies to search for non-thermal emission due to relativistic particle acceleration. In this paper, we present observations of two WR stars, WR 114 and WR 142, in Band 4 (550-950 MHz) and Band 5 (1050-1450 MHz) using the upgraded Giant Meterwave Radio Telescope (uGMRT). Neither star is detected in the observed frequency bands, nor extended emission associated with them. The upper limit to the free-free radio emission from the stellar wind enables us to constrain the mass-loss rate of WR 114 to $\lesssim \rm 10^{-5}\,M_{ \odot}\,yr^{-1}$; this is a factor three smaller than previously estimated using spectroscopic modelling. If we further assume that the WR stars are binaries, the non-detection of synchrotron emission from the putative wind collision region implies that the stars are either in very wide binary systems away from periastron, or that the stars are in close binary systems with an orbital separation $<70$ AU for WR 114 and $<20$ AU for WR 142. The non-detection of low-frequency radio emission from these two systems thus provides evidence that narrows their nature, though it does not rule them out as bonafide particle-accelerating colliding-wind binaries.

We employ our recent model of the cosmic-ray (CR) halo by Chernyshov et al. (2022) to compute the Galactic spectra of stable and unstable secondary nuclei. In this model, confinement of the Galactic CRs is entirely determined by the self-generated Alfvenic turbulence whose spectrum is controlled by nonlinear Landau damping. We analyze the physical parameters affecting propagation characteristics of CRs, and estimate the best set of free parameters providing accurate description of available observational data. We also show that agreement with observations at lower energies may be further improved by taking into account the effect of ion-neutral damping which operates near the Galactic disk.

TeV-range cosmic-ray electrons and positrons (CREs) have been directly measured in the search for new physics or unknown astrophysical sources. CREs can inverse-Compton scatter solar photons and boost their energies into gamma-ray bands. Any potential CRE excess would enhance the resultant inverse Compton emission spectrum in the relevant energy range, offering a new window to verify the measured CRE spectrum. In this paper, we show that an excess in the TeV range of the CRE spectrum, such as the one indicated by the DAMPE experiment, can induce a characteristic solar gamma-ray signal. Accounting for contamination from extragalactic gamma-ray backgrounds (EGB), we forecast the DAMPE feature is testable ($\gtrsim 4 \sigma$) with a $\sim 10^{5}\,\mathrm{m}^2\,{\rm yr}$ exposure in the off-disk direction. This can be achieved by long-exposure observations of water Cherenkov telescopes, such as LHAASO (7.2 years) and HAWC (25.9 years).

The evolution of the five largest satellites of Uranus during the crossing of the 5/3 mean motion resonance between Ariel and Umbriel is strongly affected by chaotic motion. Studies with numerical integrations of the equations of motion and analysis of Poincar\'e surface sections provided helpful insights to the role of chaos on the system. However, they lack of a quantification of this chaos in the phase-space. Here, we construct stability maps using the frequency analysis method. We determine that for low energies (small eccentricity and/or inclinations), the phase-space is mainly stable. As the energy increases, the chaotic regions replace the stable motion, until only small, localized libration regions remain stable.

G\,29$-$38 (TIC~422526868) is one of the brightest ($V=13.1$) and closest ($d = 17.51$\,pc) pulsating white dwarfs with a hydrogen-rich atmosphere (DAV/ZZ Ceti class). It was observed by the {\sl TESS} spacecraft in sectors 42 and 56. The atmosphere of G~29$-$38 is polluted by heavy elements that are expected to sink out of visible layers on short timescales. The photometric {\sl TESS} data set spans $\sim 51$ days in total, and from this, we identified 56 significant pulsation frequencies, that include rotational frequency multiplets. In addition, we identified 30 combination frequencies in each sector. The oscillation frequencies that we found are associated with $g$-mode pulsations, with periods spanning from $\sim$ 260 s to $\sim$ 1400 s. We identified %three distinct rotational frequency triplets with a mean separation $\delta \nu_{\ell=1}$ of 4.67 $\mu$Hz and a quintuplet with a mean separation $\delta \nu_{\ell=2}$ of 6.67 $\mu$Hz, from which we estimated a rotation period of about $1.35 \pm 0.1$ days. We determined a constant period spacing of 41.20~s for $\ell= 1$ modes and 22.58\,s for $\ell= 2$ modes. We performed period-to-period fit analyses and found an asteroseismological model with $M_{\star}/M_{\odot}=0.632 \pm 0.03$, $T_{\rm eff}=11\, 635\pm 178$ K, and $\log{g}=8.048\pm0.005$ (with a hydrogen envelope mass of $M_{\rm H}\sim 5.6\times 10^{-5}M_{\star}$), in good agreement with the values derived from spectroscopy. We obtained an asteroseismic distance of 17.54 pc, which is in excellent agreement with that provided by {\sl Gaia} (17.51 pc).

We use the three-dimensional Monte Carlo radiative transfer code HDUST to model Be stars where the disc is tilted from the equatorial plane of the star. We compute 128 models across 4 spectral types, B0, B2, B5 and B8, tilting the disc by $0^o$, $10^o$, $20^o$, and $40^o$, while varying disc density according to spectral type. We also compute every model for an average and high stellar rotation rate. We first discuss non-tilted disc temperatures and show its non-linear dependence on stellar and disc parameters. We find that tilting the disc minimally affects the density-weighted average disc temperature, but tilting does create a temperature asymmetry in disc cross sections, which is more pronounced for a faster rotation rate. We also investigate the effect tilting has on $V$-band magnitude, polarization, and the H$\alpha$ line. Tilting the disc does affect these observables, but the changes are entirely dependent on the position of the observer relative to the direction of tilt. We find the observables that distinguish tilting from a change in density or geometry are the H$\alpha$ line shape, where it can transition between single-peaked and double-peaked, and the polarization position angle, whose value is dependent on the projected major elongation axis of the disc on the sky. We also present one early and one late-type model with warped discs. We find their temperature structure varies a small amount from the uniformly tilted models, and the different observables correspond to different tilt angles, consistent with their expected volume of origin within the disc.

Published near-infrared spectra of the four largest classical Uranian satellites display the presence of discrete deposits of CO$_2$ ice, along with subtle absorption features around 2.2 $\mu$m. The two innermost satellites, Miranda and Ariel, also possess surfaces heavily modified by past endogenic activity. Previous observations of the smallest satellite, Miranda, have not detected the presence of CO$_2$ ice, and a report of an absorption feature at 2.2 $\mu$m has not been confirmed. An absorption feature at 2.2 $\mu$m could result from exposed or emplaced NH$_3$- or NH$_4$-bearing species, which have a limited lifetime on Miranda's surface, and therefore may imply that Miranda's internal activity was relatively recent. In this work, we analyzed near-infrared spectra of Miranda to determine whether CO$_2$ ice and the 2.2-$\mu$m feature are present. We measured the band area and depth of the CO$_2$ ice triplet (1.966, 2.012, and 2.070 $\mu$m), a weak 2.13-$\mu$m band attributed to CO$_2$ ice mixed with H$_2$O ice, and the 2.2-$\mu$m band. We confirmed a prior detection of a 2.2-$\mu$m band on Miranda, but we found no evidence for CO$_2$ ice, either as discrete deposits or mixed with H$_2$O ice. We compared a high signal-to-noise spectrum of Miranda to synthetic and laboratory spectra of various candidate compounds to shed light on what species may be responsible for the 2.2-$\mu$m band. We conclude that the 2.2-$\mu$m absorption is best matched by a combination of NH$_3$ ice with NH$_3$-hydrates or NH$_3$-H$_2$O mixtures. NH$_4$-bearing salts like NH$_4$Cl are also promising candidates that warrant further investigation.

We examine the effects of a stiff pre-recombination era on the present day's energy spectrum of the primordial gravitational waves. If the background total equation of state parameter at the pre-recombination era is described by a kination era one, this directly affects the modes with characteristic wavenumbers which reenter the Hubble horizon during this stiff era. The stiff era causes a broken-power-law effect on the energy spectrum of the gravitational waves. We use two approaches, one model agnostic and a specific model that can realize this scenario. In all cases, the inflationary era can be realized either by some theory leading to a standard red-tilted tensor spectral index or by some theory which has a mild tensor spectral index $n_{\mathcal{T}}=0.17-3.7$ like an Einstein-Gauss-Bonnet theory. For the model agnostic scenario case, the NANOGrav signal can be explained by the stiff pre-recombination era combined with an inflationary era with a mild blue-tilted tensor spectral index $n_{\mathcal{T}}=3.7$ and a low-reheating temperature $T_R\sim 0.1$GeV. In the same case, the red-tilted inflationary theory signal can be detectable by the future LISA, BBO and DECIGO experiments. The model dependent approach is based on a Higgs-axion model which can yield multiple deformations of the background total equation of state parameter, causing multiple broken-power-law behaviors occurring in various eras before and after the recombination era. In this case, the NANOGrav signal is explained by this model in conjunction with an inflationary era with a really mild blue-tilted tensor spectral index $n_{\mathcal{T}}=0.17$ and a low-reheating temperature $T_R\sim 20\,$GeV. In this case, the signal can be detectable by the future Litebird experiment, which is a very characteristic pattern in the tail of the primordial gravitational wave energy spectrum.

Gas rich galaxies are typically star-forming. We make use of HI-MaNGA, a program of HI follow-up for the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey of the Sloan Digital Sky Surveys to construct a sample of unusual neutral hydrogen (HI, 21cm) rich galaxies which have low Star Formation Rates (SFRs); using infra-red color from the Wide-field Infrared Survey Explorer (WISE) as a proxy for specific SFR. Out of a set of 1575 MaNGA galaxies with HI-MaNGA detections, we find 83 (5%) meet our selection criteria to be HI rich with low SFR. We construct two stellar mass-matched control samples: HI rich galaxies with typical SFR (High SF Control) and HI poor galaxies with low SFR (Low HI Control). We investigate the properties of each of these samples, comparing physical parameters such as ionization state maps, stellar and ionized gas velocity and dispersion, environment measures, metallicity, and morphology to search for the reasons why these unusual HI rich galaxies are not forming stars. We find evidence for recent external accretion of gas in some galaxies (via high counter-rotating fractions), along with some evidence for AGN feedback (from a high cLIER and/or red geyser fraction), and bar quenching (via an enhanced strong bar fraction). Some galaxies in the sample are consistent with simply having their HI in a high angular momentum, large radius, low density disc. We conclude that no single physical process can explain all HI rich, low SFR galaxies.

Mountains or non-axisymmetric deformations of rotating neutron stars (NS) efficiently radiate gravitational waves (GW). We consider analogies between NS mountains and surface features of solar system bodies. Both NS and moons such as Europa or Enceladus have thin crusts over deep oceans while Mercury has a thin crust over a large metallic core. Thin sheets may wrinkle in universal ways. Europa has linear features, Enceladus has ``Tiger" stripes, and Mercury has lobate scarps. NS may have analogous features. The innermost inner core of the Earth is anisotropic with a shear modulus that depends on direction. If NS crust material is also anisotropic this will produce an ellipticity, when the crust is stressed, that grows with spin frequency. This yields a breaking index (log derivative of spin down rate) very different from $n=5$ and could explain the maximum spin observed for neutron stars and a possible minimum ellipticity of millisecond pulsars.

We previously reported a rare class of variable star light curves isolated from a sample of 4.7 million candidate variables from the ATLAS survey. Dubbed `UCBH' light curves, they have broad minima and narrow, symmetrical maxima, with typical periods of 1-10 days and amplitudes of 0.05--0.20 mag. They maintain constant amplitude, shape, and phase coherence over multiple years, but do not match any known class of pulsating variables. A localized bright spot near the equator of a rotating star will produce a UCBH-type light curve for most viewing geometries. Most stars that exhibit rotational variability caused primarily by a single bright feature should therefore appear as UCBH stars, although a rotating bright spot is not the only thing that could produce a UCBH-type lightcurve. We have spectroscopically investigated fourteen UCBH stars and found ten of them to be Ap/Bp stars: A-type or B-type stars with greatly enhanced photospheric abundances of specific heavy elements. Rotationally variable Ap/Bp stars are referred to as $\alpha^2$ CVn variables. Most ATLAS UCBH stars are therefore $\alpha^2$ CVn stars, although only a minority of $\alpha^2$ CVn stars in the literature have UCBH light curves. The fact that $\alpha^2$ CVn stars dominate the UCBH class suggests that lone bright spots with sufficient size and contrast develop more readily on Ap/Bp stars than on any other type. The $\alpha^2$ CVn UCBH stars may be characterized by a specific magnetic field topology, making them intriguing targets for future Zeeman-Doppler imaging.

Coronagraphic instruments provide a great chance of enabling high contrast spectroscopy for the pursuit of finding a habitable world. Future space telescope coronagraph instruments require high performing focal plane masks in combination with precise wavefront sensing and control techniques to achieve dark holes for planet detection. Several wavefront control algorithms have been developed in recent years that might vary in performance depending on the coronagraph they are paired with. This study compares 3 model-free and model-based algorithms when coupled with either a Vector (VVC) or a Scalar (SVC) Vortex Coronagraph mask in the same laboratory conditions: Pairwise Probing with Electric Field Conjugation, the Self-Coherent Camera with Electric Field Conjugation, and Implicit Electric Field Conjugation. We present experimental results from the In-Air Coronagraph Testbed (IACT) at JPL in narrowband and broadband light, comparing the pros and cons of each of these wavefront sensing and control algorithms with respect to their potential for future space telescopes.

The utilization of a 6.5m monolithic primary mirror in a compact three-mirror anastigmat (TMA) telescope design offers unprecedented capabilities to accommodate various next generation science instruments. This design enables the rapid and efficient development of a large aperture telescope without segmented mirrors while maintaining a compact overall form factor. With its exceptional photon collection area and diffraction-limited resolving power, the TMA design is ideally suited for both the ground and space active/adaptive optics concepts, which require the capture of natural guide stars within the field of view for wavefront measurement to correct for misalignments and shape deformation caused by thermal gradients. The wide field of view requirement is based on a statistical analysis of bright natural guide stars available during observation. The primary mirror clear aperture, compactness requirement, and detector pixel sizes led to the choice of TMA over simpler two-mirror solutions like Ritchey-Chretien (RC) telescopes, and the TMA design offers superior diffraction-limited performance across the entire field of view. The standard conic surfaces applied to all three mirrors (M1, M2, and M3) simplify the optical fabrication, testing, and alignment process. Additionally, the TMA design is more tolerant than RC telescopes. Stray light control is critical for UV science instrumentation, and the field stop and Lyot stop are conveniently located in the TMA design for this purpose.

The detection and characterization of Earth-like exoplanets around Sun-like stars is a primary science motivation for the Habitable Worlds Observatory. However, the current best technology is not yet advanced enough to reach the 10^-10 contrasts at close angular separations and at the same time remain insensitive to low-order aberrations, as would be required to achieve high-contrast imaging of exo-Earths. Photonic technologies could fill this gap, potentially doubling exo-Earth yield. We review current work on photonic coronagraphs and investigate the potential of hybridized designs which combine both classical coronagraph designs and photonic technologies into a single optical system. We present two possible systems. First, a hybrid solution which splits the field of view spatially such that the photonics handle light within the inner working angle and a conventional coronagraph that suppresses starlight outside it. Second, a hybrid solution where the conventional coronagraph and photonics operate in series, complementing each other and thereby loosening requirements on each subsystem. As photonic technologies continue to advance, a hybrid or fully photonic coronagraph holds great potential for future exoplanet imaging from space.

New development approaches, including launch vehicles and advances in sensors, computing, and software, have lowered the cost of entry into space, and have enabled a revolution in low-cost, high-risk Small Satellite (SmallSat) missions. To bring about a similar transformation in larger space telescopes, it is necessary to reconsider the full paradigm of space observatories. Here we will review the history of space telescope development and cost drivers, and describe an example conceptual design for a low cost 6.5 m optical telescope to enable new science when operated in space at room temperature. It uses a monolithic primary mirror of borosilicate glass, drawing on lessons and tools from decades of experience with ground-based observatories and instruments, as well as flagship space missions. It takes advantage, as do large launch vehicles, of increased computing power and space-worthy commercial electronics in low-cost active predictive control systems to maintain stability. We will describe an approach that incorporates science and trade study results that address driving requirements such as integration and testing costs, reliability, spacecraft jitter, and wavefront stability in this new risk-tolerant "LargeSat" context.

Several different methods are regularly used to infer the properties of the neutral interstellar medium (ISM) using atomic hydrogen (H I) 21cm absorption and emission spectra. In this work, we study various techniques used for inferring ISM gas phase properties, namely the correlation between brightness temperature and optical depth $(T_B(v)$, $\tau(v))$ at each channel velocity $(v)$, and decomposition into Gaussian components, by creating mock spectra from a 3D magnetohydrodynamic simulation of a two-phase, turbulent ISM. We propose a physically motivated model to explain the $T_B(v)-\tau(v)$ distribution and relate the model parameters to properties like warm gas spin temperature and cold cloud length scales. Two methods based on Gaussian decomposition -- using only absorption spectra and both absorption and emission spectra -- are used to infer the column density distribution as a function of temperature. In observations, such analysis reveals the puzzle of large amounts (significantly higher than in simulations) of gas with temperature in the thermally unstable range of $\sim200\mathrm{\ K}$ to $\sim2000\mathrm{\ K}$ and a lack of the expected bimodal (two-phase) temperature distribution. We show that, in simulation, both methods are able to recover the true gas distribution till temperatures $\lesssim2500\mathrm{\ K}$ (and the two-phase distribution in general) reasonably well. We find our results to be robust to a range of effects such as noise, varying emission beam size, and simulation resolution. This shows that the observational inferences are unlikely to be artifacts, thus highlighting a tension between observations and simulations. We discuss possible reasons for this tension and ways to resolve it.

This document collects the contributions of the KM3NeT collaboration to the ICRC2023 conference, held from July 26 to August 3, 2023, in Nagoya, Japan. KM3NeT submitted 38 contributions to ICRC2023, on neutrino- and multimessenger astronomy, neutrino oscillation physics, cosmic ray physics, searches for dark matter and exotics, calibration, technical detector descriptions, and art. Proceedings are published in Proceedings of Science.

The development of techniques to measure accurately black hole spins is crucial to study the physics and astrophysics of these objects. X-ray reflection spectroscopy is currently the most popular method to estimate the spins of accreting black holes; so far it has provided a spin measurement of about 40 stellar-mass black holes in X-ray binaries and 40 supermassive black holes in active galactic nuclei. The relativistic precession model is another method to measure the spins of stellar-mass black holes: it requires the measurement of the frequencies of three simultaneous quasi-periodic oscillations and can potentially provide very precise estimates of the black hole mass and spin. However, the two methods do not seem to provide consistent results when applied to the same sources, which questions the actual reliability and accuracy of these measurements. Recently, the relativistic precession model has been applied to infer the spin of the black hole in XTE J1859+226 (Motta et al. 2022). The authors found $a_* = 0.149 \pm 0.005$ (68% CL) and there are no other spin measurements of this source. We looked for archived RXTE observations of XTE J1859+226 with blurred reflection features and found 23 spectra suitable for measuring the black hole spin. We employed two different models with relxill and relxillD and obtained a higher spin value from all these fits. From the further simultaneous fitting performed on two different subsets of the total set of 23 spectra, we infer the black hole spin parameter $a_* = 0.986^{+0.001}_{-0.004}$ and $a_* =0.987 \pm 0.003$ (90% CL, statistical) with relxill and relxillD for the first set and, $a_*$ = $0.981_{-0.007}^{+0.006}$ and $0.982_{-0.007}^{+0.006}$ (90% CL, statistical) for the second set, respectively. This clearly confirms the discrepancy between the black hole spin measurements inferred from the two techniques.

We present a comprehensive data set of supernova (SN) 2016adj located within the central dust lane of Centaurus A. SN 2016adj is significantly reddened and after correcting the peak apparent $B$-band magnitude ($m_B = 17.48\pm0.05$) for Milky Way reddening and our inferred host-galaxy reddening parameters (i.e., $R_{V}^{host} = 5.7\pm0.7$ and $A_{V}^{host} = 6.3\pm0.2$), we estimate it reached a peak absolute magnitude of $M_B \sim -18$. Detailed inspection of the optical/NIR spectroscopic time-series reveals a carbon-rich SN Ic and not a SN Ib/IIb as previously suggested in the literature. The NIR spectra shows prevalent carbon-monoxide formation occurring already by +41 days past $B$-band maximum, which is $\approx 11$ days earlier than previously reported in the literature for this object. Interestingly around two months past maximum, the NIR spectrum of SN~2016adj begins to exhibit H features, with a +97~d medium resolution spectrum revealing both Paschen and Bracket lines with absorption minima of $\sim 2000$ km/s, full-width-half-maximum emission velocities of $\sim 1000$ km/s, and emission line ratios consistent with a dense emission region. We speculate these attributes are due to circumstellar interaction (CSI) between the rapidly expanding SN ejecta and a H-rich shell of material formed during the pre-SN phase. A bolometric light curve is constructed and a semi-analytical model fit suggests the supernova synthesized 0.5 solar masses of $^{56}$Ni and ejected 4.2 solar masses of material, though these values should be approached with caution given the large uncertainties associated with the adopted reddening parameters, possible CSI contamination, and known light echo emission. Finally, inspection of Hubble Space Telescope archival data yielded no progenitor detection.

This study presents structural and fundamental astrophysical parameters of poorly studied open cluster NGC 2509. We used the third photometric and astrometric data release of the Gaia (Gaia DR3) to perform analyses. By taking into account the Gaia DR3 astrometric data, we calculated the membership probabilities of the stars in the region of NGC 2509. As a result of the membership analysis, 244 stars with membership probabilities $P\geq50$% were determined as the physical members of the cluster. The colour excess, distance and age were obtained simultaneously by fitting solar metallicity PARSEC isochrones to $G\times G_{\rm BP}-G_{\rm RP}$ colour-magnitude diagram. We considered the most likely cluster member stars during the fitting procedure and calculated the colour excess, distance and age of the NGC 2509 as $E(G_{\rm BP}-G_{\rm RP})=0.100\pm0.015$ mag, $d=2518\pm667$ pc and $t=1.5\pm 0.1$ Gyr, respectively.

The formation mechanism of high-concentration dwarf galaxies is still a mystery. We perform a comparative study of the intrinsic shape of nearby low-mass galaxies with different stellar concentration. The intrinsic shape is parameterized by the intermediate-to-major axis ratios B/A and the minor-to-major axis ratios C/A of triaxial ellipsoidal models. Our galaxies ($10^{7.5} M_\odot$ < $M_\star$ < $10^{10.0} M_\odot$) are selected to have spectroscopic redshift from SDSS or GAMA, and have broadband optical images from the HSC-SSP Wide layer survey. The deep HSC-SSP images allow to measure the apparent axis ratios $q$ at galactic radii beyond the central star-forming area of our galaxies. We infer the intrinsic axis ratios based on the $q$ distributions. We find that 1) our galaxies have typical intrinsic shape similarly close to be oblate ($\mu_{B/A}$ $\sim$ 0.9--1), regardless of the concentration, stellar mass, star formation activity, and local environment (being central or satellite); 2) galaxies with the highest concentration tend to have intrinsic thickness similar to or (in virtually all cases) slightly thinner (i.e. smaller mean $\mu_{C/A}$ or equivalently lower triaxiality) than ordinary galaxies, regardless of other properties explored here. This appears to be in contrast with the expectation of the classic merger scenario for high-concentration galaxies. Given the lack of a complete understanding of dwarf-dwarf merger, we cannot draw a definite conclusion about the relevance of mergers in the formation of high-concentration dwarfs. Other mechanisms such as halo spin may also play important roles in the formation of high-concentration dwarf galaxies.

Heartbeat stars (HBSs) with tidally excited oscillations (TEOs) are ideal laboratories for studying the effect of equilibrium and dynamical tides. However, studies of TEOs in Kepler HBSs are rare due to {the need for better modeling of the equilibrium tide in light curves}. We revisit the HBSs reported by Li et al. {and study the TEOs in these HBSs based on the derived orbital parameters that could express the equilibrium tide.} We also compile a set of analysis procedures to examine the harmonic and anharmonic TEOs in their Fourier spectrum. The TEOs of 45 HBSs (excluding eight systems studied in previous works) have been determined and presented. 19 of them show prominent TEOs (the signal-to-noise ratio of the harmonics $S/N \ge 10$). The relation between the orbital eccentricities and the harmonic number of the TEOs shows a positive correlation. The relation between the orbital periods and the harmonic number also shows a positive correlation. Furthermore, the distribution of HBSs with TEOs in the Hertzsprung-Russell (H-R) diagram shows that TEOs are more visible in hot stars with surface temperatures $T$ $\gtrsim$ 6500 K. These samples may also be valuable targets for future studies of the effect of tidal action in eccentric orbits.

Nitrogen, carbon monoxide, and methane are key materials in the far outer Solar System where their high volatility enables them to sublimate, potentially driving activity at very low temperatures. Knowledge of their vapor pressures and latent heats of sublimation at relevant temperatures is needed to model the processes involved. We describe a method for using a quartz crystal microbalance to measure the sublimation flux of these volatile ices in the free molecular flow regime, accounting for the simultaneous sublimation from and condensation onto the quartz crystal to derive vapor pressures and latent heats of sublimation. We find vapor pressures to be somewhat lower than previous estimates in literature, with carbon monoxide being the most discrepant of the three species, almost an order of magnitude lower than had been thought. These results have important implications across a variety of astrophysical and planetary environments.

James Webb Space Telescope's NIRSpec infrared imaging spectrometer observed the outer solar system dwarf planets Eris and Makemake in reflected sunlight at wavelengths spanning 1 through 5 microns. Both objects have high albedo surfaces that are rich in methane ice, with a texture that permits long optical path lengths through the ice for solar photons. There is evidence for N2 ice absorption around 4.2 um on Eris, though not on Makemake. No CO ice absorption is seen at 4.67 um on either body. For the first time, absorption bands of two heavy isotopologues of methane are observed at 2.615 um (13CH4), 4.33 um (12CH3D), and 4.57 um (12CH3D). These bands enable us to measure D/H ratios of (2.5 +/- 0.5) x 10-4 and (2.9 +/- 0.6) x 10-4, along with 13C/12C ratios of 0.012 +/- 0.002 and 0.010 +/- 0.003 in the surface methane ices of Eris and Makemake, respectively. The measured D/H ratios are much lower than that of presumably primordial methane in comet 67P/Churyumov-Gerasimenko, but they are similar to D/H ratios in water in many comets and larger outer solar system objects. This similarity suggests that the hydrogen atoms in methane on Eris and Makemake originated from water, indicative of geochemical processes in past or even ongoing hot environments in their deep interiors. The 13C/12C ratios are consistent with commonly observed solar system values, suggesting no substantial enrichment in 13C as could happen if the methane currently on their surfaces was the residue of a much larger inventory that had mostly been lost to space. Possible explanations include geologically recent outgassing from the interiors as well as processes that cycle the surface methane inventory to keep the uppermost surfaces refreshed.

we review the Baryon Acoustic Oscillations from Integrated Neutral Gas Observations (BINGO) telescope, an international collaboration, led by Brazil and China, aiming to explore the Universe history through integrated post-reionization 21cm signals and fast radio emissions. For identifying individually fast radio sources, the Advanced Bingo Dark Universe Studies (ABDUS) project has been proposed and developed and will combine the current BINGO construction with the main single-dish telescope and stations of phased-array and outrigger.

Dynamic fibrils (DFs) are commonly observed chromospheric features in solar active regions. Recent observations from the Extreme Ultraviolet Imager (EUI) aboard the Solar Orbiter have revealed unambiguous signatures of DFs at the coronal base, in extreme ultraviolet (EUV) emission. However, it remains unclear if the DFs detected in the EUV are linked to their chromospheric counterparts. Simultaneous detection of DFs from chromospheric to coronal temperatures could provide important information on their thermal structuring and evolution through the solar atmosphere. In this paper, we address this question by using coordinated EUV observations from the Atmospheric Imaging Assembly (AIA), Interface Region Imaging Spectrograph (IRIS), and EUI to establish a one-to-one correspondence between chromospheric and transition region DFs (observed by IRIS) with their coronal counterparts (observed by EUI and AIA). Our analysis confirms a close correspondence between DFs observed at different atmospheric layers, and reveals that DFs can reach temperatures of about 1.5 million Kelvin, typical of the coronal base in active regions. Furthermore, intensity evolution of these DFs, as measured by tracking them over time, reveals a shock-driven scenario in which plasma piles up near the tips of these DFs and, subsequently, these tips appear as bright blobs in coronal images. These findings provide information on the thermal structuring of DFs and their evolution and impact through the solar atmosphere.

The Habitable Worlds Observatory mission will require coronagraphs capable of achieving contrasts of 1e-10 to detect exo-Earths. The choice of coronagraph depends on finding a solution that is achromatic within a 20\% bandwidth, insensitive to low order aberrations and polarization independent. We present two scalar vortex phase mask designs which employ a Roddier phase dimple and a dual zone phase dimple to improve the achromatic performance by addressing the chromatic stellar leakage not handled by the vortex. We show that using these dimples, it is possible to substantially improve the broadband contrast performance of existing scalar vortex phase masks.

We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M$_\odot$ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the strength of the magnetic fields in the core, we compute four models. Our results show that the MHD outflows are produced only for models (two out of four), to which rapid rotation and strong magnetic fields are initially imposed. Seen from the direction perpendicular to the rotational axis, a characteristic waveform is obtained exhibiting a monotonic time increase in the wave amplitude. As previously identified, this stems from the propagating MHD outflows along the axis. We show that the GW amplitude from anisotropic neutrino emission becomes more than one order-of-magnitude bigger than that from the matter contribution, whereas seen from the rotational axis, both of the two components are in the same order-of-magnitudes. Due to the memory effect, the frequency of the neutrino GW from our full-fledged 3D-MHD models is in the range less than $\sim$10Hz. Toward the future GW detection for a Galactic core-collapse supernova, if driven by the MR mechanism, the planned next-generation detector as DECIGO is urgently needed to catch the low-frequency signals.

The EUI instrument on the Solar Orbiter spacecraft has obtained the most stable, high-resolution images of the solar corona from its orbit with a perihelion near 0.4 AU. A sequence of 360 images obtained at 17.1 nm, between 25-Oct-2022 19:00 and 19:30 UT is scrutinized. One image pixel corresponds to 148 km at the solar surface. The widely-held belief that the outer atmosphere of the Sun is in a continuous state of magnetic turmoil is pitted against the EUI data. The observed plasma variations appear to fall into two classes. By far the dominant behavior is a very low amplitude variation in brightness (1%) in the coronal loops, with larger variations in some footpoint regions. No hints of observable changes in magnetic topology are associated with such small variations. The larger amplitude, more rapid, rarer and less-well organized changes are associated with flux emergence. It is suggested therefore that while magnetic reconnection drives the latter, most of the active corona is heated with no evidence of a role for large-scale (observable) reconnection. Since most coronal emission line widths are subsonic, the bulk of coronal heating, if driven by reconnection, can only be of tangentially discontinuous magnetic fields, with angles below about $0.5c_S/c_A \sim 0.3\beta$, with $\beta$ the plasma beta parameter ($\sim 0.01)$, and $c_S$ and $c_A$ sound and Alfv\'en speeds. If heated by multiple small flare-like events, then these must be $\lesssim 10^{21}$ erg, i.e. pico-flares. But processes other than reconnection have yet to be ruled out, such as viscous dissipation, which may contribute to the steady heating of coronal loops over active regions.

Several aspects of the magnetospheric physics of magnetars are summarized, including: GeV and hard X-ray emissions of magnetars, timing behaviors during magnetar outburst (soft X-ray observations), optical/IR observations of magnetars, radio emission of magnetars, and accreting magnetars. A unified picture for pulsars and magnetars are adopted, especially wind braking of magnetars, magnetar+ fallback disk systems, twisted dipole magnetic field, and accreting low magnetic field magnetars etc. It is pointed out that magnetars are related to a broad range of astrophysical phenomena.

On 2019 April 25, the LIGO/Virgo Scientific Collaboration detected a compact binary coalescence, GW190425. Under the assumption of the binary neutron star (BNS), the total mass of $3.4^{+0.3}_{-0.1}\, M_\odot$ lies five standard deviations away from the known Galactic population mean. In the standard common envelope scenario, the immediate progenitor of GW190425 is a close binary system composed of an NS and a He-rich star. With the detailed binary evolutionary modeling, we find that in order to reproduce GW190425-like events, super-Eddington accretion (e.g., $1,000\,\dot{M}_{\rm Edd}$) from a He-rich star onto the first-born NS with a typical mass of 1.33 $M_\odot$ via stable Case BB mass transfer (MT) is necessarily required. Furthermore, the immediate progenitors should potentially have an initial mass of $M_{\rm ZamsHe}$ in a range of $3.0-3.5$ $M_\odot$ and an initial orbital period of $P_{\rm init}$ from 0.08 days to 0.12 days, respectively. The corresponding mass accreted onto NSs via stable Case BB MT phase varies from $0.70\, M_\odot$ to $0.77\, M_\odot$. After the formation of the second-born NS, the BNSs are expected to be merged due to gravitational wave emission from $\sim$ 11 Myr to $\sim$ 190 Myr.

The shear flow influences the stability of magnetohydrodynamic (MHD) waves. In the presence of a dissipation mechanism, flow shear may induce a MHD wave instability below the threshold of the Kelvin-Helmholtz instability (KHI), which is called dissipative instability (DI). This phenomenon is also called negative energy wave instability (NEWI) because it is closely related to the backward wave which has negative wave energy. Considering viscosity as a dissipation mechanism, we derive an analytical dispersion relation for the slow sausage modes in a straight cylinder with a discontinuous boundary. It is assumed that the steady flow is inside and dynamic and bulk viscosities are outside the circular flux tube under photospheric condition. When the two viscosities are weak, it is found that for the slow surface mode, the growth rate is proportional to the axial wavenumber and flow shear, consistent with in the incompressible limit. For a slow body mode, the growth rate has a peak at certain axial wavenumber and its order of magnitude is similar to surface mode. The linear relationship between the growth rate and the dynamic viscosity established in the incompressible limit develops nonlinearly when the flow shear and/or the two viscosities are sufficiently strong.

Context. We study the relation between the known binary fraction and spectral absorption feature index to judge whether (and potentially which) spectral absorption feature indices are suitable for determining the binary fraction. Aims. The determination of the binary fraction is important in studies of binary star formation, evolutionary population synthesis models, and other works. The number of binary stars is difficult to determine for nearly all stellar systems because the individual stars are need to be resolved photometrically or spectroscopically. By comparison, their integrated spectra or spectral absorption feature indices are relatively easy to obtain. Results. We find that the low-resolution (15\,\AA) spectrum is not suitable for this study and the binary fraction type would affect the results: $f$($q$>0.5) and $f$(tot)$^{\rm mc}$ exhibit better correlations with the spectral absorption feature index than $f$(tot)$^{\rm mf}$ and the difference in metallicity would significantly affect the above relationship. % Finally, to eliminate the effects of metallicity, age, and dynamical evolution, we only used those GCs with multiple spectra observed among different regions. % We find that OIII-1, OIII-2, H$_{\rm \gamma F}$, H$_{\rm \delta F}$, H$_{\rm \gamma A}$, H$_{\rm \delta A}$, H$_{\rm \beta}$, Ca4455, C$_2$4668, and TiO$_1$ indices have strong correlations with binary fraction. % The two OIII indices are the most sensitive to the binary fraction, followed by four Balmer indices -- the two narrower central bandpass Balmer indices ($\sim$20\AA, F-definition) are more sensitive than the wider two ($\sim$40\AA, A-definition) and, lastly, the Ca4455, C$_2$4668, and TiO$_1$ indices.

We report the results of fitting Insight-HXMT data to the black hole X-ray binary MAXI J1348-430, which was discovered on January 26th, 2019, with the Gas Slit Camera (GSC) on-board MAXI. Several observations at the beginning of the first burst were selected, with a total of 10 spectra. From the residuals of fits using disk plus power law models, X-ray reflection signatures were clearly visible in some of these observations. We use the state-of the-art relxill series reflection model to fit six spectra with distinct reflection signatures and a joint fit to these spectra. In particular, we focus on the results for the black hole spin values. Assuming Rin = RISCO, the spin parameter is constrained to be 0.82+0.04-0.03 with 90% confidence level (statistical only).

We report the detection of millihertz quasi-periodic oscillations ($\mathrm{mHz}$ QPOs) from the low-mass X-ray binary XTE 1701$-$462. The discovery came from a search of the legacy data set of the Rossi X-ray Timing Explorer, in order to detect the periodic signals in all observations of sources exhibiting thermonuclear bursts. We found that $47$ out of $860$ observations of XTE 1701$-$462; covering the 2006--7 outburst exhibits signals with a significance above the detection threshold, which was determined separately for each observation via a Monte Carlo approach. We chose the four strongest candidates, each with maximum power exceeding $4\sigma$ of the simulated wavelet noise power distribution, to demonstrate the properties of the QPOs. The frequencies of the signals in the four observations are $\sim 3.5\;\text{to}\;5.6\; \mathrm{mHz}$, and the fractional R.M.S. amplitudes vary between $0.74 \pm 0.05\,\%$ and $3.54 \pm 0.04\,\%$. Although previously reported signals in other sources typically disappear immediately before a burst, we do not observe this behaviour in XTE 1701$-$462. Instead, we found that the QPOs and bursts occurred in separate accretion regimes. When the persistent luminosity dropped near the end of the outburst, the source showed bursts and no QPOs were detected, which is the behaviour predicted by theory for the transition from stable to unstable burning. On the basis of this new detection, we reassess the cases for identifying these $\mathrm{mHz}$ QPOs in this and other sources as arising from marginally stable burning.

Modified Newtonian Dynamics (MOND) is one of the most popular alternative theories of dark matter to explain the missing mass problem in galaxies. Although it remains controversial regarding MOND as a fundamental theory, MOND phenomenology has been shown to widely apply in different galaxies, which gives challenges to the standard $\Lambda$ cold dark matter model. In this article, we derive analytically the galactic rotation curve gradient in the MOND framework and present a rigorous analysis to examine the MOND phenomenology in our Galaxy. By assuming a benchmark baryonic disk density profile and two popular families of MOND interpolating functions, we show for the first time that the recent discovery of the declining Galactic rotation curve in the outer region ($R \approx 17-23$ kpc) can almost rule out the MOND phenomenology at more than $5\sigma$. This strongly supports some of the previous studies claiming that MOND is neither a fundamental theory nor a universal description of galactic properties.

Context. Several bright emission line regions have been observed in the central 100 parsecs of the active galaxy NGC 1068. Aims. We aim to determine the properties and ionization mechanism of three regions of NGC 1068: the nucleus (B) and two clouds located at 0.3" and 0.7" north of it (C and D). Methods. We combined SPHERE (0.95 - 1.65 um) and SINFONI (1.5 - 2.45 um) spectra for the three regions B, C, and D. We compared these spectra to several CLOUDY photoionization models and to the MAPPINGS III Library of Fast Radiative Shock Models. Results. The emission line spectra of the three regions are almost identical to each other and contribute to most of the emission line flux in the nuclear region. The emitting media contain multiple phases, the most luminous of which have temperatures ranging from 104.8 K to 106 K. Central photoionization models can reproduce some features of the spectra, but the fast radiative shock model provides the best fit to the data. Conclusions. The similarity between the three regions indicates that they belong to the same class of objects. Based on our comparisons, we conclude that they are shock regions located where the jet of the active galactic nucleus impacts massive molecular clouds.

We study gas and metal outflow from massive galaxies in protocluster regions at $z=3-9$ by using the results of the FOREVER22 simulation project. Our simulations contain massive haloes with $M_{\rm h} \gtrsim 10^{13}~\rm M_{\odot}$, showing high star formation rates of $> 100~\rm M_{\odot}~yr^{-1}$ and hosting supermassive black holes with $M_{\rm BH} \gtrsim 10^{8}~\rm M_{\odot}$. We show that the mass loading factor ($\eta_{\rm M}$) sensitively depends on the halo mass and it is $\eta_{\rm M} = 1.2~(9.2)$ for $M_{\rm h} = 10^{13}~(10^{11})~\rm M_{\odot}$. Once the halo mass exceeds $\sim 10^{12.5}~\rm M_{\odot}$, the outflow velocity of the gas rapidly decreases near a virial radius, and the gas returns to a galactic centre finally as a fountain flow. Also, the metal inflow and outflow rates sensitively depend on the halo mass and redshift. At $z=3$, the inflow rate becomes larger than the outflow one if $M_{\rm h} \gtrsim 10^{13.0}~\rm M_{\odot}$. Thus, we suggest that massive haloes cannot be efficient metal enrichment sources beyond virial radii that will be probed in future observations, e.g., studies of metal absorption lines with the Prime Focus Spectrograph on the Subaru telescope.

We present an investigation into the spectroscopic properties of non-equilibrium photoionization processes operating in a time-evolving mode. Through a quantitative comparison between equilibrium and time-evolving models, we find that the time-evolving model exhibits a broader distribution of charge states compared to the equilibrium model, accompanied by a slight shift in the peak ionization state depending on the source variability and gas density. The time-evolving code, tpho in SPEX, has been successfully employed to analyze the spectral properties of warm absorbers in the Seyfert galaxy NGC 3783. The incorporation of variability in the tpho model improves the fits of the time-integrated spectra, providing more accurate descriptions to the average charge states of several elements, in particular for Fe which is peaked around Fe XIX. The inferred densities and distances of the relevant X-ray absorber components are estimated to be approximately a few 1E11 per cubic meter and less than 1 pc, respectively. Furthermore, the updated fit suggests a potential scenario in which the observed absorbers are being expelled from the central AGN at the escape velocities. This implies that these absorbers might not play a significant role in the AGN feedback mechanism.

Growing evidence in support of a connection between Active Galactic Nuclei (AGN) activity and the Ram-Pressure Stripping (RPS) phenomenon has been found both observationally and theoretically in the past decades. In this work, we further explore the impact of RPS on the AGN activity by estimating the gas-phase metallicity of nuclear regions and the mass-metallicity relation of galaxies at $z \leq$ 0.07 and with stellar masses $\log {\rm M}_* / {\rm M}_\odot \geq 9.0 $, either experiencing RPS or not. To measure oxygen abundances, we exploit Integral Field Spectroscopy data from the GASP and MaNGA surveys, photoionization models generated with the code CLOUDY and the code Nebulabayes to compare models and observations. In particular, we build CLOUDY models to reproduce line ratios induced by photoionization from stars, AGN, or a contribution of both. We find that the distributions of metallicity and [O III]$\lambda$5007 luminosity of galaxies undergoing RPS are similar to the ones of undisturbed galaxies. Independently of the RPS, we do not find a correlation between stellar mass and AGN metallicity in the mass range $\log {\rm M}_* / {\rm M}_\odot \geq 10.4$, while for the star-forming galaxies we observe the well-known mass-metallicity relation (MZR) between $ 9.0 \leq \log \ {\rm M}_* /{\rm M}_\odot \leq 10.8$ with a scatter mainly driven by the star-formation rate (SFR) and a plateau around $\log {\rm M}_* / {\rm M}_\odot \sim 10.5$. The gas-phase metallicity in the nuclei of AGN hosts is enhanced with respect to those of SF galaxies by a factor of $\sim$ 0.05 dex regardless of the RPS.

Some supernovae (SNe), such as Type IIn SNe, are powered by collision of the SN ejecta with a dense circumstellar matter (CSM). Their emission spectra show characteristic line shapes of combined broad emission and narrow P-Cyg lines, which should closely relate to the CSM structure and the mass-loss mechanism that creates the dense CSM. We quantitatively investigate the relationship between the line shape and the CSM structure by Monte Carlo radiative transfer simulations, considering two representative cases of dense CSM formed by steady and eruptive mass loss. Comparing the H$\alpha$ emission between the two cases, we find that a narrow P-Cyg line appears in the eruptive case while it does not appear in the steady case, due to the difference in the velocity gradient in the dense CSM. We also reproduce the blue-shifted photon excess observed in some SNe IIn, which is formed by photon transport across the shock wave and find the relationship between the velocity of the shocked matter and the amount of the blue shift of the photon excess. We conclude that the presence or absence of narrow P-Cyg lines can distinguish the mass loss mechanism, and suggest high-resolution spectroscopic observations with $\lambda/ \Delta \lambda \gtrsim 10^4$ after the light curve peak for applying this diagnostic method.

Current simulations of air showers produced by ultra-high energy cosmic rays (UHECRs) do not satisfactorily describe recent experimental data, particularly when looking at the muonic shower component relative to the electromagnetic one. Discrepancies can be seen in both average values and on an individual shower-by-shower basis. It is thought that the muonic part of the air showers isn't accurately represented in simulations, despite various attempts to boost the number of muons within standard hadronic interaction physics. In this study, we investigate whether modifying the final state of events created with Sibyll~2.3d in air shower simulations can achieve a more consistent description of the muon content observed in experimental data. We create several scenarios where we separately increase the production of baryons, $\rho^0$, and strange particles to examine their impact on realistic air shower simulations. Our results suggest that these ad-hoc modifications can improve the simulations, providing a closer match to the observed muon content in air showers. One side-effect of the increased muon production in the considered model versions is a smaller difference in the predicted total muon numbers for proton and iron showers. However, more research is needed to find out whether any of these adjustments offers a realistic solution to the mismatches seen in data, and to identify the precise physical process causing these changes in the model. We hope that these modified model versions will also help to develop improved machine-learning analyses of air shower data and to estimate sys.{} uncertainties related to shortcomings of hadronic interaction models.

In this paper, we present an updated version of our model (KYNXiltr) which considers thermal reverberation of a standard Novikov-Thorne accretion disc illuminated by an X-ray point-like source. Previously, the model considered only two cases of black hole spins, and assumed a colour correction factor $f_{\rm col} = 2.4$. Now, we extend the model to any spin value and colour correction. In addition, we consider two scenarios of powering the X-ray corona, either via accretion, or external to the accretion disc. We use KYNXiltr to fit the observed time lags obtained from intense monitoring of four local Seyfert galaxies (NGC 5548, NGC 4395, Mrk 817, and Fairall 9). We consider various combinations of black hole spin, colour correction, corona height, and fraction of accretion power transferred to the corona. The model fits well the overall time-lags spectrum in these sources (for a large parameter space). For NGC 4593 only, we detect a significant excess of delays in the U-band. The contribution of the diffuse BLR emission in the time-lags spectrum of this source is significant. It is possible to reduce the large best-fitting parameter space by combining the results with additional information, such as the observed Eddington ratio and average X-ray luminosity. We also provide an update to the analytic expression provided by Kammoun et al., for an X-ray source that is not powered by the accretion process, which can be used for any value of colour correction, and for two values of the black hole spin (0 and 0.998).

Context. Due to the presence of magnetic fields, protostellar jets or outflows are a natural consequence of accretion onto protostars. They are expected to play an important role for star and protoplanetary disk formation. Aims. We aim to determine the influence of outflows on star and protoplanetary disk formation in star forming clumps. Methods. Using RAMSES, we perform the first magnetohydrodynamics calculation of massive star-forming clumps with ambipolar diffusion, radiative transfer including the radiative feedback of protostars and protostellar outflows while systematically resolving the disk scales. We compare it to a model without outflows. Results. We find that protostellar outflows have a significant impact on both star and disk formation. They provide additional turbulent and magnetic support to the clump, with typical velocities of a few 10 km/s, impact the disk temperatures, and reduce the accretion rate onto the protostars. While they promote a more numerous stellar population, we do not find that they control the mass scale of the stellar IMF. We find, however, that they have an influence on the high-mass end and the shape of the stellar IMF. Conclusions. Protostellar outflows appear to have a significant influence on both star and disk formation and should therefore be included in realistic simulations of star-forming environments.

The spatial metallicity distribution of star clusters in the Small Magellanic Cloud (SMC) has recently been found to correlate as a V-shaped function with the semi-major axis of an elliptical framework proposed to assume a projected galaxy flattening. We report results on the impact that the use of such a framework can produce on our understanding of the SMC formation and its chemical enrichment. We show that clusters with similar semi-major axes are placed at a very different distances from the SMC centre. The recently claimed bimodal metallicity distribution of clusters projected on the innermost SMC regions and the V-shaped metallicity gradient fade away when actual distances are used. Although a large dispersion prevails, clusters older than $\sim$ 1 Gyr exhibit a shallow metallicity gradient, caused by slightly different spatial distributions of clusters younger and older than $\sim$ 4 Gyr; the former being more centrally concentrated and having a mean metallicity ([Fe/H]) $\sim$ 0.15 dex more metal-rich than that of older clusters. This metallicity gradient does not show any dependence with the position angle, except for clusters placed beyond 11 kpc, which are located in the eastern side of the galaxy.

Prompt emissions from TeV blazars pair produce off the extragalactic background light and the highly energetic resulting pair beams then cascade through inverse Compton scattering to give rise to secondary gamma-rays. Such reprocessed cascade emission that can be associated with individual blazar sources has not been detected thus far. The absence of pair halos around these sources, along with the non-observation of isotropic gamma-ray background excess, seems to suggest that collective plasma effects, such as beam-plasma instabilities, can play a crucial role in alleviating this GeV-TeV tension by transferring the energy from the pair beams into the background plasma of the intergalactic medium (IGM). This has profound implications not only for TeV astrophysics, but also the strength of the intergalactic magnetic field and properties of dark matter (DM). A direct consequence of the instability losses and IGM heating is the modification of thermal history at late times, which suppresses structure formation particularly in baryonically underdense regions, potentially holding a clue towards resolving the small-scale crisis in cosmology. In a blazar-heated universe, the observation of dwarf galaxies and Lyman-$\alpha$ measurements present a favoured mass range for DM candidates such as light axion-like particles.

De-excitation gamma-ray lines, produced by nuclei colliding with protons, provide information about astrophysical environments where particles have kinetic energies of $10-100$ MeV per nucleon. In general, such environments can be categorized into two types: the interaction between non-thermal MeV cosmic rays and ambient gas, and the other is thermal plasma with a temperature above a few MeV. In this paper, we focus on the latter type and investigate the production of de-excitation gamma-ray lines in very hot thermal plasma, especially the dependence of the line profile on the plasma temperature. We have calculated the line profile of prompt gamma rays from $^{12}$C and $^{16}$O and found that when nuclei have a higher temperature than protons, gamma-ray line profiles can have a complex shape unique to each nucleus species. This is caused by anisotropic gamma-ray emission in the nucleus rest frame. We propose that the spectroscopy of nuclear de-excitation gamma-ray lines may enable to probe energy distribution in very hot astrophysical plasmas. This diagnostics can be a new and powerful technique to investigate the physical state of a two-temperature accretion flows onto a black hole, especially the energy distributions of the protons and nuclei, which are difficult to access for any other diagnostics.

The distribution and evolution of photospheric magnetic field in sunspots, plages and network, and variations in their relative flux content, play key roles in radial velocity (RV) fluctuations observed in Sun-as-a-star spectra. Differentiating and disentangling such magnetic contributions to RVs help in building models to account for stellar activity signals in high precision RV exoplanet searches. In this work, as earlier authors, we employ high-resolution images of the solar magnetic field and continuum intensities from SDO/HMI to understand the activity contributions to RVs from HARPS-N solar observations. Using well observed physical relationships between strengths and fluxes of photospheric magnetic fields, we show that the strong fields (spots, plages and network) and the weak internetwork fields leave distinguishing features in their contributions to the RV variability. We also find that the fill-factors and average unsigned magnetic fluxes of different features correlate differently with the RVs and hence warrant care in employing either of them as a proxy for RV variations. In addition, we examine disk averaged UV intensities at 1600 \r{A} and 1700 \r{A} wavelength bands imaged by SDO/AIA and their performances as proxies for variations in different magnetic features. We find that the UV intensities provide a better measure of contributions of plage fields to RVs than the Ca II H-K emission indices, especially during high activity levels when the latter tend to saturate.

Even if their detection is for now challenging, observation of small terrestrial planets will be easier in a near future thanks to continuous improvements of detection and characterisation instruments. In this quest, climate modeling is a key step to understand their characteristics, atmospheric composition and possible history. If a surface water reservoir is present on such a terrestrial planet, an increase in insolation may lead to a dramatic positive feedback induced by water evaporation: the runaway greenhouse. The resulting rise of global surface temperature leads to the evaporation of the entire water reservoir, separating two very different population of planets: 1) temperate planets with a surface water ocean and 2) hot planets with a puffed atmosphere dominated by water vapor. In this work we use a 3D General Circulation Model (GCM), the Generic-PCM, to study the runaway greenhouse transition, linking temperate and post-runaway states. Our simulations are made of two steps. First, assuming initially a liquid surface ocean, an evaporation phase which enriches the atmosphere in water vapor. Second, when the ocean is considered entirely evaporated, a dry transition phase for which the surface temperature increases dramatically. Finally, it converges on a hot and stable post-runaway state. By describing in detail the evolution of the climate during these two steps, we show a rapid transition of the cloud coverage and of the wind circulation from the troposphere to the stratosphere. By comparing our result to previous studies using 1D models, we discuss the effect of intrinsically 3D processes such as the global dynamics and the clouds, keys to understand the runaway greenhouse. We also explore the potential reversibility of the runaway greenhouse, limited by its radiative unbalance.

Coronal mass ejections (CMEs) are large eruptions from the Sun that propagate through the heliosphere after launch. Observational studies of these transient phenomena are usually based on 2D images of the Sun, corona, and heliosphere (remote-sensing data), as well as magnetic field, plasma, and particle samples along a 1D spacecraft trajectory (in-situ data). Given the large scales involved and the 3D nature of CMEs, such measurements are generally insufficient to build a comprehensive picture, especially in terms of local variations and overall geometry of the whole structure. This White Paper aims to address this issue by identifying the data sets and observational priorities that are needed to effectively advance our current understanding of the structure and evolution of CMEs, in both the remote-sensing and in-situ regimes. It also provides an outlook of possible missions and instruments that may yield significant improvements into the subject.

We present the median surface brightness profiles of diffuse Ly$\alpha$ haloes (LAHs) around star-forming galaxies by stacking 155 spectroscopically confirmed Ly$\alpha$ emitters (LAEs) at 3<z<4 in the MUSE Extremely Deep Field (MXDF), with median Ly$\alpha$ luminosity $\mathrm{L_{Ly\alpha} \approx 10^{41.1} erg\,s^{-1}}$. After correcting for a systematic surface brightness offset we identified in the datacube, we detect extended Ly$\alpha$ emission out to a distance of 270 kpc. The median Ly$\alpha$ surface brightness profile shows a power-law decrease in the inner 20 kpc, and a possible flattening trend at larger distance. This shape is similar for LAEs with different Ly$\alpha$ luminosities, but the normalisation of the surface brightness profile increases with luminosity. At distances larger than 50 kpc, we observe strong overlap of adjacent LAHs, and the Ly$\alpha$ surface brightness is dominated by the LAHs of nearby LAEs. We find no clear evidence of redshift evolution of the observed Ly$\alpha$ profiles when comparing with samples at 4<z<5 and 5<z<6. Our results are consistent with a scenario in which the inner 20 kpc of the LAH is powered by star formation in the central galaxy, while the LAH beyond a radius of 50 kpc is dominated by photons from surrounding galaxies.

Near-infrared bandpasses on spaceborne observatories diverge from their ground-based counterparts as they are free of atmospheric telluric absorption. Available transformations between respective filter systems in the literature rely on theoretical stellar atmospheres, which are known to have difficulties reproducing observed spectral energy distributions of cool giants. We present new transformations between the 2MASS $JHK_S$ and HST WFC3/IR F110W, F125W, & F160W photometric systems based on synthetic photometry of empirical stellar spectra from four spectral libraries. This sample comprises over 1000 individual stars, which together span nearly the full HR diagram and sample stellar populations from the solar neighborhood out to the Magellanic Clouds, covering a broad range of ages, metallicities, and other relevant stellar properties. In addition to global color-dependent transformations, we examine band-to-band differences for cool, luminous giant stars in particular, including multiple types of primary distance indicators.

Late-time spectra of Type Ia supernovae (SNe Ia) are important in clarifying the physics of their explosions, as they provide key clues to the inner structure of the exploding white dwarfs. We examined late-time optical spectra of 36 SNe Ia, including five from our own project (SNe 2019np, 2019ein, 2021hpr, 2021wuf, and 2022hrs), with phase coverage of $\sim 200$ to $\sim 400$ days after maximum light. At this late phase, the outer ejecta have become transparent and the features of inner iron-group elements emerge in the spectra. Based on multicomponent Gaussian fits and reasonable choices for the pseudocontinuum around Ni and Fe emission features, we get reliable estimates of the Ni to Fe ratio, which is sensitive to the explosion models of SNe Ia. Our results show that the majority (about 67%) of our SNe Ia are more consistent with the sub-Chandrasekhar-mass (i.e., double-detonation) model, although they could be affected by evolutionary or ionisation effects. Moreover, we find that the Si II $\lambda$6355 velocity measured around the time of maximum light tends to increase with the Ni to Fe ratio for the subsample with either redshifted or blueshifted nebular velocities, suggesting that progenitor metallicity might play an important role in accounting for the observed velocity diversity of SNe Ia.

High-precision exoplanet eclipse light curves, like those possible with JWST, enable flux and temperature mapping of exoplanet atmospheres. These eclipse maps will have unprecedented precision, providing an opportunity to constrain current theoretical predictions of exoplanet atmospheres. However, eclipse mapping has unavoidable mathematical limitations because many map patterns are unobservable. This ``null space'' has implications for making comparisons between predictions from general circulation models (GCMs) and the observed planet maps, and, thus, affects our understanding of the physical processes driving the observed maps. We describe the eclipse-mapping null space and show how GCM forward models can be transformed to their observable modes for more appropriate comparison with retrieved eclipse maps, demonstrated with applications to synthetic data of an ultra-hot Jupiter and a cloudy warm Jupiter under JWST-best-case- and extreme-precision observing scenarios. We show that the effects of the null space can be mitigated and manipulated through observational design, and JWST exposure times are short enough to not increase the size of the null space. Furthermore, we show the mathematical connection between the null space and the ``eigenmapping'' method, demonstrating how eigenmaps can be used to understand the null space in a model-independent way. We leverage this connection to incorporate null-space uncertainties in retrieved maps, which increases the uncertainties to now encompass the ground truth for synthetic data. The comparisons between observed maps and forward models that are enabled by this work, and the improved eclipse-mapping uncertainties, will be critical to our interpretation of multidimensional aspects of exoplanets in the JWST era.

Dwarf planets Eris and Makemake have surfaces bearing methane ice of unknown origin. D/H ratios were recently determined from James Webb Space Telescope (JWST) observations of Eris and Makemake (Grundy et al., submitted), giving us new clues to decipher the origin of methane. Here, we develop geochemical models to test if the origin of methane could be primordial, derived from CO$_2$ or CO ("abiotic"), or sourced by organics ("thermogenic"). We find that primordial methane is inconsistent with the observational data, whereas both abiotic and thermogenic methane can have D/H ratios that overlap the observed ranges. This suggests that Eris and Makemake either never acquired a significant amount of methane during their formation, or their original inventories were removed and then replaced by a source of internally produced methane. Because producing abiotic or thermogenic methane likely requires temperatures in excess of ~150{\deg}C, we infer that Eris and Makemake have rocky cores that underwent substantial radiogenic heating. Their cores may still be warm/hot enough to produce methane. This heating could have driven hydrothermal circulation at the bottom of an ice-covered ocean to generate abiotic methane, and/or metamorphic reactions involving accreted organic matter could have occurred in response to heating in the deeper interior, generating thermogenic methane. Additional analyses of thermal evolution model results and predictions from modeling of D-H exchange in the solar nebula support our findings of elevated subsurface temperatures and a lack of primordial methane on Eris and Makemake. It remains an open question whether their D/H ratios may have evolved subsequent to methane outgassing. Recommendations are given for future activities to further test proposed scenarios of abiotic and thermogenic methane production on Eris and Makemake, and to explore these worlds up close.

The search for habitable environments and biomarkers in exoplanetary atmospheres is the holy grail of exoplanet science. The detection of atmospheric signatures of habitable Earth-like exoplanets is challenging due to their small planet-star size contrast and thin atmospheres with high mean molecular weight. Recently, a new class of habitable exoplanets, called Hycean worlds, has been proposed, defined as temperate ocean-covered worlds with H2-rich atmospheres. Their large sizes and extended atmospheres, compared to rocky planets of the same mass, make Hycean worlds significantly more accessible to atmospheric spectroscopy with the JWST. Here we report a transmission spectrum of the candidate Hycean world, K2-18 b, observed with the JWST NIRISS and NIRSpec instruments in the 0.9-5.2 $\mu$m range. The spectrum reveals strong detections of methane (CH4) and carbon dioxide (CO2) at 5$\sigma$ and 3$\sigma$ confidence, respectively, with high volume mixing ratios of ~1% each in a H2-rich atmosphere. The abundant CH4 and CO2 along with the non-detection of ammonia (NH3) are consistent with chemical predictions for an ocean under a temperate H2-rich atmosphere on K2-18 b. The spectrum also suggests potential signs of dimethyl sulfide (DMS), which has been predicted to be an observable biomarker in Hycean worlds, motivating considerations of possible biological activity on the planet. The detection of CH4 resolves the long-standing missing methane problem for temperate exoplanets and the degeneracy in the atmospheric composition of K2-18 b from previous observations. We discuss possible implications of the findings, open questions, and future observations to explore this new regime in the search for life elsewhere.

We report multi-frequency VLBI studies of the sub-parsec scale structure of the two-sided jet in the nearby radio galaxy NGC 4261. Our analyses include new observations using the Source Frequency Phase Referencing technique with the Very Long Baseline Array at 44 and 88 GHz, as well as archival data at 15 and 43 GHz. Our results show an extended double-sided structure at 43/44 GHz and provide a clear image of the nuclear region at 88 GHz, showing a core size of $\sim$0.09 mas and a brightness temperature of $\sim1.3\times10^{9}$ K. Proper motions are measured for the first time in the two-sided jet, with apparent speeds ranging from $0.31\pm0.14\,c$ to $0.59\pm0.40\,c$ in the approaching jet and $0.32\pm0.14\,c$ in the receding jet. The jet-to-counter-jet brightness ratio allows us to constrain the viewing angle to between $\sim54^{\circ}$ and $84^{\circ}$ and the intrinsic speed to between $\sim0.30\,c$ and $0.55\,c$. We confirm the parabolic shape of the upstream jet on both sides of the central engine, with a power-law index of $0.56\pm0.07$. Notably, the jet collimation is found to be already completed at sub-parsec scales, with a transition location of about 0.61 pc, which is significantly smaller than the Bondi radius of 99.2 pc. This behavior can be interpreted as the initial confinement of the jet by external pressure from either the geometrically thick, optically thin advection-dominated accretion flows (ADAF) or the disk wind launched from it. Alternatively, the shape transition may also be explained by the internal flow transition from a magnetically dominated to a particle-dominated regime.

Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4-7 per cent).

Horizon-scale observations of the jetted active galactic nucleus M87 are compared with simulations spanning a broad range of dissipation mechanisms and plasma content in three-dimensional general relativistic flows around spinning black holes. Observations of synchrotron radiation from radio to X-ray frequencies can be compared with simulations by adding prescriptions specifying the relativistic electron-plus-positron distribution function and associated radiative transfer coefficients. A suite of time-varying simulations with various spins and plasma magnetizations is chosen to represent distinct possibilities for the M87 jet/accretion flow/black hole (JAB) system. We then input turbulent heating and equipartition-based emission prescriptions (and piecewise combinations thereof) in the time-dependent 3D simulations, in which jet morphology, polarization and variation are "observed" and compared with real observations so as to try to infer the rules that govern the polarized emissivity. The models in this paper support a magnetically arrested disk (MAD) with several possible spin/emission model combinations supplying the jet in M87, whose inner jet and black hole shadow have been observed down to the photon ring at 230 GHz by the Event Horizon Telescope (EHT). We also show that some MAD cases that are dominated by intrinsic circular polarization have near-linear V/I dependence on unpaired electron or positron content while SANE polarization exhibits markedly greater positron-dependent Faraday effects -- future probes of the SANE/MAD dichotomy and plasma content with the EHT. This is the second work in a series also applying the "observing" simulations methodology to near-horizon regions of supermassive black holes in Sgr A* and 3C 279.

One of the most exciting and pressing issues in cosmology today is the discrepancy between some measurements of the local Hubble constant and other values of the expansion rate inferred from the cosmic microwave background (CMB) radiation. Resolving these differences holds the potential for the discovery of new physics beyond the standard model of cosmology: Lambda Cold Dark Matter (LCDM), a successful model that has been in place for more than 20 years. Given both the fundamental significance of this outstanding discrepancy, and the many-decades-long effort to increase the accuracy of the extragalactic distance scale, it is critical to demonstrate that the local measurements are convincingly free from residual systematic errors. We review the progress over the past quarter century in measurements of the local value of the Hubble constant, and discuss remaining challenges. Particularly exciting are new data from the James Webb Space Telescope (JWST), that is delivering high-resolution near-infrared imaging data to both test for and to address directly several of the systematic uncertainties that have historically limited the accuracy of extragalactic distance scale. We present an overview of our new JWST program to observe Cepheids, TRGB and JAGB stars. For the first galaxy in our program, NGC 7250, the high-resolution JWST images demonstrate that many of the Cepheids observed with the Hubble Space Telescope (HST) are significantly crowded by nearby neighbors. Avoiding the more significantly crowded variables, the scatter in the JWST near-infrared (NIR) Cepheid period-luminosity relation is decreased by a factor of two compared to those from HST, illustrating the power of JWST for improvements to local measurements of Ho. Ultimately, these data will either confirm the standard model, or provide robust evidence for the inclusion of additional new physics.

The conditions required for massive star formation are debated, particularly whether massive stars must form in conjunction with massive clusters. Some authors have advanced the view that stars of any mass (below the total cluster mass) can form in clusters of any mass with some probability (random sampling). Others pointed out that the scatter in the determinations of the most massive star mass for a given cluster mass was consistent with the measurement error, such that the mass of the most massive star was determined by the total cluster mass (optimal sampling). Here we investigate the relation between cluster mass (M\textsubscript{ecl}) and the maximum stellar mass (M\textsubscript{max}) using a suite of SPH simulations. Varying cloud mass and turbulence random seed results in a range of cluster masses which we compare with their respective maximum star masses. We find that more massive clusters will have, on average, higher mass stars with this trend being steeper at lower cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.31}$ for $M\textsubscript{ecl}<500M\,_{\odot}$) and flattening at higher cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.11}$ for $M\textsubscript{ecl}>500M\,_{\odot}$). This rules out purely stochastic star formation in our simulations. Significant scatter in the maximum masses with identical initial conditions also rules out the possibility that the relation is purely deterministic (that is that a given cluster mass will result in a specific maximum stellar mass). In conclusion our simulations disagree with both random and optimal sampling of the initial mass function.

The auto-correlation function of the Lyman-$\alpha$ (Ly$\alpha$) forest flux from high-z quasars can statistically probe all scales of the intergalactic medium (IGM) just after the epoch of reionization. The thermal state of the IGM, which is determined by the physics of reionization, sets the amount of small-scale power seen in the \lya forest. To study the sensitivity of the auto-correlation function to the thermal state of the IGM, we compute the auto-correlation function from cosmological hydrodynamical simulations with semi-numerical models of the thermal state of the IGM. We create mock data sets of 20 quasars to forecast constraints on $T_0$ and $\gamma$, which characterize a tight temperature-density relation in the IGM, at $5.4 \leq z \leq 6$. At $z = 5.4$ we find that an ideal data set constrains $T_0$ to 29\% and $\gamma$ to 9\%. In addition, we investigate four realistic reionization scenarios that combine temperature and ultra-violet background (UVB) fluctuations at $z = 5.8$. We find that, when using mock data generated from a model that includes temperature and UVB fluctuations, we can rule out a model with no temperature or UVB fluctuations at $>1\sigma$ level 50.5\% of the time.

The evolution of protoplanetary disks is regulated by an interplay of several processes, either internal to the system or related to the environment. As most of the stars and planets have formed in massive stellar clusters, studying the effects of UV radiation on disk evolution is of paramount importance. Here we test the impact of external photoevaporation on the evolution of disks in the $\sigma$ Orionis cluster by conducting the first combined large-scale UV to IR spectroscopic and mm-continuum survey of this region. We study a sample of 50 targets located at increasing distances from the central, OB system $\sigma$ Ori. We combine new VLT/X-Shooter spectra with new and previously published ALMA measurements of disk dust and gas fluxes and masses. We confirm the previously found decrease of $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. This is particularly evident when considering the disks around the more massive stars ($\ge$ 0.4 $M_{\odot}$), where those located in the inner part ($<$ 0.5 pc) have $M_{\rm dust}$ about an order of magnitude lower than the more distant ones. About half of the sample is located in the region of the $\dot{M}_{\rm acc}$ vs $M_{\rm disk}$ expected by models of external photoevaporation, namely showing shorter disk lifetimes. These are observed for all targets with projected separation from $\sigma$ Ori $<$ 0.5 pc, proving that the presence of a massive stellar system affects disk evolution. External photoevaporation is a viable mechanism to explain the observed shorter disk lifetimes and lower $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. Follow-up observations of the low stellar mass targets are crucial to confirm the dependence of the external photoevaporation process with stellar host mass. This work confirms that the effects of external photoevaporation are significant down to impinging radiation as low as $\sim 10^{4}$ G$_0$.

Future constraints of cosmological parameters from Type Ia supernovae (SNe Ia) will depend on the use of photometric samples, those samples without spectroscopic measurements of the SNe Ia. There is a growing number of analyses that show that photometric samples can be utilised for precision cosmological studies with minimal systematic uncertainties. To investigate this claim, we perform the first analysis that combines two separate photometric samples, SDSS and Pan-STARRS, without including a low-redshift anchor. We evaluate the consistency of the cosmological parameters from these two samples and find they are consistent with each other to under $1\sigma$. From the combined sample, named Amalgame, we measure $\Omega_M = 0.328 \pm 0.024$ with SN alone in a flat $\Lambda$CDM model, and $\Omega_M = 0.330 \pm 0.018$ and $w = -1.016^{+0.055}_{-0.058}$ when combining with a Planck data prior and a flat $w$CDM model. These results are consistent with constraints from the Pantheon+ analysis of only spectroscopically confirmed SNe Ia, and show that there are no significant impediments to analyses of purely photometric samples of SNe Ia.

Despite a growing sample of precisely measured stellar rotation periods and ages, the strength of magnetic braking and the degree of departure from standard (Skumanich-like) spindown have remained persistent questions, particularly for stars more evolved than the Sun. Rotation periods can be measured for stars older than the Sun by leveraging asteroseismology, enabling models to be tested against a larger sample of old field stars. Because asteroseismic measurements of rotation do not depend on starspot modulation, they avoid potential biases introduced by the need for a stellar dynamo to drive starspot production. Using a neural network trained on a grid of stellar evolution models and a hierarchical model-fitting approach, we constrain the onset of weakened magnetic braking. We find that a sample of stars with asteroseismically-measured rotation periods and ages is consistent with models that depart from standard spindown prior to reaching the evolutionary stage of the Sun. We test our approach using neural networks trained on model grids produced by separate stellar evolution codes with differing physical assumptions and find that the choices of grid physics can influence the inferred properties of the braking law. We identify the normalized critical Rossby number ${\rm Ro}_{\rm crit}/{\rm Ro}_\odot = 0.91\pm0.03$ as the threshold for the departure from standard rotational evolution. This suggests that weakened magnetic braking poses challenges to gyrochronology for roughly half of the main sequence lifetime of sun-like stars.

This article aims at comparing gravitational wave memory effect in a Schwarzschild spacetime with that of other compact objects with static and spherically symmetric spacetime, with the purpose of proposing a procedure for differentiating between various compact object geometries. We do this by considering the relative evolution of two nearby test geodesics with in different backgrounds in the presence and absence of a gravitational wave pulse and comparing them. Memory effect due to a gravitational wave would ensure that there is a permanent effect on each spacetime and the corresponding geodesic evolution, being metric dependent, would display distinct results in each case. For a complete picture, we have considered both displacement and velocity memory effect in each geometry.

Most tests of general relativity with gravitational-wave observations rely on inferring the degree to which a signal deviates from general relativity in conjunction with the astrophysical parameters of its source, such as the component masses and spins of a compact binary. Due to features of the signal, measurements of these deviations are often highly correlated with the properties of astrophysical sources. As a consequence, prior assumptions about astrophysical parameters will generally affect the inferred magnitude of the deviations. Incorporating information about the underlying astrophysical population is necessary to avoid biases in the inference of deviations from general relativity. Current tests assume that the astrophysical population follows an unrealistic fiducial prior chosen to ease sampling of the posterior -- for example, a prior flat in component masses -- which is is inconsistent with both astrophysical expectations and the distribution inferred from observations. We propose a framework for fortifying tests of general relativity by simultaneously inferring the astrophysical population using a catalog of detections. Although this method applies broadly, we demonstrate it concretely on massive graviton constraints and parameterized tests of deviations to the post-Newtonian phase coefficients. Using observations from LIGO-Virgo-KAGRA's third observing run, we show that concurrent inference of the astrophysical distribution strengthens constraints and improves overall consistency with general relativity. We provide updated constraints on deviations from the theory, finding that, upon modeling the astrophysical population, the 90\%-credible upper limit on the mass of the graviton improves by $25\%$ to $m_g \leq 9.6 \times 10^{-24}\, \mathrm{eV}/c^2$ and the inferred population-level post-Newtonian deviations move ${\sim} 0.4 \sigma$ closer to zero.

Solar flare prediction is a central problem in space weather forecasting and recent developments in machine learning and deep learning accelerated the adoption of complex models for data-driven solar flare forecasting. In this work, we developed an attention-based deep learning model as an improvement over the standard convolutional neural network (CNN) pipeline to perform full-disk binary flare predictions for the occurrence of $\geq$M1.0-class flares within the next 24 hours. For this task, we collected compressed images created from full-disk line-of-sight (LoS) magnetograms. We used data-augmented oversampling to address the class imbalance issue and used true skill statistic (TSS) and Heidke skill score (HSS) as the evaluation metrics. Furthermore, we interpreted our model by overlaying attention maps on input magnetograms and visualized the important regions focused on by the model that led to the eventual decision. The significant findings of this study are: (i) We successfully implemented an attention-based full-disk flare predictor ready for operational forecasting where the candidate model achieves an average TSS=0.54$\pm$0.03 and HSS=0.37$\pm$0.07. (ii) we demonstrated that our full-disk model can learn conspicuous features corresponding to active regions from full-disk magnetogram images, and (iii) our experimental evaluation suggests that our model can predict near-limb flares with adept skill and the predictions are based on relevant active regions (ARs) or AR characteristics from full-disk magnetograms.

Finite-energy particles in free fall can collide with diverging center-of-mass energy near rapidly rotating black holes. What are the most salient observational signatures of this remarkable geometric effect? Here we revisit the problem from the standpoint of the near-horizon extreme Kerr geometry, where these collisions naturally take place. It is shown that the ingoing particle kinematics admits a simple, universal form. Given a scattering cross section, determination of emission properties is reduced to evaluation of particular integrals on the sky of a near-horizon orbiting particle. We subsequently apply this scheme to the example of single-photon bremsstrahlung, substantiating past results which indicate that ejected particles are observable, but their energies are bounded by the rest masses of the colliding particles. Our framework is readily applicable for any scattering process.

The Southern Wide-field Gamma-ray Observatory (SWGO) is an R&D project to plan and design the next observatory to detect gamma rays in the Southern hemisphere. The experiment, planned to be placed at an altitude greater than 4400 m, is primarily based on water Cherenkov detectors units and is expected to measure gamma rays from a few hundred GeV up to the PeV scale. SWGO will complement CTA and the existing ground-based particle detectors of the Northern Hemisphere, namely HAWC and LHAASO, having a rich science programme. The collaboration is highly invested in evaluating different detector and array configurations, prototyping, and site search. In this presentation, I shall present an overview of the project's activities, achievements and future plans.

We consider F-term hybrid inflation and supersymmetry breaking in the context of a model which largely respects a global U(1) R symmetry. The Kaehler potential parameterizes the Kaehler manifold with an enhanced U(1)x(SU(1,1)/U(1)) symmetry, where the scalar curvature of the second factor is determined by the achievement of a supersymmetry-breaking de Sitter vacuum without ugly tuning. The magnitude of the emergent soft tadpole term for the inflaton can be adjusted in the range (1.2-460) TeV -- increasing with the dimensionality of the representation of the waterfall fields -- so that the inflationary observables are in agreement with the observational requirements. The mass scale of the supersymmetric partners turns out to lie in the region (0.09-253) PeV which is compatible with high-scale supersymmetry and the results of LHC on the Higgs boson mass. The mu parameter can be generated by conveniently applying the Giudice-Masiero mechanism and assures the out-of-equilibrium decay of the R saxion at a low reheat temperature Trh<~163 GeV.

The discovery of universe's late-time acceleration and dark energy has overseen a great deal of research into cosmological singularities and particularly future singularities. Perhaps the most extreme of such singlarities is the big rip, which has propelled a lot of work into ways of moderating it or seeking out alternatives to it and two such alternatives to the big rip are the Little rip and Pseudo rip. Another possibility to consider the far future of the universe is through bounce cosmologies, which presents its own interesting ideas. So in this work we investigate the Little rip, Pseudo rip and Bounce cosmology in non-standard cosmological backgrounds with a generalized equation of state in the presence of a viscous fluid. In particular we discuss about Chern-Simons cosmology and the RS-II Braneworld and discuss how the exotic and non-conventional nature of gravity in such cosmologies affect universal evolution in these scenarios. We find out that there are very significant differences in the behaviour of such cosmic scenarios in these universes in comparison to how they appear in the simple general relativistic universe.

The LIGO-Virgo analysis of the signals from compact binary mergers observed so far have assumed in-vacuum isolated binary systems, neglecting the potential presence of astrophysical environments. Non-trivial environments may alter gravitational-wave emission, leaving imprints that can be observable via a characteristic dephasing of the emitted signal with respect to the vacuum scenario. We present here the first investigation of environmental effects on the events of the first gravitational-wave catalog (GWTC-1) by LIGO-Virgo. We include the effects of accretion and dynamical friction through a post-Newtonian deformation of the inspiral part of the waveform relative to the vacuum one. We find no evidence for the presence of environmental effects in GWTC-1. Most of the events decisively exclude the scenario of dynamical fragmentation of massive stars as their possible formation channel. Our analysis of GW170817 results in the upper bound on the medium density of $\lesssim 21\: \text{g/cm}^3$. We find that environmental effects can substantially bias the recovered parameters in the vacuum model, even when they are not detectable. Our results forecast that the future 2030s detectors Einstein Telescope and B-DECIGO will be able to probe the environmental effects of accretion disk and superradiant boson clouds on compact binaries.

In the efforts of searching for dark matter, gravitational wave interferometers have been recently proposed as a promising probe. These highly sensitive instruments are potentially be able to detect the interactions of dark matter with the detectors. In this work, we explored the possibilities of discovering topological dark matter with LIGO detectors. We analyzed domain walls consisting of axion-like dark matter passing through Earth, leaving traces in multiple detectors simultaneously. Considering dark matter interactions with the light in the interferometer and with the beamsplitter, we performed the first analysis of the topological dark matter with the gravitational-wave strain data. We examined whether astrophysically unexpected triggers could be explained by domain wall passages. We found that all of the binary black hole mergers we analyzed favored the binary black hole merger hypothesis rather than the domain wall hypothesis, with the closest being GW190521. Moreover we found that some of topological dark matter signals can be caught by binary black hole searches. Finally, we found that glitches in the data can inevitably limit the dark matter searches for certain parameters. These results are expected to guide the future searches and analyses.

A massive spin-2 field can grow unstably around a black hole, potentially probing the existence of such fields. In this work, we use time-domain evolutions to study such instabilities. Considering the linear regime by solving the equations generically governing a massive tensor field on the background of a Kerr black hole, we find that black hole spin significantly increases the mass range, and the growth rate, of the axisymmetric (azimuthal number $m=0$) instability, which takes the form of the Gregory-Laflamme black string instability for zero spin. We also consider the superradiant unstable modes with $1 \leq m \leq 3$, extending previous results to higher spin-2 masses, black hole spins, and azimuthal numbers. We find that the superradiant modes grow slower than the $m=0$ modes, except for a narrow range of high spins and masses, with $m=1$ and 2 requiring a dimensionless black hole spin of $a_{\rm BH}\gtrsim 0.95$ to be dominant. Thus, in most of the parameter space, the backreaction of the $m=0$ instability must be taken into account when using black holes to constrain massive spin-2 fields. As a simple model of this, we consider nonlinear evolutions in quadratic gravity, in particular Einstein-Weyl gravity. We find that, depending on the initial perturbation, the black hole may approach zero mass with the curvature blowing up in finite time, or can saturate at a larger mass with a surrounding cloud of the ghost spin-2 field.

We investigate the cosmological collider (CC) signal arising from the tree-level exchange of a scalar spectator particle with a non-Bunch Davies (BD) initial state. We decompose the inflaton correlators into seed integrals, which we compute analytically by solving the bootstrap equations. We show that the non-BD initial state eliminates the Hubble scale Boltzmann suppression $e^{-\pi m /H}$ that usually affects the CC signal. Consequently, in this scenario, the CC can probe an energy scale much higher than the inflationary Hubble scale $H$.

Despite its elegance, the theory of General Relativity is subject to experimental, observational, and theoretical scrutiny to arrive at tighter constraints or an alternative, more preferred theory. In alternative gravity theories, the macroscopic properties of neutron stars, such as mass, radius, tidal deformability, etc. are modified. This creates a degeneracy between the uncertainties in the equation of state (EoS) and gravity since assuming a different EoS can be mimicked by changing to a different theory of gravity. We formulate a hierarchical Bayesian framework to simultaneously infer the EoS and gravity parameters by combining multiple astrophysical observations. We test this framework for a particular 4D Horndeski scalar-tensor theory originating from higher-dimensional Einstein-Gauss-Bonnet gravity and a set of 20 realistic EoS and place improved constraints on the coupling constant of the theory with current observations. Assuming a large number of observations with upgraded or third-generation detectors, we find that the $A+$ upgrade could place interesting bounds on the coupling constant of the theory, whereas with the LIGO Voyager upgrade or the third-generation detectors (Einstein Telescope and Cosmic Explorer), the degeneracy between EoS and gravity could be resolved with high confidence, even for small deviations from GR.

Gravitational wave (GW) predictions of cosmological phase transitions are almost invariably evaluated at either the nucleation or percolation temperature. We investigate the effect of the transition temperature choice on GW predictions, for phase transitions with weak, intermediate and strong supercooling. We find that the peak amplitude of the GW signal varies by a factor of a few for weakly supercooled phase transitions, and by an order of magnitude for strongly supercooled phase transitions. The variation in amplitude for even weakly supercooled phase transitions can be several orders of magnitude if one uses the mean bubble separation, while the variation is milder if one uses the mean bubble radius instead. We also investigate the impact of various approximations used in GW predictions. Many of these approximations introduce at least a 10% error in the GW signal, with others introducing an error of over an order of magnitude.

From investigating molecular vibrations to observing galaxies, terahertz technology has found extensive applications in research and development over the past three decades. Terahertz time-domain spectroscopy and imaging have experienced significant growth and now dominate spectral observations ranging from 0.1 to 10 THz. However, the lack of standardised protocols for data processing, dissemination, and archiving poses challenges in collaborating and sharing terahertz data between research groups. To tackle these challenges, we present the dotTHz project, which introduces a standardised terahertz data format and the associated open-source tools for processing and interpretation of dotTHz files. The dotTHz project aims to facilitate seamless data processing and analysis by providing a common framework. All software components are released under the MIT licence through GitHub repositories to encourage widespread adoption, modification, and collaboration. We invite the terahertz community to actively contribute to the dotTHz project, fostering the development of additional tools that encompass a greater breadth and depth of functionality. By working together, we can establish a comprehensive suite of resources that benefit the entire terahertz community.

We study the interaction of several types of static straight cosmic strings, including local strings, global strings, and bosonic superconducting strings with and without magnetic currents. First, we evaluate the interaction energy of two widely separated cosmic strings using the point source formalism and show that the most dominant contribution to the interaction energy comes from the excitation of the lightest mediator particles in a underlying theory. The interaction energy at arbitrary separation distances is then analyzed numerically by the gradient flow method. It turns out that an additional scalar field introduced in the bosonic superconducting string becomes an additional source of attraction. For such a bosonic superconducting string, we find that a string with two winding numbers is energetically favorable compared to two strings with a single winding number in a certain parameter region. Our analysis reveals that a phase structure of bosonic superconducting strings is richer than that of local and global strings and that the formation of bound states at intersections of bosonic superconducting strings is favored.

We investigate the linear and nonlinear evolution of the ion-acoustic instability in a collisionless plasma via two-dimensional (2D2V) Vlasov-Poisson numerical simulations. We initialize the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding super-critical initial conditions that are physically unrealizable. The nonlinear evolution of ion-acoustic turbulence (IAT) is characterized in detail, including the particles' distribution functions, particle heating, (two-dimensional) wave spectrum, and the resulting anomalous resistivity. An important result is that no steady saturated nonlinear state is ever reached in our simulations: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived. Electron-acoustic waves (EAWs) are triggered during the late nonlinear evolution of the system, caused by strong modifications to the particle distribution induced by IAT.

We provide accurate universal relations that allow to estimate the moment of inertia $I$ and the ratio of kinetic to gravitational binding energy $T/W$ of uniformly rotating neutron stars from the knowledge of mass, radius, and moment of inertia of an associated non-rotating neutron star. Based on these, several other fluid quantities can be estimated as well. Astrophysical neutron stars rotate to varying degrees and although rotational effects may be neglected in some cases, not modeling them will inevitably introduce bias when performing parameter estimation. This is especially important for future, high-precision measurements coming from electromagnetic and gravitational wave observations. The proposed universal relations facilitate computationally cheap EOS inference codes that permit the inclusion of observations of rotating neutron stars. To demonstrate this, we deploy them into a recent Bayesian framework for equation of state parameter estimation that is now valid for arbitrary, uniform rotation. Our inference results are robust up to around percent level precision for the generated neutron star observations, consisting of the mass, equatorial radius, rotation rate, as well as co- and counter-rotating $f$-mode frequencies, that enter the framework as data.