Papers

The Zwicky Transient Facility Data Release 2 (ZTF DR2) includes a total of 3,628 Type Ia supernovae (SNe~Ia), providing the largest and most complete sample of spectroscopically confirmed SNe~Ia at low redshift to date. In this paper, we present a photometric and spectroscopic analysis of 124 subluminous SNe~Ia, the largest sample of spectroscopically classified subluminous SNe~Ia observed with a single instrument, comprising 87 91bg-like, 12 86G-like, 18 04gs-like, and 7 02es-like events. We complement the published DR2 SALT2 light-curve parameters with new parameters obtained using template-based fits from SNooPy. Expansion velocities and pseudo-equivalent widths pEW of key spectral features are measured using Spextractor, and spectral averages are constructed for each subluminous subtype, binned by phase. We also analyze the host galaxy environments, both global and local, in terms of $g - z$ color, stellar mass, and directional light radius $d_{DLR}$. We find that all subluminous SNe~Ia (except the 02es-like subtype) are intrinsically red. This is evident by separating extrinsic from intrinsic color components. Since SALT2 is not trained on subluminous SNe~Ia, it compensates for their redder colors by inflating the $c$ parameter, thus extending the luminosity-width relation to negative values of x1. As expected, all subluminous SNe~Ia fall within the Cool region of the Branch et al. (2006) diagram, with the exception of 02es-like events, which show lower Si II 5972 pEW values. All subluminous subtypes tend to occur in more massive, redder host galaxies, and in the reddest local environments. Notably, 91bg- and 86G-like SNe~Ia explode at significantly larger normalized galactocentric distances. Finally, we identify the $pEW$ of the blended Ti II+Si II+Mg II absorption feature at 4300~A, along with s_BV, as robust and sufficient indicators for subclassifying subluminous SNe~Ia.

High-energy blazar light curves, in X-rays and beyond, have historically preferred a log-normal flux distribution, signifying multiplicative processes either in the jet itself or due to connection(s) with accretion. Here we present 18 year archival Fermi-LAT light curves (0.1-100 GeV) of the flat spectrum radio quasar (FSRQ) CTA 102 from August 2008 to November 2025, which underwent a huge flare in 2017, with a $\sim$ factor of 100 jump in $\gamma$-ray flux, along with similar flaring in X-rays. Our statistical analyses confirm that neither the pre nor the post-flare total GeV light curves follow a strictly log-normal distribution. Instead, we observe a statistically significant reduction in skewness from the pre to the post-flare light curves, which implies the blazar transitioned from an energetic state with frequent flaring to a more plateaued state with occasional flaring. We further find that this state transition can be explained through magnetic relaxation, where many reconnection events caused the 2017 flare, after which the magnetic field was ordered and its energy reached a minimum. To explain this further, we use a Monte Carlo simulation of a modified minijets-in-a-jet model where GeV flares are produced only when a maximum number of minijets move toward the broad line region and towards the line of sight, in the context of an external Compton model. The flux distributions (both observed and simulated) could be fit by a modified log-normal power-law distribution, implying our minijets model can reproduce the GeV flares in CTA 102 as well as their flux distributions.

We construct a gravitational open extension of the effective field theory of inflation in the Schwinger-Keldysh framework. While physical symmetries allow many open operators in the Schwinger-Keldysh action, most of them overconstrain the equations of motion, yielding inconsistent dynamics. We identify the minimal open operators compatible with propagating scalar and tensor modes and build the gravitational action, recovering dissipative models of inflation.

In this work, we carry out a new spectral reanalysis of NuSTAR and Suzaku observations of the disk reflection spectra in the stellar-mass black hole X-ray binary Cyg~X-1. We compare three types of models: a broken power-law disk emissivity profile with no assumption about the coronal shape used in the previous work of the same observations, a compact lamppost corona, and an extended disk-like corona motivated by recent X-ray polarization results. Our goal is to measure the systematic uncertainties caused by the assumed geometry, with a focus on key parameters such as the black hole spin and the inclination of the inner accretion disk. We find that the disk-like corona gives a fit that is statistically similar to the broken power-law and lamppost models, but it leads to more physically reasonable results, such as a lower inclination angle of about $30^{\circ}$. By using a variable disk density model, we measure the disk density to be $n_{\rm e}\approx10^{20}$\,cm$^{-3}$, which is similar to earlier results. While the extended corona model infers a wider allowed parameter space for black hole spin and the inner radius of the disk-shaped coronal region, this reflects the additional physical freedom of the model. Even so, the disk-like corona remains a strong and physically well-motivated candidate for explaining the X-ray emission from Cyg~X-1.

