Papers for Wednesday, Jan 08 2025

Using a recent Chandra ACIS observation of the eclipsing cataclysmic variable (CV) V1460 Her (2MASS J16211735+4412541) with a fast-rotating ($P_{\rm spin} = 38.9$~s) white dwarf, we estimated a flux $F_{0.5-7 {\rm keV}}^{\rm abs} = (1.9\pm 0.5)\times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, a factor of $\sim7$ lower than found from previous Swift XRT observations $\approx8$ years ago, when the CV was quiescent. The drop in the flux suggests a corresponding drop in the accretion rate, and the resulting intrinsic luminosity $L_{0.5-7 {\rm keV}} \sim 1.6\times 10^{29}$ erg s$^{-1}$ places V1460 Her among the lowest-luminosity magnetic CVs known, in a state of very low accretion.

Wet extreme mass-ratio inspirals (wet EMRIs), which arise from stellar-mass black holes (sBHs) inspiral into supermassive black holes (SMBHs) within the gas-rich environments of Active Galactic Nuclei (AGN), are primary sources of gravitational waves (GWs) for space-borne detectors like LISA, TianQin, and Taiji. Unlike "dry EMRIs", which form through gravitational scattering in nuclear star clusters, wet EMRIs are naturally accompanied by interactions with accretion disks, offering rich multi-messenger science opportunities. They are distinct in generating transient electromagnetic (EM) signals, such as quasi-periodic eruptions (QPEs), which serve as valuable probes of accretion disk physics and SMBH environments. Their GW signals provide an unprecedented precision of the order of $O(10^{-4}\sim 10^{-6})$ in measuring SMBH mass and spin, enabling the calibration of traditional EM techniques and offering insights into jet formation models. Additionally, wet EMRIs serve as bright and dark sirens for cosmology, facilitating percent-level precision measurements of Hubble parameter through AGN host identification or statistical association. These systems hold immense potential for advancing our understanding of black hole dynamics, accretion physics, and cosmology.

The spectroscopic performance of an X-ray microcalorimeter is compromised at high count rates. In this study, we utilize the Resolve X-ray microcalorimeter onboard the XRISM satellite to examine the effects observed during high count rate measurements and propose modeling approaches to mitigate them. We specifically address the following instrumental effects that impact performance: CPU limit, pile-up, and untriggered electrical cross talk. Experimental data at high count rates were acquired during ground testing using the flight model instrument and a calibration X-ray source. In the experiment, data processing not limited by the performance of the onboard CPU was run in parallel, which cannot be done in orbit. This makes it possible to access the data degradation caused by limited CPU performance. We use these data to develop models that allow for a more accurate estimation of the aforementioned effects. To illustrate the application of these models in observation planning, we present a simulated observation of GX 13+1. Understanding and addressing these issues is crucial to enhancing the reliability and precision of X-ray spectroscopy in situations characterized by elevated count rates.

Modeling noise in gravitational-wave observatories is crucial for accurately inferring the properties of gravitational-wave sources. We introduce a transdimensional Bayesian approach to characterise the noise in ground-based gravitational-wave observatories using the Bayesian inference software $\texttt{Bilby}$. The algorithm models broadband noise with a combination of power laws; narrowband features with Lorentzians; and shapelets to capture any additional features in the data. We show that our noise model provides a significantly improved fit of the LIGO and Virgo noise amplitude spectral densities compared to currently available noise fits obtained with on-source data segments. We perform astrophysical inference on well-known events in the third Gravitational-Wave Transient Catalog using our noise model and observe shifts of up to $7\%$ in the $90\%$ boundaries of credible intervals for some parameters. We discuss plans to deploy this framework systematically for gravitational-wave inference along with possible areas of improvement.

TOPCAT and STILTS are related packages for desktop analysis of tabular data, presenting GUI and command-line interfaces respectively to much of the same functionality. This paper presents features in TOPCAT that facilitate use of STILTS.

Numerous decameter-sized asteroids have been observed impacting Earth as fireballs. These objects can have impact energies equivalent to hundreds of kilotons of TNT, posing a hazard if they impact populated areas. Previous estimates of meteoroid flux using fireball observations have shown an Earth impact rate for decameter-size objects of about once every $2$-$3$ years. In contrast, telescopic estimates of the near-Earth asteroid population predict the impact rate of such objects to be of order $20$-$40$ years, an order-of-magnitude difference. While the cause of this discrepancy remains unclear, tidal disruption of a larger near-Earth body has been proposed as an explanation for these excess decameter-sized impactors. The release in 2022 of previously classified United States Government (USG) satellite sensor data for fireball events has provided a wealth of new information on many of these impacts. Using this newly available USG sensor data, we present the first population-level study characterizing the orbital and dynamical properties of 14 decameter-sized Earth impactors detected by USG sensors since 1994, with a particular focus on searching for evidence of tidal disruption as the cause of the impact rate discrepancy. We find there is no evidence for recent ($\lesssim 10^4$ years) tidal disruption and weak evidence for longer-term tidal disruption in the decameter impactor population, but that the latter conclusion is limited by small number statistics. We also investigate the origins of both the impactor and near-Earth asteroid populations of decameter-sized objects from the main asteroid belt. We find that both populations generally originate from the same source regions: primarily from the $\nu_6$ secular resonance ($\sim70$%) with small contributions from the Hungaria group ($\sim20$%) and the 3:1 Jupiter mean-motion resonance ($\sim10$%).

I hypothesize a physical explanation for the "Little Red Dots" (LRDs) discovered by the James Webb Space Telescope (JWST). The first star formation in the universe occurs in dense clusters, some of which may undergo runaway collapse and form an intermediate mass black hole. This process would appear as a very dense stellar system, with recurring tidal disruption events (TDEs) as stellar material is accreted by the black hole. Such a system would be compact, UV-emitting, and exhibit broad H-alpha emission. If runaway collapse is the primary mechanism for forming massive black hole seeds, this process could be fairly common and explain the large volume densities of LRDs. In order to match the predicted number density of runaway collapse clusters, the tidal disruption rate must be on the order of 10^-4 per year. A top-heavy stellar initial mass function may be required to match observations without exceeding the predicted LambdaCDM mass function. The TDE LRD hypothesis can be verified with followup JWST observations looking for TDE-like variability.

Long Period radio Transients (LPTs) are a mysterious new class of radio transients pulsating on periods of minutes to hours. So far, eight LPTs have been discovered predominantly at low Galactic latitudes, yet their nature remains unknown. Here, I present the first phase-resolved optical spectroscopy of the 2.9-h LPT GLEAM-X J0704-37, acquired with the 10-m Keck I telescope. Radial velocity (RV) shifts of $189\pm 3 \textrm{km s}^{-1}$ of an M5-type star in a binary system are detected on a period nearly equal to the radio period. Weak H$\alpha$ emission is also present, with some of it possibly originating from outside of the M dwarf. Based on the RV amplitude, and assuming a typical M dwarf mass, the companion mass must be $M \geq 0.22 M_\odot$. Calibrating the spectra with space-based \textit{Gaia} photometry reveals that the system is nearly four times closer than previously reported, at $d \approx 400$ pc, suggesting that more systems could be nearby and amenable to optical characterization. The optical spectrum between 3500-10,000 Angstrom is well modeled by a binary comprised of a massive white dwarf (WD; $T_\textrm{eff}\approx$7,300 K, $M\approx0.8-1.0M_\odot$) and M dwarf ($T_\textrm{eff}\approx$3,000 K, $M\approx0.14M_\odot$). Radio pulses arrive when the WD is at nearly maximum blueshift and the M dwarf at nearly maximum redshift, in contrast to what has been reported in a similar LPT, ILT J1101+5521. GLEAM-X J0704-37 is now the second LPT with an orbital period nearly equal to the radio period, hinting at two classes of LPTs: ``long LPTs'' ($P\gtrsim$78 min) associated with WD + M dwarf binary orbits, and ``short LPTs'' ($P\lesssim$78 min) related to WD or neutron star spins. This work demonstrates that precise localization of LPTs, which enables optical follow-up, will be key in uncovering the mechanism(s) that power this new class of phenomenon.

Hydrogen-poor supernovae (SNe) of Type Ibc are explosions of massive stars that lost their hydrogen envelopes, typically due to interactions with a binary companion. We consider the case where the natal kick imparted to the neutron star (NS) or black hole (BH) remnant brings the compact object to a collision with a main-sequence companion, eventually leading to full tidal disruption of the companion. Subsequently, super-Eddington accretion onto the NS/BH launches a powerful, fast wind which collides with the SN ejecta and efficiently converts the kinetic energy of the wind into radiation. The radiation is reprocessed by the surrounding ejecta into a luminous ($\sim 10^{44}$ erg s$^{-1}$ at peak), days to months-long transient with optical peaks from $-19$ to $-21$ mag, comparable to (super)luminous Type Ibc SNe and fast blue optical transients (FBOTs) like AT2018cow. From a Monte-Carlo analysis we estimate the fraction of tidal disruptions following SNe in binaries to be $\sim 0.1$--$1$\%, roughly compatible with the event rates of these luminous SNe. At the broad-brush level, our model reproduces the multi-wavelength and spectral observations of FBOTs, and has the potential to explain peculiar features seen in some (super)luminous SNe which are difficult to reproduce by the conventional magnetar spindown mechanism, such as late-time hydrogen lines, bumpy light curves, and pre-peak excess.

Star-Planet Magnetic Interactions (SPMI) have been proposed as a mechanism for generating stellar hot-spots with energy outputs on the order of $10^{19-21}$ watts. This interaction is primarily believed to be mediated by Alfvén waves propagating towards the star. The stellar atmosphere dictates where and how much of this incoming energy can actually be deposited as heat. The stellar transition region separating the chromosphere from the corona of cool stars gives rise to a significant variation of the Alfvén speed over a short distance, and therefore a reflection of the Alfvén waves at the transition region is naturally expected. We aim to characterize the efficiency of energy transfer due to SPMI by quantifying a frequency dependent reflection of the wave energy at the stellar transition region and its transmission to the stellar chromosphere. Magnetohydrodynamic simulations are employed to model the frequency-dependent propagation of Alfvén waves through a realistic background stellar wind profile. The transmission efficiency as a function of the wave frequency is quantified, and further analysis is conducted to characterize the overall energy transfer efficiency of SPMI in several candidate systems where chromospheric hotspots have been tentatively detected. Low-frequency waves experience greater reflection compared to high-frequency waves, resulting in reduced energy transfer efficiency for lower frequencies. Conversely, the parametric decay instability of Alfvén waves substantially diminishes the energy transfer efficiency at higher frequencies. As a result, there exists a frequency range where energy transfer is most efficient. A significant fraction of the Alfvén wave energy is reflected at the stellar transition region, and in most realistic scenarios, the transmission efficiency to the chromosphere is found to be approximately 10 percent.

