Papers for Tuesday, May 09 2023

Recently, the ANITA collaboration announced the detection of new, unsettling upgoing Ultra-High-Energy (UHE) events. Understanding their origin is pressing to ensure success of the incoming UHE neutrino program. In this work, we study their internal consistency and the implications of the lack of similar events in IceCube. We introduce a generic, simple parametrization to study the compatibility between these two observatories in Standard Model-like and Beyond Standard Model scenarios: an incoming flux of particles that interact with Earth nucleons with cross section $\sigma$, producing particle showers along with long-lived particles that decay with lifetime $\tau$ and generate a shower that explains ANITA observations. We find that the ANITA angular distribution imposes significant constraints, and when including null observations from IceCube only $\tau \sim 10^{-3}$ - $10^{-2} \,\mathrm{s}$ and $\sigma \sim 10^{-33}$ - $10^{-32}\,\mathrm{cm^2}$ can explain the data. This hypothesis is testable with future IceCube data. Finally, we discuss a specific model that can realize this scenario. Our analysis highlights the importance of simultaneous observations by high-energy optical neutrino telescopes and new UHE radio detectors to uncover cosmogenic neutrinos or discover new physics.

The community agrees that Type Ia supernovae arise from Carbon/Oxygen white dwarfs undergoing thermonuclear runaway. However, the full progenitor system and the process that prompts the white dwarf to explode remain unknown. Most current models suggest that the white dwarf explodes because of interaction with a binary companion which may survive the process and remain within the resulting remnant of the exploded star. Furthermore, both the pre-supernova interaction process and the explosion of the primary are expected to imprint a significant departure from ordinary stellar radii and temperatures onto the secondary, making the star identifiable against the unrelated stellar population. Identification of a surviving companion inside an SN Ia remnant might confirm a specific corresponding SN Ia progenitor channel based on the identity of the companion. We conducted a surviving companion search of the Type Ia remnant SNR 0509-67.5 based in the Large Magellanic Cloud. The well-constrained distance to and foreground extinction of the Large Magellanic Cloud allow for Bayesian inference of stellar parameters with low correlation and uncertainties. We present a deep catalog of fully characterized stars interior to SNR 0509-67.5 with radii, effective temperatures, and metallicities inferred using combined Hubble Space Telescope photometric observations across multiple visits. We then compile a list of surviving companion models appropriate for the age of the remnant (roughly 400 years after the explosion). We compare these predictions with the inferred stellar parameters and conclude that none of the stars are consistent with the predicted signatures of a surviving companion.

Optical emission-line ratios are traditionally used to estimate gas metallicities from observed galaxy spectra. While such estimators have been calibrated primarily at low redshift, they are commonly used to study high-redshift galaxies, where their applicability may be questioned. We use comprehensive emission-line catalogues of galaxies from the IllustrisTNG simulation including ionization by stars, active nuclei and shocks to reassess the calibrations of both optical and ultraviolet metallicity estimators at redshifts $0 \geq z \geq 8$. For present-day galaxies, the predicted optical-line calibrations are consistent with previously published ones, while we find different ultraviolet-line ratios, such as HeII$\lambda$1640/CIII]$\lambda$1908, to provide powerful metallicity diagnostics. At fixed metallicity, most emission-line ratios are predicted to strongly increase or decrease with redshift (with the notable exception of N2O2=[NII]$\lambda$6584/[OII]$\lambda$3727), primarily because of a change in ionization parameter. The predicted dependence of R3=[OIII]$\lambda$5007/H$\beta$ and R23=([OII]$\lambda$3727+[OIII]$\lambda\lambda$4959,5007)/H$\beta$, and to a slightly lesser extent R2=[OII]$\lambda$3727/H$\beta$ and O32=[OIII]$\lambda$5007/[OII]$\lambda$3727, on O abundance for galaxies at $z=4-8$ agrees remarkably well with T$_e$-based measurements in 14 galaxies observed with JWST. This success leads us to provide new calibrations of optical and ultraviolet metallicity estimators specifically designed for galaxies at z $>$ 4, to guide interpretations of future, high-redshift spectroscopic surveys. We further demonstrate that applying classical z = 0 calibrations to high-redshift galaxies can bias O-abundance estimates downward by up to 1 dex, leading to the conclusion of a stronger evolution of the mass-metallicity relation than the actual one.

We report on the 2019 XMM-Newton+NuSTAR monitoring campaign of the Seyfert galaxy NGC 2992, observed at one of its highest flux levels in the X-rays. The time-averaged spectra of the two XMM-Newton orbits show Ultra Fast Outflows (UFOs) absorbing structures above 9 keV with $> 3 \sigma$ significance. A detailed investigation of the temporal evolution on a $\sim$5 ks time scale reveals UFO absorption lines at a confidence level $>$95% (2$\sigma$) in 8 out of 50 XMM-Newton segments, estimated via Monte Carlo simulations. We observe a wind variability corresponding to a length scale of 5 Schwarzschild radii $r_S$. Adopting the novel Wind in the Ionised Nuclear Environment (WINE) model, we estimate the outflowing gas velocity ($v=0.21-0.45 c$), column density ($N_H=4-8\cdot 10^{24} cm^{-2}$) and ionisation state ($\log(\xi_0/erg\ cm\ s^{-1})=3.7-4.7$), taking into account geometrical and special relativity corrections. These parameters lead to instantaneous mass outflow rates $\dot{M}_{out}\simeq0.3-0.8 M_{\odot} yr^{-1}$, with associated outflow momentum rates $\dot{p}_{out}\simeq 20-90 L_{Bol}/c$ and kinetic energy rates $\dot{E}_K \simeq 2-25 L_{Bol}$. We estimate a wind duty cycle $\approx$ 12% and a total mechanical power $\approx$ 2 times the AGN bolometric luminosity, suggesting the wind may drive significant feedback effects between the AGN and the host galaxy. Notably, we also provide an estimate for the wind launching radius and density $\approx 5 r_S, 10^{11} {cm}^{-3}$, respectively.

Cool star rotation periods have become an important tool in determining ages of open clusters. We aim to estimate the age of the open cluster NGC 2281 based on the rotational properties of its low-mass members. Previous age estimates for this open cluster range from 275 Myr to 630 Myr. Based on an eight month-long photometric time series obtained at the 1.2 m robotic STELLA telescope in Tenerife, we measured rotation periods for 126 cool star members (70% of the observed members) of NGC 2281. The large set of rotation periods allows us to construct a rich colour-period diagram for NGC 2281 with very few outliers above the slow rotator sequence. We identify an evolved fast rotator sequence which can be used to accurately age date the open cluster relative to other open clusters. Comparisons with M37 and M48 show that all three open clusters are roughly coeval, and we estimate the age of NGC 2281 to be $435\pm50$ Myr. Through comparisons with the younger NGC 3532 and the older Praesepe, we determine the spin down rates of mid-K and early-M fast rotators to be significantly lower than for early-K stars. We suspect that the spin down of early-K fast rotators might be governed by an additional mass dependence. Finally, we show the path towards an empirical description of the evolved fast rotator sequences in open clusters.

Ram pressure stripping of satellite galaxies is thought to be a ubiquitous process in galaxy clusters, and a growing number of observations reveal satellites at different stages of stripping. However, in order to determine the fate of any individual galaxy, we turn to predictions from either simulations or analytic models. It is not well-determined whether simulations and analytic models agree in their predictions, nor the causes of disagreement. Here we investigate the ram pressure stripping of galactic gas in the reference EAGLE hydrodynamical cosmological simulation, and compare the results to predictions from analytic models. We track the evolution of galaxies with stellar mass $M_{*} > 10^{9} \mathrm{M_{\odot}}$ and initial bound gas mass $M_{\mathrm{gas}} > 10^{9} \mathrm{M_{\odot}}$ that fall into galaxy clusters ($M_{\mathrm{200c}} > 10^{14} \mathrm{M_{\odot}}$) between $z = 0.27$ and $z = 0$. We divide each galaxy into its neutral gas disk and hot ionized gas halo and compare the evolution of the stripped gas fraction in the simulation to that predicted by analytic formulations for the two gas phases, as well as to a toy model that computes the motions of gas particles under the combined effects of gravity and a spatially uniform ram pressure. We find that the deviations of the analytic models from the simulation are primarily the result of opposing effects: the enhancement of ram pressure stripping by stellar feedback, and the suppression of stripping due to compaction of the galactic gas.

Pairs of galaxies hosting active galactic nuclei (AGN) are powerful probes of merger-driven supermassive black hole (SMBH) growth as they can resolve individual AGN and trace mergers over a large range of physical separations. To exploit this on a large scale for the first time for both obscured and unobscured AGN, we use photometric redshifts of AGN selected by the Wide-field Infrared Survey Explorer (WISE) to find probabilistic pairs (<100 kpc separations) across the sky, along with a comparison sample of inactive galaxy pairs. Our final sample of integrated pair probabilities yields 198 AGN-AGN pairs (dual AGN) and 2767 AGN-galaxy pairs (offset AGN) with uniformly measured AGN and host galaxy physical properties. We find the fraction of galaxy pairs hosting WISE AGN is dominated by offset AGN and significantly elevated above that of inactive galaxies for large host stellar masses. We show how the AGN merger fraction directly increases with AGN extinction for both offset and dual AGN, with up to ~40% of heavily obscured AGN found in galaxy pairs. Elevated AGN merger fractions coincide with increased host specific star formation rates that suggest merger-driven co-evolution of galaxies and SMBHs. Among dual AGN, the most rapid SMBH growth may occur within the less massive galaxy. Relative to stochastic mechanisms, mergers produce an excess of AGN at increasingly smaller separations, especially for obscured AGN (up to a factor of ~5), and augmented by correlated triggering. Finally, this excess is stronger than for lower luminosity optically-selected AGN, regardless of AGN obscuration level.

Large-scale astrophysics datasets present an opportunity for new machine learning techniques to identify regions of interest that might otherwise be overlooked by traditional searches. To this end, we use Classification Without Labels (CWoLa), a weakly-supervised anomaly detection method, to identify cold stellar streams within the more than one billion Milky Way stars observed by the Gaia satellite. CWoLa operates without the use of labeled streams or knowledge of astrophysical principles. Instead, we train a classifier to distinguish between mixed samples for which the proportions of signal and background samples are unknown. This computationally lightweight strategy is able to detect both simulated streams and the known stream GD-1 in data. Originally designed for high-energy collider physics, this technique may have broad applicability within astrophysics as well as other domains interested in identifying localized anomalies.

From a purely photometric perspective galaxies are generally decomposed into a bulge+disc system, with bulges being dispersion-dominated and discs rotationally-supported. However, recent observations have demonstrated that such framework oversimplifies complexity, especially if one considers galaxy kinematic.To address this issue we introduced with the GPU-based code \textsc{bang} a novel approach that employs analytical potential-density pairs as galactic components, allowing for a computationally fast, still reliable fit of the morphological and kinematical properties of galaxies. Here we apply \textsc{bang} to the SDSS-MaNGA survey, estimating key parameters such as mass, radial extensions, dynamics, for both bulges and discs of +10,000 objects. We test our methodology against a smaller subsample of galaxies independently analysed with an orbit-based algorithm, finding agreement in the recovered total stellar mass. We also manage to reproduce well-established scaling relations, demonstrating how a proper dynamical modelling can result in tighter correlations and provide corrections to standard approaches. Finally, we propose a more general way of decomposing galaxies into "hot" and "cold" components, showing a correlation with orbit-based approaches and visually determined morphological type. Unexpected tails in the "hot-to-total" mass-ratio distribution are present for galaxies of all morphologies, possibly due to visual morphology misclassifications.

Over the past decade, observations of relativistic outflows from outbursting X-ray binaries in the Galactic field have grown significantly. In this work, we present the first detection of moving and decelerating radio-emitting outflows from an X-ray binary in a globular cluster. MAXI J1848-015 is a recently discovered transient X-ray binary in the direction of the globular cluster GLIMPSE-C01. Using observations from the VLA, and a monitoring campaign with the MeerKAT observatory for 500 days, we model the motion of the outflows. This represents some of the most intensive, long-term coverage of relativistically moving X-ray binary outflows to date. We use the proper motions of the outflows from MAXI J1848-015 to constrain the component of the intrinsic jet speed along the line of sight, $\beta_\textrm{int} \cos \theta_\textrm{ejection}$, to be $=0.19\pm0.02$. Assuming it is located in GLIMPSE-C01, at 3.4 kpc, we determine the intrinsic jet speed, $\beta_\textrm{int}=0.79\pm0.07$, and the inclination angle to the line of sight, $\theta_\textrm{ejection}=76^\circ\pm2^{\circ}$. This makes the outflows from MAXI J1848-015 somewhat slower than those seen from many other known X-ray binaries. We also constrain the maximum distance to MAXI J1848-015 to be $4.3$ kpc. Lastly, we discuss the implications of our findings for the nature of the compact object in this system, finding that a black hole primary is a viable (but as-of-yet unconfirmed) explanation for the observed properties of MAXI J1848-015. If future data and/or analysis provide more conclusive evidence that MAXI J1848-015 indeed hosts a black hole, it would be the first black hole X-ray binary in outburst identified in a Galactic globular cluster.

