Papers for Monday, Nov 27 2023

The properties of rotating neutron stars are investigated using eight equations of state (EOSs). We also study the relations between various observables corresponding to different angular velocities for all those EOSs. All of these EoSs lead to non-rotating compact stars with maximum masses between 1.8 to 2.25 $M_{\odot}$. We calculate the moment of inertia and studied its variation with mass and the relation between central energy density and angular momentum. We compare our results with the observational findings from the most massive pulsar PSR J0740+6620 and the heaviest secondary component in the black hole - neutron star merger GW190814. It is noted that the secondary compact object of GW190814 having mass $\sim 2.6 M_{\odot}$ might be explained as a rapidly rotating neutron star (NS) with frequency larger than 1000 Hz.

To keep current global warming below 1.5{\deg}C compared with the pre-industrial era, measures must be taken as quickly as possible in all spheres of society. Astronomy must also make its contribution. In this proceeding, and during the workshop to which it refers, different levers of actions are discussed through various examples: individual efforts, laboratory-level actions, impact evaluation and mitigation in major projects, institutional level, and involvement through collectives.

Quasar absorption spectra measurements suggest that reionization proceeded rapidly, ended late at $z \sim 5.5$, and was followed by a flat evolution of the ionizing background. Simulations that can reproduce this behavior often rely on a fine-tuned galaxy ionizing emissivity, which peaks at $z \sim 6 - 7$ and drops by a factor of $1.5-2.5$ by $z \sim 5$. This is puzzling since the abundance of galaxies has been observed to grow monotonically during this period. Explanations for this include effects such as dust obscuration of ionizing photon escape and feedback due to photo-heating of the IGM. We explore the possibility that this drop in emissivity is instead an artifact of one or more modeling deficiencies in reionization simulations. These include possibly incorrect assumptions about the ionizing spectrum and/or inaccurate modeling of the clumpiness of the IGM. Our results suggest that the need for a drop could be alleviated if simulations are underestimating the IGM opacity from massive, star-forming halos. Other potential modeling issues either have a small effect or require a steeper drop when remedied. We construct an illustrative model in which the emissivity is nearly flat at the end of reionization, evolving only $\sim 0.05$ dex at $5 < z < 7$. More realistic scenarios, however, require a $\sim 0.1-0.3$ dex drop. We also study the evolution of the Ly$\alpha$ effective optical depth distribution in these scenarios and compare them to recent measurements. We find models that feature a hard ionizing spectrum and/or are driven by faint, low-bias sources can most easily reproduce the mean transmission and optical depth distribution of the forest simultaneously. Lastly, we show that the reduced speed of light approximation and low spatial resolution in the forest can lead to erroneous conclusions about the end of reionization.

In this work, we analyse the density profiles of subhaloes with masses $M_\mathrm{sh} \geq 1.4 \times 10^8$ M$_\odot$ in the TNG50 simulation, with the aim of including baryonic effects. We evaluate the performance of frequently used models, such as the standard NFW, the Einasto, and a smoothly truncated version of the NFW profile. We find that these models do not perform well for the majority of subhaloes, with the NFW profile giving the worst fit in most cases. This is primarily due to mismatches in the inner and outer logarithmic slopes, which are significantly steeper for a large number of subhaloes in the presence of baryons. To address this issue, we propose new three-parameter models and show that they significantly improve the goodness of fit independently of the subhalo's specific properties. Our best-performing model is a modified version of the NFW profile with an inner log-slope of -2 and a variable truncation that is sharper and steeper than the slope transition in the standard NFW profile. Additionally, we investigate how both the parameter values of the best density profile model and the average density profiles vary with subhalo mass, $V_\mathrm{max}$, distance from the host halo centre, baryon content and infall time, and we also present explicit scaling relations for the mean parameters of the individual profiles. The newly proposed fit and the scaling relations are useful to predict the properties of realistic subhaloes in the mass range $10^8$ M$_\odot$ $\leq M_\mathrm{sh} \leq$ $10^{13}$ M$_\odot$ that can be influenced by the presence of baryons.

The gravitational lensing signal produced by a galaxy or a galaxy cluster is determined by its total matter distribution, providing us with a way to directly constrain their dark matter content. State-of-the-art numerical simulations successfully reproduce many observed properties of galaxies and can be used as a source of mock observations and predictions. Many gravitational lensing studies aim at constraining the nature of dark matter, discriminating between cold dark matter and alternative models. However, many past results are based on the comparison to simulations that did not include baryonic physics. Here we show that the presence of baryons can significantly alter the predictions: we look at the structural properties (profiles and shapes) of elliptical galaxies and at the inner density slope of subhaloes. Our results demonstrate that future simulations must model the interplay between baryons and alternative dark matter, to generate realistic predictions that could significantly modify the current constraints.

The shell of the classical nova V5668 Sgr was resolved by ALMA at the frequency of 230 GHz 927 days after eruption, showing that most of the continuum bremsstrahlung emission originates in clumps with diameter smaller than $10^{15}$ cm. Using VLA radio observations, obtained between days 2 and 1744 after eruption, at frequencies between 1 and 35 GHz, we modeled the nova spectra, assuming first that the shell is formed by a fixed number of identical clumps, and afterwards with the clumps having a power law distribution of sizes, and were able to obtain the clump's physical parameters (radius, density and temperature). We found that the density of the clumps decreases linearly with the increase of the shell's volume, which is compatible with the existence of a second media, hotter and thinner, in pressure equilibrium with the clumps. We show that this thinner media could be responsible for the emission of the hard X-rays observed at the early times of the nova eruption, and that the clump's temperature evolution follows that of the super-soft X-ray luminosity. We propose that the clumps were formed in the radiative shock produced by the collision of the fast wind of the white dwarf after eruption, with the slower velocity of the thermonuclear ejecta. From the total mass of the clumps, the observed expansion velocity and thermonuclear explosion models, we obtained an approximate value of 1.25 M$_{\odot}$ for the mass of the white dwarf, a central temperature of $10^7$ K and an accretion rate from the secondary star of $10^{-9}-10^{-8}$ M$_{\odot}$ y$^{-1}$.

For magnetic knots and links in plasmas we introduce an internal twist and study their dynamical behavior in numerical simulations. We use a set of helical and non-helical configurations and add an internal twist that cancels the helicity of the helical configurations or makes a non-helical system helical. These fields are then left to relax in a magnetohydrodynamic environment. In line with previous works we confirm the importance of magnetic helicity in field relaxation. However, an internal twist, as could be observed in coronal magnetic loops, does not just add or subtract helicity, but also introduces an alignment of the magnetic field with the electric current, which is the source term for helicity. This source term is strong enough to lead to a significant change of magnetic helicity, which for some cases leads to a loss of the stabilizing properties expressed in the realizability condition. Even a relatively weak internal twist in these magnetic fields leads to a strong enough source term for magnetic helicity that for various cases even in a low diffusion environment we observe an inversion in sign of magnetic helicity within time scales much shorter than the diffusion time. We conclude that in solar and stellar fields an internal twist does not automatically result in a structurally stable configuration and that the alignment of the magnetic field with the electric current must be taken into account.

Smoothed particle magnetohydrodynamics has reached a level of maturity that enables the study of a wide range of astrophysical problems. In this review, the numerical details of the modern SPMHD method are described. The three fundamental components of SPMHD are methods to evolve the magnetic field in time, calculate accelerations from the magnetic field, and maintain the divergence-free constraint on the magnetic field (no monopoles). The connection between these three requirements in SPMHD will be highlighted throughout. The focus of this review is on the methods that work well in practice, with discussion on why they work well and other approaches do not. Numerical instabilities will be discussed, as well as strategies to overcome them. The inclusion of non-ideal MHD effects will be presented. A prospective outlook on possible avenues for further improvements will be discussed.

Recently, GRB 221009A, known as the brightest of all time (BOAT), has been observed across an astounding range of $\sim 18$ orders of magnitude in energy, spanning from radio to VHE bands. Notably, the Large High Altitude Air Shower Observatory (LHAASO) recorded over $60000$ photons with energies exceeding $0.2\rm~TeV$, including the first-ever detection of photons above $10\rm~TeV$. However, explaining the observed energy flux evolution in the VHE band alongside late-time multi-wavelength data poses a significant challenge. Our approach involves a two-component structured jet model, consisting of a narrow core dominated by magnetic energy and a wide jet component dominated by matter. We show that the combination of the forward shock electron synchrotron self-Compton emission from both jets and reverse shock proton synchrotron emission from the wide jet could account for both the energy flux and spectral evolution in the VHE band, and the early TeV lightcurve may be influenced by prompt photons which could explain the initial steep rising phase. We noticed the arrival time of the highest energy photons detected by LHAASO-KM2A coincident with the peak of the reverse shock proton synchrotron emission, especially a minor flare occurring about $\sim500-800$ seconds after the trigger, coinciding with the observed spectral hardening and arrival time of the $\sim 13\rm~TeV$ photons detected by LHAASO. These findings imply that the GRB reverse shock may serve as a potential accelerator of ultra-high-energy cosmic rays, a hypothesis that could be tested through future multimessenger observations.

We obtained a 108-ks Chandra X-ray Observatory (CXO) observation of PSR J1849-0001 and its PWN coincident with the TeV source HESS J1849-000. By analyzing the new and old (archival) CXO data we resolved the pulsar from the PWN, explored the PWN morphology on arcsecond and arcminute scales, and measured the spectra of different regions of the PWN. Both the pulsar and the compact inner PWN spectra are hard with power-law photon indices of $1.20\pm0.07$ and $1.36\pm0.14$, respectively. The jet-dominated PWN with a relatively low luminosity, the lack of gamma-ray pulsations, the relatively hard and non-thermal spectrum of the pulsar, and its sine-like pulse profile, indicate a relatively small angle between the pulsar's spin and magnetic dipole axis. In this respect, it shares similar properties with a few other so-called MeV pulsars. Although the joint X-ray and TeV SED can be roughly described by a single-zone model, the obtained magnetic field value is unrealistically low. A more realistic scenario is the presence of a relic PWN, no longer emitting synchrotron X-rays but still radiating in TeV via inverse-Compton upscattering. We also serendipitously detected surprisingly bright X-ray emission from a very wide binary whose components should not be interacting.

We describe {\it JWST}/NIRSpec prism measurements of Ly$\alpha$ emission in $z\gtrsim 5$ galaxies. We identify Ly$\alpha$ detections in 10 out of 69 galaxies with robust rest-optical emission line redshift measurements at $5\leq z<7$ in the CEERS and DDT-2750 observations of the EGS field. Galaxies at $z\simeq 6$ with faint continuum (F150W $=$ 27--29 mag) are found with extremely large Ly$\alpha$ equivalent widths (ranging up to 286 A). Likely Ly$\alpha$ detections are also seen in two new $z>7$ galaxies ($z=$ 7.49 and 7.17) from the second epoch of CEERS observations, both of which show large Ly$\alpha$ equivalent widths that likely indicate significant transmission through the IGM. We measure high Ly$\alpha$ escape fractions in the 12 Ly$\alpha$ emitters in our sample (median 0.28), two of which show $f_{\rm esc}^{ {\rm Ly}\alpha}$ near unity ($>0.80$). We find that $50_{-11}^{+11}$% of $z\simeq 6$ galaxies with [OIII]+H$\beta$ EW $>$ 1000 A have $f_{\rm esc}^{ {\rm Ly}\alpha}$ $>0.2$, consistent with the fractions found in lower-redshift samples with matched [OIII]+H$\beta$ EWs. While uncertainties are still significant, we find that only $10_{-5}^{+9}$% of $z>7$ galaxies with similarly strong rest optical emission lines show such large $f_{\rm esc}^{ {\rm Ly}\alpha}$, as may be expected if IGM attenuation of Ly$\alpha$ increases towards higher redshifts. We identify photometric galaxy overdensities near the $z\gtrsim 7$ Ly$\alpha$ emitters, potentially providing the ionizing flux necessary to create large ionized sightlines that facilitate Ly$\alpha$ transmission. Finally, we investigate the absence of Ly$\alpha$ emission in a comparable (and spectroscopically confirmed) galaxy overdensity at $z=7.88$ in the Abell 2744 field, discussing new prism spectra of the field obtained with the UNCOVER program.

We discuss the chosen concepts, detailed design, implementation and calibration of the 85-electrode adaptive optics (AO) system of the Swedish 1-meter Solar Telescope (SST). The AO system includes a 52 mm diameter monomorph deformable mirror with 34 mm pupil diameter, and an Intel PC workstation that performs the heavy image processing associated with cross correlations and real-time control at 2 kHz update rate with very low latency. The AO system is unusual by using a combination of a monomorph mirror with a Shack-Hartmann (SH) wavefront sensor (WFS), and uses a second high-resolution SH microlens array to aid the DM characterization, calibration, and modal control. The computer and software continue the successful implementation since 1995 of earlier generations of correlation tracker and AO systems at SST and its predecessor SVST by relying entirely on work station technology and an extremely efficient algorithm for implementing cross correlations with the large field-of-view of the WFS. We describe critical aspects of the design, calibrations, software and functioning of the AO system. The exceptionally high performance is testified through the highest Strehl ratio (as inferred from the measured granulation contrast) of existing meter-class solar telescopes, and in particular the unparalleled image quality achieved at wavelengths below 400 nm. We expect that some aspects of this AO system may be of interest also outside the solar community.

We present the hydrodynamic model of Type IIP SN 2020jfo with the unusually short (nearly 60 days) light curve plateau. The model suggests the explosion of about 8 Msun red supergiant that ejected nearly 6 Msun with the energy of about 0.8x10^51 erg. The presupernova wind density turns out highest among known SNe IIP. Yet the presupernova was not embedded into a very dense confined circumstellar shell that is a feature of some Type IIP supernovae, so the circumstellar interaction in close environment does not contribute noticeably to the initial (about 10 days) bolometric luminosity. Despite uncommon appearance SN 2020jfo turns out similar to SN 1970G in the V-band light curve, photospheric velocities, and, possibly, luminosity as well.

