Papers for Tuesday, Jul 25 2023

Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we present a new numerical scheme for the first two resistive terms, which we implement in the moving-mesh code AREPO using the single-fluid approximation combined with a new gradient estimation technique based on a least-squares fit per interface. Through various test calculations including the diffusion of a magnetic peak, the structure of a magnetic C-shock, and the damping of an Alfv\'en wave, we show that we can achieve an accuracy comparable to the state-of-the-art code ATHENA++. We apply the scheme to the linear growth of the magnetorotational instability and find good agreement with the analytical growth rates. By simulating the collapse of a magnetised cloud with constant magnetic diffusion, we show that the new scheme is stable even for large density contrasts. Thanks to the Lagrangian nature of the moving mesh method the new scheme is thus well suited for intended future applications where a high resolution in the dense cores of collapsing protostellar clouds needs to be achieved. In a forthcoming work we will extend the scheme to the Hall effect.

We present 307 Type Ia supernova (SN) light curves from the first four years of the TESS mission. We use this sample to characterize the shapes of the early time light curves, measure the rise times from first light to peak, and search for companion star interactions. Using simulations, we show that light curves must have noise $<$10% of the peak to avoid biases in the early time light curve shape, restricting our quantitative analysis to 74 light curves. We find that the mean power law index $t^{\beta_1}$ of the early time light curves is 1.83$\pm$ 0.57 and the mean rise time to peak is 15.7 $\pm$ 3.5 days. We also estimate the underlying population distribution and find a Gaussian component with mean $\beta_1 = 2.29$, width 0.34, and a tail extending to values less than 1.0. We use model comparison techniques to test for the presence of companion interactions. In contrast to recent results in the literature, we find that the data can rarely distinguish between models with and without companion interactions, and caution is needed when claiming detections of early time flux excesses. Nevertheless, we find three high-quality SN light curves that tentatively prefer the addition of a companion interaction model, but the statistical evidence is not robust. We also find two SNe that disfavor the addition of a companion interaction model to a curved power law model. Taking the 74 SNe together, we calculate 3$\sigma$ upper limits on the presence of companion signatures to control for orientation effects that can hide companions in individual light curves. Our results rule out common progenitor systems with companions having Roche lobe radii $>$ 31 R$_{\odot}$ (99.9% confidence level) and disfavor companions having Roche lobe radii $>$ 10 R$_{\odot}$ (95% confidence level). Lastly, we discuss the implications of our results for the intrinsic fraction of single degenerate progenitor systems.

Early results from the JWST-MIRI guaranteed time programs on protostars (JOYS) and disks (MINDS) are presented. Thanks to the increased sensitivity, spectral and spatial resolution of the MIRI spectrometer, the chemical inventory of the planet-forming zones in disks can be investigated with unprecedented detail across stellar mass range and age. Here data are presented for five disks, four around low-mass stars and one around a very young high-mass star. The mid-infrared spectra show some similarities but also significant diversity: some sources are rich in CO2, others in H2O or C2H2. In one disk around a very low-mass star, booming C2H2 emission provides evidence for a ``soot'' line at which carbon grains are eroded and sublimated, leading to a rich hydrocarbon chemistry in which even di-acetylene (C4H2) and benzene (C6H6) are detected (Tabone et al. 2023). Together, the data point to an active inner disk gas-phase chemistry that is closely linked to the physical structure (temperature, snowlines, presence of cavities and dust traps) of the entire disk and which may result in varying CO2/H2O abundances and high C/O ratios >1 in some cases. Ultimately, this diversity in disk chemistry will also be reflected in the diversity of the chemical composition of exoplanets.

The physical processes which make a galaxy a Lyman Alpha Emitter have been extensively studied for the past 25 years. However, the correlations between physical and morphological properties of galaxies and the strength of the Ly$\alpha$ emission line are still highly debated. Therefore, we investigate the correlations between the rest-frame Ly$\alpha$ equivalent width and stellar mass, star formation rate, dust reddening, metallicity, age, half-light semi-major axis, S\'ersic index and projected axis ratio in a sample of 1578 galaxies in the redshift range $2 \leq z \leq 7.9$ from the GOODS-S, UDS and COSMOS fields. From the large sample of Ly$\alpha$ emitters (LAEs) in the dataset we find that LAEs are typically common main sequence star forming galaxies which show stellar mass $ \leq 10^9 \text{M}_{\odot}$, star formation rate $ \leq 10^{0.5} \text{M}_{\odot}/\text{yr}$, $E(B-V) \leq 0.2$ and half-light semi-major axis $\leq 1 \text{kpc}$. Building on these findings we develop a new method based on Random Forest (i.e. a Machine Learning classifier) in order to select galaxies which have the highest probability of being Ly$\alpha$ emitters. When applied to a population in the redshift range $z \in [2.5, 4.5]$, our classifier holds a $(80 \pm 2)\%$ accuracy and $(73 \pm 4)\%$ precision. At higher redshifts ($z \in [4.5, 6]$), we obtain a $73\%$ accuracy and a $80\%$ precision. These results highlight it is possible to overcome the current limitations in assembling large samples of LAEs by making informed predictions that can be used for planning future large scale spectroscopic surveys.

We present new astrometric and polarimetric observations of flares from Sgr A* obtained with GRAVITY, the near-infrared interferometer at ESO's Very Large Telescope Interferometer (VLTI), bringing the total sample of well-covered astrometric flares to four and polarimetric ones to six, where we have for two flares good coverage in both domains. All astrometric flares show clockwise motion in the plane of the sky with a period of around an hour, and the polarization vector rotates by one full loop in the same time. Given the apparent similarities of the flares, we present a common fit, taking into account the absence of strong Doppler boosting peaks in the light curves and the EHT-measured geometry. Our results are consistent with and significantly strengthen our model from 2018: We find that a) the combination of polarization period and measured flare radius of around nine gravitational radii ($9 R_g \approx 1.5 R_{ISCO}$, innermost stable circular orbit) is consistent with Keplerian orbital motion of hot spots in the innermost accretion zone. The mass inside the flares' radius is consistent with the $4.297 \times 10^6 \; \text{M}_\odot$ measured from stellar orbits at several thousand $R_g$. This finding and the diameter of the millimeter shadow of Sgr A* thus support a single black hole model. Further, b) the magnetic field configuration is predominantly poloidal (vertical), and the flares' orbital plane has a moderate inclination with respect to the plane of the sky, as shown by the non-detection of Doppler-boosting and the fact that we observe one polarization loop per astrometric loop. Moreover, c) both the position angle on sky and the required magnetic field strength suggest that the accretion flow is fueled and controlled by the winds of the massive, young stars of the clockwise stellar disk 1-5 arcsec from Sgr A*, in agreement with recent simulations.

Compact groups of galaxies (CGs) show members with morphological disturbances, mainly products of galaxy-galaxy interactions, thus making them ideal systems to study galaxy evolution, in high-density environment. To understand how this environment affects the properties of galaxies, we select a sample of 340 CGs in the Stripe 82 region, for a total of 1083 galaxies, and a sample of 2281 field galaxies as a control sample. By performing a multi-wavelength morphological fitting process using S-PLUS data, we divide our sample into early-type (ETG), late-type (LTG), and transition galaxies using the r-band S\'ersic index and the colour (u-r). We find a bimodal distribution in the plane of the effective radius-S\'ersic index, where a secondary "peculiar" galaxy population of smaller and more compact galaxies is found in CGs, which is not observed in the control sample. This indicates that galaxies are undergoing a morphological transformation in CGs. In addition, we find significant statistical differences in the distribution of specific Star Formation Rate (sSFR) when we compare both environments for LTGs and ETGs. We also find a higher fraction of quenched galaxies and a lower median sSFR in CGs than in the control sample, suggesting the existence of environmental effects favoring the cessation of star formation, regardless of galaxy type. Our results support the notion that CGs promote morphological and physical transformations, highlighting their potential as ideal systems for galaxy pre-processing.

Current methods of identifying the ionizing source of nebular emission in galaxies are well defined for the era of single fiber spectroscopy, but still struggle to differentiate the complex and overlapping ionization sources in some galaxies. With the advent of integral field spectroscopy, the limits of these previous classification schemes are more apparent. We propose a new method for distinguishing the ionizing source in resolved galaxy spectra by use of a multi-dimensional diagnostic diagram that compares emission line ratios with velocity dispersion on a spaxel by spaxel basis within a galaxy. This new method is tested using the SAMI Galaxy Survey Data Release 3, which contains 3068 galaxies at z $<$ 0.12. Our results are released as ionization maps available alongside the SAMI DR3 public data. Our method accounts for a more diverse range of ionization sources than the standard suite of emission line diagnostics; we find 1433 galaxies with significant contribution from non-star-forming ionization using our improved method as compared to 316 galaxies identified using only emission line ratio diagnostics. Within these galaxies, we further identify 886 galaxies hosting unique signatures inconsistent with standard ionization by H2 regions, AGN, or shocks. These galaxies span a wide range of masses and morphological types and comprise a sizable portion of the galaxies used in our sample. With our revised method, we show that emission line diagnostics alone do not adequately differentiate the multiple ways to ionize gas within a galaxy.

We quantify the cosmological spread of baryons relative to their initial neighboring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighboring distribution owing to chaotic gravitational dynamics on spatial scales comparable to their host dark matter halo. In contrast, gas in hydrodynamic simulations spreads much further from the initial neighboring dark matter owing to feedback from supernovae (SNe) and Active Galactic Nuclei (AGN). We show that large-scale baryon spread is very sensitive to model implementation details, with the fiducial \textsc{SIMBA} model spreading $\sim$40\% of baryons $>$1\,Mpc away compared to $\sim$10\% for the IllustrisTNG and \textsc{ASTRID} models. Increasing the efficiency of AGN-driven outflows greatly increases baryon spread while increasing the strength of SNe-driven winds can decrease spreading due to non-linear coupling of stellar and AGN feedback. We compare total matter power spectra between hydrodynamic and paired $N$-body simulations and demonstrate that the baryonic spread metric broadly captures the global impact of feedback on matter clustering over variations of cosmological and astrophysical parameters, initial conditions, and galaxy formation models. Using symbolic regression, we find a function that reproduces the suppression of power by feedback as a function of wave number ($k$) and baryonic spread up to $k \sim 10\,h$\,Mpc$^{-1}$ while highlighting the challenge of developing models robust to variations in galaxy formation physics implementation.

