Papers for Thursday, Nov 14 2024

A core-corona decomposition of compact (neutron) star models is compared to recent NICER data of masses and radii. It is in particular interesting to capture the outlier XTE~J1814-338. Instead of integrating the TOV equations from the center to surface, we follow here another pathway by accommodating all uncertainties of the equation(s) of state (EoS) at supra-nuclear density or/and an unknown dark matter admixture in a parameterization of the core by its radius $r_x$, the included mass $m_x$ and the pressure $p_x$ at $r_x$. The corona, which may be dubbed also envelope or halo or outer crust, is assumed to be of standard-model matter where the EoS is supposed to be faithfully known.

The paper explores the use of various machine learning methods to search for heterogeneous or atypical structures on astronomical maps. The study was conducted on the maps of the cosmic microwave background radiation from the Planck mission obtained at various frequencies. The algorithm used found a number of atypical anomalous structures in the actual maps of the Planck mission. This paper details the machine learning model used and the algorithm for detecting anomalous structures. A map of the position of such objects has been compiled. The results were compared with known astrophysical processes or objects. Future research involves expanding the dataset and applying various algorithms to improve the detection and classification of outliers.

Cosmological data collected on a sphere, such as CMB anisotropies, are typically represented by the spherical harmonic coefficients, denoted as $a_{\ell m}$. The angular power spectrum, or $C_\ell$, serves as the fundamental estimator of the variance in this data. Alternatively, spherical data and their variance can also be characterized using Multipole Vectors (MVs) and the Fréchet variance. The vectors that minimize this variance, known as Fréchet Vectors (FVs), define the center of mass of points on a compact space, making them highly sensitive to small displacements of these points. This sensitivity makes FVs excellent indicators of statistical correlations between different multipoles. We demonstrate this using both simulations and real data. Through simulations, we show that FVs enable a blind detection and reconstruction of the location associated with a mock Cold Spot anomaly introduced in an otherwise isotropic sky. Applying this to the 2018 Planck maps, we implement several improvements on previous model-independent tests of Gaussianity and statistical isotropy, down to arc-minute scales. When compared with simulated maps that incorporate masking and anisotropic noise, for $2 \leq\ell \leq 1500$, while Planck's MVs appear consistent with these hypotheses, the corresponding FVs reject them with significances between 5.2 and $8.3\sigma$, depending on the component separation method.

Ultraviolet flare emission can drive photochemistry in exoplanet atmospheres and even serve as the primary source of uncertainty in atmospheric retrievals. Additionally, flare energy budgets are not well-understood due to a paucity of simultaneous observations. We present new near-UV (NUV) and optical observations of flares from three M dwarfs obtained at 20 s cadence with Swift and TESS, along with a re-analysis of flares from two M dwarfs in order to explore the energy budget and timing of flares at NUV--optical wavelengths. We find a 9000 K blackbody underestimates the NUV flux by $\geq$2$\times$ for 54$\pm$14% of flares and 14.8$\times$ for one flare. We report time lags between the bands of 0.5--6.6 min and develop a method to predict the qualitative flare shape and time lag to 36$\pm$30% accuracy. The scatter present in optical-NUV relations is reduced by a factor of 2.0$\pm$0.6 when comparing the total NUV energy with the TESS energy during the FWHM duration due to the exclusion of the $T_\mathrm{eff}\approx$5000 K tail. We show the NUV light curve can be used to remove flares from the optical light curve and consistently detect planets with 20% smaller transits than is possible without flare detrending. Finally, we demonstrate a 10$\times$ increase in the literature number of multi-wavelength flares with the Early eVolution Explorer (EVE), an astrophysics Small Explorer concept to observe young clusters with simultaneous NUV and optical bands in order to detect young planets, assess their photochemical radiation environments, and observe accretion.

The Galactic bulge hosts the Milky Way's oldest stars, possibly coming from disrupted globular clusters (GCs) or the bulge's primordial building blocks, making these stars witnesses to the Galaxy's early chemical enrichment. The Galactic bar currently dominates the bulge's region, altering the orbits of objects formed before its formation and complicating the trace of the field stars' original clusters. Here, we present the discovery of a fossil record of this evolution, SOS1 -- a star trapped in the bar, exhibiting significant enhancements in nitrogen, sodium, and aluminum, typical of second-generation GC stars. SOS1 also shows an s-process Ce enhancement, suggesting an old age and early enrichment by fast-rotating massive stars in the Galaxy's earliest phases. With the purpose of finding the SOS1's parent GC, we derive its precise chemodynamical properties by combining high-precision proper motions from Gaia with APOGEE detailed chemical abundances. Our analysis suggests that SOS1 was possibly stripped from the GC Terzan 5 by the Galactic bar's gravitational influence approximately 350 Myr ago. We also found chemical similarities suggesting that SOS1 belonged to the most metal-poor, ancient, and peripheral stellar population of Terzan 5. These results not only support the hypothesis that Terzan 5 is a remnant of a primordial building block of the Galactic bulge, but also suggest this cluster continues losing stars to the bar. Our method highlights how powerful the use of chemodynamical properties in the Gaia era is for tracing the Galaxy's evolutionary history.

Compact binaries containing hot subdwarfs and white dwarfs have the potential to evolve into a variety of explosive transients. These systems could also explain hypervelocity runaway stars such as US 708. We use the detailed binary evolution code MESA to evolve hot subdwarf and white dwarf stars interacting in binaries. We explore their evolution towards double detonation supernovae, helium novae, or double white dwarfs. Our grid of 3120 models maps from initial conditions such as orbital period and masses of hot subdwarf and white dwarf to these outcomes. The minimum amount of helium required to ignite the helium shell that leads to a double detonation supernova in our grid is $\approx 0.05 \, \mathrm{M_{\odot}}$, likely too large to produce spectra similar to normal type Ia supernovae, but compatible with inferred helium shell masses from some observed peculiar type I supernovae. We also provide the helium shell masses for our double white dwarf systems, with a maximum He shell mass of $\approx 0.18\,\mathrm{M_{\odot}}$. In our double detonation systems, the orbital velocity of the surviving donor star ranges from $\approx 450 \, \mathrm{km\,s^{-1}}$ to $\approx 1000 \, \mathrm{km\,s^{-1}}$. Among the surviving donors, we also estimate the runaway velocities of proto-white dwarfs, which have higher runaway velocities than hot subdwarf stars of the same mass. Our grid will provide a first-order estimate of the potential outcomes for the observation of binaries containing hot subdwarfs and white dwarfs from future missions like Gaia, LSST, and LISA.

Over the past decades, a population of galaxies invisible in optical/near-infrared, but bright at longer wavelengths, have been identified through color selections. These so-called optically faint/dark galaxies are considered to be massive quiescent galaxies or highly dust-attenuated galaxies. Having the entire galaxy obscured by dust, however, is likely an extreme case of the much more common occurrence of optically thin and thick absorption coexisting in the same system. With the power of JWST imaging, we are able to spatially resolve massive galaxies at z~3, accurately model their spectral energy distributions, and identify candidate optically thick substructures. We target galaxies with log(M*/Msun)>10.3 and 2.5<z<3.5, and get 486 galaxies in CEERS and PRIMER fields. Based on excess NIR luminosity, we identify 162 galaxies (~33\% of the parent sample) as candidate hosts of optically thick substructures. We then carry out spatially resolved SED modeling to explore the physical properties of those dark substructures and estimate the amount of optically thick obscuration. We find that optically thick dust is ubiquitous in normal massive galaxies with a wide variety of SFR and morphology. 10-20\% of the stellar mass/SFR are unaccounted for in our selected galaxies, and the fraction is insensitive to stellar mass or SFR. The dark substructures are generally dustier than the rest of the galaxies and are irregularly distributed, arguing against an obscured AGN as the source of the NIR excess. A correlation between the obscured luminosity and the presence of a recent starburst in the past <100 Myr is also observed.

We present the results of Chandra and XMM-Newton X-ray imaging and spatially-resolved spectroscopy, as well as new MUSTANG2 90~GHz observations of the thermal Sunyaev-Zeldovich from MACS J0717.5+3745, an intermediate redshift ($z = 0.5458$) and exceptionally massive ($3.5 \pm 0.6 \times 10^{15}~\rm M_\odot$) Frontier Fields cluster experiencing multiple mergers and hosting an apparent X-ray bright large scale structure filament. Thermodynamical maps are produced from Chandra, XMM-Newton, and ROSAT data using a new method for modelling the astrophysical and instrumental backgrounds. The temperature peak of $24 \pm 4$ keV is also the pressure peak of the cluster and closely correlates spatially with the Sunyaev-Zeldovich peak from the MUSTANG2 data. The cluster center hosts shock fronts to the north and south, for which we report lower limits for the shock Mach numbers of $M = 1.6 \pm 0.4$ and $M = 1.9 \pm 0.3$, respectively. Bayesian X-ray Analysis methods were used to disentangle different projected spectral signatures for the filament structure, with Akaike and Bayes criteria being used to select the most appropriate model to describe the various temperature components. We report an X-ray filament temperature of $2.9_{-0.3}^{+0.5}$ keV and a density $(1.60\pm0.05)\times10^{-4}\,{\rm cm^{-3}}$, corresponding to an overdensity of 150 relative to the critical density of the Universe. We estimate the hot gas mass of the filament to be $\sim4.4\times10^{12}~\rm M_\odot$, while its total projected weak lensing measured mass is $\sim 6.8 \pm 2.7 \times 10^{13}~\rm M_\odot$, indicating a hot baryon fraction of 4-10\%.

We characterise the co-evolution of radio-loud AGN and their galaxies by mapping the dependence of radio-loud AGN activity on stellar mass and star-formation rate (SFR) across cosmic time (out to $z \sim 1.5$). Deep LOFAR radio observations are combined with large galaxy samples to study the incidence of radio-loud AGN across the galaxy population; the AGN are further split into low-excitation radio galaxies (LERGs) and high-excitation radio galaxies (HERGs). We find that LERG activity occurs over a wide range of SFRs, whereas HERGs are typically found in galaxies with ongoing star formation. The LERGs are then split based on their SFRs relative to the main sequence, across redshift. Within quiescent galaxies, LERG activity shows a steep stellar mass dependence with the same normalisation across the past $\sim$ 10 Gyr; this indicates that hot gas fuels LERGs in quiescent galaxies across cosmic time. In massive galaxies ($\log_{10}(M/\rm{M_{\odot}}) \gtrsim 11$), the incidence of LERGs is roughly constant across the galaxy population, suggesting that LERGs in massive galaxies may be fuelled by hot gas regardless of the star-formation activity. At lower masses, however, LERG activity is significantly more enhanced (by a factor of up to 10) in star-forming galaxies compared to quiescent galaxies; this suggests that an additional fuelling mechanism, likely associated with cold gas, may fuel the LERGs in galaxies with higher SFRs. We find that HERGs typically accrete above 1 per cent of the Eddington-scaled accretion rate, and the LERGs typically accrete below this level.

We report a novel pilot project to characterise intra-night optical variability (INOV) of an extremely rare type of quasar, which has recently been caught in the act of transiting from a radio-quiet to radio-loud state, on a decadal time scale. Such rare transitions may signify a recurrence, or conceivably the first switch-on of jet activity in optically luminous quasars. The newly formed jet could well be jittery and unsteady, both in power and direction. The optically brightest among such radio-state transition candidates, the quasar J0950+5128 ($z = 0.2142$), was monitored by us with dense sampling in the R-band, during 2020-21 in 6 sessions, each lasting $>$ 4 hours. This is the first attempt to characterise the INOV properties associated with this recently discovered, extremely rarely observed phenomenon of quasar radio-state transition. The non-detection of INOV in any of the 6 sessions, down to the 1-2% level, amounts to a lack of evidence for a blazar-like optical activity, $\sim$ 2 years after its transition to radio-loud state was found. The only INOV feature detected in J0950+5128 during our observational campaign was a $\sim$ 0.15-mag spike lasting < 6 minutes, seen at 13.97 UT on 18-March-2021. We also report the available optical light curves of this quasar from the Zwicky Transient Facility (ZTF) survey, which indicate that it had experienced a phase of INOV activity around the time its transition to the radio-loud state was detected, however that phase did not sustain until the launch of our INOV campaign $\sim$ 2 years later.

