Papers for Tuesday, Oct 17 2023

We report galaxy MACS0416-Y3 behind the lensing cluster MACSJ0416.1--2403 as a tentative rotating disk at $z=8.34$ detected through its [OIII]$\lambda5007$ emission in JWST NIRCam wide-field slitless spectroscopic observations. The discovery is based on our new grism dynamical modeling methodology for JWST NIRCam slitless spectroscopy, using the data from ``Median-band Astrophysics with the Grism of NIRCam in Frontier Fields'' (MAGNIF), a JWST Cycle-2 program. The [OIII]$\lambda5007$ emission line morphology in grism data shows velocity offsets compared to the F480M direct imaging, suggestive of rotation. Assuming a geometrically thin disk model, we constrain the rotation velocity of $v_{\rm rot}=58^{+53}_{-35}$ km s$^{-1}$ via forward modeling of the two-dimensional (2D) spectrum. We obtain the kinematic ratio of $v_{\rm rot}/\sigma_v=1.6^{+1.9}_{-0.9}$, where $\sigma_v$ is the velocity dispersion, in line with a quasi-stable thin disk. The resulting dynamical mass is estimated to be $\log(M_{\rm dyn}/M_{\odot})=8.4^{+0.5}_{-0.7}$. If the rotation confirmed, our discovery suggests that rotating gaseous disks may have already existed within 600 million years after Big Bang.

Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code Cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (SHAP) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behavior of cooling and heating functions.

Emission line galaxies (ELGs) are now the preeminent tracers of large-scale structure at z>0.8 due to their high density and strong emission lines, which enable accurate redshift measurements. However, relatively little is known about ELG evolution and the ELG-halo connection, exposing us to potential modeling systematics in cosmology inference using these sources. In this paper, we propose a physical picture of ELGs and improve ELG-halo connection modeling using a variety of observations and simulated galaxy models. We investigate DESI-selected ELGs in COSMOS data, and infer that ELGs are rapidly star-forming galaxies with a large fraction exhibiting disturbed morphology, implying that many of them are likely to be merger-driven starbursts. We further postulate that the tidal interactions from mergers lead to correlated star formation in central-satellite ELG pairs, a phenomenon dubbed "conformity." We argue for the need to include conformity in the ELG-halo connection using galaxy models such as IllustrisTNG, and by combining observations such as the DESI ELG auto-correlation, ELG cross-correlation with Luminous Red Galaxies (LRGs), and ELG-cluster cross-correlation. We also explore the origin of conformity using the UniverseMachine model and elucidate the difference between conformity and the well-known galaxy assembly bias effect.

In this study, we compare the structural parameters of Ultra-Diffuse Galaxies (UDGs) to those of other dwarf galaxies and investigate whether UDGs form a distinct population. We observed deep $u'$-, $g'$-, and $r'$-band images (maximum limiting surface brightness [3$\sigma$, $10"\times10"$] u' and g': $\mathrm{\approx 30\,mag\,arcsec^{-2}}$; r': $\mathrm{\approx 29\,mag\,arcsec^{-2}}$) of Abell 1656 (Coma cluster) and Abell 262 with the Wendelstein Wide Field Imager at the 2.1m Fraunhofer telescope on the Wendelstein Observatory. We measure $u'-g'$ and $g'-r'$ colors and structural parameters using parametric fitting of tens of thousands of potential UDGs and other dwarf galaxies. Cluster members are identified and separated from diffuse background galaxies based on red sequence membership and location in the $u'-g'$ vs. $g'-r'$ color-color diagram. We find 11 UDGs in Abell 262 and 48 UDGs in Abell 1656. The latter is 6 times more than van Dokkum et al. found in the overlapping region. By comparing the structural parameters of UDGs to non-UDGs in our sample and to spheroidals from the literature, we do not find any separation in all tested parameter spaces. Instead, UDGs form the diffuse end of the already well-known spheroidal population and slightly extend it. Furthermore, we find that the UDG definition used by Koda et al. and Yagi et al. mainly extends the definition by van Dokkum et al. toward ordinary spheroidals.

This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample N=50,000 extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties and escape information are provided for all simulated stars. Stellar orbits are integrated in seven different static and time-varying Milky Way gravitational potential models where the structure of the disc, perturbations from the Large Magellanic Cloud and the mass and sphericity of the Milky Way's dark matter halo are all investigated. We find that the action coordinates of the mock extra-tidal stars are largely Galactic model independent, where minor offsets and broadening of the distributions between models are likely due to interactions with substructure. Importantly, we also report the first evidence for stellar stream contamination by globular cluster core stars and binaries for clusters with pericentre radii larger than five kiloparsecs. Finally, we provide a quantitative tool that uses action coordinates to match field stars to host clusters with probabilities. Ultimately, combining data from the GEMS catalogue with information of observed stars will allow for association of extra-tidal field stars with any Galactic globular cluster; a requisite tool for understanding population-level dynamics and evolution of clusters in the Milky Way.

Most ultra hot Jupiters (UHJs) show evidence of temperature inversions, in which temperature increases with altitude over a range of pressures. Temperature inversions can occur when there is a species that absorbs the stellar irradiation at a relatively high level of the atmospheres. However, the species responsible for this absorption remains unidentified. In particular, the UHJ KELT-20b is known to have a temperature inversion. Using high resolution emission spectroscopy from LBT/PEPSI we investigate the atomic and molecular opacity sources that may cause the inversion in KELT-20b, as well as explore its atmospheric chemistry. We confirm the presence of Fe I with a significance of 17$\sigma$. We also report a tentative $4.3\sigma$ detection of Ni I. A nominally $4.5\sigma$ detection of Mg I emission in the PEPSI blue arm is likely in fact due to aliasing between the Mg I cross-correlation template and the Fe I lines present in the spectrum. We cannot reproduce a recent detection of Cr I, while we do not have the wavelength coverage to robustly test past detections of Fe II and Si I. Together with non-detections of molecular species like TiO, this suggests that Fe I is likely to be the dominant optical opacity source in the dayside atmosphere of KELT-20b and may be responsible for the temperature inversion. We explore ways to reconcile the differences between our results and those in literature and point to future paths to understand atmospheric variability.

A recent article on high-resolution 86 GHz observations with the Global Millimeter VLBI Array, the phased Atacama Large Millimeter/submillimeter Array, and the Greenland Telescope describes the detection of a limb-brightened cylindrical jet, $25 \mu\rm{as}< z< 100 \mu\rm{as}$, where $z$ is the axial displacement from the supermassive black hole in the sky plane. It was shown to be much wider and much more collimated than 2D simulations of electromagnetic (Blandford-Znajek) jets from the event horizon predicted. This was an unanticipated discovery. The claimed detection of a jet connected to the accretion flow provides a direct observational constraint on the geometry and physics of the jet launching region for the first time in any black hole jetted system. This landmark detection warrants further analysis. This Letter focuses on the most rudimentary properties, the shape and size of the source of the detected jet emission, the determination of which is not trivial due to line-of-sight effects. Simple thick-walled cylindrical shell models for the source were analyzed to constrain the thickness of the jet wall. The analysis indicates a tubular jet source with a radius $R\approx 144 \mu\rm{as}\approx 38M$ and that the tubular jet walls have a width $W \approx 36\mu\rm{as} \approx 9.5 M$, where $M$ is the geometrized mass of the black hole (a volume comparable to that of the interior cavity). The observed cylindrical jet connects continuously to the highly limb-brightened jet (previously described as a thick-walled tubular jet) that extends to $z> 0.65$ mas, and the two are likely in fact the same outflow (i.e., from the same central engine).

This letter reports the first JWST spectroscopy of a white dwarf debris disk, giving a preliminary assessment of the salient features, and recommendations for future observations. The polluted and dusty star WD 0145+234 experienced a major collisional event in 2018, accompanied by an infrared outburst, and subsequently a gradual decrease in thermal emission. Time-series NIRSPEC observations demonstrate that the circumstellar disk is returning to a quiescent state with a T~1000 K infrared excess similar to the bulk of known dusty white dwarfs. MIRI spectroscopy reveals a 9-12 micron solid-state emission feature consistent with silicate minerals as observed in debris disks observed with Spitzer IRS. The strength and morphology of the silicate feature appear unchanged relative to the continuum in spectra taken over a year apart, consistent with steady-state collisional evolution of the circumstellar debris. A tentative emission feature around 7 microns may be due to carbonates, and if confirmed would indicate aqueous alteration in the parent body.

We present the first version release of SESAMME, a public, Python-based full spectrum fitting tool for Simultaneous Estimates of Star-cluster Age, Metallicity, Mass, and Extinction. SESAMME compares an input spectrum of a star cluster to a grid of stellar population models with an added nebular continuum component, using Markov chain Monte Carlo (MCMC) methods to sample the posterior probability distribution in four dimensions: cluster age, stellar metallicity $Z$, reddening $E(B-V)$, and a normalization parameter equivalent to a cluster mass. SESAMME is highly flexible in the stellar population models that it can use to model a spectrum; our testing and initial science applications use both BPASS and Starburst99. We illustrate the ability of SESAMME to recover accurate ages and metallicities even at a moderate signal-to-noise ratio (S/N ~3 - 5 per wavelength bin) using synthetic, noise-added model spectra of young star clusters. Finally, we test the consistency of SESAMME with other age and metallicity estimates from the literature using a sample of HST/COS far-UV spectra towards young, massive clusters in M83 and NGC 1313. We find that, on the whole, SESAMME infers star cluster properties that are consistent with the literature in both low- and high-metallicity environments.

Using the photometric population prediction method Red Dragon, we characterize the Red Sequence (RS) and Blue Cloud (BC) of DES galaxies in the COSMOS field. Red Dragon (RD) uses a redshift-evolving, error-corrected Gaussian mixture model to detail the distribution of photometric colors, smoothly parameterizing the two populations with relative weights, mean colors, intrinsic scatters, and inter-color correlations. This resulting fit of RS and BC yields RS membership probabilities $P_{\rm RS}$ for each galaxy. Even when training on only DES main bands $griz$, RD selects the quiescent population (galaxies with $\lg {\rm sSFR \cdot yr} < -11$) with $\gtrsim 90\%$ balanced accuracy out to $z=2$; augmenting with extended photometry from VIRCAM improves this accuracy to $\sim 95\%$ out to $z=3$. We measure redshift evolution of sSFR and galactic age in several stellar mass bins, finding that the BC is consistently more star-forming (by $\gtrsim 1~{\rm dex}$) and typically younger (by $\gtrsim 1~{\rm Gyr}$) than the RS (up to $z \sim 1.4$). This characterization of both RS and BC as functions of redshift and stellar mass improves our understanding of both populations and opens the door to more precise galaxy population characterization in future deep optical and IR systems.