A significant amount of star formation occurs and has occurred in environments unlike the solar neighbourhood. The majority of stars formed closer to the peak of the cosmic star formation rate (z > 1.3) and a great deal of star formation presently occurs in the central molecular zone (CMZ) of the Galaxy. These environments are unified by the presence of a high interstellar radiation field (ISRF) and a high cosmic ray ionisation rate (CRIR). Numerical studies of stellar birth typically neglect this fact, and those that do not have thus far been limited in scope. In this work we present the first comprehensive analysis of hydrodynamical simulations of star formation in extreme environments where we have increased the ISRF and CRIR to values typical of the CMZ and starburst galaxies. We note changes in the fragmentation behaviour on both the core and stellar system scale, leading to top-heavy core and stellar system mass functions in high ISRF/CRIR clouds. Clouds fragment less on the core scale, producing fewer but more massive cores. Conversely, the cores fragment more intensely and produce richer clusters of stellar systems. We present a picture where high ISRF/CRIR clouds fragment less on the scale of cores and clumps, but more on the scale of stellar systems. The change in fragmentation behaviour subsequently changes the mass function of the stellar systems that form through enhanced accretion rates.

Recent JWST observations have revealed the prevalence of spiral structures at $z > 1$. Unlike in the local Universe, the origin and the consequence of spirals at this epoch remain unexplored. We use public JWST/NIRCam data from the COSMOS-Web survey to map spiral structures in eight massive ($> 10^{10.5}\,\rm M_{\odot}$) star-forming galaxies at $z_{\rm spec} \sim 1.5$. We present a method for systematically quantifying spiral arms at $z>1$, enabling direct measurements of flux distributions. Using rest-frame near-IR images, we construct morphological models accurately tracing spiral arms. We detect offsets ($\sim 0.2 - 0.8\,\rm kpc$) between the rest-frame optical and near-IR flux distributions across most arms. Drawing parallels to the local Universe, we conclude that these offsets reflect the presence of density waves. For nine out of eighteen arms, the offsets indicate spiral shocks triggered by density waves. Five arms have offsets in the opposite direction and are likely associated with tidal interactions. For the remaining cases with no detected offsets, we suggest that stochastic 'clumpy' star formation is the primary driver of their formation. In conclusion, we find a multi-faceted nature of spiral arms at $z > 1$, similar to that in the local Universe.

Star-forming clumps have been found to significantly influence the star formation of gas-rich $z>1$ galaxies. Using public data from JWST/NIRCam (COSMOS-Web) and ALMA (FMOS-COSMOS), we study a sample of 32 massive ($>10^{10.5}\,\rm M_{\odot}$) main-sequence galaxies at $z_{\rm spec}\sim1.5$ with $\sim0.3\,\rm kpc$ resolution. We create composite morphological models consisting of bulge, disk, and clumps to fully 'deconstruct' the galaxy images. With the resulting measurements of the flux and size of these components, we find the following: (I)The combined contribution of clumps is $1-30\%$ towards the net star formation rate (SFR) of the host while contributing $1-20\%$ to its stellar mass. The clumps show a correlation between their stellar mass and SFR, but have an increased specific-SFR (sSFR) relative to the star-forming main sequence, with offsets ranging from $0\lesssim\Delta\log\rm sSFR\lesssim 0.4$. They feature star formation surface densities of $10^{-2}-10^{2}\,\rm M_{\odot}/yr/kpc^{2}$, consistent with values observed in local star-forming and starburst galaxies. (II)The clumps span a large range of characteristic sizes ($r_{e}\sim0.1-1\,\rm kpc$) and stellar masses ($\sim 10^{8.0-9.5}\,\rm M_{\odot}$). We estimate a mass-size relation ($r_{e}\propto\rm M_{\star}^{\,0.52\pm0.07}$) along with a stellar mass function (slope, $\alpha=-1.85\pm 0.19$), both suggesting a hierarchical nature similar to that expected in star-forming regions in local galaxies. (III)Our measurements agree with the properties of stellar clumps in $z\gtrsim1$ lensed systems, bridging the gap between lensed and unlensed studies by detecting structures at sub-kpc scales.(IV)Clumps are found to be preferentially located along spiral features visible primarily in the residual rest-frame near-IR images. In conclusion, we present an observation-based, coherent picture of star-forming clumps at $z>1$.

Thorne-Żytkow Objects (TŻOs) have been predicted to form when a neutron star is engulfed by a diffuse, convective giant envelope. Accretion onto a neutron star at a rate that is larger than $10^{-4}\, M_\odot$ yr$^{-1}$ is expected to lead to significant emission of neutrinos of all flavors with energy of 1-100 MeV. Since the neutrino signal is expected to largely vary in time (from milliseconds to thousands of years), we outline detection strategies tailored to the signal duration. We find that neutrino detection from TŻOs up to the Small Magellanic Cloud is within the reach of current- and next-generation neutrino observatories, such as Super- and Hyper-Kamiokande, the IceCube Neutrino Observatory, and JUNO. Interestingly, if targeted searches for neutrinos from TŻO candidates (e.g. VX Sgr in our Galaxy as well as HV 2112 and HV 11417 in the Small Magellanic Cloud) should lead to positive results, neutrinos could positively identify the nature of such sources and their accretion rate. Our findings should serve as motivation for establishing dedicated searches for neutrino emission from TŻOs. This is especially timely since it is challenging to detect TŻOs via electromagnetic radiation unambiguously, and the TŻO gravitational wave signal could be probed with next-generation detectors for sources within our Galaxy only.

We present an analysis of Hubble Space Telescope (HST) and XMM-Newton data of the tidal disruption event (TDE) candidate and quasi-periodic eruption (QPE) source GSN 069. Using ultraviolet (UV) and optical images at HST resolution, we show that GSN 069's emission consists of a point source superimposed on a diffuse stellar component. The latter accounts for $\leq 5\%$ of the UV emission in the inner 0.5"$\times$0.5" region, while the luminosity of the former cannot be attributed to stars. Analyzing the 2014/2018 HST UV spectra, we show that to leading order the intrinsic spectral shape is $\nu L_{\nu} \propto \nu^{4/3}$, with $\sim 10\%$ far UV flux variability between epochs. The contemporaneous X-ray and UV spectra can be modeled self-consistently in a thin disk framework. At observed epochs, the disk had an outer radius ($R_{\rm out}$) of $\mathcal{O}(10^3 R_{\rm g})$, and shows both cooling and expansion over the four years. Incorporating relativistic effects via numerical ray tracing, we constrain the disk inclination angle ($i$) to be $30^\circ \lesssim i \lesssim 65^\circ $ and identify a narrow region of spin-inclination parameter space that describes the observations. These findings confirm that GSN 069 hosts a compact, viscously expanding accretion disk likely formed after a TDE. Implications for QPE models are: (i) No published disk instability model can explain the disk's stability in 2014 (no QPEs) and its instability in 2018 (QPEs present); (ii) While the disk size in 2018 allows for orbiter/disk interactions to produce QPEs, in 2014 the disk was already sufficiently extended, yet no QPEs were present. These findings pose challenges to existing QPE models.

X-ray quasi-periodic eruptions (QPEs) are a novel mode of variability in nearby galactic nuclei whose origin remains unknown. Their multi-wavelength properties are poorly constrained, as studies have focused almost entirely on the X-ray band. Here we report on time-resolved, coordinated Hubble Space Telescope far ultraviolet and XMM-Newton X-ray observations of the shortest period X-ray QPE source currently known, eRO-QPE2. We detect a bright UV point source ($L_{\rm FUV} \approx {\rm few} \times 10^{41}$ erg s$^{-1}$) that does not show statistically significant variability between the X-ray eruption and quiescent phases. This emission is unlikely to be powered by a young stellar population in a nuclear stellar cluster. The X-ray-to-UV spectral energy distribution can be described by a compact accretion disk ($R_{\rm out} = 343^{+202}_{-138} \ R_{\rm g}$). Such compact disks are incompatible with typical disks in active galactic nuclei, but form naturally following the tidal disruption of a star. Our results rule out models (for eRO-QPE2) invoking i) a classic AGN accretion disk and ii) no accretion disk at all. For orbiter models, the expected radius derived from the timing properties would naturally lead to disk-orbiter interactions for both quasi-spherical and eccentric trajectories. We infer a black hole mass of log$_{10}(M_{\rm BH}) = 5.9 \pm 0.3$ M$_{\odot}$ and Eddington ratio of 0.13$^{+0.18}_{-0.07}$; in combination with the compact outer radius this is inconsistent with existing disk instability models. After accounting for the quiescent disk emission, we constrain the ratio of X-ray to FUV luminosity of the eruption component to be $L_{\rm X} / L_{\rm FUV} > 16-85$ (depending on the intrinsic extinction).

We present multi-wavelength analysis of ZTF23abelseb (AT 2023sva), an optically discovered fast-fading ($\Delta m_r = 2.2$ mag in $\Delta t = 0.74 $ days), luminous ($M_r \sim -30.0$ mag) and red ($g-r = 0.50$ mag) transient at $z = 2.28$ with accompanying luminous radio emission. AT 2023sva does not possess a $\gamma$-ray burst (GRB) counterpart to an isotropic equivalent energy limit of $E_{\rm{\gamma, \, iso}} < 1.6 \times 10^{52}$ erg, determined through searching $\gamma$-ray satellite archives between the last non-detection and first detection, making it the sixth example of an optically-discovered afterglow with a redshift measurement and no detected GRB counterpart. We analyze AT 2023sva's optical, radio, and X-ray observations to characterize the source. From radio analyses, we find the clear presence of strong interstellar scintillation (ISS) 72 days after the initial explosion, allowing us to place constraints on the source's angular size and bulk Lorentz factor. When comparing the source sizes derived from ISS of orphan events to those of the classical GRB population, we find orphan events have statistically smaller source sizes. We also utilize Bayesian techniques to model the multi-wavelength afterglow. Within this framework, we find evidence that AT 2023sva possesses a shallow power-law structured jet viewed slightly off-axis ($\theta_{\rm{v}} = 0.07 \pm 0.02$) just outside of the jet's core opening angle ($\theta_{\rm{c}} = 0.06 \pm 0.02$). We determine this is likely the reason for the lack of a detected GRB counterpart, but also investigate other scenarios. AT 2023sva's evidence for possessing a structured jet stresses the importance of broadening orphan afterglow search strategies to a diverse range of GRB jet angular energy profiles, to maximize the return of future optical surveys.

We investigate the chemical evolution of the Milky Way disc exploring various schemes of recent (last several Gyr) star formation episodes, as reported in recent observational works. We use a semi-analytical model with parametrized radial migration and we introduce gaussian star formation episodes constrained by those recent observations. We find significant impact of the star formation episodes on several observables, like the local age-metallicity and [alpha/Fe] vs metallicity relations, as well as the local stellar metallicity distribution or the existence of young [alpha/Fe]-rich stars. Moreover, we show that the recently found "wiggly" behaviour of the disk abundance gradient with age can be interpreted in terms of either star formation or infall episodes.

The solar system planetary architecture has been proposed to be consistent with the terrestrial and giant planets forming from material rings at ~1 au and ~5 au, respectively. Here, we show that super-Earths and mini-Neptunes may share a similar formation pathway. In our simulations conducted with a disk alpha-viscosity of 4e-3, super-Earths accrete from rings of rocky material in the inner disk, growing predominantly via planetesimal accretion. Mini-Neptunes primarily originate from rings located beyond the water snowline, forming via pebble accretion. Our simulations broadly match the period-ratio distribution, the intra-system size uniformity, and the planet multiplicity distribution of exoplanets. The radius valley constrains the typical total mass available for rocky planet formation to be less than 3-6 Earth masses. Our results predict that planets at ~1 au in systems with close-in super-Earths and mini-Neptunes are predominantly water-rich. Though relatively uncommon, at ~1% level, such systems might also host rocky Earth-sized planets in the habitable zone that underwent late giant impacts, akin to the Moon-forming event.