In this work, we investigate how the choice of initial vector potential and plasma parameters influences the development of accretion columns and jet formation in magnetized accretion flows. Using general relativistic magnetohydrodynamic simulations, we explore two different configurations of the vector potential $A_{\phi}$ and three plasma beta values $\beta$ = 50, 100, 500. We analyze how variations in the poloidal magnetic field strength and plasma magnetization affect magnetic flux accumulation near the black hole and the subsequent growth of the accretion column. Our results highlight the dependence of jet launching efficiency and accretion dynamics on the initial magnetic field topology and plasma beta, offering insight into the conditions that favor magnetically arrested disk or standard and normal evolution states.

The disc instability model successfully reproduces many of the observed properties of cataclysmic variables. However, additional ingredients such as mass-transfer variations, disc irradiation, stream-disc overflow, or inner-disc truncation must be included to explain certain systems. The physics underlying these processes is often poorly constrained, and our lack of knowledge is typically absorbed into extra free parameters, much like the $\alpha$-prescription for viscosity. In this paper, I examine how each of these ingredients affects the predicted light curves and discuss the limitations that arise from the growing number of unconstrained parameters on the model's predictive power.

We revisit the Starobinsky model of inflation in light of recent data from the Atacama Cosmology Telescope (ACT), which indicates a potential preference for a slightly larger scalar spectral index $n_s$ than predicted by the standard $R^2$ scenario. We demonstrate that a natural one-parameter generalization to a quadratic model $\sim L+L^2$ in the Lovelock invariant $L=R+\frac{\alpha}{4}{\cal G}$ ($\cal G$ is the Gauss–Bonnet term), can effectively resolve this minor tension. Scalar-tensor formulation of this theory yields an Einstein-frame Starobinsky-type scalar potential augmented by Gauss–Bonnet and derivative couplings, which modify the inflationary slow-roll dynamics. We show that a non-zero coupling $\alpha$ for the Gauss-Bonnet term can shift $(n_s, r)$ along a trajectory that brings the predictions into better agreement with the ACT likelihood. We also find that $L+L^2$ gravity, in its scalar-tensor formulation, is equivalent to Higgs inflation coupled to the Gauss–Bonnet term, and belongs to the Horndeski/galileon class of modified gravities. This work establishes the quadratic $f(L)$ gravity as a compelling and physically motivated extension that preserves the successes of Starobinsky inflation while improving its fit to modern precision cosmological data.

Standard perturbative calculations of scalar-induced gravitational waves (SIGWs) have neglected nonperturbative effects in the large-amplitude regime. We develop a hybrid numerical framework to signify nonperturbative effects on the stochastic gravitational wave (GW) background sourced by primordial curvature perturbations, focusing on the acoustic channel (fluid motions). Fully general-relativistic, spherically symmetric simulations are used to extract nonperturbative sound-shell profiles from isolated curvature peaks; these profiles are then embedded into three-dimensional lattice evolutions of relativistic hydrodynamics coupled to transverse-traceless metric perturbations to compute the acoustic GW spectra. The acoustic signal has a peak frequency determined by the comoving shell thickness, and its amplitude is extremely sensitive to the mean comoving separation of peaks, scaling approximately as $R_{*c}^{-7}$. We find a robust causal low-frequency tail $\propto k^{3}$, and the nonlinear hydrodynamic interactions can enhance the ultraviolet power. Comparing with SIGWs computed perturbatively from the same real-space configuration, we show that acoustic GWs can be amplified by an order of magnitude and display a peak shifted to a lower frequency in the large-amplitude regime. These results highlight the importance of nonperturbative effects for accurate predictions of stochastic GW signals induced from primordial curvature perturbations.

