Papers for Tuesday, Sep 24 2024

Currently the only method to establish the prevalence of particles, space debris or meteoroids, sized between 1 micrometre and a few centimetres, in Earth orbit is by instruments or witness plates dedicated to in-situ detection. Derived usable datasets are remarkably scarce and generally only cover a short period of time and a single orbital region. Nonetheless, space environment models use those limited datasets as anchor points to extrapolate results to the entirety of Earth orbit, from the beginning of the space age to decades into the future. Here we present a readout of over 20 years of DEBris In orbit Evaluator 1 (DEBIE-1), an in-situ detector that was launched in October 2001, providing the longest continuous set of measurements available to date. The dataset has not been used in the generation of space environment models and hence provides a first independent source for the detection of environment changing events and for the calibration of long term evolution models.

Transient surveys routinely detect supernovae (SNe) without obvious host galaxies. To understand the demographics of these "hostless" SNe and to constrain the possible host properties, we identify 161 SNe reported to the Transient Name Server since 2016 that do not have hosts cataloged from pre-explosion wide-field galaxy surveys. Using forced aperture photometry, we detect excess flux around only 56 of these SNe. Both thermonuclear and core-collapse (CC) SNe are present in our sample. Compared to flux-limited SNe samples with known hosts, superluminous supernovae (SLSNe), particularly hydrogen-deficient SLSNe, are over-represented here relative to all other SNe types; among CC SNe, there is also a higher fraction of interacting SNe than non-interacting. On the low-luminosity side, seven SNe have host absolute magnitude upper limits fainter than M_g=-12, about 1 per cent of the Small Magellanic Cloud's luminosity; the faintest limits are close to the luminosity of globular clusters or ultra-faint dwarf galaxies (M_g~-8). Fitting multi-band forced photometry, 11 SNe have host stellar masses <10^6 Msun assuming quiescent hosts, and 13 SNe have host stellar masses <10^5 Msun assuming star-forming hosts. The spatial distribution of hostless SNe indicates that the majority are not associated with known galaxy groups and clusters, ruling out intracluster stellar light as the primary contributor of such SNe. Hostless Type Ia SNe tend to be more luminous and slow-fading than SNe Ia with known host galaxies, implying a hidden population of low-mass and star-forming hosts. We conclude that any undetected host galaxies are likely star-forming dwarfs in the field.

This study explores the formation and implications of mini-active galactic nuclei (mAGN) disks around intermediate-mass black holes (IMBHs) embedded in gas-rich globular/nuclear clusters (GCs). We examine the parameter space for stable mAGN disks, considering the influence of IMBH mass, disk radius, and gas density on disk stability. The dynamics of stars and black holes within the mAGN disk are modeled, with a focus on gas-induced migration and gas dynamical friction. These dynamical processes can lead to several potentially observable phenomena, including the alignment of stellar orbits into the disk plane, the enhancement of gravitational wave mergers (particularly IMRIs and EMRIs), and the occurrence of mili/centi-tidal disruption events (mTDEs/cTDEs) with unique observational signatures. We find that gas hardening can significantly accelerate the inspiral of binaries within the disk, potentially leading to a frequency shift in the emitted gravitational waves. Additionally, we explore the possibility of forming accreting IMBH systems from captured binaries within the mAGN disk, potentially resulting in the formation of ultraluminous X-ray sources (ULXs). The observational implications of such accreting systems, including X-ray emission, optical signatures, and transient phenomena, are discussed. Furthermore, we investigate the possibility of large-scale jets emanating from gas-embedded IMBHs in GCs. While several caveats and uncertainties exist, our work highlights the potential for mAGN disks to provide unique insights into IMBH demographics, accretion physics, and the dynamics of GCs.

We analyse high-resolution ESPaDOnS/CFHT spectra of 20 very metal-poor stars ([Fe/H]~$<-2.0$) in the solar neighbourhood (within $\sim2$ kpc) selected to be on planar orbits (maximum height of $<4$ kpc). Targets include 11 prograde and 9 retrograde stars, spanning a wide range of eccentricities ($0.20-0.95$). Their chemical abundances are consistent with those observed in the Galactic halo but show a smaller spread, with no notable difference between progrades and retrogrades. This suggests a common chemical evolution and likely a shared formation site (except for one star). In this case, chemical evolution models indicate that the formation site would have had a baryonic mass of $\sim1.4\times10^9\msun$, similar to classical dwarf galaxies. High-energy supernovae and hypernovae are needed to reproduce the [X/Fe] up to the Fe-peak, while fast-rotating massive stars and neutron star merger events explain the [X/Fe] of the neutron-capture elements. The absence of Type Ia supernova signatures suggests a star formation duration of $<1$Gyr. Cosmological zoom-in simulations support the scenario that an in-plane infall of a single system could disperse stars over a wide range of angular momenta during the early Galactic assembly. We propose that these stars originated in a proto-Galactic building block, which we name \textit{Loki}. Less likely, if progrades and retrogrades formed in two different systems, their chemical evolution must have been very similar, with a combined baryonic mass twice that of a single system. The low number of targets and their limited metallicity range prevent us to exclude if these stars share a common progenitor with other detected structures, like GSE. A comparison (primarily [$\alpha$/Fe]) with other VMPs moving in planar orbits suggests multiple systems contributed to the Galactic planar population, presenting some differences in their kinematical parameters.

Atmospheric composition of exoplanets is often considered as a probe of the planet's formation condition. How exactly the initial chemical memory may be altered from the birth to the final state of the planet, however, remains unknown. Here, we develop a simple model of pollution of planetary atmosphere by the vaporization of infalling planetesimal of varying sizes and composition (SiO$_2$ inside 1 au and H$_2$O outside 1 au), following their trajectory and thermal evolution through the upper advective and radiative layers of a sub-Neptune class planet during the late stage of disk evolution. We vary the rate of pollution by changing the solid content of the disk and by dialing the level of disk gas depletion which in turn determines the rate of planetary migration. We find that pollution by silicate grains will always be limited by the saturation limit set by the thermal state of the atmosphere. By contrast, pollution by water ice can lead to $\sim$2--4 orders of magnitude variation in the atmospheric water mass fraction depending on the solid and gas content of the disk. Both cases suggest that post-formation pollution can erase the initial compositional memory of formation. Post-formation pollution can potentially transform sub-Neptunes with H/He-dominated envelope that initially formed beyond the iceline to waterworlds (water-enriched envelope) when the disk gas is depleted by $\gtrsim$2 orders of magnitude, allowing gentle migration. We additionally discuss the expected C/O ratio profile under pollution by water and refractory carbon species.

We use a neural network model and ALMA observations of HCN and HNC to constrain the physical conditions, most notably the cosmic-ray ionization rate (CRIR, zeta), in the Central Molecular Zone (CMZ) of the starburst galaxy NGC 253. Using output from the chemical code UCLCHEM, we train a neural network model to emulate UCLCHEM and derive HCN and HNC molecular abundances from a given set of physical conditions. We combine the neural network with radiative transfer modeling to generate modeled integrated intensities, which we compare to measurements of HCN and HNC from the ALMA Large Program ALCHEMI. Using a Bayesian nested sampling framework, we constrain the CRIR, molecular gas volume and column densities, kinetic temperature, and beam-filling factor across NGC 253's CMZ. The neural network model successfully recovers UCLCHEM molecular abundances with about 3 percent error and, when used with our Bayesian inference algorithm, increases the parameter inference speed tenfold. We create images of these physical parameters across NGC 253's CMZ at 50 pc resolution and find that the CRIR, in addition to the other gas parameters, is spatially variable with zeta a few times 10^{14} s^{-1} at greater than 100 pc from the nucleus, increasing to zeta greater than 10^{-13} s^{-1} at its center. These inferred CRIRs are consistent within 1 dex with theoretical predictions based on non-thermal emission. Additionally, the high CRIRs estimated in NGC 253's CMZ can be explained by the large number of cosmic-ray-producing sources as well as a potential suppression of cosmic-ray diffusion near their injection sites.

SPT2349$-$56 is a protocluster discovered in the 2500 deg$^2$ South Pole Telescope (SPT) survey. In this paper, we study the kinematics of the galaxies found in the core of SPT2349$-$56 using high-resolution (1.55 kpc spatial resolution at $z = 4.303$) redshifted [CII] 158-$\mu$m data. Using the publicly available code 3D-Barolo, we analyze the seven far-infrared (FIR) brightest galaxies within the protocluster core. Based on conventional definitions for the detection of rotating discs, we classify six sources to be rotating discs in an actively star-forming protocluster environment, with weighted mean $V_{\mathrm{rot}}/\sigma_{\mathrm{disp}} = 4.5 \pm 1.3$. The weighted mean rotation velocity ($V_{\mathrm{rot}}$) and velocity dispersion ($\sigma_{\mathrm{disp}}$) for the sample are $ 357.1 \pm 114.7$ km s$^{-1}$ and $43.5 \pm 23.5$ km s$^{-1}$, respectively. We also assess the disc stability of the galaxies and find a mean Toomre parameter of $Q_\mathrm{T} = 0.9 \pm 0.3$. The galaxies show a mild positive correlation between disc stability and dynamical support. Using the position-velocity maps, we find that five sources further classify as disturbed discs, and one classifies as a strictly rotating disc. Our sample joins several observations at similar redshift with high $V/\sigma$ values, with the exception that they are morphologically disturbed, kinematically rotating and interacting galaxies in an extreme protocluster environment.