The origin of thick disks and their evolutionary connection with thin disks are still a matter of debate. We provide new insights into this topic by connecting the stellar populations of thick disks at redshift $z=0$ with their past formation and growth, in 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We projected each galaxy edge on, and decomposed it morphologically into two disk components, in order to define geometrically the thin and the thick disks as usually done in observations. We produced age, metallicity and [Mg/Fe] edge-on maps. We quantified the impact of satellite mergers by mapping the distribution of ex-situ stars. Thick disks are on average $\sim 3$~Gyr older, $\sim 0.25$~dex more metal poor and $\sim 0.06$~dex more [Mg/Fe]-enhanced than thin disks. Their average ages range from $\sim 6$ to $\sim 9$~Gyr, metallicities from $\sim -0.15$ to $\sim 0.1$~dex, and [Mg/Fe] from $\sim 0.12$ to $\sim 0.16$~dex. These properties are the result of an early initial in-situ formation, followed by a later growth driven by the combination of direct accretion of stars, some in-situ star formation fueled by mergers, and dynamical heating of stars. The balance between these processes varies from galaxy to galaxy. Mergers play a key role in the mass assembly of thick disks, contributing an average accreted mass fraction of $\sim 22$\% in the analyzed thick-disk dominated regions. In two galaxies, about half of the geometric thick-disk mass was directly accreted. While primordial thick disks form at high redshift in all galaxies, young metal-rich thin disks, with much lower [Mg/Fe] abundances, start to form later but at different times (higher or lower redshift) depending on the galaxy. We conclude that thick disks result from the interplay of external processes with the internal evolution of the galaxy.

Interactions of cosmic ray protons and nuclei in their sources and in the interstellar medium produce "hadronic" gamma-ray emission. Gamma-rays can also be of "leptonic" origin, i.e. originating from high-energy electrons accelerated together with protons. It is difficult to distinguish between hadronic and leptonic emission mechanisms based on gamma-ray data alone. This can be done via detection of neutrinos, because only hadronic processes lead to neutrino production. We use publicly available ten-year IceCube neutrino telescope dataset to demonstrate the hadronic nature of high-energy emission from the direction of Cygnus region of the Milky Way. We find a 3-sigma excess of neutrino events from an extended Cygnus Cocoon, with the flux comparable to the flux of gamma-rays in the multi-TeV energy range seen by HAWC and LHAASO telescopes.

Binary black holes embedded within the accretion disks that fuel active galactic nuclei (AGN) are promising progenitors for the source of gravitational wave events detected by LIGO/VIRGO. Several recent studies have shown that when these binaries form they should be highly eccentric and retrograde. However, many uncertainties remain concerning the orbital evolution of these binaries as they either inspiral towards merger or disassociate. Previous hydrodynamical simulations exploring their orbital evolution have been predominantly two-dimensional, or have been restricted to binaries on nearly circular orbits. We present the first high-resolution, three-dimensional local shearing-box simulations of both prograde and retrograde eccentric binary black holes embedded in AGN disks. We find that retrograde binaries shrink several times faster than their prograde counterparts and exhibit significant orbital eccentricity growth, the rate of which monotonically increases with binary eccentricity. Our results suggest that retrograde binaries may experience runaway orbital eccentricity growth, which may bring them close enough together at pericenter for gravitational wave emission to drive them to coalescence. Although their eccentricity is damped, prograde binaries shrink much faster than their orbital eccentricity decays, suggesting they should remain modestly eccentric as they contract towards merger. Finally, binary precession driven by the AGN disk may dominate over precession induced by the supermassive black hole depending on the binary accretion rate and its location in the AGN disk, which can subdue the evection resonance and von Ziepel-Lidov-Kozai cycles.

In this work, we analyze a power-law inflationary potential enhanced with a step that can introduce features in the primordial power spectrum. We focus on the computation of the Spectral Distortions (SD) induced by these features obtained from the inflationary dynamics. In this scenario, we explore the potential of upcoming experimental missions like PIXIE to detect the SD of the model within a power of $n = 2/3$, a power that agrees with recent tensor-to-scalar ratio constraints. The model offers insights into models with cosmological phases and different scalar field dynamics. Introducing a step in the inflaton potential leads to distinct features in the primordial power spectrum, such as oscillations and localized enhancements/suppressions at specific scales. We analyze the impact of three primary parameters$-\beta$, $\delta$, and $\phi_{\text{step}}-$on the amplitude and characteristics of the SD. The $\phi_{\rm step}$ places the onset of the oscillations in the primordial power spectrum. The $\beta$ parameter significantly influences the magnitude of the $\mu$-SD, with its increase leading to larger SD and vice versa. Similarly, the $\delta$ parameter affects the smoothness of the step in the potential, with larger values resulting in smaller SD. Our findings indicate a distinct parameter space defined by $0.02 <\delta/{\rm M_{pl}} \lesssim 0.026$, $0.10 \lesssim \beta < 0.23$, and $ 7.53 \lesssim \phi_{\rm step}/{\rm M_{pl}} \lesssim 7.55$, which produces SD potentially detectable by PIXIE. This region also corresponds to the maximum observed values of $\mu$ and $y$ SD, which in special cases are an order of magnitude larger than the expected for $\Lambda$CDM. However, we also identify parameter ranges where $\mu$ and $y$ SD may not be detectable due to the limitations of current observational technology.

In this work, we analytically calculate the spectra of primordial perturbations at the end of Type III hilltop inflation models under the slow-roll approximation. We examine the one-loop corrections of the spectra and find that those from the inflaton self-interaction are negligible. On the contrary, the loop effects from the interaction between the inflaton field and the waterfall field can be significant when the vacuum expectation value of the waterfall field is small. The implications are discussed.

The observed scaling relations between supermassive black hole masses and their host galaxy properties indicate that supermassive black holes influence the evolution of galaxies. However, the scaling relations may be affected by selection biases. We propose to measure black hole masses in a mass-limited galaxy sample including all non-detections to inprove constraints on galaxy mass - black hole mass scaling relations and test for selection bias. We use high spatial resolution spectroscopy from the Keck and Gemini telescopes, and the Jeans Anisotropic Modelling method to measure black hole masses in early type galaxies from the Virgo Cluster. We present four new black hole masses and one upper limit in our mass-selected sample of galaxies of galaxy mass (1.0-3.2) x $10^{10} M_\odot$. This brings the total measured to 11 galaxies out of a full sample of 18 galaxies, allowing us to constrain scaling relations. We calculate a lower limit for the average black hole mass in our sample of 3.7 x $10^{7} M_\odot$. This is at an average galaxy stellar mass of (1.81 +/- 0.14) x $10^{10} M_\odot $ and an average bulge mass of (1.31 +/- 0.15) x $10^{10} M_\odot$. This lower limit shows that black hole masses in early type galaxies are not strongly affected by selection biases.

The Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS), a JWST GTO program, obtained a set of unique NIRCam observations that have enabled us to significantly improve the default photometric calibration across both NIRCam modules. The observations consisted of three epochs of 4-band (F150W, F200W, F356W, and F444W) NIRCam imaging in the Spitzer IRAC Dark Field (IDF). The three epochs were six months apart and spanned the full duration of Cycle 1. As the IDF is in the JWST continuous viewing zone, we were able to design the observations such that the two modules of NIRCam, modules A and B, were flipped by 180 degrees and completely overlapped each other's footprints in alternate epochs. We were therefore able to directly compare the photometry of the same objects observed with different modules and detectors, and we found significant photometric residuals up to ~ 0.05 mag in some detectors and filters, for the default version of the calibration files that we used (jwst_1039.pmap). Moreover, there are multiplicative gradients present in the data obtained in the two long-wavelength bands. The problem is less severe in the data reduced using the latest pmap (jwst_1130.pmap as of September 2023), but it is still present, and is non-negligible. We provide a recipe to correct for this systematic effect to bring the two modules onto a more consistent calibration, to a photometric precision better than ~ 0.02 mag.

