Papers for Friday, Jun 07 2024

GWnext 2024 was a meeting held in the Kavli Institute for Astronomy and Astrophysics at Peking University in March $4^\text{th} - 8^\text{th}$, 2024. In the meeting researchers at different career stages -- with a particular focus on early career scientists -- working on the different aspects of gravitational wave (GW) astronomy gathered to discuss the current status as well as prospects of the field. The meeting was divided into three core sessions: Astrophysics, GW Theory, and Detection. Each session consisted of introductory talks and extended discussion sessions. Moreover, there was a poster session where students could present their results. In this paper, we summarize the results presented during the meeting and present the most important outcomes.

We study the emission of gravitational waves produced by the magnetosphere of magnetars. We argue that several features in the spectrum could facilitate the identification of that source. In addition, in cases of extremely large magnetic fields we demonstrate that this emission can make the braking index of such stars to be well over 3, which is the standard prediction of the magnetic dipole radiation and aligned rotator mechanisms. A similar picture arises if one focuses on the second braking index. Moreover the braking index depends on both the rotational frequency and the strength of the magnetic field in striking difference from the other mechanisms. We also show that gravitational waves can be produced by polar gap regions due to their rapid charge-discharge process that takes place in timescales from nanosec to microsec. This can provide an alternative way to probe magnetars and test the polar gap model.

The Mid-Infrared Instrument (MIRI) for JWST is supplied with a suite of imaging bandpass filters optimized for full spectral coverage in eight intermediate-width bands from 5 to 26 microns and a narrower one at 11.3 microns. This contrasts with previous infrared space telescopes, which generally have provided only two broad bands, one near 10 microns and the other near 20 microns. The expanded MIRI spectral capability provides new possibilities for detailed interpretation of survey results. This is an important feature of the instrument, on top of its great increase in sensitivity and angular resolution over any previous mission. The Systematic Mid-infrared Instrument Legacy Extragalactic Survey (SMILES) was designed to take full advantage of this capability. This paper briefly describes the history of infrared surveys that paved the way for MIRI on JWST and for our approach to designng SMILES. It illustrates the use of the observations for a broad range of science programs, and concludes with a brief summary of the need for additional surveys with JWST/MIRI.

We investigate the contribution of host galaxies to the overall Dispersion Measures (DMs) for Fast Radio Bursts (FRBs) using the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite. We calculate DMs from every star particle in the simulated L* galaxies by ray-tracing through their multi-phase interstellar medium (ISM), summing the line-of-sight free thermal electron column for all gas elements within $\pm$20 kpc of the galactic mid-plane. At $z=0$, we find average (median) host-galaxy DMs of 74 (43) and 210 (94) pc cm$^{-3}$ for older ($\gtrsim$10 Myr) and younger ($\lesssim$10 Myr) stellar populations, respectively. Inclination raises the median DM measured for older populations ($\gtrsim$10 Myr) in the simulations by a factor of $\sim$2, but generally does not affect the younger stars deeply embedded in H{\small II} regions except in extreme edge-on cases (inclination $\gtrsim 85^\circ$). In kinematically disturbed snapshots ($z = 1$ in FIRE), the average (median) host-galaxy DMs are higher: 80 (107) and 266 (795) pc cm$^{-3}$ for older ($\gtrsim$10 Myr) and younger ($\lesssim$10 Myr) stellar populations, respectively. FIRE galaxies tend to have higher DM values than cosmological simulations such as IllustrisTNG. As a result, FRB host galaxies may be closer (lower redshift) than previously inferred. Furthermore, constraining host-galaxy DM distributions may help significantly constrain FRB progenitor models.

The atmospheric response for MeV $\gamma$ rays ($\sim$ 0.1 $-$ 10 MeV) can be characterized in terms of two observed components. The first component is due to photons that reach the detector without scattering. The second component is due to photons that reach the detector after scattering one or more times. While the former can be determined in a straightforward manner, the latter is much more complex to quantify, as it requires tracking the transport of all source photons that are incident on Earth's atmosphere. The scattered component can cause a significant energy-dependent distortion in the measured spectrum, which is important to account for when making balloon-borne observations. In this work we simulate the full response for $\gamma$-ray transport in the atmosphere. We find that the scattered component becomes increasingly more significant towards lower energies, and at 0.1 MeV it may increase the measured flux by as much as a factor of $\sim2-4$, depending on the photon index and off-axis angle of the source. This is particularly important for diffuse sources, whereas the effect from scattering can be significantly reduced for point sources observed with an imaging telescope.

