Papers for Tuesday, May 06 2025

In this paper we present the latest results of our ongoing multiplicity survey of (Community) TESS Objects of Interest, using astrometry and photometry from the latest data release of the ESA Gaia mission to detect stellar companions of these stars and to characterize their properties. A total of 92 binary and two hierarchical triple star systems are identified among the 745 target stars whose multiplicity is explored in this study, all at distances of less than 500pc around the Sun. As expected for components of gravitationally bound star systems, the targets and the detected companions are at the same distance and share a common proper motion, as shown by their accurate Gaia astrometry. The companions have masses of about 0.12 to 1.6$M_\odot$. and are most frequently found in the mass range up to 0.6$M_\odot$. The companions have projected separations from the targets between about 110 and 9600au. Their frequency is highest and constant from about 300 to 800au, decreasing at larger projected separations. In addition to main sequence stars, five white dwarf companions are detected in this study, whose true nature is unveiled by their photometric properties.

We present XRISM Resolve observations centered on Hydra-A, a redshift z = 0.054 brightest cluster galaxy which hosts one of the largest and most powerful FR-I radio sources in the nearby Universe. We examine the effects of its high jet power on the velocity structure of the cluster's hot atmosphere. Hydra-A's central radio jets have inflated X-ray cavities with energies upward of $10^{61}~\rm erg$. They reach altitudes of 225 kpc from the cluster center, well beyond the atmosphere's central cooling region. Resolve's $3\times3$ arcmin field-of-view covers $190\times190$ kpc, which encompasses most of the cooling volume. We find a one dimensional atmospheric velocity dispersion across the volume of $164\pm10\,\rm{km\,s}^{-1}$. The fraction in isotropic turbulence or unresolved bulk velocity is unknown. Assuming pure isotropic turbulence, the turbulent kinetic energy is $2.5\,\%$ of the thermal energy radiated away over the cooling timescale, implying that kinetic energy must be supplied continually to offset cooling. While Hydra-A's radio jets are powerful enough to supply kinetic energy to the atmosphere at the observed level, turbulent dissipation alone would struggle to offset cooling throughout the cooling volume. The central galaxy's radial velocity is similar to the atmospheric velocity, with an offset of $-37 \pm 23\,$ km s$^{-1}$.

The radii and masses of many giant exoplanets imply their interiors each contain more than $\sim$100 $M_\oplus$ of solids. A large metal content may arise when a giant planet grows by colliding and merging with multiple $\sim$10 $M_\oplus$ solid cores. Here we show that a giant impact with a young gas giant excites long-lived seismic oscillations that can be detected photometrically. Mode lifetimes are close to the planet's Kelvin-Helmholtz time, a significant fraction of a young planet's age. Oscillation periods lie between tens of minutes to an hour, and variability amplitudes can exceed a percent for several million years. Beta Pictoris b is a young super-Jupiter known to be highly metal-enriched. If a Neptune-mass (17 $M_\oplus$) body impacted $\beta$ Pictoris b in the past $\sim$9--18 Myr, the planet could still be ringing with a percent-level photometric variability measurable with JWST.

Mass loss in massive stars is crucial to understanding how these stars evolve and explode. Despite increasing evidence indicating its importance, episodic mass loss remains poorly understood. Here we report the results of an optical spectroscopic survey of evolved massive stars in NGC 6822, IC 10, and IC 1613, conducted by the ASSESS project (Episodic Mass Loss in Evolved Massive Stars: Key to Understanding the Explosive Early Universe), which aimed to investigate the role of episodic mass loss, by targeting stars with infrared excesses indicating a dusty circumstellar environment. We assigned a spectral class to 122 unique sources, the majority of which are dusty. The rate of evolved massive stars was over 60% for the highest-priority targets. We discovered 2 blue supergiants, 1 yellow supergiant, 1 emission-line object, and confirmed 2 supernova remnant candidates, a Wolf-Rayet star, and 2 H II regions. Twenty-eight unique sources were classified as red supergiants, 21 of which are new discoveries. In IC 10, we increased the sample of spectroscopically confirmed RSGs from 1 to 17. We used the MARCS models to obtain their surface properties, most importantly the effective temperature, and spectral energy distribution fitting to obtain the stellar luminosity for 17 of them. The dusty RSGs are cooler, more luminous, more extinguished, and more evolved than the non-dusty ones, in agreement with previous findings. Investigating the optical photometric variability of the RSGs from light curves covering a period over a decade revealed that the dusty RSGs are more variable. We further highlight the very extinguished emission-line object, two RSGs that display a significant change in spectral type between two observed epochs, and four dusty K-type RSGs as candidates for having undergone episodic mass loss.

We report the discovery of a Lyman $\alpha$ emitter (LAE) candidate in the immediate foreground of the quasar PSO J158-14 at $z_{\rm QSO}=6.0685$ at a projected distance $\sim29\ {\rm pkpc}$ that is associated with an extremely metal-poor absorption system. This system was found in archival observations of the quasar field with the Very Large Telescope/Multi-Unit Spectroscopic Explorer (VLT/MUSE) and was previously missed in searches of absorption systems using quasar absorption line spectroscopy as it imparts no detectable metal absorption lines on the background quasar spectrum. The detected Ly$\alpha$ emission line at a redshift of $z_{\rm LAE}=6.0323$ is well aligned with the outer edge of the quasar's proximity zone and can plausibly cause its observed damping wing if it is associated with a proximate sub-damped Ly$\alpha$ absorption system with a column density of $\log {N_{\rm HI} / {\rm cm}^{-2}} \approx 19.8$. A $>10$ hour medium-resolution spectrum of the quasar observed with the Magellan/Folded-port InfraRed Echellette (FIRE) and VLT/X-Shooter spectrographs reveals a metallicity constraint of ${\rm [Z/H]} < -3$. Such low metallicity makes this system an extremely metal-poor galaxy candidate and provides an exciting site to study possible signatures of Population III stars.

Galaxy groups and clusters are one of the best probes of structure formation and growth in a cosmological context. Most of their baryonic component is dominated by the intracluster medium (ICM), whose thermodynamical properties serve as indicators of the halo's dynamical state and can be used for the halo mass determination in the self-similar scenario. However, baryonic processes, such as AGN feedback and gas cooling, may affect the global properties of the ICM, especially in the group regime. These effects might lead to deviations from self-similar predictions in galaxy groups' scaling relations, while they remain in place for massive galaxy clusters. Additionally, the low-mass end of the scaling relations, ranging from $10^{13}$ to $10^{14} M_\odot$, remains unclear and poorly populated, as current X-ray surveys detect only the brightest groups. Here, we present the Mass-Temperature relation across the full mass range, from massive clusters to low-mass groups ($10^{13}M_\odot$), as observed by eROSITA. Using spectral stacking from eROSITA eRASS1 data for optically selected galaxy groups, we find that, in the lower mass range, galaxy groups follow the power-law relation known for galaxy clusters. We further validate these results by conducting the same stacking procedure on mock eRASS:4 data using the Magneticum hydrodynamical simulation. This indicates that AGN feedback is more likely to affect the distribution of baryons in the intragroup medium rather than the overall halo gas temperature. No significant changes in the Mass-Temperature relation slope suggest that temperature can serve as a reliable mass proxy across the entire mass range. This validates the use of temperature-derived masses, particularly in cosmological studies, significantly broadening the mass range and enabling applications such as improving the cluster mass function studies and cosmological parameter estimates.

As dark matter appears to comprise most of the Galactic mass, some of it may accumulate in the cores of stars, thereby making the Sun a laboratory for constraining various dark matter theories. We consider the effects on the solar structure arising from a general class of macroscopic dark matter candidates that include strange quark matter, compact dark objects, and others. We calibrate standard solar evolution models (i.e., models that reproduce the mass, luminosity, radius, and metallicity of the Sun at its present age) with variable compact dark core masses ranging from $10^{-8}$ to $10^{-2}~\rm{M}_\odot$ and assess their properties. We find that the weakest constraints come from solar neutrino flux measurements, which only rule out the most massive dark core comprising at least $\sim 1\%$ of the total solar mass. The Sun's acoustic oscillations impose stronger constraints, probing masses down to $\sim 10^{-5}~\rm{M}_\odot$. We find that a model with a $10^{-3}~\rm{M}_\odot$ dark core appears to improve the agreement with helioseismic observations. We nevertheless caution against interpreting this as evidence for dark matter in the solar interior, and suggest plausible effects that the dark core may instead be emulating. Finally, we show that future measurements of solar $g$~modes may constrain dark core masses down to $10^{-7}~\rm{M}_\odot$.

The evolution of binary stellar systems involves a wide range of physical processes, many of which are not yet well understood. We aim to build a general-purpose algorithm based on inverse population synthesis techniques, able to reconstruct the past history of binary systems. This algorithm will be applied to a sample of eclipsing binaries, aiming to ascertain their progenitors and past histories. Once validated, it was applied to a sample 30 white dwarf plus main-sequence eclipsing binaries observed by the Zwicky Transient Facility survey. We determined the input space parameters of the progenitors for the 30 eclipsing binary systems to which the algorithm was applied. These parameters included the initial primary and secondary masses, the orbital separation and eccentricity, the common-envelope efficiency ($\alpha_{\rm CE}$), and the age at which the system was formed. Furthermore, the analysis of the global properties revealed some important features: a mild anticorrelation between the common-envelope efficiency parameter and the secondary mass, the absence of a universal value of $\alpha_{\rm CE}$ along with no need for internal energy, although in the low-mass regime, the high values of $\alpha_{\rm CE}$ suggest a possible contribution, and an initial thermalized eccentricity distribution. Although a strong degeneracy among the input parameters exists in the reconstruction of post-common envelope binary systems, the high accuracy obtained for the eclipsing-binary systems analyzed here has allowed our algorithm to make a reasonable determination of the initial parameters without the need to include external constraints. The global properties found here so far, can be substantially improved when analyzing a future volume-complete sample.

The discovery of inhabited exoplanets hinges on identifying biosignature gases. JWST is revealing potential biosignatures in exoplanet atmospheres, though their presence is yet to provide strong evidence for life. The central challenge is attribution: how to confidently identify biogenic sources while ruling out, or deeming unlikely, abiotic explanations? Attribution is particularly difficult for individual planets, especially regarding system-scale stochastic processes that could set atmospheric conditions. To address this, we here propose a comparative multi-planet approach: comparing atmospheric compositions across multiple planets within a system to empirically define the 'abiotic baseline'. This baseline serves as a reference point for biosignatures, and enables marginalisation over inaccessible, shared abiotic parameters. This is possible because planets within a system are linked by their birth in the same natal disk, having been irradiated by the same evolving star, and having a related dynamical history. Observations aligning with the abiotic baseline, where the locally informed abiotic models demonstrate high out-of-sample predictive accuracy, are likely non-biological. Deviations from the baseline -- potentially biotic anomalies -- suggest an alternative origin. We present the application of Bayesian leave-one-out cross-validation to evaluate the performance of geochemical- and biogeochemical-climate models in explaining these anomalies, using the expected log pointwise predictive density as a diagnostic. When biogeochemical models outperform their abiotic counterparts, the anomaly may be shaped by life, and constitutes a comparative biosignature. If both models perform poorly, the anomaly is flagged as an "unknown unknown" -- a signature of either unrecognised abiotic chemistry, or life as we don't yet know it.

We investigate the observational constraints on $\alpha$-attractor inflationary models and their post-inflationary reheating dynamics in light of the latest CMB data from ACT DR6 combined with Planck18, BICEP/$Keck$ 2018, and DESI (collectively denoted P-ACT-LB-BK18). Focusing on both E and T type attractor potentials, we analyze how inflationary observables namely the scalar spectral index $n_s$ and the tensor-to-scalar ratio $r$ are indirectly influenced by reheating parameters such as the reheating temperature $T_{\text{RH}}$, the inflaton equation-of-state $w_\phi$, and the inflaton's couplings to Starndard Model particles. We incorporate indirect constraints from the overproduction of primordial gravitational waves (PGWs), particularly via $\Delta N_{\rm eff}$ bounds on BBN, which become significant for stiff post-inflationary dynamics. Our results show that E-models permit a wide range of reheating scenarios, including matter-like reheating ($w_\phi = 0$), while T-models are viable only for $w_\phi \gtrsim 0.44$. We derive bounds on inflaton-Standard Model couplings for both decay and scattering channels and identify parameter regimes compatible with recent ACT data for successful reheating. These findings establish a robust connection between inflationary theory, thermal history, and particle phenomenology, and provide predictive targets for future CMB missions.

