Papers for Thursday, Feb 06 2025

Using a dataset of stars from a barred galaxy that provides instantaneous positions and velocities in the sky, we propose a method to identify its key structures, such as the central bar and spiral arms, and to determine the bar pattern speed in the short term. This is accomplished by applying combinatorial and topological data analysis techniques to the available star-state information set. To validate the methodology, we apply it to two scenarios: a test particle simulation and a N-body simulation, both evaluated when the bar is fully formed. In both cases, the results are robust and consistent. We then use these outcomes to calibrate a classical dynamical model and compute related structures, such as equilibrium points and invariant manifolds, which we verify align with the galaxy's spiral arms.

The field of the Search for Extraterrestrial Intelligence (SETI) searches for ``technosignatures'' that could provide the first detection of life beyond Earth through the technology that an extraterrestrial intelligence (ETI) may have created. Any given SETI survey, if no technosignatures are detected, should set upper limits based on the kinds of technosignatures it should have been able to detect; the sensitivity of many SETI searches requires that their target sources (e.g., Dyson spheres or Kardashev II/III level radio transmitters) emit with power far exceeding the kinds of technology humans have developed. In this paper, we instead turn our gaze Earthward, minimizing the axis of extrapolation by only considering transmission and detection methods commensurate with an Earth-2024 level. We evaluate the maximum distance of detectability for various present-day Earth technosignatures -- radio transmissions, atmospheric technosignatures, optical and infrared signatures, and objects in space or on planetary surfaces -- using only present-day Earth instruments, providing one of the first fully cross-wavelength comparisons of the growing toolbox of SETI techniques. In this framework, we find that Earth's space-detectable signatures span 13 orders of magnitude in detectability, with intermittent, celestially-targeted radio transmission (i.e., planetary radar) beating out its nearest non-radio competitor by a factor of $10^3$ in detection distance. This work highlights the growing range of ways that exoplanet technosignatures may be expressed, the growing complexity and visibility of the human impact upon our planet, and the continued importance of the radio frequencies in SETI.

Cosmological simulations provide much of the theoretical framework within which we interpret extragalactic observations. However, even if a given simulation reproduces the integrated properties of galaxies well, it may not reproduce the detailed structures of individual galaxies. Comparisons between the 2D light distributions of simulated and observed galaxies -- particularly in the dwarf regime, where key processes like tidal perturbations and baryonic feedback most strongly influence galaxy structure -- thus provide an additional valuable test of the simulation's efficacy. We compare scaling relations derived from mock observations of simulated galaxies, drawn from the two largest halos in the high-resolution NewHorizon cosmological simulation, with galaxies in the Fornax cluster. While Fornax is significantly more massive than either group, it is the lowest-mass cluster in the local Universe, and contains a well-studied population of spatially resolved dwarfs, hence serves as a useful benchmark. Per unit stellar mass, NewHorizon dwarfs are systematically larger in half-light radius, much fainter in surface brightness, and bluer in colour than their Fornax counterparts, albeit with similar light profile shapes. We discuss potential reasons for these discrepancies, including environmental effects, baryonic feedback, resolution, or couplings of these factors. As observations of dwarfs outside of the local Universe become more plentiful through on-going or up-coming surveys such as Euclid and LSST, 2D comparisons such as these, where properties are measured in the same way across both simulations and observations, can place strong constraints on processes that alter the spatial distribution of baryons in galaxies.

This study explores the dynamical impact of cosmic rays (CRs) in Milky Way-like galaxies using the Rhea simulation suite. Cosmic rays, with their substantial energy density, influence the interstellar medium (ISM) by supporting galactic winds, modulating star formation, and shaping ISM energetics. The simulations incorporate a multi-phase ISM, self-consistent CR transport in the advection-diffusion approximation, and interactions with magnetic fields to study their effect on galaxy evolution. Key findings reveal that CRs reduce star formation rates, and drive weak but sustained outflows with mass loading factors of $\sim0.2$, transporting a substantial fraction (20%-60%) of the injected CR energy. These CR-driven outflows are launched not just from the galactic center but across the entire disk, illustrating their pervasive dynamical influence. Galactic disks supported by CRs exhibit broader vertical structures compared to magnetic-field-dominated setups, though the scale heights are similar. CR feedback enhances magnetic flux transport to the circumgalactic medium (CGM), leading to a magnetically enriched CGM with field strengths of $\sim0.5\mu\mathrm{G}$ while reducing gas temperatures to $\lesssim10^5\,\mathrm{K}$. The CR energy is relatively smoothly distributed in the disk, with gradient lengths exceeding the typical size of molecular clouds, indicating that the CR behavior is not adiabatic.

We explore interacting dark matter (DM) models that allow DM and baryons to scatter off of each other with a cross section that scales with relative particle velocity. Using the effective field theory of large-scale structure, we perform the first analysis of BOSS full-shape galaxy clustering data for velocity-dependent DM-baryon interactions. We determine that while the addition of BOSS full-shape data visibly modifies the shape of the posterior distribution, it does not significantly alter the 95% confidence level intervals for the interaction cross section obtained from an analysis of the cosmic microwave (CMB) anisotropy from Planck measurements alone. Moreover, in agreement with previous findings, we note that the DM-baryon interacting model presents a good fit to both large-scale structure (LSS) data and CMB data and alleviates the $S_8$ tension between the two data sets. After combining LSS and CMB data with weak lensing data from the Dark Energy Survey, we find a $\gtrsim2\sigma$ preference for non-zero interactions between DM and baryons in a velocity-independent model. We also explore a scenario where only a fraction of DM undergoes scattering with baryons; we find a similar $\gtrsim2\sigma$ preference for the presence of interactions. Our results suggest that a suppression of the linear matter power spectrum at small scales may be needed to resolve certain discrepancies between LSS and CMB data that are found in the cold DM (CDM) scenario.

The James Webb Space Telescope (JWST) is revolutionizing our understanding of the Universe by unveiling faint, near-infrared dropouts previously beyond our reach, ranging from exceptionally dusty sources to galaxies up to redshift $z \sim 14$. In this paper, we identify F200W-dropout objects in the Cosmic Evolution Early Release Science (CEERS) survey which are absent from existing catalogs. Our selection method can effectively identify obscured low-mass ($\log \text{M}_* \leq 9$) objects at $z \leq 6$, massive dust-rich sources up to $z \sim 12$, and ultra-high-redshift ($z > 15$) candidates. Primarily relying on NIRCam photometry from the latest CEERS data release and supplementing with Mid-Infrared/(sub-)mm data when available, our analysis pipeline combines multiple SED-fitting codes, star formation histories, and CosMix - a novel tool for astronomical stacking. Our work highlights three $2<z<3$ dusty dwarf galaxies which have larger masses compared to the typical dusty dwarfs previously identified in CEERS. Additionally, we reveal five faint sources with significant probability of lying above $z>15$, with best-fit masses compatible with $\Lambda$CDM and a standard baryons-to-star conversion efficiency. Their bi-modal redshift probability distributions suggest they could also be $z<1.5$ dwarf galaxies with extreme dust extinction. We also identify a strong line emitter galaxy at $z \sim 5$ mimicking the near-infrared emission of a $z \sim 13$ galaxy. Our sample holds promising candidates for future follow-ups. Confirming ultra high-redshift galaxies or lower-z dusty dwarfs will offer valuable insights into early galaxy formation, evolution with their central black holes and the nature of dark matter, and/or cosmic dust production mechanisms in low-mass galaxies, and will help us to understand degeneracies and contamination in high-z object searches.

In this paper, we propose an Alcock-Paczyński (AP) test to constrain cosmology using HII bubbles during the Epoch of Reionization. Similarly to cosmic voids, a stack of HII bubbles is spherically symmetric because ionizing fronts propagate isotropically on average (even if individual bubbles may not be spherical), making them standard spheres to be used in an AP test. Upcoming 21-cm observations, from the Square Kilometer Array (SKA) for instance, will contain tomographic information about HII regions during reionization. However, extracting the bubbles from this signal is made difficult because of instrumental noise and foreground systematics. Here, we use a neural network to reconstruct neutral-fraction boxes from the noisy 21-cm signal, from which we extract bubbles using a watershed algorithm. We then run the purely geometrical AP test on these stacks, showing that a SKA-like experiment will be able to constrain the product of the angular-diameter distance $D_{\text{A}}$ and Hubble parameter $H$ at reionization redshifts with $\sim 2\%$ precision, robustly to astrophysical and cosmological uncertainties within the models tested here. This AP test, whether performed on 21-cm observations or other large surveys of ionized bubbles, will allow us to fill the knowledge gap about the expansion rate of our Universe at reionization redshifts.

Rapidly rotating classical OBe stars have been proposed as the products of binary interactions, and the fraction of Be stars with compact companions implies that at least some are. However, to constrain the interaction physics spinning up the OBe stars, a large sample of homogeneously analysed OBe stars with well-determined binary characteristics and orbital parameters are required. We investigate the multiplicity properties of a sample of 18 Oe, 62 Be, and two Of?p stars observed within the BLOeM survey in the Small Magellanic Cloud. We analyse the first nine epochs of spectroscopic observations obtained over approximately three months in 2023. Radial velocities (RVs) of all stars are measured. Applying commonly-used binarity criteria we classify objects as binaries, binary candidates, and apparently single (RV stable) objects. We further inspect the spectra for double-lined spectroscopic binaries and cross-match with catalogues of X-ray sources and photometric binaries. We classify 14 OBe stars as binaries, and an additional 11 as binary candidates. The two Of?p stars are apparently single. Two more objects are most likely currently interacting binaries. Without those, the observed binary fraction for the OBe sample (78 stars) is f_OBe_obs=0.18+/-0.04 (f_obs_cand=0.32+/-0.05 including candidates). This fraction is less than half of that measured for OB stars in BLOeM. Combined with the lower fraction of SB2s, this suggests that OBe stars have indeed fundamentally different binary properties than OB stars. We find no evidence for OBe binaries with massive compact companions, in contrast to expectations from binary population synthesis. Our results support the binary scenario as an important formation channel for OBe stars, as post-interaction binaries may have been disrupted or the stripped companions of OBe stars are harder to detect.

