Papers for Monday, Mar 24 2025

Accurate identification of solar radio bursts (SRBs) is essential for advancing research in solar physics and predicting space weather. However, the majority of current studies mainly concentrate on detecting whether SRBs are present or absent, often focusing on only one particular type of burst. Moreover, the neural network models used for SRB detection are typically complex, involving a large number of parameters, which results in slower processing speeds. This study establishes a dataset encompassing Type II, Type III, Type IIIs, Type IV, and Type V SRBs collected from e-CALLISTO, including 8,752 SRB spectrum images and achieving annotations for 10,822 SRBs. We propose a multi-category SRB detection model based on task-aligned one-stage object detection (TOOD). TOOD can solve the problem of inconsistent predictions in classification and localization tasks, and it improves the detection recall rate. This model aligns classification and localization tasks and optimizes the neck network by incorporating a channel attention mechanism. This model achieves higher recall and accuracy with fewer parameters. This model can accurately detect five types of SBRs. The experimental results show that the model achieved an accuracy of 79.9\% (AP50) and a recall rate of 95.1\% on the SBRs dataset. A higher recall rate than other models means fewer SRBs are missed in automatic detection. The model we propose has the potential to make a substantial impact on solar physics research and space weather studies. Additionally, the findings in this paper could provide valuable insights for processing other small-sample astronomical this http URL source code and data is available at this https URL.

During the 2024's quinquennial scientific roadmap of CNES, a specific group worked on setting recommendations to decrease the environmental footprint of space science activities. This correspondence to Nature Astronomy highlights the efforts of the french space research to move towards sustainability. It relies on two complementary methods: decarbonisation and frugality.

JWST has uncovered large numbers of compact galaxies at high redshift with broad hydrogen/helium lines. These include the enigmatic population known as "little red dots" (LRDs). Their nature is debated, but they are thought to be powered by supermassive black holes (SMBHs) or intense star formation. They exhibit unusual properties for SMBHs, such as black holes that are overmassive for their host galaxies and extremely weak X-ray and radio emission. Using the highest-quality JWST spectra, we show here that the lines are broadened by electron scattering with a narrow intrinsic line core. The data require high electron column densities and compact sizes (light days), which, when coupled with their high luminosities can only be explained by SMBH accretion. The narrow intrinsic cores of the lines imply upper limits on the black hole masses of $10^{5-7}$ $M_{\odot}$, two orders of magnitude lower than previous estimates. These are among the lowest mass SMBHs known at high redshift and suggest that this is a population of young, rapidly growing SMBHs. They are enshrouded in a dense cocoon of ionized gas, probably related to their youth, from which they are accreting close to the Eddington limit. Reprocessed nebular emission from the dense cocoon dominates the optical spectrum, explaining most LRD spectral characteristics and helping to suppress radio and X-ray emission.

The abundance and distribution of sulfur-bearing molecules in protoplanetary disks directly influences the composition and potential habitability of nascent planets in addition to providing powerful probes of the physical gas conditions in the disks themselves. Here, we present new and archival ALMA and SMA observations of CS and H$_2$CS, and their C$^{34}$S and H$_2$C$^{34}$S isotopologues, at high-angular resolution (${\approx}$0."2-0."4; 20-40 au) in the HD 163296 disk, which reveal a central cavity and multi-ringed emission structure. These observations comprise the most comprehensive, multi-line CS data in a planet-forming disk to date, spanning a wide range of excitation conditions from E$_{\rm{u}}$=7.1 K to 129.3 K, and include new detections of C$^{34}$S, H$_2$CS, and H$_2$C$^{34}$S in this system. Using these data, we derive spatially-resolved rotational temperature and column density profiles for all species. We find a column density ratio N(H$_2$CS)/N(CS) $\approx$ 0.5, which is comparable to that of the similar MWC 480 disk and suggests that organic sulfur compounds may constitute a large fraction of the volatile sulfur reservoir in disks around Herbig stars generally. We derive $^{32}$S/$^{34}$S ratios of ${\approx}$5 (CS/C$^{34}$S) and ${\approx}$2 (H$_2$CS/H$_2$C$^{34}$S) based on disk-averaged and spatially-resolved analyses. Both values are consistent across these two pairs of optically-thin molecules and are well-below the expected ISM ratio of ${\approx}$22, suggesting significant sulfur fractionation. We also constrain the CS emitting layer ($z/r\lesssim 0.1$) using the vertical separations of the disk surfaces in the channel maps and based on the known 2D gas structure of the HD 163296 disk combined with our excitation analysis.

The extreme energy of the KM3-230213A event could transform our understanding of the most energetic sources in the Universe. However, it also reveals an inconsistency between the KM3NeT detection and strong IceCube constraints on the ultra-high energy neutrino flux. The most congruous explanation for the KM3NeT and IceCube data requires KM3-230213A to be produced by a (potentially transient) source fortuitously located in a region where the KM3NeT acceptance is maximized. In hadronic models of ultra-high-energy neutrino production, such a source would also produce a bright {\gamma}-ray signal, which would cascade to GeV--TeV energies due to interactions with extragalactic background light. We utilize the {\gamma}-Cascade package to model the spectrum, spatial extension, and time-delay of such a source, and scan a region surrounding the KM3NeT event to search for a consistent {\gamma}-ray signal. We find no convincing evidence for a comparable \textit{Fermi}-LAT source and place constraints on a combination of the source redshift and the intergalactic magnetic field strength between the source and Earth.

The baryonic Tully Fisher relation (bTFR) provides an empirical connection between baryonic mass and dynamical mass (measured by the maximum rotation velocity) for galaxies. Due to the impact of baryonic feedback in the shallower potential wells of dwarf galaxies, the bTFR is predicted to turn down at low masses from the extrapolated power-law relation at high masses. The low-mass end of the bTFR is poorly constrained due to small samples and difficulty in connecting the galaxy's gas kinematics to its dark matter halo. Simulations can help us understand this connection and interpret observations. We measure the bTFR with 66 dwarf galaxies from the Marvel-ous and Marvelous Massive Dwarfs hydrodynamic simulations. Our sample has M$_\star = 10^6-10^9$ M$_\odot$, and is mostly gas dominated. We compare five velocity methods: V$_\text{out,circ}$ (spatially resolved mass-enclosed), V$_\text{out,mid}$ (spatially resolved midplane gravitational potential), and unresolved HI linewidths at different percentages of the peak flux (W$_\text{10}$, W$_\text{20}$, and W$_\text{50}$). We find an intrinsic turndown in the bTFR for maximum halo speeds $\lesssim 50$ km s$^{-1}$ (or total baryonic mass, M$_\text{bary}\lesssim 10^{8.5}$ M$_\odot$). We find that observing HI in lower-mass galaxies to the conventional surface density limit of 1M$_\odot$pc$^{-2}$ is not enough to detect a turndown in the bTFR; none of the HI velocity methods (spatially resolved or unresolved) recover the turndown, and we find bTFR slopes consistent with observations of higher-mass galaxies. However, we predict that the turndown can be recovered by resolved rotation curves if the HI limit is $\lesssim 0.08$ M$_\odot$ pc$^{-2}$, which is within the sensitivity of current HI surveys like FEASTS and MHONGOOSE.

The sub-Neptune planets have no solar system analogues, and their low bulk densities suggest thick atmospheres containing degenerate quantities of volatiles and H/He, surrounding cores of unknown sizes. Measurements of their atmospheric composition can help break these degeneracies, but many previous studies at low spectral resolution have largely been hindered by clouds or hazes, returning muted spectra. Here, we present the first comprehensive study of a short-period sub-Neptune using ground-based, high-resolution spectroscopy, which is sensitive to the cores of spectral lines that can extend above potential high altitude aerosol layers. We observe four CRIRES+ $\textit{K}$-band transits of the warm sub-Neptune GJ 3090 b (T$_{\text{eq}}$ = 693$\pm$18 K) which orbits an M2V host star. Despite the high quality data and sensitivity to CH$_4$, H$_2$O, NH$_3$, and H$_2$S, we detect no molecular species. Injection-recovery tests are consistent with two degenerate scenarios. First, GJ 3090 b may host a highly metal-enriched atmosphere with $>\,$150 Z$_{\odot}$ and mean molecular weight $>\,$7.1 g mol$^{-1}$, representing a volatile dominated envelope with a H/He mass fraction $x_{\text{H/He}} < 33\%$, and an unconstrained aerosol layer. Second, the data are consistent with a high altitude cloud or haze layer at pressures $<\,$3.3$\times$10$^{-5}~$bar, for any metallicity. GJ 3090 b joins the growing evidence to suggest that high metallicity atmospheres and high altitude aerosol layers are common within the warm (500$~<T_{\text{eq}}<~$800 K) sub-Neptune population. We discuss the observational challenges posed by the M-dwarf host star, and suggest observing strategies for transmission spectroscopy of challenging targets around M-dwarfs for existing and ELT instrumentation.

