Papers for Wednesday, May 15 2024

We explore the formation of the intragroup light (IGL) and intracluster light (ICL), representing diffuse lights within groups and clusters, since $z=1.5$. For this, we perform multi-resolution cosmological N-body simulations using the ``galaxy replacement technique" (GRT) and identify the progenitors in which the diffuse light stars existed when they fell into the groups or clusters. Our findings reveal that typical progenitors contributing to diffuse lights enter the host halo with the massive galaxies containing a stellar mass of $10 < \log M_{\rm{gal}}~[M_{\odot}]< 11$, regardless of the mass or dynamical state of the host halos at $z=0$. In cases where the host halos are dynamically unrelaxed or more massive, diffuse lights from massive progenitors with $\log M_{\rm{gal}}~[M_{\odot}]> 11$ are more prominent, with over half of them already pre-processed before entering the host halo. Additionally, we find that the main formation mechanism of diffuse lights is the stripping process of satellites, and a substantial fraction ($40-45\%$) of diffuse light stars is linked to the merger tree of the BCG. Remarkably, all trends persist for groups and clusters at higher redshifts. The fraction of diffuse lights in the host halos with a similar mass decreases as the redshift increases, but they are already substantial at $z=1.5$ ($\sim10\%$). However, it's crucial to acknowledge that detection limits related to the observable radius and faint-end surface brightness may obscure numerous diffuse light stars and even alter the main formation channel of diffuse lights.

A phase shift in the acoustic oscillations of cosmic microwave background (CMB) spectra is a characteristic signature for the presence of non-photon radiation propagating differently from photons, even when the radiation couples to the Standard Model particles solely gravitationally. It is well-established that compared to the presence of free-streaming radiation, CMB spectra shift to higher $\ell$-modes in the presence of self-interacting non-photon radiation such as neutrinos and dark radiation. In this study, we further demonstrate that the scattering of non-photon radiation with dark matter can further amplify this phase shift. We show that when the energy density of the interacting radiation surpasses that of interacting dark matter around matter-radiation equality, the phase shift enhancement is proportional to the interacting dark matter abundance and remains insensitive to the radiation energy density. Given the presence of dark matter-radiation interaction, this additional phase shift emerges as a generic signature of models featuring an interacting dark sector or neutrino-dark matter scattering. Using neutrino-dark matter scattering as an example, we numerically calculate the amplified phase shift and offer an analytical interpretation of the result by modeling photon and neutrino perturbations with coupled harmonic oscillators. This framework also explains the phase shift contrast between self-interacting and free-streaming neutrinos. Fitting models with neutrino-dark matter or dark radiation-dark matter interactions to CMB and large-scale structure data, we validate the presence of the enhanced phase shift, affirmed by the linear dependence observed between the preferred regions of the sound horizon angle $\theta_s$ and interacting dark matter abundance.

We present estimates of the ultraviolet (UV) and Lyman continuum flux density contributed by galaxies of luminosities from $M_{\rm UV}\approx -25$ to $M_{\rm UV}=-4$ at redshifts $5\leq z\leq 10$ using a galaxy formation model that reproduces properties of local dwarf galaxies down to the luminosities of the ultra-faint satellites. We characterize the UV luminosity function (LF) of galaxies and their abundance as a function of the ionizing photon emission rate predicted by our model and present accurate fitting functions describing them. Although the slope of the LF becomes gradually shallower with decreasing luminosity due to feedback-driven outflows, the UV LF predicted by the model remains quite steep at the luminosities $M_{\rm UV}\lesssim -14$. After reionization, the UV LF flattens at $M_{\rm UV}\gtrsim -12$ due to UV heating of intergalactic gas. However, before reionization, the slope of the LF remains steep and approximately constant from $M_{\rm UV}\approx -14$ to $M_{\rm UV}=-4$. We show that for a constant ionizing photon escape fraction the contribution of faint galaxies with $M_{\rm UV}>-14$ to the UV flux and ionizing photon budget is $\approx 40-60\%$ at $z>7$ and decreases to $\approx 20\%$ at $z=6$. Before reionization, even ultra-faint galaxies of $M_{\rm UV}>-10$ contribute $\approx 10-25\%$ of ionizing photons. If the escape fraction increases strongly for fainter galaxies, the contribution of $M_{\rm UV}>-14$ galaxies before reionization increases to $\approx 60-75\%$. Our results imply that dwarf galaxies fainter than $M_{\rm UV}=-14$, beyond the James Webb Space Telescope limit, contribute significantly to the UV flux density and ionizing photon budget before reionization alleviating requirements on the escape fraction of Lyman continuum photons.

We present a multi-frequency (144 MHz - 9 GHz) and multi-scale (5 pc - 50 kpc) investigation of the central radio galaxy in RBS 797, by means of JVLA, LOFAR (with international stations), e-Merlin, VLBA and EVN data. We investigate the morphological and spectral properties of the radio lobes, the jets, and the active core. We confirm the co-spatiality of the radio lobes with the four perpendicular X-ray cavities (see arXiv:2111.03679). The radiative ages of the E-W lobes ($31.4\pm6.6$ Myr) and of the N-S lobes ($32.1\pm9.9$ Myr) support a coeval origin of the perpendicular outbursts, that also have similar active phase duration ($\sim$12 Myr). For the inner N-S jets (on scales of $\leq10$ kpc), we (a) confirm the S-shaped jet morphology; (b) show the presence of two hotspots per jet with a similar spectral index; (c) estimate the age of the twisting jets to be less than $\sim8$ Myr. Based on these results, we determine that jet precession, with period $\sim$9 Myr, half-opening angle $\sim$24$^{\circ}$ and jet speed $\sim$0.01$c$, can explain the properties of the N-S jets. We also find that the synchrotron injection index has steepened from the large, older outbursts ($\Gamma\sim0.5$) to the younger S-shaped jets ($\Gamma\sim0.9$), possibly due to a transition from an FR I-like to an FR II-like activity. The VLBI data reveal a single, compact core at the heart of RBS 797, surrounded by extended radio emission whose orientation depends on the spatial scale sampled by the data. We explore several engine-based scenarios to explain these results. Piecing together the available evidence, we argue that RBS 797 likely hosts (or hosted) binary active SMBHs. This is still consistent with the detection of a single component in the VLBI data, since the predicted separation of the binary SMBHs ($\leq$0.6 pc) is an order of magnitude smaller than the resolution of the available radio data (5 pc).

Discovery of pulsations from a number of ultra-luminous X-ray (ULX) sources proved that accretion onto neutron stars can produce luminosities exceeding the Eddington limit by a couple of orders of magnitude. The conditions necessary to achieve such high luminosities as well as the exact geometry of the accretion flow in the neutron star vicinity are, however, a matter of debate. The pulse phase-resolved polarization measurements that became possible with the launch of the IXPE can be used to determine the pulsar geometry and its orientation relative to the orbital plane. They provide an avenue to test different theoretical models of ULX pulsars. In this paper we present the results of three IXPE observations of the first Galactic ULX pulsar Swift J0243.6+6124 during its 2023 outburst. We find strong variations of the polarization characteristics with the pulsar phase. The average polarization degree increases from about 5% to 15% as the flux dropped by a factor of three in the course of the outburst. The polarization angle (PA) as function of the pulsar phase shows two peaks in the first two observations, but changes to a characteristic sawtooth pattern in the remaining data set. This is not consistent with a simple rotating vector model. Assuming the existence of an additional constant polarized component, we were able to fit the three observations with a common rotating vector model and obtain constraints on the pulsar geometry. In particular, we find the pulsar angular momentum inclination with respect to the line-of-sight of 15-40 deg, the magnetic obliquity of 60-80 deg, and the pulsar spin position angle of -50 deg, which differs from the constant component PA of about 10 deg. Combining these X-ray measurements with the optical PA, we find evidence for a 30 deg misalignment between the pulsar spin and the binary orbital axis.

Primordial black holes (PBHs) are hypothetical objects that could have originated from density fluctuations in a very early phase of our Universe. Recent observations restrict the masses that such PBHs could have, if they are to constitute all of dark matter today: $10^{17} \, {\rm g} \leq m \leq 10^{23} \, {\rm g}$. With such low masses, general relativity predicts that the corresponding radii for the PBHs would be atomic or subatomic in size. When captured by a star, such a tiny PBH could exhibit an orbit completely or partially inside the body of the star, without significantly changing its mass for quite a long time. Here we examine the possible trajectories of a PBH that is captured by a Sun-like star. When in motion in the interior of the star, the amount of stellar mass that effectively interacts with the PBH will be a function of its distance to the center of the star. As a consequence, a strong effect on the shape of the orbits emerges, leading to PBH trajectories that could be open or closed, and exhibiting a rich variety of patterns.