Ultrahot Jupiters (UHJs), being the hottest class of exoplanets known, provide a unique laboratory for testing atmospheric interactions with internal planetary magnetic fields at a large range of temperatures. Thermal ionization of atmospheric species on the dayside of these planets results in charged particles becoming embedded in the planet's mostly neutral wind. The charges will resist flow across magnetic field lines as they are dragged around the planet and ultimately alter the circulation pattern of the atmosphere. We model this process to study this effect on high resolution emission and transmission spectra in order to identify observational signatures of the magnetic circulation regime that exist across multiple UHJs. Using a state-of-the-art kinematic MHD/active drag approach in a 3D atmospheric model, we simulate three different ultrahot Jupiters with and without magnetic effects. We post-process these models to generate high resolution emission and transmission spectra and explore trends in net Doppler shift as a function of phase. In emission spectra, we find that the net Doppler shift before and after secondary eclipse can be influenced by the presence of magnetic drag and wavelength choice. Trends in transmission spectra show our active drag models consistently produce a unique shape in their Doppler shift trends that differs from the models without active drag. This work is a critical theoretical step to understanding how magnetic fields shape the atmospheres of UHJs and provides some of the first predictions in high resolution spectroscopy for observing these effects.

SNIa and baryonic acoustic oscillations cosmic probes require performing measurements on astrophysical objects. Such objects are mostly located in overdense regions of the universe, where matter has grouped together. Measurements are thus biased in the sense that they collect information on a particular part of the universe only. Also, given the fact that the observations present a dynamics corresponding to an overall density evolving over time according to what we would expect for the density in overdense regions, we may wonder if the these cosmic probes measure the dynamics of such regions instead of the one of the global universe. In order to investigate this possibility, we improve here the two-regions model that was developed recently, to take into account the existence of overdense and underdense regions. Applying this model on SNIa and BAO cosmic probes, we highlight the reasons for which in practice we observe an accelerated expansion of the universe, involving an apparent cosmological constant.

Radio surveys typically sample extragalactic sources in higher redshift regimes than is typical for optical surveys, resulting in many radio sources not having a detected optical counterpart. Over the next decade the Legacy Survey of Space and Time (LSST) will be performing the deepest ($i < 26.4\,$mag) wide-area optical survey to date increasing the fraction of radio sources for which we have optical data. In this Research Note we use the Hyper Suprime-Cam survey to analyse how the fraction of radio sources in the Very Large Array Sky Survey (VLASS) with optical detections varies as a function of $i$-band magnitude and extrapolate to predict the number of optical counterparts we expect LSST to detect. Assuming a final VLASS point source depth of $S_{3\,\text{GHz}}\lesssim350\,\mu$Jy, we expect LSST to identify optical counterparts to $\sim 10^6$ radio sources in VLASS.

We present age estimates for over 8100 sub-giant branch (SGB) stars in Omega Centauri ($\omega$ Cen) to study its star formation history. Our large data set, which combines multi-wavelength HST photometry with MUSE metallicities, provides an unprecedented opportunity to measure individual stellar ages. We do this by fitting each star's photometry and metallicity with theoretical isochrones, that are embedded with an empirical [C+N+O]-[Fe/H] relation specifically for $\omega$ Cen. The bulk of the stars have ages between 13 and 10 Gyr, with the mean stellar age being 12.08$\pm$0.01 Gyrs and the median age uncertainty being 0.68 Gyrs. From these ages we construct the most complete age-metallicity relation (AMR) for $\omega$ Cen to date. We find that the mean age of stars decreases with increasing metallicity and find two distinct streams in the age-metallicity plane, hinting at different star formation pathways. We derive an intrinsic spread in the ages of 0.75$\pm$0.01 Gyr for the whole cluster, with the age spread showing a clear increase with metallicity. We verify the robustness of our age estimations by varying isochrone parameters and constraining our systematics. We find the C+N+O relation to be the most critical consideration for constraining the AMR. We also present the SGB chromosome map with age information. In the future, these stellar ages could be combined with chemical abundances to study age differences in subpopulations, and uncover the chemical evolution history of this massive nuclear star cluster.

Diffuse synchrotron emission in the form of radio halos and radio relics probe the existence of relativistic electrons and magnetic fields in galaxy clusters. These nonthermal components are generated from the dissipation of kinetic energy released by turbulence and shocks injected in the intracluster medium (ICM) during the large-scale structure formation process. By using the deepest images ever obtained on a galaxy cluster at low-frequency (72 h LOFAR-HBA + 72 h LOFAR-LBA), in arXiv:2211.01493 we provided an unprecedented view of the distribution of relativistic electrons and magnetic fields in the far outskirts of Abell 2255. In particular, we observed pervasive radio emission that fills the entire cluster volume and extends up to the cluster virial radius, reaching a maximum projected linear size of 5 Mpc. By combining radio and X-ray observations with advanced numerical simulations, we estimated the magnetic field and energy budget associated to turbulent motions at such large distances from the cluster center. Our results suggest an efficient transfer of kinetic energy into nonthermal components in the extremely diluted cluster outskirts. In the past two years, the total LOFAR-HBA observation time on Abell 2255 has increased to 336 hours. The analysis of this ultra-deep dataset aims to further advance our understanding of relativistic electrons and magnetic fields in cluster peripheries.

At sufficiently high pressures (~Mbar) and low temperatures (~1e3-1eK), hydrogen and helium become partly immiscible. Interpretations of Jupiter and Saturn's magnetic fields favor the existence of a statically stable layer near the Mbar pressure level. From experimental and computational data for the hydrogen-helium phase diagram we find that moist convection and diffusive convection are inhibited, implying a stable helium rain layer in both Jupiter and Saturn. However, we find a significant difference in terms of structure and evolution: In Jupiter, helium settling leads to a stable yet super-adiabatic temperature gradient that is limited by conductive heat transport. The phase separation region should extend only a few tens of kilometers instead of thousands in current-day models, and be characterized by a sharp increase of the temperature of about 500K for standard phase separation diagrams. In Saturn, helium rains occurs much deeper, implying a larger helium flux relative to planetary mass. We find that the significant latent heat associated with helium condensation implies that a large fraction, perhaps close to 100%, of the planet's intrinsic heat flux may be locally transported by the sinking helium droplets. This implies that Saturn may possess a much more extended helium-rain region. This also accounts, at least qualitatively, for the differences in strength and characteristics of the magnetic fields of the two planets. Dedicated models of magnetic field generation in both planets may offer observational constraints to further refine these findings.

We revisit the small field double-well inflationary model and investigate its observational viability in light of the current Cosmic Microwave Background data. In particular, considering scenarios with minimal and nonminimal coupling between the inflaton field and the Ricci scalar, we perform a Monte Carlo Markov chain analysis to probe the model's parameter space. We also investigate the consequences of the cosmological results in the canonical type-I seesaw mechanism context and obtain constraints on the vacuum expectation value of the inflaton field, together with the amplitude of the coupling to gravity in the nonminimal case. We employ a Bayesian procedure to compare the model's predictions with the Starobinsky inflationary scenario and find a strong statistical preference for the latter against the minimal and nonminimal coupled double-well scenario.

We introduce the construction of polarized intensity cubes $P$(RA, Dec, $\Phi$) and their visualization as movies, as a powerful technique for interpreting Faraday structure. $P$ is constructed from maps of peak polarized intensity P(RA, Dec) with their corresponding Faraday depth maps $\Phi$(RA, Dec). We illustrate the extensive scientific potential of such visualizations with a variety of science use cases from ASKAP and MeerKAT, presenting models that are consistent with the data but not necessarily unique. We demonstrate how one can, in principle, distinguish between cube structures which originate from unrelated foreground screens from those due to magnetized plasmas local to the emitting source. Other science use cases illustrate how variations in the local $n_e$ $B$, and line-of-sight distance to the synchrotron emitting regions can be distinguished using Faraday rotation. We show, for the first time, how the line-of-sight orientation of AGN jets can be determined. We also examine the case of M87 to show how internal jet magnetic field configurations can be identified, and extend earlier results. We recommend using this technique to re-evaluate all previous analyses of polarized sources that are well-resolved both spatially and in Faraday depth. Recognizing the subjective nature of interpretations at this early stage, we also highlight the need and utility for further scientific and technical developments.

Recently the Nereides nebula was discovered through deep optical emission line observations and was classified as a supernova remnant (SNR) candidate, G107.7-5.1. Since very little is known about this SNR we look at several archival data sets to better understand the environment and properties of the object. We present a detailed analysis of the gamma-ray emission detected by the \textit{Fermi} Large Area Telescope in the region of the nebula. A model of the non-thermal emission is presented which allows us to derive the particle distribution responsible for the gamma rays. We also use molecular gas and atomic hydrogen observations to try to constrain the source age and distance. An extended (~2deg) GeV source coincident with the location of the nebula is found. The non-thermal emission has a hard spectrum and is detected up to about 100 GeV, confirming the SNR nature of this object. The GeV properties of G107.7-5.1 are similar to those of other SNRs such as G150.3+4.5, and likely expands in a relatively low-density medium. The Nereides nebula is one more example of a growing population of dim SNRs detected at high energies. A simple leptonic model is able to account for the gamma-ray emission. Standard SNR evolutionary models constrain the age to be in the 10-50 kyr range, which is consistent with estimates of the maximum particle energy obtained from GeV observations. However, more detailed observations of the source should be carried out to better understand its properties.

Accurate solar flare prediction is crucial due to the significant risks that intense solar flares pose to astronauts, space equipment, and satellite communication systems. Our research enhances solar flare prediction by utilizing advanced data preprocessing and classification methods on a multivariate time series-based dataset of photospheric magnetic field parameters. First, our study employs a novel preprocessing pipeline that includes missing value imputation, normalization, balanced sampling, near decision boundary sample removal, and feature selection to significantly boost prediction accuracy. Second, we integrate contrastive learning with a GRU regression model to develop a novel classifier, termed ContReg, which employs dual learning methodologies, thereby further enhancing prediction performance. To validate the effectiveness of our preprocessing pipeline, we compare and demonstrate the performance gain of each step, and to demonstrate the efficacy of the ContReg classifier, we compare its performance to that of sequence-based deep learning architectures, machine learning models, and findings from previous studies. Our results illustrate exceptional True Skill Statistic (TSS) scores, surpassing previous methods and highlighting the critical role of precise data preprocessing and classifier development in time series-based solar flare prediction.