Using multi-viewpoint observations from STEREO and SOHO during three solar cycles from 23 to 25, we study the magnetic flux rope (MFR) structures of coronal mass ejections (CMEs) near the Sun and magnetic clouds (MCs) at 1au. The study aims to investigate two phenomena: 1) the occurrence rate of CMEs near Hale sector boundaries (HBs) and 2) solar-cycle variation of MFR axial orientations in CMEs and MCs. Our preliminary results include: 1) the axes of MFRs in cycle 25 present a systematic northward orientation, which is the same as in cycle 23 but opposite to cycle 24; 2) the majority of the MFRs occurred near HBs (within 30 degrees) and some exceptional events occurred at non-HBs; 3) the axial fields in MCs present a similar north-south orientation, which changes from cycle to cycle. We discuss the implication of solar cycle variations of MFR axial orientations for space weather forecasts.

We report the results from a study of two massive ($M_{500c} > 6.0 \times 10^{14} M_{\odot}$) strong lensing clusters selected from the South Pole Telescope cluster survey for their high Einstein radius ($R_E > 40''$), SPT-CLJ2325$-$4111 and SPT-CLJ0049$-$2440. Ground-based and shallow HST imaging indicated extensive strong lensing evidence in these fields, with giant arcs spanning 18\arcsec\ and 31\arcsec, respectively, motivating further space-based imaging followup. Here, we present multiband HST imaging and ground-based Magellan spectroscopy of the fields, from which we compile detailed strong lensing models. The lens models of SPT-CL\,J2325$-$4111 and SPT-CL\,J0049$-$2440 were optimized using 9, and 8 secure multiple-imaged systems with a final image-plane rms of 0\farcs63 and 0\farcs73, respectively. From the lensing analysis, we measure the projected mass density within 500~kpc of $M(<500 ~{\rm kpc}) = 7.30\pm0.07 \times 10^{14}$$M_{\odot}$, and $M(<500 ~{\rm kpc})=7.12^{+0.16}_{-0.19}\times 10^{14}$ $M_{\odot}$ for these two clusters, and a sub-halos mass ratio of $0.12\pm{0.01}$ and $0.21^{+0.07}_{-0.05}$, respectively. Both clusters produce a large area with high magnification ($\mu\geq 3$) for a source at $z=9$, $A^{lens}_{| \mu | \geq 3 }=4.93^{+0.03}_{-0.04} arcmin^2$, and $A^{lens}_{| \mu | \geq 3 }=3.64^{+0.14}_{-0.10} arcmin^2$ respectively, placing them in the top tier of strong lensing clusters. We conclude that these clusters are spectacular sightlines for further observations that will reduce the systematic uncertainties due to cosmic variance. This paper provides the community with two additional well-calibrated cosmic telescopes, as strong as the Frontier Fields, suitable for studies of the highly magnified background Universe.

In this paper, we investigate the effects of varying bulk viscosity coefficients $\zeta(t)=\zeta_{0}+\zeta_{1}H$ on cosmic evolution within the framework of $f(T)$ teleparallel gravity. We focus on two cases: (i) $\zeta_{1} \neq0$ and (ii) $\zeta_{1} =0$, deriving the Hubble parameter $H$ as a function of redshift $z$ using a linear $f(T)$ model ($f(T) = \alpha T$ where $\alpha \neq 0$). Using the combined $H(z)+Pantheon^{+}+BAO$ dataset, we obtain observational constraints on model parameters. For Case I ($\zeta_1 \neq 0$), best-fit values are $H_0=60.0^{+2.0}_{-1.9}$ km/s/Mpc, $\alpha=1.01^{+0.10}_{-0.098}$, $\zeta_0=40.1^{+1.9}_{-2.0}$, and $\zeta_1=0.123^{+0.093}_{-0.088}$, while for Case II ($\zeta_1 = 0$), they are $H_0=67.5^{+1.3}_{-1.3}$ km/s/Mpc, $\alpha=0.94^{+0.14}_{-0.13}$, and $\zeta_0=34.7^{+2.0}_{-2.0}$. The analysis reveals a transition in the deceleration parameter, indicating a shift from deceleration to acceleration of the universe's expansion, with present-day values of $q_{0} \approx -0.49$ and $q_{0} \approx -0.32$ for the respective cases. The jerk parameter $j(z)$ and effective EoS for the cosmic viscous fluid also support the cosmic acceleration, with trajectories aligning with the quintessence scenario. These findings underscore the potential of our $f(T)$ model dominated by bulk viscous matter in explaining cosmic acceleration.

We present a novel method that enables us to estimate the acceleration of individual millisecond pulsars (MSPs) using only their spin period and its time derivative. For our binary MSP sample, we show that one can obtain an empirical calibration of the magnetic braking term that relies only on observed quantities. We find that such a model for magnetic braking is only valid for MSPs with small surface magnetic field strengths ($<3\times10^8$ G) and large characteristic ages ($>$ 5 Gyr). With this method we are able to effectively double the number of pulsars with line-of-sight acceleration measurements, from 28 to 54 sources. This expanded dataset leads to an updated measurement of the total density in the midplane, which we find to be $\rho_0$ = 0.086 $\pm$ 0.001 stat. $\pm$ 0.006 sys. M$_\odot$/pc$^3$, and an updated measurement of the local dark matter density, which we calculate to be $\rho_{0,\mathrm{DM}}$ = 0.014 $\pm$ 0.005 stat. $\pm$ 0.006 sys. M$_\odot$/pc$^3$ (0.53 $\pm$ 0.30 GeV/cm$^3$). This updated value for $\rho_{0,\mathrm{DM}}$ is in good agreement with literature values derived from kinematic estimates. We show that each new acceleration measurement improves the precision of $\rho_{0,\mathrm{DM}}$ by the same amount as roughly 10$^5$ stars. The pulsar accelerations are very asymmetric above and below the disk; we show that the shape and size of this asymmetry can be largely explained by the north-south asymmetry of disk star counts and the offset in the Milky Way disk and halo centers of mass due to the Large Magellanic Cloud.

We present 6 GHz e-MERLIN observations of 42 $z<0.2$ type 1 and type 2 mostly radio-quiet quasars ($L_{\rm[OIII]}\gtrsim10^{42}$ erg s$^{-1}$; $L_{\rm AGN}\gtrsim10^{45}$ erg s$^{-1}$) from the Quasar Feedback Survey. The nature and origin of radio emission in these types of sources is typically ambiguous based on all-sky, low-resolution surveys. With e-MERLIN, we investigate radio emission on sub-kiloparsec scales ($\sim$10s-100s pc). We find 37/42 quasars are detected, with a diversity of radio morphologies, including compact cores, knots and extended jet-like structures, with sizes of 30-540 pc. Based on morphology and brightness temperature, we classify 76 per cent of the quasars as radio-AGN, compared to the $\sim$57 per cent identified as radio-AGN at the $\sim$1-60 kpc scales probed in prior studies. Combining results from e-MERLIN and the Very Large Array, 86 per cent reveal a radio-AGN. On average, $\sim$60 per cent of the total radio flux is resolved away in the e-MERLIN maps, and is likely dominated by jet-driven lobes and outflow-driven shocks. We find no significant differences in measured radio properties between type 1 and type 2 quasars, and estimate sub-relativistic jet speeds of $\sim$0.2-0.3c and modest jet powers of $P_\mathrm{jet} \approx \times$10$^{43}$ erg s$^{-1}$ for the few targets, where these measurements were possible. These quasars share characteristics with compact radio-selected populations, and the global radio emission likely traces strong interactions between the AGN (jets/outflows) and their host galaxy ISM from 10s parsec to 10s kiloparsec scales.

The Maunder Minimum (MM) was a period of prolonged solar activity minimum between 1645 and 1715. Several works have identified a significant number of problematic spotless days in the MM included in existing databases. We have found a list of exact spotless (in the second half of 1709) and spot days (January and August 1709) provided by Johann Heinrich Muller. We computed the most probable value and upper/lower limits of the active day fraction (ADF) from Muller's data using the hypergeometrical probability distribution. Our sample is not strictly random because Muller recorded observations in consecutive days when he observed sunspots. Therefore, our result represents an upper threshold of solar activity for 1709. We compared this result with annual values of the ADF calculated for the Dalton Minimum and the most recent solar cycles. We concluded that it was less active than most years both in the Dalton Minimum and in the most recent solar cycles. Therefore, the solar activity level estimated in this work for 1709 represents robust evidence of low solar activity levels in the MM.

We explore a scenario in which dark matter is a massive bosonic field, arising solely from quantum fluctuations generated during inflation. In this framework, dark matter exhibits primordial isocurvature perturbations with an amplitude of ${\cal O}(1)$ at small scales that are beyond the reach of current observations such as those from the CMB and large-scale structure. We derive an exact transfer function for the dark matter field perturbations during the radiation dominated era. Based on this result, we also derive approximate expressions of the transfer function in some limiting cases where we confirm that the exact transfer function reproduces known behaviors. Assuming a monochromatic initial power spectrum, we use the transfer function to identify the viable parameter space defined by the dark matter mass and the length scale of perturbations. A key prediction of this scenario is copious formation of subsolar mass dark matter halos at high redshifts. Observational confirmation of a large population of such low-mass halos will support for the hypothesis that dark matter originated purely from inflationary quantum fluctuations.

We study spatially-flat dynamical dark energy parametrizations, $w(z)$CDM, with redshift-dependent dark energy equation of state parameter $w(z)$ expressed using three different quadratic and other polynomial forms (as functions of $1-a$, where $a$ is the scale factor), without and with a varying cosmic microwave background (CMB) lensing consistency parameter $A_L$. We use Planck CMB anisotropy data (P18 and lensing) and a large, mutually-consistent non-CMB data compilation that includes Pantheon+ type Ia supernova, baryon acoustic oscillation (BAO), Hubble parameter ($H(z)$), and growth factor ($f\sigma_8$) measurements, but not recent DESI BAO data. The six $w(z)$CDM ($+A_L$) parametrizations show higher consistency between the CMB and non-CMB data constraints compared to the XCDM ($+A_L$) and $w_0 w_a$CDM ($+A_L$) cases. Constraints from the most-restrictive P18+lensing+non-CMB data compilation on the six $w(z)$CDM ($+A_L$) parametrizations indicate that dark energy dynamics is favored over a cosmological constant by $\gtrsim 2\sigma$ when $A_L = 1$, but only by $\gtrsim 1\sigma$ when $A_L$ is allowed to vary (and $A_L>1$ at $\sim2\sigma$ significance). Non-CMB data dominate the P18+lensing+non-CMB compilation at low $z$ and favor quintessence-like dark energy. At high $z$ P18+lensing data dominate, favoring phantom-like dark energy with significance from $1.5\sigma$ to $2.9 \sigma$ when $A_L = 1$, and from $1.1\sigma$ to $1.8\sigma$ when $A_L$ varies. These results suggest that the observed excess weak lensing smoothing of some of the Planck CMB anistropy multipoles is partially responsible for the $A_L = 1$ cases $\gtrsim 2\sigma$ evidence for dark energy dynamics over a cosmological constant.