We conduct a comprehensive search for transiting exomoons and exosatellites within 44 archival Spitzer light curves of 32 substellar worlds with estimated masses ranging between 3-30M$_{\rm Jup}$. This sample's median host mass is 16M$_{\rm Jup}$, inclusive of 14 planetary-mass objects, among which one is a wide-orbit exoplanet. We search the light curves for exosatellite signatures and implement a transit injection-recovery test, illustrating our survey's capability to detect $>$0.7R$_{\oplus}$ exosatellites. Our findings reveal no substantial ($>$5$\sigma$) evidence for individual transit events. However, an unusual fraction of light curves favor the transit model at the 2-3$\sigma$ significance level, with fitted transit depths consistent with terrestrial-sized (0.7-1.6R$_{\oplus}$) bodies. Comparatively, fewer than 2.2% of randomly generated normal distributions from an equivalent sample size exhibit a similar prevalence of outliers. Should one or two of these outliers represent a real exosatellite transit, it would imply an occurrence rate of $\eta = 0.61^{+0.49}_{-0.34}$ short-period terrestrial exosatellites per system, consistent with the known occurrences rates for both solar system moons and mid M-dwarf exoplanets. We explore alternative astrophysical interpretations for these outliers, underscoring that transits are not the only plausible explanation. For orbital periods $<$0.8 days, the typical duration of the light curves, we constrain the occurrence rate of sub-Neptunes to $\eta<$0.35 (95% confidence) and, if none of the detected outlier signals are real, the occurrence rate of terrestrial ($\sim$Earth-sized) exosatellites to $\eta<$0.51 (95% confidence). Forthcoming JWST observations of substellar light curves will enable detection of sub-Io-sized exosatellites, allowing for much stronger constraints on this exosatellite population.

Magnetic fields likely play an important role in the formation of young protostars. Multiscale and multiwavelength dust polarization observations can reveal the inferred magnetic field from scales of the cloud to core to protostar. We present continuum polarization observations of the young protostellar triple system IRAS 16293-2422 at 89 $\mu$m using HAWC+ on SOFIA. The inferred magnetic field is very uniform with an average field angle of 89$^\circ\pm$23$^\circ$ (E of N), which is different from the $\sim$170$^\circ$ field morphology seen at 850 $\mu$m at larger scales (> 2000 au) with JCMT POL-2 and at 1.3 mm on smaller scales (< 300 au) with ALMA. The HAWC+ magnetic field direction is aligned with the known E-W outflow. This alignment difference suggests that the shorter wavelength HAWC+ data is tracing the magnetic field associated with warmer dust likely from the outflow cavity, whereas the longer wavelength data are tracing the bulk magnetic field from cooler dust. Also, we show in this source the dust emission peak is strongly affected by the observing wavelength. The dust continuum peaks closer to source B (northern source) at shorter wavelengths and progressively moves toward the southern A source with increasing wavelength (from 22 $\mu$m to 850 $\mu$m).

Nebular spectral lines provide insight into the properties of the interstellar medium and the ionizing radiation within galaxies. The presence of high-energy ionization lines, such as He II, indicates the existence of hard ionizing photons, because the second ionization energy of helium is high (54 eV). The enigma surrounding the origin of these lines observed in star-forming galaxies persists, as stellar ionization cannot account for such high energy emission. This paper proposes that supersoft X-ray sources (SSS) may produce these high energy ionization lines in star-forming galaxies. We model the spectra of such sources using blackbody radiation and then add the blackbody to stellar population spectra to represent the overall spectra of galaxies. We then use a photoionization model to predict the resulting emission lines and compare the contribution of SSS to the observation of highly ionized lines in star-forming galaxies, both at low and high redshifts. We find that incorporating a blackbody with temperatures between kT = 20-100eV can boost the He II emission line ratio to the observed level. The blackbody temperature range aligns with the observed temperatures of the SSSs. The number of SSSs in spiral galaxies listed in Chandra catalogues, and our estimates of the total population, confirms that SSSs are promising candidates for the source of the He II ionization.

Effective temperatures, surface gravities, and iron abundances were derived for 109 stars observed by the K2 mission using equivalent width measurements of Fe I and Fe II lines. Calculations were carried out in LTE using Kurucz model atmospheres. Stellar masses and radii were derived by combining the stellar parameters with Gaia DR3 parallaxes, V-magnitudes, and isochrones. The derived stellar and planetary radii have median internal precision of 1.8%, and 2.3%, respectively. The radius gap near $\rm R_{planet}\sim 1.9 R_\oplus$ was detected in this K2 sample. Chromospheric activity was measured from the Ca II H and K lines using the Values of $\log R^\prime_{\rm HK}$ were investigated as a function of stellar rotational period (P$_{rot}$) and we found that chromospheric activity decreases with increasing P$_{rot}$, although there is a large scatter in $\log R^\prime_{\rm HK}$ ($\sim$0.5) for a given P$_{rot}$. Activity levels in this sample reveal a paucity of F & G dwarfs with intermediate activity levels (Vaughan-Preston gap). The effect that stellar activity might have on the derivation of stellar parameters was investigated by including magnetically-sensitive Fe I lines in the analysis and we find no significant differences between parameters with and without magnetically-sensitive lines, although the more active stars ($\log R^\prime _{\rm HK}>-5.0$) exhibit a larger scatter in the differences in $T_{\rm eff}$ and [Fe/H].

Despite being studied for many years, the structure of collisionless shocks is still not fully determined. Such shocks are known to be accelerators of cosmic rays, which, in turn, modify the shock structure. The shock width $\lambda$ is known to be connected to the cosmic rays (CRs) spectral index, $a$. Here, we use an instability analysis to derive the shock width in the presence of CRs. We obtain an analytical expression connecting the shock width to the CRs index and to the fraction of upstream particles that are accelerated. We find that when this fraction becomes larger than $\sim$~30\%, a new instability becomes dominant. The shock undergoes a transition where its width increases by a factor $\sim 8- 10$, and the CRs acceleration effectively ends. Our analysis is valid for strong, non-relativistic and unmagnetized shocks. We discuss the implication of these results to the expected range of CRs spectra and flux observed, and on the structure of non-relativistic collisionless shocks.

Astrophysical observations, theoretical models, and terrestrial experiments probe different regions of neutron star (NS) interior. Therefore, it is essential to consistently combine the information from these sources. This analysis requires multiple evaluations of Tolman Oppenheimer Volkoff equations which can become computationally expensive with a large number of observations. Further, multi-messenger astronomy requires rapid NS characterization via gravitational waves for efficient electromagnetic follow-up. In this work, we develop a novel neural network-based map from the EoS curve to the mass and radius of cold non-rotating NS. We estimate a speed-up of an order of magnitude when compared with the state-of-the-art RePrimAnd solver and an average error of 1e-3 when calculating the mass and radius of the neutron star. Additionally, we also develop neural network solvers for obtaining EoS curves from a physics conforming EoS model, FRZ$\chi_{1.5}$. We utilize this efficient continuous map to measure the sensitivity of model parameters of FRZ$\chi_{1.5}$ towards mass and radius. We show that 8 out of 18 parameters of this model are sensitive by at least three orders of magnitude higher than the remaining 10 parameters. This information will be useful in further speeding up, as well as probing the crucial parameter space, in the parameter estimation from astrophysical observations using this physics-conforming EoS model.

Lithium (Li) abundance is an age indicator for G, K, and M stellar types, as its abundance decreases over time for these spectral types. However, despite the observational efforts made over the past few decades, the role of rotation, activity, and metallicity in the depletion of Li is still unclear. We have investigated how Li depletion is affected by rotation and metallicity in G and K members of the roughly Pleiades-aged open cluster M35. To do so, we have collected a sample of 165 candidate members observed with the WIYN/Hydra spectrograph. In addition, we have taken advantage of three previous spectroscopic studies of Li in M35. As a result, we have collected a final sample of 396 stars which we have classified as members and non-members of the cluster. We have measured iron abundances, Li equivalent widths, and Li abundances for the 110 M35 members added to the existing sample by this study. Finally, rotation periods for cluster members have been obtained from the literature or derived from Zwicky Transient Facility light curves. As a result, we have confirmed that fast G and K rotators are Li-rich in comparison with slow rotators of similar effective temperature. Furthermore, while we derived subsolar metallicity for M35 from our spectra, the distribution of Li in this cluster is similar to those observed for the Pleiades and M34, which have solar metallicity and slightly different ages. In addition, we have shown that an empirical relationship proposed to remove the contribution of the Fe I line at 670.75 nm to the blended feature at 670.78 nm overestimates the contribution of this iron line for M35 members. We conclude that a 0.2-0.3 dex difference in metallicity makes little difference in the Li distributions of open clusters with ages between 100 and 250 Myr.

The expanding solar wind plasma ubiquitously exhibits anisotropic non-thermal particle velocity distributions. Typically, proton Velocity Distribution Functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. This phenomenon gives rise to VDFs that resemble a "hammerhead". In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam populations that drift with respect to each other. This configuration triggers a proton-beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing significant perpendicular heating at the expense of the parallel temperature. Furthermore, the proton perpendicular heating induces a hammerhead-like shape in the resulting VDF. Our results suggest that this mechanism probably contribute to producing the observed hammerhead distributions.