The Milky Way's (MW) inner stellar halo contains an [Fe/H]-rich component with highly eccentric orbits, often referred to as the "last major merger." Hypotheses for the origin of this component include Gaia-Sausage/Enceladus (GSE), where the progenitor collided with the MW proto-disk 8-11 Gyr ago, and the Virgo Radial Merger (VRM), where the progenitor collided with the MW disk within the last 3 Gyr. These two scenarios make different predictions about observable structure in local phase space, because the morphology of debris depends on how long it has had to phase mix. The recently-identified phase-space folds in Gaia DR3 have positive caustic velocities, making them fundamentally different than the phase-mixed chevrons found in simulations at late times. Roughly 20% of the stars in the prograde local stellar halo are associated with the observed caustics. We compare the observed phase-space distribution to a time-series of FIRE-2 Latte simulations of a GSE-like merger, using a quantitative metric (2D causticality) that measures how phase-mixed a given distribution is. We find that the observed local phase-space distribution best matches the simulated data within 1 Gyr after the merger collides with the host galaxy disk, and certainly not later than 3 Gyr after collision. This is further evidence that the progenitor of the "last major merger" did not collide with the MW proto-disk at early times, as is thought for the GSE, but instead collided with the MW disk within the last few Gyr, consistent with the body of work surrounding the VRM.

Gaussian decomposition of HI4PI profiles was done for two regions along high velocity Complex M, the first with corresponding H-alpha emission and the second without. In the former, we detected the Critical Ionization Velocity (CIV) signature, a line width of 25 km/s, of H-alpha. In low-velocity gas, the CIV effect ionizes the atoms and produces the signatures of He (34 km/s), CNO (13.4 km/s), and the metals (around 6 km/s). In anomalous-velocity gas, the CIV effect ionizes the ions and produces the signatures of ionized CNO (20 km/s) and-or H-alpha (25.5 km/s), as well as the CIV signatures of metals around 6 km/s. These distinct CIV signatures indicate that the physical conditions are different in each velocity regime. The particles in the low-velocity gas cloud are initially neutral, the temperatures can be cold (50-100 K), and the abundances can be cosmic. The densities need to be high enough to impose regular recombination. The anomalous-velocity gas, on the other hand, initially contains a significant population of ionized and excited particles. The temperature and abundance values may be similar to those that characterize low-density gas, but the densities must be low enough so that the plasma is essentially collision-less, similar to the conditions in the solar wind.

We show that the inferred properties of the meteor CNEOS-2014-01-08 (IM1) can be naturally explained by tidal disruption of rocky planets on highly eccentric orbits around the most common stars, M-dwarfs. Rocky planets can develop a high orbital eccentricity as a result of a secular torque from an outer giant planet or a binary companion. Melting of their rocky crust during many tight periapse passages may result in additional elemental differentiation, leading to enhanced abundances of beryllium, lanthanum and uranium in the crust, as inferred from the composition of unique spherules along IM1's path. The excess energy imparted to the tidal stream in highly eccentric orbits naturally accounts for the inferred IM1 speed of $\sim 60~{\rm km~s^{-1}}$ relative to the Local Standard of Rest. Finally, the tidal disruption of $\sim 10M_\oplus$ reservoirs of rocky material around M-dwarfs accounts for an impact rate of 500kg meteors on Earth of once per decade, consistent with IM1.

The primary goal of the search for extraterrestrial intelligence (SETI) is to gain an understanding of the prevalence of technologically advanced beings (organic or inorganic) in the Galaxy. One way to approach this is to look for technosignatures: remotely detectable indicators of technology, such as temporal or spectral electromagnetic emissions consistent with an artificial source. With the new Commensal Open-Source Multimode Interferometer Cluster (COSMIC) digital backend on the Karl G. Jansky Very Large Array (VLA), we aim to conduct a search for technosignatures that is significantly larger, more sensitive, and more efficient than previously attempted. The COSMIC system is currently operational on the VLA, recording data, and designed with the flexibility to provide user-requested modes. This paper describes the hardware system design, the current software pipeline, and plans for future development.

This paper presents the identification of a new transient ULX candidate (ULX-3) with reaching a peak luminosity of ~ 4e39 erg/s in NGC 4254 by using archival Chandra, Swift X-Ray Telescope (Swift/XRT), Hubble Space Telescope (HST), and James Webb Space Telescope (JWST) observations. From precise astrometric calculations, unique optical, near-infrared (NIR) and mid-infrared (mid-IR) counterparts were found. The spectral energy distribution (SED) and color-magnitude diagrams (CMDs) of counterparts of the new ULX candidate were plotted to constrain the nature of the possible donor star. Evidence of a circumbinary disk was found from its SED with two blackbody temperatures of 1000 and 200 K. Moreover, according to the X-ray hardness ratios, ULX-3 exhibits very hard to very soft transitions as seen in some high-mass X-ray binaries (HMXBs) with Be-star donors Moreover, ULX-3 varies by more than two orders of magnitude in the 0.3-10 keV energy band as seen in typical transient ULXs.

The intracluster medium of the Perseus Cluster exhibits spiral-shaped X-ray surface brightness discontinuities known as ``cold fronts'', which simulations indicate are caused by the sloshing motion of the gas after the passage of a subcluster. Recent observations of Perseus have shown that these fronts extend to large radii. In this work, we present simulations of the formation of sloshing cold fronts in Perseus using the AREPO magnetohydrodynamics code, to produce a plausible scenario for the formation of the large front at a radius of 700 kpc. Our simulations explore a range of subcluster masses and impact parameters. We find that low-mass subclusters cannot generate a cold front that can propagate to such a large radius, and that small impact parameters create too much turbulence which leads to the disruption of the cold front before it reaches such a large distance. Subclusters which make only one core passage produce a stable initial front that expands to large radii, but without a second core passage of the subcluster other fronts are not created at a later time in the core region. We find a small range of simulations with subclusters with mass ratios of $R \sim 1:5$ and initial impact parameter of $\theta~\sim~20-25^\circ$ which not only produce the large cold front but a second set in the core region at later times. These simulations indicate that the ``ancient' cold front is $\sim$6-8.5 Gyr old. For the simulations providing the closest match with observations, the subcluster has completely merged into the main cluster.

We present continuum observations from the Atacama Large Millimeter/submillimeter Array (ALMA) of 10 high-redshift ($2.2 \le z \le 2.7$) ultraluminous quasars (QSOs) and constrain the presence of hot, ionized, circum-galactic gas in a stacking analysis. We measure a Compton-y parameter profile with a peak value of $(1.7 \pm 1.1) \times 10^{-6}$ at a radius of $\sim50$ kpc. We compare our stacked observations to active galactic nucleus (AGN) feedback wind models and generalized Navarro-Frenk-White (gNFW) pressure profile models to constrain the wind luminosity and halo mass of the stacked QSOs. Our observations constrain the observed stack's halo mass to $<1\times 10^{13}M_{\odot}$ and the stack's feedback wind power $<1\times 10^{12}L_{\odot}$, which is $<1$% of the bolometric luminosity of the quasar.

Can information theory provide insights into whether exoplanets are habitable? Here we apply information theory to a range of simulated exoplanet transmission spectra as a diagnostic tool to search for potential signatures of life on Earth-analog planets. We test the algorithms on three epochs of evolution for Earth-like planets orbiting a range of host stars. The James Webb Space Telescope and upcoming ground- and space-based missions promise to achieve sufficient high-resolution data that information theory can be applied to assess habitability. This approach provides a framework and a tool for observers to assess whether an exoplanet shows signs of habitability.

Phosphorus (P) is one of the important elements for the formation of life and plays a crucial role in several biochemical processes. Recent spectral line surveys have confirmed the existence of P-bearing molecules, especially PN and PO, in the star-formation regions, but their formation mechanisms are poorly understood. The P-bearing molecule phosphorus nitride (PN) is detected in several star-forming regions, but this molecule has been poorly studied at high gas densities ($\geq$10$^{6}$ cm$^{-3}$) hot molecular cores. In this article, we present the detection of the rotational emission line of PN with transition J = 3$-$2 towards the hot molecular cores G10.47+0.03 and G31.41+0.31, using the Atacama Compact Array (ACA). The estimated column densities of PN for G10.47+0.03 and G31.41+0.31 using the local thermodynamic equilibrium (LTE) model are (3.60$\pm$0.2)$\times$10$^{13}$ cm$^{-2}$ and (9.10$\pm$0.1)$\times$10$^{12}$ cm$^{-2}$ with an excitation temperature of 150$\pm$25 K. The fractional abundance of PN relative to H$_{2}$ is 2.76$\times$10$^{-10}$ for G10.47+0.03 and 5.68$\times$10$^{-11}$ for G31.41+0.031. We compute the two-phase warm-up chemical model of PN to understand the chemical evolution in the environment of hot molecular cores. After chemical modelling, we claim that PN is created in the gas phase via the neutral-neutral reaction between PO and N in the warm-up stage. Similarly, PN is destroyed via the ion-neutral reaction between H$_{3}$O$^{+}$ and PN.

Magnetic activity is a ubiquitous feature of stars with convective outer layers, with implications from stellar evolution to planetary atmospheres. Investigating the mechanisms responsible for the observed stellar activity signals from days to billions of years is important in deepening our understanding of the spatial configurations and temporal patterns of stellar dynamos, including that of the Sun. In this paper, we focus on three problems and their possible solutions. We start with direct field measurements and show how they probe the dependence of magnetic flux and its density on stellar properties and activity indicators. Next, we review the current state-of-the-art in physics-based models of photospheric activity patterns and their variation from rotational to activity-cycle timescales. We then outline the current state of understanding in the long-term evolution of stellar dynamos, first by using chromospheric and coronal activity diagnostics, then with model-based implications on magnetic braking, which is the key mechanism by which stars spin down and become inactive as they age. We conclude by discussing possible directions to improve the modeling and analysis of stellar magnetic fields.