We created new reddening maps and derived new extinction laws from visual to near-infrared passbands using improved RR~Lyrae period-absolute magnitude-metallicity relations, thus enabling distance estimates for individual bulge RR~Lyrae variables. The extinction law is most uniform in RIK and RJK and the distances to individual RR~Lyrae based on these colors are determined with an accuracy six and four percent, respectively. Using only the near-infrared passbands for distance estimation we inferred the distance to the Galactic center equal to djk = 8.2 +- 0.001(stat) +- 0.53(sys)pc after geometrical correction. We show that variations in the extinction law toward the Galactic bulge can mimic a barred spatial distribution in the bulge RR~Lyrae star population in visual passbands. This arises from a gradient in extinction differences along Galactic longitudes and latitudes, which can create the perception of the Galactic bar, particularly when using visual passband-based distances. A barred angle in the RR~Lyrae spatial distribution disappears when near-infrared passband-based distances are used, as well as when reddening law variations are incorporated in visual passband-based distances. The prominence of the bar, as traced by RR~Lyrae stars, depends on their metallicity, with metal-poor RR~Lyrae stars ([Fe/H]<-1.0dex) showing little to no tilt with respect to the bar. Metal-rich ([Fe/H]>-1.0dex) RR~Lyrae stars do show a barred/bulge signature in spatial properties derived using near-infrared distances, with an angle {\iota} = 18 +- 5deg, consistent with previous bar measurements from the literature. The 5D kinematic analysis, primarily based on transverse velocities, indicates a rotational lag in RR~Lyrae stars compared to red clump giants. Despite variations in the extinction law, our kinematic conclusions are robust across different distance estimation methods.

Pulsars are often lauded for their (relative) rotational and radio emission stability over long time scales. However, long-term observing programmes are identifying an increasing number that deviate from this preconceived notion. Using Gaussian process regression and Bayesian inference techniques, we investigated the emission and rotational stability of 259 radio pulsars that have been monitored using Murriyang, the Parkes 64 m radio telescope, over the past three decades. We found 238 pulsars display significant variability in their spin-down rates, 52 of which also exhibit changes in profile shape. Including 23 known state-switching pulsars, this represents the largest catalogue of variable pulsars identified to date and indicates that these behaviours are ubiquitous among the wider population. The intensity of spin-down fluctuations positively scales with increasing pulsar spin-down rate, with only a marginal dependence on spin-frequency. This may have substantial implications for ongoing searches for gravitational waves in the ensemble timing of millisecond pulsars. We also discuss challenges in explaining the physical origins of quasi-periodic and transient profile/spin-down variations detected among a subset of our pulsars.

Gravitational microlensing provides a unique opportunity to probe the mass distribution of stars, black holes, and other objects in the Milky Way. Population simulations are necessary to interpret results from microlensing surveys. The contribution from binary objects is often neglected or minimized in analysis of observations and simulations despite the high percentage of binary systems and microlensing's ability to probe binaries. To simulate the population effects we added multiple systems to Stellar Population Interface for Stellar Evolution and Atmospheres (SPISEA), which simulates stellar clusters. We then inject these multiples into Population Synthesis for Compact-object Lensing Events (PopSyCLE), which simulates Milky Way microlensing surveys. When making OGLE observational selection criteria, we find that 55% of observed microlensing events involve a binary system. Specifically, 14.5% of events have a multiple-lens and a single source, 31.7% have a single lens and a multiple-source, and 8.8% have a multiple-lens and a multiple-source. The majority of these events have photometric lightcurves that appear single and are fit well by a single-lens, single-source model. This suggests that binary source and binary lens-binary source models should be included more frequently in event analysis. The mean Einstein crossing time shifts from 19.1 days for single events only to 21.3 days for singles and multiple events, after cutting binary events with multiple peaks. The Einstein crossing time distribution of singles and single-peaked multiple events is better aligned with observed distributions from OGLE (arXiv:1707.07634) than singles alone, indicating that multiple systems are a significant missing piece between simulations and reality.