Context. Jupiter-family comets (JFCs), which originate from the Kuiper belt and scattered disk, exhibit low-inclination and chaotic trajectories due to close encounters with Jupiter. Despite their typically short incursions into the inner solar system, a notable number of them are on Earth-crossing orbits, with fireball networks detecting many objects on ``JFC-like'' (2 < TJ < 3) orbits. Aims. This investigation aims to examine the orbital dynamics of JFCs and comet-like fireballs over 10,000 yr timescales, focusing on the trajectories and stability of these objects in the context of gravitational interactions within the solar system. Methods. We employed an extensive fireball dataset from the Desert Fireball Network (DFN), European Fireball Network (EFN), Fireball Recovery and InterPlanetary Observation Network (FRIPON), and Meteorite Observation and Recovery Project (MORP), alongside telescopically observed cometary ephemeris from the NASA HORIZONS database. The study integrates 646 fireball orbits with 661 JFC orbits for a comparative analysis of their orbital stability and evolution. Results. The analysis confirms frequent Jupiter encounters among most JFCs, inducing chaotic orbital behavior with limited predictability and short Lyapunov lifetimes (about 120 years), underscoring Jupiter`s significant dynamical influence. In contrast, ``JFC-like'' meteoroids detected by fireball networks largely exhibit dynamics divergent from genuine JFCs, with 79-92% on ``JFC-like'' orbits shown not to be prone to frequent Jupiter encounters; in particular, only 1-5% of all fireballs detected by the four networks exhibit dynamics similar to that of actual JFCs. In addition, 22% (16 of 72) of near-Earth JFCs are on highly stable orbits, suggesting a potential main belt origin for some of the bodies.

The observation of a gamma-ray burst (GRB) associated with a supernova (SN) coincides remarkably with the energy output from a binary system comprising a very massive carbon-oxygen (CO) core and an associated binary neutron star (NS) by the Binary-Driven Hypernova (BdHN) model. The dragging effect in the late evolution of such systems leads to co-rotation, with binary periods on the order of minutes, resulting in a very fast rotating core and a binary NS companion at a distance of $\sim 10^5$ km. Such a fast-rotating CO core, stripped of its hydrogen and helium, undergoes gravitational collapse and, within a fraction of seconds, leads to a supernova (SN) and a newly born, fast-spinning neutron star ($\nu$NS), we name the emergence of the SN and the $\nu$NS as the SN-rise and $\nu$NS-rise. Typically, the SN energies range from $10^{51}$ to $10^{53}$ erg. We address this issue by examining 10 cases of Type-I BdHNe, the most energetic ones, in which SN accretion onto the companion NS leads to the formation of a black hole (BH). In all ten cases, the energetics of the SN events are estimated, ranging between $0.18$ and $12 \times 10^{52}$ erg. Additionally, in all 8 sources at redshift $z$ closer than $4.61$, a clear thermal blackbody component has been identified, with temperatures between $6.2$ and $39.99$ keV, as a possible signature of pair-driven SN. The triggering of the X-ray afterglow induced by the $\nu$NS-rise are identified in three cases at high redshift where early X-ray observations are achievable, benefits from the interplay of cosmological effects.

The process of cometary activity continues to pose a challenging question in cometary science. The activity modeling of comet 67P/Churyumov-Gerasimenko, based on data from the Rosetta mission, has significantly enhanced our comprehension of cometary activity. But thermophysical models have difficulties in simultaneously explaining the production rates of various gas species and dust. It has been suggested that different gas species might be responsible for the ejection of refractory material in distinct size ranges. This work focuses on investigating abundance and the ejection mechanisms of large ($\gtrsim$ 1 cm) aggregates from the comet nucleus. We aim to determine their properties and map the distribution of their source regions across the comet surface. This can place constraints on activity models for comets. We examined 189 images acquired at five epochs by the OSIRIS/NAC instrument. Our goal was to identify bright tracks produced by individual aggregates as they traversed the camera field of view. We generated synthetic images based on the output of dynamical simulations involving various types of aggregates. By comparing these synthetic images with the observations, we determine the characteristics of the simulated aggregates that most closely resembled the observations. We identified over 30000 tracks present in the OSIRIS images, derived constraints on the characteristics of the aggregates and mapped their origins on the nucleus surface. The aggregates have an average radius of $\simeq5$ cm, and a bulk density consistent with that of the comet's nucleus. Due to their size, gas drag exerts only a minor influence on their dynamical behavior, so an initial velocity is needed in order to bring them into the camera field of view. The source regions of these aggregates are predominantly located near the boundaries of distinct terrains on the surface.

We perform a long-term simulation of star and disk formation using three-dimensional non-ideal magnetohydrodynamics. The simulation starts from a prestellar cloud and proceeds through the long-term evolution of the circumstellar disk until $\sim 1.5\times10^5$ yr after protostar formation. The disk has size $\lesssim 50$ au and little substructure in the main accretion phase because of the action of magnetic braking and the magnetically-driven outflow to remove angular momentum. The main accretion phase ends when the outflow breaks out of the cloud, causing the envelope mass to decrease rapidly. The outflow subsequently weakens as the mass accretion rate also weakens. While the envelope-to-disk accretion continues, the disk grows gradually and develops transient spiral structures due to gravitational instability. When the envelope-to-disk accretion ends, the disk becomes stable and reaches a size $\gtrsim 300$ au. In addition, about 30% of the initial cloud mass has been ejected by the outflow. A significant finding of this work is that after the envelope dissipates, a revitalization of the wind occurs, and there is mass ejection from the disk surface that lasts until the end of the simulation. This mass ejection (or disk wind) is generated since the magnetic pressure significantly dominates both the ram pressure and thermal pressure above and below the disk at this stage. Using the angular momentum flux and mass loss rate estimated from the disk wind, the disk dissipation timescale is estimated to be $\sim10^6$ yr.

Monitoring the evolution of the anthropogenic light emissions is a priority task in light pollution research. Among the complementary approaches that can be adopted to achieve this goal stand out those based on measuring the direct radiance of the sources at ground level or from low Earth orbit satellites, and on measuring the scattered radiance (known as artificial night sky brightness or skyglow) using networks of ground-based sensors. The terrestrial atmosphere is a variable medium interposed between the sources and the measuring instruments, and the fluctuation of its optical parameters sets a lower limit for the actual source emission changes that can be confidently detected. In this paper we analyze the effect of the fluctuations of the molecular and aerosol optical depths. It is shown that for reliably detecting changes in the anthropogenic light emissions of order ~1% per year, the inter-annual variability of the annual means of these atmospheric parameters in the measurement datasets must be carefully controlled or efficiently corrected for.

We present the projected rotational velocity and molecular abundances for HD 33632 Ab obtained via Keck Planet Imager and Characterizer high-resolution spectroscopy. HD 33632 Ab is a nearby benchmark brown dwarf companion at a separation of $\sim$20 au that straddles the L/T transition. Using a forward-modeling framework with self-consistent substellar atmospheric and retrieval models for HD 33632 Ab, we derive a projected rotational velocity of 53 $\pm$ 3 km/s and water plus carbon monoxide mass fractions of log CO = $-$2.3 $\pm$ 0.3 and log H$_2$O = $-$2.7 $\pm$ 0.2. The inferred carbon-to-oxygen ratio (C/O = 0.58 $\pm$ 0.14), molecular abundances, and metallicity ([C/H] = 0.0 $\pm$ 0.2 dex) of HD 33632 Ab are consistent with its host star. Although detectable methane opacities are expected in L/T transition objects, we did not recover methane in our KPIC spectra, partly due to the high $v\sin{i}$ and to disequilibrium chemistry at the pressures we are sensitive to. We parameterize the spin as the ratio of rotation over break-up velocity, and compare HD 33632 Ab to a compilation of >200 very low-mass objects (M$\lesssim$0.1 M$_{\odot}$) that have spin measurements in the literature. There appears to be no clear trend for the isolated field low-mass objects versus mass, but a tentative trend is identified for low-mass companions and directly imaged exoplanets, similar to previous findings. A larger sample of close-in gas giant exoplanets and brown dwarfs will critically examine our understanding of their formation and evolution through rotation and chemical abundance measurements.

We explore the influence of radio-mode feedback on the properties of the cool circumgalactic medium (CGM). To this end, we assemble a statistical sample of approximately 30,000 radio galaxies with background quasars by combining optical spectroscopic measurements of luminous red galaxies (LRGs) and quasars from the year 1 dataset of Dark Energy Spectroscopic Instrument (DESI) and radio sources from the LOw-Frequency ARray Two-metre Sky Survey (LoTSS) DR2 catalog and the Very Large Array Sky Survey (VLASS) quick look catalog. Galaxies with similar optical properties but with no radio counterparts in LoTSS and VLASS are selected as the control group. We measure the cool CGM properties of radio galaxies and their control samples traced by MgII absorption lines, including covering fraction, rest equivalent width, and gas kinematics. Our results show no significant difference in the properties of gas around radio galaxies and their control sample, indicating that the operating radio-mode feedback of massive galaxies does not produce detectable effects on the properties of the cool CGM. Finally, we show that the CGM of radio galaxies contain a non-negligible amount of cool gas with approximately 10^10 solar masses. This abundance can place a stringent constraint on the radio-mode feedback models.

Fast radio bursts (FRBs) are cosmological radio transients with millisecond durations and extremely high brightness temperatures. One FRB repeater, FRB 180916.J0158+65 (FRB 180916B), was confirmed to appear 16.35-day periodic activities with 5-day activity window. Another FRB repeater, FRB 121102, and two soft gamma-ray repeaters (SGRs), SGR 1935+2154 and SGR 1806-20, also show possible periodic activities. These periodicities might originate from the precession process of young magnetars due to the anisotropic pressure from the inner magnetic fields as proposed in the literature. In this work, we analyze a self-consistent model for the rotation evolution of magnetars and obtain the evolutions of magnetar precession and obliquity. We find that if the FRB repeaters and the SGRs with (possible) periodic activities originate from the magnetar precession, their ages would be constrained to be hundreds to tens of thousands of years, which is consistent with the typical ages of magnetars. Assuming that the FRB emission is beaming in the magnetosphere as proposed in the literature, we calculate the evolution of the observable probability and the duty cycle of the active window period. We find that for a given magnetar the observable probability increases with the magnetar age in the early stage and decreases with the magnetar age in the later stage, meanwhile, there are one or two active windows in one precession period if the emission is not perfectly axisymmetric with respect to the deformation axis of a magnetar, which could be tested by the future observation for repeating FRB sources.

We present results from three-dimensional, magnetohydrodynamic, core-collapse simulations of sixteen progenitors following until 0.5 s after bounce. We use non-rotating solar-metallicity progenitor models with zero-age main-sequence mass between 9 and 24 $M_{\odot}$. The examined progenitors cover a wide range of the compactness parameter including a peak around $23 M_{\odot}$. We find that neutrino-driven explosions occur for all models within 0.3 s after bounce. We also find that the properties of the explosions and the central remnants are well correlated with the compactness. Early shock evolution is sensitive to the mass accretion rate onto the central core, reflecting the density profile of the progenitor stars. The most powerful explosions with diagnostic explosion energy $E_{\rm exp} \sim 0.75 \times 10^{51}$ erg are obtained by 23 and 24 $M_{\odot}$ models, which have the highest compactness among the examined models. These two models exhibit spiral SASI motions during 150-230 ms after bounce preceding a runaway shock expansion and leave a rapidly rotating neutron star with spin periods $\sim 50$ ms. Our models predict the gravitational masses of the neutron star ranging between $1.22 M_{\odot}$ and $1.67 M_{\odot}$ and their spin periods 0.04-4 s. The number distribution of these values roughly matches observation. On the other hand, our models predict small hydrodynamic kick velocity (15-260 km/s), although they are still growing at the end of our simulations. Further systematic studies, including rotation and binary effects, as well as long-term simulations up to several seconds, will enable us to explore the origin of various core-collapse supernova explosions.

We present six nearly full-sky maps made from data taken by radiometers on the Juno satellite during its 5-year flight to Jupiter. The maps represent integrated emission over $\sim 4\%$ passbands spaced approximately in octaves between 600 MHz and 21.9 GHz. Long time-scale offset drifts are removed in all bands, and, for the two lowest frequency bands, gain drifts are also removed from the maps via a self-calibration algorithm similar to the NPIPE pipeline used by the Planck collaboration. We show that, after this solution is applied, residual noise in the maps is consistent with thermal radiometer noise. We verify our map solutions with several consistency tests and end-to-end simulations. We also estimate the level of pixelization noise and polarization leakage via simulations.

Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. We intend to understand which type of regions are probed by sulphur-bearing species and which process(es) dominate(s) the release of sulphur into the gas phase. We used the high-angular resolution (1.6" or 27 pc) observations from the ALCHEMI ALMA Large Program to image several sulphur-bearing species towards the central molecular zone (CMZ) of NGC 253. We performed local thermodynamic equilibrium (LTE) and non-LTE large velocity gradient (LVG) analyses to derive the physical conditions of the gas in which S-bearing species are emitted, and their abundance ratios across the CMZ. Finally, we compared our results with previous ALCHEMI studies and a few selected Galactic environments. We found that not all sulphur-bearing species trace the same type of gas: strong evidence indicates that H2S and part of the emission of OCS, H2CS, and SO, are tracing shocks whilst part of SO and CS emission rather trace the dense molecular gas. For some species, such as CCS and SO2, we could not firmly conclude on their origin of emission. The present analysis indicates that the emission from most sulphur-bearing species throughout the CMZ is likely dominated by shocks associated with ongoing star formation. In the inner part of the CMZ where the presence of super star clusters was previously indicated, we could not distinguish between shocks or thermal evaporation as the main process releasing the S-bearing species.

We develop a model-independent approach to lagrangian perturbation theory for the large scale structure of the universe. We focus on the displacement field for dark matter particles, and derive its most general structure without assuming a specific form for the equations of motion, but implementing a set of general requirements based on symmetry principles and consistency with the perturbative approach. We present explicit results up to sixth order, and provide an algorithmic procedure for arbitrarily higher orders. The resulting displacement field is expressed as an expansion in operators built up from the linear density field, with time-dependent coefficients that can be obtained, in a specific model, by solving ordinary differential equations. The derived structure is general enough to cover a wide spectrum of models beyond $\Lambda$CDM, including modified gravity scenarios of the Hordenski type and models with multiple dark matter species. This work is a first step towards a complete model-independent lagrangian forward model, to be employed in cosmological analyses with power spectrum and bispectrum, other summary statistics, and field-level inference.

Context. A necessary ingredient in understanding the star formation history of a young cluster is knowledge of the extinction towards the region. This has typically been done by making use of the colour-difference method with photometry, or similar methods utilising the colour-colour diagram. Aims. The colour-excess can be independently determined by studying the decrements of the recombination lines produced by the nebular gas. Having access to many recombination lines from the same spectral series removes the need of adopting an extinction curve. Methods. Using the Micro-Shutter Assembly (MSA) on board the Near InfraRed Spectrograph (NIRSpec), multi-object spectroscopy was performed, yielding 600 nebular spectra from the Galactic massive star formation region NGC 3603. The recombination line intensity ratios were used to determine independent values of colour-excess. Results. The extinction characteristics of NGC 3603 are similar to other Galactic HII regions like Orion, as well as starburst regions such 30 Doradus in the Large Magellanic Cloud, in that we find a relatively large value of total-to-selective extinction of 4.8 $\pm$ 1.06, larger than the Galactic average of 3.1. We find a typical value for the colour-excess of 0.64 $\pm$ 0.27, significantly lower than values determined in previous studies. We also present a stacked nebular spectrum with a typical continuum S/N = 70. This spectrum highlights the recombination lines of the HII region, several s-process elements such as Kr III and Se IV, and molecular H2 emission lines. Conclusions. Using ratios of hydrogen recombination lines, we calculated the total-to-selective extinction, colour-excess and visual extinction for > 200 lines of sight across NGC 3603. An extinction curve with a total-to-selective extinction of 4.8 $\pm$ 1.06 was found, corresponding to a colour-excess of 0.64 $\pm$ 0.27.

Planets form from protoplanetary discs of dust and gas that surround stars younger than a few million years. The properties of these discs dictate how planets grow, determining the nature, architecture and composition of the observed exoplanet population. The nature of (turbulent) angular momentum transport through the disc in particular is critical for planet growth and migration. While the dominant physical processes in disc evolution remain unknown, theoretical studies widely assume an isolated star-disc system. Here we challenge this conventional assumption by showing that, far from being isolated systems, discs in typical star forming environments are constantly replenished by capturing new material from the interstellar medium (ISM). We show that this in-fall alone explains observed disc masses, outer radii and stellar accretion rates as a function of time and stellar mass. ISM capture is also capable of driving disc turbulence corresponding to viscosity in the range $\alpha_\mathrm{SS} \sim 10^{-5}{-}10^{-1}$, as observationally inferred. We find that $20{-}70$~percent of discs are composed of material captured in the most recent half of their life-time, implying their properties are not direct probes of internal disc physics. Our results suggest that planet formation is driven by the turbulent cascade of energy from galactic-scales down to the protoplanetary disc and that recent evidence of in-fall from the ISM onto mature protoplanetary discs are part of an important, general process. This represents a far-reaching shift in our understanding of every stage of planet formation, from dust aggregation to the migration of giant planets through the disc.

Light mass warm dark matter is an interesting and viable alternative to the cold dark matter paradigm. An intriguing variation of this scenario is the mass-varying dark matter model where the dark matter mass varies with time during its cosmic history. This is realized in multiple particle physics models. In this work, we study the cosmological constraints on such a model where the dark matter mass transitions from zero to a finite value in the early Universe. In this model, the matter power spectrum exhibits power suppression below a scale that depends on the epoch of transition, and the angular power spectrum of the cosmic microwave background show a distinctive phase shift. We use the latest cosmic microwave background and the weak lensing data to place lower limit on the transition redshift and ease the $S_8$ tension, unlike the warm dark matter model. This analysis also facilitates a marginal detection of the dark matter (DM) mass. Our findings reveal that while Planck data alone reduces the $S_8$ tension to approximately $2\sigma$, it does not sufficiently constrain the DM mass. However, when combined with the $S_8$ measurement from KIDS1000+BOSS+2dfLenS, the tension significantly decreases to roughly $1.3\sigma$, and we observe the detection of a DM mass at $41.7^{+7.81}_{-27.5}\,\mathrm{eV}$. Further analysis incorporating a combined data set from ACT and weak lensing results in an even more pronounced reduction in the tension to approximately $0.4\sigma$, alongside a higher detected mass of $51.2^{+16}_{-33.5}\,\mathrm{eV}$. We also find a better fit to the combined data compared to the $\Lambda$CDM model.

We present the results from a long term X-ray analysis of Mrk 279 during the period 2018-2020. We use data from multiple missions - AstroSat, NuSTAR and XMM-Newton, for the purpose. The X-ray spectrum can be modelled as a double Comptonisation along with the presence of neutral Fe K${\alpha}$ line emission, at all epochs. We determined the source's X-ray flux and luminosity at these different epochs. We find significant variations in the source's flux state. We also investigated the variations in the source's spectral components during the observation period. We find that the photon index and hence the spectral shape follow the variations only over longer time periods. We probe the correlations between fluxes of different bands and their photon indices, and found no significant correlations between the parameters.

Rapidly rotating neutron stars are similar to highly flattened ellipsoids. Observed spectra of flattened stars must exhibit effects of non spherical shape and the gravitational darkening. We examined in detail the influence of both effects on the observed central energies and profiles of lines of highly ionized iron, FeXXV and FeXXVI. We note that the gravitational darkening effect do not change the central energy of lines. What is important, spectra of neutron stars that rotate with different frequencies and are seen at various inclinations angles differ significantly. The appearance and the depth of lines strongly depends on the parameters like the inclination angle of the star or the frequency of the star rotation. In this paper we clearly show that the gravitational darkening effect should be included in realistic model of the atmospheres of the neutron stars.

Episodic accretion in protostars leads to luminosity outbursts that end up heating their surroundings. This rise in temperature pushes the snow lines back, enabling the desorption of chemical species from dust grain surfaces, which may significantly alter the chemical history of the accreting envelope. However, a limited number of extensive chemical surveys of eruptive young stars have been performed thus far. In the present study, we carry out a large spectral survey of the binary Class I protostar L1551 IRS 5, known to be a FUor-like object, in the 3mm and 2mm bands with the IRAM-30m telescope. As a result, we detected more than 400 molecular lines. The source displays a great chemical richness with the detection of 75 species, including isotopologues. Among these species, there are 13 hydrocarbons, 25 N-bearing species, 30 O-bearing species, 15 S-bearing species, 12 deuterated molecules, and a total of 10 complex organic molecules (l-C4H2, CH3CCH, CH2DCCH, CH3CHO, CH3CN, CH3OCH3, CH3OCHO, CH3OH, CH2DOH, and HC5N). With the help of local thermodynamic equilibrium (LTE) and non-LTE models, we determined the column densities of most molecules as well as excitation and kinetic temperatures. While most of those molecules trace the cold envelope (. 20 K), the OCS and CH3OH emission arise from the warm (> 100 K) innermost (< 2'' ) regions. We compared the chemical inventory of L1551 IRS 5 and its column density ratios, including isotopic ratios, with other protostellar sources. A broad chemical diversity is seen among Class I objects. More observations with both single-dish telescopes and interferometers are needed to characterize the diversity in a larger sample of protostars, while more astrochemical models would help explain this diversity, in addition to the impact of luminosity outbursts on the chemistry of protostellar envelopes.

Cosmological inflation is a popular paradigm for understanding Cosmic Microwave Background Radiation (CMBR); however, it faces many conceptual challenges. An alternative mechanism to inflation for generating an almost scale-invariant spectrum of perturbations is a \emph{bouncing cosmology} with an initial matter-dominated contraction phase, during which the modes corresponding to currently observed scales exited the Hubble radius. Bouncing cosmology avoids the initial singularity but has fine-tuning problems. Taking an \emph{agnostic view} of the two early-universe paradigms, we propose a quantum measure -- Dynamical Fidelity Susceptibility (DFS) of CMBR -- that distinguishes the two scenarios. Taking two simple models with the same power-spectrum, we explicitly show that DFS behaves differently for the two scenarios. We discuss the possibility of using DFS as a distinguisher in the upcoming space missions.

We study the mass distribution of galaxy clusters in Milgromian dynamics, or modified Newtonian dynamics (MOND). We focus on five galaxy clusters from the X-COP sample, for which high-quality data are available on both the baryonic mass distribution (gas and stars) and internal dynamics (from the hydrostatic equilibrium of hot gas and the Sunyaev-Zeldovich effect). We confirm that galaxy clusters require additional `missing matter' in MOND, although the required amount is drastically reduced with respect to the non-baryonic dark matter in the context of Newtonian dynamics. We studied the spatial distribution of the missing matter by fitting the acceleration profiles of the clusters with a Bayesian method, finding that a physical density profile with an inner core and an outer $r^{-4}$ decline (giving a finite total mass) provide good fits within $\sim$1 Mpc. At larger radii, the fit results are less satisfactory but the combination of the MOND external field effect and hydrostatic bias (quantified as 10$\%$-40$\%$) can play a key role. The missing mass must be more centrally concentrated than the intracluster medium (ICM). For relaxed clusters (A1795, A2029, A2142), the ratio of missing-to-visible mass is around $1-5$ at $R\simeq200-300$ kpc and decreases to $0.4-1.1$ at $R\simeq2-3$ Mpc, showing that the total amount of missing mass is smaller than or comparable to the ICM mass. For clusters with known merger signatures (A644 and A2319), this global ratio increases up to $\sim$5 but may indicate out-of-equilibrium dynamics rather than actual missing mass. We discuss various possibilities regarding the nature of the extra mass, in particular `missing baryons' in the form of pressure-confined cold gas clouds with masses of $<10^5$ M$_\odot$ and sizes of $< 50$ pc.