We investigated the correlation between intensities of the $^{12}$CO and $^{13}$CO ($J=1$-0) lines toward the Galactic giant molecular clouds (GMCs) W51A, W33, N35-N36 complex, W49A, M17SW, G12.02-00.03, W43, and M16 using the FUGIN (FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45-m telescope) CO line data. All the GMCs show intensity saturation in the $^{12}$CO line when the brightness temperature of $^{13}$CO is higher than a threshold temperature of about $\sim 5$ K. We obtained high-resolution ($\sim 20"$) distribution maps of the $X_{\rm CO}$ factor ($X_{\rm CO, iso}$) in individual GMCs using correlation diagrams of the CO isotopologues. It is shown that $X_{\rm CO, iso}$ is variable in each GMC within the range of $X_{\rm CO, iso} \sim (0.9 {\rm -} 5) \times 10^{20}$ cm$^{-2}$ (K km s$^{-1})^{-1}$. Despite the variability in the GMCs, the average value among the GMCs is found to be nearly constant at $X_{\rm CO, iso} = (2.17 \pm 0.27) \times 10^{20}$ cm$^{-2}$ (K km s$^{-1})^{-1}$, which is consistent with that from previous studies in the Milky Way.

Magnetars are highly-magnetised rotating neutron stars that are predominantly observed as high-energy sources. Six of this class of neutron star are known to also emit radio emission, and magnetars are, thus, a favoured model for the origin for at least some of the Fast Radio Bursts (FRBs). If magnetars, or neutron stars in general, are indeed responsible, sharp empirical constraints on the mechanism producing radio emission are required. Here we report on the detection of polarised quasi-periodic sub-structure in the emission of all well-studied radio-detected magnetars. A correlation previously seen, relating sub-structure in pulsed emission of radio emitting neutron stars to their rotational period, is extended, and shown to now span more than six of orders of magnitude in pulse period. This behaviour is not only seen in magnetars but in members of all classes of radio-emitting rotating neutron stars, regardless of their evolutionary history, their power source or their inferred magnetic field strength. If magnetars are responsible for FRBs, it supports the idea of being able to infer underlying periods from sub-burst timescales in FRBs.

Dark matter dominates the matter budget of the universe but its nature is unknown. Deviations from the standard model, where dark matter clusters with the same gravitational strength as baryons, and has the same pressureless equation of state as baryons, can be tested by cosmic growth measurements. We take a model independent approach, allowing deviations in bins of redshift, and compute the constraints enabled by ongoing cosmic structure surveys through redshift space distortions and peculiar velocities. These can produce constraints at the $3-14\%$ level in four independent redshift bins over $z=[0,4]$.

An association of GW190425 and FRB 20190425A had been claimed recently. Given the $\sim 2.5$ hour delay of the occurrence of FRB 20190425A, a uniformly rotating supramassive magnetar remnant is favored. The required maximum gravitational mass of the nonrotating neutron star (NS) is $M_{\rm TOV}\approx 2.77M_\odot$, which is strongly in tension with the low $M_{\rm TOV}\approx 2.25M_\odot$ obtained in current equation of state (EOS) constraints incorporating perturbative quantum chromodynamics (pQCD) information. However, the current mass-radius and mass-tidal deformability measurements of NSs alone do not convincingly exclude the high $M_{\rm TOV}$ possibility. By performing EOS constraints with mock measurements, we find that with a $2\%$ determination for the radius of PSR J0740+6620-like NS it is possible to distinguish between the low and high $M_{\rm TOV}$ scenarios. We further explore the prospect to resolve the issue of the appropriate density to impose the pQCD constraints with future massive NS observations or determinations of $M_{\rm TOV}$ and/or $R_{\rm TOV}$. It turns out that measuring the radius of a PSR J0740+6620-like NS is insufficient to probe the EOSs around 5 nuclear saturation density, where the information from pQCD becomes relevant. The additional precise $M_{\rm TOV}$ measurements, anyhow, could help. Indeed, supposing the central engine of GRB 170817A is a black hole formed via the collapse of a supramassive NS, the resulting $M_{\rm TOV}\approx 2.2M_\odot$ considerably softens the EOS at the center of the most massive NS, which is in favor of imposing the pQCD constraint at density beyond the one achievable in the NSs.

RR Lyrae stars play a central role in tracing phase-space structures within the Milky Way because they are easy to identify, are relatively luminous, and are found in large numbers in the Galactic bulge, disk, and halo. In this work, we present a new set of spectroscopic metallicity calibrations that use the equivalent widths of the Ca II K and Balmer H-gamma and H-delta lines to calculate metallicity values from low-resolution spectra. This builds on an earlier calibration from Layden by extending the range of equivalent widths which map between Ca II K and the Balmer lines. We have developed the software rrlfe to apply this calibration to spectra in a consistent, reproducible, and extensible manner. This software is open-source and available to the community. The calibration can be updated with additional datasets in the future.

We study zero and finite temperature static Bose-Einstein condensate (BEC) stars in the combined Rastall-Rainbow (RR) theory of gravity by considering different BEC equation of states (EoSs). We obtain the global properties of BEC stars by solving the modified Tolman-Oppenheimer-Volkoff equations with values of Rastall parameter $\kappa$ and Rainbow function $\Sigma$ chosen accordingly to get the results in theories of Rastall, Rainbow and RR. We observe that the parameter $\kappa$ has negligible effect on the maximum mass of the stars considered, whereas $\Sigma$ alters it significantly, and increasing the value of $\kappa$ beyond a certain limit results in unstable solutions for any value of $\Sigma$. We report that the inclusion of temperature in our analysis expands the parameter space by including more values of $\kappa$. However, temperature has negligible effect on the maximum mass of the stellar profiles in all the three theories. We find that the maximum masses and radii of the stars within RR theory can have good agreement with the observational data on pulsars for all the EoSs considered and in particular, the Colpi-Wasserman-Shapiro EoS, which was ruled out in General Relativity (GR). We also find that, in contrast to the results of GR, BEC stars consistent with observations can be realised in the RR theory with smaller bosonic self-interaction strength.

Narrow-line Seyfert 1 (NLSy1) galaxies are an enigmatic class of active galactic nuclei (AGN) that exhibit peculiar multiwavelength properties across the electromagnetic spectrum. For example, these sources have allowed us to explore the innermost regions of the central engine of AGN using X-ray observations and have also provided clues about the origin of relativistic jets considering radio and gamma-ray bands. Keeping in mind the ongoing and upcoming wide-field, multi-frequency sky surveys, we present a new catalogue of NLSy1 galaxies. This was done by carrying out a detailed decomposition of >2 million optical spectra of quasars and galaxies from the Sloan Digital Sky Survey Data Release 17 (SDSS-DR17) using the publicly available software "Bayesian AGN Decomposition Analysis for SDSS Spectra". The catalogue contains 22656 NLSy1 galaxies which is more than twice the size of the previously identified NLSy1s based on SDSS-DR12. As a corollary, we also release a new catalogue of 52273 broad-line Seyfert 1 (BLSy1) galaxies. The estimated optical spectral parameters and derived quantities confirm the previously known finding of NLSy1 galaxies being AGN powered by highly accreting, low-mass black holes. We conclude that this enlarged sample of NLSy1 and BLSy1 galaxies will enable us to explore the low-luminosity end of the AGN population by effectively utilizing the sensitive, high-quality observations delivered by ongoing/upcoming wide-field sky surveys. The catalogue has been made public at https://www.ucm.es/blazars/seyfert

We present the first direct constraints on a Degenerate Higher Order Scalar Tensor (DHOST) inflation model using the Planck 2018 Cosmic Microwave Background (CMB) results on non-Gaussianities. We identify that the bispectrum consists of a fixed contribution following from the power spectrum and a linear combination of terms depending on five free parameters defining the cubic perturbations to the DHOST model. The former peaks in the squeezed limit, while the latter is maximised in the equilateral limit. We directly confront the model predictions to the CMB bispectrum statistics via the public code CMB-BEST and marginalize over the free parameters. We explicitly show that there are viable DHOST inflationary models satisfying both power spectrum and bispectrum constraints from Planck. However, rather surprisingly, the constraints exclude certain models at the $6\sigma$-level even though they pass the conventional fudge factor tests. In this case and despite having a handful of free parameters, the model's large squeezed bispectrum cannot be cancelled out without introducing a large bispectrum in other limits which are strongly constrained by Planck's non-detection of primordial non-Gaussianity. We emphasize that first-order approximations such as fudge factors, albeit commonly used in the literature, may be misleading and provide weaker constraints. A proper analysis of the constraints from Planck requires a more robust approach, such as the one provided by the CMB-BEST code.

We assert that non-linear features of fast neutrino-flavor conversion (FFC) can be qualitatively different between core-collapse supernovae (CCSNe) and binary neutron star mergers (BNSMs). This argument arises from recent global FFC simulations in BNSM, in which fast flavor swap (FFS) emerges in very narrow spatial regions, whereas neutrinos in CCSN tend to evolve towards flavor equipartition. In this {\it Letter}, we provide the physical mechanism of FFS based on a colliding neutrino beam model. Neutrinos/antineutrinos can undergo FFS when they propagate in ambient neutrino gas that propagates in the opposite direction and also has the opposite sign of ELN-XLN, where ELN and XLN denote electron- and heavy-leptonic neutrino number, respectively. Such environments can be naturally realized in BNSMs, whereas they are unlikely in CCSNe unless the neutrino sphere is strongly deformed aspherically. Our study exhibits the diversity of non-linear dynamics of FFC.

KS 1947+300 is a Be/X-ray binary. Despite its nearly circular orbit, it displays both giant and regular less intense X-ray outbursts. According to the viscous decretion disk model, such low eccentric binaries should not show periodic outbursts. We have studied the long-term optical variability of KS 1947+300 and its relationship with X-ray activity. Our objective is to investigate the origin of this variability. KS 1947+300 exhibits changes in brightness and H$\alpha$ emission on time scales from months to years. The optical and IR variability shows small amplitude changes during the active X-ray state, and a long, smooth decrease during the quiescent state. The fact that the amplitude of variability increased with wavelength suggests that the long-term decrease of the optical emission is due to the weakening of the circumstellar disk. Structural changes in the disk may also be the origin of the periodic signals with periods $\sim200$ days detected in the ZTF-{\it g} and {\it r} band light curves. We speculate that these changes are related to the mechanism that ejects matter from the photosphere of the Be star into the disk. We have also studied the X-ray variability that manifested as a series of type I outbursts after the two giant outburst in 2000 and 2013 and found that the intensity and peak orbital phase differ from outburst to outburst. The type I outbursts in KS 1947+300 are not strictly periodic. This irregularity could result from precession of the interacting points between the neutron star and the disk, namely the disk apastron and the two nodes of the disk. The long-term changes in optical continuum and line emission and the X-ray variability patterns are attributed to an evolving and distorted decretion disk.

We report the time resolved spectroscopy result from two observations of Cen X-3, over one binary orbit with ASTROSAT and two binary orbit with NuSTAR. NuSTAR covered two intensity states where the light curve showed transition in count rate from first to second binary orbit by a factor of $\sim$ 3. A phenomenological model comprising of partially absorbed powerlaw with smoothed high energy cutoff, cyclotron absorption $\sim$ 24 keV and 6.4 keV iron emission gave good fit for ASTROSAT observation. NuSTAR spectra required two additional emission components, a broad one $\sim$ 5.7 keV and a narrow one $\sim$ 6.9 keV. A weak secondary absorption feature at $\sim$ 11.6 keV and $\sim$ 14.5 keV was seen in the residuals of the spectral fit for ASTROSAT and NuSTAR data respectively. The secondary absorption energy showed no correlation with the cutoff energy. Its strength varied within 0.1 to 0.6 keV with its width $\sim$ 1.6 keV. Its energy and optical depth showed linear positive correlation with the fundamental cyclotron line energy and depth respectively. The cyclotron line energy showed anti-correlation to flux described by a powerlaw with negative index and the secondary absorption also showed similar trend to flux. Depth of secondary absorption was $\sim$ 45 \% and centroid energy was $\sim$ 54 \% of fundamental. Depth and energy ratio of secondary to fundamental lied within 2$\sigma$ deviation from 0.5. We suggest this secondary absorption to be a redshifted dipolar cyclotron resonance feature exhibiting sub-harmonic behaviour.

In this work, we present the detection of twelve doublets with quantum numbers of N=12-11 to N=17-16 of the v11 vibrationally excited state of C6H towards TMC-1. This marks the first time that an excited vibrational state of a molecule has been detected in a cold starless core. The data are part of the QUIJOTE line survey gathered with the Yebes 40m radio telescope. The line intensities have been aptly reproduced with a rotational temperature of 6.2 +/- 0.4K and a column density of (1.2+/-0.2)e11 cm-2. We also analysed the ground state transitions of C6H, detecting fourteen lines with quantum numbers of J = 23/2-21/2 to J = 35/2 for each of the two 2Pi_3/2 and 2Pi_1/2 ladders. It is not possible to model the intensities of all the transitions of the ground state simultaneously using a single column density. We considered the two ladders as two different species and found that the rotational temperature is the same for both ladders, Trot(2Pi_3/2)=Trot(2Pi_1/2)=6.2+/-0.2, achieving a result that is comparable to that of the v11 state. The derived column densities are N(2Pi_3/2)(6.2+/-0.3)e12cm-2 and N(2Pi_1/2)=(8.0+/-0.4)e10cm-2. The fraction of C6H molecules in its 2Pi_3/2, 2Pi_1/2, and v11 states is 96.8 %, 1.3 %, and 1.9 %, respectively. Finally, we report that this vibrational mode has also been detected towards the cold cores Lupus-1A and L1495B, as well as the low-mass star forming cores L1527 and L483, with fractions of C6H molecules in this mode of 3.8%, 4.1%, 14.8%, and 6%, respectively.

Blazars, a subclass of active galactic nuclei (AGN), are known to be bright $\gamma$-ray sources, frequently exhibiting active (flaring) periods. The blazar PKS~2155-304 is a high synchrotron-peaked BL Lac object located at redshift $z=0.116$. On 2006 July 28, an extremely remarkable outburst of VHE $\gamma$-ray emission from this blazar was reported by the H.E.S.S. experiment, with an average flux more than 10 times the low-state level. The variability timescale of this extraordinary flare was as short as approximately 200~s. In order to guarantee the transparency of the emission region for TeV photons, the fast variability demands an extremely high Doppler factor $\delta_{\rm D}>50$ of the jet within the classical one-zone model, leading to the so-called ``Doppler factor crisis''. Here we demonstrate that the stochastic dissipation model, which is a multi-blob scenario for blazars, can self-consistently explain the giant TeV flares of PKS~2155-304 and the low-state emission before and after the flares, in terms of both multi-wavelength spectral and variability characteristics. The required Doppler factor in this model can be as low as 20, which is a reasonable and typical value for blazar jets. The obtained model parameters may shed some light on the physical properties of the relativistic jet.

Cloud cover plays a pivotal role in assessing observational conditions for astronomical site-testing. Except for the fraction of observing time, its fragmentation also wields a significant influence on the quality of nighttime sky clarity. In this article, we introduce the function Gamma, designed to comprehensively capture both the fraction of available observing time and its continuity. Leveraging in situ measurement data gathered at the Muztagh-ata site between 2017 and 2021, we showcase the effectiveness of our approach. The statistical result illustrates that the Muztagh-ata site affords approximately 122 nights of absolute clear and 205 very good nights annually, corresponding to Gamma greater than or equal 0.9 and Gamma greater than or equal 0.36 respectively.

Accurately measuring the Hubble parameter is vital for understanding the expansion history and properties of the universe. In this paper, we propose a new method that supplements the covariance between redshift pairs to improve the reconstruction of the Hubble parameter using the OHD dataset. Our approach utilizes a cosmological model-independent radial basis function neural network (RBFNN) to describe the Hubble parameter as a function of redshift effectively. Our experiments show that this method results in a reconstructed Hubble parameter of $H_0 = 67.1\pm9.7~\mathrm{km~s^{-1}~Mpc^{-1}}$ , which is more noise-resistant and fits better with the $\Lambda$CDM model at high redshifts. Providing the covariance between redshift pairs in subsequent observations will significantly improve the reliability and accuracy of Hubble parametric data reconstruction. Future applications of this method could help overcome the limitations of previous methods and lead to new advances in our understanding of the universe.

SN 2019neq was a very fast evolving superluminous supernova. At a redshift z=0.1059, its peak absolute magnitude was -21.5+/-0.2 mag in g band. In this work, we present data and analysis from an extensive spectrophotometric follow-up campaign using multiple observational facilities. Thanks to a nebular spectrum of SN 2019neq, we investigated some of the properties of the host galaxy at the location of SN 2019neq and found that its metallicity and specific star formation rate are in a good agreement with those usually measured for SLSNe-I hosts. We then discuss the plausibility of the magnetar and the circumstellar interaction scenarios to explain the observed light curves, and interpret a nebular spectrum of SN 2019neq using published SUMO radiative-transfer models. The results of our analysis suggest that the spindown radiation of a millisecond magnetar with a magnetic field B~6e14 G could boost the luminosity of SN 2019neq.

Coronal mass ejections (CMEs) are solar eruptions of plasma and magnetic fields that significantly impact Space Weather, causing disruptions in technological systems and potential damage to power grids when directed towards Earth. Traditional coronagraphs along the Sun-Earth line struggle to precisely track the early evolution of Earth-directed CMEs. Coronal dimmings, localized reductions in extreme-ultraviolet (EUV) and soft X-ray emissions, are key indicators of CMEs in the low corona, resulting from mass loss and expansion during the eruption. This study introduces a novel method, DIRECD (Dimming InfeRred Estimate of CME Direction), to estimate the early propagation direction of CMEs based on the expansion of coronal dimmings. The approach involves 3D simulations of CMEs using a geometric cone model, exploring parameters like width, height, source location, and deflection from the radial direction. The dominant direction of dimming evolution is then determined, and an inverse problem is solved to reconstruct an ensemble of CME cones at various heights, widths, and deflections. By comparing the CME orthogonal projections onto the solar sphere with the dimming geometry, the 3D CME direction is derived.Validated through case studies on October 1, 2011, and September 6, 2011, the DIRECD method reveals the early propagation directions of CMEs. The CME on October 1, 2011, predominantly expands towards the South-East, while the CME on September 6, 2011, inclines towards the North-West. These findings align with previous studies using multi-viewpoint coronagraphic observations. The study demonstrates the utility of coronal dimming information for early CME direction estimation, providing valuable data for space weather forecasting and mitigating potential adverse impacts on Earth before observation in coronographs' field-of-view.

We develop a spherical self-similar model for the formation of a galaxy through gas collapsing in an isolated self-gravitating dark matter halo. We improve upon the existing literature on self-similar collapse in two ways. First, we include the effects of radiative cooling and the formation of a pseudo-disk at the center of collapse, in a parametrised manner. More importantly, for the first time we solve for the evolution of gas and dark matter simultaneously and self-consistently using a novel iterative approach. As a result, our model produces shell trajectories of both gas and dark matter that qualitatively agree with the results of full hydrodynamical simulations. We discuss the impact of various ingredients such as the accretion rate, gas equation of state, disk radius and cooling rate amplitude on the evolution of the gas shells. The self-consistent evolution of gas and dark matter allows us to study the response of the dark matter trajectories to the presence of collapsing gas, an effect that has gained increasing importance recently in the context of precision estimates of small-scale statistics like the matter power spectrum. Our default configuration produces a relaxation relation in qualitative agreement with that seen in cosmological hydrodynamical simulations, and further allows us to easily study the impact of the model ingredients mentioned above. As an initial application, we vary one ingredient at a time and find that the accretion rate and gas equation of state have the largest impact on the relaxation relation, while the cooling amplitude plays only a minor role. Our model thus provides a convenient framework to rapidly explore the coupled nonlinear impact of multiple astrophysical processes on the mass and velocity profiles of dark matter in galactic halos. (Abridged)

We made statistical analysis of the Fermi GBM and Swift BAT observational material, accumulated over 15 years. We studied how GRB parameters (T$_{90}$ duration, fluence, peak flux) that were observed by only one satellite differ from those observed by both. In the latter case, it was possible to directly compare the values of the parameters that both satellites measured. The GRBs measured by both satellites were identified using the knn() k-nearest neighbour algorithm in the FNN library of the R statistical package. In the parameter space we determined the direction in which the jointly detected GRBs differ most from those detected by only one of the instruments using the lda() in MASS library of R. To get the strength of the relationship between the parameters obtained from the GBM and BAT, a canonical correlation was performed using the cc() procedure in the CCA library in R. The GBM and BAT T$_{90}$ distributions were fitted with a linear combination of lognormal functions. The optimal number of such functions required for fit is two for GBM and three for BAT. Contrary to the widely accepted view, we found that the number of lognormal functions required for fitting the observed distribution of GRB durations does not allow us to deduce the number of central engine types responsible for GRBs.

The transport of angular momentum (AM) and chemical elements within evolving stars remains poorly understood. Recent observations showed that the radiative cores of low mass main sequence stars and red giants rotate orders of magnitude slower than classical stellar evolution models predictions and that their surface light elements abundances are too small. Magnetohydrodynamic (MHD) turbulence can enhance the transport in radiative stellar interiors but its efficiency is still largely uncertain. Here we explore the transport of AM and chemical elements due to azimuthal magnetorotational instability (AMRI) using 3D MHD direct numerical simulations in a spherical shell. First, we provide evidence of AMRI in the parameter regime expected from local and global linear stability studies and then we analyze its nonlinear evolution. For unstratified flow, we observe AMRI-driven dynamo action at values of the magnetic Prandtl number Pm in the range $0.6-1$, the smallest ever reported in a global setup. When considering stable stratification (under the Boussinesq approximation) at $\text{Pm}=1$, the turbulence is instead transitional and becomes less homogeneous and isotropic upon increasing buoyancy effects. We find that the transport of AM occurs radially outwards and is dominated by the Maxwell stresses when stratification is large enough. The associated turbulent viscosity decreases when buoyancy effects strengthen and scales with the square root of the ratio of the rotation rate to the Brunt-V\"ais\"al\"a frequency, predicting that, for example, red giant cores may reach a state of uniform rotation in a few thousand years. A passive scalar allows us to study the transport of chemical elements. The chemical turbulent diffusion coefficient scales with stratification similarly to the turbulent viscosity but is lower in amplitude as suggested by recent stellar evolution models of low mass stars.

The Vera C. Rubin Observatory is preparing to execute the most ambitious astronomical survey ever attempted, the Legacy Survey of Space and Time (LSST). Currently the final phase of construction is under way in the Chilean Andes, with the Observatory's ten-year science mission scheduled to begin in 2025. Rubin's 8.4-meter telescope will nightly scan the southern hemisphere collecting imagery in the wavelength range 320-1050 nm covering the entire observable sky every 4 nights using a 3.2 gigapixel camera, the largest imaging device ever built for astronomy. Automated detection and classification of celestial objects will be performed by sophisticated algorithms on high-resolution images to progressively produce an astronomical catalog eventually composed of 20 billion galaxies and 17 billion stars and their associated physical properties. In this article we present an overview of the system currently being constructed to perform data distribution as well as the annual campaigns which reprocess the entire image dataset collected since the beginning of the survey. These processing campaigns will utilize computing and storage resources provided by three Rubin data facilities (one in the US and two in Europe). Each year a Data Release will be produced and disseminated to science collaborations for use in studies comprising four main science pillars: probing dark matter and dark energy, taking inventory of solar system objects, exploring the transient optical sky and mapping the Milky Way. Also presented is the method by which we leverage some of the common tools and best practices used for management of large-scale distributed data processing projects in the high energy physics and astronomy communities. We also demonstrate how these tools and practices are utilized within the Rubin project in order to overcome the specific challenges faced by the Observatory.

The James Webb Space Telescope was not even launched yet when the Astro2020 Decadal Survey American report recommended the development of what is now called the Habitable World Observatory, also mentioned by the Voyage 2050 European report. This future space telescope, at 11 billions dollars and at least 6 m diameter, should allow, around 2040, the characterization of at least 25 exoplanets similar to Earth and orbiting around main sequence stars, with the hope of discovering one where life could have developed. This objective represents a technological challenge since it requires the design of spectro-imagers able to access very high contrasts (10^-8-10^-10) at low angular separations (smaller than 100 mas). This proceeding and the talk it is associated to address various obstacles that remain to be overcome in order to one day allow HabWorld to reach its ultimate performance.

We investigate the structure of admixed neutron stars with a regular hadronic component and a fraction of fermionic self-interacting dark matter. Using two limiting equations of state for the dense baryonic interior, constructed from piecewise generalised polytropes, and an asymmetric self-interacting fermionic dark component, we analyse different scenarios of admixed neutron stars depending on the mass of dark fermions $m_\chi$, interaction mediators $m_\phi$, and self-interacting strengths $g$. We find that the contribution of dark matter to the masses and radii of neutron stars leads to tension with mass estimates of the pulsar J0453+1559, the least massive neutron star, and with the constraints coming from the GW170817 event. We discuss the possibilities of constraining dark matter model parameters $g$ and $y \equiv m_\chi/m_\phi$, using current existing knowledge on neutron star estimations of mass, radius, and tidal deformability, along with the accepted cosmological dark matter freeze-out values and self-interaction cross-section to mass ratio, $\sigma_\mathrm{SI}/m_\chi$, fitted to explain Bullet, Abell, and dwarf galaxy cluster dynamics. By assuming the most restrictive upper limit, $\sigma_\mathrm{SI}/m_\chi < 0.1$ cm$^2$/g, along with dark matter freeze-out range values, the allowed $g$-$y$ region is $0.01 \lesssim g \lesssim 0.1$, with $0.5 \lesssim y \lesssim 200$. For the first time, the combination of updated complementary restrictions is used to set constraints on self-interacting dark matter.

We analyze the validity of optical diagnostic diagrams relying on emission-lines ratios and in the context of classifying Active Galactic Nuclei (AGNs) according to the cosmic metallicity evolution in the redshift range 0 < z < 11.2. In this regard, we fit the results of chemical evolution models (CEMs) to the radial gradients of the N/O abundances ratio derived through direct estimates of electron temperatures (Te-method) in a sample of four local spiral galaxies. This approach allows us to select representative CEMs and extrapolate the radial gradients to the nuclear regions of the galaxies in our sample, inferring in this way the central N/O and O/H abundances. The nuclear abundance predictions for theoretical galaxies from the selected CEMs, at distinct evolutionary stages, are used as input parameters in AGN photoionization models built with the Cloudy code. We found that standard BPT diagnostic diagrams are able to classify AGNs with oxygen abundances 12+logO/H > 8.0 [(Z/Zsolar) > 0.2) preferably found at redshift z > 4. On the other hand, the HeII4685/Hbeta versus [N II]6584/Halpha diagram produces a reliable AGN classification independent of the evolutionary stage of these objects.

Communications to and from a spacecraft undertaking launch-landing interstellar travel at near light speed faces significant challenges. Photon-based communication is significantly impacted by large photon propagation delay and relativistic time dilation. The timing of communications by photon transfer, as measured specifically by local clocks at origin and destination and aboard spacecraft, is analyzed and illustrated for concrete mission scenarios. These include a spacecraft experiencing indefinite constant self-acceleration, and a launch-landing mission, in which a spacecraft experiences constant acceleration for the first half of its cruise phase and a like deceleration for the second half. The origin and destination are assumed to be at rest within a common inertial frame with a wide range of fixed distances separating them. Several typical communication modes are considered, including one-way messaging, two-way message query with an expected response, and the one-way streaming of long program material such as a podcast or video. The local-clock relative timing experienced by the communicating entities including clock images (relation of transmit and receive clocks in one-way communication), the query-response latency (the elapsed time between a query message and reception of a message in response), and the time warping of a streaming program (nonlinear stretching or shrinking of the time axis) are included. In particular, large query-response latency, except for a short interval following launch or before landing, is a severe limit on remote control and social interaction. When photons must travel in the same direction as the spacecraft, communication blackouts strongly limit the periods of time during which communication is possible, and restrict the opportunities for both one-way and two-way communication.

As gravitational wave detectors improve, observations of black hole (BH) mergers will provide the joint distribution of their masses and spins. This will be a critical benchmark to validate formation scenarios. Merging binary BHs formed through isolated binary evolution require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by their surface critical rotation, restricting their angular momentum content at core-collapse. We use stripped star models at low metallicities ($Z_\odot/10$, $Z_\odot/50$ and $Z_\odot/250$) to determine the spins of BHs produced by critically rotating stellar progenitors. To study how such progenitors can arise, we consider their formation through chemically homogeneous evolution (CHE). We use a semianalytical model to study the final spins of CHE binaries, and compare our results against available detailed population synthesis models. We find that above BH masses of $\simeq 25M_\odot$, the dimensionless spin of critically rotating stripped stars ($a = Jc/(GM^2$)) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by binary evolution. CHE can produce binaries where both BHs hit this limit, producing a pile-up at the boundary of the excluded region. Highly spinning BHs arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary BHs in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($\chi_\mathrm{eff}<0.5$) which merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high spin population and the existence of a limit set by critical rotation.

Simultaneous in situ measurements of coronal mass ejections (CMEs), including both plasma and magnetic field, by two spacecraft in radial alignment have been extremely rare. Here, we report on one such CME measured by Solar Orbiter (SolO) and Wind on 2021 November 3--5, while the spacecraft were radially separated by a heliocentric distance of 0.13 au and angularly by only 2.2{\deg}. We focus on the magnetic cloud (MC) part of the CME. We find notable changes in the R and N magnetic field components and in the speed profiles inside the MC between SolO and Wind. We observe a greater speed at the spacecraft further away from the Sun without any clear compression signatures. Since spacecraft are close to each other and computing fast magnetosonic wave speed inside the MC we rule out temporal evolution as the reason on the observed differences suggesting that spatial variations over 2.2{\deg} of the MC structure are at the heart of the observed discrepancies. Moreover, using shock properties at SolO, we forecast an arrival time 2h30 too late for a shock that is just 5h31 away hours from Wind. Predicting the north-south component of the magnetic field at Wind from SolO measurements leads to a relative error of 55 %. These results show that even angular separations as low as 2.2{\deg} (or 0.03 au in arclength) between spacecraft can have a large impact on the observed CME properties, rising up the issue of the resolutions of current CME models and potentially affecting our forecasting capabilities.

Observations of planetary material polluting the atmospheres of white dwarfs are an important probe of the bulk composition of exoplanetary material. Medium- and high-resolution optical and ultraviolet spectroscopy of seven white dwarfs with known circumstellar dust and gas emission are presented. Detections or meaningful upper limits for photospheric absorption lines are measured for: C, O, Na, S, P, Mg, Al, Si, Ca, Ti, Cr, Fe, and Ni. For 16 white dwarfs with known observable gaseous emission discs (and measured photospheric abundances), there is no evidence that their accretion rates differ, on average, from those without detectable gaseous emission. This suggests that, typically, accretion is not enhanced by gas drag. At the effective temperature range of the white dwarfs in this sample (16,000-25,000K) the abundance ratios of elements are more consistent than absolute abundances when comparing abundances derived from spectroscopic white dwarf parameters versus photometric white dwarf parameters. Crucially, this highlights that the uncertainties on white dwarf parameters do not prevent white dwarfs from being utilised to study planetary composition. The abundances of oxygen and silicon for the three hydrogen-dominated white dwarfs in the sample with both optical and ultraviolet spectra differ by 0.62 dex depending on if they are derived from the optical or ultraviolet spectra. This optical/ultraviolet discrepancy may be related to differences in the atmospheric depth of line formation; further investigations into the white dwarf atmospheric modelling are needed to understand this discrepancy.

The Steward Observatory LEO Satellite Photometric Survey is a comprehensive observational survey to characterize the apparent brightness of the Starlink and OneWeb low Earth orbit satellites and evaluate the potential impact on astronomy. We report the results of over 16,000 independent measurements of nearly 2800 individual satellites. In addition to photometry, we also measured the astrometric position of each satellite and evaluated the accuracy of predicting satellite position with the available two-line element sets. The apparent brightness of a satellite seen in the sky is not constant and depends on the Sun-satellite-observer geometry. To capture this, we designed the survey to create an all-geometries set of measurements to fully characterize the brightness of each population of satellites as seen in the sky. We visualize the data with sky-plots that show the correlation of apparent brightness with on-sky position and relative Sun-satellite-observer geometry. The sky-plots show where in the sky the satellites are brightest. In addition to visual magnitudes, we also present two new metrics: the expected photon flux and the effective albedo. The expected photon flux metric assesses the potential impact on astronomy sensors by predicting the flux for a satellite trail in an image from a theoretical 1 m class telescope and sensor. The effective albedo metric assesses where a satellite is more reflective than baseline, which ties to the physical structure of the satellite and indicates the potential for brightness-reducing design changes. We intend to use this methodology and resulting data to inform the astronomy community about satellite brightness.

A detailed study of stellar populations in Milky Way (MW) satellite galaxies remains an observational challenge due to their faintness and fewer spectroscopically confirmed member stars. We use unsupervised machine learning methods to identify new members for nine nearby MW satellite galaxies using Gaia data release-3 (Gaia DR3) astrometry and the Dark Energy Survey (DES) and the DECam Local Volume Exploration Survey (DELVE) photometry. Two density-based clustering algorithms, DBSCAN and HDBSCAN, have been used in the four-dimensional astrometric parameter space to identify member stars belonging to MW satellite galaxies. Our results indicate that we can recover more than 80% of the known spectroscopically confirmed members in most of the satellite galaxies and also reject 95-100% of spectroscopic non-members. We have also added many new members using this method. We compare our results with previous studies that also use photometric and astrometric data and discuss the suitability of density-based clustering methods for MW satellite galaxies

PRAIA - Package for the Reduction of Astronomical Images Automatically - is a suite of photometric and astrometric tasks designed to cope with huge amounts of heterogeneous observations with fast processing, no human intervention, minimum parametrization and yet maximum possible accuracy and precision. It is the main tool used to analyse astronomical observations by an international collaboration involving Brazilian, French and Spanish researchers under the Lucky Star umbrella for Solar System studies. Here, we focus on the concepts of differential aperture photometry and digital coronagraphy underneath PRAIA, used in the reduction of stellar occultations, rotational light curves, mutual phenomena and natural satellite observations. We highlight novelties developed by us and never before reported in the literature, which significantly enhance the precision and automation of photometry and digital coronagraphy, such as: a) PRAIA's pixelized aperture photometry (PCAP); b) fully automatic object detection and aperture determination (BOIA); c) better astrometry improving the aperture and coronagraphy centre, including the new Photogravity Center Method besides circular and elliptical Gaussian and Lorentzian generalized profiles; d) coronagraphy of faint objects close to bright ones and vice-versa; e) use of elliptical rings for the coronagraphy of elongated profiles; f) refined quartile ring statistics; g) multiprocessing image capabilities for faster computation speed. We give examples showing the photometry performance, discuss the advantages of PRAIA over other popular packages, and point out the uniqueness of its digital coronagraphy in comparison with other coronagraphy tools. Besides Solar System works, PRAIA can also be used in the differential photometry and digital coronagraphy of any astrophysical observations. PRAIA codes are publicly available at: https://ov.ufrj.br/en/PRAIA/.