We present an independent validation and comprehensive re-calibration of S-PLUS Ultra-Short Survey (USS) DR1 12-band photometry using about 30,000--70,000 standard stars from the BEst STar (BEST) database. We identify spatial variation of zero-point offsets, up to 30--40\,mmag for blue filters ($u$, $J0378$, $J0395$) and 10\,mmag for others, predominantly due to the higher uncertainties of the technique employed in the original USS calibration. Moreover, we detect large- and medium-scale CCD position-dependent systematic errors, up to 50\,mmag, primarily caused by different aperture and flat-field corrections. We then re-calibrate the USS DR1 photometry by correcting the systematic shifts for each tile using second-order two-dimensional polynomial fitting combined with a numerical stellar flat-field correction method. The re-calibrated results from the XPSP and the SCR standards are consistent within 6\,mmag in the USS zero-points, demonstrating both the typical precision of re-calibrated USS photometry and a sixfold improvement in USS zero-point precision. Further validation using SDSS and Pan-STARRS1, as well as LAMOST DR10 and Gaia photometry, also confirms this precision for the re-calibrated USS photometry. Our results clearly demonstrate the capability and the efficiency of the BEST database in improving calibration precision to the milli-magnitude level for wide-field photometric surveys. The re-calibrated USS DR1 photometry is publicly available (\href{this https URL}{doi: https://doi.org/10.12149/101503}).

Contemporary real-time RFI mitigation is carried out at different stages primarily using regulatory and technical approaches. Regulatory approaches include spectrum management, radio quiet zones, and ensuring protection from self-generated RFI. The technical approaches include mitigation RFI in the analog and RF frontend systems, digital signal processing systems, and offline systems. As is known, the signal received by a radio telescope is a combination of the contributions from the astronomical signal and a combination of system and sky background noise. RFI has an additive effect on the signal received by the radio telescope. The key distinguishing properties of RFI are that it is generally stronger than the signal and non-Gaussian. The capability of signal processing receiver systems has grown manifold with the advent of high-speed signal processing platforms like Field Programmable Gate Arrays (FPGA)and Graphics Processing Unit (GPU). This has enabled the development of different signal-processing techniques for real-time RFI mitigation algorithms. This document provides an overview of contemporary techniques while focusing on the implementation of the same in specific radio telescopes.

Young planetary nebulae (PNe) are characterized by their hot central stars and the presence of abundant neutral and molecular components, which result from significant mass loss during the asymptotic giant branch (AGB) phase of stellar evolution. Far-UV \ion{He}{2}$\lambda$1025 line photons produced near the central star can undergo Raman scattering by hydrogen atoms, creating a broad emission feature centered at $\sim$ 6545~Å. We conducted high-resolution spectroscopy of 12 young PNe from April 2019 to March 2020 using the Bohyunsan Observatory Echelle Spectrograph ({\it BOES}). Building on the study by Choi and Lee, who identified Raman-scattered \ion{He}{2} at 6545~Å in NGC~6881 and NGC~6886, we report new detections of this feature in NGC~6741 and NGC~6884. Profile fitting reveals that the velocity of the \ion{H}{1} component relative to the \ion{He}{2} emission region ranges from $26-33~{\rm km~s^{-1}}$ in these PNe. Using photoionization modeling, we estimate the line flux of \ion{He}{2}$\lambda$1025 and derive Raman conversion efficiencies of 0.39, 0.21, 0.24, and 0.07 for NGC~6881, NGC~6741, NGC~6886, and NGC~6884, respectively. These results, combined with radiative transfer modeling, suggest the presence of \ion{H}{1} components with masses around $10^{-2}~M_\odot$, moving outward from the central \ion{He}{2} emission region at speeds characteristic of the slow stellar wind from a mass-losing giant star.

Solar flares release magnetic energy through reconnection, accelerating electrons into nonthermal velocity distributions, including crescent-shaped electron populations. These energetic electron distributions are crucial in driving instabilities which can lead to distinct electromagnetic emissions. This study investigates the emission properties of crescent-shaped electron velocity distribution functions (EVDFs) under different frequency ratios ($\omega_{pe}/\Omega_{ce}$), critical for understanding plasma conditions in various astrophysical environments, by comparing the emissions and intensities of waves among different cases. Here, we study and analyze three distinct frequency ratio conditions (2.2, 10, and 1, designated as cases A, B, and C, respectively). We found that the beam-Langmuir (BL) and upper-hybrid (UH) modes can be efficiently excited, leading to further plasma emissions in different cases. Our study reveals that the fundamental (O/F) emission can reach a maximum value of $\sim$$10^{-4} E_{\mathrm{k}0}$, while the harmonics (H) can extend to $\sim$$1.5 \times 10^{-5} E_{\mathrm{k}0}$ depending on the frequency ratio of the environment. The intensity of the fundamental mode exceeds previous findings for pure-ring, beam, and ring-beam distributions, highlighting the impact of crescent-shaped electron velocity distributions on wave excitation and emission processes. This effect is notably influenced by different frequency ratios, offering new insights into the way that nonthermal electron distributions affect the plasma emission process.

We present broadband X-ray observations of the High Mass X-ray Binary (HMXB) pulsar SMC X-2, using concurrent NuSTAR and NICER observations during its 2022 outburst. The source is found to be spinning with a period of 2.37281(3) s. We confirm the existence of cyclotron resonant scattering feature (CRSF) at 31 keV in addition to the iron emission line in the X-ray continuum of the source. Spectral analysis performed with the physical bulk and thermal Comptonization model indicates that the bulk Comptonization dominates the thermal Comptonization. Using phase-resolved spectroscopy, we have investigated the variations of the spectral parameters relative to pulse phase that may be due to the complex structure of magnetic field of the pulsar or the impact of the emission geometry. It is observed that the spectral parameters exhibit significant variabilities relative to the pulsed phase. Time-resolved spectroscopy is employed to examine the evolution of the continuum and changes in the spectral characteristics. Measurements of luminosity along with variations in cyclotron line energy and photon index suggest that the source may be accreting in the super-critical regime.

The existence of a strange quark star (QS) predicted in the Bodmer-Witten hypothesis has been a matter of debate. The combustion from a neutron star to a strange QS in its accreted process in a low-mass X-ray binary is proposed to be a scenario that generates gamma-ray bursts (GRBs); the baryon contamination of the outflow is very low and mainly from the masses of crusts ($M_{\rm crust}$) of QSs. A special subset of GRBs detected in the past 16 years are collected and used to estimate $M_{\rm crust}$ under this assumption of QSs as central engines. Correspondingly, $M_{\rm crust}$ is calculated in the frameworks of several models for cold dense quark matter (MIT bag model and Nambu-Jona-Lasino model with or without the impacts from the formation of color superconducting condensates being considered), for comparison with the observation. In conclusion, we find that the GRB samples have so far failed to provide positive support for this hypothesis.

In this work, we analyze the most recent short gamma-ray burst (sGRB) sample detected by the \emph{Fermi} satellite to reassess the sGRB luminosity function and formation rate. Using the empirical luminosity correlation, we first determine the pseudo redshifts of 478 sGRBs. Then, we use the maximum likelihood method to constrain the luminosity function and formation rate of sGRBs under various delay-time distribution models, finding the power-law delay model marginally preferred over the Gaussian and lognormal delay models based on the Akaike Information Criterion. The local formation rate of sGRBs is $1.73_{-0.45}^{+0.60}$ $\mathrm{Gpc^{-3}\,yr^{-1}}$, largely independent of the adopted delay-time distribution model. Additionally, we investigate the potential for joint detection of sGRBs and their gravitational wave (GW) counterparts from binary neutron star mergers using both current and future GRB and GW facilities. For sGRB detection, we consider three existing satellites: \emph{Fermi}, the Space-based multi-band astronomical Variable Objects Monitor (\emph{SVOM}), and the Einstein Probe (\emph{EP}). For GW detection, we examine two International GW Networks (IGWN): a four-detector network consisting of LIGO Hanford, Livingston, Virgo, and Kagra (IGWN4) and an upcoming five-detector network that includes these four detectors plus LIGO India (IGWN5). Our results indicate that for different delay-time distribution models, the joint sGRB and GW detection rates for \emph{Fermi}, \emph{SVOM}, and \emph{EP} with IGWN4 (IGWN5) lie within 0.09--0.31 $\mathrm{yr^{-1}}$ (0.55--1.98 $\mathrm{yr^{-1}}$), 0.03--0.11 $\mathrm{yr^{-1}}$ (0.26--0.80 $\mathrm{yr^{-1}}$), and 0.01--0.04 $\mathrm{yr^{-1}}$ (0.09--0.27 $\mathrm{yr^{-1}}$), respectively.

Theoretical predictions of element yields from the rapid neutron capture (r-) process are subject to large uncertainties due to incomplete knowledge of nuclear properties and approximative hydrodynamical modeling of matter ejection. A major source of uncertainty in determining ejecta composition and radioactive decay heat is the lack of nuclear mass data for exotic neutron-rich nuclei produced during neutron irradiation. We examine both model (systematic) and parameter (statistical) uncertainties affecting nuclear mass predictions and their impact on r-process nucleosynthesis, and consequently, the composition of neutron star merger ejecta. To estimate the effect of model uncertainties, we consider five nuclear mass models that accurately describe known masses. We also use a backward-forward Monte Carlo method to estimate uncorrelated uncertainties from local variations in model parameters, constraining them to experimentally known masses before propagating them to unknown masses of neutron-rich nuclei. These mass uncertainties are then applied to a 1.38-1.38 M$_{\odot}$ neutron star merger model, considering a wide range of ejecta trajectories. We find that uncorrelated parameter uncertainties lead to ejected abundance uncertainties of 20% up to A $\simeq$ 130, 40% between A=150 and 200, with peaks around A $\simeq$ 140 and A $\simeq$ 203, leading to deviations of 100-300%. While correlated model uncertainties generally exceed parameter uncertainties for most nuclei, both have a significant impact on heavy element production. Overall, improvements in nuclear models are essential to reducing uncertainties in r-process predictions. Both correlated model uncertainties and coherent determination of parameter uncertainties are crucial for sensitivity analysis in r-process nucleosynthesis.

In the context of the gravity-thermodynamics conjecture based on the Tsallis entropy-area relation, we investigate the standard $\Lambda$CDM cosmology and examine some potential deviations from it. Utilizing recent updates from geometrical datasets, including the DESI BAO measurements (2024), Planck CMB anisotropy measurements (2018), and the Pantheon+ catalogue for SNIa (2022), we conduct a thorough analysis via the window of Tsallis cosmology. In the first step, our analysis reveals no significant deviation from the standard model when using DESI BAO data alone, Pantheon+ data alone, or a combination of both. In the next step, we incorporate all datasets by adding the CMB data to our analysis, indicating a potential deviation from the standard model within the framework of Tsallis cosmology. Imposing the Planck prior to the sound horizon at the baryon drag epoch, we observe support for the standard model and consistency between our constraints on the Hubble constant and the Planck value. Finally, we compare the Tsallis and $\Lambda$CDM cosmologies using the Akaike Information Criterion (AIC).

Coronal holes are thought to be composed of relatively broad columnar structures known as plumes. Here we demonstrate that the plumes (and inter-plumes) in polar coronal holes are composed of fine-scale filamentary structure, with average scales of 2-10$^{\arcsec}$. The fine structure is the off-limb analogue of the previously found 'plumelets' of \cite{Uritsky_2021}. The off-limb observations enable an examination of the fine-structure without the influence of the underlying atmosphere along the line of sight. Hence, we show that the fine-scale structure is present at least until the edge of the field of view of the Solar Dynamics Observatory. The fine structure is found to have spatial distribution that follows a $k^{-1}$ power law perpendicular to the inferred magnetic field direction. For a small sample of the fine structure, the cross-sectional profiles are measured as a function of height. In some cases, the measurements indicate that the fine structure expands super-radially, consistent with existing models of polar field expansion and the expansion of the plumes. We discuss the implications of the presence of the fine structure with respect to understanding wave propagation in the coronal holes and their contribution to powering the solar wind.

Context. The Spectrometer/Telescope for Imaging X-Rays (STIX) onboard Solar Orbiter was designed to observe solar flares in the X-ray range of 4-150 keV, providing spectral, temporal and spatial information. Besides 30 imaging detectors, STIX has two additional detectors, the coarse flare locator (CFL) and the background (BKG) detector. Flares observed from Earth are classified using their peak X-ray flux observed by the GOES satellites. Roughly half of all flares observed by STIX are located on the backside of the Sun. These flares lack a GOES-class classification. Aims. In this paper, we describe the calibration of the BKG detector aperture sizes. Using the calibrated measurements of the BKG detector, we explore the relationship between the peak flux for flares jointly observed by STIX and GOES. This allows us to estimate the GOES flare classes of backside flares using STIX measurements. Methods. We looked at the 500 largest flares observed by both STIX and GOES in the time range Feb. 21 to Apr. 23. Aperture size calibration is done by comparing 4-10 keV counts of the BKG detector with the CFL measurements. In a second step, we correlate the calibrated STIX BKG peak flux with the GOES peak flux for individual flares. Results. We calibrated the BKG detector aperture sizes of STIX. Further, we showed that for the larger flares a close power law fit exists between the STIX BKG and GOES peak flux with a Pearson correlation coefficient of 0.97. This correlation provides a GOES proxy with a one sigma uncertainty of 11%. We were able to show that the BKG detector can reliably measure a broad range of GOES flare classes from roughly B5 up to at least X85 (assuming a radial distance of 1AU), making it an interesting detector-concept for future space weather missions. The largest flare observed by STIX to date is an estimated X16.5 $\pm$ 1.8 backside flare on the 20 Mai 2024.

It has been recently shown that cosmological models with scale-dependent primordial non-Gaussianities (sPNG) could provide a possible path to solve current cosmic tensions. Moreover, it has been pointed out that some of these models might mimic the effects of Warm Dark Matter (WDM) for several observables at low redshift. Here, we confirm the qualitative similarity of the matter power spectrum for sPNG and WDM models, but also point out differences in the halo mass function and void size function. We then jointly simulate WDM and sPNG together. Such simulations allow us to demonstrate that the joint impact of WDM and sPNG is close to the linear superposition of their respective effects at low redshift, at the percent level. We finally propose a model with mixed hot and cold dark matter together with sPNG, that reproduces the $\Lambda$CDM power spectrum at redshifts $z \leq 3$ but is still distinct in terms of halo statistics.

T Tauri stars and their discs are crucial for understanding stellar evolution and the formation of planets in low-mass systems. These stars exhibit significant variability, notably emitting intense X-ray flares due to magnetic reconnection events. During these events, magnetic energy is converted into kinetic energy of particles. Some of these particles then heat the plasma of the underlying chromosphere, emitting the observed X-rays. A portion of particles are thought to escape the chromosphere to interact with the surrounding circumstellar environment. The question is what impact do particles produced by magnetic reconnection events have on the discs of young stars. The complex characteristics of protoplanetary discs around T Tauri stars require an interdisciplinary strategy to enhance our understanding of these objects. This thesis contribute to establish a framework combining observational methodologies, chemical and dynamic models of protoplanetary discs, and the mechanics of acceleration and transport of energetic particles. We indroduce the role of ionisation in disc dynamics, including its sources. Not only from standard sources like stellar radiation and galactic cosmic rays but also from non-thermal ionisation due to magnetic reconnection events. We then examine energetic particles accelerated during magnetic reconnection events in T Tauri flares as an alternative source of ionisation, requiring the construction of a theoretical model based on solar flares. Then, we present a study revealing that particles from magnetic reconnection events could significantly contribute to the ionisation of the inner disc. Finally, we present a complementary study accounting for temporal effects, showing that considering these particles could increase the ionisation rate as well as viscosity, accretion rate, volumetric heating rate, and chemical complexity of inner protoplanetary discs.

Context. Recent developments in exoplanetary research highlight the importance of Love numbers in understanding their internal dynamics, formation, migration history and their potential habitability. Love numbers represent crucial parameters that gauge how exoplanets respond to external forces such as tidal interactions and rotational effects. By measuring these responses, we can gain insights into the internal structure, composition, and density distribution of exoplanets. The rate of apsidal precession of a planetary orbit is directly linked to the second-order fluid Love number, thus we can gain valuable insights into the mass distribution of the planet. Aims. In this context, we aim to re-determine the orbital parameters of WASP-43b-in particular, orbital period, eccentricity, and argument of the periastron-and its orbital evolution. We study the outcomes of the tidal interaction with the host star:whether tidal decay and periastron precession are occurring in the system. Method. We observed the system with HARPS, whose data we present for the first time, and we also analyse the newly acquired JWST full-phase light curve. We fit jointly archival and new radial velocity and transit and occultation mid-times, including tidal decay, periastron precession and long-term acceleration in the system. Results. We detected a tidal decay rate of \dotP_a=(-1.99pm0.50) and a periastron precession rate of \dotomega=(0.1851+0.0070-0.0077)=(0.1727+0.0083-0.0089)deg/d=(621.72+29.88-32.04)arcsec/d. This is the first time that both periastron precession and tidal decay are simultaneously detected in an exoplanetary system. The observed tidal interactions can neither be explained by the tidal contribution to apsidal motion of a non-aligned stellar or planetary rotation axis nor by assuming non-synchronous rotation for the planet, and a value for the planetary Love number cannot be derived. [...]

Context. Blazars are a distinct subclass of active galactic nuclei (AGN), known for their fast variability, high polarization, and intense emission across the electromagnetic spectrum, from radio waves to gamma rays. Gamma-ray blazar candidates of uncertain type (BCU) are an ongoing challenge in gamma-ray astronomy due to difficulties in classification and redshift determination. Aims. This study continues an optical spectroscopic campaign aimed at identifying the characteristics of BCUs to improve classification and redshift estimates, particularly focusing on low-synchrotron-peak sources. Methods. We conducted a detailed analysis of optical spectroscopic data for a sample of 21 low-synchrotron-peak BCUs plus one bl lac with contradictory results in the literature, using the 3.58-m Telescopio Nazionale Galileo (TNG, La Palma, Spain). Results. Our analysis identifies 14 out of the 21 BCUs as flat-spectrum radio quasars (FSRQs), demonstrating the effectiveness of our selection criteria. Notably, four FSRQs have redshifts exceeding 1, including 4FGL J2000.0+4214 at z = 2.04. Six sources are classified as bl lacs, with one of them, 4FGL J0746.5-0719, showing a featureless spectrum in this work despite previously exhibiting strong lines, suggesting it may be a changing-look blazar. One source remains classified as a BCU due to a noisy spectrum. Additionally, we observed a bl lac object, 4FGL J1054.5+2211, due to inconsistent redshift estimates in the literature, but we could not confirm any redshift due to its featureless spectrum. Our findings provide insights into the classification and redshift estimation of blazar candidates, emphasizing the need for continued spectroscopic monitoring.

Context. The detection of the He I 10830 A triplet in exoplanet atmospheres has opened a new window for probing planetary properties, including atmospheric escape. Unlike Lyman alpha, the triplet is less affected by ISM absorption. Sufficient XUV stellar irradiation may trigger the formation of the He I triplet via photoionization and posterior recombination processes in the planet atmospheres. Only a weak trend between stellar XUV and the planetary He I strength has been observed so far. Aims. We aim to confirm this mechanism for producing the He I absorption in exoplanetary atmospheres by examining a sample of planetary systems. Methods. We obtained homogeneous measurements of the planetary He I line EW and consistently computed the stellar XUV ionizing irradiation. We first derived new coronal models for the planet-host stars. We used updated data from the X-exoplanets database, archival X-ray spectra of M-type stars (including AU Mic and Proxima Cen), and new XMM-Newton X-ray data obtained for the CARMENES project. These data were complemented at longer wavelengths with publicly available HST, FUSE, and EUVE spectra. A total of 75 stars are carefully analyzed to obtain a new calibration between X-ray and EUV emission. Results. Two distinct relationships between stellar X-ray emission (5-100 A) and EUV_H (100-920 A) or EUV_He (100-504 A) radiation are obtained to scale the emission from late-type stellar coronae. A total of 48 systems with reported planetary He I 10830 A studies, exhibit a robust relationship between the planetary He I feature and the ionizing XUV_He received by the planet, corrected by stellar and planetary radii, and the planet's gravitational potential. Some outliers could be explained by a different atmospheric composition or the lack of planetary gaseous atmospheres. This relation may be used to predict the He I 10830 A absorption in exoplanet atmospheres.

Context. The Sun is a privileged laboratory of stellar evolution, thanks to the quality and complementary nature of available constraints. Using these observations, we are able to draw a detailed picture of its internal structure and dynamics which form the basis of the successes of solar modelling. Amongst such constraints, the depletion of lithium and beryllium are key tracers of the required efficiency and extent of macroscopic mixing just below the solar convective envelope. Thanks to revised determinations of these abundances, we may use them in conjunction with other existing spectroscopic and helioseismic constraints to study in detail the properties of macroscopic transport. Aims. We aim at constraining the efficiency of macroscopic transport at the base of the convective envelope and determining the compatibility of the observations with a suggested candidate linked with the transport of angular momentum in the solar radiative interior. Methods. We use recent spectroscopic observations of lithium and beryllium abundance and include them in solar evolutionary model calibrations. We test the agreement of such models in terms of position of the convective envelope, helium mass fraction in convective zone, sound speed profile inversions and neutrino fluxes. Results. We constrain the required efficiency and extent of the macroscopic mixing at the base of the solar convective envelope, finding that a power law of density with an index n between 3 and 6 would reproduce the data, with efficiencies at the base of the envelope of about 6000 cm2 /s, depending on the value of n. We also confirm that macroscopic mixing worsens the agreement with neutrino fluxes and that the current implementations of the magnetic Tayler instability are unable to explain the observations.

NGC 7419 is a young open cluster notable for hosting five Red Supergiants and a much higher abundance of Classical Be stars (CBe) than typical open clusters. We perform a membership analysis using GAIA DR3 data and machine learning techniques like Gaussian Mixture Models (GMM) and Random Forest (RF) and determine the cluster's mean distance to be 3.6 $\pm$ 0.7 kpc. We identify 499 GAIA-based members with a mass above $\sim$ 1.5 M$\odot$, and estimate the cluster's age to be $21.1 ^{+1.6}_{-0.6}$ Myr. Using our revised $H\alpha$ excess-based analysis, we find 42 CBe stars containing many known CBe stars, bringing the total number of CBe stars in NGC 7419 to 49 and the fraction of CBe to (B+CBe) members to 12.2 %. We investigate the variability of the candidate members from ZTF and NEOWISE data using standard deviation, median absolute deviation, and Stetson Index (J), and their periodicity using the Generalized Lomb Scargle Periodogram variability. We find that 66 % of CBe stars are variable: 23 % show periodic signals from pulsation or rotation, 41 \% exhibit variability from disk dynamics or binarity, and 14 % have long-term variations due to disk dissipation/formation. We also find that all pulsating CBe stars are early-type, while 50 \% of stars with long-term variations are early-type, and the other 50 % are mid-type. Our results agree with previous findings in the literature and confirm that CBe stars display variability through multiple mechanisms across different timescales.

In recent years, the interaction between dark matter (DM) and dark energy has become a topic of interest in cosmology. Interacting dark matter-dark energy (IDE) models have a substantial impact on the formation of cosmological large-scale structures, which serve as the background for DM halo evolution. This impact can be examined through the shape and spin orientation of halos in numerical simulations incorporating IDE effects. In our work, we use the N-body simulation pipeline ME-GADGET to simulate and study the halo spin and orientation in IDE models. We found that in models where DM transfers into DE (IDE I), the alignment of halo shapes with the surrounding tidal field is enhanced, while the alignment of halo spins with the tidal field is decreased compared to {\Lambda}CDM. Conversely, in models where DE transfers into DM (IDE II), the opposite occurs. We have provided fitted functions to describe these alignment signals. Our study provides the foundation for more accurate modeling of observations in the future such as China Space Station Telescope.

Interior models of gas giants in the Solar System traditionally assume a fully convective molecular hydrogen envelope. However, recent observations from the Juno mission suggest a possible depletion of alkali metals in Jupiter's molecular hydrogen envelope, indicating that a stable radiative layer could exist at the kilobar level. Recent studies propose that deep stable layers help reconcile various Jupiter observations, including its atmospheric water and CO abundances and the depth of its zonal winds. However, opacity tables used to infer stable layers are often outdated and incomplete, leaving the precise molecular hydrogen envelope composition required for a deep radiative zone uncertain. In this paper, we determine atmospheric compositions that can lead to the formation of a radiative zone at the kilobar level in Jupiter and Saturn today. We computed radiative opacity tables covering pressures up to $10^5$ bar, including the most abundant molecules present in the gas giants of the Solar System, as well as contributions from free electrons, metal hydrides, oxides, and atomic species, using the most up-to-date line lists published in the literature. These tables were used to calculate Rosseland-mean opacities for the molecular hydrogen envelopes of Jupiter and Saturn, which were then compared to the critical mean opacity required to maintain convection. We find that the presence of a radiative zone is controlled by the existence of K, Na, and NaH in the atmosphere of Jupiter and Saturn. For Jupiter, the elemental abundance of K and Na must be less than $\sim 10^{-3}$ times solar to form a radiative zone. In contrast, for Saturn, the required abundance for K and Na is below $\sim 10^{-4}$ times solar.

The presence and nature of low-frequency (0.1-10~mHz) Alfvénic waves in the corona has been established over the last decade, with many of these results coming from coronagraphic observations of the infrared Fe XIII line. The Cryo-NIRSP instrument situated at DKIST has recently begun acquiring science quality data of the same Fe XIII line, with at least a factor of 9 improvement in spatial resolution, a factor 30 increase in temporal resolution and an increase in signal-to-noise, when compared to the majority of previously available data. Here we present an analysis of 1~s cadence sit-and-stare data from Cryo-NIRSP, examining the Doppler velocity fluctuations associated with the Fe XIII 1074~nm coronal line. We are able to confirm previous results of Alfvénic waves in the corona as well as explore a new frequency regime. The data reveals that the power law behaviour of the Doppler velocity power spectrum extends to higher frequencies. This result appears to challenge some models of photospheric-driven Alfvénic waves that predict a lack of high frequency wave power in the corona due to strong chromospheric damping. Moreover, the high-frequency waves do not transport as much energy as their low-frequency counterparts, with less time-averaged energy per frequency interval. We are also able to confirm the incompressible nature of the fluctuations with little coherence between the line amplitude and Doppler velocity time-series.

It is known for long that the observed mass surface density of cored dark matter (DM) halos is approximately constant, independently of the galaxy mass (i.e., rhoc X rc simeq constant}, with rhoc and rc the central volume density and the radius of the core, respectively). Here we review the evidence supporting this empirical fact as well as its theoretical interpretation. It seems to be an emergent law resulting from the concentration-halo mass relation predicted by the current cosmological model, where the DM is made of collisionless cold DM particles (CDM). We argue that the prediction rhoc X rc simeq constant is not specific to this particular model of DM but holds for any other DM model (e.g., self-interacting) or process (e.g., stellar or AGN feedback) that redistributes the DM within halos conserving its CDM mass. In addition, the fact that rhoc X rc simeq constant is shown to allow the estimate of the core DM mass and baryon fraction from stellar photometry alone, particularly useful when the observationally-expensive conventional spectroscopic techniques are unfeasible.

We present a detailed analysis of the high-mass binary system V1216 Sco, an eclipsing Algol-type binary hosting a $\beta$ Cephei pulsator, with an orbital period of 3.92 days. This system was analyzed using TESS photometry and high-resolution spectroscopy from SALT HRS to investigate its orbital parameters, stellar properties, and evolutionary history. The TESS light curve, comprising over 12000 data points, revealed five independent pulsation frequencies within the beta Cephei range of 5 to 7 d$^{-1}$. Spectroscopic analysis provided radial velocities and disentangled atmospheric parameters, enabling precise orbital and evolutionary modeling. The system features a primary star of 11.72 M$_{\odot}$ and a secondary of 4.34 M$_{\odot}$. The secondary star has a radius near its Roche lobe, indicating recent or ongoing mass transfer. Evolutionary modeling with MESA-binary suggests that V1216 Sco underwent a case A mass transfer scenario, where mass transfer began while the donor was still on the main sequence. The system's evolutionary models indicate an age of 15 to 30 million years, highlighting the significant impact of binary interactions on stellar evolution. This study underscores the value of combining observational and theoretical approaches to understanding complex systems like V1216 Sco, emphasizing the role of mass transfer in shaping binary star evolution.

Fanaroff-Riley (FR) type 0 radio galaxies are a subclass of radio-loud active galactic nuclei (AGN) that lack extended kpc-scale jets, different from the classical FRI and FRII radio galaxies. They constitute the most abundant population of radio galaxies in the local Universe (z<0.1), yet remain largely unexplored. VLBI observations of a limited number of FR0s demonstrated that their central supermassive black hole (SMBH) are able to lunch mostly two-sided jets, with mildly relativistic bulk speed. In this work, we highlight the need of further high-resolution radio observations to probe the jet structures of these compact radio galaxies, by showing exploratory results of our EVN+eMERLIN observation campaign of FR0s. A preliminary analysis of these recent data reveals a possible change of the jet direction at different scales. We shortly discuss their physical conditions to explain the observed jet compactness, stressing the role of the SMBH spin vector in shaping their radio morphology

Solar flares stronger than X10 (S-flares, >X10) are the highest class flares which significantly impact on the Sun's evolution and space weather. Based on observations of Geostationary Orbiting Environmental Satellites (GOES) at soft X-ray (SXR) wavelength and the daily sunspot numbers (DSNs) since 1975, we obtained some interesting and heuristic conclusions: (1) Both S-flares and the more powerful extremely strong flares (ES-flares, >X14.3) mostly occur in the late phases of solar cycles and low-latitude regions on the solar disk; (2) Similar to X-class flares, the occurrence of S-flares in each solar cycle is somewhat random, but the occurrence of ES-flares seems to be dominated by the mean DSN (Vm) and its root-mean-square deviation during the valley phase (Vd) before the cycle: the ES-flare number is strongly correlated with Vd, and the occurrence time of the first ES-flare is anti-correlated with Vd and Vm. These facts indicate that the higher the Vm and Vd, the stronger the solar cycle, the more the ES-flares and the earlier they occurred. We proposed that the Sun may have a low-latitude active zone (LAZ), and most ES-flares are generated from the interaction between LAZ and the newly emerging active regions. The correlations and the linear regression functions may provide an useful method to predict the occurrence of ES-flares in an upcoming solar cycle, which derives that solar cycle 25 will have about 2 ES-flares after the spring of 2027.

Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multi-dimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localised flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20M_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for timescales typically considered by two-dimensional numerical simulations.

The black hole X-ray binary GRS 1915+105 was bright for 26 years since its discovery and is well-known for its disk instabilities, quasi-periodic oscillations, and disk wind signatures. We report a long-term spectral-timing tracing of this source from mid-2017 until the onset of the "obscured state", based on the complete data from the Neutron Star Interior Composition Explorer (NICER) and the Insight--Hard X-ray Modulation Telescope (HXMT), whose hard coverage decisively informs the modeling at lower energies. In the soft state predating 2018, we observed highly ionized winds. However, in the hard state shortly before transitioning into the "obscured state" on May 14, 2019 (MJD 58617), the winds exhibited a discernible reduction in ionization degree ($\log \xi$), decreasing from above 4 to approximately 3. Our analysis involves the measurement of the frequencies of the quasi-periodic oscillations and the estimation of the properties of the ionized winds and the intensities of different spectral components through spectroscopy during the decay phase. We delve into the origin of these infrequently observed warm outflows in the hard state. It is found that the launching radius of the winds in the hard decay phase is similar to that in the soft state, indicating the launching mechanism of those winds in both states is likely the same. The presence of the ionized winds is preferentially dependent on the periphery of the accretion disk, but not directly related to the corona activities in the center of the binary system.

We present the most extensive set to date of high-quality RR Lyrae light curve templates in the griz bands, based on time-series observations of the Dark Energy Camera Plane Survey (DECaPS) East field, located in the Galactic bulge at coordinates (RA, DEC)(J2000) = (18:03:34, -29:32:02), obtained with the Dark Energy Camera (DECam) on the 4-m Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). Our templates, which cover both fundamental-mode (RRab) and first-overtone (RRc) pulsators, can be especially useful when there is insufficient data for accurately calculating the average magnitudes and colors, hence distances, as well as to inform multi-band light curve classifiers, as will be required in the case of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). In this paper, we describe in detail the procedures that were adopted in producing these templates, including a novel approach to account for the presence of outliers in photometry. Our final sample comprises 136 RRab and 144 RRc templates, all of which are publicly available. Lastly, in this paper we study the inferred Fourier parameters and other light curve descriptors, including rise time,skewness, and kurtosis, as well as their correlations with the pulsation mode, period, and effective wavelength.