We present a comprehensive study of the 2021 outburst of GX 339-4 using AstroSat observations in the hard-intermediate (HIMS) and soft-intermediate states (SIMS). Spectral and timing analyses across these states suggest that during the SIMS, unabsorbed flux (0.1-3 keV), inner disc temperature, and "apparent" inner disc radius do not change, suggesting the stability of the disc. In the SIMS, the photon index decreases from 2.1 to 1.7, indicating spectral hardening. The power density spectra (PDS) suggest the presence of quasi-periodic oscillations (QPOs) in the HIMS and SIMS. The QPO frequency evolves from 0.1 Hz to 0.2 Hz in the HIMS, and further to 5.7 Hz in the SIMS. We also observe a decrease in QPO frequency from 5.7 Hz to 4.5 Hz during the SIMS. We discuss the evolution of the QPO, fractional root mean square (rms) amplitude, and time-lag spectra. We discover that variations in disc normalization, disc temperature, and coronal heating rate can reproduce the observed rms and lag spectra with a time delay between them.

$\delta$ Scuti stars in binary or multiple systems serve as crucial probes for studying stellar pulsation and evolution. However, many such systems are not ideal for asteroseismology due to uncertainties in mass transfer with close companions and the challenges of dynamically measuring all components' physical properties. The triple system DG~Leo, comprising an inner binary and a distant $\delta$ Scuti star, is an ideal target due to its well-separated pulsator. By combining new TESS photometry with archival spectroscopy, our dynamical analysis shows that the system's three components share similar masses, radii, and luminosities within errors, occupying coincident Hertzsprung–Russell diagram positions, indicative of coeval evolution. By fitting seven observed $\delta$ Scuti frequencies through asteroseismic modeling with dynamically constrained theoretical grids, we simultaneously trace the pulsating star's evolution and constrain the triple system's evolutionary stage, with the derived fundamental parameters showing consistency with the dynamical solutions. Our analysis reveals that all three components of DG~Leo are in the post-main-sequence phase, with a system age of $0.7664^{+0.1402}_{-0.1258}$~Gyr. Additionally, the $\delta$ Scuti component shows multiple non-radial modes with significant mixed-character frequencies, providing precise constraints on its convective core extent ($R_{\mathrm{cz}}/R = 0.0562^{+0.0137}_{-0.0021}$).

The generation of helical magnetic fields and the associated chiral asymmetry via the chiral anomaly is a generic feature in pseudoscalar inflation. In the presence of a Chern–Simons coupling between the inflaton and a U(1) gauge field, the homogeneous evolution of the inflaton induces a tachyonic instability in one circular polarization of the gauge field, resulting in the production of helical magnetic fields. In this work, we show that, in the case of a gauged lepton flavor symmetry, U(1)$_{L_i-L_j}$, this mechanism can lead to the generation of a sizable lepton asymmetry. In a simple setup, however, the resulting lepton asymmetry is typically too small to have an observational consequences, even setting aside constraints from baryon overproduction via sphaleron processes, due to the backreaction of the produced gauge fields and fermions on the inflationary dynamics. We demonstrate that this limitation can be overcome by implementing a mechanism to suppress fermion production during inflation. As a result, a much larger lepton asymmetry can be generated from the subsequent decay of magnetic helicity. Remarkably, for the gauged U(1)$_{L_\mu-L_\tau}$ symmetry, the generated asymmetry can be sufficiently large to suppress the primordial helium abundance, as may be inferred from recent cosmic microwave background observations by ACT.

We investigate axion-photon conversion in stochastic magnetic fields, focusing on the evolution of the photon intensity and polarizations induced by conversion into axions. Assuming Gaussian magnetic fields characterized by the power spectra of their helical/non-helical components, we express the expectation values and variances of the photon intensity and linear/circular polarizations after conversion in terms of these spectra. We find nontrivial dependencies of these statistical quantities on the characteristic magnetic field correlation length, the propagation distance, and the axion mass. Moreover, we find that nontrivial polarizations emerge even if the photons are initially unpolarized, that the variances of these observables become suppressed in specific frequency regions, and that a peak structure arises in the expectation value of the circular polarization in the presence of statistically helical magnetic fields. We also point out consistency relations among these statistical quantities that hold independently of the specific forms of the magnetic field power spectra.