We have assembled a sample of $\sim$8200 stars with spectral types F5V-M5V, all having directly measured X-ray luminosities from eROSITA and rotation periods from TESS, and having empirically estimated ages via their membership in stellar clusters and groups identified in Gaia astrometry (ages 3-500 Myr). This is the largest such study sample yet assembled for the purpose of empirically constraining the evolution of rotationally driven stellar X-ray activity. We observe rotation-age-activity correlations that are qualitatively as expected: stars of a given spectral type spin down with age and they become less X-ray active as they do so. We provide simple functional representations of these empirical relationships that predict X-ray luminosity from basic observables to within 0.3 dex. Interestingly, we find that the rotation-activity relationship is far simpler and more monotonic in form when expressed in terms of stellar angular momentum instead of rotation period. We discuss how this finding may relate to the long-established idea that rotation-activity relationships are mediated by stellar structure (e.g., convective turnover time, surface area). Finally, we provide an empirical relation that predicts stellar angular momentum from basic observables, and without requiring a direct measurement of stellar rotation, to within 0.5 dex.

We report the study of a huge optical intraday flare on November 12, 2021, at 2 am UT, in the blazar OJ287. In the binary black hole model it is associated with an impact of the secondary black hole on the accretion disk of the primary. Our multifrequency observing campaign was set up to search for such a signature of the impact, based on a prediction made eight years earlier. The first I-band results of the flare have already been reported by \cite{2024ApJ...960...11K}. Here we combine these data with our monitoring in the R-band. There is a big change in the R-I spectral index by $1.0\pm0.1$ between the normal background and the flare, suggesting a new component of radiation. The polarization variation during the rise of the flare suggests the same. The limits on the source size place it most reasonably in the jet of the secondary black hole. We then ask why we have not seen this phenomenon before. We show that OJ287 was never before observed with sufficient sensitivity on the night when the flare should have happened according to the binary model. We also study the probability that this flare is just an oversized example of intraday variability, using the Krakow-dataset of intense monitoring between 2015 and 2023. We find that the occurrence of a flare of this size and rapidity is unlikely. In the Appendix, we give the full orbit-linked historical light curve of OJ287 as well as the dense monitoring sample of Krakow.

