Papers for Thursday, Mar 14 2024

Photometric classifications of supernova (SN) light curves have become necessary to utilize the full potential of large samples of observations obtained from wide-field photometric surveys, such as the Zwicky Transient Facility (ZTF) and the Vera C. Rubin Observatory. Here, we present a photometric classifier for SN light curves that does not rely on redshift information and still maintains comparable accuracy to redshift-dependent classifiers. Our new package, Superphot+, uses a parametric model to extract meaningful features from multiband SN light curves. We train a gradient-boosted machine with fit parameters from 6,061 ZTF SNe that pass data quality cuts and are spectroscopically classified as one of five classes: SN Ia, SN II, SN Ib/c, SN IIn, and SLSN-I. Without redshift information, our classifier yields a class-averaged F1-score of 0.61 +/- 0.02 and a total accuracy of 0.83 +/- 0.01. Including redshift information improves these metrics to 0.71 +/- 0.02 and 0.88 +/- 0.01, respectively. We assign new class probabilities to 3,558 ZTF transients that show SN-like characteristics (based on the ALeRCE Broker light curve and stamp classifiers), but lack spectroscopic classifications. Finally, we compare our predicted SN labels with those generated by the ALeRCE light curve classifier, finding that the two classifiers agree on photometric labels for 82 +/- 2% of light curves with spectroscopic labels and 72% of light curves without spectroscopic labels. Superphot+ is currently classifying ZTF SNe in real time via the ANTARES Broker, and is designed for simple adaptation to six-band Rubin light curves in the future.

Despite being energetically important, the effect of cosmic rays on the dynamics of the interstellar medium (ISM) is assumed to be negligible because the cosmic ray energy diffusion coefficient parallel to the magnetic field is relatively large. Using numerical simulations, we explore how variation of the cosmic ray diffusion coefficient as a function of gas temperature could impact the dynamics of the ISM. We create a two-zone model of cosmic ray transport, reflecting the strong damping of the small scale magnetic field fluctuations, which scatter the cosmic rays, in a gas with low ionization. The variable diffusion coefficient allows more cold gas to form. However, setting the diffusion coefficient at a critical value in the warm phase allows the cosmic rays to adjust the kinetic energy cascade. Specifically, we show the slope of the cascade changes for motion perpendicular to the mean magnetic field, whereas kinetic energy parallel to the magnetic field is reduced equally across inertial scales. We show that cosmic ray energization (or reacceleration) comes at the expense of total radiated energy generated during the formation of a cold cloud. We also show that our two-zone model of cosmic ray transport is capable of matching estimates of the grammage for some paths through the simulation, but full comparison of the grammage requires simulating turbulence in a larger volume.

We describe features of the X-ray: Generate and Analyse (XGA) open-source software package that have been developed to facilitate automated hydrostatic mass ($M_{\rm hydro}$) measurements from XMM X-ray observations of clusters of galaxies. This includes describing how XGA measures global, and radial, X-ray properties of galaxy clusters. We then demonstrate the reliability of XGA by comparing simple X-ray properties, namely the X-ray temperature and gas mass, with published values presented by the XMM Cluster Survey (XCS), the Ultimate XMM eXtragaLactic survey project (XXL), and the Local Cluster Substructure Survey (LoCuSS). XGA measured values for temperature are, on average, within 1% of the values reported in the literature for each sample. XGA gas masses for XXL clusters are shown to be ${\sim}$10% lower than previous measurements (though the difference is only significant at the $\sim$1.8$\sigma$ level), LoCuSS $R_{2500}$ and $R_{500}$ gas mass re-measurements are 3% and 7% lower respectively (representing a 1.5$\sigma$ and 3.5$\sigma$ difference). Like-for-like comparisons of hydrostatic mass are made to LoCuSS results, which show that our measurements are $10{\pm}3%$ ($19{\pm}7%$) higher for $R_{2500}$ ($R_{500}$). The comparison between $R_{500}$ masses shows significant scatter. Finally, we present new $M_{\rm hydro}$ measurements for 104 clusters from the SDSS DR8 redMaPPer XCS sample (SDSSRM-XCS). Our SDSSRM-XCS hydrostatic mass measurements are in good agreement with multiple literature estimates, and represent one of the largest samples of consistently measured hydrostatic masses. We have demonstrated that XGA is a powerful tool for X-ray analysis of clusters; it will render complex-to-measure X-ray properties accessible to non-specialists.

Studying and understanding the physical and chemical processes that govern hot Jupiters gives us insights on the formation of these giant planets. Having a constraint on the molecular composition of their atmosphere can help us pinpoint their evolution timeline. Namely, the metal enrichment and carbon-to-oxygen ratio can give us information about where in the protoplanetary disk a giant planet may have accreted its envelope, and subsequently, indicate if it went through migration. Here we present the first analysis of the atmosphere of the hot Jupiter HIP 65Ab. Using near-infrared high-resolution observations from the IGRINS spectrograph, we detect H$_2$O and CO absorption in the dayside atmosphere of HIP 65Ab. Using a high-resolution retrieval framework, we find a CO abundance of log(CO) = $-3.85^{+0.33}_{-0.36}$, which is slightly under abundant with expectation from solar composition models. We also recover a low water abundance of log(H$_2$O) = $-4.42\pm{0.18}$, depleted by 1 order of magnitude relative to a solar-like composition. Upper limits on the abundance of all other relevant major carbon- and oxygen-bearing molecules are also obtained. Overall, our results are consistent with a sub-stellar metallicity but slightly elevated C/O. Such a composition may indicate that HIP 65Ab accreted its envelope from beyond the water snowline and underwent a disk-free migration to its current location. Alternatively, some of the oxygen on HIP 65Ab could be condensed out of the atmosphere, in which case the observed gas-phase abundances would not reflect the true bulk envelope composition.

White dwarfs with a F, G or K type companion represent the last common ancestor for a plethora of exotic systems throughout the galaxy, though to this point very few of them have been fully characterised in terms of orbital period and component masses, despite the fact several thousand have been identified. Gaia data release 3 has examined many hundreds of thousands of systems, and as such we can use this, in conjunction with our previous UV excess catalogues, to perform spectral energy distribution fitting in order to obtain a sample of 206 binaries likely to contain a white dwarf, complete with orbital periods, and either a direct measurement of the component masses for astrometric systems, or a lower limit on the component masses for spectroscopic systems. Of this sample of 206, four have previously been observed with Hubble Space Telescope spectroscopically in the ultraviolet, which has confirmed the presence of a white dwarf, and we find excellent agreement between the dynamical and spectroscopic masses of the white dwarfs in these systems. We find that white dwarf plus F, G or K binaries can have a wide range of orbital periods, from less than a day to many hundreds of days. A large number of our systems are likely post-stable mass transfer systems based on their mass/period relationships, while others are difficult to explain either via stable mass transfer or standard common envelope evolution.