Aiming to capture the formation and eruption of flux ropes (FRs) in the source active regions (ARs), we simulate the coronal magnetic field evolution of the AR 11429 employing the time-dependent magneto-friction model (TMF). The initial field is driven by electric fields that are derived from time-sequence photospheric vector magnetic field observations by invoking ad-hoc assumptions. The simulated magnetic structure evolves from potential to twisted fields over the course of two days, followed by rise motion in the later evolution, depicting the formation of FR and its slow eruption later. The magnetic configuration resembles an inverse S-sigmoidal structure, composed of a potential field enveloping the inverse J-shaped fields that are shared past one another and a low lying twisted field along the major PIL. To compare with observations, proxy emission maps based on averaged current density along the field lines are generated from the simulated field. These emission maps exhibit a remarkable one-to-one correspondence with the spatial characteristics in coronal EUV images, especially the filament-trace supported by the twisted magnetic field in the south-west subregion. Further, the topological analysis of the simulated field reveals the co-spatial flare ribbons with the quasi-separatrix layers, which is consistent with the standard flare models; therefore, the extent of the twist and orientation of the erupting FR is indicated to be the real scenario in this case. The TMF model simulates the coronal field evolution, correctly capturing the formation of the FR in the observed time scale and the twisted field generated from these simulations serve as the initial condition for the full MHD simulations.

The Observatorio del Roque de los Muchachos will host the northern site of the Cherenkov Telescope Array Observatory (CTAO), in an area about 200 m below the mountain rim, where the optical telescopes are located. The site currently hosts the MAGIC Telescopes, which have gathered a unique series of 20 years of weather data. We use advanced profile likelihood methods to determine seasonal cycles, the occurrence of weather extremes, weather downtime, and long-term trends correctly taking into account data gaps. The fractality of the weather data is investigated by means of multifractal detrended fluctuation analysis. The data are published according to the Findable, Accessible, Interoperable, and Reusable (FAIR) principles. We find that the behaviour of wind and relative humidity show significant differences compared to the mountain rim. We observe an increase in temperature of $0.55\pm0.07\mathrm{(stat.)}\pm0.07\mathrm{(syst.)}^\circ C$/decade, the diurnal temperature range of $0.13\pm0.04\mathrm{(stat.)}\pm0.02\mathrm{(syst.)}^\circ C$/decade (accompanied by an increase of seasonal oscillation amplitude of $\Delta C_m=0.29\pm0.10\mathrm{(stat.)}\pm0.04\mathrm{(syst.)}^\circ C$/decade) and relative humidity of $4.0\pm0.4\mathrm{(stat.)}\pm1.1\mathrm{(syst.)}$%/decade, and a decrease in trade wind speeds of $0.85\pm0.12\mathrm{(stat.)}\pm0.07\mathrm{(syst.)}$(km/h)/decade. The occurrence of extreme weather, such as tropical storms and long rains, remains constant over time. We find a significant correlation of temperature with the North Atlantic Oscillation Index and multifractal behaviour of the data. The site shows a weather-related downtime of 18.5%-20.5%, depending on the wind gust limits employed. No hints are found of a degradation of weather downtime under the assumption of a linear evolution of environmental parameters over time.

We present a comprehensive analysis of the X-ray spectral properties of 198 newly identified active galactic nuclei (AGNs), leveraging archival data from the {\it Chandra} X-ray Observatory. All these AGNs exhibit a powerlaw spectral signature spanning a broad energy range of $0.5-7.0$ keV, characterized by the photon index ($\Gamma$) values ranging from $0.3^{+0.16}_{-0.14}$ to $2.54^{+0.14}_{-0.13}$. Particularly, 76 of these AGNs display discernible levels of intrinsic absorption, after considering the Galactic absorption. The column densities associated with this local absorption ($n_{\rm H}^{\rm local}$) are within a range of $\sim 10^{19} - 10^{22}\ {\rm cm^{-2}}$. We study the cosmological evolution of AGNs using the variation of $n_{\rm H}^{\rm local}$ and $\Gamma$ with their estimated redshift. The intrinsic spectral signature did not reveal any significant cosmological evolution; however, a deficit of hard sources at high redshift is possibly intrinsic. Our sample covers several decades of broadband intrinsic luminosity ($L_{\rm B}^{\rm intr}$) ranging from $4.59^{+0.41}_{-0.41} \times 10^{42}$ to $2.4^{+0.12}_{-0.12} \times 10^{46}\, {\rm erg~s}^{-1}$ with peak at 1.84 redshift. We also investigate the hardness-luminosity diagram (HLD) to further probe the AGNs. We conduct a sanity check by applying our findings to known AGNs, and the results are consistent with our observations.

Since the first direct detection of gravitational waves by the LIGO--Virgo collaboration in 2015, the size of the gravitational-wave transient catalog has grown to nearly 100 events, with more than as many observed during the ongoing fourth observing run. Extracting astrophysical/cosmological information from these observations is a hierarchical Bayesian inference problem. GWPopulation is designed to provide simple-to-use, robust, and extensible tools for hierarchical inference in gravitational-wave astronomy/cosmology. It has been widely adopted for gravitational-wave astronomy, including producing flagship results for the LIGO-Virgo-KAGRA collaborations. While designed to work with observations of compact binary coalescences, GWPopulation may be available to a wider range of hierarchical Bayesian inference problems.

Accretion physics studies the process of gravitational capture of ambient matter by massive stars. The background processes are very challenging to observe and measure due to the extreme conditions in these systems. Numerical simulations play a crucial role in accretion physics because they provide the only practical method to model the complex processes occurring in accretion disks. In this review, we outline different branches of numerical simulations, such as hydrodynamic simulations, magnetohydrodynamic simulations, and Monte-Carlo simulations, and their methodology, and we discuss possible implications for modeling accretion physics around black holes, neutron stars, and protoplanetary disks.

Growing observational evidence suggests that enhanced mass loss from the progenitors of core-collapse supernovae (SNe) is common during $\sim1$ yr preceding the explosion, creating an optically thick circum-stellar medium (CSM) shell at $\sim10^{14.5}$ cm radii. We show that if such mass loss is indeed common, then the breakout of the SN shock through the dense CSM shell produces a neutrino flux that may account for a significant fraction of the observed $\gtrsim10$ TeV neutrino background. The neutrinos are created within a few days from the explosion, during and shortly after the shock breakout, which produces also large UV (and later X-ray) luminosity. The compact size and large UV luminosity imply a pair production optical depth of $\sim10^4$ for $>100$ GeV photons, naturally accounting for the lack of a high-energy gamma-ray background accompanying the neutrino background. SNe producing $>1$ neutrino event in a 1 km$^2$ detector are expected at a rate of $\lesssim0.1$/yr. A quantitative theory describing the evolution of the electromagnetic spectrum during a breakout, as the radiation-mediated shock is transformed into a collisionless one, is required to enable (i) using data from upcoming surveys that will systematically detect large numbers of young, $<1$ d old SNe, to determine the pre-explosion mass loss history of the SN progenitor population, and (ii) a quantitative determination of the neutrino luminosity and spectrum.

We present the first results of an extensive spectroscopic survey of directly imaged planet host stars. The goal of the survey is the measurement of stellar properties and abundances of 15 elements (including C, O, and S) in these stars. In this work, we present the analysis procedure and the results for an initial set of five host stars, including some very well-known systems. We obtain C/O ratios using a combination of spectral modeling and equivalent width measurements for all five stars. Our analysis indicates solar C/O ratios for HR 8799 (0.59 $\pm$ 0.11), 51 Eri (0.54 $\pm$ 0.14), HD 984 (0.63 $\pm$ 0.14), and GJ 504 (0.54 $\pm$ 0.14). However, we find a super-solar C/O (0.81 $\pm$ 0.14) for HD 206893 through spectral modeling. The ratios obtained using the equivalent width method agree with those obtained using spectral modeling but have higher uncertainties ($\sim$0.3 dex). We also calculate the C/S and O/S ratios, which will help us to better constrain planet formation, especially once planetary sulfur abundances are measured using JWST. Lastly, we find no evidence of highly elevated metallicities or abundances for any of our targets, suggesting that a super metal-rich environment is not a prerequisite for large, widely separated gas planet formation. The measurement of elemental abundances beyond carbon and oxygen also provides access to additional abundance ratios, such as Mg/Si, which could aid in further modeling of their giant companions.

Large near-future galaxy surveys offer sufficient statistical power to make our cosmology analyses data-driven, limited primarily by systematic errors. Understanding the impact of systematics is therefore critical. We perform an end-to-end analysis to investigate the impact of some of the systematics that affect large-scale structure studies by doing an inference analysis using simulated density maps with various systematics; these include systematics caused by photometric redshifts (photo-$z$s), Galactic dust, structure induced by the telescope observing strategy and observing conditions, and incomplete covariance matrices. Specifically, we consider the impacts of incorrect photo-$z$ distributions (photometric biases, scatter, outliers; spectroscopic calibration biases), dust map resolution, incorrect dust law, selecting none or only some contaminant templates for deprojection, and using a diagonal covariance matrix instead of a full one. We quantify the biases induced by these systematics on cosmological parameter estimation using tomographic galaxy angular power spectra, with a focus on identifying whether the maximum plausible level of each systematic has an adverse impact on the estimation of key cosmological parameters from a galaxy clustering analysis with Rubin Observatory Legacy Survey of Space and Time (LSST). We find photo-$z$ systematics to be the most pressing out of the systematics investigated, with spectroscopic calibration biases leading to the greatest adverse impact while helpfully being flagged by a high $\chi^2$ value for the best fit model. Larger-than-expected photo-$z$ scatter, on the other hand, has a significant impact without necessarily indicating a poor fit. In contrast, in the analysis framework used in this work, biases from observational systematics and incomplete covariance matrices are comfortably subdominant.