Y dwarfs, the coolest known spectral class of brown dwarfs, overlap in mass and temperature with giant exoplanets, providing unique laboratories for studying low-temperature atmospheres. However, only a fraction of Y dwarf candidates have been spectroscopically confirmed. We present Keck/NIRES near-infrared spectroscopy of the nearby ($d \approx 6-8$ pc) brown dwarf CWISE J105512.11+544328.3. Although its near-infrared spectrum aligns best with the Y0 standard in the $J$-band, no standard matches well across the full $YJHK$ wavelength range. The CWISE J105512.11+544328.3 NH$_3$-$H$ = 0.427 $\pm$ 0.0012 and CH$_4$-$J$ = 0.0385 $\pm$ 0.0007 absorption indices and absolute Spitzer [4.5] magnitude of 15.18 $\pm$ 0.22 are also indicative of an early Y dwarf rather than a late T dwarf. CWISE J105512.11+544328.3 additionally exhibits the bluest Spitzer [3.6]$-$[4.5] color among all spectroscopically confirmed Y dwarfs. Despite this anomalously blue Spitzer color given its low luminosity, CWISE J105512.11+544328.3 does not show other clear kinematic or spectral indications of low metallicity. Atmospheric model comparisons yield a log(g) $\le$ 4.5 and $T_{\rm eff} \approx 500 \pm 150$ K for this source. We classify CWISE J105512.11+544328.3 as a Y0 (pec) dwarf, adding to the remarkable diversity of the Y-type population. JWST spectroscopy would be crucial to understanding the origin of this Y dwarf's unusual preference for low-gravity models and blue 3-5 $\mu$m color.

This work shows that a class of astrodynamics problems subject to mission constraints can be efficiently solved using the Theory of Functional Connections (TFC) mathematical framework by a specific change of coordinates. In these problems, the constraints are initially written in non-linear and coupled mathematical forms using classical rectangular coordinates. The symmetries of the constrained problem are used to select a new system of coordinates that transforms the non-linear constraints into linear. This change of coordinates is also used to isolate the components of the constraints. This way the TFC technique can be used to solve the ordinary differential equations governing orbit transfer problems subject to mission constraints. Specifically, this paper shows how to apply the change of coordinates method to the perturbed Hohmann-type and the one-tangent burn transfer problems.

In this study, we employ a convolutional neural network to classify gravitational waves originating from core-collapse supernovae. Training was conducted using spectrograms derived from three-dimensional numerical simulations of waveforms, which were injected onto real noise data from the third observing run of both Advanced LIGO and Advanced Virgo. To gain insights into the model's decision-making process, we apply class activation mapping techniques to visualize the regions in the input image that are significant for the model's prediction. Our model distinguished between 9 different waveforms and noise with an accuracy of 98.4% at 1 kpc. Visualization through class activation mapping revealed that the model's predictions predominantly rely on specific features within the input spectrograms, namely the g-mode and low-frequency modes.

Atmospheric turbulence profile plays an important role in designing and operating adaptive optics (AO) systems with multiple laser guide stars. To obtain representative free atmospheric profiles and resolved ground layer profiles for future AO systems at the Subaru telescope, we are conducting the SHARPEST (Shack-Hartmann Atmospheric tuRbulence Profiling Experiment at the Subaru Telescope) project. In this project, we develop a turbulence profiler comprising two Shack-Hartmann (SH) sensors to observe a pair of bright stars through the Subaru telescope with high spatial sampling by 2 cm subapertures. We perform two analyses on the SH spot data: variance analysis on the spot scintillation for free atmospheric profiles, and on the spot slope for ground layer profiles. This paper introduces the initial results of free atmospheric profiles as well as total seeing values and wind profiles obtained by the first two engineering runs. The free atmospheric profiles reconstructed by the two independent SH sensors show good agreement. The results are also consistent with simultaneous measurements by another profiler except for turbulence strength at ~1 km, which is explained by an overestimation problem of scintillation-based profilers. Measured total seeing values are also smaller than the simultaneous measurements, possibly due to the difference in ground layer turbulence between the two sites. The wind profiles show good consistency with the direct measurements by a rawinsonde. Through this study, we establish a method to constrain the free atmospheric profile, the total seeing, and the wind profile by analysing data from a single SH sensor with fine subapertures.

We consider the possibility that aeolian (wind blown) processes occur on small, 10 to 100~km diameter, planetesimals when they were embedded in the protosolar nebula. Drag from a headwind within a protostellar disk is sufficiently large to loft cm and smaller sized particles off the surface of a 10 km diameter asteroid in the inner solar system (at a few AU), and micron sized particles off the surface of a 10 km diameter object in the Transneptunian region. The headwind is sufficiently strong to overcome surface cohesion in the inner solar system, but not in the outer solar system. However, in the outer solar system, surface particles can be redistributed or escape due to impacts from particles that are in the protosolar disk's wind. Based on scaling crater ejecta, we estimate that impacts from particles in the headwind will lead to erosion of mass rather than accretion for planetesimals below about 10 km in diameter. The erosion limit is independent of material strength but proportional to the wind velocity. We explore the sensitivity of splash particle trajectories to particle size, headwind velocity and Reynolds number. Winds from a protostellar disk could account for Kuiper Belt Object (486958) Arrokoth's smooth undulating terrain but only during an epoch of high particle flux and low wind velocity. These conditions could have been present during and just after coalescence of Arrokoth's building blocks.

Some fraction of the diffuse photon background is supposed to be linked to high-energy neutrinos by astrophysical mechanisms of production and electromagnetic cascades. This article presents a simulation study of axion-like particles (ALPs) implications for that component, exploiting transport equations. Alternations of that spectrum due to ALP-photon conversion in the intergalactic magnetic field (IGMF) in the cases of various ALP parameters and mixing regimes at sources are studied. The results indicate considerable influence of IGMF-conversion on the ALP-photon flux even in the case of inverse ALP-photon coupling constant $M$ equal to $10^{11.5}$ GeV and some residual effects in the case of $M=10^{12}$ GeV. Furthermore, the scenario shows another aspect of a complex multimessenger interplay between IceCube and Fermi data, to a certain extent relieving the tension between them.

The Atacama Large Millimeter/submillimeter Array (ALMA) was used in 2021 to image the carbon-rich evolved star R Lep in Bands 8-10 (397-908 GHz) with baselines up to 16 km. The goal was to validate the calibration, using band-to-band (B2B) phase referencing with a close phase calibrator J0504-1512, 1.2 deg from R Lep in this case, and the imaging procedures required to obtain the maximum angular resolution achievable with ALMA. Images of the continuum emission and the HCN maser line at 890.8 GHz, from the J=10-9 transition between the (1110) and (0400) vibrationally excited states, achieved angular resolutions of 13, 6, and 5 mas in Bands 8-10, respectively. Self-calibration (self-cal) was used to produce ideal images as to compare with the B2B phase referencing technique. The continuum emission was resolved in Bands 9 and 10, leaving too little flux for self-cal of the longest baselines, so these comparisons are made at coarser resolution. Comparisons showed that B2B phase referencing provided phase corrections sufficient to recover 92%, 83%, and 77% of the ideal image continuum flux densities. The HCN maser was sufficiently compact to obtain self-cal solutions in Band 10 for all baselines (up to 16 km). In Band 10, B2B phase referencing as compared to the ideal images recovered 61% and 70% of the flux density for the HCN maser and continuum, respectively.

Novae are luminous explosions in close binaries which host a white dwarf and a companion donor star. They are triggered by a thermonuclear runaway when the white dwarf accretes a critical amount of matter from the secondary. Though novae are established as high-energy gamma-ray emitters through observations by the Fermi Large Area Telescope (LAT), the origin of the gamma-ray emission, whether it is hadronic or leptonic, had been under intense debate until very recently. RS Ophiuchi (RS Oph) is a well-known recurrent symbiotic nova with a recurrence time scale of 15 years. The most recent outburst of RS Oph in 2021 brought the first detection of very-high-energy (VHE) gamma rays from a nova ever. The first Large-Sized Telescope prototype (LST-1) of the Cherenkov Telescope Array observed this historic event along with H.E.S.S. and MAGIC. The LST-1 observations in the first days after the burst onset show a clear VHE gamma-ray signal from RS Oph. The low energy threshold of LST-1 allows us to reconstruct the RS Oph gamma-ray spectrum down to $\sim$30 GeV, providing the best connection of the VHE gamma-ray data to the Fermi LAT energy range. The results from the analysis of the LST-1 observations are consistent with those obtained with H.E.S.S. and MAGIC, and also support a hadronic origin for the observed gamma-ray fluxes. In this contribution, we will present the analysis results of the LST-1 observations of the 2021 outburst of RS Oph.

Many massive quiescent galaxies have been discovered at $z>2$ thanks to multi-wavelength deep and wide surveys, however, substantial deep near-infrared spectroscopic observations are needed to constrain their star-formation histories statistically. Here, we present a new technique to select quiescent galaxies with a short quenching timescale ($\leq0.1$ Gyr) at $z\sim2$ photometrically. We focus on a spectral break at $\sim1600$ \AA~that appears for such fast-quenching galaxies $\sim1$ Gyr after quenching when early A-type stars go out, but late A-type stars still live. This spectral break at $z\sim2$ is similar to a Lyman break at $z\sim4$. We construct a set of color criteria for $z\sim2$ fast-quenching galaxies on $g-r$ vs. $r-i$ and $i-J$ vs. $J-H$ or $\rm i-[3.6]$ vs. $\rm [3.6]-[4.5]$ color diagrams, which are available with the existing and/or future wide imaging surveys, by simulating various model galaxy spectra and test their robustnesses using the COSMOS2020 catalog. Galaxies with photometric and/or spectroscopic redshifts $z\sim2$ and low specific star formation rates are successfully selected using these colors. The number density of these fast-quenching galaxy candidates at $z\sim2$ suggests that massive galaxies not so far above the star-formation main sequence at $z=3-4$ should be their progenitors.

