Papers for Tuesday, Apr 19 2022

Metallicity estimates of circumgalactic gas based on absorption line measurements typically require photoionization modeling to account for unseen ionization states. We explore the impact of uncertainties in the extreme ultraviolet background (EUVB) radiation on such metallicity determinations for the z < 1 circumgalactic medium (CGM). In particular, we study how uncertainties in the power-law slope of the EUV radiation, alpha_EUVB, from active galactic nuclei affect metallicity estimates in a sample of 34 absorbers with HI column densities between 15.25 < log (NHI / cm^-2) < 17.25 and measured metal ion column densities. We demonstrate the sensitivity of metallicity estimates to changes in the EUV power-law slope of active galactic nuclei, alpha_EUVB, at low redshift (z<1), showing derived absorber metallicities increase on average by approximately 0.3 dex as the EUV slope is hardened from alpha_EUVB = -2.0 to -1.4. We use Markov Chain Monte Carlo sampling of photoionization models with alpha_EUVB as a free parameter to derive metallicities for these absorbers. The current sample of absorbers does not provide a robust constraint on the slope alpha_EUVB itself; we discuss how future analyses may provide stronger constraints. Marginalizing over the uncertainty in the slope of the background, we find the average uncertainties in the metallicity determinations increase from 0.08 dex to 0.14 dex when switching from a fixed EUVB slope to one that freely varies. Thus, we demonstrate that EUVB uncertainties can be included in ionization models while still allowing for robust metallicity inferences.

We present results from extensive broadband follow-up of GRB 210204A over the period of thirty days. We detect optical flares in the afterglow at 7.6 x 10^5 s and 1.1 x 10^6 s after the burst: the most delayed flaring ever detected in a GRB afterglow. At the source redshift of 0.876, the rest-frame delay is 5.8 x 10^5 s (6.71 d). We investigate possible causes for this flaring and conclude that the most likely cause is a refreshed shock in the jet. The prompt emission of the GRB is within the range of typical long bursts: it shows three disjoint emission episodes, which all follow the typical GRB correlations. This suggests that GRB 210204A might not have any special properties that caused late-time flaring, and the lack of such detections for other afterglows might be resulting from the paucity of late-time observations. Systematic late-time follow-up of a larger sample of GRBs can shed more light on such afterglow behaviour. Further analysis of the GRB 210204A shows that the late time bump in the light curve is highly unlikely due to underlying SNe at redshift (z) = 0.876 and is more likely due to the late time flaring activity. The cause of this variability is not clearly quantifiable due to the lack of multi-band data at late time constraints by the bad weather conditions. The flare of GRB 210204A is the latest flare detected to date.

Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present the first measurement of the evolution of the radio luminosity function (LF) of LERGs out to $z\sim2.5$; this shows relatively mild evolution. We split the LERGs into those hosted by quiescent and star-forming galaxies, finding a new dominant population of LERGs hosted by star-forming galaxies at high redshifts. The incidence of LERGs in quiescent galaxies shows a steep dependence on stellar-mass out to $z \sim1.5$, consistent with local Universe measurements of accretion occurring from cooling of hot gas haloes. The quiescent-LERGs dominate the LFs at $z<1$, showing a strong decline in space density with redshift, tracing that of the available host galaxies, while there is an increase in the characteristic luminosity. The star-forming LERG LF increases with redshift, such that this population dominates the space densities at most radio-luminosities by $z \sim 1$. The incidence of LERGs in star-forming galaxies shows a much weaker stellar-mass dependence, and increases with redshift, suggesting a different fuelling mechanism compared to their quiescent counterparts, potentially associated with the cold gas supply present in the star-forming galaxies.

We present realistic expectations for the number and properties of neutron star binary mergers to be detected as multi-messenger sources during the upcoming fourth observing run (O4) of the LIGO-Virgo-KAGRA gravitational wave (GW) detectors, with the aim of providing guidance for the optimization of observing strategies. Our predictions are based on a population synthesis model which includes the GW signal-to-noise ratio, the kilonova (KN) optical and near-infrared light curves, the relativistic jet gamma-ray burst (GRB) prompt emission peak photon flux, and the afterglow light curves in radio, optical and X-rays. Within our assumptions, the rate of GW events to be confidently detected during O4 is $7.7^{+11.9}_{-5.7}$ yr$^{-1}$, 78% of which will produce a KN, and a lower 52% will also produce a relativistic jet. The typical depth of current optical electromagnetic search and follow up strategies is still sufficient to detect most of the KNe in O4, but only for the first night or two. The prospects for detecting relativistic jet emission are not promising. While closer events (within z<0.02) will likely still have a detectable cocoon shock breakout, most events will have their GRB emission (both prompt and afterglow) missed unless seen under a favorably small viewing angle. This reduces the fraction of events with detectable jets to 2% (prompt emission, serendipitous) and 10% (afterglow, deep radio monitoring), corresponding to detection rates of $0.17^{+0.26}_{-0.13}$ and $0.78^{+1.21}_{-0.58}$ yr$^{-1}$, respectively. When considering a GW sub-threshold search triggered by a GRB detection, our predicted rate of GW+GRB prompt emission detections increases up to a more promising $0.75^{+1.16}_{-0.55}$ yr$^{-1}$.

The analysis of the APOGEE data suggests the existence of a clear distinction between two sequences of disc stars in the [$\alpha$/Fe] vs. [Fe/H] abundance ratio space. We aim to test if the two-infall chemical evolution models designed to reproduce these two sequences in the solar neighbourhood are also capable to predict the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Along with the predicted chemical composition of SSPs born at different Galactic times in the solar vicinity, we provide their maximum vertical height |zmax| above the Galactic plane computed assuming the relation between the vertical action and stellar age in thin disc stars. The predicted vertical distribution of the [Mg/Fe] abundance ratio is in agreement with the one observed combining the APOGEE DR16 data and the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-$\alpha$ sequence stars). Including the high-$\alpha$ disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z| < 2 kpc. However, our model predicts a too flat growth of the |zmax| as a function of [Mg/Fe] for high-$\alpha$ objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are: i) the data sample with |z| < 2 kpc is more likely contaminated by halo stars, causing the median values to be kinematically hotter, ii) external perturbations such as minor mergers could have heated up the disc, and the heating of the orbits cannot be modelled by only scattering processes. Assuming for the data a disc dissection based on chemistry, the observed |zmax| distributions for high-$\alpha$ and low-$\alpha$ sequences are in good agreement with our model predictions if we consider in the calculation the errors in the vertical action estimates.