AT 2020mot is a typical UV/optical tidal disruption event (TDE) with no radio or X-ray signatures in a quiescent host. We find an i-band excess and re-brightening along the decline of the light curve which could be due to two consecutive dust echoes from a TDE. We model our observations following van Velzen et al. (2016) and find that the near-infrared light curve can be explained by concentric rings of thin dust within $\sim$0.1 parsecs of a 6e6 M$_{\odot}$ supermassive black hole (SMBH), among the smallest scales at which dust has been inferred near SMBHs. We find dust covering factors of order fc $\leq$ 2%, much lower than found for dusty tori of active galactic nuclei. These results highlight the potential of TDEs for uncovering the environments around black holes when including near-infrared observations in high-cadence transient studies.

Vortex flows have been found in the photosphere, chromosphere and low corona in observations and simulations. It has been suggested that vortices play an important role for channeling energy and plasma into the corona, but the impact of vortex flows on the corona has not directly been studied in a realistic setup. We investigate the role vortices play for coronal heating using high resolution simulations of coronal loops. The vortices are not artificially driven, but arise self-consistently from magnetoconvection. We perform 3D resistive MHD simulations with the MURaM code. Studying an isolated coronal loop in a Cartesian geometry allows us to resolve the structure of the loop interior. We conduct a statistical analysis to determine vortex properties as a function of height from the chromosphere into the corona. We find that the energy injected into the loop is generated by internal coherent motions within strong magnetic elements. A significant part of the resulting Poynting flux is channeled through the chromosphere in vortex tubes forming a magnetic connection between the photosphere and corona. Vortices can form contiguous structures that reach up to coronal heights, but in the corona itself the vortex tubes get deformed and eventually lose their identity with increasing height. Vortices show increased upward directed Poynting flux and heating rate both in the chromosphere and corona, but their effect becomes less pronounced with increasing height. While vortices play an important role for the energy transport and structuring in the chromosphere and low corona, their importance higher up in the atmosphere is less clear since the swirls are less distinguishable from their environment. Vortex tubes reaching the corona show a complex relationship with the coronal emission.

We present high-cadence ultraviolet through near-infrared observations of the Type Ia supernova (SN Ia) 2023bee in NGC~2708 ($D = 32 \pm 3$ Mpc), finding excess flux in the first days after explosion relative to the expected power-law rise from an expanding fireball. This deviation from typical behavior for SNe Ia is particularly obvious in our 10-minute cadence $TESS$ light curve and $Swift$ UV data. Compared to a few other normal SNe Ia with detected early excess flux, the excess flux in SN 2023bee is redder in the UV and less luminous. We present optical spectra of SN 2023bee, including two spectra during the period where the flux excess is dominant. At this time, the spectra are similar to those of other SNe Ia but with weaker Si II, C II and Ca II absorption lines, perhaps because the excess flux creates a stronger continuum. We compare the data to several theoretical models that have been proposed to explain the early flux excess in SNe Ia. Interaction with either a nearby companion star or close-in circumstellar material is expected to produce a faster evolution than seen in the data. Radioactive material in the outer layers of the ejecta, either from a double detonation explosion or simply an explosion with a $^{56}$Ni clump near the surface, can not fully reproduce the evolution either, likely due to the sensitivity of early UV observable to the treatment of the outer part of ejecta in simulation. We conclude that no current model can adequately explain the full set of observations. We find that a relatively large fraction of nearby, bright SNe Ia with high-cadence observations have some amount of excess flux within a few days of explosion. Considering potential asymmetric emission, the physical cause of this excess flux may be ubiquitous in normal SNe Ia.

Planetary debris disks around other stars are analogous to the Asteroid and Kuiper belts in the Solar System. Their structure reveals the configuration of small bodies and provides hints for the presence of planets. The nearby star Fomalhaut hosts one of the most prominent debris disks, resolved by HST, Spitzer, Herschel, and ALMA. Images of this system at mid-infrared wavelengths using JWST/MIRI not only show the narrow Kuiper-Belt-analog outer ring, but also that (1) what was thought from indirect evidence to be an asteroid-analog structure is instead broad, extending outward into the outer system; (2) there is an intermediate belt, probably shepherded by an unseen planet. The newly discovered belt is demarcated by an inner gap, located at ~ 78 au, and it is misaligned relative to the outer belt. The previously known collisionally generated dust cloud, Fomalhaut b, could have originated from this belt, suggesting increased dynamical stirring and collision rates there. We also discovered a large dust cloud within the outer ring, possible evidence of another dust-creating collision. Taken together with previous observations, Fomalhaut appears to be the site of a complex and possibly dynamically active planetary system.

We have measured line widths in active region coronal loops in order to determine whether the non-thermal broadening is anisotropic with respect to the magnetic field direction. These non-thermal velocities are caused by unresolved fluid motions. Our analysis method combines spectroscopic data and a magnetic field extrapolation. We analyzed spectra from the Extreme Ultraviolet Imaging Spectrometer on Hinode. A differential emission measure analysis showed that many spectral lines that are commonly considered to be formed in the active region have a substantial contribution from the background quiet Sun. From these spectra we identified lines whose emission was dominated by the active region loops rather than background sources. Using these lines, we constructed maps of the non-thermal velocity. With data from the Helioseismic Magnetic Imager on the Solar Dynamics Observatory and the Coronal Modeling System nonlinear force-free magnetic field reconstruction code, we traced several of the magnetic field lines through the active region. Comparing the spectroscopic and magnetic data, we looked for correlations of non-thermal velocity with the viewing angle between the line of sight and the magnetic field. We found that non-thermal velocities show a weak anti-correlation with the viewing angle. That is, the tendency is for the non-thermal velocity to be slightly larger in the parallel direction. This parallel broadening may be due to acoustic waves or unresolved parallel flows.

Accreting white dwarfs in Cataclysmic variables (CVs) show short-period (tens of minutes) brightness variations that are consistent with non-radial oscillations similar to gravity (g) modes observed in isolated white dwarfs (WDs). GW Librae, a dwarf nova, was the first CV in which non-radial oscillations were observed and continues to be the best studied accreting WD displaying these pulsations. Unlike isolated WDs, accreting WDs rotate rapidly, with spin periods comparable to or shorter than typical low-order oscillation periods. Accreting WDs also have a different relationship between their interior temperature and surface temperature. The surface temperature of an accreting WD varies on a months to year timescale between dwarf novae accretion events, allowing study of how this temperature change effects g-mode behavior. Here we show results from adiabatic seismological calculations for accreting WDs, focusing on low-order ($\ell=1$) modes. We demonstrate how g-modes vary in response to temperature changes in the subsurface layers due to a dwarf nova accretion event. These calculations include rotation non-perturbatively, required by the high spin rate. We discuss the thermal history of these accreting WDs, and compare the seismological properties with and without rotation. Comparison of $g$-mode frequencies to observed objects may allow inference of features of the structure of the WD such as mass, surface abundance, accretion history, and more. The variation of mode frequencies during cooling after an outburst provides a novel method of identifying modes.

Massive galaxies are expected to grow through different transformative evolutionary phases where high-redshift starburst galaxies and quasars are examples of such phases. The physical mechanisms driving these phases include companion galaxy interactions, active galactic nuclei feedback, and magnetic fields. Our aim is to characterize the physical properties and the environment of the submillimeter galaxy AzTEC-3 at z = 5.3 and the lensed quasar BRI 0952-0115 at z = 4.4, to set a limit on the polarization properties, as well as placing both in the broader context of galaxy evolution. We used full polarization, sub-arcsecond-resolution, ALMA band-7 observations of both BRI 0952-0115 and AzTEC-3 and detect [CII] line emission towards both galaxies, along with companions in each field. We present an updated gravitational lensing model for BRI 0952-0115. We present infrared luminosities, star-formation rates, and [CII] line to infrared luminosity ratios for each source. The [CII] emission line profile for both BRI 0952-0115 and AzTEC-3 exhibit a broad, complex morphology, indicating the possible presence of outflows. We present evidence of a 'gas bridge' between AzTEC-3 and a companion source. Using a simple dynamical mass estimate for the sources, we suggest that both systems are undergoing minor or major mergers. No polarization is detected for the [CII], placing an upper limit below that of theoretical predictions. Our results show that high-velocity wings are detected, indicating possible signs of massive outflows; however, the presence of companion galaxies can affect the final interpretation. Furthermore, the results provide additional evidence in support of the hypothesis that massive galaxies form in overdense regions, growing through interactions. Finally, strong, ordered magnetic fields are unlikely to exist at the kiloparsec scale in the two studied sources.