We analyse jointly phase-curve observations of the ultra-hot Jupiter WASP-18 b from the visible to the mid-infrared, using unpublished data from CHEOPS, TESS and Spitzer. We aim to characterise the planetary atmosphere with a consistent view over the large wavelength range covered using GCMs and retrieval analyses, and including JWST data. We obtain new ephemerides with unprecedented precisions of 1 second and 1.4 millisecond on the time of inferior conjunction and orbital period, respectively. We compute a planetary radius of $R_p = 1.1926 \pm 0.0077 R_J$ with a precision of 0.65%, or 550 km. Through a timing inconsistency with JWST, we discuss and confirm orbital eccentricity ($e = 0.00852 \pm 0.00091$), and constrain the argument of periastron to $\omega = 261.9^{+1.3}_{-1.4}$ deg. We show that the large dayside emission implies the presence of magnetic drag and super-solar metallicity. We find a steep thermally-inverted gradient in the planetary atmosphere, which is common for UHJs. We detect the presence of strong CO emission lines at 4.5 $\mu$m from an excess of dayside brightness in the Spitzer/IRAC/Ch2 passband. Using these models to constrain the reflected contribution in the CHEOPS passband, we derive an extremely low geometric albedo $A_g^\text{CHEOPS} = 0.027 \pm 0.011$.

The transfer of a star's angular momentum to its atmosphere is a topic of considerable and wide-ranging interest in astrophysics. This letter considers the effect of kinetic and magnetic turbulence on the solar wind's angular momentum. The effects are quantified in a theoretical framework that employs Reynolds-averaged mean field magnetohydrodynamics, allowing for fluctuations of arbitrary amplitude. The model is restricted to the solar equatorial (\(r-\phi\)) plane with axial symmetry, which permits the effect of turbulence to be expressed in analytical form as a modification to the classic Weber & Davis (1967) theory, dependent on the \(r,\phi\) shear component of the Reynolds stress tensor. A solar wind simulation with turbulence transport modeling and Parker Solar Probe observations at the Alfvén surface are employed to quantify this turbulent modification to the solar wind's angular momentum, which is found to be ~ 3% - 10% and tends to be negative. Implications for solar and stellar rotational evolution are discussed.

The sub-Jovian exoplanet WASP-107b ranks among the best-characterized low-density worlds, featuring a Jupiter-like radius and a mass that lies firmly in the sub-Saturn range. Recently obtained JWST spectra reveal significant methane depletion in the atmosphere, indicating that WASP-107b's envelope has both a high metallicity and an elevated internal heat flux. Together with a detected non-zero orbital eccentricity, these data have been interpreted as evidence of tidal heating. However, explaining the observed luminosity with tidal dissipation requires an unusually low tidal quality factor of $Q \sim 100$. Moreover, we find that secular excitation by the RV-detected outer companion WASP-107c, generally cannot sustain WASP-107b's eccentricity in steady state against tidal circularization. As an alternative explanation, we propose that Ohmic dissipation -- generated by interactions between zonal flows and the planetary magnetic field in a partially ionized atmosphere -- maintains the observed thermal state. Under nominal assumptions for the field strength, atmospheric circulation, and ionization chemistry, we show that Ohmic heating readily accounts for WASP-107b's inflated radius and anomalously large internal entropy.

This study presents the results of ALMA band 6 archival data of G328.24$-$0.55, with the aim to pin down the physical and kinematic properties of young stellar objects (YSOs) in G328.24$-$0.55 star forming region. The dust continuum image reveals 5 protostellar objects (MM1a, MM1b, MM1c, MM2 and MM3), with MM1a dominating the region. The dust continuum peaks do not coincide with the strongest radio continuum peak previously detected in the region in a MeerKAT observation, but coincide with the weaker MeerKAT peak. The dust continuum objects are associated with faint unresolved infrared emission. We detected 70, 49, 26, 7 and 8 molecular transitions toward MM1a, MM1b, MM1c, MM2 and MM3, respectively. This variation in the number of detected molecular transitions supports different excitation conditions in these objects. The excitation temperatures estimated toward MM1a, MM1b and MM1c are $\sim$ 183, 168 and 110\,K, respectively. MM2 and MM3 lack multiple transitions of molecular lines to determine their excitation temperatures. The masses of MM1a, MM1b, MM1c, MM2 and MM3 were calculated to be 23.2, 16.1, 12.0, 9.8 and 14.9$M_{\odot}$, respectively. The velocity gradient of CH$_{3}$OH ($10_{2,8}-9_{3,7}$) emission traces a rotating structure, probably an envelope of gas around MM1a. Bipolar outflow traced by CO emission is seen towards MM1a. The properties of MM1a clearly point to the existence of a massive protostellar object that is still undergoing accretion and outflow in its early formative stage.

Galaxy clusters trace the most massive dark matter haloes, whose shapes and orientations reflect the imprint of the cosmic large-scale tidal field. This paper introduces the Subhalo-based Halo Alignment and Projected Ellipticity (SHAPE) technique, which reconstructs cluster halo shapes from the projected distribution of subhaloes, providing a novel approach to investigate intrinsic alignment (IA) correlations between cluster halo shapes and the surrounding density field. We measure halo shapes and orientations using different line-of-sight projection depths and find that, with modest projection depths, the shapes and orientations recovered by SHAPE show good agreement with those measured directly from the simulation particles. Using these SHAPE-derived shapes, we compute IA correlation functions from N-body simulations in both real and redshift space. The IA correlation multipoles exhibit features consistent with baryon acoustic oscillations around 100 Mpc/h and show redshift-space distortion (RSD) effects that agree well with predictions from a non-linear alignment model incorporating RSD. We further demonstrate that the structure growth rate parameter can be robustly estimated without bias from these IA correlations, providing a new avenue for cosmological parameter estimation. Expanding the IA correlations in an associated Legendre basis yields results consistent with those from the standard Legendre expansion, but with improved statistical significance. These results suggest that SHAPE may enhance cosmological parameter constraints in future galaxy surveys.

The Dark Energy Spectroscopic Instrument (DESI) survey uses an automatic spectral classification pipeline to classify spectra. QuasarNET is a convolutional neural network used as part of this pipeline originally trained using data from the Baryon Oscillation Spectroscopic Survey (BOSS). In this paper we implement an active learning algorithm to optimally select spectra to use for training a new version of the QuasarNET weights file using only DESI data, specifically to improve classification accuracy. This active learning algorithm includes a novel outlier rejection step using a Self-Organizing Map to ensure we label spectra representative of the larger quasar sample observed in DESI. We perform two iterations of the active learning pipeline, assembling a final dataset of 5600 labeled spectra, a small subset of the approx 1.3 million quasar targets in DESI's Data Release 1. When splitting the spectra into training and validation subsets we meet or exceed the previously trained weights file in completeness and purity calculated on the validation dataset with less than one tenth of the amount of training data. The new weights also more consistently classify objects in the same way when used on unlabeled data compared to the old weights file. In the process of improving QuasarNET's classification accuracy we discovered a systemic error in QuasarNET's redshift estimation and used our findings to improve our understanding of QuasarNET's redshifts.

Neutron star (NS) mergers, including both binary NS mergers and black hole - NS mergers, are multi-messenger sources detectable in both gravitational waves (GWs) and electromagnetic (EM) radiation. The expected EM emission signatures depend on the source's progenitor, merger remnant, and observer's line of sight (LoS). Widely discussed EM counterparts of NS mergers have been focused in the gamma-ray (in terms of short-duration gamma-ray bursts) and optical (in terms of kilonova) bands. In this paper, we show that X-ray emission carries unique post-merger signatures that are inaccessible in other EM bands and plays an important role in understanding the physics and geometry of these mergers. We consider several binary progenitor and central engine models and investigate X-ray emission signatures from the prompt phase immediately after the merger to the afterglow phase extending years later. For the prompt phase, we devise a general method for computing phenomenological X-ray light curves and spectra for structured jets viewed from any LoS, which can be applied to X-ray observations of NS mergers to constrain the geometry. The geometric constraints can in turn be used to model the afterglow and estimate a peak time and flux -- to preemptively determine afterglow characteristics would be monumental for follow-up observation campaigns of future GW sources. Finally, we provide constraints on the time window for X-ray counterpart searches of NS mergers across a range of luminosity distances and detector sensitivities.

We use the cosmological hydrodynamical simulation TNG50 to study the galaxy mass-morphology relation, as measured by the rotational support of the stellar component of simulated galaxies. For isolated galaxies with stellar mass in the range $8<\log(M_{*}/M_{\odot}) < 11$, rotational support increases with $M_*$, from dispersion-supported spheroidal dwarfs to massive galaxies with prominent rotationally supported discs. Our results indicate that this correlation arises from the spatial distribution of star formation in TNG50 galaxies, which occurs primarily in two distinct regions: an unresolved, non-rotating central baryonic clump $(r \lesssim 1~\mathrm{kpc})$ and a rotationally supported outer disc, separated by a quiescent region. The importance of the inner clump increases with decreasing $M_*$; it makes up more than 80\% of all stars in dwarfs. This explains why dwarfs have less rotational support than massive galaxies, as well as why all dwarfs have similar stellar half-mass radii regardless of $M_*$. It also explains why dwarfs in TNG50 appear to form outside-in, as star formation in the dominant inner clump moves progressively inward. The clump-disc segregation of star formation in TNG50 galaxies is probably numerical in origin. Inner clumps are formed by the accumulation of low-angular momentum gas supported by the equation of state introduced to prevent artificial fragmentation. The decoupled-wind feedback implementation in TNG50 helps to preserve the clumps but disrupts disc formation in its immediate surroundings. This hinders the formation of discs in (dwarf) galaxies whose sizes are not substantially larger than the clump, but has little effect on the larger discs of more massive systems. Our results argue for caution when interpreting the dependence on stellar mass of TNG50 galaxy morphologies, or the evolution of galaxy sizes, especially at the dwarf end.

Active Galactic Nuclei (AGNs) are among the most luminous objects in the universe, making them valuable probes for studying galaxy evolution. However, understanding how AGN properties evolve over cosmic time remains a fundamental challenge. This study investigates whether AGNs at low redshift (nearby) can serve as proxies for their high-redshift (distant) counterparts by identifying spectral 'doppelgängers', AGNs with remarkably similar emission line properties despite being separated by vast cosmic distances. We analyze key spectral features of bona fide AGNs using the Sloan Digital Sky Survey's Data Release 16, including continuum and emission lines: Nitrogen (N V), Carbon (C IV), Magnesium (Mg II), Hydrogen-beta (H$\beta$), and Iron (Fe II - optical and UV) emission lines. We incorporated properties such as equivalent width, velocity dispersion in the form of full width at half maximum (FWHM), and continuum luminosities (135nm, 300nm, and 510nm) closest to these prominent lines. Our initial findings suggest the existence of multiple AGNs with highly similar spectra, hinting at the possibility that local AGNs may indeed share intrinsic properties with high-redshift ones. We showcase here one of the better candidate pairs of AGNs resulting from our analyses.

The Hubble tension has become one of the central problems in cosmology. In this work, we determine the Hubble constant $H_0$ and sound horizon $r_d$ by using the combination of Baryon Acoustic Oscillations (BAOs) from DESI surveys, time-delay lensed quasars from H0LiCOW collaborations and the Pantheon supernovae observations. We consider two cosmological approaches, i.e., Taylor series and Padé polynomials, to avoid cosmological dependence. The reason for using this combination of data is that the absolute distance provided by strong gravitational lensing helps anchor the relative distance of BAO, and supernovae provide a robust history of universe evolution. {{Combining the 6 time-delay distance (6$D_{\Delta t}$) plus 4 angular diameter distance to the deflector (4$D_d$) measurements of time-delay lensed quasars,}} the BAO and the type Ia of supernovae (SNe Ia) datasets, we obtain a model-independent result of $r_d = 138.2_{-3.9}^{+3.3}$ Mpc and $H_0 = 72.9^{+1.8}_{-1.8}$ ${\mathrm{~km~s^{-1}~Mpc^{-1}}}$ for the Taylor series cosmography and $r_d = 137.0_{-3.7}^{+3.2}$ Mpc and $H_0 = 73.1_{-1.7}^{+1.8}$ ${\mathrm{~km~s^{-1}~Mpc^{-1}}}$ for the Padé polynomials cosmography. The determination of $r_d$ and $H_0$ prefers larger $H_0$ and smaller $r_d$ than Planck data under the assumption of flat-$\Lambda$CDM model. However, the values of $H_0$ are consistent with the $H_0$ determination from SH0ES collaboration.