Ultracool dwarfs represent a significant proportion of stars in the Milky Way,and deep samples of these sources have the potential to constrain the formation history and evolution of low-mass objects in the Galaxy. Until recently, spectral samples have been limited to the local volume (d<100 pc). Here, we analyze a sample of 164 spectroscopically-characterized ultracool dwarfs identified by Aganze et al. (2022) in the Hubble Space Telescope WFC3 Infrared Spectroscopic Parallel (WISP) Survey and 3D-HST. We model the observed luminosity function using population simulations to place constraints on scaleheights, vertical velocity dispersions and population ages as a function of spectral type}. Our star counts are consistent with a power-law mass function and constant star formation history for ultracool dwarfs, with vertical scaleheights 249$_{-61}^{+48}$ pc for late M dwarfs, 153$_{-30}^{+56}$ pc for L dwarfs, and 175$_{-56}^{+149}$ pc for T dwarfs. Using spatial and velocity dispersion relations, these scaleheights correspond to disk population ages of 3.6$_{-1.0}^{+0.8}$ for late M dwarfs, 2.1$_{-0.5}^{+0.9}$ Gyr for L dwarfs, and 2.4$_{-0.8}^{+2.4}$ Gyr for T dwarfs, which are consistent with prior simulations that predict that L-type dwarfs are on average a younger and less dispersed population. There is an additional 1-2 Gyr systematic uncertainty on these ages due to variances in age-velocity relations. We use our population simulations to predict the UCD yield in the JWST PASSAGES survey, a similar and deeper survey to WISPS and 3D-HST, and find that it will produce a comparably-sized UCD sample, albeit dominated by thick disk and halo sources.

We present new period-$\phi_{31}$-[Fe/H] relations for first overtone RRL stars (RRc), calibrated over a broad range of metallicities ($-2.5 < \textrm{[Fe/H]}< 0.0$) utilizing the largest currently available set of Galactic halo field RRL with homogeneous spectroscopic metallicities. Our relations are defined in the optical (ASAS-SN $V$-band) and, inaugurally, in the infrared (WISE $W1$ and $W2$ bands). Our $V$-band relation can reproduce individual RRc spectroscopic metallicities with a dispersion of 0.30 dex over the entire metallicity range of our calibrator sample (an RMS smaller than what we found for other relations in literature including non-linear terms). Our infrared relation has a similar dispersion in the low and intermediate metallicity range ($\textrm{[Fe/H]} < -0.5$) but tends to underestimate the [Fe/H] abundance around solar metallicity. We tested our relations by measuring both the metallicity of the Sculptor dSph and a sample of Galactic globular clusters, rich in both RRc and RRab stars. The average metallicity we obtain for the combined RRL sample in each cluster is within $\pm 0.08$ dex of their spectroscopic metallicities. The infrared and optical relations presented in this work will enable deriving reliable photometric RRL metallicities in conditions where spectroscopic measurements are not feasible; e.g., in distant galaxies or reddened regions (observed with upcoming Extremely Large Telescopes and the James Webb Space Telescope), or in the large sample of new RRL that will be discovered in large-area time-domain photometric surveys (such as LSST and the Roman space telescope).

The highest critical mass of neutron stars (NSs) was reviewed in the context of equation of state and observationsl results/ It was predicted that the maximum NS mass (MNS) exists in the range MNS ~ 2.2-2.9 M_Sun. However, recent observations of gravitationsl waves and other studies had suggested the higher mass limit of NSs, MNS ~ 3/2 M_Sun. The NS mass up to the value of MNS ~ 2 M_Sun is well understood, and with such a mass value it was meaning ful to discuss the "mass gap" ("m-gap") between the NS and black hole (BH) collapsars.The "m-gap" exists in between the highest mass of NS and the lowest mass of BH collapsars (Mm-gap ~ 2-5 M_Sun). in the mass distribution, the maximim population of NSs and BHs is located at MNS = 1.4 M_Sun and MBH = 6.7 M_Sun, respectively. However, recent ofservational results predicted filling the "m-gap" by the compact objects. In this paper, the concept of gravidynamics was reported to resolve the problem of peak likelihood value of gravitational mass at Mpeak = 6.7 M_Sun and the "m-gap" (Mm-gap ~ 2-5 M_Sun). This concept was based on a non-metric scalar-tensor model of gravitational interaction with localizable field energy. The gravidynamics model shows the total mass (MQ) of a compact relativistic object filled with matter of quark-gluon plasma of the radius r* = GMQ/c2 ~ 10 km, consistent with the "m-gap". It was conceptualized that the total measurable gravitational mass of such an extremely dense object consists of both matter and field, which is described by scalar-tensor components. This model is also useful for predicting the collapsars within the "m-gap".

We explore the possibility of multi-parametric resonances from time varying sound speed during cosmological inflation. In particular, we fix our set-up to the simpler case beyond a single oscillation model already explored in literature: two sinusoidal harmonics around a constant sound speed equal to one. We find that, within the perturbative regime, except for some certain extreme corners of the parameter space, the primordial density spectrum is characterized by two groups of amplified peaks centered around two critical oscillatory frequencies of the sound speed. As a general result, we show that the energy spectrum of the secondary induced GWs from the inflationary era has a single major broad peak, whereas the one from the radiation dominated phase consists of one/two principle peak-like configuration(s) for relatively small/large ratio of two oscillatory frequencies. The GW relic stochastic backgrounds carry a gravitational memory of the parametric resonances during inflation. GW signals from double sound speed resonances can be tested in complementary channels from Pulsar-timing radio-astronomy, space and terrestrial GW interferometers.

We evaluate the potential for gravitational-wave (GW) detection in the frequency band from 10 nHz to 1 $\mu$Hz using extremely high-precision astrometry of a small number of stars. In particular, we argue that non-magnetic, photometrically stable hot white dwarfs (WD) located at $\sim$ kpc distances may be optimal targets for this approach. Previous studies of astrometric GW detection have focused on the potential for less precise surveys of large numbers of stars; our work provides an alternative optimization approach to this problem. Interesting GW sources in this band are expected at characteristic strains around $h_c \sim 10^{-17} \times \left(\mu\text{Hz}/f_{\text{GW}}\right)$. The astrometric angular precision required to see these sources is $\Delta \theta \sim h_c$ after integrating for a time $T \sim 1/f_{\text{GW}}$. We show that jitter in the photometric center of WD of this type due to starspots is bounded to be small enough to permit this high-precision, small-$N$ approach. We discuss possible noise arising from stellar reflex motion induced by orbiting objects and show how it can be mitigated. The only plausible technology able to achieve the requisite astrometric precision is a space-based stellar interferometer. Such a future mission with few-meter-scale collecting dishes and baselines of $\mathcal{O}(100\text{ km})$ is sufficient to achieve the target precision. This collector size is broadly in line with the collectors proposed for some formation-flown, space-based astrometer or optical synthetic-aperature imaging-array concepts proposed for other science reasons. The proposed baseline is however somewhat larger than the km-scale baselines discussed for those concepts, but we see no fundamental technical obstacle to utilizing such baselines. A mission of this type thus also holds the promise of being one of the few ways to access interesting GW sources in this band.

The Solar Ultraviolet Imaging Telescope (SUIT) is an instrument onboard the Aditya-L1 mission of ISRO that will measure and monitor the solar radiation emitted in the near-ultraviolet wavelength range (200-400 nm). SUIT will simultaneously map the photosphere and the chromosphere of the Sun using 11 filters sensitive to different wavelengths and covering different heights in the solar atmosphere and help us understand the processes involved in the transfer of mass and energy from one layer to the other. SUIT will also allow us to measure and monitor spatially resolved solar spectral irradiance that governs the chemistry of oxygen and ozone in the stratosphere of Earth's atmosphere. This is central to our understanding of the Sun climate relationship.

