Search

The first few hours of a supernova contain significant information about the progenitor system. The most modern wide-field surveys that scan the sky repeatedly every few days can discover all kinds of transients in those early epochs. At such times, some progenitor footprints may be visible, elucidating critical explosion parameters and helping to distinguish between leading explosion models. A dedicated spectroscopic classification programme using the optical spectrograph OSIRIS mounted to the Gran Telescopio de Canarias was set up to try to obtain observations of supernova at those early epochs. With the time awarded, we obtained spectra for 10 SN candidates, which we present here. Half of them were thermonuclear SNe, while the other half were core-collapse SNe. Most (70\%) were observed within the first six days of the estimated explosion, with two being captured within the first 48 hours. We present a characterization of the spectra, together with other public ancillary photometry from ZTF and ATLAS. This programme shows the need for an accompanying rapid-response spectroscopic programme to existing and future deep photometric wide-field surveys located at the right longitude to be able to trigger observations in a few hours after the discovery of the supernova candidate. Both the future La Silla Southern Supernova Survey (LS4) and the Legacy Survey of Space and Time (LSST) both located in Chile will be providing discovery and follow up of most of the transients in the southern hemisphere. This paper demonstrates that with a rapid spectroscopic programme and stringent triggering criteria, obtaining a sample of SN with spectra within a day of the explosion is possible., Comment: 20 pages, 16 figures. Submitted to A&A

Massive star formation involves significant ionization in the innermost regions near the central object, such as gravitationally trapped H II regions, jets, ionized disks, or winds. Resolved observations of the associated continuum and recombination line emission are crucial for guiding theory. The next-generation Very Large Array (ngVLA) will enable unprecedented observations of thermal emission with 1 mas resolution, providing a new perspective on massive star formation at scales down to a few astronomical units at kiloparsec distances. This work presents synthetic interferometric ngVLA observations of the free-free continuum (93-GHz band), $\mathrm{H41\alpha}$, and $\mathrm{H38\alpha}$ recombination lines from ionized jets and disks around massive protostars. Using the sf3dmodels Python package, we generate gas distributions based on analytical models, which are then processed through the RADMC-3D radiative transfer code. Our results indicate that the ngVLA can easily resolve, both spatially and spectrally, the ionized jet from a 15 $\mathrm{M_\odot}$ protostar at 700 pc, distinguishing between collimated jets and wide-angle winds, and resolving their launching radii, widths, and any substructure down to a few astronomical units. Detailed studies of radio jets launched by massive protostars will be feasible up to distances of $\sim 2$ kpc. Furthermore, ngVLA will be able to study in detail the ionized disks around massive ($> 10~\mathrm{M_\odot}$) protostars up to distances from 4 to 12 kpc, resolving their kinematics and enabling the measurement of their central masses across the Galaxy. These observations can be conducted with on-source integrations of only a few hours., Comment: Accepted for publication in ApJ

Network simulation is pivotal in network modeling, assisting with tasks ranging from capacity planning to performance estimation. Traditional approaches such as Discrete Event Simulation (DES) face limitations in terms of computational cost and accuracy. This paper introduces RouteNet-Gauss, a novel integration of a testbed network with a Machine Learning (ML) model to address these challenges. By using the testbed as a hardware accelerator, RouteNet-Gauss generates training datasets rapidly and simulates network scenarios with high fidelity to real-world conditions. Experimental results show that RouteNet-Gauss significantly reduces prediction errors by up to 95% and achieves a 488x speedup in inference time compared to state-of-the-art DES-based methods. RouteNet-Gauss's modular architecture is dynamically constructed based on the specific characteristics of the network scenario, such as topology and routing. This enables it to understand and generalize to different network configurations beyond those seen during training, including networks up to 10x larger. Additionally, it supports Temporal Aggregated Performance Estimation (TAPE), providing configurable temporal granularity and maintaining high accuracy in flow performance metrics. This approach shows promise in improving both simulation efficiency and accuracy, offering a valuable tool for network operators., Comment: 13 pages, 11 figures