Hydrogen, the most abundant element in the Universe, has traditionally been used to investigate astrophysical processes within and around our own Galaxy. In its chemically neutral, atomic form (known as HI in the astronomical literature), it has tremendous potential today as a tool for precision cosmology and testing theories of fundamental physics. Cosmological HI is accessed through two of its main spectral lines: the Lyman-$\alpha$, with a rest wavelength of 1216 $Å$, in the ultraviolet and visible part of the spectrum, and the 21-cm, which manifests in the radio frequency band. A plethora of radio telescopes worldwide are focused on detecting the faint 21 cm signal from the dark ages and Cosmic Dawn, some of the earliest epochs of the Universe. This chapter will describe the formalism for doing cosmology with HI, the recent results from the facilities and their prospects for studying the evolution of the Universe.

We present several nonlinear wavefront sensing techniques for few-mode sensors, all of which are empirically calibrated and agnostic to the choice of wavefront sensor. The first class of techniques involves a straightforward extension of the linear phase retrieval scheme to higher order; the resulting Taylor polynomial can then be solved using the method of successive approximations, though we discuss alternate methods such as homotopy continuation. In the second class of techniques, a model of the WFS intensity response is created using radial basis function interpolation. We consider both forward models, which map phase to intensity and can be solved with nonlinear least-squares methods such as the Levenberg-Marquardt algorithm, as well as backwards models which directly map intensity to phase and do not require a solver. We provide demonstrations for both types of techniques in simulation using a quad-cell sensor and a photonic lantern wavefront sensor as examples. Next, we demonstrate how the nonlinearity of an arbitrary sensor may studied using the method of numerical continuation, and apply this technique both to the quad-cell sensor and a photonic lantern sensor. Finally, we briefly consider the extension of nonlinear techniques to polychromatic sensors.

We investigate how magnetic field variations around accreting black holes on event horizon scales affect the morphology of magnetically-driven jet on larger scales. By performing radiative transfer calculations on general relativistic magnetohydrodynamics simulations, we find that temporal variation in the magnetic flux on the event horizon and the jet power are imprinted on the variability of jet width up to several hundred gravitational radii. When the magnetic flux around the black hole drops and then rises, the jet initially narrows or becomes truncated, then widens, creating a thin-thick pattern that propagates down the jet. This suggests that extended jet observations can provide a history record of horizon-scale magnetic field dynamics, and conversely, upcoming changes in the jet image can be predicted from direct observation of the magnetized accreting plasma near the black hole. Furthermore, the pattern of jet width variations shows acceleration up to the relativistic regime as it moves away from the black hole, aligning with plasma bulk motion. We also find in time-averaged images that both the bulk plasma motion and magnetic field configuration in the jet-launching region, which are sensitive to black hole spin, shape diverse features through relativistic beaming and aberration. Higher black hole spins result in more poloidal bulk motion and toroidal magnetic fields, leading to more symmetric jet images and linear polarization patterns. These results suggest a new method for testing the magnetically arrested disk model and the Blandford-Znajek process, and for determining the black hole spin through observations bridging horizon and jet-launching scales.

We investigate the formation history of intrahalo light (IHL) using the high-resolution (~1 kpc), large-scale (~Gpc) cosmological hydrodynamical simulation, Horizon Run 5 (HR5). IHL particles are identified by carefully considering both their binding energies and positions with respect to the tidal radii of individual galaxies. By analyzing more than 1,200 galaxy groups and clusters with $\geq 10^{13} M_{\odot}$ and tracing their individual IHL particles back in time, we classify the origin of each IHL particle at each epoch based on the status of the originating galaxy into three categories: brightest halo galaxy (BHG) formation/merger, satellite galaxy stripping, and pre-processing. Our study reveals that the IHL production through BHG formation/merger is the predominant production channel, contributing over 60\% of the total IHL mass across all redshifts. The second most significant IHL production channel is pre-processing, providing more than 20\% in the final HR5 snapshot. Stripping is negligible at $z>4$ but becomes gradually more important as halos mature at $z<4$. Finally, we verify that IHL production through the disruption of dwarf galaxies and in-situ formation is negligible, contributing less than ~3\% and ~0.5\% to the total IHL production, respectively.

Coupled-mode theory (CMT) is a powerful tool for simulating near-harmonic systems. In telecommunications, variations of the theory have been used extensively to study waveguides, both analytically and through numerical modelling. Analogous mathematical techniques to the CMT are also widely used in quantum mechanics. The purpose of this work is to collect different formulations of the CMT and their underlying connections to quantum mechanical techniques, and to showcase their utility in modelling slowly varying waveguides including directional couplers and photonic lanterns. My choice of example waveguides is motivated by the astronomical applications of such devices in starlight nulling, wavefront sensing, and high-resolution spectroscopy. I first provide a brief review of the standard form of the CMT, applicable for waveguides with fixed eigenmodes. Next, I show that the CMT also applies for slowly varying waveguides, and demonstrate the close relation between the CMT and several well-known approximation methods from quantum mechanics, as well as concepts like geometric phase. Finally, I present a verification of my analysis, in the form of the numerical package cbeam.

Binaries are known to play a key role in the mass loss and dynamical environments of evolved stars. Stellar and sub-stellar companion interactions produce complex wind morphologies including rotating/expanding disks, bipolar outflows, and spiral wind patterns; however, the connection between these many structures and the gas phase chemistry they harbor is not well-constrained. To expand the sample of chemical inventories in interacting systems, we present a detailed spectroscopic case study of the binary C-rich Asymptotic Giant Branch (AGB) star V Hya. Using spatially resolved ALMA observations at Bands 3, 6 and 7, we characterize the rotational emission lines and distributions of molecules in its surrounding disk undergoing dynamical expansion (DUDE). We detect emission from over 15 molecules and isotopologues toward this source, and present resolved maps for the brightest tracers of carbonaceous chemistry (e.g. CCH, C4H, HC5N, HNC, CH3CN). Employing LTE and non-LTE models of emission from the DUDE, we estimate the abundance distributions for optically thin species, and compare them with prototypical carbon-rich AGB envelopes. We find that the average abundances of detected species are within a factor of ${\sim}5$ from sources with similar mass-loss rates; however, the distribution of daughter species in V Hya is much more compact, with carbon chain species (CCH, C4H, HC3N) appearing with abundances $>$10$^{-7}$ even in the innermost sampled regions (200 au) of the disk.

In this work, we devise a model for estimating UV and optical line emission (i.e., CIII] $1909$A, H$\beta$, [OIII] $5007$A, H$\alpha$, [NII] $6583$A) tracing HII regions in the interstellar medium (ISM) of galaxies at $z\sim4-6$ from the ALMA Large Programme ALPINE. The aim is to investigate the impact of binary stars in the stellar population along with an abrupt quenching in the Star Formation History (SFH) on line emission. This is crucial for understanding the ISM's properties in early galaxies and identifying new star formation tracers in high-$z$ galaxies. The model simulates HII+Photodissociation Region (PDR) complexes through radiative transfer in 1D slabs, characterized by gas density ($n$), ionisation parameter ($U$), and metallicity ($Z$). It considers: (a) heating from star formation (SF), simulated with Starburst99 and BPASS to quantify binary stars impact; (b) constant, exponentially declining, and quenched SFH scenarios. For each galaxy, we select theoretical ratios from CLOUDY models between [CII] line emission, tracing PDRs, and nebular lines from HII regions, using these to derive expected optical/UV lines from observed [CII]. We find binary stars strongly impact line emission post-quenching, keeping UV photon flux higher for longer, maintaining free electron temperature and ionised column density in HII regions up to 5 Myr after quenching. We constrain ISM properties of our subsample, finding a low ionisation parameter $\log U{\approx}-3.8\pm 0.2$ and moderate/high densities $\log(n/\rm cm^{-3}){\approx}2.9\pm 0.6$. Finally, we derive UV/optical line luminosities-SFR relations for different burstiness parameters ($k_s$). In the fiducial BPASS model, relations show negligible SFH dependence but depend on $k_s$, while in the SB99 case, dependence is on SFH. We propose their use for characterising the burstiness of high-$z$ galaxies.

In 2018 the EDGES experiment claimed the first detection of the global cosmic 21cm signal, which featured an absorption trough centered around $z \sim 17$ with a depth of approximately -500mK. This amplitude is deeper than the standard prediction (in which the radio background is determined by the cosmic microwave background) by a factor of two and potentially hints at the existence of a radio background excess. While this result was obtained by fitting the data with a phenomenological flattened-Gaussian shape for the cosmological signal, here we develop a physical model for the inhomogeneous radio background sourced by the first galaxies hosting population III stars. Star formation in these galaxies is quenched at lower redshifts due to various feedback mechanisms, so they serve as a natural candidate for the excess radio background hinted by EDGES, without violating present day measurements by ARCADE2. We forward-model the EDGES sky temperature data, jointly sampling our physical model for the cosmic signal, a foreground model, and residual calibration errors. We compare the Bayesian evidences obtained by varying the complexity and prior ranges for the systematics. We find that the data is best explained by a model with seven log-polynomial foreground terms, and that it requires calibration residuals. Interestingly, the presence of a cosmic 21cm signal with a non-standard depth is decisively disfavored. This is contrary to previous EDGES analysis in the context of extra radio background models, serving as a caution against using a ''pseudo-likelihood'' built on a model (flattened Gaussian) that is different from the one being used for inference. We make our simulation code and associated emulator publicly-available.

The standard $\Lambda$CDM cosmological model informed by cosmic microwave background (CMB) anisotropies makes a precise prediction for the growth of matter density fluctuations over cosmic time on linear scales. A variety of cosmological observables offer independent and complementary ways of testing this prediction, but results have been mixed, with many constraints on the amplitude of structure $S_8$ being 2-3$\sigma$ lower than the expectation from Planck primary CMB anisotropies. It is currently unclear whether these discrepancies are due to observational systematics, non-linearities and baryonic effects or new physics. We review how gravitational lensing of the CMB has and will continue to provide insights into this problem, including through tomographic cross-correlations with galaxy surveys over cosmic time.

Context: Understanding galaxy evolution in dense environments, particularly proto-clusters, is crucial for studying mechanisms driving star formation and quenching. Aims: This study examines how two proto-cluster over-densities at 3 < z < 4 impact star formation rate (SFR), stellar mass, and morphology, focusing on quenched galaxies. Methods: We identified proto-cluster over-densities in the Chandra Deep Field South (CDFS) and Ultra Deep Survey (UDS) regions of the VANDELS survey. Using spectral energy distribution analysis, Bayesian methods (BEAGLE and BAGPIPES) helped derive best-fit parameters and U-V and V-J rest-frame colours (UVJ), classifying galaxies as quenched or star-forming based on UVJ diagrams and specific star formation rates (sSFR). TNG300 simulations aided interpretation. Results: Two of 13 proto-cluster over-densities host quenched galaxies with red U-V colours, low sSFR, and properties like massive passive galaxies. These quenched members are redder, older, more massive, and more compact. The highest-density peaks at z=3.55 and z=3.43 have dark matter halo masses consistent with proto-clusters and host AGNs, with five and three AGNs, respectively. Compared to field galaxies, these quenched members are in denser environments. TNG300 simulations suggest proto-clusters with quenched galaxies at high redshift evolve to contain more passive galaxies by z=1. Conclusions: The over-densities host massive quenched galaxies and AGNs in their densest peaks. Simulations reveal that sSFR for passive galaxies in proto-clusters was high at z=6, with median mass growth rates of 96% from z=6 to z=3. Conditions for mass assembly likely involve galaxy interactions and high gas accretion in dense environments. Black hole growth and AGN feedback appear to drive quenching at z=3, aligning with the properties of quenched galaxies observed in our study.

We present the discovery of TOI-6303b and TOI-6330b, two massive transiting super-Jupiters orbiting a M0 and a M2 star respectively, as part of the Searching for GEMS survey. These were detected by TESS and then confirmed via ground-based photometry and radial velocity observations with the Habitable-zone Planet Finder (HPF). TOI-6303b has a mass of 7.84 +/- 0.31 MJ, a radius of 1.03 +/- 0.06 RJ , and an orbital period of 9.485 days. TOI-6330b has a mass of 10.00 +/- 0.31 MJ , a radius of 0.97 +/- 0.03 RJ , and an orbital period of 6.850 days. We put these planets in context of super-Jupiters around M-dwarfs discovered from radial-velocity surveys, as well as recent discoveries from astrometry. These planets have masses that can be attributed to two dominant planet formation mechanisms - gravitational instability and core-accretion. Their masses necessitate massive protoplanetary disks that should either be gravitationally unstable, i.e. forming through gravitational instability, or be amongst some of the most massive protoplanetary disks to form objects through core-accretion. We also discuss the eccentricity distribution of these objects, as a potential indicator of their formation and evolutionary mechanisms.