In the Brief History of Time Stephen Hawking was pessimistic about astronomers detecting primordial black holes (PBHs). He would not be the only distinguished scientist to underestimate the extraordinary power of new technology. In a related area Albert Einstein published the equations for microlensing, but wrote off their practicality. Perhaps they meant "during my lifetime." The amazing properties of PBHs, however, validate heroic efforts to detect them. If they exist, their niches in our current history of time include supplying dark matter to bind galaxies, offering a solution for the Hubble tension, and, as supermassive black holes, giving us quasars as far as the eye can see. This Research Note describes a search for PBHs in the Gaia archive. In spite of the high density of local dark matter, it was unsuccessful. Microlensing with the Rubin telescope is the tool at our disposal to open the asteroid window for PBH.

We present high-resolution observations made with the Australia Telescope Compact Array (ATCA) at 5.5 GHz and 9.0 GHz of a sample of twenty-eight massive protostars. These protostars were initially detected at radio wavelengths in the MeerKAT Galactic Plane Survey at 1.3 GHz, where they were unresolved with an angular resolution of approximately 8 arcseconds. The resolution of the ATCA observations at 5.5 GHz and 9.0 GHz are significantly higher at, 2.0 and 1.2 arc seconds, respectively. The highest angular resolution corresponds to a linear resolution of 1920 AU for our nearest object. This improvement in resolution enabled us to resolve the components of extended emission and to more accurately determine the nature of the objects. Approximately 80% of the detections were classified as jets with 68% of the jets found to be associated with non-thermal emission.

We present an investigation of the L-band emission from known massive young stellar objects (MYSOs) in the SARAO MeerKAT Galactic Plane Survey to search for non-thermal radio emitters in the sample. A total of 398 massive protostars, identified from the Red MSX Source (RMS) survey, are located within the survey region. Among these, 162 fields that host the protostars are isolated from nearby bright HII regions, allowing for the study of any ionized jets present. Seventy-one of these fields have jets with five-sigma detections or higher, corresponding to a detection rate of 44%. The MeerKAT fluxes of the detections, together with the upper limits of the non-detections and any other fluxes from previous observations, were used to estimate the spectral indices of the jets, and to search for the presence of non-thermal radiation. In cases where a source manifests as single in a given observation but is resolved into multiple components in observations of higher resolutions, the sum of the fluxes of the resolved components was used in estimating the indices. Any e!ects from missing flux in higher-resolution observations were incorporated into the index uncertainties. The spectral indices of the sample show that at least 50% of the jets emit non-thermal radiation. Additionally, the spectral energy distribution (SED) of some of the sources, as well as their radio luminosities exhibit evidence of non-thermal emission, especially in extended sources.

In this paper, we carry out multiwavelength and multiview observations of the eruption of an intermediate prominence originating from the farside of the Sun on 2023 March 12. The southeast footpoint of the prominence is located in active region (AR) 13252. The eruption generates a B7.8 class flare and a partial halo coronal mass ejection (CME). The prominence takes off at 02:00 UT and accelerates for nearly three hours. Rotation of the southeast leg of the prominence in the counterclockwise direction is revealed by spectroscopic and imaging observations. The apex of the prominence changes from a smooth loop to a cusp structure during the rising motion and the northwest leg displays a drift motion after 04:30 UT, implying a writhing motion. Hence, the prominence eruption is most likely triggered by ideal kink instability. For the first time, we apply the Graduated Cylindrical Shell (GCS) modeling in three-dimensional reconstruction and tracking of the prominence for nearly two hours. Both the source region (110$\degr$E, 43$\degr$N) and northwest footpoint (162$\degr$E, 44$\degr$N) are located. The edge-on and face-on angular widths of the prominence are $\sim$6$\degr$ and $\sim$86$\degr$, respectively. The axis has a tilt angle of $\sim$70$\degr$ with the meridian. The heliocentric distance of the prominence leading edge increases from $\sim$1.26\,$R_{\sun}$ to $\sim$2.27\,$R_{\sun}$. The true speed of the CME increases from $\sim$610 to $\sim$849 km s$^{-1}$.

In this paper, the minimum orbit intersection distances (MOIDs) of near-Earth asteroids (NEAs) over the next 200 years were computed and analyzed in detail. It was shown that the MOID of a NEA relative to the Earth-Moon barycenter (EMB) is usually a superior metric for predicting a potential impact than that relative to the Earth. Subsequently, a novel MOID Evolution Index (MEI) spanning from 0.0 to 9.9 was proposed and the orbits of NEAs are classified into 100 distinct categories by considering the variations of the MOID over time, which is useful for quickly screening and prioritizing hazardous asteroids for future research. Furthermore, it was demonstrated that a linear fitting to the MOID evolution provides a simple yet valid approach for most of the NEAs, which is useful for quickly estimating the MOID value without the need to perform an orbit propagation. As a result, a scheme with several parameters was proposed to characterize the MOID variations as well as the relative position information of the critical points along the orbits associated with the minimum distances. A database incorporating these parameters and the MEI values was therefore established for the cataloged NEAs, enabling the derivation of statistically constrained upper bounds for secular MOID drift rate as function of the semi-major axes. Finally, some special orbital configurations and dynamical mechanisms that may lead to a large deviation from the linear fit or multiple orbit crossings were also investigated, indicating the intricate nature in the patterns of MOID evolution for some NEAs.

Millisecond pulsars (MSPs) are prominent GeV $\gamma$-ray emitters, but their emission mechanisms, especially during off-pulse phases and at high energies, remain debated. We analyze 38 MSPs using approximately 15 years of $Fermi$-LAT data. Employing high-precision ephemerides and Bayesian Blocks algorithm, we detect significant off-pulse emission from 15 MSPs. Spectral analysis reveals significant cutoffs in 10 of these, favoring a magnetospheric origin. The remaining 5 likely also have magnetospheric emission with cutoffs below current sensitivity. Phase-resolved spectral analysis for the 15 MSPs with off-pulse flux shows a significant correlation between cutoff energy ($E_{\text{cut}}$) and average photon counts in 11 sources. We introduce a phase-resolved pseudo-luminosity ($L$) and discover a novel, statistically significant linear correlation between $\log_{10}(L)$ and $\log_{10}(E_{\text{cut}})$ across their phase bins, with the best-fit slope $\alpha = 2.31^{+0.22}_{-0.25}$. Furthermore, we detect significant pulsed emission above 10 GeV from 9 MSPs with off-pulse emission and 7 others (totaling 16 MSPs and tripling previous cataloged detections). Emission above 25 GeV is confirmed for PSR J0614$-$3329 and newly detected for PSR J1536$-$4948, both of which also exhibit significant off-pulse emission. Their highest-energy photons ($\sim$61 GeV and $\sim$57 GeV, respectively) coincide with pulse peaks. The co-existence of significant off-pulse emission and high-energy pulsed emission challenges standard Outer Gap models, while the observed $\log_{10}(L)$-$\log_{10}(E_{\text{cut}})$ correlation may provide new evidence supporting the equatorial current sheet scenario.

We introduce RVSNUpy, a new Python package designed to measure spectroscopic redshifts. Based on inverse-variance weighted cross-correlation, RVSNUpy determines the redshifts by comparing observed spectra with various rest-frame template spectra. We test the performance of RVSNUpy based on ~ 6000 objects in the HectoMAP redshift survey observed with both SDSS and MMT/Hectospec. We demonstrate that a slight redshift offset (~ 40 km/s) between SDSS and MMT/Hectospec measurements reported from previous studies results from the small offsets in the redshift template spectra used for SDSS and Hectospec reductions. We construct the universal set of template spectra, including empirical SDSS template spectra, carefully calibrated to the rest frame. Our test for the HectoMAP objects with duplicated observations shows that RVSNUpy with the universal template spectra yields the homogeneous redshift from the spectra obtained with different spectrographs. We highlight that RVSNUpy is a powerful redshift measurement tool for current and future large-scale spectroscopy surveys, including A-SPEC, DESI, 4MOST, and Subaru/PFS.

NGC2261 is a reflection nebula illuminated by the young star R Monocerotis. Objects moving near the star occasionally cast shadows on the nebula, giving rise to its alternative name: Hubble's Variable Nebula. For 7 years since Spring 2017 robotic telescopes have been used to compile a roughly twice-weekly record of changes in the object. The results, over 1000 images at separate epochs, have been compiled into a movie. This shows that, as well as the large scale but infrequent variability for which it is famous, the nebula is continually traversed by low level `ripples' of light and dark. These record changes in the light output from R Mon and analysis of their progress indicates that the reflecting material takes the form of a thin (<3x10^16cm) screen whose shape resembles a half paraboloid, rooted at the star and bowed towards us. The brightness of the screen in Herschel far-IR maps indicates a density n_H> 1.7x10^5cm^-3 and CO observations show the material is moving towards us at a few km/s relative to the rest cloud, consistent it with being a dense shell of material displaced by R~Mon's outflow. The results demonstrate the value of studying such objects in the time domain, and are a glimpse of what will be achieved by instruments like the Zwicky Transient Facility and Vera Rubin Observatory.

This study involves the use of integral field spectroscopy (IFS) data from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) to investigate whether the morphology influences the environmental dependence of galaxies' colours and the colour-magnitude planes. The galaxies are classified into six morphologies (Elliptical, Lenticular, Early-type, Intermediate-type, Late-type spirals and Irregular) and further in field and group environments. The distributions of colours ($B-V, B-R, u-g, g-r, r-i$ and $i-z$) are compared between field and group environments and then the colour-magnitude planes are analysed. It is observed that Intermediate and Late-types spirals preferentially exist in field environments while Early-type spirals exist in groups. The colours and colour-magnitude planes of Elliptical, Lenticular, Early-type and Intermediate-type spirals depend on the environment while for the Late-type and Irregular galaxies, their dependence on the environment is very weak. The study concludes that the dependence of colours and colour-magnitude planes on the environment is influenced by morphology.

Blazars, supermassive black hole systems (SMBHs) with highly relativistic jets aligned with the line of sight, are the most powerful long-lived emitters of electromagnetic emission in the Universe. We report here on a radio to gamma-ray multiwavelength campaign on the blazar BL Lacertae with unprecedented polarimetric coverage from radio to X-ray wavelengths. The observations caught an extraordinary event on 2023 November 10-18, when the degree of linear polarization of optical synchrotron radiation reached a record value of 47.5%. In stark contrast, the Imaging X-ray Polarimetry Explorer (IXPE) found that the X-ray (Compton scattering or hadron-induced) emission was polarized at less than 7.4% (3sigma confidence level). We argue here that this observational result rules out a hadronic origin of the high energy emission, and strongly favors a leptonic (Compton scattering) origin, thereby breaking the degeneracy between hadronic and leptonic emission models for BL Lacertae and demonstrating the power of multiwavelength polarimetry to address this question. Furthermore, the multiwavelength flux and polarization variability, featuring an extremely prominent rise and decay of the optical polarization degree, is interpreted for the first time by the relaxation of a magnetic "spring" embedded in the newly injected plasma. This suggests that the plasma jet can maintain a predominant toroidal magnetic field component parsecs away from the central engine.

Cosmic rays and cosmic ray induced photons are vital components of chemical evolution in areas of interstellar medium that are impenetrable by external ultraviolet radiation. However, rates of reactions with cosmic ray induced photons used in astrochemical models were calculated for molecular clouds and can be different in protoplanetary discs, where dust grows up to larger sizes. Using ANDES astrochemical model, we study how an increase in both upper dust size and rates of reactions with cosmic ray induced photons can influence species abundances in protoplanetary discs. We show that the increase in these reactions' rates has a significant impact on the ice mass fraction in area between 2 and 20 au but has little impact on ionisation degree in disc.

NGC 6528 and NGC 6553 are among the most metal-rich globular clusters in the Galactic bulge. They represent the upper end of the chemical enrichment in the Galaxy, and can inform on the processes of cluster formation and enrichment. We aim to refine the fundamental parameters of NGC 6528 and NGC 6553, based on proper motion-corrected Hubble Space Telescope WFC3 and ACS photometries. Methods. In order to derive the fundamental parameters age, distance, reddening, and the total-to-selective absorption coefficient, we employed a Bayesian isochrone fitting. Age and metallicity are mainly constrained by the turn-off morphology, thanks to the unprecedented quality of the proper-motion cleaned photometry. The two clusters show remarkably similar Colour-Magnitude Diagrams. We derived an age of 11+-0.5 Gyr with a solar metallicity for both clusters. The reddening for NGC 6528 and NGC 6553 is E(B-V) = 0.63 and 0.76 and the distances from the Sun are d = 7.85 and 5.1 kpc, respectively, recalling that distances strictly depend on the adopted total-to-selective absorption parameter. The age of these metal-rich clusters is about 2 Gyr younger than the moderately metal-poor bulge clusters. The ages and metallicities are remarkably identical to that of the bulk of bulge field stars.

