Papers for Tuesday, Apr 11 2023

Reconstructing the initial conditions of the universe is a key problem in cosmology. Methods based on simulating the forward evolution of the universe have provided a way to infer initial conditions consistent with present-day observations. However, due to the high complexity of the inference problem, these methods either fail to sample a distribution of possible initial density fields or require significant approximations in the simulation model to be tractable, potentially leading to biased results. In this work, we propose the use of score-based generative models to sample realizations of the early universe given present-day observations. We infer the initial density field of full high-resolution dark matter N-body simulations from the present-day density field and verify the quality of produced samples compared to the ground truth based on summary statistics. The proposed method is capable of providing plausible realizations of the early universe density field from the initial conditions posterior distribution marginalized over cosmological parameters and can sample orders of magnitude faster than current state-of-the-art methods.

The shrinking of a binary orbit driven by the interaction with a gaseous circumbinary disc, initially advocated as a potential way to catalyze the binary merger, has been recently debated in the case of geometrically thick (i.e. with $H/R\gtrsim 0.1$) discs. However, a clear consensus is still missing mainly owing to numerical limitations, such as fixed orbit binaries or lack of resolution inside the cavity carved by the binary in its circumbinary disc. In this work, we asses the importance of evolving the binary orbit by means of hydrodynamic simulations performed with the code {\sc gizmo} in meshless-finite-mass mode. In order to model the interaction between equal mass circular binaries and their locally isothermal circumbinary discs, we enforce hyper-Lagrangian resolution inside the cavity. We find that fixing the binary orbit ultimately leads to an overestimate of the gravitational torque that the gas exerts on the binary, and in an underestimate of the torque due to the accretion of material onto the binary components. Furthermore, we find that the modulation of the accretion rate on the binary orbital period is strongly suppressed in the fixed orbit simulation, while it is clearly present in the live binary simulations. This has potential implications for the prediction of the observable periodicities in massive black hole binary candidates.