By utilising Gaia EDR3 photometric/astrometric data, we studied the dynamical evolution from the obtained astrophysical, structural and dynamical parameters of the open clusters (OCs), Berkeley\,10 (Be~10), Berkeley\,81 (Be~81), Berkeley\,89 (Be~89), and Ruprecht\,135 (Ru~135). The Gaia EDR3 photometric distances from the isochrone fitting method are smaller than the ones of Gaia EDR2. The relaxation times of four OCs are smaller than their ages, in this regard, they are dynamically relaxed. Their steep overall mass function slopes mean that their low mass stars outnumber their massive ones. Their large $\tau$/relatively small $t_{rlx}$ values imply an advanced mass segregation. Therefore, they seem to have lost their low-mass stars much to the field. Be~89's outer parts indicate an expansion with time. However, Be~10 and Be~81 show the relatively shrinkage core/cluster radii due to dynamical evolution. Ru~135 (0.9~Gyr) may have a primordial origin, instead of shrinking in size and mass with time. Be~89's tidal radius is less than its cluster radius. This means that its member stars lie within its tidal radius, in the sense it is gravitationally bound to the cluster. For the rest OCs, the cluster members beyond their tidal radii are gravitationally unbound to the clusters, which are more influenced by the potential of the Galaxy.

Microflares, one of small-scale solar activities, are believed to be caused by magnetic reconnection. Nevertheless, their three-dimensional (3D) magnetic structures, thermodynamic structures, and physical links to the reconnection have been unclear. In this Letter, based on high-resolution 3D radiative magnetohydrodynamic simulation of the quiet Sun spanning from the upper convection zone to the corona, we investigate 3D magnetic and thermodynamic structures of three homologous microflares. It is found that they originate from localized hot plasma embedded in the chromospheric environment at the height of 2--10 Mm above the photosphere and last for 3--10 minutes with released magnetic energy in the range of $10^{27}-10^{28}$ erg. The heated plasma is almost co-spatial with the regions where the heating rate per particle is maximal. The 3D velocity field reveals a pair of converging flows with velocities of tens of km s$^{-1}$ toward and outflows with velocities of about 100 km s$^{-1}$ moving away from the hot plasma. These features support that magnetic reconnection plays a critical role in heating the localized chromospheric plasma to coronal temperature, giving rise to observed microflares. The magnetic topology analysis further discloses that the reconnection region is located near quasi-separators where both current density and squashing factors are maximal although the specific topology may vary from tether-cutting to fan-spine-like structure.

We present a 3D mapping of the Wolf-Rayet (WR) nebula M1-67 around WR124. We obtained high-resolution San Pedro M\'{a}rtir (SPM) Manchester Echelle Spectrograph (MES) observations along 17 long-slit positions covering all morphological features in M1-67. We are able to unveil the true morphology of M1-67 and its kinematics by interpreting the SPM MES observations by means of the 3D modelling tool for Astrophysics SHAPE. Our SHAPE model that best reproduces the SPM MES data includes three concentric bipolar structures composed by a hollow ellipsoidal structure and a torus. In addition, the model requires the presence of expanding jets and broken blisters in order to reproduce specific spectral features. Our results are consistent with the idea that M1-67 and its progenitor star WR124 have formed through a common envelope scenario that occurred 11.8$^{+4.6}_{-0.8}$ kyr ago. Our bipolar model strongly questions previous suggestions of the presence of a bow shock structure surrounding M1-67. We interpret that the bright structures detected in the spectra extracted from the central regions are produced by wind compression at the receding region of the innermost structure in M1-67. Furthermore, WR124 is moving through a low-density region above the Galactic plane that has negligibly affected the formation history of M1-67.

The origin of strong magnetic fields in white dwarfs has been a puzzle for decades. Recently, a dynamo mechanism operating in rapidly rotating and crystallizing white dwarfs has been suggested to explain the occurrence rates of strong magnetic fields in white dwarfs with close low-mass main sequence star companions. Here we investigate whether the same mechanism may produce strong magnetic fields in close double white dwarfs. The only known strongly magnetic white dwarf that is part of a close double white dwarf system, the magnetic component of NLTT 12758, is rapidly rotating and likely crystallizing and therefore the proposed dynamo mechanism represents an excellent scenario for the origin of its magnetic field. Presenting a revised formation scenario for NLTT 12758, we find a natural explanation for the rapid rotation of the magnetic component. We furthermore show that it is not surprising that strong magnetic fields have not been detected in all other known double white dwarfs. We therefore conclude that the incidence of magnetic fields in close double white dwarfs supports the idea that a rotation and crystallization driven dynamo plays a major role in the generation of strong magnetic fields in white dwarfs.

Time series datasets often have missing or corrupted entries, which need to be ignored in subsequent data analysis. For example, in the context of space physics, calibration issues, satellite telemetry issues, and unexpected events can make parts of a time series unusable. Various approaches exist to tackle this problem, including mean/median imputation, linear interpolation, and autoregressive modeling. Here we study the utility of artificial neural networks (ANNs) to predict statistics, particularly second-order structure functions, of turbulent time series concerning the solar wind. Using a dataset with artificial gaps, a neural network is trained to predict second-order structure functions and then tested on an unseen dataset to quantify its performance. A small feedforward ANN, with only 20 hidden neurons, can predict the large-scale fluctuation amplitudes better than mean imputation or linear interpolation when the percentage of missing data is high. Although, they perform worse than the other methods when it comes to capturing both the shape and fluctuation amplitude together, their performance is better in a statistical sense for large fractions of missing data. Caveats regarding their utility, the optimisation procedure, and potential future improvements are discussed.

We have statistically analyzed 379 radio-loud (RL) CMEs and their associated flares during the period 1996 - 2019 covering both solar cycles (SC) 23 and 24. We classified them into two sets of populations based on the observation period: i) 235 events belong to SC 23 (August 1996 - December 2008) and ii) 144 events belong to SC 24 (January 2009 - December 2019). The average residual acceleration of RL CMEs in SC 24 (--17.39 $\pm$ 43.51 m s$^{-2}$) is two times lower than that of the RL CMEs in SC 23 (--8.29 $\pm$ 36.23 m s$^{-2}$), which means that deceleration of RL CMEs in SC 24 is twice as fast as in SC 23. RL CMEs of SC 23 (1443 $\pm$ 504 km s$^{-1}$; 13.82 $\pm$ 7.40 \emph{R}$_{\circledcirc}$) reach their peak speed at higher altitudes than RL CMEs of SC 24 (1920 $\pm$ 649 km s$^{-1}$; 12.51 $\pm$ 7.41 \emph{R}$_{\circledcirc}$).We also observed that the mean apparent widths of RL CMEs in SC 23 are less than in SC 24which is statistically significant. SC 23 has a lower average CME nose height (3.85 \emph{R}$_{\circledcirc}$) at the start time of DH type II bursts than that of SC 24 (3.46 \emph{R}$_{\circledcirc}$). The starting frequencies of DH type II bursts associated with RL CMEs for SC 24 are significantly larger (formed at lower heights) than that of SC 23. We found that there is a good correlation between the drift rates and the mid-frequencies of DH type II radio bursts for both the solar cycles (\emph{R} = 0.80, $\epsilon$ = 1.53). Most of the RL CMEs kinematics and their associated solar flare properties are found similar for SC 23 and SC 24. We concluded that the reduced total pressure in the heliosphere for SC 24 enables RL CMEs to expand wider and decelerate faster, resulting in DH type II radio emissions at lower heights than SC 23.

