The galaxy integrated star-formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) has been proposed as a valuable diagnostic of the mass accumulation in galaxies as being more tightly related to the physics of star-formation (SF) and stellar feedback than other SF indicators. In this paper, we assemble a statistical sample of 230 galaxies observed with JWST in the GLASS and CEERS spectroscopic surveys to estimate Balmer line based dust attenuations and SFRs, and UV rest-frame effective radii. We study the evolution of galaxy SFR and $\Sigma_{\rm SFR}$ in the first 1.5 Billion years of our Universe, finding that $\Sigma_{\rm SFR}$ is mildly increasing with redshift with a linear slope of $0.16 \pm 0.06$. We also explore the dependence of SFR and $\Sigma_{\rm SFR}$ on stellar mass, showing that a SF 'Main-Sequence' and a $\Sigma_{\rm SFR}$ `Main-Sequence' are in place out to z=10, with a similar slope compared to the same relations at lower redshifts. We find that the specific SFR (sSFR) and $\Sigma_{\rm SFR}$ are correlated with the [OIII]5007/[OII]3727 ratio and with indirect estimates of the escape fraction of Lyman continuum photons, hence they likely play an important role in the evolution of ionization conditions and in the escape of ionizing radiation. We also search for spectral outflow signatures in a subset of galaxies observed at high resolution, finding an outflow incidence of $2/11$ ($=20\%^{32\%}_{9\%}$) at z<6, but no evidence at z>6 (<26\%). Finally, we find a positive correlation between A$_V$ and $\Sigma_{\rm SFR}$, and a flat trend as a function of sSFR, indicating that there is no evidence of a drop of A$_V$ in extremely star-forming galaxies between z=4 and 10. This might be at odds with a dust-clearing outflow scenario, which might instead take place at redshifts $z\geq 10$, as suggested by some theoretical models.

With the increasing volume of astronomical data generated by modern survey telescopes, automated pipelines and machine learning techniques have become crucial for analyzing and extracting knowledge from these datasets. Anomaly detection, i.e. the task of identifying irregular or unexpected patterns in the data, is a complex challenge in astronomy. In this paper, we propose Multi-Class Deep Support Vector Data Description (MCDSVDD), an extension of the state-of-the-art anomaly detection algorithm One-Class Deep SVDD, specifically designed to handle different inlier categories with distinct data distributions. MCDSVDD uses a neural network to map the data into hyperspheres, where each hypersphere represents a specific inlier category. The distance of each sample from the centers of these hyperspheres determines the anomaly score. We evaluate the effectiveness of MCDSVDD by comparing its performance with several anomaly detection algorithms on a large dataset of astronomical light-curves obtained from the Zwicky Transient Facility. Our results demonstrate the efficacy of MCDSVDD in detecting anomalous sources while leveraging the presence of different inlier categories. The code and the data needed to reproduce our results are publicly available at this https URL.

Quasi-periodic eruptions (QPEs) are recurring X-ray bursts originating from the vicinity of supermassive black holes, but their driving mechanisms remain under debate. This study analyzes new NICER observations of QPEs in Ansky (a transient event in the nucleus of the galaxy SDSS J1335+0728), taken between January and June 2025. By examining flare durations, peak-to-peak recurrence times, and profiles, we compare the 2025 data with those from 2024 to investigate changes in energy, timescales, and flare shapes. The 2025 QPEs are found to be four times more energetic, with recurrence times of approximately 10 days and flare durations ranging from 2.5 to 4 days, making them both about twice as long as in 2024. Additionally, the flare profiles have become more asymmetric, showing longer decays. We explore different theoretical scenarios to explain the observed properties of the QPEs in Ansky, including evolving stream-disk interactions in an extreme mass-ratio inspiral (EMRI) system as a potential mechanism behind the observed changes in recurrence time and energetics, while also considering alternative models based on mass transfer and accretion disk instabilities. Continued observational efforts will be crucial for unveiling the nature of Ansky.

We apply pre-trained Vision Transformers (ViTs), originally developed for image recognition, to the analysis of astronomical spectral data. By converting traditional one-dimensional spectra into two-dimensional image representations, we enable ViTs to capture both local and global spectral features through spatial self-attention. We fine-tune a ViT pretrained on ImageNet using millions of spectra from the SDSS and LAMOST surveys, represented as spectral plots. Our model is evaluated on key tasks including stellar object classification and redshift ($z$) estimation, where it demonstrates strong performance and scalability. We achieve classification accuracy higher than Support Vector Machines and Random Forests, and attain $R^2$ values comparable to AstroCLIP's spectrum encoder, even when generalizing across diverse object types. These results demonstrate the effectiveness of using pretrained vision models for spectroscopic data analysis. To our knowledge, this is the first application of ViTs to large-scale, which also leverages real spectroscopic data and does not rely on synthetic inputs.

Dipper stars are extrinsically variable stars with deep dimming events due to extended, often dusty, structures produced by a wide range of mechanisms such as collisions, protoplanetary evolution or stellar winds. ASAS-SN has discovered 12 dipper-like objects as part of its normal operations. Here we systematically search the $\sim 5.1$ million ASAS-SN targets with $13

The quasar main sequence (QMS), characterized by the Eigenvector 1 (EV1), serves as a unifying framework for classifying type-1 active galactic nuclei (AGNs) based on their diverse spectral properties. Although a fully self-consistent physical interpretation has long been lacking, our physically motivated 2.5D FRADO (Failed Radiatively Accelerated Dusty Outflow) model naturally predicts that the Eddington ratio ($\dot{m}$) is the primary physical driver of the QMS, with black hole mass ($M_{\rm BH}$) and inclination ($i$) acting as secondary contributors. We employed a dense grid of FRADO simulations of the geometry and dynamics of the broad-line region (BLR), covering a representative range of $M_{\rm BH}$ and $\dot{m}$. For each simulation, we computed the full width at half maximum (FWHM) of the H$\beta$ line under different $i$. The resulting FWHM--$\dot{m}$ diagram closely resembles the characteristic trend observed in the EV1 parameter space. This establishes the role of $\dot{m}$ as the true proxy for the Fe II strength parameter ($R_{\rm Fe}$), and vice versa. Our results suggest that $\dot{m}$ can be regarded as the sole underlying physical tracer of $R_{\rm Fe}$ and should therefore scale directly with it. The $M_{\rm BH}$ accounts for the virial mass-related scatter in FWHM, while $i$ acts as a secondary driver modulating $R_{\rm Fe}$ and FWHM for a given $\dot{m}$ and $M_{\rm BH}$.

We present new frontiers in the modelling of the spectral energy distributions (SED) of active galaxies by introducing the radio-to-X-ray fitting capabilities of the publicly available Bayesian code AGNfitter. The new code release, called AGNfitter-rx, models the broad-band photometry covering the radio, infrared (IR), optical, ultraviolet (UV) and X-ray bands consistently, using a combination of theoretical and semi-empirical models of the AGN and host galaxy emission. This framework enables the detailed characterization of four physical components of the active nuclei: the accretion disk, the hot dusty torus, the relativistic jets/core radio emission, and the hot corona; alongside modeling three components within the host galaxy: stellar populations, cold dust, and the radio emission from the star-forming regions. Applying AGNfitter-rx to a diverse sample of 36 AGN SEDs at z<0.7 from the AGN SED ATLAS, we investigate and compare the performance of state-of-the-art torus and accretion disk emission models on fit quality and inferred physical parameters. We find that clumpy torus models that include polar winds and semi-empirical accretion disk templates including emission line features significantly increase the fit quality in 67% of the sources, by effectively reducing by $2\sigma$ fit residuals in the $1.5-5 \mu \rm m$ and $0.7 \mu \rm m$ regimes.We demonstrate that, by applying AGNfitter-rx on photometric data, we are able to estimate inclination and opening angles of the torus, consistent with spectroscopic classifications within the AGN unified model, as well as black hole mass estimates in agreement with virial estimates based on H$\alpha$. The wavelength coverage and the flexibility for the inclusion of state-of-the-art theoretical models make AGNfitter-rx a unique tool for the further development of SED modelling for AGNs in present and future radio-to-X-ray galaxy surveys.

Dust attenuation in galaxies has often been used as a proxy for the extinction of point sources, such as supernovae, even though this approach ignores fundamental differences between the two cases. We present an analysis of the impact of geometric effects and scattering within dusty media on recovered galaxy dust properties. We use SKIRT, a radiative transfer code, to simulate observations of point sources embedded in dust clouds, as well as spiral and elliptical galaxies. We examine various galaxy morphologies, inclinations, and instrument apertures. We find that in galaxies the scattering of light into the line of sight and the presence of sources at different depths within the galaxy make attenuation fundamentally different from extinction. For a medium with intrinsic extinction slope Rv=3.068, we recover effective attenuation slopes Rv_e ranging from 0.5 to 7, showing that the two quantities are not analogous, even for local resolved observations. We find that Rv_e greatly depends on dust density, galaxy morphology, and inclination, the latter being the most significant. A single simulated galaxy, viewed from different angles, can reproduce the well-known relation between attenuation strength Av_e and Rv_e observed for star-forming galaxy samples. An increase in dust density leads to higher Rv_e across all inclinations, which, assuming a correlation between stellar mass and dust density, explains the increase in Rv_e with mass observed in star-forming galaxies. However, we are unable to explain the differences in Rv_e between star-forming and quiescent high-mass galaxies. We conclude that highly attenuated regions of simulated face-on galaxies yield Rv_e within 10% of the intrinsic extinction slope of the medium, allowing for the distinction of different dust types. For edge-on spirals, however, the median Rv_e for low Av_e regions appears to better approximate the extinction slope.

We present a census of the Compton-thick (CT) active galactic nucleus (AGN) population and the column density ($N_{\rm{H}}$) distribution of AGN in our cosmic backyard using a mid-infrared selected AGN sample within 15 Mpc. The column densities are measured from broadband X-ray spectral analysis, mainly using data from $\textit{Chandra}$ and $\textit{NuSTAR}$. Our sample probes AGN with intrinsic 2-10 keV luminosities of $L_{\rm 2-10, int} = 10^{37}$-$10^{43}$ erg s$^{-1}$, reaching a parameter space inaccessible to more distant samples. We directly measure a 32$^{+30}_{-18}\%$ CT AGN fraction and obtain an $N_{\rm{H}}$ distribution that agrees with that inferred by the $\textit{Swift}$-BAT survey. Restricting the sample to the largely unexplored domain of low-luminosity AGN with $L_{\rm 2-10, int}$ $\leq$ $10^{42}$ erg s$^{-1}$, we found a CT fraction of 19$^{+30}_{-14}\%$, consistent with those observed at higher luminosities. Comparing the host-galaxy properties between the two samples, we find consistent star formation rates, though the majority of our galaxy have lower stellar masses (by $\approx 0.3$ dex). In contrast, the two samples have very different black hole mass ($M_{\rm BH}$) distributions, with our sample having $\approx$1.5 dex lower mean mass ($M_{\rm BH}$$\sim$10$^{6}$$M_\odot$). Additionally, our sample contains a significantly higher number of LINERs and H$_{\rm{II}}$-type nuclei. The Eddington ratio range probed by our sample, however, is the same as $\textit{Swift}$-BAT, although the latter dominates at higher accretion rates, and our sample is more evenly distributed. The majority of our sample with $\lambda_{\rm Edd} \ge$ 10$^{-3}$ tend to be CT, while those with $\lambda_{\rm Edd} <$10$^{-3}$ are mostly unobscured or mildly obscured.

We present our sixth work in a series dedicated to variability studies of active galactic nuclei (AGN) based on the survey of the COSMOS field by the VLT Survey Telescope (VST). Its 54 r-band visits over 3.3 yr and single-visit depth of 24.6 r-band mag make this dataset a valuable scaled-down version that can help forecast the performance of the Rubin Observatory Legacy Survey of Space and Time (LSST). This work is centered on the analysis of the structure function (SF) of VST-COSMOS AGN, investigating possible differences in its shape and slope related to how the AGN were selected, and explores possible connections between the ensemble variability of AGN and black-hole mass, accretion rate, bolometric luminosity, redshift, and obscuration of the source. Given its features, our dataset opens up the exploration of samples ~2 mag fainter than most of the literature to date. We identify several samples of AGN - 677 in total - obtained by a variety of selection techniques which partly overlap. Our analysis compares results for the various samples. We split each sample in two based on the median of the physical property of interest, and analyze differences in the shape and slope of the SF, and possible causes. While the shape of the SF does not change with depth, it is highly affected by the type of AGN (unobscured/obscured) included in the sample. Where a linear region can be identified, we find that the variability amplitude anticorrelates with accretion rate and bolometric luminosity, consistent with previous literature on the topic, while no dependence on black-hole mass emerges from this study. With its longer baseline and denser and more regular sampling, the LSST will allow an improved characterization of the SF and its dependencies on the mentioned physical properties over much larger AGN samples.

Planet formation models suggest that the formation of giant planets is significantly harder around low-mass stars, due to the scaling of protoplanetary disc masses with stellar mass. The discovery of giant planets orbiting such low-mass stars thus imposes strong constraints on giant planet formation processes. Here, we report the discovery of a transiting giant planet orbiting a $0.207 \pm 0.011 M_{\odot}$ star. The planet, TOI-6894 b, has a mass and radius of $M_P = 0.168 \pm 0.022 M_J (53.4 \pm 7.1 M_{\oplus})$ and $R_P = 0.855 \pm 0.022 R_J$, and likely includes$12 \pm 2 M_{\oplus}$ of metals. The discovery of TOI-6894 b highlights the need for a better understanding of giant planet formation mechanisms and the protoplanetary disc environments in which they occur. The extremely deep transits (17% depth) make TOI-6894 b one of the most accessible exoplanetary giants for atmospheric characterisation observations, which will be key for fully interpreting the formation history of this remarkable system and for the study of atmospheric methane chemistry.

We present optical photometric and spectroscopic observations of the peculiar Type Ia supernova ASASSN-20jq/SN 2020qxp. It is a low-luminosity object with a peak absolute magnitude of $M_B=-17.1\pm0.5$ mag. Despite its low luminosity, its post-peak light-curve decline rate ($\Delta m_{15}(B)=1.35\pm0.09$ mag) and color-stretch parameter (sBV>0.82) are similar to normal SNe Ia, making it an outlier in the luminosity-width and luminosity-color-stretch relations. Early light curves suggest a "bump" during the first 1.4 days of explosion. ASASSN-20jq synthesized a low radioactive $^{56}$Ni mass of $0.09\pm0.01M_\odot$. Near-maximum light spectra reveal strong Si II absorption lines, indicating a cooler photosphere than normal SNe Ia, but lack Ti II absorption lines. Unusually strong O I $\lambda$7773 and Ca II near-infrared triplet absorption features are present. Nebular spectra show a strong, narrow forbidden [Ca II] $\lambda\lambda$7291,7324 doublet emission, rarely seen in SNe Ia except in some Type Iax events. Marginal detection of [O I] $\lambda\lambda$6300,6364 doublet emission, which is extremely rare, is observed. Both [Ca II] and [O I] lines are redshifted by $\sim2000$ km/s. A strong [Fe II] $\lambda$7155 emission line with a tilted-top profile, identical to the [Fe II] $\lambda$16433 profile, is also observed. These asymmetric [Fe II] profiles and redshifted [Ca II] and [O I] emissions suggest a high central density white dwarf progenitor undergoing an off-center delayed-detonation explosion mechanism, producing roughly equal amounts of $^{56}$Ni in deflagration and detonation phases. This distinguishes ASASSN-20jq from normal and subluminous SNe Ia. ASASSN-20jq's light curve and spectra do not align with any single SNe Ia subclass but show similarities to 2002es-like objects. Thus, we add it as an extreme candidate within the heterogeneous parameter space of 2002es-like SNe Ia.

Context: The VISTA Variables in the Via Lactea (VVV) and its extension (VVVX) are near-infrared surveys mapping the Galactic bulge and adjacent disk. These data have enabled the discovery of numerous star clusters obscured by high and spatially variable extinction. Most previous searches relied on visual inspection of individual tiles, which is inefficient and biased against faint or low-density systems. Aims: We aim to develop an automated, homogeneous algorithm for systematic cluster detection across different surveys. Here, we apply our method to VVVX data covering low-latitude regions of the Galactic bulge and disk, affected by extinction and crowding. Methods: We introduce the Consensus-based Algorithm for Nonparametric Detection of Star Clusters (CANDiSC), which integrates kernel density estimation, the Density-Based Spatial Clustering of Applications with Noise (DBSCAN), and nearest-neighbour density estimation within a consensus framework. A stellar overdensity is classified as a candidate if identified by at least two of these methods. We apply CANDiSC to 680 tiles in the VVVX PSF photometric catalogue, covering approximately 1100 square degrees. Results: We detect 163 stellar overdensities, of which 118 are known clusters. Cross-matching with recent catalogues yields five additional matches, leaving 40 likely new candidates absent from existing compilations. The estimated false-positive rate is below 5 percent. Conclusions: CANDiSC offers a robust and scalable approach for detecting stellar clusters in deep near-infrared surveys, successfully recovering known systems and revealing new candidates in the obscured and crowded regions of the Galactic plane.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

WASP-4 b is a hot Jupiter exhibiting a decreasing orbital period, prompting investigations into potential mechanisms driving its evolution. We analyzed 173 transit light curves, including 37 new observations, and derived mid-transit timings with EXOFAST, forming the most extensive TTV dataset for this system. Adding 58 literature timings and removing unreliable data, we constructed a TTV diagram with 216 points. Our analysis considered linear, quadratic, and apsidal motion models, with the quadratic model proving to be significantly superior in all model comparison statistics. We found no significant periodic signals in the data. The quadratic model allows us to infer a tidal quality factor of Q' ~ 80,000 from the orbital decay rate if this is due to stellar tides. Theoretical considerations indicate that such efficient dissipation is possible due to internal gravity waves in the radiative core of WASP-4, but only in our models with a more evolved host star, possibly near the end of its main-sequence lifetime, and with a larger radius than the observed one. Our main-sequence models produce only about a third of the required dissipation (Q' ~ 200,000 - 500,000). Therefore, the observed orbital decay can only be explained by a slightly larger or more evolved host, resembling the case for WASP-12. Our findings highlight the need for further stellar modeling and improvement in our current understanding of tidal dissipation mechanisms driving orbital decay in close-in exoplanetary systems.

Strong gravitational magnifications enable to detect faint background sources, resolve their internal structures, and even identify individual stars in distant galaxies. Highly magnified individual stars allow various applications, including studies of stellar populations in distant galaxies and constraining dark matter structures in the lensing plane. However, these applications have been hampered by the small number of individual stars observed, as typically one or a few stars are identified from each distant galaxy. Here, we report the discovery of more than 40 microlensed stars in a single galaxy behind Abell 370 at redshift of 0.725 when the Universe was half of its current age (dubbed ``the Dragon arc''), using James Webb Space Telescope (JWST) observations with the time-domain technique. These events are found near the expected lensing critical curves, suggesting that these are magnified stars that appear as transients from intracluster stellar microlenses. Through multi-wavelength photometry, we constrain stellar types and find that many of them are consistent with red giants/supergiants magnified by factors of hundreds. This finding reveals an unprecedented high occurrence of microlensing events in the Dragon arc, and proves that {\it JWST}'s time-domain observations open up the possibility of conducting statistical studies of high-redshift stars.

Context. The classification of active galactic nuclei (AGNs) is a challenge in astrophysics. Variability features extracted from light curves offer a promising avenue for distinguishing AGNs and their subclasses. This approach would be very valuable in sight of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Aims. Our goal is to utilize self-organizing maps (SOMs) to classify AGNs based on variability features and investigate how the use of different subsets of features impacts the purity and completeness of the resulting classifications. Methods. We derived a set of variability features from light curves, similar to those employed in previous studies, and applied SOMs to explore the distribution of AGNs subclasses. We conducted a comparative analysis of the classifications obtained with different subsets of features, focusing on the ability to identify different AGNs types. Results. Our analysis demonstrates that using SOMs with variability features yields a relatively pure AGNs sample, though completeness remains a challenge. In particular, Type 2 AGNs are the hardest to identify, as can be expected. These results represent a promising step toward the development of tools that may support AGNs selection in future large-scale surveys such as LSST.

Some supernovae such as pair-instability supernovae are predicted to have the duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2m Subaru telescope. HSC is a wide-field (a 1.75 deg2 field-of-view) camera and it can efficiently conduct transient surveys. We observed the same 1.75 deg2 field repeatedly using the g, r, i, and z band filters with the typical depth of 26 mag for 4 seasons (from late 2016 to early 2020). Using these data, we searched for transients lasting for more than a year. Two supernovae were detected in 2 continuous seasons, one supernova was detected in 3 continuous seasons, but no transients lasted for all 4 seasons searched. The discovery rate of supernovae lasting for more than a year with the typical limiting magnitudes of 26 mag is constrained to be 1.4^{+1.3}_{-0.7}(stat.)^{+0.2}_{-0.3}(sys.) events deg-2 yr-1. All the long-lasting supernovae we found are likely Type IIn supernovae and our results indicate that about 40% of Type IIn supernovae have long-lasting light curves. No plausible pair-instability supernova candidates lasting for more than a year are discovered. By comparing the survey results and survey simulations, we constrain the luminous pair-instability supernova rate up to z ~ 3 should be of the order of 100 Gpc-3 yr-1 at most, which is 0.01 - 0.1 per cent of the core-collapse supernova rate.