The measurement of pulsar pulse times-of-arrival (ToAs) is a crucial step in detecting low-frequency gravitational waves. To determine ToAs, we can use template-matching to compare each observed pulse profile with a standard template. However, using different combinations of templates and template-matching methods (TMMs) without careful consideration may lead to inconsistent results. In pulsar timing array (PTA) experiments, distinct ToAs from the same observations can be obtained, due to the use of diverse templates and TMMs. In other words, employing diverse approaches can yield different timing results and would thus have a significant impact on subsequent gravitational wave searches. In this paper, we examine several commonly used combinations to analyze their effect on pulse ToAs. we evaluate the potential impact of template and TMM selection on thirteen typical millisecond pulsars within the European PTA. We employ pulsar timing methods, specifically the root mean square and reduced chi-square $\chi_r^2$ of the residuals of the best timing solution to assess the outcomes. Additionally, we evaluate the system-limited noise floor (SLNF) for each pulsar at various telescopes operating around 1.4~GHz using frequency-resolved templates. Our findings suggest that utilizing data-derived and smoothed templates in conjunction with the Fourier-domain with Markov-chain Monte Carlo (FDM) TMM is generally the most effective approach, though there may be exceptions that require further attention. Furthermore, we determine that pulse phase jitter noise does not significantly limit the current precision of the European PTA's timing, as jitter levels derived from other studies are much smaller than the SLNF.