We present measurements and simulations of the polarization purity of leaky lens-antenna coupled microwave Kinetic Inductance Detectors (KIDs) at 1.5 THz. We find the integrated cross-polarization level to be at -21.5 dB for 1 f\#$\lambda$ spatial sampling. The measurements agree well with the theoretical description which is based on a combination of in-transmission simulation of the antenna feed, and an in-reception analysis of the antenna-KID system. Combined with the measured noise equivalent power of 5--7$\times$10$^{-20}$ W/$\sqrt{\mathrm{Hz}}$, these detectors are excellent candidates for large scale and high performance imaging polarimetric instruments.

Mrk 421 was in its most active state around early 2010, which led to the highest TeV gamma-ray flux ever recorded from any active galactic nuclei. We aim to characterize the multiwavelength behavior during this exceptional year for Mrk 421, and evaluate whether it is consistent with the picture derived with data from other less exceptional years. We investigated the period from November 5, 2009, (MJD 55140) until July 3, 2010, (MJD 55380) with extensive coverage from very-high-energy (VHE; E$\,>\,$100$\,$GeV) gamma rays to radio with MAGIC, VERITAS, Fermi-LAT, RXTE, Swift, GASP-WEBT, VLBA, and a variety of additional optical and radio telescopes. We investigated the variability and correlation behavior among different energy bands in great detail. We find the strongest variability in X-rays and VHE gamma rays, and PSDs compatible with power-law functions. We observe strong correlations between X-rays and VHE gamma rays. We also report a marginally significant positive correlation between high-energy (HE; E$\,>\,$100$\,$MeV) gamma rays and the ultraviolet band. We detected marginally significant correlations between the HE and VHE gamma rays, and between HE gamma rays and the X-ray, that disappear when the large flare in February 2010 is excluded from the correlation study. The activity of Mrk 421 also yielded the first ejection of features in the VLBA images of the jet of Mrk 421. Yet the large uncertainties in the ejection times of these radio features prevent us from firmly associating them to the specific flares recorded during the campaign. We also show that the collected multi-instrument data are consistent with a scenario where the emission is dominated by two regions, a compact and extended zone, which could be considered as a simplified implementation of an energy-stratified jet as suggested by recent IXPE observations.