Individual stars in the Milky Way (MW) and its satellites have been shown to trace galaxy stellar mass dependent sequences in the $\alpha$-abundance ([$\alpha$/Fe]) vs metallicity ([Fe/H]) plane. Testing the universality of such sequences has been elusive as deep absorption-line spectra required for [$\alpha$/Fe] and [Fe/H] measurements beyond the local group are mostly limited to integrated light from nearby, relatively high-mass, early-type galaxies. However, analogous to [$\alpha$/Fe] vs [Fe/H] for stars, we now have log(O/Ar) vs 12+log(Ar/H) for the integrated nebular light of star-forming galaxies (SFGs). From Sloan-Digital Sky Survey (SDSS) observations of $\sim3000$ SFGs out to z$\sim0.3$, where we directly determined O & Ar abundances, we obtain for the first time the distribution of an ensemble of SFGs in the log(O/Ar) vs 12+log(Ar/H) plane. We show that higher (<M$\rm_{*}$>$\sim2.6\times10^9$M$_{\odot}$) and lower mass (<M$\rm_{*}$>$\sim1.7\times10^7$M$_{\odot}$) SFGs clearly trace distinct mass dependent sequences in this plane, qualitatively consistent with the mass dependence of chemical enrichment sequences observed for the stars in the MW and its satellites. Such sequences are consistent with expectations from galaxy chemical evolution (GCE) models that are driven primarily by the interplay of core-collapse and Type Ia supernovae.

Simultaneous optical and ionosonde detections of meteors offer a great opportunity to measure the transient physical properties of the meteor's ionization trail. One of the key parameters of the ionization trail is its true geometric height above the surface of the Earth, which can be determined from the trajectory of the optical meteors by multi-site measurements. During the peak of the 2019 Geminid shower we looked for optical meteors among the Konkoly Meteor Observatory Network data taken at ELTE Gothard Observatory, Szombathely, and contemporaneous meteor signals detected by the DPS-4D type ionosonde (Digisonde) operating at Széchenyi István Geophysical Observatory (SZIGO), Nagycenk, Hungary. From the 2 simultaneous detections we inferred the atmospheric height of the meteors using the method of intersecting planes. The results, about 90 and 85 km, are consistent with the detection limit (above 80 km) of the Digisonde.

We present the first results from long-term high-cadence spectroscopic monitoring of 15 PG quasars with relatively strong Fe II emission as a part of a broader reverberation mapping campaign performed with the Calar Alto Observatory 2.2m telescope. The $V$-band, 5100 Å continuum, and H$\beta$ broad emission line light curves were measured for a set of quasars for between dozens to more than a hundred epochs from May 2017 to July 2020. Accurate time lags between the variations of the H$\beta$ broad line fluxes and the optical continuum strength are obtained for all 15 quasars, ranging from $17.0_{-3.2}^{+2.5}$ to $95.9_{-23.9}^{+7.1}$ days in the rest frame. The virial masses of the central supermassive black holes are derived for all 15 quasars, ranging between $0.50_{-0.19}^{+0.18}$ and $19.17_{-2.73}^{+2.98}$ in units of $10^7 M_\odot$. For 11 of the objects in our sample, this is the first reverberation analysis published. Of the rest, two objects have been the subject of previous reverberation studies, but we determine time lags for these that are only half as long as found in the earlier investigations, which had only been able to sample much more sparsely. The remaining two objects have previously been monitored with high sampling rates. Our results here are consistent with the earlier findings in the sense that the time lag and the line width vary inversely consistent with virialization.

We report the results of a long-duration high-cadence reverberation mapping campaign of a second batch of 11 PG quasars using the 2.2m telescope at the Calar Alto Observatory. This follows a similar earlier study of another sample of 15 objects reported by Hu et al. (2021). Among the 11 PG quasars, 8 objects have the H$\beta$ time lags measured for the first time, while the other 3 objects were observed in previous campaigns, but only had highly uncertain H$\beta$-lag measurements. Long-term light curves are presented of photometric $V$-band, spectroscopic 5100 Å continuum, and the H$\beta$ emission line, lasting for $\sim$3--6 years with a cadence of $\sim$6--14 days. Accurate H$\beta$ time lags ranging from $\sim$20 to 150 days in the rest frame are obtained. The estimated virial masses of the central supermassive black holes range from $\sim$(3--300)$\times10^7 M_\odot$. Combining these results with those reported in Hu et al. (2021), we now have 26 PG quasars, with representative properties, having reliable H$\beta$ time-lag measurements from our long-duration high-cadence campaign. A tentative fit to the relation between the H$\beta$ time lag and the continuum luminosity for these 26 objects gives a slope of 0.53.

The transformation of radiation signals (e.g., photon occupation number and integrated intensity) between moving frames is a common task is physics, astrophysics and cosmology. Here we show that the required boost operator, relating the frequency-dependent spin-weighted spherical harmonic coefficients of the considered observable between the frames, is directly given by the aberration kernel with the Doppler weight parameter being replaced by a differential operator. The aberration kernel has been previously studied in great detail, meaning that this simplification allows us to directly compute the boost operator using the expressions of the aberration kernel. As a preparatory step, we generalize the differential equation that determines the aberration kernel to general Doppler weight. This avoids the intermediate step of Doppler-weight raising and lowering operations in computations of the boost operator. We then clarify all the properties of the boost operator (e.g., raising and lowering operations, symmetries and commutation relations) and derive a formal operator differential equation for the boost operator. This differential equation allows us to quickly generate the boost operator for which we give exact expressions up to second order in v/c. For illustration, we then apply the boost operator to transformations of the cosmic microwave background (CMB), validating that measurements of the lowest CMB multipoles do not allow determining the amplitude of the primordial CMB dipole. We also derive the kinematic corrections to the Thomson scattering process (to all orders in v/c), giving explicit expressions up to second order in v/c, showcasing an application of the boost operator in radiative transfer problems.

X-ray transients (XTs) driven by newborn magnetars from mergers of neutron star binaries (NSBs) were occasionally detected in the narrow-field {\it Chandra} Deep Field-South survey (CDF-S) and the {\it Swift}/XRT observations of short gamma-ray bursts (sGRBs). Quantifying their event rate density (ERD) and luminosity function (LF) is critical for understanding NSB coalescence and magnetar formation. Utilizing population synthesis calculations incorporating various equations of state (EoS), we derive a local ERD of $\sim 300\,{\rm Gpc^{-3}\,yr^{-1}}$ and a redshift-dependent ERD profile peaking at $z=1.81$ followed by rapid decline beyond $z \sim 4$. Constructing an XT sample based on CDF-S and {\it Swift} observations, we characterize the LF by a single power-law function at $L \leq 4.75 \times 10^{46}\;{\rm erg\;s^{-1}}$ with a slope of $-1.03$, following by a broken power-law function in which the break luminosity is $L_{\rm b} = 4.38 \times 10^{47}\;{\rm erg\;s^{-1}}$ and the slopes are $-0.28$ and $-1.66$. Based on the ERD and the LF, we estimate that the {\it Einstein Probe} ({\it EP}) detection rate is $\sim 31\;{\rm yr^{-1}}$, adopting a conservative threshold flux of $10^{-9}\;{\rm erg}\;{\rm s^{-1}}$, an luminosity range of $L \in [2\times 10^{44},2\times 10^{49}]\;{\rm erg\;s^{-1}}$, and a correction for jet opening angle of $\sim 16^{\circ}$. This detection rate is consistent with the {\it EP} observations during its first-year operation. It is important to note that our estimation is subject to uncertainties arising from the LF derivation. Future {\it EP} observations of these XT events will be crucial in reducing these uncertainties.

The new method for aerosol measurement using wide-field stellar photometry was originally developed for B-filter data. The dependence of VAOD (vertical aerosol optical depth) on wavelength can be used to understand the physical characteristics of aerosol. The method has been extended for V and R filters using the Gaia stellar photometric catalog as the reference. We calculated the coefficients of the equation describing molecular and aerosol extinction of the atmosphere. For data in a particular filter, they depend on the telescope's location. Light towards shorter wavelengths has higher extinction. The VAOD, calculated in B, V, and R filters does not show monotonous dependence on wavelength, and there is increased dependence on the photometric methods for V and R filters. We studied the possible reasons for this unexpected behavior.

The Chinese Space Station Telescope (CSST) is China's upcoming next-generation ultraviolet and optical survey telescope, with imaging resolution capabilities comparable to the Hubble Space Telescope (HST). In this study, we utilized a comprehensive sample of 3,679 CSST realistic mock galaxies constructed from HST CANDELS/GOODS-North deep imaging observations, with stellar masses $\log\left(M_{*} / M_{\odot}\right) > 9.0$ and redshifts $z < 2$. We evaluate the detection capabilities of CSST surveys and the accuracy in measuring the non-parametric morphological indicators ($C$, $A$, $Gini$, $M_{\rm 20}$, $A_{\rm O}$, $D_{\rm O}$) of galaxies. Our findings show that in terms of galaxy detection capabilities, CSST's deep field surveys can achieve the same level as HST's deep field observations; however, in wide-field surveys, CSST exhibits a significant deficiency in detecting high-redshift, low-mass, low-surface-brightness galaxies. Regarding the measurement of galaxy morphology, CSST's deep field surveys achieve high accuracy across all indicators except for the asymmetry indicator ($A$), whereas its wide-field surveys suffer from significant systematic biases. We thus provide simple correction functions to adjust the non-parametric morphological indicators obtained from CSST's wide-field and deep-field observations, thereby aligning CSST measurements with those from HST. This adjustment enables the direct application of non-parametric morphological classification methods originally developed for HST data to galaxies observed by CSST.

We generalize the magnetically enhanced radiative torque (MRAT) alignment theory for general astrophysical environments described by a dimensionless parameter $U/(n_{1}T_{2})$ with $U$ local radiation strength, $n_{1}=n_{\rm H}/(10{\rm cm}^{-3})$ the hydrogen density, and $T_{2}=T_{\rm gas}/100{\rm K}$ the gas temperature. We first derive the critical magnetic relaxation $\delta_{\rm mag,cri}$ required to produce high-J attractors for different RAT models and local conditions and find that $\delta_{\rm mag,cri}$ must be larger for stronger radiation fields. We then numerically study the grain alignment and rotational disruption by the MRAT mechanism taking into account gas collisions and magnetic fluctuations. We find that, for the collision-dominated (CD) regime ($U/(n_{1}T_{2})\leq 1$), collisional and magnetic excitations can slowly transport large grains from low-J rotation to high-J attractors, leading to the perfect {\it slow alignment} within $\sim 10-100$ damping times due to MRATs. However, for the radiation-dominated (RD) regime ($U/(n_{1}T_{2})>1$), only a fraction of grains can be fast aligned at high-J attractors by MRATs, and the majority of grains are trapped at low-J rotation due to strong radiative torques, a new effect we term {\it radiative torque (RAT) trapping}. For extreme radiation fields of $U/(n_{1}T_{2})>10^{4}$, the efficiency of magnetic relaxation on grain alignment is suppressed, and grains only have fast alignment and disruption purely determined by RATs. We quantified the fraction of grains with fast alignment at high-J attractors, $f_{\rm high-J}^{\rm fast}$, for different RAT models, magnetic relaxation, and $U/(n_{1}T_{2})$, and found that the maximum $f_{\rm high-J}^{\rm fast}$ can reach $45\%$ by MRATs and $22\%$ by RATs.

Recent DESI DR2 baryon acoustic oscillation (BAO) measurement shows inconsistency with the cosmic microwave background (CMB) observation in $\Lambda$CDM. We used the Bayesian suspiciousness and goodness-of-fit tension metrics to quantify the significance of the tension. Both tension metrics consistently confirm a $\sim2\sigma$ tension between DESI BAO and CMB. We identified the key assumptions of $\Lambda$CDM underlying the tension between BAO and CMB. If the tension is physical, one or multiple of the assumptions might need to be broken by the potential new physics addressing the tension. In particular, the analysis shows that if dynamical dark energy is the only new physics behind the tension, phantom crossing is required to restore concordance between CMB and BAO.

We explore the impact of finite-temperature quantum gravity effects on cosmological parameters, particularly the cosmological constant $\Lambda$, by incorporating temperature-dependent quantum corrections into the Hubble parameter. For that purpose, we modify the Cosmic Linear Anisotropy Solving System. We introduce new density parameters, $\Omega_{\Lambda_2}$ and $\Omega_{\Lambda_3}$, arising from finite-temperature quantum gravity contributions, and analyze their influence on the cosmic microwave background power spectrum using advanced machine learning techniques, including artificial neural networks and stochastic optimization. Our results reveal that $\Omega_{\Lambda_3}$ assumes a negative value, consistent with dimensional regularization in renormalization and that the presence of $\Omega_{\Lambda_2}$ as well as $\Omega_{\Lambda_3}$ significantly enhances model accuracy. Numerical analyses demonstrate that the inclusion of these parameters improves the fit to 2018 Planck data, suggesting that finite-temperature quantum gravity effects play a non-negligible role in cosmological evolution. Although the Hubble tension persists, our findings highlight the potential of quantum gravitational corrections in refining cosmological models and motivate further investigation into higher-order thermal effects and polarization data constraints.