The study examines the heating profile of hot solar transition region loops, particularly focusing on transient brightenings observed in IRIS 1400Å slit-jaw images. The findings challenge the adequacy of simplistic, singular heating mechanisms, revealing that the heating is temporally impulsive and requires a spatially complex profile with multiple heating scales. A forward modeling code is utilized to generate synthetic IRIS emission spectra of these loops based on HYDRAD output, confirming that emitting ions are out of equilibrium. The modeling further indicates that density-dependent dielectronic recombination rates must be included to reproduce the observed line ratios. Collectively, this evidence substantiates that the loops are subject to impulsive heating and that the components of the transiently brightened plasma are driven far from thermal equilibrium. Heating events such as these are ubiquitous in the transition region and the analysis described above provides a robust observational diagnostic tool for characterizing the plasma.

The linear relationship between pulsar micro-pulse widths and rotation period is consistent with the existence of a physical length L on the neutron-star surface and seen on the observer arc of transit across the polar cap. Within the ion-proton model it is the width of the minimum area of surface that can support the critical growth rate of the unstable two-beam Langmuir mode that is the source of the radio emission.

We use spectro-polarimetric data recorded by Hinode to analyze the magnetic field configuration of a part of a sunspot (AR 10923) where a bundle of penumbral filaments are intruding into its umbra. We want to explore the role of the sunspot magnetic configuration in the formation and kinematics of the fine-structures, such as umbral dots and light bridges, inside the sunspot umbra. Both direct inferences from polarization Stokes profiles and the inversion results using the SIR code indicate an aligned magnetic field configuration in the umbra where moving umbral dots are easily formed at the leading edges of the rapidly intruding penumbral filaments. We suggest that the magnetic field topology is rearranged leading to the observed aligned magnetic field lines via magnetic reconnection process by which a part of the magnetic energy is converted into thermal and kinetic energy. This new configuration causes the umbral fine-structures to form easily and more frequently.

NASA's Meteoroid Engineering Model (MEM) is designed to provide aerospace engineers with an accurate description of potentially hazardous meteoroids. It accepts a spacecraft trajectory as input and its output files describe the flux, speed, directionality, and density of microgram- to gram-sized meteoroids relative to the provided trajectory. MEM provides this information at a fairly fine level of detail in order to support detailed risk calculations. However, engineers and scientists in the very early planning stages of a mission may not yet have developed a trajectory or acquired the tools to analyze environment data. Therefore, we have developed an online library of sample MEM runs that allow new users or overloaded mission planners to get a quick feel for the characteristics of the meteoroid environment. This library provides both visualizations of these runs and input files that allow users to replicate them exactly. We also discuss the number of state vectors needed to obtain an accurate representation of the environment encountered along our sample trajectories, and outline a process for verifying that any given trajectory is adequately sampled.

We investigate the alignment mechanism between the magnetic field and interstellar clouds formed by the collision of warm atomic gas. We find that the magnetic field, initially oriented parallel to the flow, is perturbed by a fast MHD shock, which amplifies magnetic field fluctuations parallel to the shock front. Behind the shock, the compressive downstream velocity field further amplifies the magnetic field component parallel to the shock front. This mechanism causes the magnetic field to become increasingly parallel to the dense layer, and the development of a shear flow around the latter. Furthermore, the bending-mode perturbations on the dense layer are amplified by the non-linear thin-shell instability (NTSI), stretching the density structures formed by the thermal instability, and rendering them parallel to the bent field lines. By extension, we suggest that a tidal stretching velocity gradient such as that produced in gas infalling into a self-gravitating structure must straighten the field lines along the accretion flow, orienting them perpendicular to the density structures. We also find that the upstream superalfvénic regime transitions to a transalfvénic regime between the shock and the condensation front, and then to a subalfvénic regime inside the condensations. Our results provide a feasible physical mechanism for the observed transition from parallel to perpendicular relative orientation of the magnetic field and the density structures as the density structures become increasingly dominated by self-gravity.

Deep learning (DL) techniques are a promising approach among the set of methods used in the ever-challenging determination of stellar parameters in M dwarfs. In this context, transfer learning could play an important role in mitigating uncertainties in the results due to the synthetic gap (i.e. difference in feature distributions between observed and synthetic data). We propose a feature-based deep transfer learning (DTL) approach based on autoencoders to determine stellar parameters from high-resolution spectra. Using this methodology, we provide new estimations for the effective temperature, surface gravity, metallicity, and projected rotational velocity for 286 M dwarfs observed by the CARMENES survey. Using autoencoder architectures, we projected synthetic PHOENIX-ACES spectra and observed CARMENES spectra onto a new feature space of lower dimensionality in which the differences between the two domains are reduced. We used this low-dimensional new feature space as input for a convolutional neural network to obtain the stellar parameter determinations. We performed an extensive analysis of our estimated stellar parameters, ranging from 3050 to 4300 K, 4.7 to 5.1 dex, and -0.53 to 0.25 dex for Teff, logg, and [Fe/H], respectively. Our results are broadly consistent with those of recent studies using CARMENES data, with a systematic deviation in our Teff scale towards hotter values for estimations above 3750 K. Furthermore, our methodology mitigates the deviations in metallicity found in previous DL techniques due to the synthetic gap. We consolidated a DTL-based methodology to determine stellar parameters in M dwarfs from synthetic spectra, with no need for high-quality measurements involved in the knowledge transfer. These results suggest the great potential of DTL to mitigate the differences in feature distributions between the observations and the PHOENIX-ACES spectra.

Characterizing the exchange of water between the Martian atmosphere and the (sub)surface is a major challenge for understanding the mechanisms that regulate the water cycle. Here we present a new dataset of water ice detected on the Martian surface with the Thermal Emission Imaging System (THEMIS). The detection is based on the correlation between bright blue-white patterns in visible images and a temperature measured in the infrared that is too warm to beassociated with CO2 ice and interpreted instead as water ice. Using this method, we detect ice down to 21.4°S, 48.4°N, on the pole-facing slopes at mid-latitudes, and on any surface orientation poleward of 45° latitude. Water ice observed with THEMIS is most likely seasonal rather than diurnal. Our dataset is consistent with near-infrared spectroscopic data predictions by the Mars Planetary Climate Model. The water frost average temperature is 170 K, and the maximum temperature measured is 243 K, lower than the water ice melting point. We show that the melting of pure water ice on the surface is unlikely due to cooling by latent heat during its sublimation. However, 243 THEMIS images show frosts that are hot enough to form brines if salts are present on the surface. The water vapor pressure at the surface, calculated from the ice temperature, indicates a dry atmosphere in early spring, during the recession of the CO2 ice cap. When it sublimes, the frost acts as a vapor source that is wetter than the near-surface atmosphere, which stabilizes the subsurface ice.

We present the design of a portable coronagraph, CATEcor, that incorporates a novel "shaded truss" style of external occultation and serves as a proof-of-concept for that family of coronagraphs. The shaded truss design style has the potential for broad application in various scientific settings. We conceived CATEcor itself as a simple instrument to observe the corona during the darker skies available during a partial solar eclipse, or for students or interested amateurs to detect the corona under ideal non-eclipsed conditions. CATEcor is therefore optimized for simplicity and accessibility to the public. It is implemented using an existing dioptric telescope and an adapter rig that mounts in front of the objective lens, restricting the telescope aperture and providing external occultation. The adapter rig, including occulter, is fabricated using fusion deposition modeling (FDM; colloquially "3D printing"), greatly reducing cost. The structure is designed to be integrated with moderate care and may be replicated in a university or amateur setting. While CATEcor is a simple demonstration unit, the design concept, process, and trades are useful for other more sophisticated coronagraphs in the same general family, which might operate under normal daytime skies outside the annular-eclipse conditions used for CATEcor.

Observations of velocity dispersions of galactic structures over a wide range of scales point to the existence of a universal acceleration scale $a_0\sim 10^{-10}$ m/s$^2$. Focusing on the fuzzy dark matter paradigm, which proposes ultralight dark matter with mass around $10^{-22}$ eV and de Broglie wavelength $\lambda\sim {\rm few}\times10^{2}$ parsecs, we highlight the emergence of the observed acceleration scale from quantum effects in a fluid-like description of the dark matter dynamics. We then suggest the possibility of a natural connection between the acceleration scale and dark energy within the same paradigm.