We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in our universe, and therefore to reside in proto-clusters encompassing an abundance of galaxies in close vicinity. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detected [OIII]-emitting galaxies in the quasars' environments, ranging between a density enhancement of $\delta>100$ within a $2$ cMpc radius - one of the largest proto-clusters during the Epoch of Reionization discovered to date - to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy auto-correlation function, we infer a correlation length of quasars at $\langle z\rangle=6.25$ of $r_0^{\rm QQ}=21.3^{+2.7}_{-2.6}~{\rm cMpc}\,h^{-1}$, while we obtain a correlation length of the [OIII]-emitting galaxies of $r_0^{\rm GG}=4.2\pm0.1~{\rm cMpc}\,h^{-1}$. By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars' host dark matter halos to be $\log_{10}(M_{\rm halo, min}/M_\odot)=12.30\pm0.14$ (and $\log_{10}(M_{\rm halo, min}^{\rm [OIII]}/M_\odot) = 10.72\pm0.03$ for the [OIII]-emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is $f_{\rm duty}\ll1$. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates.

Based on our Chandra imaging-spectroscopic observations, we present the latest evolution of the X-ray remnant of SN 1987A. Recent changes in the electron temperatures and volume emission measures suggest that the blast wave in SN 1987A is moving out of the dense inner ring structure, also called the equatorial ring (ER). The 0.5-2.0 keV X-ray light curve shows a linearly declining trend (by $\sim$4.5 % yr$^{-1}$) between 2016 and 2020, as the blast wave heats the hitherto unknown circumstellar medium (CSM) outside the ER. While the peak X-ray emission in the latest 0.3-8.0 keV image is still within the ER, the radial expansion rate in the 3.0-8.0 keV images suggests an increasing contribution of the X-ray emission from less dense CSM since 2012, at least partly from beyond the ER. It is remarkable that, since 2020, the declining soft X-ray flux has stabilized around $\sim$7 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, which may signal a contribution from the reverse-shocked outer layers of ejecta as predicted by the 3-D magneto-hydrodynamic (MHD) models. In the latest ACIS spectrum of supernova remnant (SNR) 1987A in 2022 we report a significant detection of the Fe K line at $\sim$6.7 keV, which may be due to changing thermal conditions of the X-ray emitting CSM and/or the onset of reverse shock interactions with the Fe-ejecta.

Evidence for the majority of the supermassive black holes in the local universe has been obtained dynamically from stellar motions with the Schwarzschild orbit superposition method. However, there have been only a handful of studies using simulated data to examine the ability of this method to reliably recover known input black hole masses $M_{BH}$ and other galaxy parameters. Here we conduct a comprehensive assessment of the reliability of the triaxial Schwarzschild method at $\textit{simultaneously}$ determining $M_{BH}$, stellar mass-to-light ratio $M^{*}/L$, dark matter mass, and three intrinsic triaxial shape parameters of simulated galaxies. For each of 25 rounds of mock observations using simulated stellar kinematics and the $\texttt{TriOS}$ code, we derive best-fitting parameters and confidence intervals after a full search in the 6D parameter space with our likelihood-based model inference scheme. The two key mass parameters, $M_{BH}$ and $M^{*}/L$, are recovered within the 68% confidence interval, and other parameters are recovered between 68% and 95% confidence intervals. The spatially varying velocity anisotropy of the stellar orbits is also well recovered. We explore whether the goodness-of-fit measure used for galaxy model selection in our pipeline is biased by variable complexity across the 6D parameter space. In our tests, adding a penalty term to the likelihood measure either makes little difference, or worsens the recovery in some cases.

We perform a comprehensive CO study toward the Monoceros OB1 (Mon OB1) region based on the MWISP survey at an angular resolution of about $50''$. The high-sensitivity data, together with the high dynamic range, shows that molecular gas in the $\rm 8^{\circ}\times4^{\circ}$ region displays complicated hierarchical structures and various morphology (e.g., filamentary, cavity-like, shell-like, and other irregular structures). Based on Gaussian decomposition and clustering for $\mathrm{^{13}CO}$ data, a total of 263 $\mathrm{^{13}CO}$ structures are identified in the whole region, and 88\% of raw data flux is recovered. The dense gas with relatively high column density from the integrated CO emission is mainly concentrated in the region where multiple $\rm ^{13}CO$ structures are overlapped. Combining the results of 32 large $\mathrm{^{13}CO}$ structures with distances from Gaia DR3, we estimate an average distance of $\rm 729^{+45}_{-45}~pc$ for the GMC complex. The total mass of the GMC Complex traced by $\mathrm{^{12}CO}$, $\mathrm{^{13}CO}$, and $\mathrm{C^{18}O}$ are $1.1\times10^5~M_\odot$, $4.3\times10^4~M_\odot$, and $8.4\times10^3~M_\odot$, respectively. The dense gas fraction shows a clear difference between Mon OB1 GMC East (12.4\%) and Mon OB1 GMC West (3.3\%). Our results show that the dense gas environment is closely linked to the nearby star-forming regions. On the other hand, star-forming activities have a great influence on the physical properties of the surrounding molecular gas (e.g., greater velocity dispersion, higher temperatures, and more complex velocity structures, etc.). We also discuss the distribution/kinematics of molecular gas associated with nearby star-forming activities.

Hot Jupiters are expected to form far from their host star and move toward close-in, circular orbits via a smooth, monotonic decay due to mild and constant tidal dissipation. Yet, three systems have recently been found exhibiting planet-induced stellar pulsations suggesting unexpectedly strong tidal interactions. Here we combine stellar evolution and tide models to show that dynamical tides raised by eccentric gas giants can give rise to chains of resonance locks with multiple modes, enriching the dynamics seen in single-mode resonance locking of circularized systems. These series of resonance locks yield orders-of-magnitude larger changes in eccentricity and harmonic pulsations relative to those expected from a single episode of resonance locking or nonresonant tidal interactions. Resonances become more frequent as a star evolves off the main sequence providing an alternative explanation to the origin of some stellar pulsators and yielding the concept of "dormant migrating giants". Evolution trajectories are characterized by competing episodes of inward/outward migration and spin-up/-down of the star which are sensitive to the system parameters, revealing a new challenge in modeling migration paths and in contextualizing the observed populations of giant exoplanets and stellar binaries. This sensitivity however offers a new window to constrain the stellar properties of planetary hosts via tidal asteroseismology.

Strongly lensed Type Ia supernovae (LSNe Ia) are a promising probe to measure the Hubble constant ($H_0$) directly. To use LSNe Ia for cosmography, a time-delay measurement between the multiple images, a lens-mass model, and a mass reconstruction along the line of sight are required. In this work, we present the machine learning network LSTM-FCNN which is a combination of a Long Short-Term Memory Network (LSTM) and a fully-connected neural network (FCNN). The LSTM-FCNN is designed to measure time delays on a sample of LSNe Ia spanning a broad range of properties, which we expect to find with the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST) and for which follow-up observations are planned. With follow-up observations in $i$ band (cadence of one to three days with a single-epoch $5\sigma$ depth of 24.5 mag), we reach a bias-free delay measurement with a precision around 0.7 days over a large sample of LSNe Ia. The LSTM-FCNN is far more general than previous machine learning approaches such as the Random Forest (RF), where a RF has to be trained for each observational pattern separately, and yet the LSTM-FCNN outperforms the RF by a factor of roughly three. Therefore, the LSTM-FCNN is a very promising approach to achieve robust time delays in LSNe Ia, which is important for a precise and accurate constraint on $H_0$

In this study, we compare the radio, optical and environmental properties of GRGs with those of a control sample of smaller RGs we found in the three LOw-Frequency ARray (LOFAR) deep fields, namely the Bootes, ELAIS-N1, Lockman Hole, for a total area of about 95 deg^2. We inspected the LOFAR deep fields and created a catalogue of 1609 extended radio galaxies (ERGs). By visual inspection, we identified their host galaxies and spectroscopically or photometrically classified 280 of these as GRGs. We studied their properties, such as their accretion state, stellar mass and star formation rate (SFR) using deep optical and infrared survey data. Moreover, we explored the environment in terms of the surface number density of neighbouring galaxies within these surveys. Integrated flux densities and radio luminosities were also determined for a subset of ERGs through available survey images at 50, 150, 610, and 1400 MHz to compute integrated spectral indices. Considering the fraction of GRGs displaying an FRII morphology alongside the host galaxy properties, we suggest that GRGs consistently possess sufficient power to overcome jet frustration caused by the interstellar medium. Moreover, clear differences emerge in the environmental densities between GRGs and smaller RGs, using the number of neighbouring galaxies within 10 Mpc from the host galaxy as a proxy. GRGs preferentially reside in sparser environments compared to their smaller counterparts. In particular, only 3.6% of the GRGs reside within a 3D comoving distance of 5 Mpc from a previously reported galaxy cluster. We found that larger sources exhibit steeper integrated spectral indices, suggesting that GRGs are late-stage versions of RGs. These results suggest that GRGs are amongst the oldest radio sources with the most stable nuclear activity that reside in sparse environments.

The Sun experienced a period of unprecedented activity during the 20th century, now called the Modern Maximum (MM). The decay of the MM after cycle 19 has changed the Sun, the heliosphere, and the planetary environments in many ways. However, studies disagree on whether this decay has proceeded synchronously in different solar parameters or not. One key issue is if the relation between two long parameters of solar activity, the sunspot number and the solar 10.7cm radio flux, has remained the same during this decay. A recent study argues that there is an inhomogeneity in the 10.7cm radio flux in 1980, which leads to a step-like jump ("1980 jump") in this relation. Here we show that the relation between sunspot number and 10.7cm radio flux varies in time, not due to an inhomogeneous radio flux but due to physical changes in the solar atmosphere. We used radio fluxes at four different wavelengths measured in Japan, and studied their long-term relation with the sunspot number and the 10.7cm radio flux. We also used two other solar parameters, the MgII index and the number of active regions. We find that the 1980 jump is only the first of a series of 1-2-year "humps" that mainly occur during solar maxima. All radio fluxes increase with respect to the sunspot number from the 1970s to 2010s. These results reestablish the 10.7cm flux as a homogeneous measure of solar activity. The fluxes of the longer radio waves are found to increase with respect to the shorter waves, which suggests a long-term change in the solar radio spectrum. We also find that the MgII index and the number of active regions also increased with respect to the sunspot number, further verifying the difference in the long-term evolution in chromospheric and photospheric parameters. Our results provide evidence for important structural changes in solar magnetic fields and the solar atmosphere during the decay of the MM.

Yellow supergiants (YSGs) are rare and poorly understood, and studying them is critical to constraining massive star evolution. We obtained flux-calibrated Magellan Inamori Kyocera Echelle (MIKE) high-resolution spectra of 40 YSGs in the Large Magellanic Cloud (LMC); this sample likely contains post-red supergiants (RSGs). Fitting these data with ATLAS9 model atmospheres, we determined fundamental parameters for these stars. We measure the first spectroscopic luminosities for YSGs above 20 $M_\odot$, providing us a novel probe of the luminosity-to-mass ratio. Many stars in our sample appear to have anomalously high surface gravities, despite being confirmed LMC supergiants. We manually inspected our data finding evidence for binary companions and ongoing mass loss. Our work demonstrates the valuable role of high-resolution spectroscopy in interpreting the evolutionary status of cool supergiants.

Stellar abundance measurements are subject to systematic errors that induce extra scatter and artificial correlations in elemental abundance patterns. We derive empirical calibration offsets to remove systematic trends with surface gravity $\log(g)$ in 17 elemental abundances of 288,789 evolved stars from the SDSS APOGEE survey. We fit these corrected abundances as the sum of a prompt process tracing core-collapse supernovae and a delayed process tracing Type Ia supernovae, thus recasting each star's measurements into the amplitudes $A_{\text{cc}}$ and $A_{\text{Ia}}$ and the element-by-element residuals from this two-parameter fit. As a first application of this catalog, which is $8\times$ larger than that of previous analyses that used a restricted $\log(g)$ range, we examine the median residual abundances of 14 open clusters, nine globular clusters, and four dwarf satellite galaxies. Relative to field Milky Way disk stars, the open clusters younger than 2 Gyr show $\approx 0.1-0.2$ dex enhancements of the neutron-capture element Ce, and the two clusters younger than 0.5 Gyr also show elevated levels of C+N, Na, S, and Cu. Globular clusters show elevated median abundances of C+N, Na, Al, and Ce, and correlated abundance residuals that follow previously known trends. The four dwarf satellites show similar residual abundance patterns despite their different star formation histories, with $\approx 0.2-0.3$ dex depletions in C+N, Na, and Al and $\approx 0.1$ dex depletions in Ni, V, Mn, and Co. We provide our catalog of corrected APOGEE abundances, two-process amplitudes, and residual abundances, which will be valuable for future studies of abundance patterns in different stellar populations and of additional enrichment processes that affect galactic chemical evolution.