The recently discovered bright transient black hole candidate Swift J1727.8-1613 is studied in a broad energy range ($0.5-79$ keV) using combined NICER and NuSTAR data on 29 August 2023. A promonient type-C Quasi-Periodic Oscillation (QPO) at $0.89 \pm 0.01$ Hz with its harmonic was observed in NICER data of $0.5-10$ keV. Interestingly, the harmonic becomes weaker in the lower energy bands ($0.5-1$ & $1-3$ keV). We also report the first detection of a soft time-lag of $0.014 \pm 0.001$ s at the QPO frequency between harder ($3-10$ kev) and softer ($0.5-3$ keV) band photons observed with the NICER/XTI instrument. This indicates that the inclination of the accretion disk in the binary system might be high. From the detailed spectral analysis with the relxill reflection model, we found the disk inclination angle of source to be $\sim 85^\circ$. We discuss how the accretion flow configuration inferred from spectral analysis can help us to understand the origin of QPOs and soft lag in this source.

Large Array of imaging atmospheric Cherenkov Telescope (LACT) is an array of 32 Cherenkov telescopes with 6-meter diameter mirrors to be constructed at the LHAASO site. In this work, we present a study on the layout optimization and performance analysis of LACT. We investigate two observation modes: large zenith angle observations for ultra-high energy events and small zenith angle observations for lower energy thresholds. For large zenith angles (60°), simulations show that an 8-telescope subarray can achieve an effective area of $3 ~\rm km^2$ and excellent angular resolution. For small zenith angles, we optimize the layout of 4-telescope cells and the full 32-telescope array. The threshold of the full array is about $200~\rm GeV$, which is particularly crucial for studying transient phenomena, including gamma-ray bursts (GRBs) and active galactic nuclei (AGNs). This study provides important guidance for the final LACT layout design and performance estimates under different observational conditions, demonstrating LACT's potential for deep observations of ultra-high energy \gray sources and morphological studies of PeVatrons, as well as time-domain \gray astronomy.

HOPS 87 is a Class 0 protostellar core known to harbor an extremely young bipolar outflow and a hot corino. We report the discovery of localized, chemically rich regions near the bases of the two-lobe bipolar molecular outflow in HOPS 87 containing molecules such as H$_2$CO, $^{13}$CS, H$_2$S, OCS, and CH$_3$OH, the simplest complex organic molecule (COM). The locations and kinematics suggest that these localized features are due to jet-driven shocks rather than being part of the hot corino region encasing the protostar. The COM compositions of the molecular gas in these jet-localized regions are relatively simpler than those in the hot corino zone. We speculate that this simplicity is due to either the liberation of ice with a less complex chemical history or the effects of shock chemistry. Our study highlights the dynamic interplay between the protostellar bipolar outflow, disk, inner core environment, and the surrounding medium, contributing to our understanding of molecular complexity in solar-like young stellar objects.

In general, large solar flares are more efficient at accelerating high-energy electrons than microflares. Nonetheless, we sometimes observe microflares that accelerate electrons to high energies. We statistically characterize 39 microflares with strikingly hard spectra in the hard X-ray (HXR) range, which means that they are efficient in accelerating high-energy electrons. We refer to these events as "hard microflares." The statistical analysis is built upon spectral and imaging information from STIX, combined with EUV and magnetic field maps from SDO. The key observational result is that all hard microflares in this dataset have one of the footpoints rooted directly within a sunspot (either in the umbra or the penumbra). This clearly indicates that the underlying magnetic flux densities are large. For the events with the classic two-footpoints morphology, the absolute value of the mean line-of-sight magnetic flux density (and vector magnetic field strength) at the footpoint rooted within the sunspot ranges from 600 to 1800 G (1500 to 2500 G), whereas the outer footpoint measures from 10 to 200 G (100 to 400 G), therefore about 10 times weaker. Approximately 78% of hard microflares, which exhibited two HXR footpoints, have similar or even stronger HXR flux from the footpoint rooted within the sunspot. This contradicts the magnetic mirroring scenario. In addition, about 74% of the events could be approximated by a single loop geometry, demonstrating that hard microflares typically have a relatively simple morphology. We conclude that all hard microflares are rooted in sunspots, which implies that the magnetic field strength plays a key role in efficiently accelerating high-energy electrons, with hard HXR spectra associated with strong fields. This key result will allow us to further constrain our understanding of the electron acceleration mechanisms in flares and space plasmas.

During the fourth observing run, the LIGO-Virgo-KAGRA Collaboration reported the detection of a coalescing compact binary (GW230529$_{-}$181500) with component masses estimated at $2.5-4.5\, M_\odot$ and $1.2-2.0\, M_\odot$ with 90\% credibility. Given the current constraints on the maximum neutron star (NS) mass, this event is most likely a lower-mass-gap (LMG) black hole-neutron star (BHNS) binary. The spin magnitude of the BH, especially when aligned with the orbital angular momentum, is critical in determining whether the NS is tidally disrupted. An LMG BHNS merger with a rapidly spinning BH is an ideal candidate for producing electromagnetic counterparts. However, no such signals have been detected. In this study, we employ a detailed binary evolution model, incorporating new dynamical tide implementations, to explore the origin of BH spin in an LMG BHNS binary. If the NS forms first, the BH progenitor (He-rich star) must begin in orbit shorter than 0.35 days to spin up efficiently, potentially achieving a spin magnitude of $\chi_{\rm BH} > 0.3$. Alternatively, if a non-spinning BH (e.g., $M_{\rm BH} = 3.6\, M_\odot$) forms first, it can accrete up to $\sim 0.2\, M_\odot$ via Case BA mass transfer (MT), reaching a spin magnitude of $\chi_{\rm BH} \sim 0.18$ under Eddington-limited accretion. With a higher Eddington accretion limit (i.e., 10.0 $\Dot{M}_{\rm Edd}$), the BH can attain a significantly higher spin magnitude of $\chi_{\rm BH} \sim\,0.65$ by accreting approximately $1.0\, M_\odot$ during Case BA MT phase.

Standard definitions of habitability assume that life requires the presence of planetary gravity wells to stabilize liquid water and regulate surface temperature. Here the consequences of relaxing this assumption are evaluated. Temperature, pressure, volatile loss, radiation levels and nutrient availability all appear to be surmountable obstacles to the survival of photosynthetic life in space or on celestial bodies with thin atmospheres. Biologically generated barriers capable of transmitting visible radiation, blocking ultraviolet, and sustaining temperature gradients of 25-100 K and pressure differences of 10 kPa against the vacuum of space can allow habitable conditions between 1 and 5 astronomical units in the solar system. Hence ecosystems capable of generating conditions for their own survival are physically plausible, given the known capabilities of biological materials on Earth. Biogenic habitats for photosynthetic life in extraterrestrial environments would have major benefits for human life support and sustainability in space. Because the evolution of life elsewhere may have followed very different pathways from on Earth, living habitats could also exist outside traditional habitable environments around other stars, where they would have unusual but potentially detectable biosignatures.

The early time observations of Type Ia supernovae (SNe Ia) play a crucial role in investigating and resolving longstanding questions about progenitor stars and the explosion mechanisms of these events. Colors of supernovae (SNe) in the initial days after the explosion can help differentiate between different types of SNe. However, the use of true color information to identify SNe Ia at the early-time explosion is still in its infancy. The Multi-channel Photometric Survey Telescope (Mephisto) is a photometric survey telescope equipped with three CCD cameras, capable of simultaneously imaging the same patch of sky in three bands (\emph{u, g, i} or \emph{v, r, z}), yielding real-time colors of astronomical objects. In this paper, we introduce a new time-series classification tool named Mephisto Early Supernovae Ia Rapid Identifier (\emph{\texttt{Mesiri}}), which for the first time, utilizes real-time color information to distinguish early-time SNe Ia from core-collapse supernovae (CCSNe). \emph{\texttt{Mesiri}} is based on the deep learning approach and can achieve an accuracy of $96.75\pm0.79$\%, and AUC of $98.87\pm0.53$\% in case of single epoch random observation before the peak brightness. These values reach towards perfectness if additional data points on several night observations are considered. The classification with real-time color significantly outperforms that with pseudo-color, especially at the early time, i.e., with only a few points of observations. The architecture of BiLSTM shows the best performance than the others that have been tested in this work.

We investigate the potential of using gravitational wave (GW) signals from rotating core-collapse supernovae to probe the equation of state (EOS) of nuclear matter. By generating GW signals from simulations with various EOSs, we train machine learning models to classify them and evaluate their performance. Our study builds on previous work by examining how different machine learning models, parameters, and data preprocessing techniques impact classification accuracy. We test convolutional and recurrent neural networks, as well as six classical algorithms: random forest, support vector machines, naive Bayes, logistic regression, k-nearest neighbors, and eXtreme gradient boosting. All models, except naive Bayes, achieve over 90 per cent accuracy on our dataset. Additionally, we assess the impact of approximating the GW signal using the general relativistic effective potential (GREP) on EOS classification. We find that models trained on GREP data exhibit low classification accuracy. However, normalizing time by the peak signal frequency, which partially compensates for the absence of the time dilation effect in GREP, leads to a notable improvement in accuracy.