The evolutionary scenario of early-type nitrogen-sequence Wolf-Rayet (WNE) stars predicts a slowly rotating subclass that typically forms after the red supergiant (RSG) phase. Their slow rotation rates are attributed to stellar winds that remove angular momentum transferred outward during core contraction. We incorporate improved prescriptions for internal gravity waves and the magnetic Tayler instability into single massive star evolution models. Our simulations successfully produce slowly rotating WNE stars and determine optimal parameters for both mechanisms ($A \ge 10$ for internal gravity waves (IGWs), $\alpha = 0.01$ for revised Tayler instability (TSF)). The results demonstrate that the efficiency of angular momentum transfer in massive stars is significantly enhanced compared to low-mass stars, both processes can self-consistently explain the slow rotation of WNE stars, confirming their efficiency in angular momentum redistribution and providing crucial theoretical support for the existence of this predicted stellar population.

We present the quantum-mechanical calculations of the dimensionless sensitivity coefficients Q to small changes in the fundamental physical constant mu = m_e/m_p - the electron-to-proton mass ratio - for a number of low-frequency (1-50 GHz) transitions of the acetaldehyde (CH3CHO) molecule. The calculations show that Q varies in the range from 0.62 to 3.61. An example of the practical use of the CH3CHO and CH3OH lines tracing the same regions in three molecular clouds, located at large galactocentric distances (D_GC ~ 8 kpc) is considered. This results in a limit on the mu variations of Delta mu/mu = (0.1 +/- 0.4)*10^(-7) which is in line with previously obtained most stringent upper limits on changes in mu based on other molecules and methods. The limit obtained restricts hypothetical violations of the Einstein principle of the local position invariance at the level of 4*10^(-8) in the Galactic disk at large galactocentric distances.

Collisions between hot hydrogen atoms and CO$_2$ play a central role in energy transfer and atmospheric escape in CO$_2$-rich planetary atmospheres. We present quantum mechanical $j_z$-conserving coupled-states calculations of state-resolved cross sections for H/D–CO$_2$ collisions at energies up to 5~eV, benchmarked to within 7% of close-coupling results. Scattering is strongly forward-peaked, yielding momentum-transfer cross sections substantially smaller than commonly assumed: mass-scaling from O/C–CO$_2$ systems overestimates H–CO$_2$ total cross sections by factors of 30–45, while existing empirical fits underestimate the low-energy regime by up to $\sim$45%. Isotopic substitution (H/D) produces energy-dependent differences of up to 35% at $E<0.1$~eV, invalidating uniform scaling approaches for D/H fractionation. Maxwellian-averaged rate coefficients derived from our cross sections are significantly smaller than mass-scaled values, implying reduced H–CO$_2$ energy transfer efficiency. In atmospheric escape modelling, these revisions can shift Martian exobase altitudes by 10–20~km, leading to order-unity changes in thermal escape rates, and have implications for hydrogen loss in early CO$_2$-dominated planetary atmospheres. Our results provide essential quantum-mechanical inputs for revisiting atmospheric evolution scenarios on Mars, early Earth, and CO$_2$-rich exoplanets.

In this work, we present a comprehensive semiclassical analysis of black hole radiation in a spatially flat FRW Universe for two fundamental nonclassical states: the Squeezed Number State (SNS) and the Coherent Squeezed Number State (CSNS). Unlike thermally modified earlier studies, SNS and CSNS constitute fully non-thermal, number-state-dependent quantum configurations. By embedding these states within the framework of semiclassical theory of gravity, we derive state-resolved expressions for the Hawking temperature, entropy variation, and corresponding mass loss of an evaporating black hole. The influence of the squeezing parameter $\rho$ and number state parameter $n$ on Hawking emission is examined through a series of analytical results supported by twelve detailed plots. The analysis reveals that the Hawking temperature exhibits monotonic growth with increasing $\rho$ and $n$, thereby elevating the effective temperature experienced at the black hole horizon. The entropy variations $\Delta \mathbb{S}_{\mathrm{SNS}}$ and $\Delta \mathbb{S}_{\mathrm{CSNS}}$ show strong nonlinear enhancement, especially at moderate and large squeezing values. Overall, the study extends earlier thermal squeezed-state approaches to a fully number-state-resolved framework, highlighting the sensitivity of Hawking emission to nonclassical quantum configurations. These findings contribute a new perspective on gravitational particle creation in cosmological settings.