Radial migration is an important dynamical effect that has reshaped the Galactic disc, but its origin has yet to be elucidated. In this work, we present evidence that resonant dragging by the corotation of a decelerating bar could be the main driver of radial migration in the Milky Way disc. Using a test particle simulation, we demonstrate this scenario explains the two distinct age-metallicity sequences observed in the solar vicinity: the plateauing upper sequence is interpreted as stars dragged outwards by the expanding corotation of the decelerating bar and the steeper lower sequence as stars formed locally around the solar circle. The upper migrated sequence dominates at guiding radii around the current corotation radius of the bar, $R\sim7\,\mathrm{kpc}$, but rapidly dies away beyond this where the mechanism cannot operate. This behaviour naturally explains the radial dependence of the $\mathrm{[\alpha/Fe]}$-bimodality, in particular the truncation of the high-$\mathrm{[\alpha/Fe]}$ disc beyond the solar circle. Under our proposed radial migration scenario, we constrain the Milky Way bar's pattern speed evolution using the age-metallicity distribution of stars currently trapped at corotation. We find the bar likely formed with an initial pattern speed of $60-100$ km s$^{-1}$ kpc$^{-1}$ and began decelerating $6-8$ Gyr ago at a rate $-\dot{\Omega}/\Omega^2\sim0.0025-0.0040$ (where the quoted ranges include systematic uncertainties).

Given the uncertain evolutionary status of blue supergiant stars, their multiplicity properties hold vital clues to better understand their origin and evolution. As part of The Binarity at LOw Metallicity (BLOeM) campaign in the Small Magellanic Cloud we present a multi-epoch spectroscopic survey of 128 supergiant stars of spectral type B5--F5, which roughly correspond to initial masses in the range 6 to 30 solar masses. The observed binary fraction for the B5-9 supergiants is 25+/-6 % (10+/-4 %) and 5+/-2 % (0 %) for the A-F stars, using a radial velocity (RV) variability threshold of 5 kms (10 kms) as a criterion for binarity. Accounting for observational biases we find an intrinsic multiplicity fraction of less than 18% for the B5-9 stars and 8$^{+9}_{-7}$% for the AF stars, for the orbital periods up to 10$^{3.5}$day and mass-ratios (q) in the range 0.1 < q < 1. The large stellar radii of these supergiant stars prevent short orbital periods but we demonstrate that this effect alone cannot explain our results. We assess the spectra and RV time series of the detected binary systems and find that only a small fraction display convincing solutions. We conclude that the multiplicity fractions are compromised by intrinsic stellar variability such that the true multiplicity fraction may be significantly smaller. Our main conclusions from comparing the multiplicity properties of the B5-9 and AF supergiants to that of their less evolved counterparts is that such stars cannot be explained by a direct evolution from the main sequence. Furthermore, by comparing their multiplicity properties to red supergiant stars we conclude that the AF supergiant stars are neither progenitors nor descendants of red supergiants.

The Milky Way is a complex ecosystem, for which we can obtain detailed observations probing the physical mechanisms determining the interstellar medium. For a detailed comparison with observations, and to provide theories for missing observables, we need to model the Milky Way as closely as possible. However, details of the Galactic structure are not fully defined by observations, raising the need for more generalized models. With the Rhea simulations we present a set of Milky Way like simulations, containing detailed physics of the interstellar medium, as well as star formation and stellar feedback. We conduct two simulations that differ in the gravitational potential: one fitted to several structural details derived from observations, the other just reproducing the most basic quantities. We find little difference in the overall morphology except for the bar region, which funnels gas towards the Galactic inner region and therefore prevents quenching in the center. Despite differences with galacto-centric radius, the global star formation rate is almost identical in both setups. A spiral arm potential does not influence properties of groups of formed stars. A bar potential, however, lowers size and formation time of those groups. We therefore conclude for a spiral arm potential to have little influence on star formation in the Galaxy, except for producing long-lived spiral structures instead of transient ones. A Galactic bar potential has noticeable influence on star formation mainly within the innermost 2.5kpc.

Central bars and spirals are known to strongly impact the evolution of their host galaxies, both in terms of dynamics and star formation. Their typically different pattern speeds cause them to regularly overlap, which induces fluctuations in bar parameters. In this paper, we analyze both numerical simulations of disk galaxies and observational data to study the effect of bar-spiral physical overlap on stellar radial migration and star formation in the bar vicinity, as a function of time and galactic azimuth. We study three different numerical models, two of which are in a cosmological context, as well as APOGEE DR17 data and the WISE catalog of Galactic HII regions. We find that periodic boosts in stellar radial migration occur when the bar and spiral structure overlap. This mechanism causes net inward migration along the bar leading side, while stars along the bar trailling side and minor axis are shifted outward. The signature of bar-spiral induced migration is seen between the bar's inner Lindbald resonance and well outside its corotation, beyond which other drivers take over. We also find that, in agreement with simulations, APOGEE DR17 stars born at the bar vicinity (mostly metal-rich) can migrate out to the solar radius while remaining on cold orbits. For the Milky Way, 13% of stars in the solar vicinity were born inside the bar, compared to 5-20% in the simulations. Bar-spiral reconnections also result in periodic starbursts at the bar ends with an enhancement of up to a factor of 4, depending on the strength of the spiral structure. Similarly to the migration bursts, these do not always happen simultaneously at the two sides of the bar, hinting at the importance of odd spiral modes. Data from the WISE catalog suggest this phhenomenon is also relevant in our own Galaxy.

Dwarf galaxies dominate the galaxy number density, making them critical to our understanding of galaxy evolution. However, typical dwarfs are too faint to be visible outside the very local Universe in past surveys like the SDSS, which offer large footprints but are shallow. Dwarfs in such surveys have relatively high star formation rates, which boost their luminosity, making them detectable in shallow surveys, but also biased and potentially unrepresentative of dwarfs as a whole. Here, we use deep data to perform an unbiased statistical study of ~7,000 nearby (z<0.25) dwarfs (10^8 MSun < M < 10^9.5 MSun) in the COSMOS field which, at these redshifts, is a relatively low-density field. At z~0.05, ~40 per cent of dwarfs in low-density environments are red/quenched, falling to ~30 per cent by z~0.25. Red dwarfs reside closer to nodes, filaments and massive galaxies. Proximity to a massive galaxy appears to be more important in determining whether a dwarf is red, rather than simply its distance from nodes and filaments or the mean density of its local environment. Interestingly, around half of the red dwarfs reside outside the virial radii of massive galaxies and around a third of those also inhabit regions in the lower 50 per cent in density percentile (i.e. regions of very low ambient density). Around half of the red dwarf population is, therefore, quenched by mechanisms unrelated to environment, which are likely to be internal processes such as stellar and AGN feedback.

Runaway stars with peculiar high velocities can generate stellar bow shocks. Only a few bow shocks show clear radio emission. Our goal is to identify and characterize new stellar bow shocks around O and Be runaway stars in the infrared (IR), and to study their possible radio emission and nature. Our input data is a catalog of O and Be runaways compiled using Gaia DR3. We used WISE IR images to search for bow shocks around these runaways, Gaia DR3 data to determine the actual motion of the runaway stars corrected for interstellar medium (ISM) motion caused by Galactic rotation, and archival radio data to search for emission signatures. We finally explored the radio detectability of these sources under thermal and nonthermal scenarios. We found 9 new stellar bow shock candidates, 3 new bubble candidates, and 1 intermediate structure candidate. One of them is an in situ bow shock candidate. We also found 17 already known bow shocks in our sample, though we discarded one, and 62 miscellaneous sources showing some IR emission around the runaways. We geometrically characterized the sources in IR using the WISE-4 band and estimated the ISM density at the bow shock positions, obtaining median values of ~6 and ~4 cm$^{-3}$ using 2D and 3D peculiar velocities. Most of the new discovered bow shocks come from new runaway discoveries. Within our samples we found that ~24% of the O-type runaway stars show bow shocks, while this decreases to ~3% for Be-type runaway stars. Two bow shocks present radio emission but not as clear counterparts, and two others show hints of radio emission. The physical scenarios indicate that two sources could still be compatible with nonthermal radio emission. The new sample of O and Be runaway stars allowed us to discover both new stellar bow shocks and bubbles. Their geometrical characterization can be used to assess the physical scenario of the radio emission. (Abridged)

Many astrophysical explosions, such as type Ia supernovae, classical novae, and X-ray bursts, are dominated by thermonuclear runaway. To model these processes accurately, one must evolve nuclear reactions concurrently with hydrodynamics. We present an application of the moving mesh technique to this field of computation with the aim of explicitly testing the advantages of the method against the fixed mesh case. By way of traditional Strang splitting, our work couples a 13 isotope nuclear reaction network to a 1D moving mesh, Cartesian geometry hydrodynamics code. We explore three reacting problems: an acoustic pulse, a burning shock, and an advecting deflagration. Additionally using the shock jump conditions, we semi-analytically solve the burning shock problem under the assumption of quick, complete burning with the hope of establishing a useful and easy to set-up test problem. Strong moving mesh advantages are found in advecting, deflagrating flame fronts, where the technique dramatically reduces numerical diffusion that would otherwise lead to very fast artificial deflagration.

Foundational models have emerged as a powerful paradigm in deep learning field, leveraging their capacity to learn robust representations from large-scale datasets and effectively to diverse downstream applications such as classification. In this paper, we present Astromer 2 a foundational model specifically designed for extracting light curve embeddings. We introduce Astromer 2 as an enhanced iteration of our self-supervised model for light curve analysis. This paper highlights the advantages of its pre-trained embeddings, compares its performance with that of its predecessor, Astromer 1, and provides a detailed empirical analysis of its capabilities, offering deeper insights into the model's representations. Astromer 2 is pretrained on 1.5 million single-band light curves from the MACHO survey using a self-supervised learning task that predicts randomly masked observations within sequences. Fine-tuning on a smaller labeled dataset allows us to assess its performance in classification tasks. The quality of the embeddings is measured by the F1 score of an MLP classifier trained on Astromer-generated embeddings. Our results demonstrate that Astromer 2 significantly outperforms Astromer 1 across all evaluated scenarios, including limited datasets of 20, 100, and 500 samples per class. The use of weighted per-sample embeddings, which integrate intermediate representations from Astromer's attention blocks, is particularly impactful. Notably, Astromer 2 achieves a 15% improvement in F1 score on the ATLAS dataset compared to prior models, showcasing robust generalization to new datasets. This enhanced performance, especially with minimal labeled data, underscores the potential of Astromer 2 for more efficient and scalable light curve analysis.