Near-solar mass black holes (BHs) could have been involved in the two recent gravitational wave events, GW190425 and GW190814. Since such a low mass BH cannot be formed via stellar evolution, a model has been proposed based on the core collapse of a neutron star initiated by a certain number of dark matter (DM) particles. In this process, the accumulated DM particles collapse to form a tiny BH inside the neutron star, and the entire neutron star is transmuted into a BH after a certain time due to the accretion of matter by the endoparasitic BH from its host. Here, we argue that, depending on the initial conditions, a dark core collapse could give rise to either a BH or a naked singularity. For example, if the accumulated cloud of DM particles in the core of a neutron star can be modeled as an anisotropic fluid and it fulfils the criterion for collapse, an endoparasitic naked singularity could form instead of an endoparasitic BH. Immediately after its formation, the naked singularity should begin accreting matter from the host neutron star, thus eventually transmuting the entire host into a near-solar mass, relatively slowly-spinning naked singularity. We also propose a general technique to constrain the DM particle--neutron scattering cross section using the lack of pulsars near the Galactic centre and assuming that these missing pulsars have already been transmuted into BHs and/or naked singularities. Thus, the missing pulsars also indicate the existence of many such singularities near the Galactic center.

Previous X-ray works on Mrk 1239 have revealed a complex Narrow Line Seyfert 1 (NLS1) that exhibits substantial absorption and strong emission from both collisional (CIE) and photoionized (PIE) plasmas. Here, we report on deep-pointed observations with $XMM{\rm -}Newton$ and $NuSTAR$, along with $Swift$ monitoring, to understand the $0.3-30$ keV continuum emission and the central engine geometry. A strong X-ray flare, where the AGN brightens by a factor of five in $\sim30$ ks, is captured between $4-30$ keV and can be attributed to a brightening of the primary continuum. However, the lack of any variability below $\sim3$ keV on long- or short-time scales requires complete absorption of the AGN continuum with a neutral medium of column density $\sim 10^{23.5}{\rm cm}^{-2}$. The timing and spectral properties are consistent with a blurred reflection interpretation for the primary emission. The variability and presence of a Compton hump disfavours ionized partial covering. The neutral absorber, if outflowing, could be crashing into the surrounding medium and ISM to produce the low-energy continuum and CIE. Scattered emission off the inner torus could produce the PIE. The intricate scenario is demanded by the data and highlights the complexity of the environment that is normally invisible when overwhelmed by the AGN continuum. Objects like Mrk 1239 serve as important sources for unveiling the interface between the AGN and host galaxy environments.

Extreme-mass-ratio inspirals consist of binary systems of compact objects, with orders of magnitude differences in their masses, in the regime where the dynamics are driven by gravitational wave emission. The unique nature of extreme-mass-ratio inspirals facilitates the exploration of diverse and stimulating tests related to various aspects of black holes and General Relativity. These tests encompass investigations into the spacetime geometry, the dissipation of the binary system energy and angular momentum, and the impact of the intrinsic non-linear effects of the gravitational interaction. This review accounts for the manifold opportunities for testing General Relativity, ranging from examinations of black hole properties to cosmological and multimessenger tests.

In the present study, we generalize the possible ghost field configurations within the framework of $k$-essence theory to the Simpson-Visser metric area function $\Sigma^2=x^2+a^2$. Our analysis encompasses field configurations for the region-defined metric function $dA_\pm$ as well as the general solution that asymptotically behaves as Schwarzschild-de Sitter for $x\to-\infty$. Specifically, we investigate two scalar field configurations and define the associated potential for each one. Through rigorous calculations, we verify that all equations of motion are satisfied. Notably, our findings indicate that even when proposing new configurations of ghost scalar fields, the energy conditions remain unchanged. This result serves to validate the wormhole solutions obtained in previous studies.

In this work, we have reconstructed the extended $f(\mathcal{P})$ cubic gravity and symmetric $f(\mathcal{Q})$ teleparallel gravity from $(m,n)$-type Barrow Holographic Dark Energy (BHDE) and find the unknown functions $f(\mathcal{P})$ and $f(\mathcal{Q})$ in terms of $\mathcal{P}$ and $\mathcal{Q}$ by taking the universe to be flat, homogeneous and isotropic. We then analyzed the behavior and stability of each model for the entire stages of the evolution of the universe by studying several important parameters such as the deceleration parameter, equation of state (EoS) parameter $\omega_{DE}$, square of the speed of sound $v_s^2$. Apart from this, we have studied the cosmographic behavior by plotting the jerk parameter, snap parameter, and lerk parameter against the redshift. We have also examined the $\omega'_{DE}-\omega_{DE}$ phase plane and $(r,s^*)$, $(r,q)$ statefinder parameters that provide valuable insights into the dynamics of the universe and the distinctive features of the dark energy. All these analyses pointed out that our model can produce a universe going through an accelerated expansion with the quintessence type dark energy.

Einstein's Field Equations have proven applicable across many scales, from black holes to cosmology. Even the mysterious Cosmological Constant found a physical interpretation in the so-called ``dark energy'' causing the accelerated cosmic expansion as inferred from multiple observables. Yet, we still lack a material source for this dark fluid. Probing the local universe to find it yields complementary information to the one from the cosmic microwave background. Could dark energy be sourced by super-extremal charged black holes? Contrary to intuition, such objects could exist with only weak observational signatures. The latter are introduced here to outline how sky surveys can identify individual candidates which challenge Cosmic Censorship on the one hand but may explain the physical origin of the Cosmological Constant on the other.

Cosmic Microwave Background (CMB) photons can undergo resonant conversion into axions in the presence of magnetized plasma distributed inside non-linear large-scale structure (LSS). This process leads to axion-induced patchy screening: secondary temperature and polarization anisotropies with a characteristic non-blackbody frequency dependence that are strongly correlated with the distribution of LSS along our past light cone. We compute the axion-induced patchy screening contribution to two- and three- point correlation functions that include CMB anisotropies and tracers of LSS within the halo model. We use these results to forecast the sensitivity of existing and future surveys to photon-axion couplings for axion masses between $2\times 10^{-13}$ eV and $3\times 10^{-12}$ eV, using a combination of empirical estimates from Planck data of the contribution from instrumental noise and foregrounds as well as modeled contributions on angular scales only accessible with future datasets. We demonstrate that an analysis using Planck and the unWISE galaxy catalogue would be complementary to the most sensitive existing astrophysical axion searches, probing couplings as small as $3\times 10^{-12} \, {\rm GeV}^{-1}$, while observations from a future survey such as CMB-S4 could extend this reach by almost an additional order of magnitude.

We treat the guiding-center dynamics in a non-uniform external Maxwell field using a manifestly Lorentz covariant action principle which easily reproduces the known Vandervoort equations of motion. We derive the corresponding kinetic theory and ideal hydrodynamic theory. In contrast to conventional five-equation hydrodynamics, the guiding-center hydrodynamics needs only three equations due to a constraint on the motion across magnetic field. We argue that this hydrodynamics applies more generally than the kinetic theory, e.g., for strongly-coupled quark-gluon plasma.

We derive a novel BPS bound from chiral perturbation theory minimally coupled to electrodynamics at finite isospin chemical potential. At a critical value of the isospin chemical potential, a system of three first-order differential field equations (which implies the second-order field equations) for the gauge field and the hadronic profile can be derived from the requirement to saturate the bound. These BPS configurations represent magnetic multi-vortices with quantized flux supported by a superconducting current. The corresponding topological charge density is related to the magnetic flux density, but is screened by the hadronic profile. Such a screening effect allows to derive the maximal value of the magnetic field generated by these BPS magnetic vortices, being $B_{max}=2,04 \times 10^{14}G$. The solution for a single BPS vortex is discussed in detail, and some physical consequences, together with the comparison with the magnetic vortices in the Ginzburg-Landau theory at critical coupling, are described.

We consider the single-parameter $\mathcal{R}+ c\mathcal{R}^2$ gravitational action and use constraints from astrophysics and the laboratory to derive a natural relation between the coefficient $c$ and the value of the cosmological constant. We find that the renormalisation of $c$ from the energy of the inflationary phase to the infrared, where the acceleration of the expansion of the Universe takes place, is correlated with the evolution of the vacuum energy. Our results suggest that the coefficient of the $\mathcal{R}^2$ term may provide an unexpected bridge between high-energy physics and cosmological phenomena such as inflation and dark energy.

