Papers for Wednesday, Sep 25 2024

Relativistic outflows emanating from active galactic nuclei can extend up to kiloparsec scales in length, displaying a variety of complex morphologies. This study explores the intricate morphologies of such relativistic jets, mainly focusing on creating a bridge between magnetic instabilities in jets with observational signatures from complex radio galaxies. In particular, we aim to study the role of dynamical instabilities in forming distinctive morphological features by employing 3D relativistic magnetohydrodynamic (RMHD) simulations of rotating jets. Our simulations have further used the hybrid Eulerian-Lagrangian framework of the PLUTO code and generated the synthetic synchrotron emission and polarisation maps to compare with the observed signatures. Our analysis based on simulations of a continuously injected jet suggests that current-driven instabilities, notably the $|m|=1$ mode, generate rib-like structures that are seen in some of the recent radio galaxies using MeerKat, e.g. MysTail. In our contrasting simulations of the restarted jet, the kink-instability driven rib-like structures were formed relatively near the nozzle. In both cases, the jet dissipates its pre-existing magnetic energy through these instabilities, transitioning to a more kinetic energy dominant state. The turbulent structures resulting from this dissipation phase are filamentary and resemble the tethers as observed for the case of MysTail. This pilot study essentially provides a plausible qualitative explanation by bridging simulations of kink instability to produce synthetic radio features resembling the observed complex radio morphology of MysTail.

The recent development of statistical methods that can distinguish between stellar activity and dynamical signals in radial velocity (RV) observations has facilitated the discovery and characterization of planets orbiting young stars. One such technique, Gaussian process (GP) regression, has been regularly employed to improve the detection of a growing number of planets, but the impact of this model for mitigating stellar activity has not been uniformly analyzed for a large sample with real observations. The goal of this study is to investigate how GPs can affect the inferred parameters of RV-detected planets. We homogeneously analyze how two commonly adopted GP frameworks, a GP trained on RVs alone and a GP pretrained on photometry and then applied to RVs, can influence the inferred physical and orbital parameters compared to a traditional Keplerian orbit fit. Our sample comprises 17 short-period giant planets orbiting stars that exhibit a broad range of activity levels. We find that the decision to adopt GPs, as well as the choice of GP framework, can result in variations of inferred parameters such as minimum planet mass and eccentricity by up to 67% and 95%, respectively. This implies that the method for modeling stellar activity in RVs of young planet-hosting stars can have widespread ramifications on the interpretation of planet properties including their masses, densities, circularization timescales, and tidal quality factors. When mitigating stellar activity with GPs, we recommend carrying out comparative tests between different models to assess the sensitivity of planet physical and orbital parameters to these choices.

We report on the first self-consistent multidimensional particle-in-cell numerical simulations of non-homogeneous pair discharges in polar caps of rotation-powered pulsars. By introducing strong inhomogeneities in the initial plasma distribution in our simulations, we analyze the degree of self-consistently emerging coherence of discharges across magnetic field lines. In 2D, we study discharge evolution for a wide range of physical parameters and boundary conditions corresponding to both the absent and free escape of charged particles from the surface of a neutron star. We also present the results of the first 3D simulations of discharges in a polar cap with a distribution of the global magnetospheric current appropriate for a pulsar with $60^{\circ}$ inclination angle. For all parameters, we find the coherence scale of pair discharges across magnetic field lines to be of the order of gap height. We also demonstrate that the popular "spark" model of pair discharges is incompatible with the universally adopted force-free magnetosphere model: intermittent discharges fill the entire zone of the polar cap that allows pair cascades, leaving no space for discharge-free regions. Our findings disprove the key assumption of the spark model about the existence of isolated distinct discharge columns.

The clusters Abell 2061 and Abell 2067 in the Corona Borealis supercluster have been studied at different radio frequencies and are both known to host diffuse radio emission. The aim of this work is to investigate the radio emission in between them, suggested by low resolution observations. We analyse deep LOFAR HBA observations at 144 MHz to follow up on the possible intercluster filament suggested by previous 1.4 GHz observations. We investigate the radial profiles and the point-to-point surface brightness correlation of the emission in A2061 with radio and Xray observations, to describe the nature of the diffuse emission. We report the detection of diffuse radio emission on 800 kpc scale, more extended than previously known, reaching beyond the radio halo in A2061 towards A2067 and over the separation outside the two clusters R500 radii. We confirm the presence of a radio halo in A2061, while do not find evidence of diffuse emission in A2067. The surface brightness profile from the centre of A2061 shows an excess of emission with respect to the azimuthally averaged radio halo profile and X-ray background. We explore three different dynamical scenario to explain the nature of the diffuse emission. We analyse a trail of emission of 760 kpc between the radio halo and radio relic in A2061. This pre merger system closely resembles the two other cluster pairs where radio bridges connecting the radio halos on Mpc scales have been detected. The diffuse emission extends beyond each cluster R500 radius but in this unique case, the absence of the radio halo in A2067 is likely the reason for the observed 'gap' between the two systems. However, the point-to-point correlation results are challenging to explain. The classification of the emission remains unclear, and detailed spectral analysis and further Xray observations are required to understand the origin of the diffuse emission.

Wide, highly eccentric ($e>0.9$) compact binaries can naturally arise as progenitors of gravitational wave (GW) mergers. These systems are expected to have a significant population in the mHz band, with their GW signals characterized by ``repeated bursts" emitted upon each pericenter passage. In this study, we show that the detection of mHz GW signals from highly eccentric stellar mass binaries in the local universe can strongly constrain their orbital parameters. Specifically, it can achieve a relative measurement error of $\sim 10^{-6}$ for orbital frequency and $\sim 1\%$ for eccentricity (as $1-e$) in most of the detectable cases. On the other hand, the binary's mass ratio, distance, and intrinsic orbital inclination may be less precisely determined due to degeneracies in the GW waveform. We also perform mock LISA data analysis to evaluate the realistic detectability of highly eccentric compact binaries. Our results show that highly eccentric systems could be efficiently identified when multiple GW sources and stationary Gaussian instrumental noise are present in the detector output. This work highlights the potential of extracting the signal of ``bursting'' LISA sources to provide valuable insights into their orbital evolution, surrounding environment, and formation channels.

Dwarf spheroidal galaxies are known to be dominated by old stellar populations. This has led to the assumption that their gas-rich progenitors lost their gas during their infall in the Milky Way (MW) halo at distant look-back times. Here, we report a discovery of a tiny but robustly detected population of possibly young ($\sim$ 1 Gyr old) and intermediate-mass ($\rm 1.8 M_{\odot} \le M < 3 M_{\odot}$) stars in MW dwarf spheroidal galaxies. This was established on the basis of their positions in color-magnitude diagrams, after filtering out the bulk of the foreground MW using Gaia DR3 proper motions. We have considered the possibility that this population is made of evolved blue stragglers. For Sculptor, it seems unlikely, because 95.5% of its stars are older than 8 Gyr, leading to masses smaller than 0.9 M$_{\odot}$. This would only allow blue straggler masses of less than 1.8 M$_{\odot}$, which is much lower than what we observed. Alternatively, it would require the merger of three turnoff stars, which appears even more unlikely. On the other hand, the recent Gaia proper motion measurements of MW dwarf galaxies infer their low binding energies and large angular momenta, pointing to a more recent, $\le$ 3 Gyr, infall. Although the nature of the newly discovered stars still needs further confirmation, we find that they are consistent with the recent infall of the dwarf galaxies into the MW halo, when star formation occurred from the ram pressurization of their gas content before its removal by the hot Galactic corona. The abundance of this plausibly young population of stars is similar to the expectations drawn from hydrodynamical simulations. These results point to a novel origin for MW dwarf spheroidal galaxies.

We study the properties of the Luminous Red Galaxies (LRGs) selected from the 4th data release of the Kilo Degree Survey (KiDS-1000) via galaxy-galaxy lensing of the background galaxies from KiDS-1000. We use a halo model formalism to interpret our measurements and obtain estimates of the halo masses and the satellite fractions of the LRGs, resulting in halo masses $2.7 \times 10^{12} h^{-1} {\rm M}_{\odot}<M_{\rm h}< 2.6 \times 10^{13} h^{-1} {\rm M}_{\odot}$. We study the strength of intrinsic alignments (IA) using the position-shape correlations as a function of LRG luminosity, where we use a double power law to describe the relation between luminosity and halo mass to allow a comparison with previous work. Here, we directly link the observed IA of the (central) galaxy to the mass of the hosting halo, which is expected to be a fundamental quantity in establishing the alignment. We find that the dependence of the IA amplitude on halo mass is described well by a single power law, with amplitude $A = 5.74\pm{0.32}$ and slope $\beta_M = 0.44\pm{0.04}$, in the range $1.9 \times 10^{12}h^{-1} {\rm M}_{\odot}<M_{\rm h}<3.7 \times 10^{14} h^{-1} {\rm M}_{\odot}$. We also find that both red and blue galaxies from the source sample associated with the LRGs are oriented randomly with respect to the LRGs, although our detection significance is limited by the uncertainty in our photometric redshifts.

Recent observations by the Event Horizon Telescope (EHT) of supermassive black holes M87* and Sgr A* offer valuable insights into their spacetime properties and astrophysical conditions. Utilizing a library of model images (~2 million for Sgr A*) generated from general-relativistic magnetohydrodynamic (GRMHD) simulations, limited and coarse insights on key parameters such as black hole spin, magnetic flux, inclination angle, and electron temperature were gained. The image orientation and black hole mass estimates were obtained via a scoring and an approximate rescaling procedure. Lifting such approximations, probing the space of parameters continuously, and extending the parameter space of theoretical models is both desirable and computationally prohibitive with existing methods. To address this, we introduce a new Bayesian scheme that adaptively explores the parameter space of ray-traced, GRMHD models. The general relativistic radiative transfer code \ipole is integrated with the EHT parameter estimation tool THEMIS. The pipeline produces a ray-traced model image from GRMHD data, computes predictions for VLBI observables from the image for a specific VLBI array configuration and compares to data thereby sampling the likelihood surface via an MCMC scheme. At this stage we focus on four parameters: accretion rate, electron thermodynamics, inclination, and source position angle. Our scheme faithfully recovers parameters from simulated VLBI data and accommodates time-variabibility via an inflated error budget. We highlight the impact of intrinsic variability on model fitting approaches. This work facilitates more informed inferences from GRMHD simulations and enables expansion of the model parameter space in a statistically robust and computationally efficient manner.