Unsupervised learning, a branch of machine learning that can operate on unlabelled data, has proven to be a powerful tool for data exploration and discovery in astronomy. As large surveys and new telescopes drive a rapid increase in data size and richness, these techniques offer the promise of discovering new classes of objects and of efficient sorting of data into similar types. However, unsupervised learning techniques generally require feature extraction to derive simple but informative representations of images. In this paper, we explore the use of self-supervised deep learning as a method of automated representation learning. We apply the algorithm Bootstrap Your Own Latent (BYOL) to Galaxy Zoo DECaLS images to obtain a lower dimensional representation of each galaxy. We briefly validate these features using a small supervised classification problem. We then move on to apply an automated clustering algorithm, demonstrating that this fully unsupervised approach is able to successfully group together galaxies with similar morphology. The same features prove useful for anomaly detection, where we use the framework astronomaly to search for merger candidates. Finally, we explore the versatility of this technique by applying the exact same approach to a small radio galaxy dataset. This work aims to demonstrate that applying deep representation learning is key to unlocking the potential of unsupervised discovery in future datasets from telescopes such as the Vera C. Rubin Observatory and the Square Kilometre Array.

Pulsar timing arrays seek and study gravitational waves (GWs) through the angular two-point correlation function of timing residuals they induce in pulsars. The two-point correlation function induced by the standard transverse-traceless GWs is the famous Hellings-Downs curve, a function only of the angle between the two pulsars. Additional polarization modes (vector/scalar) that may arise in alternative-gravity theories have different angular correlation functions. Furthermore, anisotropy, linear, or circular polarization in the stochastic GW background gives rise to additional structure in the two-point correlation function that cannot be written simply in terms of the angular separation of the two pulsars. In this paper, we provide a simple formula for the most general two-point correlation function -- or overlap reduction function (ORF) -- for a gravitational-wave background with an arbitrary polarization state, possibly containing anisotropies in its intensity and polarization (linear/circular). We provide specific expressions for the ORFs sourced by the general-relativistic transverse-traceless GW modes as well as vector (or spin-1) modes that may arise in alternative-gravity theories.