We present the first study dedicated to measuring the timescales for black hole accretion and jet launch in a quasar at the edge of Reionization, PSO J352.4034-15.3373 at z = 5.832 $\pm$ 0.001. Previous work presented evidence of the strong radio synchrotron emission from the jet affecting the host galaxy dust-dominated continuum emission at $\nu_{\rm rest}=683$ GHz ($\nu_{\rm obs}=100$ GHz), implying a break in the synchrotron spectrum. In this work, we present quasi-simultaneous observations at 1.5\, GHz - 42\,GHz with the Karl G. Jansky Very Large Array (VLA), and derive a frequency break at $\nu^{\rm break}_{\rm rest} = 196.46$ GHz ($\nu^{\rm break}_{\rm obs} = 28.76$ GHz). Modeling these observations, we calculate the jet spectral aging from the cooling of electrons to be $t_{\mathrm{spec}}\sim 580$ yr. From this measurement, we approximate the dynamical age $t_{\mathrm{dyn}}$ to be $\sim2,000$ yr, implying a recent jet ejection. We compare the jet timescale to the quasar's lifetime ($t_{\mathrm{Q}}$) that indicates the duration of the latest black hole accretion event and is derived from the proximity zone size in the rest-UV/optical spectrum. However, a ghostly Damped Ly$\alpha$ (DLA) system affects this measurement yielding an upper limit of $t_{\mathrm{Q}} \lesssim 10^4$ yr, consistent with the jet lifetime and indicative of a young quasar. This suggests that the triggering of a UV-bright quasar phase may occur within comparable timescales as the launch of a relativistic radio jet. Therefore, we may be witnessing an early stage of black hole and jet interactions in a quasar during the first gigayear of the universe.

Radiation is crucial not only for observing astrophysical objects, but also for transporting energy and momentum. However, accurate on-the-fly radiation transport in astrophysical simulations is challenging and computationally expensive. Here we introduce AREPO-IDORT (Implicit Discrete Ordinates Radiation Transport), a scheme coupled to the explicit magnetohydrodynamic (MHD) solver in the 3D moving-mesh code AREPO. The discrete ordinates scheme means we directly solve for the specific intensities along discrete directions. We solve the time-dependent relativistic radiation transport equation via an implicit Jacobi-like iterative finite-volume solver, which overcomes the small radiation time-steps needed by explicit methods. Compared to commonly-used moment-based methods, e.g. flux-limited diffusion or M1 closure, this scheme has the advantage of correctly capturing the direction of radiation in both optically-thick and thin regions. It is based on the scheme by Jiang 2021 for the adaptive mesh refinement code ATHENA++, but we generalize the scheme to support (1) an unstructured moving-mesh, (2) local time-stepping, and (3) general equations of state. We show various test problems that commonly-used moment-based methods fail to reproduce accurately. To apply the scheme to a real astrophysics problem, we show the first global 3D radiation hydrodynamic simulation of the entire convective envelope of a red supergiant star. (abridged) For this problem, the radiation module only takes less than half of the total computational cost. Our current scheme assumes grey radiation, is first-order accurate in both time and space (abridged). We expect our scheme will enable more accurate multi-scale radiation MHD simulations involving supersonic bulk motions, ranging from planet formation in protoplanetary disks, stars and associated transients, to accretion flows near black holes.

Understanding the creation of relativistic jets originating from active galactic nuclei, require a thorough understanding of the accompanying plasma instabilities. Our high sensitivity, high resolution, global very long baseline interferometry observations of the jet in the radio galaxy 3C 84 enable us to study its inner morphology, which resembles a thread-like pattern. We find that this pattern can be described by a Kelvin-Helmholtz instability, consisting of four instability modes. Our model favours a jet described by a Mach number of $M_\textrm{j} = 5.0\pm1.7$ and a sound speed of $\alpha_\textrm{j} = 0.14\pm0.06$. With it, we are able to describe the internal structure of 3C 84 and to tentatively connect the origin of the instability to accretion disc activity.

The physico-chemical properties of the planetary nebula (PN) NGC 3242 are investigated in both 1D and 2D, using Integral Field Unit (IFU) data. This PN has a complex morphology with multiple shells and contains a pair of structures with a lower degree of ionization compared to the main nebular components. These structures are known as low ionization structures (LISs), and their origin is still a mystery. With the capabilities provided by IFU spectroscopy, we aim to gain a better understanding of the behavior of nebular properties in the LISs. Data from the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT) were used in order to perform a spatially resolved physico-chemical analysis of NGC 3242 both in 2D, through the analysis of emission line maps, and in 1D, simulating long-slit spectroscopy, with pseudo-slits. Through the deeper investigation of MUSE data, we detect new structures perpendicular to the pair of LISs of NGC 3242, which are mainly seen in the light of [S III] and [N II]. In addition, two arc-like structures are revealed. Moreover, an inner jet-like structure is found through its [Fe III] emission. The interaction of the jet with the rim may be related to the formation of knots and blobs. The higher value of Te, is estimated from the [S III] diagnostic lines, followed by Te ([N II]), Te(H I) and finally Te (He I). In all cases, Te is higher at the inner nebular structures. Regarding electron density, ne, is lower at the LISs, while an increase is observed at the nebular rim. Diagnostic diagrams confirm that NGC 3242 is a highly ionized nebula. Moreover, the MUSE data unveiled for the first time in this PN, the atomic line [C I] {\lambda}8727, primarily emitted from the LISs. This finding suggests that these structures may consist of a molecular core surrounded by neutral and ionized gas

One of the mechanisms responsible for disk formation around the black holes is Bondi-Hoyle-Lyttleton (BHL) this http URL fact that BHL accretion can be interrupted by various astrophysical phenomena, such as stellar winds or astrophysical jets, makes it crucial to study the behavior of shock cones formed by BHL accretion around black holes once the accretion process is halted. Investigating the new plasma structures that emerge in these scenarios can provide insights into observational results. In this context, a new plasma structure forming around the Kerr black hole has been numerically modeled as a function of the black hole spin parameter and the asymptotic velocity of BHL accretion. The numerical analysis revealed that high spin (a/M=0.9) and supersonic flow ( M > 1) are necessary conditions for low-frequency quasi-periodic oscillations (LFQPOs) formation. On the other hand, the fundamental mode of the high-frequency quasi-periodic oscillations (HFQPOs) are found to be independent of both the black hole spin and asymptotic velocity and are instead governed by general relativistic effects. Additionally, the study demonstrated that for 3:2 and 2:1 resonance states to form, nonlinear couplings needs to be occurred when the black hole rotates rapidly. These couplings produce harmonic frequencies, providing an explanation for the observed quasi-periodic oscillation (QPO) resonances in black hole binaries. These findings align with precession models and nonlinear resonance models, both of which play a crucial role in QPO generation. Finally, the LFQPOs and HFQPOs obtained from numerical simulations are consistent with the observed QPO frequencies in the microquasars GRS 1915+105 and XTE J1550-564, as well as in the AGN REJ1034+396, which harbors a supermassive black hole at its center.

Calibration of instrumental polarization is critical for measuring polarized radio emissions from astrophysical sources to extract the magnetic field information in astrophysical, heliospheric, and terrestrial plasmas. At meter wavelengths, calibration of radio polarimetric observations is particularly challenging because of the scarcity of bright polarized sources due to significant Faraday depolarization. Here, we present a novel formalism for polarization calibration using an unpolarized sky model. The formalism is specifically designed for wide-field, low-frequency instruments like the Murchison Widefield Array (MWA), the LOw Frequency ARray (LOFAR), New Extension in Nançay Upgrading LoFAR (NenuFAR), Owens Valley Radio Observatory - Long Wavelength Array (OVRO-LWA), low-frequency telescope of the Square Kilometre Array Observatory (SKAO-low), etc. By leveraging the apparent polarization of the unpolarized sky induced by the polarized primary beam of the radio telescope, this method avoids dependence on bright polarized calibrators. It is also immune to ionospheric Faraday rotation. The validation of the approach via MWA observations confirms the accuracy of the method. This formalism provides a robust framework for low-frequency polarization calibration. It addresses the longstanding calibration challenges and advances the field of low-frequency polarimetry by enabling polarization studies of astrophysical radio sources.

In deep, ground-based imaging, about 15%-30% of object detections are expected to correspond to two or more true objects - these are called ``unrecognized blends''. We use Machine Learning algorithms to detect unrecognized blends in deep ground-based photometry using only catalog-level information: colors, magnitude, and size. We compare the performance of Self Organizing Map, Random Forest, k-Nearest Neighbors, and Anomaly Detection algorithms. We test all algorithms on 9-band ($uBVri^{+}z^{++}YJH$) and 1-size (flux_radius in $\textit{i}$-band) measurements of the ground-based COSMOS catalog, and use COSMOS HST data as the truth for unrecognized blend. We find that 17% of objects in the ground-based COSMOS catalog are unrecognized blends. We show that some unrecognized blends can be identified as such using only catalog-level information; but not all blends can be easily identified. Nonetheless, our methods can be used to improve sample purity, and can identify approximately 30% to 80% of unrecognized blends while rejecting 10% to 50% of all detected galaxies (blended or unblended). The results are similar when only optical bands ($uBVri^{+}z^{++}$) and the size information is available. We also investigate the ability of these algorithms to remove photo-z outliers (identified with spectroscopic redshifts), and find that algorithms targeting color outliers perform better than algorithms targeting unrecognized blends. Our method can offer a cleaner galaxy sample with lower blending rates for future cosmological surveys such as the Legacy Survey of Space and Time (LSST), and can potentially improve the accuracy on cosmological parameter constraints at a moderate cost of precision.