Making precise measurements of pulsar dispersion measures (DMs) and applying suitable corrections for them is amongst the major challenges in high-precision timing programmes such as pulsar timing arrays (PTAs). While the advent of wide-band pulsar instrumentation can enable more precise DM measurements and thence improved timing precision, it also necessitates doing careful assessments of frequency-dependent (chromatic) DMs that was theorised by Cordes et al. (2016). Here we report the detection of such an effect in broadband observations of the millisecond pulsar PSR J2241-5236, a high-priority target for current and future PTAs. The observations were made contemporaneously using the wide-band receivers and capabilities now available at the Murchison Widefield Array (MWA), the upgraded Giant Metrewave Radio Telescope (uGMRT) and the Parkes telescopes, thus providing an unprecedentedly large frequency coverage from 80 MHz to 4 GHz. Our analysis shows the measurable changes in DM that scale with the observing frequency ($\nu$) as $\rm \delta DM \sim \nu^{2.5 \pm 0.1}$. We discuss the potential implications of such a frequency dependence in the measured DMs, and the likely impact on the timing noise budget, and comment on the usefulness of low-frequency observations in advancing PTA efforts.

The search for the physical mechanism underlying the observational evidence for the acceleration of the recent universe is a compelling goal of modern fundamental cosmology. Here we quantitatively study a class of homogeneous and isotropic cosmological models in which the matter side of Einstein's equations includes, in addition to the canonical term, a term proportional to the trace of the energy-momentum tensor, $T=\rho-3p$, and constrain these models using low redshift background cosmology data. One may think of these models as extensions of general relativity with a nonlinear matter Lagrangian, and they can be studied either as phenomenological extensions of the standard $\Lambda$CDM model, containing both matter and a cosmological constant, or as direct alternatives to it, where there is no cosmological constant but the additional terms would have to be responsible for accelerating the universe. Overall, our main finding is that parametric extensions of $\Lambda$CDM are tightly constrained, with additional model parameters being constrained to their canonical behaviours to within one standard deviation, while alternative models in this class (which do not have a $\Lambda$CDM limit) are ruled out. This provides some insight on the level of robustness of the $\Lambda$CDM model and on the parameter space still available for phenomenological alternatives and extensions.

Relativistic hadronic plasmas can become under certain conditions supercritical, abruptly and efficiently releasing the energy stored in protons through photon outbursts. Past studies have tried to relate the features of such hadronic supercriticalities (HSC) to the phenomenology of Gamma-Ray Burst (GRB) prompt emission. In this work we investigate, for the first time, HSC in adiabatically expanding sources. We examine the conditions required to trigger HSC, study the role of expansion velocity, and discuss our results in relation to GRB prompt emission. We find multi-pulse light curves from slowly expanding regions ($u_{\rm exp}\lesssim 0.01 c)$ that are a manifestation of the natural HSC quasi-periodicity, while single-pulse light curves with a fast rise and slow decay are found for higher velocities. The formation of the photon spectrum is governed by an in-source electromagnetic cascade. The peak photon energy is $\sim 1$ MeV ($\sim 1$ GeV) for maximum proton energies $\sim 1-10$ PeV ($1-10$ EeV) assuming a jet Lorentz factor 100. Peak $\gamma$-ray luminosities are in the range $10^{49}-10^{52}$ erg s$^{-1}$, with the MeV-peaked spectra being $\sim 100-300$ times more luminous than their GeV-peaked analogues. HSC bursts peaking in the MeV are also copious $\sim 10$ TeV neutrino emitters, with an all-flavour fluence $\sim 10\%$ of the $\gamma$-ray one. The hypothesis that typical long-duration GRBs are powered by HSC could be tested in the near future with more sensitive neutrino telescopes like IceCube-Gen2.

Co-orbital planets (in a $1:1$ mean motion resonance) can be formed within a Laplace resonance chain. Here, we develop a secular model to study the dynamics of the resonance chain $p:p:p+1$, where the co-orbital pair is in a first-order mean motion resonance with the outermost third planet. Our model takes into account tidal dissipation through the use of a Hamiltonian version of the constant time-lag model, which extends the Hamiltonian formalism of the point-mass case. We show the existence of several families of equilibria, and how these equilibria extend to the complete system. In one family, which we call the main branch, a secular resonance between the libration frequency of the co-orbitals and the precession frequency of the pericentres has unexpected dynamical consequences when tidal dissipation is added. We report the existence of two distinct mechanisms that make co-orbital planets much more stable within the $p:p:p+1$ resonance chain rather than outside it. The first one is due to negative real parts of the eigenvalues of the linearised system with tides, in the region of the secular resonance mentioned above. The second one comes from non-linear contributions of the vector field and it is due to eccentricity damping. These two stabilising mechanisms increase the chances of a still-to-come detection of exoplanets in the co-orbital configuration.

During cosmic dawn, the thermal history of the universe is well studied, and a study of this era can give us some of the most useful insight into the universe before the recombination epoch. Its precise modeling and future high-precision measurements will be a valuable tool for determining the thermal history of the universe. In the present work, we study the thermal and ionization history of IGM in the presence of decaying magnetic fields via ambipolar and turbulent decay, Baryon-Dark matter (BDM) interaction, including the DM decay/annihilation. The BDM interaction cross-sections considered are of the form $\sigma=\sigma_0 v^{n}$, where $n=-2$ and $n=-4$. In this work, we show that in the current scenario, the decay/annihilation of the DM particles have a considerable impact on the temperature and ionization histories at low redshift. With the addition of the concept of fractional interaction, which states that if a fraction of the DM particles interacts with the baryons, the temperature and ionization fraction of the baryons show a strong dependence on the percentage of DM particles interacting with the baryons. We have also studied the interesting consequences of the present scenario on the thermal Sunyaev-Zeldovich (tSZ) effect. We show that the highest value of the absolute value of the mean $y-$parameter in the current DM decay/annihilation scenario is well within the values derived from experimental data such as PLANCK, FIRAS, and PIXIE. Later we calculate the bound on the ordinary magnetic fields originating from the Dark photons.

The universe's biggest galaxies have both vast atmospheres and supermassive central black holes. This article reviews how those two components of a large galaxy couple and regulate the galaxy's star formation rate. Models of interactions between a supermassive black hole and the large-scale atmosphere suggest that the energy released as cold gas clouds accrete onto the black hole suspends the atmosphere in a state that is marginally stable to formation of cold clouds. A growing body of observational evidence indicates that many massive galaxies, ranging from the huge central galaxies of galaxy clusters down to our own Milky Way, are close to that marginal state. The gas supply for star formation within a galaxy in such a marginal state is closely tied to the central velocity dispersion (sigma_v) of its stars. We therefore explore the consequences of a model in which energy released during blackhole accretion shuts down star formation when sigma_v exceeds a critical value determined by the galaxy's supernova heating rate.