The launch of JWST has ushered in a new era of high precision infrared astronomy, allowing us to probe nearby white dwarfs for cold dust, exoplanets, and tidally heated exomoons. While previous searches for these exoplanets have successfully ruled out companions as small as 7-10 Jupiter masses, no instrument prior to JWST has been sensitive to the likely more common sub-Jovian mass planets around white dwarfs. In this paper, we present the first multi-band photometry (F560W, F770W, F1500W, F2100W) taken of WD 2149+021 with the Mid-Infrared Instrument (MIRI) on JWST. After a careful search for both resolved and unresolved planets, we do not identify any compelling candidates around WD 2149+021. Our analysis indicates that we are sensitive to companions as small as ~0.34 MJup outwards of 1."263 (28.3 au) and ~0.64 MJup at the innermost working angle (0."654, 14.7 au) with 5 sigma confidence, placing significant constraints on any undetected companions around this white dwarf. The results of these observations emphasize the exciting future of sub-Jovian planet detection limits by JWST, which can begin to constrain how often these planets survive their host stars evolution.

We argue that the broad-band observations of the brightest-of-all-time (BOAT) GRB 221009A reveal a physical picture involving two jet components: a narrow ($\sim 0.6$ degree half opening angle) pencil-beam jet that has a Poynting-flux-dominated jet composition, and a broader matter-dominated jet with an angular structure. We discuss various observational evidence that supports such a picture. To treat the problem, we develop an analytical structured jet model for both forward and reverse shock emission from the matter dominated structured jet wing during the deceleration phase. We discuss the physical implications of such a two-component jet configuration for this particular burst and for GRBs in general. We argue that some bright X-ray flares could be similar narrow jets viewed slightly outside the narrow jet cone and that narrow jets may exist in many more GRBs without being detected.

We build on our previous work involving smoothed particle hydrodynamic simulations of Be stars, by using the model that exhibited disc tearing as input into the three-dimensional Monte Carlo radiative transfer code HDUST to predict observables from a variety of viewing angles throughout the disc tearing process. We run one simulation at the start of each orbital period from 20 to 72 orbital periods, which covers two complete disc tearing events. The resulting trends in observables are found to be dependent on the relative position of the observer and the tearing disc. The $\rm H\alpha$ equivalent width, $V$ magnitude, and polarization can all increase or decrease in any combination depending on the viewpoint of the observer. The $\rm H\alpha$ line profile also displays changes in strength and peak separation throughout the tearing process. We show how the outer disc of the torn system can have a large effect on the $\rm H\alpha$ line profile, and also contributes to a wavelength-dependent polarization position angle, resulting in a similar sawtooth shape to the polarization percentage. Finally, we compare our predictions to Pleione (28 Tau) where evidence has suggested that a disc tearing event has occurred in the past. We find that our tearing disc model can broadly match the trends seen in Pleione's observables, as well as produce the two-component $\rm H\alpha$ lines observed in Pleione. This is the strongest evidence, thus far, of Pleione's disc having indeed experienced a tearing event.

A planet's Lyman-{\alpha} (Ly{\alpha}) emission is sensitive to its thermospheric structure. Here, we report joint Hubble Space Telescope (HST) and Cassini cross-calibration observations of the Saturn Ly{\alpha} emission made two weeks before the Cassini grand finale. To investigate the long-term Saturn Ly{\alpha} airglow observed by different ultraviolet instruments, we cross-correlate their calibration, finding that while the official Cassini/UVIS sensitivity should be lowered by ~75%, the Voyager 1/UVS sensitivities should be enhanced by ~20% at the Ly{\alpha} channels. This comparison also allowed us to discover a permanent feature of the Saturn disk Ly{\alpha} brightness that appears at all longitudes as a brightness excess (Ly{\alpha} bulge) of ~30% (~12{\sigma}) extending over the latitude range ~5-35N compared to the regions at equator and ~60N. This feature is confirmed by three distinct instruments between 1980 & 2017 in the Saturn north hemisphere. To analyze the Ly{\alpha} observations, we use a radiation transfer (RT) model of resonant scattering of solar and interplanetary Ly{\alpha} photons, and a latitude-dependent photochemistry model of the upper atmosphere constrained by occultation and remote-sensing observations. For each latitude, we show that the Ly{\alpha} observations are sensitive to the temperature profile in the upper stratosphere and lower thermosphere, thus providing useful information in a region of the atmosphere that is difficult to probe by other means. In the Saturn Ly{\alpha} bulge region, at latitudes between ~5 to ~35{\deg}, the observed brightening and line broadening support seasonal effects, variation of the temperature vertical profile, and potential superthermal atoms that require confirmation.

Baryon Acoustic Oscillation (BAO) observations offer a robust method for measuring cosmological expansion. However, the BAO signal in a sample of galaxies can be diluted and shifted by interlopers - galaxies that have been assigned the wrong redshifts. Because of the slitless spectroscopic method adopted by the Roman and Euclid space telescopes, the galaxy samples resulting from single line detections will have relatively high fractions of interloper galaxies. Interlopers with a small displacement between true and false redshift have the strongest effect on the measured clustering. In order to model the BAO signal, the fraction of such interlopers and their clustering need to be accurately known. We introduce a new method to self-calibrate these quantities by shifting the contaminated sample towards or away from us along the line of sight by the interloper offset, and measuring the cross-correlations between these shifted samples. The contributions from the different components are shifted in scale in this cross-correlation compared to the auto-correlation of the contaminated sample, enabling the decomposition and extraction of the component terms. We demonstrate the application of the method using numerical simulations and show that an unbiased BAO measurement can be extracted. Unlike previous attempts to model the effects of contaminants, self-calibration allows us to make fewer assumptions about the form of the contaminants such as their bias.

Alternative cosmological models have been proposed to alleviate the tensions reported in the concordance cosmological model, or to explain the current accelerated phase of the universe. One way to distinguish between General Relativity and modified gravity models is using current astronomical data to measure the growth index $\gamma$, a parameter related to the growth of matter perturbations, which behaves differently in different metric theories. We propose a model independent methodology for determining $\gamma$, where our analyses combine diverse cosmological data sets, namely $\{ f(z_i) \}$, $\{ [f\sigma_8](z_i) \}$, and $\{ H(z_i) \}$, and use Gaussian Processes, a non-parametric approach suitable to reconstruct functions. This methodology is a new consistency test for $\gamma$ constant. Our results show that, for the redshift interval $0 < z < 1$, $\gamma$ is consistent with the constant value $\gamma = 0.55$, expected in General Relativity theory, within $2 \sigma$ confidence level (CL). Moreover, we find $\gamma(z=0)$ = $0.311 \pm 0.229$ and $\gamma(z=0) = 0.336 \pm 0.107$ for the reconstructions using the $\{ f(z_i) \}$ and $\{ [f\sigma_8](z_i) \}$ data sets, respectively, values that also agree at a 2$\sigma$ CL with $\gamma = 0.55$. Our main results are a third way in light of the current discussion in literature that point out to some possible evidence for the growth index evolution.

A quarter of a century has passed since the observing technique of integral field spectroscopy (IFS) was first applied to planetary nebulae (PNe). Progress after the early experiments was relatively slow, mainly because of the limited field-of-view (FoV) of first generation instruments.With the advent of MUSE at the ESO Very Large Telescope, this situation has changed. MUSE is a wide field-of-view, high angular resolution, one-octave spanning optical integral field spectrograph with high throughput. Its major science mission has enabled an unprecedented sensitive search for Ly{\alpha} emitting galaxies at redshift up to z=6.5. This unique property can be utilized for faint objects at low redshift as well. It has been demonstrated that MUSE is an ideal instrument to detect and measure extragalactic PNe with high photometric accuracy down to very faint magnitudes out to distances of 30 Mpc, even within high surface brightness regions of their host galaxies. When coupled with a differential emission line filtering (DELF) technique, MUSE becomes far superior to conventional narrow-band imaging, and therefore MUSE is ideal for accurate Planetary Nebula Luminosity Function (PNLF) distance determinations. MUSE enables the PNLF to become a competitive tool for an independent measure of the Hubble constant, and stellar population studies of the host galaxies that present a sufficiently large number of PNe.

Cosmic rays are energetic charged particles from extraterrestrial sources, with the highest energy events thought to come from extragalactic sources. Their arrival is infrequent, so detection requires instruments with large collecting areas. In this work, we report the detection of an extremely energetic particle recorded by the surface detector array of the Telescope Array experiment. We calculate the particle's energy as 244 +- 29 (stat.) +51,-76 (syst.) exa-electron volts (~40 joules). Its arrival direction points back to a void in the large-scale structure of the Universe. Possible explanations include a large deflection by the foreground magnetic field, an unidentified source in the local extragalactic neighborhood or an incomplete knowledge of particle physics.

Abundance determinations in planetary nebulae (PNe) are crucial for understanding stellar evolution and the chemical evolution of the host galaxy. We discuss the complications involved when the presence of a metal-rich phase is suspected in the nebula. We demonstrate that the presence of a cold region emitting mainly metal recombination lines necessitates a detailed treatment to obtain an accurate assessment of the enrichment of this cold gas phase.

Tidal streams and stellar shells are naturally formed in galaxy interactions and mergers. The Giant Stellar Stream (GSS), the North-East (NE), and Western (W) stellar shelves observed in Andromeda galaxy (M31) are examples of these structures and were formed through the merger of M31 and a satellite galaxy. Recent observational papers have provided strong evidence that the shells and GSS originate from a single progenitor. In this paper, we investigate the formation of these two stellar shelves and the detailed nature of their relationship to the GSS. We present numerical simulations of tidal disruption of a satellite galaxy assuming that it is a progenitor of the GSS and the shell system. We represent the progenitor as a dwarf spheroidal galaxy with the stellar mass of $10^{9} M_{\odot}$ and evolve its merger with M31 for 3 Gyrs to reproduce the chemodynamical properties of the NE and W shelves. We find that an initial metallicity of the progenitor with a negative radial gradient of $\Delta$ FeH = -0.3 $\pm$ 0.2, successfully reproduces observed metallicities of the NE, W shelves, and the GSS, showing that all these structures can originate from the same merger event.

We present the Bidimensional Exploration of the warm-Temperature Ionised gaS (BETIS) project, designed for the spatial and spectral study of the Diffuse Ionised Gas (DIG) in a selection of 265 spiral galaxies observed with MUSE, of different morphologies and an average spatial resolution of 274 pc (z<0.06). Our main objective is to study the ionisation mechanisms at play within the DIG by analysing the distribution of various species in the optical spectra. We introduce a new methodology for defining the DIG, employing an innovative adaptive binning technique on the observed datacube based on the spectroscopic signal-to-noise ratio (SN) of the [S II] line, in order to increase the SN of the rest of the lines. Subsequently, we create a DIG mask by subtracting the emission associated with both bright and faint H II regions. We also examine the suitability of using H{\alpha} equivalent width (EW(H{\alpha})) as a proxy for defining the DIG and its associated ionisation regime. Notably, for EW(H{\alpha})< 3{\AA}, the measured value is contingent on the chosen population synthesis technique performed. Our analysis in a showcase subsample reveals a consistent cumulative DIG fraction across them, averaging around 40%-70%. The average radial distribution of the main line ratios are enhanced in the DIG regimes (up to 0.2 dex), and follows similar trends between the DIG and the H II regions, as well as H{\alpha} surface brightness, indicating a correlation between the ionisation of these species in both the DIG and the H II regions. The DIG loci in the BPT diagrams are found within the line ratios that correspond to photoionisation due to the star formation. There is an offset corresponding to ionisation due to fast shocks, primarily due to the contribution of the two Seyfert galaxies in our subsample, which closely align with models of ionisation from fast shocks, thus mimicking the DIG emission.