We show that the behavior of the cosmic ray electron spectrum in the TeV energy band near the Earth is dominated by gluon condensation and anomalous electron/positron pair-production in Cygnus X.

We have conducted a planetary radial velocity measurement of the ultra-hot Jupiter WASP-121b using JWST NIRSpec phase curve data. Our analysis reveals the Doppler shift of the planetary spectral lines across the full orbit, which shifts considerably across the detector ($\sim$ 10 pixels). Using cross-correlation techniques, we have determined an overall planetary velocity amplitude of $K_{\rm p}=215.7\pm1.1$ km/s, which is in good agreement with the expected value. We have also calculated the dynamical mass for both components of the system by treating it as an eclipsing double-line spectroscopic binary, with WASP-121A having a mass of M$_{\star}$=1.330 $\pm$ 0.019 M$_{\odot}$, while WASP-121b has a mass of M$_{\rm p}$= 1.170 $\pm$ 0.043 M$_{\rm Jup}$. These dynamical measurements are $\sim3\times$ more precise than previous estimates and do not rely on any stellar modeling assumptions which have a $\sim$5\% systematic floor mass uncertainty. Additionally, we used stellar evolution modeling constrained with a stellar density and parallax measurement to determine a precise age for the system, found to be 1.11 $\pm$ 0.14 Gyr. Finally, we observed potential velocity differences between the two NIRSpec detectors, with NRS1 lower by 5.5$\pm$2.2 km/s. We suggest that differences can arise from day/night asymmetries in the thermal emission, which can lead to a sensitivity bias favoring the illuminated side of the planet, with planetary rotation and winds both acting to lower a measured $K_{\rm P}$. The planet's rotation can account for 1 km/s of the observed velocity difference, with 4.5$\pm$2.2 km/s potentially attributable to vertical differences in wind speeds.

In contemporary astronomy and astrophysics (A&A), the integration of high-performance computing (HPC), big data analytics, and artificial intelligence/machine learning (AI/ML) has become essential for advancing research across a wide range of scientific domains. These tools are playing an increasingly pivotal role in accelerating discoveries, simulating complex astrophysical phenomena, and analyzing vast amounts of observational data. For India to maintain and enhance its competitive edge in the global landscape of computational astrophysics and data science, it is crucial for the Indian A&A community to fully embrace these transformative technologies. Despite limited resources, the expanding Indian community has already made significant scientific contributions. However, to remain globally competitive in the coming years, it is vital to establish a robust national framework that provides researchers with reliable access to state-of-the-art computational resources. This system should involve the regular solicitation of computational proposals, which can be assessed by domain experts and HPC specialists, ensuring that high-impact research receives the necessary support. By building such a system, India can cultivate the talent, infrastructure, and collaborative environment necessary to foster world-class research in computational astrophysics and data science.

Context. The supernova remnant (SNR) W44 and its surroundings are a prime target for studying the acceleration of cosmic rays (CRs). Several previous studies established an extended gamma-ray emission that is set apart from the radio shell of W44. This emission is thought to originate from escaped high-energy CRs that interact with a surrounding dense molecular cloud complex. Aims. We present a detailed analysis of Fermi-LAT data with an emphasis on the spatial and spectral properties of W44 and its surroundings. We also report the results of the observations performed with the MAGIC telescopes of the northwestern region of W44. Finally, we present an interpretation model to explain the gamma-ray emission of the SNR and its surroundings. Methods. We first performed a detailed spatial analysis of 12 years of Fermi-LAT data at energies above 1 GeV, in order to exploit the better angular resolution, while we set a threshold of 100MeV for the spectral analysis. We performed a likelihood analysis of 174 hours of MAGIC data above 130 GeV using the spatial information obtained with Fermi-LAT. Results. The combined spectra of Fermi-LAT and MAGIC, extending from 100MeV to several TeV, were used to derive constraints on the escape of CRs. Using a time-dependent model to describe the particle acceleration and escape from the SNR, we show that the maximum energy of the accelerated particles has to be ' 40 GeV. However, our gamma-ray data suggest that a small number of lower-energy particles also needs to escape. We propose a novel model, the broken-shock scenario, to account for this effect and explain the gamma-ray emission.

Direct exo-Earth imaging is a key science goal for astronomy in the next decade. This ambitious task imposes a target contrast of ~10^-7 at wavelengths from I to J-band. In our prior study, we determined that polarization aberrations can limit the achievable contrast to 10^-5 to 10^-6 in the infrared. However, these results assumed a perfect coronagraph coupled to a telescope with an ideal coating on each of the mirrors. In this study we seek to understand the influence of polarization aberrations from segment-to-segment coating variations on coronagraphy and polarimetry. We use the Poke open-source polarization ray tracing package to compute the Jones pupil of each GSMT with spatially-varying coatings applied to the segments. The influence of the resultant polarization aberrations is simulated by propagating the Jones pupil through physical optics models of coronagraphs using HCIPy. After applying wavefront control from an ideal adaptive optics system, we determine that the segment-to-segment variations applied limit the performance of coronagraphy to a raw contrast of approximately 10^-8 in I-band, which is 2-3 orders of magnitude lower the target performance for high-contrast imaging systems on the ground. This is a negligible addition to the nominal polarization aberrations for ground-based systems. We further observe negligible degradation in polarimetric imaging of debris disks from segment-to-segment aberrations above and beyond the impact of nominal polarization aberration.

The ringed disk around HL Tau stands out as the iconic signature of planet formation, but the origin of the substructures is still debated. The HL Tau system also drives a powerful bipolar wind, and we analyze its outermost component traced by CO emission, to determine the relationship of the flow with the disk and its substructures. We use ALMA observations of the ${}^{12}$CO (2-1) line at 1.3 mm, with 0.2 km/s and ~ 0.28" resolution, conducted within the ALMA-DOT project. The channel maps and position-velocity diagrams show a rich structure of concatenated bubble- and arc-shaped features, whose size and distance from the source increase with velocity. The superposition of the features generates the apparent conical shape. The tomographic reconstruction of the morphology and kinematics of the red-shifted lobe suggests the presence of distinct nested shells having higher velocity and faster acceleration going toward the axis, rotating in the same sense of the disk. Such configuration can be justified by different classes of models. In this paper we compare the derived wind parameters with the predictions of magnetohydrodynamic (MHD) disk winds. Under this hypothesis, the launch radii of the three outermost shells are found to be at 55, 67 and 86 au from the star, coincident with regions of enhanced gas density in the disk. The wind may be capable of removing angular momentum from the outer disk, and we derive a magnetic lever arm of $\lambda \sim 4 - 5$, higher than that commonly adopted for MHD winds from these regions. Interpretations are discussed. The arrangement of the wind in nested shells with brighter emission rooted in rings of enhanced gas density could support the results of non-ideal MHD simulations according to which magnetic instabilities can generate the disk ring-gap system with a connected layered wind, alternatively to the action of yet undetected protoplanets.