(Abridged) Recent molecular surveys have revealed a rich gas organization of sonic-like fibers in all kind of environments prior to the formation of low- and high-mass stars. This paper introduces the EMERGE project aiming to investigate whether complex fiber arrangements could explain the origin of high-mass stars and clusters. We analyzed the EMERGE Early ALMA Survey including 7 star-forming regions in Orion (OMC-1/2/3/4 South, L1641N, NGC2023, and Flame Nebula) homogeneously surveyed in both molecular lines (N$_2$H$^+$ J=1-0, HNC J=1-0, plus HC3N J=10-9) and 3mm-continuum using a combination of interferometric ALMA mosaics and IRAM-30m single-dish (SD) maps. Based on our low-resolution (SD) observations, we describe the global properties of our sample covering a wide range of physical conditions including low-, intermediate, and high-mass star-forming regions in different evolutionary stages. Their comparison with ancillary YSO catalogs denotes N$_2$H$^+$ as the best proxy for the dense, star-forming gas in our targets showing a constant star formation efficiency and a fast time evolution of <1 Myr. While apparently clumpy and filamentary in our SD data, all targets show a much more complex fibrous substructure at the enhanced resolution of our ALMA+IRAM-30m maps. A large number of filamentary features at sub-parsec scales are clearly recognized in the high-density gas traced by N$_2$H$^+$ directly connected to the formation of individual protostars. This complex gas organization appears to extend further into the more diffuse gas traced by HNC. This paper presents the EMERGE Early ALMA survey including a first data release of continuum maps and spectral products for this project to be analysed in future papers of this series. A first look at these results illustrates the need of advanced data combination techniques to investigate the intrinsic multi-scale, gas structure of the ISM.

The search for dual supermassive black holes (SMBHs) is of immense interest in modern astrophysics. Galaxy mergers may be an important route to fuel and to produce SMBH pairs. Actively accreting SMBH pairs can be observed as a dual quasar, which are vital probes of SMBH growth. Gaia observations have enabled a novel technique to systematically search for such dual quasars at previously unreachable sub-kpc scales, based on the small jitters of the light centroid as the two quasars vary stochastically. Here we present the first detailed study of a 0.46'', 3.8 kpc separation, VODKA-selected dual quasar, J0749+2255, at $z=2.17$ using JWST/NIRSpec integral field unit spectroscopy. This is one of the most distant, small separation dual quasars identified today. Dual quasars at cosmic noon are not well characterized. We detect the faint ionized gas of the host galaxy, best traced by the narrow \ha\ emission. Line ratio diagnostics show a mix of ionization from the two quasars and intense star formation. The spatially-resolved spectra of the two quasars suggest that they have very similar black hole properties (two $M_{BH}\sim 10^9\ \textrm{M}_{\odot}$ with large Eddington ratio reaching $L/L_{Edd}\sim0.2$) hinting at the possible synchronized growth and accretion from the same gas supply. Surprisingly, the ionized gas kinematics suggest an extended, rotating disk rather than a disturbed system that would be expected in a major gas-rich galaxy merger. While it is unclear if J0749+2255 is representative of the dual quasar evolution, the observations with JWST revealed a major puzzle. It would be interesting to see what observations of other dual quasars will show.

This paper presents the results of a systematic study of projection biases in the Weak Lensing analysis of cosmic shear and the combination of galaxy clustering and galaxy-galaxy lensing using data collected during the first-year of running the Dark Energy Survey experiment. The study uses $\Lambda$CDM as the cosmological model and two-point correlation functions for the WL analysis. The results in this paper show that, independent of the WL analysis, projection biases of more than $1\sigma$ exist, and are a function of the position of the true values of the parameters $h$, $n_{s}$, $\Omega_{b}$, and $\Omega_{\nu}h^{2}$ with respect to their prior probabilities. For cosmic shear, and the combination of galaxy clustering and galaxy-galaxy lensing, this study shows that the coverage probability of the $68.27\%$ credible intervals ranges from as high as $93\%$ to as low as $16\%$, and that these credible intervals are inflated, on average, by $29\%$ for cosmic shear and $20\%$ for the combination of galaxy clustering and galaxy-galaxy lensing. The results of the study also show that, in six out of nine tested cases, the reduction in error bars obtained by transforming credible intervals into confidence intervals is equivalent to an increase in the amount of data by a factor of three.

The 2003 October 28 (X17.2) eruptive flare was a unique event. The coronal electric field and the {\pi}-decay {\gamma}-ray emission flux had the highest values ever inferred in solar flares. This study reveals physical links between the magnetic reconnection process, the energy release, and the acceleration of electrons and ions to high energies in the chain of the magnetic energy transformations in the impulsive phase of the solar flare. The global reconnection rate and the local reconnection rate are calculated from flare ribbon separation in H{\alpha} filtergrams and photospheric magnetic field maps. Available results of INTEGRAL and CORONAS-F/SONG observations are combined with Konus-Wind data to quantify time behavior of electron and proton acceleration. Prompt {\gamma}-ray lines and delayed 2.2 MeV line temporal profiles observed with Konus-Wind and INTEGRAL/SPI used to detect and quantify the nuclei with energies of 10-70 MeV. The global and local reconnection rates reach their peaks at the end of the main rise phase of the flare. The spectral analysis of the high-energy {\gamma}-ray emission revealed a close association between the acceleration process efficiency and the reconnection rates. High-energy bremsstrahlung continuum and narrow {\gamma}-ray lines were observed in the main rise phase. In the main energy release phase, the upper energy of the bremsstrahlung spectrum was significantly reduced and the pion-decay {\gamma}-ray emission appeared abruptly. We discuss the reasons why the change of the acceleration regime occurred along with the large-scale magnetic field restructuration of this flare. We argue that the main energy release and proton acceleration up to subrelativistic energies began just when the reconnection rate was going through the maximum, i.e., after a major change of the flare topology.

Diffusion coefficients of energetic charged particles in turbulent magnetic fields are a fundamental aspect of diffusive transport theory but remain incompletely understood. In this work, we use quasi-linear theory to evaluate the spatial variation of the parallel diffusion coefficient $\kappa_\parallel$ from the measured magnetic turbulence power spectra in the inner heliosphere. We consider the magnetic field and plasma velocity measurements from Parker Solar Probe made during Orbits 5-13. The parallel diffusion coefficient is calculated as a function of radial distance from 0.062 to 0.8 AU, and the particle energy from 100 keV to 1GeV. We find that $\kappa_\parallel$ increases exponentially with both heliocentric distance and energy of particles. The fluctuations in $\kappa_\parallel$ are related to the episodes of large-scale magnetic structures in the solar wind. By fitting the results, we also provide an empirical formula of $\kappa_{\parallel}=(5.16\pm1.22) \times 10^{18} \: r^{1.17\pm0.08} \: E^{0.71\pm 0.02} \; (cm^2/s)$ in the inner heliosphere which can be used as a reference in studying the transport and acceleration of solar energetic particles as well as the modulation of cosmic rays.

Recently discovered inertial waves, observed on the solar surface, likely extend to the deeper layers of the Sun. Utilizing helioseismic techniques, we explore these motions, allowing us to discern inertial-mode eigenfunctions in both radial and latitudinal orientations. We analyze $8$ years of space-based observations ($2010 - 2017$) taken by the Helioseismic and Magnetic Imager (HMI) onboard the Solar dynamic observatory (SDO) using normal-mode coupling. Coupling between same and different-degree acoustic modes and different frequency bins are measured in order to capture the various length scales of inertial modes. We detect inertial modes at high latitude with azimuthal order $t=1$ and frequency $\sim -80$ nHz. This mode is present in the entire convection zone. The presence of Rossby modes may be seen down to a depth of $\sim 0.83R_\odot$ and the Rossby signal is indistinguishable from noise below that depth for high azimuthal order. We find that the amplitudes of these modes increase with depth down to around $0.92 R_\odot$ and decrease below that depth. We find that the latitudinal eigenfunctions of Rossby modes deviate from sectoral spherical harmonics if we use a similar approach as adopted in earlier studies. We found that spatial leakage and even pure noise in the measurements of non-sectoral components can also explain the above-mentioned characteristics of the latitudinal eigenfunctions. This realization underscores the necessity for careful interpretation when considering the latitudinal eigenfunctions of Rossby modes. Exploring the depth-dependent characteristics of these modes will enable us to capture interior dynamics distinctly, separate from p-mode seismology.

We present SN 2023zaw $-$ a sub-luminous ($\mathrm{M_r} = -16.7$ mag) and rapidly-evolving supernova ($\mathrm{t_{1/2,r}} = 4.9$ days), with the lowest nickel mass ($\approx0.002$ $\mathrm{M_\odot}$) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He I and Ca NIR emission lines with velocities of $\sim10\ 000 - 12\ 000$ $\mathrm{km\ s^{-1}}$. The late-time spectra show prominent narrow He I emission lines at $\sim$1000$\ \mathrm{km\ s^{-1}}$, indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of $\approx0.2$ $\mathrm{M_\odot}$ and an envelope radius of $\approx50$ $\mathrm{R_\odot}$. The extremely low nickel mass and low ejecta mass ($\approx0.5$ $\mathrm{M_\odot}$) suggest an ultra-stripped SN, which originates from a mass-losing low mass He-star (ZAMS mass $<$ 10 $\mathrm{M_\odot}$) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass ($< 0.005$ $\mathrm{M_\odot}$) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys.

Historically, metallicity profiles of galaxies have been modelled using a radially symmetric, two-parameter linear model, which reveals that most galaxies are more metal-rich in their central regions than their outskirts. However, this model is known to yield inaccurate results when the point-spread function (PSF) of a telescope is large. Furthermore, a radially symmetric model cannot capture asymmetric structures within a galaxy. In this work, we present an extension of the popular forward-modelling python package LENSTRONOMY, which allows the user to overcome both of these obstacles. We demonstrate the new features of this code base through two illustrative examples on simulated data. First, we show that through forward modelling, LENSTRONOMY is able to recover accurately the metallicity gradients of galaxies, even when the PSF is comparable to the size of a galaxy, as long as the data is observed with a sufficient number of pixels. Additionally, we demonstrate how LENSTRONOMY is able to fit irregular metallicity profiles to galaxies that are not well-described by a simple surface brightness profile. This opens up pathways for detailed investigations into the connections between morphology and chemical structure for galaxies at cosmological distances using the transformative capabilities of JWST. Our code is publicly available and open source, and can also be used to model spatial distributions of other galaxy properties that are traced by its surface brightness profile.