Pyramid wavefront sensors (PWFSs) are the preferred choice for current and future extreme adaptive optics (XAO) systems. Almost all instruments use the PWFS in its modulated form to mitigate its limited linearity range. However, this modulation comes at the cost of a reduction in sensitivity, a blindness to petal-piston modes, and a limit to the sensor's ability to operate at high speeds. Therefore, there is strong interest to use the PWFS without modulation, which can be enabled with nonlinear reconstructors. Here, we present the first on-sky demonstration of XAO with an unmodulated PWFS using a nonlinear reconstructor based on convolutional neural networks. We discuss the real-time implementation on the Magellan Adaptive Optics eXtreme (MagAO-X) instrument using the optimized TensorRT framework and show that inference is fast enough to run the control loop at >2 kHz frequencies. Our on-sky results demonstrate a successful closed-loop operation using a model calibrated with internal source data that delivers stable and robust correction under varying conditions. Performance analysis reveals that our smart PWFS achieves nearly the same Strehl ratio as the highly optimized modulated PWFS under favorable conditions on bright stars. Notably, we observe an improvement in performance on a fainter star under the influence of strong winds. These findings confirm the feasibility of using the PWFS in its unmodulated form and highlight its potential for next-generation instruments. Future efforts will focus on achieving even higher control loop frequencies (>3 kHz), optimizing the calibration procedures, and testing its performance on fainter stars, where more gain is expected for the unmodulated PWFS compared to its modulated counterpart.