Carbon-rich Wolf-Rayet binaries are a prominent source of carbonaceous dust that contribute to the dust budget of galaxies. The "textbook" example of an episodic dust producing WR binary, WR140 (HD193793), provides us with an ideal laboratory for investigating the dust physics and kinematics in an extreme environment. This study is among the first to utilize two separate JWST observations, from Cycle 1 ERS (July 2022) and Cycle 2 (Sept. 2023), to measure WR140's dust kinematics and confirm its morphology. To measure the proper motions and projected velocities of the dust shells, we performed a novel PSF subtraction to reduce the effects of the bright diffraction spikes and carefully aligned the Cycle 2 to the Cycle 1 images. At 7.7 $\mu$m, through the bright feature common to 16 dust shells (C1), we find an average dust shell proper motion of $390\pm29$ mas yr$^{-1}$, which equates to a projected velocity of $2714\pm188$ km s$^{-1}$ at a distance of 1.64 kpc. Our measured speeds are constant across all visible shells and consistent with previously reported dust expansion velocities. Our observations not only prove that these dusty shells are astrophysical (i.e., not associated with any PSF artifact) and originate from WR140, but also confirm the "clumpy" morphology of the dust shells, in which identifiable substructures within certain shells persist for at least 14 months from one cycle to the next. These results support the hypothesis that clumping in the wind collision region is required for dust production in WR binaries.

Three uniquely powerful solar and heliospheric facilities are now operational at the same time. The US National Science Foundation's Daniel K Inouye Solar Telescope, NASA's Parker Solar Probe, and ESA's Solar Orbiter each represent frontiers in space science, and each pursue richly tailored science missions. At the intersection of these missions, though, lie unparalleled opportunities for multi-vantage point science. This symbiotic relationship is especially pronounced during PSP's perihelia and Solar Orbiter remote science windows. As the most advanced solar polarimeter ever built, DKIST strengthens many of the multi-facility use cases by opening new diagnostic windows into solar magnetism -- spanning the photosphere, chromosphere, and corona -- at unprecedented spatial, spectral, and temporal resolution. In this article, we report recent efforts to maximize the scientific potential of coordinated DKIST, PSP, and Solar Orbiter observations. Existing DKIST data from coordinated observations with Solar Orbiter and PSP are highlighted alongside some first investigations of these data.

The halo model offers a framework for investigating galaxy clustering, and for understanding the growth of galaxies and the distribution of galaxies of different types. Here, we use the halo model to study the small-scale clustering and halo occupation distribution (HOD) of the unWISE-Blue galaxy sample, an infrared-selected sample of $\sim$100 million galaxies across the entire extragalactic sky at $z\sim 0.5$ $-$ similar redshifts to the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample. Although the photometric unWISE galaxies cannot be easily split in redshift, we use their cross-correlation with the BOSS CMASS sample to tomographically probe the HOD of the unWISE galaxies at $0.45 < z < 0.75$. To do so, we develop a new method for applying the halo model to cross-correlations between a photometric sample and a spectroscopic sample in narrow redshift bins, incorporating halo exclusion, post-Limber corrections, and redshift-space distortions. We reveal strong evolution in the CMASS HOD, and modest evolution in the unWISE-Blue HOD. For unWISE-Blue, we find that the average bias and mean halo mass drop from $b = 1.6$ and $\log_{10}(M_{\mathrm{h}}/M_{\odot}) \sim 13.4$ at $z \sim 0.5$ to $b = 1.4$ and $\log_{10}(M_{\mathrm{h}}/M_{\odot}) \sim 13.1$ at $z \sim 0.7$, and that the satellite fraction drops modestly from $\sim$20% to $\sim$10% in the same range. These results are useful for creating mock samples of the unWISE-Blue galaxies. Furthermore, the techniques developed to obtain these results are applicable to other tomographic cross-correlations between photometric samples and narrowly-binned spectroscopic samples, such as clustering redshifts.

We present a high-power continuous-wave (CW) lunar laser ranging (LLR) technique that has the potential to significantly improve Earth--Moon distance measurements. Using a 1 kW CW laser at 1064 nm and a 1 m-aperture telescope as an example, we develop a detailed link budget and analyze the prevailing noise sources to assess system performance when ranging to next-generation ~10 cm corner-cube retroreflectors (CCRs). Unlike legacy arrays, these smaller CCRs are designed to yield lower intrinsic range errors, yet their reduced reflective area results in lower photon return rates, posing challenges for pulsed LLR systems. The photon-rich CW approach, by providing continuous high-power illumination, overcomes this limitation, reducing shot noise and enabling sustained millimeter-level ranging with a pathway to sub-0.1 mm precision. Furthermore, by alternating measurements between widely separated lunar reflectors, differential LLR mitigates common-mode station errors to achieve tens-of-micrometer precision, limited primarily by uncorrelated atmospheric turbulence. This scalable approach -- integrating high-power CW lasers, narrowband filtering, and rapid atmospheric turbulence averaging -- enables next-generation gravitational tests, precision lunar geodesy, and improved lunar reference frames in support of planetary exploration.

Coherent radio emission mechanism of solar radio bursts is one of the most complicated and controversial topics in solar physics. To clarify the mechanism(s) of different types of solar radio bursts, (radio) wave excitation by energetic electrons in homogeneous plasmas has been widely studied via particle-in-cell (PIC) code numerical simulations. The solar corona is, however, inhomogeneous over almost all spatial scales. Inhomogeneities of the plasma could influence the emission properties of solar radio bursts. In this paper, we, hence, investigate effects of inhomogeneity (in the magnetic field, plasma density and temperature) of plasmas in the solar corona on radio wave emission by ring-beam distributed energetic electrons utilizing 2.5-dimensional PIC simulations. Both the beam and electron cyclotron maser (ECM) instabilities could be triggered with the presence of the energetic ring-beam electrons. The resultant spectrum of the excited electromagnetic waves presents a zebra-stripe pattern in the frequency space. The inhomogeneous density or temperature in plasmas influences the frequency bandwidth and location of these excited waves. Our results can, hence, help to diagnose the plasma properties at the emission sites of solar radio bursts. Applications of our results to the solar radio bursts with zebra-stripe pattern are discussed.

The Laser Interferometer Space Antenna (LISA) will soon detect gravitational waves (GWs) emitted by massive black hole (MBH) mergers. Some theoretical models have predicted transient electromagnetic (EM) emission from these mergers, enabling the association of LISA GW sources with their EM counterparts via telescope follow-up. However, the number of unrelated EM transients that might contaminate telescope searches for the true transient counterparts of LISA MBH mergers is unknown. We investigate the expected numbers of unrelated EM transients that will coincide with simulated LISA localization volumes of MBH mergers, as a function of the merger total mass and redshift. We find that the number of potential contaminants in LISA localization volumes drops to unity for mergers at $z \lesssim 0.8$ and at 1 hour before coalescence. After coalescence, the parameter space corresponding to a maximum of one potential contaminant expands to $z \lesssim 1.5$. In contrast, if the redshifts for all transients detected in LISA sky localization regions are not available, the number of potential contaminants increases by an average factor of $\sim100$, and never drops below unity. Overall, we expect the average number of contaminating transients in telescope follow-up of LISA MBH mergers to be non-negligible, especially without redshift information for the detected transients. We recommend that endeavors designing follow-up strategies of LISA events should focus on: (1) building large redshift catalogs for host galaxies, (2) developing robust real-time transient classification algorithms, (3) and coordinating telescope resources to obtain redshifts for candidate transient EM counterparts in a timely manner.

Context. Multiple systems are common in field stars, and the frequency is found to be higher in early evolutionary stages. Thus, the study of young multiple systems during the embedded stages is key to have a comprehensive understanding of star formation. In particular, how material accretes from the large-scale envelope into the inner region and how this flow interacts with the system physically and chemically has not yet been well characterized observationally. Aims. We aim to provide a snapshot of the forming protobinary system SVS13A, consisting of VLA4A and VLA4B. This includes clear pictures of its kinematic structures, physical conditions, and chemical properties. Methods. We conducted ALMA observations toward SVS13A targeting CH3CN and CH3-13CN J=12-11 K-ladder line emission with a high spatial resolution of ~30 au at a spectral resolution of ~0.08 km s-1 Results. We perform LTE radiative transfer models to fit the spectral features of the line emission. We find the two-layer LTE radiative model including dust absorption is essential to interpret the CH3CN and CH3-13-CN line emission. We identify two major and four small kinematic components, and derive their physical and chemical properties. Conclusions. We find a possible infalling signature toward the bursting secondary source VLA4A, which may be fed by an infalling streamer from the large-scale envelope. The mechanical heating in the binary system, as well as the infalling shocked gas, likely play a role in the thermal structure of the protobinary system. By accumulating mass from the streamer, the system might have experienced a gravitationally unstable phase before the accretion outburst. Finally, the derived CH3CN/CH3-13-CN ratio is lower than the canonical ratio in the ISM and is different between VLA4A and VLA4B.

Before a binary system enters into a common envelope (CE) phase, accretion from the primary star onto the companion star through Roche Lobe overflow (RLOF) will lead to the formation of an accretion disk, which may generate jets. Accretion before and during the CE may alter the outcome of the interaction. Previous studies have considered different aspects of this physical mechanism. Here we study the properties of an accretion disk formed via 3D hydrodynamic simulations of the RLOF mass transfer between a 7 M$_\odot$, red supergiant star and a 1.4 M$_\odot$, neutron star companion. We simulate only the volume around the companion for improved resolution. We use a 1D implicit MESA simulation of the evolution of the system during 30,000 years between the on-set of the RLOF and the CE to guide the binary parameters and the mass-transfer rate, while we simulate only 21 years of the last part of the RLOF in 3D using an ideal gas isothermal equation of state. We expect that a pre-CE disk under these parameters will have a mass of $\sim 5\times 10^{-3}$ M$_\odot$ and a radius of $\sim$40 R$_\odot$ with a scale height of $\sim$5 R$_\odot$. The temperature profile of the disk is shallower than that predicted by the formalism of Shakura and Sunyaev, but more reasonable cooling physics would need to be included. We stress test these results with respect to a number of physical and numerical parameters, as well as simulation choices, and we expect them to be reasonable within a factor of a few for the mass and 15% for the radius. We also contextualize our results within those presented in the literature, in particular with respect to the dimensionality of simulations and the adiabatic index. We consider what properties of magnetic fields and jets may be supported by our disk and discuss prospects for future work.

