Papers for Tuesday, Sep 10 2024

We propose a new T, Y dwarf search model using machine learning (ML), called the "BROWn Dwarf Image Explorer (BROWDIE). Brown dwarfs (BD) are estimated to make up 25 percent of all celestial objects in the Galaxy, yet only a small number have been thoroughly studied. Homogeneous and complete samples of BDs are essential to advance the studies. However, due to their faintness, conducting spectral studies of BDs can be challenging. T\&Y brown dwarfs, a redder and fainter subclass of BDs, are even harder to detect. As a result, only a few T\&Y dwarfs have been extensively studied. Numerous attempts, including ML using various color band observations, have been made to identify BDs based on their colors. However, those models often require a large number of color observations, which can be a limitation. This study implemented an ML model by utilizing data from the J, H, and K photometry of the UKIRT Infrared Deep Sky Surveys (UKIDSS) simultaneously, effectively distinguishing celestial objects similar to BDs. The BROWDIE model was trained using UKIDSS data, and based on the model, 118 T dwarfs and 14 Y dwarfs were found in the UKIDSS DR11PLUS LAS L4 zone.

The ability to accurately discern active massive black holes (BHs) in nearby dwarf galaxies is paramount to understanding the origins and processes of "seed" BHs in the early Universe. We present Chandra X-ray Observatory observations of a sample of three local dwarf galaxies (M$_{*}$ $\leqslant 3 \times 10^{9}$ M$_\odot$, z $\leqslant$ 0.15) previously identified as candidates for hosting active galactic nuclei (AGN). The galaxies were selected from the NASA-Sloan Atlas (NSA) with spatially coincident X-ray detections in the eROSITA Final Equatorial Depth Survey (eFEDS). Our new Chandra data reveal three X-ray point sources in two of the target galaxies with luminosities between log(L$_{\rm \text{2-10 keV}}$ [erg s$^{-1}$]) = 39.1 and 40.4. Our results support the presence of an AGN in these two galaxies and a ULX in one of them. For the AGNs, we estimate BH masses of $M_{\rm BH} \sim 10^{5-6} M_\odot$ and Eddington ratios on the order of $\sim 10^{-3}$.

Globular clusters (GCs) are sensitive tracers of galaxy assembly histories but interpreting the information they encode is challenging because mergers are thought to promote both the formation and disruption of GCs. We use simulations with controlled merger histories to examine the influence of merger mass ratio on the GC population of a present-day $L^\ast$ galaxy, using the genetic modification technique to adjust the initial conditions of a galaxy that experiences major mergers at $z = 1.7$ and $z = 0.77$ (ORGANIC case), so the later merger has twice its original mass ratio (ENHANCED case), or is prevented from occurring (SUPPRESSED case). We evolve the three realizations with E-MOSAICS, which couples sub-grid star cluster formation and evolution models to the EAGLE galaxy formation model. Relative to the ORGANIC case, the mass of surviving GCs is elevated (reduced) in the ENHANCED (SUPPRESSED) case, indicating that major mergers promote a net boost to the GC population. The boost is clearly quantified by the GC specific mass, $S_{\rm M}$, because it is sensitive to the number of the most massive GCs, whose long characteristic disruption timescales enable them to survive their hostile natal environments. In contrast, the specific frequency, $T_{\rm N}$, is insensitive to assembly history because it primarily traces low-mass GCs that tend to be disrupted soon after their formation. The promotion of GC formation and disruption by major mergers imprints a lasting and potentially observable signature: an elevated mass fraction of field stars in the galaxy's stellar halo that were born in star clusters.

If primordial black holes (PBHs) of asteroidal mass ($M_{\rm PBH}\in[10^{17},10^{23}]$ g) make up the entire dark matter they could be detectable through their gravitational influence in the solar system. In this work we study the perturbations that PBHs induce on the orbits of planets. Detailed numerical simulations of the solar system, embedded in a halo of primordial black holes are performed. We show that the perturbations are too small to be directly detectable with current data, challenging recent results that have ruled out PBHs as a dark matter candidate. Using the Earth-Mars system as an observational probe, we estimate that an improvement in the measurement accuracy by more than an order of magnitude is required to detect the gravitational influence of PBHs in the solar system in the foreseeable future.

The wealth of high-quality observational data from the epoch of reionization that will become available in the next decade motivates further development of modeling techniques for their interpretation. Among the key challenges in modeling reionization are (1) its multi-scale nature, (2) the computational demands of solving the radiative transfer (RT) equation, and (3) the large size of reionization's parameter space. In this paper, we present and validate a new RT code designed to confront these challenges. FlexRT (Flexible Radiative Transfer) combines adaptive ray tracing with a highly flexible treatment of the intergalactic ionizing opacity. This gives the user control over how the intergalactic medium (IGM) is modeled, and provides a way to reduce the computational cost of a FlexRT simulation by orders of magnitude while still accounting for small-scale IGM physics. Alternatively, the user may increase the angular and spatial resolution of the algorithm to run a more traditional reionization simulation. FlexRT has already been used in several contexts, including simulations of the Lyman-$\alpha$ forest of high-$z$ quasars, the redshifted 21cm signal from reionization, as well as in higher resolution reionization simulations in smaller volumes. In this work, we motivate and describe the code, and validate it against a set of standard test problems from the Cosmological Radiative Transfer Comparison Project. We find that FlexRT is in broad agreement with a number of existing RT codes in all of these tests. Lastly, we compare FlexRT to an existing adaptive ray tracing code to validate FlexRT in a cosmological reionization simulation.

Due to high-cadence automated surveys, we can now detect and classify supernovae (SNe) within a few days after explosion, if not earlier. Early-time spectra of young SNe directly probe the outermost layers of the ejecta, providing insights into the extent of stripping in the progenitor star and the explosion mechanism in the case of core-collapse supernovae. However, many SNe show overlapping observational characteristics at early time, complicating the early-time classification. In this paper, we focus on the study and classification of Type Ib supernovae (SNe Ib), which are a subclass of core-collapse supernovae that lack strong hydrogen lines but show helium lines in their spectra. Here we present a spectral dataset of 8 SNe Ib, chosen to have at least 3 pre-maximum spectra, which we call early spectra. Our dataset was obtained mainly by the the Las Cumbres Observatory (LCO) and consists of a total of 82 optical photospheric spectra, including 38 early spectra. This data set increases the number of published SNe Ib with at least three early spectra by ~60%. For our classification efforts, we use early spectra in addition to spectra taken around maximum light. We also convert our spectra into SN Identification (SNID) templates and make them available to the community for easier identification of young SNe Ib. Our data set increases the number of publicly available SNID templates of early spectra of SNe Ib by ~43%. Almost half of our sample has SN types that change over time or are different from what is listed on the Transient Name Server (TNS). We discuss the implications of our dataset and our findings for current and upcoming SN surveys and their classification efforts.

Context: HR2562B is a planetary-mass companion located 0.56arcsec (19au) from its host star. It is one of a few L/T transitional objects orbiting a young star. This companion provides insight into the evolution of young objects in the L/T transition. However, its key physical properties, such as Teff and mass, remain poorly constrained, with large uncertainties (34% for Teff, 22% for log(g)) based on near-infrared observations alone. Aims: We aim to refine these uncertainties, especially for Teff (1200-1700K) and log(g) (4-5), using new MIR data from the JWST/MIRI filters (10.65, 11.40, and 15.50 microns), and better understand the companion's chemical composition and its role in the L/T transition. Methods: MIRI data were processed using reference star differential imaging, revealing HR2562B at high S/N (16) in all 3 filters. We used 2 atmospheric models, ATMO and ExoREM, to fit the SED, combining MIR and NIR datasets. Additionally, we used CMD with brown dwarfs to explore the chemical composition of HR2562B's atmosphere and compare it to another L/T transition object, VHS1256b. Results: Our analysis improved the temperature precision (Teff=1255+-15K) by 6x compared to previous estimates. We also narrowed its luminosity to -4.69+-0.01 dex. Surface gravity remains uncertain (4.4-4.8), and its mass is estimated between 8 and 18.5Mj, depending on modeling and astrometry. Sensitivity analysis revealed the ability to detect objects between 2-5Mj at 100au. Conclusions: HR2562B likely has a near cloud-free atmosphere, with the ATMO model fitting better than ExoREM. Silicate absorption features are weak, requiring further spectroscopic observations. While HR2562B and VHS1256b share similarities, they are in different evolutionary stages, making HR2562B key to understanding young objects in the L/T transition. It is likely a planetary-mass companion, suggesting a reclassification as HR2562b.

Lacking terrestrial experimental data, our best constraints on the behavior of matter at high densities up to and above nuclear density arise from observations of neutron stars. Current constraints include those based on measurements of stellar masses, radii, and tidal deformabilities. Here we explore how orbits of primordial black holes - should they exist - inside neutron stars could provide complementary constraints on the nuclear equation of state (EOS). Specifically, we consider a sample of candidate EOSs, construct neutron star models for these EOSs, and compute orbits of primordial black holes inside these stars. We discuss how the pericenter advance of eccentric orbits, i.e. orbital precession, results in beat phenomena in the emitted gravitational wave signal. Observing this beat frequency could constrain the nuclear EOS and break possible degeneracies arising from other constraints, as well as provide information about the host star.

Young stellar objects are thought to commonly undergo sudden accretion events that result in a rise in bolometric luminosity. These outbursts likely coincide with the onset of planet formation, and could impact the formation of planets. The reason behind this dramatic enhancement of accretion is an active area of research, and the mass of the system is a critical parameter. Using Northern Extended Millimeter Array, we survey five outbursting sources (three FU Ori, one EX Or, one 'peculiar' source) with the primary goal of determining the system's mass using an optically thin line of CO. We estimate the mass of a central region for each object that using both continuum emission and C17O J=2-1. The C17O emission likely includes both disk and inner envelope material, thus acts as an upper limit on the disk mass, ranging from 0.33-3.4 Msun for our sources. These derived masses suggest that the inner approx. 1000 au contains enough mass along the line of sight for these sources to be gravitationally unstable.

Our understanding of the assembly timeline of the Milky Way has been transforming along with the dramatic increase in astrometric and spectroscopic data available over the past several years. Many substructures in chemo-dynamical space have been discovered and identified as the remnants of various galactic mergers. To investigate the timeline of these mergers we select main sequence turn off & subgiant stars (MSTOs) from the H3 survey, finding members in seven metal poor components of the halo: GSE, the Helmi Streams, Thamnos, Sequoia, Wukong/LMS-1, Arjuna, and I'itoi. We also select out the metal poor in situ disk to facilitate comparison to the evolution of the Milky Way itself at these early epochs. We fit individual isochrone ages to the MSTOs in each of these substructures and use the resulting age distributions to infer simple star formation histories. For GSE we resolve an extended star formation history that truncates $\approx10$ Gyr ago, as well as a clear age -- metallicity relation. From this age distribution and measured star formation history we infer that GSE merged with the Milky Way at a time $9.5-10.2$ Gyr ago, in agreement with previous estimates. We infer that the other mergers occurred at various times ranging from $9-13$ Gyr ago, and that the metal poor component of the disk built up within only a few billion years. These results reinforce the emerging picture that both the disk and halo of the Milky Way experienced a rapid assembly.

FU Ori and EX Lup type objects present natural experiments for understanding a critical stage in the star and planet formation process. These objects offer insight into the diversity of molecules available to forming planetary systems due to a sudden increase in accretion and central luminosity causes the disk and surrounding material to increase in temperature. This allows for volatiles to sublimate off of grains and exist in the gas-phase for tens to hundreds of years post initial outburst. While this dynamic stage may be common for solar-type protostars, observations of the chemical impact of these bursts are rare. In this article, we present observations from the NOrthern Extended Millimeter Array (NOEMA) of five Young Stellar Objects (YSOs) that have undergone outbursts within the past 100 years and catalog the volatile chemistry found within approx 1000 au of the YSO. Only one source clearly shows a line rich spectra with >11 molecules detected including complex organics and water, as is an expected spectra signature for a post-outburst source. This source is V1057 Cyg, and we present it as the northern analog to the well studied and molecule-rich FU Ori source, V883 Ori. Our conclusions on the chemical inventory of the other four sources in our sample are sensitivity limited, as V1057 Cyg contains the highest disk/envelope gas mass.

We present a study of the white dwarf (WD) cooling sequence (CS) in the globular cluster (GC) Omega Centauri, the primary goal of a dedicated Hubble Space Telescope (HST) programme. Our analysis has revealed that the peak at the termination of the WD CS is located at $m_{\rm F606W}$=30.1$\pm$0.2 (equivalent to $V$$\sim$31). The brighter part of Omega Centauri's WD CS is consistent with the presence of massive He-core WDs, in agreement with previous HST analyses with ultraviolet and blue filters. Comparative analyses of the WD luminosity function (LF) with theoretical counterparts have shown that a single-age population for the cluster is compatible with the data. However, an analysis of just the WD LF cannot entirely exclude the possibility of an age range, due to uncertainties in the present-day WD mass function, with a star formation history potentially spanning up to 5 billion years, predominantly comprising stars about 13 Gyr old, and with just a minority potentially as young as 8 Gyr. This underscores the need for global spectroscopic and photometric investigations that include simultaneously the WD populations together with the previous evolutionary phases to fully understand the cluster's diverse chemical compositions and ages.

While elementary particles are the favored candidate for the elusive dark matter, primordial black holes (PBHs) have also been considered to fill that role. Gravitational microlensing is a very well-suited tool to detect and measure the abundance of compact objects in galaxies. Previous studies based on quasar microlensing exclude a significant presence of substellar to intermediate-mass BHs ($\lesssim 100\,\mathrm{M}_\odot$). However, these studies were based on a spatially uniform distribution of BHs while, according to current theories of PBHs formation, they are expected to appear in clusters. We study the impact of clustering in microlensing flux magnification finding that at large scales clusters act like giant pseudo-particles, strongly affecting the emission coming from the Broad Line Region, which can no longer be used to define the zero microlensing baseline. As an alternative, we set this baseline from the intrinsic magnification ratios of quasar images predicted by macro lens models and compare them with the observed flux ratios in emission lines, infrared (IR), and radio. The (magnitude) differences are the flux-ratio anomalies attributable to microlensing, which we estimate for 35 image pairs corresponding to 12 lens systems. A Bayesian analysis indicates that the observed anomalies are incompatible with the existence of a significant population of clustered PBHs. Furthermore, we find that more compact clusters exhibit a stronger microlensing impact. Consequently, we conclude that clustering makes the existence of a significant population of BHs in the substellar to intermediate mass range even more unlikely.