Protoplanetary disks exhibit a rich variety of substructure in millimeter continuum emission, often attributed to unseen planets. As these planets carve gaps in the gas, dust particles can accumulate in the resulting pressure bumps, forming bright features in the dust continuum. We investigate the role of dust dynamics in the gap-opening process with 2D radiation hydrodynamics simulations of planet--disk interaction and a two-population dust component modeled as a pressureless fluid. We consider the opacity feedback and backreaction due to drag forces as mm grains accumulate in pressure bumps at different stages of dust growth. We find that dust dynamics can significantly affect the resulting substructure driven by the quasi-thermal-mass planet with $M_p/M_\star=10^{-4}$. Opacity feedback causes nonaxisymmetric features to become more compact in azimuth, whereas the drag-induced backreaction tends to dissolve nonaxisymmetries. For our fiducial model, this results in multiple concentric rings of dust rather than the expected vortices and corotating dust clumps found in models without dust feedback. A higher coagulation fraction disproportionately enhances the effect of dust opacity feedback, favoring the formation of crescents rather than rings. Our results suggest that turbulent diffusion is not always necessary to explain the rarity of observed nonaxisymmetric features, and that incorporating dust dynamics is vital for interpreting the observed substructure in protoplanetary disks. We also describe and test the implementation of the publicly-available dust fluid module in the PLUTO code.

We present a study of OVI-bearing gas around 247 low-mass (median log(M*/Msun)~8.7) galaxies at low redshifts (0.1 < z < 0.7) using background quasars as part of the MUSE Quasar-fields Blind Emitters Survey (MUSEQuBES). We find that the average OVI column density, ${\rm log}_{10}<N({\rm OVI})/{\rm cm}^{-2}>$ = $14.14^{+0.09}_{-0.10}$, measured within the virial radius for our sample, is significantly lower than for L_* galaxies. Combining 253 star-forming galaxies (mostly more massive) from the literature with 176 star-forming galaxies from MUSEQuBES, we find that both <N(OVI)> and the average covering fraction peak at log(M*/Msun)~9.5. The virial temperature corresponding to this stellar mass is ideal for OVI production via collisional ionization. However, we argue that photoionization and/or non-equilibrium processes are necessary to produce the OVI associated with low-mass, dwarf galaxies (log(M*/Msun)<9). The average OVI mass within the virial radius of dwarf galaxies is measured to be $10^{5.2_{-0.1}^{+0.1}}$ Msun. The characteristic normalized impact parameter at which the OVI covering fraction drops to half of its peak value is the largest (~1.1) for galaxies with stellar mass log(M*/Msun)~9.5. We report the presence of a highly ionized metal floor with ${\rm log}_{10}(N({\rm OVI})/{\rm cm}^{-2}) = 13.2$ outside the virial radius of dwarf galaxies inferred from median spectral stacking.

We present our findings from the first long X-ray observation of the hard X-ray source Swift J0826.2-7033 with XMM-Newton, which has shown characteristics of magnetic accretion. The system appears to have a long orbital period (~7.8 h) accompanied by short timescale variabilities, which we tentatively interpret as the spin and beat periods of an intermediate polar. These short- and long-timescale modulations are energy-independent, suggesting that photoelectric absorption does not play any role in producing the variabilities. If our suspected spin and beat periods are true, then Swift J0826.2-7033 accretes via disc-overflow with an equal fraction of accretion taking place via disc and stream. The XMM-Newton and Swift-BAT spectral analysis reveals that the post-shock region in Swift J0826.2-7033 has a multi-temperature structure with a maximum temperature of ~43 keV, which is absorbed by a material with an average equivalent hydrogen column density of ~1.6$\times$10$^{22}$ cm$^{-2}$ that partially covers ~27% of the X-ray source. The suprasolar abundances, with hints of an evolved donor, collectively make Swift J0826.2-7033 an interesting target, which likely underwent a thermal timescale mass transfer phase.

In this paper we analyse the spectro-photometric properties of the Type II supernova \sn, exploded at a distance of $19.9\,\rm{Mpc}$, in NGC~3206. Its early spectra are characterised by narrow high-ionisation emission lines, often interpreted as signatures of ongoing interaction between rapidly expanding ejecta and a confined dense circumstellar medium. However, we provide a model of the bolometric light curve of the transient that does not require sources of energy different than the H recombination and radioactive decays. Our model can reproduce the bolometric light curve of SN~2024bch adopting an ejected mass of $M_{bulk}\simeq5$\msun~surrounded by an extended envelope of only 0.2\msun~with an outer radius $R_{env}=7.0\times10^{13}\,\rm{cm}$. An accurate modelling focused on the radioactive part of the light curve, which accounts for incomplete $\gamma-$ray trapping, gives a $^{56}\rm{Ni}$ mass of 0.048\msun. We propose narrow lines to be powered by Bowen fluorescence induced by scattering of \ion{He}{II} Ly$\alpha$ photons, resulting in the emission of high-ionisation resonance lines. Simple light travel time calculations based on the maximum phase of the narrow emission lines place the inner radius of the H-rich, un-shocked shell at a radius $\simeq4.4\times10^{15}\,\rm{cm}$, compatible with an absence of ejecta-CSM interaction during the first weeks of evolution. Possible signatures of interaction appear only $\sim69\,\rm{days}$ after explosion, although the resulting conversion of kinetic energy into radiation does not seem to contribute significantly to the total luminosity of the transient.

We present a detailed study of the kinematics of OVI-bearing gas around 60 low-mass (median log(M*/Msun)~8.9) galaxies at low redshift (0.1 < z < 0.7) using background quasars (median impact parameter $\approx115$ kpc) as part of the MUSE Quasar-fields Blind Emitters Survey (MUSEQuBES). We find that the majority of the OVI absorbers detected within the virial radius have line-of-sight velocities smaller than the escape velocities and are thus consistent with being gravitationally bound, irrespective of the halo mass. However, the fraction of such absorbers declines at larger impact parameters. The Doppler $b$ parameter and the velocity width ($\Delta v_{90}$) of the OVI absorbers exhibit large scatter inside the virial radius of the host galaxies, but the scatter declines sharply at impact parameter $D \gtrsim 2R_{\rm vir}$. For high-mass galaxies (log(M*/Msun)>9), OVI absorption displays a larger kinematic spread, quantified by the pixel-velocity two-point correlation function (TPCF). However, this difference disappears for the isolated galaxies when the pixel velocities are scaled by the galaxy's circular velocity. We do not find any significant difference between the TPCF of isolated and group galaxies when the stellar mass is controlled for. A significant fraction of groups (4/6) with four or more member galaxies do not show any detectable OVI absorption, likely due to the passive nature of nearest galaxies.

We describe the spectrophotometric calibration of the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrometer (MRS) aboard the James Webb Space Telescope (JWST). This calibration is complicated by a time-dependent evolution in the effective throughput of the MRS; this evolution is strongest at long wavelengths, approximately a factor of 2 at 25um over the first two years of the mission. We model and correct for this evolution through regular observations of internal calibration lamps. Pixel flatfields are constructed from observations of the infrared-bright planetary nebula NGC 7027, and photometric aperture corrections from a combination of theoretical models and observations of bright standard stars. We tie the 5--18um flux calibration to high signal/noise (S/N; ~ 600-1000) observations of the O9 V star 10 Lacertae, scaled to the average calibration factor of nine other spectrophotometric standards. We calibrate the 18--28um spectral range using a combination of observations of main belt asteroid 515 Athalia and the circumstellar disk around young stellar object SAO 206462. The photometric repeatability is stable to better than 1% in the wavelength range 5--18um, and the S/N ratio of the delivered spectra is consistent between bootstrapped measurements, pipeline estimates, and theoretical predictions. The MRS point-source calibration agrees with that of the MIRI imager to within 1% from 7 to 21um and is approximately 1% fainter than prior Spitzer observations, while the extended source calibration agrees well with prior Cassini/CIRS and Voyager/IRIS observations.

The spin of intermediate-mass black holes (IMBHs) growing through repeated black hole mergers in stellar clusters statistically asymptotes to zero. Putative observations of IMBHs with dimensionless spin parameter $\chi\gtrsim 0.6$ would require a phase of coherent gas accretion to spin up the black hole. We estimate the amount of gas necessary to produce a given IMBH spin. If the observed IMBH mass and spin are $M\gtrsim 1000~M_\odot$ and $\chi\gtrsim 0.6$, respectively, the IMBH must have coherently accreted at least $\sim 100~M_\odot$ of gas. In this scenario, as long as the spin is not maximal, the IMBH can only accrete at most half of its mass. Our estimates can constrain the relative contribution of accretion and mergers to the growth of IMBHs in dense stellar environments.

The galaxy cluster Abell 496 has been extensively studied in the past for the clear sloshing motion of the hot gas on large scales, but the interplay between the central radio galaxy and the surrounding cluster atmosphere is mostly unexplored. We present a dedicated radio, X-ray, and optical study of Abell 496 aimed at investigating this connection. We use deep radio images obtained with the Giant Metrewave Radio Telescope at 150, 330 and 617 MHz, Very Large Array at 1.4 and 4.8 GHz, and VLA Low Band Ionosphere and Transient Experiment at 340 MHz, with angular resolutions ranging from 0.''5 to 25''. Additionally, we use archival Chandra and Very Large Telescope MUSE observations. The radio images reveal three distinct periods of jet activity: an ongoing episode on sub-kpc scales with an inverted radio spectrum; an older episode that produced lobes on scales $\sim$20 kpc which now have a steep spectral index ($\alpha = 2.0 \pm 0.1$); and an oldest episode that produced lobes on scales of $\sim$50 - 100 kpc with an ultra-steep spectrum ($\alpha = 2.7 \pm 0.2$). Archival Chandra X-ray observations show that the older and oldest episodes have excavated two generations of cavities in the hot gas of the cluster. The outermost X-ray cavity has a clear mushroom-head shape, likely caused by its buoyant rise in the cluster's potential. Cooling of the hot gas is ongoing in the innermost 20 kpc, where H$\alpha$-bright warm filaments are visible in VLT-MUSE data. The H$\alpha$-filaments are stretched towards the mushroom-head cavity, which may have stimulated ICM cooling in its wake. We conclude by commenting on the non-detection of a radio mini-halo in this vigorously sloshing, but low-mass, galaxy cluster.