The Solar Neutrino and Astro-Particle PhYsics (SNAPPY) Cubesat is expected to launch in 2025 and it will carry into a polar orbit a prototype test detector for solar neutrino background studies while over the Earth's poles for the neutrino Solar Orbiting Laboratory future project ($\nu$SOL). During this flight it is possible to do other science measurements. One of these is an improved study of the solar wind particles through improved particle identification and energy measurements. This study aimed to understand how well could the solar wind particles be identified using the planned detector but instead of using the veto array as an anti-coincidence it would be used as a $\Delta$E energy sampling of a phoswich particle ID system.

Owing to their rapid cooling rate and hence loss-limited propagation distance, cosmic-ray electrons and positrons (CRe) at very high energies probe local cosmic-ray accelerators and provide constraints on exotic production mechanisms such as annihilation of dark matter particles. We present a high-statistics measurement of the spectrum of CRe candidate events from 0.3 to 40 TeV with the High Energy Stereoscopic System (H.E.S.S.), covering two orders of magnitude in energy and reaching a proton rejection power of better than $10^{4}$. The measured spectrum is well described by a broken power law, with a break around 1 TeV, where the spectral index increases from $\Gamma_1 = 3.25$ $\pm$ 0.02 (stat) $\pm$ 0.2 (sys) to $\Gamma_2 = 4.49$ $\pm$ 0.04 (stat) $\pm$ 0.2 (sys). Apart from the break, the spectrum is featureless. The absence of distinct signatures at multi-TeV energies imposes constraints on the presence of nearby CRe accelerators and the local CRe propagation mechanisms.

Free electrons in the interstellar medium refract and diffract radio waves along multiple paths, resulting in angular and temporal broadening of radio pulses that limits pulsar timing precision. We determine multifrequency, multi-epoch scattering times for the large dispersion measure millisecond pulsar J1903+0327 by developing a three component model for the emitted pulse shape that is convolved with a best fit pulse broadening function (PBF) identified from a family of thin-screen and extended-media PBFs. We show that the scattering time, $\tau$, at a fiducial frequency of 1500 MHz changes by approximately 10% over a 5.5yr span with a characteristic timescale of approximately 100 days. We also constrain the spectral index and inner scale of the wavenumber spectrum of electron density variations along this line of sight. We find that the scaling law for $\tau$ vs. radio frequency is strongly affected by any mismatch between the true and assumed PBF or between the true and assumed intrinsic pulse shape. We show using simulations that refraction is a plausible cause of the epoch dependence of $\tau$, manifesting as changes in the PBF shape and $1/e$ time scale. Finally, we discuss the implications of our scattering results on pulsar timing including time of arrival delays and dispersion measure misestimation.

Ever since they were first detected in the interstellar medium, the radio wavelength (3.3 GHz) hyperfine-structure splitting transitions in the rotational ground state of CH have been observed to show anomalous excitation. Astonishingly, this behaviour has been uniformly observed towards a variety of different sources probing a wide range of physical conditions. While the observed level inversion can be explained globally by a pumping scheme involving collisions, a description of the extent of 'over-excitation' observed in individual sources requires the inclusion of radiative processes, involving transitions at higher rotational levels. Therefore, a complete description of the excitation mechanism in the CH ground state, observed towards individual sources entails observational constraints from the rotationally excited levels of CH and in particular that of its first rotationally excited state. Given the limited detections of these lines, the objective of this work is to characterise the physical and excitation properties of the rotationally excited lines of CH near 700 MHz, and investigate their influence on the pumping mechanisms of the ground-state lines of CH. This work presents the first interferometric search for the rotationally excited lines of CH near 700 MHz carried out using the uGMRT array and jointly models the physical and excitation conditions traced by lines from both the ground and first rotationally excited states of CH.

Context. Solar energetic particles (SEPs) are detected in interplanetary space in association with flares and coronal mass ejections (CMEs) at the Sun. The magnetic connection between the observing spacecraft and the solar active region (AR) source of the event is a key parameter in determining whether SEPs are observed and the properties of the particle event. Aims. We investigate whether an east-west asymmetry in the detection of SEP events is present in observations and discuss its possible link to corotation of magnetic flux tubes with the Sun. Methods. We used a published dataset of 239 CMEs recorded between 2006 and 2017 and having source regions both on the front side and far side of the Sun as seen from Earth. We produced distributions of occurrence of in-situ SEP intensity enhancements associated with the CME events, versus \Delta \phi, the separation in longitude between the source active region and the magnetic footpoint of the observing spacecraft based on the nominal Parker spiral. We focused on protons of energy >10 MeV measured by the STEREO A, STEREO B and GOES spacecraft at 1 au. We also considered the occurrence of 71-112 keV electron events detected by MESSENGER between 0.31 and 0.47 au. Results. We find an east-west asymmetry in the detection of >10 MeV proton events and of 71-112 keV electron events. For protons, observers for which the source AR is on the east side of the spacecraft footpoint and not well connected (-180 < \Delta \phi < -40) are 93% more likely to detect an SEP event compared to observers with +40 < \Delta \phi < +180. The asymmetry may be a signature of corotation of magnetic flux tubes with the Sun, given that for events with \Delta \phi < 0 corotation sweeps the particle-filled flux tubes towards the observing spacecraft, while for \Delta \phi > 0 it takes them away from it.

We analyse continuum and molecular emission, observed with ALMA, from the dust-enshrouded intermediate-mass AGB star OH 30.1 -0.7. We find a secondary peak in the continuum maps, "feature B", separated by 4.6" from the AGB star, which corresponds to a projected separation of $1.8 \times 10^4$ au, placing a lower limit on the physical separation. This feature is most likely composed of cold dust and is likely to be ejecta associated with the AGB star, though we cannot rule out that it is a background object. The molecular emission we detect includes lines of CO, SiS, CS, SO$_2$, NS, NaCl, and KCl. We find that the NS emission is off centre and arranged along an axis perpendicular to the direction of feature B, indicative of a UV-emitting binary companion (e.g. a G-type main sequence star or hotter), perhaps on an eccentric orbit, contributing to its formation. However, the NaCl and KCl emission constrain the nature of that companion to not be hotter than a late B-type main sequence star. We find relatively warm emission arising from the inner wind and detect several vibrationally excited lines of SiS (3 = 1), NaCl (up to 3 = 4) and KCl (up to 3 = 2), and emission from low energy levels in the mid to outer envelope, as traced by SO$_2$. The CO emission is abruptly truncated around 3.5" or 14,000 au from the continuum peak, suggesting that mass loss at a high rate may have commenced as little as 2800 years ago.

The cosmic velocity field is an unbiased probe of the total matter distribution but is challenging to measure directly at intermediate and high redshifts. The large-scale velocity field imprints a signal in the cosmic microwave background (CMB) through the kinetic Sunyaev-Zeldovich (kSZ) effect. We perform the first 3d reconstruction of the large-scale velocity field from the kSZ effect by applying a quadratic estimator to CMB temperature maps and the 3d positions of galaxies. We do so by combining CMB data from the fifth data release of the Atacama Cosmology Telescope (in combination with Planck) and a spectroscopic galaxy sample from the Sloan Digital Sky Survey. We then measure the galaxy-velocity cross-power spectrum and detect the presence of the kSZ signal at a signal-to-noise ratio of 7.2$\sigma$. Using this galaxy-velocity cross-correlation alone, we constrain the amplitude of local primordial non-Gaussianity finding $f_{\rm NL}=-90^{+210}_{-350}$. This pathfinder measurement sets the stage for joint galaxy-CMB kSZ constraints to significantly enhance the $f_{\rm NL}$ information obtained from galaxy surveys through sample variance cancellation.

In weak magnetic fields ($\lesssim 50 \,\mbox{G}$), parallel and perpendicular viscosities, mainly from neutrals, may exceed magnetic diffusivities (Ohm, Hall, ambipolar) in the middle and upper chromosphere. Ion-driven gyroviscosity may dominate in the upper chromosphere and transition region. In strong fields ($\gtrsim 100\, \mbox{G}$), viscosities primarily exceed diffusivities in the upper chromosphere and transition region. Parallel and perpendicular viscosities, being similar in magnitude, dampen waves and potentially compete with ambipolar diffusion in plasma heating, potentially inhibiting Hall and ambipolar instabilities when equal. The perpendicular viscosity tensor has two components, $\nu_1$ and $\nu_2$, which differ slightly and show weak dependence on ion magnetization. Their differences, combined with shear, may destabilize waves, though magnetic diffusion introduces a cutoff for this instability. In configurations with a magnetic field $\bf{B}$ having vertical ($b_z=B_z/|\bf{B}|$) and azimuthal ($b_y=B_y/|\bf{B}|$) components, and a wavevector $\bf{k}$ with radial ($\kx=k_x/|\bf{k}|$) and vertical ($\kz=k_z/|\bf{k}|$) components, parallel viscosity and Hall diffusion can generate the viscous-Hall instability. Gyroviscosity further destabilizes waves in the upper regions. These findings indicate that the solar atmosphere may experience various viscous instabilities, revealing complex interactions between viscosity, magnetic fields, and plasma dynamics across different atmospheric regions.

We investigate when and how the relations of galaxy morphology and star forming activity with clustercentric radius become evident in galaxy clusters. We identify 162 galaxy clusters with total mass $M_{\rm tot}^{\rm cl} > 5 \times 10^{13} {\rm M}_\odot$ at $z = 0.625$ in the Horizon Run 5 (HR5) cosmological hydrodynamical simulation and study how the properties of the galaxies with stellar mass $M_\ast > 5 \times 10^9 {\rm M}_\odot$ near the cluster main progenitors have evolved in the past. Galaxies are classified into disk, spheroid, and irregular morphological types according to the asymmetry and Sersic index of their stellar mass distribution. We also classify galaxies into active and passive ones depending on their specific star-formation rate. We find that the morphology-clustercentric radius relation (MRR) emerges at $z \simeq 1.8$ as the fraction of spheroidal types exceeds 50% in the central region ($d \lesssim 0.1 R_{200}$). Galaxies outside the central region remain disk-dominated. Numerous encounters between galaxies in the central region seem to be responsible for the morphology transformation from disks to spheroids. We also find that the star formation activity-clustercentric radius relation emerges at an epoch different from that of MRR. At $z\simeq0.8$, passive galaxies start to dominate the intermediate radius region ($0.1\lesssim d/R_{200} \lesssim0.3$) and this "quenching region" grows inward and outward thereafter. The region dominated by early-type galaxies (spheroids and passive disks) first appears at the central region at $z\simeq 1.8$, expands rapidly to larger radii as the population of passive disks grows in the intermediate radii, and clusters are dominated by early types after $z\simeq 0.8$.

We investigate the evolution of a changing-look blazar (CLB) on long timescales and expect to trace the state change of a CLB. Three morphological types, including a flat spectrum radio quasar (FSRQ) state, transition state, and BL Lacertae (BL Lac) state are classified according to the criteria proposed by analyzing the relationship between the equivalent width of the emission line and the $\gamma$-ray photon spectral index $\Gamma_{\gamma}$. The multiwavelength light curves and spectral energy distributions corresponding to different epochs are obtained. The efforts found that $\Gamma_{\gamma}$ satisfy the relationships with $\Gamma_{\gamma} $ $\gtrsim$ 2.2 for the FSRQ state, $2.0 < \Gamma_{\gamma} < 2.2$ for the transition state, and $\Gamma_{\gamma}$ $\lesssim$ 2.0 for the BL Lac state. We apply the criteria to the photon spectrum evolution of CLB OQ 334 during MJD 58678 - 60387. The evolution is subdivided into five FSRQ states, nine transition states, and four BL Lac states. Moreover, we use the model spectra parameters of each state epoch to test the reliability of subdivided morphological types. The result shows that: (1) the accretion rate parameter is consistent with our earlier research; and (2) there is an increasing trend in the epochs of the BL Lac states, even if there is not an obvious decreasing trend in epochs of the FSRQ states. We issue that strong evidence that a CLB is an especial epoch in the evolution of blazars that could be obtained from the oscillation phenomenon in the CLB evolution.