We present the first characterization of the statistical relationship between a large sample of novae in M31 and their progenitor stellar populations in the form of a delay time distribution. To this end, we leverage the spatially resolved stellar age distribution of the M31 disk derived from deep HST photometry by the Panchromatic Hubble Andromeda Treasury (PHAT) survey and a large catalog of novae in M31. Our delay time distribution has two statistically significant detections: one population of nova progenitors, ages between 2 and 3.2 Gyr, with an unnormalized rate of ($3.7^{+6.8}_{-3.5} \pm 2.1) \cdot 10^{-9}$ events / $M_{\odot}$, and another of ages between 7.9 Gyr and the age of the Universe with ($4.8^{+1.0}_{-0.9} \pm 0.2) \cdot 10^{-9}$ events / $M_{\odot}$ (uncertainties are statistical and systematic, respectively). Together with the upper limits we derive at other time bins, these detections are consistent with either a constant production efficiency or a higher production efficiency of novae at earlier delay times., Comment: 14 pages, 6 figures. Submitted to The Astrophysical Journal, currently under review

We develop a method for constraining the potential of the Milky Way using stellar streams with a known progenitor. The method expresses the stream in angle-action coordinates and integrates the orbits of the stars backwards in time to obtain the stripping point positions of the stream stars relative to the cluster. In the potential that generated the stream, the stars return approximately to the cluster centre. In a different potential, they are redirected to different locations. The free parameters of the model are estimated by maximising the degree of clustering of the stripping point distribution. We test this method with the stellar stream of the globular cluster M68 (NGC 4590). We use an N-body code to simulate the stream and generate a realistic star sample using a model of the Gaia selection function. We also simulate the expected observational uncertainties, and estimate the heliocentric distances and radial velocities of the stream stars from the cluster orbit. Using this sample of stars, we recover the values of four free parameters characterising the potential of the disc and the dark halo to an accuracy of about 10 per cent of the correct values. We show that this accuracy is improved up to about 2 per cent using the expected end-of-mission Gaia data. In addition, we obtain a strong correlation between the mass of the disc and the dark halo axis ratio, which is explained by the fact that the stream flows close and parallel to the disc plane., Comment: 27 pages, 16 figures

Brown dwarfs are the bridge between low-mass stars and giant planets. One way of shedding light on their dominant formation mechanism is to study them at the earliest stages of their evolution, when they are deeply embedded in their parental clouds. Several works have identified pre- and proto-brown dwarfs candidates using different observational approaches. The aim of this work is to create a database with all the objects classified as very young substellar candidates in the litearature in order to study them in an homogeneous way. We have gathered all the information about very young substellar candidates available in the literature until 2020. We have retrieved their published photometry from the optical to the centimeter regime, and we have written our own codes to derive their bolometric temperatures and luminosities, and their internal luminosities. We have also populated the database with other parameters extracted from the literature, like e.g. the envelope masses, their detection in some molecular species, and presence of outflows. The result of our search is the SUCANES database, containing 174 objects classified as potential very young substellar candidates in the literature. We present an analysis of the main properties of the retrieved objects. Since we have updated the distances to several star forming regions, this has allowed us to reject some candidates based on their internal luminosities. We have also discussed the derived physical parameters and envelope masses for the best substellar candidates isolated in SUCANES. As an example of a scientific exploitation of this database, we present a feasibility study for the detection of radiojets with upcoming facilities: the ngVLA and the SKA interferometers. The SUCANES database is accessible through a Graphical User Interface and it is open to any potential user., Comment: 17 pages, 7 figures

The landscape of computing technologies is changing rapidly, straining existing software engineering practices and tools. The growing need to produce and maintain increasingly complex multi-architecture applications makes it crucial to effectively accelerate and automate software engineering processes. At the same time, artificial intelligence (AI) tools are expected to work hand-in-hand with human developers. Therefore, it becomes critical to model the software accurately, so that the AI and humans can share a common understanding of the problem. In this contribution, firstly, an in-depth overview of these interconnected challenges faced by modern software engineering is presented. Secondly, to tackle them, a novel architecture based on the emerging WebAssembly technology and the latest advancements in neuro-symbolic AI, autonomy, and knowledge graphs is proposed. The presented system architecture is based on the concept of dynamic, knowledge graph-based WebAssembly Twins, which model the software throughout all stages of its lifecycle. The resulting systems are to possess advanced autonomous capabilities, with full transparency and controllability by the end user. The concept takes a leap beyond the current software engineering approaches, addressing some of the most urgent issues in the field. Finally, the efforts towards realizing the proposed approach as well as future research directions are summarized.