Ultra-massive white dwarf stars are currently being discovered at a considerable rate, thanks to surveys such as the {\it Gaia} space mission. These dense and compact stellar remnants likely play a major role in type Ia supernova explosions. It is possible to probe the interiors of ultra-massive white dwarfs through asteroseismology. In the case of the most massive white dwarfs, General Relativity could affect their structure and pulsations substantially. In this work, we present results of relativistic pulsation calculations employing relativistic ultra-massive ONe-core white dwarf models with hydrogen-rich atmospheres and masses ranging from $1.29$ to $1.369 M_{\odot}$ with the aim of assessing the impact of General Relativity on the adiabatic gravity ($g$)-mode period spectrum of very-high mass ZZ Ceti stars. Employing the relativistic Cowling approximation for the pulsation analysis, we find that the critical buoyancy (Brunt-V\"ais\"al\"a) and acoustic (Lamb) frequencies are larger for the relativistic case, compared to the Newtonian case, due to the relativistic white dwarf models having smaller radii and higher gravities for a fixed stellar mass. In addition, the $g$-mode periods are shorter in the relativistic case than in the Newtonian computations, with relative differences of up to $\sim 50$ \% for the highest-mass models ($1.369 M_{\odot}$) and for effective temperatures typical of the ZZ Ceti instability strip. Hence, the effects of General Relativity on the structure, evolution, and pulsations of white dwarfs with masses larger than $\sim 1.29 M_{\odot}$ cannot be ignored in the asteroseismological analysis of ultra-massive ZZ Ceti stars.

We present observational evidence of the impact of triaxiality on radial profiles that extend to 40~Mpc from galaxy cluster centres in optical measurements. We perform a stacked profile analysis from a sample of thousands of nearly relaxed galaxy clusters from public data releases of the Dark Energy Survey (DES) and the Dark Energy Camera Legacy Survey (DECaLS). Using the central galaxy elliptical orientation angle as a proxy for galaxy cluster orientation, we measure cluster weak lensing and excess galaxy density axis-aligned profiles, extracted along the central galaxy's major or minor axes on the plane-of-the-sky. Our measurements show a $\gtrsim2-3\sigma$ difference per radial bin between the normalized axis-aligned profiles. The profile difference between each axis-aligned profile and the azimuthally averaged profile ($\sim\pm10-20\%$ along major/minor axis) appears inside the clusters ($\sim0.4$ Mpc) and extends to the large-scale structure regime ($\sim10-20$ Mpc). The magnitude of the difference appears to be relatively insensitive to cluster richness and redshift, and extends further out in the weak lensing surface mass density than in the galaxy overdensity. Looking forward, this measurement can easily be applied to other observational or simulation datasets and can inform the systematics in cluster mass modeling related to triaxiality. We expect imminent upcoming wide-area deep surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), to improve our quantification of optical signatures of cluster triaxiality.

Merging of galaxy clusters are some of the most energetic events in the Universe, and they provide a unique environment to study galaxy evolution. We use a sample of 84 merging and relaxed SPT galaxy clusters candidates, observed with the Dark Energy Camera in the $0.11<z<0.88$ redshift range, to build colour-magnitude diagrams to characterize the impact of cluster mergers on the galaxy population. We divided the sample between relaxed and disturbed, and in two redshifts bin at $z = 0.55$. When comparing the high-z to low-z clusters we find the high-z sample is richer in blue galaxies, independently of the cluster dynamical state. In the high-z bin we find that disturbed clusters exhibit a larger scatter in the Red Sequence, with wider distribution and an excess of bluer galaxies compared to relaxed clusters, while in the low-z bin we find a complete agreement between the relaxed and disturbed clusters. Our results support the scenario in which massive cluster halos at $z<0.55$ galaxies are quenched as satellites of another structure, i.e. outside the cluster, while at $z \geq 0.55$ the quenching is dominated by in-situ processes.

We study the effect an atmosphere has on pebble orbits and spin build-up on a planet's surface during pebble accretion in the extreme case of a static atmosphere. We numerically integrate the equations of motion of pebbles in a planar, global frame with a planet, a central star and gas from a protoplanetary disc. An adiabatic atmosphere is then placed around the planet, and the spin deposited onto the planet's surface is measured. These simulations are evaluated for different distances to the star, Stokes numbers, and planet masses. Pebble feedback to the gas is not taken into account. We find that a static atmosphere dampens the spin the planet's surface receives by absorbing part of the angular momentum of the pebbles and circularising their orbits. This could prevent the excessive spin values predicted in some 3D pebble accretion simulations without an atmosphere. For planets larger than 0.5 Earth masses, a stationary atmosphere absorbs all angular momentum, leaving no spin for the surface. Significant quantities of angular momentum are stored in the inner and intermediate atmosphere ($<0.3$ Bondi radii). Depending on the atmospheric and disc model, this spin could be transported either to the disc through atmospheric recycling or to the planet through drag between the surface and the atmosphere. Further research is required to quantify the spin transfer within the atmosphere.

We investigated the molecular environment of AFGL 2298, an obscured Galactic Luminous blue variable (LBV) star which hosts a highly structured circumstellar environment with hints of multiple mass-loss events in the last few $10^4$ a. We present spectral line observations of AFGL 2298 at 1 and 3 mm performed with the IRAM 30m radio telescope. Furthermore, we report the detection of several carbon- and nitrogen- bearing species (CO, $^{13}$CO, C$^{18}$O, C$^{17}$O, HCO$^+$, HCN, HNC, H$^{13}$CO$^+$, CN, N$_2$H$^+$, and C$_2$H) in the surroundings of AFGL 2298. In addition, we identified three velocity components that clearly stand out from the Galactic background. The morphology, kinematics, masses and isotopic ratios, together with a comparative study of the fractional abundances, lead us to suggest that two of these components (36 and 70 km/s) have a stellar origin. The other component (46 km/s) most likely traces swept-up interstellar material, probably harbouring also a photon-dominated region. The first inventory of the circumstellar molecular gas around AFGL 2298 is provided. The results are compatible with the hypothesis of former mass-loss events, produced before the one that created the infrared nebula. There are chemical hints of the presence of ejected stellar material, and also swept up gas. These findings will help to better understand the mass-loss history of this class of evolved massive stars, which heavily influence the overall chemical evolution of the Galaxy.

The heaviest black holes discovered through gravitational waves have masses that are difficult to explain with current standard stellar models. This discrepancy may be due to a series of hierarchical mergers, where the observed black holes are themselves the products of previous mergers. Here we present a method to estimate the masses and spins of previous generations of black holes based on the masses and spins of black holes in a binary. Examining the merger GW190521, we find that assuming black hole spins that are consistent with those of merger remnants will alter the reconstructed ancestral spins when compared to results with uninformed priors. At the same time, the inclusion of black hole spins does not significantly affect the mass distributions of the ancestral black holes.

We present positions, proper motions, and near-infrared photometry for 966 known objects with spectral types later than M observed as part of the the UKIRT Hemisphere Survey (UHS). We augment the photometry and astrometry from UHS with information from Gaia DR3, Pan-STARRS DR2, and CatWISE 2020 to produce a database of homogeneous photometry and astrometry for this sample. The multi-epoch survey strategy of UHS allows us to determine proper motions for most sources, with a median proper motion uncertainty of $\sim$3.6 mas yr$^{-1}$. Our UHS proper motion measurements are generally in good agreement with those from Gaia DR3, Pan-STARRS, and CatWISE 2020, with UHS proper motions typically more precise than those from CatWISE 2020 and Pan-STARRS but not Gaia DR3. We critically analyze publicly available spectra for 406 members of this sample and provide updated near-infrared spectral types for $\sim$100 objects. We determine typical colors as a function of spectral type and provide absolute magnitude vs. spectral type relations for UHS $J$- and $K$-band photometry. Using newly determined proper motions, we highlight several objects of interest, such as objects with large tangential velocities, widely separated co-moving companions, and potential members of young nearby associations.

Interstellar objects (ISOs) are small bodies that can travel through our solar system from other star systems. When present in our solar system, they represent an opportunity to study the properties and origins of these objects, as well as the potential for cross-pollination of material between star systems. With current propulsion technology, rendezvous with these objects is likely infeasible, and thus the maximum science return results from a rapid response flyby and impactor. However, while trajectories to ISOs may be feasible, their potentially high ephemeris uncertainties and high-speed hyperbolic orbits present significant challenges to navigation. In this paper we assess these challenges by modeling the uncertainties of reachable synthetic ISOs as a function of time, as derived by measurements from ground observatories and an approaching spacecraft. From these uncertainties we derive the final delivery accuracy of fast flyby spacecraft to the ISO and required statistical delta-v for navigation. We find that these two challenges can lead to hundreds of meters-per-second or even kilometers-per-second of required statistical delta-v for navigation, reduce delivery accuracy to hundreds of kilometers, and make autonomous navigation a requirement.

The Keck Planet Imager and Characterizer (KPIC) is an instrument at the Keck II telescope that enables high-resolution spectroscopy of directly imaged exoplanets and substellar companions. KPIC uses single-mode fibers to couple the adaptive optics system to Keck's near-infrared spectrometer (NIRSPEC). However, KPIC's sensitivity at small separations is limited by the leakage of stellar light into the fiber. Speckle nulling uses a deformable mirror to destructively interfere starlight with itself, a technique typically used to reduce stellar signal on a focal-plane imaging detector. We present the first on-sky demonstration of speckle nulling through an optical fiber with KPIC, using NIRSPEC to collect exposures that measure speckle phase for quasi-real-time wavefront control while also serving as science data. We repeat iterations of measurement and correction, each using at least 5 exposures. We show a decrease in the on-sky leaked starlight by a factor of 2.6 to 2.8 in the targeted spectral order, at a spatial separation of 2.0 {\lambda}/D in K-band. This corresponds to an estimated factor of 2.6 to 2.8 decrease in the required exposure time to reach a given SNR, relative to conventional KPIC observations. The performance of speckle nulling is limited by instability in the speckle phase: when the loop is opened, the null-depth degrades by a factor of 2 on the timescale of a single phase measurement, which would limit the suppression that can be achieved. Future work includes exploring gradient-descent methods, which may be faster and thereby able to achieve deeper nulls. In the meantime, the speckle nulling algorithm demonstrated in this work can be used to decrease stellar leakage and improve the signal-to-noise of science observations.

We determine the efficacy of the kinematic Sunyaev-Zel'dovich signal extraction pipeline, using pairwise kSZ measurements, in recovering unbiased estimates of the signal and inference of the associated optical depth. We consider the impact of cluster co-alignments along the line of sight, the modeling of baryonic clustering, and the presence of diffuse gas, as well as instrument beam convolution and noise. We demonstrate that two complementary approaches, aperture photometry, and a matched filter, can be used to recover an unbiased estimate of the cluster kSZ signal and the associated optical depth. Aperture photometry requires a correction factor accounting for the subtraction of signal in the annulus while the matched filter requires a tuning of the signal template profile. We show that both of these can be calibrated from simulated survey data. The optical depth estimates are also consistent with those inferred from stacked thermal SZ measurements. We apply the approaches to the publicly available Atacama Cosmology Telescope (ACT) data. The techniques developed here provide a promising method to leverage upcoming kSZ measurements, from ACT, Simons Observatory, CCAT, and CMB-S4 with spectroscopic galaxy surveys from DESI, Euclid, and Roman, to constrain cosmological properties of the dark energy, gravity, and neutrino masses.

In this study, we propose a reparameterization of a specific viable $f(R)$ gravity model to represent it as a perturbation of the $\Lambda$CDM model. The $f(R)$ gravity model under consideration includes two parameters, $b$ and $n$, which control how close the proposed model can be to $\Lambda$CDM, allowing for arbitrary proximity. Furthermore, it is shown that the Hu-Sawicki (HS) model is a limiting case of this reparameterized model. Following the existing literature, we also derive an analytical approximation for the expansion rate $H(z)$, which shows an excellent agreement between this analytical approximation and the numerical solution over a wide range of redshifts for realistic values of the deviation parameter $b$. By appropriately selecting values for the model parameters, we plot the cosmological parameters $w_{\rm{DE}}$, $w_{\rm{eff}}$, $\Omega_{\rm{DE}}$, and $H(z)$, as well as the statefinder quantities $q$, $j$, $s$, and $Om(z)$. We find that their present values (at $z=0$) are consistent with the observations from Planck 2018 and the values predicted by the $\Lambda$CDM model. It is important to note that the examined cosmological and statefinder parameters do not exhibit significant oscillations of effective dark energy, which could lead to singular and unphysical solutions at high redshifts. This anomalous behavior has been avoided here by utilizing the approximate analytical solution for $H(z)$. Additionally, we conduct a detailed analysis of the evolution of matter density perturbations within the introduced $f(R)$ gravity model. The results demonstrate that this viable $f(R)$ gravity model is practically indistinguishable from the $\Lambda$CDM model at the background level.

We report a broad-band analysis of a Galactic X-ray source, CXOGBS J174517.0-321356 (J1745), with a 614-second periodicity. Chandra discovered the source in the direction of the Galactic Bulge. Gong (2022) proposed J1745 was either an intermediate polar (IP) with a mass of ~1 $M_{\odot}$, or an ultra-compact X-ray binary (UCXB). By jointly fitting XMM-Newton and NuSTAR spectra, we rule out a UCXB origin. We have developed a physically realistic model that considers finite magnetosphere radius, X-ray absorption from the pre-shock region, and reflection from the WD surface to determine the IP properties, especially its WD mass. To assess systematic errors on WD mass measurement, we consider a broad range of specific accretion rates ($\dot{m}$ = 0.6 - 44 g\cm$^2$\s) based on the uncertain source distance (d = 3-8 kpc) and fractional accretion area (f = 0.001-0.025). Our model properly implements the fitted accretion column height in the X-ray reflection model and accounts for the underestimated mass accretion rate due to the (unobserved) soft X-ray blackbody and cyclotron cooling emissions. We found that the lowest accretion rate of $\dot{m}$ = 0.6 g\cm$^2$\s, which corresponds to the nearest source distance and maximum f value, yield the WD mass of $(0.92\pm0.08) M_{\odot}$. However, if the accretion rate is $\dot{m}$ > ~3 g\cm$^2$\s, the WD mass is robustly measured to be $(0.81\pm0.06) M_{\odot}$, nearly independent of $\dot{m}$. The derived WD mass range is consistent with the mean WD mass of nearby IPs. Assuming spin equilibrium between the WD and accretion disk, we constrained the WD magnetic field to B > ~7 MG, indicating that it could be a highly magnetized IP. Our analysis presents the most comprehensive methodology for constraining the WD mass and B-field of an IP by consolidating the effects of cyclotron cooling, finite magnetospheric radius, and accretion column height.

When cool clouds are ram-pressure accelerated by a hot supersonic galactic wind, some of the clouds may be shredded by hydrodynamical instabilities and incorporated into the hot flow. Recent one-dimensional steady-state calculations show how cool cloud entrainment directly affects the bulk thermodynamics, kinematics, and observational characteristics of the hot gas. In particular, mass-loading decelerates the hot flow and changes its entropy. Here, we investigate the stability of planar and spherical mass-loaded hot supersonic flows using both perturbation analysis and three-dimensional time-dependent radiative hydrodynamical simulations. We show that mass-loading is stable over a broad range of parameters and that the 1D time-steady analytic solutions exactly reproduce the 3D time-dependent calculations, provided that the flow does not decelerate sufficiently to become subsonic. For higher values of the mass-loading, the flow develops a sonic point and becomes thermally unstable, rapidly cooling and forming elongated dense cometary filaments. We explore the mass-loading parameters required to reach a sonic point and the radiative formation of these filaments. For certain approximations, we can derive simple analytic criteria. In general a mass-loading rate similar to the initial mass outflow rate is required. In this sense, the destruction of small cool clouds by a hot flow may ultimately spontaneously generate fast cool filaments, as observed in starburst superwinds. Lastly, we find that the kinematics of filaments is sensitive to the slope of the mass-loading function. Filaments move faster than the surrounding wind if mass-loading is over long distances whereas filaments move slower than their surroundings if mass-loading is abrupt.

Classifying Active Galactic Nuclei (AGN) is a challenge, especially for BL Lac Objects (BLLs), which are identified by their weak emission line spectra. To address the problem of classification, we use data from the 4th Fermi Catalog, Data Release 3. Missing data hinders the use of machine learning to classify AGN. A previous paper found that Multiple Imputation by Chain Equations (MICE) imputation is useful for estimating missing values. Since many AGN have missing redshift and the highest energy, we use data imputation with MICE and K-nearest neighbor (kNN) algorithm to fill in these missing variables. Then, we classify AGN into the BLLs or the Flat Spectrum Radio Quasars (FSRQs) using the SuperLearner, an ensemble method that includes several classification algorithms like logistic regression, support vector classifiers, Random Forests, Ranger Random Forests, multivariate adaptive regression spline (MARS), Bayesian regression, Extreme Gradient Boosting. We find that a SuperLearner model using MARS regression and Random Forests algorithms is 91.1% accurate for kNN imputed data and 91.2% for MICE imputed data. Furthermore, the kNN-imputed SuperLearner model predicts that 892 of the 1519 unclassified blazars are BLLs and 627 are Flat Spectrum Radio Quasars (FSRQs), while the MICE-imputed SuperLearner model predicts 890 BLLs and 629 FSRQs in the unclassified set. Thus, we can conclude that both imputation methods work efficiently and with high accuracy and that our methodology ushers the way for using SuperLearner as a novel classification method in the AGN community and, in general, in the astrophysics community.

We present a comprehensive study of the rotational and emission properties of PSR J0738$-$4042 using a combination of observations taken by the Deep Space Network, Hartebeesthoek, Parkes (Murriyang) and Molonglo observatories between 1972 and 2023. Our timing of the pulsar is motivated by previously reported profile/spin-down events that occurred in September 2005 and December 2015, which result in an anomalously large braking index of $n = 23300 \pm 1800$. Using a Gaussian process regression framework, we develop continuous models for the evolution of the pulsar spin-down rate ($\dot{\nu}$) and profile shape. We find that the pulse profile variations are similar regardless of radio observing frequency and polarisation. Small-scale differences can be ascribed to changes in the interstellar medium along the line of sight and frequency-dependent changes in magnetospheric radio emission height. No new correlated spin-down or profile events were identified in our extended dataset. However, we found that the disappearance of a bright emission component in the leading edge of archival profiles between 1981-1988 was not associated with a substantial change in $\dot{\nu}$. This marks a notable departure from the previous profile/spin-down events in this pulsar. We discuss the challenges these observations pose for physical models and conclude that interactions between the pulsar and in-falling asteroids or a form of magnetospheric state-switching with a long periodicity are plausible explanations.

[Abridged] It is expected that the ongoing and future space-borne planet survey missions including TESS, PLATO, and Earth 2.0 will detect thousands of small to medium-sized planets via the transit technique, including over a hundred habitable terrestrial rocky planets. To conduct a detailed study of these terrestrial planets, particularly the cool ones with wide orbits, the exoplanet community has proposed various follow-up missions. The currently proposed ESA mission ARIEL is capable of characterization of planets down to warm super-Earths mainly using transmission spectroscopy. The NASA 6m UV/Opt/NIR mission proposed in the Astro2020 Decadal Survey may further tackle down to habitable rocky planets, and is expected to launch around 2045. In the meanwhile, China is funding a concept study of a 6-m class space telescope named Tianlin (A UV/Opt/NIR Large Aperture Space Telescope) that aims to start its operation within the next 10-15 years and last for 5+ years. Tianlin will be primarily aimed to the discovery and characterization of rocky planets in the habitable zones (HZ) around nearby stars and to search for potential biosignatures mainly using the direct imaging method. Transmission and emission spectroscopy at moderate to high resolution will be carried out as well on a population of exoplanets to strengthen the understanding of the formation and evolution of exoplanets. It will also carry out in-depth studies of the cosmic web and early galaxies, and constrain the nature of the dark matter and dark energy. We describe briefly the primary scientific motivations and main technical considerations based on our preliminary simulation results. We find that a monolithic off-axis space telescope with a primary mirror diameter larger than 6m equipped with a high contrast chronograph can identify water in the atmosphere of a habitable-zone Earth-like planet around a Sun-like star.

Recent radiation-magnetohydrodynamic simulations of active galactic nuclei predict the presence of the disk winds, which may get unstable and turn into fragmented clumps far from the central black hole. These inner winds and the outer clumps may be observed as the ultrafast outflows (UFOs) and the partial absorbers, respectively. However, it is challenging to observationally constrain their origins because of the complicated spectral features and variations. To resolve such degeneracies of the clumpy absorbers and other components, we developed a novel ``spectral-ratio model fitting'' technique that estimates the variable absorbing parameters from the ratios of the partially absorbed spectra to the non-absorbed one, canceling the complex non-variable spectral features. We applied this method to the narrow-line Seyfert 1 galaxy \iras observed by \xmm in 2016 for $\sim$1.5 Ms. As a result, we found that the soft spectral variation is mostly caused by changes in the partial covering fraction of the mildly-ionized clumpy absorbers, whose outflow velocities are similar to those of the UFO ($\sim$0.2--0.3 $c$). Furthermore, the velocities of the clumpy absorbers and UFOs increase similarly with the X-ray fluxes, consistent with the change in the UV-dominant continuum flux. We also discovered a striking correlation between the clump covering fraction and the equivalent width of the UFO absorption lines, which indicates that increasing the outflow in the line-of-sight lead to more prominent UFOs and more partial absorption. These findings strongly suggest that the clumpy absorbers and the UFO share the same origin, driven by the same UV-dominant continuum radiation.

With new calibration data, thermal emission from the jet of radio galaxy M87 is studied with about 700 ks archival data with Chandra. For nucleus, HST-1, knot D, X-ray energy spectra is well fitted with a power law. However, For knot A, a power law model is rejected with a high significance and an X-ray energy spectra is well fitted with a combination model of a power law and an apec model of 0.2 keV and a metal abundance 0.00. Thermal emission from knot A is confirmed.

Terrestrial and sub-Neptune planets are expected to form in the inner ($<10~$AU) regions of protoplanetary disks. Water plays a key role in their formation, although it is yet unclear whether water molecules are formed in-situ or transported from the outer disk. So far Spitzer Space Telescope observations have only provided water luminosity upper limits for dust-depleted inner disks, similar to PDS 70, the first system with direct confirmation of protoplanet presence. Here we report JWST observations of PDS 70, a benchmark target to search for water in a disk hosting a large ($\sim54~$AU) planet-carved gap separating an inner and outer disk. Our findings show water in the inner disk of PDS 70. This implies that potential terrestrial planets forming therein have access to a water reservoir. The column densities of water vapour suggest in-situ formation via a reaction sequence involving O, H$_2$, and/or OH, and survival through water self-shielding. This is also supported by the presence of CO$_2$ emission, another molecule sensitive to UV photodissociation. Dust shielding, and replenishment of both gas and small dust from the outer disk, may also play a role in sustaining the water reservoir. Our observations also reveal a strong variability of the mid-infrared spectral energy distribution, pointing to a change of inner disk geometry.

We use a full general relativistic framework to study the self-similar expansion of bubbles of the stable phase into a flat Friedmann-Lema\^itre-Robertson-Walker Universe in a first order phase transition in the early Universe. With a simple linear barotropic equation of state in both phases, and in the limit of a phase boundary of negligible width, we find that self-similar solutions exist, which are qualitatively similar to the analogous solutions in Minkowski space, but with distinguishing features. Rarefaction waves extend to the centre of the bubble, while spatial sections near the centre of the bubble have negative curvature. Gravitational effects redistribute the kinetic energy of the fluid around the bubble, and can change the kinetic energy fraction significantly. The kinetic energy fraction of the gravitating solution can be enhanced over the analogous Minkowski solution by as much as $\mathcal{O}(1)$, and suppressed by a factor as larger as $\mathcal{O}(10)$ in case of fast detonations. The amount of negative spatial curvature at the centre of the bubble is of the same order of magnitude of the naive expectation based on considerations of the energy density perturbation in Minkowski solutions, with gravitating deflagrations less negatively curved, and detonations more. We infer that general relativistic effects might have a significant impact on accurate calculations of the gravitational wave power spectrum when the bubble size becomes comparable to the cosmological Hubble radius, affecting the primary generation from the fluid shear stress, and inducing secondary generation by scalar perturbations.

Binary black hole (BBH) systems detected via gravitational-wave (GW) emission are a recently opened astrophysical frontier with many unknowns and uncertainties. Accurate reconstruction of the binary distribution with as few assumptions as possible is desirable for inference on formation channels and environments. Most population analyses have, though, assumed a power law in binary mass ratio $q$, and/or assumed a universal $q$ distribution regardless of primary mass. Kernel density estimation (KDE)-based methods allow us to dispense with such assumptions and directly estimate the joint binary mass distribution. We deploy a self-consistent iterative method to estimate this full BBH mass distribution, finding local maxima in primary mass consistent with previous investigations and a secondary mass distribution with a partly independent structure, inconsistent with both power laws and with a constant function of $q$. We find a weaker preference for near-equal mass binaries than in most previous investigations; instead, the secondary mass has its own "spectral lines" at slightly lower values than the primary, and we observe an anti-correlation between primary and secondary masses around the ~$10M_\odot$ peak.

Real extreme/intermediate mass ratio inspiral(E/IMRI) systems are likely to contain large accretion disks which could be as massive as the central supermassive black hole. Therefore, contrary to its ideal model, a real E/IMRI system contains a third important component: the accretion disk. We study the influence of these disks on the emitted GW profile and its detectability through proposed LISA observation. We use a semi-relativistic formalism in the Kerr background (Gair & Glampedakis 2006; Barausse & Rezzolla 2008) for the case of transonic accretion flow which is a potential candidate to describe the accretion flows around AGN. The hydrodynamic drag of the disks modified the motion of the companion as a result the emitted wave changes in amplitude and phase. We found that these changes are detectable through the last few years of observation by LISA (in some cases as small as six months) for EMRIs residing within 3 GPc from the detector and for the accretion rate of the primary black hole of the order of $\dot{M}=1 \dot{M}_{Edd}$. These choices of parameter values are consistent with real systems. The drag effect and hence the detectability of the emitted GW is sensitive to the hydrodynamical model of the disk. Therefore such observations will help one to identify the nature of the accretion flow and verify various paradigms of accretion physics.

The in-ice or in-water Cherenkov neutrino telescope such as IceCube has already proved its power in measuring the Glashow resonance by searching for the bump around $E^{}_{\rm \nu} = 6.3~{\rm PeV}$ arising from the $W$-boson production. In the next few decades, there are many proposals that observe cosmic tau neutrinos with extensive air showers, also known as tau neutrino telescopes. As has been recognized, the air shower telescope is in principle sensitive to the Glashow resonance via the channel $W \to \tau \nu^{}_{\tau}$ followed by the tau decay in the air. However, with a thorough numerical analysis we have identified several limitations for those telescopes on hunting the resonance. If ultrahigh-energy neutrinos are dominantly produced from the meson decay, it will be statistically difficult for a rather advanced proposal, such as TAMBO with a geometric area around $500~{\rm km^2}$, to discriminate the Glashow resonance induced by $\overline{\nu}^{}_{e}$ from the intrinsic $\nu^{}_{\tau}/\overline{\nu}^{}_{\tau}$ background. The discovery significance is only around $1\sigma$ considering the flux parameters measured by IceCube as the input. Nevertheless, the significance will be improved to $90\%$ if PeV neutrinos mainly originate from the neutron decay, which is, however, thought to be only a subdominant neutrino source. The presence of new physics can also increase the significance. Compared to the in-ice or in-water telescope, the challenge for the Glashow resonance search is ascribed to several factors: (i) a suppressed branching ratio of $11\%$ for the decay $W \to \tau \nu^{}_{\tau}$; (ii) the smearing effect and the reduced acceptance because the daughter neutrino takes away $\langle y \rangle \sim 75\%$ of the energy from the $W$ decay; (iii) a large attenuation effect for Earth-skimming neutrinos with the resonance.

The PLATO mission is scheduled for launch in 2026. This study aims to estimate the number of exoplanets that PLATO can detect as a function of planetary size and period, stellar brightness, and observing strategy options. Deviations from these estimates will be informative of the true occurrence rates of planets, which helps constraining planet formation models. For this purpose, we developed the Planet Yield for PLATO estimator (PYPE), which adopts a statistical approach. We apply given occurrence rates from planet formation models and from different search and vetting pipelines for the Kepler data. We estimate the stellar sample to be observed by PLATO using a fraction of the all-sky PLATO stellar input catalog (PIC). PLATO detection efficiencies are calculated under different assumptions that are presented in detail in the text. The results presented here primarily consider the current baseline observing duration of four years. We find that the expected PLATO planet yield increases rapidly over the first year and begins to saturate after two years. A nominal (2+2) four-year mission could yield about several thousand to several tens of thousands of planets, depending on the assumed planet occurrence rates. We estimate a minimum of 500 Earth-size (0.8-1.25 RE) planets, about a dozen of which would reside in a 250-500d period bin around G stars. We find that one-third of the detected planets are around stars bright enough (V $\leq 11$) for RV-follow-up observations. We find that a three-year-long observation followed by 6 two-month short observations (3+1 years) yield roughly twice as many planets as two long observations of two years (2+2 years). The former strategy is dominated by short-period planets, while the latter is more beneficial for detecting earths in the habitable zone.

In this followup analysis, we update previous constraints on the Transitional Planck Mass (TPM) modified gravity model using the latest version of EFTCAMB and provide new constraints using SPT and Planck anisotropy data along with Planck CMB lensing, BAO, SNe Ia, and an $H_0$ prior from local measurements. We find that large shifts in the Planck mass lead to large suppression of power on small scales that is disfavored by both SPT and Planck. Using only SPT TE-EE data, this suppression of power can be compensated for by an upward shift of the scalar index to $n_s = 1.003 \pm 0.016$ resulting in $H_0 = 71.94^{+0.86}_{-0.85}$ kms$^{-1}$Mpc$^{-1}$ and a $\sim7\%$ shift in the Planck mass. Including Planck TT $\ell \leq 650$ and Planck TE-EE data restricts the shift to be $<5\%$ at $2\sigma$ with $H_0 = 70.65 \pm 0.66$ kms$^{-1}$Mpc$^{-1}$. Excluding the $H_0$ prior, SPT and Planck data constrain the shift in the Planck mass to be $<3\%$ at $2\sigma$ with a best-fit value of $0.04\%$, consistent with the $\Lambda$CDM limit. In this case $H_0 = 69.09^{+0.69}_{-0.68}$ kms$^{-1}$Mpc$^{-1}$, which is partially elevated by the dynamics of the scalar-field in the late universe. This differs from EDE models that prefer higher values of $H_0$ when high $\ell$ Planck TT data are excluded. We additionally constrain TPM using RSD data from BOSS DR 12 and cosmic shear, galaxy-galaxy lensing, and galaxy clustering data from DES Y1 finding both disfavor transitions close to recombination, but earlier Planck mass transitions are allowed.

We investigate how a GHz radio burst emitted near a magnetar propagates through its magnetosphere at radii $r=10^7$-$10^9$ cm. Bursts propagating near the magnetic equator behave as magnetohydrodynamic (MHD) waves if they have luminosity $L\gg 10^{40}$ erg/s. The waves develop plasma shocks in each oscillation and dissipate at $r\sim 3 \times 10^8 L_{42}^{-1/4}$ cm. GHz waves with lower $L$ or propagation directions closer to the magnetic axis do not obey MHD. Instead, they interact with individual particles, which requires a kinetic description. The kinetic interaction quickly accelerates particles to Lorentz factors $10^4$-$10^5$ at the expense of the wave energy, which again results in strong damping of the wave. In either regime of wave propagation, MHD or kinetic, the magnetosphere acts as a pillow absorbing the GHz burst and re-radiating the absorbed energy in X-rays. We conclude that a GHz source confined in the inner magnetosphere would be blocked by the outer magnetosphere at practically all relevant luminosities and viewing angles. This result constrains the origin of observed fast radio bursts (FRBs). We argue that observed FRBs come from magnetospheric explosions ejecting powerful outflows.

Ground level enhancements (GLEs) are extreme solar energetic particle (SEP) events that are of particular importance in space weather. In solar cycle 24, two GLEs were recorded on 2012 May 17 (GLE 71) and 2017 September 10 (GLE 72), respectively, by a range of advanced modern instruments. Here we conduct a comparative analysis of the two events by focusing on the effects of large-scale magnetic field configuration near active regions on particle acceleration and release. Although the active regions both located near the western limb, temporal variations of SEP intensities and energy spectra measured in-situ display different behaviors at early stages. By combining a potential field model, we find the CME in GLE 71 originated below the streamer belt, while in GLE 72 near the edge of the streamer belt. We reconstruct the CME shock fronts with an ellipsoid model based on nearly simultaneous coronagraph images from multi-viewpoints, and further derive the 3D shock geometry at the GLE onset. The highest-energy particles are primarily accelerated in the shock-streamer interaction regions, i.e., likely at the nose of the shock in GLE 71 and the eastern flank in GLE 72, due to quasi-perpendicular shock geometry and confinement of closed fields. Subsequently, they are released to the field lines connecting to near-Earth spacecraft when the shocks move through the streamer cusp region. This suggests that magnetic structures in the corona, especially shock-streamer interactions, may have played an important role in the acceleration and release of the highest-energy particles in the two events.

We have modeled the diffuse background at the Galactic Poles in the far-ultraviolet (FUV: 1536 \AA) and the near-ultraviolet (NUV: 2316 \AA). The background is well-fit using a single-scattering dust model with an offset representing the extragalactic light plus any other contribution to the diffuse background. We have found a dust albedo of 0.35 -- 0.40 (FUV) and 0.11 -- 0.19 in the NGP ($b > 70^{\circ}$) and 0.46 -- 0.56 (FUV) and 0.31 -- 0.33 (NUV) in the SGP ($b < 70^{\circ}$. The differences in the albedo may reflect changes in the dust-to-gas ratio over the sky or in the dust distribution. We find offsets at zero-reddening of 273 -- 286 and 553 -- 581 photons cm$^{-2}$ s$^{-1}$ sr$^{-1}$ \AA$^{-1}$ in the FUV and NUV, respectively, in the NGP with similar values in the SGP.

The recent detection of a stochastic background of gravitational waves in the nHz band by pulsar-timing array (PTA) experiments has shed new light on the formation and evolution of massive black hole binaries with masses $\sim 10^8$-$10^9 M_\odot$. The PTA data are consistent with a population of such binaries merging efficiently after the coalescence of their galactic hosts, and presenting masses slightly larger than previously expected. This momentous discovery calls for investigating the prospects of detecting the smaller ($\sim 10^5$-$10^7 M_\odot$) massive black hole binaries targeted by the Laser Interferometer Space Antenna (LISA). By using semi-analytic models for the formation and evolution of massive black hole binaries calibrated against the PTA results, we find that LISA will observe at least a dozen and up to thousands of black hole binaries during its mission duration. The minimum number of detections rises to $\sim 70$ if one excludes models that only marginally reproduce the quasar luminosity function at $z=6$. We also assess LISA's parameter estimation capabilities with state-of-the-art waveforms including higher modes and realistic instrumental response, and find that the masses, sky position, and distance will typically be estimated to within respectively 1%, 10 square degrees, and 10% for the detected systems.

Ground-based cosmic ray experiments detect cosmic ray mainly by measuring the longitudinal and lateral distribution of secondary particles produced in the extensive air shower (EAS). The EAS of cosmic ray in the knee energy region is simulated via CORSIKA software. Several simulation samples with different energy, composition and zenith angles were carried out to understand the longitudinal development of electron, muon and Cherenkov light in EAS. All the results presented were obtained assuming an observation plane at an altitude of 4400 m a.s.l. The differences of longitudinal development between electron and Cherenkov light were studied, and the reconstruction uncertainty of shower maximum for electron from Cherenkov light was estimated to be 10-15g/cm$^2$ for nuclei above 1 PeV. The performances of energy measurement and the composition discrimination ability based on longitudinal development were studied and compared with that from lateral distribution. It was found that number of electron per depth at its shower maximum has the smallest shower-to-shower fluctuations, but the shower-to-shower fluctuations of electron density measured at observation level was very close to it when the appropriate zenith angle was employed. The shower-to-shower fluctuations of shower maximum for electron is 50-55 g/cm$^2$ for proton, and 20-25 g/cm$^2$ for iron, but the composition discrimination ability between nuclei from muon density measured at observation level is much better than the shower maximum variable from longitudinal development. The hadronic model dependencies of the longitudinal development and lateral distribution were also discussed.

The second moment of the stellar velocity within the effective radius, denoted by $\sigma_{\rm e}^2$, is a crucial quantity in galaxy studies as it provides insight into galaxy properties and their mass distributions. However, large spectroscopic surveys typically do not measure $\sigma_{\rm e}$ directly, instead providing $\sigma_{\rm aper}$, the second moment of the stellar velocity within a fixed fiber aperture. In this paper, we derive an empirical aperture correction formula, given by $\sigma_{\rm aper}/\sigma_{\rm e}=(R_{\rm aper}/R_{\rm e})^{\alpha}$, using spatially resolved stellar kinematics extracted from approximately 10,000 Sloan Digital Sky Survey-Mapping Nearby Galaxies at Apache Point Observatory (SDSS-MaNGA) integral field unit observations. Our analysis reveals a strong dependence of $\alpha$ on the $r$-band absolute magnitude $M_{\rm r}$, $g-i$ color, and Sersic index $n_{\rm Ser}$, where $\alpha$ values are lower for brighter, redder galaxies with higher Sersic indices. Our results demonstrate that the aperture correction derived from previous literature on early-type galaxies cannot be applied to predict the aperture corrections for galaxies with intermediate Sersic indices. We provide a lookup table of $\alpha$ values for different galaxy types, with parameters in the ranges of $-18>M_{\rm r}>-24$, $0.4<g-i<1.6$, and $0<n_{\rm Ser}<8$. A Python script is provided to obtain the correction factors from the lookup table.

Anisotropic diffusion is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and the interplay of cosmic rays with the Galactic magnetic field. This diffusion term contributes to the highly stiff nature of the cosmic ray transport equation. To conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators (namely, Crank-Nicolson (CN)) have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods to treat the linear anisotropc diffusion equation in the presence of advection and time-independent and time-dependent sources. These methods allow us to take even larger step sizes that can substantially speed-up the simulations whilst generating highly accurate solutions. In or subsequent work, we will use these exponential solvers in the Picard code to study anisotropic cosmic ray diffusion and we will consider additional physical processes such as continuous momentum losses and reacceleration.

Magnetic fields are significant in the structure and evolution of stars. We present a comprehensive catalogue of 1784 known magnetic stars, detailing their identifications, HD numbers, precise locations, spectral types, and averaged quadratic effective magnetic fields among other important information. The group comprises 177 O-type stars, 551 B-type stars, 520 A-type stars, 91 F-type stars, 53 G-type stars, 61 K-type stars, 31 M-type stars, and an additional 300 stars whose spectral classification remains indeterminate. Our analysis examines the statistical properties of these magnetic stars. The relative integrated distribution function and number distribution function for all magnetic stars of the same spectral type can be effectively approximated using an exponential function of the averaged quadratic effective magnetic field. The analysis further reveals that A and B-type stars possess the strongest mean magnetic fields, indicating an easier detection of their magnetic fields.

Observations from space missions have allowed significant progress in many scientific domains due to the absence of atmospheric noise contributions and having uninterrupted data sets. In the context of asteroseismology, this has been extremely beneficial because many oscillation frequencies with small amplitudes, not observable from the ground, can be detected. One example of this success is the large number of hybrid delta Sct-gamma Dor stars discovered. These stars have radial and non-radial p- and g-modes simultaneously excited to an observable level allowing us to probe both the external and near-to-core layers of the star. We analyse the light curve of hybrid delta Sct-gamma Dor star CoRoT ID 102314644 and characterise its frequency spectrum. We detected 29 gamma Dor type frequencies in the range [0.32-3.66] cycles per day (c/d) and a series of 6 equidistant periods with a mean period spacing of DeltaPi=1612 s. In the delta Sct domain we found 38 frequencies in the range 8.63-24.73 c/d and a quintuplet centred on the frequency p_1=11.39 c/d and derived a possible rotational period of 3.06 d. The frequency analysis of this object suggests the presence of spots at the stellar surface, nevertheless we could not dismiss the possibility of a binary system. The initial modelling of the frequency data along with external constraints has allowed us to refine its astrophysical parameters giving a mass of approximately 1.75 solar masses, a radius of 2.48 solar radii and an age of 1241 Myr. The observed period spacing, a p-mode quintuplet, the possible rotation period and the analysis of the individual frequencies provide important input constraints for the understanding of different transport phenomena in A-F-type stars.[abridged]

Theoretical physical-chemical models for the formation of planetary systems depend on data quality for the Sun's composition, that of stars in the solar neighbourhood, and of the estimated "pristine" compositions for stellar systems. The effective scatter and the observational uncertainties of elements within a few hundred parsecs from the Sun, even for the most abundant metals like carbon, oxygen and silicon, are still controversial. Here we analyse the stellar production and the chemical evolution of key elements that underpin the formation of rocky (C, O, Mg, Si) and gas/ice giant planets (C, N, O, S). We calculate 198 galactic chemical evolution (GCE) models of the solar neighbourhood to analyse the impact of different sets of stellar yields, of the upper mass limit for massive stars contributing to GCE ($M_{\rm up}$) and of supernovae from massive-star progenitors which do not eject the bulk of the iron-peak elements (faint supernovae). Even considering the GCE variation produced via different sets of stellar yields, the observed dispersion of elements reported for stars in the Milky Way disk is not reproduced. Among others, the observed range of super-solar [Mg/Si] ratios, sub-solar [S/N], and the dispersion of up to 0.5 dex for [S/Si] challenge our models. The impact of varying $M_{\rm up}$ depends on the adopted supernova yields. Thus, observations do not provide a constraint on the M$_{\rm up}$ parametrization. When including the impact of faint supernova models in GCE calculations, elemental ratios vary by up to 0.1-0.2 dex in the Milky Way disk; this modification better reproduces observations.

We investigate the alignment of galaxy and halo orientations using the TNG300-1 hydrodynamical simulation. Our analysis reveals that the distribution of the 2D misalignment angle $\theta_{\rm{2D}}$ can be well described by a truncated shifted exponential (TSE) distribution with only {\textit{one}} free parameter across different redshifts and galaxy/halo properties. We demonstrate that the galaxy-ellipticity (GI) correlations of galaxies can be reproduced by perturbing halo orientations with the obtained $\theta_{\rm{2D}}$ distribution, with only a small bias ($<3^{\circ}$) possibly arising from unaccounted couplings between $\theta_{\rm{2D}}$ and other factors. We find that both the 2D and 3D misalignment angles $\theta_{\rm{2D}}$ and $\theta_{\rm{3D}}$ decrease with ex situ stellar mass fraction $F_{\rm{acc}}$, halo mass $M_{\rm{vir}}$ and stellar mass $M_{*}$, while increasing with disk-to-total stellar mass fraction $F_{\rm{disk}}$ and redshift. These dependences are in good agreement with our recent observational study based on the BOSS galaxy samples. Our results suggest that $F_{\rm{acc}}$ is a key factor in determining the galaxy-halo alignment. Grouping galaxies by $F_{\rm{acc}}$ nearly eliminates the dependence of $\theta_{\rm{3D}}$ on $M_{\rm{vir}}$ for all three principle axes, and also reduces the redshift dependence. For $\theta_{\rm{2D}}$, we find a more significant redshift dependence than for $\theta_{\rm{3D}}$ even after controlling $F_{\rm{acc}}$, which may be attributed to the evolution of galaxy and halo shapes. Our findings present a valuable model for observational studies and enhance our understanding of galaxy-halo alignment.

We combine mid-infrared diagnostics obtained from integral-field unit observations taken with MIRI/MRS on JWST with cold molecular gas information derived from ALMA observations of CO(1-0) emission to investigate the star formation rate and efficiency within the central $\sim$ 1.5 kpc$\times$1.3 kpc region of the Seyfert 1 galaxy NGC 7469 on $\sim$ 100 pc scales. The active nucleus leaves a notable imprint on its immediate surroundings by elevating the temperature of the warm molecular gas, driving an ionized gas outflow on sub-kpc scales, and selectively destroying small dust grains. These effects, nevertheless, have relatively little impact on the cold circumnuclear medium or its ability to form stars. Most of the star formation in NGC 7469 is confined to a clumpy starburst ring, but the star formation efficiency remains quite elevated even for the nuclear region that is most affected by the active nucleus.

In this paper, we seek optimal solutions for a transfer from a parking orbit around the Moon to a halo orbit around $L_2$ of the Earth-Moon system, by applying a single maneuver and exploiting the stable invariant manifold of the hyperbolic parking solution at arrival. For that, we propose an optimization problem considering as variables both the orbital characteristics of a parking solution around the Moon, (namely, its Keplerian elements) and the characteristics of a transfer trajectory guided by the stable manifold of the arrival Halo orbit. The problem is solved by a nonlinear programming method (NLP), aiming to minimize the cost of $\Delta V$ to perform a single maneuver transfer, within the framework of the Earth-Moon system of the circular restricted three-body problem. Results with low $\Delta V$ and suitable time of flight show the feasibility of this kind of transfer for a Cubesat.

Radiative turbulent mixing layers are expected to form pervasively at the phase boundaries in multiphase astrophysical systems. This inherently small scale structure is dynamically crucial because it directly regulates the mass, momentum and energy exchanges between adjacent phases. Previous studies on hydrodynamic turbulent mixing layers have revealed the interactions between cold and hot phases in the context of the circumgalactic medium, offering important insight into the fate of cold clouds traveling through hot galactic winds. However, the role of magnetic field has only been sparsely investigated. We perform a series of 3D magnetohydrodynamics (MHD) simulations of such mixing layers in the presence of weak to modest background magnetic field. We find that due to field amplification, even relatively weak background magnetic fields can significantly reduce the surface brightness and inflow velocity of the hot gas in the mixing layer. This reduction is attributed to a combination of magnetic pressure support and direct suppression of turbulent mixing, both of which alter the phase structures. Our results are largely independent of thermal conduction and converged with resolution, offering insights on the survival of cold gas in multiphase systems.

The high-energy diffuse gamma-ray emission and neutrino emission are expected from the Galactic plane, generated by hadronuclear interactions between cosmic rays (CR) and interstellar medium (ISM). Therefore, measurements of these diffuse emissions will provide important clues on the origin and nature of Galactic CRs. Comparing the latest observations of LHAASO and IceCube on the diffuse Galactic gamma-ray and neutrino emissions respectively, we suggest that the diffuse gamma-ray emission at multi-TeV energies contains a considerable contribution of a leptonic component. By modelling the gamma-ray halos powered by middle-aged pulsars in our Galaxy with taking into account the magnetic field configuration and the interstellar radiation field in the Galaxy, we demonstrate that the collective contribution of pulsar halos can account for the excess in the measured diffuse gamma-ray emission with respect to the predicted flux from CR-ISM interactions. The resulting one-dimensional profile along the Galactic longitude is also consistent with the observation.

In order to advance our understanding of the dynamic interactions between coronal mass ejections (CMEs) and the magnetized solar wind, we investigate the impact of magnetic erosion on the well-known aerodynamic drag force acting on CMEs traveling faster than the ambient solar wind. In particular, we start by generating empirical relationships for the basic physical parameters of CMEs that conserve their mass and magnetic flux. Furthermore, we examine the impact of the virtual mass on the equation of motion by studying a variable-mass system. We next implement magnetic reconnection into CME propagation, which erodes part of the CME magnetic flux and outer-shell mass, on the drag acting on CMEs, and we determine its impact on their time and speed of arrival at 1 AU. Depending on the strength of the magnetic erosion, the leading edge of the magnetic structure can reach near-Earth space up to $\approx$ three hours later, compared to the non-eroded case. Therefore, magnetic erosion may have a significant impact on the propagation of fast CMEs and on predictions of their arrivals at 1 AU. Finally, the modeling indicates that eroded CMEs may experience a significant mass decrease. Since such a decrease is not observed in the corona, the initiation distance of erosion may lie beyond the field-of-view of coronagraphs (i.e. 30 $\mathrm{R_{\odot}}$).

Cepheids have long been used as standard candles to determine distances around the Milky Way and to nearby galaxies. A discrepancy still remains for Hubble Constant determinations using Cepheids vs. the cosmic microwave background or calibrations to the tip of the red-giant branch. Therefore, refinement of Cepheid period-luminosity relations continues to be an active topic of research. Cepheids are also important laboratories for testing stellar physics. This paper explores outstanding questions in Cepheid evolution and pulsation modeling. We examine the discrepancy between Cepheid masses determined from pulsation properties and binary orbital dynamics and those determined using stellar evolution models. We review attempts to resolve the discrepancy by including rotation, convective overshooting, and mass loss. We review the impact of uncertainties in nuclear reaction rates on Cepheid evolution and the extent of blue loops in the Hertzsprung-Russell diagram. We consider implications for Cepheids of stellar opacity revisions suggested in light of findings for the Sun and other types of variable stars. We apply the 1-D open-source MESA stellar evolution code and the MESA radial stellar pulsation (RSP) nonlinear hydrodynamics code to investigate changes in input physics for Cepheid models. We touch on progress in 2-D and 3-D stellar modeling applied to Cepheids. Additional areas in which Cepheid models are being tested against observations include: predicting the edges of the Cepheid pulsation instability strip; predicting period-change rates and implications for instability strip crossings; explaining period and amplitude modulations and periodicities that may be non-radial pulsation modes; discovering what can be learned from Cepheid observations in X-ray, ultraviolet, and radio wavelengths. We also show a few examples of Cepheid light curves from NASA TESS photometry.

NEID is a high-resolution red-optical precision radial velocity (RV) spectrograph recently commissioned at the WIYN 3.5 m telescope at Kitt Peak National Observatory, Arizona, USA. NEID has an extremely stable environmental control system, and spans a wavelength range of 380 to 930 nm with two observing modes: a High Resolution (HR) mode at R $\sim$ 112,000 for maximum RV precision, and a High Efficiency (HE) mode at R $\sim$ 72,000 for faint targets. In this manuscript we present a detailed description of the components of NEID's optical fiber feed, which include the instrument, exposure meter, calibration system, and telescope fibers. Many parts of the optical fiber feed can lead to uncalibratable RV errors, which cannot be corrected for using a stable wavelength reference source. We show how these errors directly cascade down to performance requirements on the fiber feed and the scrambling system. We detail the design, assembly, and testing of each component. Designed and built from the bottom-up with a single-visit instrument precision requirement of 27 $\textrm{cm~s}^{-1}$, close attention was paid to the error contribution from each NEID subsystem. Finally, we include the lab and on-sky tests performed during instrument commissioning to test the illumination stability, and discuss the path to achieving the instrumental stability required to search for a true Earth twin around a Solar-type star.

Nomadic worlds, i.e., objects not gravitationally bound to any star(s), are of great interest to planetary science and astrobiology. They have garnered attention recently due to constraints derived from microlensing surveys and the recent discovery of interstellar planetesimals. In this paper, we roughly estimate the prevalence of nomadic worlds with radii of $100\,\mathrm{km} \lesssim R \lesssim 10^4\,\mathrm{km}$. The cumulative number density $n_>\left(>R\right)$ appears to follow a heuristic power law given by $n_> \propto R^{-3}$. Therefore, smaller objects are probably much more numerous than larger rocky nomadic planets, and statistically more likely to have members relatively close to the inner Solar system. Our results suggest that tens to hundreds of planet-sized nomadic worlds might populate the spherical volume centered on Earth and circumscribed by Proxima Centauri, and may thus comprise closer interstellar targets than any planets bound to stars. For the first time, we systematically analyze the feasibility of exploring these unbounded objects via deep space missions. We investigate what near-future propulsion systems could allow us to reach nomadic worlds of radius $> R$ in a $50$-year flight timescale. Objects with $R \sim 100$ km are within the purview of multiple propulsion methods such as electric sails, laser electric propulsion, and solar sails. In contrast, nomadic worlds with $R \gtrsim 1000$ km are accessible by laser sails (and perhaps nuclear fusion), thereby underscoring their vast potential for deep space exploration.

One of the key open questions in Cosmology is the nature of the sources that completed the cosmological hydrogen Reionization at z~5.2. High-z primeval galaxies have been long considered the main drivers for Reionization, with a minor role played by high-z AGN. However, in order to confirm this scenario, it is fundamental to measure the photo-ionization rate produced by active SMBHs close to the epoch of Reionization. Given the pivotal role played by spectroscopically complete observations of high-z QSOs, in this paper we present the first results of the RUBICON (Reionizing the Universe with BrIght COsmological Nuclei) survey. It consists of a color selected sample of bona-fide z~5 QSO candidates from the Hyper Suprime-Cam Subaru Strategic Survey. Our QSO candidates have been validated both by photometric redshifts based on SED fitting and by spectroscopic redshifts, confirming that they lie at 4.5<z_spec<5.2. A relatively large space density of QSOs (Phi~1.4x10^-8 cMpc^-3) is thus confirmed at z~5 and M1450~-27, consistent with a pure density evolution of the AGN luminosity function from z=4 to z=5, with a mild density evolution rate of 0.25 dex. This indicates that AGN could play a non-negligible role in the cosmic Reionization. The Rubicon of Reionization has been crossed.

The fragmentation of gas to form stars in molecular clouds is intrinsically linked to the turbulence within them. These internal motions are set at the birth of the cloud and may vary with galactic environment and as the cloud evolves. In this paper, we introduce a new suite of 15 high-resolution molecular cloud simulations using the moving mesh code AREPO, to investigate the role of different decaying turbulent modes (mixed, compressive and solenoidal) and Virial ratios on the evolution of a $10^4\mathrm{M}_{\odot}$ molecular cloud. We find that diffuse regions maintain a strong relic of the initial turbulent mode, whereas the initial gravitational potential dominates dense regions. Solenoidal seeded models thus give rise to a diffuse cloud with filament-like morphology, and an excess of brown dwarf mass fragments. Compressive seeded models have an early onset of star-formation, cluster-like morphologies and a higher accretion rate, along with overbound clouds, compared to other simulations. Filaments identified using DisPerSE, and analyzed through a new Python toolkit we develop and make publicly available with this work called FIESTA, show no clear trend in lengths, masses and densities between initial turbulent modes. Overbound clouds, however, produce more filaments and thus have more mass in filaments. The hubs formed by converging filaments are found to favour star-formation, with surprisingly similar mass distributions independent of the number of filaments connecting the hub.

Relativistic electrons are an essential component in many astrophysical sources, and their radiation may dominate the high-energy bands. Inverse Compton (IC) emission is the radiation mechanism that plays the most important role in these bands. The basic properties of IC, such as the total and differential cross sections, have long been studied; the properties of the IC emission depend strongly not only on the emitting electron distribution but also on the properties of the target photons. This complicates the phenomenological studies of sources, where target photons are supplied from a broad radiation component. We study the spectral properties of IC emission generated by a power-law distribution of electrons on a power-law distribution of target photons. We approximate the resulting spectrum by a broken-power-law distribution and show that there can be up to three physically motivated spectral breaks. If the target photon spectrum extends to sufficiently low energies, $\varepsilon_{\mathrm{min}}< m_e^2c^4/E_{\mathrm{max}}$ ($m_e$ and $c$ are electron mass and speed of light, respectively; $\varepsilon_{\mathrm{min}}$ and $E_{\mathrm{max}}$ are the minimum/maximum energies of target photons and electrons, respectively), then the high energy part of the IC component has a spectral slope typical for the Thomson regime with an abrupt cutoff close to $E_{\mathrm{max}}$. The spectra typical for the Klein-Nishina regime are formed above $m_e^2c^4/\varepsilon_{\mathrm{min}}$. If the spectrum of target photons features a cooling break, i.e., a change of the photon index by $0.5$ at $\varepsilon_{\mathrm{br}}$, then the transition to the Klein-Nishina regime proceeds through an intermediate change of the photon index by $0.5$ at $m_e^2c^4/\varepsilon_{\mathrm{br}}$.

Gravitational lensing provides unique insights into astrophysics and cosmology, including the determination of galaxy mass profiles and constraining cosmological parameters. We present spectroscopic confirmation and lens modeling of the strong lensing system DESI-253.2534+26.8843, discovered in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys data. This system consists of a massive elliptical galaxy surrounded by four blue images forming an Einstein Cross pattern. We obtained spectroscopic observations of this system using the Multi Unit Spectroscopic Explorer (MUSE) on ESO's Very Large Telescope (VLT) and confirmed its lensing nature. The main lens, which is the elliptical galaxy, has a redshift of $z_{L1} = 0.636\pm 0.001$, while the spectra of the background source images are typical of a starburst galaxy and have a redshift of $z_s = 2.597 \pm 0.001$. Additionally, we identified a faint galaxy foreground of one of the lensed images, with a redshift of $z_{L2} = 0.386$. We employed the GIGA-Lens modeling code to characterize this system and determined the Einstein radius of the main lens to be $\theta_{E} =2.520{''}_{-0.031}^{+0.032}$, which corresponds to a velocity dispersion of $\sigma$ = 379 $\pm$ 2 km s$^{-1}$. Our study contributes to a growing catalog of this rare kind of strong lensing systems and demonstrates the effectiveness of spectroscopic integral field unit observations and advanced modeling techniques in understanding the properties of these systems.

Tremendous information is hidden in the light curve of a gamma-ray burst (GRB). Based on CGRO/BATSE data, Hakkila (2021) found a majority of GRBs can be characterized by a smooth, single-peaked component superposed with a temporally symmetrical residual structure, i.e., a mirror feature for the fast varying component. In this study, we conduct a similar analysis on the same data, as well as on Fermi/GBM data. We got a similar conclusion that most GRBs have this symmetrical fast varying component. Further more, we chose an alternative model to characterize the smooth component and used a three-parameter model to identify the residual, i.e., the fast component. By choosing 226 BATSE GRBs based on a few criteria, we checked the time symmetrical feature and time translational feature for the fast components and found the ratio is roughly 1:1. We propose that both features could come from the structure of the ejected shells. Future SKA might be able to observe the early radio emission from the collision of the shells.

Early results from the James Webb Space Telescope (JWST) observations have hinted at two traces beyond the standard cosmological framework. One is the extraordinarily high stellar masses and their density at $z=7.5\sim9.1$, another is the unexpected abundance of ultra-violet (UV) bright galaxies at $z\ge10$. Nevertheless, both pieces of evidence are not statistically robust, yet. In this work, we construct rest-frame UV luminosity functions (LFs) based on a general formation model for these high-redshift galaxy candidates, since UV LFs always carry the information of stellar formation efficiency (SFE), initial mass function (IMF), dust attenuation, and other crucial elements for galaxy evolution. By updating the massive galaxies candidates with spectroscopic observations and exploring the parameter space of SFE, we are able to reasonably explain the cumulative stellar mass density within the redshift range of $7.5\sim9.1$, with only one galaxy exhibiting unusual characteristics. We also reveal a potential non-monotonic trend of SFE with the increasing redshift. At higher redshift ($z\sim13$), bright UV LFs can be well fitted with non-dust attenuation or Top-heavy IMF for Population III stars. The Population III star scenario can also naturally account for the possible dip of SFE at $z\sim9$.

We present deep Hubble Space Telescope images taken to examine the ejecta from the DART spacecraft impact into asteroid Dimorphos. The images reveal an extensive population of co-moving boulders, the largest of which is about 7 m in diameter (geometric albedo 0.15 assumed). Measurements of 37 boulders show a mean sky-plane velocity dispersion of 0.30+/-0.03 m/s, only slightly larger than the 0.24 m/s gravitational escape velocity from the Didymos/Dimorphos binary system. The total boulder mass, 5e6 kg (density 2200 kg/m3 assumed), corresponds to about 0.1 percent of the mass of Dimorphos and the boulders collectively carry about 3e-5 of the kinetic energy delivered by the DART spacecraft impact. The sky-plane distribution of the boulders is asymmetric, consistent with impact into an inhomogeneous, likely rubble-pile, body. Surface boulder counts on Didymos show that the observed boulder swarm could be ejected from as little as 2 percent of the surface of Dimorphos (for example a circular crater at the impact point about 50 m in diameter). The large, slow-moving boulders are potential targets to be investigated in-situ by the upcoming ESA HERA mission.

We present an incremental version (4FGL-DR4, for Data Release 4) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first 14 years of science data in the energy range from 50 MeV to 1 TeV, it uses the same analysis methods as the 4FGL-DR3 catalog did for 12 years of data, with only a few improvements. The spectral parameters, spectral energy distributions, light curves and associations are updated for all sources. We add four new extended sources and modify two existing ones. Among the 6658 4FGL-DR3 sources, we delete 14 and change the localization of 10, while 26 are newly associated and two associations were changed. We add 546 point sources, among which 8 are considered identified and 228 have a plausible counterpart at other wavelengths. Most are just above the detection threshold, and 14 are transient sources below the detection threshold that can affect the light curves of nearby sources.

We perform a structural decomposition of galaxies identified in three cosmological hydrodynamical simulations by applying Gaussian Mixture Models (GMMs) to the kinematics of their stellar particles. We study the resulting disc, bulge, and intra-halo light (IHL) components of galaxies whose host dark matter haloes have virial masses in the range $M_{200}=10^{11}$-- $10^{15}\,{\rm M_\odot}$. Our decomposition technique isolates galactic discs whose mass fractions, $f_{\rm disc}$, correlate strongly with common alternative morphology indicators; for example, $f_{\rm disc}$ is approximately equal to $\kappa_{{\rm co}}$, the fraction of stellar kinetic energy in co-rotation. The primary aim of our study, however, is to characterise the IHL of galaxies in a consistent manner and over a broad mass range, and to analyse its properties from the scale of galactic stellar haloes up to the intra-cluster light. Our results imply that the IHL fraction, $f_{\rm IHL}$, has appreciable scatter and is strongly correlated with galaxy morphology: at fixed stellar mass, the IHL of disc galaxies is typically older and less massive than that of spheroids. Above $M_{200}\approx 10^{13}\,{\rm M_\odot}$, we find, on average, $f_{\rm IHL}\approx 0.45$, albeit with considerable scatter. The transition radius beyond which the IHL dominates the stellar mass of a galaxy is roughly $30\,{\rm kpc}$ for $M_{200}\lesssim 10^{12.8}\,{\rm M_\odot}$, but increases strongly towards higher masses. However, we find that no alternative IHL definitions -- whether based on the ex-situ stellar fraction, or the stellar mass outside a spherical aperture -- reproduce our dynamically-defined IHL fractions.

Recently, LHAASO reported the detection of brightest-of-all-time GRB 221009A, revealing the early onset of a TeV afterglow. However, there is no evidence of afterglow emission at such early time at other wavelengths. Here we report the discovery of a hidden afterglow component during the prompt emission phase with Fermi Gamma-Ray Burst Monitor (GBM) observations. We analyze the spectral evolution of the X-ray/$\gamma$-ray emission of GRB 221009A measured by GBM during the dips of two prompt emission pulses (i.e., intervals $T_{0}+[300-328]\rm~s$ and $T_{0}+[338-378]\rm~s$, where $T_0$ is the GBM trigger time). We find that the spectra at the dips transit from the Band function to a power-law function, indicating a transition from the prompt emission to the afterglow. After $\sim T_{0}+ 660 \rm~s$, the spectrum is well described by a power-law function and the afterglow becomes dominant. Remarkably, the underlying afterglow emission at the dips smoothly connect with the afterglow after $\sim T_{0}+ 660 \rm~s$. The entire afterglow emission measured by GBM can be fitted by a power-law function $F\sim t^{-0.95\pm0.05}$, where $t$ is the time since the first main pulse at $T^*=T_0+226~{\rm s}$, consistent with the TeV afterglow decay measured by LHAASO. The start time of this power-law decay indicates that the afterglow peak of GRB 221009A should be earlier than $T_{0}+300 \rm ~s$. We also test the possible presence of a jet break in the early afterglow light curve, finding that both the jet break model and single power-law decay model are consistent with the GBM data. The two models can not be distinguished with the GBM data alone because the inferred jet break time is quite close to the end of GBM observations.

Pulsar timing arrays (PTAs) provide a way to detect gravitational waves (GWs) at nanohertz frequencies. To ensure the detection of GWs, observational data must exhibit the Hellings-Downs angular correlation. It is also known that PTAs can probe ultralight dark matter. In this paper, we consider possible contamination of the Hellings-Downs angular correlation by the ultralight dark matter. We find that ultralight vector dark matter can give rise to the deformation of the Hellings-Downs correlation curve. Thus, the Hellings-Downs correlation curve could contain information on ultralight dark matter with a spin.

The heating mechanisms of solar white-light flares remain unclear. We present an X1.0 white-light flare on 2022 October 2 (SOL2022-10-02T20:25) observed by the Chinese \ha\ Solar Explorer (CHASE) that provides two-dimensional spectra in the visible light for the full solar disk with a seeing-free condition. The flare shows a prominent enhancement of $\sim$40\% in the photospheric \fe\ line at 6569.2 \AA, and the nearby continuum also exhibits a maximum enhancement of $\sim$40\%. For the continuum near the \fe\ line at 6173 \AA\ from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO), it is enhanced up to $\sim$20\%. At the white-light kernels, the \fe\ line at 6569.2 \AA\ has a symmetric Gaussian profile that is still in absorption and the H$\alpha$ line at 6562.8 \AA\ displays a very broad emission profile with a central reversal plus a red or blue asymmetry. The white-light kernels are co-spatial with the microwave footpoint sources observed by the Expanded Owens Valley Solar Array (EOVSA) and the time profile of the white-light emission matches that of the hard X-ray emission above 30 keV from the Gamma-ray Burst Monitor (GBM) on Fermi. These facts indicate that the white-light emission is qualitatively related to a nonthermal electron beam. We also perform a radiative hydrodynamic simulation with the electron beam parameters constrained by the hard X-ray observations from Fermi/GBM. The result reveals that the white-light enhancement cannot be well explained by a pure electron-beam heating together with its induced radiative backwarming but may need additional heating sources such as \alfven\ waves.

For binary systems with an unseen primary and a luminous secondary, the astrometric wobble of the secondary could be used to study the primary. With Gaia, it is possible to measure the mass of the black hole or neutron star with a luminous companion (hereafter BH/NS-LC). Our aim is to provide a method for predicting Gaia's ability in measuring the mass of BH/NS-LCs. We also tried to estimate the number of solvable BH/NS-LCs using Gaia. We used a realistic Markov chain Monte Carlo simulation of mock Gaia observations to obtain a relation between the uncertainty of mass measurement of the primary in BH/NS-LCs with the observable variables of the secondary astrometric orbit. Furthermore, we used the MOBSE code to evolve a Galactic BH/NS-LC sample with a combined Milky Way model. Our relation is applied to this sample to estimate the number of solvable BH/NS-LCs. We derived a good relation between the mass uncertainty and the binary parameters. For the first time, we show the quantitive influence of the period P, inclination i, eccentricity e, and ecliptic latitude $\beta$ to the mass measurement. Our results suggest that $48^{+7}_{-7}$ BH-LCs and $102^{+11}_{10}$ NS-LCs are solvable during a 5 yr Gaia mission. We also give the distribution of the distance and apparent magnitude of the Gaia solvable BH/NS-LCs. This solvable sample would be increased by additional spectroscopic data or a prolonged Gaia mission. The mass uncertainty relation could be used in future simulations of BH/NS-LCs observed by Gaia. The prediction of the solvable BH/NS-LCs is not only influenced by the process in generating the Galactic BH/NS-LC sample, but is also affected by our uncertainty relation. In particular, the relations of parameters such as $[P, e, i, \beta]$ are very useful to correct the selection effect in the statistic results of the future BH/NS-LC sample observed by Gaia.

The presence of pulsating stars in eclipsing binary systems (EBs) makes these objects significant since they allow us to investigate the stellar interior structure and evolution. Different types of pulsating stars could be found in EBs such as Delta Scuti variables. Delta Scuti stars in EBs have been known for decades and the increasing number of such systems is important for understanding pulsational structure. Hence, in this study, a research was carried out on the southern TESS field to discover new Delta Scuti stars in EBs. We produced an algorithm to search for detached and semi-detached EBs considering three steps; the orbital period (P$_{orb}$)'s harmonics in the Fourier spectrum, skewness of the light curves, and classification of \textsc{UPSILON} program. If two of these steps classify a system as an EB, the algorithm also identifies it as an EB. The TESS pixel files of targets were also analyzed to see whether the fluxes are contaminated by other systems. No contamination was found. We researched the existence of pulsation through EBs with a visual inspection. To confirm Delta Scuti-type oscillations, the binary variation was removed from the light curve, and residuals were analyzed. Consequently, we identified 42 Delta Scuti candidates in EBs. The P$_{orb}$, $L$, and M$_{V}$ of systems were calculated. Their positions on the H-R diagram and the known orbital-pulsation period relationship were analyzed. We also examined our targets to find if any of them show frequency modulation with the orbital period and discovered one candidate of tidally tilted pulsators.

We analyze a Neutron Star Interior Composition Explorer (NICER) observation of the black hole X-ray binary MAXI J1820+070 during a transition from type-C to type-B quasi-periodic oscillations (QPOs). We find that below ~2 keV, for the type-B QPOs the rms amplitude is lower and the magnitude of the phase lags is larger than for the type-C QPOs. Above that energy, the rms and phase-lag spectra of the type-B and type-C QPOs are consistent with being the same. We perform a joint fit of the time-averaged spectra of the source, and the rms and phase-lag spectra of the QPOs with the time-dependent Comptonization model vkompth to study the geometry of the corona during the transition. We find that the data can be well-fitted with a model consisting of a small and a large corona that are physically connected. The sizes of the small and large coronae increase gradually during the type-C QPO phase whereas they decrease abruptly at the transition to type-B QPO. At the same time, the inner radius of the disc moves inward at the QPO transition. Combined with simultaneous radio observations showing that discrete jet ejections happen around the time of the QPO transition, we propose that a corona that expands horizontally during the type-C QPO phase, from ~10^{4} km (~800 Rg) to ~10^{5} km (~8000 Rg) overlying the accretion disc, transforms into a vertical jet-like corona extending over ~10^{4} km (~800 Rg) during the type-B QPO phase.

Using observations with e-MERLIN and the VLA, together with archival data from ALMA, we obtain high-resolution radio images of two binary YSOs: L1551 IRS 5 and L1551 NE, covering a wide range of frequencies from 5 - 336 GHz, and resolving emission from the radio jet on scales of only ~15 au. By comparing these observations to those from a previous epoch, it is shown that there is a high degree of variability in the free-free emission from the jets of these sources. In particular, the northern component of L1551 IRS 5 shows a remarkable decline in flux density of a factor of ~5, suggesting that the free-free emission of this source has almost disappeared. By fitting the spectra of the sources, the ionised mass-loss rates of the jets are derived and it is shown that there is significant variability of up to a factor of ~6 on timescales of ~20 years. Using radiative transfer modelling, we also obtained a model image for the jet of the southern component of L1551 IRS 5 to help study the inner region of the ionised high-density jet. The findings favour the X-wind model launched from a very small innermost region.

The past 10 years have seen tremendous progress in our capability of testing General Relativity in the strong field regime with black hole observations. 10 years ago, the theory of General Relativity was almost completely unexplored in the strong field regime. Today, we have gravitational wave data of the coalescence of stellar-mass black holes, radio images of the supermassive black holes SgrA$^*$ and M87$^*$, and high-quality X-ray data of stellar-mass black holes in X-ray binaries and supermassive black holes in active galactic nuclei. In this manuscript, we will review current efforts to test General Relativity with black hole X-ray data and we will provide a detailed description of the public codes available on ABHModels.

Early data from the James Webb Space Telescope (JWST) has uncovered the existence of a surprisingly abundant population of very massive galaxies at extremely high redshift, which are hard to accommodate within the standard $\Lambda$CDM cosmology. We explore whether the JWST observations may be pointing towards more complex dynamics in the dark energy (DE) sector. Motivated by the ubiquity of anti-de Sitter vacua in string theory, we consider a string-inspired scenario where the DE sector consists of a negative cosmological constant (nCC) and a evolving component with positive energy density on top, whose equation of state is allowed to cross the phantom divide. We show that such a scenario can drastically alter the growth of structure compared to $\Lambda$CDM, and accommodate the otherwise puzzling JWST observations if the dynamical component evolves from the quintessence-like regime in the past to the phantom regime today: in particular, we demonstrate that the presence of a nCC (which requires a higher density for the evolving component) plays a crucial role in enhancing the predicted cumulative comoving stellar mass density. Our work reinforces the enormous potential held by observations of the abundance of high-$z$ galaxies in probing cosmological models and new fundamental physics, including string-inspired ingredients.

The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy neutrinos, cosmic rays, and gamma rays, with energies above 100 PeV. GRAND targets the radio signals emitted by extensive air showers induced by the interaction of ultra-high-energy particles in the atmosphere, using an array of 200,000 radio antennas split into sub-arrays deployed worldwide. GRANDProto13 (GP13) is a 13-antenna demonstrator array deployed in February 2023 in the Gansu province of China, as a precursor for GRANDProto300, which will validate the detection principle of the GRAND experiment. Its goal is to measure the radio background present at the site, validate the design of the detection units and develop an autonomous radio trigger for air showers. We will describe GP13 and its operation, and show preliminary results on noise monitoring.

Numerous exoplanets with masses ranging from Earth to Neptune and radii larger than Earth have been found through observations. These planets possess atmospheres that range in mass fractions from 1% to 30%, reflecting the diversity of atmospheric mass fractions. Such diversities are supposed to be caused by differences in the formation processes or evolution. Here we consider head-on giant impacts onto planets causing atmosphere losses in the later stage of their formation. We perform smoothed particle hydrodynamic simulations to study the impact-induced atmosphere loss of young super-Earths with 10%-30% initial atmospheric mass fractions. We find that the kinetic energy of the escaping atmosphere is almost proportional to the sum of the kinetic impact energy and self-gravitational energy released from the merged core. We derive the relationship between the kinetic impact energy and the escaping atmosphere mass. The giant impact events for planets of comparable masses are required in the final stage of the popular scenario of rocky planet formation. We show it results in a significant loss of the atmosphere, if the impact is a head-on collision with comparable masses. This latter fact provides a constraint on the formation scenario of rocky planets with substantial atmospheres.

Context: High-spatial resolution Atacama Large Millimeter/submillimeter Array (ALMA) data have revealed a plethora of substructures in protoplanetary disks. Some of those features are thought to trace the formation of embedded planets. One example is the gas and dust that accumulated in the co-orbital Lagrangian regions $L_4$/$L_5$, which were tentatively detected in recent years and might be the pristine material for the formation of Trojan bodies. Aims: This work is part of the TROY project, whose ultimate goal is to find robust evidence of exotrojan bodies and study their implications in the exoplanet field. Here, we focus on the early stages of the formation of these bodies by inspecting the iconic system PDS 70, the only confirmed planetary system in formation. Methods: We reanalyzed archival high-angular resolution Band 7 ALMA observations from PDS 70 by doing an independent imaging process to look for emission in the Lagrangian regions of the two detected gas giant protoplanets, PDS 70 b and c. We then projected the orbital paths and visually inspected emission features at the regions around the $L_4$/$L_5$ locations as defined by $\pm$ 60$^{\circ}$ in azimuth from the planet position. Results: We found emission at a $\sim$4-$\sigma$ level ($\sim$6-$\sigma$ when correcting from a cleaning effect) at the position of the $L_{5}$ region of PDS 70 b. This emission corresponds to a dust mass in a range of 0.03- 2 M$_{Moon}$, which potentially accumulated in this gravitational well. Conclusions: The tentative detection of the co-orbital dust trap that we report requires additional observations to be confirmed. We predict that we could detect the co-orbital motion of PDS 70 b and the dust presumably associated with $L_5$ by observing again with the same sensitivity and angular resolution as early as February 2026.

One key challenge in detecting 21 cm cosmological signal at z > 6 is to separate the cosmological signal from foreground emission. This can be studied in a power spectrum space where the foreground is confined to low delay modes whereas the cosmological signal can spread out to high delay modes. When there is a calibration error, however, chromaticity of gain errors propagates to the power spectrum estimate and contaminates the modes for cosmological detection. The Hydrogen Epoch of Reionization Array (HERA) employs a high-precision calibration scheme using redundancy in measurements. In this study, we focus on the gain errors induced by nonredundancies arising from feed offset relative to the HERA's 14 meter parabolic dish element, and investigate how to mitigate the chromatic gain errors using three different methods: restricting baseline lengths for calibration, smoothing the antenna gains, and applying a temporal filter prior to calibration. With 2 cm/2 degree perturbations for translation/tilting motions, a level achievable under normal HERA operating conditions, the combination of the baseline cut and temporal filtering indicates that the spurious gain feature due to nonredundancies is significantly reduced, and the power spectrum recovers the clean foreground-free region. We found that the mitigation technique works even for large feed motions but in order to keep a stable calibration process, the feed positions need to be constrained to 2 cm for translation motions and 2 degree for tilting offset relative to the dish's vertex.

Strong gravitationally lensed supernovae (glSNe) are a powerful probe to obtain a measure of the expansion rate of the Universe, but they are also extremely rare. To date, only two glSNe with multiple images strongly lensed by galaxies have been found, but their short time delays make them unsuitable for cosmography. We simulate a realistic catalogue of lensed supernovae and study the characteristics of the population of detectable systems for different surveys. Our simulations show that the properties of glSNe in shallow surveys (such as the Zwicky Transient Facility; ZTF) are determined by the need for large magnifications, which favours systems of four images with short time delays and low image separations. This picture is consistent with the properties of iPTF16geu and SN~Zwicky, but is not representative of the population found in deeper surveys, which are limited by the volume of the Universe that is strongly lensed. In our simulations of the Legacy Survey of Space and Time (LSST), glSNe show longer time delays and greater angular separations. Of these systems in LSST, 35\% will allow for time-delay measurements with a precision of 10\% or better. In the 10 years of the survey LSST should be able to find $\approx$ 180 systems, of which 60 will be suited for cosmography enabling a $\approx 1.5 \%$ precision $H_0$ measurement with LSST glSNe.

Primordial black hole (PBH) binaries forming in the early Universe may contribute to the merger events observed by the LIGO-Virgo-KAGRA collaborations. Moreover, the inferred merger rate constraints the fraction of PBH with masses $m \sim 10 \, M_{\odot}$ in the dark matter (DM) to $f_{PBH} \lesssim 10^{-3}$. This constraint assumes that after the formation of PBH binaries, they do not get destroyed or their parameters are not perturbed until the merger. However, PBHs themselves contribute to the formation of early DM structures in which the interactions between PBHs take place actively. This leads to the fact that the binaries can be perturbed in such a way that their lifetime becomes longer than the Hubble time $t_H$. In this work, we consider the effect of the initial spatial Poisson distribution of PBHs on the structure formation at the high redshifts $z \gtrsim 10$. Next, we explore the evolution of such halos due to the interaction of PBHs with each other and with DM particles. We show that the early halos evolve on timescales much shorter than the age of the Universe. Furthermore, for fractions of PBHs $f_{PBH} < 1$, the internal dynamics of a halo is significantly accelerated due to the dynamical friction of PBHs against DM particles. As a result, a significant fraction of binaries will be perturbed in such structures, and the gravitational waves constraints on PBHs with masses $m \sim 10 \, M_{\odot}$ can be weakened to $f_{PBH} \sim 0.1$.

CMB-S4, the next-generation ground-based cosmic microwave background (CMB) observatory, will provide detailed maps of the CMB at millimeter wavelengths to dramatically advance our understanding of the origin and evolution of the universe. CMB-S4 will deploy large and small aperture telescopes with hundreds of thousands of detectors to observe the CMB at arcminute and degree resolutions at millimeter wavelengths. Inflationary science benefits from a deep delensing survey at arcminute resolutions capable of observing a large field of view at millimeter wavelengths. This kind of survey acts as a complement to a degree angular resolution survey. The delensing survey requires a nearly uniform distribution of cameras per frequency band across the focal plane. We present a large-throughput, large-aperture (5-meter diameter) freeform three-mirror anastigmatic telescope and an array of 85 cameras for CMB observations at arcminute resolutions, which meets the needs of the delensing survey of CMB-S4. A detailed prescription of this three-mirror telescope and cameras is provided, with a series of numerical calculations that indicate expected optical performance and mechanical tolerance.

We study the time evolution of the mutual information between the mass distributions in spatially separated but casually connected regions in an expanding universe. The evolution of the mutual information is primarily determined by the configuration entropy rate which depends on the dynamics of the expansion and the growth of the density perturbations. The joint entropy between the distributions from the two regions plays a negligible role in such evolution. The mutual information decreases with time in a matter dominated Universe whereas it stays constant in a $\Lambda$-dominated Universe. The $\Lambda$CDM model and some other models of dark energy predict a minimum in the mutual information beyond which the dark energy dominates the dynamics of the Universe. The mutual information may have deeper connections to the dark energy and the accelerated expansion of the universe.

We present the preliminary timing analysis of confirmed intermediate polar UU Col and possible intermediate polar Swift J0939.7-3224 in the optical band with the help of long-term, high-cadence continuous photometry from Transiting Exoplanet Survey Satellite (TESS). For UU Col, we revise previously reported orbital and spin periods as 3.464 $\pm$ 0.005 h and 863.74 $\pm$ 0.08 s, respectively. Using the second harmonic of the beat frequency, the beat period is estimated as $\sim$928 s. These findings indicate that UU Col is a disc-fed dominated disc-overflow accretor. For J0939, we establish the spin period as 2671.8 $\pm$ 0.8 s and refine the provisionally suggested orbital period as 8.49 $\pm$ 0.03 h. The absence of beat frequency in J0939 signifies that it might be a pure disc-fed accretor; however, an X-ray study of this source will help to understand its true nature.

This paper is devoted to studying the observational signatures modified by Bardeen black hole via shadow and strong lensing observations. Influence of the modified Bardeen black hole parameters q, g, and the parameter $\mu$ on the shadow radius of the black hole have been investigated numerically and graphically. Recently, EHT collaboration observed the image and shadow of supermassive black holes $M87^*$ and $SgrA^*$ where the shadow angular diameter $\theta_d=42\pm3$ for $M87^*$ and $\theta_d=51.8\pm2.3$ for $SgrA^*$. The modified black hole parameters q and $\mu$ for the fixed value of g have been constrained by the EHT collaboration data for the angular shadow diameter of $M87^*$ and $SgrA^*$. It has been observed that the constrain ranges of the parameters $\mu$ and $q$ of modified Bardeen black hole as $-0.89\leq \mu/8M^2 \leq 0.4$ and $0\leq |q|\leq 0.185$ for $M87^*$; and $-1.38\leq \mu/8M^2 \leq 0.1$ and $0\leq |q|\leq 0.058$ for $SgrA^*$, keeping the fixed value $g/2M=0.2$. Modified Bardeen black holes with the additional parameters $\mu$,$g$ and $q$ besides the mass M of the black hole as the supermassive black holes $M87^*$ and $SgrA^*$; and it is observed that to be a viable astrophysical black hole candidate. Furthermore, Gravitational lensing in the strong field limit for modified Bardeen black hole has been investigated numerically as well as graphically and compared to the other ordinary astrophysical black hole such as Schwarzschild ($\mu=\&q=0$) and regular Bardeen ($\mu=0$) black hole.

Recent work has shown that axions can be efficiently produced via non-stationary pair plasma discharges in the polar cap region of pulsars. Here, we point out that for axion masses $10^{-9} \, {\rm eV} \lesssim m_a \lesssim 10^{-4} \, \rm eV$, a sizable fraction of the sourced axion population will be gravitationally confined to the neutron star. These axions accumulate over astrophysical timescales, thereby forming a dense `axion cloud' around the star. We argue that the existence of such a cloud, with densities reaching and potentially exceeding $\mathcal{O}(10^{22}) \, {\rm GeV \, cm^{-3}}$, is a generic expectation across a wide range of parameter space. For axion masses $m_a \gtrsim 10^{-7} \, \rm eV$, energy is primarily radiated from the axion cloud via resonant axion-photon mixing, generating a number of distinctive signatures that include: a sharp line in the radio spectrum of each pulsar (located at the axion mass, and with a percent-level width), and transient events arising from the reconfiguration of charge densities in the magnetosphere. While a deeper understanding of the systematic uncertainties in these systems is required, our current estimates suggest that existing radio telescopes could improve sensitivity to the axion-photon coupling by more than an order of magnitude.

With axions now a primary candidate for dark matter, understanding their indirect astrophysical signatures is of paramount importance. Key to this is the production of photons from axions in magnetised astrophysical plasmas. While simple formulae for axion-photon mixing in 1D have been sketched several decades ago, there has recently been renewed interest in robust calculations for this process in arbitrary 3D plasmas. These calculations are vital for understanding, amongst other things, the radio production from axion dark matter conversion in neutron stars, which may lead to indirect axion dark matter detection with current telescopes or future searches, e.g., by the SKA. In this paper, we derive the relevant transport equations in magnetised plasmas. These equations describe both the production and propagation of photons in an arbitrary 3D medium due to the resonant conversion of axions into photons. They also fully incorporate the refraction of photons, and we find no evidence for a conjectured phenomenon of de-phasing. Our result is free of divergences which plagued previous calculations, and our kinetic theory description provides a direct link between ray tracing and the production mechanism. These results mark an important step toward solving one of the major open questions concerning indirect searches of axions in recent years, namely how to compute the photon production rate from axions in arbitrary 3D plasmas.

Dark matter direct (and indirect) detection experiments usually can only determine a specific combination of a power of the coupling and the dark matter density. This is also true for axion haloscopes which are sensitive to the product $g^{2}_{a\gamma\gamma}\rho_{\rm DM}$, the combination of axion-photon coupling squared and the dark matter density. In this note we show, that in the lucky case when we intersect with a so-called axion minicluster of a suitable size, we can utilize the spectral information available in haloscopes to determine the gravitational potential of the minicluster. We can then use this to measure separately the coupling and the density of the minicluster.

This paper introduces significant improvements to the GravAD pipeline, a Python-based system for gravitational wave detection. These advancements address the complexities of waveform templates and the heightened sensitivity of LIGO detectors. By integrating simulated signals and utilising optimisation techniques, GravAD exhibits increased performance, efficiency, and accuracy in processing GW data. Notable advancements include a reduction in required templates, leading to more efficient detection and freeing computational resources for further research. This pipeline also applies adaptive termination procedures for resource optimisation, enhancing GW detection speed and precision. The paper emphasises the importance of robust, efficient tools in gravitational wave data analysis, particularly given the finite nature of computational resources. Acknowledging system limitations such as dependency on ripple python library capabilities and suggests future enhancements in waveform generation and differentiation.

GW190521 is a short-duration, low-frequency gravitational-wave signal in the LIGO-Virgo catalogue. The signal is consistent with the ringdown and possibly some of the inspiral-merger of an intermediate-mass binary black-hole coalescence. We find that previous models of the quasinormal mode spectrum in the ringdown of GW190521 give remnant mass and spin estimates which are not fully consistent with those of many inspiral-merger-ringdown waveforms. In our own analysis, we find that ringdown models which include both the angular ${l=2}$, ${m=1}$ and ${l=m=2}$ fundamental quasinormal modes are in full agreement with most inspiral-merger-ringdown waveforms, and in particular with the numerical relativity surrogate NRSur7dq4. We also find some support for including the ${l=3}$, ${m=2}$ fundamental quasinormal mode in our fits, building on Capano et al.'s findings regarding a higher-frequency subdominant mode. We propose an interpretation of our GW190521 ringdown model that links precession to the excitation of ${l\neq m}$ quasinormal modes, but we do not rule out eccentricity or other interpretations.

We investigate the cosmological consequences of light sterile neutrinos with altered dispersion relations (ADRs) and couplings to an ultra-light, axion-like scalar field. In particular we study the impact on the number of additional, light, fermionic degrees of freedom and primordial nucleosynthesis. While the ADR leads to a new potential term in the Hamiltonian, the coupling to the scalar field results in a time dependent, effective mass contribution. We solve the quantum kinetic equations (QKEs) for the neutrino density matrix and find that in certain parameter regions both new physics effects can individually yield a suppressed population of sterile neutrino species and the correct observed amount of helium in nucleosynthesis. Combining both effects opens up new patches of parameter space excluded by experimental bounds applying to models featuring only one of the effects.

We present a Bayesian parameter estimation progress to infer the stellar mass binary black hole properties by TianQin, LISA, and TianQin+LISA. Two typical Stellar-mass Black Hole Binary systems, GW150914 and GW190521 are chosen as the fiducial sources. In this work, we establish the ability of TianQin to infer the parameters of those systems and first apply the full frequency response in TianQin's data analysis. We obtain the parameter estimation results and explain the correlation between them. We also find the TianQin+LISA could marginally increase the parameter estimation precision and narrow the $1\sigma$ area compared with TianQin and LISA individual observations. We finally demonstrate the importance of considering the effect of spin when the binaries have a non-zero component spin and great derivation will appear especially on mass, coalescence time and sky location.

In this work, we propose a meta-modelling technique to nuclear matter on the basis of a relativistic density functional with density-dependent couplings. Identical density dependence for the couplings both in the isoscalar and isovector sectors is employed. We vary the coupling parameters of the model to capture the uncertainties of the empirical nuclear matter parameters at saturation. Then, we construct a large ensemble of unified equations of state in a consistent manner both for clusterized and uniform matter in $\beta$-equilibrium at zero temperature. Finally, we calculate neutron star properties to check the consistency with astrophysical observations within a Bayesian framework. Out of the different sets of astrophysical data employed, constraint on tidal deformability from the GW170817 event was found to be the most stringent in the posteriors of different neutron star properties explored in the present study. We demonstrate in detail the impact of the isovector incompressibility ($K_{sym}$) on high-density matter that leads to a considerable variation in the composition of neutron star matter. A couple of selected models with extreme values of $K_{sym}$, which satisfy various modern nuclear physics and neutron star astrophysics constraints, are uploaded in the \textsc{CompOSE} \cite{Typel:2013rza} database for use by the community.

Recently, pulsar timing array (PTA) experiments reported the observation of a stochastic gravitational wave (GW) background in the nanohertz range frequency band. We show that such signal can be originated from a cosmological first-order phase transition (PT) within a well-motivated heavy (visible) QCD axion model. Considering the Peccei-Quinn symmetry breaking at the TeV scale in the scenario, we find a supercooled PT, in the parameter space of the model, prolonging the PT with the reheating temperature at the GeV scale.

Collisionless plasma shocks are a common feature of many space and astrophysical systems and are sources of high-energy particles and non-thermal emission, channeling as much as 20\% of the shock's energy into non-thermal particles. The generation and acceleration of these non-thermal particles have been extensively studied, however, how these particles feed back on the shock hydrodynamics has not been fully treated. This work presents the results of self-consistent hybrid particle-in-cell simulations that show the effect of self-generated non-thermal particle populations on the nature of collisionless, quasi-parallel shocks. They contribute to a significant heat flux density upstream of the shock. Non-thermal particles downstream of the shock leak into the upstream region, taking energy away from the shock. This increases the compression ratio, slows the shock down, and flattens the non-thermal population's spectral index for lower Mach number shocks. We incorporate this into a revised theory for the Rankine-Hugoniot jump conditions that include this effect and it shows excellent agreement with simulations. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies of predictions of arrival times of coronal mass ejections and the conflicting radio and x-ray observations of intracluster shocks. These effects will likely need to be included in fluid modeling to accurately predict shock evolution.

Extreme Mass Ratio Inspirals (EMRIs) are one of the key sources for future space-based gravitational wave interferometers. Measurements of EMRI gravitational waves are expected to determine the characteristics of their sources with sub-percent precision. However, their waveform generation is challenging due to the long duration of the signal and the high harmonic content. Here, we present the first ready-to-use Schwarzschild eccentric EMRI waveform implementation in the frequency domain for use with either graphics processing units (GPUs) or central processing units (CPUs). We present the overall waveform implementation and test the accuracy and performance of the frequency domain waveforms against the time domain implementation. On GPUs, the frequency domain waveform takes in median $0.044$ seconds to generate and is twice as fast to compute as its time domain counterpart when considering massive black hole masses $\geq 2 \times 10^6 \,{\rm M_\odot}$ and initial eccentricities $e_0 > 0.2$. On CPUs, the median waveform evaluation time is $5$ seconds, and it is five times faster in the frequency domain than in the time domain. Using a sparser frequency array can further speed up the waveform generation, reaching up to $ 0.3$ seconds. This enables us to perform, for the first time, EMRI parameter inference with fully relativistic waveforms on CPUs. Future EMRI models which encompass wider source characteristics (particularly black hole spin and generic orbit geometries) will require significantly more harmonics. Frequency-domain models will be essential analysis tools for these astrophysically realistic and important signals.

The Dark Matter Particle Explorer (DAMPE) is a satellite-borne detector designed to measure high energy cosmic-rays and $\gamma$-rays. As a key sub-detector of DAMPE, the Bismuth Germanium Oxide (BGO) imaging calorimeter is utilized to measure the particle energy with a high resolution. The nonlinear fluorescence response of BGO for large ionization energy deposition, known as the quenching effect, results in an under-estimate of the energy measurement for cosmic-ray nuclei. In this paper, various models are employed to characterize the BGO quenching factors obtained from the experimental data of DAMPE. Applying the proper quenching model in the detector simulation process, we investigate the tuned energy responses for various nuclei and compare the results based on two different simulation softwares, i.e. GEANT4 and FLUKA. The BGO quenching effect results in a decrease of the measured energy by approximately $2.5\%$ ($5.7 \%$) for carbon (iron) at $\sim$10 GeV/n and $<1\%$ above 1 TeV/n, respectively. Accordingly, the correction of the BGO quenching effect leads to an increase of the low-energy flux measurement of cosmic-ray nuclei.

We perform an improvement in a van der Waals-type model by including it effects of short-range correlations (SRC). Attractive and repulsive parts of the nucleon-nucleon interaction are assumed to be density-dependent functions, more specifically, we adopt the Carnahan-Starling (CS) method for the latter, and a suitable expression for the former in order to reproduce the structure of the Clausius (C) real gas model. The parametrizations of the resulting model, named as CCS-SRC model, are shown to be capable of reproducing the flow constraint at the high-density regime of symmetric nuclear matter for incompressibility values inside the range of $K_0=(240\pm 20)$ MeV. In the context of stellar matter, our findings point out a good agreement of the CCS-SRC model with recent astrophysical observational data, namely, mass-radius contours and dimensionless tidal deformability regions and values, coming from gravitational waves data related to the GW170817 and GW190425 events, and from the NASA's Neutron star Interior Composition Explorer (NICER) mission. Furthermore, the values for the symmetry energy slope of the model ($L_0$) are in agreement with a recent range found for this quantity, claimed to be consistent with results reported by the updated lead radius experiment (PREX-2) collaboration. In this case, higher values of $L_0$ are favored, while the opposite scenario does not allow simultaneous compatibility between the model and the astrophysical data.

Halo dark matter (DM) particles could lose energy due to the scattering off nuclei within the Earth before reaching the underground detectors of DM direct detection experiments. This Earth shielding effect can result in diurnal modulation of the DM-induced recoil event rates observed underground due to the self-rotation of the Earth. For electron recoil signals from DM-electron scatterings, the current experimental constraints are very stringent such that the diurnal modulation cannot be observed for halo DM. We propose a novel type of diurnal modulation effect: diurnal modulation in electron recoil signals induced by DM-nucleon scattering via the Migdal effect. We set so far the most stringent constraints on DM-nucleon scattering cross section via the Migdal effect for sub-GeV DM using the S2-only data of PandaX-II and PandaX-4T with improved simulations of the Earth shielding effect. Based on the updated constraints, we show that the Migdal effect induced diurnal modulation of electron events can still be significant in the low energy region, and can be probed by experiments such as PandaX-4T in the near future.

We present the design, control system, and noise analysis of a 6-axis seismometer comprising a mass suspended by a single fused silica fibre. We utilise custom-made, compact Michelson interferometers for the readout of the mass motion relative to the table and successfully overcome the sensitivity of existing commercial seismometers by over an order of magnitude in the angular degrees of freedom. We develop the sensor for gravitational-wave observatories, such as LIGO, Virgo, and KAGRA, to help them observe intermediate-mass black holes, increase their duty cycle, and improve localisation of sources. Our control system and its achieved sensitivity makes the sensor suitable for other fundamental physics experiments, such as tests of semiclassical gravity, searches for bosonic dark matter, and studies of the Casimir force.

We investigate the production process of induced gravity waves due to large scalar fluctuations in the paradigm of quintessential inflation. We numerically solve the Mukhanov-Sasaki equation for different sets of parameters to obtain the power spectra. We demonstrate that the induced gravity wave signal generated in this framework can falls within the region of the NANOGrav data for chosen values of model parameters. We show that there is an allowed region of parameter space where the effect shifts to high frequency regime relevant to LISA and other available sensitivities.

In this work, we numerically investigate and visually illustrate the dynamical properties of the dissipative spin-orbit problem such as the co-existence of multiple periodic and quasi-periodic attractors, and the complexity of the corresponding basins of attraction. Our model is composed by a triaxial satellite (planet) orbiting a planet (star) in a fixed Keplerian orbit with zero obliquity. A dissipative tidal torque that is proportional to its rotational angular velocity is assumed to be acting on the satellite. We use Hyperion as a toy model to characterize the methodology used, since this system has a very rich conservative dynamical scenario, and we later apply our methodology to the Moon and Mercury. Our results show that the basins of attraction may possess an intricate structure in all cases which changes with the orbital eccentricity, and that the Gibbs entropy is a good measure on how dominant one basin is over the others in the phase space.