The emission and escape of Lyman-$\alpha$ photons from star-forming galaxies is determined through complex interactions between the emitted photons and a galaxy's interstellar and circumgalactic gas, causing Lyman-$\alpha$ emitters (LAEs) to commonly appear not as point sources but in spatially extended halos with complex spectral profiles. We develop a 3D spatial-spectral model of Lyman-$\alpha$ halos (LAHs) to replicate LAH observations in integral field spectroscopic studies, such as those made with VLT/MUSE. The profile of this model is a function of 6 key halo properties: the halo- and compact-source exponential scale lengths ($r_{sH}$ and $r_{sC}$), the halo flux fraction ($f_H$), the compact component ellipticity ($q$), the spectral line width ($\sigma$), and the spectral line skewness parameter ($\gamma$). Placing a series of model LAHs into datacubes reflecting observing conditions in the MUSE UDF-Mosaic survey, we test their detection recoverability and determine that $\sigma$, $r_{sH}$, and $f_H$ are expected to have the most significant effect on the detectability of the overall LAH at a given central wavelength and intrinsic line luminosity. We develop a general selection function model spanning a grid of these halo parameters, and with a sample of 145 UDF-Mosaic LAHs with measured halo properties, we derive completeness-corrected, intrinsic distributions of the values of $\sigma$, $r_{sH}$, and $f_H$ for $3<z<5$ LAHs. We present best-fit functional forms of the distributions, and a $\sigma$ distribution corrected for instrumental line-spread function (LSF) broadening, and thereby show the physical line spread distribution of the intrinsic population. Finally, we discuss implications of these distributions for Ly$\alpha$ emission through the circumgalactic medium, finding that observations undercount LAHs with extended halo scale lengths compared to the intrinsic population.

We present cogsworth, an open-source Python tool for producing self-consistent population synthesis and galactic dynamics simulations. With cogsworth one can (1) sample a population of binaries and star formation history, (2) perform rapid (binary) stellar evolution, (3) integrate orbits through the galaxy and (4) inspect the full evolutionary history of each star or compact object, as well as their positions and kinematics. We include the functionality for post-processing hydrodynamical zoom-in simulations as a basis for galactic potentials and star formation histories to better account for initial spatial stellar clustering and more complex potentials. Alternatively, several analytic models are available for both the potential and star formation history. cogsworth can transform the intrinsic simulated population to an observed population through the joint application of dust maps, bolometric correction functions and survey selection functions. We provide a detailed explanation of the functionality of cogsworth and demonstrate its capabilities through a series of use cases: (1) We predict the spatial distribution of compact objects and runaways in both dwarf and Milky-Way-like galaxies, (2) using a star cluster from a hydrodynamical simulation, we show how supernovae can change the orbits of stars in several ways, and (3) we predict the separation of disrupted binary stellar companions on the sky and create a synthetic Gaia colour-magnitude diagram. We also discuss some current limitations and plans for future developments. We designed cogsworth and its online documentation to provide a powerful tool for constraining binary evolution, but also a flexible and accessible resource for the entire community.

The Zernike wavefront sensor (ZWFS) stands out as one of the most sensitive optical systems for measuring the phase of an incoming wavefront, reaching photon efficiencies close to the fundamental limit. This quality, combined with the fact that it can easily measure phase discontinuities, has led to its widespread adoption in various wavefront control applications, both on the ground but also for future space-based instruments. Despite its advantages, the ZWFS faces a significant challenge due to its extremely limited dynamic range, making it particularly challenging for ground-based operations. To address this limitation, one approach is to use the ZWFS after a general adaptive optics (AO) system; however, even in this scenario, the dynamic range remains a concern. This paper investigates two optical configurations of the ZWFS: the conventional setup and its phase-shifted counterpart, which generates two distinct images of the telescope pupil. We assess the performance of various reconstruction techniques for both configurations, spanning from traditional linear reconstructors to gradient-descent-based methods. The evaluation encompasses simulations and experimental tests conducted on the Santa cruz Extreme Adaptive optics Lab (SEAL) bench at UCSC. Our findings demonstrate that certain innovative reconstruction techniques introduced in this study significantly enhance the dynamic range of the ZWFS, particularly when utilizing the phase-shifted version.

Transit timing variation (TTV) provides rich information about the mass and orbital properties of exoplanets, which are often obtained by solving an inverse problem via Markov Chain Monte Carlo (MCMC). In this paper, we design a new data-driven approach, which potentially can be applied to problems that are hard to traditional MCMC methods, such as the case with only one planet transiting. Specifically, we use a deep learning approach to predict the parameters of non-transit companion for the single transit system with transit information (i.e., TTV, and Transit Duration Variation (TDV)) as input. Thanks to a newly constructed \textit{Transformer}-based architecture that can extract long-range interactions from TTV sequential data, this previously difficult task can now be accomplished with high accuracy, with an overall fractional error of $\sim$2\% on mass and eccentricity.

Context. On 2022 January 20, the Energetic Particle Detector (EPD) on board Solar Orbiter measured a solar energetic particle (SEP) event showing unusual first arriving particles from the anti-Sun direction. Near-Earth spacecraft separated 17° in longitude to the west from Solar Orbiter measured classic antisunward-directed fluxes. STEREO-A and MAVEN, separated 18° to the east and 143° to the west from Solar Orbiter respectively, also observed the event, suggesting that particles spread over at least 160° in the heliosphere. Results. Solar Orbiter was embedded in a MC erupting on 16 January from the same active region as the one related to the SEP event on 20 January. The SEP event is related to a M5.5 flare and a fast CME-driven shock of 1433 km/s, which injected particles within and outside the MC. The hard SEP spectra, the presence of a Type II radio burst, and the co-temporal Type III radio bursts being observed from 80 MHz that seems to emanate from the Type II, points to the shock as the relevant accelerator of the particles. Conclusions. The detailed analysis of the SEP event strongly suggest that the energetic particles are injected mainly by a CME-driven shock into and outside of a previous MC present in the heliosphere at the time of the particle onset. The sunward propagating SEPs measured by Solar Orbiter are produced by the injection of particles along the longer (western) leg of the MC still connected to the Sun at the time of the release of the particles. The determined electron propagation path length inside the MC is around 30% longer than the estimated length of the loop leg of the MC itself (based on the graduated cylindrical shell model) consistent with a low number of field line rotations.

Explosive energy release in the solar atmosphere is driven magnetically, but mechanisms triggering the onset of the eruption remain in debate. In the case of flares and CMEs, ideal or non-ideal instabilities usually occur in the corona, but direct observations and diagnostics there are difficult to obtain. To overcome this difficulty, we analyze observational signatures in the upper chromosphere or transition region, in particular, brightenings and dimmings at the feet of coronal magnetic structures. In this paper, we examine the time evolution of spatially resolved light curves in two eruptive flares, and identify a variety of tempo-spatial sequences of brightenings and dimmings, such as dimming followed by brightening, and dimming preceded by brightening. These brightening-dimming sequences are indicative of the configuration of energy release in the form of plasma heating or bulk motion. We demonstrate the potential of using these analyses to diagnose properties of magnetic reconnection and plasma expansion in the corona during the early stage of the eruption.

Abridged: Very Low Luminosity Objects (VeLLOs) are deeply embedded, and extremely faint objects and are thought to be in the quiescent phase of the episodic accretion process. They fill an important gap in our understanding of star formation. The VeLLO in the isolated DC3272+18 cloud has undergone an outburst, and is thus an ideal target for investigating the chemical inventory in the gas phase of an object of its type. Observations with the Atacama Pathfinder EXperiment (APEX) in four spectral windows in the frequency range of 213.6--272.4~GHz have been carried out to identify molecules that can be directly linked to the past outburst, utilize the line fluxes, column densities, and the abundance ratios of the detected species to characterize the different physical components of the VeLLO, and probe for the presence of complex organic molecules. Nitric oxide (NO) is detected for the first time in a source of this type, and its formation could be induced by the sublimation of grain-surface species during the outburst. A pathway to form NO directly in the gas phase is from the photodissociation products created after the sublimation of H$_2$O and NH$_3$ from the ices. While the present time water snowline has likely retreated to pre-outburst small radius, the volatile NO species is still extensively present in the gas phase, as evident by its high column density relative to methanol in the observations. This suggests that NO could be potentially used to trace the water snowline in outbursting sources. In order to rule out non-thermal desorption processes that could also have led to the formation of NO, this proposition has to be verified with future observations at higher spatial resolution, and by searching for NO in additional targets.

We present a complete analysis of Fermi Large Area Telescope (LAT) data of GRB 221009A, the brightest Gamma-Ray Burst (GRB) ever detected. The burst emission above 30 MeV detected by the LAT preceded by 1 s the low-energy (< 10 MeV) pulse that triggered the Fermi Gamma-Ray Burst Monitor (GBM), as has been observed in other GRBs. The prompt phase of GRB 221009A lasted a few hundred seconds. It was so bright that we identify a Bad Time Interval (BTI) of 64 seconds caused by the extremely high flux of hard X-rays and soft gamma rays, during which the event reconstruction efficiency was poor and the dead time fraction quite high. The late-time emission decayed as a power law, but the extrapolation of the late-time emission during the first 450 seconds suggests that the afterglow started during the prompt emission. We also found that high-energy events observed by the LAT are incompatible with synchrotron origin, and, during the prompt emission, are more likely related to an extra component identified as synchrotron self-Compton (SSC). A remarkable 400 GeV photon, detected by the LAT 33 ks after the GBM trigger and directionally consistent with the location of GRB 221009A, is hard to explain as a product of SSC or TeV electromagnetic cascades, and the process responsible for its origin is uncertain. Because of its proximity and energetic nature, GRB 221009A is an extremely rare event.

Stellar age estimates are often calculated by interpolating a star's properties in a grid of models. However, different model grids will give different ages for the same star. We used the open cluster M67 to compare four different model grids: DSEP, GARSTEC, MIST, and YREC. Across all model grids, age estimates for main sequence stars were consistently higher than the accepted age of M67, while age estimates for red giant stars were lower. We compared model-generated age and mass values to external constraints as an additional test of the reliability of each model grid. For stars near solar age and metallicity, we recommend using the DSEP model grid to estimate the ages of main sequence stars and the GARSTEC model grid for red giant stars.

Sub-parsec supermassive black hole (SMBH) binaries are expected to be common in active galactic nuclei (AGN), as a result of the hierarchical build-up of galaxies via mergers. While direct evidence for these compact binaries is lacking, a few hundred candidates have been identified, most based on the apparent periodicities of their optical light-curves. Since these signatures can be mimicked by AGN red-noise, additional evidence is needed to confirm their binary nature. Recurring self-lensing flares (SLF), occurring whenever the two BHs are aligned with the line of sight within their Einstein radii, have been suggested as additional binary signatures. Furthermore, in many cases, lensing flares are also predicted to contain a "dip", whenever the lensed SMBH's shadow is comparable in angular size to the binary's Einstein radius. This feature would unambiguously confirm binaries and additionally identify SMBH shadows that are spatially unresolvable by high-resolution VLBI. Here we estimate the number of quasars for which these dips may be detectable by LSST, by extrapolating the quasar luminosity function to faint magnitudes, and assuming that SMBH binaries are randomly oriented and have mass-ratios following those in the Illustris simulations. Under plausible assumptions about quasar lifetimes, binary fractions, and Eddington ratios, we expect tens of thousands of detectable flares, of which several dozen contain measurable dips.

Parametric strong lensing studies of galaxy clusters often display "misleading features". This is the case in the galaxy cluster Abell 370. Using strong lensing techniques, it has been described parametrically by a four dark matter clumps model and galaxy scale perturbers, as well as a significant external shear component, which physical origin remains a challenge. The dark matter distribution features a mass clump with no stellar counterpart and a significant offset between one of the dark matter clumps and its associated stellar counterpart. In this paper, based on BUFFALO data, we begin by revisiting this mass model. We find a four dark matter clumps solution which does not require any external shear and provides a slightly better RMS compared to previous models. Investigating further this new solution, we present a class of models which can accurately reproduce the strong lensing data, but whose parameters for the dark matter component are poorly constrained. We then develop a model where each large scale dark matter component must be associated with a stellar counterpart. This model is unable to reproduce the observational constraints with an RMS smaller than 2.3", and the parameters describing this dark matter component remain poorly constrained. Examining the total projected mass maps, we find a good agreement between the total mass and the stellar distribution, both being bimodal. We interpret the "misleading features" of the four dark matter clumps mass model and the failure of the three dark matter clumps mass model as being symptomatic of the lack of realism of a parametric description of the dark matter distribution, and encourage caution and criticism on the outputs of parametric strong lensing modelling. We briefly discuss the implications of our results for using Abell 370 as a gravitational telescope.

A successful theory of star formation should predict the number of objects as a function of their mass produced through star-forming events. Previous studies in star-forming regions and the solar neighborhood identify a mass function increasing from the hydrogen-burning limit down to about 10 M$_{J}$. Theory predicts a limit to the fragmentation process, providing a natural turnover in the mass function down to the opacity limit of turbulent fragmentation thought to be 2-10 M$_{J}$. Programs to date have not been sensitive enough to probe the hypothesized opacity limit of fragmentation. Here we present the first identification of a turnover in the initial mass function below 12 M$_{J}$ within NGC 2024, a young star-forming region. With JWST/NIRCam deep exposures across 0.7-5 {\mu}m, we identified several free floating objects down to ~ 3 M$_{J}$ with sensitivity to 0.5 M$_{J}$. We present evidence for a double power law model increasing from about 60 M$_{J}$ to roughly 12 M$_{J}$, consistent with previous studies, followed by a decrease down to 0.5 M$_{J}$. Our results support the predictions of star and brown dwarf formation theory, identifying the theoretical turnover in the mass function and suggest the fundamental limit of turbulent fragmentation near 3 M$_{J}$.