The All-Sky Automated Survey for Supernovae (ASAS-SN) began observing in late-2011 and has been imaging the entire sky with nightly cadence since late 2017. A core goal of ASAS-SN is to release as much useful data as possible to the community. Working towards this goal, in 2017 the first ASAS-SN Sky Patrol was established as a tool for the community to obtain light curves from our data with no preselection of targets. Then, in 2020 we released static V-band photometry from 2013--2018 for 61 million sources. Here we describe the next generation ASAS-SN Sky Patrol, Version 2.0, which represents a major progression of this effort. Sky Patrol 2.0 provides continuously updated light curves for 111 million targets derived from numerous external catalogs of stars, galaxies, and solar system objects. We are generally able to serve photometry data within an hour of observation. Moreover, with a novel database architecture, the catalogs and light curves can be queried at unparalleled speed, returning thousands of light curves within seconds. Light curves can be accessed through a web interface (this http URL) or a Python client (https://asas-sn.ifa.hawaii.edu/documentation). The Python client can be used to retrieve up to 1 million light curves, generally limited only by bandwidth. This paper gives an updated overview of our survey, introduces the new Sky Patrol, and describes its system architecture. These results provide significant new capabilities to the community for pursuing multi-messenger and time-domain astronomy.

In gravitational lensing, central images in quads can serve as a powerful probe of the inner regions of lens galaxies. The presence of an offset central supermassive black hole (SMBH) has the potential to distort the time-delay surface in a way such that 3 central images form: a strongly de-magnified image near the SMBH, and two less de-magnified (and potentially observable) images at a central maximum and saddle point. Using a quad lens macro model, we simulate the constraints that could be placed on various lens galaxy parameters based on their central images probability of detection or non-detection. Informed by a recent low-redshift distribution of off-nucleus AGN, we utilize Bayesian inference to constrain the mean SMBH off-nucleus distance and galactic core radius for a sample of 6 quads. In general, we find that a detection of the central image in any quad would favor larger SMBH off-nucleus distances and galaxy core sizes. Assuming a linear relationship between core radii and velocity dispersion $r_c = b\sigma$, these results similarly imply strong constraints on $b$, where the likely case of a central image non-detection in each quad constraining $b$ to $3.11^{+2.72}_{-2.26} \times 10^{-4}$ kpc km$^{-1}$ s. Our results show that tight constraints on lens galaxy parameters can be made regardless of a detection or non-detection of a central image. Therefore, we recommend observational searches for the central image, possibly using our suggested novel detection technique in UV filters, to formalize stronger constraints on lens galaxy parameters.

Using high-resolution Romulus simulations, we explore the origin and evolution of the circumgalactic medium (CGM) in the zone 0.1 $\leq \mathrm{R}/\mathrm{R}_\mathrm{500} \leq$ 1 around massive central galaxies in group-scale halos. We find that the CGM is both multiphase and highly dynamic. Investigating the dynamics, we identify seven patterns of evolution. We show that these are robust and detected consistently across various conditions. There are two pathways by which the gas cools: (1) filamentary cooling inflows and (2) condensations forming from rapidly cooling density perturbations. In our cosmological simulations, the perturbations are mainly seeded by orbiting substructures. We find that condensations can form even when the median $t_\mathrm{cool} / t_\mathrm{ff}$ of the X-ray emitting gas is above the canonical threshold of 10 or 20. Strong amplitude perturbations can provoke runaway cooling regardless of the state of the background gas. We also find perturbations whose local $t_\mathrm{cool} / t_\mathrm{ff}$ ratios drop below the threshold but which do not condense. Rather, the ratios fall to some minimum value and then bounce. These are weak perturbations that are temporarily swept up in satellite wakes and carried to larger radii. Their $t_\mathrm{cool} / t_\mathrm{ff}$ ratios decrease because the denominator ($t_\mathrm{ff}$) is increasing, not because the numerator ($t_\mathrm{cool}$) is decreasing. For structures forming hierarchically, our study highlights the challenge of using a simple threshold argument to infer the CGM's evolution. It also highlights that the median hot gas properties are suboptimal determinants of the CGM's state and dynamics. Realistic CGM models must factor in the effects and after-effects of mergers and orbiting satellites, along with the CGM's heating and cooling cycles.

We introduce Mahakala, a Python-based, modular, radiative ray-tracing code for curved space-times. We employ Google's JAX framework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr metrics. JAX also enables Mahakala to run in parallel on both CPUs and GPUs and achieve speeds comparable to C-based codes. Mahakala natively uses the Cartesian Kerr-Schild coordinate system, which avoids numerical issues caused by the "pole" of spherical coordinates. We demonstrate Mahakala's capabilities by simulating the 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature of Mahakala allows us to easily quantify the relative contribution of different regions of the flow to image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole. We also quantify the relative contribution of the disk, forward jet, and counter jet to 1.3 mm images.

We present the first elemental abundance measurements of the K dwarf (K7V) exoplanet-host star WASP-107 using high-resolution (R = 45,000), near-infrared (H- and K-band) spectra taken from Gemini-S/IGRINS. We use the previously determined physical parameters of the star from the literature and infer the abundances of 15 elements: C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni, all with precision < 0.1 dex, based on model fitting using MARCS model atmospheres and the spectral synthesis code Turbospectrum. Our results show near-solar abundances and a carbon-to-oxygen ratio (C/O) of 0.50 (+/-0.10), consistent with the solar value of 0.54 (+/-0.09). The orbiting planet, WASP-107b, is a super Neptune with a mass in the Neptune regime (= 1.8 M_Nep) and a radius close to Jupiter's (= 0.94 R_Jup). This planet is also being targeted by four JWST Cycle 1 programs in transit and eclipse, which should provide highly precise measurements of atmospheric abundances. This will enable us to properly compare the planetary and stellar chemical abundances, which is essential in understanding the formation mechanisms, internal structure, and chemical composition of exoplanets. Our study is a proof-of-concept that will pave the way for such measurements to be made for all JWST's cooler exoplanet-host stars.

The Event Horizon Telescope (EHT) has produced images of M87* and Sgr A*, and will soon produce time sequences of images, or movies. In anticipation of this, we describe a technique to measure a rotation rate, or pattern speed $\Omega_p$, from movies using an autocorrelation technique. We validate the technique on Gaussian random field models with a known rotation rate and apply it to a library of synthetic images of Sgr A* based on general relativistic magnetohydrodynamics (GRMHD) simulations. We predict that EHT movies will have $\Omega_p \approx 1$ degree per $\mathrm{GMc^{-3}}$, which is of order $15\%$ of the Keplerian orbital frequency in the emitting region. We can plausibly attribute the slow rotation seen in our models to the pattern speed of inward-propagating spiral shocks. We also find that $\Omega_p$ depends strongly on inclination. Application of this technique will enable us to compare future EHT movies with the clockwise rotation of Sgr A* seen in near-infrared flares by GRAVITY. Pattern speed analysis of future EHT observations of M87* and Sgr A* may also provide novel constraints on black hole inclination and spin, as well as an independent measurement of black hole mass.

Characterisation of atmospheres undergoing photo-evaporation is key to understanding the formation, evolution, and diversity of planets. However, only a few upper atmospheres that experience this kind of hydrodynamic escape have been characterised. Our aim is to characterise the upper atmospheres of the hot Jupiters HAT-P-32 b and WASP-69 b, the warm sub-Neptune GJ 1214 b, and the ultra-hot Jupiter WASP-76 b through high-resolution observations of their HeI triplet absorption. In addition, we also reanalyse the warm Neptune GJ 3470 b and the hot Jupiter HD 189733 b. We used a spherically symmetric 1D hydrodynamic model coupled with a non-local thermodynamic equilibrium model. Comparing synthetic absorption spectra with observations, we constrained the main parameters of the upper atmosphere of these planets and classify them according to their hydrodynamic regime. Our results show that HAT-P-32 b photo-evaporates at (130$\pm$70)$\times$10$^{11}$ gs$^{-1}$ with a hot (12 400$\pm$2900 K) upper atmosphere; WASP-69 b loses its atmosphere at (0.9$\pm$0.5)$\times$10$^{11}$ gs$^{-1}$ and 5250$\pm$750 K; and GJ 1214 b, with a relatively cold outflow of 3750$\pm$750 K, photo-evaporates at (1.3$\pm$1.1)$\times$10$^{11}$ gs$^{-1}$. For WASP-76 b, its weak absorption prevents us from constraining its temperature and mass-loss rate significantly; we obtained ranges of 6000-17 000\,K and 23.5$\pm$21.5$\times$10$^{11}$ gs$^{-1}$. Our reanalysis of GJ 3470 b yields colder temperatures, 3400$\pm$350 K, but practically the same mass-loss rate as in our previous results. Our reanalysis of HD 189733 b yields a slightly higher mass-loss rate, (1.4$\pm$0.5)$\times$10$^{11}$ gs$^{-1}$, and temperature, 12 700$\pm$900 K compared to previous estimates. Our results support that photo-evaporated outflows tend to be very light.

We reanalyze near-infrared spectra of the young extrasolar giant planet 51 Eridani b which was originally presented in (Macintosh et al. 2015) and (Rajan et al. 2017) using modern atmospheric models which include a self-consistent treatment of disequilibrium chemistry due to turbulent vertical mixing. In addition, we investigate the possibility that significant opacity from micrometeors or other impactors in the planet's atmosphere may be responsible for shaping the observed spectral energy distribution (SED). We find that disequilibrium chemistry is useful for describing the mid-infrared colors of the planet's spectra, especially in regards to photometric data at M band around 4.5 $\mu$m which is the result of super-equilibrium abundances of carbon monoxide, while the micrometeors are unlikely to play a pivotal role in shaping the SED. The best-fitting, micrometeroid-dust-free, disequilibrium chemistry, patchy cloud model has the following parameters: effective temperature $T_\textrm{eff} = 681$ K with clouds (or without clouds, i.e. the grid temperature $T_\textrm{grid}$ = 900 K), surface gravity $g$ = 1000 m/s$^2$, sedimentation efficiency $f_\textrm{sed}$ = 10, vertical eddy diffusion coefficient $K_\textrm{zz}$ = 10$^3$ cm$^2$/s, cloud hole fraction $f_\textrm{hole}$ = 0.2, and planet radius $R_\textrm{planet}$ = 1.0 R$_\textrm{Jup}$.

The propagation of active galactic nucleus jets depends both on the environment into which they propagate and on their internal structure. To test the impact that different magnetic topologies have on the observable properties of radio galaxies on kpc scales, we conducted a series of magneto-hydrodynamic simulations of jets injected with different magnetic field configurations propagating into a gaseous atmosphere modeled on the Perseus cluster. The simulations show that the structure of the field affects the collimation and propagation of the jets on cluster scales and thus the morphology of the radio lobes inflated by the jets, due both to magnetic collimation and the development of dynamical instabilities in jets with different magnetic topologies. In all cases, the simulations show a distinct reversal of the sychrotron spectral age gradient in the radio lobes about a dynamical time after the jets turn off due to large scale circulation inside the radio lobe, driven primarily by buoyancy, which could provide a way to constrain the age of radio sources in cluster environments without the need for detailed spectral modeling and thus constrain the radio mode feedback efficiency. We suggest a robust diagnostic to search for such age gradients in multi-frequency radio data.

We inspect 4 microlensing events KMT-2021-BLG-1968, KMT-2021-BLG-2010, KMT-2022-BLG-0371, and KMT-2022-BLG-1013, for which the light curves exhibit partially covered short-term central anomalies. We conduct detailed analyses of the events with the aim of revealing the nature of the anomalies. We test various models that can give rise to the anomalies of the individual events including the binary-lens (2L1S) and binary-source (1L2S) interpretations. Under the 2L1S interpretation, we thoroughly inspect the parameter space to check the existence of degenerate solutions, and if they exist, we test the feasibility of resolving the degeneracy. We find that the anomalies in KMT-2021-BLG-2010 and KMT-2022-BLG-1013 are uniquely defined by planetary-lens interpretations with the planet-to-host mass ratios of $q\sim 2.8\times 10^{-3}$ and $\sim 1.6\times 10^{-3}$, respectively. For KMT-2022-BLG-0371, a planetary solution with a mass ratio $q\sim 4\times 10^{-4}$ is strongly favored over the other three degenerate 2L1S solutions with different mass ratios based on the $\chi^2$ and relative proper motion arguments, and a 1L2S solution is clearly ruled out. For KMT-2021-BLG-1968, on the other hand, we find that the anomaly can be explained either by a planetary or a binary-source interpretation, making it difficult to firmly identify the nature of the anomaly. From the Bayesian analyses of the identified planetary events, we estimate that the masses of the planet and host are $(M_{\rm p}/M_{\rm J}, M_{\rm h}/M_\odot) = (1.07^{+1.15}_{-0.68}, 0.37^{+0.40}_{-0.23})$, $(0.26^{+0.13}_{-0.11}, 0.63^{+0.32}_{-0.28})$, and $(0.31^{+0.46}_{-0.16}, 0.18^{+0.28}_{-0.10})$ for KMT-2021-BLG-2010L, KMT-2022-BLG-0371L, and KMT-2022-BLG-1013L, respectively.

The Andromeda galaxy (M31) is our closest neighbouring spiral galaxy, making it an ideal target for studying the physics of the interstellar medium in a galaxy very similar to our own. Using new observations of M31 at 4.76GHz by the C-Band All-Sky Survey (C-BASS), and all available radio data at $1^\circ$ resolution, we produce the integrated spectrum and put new constraints on the synchrotron spectral index and anomalous microwave emission (AME) from M31. We use aperture photometry and spectral modelling to fit for the integrated spectrum of M31, and subtract a comprehensive model of nearby background radio sources. The AME in M31 is detected at $3\sigma$ significance with a peak near 30GHz and flux density $0.27\pm0.09$Jy. The synchrotron spectral index of M31 is flatter than our own Galaxy at $\alpha = -0.66 \pm 0.03$ with no strong evidence of spectral curvature. The emissivity of AME, averaged over the total emission from M31 is lower than typical AME sources in our Galaxy, implying that AME is not uniformly distributed throughout M31 and instead is likely confined to sub-regions -- this will need to be confirmed using future higher resolution observations around 20--30GHz.

To address the problem of the low accuracy of transverse velocity field measurements for small targets in high-resolution solar images, we proposed a novel velocity field measurement method for high-resolution solar images based on PWCNet. This method transforms the transverse velocity field measurements into an optical flow field prediction problem. We evaluated the performance of the proposed method using the Ha and TiO datasets obtained from New Vacuum Solar Telescope (NVST) observations. The experimental results show that our method effectively predicts the optical flow of small targets in images compared with several typical machine- and deep-learning methods. On the Ha dataset, the proposed method improves the image structure similarity from 0.9182 to 0.9587 and reduces the mean of residuals from 24.9931 to 15.2818; on the TiO dataset, the proposed method improves the image structure similarity from 0.9289 to 0.9628 and reduces the mean of residuals from 25.9908 to 17.0194. The optical flow predicted using the proposed method can provide accurate data for the atmospheric motion information of solar images. The code implementing the proposed method is available on https://github.com/lygmsy123/transverse-velocity-field-measurement.

Gaussian Process (GP) has gained much attention in cosmology due to its ability to reconstruct cosmological data in a model-independent manner. In this study, we compare two methods for GP kernel selection: Approximate Bayesian Computation (ABC) Rejection and nested sampling. We analyze three types of data: cosmic Chronometer data (CC), Type Ia Supernovae (SNIa), and Gamma Ray Burst (GRB), using five kernel functions. To evaluate the differences between kernel functions, we assess the strength of evidence using Bayes factors. Our results show that, for ABC Rejection, the Mat\'ern kernel with $\nu$=5/2 (M52 kernel) outperformes the commonly used Radial Basis Function (RBF) kernel in approximating all three datasets. Bayes factors indicate that the M52 kernel typically supports the observed data better than the RBF kernel, but with no clear advantage over other alternatives. However, nested sampling gives different results, with the M52 kernel losing its advantage. Nevertheless, Bayes factors indicate no significant dependence of the data on each kernel.

We examine the time delay between radio and X-ray and between narrow radio frequency flares in Sagittarius A* (Sgr A*), from analyses of the synchrotron, bremsstrahlung and monochromatic luminosity curves. Using the results of 2D relativistic radiation magnetohydrodynamic (MHD) simulations based on the shock oscillation model, we find three types of time delay between the synchrotron and bremsstrahlung emissions: Type A with a time delay of 2 -- 3 h on the shock descending branch, Type B with no time delay and Type C with an inverse time delay of 0.5 -- 1 h on the shock ascending branch. The time delays in Types A and C are interpreted as a transit time of Alfv\'{e}n and acoustic waves between both emission dominant regions, respectively. The delay times between 22 and 43 GHz flares and between 8 and 10 GHz flares are $\sim$ 13 -- 26 min and 13 min, respectively, while the inverse delay also occurs dependently on the shock location branch. These time delays between the narrow radio bands are interpreted as the transit time of the acoustic wave between the frequency-dependent effective radii $R_{\tau_{\rm \nu=1}}$, at which the optical depth $\tau_{\rm \nu}$ at the accretion disc surface becomes $\sim$ unity. The shock oscillation model explains well the observed delay times of 0.5 -- 5 h between radio and X-ray, 20 -- 30 min between 22 and 43 GHz and $\sim$ 18 min between 8 and 10 GHz in Sgr A*.

In this work, we report the presence of rapid intra-night optical variations in both -- flux and polarization of the blazar BL Lacertae during its unprecedented 2020--2021 high state of brightness. The object showed significant flux variability and some color changes, but no firmly detectable time delays between the optical bands. The linear polarization was also highly variable in both -- polarization degree and angle (EVPA). The object was observed from several observatories throughout the world, covering in a total of almost 300 hours during 66 nights. Based on our results, we suggest, that the changing Doppler factor of an ensemble of independent emitting regions, travelling along a curved jet that at some point happens to be closely aligned with the line of sight can successfully reproduce our observations during this outburst. This is one of the most extensive variability studies of the optical polarization of a blazar on intra-night timescales.

Fragmentation is a key step in the process of transforming clouds (and their substructures such as filaments, clumps, and cores) into protostars. The thermal gas pressure and gravitational collapse are believed to be the primary agents governing this process, referred to as the thermal Jeans fragmentation. However, the contributions of other factors (such as magnetic fields and turbulence) to the fragmentation process remain less explored. In this work, we have tested possible fragmentation mechanisms by estimating the mean core mass and mean inter-core separation of the B213 filament. We have used the $\sim$14" resolution James Clerk Maxwell Telescope (JCMT) Submillimetre Common-User Bolometer Array 2 (SCUBA-2)/POL-2 850 $\mu$m dust continuum map and combined it with a Planck 850 $\mu$m map and Herschel data. We find that in addition to the thermal contribution, the presence of ordered magnetic fields could be important in the fragmentation of the B213 filament.

We present ALMA CO(2-1) observations of the Seyfert 2 galaxy NGC 3281 at $\sim$ 100 pc spatial resolution. This galaxy was previously known to present a bi-conical ionised gas outflow extending to 2 kpc from the nucleus. The analysis of the CO moment and channel maps, as well as kinematic modelling reveals two main components in the molecular gas: one rotating in the galaxy plane and another outflowing and extending up to $\sim$ 1.8 -- 2.6 kpc from the nucleus, partially encasing the ionised component. The mass of the outflowing molecular gas component is $M_{\mathrm{mol},\mathrm{out}}$ = $(2.5\pm1.6){\times}10^{6}$ $\rm{M_{\odot}}$, representing $\sim$ 1.7 -- 2 % of the total molecular gas seen in emission within the inner 2.3 kpc. The corresponding mass outflow rate and power are $\dot{M}_{\mathrm{mol},\mathrm{out}}$ = 0.12 -- 0.72 $\rm{M_{\odot} yr^{-1}}$ and $\dot{E}_{\mathrm{mol},\mathrm{out}}$ = (0.045 -- 1.6) ${\times} 10^{40}$ $\rm{erg s^{-1}}$, which translates to a kinetic coupling efficiency with the AGN power of only $10^{-4}$ -- 0.02 %. This value reaches up to 0.1 % when including both the feedback in the ionised and molecular gas, as well as considering that only part of the energy couples kinetically with the gas. Some of the non-rotating CO emission can also be attributed to inflow in the galaxy plane towards the nucleus. The similarity of the CO outflow -- encasing the ionised gas one and the X-ray emission -- to those seen in other sources, suggests that this may be a common property of galactic outflows.

The gamma-ray burst GRB 211211A and its associated kilonova-like emission were reported recently. A significant difference between this association event and GRB 170817A/AT 2017gfo is that GRB 211211A has a very long duration. In this paper, we show that this association event may arise from a neutron star$-$white dwarf (NS$-$WD) merger if a magnetar leaves finally in the central engine. Within the NS$-$WD merger, the main burst of GRB 211211A could be produced by magnetic bubble eruptions from toroidal magnetic field amplification of the pre-merger NS. This toroidal field amplification can be induced by the runaway accretion from the WD debris disc if the disc is in low initial entropy and efficient wind. While the extended emission of GRB 211211A is likely involved with magnetic propelling. The observed energetics and duration of the prompt emission of GRB 211211A can be fulfilled in comparison with those of accretion in hydrodynamical thermonuclear simulation, as long as the WD has a mass $\gtrsim1M_{\odot}$. Moreover, if the X-ray plateau in GRB afterglows is due to the magnetar spin-down radiation, GRB optical afterglows and kilonova-like emission can be well jointly modeled combining the standard forward shock with the radioactive decay power of $^{56}{\rm Ni}$ adding a rotational power input from the post-merger magnetar.

We consider Galileon inflation in the Effective Field Theory (EFT) framework and examine the possibility for PBH formation during slow roll (SR) to ultra slow roll (USR) transitions. We show that loop corrections to the power spectrum, in this case, do not impose additional constraints on the masses of PBHs produced. We indicate that the remarkable non-renormalization property of Galileon due to generalized shift symmetry is responsible for protecting PBH formation from quantum loop corrections.

NGC 4151, a nearby Seyfert galaxy, has recently been reported to emit gamma rays in the GeV range, posing an intriguing astrophysical mystery. The star formation rate of NGC 4151 is too low to explain the observed GeV flux, but the galaxy is known for its coronal activity in X-ray and jet activity in radio. We propose that either the combination of these two activities or the jet activity alone can account for the gamma-ray spectrum. An energy-dependent variability search would allow us to distinguish between the two scenarios, as the coronal component can only contribute at energies of $\lesssim1$ GeV. Our analysis also indicates that the expected neutrino flux from the coronal component is likely to be undetectable by IceCube.

Engine-driven explosions with continuous energy input from the central system have been suggested for supernovae (SNe) associated with a Gamma-Ray Burst (GRB), super-luminous SNe (SLSNe), and at least a fraction of broad-lined SNe Ic (SNe Ic-BL) even without an associated GRB. In the present work, we investigate observational consequences in this scenario, focusing on the case where the energy injection is sufficiently brief, which has been suggested for GRB-SNe. We construct a simplified, spherical ejecta model sequence taking into account the major effects of the central engine; composition mixing, density structure, and the outermost ejecta velocity. Unlike most of the previous works for GRB-SNe, we solve the formation of the photosphere self-consistently, with which we can predict the photometric and spectroscopic observables. We find that these ejecta properties strongly affect their observational appearance in the initial phase (~ a week since the explosion), highlighted by blended lines suffering from higher-velocity absorptions for the flatter density distribution and/or higher outermost ejeca velocity. This behaviour also affects the multi-band light curves in a non-monotonic way. Prompt follow-up observations starting immediately after the explosion thus provides key diagnostics to unveil the nature of the central engine behind GRB-SNe and SNe Ic-BL. For SN 2017iuk associated with GRB 171205A these diagnosing observational data are available, and we show that the expected structure from the engine-driven explosion, i.e., a flat power-law density structure extending up to >~ 100,000 km/s, can explain the observed spectral evolution reasonably well.

This article provides a review on X-ray pulsar-based navigation (XNAV). The review starts with the basic concept of XNAV, and briefly introduces the past, present and future projects concerning XNAV. This paper focuses on the advances of the key techniques supporting XNAV, including the navigation pulsar database, the X-ray detection system, and the pulse time of arrival estimation. Moreover, the methods to improve the estimation performance of XNAV are reviewed. Finally, some remarks on the future development of XNAV are provided.

Electron(positron), proton and nuclei can be accelerated to very high energy by local supernova remnants (SNR). The famous excesses of electron and proton (nuclei) potentially come from such kind of local sources. Recently, the DAMPE experiment measured the electron spectrum (including both electrons and positrons) of cosmic rays with high-accuracy. It provides an opportunity to further explore the excess of electrons. According to the gluon condensation (GC) theory, once GC occurs, huge number of gluons condense at a critical momentum, and the production spectrum of electron and proton showing typical GC characteristics. There are exact correlations between the electron and proton spectrum from a same GC process. It is possible to interpret the power-law break of cosmic rays in view of GC phenomenon, and predict one from another based on the relations between electron and proton spectrum. In this work, we point out the potential existence of a second excess in the electron spectrum, the characteristics of this excess is derived from experimental data of proton. We hope that the future DAMPE experiments will confirm the existence of this second excess and support the result of GC model.

Observations show that there are twisted magnetic flux tubes and plasma flow throughout the solar atmosphere. The main purpose of this work is to obtain the damping rate of sausage modes in the presence of magnetic twist and plasma flow. We obtain the dispersion relation for sausage modes in slow continuity in an inhomogeneous layer under the conditions of magnetic pores, then we solve it numerically. For the selected density profile, the magnetic field, and the plasma flow as a function of radius across the inhomogeneous layer, we show that the effect of the twisted magnetic field on the resonance absorption at low speed of the plasma flow is greater than one at high speed.

Active galactic nuclei (AGN) are believed to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The variability of an AGN's brightness over time can reveal important information about the physical properties of the underlying black hole. The temporal variability is believed to follow a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). With upcoming wide-field surveys set to observe 100 million AGN in multiple bandpass filters, there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are well-suited for modeling AGN time series data, as they can explicitly capture the underlying stochastic dynamics. In this work, we modify latent SDEs to jointly reconstruct the unobserved portions of multivariate AGN light curves and infer their physical properties such as the black hole mass. Our model is trained on a realistic physics-based simulation of ten-year AGN light curves, and we demonstrate its ability to fit AGN light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multi-output Gaussian process regression baseline.

Core-collapse supernovae are among the astrophysical sources of gravitational waves that could be detected by third-generation gravitational-wave detectors. Here, we analyze the gravitational-wave strain signals from two- and three-dimensional simulations of core-collapse supernovae generated using the code F{\sc{ornax}}. A subset of the two-dimensional simulations has non-zero core rotation at the core bounce. A dominant source of time changing quadrupole moment is the $l=2$ fundamental mode ($f-$ mode) oscillation of the proto-neutron star. From the time-frequency spectrogram of the gravitational-wave strain we see that, starting $\sim 400$ ms after the core bounce, most of the power lies within a narrow track that represents the frequency evolution of the $f-$mode oscillations. The $f-$mode frequencies obtained from linear perturbation analysis of the angle-averaged profile of the protoneutron star corroborate what we observe in the spectrograms of the gravitational-wave signal. We explore the measurability of the $f-$mode frequency evolution of protoneutron star for a supernova signal observed in the third-generation gravitational-wave detectors. Measurement of the frequency evolution can reveal information about the masses, radii, and densities of the proto-neutron stars. We find that if the third generation detectors observe a supernova within 10 kpc, we can measure these frequencies to within $\sim$90\% accuracy. We can also measure the energy emitted in the fundamental $f-$mode using the spectrogram data of the strain signal. We find that the energy in the $f-$mode can be measured to within 20\% error for signals observed by Cosmic Explorer using simulations with successful explosion, assuming source distances within 10 kpc.

Cross-correlation spectroscopy is an invaluable tool in the study of exoplanets. However, aliasing between spectral lines makes it vulnerable to systematic biases. This work strives to constrain the aliases of the cross-correlation function to provide increased confidence in the detections of elements in the atmospheres of ultra-hot Jupiters (UHJs) observed with high-resolution spectrographs. We use a combination of archival transit observations of the UHJ KELT-9 b obtained with the HARPS-N and CARMENES spectrographs and show that it is possible to leverage each instrument's strengths to produce robust detections at substantially reduced signal-to-noise. Aliases that become present at low signal-to-noise regimes are constrained through a linear regression model. We confirm previous detections of H I, Na I, Mg I, Ca II, Sc II, Ti II, Cr II, Fe I, and Fe II, and detect eight new species Ca I, Cr I, Ni I, Sr II, Tb II at the 5$\sigma$ level and Ti I, V I, Ba II above the 3$\sigma$ level. Ionised terbium (Tb II) has never before been seen in an exoplanet atmosphere. We further conclude that a 5$\sigma$ threshold may not provide a reliable measure of confidence when used to claim detections, unless the systematics in the cross-correlation function caused by aliases are taken into account.

This paper presents photometric, astrometric, and kinematic analyses of the open clusters NGC 189, NGC 1758 and NGC 7762 based on CCD UBV photometric and Gaia Data Release 3 (DR3) data. According to membership analyses, we identified 32, 57 and 106 most probable member stars with membership probabilities $P\geq 0.5$ in NGC 189, NGC 1758 and NGC 7762, respectively. The color excesses and photometric metallicities of each cluster were determined separately using UBV two-color diagrams. The color excess $E(B-V)$ is $0.590 \pm 0.023$ mag for NGC 189, $0.310 \pm 0.022$ mag for NGC 1758 and $0.640 \pm 0.017$ mag for NGC 7762. The photometric metallicity [Fe/H] is $-0.08 \pm 0.03$ dex for both NGC 189 and NGC 1758, and $-0.12 \pm 0.02$ dex for NGC 7762. Distance moduli and ages of the clusters were obtained by comparing PARSEC isochrones with the color-magnitude diagrams constructed from UBV and Gaia photometric data. During this process, we kept as constant color excess and metallicity for each cluster. The estimated isochrone distance is $1201 \pm 53$ pc for NGC 189, $902 \pm 33$ pc for NGC 1758 and $911 \pm 31$ pc for NGC 7762. These are compatible with the values obtained from trigonometric parallax. Ages of the clusters are $500\pm 50$ Myr, $650\pm 50$ Myr and $2000\pm 200$ Myr for NGC 189, NGC 1758 and NGC 7762, respectively. Galactic orbit integration of the clusters showed that NGC 1758 completely orbits outside the solar circle, while NGC 189 and NGC 7762 enter the solar circle during their orbits.

We present early spectral observations of the very slow Galactic nova Gaia22alz, over its gradual rise to peak brightness that lasted 180 days. During the first 50 days, when the nova was only 3--4 magnitudes above its normal brightness, the spectra showed narrow (FWHM $\approx$ 400 km s$^{-1}$) emission lines of H Balmer, He I, He II, and C IV, but no P Cygni absorption. A few weeks later, the high-excitation He II and C IV lines disappeared, and P Cygni profiles of Balmer, He I, and eventually Fe II lines emerged, yielding a spectrum typical of classical novae before peak. We propose that the early spectra of Gaia22alz are produced in the white dwarf's envelope or accretion disk, reprocessing X-ray and ultraviolet emission from the white dwarf after a dramatic increase in the rate of thermonuclear reactions, during a phase known as the ``early X-ray/UV flash''. If true, this would be one of the rare times that the optical signature of the early X-ray/UV flash has been detected. While this phase might last only a few hours in other novae and thus be easily missed, it was possible to detect in Gaia22alz due to its very slow and gradual rise and thanks to the efficiency of new all-sky surveys in detecting transients on their rise. We also consider alternative scenarios that could explain the early spectral features of Gaia22alz and its unusually slow rise.

A substantial number of ultra-high redshift (8 < z < 17) galaxy candidates have been detected with JWST, posing the question: are these observational results surprising in the context of current galaxy formation models? We address this question using the well-established Santa Cruz semi-analytic models, implemented within merger trees from the new suite of cosmological N-body simulations GUREFT, which were carefully designed for ultra-high redshift studies. Using our fiducial models calibrated at z=0, we present predictions for stellar mass functions, rest-frame UV luminosity functions, and various scaling relations. We find that our (dust-free) models predict galaxy number densities at z~11 (z~ 13) that are an order of magnitude (a factor of ~30) lower than the observational estimates. We estimate the uncertainty in the observed number densities due to cosmic variance, and find that it leads to a fractional error of 30-70% at z=11 (25-150% at z=13) for the available observed fields. We explore which processes in our models are most likely to be rate-limiting for the formation of luminous galaxies at these early epochs, considering the halo formation rate, gas cooling, star formation, and stellar feedback, and conclude that it is mainly efficient stellar-driven winds. We find that a modest boost of a factor of ~4 to the UV luminosities, which could arise from a top-heavy stellar initial mass function characteristic of Pop III stars, would bring our current models into agreement with the observations.

Due to the prevalence of magnetic fields in astrophysical environments, magnetohydrodynamic (MHD) simulation has become a basic tool for studying astrophysical fluid dynamics. To further advance the precision of MHD simulations, we have developed a new simulation code that solves ideal adiabatic or isothermal MHD equations with high-order accuracy. The code is based on the finite-difference weighted essentially non-oscillatory (WENO) scheme and the strong stability-preserving Runge-Kutta (SSPRK) method. Most of all, the code implements a newly developed, high-order constrained transport (CT) algorithm for the divergence-free constraint of magnetic fields, completing its high-order competence. In this paper, we present the version in Cartesian coordinates, which includes a fifth-order WENO and a fourth-order five-stage SSPRK, along with extensive tests. With the new CT algorithm, fifth-order accuracy is achieved in convergence tests involving the damping of MHD waves in three-dimensional space. And substantially improved results are obtained in magnetic loop advection and magnetic reconnection tests, indicating a reduction in numerical diffusivity. In addition, the reliability and robustness of the code, along with its high accuracy, are demonstrated through several tests involving shocks and complex flows. Furthermore, tests of turbulent flows reveal the advantages of high-order accuracy, and show the adiabatic and isothermal codes have similar accuracy. With its high-order accuracy, our new code would provide a valuable tool for studying a wide range of astrophysical phenomena that involve MHD processes.

We present the X-ray spectral and temporal analysis of the black hole X-ray transient Swift J1658.2--4242 observed by {\it AstroSat}. Three epochs of data have been analysed using the JeTCAF model to estimate the mass accretion rates and to understand the geometry of the flow. The best-fit disc mass accretion rate ($\dot m_d$) varies between $0.90^{+0.02}_{-0.01}$ to $1.09^{+0.04}_{-0.03}$ $\dot M_{\rm Edd}$ in these observations, while the halo mass accretion rate changes from $0.15^{+0.01}_{-0.01}$ to $0.25^{+0.02}_{-0.01}$ $\dot M_{\rm Edd}$. We estimate the size of the dynamic corona, that varies substantially from $64.9^{+3.9}_{-3.1}$ to $34.5^{+2.0}_{-1.5}$ $r_g$ and a moderately high jet/outflow collimation factor stipulates isotropic outflow. The inferred high disc mass accretion rate and bigger corona size indicate that the source might be in the intermediate to soft spectral state of black hole X-ray binaries. The mass of the black hole estimated from different model combinations is $\sim 14 M_\odot$. In addition, we compute the quasi-periodic oscillation (QPO) frequencies from the model-fitted parameters, which match the observed QPOs. We further calculate the binary parameters of the system from the decay profile of the light curve and the spectral parameters. The estimated orbital period of the system is $4.0\pm0.4$ hr by assuming the companion as a mid or late K-type star. Our analysis using the JeTCAF model sheds light on the physical origin of the spectro-temporal behaviour of the source, and the observed properties are mainly due to the change in both the mass accretion rates and absorbing column density.

Modern lunar-planetary ephemerides are numerically integrated on the observational timespan of more than 100 years (with the last 20 years having very precise astrometrical data). On such long timespans, not only finite difference approximation errors, but also the accumulating arithmetic roundoff errors become important because they exceed random errors of high-precision range observables of Moon, Mars, and Mercury. One way to tackle this problem is using extended-precision arithmetics available on x86 processors. Noting the drawbacks of this approach, we propose an alternative: using double-double arithmetics where appropriate. This will allow to use only double precision floating-point primitives which have ubiquitous support.

Jupiter's weather layer exhibits long-term and quasi-periodic cycles of meteorological activity that can completely change the appearance of its belts and zones. There are cycles with intervals from 4 to 9 years, dependent on the latitude, which were detected in 5$\mu$m radiation, which provides a window into the cloud-forming regions of the troposphere; however, the origin of these cycles has been a mystery. Here we propose that magnetic torsional oscillations/waves arising from the dynamo region could modulate the heat transport and hence be ultimately responsible for the variability of the tropospheric banding. These axisymmetric waves are magnetohydrodynamic waves influenced by the rapid rotation, which have been detected in Earth's core, and have been recently suggested to exist in Jupiter by the observation of magnetic secular variations by Juno. Using the magnetic field model JRM33, together with the density distribution model, we compute the expected speed of these waves. For the waves excited by variations in the zonal jet flows, their wavelength can be estimated from the width of the alternating jets, yielding waves with a half period of 3.2-4.7 years in 14-23$^\circ$N, consistent with the intervals with the cycles of variability of Jupiter's North Equatorial Belt and North Temperate Belt identified in the visible and infrared observations. The nature of these waves, including the wave speed and the wavelength, is revealed by a data-driven technique, dynamic mode decomposition, applied to the spatio-temporal data for 5$\mu$m emission. Our results imply that exploration of these magnetohydrodynamic waves may provide a new window to the origins of quasi-periodic patterns in Jupiter's tropospheric clouds and to the internal dynamics and the dynamo of Jupiter.

One of the most important questions in cosmology is concerning the fundamental nature of dark matter (DM). DM could consist of spinless particles of very small mass i.e. $m \sim 10^{-22}\ \text{eV}$. This kind of ultralight dark matter (ULDM) would form cored density profiles (called "solitons") at the centre of galaxies. In this context, recently it has been argued that (a) there exists a power law relation between the mass of the soliton and mass of the surrounding halo called the Soliton-Halo (SH) relation, and, (b) the requirement of satisfying observed galactic rotation curves as well as SH relations is so stringent that ULDM is disfavoured from comprising $100\%$ of the total cosmological dark matter. In this work, we revisit these constraints for ULDM particles with non-negligible quartic self-interactions. Using a recently obtained soliton-halo relation which takes into account the effect of self-interactions, we present evidence which suggests that, for $m = 10^{-22}\ \text{eV}$, the requirement of satisfying both galactic rotation curves as well as SH relations can be fulfilled with repulsive self-coupling $\lambda \sim \mathcal{O}(10^{-90})$.

We select a disk-like galaxy sample with observations of the $HI$, $H_{2}$ and dust from Herschel Reference Survey (HRS), and derive inner HI masses within the optical radius. We find that the inner gas-to-dust ratio is almost independent of gas-phase metallicity, and confirm that the inner gas mass ($HI$+$H_{2}$) shows tighter relationship with dust mass and monochromatic 500 $\mu m$ luminosity than the integral gas mass. It supports that dust is more closely associated with co-spatial cold gas than the overall cold gas. Based on the newly calibrated relationship between inner gas mass and dust mass, we predict dust masses for disk-dominated galaxies from the xCOLD GASS sample. The predicted dust masses show scaling relations consistent with fiducial ones in the literature, supporting their robustness. Additionally, we find that at a given dust mass and star formation rate (SFR), the galactic WISE W3 luminosities show significant dependence on the [NII] luminosity and the stellar mass surface density. Such dependence highlights the caveat of using the W3 luminosity as integral SFR indicator, and is consistent with findings of studies which target star-forming regions in more nearby galaxies and accurately derive dust masses based on mapping-mode spectroscopy.

Context. Massive substellar companions orbiting active low-mass stars are rare. They, however, offer an excellent opportunity to study the main mechanisms involved in the formation and evolution of substellar objects. Aims. We aim to unravel the physical nature of the transit signal observed by the TESS space mission on the active M dwarf TOI-5375. Methods. We analysed the available TESS photometric data as well as high-resolution (R $\sim$ 115000) HARPS-N spectra. We combined these data to characterise the star TOI-5375 and to disentangle signals related to stellar activity from the companion transit signal in the light-curve data. We ran an MCMC analysis to derive the orbital solution and apply state-of-the-art Gaussian process regression to deal with the stellar activity signal. Results. We reveal the presence of a companion in the brown dwarf / very-low-mass star boundary orbiting around the star TOI-5375. The best-fit model corresponds to a companion with an orbital period of 1.721564 $\pm$ 10$^{\rm -6}$ d, a mass of 77 $\pm$ 8 $M_{\rm J}$ and a radius of 0.99 $\pm$ 0.16 $R_{\rm J}$. We derive a rotation period for the host star of 1.9692 $\pm$ 0.0004 d, and we conclude that the star is very close to synchronising its rotation with the orbital period of the companion.

Kuiper belt objects, such as Arrokoth, the probable progenitors of short-period comets, formed and evolved at large heliocentric distances, where the ambient temperatures appear to be sufficiently low for preserving volatile ices. By detailed numerical simulations, we follow the long-term evolution of small bodies, composed of amorphous water ice, dust, and ices of other volatile species that are commonly observed in comets. The heat sources are solar radiation and the decay of short-lived radionuclides. The bodies are highly porous and gases released in the interior flow through the porous medium. The most volatile ices, CO and CH$_4$ , are found to be depleted down to the center over a time scale on the order of 100 Myr. Sublimation fronts advance from the surface inward, and when the temperature in the inner part rises sufficiently, bulk sublimation throughout the interior reduces gradually the volatile ices content until they are completely lost. All the other ices survive, which is compatible with data collected by New Horizons on Arrokoth, showing the presence of methanol, and possibly, H$_2$O, CO$_2$, NH$_3$ and C$_2$H$_6$, but no hypervolatiles. The effect of short-lived radionuclides is to increase the sublimation equilibrium temperatures and reduce volatile depletion times. We consider the effect of the bulk density, abundance ratios and heliocentric distance. At 100~au, CO is depleted, but CH$_4$ survives to present time, except for a thin outer layer. Since CO is abundantly detected in comets, we conclude that the source of highly volatile species in active comets must be gas trapped in amorphous ice.

The precise estimation of the statistical errors and accurate removal of the systematical errors are the two major challenges for the stage IV cosmic shear surveys. We explore their impact for the China Space-Station Telescope (CSST) with survey area $\sim17,500\deg^2$ up to redshift $\sim4$. We consider statistical error contributed from Gaussian covariance, connected non-Gaussian covariance and super-sample covariance. We find the super-sample covariance can largely reduce the signal-to-noise of the two-point statistics for CSST, leading to a $\sim1/3$ loss in the figure-of-merit for the matter clustering properties ($\sigma_8-\Omega_m$ plane) and $1/6$ in the dark energy equation-of-state ($w_0-w_a$ plane). We further put requirements of systematics-mitigation on: intrinsic alignment of galaxies, baryonic feedback, shear multiplicative bias, and bias in the redshift distribution, for an unbiased cosmology. The $10^{-2}$ to $10^{-3}$ level requirements emphasize strong needs in related studies, to support future model selections and the associated priors for the nuisance parameters.

High mass accretion rate onto strongly magnetised neutron stars results in the appearance of accretion columns supported by the radiation pressure and confined by the strong magnetic field of a star. At mass accretion rates above $\sim 10^{19}\,{\rm g\,s^{-1}}$, accretion columns are expected to be advective. Under such conditions, a noticeable part of the total energy release can be carried away by neutrinos of a MeV energy range. Relying on a simple model of the neutrino luminosity of accreting strongly magnetised neutron stars, we estimate the neutrino energy fluxes expected from six ULX pulsars known up to date and three brightest Be X-ray transits hosting magnetised neutron stars. Despite the large neutrino luminosity expected in ULX pulsars, the neutrino energy flux from the Be X-ray transients of our Galaxy, SMC and LMC is dominant. However, the neutrino flux from the brightest X-ray transients is estimated to be below the isotropic background by two orders of magnitude at least, which makes impossible direct registration of neutrino emission from accreting strongly magnetised neutron stars nowadays.

Anisotropic diffusion is one of the potential interpretations for the morphology of the Geminga pulsar halo. It interprets the observed slow-diffusion phenomenon through a geometric effect, assuming the mean magnetic field direction around Geminga is closely aligned with the line of sight toward it. However, this direction should not extend further than the correlation length of the turbulent magnetic field $L_c$, which could be $100$ pc or less. We first revisit the $L_c=\infty$ scenario and show that the halo asymmetry predicted by this scenario is mainly contributed by the electrons located beyond the ``core" section around Geminga, which has a length of $100$ pc. Then, considering the directional variation of the magnetic field beyond the core section, we take one magnetic field configuration as an example to investigate the possible halo morphology. The predicted morphology has some different features compared to the $L_c=\infty$ scenario. The current experiments may already be able to test these features. In addition, we use a semi-analytical method to solve the anisotropic propagation equation, which offers significant convenience compared to numerical approaches.

The PRime-focus Infrared Microlensing Experiment (PRIME) will be the first to conduct a dedicated near infrared (NIR) microlensing survey by using a 1.8m telescope with a wide field of view of 1.45 ${\rm deg^{2}}$ at the South African Astronomical Observatory (SAAO). The major goals of the PRIME microlensing survey are to measure the microlensing event rate in the inner Galactic bulge to help design the observing strategy for the exoplanet microlensing survey by the {\it Nancy Grace Roman Space Telescope} and to make a first statistical measurement of exoplanet demographics in the central bulge fields where optical observations are very difficult owing to the high extinction in these fields. Here we conduct a simulation of the PRIME microlensing survey to estimate its planet yields and determine the optimal survey strategy, using a Galactic model optimized for the inner Galactic bulge. In order to maximize the number of planet detections and the range of planet mass, we compare the planet yields among four observation strategies. Assuming {the \citet{2012Natur.481..167C} mass function as modified by \citet{2019ApJS..241....3P}}, we predict that PRIME will detect planetary signals for $42-52$ planets ($1-2$ planets with $M_p \leq 1 M_\oplus$, $22-25$ planets with mass $1 M_\oplus < M_p \leq 100 M_\oplus$, $19-25$ planets $100 M_\oplus < M_p \leq 10000 M_\oplus$), per year depending on the chosen observation strategy.

A minor population of antistars in galaxies has been predicted by some of non-standard models of baryogenesis and nucleosynthesis in the early Universe, and their presence is not yet excluded by the currently available observations. Detection of an unusually high abundance of antinuclei in cosmic rays can probe the baryogenesis scenarios in the early Universe. Recent report of the \textit{AMS-02} collaboration on the tentative detection of a few antihelium nuclei in GeV cosmic rays provided a great hope on the progress in this issue. We discuss possible sources of antinuclei in cosmic rays from antistars which are predicted in a modified Affleck-Dine baryogenesis scenario by Dolgov and Silk (1993). The model allows us to estimate the expected fluxes and isotopic content of antinuclei in the GeV cosmic rays produced in scenarios involving antistars. We show that the flux of antihelium CRs reported by the \textit{AMS-02} experiment can be explained by Galactic anti-nova outbursts, thermonuclear anti-SN Ia explosions, a collection of flaring antistars or an extragalactic source with abundances not violating existing gamma-ray and microlensing constraints on the antistar population.

We have examined inter-night variability of K2-discovered Dippers that are not close to being viewed edge-on, as determined from previously-reported ALMA images, using the SpeX spectrograph and the NASA Infrared Telescope facility (IRTF). The three objects observed were EPIC 203850058, EPIC 205151387, and EPIC 204638512 (2MASS J16042165-2130284). Using the ratio of the fluxes between two successive nights, we find that for EPIC 204638512 and EPIC 205151387, we find that the properties of the dust differ from that seen in the diffuse interstellar medium and denser molecular clouds. However, the grain properties needed to explain the extinction does resemble those used to model the disks of many young stellar objects. The wavelength-dependent extinction models of both EPIC 204638512 and EPIC 205151387 includes grains at least 500 microns in size, but lacks grains smaller than 0.25 microns. The change in extinction during the dips, and the timescale for these variations to occur, imply obscuration by the surface layers of the inner disks. The recent discovery of a highly mis-inclined inner disk in EPIC 204638512 is suggests that the variations in this disk system may point to due to rapid changes in obscuration by the surface layers of its inner disk, and that other face-on Dippers might have similar geometries. The He I line at 1.083 microns in EPIC 205151387 and EPIC 20463851 were seen to change from night to night, suggesting that we are seeing He I gas mixed in with the surface dust.

Long gamma-ray bursts (GRBs), which signify the end-life collapsing of very massive stars, are produced by extremely relativistic jets colliding into circumstellar medium. Huge energy is released both in the first few seconds, namely the internal dissipation phase that powers prompt emissions, and in the subsequent self-similar jet-deceleration phase that produces afterglows observed in broad-band electromagnetic spectrum. However, prompt optical emissions of GRBs have been rarely detected, seriously limiting our understanding of the transition between the two phases. Here we report detection of prompt optical emissions from a gamma-ray burst (i.e. GRB 201223A) using a dedicated telescope array with a high temporal resolution and a wide time coverage. The early phase coincident with prompt {\gamma}-ray emissions show a luminosity in great excess with respect to the extrapolation of {\gamma}-rays, while the later luminosity bump is consistent with onset of the afterglow. The clearly detected transition allows us to differentiate physical processes contributing to early optical emissions and to diagnose the composition of the jet

We present an identification of dust-attenuated star-forming galactic-disk substructures in a typical star-forming galaxy (SFG), UDF2, at $z = 2.696$. To date, substructures containing significant buildup of stellar mass and actively forming stars have yet to be found in typical (i.e., main-sequence) SFGs at $z > 2$. This is due to the strong dust attenuation common in massive galaxies at the epoch and the scarcity of high-resolution, high-sensitivity extinction-independent imaging. To search for disk substructures, we subtracted the central stellar-mass disk from the JWST/NIRCam rest-frame 1.2 $\mu$m image ($0.13''$ resolution) and subtracted, in the visibility plane, the central starburst disk from ALMA rest-frame 240 $\mu$m observations ($0.03''$ resolution). The residual images revealed substructures at rest-frame 1.2 $\mu$m co-located with those found at rest-frame 240 $\mu$m, $\simeq 2$ kpc away from the galactic center. The largest substructure contains $\simeq20$% of the total stellar mass and $\simeq5$% of the total SFR of the galaxy. While UDF2 exhibits a kinematically-ordered velocity field of molecular gas consistent with a secularly evolving disk, more sensitive observations are required to characterize the nature and the origin of this substructure (spiral arms, minor merger, or other types of disk instabilities). UDF2 resides in an overdense region ($N \geqslant 4$ massive galaxies within 70 kpc projected distance at $z=2.690-2.697$) and the substructures may be associated with interaction-induced instabilities. Importantly, a statistical sample of such substructures identified with JWST and ALMA could play a key role in bridging the gap between the bulge-forming starburst and the rest of the galaxy.

The formation of the first supermassive black holes is expected to have occurred in some most pronounced matter and galaxy overdensities in the early universe. We have conducted a sub-mm wavelength continuum survey of 54 $z\sim6$ quasars using the Submillimeter Common-User Bolometre Array-2 (SCUBA2) on the James Clerk Maxwell Telescope (JCMT) to study the environments around $z \sim 6$ quasars. We identified 170 submillimeter galaxies (SMGs) with above 3.5$\sigma$ detections at 450 or 850 \um\, maps. Their FIR luminosities are 2.2 - 6.4 $\times$ 10$^{12} L_{\odot}$, and star formation rates are $\sim$ 400 - 1200 M$_{\odot}$ yr$^{-1}$. We also calculated the SMGs differential and cumulative number counts in a combined area of $\sim$ 620 arcmin$^2$. To a $4\sigma$ detection (at $\sim$ 5.5 mJy), SMGs overdensity is $0.68^{+0.21}_{-0.19}$($\pm0.19$), exceeding the blank field source counts by a factor of 1.68. We find that 13/54 quasars show overdensities (at $\sim$ 5.5 mJy) of $\delta_{SMG}\sim$ 1.5 - 5.4. The combined area of these 13 quasars exceeds the blank field counts with the overdensity to 5.5 mJy of \dsmg $\sim$ $2.46^{+0.64}_{-0.55}$($\pm0.25$) in the regions of $\sim$ 150 arcmin$^2$. However, the excess is insignificant on the bright end (e.g., 7.5 mJy). We also compare results with previous environmental studies of Lyman alpha emitters (LAEs) and Lyman-Break Galaxies (LBGs) on a similar scale. Our survey presents the first systematic study of the environment of quasars at $z\sim6$. The newly discovered SMGs provide essential candidates for follow-up spectroscopic observations to test whether they reside in the same large-scale structures as the quasars and search for protoclusters at an early epoch.

Using both ground-based transit photometry and high-precision radial velocity (RV) spectroscopy, we confirm the planetary nature of TOI-3785 b. This transiting Neptune orbits an M2-Dwarf star with a period of ~4.67 days, a planetary radius of 5.14 +/- 0.16 Earth Radii, a mass of 14.95 +4.10, -3.92 Earth Masses, and a density of 0.61 +0.18, -0.17 g/cm^3. TOI-3785 b belongs to a rare population of Neptunes (4 Earth Radii < Rp < 7 Earth Radii) orbiting cooler, smaller M-dwarf host stars, of which only ~10 have been confirmed. By increasing the number of confirmed planets, TOI-3785 b offers an opportunity to compare similar planets across varying planetary and stellar parameter spaces. Moreover, with a high transmission spectroscopy metric (TSM) of ~150 combined with a relatively cool equilibrium temperature of 582 +/- 16 K and an inactive host star, TOI-3785 b is one of the more promising low-density M-dwarf Neptune targets for atmospheric follow-up. Future investigation into atmospheric mass loss rates of TOI-3785 b may yield new insights into the atmospheric evolution of these low-mass gas planets around M-dwarfs.

The CO-to-H$_2$ conversion factor ($\alpha_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $\alpha_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NGC 3627 and NGC 4321 tracing $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales. Our multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of $\alpha_\rm{CO}$. The central 300 pc nuclei in both galaxies show strong enhancement of temperature $T_\rm{k}>100$ K and density $n_\rm{H_2}>10^3$ cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower $\alpha_\rm{CO}$ than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of $\alpha_\rm{CO}$ with low-J $^{12}$CO optical depths ($\tau_\rm{CO}$), as well as an anti-correlation with $T_\rm{k}$. The $\tau_\rm{CO}$ correlation explains most of the $\alpha_\rm{CO}$ variation in the three galaxy centers, whereas changes in $T_\rm{k}$ influence $\alpha_\rm{CO}$ to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with $\tau_\rm{CO}$ variation in these centers, and thus they are useful observational indicators for $\alpha_\rm{CO}$ variation. We also test current simulation-based $\alpha_\rm{CO}$ prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations.

Wave dark matter (DM) represents a class of the most representative DM candidates. Due to its periodic perturbation to spacetime, the wave DM can be detected with a galactic interferometer - pulsar timing array (PTA). We perform in this Letter a first analysis of applying the $\gamma$-ray PTA to detect the wave DM, with the data of Fermi Large Area Telescope (Fermi-LAT). Despite the limitation in statistics, the $\gamma$-PTA demonstrates a promising sensitivity potential for a mass $\sim 10^{-23}-10^{-22}$ eV. We show that the upper limits not far from those of the dedicated radio-PTA projects can be achieved. Particularly, we have fulfilled an analysis to cross-correlate the pulsar data, which has been essentially missing so far in real data analysis but is known to be crucial for identifying the nature of potential signals, with the Fermi-LAT data of two pulsars.

Feedback likely plays a crucial role in resolving discrepancies between observed and theoretical predictions of dwarf galaxy properties. Stellar feedback was once believed to be sufficient to explain these discrepancies, but it has thus far failed to fully reconcile theory and observations. The recent discovery of energetic galaxy-wide outflows in dwarf galaxies hosting Active Galactic Nuclei (AGN) suggests that AGN feedback may have a larger role in the evolution of dwarf galaxies than previously suspected. In order to assess the relative importance of stellar versus AGN feedback in these galaxies, we perform a detailed Keck/KCWI optical integral field spectroscopic study of a sample of low-redshift star-forming (SF) dwarf galaxies that show outflows in ionized gas in their SDSS spectra. We characterize the outflows and compare them to observations of AGN-driven outflows in dwarfs. We find that SF dwarfs have outflow components that have comparable widths (W$_{80}$) to those of outflows in AGN dwarfs, but are much less blue-shifted, indicating that SF dwarfs have significantly slower outflows than their AGN counterparts. The outflows in SF dwarfs are spatially resolved and significantly more extended than those in AGN dwarfs. The mass loss rates, momentum and energy rates of SF-driven outflows are much lower than those of AGN-driven outflows. Our results indicate that AGN feedback in the form of gas outflows may play an important role in dwarf galaxies and should be considered along with SF feedback in models of dwarf galaxy evolution.

In this work, we study the images of a Kerr black hole (BH) immersed in uniform magnetic fields, illuminated by the synchrotron radiation of charged particles in the jet. We particularly focus on the spontaneously vortical motions (SVMs) of charged particles in the jet region and investigate the polarized images of electromagnetic radiations from the trajectories along SVMs. We notice that there is a critical value $\omega_c$ for charged particle released at a given initial position and subjected an outward force, and once $|qB_0/m|=|\omega_B|>|\omega_c|$ charged particles can move along SVMs in the jet region. We obtain the polarized images of the electromagnetic radiations from the trajectories along SVMs. Our simplified model suggests that the SVM radiations can act as the light source to illuminate the BH and form a photon ring structure.

Supermassive black hole binaries are promising sources of low-frequency gravitational waves (GWs) and bright electromagnetic emission. Pulsar timing array searches for resolved binaries are complex and computationally expensive and so far limited to only a few sources. We present an efficient approximation that empowers large-scale targeted multi-messenger searches by neglecting GW signal components from the pulsar term. This Earth-term approximation provides similar constraints on the total mass and GW frequency of the binary, yet is $>100$ times more efficient.

We provide a comprehensive discussion of the Everpresent $\Lambda$ cosmological model arising from fundamental principles in causal set theory and unimodular gravity. In this framework the value of the cosmological constant ($\Lambda$) fluctuates, in magnitude and in sign, over cosmic history. At each epoch, $\Lambda$ stays statistically close to the inverse square root of the spacetime volume. Since the latter is of the order of $H^2$ today, this provides a way out of the cosmological constant puzzle without fine tuning. Our discussion includes a review of what is known about the topic as well as new motivations and insights supplementing the original arguments. We also study features of a phenomenological implementation of this model, and investigate the statistics of simulations based on it. Our results show that while the observed values of $H_0$ and $\Omega_\Lambda^0$ are not typical outcomes of the model, they can be achieved through a modest number of simulations. We also confirm some expected features of $\Lambda$ based on this model, such as the fact that it stays statistically close to the value of the total ambient energy density (be it matter or radiation dominated), and that it is likely to change sign roughly every Hubble timescale.

Recently there has been an interesting revival of the idea to use large extra dimensions to address the dark energy problem, exploiting the (true) observation that towers of states with masses split, by $M^2_N = f(N) m^2,$ with $f$ an unbounded function of the integer $N$, sometimes contribute to the vacuum energy only an amount of order $m^D$ in $D$ dimensions. It has been argued that this fact is a consequence of swampland conjectures and may require a departure from Effective Field Theory (EFT) reasoning. We test this claim with calculations for Casimir energies in extra dimensions. We show why the domain of validity for EFTs ensures that the tower spacing scale $m$ is always an upper bound on the UV scale for the lower-energy effective theory; use of an EFT with a cutoff part way up a tower is not a controlled approximation. We highlight the role played by the sometimes-suppressed contributions from towers in extra-dimensional approaches to the cosmological constant problem, old and new, and point out difficulties encountered in exploiting it. We compare recent swampland realizations of these arguments with earlier approaches using standard EFT examples, discussing successes and limitations of both.

Unlike F(R) gravity, pure metric F(T) gravity in the vacuum dominated era, ends up with an imaginary action and is therefore not feasible. This eerie situation may only be circumvented by associating a scalar field, which can also drive inflation in the very early universe. We show that, despite diverse claims, F(T) theory admits Noether symmetry only in the pressure-less dust era in the form F(T) proportional to the nth power of T, n being odd integers. A suitable form of F(T), admitting a viable Friedmann-like radiation dominated era, together with early deceleration and late-time accelerated expansion in the pressure-less dust era, has been proposed.

In General Relativity, the gravitational field of an electrically charged, non-rotating, spherically symmetric body is described by the Reissner-Nordstrom (RN) metric. In the naked-singularity regime, a general property of this metric is the existence of a radius, known as the zero-gravity radius, where a test particle would remain at rest. As a consequence of repulsive gravity there is no circular orbit inside this radius. A part of any quasi-stable structure must necessarily lie outside of it. Assuming the compact source Sgr A* at the galactic center may be a naked singularity in RN metric, we provide constraints on the electric charge-to-mass ratio Q/M based on different observations. The compari12 son between the Event Horizon Telescope (EHT) observations and the space-time zero-gravity radius provides the most conservative limit on the charge of Sgr A* to be Q/M < 2.32. Therefore, a charged naked singularity respecting this charge-to-mass constraint is indeed consistent with the current EHT observations.

Atomic spectroscopy is used to search for the space-time variation of fundamental constants which may be due to an interaction with scalar and pseudo-scalar (axion) dark matter. In this letter, we study the effects which are produced by the variation of the nuclear radius and electric quadrupole moment. The sensitivity of the electric quadrupole hyperfine structure to both the variation of the quark mass and the effects of dark matter exceeds that of the magnetic hyperfine structure by 1-2 orders of magnitude. Therefore, the measurement of the variation of the ratio of the electric quadrupole and magnetic dipole hyperfine constants is proposed. The sensitivity of the optical clock transitions in the Yb$^+$ ion to the variation of the nuclear radius allows us to extract, from experimental data, limits on the variation of the hadron and quark masses, the QCD parameter $\theta$ and the interaction with axion and scalar dark matter.

The aim of this paper is to highlight the challenges and potential gains surrounding a coherent description of physics from the high-energy scales of inflation down to the lower energy scales probed in particle-physics experiments. As an example, we revisit the way inflation can be realised within an effective Minimal Supersymmetric Standard Model (eMSSM), in which the $LLe$ and $udd$ flat directions are lifted by the combined effect of soft-supersymmetric-breaking masses already present in the MSSM, together with the addition of effective non-renormalizable operators. We clarify some features of the model and address the question of the one-loop Renormalization Group improvement of the inflationary potential, discussing its impact on the fine-tuning of the model. We also compare the parameter space that is compatible with current observations (in particular the amplitude, $A_{\scriptscriptstyle{\mathrm{S}}}$, and the spectral index, $n_{\scriptscriptstyle{\mathrm{S}}}$, of the primordial cosmological fluctuations) at tree level and at one loop, and discuss the role of reheating. Finally we perform combined fits of particle and cosmological observables (mainly $A_{\scriptscriptstyle{\mathrm{S}}}$, $n_{\scriptscriptstyle{\mathrm{S}}}$, the Higgs mass, and the cold-dark-matter energy density) with the one-loop inflationary potential applied to some examples of dark-matter annihilation channels (Higgs-funnel, Higgsinos and A-funnel), and discuss the status of the ensuing MSSM spectra with respect to the LHC searches.

We study the impact of time-dependent solar cycles in the atmospheric neutrino rate at DUNE and Hyper-Kamiokande (HK), focusing in particular on the flux below 1 GeV. Including the effect of neutrino oscillations for the upward-going component that travels through the Earth, we find that across the solar cycle the amplitude of time variation is about $\pm5\%$ at DUNE, and $\pm 1\%$ at HK. At DUNE, the ratio of up/down-going events ranges from 0.45 to 0.85, while at HK, it ranges from 0.75 to 1.5. Over the 11-year solar cycle, we find that the estimated statistical significance for observing time modulation of atmospheric neutrinos is $4.8\sigma$ for DUNE and $2.0\sigma$ for HK. Flux measurements at both DUNE and HK will be important for understanding systematics in the low-energy atmospheric flux as well as for understanding the effect of oscillations in low-energy atmospheric neutrinos.

If multiple thermal weakly interacting massive particle (WIMP) dark matter candidates exist, then their capture and annihilation dynamics inside a massive stars such as Sun could change from conventional method of study. With a simple correction to time evolution of dark matter (DM) number abundance inside the Sun for multiple dark matter candidates, significant changes in DM annihilation flux depending on annihilation, direct detection cross-section, internal conversion and their contribution to relic abundance are reported in present work.