As the most ancient branch of astronomy, astrometry has been developed for thousands of years. However, it has only recently become possible to utilize astrometry for the detection of exoplanets. Gaia, an astrometric surveyor of 1 billion stars, is capable of measuring the position of stars with a precision as high as 20 $\mu$as. Gaia is expected to discover more than 10,000 exoplanets by the end of its mission, surpassing the productivity of most exoplanet surveys. In this chapter, I will introduce different techniques used to achieve high-precision astrometry. Subsequently, I will explore how both relative and absolute astrometry can be employed to detect exoplanets. Finally, I will present the detection limit of the Gaia astrometric survey.

We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-$\alpha$ (Ly-$\alpha$) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of 205,012 quasar spectra in the redshift range 2 <= z <= 4 at an effective redshift z=2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp(ikr), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artifacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line-of-sight of Ly-$\alpha$ skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly-$\alpha$ skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range 0.02 <= k <= 0.35 h/Mpc at the 1-2 sigma level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly-$\alpha$ forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO and 4MOST.

We measure the evolution of the rest-frame $NUV-V$ colors for early-type galaxies in clusters at $0<z<1.1$ using data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), CFHT Large Area U-band Deep Survey (CLAUDS) and local SDSS clusters observed with GALEX. Our results show that there is an excess in the ultraviolet spectrum in most quiescent galaxies (compared to the expectations from models fitting their optical/infrared colors and spectra) below $z\sim0.6$, beyond which the excess UV emission fades rapidly. This evolution of the UV color is only consistent with the presence of a highly evolved, hot horizontal branch sub-population in these galaxies (amongst the majority cool and optically bright stars), comprising on average 10\% of the total stellar mass and forming at $z>3$. The blue UV colors of early-type galaxies at low-intermediate redshifts are likely driven by this sub-population being enriched in helium up to $\sim44\%$. At $z>0.8$ (when the extra UV component has not yet appeared) the data allows us to constrain the star formation histories of galaxies by fitting models to the evolution of their UV colors: we find that the epoch at which the stellar populations formed ranges between $3<z_{form}<10$ (corresponding to $0.5-2.2$ Gyrs after the Big Bang) with a star-formation e-folding timescale of $\tau=0.35-0.7$ Gyr, suggesting that these galaxies formed the majority of stars at very high redshift, with a brief yet intense burst of star-formation activity. The star formation history and chemical evolution of early-type galaxies resemble those of globular clusters, albeit on much larger scales.

We have translated the results of $N$-body simulations of one barred model into the language of action variables and frequencies. Using this language, we analysed the behaviour of all orbits in the model on a large time scale at the stage of a mature bar. We show that the orbits join the bar while preserving their adiabatic invariant, which takes into account the 3D structure of the orbits. This allows us to apply the concept of the Lynden-Bell derivative for each of these orbits and trace how the sign of the derivative changes, i.e. how asynchronous changes in angular momentum $L_z$ and orbital precession rate $\Omega_\mathrm{pr}$ (normal orbital mode) change to synchronous (abnormal mode). The transition to the abnormal mode occurs when $\Omega_\mathrm{pr}$ reaches the angular velocity of the pattern $\Omega_\mathrm{p}$, after which the orbit becomes stuck in the bar trap. All this happens against the background of secular changes in actions ($L_z$ decreases, $J_\mathrm{R}$ and $J_z$ increase). At the same time, corotation particles near two stable Lagrange points are also subject to secular changes in their actions. They increase $L_z$ and drift to the periphery, shifting corotation outwards. We also show that a change in the orbital mode from normal to abnormal and the trapping of orbits in a bar is possible only when the bar speed decreases with time, regardless of what is causing the bar to slow down. Our findings clarify and expand the picture of bar formation and evolution in numerical models.

We study the path of light rays passing near a massive object, in the context of the scale invariant equation of the geodesics first obtained by Dirac (1973). Using the exterior Schwarzschild solution for the metric, we derive the complete equations of the geodesics in the scale invariant context. We find that scale invariance introduces two additional terms to the Einstein term producing the deflection angle and that can potentially act over cosmological distances. Numerical integration of the scale-invariant geodesics, for the specific case of the z_L=1.94 lens galaxy in the extreme system JWST-ER1 (van Dokkum et al. 2023) shows that the two additional terms introduce only negligible effects, typically 1E-06 of the Einstein term. We conclude that the lensing deflection angle derived in Einstein's General Relativity is essentially independent of the scale invariant effects and that the photon's geodesics remain unchanged. We also explore the possible origin of the differences in the mass estimates from lensing and photometry in JWST-ER1 and in the SLACS galaxies, differences which appear larger at higher redshifts.

Observations performed with high-resolution imaging techniques revealed the existence of shadows in circumstellar disks that can be explained by the misalignment of an inner with respect to an outer disk. The cause of misalignment, however, is still debated. In this study, we investigate the feasibility of observing shadows induced by one prominent scenario that may lead to misalignment, which involves the late infall of material onto a protostellar system. In particular, we use previously performed hydrodynamical simulations of such events, and generate flux maps in the visible, near-infrared, submillimeter, and millimeter wavelength range using Monte Carlo radiative transfer. Based on that, we derive synthetic observations of these systems performed with the instruments SPHERE/VLT and ALMA, which we use as a basis for our subsequent analysis. We find that near-infrared observations with SPHERE are particularly well suited for detecting shadows via direct imaging alongside other features such as gaps, arcs, and streamers. On the contrary, performing a shadow detection based on reconstructed ALMA observations is very challenging due to the high sensitivity that is required for this task. Thus, in cases that allow for a detection, sophisticated analyses may be needed, for instance by the utilization of carefully constructed azimuthal profiles, aiding the search for potentially shallow shadows. Lastly, we conclude that late infall-induced disk misalignment offers a plausible explanation for the emergence of shadows that are observed in various systems.

To study the chemical evolution across cosmic epochs, we investigate Ne, S, Cl, and Ar abundance patterns in the COS Legacy Archive Spectroscopic SurveY (CLASSY). CLASSY comprises local star-forming galaxies (0.02 < z < 0.18) with enhanced star-formation rates, making them strong analogues to high-z star-forming galaxies. With direct measurements of electron temperature, we derive accurate ionic abundances for all elements and assess ionization correction factors (ICFs) to account for unseen ions and derive total abundances. We find Ne/O, S/O, Cl/O, and Ar/O exhibit constant trends with gas-phase metallicity for 12+log(O/H) < 8.5 but significant correlation for Ne/O and Ar/O with metallicity for 12+log(O/H) > 8.5, likely due to ICFs. Thus, applicability of the ICFs to integrated spectra of galaxies could bias results, underestimating true abundance ratios. Using CLASSY as a local reference, we assess the evolution of Ne/O, S/O, and Ar/O in galaxies at z>3, finding no cosmic evolution of Ne/O, while the lack of direct abundance determinations for S/O and Ar/O can bias the interpretation of the evolution of these elements. We determine the fundamental metallicity relationship (FMR) for CLASSY and compare to the high-redshift FMR, finding no evolution. Finally, we perform the first mass-neon relationship analysis across cosmic epochs, finding a slight evolution to high Ne at later epochs. The robust abundance patterns of CLASSY galaxies and their broad range of physical properties provide essential benchmarks for interpreting the chemical enrichment of the early galaxies observed with the JWST.

Modeling the surface brightness distribution of stars is of prime importance to interpret observations. Nevertheless, this remains quite challenging for cool stars as it requires one to model the MHD turbulence that develops in their convective envelope. In Paper I, the effect of the Coriolis acceleration on the surface heat flux has been studied by means of hydrodynamic simulations. In this paper, we aim to investigate the additional effect of dynamo magnetic fields. We focus on an envelope thickness that is representative of either a $\sim0.35~M_\odot$ M dwarf, a young red giant star or a pre-main sequence star. We performed a parametric study using numerical MHD simulations of anelastic convection in thick rotating spherical shells. For each model, we computed the mean surface distribution of the heat flux, and examined the leading-order effect of the magnetic field on the obtained latitudinal luminosity profile. We identify three different regimes. Close to the onset of convection, while the first unstable modes tend to convey heat more efficiently near the equator, magnetic fields are shown to generally enhance the mean heat flux close to the polar regions (and the tangent cylinder). By progressively increasing the Rayleigh number, the development of a prograde equatorial jet was previously shown to make the equator darker when no magnetic field is taken into account. For moderate Rayleigh numbers, magnetic fields can instead inverse the mean pole-equator brightness contrast (which means going from a darker to a brighter equator when a dynamo sets in) and finally induce a similar regime to that found close to the onset of convection. For more turbulent models with larger Rayleigh numbers, magnetic fields alternatively tend to smooth out the brightness contrast. This general behavior is shown to be related to the quenching of the surface differential rotation by magnetic fields.

The origin of the diffuse gamma-ray background in the range from hundreds keV to several MeV is not known conclusively. From current models and observations it is believed that, at least partially, this background is formed by blazars and remnants of supernovae (SN) of type Ia in distant galaxies. However, these contributions are not sufficient to reproduce the observed level of the signal. In this work we propose another source which could contribute to this background, namely the jets of active galactic nuclei (AGN). The composition of jets is not known, but there are observational hints that the fraction of positrons there is substantial. Positrons are partially evacuated to the intergalactic medium and partially mix with the circumgalactic medium and annihilate there comparatively quickly. Using the AGN luminosity function, we estimated the positron production rate and the contribution of the positron annihilation to the cosmic background below 511 keV. We also estimated the analogous contribution from positron annihilation within SN Ia remnants in distant galaxies. The contribution of AGNs is estimated to be a factor of 5 - 10 smaller than the observed background intensity, and the contribution from SNe is yet smaller by one order of magnitude. Nevertheless, the contribution of AGNs appeared to be larger than the contribution of blazars estimated from Swift-BAT and Fermi-LAT observations. The main uncertainty in our model is the fraction of positrons remaining in the circumgalactic medium which makes our estimation an upper limit.