Circumstellar OH maser lines are useful for studying the dynamics of the circumstellar envelope (CSE) around evolved stars. This study aims to identify CSEs around cold stars bf, which exhibit deviations from the spherical expansion, by comparing the velocity ranges of the OH main lines (1665/1667 MHz) with those of the satellite line (1612 MHz), using a database of circumstellar OH maser sources. We performed this comparison for 377 circumstellar OH maser sources. In addition, using infrared two-color diagrams, we examined the evolutionary stages and infrared properties of objects showing velocity excess (velocity excess means the detection of the main lines outside the velocity range of the satellite line). A periodicity analysis of the WISE light curves was also carried out. As a result of the velocity range comparison, eight circumstellar OH maser sources were found to exhibit velocity excess. The infrared colors of these objects match those of post-AGB stars. Periodic variations were observed in the WISE light curves of five of these eight objects. The results suggest that examining velocity excess of the main lines relative to the satellite line is scientifically significant because mainline masers probe the CSE dynamics over a broader range of evolutionary stages compared to the 22.235 GHz H${_2}$O maser line. Additionally, during the post-AGB phase, the emission regions of the mainline and 22.235 GHz H${_2}$O masers may overlap in a CSE, whereas they originate from different regions during the AGB phase.

We present a preliminary laboratory test of a setup designed to measure Hanbury Brown and Twiss-type intensity correlations from a chaotic light source using five spectral channels simultaneously. After averaging the zero-delay correlation peaks from all channels, we obtain an improvement of the signalto-noise ratio fairly consistent with theory. The goal is to demonstrate the feasibility and scalability of this technique to improve the sensitivity of stellar intensity interferometry using optical telescopes.

Inertial waves in convective regions of stars exhibit topological properties linked to a Chern number of 1. The first of these is a unique, unidirectional, prograde oscillation mode within the cavity, which propagates at arbitrarily low frequencies for moderate azimuthal wavenumbers. The second one are phase singularities around which the phase winds in Fourier space, with winding numbers of $\pm 1$ depending on the hemisphere. Phase winding is a collective effect over waves propagating in all directions that is strongly robust to noise. This suggests a topology-based method for wave detection in noisy observational data.

M-dwarfs show frequent flares and associated coronal mass ejections (CMEs) may significantly impact close-in habitable planets. M-dwarf flares sometimes show red/blue asymmetries in the H$\alpha$ line profile, suggesting prominence eruptions as an early stage of CMEs. However, their high-time-cadence observations are limited. We conducted spectroscopic monitoring observations of the active M-dwarf YZ Canis Minoris with $\sim$1 minute time cadence using the Seimei telescope, simultaneously with the optical photometric observations by Transiting Exoplanet Survey Satellite. We detected 27 H$\alpha$ flares with H$\alpha$ energies ranging from 1.7 $\times$ 10$^{29}$ to 3.8 $\times$ 10$^{32}$ erg and durations from 8 to 319 minutes. Among them, we identified 3 blue asymmetry and 5 red asymmetry events based on criteria using the Bayesian Information Criterion. The maximum velocity of the blue- and red-shifted components ranges from 250 to 450 km s$^{-1}$ and 190 to 400 km s$^{-1}$, respectively. The duration and time evolution show variety, and in particular, we discovered rapid, short-duration blue/red asymmetry events with the duration of 6--8 minutes. Among the 8 blue/red asymmetry events, two blue and one red asymmetry events are interpreted as prominence eruptions because of their fast velocity and time evolution. Based on this interpretation, the lower limit of occurrence frequency of prominence eruptions can be estimated to be $\sim$1.1 events per day. Our discovery of short-duration events suggests that previous studies with low time cadence may have missed these events, potentially leading to an underestimation of the occurrence frequency of prominence eruptions/CMEs.

Understanding the circumgalactic medium (CGM) distribution of galaxies is the key to revealing the dynamical exchange of materials between galaxies and their surroundings. In this work, we use DESI EDR dataset to investigate the cool CGM of galaxies ($0.3<z<1.7$) with stacking the spectra of background QSOs to obtain Mg II absorption of foreground galaxies. The equivalent width of Mg II absorption strongly correlates to stellar mass with EW(Mg II) $\propto M_{*}^{0.5}$ for star-forming galaxies with $\log M_{*}/M_{\odot} < 10$, but is independent with mass for galaxies above this mass. At given stellar mass, EW(Mg II) is larger for galaxies of higher star formation rate with impact parameter less than $50$ kpc, while showing little dependence on galaxy size. By studying the dependence on azimuthal angle, we find EW(Mg II) is strongest at the direction near the minor axis for star-forming galaxies with $\log M_{*}/M_{\odot} < 10.0$, while no dependence on azimuthal angle is seen for luminous red galaxies. This indicates that the outflow associated with star formation enhances the Mg II absorption. However, for galaxies with $\log M_{*}/M_{\odot} > 10.0$, the EW(Mg II) at the minor axis is largely suppressed with respect to low mass galaxies. This suggests that the competing processes, such as stellar feedback and gravity, play a key role in shaping the distribution of outflowing gas.

The Tayler instability (TI) of toroidal magnetic fields is a candidate mechanism for driving turbulence, angular momentum transport, and dynamo action in stellar radiative zones. Recently Skoutnev & Beloborodov (2024) revisited the linear stability analysis of a toroidal magnetic field in a rotating and stably stratified fluid. In this paper, we extend the analysis to include both thermal and compositional stratification, allowing for general application to stars. We formulate an analytical instability criterion for use as a ``toggle switch" in stellar evolution codes. It determines when and where in a star the Tayler-Spruit dynamo may be responsible for angular momentum transport. We implement such a ``toggle switch" in the MESA stellar evolution code and map out the stability of each mode of the TI on a grid of stellar evolution models. In evolved lower mass stars, the TI becomes suppressed in compositionally stratified regions, undermining the Tayler-Spruit dynamo as an explanation of the core-envelope coupling. In higher mass stars, the TI is active throughout the radiative zones, but at different wavenumbers than previously expected.

Determining the spatial curvature $\Omega_K$ of the Universe has long been crucial in cosmology. In practice, this effort is often entangled with assumptions of dark energy. A combination of distance ($D_{\rm M}$, $D_{\rm L}$) and expansion rate ($H(z)$) measurements can break this degeneracy. However, fitting against discrete data points requires parameterizations of distance and expansion rate as functions of redshifts, which often induces cosmological model dependence. In this work, we propose a new dark energy model-independent parameterization of the cosmological comoving radial distance $\chi$. Fitting data combining distance ($D_{\rm M}$, $D_{\rm L}$) and Hubble parameter (or equivalently $D_H$) measurements, we are then able to obtain $\Omega_K$ in a dark energy model-independent manner. We test this parameterization and the associated fitting scheme with mock data generated with a wide range of fiducial dark energy equations of state ($-1.3<w<1.3$), finding that the best-fit $\Omega_K$ is always unbiased. Then we combine SDSS Baryon Acoustic Oscillation (BAO), Pantheon+ sample of Type Ia Supernovae (SNe Ia), and Observational Hubble Data (OHD) to constrain $\Omega_K$. We find a flat universe with $\Omega_K=-0.01\pm 0.09$. Most constraining power is contributed by SDSS BAO, with the BAO-alone constraint $\Omega_K=-0.03 \pm 0.10$. When replacing SDSS BAO with DESI year-one BAO measurement, we obtain $\Omega_K=0.06 \pm 0.08$. With the full DESI BAO data alone, we forecast $\sigma(\Omega_K)\sim 0.03$. Our result verifies the flatness of the universe free of dark energy modeling, and the proposed parameterization would be useful for future investigation of $\Omega_K$ and other parameters of interest, such as the horizon radius.

In astrophysical environments such as core-collapse supernovae (CCSNe) and binary neutron star mergers (BNSMs), neutrinos potentially experience substantial flavor mixing due to the refractive effects of neutrino self-interactions. Determining the survival probability of neutrinos in asymptotic states is paramount to incorporating flavor conversions' effects in the theoretical modeling of CCSN and BNSM. Some phenomenological schemes have shown good performance in approximating asymptotic states of fast neutrino-flavor conversions (FFCs), known as one of the collective neutrino oscillation modes induced by neutrino self-interactions. However, a recent study showed that they would yield qualitatively different asymptotic states of FFC if the neutrino number is forced to evolve. It is not yet fully understood why the canonical phenomenological models fail to predict asymptotic states. In this paper, we perform detailed investigations through numerical simulations and then provide an intuitive explanation with a quasi-homogeneous analysis. Based on the analysis, we propose a new phenomenological model, in which the quasi-steady evolution of FFCs is analytically determined. The model also allows us to express the convolution term of spatial wave number as a concise form, which corresponds to useful information on analyses for the non-linear feedback from small-scale flavor conversions to large-scale ones. Our model yields excellent agreement with numerical simulations, which lends support to our interpretation.

In this paper we present three different applications, based on deep learning methodologies, that we are developing to support the scientific analysis conducted within the ASKAP-EMU and MeerKAT radio surveys. One employs instance segmentation frameworks to detect compact and extended radio sources and imaging artefacts from radio continuum images. Another application uses gradient boosting decision trees and convolutional neural networks to classify compact sources into different astronomical classes using combined radio and infrared multi-band images. Finally, we discuss how self-supervised learning can be used to obtain valuable radio data representations for source detection, and classification studies.

The jet may compress an interstellar medium and soft X-ray may be absorbed outside the jet by a compressed interstellar medium and thermal emission heated by a shock is expected inside the jet. I analyze X-ray image of surroundings of the jet of M87 by Chandra. There is an original work for image analysis. In this time, I select data for image by removal of pile up event completely. I confirm a dip in soft X-ray outside the jet between knot E and knot F with 260 ks archival data. I analyze X-ray energy spectra for knot E and knot F with 400 ks archival data by Chandra in order to search for thermal emission. However, X-ray energy spectra for knot E and knot F are well described with synchrtron emission.

In this paper, we developed a spectral emulator based on the Mapping Nearby Galaxies at Apache Point Observatory Stellar Library (MaStar) and a grouping optimization strategy to estimate effective temperature (T_eff), surface gravity (log g), metallicity ([Fe/H]) and the abundance of alpha elements with respect to iron ([alpha/Fe]) for O-M-type stars within the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution spectra. The primary aim is to use a rapid spectral-fitting method, specifically the spectral emulator with the grouping optimization strategy, to create a comprehensive catalog for stars of all types within LAMOST, addressing the shortcomings in parameter estimations for both cold and hot stars present in the official LAMOST AFGKM-type catalog. This effort is part of our series of studies dedicated to establishing an empirical spectral library for LAMOST. Experimental results demonstrate that our method is effectively applicable to parameter prediction for LAMOST, with the single-machine processing time within $70$ hr. We observed that the internal error dispersions for T_eff, log g, [Fe/H], and [alpha/Fe] across different spectral types lie within the ranges of $15-594$ K, $0.03-0.27$ dex, $0.02-0.10$ dex, and $0.01-0.04$ dex, respectively, indicating a good consistency. A comparative analysis with external data highlighted deficiencies in the official LAMOST catalog and issues with MaStar parameters, as well as potential limitations of our method in processing spectra with strong emission lines and bad pixels. The derived atmospheric parameters as a part of this work are available at this https URL .

In comparing the two alternative explosion mechanisms of core-collapse supernovae (CCSNe), I examine recent three-dimensional (3D) hydrodynamical simulations of CCSNe in the frame of the delayed-neutrino explosion mechanism (neutrino mechanism) and argue that these valuable simulations show that neutrino heating can supply a non-negligible fraction of the explosion energy but not the observed energies, hence cannot be the primary explosion mechanism. In addition to the energy crisis, the neutrino mechanism predicts many failed supernovae that are not observed. The most challenging issue of the neutrino mechanism is that it cannot account for point-symmetric morphologies of CCSN remnants, many of which were identified in 2024. These contradictions with observations imply that the neutrino mechanism cannot be the primary explosion mechanism of CCSNe. The alternative jittering-jets explosion mechanism (JJEM) seems to be the primary explosion mechanism of CCSNe; neutrino heating boosts the energy of the jittering jets. Even if some simulations show explosions of stellar models (but usually with energies below observed), it does not mean that the neutrino mechanism is the explosion mechanism. Jittering jets, which simulations do not include, can explode the core before the neutrino heating process does. Morphological signatures of jets in many CCSN remnants suggest that jittering jets are the primary driving mechanism, as expected by the JJEM.