We develop the framework of Linear Simulation-based Inference (LSBI), an application of simulation-based inference where the likelihood is approximated by a Gaussian linear function of its parameters. We obtain analytical expressions for the posterior distributions of hyper-parameters of the linear likelihood in terms of samples drawn from a simulator, for both uniform and conjugate priors. This method is applied sequentially to several toy-models and tested on emulated datasets for the Cosmic Microwave Background temperature power spectrum. We find that convergence is achieved after four or five rounds of $\mathcal{O}(10^4)$ simulations, which is competitive with state-of-the-art neural density estimation methods. Therefore, we demonstrate that it is possible to obtain significant information gain and generate posteriors that agree with the underlying parameters while maintaining explainability and intellectual oversight.

The mass accretion process controls pre-main-sequence evolution, although its intrinsic instability has yet to be fully understood, especially towards the protostellar stage. In this work, we have undertaken a thorough examination of the mid-infrared variability of Spitzer-selected YSOs in the Cygnus-X star-forming region over the last decade, using the NEOWISE time series. This work compares two groups of young stars: embedded Class I objects, and the more evolved flat-spectrum/Class II sources. We report on 48 candidate eruptive variables within these groups, including 14 with characteristics that resemble the photometric behaviour of FUors. We also include an additional 20 YSOs, which are of a less certain categorisation. We find the candidate FUors to be an order of magnitude more common among the younger Class I systems than more evolved objects. A large number of the identified short-duration eruptive YSOs display mid-infrared colour behaviour that is redder-when-brighter, which contrasts with optically bright outbursts seen in YSOs. Finally, we note the unusual long-term rising behaviours of four Class I YSOs, with rise timescales longer than five years, which is far slower than 6-12 month timescale for the majority of optically discovered FUors. Additionally, our broader investigation of MIR variability for embedded class I YSOs shows that there is a higher incidence of high amplitude variability for these stars, than is seen in class II sources. This holds true for all variable class I YSOs, not just the eruptive sources.

Using the one-dimensional numerical code MESA, we simulate mass accretion at very high rates onto massive main sequence stars, M=30, 60, 80 Mo, and find that these stars can accrete up to 10% of their mass without expanding much if we consider a simultaneous mass removal by jets. In this jetted-mass-removal accretion scenario, the accretion is through an accretion disk that launches jets. When the star expands due to rapid mass accretion, it engulfs the inner zones of the accretion disk and the jets it launches. We assume that these jets remove the outer layers of the envelope. We mimic this in the one-dimensional numerical code by alternating mass addition and mass removal parts. We add mass and energy, the accretion energy, to the outer layers of the envelope, leading to rapid stellar expansion. When the star expands by a few tens of percent, we stop mass addition and start mass removal until the star returns to its initial radius. We also show that the density of the accretion disk is larger than the density of the outer layers of the inflated envelope, allowing the disk to launch jets inside the outer inflated envelope layers. Our results show that main sequence stars can accrete mass at high rates while maintaining the deep potential well, as some models of eruptive systems require, e.g., some luminous red novae, the grazing envelope evolution, and the 1837-1856 Great Eruption of Eta Carinae.

Studying the isotopic composition of cosmic-rays (CRs) provides crucial insights into the galactic environment and helps improve existing propagation models. Special attention is given to the secondary-to-primary ratios of light isotopic components in CRs, as these measurements can offer complementary data compared to traditional secondary-to-primary ratios like B/C. Recently, a precision measurement of the Deuterium (D) abundance in CR in the 2-21 GV rigidity range provided by the AMS02 experiment unexpectedly detected an excess of D with respect to its expected secondary nature, opening the field for new measurements at high rigidity to determine how the spectrum evolves and whether there is confirmation of a primary or primary plus secondary origin. While there are theoretical models that attempt to explain this excess, the experimental uncertainties on D production cross-sections and on CR propagation models remain significant, and only new and precise measurements can dissipate existing doubts. In this work we review the current experimental scenario and we propose a dedicated experiment able to extend the D abundance measurement up to 100 GeV/nucl without the need of a magnetic spectrometer, using a multiple scattering based technique for the measurement of particle momentum. The expected performances of the proposed detector were assessed through a dedicated simulation using the GEANT4 package, and its role in the current particle physics scenario is discussed.

We investigate AGN feedback from an intermediate-mass black hole at the center of a dwarf spheroidal galaxy, by performing isolated galaxy simulations using a modified version of the GADGET-3 code. We consider Leo II (PGC 34176) in the Local Group as our simulation reference model. Beginning with black hole seeds ranging from $10^3$ to $10^6$ M$_{\odot}$, our simulations focus on comparing stellar-only feedback with AGN+stellar/SN feedback over 13.7 Gyr of galactic evolution. Our results indicate that a low-mass AGN in a dwarf galaxy influences the star formation history under specific physical conditions. While AGN feedback is generally negative on star formation, instances of positive feedback were also identified. Despite measurable effects on the evolution of the dwarf host galaxy, black hole seeds exhibited only marginal growth. We tested several physical scenarios as modified models in our simulations, primarily concerning the dynamics of the central black holes, which may wander within dwarf galaxies rather than being centrally located. However, none of these adjustments significantly impacted the growth of the black hole seeds. This suggests that intermediate-mass black holes may struggle to achieve higher masses in isolated environments, with mergers and interactions likely playing crucial roles in their growth. Nevertheless, AGN feedback exhibited non-negligible effects in our simulated dwarf spheroidal galaxies, despite the assumed dominant role of stellar feedback in the low-mass regime.

Over the past decade, observations of evaporating exoplanets have become increasingly common, driven by the discovery of the near-infrared helium-triplet line as a powerful probe of atmospheric escape. This process significantly influences the evolution of exoplanets, particularly those smaller than Jupiter. Both theoretical and observational studies have aimed to determine how efficiently exoplanets convert their host star's X-ray and ultraviolet (XUV) radiation into atmospheric mass loss. In this study, we employ the open-source atmospheric escape model p-winds to systematically analyze all publicly available helium triplet spectroscopic detections related to exoplanetary atmospheric escape. Our findings indicate that the retrieved outflows strongly depend on the ratio of XUV flux to planetary density ($F_{\text{XUV}}/\rho_p$), supporting the theoretical framework of energy-limited mass loss. We constrain population-level photoevaporative efficiencies to $0.34 \pm 0.13$ and $0.75 \pm 0.21$ for hydrogen-helium fractions of $0.90$ and $0.99$, respectively. These results offer new insights into exoplanetary atmospheric evolution and will aid future studies on exoplanet population demographics.

GUT phase transition during inflation solves the GUT monopole problem on one hand. On the other hand, gravitational waves (GWs) from such a phase transition, if it is first-order, are redshifted and deformed, and might be observed today in GW observatories. We review the formalism of inflated GWs and derive the general deformation function between inflated and uninflated GW spectra in the instant-source application. It is valid for any e-folding number of instant source. Applying the formalism to GUT phase transition, we find that the e-folding number at 15 or 25 can shift the GWs to 10 Hz or mHz hands, respectively, which might be tested in the future ground-based or space-based interferometers. We further generalise the discussion to inflated GWs via phase transition below the GUT scale. It is worth mentioning that, due to the deformation of the spectrum, the peak of inflated GWs is not simply a redshift of the peak of uninflated GWs.

This paper explicitly develops the three-field cosmological perturbation theory with a flat field space. We solve the background and perturbation equations numerically for three different cases. First, to check the consistency of the three-field formalism, we investigate an effective two-field model motivated by the two-block case of the multi-giant vacua matrix Inflation model. Then we investigate a completely three-field case without any direct interaction between different fields, and finally, a three-field case containing direct interactions. The power spectra of the curvature perturbations in all cases are computed numerically, and the effects of rapid turn in the power spectrum are highlighted.

Cosmic strings are powerful witnesses to cosmic events including any period of early matter domination. If such a period of matter domination was catalysed by metastable, long-lived particles, then there will be complementary signals to ascertain the nature of dark sector in experiments detecting primordial features in the gravitational wave (GW) power spectrum and laboratory searches for long-lived particles. We give explicit examples of global and local U(1) gauge extended dark sectors to demonstrate such a complementarity as the union of the two experiments reveals more information about the dark sector than either experiment. Demanding that Higgs-portal long-lived scalar be looked for, in various experiments such as DUNE, FASER, FASER-II, MATHUSLA, SHiP, we identify the parameter space which leads to complementary observables for GW detectors such as LISA and ET.

In some scenarios for the early universe, non-relativistic thermal dark matter chemically decouples from the thermal environment once the temperature drops below the Hubble rate. The value at which the energy density freezes out depends on the underlying model. In a simple setting, we provide a comprehensive study of heavy fermionic dark matter interacting with the light degrees of freedom of a dark thermal sector whose temperature $T$ decreases from an initial value close to the freeze-out temperature. Different temperatures imply different hierarchies of energy scales. By exploiting the methods of non-relativistic effective field theories at finite $T$, we systematically determine the thermal and in-vacuum interaction rates. In particular, we address the impact of the Debye mass on the observables and ultimately on the dark matter relic abundance. We numerically compare the corrections to the present energy density originating from the resummation of Debye mass effects with the corrections coming from a next-to-leading order treatment of the bath-particle interactions. We observe that the fixed-order calculation of the inelastic heavy-light scattering at high temperatures provides a larger dark matter depletion, and hence an undersized yield for given benchmark points in the parameter space, with respect to the calculation where Debye mass effects are resummed.

The transport of energetic particles is intimately related to the properties of plasma turbulence, a ubiquitous dynamic process that transfers energy across a broad range of spatial and temporal scales. However, the mechanisms governing the interactions between plasma turbulence and energetic particles remain incompletely understood. Here we present comprehensive observations from the upstream region of a quasi-perpendicular interplanetary (IP) shock on 2004 January 22, using data from four Cluster spacecraft to investigate the interplay between turbulence dynamics and energetic particle transport. Our observations reveal a transition in energetic proton fluxes from exponential to power-law decay with increasing distance from the IP shock. This result provides possible observational evidence of a shift in transport behavior from normal diffusion to superdiffusion. This transition correlates with an increase in the time ratio from $\tau_s/\tau_{c}<1$ to $\tau_s/\tau_{c}\gg1$, where $\tau_s$ is the proton isotropization time, and $\tau_{c}$ is the turbulence correlation time. Additionally, the frequency-wavenumber distributions of magnetic energy in the power-law decay zone indicate that energetic particles excite linear Alfvén-like harmonic waves through gyroresonance, thereby modulating the original turbulence structure. These findings provide valuable insights for future studies on the propagation and acceleration of energetic particles in turbulent astrophysical and space plasma systems.