We present the first statistical investigation of spatially resolved emission-line properties in a sample of 63 low-mass galaxies at $4\leq z<10$, using JWST/NIRSpec MSA data from the JWST Advanced Deep Extragalactic (JADES) survey focusing on deep, spatially resolved spectroscopy in the GOODS-S extragalactic field. By performing a stacking of the 2D spectra of the galaxies in our sample, we find an increasing or flat radial trend with increasing radius for [OIII]$\lambda5007$/H$\beta$ and a decreasing one for [NeIII]$\lambda3869$/[OII]$\lambda3727$ (3--4 $\sigma$ significance). These results are still valid when stacking the sample in two redshift bins (i.e., $4\leq z<5.5$ and $5.5\leq z<10$). The comparison with star-formation photoionization models suggests that the ionization parameter increases by $\sim 0.5$ dex with redshift. We find a tentative metallicity gradient that increases with radius (i.e., 'inverted') in both redshift bins. Moreover, our analysis reveals strong negative gradients for the equivalent width of \Hbeta (7$\sigma$ significance). This trend persists even after removing known AGN candidates, therefore, it is consistent with a radial gradient primarily in stellar age and secondarily in metallicity. Taken all together, our results suggest that the sample is dominated by active central star formation, with possibly inverted metallicity gradients sustained by recent episodes of accretion of pristine gas or strong radial flows. Deeper observations and larger samples are needed to confirm these preliminary results and to validate our interpretation.

Ultraviolet (UV) processing in the insterstellar medium (ISM) induces the dehydrogenation of hydrocarbons. Aliphatics, including alkanes, are present in different interstellar environments, being prevalently formed in evolved stars; thus, the dehydrogenation by UV photoprocessing of alkanes plays an important role in the chemistry of the ISM, leading to the formation of unsaturated hydrocarbons and eventually to aromatics, the latter ubiquitously detected in the ISM. Here, through combined experimental results and \textit{ab-initio} calculations, we show that UV absorption (mainly at the Ly-$\alpha$ emission line of hydrogen at 121.6 nm) promotes an alkane to an excited Rydberg state from where it evolves towards fragmentation inducing the formation of olefinic C=C bonds, which are necessary precursors of aromatic hydrocarbons. We show that photochemistry of aliphatics in the ISM does not primarily produce direct hydrogen elimination but preferential C-C photocleavage. Our results provide an efficient synthetic route for the formation of unsaturated aliphatics, including propene and dienes, and suggest that aromatics could be formed in dark clouds by a bottom-up mechanism involving molecular fragments produced by UV photoprocessing of aliphatics.

Context. The third Gaia Data Release, which includes BP/RP spectra for 219 million sources, has opened a new window in the exploration of the chemical history and evolution of the Milky Way. The wealth of information encapsulated in these data is far greater than their low resolving power (R=50) at first glance would suggest, as shown in many studies. We zero in on the use of this data for the purpose of the detection of ''new'' metal-poor stars, which are hard to find yet essential for understanding - among other - several aspects of the origin of the Galaxy, star formation and the creation of the elements. Aims. We strive to refine a metal-poor candidate selection method which was developed with simulated Gaia BP/RP spectra, with an ultimate objective of providing the community with both a recipe to select stars for medium/high resolution observations and a catalogue of stellar metallicities. Methods. We used a datased comprised of GALAH DR3 and SAGA database stars in order to verify and adjust to real world data our selection method. For that purpose, we used dereddening as a mean to tackle the issue of extinction, and then we applied our fine-tuned method to select metal-poor candidates, which we thereafter observed and analysed. Results. We were able to infer metallicities for GALAH DR3 and SAGA stars - with color excesses up to E(B-V)<1.5 - with an uncertainty of 0.36 dex, which is good enough for the purpose of identifying new metal-poor stars. Further, we selected 26 metal-poor candidates - via our method - for observations. As spectral analysis showed, 100% of them had [Fe/H]<-2.0, 57% had [Fe/H]<-2.5 and 8% had [Fe/H]<-3.0. We inferred metallicities for these stars with an uncertainty of 0.31 dex, as was proven when comparing to the spectroscopic [Fe/H]. Finally, we assembled a catalogue of metallicities for 10 861 062 stars.

We offer a new way to look at the Large Magellanic Cloud through stellar mono-abundance and mono-age-mono-abundance maps. These maps are based on $\gtrsim 500\,000$ member stars with photo-spectroscopic [M/H] and age estimates from Gaia DR3 data, and they are the first area-complete, metallicity- and age-differentiated stellar maps of any disk galaxy. Azimuthally averaged, these maps reveal a surprisingly simple picture of the Milky Way's largest satellite galaxy. For any [M/H] below -0.1 dex, the LMC's radial profile is well described by a simple exponential, but with a scale length that steadily shrinks towards higher metallicities, from nearly 2.3~kpc at [M/H]$=-1.8$ to only 0.75~kpc at [M/H]$=-0.25$. The prominence of the bar decreases dramatically with [M/H], making it barely discernible at [M/H]$\lesssim -1.5$. Yet, even for metal-rich populations, the bar has little impact on the azimuthally averaged profile of the mono-abundance components. Including ages, we find that the scale length is a greater function of age than of metallicity, with younger populations far more centrally concentrated. At old ages, the scale length decreases with increasing metallicity; at young ages, the scale-length is independent of metallicity. These findings provide quantitative support for a scenario where the LMC built its stellar structure effectively outside in.

I point out similarities between point-symmetric X-ray morphologies in cooling flow groups and clusters of galaxies, which are observed to be shaped by jets, and point-symmetric morphologies of eight core-collapse supernova (CCSN) remnants. I use these similarities to strengthen the jittering jet explosion mechanism (JJEM) of CCSNe, which predicts that the last pairs of jets to be launched by the newly born neutron star might shape some CCSN remnants to point-symmetric morphology. The point-symmetric morphologies in both types of objects are composed of two or more pairs of opposite bubbles (cavities), nozzles, some clumps, small protrusions (termed ears), and rims. The typically large volume of a CCSN remnant shaped by jets implies that the shaping jets carry an energy comparable to that of the ejecta, which in turn implies that jets exploded the remnant's massive star progenitor. The morphological similarities studied here add to the similarity of CCSN remnants, not only point-symmetric ones, to planetary nebulae shaped by jets. Together, these similarities solidify the JJEM as the main explosion mechanism of CCSNe. I consider the identification of point-symmetry in CCSNe, as expected by jet-shaping in the JJEM, to be the most severe challenge to the competing neutrino-driven explosion mechanism. I reiterate my earlier claim, but in a more vocal voice, that the main explosion mechanism of CCSNe is the JJEM.

Recent models of solar system formation suggest that a dynamical instability among the giant planets happened within the first 100 Myr after disk dispersal, perhaps before the Moon-forming impact. As a direct consequence, a bombardment of volatile-rich impactors may have taken place on Earth before internal and atmospheric reservoirs were decoupled. However, such a timing has been interpreted to potentially be at odds with the disparate inventories of Xe isotopes in Earth's mantle compared to its atmosphere. This study aims to assess the dynamical effects of an Early Instability on the delivery of carbonaceous asteroids and comets to Earth, and address the implications for the Earth's volatile budget. We perform 20 high-resolution dynamical simulations of solar system formation from the time of gas disk dispersal, each starting with 1600 carbonaceous asteroids and 10000 comets, taking into account the dynamical perturbations from an early giant planet instability. Before the Moon-forming impact, the cumulative collision rate of comets with Earth is about 4 orders of magnitude lower than that of carbonaceous asteroids. After the Moon-forming impact, this ratio either decreases or increases, often by orders of magnitude, depending on the dynamics of individual simulations. An increase in the relative contribution of comets happens in 30\% of our simulations. In these cases, the delivery of noble gases from each source is comparable, given that the abundance of 132Xe is 3 orders of magnitude greater in comets than in carbonaceous chondrites. The increase in cometary flux relative to carbonaceous asteroids at late times may thus offer an explanation for the Xe signature dichotomy between the Earth's mantle and atmosphere.

Solar radio spikes are short lived, narrow bandwidth features in low frequency solar radio observations. The timing of their occurrence and the number of spikes in a given observation is often unpredictable. The high temporal and frequency of resolution of modern radio telescopes such as NenuFAR mean that manually identifying radio spikes is an arduous task. Machine learning approaches to data exploration in solar radio data is on the rise. Here we describe a convolutional neural network to identify the per pixel location of radio spikes as well as determine some simple characteristics of duration, spectral width and drift rate. The model, which we call SpikeNet, was trained using an Nvidia Tesla T4 16GB GPU with ~100000 sample spikes in a total time of 2.2 hours. The segmentation performs well with an intersection over union in the test set of ~0.85. The root mean squared error for predicted spike properties is of the order of 23%. Applying the algorithm to unlabelled data successfully generates segmentation masks although the accuracy of the predicted properties is less reliable, particularly when more than one spike is present in the same 64 X 64 pixel time-frequency range. We have successfully demonstrated that our convolutional neural network can locate and characterise solar radio spikes in a number of seconds compared to the weeks it would take for manual identification.

Studies of micrometeorites in mid-Ordovician limestones and Earth's impact craters indicate that our planet witnessed a massive infall of ordinary L chondrite material 466 million years (My) ago (Heck et al. 2017, Schmieder & Kring 2020, Kenkmann 2021) that may have been at the origin of the first major mass extinction event (Schmitz et al. 2019). The breakup of a large asteroid in the main belt is the likely cause of this massive infall. In modern times, material originating from this breakup still dominates meteorite falls (>20% of all falls) (Swindle et al. 2014). Here, we provide spectroscopic observations and dynamical evidence that the Massalia collisional family is the only plausible source of this catastrophic event and of the most abundant class of meteorites falling on Earth today. It is suitably located in the inner belt, at low-inclination orbits, which corresponds to the observed distribution of L-chondrite-like near-Earth objects (NEOs) and of interplanetary dust concentrated at 1.4 degrees (Sykes 1990, Reach et al. 1997).

Understanding the origin of bright shooting stars and their meteorite samples is among the most ancient astronomy-related questions that at larger scales has human consequences [1-3]. As of today, only ${\sim}\,6\%$ of meteorite falls have been firmly linked to their sources (Moon, Mars, and asteroid (4) Vesta [4-6]). Here, we show that ${\sim}\,70\%$ of meteorites originate from three recent breakups of $D > 30\,{\rm km}$ asteroids that occurred 5.8, 7.5 and less than ${\sim}\,40$ million years ago. These breakups, including the well-known Karin family [7], took place in the prominent yet old Koronis and Massalia families and are at the origin of the dominance of H and L ordinary chondrites among meteorite falls. These young families distinguish themselves amidst all main belt asteroids by having a uniquely high abundance of small fragments. Their size-frequency distribution remains steep for a few tens of millions of years, exceeding temporarily the production of metre-sized fragments by the largest old asteroid families (e.g., Flora, Vesta). Supporting evidence includes the existence of associated dust bands [8-10], the cosmic-ray exposure ages of H-chondrite meteorites [11,12], or the distribution of pre-atmospheric orbits of meteorites [13-15].