The mass range of observed black holes extends from stellar-mass to supermassive scales, yet the existence of objects in the intermediate-mass range of $10^{2} - 10^{5} \text{M}_{\odot}$ remains unconfirmed. Black holes are suspected to compress the surrounding dark matter distribution, forming a ``spike''. If dark matter is self-annihilating, the spike could produce gamma-ray emission sufficiently luminous to be detected. This work aims to estimate the number of expected unmerged intermediate-mass black holes in a Milky Way-like galaxy that could form such spikes. These intermediate-mass black holes are assumed to have formed from the collapse of high-mass Population III stars, such that the resulting merger rate is constrained by observations of gravitational wave emission. It is furthermore estimated to what extent the progenitor Population III stars contribute to the extragalactic background light. The Population III stars are simulated and tracked using the A-SLOTH semi-analytical simulation code and the resulting number of intermediate-mass black holes is constrained by applying the Population III binary black hole merger rate to an effective volume determined from the Population III star formation rate. In this framework, $\sim 130$ unmerged IMBHs from Population III stars are expected to reside in a Milky Way-like galaxy. The contribution of their progenitors to the EBL in the near-infrared is less than $10^{-3} \text{nW} \text{m}^{-2} \text{sr}^{-1}$, well below previous estimates.

We utilized the Edelson and Malkan (2012) and Stern et al. (2012) selection techniques and other methods to identify AGN candidates that were monitored during the Kepler prime and K2 missions. Subsequent to those observations, we obtained 125 long-slit optical spectra with the Lick 3-m telescope, 58 spectra with the Palomar 5-m telescope, and three with the Keck 10-m telescope to test these identifications. Of these 186 AGN candidates, 105 were confirmed as Type 1 AGN and 35 as Type 2 AGN, while the remaining 46 were found to have other identifications (e.g., stars and normal galaxies). This indicated an overall reliability of about 75%, while the two main methods had much higher reliability, 87%-96%. The spectra indicated redshifts out to z = 3.4. Then, we examined the AGN sample properties through the Baldwin, Phillips & Terlevich diagram and compared the AGN's spectral energy distributions (SEDs) with those from the literature. We found that our sample yielded the same AGN population as those identified through other methods, such as optical spectroscopy.

Massive stars are extremely influential objects and multiplicity has the potential to dictate how they evolve over time by causing dynamical interactions, common-envelope evolution, mergers and more. While O-star multiplicity has been studied over a broad separation range, studies of B star multiplicity are lacking despite the fact that they dominate the production of core-collapse supernovae and neutron stars. We analyse interferometric data at high-angular resolution taken with PIONIER/VLTI for a sample of 32 B stars. Using parametric modelling of the closure phases and visibilities, we determine best-fit models to each of the systems and investigate whether each source is best represented by a single star or a higher-order system. Detection limits were calculated for companions to determine if they were significant. We then combined our findings from the interferometric data with results from a literature search to see if other companions were reported at different separation ranges. We find that, within the interferometric range, 72$\pm$8\% of the B stars are resolved as multiple systems. The most common type of system is a triple system, followed by binary systems, presumably single stars and then higher-order systems. The interferometric companion fraction derived for the sample is 1.88$\pm$0.24. When accounting for the spectroscopic companions confirmed in the literature and wide companions inferred from Gaia data in addition to companions we find with interferometry, multiplicity and companion fractions of 0.88$\pm$0.06 and 2.31$\pm$0.27 respectively are obtained for our sample. The number of triple systems in particular increases significantly when accounting for spectroscopic companions, suggesting that binarity and higher-order multiplicity are integral to the evolution of B stars, as they are for O stars.

The abundance of dust within galaxies directly influences their evolution. Contemporary models attempt to match this abundance by simulating the processes of dust creation, growth, and destruction. While these models are accurate, they require refinement, especially at earlier epochs. This study aims to compare simulated and observed datasets and identify discrepancies between the two, providing a basis for future improvements. We utilise simulation data from the SIMBA cosmological simulation suite and observed data from the Galaxy and Mass Assembly (GAMA), a subset of the Cosmic Evolution Survey (G10-COSMOS), and the Hubble Space Telescope (3D-HST). We selected galaxies in the observed and simulated data in a stellar mass range of $(10^{8.59} < M_\odot < 10^{11.5})$ and at redshift bins centering around $(z = 0.0)$, $(z = 0.1)$, $(z = 0.5)$, $(z = 1.0)$, and $(z = 1.5)$ in a homogeneous dust mass range $((10^{6} < M_D [M_\odot] < 10^{9}))$. Our results show notable deviations between SIMBA and observed data for dust-poor and rich galaxies, with strong indications that differences in galaxy populations and SIMBA limitations are the underlying cause rather than the dust physics implemented in SIMBA itself.

We aim to detect activity cycles in young main-sequence stars, analogous to the 11-year solar cycle, using combined photometric survey data. This research will enhance our understanding of how cycle periods relate to rotation rates in fast-rotating stars. We measured activity cycles for 138 G-K-type main sequence stars using combined time-series photometry spanning ~14 years. The first set of 70 stars used data from Kepler Full Frame Images (FFIs)-ASAS-SN-ZTF, and the second set of 68 stars used data from Kepler-FFIs-ZTF. Additionally, we measured the activity cycles for 25 RS CVn candidates. For our sample, we analyzed the correlation or anti-correlation between flux variations and photospheric activity, which arises due to the presence of faculae or starspots. We identified fast-rotating K-type stars that are faculae-dominated by tracking spot or faculae evolution in Kepler RMS data. Our findings reveal that fast-rotating G-K-type stars show no strong correlation between cycle length and rotation period. Previous studies have identified active and inactive branches in the cycle-rotation diagram. However, we find that G-K-type stars do not show a clear trend aligning with the active branch, with 34 per cent of stars falling within the intermediate region between the two branches, where our Sun resides. Our results highlight that the proposed distinction between the two branches may not be as definitive as previously thought, particularly regarding the active branch. Furthermore, we also detected 23 per cent of young Sun-like stars in the intermediate region, where our Sun is located, implying that our Sun may not be unique in this regard.

The Higgs inflation model with nonminimal coupling, while disfavored by the 1$\sigma$ region of the latest Atacama Cosmology Telescope (ACT) observational data, can be reconciled with the ACT data by incorporating the effects of reheating. In this paper, we consider reheating with a constant equation of state $w_{re}$. To simultaneously satisfy the ACT data and ensure that the temperature at the end of reheating is above the threshold required for Big Bang Nucleosynthesis, we find that the equation of state must satisfy $w_{re} \gtrsim 2/3$. For the special cases of $w_{ re} = 2/3$ and $w_{ re} = 1$, the number of e-folds during reheating must lie within the ranges $31.6 \leq N_{ re} \leq 32.8$ and $15.8 \leq N_{ re} \leq 27.3$, respectively. Our findings suggest that by considering reheating, a wide range of inflationary models, such as $R^2$ inflation, hilltop inflation, E-model inflation, and T-model inflation, can also be made consistent with the ACT observational data.

The formation of complex organic molecules in the interstellar medium (ISM) is central to astrochemistry and prebiotic chemistry, as these species may act as precursors to biomolecules essential for life. Although glyceraldehyde (HOCH2CH(OH)C(O)H, GCA) has not yet been detected in the ISM, the presence of structurally related compounds in various astronomical environments suggests that it may form under interstellar conditions. In this study, we employed the automated reaction discovery tool AutoMeKin to systematically explore the gas-phase chemical reaction networks (CRNs) of C3H6O3 (GCA), C3H7O3 (a hydrogenated analog), and C2H5O2. Reaction pathways were characterized at the wB97XD/Def2-TZVPP level of theory, and rate coefficients for key processes were computed using the competitive canonical unified statistical (CCUS) model, which accounts for multiple dynamic bottlenecks. Our analysis revealed several barrierless pathways leading to GCA or to GCA and a leaving group. Notably, the reaction between glyoxal (HCOHCO) and the HOCHCH2OH radical, though neither has yet been detected in the ISM, was found to efficiently produce GCA and a formyl radical. However, aside from the aforementined exception, most GCA formation channels result in highly vibrationally excited intermediates that are more likely to undergo rapid unimolecular decomposition than to be stabilized by radiative emission under typical ISM conditions. These results suggest that while gas-phase GCA formation is chemically feasible, it is likely transient and difficult to detect directly. In contrast, alternative products such as formaldehyde, glycolaldehyde, and (Z)-ethene-1,2-diol dominate many pathways and align better with current astronomical observations.

Tidal disruption events (TDEs) are luminous black hole (BH) transient sources, which are detected mainly in X-ray and optical bands. It is generally believed that the X-ray emission in TDEs is produced by an accretion disc formed as the stellar debris accreted onto the central BH. The origin of the optical emission is not determined, but could be explained by the `reprocessing' model with the X-ray emission reprocessed into optical band by a surrounding optically thick envelope or outflow. In this paper, we performed radiation hydrodynamic simulations of super-Eddington accretion flow with Athena++ code in the environment of TDEs, i.e., injecting a continuous mass flow rate at the circularization radius in the form of $\dot M_{\rm inject} \propto t^{-5/3}$ for the mass supply rate. We show that a significant fraction of the matter in the accretion inflow are blowed off forming outflow, and the properties of the outflow are viewing-angle dependent. We further calculate the emergent spectra of such an inflow/outflow system for different viewing angles with the method of Monto Carlo radiative transfer. Based on the emergent spectra, we show that the observed features of TDEs, such as the X-ray and optical luminosities, the blackbody temperature of X-ray and optical emission and the corresponding emission radii, the ratio of X-ray luminosity to optical luminosity, as well as the evolution of these quantities can be explained in the framework of viewing-angle effect of super-Eddington accretion around a BH.

We investigate the impact of early-time initial conditions on nonlinear structure formation and evolution within the framework of the semi-analytical Excursion Set Theory (EST). Our analysis reveals that adding a Gaussian bump to the initial curvature power spectrum at small scales enhances the abundance of massive halos while sharply reducing the number of small-mass halos, and consequently, satellite galaxies. Moreover, this modification increases the frequency of major mergers while suppressing high-mass-ratio minor mergers. These features may offer resolutions to the missing satellite and Too Big to Fail (TBTF) problems. In underdense regions -- voids -- the same modifications increase the likelihood of finding massive halos embedded in voids while similarly decreasing the small-halo population, and consequently, faint galaxies. This behavior suggests a potential solution to the void phenomenon, in which embedded halos, despite being too massive, were too rare to be noticed. More precisely, our results indicate that an excess of massive structures emerges at mass scales near the center of the Gaussian bump: $k_* = 1.85 \,\rm{h/Mpc}$ and $k_* = 3.95 \,\rm{h/Mpc}$. These scales correspond to mass scales of $M_* = 10^{11}$ and $M_* = 10^{10}$, respectively. This modification extends up to two orders of magnitude in higher mass scales, while reducing the abundance of halos below $M_*$ by two to three orders of magnitude. Additionally, we find that evolutionary conditions, halo-in-halo, and particularly halo-in-void statistics serve as more sensitive and complementary probes for differentiating among cosmological models.

Supernova remnants (SNRs) are often considered as the main sites of acceleration of cosmic rays in our Galaxy, possibly up to the knee. However, their ability to accelerate particles to reach PeV energies is questionable and lacks observational evidence. Theoretical predictions suggest that only a small subclass of very young SNRs evolving in dense environments could potentially satisfy the necessary conditions to accelerate particles to PeV energies. Most such theoretical investigations are carried out either in the standard interstellar medium or in the wind of the progenitor. Since most core collapse supernovae occur in star clusters, it is important to extend such investigation to SNRs taking place in a star cluster. In this work we focus on a SNR shock propagating in the collective wind of a compact star cluster, and we study the acceleration process as a function time, with special emphasis on the maximum energy of accelerated particles. Using both analytic and numerical approaches we investigate the spectrum of accelerated particles and maximum achievable energy in the case of pre-existing turbulence in the collective wind and self-generated magnetic perturbations. We find that similar to isolated SNRs, acceleration to PeV energies is plausible only for extreme conditions achievable only in a small subset of SNRs.

The goal of the MeerKAT radio telescope's pulsar timing array programme (MPTA) is the detection of gravitational waves (GWs) of nanohertz frequencies. Evidence for such a signal was recently announced by the MPTA and several other pulsar timing array (PTA) consortia. Given an array of pulsars and an observation strategy, we consider whether small adjustments to the observing schedule can provide gains in signal-to-noise ratio (S/N) for a stochastic GW background signal produced by a population of massive black hole binaries. Our approach uses a greedy algorithm to reallocate available integration time between pulsars in the array. The overall time dedicated to MPTA observing is kept constant so that there is only minimal disruption to the current observation strategy. We assume a GW signal consistent with those reported. For the sake of demonstrating our method, we also make several simplifying assumptions about the noise properties of the pulsars in the MPTA. Given these assumptions, we find that small adjustments to the observing schedule can provide an increased S/N by $\approx 20\%$ for a $10\,{\rm yr}$ PTA lifespan.

Observations have shown the presence of the Kelvin-Helmholtz instability (KHi) in solar prominences. Effects due to partial ionization of the prominence plasma may influence the KHi onset. We study the triggering of the KHi in an interface model that consists of a partially ionized prominence region and a fully ionized coronal region, with a uniform magnetic field parallel to the interface. There is a longitudinal flow in the prominence region. The plasma is compressible and the role of ambipolar diffusion, which accounts for collisions between charges and neutrals, is taken into account in the prominence plasma. We derive the dispersion relation of linear perturbations on the interface and analyze some limit cases analytically. Numerical results are obtained for realistic prominence parameters. We find that compressibility and gas pressure are important in determining the unstable flow velocities, specially in the range of sub-Alfvénic flows that are consistent with the observations. The ambipolar diffusion has a generally destabilizing influence and reduces the threshold flow velocity for the KHi onset.

Water vapour is delivered to Saturn's stratosphere by Enceladus' plumes and subsequent diffusion in the planet system. It is expected to condense into a haze in the middle stratosphere. The hot stratospheric vortex (the `beacon') that formed as an aftermath of Saturn's Great Storm of 2010 significantly altered the temperature, composition, and circulation in Saturn's northern stratosphere. Previous photochemical models suggested haze sublimation and vertical winds as processes likely to increase the water vapour column density in the beacon. We aim to quantify the temporal evolution of stratospheric water vapour in the beacon during the storm. We mapped Saturn at 66.44 and 67.09 $\mu$m on seven occasions from July 2011 to February 2013 with the PACS instrument of the Herschel Space Observatory. The observations probe the millibar levels, at which the water condensation region was altered by the warmer temperatures in the beacon. Using radiative transfer modelling, we tested several empirical and physically based models to constrain the cause of the enhanced water emission found in the beacon. The observations show an increased emission in the beacon that cannot be reproduced only accounting for the warmer temperatures. An additional source of water vapour is thus needed. We find a factor (7.5$\pm$1.6) increase in the water column in the beacon compared to pre-storm conditions using empirical models. Combining our results with a cloud formation model for July 2011, we evaluate the sublimation contribution to 45-85% of the extra column derived from the water emission increase in the beacon. The observations confirm that the storm conditions enhanced the water abundance at the millibar levels because of haze sublimation and vertical winds in the beacon. Future work on the haze temporal evolution during the storm will help to better constrain the sublimation contribution over time.