As part of the GLOSTAR (GLObal view of STAR formation in the Milky Way) survey, we present the high-resolution continuum source catalog for the regions (l = 2-28, 36-40, 56-60, &|b|<1.0), observed with the Karl G. Jansky Very Large Array (VLA) in its B-configuration. The continuum images are optimized to detect compact sources on angular scales up to 4", and have a typical noise level of 1sigma ~ 0.08mJy/beam for an angular resolution of 1", which makes GLOSTAR currently the highest resolution as well as the most sensitive radio survey of the northern Galactic plane at 4-8GHz. We extracted 13354 sources above a threshold of 5sigma and 5437 sources above 7sigma that represent the high-reliability catalog. We determined the in-band spectral index (alpha) for the sources in the 7sigma-threshold catalog. The mean value is alpha=-0.6, which indicates that the catalog is dominated by sources emitting non-thermal radio emission. We identified the most common source types detected in radio surveys: 251 HII region candidates (113 new), 282 planetary nebulae (PNe) candidates (127 new), 784 radio star candidates (581 new), and 4080 extragalactic radio source candidates (2175 new). A significant fraction of HII regions and PNe candidates have alpha<-0.1 indicating that these candidates could contain radio jets, winds or outflows from high-mass and low-mass stellar objects. We identified 245 variable radio sources by comparing the flux densities of compact sources from the GLOSTAR survey and the Co-Ordinated Radio `N' Infrared Survey for High-mass star formation (CORNISH), and find that most of them are infrared quiet. The catalog is typically 95% complete for point sources at a flux density of 0.6 mJy (i.e. typical 7sigma level) and the systematic positional uncertainty is <= 0.1". The GLOSTAR data and catalogs are available online at https://glostar.mpifr-bonn.mpg.de.

Various particle accelerators operate in the space plasmas, filling the Galaxy with high energy particles, primary cosmic rays. Reaching the atmosphere of the earth, these particles originate extensive air showers consisting of millions of elementary particles, secondary cosmic rays, covering several large areas on the ground. During thunderstorms, strong electric fields modulate the energy spectra of secondary particles and, consequently, originate short and long particle bursts. Impulse amplifications of particle fluxes, the thunderstorm ground enhancements, TGEs manifest themselves as peaks in the time series of count rates of particle detectors, which coincide with thunderstorms, during which free electrons are accelerating and multiplied, forming electron gamma ray avalanches. Thus, electron accelerators emerging in thunderous atmospheres can significantly alter the frequency of EAS triggers.

The first science image released by the JWST reveals numerous galaxies in the distant background of the galaxy cluster SMACS J0723.3-7327. Some have claimed redshifts of up to $z \simeq 20$, challenging standard cosmological models for structure formation. Here, we present a lens model for SMACS J0723.3-7327 anchored on five spectroscopically-confirmed systems at $1.38 \leq z \leq 2.21$ that are multiply lensed, along with twelve other systems with proposed image counterparts sharing common colours, spectral energy distributions, and morphological features, but having unknown redshifts. Constrained only by their image positions and, where available, redshifts, our lens model correctly reproduces the positions and correctly predicts the morphologies and relative brightnesses of all these image counterparts, as well as providing geometrically-determined redshifts spanning $1.4 \lesssim z \lesssim 6.7$ for the twelve candidate multiply-lensed galaxies lacking spectroscopic measurements. From this lens model, we create a lens finder map that defines regions over which galaxies beyond a certain redshift are predicted to be multiply lensed. Applying this map to three galaxies claimed to be at $10 \lesssim z \lesssim 20$, we find no image counterparts at locations (with an uncertainty of $\sim$$0.^{\prime\prime}5$) where they ought to be sufficiently magnified to be detectable - suggesting instead that these galaxies lie at $z \lesssim 1.7-3.2$. In lieu of spectroscopy, the creation of reliable lens finder maps for cluster fields are urgently needed to test and constrain redshifts inferred from photometry for a rapidly increasing number of candidate high-$z$ galaxies found with the JWST.

In this paper, we investigate the momentum coupling between early dark energy (EDE) and cold dark matter to alleviate cosmological tensions. EDE has exhibited promising efficacy in addressing the Hubble tension, but it exacerbates the large-scale structure tension. We consider the interaction between EDE and cold dark matter, introducing a pure momentum exchange between them to alleviate the large-scale structure tension introduced by the EDE model. We find that this coupling model is consistent with the EDE model, yielding a higher value for $H_0$, which can resolve the Hubble tension. Additionally, the new model exhibits a suppressive effect on structure growth, contributing to the alleviation of the large-scale structure tension. By utilising the Markov Chain Monte Carlo method and incorporating various cosmological data, the coupling model constrains the best-fit values for $H_0$ to be $72.23$ km/s/Mpc and for $S_8$ to be 0.8192. Compared to the $\Lambda$CDM model, the new models have not fully resolved the large-scale structure tension. However, in contrast to the best-fit value of 0.8316 for $S_8$ obtained from the EDE model, the new model alleviates the negative impact of the EDE model.

TESS observations of pulsating hot main sequence stars paint a very different picture from what is currently accepted. There are large numbers of delta Scuti (DSCT) stars hotter than the theoretical hot edge of the instability strip, continuing to what appear to be DSCT stars of mid-B type (historically known as Maia variables). The frequencies of maximum amplitude in DSCT stars are in poor agreement with unstable frequencies from the models. There is a well-defined upper envelope in the frequencies of maximum amplitude as a function of effective temperature for DSCT and MAIA stars which requires an explanation. The gamma Doradus (GDOR) stars should be regarded as DSCT stars with suppressed high frequencies rather than a separate class. They are found mostly among the cool DSCT stars, but occur throughout the DSCT instability strip and as early A-type stars, where they merge with the SPB variables. The mixture of DSCT and GDOR stars throughout the instability strip is one example of the unexplained large variation of frequency patterns in DSCT stars. The location of beta Cephei stars in the H-R diagram agrees quite well with the models, but the observed frequencies are generally higher than predicted. There are no discernible boundaries between the traditional classes of pulsating stars. Existing pulsation models do not describe the observations at all well and a re-evaluation is required.

After 28 years from the conception of the pyramid WFS several new kind of devices able to convert wavefront shape into some sort of different illumination on a detector have been conceived. While, suspending momentarily any kind of modesty, I claim credit for being among the few that contributed to show at the time that there could be much more than just a lenslet array I continued -- with alternating successes -- to conceive other types of such devices.

After the second data release of Gaia, the number of new globular cluster candidates has increased importantly. However, most of them need to be properly characterised, both spectroscopically and photometrically, by means of radial velocities, metallicities, and deeper photometric observations. Our goal is to provide an independent confirmation of the cluster nature of Gran 4, a recently discovered globular cluster, with follow-up spectroscopic observations. The derived radial velocity for individual stars, coupled with proper motions, allows us to isolate cluster members from field stars, while the analysis of their spectra allows us to derive metallicities. By including in the analysis the recently confirmed clusters Gran 1, 2, 3, and 5, we aim to completely characterise the sample presented in Gran et al. 2022. Using Gaia DR3 and VVV catalogue data and MUSE@VLT observations, we performed a selection of cluster members based on their proper motions, radial velocities and their position in colour-magnitude diagrams. Furthermore, full spectral synthesis was performed on the cluster members, extracting surface parameters and metallicity from MUSE spectra. Finally, a completeness estimation was performed on the total globular cluster population of the Milky Way. We confirm the nature of Gran 4, a newly discovered globular cluster behind the Galactic bulge, with a mean radial velocity of ${\rm RV} = -265.28 \pm 3.92$ km s$^{-1}$ and a mean metallicity of ${\rm [Fe/H] = -1.72 \pm 0.32}$ dex. Additionally, independent measurements of the metallicities were derived for Gran 1, 2, 3, and 5. We also revise the observational lower mass limit for a globular cluster to survive in the bulge/disk environment. We estimate that $\sim 12-26$ globular clusters have still to be discovered on the other side of the Galaxy (i.e., behind the bulge/bar/disk), up to 20 kpc.

The JWST NIRSpec integral field unit (IFU) presents a unique opportunity to observe directly imaged exoplanets from 3-5um at moderate spectral resolution (R~2,700) and thereby better constrain the composition, disequilibrium chemistry, and cloud properties of their atmospheres. In this work, we present the first NIRSpec IFU high-contrast observations of a substellar companion that requires starlight suppression techniques. We develop specific data reduction strategies to study faint companions around bright stars, and assess the performance of NIRSpec at high contrast. First, we demonstrate an approach to forward model the companion signal and the starlight directly in the detector images, which mitigates the effects of NIRSpec's spatial undersampling. We demonstrate a sensitivity to planets that are 2e-6 fainter than their stars at 1'', or 2e-5 at 0.3''. Then, we implement a reference star point spread function (PSF) subtraction and a spectral extraction that does not require spatially and spectrally regularly sampled spectral cubes. This allows us to extract a moderate resolution (R~2,700) spectrum of the faint T-dwarf companion HD~19467~B from 2.9-5.2um with signal-to-noise ratio (S/N)~10 per resolution element. Across this wavelength range, HD~19467~B has a flux ratio varying between 1e-5-1e-4 and a separation relative to its star of 1.6''. A companion paper by Hoch et al. more deeply analyzes the atmospheric properties of this companion based on the extracted spectrum. Using the methods developed here, NIRSpec's sensitivity may enable direct detection and spectral characterization of relatively old (~1Gyr), cool (~250K), and closely separated (~3-5au) exoplanets that are less massive than Jupiter.

One of the hottest questions in the cosmology of self-interacting dark matter (SIDM) is whether scatterings can induce detectable core-collapse in halos by the present day. Because gravitational tides can accelerate core-collapse, the most promising targets to observe core-collapse are satellite galaxies and subhalo systems. However, simulating small subhalos is computationally intensive, especially when subhalos start to core-collapse. In this work, we present a hierarchical framework for simulating a population of SIDM subhalos, which reduces the computation time to linear order in the total number of subhalos. With this method, we simulate substructure lensing systems with multiple velocity-dependent SIDM models, and show how subhalo evolution depends on the SIDM model, subhalo mass and orbits. We find that an SIDM cross section of $\gtrsim 200$ cm$^2$/g at velocity scales relevant for subhalos' internal heat transfer is needed for a significant fraction of subhalos to core-collapse in a typical lens system at redshift $z=0.5$, and that core-collapse has unique observable features in lensing. We show quantitatively that core-collapse in subhalos is typically accelerated compared to field halos, except when the SIDM cross section is non-negligible ($\gtrsim \mathcal{O}(1)$ cm$^2$/g) at subhalos' orbital velocities, in which case evaporation by the host can delay core-collapse. This suggests that substructure lensing can be used to probe velocity-dependent SIDM models, especially if line-of-sight structures (field halos) can be distinguished from lens-plane subhalos. Intriguingly, we find that core-collapse in subhalos can explain the recently reported ultra-steep density profiles of substructures found by lensing with the \emph{Hubble Space Telescope}

We study mass ejection from a binary neutron star merger producing a long-lived massive neutron star remnant with general-relativistic neutrino-radiation hydrodynamics simulations. In addition to outflows generated by shocks and tidal torques during and shortly after the merger, we observe the appearance of a wind driven by spiral density waves in the disk. This spiral-wave-driven outflow is predominantly located close to the disk orbital plane and have a broad distribution of electron fractions. At higher latitudes, a high electron-fraction wind is driven by neutrino radiation. The combined nucleosynthesis yields from all the ejecta components is in good agreement with Solar abundance measurements.

Context. The advent of multi-dimensional solar flare simulations has led to numerous investigations of coronal flows and new physical insights. These studies have not yet included detailed analysis of the lower atmospheric responses such as down-flowing chromospheric compressions and chromospheric evaporation processes. Aims. In this work, we present an analysis of multi-dimensional flare simulations, including analysis of chromospheric up-flows and down-flows that help to elucidate multi-dimensional effects. We also provide important groundwork for comparing 1D and multi-dimensional models, with the aim that future multi-dimensional simulations can include detailed field-aligned physical processes. Methods. A localized anomalous resistivity initiates magnetic reconnection, which drives the evolution of a standard solar flare model. We vary the background magnetic field strength, to produce four flare simulations that cover a large span of observationally reported solar flare strengths. Chromospheric energy fluxes, and energy density maps are used to analyse the transport of energy from the corona to the lower atmosphere, and the resultant evolution of the flare. Quantities traced along 1D field-lines allow for detailed comparison with 1D evaporation models. We highlight the similarities, stressing deficiencies from simplified physics along these 1D flux tubes, and crucial effects that enter by multi-dimensional effects.

In this paper, we present the discovery of a dual AGN in a $\sim$14:1 minor merger between the galaxy SDSS J125417.98+274004.6 and its unnamed dwarf satellite. We calculated stellar masses of the primary and secondary galaxy to be 3.8$\times$10$^{10}$M$_{\odot}$ and 2.7$\times$10$^{9}$M$_{\odot}$, respectively. We used archival Chandra X-ray observations to assess AGN presence. We found that both AGN have comparable luminosities of $\sim$2$\times$10$^{42}$ erg s$^{-1}$, with the secondary AGN being more likely to be the dominant one. The galaxies are in the early stages of the merger and are connected by a tidal bridge. Computational works suggest that the secondary AGN should experience a brief but intensive period of Eddington-limit approaching accretion during the early stages of the merger. During the merger, the secondary black hole can increase its mass by a factor of ten. SDSS J125417.98+274004.6 is the first known dual AGN in an early-stage minor merger with a comparably or more luminous secondary AGN. As such, it will be of great value for future studies of merger-triggered accretion and black hole growth.

CORSIKA 8 is a new framework for air shower simulations implemented in modern C++17, based on past experience with existing codes like CORSIKA 7. The flexible and modular structure of the project allows the development of independent modules that can produce a fully customizable air shower simulation. The radio module in particular is designed to treat the signal propagation and electric field calculation to each antenna in an autonomous and flexible way. It provides the possibility to simulate simultaneously the radio emission calculated with two independent time-domain formalisms, the "Endpoint formalism" as implemented in CoREAS and the "ZHS" algorithm as ported from ZHAireS. Future development for the simulation of radio emission from particle showers in complex scenarios, for example cross-media showers penetrating from air into ice, can build on the existing radio module, re-using the establishes interfaces. In this work, we will present the design and implementation of the radio module in CORSIKA 8, and show a direct comparison of radio emission from air showers simulated with CORSIKA 8, CORSIKA 7 and ZHAireS.

The North Polar Spur (NPS) is a giant structure that is clearly visible in both radio and X-ray all-sky maps. We analyzed broadband radio observations covering a range between 22 MHz and 70 GHz to systematically analyze the thermal/non-thermal emissions associated with the NPS. We demonstrate that the radio emission of the NPS comprises synchrotron, free-free, and dust emission; however, synchrotron emissions dominate over other emissions, especially at high galactic-latitudes. Moreover, the synchrotron spectra exhibit a power-law behavior with $N(\gamma)\propto\gamma^{-s}$ ($s\simeq1.8-2.4$) up to a few GHz moderated by a turnover at $\nu_{\rm brk} \simeq 1$ GHz, above which the spectral index $s$ decrease by one. Assuming that the turnover is due to the electrons cooled by synchrotron radiation before escaping (or advecting) from the emission region, the magnetic field strength can be estimated to be $B\sim 8 \rm\mu G$ if the NPS is a distant structure that is near the Galactic Center (GC). However, an unreasonably strong $B\sim 114\rm\mu G$ is required if the NPS is near the local supernova remnant (SNR). The corresponding non-thermal energy stored in the NPS is $E_{\rm n/th}\simeq 4.4\times 10^{55}$ erg in the GC scenario, whereas $E_{\rm n/th}\simeq 4.1\times 10^{52}$ erg is difficult to explain with a single local SNR. We also estimated the gamma-ray emission associated with the NPS through inverse Comptonization of the cosmic microwave background (CMB), which peaks at 100 - 1000 keV with a flux of $\nu F_{\nu}\sim 10^{-9}$ $\rm erg\,cm^{-2}s^{-1}sr^{-1}$ in the GC model, and may be a good candidate for detection by future X-ray/gamma-ray observatories.

The Hubble constant problem is that the values of Hubble constant from the observation of cosmic microwave background assuming the LambdaCDM model disagrees with the values from direct measurements. This problem suggests some new physics beyond the LambdaCDM model. Typically there are two ways of reconciliation: one is the realization of smaller value of sound horizon at recombination, and the other is the modification of the way of expansion of the universe after recombination. In this letter we examine the latter possibility by comparing two typical phenomenological dark energy models with the distance-redshift relation provided by Pantheon catalogue of supernova observations and galaxy surveys by BOSS and eBOSS collaborations. Though these phenomenological dark energy models globally fit observations better than the LambdaCDM model, they are strongly disfavored by the distance-redshift relation as almost the same level as the LambdaCDM model defined by cosmic microwave background observations. The distance-redshift relation strongly suggests some new physics which realize smaller value of sound horizon at recombination.

We present the host galaxies of four apparently non-repeating fast radio bursts (FRBs), FRBs 20181223C, 20190418A, 20191220A, and 20190425A, reported in the first Canadian Hydrogen Intensity Mapping Experiment (CHIME/FRB) catalog. Our selection of these FRBs is based on a planned hypothesis testing framework where we search all CHIME/FRB Catalog-1 events that have low extragalactic dispersion measure (< 100 pc cm$^{-3}$), with high Galactic latitude (|b| > 10$\deg$) and saved baseband data. We associate the selected FRBs to galaxies with moderate to high star-formation rates located at redshifts between 0.027 and 0.071. We also search for possible multi-messenger counterparts, including persistent compact radio and gravitational wave (GW) sources, and find none. Utilizing the four FRB hosts from this study along with the hosts of 14 published local Universe FRBs (z < 0.1) with robust host association, we conduct an FRB host demographics analysis. We find all 18 local Universe FRB hosts in our sample to be spirals (or late-type galaxies), including the host of FRB 20220509G, which was previously reported to be elliptical. Using this observation, we scrutinize proposed FRB source formation channels and argue that core-collapse supernovae are likely the dominant channel to form FRB progenitors. Moreover, we infer no significant difference in the host properties of repeating and apparently non-repeating FRBs in our local Universe FRB host sample. Finally, we find the burst rates of these four apparently non-repeating FRBs to be consistent with those of the sample of localized repeating FRBs observed by CHIME/FRB. Therefore, we encourage further monitoring of these FRBs with more sensitive radio telescopes.

TianQin (TQ) project plans to deploy three satellites in space around the Earth to measure the displacement change of test masses caused by gravitational waves via laser interferometry. The requirement of the acceleration noise of the test mass is on the order of $10^{-15}~\,{\rm m}\,{\rm s}^{-2}\,{\rm Hz}^{-1/2}$ in the sensitive frequency range of TQ, %the extremely precise acceleration measurement requirements make it necessary to investigate acceleration noise due to space magnetic fields. which is so stringent that the acceleration noise caused by the interaction of the space magnetic field with the test mass needs to be investigated. In this work, by using the Tsyganenko model, a data-based empirical space magnetic field model, we obtain the magnetic field distribution around TQ's orbit spanning two solar cycles in 23 years from 1998 to 2020. With the obtained space magnetic field, we derive the distribution and amplitude spectral densities (ASDs) of the acceleration noise of TQ in 23 years. Our results reveal that the average values of the ratio of the acceleration noise cauesd by the space magnetic field to the requirements of TQ at 1 mHz ($R_{\rm 1mHz}$) and 6 mHz ($R_{\rm 6mHz}$) are 0.123$\pm$0.052 and 0.027$\pm$0.013, respectively. The occurence probabilities of $R_{\rm 1mHz}>0.2$ and $>0.3$ are only 7.9\% and 1.2\%, respectively, and $R_{\rm 6mHz}$ never exceeds 0.2.

The population of optical transients evolving within a time-scale of a few hours or a day (so-called Fast Optical Transients, FOTs) has recently been debated extensively. In particular, our understanding of extragalactic FOTs and their rates is limited. We present a search for extragalactic FOTs with the Tomo-e Gozen high-cadence survey. Using the data taken from 2019 August to 2022 June, we obtain 113 FOT candidates. Through light curve analysis and cross-matching with other survey data, we find that most of these candidates are in fact supernovae, variable quasars, and Galactic dwarf novae, that were partially observed around their peak brightness. We find no promising candidate of extragalactic FOTs. From this non-detection, we obtain upper limits on the event rate of extragalactic FOTs as a function of their time-scale. For a very luminous event (absolute magnitude M<-26 mag), we obtain the upper limits of 4.4 x 10^{-9} Mpc^{-3} yr^{-1} for a time-scale of 4 h, and 7.4 x 10^{-10} Mpc^{-3} yr^{-1} for a time-scale of 1 d. Thanks to our wide (although shallow) surveying strategy, our data are less affected by the cosmological effects, and thus, give one of the more stringent limits to the event rate of intrinsically luminous transients with a time-scale of < 1 d.

Aerocapture uses atmospheric drag to decelerate spacecraft and achieve orbit insertion. One of the significant risks associated with aerocapture is the uncertainty in the atmospheric density, particularly for outer planets. The paper performs a comparative study of the atmospheric uncertainties and provides design rules for aerocapture missions. The atmospheres of Venus, Mars, and Titan are well-characterized for engineering purposes. At the altitude ranges relevant for aerocapture, the 3$\sigma$ density variation is approximately $\pm$30%, $\pm$50%, $\pm$30% for Venus, Mars, and Titan respectively. With no in-situ data, the atmospheres of Uranus and Neptune are not as well characterized as the other bodies. For both Uranus and Neptune, the GRAM suite provides a 3$\sigma$ density variation of approximately $\pm$30% for the relevant altitude ranges which is considered an optimistic estimate. Until in-situ data from an atmospheric probe becomes available, a more conservative global min-max estimate is recommended to accommodate the worst-case scenario. The study presents a graphical method for selection of the optimal entry flight path angle when considering the atmospheric uncertainties to ensure the on-board guidance is given the best possible initial state for targeting the desired exit state post aerocapture.

We report the detection of a $\gamma$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $\Gamma (E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $\gamma$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $\gamma$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

Ultralight axion (ULA) can be one of the potential candidates for dark matter. The extremely low mass of the ULA can lead to a de Broglie wavelength around galaxy-size, and result in a suppression of the structure growth at small scales. This allows us to study its properties via galaxy surveys for probing cosmic structure formation. In this work, we forecast the constraint on the ULA particle mass $m_{\text{a}}$ and relative fraction to dark matter $f_{\text{a}} = \Omega_{\text{a}}/\Omega_{\text{d}}$ for the forthcoming Stage IV space-based optical survey equipment $\it{CSST}$ (China Space Station Telescope). We focus on the $\it{CSST}$ cosmic shear and galaxy clustering photometric surveys, and forecast the measurements of shear, galaxy, and galaxy-galaxy lensing power spectra (i.e. 3$\times$2pt). The effects of neutrino, baryonic feedback, and uncertainties of intrinsic alignment, shear calibration, galaxy bias, and photometric redshift are also included in the analysis. After performing a joint constraint on all the cosmological and systematical parameters using Markov Chain Monte Carlo (MCMC) technique, we obtain a lower limit of the ULA particle mass $\text{log}_{10}(m_{\text{a}}/\text{eV}) \geqslant -22.5$ and an upper limit of the ULA fraction $f_{\text{a}} \leqslant 0.83$ at 95\% confidence level, and $\text{log}_{10}(m_{\text{a}}/\text{eV}) \geqslant -21.9$ with $f_{\text{a}} \leqslant 0.77$ when ignoring the baryonic feedback. We find that the CSST photometric surveys can improve the constraint on the ULA mass by about one order of magnitude, compared to the current constraints using the same kind of observational data.

Compilation of papers presented by the GAPS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan.

Context. Binary Cepheids play an important role in investigating the calibration of the classical Cepheid period-luminosity relationship. Therefore a thorough study of individual Cepheids belonging to binary systems is necessary. Aims. Our aim is to determine the orbit of the binary system V1344 Aql using newly observed and earlier published spectroscopic and photometric data. Methods. We collected new radial velocity observations using medium resolution (${R \approx 11000}$ and ${R \lessapprox 20000}$) spectrographs and we updated the pulsation period of the Cepheid based on available photometric observations using $O-C$ diagram. Separating the pulsational and orbital radial velocity variations for each observational season (year), we determined the orbital solution for the system using $\chi^2$ minimisation. Results. The updated pulsation period of the Cepheid estimated for the epoch of HJD 2458955.83 is 7.476826 days. We determined orbital elements for the first time in the literature. The orbital period of the system is about 34.6 years, with an eccentricity e = 0.22.

The vertical phase-space spiral (snail) is a direct sign of dis-equilibrium of Milky Way's disc. Nevertheless, the wrapping of the phase snail contains information of the vertical potential. We propose a novel method to measure the vertical potential utilizing the intersections between the snail and $z$/$V_{z}$ axes, for which we know the maximum vertical heights ($Z_{max}$) or the maximum vertical velocities ($V_{z,max}$). Using a refined linear interpolation method, we directly obtain $(Z_{max},\ \frac{1}{2}V_{z,max}^{2})$ for these snail intersections to constrain the vertical potential profile empirically. Our method is model independent since no assumptions about the snail shape or the vertical potential have been made. Although the snail binned by the guiding center radius ($R_{g}$) is more prominent, it traces a vertical potential shallower than that of the snail binned by the same Galactocentric radius ($R$). We apply an empirical method to correct this difference. We measure the snail intersections in several $R_{g}$ bins within $7.5< R_{g} < 11.0$ kpc for Gaia DR3, and apply the interpolation method to deduce the potential values at several vertical heights. The potential at the snail intersections, as well as the following mass modeling are consistent with the popular Milky Way potentials in the literature. For the $R_{g}$-binned phase snail in the Solar neighborhood, the mass modeling indicates a local dark matter density of $\rho_{\rm dm}= 0.0150\pm0.0031$ $\rm M_{\odot}\,pc^{-3}$, consistent with previous works. Our method could be applied to larger radial ranges in future works, to provide independent and stronger constraints on the Milky Way's potential.

We present a comprehensive study of the molecular gas properties of 17 Type 2 quasars at $z <$ 0.2 from the Quasar Feedback Survey (L$_{[OIII]}$ > $10^{42.1}$ $\rm ergs^{-1}$), selected by their high [OIII] luminosities and displaying a large diversity of radio jet properties, but dominated by LIRG-like galaxies. With these data, we are able to investigate the impact of AGN and AGN feedback mechanisms on the global molecular interstellar medium. Using APEX and ALMA ACA observations, we measure the total molecular gas content using the CO(1-0) emission and homogeneously sample the CO spectral line energy distributions (SLEDs), observing CO transitions (J$_{up}$ = 1, 2, 3, 6, 7). We observe high $r_{21}$ ratios (r$_{21}$ = L'$_{CO(2-1)}$/L'$_{CO(1-0)}$) with a median $r_{21}$ = 1.06, similar to local (U)LIRGs (with $r_{21}$ $\sim$ 1) and higher than normal star-forming galaxies (with r$_{21}$ $\sim$ 0.65). Despite the high $r_{21}$ values, for the 7 targets with the required data we find low excitation in CO(6-5) & CO(7-6) ($r_{61}$ and $r_{62}$ < 0.6 in all but one target), unlike high redshift quasars in the literature, which are far more luminous and show higher line ratios. The ionised gas traced by [OIII] exhibit systematically higher velocities than the molecular gas traced by CO. We conclude that any effects of quasar feedback (e.g. via outflows and radio jets) do not have a significant instantaneous impact on the global molecular gas content and excitation and we suggest that it only occurs on more localised scales.

Detailed analysis of kinematics of the Milky Way disk in the solar neighborhood based on the GAIA DR3 catalog reveals the existence of peculiarities in the stellar velocity distribution perpendicular to the galactic plane. We study the influence of resonances -- the outer Lindblad resonance and the outer vertical Lindblad resonance -- of a rotating bar with stellar oscillations perpendicular to the plane of the disk, and their role in shaping the spatial and the velocity distributions of stars. We find that the $Z$ and $V_Z$ distributions of stars with respect to $L_Z$ are affected by the outer Lindblad resonance. The existence of bar resonance with stellar oscillations perpendicular to the plane of the disk is demonstrated for a long (large semi-axis 5 kpc) and fast rotating bar with $\Omega_{b}= 60.0$ $km~s^{-1}~kpc^{-1}$. We show also that, in the model with the long and fast rotating bar, some stars in the 2:1 OLR region deviate far from their original places, entering the bar region. A combination of resonance excitation of stellar motions at the 2:1 OLR region together with strong interaction of the stars with the bar potential leads to the formation of the group of 'escapees', i.e., stars that deviate in $R$ and $Z$ -- directions at large distances from the resonance region. Simulations, however, do not demonstrate any noticeable effect on $V_Z$-distribution of stars in the solar neighborhood

We use multiband archival HST observations to measure the photometric and structural parameters of the M81 globular cluster that hosts the Fast Radio Burst FRB 20200120E. Our best-fitting King model has an effective radius $r_h = 3.06$ pc with a moderate King model concentration of $c = 53$, and an inferred core radius of 0.81 pc. We revisit the exact astrometric location of the FRB within the cluster, and find that FRB 20200120E is located 1.92 pc from the center, but within the projected half-light radius. We estimate the relative encounter rate of the FRB host, along with the corresponding rates of 210 other globular clusters in M81, and compare these values with the encounter rates of Galactic globular clusters. The FRB resides in a globular cluster with an encounter rate that is moderately higher than the median stellar encounter rate in our two comparison samples. While the estimated encounter rate of the FRB host cluster (e.g., $\sim50\%$ of a cluster like 47 Tuc) is sufficient to allow the possibility that the FRB formed dynamically, our results do not place strong constraints on this scenario due to the limitations of the available HST data and the possible systematic uncertainties and selection effects in the comparison data.

HD 139139 (a.k.a. 'The random transiter') is a star that exhibited enigmatic transit-like features with no apparent periodicity in K2 data. The shallow depth of the events ($\sim$200 ppm -- equivalent to transiting objects with radii of $\sim$1.5 R$_\oplus$ in front of a Sun-like star), and their non-periodicity, constitutes a challenge for the photometric follow-up of this star. The goal of this study is to confirm with independent measurements the presence of shallow, non-periodic transit-like features on this object. We performed observations with the CHEOPS satellite, for a total accumulated time of 12.75 d, distributed in visits of roughly 20 h, in two observing campaigns, in years 2021 and 2022. The precision of the data is sufficient to detect 150 ppm features with durations longer than 1.5 h. We use the duration and times of the events seen in the K2 curve to estimate how many should have been detected in our campaigns, under the assumption that their behaviour during the CHEOPS observations is the same as in the K2 data of 2017. We do not detected events with depths larger than 150 ppm in our data set. If the frequency, depth, and duration of the events were the same as in the K2 campaign, we estimate the probability of having missed all events due to our limited observing window as 4.8 %. We suggest three different scenarios: 1) our observing window was not long enough, and the events were missed with the estimated 4.8 % probability, 2) the events recorded in K2 observations were time-critical, and the mechanism producing them was either not active in the two 2021 and 2022 campaigns, or producing shallower events under our detectability level, or 3) the enigmatic events in K2 data are the result of an unidentified and infrequent instrumental noise in the original data set or its data treatment.

We explore the influence of the Milky Way galaxy's chemical evolution on the formation, structure, and habitability of the Solar system. Using a multi-zone Galactic Chemical Evolution (GCE) model, we successfully reproduce key observational constraints, including the age-metallicity ([Fe/H]) relation, metallicity distribution functions, abundance gradients, and [X/Fe] ratio trends for critical elements involved in planetary mineralogy, including C, O, Mg, and Si. Our GCE model suggests that the Sun formed in the inner Galactic disc, $R_{\rm birth,\odot}\approx 5$ kpc. We also combined a stoichiometric model with the GCE model to examine the temporal evolution and spatial distribution of planet building blocks (PBBs) within the Milky Way galaxy, revealing trends in the condensed mass fraction ($f_{\rm cond}$), iron-to-silicon mass fraction ($f_{\rm iron}$), and water mass fraction ($f_{\rm water}$) over time and towards the inner Galactic disc regions. Specifically, our model predicts a higher $f_{\rm cond}$ in the protoplanetary disc within the inner regions of the Milky Way galaxy, as well as an increased $f_{\rm iron}$ and a decreased $f_{\rm water}$ in the inner regions. Based on these findings, we discuss the potential impact of the Sun's birth location on the overall structure and habitability of the Solar System.

The X-ray source CXOU J163802.6-471358is thought to be a pulsar wind nebula (PWN), as it shows an extended, $\approx 40$ arcsec trail from a compact source. Here we present GMRT observations of this source at 330 and 1390 MHz, which reveal a remarkable linear radio trail $\approx 90$ arcsec in extent. Although the radio trail points back to the supernova remnant (SNR) G338.1+0.4, $\approx 50$ arcmin from ,CXOU J163802.6-471358 associating it with this remnant would require a very large velocity for the pulsar. There are no known galactic SNRs close to the PWN and radio trail. No pulsar has yet been identified in CXOU J163802.6-471358, but if one could be found, this would allow more quantitative studies of the PWN and radio trail to be made.

Assuming a slow-roll inflationary model where conformal invariance of the Maxwell action is broken via a non-minimal kinetic coupling term, we investigate the non-Gaussian three-point cross-correlation function between the primordial curvature perturbation and the primordial magnetic field, under a fairly general choice of initial vacua for both the scalar and the gauge field sectors. Among the possible triangular configurations of the resulting cross-bispectrum, we find that the squeezed limit leads to local-type non-Gaussianity allowing a product form decomposition in terms of the scalar and magnetic power spectra, which is a generic result independent of any specific choice of the initial states. We subsequently explore its detection prospects in the CMB via correlations between pre-recombination $\mu$-type spectral distortions and temperature anisotropies, sourced by such a primordial cross-correlation. Our analysis with several proposed next-generation CMB missions forecasts a low value of the signal-to-noise ratio (SNR) for the $\mu T$ spectrum if both the vacua are assumed to be pure Bunch-Davies. On the contrary, the SNR may be enhanced significantly for non-Bunch-Davies initial states for the magnetic sector within allowed bounds from current CMB data.

To understand the chemistry of sulphur (S) in the interstellar medium, models need to be tested by observations of S-bearing molecules in different physical conditions. We analyse observations obtained with the IRAM 30m telescope towards 15 well-known cores classified in the three main evolutionary stages of the high-mass star-formation process: high-mass starless cores,high-mass protostellar objects, and ultracompact HII regions. We detected rotational lines of SO, SO+, NS, C34S, 13CS, SO2, CCS, H2S, HCS+, OCS, H2CS, and CCCS. We also analyse for the first time lines of the NO molecule to complement the analysis. From a local thermodynamic equilibrium approach, we have derived column densities of each species and excitation temperatures. Based on a statistical analysis on the line widths and the excitation temperatures, we find that: NS, C34S, 13CS, CCS, and HCS+ trace cold, quiescent, and likely extended material; OCS, and SO2 trace warmer, more turbulent, and likely denser and more compact material; SO and perhaps SO+ trace both quiescent and turbulent material depending on the target. The abundances of SO, SO2, and H2S show the strongest positive correlations with the kinetic temperature, believed to be an evolutionary indicator. Moreover, the sum of all molecular abundances show an enhancement of gaseous S from the less evolved to the more evolved stages. These trends could be due to the increasing amount of S sputtered from dust grains owing to the increasing protostellar activity with evolution. The average abundances in each evolutionary group increase especially in the oxygen-bearing molecules, perhaps due to the increasing abundance of atomic oxygen with evolution owing to photodissociation of water in gas phase. Our observational work represents a test-bed for theoretical studies aimed at modelling the chemistry of sulphur during the evolution of high-mass star-forming cores

While looking for recently discovered distant Cepheids with an interesting chemical composition, we noticed one star (OGLE GD-CEP-1353) with extremely large equivalent widths of spectral lines of heavy elements. The aim of this work is to perform an abundance analysis, and to find a possible explanation for the found chemical anomaly. Quantitative analysis of the equivalent widths and synthetic spectrum synthesis were used to derive abundances in this star. Both local and nonlocal thermodynamic equilibrium (LTE and NLTE) approximations were used in our analysis. Abundances of 28 chemical elements from carbon to thorium were derived. While light and iron peak elements show abundances typical for distant Cepheids (located in the outer disk), the s-process elements are overabundant about one dex. r-process elements are slightly less overabundant. This makes the star a unique Cepheid of our Galaxy.

Continuing from our prior work (Alhassan et al. 2022), where a single detector data of the Einstein Telescope (ET) was evaluated for the detection of binary black hole (BBHs) using deep learning (DL). In this work we explored the detection efficiency of BBHs using data combined from all the three proposed detectors of ET, with five different lower frequency cutoff (Flow): 5 Hz, 10 Hz, 15 Hz, 20 Hz and 30 Hz, and the same previously used SNR ranges of: 4-5, 5-6, 6-7, 7-8 and >8. Using ResNet model (which had the best overall performance on single detector data), the detection accuracy has improved from 60%, 60.5%, 84.5%, 94.5% and 98.5% to 78.5%, 84%, 99.5%, 100% and 100% for sources with SNR of 4-5, 5-6, 6-7, 7-8 and >8 respectively. The results show a great improvement in accuracy for lower SNR ranges: 4-5, 5-6 and 6-7 by 18.5%, 24.5%, 13% respectively, and by 5.5% and 1.5% for higher SNR ranges: 7-8 and >8 respectively. In a qualitative evaluation, ResNet model was able to detect sources at 86.601 Gpc, with 3.9 averaged SNR (averaged SNR from the three detectors) and 13.632 chirp mass at 5 Hz. It was also shown that the use of the three detectors combined data is appropriate for near-real-time detection, and can be significantly improved using more powerful setup.

The indirect search for dark matter is typically restricted to individual photon bands and instruments. In the context of multiwavelength observations, finding a weak signal in large fore- and backgrounds at only one wavelength band is hampered by systematic uncertainties dominating the signal strength. Dark matter particle annihilation is producing Standard Model particles of which the prompt photon emission is searched for in many studies. However, also the secondary emission of charged particles from dark matter annihilation in the TeV range results in comparable or even stronger fluxes in the GHz-GeV range. In this study, we calculate the prompt and secondary emission of a scotogenic WIMP with a mass of $1\,\mathrm{TeV}$ in 27 dwarf galaxies of the Milky Way. For the secondary emission, we include Inverse Compton scattering, bremsstrahlung, and synchrotron radiation, which results in a "triple hump" structure characteristic for only dark matter and no other astrophysical source. In order to determine the best candidates for multi-instrument analyses, we estimate the diffuse emission component of the Milky Way itself, including its own dark matter halo from the same scotogenic WIMP model. We find signal-to-background ratios of individual sources on the order of $10^{-3}$-$10^{-2}$ across X- to $\gamma$-rays assuming $J$-factors for the cold dark matter distribution inferred from observations and no additional boosting due to small-scale clumping. We argue that a joint multi-wavelength analysis of all nearby galaxies as well as the extension towards the Cosmic Gamma-ray Background is required to disentangle possible dark matter signals from astrophysical back- and foregrounds.

We present the results obtained from the spectral studies of black hole X-ray binary GX~339--4 using \astrosat~ observations during its 2021 outburst. \astrosat~ observed the source in the intermediate state for $\sim600$ ks. The combined spectra of SXT and LAXPC in the $0.7-25$ keV energy range are studied with phenomenological and physical models. The spectral study reveals a receding disc and a contracting corona during the observation period. The outflow rate is found to be increased though the accretion rates did not vary during the observation period. The X-ray flux decreases as the disc recedes and the spectrum becomes hard. At the same time, the Comptonized flux decreases with increasing fraction of thermal emission. This could be plausible that episodic jet ejection modified the corona and reduced Comptonized flux. An iron emission line at 6.4 keV is observed in the spectra of all the orbits of observation. We find that the equivalent width of the iron emission line correlates with the photon index, indicating a decrease in the reflection strength as the spectrum becomes hard. We observe that the disc flux does not follow $F_{\rm DBB}-T^{4}$ relation.

Context. Extended gamma-ray TeV emission (TeV halos) around middle-aged pulsars has been detected. A proposed model to explain these TeV halos is that electrons from a degree-wide Pulsar Wind Nebula (PWN) get up-scattered by cosmic microwave background photons through inverse Compton processes. However, no X-ray degree-wide faint diffuse PWNe have been found around these middle-aged pulsars in previous X-ray observations. Aims. We have performed a search for degree wide PWNe around Geminga, PSR B0656+14, B0540+23, J0633+0632, and J0631+1036, using data from the first four consecutive Spectrum Roentgen Gamma/eROSITA all-sky surveys. In order to better understand the mechanisms underlying the formation of TeV halos, we investigated the magnetic field strength in the degree wide neighbourhood of those pulsars. Results. We did not detect degree-wide diffuse emission around Geminga, PSR B0656+14, B0540+23, J0633+0632, and J0631+1036, which can be attributed to being powered by the rotation-powered pulsars. Indeed, a close inspection of the data shows that the pulsars of interest are all embedded in diffuse emission from supernova remnants like the Monogem Ring or the Rosetta Nebula, while PSR B0540+23 is located ~2.5 degrees away from the bright Crab pulsar, which shines out the eROSITA point-spread function up to the position of PSR B0540+23 and thus reduced the sensitivity to search for degree wide bright diffuse X-ray emission strongly. Conclusions. Despite the non-detection of any degree-wide PWN surrounding the analysed pulsars, we set flux upper limits to provide useful information on magnetic field strength and its spatial distribution around those pulsars, providing additional constraints to the proposed theory for the formation of TeV halos around pulsars.

The active asteroid 311P is one of the two targets of a planned Chinese asteroid exploration mission Tianwen-2. During 2013, 311P experienced several mass-loss events and exhibited multiple comet-like tails. Here we analyze the morphology and surface brightness of the tails to investigate the dust environment around the nucleus and mechanism of activities. We enhance the features of the tails using image processing techniques to obtain information about the morphology of the tails, and fit processed images to the syndyne-synchrone diagrams. The fitting results give estimations of the upper limits of the durations ($2\sim8$ days) of eruptions and the dust size range ($0.006\sim38.9$ mm) in the tails. The results of surface photometry performed for each dust tail show that the brightness distribution index of each tail ranged from approximately -1.81 to 0 and the dust size distribution indices of 311P's tails ranged from -2.29 to -1.45. The quantity of particles in each tail ranged from 0.5 to 8 $\times10^6$\,kg, which leads to a total dust-loss quantity of $3.0\times10^7$\,kg and a mass loss rate of 1.59 kg s$^{-1}$. Sublimation, continuous impacts or tidal forces of planets are excluded as explanations for 311P's activities, and rotational instability remains a possible activation cause without strong evidence against it.

The large side aperture of the NuSTAR telescope for unfocused photons (so-called stray light) is a known source of rich astrophysical information. To support many studies based on the NuSTAR stray light data, we present a fully automatic method for determining detector area suitable for background analysis and free from any kind of focused X-ray flux. The method's main idea is `a trous' wavelet image decomposition, capable of detecting structures of any spatial scale and shape, which makes the method of general use. Applied to the NuSTAR data, the method provides a detector image region with the highest possible statistical quality, suitable for the NuSTAR stray light studies. We developed an open-source Python nuwavdet package, which implements the presented method. The package contains subroutines to generate detector image region for further stray light analysis and/or to produce a list of detector bad-flagged pixels for processing in the NuSTAR Data Analysis Software for conventional X-ray analysis.