Low Level RF (LLRF) control systems of linear accelerators (LINACs) are typically implemented with heterodyne based architectures, which have complex analog RF mixers for up and down conversion. The Gen 3 Radio Frequency System-on-Chip (RFSoC) device from AMD Xilinx integrates data converters with maximum RF frequency of 6~GHz. This enables direct RF sampling of C-band LLRF signal typically operated at 5.712 GHz without any analogue mixers, which can significantly simplify the system architecture. The data converters sample RF signals in higher order Nyquist zones and then up or down convert digitally by the integrated data path in RFSoC. The closed-loop feedback control firmware implemented in FPGA integrated in RFSoC can process the base-band signal from the ADC data path and calculate the updated phase and amplitude to be up-mixed by the DAC data path. We have developed a C-band LLRF control RFSoC platform with direct RF sampling, which targets Cool Copper Collider (C3) and other C or S band LINAC research and development projects. In this paper, the architecture of the platform will be described. We have optimized the configuration of the data converter and characterized performance of them with RF pulses. The test results for some of the key performance parameters for the LLRF platform with our custom solid-state amplifier, such as phase and amplitude stability, will be discussed in this paper.

Extracting the rotational energy from a Kerr black hole (BH) is one of the crucial topics in relativistic astrophysics. Here, we give special attention to the Penrose ballistic process based on the fission of a massive particle $\mu_0$ into two particles $\mu_1$ and $\mu_2$, occurring in the ergosphere of a Kerr BH. Bardeen et al. indicated that for the process to occur, some additional "{\it hydrodynamical forces or superstrong radiation reactions}" were needed. This idea was further expanded by Wald and Chandrasekhar. This animosity convinced T. Piran and collaborators to move from a simple three-body system characterizing the original Penrose process to a many-body system. This many-body approach was further largely expanded by others, introducing additional processes, some questionable in their validity. In this letter, we return to the simplest original Penrose process and show that the solution of the equations of motion, imposing the turning point condition on their trajectories, leads to the rotational energy extraction from the BH expected by Penrose. The efficiency of energy extraction by a single process is precisely quantified for three different single decay processes occurring respectively at $r=1.2 M$, $r=1.5 M$, and $r=1.9 M$. An interesting repetitive model has been proposed by Misner, Thorne \& Wheeler. Indeed, it would appear that a repetitive sequence of $246$ decays of the above injection process at $r=1.2 M$ and the corresponding ones at $r=1.5 M$ and $r=1.9 M$ could extract $100\%$ of the rotational energy of the BH. The accompanying paper shows that accounting for the irreducible mass introduces a non-linear approach that avoids violating energy conservation, leading to new energy extraction processes and demonstrating the impossibility of extracting the whole BH rotational energy by a sequence of Penrose processes.

The recent observations of the galactic center of the M87 galaxy have made the field of observing black holes and calculating its shadow much more intriguing. Approaching the question of calculating shadows, many approximations are made in order to simplify the equations which makes the considered case less realistic. Understanding the shadow of different singularities under the influence of magnetic field is of more importance astrophysically as the accreting matter around the singularity would generate electromagnetic fields as it would be in plasma state due to the high tidal effects. Here, we use Ernst technique to immerse spacetimes in uniform, sourceless and asymptotic magnetic field. Later, we compare the effective potential of null geodesics in magnetized and non-magnetized cases. This study would be helpful in understanding the M87 shadow and the forthcoming image of shadow of SagA*.

State estimation of nonlinear dynamical systems has long aimed to balance accuracy, computational efficiency, robustness, and reliability. The rapid evolution of various industries has amplified the demand for estimation frameworks that satisfy all these factors. This study introduces a neuromorphic approach for robust filtering of nonlinear dynamical systems: SNN-EMSIF (spiking neural network-extended modified sliding innovation filter). SNN-EMSIF combines the computational efficiency and scalability of SNNs with the robustness of EMSIF, an estimation framework designed for nonlinear systems with zero-mean Gaussian noise. Notably, the weight matrices are designed according to the system model, eliminating the need for a learning process. The framework's efficacy is evaluated through comprehensive Monte Carlo simulations, comparing SNN-EMSIF with EKF and EMSIF. Additionally, it is compared with SNN-EKF in the presence of modeling uncertainties and neuron loss, using RMSEs as a metric. The results demonstrate the superior accuracy and robustness of SNN-EMSIF. Further analysis of runtimes and spiking patterns reveals an impressive reduction of 85% in emitted spikes compared to possible spikes, highlighting the computational efficiency of SNN-EMSIF. This framework offers a promising solution for robust estimation in nonlinear dynamical systems, opening new avenues for efficient and reliable estimation in various industries that can benefit from neuromorphic computing.

The $\gamma$ ray emission originating from in-flight annihilation (IA) of positrons is a powerful observable for constraining high-energy positron production from exotic sources. By comparing diffuse $\gamma$ ray observations of INTEGRAL, COMPTEL and EGRET to theoretical predictions, we set the most stringent constraints on electrophilic feebly interacting particles (FIPs), thereby proving IA as a valuable probe of new physics. In particular, we extensively discuss the case of MeV-scale sterile neutrinos, where IA sets the most stringent constraints, excluding $|U_{\mu4}|^{2} \gtrsim 10^{-13}$ and $|U_{\tau4}|^{2} \gtrsim 2\times 10^{-13}$ for sterile neutrinos mixed with $\mu$ and $\tau$ neutrinos respectively. These constraints improve existing limits by more than an order of magnitude. We briefly discuss the application of these results to a host of exotic positron sources such as dark photons, axion-like particles, primordial black holes (PBHs) and sub-GeV dark matter (DM).

Cherenkov-type particle detectors or scintillators use as a fundamental element photomultiplier tubes, whose efficiency decreases when subjected to the Earth's magnetic field. This work develops a geomagnetic field compensation system based on coils for large scale cylindrical detectors. The effect of different parameters such as the size of the detector, the distance between coils or the magnetic field strength on the compensation using a basic coil system composed of circular and rectangular coils is studied. The addition of coils of very specific geometry and position to the basic configuration is proposed in order to address the compensation in the areas of the detector where it is more difficult to influence, in order to minimize the loss of efficiency. With such improvement, in the considered simulated system, more than 99.5% of the photomultiplier tubes in the detector experience an efficiency loss of less than 1% due to the effect of the magnetic fields.

Many galaxies contain supermassive black holes (SMBHs), whose formation and history raise many puzzles. Pulsar timing arrays have recently discovered a low-frequency cosmological "hum" of gravitational waves that may be emitted by SMBH binary systems, and the JWST and other telescopes have discovered an unexpectedly large population of high-redshift SMBHs. We argue that these two discoveries may be linked, and that they may enhance the prospects for measuring gravitational waves emitted during the mergers of massive black holes, thereby opening the way towards resolving many puzzles about SMBHs as well as providing new opportunities to probe general relativity.

I review recent developments in the field of dark photons-here taken to be U(1) gauge bosons with mass less than the Z-including both kinetically mixed vectors and those that couple to anomaly-free U(1)'s. Distinctions between Higgs and Stueckelberg masses are highlighted, with discussion of swampland constraints, UV completions, and new experimental search strategies.

In this work, a new fast-roll (FR) mechanism to generate primordial black holes (PBHs) and their coeval gravitational waves (GWs) in generalized non-canonical natural inflation is introduced. In this model, choosing a suitable function for non-canonical mass scale parameter $M(\phi)$ gives rise to produce a cliff-like region in field evolution path. When inflaton rolls down the steep cliff, its kinetic energy during the FR stage is amplified in comparison with slow-roll (SR) stage. Hence, seeds of PBHs production are born in this transient FR stage. Depending on the position of the cliff, appropriate cases of PBHs for explaining total dark matter (DM), microlensing effects, LIGO-VIRGO events and NANOGrav 15 year data can be formed. The density spectrum of GWs related to one case of the model lies in the NANOGrav 15 year domain and behaves like $\Omega_{\rm GW_0}\sim f^{5-\gamma}$. The spectral index $\gamma=3.42$ for this case satisfies the NANOGrav 15 year constraint. Moreover, regarding reheating considerations, it is demonstrated that PBHs are born in radiation-dominated (RD) era. Furthermore, viability of the model in light of theoretical swampland criteria and observational constraints on cosmic microwave background (CMB) scales are illustrated.

We explore the fermion oscillation in a degenerate environment. The direct consequence is introducing a Pauli blocking factor $1 - f_i$, where $f_i$ is the phase space distribution function, for each intermediate mass eigenstate during propagation. It is then much easier for a state with larger existing fraction or density to oscillate into other states with less degeneracy while the reversed process is not enhanced. This can significantly modify the oscillation behaviors. We apply this degenerate fermion oscillation to a concrete scenario of neutron-antineutron oscillation in neutron star. It turns out antineutrons receive a standing fraction to annihilate with the environmental neutrons. The subsequent neutron star heating can put an extremely stringent bound on the baryon number violating cross mass term between neutron and antineutron.

With the direct discovery of gravitational waves, black holes have regain interest in the recent years. In particular primordial black holes (PBHs), which originate from the very early Universe, may constitute (at least in part) dark matter. The possibility that dark matter is made of black holes is particularly appealing, and multi-messenger searches are important to probe this hypothesis. In this paper I will discuss the concept of primordial black holes, their origins, their characteristics and the current constraints. In addition I will explain that the study of black holes is of utmost interest since they may constitute portals to new physics and to quantum gravity.