We report on the cryogenic properties of a low-contraction silicon-aluminum composite, namely Japan Fine Ceramics SA001, to use as a packaging structure for cryogenic silicon devices. SA001 is a silicon--aluminum composite material (75% silicon by volume) and has a low thermal expansion coefficient ($\sim$1/3 that of aluminum). The superconducting transition temperature of SA001 is measured to be 1.18 K, which is in agreement with that of pure aluminum, and is thus available as a superconducting magnetic shield material. The residual resistivity of SA001 is 0.065 $\mathrm{\mu \Omega m}$, which is considerably lower than an equivalent silicon--aluminum composite material. The measured thermal contraction of SA001 immersed in liquid nitrogen is $\frac{L_{293\mathrm{K}}-L_{77\mathrm{K}}}{L_{293\mathrm{K}}}=0.12$%, which is consistent with the expected rate obtained from the volume-weighted mean of the contractions of silicon and aluminum. The machinability of SA001 is also confirmed with a demonstrated fabrication of a conical feedhorn array, with a wall thickness of 100 $\mathrm{\mu m}$. These properties are suitable for packaging applications for large-format superconducting detector devices.

To recover the 21 cm hydrogen line, it is essential to separate the cosmological signal from the much stronger foreground contributions at radio frequencies. The BINGO radio telescope is designed to measure the 21 cm line and detect BAOs using the intensity mapping technique. This work analyses the performance of the GNILC method, combined with a power spectrum debiasing procedure. The method was applied to a simulated BINGO mission, building upon previous work from the collaboration. It compares two different synchrotron emission models and different instrumental configurations, in addition to the combination with ancillary data to optimize both the foreground removal and recovery of the 21 cm signal across the full BINGO frequency band, as well as to determine an optimal number of frequency bands for the signal recovery. We have produced foreground emissions maps using the Planck Sky Model, the cosmological Hi emission maps are generated using the FLASK package and thermal noise maps are created according to the instrumental setup. We apply the GNILC method to the simulated sky maps to separate the Hi plus thermal noise contribution and, through a debiasing procedure, recover an estimate of the noiseless 21 cm power spectrum. We found a near optimal reconstruction of the Hi signal using a 80 bins configuration, which resulted in a power spectrum reconstruction average error over all frequencies of 3%. Furthermore, our tests showed that GNILC is robust against different synchrotron emission models. Finally, adding an extra channel with CBASS foregrounds information, we reduced the estimation error of the 21 cm signal. The optimisation of our previous work, producing a configuration with an optimal number of channels for binning the data, impacts greatly the decisions regarding BINGO hardware configuration before commissioning.

Fast radio bursts (FRBs) are mysterious cosmic sources emitting millisecond-duration radio bursts\cite{Lorimer2007}. Although several hundreds FRBs have been discovered\cite{TheCHIME/FRBCollaboration2021}, their physical nature and central engine are still unclear\cite{Cordes2019,Petroff2019,Zhang2020,Xiao2021}. The variations of Faraday rotation measure and dispersion measure, due to local environment, are crucial clues to understand their physical nature\cite{Michilli2018,Anna-Thomas2022}. The recent observations on the rotation measure of FRB 20201124A show a significant variation on a daytime scale\cite{FAST21}. The variable rotation measure demonstrates that FRB 20201124A is in a dynamic magneto-ionic environment\cite{FAST21}. Intriguingly, the oscillation of rotation measure supports that the local contribution can change sign, which indicates the direction changes of the magnetic field along the line of sight. Here we propose that this FRB resides in a binary system containing a magnetar and a Be star with a decretion disk. When the magnetar approaches the periastron, the propagation of radio waves through the disk of the Be star naturally leads to the observed varying rotation measure, depolarization, large scattering timescale, and Faraday conversion. Searching for FRB signals from Be/X-ray binaries is encouraging.

AliCPT-1 is the first Chinese CMB experiment aiming for high precision measurement of Cosmic Microwave Background B-mode polarization. The telescope, currently under deployment in Tibet, will observe in two frequency bands centered at 90 and 150 GHz. We forecast the CMB lensing reconstruction, lensing-galaxy as well as lensing-CIB (Cosmic Infrared Background) cross correlation signal-to-noise ratio (SNR) for AliCPT-1. We consider two stages with different integrated observation time, namely "4 module*yr" (first stage) and "48 module*yr" (final stage). For lensing reconstruction, we use three different quadratic estimators, namely temperature-only, polarization-only and minimum-variance estimators, using curved sky geometry. We take into account the impact of inhomogeneous hit counts as well as of the mean-field bias due to incomplete sky coverage. In the first stage, our results show that the 150 GHz channel is able to measure the lensing signal at $15\sigma$ significance with the minimum-variance estimator. In the final stage, the measurement significance will increase to $31\sigma$. We also combine the two frequency data in the harmonic domain to optimize the SNR. Our result show that the coadding procedure can significantly reduce the reconstruction bias in the multiple range l>800. Thanks to the high quality of the polarization data in the final stage of AliCPT-1, the EB estimator will dominate the lensing reconstruction in this stage. We also estimate the SNR of cross-correlations between AliCPT-1 CMB lensing and other tracers of the large scale structure of the universe. For its cross-correlation with DESI galaxies/quasars, we report the cross-correlation SNR = 10-20 for the 4 redshift bins at 0.05<z<2.1. In the first stage, the total SNR is about $32$. In the final stage, the lensing-galaxy cross-correlation can reach SNR=52.

Cosmic microwave background (CMB) spectral distortion from Rayleigh scattering is calculated for the first time in rigorous second-order cosmological perturbation theory. The new spectral distortion is sensitive to acoustic dissipation at $10^{-2}<k{\rm Mpc}/h<1$, which slightly extends the scale constrained by the CMB anisotropies. The spectral shape is different from either temperature perturbations or any other traditional spectral distortions from Compton scattering, such as $y$ and $\mu$. The new spectral distortion is not formed in the late Universe, unlike the thermal Sunyaev-Zel'dovich effect degenerated with the primordial $y$ distortions since photons must be hot for Rayleigh scattering. Therefore, ideal measurements can distinguish the signal from the other effects and extract new information during recombination. Assuming cosmological parameters consistent with the recent CMB anisotropy measurements, we find the new spectral distortion is $6.5\times 10^{-3}$Jy/str, which is one order of magnitude smaller than the currently proposed target sensitivity range of voyage 2050.

The interactions between a relativistic magnetized collisionless shock and dense clumps have been expected to play a crucial role on the magnetic field amplification and cosmic-ray acceleration. We investigate this process by two-dimensional Particle-In-Cell (PIC) simulations for the first time, where the clump size is much larger than the gyroradius of downstream particles. We also perform relativistic magnetohydrodynamic (MHD) simulations for the same condition to see the kinetic effects. We find that particles escape from the shocked clump along magnetic field lines in the PIC simulations, so that the vorticity is lower than that in the MHD simulations. Moreover, in both the PIC and MHD simulations, the shocked clump quickly decelerates because of the Lorentz contraction. Owing to the escape and the deceleration, the shocked clump cannot amplify the downstream magnetic field in relativistic collisionless shocks. This large-scale PIC simulation opens a new window to understand large-scale behaviors in collisionless plasma systems.