We describe a machine-learning-based surrogate model for reproducing the Bayesian posterior distributions for exoplanet atmospheric parameters derived from transmission spectra of transiting planets with typical retrieval software such as TauRex. The model is trained on ground truth distributions for seven parameters: the planet radius, the atmospheric temperature, and the mixing ratios for five common absorbers: $H_2O$, $CH_4$, $NH_3$, $CO$ and $CO_2$. The model performance is enhanced by domain-inspired preprocessing of the features and the use of semi-supervised learning in order to leverage the large amount of unlabelled training data available. The model was among the winning solutions in the 2023 Ariel Machine Learning Data Challenge.

The early afterglow of a Gamma-ray burst (GRB) can provide critical information on the jet and progenitor of the GRB. The extreme brightness of GRB 221009A allows us to probe its early afterglow in unprecedented detail. In this letter, we report comprehensive observation results of the early afterglow of GRB 221009A (from $T_0$+660 s to $T_0$+1860 s, where $T_0$ is the \textit{Insight}-HXMT/HE trigger time) in X/$\gamma$-ray energy band (from 20 keV to 20 MeV) by \textit{Insight}-HXMT/HE, GECAM-C and \textit{Fermi}/GBM. We find that the spectrum of the early afterglow in 20 keV-20 MeV could be well described by a cutoff power-law with an extra power-law which dominates the low and high energy bands respectively. The cutoff power-law $E_{\rm peak}$ is $\sim$ 30 keV and the power-law photon index is $\sim$ 1.8 throughout the early afterglow phase. By fitting the light curves in different energy bands, we find that a significant achromatic break (from keV to TeV) is required at $T_0$ + 1246$^{+27}_{-26}$ s (i.e. 1021 s since the afterglow starting time $T_{\rm AG}$=$T_0$+225 s), providing compelling evidence of a jet break. Interestingly, both the pre-break and post-break decay slopes vary with energy, and these two slopes become closer in the lower energy band, making the break less identifiable. Intriguingly, the spectrum of the early afterglow experienced a slight hardening before the break and a softening after the break. These results provide new insights into the understanding of this remarkable GRB.