We study the secular periodic evolution of quasi-periodic eruptions (QPEs) for GSN069 and eRO-QPE2 assuming that they are driven by star-disc collisions. We set up numerical simulations and compared them with the observed periodic decay of $\sim -3160\pm720$ s yr$^{-1}$ in GSN069 and $\sim -370\pm40$ s yr$^{-1}$ in eRO-QPE2. We find that: (1) Stellar mass black holes are unlikely the orbiters in these two sources, as their periodic decay are on the order of $<10$ s yr$^{-1}$; (2) A naked degenerate core (including white dwarf) is unlikely the orbiter in GSN069, as the decay is on the order of $<200$ s yr$^{-1}$. However, it is possible in eRO-QPE2, although the required surface density of the accretion disc is relatively high (e.g., $\Sigma\gtrsim10^7\sim 10^8$ g cm$^{-2}$); (3) Both the orbiters in GSN069 and eRO-QPE2 can be solar-like main-sequence stars (MSs). However, each collision can lead to gradual ablation of the stellar envelope in the order of $10^{-5}\sim 10^{-3}M_\odot$. To reproduce the observed decay while surviving for $\gtrsim 3$ yr, the surface density of the disc needs to be within a certain range. For example, given a $1M_\odot$ MS orbiter the surface density of the disc gas should be in the range of $3\times10^5\sim 2\times10^6$g cm$^{-2}$ for GSN069 or $5\times10^4\sim 10^6$ g cm$^{-2}$ for eRO-QPE2. In both of these two sources, the MS can not survive for more than $\sim 12$ yr. We expect that future observations of these two sources can help to distinguish whether the orbiters are degenerated compact objects or gaseous stars.

We describe a new method to incorporate thermonuclear heating in the envelope of accreting neutron star into long term simulations of their thermal evolution. We obtain boundary conditions for the heat exchange between the envelope and the crust based on stationary models which include nuclear burning and validate these values comparing to the results of the time-dependent code \texttt{MESA}. These simple boundary conditions allow us to explore a large parameter space. We quantify the amount of heat flowing from the envelope into the crust, or viceversa, depending on the mass accretion rate, outburst duration and duty cycle, and especially crust/core physical parameters such as impurities, crustal heating, and neutrino cooling rate.

Supernova observations imply the presence of a dense and asymmetric circumstellar environment around SN Type II progenitors, whereas the mass loss from these progenitors, namely, red supergiants, is still poorly constrained. We aim to characterise the dust and gas in the circumstellar environment of the extreme Galactic red supergiant \nmlcyg in terms of mass, morphology, and kinematics. Based on interferometric observations with NOEMA at 230 GHz we estimated dust masses and temperatures, and measured the extent and morphological complexity of the circumstellar environment. We detected two strong continuum components, amounting to an estimated total dust mass of $\sim2\times10^{-3}M_{\odot}$ located out to ~2000 AU from the star, largely beyond the dust detected at optical/infrared wavelengths. The extent of the detected CO emission supports the notion that the outflow is formed by a mass-loss rate of several $10^{-4}\,M_{\odot}$/year and that it is not primarily shaped by extreme irradiation from the Cyg OB2 cluster it has been associated with. We have detected, but not resolved, previously unseen high-velocity components close to the star. The observations reveal a very complex circumstellar morphology and we propose that some of the detected components could be the imprint of a hitherto unknown binary companion.

We show that the latest empirical constraints on cosmology from a combination of DESI, CMB and supernova data, can be accounted for if a small component of dark matter has an evolving and oscillating equation of state within $-1<w<1$. From a fundamental physics perspective, this interpretation is more appealing than an evolving phantom dark energy with $w<-1$, which violates the null energy condition.

This paper presents an alternative way of analysing Baryon Acoustic Oscillation (BAO) distance measurements via rotations to define new quantities Dperp and Dpar. These quantities allow simple tests of consistency with the Planck LCDM cosmology. The parameter Dperp is determined with negligible uncertainty from Planck under the assumption of LCDM. Comparing with measurements from the Dark Energy Spectroscopic Instrument (DESI), we find that the measurements of Dperp from Data Release 2 (DR2) move into significantly better agreement with the Planck LCDM cosmology compared to DESI Data Release 1 (DR1). The quantity in the orthogonal direction Dpar provides a measure of the physical matter density omega_m in the LCDM cosmology. The DR2 measurements of Dpar\ also come into better agreement with Planck LCDM compared to the earlier DR1 results. From the comparison of Planck and DESI BAO measurements, we find no significant evidence in support of evolving dark energy. We also investigate a rotation in the theory space of the w_0 and w_a parameterization of the dark energy equation-of-state w(z). We show that the combination of DESI BAO measurements and the CMB constrain w(z=0.5) = -0.996 pm 0.046, i.e. very close to the value expected for a cosmological constant. We present a critique of the statistical methodology employed by the DESI collaboration and argue that it gives a misleading impression of the evidence in favour of evolving dark energy. An Appendix shows that the cosmological parameters determined from the Dark Energy Survey 5 Year supernova sample are in tension with those from DESI DR2 and parameters determined by Planck.

We follow the spin vector evolutions of well resolved dark matter haloes (containing more than 300 particles) in merger tree main branches from the Millennium and Millennium-II N-body simulations, from z about 3.3 to z = 0. We find that there seems to be a characteristic plane for the spin vector evolution along each main branch. In the direction perpendicular to it, spin vectors oscillate around the plane, while within the plane, spin vectors show a coherent direction change as well as a diffusion in direction (possibly corresponds to a Gaussian white noise). This plane may reflect the geometry of surrounding large-scale structures. We also construct a simple stochastic model in which halo spin vector evolution is assumed to be driven by accretion of halo mass and angular momentum. This model can reproduce major features of the results from N-body simulations.

We investigate the effect of the environment on the infrared and radio emission of cluster galaxies during the transition epoch at 1 < z < 2 when they first start to quench consistently in the majority of galaxy clusters. We considered a sample of 129 cluster member galaxies from 11 massive clusters at a confirmed redshift of 1.0-1.8 from the IRAC Shallow Cluster Survey (ISCS), the IRAC Distant Cluster Survey (IDCS), and new 3 GHz images from the Karl G. Jansky Very Large Array (VLA). We calculated the IR-radio correlation slope parameter, q, in order to identify differences in the ratios of IR to radio of cluster galaxies and field galaxy comparison samples at different redshifts. Active galactic nuclei (AGNs) were identified and analyzed to search for any effect on the IR-radio correlation. The correlation parameter values were also compared by the Kolmogorov-Smirnov test with field galaxies. Our comparison of the IR-radio correlation in cluster galaxies to the control sample of field galaxies reveals a marginally to modestly significant difference in the correlation slope parameter at the ~2sigma-3sigma level. A split of the clusters into low-redshift (1 < z < 1.37) and high-redshift (1.37 < z < 1.8) bins indicates a more significant difference in the correlation parameter in the lower redshift cluster subsample, where widespread quenching begins. We find no difference in the IR-radio correlation between galaxies that host AGNs and non-active star-forming galaxies either. This suggests that our AGNs are overwhelmingly radio quiet and therefore do not affect the results we described above. We conclude that further investigations based on larger datasets are needed to constrain the impact of the cluster environment on the IR-radio correlation better.

A plausible explanation for the absence of primordial argon, krypton, and xenon in Titan's current atmosphere is that these gases were sequestered in clathrate hydrates during Titan's "open-ocean" phase. We examine how clathrate hydrate formation at Titan's ocean surface in its early history may have contributed to noble gas depletion in the primordial atmosphere. Starting with vapor-liquid equilibrium modeling between water and volatiles, we used a statistical thermodynamic model to determine the clathrate hydrate crust thickness needed to deplete the primordial atmosphere of noble gases. Our computations suggest that if Titan's volatile budget was delivered by icy planetesimals with a comet-like composition, its primordial atmosphere should be rich in CO$_2$ and CH$_4$, with NH$_3$ largely retained in water as ions. We show that at 273.15 K, a clathrate crust tens of kilometers thick would deplete the primordial atmosphere of xenon and krypton. The lack of primordial argon in Titan's atmosphere may result from the partial de-volatilization of its accreted materials.

Abridged: Supermassive black hole binaries (SMBHBs) separated by (sub)-pc scales represent one of the latest stages of hierarchical galaxy assembly. However, many of these objects are hidden behind large columns of gas and dust at the center of galaxies and are difficult to detect. The X-ray and UV emission in these systems are predicted to vary regularly on timescales comparable to that of the orbital period of the binary. This is the first of a series of papers where we aim to find SMBHB candidates based on quasi-periodic light curves from the soft X-ray instrument eROSITA on board the Spectrum-Roentgen-Gamma (SRG) observatory and X-ray follow-up. We searched the multi-epoch SRG/eROSITA all-sky surveys for extragalactic sources that show an `up-down-up-down' or `down-up-down-up' profile (from scan to scan) in their 0.2--2.3 keV flux light curves. We compiled a sample of 16 sources that are suitable for X-ray follow-up campaigns given their brightness and significant variability between bright and faint SRG/eROSITA flux levels. We triggered extensive Swift-XRT and NICER monitoring campaigns on the best SMBHB candidates to confirm or discard their tentative periodicities. Optical spectroscopic observations confirmed the nuclear and extragalactic nature of 15/16 objects and enabled single-epoch SMBH mass measurements and BPT classifications. Our most promising candidate, eRASSt J0530-4125, shows X-ray quasi-periodic variability with a typical time scale of one year in the observed frame. By stacking the X-ray observations of each source in our sample, we find that 14/15 sources can be modeled by a power law with a photon index ranging from $\Gamma\sim1.8-2.8$. Based on our selection, we estimate an optimistic upper limit on the fraction of SMBHB candidates to be $\sim 0.05$ per galaxy. We emphasize that further observational evidence is needed to confirm the SMBHB nature of our sources.