The Quick Look data products from the Very Large Array Sky Survey (VLASS) contain widespread imaging artifacts arising from the simplified imaging algorithm used in their production. The catalog of double radio sources associated with active galactic nuclei (DRAGNs) found in the VLASS first epoch Quick Look release using the DRAGNhunter algorithm suffers from contamination from these artifacts. These sources contain two or three individual components, each of which can be an artifact. We train random forest models to classify these DRAGNs based on the number of artifacts they contain, ranging from zero to three artifacts. We optimize our models and mitigate the class imbalance of our dataset with judicious training set selection, and the best of our models achieves a weighted F1 score of $97.01\%^{+1.12\%}_{-1.32\%}$. Using our classifications, we produce a catalog of VLASS DRAGNs from which an estimated 99.3% complete catalog of 97.7% artifact-free sources can be extracted.

We construct a generic X-ray spectral model for the reflection component from the clumpy torus and dusty gas in the polar region (polar dusty gas) in an active galactic nucleus (AGN), designated as Inclusive spectral energy distribution Model of Polar dust And Clumpy Torus for X-ray (IMPACTX). To calculate the spectra, we utilize the Monte-Carlo based, 3-dimensional radiative transfer code SKIRT. The adopted geometry is the same as that of the IMPACT model (Ogawa et al. in prep.), consisting of a clumpy torus defined by Nenkova et al. (2008) and paraboloid-shell polar dusty gas launched at the inner radius of the torus. We discuss the dependence of the X-ray spectrum on geometrical parameters in comparison with the conventional torus-only model (XCLUMPY). As an example, we have simultaneously applied the IMPACTX and IMPACT models to the X-ray and infrared spectra of the nearby Seyfert 2 galaxy NGC 4388, and have found a solution that can well reproduce both spectra. This demonstrates the importance of using both X-ray and infrared data to constrain the nuclear structure of an AGN.

GW241011 and GW241110 both exhibit extremely asymmetric masses, high primary spins, and significant spin-orbit misalignment, which challenge the formation of first-generation binary black hole mergers formed from stellar collapse. This implies that these two gravitational wave events might originate from the hierarchical merger mechanism, with at least one of the black holes being the remnant of a previous merger. Here we investigate the origin of hierarchical mergers for GW241011 and GW241110 using the merger entropy index which measures the efficiency of entropy transfer for binary black hole mergers in general relativity. We find that GW241011 is consistent with hierarchical mergers in dense star clusters. The origin of GW241110 remains under debate due to its large distribution uncertainty, which leads to method-dependent inference and should be taken into account when interpreting this event in terms of hierarchical mergers.

We performed an optical study of the magnetic period-bouncer candidate IL Leo. Long-term photometric analysis over $\approx 20$ years reveals multiple state transitions. Modelling the ultraviolet and optical spectral energy distribution refined the white dwarf parameters, yielding a mass of $M_\textrm{wd} = 0.74 \pm 0.05 M_{\odot}$ and an effective temperature of $T_\mathrm{eff} = 12700 \pm 360$ K. We analyzed phase-resolved spectroscopy obtained with the 6-m BTA telescope and the VLT during the low state. Orbital variability of the H$\alpha$ emission, inferred from dynamical spectra and Doppler tomograms, suggests that it originates in the accretion stream. Zeeman splitting gives a mean magnetic field of $B = 40.7 \pm 0.5$ MG. Modelling two sets of cyclotron spectra determined a low-state accretion rate of $\dot{M} = (2.5 - 4.1) \times 10^{-13}~M_{\odot}$ yr$^{-1}$ and a magnetic field of $B_\mathrm{m} \approx 41$ MG near magnetic pole.