The classical model of massive-star mechanical feedback is based on effects at solar metallicity (Zsun), yet feedback parameters are very different at low metallicity. Metal-poor stellar winds are much weaker, and more massive supernova progenitors likely collapse directly to black holes without exploding. Thus, for ~0.4 Zsun we find reductions in the total integrated mechanical energy and momentum of ~40% and 75%, respectively, compared to values classically expected at solar metallicity. But in particular, these changes effectively delay the onset of mechanical feedback until ages of ~10 Myr. Feedback from high-mass X-ray binaries could slightly increase mechanical luminosity between ages 5-10 Myr, but it is stochastic and unlikely to be significant on this timescale. Stellar dynamical mechanisms remove most massive stars from clusters well before 10 Myr, which would further promote this effect; this process is exacerbated by gas retention implied by weak feedback. Delayed mechanical feedback implies that radiation feedback therefore dominates at early ages, which is consistent with the observed absence of superwinds in some extreme starbursts. This scenario may lead to higher star-formation efficiencies, multiple stellar populations in clusters, and higher Lyman continuum escape. This could explain the giant star-forming complexes in metal-poor galaxies and the small sizes of OB superbubble shells relative to their inferred ages. It could also drive modest effects on galactic chemical evolution, including on oxygen abundances. Thus, delayed low-metallicity mechanical feedback may have broad implications, including for early cosmic epochs.