The Dark Energy Spectroscopic Instrument (DESI) collaboration measured a tight relation between the Hubble constant ($H_0$) and the distance to the Coma cluster using the fundamental plane (FP) relation of the deepest, most homogeneous sample of early-type galaxies. To determine $H_0$, we measure the distance to Coma by several independent routes each with its own geometric reference. We measure the most precise distance to Coma from 12 Type Ia Supernovae (SNe Ia) in the cluster with mean standardized brightness of $m_B^0=15.712\pm0.041$ mag. Calibrating the absolute magnitude of SNe Ia with the HST distance ladder yields $D_{\textrm Coma}=98.5\pm2.2$ Mpc, consistent with its canonical value of 95--100 Mpc. This distance results in $H_0=76.5 \pm 2.2$ km/s/Mpc from the DESI FP relation. Inverting the DESI relation by calibrating it instead to the Planck+$\Lambda$CDM value of $H_0=67.4$ km/s/Mpc implies a much greater distance to Coma, $D_{\textrm Coma}=111.8\pm1.8$ Mpc, $4.6\sigma$ beyond a joint, direct measure. Independent of SNe Ia, the HST Key Project FP relation as calibrated by Cepheids, Tip of the Red Giant Branch from JWST, or HST NIR surface brightness fluctuations all yield $D_{\textrm Coma}<$ 100 Mpc, in joint tension themselves with the Planck-calibrated route at $>3\sigma$. From a broad array of distance estimates compiled back to 1990, it is hard to see how Coma could be located as far as the Planck+$\Lambda$CDM expectation of $>$110 Mpc. By extending the Hubble diagram to Coma, a well-studied location in our own backyard whose distance was in good accord well before the Hubble Tension, DESI indicates a more pervasive conflict between our knowledge of local distances and cosmological expectations. We expect future programs to refine the distance to Coma and nearer clusters to help illuminate this new, local window on the Hubble Tension.