Brown Dwarfs (BDs) are crucial objects in our understanding of both star and planet formation. However, there is still an unconcluded debate about which is the dominant formation mechanism of these objects. For this, it is mandatory to study BDs in their earliest evolutionary stages (what we call pre- and proto-BDs), comparable to the `pre-stellar' and `Class 0/I' stages well characterized for the formation of low-mass stars. In this review, the recent efforts aimed at searching, identifying and characterising pre- and proto-BD candidates in nearby star-forming regions are presented, and revised requirements for an object to be a promising proto-BD or pre-BD candidate are provided, based on a new, unexplored so far, relation between the internal luminosity and the accreted mass. By applying these requirements, a list of 67 promising proto-BD candidates is presented, along with a compilation of possible pre-BDs from the literature. Updated correlations of protostellar properties such as mass infall rate or outflow momentum rate with bolometric luminosity are provided down to the low-mass BD regime, where no significant deviations are apparent. Furthermore, the number of proto-BD candidates in different clouds of the Solar Neighborhood seem to follow the known relations of number of protostars with cloud properties. In addition, proto(star-to-BD) ratios for the different clouds are also explored, unveiling a particular underproduction of low-mass proto-BD candidates in Ophiuchus compared to Lupus and Taurus. Possible explanations for this behavior are discussed, including heating of the Ophiuchus cloud by the nearby OB stars. The overall results of this work tend to favor a star-like process for BD formation down to the planetary boundary, of about 10 Mjup, below which other mechanisms might be at work., Comment: Invited review to New Astronomy Reviews. 63 pages. In press

The energy gap of Bernal-stacked bilayer graphene can be tuned by applying a perpendicular electric field. The origin of this gap can be traced down to the breaking of its inversion symmetry by an onsite potential difference between the layers. This degree of tunability makes bilayer graphene a perfect playground for the study of the effects of electric fields, such as the crystalline field, which are developed when layers of other materials are deposited on top of it. Here, we introduce a novel device architecture allowing a simultaneous control over the applied displacement field and the crystalline alignment between two materials. Our experimental and numerical results confirm that the crystal field and electrostatic doping due to the interface reflect the 120$^{\circ}$ symmetry of the bilayer graphene/BN heterostructure and are highly affected by the commensurate state. These results provide an unique insight into the role of the twist angle in the development of internal crystal fields and intrinsic electrostatic doping in heterostructures. Our results highligth the importance of layer alignment, beyond the existence of a moir\'e superlattice, to understand and simulate the intrinsic properties of a van der Waal heterostructure.

Hub-filament systems are considered as natural sites for high-mass star formation. Kinematic analysis of the surroundings of hub-filaments is essential to better understand high-mass star formation within such systems. In this work, we present a detailed study of the massive Galactic protocluster IRAS 15394$-$5358, using continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey. The 3~mm dust continuum map reveals the fragmentation of the massive ($\rm M=843~M_{\odot}$) clump into six cores. The core C-1A is the largest (radius = 0.04~pc), the most massive ($\rm M=157~M_{\odot}$), and lies within the dense central region, along with two smaller cores ($\rm M=7~and~3~M_{\odot}$). The fragmentation process is consistent with the thermal Jeans fragmentation mechanism and virial analysis shows that all the cores have small virial parameter values ($\rm \alpha_{vir}<<2$), suggesting that the cores are gravitationally bound. The mass vs. radius relation indicates that three cores can potentially form at least a single massive star. The integrated intensity map of $\rm H^{13}CO^{+}$ shows that the massive clump is associated with a hub-filament system, where the central hub is linked with four filaments. A sharp velocity gradient is observed towards the hub, suggesting a global collapse where the filaments are actively feeding the hub. We discuss the role of global collapse and the possible driving mechanisms for the massive star formation activity in the protocluster., Comment: 23 pages, 19 figures, accepted for publication in MNRAS