Investigating the metal enrichment in the early universe helps us constrain theories about the first stars and study the ages of galaxies. The lensed galaxy MACS0647$-$JD at $z=10.17$ is the brightest galaxy known at $z > 10$. Previous work analyzing JWST NIRSpec and MIRI data yielded a direct metallicity $\rm{12+log(O/H)}=7.79\pm0.09$ ($\sim$ 0.13 $Z_\odot$) and electron density $\rm{log}(n_e / \rm{cm^{-3}}) = 2.9 \pm 0.5$, the most distant such measurements to date. Here we estimate the direct C/O abundance for the first time at $z > 10$, finding a sub-solar ${\rm log(C/O)}=-0.44^{+0.06}_{-0.07}$. This is higher than other $z>6$ galaxies with direct C/O measurements, likely due to higher metallicity. It is also slightly higher than galaxies in the local universe with similar metallicity. This may suggest a very efficient and rapid burst of star formation, a low effective oxygen abundance yield, or the presence of unusual stellar populations including supermassive stars. Alternatively, the strong CIII]${\rm {\lambda}{\lambda}}$1907,1909 emission ($14\pm 3\,{Å}$ rest-frame EW) may originate from just one of the two component star clusters JDB ($r \sim 20$ pc). Future NIRSpec IFU spectroscopic observations of MACS0647$-$JD will be promising for disentangling C/O in the two components to constrain the chemistry of individual star clusters just 460 Myr after the Big Bang.