Continuous nanohertz gravitational waves from individual supermassive black hole binaries may be detectable with pulsar timing arrays. A novel search strategy is developed, wherein intrinsic achromatic spin wandering is tracked simultaneously with the modulation induced by a single gravitational wave source in the pulse times of arrival. A two-step inference procedure is applied within a state-space framework, such that the modulation is tracked with a Kalman filter, which then provides a likelihood for nested sampling. The procedure estimates the static parameters in the problem, such as the sky position of the source, without fitting for ensemble-averaged statistics such as the power spectral density of the timing noise, and therefore complements traditional parameter estimation methods. It also returns the Bayes factor relating a model with a single gravitational wave source to one without, complementing traditional detection methods. It is shown via astrophysically representative software injections in Gaussian measurement noise that the procedure distinguishes a gravitational wave from pure noise down to a characteristic wave strain of $h_0 \approx 2 \times 10^{-15}$. Full posterior distributions of model parameters are recovered and tested for accuracy. There is a bias of $\approx 0.3$ rad in the marginalised one-dimensional posterior for the orbital inclination $\iota$, introduced by dropping the so-called `pulsar terms'. Smaller biases $\lesssim 10 \%$ are also observed in other static parameters.

We present experiments to reproduce the characteristics of core-collapse supernovae with different stellar masses and initial explosion energies in the laboratory. In the experiments, shocks are driven in 1.2 atm and 1.9 atm xenon gas by laser with energy from 1600J to 2800J on the SGIII prototype laser facility. The average shock velocities and shocked densities are obtained from experiments. Experimental results reveal that higher laser energy and lower Xe gas density led to higher shock velocity, and lower Xe gas initial density has a higher compression. Modeling of the experiments using the 2D radiation hydrodynamic codes Icefire shows excellent agreement with the experimental results and gives the temperature. These results will contribute to time-domain astrophysical systems, such as gravitational supernovae, where a strong radiative shock propagates outward from the center of the star after the core collapses.

About 25\% -50\% of white dwarfs (WDs) are found to be polluted by heavy elements. It has been argued that the pollution could be caused by the tidal disruption of an approaching planet around the WD, during which a large number of clumps would be produced and would finally fall onto the WD. The reason that the planet approaches the WD is usually believed to be due to gravitational perturbations from another distant planet or stellar companion. However, the dynamics of the perturbation and the detailed partial disruption process are still poorly understood. In this study, we present an in-depth investigation of these issues. A triple system composed of a WD, an inner orbit planet, and an outer orbit planet is considered. The inner plant would be partially disrupted periodically in the long-term evolution. Fragments generated in the process are affected by the gravitational perturbations from the remnant planet, facilitating their falling toward the WD. The mass loss rate of the inner planet depends on both its internal structure and also on the orbital configuration of the planetary system.

The impact of solar-stellar activity on planetary environments is a topic of great interest within the Sun-Earth system as well as exoplanetary systems. In particular, extreme events such as flares and coronal mass ejections have a profound effect on planetary atmospheres. In May this year, a magnetic active region on the Sun (AR 13664) -- with a size exceeding hundred times that of Earth -- unleashed a large number of high energy X-class flares and associated mass ejections. The resulting Earth impact (geomagnetic storm) on May 10-11 was the strongest in the last two decades. We perform the first comprehensive analysis of the magnetic properties of the active region that spawned these flares and identify this to be a super active region with very rare physical characteristics. We also demonstrate how the rate of energization of the system is related to the flaring process. Our work illuminates how flare productive super active regions on the Sun and stars can be identified and what are their salient physical properties. Specifically, we put AR 13664 in historical context over the cumulative period of 1874 May-2024 June. We find that AR 13664 stands at 99.95 percentile in the distribution of area over 1874 May-2024 June, and at 99.10 percentile in terms of flux content among all ARs over the period 1996 April-2024 June. Our analysis indicates that five of its magnetic properties rank highest among all ARs recorded in SHARP data series during 2010 May-2024 June by the Solar Dynamic Observatory. Furthermore, we demonstrate that AR 13664 reached its most dynamic flare productive state following a rapid rate of rise of its flare-relevant parameters and that the X-class flares it spawned were more frequent near their peak values. Our analyses establish AR 13644 to be solar super active region and provide a paradigm for investigating their flare-relevant physical characteristics.

In this paper, we use hydrodynamic zoom-in simulations of Milky Way-type haloes to explore using dust as an observational tracer to discriminate between cold and warm dark matter universes. Comparing a cold and 3.5keV warm dark matter particle model, we tune the efficiency of galaxy formation in our simulations using a variable supernova rate to create Milky Way systems with similar satellite galaxy populations while keeping all other simulation parameters the same. Cold dark matter, having more substructure, requires a higher supernova efficiency than warm dark matter to achieve the same satellite galaxy number. These different supernova efficiencies create different dust distributions around their host galaxies, which we generate by post-processing the simulation output with the POWDERDAY codebase. Analysing the resulting dust in each simulation, we find $\sim$4.5 times more dust in our cold dark matter Milky Way halos compared with warm dark matter. The distribution of dust out to R$_{200\text{c}}$ is then explored, revealing that the warm dark matter simulations are noticeably less concentrated than their cold dark matter counterparts, although differences in substructure complicate the comparison. Our results indicate that dust is a possible unique probe to test theories of dark matter.

We present a method to refine photometric redshift galaxy catalogs by comparing their color-space matching with overlapping spectroscopic calibration data. We focus on cases where photometric redshifts (photo-$z$) are estimated empirically. Identifying galaxies that are poorly represented in spectroscopic data is crucial, as their photo-$z$ may be unreliable due to extrapolation beyond the training sample. Our approach uses a self-organizing map (SOM) to project a multi-dimensional parameter space of magnitudes and colors onto a 2-D manifold, allowing us to analyze the resulting patterns as a function of various galaxy properties. Using SOM, we compare the Kilo-Degree Survey bright galaxy sample (KiDS-Bright), limited to $r<20$ mag, with various spectroscopic samples, including the Galaxy And Mass Assembly (GAMA). Our analysis reveals that GAMA under-represents KiDS-Bright at its faintest ($r\gtrsim19.5$) and highest-redshift ($z\gtrsim0.4$) ranges, however no strong trends in color or stellar mass. By incorporating additional spectroscopic data from the SDSS, 2dF, and early DESI, we identify SOM cells where photo-$z$ are estimated suboptimally. We derive a set of SOM-based criteria to refine the photometric sample and improve photo-$z$ statistics. For the KiDS-Bright sample, this improvement is modest: exclusion of the least represented 20% of the sample reduces photo-$z$ scatter by less than 10%. We conclude that GAMA, used for KiDS-Bright photo-$z$ training, is sufficiently representative for reliable redshift estimation across most of the color space. Future spectroscopic data from surveys such as DESI should be better suited for exploiting the full improvement potential of our method.

Turbulent gas motions are expected to dominate the non-thermal energy budget of the intracluster medium (ICM). The measurement of pressure fluctuations from high angular resolution Sunyaev-Zel'dovich imaging opens a new avenue to study ICM turbulence, complementary to X-ray density fluctuation measures. We develop a methodological framework designed to optimally extract information on the ICM pressure fluctuation power spectrum statistics, and publicly release the associated software named PITSZI. We apply this tool to the NIKA data of the merging cluster MACSJ0717 to measure its pressure fluctuation power spectrum at high significance, and to investigate the implications for its nonthermal content. Depending on the choice of the radial pressure model and the details of the applied methodology, we measure an energy injection scale $L_{inj} \sim 800$ kpc. The power spectrum normalization corresponds to a characteristic amplitude reaching $A(k_{peak}) \sim 0.4$. These results are are obtained assuming that MACSJ0717 can be described as pressure fluctuations on top of a single (smooth) halo, and are dominated by systematics due to the choice of the radial pressure model. Using simulations, we estimate that fitting a radial model to the data can suppress the observed fluctuations by up to 50\%, while a poorly representative radial model can induce spurious fluctuations, which we also quantify. Assuming standard scaling relations between the pressure fluctuations and turbulence, we find that MACSJ0717 presents a turbulent velocity dispersion $\sigma_v \sim 1200$ km/s, a kinetic to kinetic plus thermal pressure fraction $P_{k} / P_{k+th} \sim 20\%$, and we estimate the hydrostatic mass bias to $b_{HSE} \sim 0.3-0.4$. Our results are in excellent agreement with alternative measurements from X-ray surface brightness fluctuations, and in agreement with the fluctuations being adiabatic in nature.

Previous studies have found that localized strengthening of the f mode recedes the emergence of active regions on the Sun by one day to three days. To help interpret these observations, we have performed nonlinear simulations of convection with imposed magnetic fields at different depths. We find that the f mode is strengthened when a super-equipartition magnetic field is imposed near the top of the domain. However, neither a magnetic field of equal strength near the bottom of the domain nor an equipartition magnetic field near the top of the domain have a significant effect. Our results suggest that the magnetic precursors of active regions are present near the surface of the Sun for much longer than would be expected if active regions were formed by flux tubes rising from deep within the convection zone. Application to observations should account for the fact that the effects we observe are transient.

We use the rms and lag spectra of the type-C quasi-periodic oscillation (QPO) to study the properties of the Comptonisation region (aka corona) during the low/hard and hard-intermediate states of the main outburst and reflare of MAXI J1348$-$630. We simultaneously fit the time-averaged energy spectrum of the source and the fractional rms and phase-lag spectra of the QPO with the time-dependent Comptonization model vKompth. The data can be explained by two physically connected coronae interacting with the accretion disc via a feedback loop of X-ray photons. The best-fitting model consists of a corona of $\sim$10$^3$ km located at the inner edge of the disc and a second corona of $\sim$10$^4$ km horizontally extended and covering the inner parts of the accretion disc. The properties of both coronae during the reflare are similar to those during the low/hard state of the main outburst, reinforcing the idea that both the outburst and the reflare are driven by the same physical mechanisms. We combine our results for the type-C QPO with those from previous work focused on the study of type-A and type-B QPOs with the same model to study the evolution of the geometry of the corona through the whole outburst, including the reflare of MAXI J1348$-$630. Finally, we show that the sudden increase in the phase-lag frequency spectrum and the sharp drop in the coherence function previously observed in MAXI J1348$-$630 are due to the type-C QPO during the decay of the outburst and can be explained in terms of the geometry of the coronae.

Large-scale structure (LSS) studies in cosmology map and analyse matter in the Universe on the largest scales. Understanding the LSS can provide observational support for the Cosmological Principle (CP) and the Standard Cosmological Model ($\Lambda$CDM). In recent years, many discoveries have been made of LSSs that are so large that they become difficult to understand within $\Lambda$CDM. Reasons for this are: they potentially challenge the CP, (i.e. the scale of homogeneity); and their formation and origin are not fully understood. In this article we review two recent LSS discoveries: the Giant Arc (GA, $\sim 1$ Gpc) and the Big Ring (BR, $\sim 400$ Mpc). Both structures are in the same cosmological neighbourhood -- at the same redshift $z \sim 0.8$ and with a separation on the sky of only $\sim 12^\circ$. Both structures exceed the often-cited scale of homogeneity (Yadav+ 2010), so individually and together, these two intriguing structures raise more questions for the validity of the CP and potentially hint at new physics beyond the Standard Model. The GA and BR were discovered using a novel method of mapping faint matter at intermediate redshifts, interpreted from the MgII absorption doublets seen in the spectra of background quasars.

Recent constraints on neutron star mass and radius have advanced our understanding of the equation of state of cold dense matter. Some of them have been obtained by modeling the pulses of three millisecond X-ray pulsars observed by the Neutron Star Interior Composition Explorer (NICER). Here, we present a Bayesian parameter inference for a fourth pulsar, PSR J1231-1411, using the same technique with NICER and XMM-Newton data. When applying a broad mass-inclination prior from radio timing measurements and the emission region geometry model that can best explain the data, we find likely-converged results only when using a limited radius prior. If limiting the radius to be consistent with the previous observational constraints and equation of state analyses, we infer the radius to be $12.6 \pm 0.3$ km and the mass to be $1.04_{-0.03}^{+0.05}$ $M_\odot$, each reported as the posterior credible interval bounded by the $16\,\%$ and $84\,\%$ quantiles. If using an uninformative prior but limited between $10$ and $14$ km, we find otherwise similar results, but $R_{\mathrm{eq}} = 13.5_{-0.5}^{+0.3}$ km for the radius. In both cases, we find a non-antipodal hot region geometry where one emitting spot is at the equator or slightly above, surrounded by a large colder region, and where a non-circular hot region lies close to southern rotational pole. If using a wider radius prior, we only find solutions that fit the data significantly worse. We discuss the challenges in finding the better fitting solutions, possibly related to the weak interpulse feature in the pulse profile.

We search for a stochastic gravitational-wave background (SGWB) originating from scalar-induced gravitational waves (SIGWs) with the sound speed resonance (SSR) effect using data from Advanced LIGO and Advanced Virgo's first three observing runs. The SSR mechanism, characterized by an oscillating sound speed squared term, can induce a nonperturbative parametric amplification of specific perturbation modes during inflation, leading to enhanced primordial curvature perturbations and a significant SIGW signal. We perform a Bayesian analysis to constrain the model parameters describing the SGWB spectrum from the SSR effect. Our results show no statistically significant evidence for the presence of such a signal in the current data. Consequently, we place an upper limit of $|\tau_0| \lesssim 5.9 \times 10^3$ s at $95\%$ confidence level on the start time of the oscillation in the SSR model. These results demonstrate the capability of current gravitational wave detectors to probe inflation models through the SSR mechanism and paves the way for future searches with improved sensitivity.

The Linear Mixing Approximation (LMA) is often used in planetary models for calculating the equations of state (EoSs) of mixtures. A commonly assumed planetary composition is a mixture of rock and water. Here we assess the accuracy of the LMA for pressure-temperature conditions relevant to the interiors of Uranus and Neptune. We perform MD simulations using ab-initio simulations and consider pure-water, pure-silica, and 1:1 and 1:4 silica-water molecular fractions at temperature of 3000 K and pressures between 30 and 600 GPa. We find that the LMA is valid within a few percent (<~5%) between ~150-600 Gpa, where the sign of the difference in inferred density depends on the specific composition of the mixture. We also show that the presence of rocks delays the transition to superionic water by ~70 GPa for the 1:4 silica-water mixture. Finally, we note that the choice of electronic theory (functionals) affect the EoS and introduces an uncertainty in of the order of 10% in density. Our study demonstrates the complexity of phase diagrams in planetary conditions and the need for a better understanding of rock-water mixtures and their effect on the inferred planetary composition.

A key result of the string phenomenology framework is the two swampland conjectures (SCs), namely that the values of the scalar field excursion and gradient of the potential have to be $\mathcal{O}(1)$ in Planck units. Using the first data release from the Dark Energy Spectroscopic Instrument (DESI) survey and model-independent reconstructions of the SCs at late times, via a particular machine learning approach known as the genetic algorithms (GA), we find that the reconstructions of the SCs are, in fact, in agreement with their theoretical values of $\mathcal{O}(1)$. Specifically, the reconstruction of the second SC is in fact several sigmas away from zero, thus indicating a very steep potential in contrast to recent model-specific analyses, assuming exponential potentials. Our approach is based solely on model-independent reconstructions of cosmological observables, such as the angular diameter distance and the Hubble expansion history, thus it is readily applicable to the forthcoming data from Stage IV surveys and will help settle in the near future the issue on the observational falsification, or not, of the two conjectures.

The intrinsic alignment of galaxies is a key factor in modeling weak-lensing observations and can serve as a valuable signal for both cosmological and astrophysical studies. Modelling this signal requires understanding how galaxy shapes form, and their relations to the large-scale gravitational field -- typically encoded in the value of large-scale shape-bias parameters. In this article we contribute to this topic in three ways: (i) developing new estimators of Lagrangian shape-biases (ii) applying them to measure the shape-biases of dark-matter halos (iii) interpreting these measurements to gain insight on the process of halo-shape formation. We show that our estimators produce results consistent with previous literature, and that they possess advantages with respect to previous methods, namely that the measurement of each bias parameter is completely independent from the others, and that bias parameters can be defined for each individual object. We measure universal relations between shape-bias parameters and peak-significance, $\nu$. This relation for the first-order shape-bias parameter is linear at high $\nu$, and converges to zero at low $\nu$, which we interpret as strong evidence against the proposed scenario according to which galaxy shapes arise due to post-formation interaction with the large-scale tidal-field. We anticipate our estimators to be very useful in analyzing hydrodynamical simulations to extract physical understandings of galaxy shape formation, as well as establishing priors on the values of intrinsic-alignment biases.

Context. Quasar variability can potentially unlock crucial insights into the accretion process. Understanding how this variability is influenced by wavelength is crucial for validating and refining quasar variability models. Aims. This paper aims to enhance the understanding of the dependence of variability on rest-frame wavelength by isolating the variance on different time scales in well-defined wavelength bins and examining the Corona-heated Accretion-disk (CHAR) model. Methods. We investigated the relation between variance and rest-frame wavelength using optical g and r-band light curves from the Zwicky Transient Facility (ZTF) Data Release 15 for 5000 quasars within narrow ranges of black hole mass ($M_{BH}$) and Eddington ratio ($R_{Edd}$). A spectral model taking into account disk continuum emission, Balmer transitions and Fe II pseudo-continuum emission, and other emission lines is necessary to best interpret the variance spectrum. Results. Our analysis indicates a strong anti-correlation between median variance and rest-frame wavelength for quasars with $M_{BH} = 10^{8}$ and $R_{Edd} = 10^{-1}$ at different timescales. This anti-correlation is more pronounced at shorter timescales. The results align well with a bending power-law Power Spectrum Density (PSD) model with both the damping timescale and the high-frequency slope of the PSD depending on wavelength. The accurate predictions provided by the CHAR model on the variance spectrum across most timescales studied showcases its potential in constraining temperature variations within the accretion disk. Key words. accretion, accretion discs; galaxies:active; quasars: supermassive black holes

We select and investigate six global solar extreme ultraviolet (EUV) wave events using data from the Solar Dynamics Observatory (SDO) and the Solar and Heliospheric Observatory (SOHO). These eruptions are all on the limb but recorded as halo coronal mass ejections (CMEs) because the CME-driven shocks have expanded laterally to the opposite side. With the limb observations avoiding the projection effect, we have measured the inclination and speed of the EUV wavefront from 1.05 to 1.25 $R_\odot$. We also investigate the coupling and connection of the EUV wavefront with the CME boundary and the CME-driven shock, respectively. The major findings in the six events are: (1) the forward inclination of the primary and coronal-hole transmitted EUV wavefronts is estimated, respectively, and the origins of these inclinations and their effects on the estimate of actual wavefront speed are investigated; (2) the wavefront speed can be elevated by loop systems near the coronal base, and the average speed in the low corona has no clear correlation with the lateral expansion of the CME-driven shock in the high corona; (3) the fast magnetosonic Mach number of the wavefront is larger than unity from the coronal base; (4) the EUV wavefront is coupled with the CME driver throughout the propagation in two events; (5) after the EUV wavefront vanishes, the CME-driven shock continues traveling on the opposite side and disconnects from the EUV wavefront in four events. These results and their implications are discussed, which provide insight into the properties of global EUV waves.

We present a refined speckle-interferometric orbit of a binary system $\mu$ Cet, with the main component studied based on the analysis of photometric and spectroscopic data, obtained at the SAO RAS 6-m telescope. The object was initially classified as a giant with chemical composition anomalies. As a result of our analysis, we conclude that the star belongs to the Main Sequence, to the class of non-peculiar stars. Analysis of photometric data from the TESS mission indicates that the main component of the system belongs to the $\gamma$ Dor pulsators.

HERMES (High Energy Rapid Modular Ensemble of Satellites) Pathfinder mission aims to observe and localize Gamma Ray Bursts (GRBs) and other transients using a constellation of nanosatellites in low-Earth orbit (LEO). Scheduled for launch in early 2025, the 3U CubeSats will host miniaturized instruments featuring a hybrid Silicon Drift Detector (SDD) and GAGG:Ce scintillator photodetector system, sensitive to X-rays and gamma-rays across a wide energy range. Each HERMES payload contains 120 SDD cells, each with a sensitive area of 45 mm^2, organized into 12 matrices, reading out 60 12.1x6.94x15.0 mm^3 GAGG:Ce scintillators. Photons interacting with an SDD are identified as X-ray events (2-60 keV), while photons in the 20-2000 keV range absorbed by the crystals produce scintillation light, which is read by two SDDs, allowing event discrimination. The detector system, including front-end and back-end electronics, a power supply unit, a chip-scale atomic clock, and a payload data handling unit, fits within a 10x10x10 cm^3 volume, weighs 1.5 kg, and has a maximum power consumption of about 2 W. This paper outlines the development of the HERMES constellation, the design and selection of the payload detectors, and laboratory testing, presenting the results of detector calibrations and environmental tests to provide a comprehensive status update of the mission.

Quantifying the energy content of accelerated electron beams during solar eruptive events is a key outstanding objective that must be constrained to refine particle acceleration models and understand the electron component of space weather. Previous estimations have used in situ measurements near the Earth, and consequently suffer from electron beam propagation effects. In this study, we deduce properties of a rapid sequence of escaping electron beams that were accelerated during a solar flare on 22 May 2013 and produced type III radio bursts, including the first estimate of energy density from remote sensing observations. We use extreme-ultraviolet observations to infer the magnetic structure of the source active region NOAA 11745, and Nançay Radioheliograph imaging spectroscopy to estimate the speed and origin of the escaping electron beams. Using the observationally deduced electron beam properties from the type III bursts and co-temporal hard X-rays, we simulate electron beam properties to estimate the electron number density and energy in the acceleration region. We find an electron density (above $30\ \mathrm{keV}$) in the acceleration region of $10^{2.5}\ \mathrm{cm}^{-3}$ and an energy density of $2\times10^{-5}\ \mathrm{erg\ cm}^{-3}$. Radio observations suggest the particles travelled a very short distance before they began to produce radio emission, implying a radially narrow acceleration region. A short but plausibly wide slab-like acceleration volume of $10^{26}-10^{28}\ \mathrm{cm}^{3}$ atop the flaring loop arcade could contain a total energy of $10^{23}-10^{25}\ \mathrm{erg}$ ($\sim 100$ beams), which is comparable to energy estimates from previous studies.

Solar pores are ideal magnetic structures for wave propagation and transport of energy radially-outwards across the upper layers of the solar atmosphere. We aim to model the excitation and propagation of magnetohydrodynamic waves in a pore with a lightbridge modelled as two interacting magnetic flux tubes separated by a thin, weaker field, layer. We solve the three-dimensional MHD equations numerically and calculate the circulation as a measure of net torsional motion. We find that the interaction between flux tubes results in the natural excitation of propagating torsional Alfvén waves, but find no torsional waves in the model with a single flux tube. The torsional Alfvén waves propagate with wave speeds matching the local Alfvén speed where wave amplitude peaks.

The distribution of period ratios for 580 known two-planet systems is apparently nonuniform, with several sharp peaks and troughs. In particular, the vicinity of the 2:1 commensurability seems to have a deficit of systems. Using Monte Carlo simulations and an empirically inferred population distribution of period ratios, we prove that this apparent dearth of near-resonant systems is not statistically significant. The excess of systems with period ratios in the wider vicinity of the 2:1 resonance is significant, however. Long-term WHFast integrations of a synthetic two-planet system on a grid period ratios from 1.87 through 2.12 reveal that the eccentricity and inclination exchange mechanism between non-resonant planets represents the orbital evolution very well in all cases, except at the exact 2:1 mean motion resonance. This resonance destroys the orderly exchange of eccentricity, while the exchange of inclination still takes place. Additional simulations of the Kepler-113 system on a grid of initial inclinations show that the secular periods of eccentricity and inclination variations are well fitted by a simple hyperbolic cosine function of the initial mutual inclination. We further investigate the six known two-planet systems with period ratios within 2\% of the exact 2:1 resonance (TOI-216, KIC 5437945, Kepler-384, HD 82943, HD 73526, HD 155358) on a grid of initial inclinations and for two different initial periastron longitudes corresponding to the aligned and anti-aligned states. All these systems are found to be long-term stable except HD 73526, which is likely a false positive. The periodic orbital momentum exchange is still at work in some of these systems, albeit with much shorter cycling periods of a few years.

For years, the cosmological constant $\Lambda$ and cold dark matter (CDM) model ($\Lambda\text{CDM}$) has stood as a cornerstone in modern cosmology and serves as the predominant theoretical framework for current and forthcoming surveys. However, the latest results shown by the Dark Energy Spectroscopic Instrument (DESI), along other cosmological data, show hints in favor of an evolving dark energy. Given the elusive nature of dark energy and the imperative to circumvent model bias, we introduce a novel null test, derived from Noether's theorem, that uses measurements of the Weyl potential (the sum of the spatial and temporal distortion) at different redshifts. In order to assess the consistency of the concordance model we quantify the precision of this null test through the reconstruction of mock catalogs based on $\Lambda\text{CDM}$ using forthcoming survey data, employing Genetic Algorithms, a machine learning technique. Our analysis indicates that with forthcoming LSST-like and DESI-like mock data our consistency test will be able to rule out several cosmological models at around 4$\sigma$ and help to check for tensions in the data.

This paper presents chemical abundances of twelve elements (C, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe) for 80 FGK dwarfs in the Pleiades open cluster, which span a temperature range of $\sim$2000 K in T$_{\rm eff}$, using the high-resolution (R$\sim$22,500) near-infrared SDSS-IV/APOGEE spectra ($\lambda$1.51--1.69 \micron). Using a 1D LTE abundance analysis, we determine an overall metallicity of [Fe/H]=+0.03$\pm$0.04 dex, with the elemental ratios [$\alpha$/Fe]=+0.01$\pm$0.05, [odd-z/Fe]=-0.04$\pm$0.08, and [iron peak/Fe]=-0.02$\pm$0.08. These abundances for the Pleiades are in line with the abundances of other open clusters at similar galactocentric distances as presented in the literature. Examination of the abundances derived from each individual spectral line revealed that several of the stronger lines displayed trends of decreasing abundance with decreasing $T_{\rm eff}$. The list of spectral lines that yield abundances that are independent of $T_{\rm eff}$ are presented and used for deriving the final abundances. An investigation into possible causes of the temperature-dependent abundances derived from the stronger lines suggests that the radiative codes and the APOGEE line list we employ may inadequately model van der Waals broadening, in particular in the cooler K dwarfs.

Wolf-Rayet stars (WRs) are thought to be the final evolutionary stages of massive stars and immediate progenitors of stellar-mass black holes. Their multiplicity forms an important anchor point in single and binary population models for predicting gravitational-wave progenitors. Recent spectroscopic campaigns have suggested incompatible multiplicity fractions and period distributions for N- and C-rich Galactic WRs (WNs and WCs) at short as well as long orbital periods, in contradiction with evolutionary model predictions. In this work, we employed infrared interferometry using the $K$-band instrument GRAVITY at the VLTI to investigate the multiplicity of WRs at long periods and explore the nature of their companions. We present a survey of 39 Galactic WRs, including 11 WN, 15 WC and 13 H-rich WN (WNh) stars. We detected wide companions with GRAVITY for only four stars: WR 48, WR 89, WR 93 and WR 115. Combining with spectroscopic studies, we arrived at multiplicity fractions of $f^{\rm WN}_{\rm obs} = 0.55\pm0.15$, $f^{\rm WC}_{\rm obs} = 0.40\pm0.13$ and $f^{\rm WNh}_{\rm obs} = 0.23\pm0.12$. In addition, we also found other features in the GRAVITY dataset such as (i) a diffuse extended component in over half the WR sample; (ii) five known spectroscopic binaries resolved in differential phase data and (iii) spatially resolved winds in four stars: WR 16, WR 31a, WR 78 and WR 110. Our survey reveals a lack of intermediate (few 100s d) and long- (few years to decades) period WR systems. The 200-d peak in the period distributions of WR+OB and BH+OB binaries predicted by Case B mass-transfer binary evolution models is not seen in our data. The rich companionship of their O-type progenitors in this separation range suggest that the WR progenitor stars expand and interact with their companions, most likely through unstable mass-transfer, resulting in either a short-period system or a merger.

We present a study of a previously identified candidate ``jellyfish'' galaxy in the Abell 2744 cluster, F0083, which showed faint signs of a tidal interaction in archival imaging. We used publicly available, PSF-matched deep photometric data from the Hubble and James Webb Space Telescopes, to infer the spatially-resolved star-formation history of this galaxy. F0083 shows clear signs of ram-pressure stripping (RPS), with a recently enhanced star-formation rate (SFR) orientated towards the South-West quadrant of the stellar disc. The stellar mass surface density is heavily asymmetric, with a variation of nearly 1 dex between the Western spiral arm and the postulated tidal feature. This feature appears to contain a high proportion of older stars, ruling out RPS as the cause of this ``unwinding''. An investigation of nearby galaxies reveals the most probable interaction candidate to be a compact elliptical galaxy, 28171, of mass $\log_{10}(M_*/\rm{M}_{\odot})=8.53\pm0.04$ embedded in the tidal stream along our line-of-sight. The star-formation history of the tidal feature indicates a steep change in SFR at lookback times $t_L\lesssim1\,$Gyr. The most probable formation scenario of F0083 thus indicates a low-speed interaction with 28171, followed by RPS as the combined system approaches pericentre passage. Our results demonstrate the ability of photometric data to distinguish between these consecutive processes.

Knowledge of the global magnetic field distribution and its evolution on the Sun's surface is crucial for modeling the coronal magnetic field, understanding solar wind dynamics, computing the heliospheric open flux distribution and predicting solar cycle strength. As the far side of the Sun cannot be observed directly and high-latitude observations always suffer from projection effects, we often rely on surface flux transport simulations (SFT) to model long-term global magnetic field distribution. Meridional circulation, the large-scale north-south component of the surface flow profile, is one of the key components of the SFT simulation that requires further constraints near high latitudes. Prediction of the photospheric magnetic field distribution requires knowledge of the flow profile in the future, which demands reconstruction of that same flow at the current time so that it can be estimated at a later time. By performing Observing System Simulation Experiments, we demonstrate how the Ensemble Kalman Filter technique, when used with a SFT model, can be utilized to make ``posterior'' estimates of flow profiles into the future that can be used to drive the model forward to forecast photospheric magnetic field distribution.

The star Beta Pic is widely known for harbouring a large population of exocomets, which create variable absorption signatures in the stellar spectrum as they transit the star. While the physical and chemical properties of these objects have long remained elusive, Vrignaud et al. (2024) recently introduced the exocomet curve of growth approach, enabling, for the first time, the estimate of exocometary column densities and excitation temperature using absorption measurements in many spectral lines. Using this new tool, we present a refined study of a Beta Pic exocomet observed on December 6, 1997 with the HST. We first show that the comet's signature in FeII lines is well explained by the transit of two gaseous components, with different covering factors and opacities. Then, we show that the studied comet is detected in the lines of other species, such as NiII and CrII. These species are shown to experience similar physical conditions than FeII (same radial velocity profile, same excitation temperature), hinting that they are well-mixed. Finally, using almost 100 FeII lines rising from energy levels between 0 and 33000 cm-1, we derive the complete excitation diagram of Fe+ in the comet. The transiting gas is found to be populated at an excitation temperature of 8190+-160 K, very close to the stellar effective temperature (8052 K). Using a model of radiative and collisional excitation, we show that the observed excitation diagram is compatible with a radiative regime, associated with a close transit distance (< 0.43 au) and a low electronic density (< 1e7 cm-3). In this regime, the excitation of Fe+ is controlled by the stellar flux, and do not depend on the local electronic temperature or density. These results allow us to derive the Ni+/Fe+ and Cr+/Fe+ ratios in the December 6, 1997 comet, at 8.5 +- 0.8 x 10-2 and 1.04 +- 0.15 x 10-2 respectively, close to solar abundances.

In this chapter, we will discuss the use of Machine Learning methods for the identification and localization of cometary activity for Solar System objects in ground and in space-based wide-field all-sky surveys. We will begin the chapter by discussing the challenges of identifying known and unknown active, extended Solar System objects in the presence of stellar-type sources and the application of classical pre-ML identification techniques and their limitations. We will then transition to the discussion of implementing ML techniques to address the challenge of extended object identification. We will finish with prospective future methods and the application to future surveys such as the Vera C. Rubin Observatory.

Near-sun sky twilight observations allow for the detection of asteroid interior to the orbit of Venus (Aylos), the Earth (Atiras), and comets. We present the results of observations with the Palomar 48-inch telescope (P48)/Zwicky Transient Facility (ZTF) camera in 30 s r-band exposures taken during evening astronomical twilight from 2019 Sep 20 to 2022 March 7 and during morning astronomical twilight sky from 2019 Sep 21 to 2022 Sep 29. More than 46,000 exposures were taken in evening and morning astronomical twilight within 31 to 66 degrees from the Sun with an r-band limiting magnitude between 18.1 and 20.9. The twilight pointings show a slight seasonal dependence in limiting magnitude and ability to point closer towards the Sun, with limiting magnitude slightly improving during summer. In total, the one Aylo, (594913) 'Ayló'chaxnim, and 4 Atiras, 2020 OV1, 2021 BS1, 2021 PB2, and 2021 VR3, were discovered in evening and morning twilight observations. Additional twilight survey discoveries also include 6 long-period comets: C/2020 T2, C/2020 V2, C/2021 D2, C/2021 E3, C/2022 E3, and C/2022 P3, and two short-period comets: P/2021 N1 and P/2022 P2 using deep learning comet detection pipelines. The P48/ZTF twilight survey also recovered 11 known Atiras, one Aylo, three short-period comes, two long-period comets, and one interstellar object. Lastly, the Vera Rubin Observatory will conduct a twilight survey starting in its first year of operations and will cover the sky within 45 degrees of the Sun. Twilight surveys such as those by ZTF and future surveys will provide opportunities for discovering asteroids inside the orbits of Earth and Venus.

Galaxy interactions can trigger drastic changes in the resolved star-forming properties of their constituents, but it remains unclear as to whether those changes are discernible from secular starburst triggers. In this Letter we investigate whether or not post-merger galaxies create unique star-forming trends on a kiloparsec scale. We present radial trends in star-formation-rate (SFR) surface density ($\Sigma_{SFR}$) for 150 post-merger galaxies with moderate to extremely heightened global SFRs using observations from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We juxtapose these profiles with those of noninteracting galaxies (excluding both galaxy pairs and post-merger galaxies) with similarly enhanced global SFRs. Post-merger galaxies have a much stronger central starburst than isolated galaxies with similar global star-formation enhancements. Indeed, isolated starburst galaxies (SBs) lack a marked central enhancement and instead show a fairly uniform enhancement in $\Sigma_{SFR}$ with radius. Moreover, the difference in central star formation between post-merger galaxies and noninteracting galaxies is more radially extended and pronounced when the global enhancement in star formation is larger. We conclude that post-merger galaxies create a unique signature in their resolved star-forming properties that is distinct from secular processes that can trigger similar global SFR enhancements.

The nuclear stellar disc (NSD) is a flat and dense stellar structure at the centre of the Milky Way. Previous work has identified the presence of metal-rich and metal-poor stars in the NSD, suggesting that they have different origins. The recent publication of photometric, metallicity, proper motion, and orbital catalogues allows the NSD stellar population to be characterised with unprecedented detail. We aim to explore the proper motions and orbits of NSD stars with different metallicities to assess whether they have different origins and to better understand the metallicity distribution in the NSD. We distinguished between metal-rich and metal-poor stars by applying a Gaussian mixture model, as done in previous work, and analysed the proper motions, orbits, and spatial distribution of stars with different metallicities. We find that metal-rich stars exhibit a lower velocity dispersion, suggesting that they trace a kinematically cooler component compared to metal-poor ones. Furthermore, z-tube orbits are predominant among metal-rich stars, while chaotic/box orbits are more common among metal-poor ones. We also find that metal-rich and metal-poor stars show a similar extinction and are present throughout the analysed regions. As a secondary result, we detected a metallicity gradient in the metal-rich population with higher metallicity towards the centre of the NSD and a tentative gradient for the metal-poor stars, which is consistent with previous studies that did not distinguish between the two populations. Our results suggest that metal-rich stars trace the NSD, whereas metal-poor ones are related to the Galactic bar and probably constitute Galactic bar interlopers and/or are NSD stars that originated from accreted clusters. The detected metallicity gradients aligns with the currently accepted inside-out formation of the NSD.

We investigate the impact of primordial mass segregation on the formation and evolution of dark star clusters (DSCs). Considering a wide range of initial conditions, we conducted $N$-body simulations of globular clusters (GCs) around the Milky Way. In particular, we assume a canonical IMF for all GCs without natal kicks for supernovae remnants, namely neutron stars or black holes. Our results demonstrate that clusters with larger degrees of primordial mass segregation reach their DSC phase earlier and spend a larger fraction of their dissolution time in such a phase, compared to clusters without mass segregation. In primordially segregated clusters, the maximum Galactocentric distance that the clusters can have to enter the DSC phase is almost twice that of the clusters without primordial mass segregation. Primordially segregated clusters evolve with a higher number of stellar mass black holes, accelerating energy creation in their central regions and consequently increasing evaporation rates and cluster sizes during dark phases. The simulations reveal that aggregating heavy components at the centre doubles the time spent in the dark phase. Additionally, the study identifies potential links between simulated dark clusters and initial conditions of Milky Way globular clusters, suggesting some may transition to dark phases before dissolution. Higher primordial mass segregation coefficients amplify the average binary black hole formation rate by 2.5 times, raising higher expectations for gravitational wave emissions.