During slow-roll inflation, non-perturbative transitions can produce bubbles of metastable vacuum. These bubbles expand exponentially during inflation to super-horizon size, and later collapse into black holes when the expansion of the universe is decelerating. In a broad class of models, the inflationary fine-tuning that gives rise to small density fluctuations causes these bubbles to appear only during a time interval that is short compared to the inflationary Hubble time. As a result, despite the fact that the final mass of the black hole is exponentially sensitive to the moment bubbles form during inflation, the resulting primordial black hole mass spectrum can be nearly monochromatic. If the transition occurs near the middle of inflation, the mass can fall in the "asteroid" range 10^17-10^22 g in which all known observations are compatible with black holes comprising 100% of dark matter.

We perform a global fit to the electroweak vertices and 4-fermion operators of the standard model effective field theory in this work using $N_{\rm eff}$ from cosmological probes, as well as data sets from colliders and low-energy experiments. Assuming flavor universality for both leptons and quarks, we find $N_{\rm eff}$ from CMB-S4/HD can improve the fit for several operators by a factor of a few. Without making any flavor assumption instead, $N_{\rm eff}$ can probe a different direction compared with the well-known flat one presented in the neutrino trident process in the $4\mu$ sector. Fresh sensitivity to two extra operators in the $\tau$ sector is also identified at the 10% level. In addition, we find reducing the uncertainty in primordial helium abundance from metal-poor galaxies down below 0.2% could play an important role in deepening our understanding on the free neutron lifetime anomaly.