Take a self-similar fragmentation process with dislocation measure $\nu$ and index of self-similarity $\alpha > 0$. Let $e^{-m_t}$ denote the size of the largest fragment in the system at time $t\geq 0$. We prove fine results for the asymptotics of the stochastic process $(m_t)_{t \geq 0}$ for a broad class of dislocation measures satisfying $$\nu(1-s_1 > \delta ) = \delta^{-\theta} \ell(1/\delta),$$ for some $\theta \in (0,1)$ and $\ell:(0,\infty) \to (0,\infty)$ slowly varying at infinity. Under this regularity condition, we find that if for $t \geq 0$ we have the almost-sure convergence $\lim_{t \to \infty} (m_t - g(t)) = 0$, where $$ g(t) :=\left( \log t - (1-\theta) \log \log t +f(t) \right)/\alpha$$ and $f(t) = o(\log \log t)$ can be given explicitly in terms of $\alpha,\theta$ and $\ell(\cdot)$. We prove a similar result in the finite activity case, which corresponds roughly speaking to the setting $\theta = 0$. Our results sharpen significantly the best prior result on general self-similar fragmentation processes, due to Bertoin, which states that $m_t = (1+o(1)) \log (t)/\alpha$., Comment: are welcome!

We explore the Schmidt-Kennicutt (SK) relations and the star formation efficiency per free-fall time ($\eff$), mirroring observational studies, in numerical simulations of filamentary molecular clouds undergoing gravitational contraction. We find that {\it a)} collapsing clouds accurately replicate the observed SK relations for galactic clouds and {\it b)} the so-called efficiency per free-fall time ($\eff$) is small and constant in space and in time, with values similar to those found in local clouds. This constancy is a consequence of the similar radial scaling of the free-fall time and the internal mass in density structures with spherically-averaged density profiles near $r^{-2}$. We additionally show that {\it c)} the star formation rate (SFR) increases rapidly in time; {\it d)} the low values of $\eff$ are due to the different time periods over which $\tauff$ and $\tausf$ are evaluated, together with the fast increasing SFR, and {\it e)} the fact that star clusters are significantly denser than the gas clumps from which they form is a natural consequence of the fast increasing SFR, the continuous replenishment of the star-forming gas by the accretion flow, and the near $r^{-2}$ density profile generated by the collapse Finally, we argue that the interpretation of $\eff$ as an efficiency is problematic because its maximum value is not bounded by unity, and because the total gas mass in the clouds is not fixed, but rather depends on the environment where clouds are embedded. In summary, our results show that the SK relation, the typical observed values of $\eff$, and the mass density of clusters arise as a natural consequence of gravitational contraction., Comment: Submitted to MNRAS

We explain the large variety of star formation laws in terms of one single, simple law that can be inferred from the definition of the star formation rate and basic algebra. The resulting equation, $\SFR = \eff\ \Mcollapsing/\tauff$, although it has been presented elsewhere, is interpreted in terms of clouds undergoing collapse { rather than being turbulence-supported, an idea that different groups have pursued this century}. Under such assumption, one can explain the constancy of $\eff$, the different intra-cloud correlations observed in Milky Way's molecular clouds, as well as the resolved and unresolved extragalactic relationships between SFR and a measurement of the mass in CO, HCN, and CO+HI. We also explain why the slope of the correlation changes when the orbital time $\tauorb$ is considered instead of the free-fall time, and why estimations of the free-fall time from extragalactic observations skew the correlation, providing a false sublinear correlation. We furthermore show that the apparent nearly linear correlation between the star formation rate and the dynamical equilibrium pressure in the midplane of the galaxies, $\PDE$, is just a consequence of $\PDE$ values being dominated by the variation of the column density of molecular gas. All in all, we argue that the star formation law is driven by the collapse of cold, dense gas, which happens to be primarily molecular in the present Universe, and that the role of stellar feedback is just to shut down the star formation process, not to shape the star formation law., Comment: 17 pages, one figure. Accepted by MNRAS