We present a catalog of 71,364 point-like UV sources with SDSS photometry and GALEX FUV-NUV less or equal 0.1mag. The limit corresponds to stellar Teff greater than 15,000 to 20,000K, slightly depending on gravity but nearly reddening-independent for Milky-Way-type dust. Most sources are hot white-dwarfs (WDs) and sub-dwarfs (SDs). Comparing the SED (GALEX FUV, NUV, SDSS u,g,r,i,z) of 35,294 sources having good photometry with colors of stellar models and known objects, we identify 12,404+1871-1267 binary hot-compact stars with a cooler, less-evolved companion (with a possible 8% to 15% contamination by low-redshift QSOs), and 22,848+1267-3853 single-star candidates. Single-star counts are an upper limit because pairs of similar stars have single-star-like SED, and hot-WDs with main-sequence companions of certain types (depending on WD's radius) are missed or counted as single in the available wavelength range and selection. The catalog offers unique leverage for identifying hot WDs, elusive at longer wavelengths when a cooler, larger companion dominates optical-IR fluxes: 51% of the binary- and 20% of the single-star candidates are previously un-known objects. Gaia DR3 provides a parallax with error less than 20% for 34% of the binaries- and 45% of single-star candidates, allowing Teff, E(B-V), radius and Lbol to be derived from SED analysis. The binary-candidate sample usefully expands the overall current binary-WD census to subpopulations elusive to Gaia and to other searches. The binary fraction among this specific sample of hot-compact objects, albeit with the mentioned biases, B_f>46%, compared with that of their progenitors (>80%-50% for mass range 8-1Msun, Moe (2019)), implies a lower merging rate than found for massive stars by Sana et al (2017).

When a telescope doesn't reach a reasonable point spread function on the detector or detectable wavefront quality after initial assembly, a coarse phase alignment on-sky is crucial. Before utilizing a closed loop adaptive optics system, the observatory needs a strategy to actively align the telescope sufficiently for fine wavefront sensing. This paper presents a method of early-stage alignment using a stochastic parallel-gradient-descent (SPGD) algorithm which performs random perturbations to the optics of a three mirror anastigmat telescope design. The SPGD algorithm will drive the telescope until the wavefront error is below the acceptable range of the fine adaptive optics system to hand the telescope over. The focused spot size over the field of view is adopted as a feed parameter to the SPGD algorithm and wavefront peak-to-valley error values are monitored to directly compare our mechanical capabilities to our alignment goal of diffraction limited imaging and fine wavefront sensing.

Type Ia supernovae (SNe Ia) constitute an historical probe to derive cosmological parameters through the fit of the Hubble-Lemaître diagram, i.e. SN Ia distance modulus versus their redshift. In the era of precision cosmology, realistic simulation of SNe Ia for any survey entering in an Hubble-Lemaître diagram is a key tool to address observational systematics, like Malmquist bias. As the distance modulus of SNe Ia is derived from the fit of their light-curves, a robust simulation framework is required. In this paper, we present the performances of the simulation framework skysurvey to reproduce the the Zwicky Transient Facility (ZTF) SN Ia DR2 covering the first phase of ZTF running from April 2018 up to December 2020. The ZTF SN Ia DR2 sample correspond to almost 3000 classified SNe Ia of cosmological quality. First, a targeted simulation of the ZTF SN Ia DR2 was carried on to check the validity of the framework after some fine tuning of the observing conditions and instrument performance. Then, a realistic simulation has been run using observing ZTF logs and ZTF SN Ia DR2 selection criteria on simulated light-curves to demonstrate the ability of the simulation framework to match the ZTF SN Ia DR2 sample. Furthermore a redshift dependency of SALT2 light-curve parameters (stretch and colour) was conducted to deduce a volume limited sample, i.e. an unbiased SNe Ia sample, characterized with $z_{lim} \leq 0.06$. This volume limited sample of about 1000 SNe Ia is unique to carry on new analysis on standardization procedure with a precision never reached (those analysis are presented in companion papers).

The dispersion measures of fast radio bursts have been identified as a powerful tool for testing the zero-mass hypothesis of the photon. The classical approach treats the massive photon-induced and plasma-induced time delays as two separate phenomena. Recently, Wang et al. (2024) suggested that the joint influence of the nonzero photon mass and plasma effects should be considered, and proposed a revised time delay for massive photons propagating in a plasma medium, denoted as $\Delta t'_{m_{\gamma}} \propto \nu^{-4}$, which departures from the classical dispersion relation ($\propto \nu^{-2}$). Here we discuss the derivation presented by Wang et al. (2024) and show that the classical dispersion relation remains valid based on Proca equations.

We present the first joint spectral and imaging analysis of hard X-ray (HXR) emission from 3 microflares observed by the Nuclear Spectroscopic Telescope ARray (NuSTAR) and Solar Orbiter/Spectrometer/Telescope for Imaging X-rays (STIX). We studied 5 joint spectra from GOES A7, B1 and B6 class microflares from active region AR12765 on 2020 June 6 and 7. As these events are very bright for NuSTAR, resulting in extremely low (<1%) livetime, we introduce a pile-up correction method. All five joint spectra were fitted with an isothermal model finding temperatures in the 9-11 MK range. Furthermore, three joint spectra required an additional non-thermal thick-target model finding non-thermal powers of $10^{25}$-$10^{26}$ erg s$^{-1}$. All the fit parameters were within the ranges expected for HXR microflares. The fit results give a relative scaling of STIX and NuSTAR mostly between 6-28% (one outlier at 52%) suggesting each instrument are well calibrated. In addition to spectral analysis, we performed joint HXR imaging of the June 6 and one of the June 7 microflares. In NuSTAR's field of view (FOV), we observed two separate non-thermal sources connected by an elongated thermal source during the June 6 microflares. In STIX's FOV (44 degrees W with respect to NuSTAR), we imaged thermal emission from the hot flare loops which when reprojected to an Earth viewpoint matches the thermal sources seen with NuSTAR and in the hotter EUV channels with the Solar Dynamic Observatory's Atmospheric Imaging Assembly.

Millisecond pulsars (MSPs) are laboratories for stellar evolution, strong gravity, and ultra-dense matter. Although MSPs are thought to originate in low-mass X-ray binaries (LMXBs), approximately 27% lack a binary companion, while others are found in systems with large orbital eccentricities. Understanding how these systems form may provide insight into the internal properties of neutron stars (NSs). We study the formation of a twin compact star through rapid first-order phase transitions in NS cores due to mass accretion in LMXBs. We investigate whether this mechanism, possibly coupled with secondary kick effects such as neutrino or electromagnetic rocket effects, may leave an observable long-lasting imprint on the orbit. We simulate mass accretion in LMXBs consisting of a NS and a low-mass main sequence companion, following the evolution of the NS mass, radius, and spin until a strong phase transition is triggered. For the NS structure, we assume a multipolytrope equation-of-state that allows for a sharp phase transition from hadronic to quark matter and satisfies observational constraints. We find that in compact binaries with relatively short pre-Roche lobe overflow orbital periods, an accretion-induced phase transition may occur during the LMXB phase. In contrast, in systems with wider orbits, this transition may take place during the spin-down phase, forming an eccentric binary MSP. If the transition is accompanied by a secondary kick (w > 20 km/s), the binary is likely to be disrupted, forming an isolated MSP or reconfigured to an ultra-wide orbit. Our results suggest that accretion in LMXBs provides a viable path for forming twin compact stars, potentially leaving an observable imprint on the orbit. The eccentricity distribution of binary MSPs with long (> 50 d) orbital periods could provide constraints on first-order phase transitions in dense nuclear matter.

We present a multi-wavelength spectral study of NGC 4151 based on five epochs of simultaneous AstroSat observations in the near ultra-violet (NUV) to hard X-ray band ($\sim 0.005-80$ keV) during $2017 - 2018$. We derived the intrinsic accretion disk continuum after correcting for internal and Galactic extinction, contributions from broad and narrow line regions, and emission from the host galaxy. We found a bluer continuum at brighter UV flux possibly due to variations in the accretion disk continuum or the UV reddening. We estimated the intrinsic reddening, $E(B-V) \sim 0.4$, using high-resolution HST/STIS spectrum acquired in March 2000. We used thermal Comptonization, neutral and ionized absorption, and X-ray reflection to model the X-ray spectra. We obtained the X-ray absorbing neutral column varying between $N_H \sim 1.2-3.4 \times 10^{23} cm^{-2}$, which are $\sim 100$ times larger than that estimated from UV extinction, assuming the Galactic dust-to-gas ratio. To reconcile this discrepancy, we propose two plausible configurations of the obscurer: (a) a two-zone obscurer consisting of dust-free and dusty regions, divided by the sublimation radius, or (b) a two-phase obscurer consisting of clumpy, dense clouds embedded in a low-density medium, resulting in a scenario where a few dense clouds obscure the compact X-ray source substantially, while the bulk of UV emission arising from the extended accretion disk passes through the low-density medium. Furthermore, we find a positive correlation between X-ray absorption column and $NUV-FUV$ color and UV flux, indicative of enhanced winds possibly driven by the 'bluer-when-brighter' UV continuum.

We re-consider the already published phenomenon: the blue shift of diffuse interstellar bands, observed in spectra of HD34078 (AE Aur) and members of the Sco OB1 association, in particular HD152233. We have analyzed 29 diffuse bands. Part of them, already proven as blue-shifted in our earlier study, are now confirmed using another instrument: this time the 6.5m Clay telescope equipped with the MIKE spectrograph. The high signal-to-noise ratio (over 600) of our spectra allowed us to reveal even small small scale {\bf displacements} of positions (both: blue and red-shifts) of diffuse bands along the considered lines of sight. In some cases the magnitude of deviation exceeds 10 km\,s$^{-1}$. Also, we prove that profiles of many diffuse bands in spectra of HD34078 suffer significant broadening. The origin of the observed phenomena is discussed.

Gamma-ray binary LS I+61$^{\circ}$ 303 consists of a neutron star orbiting around a Be star with a period of $P_{\rm orb}\simeq26.5\ {\rm d}$. Apart from orbital modulations, the binary shows long-term flux variations with a superorbital period of $P_{\rm sup}\simeq4.6\ {\rm yrs}$ as seen in nearly all wavelengths. The origin of this superorbital modulation is still not well understood. Under the pulsar wind-stellar outflow interaction scenario, we propose that the superorbital modulations of LS I+61$^{\circ}$ 303 could be caused by the precession of the Be disk. Assuming X-rays arise from synchrotron radiation of the intrabinary shock, we develop an analytical model to calculate expected flux modulations over the orbital and superorbital phases. The asymmetric two-peak profiles in orbital light curves and sinusoidal-like long-term modulations are reproduced under the precessing disk scenario. The observed orbital phase drifting of the X-ray peak and our fitting of long-term X-ray data indicate that the neutron star is likely orbiting around the star with a small eccentricity and periastron phase around $\Phi_{\rm p}\sim0.6$. We compare the Corbet diagrams of LS I+61$^{\circ}$ 303 with other Be/X-ray binaries and the linear correlation in the $P_{\rm sup}-P_{\rm orb}$ diagram suggests that the precession of the Be disk in LS I+61$^{\circ}$ 303 is induced by the tidal torque of its neutron star companion.

We have used the Dark Energy Camera (DECam) on the CTIO Blanco 4-m telescope to perform a new emission-line survey of the Large Magellanic Cloud (LMC) using narrow-band H-alpha and [SII] filters in addition to a continuum band for use in creating pure emission-line images. We refer to this new survey as DeMCELS, to distinguish it from the earlier Magellanic Cloud Emission Line Survey (MCELS) that has been in service for nearly 25 years. DeMCELS covers $\sim 54$ degrees$^{2}$, encompassing most of the bright optical disk of the LMC. With DECam's pixel size of only 0.27", our DeMCELS survey provides a seeing-limited improvement of 3-5 times over MCELS and is comparable in depth, with surface brightness limits of 3.3E-17 erg cm$^{-2}$ s$^{-1}$ arcsec$^{-2}$ in H-alpha and 2.9E-17 erg cm$^{-2}$ s$^{-1}$ arcsec$^{-2}$ in H-alpha and [SII], respectively. DeMCELS provides detailed morphological information on nebulae of all scales, from the largest supershells to individual [HII] regions and supernova remnants, to bubbles of emission surrounding individual stars, and even to faint structures in the diffuse ionized gas of the LMC. Many complex regions of emission show significant variations in the ratio of [SII] to H-alpha, a sign of the mixture of shocks from stellar winds and/or supernovae with photoionization by embedded hot, young stars. We present the details of the observing strategy and data processing for this survey, and show selected results in comparison with previous data. A companion project for the Small Magellanic Cloud is in progress and will be reported separately. We are making these new data available to the community at large via the NOIRLab's Data Lab site.

We have analyzed the prompt and afterglow characteristics of the intermediate luminosity burst ``GRB 210210A". Our prompt emission analysis indicates that GRB 210210A is among the softest long GRBs detected by the Swift-BAT. The time-integrated prompt emission spectrum of GRB 210210A is aptly described by a power law function with an exponential cutoff. The spectral peak energy (E$_{p,z}$) in rest-frame and the E$_{\rm \gamma, iso}$ for this GRB marginally satisfy the 2$\sigma$ Amati correlation, a common feature observed in low/intermediate luminosity GRBs. Notably, an early bump is observed in the Swift-XRT light curve (a rare feature); the optical afterglow light curve, on the other hand, appears to follow a power law decay. However, due to the lack of sufficient early optical observations, we cannot completely rule out the possibility of an early bump in the optical light curve. For the bump observed in the early X-ray light curve, we calculated parameters such as peak time, rise time, decay time, and bulk Lorentz factor ($\Gamma_{0}$ $\sim$ 156), which perfectly satisfy the correlation between the parameters of the onset of the afterglow in GRBs. Both the optical and X-ray (including our observations) light curves exhibit a chromatic break in the late afterglow. Based on the prompt and afterglow parameters, we confirm that the intermediate luminosity GRB 210210A favors a collapsar scenario and is possibly powered by a magnetar.

This paper presents a catalog of the energy distribution in the spectra of 263 stars in the wavelength range from 0.4 to 100 {\mu}m, which are at late stages of evolution and have been observed by the ISO space observatory. For each object in the catalog, estimates of the observed bolometric fluxes were derived from smoothed energy distribution curves. The catalog is available at this https URL\_lss\_sed both as a table and in machine-readable format. It is shown that for the specified sample of objects their ISO SWS spectra in the range 2.4-45 {\mu}m only 60% of cases correspond to the general shape of the continuum, and can be used without recalibration. A selection of carbon stars, accessible for the infrared observations from the MSU observatories has been made. For some of them the first brightness estimates in the K, L, and M bands were obtained with the new IR camera of the 2.5-m telescope of CMO.

The popularity of pulsating stars resides in their capacity of determining several crucial and relevant parameters such as heliocentric distances, ages, metallicity gradients and reddening. RR Lyrae stars are old stellar tracers and have been detected in nearly all nearby galaxies that have been searched for these stars, with just a few exceptions of very low mass dwarfs. Less common but also of great importance are Anomalous Cepheids, indicators of either old or intermediate-age population, depending on their stellar origin. Classical Cepheids are only found within young stellar populations, and because of their brighter absolute magnitudes, they can be detected in galaxies farther than the Local Group. This paper presents a concise review built upon the aforementioned pulsating stars in Local Group dwarf galaxies and some of their applications to infer important properties of their host galaxies.

This paper examines nucleosynthesis in a low-mass neutron star crust that loses mass due to accretion in a close binary system and, reaching a hydrodynamically unstable configuration explodes. The r-process proceeds mainly in the inner crust. Nucleosynthesis in the outer crust is an explosive process with a sharp increase in temperature caused by an outward-propagating shockwave (shock-wave nucleosynthesis). The number of heavy elements produced in a low-mass neutron star crust during the explosion is approximately equals 0.041 solar masses, which exceeds the number of heavy elements ejected as jets in the neutron star merger scenario.

The origins of lenticular galaxies (S0s) can be classified into two main categories: ``minor mergers" in low-density environments (LDEs) and ``faded spirals" in high-density environments (HDEs). The transitional phase in the evolution of S0s, namely, star-forming lenticular galaxies (SFS0s), can serve as an important probe for analyzing the complex processes involved in the transformation between different galaxy types and the quenching of star formation (SF). We attempt to find the impact of different environments on the global properties and spatially resolved quantities of SFS0s. We selected 71 SFS0s from the SDSS-IV MaNGA Survey, comprising 23 SFS0s in HDEs (SFS0s$\_$HE) and 48 SFS0s in LDEs (SFS0s$\_$LE). We examined the effects of the environment, by studying the global properties, concentration index, and radial profiles of the derived quantities. The varied environments of SFS0s do not lead to any significant difference in global properties (e.g., S$\acute{\rm e}$rsic index). By calculating $CI_{\rm H_{\alpha}/cont}$, we observe that different environments may cause varying concentrations of SF. Specifically, SFS0s$\_$LE, affected by external gas mergers or inflow, exhibit a more centrally concentrated SF (i.e., larger $CI_{\rm H_{\alpha}/cont}$). This trend is further supported by $CI_{\rm SFR, H_{\alpha}}$, which only considers the gas disk of the galaxy. This observation is aligned with the observed shrinking of gas disks in galaxies affected by ram-pressure stripping in HDEs. Furthermore, their $\Sigma_{\rm SFR}$ or resolved sSFR are comparable. On average, SFS0s$\_$LE display significantly higher values for both quantities. Finally, the observed D$_{\rm n}4000$ and gas-phase metallicity gradient correspond well to their assumed origins. However, we did not find a significantly lower gas-phase metallicity in SFS0s$\_$LE. Abridged

Swift J1858.6$-$0814 (hereafter J1858) is a transient neutron star low-mass X-ray binary (NS LMXB). There is controversy regarding its donor mass derived from observations and theoretical calculations. In this paper, we adopt seven magnetic braking (MB) prescriptions suggested in the literature and different metallicity $Z$ to simulate the evolution of the LMXB. Our results show that, employing the MB model proposed by \citet{2012ApJ...746...43R} ("rm12"), the Convection And Rotation Boosted ("carb") model \citep{2019ApJ...886L..31V}, as well as the Intermediate ("inter") and Convection-boosted ("cboost") models in \citet{2019MNRAS.483.5595V} can match (part of) the observational parameters of J1858 well. We then apply our method to other observed LMXBs and find that the "rm12" and "inter" MB laws are most promising in explaining transient LMXBs. In comparison, the simulations with the "cboost" and "carb" MB laws are more inclined to reproduce persistent LMXBs and ultra-compact X-ray binaries (UCXBs), respectively. Our results, though subject to computational and/or observational bias, show that it is challenging to find a unified MB law that applies to the NS LMXB sub-populations simultaneously, indicating our lack of understanding of the true MB law. In addition, we explore the influence of various MB laws on the magnitude of the bifurcation periods in LMXBs.

The detection of gravitational waves (GWs) from binary neutron stars (BNSs) with possible telescope follow-ups opens a window to ground-breaking discoveries in the field of multi-messenger astronomy. With the improved sensitivity of current and future GW detectors, more BNS detections are expected in the future. Therefore, enhancing low-latency GW search algorithms to achieve rapid speed, high accuracy, and low computational cost is essential. One innovative solution to reduce latency is the use of machine learning (ML) methods embedded in field-programmable gate arrays (FPGAs). In this work, we present a novel \texttt{WaveNet}-based method, leveraging the state-of-the-art ML model, to produce early-warning alerts for BNS systems. Using simulated GW signals embedded in Gaussian noise from the Advanced LIGO and Advanced Virgo detectors' third observing run (O3) as a proof-of-concept dataset, we demonstrate significant performance improvements. Compared to the current leading ML-based early-warning system, our approach enhances detection accuracy from 66.81\% to 76.22\% at a 1\% false alarm probability. Furthermore, we evaluate the time, energy, and economical cost of our model across CPU, GPU, and FPGA platforms, showcasing its potential for deployment in real-time gravitational wave detection pipelines.

The circular velocity function (CVF) of galaxies is a fundamental test of the $\Lambda$ Cold Dark Matter (CDM) paradigm as it traces the variation of galaxy number densities with circular velocity ($v_{\rm{circ}}$), a proxy for dynamical mass. Previous observational studies of the CVF have either been based on \ion{H}{I}-rich galaxies, or encompassed low-number statistics and probed narrow ranges in $v_{\rm{circ}}$. We present a benchmark computation of the CVF between $100-350\ \rm{km\ s^{-1}}$ using a sample of 3527 nearby-Universe galaxies, representative for stellar masses between $10^{9.2}-10^{11.9} \rm{M_{\odot}}$. We find significantly larger number densities above 150 $\rm{km\ s^{-1}}$ compared to results from \ion{H}{I} surveys, pertaining to the morphological diversity of our sample. Leveraging the fact that circular velocities are tracing the gravitational potential of halos, we compute the halo mass function (HMF), covering $\sim$1 dex of previously unprobed halo masses ($10^{11.7}-10^{12.7} \rm{M_{\odot}}$). The HMF for our sample, representative of the galaxy population with $M_{200}\geqslant10^{11.35} \rm{M_{\odot}}$, shows that spiral morphologies contribute 67 per cent of the matter density in the nearby Universe, while early types account for the rest. We combine our HMF data with literature measurements based on \ion{H}{I} kinematics and group/cluster velocity dispersions. We constrain the functional form of the HMF between $10^{10.5}-10^{15.5} \rm{M_{\odot}}$, finding a good agreement with $\Lambda$CDM predictions. The halo mass range probed encompasses 72$\substack{+5 \\ -6}$ per cent ($\Omega_{\rm{M,10.5-15.5}} = 0.227 \pm 0.018$) of the matter density in the nearby Universe; 31$\substack{+5 \\ -6}$ per cent is accounted for by halos below $10^{12.7}\rm{M_{\odot}}$ occupied by a single galaxy.

Stability of Hilda Asteroids in the solar system around the 3:2 resonance point is analyzed in terms of the Sun-Jupiter-asteroid elliptic restricted three-body problem. We show that the Hamiltonian of the system is well-approximated by a single-resonance Hamiltonian around the 3:2 resonance. This implies that orbits of the Hilda asteroids are approximately integrable, thus their motion is stable. This is in contrast to other resonances such as the 3:1 and the 2:1 resonances at which Kirkwood gaps occur. Indeed, around the 3:1 and the 2:1 resonances, the Hamiltonians are approximated by double-resonance Hamiltonians that are nonintegrable and thus indicate chaotic motions. By a suitable canonical transformation, we reduce the number of degrees of freedom for the system and derive a Hamiltonian which has two degrees of freedom. As a result, we can analyze the stability of the motion by constructing Poincare surface of section.

In this contribution we report the on-going progresses of the project FATE, an operational automatic forecast system conceived to deliver forecasts of a set of astroclimatic and atmospheric parameters having the aim to support the science operations (i.e. the Service Mode) at the Very Large Telescope. The project has been selected at conclusion of an international open call for tender opened by ESO and it fits with precise technical specifications. In this contribution we will present the ultimate goals of this service once it will be integrated in the VLT operations, the forecasts performances at present time and the state of the art of the project. FATE is supposed to draw the roadmap towards the optical turbulence forecast for the ELT.

In this contribution we present preliminary results of a study applied to the Observatories of Roque de Los Muchachos (La Palma) and Teide (Tenerife) in Canary Islands aiming to investigate the possibility to implement an automatic system for the optical turbulence forecasting for the European Solar Telescope (EST) telescope. The study has been carried out in the context of the SOLARNET project and the two mentioned sites were the pre-selected sites for EST. This analysis aimed to investigate the possibility to extend the methodology of the forecast of the optical turbulence developed by our team and performed on top-class ground-based telescopes dedicated to night time observations such as ALTA (@ LBT) and FATE (@ VLT) to the day-time regime. As an ancillary output our very preliminary analysis concludes, that the two sites of Roque de Los Muchachos Observatory (ORM) and Teide Observatory (TO) show comparable characteristics during the day time. Considering that the site of EST has been already identified to be at ORM this can be considered a very useful information from a scientific point of view.

Full-Stokes polarimetric datasets, originating from slit-spectrograph or narrow-band filtergrams, are routinely acquired nowadays. The data rate is increasing with the advent of bi-dimensional spectropolarimeters and observing techniques that allow long-time sequences of high-quality observations. There is a clear need to go beyond the traditional pixel-by-pixel strategy in spectropolarimetric inversions by exploiting the spatiotemporal coherence of the inferred physical quantities. We explore the potential of neural networks as a continuous representation of the physical quantities over time and space (also known as neural fields), for spectropolarimetric inversions. We have implemented and tested a neural field to perform the inference of the magnetic field vector (approach also known as physics-informed neural networks) under the weak-field approximation (WFA). By using a neural field to describe the magnetic field vector, we can regularize the solution in the spatial and temporal domain by assuming that the physical quantities are continuous functions of the coordinates. We investigated the results in synthetic and real observations of the Ca II 8542 A line. We also explored the impact of other explicit regularizations, such as using the information of an extrapolated magnetic field, or the orientation of the chromospheric fibrils. Compared to the traditional pixel-by-pixel inversion, the neural field approach improves the fidelity of the reconstruction of the magnetic field vector, especially the transverse component. This implicit regularization is a way of increasing the effective signal-to-noise of the observations. Although it is slower than the pixel-wise WFA estimation, this approach shows a promising potential for depth-stratified inversions, by reducing the number of free parameters and inducing spatio-temporal constraints in the solution.

Lenticular galaxies (S0s) in the local universe are generally absent of recent star formation and lack molecular gas. In this paper, we investigate one massive ($M_*$$\sim$5$\times10^{10}$ M$_\odot$) star-forming S0, PGC 39535, with the Northern Extended Millimeter Array (NOEMA). Using optical data from SDSS-IV MaNGA survey, we find star formation mainly concentrates in the central region of PGC 39535. The total star formation rate estimated using extinction-corrected H$\alpha$ flux is 1.57 M$_\odot$ yr$^{-1}$. Results of NOEMA observation suggest that the molecular gas mainly concentrates in the central regions as a gaseous bar and a ring-like structure, and shows similar kinematics as the stellar and ionized gas components. The total molecular gas mass estimated from CO(1-0) is (5.42$\pm$1.52)$\times$10$^{9}$ M$_{\odot}$. We find PGC 39535 lies on the star-forming main sequence, but falls below Kennicutt-Schmidt relation of spiral galaxies, suggesting that the star formation efficiency may be suppressed by the massive bulge. The existence of a second Gaussian component in the CO spectrum of the central region indicates possible gas flows. Furthermore, our analyses suggest that PGC 39535 resides in the center of a massive group and the derived star formation history indicates it may experience a series of gas-rich mergers over the past 2$\sim$7 Gyr.

Feedback is widely recognized as an essential condition for Lyman continuum (LyC) escape in star-forming galaxies. However, the mechanisms by which galactic outflows clear neutral gas and dust remain unclear. In this paper, we model the Mg II 2796Å, 2804Å absorption + emission lines in 29 galaxies taken from the Low-z LyC Survey (LzLCS) to investigate the impact of (radiation + mechanical) feedback on LyC escape. Using constraints on Mg$^+$ and photoionization models, we map the outflows' neutral hydrogen content and predict $f_{esc}^{LyC}$ with a multiphase wind model. We measure mass, momentum, and energy loading factors for the neutral winds, which carry up to 10% of the momentum and 1% of the energy in SFR-based deposition rates. We use SED template fitting to determine the relative ages of stellar populations, allowing us to identify radiation feedback dominant systems. We then examine feedback related properties (stellar age, loading factors, etc.) under conditions that optimize feedback efficiency, specifically high star formation rate surface density and compact UV half-light radii. Our findings indicate that the strongest leakers are radiation feedback dominant, lack Mg II outflows, but have extended broad components in higher ionization lines like [O III] 5007Å, as observed by Amorín et al. (2024). In contrast, galaxies experiencing supernovae feedback typically exhibit weaker $f_{esc}^{LyC}$ and show evidence of outflows in both Mg II and higher ionization lines. We attribute these findings to rapid or "catastrophic" cooling in the radiation-dominant systems, which, given the low metallicities in our sample, are likely experiencing delayed supernovae.

The neutral hydrogen (HI) power spectrum, measured from intensity fluctuations in the 21-cm background, offers insights into the large-scale structures (LSS) of our Universe in the post-reionization era (redshift $z<6$). A significant amount of HI is expected to reside in low- and intermediate-density environments, but the power spectrum mainly captures information from high-density regions. To more fully extract the information contained in the HI field, we investigate the use of a marked power spectrum statistic. Here, the power spectrum is effectively re-weighted using a non-linear mark function which depends on the smoothed local density, such that low- or high-density regions are up- or down-weighted. This approach may also capture information on some higher-order statistical moments of the field. We model the HI distribution using semi-numerical simulations and for the first time study the marked HI power spectrum, across $1 \leq z \leq 5$. Our analysis indicates that there is considerable evolution of the HI field during the post-reionization era. Over a wide range of length scales (comoving wave numbers $0.05\leq k \leq 1.0$ Mpc$^{-1}$) we expectedly find that the HI evolves slowly at early times, but more rapidly at late times. This evolution is not well-captured by the power spectrum of the standard (unmarked) HI field. We also study how the evolution of the HI field depends on the chosen smoothing scale for the mark, and how this affects the marked power spectrum. We conclude that the information about the HI content at low and intermediate densities is important for a correct and consistent analysis of HI content and evolution based on the 21-cm background. The marked power spectrum can thus provide a less biased statistic for parameter constraints than the normal power spectrum.

Many compact objects (black holes and neutron stars) exist in binaries. These binaries are normally discovered through their interactions, either from accretion as an X-ray binary or collisions as a gravitational wave source. However, the majority of compact objects in binaries should be non-interacting. Recently proposed discoveries have used radial velocities of a bright star (main sequence or evolved) that are indicative of a massive but dark companion, which is inferred to be a compact object. Unfortunately, this burgeoning new field has been hindered by false positives, including the ``Unicorn'' (V723 Mon) which was initially believed to be a red giant/black hole binary before being refuted. In this work, we investigate the evolution of stellar binary populations over time, using the binary evolution code COSMIC to simulate binary populations and determine the probability of a candidate object being either a ``true Unicorn'' (actual compact objects in binaries) or a false positive. We find that main sequence stars have a higher true Unicorn probability than red giants or naked helium stars (an exposed core of an evolved star), particularly if the companion is more massive and is >3 times less luminous than the MS star. We also find that a top-heavy initial mass function raises the true Unicorn probability further, that super-solar metallicity reduces the probability, and that most true Unicorns are found at periods <100 days. Finally, we find that a significant fraction of true Unicorns do not evolve into x-ray binaries during the age of the universe.

Most galaxies follow well-defined scaling relations of metallicity and stellar mass; however, some outliers at the low mass end of the observed galaxy population exhibit unusually high metallicity for their mass. Understanding how these objects get to be so metal-rich is vital for understanding the role of feedback in galaxy formation. Using the TNG50 simulation, we explore the origins of this phenomenon. We identify 227 metal-rich, Compact Stellar Systems (CSSs) that deviate significantly from this scaling relation. These CSSs are satellites located in the vicinity of massive host galaxies, with stellar masses ranging from $10^{8} M_{\odot}$ to $10^{10} M_{\odot}$ (including six systems that are close analogs of the M31-M32 system). Contrary to the previously assumed scenario that such objects are predominantly products of tidal stripping, our results suggest a more prevalent role for ram pressure in their formation. Indeed, 76\% (173) of these CSSs are formed through a burst of star formation occurring around the time of the first pericentric passage, typically at redshifts $z\lesssim1$, aided by strong ram pressure and tidal forces. The high ram pressure, resulting from the CSSs' rapid motion near the halo center, facilitates metal enrichment, producing high-metallicity CSSs by confining the metal-rich gas from bursty star formation, which leads to distinct stellar populations characterized by enhanced metallicity as well as high $\alpha$-abundance. Only the remaining 24\% (54) of metal-rich CSSs are generated through the tidal stripping of massive progenitors. Our results further indicate that M32 is more likely to have formed through intense star formation events rather than through gradual, tidal stripping, thereby providing crucial insights into the nature of low mass, compact galaxy formation.

The petaling effect, induced by pupil fragmentation from the telescope spider, drastically affects the performance of high contrast instruments by inducing core splitting on the PSF. Differential piston/tip/tilt aberrations within each optically separated fragment of the pupil are poorly measured by commonly used Adaptive Optics (AO) systems. We here pursue a design of dedicated low-order wavefront sensor -- or petalometers -- to complement the main AO. Interferometric devices sense differential aberrations between fragments with optimal sensitivity; their weakness though is their limitation to wrapped phase measurements. We show that by combining multiple spectral channels, we increase the capture range for petaling aberrations beyond several microns, enough to disambiguate one-wave wrapping errors made by the main AO system. We propose here to implement a petalometer from the multi-wavelength imaging mode of the VAMPIRES visible-light instrument, deployed on SCExAO at the Subaru Telescope. The interferometric measurements obtained in four spectral channels through a 7 hole non-redundant mask allow us to effiiently reconstruct diffierential piston between pupil petals.

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

This paper investigates the secular motion of a massless asteroid within the framework of the double-averaged elliptic restricted three-body problem. By employing Poincaré variables, we analyze the stability properties of asteroid orbits in the presence of planetary perturbations. Our study reveals that periodic orbits identified in the planar configuration maintain stability in the spatial perturbed problem across a wide range of parameter values. These findings, supported by numerical simulations, contribute to a deeper understanding of asteroid dynamics and have implications for studying exoplanetary systems with highly eccentric host stars.

The (re)analysis of data on the X-ray emitting pulsars PSR J0030+0451 and PSR J0740+6620, as well as new results on PSR J0437-4715, are confronted with the predictions of the equation of state (EoS) models allowing for strong first-order phase transition for the mass-radius ($M$-$R$) diagram. We use models that are based on a covariant density functional (CDF) EoS for nucleonic matter at low densities and a quark matter EoS, parameterized by the speed of sound, at higher densities. To account for the variations in the ellipses for PSR J0030+0451 obtained from different analyses, we examined three scenarios to assess their consistency with our models, focusing particularly on the potential formation of twin stars. We found that in two scenarios, where the ellipses for PSR J0030+0451 and PSR J0437-4715 with masses close to the canonical mass $\sim 1.4\,M_{\odot}$ are significantly separated, our models allow for the presence of twin stars as a natural explanation for potential differences in the radii of these stars.

The Small-Correlated-Against-Large Estimator (SCALE) for small-scale lensing of the cosmic microwave background (CMB) provides a novel method for measuring the amplitude of CMB lensing power without the need for reconstruction of the lensing field. In our previous study, we showed that the SCALE method can outperform existing reconstruction methods to detect the presence of lensing at small scales ($\ell \gg 3000$). Here we develop a procedure to include information from SCALE in cosmological parameter inference. We construct a precise neural network emulator to quickly map cosmological parameters to desired CMB observables such as temperature and lensing power spectra and SCALE cross spectra. We also outline a method to apply SCALE to full-sky maps of the CMB temperature field, and construct a likelihood for the application of SCALE in parameter estimation. SCALE supplements conventional observables such as the CMB power spectra and baryon acoustic oscillations in constraining parameters that are sensitive to the small-scale lensing amplitude such as the neutrino mass $m_\nu$. We show that including estimates of the small-scale lensing amplitude from SCALE in such an analysis provides enough constraining information to measure the minimum neutrino mass at $4\sigma$ significance in the scenario of minimal mass, and higher significance for higher mass. Finally, we show that SCALE will play a powerful role in constraining models of clustering that generate scale-dependent modulation to the distribution of matter and the lensing power spectrum, as predicted by models of warm or fuzzy dark matter.

The polarization of the Mg II h & k resonance lines is the result of the joint action of scattering processes and the magnetic field induced Hanle, Zeeman, and magneto-optical effects, thus holding significant potential for the diagnostic of the magnetic field in the solar chromosphere. The Chromospheric LAyer Spectro-Polarimeter sounding rocket experiment, carried out in 2019, successfully measured at each position along the 196 arcsec spectrograph slit the wavelength variation of the four Stokes parameters in the spectral region of this doublet around 280 nm, both in an active region plage and in a quiet region close to the limb. We consider some of these CLASP2 Stokes profiles and apply to them the recently-developed HanleRT Tenerife Inversion Code, which assumes a one-dimensional model atmosphere for each spatial pixel under consideration (i.e., it neglects the effects of horizontal radiative transfer). We find that the non-magnetic causes of symmetry breaking, due to the horizontal inhomogeneities and the gradients of the horizontal components of the macroscopic velocity in the solar atmosphere, have a significant impact on the linear polarization profiles. By introducing such non-magnetic causes of symmetry breaking as parameters in our inversion code, we can successfully fit the Stokes profiles and provide an estimation of the magnetic field vector. For example, in the quiet region pixels, where no circular polarization signal is detected, we find that the magnetic field strength in the upper chromosphere varies between 1 and 20 gauss.

Most star formation in the local Universe occurs in spiral galaxies, but their origin remains an unanswered question. Various theories have been proposed to explain the development of spiral arms, each predicting different spatial distributions of the interstellar medium. This study maps the star formation rate (SFR) and gas-phase metallicity of nine spiral galaxies with the TYPHOON survey to test two dominating theories: density wave theory and dynamic spiral theory. We discuss the environmental effects on our galaxies, considering reported environments and merging events. Taking advantage of the large field of view covering the entire optical disk, we quantify the fluctuation of SFR and metallicity relative to the azimuthal distance from the spiral arms. We find higher SFR and metallicity in the trailing edge of NGC~1365 (by 0.117~dex and 0.068~dex, respectively) and NGC~1566 (by 0.119~dex and 0.037~dex, respectively), which is in line with density wave theory. NGC~2442 shows a different result with higher metallicity (0.093~dex) in the leading edge, possibly attributed to an ongoing merging. The other six spiral galaxies show no statistically significant offset in SFR or metallicity, consistent with dynamic spiral theory. We also compare the behaviour of metallicity inside and outside the co-rotation radius (CR) of NGC~1365 and NGC~1566. We find comparable metallicity fluctuations near and beyond the CR of NGC~1365, indicating gravitational perturbation. NGC~1566 shows the greatest fluctuation near the CR, in line with the analytic spiral arms. Our work highlights that a combination of mechanisms explains the origin of spiral features in the local Universe.

Globular clusters offer a powerful way to test the properties of stellar populations and the late stages of low-mass stellar evolution. In this paper we study oscillating giant stars and overtone RR Lyrae-type pulsators in the nearest globular cluster, M4, with the help of high-precision, continuous light curves collected by the Kepler space telescope in the K2 mission. We determine the frequency composition of five RRc stars and model their physical parameters with a grid of linear pulsation models. We are able, for the first time, to compare seismic masses of RR Lyrae stars directly to the masses of the very similar red horizontal branch stars in the same stellar population, independently determined from asteroseismic scaling relations. We find a close match, with an average seismic mass of $0.651\pm0.028\,M_\odot$ for RR Lyrae stars and $0.657\pm0.034\,M_\odot$ for red horizontal-branch stars. While the validity of our RR Lyrae masses still relies on the similarity of neighboring horizontal branch subgroups, this result strongly indicates that RRc stars may indeed exhibit high-degree, $l = 8$ and 9 non-radial modes, and modeling these modes can provide realistic mass estimates. We also determine the He content of the cluster to be $Y = 0.266\pm 0.008$, and compare the seismic masses for our sample of RR Lyrae to theoretical mass relations and highlight the limitations of these relations.

We develop a relativistic scenario of fast magnetic reconnection process, for general magnetohydrodynamical plasmas around Kerr black holes. Generalizing the Petschek model, we study various properties of the reconnection layer in distinct configurations. When current sheet forms in the zero-angular-momentum (ZAMO) frame which corotates with the black hole, the reconnection rate for both radial and azimuthal configurations is decreased by spacetime curvature. However, when the current sheet forms in a non-ZAMO frame, which rotates either faster or slower than the black hole, detail analysis establishes that for any given slow rotations (subrelativistic at most) and mildly relativistic inflow, the ZAMO observer will find asymmetric reconnection rates for radial configuration: it is decreased on one side of the current sheet and is increased on the other side in comparison to the unrotation limit. This is valid to both the Sweet-Parker and the Petschek scenario. The results clarify the effects of rotation on the reconnection layer in the laboratory frame in the flat spacetime limit.

This study applied machine learning models to estimate stellar rotation periods from corrected light curve data obtained by the NASA Kepler mission. Traditional methods often struggle to estimate rotation periods accurately due to noise and variability in the light curve data. The workflow involved using initial period estimates from the LS-Periodogram and Transit Least Squares techniques, followed by splitting the data into training, validation, and testing sets. We employed several machine learning algorithms, including Decision Tree, Random Forest, K-Nearest Neighbors, and Gradient Boosting, and also utilized a Voting Ensemble approach to improve prediction accuracy and robustness. The analysis included data from multiple Kepler IDs, providing detailed metrics on orbital periods and planet radii. Performance evaluation showed that the Voting Ensemble model yielded the most accurate results, with an RMSE approximately 50\% lower than the Decision Tree model and 17\% better than the K-Nearest Neighbors model. The Random Forest model performed comparably to the Voting Ensemble, indicating high accuracy. In contrast, the Gradient Boosting model exhibited a worse RMSE compared to the other approaches. Comparisons of the predicted rotation periods to the photometric reference periods showed close alignment, suggesting the machine learning models achieved high prediction accuracy. The results indicate that machine learning, particularly ensemble methods, can effectively solve the problem of accurately estimating stellar rotation periods, with significant implications for advancing the study of exoplanets and stellar astrophysics.

Filaments are believed to play a key role in high-mass star formation. We present a systematic study of the filaments and their hosting clumps in the G35 molecular complex using JCMT SCUBA-2 850 $\micron$ continuum data. We identified five clouds in the complex and 91 filaments within them, some of which form 10 hub-filament systems (HFSs), each with at least 3 hub-composing filaments. We also compiled a catalogue of 350 dense clumps, 183 of which are associated with the filaments. We investigated the physical properties of the filaments and clumps, such as mass, density, and size, and their relation to star formation. We find that the global mass-length trend of the filaments is consistent with a turbulent origin, while the hub-composing filaments of high line masses ($m_{\rm l}\,>$\,230\,$\mathrm{M_{\odot}~pc^{-1}}$) in HFSs deviate from this relation, possibly due to feedback from massive star formation. We also find that the most massive and densest clumps (R\,$>$\,0.2\,pc, M\,$>35\,\mathrm{M_{\odot}}$, $\mathrm{\Sigma}>\,0.05\,\mathrm{g~cm^{-2}}$) are located in the filaments and in the hubs of HFS with the latter bearing a higher probability of occurrence of high-mass star-forming signatures, highlighting the preferential sites of HFSs for high-mass star formation. We do not find significant variation in the clump mass surface density across different evolutionary environments of the clouds, which may reflect the balance between mass accretion and stellar feedback.

This contribution summarizes the main activities and objectives of the outreach project Astroaccesible, whose main aim is to carry the teaching and diffusion of astronomy among all kinds of collectives, focusing on blind and visually impaired (BVI) people. This project is led by a blind astronomer and aims to use a variety of resources based on different sensory channels, avoiding limiting the transmission of concepts to visual perception. This principle favors inclusion and benefits everyone, as the information is not presented using just one channel. This strategy is especially convenient for the nowadays typical data acquisition, where a variety of sources of information, not solely based on the collection of different spectral domains of electromagnetic radiation, is used. Moreover, the study of new multi-messenger astronomy could be much better understood using a multi-messenger teaching approach, favoring inclusion, motivation, and creativity.

Recent measurements of the cosmic-ray electron plus positron spectrum by several experiments have confirmed the presence of a break at $\sim\,1$ TeV. The origin of the break is still not clearly understood. In this work, we explore different possibilities for the origin which include an electron source spectrum with a broken power-law, a power-law with an exponential or super-exponential cut-offs and the absence of potential nearby cosmic-ray sources. Based on the observed electron plus positron data from the DAMPE and the H.E.S.S experiments, and considering supernova remnants as the main sources of cosmic rays in the Galaxy, we find statistical evidence in favour of the scenario with a broken power-law source spectrum with the best-fit source parameters obtained as $\Gamma=2.39$ for the source spectral index, $E_0\approx 1.6$ TeV for the break energy and $f=1.59\times 10^{48}$ ergs for the amount of supernova kinetic energy injected into cosmic-ray electrons. Such a power-law break in the spectrum has been predicted for electrons confined inside supernova remnants after acceleration via diffusive shock acceleration process, and also indicated by the multi-wavelength study of supernova remnants. All these evidences have shown that the observed spectral break provides a strong indication of a direct link between cosmic-ray electrons and their sources. Our findings further show that electrons must undergo spectral changes while escaping the source region in order to reconcile the difference between the spectral index of electrons observed inside supernova remnants and that obtained from Galactic cosmic-ray propagation studies.

Accreting neutron stars (NSs) are expected to emit a redshifted 2.2 MeV line due to the capture of neutrons produced through the spallation processes of $^4$He and heavier ions in their atmospheres. Detecting this emission would offer an independent method for constraining the equation of state of NSs and provide valuable insights into nuclear reactions occurring in extreme gravitational and magnetic environments. Typically, a higher mass accretion rate is expected to result in a higher 2.2 MeV line intensity. However, when the mass accretion rate approaches the critical threshold, the accretion flow is decelerated by the radiative force, leading to a less efficient production of free neutrons and a corresponding drop in the flux of the spectral line. This makes the brightest X-ray pulsars unsuitable candidates for gamma-ray line detection. In this work, we present a theoretical framework for predicting the optimal X-ray luminosity required to detect a redshifted 2.2 MeV line in a strongly magnetized NS. As the INTEGRAL mission nears its conclusion, we have undertaken a thorough investigation of the SPI data of this line in a representative sample of accreting NSs. No redshifted 2.2 MeV line was detected. For each spectrum, we have determined the 3-sigma upper limits of the line intensity, assuming different values of the line width. Our findings suggest that advancing our understanding of the emission mechanism of the 2.2 MeV line, as well as the accretion flow responsible for it, will require a substantial increase in sensitivity from future MeV missions. For example, for a bright X-ray binary such as Sco X-1, we would need at least a 3-sigma line point source sensitivity of ~1E-6 ph/cm^2/s, that is, about two orders of magnitude better than that currently achieved. [Abridged]

We report new X-ray results from the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL), Swift, Chandra, and XMM-Newton observations of the hitherto poorly studied unidentified X-ray transient IGR J17419-2802. We studied in detail the temporal, spectral, and energetic properties of three hard X-ray outbursts detected above 20 keV by INTEGRAL. They are all characterized by an average X-ray luminosity of 3$\times$10$^{35}$~erg~s$^{-1}$ and a constrained duration of a few days. This marks a peculiarly short and faint X-ray transient nature for IGR J17419-2802. From archival unpublished soft X-ray observations, we found that the source spends most of the time undetected at very low X-ray fluxes (down to $<4.7\times10^{-14}$ erg cm $^{-2}$ s$^{-1}$) for a dynamic range >2,000 when in outburst. We provided an accurate arcsecond-sized source error circle. Inside it, we pinpointed the best candidate near-infrared counterpart whose photometric properties are compatible with a late-type spectral nature. Based on our new findings, we suggest that IGR J17419-2802 is a new member of the very faint X-ray transients (VFXTs) class. Detailed investigations of VFXT outbursts above 20 keV are particularly rare. In this respect, our reported INTEGRAL outbursts are among the best studied to date; in particular, their constrained duration of a few days is among the shortest ever measured for a VFXT.