Transition disks with small inner dust cavities are interesting targets for the study of disk clearing mechanisms. Such disks have been identified through a deficit in the infrared part of their SED, but spatially resolved millimeter imaging is required to confirm the presence of an inner dust cavity. We use high-resolution ALMA observations of 30 mas resolution in Band 6 continuum and $^{12}$CO 2--1 emission of 10 transition disk candidates in the Lupus star forming region, in order to confirm the presence of inner dust cavities and infer the responsible mechanism. The continuum data are analyzed using visibility modeling and the SEDs are compared with radiative transfer models. Out of the six transition disk candidates selected from their SED, only one disk revealed an inner dust cavity of 4 au in radius. Three of the other disks are highly inclined, which limits the detectability of an inner dust cavity but it is also demonstrated to be the possible cause for the infrared deficit in their SED. The two remaining SED-selected disks are very compact, with dust radii of only $\sim$3 au. From the four candidates selected from low-resolution images, three new transition disks with large inner cavities $>$20 au are identified, bringing the total number of transition disks with large cavities in Lupus to 13. SED-selected transition disks with small cavities are biased towards highly inclined and compact disks, which casts doubt on the use of their occurrence rates in estimating dispersal timescales of photoevaporation. Using newly derived disk dust masses and radii, we re-evaluate the size-luminosity and $M_{\rm dust}-M_{\rm star}$ relations. These relations can be understood if the bright disks are dominated by disks with substructure whereas faint disks are dominated by drift-dominated disks. (Abridged)

Theories on the origins of life propose that early cell membranes were synthesized from amphiphilic molecules simpler than phospholipids such as fatty alcohols. The discovery in the interstellar medium (ISM) of ethanolamine, the simplest phospholipid head group, raises the question whether simple amphiphilic molecules are also synthesized in space. We investigate whether precursors of fatty alcohols are present in the ISM. For this, we have carried out a spectral survey at 7, 3, 2 and 1 mm toward the Giant Molecular Cloud G+0.693-0.027 located in the Galactic Center using the IRAM 30m and Yebes 40m telescopes. Here, we report the detection in the ISM of the primary alcohol n-propanol (in both conformers Ga-n-C3H7OH and Aa-n-C3H7OH), a precursor of fatty alcohols. The derived column densities of n-propanol are (5.5+-0.4)x10^13 cm^-2 for the Ga conformer and (3.4+-0.3)x10^13 cm^-2 for the Aa conformer, which imply molecular abundances of (4.1+-0.3)x10^-10 for Ga-n-C3H7OH and of (2.5+-0.2)x10^-10 for Aa-n-C3H7OH. We also searched for the AGa conformer of n-butanol (AGa-n-C4H9OH) without success yielding an upper limit to its abundance of <4.1x10^-11. The inferred CH3OH:C2H5OH:C3H7OH:C4H9OH abundance ratios go as 1:0.04:0.006:<0.0004 toward G+0.693-0.027, i.e. they decrease roughly by one order of magnitude for increasing complexity. We also report the detection of both syn and anti conformers of vinyl alcohol, with column densities of (1.11+-0.08)x10^14 cm^-2 and (1.3+-0.4)x10^13 cm^-2, and abundances of (8.2+-0.6)x10^-10 and (9.6+-3.0)x10^-11, respectively. The detection of n-propanol, together with the recent discovery of ethanolamine in the ISM, opens the possibility that precursors of lipids according to theories of the origin of life, could have been brought to Earth from outer space.

Primordial gravitational waves with extremely low frequency are expected to origin from inflation in the early Universe. The detection of such kind of gravitational waves is of great significance to verify the inflationary theory and determine the energy scale of inflation. Due to the foreground contamination from dust in our Milky Way galaxy, the traditional method using B-mode polarization faces challenges. In this work, we propose an alternative way of detection by investigating the effect of primordial gravitational wave with extremely low frequency on a gravitational lens system with a non-aligned source-deflector-observer configuration. The results show that, with a series of chosen parameters, gravitational lens system with perturbation from extremely low frequency primordial gravitational wave could induce time delay which could deviate from the time delay deduced from theoretical model as much as about sixty percent, meaning that gravitational lens system with a non-aligned configuration could serve as a potential long-base-line detector of extremely low frequency primordial gravitational wave.

Planetary systems are thought to be born in protoplanetary disks. Isotope ratios are a powerful tool for investigating the material origin and evolution from molecular clouds to planetary systems via protoplanetary disks. However, it is challenging to measure the isotope (isotopologue) ratios, especially in protoplanetary disks, because the emission lines of major species are saturated. We developed a new method to overcome these challenges by using optically thin line wings induced by thermal broadening. As a first application of the method, we analyzed two carbon monoxide isotopologue lines, $^{12}$CO $3-2$ and $^{13}$CO $3-2$, from archival observations of a protoplanetary disk around TW Hya with the Atacama Large Millimeter/sub-millimeter Array. The $^{12}$CO/$^{13}$CO ratio was estimated to be ${ 20\pm5}$ at disk radii of ${ 70-110}$ au, which is significantly smaller than the value observed in the local interstellar medium, $\sim69$. It implies that an isotope exchange reaction occurs in a low-temperature environment with $\rm C/O>1$ . In contrast, it is suggested that $^{12}$CO/$^{13}$CO is higher than $\sim{ 84}$ in the outer disk ($r > { 130}$ au), which can be explained by the difference in the binding energy of the isotopologues on dust grains and the CO gas depletion processes. Our results imply that the gas-phase $^{12}$CO/$^{13}$CO can vary by a factor of ${ > 4}$ even inside a protoplanetary disk, and therefore, can be used to trace material evolution in disks.

There exists much uncertainty surrounding interstellar grain-surface chemistry. One of the major reaction mechanisms is grain-surface diffusion for which the the binding energy parameter for each species needs to be known. However, these values vary significantly across the literature which can lead to debate as to whether or not a particular reaction takes place via diffusion. In this work we employ Bayesian inference to use available ice abundances to estimate the reaction rates of the reactions in a chemical network that produces glycine. Using this we estimate the binding energy of a variety of important species in the network, by assuming that the reactions take place via diffusion. We use our understanding of the diffusion mechanism to reduce the dimensionality of the inference problem from 49 to 14, by demonstrating that reactions can be separated into classes. This dimensionality reduction makes the problem computationally feasible. A neural network statistical emulator is used to also help accelerate the Bayesian inference process substantially. The binding energies of most of the diffusive species of interest are found to match some of the disparate literature values, with the exceptions of atomic and diatomic hydrogen. The discrepancies with these two species are related to limitations of the physical and chemical model. However, the use of a dummy reaction of the form H + X -> HX is found to somewhat reduce the discrepancy with the binding energy of atomic hydrogen. Using the inferred binding energies in the full gas-grain version of UCLCHEM results in almost all the molecular abundances being recovered.

A typical inflatable reflector for space application consists of two thin membranes with a parabolic shape. It is critical to understand the interaction of the inflatable and the micrometeoroid environment to which it is exposed. This interaction leads to a series of penetrations of the inflatable membrane on entrance and exit of the impacting particle, creating a pathway for gas escape. To increase the fidelity of the damage expected, we examine the literature for descriptions of micrometeoroid fragmentation and present a theoretical formulation for the damage caused by an impacting particle to the entrance and exit membranes. This theory is compared to an initial set of hyper-velocity tests for micrometeoroid-sized particles on thin film membranes. We use the results of these tests to produce a predictive model. This model is applied to estimate the damage rate near the 1 AU location and output predictions for the effectiveness of a micrometeoroid shield to reduce the damage on the lenticular and effectively optimize its lifetime. Lastly, we apply the kinetic theory of gasses to develop expressions for the expenditure of gas over a specified mission lifetime due to penetrations. Although this paper examines the specific case of an inflated lenticular protected by extra membrane layers, our predictive model can be applied for any gossamer structure composed of polyimide membranes.