We provide a detailed comparison between the ``gravoturbulent'' (GT) and ``global hierarchical collapse'' (GHC) models for molecular clouds and star formation, their respective interpretations of the observational data, the features they share, and suggested tests and observations to discern between them. Also, we clarify common misconceptions in recent literature about the global and hierarchical nature of the GHC scenario, and briefly discuss the evolution of some aspects of both models toward convergence. GT assumes that molecular clouds and their substructures are in approximate virial equilibrium and are in a near-stationary state, interprets the linewidth-derived nonthermal motions exclusively as turbulence, which provides additional pressure against self-gravity. Conversely, GHC assumes that most star-forming molecular clouds and their substructures are part of a continuous gravitationally-driven flow, each accreting from their parent structure. Thus, the clouds and their star formation rate evolve in time. GHC interprets nonthermal motions as a mixture of infall and turbulent components, with the relative importance of the former increasing as the objects become denser and/or more massive. Tests that may provide clues to distinguishing between GT and GHC must take into account that the innermost parts of globally gravitationally bound structures may not locally appear bound, and thus the binding may have to be searched for at the largest scale of the structure., Comment: Submitted to MNRAS. Comments welcome. Animation for Fig. 3 available at https://www.irya.unam.mx/gente/g.gomez/publica/24-GT_vs_GHC-f3.mpg

The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of Crutcher et al. (2010) with a large body of magnetic data that includes 700 observations with the Davis-Chandrasekhar-Fermi method. By using a new multiparameter Bayesian analysis we present a new, more general, time-averaged observational relation: B \propto n^{0.27 \pm 0.017} for n \leq n_0 and B \propto n^{0.54 \pm 0.18} for n \geq n_0, with n_0 = 924^(+145-144) cm^-3. We perform a separate analysis on 19 numerical magnetohydrodynamics simulations that cover a wide range of scales, resolutions, initial conditions, and completed with a variety of codes: arepo, flash, pencil, and ramses. The power law exponents derived from the simulations depend on several physical factors including: dynamo effects, time scales, turbulence, and the initial seed field strength. In particular, early-time simulations where the density, velocity and magnetic fields are unevolved do not match the observational scalings. The simulations that trace the observed density range best, the evolved dwarf galaxy and Milky Way like galaxy simulations, settle into a near consistent exponent of = 0.5 in the dense gas, with variability in the diffuse gas exponent., Comment: 23 figures, 28 pages, submitted to MNRAS, Comments welcome

Nanomechanical measurements of minimally twisted van der Waals materials remained elusive despite their fundamental importance for device realisation. Here, we use Ultrasonic Force Microscopy (UFM) to locally quantify the variation of out-of-plane Young's modulus in minimally twisted double bilayer graphene (TDBG). We reveal a softening of the Young's modulus by 7\% and 17\% along single and double domain walls, respectively. Our experimental results are confirmed by force-field relaxation models. This study highlights the strong tunability of nanomechanical properties in engineered twisted materials, and paves the way for future applications of designer 2D nanomechanical systems.

Transformation optics applied to low frequency magnetic systems has been recently implemented to design magnetic field concentrators and cloaks with superior performance. Although this achievement has been amply demonstrated theoretically and experimentally in bulk 3D macrostructures, the performance of these devices at low dimensions remains an open question. In this work, we numerically investigate the non-monotonic evolution of the gain of a magnetic metamaterial field concentrator as the axial dimension is progressively shrunk. In particular, we show that in planar structures the role played by the diamagnetic components becomes negligible, whereas the paramagnetic elements increase their magnetic field channeling efficiency. This is further demonstrated experimentally by tracking the gain of superconductor-ferromagnet concentrators through the superconducting transition. Interestingly, for thicknesses where the diamagnetic petals play an important role for the concentration gain, they also help to reduce the stray field of the concentrator, thus limiting the perturbation of the external field (invisibility). Our findings establish a roadmap and set clear geometrical limits for designing low dimensional magnetic field concentrators.