We present highly resolved and sensitive imaging of the five nearby massive spiral galaxies (with rotation velocities $\rm > 300 km s^{-1}$) observed by the UltraViolet Imaging Telescope onboard India's multi-wavelength astronomy satellite ASTROSAT, along with other archival observations. These massive spirals show a far-ultraviolet star formation rate in the range of $\sim$$\rm 1.4$$-$$\rm13.7 M_{\odot} yr^{-1}$ and fall in the `Green Valley' region with a specific star formation rate within $\sim$$\rm10^{-11.5}$$-$$\rm10^{-10.5}yr^{-1}$. Moreover, the mean star formation rate density of the highly resolved star-forming clumps of these objects are in the range $\rm 0.011$$-$$\rm 0.098 M_{\odot} yr^{-1}kpc^{-2}$, signifying localised star formation. From the spectral energy distributions, under the assumption of a delayed star formation model, we show that the star formation of these objects had peaked in the period of $\sim$ $\rm0.8$$-$$\rm2.8$ Gyr after the `Big Bang' and the object that has experienced the peak sooner after the `Big Bang' show relatively less star-forming activity at $\rm\it{z}$$\sim$0 and falls below the main-sequence relation for a stellar content of $\rm \gtrsim 10^{11} M_{\odot}$. We also show that these objects accumulated much of their stellar mass in the early period of evolution with $\sim31$$-$$42$ per cent of the total stellar mass obtained in a time of $(1/16)$$-$$(1/5)^{\rm th}$ the age of the Universe. We estimate that these massive objects convert their halo baryons into stars with efficiencies falling between $\sim 7-31$ percent.

Using 3D particle-in-cell simulation, we characterize energy conversion, as a function of guide magnetic field, in a thin current sheet in semirelativistic plasma, with relativistic electrons and subrelativistic protons. There, magnetic reconnection, the drift-kink instability (DKI), and the flux-rope kink instability all compete and interact in their nonlinear stages to convert magnetic energy to plasma energy. We compare fully 3D simulations with 2D in two different planes to isolate reconnection and DKI effects. In zero guide field, these processes yield distinct energy conversion signatures: ions gain more energy than electrons in 2Dxy (reconnection), while the opposite is true in 2Dyz (DKI), and the 3D result falls in between. The flux-rope instability, which occurs only in 3D, allows more magnetic energy to be released than in 2D, but the rate of energy conversion in 3D tends to be lower. Increasing the guide magnetic field strongly suppresses DKI, and in all cases slows and reduces the overall amount of energy conversion; it also favors electron energization through a process by which energy is first stored in the motional electric field of flux ropes before energizing particles. Understanding the evolution of the energy partition thus provides insight into the role of various plasma processes, and is important for modeling radiation from astrophysical sources such as accreting black holes and their jets.

We present sensitive, high angular resolution Jansky Very Large Array observations made in 2014 at 1.50 GHz toward the field of Tycho's supernova remnant. We detect a total of 36 compact sources in a field with radius of 13 arcmin. This number is consistent with the expected number of background sources. We use older observations made with the classic Very Large Array to compare with the 2014 observations and search for sources showing large proper motions that could be related to the donor companion of the exploding white dwarf that produced the supernova in 1572. The comparison of the positions for the two sets of observations does not show sources with large proper motions and supports the conclusion that all sources detected are extragalactic and unrelated to the supernova field.

Radio Frequency Interference (RFI) detection and mitigation is critical for enabling and maximising the scientific output of radio telescopes. The emergence of machine learning methods capable of handling large datasets has led to their application in radio astronomy, particularly in RFI detection. Spiking Neural Networks (SNNs), inspired by biological systems, are well-suited for processing spatio-temporal data. This study introduces the first application of SNNs to an astronomical data-processing task, specifically RFI detection. We adapt the nearest-latent-neighbours (NLN) algorithm and auto-encoder architecture proposed by previous authors to SNN execution by direct ANN2SNN conversion, enabling simplified downstream RFI detection by sampling the naturally varying latent space from the internal spiking neurons. We evaluate performance with the simulated HERA telescope and hand-labelled LOFAR dataset that the original authors provided. We additionally evaluate performance with a new MeerKAT-inspired simulation dataset. This dataset focuses on satellite-based RFI, an increasingly important class of RFI and is, therefore, an additional contribution. Our SNN approach remains competitive with the original NLN algorithm and AOFlagger in AUROC, AUPRC and F1 scores for the HERA dataset but exhibits difficulty in the LOFAR and MeerKAT datasets. However, our method maintains this performance while completely removing the compute and memory-intense latent sampling step found in NLN. This work demonstrates the viability of SNNs as a promising avenue for machine-learning-based RFI detection in radio telescopes by establishing a minimal performance baseline on traditional and nascent satellite-based RFI sources and is the first work to our knowledge to apply SNNs in astronomy.