We present a study on the incidence of major mergers and their impact on the triggering of nuclear activity in 47 type I and 236 type II optically-selected AGN from the MaNGA DR15 sample. From an estimate of non-parametric image predictors ($Gini$, M$_{20}$, concentration (C), asymmetry (A), clumpiness (S), S\'ersic index (n), and shape asymmetry($A_S$)) using the SDSS images, in combination with a Linear Discriminant Analysis Method, we identified major mergers and merger stages. We reinforced our results by looking for bright tidal features in our post-processed SDSS and DESI legacy images. We find a statistically significant higher incidence of major mergers of 29\% $\pm$ 3\% in our type I+II AGN sample compared to 22\% $\pm$ 0.8\% for a non-AGN sample matched in redshift, stellar mass, color and morphological type, finding also a prevalence of post-coalescence (51\% $\pm$ 5\%) over pre-coalescence (23\% $\pm$ 6\%) merger stages. The levels of AGN activity among our massive major mergers are similar to those reported in other works using \oiii\ tracers. However, similar levels are produced by our AGN-galaxies hosting stellar bars, suggesting that major mergers are important promoters of nuclear activity but are not the main nor the only mechanism behind the AGN triggering. The tidal strength parameter $Q$ was considered at various scales looking for environmental differences that could affect our results on the merger incidence, finding non-significant differences. Finally, the H-H$\beta$ diagram could be used as an empirical predictor for the flux coming from an AGN source, useful to correct photometric quantities in large AGN samples emerging from surveys.

We study three astrophysical/cosmological consequences of nonminimal couplings to gravity in wavelike vector dark matter. In the nonrelativistic limit, the nonminimal coupling with the lowest mass dimension leads to effective self-interactions that affect the mass-radius relation of vector solitons, growth of linear perturbations during structure formation, and the speed of gravitational waves (GWs). Based on the success of cold dark matter on large-scale perturbations and the current limits on GW speed, we constrain the dark matter mass and nonminimal coupling strength to be within the range $|\xi_1| / m^2 \ll 10^{50} \mathrm{eV^{-2}}$ and $-3\times 10^{46} \mathrm{eV^{-2}} \lesssim \xi_2 / m^2 \lesssim 8 \times 10^{48} \mathrm{eV^{-2}}$.

Observations of substructure in protoplanetary disks have largely been limited to the brightest and largest disks, excluding the abundant population of compact disks which are likely sites of planet formation. Here, we reanalyze ~0.1'', 1.33 mm ALMA continuum observations of 12 compact protoplanetary disks in the Taurus star-forming region. By fitting visibilities directly, we identify substructures in 6 of the 12 compact disks. We then compare the substructures identified in the full Taurus sample of 24 disks in single star systems and the ALMA DSHARP survey, differentiating between compact (R_eff,90% < 50 au) and extended (R_eff,90% > 50 au) disk sources. We find that substructures are detected at nearly all radii in both small and large disks. Tentatively, we find fewer wide gaps in intermediate-sized disks with R_eff,90% between 30 and 90 au. We perform a series of planet-disk interaction simulations to constrain the sensitivity of our visibility-fitting approach. Under an assumption of planet-disk interaction, we use the gap widths and common disk parameters to calculate potential planet masses within the Taurus sample. We find that the young planet occurrence rate peaks near Neptune masses, similar to the DSHARP sample. For 0.01 $M_J/M_\odot$ $\lesssim$ $M_p/M_*$ $\lesssim$ 0.1 $M_J/M_\odot$, the rate is 17.4$\pm$8.3%; for 0.1 $M_J/M_\odot$ $\lesssim$ $M_p/M_*$ $\lesssim$ 1 $M_J/M_\odot$, it is 27.8$\pm$8.3%. Both of them are consistent with microlensing surveys. For gas giants more massive than 5 $M_J$, the occurrence rate is 4.2$\pm$4.2%, consistent with direct imaging surveys.

Gamma-ray bursts (GRBs) can be classified with their linearly dependent parameters alongside the standard $T_{90}$ distribution. The Generalized linear mixture model(GLM) identifies the number of linear dependencies in a two-parameter space. Classically, GRBs are classified into two classes by the presence of bimodality in the histogram of T$_{90}$. However, additional classes and sub-classes of GRBs are fascinating topics to explore. In this work, we investigate the GRBs classes in the $ T_{90} {-}T_{50}$ plane using the Generalized Linear Models(GLM) for \textit{Fermi} GBM and BATSE catalogs. This study shows five linear features for the \textit{Fermi} GBM catalog and four linear features for the BATSE catalog, directing towards the possibility of more than two GRB classes.

The weak gravitational lensing is a powerful tool in modern cosmology. To accurately measure the weak lensing signal, one has to control the systematic bias to a small level. One of the most difficult problems is how to correct the smearing effect of the Point Spread Function (PSF) on the shape of the galaxies. The chromaticity of PSF for a broad-band observation can lead to new subtle effects. Since the PSF is wavelength dependent and the spectrum energy distributions between stars and galaxies are different, the effective PSF measured from the star images will be different from that smears the galaxies. Such a bias is called colour bias. We estimate it in the optical bands of the Chinese Space Station Survey Telescope from simulated PSFs, and show the dependence on the colour and redshift of the galaxies. Moreover, due to the spatial variation of spectra over the galaxy image, there exists another higher-order bias, colour gradient bias. Our results show that both colour bias and colour gradient bias are generally below $0.1$ percent in CSST. Only for small-size galaxies, one needs to be careful about the colour gradient bias in the weak lensing analysis using CSST data.

Radio-loud active galactic nuclei (RLAGNs) are rare among AGN populations. Lacking high-resolution and high-frequency observations, their structure and evolution stages are not well understood at high redshifts. In this work, we report ALMA 237 GHz continuum observation at $0.023''$ resolution and VLA 44 GHz continuum observation at $0.08''$ resolution of the radio continuum emission from a high-redshift radio and hyper-luminous infrared galaxy at $z=1.92$. The new observations confirm the South-East (SE) and North-West (NW) hotspots identified by previous low-resolution VLA observations at 4.7 and 8.2 GHz and identify a radio core undetected in all previous observations. The SE hotspot has a higher flux density than the NW one does by a factor of 6, suggesting that there can be a Doppler boosting effect in the SE one. In this scenario, we estimate the advance speed of the jet head, ranging from $\sim$0.1c -- 0.3c, which yields a mildly relativistic case. The projected linear distance between the two hotspots is $\sim13$ kpc, yielding a linear size ($\leq20$ kpc) of a Compact-Steep-Spectrum (CSS) source. Combined with new \black{high-frequency ($\nu_\text{obs}\geq44$ GHz) and archived low-frequency observations ($\nu_\text{obs}\leq8.2$ GHz)}, we find that injection spectra of both NW and SE hotspots can be fitted with a continuous injection (CI) model. Based on the CI model, the synchrotron ages of NW and SE hotspots have an order of $10^5$ yr, consistent with the order of magnitude $10^3 - 10^5$ yr observed in CSS sources associated with radio AGNs at an early evolution stage. The CI model also favors the scenario in which the double hotspots have experienced a quiescent phase, suggesting that this RLAGN may have transient or intermittent activities.

Mass accretion is a fundamental process for the growth of supermassive black holes and activating the central engine. However, detailed accretion properties have not been observationally identified at the central ~10 parsec of active galaxies due to its compactness. Here we for the first time ever report the direct detection of parsec-scale (i.e., 0.01% scale of the host galaxy) dense molecular inflow in the active nucleus of the Circinus galaxy. Only a tiny portion (< 3%) of this inflow is consumed in the actual black hole growth but a bulk portion is blown-out by multiphase outflows. The dense gas disk is gravitationally unstable and drives accretion down to the central ~1 parsec, but another process will be required for the final subparsec-scale accretion.

Close-in, sub-stellar companions to $\delta$ Scuti type stars present a highly suitable testbed for examining how planetary-mass objects can influence stellar pulsations. We aim to constrain the mass of HD 31221 b, probe its atmosphere, and demonstrate how it affects the pulsational pattern of its host, HD 31221. We made use of the available data from the short-cadence Transiting Exoplanet Survey Satellite (TESS). We modeled the nine observed transits and the out-of-phase variations, including Doppler beaming, ellipsoidal variations, and the reflection effect. We also incorporated ground-based photometry from the MuSCAT2 imager installed at the 1.52 m Telescopio Carlos Sanchez in the Teide Observatory, Spain. We found HD 31221 b to have an orbital period of $4.66631 \pm 0.00011$ days, with a radius of $1.32 \pm 0.14$ R$_J$ and a mass of $11.5 \pm 10.3$ M$_J$ (from the ellipsoidal effect), making it consistent with either a brown dwarf or a giant planet. As HD 31221 is a rapid rotator ($v \sin I_\star = 175.31 \pm 1.74$ km s$^{-1}$), we deduced the spin-orbit misalignment to be $\lambda = -121.6 \pm 14.4^\circ$ and $I_\star = 55.9 \pm 11.3^\circ$. The phase curve is dominated by the reflection effect, with a geometric albedo of $1.58 \pm 0.50$. We also found evidence that HD 31221 is a $\delta$ Scuti/$\gamma$ Doradus hybrid pulsator. There are three cases for which the $3$rd, $85$th, and $221$st orbital harmonics almost exactly coincide with peaks in the Fourier spectrum of the star, hinting at tidally perturbed stellar oscillations. HD 31221 b is the third substellar object that is found to be disrupting the pulsations of its host, following HAT-P-2 and WASP-33. Additional photometric observations by CHEOPS and/or PLATO can be used to further constrain its mass and provide a more in-depth analysis of its atmosphere.

Astrophysical events that occur in active galactic nucleus (AGN) disks are believed to differ significantly from the ordinary in the interstellar medium. We show that stars located in the outer region of the AGN disk would explode near the original migration starting points instead of being accreted by the central supermassive black hole due to the effect of viscosity. AGN disks provide a dense environment for supernova (SN) explosions, which inevitably involve ejecta-disk interactions. In this paper, we investigate the light curves (LCs) of core-collapse SN exploded in AGN disks. In addition to the fundamental energy source of $^{56} \mathrm{Ni}$--$^{56} \mathrm{Co}$--$^{56} \mathrm{Fe}$ decay reaction powering the SN LCs, the forward-reverse shock produced during interactions may contribute significantly to the observed flux. If the stellar winds manage to create a cavity surrounded by a shell near the star before the SN explosion, the ejecta-winds-disk configurations are expected. We present various SN LCs from different types of progenitors and find that the SN LCs are dominated by the radiation of ejecta-disk interaction-induced shocks. The resulting SNe in the AGN disk is a promising transient source for UV and optical band detection by the Neil Gehrels Swift Observatory (Swift), the Ultraviolet Explorer (UVEX) and wide field survey telescopes such as Ultraviolet Transient Astronomy Satellite (ULTRASAT), Wide Field Survey Telescope (WFST) and Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory. These detections could aid in the investigation of AGN discs and the associated high-energy transient occurrences.

The second generation of stars in the GCs of the MW exhibit unusually high N, Na, or Al, compared to typical Galactic halo stars at similar metallicities. The halo field stars enhanced with such elements are believed to have originated in disrupted GCs or escaped from existing GCs. We identify such stars in the metallicity range -3.0 < [Fe/H] < 0.0 from a sample of ~ 36,800 giant stars observed in the SDSS and LAMOST survey, and present their dynamical properties. The N-rich population and N-normal population among our giant sample do not exhibit similarities in either in their metallicity distribution function or dynamical properties. We find that, even though the MDF of the NRP looks similar to that of the MW's GCs in the range of [Fe/H] < -1.0, our analysis of the dynamical properties does not indicate similarities between them in the same metallicity range, implying that the escaped members from existing GCs may account for a small fraction of our N-rich stars, or the orbits of the present GCs have been altered by the dynamical friction of the MW. We also find a significant increase in the fraction of N-rich stars in the halo field in the very metal-poor (VMP; [Fe/H] < -2.0) regime, comprising up to ~ 20% of the fraction of the N-rich stars below [Fe/H] = -2.5, hinting that partially or fully destroyed VMP GCs may have in some degree contributed to the Galactic halo. A more detailed dynamical analysis of the NRP reveals that our sample of N-rich stars do not share a single common origin. Although a substantial fraction of the N-rich stars seem to originate from the GCs formed in situ, more than 60% of them are not associated with those of typical Galactic populations, but probably have extragalactic origins associated with GSE, Sequoia, and Sagittarius dwarf galaxies, as well as with presently unrecognized progenitors.

The power spectrum and bispectrum of dark matter tracers are key and complementary probes of the Universe. Next-generation surveys will deliver good measurements of the bispectrum, opening the door to improved cosmological constraints and the breaking of parameter degeneracies, from the combination of the power spectrum and bispectrum. Multi-tracer power spectra have been used to suppress cosmic variance and mitigate the effects of nuisance parameters and systematics. We present a bispectrum multi-tracer formalism that can be applied to next-generation survey data. Then we perform a simple Fisher analysis to illustrate qualitatively the improved precision on primordial non-Gaussianity that is expected to come from the bispectrum multi-tracer. In addition, we investigate the parametric dependence of conditional errors from multi-tracer power spectra and multi-tracer bispectra, on the differences between the biases and the number densities of two tracers. Our results suggest that optimal constraints arise from maximising the ratio of number densities, the difference between the linear biases, the difference between the quadratic biases, and the difference between the products $b_1\,b_\Phi$ for each tracer, where $b_\Phi$ is the bias for the primordial potential.

The link between black holes and star formation allows us to draw a connection between black holes and the places and times extraterrestrial intelligences (ETIs) had a greater chance of emerging. Within the context of the gap paradigm for black holes, we show that denser cluster environments that led to gas rich mergers and copious star formation were places less compatible on average with the emergence of ETIs compared to isolated elliptical galaxies by almost two orders of magnitude. The probability for ETIs peaked in these isolated environments around 6 billion years ago and cosmic downsizing shifted the likelihood of ETIs emerging to galaxies with weak black hole feedback, such as in spiral galaxies, at late times.

OH+ absorption is a powerful tracer of inflowing and outflowing gas in the predominantly atomic diffuse and turbulent halo surrounding galaxies. In this letter, we present observations of OH+(1_1-1_0), CO(9-8) and the underlying dust continuum in 5 strongly lensed z~2-4 QSOs, using ALMA to detect outflowing neutral gas. Blue-shifted OH+ absorption is detected in 3/5 QSOs and tentatively detected in a 4th. Absorption at systemic velocities is also detected in one. OH+ emission is observed in 3/5 QSOs at systemic velocities and CO(9-8) is detected in all 5 QSOs at high S/N, providing information on the dense molecular gas within the host galaxy. We compare our sample to high-z far-infrared (FIR) luminous star-forming and active galaxies from the literature. We find no difference in OH+ absorption line properties between active and star-forming galaxies with both samples following the same optical depth-dust temperature relation, suggesting that these observables are driven by the same mechanism in both samples. Similarly, star-forming and active galaxies both follow the same OH+ emission-FIR relation. Obscured QSOs display broader (>800 km/s) emission than the unobscured QSOs and all but one of the high-z star-forming galaxies, likely caused by the warm molecular gas reservoir obscuring the accreting nucleus. Broader CO(9-8) emission (>500 km/s) is found in obscured versus unobscured QSOs, but overall cover a similar range in line widths as the star-forming galaxies and follow the CO(9-8)-FIR luminosity relation found in low-z galaxies. We find that outflows traced by OH+ are only detected in extreme star-forming galaxies (broad CO emission) and in both types of QSOs, which, in turn, display no red-shifted absorption. This suggests that diffuse neutral outflows in galaxy halos may be associated with the most energetic evolutionary phases leading up to and following the obscured QSO phase.

The Granger causality is an econometric test for determining whether one time series is useful for forecasting another one with a certain Granger lag. Here, the light curves in the 0.3-1 keV (reflection dominated, soft) and 1.2-5 keV (continuum dominated, hard) bands of Active Galactic Nuclei (AGNs) are produced, and the Granger lags are estimated and compared to the traditional lag-frequency spectra. We find that the light curves in the hard band Granger-cause (lead) those in the soft band, whereas the obtained Granger lags could be interpreted as the intrinsic reverberation lags. Then, we extract the Granger-lag profiles from 14 XMM-Newton observations of IRAS 13224-3809, and find that the lags are significant in 12 observations. The majority of the obtained Granger (intrinsic) lags are ~ 200-500 s. With the IRAS 13224-3809 mass of 2 $\times 10^6$ $M_{\odot}$, these lags correspond to the true light-travel distance of ~ 20-50 $r_g$. Hence, by assuming a lamp-post geometry and a face-on disc, this places the corona at ~ 10-25 $r_{g}$ above the central black hole. Moreover, multiple Granger lags consisting of the small and large lags of < 500 s and > 1000 s are detected in 4 observations. This means that the corona height can significantly evolve from ~ 10-25 $r_{g}$ to ~ 55 $r_{g}$, or vice versa, towards the end of the observation. Granger lags are a promising way to measure the intrinsic lags, and provide evidence of coronal height variability within each individual observation.

A sample of 895 s-process-rich candidates has been found among the 454180 giant stars surveyed by LAMOST at low spectral resolution (R~1800). In a previous study, taking advantage of the higher resolution (R~86 000) offered by the the HERMES-Mercator spectrograph, we performed the re-analysis of 15 among the brightest stars of this sample. Among these 15 program stars, having close-to-solar metallicities, 11 showed mild to strong heavy element overabundances. The nucleosynthesis process(es) at the origin of these overabundances were however not questioned in our former study. We derive the abundances in s- and r-process elements of the 15 targets in order to investigate whether some stars also show an i-process signature, as sometimes found in their lower metallicity counterparts (namely, the Carbon-Enhanced Metal-Poor (CEMP)-rs stars). Abundances are derived from the high-resolution HERMES spectra for Pr, Nd, Sm, and Eu, using the TURBOSPECTRUM radiative transfer LTE code with MARCS model atmospheres. Using the new classification scheme proposed in our recent study we find that two stars show overabundances in both s- and r-process elements well above the level expected from the Galactic chemical evolution, an analogous situation to the one of CEMP-rs stars at lower metallicities. We compare the abundances of the most enriched stars with the nucleosynthetic predictions from the STAREVOL stellar evolutionary code and find abundances compatible with an i-process occurring in AGB stars. Despite a larger number of heavy elements to characterize the enrichment pattern, the limit between CEMP-s and CEMP-rs stars remains fuzzy. It is however interesting to note that an increasing number of extrinsic stars are found to have abundances better reproduced by an i-process pattern even at close-to-solar metallicities.

We found an active state lasting for ~200 d in the AM CVn star NSV 1440 in 2022. During this state, the object reached a magnitude of 16.5, 2.0-2.5 mag above quiescence, and showed a number of superposed normal outbursts. Such an active state was probably brought either by an enhanced mass-transfer from the secondary or increased quiescent viscosity of the accretion disk. These possibilities are expected to be distinguished by an observation of the interval to the next superoutburst. We also found that the brightness and the course toward the end of the event were similar to the post-superoutburst fading tail in 2021. The mechanism producing the 2022 active state and post-superoutburst fading tails in AM CVn stars may be the same, and the present finding is expected to clarify the nature of these still poorly understood fading tails in AM CVn stars, and potentially of the corresponding phenomenon in hydrogen-rich WZ Sge stars. We also note that the faint, long "superoutbursts" in long-period AM CVn stars claimed in the past were not true outbursts powered by disk instability, but were more likely phenomena similar to the 2022 active state in NSV 1440.

Accretion of mass onto a white dwarf (WD) in a binary system can lead to stellar explosions. If a WD accretes from stellar wind of a distant evolved giant in a symbiotic binary, it can undergo occasional outbursts in which it brightens by several magnitudes, produces a low- and high-velocity mass-outflow, and, in some cases, ejects bipolar jets. In this paper, we complement the current picture of these outbursts by the transient emergence of a neutral region in the orbital plane of symbiotic binaries consisting of wind from the giant. We prove its presence by determining H$^0$ column densities ($N_{\rm H}$) in the direction of the WD and at any orbital phase of the binary by modeling the continuum depression around the Ly$\alpha$ line caused by Rayleigh scattering on atomic hydrogen for all suitable objects, i.e., eclipsing symbiotic binaries, for which a well-defined ultraviolet spectrum from an outburst is available. The $N_{\rm H}$ values follow a common course along the orbit with a minimum and maximum of a few times $10^{22}$ and $10^{24}$ cm$^{-2}$ around the superior and inferior conjunction of the giant, respectively. Its asymmetry implies an asymmetric density distribution of the wind from the giant in the orbital plane with respect to the binary axis. The neutral wind is observable in the orbital plane due to the formation of a dense disk-like structure around the WD during outbursts, which blocks ionizing radiation from the central burning WD in the orbital plane.

The evolutionary sequence for high-mass star formation starts with massive starless clumps that go on to form protostellar, young stellar objects and then compact HII regions. While there are many examples of the three later stages, the very early stages have proved to be elusive. We follow-up a sample of 110 mid-infrared dark clumps selected from the ATLASGAL catalogue with the IRAM telescope in an effort to identify a robust sample of massive starless clumps. We have used the HCO+ (1-0) and HNC (1-0) transitions to identify clumps associated with infall motion and the SiO (2-1) transition to identity outflow candidates. We have found blue asymmetric line profile in 65% of the sample, and have measured the infall velocities and mass infall rates (0.6-$36 \times 10^{-3}$ Msun/yr) for 33 of these clumps. We find a trend for the mass infall rate decreasing with an increase of bolometric luminosity to clump mass i.e. star formation within the clumps evolves. Using the SiO 2-1 line, we have identified good outflow candidates. Combining the infall and outflow tracers reveals that 67% of quiescent clumps are already undergoing gravitational collapse or are associated with star formation; these clumps provide us with our best opportunity to determined the initial conditions and study the earliest stages of massive star formation. Finally, we provide an overview of a systematic high-resolution ALMA study of quiescent clumps selected that allows us to develop a detailed understanding of earliest stages and their subsequent evolution.

We analyzed the 3-mm wavelength spectral line survey of 408 ATLASGAL clumps observed with the IRAM 30m-telescope, focusing on the class I methanol masers with frequencies near 84, 95 and 104.3 GHz. We detect narrow, maser-like features towards 54, 100 and 4 sources in the maser lines near 84, 95 and 104.3 GHz, respectively. Among them, fifty 84 GHz masers, twenty nine 95 GHz masers and four rare 104.3 GHz masers are new discoveries. The new detections increase the number of known 104.3 GHz masers from 5 to 9. The 95 GHz class I methanol maser is generally stronger than the 84 GHz maser counterpart. We find 9 sources showing class I methanol masers but no SiO emission, indicating that class I methanol masers might be the only signpost of protostellar outflow activity in extremely embedded objects at the earliest evolutionary stage. Class I methanol masers that are associated with sources that show SiO line wings are more numerous and stronger than those without such wings. The total integrated intensity of class I methanol masers is well correlated with the integrated intensity and velocity coverage of the SiO (2--1) emission. The properties of class I methanol masers are positively correlated with the bolometric luminosity, clump mass, peak H$_2$ column density of their associated clumps but uncorrelated with the luminosity-to-mass ratio, dust temperature, and mean H$_2$ volume density. We suggest that the properties of class I masers are related to shocks traced by SiO. Based on our observations, we conclude that class I methanol masers at 84 and 95 GHz can trace a similar evolutionary stage as H$_2$O maser, and appear prior to 6.7 and 12.2 GHz methanol and OH masers. Despite their small number, the 104.3 GHz class I masers appear to trace a short and more evolved stage compared to the other class I masers. [abridged]

Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the polarized sky patches of the size of around ten-by-ten arcminutes, at a distance of up to 20 degrees from the Moon, on bright-Moon nights as a wide-field linear polarimetric flat source. We observed 19 patches of the sky adjacent to the bright-Moon with the RoboPol instrument in the SDSS-r broadband filter. These were observed on five nights within two days of the full-Moon across two RoboPol observing seasons. We find that for 18 of the 19 patches, the uniformity in the measured normalized Stokes parameters $q$ and $u$ is within 0.2 %, with 12 patches exhibiting uniformity within 0.07 % or better for both $q$ and $u$ simultaneously, making them reliable and stable wide-field linear polarization flats. We demonstrate that the sky on bright-Moon nights is an excellent wide-field linear polarization flat source. Various combinations of the normalized Stokes parameters $q$ and $u$ can be obtained by choosing suitable locations of the sky patch with respect to the Moon

Galaxy number counts in the $K$-, $H$-, $I$-, $R$-, $B$- and $U$-bands from the Durham Extragalactic Astronomy and Cosmology catalogue could be well-fitted over their whole range using luminosity function (LF) parameters derived from the SDSS at the bright region and required only modest luminosity evolution with the steepening of the LF slope ($\alpha$), except for a sudden steep increase in the $B$-band and a less steep increase in the $U$-band at faint magnitudes that required a starburst evolutionary model to account for the excess faint number counts. A cosmological model treating Hubble expansion as an Einstein curvature required less correction at faint magnitudes than a standard $\Lambda$CDM model, without requiring dark matter or dark energy. Data from DR17 of the SDSS in the $g$, $i$, $r$, $u$ and $z$ bands over two areas of the sky centred on the North Galactic Cap (NGC) and above the South Galactic Cap (SGC), with areas of 5954 and 859 sq. deg., respectively, and a combined count of 622,121 galaxies, were used to construct bright galaxy number counts and galaxy redshift/density plots within the limits of redshift $\leq0.4$ and mag $\leq20$. Their comparative densities confirmed an extensive void in the Southern sky with a deficit of 26\% out to a redshift $z$$\leq$0.15. Although not included in the number count data set because of its incompleteness at fainter magnitudes, extending the SDSS redshift-number count survey to fainter and more distant galaxies with redshift $\leq1.20$ showed a secondary peak in the number counts with many QSOs, bright X-ray and radio sources, and evolving irregular galaxies with rapid star formation rates. This sub-population at redshifts of 0.45--0.65 may account for the excess counts observed in the $B$-band.

We present first results from the KMOS Ultra-deep Rotation Velocity Survey (KURVS), aimed at studying the outer rotation curves shape and dark matter content of 22 star-forming galaxies at $z\sim1.5$. These galaxies represent `typical' star-forming discs at $z \sim 1.5$, being located within the star-forming main sequence and stellar mass-size relation with stellar masses $9.5\leqslant$log$(M_{\star}/\mathrm{M_{\odot}})\leqslant11.5$. We extract individual rotation curves out to 4 times the effective radius, on average, or $\sim 10-15$ kpc. Most rotation curves are flat or rising between three- and six-disc scale radii. Only three objects with dispersion-dominated dynamics ($v_{\rm rot}/\sigma_0\sim0.2$) have declining outer rotation curves at more than 5$\sigma$ significance. After accounting for seeing and pressure support, the nine rotation-dominated discs with $v_{\rm rot}/\sigma_0\geqslant1.5$ have average dark matter fractions of $50 \pm 20\%$ at the effective radius, similar to local discs. Together with previous observations of star-forming galaxies at cosmic noon, our measurements suggest a trend of declining dark matter fraction with increasing stellar mass and stellar mass surface density at the effective radius. Simulated EAGLE galaxies are in quantitative agreement with observations up to log$(M_{\star}R_{\rm eff}^{-2}/\mathrm{M_{\odot}kpc^{-2}}) \sim 9.2$, and over-predict the dark matter fraction of galaxies with higher mass surface densities by a factor of $\sim 3$. We conclude that the dynamics of typical rotationally-supported discs at $z \sim 1.5$ is dominated by dark matter from effective radius scales, in broad agreement with cosmological models. The tension with observations at high stellar mass surface density suggests that the prescriptions for baryonic processes occurring in the most massive galaxies (such as bulge growth and quenching) need to be reassessed.

In the linear theory, the galaxy angular momentum vectors which originate from the initial tidal interactions with surrounding matter distribution intrinsically develop perpendicular alignments with the directions of maximum matter compression, regardless of galaxy mass. In simulations, however, the galaxy spins exhibit parallel alignments in the mass-range lower than a certain threshold, which depends on redshift, web type, and background cosmology. We show that the observed three dimensional spins of the spiral galaxies located on the void surfaces from the Sloan Digital Sky Survey indeed transit from the perpendicular to the parallel alignments with the directions toward the void centers at the threshold zone, $9.51\le\log (M_{\star}/h^{-1}\,M_{\odot})\le10.03$. The null hypothesis of no spin transition is rejected at the 99.9% confidence level by the Kolmogorov-Smirnov test. This study presents a first direct observational evidence for the occurrence of the mass-dependent spin transition of the real galaxies with respect to the non-filamentary structures of the cosmic web, opening a way to constrain the initial conditions of the early universe by measuring the spin transition threshold.

Open clusters (OCs) are infrequent survivors of embedded clusters gestated in molecular clouds. Up to now, little is known about the initial conditions for the formation of OCs. Here, we studied this issue using high-precision astrometric parameters provided by Gaia data release 3. The statistics show that the peculiar motion velocities of OCs vary little from infancy to old age, providing a remarkable opportunity to use OCs to trace their progenitors. Adopting a dynamical method, we derived the masses of the progenitor clumps where OCs were born, which have statistical characteristics comparable to previously known results for clumps observed in the Galaxy. Moreover, the masses of the progenitor clumps of OCs indicate they should be capable of gestating massive O-type stars. In fact, after inspecting the observed OCs and O-type stars, we found that there are many O-type stars in OCs. The destructive stellar feedback from O-type stars may disintegrate the vast majority of embedded clusters, and only those sufficiently dense ones can survive as OCs.

The 4U 1820-30 is a ultra-compact low mass X-ray binary (LMXB) near the center of the globular cluster NGC 6624. Its negative orbital period derivative, observed from the phase evolution of its sinusoidal-like orbital variation, contradicts the positive value obtained from the theoretical prediction. In this paper, we present the analysis of the 4U 1820-30 orbital modulation from light curves obtained from the {\it Neutron star Interior Composition ExploreR (NICER)} observations from 2017 to mid 2022. Combined with historical records, the orbital derivative is measured from the orbital phase evolution between 1976 and 2002 is $\dot P /P =(-5.21 \pm 0.13) \times 10^{-8}$ yr$^{-1}$. No significant second order orbital period derivative is detected with a 2$\sigma$ upper limit of $|\ddot P|<5.48 \times 10^{-22}$ s s$^{-2}$. We discuss the possible intrinsic orbital period derivative of 4U 1820-30 and suggest that this binary system may have a significant mass outflow similar to some other LMXBs. In addition, a periodic modulation with a period of $691.6 \pm 0.7$ s, which is consistent with the superhump period discovered in the far ultraviolet band of the {\it Hubble Space Telescope}, was also detected in in the X-ray light curves collected by {\it NICER}. We conclude that such modulation is probably caused by a period of $0.8 \pm 0.1$ day apsidal precession of accretion disk similar to the SU UMa type dwarf novae and some LMXBs. However we cannot exclude the possibility that it is induced by a hierarchical third star orbiting around the binary system.

This is a technical report for band 4 (550-900 MHz) polarization data with the upgraded GMRT (uGMRT). The report describes the band 4 polarization data analysis procedure and includes notes for observers who are planning polarization observations with the uGMRT. A few pipelines that are currently being used and tested by astronomers at NCRA are discussed as well.

The recent observation of a common red-noise process in pulsar timing arrays (PTAs) suggests that the detection of nanohertz gravitational waves might be around the corner. However, in order to confidently attribute this red process to gravitational waves, one must observe the Hellings-Downs curve -- the telltale angular correlation function associated with a gravitational-wave background. This effort is complicated by the complex modelling of pulsar noise. Without proper care, mis-specified noise models can lead to false-positive detections. Background estimation using bootstrap methods such as sky scrambles and phase shifts, which use the data to characterize the noise, are therefore important tools for assessing significance. We investigate the ability of current PTA experiments to estimate their background with "quasi-independent" scrambles -- characterized by a statistical "match" below the fiducial value: $|M|<0.1$. We show that sky scrambling is affected by "saturation" after ${\cal O}(10)$ quasi-independent realizations; subsequent scrambles are no longer quasi-independent. We show phase scrambling saturates after ${\cal O}(100)$ quasi-independent realizations. With so few independent scrambles, it is difficult to make reliable statements about the $\gtrsim 5 \sigma$ tail of the null distribution of the detection statistic. We discuss various methods by which one may increase the number of independent scrambles. We also consider an alternative approach wherein one re-frames the background estimation problem so that the significance is calculated using statistically \textit{dependent} scrambles. The resulting $p$-value is in principle well-defined but may be susceptible to failure if assumptions about the data are incorrect.

The coupled, time-dependent Gross-Pitaevskii and Ginzburg-Landau equations are solved simultaneously in three dimensions to investigate the equilibrium state and far-from-equilibrium, spin-down dynamics of an interpenetrating neutron superfluid and proton type-II superconductor, as an idealized description of the outer core of a neutron star. The simulations generalize previous calculations without the time-dependent Ginzburg-Landau equation, where proton feedback is absent. If the angle $\theta$ between the rotation and magnetic axes does not equal zero, the equilibrium state consists of geometrically complicated neutron vortex and proton flux-tube tangles, as the topological defects pin to one another locally but align with different axes globally. During spin-down, new types of motion are observed. For $\theta = 0$, entire vortices pair rectilinearly with flux tubes and move together while pinned. For $\theta \neq 0$, vortex segments pair with segments from one or more flux tubes, and the paired segments move together while pinned. The degree to which proton feedback impedes the deceleration of the crust is evaluated as a function of $\theta$ and the pinning strength, $\eta$. Key geometric properties of vortex-flux-tube tangles, such as filament length, mean curvature, and polarity are analysed. It is found that proton feedback smooths the deceleration of the crust, reduces the rotational glitch sizes, and stabilizes the vortex tangle dynamics. The dimensionless control parameters in the simulations are mutually ordered to match what is expected in a real neutron star, but their central values and dynamics ranges differ from reality by many orders of magnitude due to computational limitations.

The Cassiopeia A supernova remnant has a complex structure, manifesting the multidimensional nature of core-collapse supernova explosions. To further understand this, we carried out near-infrared multi-object spectroscopy on the ejecta knots located in the "northeastern (NE) jet" and the "Fe K plume" regions, which are two distinct features in the outer eastern area of the remnant. Our study reveals that the knots exhibit varying ratios of [S II] 1.03 $\mu$m, [P II] 1.189 $\mu$m, and [Fe II] 1.257 $\mu$m lines depending on their locations within the remnant, suggesting regional differences in elemental composition. Notably, the knots in the NE jet are mostly 'S-rich' with weak or no [P II] lines, implying that they originated below the explosive Ne burning layer, consistent with the results of previous studies. We detected no ejecta knots exhibiting only [Fe II] lines in the NE jet area that are expected in the jet-driven SN explosion model. Instead, we discovered a dozen 'Fe-rich' knots in the Fe K plume area. We propose that they are dense knots produced by a complete Si burning with $\alpha$-rich freezeout in the innermost region of the progenitor and ejected with the diffuse X-ray emitting Fe ejecta but decoupled after crossing the reverse shock. In addition to these metal-rich ejecta knots, several knots emitting only He I 1.083 $\mu$m lines were detected, and their origin remains unclear. We also detected three extended H emission features of circumstellar or interstellar origin in this area and discuss its association with the supernova remnant.

The Interstellar Neutral Helium (ISN He) is an important source of information on the physical state of the Local Interstellar Medium. Radiation pressure acting on the neutral helium atoms in the heliosphere has always been neglected, its effect has been considered insignificant compared to gravitational force. The most advanced numerical models of ISN He take into account more and more subtle effects, therefore it is important to check if the effect of radiation pressure is still negligible. In this paper, we use the most up-to-date version of the Warsaw Test Particle Model (WTPM) to calculate the expected helium distribution in the heliosphere, and simulate the flux of ISN He observed by the Interstellar Boundary Explorer (IBEX) and in the future by the Interstellar Mapping and Acceleration Probe (IMAP). We compare results calculated with and without radiation pressure during low and high solar activity. The results show that in the analysis of IBEX-Lo observations the radiation pressure acting on typical helium causes flux differences at a level of 1-4% and is comparable to the observational errors. For the more sensitive IMAP-Lo instrument, there are some regions in the considered observations configurations where radiation pressure causes potentially statistically significant changes in the calculated fluxes. The effect can be up to 9% for the indirect beam and is likely to be higher than the estimated errors. Therefore, we claim that in the future analysis of the IMAP-Lo observations radiation pressure acting on ISN He should be considered.

Modifications in the law of gravity can leave signatures on the large-scale structure, especially on the non-linear and quasi-linear scales. In this work, we propose that the four Minkowski functionals (MFs), which fully encapsulate the morphological information of the large-scale structure, are a powerful tool to constrain modified gravity(MG). With the assistance of N-body simulations, we quantify the MFs' constraining power on the Hu-Sawicki $f(R)$ gravity model by the Fisher matrix technique. The cosmic variance from the finite simulation volume is discussed. By measuring the MFs from the density field of dark matter with different smoothing scales and at different redshifts, we systematically evaluate how much information can be extracted from the MFs in these situations, and try to combine them to improve the constraint. Furthermore, we study the MG signatures imprinted in the MFs of the large-scale structures (LSS) when the biased tracer -- dark matter halo -- is used. And find the response to the fifth force is different from the matter. We also study how the constraints on modified gravity parameters change when using biased tracers of LSS. We expect that the MFs of LSS can provide us with a stringent constraint on modified gravity in the near future.

We construct a catalogue of stellar masses and ages for 696,680 red giant branch (RGB) stars, 180,436 primary red clump (RC) stars, and 120,907 secondary RC stars selected from the LAMOST\,DR8. The RGBs, primary RCs, and secondary RCs are identified with the large frequency spacing ($\Delta \nu$) and period spacing ($\Delta P$), estimated from the LAMOST spectra with spectral SNRs $> 10$ by the neural network method supervised with the seismologic information from LAMOST-Kepler sample stars. The purity and completeness of both RGB and RC samples are better than 95\% and 90\%, respectively. The mass and age of RGBs and RCs are determined again with the neural network method by taking the LAMOST-Kepler giant stars as the training set. The typical uncertainties of stellar mass and age are, respectively, 10\% and 30\% for the RGB stellar sample. For RCs, the typical uncertainties of stellar mass and age are 9\% and 24\%, respectively. The RGB and RC stellar samples cover a large volume of the Milky Way (5 $< R < 20$\,kpc and $|Z| <$\,5\,kpc), which are valuable data sets for various Galactic studies.

The motion of solid particles embedded in gaseous protoplanetary disks is influenced by turbulent fluctuations. Consequently, the dynamics of moderately to weakly coupled solids can be distinctly different from the dynamics of the gas. Additionally, gravitational perturbations from an embedded planet can further impact the dynamics of solids. In this work, we investigate the combined effects of turbulent fluctuations and planetary dust stirring in a protoplanetary disk on three-dimensional dust morphology and on synthetic ALMA continuum observations. We carry out 3D radiative two-fluid (gas+1-mm-dust) hydrodynamic simulations in which we explicitly model the gravitational perturbation of a Jupiter-mass planet. We derived a new momentum-conserving turbulent diffusion model that introduces a turbulent pressure to the pressureless dust fluid to capture the turbulent transport of dust. The model implicitly captures the effects of orbital oscillations and reproduces the theoretically predicted vertical settling-diffusion equilibrium. We find a Jupiter-mass planet to produce distinct and large-scale three-dimensional flow structures in the mm-size dust, which vary strongly in space. We quantify these effects by locally measuring an effective vertical diffusivity (equivalent alpha) and find azimuthally averaged values in a range $\delta_\mathrm{eff}\sim5\cdot 10^{-3} - 2\cdot 10^{-2}$ and local peaks at values of up to $\delta_\mathrm{eff}\sim3\cdot 10^{-1}$. In synthetic ALMA continuum observations of inclined disks, we find effects of turbulent diffusion to be observable, especially at disk edges, and effects of planetary dust stirring in edge-on observations.

Many galaxy clusters show diffuse cluster-scale emission in the form of radio halos, showing that magnetic fields and relativistic electrons are mixed in with the intra-cluster medium (ICM). There is general agreement that the origin of radio halos is connected to turbulence, generated during cluster mergers. Statistical studies of large samples of galaxy clusters in the radio band have the potential to unveil the connection between the properties of radio halos and the mass and dynamics of the host clusters. Previous studies have been limited to massive clusters and based on a small number of radio halos. The aim of this paper is to investigate the scaling relation between the radio power of radio halos and the mass of the host clusters at low frequencies and down to lower cluster masses. We analysed the clusters from the second catalogue of Planck Sunyaev Zel'dovich sources that lie within the 5634 sq deg covered by the second Data Release of the LOFAR Two-meter Sky Survey. We derived the correlation between the radio power and the mass of the host clusters and we investigated the distribution of clusters without radio halos with respect to the correlation. We use X-ray observations to classify the dynamical state of clusters and investigate its role on the power of radio halos. We found a correlation between the power of radio halos at 150 MHz and the mass of the host clusters down to 3e14 Msun. This correlation has a large scatter, part of which can be attributed to the different dynamical states of host clusters. We used two statistical test to show that the distribution of clusters with and without (upper limits) radio halos in the mass-radio power diagram is not compatible with a single correlation and that it is also not compatible with clusters being uniformly distributed below an upper envelope constituted by the correlation.

Evidence of efficient acceleration of cosmic rays in massive young stellar objects has been recently reported. Among these massive protostars, S255 NIRS 3 for which extreme flaring events associated with radio jets have been detected, is one of the best objects to test this hypothesis. We search for gamma-ray emission associated with this object in Fermi-LAT data and inspect the gas content in different molecular lines using the MWISP survey. A GeV source dubbed 4FGL J0613.1+1749c lies on top of the MYSO region, where two filamentary ~10 pc CO structures extend along the same direction of the sub-parsec radio jets. We investigate the spectrum, morphology, and light curve of the gamma-ray source and compare it with the theoretical emission expected from hadronic and leptonic populations accelerated in the radio jets. We argue that the gamma-ray source could be powered by particles accelerated in the S255 NIRS 3 jets, radiating via Bremsstrahlung or proton-proton interaction, and with a synchrotron component shinning in radio from primary or secondary electrons in the case of a leptonic or hadronic population.

Dwarf galaxies have been extensively studied in the Local Group, in nearby groups, and selected clusters, giving us a robust picture of their global stellar and dynamical properties in particular locations in the Universe. Intense study of these properties has revealed correlations between them, including the well known universal stellar mass-metallicity relation. However, since dwarfs play a role in a vast range of different environments, much can be learned about galaxy formation and evolution through extending the study of these objects to various locations. We present MUSE spectroscopy of a sample of 56 dwarf galaxies as a follow-up to the MATLAS survey in low-to-moderate density environments beyond the Local Volume. The dwarfs have stellar masses in the range of $M_{*}/M_{\odot}$ = 10$^{6.1}$-10$^{9.4}$ and show a distance range of D = 14-148 Mpc, the majority (75%) of which are located in the range targeted by the MATLAS survey (10-45 Mpc). We thus report a 75% (79% for dwarf ellipticals) success rate for the semi-automatic identification of dwarf galaxies in the MATLAS survey on the here presented subsample. Using pPXF full spectrum fitting, we determine their line-of-sight velocity and can match the majority of them with their massive host galaxy. Close inspection of their spectra reveals that ~30% show clear emission lines and thus star formation activity. We estimate their stellar population properties (age and metallicity) and compare our results with other works investigating Local Volume and cluster dwarf galaxies. We find that the dwarf galaxies presented in this work show a systematic offset from the stellar mass-metallicity relation towards lower metallicities at the same stellar mass. A similar deviation is present in other works in the stellar mass range probed in this work and might be attributed to the use of different methodologies for deriving the metallicity.

In the context of a project aimed at characterizing the properties of the so-called Bulge Fossil Fragments (the fossil remnants of the bulge formation epoch), here we present the first determination of the metallicity distribution of Liller 1. For a sample of 64 individual member stars we used ESO- MUSE spectra to measure the equivalent width of the CaII triplet and then derive the iron abundance. To test the validity of the adopted calibration in the metal-rich regime, the procedure was first applied to three reference bulge globular clusters (NGC 6569, NGC 6440, and NGC 6528). In all the three cases, we found single-component iron distributions, with abundance values fully in agreement with those reported in the literature. The application of the same methodology to Liller 1 yielded, instead, a clear bimodal iron distribution, with a sub-solar component at $\text{[Fe/H]}= -0.48\,$dex ($\sigma = 0.22$) and a super-solar component at $\text{[Fe/H]}= +0.26\,$dex ($\sigma = 0.17$). The latter is found to be significantly more centrally concentrated than the metal-poor population, as expected in a self-enrichment scenario and in agreement with what found in another bulge system, Terzan 5. The obtained metallicity distribution is astonishingly similar to that predicted by the reconstructed star formation history of Liller 1, which is characterized by three main bursts and a low, but constant, activity of star formation over the entire lifetime. These findings provide further support to the possibility that, similar to Terzan 5, also Liller 1 is a Bulge Fossil Fragment.

The ${\rm H_2NC}$ radical is the high-energy metastable isomer of ${\rm H_2CN}$ radical, which has been recently detected for the first time in the interstellar medium towards a handful of cold galactic sources, besides a warm galaxy in front of the PKS 1830-211 quasar. These detections have shown that the ${\rm H_2CN}$/${\rm H_2NC}$ isomeric ratio, likewise the HCN/HNC ratio, might increase with the kinetic temperature ($T_{\rm kin}$), but the shortage of them in warm sources still prevents us to confirm this hypothesis and shed light about their chemistry. In this work, we present the first detection of ${\rm H_2CN}$ and ${\rm H_2NC}$ towards a warm galactic source, the G+0.693-0.027 molecular cloud (with $T_{\rm kin} > 70 \, {\rm K}$), using IRAM 30m observations. We have detected multiple hyperfine components of the $N_{K_\text{a}K_\text{c}} = 1_{01} - 0_{00}$ and $2_{02} - 1_{01}$ transitions. We derived molecular abundances with respect to ${\rm H_2}$ of (6.8$\pm$1.3)$\times 10^{-11}$ for ${\rm H_2CN}$ and of (3.1$\pm$0.7)$\times 10^{-11}$ for ${\rm H_2NC}$, and a ${\rm H_2CN}$/${\rm H_2NC}$ abundance ratio of 2.2$\pm$0.5. These detections confirm that the ${\rm H_2CN}$/${\rm H_2NC}$ ratio is $\gtrsim$2 for sources with $T_{\rm kin} > 70 \, {\rm K}$, larger than the $\sim$1 ratios previously found in colder cores ($T_{\rm kin}\sim10 \, {\rm K}$). This isomeric ratio dependence with temperature cannot be fully explained with the currently proposed gas-phase formation and destruction pathways. Grain surface reactions, including the ${\rm H_2NC} \rightarrow {\rm H_2CN}$ isomerization, deserve consideration to explain the higher isomeric ratios and ${\rm H_2CN}$ abundances observed in warm sources, where the molecules can be desorbed into the gas phase through thermal and/or shock-induced mechanisms.

A gravitational anomaly is found at weak gravitational acceleration $g_{\rm{N}} < 10^{-9}$ m s$^{-2}$ from analyses of the dynamics of wide binary stars selected from the Gaia EDR3 database that have accurate distances, proper motions, and reliably inferred stellar masses. Implicit high-order multiplicities are required and the multiplicity fraction is calibrated so that binary internal motions agree statistically with Newtonian dynamics at a high enough acceleration of $10^{-8}$ m s$^{-2}$. The observed sky-projected motions and separation are deprojected to the three-dimensional relative velocity $v$ and separation $r$ through a Monte Carlo method, and a statistical relation between the Newtonian acceleration $g_{\rm{N}} \equiv GM/r^2$ (where $M$ is the total mass of the binary system) and a kinematic acceleration $g \equiv v^2/r$ is compared with the corresponding relation predicted by Newtonian dynamics. The empirical acceleration relation at $< 10^{-9}$ m s$^{-2}$ systematically deviates from the Newtonian expectation. A gravitational anomaly parameter $\delta_{\rm{obs-newt}}$ between the observed acceleration at $g_{\rm{N}}$ and the Newtonian prediction is measured to be: $\delta_{\rm{obs-newt}}= 0.034\pm 0.007$ and $0.109\pm 0.013$ at $g_{\rm{N}}\approx10^{-8.91}$ and $10^{-10.15}$ m s$^{-2}$, from the main sample of 26,615 wide binaries within 200 pc. These two deviations in the same direction represent a $10\sigma$ significance. The deviation represents a direct evidence for the breakdown of standard gravity at weak acceleration. At $g_{\rm{N}}=10^{-10.15}$ m s$^{-2}$, the observed to Newton predicted acceleration ratio is $g_{\rm{obs}}/g_{\rm{pred}}=10^{\sqrt{2}\delta_{\rm{obs-newt}}}=1.43\pm 0.06$. This systematic deviation agrees with the boost factor that the AQUAL theory predicts for kinematic accelerations in circular orbits under the Galactic external field.

We present a study of rotation across 30 square degrees of the Orion Star-forming Complex, following a $\sim$200 d photometric monitoring campaign by the Next Generation Transit Survey (NGTS). From 5749 light curves of Orion members, we report periodic signatures for 2268 objects and analyse rotation period distributions as a function of colour for 1789 stars with spectral types F0$-$M5. We select candidate members of Orion using $\textit{Gaia}$ data and assign our targets to kinematic sub-groups. We correct for interstellar extinction on a star-by-star basis and determine stellar and cluster ages using magnetic and non-magnetic stellar evolutionary models. Rotation periods generally lie in the range 1$-$10 d, with only 1.5 per cent of classical T Tauri stars or Class I/II young stellar objects rotating with periods shorter than 1.8 d, compared with 14 per cent of weak-line T Tauri stars or Class III objects. In period$-$colour space, the rotation period distribution moves towards shorter periods among low-mass (>M2) stars of age 3$-$6 Myr, compared with those at 1$-$3 Myr, with no periods longer than 10 d for stars later than M3.5. This could reflect a mass-dependence for the dispersal of circumstellar discs. Finally, we suggest that the turnover (from increasing to decreasing periods) in the period$-$colour distributions may occur at lower mass for the older-aged population: $\sim$K5 spectral type at 1$-$3 Myr shifting to $\sim$M1 at 3$-$6 Myr.

Based on the spectral decomposition through a code of PrepSpec, the light curves (spanning 6.5 years in the observed frame) of the broad-line Balmer decrement, i.e., the flux ratio of the broad \ha to the broad \hb line, are calculated for a sample of 44 Sloan Digital Sky Survey reverberation-mapped quasars ($z<0.53$). It is found that the logarithm of the mean broad-line Balmer decrement is 0.62 with a standard deviation of 0.15 dex. The relations between the mean Balmer decrement and the SMBH accretion properties (the luminosity, black hole mass, Eddington ratio, accretion rate) are investigated and no obvious correlations are found. It is found that there are 27 quasars ($61\%$) showing strong negative correlations between the Balmer decrement variance and the continuum variance, i.e., the Balmer decrement would be smaller with larger continuum flux. Assuming that the dust obscuration leads to the variance in the Balmer decrement and the continuum, an expected slope is $-1/3$, which is not consistent with most of measured slopes. Using the interpolated cross-correlation function, the time delays between the inverse Balmer decrement and the continuum are measured for 14 quasars with the maximum correlation coefficient larger the 0.6. It suggests that the size corresponding to the Balmer decrement lag extends from the BLR size to the torus size.

The substantial populations of massive quiescent galaxies at $z\ge3$ challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at $z> 3$. As a result, studies report that as much as 70\% of quiescent galaxies at $z> 3$ may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at $3\lesssim z \lesssim 6$ using JWST NIRCam imaging data. We use empirically-constrained galaxy SED templates to define a region in the $F277W-F444W$ vs. $F150W-F277W$ color plane that appears unique in capturing quiescent galaxies at $z> 3$ and minimizes contamination from other red galaxy populations. We apply this color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and filter out $> 99\%$ of sources. We identify 44 candidate $z\gtrsim3$ quiescent galaxies and derive volume density estimates of $n\sim1-4\times10^{-5}$ Mpc$^{-3}$ at $3< z< 5$, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam's wavelength coverage and sensitivity, this technique provides an efficient filter that aids in the search for these rare galaxies.

The dynamics of repeating fast radio bursts (FRBs) are driven by their physical nature and central engine, however their event rate, energy distribution and temporal occurrence behaviour are still remain uncertain due to the server lack of information of bursts. Recently, the available of high-frequency observation data for the Five-hundred-meter Aperture Spherical radio Telescope (FAST) has made it possible to statistically study the temporal occurrence on timescales from several milliseconds to over several thousand seconds. In this research we studied both the FRB121102 and FRB20201124A temporal occurrence and report here a statistical result about the behaviour of the waiting time (or recurrence-time) between successive bursts. The results exhibit novel scaling and universality which have not reported in the field yet. Specifically, we find the scaling law for FRBs recurrence-time distribution which is a clear indication of the importance of correlations in the structure of its physical nature and central engine. The scaling relationships were observed for time scales spanning three orders of magnitude. Given that they are sharing the same scaling law between two repeating FRBs, we infer that the scaling law of waiting time distribution should acts as a indicator which provides insights into the physical nature and the development of the central engine model.

The proto-clusters in the epoch of reionization, traced by galaxies overdensity regions, are ideal laboratories to study the process of stellar assembly and cosmic reionization. We present the spectroscopic confirmation of the core of the most distant proto-cluster at $z = 7.88$, A2744-z7p9OD, with the James Webb Space Telescope NIRSpec integral field unit spectroscopy. The core region includes as many as 4 galaxies detected in [O III] 4960 A and 5008 A in a small area of $\sim 3" \times 3"$, corresponding to $\sim$ 11 kpc $\times$ 11 kpc. Three member galaxies are also tentatively detected in dust continuum in ALMA Band 6, which is consistent with their red ultraviolet continuum slopes, $\beta \sim -1.3$. The member galaxies have stellar masses in the range of log($M_{*}/M_{\rm \odot}$) $\sim 7.6-9.2$ and star formation rates of $\sim 3-50$ $M_{\rm \odot}$ yr$^{-1}$, showing a diversity in their properties. FirstLight cosmological simulations reproduce the physical properties of the member galaxies including the stellar mass, [OIII] luminosity, and dust-to-stellar mass ratio, and predict that the member galaxies are on the verge of merging in a few to several tens Myr to become a large galaxy with $M_{\rm *}\sim 6\times10^{9} M_{\rm \odot}$. The presence of a multiple merger and evolved galaxies in the core region of A2744-z7p9OD indicates that environmental effects are already at work 650 Myr after the Big Bang.

Short-lived radionuclides, such as 26Al and 60Fe, are tracers of star formation. Therefore, their abundances can unravel the recent star formation history of the host galaxy. In view of future gamma-ray surveys, we predict the masses and fluxes of these two elements in the Large Magellanic Cloud (LMC) using new chemical evolution models. Our best model reproduces the abundance patterns of alpha/Fe ratios, the gas mass, the average metallicity, the present time supernova and nova rates observed in LMC. We show three main results: i) the best model for the LMC suggests a star formation rate very mild at the beginning with a recent burst, and a Salpeter-like initial mass function. ii) The predicted mass of 26Al is 0.33 M_Sun, 2/3 produced by massive stars and 1/3 by novae. iii) The predicted mass of 60Fe is 0.44 M_Sun, entirely produced by massive stars. This result suggests a larger fraction of 60Fe, at variance with the Milky Way. The explanation for this lies in the adopted initial mass function, that for the LMC contains more massive stars than for the Milky Way. These predictions can be useful for the COSI-SMEX mission planned for launch in 2027. The expected gamma-ray line fluxes for the 1.809 MeV line of 26Al and the 1.173 and 1.332 MeV lines of 60Fe are in the range of (0.2-2.7)x10^-6 ph cm^-2 s^-1 and (0.7-2.8)x10^-7 ph cm^-2 s^-1, respectively. This new instrument could have the sensitivity to detect the upper end of the predicted 1.8 MeV flux within its nominal two-year mission.

Several short-chain fatty acids and their corresponding potential existing hydrated forms are important molecules in interstellar space. Their structures were optimized with twelve different computational methods. The dipole moments and the spectral constants, including rotational constants and centrifugal distortion constants were obtained. According to the benchmark study, revDSD-PBEP86-D3(BJ) is the most suitable method that was selected for rotational calculation. Symmetry-adapted perturbation theory was used to study the strength and composition of the interaction between acids and water in clusters. The possibility of its existing under the low-temperature and low-pressure conditions was confirmed by calculating of binding free energy. Furthermore, ab initio molecular dynamics simulations were used to investigate whether the internal rotations of acids could be observed. The 3-fold splitting from the predicted high-resolution microwave rotational spectra of the acetic acid monohydrate at different temperatures perfectly proved the accuracy of the simulations.

We applied a scientific approach to the problem of the effective area cross-calibration of the XMM-Newton EPIC instruments. Using a sample of galaxy clusters observed with the XMM-Newton EPIC, we quantified the effective area cross-calibration bias between the EPIC instruments as implemented in the public calibration data base on November 2021 in the 0.5-6.1 keV energy band. We invested significant efforts in controlling and minimising the systematic uncertainties of the cross-calibration bias below 1%. The statistical uncertainties are similar and thus we can reliably measure effects at 1% level. The effective area cross-calibration in the 0.5-6.1 keV band between MOS and pn is biased at a substantial level. The MOS/pn bias is systematic suggesting that MOS (pn) effective area may be calibrated too low (high), by $\sim$3-27% on average depending on the instrument and energy band. The excellent agreement of the energy dependencies (i.e. shapes) of the effective area of MOS2 and pn suggest that they are correctly calibrated within $\sim$1% in the 0.5-4.5 keV band. Comparison with an independent data set of point sources (3XMM) confirms this. The cluster sample indicates that the MOS1/pn effective area shape cross-calibration has an approximately linear bias amounting to $\sim$10% in maximum in the 0.5-4.5 keV band. The effective area cross-calibration of XMM-Newton/EPIC instruments on November 2021 in the 0.5-4.5 keV band is in a relatively good shape. However, the cluster-to-cluster rms scatter of the bias is substantial compared to the median bias itself. Thus, a statistically robust implementation of the cross-calibration uncertainties to a scientific analysis of XMM-Newton/EPIC data should include the propagation of the scatter to the best-fit parameters, instead of a simple average bias correction of the effective area.

We present magnetohydrodynamic simulations of a jet-wind interaction in a galaxy cluster and the radio to gamma-ray and the neutrino emissions from this "head-tail galaxy". Our simulation follows the evolution of cosmic-ray (CR) particle spectra with energy losses and the stochastic turbulence acceleration. We find that the reacceleration is essential to explain the observed radio properties of head-tail galaxies, in which the radio flux and spectral index do not drastically change. Our models suggest that hard X-ray emissions can be detected around the head-tail galaxy in the Perseus cluster by the hard X-ray satellites, such as FORCE, and it will potentially constrain the acceleration efficiency. We also explore the origin of the collimated synchrotron threads, which are found in some head-tail galaxies by recent high-quality radio observations. Thin and elongated flux tubes, connecting the two tails, are formed by strong backflows at an early phase. We find that these threads advect with the wind for over 300 Myr without disrupting. The radio flux from the flux tubes is much lower than the typical observed flux. An efficient CR diffusion process along the flux tubes, however, may solve this discrepancy.

Gas-phase molecules in cometary atmospheres (comae) originate primarily from (1) outgassing by the nucleus, (2) sublimation of icy grains in the near-nucleus coma, and (3) coma (photo-)chemical processes. However, the majority of cometary gases observed at radio wavelengths have yet to be mapped, so their production/release mechanisms remain uncertain. Here we present observations of six molecular species towards comet 46P/Wirtanen, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) during the comet's unusually close (~0.1 au) approach to Earth in December 2018. Interferometric maps of HCN, CH$_3$OH, CH$_3$CN, H$_2$CO, CS and HNC were obtained at an unprecedented sky-projected spatial resolution of up to 25 km, enabling the nucleus and coma sources of these molecules to be accurately quantified. The HCN, CH$_3$OH and CH$_3$CN spatial distributions are consistent with the production from direct outgassing from (or very near to) the nucleus, with a significant proportion (~50%) of the observed CH$_3$OH originating from sublimation of icy grains in the near-nucleus coma. On the other hand, H$_2$CO, CS and HNC originate primarily from distributed coma sources. The HCN, CH$_3$OH and HNC abundances in 46P are consistent with the average values previously observed in comets, whereas the H$_2$CO, CH$_3$CN and CS abundances are relatively low.

Efficient exploration of parameter spaces is crucial to extract physical information about the Epoch of Reionization from various observational probes. To this end, we propose a fast technique based on Gaussian Process Regression (GPR) training applied to a semi-numerical photon-conserving reionization model, SCRIPT. Our approach takes advantage of the numerical convergence properties of SCRIPT and constructs a training set based on low-cost, coarse-resolution simulations. A likelihood emulator is then trained using this set to produce results in approximately two orders of magnitude less computational time than a full MCMC run, while still generating reasonable 68% and 95% confidence contours. Furthermore, we conduct a forecasting study using simulated data to demonstrate the applicability of this technique. This method is particularly useful when full MCMC analysis is not feasible due to expensive likelihood computations.

Neutron stars might have multipole magnetic fields as implied by recent observations of pulsars. The presence of the quadrupole field might have an effect on the interaction between the disc and the neutron star depending on the location of the inner radius of the disc and the strength of the quadrupole field. For a quadrudipole stellar field, we calculate the toroidal field generated within the disc, the magnetospheric radius and the torque exerted onto the star. Also, we deduce the effect of the rotation of the star on the magnetospheric radius which is relevant even for pure dipole magnetic fields.

Polarization phase-curve measurements provide a unique constraint on the surface properties of asteroids that are complementary to those from photometry and spectroscopy, and have led to the identification of the ``Barbarian'' asteroids as a class of objects with highly unusual surfaces. We present new near-infrared polarimetric observations of six Barbarian asteroids obtained with the WIRC+Pol instrument on the Palomar Hale telescope. We find a dramatic change in polarimetric behavior from visible to near-infrared for these objects, including a change in the polarimetric inversion angle that is tied to the index of refraction of the surface material. Our observations support a two-phase surface composition consisting of high albedo, high index of refraction inclusions with a small optical size scale embedded in a dark matrix material more closely related to C-complex asteroids. These results are consistent with the interpretation that the Barbarians are remnants of a population of primitive bodies that formed shortly after CAIs. Near-infrared polarimetry provides a direct test of the constituent grains of asteroid surfaces.

The catalog of km-sized near-Earth objects (NEOs) is nearly complete. Typical impact monitoring analyses search for possible impacts over the next 100 years and none of the km-sized objects represent an impact threat over that time interval. Assessing the impact risk over longer time scales is a challenge since orbital uncertainties grow. To overcome this limitation we analyze the evolution of the Minimum Orbit Intersection Distance (MOID), which bounds the closest possible encounters between the asteroid and the Earth. The evolution of the MOID highlights NEOs that are in the vicinity of the Earth for longer periods of time, and we propose a method to estimate the probability of a deep Earth encounter during these periods. This metric is used to rank the km-sized catalog in terms of their long-term impact hazard to identify targets of potential interest for additional observation and exploration.

The paper presents the first results of the ongoing spectropolarimetric monitoring of magnetic fields of stars, whose chemically peculiar nature has been previously revealed with the 1-m SAO RAS telescope. We selected the sample candidates using the photometric data of the Kepler and TESS space missions. The efficiency of the method of searching for new CP stars based on photometric light curves has been confirmed. We present the magnetic field measurements and estimate the atmospheric parameters of the objects under study.

The Red Nova ZTF SLRN-2020 is the third transient event with properties that are compatible with the merger of a planet with a main sequence (or close to) star on a dynamical timescale. While the two first transient events occurred in young systems, ZTF SLRN-2020 occurred in an old system. Nonetheless, I show that the three star-planet intermediate luminosity optical transients (ILOTs, also termed Red Novae) occupy the same area in the energy-time diagram of ILOTs. Based on models for ILOTs that are power by stellar binary interaction I suggest that the planet in ZTF SLRN-2020 launched jets at about its escape speed before it was engulfed by the star. Interestingly, the escape speed from the planet is similar to the orbital speed of the planet. This leads to an outflow with a very low terminal velocity, much below the escape velocity from the star, and in concentration around ~45 degrees to the equatorial plane. As well, the planet might have lost back some of the accreted mass just before engulfment, forming an accretion disk around the star. This disk might have launched jets during the main outburst of the event. The jets form a bipolar expanding nebula.

I show that exoplanets can be used to test origins scenarios. Origins scenarios start with certain initial conditions, proceed via a network of chemical reactions and, if successful, result in a chemistry that is closer to a living system than the initial conditions. Exoplanet environments can be applied to test each of these three aspects of origins scenarios. I show what tests can be applied to the UV-driven cyanosulfidic scenario and how the application of some of these tests has already falsified certain versions of this scenario. Testing initial conditions has replaced certain reactants with others and has affected the overall chemical network underlying the cyanosulfidic scenario. The sequence of reactions the scenario invokes provide a predicted upper limit on the ubiquity of life in the universe that has ample room for improvement. The outcome of the experiments in different environments is part of a predicted distribution of biosignature detections that can be compared to future observed distributions.

Small planets located at the lower mode of the bimodal radius distribution are generally assumed to be composed of iron and silicates in a proportion similar to that of the Earth. However, recent discoveries are revealing a new group of low-density planets that are inconsistent with that description. We intend to confirm and characterize the TESS planet candidate TOI-244.01, which orbits the bright ($K$ = 7.97 mag), nearby ($d$ = 22 pc), and early-type (M2.5 V) M-dwarf star GJ 1018 with an orbital period of 7.4 days. We used Markov Chain Monte Carlo methods to model 57 precise radial velocity measurements acquired by the ESPRESSO spectrograph together with TESS photometry and complementary HARPS data. We find TOI-244 b to be a super-Earth with a radius of $R_{\rm p}$ = 1.52 $\pm$ 0.12 $\rm R_{\oplus}$ and a mass of $M_{\rm p}$ = 2.68 $\pm$ 0.30 $\rm M_{\oplus}$. These values correspond to a density of $\rho$ = 4.2 $\pm$ 1.1 $\rm g \cdot cm^{-3}$, which is below what would be expected for an Earth-like composition. We find that atmospheric loss processes may have been efficient to remove a potential primordial hydrogen envelope, but high mean molecular weight volatiles such as water could have been retained. Our internal structure modeling suggests that TOI-244 b has a $479^{+128}_{-96}$ km thick hydrosphere over a 1.17 $\pm$ 0.09 $\rm R_{\oplus}$ solid structure composed of a Fe-rich core and a silicate-dominated mantle compatible with that of the Earth. On a population level, we find two tentative trends in the density-metallicity and density-insolation parameter space for the low-density super-Earths, which may hint at their composition. With a 8$\%$ precision in radius and 12$\%$ precision in mass, TOI-244 b is among the most precisely characterized super-Earths, which, together with the likely presence of an extended hydrosphere, makes it a key target for atmospheric observations.

We study the formation of primordial black holes (PBHs) in strongly supercooled first-order phase transitions. The mechanism is based on the presence of remnants dominated by the false vacuum that scale slower with the expansion of the Universe than their surroundings where this energy was already converted into radiation. We compute the PBH formation from these remnants including the contribution from the false vacuum and the bubble walls, by estimating the collapse using the hoop conjecture and by considering both regions collapsing immediately when entering the horizon and sub-horizon regions that collapse as their compactness grows. We show that for exponential bubble nucleation rate, $\Gamma \propto e^{\beta t}$, the primordial black hole formation implies $\beta/H \gtrsim 3.8$, where $H$ denotes the Hubble rate, if the potential energy of the false vacuum is $\Delta V \lesssim (10^{12} {\rm GeV})^4$, as otherwise a too large abundance of long-lived PBHs forms. The observed dark matter abundance can be formed in asteroid mass PBHs if $\beta/H \simeq 3.8$ and $10^5 {\rm GeV} \lesssim \Delta V^{1/4} \lesssim 10^8 {\rm GeV}$. Finally, we consider also the effect of the second order correction to the exponential nucleation rate showing that the PBH abundance is mainly determined by the average radius of the true vacuum bubbles.

Recent studies reveal that more than a dozen of white dwarfs displaying near-perfect blackbody spectra in the optical range have been lurking in the Sloan Digital Sky Survey catalog. We point out that, in a way analogous to the Cosmic Microwave Background, these stars serve as excellent testbeds for new physics. Specifically, we show how their observed lack of spectral distortions translates into limits on the parameter space of axions with electromagnetic coupling. The prospects for future improvements are also discussed.

Axion-like particles (ALPs) are leading candidates to explain the dark matter in the universe. Their production via the misalignment mechanism has been extensively studied for cosine potentials characteristic of pseudo-Nambu-Goldstone bosons. In this work we investigate ALPs with non-periodic potentials, which allow for large misalignment of the field from the minimum. As a result, the ALP can match the relic density of dark matter in a large part of the parameter space. Such potentials give rise to self-interactions which can trigger an exponential growth of fluctuations in the ALP field via parametric resonance, leading to the fragmentation of the field. We study these effects with both Floquet analysis and lattice simulations. Using the Press-Schechter formalism, we predict the halo mass function and halo spectrum arising from ALP dark matter. These halos can be dense enough to produce observable gravitational effects such as astrometric lensing, diffraction of gravitational wave signals from black hole mergers, photometric microlensing of highly magnified stars, perturbations of stars in the galactic disk or stellar streams. These effects would provide a probe of dark matter even if it does not couple to the Standard Model. They would not be observable for halos predicted for standard cold dark matter and for ALP dark matter in the standard misalignment mechanism. We determine the relevant regions of parameter space in the (ALP mass, decay constant)-plane and compare predictions in different axion fragmentation models.

I compute the rate of change of mass and angular momentum of a black hole, namely tidal heating, in an eccentric orbit. The change is caused due to the tidal field of the orbiting companion. I compute the result for both the spinning and non-spinning black holes in the leading order of the mean motion, namely $\xi$. I demonstrate that the rates get enhanced significantly for nonzero eccentricity. Since eccentricity in a binary evolves with time I also express the results in terms of an initial eccentricity and azimuthal frequency $\xi_{\phi}$. In the process, I developed a prescription that can be used to compute all physical quantities in a series expansion of initial eccentricity, $e_0$. This result was only known in the leading order while ignoring the contribution of the spin on the eccentricity evolution. Although the eccentricity evolution result still ignores the spin effect in the current work, the prescription can be used to compute higher-order corrections of initial eccentricity post-leading order. Using this result I computed the rate of change of mass and angular momentum of a black hole in terms of initial eccentricity and azimuthal frequency up to $\mathcal{O}(e_0^2)$.

We describe cosmological solutions of the holographic dilaton with the aim of exploring alternatives to the commonly studied thermal Randall-Sundrum phase transition. It is well known that the thermal transition is typically strongly first order, with the requirement of a perturbative 5D gravity theory obstructing completion of the transition. This thermal transition corresponds to nucleation of an infrared brane through the surface of an AdS-Schwarzschild horizon. The approach we study instead invokes an early epoch in which the cosmology is fully 5-dimensional, with highly relativistic brane motion, and with Rindler horizons obscuring the infrared brane at early times. Our approach corresponds, via AdS/CFT, to a non-equilibrium approach to the conformal phase transition. We comment on a class of initial conditions that generically leads to completion of the phase transition without sacrificing perturbativity of the 5D theory.

We study a simple extension of the Standard Model featuring a dark sector with an ultralight pseudo Nambu-Goldstone boson as dark matter candidate. We focus on the mass range $\mathcal{O}(10^{-20} - 10^{-10})$ eV, relevant for strong gravity applications, and explore its production and evolution in the early Universe. The model is formulated in such a way that dark matter does not couple directly to photons or other Standard Model particles avoiding some of the most stringent cosmological bounds related to axion-like particles. In this work, two different scenarios are considered depending on whether dark matter is produced in a pre-inflationary or post-inflationary regime. We also discuss the effect from emergent topological defects such as cosmic strings and domain walls, and estimate the spectrum of stochastic gravitational waves produced by their decay, enabling to test the model at current and future gravitational-wave experiments.

We provide a pedagogical review article on fundamentals and applications of the quantum dynamics in strong electromagnetic fields in QED and QCD. The fundamentals include the basic picture of the Landau quantization and the resummation techniques applied to the class of higher-order diagrams that are enhanced by large magnitudes of the external fields. We then discuss observable effects of the vacuum fluctuations in the presence of the strong fields, which consist of the interdisciplinary research field of nonlinear QED. We also discuss extensions of the Heisenberg-Euler effective theory to finite temperature/density and to non-Abelian theories with some applications. Next, we proceed to the paradigm of the dimensional reduction emerging in the low-energy dynamics in the strong magnetic fields. The mechanisms of superconductivity, the magnetic catalysis of the chiral symmetry breaking, and the Kondo effect are addressed from a unified point of view in terms of the renormalization-group method. We provide an up-to-date summary of the lattice QCD simulations in magnetic fields for the chiral symmetry breaking and the related topics as of the end of 2022. Finally, we discuss novel transport phenomena induced by chiral anomaly and the axial-charge dynamics. Those discussions are supported by a number of appendices.

We investigate gravitational lensing in the strong deflection regime by loop quantum gravity (LQG)-motivated rotating black hole (LMRBH) metrics with an additional parameter $l$ besides mass $M$ and rotation $a$. The LMRBH spacetimes are regular everywhere, asymptotically encompassing the Kerr black hole as a particular case and, depending on the parameters, describe black holes with one horizon only (BH-I), black holes with an event horizon and a Cauchy horizon (BH-II), black holes with three horizons (BH-III), or black holes with no horizons (NH) spacetime. It turns out that as the LQG parameter $l$ increases, the unstable photon orbit radius $x_{ps}$, the critical impact parameter $u_{ps}$, the deflection angle $\alpha_D(\theta)$ and angular position $\theta_{\infty}$ also increases. Meanwhile, the angular separation $s$ decreases, and relative magnitude $r_{mag}$ increases with increasing $l$ for prograde motion but they show opposite behaviour for the retrograde motion. For Sgr A*, the angular position $\theta_{\infty}$ is $\in$ (16.404, 39.8044) $\mu$as, while for M87* $\in$ (12.3246, 29.9057) $\mu$as. The angular separation $s$ is ranging $\in$ (0.008306-0.37573) $\mu$as for Sgr A* and $\in$ (0.00624-0.282295) $\mu$as for M87*. The relative magnitude $r_{mag}$ $\in$ (0.04724, 1.53831). We estimate the time delay between the first and second relativistic images using twenty supermassive galactic centre black holes as lenses. Our analysis concludes that, within the $1 \sigma$ region, a significant portion of the BH-I and BH-II and for a small portion of BH-III parameter space agrees with the EHT results of M87* and Sgr A* whereas NH is completely ruled out. We discover that the EHT results of Sgr A* place more stringent limits on the parameter space of LMRBH black holes than those established by the EHT results of M87*.

A supermassive binary black-hole candidate SDSS J1430+2303 reported recently motivates us to investigate an imminent binary of supermassive black holes as potential gravitational wave source, the radiated gravitational waves at the end of the merger are shown to be in the band of space-borne detectors. We provide a general analysis on the required detecting sensitivity needed for probing such type gravitational wave sources and make a full discussion by considering two typically designed configurations of space-borne antennas. If a source is so close, it is possible to be detected with Taiji pathfinder-plus which is proposed to be an extension for the planned Taiji pathfinder by just adding an additional satellite to the initial two satellites. The gravitational wave detection on such kind of source enables us to explore the properties of supermassive black holes and the nature of gravity.

The space-based laser interferometers, LISA, Taiji and TianQin, are targeting to observe milliHz gravitational waves (GWs) in the 2030s. The joint observations from multiple space-based detectors yield significant advantages. In this work, we recap the studies and investigations for the joint space-based GW detector networks to highlight: 1) the high precision of sky localization for the massive binary black hole (BBH) coalescences and the GW sirens in the cosmological implication, 2) the effectiveness to test the parity violation in the stochastic GW background observations, 3) the efficiency of subtracting galactic foreground, 4) the improvement in stellar-mass BBH observations. We inspect alternative networks by trading off massive BBH observations and stochastic GW background observation.

Different astrophysical methods can be combined to detect possible deviations from General Relativity. In this work, we consider a class of $f(R)$ gravity models selected by the existence of Noether symmetries. In this framework, it is possible to determine a set of static and spherically symmetric black hole solutions, encompassing small departures from the Schwarzschild geometry. In particular, when gravity is the only dominating interaction, we exploit the ray-tracing technique to reconstruct the image of a black hole, the epicyclic frequencies, and the black hole shadow profile. Moreover, when matter dynamics is also affected by an electromagnetic radiation force, we take into account the general relativistic Poynting-Robertson effect. In light of the obtained results, the proposed strategy results to be robust and efficient: on the one hand, it allows to investigate gravity from strong to weak field regimes; on the other hand, it is capable of detecting small departures from General Relativity, depending on the current observational sensitivity.

Binary neutron stars (BNS) and neutron star-black hole (NSBH) binaries are one of the most promising gravitational wave (GW) sources to probe matter effects. Upcoming observing runs of LIGO-Virgo-KAGRA detectors and future third generation detectors like Einstein Telescope and Cosmic Explorer will allow the extraction of detailed information on these matter effects from the GW signature of BNS and NSBH systems. One subtle effect which may be helpful to extract more information from the detection of compact binary systems is the non-linear memory. In this work, we investigate the observational consequences of gravitational wave non-linear memory in the presence of matter effects. We start by quantifying the impact of non-linear memory on distinguishing BNS mergers from binary black holes (BBH) or NSBH mergers. We find that for the third generation detectors, the addition of non-linear memory to the GW signal model expands the parameter space where BNS signals become distinguishable from the BBH and NSBH signals. Using numerical relativity simulations, we also study the non-linear memory generated from the post-merger phase of BNS systems. We find that it does not show a strong dependence on the equation of state of the NS. However, the non-linear memory from the BNS post-merger phase is much lower than the one from BBH systems of the same masses. Furthermore, we compute the detection prospects of non-linear memory from the post-merger phase of NS systems by accumulating signal strength from a population of BNS mergers for the current and future detectors. Finally, we discuss the impact of possible linear memory from the dynamical ejecta of BNS and NSBH systems and its signal strength relative to the non-linear memory. We find that linear memory almost always has a much weaker effect than non-linear memory.

This paper responds to suggestions that the standard approach to collective neutrino oscillations leaves out potentially important quantum many-body correlations. Arguments in favor of this idea have been based on calculations that, on close scrutiny, offer no evidence either way. Inadequacies of the usual quantum-kinetic formalism are not currently supported by the literature.