Neutrinos are guaranteed to be observable from the next galactic supernova (SN). Optical light and gravitational waves are also observable, but may be difficult to observe if the location of the SN in the galaxy or the details of the explosion are unsuitable. The key to observing the next supernova is to first use neutrinos to understand various physical quantities and then link them to other signals. In this paper, we present Monte Carlo sampling calculations of neutrino events from galactic supernova explosions observed with Super-Kamiokande. The analytical solution of neutrino emission, which represents the long-term evolution of neutrino-light curve from supernovae, is used as a theoretical template. It gives the event rate and event spectrum through inverse beta decay interactions with explicit model parameter dependence. Parameter estimation is performed on these simulated sample data by fitting least squares using the analytical solution. The results show that the mass, radius and total energy of a remnant neutron star produced by a SN can be determined with an accuracy of $\sim 0.1M_\odot$, $\sim 1$ km, and $\sim 10^{51}$ erg, respectively, for a galactic SN at 8 kpc.

Short rise times of Fast Blue Optical Transients (FBOTs) require very light ejected envelopes, $M_{ej} \leq 10^{-1} M_\odot$, much smaller than of a typical supernova. Short peak times also mean that FBOTs should be hydrodynamically, not radioactively powered. The detection by Chandra of X-ray emission in AT2020mrf of $L_X \sim 10^{42} $ erg s$^{-1}$ after 328 days implies total, overall dominant, X-ray energetics at the Gamma Ray Bursts (GRBs) level of $\sim 6 \times 10^{49}$ erg. We further develop a model of Lyutikov & Toonen (2019), whereby FBOTs are the results of a late accretion induced collapse (AIC) of the product of super-Chandrasekhar double white dwarf (WD) merger between ONeMg WD and another WD. Small ejecta mass, and the rarity of FBOTs, result from the competition between mass loss from the merger product to the wind, and ashes added to the core, on time scale of $\sim 10^3-10^4$ years. FBOTs occur only when the envelope mass before AIC is $\leq 10^{-1} M_\odot$. FBOTs proper come from central engine-powered radiation-dominated forward shock as it propagates through ejecta. FBOTs' duration is determined by the diffusion time of photons produced by the NS-driven forward shock within the expanding ejecta. All the photons produced by the central source deep inside the ejecta escape almost simultaneously, producing a short bright event. The high energy emission is generated at the highly relativistic and highly magnetized termination shock, qualitatively similar to Pulsar Wind Nebulae. The X-ray bump observed in AT2020mrf by SRG/eROSITA, predicted by Lyutikov & Toonen (2019), is coming from the break-out of the engine-powered shock from the ejecta into the preceding wind. The model requires total energetics of just few $\times 10^{50}$ ergs, slightly above the observed X-rays. We predict that the system is hydrogen poor.

We report g-r and r-i new colors for 21 Saturn Irregular Satellites, among them, 4 previously unreported. This is the highest number of Saturn Irregular satellites reported in a single survey. These satellites were measured by "stacking" their observations to increase their signal without trailing. This work describes a novel processing algorithm that enables the detection of faint sources under significant background noise and in front of a severely crowded field. Our survey shows these new color measurements of Saturn Irregular Satellites are consistent with other Irregular Satellites populations as found in previous works and reinforcing the observation that the lack of ultra red objects among the irregular satellites is a real feature that separates them from the trans-Neptunian objects (their posited source population).

The stellar halo of the Milky Way is known to have a highly lumpy structure due to the presence of tidal debris and streams accreted from the satellite galaxies. The abundance and distribution of these substructures can provide a wealth of information on the assembly history of the Milky Way. We use some information-theoretic measures to study the anisotropy in a set of Milky Way-sized stellar halos from the Bullock & Johnston suite of simulations. Our analysis shows that the radial anisotropy in each stellar halo increases with the distance from its centre and eventually plateaus out beyond a certain radius. All the stellar halo have a very smooth structure within a radius of $\sim 50$ kpc and a highly anisotropic structure in the outskirts. At a given radius, the polar and azimuthal anisotropies in a stellar halo have two distinct components: (i) a nearly unvarying component from the smooth distribution of stars and (ii) a highly fluctuating component possibly from the substructures. We destroy the substructures and any non-spherical shape of the halo by randomizing the polar and azimuthal coordinates of the stellar particles while keeping their radial distances fixed. We observe that the fluctuating part of the anisotropy is completely eliminated, and the unvarying part of the anisotropy is significantly reduced after the sphericalization. A comparison between the original halo and its sphericalized versions reveals that the unvarying part of the anisotropy originates from the discreteness noise and the non-spherical shape of the halo whereas the substructures contribute to the fluctuating part. We show that such distinction between the anisotropies can constrain the shape of the stellar halo and its substructures. Finally, we map the distribution of the individual substructures by combining the polar and azimuthal anisotropy profiles of the halo at different radial distances.

Odd Radio Circles or ORCs are recently discovered low surface brightness diffuse radio sources, whose progenitors and astrophysical processes responsible for their origins are presently debated. Some ORCs appear to be hosted in distant galaxies and some are host-less. Two plausible scenarios consider ORCs as either nearby supernovae remnants of sizes a few hundred pc in the intragroup medium of the local group of galaxies or alternatively shocked halos of sizes a few hundred kpc around distant galaxies. The input energy for shocks, required to create ORCs of a few hundred kpc size is estimated to be $10^{56} - 10^{59}$ erg. It is shown here that the energy released in shocks with a rate of $10^{41} - 10^{44}$ erg s$^{-1}$ due to multiple ($10^{3} - 10^{6}$) tidal disruption events in a short span of time (a few tens of Myr) in a merger galaxy system, hosting a massive black hole can supply required energy to generate ORCs around distant galaxies. The most plausible and abundant hosts for ORCs can be post-starburst galaxies at intermediate redshifts. The presently observed dominance of tidal disruption events in post-starburst galaxies, redshift evolution and environments of post-starburst galaxies at intermediate redshifts favourably support the current observations.

We present an exact $\gamma=5/3$ spherical accretion model which modifies the Bondi boundary condition of $\rho \to const.$ as $r\to \infty$ to one where $\rho \to 0$ as $r \to \infty$. This change allows for simple power law solutions on the density and infall velocity fields, ranging from a cold empty free-fall condition where pressure tends to zero, to a hot hydrostatic equilibrium limit with no infall velocity. As in the case of the Bondi solution, a maximum accretion rate appears. As in the $\gamma=5/3$ case of the Bondi solution, no sonic radius appears, this time however, because the flow is always characterised by a constant Mach number. This number equals 1 for the case of the maximum accretion rate, diverges towards the cold empty state, and becomes subsonic towards the hydrostatic equilibrium limit. Deviations from sphericity are then explored through an analytic perturbative analysis. The perturbed solution yields a rich phenomenology through density and radial velocity fields in terms of Legendre polynomials, which we begin to explore for simple angular velocity boundary conditions having zeros on the plane and pole. Infall/outflow solutions appear for the sub-sonic parameter region, closely resembling the outputs of recent numerical experiments. The strong density gradients in these cases result in significant pressure gradients which accelerate the polar outflows to many times the local escape velocity, well within the validity range of the perturbative analysis. Our results could complement our understanding of the origin of outflows in a variety of astrophysical settings surrounding infall situations, through purely hydrodynamical physics.