At present, there is practically no doubt that general relativity is closely related to gravity. Moreover, after the work of Jacobson, Padmanabhan and others, it became clear that a thermodynamic interpretation of Einstein's relativistic equations is possible. On the other hand, we are witnessing the conceptual problems of the SCM (the problem of the cosmological constant, the problem of coincidences) and many years of futile attempts to directly fix the main components of the model (dark energy and dark matter). The combination of these two factors gave rise to a natural desire, at least at the phenomenological level, to build a cosmological model that represents the synthesis of gravity and thermodynamics and does not include components of an unknown nature. It is this model\`uentropic cosmology\`uthat is considered in this review. We have set as our goal, omitting the details that can be found in the references given, to present the conceptual foundations of the model.

Whether astrophysical black holes (BHs) can have charge is a question to be addressed by observations. In the era of gravitational wave (GW) astronomy, one can constrain the charge of a merged BH remnant using the merger-ringdown signal of the GW data. Extending earlier studies, we analyze five GW events in GWTC-3, assuming Kerr-Newman BHs. Our results show no strong evidence for a charged BH, and give a limit on the charge-to-mass-ratio $Q<0.37$ at $90\%$ credible level (CL). Due to the charge-spin degeneracy in the waveform and the limited signal-to-noise ratios (SNRs), it is challenging for LIGO/Virgo/KAGRA observations to provide better constraints. We further simulate data for the Einstein Telescope (ET), where SNRs can be as large as $\sim270$ in the ringdown signal. These simulated events allow us to consider the 220, 221, and 330 ringdown modes altogether, which can help break the charge-spin degeneracy. The analysis of our simulation shows that ET can improve the constraints on the charge-to-mass-ratio to $Q \lesssim 0.2$ at $90\%$ CL with ringdown signals.

In this work, we investigate the electromagnetic energy released by astrophysical black holes within the Kerr-Taub-NUT solution, which describes rotating black holes with a nonvanishing gravitomagnetic charge. In our study, we consider the black holes in the X-ray binary systems GRS 1915+105, GRO J1655-40, XTE J1550-564, A0620-00, H1743-322, and GRS 1124-683. We show that the Kerr-Taub-NUT spacetime can explain the radiative efficiency of these sources inferred from the continuum fitting method (CFM). We also show that, in the framework of the Blandford-Znajeck mechanism, it is possible to reproduce the observed jet power. We unify the results of the two analyses for the selected objects to get more stringent constraints on the spacetime parameters. We show that, as in the case of the Kerr spacetime, the Kerr-Taub-NUT solution cannot simultaneously explain the observed jet power and radiative efficiency of GRS 1915+105.

We propose a novel method for distinguishing the spin of ultralight dark matter (ULDM) using interferometric gravitational wave detectors. ULDM can be a bosonic field of spin-0, 1, or 2, and each induces distinctive signatures in signals. We find that the finite-time traveling effect causes a dominant signal for spin-0 and spin-1 ULDM, while not for spin-2. By using overlap reduction functions (ORF) of multiple detectors, we can differentiate between the spins of ULDM. Furthermore, we point out that the current constraint on the coupling constant of spin-1 ULDM to baryons becomes 30 times weaker when the finite-time light-travel effect on the ORF is taken into account.