Mid-infrared (mid-IR) observations are crucial for understanding galaxy evolution, tracing star formation, and active galactic nuclei (AGN) activity via dust emission. This work presents mid-IR galaxy luminosity functions (LFs) at $0.5 < z < 6$, derived from the JWST Systematic Mid-infrared Instrument Legacy Extragalactic Survey (SMILES) program. We combine 8 MIRI bands ($5-25$ $\mu$m) of SMILES and archival 23-band HST+JWST NIRCam photometry to construct an extensive catalog containing 2,813 galaxies with sub-$\mu$Jy level completeness in the mid-IR. We obtain monochromatic (in 5.6, 7.7, 10, 12.8, 15, 18, 21, and 25.5 $\mu$m), $L_{\rm IR}$, and AGN LFs, with a limiting luminosity down to $10^{9.5}$ $L_\odot$ at $z=0.5-1.0$, $\sim10^{10.5}$ $L_\odot$ at $z=2.0-4.0$, and to $\sim10^{11}$ $L_\odot$ at $z=4.0-6.0$. With the unprecedented sensitivity and resolution of JWST, here we better constrain the faint-end slope and its evolution of the mid-IR LFs, quantifying the dusty cosmic star formation and black hole accretion histories out to $z \sim 5$. These results provide essential insights to refine our understanding of the obscured star formation and galaxy-AGN co-evolution over cosmic time.

Star clusters have long been central to the study of stellar evolution due to their chemically and chronologically homogeneous populations. Asteroseismology, the analysis of stellar oscillations and pulsations, provides precise information about properties such as masses, radii, and ages of stars in the field. However, these stars lack calibration to an absolute scale, and so this project seeks to utilize the data from NASA's TESS mission and ESA's Gaia mission to identify additional cluster stars suitable for asteroseismic analysis and calibration. In this work we analyze 14 stars belonging to 3 well-populated clusters, 5 additional stars that are the only detected oscillators in their respective clusters, and 3 detected oscillators of unknown cluster membership. By significantly expanding the number of clusters with measured oscillating giants, this project increases the opportunity for cross-validation between classical stellar models and asteroseismic methods, allowing for improvements in both calibration techniques and age estimations across the galaxy.

The pair-instability supernova (PISN) mechanism predicts a mass gap in the black hole population, where no stellar-origin black holes are expected to form. However, several binary black hole (BBH) merger events exhibit component masses that appear to lie within the PISN mass gap. If a gravitational-wave (GW) source is lensed, neglecting lensing effects leads to an underestimation of the luminosity distance and hence the redshift, resulting in an overestimation of the source-frame masses and potentially placing them within the mass-gap region. In this work, we analyze two representative events, GW190521 and GW231123. We establish a direct mapping between the lensing magnification and the fraction of posterior samples with source-frame masses below the mass-gap boundary. Adopting a lower bound of $65\,M_{\odot}$, we find that the magnifications required for $90\%$ of the posterior samples to fall below this boundary are $\mu=444$ for GW231123 and $\mu=39$ for GW190521. At these magnifications, the corresponding source-frame masses of the primary black hole are $52^{+10}_{-10}\,M_{\odot}$ and $42^{+19}_{-14}\,M_{\odot}$, with lensed source redshifts of $2.3^{+0.8}_{-0.5}$ and $2.7^{+1.5}_{-1.2}$, respectively. These results provide a quantitative framework for assessing the lensing hypothesis as a possible explanation for BBH mergers observed within the PISN mass gap, and highlight that the extreme magnifications required for GW231123 may challenge the astrophysical plausibility of simple strong-lensing interpretations.

Determining the relationship between nanoflare energies and their delays is the key for understanding the physical mechanism of the events and the plasma response. Nanoflares analyzed in this study were generated self-consistently via prescribed photospheric motions in a 3D multi-strand simulation of a subset of active region magnetic flux. Energies and durations were quantified using three distinct methods. In this study, we investigated the correlation between nanoflare energies (E) and delays ($\tau_D$) using two non-parametric, rank-based statistical tests. Across all methods, results consistently show little to no correlation. This is further supported by the distribution of the exponent $\alpha$ in the assumed relation $E \propto \tau_D^\alpha$, which peaks near zero, and by broad delay distributions within fixed energy bins. These findings are irrespective of whether delays are correlated with the energy of the preceding or subsequent event. They also hold for a subset of high-energy nanoflares. The absence of correlation suggests that nanoflare onset is not solely determined by a critical value of magnetic stress and may involve triggering by other events, perhaps related to a locally complex topology.