Pulsars are spinning neutron stars with very regular periods. These pulsars have, however, had instances where they exhibit a change in their periods. Older theories have shown that older pulsars have a tendency to skip and speed up. Newer theories have been created, due to the discovery that younger X-ray pulsars exhibit the same skips. The older theories explain that the core of the pulsar is a superfluid with a differential rotation and the core will occasionally exhibit solid properties to catch the crust of the pulsar and speed it up. The newer quantum mechanical theory states that quantum particle packets, called the strange nuggets, slam into the side of the pulsar to add angular momentum to the pulsar and then release it later.

Recent experiments searching for sub-GeV/$c^2$ dark matter have observed event excesses close to their respective energy thresholds. Although specific to the individual technologies, the measured excess event rates have been consistently reported at or below event energies of a few-hundred eV, or with charges of a few electron-hole pairs. In the present work, we operated a 1-gram silicon SuperCDMS-HVeV detector at three voltages across the crystal (0 V, 60 V and 100 V). The 0 V data show an excess of events in the tens of eV region. Despite this event excess, we demonstrate the ability to set a competitive exclusion limit on the spin-independent dark matter--nucleon elastic scattering cross section for dark matter masses of $\mathcal{O}(100)$ MeV/$c^2$, enabled by operation of the detector at 0 V potential and achievement of a very low $\mathcal{O}(10)$ eV threshold for nuclear recoils. Comparing the data acquired at 0 V, 60 V and 100 V potentials across the crystal, we investigated possible sources of the unexpected events observed at low energy. The data indicate that the dominant contribution to the excess is consistent with a hypothesized luminescence from the printed circuit boards used in the detector holder.

Implementation of many statistical methods for large, multivariate data sets requires one to solve a linear system that, depending on the method, is of the dimension of the number of observations or each individual data vector. This is often the limiting factor in scaling the method with data size and complexity. In this paper we illustrate the use of Krylov subspace methods to address this issue in a statistical solution to a source separation problem in cosmology where the data size is prohibitively large for direct solution of the required system. Two distinct approaches are described: one that uses the method of conjugate gradients directly to the Kronecker-structured problem and another that reformulates the system as a Sylvester matrix equation. We show that both approaches produce an accurate solution within an acceptable computation time and with practical memory requirements for the data size that is currently available.

Astronomical interest within the current Colombian territory has its roots in the Botanical Expedition of the New Kingdom of Granada, which stimulated the creation of an astronomical observatory in 1803, the first one established in the New World to pursue systematic observations and meteorological studies. After the death in 1816 of its first director, Francisco Jos\'e de Caldas, during the convulsive independence period, no major astronomical observations were made for decades, with few exceptions. In this work we delve into the contributions of the astronomer Jos\'e Mar\'ia Gonz\'alez Benito, the main reactivator of the National Astronomical Observatory of Colombia in the second half of the 19th century, pointing out his pioneering efforts that put worldwide attention to it, and to his own private observatory making him one of the most committed figures to the development of astronomical sciences in the country and the most renowned Colombian in the international astronomical research scene of his time.

We propose a type II seesaw model for light Dirac neutrinos to provide an explanation for the recently reported anomaly in W boson mass by the CDF collaboration with $7\sigma$ statistical significance. In the minimal model, the required enhancement in W boson mass is obtained at tree level due to the vacuum expectation value of a real scalar triplet, which also plays a role in generating light Dirac neutrino mass. Depending upon the couplings and masses of newly introduced particles, we can have thermally or non-thermally generated relativistic degrees of freedom $\Delta N_{\rm eff}$ in the form of right handed neutrinos which can be observed at future cosmology experiments. Extending the model to a radiative Dirac seesaw scenario can also accommodate dark matter and lepton anomalous magnetic moment.

We investigate two Type-IIa Minimally Modified Gravity theories, namely VCDM and Cuscuton theories. We confirm that all acceptable Cuscuton solutions are always solutions for VCDM theory. However, the inverse does not hold. We find that VCDM allows for the existence of exact General Relativity (GR) solutions with or without the presence of matter fields and a cosmological constant. We determine the conditions of existence for such GR-VCDM solutions in terms of the trace of the extrinsic curvature and on the fields which define the VCDM theory. On the other hand, for the Cuscuton theory, we find that the same set of exact GR solutions (such as Schwarzschild and Kerr spacetimes) is not compatible with timelike configurations of the Cuscuton field and therefore cannot be considered as acceptable solutions. Nonetheless, in Cuscuton theory, there could exist solutions which are not the same but close enough to GR solutions. We also show the conditions to determine intrinsic-VCDM solutions, i.e. solutions which differ from GR and do not belong to the Cuscuton model. We finally show that in cosmology a mapping between VCDM and the Cuscuton is possible, for a generic form of the VCDM potential. In particular, we find that for a quadratic potential in VCDM theory, this mapping is well defined giving an effective redefinition of the Planck mass for the cosmological background solutions of both theories.

We develop a model for an interstellar communication network that is composed of relay nodes that transmit diffraction-limited beams of photons. We provide a multi-dimensional rationale for such a network of communication in lieu of interstellar beacons. We derive a theoretical expression for the bit rate of communication based on fundamental physics, constrained by the energy available for photons and the diffraction of the beam that dilutes the information by the inverse square law. We find that meter-scale probes are severely limited in their bit rate, under 1 Gbps, over distances of a light year. However, that bit rate is proportional to the 4th power of the size of the optics that transmit and receive the photons, and inversely proportional to the square of the distance between them, thus favoring large optics and short separations between nodes. The optimized architecture of interstellar communication consists of a network of nodes separated by sub-light-year distances and strung out between neighboring stars.

The millisecond pulsars, old-recycled objects spinning with high frequency $\mathcal{O}$(kHz) sustaining the deformation from their spherical shape, may emit gravitational-waves (GW). These are one of the potential candidates contributing to the anisotropic stochastic gravitational-wave background (SGWB) observable in the ground-based GW detectors. Here, we present the results from a likelihood-based targeted search for the SGWB due to millisecond pulsars in the Milky Way, by analyzing the data from the first three observing runs of Advanced LIGO and Advanced Virgo detector. We assume that the shape of SGWB power spectra and the sky distribution is known a priori from the population synthesis model. The information of the ensemble source properties, i.e., the in-band number of pulsars, $N_{obs}$ and the averaged ellipticity, $\mu_\epsilon$ is encoded in the maximum likelihood statistic. We do not find significant evidence for the SGWB signal from the considered source population. The best Bayesian upper limit with $95\%$ confidence for the parameters are $N_{obs}\leq8.8\times10^{4}$ and $\mu_\epsilon\leq1.1\times10^{-7}$, which is comparable to the bounds on mean ellipticity with the GW observations of the individual pulsars. Finally, we show that for the plausible case of $N_{obs}=40,000$, with the one year of observations, the one-sigma sensitivity on $\mu_\epsilon$ might reach $10^{-8}$ and $2.7\times10^{-9}$ for the second-generation detector network having A+ sensitivity and third-generation detector network respectively.