The reionisation of the intergalactic medium (IGM) was driven by the first stars, galaxies, and accreting black holes. However, the relative importance of these sources and the efficiency by which ionising photons escape into the IGM remain poorly understood. Most reionisation modelling frameworks assume idealised, isotropic emissions. We investigate this assumption by examining a suite of simulations incorporating directed, anisotropic photon emissions. We find that such anisotropic emissions of ionising photons yield a different reionisation geometry compared to the standard, isotropic, case. During the early stages of reionisation (when less than 30 per cent of the Universe is ionised), simulations with narrow photon leakage channels produce smaller ionised bubbles on average. However, these bubbles grow to similar sizes during the middle stages of reionisation. This anisotropy not only produces a distinctive evolution of the size distribution of the ionised regions, but also imprints a feature onto the spherically averaged power spectra of the 21-cm signal throughout reionisation. We observe a suppression in power by about 10-40 per cent at scales corresponding to wavenumbers $k = 0.1-1 \, h \, \mathrm{Mpc}^{-1}$, corresponding to the range in which current radio interferometers are most likely to measure the power spectrum. The simulation with the narrowest channel of ionisation emission shows the strongest suppression. However, this anisotropic emission process does not introduce any measurable anisotropy in the 21-cm signal.

Forced magnetic reconnection is triggered by external perturbations, which are ubiquitous in the solar corona. This process plays a crucial role in the energy release during solar transient events, which are often associated with electric current sheets (CSs). The CSs can often disintegrate through the development of the tearing instability, which may lead to the formation of plasmoids in the non-linear phase of evolution. However, the complexity of the dynamics, and the magnetic and thermodynamic evolution due to the coalescence of the plasmoids are not fully understood. We used a resistive magnetohydrodynamic simulation of a 2.5D current layer embedded in a stratified medium in the solar corona, incorporating field-aligned thermal conduction. Multiple levels of adaptive mesh-refined grids are used to resolve the fine structures that result during the evolution of the system. The instability in the CS is triggered by imposing impulsive velocity perturbations concentrated at three different locations in the upper half along the CS plane; this leads to the formation of plasmoids and their later coalescence. We demonstrate that a transition from purely 2D reconnection to 2D reconnection with guide field takes place at the interface between the plasmoids as the latter evolve from the pre-merger to the merged state. The small-scale, short-lived, and collimated outflows during the merging process share various physical properties with the recently discovered nanojets. The subsequent thermodynamic change within and outside the merged plasmoid region is governed by the combined effect of Ohmic heating, thermal conduction, and expansion/contraction of the plasma. Our results imply that impulsive perturbations in coronal CSs can be the triggering agents for plasmoid coalescence, which leads to the subsequent magnetic, and thermodynamic change in and around the CS.

Carbon plays key roles in the InterStellar Medium (ISM) -- as a constituent of dust, as the carrier of the dominant far-infrared cooling line, and as a component of various important molecules. But despite this, there are very few measurements of the abundance and depletion of carbon in the diffuse ISM. As with other elements, these measurements are traditionally performed in the ultraviolet. But for carbon, such measurements are extremely difficult, and less than 20 have been reported in the literature to date. Here, we present a novel method of measuring the abundance and depletion of carbon in the diffuse ISM: by observing absorption of the 158 $\mu$m [CII] line in the far-infrared. We present a catalog of 432 candidate sightlines that use bright nearby galaxies as background sources, and predict the [CII] absorption expected towards each. We conducted a pilot study using SOFIA, targeting sightlines towards the galaxies IC342 and Circinus. We report a potential detection of Galactic [CII] absorption along the IC342 sightline, although it requires disentangling [CII] emission from IC342 itself. The Circinus sightline had an insufficiently stable instrumental baseline to allow a detection. This SOFIA study informs the prospects for [CII] absorption measurements with future facilities. To that end, we explore the potential for four proposed future FIR telescopes -- PRIMA, FIRSST, SALTUS, and Origins -- to detect [CII] absorption. We find that all four facilities would be able to detect [CII] absorption along a significant number of sightlines.

To investigate the impact of winds and low-to-moderate power jets on the cold molecular gas reservoirs of AGN, we present high angular resolution ALMA CO(2-1) and CO(3-2) observations of a sample of six type-2 quasars (QSO2s) at z$\sim$0.1 from the Quasar Feedback (QSOFEED) sample. Spatially resolved molecular line ratio maps, defined as $R_{32}=L'_{CO(3-2)}/L'_{CO(2-1)}$, and kinematic modelling were used to constrain changes in gas excitation and to identify gas outflows, respectively. We find that the molecular outflows are co-spatial with regions with $R_{32}$>1, indicating enhanced temperature relative to the discs and the presence of optically thin gas in the outflows. We find mass outflow rates of 5$<\dot{M}_{out}<$150$M_\odot$/yr, much lower than those expected from their AGN luminosities of $10^{45.5-46}$erg/s. The outflow kinetic energies might be driven by the combined action of jets and winds/radiation pressure, with radiative coupling efficiencies ($\epsilon_{AGN}=\dot{E}_{out}/L_{bol}$) ranging from $10^{-6}<\epsilon_{AGN}<10^{-4}$ and jet coupling efficiencies ($\epsilon_{jet}=\dot{E}_{out}/P_{jet}$) from $10^{-3}<\epsilon_{AGN}<10^{-2}$. A linear regression including the six QSO2s follows the locus of $\epsilon_{jet}\sim$0.1\%. Our results provide evidence that AGN-driven jets/winds disturb the molecular gas kinematics and excitation within the central kpc of the galaxies. The coupling between compact jets and the ISM might be relevant to AGN feedback, even in the case of radio-quiet galaxies, which are more representative of the AGN population. Finally, we find that the warm and cold molecular gas phases seem to be tracing the same outflow, with the main distinction between them being the mass they carry, while the warm ionized outflows do not seem to be another face of the same outflow, as they show different orientation, velocity, and radius.

We report discoveries of bow-shock nebulae, seen in Halpha and [O III] 5007, around two cataclysmic variables (CVs): LS Pegasi and ASASSN-V J205457.73+515731.9 (hereafter ASASJ2054). Additionally, both stars lie near the edges of faint extended Halpha-emitting nebulae. The orientations of the bow shocks are consistent with the directions of the objects' proper motions. The properties of LS Peg and ASASJ2054, and of their nebulae, are remarkably similar to those of SY Cancri, which we described in a recent paper; SY Cnc is a CV likewise associated with a bow shock and an off-center Halpha nebula. These objects join V341 Arae and BZ Camelopardalis, CVs that are also accompanied by similar nebulae. All five stars belong to the nova-like variable (NLV) subclass of CVs, characterized by luminous optically thick accretion disks that launch fast winds into the surrounding space. We suggest that the bow shocks and nebulae result from chance encounters of the NLVs with interstellar gas clouds, with the stars leaving in their wakes Stromgren zones that are recombining after being photoionized by the CVs' ultraviolet and X-ray radiation. Our discoveries illustrate the power of small telescopes equipped with modern instrumentation, and used to accumulate extremely long exposure times, for the detection of very low-surface-brightness nebulae.

Galaxy clusters are powerful probes of cosmology, and the halo mass function (HMF) serves as a fundamental tool for extracting cosmological information. Previous calibrations of the HMF in dynamical dark energy (DE) models either assumed a homogeneous DE component or a fixed sound speed of unity, which strongly suppresses DE perturbations. We extend the HMF calibration to clustering dark energy (CDE) models by allowing for a sound speed $(c_{\rm s})$ value different than unity. This generalization enables a broader description of the impact of DE perturbations on structure formation. Our approach builds upon the DUCA simulation suite that accounts for DE at the background and perturbative levels. We present an HMF calibration based on introducing an effective peak height while maintaining the multiplicity function as previously calibrated. The effective peak height is written as a function of the peak height computed using the matter power spectrum of the homogeneous DE case, but it is modulated by the amplitude of DE and matter perturbations on the non-homogeneous case at the turnaround. The model depends on one single parameter, which we calibrate using $N$-body simulations, following a Bayesian approach. The resulting HMF model achieves sub-percent accuracy over a wide range of $c_{\rm s}$ values. Our analysis reveals that, although the overall impact of CDE on halo abundances remains modest (typically a few percent), the effects are more pronounced in non-phantom DE scenarios. Our model qualitatively agrees with predictions based on the spherical collapse model, but predicts a significantly lower impact for low $c_{\rm s}$. Our results underscore the need for more precise modeling of CDE's nonlinear regime. Numerical simulations and theoretical approaches must be advanced to capture the complex interplay between DE perturbations and matter fully.

OCARS (Optical Characteristics of Astrometric Radio Sources) is a compiled catalog of various additional data associated with astrometric radio sources whose coordinates have been determined from very long baseline interferometry (VLBI) observations. It contains source coordinates, object type, redshift, optical and near-infrared magnitudes. Until now, OCARS source coordinates were simply copied from input catalogs and, as a result, were systematically inhomogeneous. This work is the first attempt to obtain a unified set of radio source coordinates aligned to the International Celestial Reference Frame (ICRF), more specically to the third ICRF release, ICRF3. Comparison of the source coordinates in the old OCARS version as of December 2024 and the new OCARS version as of March 2025 with the ICRF3-SX catalog using the vector spherical harmonics (VSH) technique showed almost complete elimination of systematic errors in new OCARS positions relative to the ICRF3 frame.

The neutron star zoo comprises several sub-populations that range from energetic magnetars and thermally emitting X-ray neutron stars to radio-emitting pulsars. Despite studies over the last five decades, it has been challenging to obtain a clear physical link between the various populations of neutron stars, vital to constrain their formation and evolutionary pathways. Here we report the detection of a burst of coherent radio emission from a known radio-quiet, thermally emitting neutron star 2XMM J104608.7$-$594306in the Carina Nebula. The burst has a distinctive sharp rise followed by a decay made up of multiple components, which is unlike anything seen from other radio-emitting neutron stars. It suggests an episodic event from the neutron star surface, akin to transient radio emission seen from magnetars. The radio burst confirms that the X-ray source is a neutron star and suggests a new link between these apparently radio-quiet X-ray emitting sources and other transient or persistent radio-emitting neutron stars. It also suggests that a common physical mechanism for emission might operate over a range of magnetic field strengths and neutron star ages. We propose that 2XMM J104608.7$-$594306 straddles the boundary between young, energetic neutron stars and their evolved radio-emitting cousins and may bridge these two populations. The detection of such a radio burst also shows that other radio-quiet neutron stars may also emit such sporadic radio emission that has been missed by previous radio surveys and highlights the need for regular monitoring of this unique sub-population of neutron stars.

The detection and characterization of habitable planets around nearby stars persist as one of the foremost objectives in contemporary astrophysics. This work investigates the synergistic integration of astrometric and direct imaging techniques by capitalizing on the complementary capabilities of the Closeby Habitable Exoplanet Survey (CHES) and Habitable Worlds Observatory (HWO). Planetary brightness and position vary over time due to phase effects and orbital architectures, information that can be precisely provided by CHES's astrometric measurements. By combining the precise orbital constraints from CHES with the imaging capabilities of HWO, we evaluate the improvements in detection efficiency, signal-to-noise ratio and overall planet yield. Completeness is quantified as the fraction of injected planets that are successfully detected, while yields are estimated for various scenarios using terrestrial planet occurrence rates derived from the Kepler dataset. Our results indicate that prior astrometric data significantly enhance detection efficiency. Under the adopted detection limit, our analysis indicates that prior CHES observations can increase completeness by approximately 10% and improve detection efficiency by factors ranging from two to thirty. The findings underscore the importance of interdisciplinary approaches in the search for and characterization of habitable worlds.