We study the dependence of the physical and observable properties of the CGM on its halo mass. We analyse 22 simulations from the Auriga suite of high resolution cosmological `zoom-in' simulations at $z=0$ with halo masses $10^{10}~\text{M}_{\odot}\leq\text{M}_{\mathrm{200c}}\leq10^{12}~\text{M}_{\odot}$. We find a larger scatter in temperature and smaller scatter in metallicity in more massive haloes. The scatter of temperature and metallicity as a function of radius increases out to larger radii. The median and scatter of the volume-weighted density and mass-weighted radial velocity show no significant dependence on halo mass. Our results highlight that the CGM is more multiphase in haloes of higher mass. We additionally investigate column densities for HI and the metal ions CIV, OVI, MgII and SiII as a function of stellar mass and radius. We find the HI and metal ion column densities increase with stellar mass at sufficiently large radii ($R\gtrsim{0.2}$R$_{\mathrm{200c}}$). We find good agreement between our HI column densities and observations outside $20$% of the virial radius and overpredict within $20$%. MgII and SiII are similarly overpredicted within $20$% of the virial radius, but drop off steeply at larger radii. Our OVI column densities underpredict observations for stellar masses between $10^{9.7}~\text{M}_{\odot}\leq\text{M}_{\star}<10^{10.8}~\text{M}_{\odot}$ with reasonable agreement at $10^{10.8}~\text{M}_{\odot}$. CIV column densities agree with observational detections above a halo mass of $10^{9.7}~\text{M}_{\odot}$. We find that OVI (MgII) traces the highest (lowest) temperatures, and lowest (highest) density and metallicity. OVI and CIV are photo-ionized (collisionally ionized) at low (high) halo masses with a transition to higher temperatures at $10^{11}~\text{M}_{\odot}$. However, there is no clear trend for the radial velocity of the ions.

Stellar mass black holes in the disks around active galactic nuclei (AGN) are promising sources for gravitational wave detections by LIGO/VIRGO. Recent studies suggest this environment fosters the formation and merger of binary black holes. Many of these studies often assumed a simple, laminar AGN disk without magnetic fields or turbulence. In this work, we present the first 3D magnetohydrodynamical simulations of circum-single disks around isolated and binary black holes in strongly magnetized, stratified accretion disks with turbulence driven by magneto-rotational instability. We simulated three scenarios with varying initial net-vertical magnetic field strengths: weak, intermediate, and strong. Our results show that weakly magnetized models produce circum-single disks aligned with the AGN disk's equatorial plane, similar to past hydrodynamic simulations. However, intermediate and strong magnetic fields result in randomly misaligned disks, contingent upon the availability of local ambient angular momentum within turbulent regions. Our findings emphasize the significant impact of ambient gas in the AGN disk on the inclination of circum-single disks, linked to magnetically induced inhomogeneity and angular momentum during disk formation. The presence of misaligned disks, both in single and binary black hole systems, could have profound implications for the long-term evolution of black hole spin and the inclination of the disk at the horizon scale.

In this work we revisit $\alpha$-Starobinsky inflation, also know as $E$-model, in the light of current CMB and LSS observations. The inflaton potential in the Einstein frame for this model contains a parameter $\alpha$ in the exponential, which alters the predictions for the scalar and tensor power spectra of Starobinsky inflation. We obtain these power spectra numerically without using slow-roll approximation and perform MCMC analysis to put constraints on parameters $M$ and $\alpha$ from Planck-2018, BICEP/Keck (BK18) and other LSS observations. We consider general reheating scenario by varying the number of e-foldings during inflation, $N_{pivot}$, along with the other parameters. We find $\log_{10}\alpha = 0.0^{+1.6}_{-5.6}$, $\log_{10}M= -4.91^{+0.69}_{-2.7}$ and $N_{pivot} = 53.2^{+3.9}_{-5}$ with $95\%$ C. L.. This implies that the present CMB and LSS observations are insufficient to constrain the parameter $\alpha$. We also find that there is no correlation between $N_{pivot}$ and $\alpha$.

Blazars exhibit multiwavelength variability, a phenomenon whose underlying mechanisms remain elusive. This study investigates the origin of such variability through leptonic blazar emission simulations, focusing on stochastic fluctuations in environmental parameters. By analyzing the spectral indices of the power spectral densities of the variability, we assess their relationship with the underlying fluctuations. Our findings reveal that the variability spectral indices remain almost independent of the variations responsible for their emergence. This suggests a complex interplay of factors contributing to the observed multiwavelength variability in blazars.

The X-ray Integral Field Unit (X-IFU) is an instrument of ESA's future NewAthena space observatory, with the goal to provide high-energy resolution ($<$ 4 eV at X-ray energies up to 7 keV) and high-spatial resolution (9") spectroscopic imaging over the X-ray energy range from 200 eV to 12 keV, by means of an array of about 1500 transition-edge sensors (TES) read out via SQUID time-division multiplexing (TDM). In 2022, to aid in the transfer of TDM readout technology from the laboratory toward flight hardware, our team commissioned a new TDM-based laboratory test-bed at SRON. This setup hosts an array of $75\times 75\ \mu$m$^2$ TESs that are read out via 2-column $\times$ 32-row TDM. A system component that is critical to high-performance operation is the wiring harness that connects the room-temperature electronics to the cryogenic readout componentry. In November 2023, we implemented a re-designed flex harness, which in the SRON test-bed has a length close to what is envisioned for the X-IFU flight harness. We report here on our characterization of the TDM system with the new flex harness, which allowed the system to achieve a co-added energy resolution at a level of 2.7~eV FWHM at 6~keV via 32-row readout. In addition, we provide an outlook on the upcoming integration of TDM readout into the X-IFU Focal-Plane Assembly Development Model.

GHOST is a newly operational optical fiber-fed high-resolution spectrograph at the Gemini South 8.1m telescope. It currently offers the choice of two resolution modes captured by one (or two) input IFUs with a FOV of 1.2'' and a spectral resolving power of 56,000 and 76,000 for the unbinned CCDs. At the high-resolution mode, one can also instigate a simultaneous ThXe calibration lamp, which along with a simultaneous pseudo-slit profile constructed from reformatting the input IFU image will allow for precision radial velocity measurements. Here we talk about the proposed roadmap towards full queue operations, potential upgrades, and the error terms contributing to the final on-sky RV precision, which is estimated to be in the 1-10 m s$^{-1}$ range.

Magnetic wind braking drives the spin-down of low-mass stars and the evolution of most interacting binary stars. A magnetic braking prescription that was claimed to reproduce both the period distribution of cataclysmic variables (CVs) and the evolution of the rotation rates of low-mass stars is based on a relation between the angular momentum loss rate and magnetic field complexity. The magnetic braking model based on field complexity has been claimed to predict a detached phase that could explain the observed period gap in the period distribution of CVs but has never been tested in detailed models of CV evolution. Here we fill this gap. We incorporated the suggested magnetic braking law in MESA and simulated the evolution of CVs for different initial stellar masses and initial orbital periods. We find that the prescription for magnetic braking based on field complexity fails to reproduce observations of CVs. The predicted secondary star radii are smaller than measured, and an extended detached phase that is required to explain the observed period gap (a dearth of non-magnetic CVs with periods between ${\sim}2$ and ${\sim}3$ hours) is not predicted. Proposed magnetic braking prescriptions based on a relation between the angular momentum loss rate and field complexity are too weak to reproduce the bloating of donor stars in CVs derived from observations and, in contrast to previous claims, do not provide an explanation for the observed period gap. The suggested steep decrease in the angular momentum loss rate does not lead to detachment. Stronger magnetic braking prescriptions and a discontinuity at the fully convective boundary are needed to explain the evolution of close binary stars that contain compact objects. The tension between braking laws derived from the spin-down of single stars and those required to explain CVs and other close binaries containing compact objects remains.

On large scales, measurements of the Lyman-$\alpha$ forest offer insights into the expansion history of the Universe, while on small scales, these impose strict constraints on the growth history, the nature of dark matter, and the sum of neutrino masses. This work introduces ForestFlow, a cosmological emulator designed to bridge the gap between large- and small-scale Lyman-$\alpha$ forest analyses. Using conditional normalizing flows, ForestFlow emulates the 2 Lyman-$\alpha$ linear biases ($b_\delta$ and $b_\eta$) and 6 parameters describing small-scale deviations of the 3D flux power spectrum ($P_\mathrm{3D}$) from linear theory. These 8 parameters are modeled as a function of cosmology $\unicode{x2013}$ the small-scale amplitude and slope of the linear power spectrum $\unicode{x2013}$ and the physics of the intergalactic medium. Thus, in combination with a Boltzmann solver, ForestFlow can predict $P_\mathrm{3D}$ on arbitrarily large (linear) scales and the 1D flux power spectrum ($P_\mathrm{1D}$) $\unicode{x2013}$ the primary observable for small-scale analyses $\unicode{x2013}$ without the need for interpolation or extrapolation. Consequently, ForestFlow enables for the first time multiscale analyses. Trained on a suite of 30 fixed-and-paired cosmological hydrodynamical simulations spanning redshifts from $z=2$ to $4.5$, ForestFlow achieves $3$ and $1.5\%$ precision in describing $P_\mathrm{3D}$ and $P_\mathrm{1D}$ from linear scales to $k=5\,\mathrm{Mpc}^{-1}$ and $k_\parallel=4\,\mathrm{Mpc}^{-1}$, respectively. Thanks to its parameterization, the precision of the emulator is also similar for both ionization histories and two extensions to the $\Lambda$CDM model $\unicode{x2013}$ massive neutrinos and curvature $\unicode{x2013}$ not included in the training set. ForestFlow will be crucial for the cosmological analysis of Lyman-$\alpha$ forest measurements from the DESI survey.

We analyze JWST MIRI/MRS observations of the infrared PAH bands in the nuclear and circumnuclear regions of local AGN from the GATOS Survey. In this work, we examine the PAH properties in the circumnuclear regions of AGN and AGN-outflows, and compare them to those in star-forming regions and the innermost regions of AGN. This study employs 4.9-28.1 micron sub-arcsecond angular resolution data to investigate the properties of PAH in three nearby sources (DL~30-40 Mpc). Our findings align with previous JWST studies, showing that the central regions of AGN show a larger fraction of neutral PAH molecules (i.e. elevated 11.3/6.2 and 11.3/7.7 PAH ratios) compared to star-forming galaxies. We find that the AGN might affect not only the PAH population in the innermost region but also in the extended regions up to ~kpc scales. By comparing our observations to PAH diagnostic diagrams, we find that, in general, regions located in the projected direction of the AGN-outflow occupy similar positions on the PAH diagnostic diagrams as those of the innermost regions of AGN. Star-forming regions that are not affected by the AGN in these galaxies share the same part of the diagram as Star-forming galaxies. We examine the potential of the PAH-H2 diagram to disentangle AGN versus star-forming activity. Our results suggest that in Sy-like AGN, illumination and feedback from the AGN might affect the PAH population at nuclear and kpc scales, in particular, the ionization state of the PAH grains. However, PAH sizes are rather similar. The carriers of the ionized PAH bands (6.2 and 7.7 micron) are less resilience than those of neutral PAH bands (11.3 micron), which might be particularly important for strongly AGN-host coupled systems. Therefore, caution must be applied when using PAH bands as star-formation rate indicators in these systems even at kpc scales, with the ionized ones being more affected by the AGN.

Marked power spectra provide a computationally efficient way to extract non-Gaussian information from the matter density field using the usual analysis tools developed for the power spectrum without the need for explicit calculation of higher-order correlators. In this work, we explore the optimal form of the mark function used for re-weighting the density field, to maximally constrain cosmology. We show that adding to the mark function or multiplying it by a constant leads to no additional information gain, which significantly reduces our search space for optimal marks. We quantify the information gain of this optimal function and compare it against mark functions previously proposed in the literature. We find that we can gain around $\sim2$ times smaller errors in $\sigma_8$ and $\sim4$ times smaller errors in $\Omega_m$ compared to using the traditional power spectrum alone, an improvement of $\sim60\%$ compared to other proposed marks when applied to the same dataset.

Aim: The cold molecular gas mass is one of the crucial, yet challenging parameters in galaxy evolution studies. Here, we introduce a new calibration for estimating molecular gas masses using mid-infrared (MIR) photometry. This topic is timely, as JWST now allows us to detect the MIR emission of typical main-sequence galaxies across a wide range of masses and star formation rates with modest time investments. This Letter highlights the strong synergy between ALMA and JWST for studies of dust and gas at cosmic noon. Methods: We combine a sample of 14 main sequence galaxies at z=1-3 with robust CO detections and multi-band MIR photometry, along with a literature sample at z=0-4 with CO and PAH spectroscopy, to study the relationship between PAH, CO(1-0), and total IR luminosities. PAH luminosities are derived from modeling rest-frame UV to sub-mm data. The new z=1-3 sample extends previous high-z studies to about an order-of-magnitude lower PAH and CO luminosities, into the regime of local starbursts for the first time. Results: The PAH-to-CO luminosity ratio remains constant across a wide range of luminosities, for various galaxy types, and throughout the explored redshift range. In contrast, the PAH-to-IR and CO-to-IR luminosity ratios deviate from a constant value at high L(IR). The intrinsic scatter in the L(PAH)-L'(CO) relation is 0.21 dex, with a median of 1.40, and a power-law slope of $1.07 \pm 0.04$. Both the PAH-IR and CO-IR relations are sub-linear. Given the tight and uniform PAH-CO relation over ~3 orders of magnitude, we provide a recipe to estimate the cold molecular gas mass of galaxies from PAH luminosities, with a PAH-to-molecular gas conversion factor of $\alpha_{\rm PAH7.7} = (3.08 \pm 1.08)(4.3/\alpha_{\rm CO})\,M_{\odot}/L_{\odot}$. This method opens a new window to explore the gas content of galaxies beyond the local Universe using multi-wavelength JWST/MIRI imaging.

46P/Wirtanen is a Jupiter-family comet, probably originating from the Solar System's Kuiper belt, that now resides on a 5.4 year elliptical orbit. During its 2018 apparition, comet 46P passed unusually close to the Earth (within 0.08 au), presenting an outstanding opportunity for close-up observations of its inner coma. Here we present observations of HCN, H$^{13}$CN and HC$^{15}$N emission from 46P using the Atacama Compact Array (ACA). The data were analyzed using the SUBLIME non-LTE radiative transfer code to derive $^{12}$C/$^{13}$C and $^{14}$N/$^{15}$N ratios. The HCN/H$^{13}$CN ratio is found to be consistent with a lack of significant $^{13}$C fractionation, whereas the HCN/HC$^{15}$N ratio of $68\pm27$ (using our most conservative $1\sigma$ uncertainties), indicates a strong enhancement in $^{15}$N compared with the solar and terrestrial values. The observed $^{14}$N/$^{15}$N ratio is also significantly lower than the values of $\sim140$ found in previous comets, implying a strong $^{15}$N enrichment in 46P's HCN. This indicates that the nitrogen in Jupiter-family comets could reach larger isotopic enrichments than previously thought, with implications for the diversity of $^{14}$N/$^{15}$N ratios imprinted into icy bodies at the birth of the Solar System.

Recent large-scale structure (LSS) surveys have revealed a persistent tension in the value of $S_8$ compared to predictions from the standard cosmological model. This tension may suggest the need for new physics beyond the standard model, but an accurate characterisation of baryonic effects is essential to avoid biases. Although some studies indicate that baryonic effects are too small to resolve this tension, others propose that more aggressive feedback mechanisms could reconcile differences between cosmic microwave background (CMB) measurements and low-redshift LSS observations. In this paper, we investigate the role of baryonic effects in alleviating the $S_8$ tension. We extend the SP(k) model (Salcido et al. 2023), which was trained on hundreds of cosmological hydrodynamical simulations to map the suppression of the matter power spectrum to the baryon fraction in groups and clusters, to predict the required baryon fraction for a given $P(k)$ suppression. We then compare predictions from recent cosmic shear (weak lensing) analyses with the latest baryon budget measurements from X-ray and weak gravitational lensing studies. Our findings show that studies marginalising over baryonic effects while fixing cosmological parameters to a Planck-like cosmology predict strong $P(k)$ suppression and baryon fractions that are much lower than existing low-redshift baryon budget estimates of galaxy groups and clusters. Conversely, most studies that marginalise over both cosmological parameters and baryonic effects imply baryon fractions that are consistent with observations but lower values of $S_8$ than inferred from the CMB. Unless the observed baryon fractions are biased high by a factor of several, these results suggest that a mechanism beyond baryonic physics alone is required to modify or slow down the growth of structure in the universe in order to resolve the $S_8$ tension.

We report the discovery of the Fe K line emission at $\sim6.62^{+0.06}_{-0.06}$ keV with a width of $\sim0.19^{+0.05}_{-0.05}$ keV using two epochs of {\it Chandra} archival data from the nucleus of the galaxy 4C+37.11, which is known to host a binary supermassive black hole (BSMBH) system where the SMBHs are separated by $\sim7$ mas or $\sim$ 7pc. Our study reports the first detection of the Fe K line from a known binary AGN, and has an F-statistic value of 20.98 and probability $2.47\times 10^{-12}$. Stacking of two spectra reveals another Fe K line component at $\sim7.87^{+0.19}_{-0.09}$ keV. Different model scenarios indicate that the lines originate from the combined effects of accretion disk emission and circumnuclear collisionally ionized medium. The observed low column density favors the gas-poor merger scenario, where the high temperature of the hot ionized medium may be associated with the shocked gas in the binary merger and not with star formation activity. The estimated total BSMBH mass and disk inclination are $\sim1.5\times10^{10}$ M$_\odot$ and $\gtrsim75^\circ$, indicating that the BSMBH is probably a high inclination system. The spin parameter could not be tightly constrained from the present data sets. Our results draw attention to the fact that detecting the Fe K line emissions from BSMBHs is important for estimating the individual SMBH masses, and the spins of the binary SMBHs, as well as exploring their emission regions.

This study introduces PI-AstroDeconv, a physics-informed semi-supervised learning method specifically designed for removing beam effects in astronomical telescope observation systems. The method utilizes an encoder-decoder network architecture and combines the telescope's point spread function or beam as prior information, while integrating fast Fourier transform accelerated convolution techniques into the deep learning network. This enables effective removal of beam effects from astronomical observation images. PI-AstroDeconv can handle multiple PSFs or beams, tolerate imprecise measurements to some extent, and significantly improve the efficiency and accuracy of image deconvolution. Therefore, this algorithm is particularly suitable for astronomical data processing that does not rely on annotated data. To validate the reliability of the algorithm, we used the SKA Science Data Challenge 3a datasets and compared it with the CLEAN deconvolution method at the 2-D matter power spectrum level. The results demonstrate that our algorithm not only restores details and reduces blurriness in celestial images at the pixel level but also more accurately recovers the true neutral hydrogen power spectrum at the matter power spectrum level.

The accretion flow onto a neutron star will be impacted due to irradiation by a Type I X-ray burst. The burst radiation exerts Poynting-Robertson (PR) drag on the accretion disk, leading to an enhanced mass accretion rate. Observations of X-ray bursts often find evidence that the normalization of the disk-generated persistent emission (commonly denoted by the factor $f_a$) increases during a burst, and changes in $f_a$ have been used to infer the evolution in the mass accretion rate due to PR drag. Here, we examine this proposed relationship between $f_a$ and mass accretion rate enhancement using time-resolved data from simulations of accretion disks impacted by Type I X-ray bursts. We consider bursts from both spinning and non-spinning neutron stars and track both the change in accretion rate due to PR grad and the disk emission spectra during the burst. Regardless of the neutron star spin, we find that $f_a$ strongly correlates with the disk temperature and only weakly follows the mass accretion rate (the Pearson correlation coefficients are $\leq 0.63$ in the latter case). Additionally, heating causes the disk to emit at higher energies, reducing its contribution to a soft excess. We conclude that $f_a$ cannot accurately capture the mass accretion rate enhancement and is rather a tracer of the disk temperature.

Conventional planet formation theories predict a paucity of massive planets around small stars, especially very low-mass ($0.1 - 0.3 \ M_{\odot}$) mid-to-late M dwarfs. Such tiny stars are expected to form planets of terrestrial sizes, but not much bigger. However, this expectation is challenged by the recent discovery of LHS 3154 b, a planet with period of 3.7 days and minimum mass of $13.2 \ M_{\oplus}$ orbiting a $0.11 \ M_{\odot}$ star. Here, we propose that close-in Neptune-mass planets like LHS 3154 b formed through an anomalous series of mergers from a primordial compact system of super-Earths. We perform simulations within the context of the "breaking the chains" scenario, in which super-Earths initially form in tightly-spaced chains of mean-motion resonances before experiencing dynamical instabilities and collisions. Planets as massive and close-in as LHS 3154 b ($M_p \sim 12 - 20 \ M_{\oplus}$, $P < 7$ days) are produced in $\sim$1% of simulated systems, in broad agreement with their low observed occurrence. These results suggest that such planets do not require particularly unusual formation conditions but rather are an occasional byproduct of a process that is already theorized to explain compact multi-planet systems. Interestingly, our simulated systems with LHS 3154 b-like planets also contain smaller planets at around $\sim 30$ days, offering a possible test of this hypothesis.

NEO Surveyor will detect asteroids and comets using mid-infrared thermal emission, however ground-based followup resources will require knowledge of the expected visible light brightness in order to plan characterization observations. Here we describe the range of visual-to-infrared colors that the NEOs detected by Surveyor will span, and demonstrate that for objects that have no previously reported Visual band observations, estimates of the Johnson Visual-band brightness based on infrared flux alone will have significant uncertainty. Incidental or targeted photometric followup of objects discovered by Surveyor enables predictions of the fraction of reflected light visible and near-infrared wavelengths, supporting additional detailed characterization.

Context. Inferring spectral parameters from X-ray data is one of the cornerstones of high-energy astrophysics, and is achieved using software stacks that have been developed over the last twenty years and more. However, as models get more complex and spectra reach higher resolutions, these established software solutions become more feature-heavy, difficult to maintain and less efficient. Aims. We present jaxspec, a Python package for performing this task quickly and robustly in a fully Bayesian framework. Based on the JAX ecosystem, jaxspec allows the generation of differentiable likelihood functions compilable on core or graphical process units (resp. CPU and GPU), enabling the use of robust algorithms for Bayesian inference. Methods. We demonstrate the effectiveness of jaxspec samplers, in particular the No U-Turn Sampler, using a composite model and comparing what we obtain with the existing frameworks. We also demonstrate its ability to process high-resolution spectroscopy data and using original methods, by reproducing the results of the Hitomi collaboration on the Perseus cluster, while solving the inference problem using variational inference on a GPU. Results. We obtain identical results when compared to other softwares and approaches, meaning that jaxspec provides reliable results while being $\sim 10$ times faster than existing alternatives. In addition, we show that variational inference can produce convincing results even on high-resolution data in less than 10 minutes on a GPU. Conclusions. With this package, we aim to pursue the goal of opening up X-ray spectroscopy to the existing ecosystem of machine learning and Bayesian inference, enabling researchers to apply new methods to solve increasingly complex problems in the best possible way. Our long-term ambition is the scientific exploitation of the data from the newAthena X-ray Integral Field Unit (X-IFU).

The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon-oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. We perform three-dimensional, smoothed-particle-hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state (EOS), with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS to the critical mass for gravitational collapse into a black hole (BH) and its occurrence time, $t_\textrm{col}$. We include a non-zero angular momentum of the NSs and find that $t_\textrm{col}$ ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about five minutes orbital period), promptly form a BH and release $\gtrsim 10^{52}$ erg. They form NS-BH binaries with tens of kyr merger timescale by gravitational-wave emission. BdHNe II and III do not form BHs, release $\sim 10^{50}$-$10^{52}$ erg and $\lesssim 10^{50}$ erg. They form NS-NS with a wider range of merger timescales. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a non-rotating NS. Magnetic braking by a $10^{13}$ G field can delay BH formation, leading to BH-BH or NS-BH of tens of kyr merger timescale.

High-resolution solar observations have revealed the existence of small-scale vortices, as seen in chromospheric intensity maps and velocity diagnostics. Frequently, these vortices have been observed near magnetic flux concentrations, indicating a link between swirls and the evolution of the small-scale magnetic fields. Vortices have also been studied with magneto-hydrodynamic (MHD) numerical simulations of the solar atmosphere, revealing their complexity, dynamics, and magnetic nature. In particular, it has been proposed that a rotating magnetic field structure driven by a photospheric vortex flow at its footprint produces the chromospheric swirling plasma motion. We present a complete and comprehensive description of the time evolution of a small-scale magnetic flux concentration interacting with the intergranular vortex flow and affected by processes of intensification and weakening of its magnetic field. In addition, we study the chromospheric dynamics associated with the interaction, including the analysis of a chromospheric swirl and an impulsive chromospheric jet.

Upcoming large multiwavelength photometric surveys will provide a leap in our understanding of small body populations, among other fields of modern astrophysics. Serendipitous observations of small bodies in different orbital locations allow us to study diverse phenomena related to how their surfaces scatter solar light. In particular, multiple observations of the same object in different epochs permit us to construct their phase curves to obtain absolute magnitudes and phase coefficients. In this work, we tackle a series of long-used relationships associating these phase coefficients with the taxa of small bodies and suggest that some may need to be revised in the light of large-number statistics.

The Roman Coronagraph Instrument will be the first space facility equipped with deformable mirrors (DMs). These will lead to reach a contrast of $10^{-8}$ or better in a dark hole between $3-9 \lambda/D$. Post-processing techniques play an important role in increasing the contrast limits. Our work investigates how DMs can be used to calibrate the instrument response to controlled wavefront error maps and to improve the post-processing performance. To this goal, we are developing a simulation pipeline, CAPyBARA, that includes both a propagation model of the Coronagraph and a post-processing module and produces starlight subtracted images of a science target. This pipeline will allow us to investigate alternative observing strategies and test their performance for the Roman Coronagraph. Here we present the first version of the simulator: it currently reproduces the optical propagation, which consists in the hybrid Lyot coronagraph optical structure and dark-hole digging technique (Electric Field Conjugation coupled with $\beta$-bumping), the environment (quasi-static aberration) and the post-processing. With it, we mimic a Coronagraph Instrument observing sequence, which consists in first acquiring reference star data before slewing to the scientific target, and we investigate how the evolution of quasi-static aberrations deteriorate the contrast limit in the dark hole. We simulate a science target with planets at high contrast with their star and we perform a first post-processing analysis with classical subtraction techniques. Here we present the CAPyBARA simulator, as well as some first results. The next step will be to generate PSF libraries by injecting pre-calibrated probes on the DMs (in open loop) during the reference star acquisition and compute a PCA model. Later, we will compare the performance gain obtained with the modulated-DM reference library over standard approaches (RDI).

Comet C/2017 K2 (Pan-STARRS) provided a rare opportunity to investigate the evolution of coma composition and outgassing patterns over a transitional heliocentric distance (Rh) range where activity drivers in comets are thought to change from "hypervolatile" (CO, CH$_4$, C$_2$H$_6$, and/or CO$_2$)-dominated to H2O-dominated. We performed high-resolution, cross-dispersed, near-infrared spectroscopy of C/2017 K2 with iSHELL at the NASA Infrared Telescope Facility (IRTF) and NIRSPEC at the W. M. Keck 2 Observatory. We report gas rotational temperatures (Trot) and molecular production rates (Q; mol/s) or upper limits for the hypervolatile species CH$_4$, CO, and C$_2$H$_6$, together with less volatile ices (CH$_3$OH, H$_2$O, HCN, C$_2$H$_2$, NH$_3$, and OCS) over a range of pre-perihelion distances, Rh= 3.15 - 2.35 au. We also report (or stringently constrain) abundance ratios (mixing ratios) of the targeted species with respect to CO, C$_2$H$_6$, and (when detected) H$_2$O. All volatiles were enriched relative to water in C/2017 K2 when compared to their mean values among Oort Cloud comets, whereas abundances relative to C2H6 were consistent with their average values from other long-period comets.

We present the luminosity-halo mass relations of satellite (sLHMRs) galaxies in the SDSS redMaPPer cluster catalogue and the effects of the dense cluster environment on subhalo mass evolution. We use data from the Subaru Hyper Suprime-Cam survey Year-3 catalogue of galaxy shapes to measure the weak lensing signal around these satellites. This signal serves as a probe of the matter distribution around the satellites, thereby providing the masses of their associated subhalos. We bin our satellites based on physical observable quantities such as their luminosity or the host cluster's richness, combined with their cluster-centric radial separations. Our results indicate that although more luminous satellites tend to reside in more massive halos, the sLHMRs depend on the distance of the satellite from the cluster centre. Subhalos near the cluster centre (within $<0.3 h^{-1}Mpc$) are stripped of mass. Consequently, the ratio of subhalo mass to luminosity decreases near the cluster centre. For low luminosity galaxies ($L < 10^{10} h^{-2}L_{\odot}$), the lack of evidence of increasing subhalo masses with luminosity shows the impact of tidal stripping. We also present stellar-to-subhalo mass relations (sSHMRs) for our satellite sample evolving at different cluster-centric separations. Inferred sSHMRs in the outer radial bin appear to match that observed for the field galaxies. We show that the sSHMRs from the mock-redMaPPer run on galaxy catalogues generated by the empirical UniverseMachine galaxy formation model are in good agreement with our observational results. Satellites, when binned based on the host cluster's richness, show very little dependence of the subhalo mass on the richness.

Directly imaging Earth-like exoplanets within habitable zones is challenging because faint signals can be obscured by exozodiacal dust, analogous to our solar system's zodiacal dust. This dust scatters starlight, creating a bright background noise. This paper introduces Toy Coronagraph, a Python package designed to quantify the impact of this dust on exoplanet detection. It takes circularly symmetric disk images point spread functions (PSFs), and exoplanet orbital parameters as input, generating key metrics like contrast curves, signal-to-noise ratios, and dynamic visualizations of exoplanet motion under the dust background. The package also provides tools for generating vortex coronagraph PSFs and includes example disk images. Toy Coronagraph empowers researchers to understand exozodiacal dust, develop mitigation strategies, and optimize future telescope designs and mission time, ultimately advancing the search for potentially habitable worlds. Future work will focus on handling non-circularly symmetric inputs, incorporating realistic noise models, and estimating exoplanet yield rates for future space telescope missions.

The radial distribution of gas within galactic haloes is connected to the star formation rate and the nature of baryon-driven feedback processes. Using six variants of the hydrodynamic simulation Simba, we study the impact of different stellar/AGN feedback prescriptions on the gas density profiles of haloes in the total mass range $10^{11} \, \mathrm{M}_{\odot} < M_{\mathrm{200c}} < 10^{14} \, \mathrm{M}_{\odot}$ and redshift interval $0<z<4$. We find that the radial profiles are well represented by a power law and that, for a fixed total halo mass, the slope and amplitude of such power law are generally weakly dependent on redshift. Once AGN-driven jets are activated in the simulation, the gas density profile of haloes with $M_{\rm 200c} \gtrsim 10^{13} \, \rm M_{\odot}$ declines more gently with radial distance. We argue that this distinctive feature could be exploited with current observations to discriminate amongst the predictions of the different feedback models. We introduce a universal fitting formula for the slope and amplitude of the gas density profile as a function of total halo mass and redshift. The best-fit functions are suitable for all feedback variants considered, and their predictions are in excellent agreement with the numerical results. We provide the values of all fit parameters, making our fitting formula a versatile tool to mimic the effect of Simba feedback models onto N-body simulations and semi-analytical models of galaxy formation. Our results can also aid observational estimates of the gas mass within haloes that assume a specific slope for the underlying gas density profile.

The preparation of a telescope observation time proposal is a recurring activity in observational astronomy. It is a necessary investment of time and effort to obtain data to advance a research topic. Therefore, the success of an observation proposal is a condition for progress in observational research. This guide was created to offer a straightforward, practical, and comprehensive resource for applicants who are preparing a Gemini Observatory proposal. It reviews the fundamentals of an observation proposal, including its content and evaluation criteria, to help applicants organize their submissions effectively and improve specific aspects of their presentations. This manuscript is inspired by the recommendations of the User Advisory Council established in the documents ''Criterios de evaluación de propuestas por parte del NTAC'' (National Time Allocation Committee), ''Consideraciones básicas para la presentación de propuestas Gemini'', and in experiences of members of the NTAC.

In this study, we introduce a novel approach aimed at addressing the longstanding baryon-anti-baryon asymmetry conundrum. Our proposed mechanism suggests that baryon numbers were generated during the inflationary epoch through the dynamics of the inflaton field coupled with an explicit baryon number violating interaction. Notably, during inflation, it is possible to halt the baryon number generation process via a symmetry restoration phase transition. We elucidate that prior to this phase transition, baryon numbers could be synthesized and preserved within classical field configurations. Subsequently, following the phase transition, these baryon numbers were liberated as particles. Crucially, we demonstrate that this mechanism of baryon number production is intricately linked with significant cosmological collider signals and gravitational wave (GW) signals, offering a compelling framework to explore the origins of the universe's matter-antimatter asymmetry.

We measure the speed of light with current observations, such as Type Ia Supernova, galaxy ages, radial BAO mode, as well as simulations of future redshift surveys and gravitational waves as standard sirens. By means of a Gaussian Process reconstruction, we find that the precision of such measurements can be improved from roughly 6\% to 1.5-2\%, in light of these forthcoming observations. This result demonstrates that we will be able to perform a cosmological measurement of a fundamental physical constant with unprecedented precision, which will help us underpinning if its value is truly consistent with local measurements, as predicted by the standard model of Cosmology.

The Gemini South telescope is now equipped with a new high-resolution spectrograph called GHOST (the Gemini High-resolution Optical SpecTrograph). This instrument provides high-efficiency, high-resolution spectra covering 347-1060 nm in a single exposure of either one or two targets simultaneously, along with precision radial velocity spectroscopy utilizing an internal calibration source. It can operate at a spectral element resolving power of either 76000 or 56000, and can reach a SNR$\sim$5 in a 1hr exposure on a V$\sim$20.8 mag target in median site seeing, and dark skies (per resolution element). GHOST was installed on-site in June 2022, and we report performance after full integration to queue operations in November 2023, in addition to scientific results enabled by the integration observing runs. These results demonstrate the ability to observe a wide variety of bright and faint targets with high efficiency and precision. With GHOST, new avenues to explore high-resolution spectroscopy have opened up to the astronomical community. These are described, along with the planned and potential upgrades to the instrument.