The Cosmological Principle is a cornerstone of the standard model of cosmology and shapes how we view the Universe and our place within it. It is imperative, then, to devise multiple observational tests which can identify and quantify possible violations of this foundational principle. One possible method of probing large-scale anisotropies involves the use of weak gravitational lensing. We revisit this approach in order to analyse the imprint of late-time anisotropic expansion on cosmic shear. We show that the cross-correlation of shear $E$- and $B$-modes on large scales can be used to constrain the magnitude (and possibly direction) of anisotropic expansion. We estimate the signal to noise for multipoles $10\lesssim \ell\lesssim 100$ that is achievable by a Euclid-like survey. Our findings suggest that such a survey could detect the $E$-$B$ signal for reasonable values of the late-time anisotropy parameter.

The nature of dark matter is poorly constrained on subgalactic scales. Alternative models to cold dark matter, such as warm dark matter or self-interacting dark matter, could produce very different dark haloes on these scales. One of the few known dark haloes smaller than a galaxy was discovered in the triple source plane strong lens system J0946+1006. Previous studies have found that this structure is much more concentrated than expected in $\Lambda$CDM, but have assumed the dark halo is at the same redshift as the main deflector ($z_{\rm main}=0.222$). In this paper, we fit for the redshift of this dark halo. We reconstruct the first two sources in the system using a forward modelling approach, allowing for additional complexity from multipole perturbations. We find that the perturber redshift is $z_{\rm halo} = {0.207}^{+0.019}_{-0.019}$, and lower bounds on the evidence strongly prefer a subhalo over a line-of-sight structure. Whilst modelling both background sources does not improve constraints on the redshift of the subhalo, it breaks important degeneracies affecting the reconstruction of multipole perturbations. We find that the subhalo is a more than $5\sigma$ outlier from the $\Lambda$CDM $v_{\rm max}$-$r_{\rm max}$ relation and has a steep profile with an average slope of $\gamma_{\rm 2D} = {-1.81}^{+0.15}_{-0.11}$ for radii between $0.75-1.25$ kpc. This steep slope might indicate dark matter self-interactions causing the subhalo to undergo gravothermal collapse; such collapsed haloes are expected to have $\gamma_{\rm 2D} \approx -2$.

The phase evolution of gravitational waves (GWs) can be modulated by the astrophysical environment surrounding the source, which provides a probe for the origin of individual binary black holes (BBHs) using GWs alone. We here study the evolving phase of the GW waveform derived from a large set of simulations of BBH mergers forming in dense stellar clusters through binary-single interactions. We uncover that a well-defined fraction of the assembled eccentric GW sources will have a notable GW phase shift induced by the remaining third object. The magnitude of the GW phase shift often exceeds conservative analytical estimates due to strong 3-body interactions, which occasionally results in GW sources with clearly shifted and perturbed GW waveforms. This opens up promising opportunities for current and future GW detectors, as observing such a phase shift can identify the formation environment of a BBH, as well as help to characterise the local properties of its surrounding environment.

Axion-photon oscillation effect provides a possible explanation for the presence of very-high-energy (VHE) $\gamma$-ray signals from distant sources. In this work, we propose a model-dependent method to select possible sources that may give sufficient constraints on the axion parameters. We investigate such effect in the spectra of active galactic nuclei (AGN) B2 2234+28A and 3C 454.3 based on data obtained from Fermi Large Area Telescope (Fermi-LAT) and MAGIC U.L. We utilize the Markov Chain Monte Carlo method to fit the axion parameters, yielding a result of $g_{a\gamma}=3.05^{+0.51}_{-0.31} \times 10^{-11}$ GeV$^{-1}$ for the axion-photon coupling strength and $m_{a}=5.25^{+2.35}_{-2.65} \times 10^{-8} $ eV for the axion mass. We also perform 95\% confidence level (CL) constraints to set an upper limit for $g_{a\gamma}$.

$F(R)$ models for dark energy generally exhibit a weak curvature singularity, which can be cured by adding an $R^2$ term. This correction allows for a unified description of primordial and late-time accelerated expansions. However, most existing models struggle to achieve this, as they become unstable over certain negative ranges of the Ricci scalar, where either the first or second derivative of $F(R)$ turns negative. These instabilities may disrupt the post-inflationary evolution when the Ricci scalar oscillates about the vacuum state after the $R^2$ inflation. In this work, we introduce a new model-building to guarantee global stability, i.e., the first and second derivatives are positive for all real Ricci scalars. By extending the idea from Appleby and Battye, we demonstrate that viable models can be constructed by imposing a positive, bounded first derivative of $F(R)$ with a sigmoid shape. As examples, we first reformulate and generalize the original Appleby-Battye model. Then, we propose a new dark energy model, which successfully explains the acceleration of cosmic expansion and passes local gravity tests.

We present JWST/NIRSpec integral field spectroscopic (IFS) observations of the \Lyalpha emitter CR7 at z ~ 6.6, observed as part of the GA-NIFS program. Using low-resolution PRISM (R ~ 100) data, we confirm a bright \Lyalpha emitter, and a diffuse \Lyalpha halo extending up to 3 kpc from the peak of ionized emission, both of them associated to the most massive, UV bright galaxy in the system (CR7-A). We confirm the presence of two additional UV-bright satellites (CR7-B and CR7-C) detected at projected distances of 6.4 and 5.2 kpc from the primary source. We perform SED fitting of the low-resolution data and revealed an inverted star formation history between two satellites at early epochs and a spatially resolved anti-correlation of the gas-phase metallicity and the star formation rate density, likely driven by the gas exchange among the satellites, favouring the merger scenario for CR7. From the high-resolution G395H (R ~ 2700) data, we discover at least three additional companions mainly traced by the \OIIIL emission line, although they are not detected in continuum. We disentangle the kinematics of the system and reveal extended ionised emission linking the main galaxy and the satellites. We spatially resolve the \OIIIL, \OIII[4363], and \Hgamma emission lines and use a diagnostic diagram tailored to high-z systems to reveal tentative evidence of AGN ionisation across the main galaxy (CR7-A) and the N-E companion (CR7-B). Moreover, we detect an unresolved blue-shifted outflow from one of the satellites and present first evidence for a redshifted outflow from the main galaxy. Finally, we compute resolved electron temperature (T$_e \sim 1.6 \times 10^4$ K) and metallicity maps (log(Z/\zsun) from --0.8 to --0.5), and provide insights on how the physical properties of the system evolved at earlier epochs.

Jets launched by supermassive black holes transport relativistic leptons, magnetic fields, and atomic nuclei from the centres of galaxies to their outskirts and beyond. These outflows embody the most energetic pathway by which galaxies respond to their Cosmic Web environment. Studying black hole feedback is an astrophysical frontier, providing insights on star formation, galaxy cluster stability, and the origin of cosmic rays, magnetism, and heavy elements throughout the Universe. This feedback's cosmological importance is ultimately bounded by the reach of black hole jets, and could be sweeping if jets travel far at early epochs. Here we present the joint LOFAR-uGMRT-Keck discovery of a black hole jet pair extending over $7$ megaparsecs -- the largest galaxy-made structure ever found. The outflow, seen $7.5$ gigayears into the past, spans two-thirds of a typical cosmic void radius, thus penetrating voids at ${\sim}95\%$ probability. This system demonstrates that jets can avoid destruction by magnetohydrodynamical instabilities over cosmological distances, even at epochs when the Universe was 15 to 7 times denser than it is today. Whereas previous record-breaking outflows were powered by radiatively inefficient active galactic nuclei, this outflow is powered by a radiatively efficient active galactic nucleus, a type common at early epochs. If, as implied, a population of early void-penetrating outflows existed, then black hole jets could have overwritten the fields from primordial magnetogenesis. This outflow shows that energy transport from supermassive black holes operates on scales of the Cosmic Web and raises the possibility that cosmic rays and magnetism in the intergalactic medium have a non-local, cross-void origin.

The Galactica simulation database is a platform designed to assist computational astrophysicists with their open science approach based on FAIR (Findable, Accessible, Interoperable, Reusable) principles. It offers the means to publish their numerical simulation projects, whatever their field of application or research theme and provides access to reduced datasets and object catalogs online. The application implements the Simulation Datamodel IVOA standard. To provide the scientific community indirect access to raw simulation data, Galactica can generate, on an "on-demand" basis, custom high-level data products to meet specific user requirements. These data products, accessible through online WebServices, are produced remotely from the raw simulation datasets. To that end, the Galactica central web application communicates with a high-scalability ecosystem of data-processing servers called Terminus by means of an industry-proven asynchronous task management system. Each Terminus node, hosted in a research institute, a regional or national supercomputing facility, contributes to the ecosystem by providing both the storage and the computational resources required to store the massive simulation datasets and post-process them to create the data products requested on Galactica, hence guaranteeing fine-grained sovereignty over data and resources. This distributed architecture is very versatile, it can be interfaced with any kind of data-processing software, written in any language, handling raw data produced by every type of simulation code used in the field of computational astrophysics. Its generality and versatility, together with its excellent scalability makes it a powerful tool for the scientific community to disseminate numerical models in astrophysics in the exascale era.

By observing binary black hole (BBH) mergers out to the edge of the Universe, next-generation (XG) ground-based gravitational-wave (GW) detectors like Cosmic Explorer and Einstein Telescope will map the BBH merger rate across all of cosmic history. This merger rate traces the formation rate of their progenitor stars convolved with a delay time distribution. Given theoretically-motivated priors on the delay time distribution, we show how XG observations can measure the BBH progenitor formation rate, probing the star formation rate (SFR) up to $z > 15$. However, the progenitor formation rate does not directly give a measurement of the SFR, but rather a combination of the SFR and its metallicity distribution as a function of redshift. Fortunately, the metallicity-dependence of BBH formation likely varies as a function of BBH mass and/or formation channel. We find that if different BBH subpopulations with distinct metallicity biases can be identified, comparing their rates as a function of redshift yields a simultaneous measurement of the SFR and its metallicity distribution. Given optimistic theoretical priors and one year of observation, this may provide a $\sim10\%$ measurement of the SFR at its peak and a 0.2 dex (0.7 dex) measurement of the median metallicity out to $z = 10$ ($z = 15$) at 90\% credibility, although the uncertainties scale with theoretical uncertainties on BBH delay times and formation efficiencies.

The multiple impact hypothesis proposes that the Moon formed through a series of smaller collisions, rather than a single giant impact. This study advances our understanding of this hypothesis, as well as moon collisions in other contexts, by exploring the implications of these smaller impacts, employing a novel methodological approach that combines self-consistent initial conditions, hybrid hydrodynamic/N-body simulations, and the incorporation of material strength. Our findings challenge the conventional assumption of perfect mergers in previous models, revealing a spectrum of collision outcomes including partial accretion and mass loss. These outcomes are sensitive to collision parameters and Earth's tidal influence, underscoring the complex dynamics of lunar accretion. Importantly, we demonstrate that incorporating material strength is important for accurately simulating moonlet-sized impacts. This inclusion significantly affects fragmentation, tidal disruption, and the amount of material ejected or accreted onto Earth, ultimately impacting the Moon's growth trajectory. By accurately modeling diverse collision outcomes, our hybrid approach provides a powerful new framework for understanding the Moon's formation. We show that most collisions (~90%) do not significantly erode the largest moonlet, supporting the feasibility of lunar growth through accretion. Moreover, we revise previous estimates of satellite disruption, suggesting a higher survival rate and further bolstering the multiple-impact scenario.

Changing-look active galactic nuclei (CLAGNs) show the appearance and disappearance of broad emission lines in their UV/optical spectra on timescales of months to decades. Here, we investigate how CL transitions depend on several AGN parameters such as accretion rate, obscuration properties and black hole mass. We study a sample of 20 nearby optically-identified CLAGNs from the BAT AGN Spectroscopic Survey (BASS), using quasi-simultaneous optical and X-ray observations taken in the last $\sim 40$ years. We find that for all CLAGNs, the transition is accompanied by a change in Eddington ratio. The CL transitions are not associated with changes in the obscuration properties of the AGN. CLAGNs are found to have a median Eddington ratio lower than the AGNs in the BASS sample in which CL transitions were not detected. The median of the transition Eddington ratio (Eddington ratio at which AGN changes its state) is found to be $\sim 0.01$ for type 1 $\leftrightarrow$ 1.8/1.9/2 transition, which is consistent with the hard $\leftrightarrow$ soft state transition in black hole X-ray binaries. Most CL events are constrained to occur within 3-4 years, which is considerably shorter than the expected viscous timescale in AGN accretion disk. The transitions of the optical CLAGNs studied here are likely associated to state changes in the accretion flow, possibly driven by disk-instability.