The emergence of magnetic flux from the deep convection zone plays an important role in the solar magnetism, such as the generation of active regions and triggering of various eruptive phenomena, including jets, flares, and coronal mass ejections. To investigate the effects of magnetic twist on flux emergence, we performed numerical simulations of flux tube emergence using the radiative magnetohydrodynamic code R2D2, and conducted a systematic survey on the initial twist. Specifically, we varied the twist of the initial tube both positively and negatively from zero to twice the critical value for kink instability. As a result, regardless of the initial twist, the flux tube was lifted by the convective upflow and reached the photosphere to create sunspots. However, when the twist was too weak, the photospheric flux was quickly diffused and not retained long as coherent sunspots. The degree of magnetic twist measured in the photosphere conserved the original twist relatively well, and was comparable to actual solar observations. Even in the untwisted case, a finite amount of magnetic helicity was injected into the upper atmosphere because the background turbulence added helicity. However, when the initial twist exceeded the critical value for kink instability, the magnetic helicity normalized by the total magnetic flux was found to be unreasonably larger than the observations, indicating that the kink instability of the emerging flux tube may not be a likely scenario for the formation of flare-productive active regions.

The correction map method means extended phase-space algorithm with correction map. In our research, we have developed a correction map method, specifically the dissipated correction map method with trapezoidal rule, for numerical simulations of gravitational waves from spinning compact binary systems. This new correction map method, denoted as $CM3$, has shown remarkable performance in various simulation results, such as phase space distance, dissipated energy error, and gravitational waveform, closely resembling the high-order precision implicit Gaussian algorithm. When compared to the previously used midpoint map which denoted as $C_2$, the $CM3$ consistently exhibits a closer alignment with the highly accurate Gaussian algorithm in waveform evolution and orbital trajectory analysis. Through detailed comparisons and analyses, it is evident that $CM3$ outperforms other algorithms, including $CM2$ and $C_2$ mentioned in this paper, in terms of accuracy and precision in simulating spinning compact binary systems. The incorporation of the trapezoidal rule and the optimization with a scale factor $\gamma$ have significantly enhanced the performance of $CM3$, making it a promising method for future numerical simulations in astrophysics. With the groundbreaking detection of gravitational waves by the LIGO/VIRGO collaboration, interest in this research domain has soared. Our work contributes valuable insights for the application of matched filtering techniques in the analysis of gravitational wave signals, enhancing the precision and reliability of these detection.

Cosmic rays may be dynamically very important in driving large-scale galactic winds. Edge-on galaxies give us an outsider's view of the radio halo, which shows the presence of extra-planar cosmic-ray electrons and magnetic fields. We present a new radio continuum imaging study of the nearby edge-on galaxy NGC 4217 in order to study the distribution of extra-planar cosmic rays and magnetic fields. We both observe with the Jansky Very Large Array (JVLA) in the S-band (2-4 GHz) and with LOw Frequency ARray (LOFAR) at 144 MHz. We measure vertical intensity profiles and exponential scale heights. We re-image both JVLA and LOFAR data at matched angular resolution in order to measure radio spectral indices between 144 MHz and 3 GHz. Confusing point-like sources were subtracted prior to imaging. Intensity profiles are then fitted with cosmic-ray electron advection models, where we use an isothermal wind model that is driven by a combination of pressure from the hot gas and cosmic rays. We discover a large-scale radio halo on one (northwestern) side of the galactic disc. The morphology is reminiscent of a bubble extending up to 20 kpc away from the disc. We find spectral ageing in the bubble which allows us to measure advection speeds of the cosmic-ray electrons accelerating from 300 to 600 $\rm km\, s^{-1}$ . Assuming energy equipartition between the cosmic rays and the magnetic field, we estimate the bubble can be inflated by a modest 10 per cent of the kinetic energy injected by supernovae over its dynamical time-scale of 35 Myr. While no active galactic nucleus (AGN) has been detected, such activity in the recent past cannot be ruled out. Non-thermal bubbles with sizes of tens of kiloparsec may be a ubiquitous feature of star-forming galaxies showing the influence of feedback. To determine possible contributions by AGN feedback, will require deeper observations.

The relative contribution of emission from stellar sources and accretion onto supermassive black holes to reionization has been brought into focus again by the apparent high abundance of faint AGN at $4\lesssim z\lesssim11$ uncovered by JWST. We investigate the contribution of these faint AGN to hydrogen and the early stages of helium reionization using the GPU-based radiative transfer code ATON-HE by post-processing a cosmological hydrodynamical simulation from the SHERWOOD-RELICS suite of simulations. We study four models of reionization: two previously studied galaxy-only late-end reionization models and two new models -- a QSO-assisted model and a QSO-only model. In the QSO-assisted model, 1\% of the haloes host AGN with a 10~Myr lifetime, and the AGN luminosities are scaled such that the AGN contribution to the hydrogen-ionizing emissivity is 20\% of that contributed by galaxies. In the QSO-only model, quasars account for all the hydrogen-ionizing emissivity, with 10\% of the haloes hosting AGN, each with a 10 Myr lifetime. All models are calibrated to the observed mean Lyman-$\alpha$ forest transmission at $5\lesssim z\lesssim6.2$. We find that the QSO-assisted model requires an emissivity factor of $1.8$ lower than the galaxy-only models towards the end of reionization and fits the observed distribution of the Lyman-$\alpha$ optical depths well. Our QSO-only model is inconsistent with the observed Lyman-$\alpha$ optical depths distribution. It also results in too high IGM temperatures at $z\lesssim 5$ due to an early onset of HeII reionization unless the escape fraction of HeII-ionizing photons is assumed to be low. Our results suggest that a modest contribution to reionization by faint AGN is in good agreement with the Lyman-$\alpha$ forest data. In contrast, a scenario dominated by faint AGN appears difficult to reconcile with these observations.