The radiation physics of fast radio bursts (FRBs) remains an open question. Current observations have discovered that narrowly-banded bursts of FRB 20201124A are active in 0.4-2 GHz and their spectral peak frequency ($\nu^{\rm obs}_{p}$) are mostly toward $\sim 1$ GHz. Utilizing a sample of 1268 bursts of FRB 20201124A detected with the FAST telescope, we show that the $1\sigma$ spectral regime of 71.4\% events (in-band bursts) is within the FAST bandpass. Their intrinsic burst energies ($E^{\rm obs}_{\rm BWe}$) and spectral widths ($\sigma_s^{\rm obs}$) are well measured by fitting the spectral profile with a Gaussian function. The derived $E^{\rm obs}_{\rm BWe}$ and $\sigma_s^{\rm obs}$ distributions are log-normal and centering at $\log E^{\rm obs}_{\rm BWe}/{\rm erg}=37.2~ (\sigma=0.76)$ and $\log \sigma_s^{\rm obs}/{\rm GHz}=-1.16~ (\sigma=0.17)$. Our Monte Carlo simulation analysis infers its intrinsic $\nu_p$ distribution as a normal function centered at $\nu_{p,c}=1.16$ GHz ($\sigma=0.22$) and its intrinsic energy function as $\Phi(E)\propto E^{-0.60}e^{-E/E_c}$ with $E_c=9.49 \times 10^{37}$ erg. We compare these results with that of typical repeating FRBs 20121102A and 20190520B that are active over a broad frequency range at several specific frequencies and discuss possible observational biases on the estimation of the event rate and energy function. Based on these results, we argue that FRB 20201124A likely occurs in a fine-tuned plasma for maser radiations at a narrow frequency range, while FRB 20121102A and FRB 20190520B could involve clumpy plasma conditions that make maser emission around several specific frequencies in a broad range.

Dynamo models of stellar magnetic fields for partly and fully convective stars are guided by observational constraints. Zeeman-Doppler imaging has revealed a variety of magnetic field geometries and, for fully convective stars in particular, a dichotomy: either strong, mostly axisymmetric, and dipole-dominated or weak, non-axisymmetric, and multipole-dominated. This dichotomy is explained by dynamo bistability or by long-term magnetic cycles, but there is no definite conclusion on the matter. We analysed optical spectropolarimetric data sets collected with ESPaDOnS and Narval between 2005 and 2016, and near-infrared SPIRou data obtained between 2019 and 2022 for three active M dwarfs with masses between 0.1 and 0.6 MSun: EV Lac, DS Leo, and CN Leo. We looked for changes in time series of longitudinal magnetic field, width of unpolarised mean-line profiles, and large-scale field topology as retrieved with principal component analysis and Zeeman-Doppler imaging. We retrieved pulsating (EV Lac), stable (DS Leo), and sine-like (CN Leo) long-term trends in longitudinal field. The width of near-infrared mean-line profiles exhibits rotational modulation only for DS Leo, whereas in the optical it is evident for both EV Lac and DS Leo. The line width variations are not necessarily correlated to those of the longitudinal field, suggesting complex relations between small- and large-scale field. We also recorded topological changes: a reduced axisymmetry for EV Lac and a transition from toroidal- to poloidal-dominated regime for DS Leo. For CN Leo, the topology remained dipolar and axisymmetric, with only an oscillation in field strength. Our results show a peculiar evolution of the magnetic field for each M dwarf, confirming that M dwarfs with distinct masses and rotation periods can undergo magnetic long-term variations, and suggesting a variety of cyclic behaviours of their magnetic fields.

Primordial tensor modes can induce Cosmic Microwave Background spectral distortions during horizon re-entry. We investigate a specific mechanism proposed for this purpose, characterized by the coupling of an SU(2) gauge field to an axion undergoing a momentary stage of rapid evolution during inflation. Examining also the scalar perturbations produced by this model, we find that spectral distortions from the scalar modes significantly dominate those arising from the tensors. This holds true also for an earlier version of the model based on a U(1) gauge field. The scalar-induced distortions might be observed in future experiments, and the current COBE/FIRAS constraints already limit the parameter space of these models. Additionally, we find that delaying the onset of fast roll in the SU(2) scenario (to enhance the modes at the scales relevant for spectral distortions, while respecting the CMB constraints at larger scales) poses a greater challenge compared to the U(1) case. We propose a way to control the axion speed by varying the size of its coupling to the gauge fields.

Context. Future X-ray observatories with high spectral resolution and imaging capabilities will enable measurements and mappings of emission line shifts in the intracluster medium (ICM). Such direct measurements can serve as unique probes of turbulent motions in the ICM. Determining the level and scales of turbulence will improve our understanding of the galaxy cluster dynamical evolution and assembly, together with a more precise evaluation of the non thermal support pressure budget. This will allow for more accurate constraints to be placed on the masses of galaxy clusters, among other potential benfits. Aims. In this view, we implemented the methods presented in the previous instalments of our work to characterize the turbulence in the ICM in a feasibility study with the X-IFU on board the future European X-ray observatory, Athena. Methods. From idealized mock observations of a toy model cluster, we reconstructed the second-order structure function built with the observed velocity field to constrain the turbulence. We carefully accounted for the various sources of errors to derive the most realistic and comprehensive error budget within the limits of our approach. With prior assumptions on the dissipation scale and power spectrum slope, we constrained the parameters of the turbulent power spectrum model through the use of MCMC sampling. Results. With favourable assumptions, we were able to retrieve the injection scale, velocity dispersion, and power spectrum slope, with 1sigma uncertainties better than ~15% of the input values. We demonstrated the efficiency of our carefully set framework to constrain the turbulence in the ICM from high-resolution X-ray spectroscopic observations, paving the way for more in-depth investigation of the optimal required observing strategy within a more restrictive observational setup with the future X-IFU instrument.

Earth-based Very Long Baseline Interferometry (VLBI) has made rapid advances in imaging black holes. However, due to the limitations imposed on terrestrial VLBI by the Earth's finite size and turbulent atmosphere, it is imperative to have a space-based component in future VLBI missions. Herein, this paper investigates the effect of Earth's oblateness, also known as the $J_{2}$ effect, on orbiters in Earth-Space and Space-Space VLBI. The paper provides an extensive discussion on how the $J_{2}$ effect can directly impact orbit selection for black hole observations and how through informed choices of orbital parameters, the effect can be used to the mission's advantage, a fact that has not been addressed in existing space-VLBI investigations. We provide a comprehensive study of how the orbital parameters of several current space VLBI proposals will vary specifically due to the $J_{2}$ effect. For black hole accretion flow targets of interest, we have demonstrated how the $J_{2}$ effect leads to modest increase in shorter baseline coverage, filling gaps in the $(u,v)$ plane. Subsequently, we construct a simple analytical formalism that allows isolation of the impact of the $J_{2}$ effect on the $(u,v)$ plane without requiring computationally intensive orbit propagation simulations. By directly constructing $(u,v)$ coverage using the $J_{2}$ affected and invariant equations of motion, we obtain distinct coverage patterns for M87* and SgrA* that show extremely dense coverage on short baselines as well as long term orbital stability on longer baselines.

Fast Radio Bursts (FRBs) are a sensitive probe of the electron distribution in both the large-scale structure and their host galaxies through the dispersion measure (DM) of the radio pulse. Baryonic feedback models are crucial for modelling small scales for ongoing cosmological surveys that are expected to change the electron distribution in galaxies in a way that can be probed by FRB observations. In this paper, we explore the impact of baryonic feedback on FRB hosts using numerical simulations and make a detailed study of the host galaxy dispersion as a function of redshift, galaxy type, feedback model and how these properties vary in independent simulation codes. We find that the host galaxy dispersion varies dramatically between different implementations of baryonic feedback, allowing FRBs with host identification to be a valuable probe of feedback physics and thus provide necessary priors for upcoming analysis of the statistical properties of the large-scale structure. We further find that any dependency on the exact location of events within the halo is small. While there exists an evolution of the dispersion measure with redshift and halo mass, it is largely driven by varying star formation rates of the halo. Spectral information from FRB hosts can therefore be used to put priors on the host galaxy dispersion measure, and FRBs can be used to distinguish between competing models of baryonic feedback in future studies.

We provide an analytical description of the galaxy bispectrum covariance and the power spectrum-bispectrum cross-covariance in redshift space that captures the dominant non-Gaussian contributions. The Gaussian prediction for the variance of the halo bispectrum monopole significantly underestimates numerical estimates particularly for squeezed triangles, that is bispectrum triangular configurations where one side is much smaller than the other two, whereas the effect is relatively less important when considering the quadrupole. We propose an expression for the missing non-Gaussian contribution valid in the squeezed limit that requires an accurate modeling of the bispectrum alone. We validate our model against the numerical covariance estimated from a large suite of mock catalogs and find that it accurately predicts the variance as well as the dominant off-diagonal terms. We also present an expression for the cross-covariance between power spectrum and bispectrum multipoles and likewise find it to provide a good description of the numerical results.

The discrepancies in different measurements of the lifetime of isolated neutrons could be resolved by considering an extra neutron decay channel into dark matter, with a branching ratio of the order of $O(1$\%). Although the decay channel into a dark fermion $\chi$ plus visible matter has been already experimentally excluded, a dark decay with either a scalar or dark photon remains still a possibility. In particular, a model with a fermion mass $m_\chi\approx 1$ GeV and a scalar $m_\phi \approx O(\rm{MeV})$ could provide not only the required branching ratio to explain the anomaly but also a good dark matter (DM) candidate with the right thermal abundance today. Although the interaction DM-neutron will affect the formation of neutron stars, the combined effect of the dark matter self-interactions mediated by the light scalar and an effective repulsive interaction with the neutrons induced by the scalar-Higgs coupling would allow heavy enough neutron stars. The combined constraints from neutron lifetime, dark matter abundance, neutron star and Higgs physics, and Big Bang Nucleosynthesis, restrict the light scalar mass to the range $2 m_e < m_\phi < 2 m_e + 0.0375$ MeV.