In this study we use MUSE Integral Field Spectroscopy (IFS), along with multi-line diagnostics, for the optical identification of Supernova Remnants (SNRs) in the galaxy NGC 7793. We find in total 238 SNR candidates, 225 of them new identifications, increasing significantly the number of known SNRs in this galaxy. The velocity dispersion of the candidate SNRs was calculated, giving a mean value of 27 $ km\,s^{ -1} $. We construct the H$\rm \alpha$, [S II], [O III], and [S II] - H$\rm \alpha$ luminosity functions, and for the first time, the [N II], [N II] - H$\rm \alpha$, [N II] - [S II], [O III] - [S II], and [O III] - [N II] luminosity functions of the candidate SNRs. Shock models, along with the observed multi-line information were used, in order to estimate shock velocities. The $\sim$ 65% of the SNRs present velocities < 300 $ km\,s^{ -1} $. There is a clear correlation between shock velocity and [O III]/H$\rm \beta$ ratio, and a less clear but still evident correlation in the relation between shock velocity and the [S II]/H$\rm \alpha$, [N II]/H$\rm \alpha$ ratios. We also use the [S II]6716/31 ratio of the SNR candidates to calculate their post-shock density, assuming different temperatures. The median value of the density of our sample is $\sim$ 80 $ cm^{ -3} $, for a temperature of T = $10^4\,K$. No correlation between shock velocity and density, or density and SNRs with [S II]/H$\rm \alpha$ > 0.4 and [S II]/H$\rm \alpha$ < 0.4 is observed., Comment: Accepted for publication in MNRAS. 18 pages (main paper)

Rapid progress in information technologies has spurred the need for innovative memory concepts, for which advanced data-processing methods and tailor-made materials are required. Here we introduce a previously unexplored nanoscale magnetic object: an analog magnetic vortex controlled by electric-field-induced ion motion, termed magneto-ionic vortex or "vortion". This state arises from paramagnetic FeCoN through voltage gating and gradual N3-ion extraction within patterned nanodots. Unlike traditional vortex states, vortions offer comprehensive analog adjustment of key properties such as magnetization amplitude, nucleation/annihilation fields, or coercivity using voltage as an energy-efficient tuning knob. This manipulation occurs post-synthesis, obviating the need for energy-demanding methods like laser pulses or spin-torque currents. By leveraging an overlooked aspect of N3-magneto-ionics -- planar ion migration within nanodots -- precise control of the magnetic layer's thickness is achieved, which enables reversible transitions among paramagnetic, single-domain, and vortion states, offering future prospects for analog computing, multi-state data storage, or brain-inspired devices.

A number of works infer radial temperature profiles of envelopes surrounding young stellar objects using several rotational transitions in a pixel-by-pixel or azimuthally-averaged basis. However, in many cases the assumption that the rotational temperature is constant along the line of sight is made, while this is not the case when a partially resolved envelope, assumed to be spherically symmetric, is used to obtain values of temperature for different projected radii. This kind of analysis (homogeneous analysis) is intrinsically inconsistent. By using a spherical envelope model to interpret NH3 (1, 1) and (2, 2) observations, we tested how robust it is to infer radial temperature profiles of an envelope. The temperature and density of the model envelope are power laws of radius, but the density can be flat for an inner central part. The homogeneous analysis was applied to obtain radial temperature profiles, and resulted that for small projected radii, where the optical depth of the lines is high, the homogeneous temperature can be much higher than the actual envelope temperature. In general, for larger projected radii, both the temperature and the temperature power-law index can be underestimated by as much as 40%, and 0.15, respectively. We applied this study to the infrared dark cloud G14.225-0.506 for which the radial temperature profile was previously derived from the dust emission at submillimeter wavelengths and the spectral energy distribution. As expected, the homogeneous analysis underestimated both the temperature and the temperature power-law index., Comment: Accepted for publication in MNRAS

Only a handful of massive starless core candidates have been discovered so far, but none of them have been fully confirmed. Within the MM1 clump in the filamentary infrared dark cloud G34.43+0.24 that was covered by the ALMA-ATOMS survey at Band 3 ($\sim2\arcsec$, 6000\,au) and the ALMA-QUARKS survey at Band 6 ($\sim 0.3\arcsec$, 900\,au), two prestellar core candidates MM1-C and E1 with masses of 71 and 20 \solarmass~and radii of 2100--4400\,au were discovered. The two cores show no obvious sign of star-formation activities. In particular, MM1-C is a very promising massive prestellar core candidate with a total gas mass of 71\,\solarmass. Within MM1-C, we detected two extremely dense substructures, C1 and C2, as characterized by their high densities of $\rm n_{H_2}\sim 10^{8-9} cm^{-3}$. Moreover, evidence of further fragmentation in C2 was also revealed. We have detected the primordial fragmentation in the earliest stage of massive star formation, and we speculate that MM1-C would be the birthplace of a massive multiple system. However, we cannot fully rule out the possibility that the massive prestellar core MM1-C will just form a cluster of low-mass stars if it undergoes further fragmentation., Comment: 12 pages, 6 figures