How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core-collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios ($q=M_2/M_1$ = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion induced collapse (AIC), rapidly rotating supernova (RSN), black hole formation, and gamma ray burst (GRB). We then conduct Binary Population Synthesis calculations for 12 models, varying the initial parameter distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate $R_{\rm SN} = (1.14$--$1.57) \times10^{-2} \, {\rm yr}^{-1}$ which is consistent with observations. We find the rates of AIC, RSN, and GRB to be $\sim 1/100$ the rate of regular supernovae. Our estimated GRB rates are higher than the observed long GRB rate, but very close to the low luminosity GRB rate. Furthering binary modeling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars.

Recently, several new solar radio telescopes have been put into operation and provided spectral-imaging observations with much higher resolutions in decimeter (dm) and centimeter (cm) wavelengths. These telescopes include the Mingantu Spectral Radioheliograph (MUSER, at frequencies of 0.4 - 15 GHz), the Expanded Owens Valley Solar Array (EOVSA, at frequencies of 1 - 18 GHz), and the Siberian Radio Heliograph (SRH, at frequencies of 3 - 24 GHz). These observations offer unprecedented opportunities to study solar physics and space weather, especially to diagnose the coronal magnetic fields, reveal the basic nature of solar eruptions and the related non-thermal energy release, particle accelerations and propagation, and the related emission mechanisms. These results might be the important input to the space weather modeling for predicting the occurrence of disastrous powerful space weather events. In order to provide meaningful reference for other solar physicists and space weather researchers, this paper mainly focus on discussing the potential scientific problems of solar radio spectral-imaging observations in dm-cm wavelengths and its possible applications in the field of space weather. These results will provide a helpful reference for colleagues to make full use of the latest and future observation data obtained from the above solar radio telescopes.

Binaries containing compact objects, if viewed close to edge on, can produce periodic brightening events under certain conditions on the masses, radii, and binary separation. Such flares are caused by one object gravitational lensing the other, in what is known as self-lensing flares. We present a simulation tool that efficiently reproduces the main features of self-lensing flares and facilitates a detection sensitivity analysis for various sky surveys. We estimate the detection prospects for a handful of representative surveys when searching for systems of either two white dwarfs, or a white dwarf with other compact objects, i.e., neutron stars and black holes. We find only a marginal ability to detect such systems in existing surveys. However, we estimate many such systems could be detectable by surveys in the near future, including the Vera Rubin observatory. We provide a quantitative analysis of the detectability of double-compact object self-lensing flares across the landscape of system parameters, and a qualitative discussion of survey and followup approaches to distinguish such flares from confounding events, such as stellar flares, satellite glints, and cosmic rays. We estimate 0.3, 3 and 247 double white dwarf systems could be detected by TESS, ZTF, and LSST, respectively. A similar number of systems with a neutron star or black hole companion could be detected, but we caution that the number densities of such binaries is model dependent and so our detection estimates. Such binaries can be used to constrain models of the end states of binary evolution.

Shock breakout emission is light that arises when a shockwave, generated by core-collapse explosion of a massive star, passes through its outer envelope. Hitherto, the earliest detection of such a signal was at several hours after the explosion, though a few others had been reported. The temporal evolution of early light curves should reveal insights into the shock propagation, including explosion asymmetry and environment in the vicinity, but this has been hampered by the lack of multiwavelength observations. Here we report the instant multiband observations of a type II supernova (SN 2023ixf) in the galaxy M101 (at a distance of 6.85+/-0.15 Mpc), beginning at about 1.4 hours after the explosion. The exploding star was a red supergiant with a radius of 450 solar radii. The light curves evolved rapidly, on timescales of 1-2 hours, and appeared unusually fainter and redder than predicted by models within the first few hours, which we attribute to an optically thick dust shell before it was disrupted by the shockwave. We infer that the breakout and perhaps the distribution of the surrounding dust were not spherically symmetric.

Current cosmological tensions show that it is crucial to test the predictions from the canonical $\Lambda$CDM paradigm at different cosmic times. One very appealing test of structure formation in the universe is the growth rate of structure in our universe $f$, usually parameterized via the growth index $\gamma$, with $f\equiv \Omega_m(a)^\gamma$ and $\gamma \simeq 0.55$ in the standard $\Lambda$CDM case. Recent studies have claimed a suppression of the growth of structure from a variety of cosmological observations, characterized by $\gamma>0.55$. By employing different self-consistent growth parameterizations schemes, we show here that $\gamma<0.55$, obtaining instead \emph{an enhanced growth of structure today}. This preference reaches the $3\sigma$ significance using Cosmic Microwave Background observations, Supernova Ia and Baryon Acoustic Oscillation measurements. The addition of Cosmic Microwave Background lensing data relaxes such a preference to the $2\,\sigma$ level, since a larger lensing effect can always be compensated with a smaller structure growth, or, equivalently, with $\gamma>0.55$. We have also included the lensing amplitude $A_{\rm L}$ as a free parameter in our data analysis, showing that the preference for $A_{\rm L}>1$ still remains, except for some particular parameterizations when lensing observations are included. We also not find any significant preference for a multipole dependence of $A_{\rm L}$. To further reassess the effects of a non-standard growth, we have computed by means of N-body simulations the dark matter density fields, the dark matter halo mass functions and the halo density profiles for different values of $\gamma$. Future observations from the Square Kilometer Array, reducing by a factor of three the current errors on the $\gamma$ parameter, could finally settle the issue.

Upcoming galaxy surveys provide the necessary sensitivity to measure gravitational redshift, a general relativistic effect that generates a dipole in galaxy clustering data. Here, we study the constraining power of gravitational redshift within the framework of the effective theory of interacting dark energy. This formalism describes linear cosmological perturbations in scalar-tensor theories of gravity with a limited number of free functions, and allows each particle species to be coupled differently to the gravitational sector. In this work, we focus on Horndeski theories with a non-minimal coupling of dark matter to the scalar degree of freedom, yielding a breaking of the weak equivalence principle for this cosmic component, a scenario that is yet untested. We show that the dipole generated by gravitational redshift significantly breaks degeneracies and tightens the constraints on the parameters of the effective theory compared to the standard redshift-space distortion analysis solely based on the even multipoles in the galaxy correlation function, with an improvement of up to $\sim 50\%$. We make the Python package \texttt{EF-TIGRE} (\textit{Effective Field Theory of Interacting dark energy with Gravitational REdshift}) developed for this work publicly available.

Solar Orbiter's four in-situ instruments have recorded numerous energetic electron events at heliocentric distances between 0.5 and 1au. We analyse energetic electron fluxes, spectra, pitch angle distributions, associated Langmuir waves, and type III solar radio bursts for 3 events to understand what causes modifications in the electron flux and identify the origin and characteristics of features observed in the electron spectrum. We investigate what electron beam properties and solar wind conditions are associated with Langmuir wave growth and spectral breaks in the electron peak flux as a function of energy. We observe velocity dispersion and quasilinear relaxation in the electron flux caused by the resonant wave-particle interactions in the deca-keV range, at the energies at which we observe breaks in the electron spectrum, co-temporal with the local generation of Langmuir waves. We show, via the evolution of the electron flux at the time of the event, that these interactions are responsible for the spectral signatures observed around 10 and 50keV, confirming the results of simulations by Kontar & Reid (2009). These signatures are independent of pitch angle scattering. Our findings highlight the importance of using overlapping FOVs when working with data from different sensors. In this work, we exploit observations from all in-situ instruments to address, for the first time, how the energetic electron flux is modified by the beam-plasma interactions, and results into specific features to appear in the local spectrum. Our results, corroborated with numerical simulations, can be extended to a wider range of heliocentric distances.

Cosmologists aim to uncover the underlying cosmological model governing the formation and evolution of the Universe. One approach is through studying the large-scale structure (LSS) traced by galaxy redshift surveys. In this paper, we explore clustering and covariance errors of BOSS and eBOSS surveys in configuration and Fourier space with a new generation of galaxy lightcones. We create 16 lightcones using the UCHUU simulation: a 2$h^{-1}$Gpc $N$-body simulation tracking 2.1 trillion dark matter particles within a Planck-$\Lambda$CDM cosmology. Simulation's (sub)halos are populated with Luminous red galaxies (LRGs) using the subhalo abundance matching. For estimating covariance errors, we generate 5,040 GLAM-UCHUU LRG lightcones based on GLAM $N$-body simulations. LRGs are included using halo occupation distribution. Our simulated lightcones reproduce BOSS/eBOSS clustering statistics on scales from redshifts 0.2 to 1.0, from 2 to 150$h^{-1}$Mpc, and from 0.005 to 0.7$h$Mpc$^{-1}$, in configuration and Fourier space, respectively. We analyse stellar mass and redshift effects on clustering and bias, revealing consistency with data and noting an increasing bias factor with redshift. Our investigation leads us to the conclusion that the Planck-$\Lambda$CDM cosmology accurately explains the observed LSS. Furthermore, we compare our GLAM-UCHUU LRG lightcones with MD-PATCHY and EZMOCK, identifying large deviations from observations within 20$h^{-1}$Mpc. We examine covariance matrices, finding that our data estimated errors are higher than those previously reported, carrying significant implications for cosmological parameter inferences. Lastly, we explore cosmology's impact on galaxy clustering. Our results suggest that, given the current level of uncertainties, we are unable to distinguish models with and without massive neutrino effects on LSS.

The geographic South Pole provides unique opportunities to study cosmic particles in the Southern Hemisphere. It represents an optimal location to deploy large-scale neutrino telescopes in the deep Antarctic ice, such as AMANDA or IceCube. In both cases, the presence of an array, constructed to observe extensive air showers, enables hybrid measurements of cosmic rays. While additional neutron monitors can provide information on solar cosmic rays, large detector arrays, like SPASE or IceTop, allow for precise measurements of cosmic rays with energies above several $100\,\rm{TeV}$. In coincidence with the signals recorded in the deep ice, which are mostly due to the high-energy muons produced in air showers, this hybrid detector setup provides important information about the nature of cosmic rays. In this review, we will discuss the historical motivation and developments towards measurements of cosmic rays at the geographic South Pole and highlight recent results reported by the IceCube Collaboration. We will emphasize the important contributions by Thomas K. Gaisser and his colleagues that ultimately led to the rich Antarctic research program which today provides crucial insights into cosmic-ray physics.

We present a systematic search for periodic X-ray sources in 10 Galactic globular clusters (GCs) utilizing deep archival Chandra observations. By applying the Gregory-Loredo algorithm, we detect 28 periodic signals among 27 independent X-ray sources in 6 GCs, which include 21 newly discovered ones in the X-ray band. The remaining 4 GCs exhibit no periodic X-ray sources, mainly due to a relatively lower sensitivity of the data. Through analysis of their X-ray timing and spectral properties, complemented with available optical and ultraviolet information, we identify 21 of these periodic sources as cataclysmic variables (CVs). Combining with 11 periodic CVs in 47 Tuc similarly identified in the X-ray band, we compile the most comprehensive sample to date of GC CVs with a probable orbital period. The scarcity of old, short-period CVs in GCs compared to the Galactic inner bulge and solar neighborhood, can be attributed to both a selection effect favouring younger, dynamically-formed systems and the hindrance of CV formation through primordial binary evolution by stellar dynamical interactions common to the GC environment. Additionally, we identify a significant fraction of the GC CVs, most with an orbital period below or within the CV period gap, as probable magnetic CVs, but in the meantime there is a deficiency of luminous intermediate polars in the GC sample compared to the solar neighborhood.

The highly-filamented nature of the coronal plasma significantly influences dynamic processes in the corona such as magnetohydrodynamic waves and oscillations. Fast magnetoacoustic waves, guided by coronal plasma non-uniformities, exhibit strong geometric dispersion, forming quasi-periodic fast-propagating (QFP) wave trains. QFP wave trains are observed in extreme-ultraviolet imaging data and indirectly in microwaves and low-frequency radio, aiding in understanding the magnetic connectivity, energy, and mass transport in the corona. However, measuring the field-aligned group speed of QFP wave trains, as a key parameter for seismological analysis, is challenging due to strong dispersion and associated rapid evolution of the wave train envelope. We demonstrate that the group speed of QFP wave trains formed in plane low-$\beta$ coronal plasma non-uniformities can be assessed through the propagation of the wave train's effective centre of mass, referred to as the wave train's centroid speed. This centroid speed, as a potential observable, is shown empirically to correspond to the group speed of the most energetic Fourier harmonic in the wave train. The centroid speed is found to be almost insensitive to the waveguide density contrast with the ambient corona, and to vary with the steepness of the transverse density profile. The discrepancy between the centroid speed as the group speed measure and the phase speed at the corresponding wavelength is shown to reach 70\%, which is crucial for the energy flux estimation and interpretation of observations.

Scalar induced gravitational waves contribute to the cosmological gravitational wave background. They can be related to the primordial density power spectrum produced towards the end of inflation and therefore are a convenient new tool to constrain models of inflation. These waves are sourced by terms quadratic in perturbations and hence appear at second order in cosmological perturbation theory. While the focus of research so far was on purely scalar source terms we also study the effect of including first order tensor perturbations as an additional source. This gives rise to two additional source terms: a term quadratic in the tensor perturbations and a cross term involving mixed scalar and tensor perturbations. We present full analytical expressions for the spectral density of these new source terms and discuss their general behaviour. To illustrate the generation mechanism we study two toy models containing a peak on small scales. For these models we show that the scalar-tensor contribution becomes non-negligible compared to the scalar-scalar contribution on smaller scales. We also consider implications for future gravitational wave surveys.

We present imaging and spectroscopic diagnostics of a long filament during its formation with the observations from the Chinese H$\alpha$ Solar Explorer and Solar Dynamics Observatory. The seed filament first appeared at about 05:00 UT on 2022 September 13. Afterwards, it grew gradually and connected to another filament segment nearby, building up a long filament at about 20:00 UT on the same day. The CHASE H$\alpha$ spectra show an obvious centroid absorption with mild broadening at the main spine of the long filament, which is interpreted as the evidence of filament material accumulation. More interestingly, near the footpoints of the filament, persistent redshifts have been detected in the H$\alpha$ spectra during the filament formation, indicating continuous drainage of filament materials. Furthermore, through inspecting the extreme ultraviolet images and magnetograms, it is found that EUV jets and brightenings appeared repeatedly at the junction of the two filament segments, where opposite magnetic polarities converged and canceled to each other continuously. These results suggest the occurrence of intermittent magnetic reconnection that not only connects magnetic structures of the two filament segments but also supplies cold materials for the filament channel likely by the condensation of injected hot plasma, even though a part of cold materials fall down to the filament footpoints at the same time.

Global 21-cm experiments are built to study the evolution of the Universe between the cosmic dawn and the epoch of reionisation. FlexKnot is a function parameterised by freely moving knots stringed together by splines. Adopting the FlexKnot function as the signal model has the potential to separate the global 21-cm signal from the foregrounds and systematics while being capable of recovering the crucial features given by theoretical predictions. In this paper, we implement the FlexKnot method by integrating twice over a function of freely moving knots interpolated linearly. The function is also constrained at the lower frequencies corresponding to the dark ages by theoretical values. The FlexKnot model is tested in the framework of the realistic data analysis pipeline of the REACH global signal experiment using simulated antenna temperature data. We demonstrate that the FlexKnot model performs better than existing signal models, e.g. the Gaussian signal model, at reconstructing the shape of the true signals present in the simulated REACH data, especially for injected signals with complex structures. The capabilities of the FlexKnot signal model is also tested by introducing various systematics and simulated global signals of different types. These tests show that four to five knots are sufficient to recover the general shape of most realistic injected signals, with or without sinusoidal systematics. We show that true signals whose absorption trough of amplitude between 120 to 450 mK can be well recovered with systematics up to about 50 mK.

Close binary central stars of planetary nebulae (PNe) must have formed through a common envelope evolution during the giant phase experienced by one of the stars. Transfer of the angular momentum from the binary system to the envelope leads to the shortening of the binary separations from the radius of red giant to the radius of the order of few tenths of AU. Thus, close binary central stars of planetary nebulae are laboratories to study the common envelope phase of evolution. The close binary fraction in the Galaxy has been measured in various sky surveys, but the close binary fraction is not yet well constrained for the Magellanic Clouds, and our results may help the study of common envelope evolution in low-metallicity environments. This paper presents a continuation of our study of variability in the Magellanic Cloud planetary nebulae on the basis of data from the OGLE survey. Previously, we had analysed the OGLE data in the Small Magellanic Cloud. Here, the study is extended to the Large Magellanic Cloud (LMC). In this paper we search for close binary central stars with the aim to constrain the binary fraction and period distribution in the LMC. We identified 290 counterparts of PNe in the LMC in the I-band images from the OGLE-III and OGLE-IV surveys. However, the light curves of ten objects were not accessible in the OGLE database, and thus we analysed the time series photometry of 280 PNe. In total, 32 variables were found, but 5 of them turned out to be foreground objects. Another 18 objects show irregular or regular variability that is not attributable to the binarity of their central stars. Their status and the nature of their variability will be verified in the follow-up paper. Nine binary central stars of PNe with periods between 0.24 and 23.6 days were discovered. The obtained fraction for the LMC PNe is 3.3^(+2.6)_(-1.6)% without correcting for incompleteness.

We explore the possibility that the N-rich young proto-galaxy GN-z11 recently observed at z=10.6 by the James Webb Space Telescope is the result of the formation of second generation stars from pristine gas and Asymptotic Giant Branch (AGB) ejecta in a massive globular cluster or nuclear star cluster. We show that a second generation forming out of gas polluted by the ejecta of massive AGB stars and mixed with gas having a standard composition accounts for the unusually large N/O in the GN-z11 spectrum. The timing of the evolution of massive (4-7.5M$_{\odot}$) AGBs also provides a favourable environment for the growth of a central stellar mass black hole to the Active Galactic Nucleus stage observed in GN-z11. According to our model the progenitor system was born at an age of the Universe of $\simeq 260 - 380$Myr, well within the pre-reionization epoch.

In this work, we present a study of the void lensing signal or the excess surface mass density (ESMD) around cosmic voids. First, we propose a new void-finder algorithm that is designed to capture the ESMD around voids. We compare our algorithm applied to projected slices with the ZOBOV void finder and find significantly deeper weak-lensing profiles for voids defined by our algorithm in the context of a realistic galaxy mock. Then we test the consistency between the measurements of the ESMD as measured through the shear of background galaxies and directly calculated through the dark matter density profiles of the same voids. We found inconsistencies for voids with diameter $\geq 100h^{-1}\mathrm{Mpc}$ along the line-of-sight, but the consistency holds for smaller voids, meaning that we are indeed probing the underlying dark matter field by measuring the shear around these voids. Moreover, we show that voids found in the projected slices, which are highly sensitive to lensing, are correlated to $3$D voids exhibiting intrinsic alignments between them.

In the beyond Standard Model (BSM) scenarios, the possibility of neutrinos decaying into a lighter state is still challenging to probe for the new-generation neutrino experiments. The observation of high-energy astrophysical neutrinos by IceCube opens up a new avenue for studying neutrino decay. In this work, we investigate a novel scenario of invisible neutrino decay to axion-like particles (ALPs). These ALPs propagate unattenuated and reconvert into gamma rays in the magnetic field of the Milky Way. This is complementary to the previously done studies where gamma rays produced at the source are used to investigate the ALP hypothesis. We exploit the Fermi-LAT and IceCube observations of NGC 1068 to set constraints on the ALP parameters. Being a steady source of neutrinos, it offers a better prospect over transient sources. We obtain 95% confidence level (C.L.) upper limits on the photon-ALP coupling constant $g_{a\gamma}\lesssim 1.37 \times 10^{-11}$ GeV$^{-1}$ for ALP masses $m_{a} \leq 2 \times 10^{-9}$ eV. Our results are comparable to previous upper limits obtained using the GeV to sub-PeV gamma-ray observations.

The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, particularly those due to planet-disk interactions. The multi-wavelength surface brightness, the disk color and geometry permit identification of any asymmetries such as warps or disk offsets from the central star. We find that nineteen of the disks in this sample exhibit asymmetries in surface brightness, disk color, disk geometry, or a combination of the three, suggesting that for this sample, perturbations, as seen in scattered light, are common. The relationship between these perturbations and potential planets in the system are discussed. We also explore correlations among stellar temperatures, ages, disk properties, and observed perturbations. We find significant trends between the vertical aspect ratio and the stellar temperature, disk radial extent, and the dust grain size distribution power-law, $q$. We also confirm a trend between the disk color and stellar effective temperature, where the disk becomes increasingly red/neutral with increasing temperature. Such results have important implications on the evolution of debris disk systems around stars of various spectral types.

Very recent work has identified a new satellite galaxy, Ursa Major III/UNIONS I, which is the faintest such system ever observed. Dynamical considerations indicate that if the system is in equilibrium, it is likely to be highly dark matter dominated. This, in combination with its proximity, predicts that it may be the preeminent dwarf spheroidal galaxy target for dark matter indirect detection searches. We utilize 15 years of Fermi-LAT data to search for $\gamma$-ray emission from Ursa Major III. Finding no excess, we set strong constraints on dark matter annihilation. Intriguingly, if the high J-factor of Ursa Major III is confirmed, standard thermal dark matter annihilation to $b\bar{b}$ final states would be ruled out for dark matter masses up to 4 TeV. The discovery of Ursa Major III, combined with recent tentative measurements of other high J-factor systems, suggests the exciting possibility that near-future data could produce transformative constraints on thermal dark matter.

Recently the Telescope Array collaboration reported an observation of cosmic ray event with very high energy 244 EeV ($2.44 \times 10^{20}$ eV). Importantly, the event is hard to correlate with the matter distribution in the local Universe, even after taking into account deflections in magnetic fields. This implies that the event is likely a nucleus with a large charge. An attenuation length of the nucleus of such a high energy in intergalactic space is quite small, therefore its source should be relatively close to our Galaxy. Using these arguments we derive a new upper bound on a distance to the closest ultra-high energy cosmic ray (UHECR) source and a lower bound on the general UHECR source number density. The distance to the closest source should not exceed 5 Mpc at 95% C.L. and the 95% C.L. lower-bound on the sources number density is $\rho > 1.0 \times 10^{-4}$ Mpc$^{-3}$. The number density of UHECR sources emitting heavy nuclei is constrained for the first time.

We extend the analysis of a physical model within the standard cosmology that robustly predicts a high star-formation efficiency (SFE) in massive galaxies at cosmic dawn due to feedback-free starbursts (FFB). It implies an excess of bright galaxies at z>~10 compared to the standard models based on the low SFE at later epochs, an excess that is indicated by JWST observations. Here we provide observable predictions based on the analytic FFB scenario, to be compared with simulations and JWST observations. We use the model to approximate the SFE as a function of redshift and mass, assuming a maximum SFE of 0.2~1 in the FFB regime. From this, we derive the evolution of the galaxy mass and luminosity functions as well as the evolution of stellar and star-formation densities. We then predict the star-formation history (SFH), galaxy sizes, outflows, gas fractions, metallicities and dust attenuation, all as functions of mass and redshift in the FFB regime. The major distinguishing feature is the occurrence of FFBs above a mass threshold that declines with redshift. The luminosities and star formation rates in bright galaxies are predicted to be in excess of extrapolations of standard empirical models and cosmological simulations, an excess that grows from z~9 to higher redshifts. The FFB phase of ~100Myr is predicted to show a characteristic SFH that fluctuates on a timescale of ~10Myr. The stellar systems are compact (Re~0.3kpc at z~10 and declining with z). The galactic gas consists of a steady wind driven by supernovae from earlier generations, with high outflow velocities (FWHM~1400-6700km/s), low gas fractions (<0.1), low metallicities (<~0.1solar), and low dust attenuation ($A_{UV}$~0.5 at z~10 and declining with z). Tentative comparisons with current JWST observations are made for preliminary impression, to be performed quantitatively when more complete reliable data become available.

We propose a new regime of minimal QCD axion dark matter that lies between the pre- and post-inflationary scenarios, such that the Peccei-Quinn (PQ) symmetry is restored only on sufficiently large spatial scales. This leads to a novel cosmological evolution, in which strings and domain walls re-enter the horizon and annihilate later than in the ordinary post-inflationary regime, possibly even after the QCD crossover. Such dynamics can occur if the PQ symmetry is restored by inflationary fluctuations, i.e. the Hubble parameter during inflation $H_I$ is larger than the PQ breaking scale $f_a$, but it is not thermally restored afterwards. Solving the Fokker-Planck equation, we estimate the number of inflationary e-folds required for the PQ symmetry to be, on average, restored. Moreover, we show that, in the large parts of parameter space where the radial mode is displaced from the minimum by de Sitter fluctuations, a string network forms due to the radial mode oscillating over the top of its potential after inflation. In both cases we identify order one ranges in $H_I/f_a$ and in the quartic coupling $\lambda$ of the PQ potential that lead to the late-string dynamics. In this regime the cosmological dark matter abundance can be reproduced for axion decay constants as low as the astrophysical constraint $O(10^8)$ GeV, corresponding to axion masses up to $10^{-2}~{\rm eV}$, and with miniclusters with masses as large as $O(10)M_\odot$.

Scale-invariant extensions of the electroweak theory are not only attractive because they can dynamically generate the weak scale, but also due to their role in facilitating supercooled cosmological first-order phase transitions. We study the minimal scale-invariant $U(1)_{\rm D}$ extension of the standard model and find that Primordial Black Holes (PBHs) can be abundantly produced. The PBH mass is bounded from above by the moon mass due to QCD catalysis limiting the amount of supercooling. Lunar-mass PBHs, which are produced for dark Higgs vev $v_\phi\simeq 20~\rm TeV$, correspond to the best likelihood to explain the HSC lensing anomaly. For $v_\phi\gtrsim 400~\rm TeV$, the model can explain hundred per cent of dark matter. At even larger hierarchy of scales, it can contribute to the $511~\rm keV$ line. While the gravitational wave (GW) signal produced by the HSC anomaly interpretation is large and detectable by LISA above astrophysical foreground, the dark matter interpretation in terms of PBHs can not be entirely excluded by GW detection due to the dilution of the signal by the entropy injected during the decay of the long-lived $U(1)_{\rm D}$ scalar. This long lifetime is a natural consequence of the large hierarchy of scales.

Tensions in measurements of neutron and kaon weak decays, such as of the neutron lifetime, may speak to the existence of new particles and dynamics not present in the Standard Model (SM). In scenarios with dark sectors, particles that couple feebly to those of the SM appear. We offer a focused overview of such possibilities and describe how observations of neutron stars, that probe either their structure or dynamics, limit them. In realizing these constraints we highlight how the assessment of particle processes within dense baryonic matter impacts the emerging picture -- and we emphasize both the flavor structure of the constraints and their broader connections to cogenesis models of dark matter and baryogenesis.

Compact objects such as neutron stars and white dwarfs can source axion-like particles and QCD axions due to CP-violating axion-fermion couplings. The magnitude of the axion field depends on the stellar density and on the strength of the axion-fermion couplings. We show that even CP-violating couplings one order of magnitude smaller than existing constraints source extended axion field configurations. For axion-like particles, the axion energy is comparable to the magnetic energy in neutron stars with inferred magnetic fields of the order of $10^{13}$ G, and exceeds by more than one order of magnitude the magnetic energy content of white dwarfs with inferred fields of the order of $10^{4}$ G. On the other hand, the energy stored in the QCD axion field is orders of magnitude lower due to the smallness of the predicted CP-violating couplings. It is shown that the sourced axion field can polarize the photons emitted from the stellar surface, and stimulate the production of photons with energies in the radio band.

The atmospheric convective boundary layer (CBL) consists in three basic parts: (i) the surface layer unstably stratified and dominated by small-scale turbulence of very complex nature; (ii) the CBL core dominated by the energy-, momentum- and mass-transport of semi-organised structures (large-scale circulations), with a small contribution from small-scale turbulence produced by local structural shears; (iii) turbulent entrainment layer at the upper boundary, characterised by essentially stable stratification with negative (downward) turbulent flux of potential temperature. The energy- and flux budget (EFB) theory developed previously for atmospheric stably-stratified turbulence and the surface layer in atmospheric convective turbulence is extended to the CBL core using budget equations for turbulent energies and turbulent fluxes of buoyancy and momentum. For the CBL core, we determine global turbulent characteristics (averaged over the entire volume of the semi-organised structure) as well as kinetic and thermal energies of the semi-organised structures as the functions of the aspect ratio of the semi-organised structure and the scale separation parameter between the vertical size of the structures and the integral scale of turbulence. The obtained theoretical relationships are potentially useful in modelling applications in the atmospheric convective boundary-layer.

In this paper, we explore the behavior of light ray trajectories in the exterior geometry of a rotating black hole within the Rastall theory of gravity, which is surrounded by an inhomogeneous anisotropic electronic cold plasma. By specifying the plasma's frequency profile, we derive fully analytical solutions for the temporal evolution of spacetime coordinates using elliptic integrals and Jacobi elliptic functions. These solutions illustrate various possible orbits. Throughout the study, we compare results with the vacuum case, emphasizing the influence of the plasma. Additionally, we utilize the analytical solutions to establish the lens equation for the considered spacetime. The investigation also addresses the significance of spherical photon orbits in critical trajectories, by presenting several examples.

Following our previous work, we continue to explore gravitational dark matter production during the minimal preheating caused by inflaton self-resonance. In this situation there is only one dimensionless index parameter $n$ characterizing the inflation potential after the end of inflation, which leads to a robust prediction on the gravitational dark matter relic abundance. Using lattice method to handle the non-perturbative evolutions of relevant quantities during the inflaton self-resonance, we derive the gravitational dark matter relic abundance arising from both the inflaton condensate and fluctuation annihilation. While being absent for $n=2$, the former one can instead dominate over the later one for $n=4,6$. Our results show that gravitational dark matter mass of $1.04~(2.66)\times 10^{14}$ GeV accommodates the observed value of dark matter relic abundance for $n=4$ (6).

There are life forms in space and the ancestor of Earth life came from interstellar space according to models like panspermia. Naturally, life may also exist in molecular clouds. Here the author discusses the possibility of methanogenic life in Molecular Cloud with methane as the final metabolism product. According to the calculations, it is easy to see that the chemical reaction through methanogenesis can release sufficient free energy. If methanogenic life exist in the pre-solar nebula, then they may be the ancestor of Earth's life and there are already some tentative evidences by several molecular biology studies.

Coupling of axions or axion-like particles (ALPs) with photons may lead to photons escaping optically opaque regions by oscillating into ALPs. This phenomenon may be viewed as the Light Shining through Wall (LSW) scenario. While this LSW technique has been used previously in controlled laboratory settings to constrain the ALP-photon coupling ($g_{a\gamma}$), we show that this can also be applied in astrophysical environments. We find that obscured magnetars in particular are excellent candidates for this purpose. A fraction of photons emitted by the magnetar may convert to ALPs in the magnetar neighborhood, cross the large absorption column densities, and convert back into photons due to the interstellar magnetic field. Comparing the observed flux with the estimated intrinsic flux from the magnetar, we can constrain the contribution of this process, and hence constrain $g_{a\gamma}$. The effects of resonant conversion near the magnetar as well as ALP-photon oscillations in the interstellar medium are carefully considered. Taking a suitable magnetar candidate PSR J1622-4950, we find that the ALP-photon coupling can be constrained at $g_{a\gamma} \lesssim (10^{-10} - 10^{-11})$ GeV$^{-1}$ for low mass axions ($m_a \lesssim 10^{-12}$ eV). Our study reveals the previously unrealized potential for employing the LSW technique for obscured magnetars for probing and constraining ALP-photon couplings.

We discuss the standard Lane-Emden formalism as well as the one related to the slowly rotating objects. It is preceded by a brief introduction of different forms of the polytropic equation of state. This allows to study a wide class of astrophysical objects in the framework of a given theory of gravity, as demonstrated in a few examples. We will discuss light elements burning processes and cooling models in stars and substellar objects with the use of the Lane-Emden formalism.

Context. The mechanisms regulating the transport and energization of charged particles in space and astrophysical plasmas are still debated. Plasma turbulence is known to be a powerful particle accelerator. Large-scale structures, including flux ropes and plasmoids, may contribute to confine particles and lead to fast particle energization. These structures may also modify the properties of the turbulent nonlinear transfer across scales. Aims. We investigate how large-scale flux ropes are perturbed and, simultaneously, influence the nonlinear transfer of turbulent energy towards smaller scales. We then address how these structures affect particle transport and energization. Methods. We adopt magnetohydrodynamic simulations for perturbing a large-scale flux rope in solar-wind conditions and possibly triggering turbulence. Then, we employ test-particle methods to investigate particle transport and energization in the perturbed flux rope. Results. The large-scale helical flux rope inhibits the turbulent cascade towards smaller scales, especially if the amplitude of the initial perturbations is not large (~5%). In this case, particle transport is inhibited inside the structure. Fast particle acceleration occurs in association with phases of trapped motion within the large-scale flux rope.

We study neutrino spin oscillations while the particles scatter off a rotating black hole surrounded by a thick magnetized accretion disk. Neutrino spin precession is caused by the interaction of the neutrino magnetic moment with the magnetic field in the disk which has both toroidal and poloidal components. Our calculation of the observed neutrino fluxes, accounting for spin oscillations, are based on numerical simulations of the propagation of a great number of incoming test particles using High Performance Parallel Computing. The obtained results significantly improve our previous findings. We briefly discuss the applications for the observations of astrophysical neutrinos.

Dark matter particles could be superheavy, provided their lifetime is much longer than the age of the universe. Using the sensitivity of the Pierre Auger Observatory to ultra-high energy neutrinos and photons, we constrain a specific extension of the Standard Model of particle physics that meets the lifetime requirement for a superheavy particle by coupling it to a sector of ultra-light sterile neutrinos. Our results show that, for a typical dark coupling constant of 0.1, the mixing angle $\theta_m$ between active and sterile neutrinos must satisfy, roughly, $\theta_m \lesssim 2.5\times 10^{-6}(M_X/10^9~\mathrm{GeV})^{-2}$ for a mass $M_X$ of the dark-matter particle between $2.3\times 10^8$ and $10^{11}~$GeV.

We present a summary of the main results within the Scale Invariant Vacuum (SIV) paradigm based on the Weyl Integrable Geometry (WIG) as an extension to the standard Einstein General Relativity (EGR). After a brief review of the mathematical framework, where we also highlight the connection between the weak-field SIV equations and the notion of un-proper time parametrization within the reparametrization paradigm, we continue with the main results related to early Universe; that is, applications to inflation, Big Bang Nucleosynthesis, and the growth of the density fluctuations within the SIV. In the late time Universe the applications of the SIV paradigm are related to scale-invariant dynamics of galaxies, MOND, dark matter, and the dwarf spheroidals where one can find MOND to be a peculiar case of the SIV theory. Finally, within the recent time epoch, we highlight that some of the change in the length-of-the-day (LOD), about 0.92 cm/yr, can be accounted for by SIV effects in the Earth-Moon system.

A recent study has shown that domain walls with initial inflationary fluctuations exhibit remarkable stability against population bias due to long-range correlations, challenging the claims of prior research. In this paper, we study the dynamics of these domain walls in the presence of potential bias and show that they collapse with a lifetime several times longer than that due to thermal fluctuations. This is interpreted as a difference in the average distance between domain walls, leading us to derive a new formula for the domain wall lifetime that depends on the area parameter. In addition, we compute the spectrum of gravitational waves generated by such domain walls and find that both the peak frequency and the peak abundance are lowered in a manner that depends on the area parameter. Based on these findings, we also determine the necessary degree of vacuum degeneracy for axion domain walls to explain the isotropic cosmic birefringence.

The study of volatiles and the search for water are the primary objectives of the Luna-27 mission, which is planned to land on the south pole of the Moon in 2028. Here we present the tunable Diode Laser Spectrometer (DLS-L) that will be onboard the lander. The DLS-L will perform isotopic analysis of volatiles that are pyrolytically evolved from regolith. This article dives into the design of the spectrometer and the characterisation of isotopic signature retrieval. We look forward to expanding our knowledge of Lunar geochemistry by measuring D/H, 18O/17O/16O, 13C/12C ratios in situ, which could be the one-of-a-kind direct study of the lunar soil isotopy without sample contamination.

Many areas of science and engineering encounter data defined on spherical manifolds. Modelling and analysis of spherical data often necessitates spherical harmonic transforms, at high degrees, and increasingly requires efficient computation of gradients for machine learning or other differentiable programming tasks. We develop novel algorithmic structures for accelerated and differentiable computation of generalised Fourier transforms on the sphere $\mathbb{S}^2$ and rotation group $\text{SO}(3)$, i.e. spherical harmonic and Wigner transforms, respectively. We present a recursive algorithm for the calculation of Wigner $d$-functions that is both stable to high harmonic degrees and extremely parallelisable. By tightly coupling this with separable spherical transforms, we obtain algorithms that exhibit an extremely parallelisable structure that is well-suited for the high throughput computing of modern hardware accelerators (e.g. GPUs). We also develop a hybrid automatic and manual differentiation approach so that gradients can be computed efficiently. Our algorithms are implemented within the JAX differentiable programming framework in the S2FFT software code. Numerous samplings of the sphere are supported, including equiangular and HEALPix sampling. Computational errors are at the order of machine precision for spherical samplings that admit a sampling theorem. When benchmarked against alternative C codes we observe up to a 400-fold acceleration. Furthermore, when distributing over multiple GPUs we achieve very close to optimal linear scaling with increasing number of GPUs due to the highly parallelised and balanced nature of our algorithms. Provided access to sufficiently many GPUs our transforms thus exhibit an unprecedented effective linear time complexity.