Stellar-driven galactic winds are multiphase outflows of energy and matter connecting the interstellar and circumgalactic media (CGM) with the intergalactic medium. Galactic winds contain a hot and diffuse phase detected in X-rays, and a cold and dense phase detected via emission and absorption lines from the ions populating the outflow. The ion production within galactic winds largely depends on the background UV radiation field produced by star formation, and this in turn depends on the age of the starburst, the gas metallicity, the proximity of the outflowing gas to the central star-forming regions. Our study probes the influence of the proximity of wind-cloud systems to the UV background source, and the effects of magnetic fields on the N V ion production through the analysis of synthetic column densities and spectral lines. We utilise magnetohydrodynamical simulations to study weakly-magnetised wind-cloud systems, and extract synthetic spectral lines with Trident and yt. Our simulations indicate that magnetic fields transverse to the wind have a shielding effect on dense gas, producing broader N V absorption lines. Also, a weak (distant) UV background produces N V only in the outer cloud layers with no spectral signature, while a strong (nearby) UV background produces it in the cloud core with a narrow spectral line. Overall, transverse magnetic fields and a UV radiation at 50 kpc produce the stronger N V spectral lines.

The rovibrational level populations, and subsequent emission in various astrophysical environments, is driven by inelastic collision processes. The available rovibrational rate coefficients for water have been calculated using a number of approximations. We present a numerically exact calculation for the rovibrational quenching for all water vibrational modes due to collisions with atomic hydrogen. The scattering theory implements a quantum close-coupling (CC) method on a high level ab initio six-dimensional (6D) potential energy surface (PES). Total rovibrational quenching cross sections for excited bending levels were compared with earlier results on a 4D PES with the rigid-bender close-coupling (RBCC) approximation. General agreement between 6D-CC and 4D-RBCC calculations are found, but differences are evident including the energy and amplitude of low-energy orbiting resonances. Quenching cross sections from the symmetric and asymmetric stretch modes are provided for the first time. The current 6D-CC calculation provides accurate inelastic data needed for astrophysical modeling.

In the past decade, an asymmetry in the large-scale distribution of galaxy spin directions has been observed in data from all relevant digital sky surveys, all showing a higher number of galaxies rotating in the opposite direction relative to the Milky Way as observed from Earth. Additionally, JWST deep fields have shown that the asymmetry is clear and obvious, and can be sensed even by the naked human eye. These experiments were performed using two separate statistical methods: standard binomial distribution and simple $\chi^2$ statistics. Stiskalek \& Desmond (2024) suggested that the asymmetry in the distribution of galaxy spin directions is due to the use of binomial or $\chi^2$ statistics. Instead, they developed a new complex ad-hoc statistical method that shows random distribution in galaxy spin directions, and specifically in data from HSC. Source code for the method was also made available. The primary downside of the new method is that it is not able to identify asymmetry in the distribution of galaxy spin directions. Even when the new method is provided with synthetic data with extreme and obvious asymmetry, it still reports a null-hypothesis Universe with random distribution. That shows empirically that the method cannot sense asymmetry in the distribution of the directions of rotation of galaxies. While this further concludes that the distribution of galaxy spin direction as observed from Earth is not symmetric, it is not necessarily an indication of an anomaly in the large-scale structure. The excessive number of galaxies that rotate in the opposite direction relative to the Milky Way can also be driven by the internal structure of galaxies and the physics of galaxy rotation. The phenomenon can be related to other puzzling anomalies such the Ho tension. Data are publicly available, and no code is needed to reproduce the results since only conventional statistics is used.

We have developed a new regression technique, the maximum likelihood (ML)-based method and its variant, the KS-test based method, designed to obtain unbiased regression results from typical astronomical data. A normalizing flow model is employed to automatically estimate the unobservable intrinsic distribution of the independent variable as well as the unobservable correlation between uncertainty level and intrinsic value of both independent and dependent variables from the observed data points in a variational inference based empirical Bayes approach. By incorporating these estimated distributions, our method comprehensively accounts for the uncertainties associated with both independent and dependent variables. Our test on both mock data and real astronomical data from PHANGS-ALMA and PHANGS-JWST demonstrates that both the ML based method and the KS-test based method significantly outperform the existing widely-used methods, particularly in cases of low signal-to-noise ratios. The KS-test based method exhibits remarkable robustness against deviations from underlying assumptions, complex intrinsic distributions, varying correlations between uncertainty levels and intrinsic values, inaccuracies in uncertainty estimations, outliers, and saturation effects. We recommend the KS-test based method as the preferred choice for general applications, while the ML based method is suggested for small samples with sizes of $N < 100$. A GPU-compatible Python implementation of our methods, nicknamed ``raddest'', will be made publicly available upon acceptance of this paper.

Super-Eddington accretion onto stellar-mass compact objects powers fast outflows in Ultra-Luminous X-ray sources (ULXs). Such outflows, which can reach mildly relativistic velocities, are often observed forming bubble structures. Wind bubbles are expected to develop strong wind termination shocks, sites of great interest for diffusive shock acceleration. We develop a model of diffusive shock acceleration in the wind bubbles powered by ULXs. We find that the maximum energy in these objects can easily reach the PeV range, promoting ULX winds as a new class of PeVatrons. We specialize our model in the context of the Galactic source SS 433 and show that high-energy protons in the bubble might explain the highest energy photons (>100 TeV) and their morphology recently observed by LHAASO. We discuss the detectability of such a source in neutrinos and we analyze the possible radio counterpart of ULXs focusing on the case of W50, the nebula surrounding SS 433. We finally discuss the possible contribution of Galactic ULXs to the cosmic-ray flux at the knee concluding that their role might be substantial.

It was hypothesized in the literature that some physical parameters may be time-evolving and the astrophysical data can serve as a probe. Recently, James Webb Space Telescope (JWST) have released its early observations. In this work, we select the JWST spectroscopic observations of the high redshift ($z>7.1$) galaxies with strong [OIII] ($\lambda=4959$ Å\,and $5007$ Å\,in the rest frame) emission lines to constraint the evolution of the fine structure constant ($\alpha$). With the spectra from two galaxies at redshifts of $7.19$ and $8.47$, the deviation of $\alpha$ to its fiducial value is found to be as small as $0.44^{+8.4+1.7}_{-8.3-1.7} \times 10^{-4}$ and $-10.0^{+18+1.5}_{-18-1.5} \times 10^{-4}$, respectively (the first error is statistical and the latter is systematic). The combination of our results with the previous data reveals that $\frac{1}{\alpha} \frac{d \alpha}{dt} = 0.30^{+4.5}_{-4.5} \times 10^{-17}~{\rm yr^{-1}}$. Clearly, there is no evidence for a cosmic evolution of $\alpha$. The prospect of further constraining the time evolution of $\alpha$ is also discussed. The scalar field of dark energy is hypothesized to drive the acceleration of the universe's expansion through an interaction with the electromagnetic field. By integrating the observational data of the fine-structure constant variation, $\frac{\Delta\alpha}{\alpha}(z)$, we have established a stringent upper limit on the coupling strength between dark energy and electromagnetism. Our analysis yields $\zeta \leq 3.92 \times 10^{-7}$ at the 95\% confidence level, representing the most stringent bound to date.

The Sun's chromosphere is a critical region to understand when considering energy and mass deposition into the transition region and corona, but many of the smaller, faster events which transport a portion of this mass and energy are still difficult to observe, identify and model. Solar Spicules are small, spike-like events in the solar chromosphere that have the potential to transfer energy and mass to the transition region, but whose energetic origins are still being researched. Chromospheric spicule activity on-disk can be identified by observing temporary excursions in the red and blue wings of chromospheric emission lines. Researchers have demonstrated this in Hydrogen~Alpha (H{\alpha}, 6563 Å), Ca II (8542 Å, k 3934 Å), Mg II (h 2803 Å, k 2796 Å), and Si IV (1394 Å, 1405 Å) spectral observations, with the vast majority of identification efforts focused on lower chromospheric observations of H$\alpha$ and Ca II. Because any spicules which deposit mass and energy into the transition region must necessarily pass through the upper chromosphere, observations from this region such as Mg II or Hydrogen Lyman Alpha (Ly$\alpha$ 1216 Å) in enough quantity to perform proper statistics will be critical to fully characterizing spicules' impact on mass and energy transfer in the Sun. This research proposes a definition with numerical limits for how spicules appear in Mg II wavelengths, tunes an algorithm for automatically detecting spicules in Mg II spectral observations, and uses K Means Clustering to identify and display the full range of spicule spectrum shapes. This work will help allow statistical studies on spicules in the upper chromosphere to be as thorough as those of the lower chromosphere, allowing researchers to better understand the physical nature of spicules and their role in energy transfer and deposition in the solar atmosphere.

Rotational irregularities are one of the prominent observational features that most pulsars exhibit. These glitches, which are sudden increases in spin angular velocity, remains an open problem. In this study, we have investigated the potential role of nontrivial topological defects, specifically in the form of Nambu-goto-type CSs, and its connection to spin irregularities. Such CSs which are one-dimensional topological defects may be formed during various symmetry-breaking and phase transition scenarios and can interact with the neutron stars. In this work, we see that the appearance of such topological defects trapped within the core can lead to the coupling of the string tension with the angular velocity, leading to the abrupt rotational changes observed as pulsar glitches. We have further studied how these coupling may generate detectable gravitational waves as a mixture of continuous and burst signals. The evolution of cusps of CSs trapped within neutron stars and the neutron star's mass quadruple moment change due to rotation could produce distinctive gravitational wave signatures, well within the noise cutoff of advLIGO. Our study highlights a potential connection between topological defects, pulsar glitches, and gravitational wave emissions, offering a possible avenue for observationally testing the presence of CSs and their astrophysical effects.

The Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP) reveals complex nested morphological and kinematic features of molecular outflows through the CO (J = 2 - 1) and SiO (J = 5 - 4) emission. We characterize the jet and outflow kinematics of the ALMASOP sample in four representative sources (HOPS 10, 315, 358, and G203.21-11.20W2) through channel maps and position-velocity diagrams (PVDs) parallel and transverse to the outflow axes. The combined CO and SiO emission exhibits the coexistence of the conventional extremely-high-velocity (EHV) jets and shell-like low-velocity (LV) cavity walls and new features. More complex, nested bubble-like and filamentary structures in the images and channel maps, triangle-shaped regions near the base of the parallel PVDs, and regions composed of rhombus/oval shapes in the transverse PVDs, are also evident. Such features find natural explanations within the bubble structure of the unified model of jet, wind, and ambient medium. The reverse shock cavity is revealed on the PVD base regions, and other features naturally arise within the dynamic postshock region of magnetic interaction. The finer nested shells observed within the compressed wind region reveal previously unnoticed shocked emission between the jet and the conventional large cavity walls. These pseudopulse-produced filamentary features connect to the jet-like knotty blobs, creating an impression of episodicity in mass ejection. SiO emission is enhanced downstream of the reverse shock boundary, with jet-like excitation conditions. Combined, these observed features reveal the extended structures induced by the magnetic interplay between a jet-bearing magnetized wide-angle wind and its ambient magnetized surrounding medium.

Context. Small bodies in Earth-like orbits, the Arjunas, are good targets for scientific exploration and mining studies as they enable low-cost missions. The subset of such objects that experience recurrent temporarily captured flyby or orbiter events, also called mini-moon episodes, are among the best ranked in terms of accessibility. Only a handful of objects are known to have engaged in such a dynamical behavior. Finding and characterizing more of them may help to expand scientific and commercial research activities in space over the next few decades. Asteroid 2024 PT5 was found recently and belongs to this group of interesting objects. Aims. Here, we investigate the orbital context of 2024 PT5, and its spectral and rotational properties. Methods. We studied the short-term orbital evolution of 2024 PT5 using direct N-body simulations. We identified its spectral class from the visible reflectance spectrum and used photometric observations to derive its rotational properties. Observational data were obtained with the OSIRIS camera spectrograph at the 10.4 m Gran Telescopio Canarias and the Two-meter Twin Telescope. Results. Asteroid 2024 PT5 experiences recurrent co-orbital engagements and mini-moon events of the temporarily captured flyby type. Its visible spectrum is consistent with that of an Sv-type asteroid or perhaps lunar ejecta. Its rotational period could be less than or close to 1 h. Conclusions. The discovery of 2024 PT5 confirms that temporarily captured flybys are relatively frequent and involve objects larger than a few meters, suitable as accessible targets for scientific research activities and commercial mining ventures in space.