Over the past decades, observational evidence of circumstellar disks around massive protostars has been steadily accumulating. However, there have also been cases of non-detections in high-mass star-forming regions, leaving the role and prevalence of disks around massive protostars still uncertain. We used high-resolution (0.2") NOrthern Extended Millimeter Array (NOEMA) observations to study the 1.3 mm continuum and molecular line emission of five massive dense cores in the Cygnus-X cloud complex. Four cores host 2000-au-scale rotating structures previously identified as disk candidates in lower-resolution SMA observations, while the remaining core with no evidence for a disk serves as a comparison. With a resolution of 300 au, the 1.3 mm continuum emission reveals varying levels of fragmentation in our sample, with fragment radii ranging from 150 to 800 AU. In this work, we confirm the existence of two small, stable disks in Keplerian-like rotation at scales of 500 au out of four previously identified disk candidates from the SMA observations at coarser resolution. The lack of evidence for Keplerian disks in other disk candidates identified from the SMA data suggests that rotational signatures observed at 2000 au scales do not necessarily imply the presence of Keplerian disks at smaller scales. Therefore, higher-resolution and higher-sensitivity observations are essential to definitively identify Keplerian disks on smaller scales.

We present the Deep-learning Transient Astronomical Object (Deep-TAO), a dataset of 1,249,079 annotated images from the Catalina Real-time Transient Survey, including 3,807 transient and 12,500 non-transient sequences. Deep-TAO has been curated to provide a clean, open-access, and user-friendly resource for benchmarking deep learning models. Deep-TAO covers transient classes such as blazars, active galactic nuclei, cataclysmic variables, supernovae, and events of indeterminate nature. The dataset is publicly available in FITS format, with Python routines and Jupyter notebooks for easy data manipulation. Using Deep-TAO, a baseline Convolutional Neural Network outperformed traditional random forest classifiers trained on light curves, demonstrating its potential for advancing transient classification.

Maintaining wavefront stability while directly imaging exoplanets over long exposure times is an ongoing problem in the field of high-contrast imaging. Robust and efficient high-order wavefront sensing and control systems are required for maintaining wavefront stability to counteract mechanical and thermal instabilities. Dark zone maintenance (DZM) has been proposed to address quasi-static optical aberrations and maintain high levels of contrast for coronagraphic space telescopes. To further experimentally test this approach for future missions, such as the Habitable Worlds Observatory, this paper quantifies the differences between the theoretical closed-loop contrast bounds and DZM performance on the High-contrast Imager for Complex Aperture Telescopes(HiCAT) testbed. The quantification of DZM is achieved by traversing important parameters of the system, specifically the total direct photon rate entering the aperture of the instrument, ranging from $1.85 \times 10^6$ to $1.85 \times 10^8$ photons per second, and the wavefront error drift rate, ranging from $\sigma_{drift}$ = 0.3 - 3 $nm/\sqrt{iteration}$, injected via the deformable mirror actuators. This is tested on the HiCAT testbed by injecting random walk drifts using two Boston Micromachines kilo deformable mirrors (DMs). The parameter scan is run on the HiCAT simulator and the HiCAT testbed where the corresponding results are compared to the model-based theoretical contrast bounds to analyze discrepancies. The results indicate an approximate one and a half order of magnitude difference between the theoretical bounds and testbed results.

We introduce NEIDSpecMatch, a tool developed to extract stellar parameters from spectra obtained with the NEID spectrograph. NEIDSpecMatch is based on SpecMatch-Emp and HPFSpecMatch, which estimate stellar parameters by comparing the observed spectrum to well-characterized library spectra. This approach has proven effective for M dwarfs. Utilizing a library of 78 stellar spectra covering effective temperatures from $3000-6000$ K, NEIDSpecMatch derives key parameters, including effective temperature, metallicity, surface gravity, and projected rotational velocity. Cross-validation shows median uncertainties of $\sigma_{T_{\mathrm{eff}}} = 115\,\mathrm{K}$, $\sigma_{[\mathrm{Fe/H}]} = 0.143$, and $\sigma_{\log g} = 0.073$ across 49 orders. We showcase its application by fitting the spectrum of an M-dwarf and discuss its utility across a wide range of spectra observed with NEID. NEIDSpecMatch is pip-installable.

Using optical time series with Telescopi Joan Oró (TJO), Gaia, TESS, and NEOWISE archival data, we performed a variability study on the candidate bloated massive young stellar object (MYSO) IRAS 19520+2759. This is the first time that a bloated star candidate has been tested for the theoretically predicted periodic variability. The source is found to be variable at optical and mid-infrared wavelengths and classified as a long-period variable MYSO. The observed TJO data gives a period of the source of $\sim$ 270$\pm$40 days (in the Rc band) and $\sim$ 270$\pm$50 days (in the Ic band), which is very close to the value predicted by the theoretical Period-Luminosity relation for a bloated young star of $\sim 10^5 L\odot$. Additionally, a large period of $\sim$ 460$\pm$80 days (in the G band) and $\sim$ 440$\pm$70 (in the Rp band) is also visible in the Gaia light curve. The physical parameters of the source, such as mass, radius, and accretion rate, based on the theoretical predictions for the spherical accretion case and corresponding to a period of 270--460 days, are $\sim 24$--28$\,M\odot$, $\sim 650$--900$\,R\odot$ and $\sim (6$--$9)\times10^{-3}\,M\odot yr^{-1}$. However, these numbers are very sensitive to the effective temperatures assumed in the models. Additionally, these values strongly depend on the geometry of accretion and could significantly decrease for the case of a MYSO accreting through a disc. The observed periodic variability, the observed colour trend, and the nature of the variability are found to be consistent with the pulsational model for a bloated MYSO.

We examined dust polarisation within the planetary nebula (PN) K 3-35 using the Atacama Large Millimeter/Submillimeter Array (ALMA). This investigation aimed to identify and trace the magnetic field within the PN, as it potentially plays a crucial role in shaping this bipolar nebula. Our findings include a marginal detection of the polarised region and low fractional polarisation (peaking at 1.4%). Assuming a certain level of validity, we observed well-organised dust grains aligned along the equatorial plane of the PN, indicating a magnetic field alignment with the outflows. The limited polarisation detection at submillimeter wavelengths in this PN and others may be attributed to a pronounced optical depth. However, our K 3-35 analysis with the code DUSTY does not seem to support this idea. We also modelled the SED of K3-35, and our best-fit models included a mixture of silicates and amorphous carbon. The grains of amorphous carbon are less susceptible to alignment with the magnetic field, which could, at least partially, explain the observed low polarisation. The models presented in this article should be considered preliminary, and a more advanced approach is needed for a more complete interpretation of the results.

The mass window of ultralight axion dark matter motivated by suppressing the growth of structure on subgalactic scales, $m\sim 10^{-22}\,\mathrm{eV}$, is now severely constrained by various observation data (e.g. Lyman-$\alpha$ forest). As an attempt to reopen this mass window, we investigate an alternative ultralight dark matter candidate, the complex scalar field dark matter (SFDM). We derive the relativistic hydrodynamics of the complex SFDM in the framework of cosmological perturbation theory. Our formalism contains two novel ingredients uniquely associated with the complex SFDM model: the Eckart frame defined by the conserved Noether current, and the stiff gauge condition, $c_s^2\equiv (\delta P/\delta\rho)|_s=1$. In the Eckart frame, the complex SFDM is effectively an imperfect fluid with a dissipative energy flux, distinguishing itself from axion dark matter. The energy flux can affect the growth of density fluctuations dynamically. Meanwhile, we apply the stiff gauge condition to find new constitutive equations for the complex SFDM. We revisit the homogeneous evolution of the complex SFDM and present illustrative early-stage solutions for perturbations of the complex SFDM in a simplified setting. We demonstrate the effects of varying the model parameters on the evolution of the perturbation variables.

Unstable states of the solar coronal magnetic field structure result in various flare behaviors. In this study, we compared the confined and eruptive flares that occurred under similar magnetic circumstances in the active region 12673, on 2017 September 6, using the twist number, decay index, and height of magnetic field lines to identify observational behaviors of the flare eruption. We investigated the parameters from the magnetic field lines involved in an initial energy release, which were identified from the positions of the core of flare ribbons, i.e., flare kernels. The magnetic field lines were derived by nonlinear force-free field modeling calculated from the photospheric vector magnetic field obtained by the Solar Dynamics Observatory SDO/Helioseismic and Magnetic Imager, and flare kernels were identified from the 1600 angstrom data obtained by the SDO/Atmospheric Imaging Assembly. The twist number of all the magnetic field lines in the confined flare was below 0.6; however, the twist number in seven out of twenty-four magnetic field lines in the eruptive flare was greater than 0.6. These lines were tall. It is found that the decay index is not a clear discriminator of the confined and eruptive flares. Our study suggests that some magnetic field lines in the kink instability state may be important for eruptive flares, and that taller magnetic field lines may promote flare eruption.

There is a candidate electromagnetic counterpart to the binary black hole merger GW190521, identified as ZTF19abanrhr within AGN J124942.3 + 344929. Additionally, GW190514 is proposed as a plausible precursor merger to GW190521 within a hierarchical merger scenario. In this study, we investigate the potential association between GW190514 and GW190521 as a hierarchical triple merger associated with ZTF19abanrhr, taking into account of sky position, distance, and mass of the sources using a Bayesian criterion. Our analysis reveals that the association is favored over a random coincidence, with a log Bayes factor of 16.8, corresponding to an odds ratio of $\sim$$199:1$, assuming an astrophysical prior odds of $10^{-5}$. Notably, when accounting for the primary masses of the two gravitational wave events as potential products of mergers in the AGN formation channel, the Bayes factor increases significantly, further enhancing the preference for this association by a factor of $\sim$$10^2$, corresponding to a log Bayes factor of 21.5 and an odds ratio of $\sim$$2\times10^4:1$. Our results suggest strong evidence for the first hierarchical triple merger associated with an electromagnetic counterpart in the AGN formation channel. This work is crucial for understanding the formation mechanisms of massive black holes, the role of AGNs in hierarchical mergers, and the implications of multi-messenger astronomy.

We conducted an analysis of the continuum during the onset and initial decline phases of the 2023 outburst in transient neutron star low-mass X-ray binary Aql X$-$1 using broadband observations from the \textit{Insight-Hard X-ray Modulation Telescope (Insight-HXMT)} instrument. To determine the most appropriate model for the continuum of this outburst, we employed three models to explore the evolution of the spectral component. These observations revealed that the source transitions from the hard state to the soft state. The disk-corona and sphere-corona models both adequately described the spectra of the hard state, while the double blackbody model became preferable after the hard X-ray emission ($>$25 keV) disappeared during the state transition. In the soft state, the total emission is dominated by changes in the disk and other blackbody components. The combination of the sphere-corona model and the double blackbody model is the most suitable model for this outburst. The results suggest that as the source transitioned into the soft state, the emission from the boundary layer was enhanced, and a hot spot occurred. Notably, we identified two type-I X-ray bursts, one of which exhibited a significant hard X-ray deficit (significance $\sim$ 4.82 $\sigma$), which indicates that \textit{Insight-HXMT} has the capability to capture the evolution of the corona in a single burst.

Sound speed can be an important tool in unraveling the nature of matter that exists at the cores of neutron stars. In this study, we investigate three major classes of equations of state; monotonous, non-monotonous and discontinuous depending on the nature of the sound speed in neutron stars. The monotonous EoS refers to hadronic models, the non-monotonous refers to the quarkyonic or smooth crossover models and discontinuous refers to discontinuous first-order phase transition models. We generate a large ensemble of EoS for three classes with the model agnostic speed of sound interpolation approach. Our main aim is to check which class of EoS is most favoured by present astrophysical bounds. It is seen that although non-monotonous and discontinuous is favoured thermodynamically, the usual neutron star observations like mass-radius, and f-mode oscillation fail to provide a satisfactory result. The universal relations are also seen to be futile as they show considerable spread and significant overlaps among the different classes. The Bayesian analysis shows slight bias towards the non-monotonous model but fails to provide a decisive answer.

Three-dimensional special relativistic magnetohydrodynamic simulations are performed to investigate properties of the downstream turbulence generated by the interaction between a relativistic shock wave and multiple clumps. We analyze the properties of the downstream turbulence by performing the Helmholtz decomposition. It is shown that, in contrast to the non-relativistic shock case, the amplitude of compressive modes is comparable to that of solenoidal modes for the relativistic shock. In addition, many reflected shocks propagate in the downstream region. The strength of the compressive mode, the solenoidal mode, the reflected shock waves, and the amplified magnetic field depend on the amplitude of the upstream density fluctuations. Our simulation results suggest that the wide distribution of the ratio of the magnetic energy to the shock kinetic energy, $\epsilon_B$, in gamma-ray burst afterglows is due to the diversity of the gamma-ray burst environment. Furthermore, the inhomogeneity of density around high-energy astrophysical objects affects the spectrum of accelerated particles because the reflected shock and turbulence can inject and accelerate non-thermal particles in the shock downstream region. The probability distribution of the downstream quantities, power spectra of turbulence, and vortex generation are also analyzed and discussed in this work.

We present a comprehensive analysis of the preheating dynamics and associated gravitational wave signatures in the Higgs--$R^2$ inflationary model. Using lattice simulations, we investigate the post-inflationary evolution of the system across the parameter space, covering both the Higgs-like and $R^2$-like scenarios. We demonstrate that the efficiency of preheating is significantly dependent on the nonminimal coupling parameter $\xi$. As the $\xi$ parameter increases, moving from the $R^2$-like regime to the Higgs-like regime, we observe more efficient preheating. Through detailed numerical computations, efficient preheating is shown to lead to stronger gravitational wave production. The amplitude of the gravitational wave spectrum varies by several orders of magnitude as we move from the $R^2$-like regime to the Higgs-like regime. The resultant gravitational wave signatures can serve as a potential observational probe to distinguish between different parameter regimes of the Higgs--$R^2$ model.

Observations conducted with H.E.S.S. at high energies have led to the discovery of numerous gamma-ray sources in the Galactic plane at TeV energies. One of these sources, HESS J1832-085, has been suggested to be a pulsar wind nebula (PWN); however, its nature is not yet fully understood. In this work, we analyze Suzaku data to investigate the X-ray spectral properties of HESS J1832-085. We found that the X-ray spectra are highly absorbed and well-represented by a power-law model with a photon index of $\Gamma \sim 1.5$, and an unabsorbed X-ray flux of $F_{\rm X} \sim 0.3 \times 10^{-11}$ erg cm$^{-2}$ s$^{-1}$ in the 2-10 keV energy band. The gamma-ray flux is approximately 66 times higher than the X-ray flux. Based on our X-ray analysis, we discuss the origin of the source HESS J1832-085. We propose that the PWN scenario is possible, although several issues still need to be resolved.

Using the 3D density distribution derived from the 3D dust map of the solar neighborhood, the gravitational potential is obtained by solving the Poisson equation, from which the tidal tensor is computed. In the optimal decomposition, the external tidal tensor follows the same formalism as that of a point mass. The average tidal strength of the clouds, derived from both tidal tensor analysis and pixel-by-pixel computation, shows consistent results. The equivalent velocity dispersion of the clouds, estimated from the average tidal strength, is comparable in magnitude to the velocity dispersion measured from CO (1-0) line emission. This suggests that tidal effects from surrounding material may play a significant role in driving velocity dispersion within the clouds. Future studies should carefully consider these tidal effects in star-forming regions.

Historically, various methods have been employed to understand the origin of the elements, including observations of elemental abundances which have been compared to Galactic Chemical Evolution (GCE) models. It is also well known that 1D Local Thermodynamic Equilibrium (LTE) measurements fail to accurately capture elemental abundances. Non-LTE (NLTE) effects may play a significant role, and neglecting them leads to erroneous implications in galaxy modelling. In this paper, we calculate 3D NLTE abundances of seven key iron-peak and neutron-capture elements (Mn, Co, Ni, Sr, Y, Ba, Eu) based on carefully assembled 1D LTE literature measurements, and investigate their impact within the context of the OMEGA+ GCE model. Our findings reveal that 3D NLTE abundances are significantly higher for iron-peak elements at [Fe/H]< -3, with (for the first time ever) [Ni/Fe] and (confirming previous studies) [Co/Fe] on average reaching 0.6-0.8 dex, and [Mn/Fe] reaching -0.1 dex, which current 1D core-collapse supernova (CCSN) models cannot explain. We also observe a slightly higher production of neutron-capture elements at low metallicities, with 3D NLTE abundances of Eu being higher by +0.2 dex at [Fe/H]= -3. 3D effects are most significant for iron-peak elements in the very metal-poor regime, with average differences between 3D NLTE and 1D NLTE {reaching} up to 0.15 dex. Thus, ignoring 3D NLTE effects introduces significant biases, so including {them} should be considered whenever possible.

Weak lensing galaxy surveys are currently undergoing a dramatic revolution as the dawn of the Stage-IV surveys are upon us. Hence, ensuring that our analysis methods are as accurate and precise as the raw data is of upmost importance. This motivated the development of a new implementation of the quadratic maximum likelihood power spectrum estimation technique, the application of the theoretical uncertainties approach to mitigate baryonic feedback biases, and to re-evaluating the criterion from which binary scale cuts are derived when aiming to eliminate baryonic biases. These techniques maximise the available information from weak lensing observations while minimising potential systematic biases, and shows how this PhD thesis contributes to the advancement of weak lensing cosmology.

The peculiar motion of an observer relative to an ideal reference frame at rest with respect to the cosmic background produces boosting effects which modify and transfer at higher multipoles the frequency spectrum of the isotropic background. To mitigate the computational effort needed for accurate theoretical predictions, analytical solutions of a linear system able to evaluate the spherical harmonic expansion coefficients for (analytical or semi-analytical) background representations have been presented, and extended to generic tabulated functions potentially affected by numerical uncertainties. Owing to the dipole spectrum frequency dependence and to precise inter-frequency calibrations, it will be possible to constrain (or even detect) the tiny imprints in the background spectrum from a variety of cosmological and astrophysical processes.

We present results of numerical relativity simulations of core collapse of rotating magnetized white dwarfs (WDs) in three dimension, aiming at discussing the explosion dynamics and associate multi-messenger signals: gravitational waves (GWs), neutrinos, and electromagnetic counterparts. All WDs initiate gravitational collapse due to electron captures and then experience prompt type explosions after the proto neutron star formation. We observe the explosions dominated by a bipolar structure and the emergence of strong spiral waves in rapidly rotating models. The spiral waves facilitate to increase both the explosion energy and ejecta mass, though the final values still fall in the category of low explosion energy supernovae with small ejecta mass. The spiral waves also produce strong GWs, which may expand the horizon distance of such events against GWs up to $\sim 10$ Mpc for third-generation ground-based detectors. Additionally as an intriguing implication, we demonstrate that such accretion or merger induced collapse of WDs might be able to explain some of the rapidly evolving optical transients, such as fast blue optical transients (FBOTs), as previously suggested. Based on the simulation results together with several assumptions, we confirm that the magnetar may account for the brighter side of observed FBOTs, while a combination of ejecta-envelope interaction which can be also followed by radioactive decay of heavy elements synthesized along with the explosion might still explain the fainter branch even in the absence of magnetar formation.

We present a detailed analysis of the low-mass detached eclipsing binary system BB Persei, which contains two K-type stars in a circular orbit with a short period of 0.4856 d. We used light curves from the Transiting Exoplanet Survey Satellite, which observed BB Per in five sectors, to determine its photometric properties and a precise orbital ephemeris. The solution of the TESS light curve in Phoebe results in a detached configuration, where the temperature of the primary component was fixed to $T_1 = 5~300$ K according to Lamost, which gives us $T_2 = 5~050 \pm 50$ K for the secondary. The spectral type of the primary component was derived as K0 and the photometric mass ratio was estimated $q = 0.90$. Slow period changes on the current O-C diagram spanning the past 25 years indicate the presence of a third body orbiting the eclipsing pair with an orbital period of about 22 years. The companion could be a red dwarf of spectral type M6 - M7 with a minimal mass of about 0.1 M$_{\odot}$. The characteristics and temporal variation of the dark region on the surface of the secondary component were estimated.

We consider stably rotating highly magnetised neutron stars and glitching pulsars. We discuss the prospects for detecting continuous gravitational waves from these sources below 20 Hz with next-generation ground-based facilities such as the Einstein Telescope and Cosmic Explorer and space-based observatories such as DECIGO and Big Bang Observer. We demonstrate that these constitute interesting science targets. We use a robust sensitivity estimation method for future searches based on demonstrated performance. We show that the spin-down upper limit on the gravitational wave amplitude of more than 90% of all highly magnetised pulsars and magnetars suitable for a years-long fully coherent search, exceeds the smallest gravitational wave amplitude estimated detectable with DECIGO and Big Bang Observer. We find that the hidden magnetar candidate PSR J1852+0040 can be detected by Cosmic Explorer if it is emitting at least at 20% of its spin-down luminosity. Finally, post-glitch transient continuous gravitational waves from magnetars are an interesting target for deci-Hz detectors, with all but one of the recorded glitches giving rise to a spin-down limit signal above the smallest detectable level.

Given that secular perturbations in a binary system not only excite high orbital eccentricities but also alter the planetary orbital inclination, the classical Keplerian orbital model is no longer applicable for orbital retrieval. The combination of a dynamical model and observational data is essential for characterizing the configuration and planetary mass in close binaries. We calculate the theoretical radial velocity (RV) signal in the N-body framework and observe a drift in the RV semi-amplitude, which leads to a reduction in the $m$sin$i$ detection threshold by 20 $M_{\oplus}$, with $\sim$ 100% detection probability in the $m_1$sin$i_1$-$a_1$ parameter space. High-precision RV data with an accuracy of 1 m/s can detect such dynamical effects. For four close-in binaries-GJ 86, GJ 3021, HD 196885, and HD 41004, the deviation between the minimum mass derived from the Keplerian and N-body models is found to be $> 0.2 ~ M_{\mathrm{Jup}}$. High-precision astrometric data are also necessary to resolve the 3D orbits and true masses exoplanets. We generate astrometric simulation data with accuracies corresponding to Gaia (57.8 $\mu$as) and the Closeby Habitable Exoplanet Survey (CHES) (1 $\mu$as), respectively. Joint orbit fitting is performed for RV + Gaia and RV + CHES synergy methods. Compared with the fitting results from the astrometry-only model, the synergy models more effectively constrain the range of orbital inclinations. Using simulation data, we derive precise uncertainties for the true planetary mass, which are critical for determining the evolution of planets around binary stars and multi-planet systems.

This work analyzes the energetics of asteroid rubble piles in order to understand what asteroid morphologies should naturally arise from their formation and evolution process. In doing this, a phase diagram is developed that maps out the range of final minimum energy states that a collapsing gravitational aggregate can achieve as a function of total angular momentum and mass distribution. This is developed assuming properties associated with rubble pile asteroids, and can provide insight into the formation and subsequent evolution of contact binaries and orbital binaries in the solar system as an outcome of catastrophic disruptions. The system angular momentum is used as an independent parameter, combined with resulting minimum energy configurations as a simple function of mass morphology of the final system. The configuration of systems with an energy boosted above the minimum energy state are also considered. This paper considers an ideal case, but outlines general results that can be continued for more precise models of distributed granular media modeled using continuum models or using discrete element models.

Thanks to the advent of sensitive and broad bandwidth instrumentation, complex organic molecules (COMs) have been found in a wide variety of interstellar environments, not only in our Galaxy but also in external galaxies up to a redshift of 0.89. The detection of COMs in cold environments such as starless or prestellar cores has challenged our understanding of COM formation and new ideas are being implemented in chemical models and explored in laboratory experiments. At the protostellar stage, the advent of new interferometers such as the Atacama Large Millimeter/submillimeter Array (ALMA) has allowed the mapping of the weak emission of COMs in the protostellar envelopes and protoplanetary disks around both low-mass and high-mass protostars, pinpointing their location and revealing differentiation between the different families of molecules. In this way, thermal and non-thermal desorption mechanisms can be probed, constraining the efficiency of formation of COMs in the gas phase versus on grain surfaces. Some degree of continuity in the COM composition is found from the early to late stages of star formation, suggesting that a significant fraction of COMs are formed at the initial conditions of star formation. For extreme environments such as the Galactic Center, cosmic rays and low-velocity shocks seem to influence the COM composition of low and high-density gas components. The spectral confusion limit will be a major challenge for the detection of new COMs in future spectroscopic surveys. However, low-frequency interferometers targeting sources with low-excitation temperatures may help to overcome this limit.

The search for exoplanets is an active field in astronomy, with direct imaging as one of the most challenging methods due to faint exoplanet signals buried within stronger residual starlight. Successful detection requires advanced image processing to separate the exoplanet signal from this nuisance component. This paper presents a novel statistical model that captures nuisance fluctuations using a multi-scale approach, leveraging problem symmetries and a joint spectral channel representation grounded in physical principles. Our model integrates into an interpretable, end-to-end learnable framework for simultaneous exoplanet detection and flux estimation. The proposed algorithm is evaluated against the state of the art using datasets from the SPHERE instrument operating at the Very Large Telescope (VLT). It significantly improves the precision-recall trade-off, notably on challenging datasets that are otherwise unusable by astronomers. The proposed approach is computationally efficient, robust to varying data quality, and well suited for large-scale observational surveys.

Heartbeat (HB) stars exhibit pulses in their light curves once per orbit due to ellipsoidal distortions from strong tides at periapse. We analyze the population of HB stars in the Magellanic Clouds captured by the OGLE survey, and provide broadband spectral energy distribution fitting to estimate physical properties of the HB stars. The HB stars span a wide range of luminosities, radii, and effective temperatures. However, we find that they cluster near loci of strong tidal influence at periapse, indicating that in nearly all cases, strong tides are indeed responsible for their photometric variability. HB stars tend to populate regions away from the main sequence, where stellar evolution is particularly rapid. We examine the distribution of tidal amplitudes, and show that these can be interpreted through a simplified model of radius growth through stellar evolution and orbital circularization through linear tidal dissipation. When we compare rates of tidal dissipation, we find differences between the modeled modified tidal quality factor among hot ($T_{\rm eff}>6250$~K), $Q_\ast' \gtrsim 10^7$, and cool ($T_{\rm eff}<6250$~K), $Q_\ast'\sim 10^5$, stars, which is qualitatively consistent with models of efficient tidal dissipation in the convective envelopes of cool stars. We find that this model can reproduce the observed observed locations and amplitudes of HB stars in the Hertzsprung-Russell diagram. The hot stars, in particular, extend to amplitudes near, but not beyond, the threshold for nonlinear tidal wave breaking on stellar surfaces, suggesting a physical saturation of tidal amplitudes at this threshold.

Individually identified binary systems of very massive stars define fixed points on possible evolutionary pathways that begin with extreme star formation and end in either coalescence of compact remnants or complete disruption as pair-production supernovae. The LMC star Melnick 39 in the Tarantula Nebula is revealed to be an eccentric ($e = 0.618\pm0.014$) binary system of reasonably long period from time-series analysis of Chandra T-ReX X-ray observations. Its X-ray luminosity scales with the inverse of the binary separation, as expected for colliding-wind binaries in the adiabatic regime. The inclusion of optical time-series spectroscopy from the VLT FLAMES Tarantula Survey and archival HST spectroscopy confirms Melnick 39 as a double-lined O2.5If/WN6+O3V-III spectroscopic binary with orbital period near 648 days. We obtain a mass ratio of $q = 0.76 \pm 0.06$, and minimum dynamical masses of $105\pm11$ and $80\pm11$ $M_{solar}$ for the O2.5If/WN6 and O3V-III components, plus photometric evidence for an orbital inclination near 90 degrees. Disentangled spectroscopy allows the physical and wind properties of the primary to be determined, including $T_{\ast}$ = 44 kK, $\log L/L_{solar}$ = 6.2, $\log \dot{M}/M_{solar}$ yr$^{-1}$ = $-5.0$. Its dynamical mass agrees closely with $109 M_{solar}$ obtained from the mass-luminosity relation of very massive stars.

Aims. We study the behaviour of Cl abundance and its ratios with respect to O, S and Ar abundances in a sample of more than 200 spectra of Galactic and extragalactic H ii regions and star-forming galaxies (SFGs) of the local Universe. Methods. We use the DEep Spectra of Ionised REgions Database (DESIRED) Extended project (DESIRED-E) that comprises more than 2000 spectra of H ii regions and SFGs with direct determinations of electron temperature ($T_e$). From this database we select those spectra where it is possible to determine the Cl$^{2+}$ abundance and whose line ratios meet certain observational criteria. We calculate the physical conditions and Cl, O, S and Ar abundances in an homogeneous manner for all the spectra. We compare with results of photoionisation models to carry out an analysis of which is the most appropriate $T_e$ indicator for the nebular volume where Cl$^{2+}$ lies, proposing a scheme that improves the determination of the Cl$^{2+}$ abundance. We compare the Cl/O ratios obtained using two different ionisation correction factor (ICF) schemes. We also compare the nebular Cl/O distribution with stellar determinations. Results. Our analysis indicates that the ICF scheme proposed by Izotov et al. (2006) better reproduces the observed distributions of the Cl/O ratio. We find that the log(Cl/O) vs. 12+log(O/H) and log(Cl/Ar) vs. 12+log(Ar/H) distributions are not correlated in the whole metallicity range covered by our objects indicating a lockstep evolution of those elements. In contrast, the log(Cl/S) vs. 12+log(S/H) distribution shows a weak correlation with a slight negative slope.

This paper presents the first combined detections of CO$_2$, CO, XCN and water ices beyond the local Universe. We find gas-phase CO in addition to the solid phase CO. Our source, SSTXFLS J172458.3+591545, is a $z=0.494$ star-forming galaxy which also hosts a deeply obscured AGN. The profiles of its ice features are consistent with those of other Galactic and local galaxy sources and the implied ice mantle composition is similar to that of even more obscured sources. The ice features indicate the presence of a compact nucleus in our galaxy and allow us to place constraints on its density and temperature ($n>10^5$cm$^{-3}$ and $T=20-90K$). We infer the visual extinction towards this nucleus to be $A_V\approx6-7$. An observed plot of $\tau_{Si}$ vs. $\tau_{CO2}/\tau_{Si}$ can be viewed as a probe for both the total dustiness of a system as well as the clumpiness of the dust along the line of sight. This paper highlights the potential of using {\sl JWST} MIRI spectra to study the dust composition and geometric distribution of sources beyond the local Universe.

The strong gravitational field of a black hole bends light, forming multi-level images, yet extracting precise spacetime information from them remains challenging. In this study, we investigate how gravitational lensing leaves unique and detectable signatures in black hole movies using autocorrelation analysis. By examining the two-dimensional autocorrelation of a movie depicting a hotspot orbiting a Kerr black hole, as viewed by a near-axis observer, we identify a persistent secondary peak structure induced by gravitational lensing. Notably, these secondary peaks converge to a fixed point in the time-lag domain, largely independent of the hotspot's orbital radius. This key property suggests that combining future black hole flare observations with advanced autocorrelation analysis could effectively disentangle lensing effects from orbital dynamics, enabling direct measurement of black hole parameters. Our findings establish autocorrelation as a powerful tool for probing spacetime geometry, offering new insights into gravitational physics through time-resolved black hole images.

Context. Galactic halos host faint substructures, such as stellar streams and shells, which provide insights into the hierarchical assembly history of galaxies. To date, such features have been identified in external galaxies by visual inspection. However, with the advent of larger and deeper surveys and the associated increase in data volume, this methodology is becoming impractical. Aims. Here we aim to develop an automated method to detect low surface brightness features in galactic stellar halos. Moreover, we seek to quantify its performance when considering progressively more complex data sets, including different stellar disc orientations and redshifts. Methods. We develop the Stream Automatic Detection with Convolutional Neural Networks, SAD-CNN. This tool is trained on mock surface brightness maps obtained from simulations of the Auriga Project. The model incorporates transfer learning, data augmentation and balanced datasets to optimise its detection capabilities at surface brightness limiting magnitudes ranging from 27 to 31 mag arcsec^-2. Results. The iterative training approach, coupled with transfer learning, allowed the model to adapt to increasingly challenging datasets, achieving precision and recall metrics above 80% in all considered scenarios. The use of a well-balanced training dataset is critical for mitigating biases, ensuring that the CNN accurately distinguishes between galaxies with and without streams. Conclusions. SAD-CNN is a reliable and scalable tool for automating the detection of faint substructures in galactic halos. Its adaptability makes it well-suited for future applications, including the analysis of data from upcoming large astronomical surveys (such as LSST, JWT).

Dark matter heating in planets has been proposed as a potential probe for dark matter detection. Assuming near-equilibrium, we find that dark matter energy input raises planetary temperature and accelerates rotation. The energy distribution depends on planetary properties and external input, suggesting that previous studies have overestimated its heating effect. When dark matter density is high, planetary rotation stabilizes earlier and is primarily governed by dark matter. For Earth, our model predicts a 0.015 K temperature increase and a $10^{-8}$ rad/s angular velocity increase over 100 years.

Positive lags between the arrival time of different photon energies are commonly observed in Gamma-Ray Bursts (GRBs), where soft photons lag behind harder ones. However, some GRBs exhibit the opposite behavior. In particular, Fermi LAT observations have revealed that high-energy photons often have a delayed onset. We explore spectral lags as a tool to identify emission components, analyzing Fermi GBM and LAT Low Energy (LLE) data. Using the Discrete Correlation Function method, we compute spectral lags in four energy bands (10 keV-100 MeV) for 70 GRBs from the LLE Catalog. Lags between 10 keV and 1 MeV are mostly positive (76%), possibly due to a hard-to-soft spectral evolution. However, lags between the LLE band (30-100 MeV) and GBM (10-100 keV) vary: 40% are positive, while 37% are negative. These negative lags suggest the delayed emergence of an additional high-energy component. Spectral analysis of 56 GRBs reveals that negative lags correspond to an LLE spectral index harder than the GBM high-energy power law. LLE spectral lags can thus serve as a diagnostic tool to identify and characterize emission components, emphasizing the importance of combining temporal and spectral analyses to better understand GRB emission mechanisms.

Context: Rings around giant planets are a common feature of the solar system. Even though solar radiation pressure is known to destabilize rings by exciting the orbital eccentricity of its particles, the Centaur Chariklo (and possibly Chiron), the dwarf planet Haumea, and trans-Neptunian object Quaoar also host rings of solid material. Aims: We explore the dynamical evolution of rings around spherical Chariklo and Haumea analogs, assuming different particle sizes and tilt angles with respect to the planetary orbital plane of the ring. Methods: The ring dynamics were studied using a GPU-based N-body integrator with an 8th-order Hermite scheme for several thousand years, corresponding to 10 solar orbits. The simulations took into account the gravitational effects of the planet and the Sun, radiation pressure, and the shadow cast by the planet. Results: Two families of rings have been identified depending on the ring tilt angle. Slightly tilted rings (<=40 deg) are unstable under a critical particle size. Highly tilted rings (>=50 deg), however, show instability only for a range of particle sizes that spans 1-10 times the critical size. The planetary shadow reduces the critical size by a factor of five and extends the instability region to 0.1-10 times this newly identified critical size. Conclusions: The stabilization of highly inclined rings occurs because the plane of the ring is forced to be perpendicular to the Solar radiation. As a result, the plane of the ring rotates as the planetary bodies revolves: always facing the sun, like a celestial sunflower. Rings which are closely aligned to the orbital plane of the host planet, such as Haumea and Quaoar, presumably consist of particles with a size at least 1-4 um. However, particles in the rings which are highly tilted, like that around Chariklo and Chiron, should consist of particles <=2.5-15 um or >=60-300 um.

A knowledge of polarization properties of coherent radio waves escaping pulsar polar caps is crucial for calculating radiative transfer through the magnetosphere and for obtaining specific predictions of observable radio properties. We describe the pair cascades in the pulsar polar cap, and for the first time, determine the Stokes parameters of the escaping radio waves from first-principle kinetic simulations for a pulsar with an inclination angle of the magnetic axis 60°. Our model provides a quantitative and qualitative explanation of the observed pulsar radio powers and spectra, the pulse profiles and polarization curves, their temporal variability, the strong Stokes L and weak Stokes V polarization components, as well as the fact that linear polarization decreases with frequency and the non-existence of a radius to frequency relationship. We find that the radio emission from the polar cap can produce a diverse range of observed pulsar properties, including single or double peaked profiles. Most of the Stokes V curves from our simulations appear to be antisymmetric, but symmetric curves are also present at some viewing angles. Although the PA swing of the radiation from the polar cap can be fitted by the rotating vector model (RVM) for most viewing angles, the angles obtained from the RVM do not correspond to the angular distance of the observer from the magnetic axis. Instead, the PA is directly related to the plasma flows in the polar cap and not to the dipole geometry of the magnetic field. The observed range of other polarization features, in addition to our results, can be explained by propagation effects which are not part of the simulation. Our simulations demonstrate that pair discharges determine the majority of its typically observed properties. The usage of RVM for estimations of the magnetic field geometry from observations needs to be reevaluated.

We address the chemical link between Terzan 5 (hereafter Ter5) and the Bulge, as probed by the observed distributions of [$\alpha$/Fe] abundance ratios with varying [Fe/H] and by suitable statistical tests to evaluate their significance. We also present a comprehensive review of the kinematic and evolutionary properties of Ter5, based on all the available observational signatures and the scenarios proposed so far in the literature for the formation and evolution of Ter5, based on these observational facts and the recent modeling of its star formation and chemical enrichment history. This analysis confirms the complex nature of this massive stellar system, with robust evidences of a bulge in-situ formation and of a subsequent evolution that cannot be simply explained by a single merging/accretion event of two globulars or a globular and a giant molecular cloud, as proposed in the literature, but it requires a more complex star formation likely accompanied by some self-enrichment.

The derivation of precise stellar ages is considered the current major challenge to reconstruct the chronology of the Milky Way. Color-magnitude diagram (CMD)-fitting offers a robust alternative to individual age determinations via the derivation of dynamically evolved star formation histories (deSFH) and age-metallicity distributions (Gallart et al. 2024). Our new suite of routines, this http URL, specifically developed to analyse Gaia CMDs, produce deSFHs which are robust against sensible changes in the input parameters and extremely precise, providing an unprecedentedly detailed characterization of the successive events of star formation that, since its early evolution, have shaped the current Milky Way. Also important is the fact that, thanks to the high completeness of the Gaia photometric data, this http URL provides the actual number of stars and the mass involved in the different events of star formation. The current analysis of the deSFH for stellar populations within 100 pc of the Sun, as well as for kinematically selected stars in the thin disk, thick disk, and halo, allows us to sketch a tentative picture of Milky Way evolution. The findings indicate that star formation commenced very early in a thick disk, with a small fraction of stars having [M/H]<-0.5 forming more than 12 Gyr ago. This phase culminated in a more prominent 12 Gyr old population with [M/H]~-0.5. Approximately 11 Gyr ago, the merger with GSE triggered an intense burst of star formation, generating most of the thick disk mass and enriching its metallicity to solar levels. Subsequently, the bulk of the star formation in the thin disk started and continues with a somewhat episodic behaviour up to the present time.

Cosmic background radiation, both diffuse and discrete in nature, produced at different cosmic epochs before and after recombination, provides key information on the evolution of cosmic structures. We discuss the main classes of sources that contribute to the extragalactic background light from radio to sub-millimetre wavelenghs and the currently open question on the level of the cosmic radio background spectrum. The redshifted 21cm line signal from cosmological neutral Hydrogen during the primeval phases of cosmic structures as a probe of the cosmological reionisation process is presented, along with the route for confident detection of this signal. We then describe the basic formalism and the feasibility to study via a differential approach, based mainly on dipole analysis, the tiny imprints in the CB spectrum expected from a variety of cosmological and astrophysical processes at work during the early phases of cosmic perturbation and structure evolution. Finally, we discuss the identification of high-redshift sub-millimetre lensed galaxies with extreme magnifications in the Planck maps and their use for the comprehension of fundamental processes in early galaxy formation and evolution.

Asteroid (16) Psyche is a metal-rich body that might record an ancient coherent magnetization if some relict crust or mantle is preserved. Herein, we use magnetohydrodynamic simulations to predict (16) Psyche's field, assuming it has such relicts that were magnetized after nebula dispersal via one of two distinct pathways: i. an early solar wind-induced magnetization imparted after a larger body was impacted, forming the present-day asteroid and ii. a core dynamo magnetization imparted in an asteroid that is either presently largely intact or was a rubble pile. For pathway (i) we find the field to be predominantly dipolar and spin axis-aligned. For pathway (ii) we find the field to be either dipolar and spin axis-misaligned, or highly multipolar. Field topology and orientation may thus reveal key details of the nature and history of (16) Psyche, and our framework is broadly applicable to the study of magnetic fields from other asteroids.

We present new Hubble Space Telescope/Wide Field Camera 3 G141 grism observations for COBRA1411+3415, originally identified as a high-redshift cluster candidate in the Clusters Occupied by Bent Radio AGN (COBRA) survey using radio, infrared, and optical data. We spectroscopically identify seven cluster members within a 0.5 Mpc radius with grism redshifts in the range $1.8006 \leq z_{grism} \leq 1.8175$, consistent with COBRA1411+3415 being a high-redshift cluster with a mean redshift of $\langle z_{grism}\rangle = 1.8106 \pm 0.0006$. The detection of seven galaxies within this small redshift range is significant above the background distribution of galaxies at the level of 5$\sigma$. The line-of-sight velocity dispersion of the cluster is found to be $\sigma_{\parallel} = 701^{+347}_{-138}$ km/s with a virial mass of $M_{200} \approx 2.2^{+3.3}_{-1.3}\times 10^{14}$ M$_{\odot}$. However, the mass may be lower if the cluster is still in formation. In projected phase-space, we also identify two possible infalling members of COBRA1411+3415 and two additional structures at $z\sim 1.73$ and $z\sim 1.88$. The similar spatial distributions and small projected separation from the main cluster suggest they could be a part of the same large-scale filament and together may form a protocluster system that could eventually merge to form a single, massive cluster. COBRA1411+3415 is the highest redshift cluster to be spectroscopically confirmed using a bent, double-lobed radio source as a cluster tracer.

Recent optical surveys have found a rare population of low-$z$ ($z \sim 0.01 - 0.06$) extreme star-forming galaxies (xSFGs) that are the most metal-poor galaxies with strong emission lines, extremely high specific star-formation rate and low stellar mass. Using the Hubble Space Telescope Cosmic Origins Spectrograph, it was found that xSFGs are strong Ly$\alpha$ emitters. Their Ly$\alpha$ properties indirectly suggest that they are also strong Ly continuum (LyC) leakers. This along with several other global properties makes them similar to the recently found $z > 6$ reionization-era star-forming (SF) galaxies using the James Webb Space Telescope. Here we aim to study the radio spectral energy distribution (radio-SED) of $8$ xSFGs to understand mechanisms behind the extreme nature of SF in these galaxies, particularly their high ionisation state and possible strong LyC leakage. We present new radio continuum (RC) observations using the upgraded Giant Metrewave Radio Telescope (uGMRT) at Band-$5$ ($1060 - 1460$ MHz) together with the Karl G. Jansky Very Large Array (VLA) at S- ($2-4$ GHz), C- ($4-8$ GHz), X- ($8-12$ GHz) and Ku- ($12-18$ GHz) bands and archival LOw Frequency ARray (LOFAR) survey data between $120-168$ MHz for a few sources. The radio-SED of xSFGs is flat between $6-15$ GHz and shows a strong evidence for a turnover at lower frequencies in the range of $2-10$ GHz. They can be well described using a thermally dominated radio spectra with a free-free absorption component with a high emission measure. Thus the SF complexes in these galaxies are extremely young (pre-SNe stage; below $\sim 5$ Myrs) and the ISM is very dense suggesting a dominance of young massive star clusters. This can explain several extreme star-forming properties, along with the potential leakage of a significant amount of LyC photons in metal-poor extreme starburst galaxies. (abridged)

Most galaxies in the Universe are found in groups, which have various morphologies and dynamical states. Studying how groups evolve is an important step for our understanding in both large-scale structure formation and galaxy evolution. We analysed the system composed by two groups at z = 0.037, NGC 5098, a group dominated by a pair of elliptical galaxies, and NGC 5096, a compact system which appears to be interacting with NGC 5098. We aim to describe its current dynamical state in order to investigate how it fits in our current cosmological framework. Our analysis is based on deep Canada-France-Hawaii Telescope (CFHT/MegaCam) g and r imaging, archival Chandra X-ray data, and publicly available data of the galaxy redshift distribution. We model the surface brightness of the 12 brightest galaxies in the field-of-view and investigate the diffuse intragroup light that we detect. With a redshift sample of 112 galaxies, we study the dynamical states of both groups. We detect low surface brightness diffuse light associated with both galaxy-galaxy interactions and a possible group-group collision. The substructure we found in velocity space indicates a past interaction between both groups. This is further corroborated by the X-ray analysis. We conclude that NGC 5098 and NGC 5096 form a complex system, that may have collided in the past, producing a sloshing observed in X-rays and a large scale diffuse component of intragroup light as well as some important tidal debris.

We present an updated reconstruction of the dark energy equation of state, $w(a)$, using the newly released DESI DR2 Baryon Acoustic Oscillation (BAO) data in combination with Pantheon+ and DES5Y Type Ia supernovae measurements, respectively. Building on our previous analysis in arXiv:2503.08658, which employed a nonparametric flexknot reconstruction approach, we examine whether the evidence for dynamical dark energy persists with the improved precision of the DESI DR2 dataset. We find that while the overall qualitative structure of $w(a)$ remains consistent with our earlier findings, the statistical support for dynamical dark energy is reduced when considering DESI DR2 data alone, particularly for more complex flexknot models with higher numbers of knots. However, the evidence for simpler dynamical models, such as $w$CDM and CPL (which correspond to $n=1$ and $n=2$ knots respectively), increases relative to $\Lambda$CDM with DESI DR2 alone, consistent with previous DESI analyses. When combined with Pantheon+ data, the conclusions remain broadly consistent with our earlier work, but the inclusion of DES5Y supernovae data leads to an increase of preference for flexknot models with more than two knots, placing $w$CDM and CPL on par with $\Lambda$CDM.

Hypervelocity stars (HVSs) ejected from the Galactic Center (GC) at speeds faster than the Galactic escape velocity are useful tools to provide insight into the Milky Way's dark matter halo. However, most characterizations of HVS orbits assume static models of the Milky Way's gravitational potential. In this work, we assess the influence of the Galactic bar and the Large Magellanic Cloud (LMC) on HVS trajectories, comparing them with those from an axisymmetric potential. We simulate 28,000 HVSs ejected over the last 100 Myr and find that ignoring the bar and LMC can cause their apparent ejection location to drift by up to 100 pc. Applying two standard HVS potential fitting methods to our sample shows that they are unable to perform as designed when non-axisymmetric effects are neglected. We calculate the angle between HVS Galactocentric position and velocity and find the LMC and bar can induce a deflection angle of up to several degrees. Using mock Gaia Data Release 4 observations, however, we show that this deflection is too small in magnitude to be measured in the near future without significantly improved observational uncertainties, particularly in heliocentric distance. Our results emphasize the need to account for the bar and LMC in modeling the Galactic potential using HVSs as a tracer.

Beta Pic is a young, A5V star, known for harbouring a large number of exocomets, which frequently transit the star and produce absorption signatures. The physical and chemical properties of these exocomets can be probed by the recently introduced curve of growth approach, which enables column densities measurements in exocomets using observations in numerous spectral lines. Using this approach, we present a new study of archival spectra of Beta Pic obtained with the HST, the HARPS spectrograph, and at the Mont John University Observatory, aimed at constraining the abundance of refractory ions in Beta Pic exocomets. 29 individual objects are studied, all observed in FeII lines (used as a reference ion) and at least one other species (Ni II, Ca II, Cr II...). We find that the refractory composition of exocomets is overall stable, especially for singly ionised species, and consistent with solar abundances. This validates the use of the curve of growth approach to study exocometary composition. We also show that some ions, such as Ca II, are significantly depleted compared to solar abundances, allowing us to constrain the ionisation state in Beta Pic exocomets. We find that most refractory elements (Mg, Ni, Fe...) are split in similar fractions between their first and second ionisation states, with the exception of Ca, mostly ionized twice. A strong correlation between the Al III/Fe II ratio and radial velocity is also found, showing that the most redshifted exocomets tend to be more ionised. These results open the way for further modelling, in order to better understand the physical processes that influence the composition and shape of exocometary tails.

The detection of young transiting exoplanets represents a new frontier in our understanding of planet formation and evolution. For the population of observed close-in sub-Neptunes, two proposed formation pathways can reproduce their observed masses and radii at $\sim$Gyr ages: the "gas dwarf" hypothesis and the "water world" hypothesis. We show that a sub-Neptune's size at early ages $\lesssim 100$ Myrs is strongly dependent on the bulk mean molecular weight within its envelope. As a result, gas dwarfs and water worlds should diverge in size at early ages since the mean molecular weight of gas dwarf envelopes is predicted to be smaller than that of water worlds. We construct population models under both scenarios that reproduce Kepler demographics in the age range $\sim1-10$ Gyrs. We find tentative evidence that the gas dwarf model is more consistent with the small population of young exoplanets $< 40$ Myrs from TESS. We show that planet radius is relatively insensitive to planet mass for young, puffy sub-Neptunes, meaning that well-characterised masses are not necessarily required to exploit the effects of mean molecular weight at the population level. We confirm the predicted difference in planet size between the models is also true under mixed-envelope scenarios, in which envelopes consist of mixtures of hydrogen and steam. We highlight that transit surveys of young exoplanets should target the youngest observable stellar clusters to exploit the effects of mean molecular weight.

The gravitational wave (GW) signals from extreme mass-ratio inspirals (EMRIs), a key target for the Laser Interferometer Space Antenna (LISA), will be affected in the presence of dark matter (DM) halos. In this paper we explore whether the effects of DM are detectable by LISA within a fully relativistic framework. We model the massive EMRI component as a nonrotating black hole (BH) surrounded by a DM halo. We compute axial and polar GW fluxes for circular orbits at linear order in the mass ratio for DM density profiles with varying mass and compactness. By comparing the phase evolution with vacuum systems, we find that DM halos can induce dephasings of tens to hundreds of radians over a one-year observation period. We demonstrate that even highly diluted DM distributions can significantly affect the emitted waveforms, and that the resulting GW signals can usually be distinguished from each other. While it is important to generalize these findings to more generic orbits and to spinning BHs, our results suggest that LISA could not only reveal the presence of DM halos, but also discriminate between different halo models.

Frequency response function (FRF) measurements are widely used in Gravitational Wave (GW) detectors, e.g., for the design of controllers, calibrating signals and diagnostic problems with system dynamics. The aim of this paper is to present GraFIT: a toolbox that enables fast, inexpensive, and accurate identification of FRF measurements for GW detectors compared to the commonly used approaches, including common spectral analysis techniques. The toolbox consists of a single function to estimate the frequency response function for both open-loop and closed-loop systems and for arbitrary input and output dimensions. The toolbox is validated on two experimental case studies of the Virgo detector, illustrating more than a factor 3 reduction in standard deviation of the estimate for the same measurement times, and comparable standard deviations with up to 10 times less data for the new method with respect to the currently implemented Spectral Analysis method.

Nuclear ground state and collective excitation properties provide a means to probe the nuclear matter equation of state and establish connections between observables in finite nuclei and neutron stars. Specifically, the electric dipole polarizability, measured with high precision in various neutron-rich nuclei, serves as a robust constraint on the density dependence of the symmetry energy. In this Letter, we employ a class of relativistic energy density functionals in a twofold process: first, to link the electric dipole polarizability from recent experiments to the slope of the symmetry energy, and second, to translate this information into constraints on the tidal deformability and radii of neutron stars, in connection with multimessenger astrophysical observations from pulsars and binary neutron stars. We provide compelling evidence that the electric dipole polarizability represents a key nuclear observable to probe the neutron star properties. By significantly reducing the uncertainties in the mass-radius plane, our findings also align with recent multimessenger observations.

The primary focus of this thesis is the numerical investigation of chaos in Hamiltonian models describing charged particle orbits in plasma, star motions in barred galaxies, and orbits' diffusion in multidimensional maps. We systematically explore the interplay between magnetic and kinetic chaos in toroidal fusion plasmas, where non-axisymmetric perturbations disrupt smooth magnetic flux surfaces, generating complex particle trajectories. Using the Generalized Alignment Index (GALI) method, we efficiently quantify chaos, compare the behavior of magnetic field lines and particle orbits, visualize the radial distribution of chaotic regions, and offer GALI as a valuable tool for studying plasma physics dynamics. We also study the evolution of phase space structures in a 3D barred galactic potential, following successive 2D and 3D pitchfork and period-doubling bifurcations of periodic orbits. By employing the `color and rotation' technique to visualize the system's 4D Poincaré surface of sections, we reveal distinct structural patterns. We further investigate the long-term diffusion transport and chaos properties of single and coupled standard maps, focusing on parameters inducing anomalous diffusion through accelerator modes exhibiting ballistic transport. Using different ensembles of initial conditions in chaotic regions influenced by these modes, we examine asymptotic diffusion rates and time scales, identifying conditions suppressing anomalous transport and leading to long-term convergence to normal diffusion across coupled maps. Lastly, we perform the first comprehensive investigation into the GALI indices for various attractors in continuous and discrete-time dissipative systems, extending the method's application to non-Hamiltonian systems. A key aspect of our work involves analyzing and comparing GALIs' with Lyapunov Exponents for systems exhibiting hyperchaotic motion.

The growing energy demands of HPC systems have made energy efficiency a critical concern for system developers and operators. However, HPC users are generally less aware of how these energy concerns influence the design, deployment, and operation of supercomputers even though they experience the consequences. This paper examines the implications of HPC's energy consumption, providing an overview of current trends aimed at improving energy efficiency. We describe how hardware innovations such as energy-efficient processors, novel system architectures, power management techniques, and advanced scheduling policies do have a direct impact on how applications need to be programmed and executed on HPC systems. For application developers, understanding how these new systems work and how to analyse and report the performances of their own software is critical in the dialog with HPC system designers and administrators. The paper aims to raise awareness about energy efficiency among users, particularly in the high energy physics and astrophysics domains, offering practical advice on how to analyse and optimise applications to reduce their energy consumption without compromising on performance.