Context. Gaia18cjb is one of the Gaia-alerted eruptive young star candidates which has been experiencing a slow and strong brightening during the last 13 years, similar to some FU Orionis-type objects. Aims. The aim of this work is to derive the young stellar nature of Gaia18cjb, determine its physical and accretion properties to classify its variability. Methods. We conducted monitoring observations using multi-filter optical and near-infrared photometry, as well as near-infrared spectroscopy. We present the analysis of pre-outburst and outburst optical and infrared light curves, color-magnitude diagrams in different bands, the detection of near-IR spectral lines, and estimates of both stellar and accretion parameters during the burst. Results. The optical light curve shows an unusually long (8 years) brightening event of 5 mag in the last 13 years, before reaching a plateau indicating that the burst is still on-going, suggesting a FUor-like nature. The same outburst is less strong in the infrared light curves. The near-infrared spectra, obtained during the outburst, exhibit emission lines typical of highly accreting low-intermediate mass young stars with typical EXor features. The spectral index of Gaia18cjb SED classifies it as a Class I in the pre-burst stage and a Flat Spectrum young stellar object (YSO) during the burst. Conclusions. Gaia18cjb is an eruptive YSO which shows FUor-like photometric features (in terms of brightening amplitude and length of the burst) and EXor-like spectroscopic features and accretion rate, as V350 Cep and V1647 Ori, classified as objects in between FUors and EXors

We have used the Green Bank Telescope (GBT) to search for the OH molecule at several locations in the Smith Cloud, one of the most prominent of the high-velocity clouds that surround the Milky Way. Five positions with a high HI column density were selected as targets for individual pointings, along with a square degree around a molecular cloud detected with the Planck telescope near the tip of the Smith Cloud. Gas in the Galactic disk with similar values of $N_{HI}$ has detectable OH emission. Although we found OH at velocities consistent with the foreground Aquila molecular cloud, nothing was found at the velocity of the Smith Cloud to an rms level of 0.7 mK (T$_b$) in a 1 km $s^1$ channel. The three positions that give the strictest limits on OH are analyzed in detail. Their combined data imply a $5\sigma$ limit on $N(H_2) / N_{HI} \leq 0.03$ scaled by a factor dependent on the OH excitation temperature and background continuum $T_{ex}/(T_{ex}-T_{bg})$. There is no evidence for far-infrared emission from dust within the Smith Cloud. These results are consistent with expectations for a low-metallicity diffuse cloud exposed to the radiation field of the Galactic halo rather than a product of a galactic fountain.

Dwarf galaxies in groups of galaxies provide excellent test cases for models of structure formation. This led to a so-called small-scale crisis, including the famous missing satellite and too-big-to-fail problems. It was suggested that these two problems are solved by the introduction of baryonic physics in cosmological simulations. We test for the nearby grand spiral M83 - a Milky Way sibling - whether its number of dwarf galaxy companions is compatible with today's $\Lambda$ + Cold Dark Matter model using two methods: with cosmological simulations that include baryons, as well as with theoretical predictions from the sub-halo mass function. By employing distance measurements we recover a list of confirmed dwarf galaxies within 330 kpc around M83 down to a magnitude of $M_V =-10$. We found that both the state-of-the-art hydrodynamical cosmological simulation Illustris-TNG50 and theoretical predictions agree with the number of confirmed satellites around M83 at the bright end of the luminosity function (>10$^8$ solar masses) but underestimate it at the faint end (down to 10$^6$ solar masses) at more than 3$\sigma$ and 5$\sigma$ levels, respectively. This indicates a too-many satellites problem in $\Lambda$CDM for M83. The actual degree of tension to cosmological models is underestimated, because the number of observed satellites is incomplete due to the high contamination of spurious stars and galactic cirrus.

Multi-object spectroscopic galaxy surveys typically make use of photometric and colour criteria to select targets. Conversely, the Euclid NISP slitless spectrograph will record spectra for every source over its field of view. Slitless spectroscopy has the advantage of avoiding defining a priori a galaxy sample, but at the price of making the selection function harder to quantify. The Euclid Wide Survey aims at building robust statistical samples of emission-line galaxies with fluxes in the Halpha-NII complex brighter than 2e-16 erg/s/cm^2 and within 0.9<z<1.8. At faint fluxes, we expect significant contamination by wrongly measured redshifts, either due to emission-line misidentification or noise fluctuations, with the consequence of reducing the purity of the final samples. This can be significantly improved by exploiting Euclid photometric information to identify emission-line galaxies over the redshifts of interest. To this goal, we compare and quantify the performance of six machine-learning classification algorithms. We consider the case when only Euclid photometric and morphological measurements are used and when these are supplemented by ground-based photometric data. We train and test the classifiers on two mock galaxy samples, the EL-COSMOS and Euclid Flagship2 catalogues. Dense neural networks and support vector classifiers obtain the best performance, with comparable results in terms of the adopted metrics. When training on Euclid photometry alone, these can remove 87% of the sources that are fainter than the nominal flux limit or lie outside the range 0.9<z<1.8, a figure that increases to 97% when ground-based photometry is included. These results show how by using the photometric information available to Euclid it will be possible to efficiently identify and discard spurious interlopers, allowing us to build robust spectroscopic samples for cosmological investigations.

Galaxy velocities in clusters, rotation curves of galaxies, and "vertical" oscillations in the Milky Way currently show too high velocities with respect to the masses thought to be involved. While these velocity excesses are currently interpreted as the consequence of dark matter, it can also be naturally explained as a consequence of scale invariant theory, which rests on a very simple first principle: the addition of a new fundamental symmetry. In the present work, the case of scale invariance, present in General Relativity and Maxwell equations for the empty space without charge and current, is considered. Cosmological models predict a rapid decrease of these effects with increasing mean density up to the critical density, where they totally disappear. Starting from the scale invariant geodesic equation by Dirac (1973), for which a demonstration by an action principle is presented, a modified Newton equation is derived. The solutions of this equation are applied to clusters of galaxies, galactic rotation at different redshifts and "vertical" motions in the Milky Way. In this new framework, the convergence of theoretical predictions and observations, in different gravitational systems, epochs, mass range and spatial scales, opens interesting perspectives that deserve to be explored further.

A simple analytic description is provided of the rate of energy deposition by $\beta$-decay electrons in the homologously expanding radioactive plasma ejected in neutron star mergers, valid for a wide range of ejecta parameters -- initial entropy, electron fraction $\{s_0,Y_e\}$ and density $\rho t^3$. The formulae are derived using detailed numerical calculations following the time-dependent composition and $\beta$-decay emission spectra (including the effect of delayed deposition). The deposition efficiency depends mainly on $\rho t^3$ and only weakly on $\{s_0,Y_e\}$. The time $t_e$ at which the ratio between the rates of electron energy deposition and energy production drops to $1-e^{-1}$, is given by $t_e=t_{0e}\Big(\frac{\rho t^3}{0.5(\rho t^3)_0}\Big)^a$, where $(\rho t^3)_0=\frac{0.05M_{\odot}}{4\pi(0.2c)^3}$, $t_{0e}(s_0,Y_e)\approx17$ days and $0.4\le a(s_0,Y_e)\le0.5$. The fractional uncertainty in $t_e$ due to nuclear physics uncertainties is $\approx10\%$. The result $a\le0.5$ reflects the fact that the characteristic $\beta$-decay electron energies do not decrease with time (largely due to "inverted decay chains" in which a slowly-decaying isotope decays to a rapidly-decaying isotope with higher end-point energy). We provide an analytic approximation for the time-dependent electron energy deposition rate, reproducing the numerical results to better than $50\%$ (typically $<30\%$, well within the energy production rate uncertainty due to nuclear physics uncertainties) over a 3-4 orders-of-magnitude deposition rate decrease with time. Our results may be easily incorporated in calculations of kilonovae light curves (with general density and composition structures), eliminating the need to numerically follow the time-dependent electron spectra. Identifying $t_e$, e.g. in the bolometric light curve, will constrain the (properly averaged) ejecta $\rho t^3$.

A semi-analytic approximation is derived for the time-dependent fraction $f_\gamma(t)$ of the energy deposited by radioactive decay $\gamma$-rays in a homologously expanding plasma of general structure. An analytic approximation is given for spherically symmetric plasma distributions. Applied to Kilonovae (KNe) associated with neutron stars mergers and Type Ia supernovae, our semi-analytic and analytic approximations reproduce, with a few percent and 10% accuracy, respectively, the energy deposition rates, $\dot{Q}_\text{dep}$, obtained in numeric Monte Carlo calculations. The time $t_\gamma$ beyond which $\gamma$-ray deposition is inefficient is determined by an effective frequency-independent $\gamma$-ray opacity $\kappa_{\gamma,\text{eff}}$, $t_\gamma = \sqrt{\kappa_{\gamma,\text{eff}}\langle\Sigma\rangle t^2}$, where $\langle\Sigma\rangle\propto t^{-2}$ is the average plasma column density. For $\beta$-decay dominated energy release, $\kappa_{\gamma,\text{eff}}$ is typically close to the effective Compton scattering opacity, $\kappa_{\gamma,\text{eff}} \approx 0.025~{\rm {cm}^{2}\,g^{-1}}$ with a weak dependence on composition. For KNe, $\kappa_{\gamma,\text{eff}}$ depends mainly on the initial electron fraction $Y_e$, $\kappa_{\gamma,\text{eff}} \approx 0.03(0.05)~{\rm {cm}^{2}\,g^{-1}}$ for $Y_e \gtrsim (\lesssim) 0.25$ (in contrast with earlier work that found $\kappa_{\gamma,\text{eff}}$ larger by 1-2 orders of magnitude for low $Y_e$), and is insensitive to the (large) nuclear physics uncertainties. Determining $t_\gamma$ from observations will therefore measure the ejecta $\langle\Sigma\rangle t^2$, providing a stringent test of models. For $\langle\Sigma\rangle t^2=2\times10^{11}~{\rm g\,{cm}^{-2}\,s^2}$, a typical value expected for KNe, $t_\gamma\approx1$ d.