We report the discovery of a close binary J0606+2132 (Gaia DR3 3423365496448406272) with $P_{\rm obs}=2.77$ days containing a possible massive white dwarf or a neutron star using the LAMOST spectroscopic data. By a joint fitting of the radial velocity from LAMOST and the light curve from TESS, we derived a circular Keplerian orbit with an inclination of $i=$81.31$^{\circ}$$^{+6.26^{\circ}}_{-7.85^{\circ}}$, which is consistent with that derived from $v{\rm sin}I$. Together with the mass of the visible star, we derived the mass of the invisible object to be 1.34$^{+0.35}_{-0.40} M_{\odot}$. Spectral disentangling with the LAMOST medium-resolution spectra shows no absorption feature from an additional component, suggesting the presence of a compact object. No X-ray or radio pulsed signal is detected from ROSAT and FAST archive observations. J0606+2132 could evolve into either a Type Ia supernova or a neutron star through accretion-induced collapse if it is a white dwarf, or into an intermediate-mass X-ray binary if it is a neutron star.

Supermassive black holes (SMBHs) are observed in diverse galaxy populations across time yet a clear understanding of how they coevolve with their hosts has not been reached. Physically-motivated models of SMBH accretion and feedback vary widely between galaxy formation simulations due to the difficulty of modeling the range of scales important for galactic and SMBH processes. Here we use observational data to build an empirical model for SMBH growth. We apply observed specific accretion rate probability distributions as a function of galaxy star formation rate between $z = 0-2$ to the UniverseMachine galaxy formation model to determine SMBH accretion rates based on galaxy properties. We use observed $z = 0$ SMBH-stellar mass relations for the quiescent and star-forming populations to provide the local boundary conditions for SMBH growth histories. We then track the coevolutionary histories of galaxy stellar mass and their SMBHs backwards in time to $z = 2$. We find that the most massive SMBHs at $z = 0$ have grown very little of their total mass between $z = 0-2$, indicating early SMBH mass assembly for these systems. Conversely, lower mass SMBHs at $z = 0$ assembled their mass gradually across $z = 0-2$. This results in substantial evolution of the SMBH-stellar mass relation, shifting to higher normalization and shallower slope with increasing redshift. We find that the substantial scatter observed in the $z = 0$ SMBH-stellar mass relation results in the diversity of growth pathways found in our model, with some galaxies assembling their stellar mass before their SMBHs and others doing the opposite.

While machine-learned models are now routinely employed to facilitate astronomical inquiry, model inputs tend to be limited to a primary data source (namely images or time series) and, in the more advanced approaches, some metadata. Yet with the growing use of wide-field, multiplexed observational resources, individual sources of interest often have a broad range of observational modes available. Here we construct an astronomical multimodal dataset and propose AstroM$^3$, a self-supervised pre-training approach that enables a model to learn from multiple modalities simultaneously. Specifically, we extend the CLIP (Contrastive Language-Image Pretraining) model to a trimodal setting, allowing the integration of time-series photometry data, spectra, and astrophysical metadata. In a fine-tuning supervised setting, our results demonstrate that CLIP pre-training improves classification performance for time-series photometry, where accuracy increases from 84.6% to 91.5%. Furthermore, CLIP boosts classification accuracy by up to 12.6% when the availability of labeled data is limited, showing the effectiveness of leveraging larger corpora of unlabeled data. In addition to fine-tuned classification, we can use the trained model in other downstream tasks that are not explicitly contemplated during the construction of the self-supervised model. In particular we show the efficacy of using the learned embeddings for misclassifications identification, similarity search, and anomaly detection. One surprising highlight is the "rediscovery" of Mira subtypes and two Rotational variable subclasses using manifold learning and dimension reduction algorithm. To our knowledge this is the first construction of an $n>2$ mode model in astronomy. Extensions to $n>3$ modes is naturally anticipated with this approach.

Stochastic inflation, together with the $\Delta N$ formalism, provides a powerful tool for estimating the large-scale behaviour of primordial fluctuations. We construct a numerical code to capture the non-perturbative statistics of such fluctuations and test our code to obtain the exponential non-Gaussian tail of the curvature perturbations. We provide a numerical algorithm to compute the non-perturbative curvature power spectrum and apply it to slow-roll (SR) and ultra-slow-roll (USR) single-field models of inflation. For the USR case, we successfully reproduce the peak in the power spectrum, which agrees to a certain accuracy with the perturbative power spectrum. We highlight some important differences between non-perturbative and perturbative approaches that may suggest the inconsistency of the $\Delta N$ formalism at the transition stages between attractor and non-attractor regimes.

This paper presents the GRBSN webtool, a public facing application which hosts the most complete list of GRB-SN associations to date. In contrast to other repositories of supernova or gamma-ray burst data, this tool brings together all of the information required to study a GRB-SN association. GRBSN allows users to view and interact with plots of the data; search and filter the whole database; and download all multi-wavelength data related to a GRB-SN association, including radio, X-ray, optical/NIR photometric and spectroscopic data. The tool is fully open source and is hosted on a public GitHub repository, meaning users can upload their own data, flag missing data and suggest improvements. As the number of confirmed GRB-SN associations increases, the webtool will provide a robust framework in which to catalogue these associations and their associated data. The web application is freely available and publicly accessible at this https URL.

We show that the anisotropies in the spectrum of gravitational waves induced by scalar modes after the end of inflation in canonical, single-field models are completely determined by the tilt of the scalar and tensor power spectra. The latter contains information about anisotropies produced due to the propagation of the tensor modes in an inhomogeneous Universe, whereas the former represents the anisotropies generated at the time of production and arise only when non-Gaussian corrections to the angular power spectrum are considered. Our proof takes into account all scalar interactions in the cubic inflaton Lagrangian.

Quantifying the global reach of planetarium dome shows presents significant challenges due to the lack of standardized viewership tracking mechanisms across diverse planetarium venues. We present an analysis of the global impact of dome shows, presenting data regarding four documentary films from a single visualization lab. Specifically, we designed and administered a viewership survey of four long-running shows that contained cinematic scientific visualizations. Reported survey data shows that between 1.2 - 2.6 million people have viewed these four films across the 68 responding planetariums (mean: 1.9 million). When we include estimates and extrapolate for the 315 planetariums that licensed these shows, we arrive at an estimate of 16.5 - 24.1 million people having seen these films (mean: 20.3 million).

Terrestrial experiments that use electrons in Earth as a spin-polarized source have been demonstrated to provide strong bounds on exotic long-range spin-spin and spin-velocity interactions. These bounds constrain the coupling strength of many proposed ultralight bosonic dark-matter candidates. Recently, it was pointed out that a monopole-dipole coupling between the Sun and the spin-polarized electrons of Earth would result in a modification of the precession of the perihelion of Earth. Using an estimate for the net spin-polarization of Earth and experimental bounds on Earth's perihelion precession, interesting constraints were placed on the magnitude of this monopole-dipole coupling. Here we investigate the spin associated with Earth's electrons. We find that there are about $6 \times 10^{41}$ spin-polarized electrons in the mantle and crust of Earth oriented anti-parallel to their local magnetic field. However, when integrated over any spherically-symmetric Earth model, we find that the vector sum of these spins is zero. In order to establish a lower bound on the magnitude of the net spin along Earth's rotation axis we have investigated three of the largest breakdowns of Earth's spherical symmetry: the large low shear-velocity provinces of the mantle, the crustal composition, and the oblate spheroid of Earth. From these investigations we conclude that there are at least $5 \times 10^{38}$ spin-polarized electrons aligned anti-parallel to Earth's rotation axis. This analysis suggests that the bounds on the monopole-dipole coupling that were extracted from Earth's perihelion precession need to be relaxed by a factor of about 2000.

Next-generation low-frequency interferometers are expected to detect binary systems near supermassive black holes, where tidal effects can alter significantly the motion of the binary. This motivates a broader investigation of how external gravitational fields influence the dynamics of physical systems. In this work, we consider a charged black hole binary system subject to a gravitational tide. We first construct a stationary gravitational tide acting on a dyonic Reissner-Nordström black hole and, focusing on the extreme mass-ratio limit, we analyze the motion of a test particle. By calculating the secular Hamiltonian of the test particle, we obtain the ISCO and light ring tidal shifts in terms of explicit functions of the parameters of the binary. Our results show that tidal corrections are suppressed as the charge of the black hole increases, but they persist in the extremal limit yielding a finite contribution. This work paves the way towards studying tidal effects on other charged systems, such as topological stars.

In the present work, evaporation of a black hole immersed in a de Sitter environment is considered. Vaidya-de Sitter spacetime is used to model the process in a scenario of accelerated expansion of the Universe. The role of observers is highlighted in the development and Hayward thermodynamics for non stationary geometries is employed in the description of the compact objects. The results of the proposed dynamical model are compared with the usual description based on stationary geometries, focusing on primordial black holes (PBHs). It is found how the timescale of evaporation depends on the choice of a cosmological observer. It may differ substantially from the treatment based on stationary models for black holes. In particular, the standard assertion that there is a fixed initial mass just below $10^{15} \, \text{g} \sim 10^{-18} M_\odot$ for the PBHs which are ending their evaporation process today is imprecise, even when possible quantum corrections at the late stages are not considered. Deviations from this prediction appear when the evaporation is measured with respect to the cosmological time.

In this article, we look at the current bounds on the coupling strength of axion-like particles (ALPs) with two photons in the context of the Randall-Sundrum (RS) model. We relate the coupling strength to the compactification radius that governs the size of the extra dimension in the RS warped geometry model and show how the current bounds on the ALP can be used to derive appropriate constraints on the size of the extra fifth dimension in the RS model. We show that the resulting constraints fail to resolve the gauge hierarchy problem for light/ultralight ALPs and require a massive ALP of at least $m_{a} \gtrsim 0.1$ [GeV] to be relevant in the context of the hierarchy problem when the gauge field is in the bulk.

Core-collapse supernovae (CCSNe) are potential multimessenger events detectable by current and future gravitational wave (GW) detectors. The GW signals emitted during these events are expected to provide insights into the explosion mechanism and the internal structures of neutron stars. In recent years, several studies have empirically derived the relationship between the frequencies of the GW signals originating from the oscillations of protoneutron stars (PNSs) and the physical parameters of these stars. This study applies the Hilbert-Huang transform (HHT) [Proc. R. Soc. A 454, 903 (1998)] to extract the frequencies of these modes to infer the physical properties of the PNSs. The results exhibit comparable accuracy to a short-time Fourier transform-based estimation, highlighting the potential of this approach as a complementary method for extracting physical information from GW signals of CCSNe.

We investigate the role of hybrid and nucleonic equations of state (EOSs) within neutron star (NS) interiors using Bayesian inference to evaluate their alignment with recent observational data from NICER and LIGO-Virgo (LV) collaborations. We find that smooth hybrid EOSs are slightly favoured in explaining NS mass-radius relations, particularly for pulsars such as PSR J0030+0451 and PSR J0740+6620. However, this preference is not definitive, as gravitational wave (GW) data does not significantly differentiate between our hybrid and nucleonic models. Our analysis also reveals tensions between older NICER data and recent measurements for PSR J0437-4715, highlighting the need for more flexible EOS models. Through two sampling approaches - one fixing the hadronic EOS set and the other without fixing the same, we demonstrate that the hybrid EOS model can incorporate stiffer EOSs, resulting in a better agreement with NICER data but leading to higher tidal deformability, which is less consistent with GW observations. In some recent publications a parameter $d_c$, related to the trace anomaly and its derivative, is used to indicate the presence of deconfined quark matter. We find that our hadronic model, which does not include phase transition to deconfined matter, under the influence of imposed constraints, is able to predict values below 0.2 for $d_c$ at around five times saturation density. The hybrid model goes below this threshold at lower densities under the same conditions.