We present the Bidimensional Exploration of the warm-Temperature Ionised gaS (BETIS) project, designed for the spatial and spectral study of the diffuse ionised gas (DIG) in a selection of nearby spiral galaxies observed with the MUSE integral-field spectrograph. Our primary objective is to investigate the various ionisation mechanisms at play within the DIG. We analysed the distribution of high- and low-ionisation species in the optical spectra of the sample on a spatially resolved basis. We introduced a new methodology for spectroscopically defining the DIG, optimised for galaxies of different resolutions. We also examined the suitability of using Ha equivalent width (EWha) as a proxy for defining the DIG and its associated ionisation regime. The average radial distribution of the main line ratios are enhanced in the DIG. It follows similar trends between the DIG regime and the HII regions, as well as the Ha surface brightness, indicating a correlation between the ionisation of these species in both the DIG and the HII regions. The DIG loci in diagnostic diagrams are found within the line ratios that correspond to photoionisation due to the star formation. There is a noticeable offset correspondent to fast shocks. However, an individual diagnosis performed for each galaxy reveals that all the DIG in these galaxies can be attributed to photoionisation from star formation. The offset is primarily due to the contribution of Seyfert galaxies in our sample, which is closely aligned with models of ionisation from fast shocks and galactic outflows, thus mimicking the DIG emission. Our results indicate that galaxies exhibiting active galactic nucleus activity should be considered separately when conducting a general analysis of the DIG ionisation mechanisms, since this emission is indistinguishable from high-excitation DIG., Comment: 21 pages, 15 figures, A&A accepted

The IRDC G14.225-0.506 is associated with a network of filaments, which result in two different dense hubs, as well as with several signposts of star formation activity. The aim of this work is to study the cm continuum emission to characterize the stellar population in G14.2. We performed deep (~1.5-3 microJy) radio continuum observations at 6 and 3.6 cm using the VLA in the A configuration (~0.3''). We have also made use of observations taken during different days to study the presence of variability at short timescales. We detected a total of 66 sources, 32 in the northern region G14.2-N and 34 in the southern region G14.2-S. Ten of the sources are found to be variable. Based on their spectral index, the emission in G14.2-N is mainly dominated by non-thermal sources while G14.2-S contains more thermal emitters. Approximately 75% of the sources present a counterpart at other wavelengths. In the inner 0.4~pc region around the center of each hub, the number of IR sources in G14.2-N is larger than in G14.2-S by a factor of 4. We also studied the relation between the radio luminosity and the bolometric luminosity, finding that the thermal emission of the studied sources is compatible with thermal radio jets. For our sources with X-ray counterparts, the non-thermal emitters follow a G\"udel-Benz relation with k = 0.03. We found similar levels of fragmentation between G14.2-N and G14.2-S, suggesting that both regions are most likely twin hubs. The non-thermal emission found in the less evolved objects suggests that G14.2-N may be composed of more massive YSOs as well as being in a more advanced evolutionary stage, consistent with the filament-halo gradient in age and mass from previous works. Our results confirm a wider evolutionary sequence starting in G14.2-S as the youngest part, followed by G14.2-N, and ending with the most evolved region M17., Comment: 28 pages, 22 figures. Accepted for publication in A&A

Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions., Comment: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

The two point correlation function (2PCF) is a powerful statistical tool to measure galaxy clustering. Although 2PCF has also been used to study the clustering of stars on parsec and sub-parsec scales, its physical implication is not clear on such non-linear scales. In this study, we use the Illustris-TNG50 simulation to study the connection between the 2PCF signals of accreted halo stars and the assembly histories of Milky Way-mass galaxies. We find, in general, the 2PCF signal increases with the increase in galactocentric radii, $r$, and with the decrease in the pair separations. Galaxies which assemble late on average have stronger 2PCF signals. With $z_{1/4}$, $z_{1/2}$ and $z_{3/4}$ defined as the redshifts when galaxies accreted one-fourth, half and three-fourths of their ex-situ stellar mass at today, we find all of them show the strongest correlations with the 2PCF signals at $r\sim0.2R_{200}$. $z_{3/4}$ shows the strongest correlations at all radii than those of $z_{1/4}$ or $z_{1/2}$, as later accreted stars preserve better clusterings. However, the correlations between the 2PCF signals at different radii and the galaxy formation times all have large scatters. The 2PCFs in velocity space show weaker correlations with the galaxy formation times within $0.38R_{200}$ than real space 2PCFs, and the scatter is considerably large. Both the real and velocity space 2PCFs correlate with the assembly histories of the host dark matter halos as well. Within $0.38R_{200}$, the real space 2PCF shows stronger correlations with the galaxy formation histories than with the halo formation histories, while the velocity space 2PCFs do not show large differences. We conclude that it is difficult to use 2PCF alone to precisely predict the formation times or assembly histories of galaxies., Comment: accepted by ApJ