Advancements in gravitational-wave interferometers, particularly the next generation, are poised to profoundly impact gravitational wave astronomy and multimessenger astrophysics. A hybrid quantum algorithm is proposed to carry out quantum inference of parameters from compact binary coalescences detected in gravitational-wave interferometers. It performs quantum Bayesian Inference with Renormalization and Downsampling (qBIRD). We choose binary black hole (BBH) mergers from LIGO observatories as the first case to test the algorithm, but its application can be extended to more general instances. The quantum algorithm is able to generate corner plots of relevant parameters such as chirp mass, mass ratio, spins, etc. by inference of simulated gravitational waves with known injected parameter values with zero noise, Gaussian noise and real data, thus recovering an accuracy equivalent to that of classical Markov Chain Monte Carlo inferences. The simulations are performed with sets of 2 and 4 parameters. These results enhance the possibilities to extend our capacity to track signals from coalescences over longer durations and at lower frequencies extending the accuracy and promptness of gravitational wave parameter estimation.

In this paper, we use the holographic principle to obtain a modified metric of black holes that reproduces the exponentially corrected entropy. The exponential correction of the black hole entropy comes from non-perturbative corrections. It interprets as a quantum effect which affects black hole thermodynamics especially in the infinitesimal scales. Hence, it may affect black hole stability at the final stage. Then, we study modified thermodynamics due to the non-perturbative corrections and calculate thermodynamics quantities of several non-rotating black holes.

Scalar-tensor theories have a long history as possible phenomenological alternatives to General Relativity, but are known to potentially produce deviations from the (strong) equivalence principle in systems involving self-gravitating objects, as a result of the presence of an additional gravitational scalar field besides the tensor modes of General Relativity. We describe here a novel mechanism whereby the equivalence principle is violated for an isolated rotating neutron star, if the gravitational scalar field is changing in time far from the system. We show that the neutron star rotational period changes due to an effective coupling ("angular momentum sensitivity") to the gravitational scalar, and compute that coupling for viable equations of state for nuclear matter. We comment on the relevance of our findings for testing scalar-tensor theories and models of ultralight dark matter with pulsar timing observations, a topic that we tackle in a companion paper.

Charged leptons produced by high-energy and ultrahigh-energy neutrinos have a substantial probability of emitting prompt internal bremsstrahlung $\nu_\ell + N \rightarrow \ell + X + \gamma$. This can have important consequences for neutrino detection. We discuss observable consequences at high- and ultrahigh-energy neutrino telescopes and LHC's Forward Physics Facility. Logarithmic enhancements can be substantial (e.g.\ $\sim 20\%$) when either the charged lepton's energy, or the rest of the cascade, is measured. We comment on applications involving the inelasticity distribution including measurements of the $\nu/\bar{\nu}$ flux ratio, throughgoing muons, and double-bang signatures for high-energy neutrino observation. Furthermore, for ultrahigh-energy neutrino observation, we find that final state radiation affects flavor measurements and decreases the energy of both Earth-emergent tau leptons and regenerated tau neutrinos. Finally, for LHC's Forward Physics Facility, we find that final state radiation will impact future extractions of strange quark parton distribution functions. Final state radiation should be included in future analyses at neutrino telescopes and the Forward Physics Facility.

Over the past several years, the misidentification of SpaceX Starlink satellites as Unidentified Aerial Phenomena (UAP) by pilots and laypersons has generated unnecessary aviation risk and confusion. The many deployment and orbital evolution strategies, coupled with changing sun specular reflection angles, contribute to this gap in space situational awareness. In this paper we present a case analysis of an incident that generated multiple, corroborating reports of a UAP from five pilots on two commercial airline flights over the Pacific Ocean on August 10th, 2022. This incident included two cell phone photos and a video of an unrecognizable and possibly anomalous phenomenon. We then use supplemental two-line elements (TLEs) for the Starlink train of satellites launched that same day and Automatic Dependent Surveillance Broadcast (ADS-B) data from the flight with the photographs to reconstruct a view of these satellites from the cockpit at the time and place of the sighting. The success of this work demonstrates an approach that could, in principle, warn aviators about satellites that could be visible in unusual or novel illumination configurations, thus increasing space situational awareness and supporting aviation safety. We conclude with recommendations for governments and satellite operators to provide better a-priori data that can be used to create advisories to aviators and the public. The automated simulation of known specular reflection off constellations of satellites could also support researchers investigating sightings of unfamiliar aerial or aerospace objects as likely being from normal versus novel space events.

Shear Alfven wave parametric decay instability (PDI) provides a potential path toward significant wave dissipation and plasma heating. However, fundamental questions regarding how PDI is excited in a realistic three-dimensional (3D) open system and how critically the finite perpendicular wave scale -- as found in both the laboratory and space plasmas -- affects the excitation remain poorly understood. Here, we present the first 3D, open-boundary, hybrid kinetic-fluid simulations of kinetic Alfven wave PDI in low-beta plasmas. Key findings are that the PDI excitation is strongly limited by the wave damping present, including electron-ion collisional damping (represented by a constant resistivity) and geometrical attenuation associated with the finite-scale Alfven wave, and ion Landau damping of the child acoustic wave. The perpendicular wave scale alone, however, plays no discernible role, with different wave scales exhibiting similar instability growth. These findings are corroborated by theoretical analysis and estimates. The new understanding of 3D kinetic Alfven wave PDI physics is essential for laboratory study of the basic plasma process and may also help evaluate the relevance/role of PDI in low-beta space plasmas.

The present paper is devoted to a study of the equilibrium configurations of slowly rotating anisotropic stars in the framework of general relativity. For that purpose, we provide the equations of structure where the rotation is treated to second order in the angular velocity. These equations extend those first derived by Hartle for slowly rotating isotropic stars. As an application of the new formalism, we study the rotational properties of Bowers-Liang fluid spheres. A result of particular interest is that the ellipticity and mass quadrupole moment are negative for certain highly anisotropic configurations, thus such systems are prolate rather than oblate. Furthermore, for configurations with high anisotropy, and compactness close to their critical value, quantities like the moment of inertia, change of mass, and mass quadrupole moment approach to the corresponding Kerr black hole values, similar to other ultracompact systems like sub-Buchdahl Schwarzschild stars and analytic rotating gravastars.

Within the framework of metric-affine theories of gravity, where both the metric and connection are treated as independent variables, we consider actions quadratic in the Ricci scalar curvature coupled non-minimally to a scalar field through derivative couplings. Our analysis delves into the inflationary predictions, revealing their consistency with the latest observational constraints across a wide range of parameters. This compatibility permits adjustments such as an increase in the spectral index and a reduction in the tensor-to-scalar ratio. While we do not propose a specific reheating mechanism, our analysis demonstrates that within the quadratic model of inflation, the maximum reheating temperature can reach $\sim 3\times10^{15}\, {\rm GeV}$.

We present the results of the charge ratio ($R$) and polarization ($P^{\mu}_{0}$) measurements using the decay electron events collected from 2008 September to 2022 June by the Super-Kamiokande detector. Because of its underground location and long operation, we performed high precision measurements by accumulating cosmic-ray muons. We measured the muon charge ratio to be $R=1.32 \pm 0.02$ $(\mathrm{stat.}{+}\mathrm{syst.})$ at $E_{\mu}\cos \theta_{\mathrm{Zenith}}=0.7^{+0.3}_{-0.2}$ $\mathrm{TeV}$, where $E_{\mu}$ is the muon energy and $\theta_{\mathrm{Zenith}}$ is the zenith angle of incoming cosmic-ray muons. This result is consistent with the Honda flux model while this suggests a tension with the $\pi K$ model of $1.9\sigma$. We also measured the muon polarization at the production location to be $P^{\mu}_{0}=0.52 \pm 0.02$ $(\mathrm{stat.}{+}\mathrm{syst.})$ at the muon momentum of $0.9^{+0.6}_{-0.1}$ $\mathrm{TeV}/c$ at the surface of the mountain; this also suggests a tension with the Honda flux model of $1.5\sigma$. This is the most precise measurement ever to experimentally determine the cosmic-ray muon polarization near $1~\mathrm{TeV}/c$. These measurement results are useful to improve the atmospheric neutrino simulations.

We present a study of energy density and pressure of a free real scalar quantum field after its decoupling from a thermal bath in the spatially flat Friedman-Lemaitre-Robertson-Walker space-time by solving the Klein-Gordon equation both analytically and numerically for different predetermined scale factor functions $a(t)$. The energy density and pressure, defined by subtracting the vacuum expectation values at the decoupling time, feature corrections with respect to the classical free-streaming solution of the relativistic Boltzmann equation. We show that if the expansion rate is comparable or larger than $mc^2/\hbar$ or $KT_0/\hbar$ where $m$ is the mass and $T_0$ the decoupling temperature, both energy density and pressure gets strong quantum corrections which substantially modify their classical dependence on the scale factor $a(t)$ and drive pressure to large negative values. For a minimally coupled field with a very low mass in an expanding de Sitter universe quantum corrections are dominant driving pressure and energy density to become asymptotically constant with an equation of state $p/\varepsilon \simeq -1$, thereby mimicking a cosmological constant. For a minimally coupled massless field, quantum corrections are asymptotically dominant for any accelerated expansion.

The recent increasing interest in detecting gravitational waves (GWs) by lunar seismic measurement urges us to have a clear understanding of the response of the moon to passing GWs. In this paper, we clarify the relationship between two seemly different response functions which have been derived previously using two different methods, one taking the field-theory approach and the other using the tidal force induced by GWs. We revisit their derivation and prove, by both analytical arguments and numerical calculations, that the two response functions are equivalent. Their apparent difference can be attributed to the choice of different coordinates. Using the correct response function, we calculate the sensitivities (to GWs) of several designed lunar seismometers, and find that the sensitivity curves between $10^{-3}$ and $0.1$ Hz are much flatter than the previous calculations based on normal-mode model. Our results will help clarifying the scientific objectives of lunar GW observation, as well as provide important constraints on the design of lunar GW detectors.