We demonstrate that a broad class of modular inflation models predicts the emergence of new physics within an energy range of approximately $10^{15}\, \mathrm{GeV}$ to $10^{17} \, \mathrm{GeV}$. This prediction arises by comparing the moduli-dependent species scale with observational constraints on inflation. Specifically, we illustrate this within the context of $SL(2, \mathbb{Z})$-modular inflation models by re-expressing inflationary observables in terms of the species scale. We further discuss the implications of this approach for generic Calabi-Yau threefolds, showing that this reformulation allows us to directly constrain the fundamental parameters related to the geometry of extra dimensions, specifically the second Chern numbers.

We present a method for determining the physical parameters of a Kerr-Newman black hole through shadow observation. In a system comprising a Kerr-Newman black hole, an observer, and a light source, the relevant parameters are mass $M$, specific angular momentum $a$, electric charge $Q$, inclination angle $i$, and distance $r_o$. We consider the cases where the observer is at either a finite distance or spatial infinity. Using our method, the dimensionless parameters $(a/M, Q/M, i)$ can be determined by observing the shadow contour of the Kerr-Newman black hole from spatial infinity. We analytically prove that the shadow contour of the Kerr-Newman black hole observed from spatial infinity is unique, where uniqueness is defined as the absence of two congruent shadow contours for distinct sets of dimensionless parameter values. This method is versatile and can be applied to a range of black hole solutions with charge. Additionally, we show analytically that the shadow contour of a Kerr-Newman black hole observed from a finite distance $r_o$ is not unique, meaning that the parameters of a Kerr-Newman black hole at finite distance cannot be determined from shadow observations. This result reveals a new challenge and provides a clear direction for further research on black hole shadows.

The neutrino sector of the seesaw-modified Standard Model is investigated under the anarchy principle. The anarchy principle leading to the seesaw ensemble is studied analytically with tools of random matrix theory. The probability density function is obtained.

The Blandford-Znajek mechanism is an electromagnetic manifestation of the Penrose process that currently constitutes the best theoretical candidate to explain the launching of relativistic jets by black holes. In this talk we offer a modern review about the Blandford-Znajek mechanism and the analytic construction of black hole magnetospheres. Higher order perturbative corrections are crucial in order to produce results that are complementary to numerical simulations when the black hole is in the high-spin regime, and can potentially predict new features about the non-perturbative structure of the Blandford-Znajek theory. Moreover, we show by means of an explicit example that these perturbative corrections depend in a non-degenerate manner on the underlying theory of gravity considered, enabling one to use the BZ power emitted as a strong-gravity signature to test General Relativity against alternative theories of gravity on future horizon-scale observations.

We investigate the shadow properties of a rotating black hole with a weakly coupled global monopole charge, using a modified Newman-Janis algorithm. This study explores how this charge and rotational effects shape the black hole's shadow, causal structure, and ergoregions, with implications for distinguishing it from Kerr-like solutions. Analysis of null geodesics reveals observable features that may constrain the global monopole charge and weak coupling parameters within nonminimal gravity frameworks. Observational data from M87* and Sgr A* constrain the global monopole charge and coupling constant to $0 \leq \gamma \lesssim 0.036$ and $-0.2 \lesssim \alpha \leq 0$, respectively.

To enhance the ionization yield of liquid argon time projection chambers (LArTPC) used in dark matter and neutrino experiments it was proposed the use of dopants in LAr, with ionization energies below the scintillation threshold of Ar. While dual-phase LArTPCs have excellent sensitivity to single ionization electrons, their compatibility with photosensitive dopants is hindered by gas-phase electroluminescence photons ionizing the dopants, leading to a positive feedback loop. This can be addressed by optically decoupling the gaseous and liquid phases with a barrier that transmits electrons. A possible solution relies on the use of a pair of structures based on Gaseous Electron Multipliers(GEMs) with misaligned holes. Rather than amplifying electron signals in gas pockets within their holes, their holes will be filled with LAr and a low biasing voltage, so that incident drifting electrons are drawn into the holes but not amplified. Instead, amplification will occur in the gas phase above the structures. Its core element is a GEM-like structure machined from polyethylene naphthalate(PEN). Since PEN scintillates in the visible spectrum, the risk of increased radioactivity due to the larger mass compared to traditional wire grids is negated by the potential to veto its own radioactivity. As such, these structures may also be a useful alternative to wire grids. In this work, we report the newest developments on the production of GEM-like structures using laser-based techniques, namely the manufacture of the first batch PEN and PMMA-based GEM-like structures. This process allows low-cost, reproducible fabrication of a high volume of such structures. In addition to being a low radioactive technique, we expect that it will allow the scaling up of the production of these structures at a reduced cost. First tests indicate good electrical stability, while the performance assessment is still ongoing.

We analyse the evolution of the reduced density matrix of inflationary perturbations, coupled to a heavy entropic field via the leading-order term within the Effective Field Theory of Inflation, for two nearly de Sitter backgrounds. We perform a full quantum treatment of the open system and derive a Fokker-Planck equation to describe decoherence and the entanglement structure of the adiabatic perturbations. We find that exotic phenomena, such as recoherence and transient negative growth of entanglement entropy, appearing for the attractor solution, are absent for the non-attractor background. We comment on the relationship of these to the non-Markovian nature of the system. Finally, we generalise to the case where a few e-folds of ultra-slow roll evolution are sandwiched between phases of slow-roll inflation to find its (memory) effects on the curvature perturbation.

Recently, the temporal evolution of the angles characterizing the spatial configuration of the jet in the supermassive black hole M87$^\ast$ was measured exhibiting a precessional pattern around the hole's spin axis. It would be due to the dragging induced by the fact that the hole's external spacetime is described by the Kerr metric. Here, it is shown that the Lense-Thirring orbital precessions of a test particle moving about a rotating massive object, calculated perturbatively to the first post-Newtonian order, are able to fully reproduce all the measured features of the jet axis of M87$^\ast$. In particular, by assuming that the latter is aligned with the angular momentum of the accretion disk, modelled as an effective particle moving along a circular orbit, the condition that the predicted Lense-Thirring precessional frequency of the disk agrees with the measured value of $0.56\pm 0.02$ radians per year of the jet's one is satisfied for a range of physically meaningful values of the hole's spin parameter, close to unity, and of the effective disk radius, of the order of just over a dozen gravitational radii. Relying upon such assumptions and results, it is possible to predict that the angle between the hole's spin axis and the jet's one stays constant over the years amounting to $1.16^\circ$, in agreement with its measured value of $1.25^\circ\pm 0.18^\circ$. Furthermore, also the temporal pattern and the amplitudes of the time series of the jet's angles are reproduced by the aforementioned Lense-Thirring precessional model.

Axions and axion-like particles can be probed through gravitational waves indirectly, often referred to as "audible axions". The usual concept of audible axion relies on the coupling between the axions and the gauge fields. Here we consider an axion-like mechanism with coupling to the Nieh-Yan term. This interaction leads to the direct and efficient production of gravitational waves during the radiation-dominated era, originating from the tachyonic instability of the gravitational perturbations with the Nieh-Yan term. We calculate the energy spectral density of the chiral gravitational wave background and the comoving energy density of axion-like fields. Based on the numerical results, we explore the parameter space of axion masses and decay constants for detectable gravitational wave signals, either in pulsar timing arrays or space-based gravitational wave detections.

In this contribution, we extend the discussion about the calculation of the bulk viscosity of quark matter in the normal phase due to electroweak processes and its effect on the damping of baryon density oscillations that might occur in the coalescence of two compact stars. Employing the EoSs from the MIT bag model and perturbative quantum chromodynamics (pQCD) up to $\mathcal{O}(\alpha_s)$, we analyze our results varying densities in the range of temperatures from 0 to 10 MeV for frequencies around 1 kHz. Our estimates show that bulk viscous effects might play a relevant role during the postmerger stage if the system reaches a deconfined quark matter phase.

We investigate the formation of bound states of non-relativistic dark matter particles subject to long-range interactions through radiative capture. The initial scattering and final bound states are described by Coulomb potentials with different strengths, as relevant for non-abelian gauge interactions or theories featuring charged scalars. For bound states with generic quantum numbers $n$ and $\ell$, we provide closed-form expressions for the bound-state formation (BSF) cross sections of monopole, dipole and quadrupole transitions, and of arbitrary multipole order when $\ell=n-1$. This allows us to investigate in detail a strong enhancement of BSF that occurs for initial states in a repulsive potential. For $\ell=n-1\gg 1$, we show that the BSF cross section for each single bound state violates the perturbative unitarity bound in the vicinity of a certain critical initial velocity, and provide an interpretation in terms of a smooth matching of classical trajectories. When summing the BSF cross section over all possible bound states in the final state, this leads to a unitarity violation below a certain velocity, but within the validity range of the weakly coupled non-relativistic description. We identify an effectively strong interaction as the origin of this unitarity violation, which is caused by an "anomalously" large overlap of scattering and bound-state wave functions in Coulomb potentials of different strength.

Pulsar timing array experiments have recently found evidence for a stochastic gravitational wave (GW) background, which induces correlations among pulsar timing residuals described by the Hellings and Downs (HD) curve. Standard calculations of the HD correlation and its variance assume an isotropic background. However, for a background of astrophysical origin, we expect a higher GW spectral density in directions with higher galaxy number densities. In a companion paper, we have developed a theoretical formalism to account for the anisotropies arising from large-scale galaxy clustering, leading to a new contribution to the variance of the HD correlation. In this subsequent work, we provide numerical results for this novel effect. We consider a GW background resulting from mergers of supermassive black hole binaries, and relate the merger number density to the overdensity of galaxies. We find that anisotropies due to large-scale galaxy clustering lead to a standard deviation of the HD correlation at most at percent level, remaining well below the standard contributions to the HD variance. Hence, this kind of anisotropies in the GW source distribution does not represent a substantial contamination to the correlations of timing residuals in present and future PTA surveys. Suitable statistical methods to extract the galaxy clustering signal from PTA data will be investigated in the future.

By considering realistic equations of state (EoSs) to describe the ordinary matter of the stellar crust, in this study, we explore the effect of a dark energy core, made of Chaplygin Dark Fluid (CDF), on neutron stars (NSs). To accomplish this purpose, we solve the stellar structure equations and investigate the impact of the CDF parameters on the several macroscopic properties of NSs such as mass-radius ($M-R$) relation, and tidal deformabilities of a single star and of a binary system, the latter being of great importance when analyzing gravitational-wave signals coming from the merger of such compact objects. We also present an analysis of the radial oscillation modes for the rapid phase transition, with the aim of distinguishing regions consisting of dynamically stable stars from those of unstable ones. Specifically, our outcomes reveal that an increase in the energy density jump (controlled by a parameter $\alpha$) leads to an increase in the radial stability of the NS with a CDF core. Furthermore, our theoretical results are consistent with the observational $M-R$ measurements of millisecond pulsars from NICER data and tidal deformability constraints from the GW170817 event.

Radioactive nuclei with lifetimes on the order of millions of years can reveal the formation history of the Sun and active nucleosynthesis occurring at the time and place of its birth. Among such nuclei whose decay signatures are found in the oldest meteorites, $^{205}$Pb is a powerful example, as it is produced exclusively by slow neutron captures (the s process), with most being synthesized in asymptotic giant branch (AGB) stars. However, making accurate abundance predictions for $^{205}$Pb has so far been impossible because the weak decay rates of $^{205}$Pb and $^{205}$Tl are very uncertain at stellar temperatures. To constrain these decay rates, we measured for the first time the bound-state $\beta^-$ decay of fully ionized $^{205}$Tl$^{81+}$, an exotic decay mode that only occurs in highly charged ions. The measured half-life is 4.7 times longer than the previous theoretical estimate and our 10% experimental uncertainty has eliminated the main nuclear-physics limitation. With new, experimentally backed decay rates, we used AGB stellar models to calculate $^{205}$Pb yields. Propagating those yields with basic galactic chemical evolution (GCE) and comparing with the $^{205}$Pb/$^{204}$Pb ratio from meteorites, we determined the isolation time of solar material inside its parent molecular cloud. We find positive isolation times that are consistent with the other s-process short-lived radioactive nuclei found in the early Solar System. Our results reaffirm the site of the Sun's birth as a long-lived, giant molecular cloud and support the use of the $^{205}$Pb--$^{205}$Tl decay system as a chronometer in the early Solar System.