Several low--mass galaxy nuclei are observed to produce quasiperiodic eruptions (QPEs). Recently one of these systems, 1ES~1927+654, changed its quasiperiod drastically, from $\sim 18$ minutes to $\sim 7.1$ minutes, over a span of just two years. I suggest that this is an effect of von Zeipel -- Lidov -- Kozai (ZLK) cycles, where a more distant star orbits the QPE `binary' in which a white dwarf orbits a moderately massive central black hole. I show that in 1ES~1927+654 the white dwarf's orbital plane oscillates with angular amplitude $\simeq 71^{\circ}$ each side of the orbital plane of the distant star. This causes correlated changes of the orbital eccentricity and quasiperiod, and of the accretion luminosity driven by gravitational radiation losses. The GR luminosity has the characteristic property that it is inversely proportional to the instantaneous binary quasiperiod in all cases. The QPE system is probably just one of the effects produced by a complex infall event involving several stars. The whole system is likely to evolve rapidly, and will repay further monitoring.

Neutron stars (NS) offer an exceptional opportunity to investigate gravitational physics in extreme environments. In the context of alternative theories of General Relativity (GR), a scalar field non-minimally coupled to gravity can lead to "hairy" NS solutions with a non-zero scalar charge. As a result, binary systems containing at least one "hairy" NS can serve as probes for testing GR through observations. In this work, we analyze the gravitational wave signals from four neutron-star-black-hole (NS-BH) merger events: GW190814, GW200115, GW200105, and GW230529 examine the "hairy" NS hypothesis predicted by the luminal Horndeski theories. Our analysis constrains the scalar charge to be on the order of $10^{-3}$ to $10^{-4}$, with the most stringent limit arising from the recent GW230529 event. In addition, we discuss the implications of our results in relation to the scalarized NS solutions found in the literature.

We perform a detailed study of the cosmological constraints on the decay of a relic particle $\phi$ into neutrinos, $\phi \rightarrow \nu \bar{\nu}$, in particular those arising from the observed light-element abundances in the early Universe. We focus on the late-time disintegration of the light elements previously synthesised during BBN. Several processes are relevant, including final-state radiation associated with the decay, as well as subsequent interactions of the injected neutrinos with the thermal background neutrinos or between themselves. All processes generically contribute to the production of electromagnetic and often also hadronic material and may therefore induce late-time photodisintegration and hadrodisintegration reactions, i.e.~the destruction of light elements that have previously been formed during BBN. Here, we examine this scenario with a Monte-Carlo inspired probabilistic approach rather than Boltzmann techniques, taking into account all of these different reactions as well as their interplay. We find the resulting constraints to be very significant, covering a broad range of previously unexplored masses and lifetimes of the relic source particle.

Moist convection is a physical process where the latent heat released by condensation acts as a buoyancy source that can enhance or even trigger an overturning convective instability. Since the saturation temperature often decreases with height, condensation releases latent heat preferentially in regions of upflow. Due to this inhomogeneous heat source, moist convection may be more sensitive to changes in flow morphology, such as those induced by rotation, than dry Rayleigh-Bénard convection. In order to study the effects of rotation on flows driven by latent heat release, we present a suite of numerical simulations that solve the Rainy-Bénard equations (Vallis et al. 2019). We identify three morphological regimes: a cellular regime and a plume regime broadly analogous to those found in rotating Rayleigh Bénard convection, and a novel funnel regime that lacks a clear analog within the regimes exhibited by dry convection. We measure energy fluxes through the system and report rotational scalings of the Reynolds and moist Nusselt numbers. We find that moist static energy transport, as measured by a moist Nusselt number, is significantly enhanced in the funnel regime without a corresponding enhancement in Reynolds number, indicating that this funnel regime produces structures with more favorable correlations between the temperature and vertical velocity.

An ablative plasma shock can emanate from the interface between a cold/dense plasma and a hot/dilute ambient plasma, where the plasma mean-free-path is much longer than the temperature gradient length. The shock is driven by thermal flux from the hot plasma into the cold plasma, primarily through tail electrons mediated by an ambipolar electric field, and it propagates into the ambient hot/dilute plasma. Since the collisional mean-free-path is usually much longer than the Debye length, the ablative plasma shock is mostly collisionless, with the shock front width set by the upstream hot plasma Debye length and the shock speed by the downstream cold plasma sound speed. The shock heating of ions is extremely efficient via collisionless mixing of upstream hot ions and downstream cold ions, both of which have been converted into shock-front-bound flows accelerated by the ambipolar electric field that has a deep potential well anchored inside the shock front.

We derive a coordinate-invariant expression for the photon deflection angle in the strong deflection limit (SDL) of stationary axisymmetric spacetimes. The key logarithmic-divergence coefficient is shown to depend only on quantities locally measurable by a zero-angular-momentum observer -- curvature scalars, the circumferential radius, and the proper angular velocity. The same coefficient governs the damping rate of quasinormal modes (QNMs) in the eikonal limit, establishing a curvature-based, model-independent connection between QNMs and lensing in the SDL near rotating compact objects.

A new interpretation of Dirac singletons \cite{Dirac:1963ta}, i.e. free conformal fields in $d$ dimensions, as relativistic fields in a $d+1$-dimensional space-time with cosmological constant, that differs from the Flato-Fronsdal dipole construction in $AdS_{d+1}$ \cite{Flato:1986uh}, is proposed. The $d+1$-dimensional field is described at the level of both equations and Lagrangian. It forms an infinite-dimensional representation of the $d+1$-dimensional Lorentz group that relates fields at different space-time points. The associated well-known fact is that singleton cannot be localized at a point in ${d+1}$ dimensions, hence being unobservable via local scattering/radiation phenomena in the Standard Model ($d=3$). On the other hand, that singleton respects ${d+1}$ dimensional relativistic symmetries makes it possible to introduce its interactions with gravity and other relativistic fields in $d+1$ dimensions. It is speculated that the presence of singleton in a four-dimensional field theory with non-zero cosmological constant (dark energy) can be relevant to the dark matter phenomenon and baryon asymmetry generation.

Increasingly sensitive direct detection dark matter experiments are testing important regions of parameter space for WIMP dark matter and pushing many models to the multi-TeV regime. This brings into question the perturbativity of these models. In this context, and in light of the new limits from the LZ experiment, we investigate the status of the simplest thermal dark matter model: a singlet scalar, real or complex, coupled to the Higgs boson. We calculate the next-to-leading order (NLO) corrections to the direct detection rates as well as for the annihilations driving thermal freeze-out. For the complex case, we find that the entire perturbative region is excluded by direct detection. For the real case we find that the mass should be $\gtrsim 20\,{\rm TeV}$ at NLO, compared with the $\gtrsim 30\,{\rm TeV}$ LO limit. We highlight that a three-fold improvement on WIMP spin independent interactions can fully test the real scalar model in the perturbative regime. Finally, for both models, there is still an allowed (albeit narrow) region near the Higgs resonance where couplings are perturbative.

We present a deep neural network (DNN) accelerated Hamiltonian Monte Carlo (HMC) algorithm called DeepHMC for the inference of binary neutron star systems. The HMC is a non-random walk sampler that uses background gradient information to accelerate the convergence of the sampler. While faster converging than a random-walk sampler, in theory by a factor of the dimensionality of the problem, a known computational bottleneck for HMC algorithms is the calculation of gradients of the log-likelihood. We demonstrate that Hamiltonian trajectories based on a DNN gradients are 30 times faster than those based on the relative binning gradients, and 7000 times faster than trajectories based on a naive likelihood gradient calculation. Using the publicly available 128 second LVK data set for the binary neutron star mergers GW170817 and GW190425, we show that not only does DeepHMC produce produces highly accurate and consistent results with the LVK public data, but acquires 5000 statistically independent samples (SIS) in the $12D$ parameter space in approximately two hours on a Macbook pro for GW170817, with a cost of $<1$ second/SIS, and 2.5 days for GW190425, with a cost of $\sim25$ seconds/SIS.

We constrain the nuclear matter equation of state within the relativistic mean field model by including the isoscalar-vector and isovector-vector coupling term at a fundamental level using the Bayesian analysis. We used the nuclear saturation properties and recent astrophysical observations to constrain the dense matter equation of state. We obtained about 20000 sets of equations of states out of sampling about 60 millions sets of equations of states. All 20000 equations of states satisfy nuclear matter saturation properties at saturation densities and produces high mass neutron stars. In our findings, we find that the non-zero value of isoscalar-vector and isovector-vector coupling parameter and negative value of sigma meson self-coupling stiffen the equation of state. Our sets of equations of state produces neutron stars of mass larger than 2.5 M$_{\odot}$ to include the recent gravitational waves observation GW190419.

One in every two atoms in the Earth, Mars, and the Moon is oxygen; it is the third most abundant element in the solar system. The oxygen isotopic compositions of the terrestrial planets are different from those of the Sun and demonstrate that these planets are not direct compositional analogs of the solar photosphere. Likewise, the Suns O/Fe, Fe/Mg and Mg/Si values are distinct from those of inner solar system chondrites and terrestrial planets. These four elements (O, Fe, Mg, Si) make up about 94% by mass of the rocky planets and their abundances are determined uniquely using geophysical, geochemical and cosmochemical constraints. The rocky planets grew rapidly from planetesimals, most of which were differentiated, having a core and a mantle, before being accreted. Planetary growth in the early stages of protoplanetary disk evolution was rapid and was only partially recorded by the meteoritic record. The noncarbonaceous meteorites (NC) provide insights into the early history of the inner solar system and are used to construct a framework for how the rocky planets were assembled. NC chondrites have chondrule ages that are two to three million years younger than t_zero (the age of calcium-aluminum inclusions, CAI), documenting that chondrites are middle- to late-stage products of solar system evolution. The composition of the Earth, its current form of mantle convection, and the amount of radiogenic power that drives its engine remain controversial topics. Earths dynamics are driven by primordial and radiogenic heat sources. Measurement of the Earths geoneutrino flux defines its radiogenic power and restricts its bulk composition. Using the latest data from the KamLAND and Borexino geoneutrino experiments affirms that the Earth has 20 TW of radiogenic power and sets the proportions of refractory lithophile elements in the bulk silicate Earth at 2.5 to 2.7 times that in CI chondrites.

Overlapping signals represent one of the major data analysis challenges in next-generation gravitational wave detectors. We leverage Transformers and Normalizing Flows, state-of-the-art machine learning algorithms, to address the parameter estimation of overlapping binary black hole mergers in the Einstein Telescope (ET). Our proposed model combines a Transformer-based "Knowledge Extractor Neural Network" (KENN) with a Normalizing Flow (HYPERION) to perform rapid and unbiased inference over multiple overlapping black hole binary events. The choice of architecture leverages the strength of Transformers in capturing complex and long-range temporal structures in the strain time series data, while Normalizing Flows provide a powerful framework to sample posterior distributions. We demonstrate the effectiveness and robustness of our model over simulated gravitational wave signals, showing that it maintains the same level of accuracy regardless of the correlation level in the data. Moreover our model provides estimates of chirp mass and coalescence times within <10-20% from the true simulated value. The results obtained are promising and show how this approach might represent a first step toward a deep-learning based inference pipeline for ET and other future gravitational wave detectors.

Widespread solar energetic particle (SEP) events remain poorly understood phenomena in space weather. These events are often linked to coronal mass ejections (CMEs) and their shocks, but the mechanisms governing their global particle distribution remain debated. The 13 March 2023 event is particularly notable as a widespread SEP event associated with an exceptionally fast interplanetary shock. With speeds of up to 3000 km/s, it is one of the most extreme shocks observed in recent years. We aim to investigate whether the flanks of a wide CME-driven shock can decouple from the CME and continue propagating as freely propagating shock waves. If shocks are the primary SEP source, such a mechanism could help explain some of the widest SEP events. Using EUHFORIA, a 3D magnetohydrodynamic heliospheric model, we simulated the evolution of wide CME-driven shocks. We modified the model to allow direct shock injection at the inner boundary, upstream of the CME ejecta. Applying this to the 13 March 2023 event, we modelled two simultaneous CMEs whose shocks form a single, wide shock envelope that spans 280° in longitude. We then compared our results to in situ observations. Our simulations show that the flanks of wide CME shocks can persist as freely propagating waves beyond 2 au. For the 13 March 2023 event, the modelled shock arrival times and amplitudes of associated plasma parameters (e.g. speed and density) show good agreement with observations from various spacecraft distributed across different radial distances and longitudes. Furthermore, the combined shock structure expands into a quasi-circumsolar wave as it propagates outwards. These findings indicate that the shock flanks of fast CMEs can persist for a long time, supporting the idea that such freely propagating shock flanks play a key role in the global distribution of SEPs in widespread events.

Invoking the Klebanov-Susskind-Banks Euclidean wormhole as a bridge, we investigate the power spectrum and the entanglement between two pair-created universes. We construct a suitable global vacuum for the perturbations of the inflaton field in the Euclidean regime, which becomes a mixed state when restricted to one of the paired universes. This mixed state leads to an enhancement of the power spectrum for long-wavelength modes. In addition, entanglement between the two universes is realized by the existence of the wormhole. Thus, the power spectrum enhancement in the long-wavelength regime might be evidence of our universe being created from a Euclidean wormhole that was entangled with a partner universe, and hence our universe does not begin with a pure state.