We present an investigation of clustered stellar feedback in the form of superbubbles identified within eleven galaxies from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, at both cosmic noon (1 < z < 3) and in the local Universe. We study the spatially-resolved multiphase outflows that these supernovae drive, comparing our findings with recent theory and observations. These simulations consist of five LMC-mass galaxies and six Milky Way-mass progenitors (with a minimum baryonic particle mass of $m_{b.min} = 7100 M_{\odot}$), for which we calculate the local mass and energy loading factors on 750~pc scales from the identified outflows. We also characterize the multiphase morphology and properties of the identified superbubbles, including the `shell' of cool ($T<10^5$ K) gas and break out of energetic hot ($T>10^5$ K) gas when the shell bursts. For all galaxies, the outflow mass, momentum, and energy fluxes appear to reach their peak during the identified superbubbles, and we investigate the effects on the interstellar medium (ISM), circumgalactic medium (CGM), and subsequent star formation rates. We find that these simulations, regardless of redshift, have mass-loading factors and momentum fluxes in the cool gas that largely agree with recent observations. Lastly, we also investigate how methodological choices in measuring outflows can affect loading factors for galactic winds.

We validate the Bond et al. (2010) kinematic models for the Milky Way's disk and halo stars with Gaia Data Release 3 data. Bond et al. constructed models for stellar velocity distributions using stellar radial velocities measured by the Sloan Digital Sky Survey (SDSS) and stellar proper motions derived from SDSS and the Palomar Observatory Sky Survey astrometric measurements. These models describe velocity distributions as functions of position in the Galaxy, with separate models for disk and halo stars that were labeled using SDSS photometric and spectroscopic metallicity measurements. We find that the Bond et al. model predictions are in good agreement with recent measurements of stellar radial velocities and proper motions by the Gaia survey. In particular, the model accurately predicts the skewed non-Gaussian distribution of rotational velocity for disk stars and its vertical gradient, as well as the dispersions for all three velocity components. Additionally, the spatial invariance of velocity ellipsoid for halo stars when expressed in spherical coordinates is also confirmed by Gaia data at galacto-centric radial distances of up to 15 kpc.

The diffusion of cosmic rays (CRs) in turbulent magnetic fields is fundamental to understand various astrophysical processes. We explore the CR diffusion in the magnetic luctuations amplified by the nonlinear turbulent dynamo, in the absence of a strong mean magnetic field. Using test particle simulations, we identify three distinct CR diffusion regimes: mirroring, wandering, and magnetic moment scattering (MMS). With highly inhomogeneous distribution of the dynamo-amplified magnetic fields, we find that the diffusion of CRs is also spatially inhomogeneous. Our results reveal that lower-energy CRs preferentially undergo the mirror and wandering diffusion in the strong-field regions, and the MMS diffusion in the weak-field regions. The former two diffusion mechanisms play a more important role toward lower CR energies, resulting in a relatively weak energy dependence of the overall CR mean free path. In contrast, higher-energy CRs predominantly undergo the MMS diffusion, for which the incomplete particle gyration in strong fields has a more significant effect than the nonresonant scattering by small-scale field tangling/reversal. Compared with lower-energy CRs, they are more poorly confined in space, and their mean free paths have a stronger energy dependence. We stress the fundamental role of magnetic field inhomogeneity of nonlinear turbulent dynamo in causing the different diffusion behavior of CRs compared to that in sub-Alfvénic MHD turbulence.

We present an analysis of the JWST NIRCam and MIRI morphological properties of 80 massive ($\log_{10}(M_\ast[M_{\odot}])$=11.2$\pm$0.1) dusty star-forming galaxies at $z$$=$2.7$^{+1.2}_{-0.7}$, identified as sub-millimetre galaxies (SMGs) by ALMA, that have been observed as part of the JWST PRIMER project. To compare the structure of these massive, active galaxies to more typical star-forming galaxies, we define a sample of 850 field galaxies with matched redshifts and specific star formation rates. We identify 20$\pm$5% of the SMGs as candidate late-stage major mergers, a further 40$\pm$10% as potential minor mergers and 40$\pm$10% which have comparatively undisturbed disk-like morphologies, with no obvious massive neighbours on $\lesssim$20-30kpc (projected) scales. These rates are comparable to those for the field sample and indicate that the majority of the sub-millimetre-detected galaxies are not late-stage major mergers, but have interaction rates similar to the less-active population at $z$$\sim$2-3. Through a multi-wavelength morphological analysis, we establish that SMGs have comparable near-infrared sizes to the less active population, but exhibit lower Sérsic indices, consistent with bulge-less disks and have more structured morphologies at 2$\mu$m relative to 4$\mu$m. We find evidence for dust reddening as the origin of the morphological differences between the populations, identifying a strong correlation between the F200W$-$F444W pixel colour and the 870$\mu$m surface brightness. We conclude that SMGs and less active galaxies at the same epochs share a common disk-like structure, but the weaker bulge components of the SMGs results in a lower dynamical stability. Consequently, instabilities triggered either secularly or by minor external perturbations result in higher levels of activity (and dust content) in SMGs compared to typical star-forming galaxies. [Abridged]

The Forschungszentrum Juelich has been hosting the German part of the LOFAR archive since 2013. It is Germany's most extensive radio astronomy archive, currently storing nearly 22 petabytes (PB) of data. Future radio telescopes are expected to require a dramatic increase in long-term data storage. Here, we take stock of the current data management of the Juelich LOFAR Data Archive, describe the ingestion, the storage system, the export to the long-term archive, and the request chain. We analysed the data availability over the last 10 years and searched for the underlying data access pattern and the energy consumption of the process. We determine hardware-related limiting factors, such as network bandwidth and cache pool availability and performance, and software aspects, e.g. workflow adjustment and parameter tuning, as the main data storage bottlenecks. By contrast, the challenge in providing the data from the archive for the users lies in retrieving the data from the tape archive and staging them. Building on this analysis, we suggest how to avoid/mitigate these problems in the future and define the requirements for future even more extensive long-term data archives.

Recent stellar chemical abundance measurements of a handful of $z\sim2$ quiescent galaxies have suggested these galaxies exhibit a remarkably strong $\alpha$-enhancement compared to their local and intermediate redshift counterparts. This apparent chemical evolution following quenching suggests that even the innermost regions of massive early-type galaxies may have experienced substantial mixing of stars in mergers, challenging a purely inside-out growth model. However, larger samples are needed to determine whether a high $\alpha$-enhancement ([Mg/Fe] $\approx 0.5$) is common in $z \sim 2$ quiescent galaxies, and a comparative analysis is needed to determine whether it is consistently inferred using different stellar population synthesis models. We report age and stellar chemical abundance measurements for a sample of four gravitationally lensed quiescent galaxies at $z\sim2.1-2.65$ based on Magellan/FIRE spectroscopy. For three of these galaxies we constrain the $\alpha$-enhancement, and in two cases we measure high values comparable to earlier results when the spectra are analyzed consistently. We also find that the choice of modeling approach can exert a significant effect on the measured abundances. This model dependence can be partly, but not entirely, explained by the complex abundance patterns of $\alpha$ elements in galaxies, which has been observed at lower redshifts and in one $z \sim 2$ quiescent galaxy. Our investigation highlights the importance of independently varying abundance of $\alpha$ elements when fitting the spectra of such galaxies. Observations with JWST will soon deliver precise and spatially resolved abundances of these and other quiescent galaxies at cosmic noon, opening a new window into their evolution.

We investigate the prospects for detecting and constraining density and temperature inhomogeneities in the circumgalactic medium (CGM) using absorption measurements of metal ions. Distributions in the gas thermal properties could arise from turbulence, gas cooling from the hot phase, and mixing between the cool and hot phases. Focusing on these physically motivated models, we parameterize each with a single parameter for simplicity and provide empirical and theoretical estimates for reasonable parameter values. We then construct the probability distribution functions for each of these scenarios, calculate the effective ion fractions, and fit our models to the COS-Halos absorption measurements to infer the gas densities and metallicities. We find that the models we consider (i) produce similarly good fits to the observations with or without distributions in the gas thermal properties, and (ii) result in detectable changes in the column densities only at the boundaries of reasonable parameter values. We show that He II self-shielding can have a larger effect on the ion fractions than density and temperature fluctuations. As a result, uncertainties in cloud geometry and their spatial distribution, affecting the details of radiation transfer, may obscure the effect of inhomogeneities.

The origin of Galactic cosmic rays is still a mystery, in particular the sources and acceleration mechanism for cosmic rays with energies up to or beyond a PeV. Recently LHAASO has and H.E.S.S have shown that two gamma-ray sources associated with superbubbles created by young massive stellar clusters are likely PeVatrons. This has renewed the interest in the cosmic-ray acceleration processes in superbubbles. To study the possibility and conditions under which second-order Fermi acceleration can accelerate particles beyond PeV energies in superbubbles. An analytical equation is derived for the maximum energy a cosmic-ray particle can obtain as a function of acceleration duration. The maximum energy depends critically on the diffusion coefficient D and the Alfvén velocity, $V_A$. The analytical solutions for the acceleration time scale shows that second-order Fermi acceleration can be just as efficient as diffusive shock acceleration, when comparable relevant velocities are used-i.e. the Alfvén velocity or shock velocity. The probable values for the diffusion coefficient and Alvén speed are studied for two likely PeVatron regions, HESS J1646-458 associated with Westerlund 1, and the Cygnus Cocoon, associated Cyg OB2. It is shown that within a typical stellar cluster time scale of 1-5 Myr cosmic-rays can be accelerated to $> 10^{15}$ eV, provided that $V_{\rm A} > 300$ km/s, and the diffusion coefficient is $D \sim 10^{26}$ cm$^2$/s at 100 TeV. This suggests that second-order Fermi acceleration in superbubbles should be considered as a possible source of Galactic cosmic rays up to, or beyond a PeV.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum and CO(2-1) line emission observations toward the high-mass star formation region DR21. Five new continuum sources are found. We identify eighteen outflow streamers detected in CO emission radially arising from a common origin. The velocity spread of the outflow streamers range between $-$100 to $+$70 km s$^{-1}$. The radial velocities of each outflow roughly follow linear gradients (Hubble-Lemaitre-like expansion motions). Using the CO emission of the whole ensemble of streamers we estimate a total outflow mass of 120-210M$_{\odot}$. Additionally, we derived the dynamical age (8600 yr), momentum ($\sim10^{3}$ M$_{\odot}$ km s$^{-1}$), and kinetic energy ($\sim10^{48}$ erg) of the outflow. The morphology and kinematics presented by the CO outflow streamers confirm the presence of an explosive dispersal outflow at the heart of DR21. Five dispersal explosive outflows associated with massive star-forming regions have been confirmed in our Galaxy (Orion BN/KL, G5.89-0.39, S106-IR, IRAS16076-5134 and IRAS 12326-6245). However, their frequency of occurrence in the Galaxy and the originating nature are still uncertain.

This study is aimed to contribute to a more comprehensive understanding of the molecular hydrogen distribution in the galaxy M33 by introducing novel methods for generating high angular resolution (18.2$''$, equivalent to 75 pc) column density maps of molecular hydrogen ($N_{\rm H_2}$). M33 is a local group galaxy that has been observed with Herschel in the far-infrared wavelength range from 70 to 500 $\mum$. Previous studies have presented total hydrogen column density maps ($N_{\rm H}$), using these FIR data (partly combined with mid-IR maps), employing various methods. We first performed a spectral energy distribution fit to the 160, 250, 350, and 500 $\mum$ continuum data obtain $N_{\rm H}$, using a technique similar to one previously reported in the literature. We also use a second method which involves translating only the 250 $\mum$ map into a $N_{\rm H}$ map at the same angular resolution. An $N_{\rm H_2}$ map via each method is then obtained by subtracting the HI component. Distinguishing our study from previous ones, we adopt a more versatile approach by considering a variable emissivity index, $\beta$ and dust absorption coefficient, $\kappa_0$. This choice enables us to construct a $\kappa_0$ map, thereby enhancing the depth and accuracy of our investigation of the hydrogen column density. We address the inherent biases and challenges within both methods (which give similar results) and compare them with existing maps available in the literature. Moreover, we calculate a map of the carbon monoxide CO-to-H$_2$ conversion factor (\Xco\ factor), which shows a strong dispersion around an average value of $1.8\times10^{20}\,\cmKkms$ throughout the disk. We obtain column density probability distribution functions (N-PDFs) from the $N_{\rm H}$, $N_{\rm H_2}$, and $N_{HI}$ maps and discuss their shape, consisting of several log-normal and power-law tail components.

The Sun is not quite a perfect sphere, and its oblateness, thought to be induced through its rotation, has been measured using optical observations of its radius. Its gravitational quadrupole moment can then be deduced using solar models, or through helioseismology, and it can also be determined from measurements of its gravitational effects on Mercury's orbit. The various assessments do not agree, with the most complete and precise orbital assessments being in slight excess of other determinations. This may speak to the existence of a non-luminous disk or ring, where we also note evidence for a circumsolar dust ring within Mercury's orbit from the Solar TErrestrial RElations Observatory (STEREO) mission. Historically, too, a protoplanetary disk may have been key to reconciling the Sun's metallicity with its neutrino yield. The distribution of the non-luminous mass within Mercury's orbit can modify the relative size of the optical and orbital quadrupole moments in different ways. We develop how we can use these findings to limit a dark disk, ring, or halo in the immediate vicinity of the Sun, and we note how future orbital measurements of Mercury and near-Sun asteroids can refine these constraints.

The evolution of galaxies depends on their masses and local environments; understanding when and how environmental quenching starts to operate remains a challenge. Furthermore, studies of the high-redshift regime have been limited to massive cluster members, owing to sensitivity limits or small fields of views when the sensitivity is sufficient, intrinsically biasing the picture of cluster evolution. In this work, we use stacking to investigate the average star formation history of more than 10,000 groups and clusters drawn from the Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2). Our analysis covers near ultraviolet to far infrared wavelengths, for galaxy overdensities at $0.5 \lesssim z \lesssim 2.54$. We employ SED fitting to measure the specific star formation rates (sSFR) in four annular apertures with radii between 0 and 1000 kpc. At $z \gtrsim 1.6$, the average sSFR evolves similarly to the field in both the core and the cluster outskirts. Between $\overline{z} = 1.60$ and $\overline{z} = 1.35$, the sSFR in the core drops sharply, and continues to fall relative to the field sSFR at lower redshifts. We interpret this change as evidence that the impact of environmental quenching dramatically increases at $z \sim 1.5$, with the short time span of the transition suggesting that the environmental quenching mechanism dominant at this redshift operates on a rapid timescale. We find indications that the sSFR may decrease with increasing host halo mass, but lower-scatter mass tracers than the signal-to-noise ratio (S/N) are needed to confirm this relationship.

Gamma-ray Bursts (GRBs) are one of the most energetic phenomena in the cosmos, whose study probes physics extremes beyond the reach of laboratories on Earth. Our quest to unravel the origin of these events and understand their underlying physics is far from complete. Central to this pursuit is the rapid classification of GRBs to guide follow-up observations and analysis across the electromagnetic spectrum and beyond. Here, we introduce a compelling approach for a new and robust GRB prompt classification. Leveraging self-supervised deep learning, we pioneer a previously unexplored data product to approach this task: the GRB waterfalls.

Feedback driven massive outflows play a crucial role in galaxy evolution by regulating star formation and influencing the dynamics of surrounding media. Extracting outflow properties from spectral lines is a notoriously difficult process for a number of reasons, including the possibility that a substantial fraction of the outflow is carried by dense gas in a very narrow range in velocity. This gas can hide in spectra with insufficient resolution. Empirically motivated analysis based on the Apparent Optical Depth method, commonly used in the literature, neglects the contribution of this gas, and may therefore underestimate the true gas column density. More complex semi-analytical line transfer (e.g., SALT) models, on the other hand, allow for the presence of this gas by modeling the radial density and velocity of the outflows as power laws. Here we compare the two approaches to quantify the uncertainties in the inferences of outflow properties based on 1-D "down-the-barrel" using the UV spectra of the CLASSY galaxy sample. We find that empirical modeling may significantly underestimate the column densities relative to SALT analysis, particularly in the optically thick regime. We use simulations to show that the main reason for this discrepancy is the presence of large amount of dense material at low velocities, which can be hidden by the finite spectral resolution of the data. The SALT models in turn could over-estimate the column densities if the assumed power laws of the density profiles strong are not a property of actual outflows.

Material expelled from binary neutron star (BNS) mergers can harbor r-process nucleosynthesis and power a Kilonova (KN), both intimately related to the astrophysical conditions of the ejection. In turn such conditions indirectly depend on the equation of state (EOS) describing matter inside the neutron star. Therefore, in principle the above observables can hold valuable information on nuclear matter, as the merger gravitational wave signal already does. In this work, we consider the outcome of a set of BNS merger simulations employing different finite-temperature nuclear EOSs. The latter are obtained from a Skyrme-type interaction model where nuclear properties, such as the incompressibility and the nucleon effective mass at saturation density, are systematically varied. We post-process the ejecta using a reaction network coupled with a semi-analytic KN model, to asses the sensitivity on the input EOS of the final yields and the KN light curves. Both of them are found to be non-trivially influenced by the EOS, with the overall outcome being dominated by the heterogeneous outflows from the remnant disk, hosting a variable degree of neutron-rich material. The dynamical ejecta can be more directly related to the EOS parameters considered, however, we find its role in the yields production and the KN emission too entangled with the other ejecta components, in order to infer solid correlations. This result highlights the strong degeneracy that intervenes between the merger outcome and the behaviour of the intrinsic nuclear matter, and places itself as a limit to the employment of EOS-constraining approaches of such kind.

High-redshift ($z >3$) $\gamma$-ray blazars are rare, but they are crucial for our understanding of jet evolution, $\gamma$-ray production and propagation, and the growth of supermassive black holes in the early universe. A new analysis of Fermi-LAT data reveals a significant (5$\sigma$), spectrally soft ($\Gamma \simeq$ 3.0) $\gamma$-ray source in a specific 4-month epoch, cospatial with PKS 0201+113 ($z$ = 3.64). Monitoring of PKS 0201+113 at 15 GHz by the Owens Valley Radio Observatory 40 m Telescope from 2008 to 2023 shows a prominent flare that dominates the radio light curve. The maximum of the radio flare coincides with the $\gamma$-ray flare, strongly suggesting an association ($\textrm{p-value}=0.023$) between the $\gamma$-ray and the radio sources. PKS 0201+113 is only the third $\gamma$-ray blazar to be identified with $z> 3.5$, and it is the first such object to be identified by the detection of quasi-simultaneous $\gamma$-ray and radio flares. The jet properties of this peculiar blazar have been investigated. A detailed study of a two-zone leptonic model is presented that fits the broadband spectral energy distribution. An alternative scenario is also briefly discussed.

Fast radio bursts (FRBs) are millisecond-duration flashes with unknown origins. Its formation rate is crucial for unveiling physical origins. However, the luminosity and formation rate are degenerated when directly fitting the redshift distribution of FRBs. In contrast to previous forward-fitting methods, we use the Lynden-Bell's $c^{-}$ method to derive luminosity function and formation rate of FRBs without any assumptions. Using the non-repeating FRBs from the first CHIME/FRB catalog, we find a relatively strong luminosity evolution, and luminosity function can be fitted by a broken power-law model with a break at $1.33\times10^{41}\ \mathrm{erg}\ \mathrm{s}^{-1}$. The formation rate declines rapidly as $(1+z)^{-4.9\pm0.3}$ with a local rate $1.13\times10^4\ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$. This monotonic decrease is similar to the rate of short gamma-ray bursts. After comparing it with star formation rate and stellar mass density, we conclude that the old populations including neutron stars and black holes, are closely related to the origins of FRBs. Monte Carlo simulations are used to test our results. The distributions of mock sample are consistent with the observational data.

Some electromagnetic outbursts from the nuclei of distant galaxies have been found to repeat on months-to-years timescales, and each of these sources can putatively arise from the accretion flares generated through the repeated tidal stripping of a star on a bound orbit about a supermassive black hole (SMBH), i.e., a repeating partial tidal disruption event (rpTDE). Here we test the rpTDE model through analytical estimates and hydrodynamical simulations of the interaction between a range of stars, which differ from one another in mass and age, and an SMBH. We show that higher-mass ($\gtrsim 1 M_{\odot}$), evolved stars can survive many ($\gtrsim 10-100$) encounters with an SMBH while simultaneously losing $few \times 0.01 M_{\odot}$, resulting in accretion flares that are approximately evenly spaced in time with nearly the same amplitude, quantitatively reproducing ASASSN-14ko. We also show that the energy imparted to the star via tides can lead to a change in its orbital period that is comparable to the observed decay in the recurrence time of ASASSN-14ko's flares, $\dot{P}\simeq-0.0016$. Contrarily, lower-mass and less-evolved stars lose progressively more mass and produce brighter accretion flares on subsequent encounters for the same pericenter distances, leading to the rapid destruction of the star and cessation of flares. Such systems cannot reproduce ASASSN-14ko-like transients, but are promising candidates for recreating events such as AT2020vdq, which displayed a second and much brighter outburst compared to the first. Our results imply that the lightcurves of repeating transients are tightly coupled with stellar type.

Accretion discs are fundamental to much of astronomy. They can occur around stars both young and old, around compact objects they provide a window into the extremes of physics, and around supermassive black holes in galaxy centres they generate spectacular luminosities that can outshine the entire galaxy. However, our understanding of the inner workings of accretion discs remains far from complete. Here we revisit a conundrum in the observations of some of the simplest accreting systems; the Cataclysmic Variables (CVs). The high-accretion-rate states of (non-magnetic) CVs can be divided into the short-lived outbursts ($\sim$ a week) typical of dwarf novae (DNe) and the long-lived (and sometimes perpetual) high states of nova-like (NL) CVs. Since both sorts of high-state occur in approximately steady-state accretion discs with similar properties and accretors, we would expect them to display similar spectral energy distributions. However, previous analyses based on UV spectra from the {\it International Ultraviolet Explorer} have shown that their spectral energy distributions are different. We perform a re-analysis of the data using up to date calibrations and distance (and thus dereddening) estimates to test whether this difference persists and whether it is statistically significant over the sample. We find that it does persist and it is statistically significant. We propose routes to investigating this discrepancy further and discuss the implications this has for other accreting systems, such as X-ray binaries, Active Galactic Nuclei and protoplanetary discs.

The pion-bump structure in the gamma-ray spectrum is a direct proof for the hadronic origin of the gamma rays, and thus the decisive evidence for the acceleration of hadronic cosmic rays in astrophysical objects. However, the identification of such a spectral feature is limited by the resolution and energy coverage of current gamma-ray instruments. Furthermore, there are unavoidable bremsstrahlung emissions from secondary and primary electrons, which may dominate the gamma-ray emission below the pion-bump. Thus, the study of this gamma-ray emission component can provide unique information on the acceleration and confinement of high-energy particles. In this paper, we studied the predicted gamma-ray spectrum assuming both hadronic or leptonic origin in mid-aged supernova remnants W44, we discuss the detection potential of future MeV missions on these emissions and possible implications.

We investigate the geometry of the magnetic field towards the Radcliffe Wave, a coherent 3-kpc-long part of the nearby Local Arm recently discovered via three-dimensional dust mapping. We use archival stellar polarization in the optical and new measurements in the near-infrared to trace the magnetic field as projected on the plane of the sky. Our new observations cover the portion of the structure that is closest to the Sun, between Galactic longitudes of 122$^\circ$ and 188$^\circ$. The polarization angles of stars immediately background to the Radcliffe Wave appear to be aligned with the structure as projected on the plane of the sky. The observed magnetic field configuration is inclined with respect to the Galactic disk at an angle of 18$^\circ$. This departure from a geometry parallel to the plane of the Galaxy is contrary to previous constraints from more distant stars and polarized dust emission. We confirm that the polarization angle of stars at larger distances shows a mean orientation parallel to the Galactic disk. We discuss implications of the observed morphology of the magnetic field for models of the large-scale Galactic magnetic field, as well as formation scenarios for the Radcliffe Wave itself.

Redshift estimation and the classification of gamma-ray AGNs represent crucial challenges in the field of gamma-ray astronomy. Recent efforts have been made to tackle these problems using traditional machine learning methods. However, the simplicity of existing algorithms, combined with their basic implementations, underscores an opportunity and a need for further advancement in this area. Our approach begins by implementing a Bayesian model for redshift estimation, which can account for uncertainty while providing predictions with the desired confidence level. Subsequently, we address the classification problem by leveraging intelligent initialization techniques and employing soft voting. Additionally, we explore several potential self-supervised algorithms in their conventional form. Lastly, in addition to generating predictions for data with missing outputs, we ensure that the theoretical assertions put forth by both algorithms mutually reinforce each other.

We present a novel multimodal neural network for classifying astronomical sources in multiband ground-based observations, from optical to near infrared, to separate sources in stars, galaxies and quasars. Our approach combines a convolutional neural network branch for learning morphological features from $r$-band images with an artificial neural network branch for extracting spectral energy distribution (SED) information. Specifically, we have used 9-band optical ($ugri$) and NIR ($ZYHJK_s$) data from the Kilo-Degree Survey (KiDS) Data Release 5. The two branches of the network are concatenated and feed into fully-connected layers for final classification. We train the network on a spectroscopically confirmed sample from the Sloan Digital Sky Survey cross-matched with KiDS. The trained model achieves 98.76\% overall accuracy on an independent testing dataset, with F1 scores exceeding 95\% for each class. Raising the output probability threshold, we obtain higher purity at the cost of a lower completeness. We have also validated the network using external catalogs cross-matched with KiDS, correctly classifying 99.74\% of a pure star sample selected from Gaia parallaxes and proper motions, and 99.74\% of an external galaxy sample from the Galaxy and Mass Assembly survey, adjusted for low-redshift contamination. We apply the trained network to 27,334,751 KiDS DR5 sources with $r \leqslant 23$ mag to generate a new classification catalog. This multimodal neural network successfully leverages both morphological and SED information to enable efficient and robust classification of stars, quasars, and galaxies in large photometric surveys.

Recent discoveries of multiple long-period pulsars (periods ${\sim}10\,$s or larger) are starting to challenge the conventional notion that coherent radio emission cannot be produced by objects that are below the many theorised death lines. Many of the past pulsar surveys and software have been prone to selection effects that restricted their sensitivities towards long-period and sporadically-emitting objects. Pulsar surveys using new-generation low-frequency facilities are starting to employ longer dwell times, which makes them significantly more sensitive in detecting long-period or nulling pulsars. There have also been software advancements to aid more sensitive searches towards long-period objects. Furthermore, recent discoveries suggest that nulling may be a key aspect of the long-period pulsar population. We simulate both long-period and nulling pulsar signals, using the Southern-sky MWA Rapid Two-meter (SMART) survey data as reference, and explore the detection efficacy of popular search methods such as the fast Fourier transform (FFT), fast-folding algorithm (FFA) and single pulse search (SPS). For FFT-based search and SPS, we make use of the PRESTO implementation, and for FFA we use RIPTIDE. We find RIPTIDE's FFA to be more sensitive; however, it is also the slowest algorithm. PRESTO's FFT, although faster than others, also shows some unexpected inaccuracies in detection properties. SPS is highly sensitive to long-period and nulling signals, but only for pulses with high intrinsic signal-to-noise ratios. We use these findings to inform current and future pulsar surveys that aim to uncover a large population of long-period or nulling objects and comment on how to make optimal use of these methods in unison.

In spite of significant progress in the research of fast radio bursts (FRBs) in recent decade, their origin is still under extensive debate. Investigation on the population of FRBs can provide new insight into this interesting problem. In this paper, based on the first CHIME/FRB catalog, we construct a Bayesian framework to analyze the FRB population, with the selection effect of the CHIME telescope being properly taken into account. The energy function is modeled as the power-law with an exponential cutoff. Four redshift distribution models are considered, i.e., the star formation history (SFH) model, and three time-delayed models (Gaussian delay, log-normal delay, and power-law delay). The free parameters are simultaneously constrained using Bayesian inference method, and the Bayesian information criterion (BIC) is used in model comparison. According to BIC, the log-normal delay model fits the data best. The power-law delay model and Gaussian delay model can also give reasonable fits, although they are not as good as the log-normal delay model. However, the SFH model is strongly disfavored compared with the three time-delayed models. The energy function is tightly constrained and is almost independent of the redshift models, with the best-fitting power-law index $\alpha\approx 1.8$, and cut-off energy $\log(E_c/{\rm erg})\approx 42$. The FRB population shows on average $3\sim 5$ billion years time delay with respect to the SFH. Therefore, the hypothesis that the FRB population traces the SFH is conclusively ruled out.

Simulations are the best approximation to experimental laboratories in astrophysics and cosmology. However, the complexity, richness, and large size of their outputs severely limit the interpretability of their predictions. We describe a new, unbiased, and machine learning based approach to obtaining useful scientific insights from a broad range of simulations. The method can be used on today's largest simulations and will be essential to solve the extreme data exploration and analysis challenges posed by the Exascale era. Furthermore, this concept is so flexible, that it will also enable explorative access to observed data. Our concept is based on applying nonlinear dimensionality reduction to learn compact representations of the data in a low-dimensional space. The simulation data is projected onto this space for interactive inspection, visual interpretation, sample selection, and local analysis. We present a prototype using a rotational invariant hyperspherical variational convolutional autoencoder, utilizing a power distribution in the latent space, and trained on galaxies from IllustrisTNG simulation. Thereby, we obtain a natural Hubble tuning fork like similarity space that can be visualized interactively on the surface of a sphere by exploiting the power of HiPS tilings in Aladin Lite.

Assessing the origin of Pluto and Triton has profound implications for the bigger picture of Solar System formation and evolution. In such a context, this chapter reviews our current knowledge of the formation conditions of Pluto and Triton's constitutive building blocks in the protosolar nebula, which can be derived from their known or estimated volatile contents. Assuming that the ultravolatiles carbon monoxide and dinitrogen detected in Pluto and Triton are primordial, the presence of these molecules suggest that the two bodies accreted material originating from the vicinity of the carbon monoxide and dinitrogen icelines. Dinitrogen--rich and water--poor comets such as comet C/2016 R2 (PanSTARRS) obviously present a compositional link with Pluto and Triton, indicating that their building blocks formed in nearby regions of the protosolar nebula, despite of the variation of the water abundance among those bodies. Also, the assumption of Triton's growth in Neptune's circumplanetary disk requires that its building blocks formed at earlier epochs in the protosolar nebula, to remain consistent with its estimated composition.

As a common gravitation virialized object in the standard $\Lambda$CDM cosmology, dark matter halo connects from the large-scale structure all the way down to galaxy and star formation. However, as the nature of dark matter particles is still unclear, the smallest halo that can be formed in the universe is still unknown. Based on some simple assumptions, this paper uses the \textsc{hmf} package to investigate different halo functions used to quantify its number and mass distributions -- the halo mass function and the integrated/differential mass function (IMF/DMF) respectively. The halo mass in this study extends from the galaxy cluster to the dark matter particle mass at the GeV scale. Surprisingly, different fitting functions for the HMF are in remarkable agreement, a scatter within 2 orders of magnitude, down to dark matter particle mass, of which the halo mass spans about 80 orders of magnitude and the HMF covers over 100 orders of magnitude. The DMF reveals an interesting and consistent peak at $\sim 10^{13} \hMsun$, which implies galaxy groups have the highest contribution to the total matter mass. Furthermore, the effects of cosmology parameters on these halo functions are also examined with the most massive halos, or these halo functions at the most massive halo mass end, more sensitive to them. Different behaviours of these halo functions due to the changes in cosmology parameters can be used to break the degeneracy between them.

Unsupervised learning algorithms like self-organizing Kohonen maps are a promising approach to gain an overview among massive datasets. With UltraPINK, researchers can train, inspect, and explore self-organizing maps, whereby the toolbox of interaction possibilities grows continually. Key feature of UltraPINK is the consideration of versality in astronomical data. By keeping the operations as abstract as possible and using design patterns meant for abstract usage, we ensure that data is compatible with UltraPINK, regardless of its type, formatting, or origin. Future work on the application will keep extending the catalogue of exploration tools and the interfaces towards other established applications to process astronomical data. Ultimatively, we aim towards a solid infrastructure for data analysis in astronomy.

We report data analysis results about the outburst evolution and spectral properties during the hard state of the recently discovered X-ray transient Swift J1727.8-163 as observed by \emph{Insight}-HXMT and NuSTAR. We find that the broadband X-ray spectrum of Swift J1727.8-163 is more complex than the most typical spectral patterns of black hole X-ray binary systems, with not only a comparatively weaker reflection component but also an additional spectral continuum component, manifesting itself as a hard X-ray tail beyond the thermal Comptonization description detectable below 100 keV. This additional component can be phenomenologically well fitted by adding an extra power-law model with high energy exponential cutoff in the 2-120 keV energy band. We made an attempt to explain the broadband X-ray spectral continuum with a thermal/non-thermal hybrid plasma corona scenario , and find an ultra high compactness parameter ($l_{\rm s}\sim2000$) and a steep non-thermal electron distribution ($\Gamma_{\rm inj}>4$), suggesting the source was accreting with high Eddington rates and that the electron acceleration mechanism is not very efficient. We also present a detailed multi-epoch analysis of spectral properties using \emph{Insight}-HXMT data to investigate the evolution of the key physical properties regarding the disk and corona during the hard states. No significant variation is found with the inner disk radius and the coronal temperature during this time period, and the weak reflection and hard X-ray tail features are persistent. We discuss the physical implications of our spectral analysis results in the context of disk-corona relation, particle acceleration, and jet contribution, during the rise of a black hole X-ray binary in outburst.

Progenitor stars of long gamma-ray bursts (GRBs) could be surrounded by a significant and complex nebula structure lying at a parsec scale distance. After the initial release of energy from the GRB jet, the jet will interact with this nebula environment. We show here that for a large, plausible parameter space region, the interaction between the jet blastwave and the wind termination (reverse) shock is expected to be weak, and may be associated with a precursor emission. As the jet blast wave encounters the contact discontinuity separating the shocked wind and the shocked interstellar medium, we find that a bright flash of synchrotron emission from the newly-formed reverse shock is produced. This flash is expected to be observed at around ~100 s after the initial explosion and precursor. Such a delayed emission thus constitutes a circumburst medium (CBM) phase in a GRB, having a physically distinct origin from the preceding prompt phase and the succeeding afterglow phase. The CBM phase emission may thus provide a natural explanation to bursts observed to have a precursor followed by an intense, synchrotron-dominated main episode that is found in a substantial minority, ~10% of GRBs. A correct identification of the emission phase is thus required to infer the properties of the flow and of the immediate environment around GRB progenitors.

Evolutionary pathways of binary systems are vastly different from single stellar evolution, and thus, there is a need to quantify their frequency and diversity. Open clusters are the best test-bed to unveil the secrets of binary populations due to their coeval nature. And the availability of multi-wavelength data in recent years has been critical in characterising the binary population. NGC 752 is a solar metallicity, intermediate-age open cluster located at 460 pc. In this work, we aim to identify the optically subluminous white dwarfs in NGC 752 and identify the illusive blue lurkers by association. We used multiwavelength photometry from Astrosat/UVIT, swift/UVOT, Gaia DR3 and other archival surveys to analyse the colour-magnitude diagrams and spectral energy distributions of 37 cluster members. We detected eight white dwarfs as companions to cluster members. Four of the systems are main sequence stars with extremely low mass white dwarfs as their companions. Two are these main sequence stars are also fast rotators. The presence of low mass white dwarfs and high rotation signals a past mass transfer, and we classified the four main sequence stars as blue lurkers. The binary fraction in NGC 752 was estimated to be 50--70%, and it shows that the contribution of optically undetected stars is crucial in quantifying the present-day binary fraction.

This paper presents an investigation of spectral properties of 10 millisecond pulsars (MSPs) discovered by the uGMRT, observed from 2017-2023 using band 3 (300-500 MHz) and 4 (550-750 MHz) of uGMRT. For these MSPs, we have reported a range of spectral indices from ~0 to -4.8, while averaging the full observing band and all the observing epochs. For every MSP, we calculated the mean flux densities across 7-8 sub-bands each with approximately 25 MHz bandwidth spanning band 3 and band 4. We computed their modulation indices as well as average and maximum-to-median flux densities within each subband. Using a temporal variation of flux density we calculated the refractive scintillation time scales and estimated structure function with time lag for 8 MSPs in the sample. We note a significant temporal evolution of the in-band spectra, classified into three categories based on the nature of the best-fit power-law spectra, having single positive spectral indices, multiple broken power law, and single negative spectral indices. Additionally, indications of low-frequency turnover and a temporal variation of the turnover frequency (to the extent that turnover was observed for some of the epochs while not seen for the rest) were noted for all the MSPs. To the best of our knowledge, this is the first systematic investigation probing temporal changes in the MSP spectra as well as in turnover frequency. Future exploration with dense monitoring combined with modeling of spectra can provide vital insight into the intrinsic emission properties of the MSPs and ISM properties.

Schaefer (2024) has recently published observations of binary period derivatives $\dot P$ for 52 cataclysmic variables, and concluded that these strongly conflict with all proposed evolutionary pictures for these systems. We point out once again that using measurements of $\dot P$ is likely in practice to produce misleading evolutionary constraints in almost every case. The one identified exception to this is probably the recently-born X-ray binary SN 2022jli, because of its extremely high mass transfer rate.

Gas kinematics affect the radiative transfer and escape of hydrogen Lyman-$\alpha$ (Ly$\alpha$) emission from galaxies. We investigate this interplay empirically by relating the ionised gas kinematics of 42 galaxies in the extended Ly$\alpha$ Reference Sample with their Ly$\alpha$ escape fractions, $f_\rm{esc}$, Ly$\alpha$ equivalent widths, $\rm{EW}_\rm{Ly\alpha}$, and Ly$\alpha$ luminosities, $L_\rm{Ly\alpha}$. To this aim we use PMAS integral-field spectroscopic observations of the Balmer-$\alpha$ line. We calculate shearing velocities, $v_\rm{shear}$, and intrinsic velocity dispersions, $\sigma_0^\rm{obs}$ (empirically corrected for beam-smearing effects), as global kinematical measures for each galaxy. The sample is characterised by highly turbulent motions and more than half of the sample shows dispersion dominated kinematics. We uncover clear trends between Ly$\alpha$ observables and global kinematical statistics. We discuss statistically the importance of $v_\rm{shear}$, $\sigma_0^\rm{obs}$, and $v_\rm{shear}/\sigma_0^\rm{obs}$ for regulating the Ly$\alpha$ observables in comparison to other galaxy parameters. It emerges that $\sigma_0^\rm{obs}$ is the dominating parameter for $\rm{EW}_\rm{Ly\alpha}$ and that is as important as nebular extinction, gas covering fraction, and ionising photon production efficiency in regulating $f_\rm{esc}$. A simple scenario where the starburst age is simultaneously regulating turbulence, $\rm{EW}_\rm{Ly\alpha}$, and $f_\rm{esc}$ does not find support by our observations. However, we show that the small scale distribution of dust appears to be influenced by turbulence in some galaxies. In support of our observational result we discuss how turbulence is theoretically expected to play a significant role in modulating $f_\rm{esc}$. (abridged)

Bayesian evidence is a standard method used for comparing the ability of different models to fit available data and is used extensively in cosmology. However, since the evidence calculation involves performing an integral of the likelihood function over the entire space of model parameters this can be prohibitively expensive in terms of both CPU and time consumption. For example, in the simplest $\Lambda$CDM model and using CMB data from the Planck satellite, the dimensionality of the model space is over 30 (typically 6 cosmological parameters and 28 nuisance parameters). Even the simplest possible model requires $\mathcal{O}(10^6)$ calls to an Einstein--Boltzmann solver such as CLASS or CAMB and takes several days. Here we present calculations of Bayesian evidence using the CONNECT framework to calculate cosmological observables. We demonstrate that we can achieve results comparable to those obtained using Einstein--Boltzmann solvers, but at a minute fraction of the computational cost. As a test case, we then go on to compute Bayesian evidence ratios for a selection of slow-roll inflationary models. In the setup presented here, the total computation time is completely dominated by the likelihood function calculation which now becomes the main bottleneck for increasing computation speed.

We perform the numerical simulations of axisymmetric, relativistic, optically thin jets under the influence of the radiation field of an accretion disk. We show that starting from a very low injection velocity at the base, jets can be accelerated to relativistic terminal speeds when traveling through the radiation field. The jet gains momentum through the interaction with the radiation field. We use a relativistic equation of state for multi-species plasma, which self-consistently calculates the adiabatic index for the jet material. All the jet solutions obtained are transonic in nature. In addition to the acceleration of the jet to relativistic speeds, our results show that the radiation field also acts as a collimating agent. The jets remain well collimated under the effect of radiation pressure. We also show that if the jet starts with a rotational velocity, the radiation field will reduce the angular momentum of the jet beam.

Excess of gamma rays with a spherical morphology around the Galactic center (GC) observed in the Fermi large area telescope (LAT) data is one of the most intriguing features in the gamma-ray sky. The excess has been interpreted by annihilating dark matter as well as emission from a population of unresolved millisecond pulsars (MSPs). We use a multi-class classification of Fermi-LAT sources with machine learning to study the distribution of MSP-like sources among unassociated Fermi-LAT sources near the GC. We find that the source count distribution of MSP-like sources is comparable with the MSP explanation of the GC excess.

Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($\Delta \mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($\Delta \mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions.

With its extreme density of stars and stellar remnants, dense young massive clusters, high specific star formation rate, intense radiation field, high magnetic field strength, and properties of the interstellar medium that resemble those in high redshift galaxies and starbursts, the Galactic Centre is the most extreme environment that we can observe in detail. It is also the only nucleus of a galaxy that we can observe with a resolution of just a few milli parsecs. This makes it a crucial target to understand the physics of galactic nuclei and star formation, as well as the connection between them. It enables studies of a large number of otherwise rare objects, such as extremely massive stars and stellar remnants, at a well-defined distance, thus facilitating the interpretation of their properties. The Galactic Centre has been and is being studied intensively with the most advanced facilities. In this White Paper, we advocate for a large-area, multi-wavelength survey with the Square Kilometre Array of an area of about 1.25x0.3 deg**2 (180x40 pc**2), centered on the massive black hole Sagittarius A* and for repeated deep observations of the nuclear star cluster over a decade, which will allow the community to address multiple science problems with a single data set.

Accurate modelling of the one-to-two halo transition has long been difficult to achieve. We demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to: 1) particles \it orbiting \rm the halo; and 2) particles \it infalling \rm into the halo for the first time. Motivated by this decomposition, García [R. García et. al., MNRAS 521, 2464 (2023)] advocated for defining haloes as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo--mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterised by a single physical scale $r_{\rm h}$, which we refer to as the \it halo radius. \rm The two-halo term is described using a simple yet accurate empirical model based on the Zel'dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo--mass correlation function is accurate at the $\approx 2\%$ level for $r \in [0.1,50]\ h^{-1}\ Mpc$. The Fourier transform of our best fit model describes the halo--mass power spectrum with comparable accuracy for $k\in [0.06, 6.0]\ h\ Mpc^{-1}$.

Stellar evolution theory restricted to single stars predicts a minimum mass for core-collapse supernovae (CCSNe) of around eight solar masses; this minimum mass corresponds to a maximum age of around 45 million years for the progenitor and the coeval population of stars. Binary evolution complicates this prediction. For example, an older stellar population around 100 million years could contain stellar mergers that reach the minimum mass for core collapse. Despite this clear prediction by binary evolution, there are few, if any CCSNe associated with a distinctly older stellar population...until now. The stellar population within 150 pc of the Vela Pulsar is inconsistent with single-star evolution only; instead, the most likely solution is that the stellar population is $\ge$80 Myr old, and the brightest stars are mass gainers and/or mergers, the result of binary evolution. The evidence is as follows. Even though the main sequence is clearly dominated by a $\ge$80-Myr-old population, a large fraction of the corresponding red giants is missing. The best-fitting single-star model expects 51.5 red giants, yet there are only 22; the Poisson probability of this is $1.7 \times 10^{-6}$. In addition, there is an overabundance of bright, young-looking stars (25-30 Myrs old), yet there is not a corresponding young main sequence (MS). Upon closer inspection, the vast majority of the young-looking stars show either past or current signs of binary evolution. These new results are possible due to exquisite Gaia parallaxes and a new age-dating software called {\it Stellar Ages}.

Context: Young massive stellar clusters (YMSCs) have emerged as potential $\gamma$-ray sources, after the recent association of a dozen YMSCs with extended $\gamma$-ray emission. The large size of the detected halos, comparable to that of the wind-blown bubble expected around YMSCs, makes the $\gamma$-ray detection of individual YMSCs rather challenging. As a result, the emission from most of the Galactic YMSCs could be unresolved, thus contributing to the diffuse $\gamma$-ray radiation observed along the Galactic Plane. Aims: In this study, we estimate the possible contribution to the Galactic diffuse $\gamma$-ray emission from a synthetic population of YMSCs, and we compare it with observations obtained with different experiments, from 1 GeV to hundreds of TeV, in two regions of the Galactic Plane. Methods: As the population of galactic YMSCs is only known locally, we evaluate the contribution of $\gamma$-ray emission relying on the simulation of synthetic populations of YMSCs based on the observed properties of local clusters. We compute the $\gamma$-ray emission from each cluster assuming that the radiation is purely hadronic in nature and produced by cosmic rays accelerated at the cluster's collective wind termination shock. Results: We find that the $\gamma$-ray emission from unresolved YMSCs can significantly contribute to the observed Galactic diffuse flux, especially in the inner part of the Galaxy. The result is independent of the assumed particle transport, but an important role is played by Wolf-Rayet stars. The predicted $\gamma$-ray flux should be considered as a lower limit, given that our calculation neglects the contribution of supernovae exploding in YMSCs.

Common envelopes are thought to be the main method for producing tight binaries in the universe as the orbital period shrinks by several orders of magnitude during this phase. Despite their importance for various evolutionary channels, direct detections are rare, and thus observational constraints on common envelope physics are often inferred from post-CE populations. Population constraints suggest that the CE phase must be highly inefficient at using orbital energy to drive envelope ejection for low-mass systems and highly efficient for high-mass systems. Such a dichotomy has been explained by an interplay between convection, radiation and orbital decay. If convective transport to the surface occurs faster than the orbit decays, the CE self-regulates and radiatively cools. Once the orbit shrinks such that convective transport is slow compared to orbital decay, a burst occurs as the release of orbital energy can be far in excess of that required to unbind the envelope. With the anticipation of first light for the Rubin Observatory, we calculate light curve models for convective common envelopes and provide the time evolution of apparent magnitudes for the Rubin filters. Convection imparts a distinct signature in the light curves and lengthens the timescales during which they are observable. Given Rubin limiting magnitudes, convective CEs should be detectable out to distances of ~8 Mpc at a rate of ~0.3 per day and provide an intriguing observational test of common envelope physics.

KIC 4150611 is a high-order multiple composed of a triple system composed of the F1V primary (Aa), which is eclipsed on a 94.2d period by a tight 1.52d binary composed of two dim K/M dwarfs (Ab1, Ab2), which also eclipse each other; an 8.65d eccentric, eclipsing binary composed of two G stars (Ba, Bb); and another faint eclipsing binary composed of two stars of unknown spectral type (Ca and Cb). In addition to its many eclipses, the system is an SB3 spectroscopic multiple (Aa, Ba, and Bb) and the primary (Aa) is a hybrid pulsator. We employ a novel photometric analysis of the complicated eclipse geometry of Aa to obtain orbital and stellar properties of the triple. We acquired 51 TRES spectra at the Fred L. Whipple Observatory, calculating radial velocities and orbital elements of Aa (SB1) and the B binary (SB2). These spectra and radial velocities are used to perform spectral disentangling for Aa, Ba, and Bb. Spectral modelling is applied to the disentangled spectrum of Aa to obtain atmospheric properties. We obtain precise stellar properties of the triple, including the mass ratios (MAa/(MAb1 + MAb2) = 3.61 +/- 0.01, MAb1/MAb2 = 1.113 +/- 0.001), separation ratio (aAab/aAb1Ab2 = 21.81 +/- 0.01), orbital periods (PAab = 94.29486 +/- 0.00008d, PAb1Ab2 = 1.522248 +/- 0.000001d), and stellar radii (RAa = 1.64 +/- 0.06 Rsun, RAb1 = 0.42 +/- 0.01 Rsun, RAb2 = 0.38 +/- 0.01 Rsun). Radial velocity fitting and spectral disentangling arrive at orbital elements for Aa, Ba, and Bb in excellent agreement with each other and with previous results in the literature. Spectral modelling on the disentangled spectrum of Aa provides constraints on the effective temperature (Teff = 7280 +/- 70 K), surface gravity (log(g) = 4.14 +/- 0.18 dex), micro-turbulent velocity (vmicro = 3.61 +/- 0.19 km s-1), rotation velocity (v sin i = 127 +/- 4 km s-1), and metallicity ([M/H] = -0.23 +/- 0.06).

Motivated by the recent detection of ultra-long period radio transients, we investigate new models of coherent radio emission via low-altitude electron-positron pair production in neutron stars beyond rotationally-powered curvature radiation deathlines. We find that plastic motion (akin to 'continental drift') and qualitatively similar thermoelectric action by temperature gradients in the crusts of slowly rotating, highly magnetized neutron stars could impart mild local magnetospheric twists. Regardless of which mechanism drives twists, we find that particle acceleration initiates pair cascades across charge-starved gaps above a mild critical twist. Cascades are initiated via resonant inverse-Compton scattered photons or curvature radiation, and may produce broadband coherent radio emission. We compute the pair luminosity (maximum allowed radio luminosity) for these two channels, and derive deathlines and 'active zones' in $P-\dot{P}$ space from a variety of considerations. We find these twist-initiated pair cascades only occur for magnetar-like field strengths $B \gtrsim 10^{14}$ G and long periods: $P_{\rm RICS} \gtrsim 120 \; (T/10^{6.5} {\rm K})^{-5} \, {\rm sec}$ and $P_{\rm curv} \gtrsim 150 \; ({\rm v_{\rm pl}}/10^{3} {\, \rm cm \, yr^{-1}})^{-7/6} \, {\rm sec}$. Using a simplified geometric model, we find that plastic motion or thermoelectrically driven twists might naturally reproduce the observed luminosities, timescales, and timing signatures. We further derive 'active zones' in which rotationally-powered pair creation occurs via resonantly scattered photons, beyond standard curvature deathlines for pulsars. All cascades are generically accompanied by simultaneous (non-)thermal X-ray/UV counterparts which might be detectable with current instrumentation.

The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary disk. The multi-epoch spectroscopy yields photospheric lines to classify the star, and emission lines to study variability in the hot inner disk and to determine the mass-accretion rate. The SPHERE and ALMA observations are used to characterize the dust distribution of the small and large grains, respectively. PDS 111 is a weak-line T Tauri star with spectral type G2, exhibits strong H$\alpha$ variability and with a low mass-accretion rate of $1-5\times10^{-10}$\,M$_{\odot}$\,yr$^{-1}$. We measured an age of the system of 15.9$^{+1.7}_{-3.7}$ Myr using pre-main sequence tracks. The SPHERE observations show a strongly flaring disk with an asymmetric substructure. The ALMA observations reveal a 30 au cavity in the dust continuum emission with a low contrast asymmetry in the South-West of the disk and a dust disk mass of 45.8\,$M_\oplus$. The $^{12}$CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than suggested in literature, PDS 111 still seems relatively old; this provides insight into disk properties at an advanced stage of pre-main sequence evolution. The characteristics of this disk are very similar to its younger counterparts: strongly flaring, an average disk mass, a typical radial extent of the disk gas and dust, and the presence of common substructures. This suggests that disk evolution has not significantly changed the disk properties. These results show similarities with the "Peter Pan disks" around M-dwarfs.

The presence of young massive stars in the Galactic Centre (GC) raises questions about star formation near the black hole Sagittarius A* (Sgr A*). Additionally, the initial mass function (IMF) in this region appears different from the standard Salpeter/Kroupa law. Extreme extinction and crowding limit our understanding of the stellar population, with spectroscopic data available only for selected bright sources. We aim to improve knowledge about the distribution and IMF of young, massive stars near Sgr A*. Using intermediate band (IB) photometry, we identify candidates for massive young stars through Bayesian inference, a neural network, and a gradient-boosted trees algorithm. We obtained spectral energy distributions for 6590 stars, 1181 of which have been previously classified spectroscopically. We identify 351 stars classified as early types by all three methods, including 155 newly identified candidates. The radial density profiles for late and early-type stars fit broken power laws, with a break radius of 9.2 +- 0.6'' for early-type stars. Late-type stars show a core-like distribution around Sgr A*, while early-type stars' density increases steeply towards the black hole. We infer a top-heavy IMF of young stars near Sgr A* (R < 9''), with a power-law of 1.6 +- 0.1. At greater distances, a standard Salpeter/Kroupa IMF fits the data. IB photometry also constrains the metallicities of late-type stars, estimating metallicities for over 600 stars. The IMF variation with radial distance suggests different star formation mechanisms, with a top-heavy IMF near Sgr A* consistent with disc star formation.

The observed radiation from hot gas accreting onto a black hole depends on both the details of the flow and the spacetime geometry. The lensing behavior of a black hole produces a distinctive pattern of autocorrelations within its photon ring that encodes its mass, spin, and inclination. In particular, the time autocorrelation of the light curve is expected to display a series of peaks produced by light echoes of the source, with each peak delayed by the characteristic time lapse $\tau$ between light echoes. However, such peaks are absent from the light curves of observed black holes. Here, we develop an analytical model for such light curves that demonstrates how, even though light echoes always exist in the signal, they do not produce autocorrelation peaks if the characteristic correlation timescale $\lambda_0$ of the source is greater than $\tau$. We validate our model against simulated light curves of a stochastic accretion model ray traced with a general-relativistic code, and then fit the model to an observed light curve for Sgr A*. We infer that $\lambda_0>\tau$, providing an explanation for the absence of light echoes in the time autocorrelations of Sgr A* light curves. Our results highlight the importance for black hole parameter inference of spatially resolving the photon ring via future space-based interferometry.

We present a method to obtain rapid mass loss rates in binary systems, specifically at the onset of MT episodes. The method unifies atmospheric (underflow) and $L_1$ stream (overflow) mass rates in a single continuous procedure. The method uses averaged 3D properties of the binaries, such as effective binary potential and effective binary acceleration, to both evolve the donor and obtain properties of the matter at the $L_1$ plane. In the case of underflow, we obtain atmospheric stratification. Our method can be used for binaries with an extensive range of mass ratios, $0.01 \le q \le 100$, and can also be applied to hot donors. The considered examples show that the MT rates obtained with this revised formalism always differ from the optically thin and optically thick MT rates widely used during the computations of binary evolution.

EBLM J0608-59 / TOI-1338 / BEBOP-1 is a 12th-magnitude, F9V star in an eclipsing binary with a much fainter M-dwarf companion on a wide, eccentric orbit (P=14.6 d). The binary is orbited by two circumbinary planets: one transiting on a 95-day orbit and one non-transiting on a 215-day orbit. We have used high-precision photometry from the TESS mission combined with direct mass measurements for the two stars published recently to measure the following model-independent radii: $R_1 = 1.32 \pm 0.02 R_{\odot}$, $R_2 = 0.309 \pm 0.004 R_{\odot}$. Using $R_1$ and the parallax from Gaia EDR3 we find that this star's angular diameter is $\theta = 0.0309 \pm 0.0005$ mas. The apparent bolometric flux of the primary star corrected for both extinction and the contribution from the M-dwarf ($<0.4$%) is ${\mathcal F}_{\oplus,0} = (0.417\pm 0.005)\times10^{-9} {\rm \,erg\,cm}^{-2} {\rm \,s}^{-1}$. Hence, this F9V star has an effective temperature $T_{\rm eff,1} = 6031{\rm\,K} \pm 46{\rm \,K\,(rnd.)} \pm 10 {\rm \,K\,(sys.)}$. EBLM J0608-59 is an ideal benchmark star that can be added to the sample of such systems we are establishing for "end-to-end" tests of the stellar parameters measured by large-scale spectroscopic surveys.

We study 51 jellyfish galaxy candidates in the Fornax, Antlia, and Hydra clusters. These candidates are identified using the JClass scheme based on the visual classification of wide-field, twelve-band optical images obtained from the Southern Photometric Local Universe Survey. A comprehensive astrophysical analysis of the jellyfish (JClass > 0), non-jellyfish (JClass = 0), and independently organized control samples is undertaken. We develop a semi-automated pipeline using self-supervised learning and similarity search to detect jellyfish galaxies. The proposed framework is designed to assist visual classifiers by providing more reliable JClasses for galaxies. We find that jellyfish candidates exhibit a lower Gini coefficient, higher entropy, and a lower 2D Sérsic index as the jellyfish features in these galaxies become more pronounced. Jellyfish candidates show elevated star formation rates (including contributions from the main body and tails) by $\sim$1.75 dex, suggesting a significant increase in the SFR caused by the ram-pressure stripping phenomenon. Galaxies in the Antlia and Fornax clusters preferentially fall towards the cluster's centre, whereas only a mild preference is observed for Hydra galaxies. Our self-supervised pipeline, applied in visually challenging cases, offers two main advantages: it reduces human visual biases and scales effectively for large datasets. This versatile framework promises substantial enhancements in morphology studies for future galaxy image surveys.

In this study, we apply the Velocity Gradient Technique (VGT) to the merging Centaurus galaxy. We compare gradient maps derived from the PHANGS-ALMA survey using CO emission lines with magnetic field tracings from dust polarization data obtained via the HAWC+ instrument. Our analysis reveals a strong correspondence between the directions indicated by these two tracers across most of the galactic image. Specifically, we identify jet regions as areas of anti-alignment, consistent with previous reports that gradients tend to rotate 90 degrees in outflow regions. Statistically, we find that the alignment of magnetic fields, as revealed by polarization, is most accurate in regions with the highest signal-to-noise ratios. Our findings underscore the utility of velocity gradients as a valuable complementary tool for probing magnetic fields and dynamical processes in merging galaxies.

White dwarfs represent the end stage for 97% of stars, making precise parameter measurement crucial for understanding stellar evolution. Traditional estimation methods involve fitting spectra or photometry, which require high-quality data. In recent years, machine learning has played a crucial role in processing spectral data due to its speed, automation, and accuracy. However, two common issues have been identified. First, most studies rely on data with high signal-to-noise ratios (SNR > 10), leaving many poor-quality datasets underutilized. Second, existing machine learning models, primarily based on convolutional networks, recurrent networks, and their variants, cannot simultaneously capture both the spatial and sequential information of spectra. To address these challenges, we designed the Estimator Network (EstNet), an advanced algorithm integrating multiple techniques, including Residual Networks, Squeeze and Excitation Attention, Gated Recurrent Units, Adaptive Loss, and Monte-Carlo Dropout Layers. We conducted parameter estimation on 5,965 poor-quality white dwarf spectra (R~1800, SNR~1.17), achieving average percentage errors of 14.86% for effective temperature and 3.97% for surface gravity. These results are significantly superior to other mainstream algorithms and consistent with the outcomes of traditional theoretical spectrum fitting methods. In the future, our algorithms will be applied for large-scale parameter estimation on the Chinese Space Station Telescope and the Large Synoptic Survey Telescope.

While the secondary in a binary asteroid plays an important role in the precession of the mutual orbit, this role has not been thoroughly studied. Given the complex spin-orbit coupled dynamics in binary asteroids, we use a numerical approach to study the relationship between the secondary's shape and spin and the apsidal precession rate of the orbit. Using this approach in conjunction with observations of Didymos, we find it is likely that Dimorphos was significantly reshaped as a result of the DART impact, with its new shape more elongated than the pre-impact shape. Finally, we show that non-principal axis rotation of the secondary can lead to a chaotic evolution of the longitude of the periapsis.

We update the PlanetS catalog of transiting planets with precise and robust mass and radius measurements and use this comprehensive catalog to explore Mass-Radius (M-R) diagrams. On one hand, we propose new M-R relationships to separate exoplanets into three populations. On the other hand, we explore the transition in radius and density between Super-Earths and Sub-Neptunes around M-dwarfs and compare it with those orbiting K- and FG-dwarfs. Using Kernel Density Estimation method with a re-sampling technique, we estimate the normalized density and radius distributions, revealing connections between observations and theories on composition, internal structure, formation, and evolution of these exoplanets orbiting different spectral types. Firstly, the substantial 30% increase in the number of well-characterized exoplanets orbiting M-dwarfs compared with previous studies shows us that there is no clear gap in either composition or radius between Super-Earths and Sub-Neptunes. The "water-worlds" around M-dwarfs cannot correspond to a distinct population, their bulk density and equilibrium temperature can be interpreted by several different internal structures and compositions. The continuity in the fraction of volatiles in these planets suggests a formation scenario involving planetesimal or hybrid pebble-planetesimal accretion. Moreover, we find that the transition between Super-Earths and Sub-Neptunes appears to happen at different mass (and radii) depending on the spectral type of the star. The maximum mass of Super-Earths seems to be close to 10 M$_\oplus$ for all spectral types, but the minimum mass of Sub-Neptunes increases with the star's mass. This effect also contributes to the fading of the radius valley for M-planets compared to FGK-planets. While Sub-Neptunes are less common around M-dwarfs, smaller ones exhibit lower density than their equivalents around FGK-dwarfs.

Contemporary cosmic microwave background (CMB) experiments typically have observing bands covering the range 20 - 800 GHz. Certain science goals, including the detection of $\mu$-type distortions to the CMB spectrum and the characterization of low-frequency foregrounds, benefit from extended low-frequency coverage, but the standard CMB detector technology is not trivially adaptable to radio wavelengths. We propose using the upcoming Square Kilometer Array (SKA) as a CMB experiment, exploiting the immense raw sensitivity of SKA, in particular in single-dish mode, to measure medium-to-large-angular-scale modes of the CMB at radio wavelengths. As a worked example, we forecast the power of SKA combined with the upcoming LiteBIRD CMB space mission to constrain primordial non-Gaussianity through measurements of the correlation between anisotropies in the CMB $\mu$-distortion, temperature, and $E$-mode polarization fields. We find that adding SKA data significantly improves the constraints on $f_\textrm{nl}$, even for spatially varying low-frequency foregrounds.