The flare ribbon and an associated filament eruption are diagnosed using O iv 1401.16 A, Si iv 1402.77 A, and Mg ii k 2796.35 A spectral lines provided by IRIS. The flare ribbons have downflow (redshifts) in all these lines, and this redshift decreases from the transition region to the chromosphere. While the overlapping region (flare-ribbon+filament rise/eruption is dominated by upflows(blueshifts) in all three spectral lines. We found an extremely blueshifted Si iv profile (i.e., blueshift around -180 km/s) in the overlapping region. The mean non-thermal velocity (v_nt) in the flare ribbons is higher in O iv than Si iv. While, in the overlapping region, O iv have lower v_nt than Si iv. Note that very high v_nt around 80 km/s (in Si iv) exists in this weak B-class flare. The Mg ii k line widths are almost the same in the flare ribbon and overlapping region but, they are extremely broad than previously reported. We found double peak profiles of Si iv and O iv in the overlapping region. Most probably, one peak is due to downflow (flare ribbon) and another due to upflow (filament rise/eruption). We report a high redshift of more than 150 km/s in the weak B-class flare. In some cases, both peaks show upflows which might be the result of the superposition of two different sources, i.e., overlapping of two different velocity distributions in the line of sight.

Recent work has suggested that, during reionisation, spatial variations in the ionising radiation field should produce enhanced Ly ${\alpha}$ forest transmission at distances of tens of comoving Mpc from high-redshift galaxies. We demonstrate that the Sherwood-Relics suite of hybrid radiation-hydrodynamical simulations are qualitatively consistent with this interpretation. The shape of the galaxy--Ly ${\alpha}$ transmission cross-correlation is sensitive to both the mass of the haloes hosting the galaxies and the volume averaged fraction of neutral hydrogen in the IGM, $\bar{x}_{\rm HI}$. The reported excess Ly ${\alpha}$ forest transmission on scales r ~ 10 cMpc at $\langle z \rangle \approx 5.2$ -- as measured using C IV absorbers as proxies for high-redshift galaxies -- is quantitatively reproduced by Sherwood-Relics at z = 6 if we assume the galaxies that produce ionising photons are hosted in haloes with mass $M_{\rm h}\geq 10^{10}~h^{-1}\,{\rm M}_\odot$. However, this redshift mismatch is equivalent to requiring $\bar{x}_{\rm HI}\sim 0.1$ at $z\simeq 5.2$, which is inconsistent with the observed Ly ${\alpha}$ forest effective optical depth distribution. We speculate this tension may be partly resolved if the minimum C IV absorber host halo mass at z > 5 is larger than $M_{\rm h}=10^{10}~h^{-1}\,{\rm M}_\odot$. After reionisation completes, relic IGM temperature fluctuations will continue to influence the shape of the cross-correlation on scales of a few comoving Mpc at $4 \leq z \leq 5$. Constraining the redshift evolution of the cross-correlation over this period may therefore provide further insight into the timing of reionisation.

New concepts for observing the gravitational waves (GWs) using a detector on the Moon, such as the Lunar Gravitational-wave Antenna (LGWA), have gained increasing attention. By utilizing the Moon as a giant antenna, the LGWA is expected to detect GWs in the frequency range from 1 millihertz (mHz) to several hertz, with optimal sensitivity in the decihertz band. Despite the debated formation and evolution channel of intermediate-mass black holes (IMBHs) with masses in the range of $[10^2, 10^5]\ {\rm M_\odot}$, binary systems containing at least one IMBH are widely believed to generate GWs spanning from mHz to a few Hz, making them a key scientific target for the LGWA. We explore the detectability of IMBH binaries with the LGWA in this work. The LGWA is more sensitive to nearby binaries (i.e. with redshift $z\lesssim0.5$) with the primary mass $m_1 \in [10^4, 10^5] \ {\rm M_\odot}$, while it prefers distant binaries (i.e. $z \gtrsim 5$) with $m_1 \in [10^3, 10^4] \ {\rm M_\odot}$. Considering a signal-to-noise ratio threshold of 10, our results imply that the LGWA can detect IMBH binaries up to $z \sim \mathcal{O}(10)$. We further show that the LGWA can constrain the primary mass with relative errors $\lesssim 0.1\%$ for binaries at $z \lesssim 0.5$. Furthermore, we show that the IMBH binaries at $z \lesssim 0.1$ can be used to constrain redshift with relative errors $\lesssim 10\%$, and those with $m_1 \in [10^4, 10^5] \ {\rm M_\odot}$ can be localized by the LGWA to be within $\mathcal{O} (10)$ $\rm deg^2$.

The pulsar wind nebula CTB 87 (G74.9+1.2) is one of the sources emitting $\gamma$-rays with energies higher than 10 TeV, as measured by the Very Energetic Radiation Imaging Telescope Array System telescope (VERITAS). In this study, we undertake a reanalysis of the GeV emission from the CTB 87 region, utilising $\sim$16 years of high-energy $\gamma$-ray data collected with the Fermi Large Area Telescope. In the energy range of 0.03--1 TeV, the spectrum can be adequately described by a power-law model with an index of 1.34 $\pm$ 0.18, and the integral energy flux is calculated to be (7.25 $\pm$ 1.36) $\times$ 10$^{-13}$ erg cm$^{-2}$ s$^{-1}$. Based on the multiband data, we have employed a time-dependent model to investigate the non-thermal emission properties of CTB 87. In the model, it is assumed that particles with broken power-law energy distributions are continuously injected into the nebula. This results in multiband non-thermal emission being produced by relativistic leptons via synchrotron radiation and inverse Compton processes. Furthermore, the model suggests an energy of approximately 2.4 PeV for the most energetic particle in the nebula.

We present the stellar properties of 2908 galaxies at 0.6 < z < 1.0 from the LEGA-C survey. We emphasize the importance of high signal-to-noise, high spectral resolution spectroscopy in the inference of stellar population properties of galaxies. We estimate the galaxy properties with the SED fitting code Prospector, by fitting spectroscopy and broadband photometry together, drawn from the LEGA-C DR3 and UltraVISTA catalogs respectively. We report a positive correlation between light-weighted ages and stellar velocity dispersion ($\sigma_\star$). The trend with $\sigma_\star$ is weaker for the mass-weighted ages and stellar metallicity ($Z_\star$). On average, quiescent galaxies are characterized by high $Z_\star$, they are \sim 1.1 Gyr older, less dusty, with steeper dust attenuation slopes compared to star-forming galaxies. Conversely, star-forming galaxies are characterized by significantly higher dust optical depths and shallower (grayer) attenuation slopes. Low mass (high mass) star-forming galaxies have lower (higher) $Z_\star$, while their stellar populations are on average younger (older). A key pragmatic result of our study is that a linear-space metallicity prior is preferable to a logarithmic-space one when using photometry alone, as the latter biases the posteriors downward. Spectroscopy greatly improves stellar population measurements and is required to provide meaningful constraints on age, metallicity, and other properties. Pairing spectroscopy with photometry helps resolving the dust-age-metallicity degeneracy, yielding more accurate mass- and light-weighted ages, with ages inferred from photometry alone suffering such large uncertainties. Stellar metallicities are constrained by our spectroscopy, but precise measurements remain challenging (and impossible with photometry alone), particularly in the absence of Mg and Fe lines redward of 5000 $Å$ in the observed spectrum.

We present a detailed analysis of the kinematics of SiO maser stars around the center of the Milky Way, Sagittarius A* (Sgr A*). We used the archive data in the SiO v=1, J=2-1 emission line obtained by the Atacama Large Millimeter/Submillimeter Array (ALMA) in 2017 and 2021 (#2016.1.00940.S, PI Darling, J. and #2019.1.00292.S, PI Paine, J.). We detected 37 SiO maser stars in the channel maps and derived their angular offsets relative to Sgr A* and LSR radial velocities. We derived the proper motions of 35 stars by comparing their angular offsets in the two epochs. The proper motions of Wolf-Rayet and O star in the Nuclear Star Cluster (NSC) are reported to be rather random, except for the co-moving clusters IRS13E and IRS13N (Tsuboi et al. 2022). However, the derived proper motions of SiO maser stars do not look completely random. The proper motions of the SiO maser stars show a tendency to lie along the Galactic plane. The proper motion amplitudes of SiO maser stars are larger than the LSR velocity amplitudes. We estimated the 3D motions from the proper motions and LSR velocities. Many 3D velocities are near to or larger than the upper limit velocities for Kepler orbits around Sgr A*, whose mass is assumed to be 4x10^6 Msun. These indicate that the SiO maser stars around Sgr A* are members of the Nuclear Stellar Disk rather than the NSC.

Black hole images obtained by very long baseline interferometry (VLBI) by the Event Horizon Telescope are a new tool for testing general relativity in super-strong gravitational fields. These images demonstrated a ring-like structure which can be explained as the black hole shadow image. To date, there are no reliable methods for determining the parameters of these ring-like structures, such as diameter, width, and asymmetry. In this paper, an algorithm for determining black hole image parameters is proposed using a Gaussian asymmetric ring as an example. Using the proposed method, the diameter and asymmetry parameters of the image of a supermassive black hole in the galaxy M87$^{*}$ were estimated based on observational data obtained by the Event Horizon Telescope group.

Turbulent, relativistic nonthermal plasmas are ubiquitous in high-energy astrophysical systems, as inferred from broadband nonthermal emission spectra. The underlying turbulent nonthermal particle acceleration (NTPA) processes have traditionally been modelled with a Fokker-Planck (FP) diffusion-advection equation for the particle energy distribution. We test FP-type NTPA theories by performing and analysing particle-in-cell (PIC) simulations of turbulence in collisionless relativistic pair plasma. By tracking large numbers of particles in simulations with different initial magnetisation and system size, we first test and confirm the applicability of the FP framework. We then measure the FP energy diffusion ($D$) and advection ($A$) coefficients as functions of particle energy $\gamma m c^2$, and compare their dependence to theoretical predictions. At high energies, we robustly find $D \sim \gamma^2$ for all cases. Hence, we fit $D = D_0 \gamma^2$ and find a scaling consistent with $D_0 \sim \sigma^{3/2}$ at low instantaneous magnetisation $\sigma(t)$, flattening to $D_0 \sim \sigma$ at higher $\sigma \sim 1$. We also find that the power-law index $\alpha(t)$ of the particle energy distribution converges exponentially in time. We build and test an analytic model connecting the FP coefficients and $\alpha(t)$, predicting $A(\gamma) \sim \gamma \log \gamma$. We confirm this functional form in our measurements of $A(\gamma,t)$, which allows us to predict $\alpha(t)$ through the model relations. Our results suggest that the basic second-order Fermi acceleration model, which predicts $D_0 \sim \sigma$, may not be a complete description of NTPA in turbulent plasmas. These findings encourage further application of tracked particles and FP coefficients as a diagnostic in kinetic simulations of various astrophysically relevant plasma processes like collisionless shocks and magnetic reconnection.

Cyanopolyynes are among the largest and most commonly observed interstellar Complex Organic Molecules in star-forming regions. They are believed to form primarily in the gas-phase, but their formation routes are not well understood. We present a comprehensive study of the gas-phase formation network of cyanobutadiyne, HC$_5$N, based on new theoretical calculations, kinetics experiments, astronomical observations, and astrochemical modeling. We performed new quantum mechanics calculations for six neutral-neutral reactions in order to derive reliable rate coefficients and product branching fractions. We also present new CRESU data on the rate coefficients of three of these reactions (C$_3$N + C$_2$H$_2$, C$_2$H + HC$_3$N, CN + C$_4$H$_2$) obtained at temperatures as low as 24 K. In practice, six out of nine reactions currently used in astrochemical models have been updated in our reviewed network. We also report the tentative detection of the $^{13}$C isotopologues of HC$_5$N in the L1544 prestellar core. We derived a lower limit of $^{12}$C/$^{13}$C > 75 for the HC$_5$N isotopologues, which does not allow to bring new constraints to the HC$_5$N chemistry. Finally, we verified the impact of the revised reactions by running the GRETOBAPE astrochemical model. We found good agreement between the HC$_5$N predicted and observed abundances in cold ($\sim$10 K) objects, demonstrating that HC$_5$N is mainly formed by neutral-neutral reactions in these environments. In warm molecular shocks, instead, the predicted abundances are a factor of ten lower with respect to observed ones. In this environment possessing an higher gas ionization fraction, we speculate that the contribution of ion-neutral reactions could be significant.

The shape of the Ly-$\alpha$ transmission in the near zone of the redshift $z=5.9896$ quasar ULAS J0148$+$0600 (hereafter J0148) is consistent with a damping wing arising from an extended neutral hydrogen island in the diffuse intergalactic medium (IGM). Here we use simulations of late-ending reionisation from Sherwood-Relics to assess the expected incidence of quasars with Ly-$\alpha$ and Ly-$\beta$ absorption similar to the observed J0148 spectrum. We find a late end to reionisation at $z=5.3$ is a necessary requirement for reproducing a Ly-$\alpha$ damping wing consistent with J0148. This occurs in $\sim3$ per cent of our simulated spectra for an IGM neutral fraction $\langle x_{\rm HI}\rangle=0.14$ at $z=6$. However, using standard assumptions for the ionising photon output of J0148, the a priori probability of drawing a simulated quasar spectrum with a Ly-$\alpha$ damping wing profile and Ly-$\alpha$ near zone size that simultaneously match J0148 is very low, $p<10^{-3}$. We speculate this is because the ionising emission from J0148 is variable on timescales $t<10^{5}\rm\,yr$, or alternatively that the Ly-$\alpha$ transmission in the J0148 near zone is impacted by the transverse proximity effect from nearby star-forming galaxies or undetected quasars. We also predict the IGM temperature should be $T\sim 4\times 10^{4}\rm\,K$ within a few proper Mpc of the Ly-$\alpha$ near zone edge due to recent HI and HeII photo-heating. Evidence for enhanced thermal broadening in the Ly-$\alpha$ absorption near the damping wing edge would provide further evidence that the final stages of reionisation are occurring at $z<6$.

Recently Usov's mechanism of pair creation on the surface of compact astrophysical objects has been revisited [1] with a conclusion that the pair creation rate was previously underestimated in the literature by nearly two orders of magnitude. Here we consider an alternative hypothesis of pair creation due to a perturbation of the surface of a compact object. Radial perturbation is induced in hydrodynamic velocity resulting in a microscopic displacement of the negatively charged component with respect to the positively charged one. The result depends on the ratio between the spatial scale of the perturbation $\lambda$ and the mean free path $l$. When $\lambda\sim l$ the perturbation energy is converted into a burst of electron-positron pairs which are created in collisionless plasma oscillations at the surface; after energy excess is dissipated electrosphere returns to its electrostatic configuration. When instead $\lambda\gg l$, the perturbation is thermalized, its energy is transformed into heat, and pairs are created continuously by the heated electrosphere. We discuss the relevant astrophysical scenarios.

Radio-loud active galactic nuclei (AGN) with their jets pointed close to our line of sight constitute the majority of extragalactic $\gamma$-ray sources and significantly contribute to the radiation observed in the even higher energy regime. The upcoming Cherenkov Telescope Array (CTA) is expected to detect fainter TeV objects, leading to an anticipated increase in the proportion of non-blazar extragalactic high-energy sources. Here we present the results of our dual-frequency (1.7 and 5~GHz) European VLBI Network (EVN) and enhanced Multi Element Remotely Linked Interferometer Network (e-MERLIN) observations of two faint radio sources from the list of TeV candidate sources. They do not show signs of nuclear activity in their optical spectra, but they were hypothesized to contain faint AGN that is outshone by the host galaxy. We used the mas-scale resolution radio data to try to pinpoint the location of the compact radio emitting feature, determine its spectral index, radio power, brightness temperature and radio-X-ray luminosity ratio and thus identify the origin of the radio emission. Our results suggest that both optically passive-looking galaxies host faint compact radio-emitting AGN with steep spectra.

The collision outcomes of dust aggregates in protoplanetary disks dictate how planetesimals form. Experimental and numerical studies have suggested that bouncing collisions occurring at low impact velocities may limit aggregate growth in the disks, but the conditions under which bouncing occurs have yet to be fully understood. In this study, we perform a suite of collision simulations of moderately compact dust aggregates with various impact velocities, aggregate radii, and filling factors ranging between 0.4 and 0.5. Unlike previous simulations, we generate compact aggregates by compressing more porous ones, mimicking the natural processes through which compact aggregates form. We find that the compressed aggregates bounce above a threshold mass, which decreases with impact velocity. The threshold mass scales with impact velocity as the $-4/3$ power, consistent with the findings of previous experiments. We also find that the threshold aggregate mass for bouncing depends strongly on filling factor, likely reflecting the strong filling-factor dependence of the compressive strength of compressed aggregates. Our energy analysis reveals that nearly 90\% of the initial impact energy is dissipated during the initial compression phase, and over 70\% of the remaining energy is dissipated during the subsequent stretching phase, regardless of whether the collision results in sticking or bouncing. Our results indicate that dust aggregates with a filling factor of 0.4 cease to grow beyond 100 $\mathrm{\mu m}$ as a result of the bouncing barrier.

1LHAASO J0249+6022 is an extended very-high-energy gamma-ray source discovered by the Large High-Altitude Air Shower Observatory. Based on nearly 16.1 years of data from the Fermi Large Area Telescope, we report the probable gamma-ray emission from 1LHAASO J0249+6022 in the 0.03-1 TeV energy range. The results show that its gamma-ray spectrum can be well fitted by a single power law with an index of 1.54 $\pm$ 0.17, and integral photon flux is (4.28 $\pm$ 1.03) $\times$ 10$^{-11}$ photons cm$^{-2}$ s$^{-1}$. We also considered theoretically whether the non-thermal emission could originate from a pulsar wind nebula (PWN) scenario. Assuming that the particles injected into the nebula have a power-law distribution, the resulting spectrum from the inverse Compton scattering is consistent with the detected GeV and TeV gamma-ray fluxes. Our study shows that the PWN scenario is reasonable for 1LHAASO J0249+6022.

Knowledge of the primordial matter density field from which the large-scale structure of the Universe emerged over cosmic time is of fundamental importance for cosmology. However, reconstructing these cosmological initial conditions from late-time observations is a notoriously difficult task, which requires advanced cosmological simulators and sophisticated statistical methods to explore a multi-million-dimensional parameter space. We show how simulation-based inference (SBI) can be used to tackle this problem and to obtain data-constrained realisations of the primordial dark matter density field in a simulation-efficient way with general non-differentiable simulators. Our method is applicable to full high-resolution dark matter $N$-body simulations and is based on modelling the posterior distribution of the constrained initial conditions to be Gaussian with a diagonal covariance matrix in Fourier space. As a result, we can generate thousands of posterior samples within seconds on a single GPU, orders of magnitude faster than existing methods, paving the way for sequential SBI for cosmological fields. Furthermore, we perform an analytical fit of the estimated dependence of the covariance on the wavenumber, effectively transforming any point-estimator of initial conditions into a fast sampler. We test the validity of our obtained samples by comparing them to the true values with summary statistics and performing a Bayesian consistency test.

Recent work with JWST has demonstrated its capability to identify and chemically characterize multiple populations in globular clusters down to the H-burning limit. In this study, we explore the kinematics of multiple populations in the globular cluster 47 Tucanae by combining data from JWST, HST, and Gaia. We analyzed velocity dispersion and anisotropy profiles from the cluster center out to $\sim$10$R_h$. Our findings indicate that while 1G stars are isotropic, 2G stars are significantly radially anisotropic. These results align with the predictions of simulations of the dynamical evolution of clusters where 2G stars are initially more centrally concentrated than 1G stars. Furthermore, we subdivided the 2G population into two subpopulations: $2G_A$ and $2G_B$, with the latter being more chemically extreme. We compared their dynamical profiles and found no significant differences. For the first time, we measured the degree of energy equipartition among the multiple populations of 47 Tucanae. Overall, within the analyzed radial range ($\sim$2-4$R_h$), both populations exhibit a low degree of energy equipartition. The most significant differences between 1G and 2G stars are observed in the tangential velocity component, where 2G stars are characterized by a stronger degree of energy equipartition than 1G stars. In the radial component, the behavior of 1G and 2G stars is more variable, with differences largely dependent on radius. Finally, our analysis reveals that the ratio of rotational velocity to velocity dispersion is larger for the 2G population, while 1G stars exhibit higher skewness in their tangential proper motions, providing further evidence of differences in the kinematic properties of the 1G and 2G populations.

Low surface brightness galaxies (LSBGs) are important for understanding galaxy evolution and cosmological models. The upcoming large-scale surveys are expected to uncover a large number of LSBGs, requiring accurate automated or machine learning-based methods for their detection. We study the scope of transfer learning for the identification of LSBGs. We use transformer models divided into two categories: LSBG Detection Transformer (LSBG DETR) and LSBG Vision Transformer (LSBG ViT), trained on Dark Energy Survey (DES) data, to identify LSBGs from dedicated Hyper Suprime-Cam (HSC) observations of the Abell 194 cluster, which are two magnitudes deeper than DES. The data from DES and HSC were standardized based on pixel-level surface brightness. We used two transformer ensembles to detect LSBGs. This was followed by a single-component Sérsic model fit and a final visual inspection to filter out potential false positives and improve sample purity. We present a sample of 171 low surface brightness galaxies (LSBGs) in the Abell 194 cluster using HSC data, including 87 new discoveries. Of these, 159 were identified using transformer models, and 12 additional LSBGs were found through visual inspection. The transformer model achieved a true positive rate (TPR) of 93% in HSC data without any fine-tuning. Among the LSBGs, 28 were classified as ultra-diffuse galaxies (UDGs). The number of UDGs and the radial UDG number density suggest a linear relationship between UDG numbers and cluster mass on a log scale. UDGs share similar Sérsic parameters with dwarf galaxies and occupy the extended end of the $R_{\mathrm{eff}}-M_g$ plane, suggesting they might be an extended subpopulation of dwarf galaxies. We have demonstrated that transformer models trained on shallower surveys can be successfully applied to deeper surveys with appropriate data normalization.

As demonstrated with the Sloan Digital Sky Survey (SDSS), Pan-STARRS, and most recently with Gaia data, broadband near-UV to near-IR stellar photometry can be used to estimate distance, metallicity, and interstellar dust extinction along the line of sight for stars in the Galaxy. Anticipating photometric catalogs with tens of billions of stars from Rubin's Legacy Survey of Space and Time (LSST), we present a Bayesian model and pipeline that build on previous work and can handle LSST-sized datasets. Likelihood computations utilize MIST/Dartmouth isochrones and priors are derived from TRILEGAL-based simulated LSST catalogs from P. Dal Tio et al. The computation speed is about 10 ms per star on a single core for both optimized grid search and Markov Chain Monte Carlo methods; we show in a companion paper by K. Mrakovčić et al. how to utilize neural networks to accelerate this performance by up to an order of magnitude. We validate our pipeline, named PhotoD (in analogy with photo-z, photometric redshifts of galaxies) using both simulated catalogs and SDSS, DECam, and Gaia photometry. We intend to make LSST-based value-added PhotoD catalogs publicly available via the Rubin Science Platform with every LSST data release.

We present \texttt{SACRA-2D}, a new MPI and OpenMP parallelized, fully relativistic hydrodynamics (GRHD) code in dynamical spacetime under axial symmetry with the cartoon method using the finite-volume shock-capturing schemes for hydrodynamics. Specifically, we implemented the state-of-the-art HLLC Riemann solver and found better accuracy than the standard Total Variation Diminishing Lax-Friedrich Riemann solver. The spacetime evolves under the Baumgarte-Shapiro-Shibata-Nakamura formalism with Z4c constraint propagation. We demonstrate the accuracy of the code with some benchmark tests and excellent agreement with other codes in the literature. A wide variety of test simulations, including the head-on collision of black holes, the migration and collapse of neutron stars, and the collapse of a rotating supermassive star to a massive black hole and a disk, is also performed to show the robustness of our new code.

Interstellar material has been discovered in our Solar System, yet its origins and details of its transport are unknown. Here we present $\alpha$ Centauri as a case study of the delivery of interstellar material to our Solar System. $\alpha$ Centauri is a mature triple star system that likely harbours planets and is moving towards us with the point of closest approach approximately 28,000 years in the future. Assuming a current ejection model for the system, we find that such material can reach our Solar System and may currently be present here. The material that does reach us is mostly a product of low ($<2$ km/s) ejection velocities, and the rate at which it enters our Solar System is expected to peak around the time of $\alpha$ Centauri 's closest approach. If $\alpha$ Centauri ejects material at a rate comparable to our own Solar System, we estimate the current number of $\alpha$ Centauri particles larger than 100 m in diameter within our Oort Cloud to be $10^{6}$, and during $\alpha$ Centauri 's closest approach, this will increase by an order of magnitude. However, the observable fraction of such objects remains low as there is only a probability of $10^{-6}$ that one of them is within 10 au of the Sun. A small number ($\sim 10$) meteors greater than 100 micrometers from $\alpha$ Centauri may currently be entering Earth's atmosphere every year: this number is very sensitive to the assumed ejected mass distribution, but the flux is expected to increase as $\alpha$ Centauri approaches.

NASA's New Horizon observations yielded the most accurate measurement of the cosmic optical background (COB) intensity. The reported COB flux is $16.37\pm 1.47$ $\mathrm{nW/m^2/sr}$ at a pivot wavelength $ \lambda_{piv} = 0.608~ \mu\mathrm{m}$ observed in the range $( 0.4 ~ \mu\mathrm{m} \lesssim \lambda \lesssim 0.9 \, \mu\mathrm{m})$. After subtracting the measured intensity from the deep Hubble space telescope count, an anomalous excess flux $8.06 \pm 1.92~\mathrm{nW/m^2/sr}$ has been found. This observation could hint toward decaying dark matter producing photons. In this work, we have considered sterile neutrinos of the keV scale as well as the eV scale decaying via sterile-to-sterile transition magnetic dipole moment and active-to-sterile transition magnetic dipole moment, respectively, to explain anomalous flux. The sterile neutrinos with a mass of $\mathcal{O}(\rm keV)$ with transition magnetic moment in the range $3\times 10^{-14}\,\mu_{B} - 10^{-13}\,\mu_{B}$, and mass of $\mathcal{O}(\rm eV)$ with transition magnetic moment in the range $3.33\times 10^{-13}\,\mu_{B} - 10^{-12}\,\mu_{B}$ can successfully account for the observed anomalous intensity.

The first discovered interstellar small object, `Oumuamua (1I/2017 U1), presents unique physical properties of extremely elongated geometric shape and dual characteristics of an asteroid and a comet. These properties suggest a possible origin through tidal fragmentation, which posits that `Oumuamua was produced through intensive tidal fragmentation during a close encounter with a star or a white dwarf, resulting in its shape and ejection from its natal system. According to this mechanism, a high initial orbit eccentricity and a small pericentre of the parent body are necessary to produce `Oumuamua-like objects. To verify whether this mechanism can occur in single giant planet systems, we conduct long-term numerical simulations of systems with a low-mass ($0.5M_\odot$) host star and a giant planet in this study. We determine that an eccentric orbit ($e_\mathrm{p}\sim0.2$) and a Jupiter-mass ($M_\mathrm{p}\sim M_\mathrm{J}$) of the planet appears to be optimal to generate sufficient perturbations for the production of `Oumuamua-like objects. When the planetary semi-major axis $a_\mathrm{p}$ increases, the proportion of planetesimals ejected beyond the system $P(\mathrm{ej})$ increases accordingly, while the possibilities of ejected planetesimals undergoing stellar tidal fragmentation $P(\mathrm{tidal}|\mathrm{ej})$ remains relatively constant at $\sim0.6\%$. Focusing on stellar tidal fragmentation alone, the ratio of extremely elongated interstellar objects to all interstellar objects is $P_\mathrm{e}\sim3\%$.

This study numerically investigates the dynamics of barred spiral galaxies using 3D Ferrers bar response models. A total of 708 models were analyzed, incorporating variations in the axisymmetric potential (nucleus, bulge, disk, halo), bar length, mass, angular velocity, and disk stellar velocity dispersion. Model evaluation employed the Spearman correlation (to assess input-output relationships) and permutation feature importance in a Random Forest Regressor (to measure input variable impacts). Orbital configurations of test particles reveal the critical role of bar dynamics in shaping galaxies' morphological and kinematic properties. Key findings emphasize how bar potential influences major orbital families, affecting barred galaxies' long-term structure. These results provide deeper insights into galactic component interactions and a robust framework for understanding bar properties.

Ultra-light axions are hypothetical scalar particles that influence the evolution of large-scale structures of the Universe. Depending on their mass, they can potentially be part of the dark matter component of the Universe, as candidates commonly referred to as fuzzy dark matter. While strong constraints have been established for pure fuzzy dark matter models, the more general scenario where ultra-light axions constitute only a fraction of the dark matter has been limited to a few observational probes. In this work, we use the galaxy cluster number counts obtained from the first All-Sky Survey (eRASS1) of the SRG/eROSITA mission together with gravitational weak lensing data from the Dark Energy Survey, the Kilo-Degree Survey, and the Hyper Suprime-Cam, to constrain the fraction of ultra-light axions in the mass range $10^{-32}$ eV to $10^{-24}$ eV. We put upper bounds on the ultra-light axion relic density in independent logarithmic axion mass bins by performing a full cosmological parameter inference. We find an exclusion region in the intermediate ultra-light axion mass regime with the tightest bounds reported so far in the mass bins around $m_\mathrm{a}=10^{-27}$ eV with $\Omega_\mathrm{a} < 0.0036$ and $m_\mathrm{a}=10^{-26}$ eV with $\Omega_\mathrm{a} < 0.0084$, both at 95% confidence level. When combining with CMB probes, these bounds are tightened to $\Omega_\mathrm{a} < 0.0030$ in the $m_\mathrm{a}=10^{27}$ eV mass bin and $\Omega_\mathrm{a} < 0.0058$ in the $m_\mathrm{a}=10^{-26}$ eV mass bin, both at 95% confidence level. This is the first time that constraints on ultra-light axions have been obtained using the growth of structure measured by galaxy cluster number counts. These results pave the way for large surveys, which can be utilized to obtain tight constraints on the mass and relic density of ultra-light axions with better theoretical modeling of the abundance of halos.

Dark Scattering (DS) is an interacting dark energy model characterised by pure momentum exchange between dark energy and dark matter. It is phenomenologically interesting because it is unconstrained by CMB data and can alleviate the $S_8$ tension. We derive constraints on cosmological and DS parameters using three two-point correlation functions (3$\times$2pt) from the Dark Energy Survey third year data release (DES Y3). We then add information from the multipoles of the galaxy power spectrum combined with Baryonic Acoustic Oscillation (BAO) measurements using the twelfth data release of the Baryon Oscillation Spectroscopic Survey (BOSS DR12) and external BAO measurements. We compare results from the direct combination of the probes with the joint posterior distribution calculated with a normalising flow approach. Additionally, we run a CMB analysis with the Planck Public Release 4 (PR4) for comparison of the cosmological constraints. Overall, we find that the combination of probes allows minimising the projection effects and improves constraints without the need to include CMB information. It brings the marginalised posterior maxima closer to the corresponding best-fit values and weakens the sensitivity to the priors of the spectroscopic modelling nuisance parameters. These findings are highly relevant in light of forthcoming data of surveys like DESI, Euclid, and Rubin.

The brightest ever observed gamma ray burst GRB 221009A at redshift $z = 0.151$ was detected on October 9, 2022. Its highest energy photons have been recorded by the LHAASO collaboration up to above $12 \, \rm TeV$, and one of the at ${\cal E} = 251 \, \rm TeV$ by the Carpet-2 collaboration. Very recently, the Carpet-3 collaboration has completed the data analysis, showing that the evidence of the $251 \, {\rm TeV}$ photon is quite robust. Still, according to conventional physics photons with ${\cal E} \gtrsim 10 \, \rm TeV$ cannot be observed owing to the absorption by the extragalactic background light (EBL). Previously it has been demonstrated that an axion-like particle (ALP) with allowed parameters ensures the observability of the LHAASO photons. Here we show that the Lorentz invariance violation allows the ${\cal E} = 251 \, {\rm TeV}$ (now around 300 TeV) Carpet photon to be detected.

We present the Dark Energy Spectroscopic Instrument (DESI) Strong Lens Foundry. We discovered $\sim 3500$ new strong gravitational lens candidates in the DESI Legacy Imaging Surveys using residual neural networks (ResNet). We observed a subset (51) of our candidates using the Hubble Space Telescope (HST). All of them were confirmed to be strong lenses. We also briefly describe spectroscopic follow-up observations by DESI and Keck NIRES programs. From this very rich dataset, a number of studies will be carried out, including evaluating the quality of the ResNet search candidates and lens modeling. In this paper we present our initial effort in these directions. In particular, as a demonstration, we present the lens model for DESI-165.4754-06.0423, with imaging data from HST, and lens and source redshifts from DESI and Keck NIRES, respectively. We use a fast, fully forward-modeling Bayesian pipeline, GIGA-Lens. This is the first time a strong lens with HST data, or any high resolution imaging, has been modeled using GPUs.

In cosmic web analysis, complementary to traditional cosmological probes, the extrema (e.g. peaks and voids) two-point correlation functions (2PCFs) are of particular interest for the study of both astrophysical phenomena and cosmological structure formation. However most previous studies constructed those statistics via N-body simulations without a robust theoretical derivation from first principles. A strong motivation exists for analytically describing the 2PCFs of these local extrema, taking into account the nonlinear gravitational evolution in the late Universe. In this paper, we derive analytical formulae for the power spectra and 2PCFs of 2D critical points, including peaks (maxima), voids (minima) and saddle points, in mildly non-Gaussian weak gravitational lensing fields. We apply a perturbative bias expansion to model the clustering of 2D critical points. We successfully derive the power spectrum of weak lensing critical points up to the next-to-next-to-leading order (NNLO) in gravitational perturbation theory, where trispectrum configurations of the weak lensing field have to be included. We numerically evaluate those power spectra up to the next-to-leading order (NLO), which correspond to the inclusion of bispectrum configurations, and transform them to the corresponding 2PCFs. An exact Monte Carlo (MC) integration is performed assuming a Gaussian distributed density field to validate our theoretical predictions. Overall, we find similar properties in 2D compared to the clustering of 3D critical points previously measured from N-body simulations. Contrary to standard lensing power spectra analysis, we find distinct BAO features in the lensing peak 2PCFs due to the gradient and curvature constraints, and we quantify that non-Gaussianity makes for ~10% of the signal at quasi-linear scales which could be important for current stage-IV surveys.

Gravitational waves from merging binary neutron stars carry characteristic information about their astrophysical properties, including masses and tidal deformabilities, that are needed to infer their radii. In this study, we use Bayesian inference to quantify the precision with which radius can inferred with upgrades in the current gravitational wave detectors and next-generation observatories such as the Einstein Telescope and Cosmic Explorer. We assign evidences for a set of plausible equations of state, which are then used as weights to obtain radius posteriors. We find that prior choices and the loudness of observed signals limit the precision and accuracy of inferred radii by current detectors. In contrast, next-generation observatories can resolve the radius precisely and accurately, across most of the mass range to within $\lesssim 5\%$ for both soft and stiff equations of state. We also explore how the choice of the neutron star mass prior can influence the inferred masses and potentially affect radii measurements, finding that choosing an astrophysically motivated prior does not notably impact an individual neutron star's radius measurements.

Context. The combination between non-thermal and thermal emission in $\gamma$-ray blazars pushes them to a specific region of the mid-infrared three-dimensional color diagram, the so-called blazar locus, built based on observations performed with the Wide-field Infrared Survey Explorer. The selection of blazar candidates based on these mid-infrared colors has been extensively used in the past decade in the hunt for the counterparts of unassociated $\gamma$-ray sources observed with the Fermi Large Area Telescope and in the search for new blazars in optical spectroscopic campaigns. Aims. In this work, we provide a theoretical description of the origin of the blazar locus and show how we can reasonably reproduce it with a model consisting of only three spectral components: a log-parabola accounting for the non-thermal emission, and an elliptical host and dust torus accounting for the thermal emission. Methods. We simulate spectral energy distributions for blazars, starting with a pure log-parabola model and then increasing its complexity by adding a template elliptical galaxy and dust torus. From these simulations, we compute the mid-infrared magnitudes and corresponding colors to create our own version of the blazar locus. Results. Our modeling allows for the selection of spectral parameters that better characterize the mid-infrared emission of $\gamma$-ray blazars, such as the log-parabola curvature ($\beta < 0.04$ for 50\% of our sample) and an average spectral peak around $E_p \approx 1.5 \times 10^{-13}$ erg. We also find that the log-parabola is the main spectral component behind the observed mid-infrared blazar colors, although additional components such as a host galaxy and a dust torus are crucial to obtain a precise reconstruction of the blazar locus.

The spinning electromagnetic universe, known also as the Rotating Bertotti-Robinson(RBR) spacetime is considered as a model to represent our cosmos. The model derives from different physical considerations, such as colliding waves, throat region, and near horizon geometry of the Kerr-Newman black hole. Our interest is whether such a singularity-free spinning cosmology gives rise to a natural direction of flow, a 'chirality' for charged particles. Homochiral structures are known to be crucial for biology to start. Our concern here is cosmology rather than biology, but as in biology, the stable structures in cosmology may also rely on homochiral elements. We show the occurrence of closed timelike curves a 'la' G{ö}del. Such curves, however, seem possible only at localized cell structures, not at large scales, but according to our prescription of near horizon geometry, they arise in the vicinity of any charged, spinning black hole.

A universal theory of linear instabilities in swirling flows, occurring in both natural settings and industrial applications, is formulated. The theory encompasses a wide range of open and confined flows, including spiral isothermal flows and baroclinic flows driven by radial temperature gradients and natural gravity in rotating fluids. By employing short-wavelength local analysis, the theory generalizes previous findings from numerical simulations and linear stability analyses of specific swirling flows, such as spiral Couette flow, spiral Poiseuille flow, and baroclinic Couette flow. A general criterion, extending and unifying existing criteria for instability to both centrifugal and shear-driven perturbations in swirling flows is derived, taking into account viscosity and thermal diffusion, and guiding experimental and numerical investigations in the otherwise inaccessible parameter regimes.

Massive states produce higher derivative corrections to Einstein gravity in the infrared, which are encoded into operators of the Effective Field Theory (EFT) of gravity. These EFT operators modify the geometry and affect the tidal properties of black holes, either neutral or charged. A thorough analysis of the perturbative tidal deformation problem leads us to introduce a tidal Green function, which we use to derive two universal formulae that efficiently provide the constant and running Love numbers induced by the EFT. We apply these formulae to determine the tidal response of EFT-corrected non-spinning black holes induced by vector and tensor fields, reproducing existing results where available and deriving new ones. We find that neutral black hole Love numbers run classically for $l\geq 3$ while charged ones run for $l\geq2$. Insights from the Frobenius method and from EFT principles confirm that the Love number renormalization flow is a well-defined physical effect. We find that extremal black holes can have Love numbers much larger than neutral ones, up to ${\cal O}(1)$ within the EFT validity regime, and that the EFT cutoff corresponds to the exponential suppression of the Schwinger effect. We discuss the possibility of probing an Abelian dark sector through gravitational waves, considering a scenario in which dark-charged extremal black holes exist in the present-day Universe.

The interaction of gravitational waves (GWs) with matter is thought to be typically negligible in the Universe. We identify a possible exception in the case of resonant interactions, where GWs emitted by a background binary system such as an inspiraling supermassive black hole (SMBH) binary causes a resonant response in a stellar-mass foreground binary and the frequencies of the two systems become, and remain, synchronized. We point out that such locking is not only possible, but can significantly reduce the binaries' merger time for $\mathcal{O}(1-10^4)$ binaries in the host galaxy of the merging SMBHs of $10^{9-11}M_{\odot}$ for standard general relativity and even more either in ``wet'' SMBH mergers or in certain modified theories of gravity where the inspiral rate is reduced. This could leave an imprint on the period distribution of stellar mass binaries in post-merger galaxies which could be detectable by future GW detectors, such as LISA.

In this article, we present a description of the behaviour of shock-compressed solid materials following the Geometrical Shock Dynamics (GSD) theory. GSD has been successfully applied to various gas dynamics problems, and here we have employed it to investigate the propagation of cylindrically and spherically symmetric converging shock waves in solid materials. The analytical solution of shock dynamics equations has been obtained in strong-shock limit, assuming the solid material to be homogeneous and isotropic and obeying the Mie-Gruneisen equation of state. The non-dimensional expressions are obtained for the velocity of shock, the pressure, the mass density, the particle velocity, the temperature, the speed of sound, the adiabatic bulk modulus, and the change-in-entropy behind the strong converging shock front. The influences as a result of changes in (i) the propagation distance r from the axis or centre (r=0) of convergence, (ii) the Gruneisen parameter, and (iii) the material parameter are explored on the shock velocity and the domain behind the converging shock front. The results show that as the shock focuses at the axis or origin, the shock velocity, the pressure, the temperature, and the change-in-entropy increase in the shock-compressed titanium Ti6Al4V, stainless steel 304, aluminum 6061-T6, etc.

Planetary and stellar convection, which are compressible and turbulent, remain poorly understood. In this paper, we report numerical results on the scaling of Nusselt number ($\mathrm{Nu}$) and Reynolds number ($\mathrm{Re}$) for extreme convection. Using computationally-efficient MacCormack-TVD finite difference method, we simulate compressible turbulent convection in a two-dimensional Cartesian box up to $\mathrm{Ra} = 10^{16}$, the highest $\mathrm{Ra}$ achieved so far, and in a three-dimensional box up to $\mathrm{Ra} = 10^{11}$. We show adiabatic temperature drop in the bulk flow, leading to the Reynolds number scaling $\mathrm{Ra}^{1/2}$. More significantly, we show classical $1/3$ Nusselt number scaling: $\mathrm{Nu} \propto \mathrm{Ra}^{0.32}$ in 2D, and $\mathrm{Nu} \propto \mathrm{Ra}^{0.31}$ in 3D up to the highest $\mathrm{Ra}$.

Ball lighting (BL) has been observed for centuries. There are large number of books, review articles, and original scientific papers devoted to different aspects of BL phenomenon. Yet, the basic features of this phenomenon have never been explained by known physics. The main problem is the source which could power the dynamics of the BL. We advocate an idea that the dark matter in form of the axion quark nuggets (AQN) made of standard model quarks and gluons (similar to the old idea of the Witten's strangelets) could internally generate the required power. The corresponding macroscopically large object in form of the AQN behaves as {\it chameleon}: it does not interact with the surrounding material in dilute environment and serves as perfect cold DM candidate. However, AQN becomes strongly interacting object in sufficiently dense environment. The AQN model was invented long ago without any relation to the BL physics. It was invented with a single motivation to explain the observed similarity $\Omega_{\rm DM}\sim \Omega_{\rm visible}$ between visible and DM components. This relation represents a very generic feature of this framework, not sensitive to any parameters of the construction. However, with the same set of parameters being fixed long ago this model is capable to address the key elements of the BL phenomenology, including the source of the energy powering the BL events. In particular, we argue that the visible size of BL, its typical life time, the frequency of appearance , etc are all consistent with suggested proposal when BL represents a profound manifestation of the DM physics represented by the AQN objects. We also formulate a unique possible test which can refute or unambiguously substantiate this unorthodox proposal on nature of BL.

We demonstrate that in the presence of an external magnetic field, an uncharged classical Schwarzschild black hole moving superluminally in a dielectric with permittivity $\epsilon > 1$ produces Cherenkov emission. This is a new physical effect: classical (non-quantum) emission of electromagnetic waves by a completely charge-neutral ``particle''. The governing equations (involving General Relativity, electromagnetism, and the physics of continuous media) have no external electromagnetic source - it is the distortion of the initial electromagnetic fields by the gravity of the black hole that plays the role of a superluminally moving source. The effect relies on nonzero values of both the magnetic field and the gravitational radius, as well as on the usual Cherenkov condition on the velocity, $v/c > 1/\sqrt{\epsilon}$. Unlike Cherenkov emission by a point charge, the effective source in this case is spatially distributed, with emission generated along the single Cherenkov emission cone. The emitted spectrum is red-dominated, with power $\propto dk_z /|k_z|$ for wave numbers $|k_z| \leq 1/R_G$, where $R_G$ is the Schwarzschild radius. We comment on possible observability of this process during black hole -- neutron star mergers.

The structure formation in the local Universe is considered within the weak-field modification of General Relativity involving the cosmological constant. This approach enables to describe the dynamics of groups and clusters of galaxies, to explain the discrepancy in the observational properties of the local (late) and the global (early) Universe, i.e. the Hubble tension as a result of two flows, local and global ones, with non-equal Hubble parameters. The kinetic analysis with the modified gravitational potential involving the cosmological constant is shown to predict semi-periodical structure of filaments in the local universe. In the local scale this complements the Zeldovich pancake theory of evolution of the primordial density perturbations and of structure formation in the cosmological scale. The role of the cosmological constant is outlined in rescaling of the physical constants from one aeon to another within the Conformal Cyclic Cosmology.

The commonly used quasilinear approximation allows one to calculate the turbulent transport coefficients for the mean of a passive scalar or a magnetic field in a given velocity field. Formally, the quasilinear approximation is exact when the correlation time of the velocity field is zero. We calculate the lowest-order corrections to the transport coefficients due to the correlation time being nonzero. For this, we use the Furutsu-Novikov theorem, which allows one to express the turbulent transport coefficients in a Gaussian random velocity field as a series in the correlation time. We find that the turbulent diffusivities of both the mean passive scalar and the mean magnetic field are suppressed. Nevertheless, contradicting a previous study, we show that the turbulent diffusivity of the mean magnetic field is smaller than that of the mean passive scalar. We also find corrections to the $\alpha$ effect.

We present a description of the electronic structure of xenon within the density-functional theory formalism with the goal of accurately modeling dark matter-induced ionisation in liquid xenon detectors. We compare the calculated electronic structures of the atomic, liquid and crystalline solid phases, and find that the electronic charge density and its derivatives in momentum space are similar in the atom and the liquid, consistent with the weak interatomic van der Waals bonding. The only notable difference is a band broadening of the highest occupied $5p$ levels, reflected in the densities of states of the condensed phases, as a result of the inter-atomic interactions. We therefore use the calculated density of states of the liquid phase, combined with the standard literature approach for the isolated atom, to recompute ionisation rates and exclusion limit curves for the XENON10 and XENON1T experiments. We find that the broadening of the 5$p$ levels induced by the liquid phase is relevant only at low dark matter masses, where it increases the ionisation rate relative to that of the isolated atom. For most of the probable mass range the energies of the discrete 4$d$ and 5$s$ levels have the strongest effect on the rate. Our findings suggest a simple scheme for calculating dark matter-electron scattering rates in liquid noble gas detectors, using the calculated values for the atom weighted by the density of states of the condensed phase.

In the last decade, the ringdown community has made large strides in understanding the aftermath of binary black hole mergers through the study of numerical simulations. In this note, we introduce two flavors of fitting algorithms, that have been verified against each other, for the extraction of quasinormal mode amplitudes from ringdown waveforms - $\texttt{qnmfits}$ in Python and $\texttt{KerrRingdown}$ in Mathematica.

This work examines the Hubble constant (\(H_0\)) tension within the frameworks of perturbed \(f(R)\) gravity and perturbed \(f(R)\) gravity coupled with neutrinos, using lastest observational data. The datasets incorporate the Cosmic Microwave Background (CMB), Baryon Acoustic Oscillations (BAO), Cosmic Chronometers (CC), lensing, and Pantheon supernovae. We compare the ability of these models to bridge the discrepancy between Planck 2018 (\(H_0 = 67.4 \pm 0.5 \ \text{km/s/Mpc}\)) and the local R22 measurement (\(H_0 = 73.5 \pm 1.04 \ \text{km/s/Mpc}\)). In perturbed \(f(R)\) gravity, the derived \(H_0\) values align closely with Planck, leaving a substantial tension with R22. The inclusion of neutrino interactions introduces additional parameters that shift \(H_0\) toward higher values, reducing the tension with local measurements. Notably, the coupled model achieves a smaller residual tension compared to the standalone perturbed \(f(R)\) model, indicating that neutrino physics plays a significant role in modifying the late-time expansion dynamics. While both models provide insights into addressing the Hubble tension, the coupled \(f(R)\) gravity with neutrinos offers a more consistent alignment across the datasets.

Domain walls formed during a phase transition in a simple field theory model with $\mathbb{Z}_2$ symmetry in a periodic box have been demonstrated to annihilate as fast as causality allows and their area density scales $\propto t^{-1}$. We have performed numerical simulations of the dynamics of domain walls in the Two-Higgs Doublet Model (2HDM) where the potential has $\mathbb{Z}_2$ symmetry in two spatial dimensions. We observed significant differences with the standard case. Although the extreme long-time limit is the same for the $\approx 10^{5}$ sets of random initial configurations analysed, the percolation process is much slower due to the formation of long-lived loops. We suggest that this is due to the build up of superconducting currents on the walls which could lead ultimately to stationary configurations known as Kinky Vortons. We discuss the relevance of these findings for the production of Vortons in three spatial dimensions.