Currently the state of the art network models are based or depend on Discrete Event Simulation (DES). While DES is highly accurate, it is also computationally costly and cumbersome to parallelize, making it unpractical to simulate high performance networks. Additionally, simulated scenarios fail to capture all of the complexities present in real network scenarios. While there exists network models based on Machine Learning (ML) techniques to minimize these issues, these models are also trained with simulated data and hence vulnerable to the same pitfalls. Consequently, the Graph Neural Networking Challenge 2023 introduces a dataset of captured traffic traces that can be used to build a ML-based network model without these limitations. In this paper we propose a Graph Neural Network (GNN)-based solution specifically designed to better capture the complexities of real network scenarios. This is done through a novel encoding method to capture information from the sequence of captured packets, and an improved message passing algorithm to better represent the dependencies present in physical networks. We show that the proposed solution it is able to learn and generalize to unseen captured network scenarios., Comment: 8 pages, 4 figures

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au at 350 GHz (870 um) in continuum and line emission. 248 dense cores were identified, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early-stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion., Comment: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

In the last years there has been a substantial increase in the number of the reported massive and luminous star-forming regions with related explosive outflows thanks to the superb sensitivity and angular resolution provided by the new radio, infrared, and optical facilities. Here, we report one more explosive outflow related with the massive and bright star-forming region IRAS 12326$-$6245 using Band 6 sensitive and high angular resolution ($\sim$0.2$"$) Atacama Large Millimeter/Submillimeter Array (ALMA) observations. We find over 10 molecular and collimated well-defined streamers, with Hubble-Lemaitre like expansion motions, and pointing right to the center of a dusty and molecular shell (reported for the first time here) localized in the northern part of the UCHII region known as G301.1A. The estimated kinematic age, and energy for the explosion are $\sim$700 yrs, and 10$^{48}$ erg, respectively. Taking into account the recently reported explosive outflows together with IRAS 12326$-$6245, we estimate an event rate of once every 90 yr in our Galaxy, similar to the formation rate of massive stars., Comment: Accepted to The Astrophysical Journal Letters

We present optical, infrared, ultraviolet, and radio observations of SN 2022xkq, an underluminous fast-declining type Ia supernova (SN Ia) in NGC 1784 ($\mathrm{D}\approx31$ Mpc), from $<1$ to 180 days after explosion. The high-cadence observations of SN 2022xkq, a photometrically transitional and spectroscopically 91bg-like SN Ia, cover the first days and weeks following explosion which are critical to distinguishing between explosion scenarios. The early light curve of SN 2022xkq has a red early color and exhibits a flux excess which is more prominent in redder bands; this is the first time such a feature has been seen in a transitional/91bg-like SN Ia. We also present 92 optical and 19 near-infrared (NIR) spectra, beginning 0.4 days after explosion in the optical and 2.6 days after explosion in the NIR. SN 2022xkq exhibits a long-lived C I 1.0693 $\mu$m feature which persists until 5 days post-maximum. We also detect C II $\lambda$6580 in the pre-maximum optical spectra. These lines are evidence for unburnt carbon that is difficult to reconcile with the double detonation of a sub-Chandrasekhar mass white dwarf. No existing explosion model can fully explain the photometric and spectroscopic dataset of SN 2022xkq, but the considerable breadth of the observations is ideal for furthering our understanding of the processes which produce faint SNe Ia., Comment: 38 pages, 16 figures, accepted for publication in ApJ, the figure 15 input models and synthetic spectra are now available at https://zenodo.org/record/8379254