Tidal deformabilities are one of the observable quantities characterizing neutron stars, which are strongly associated with the stellar compactness, the ratio of the stellar mass to the radius. In addition to the tidal deformability, the quasinormal modes excited in a neutron star are also an important property for extracting information about the neutron star interior, adopting gravitational wave asteroseismology. In this study, we especially focus on the magnetic tidal deformability, which acts on the gravitational waveform from a neutron star binary merger as a higher-order effect than the electric tidal deformability, and derive the universal relations expressing the quasinormal modes, such as the fundamental ($f$-), 1st pressure ($p_1$-), and 1st spacetime ($w_1$-) modes, as a function of the magnetic tidal deformability. The universal relations derived in this study exhibit accuracy more or less comparable to those of the electric tidal deformability.

The magneto-Rayleigh-Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2-D and 3-D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

The Laser Interferometer Space Antenna (LISA), adopted by ESA and scheduled for the second half of the next decade, will drive a new revolution in the rapidly growing field of gravitational-wave astronomy, by extending GW observations into the hiterto unexplored millihertz regime. One of the key source classes of LISA is merging massive black hole binaries in the 1e4-1e7 Msun mass range detectable to very high redshifts $z\sim 15$. MBHBs lighter than 1e6 Msun can potentially be identified during their inspiral weeks before coalescence, allowing for complementary electromagnetic observations before, during, and after the two MBHs merge. This white paper aims at defining the optimal strategy to maximize the number of detected EM counterparts of LISA MBHB events.

As of late 2025 there are about 70 exoplanets that meet the formal criterion of having equilibrium temperatures allowing the presence of liquid water and about 50 of them orbit M-stars, known for their strong chromospheric activity. Most of these stars are close to the Sun and the planet-to-star mass and luminosity ratios are advantageous, allowing for a more detailed follow-up than of planets orbiting hotter and more massive stars. Many more planets orbiting late-type stars are expected to be discovered by Gaia and PLATO in the following years. However, the lingering question remains whether the UV and X-ray emission, associated with the stellar activity, allows for complex life. A comprehensive study focused on properties of flaring exoplanet hosts and their activity, on a much larger scale than these few tens (soon to become hundreds) of stars with habitable planets is called for, to answer the question if such stars can harbor habitable planets. The proposed Wide Field Survey telescope is well suited for this study.

Despite the tremendous impact of Type Ia supernovae (SNe Ia) on the field of cosmology, their underlying physics are still poorly understood. Studies have found an intriguing correlation between standardized Type Ia supernova (SN Ia) luminosities and host galaxy masses, commonly referred to as the "mass step". SNe Ia in massive galaxies appear systematically brighter than in lower-mass, star-forming hosts after standardization. However, previous analyses use host galaxy mass estimates derived largely from optical data alone and assume parametric forms for host star formation histories (SFHs), both of which are known to misestimate galaxy stellar masses. In this work, we re-examine the mass-step relation with a sample of SN Ia host galaxies complete in broadband optical (3000 Angstrom to 1 micron) and near-infrared (up to 1.8 micron) data, and in some cases including mid-infrared observations (up to 4.5 micron), using Prospector to derive nonparametric SFHs. We find that while the masses for these galaxies have indeed been underestimated, the overall trend in SN Ia luminosity versus host mass remains largely unchanged. We also uncover an environmental metallicity-dependent trend, in which low-metallicity galaxies (Z < -1 Zsol) may drive much of the observed SN Ia luminosity-mass correlation, suggesting that the chemical composition of the SN environment plays a central role in shaping the standardized brightness of SNe Ia.

We perform a Bayesian comparison between thawing quintessence and a cosmological constant, incorporating theoretically motivated priors on the phenomenological Padé-w parameters used to model thawing dynamics. We find that thawing quintessence is consistently preferred over a cosmological constant when combining BAO data from DESI DR2 and CMB data from Planck+ACT with any of the major supernova compilations, including the recently updated DES-Dovekie sample. This preference is not sensitive to our choice of prior, but it is contingent on the inclusion of supernovae in the analysis. We comment on the consistency between various information criteria and Bayesian evidence ratios, finding that the Deviance Information Criterion (DIC) tracks the Bayesian evidence more reliably than either the Akaike Information Criterion (AIC) or the Bayesian Information Criterion (BIC). Finally, we use observational likelihoods to identify which thawing trajectories are compatible with the available data, independently of theoretical priors.