We scrutinize the reported lensing anomaly of the CMB by considering several phenomenological modifications of the lensing consistency parameter, $A_{\rm L}$. Considering Planck spectra alone, we find statisically significant evidence for scale dependence (`running') of $A_{\rm L}$. We then demonstrate that the anomaly is entirely driven by Planck's low multipoles, $\ell \leq 30$. When these data points are excluded a joint analysis with several other datasets clearly favors $\Lambda$CDM over the extended $\Lambda \rm CDM+A_L$ model. Not only that the lensing anomaly and low $\ell$ anomaly of the CMB go away in this case, but also the $S_8$ tension is ameliorated, and only the Hubble tension persists.

Interaction between galaxies play a pivotal role in their evolution. Ongoing star formation in spiral galaxies can be affected by these processes. We select a sample of interacting galaxies in field environments at various interaction stages and are nearly face-on and chose galaxy pairs NGC 2207/IC 2163, NGC 4017/4016 (ARP 305) and NGC 7753/7752 (ARP 86). We use the UltraViolet Imaging Telescope (UVIT) onboard AstroSat to characterize the star-forming regions in the galaxy with a superior resolution of ~1.4". We identified and characterized star-forming regions in the UVIT images of the sample and correlated them with the neutral hydrogen (HI) distribution. We detected localized regions of enhancement in star formation surface density and distortions in the sample of galaxies. We found this consistent with the distribution of HI in the galaxy. These are possible evidence of past and ongoing interactions affecting the star formation properties in the galaxies. We then conducted a study to understand whether the observed enhancements hold true for a wider sample of interacting galaxies. We observe a moderate enhancement in the star formation rate (SFR) with the interaction class, with the maximum of 1.8 being in the merger class of galaxies. We studied the SFR enhancement for the main galaxies in our sample as a function of pair mass ratio and pair separation. We observe a strong anti-correlation between the SFR enhancement and pair mass ratio and no linear correlation between the enhancement and pair separation. This suggests that the enhancement in interaction-induced star formation may be more strongly influenced by the pair mass ratios, rather than the pair separation. We also infer that the pair separation can possibly act as a limiting parameter for the SFR enhancement.

This paper presents VARnet, a capable signal processing model for rapid astronomical timeseries analysis. VARnet leverages wavelet decomposition, a novel method of Fourier feature extraction via the Finite-Embedding Fourier Transform (FEFT), and deep learning to detect faint signals in light curves, utilizing the strengths of modern GPUs to achieve sub-millisecond single-source runtime. We apply VARnet to the NEOWISE Single-Exposure Database, which holds nearly 200 billion apparitions over 10.5 years of infrared sources on the entire sky. This paper devises a pipeline in order to extract variable candidates from the NEOWISE data, serving as a proof of concept for both the efficacy of VARnet and methods for an upcoming variability survey over the entirety of the NEOWISE dataset. We implement models and simulations to synthesize unique light curves to train VARnet. In this case, the model achieves an F1 score of $0.91$ over a 4-class classification scheme on a validation set of real variable sources present in the infrared. With $\sim2000$ points per light curve on a GPU with 22GB of VRAM, VARnet produces a per-source processing time of $<53\mu s$. We confirm that our VARnet is sensitive and precise to both known and previously undiscovered variable sources. These methods prove promising for a complete future survey of variability with WISE, and effectively showcase the power of the VARnet model architecture.

Radio detections of stellar systems provide a window onto stellar magnetic activity and the space weather conditions of extrasolar planets, information that is difficult to attain at other wavelengths. There have been recent advances observing auroral emissions from radio-bright low-mass stars and exoplanets largely due to the maturation of low-frequency radio instruments and the plethora of wide-field radio surveys. To guide us in placing these recent results in context, we introduce the foremost local analogues for the field: Solar bursts and the aurorae found on Jupiter. We detail how radio bursts associated with stellar flares are foundational to the study of stellar coronae, and time-resolved radio dynamic spectra offers one of the best prospects of detecting and characterising coronal mass ejections from other stars. We highlight the prospects of directly detecting coherent radio emission from exoplanetary magnetospheres, and early tentative results. We bridge this discussion to the field of brown dwarf radio emission, in which their larger and stronger magnetospheres are amenable to detailed study with current instruments. Bright, coherent radio emission is also predicted from magnetic interactions between stars and close-in planets. We discuss the underlying physics of these interactions and implications of recent provisional detections for exoplanet characterisation. We conclude with an overview of outstanding questions in theory of stellar, star-planet interaction, and exoplanet radio emission, and the prospects of future facilities in answering them.

Fast radio bursts (FRBs) offer unique probes of diverse cosmological phenomena due to their characteristic properties, including short duration timescale and high dispersion measure. This study investigates two distinct theoretical frameworks: the Gertsenshtein-Zel'dovich (GZ) mechanism for ultra-high-frequency gravitational waves (GWs) and fraction of dark matter in primordial mass black holes. We explore the hypothesis that ultra-high-frequency GWs could be responsible for FRB generation. Consequently, the detection of continuous GWs signal from the vicinity of an FRB by current or future detectors would disfavour merger-based FRB formation models and lend significant credence to the GZ theory, which postulates the existence of high-frequency GWs. Moreover, we examine the effects of modified gravity on the gravitational lensing of FRBs and thereby put constraints on the fraction of primordial mass black holes made up of dark matter. Our analysis suggests that modified gravity introduces a screening effect on lensing, analogous to the scattering effect by plasma on light rays. We further discuss the expected detection rates of FRBs as well as lensed FRBs with upcoming radio telescopes, primarily HIRAX.

In ultra-fast astronomical observations featuring fast transients on sub-$\mu$s time scales, the conventional Signal-to-Noise Ratio (SNR) threshold, often fixed at $5\sigma$, becomes inadequate as observational window timescales shorten, leading to unsustainably high False Alarm Rates (FAR). We provide a basic statistical framework that captures the essential noise generation processes relevant to the analysis of time series data from photon-counting detectors. In particular, we establish a protocol of defining detection limits in astronomical photon-counting experiments, such that a FAR-based criterion is preferred over the traditional SNR-based threshold scheme. We develop statistical models that account for noise sources such as dark counts, sky background, and crosstalk, and establish a probabilistic detection criterion applicable to high-speed detectors. We compare the performance of several detector technologies, including photon-counting CMOS/CCDs, SPADs, SiPMs, and PMTs, in detecting faint astronomical signals. These findings offer insights into optimizing detector choice for future ultra-fast astronomical instruments and suggest pathways for improving detection fidelity under rapid observational conditions.

The solar atmosphere is a complex environment with diverse species and varying ionization states, especially in the chromosphere, where significant ionization variations occur. This region transitions from highly collisional to weakly collisional states, leading to complex plasma state transitions influenced by magnetic strengths and collisional properties. These processes introduce numerical stiffness in multi-fluid models, imposing severe timestep restrictions on standard time integration methods. New numerical methods are essential to address these computational challenges, effectively managing the diverse timescales in multi-fluid and multi-physics models. The widely used time operator splitting technique offers a straightforward approach but requires careful timestep management to avoid stability issues and errors. Despite some studies on splitting errors, their impact on solar and stellar astrophysics is often overlooked. We focus on a Multi-Fluid Multi-Species (MFMS) model, which presents significant challenges for time integration. We propose a second-order Partitioned Implicit-Explicit Runge-Kutta (PIROCK) method that combines efficient explicit and implicit integration techniques with variable time-stepping and error control. Compared to a standard third-order explicit method and a first-order Lie splitting approach, the PIROCK method shows robust advantages in accuracy, stability, and computational efficiency. Our results reveal PIROCK's capability to solve multi-fluid problems with unprecedented efficiency. Preliminary results on chemical fractionation represent a significant step toward understanding the First-Ionization-Potential (FIP) effect in the solar atmosphere.

The wide fields of view, high sensitivity, and broad energy coverage of current X-ray and gamma-ray satellites, coupled with the high cadence observational strategy of some of them (recently Swift and Fermi) have been ideal for carrying out unprecedented studies of the variability properties of different classes of Galactic and extra-Galactic high-energy sources. These classes of objects range from nearby flaring stars to the most distant Active Galactic Nuclei (AGNs). In this paper we focus on some of the most energetic events, i.e. those powered by accretion (black hole binaries, ultra-luminous X-rays and Active Galactic Nuclei) until those leading to the first detections of GravitationalWaves (GWs), i.e. the Gamma-ray Bursts (GRBs), passing through the controversial Inter-Mediate Mass Black-Holes (IMBHs). We show the importance of X-ray and gamma-ray emission for the determination of the properties (mass, spin, inclination, height of the corona) of the compact objects residing in systems powered by accretion and the role of the LIGO and VIRGO interferometers in the case of the (compact objects) binary mergers. Gravitational waves allow the determination of the properties of the non-light-emitting compact objects (black-holes; BHs) in binary mergers (BH-BH) that otherwise would be nearly impossible. We show that for the case of neutron star (NS) mergers (NS-BH or NS-NS) electro-magnetic (EM) emission is still possible, and very powerful, being the responsible for the X-ray and gamma-ray emission of many GRBs (kilonovae). Very recently these transients have been discovered to show X-ray and gamma-ray variability patterns that could lead to very important insights into the properties of the binary mergers that originated them, opening a new view for their study.

Motivated by anomalies in cosmic microwave background observations, we investigate the implications of $f(Q, T)$ gravity in Bianchi type-I spacetime, aiming to characterize the universe's spatially homogeneous and anisotropic properties. By using a linear combination of non-metricity $Q$ and the energy-momentum tensor trace $T$, we parametrize the deceleration parameter and derive the Hubble solution, which we then impose in the Friedmann equations of $f(Q, T)$ gravity. Bayesian analysis is employed to find the best-fit values of model parameters, with $1-\sigma$ and $2-\sigma$ contour plots illustrating the constraints from observational data, including $H(z)$ data and the Pantheon+ sample. Our analysis reveals a transition from a decelerated to an accelerated expansion phase, with the present deceleration parameter indicating an accelerating universe. The energy density gradually decreases over time, approaching zero for the present and future, indicating continuous expansion. The anisotropic pressure, initially notably negative, transitions to slightly negative values, suggesting the presence of dark energy. The evolving equation of state parameter $\omega$ exhibits behavior akin to phantom energy, influenced by spacetime anisotropy. Violations of the null energy condition and the strong energy condition imply phantom-like behavior and accelerated expansion.

We present the results of the optical study of the new eclipsing polar Gaia23cer. We analyzed the brightness variability of the polar in high ($\langle r \rangle \approx 16.5\mathrm{\,mag}$) and low ($\langle r \rangle \approx 19.2\mathrm{\,mag}$) states. The system has an orbital period $P_{orb} = 102.0665 \pm 0.0015$ min and exhibits deep eclipses with a duration $\Delta t_{ecl} = 401.30 \pm 0.81$ s. The spectra have a red cyclotron continuum with the Zeeman H$\alpha$ absorption triplet forming in a magnetic field with a strength of $15.2 \pm 1.1$ MG. The source of emission lines has a high radial velocity semiamplitude ($K\approx 450$km/s) and its eclipse lags behind the eclipse of the white dwarf. The mass $M_1=0.79 \pm 0.03 M_{\odot}$ and temperature $T=11350 \pm 650 K$ of the white dwarf have been found by modelling the spectral energy distribution. The eclipse duration corresponds to a donor mass $M_2 = 0.10-0.13M_{\odot}$ and an orbital inclination $i=84.3-87.0^{\circ}$. The donor temperature was estimated to be $T\approx 2900K$ by modelling the elliptical variability and eclipse depth.

The GALEX Nearby Young Star Search (GALNYSS) yielded the identification of more than 2000 late-type stars that, based on their ultraviolet and infrared colors and pre-Gaia proper motions, are potentially of age < 200 Myr and lie within ~120 pc of Earth. We present the results of a campaign of medium- and high-resolution optical spectroscopy of 471 GALNYSS stars aimed at confirming their youth and their potential membership in nearby young stellar moving groups. We present radial velocity (RV), Li absorption, and H-$\alpha$ emission measurements for these spectroscopically observed GALNYSS stars, and assess their Li absorption and optical emission-line properties and infrared excesses. Our RV measurements are combined with literature and Gaia DR3 RV measurements and Gaia DR3 astrometry and photometry to obtain the spatiokinematics and color-magnitude positions of GALNYSS stars. We use these results to assess membership in the TW Hya, Tuc-Hor, Carina, Columba, and Argus Associations and the $\beta$-PMG and AB Dor moving groups. We have identified 132 stars as candidate members of these seven groups; roughly half of these candidates are newly identified on the basis of data presented here. At least one-third of the 132 candidates are spectroscopic and/or photometric binaries and/or have comoving (visual) binary companions in Gaia DR3. The contingent of young, low-mass stars in the solar vicinity we identify here should provide excellent subjects for future direct imaging exoplanet surveys and studies of the early evolution of low-mass stars and their planetary progeny.

(abridged) Characterizing duty cycles of recurrent phases of dormancy and activity in supermassive black holes in active galactic nuclei is crucial in understanding impact of energy released on host galaxies and their evolution. However, identifying sources in quiescent and restarted phases is challenging. Our goal is to identify and characterize a substantial sample of radio galaxies in restarted phase and explore core prominence as a signature of this activity. We expand our prior study from a $30\,\mathrm{deg^2}$ area in Lockman Hole to a larger $424\,\mathrm{deg^2}$ region in Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) field using visually selected core-dominated radio galaxies. We used 144 MHz LOFAR survey images to identify galaxies with restarting jets. By applying selection criteria, including radio core dominance, low surface brightness extended emission, spectral index properties, and morphology, we found 69 candidate restarted radio galaxies. These candidates show diverse intrinsic morphology, spanning FRI, FRII, core-with-halo, and asymmetric forms, suggesting different progenitors. Among these, nine galaxies exhibit ultra-steep spectrum extended emission combined with high radio core prominence, representing previous and current epochs of jet activity. This subset supports a model where the switch-on and switch-off mechanism occurs with fast duty cycle. The restarted candidates span radio luminosities from log$_{10}$(L$_\mathrm{144 MHz}$/$\mathrm{WHz^{-1}}$) = 23.24 to 26.80, with sizes between 88 and 1659 kpc, including 16 giant radio galaxies. Their total stellar content aligns with massive elliptical galaxies. Our findings at $z<0.4$ suggest that many restarting galaxies are not found in rich cluster environments, consistent with broader radio-galaxy population. This study confirms core prominence as an effective parameter for selecting restarted radio sources.

We investigate the kinematics and dynamics of the molecular and ionized gas in the host galaxies of three Palomar-Green quasars at low redshifts, benefiting from the archival millimeter-wave interferometric and optical integral field unit data. We study the kinematics of both cold molecular and hot ionized gas by analyzing the CO and H$\alpha$ data cubes, and construct the mass distributions of our sample through gas dynamics, utilizing a priori knowledge regarding the galaxy light distribution. We find no systematic offset between the stellar mass derived from our dynamical method and that from the broad-band photometry and mass-to-light ratio, suggesting the consistency of both methods. We then study the kinetic pressure and the weight of the interstellar medium using our dynamical mass model. By studying the relationship between kinetic pressure and gravitational pressure of the quasar host galaxies, we find an equivalence in the hydrostatic equilibrium states of ISM in the quasar host galaxies, similar to the result of gas equilibrium in normal star-forming galaxies, suggesting minimal quasar feedback. Regarding non-circular motion as indicative of quasar-driven outflows, we observe an exceptionally low coupling efficiency between molecular gas outflow and AGN bolometric luminosities. These results demonstrate the marginal influence of the central engine on the properties of cold molecular gas in quasar host galaxies.

Since the discovery of the accelerated expansion of the Universe in 1998, modified gravity (MG) theories have attracted considerable attention as alternatives to dark energy (DE). While distinguishing the effects of MG from those of DE using cosmic background expansion alone is difficult, the large-scale structure is expected to differ significantly. Among the plethora of MG models, we are particularly interested in those that introduce a scale dependence in the growth of perturbations; specifically, theories that introduce fifth forces mediated by scalar fields with a finite range accessible to cosmological probes. This is the case with $f(R)$ gravity, which is widely regarded as the most studied model in cosmology. In this work, we utilize, for the first time, the full-shape power spectrum of galaxies to constrain scale-dependent modified gravity theories. By using BOSS DR12 dataset, along with a BBN prior on $ \omega_b $ and a Planck 2018 prior on $ n_s $, we obtain an upper bound of $ |f_{R0}| < 5.89 \times 10^{-6} $ at 68% confidence level (c.l.) and $ < 1.53 \times 10^{-5} $ at 95% c.l. for the Hu-Sawicki ($ n=1 $) model. We discuss that it is highly unlikely these constraints will be significantly improved by future galaxy spectroscopic catalogs, such as DESI and Euclid.

In spite of the importance of studying the cosmic generation of heavy elements through the r-process, the detection of kilonova resulted from a merger of neutron star binaries is still a challenge task. In this paper, we show that the Visible Telescope (VT) onboard the on-going SVOM space mission is powerful for identifying kilonova candidates associated with short gamma-ray bursts (SGRBs) up to a distance of 600Mpc. A significant color variation, turn blue and then turn red, is revealed by calculating the light curves in both red and blue channels of VT by a linear combination of an afterglow and an associated kilonova. The maximum color variation is as high as $\sim0.5-1$ mag, which is far larger than the small photometry error of $\sim0.2$ mag of VT for a point source with a brightness of 23 mag. Up to a distance of 600Mpc, $\sim1-2$ kilonova candidates per year are predicted to be identified by VT.

To address the discrepancy where disk sizes exceed those predicted by standard models, we explore two extensions to disk size estimates within the UV/optical wavelength range: disk winds and color correction. We provide detailed, self-consistent derivations and analytical formulas, including those based on a power-law temperature approximation, offering efficient tools for analyzing observational data. Applying our model to four type I AGNs with intensive reverberation mapping observations, we find a shallower temperature slope ($T\propto R^{-0.66}$, compared to $R^{-3/4}$ traditionally) and a color correction factor ($f_{\rm col} \approx 1.6$), consistent with previous studies. We observe a positive correlation between accretion rate and color correction with black hole mass. However, the small sample size limits our conclusions. The strong degeneracy between the temperature slope and accretion rate suggests that incorporating flux spectra or spectral energy distributions could improve fitting accuracy. Our simulation approach offers a rapid and effective method for reverberation mapping.

The systematic errors are inevitable in Gaia published astrometric data. Lindegren et al. (L21) proposed a global recipe to correct for the GEDR3 parallax zero point offset, which did not consider the Galactic plane. The applicability of their correction model to the Galactic plane remains uncertain. We attempt to have an independent investigation into the sample dependence of the L21 correction, and its applicability to the Galactic plane. We collect various samples, including quasars, binaries, and sources with parallaxes from other surveys or methods, to validate the L21 correction, especially in the Galactic plane. We conclude that the L21 correction exhibits sample dependence, and does not apply effectively to the Galactic plane. We present a new parallax bias correction applying to the Galactic plane, offering improvements over the existing L21 correction. The correction difference between L21 and this work can go up to 0.01 mas within certain ranges of magnitude and colour. This work provides an additional recipe for users of Gaia parallaxes, especially for sources located near the Galactic plane.

In the immediate sunspots' vicinity -- their superpenumbra -- 3-minute line-of-sight (LOS) velocity oscillations dominate in the photosphere and chromosphere. Oscillations of similar periods are also registered in the transition region and lower corona above active regions. The goal of the work is to clarify whether these LOS velocity oscillations are manifestations of Alfvénic waves in the lower solar atmosphere. The study is based on three sunspots with the use of the Solar Dynamics Observatory data. Additional observations of a sunspot were carried out at ground-based Automated Solar Telescope. We use narrow-band frequency filtration (5.6--5.8 mHz) of the LOS velocity, magnetic field, and intensity signals of the Fe i 6173 A spectral line. For the analysis, we used 90-minute long series. We come to the conclusion that the 3-minute oscillations in the LOS velocity signals result from magnetoacoustic waves rather than from Alfvénic waves. However, oscillations registered in magnetic field signals indicate that Alfvénic waves may be present already in the photosphere. Further research requires simultaneous observations of LOS velocity, magnetic field strength, spectral line width, and intensity carried out at two heights of the solar atmosphere.

The Polstar small explorer concept is for an ultraviolet (UV) spectropolarimetry space telescope mission with a focus on massive star astrophysics. The instrument waveband will be from 115 nm - 286 nm for spectroscopy and 122 nm - 286 nm for polarimetry. All 4 Stokes parameters, IQUV, will be measured at a resolving power of R=20,000 (15 km/s velocity resolution). The telescope aperture will be 40 cm with an effective area of about 22 cm^2 at a reference wavelength of 150 nm. The thrust of the science goals will be to determine the astrophysics of angular momentum exchange and transport, and consequences for massive star properties and evolution. This includes the effects of rapid to critical rotation for individual stars (magnetic and non-magnetic), and the effects of mass transfer for massive binaries, including identification of stripped core stars. If selected by the NASA/SMEX program, Polstar would launch around 2031 and observe ~300 stars to achieve science goals. The mission will include a Guest Observer program to advance discovery in other areas of astrophysics.

Context: Dust grains in the interstellar medium are heated by the integrated radiation from stars in the Milky Way. A knowledge of the Local Interstellar Radiation Field (LISRF) is thus necessary to interpret observations of dust emission in the infrared and constrain (some of) the properties of interstellar grains. The LISRF representation most widely used in dust modelling still dates back to the seminal works of Mezger et al. (1982) and Mathis et al. (1983). Aims: A new version of the LISRF is presented, starting from the photometry of the Gaia Data Release 3 (DR3) and revisiting available data, among which observations from the Pioneer 10 and 11 probes. Methods: The LISRF contribution by direct starlight is estimated in the Gaia bands by summing fluxes of all stars in DR3; the LISRF is extrapolated from the optical to the ultraviolet and near-infrared, using the astrophysical parameters provided by DR3 for a subsample of Gaia stars; the correlation between dust emission at 100 um and residual diffuse emission in the Pioneer and other available maps is exploited to derive the contribution of dust-scattered starlight to the LISRF. Results: The new LISRF is significantly redder and emits ~30% more energy than the old model. The old LISRF is almost a factor two lower in the near-infrared, while in the optical it accounts only for direct starlight. For |b|<50 degree, diffuse starlight contributes on average to ~25% of the total radiation, 3x more than what predicted by literature estimates at high Galactic latitude. Conclusions: The new LISRF can modify the predicted mid-infrared dust emission beyond the uncertainties normally assumed between dust models and observational constraints; these differences should be taken into account to redefine the properties of small grains and of the carriers of the mid-infrared emission bands.

A significative fraction of high mass stars sail away through the interstellar medium of the galaxies. Once they evolved and died via a core collapse supernova, a magnetized, rotating neutron star (a pulsar) is usually their leftover. The immediate surroundings of the pulsar is the pulsar wind, which forms a nebula whose morphology is shaped by the supernova ejecta, channeled into the circumstellar medium of the progenitor star in the pre supernova time. Consequently, irregular pulsar wind nebulae display a large variety of radio appearances, screened by their interacting supernova blast wave or harboring asymmetric up down emission. Here, we present a series of 2.5 dimensional non relativistic magnetohydrodynamical simulations exploring the evolution of the pulsar wind nebulae (PWNe) generated by a red supergiant and a Wolf Rayet massive supernova progenitors, moving with Mach number M eq. to 1 and M eq. to 2 into the warm phase of the galactic plane. In such a simplified approach, the progenitors direction of motion, the local ambient medium magnetic field, the progenitor and pulsar axis of rotation, are all aligned, which restrict our study to peculiar pulsar wind nebula of high equatorial energy flux. We found that the reverberation of the termination shock of the pulsar wind nebulae, when sufficiently embedded into its dead stellar surroundings and interacting with the supernova ejecta, is asymmetric and differs greatly as a function of the past circumstellar evolution of its progenitor, which reflects into their projected radio synchrotron emission. This mechanism is particularly at work in the context of remnants involving slowly moving or very massive stars. We find that the mixing of material in plerionic core collapse supernova remnants is strongly affected by the asymmetric reverberation in their pulsar wind nebulae.

Three black holes (BHs) in wide binaries - Gaia BH1, BH2 and BH3 - were recently discovered. The likely progenitors of the BHs were massive stars that experienced a supergiant phase, reaching radii of ~1000 Rsun, before collapsing to form the BH. Such radii are difficult to accommodate with the present-day orbits of BH1 and BH2 - with semi-major axes of 1.4 and 3.7 au, respectively. In this letter, we show that the maximal radii of the supergiants are not necessarily so large, and realistic stellar evolution models, with some assumed overshooting above the convective core into the radiative stellar envelope, produce substantially smaller maximal radii. The limited expansion of supergiants is consistent with the empirical Humphreys-Davidson limit - the absence of red supergiants above an upper luminosity limit, notably lower than the highest luminosity of main-sequence stars. We propose that the evolution that led to the formation of Gaia BH1 and BH2 simply did not involve an expansion to the cool supergiant phase.

Modern disc galaxies commonly have distinct thin and thick discs that are separable in some combination of their kinematics, radial structure, chemistry and/or age. The formation mechanisms of the two discs and the timing of their onset remain open questions. To address these questions, we select edge-on galaxies from flagship JWST programs and investigate their disc structures in rest-frame, near-infrared bands. For the first time, we identify thick and thin discs at cosmological distances, dating back over 10 Gyr, and investigate their decomposed structural properties. We classify galaxies into those requires two discs (thin and thick discs) and those well fitted by a single disc. Disc radial sizes and vertical heights correlate strongly with total galaxy mass and/or disc mass, independent of cosmic time. The structure of thick discs resemble discs found in single-disc galaxies, suggesting that galaxies form a thick disc first followed by thin disc formation. The transition from single to double discs occurred around 8 Gyr ago in high-mass galaxies ($10^{9.75} - 10^{11}M_\odot$), earlier than the transition which occurred 4 Gyr ago in low-mass galaxies ($10^{9.0} - 10^{9.75}M_\odot$), indicating sequential formation proceeds in a "downsizing" manner. Toomre $Q$-regulated disc formation explains the delayed thin disc formation in low-mass galaxies, leading to the observed anti-correlation between the thick-to-thin disc mass ratio and total galaxy mass. Despite the dominant sequential formation, observations suggest that thick discs may continue to build up mass alongside their thin-disc counterparts.

Context. The spectra of fast-rotating A-type stars have strongly broadened absorption lines. This effect causes blending of the absorption lines, hindering the measurement of the abundances of the elements that are in the stellar photosphere. Aims. As the exoplanet transits across its host star, it obscures the stellar spectrum that is emitted from directly behind the planet. We aim to extract this obscured spectrum because it is less affected by rotational broadening, resolving the blending of weak lines of elements that would otherwise remain inaccessible. This allows us to more precisely measure the metal abundances in ultra-hot Jupiter systems, many of which have fast rotating host stars. Methods. We develop a novel method that isolates the stellar spectra behind the planet during a spectral time-series, and reconstructs the disc-integrated non-broadened spectrum of the host star. We have systematically tested this method with model-generated spectra of the transit of WASP-189 b across its fast-rotating A-type host star, assessing the effects of limb darkening, choice of absorption lines, signal to noise regime; and demonstrating the sensitivity to photospheric parameters ($T_{\text{eff}}$, $\log g$) and elemental abundances. We apply the method to observations by the HARPS high-resolution spectrograph. Results. For WASP-189, we obtain the metallicity and photospheric abundances for several species previously not reported in literature, Mg, Ca and Ti, with significantly improved accuracy compared to the ordinary broadened stellar spectrum. This method can be generally applied to other transiting systems in which abundance determinations via spectral synthesis are imprecise due to severe line blending. It is important to accurately determine the photospheric properties of exoplanet host stars, as it can provide further insight into the formation and evolution of the planets.

In this paper, we explore a novel framework for explaining the mass and radius relationships of observed neutron stars by considering strange stars (SSs) admixed with mirror dark matter (MDM). We develop a theoretical model that incorporates non-commutative algebra to describe the interactions between ordinary strange quark matter (SQM) and MDM, which are predicted to form compact objects that could explain recent astrophysical data, including observations of PSR J0740+6620, PSR J0030+0451, PSR J0437-4715, and the central compact object in HESS J1731-347. Notably, we demonstrate that the exotic mass-radius measurement of XTE J1814-338 can be explained by the presence of a mirror SS with an ordinary SQM core. In contrast to other explanations based on boson stars, our SS+MDM model offers a natural explanation for this system. We provide detailed mass-radius comparisons with observational data and discuss future observations that could test the predictions of our model, offering new insights into neutron star structure and the role of dark matter in compact objects.

Relativistic magnetized jets, originating near black holes, are observed to exhibit sub-structured flows. In this study, we present synthetic synchrotron emission signatures for different lines of sight and frequencies, derived from three-dimensional relativistic magneto-hydrodynamic simulations of pc-scale AGN jets. These simulations apply different injection nozzles, injecting steady, variable, and precessing jets. Extending our previous study, here, we have developed a bridge to connect jet dynamics and particle acceleration within relativistic shocks with non-thermal radiation dominant in jets. The emission is derived from Lagrangian particles - injected into the jet and following the fluid - accelerated through diffusive shock acceleration and subsequently cooled by emitting energy via synchrotron and inverse-Compton processes. Overall, the different shocks structures lead to the formation of numerous localized emission patterns - interpreted as jet knots. These knot patterns can fade or flare, also as a consequence of merging or Doppler boosting, leading to jet variability. We find knots with high-enough pattern speed supposed to be visible as superluminal motion <~5c. Synchrotron spectra of all jets reveal double-humped structures, reflecting multiple electron populations characterized by the nature of underlying shock and their age. The precessing jet is the most powerful emitter, featuring a spectrum flatter than the steady and the variable jet. The emission, although essentially governed by the acceleration through shocks, depends on the cooling history of the particle as well. Overall, the continuous re-acceleration of electrons through shocks along the jet we found, is an essential prerequisite for observing extended jet emission over large time-scales and length-scales.

CORSIKA 8 is a new framework for simulations of particle cascades in air and dense media implemented in modern C++17, based on past experience with existing codes, in particular CORSIKA 7. The flexible and modular structure of the project allows the development of independent modules that can produce a fully customizable particle shower simulation. The radio module in particular is designed to treat the electric field calculation and its propagation through complex media to each observer location in an autonomous and flexible way. It already allows for the simultaneous simulation of the radio emission calculated with two independent time-domain formalisms, the "Endpoint formalism" as previously implemented in CoREAS and the "ZHS" algorithm as ported from ZHAireS. The design acts as the baseline interface for current and future development for the simulation of radio emission from particle showers in standard and complex scenarios, such as cross-media showers penetrating from air into ice. In this work, we present the design and implementation of the radio module in CORSIKA 8, along with validation studies and a direct comparison of the radio emission from air showers simulated with CORSIKA 8, CORSIKA 7 and ZHAireS. We also present the impact of simulation details such as the step size of simulated particle tracks on radio-emission simulations and perform a direct comparison of the "Endpoints" and "ZHS" formalisms for the same underlying air showers. Finally, we present an in-depth comparison of CORSIKA 8 and CORSIKA 7 for optimum simulation settings and discuss the relevance of observed differences in light of reconstruction efforts for the energy and mass of cosmic rays.

This study aims to identify potential exoplanet signals from nearby stars with resolved debris discs. However, the high activity of many stars with debris discs limits the detection of periodic signals. Our study is constrained to a sample of 29 stars that have appropriate radial velocity data and debris disc measurements sufficient to resolve their inclination. Our results confirm and update previous findings for exoplanets around HD 10647, HD 115617, HD 69830, GJ 581, HD 22049, and HD 142091, and we identify long-term activity signals around HD 207129 and HD 202628. We utilize the inclination angles of the debris discs, assuming co-planarity between debris disc and exoplanet orbit, to determine the "disc-aligned" masses of radial velocity exoplanets in this study. The "disc-aligned" masses of HD 69830 b, HD 69830 c, and 61 Vir b suggest that they may be classified as 'hot' or 'warm' Jupiters and so might be nearby examples of planets that have undergone recent type-II disc migration.

Observations show that the X-ray emission of the accreting weakly magnetized neutron stars is polarized. Here, we develop a theoretical model, where we assume the emission of the accreting neutron star coming from the spreading layer, the extension of the boundary between the disk and the neutron star surface onto the surface. We then calculate the Stokes parameters of the emission accounting for relativistic aberration and gravitational light bending in the Schwarzschild metric. We show that regardless of the geometry, for the spreading layer, we cannot expect the polarization degree to be higher than 1.5%. Our results have implications with regard to the understanding of the X-ray polarization from weakly magnetized neutron stars observed with the Imaging X-ray Polarimetry Explorer and the future enhanced X-ray Timing and Polarimetry mission.

This study investigates the active regions of the M3.0V star G 80-21 using the observed data from the CARMENES project with synthetic spectra generated by the RH1.5D radiative transfer code. The CARMENES project aims to search for exoplanets around M dwarfs using high-resolution near-infrared and optical echelle spectrographs. By comparing the observed data and models for the chromospheric lines of H$_\alpha$ and the bluest Ca II infrared triplet line, we obtain the best-fit models for this star. The optimal fitting for the observed spectrum of G 80-21 is achieved by employing two active areas in conjunction with an inactive regions, with a calcium abundance of [Ca/H] = $-$0.4. This combination successfully fits all the observed data across varying ratios. The minor active component consistently comprises approximately 18\% of the total (ranging from 14\% to 20\%), which suggests that the minor active component is likely located in the polar regions. Meanwhile, the major active component occupies a variable proportion, ranging from 51\% to 82\%. Our method allows for the determination of the structure and size of stellar chromospheric active regions by analyzing high-resolution observed spectra.

All types of interaction of a magnetized plasma flow with an obstacle (magnetized or not) are considered, and those susceptible to produce a radio signature are identified. The role of the sub-Alfvénic or super-Alfvénic character of the flow is discussed. Known examples in the solar system are given, as well as extrapolations to star-planet plasma interactions. The dissipated power and the fraction that goes into radio waves are evaluated in the frame of the radio-magnetic scaling law, the theoretical bases and validity of which are discussed in the light of recent works. Then it is shown how radio signatures can be interpreted in the frame of the cyclotron-maser theory (developed for explaining the generation of solar system planetary auroral and satellite-induced radio emissions) for deducing many physical parameters of the system studied, including the planetary or stellar magnetic field. Recent detections of such radio signatures with new generation low-frequency radiotelescopes and future prospects are then outlined.

We present the Planetary Enriched White Dwarf Database (PEWDD), a collection of published photospheric abundances of white dwarfs accreting planetary debris alongside additional information relevant to metal-enrichment and the presence of infrared excesses, emission lines, and binary companions. PEWDD contains at the time of publishing information on 1739 white dwarfs and will be kept up-to-date with information from new publications. A total of 24 photospheric metals are recorded and are linked to accretion of exo-planetary material. The overall properties of metal-enriched white dwarfs are severely affected by observational selection effects; in particular we find that what metals are detectable strongly correlates with the effective temperature. By considering metal-enriched white dwarfs which have abundances measured by different methods, we find a spread that is comparable with the often quoted ad-hoc estimated abundance uncertainties, i.e. ~0.1-0.2dex. We draw attention to a dichotomy in the median accretion rates for metal-enriched H- and He-dominated white dwarfs, with M_{acc, H} = 7.7e7g/s and M_{acc, He} = 8.7e8g/s when extrapolating bulk compositions from bulk Earth Ca abundance. We identify 40 metal-enriched white dwarfs in binary systems and find evidence that enrichment is suppressed by binary companions within 200au.

The importance of Very Long Baseline Interferometry (VLBI) for pulsar research is becoming increasingly prominent and receiving more and more attention. In this paper, we present pathfinding pulsar observation results with the Chinese VLBI Network (CVN) incorporating the Five-hundred-meter Aperture Spherical radio Telescope (FAST). On MJD 60045 (April 11th, 2023), PSRs B0919+06 and B1133+16 were observed with the phase-referencing mode in the L-band using four radio telescopes (FAST, TianMa, Haoping and Nanshan) and correlated with the pulsar binning mode of the distributed FX-style software correlator in Shanghai. After further data processing with the NRAO Astronomical Image Processing System (AIPS), we detected these two pulsars and fitted their current positions with accuracy at the milliarcsecond level. By comparison, our results show significantly better agreement with predicted values based on historical VLBI observations than that with previous timing observations, as pulsar astrometry with the VLBI provides a more direct and model-independent method for accurately obtaining related parameters.

Extreme continuum variability in AGNs can indicate extreme changes in accretion flows onto supermassive black holes. We explore the multiwavelength nature of a continuum flare in the Seyfert LCRS B040659.9$-$385922. The all-sky X-ray surveys conducted by the eROSITA showed that its X-ray flux increased by a factor of roughly five over six months, and concurrent optical photometric monitoring with the ATLAS showed a simultaneous increase. We triggered a multiwavelength follow-up monitoring program (XMM, NICER; optical spectroscopy) to study the evolution of the accretion disk, broad-line region, and X-ray corona. During the campaign, X-ray and optical continuum flux subsided over roughly six months. We detected a soft X-ray excess near the flare peak and after it subsided, both exhibiting a power-law (nonthermal) behavior. We modeled the broadband optical/UV/X-ray spectral energy distribution at both the flare peak and post-flare times with the AGNSED model, incorporating thermal disk emission into the optical/UV and warm thermal Comptonization in the soft X-rays. Additionally, we find that the broad Heii $\lambda$4686 emission line fades significantly as the optical/UV/X-ray continuum fades, which could indicate a substantial flare of disk emission above 54 eV. We also observed a redshifted broad component in the H${\beta}$ emission line that is present during the high flux state of the source and disappears in subsequent observations. We witnessed a likely sudden strong increase in local accretion rate, which manifested itself via an increase in accretion disk emission and thermal Comptonization emission in the soft X-rays, followed by a decrease in accretion and Comptonized luminosity. The physical processes leading to such substantial variations are still an open question, and future continuous monitoring along with multi-wavelength studies will shed some light on it.

We use the James Webb Space Telescope (JWST) and its Mid-Infrared Instrument (MIRI) (5-28 um), to study the embedded HH 211 flow. We map a 0.95'x0.22' region, covering the full extent of the blue-shifted lobe, the central protostellar region, and a small portion of the red-shifted lobe. The jet driving source is not detected even at the longest mid-IR wavelengths. The overall morphology of the flow consists of a highly collimated jet, mostly molecular (H2, HD) with an inner atomic ([FeI], [FeII], [SI], [NiII]) structure. The jet shocks the ambient medium, producing several large bow-shocks, rich in forbidden atomic and molecular lines, and is driving an H2 molecular outflow, mostly traced by low-J, v=0 transitions. Moreover, 0-0 S(1) uncollimated emission is also detected down to 2"-3" (~650-1000 au) from the source, tracing a cold (T=200-400 K), less dense and poorly collimated molecular wind. The atomic jet ([FeII] at 26 um) is detected down to ~130 au from source, whereas the lack of H2 emission close to the source is likely due to the large visual extinction. Dust continuum-emission is detected at the terminal bow-shocks, and in the blue- and red-shifted jet, being likely dust lifted from the disk. The jet shows an onion-like structure, with layers of different size, velocity, temperature, and chemical composition. Moreover, moving from the inner jet to the outer bow-shocks, different physical, kinematic and excitation conditions for both molecular and atomic gas are observed. The jet mass-flux rate, momentum, and momentum flux of the warm H2 component are up to one order of magnitude higher than those inferred from the atomic jet component. Our findings indicate that the warm H2 component is the primary mover of the outflow, namely it is the most significant dynamical component of the jet, in contrast to jets from more evolved YSOs, where the atomic component is dominant.

The ionization of diffuse gas located far above the energetic midplane OB stars poses a challenge to the commonly accepted notion that radiation from OB stars is the primary ionization source for gas in galaxies. We investigated the sources of ionizing radiation, specifically leaking midplane HII regions and/or in situ hot low-mass evolved stars (HOLMES), in extraplanar diffuse ionized gas (eDIG) in a sample of eight nearby (17-52 Mpc) edge-on disk galaxies observed with the Multi Unit Spectroscopic Explorer (MUSE). We constructed a model for the photoionization of eDIG clouds and the propagation of ionizing radiation through the eDIG using subsequent runs of Cloudy photoionization code. Our model includes radiation originating both from midplane OB stars and in situ evolved stars and its dilution and processing as it propagates in the eDIG. We fit the model to the data using the vertical line ratio profiles of our sample galaxies, and find that while the ionization by in situ evolved stars is insignificant for most of the galaxies in our sample, it may be able to explain the enhanced high-ionization lines in the eDIG of the green valley galaxy ESO 544-27. Our results show that while leaking radiation from midplane HII regions is the primary ionization source for eDIG, in situ evolved stars can play a significant part in ionizing extraplanar gas in galaxies with low star forming rates.

We investigate a series of galaxy properties computed using the merger trees and environmental histories from dark matter only cosmological simulations, using the predictive semi-recurrent neural network outlined in Chittenden and Tojeiro (2023), and using stochastic improvements presented in our companion paper: Behera, Tojeiro and Chittenden (2024). We apply these methods to the dark matter only runs of the IllustrisTNG simulations to understand the effects of baryon removal, and to the gigaparsec-volume pure dark matter simulation Uchuu, to understand the effects of the lower resolution or alternative metrics for halo properties. We find that the machine learning model recovers accurate summary statistics derived from the predicted star formation and stellar metallicity histories, and correspondent spectroscopy and photometry. However, the inaccuracies of the model's application to dark simulations are substantial for low mass and slowly growing haloes. For these objects, the halo mass accretion rate is exaggerated due to the lack of stellar feedback, yet the formation of the halo can be severely limited by the absence of low mass progenitors in a low resolution simulation. Furthermore, differences in the structure and environment of higher mass haloes results in an overabundance of red, quenched galaxies. These results signify progress towards a machine learning model which builds high fidelity mocks based on a physical interpretation of the galaxy-halo connection, yet they illustrate the need to account for differences in halo properties and the resolution of the simulation.

Some X-ray binaries containing an energetic pulsar in orbit around a normal star accelerate particles to high energies in the shock cone formed where the pulsar and stellar winds collide. The magnetic field geometry in the acceleration region in such binaries is unknown. We performed the first measurement of the polarization of the X-ray synchrotron emission from a gamma-ray emitting binary system. We observed PSR B1259-63 with the Imaging X-ray Polarimetry Explorer (IXPE) during an X-ray bright phase following the periastron passage in June 2024. X-ray polarization is detected with a polarization degree of $8.3\% \pm 1.5\%$ at a significance of $5.3 \sigma$. The X-ray polarization angle is aligned with the axis of the shock cone at the time of the observation. This indicates that the predominant component of the magnetic field in the acceleration region is oriented perpendicular to the shock cone axis.

A hadron-quark phase transition (PT) may trigger supernova explosions during stellar core collapse. However, both success and failure have occurred in previous attempts to explode dying stars via this mechanism. We systematically explore the outcomes of the PT-induced collapse of proto-compact stars (PCSs), with spherically symmetric general relativistic hydrodynamic simulations and a controlled series of hybrid equations of state. Our results reveal the dependence of successful and failed explosions on the PT and quark matter characteristics. A small portion ($\sim\!0.04\%\!-\!1\%$) of the released binding energy $\Delta E_B$ transforms into the diagnostic explosion energy $E_{\rm exp,diag}$, which saturates at $\sim\!6\times10^{51}\,\rm{erg}$ near the black hole formation. We draw the phase diagrams for the possible fates of supernova explosions driven by hadron-quark PTs, where the control parameters are the onset density, energy gap of the PT, and the quark matter speed of sound. Our findings can be used to guide further investigations on PT-driven core-collapse supernovae and help identify hadron-quark PT-induced PCS collapse from future observations.

Neutron star (NS) mergers are known to produce heavy elements through rapid neutron capture (r-process) nucleosynthesis. Actinides are expected to be created solely by the r-process in the most neutron rich environments. Confirming if NS mergers provide the requisite conditions for actinide creation is therefore central to determining their origin in the Universe. Actinide signatures in kilonova (KN) spectra may yield an answer, provided adequate models are available in order to interpret observational data. In this study, we investigate actinide signatures in neutron rich merger ejecta. We use three ejecta models with different compositions and radioactive power, generated by nucleosynthesis calculations using the same initial electron fraction ($Y_e = 0.15$) but with different nuclear physics inputs and thermodynamic expansion history. These are evolved from 10 - 100 days after merger using the SUMO non-local thermodynamic equilibrium (NLTE) radiative transfer code. We highlight how uncertainties in nuclear properties, as well as choices in thermodynamic trajectory, may yield entirely different outputs for equal values of $Y_e$. We consider an actinide-free model and two actinide-rich models, and find that the emergent spectra and lightcurve evolution are significantly different depending on the amount of actinides present, and the overall decay properties of the models. We also present potential key actinide spectral signatures, of which doubly ionized $_{89}Ac$ and $_{90}Th$ may be particularly interesting as potential spectral indicators of actinide presence in KN ejecta.

Massive stars are recognized for their high degree of multiplicity, yet the mass ratio regime below 0.1 remains insufficiently explored. It is therefore unknown whether extremely low-mass (possibly substellar) companions can form and survive in the direct UV-irradiated environment of massive stars. In this paper, we discuss VLT/SPHERE IFS (0".15 - 0".85) observations of six massive O- and early B-type stars in Sco OB1 and M17 that each have a low-mass candidate companion. Two targets have companions that are brown dwarf candidates. The other four have candidate companions in the low end of the stellar mass regime ($\leq$ 0.30 M$_\odot$). For three of these, we have obtained a second epoch observation. At least two sources exhibit similar proper motion to that of their central star. However, given the expected proper motion of background objects, this does not imply certain companionship. We show how future follow-up observations of the brown dwarf candidate companions in $J$, $H$ and $L$ bands should allow for an unambiguous confirmation of their nature.

Infrared-faint white dwarfs are cool white dwarfs exhibiting significant infrared flux deficits, most often attributed to collision-induced absorption (CIA) from H$_2$-He in mixed hydrogen-helium atmospheres. We present James Webb Space Telescope (JWST) near- and mid-infrared spectra of three such objects using NIRSpec (0.6-5.3 $\mu$m) and MIRI (5-14 $\mu$m): LHS 3250, WD J1922+0233, and LHS 1126. Surprisingly, for LHS 3250, we detect no H$_2$-He CIA absorption at 2.4 $\mu$m, instead observing an unexpected small flux bump at this wavelength. WD J1922+0233 exhibits the anticipated strong absorption feature centered at 2.4 $\mu$m, but with an unexpected narrow emission-like feature inside this absorption band. LHS 1126 shows no CIA features and follows a $\lambda^{-2}$ power law in the mid-infrared. LHS 1126's lack of CIA features suggests a very low hydrogen abundance, with its infrared flux depletion likely caused by He-He-He CIA. For LHS 3250 and WD J1922+0233, the absence of a 1.2 $\mu$m CIA feature in both stars argues against ultracool temperatures, supporting recent suggestions that infrared-faint white dwarfs are warmer and more massive than previously thought. This conclusion is further solidified by Keck near-infrared spectroscopy of seven additional objects. We explore possible explanations for the unexpected emission-like features in both stars, and temperature inversions above the photosphere emerge as a promising hypothesis. Such inversions may be common among the infrared-faint population, and since they significantly affect the infrared spectral energy distribution, this would impact their photometric fits. Further JWST observations are needed to confirm the prevalence of this phenomenon and guide the development of improved atmospheric models.

The spectroscopic redshift measurement of BL Lac, a class of blazar, is challenging because its spectrum has no or weak emission lines ($\leqslant5Å$). We estimate the redshift by the photometric dropout technique for a sample of 64 blazars (59 BL Lacs and five blazar candidates of uncertainty type). Two telescopes are utilized to observe the sample: the {\it Swift} space telescope observes sources in $uvw2,\ uvm2,\ uvw1,\ u,\ b,\ v$ filters, while the ground-based telescopes SARA-CT/RM observed sources in $g',\ r,' \ i',\ z'$ filters. The photometric data are obtained using the {\it photozpy} package. We fit the photometric data by the LePhare package and report four new high-$z$ ($z>1.3$) BL Lacs at $2.03^{+0.07}_{-0.05}$, $1.84^{+0.10}_{-0.03}$, $2.04^{+0.16}_{-0.14}$, $2.93^{+0.01}_{-0.04}$ as well as upper limits for 50 sources. The work increased the number of high-$z$ BL Lacs found by this method up to 23. The high-$z$ sources are discussed in the context of the cosmic gamma-ray horizon, blazar sequence, Fermi blazar divide, and masquerading BL Lacs.

Vertical gas and dust flows in protoplanetary discs waft material above the midplane region in the presence of a protoplanet. This motion may alter the delivery of dust to the planet and its circumplanetary disc, as well as through a planetary-induced gap region and hence the inner disc chemistry. Here, we investigate the impact of a massive embedded planet on this material transport through the gap region. We use 3D global hydrodynamic simulations run using FARGO3D with gas and dust species to investigate the dust filtration and the origin of material that can make it through the gap. We find small dust particles can pass through the gap as expected from results in 2D, and that this can be considered in two parts - filtering due to the planetary-induced pressure maximum, and filtering due to accretion onto the planet. When gas accretion onto the planet is included, we find that the larger dust grains that cross the gap (i.e. those with $\mathrm{St} \sim 10^{-4}$) originate from regions near the mid-plane. We also find that dust and gas that enter the planet-carved gap region pass through the Hill sphere of the planet, where the temperature is likely to be strongly enhanced compared with the mid-plane regions from which this material originated. Considering the application of our simulations to a Jupiter-mass planet at $\sim 100\ \mathrm{AU}$, this suggests that CO ice is very likely to desorb from grains in the close proximity of the planet, without requiring any fine-tuning of the planet's location with respect to the CO snowline.

This work aims at providing some general tools for the analysis of water spectra as observed in protoplanetary disks with JWST-MIRI. We use 25 high-quality spectra from the JDISC Survey reduced with asteroid calibrators as presented in Pontoppidan et al. (2024). First, we present a spectral atlas to illustrate the clustering of water transitions from different upper level energies ($E_u$) and identify single (un-blended) lines that provide the most reliable measurements. With the atlas, we demonstrate two important excitation effects: one related to the opacity saturation of ortho-para line pairs that overlap, and the other to the sub-thermal excitation of $v=1-1$ lines scattered across the $v=0-0$ rotational band. Second, from this larger line selection we define a list of fundamental lines spanning $E_u$ from 1500 to 6000 K to develop simple line-ratio diagrams as diagnostics of temperature components and column density. Third, we report the detection of disk-rotation Doppler broadening of molecular lines, which demonstrates the radial distribution of water emission at different $E_u$ and confirms from gas kinematics a radially-extended $\approx$ 170-190 K reservoir recently suggested from the analysis of line fluxes. We also report the detection of narrow blue-shifted absorption from an inner disk wind in ro-vibrational H$_2$O and CO lines, which may be observed in disks at inclinations $> 50$ deg. We summarize these findings and tools into a general recipe that should be beneficial to community efforts to study water in planet-forming regions.

Liger is an adaptive optics (AO) fed imager and integral field spectrograph (IFS) designed to take advantage of the Keck All-sky Precision Adaptive-optics (KAPA) upgrade to the Keck I telescope. Liger adapts the design of the InfraRed Imaging Spectrograph (IRIS) for the Thirty Meter Telescope (TMT) to Keck by implementing a new imager and re-imaging optics. The performance of the imager is critical as it sequentially feeds the spectrograph and contains essential components such as the pupil wheel, filter wheel, and pupil viewing camera. We present the design and structural analysis of the Liger imager optical assembly including static, modal, and thermal simulations. We present the fabrication as well as the full assembly and characterization plan. The imager will be assembled bench-top in a clean room utilizing a coordinate-measuring machine (CMM) for warm alignment. To ensure optimal performance, the imager will be characterized in a test cryostat before integration with the full Liger instrument. This comprehensive approach to characterization ensures the precision and reliability of the imager, enhancing the observational capabilities of Liger and W.M. Keck Observatory.

Hot Jupiters are giant planets subject to intense stellar radiation. The physical and chemical properties of their atmosphere makes them the most amenable targets for the atmospheric characterization. In this paper we analyze the photometry collected during the secondary eclipses of the hot Jupiter WASP-3 b by CHEOPS, TESS and Spitzer. Our aim is to characterize the atmosphere of the planet by measuring the secondary eclipse depth in several passbands and constrain the planetary dayside spectrum. Our update of the stellar and planetary properties is consistent with previous works. The analysis of the occultations returns an eclipse depth of 92+-21 ppm in the CHEOPS passband, 83+-27 ppm for TESS and >2000 ppm in the IRAC 1-2-4 Spitzer passbands. Using the eclipse depths in the Spitzer bands we propose a set of likely emission spectra which constrain the emission contribution in the \cheops and TESS passbands to approximately a few dozens of parts per million. This allowed us to measure a geometric albedo of 0.21+-0.07 in the CHEOPS passband, while the TESS data lead to a 95\% upper limit of $\sim$0.2. WASP-3 b belongs to the group of ultra-hot Jupiters which are characterized by low Bond albedo (<0.3+-0.1), as predicted by different atmospheric models. On the other hand, it unexpectedly seems to efficiently recirculate the absorbed stellar energy, unlike similar highly irradiated planets. To explain this inconsistency, we propose that other energy recirculation mechanisms may be at play other than advection (for example, dissociation and recombination of H_2). Another possibility is that the observations in different bandpasses probe different atmospheric layers, making the atmospheric analysis difficult without an appropriate modeling of the thermal emission spectrum of WASP-3 b, which is not feasible with the limited spectroscopic data available to date.

Observations of the ultra-short period rocky exoplanet 55 Cancri e (55 Cnc e) indicate that the planet's dayside infrared radiation fluctuates by a factor of at least six on sub-weekly timescales, for unknown reasons. We propose a feedback mechanism where increased reflective clouds cool surface magma, subsequently reducing cloud formation, which may offer a potential explanation for these phenomena. In this mechanism, under less cloudy conditions, stellar radiation heats the surface magma, causing it to release more silicate vapor, which then condenses to form reflective clouds. Once formed, these clouds reduce stellar insolation at the surface, leading to surface cooling, which in turn reduces vapor supply, decreasing cloudiness. A time lag between the temperature increase of surface magma and the subsequent increase in cloudiness (likely due to lagged atmospheric transport of cloud-forming vapor) enables self-sustained oscillations in surface temperature and cloud reflectivity. These oscillations manifest as variations in both the emitted thermal radiation and the reflected stellar radiation, causing variability in secondary eclipse depths across different wavelengths without significantly affecting the transit depth. Using a simple model, we find that diverse planetary parameters can reproduce the observations. Additionally, we demonstrate that secondary eclipse depths at different wavelengths can oscillate out of phase, consistent with recent observations by the James Webb Space Telescope. Finally, we discuss observational strategies to test this proposed mechanism on 55 Cancri e. If confirmed, observable ocean-atmosphere dynamics on exoplanets would open a new window into the composition, evolution, and fate of rocky planet volatiles.

The stochastic formalism of inflation allows us to describe the scalar-field dynamics in a non-perturbative way. The correspondence between the diffusion and Schrödinger equations makes it possible to exhaustively construct analytical solutions in stochastic inflation. Those exact statistical quantities such as distribution and correlation functions have one-to-one correspondence to the exactly solvable solutions in non-relativistic quantum mechanics in terms of classical orthogonal polynomials. A class of such solutions is presented by means of isospectral Hamiltonians with an underlying symmetry called shape invariance.