Papers for Wednesday, Oct 19 2022

Recent ALMA large surveys unveiled the presence of significant dust continuum emission in star-forming galaxies at $z>4$. Unfortunately, such large programs -- i.e. ALPINE ($z\sim 5$) and REBELS ($z \sim 7$) -- only provide us with a single Far-Infrared (FIR) continuum data point for their individual targets. Therefore, high-$z$ galaxies FIR spectral energy densities (SEDs) remain mostly unconstrained, hinging on an assumption for their dust temperature ($T_{\rm d}$) in the SED fitting procedure. This introduces uncertainties in the inferred dust masses ($M_{\rm d }$), infrared luminosities ($L_{\rm IR}$), and obscured Star Formation Rate (SFR) fraction at $z > 4$. In this work we use a method that allows us to constrain $T_{\rm d}$ with a single band measurement by combining the $158\ \mathrm{\mu m}$ continuum information with the overlying [CII] emission line. We analyse the $21$ [CII] and FIR continuum detected $z\sim 5$ galaxies in ALPINE, finding a range of $T_{\rm d}=25-60\ \mathrm{K}$ and $M_{\rm d} = 0.6-25.1\ \times 10^{7}\ \mathrm{M_{\odot}}$. Given the measured stellar masses of ALPINE galaxies, the inferred dust yields are around $M_{\rm d}/M_{\star} = (0.2-8) \times 10^{-3}$, consistent with theoretical dust-production constraints. We find that $8$ out of $21$ ALPINE galaxies have $L_{\rm IR} \geq 10^{12}\ \mathrm{L_{\odot}}$, comparable to UltraLuminous IR Galaxies (ULIRGs). Relying on ultraviolet-to-optical SED fitting, the SFR was underestimated by up to $2$ orders of magnitude in $4$ of these $8$ ULIRGs-like galaxies. We conclude that these $4$ peculiar sources should be characterised by a two-phase interstellar medium structure with "spatially-segregated" FIR and ultraviolet emitting regions.

Based on their differing radio morphologies, powerful radio galaxies can be separated into the Fanaroff-Riley I (FR-I) and II (FR-II) classes. Hybrid morphology radio sources (HyMoRS) contain morphologies consistent with each type of jet on either side: a powerful, highly relativistic FR-II-like jet terminating in a hotspot on one side and an FRI-like plume on the other. HyMoRS present a unique opportunity to study the conditions which give rise to the dichotomy. Using host galaxy properties, we conduct the first multiwavelength investigation into whether orientation can explain HyMoRS morphology. Through optical spectroscopy and mid-infrared photometry, we analyze the emission characteristics, and evaluate the broad characteristics of five HyMoRS host galaxies at intermediate redshifts (0.4 < z < 1.5). The HyMoRS host galaxies in our sample have properties consistent with typical host galaxies of FR-II sources, suggesting that the observed hybrid morphologies may be caused by a dense, cluster-like environment bending FR-II jets combined with a favorable orientation which can make one side appear similar to an FR-I jet. Our results thus support the hypothesis that HyMoRS are mainly caused by environment and orientation.

Galaxy clusters accrete mass through large scale, strong, structure-formation shocks. Such a virial shock is thought to deposit fractions $\xi_e$ and $\xi_B$ of the thermal energy in cosmic-ray electrons (CREs) and magnetic fields, respectively, thus generating a leptonic virial ring. However, the expected synchrotron signal was not convincingly established until now. We stack low-frequency radio data from the OVRO-LWA around the 44 most massive, high latitude, extended MCXC clusters, enhancing the ring sensitivity by rescaling clusters to their characteristic, $R_{500}$ radii. Both high (73 MHz) and co-added low ($36\text{--}68\text{ MHz}$) frequency channels separately indicate a significant ($4\text{--}5\sigma$) excess peaked at $(2.4 \text{--} 2.6) R_{500}$, coincident with a previously stacked Fermi $\gamma$-ray signal interpreted as inverse-Compton emission from virial-shock CREs. The stacked radio signal is well fit (TS-test: $4$--$6\sigma$ at high frequency, $4$--$8\sigma$ at low frequencies, and $8$--$10\sigma$ joint) by virial-shock synchrotron emission from the more massive clusters, with $\dot{m}\xi_e\xi_B\simeq (1\text{--}4)\times 10^{-4}$, where $\dot{m}\equiv \dot{M}/(MH)$ is the dimensionless accretion rate for a cluster of mass $M$ and a Hubble constant $H$. The inferred CRE spectral index is flat, $p \simeq 2.0 \pm 0.2$, consistent with acceleration in a strong shock. Assuming equipartition or using $\dot{m}\xi_e\sim0.6\%$ inferred from the Fermi signal yields $\xi_B\simeq (2\text{--}9)\%$, corresponding to $B \simeq (0.1\text{--}0.3)~\mu\text{G}$ magnetic fields downstream of typical virial shocks.

Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are under-massive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multi-channel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fueling: constant thresholds in these properties are typically crossed within ~0.1 Hubble time of accelerated BH fueling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold Sigma1 ~ 10^9.5 Msun/kpc^2, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived, thin gas disks, as well as with virialization of the inner circumgalactic medium. The halo mass Mh ~ 10^12 Msun and stellar mass Mstar ~ 10^10.5 Msun at which BH growth accelerates correspond to ~L* galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above ~L*.

We present the Spectroscopic Classification of Astronomical Transients (SCAT) survey, which is dedicated to spectrophotometric observations of transient objects such as supernovae and tidal disruption events. SCAT uses the SuperNova Integral-Field Spectrograph (SNIFS) on the University of Hawai'i 2.2-meter (UH2.2m) telescope. SNIFS was designed specifically for accurate transient spectrophotometry, including absolute flux calibration and host-galaxy removal. We describe the data reduction and calibration pipeline including spectral extraction, telluric correction, atmospheric characterization, nightly photometricity, and spectrophotometric precision. We achieve $\lesssim 5\%$ spectrophotometry across the full optical wavelength range ($3500-9000~\r{A}$) under photometric conditions. The inclusion of photometry from the SNIFS multi-filter mosaic imager allows for decent spectrophotometric calibration ($10-20\%$) even under unfavorable weather/atmospheric conditions. SCAT obtained $\approx 650$ spectra of transients over the first 3 years of operations, including supernovae of all types, active galactic nuclei, cataclysmic variables, and rare transients such as superluminous supernovae and tidal disruption events. These observations will provide the community with benchmark spectrophotometry to constrain the next generation of hydrodynamic and radiative transfer models.

Two recent discoveries, namely PSR J0901-4046 and GLEAM-X J162759.5-523504.3 (hereafter GLEAM-X J1627), have corroborated an extant population of radio-loud periodic sources with long periods (76 s and 1091 s respectively) whose emission can hardly be explained by rotation losses. We argue that GLEAM-X J1627 is a highly-magnetized object consistent with a magnetar (an ultra long period magnetar - ULPM), and demonstrate it is unlikely to be either a magnetically or a rotationally-powered white dwarf. By studying these sources together with previously detected objects, we find there are at least a handful of promising candidates for Galactic ULPMs. The detections of these objects imply a substantial number, $N \gtrsim 12800^{+19000}_{-10100}$ and $N \gtrsim 510^{+500}_{-420}$ for PSR J0901-4046 like and GLEAM-X J1627 like objects, respectively, within our Galaxy. These source densities, as well as cooling age limits from non-detection of thermal X-rays, Galactic offsets, timing stability and dipole spindown limits, all imply the ULPM candidates are substantially older than confirmed galactic magnetars and that their formation channel is a common one. Their existence implies widespread survival of magnetar-like fields for several Myr, distinct from the inferred behaviour in confirmed Galactic magnetars. ULPMs may also constitute a second class of FRB progenitors which could naturally exhibit very long periodic activity windows. Finally, we show that existing radio campaigns are biased against detecting objects like these and discuss strategies for future radio and X-ray surveys to identify more such objects. We estimate that ${\cal O}(100)$ more such objects should be detected with SKA-MID and DSA-2000.

Spectral distortions of the cosmic microwave background (CMB) have been recognized as an important future probe of the early Universe. Existing theoretical studies primarily focused on describing the evolution and creation of average distortions, ignoring spatial perturbations in the plasma. One of the main reasons for this choice is that a treatment of the spectro-spatial evolution of the photon field deep into the primordial Universe requires solving a radiative transfer problem for the distortion signals, which in full detail is computationally challenging. Here we provide the first crucial step towards tackling this problem by formulating a new spectral discretisation of the underlying average thermalisation Green's function. Our approach allows us to convert the high-dimensional partial differential equation system (~1,000-10,000 equations) into and set of ordinary differential equations of much lower dimension (~10 equations). We demonstrate the precision of the approach and highlight how it may be further improved in the future. We also clarify the link of the observable spectral distortion parameters (e.g., mu and y) to the computational spectral basis that we use in our frequency discretisation. This reveals how several basis-dependent ambiguities can be interpreted in future CMB analysis. Even if not exact, the new Green's function discretisation can be used to formulate a generalised photon Boltzmann-hierarchy, which can then be solved with methods that are familiar from theoretical studies of the CMB temperature and polarisation anisotropies. We will carry this program out in a series of companion papers, thereby opening the path to full spectro-spatial exploration of the CMB with future CMB imagers and spectrometers.

Wide, deep, blind continuum surveys at submillimetre/millimetre (submm/mm) wavelengths are required to provide a full inventory of the dusty, distant Universe. However, conducting such surveys to the necessary depth, with sub-arcsec angular resolution, is prohibitively time-consuming, even for the most advanced submm/mm telescopes. Here, we report the most recent results from the ALMACAL project, which exploits the 'free' calibration data from the Atacama Large Millimetre/submillimetre Array (ALMA) to map the lines of sight towards and beyond the ALMA calibrators. ALMACAL has now covered 1,001 calibrators, with a total sky coverage around 0.3 deg2, distributed across the sky accessible from the Atacama desert, and has accumulated more than 1,000h of integration. The depth reached by combining multiple visits to each field makes ALMACAL capable of searching for faint, dusty, star-forming galaxies (DSFGs), with detections at multiple frequencies to constrain the emission mechanism. Based on the most up-to-date ALMACAL database, we report the detection of 186 DSFGs with flux densities down to S870um ~ 0.2mJy, comparable with existing ALMA large surveys but less susceptible to cosmic variance. We report the number counts at five wavelengths between 870um and 3mm, in ALMA bands 3, 4, 5, 6 and 7, providing a benchmark for models of galaxy formation and evolution. By integrating the observed number counts and the best-fitting functions, we also present the resolved fraction of the cosmic infrared background (CIB) and the CIB spectral shape. Combining existing surveys, ALMA has currently resolved about half of the CIB in the submm/mm regime.

We report on the changing-look nature of the active galactic nucleus (AGN) in the galaxy NGC 4156, as serendipitously discovered thanks to data acquired in 2019 at the Telescopio Nazionale Galileo (TNG) during a students' observing programme. Previous optical spectra had never shown any signatures of broad-line emission, and evidence of the AGN had come only from X-ray observations, being the optical narrow-line flux ratios unable to unambiguously denote this galaxy as a Seyfert. Our 2019 TNG data unexpectedly revealed the appearance of broad-line components in both the H$\alpha$ and H$\beta$ profiles, along with a rise of the continuum, thus implying a changing-look AGN transitioning from a type 2 (no broad-line emission) towards a (nearly) type 1. The broad-line emission has then been confirmed by our 2022 follow-up observations, whereas the rising continuum has no longer been detected, which hints at a further evolution backwards to a (nearly) type 2. The presence of broad-line components also allowed us to obtain the first single-epoch estimate of the black hole mass (log(MBH/Msun) $\sim$ 8.1) in this source. The observed spectral variability might be the result of a change in the accretion activity of NGC 4156, although variable absorption cannot be completely excluded.

We obtained ultraviolet and optical spectra for 9 M~dwarfs across a range of rotation periods to determine whether they showed stochastic intrinsic variability distinguishable from flares. The ultraviolet spectra were observed during the Far Ultraviolet M~Dwarf Evolution Survey \emph{Hubble~Space~Telescope} program using the Space Telescope Imaging Spectrograph. The optical observations were taken from the Apache Point Observatory 3.5-meter telescope using the Dual Imaging Spectrograph and from the Gemini South Observatory using the Gemini Multi-Object Spectrograph. We used the optical spectra to measure multiple chromospheric lines: the Balmer series from H$\alpha$ to H$10$ and the \ion{Ca}{2}~H and K lines. We find that after excising flares, these lines vary on the order of $1-20\%$ at minute-cadence over the course of an hour. The absolute amplitude of variability was greater for the faster rotating M~dwarfs in our sample. Among the 5 stars for which measured the weaker Balmer lines, we note a tentative trend that the fractional amplitude of the variability increases for higher order Balmer lines. We measured the integrated flux of multiple ultraviolet emission features formed in the transition region: the \ion{N}{5}, \ion{Si}{4} and \ion{C}{4} resonance line doublets, and the \ion{C}{2} and \ion{He}{2} multiplets. The signal-to-noise (S/N) ratio of the UV data was too low for us to detect non-flare variability at the same scale and time cadence as the optical. We consider multiple mechanisms for the observed stochastic variability and propose both observational and theoretical avenues of investigation to determine the physical causes of intrinsic variability in the chromospheres of M~dwarfs.

The distribution of molecules between the gas and solid phase during star and planet formation, determines the trajectory of gas and grain surface chemistry, as well as the delivery of elements to nascent planets. This distribution is primarily set by the binding energies of different molecules to water ice surfaces. We computationally estimated the binding energies of ten astrochemically relevant P-bearing species on water surface, we also validate our method for 20 species with known binding energies. We used DFT calculations (M06-2X/aug-cc-pVDZ) to calculate the energetics of molecules and water-molecule clusters (1-3 H$_2$O molecules) and from this determined the binding energy by comparing the complex and the separate molecule and cluster energies. We also explore whether these estimates can be improved by first calibrating our computational method using experimentally measured binding energies. Using the 20 reference molecules we find that the 2H$_2$O cluster size yields the best binding energy estimates and that the application of a calibration to the data may improve the results for some classes of molecules, including more refractory species. Based on these calculations we find that, small P-bearing molecules such as PH$_3$, PN, PO, HPO, PO$_2$ and POOH are relatively volatile and should desorb prior or concomitantly with water ice, while H$_2$PO, HPO$_2$, PO$_3$, PO$_2$OH can strongly bind to any hydroxylated surface, and will likely remain on the interstellar grains surface past the desorbtion of water ice. The depletion of P-carriers on grains constitute a pathway for the inclusion of Phosphorous molecules in planets and planetesimals.

The first critical fast Mach number is defined for a magnetohydrodynamic shock as the Mach number where the shock transitions from subcritical, laminar, behavior to supercritical behavior, characterized by incident ion reflection from the shock front. The ensuing upstream waves and turbulence are convected downstream leading to a turbulent shock structure. Formally this is the Mach number where plasma resistivity can no longer provide sufficient dissipation to establish a stable shock, and is characterized by the downstream flow speed becoming subsonic. We revisit these calculations, including in the MHD jump conditions terms modeling the plasma energy loss to accelerated particles and the presence of waves associated with these particles. The accelerated particle contributions make an insignificant change, but the associated waves have a more important effect. Upstream waves can be strongly amplified in intensity on passing through the shock, and represent another means of shock dissipation. The presence of such waves therefore increases the first critical fast Mach number, especially at quasi-parallel shock where wave excitation is strongest. These effects may have significance for the solar regions where shock waves accelerate particles and cause Type II and Type III radio bursts, and also contribute to the event-to-event variability of SEP acceleration.

Galaxies in the local Universe are thought to require ongoing replenishment of their gas reservoir in order to maintain the observed star formation rates. Cosmological simulations predict that such accretion can occur in both a dynamically hot and cold mode. However, until now observational evidence of the accretion required to match the observed star formation histories is lacking. This paper attempts to determine whether galaxies in the local Universe possess a significant reservoir of HI and what would be the accretion rates derived from such reservoirs. We search the vicinity of 22 nearby galaxies for isolated HI clouds or distinct streams in a systematic and automated manner. The HALOGAS observations represent one of the most sensitive and detailed HI surveys to date. These observations typically reach column density sensitivities of 10^19 cm^-2 over a 20 km/s width. We find 14 secure HI cloud candidates without an observed optical counterpart. These cloud candidates appear to be analogues to the most massive clouds detected around the Milky Way and M31. However, on average their numbers seem significantly reduced. We constrain upper limits for HI accretion in the local Universe. The average HI mass currently observed amounts to a rate of 0.05 Msun/yr with a stringent upper limit of 0.22 Msun/yr, confirming previous estimates. This is much lower than the average star formation rate in this sample. Our best estimate, based on GBT detection limits of several galaxies, suggests that another 0.04 Msun/yr could be accreted from undetected clouds and streams. These results show that in nearby galaxies HI is not being accreted at the same rate as stars are currently being formed. Our study can not exclude that other forms of gas accretion are at work. However, these observations also do not reveal extended neutral gas reservoirs around most nearby spiral galaxies.

SS433 is a Galactic microquasar with powerful outflows (double jet, accretion disk and winds) with well known orbital, precessional and nutational period. In this work we characterise different outflow parameters throughout the precessional cycle of the system. We analyse 10 NuSTAR ($3-70$ keV) observations of $\sim$30~ks that span $\sim$1.5 precessional cycles. We extract averaged spectra and model them using a combination of a double thermal jet model (bjet) and pure neutral and relativistic reflection (xillverCp and relxilllpCp) over an accretion disk. We find an average jet bulk velocity of $\beta = v/c \sim0.29$ with an opening angle of $\lesssim$6~degrees. Eastern jet kinetic power ranges from 1 to $10^{39}$~erg/s, with base "coronal" temperatures $T_o$ ranging between 14 and 18 keV. Nickel to iron abundances remain constant at $\sim$9 (within 1$\sigma$). The western to eastern jet flux ratio becomes $\sim1$ on intermediate phases, about 35% of the total precessional orbit. The $3-70$ keV total unabsorbed luminosity of the jet and disk ranges from 2 to 20 $\times$10$^{37}$~erg/s, with the disk reflection component contributing mainly to the hard $20-30$ keV excess and the stationary 6.7 keV ionized Fe line complex. At low opening angles $\Theta$ we find that the jet expands sideways following an adiabatic expansion of a gas with temperature $T_o$. Finally, the central source and lower parts of the jet could be hidden by an optically thick region of $\tau > 0.1$ and size $R\sim N_H/n_{e0}\sim1.5\times10^9$~cm$\sim$1700~$r_g$ for $M_{BH}=3~M_{\odot}$

We summarize the comet science provided by surveys. This includes surveys where the detections of comets are an advantageous benefit but were not part of the survey\'s original intent, as well as some pointed surveys where comet science was the goal. Many of the surveys are made using astrophysical and heliophysics assets. The surveys in our scope include those using ground-based as well as space-based telescope facilities. Emphasis is placed on current or recent surveys, and science that has resulted since the publication of Comets II, though key advancements made by earlier surveys (e.g. IRAS, COBE, NEAT, etc.) will be mentioned. The proportionally greater number of discoveries of comets by surveys have yielded in turn larger samples of comet populations and sub-populations for study, resulting in better defined evolutionary trends. While providing an array of remarkable discoveries, most of the survey data has been only cursorily investigated. It is clear that continuing to fund ground- and space-based surveys of large numbers of comets is vital if we are to address science goals that can give us a population-wide picture of comet properties.

If high-mass stars (>20-25 Msun) are the progenitors of stripped-envelope (SE) supernovae (SNe), their massive ejecta should lead to broad, long-duration lightcurves (LCs). Instead, literature samples of SE~SNe have reported relatively narrow LCs with ejecta masses between 1-4 Msun that favor progenitors <20-25 Msun. Working with the untargeted sample of (i)PTF SNe to better constrain their rates, we search for SE~SNe with broad LCs. Using a simple LC stretch compared to a template to measure broadness, we identified eight significantly broader Type~Ibc SNe after applying quantitative sample selection criteria. The LCs, colors, and spectra of these SNe are found to evolve more slowly relative to typical Type~Ibc SNe, proportional with the stretch. Bolometric LC modeling and their nebular spectra indicate high ejecta and nickel masses, assuming radioactive decay powering. Additionally, these objects are preferentially located in low-metallicity host galaxies with high star-formation rates, which may account for their massive progenitors, as well as their relative absence from the literature. Our study thus supports the link between broad LCs (as measured by stretch) and high-mass progenitor stars in SE~SNe with independent evidence from bolometric LC modeling, nebular spectra, host environment properties, and photometric evolution. In the first systematic search of its kind using an untargeted sample, we use the stretch distribution to identify a higher than previously appreciated fraction of SE~SNe with broad LCs (~13%). Correcting for Malmquist and LC duration observational biases, we conservatively estimate that a minimum of ~6% of SE~SNe are consistent with high-mass progenitors. This result has implications for the progenitor channels of SE~SNe, late stages of massive stellar evolution, oxygen fraction in the universe, and formation channels for stellar-mass black holes.(Abridged)

The number of terrestrial exoplanets accessible to high-contrast coronagraphic imaging with large telescopes is limited by the smallest angular offset from bright stars at which coronagraphs can observe. However, it is possible to reach inside a telescopes coronagraphic regime by employing nulling interferometry across a telescopes pupil. Indeed, cross-aperture nulling interferometry can observe significantly closer to stars than typical coronagraphs, enabling observations even within the stellar diffraction core. Identifying an optimal nulling coronagraph, i.e., one with both a very small IWA and a high throughput for exoplanet light, would thus be of great interest. A systematic examination of available nulling options has therefore been carried out, which has led to three things. The first is a topological overview that unites both multi-aperture nulling interferometers and single-aperture phase coronagraphs into a common geometrical framework. The second is a new type of phase-mask coronagraph that has emerged from a gap in this framework, called here the split-ring coronagraph. The third is a clear identification of the optimal configuration for a nulling coronagraph, which turned out to be an aperture-plane phase knife, i.e., an achromatic pi-radian phase shift applied to half the telescope pupil prior to focusing the telescopes point spread function into a single-mode fiber. The theoretical peak efficiency of the phase-knife fiber coronagraph, 35.2 percent for a circular telescope aperture, is found to be almost twice that of the next most efficient case, the vortex fiber nuller, at 19.0 percent.

The long-term dynamical evolution is a crucial point in recent planetary research. Although the amount of observational data is continuously growing and the precision allows us to obtain accurate planetary orbits, the canonical stability analysis still requires N-body simulations and phase space trajectory investigations. We propose a method for stability analysis of planetary motion based on the generalized R\'enyi entropy obtained from a scalar measurement. The radial velocity data of the central body in the gravitational three-body problem is used as the basis of a phase space reconstruction procedure. Then, Poincar\'e's recurrence theorem contributes to finding a natural partitioning in the reconstructed phase space to obtain the R\'enyi entropy. It turns out that the entropy-based stability analysis is in good agreement with other chaos detection methods, and it requires only a few tens of thousands of orbital period integration time.

Neutron star (NS) equation of state (EoS) insensitive relations or universal relations (UR) involving neutron star bulk properties play a crucial role in gravitational-wave astronomy. Considering a wide range of equations of state originating from (i) phenomenological relativistic mean field models, (ii) realistic EoS models based on different physical motivations, and (iii) polytropic EoSs described by spectral decomposition method, we update the EoS-insensitive relations involving NS tidal deformability (Multipole Love relation) and the UR between f-mode frequency and tidal deformability (f-Love relation). We analyze the binary neutron star (BNS) event GW170817 using the frequency domain TaylorF2 waveform model with updated universal relations and find that the additional contribution of the octupolar electric tidal parameter and quadrupolar magnetic tidal parameter or the change of multipole Love relation has no significant impact on the inferred NS properties. However, adding the f-mode dynamical phase lowers the 90% upper bound on $\tilde{\Lambda}$ by 16-20% as well as lowers the upper bound of NSs radii by $\sim$500m. The combined URs (multipole Love and f-Love) developed in this work predict a higher median (also a higher 90% upper bound) for $\tilde{\Lambda}$ by 6% and also predict higher radii for the binary components of GW170817 by 200-300m compared to the URs used previously in the literature. We further perform injection and recovery studies on simulated events with different EoSs in $\rm A+$ detector configuration as well as with third generation (3G) Einstein telescope. In agreement with the literature, we find that neglecting f-mode dynamical tides can significantly bias the inferred NS properties, especially for low mass NSs. However, we also find that the impact of the URs is within statistical errors.

LS I +61{\deg}303 is a high mass X-ray binary that is also catalogued as a gamma-ray binary as a result of frequent outbursts at TeV photon energies. The system has released two soft-gamma flares in the past, suggesting a magnetar interpretation for the compact primary. This inference has recently gained significant traction following the discovery of transient radio pulses, detected in some orbital phases from the system, as the measured rotation and tentative spindown rates imply a polar magnetic field strength of $B_p \gtrsim 10^{14}\,\mbox{G}$ if the star is decelerating via magnetic dipole braking. In this paper, we scrutinise magnetic field estimates for the primary in LS I +61{\deg}303 by analysing the compatibility of available data with the system's accretion dynamics, spin evolution, age limits, gamma-ray emissions, and radio pulsar activation. We find that the neutron star's age and spin evolution are theoretically difficult to reconcile unless a strong propeller torque is in operation. This torque could be responsible for the bulk of even the maximum allowed spindown, potentially weakening the inferred magnetic field by more than an order of magnitude.

Fast cosmological transients such as fast radio bursts (FRBs) and gamma-ray bursts (GRBs) represent a class of sources more compact than any other cosmological object. As such they are sensitive to significant magnification via gravitational lensing from a class of lenses which are not well-constrained by observations today. Low-mass primordial black holes are one such candidate which may constitute a significant fraction of the Universe's dark matter. Current observations only constrain their density in the nearby Universe, giving fast transients from cosmological distances the potential to form complementary constraints. Motivated by this, we calculate the effect that gravitational lensing from a cosmological distribution of compact objects would have on the observed rates of FRBs and GRBs. For static lensing geometries, we rule out the prospect that all FRBs are gravitationally lensed for a range of lens masses and show that lens masses greater than $10^{-5}M_\odot$ can be constrained with 8000 un-localised high fluence FRBs at 1.4GHz, as might be detected by the next generation of FRB-finding telescopes.

Dark matter could comprise, at least in part, primordial black holes (PBH). To test this hypothesis, we present an approach to constrain the PBH mass ($M_{\rm{PBH}}$) and mass fraction ($f_{\rm{PBH}}$) from the flux ratios of quadruply imaged quasars. Our approach uses an approximate Bayesian computation (ABC) forward modeling technique to directly sample the posterior distribution of $M_{\rm{PBH}}$ and $f_{\rm{PBH}}$, while marginalizing over the subhalo mass function amplitude, spatial distribution, and the size of the lensed source. We apply our method to 11 quadruply-imaged quasars and derive a new constraint on the intermediate-mass area of PBH parameter space $10^4 $M$_{\odot}<M_{\rm{PBH}}<10^6$M$_\odot$. We obtain an upper limit $f_{\mathrm{PBH}}<0.17$ (95\% C.L.). This constraint is independent of all other previously published limits.

We present a new formulation for numerically obtaining axisymmetric equilibrium structures of rotating stars in two spatial dimensions. With a view to apply it to the secular evolution of rotating stars, we base it on the Lagrangian description, i.e., we solve the force-balance equations to find the spatial positions of fluid elements endowed individually with a mass, specific entropy and angular momentum. The system of nonlinear equations obtained by finite-differencing the basic equations are solved with the W4 method, which is a new multi-dimensional root-finding scheme of our own devising. We augment it with a remapping scheme to avoid distortions of the Lagrangian coordinates. In this first one of a series of papers, we will give a detailed description of these methods initially. We then present the results of some test calculations, which include the construction of both rapidly rotating barotropic and baroclinic equilibrium states. We gauge their accuracies quantitatively with some diagnostic quantities as well as via comparisons with the counterparts obtained with an Eulerian code. For a demonstrative purpose, we apply the code to a toy-model cooling calculation of a rotating white dwarf.

Context. Gamma Doradus (hereafter $\gamma$~Dor) stars are gravity-mode pulsators whose periods carry information about the internal structure of the star. These periods are especially sensitive to the internal rotation and chemical mixing, two processes that are currently not well constrained in the theory of stellar evolution. Aims. We aim to identify the pulsation modes and deduce the internal rotation and buoyancy travel time for 106 $\gamma$ Dor stars observed by the Transiting Exoplanet Survey Satellite (TESS) mission in its southern continuous viewing zone (hereafter S-CVZ). We rely on 140 previously detected period-spacing patterns, that is, series of (near-)consecutive pulsation mode periods. Methods. We used the asymptotic expression to compute gravity-mode frequencies for ranges of the rotation rate and buoyancy travel time that cover the physical range in $\gamma$~Dor stars. Those frequencies were fitted to the observed period-spacing patterns by minimising a custom cost function. The effects of rotation were evaluated using the traditional approximation of rotation, using the stellar pulsation code GYRE. Results. We obtained the pulsation mode identification, internal rotation and buoyancy travel time for 60 TESS $\gamma$~Dor stars. For the remaining 46 targets, the detected patterns are either too short or contained too many missing modes for unambiguous mode identification, and longer light curves are required. For the successfully analysed stars, we found that period-spacing patterns from 1-yr long TESS light curves can constrain the internal rotation and buoyancy travel time to a precision of $\rm 0.03~d^{-1}$ and 400s, respectively, which is about half as precise as literature results based on 4-yr Kepler light curves of $\gamma$~Dor stars.

The cluster mass-richness relation (MRR) is an observationally efficient and potentially powerful cosmological tool for constraining the mean matter density of the universe and the amplitude of fluctuations using the cluster abundance technique. We derive the MRR relation using GalWCat19, a publicly available galaxy cluster catalog we created from the Sloan Digital Sky Survey-DR13 spectroscopic dataset. The MRR shows a tail at the low-richness end. Using the Illustris-TNG and mini-Uchuu cosmological numerical simulations, we demonstrate that this tail is caused by systematical uncertainties. We show that, by means of a judicious cut, identified by the use of the Hinge function, it is possible to determine a richness threshold above which the MRR is linear i.e., where cluster mass scales with richness as logM_200 = alpha + beta logN_200. We derive the MRR and show it is consistent with both sets of simulations with a slope of beta ~ 1. We use our MRR to estimate cluster masses from the GalWCat19 catalog which we then use to set constraints on omega_m and sigma_8. Utilizing the all-member MRR, we obtain constraints of omega_m = 0.31 (+0.04-0.03) and sigma_8 = 0.82 (+0.05-0.04), and utilizing the red-member MRR, we obtain omega_m = 0.31 (+0.04-0.03) and sigma_8 = 0.81 (+0.05-0.04). Our constraints on omega_m and sigma_8 are consistent and very competitive with the Planck 2018 results.

We study the impact of neutral hydrogen absorption on ${\rm H_2}$ photodissociation in protogalactic haloes exposed to soft-UV radiation. Lyman-series absorption can significantly deplete dissociating photons as line overlap with the ${\rm H_2}$ Lyman-Werner bands occurs for neutral column densities exceeding $10^{22}$ ${\rm cm^{-2}}$, but this effect has not been previously included in studies of protogalactic haloes. We use high-resolution three-dimensional hydrodynamic simulations to investigate this "HI-shielding" in three metal-free atomic cooling haloes collapsing at redshift $z \sim 10-20$. We use CLOUDY modeling to update a previous fitting formula for HI-shielding which is a better model for shielding of non-ground state ${\rm H_2}$ rovibrational populations and implement the new fit in our simulations. We find that the inclusion of HI-shielding increases the "critical flux" for suppression of ${\rm H_2}$ cooling in these haloes by $\sim 60-100$ per cent. The larger critical flux has implications in particular for the predicted numbers of candidate haloes in which "direct collapse" could seed massive ($\sim 10^5$ ${\rm M_\odot}$) black holes at $z \sim 15$.

The dividing line between gamma-ray bursts (GRBs) and ordinary stripped-envelope core-collapse supernovae (SNe) is yet to be fully understood. Observationally mapping the variety of ejecta outcomes (ultra-relativistic, mildly-relativistic or non-relativistic) in SNe of Type Ic with broad lines (Ic-BL) can provide a key test to stellar explosion models. However, this requires large samples of the rare Ic-BL events with follow-up observations in the radio, where fast ejecta can be probed largely free of geometry and viewing angle effects. Here, we present the results of a radio (and X-ray) follow-up campaign of 16 SNe Ic-BL detected by the Zwicky Transient Facility (ZTF). Our radio campaign resulted in 4 counterpart detections and 12 deep upper limits. None of the events in our sample is as relativistic as SN 1998bw and we constrain the fraction of SN 1998bw-like explosions to $< 19\%$ (3$\sigma$ Gaussian equivalent), a factor of $\approx 2$ smaller than previously established. We exclude relativistic ejecta with radio luminosity densities in between $\approx 5\times10^{27}$ erg s$^{-1}$ Hz$^{-1}$ and $\approx 10^{29}$ erg s$^{-1}$ Hz$^{-1}$ at $t\gtrsim 20$ d since explosion for $\approx 60\%$ of the events in our sample. This shows that SNe Ic-BL similar to the GRB-associated SN 1998bw, SN 2003lw, SN 2010dh, or to the relativistic SN 2009bb and iPTF17cw, are rare. Our results also exclude an association of the SNe Ic-BL in our sample with largely off-axis GRBs with energies $E\gtrsim 10^{50}$ erg. The parameter space of SN2006aj-like events (faint and fast-peaking radio emission) is, on the other hand, left largely unconstrained and systematically exploring it represents a promising line of future research.

We recently revealed that bulges and elliptical galaxies broadly define distinct, super-linear relations in the $M_{\rm bh}$-$M_{\rm *,sph}$ diagram, with the order-of-magnitude lower $M_{\rm bh}/M_{\rm *,sph}$ ratios in the elliptical galaxies due to major (disc-destroying, elliptical-building) dry mergers. Building on this, here we present a more nuanced picture. Galaxy mergers, in which the net orbital angular momentum does not cancel, can lead to systems with a rotating disc. This situation can occur with either wet (gas-rich) mergers involving one or two spiral galaxies, e.g., NGC 5128, or dry (gas-poor) collisions involving one or two lenticular galaxies, e.g., NGC 5813. The spheroid and disc masses of the progenitor galaxies and merger remnant dictate the shift in the $M_{\rm bh}$-$M_{\rm *,sph}$ and $M_{\rm bh}$-$R_{\rm e,sph}$ diagrams. We show how this explains the (previously excluded merger remnant) S\'ersic S0 galaxies near the bottom of the elliptical sequence and core-S\'ersic S0 galaxies at the top of the bulge sequence, neither of which are faded spiral galaxies. We also introduce two ellicular (ES) galaxy types, explore the location of brightest cluster galaxies and stripped `compact elliptical' galaxies in the $M_{\rm bh}$-$M_{\rm *,sph}$ diagram, and present a new merger-built $M_{\rm bh}$-$M_{\rm *,sph}$ relation which may prove helpful for studies of nanohertz gravitational waves. This work effectively brings into the fold many systems previously considered outliers with either overly massive or undermassive black holes relative to the near-linear $M_{\rm bh}$-$M_{\rm *,sph}$ `red sequence' patched together with select bulges and elliptical galaxies.

More and more experiments have identified that the energy spectra of both primary and secondary cosmic-rays exhibit a hardening above $\sim 200$ GV. Most recently, the DAMPE experiment has reported a hardening of boron-to-carbon ratio at $200$ GV. These signs call for modifications of the conventional cosmic-ray (CR) picture. In this work, we propose that the plethoric secondary cosmic rays, for example, boron, antiprotons, originate from the hadronic interactions of freshly accelerated cosmic rays with the interstellar gas near the sources. We find that secondary-to-primary ratios, for example, boron-to-carbon, boron-to-oxygen and antiproton-to-proton ratios, could be well described. The measurements of electrons and positrons could also be accounted for.

We present a search for continuous gravitational wave signals from an unidentified pulsar potentially powering HESS J1427-608, a spatially-unresolved TeV point source detected by the High Energy Stereoscopic System (H.E.S.S). The search uses a semi-coherent algorithm, which combines the maximum likelihood $\mathcal{F}-$statistic with a hidden Markov Model to efficiently detect and track quasi-monochromatic signals that wander randomly in frequency. It uses data from the second observing run of the Advanced Laser Interferometer Gravitational-Wave Observatory. Multi-wavelength observations of the H.E.S.S. source are combined with the proprieties of the population of TeV-bright pulsar wind nebulae to constrain the search parameters. We find no evidence of gravitational-wave emission from this target. We set upper limits on the characteristic wave strain $h_0^{95\%}$ (for circularly-polarised signals) at $95\%$ confidence level in sample sub-bands and interpolate it to estimate the sensitivity in the full band. We find $h_0^{95\%} = 1.3\times 10^{-25}$ near 185~Hz. The implied constraints on the ellipticity and r-mode amplitude reach $\epsilon\leq 10^{-5}$ and $\alpha \leq 10^{-3}$ at 200~Hz, respectively.

The Epoch of Reionization (EoR) remains a poorly understood cosmic era for the most part. Yet, efforts are still going on to probe and understand this epoch. We present a review of the latest developments in the techniques (especially line-intensity mapping) to study the EoR and try to highlight the contribution of the Indian community in this field. Line-emissions like [H I]$_{\rm 21cm}$, Lyman-$\alpha$, [C II]$_{\text{158}\mu\text{m}}$ and their role as tracers in probing the EoR are discussed. While the [H I]$_{\rm 21cm}$ is an excellent probe of the early IGM, the others are mainly targeted to do an unresolved and large-scale survey of the reionizing sources. Techniques to model these signals include simulations and machine learning approaches, along with the challenge to tackle foregrounds or interlopers. We also discuss synergy opportunities among the various tracers that we mention. Synergy addresses different aspects of the problem, which otherwise is difficult or impossible to tackle. They include statistics like cross-power spectrum, cross-bispectrum, and other techniques such as follow-up studies. We present updates on the relevant experiments; these include the upper limits on the [H I]$_{\rm 21cm}$ power spectrum, along with some highlights on high-redshift galaxy surveys. Finally, we highlight what can be improved further within the community: applying machine learning and simulations based on hydrodynamic and radiative-transfer techniques. Next-generation experiments also need to be conceived to address issues currently beyond our reach.

We observe the almost edge-on (i $\sim$ 90 degrees) galaxy NGC 4302 using ALMA (CO) and VLA (H I) to measure the gas disk thickness for investigating the volumetric star formation law (SFL). The recent star formation rate (SFR) is estimated based on a linear combination of IR 24 micron and H$\alpha$ emissions. The measured scale heights of CO and H I increase significantly with radius. Using the scale heights along with the vertically integrated surface densities, we derive the mid-plane volume densities of the gas ($\rho_{\rm HI}$, $\rho_{\rm H_2}$, and $\rho_{\rm gas}$ = $\rho_{\rm HI}$+$\rho_{\rm H_2}$) and the SFR ($\rho_{\rm SFR}$) and compare the volumetric SFL ($\rho_{\rm SFR} \propto \rho_{\rm gas}^n$) with the vertically integrated SFL ($\Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^N$). We find tight power-law correlations between the SFR and the gas (H I, H$_2$, and the total gas) in both volume and surface densities. The power-law indices of the total gas and H I for the volumetric SFL are noticeably smaller than the indices for the vertically integrated SFL while the H$_2$ indices for both cases are similar to each other. In terms of the star formation efficiency (SFE), we find that the molecular and total gas SFEs are roughly constant, while the atomic SFE is clearly decreasing with radius in both cases.

Dust extinction is one of the most reliable tracers of the gas distribution in the Milky Way. The near-infrared (NIR) Vista Variables in the Via Lactea (VVV) survey enables extinction mapping based on stellar photometry over a large area in the Galactic plane. We devise a novel extinction mapping approach, XPNICER, by bringing together VVV photometric catalogs, stellar parameter data from StarHorse catalogs, and previously published Xpercentile and PNICER extinction mapping techniques. We apply the approach to the VVV survey area, resulting in an extinction map that covers the Galactic disk between 295 and 350 degrees at longitude and -2 to 2 degrees at latitude, and the Galactic bulge between -10 and 5 degrees at latitude. The map has 30 arcseconds spatial resolution and it traces extinctions typically up to about 10-20 mag of visual extinction and maximally up to Av~30 mag. We compare our map to previous dust based maps, concluding that it provides a high-fidelity extinction-based map, especially in its ability to recover both the diffuse dust component of the Galaxy and moderately extincted giant molecular cloud regions. The map is especially useful as independent, extinction-based data on the Galactic dust distribution and applicable for a wide range of studies from individual molecular clouds to the studies of the Galactic stellar populations.

We explore two generic hypotheses for tracing the sources of ultra-high energy cosmic rays (UHECRs) in the Universe: star formation rate density or stellar mass density. For each scenario, we infer a set of constraints for the emission mechanisms in the accelerators, for their energetics and for the abundances of elements at escape from their environments. From these constraints, we generate sky maps above 40~EeV expected from a catalog that comprises 410,761 galaxies out to 350 Mpc and provides a near-infrared flux-limited sample to map both stellar mass and star formation rate over the full sky. Considering a scenario of intermittent sources hosted in every galaxy, we show that the main features observed in arrival directions of UHECRs can in turn constrain the burst rate of the sources provided that magnetic-horizon effects are at play in clusters of galaxies.

Asteroids having perihelion distance $q$ $<$ 1.3 AU are classified as near-Earth objects (NEOs), which are divided into different sub-groups: Vatira-class, Atira-class, Aten-class, Apollo-class, and Amor-class. 2020 $AV_2$, the first Vatira (Orbiting totally inside Venus' orbit) was discovered by the Twilight project of the Zwicky Transient Facility (ZTF) on January 4, 2020. Upon the discovery of 2020 $AV_2$, a couple of orbital studies of the short-term orbital evolution of 2020 $AV_2$ have been performed and published (e.g. de la Fuente Marcos & de la Fuente Marcos 2020; Greenstreet 2020). In this present work, we performed an assessment of the long-term orbital evolution of known near-Earth objects and known Atiras under the Yarkovsky effect by using the \textit{Mercury6} N-body code. We considered not only planetary gravitational perturbation but also the non-gravitational Yarkovsky effect. Our calculation shows that the NEOs have generally two dynamical populations, one short-lived and the other long-lived. From our calculation, the transfer probabilities of Atira-class asteroids to Vatira-class asteroids for the first transition are $\sim$13.1 $\pm$ 0.400, $\sim$13.05 $\pm$ 0.005, and $\sim$ 13.25 $\pm$ 0.450 $\%$ for different values of the Yarkovsky force (i.e. obliquity of 0, 90, and 180 deg.), respectively. It suggests that the radiation force may play some role in the long-term evolution of this asteroid population. Finally, our statistical study implicates that there should be 8.14 $\pm$ 0.133 Atira-class asteroids and 1.05 $\pm$ 0.075 Vatira-asteroids of the S-type taxonomy.

In this paper results from the optical photometric observations of the pre-main-sequence star V1180 Cas are reported. The star is a young variable associated with the dark cloud Lynds 1340, located at a distance of 600 pc from the Sun in the star forming region in Cassiopeia. V1180 Cas shows a large amplitude variability interpreted as a combination of accretion-induced and extinction-driven effects. Our data from VRI CCD photometric observations of the star are collected from September 2011 to February 2022. During our monitoring, we recorded several brightness dips with large amplitudes of up to 5 mag. (I-band). At the same time, increases in brightness over periods of several weeks have also been recorded. In this paper, we compare the photometric data obtained for V1180 Cas with observations of other low-mass pre-main sequence objects.

Measuring the rate at which the universe expands at a given time -- the 'Hubble constant' -- has been a topic of controversy since the first measure of its expansion by Edwin Hubble in the 1920's. As early as the 1970's, Sandage et de Vaucouleurs have been arguing about the adequate methodology for such a measurement. Should astronomers focus only on their best indicators, e.g., the Cepheids, and improve the precision of this measurement based on a unique object to the best possible? Or should they 'spread the risks', i.e., multiply the indicators and methodologies before averaging over their results? Is a robust agreement across several uncertain measures, as is currently argued to defend the existence of a 'Hubble crisis' more telling than a single one percent precision measurement? This controversy, I argue, stems from a misconception of what managing the uncertainties associated with such experimental measurements require. Astrophysical measurements, such as the measure of the Hubble constant, require a methodology that permits both to reduce the known uncertainties and to track the unknown unknowns. Based on the lessons drawn from the so-called Hubble crisis, I sketch a methodological guide for identifying, quantifying and reducing uncertainties in astrophysical measurements, hoping that such a guide can not only help to re-frame the current Hubble tension, but serve as a starting point for future fruitful discussions between astrophysicists, astronomers and philosophers.

The burst mode of accretion in massive star formation is a scenario linking the initial gravitational collapse of parent pre-stellar cores to the properties of their gravitationally unstable discs and of their accretion-driven bursts. In this study, we present a series of high-resolution 3D radiation-hydrodynamics numerical simulations for young massive stars formed out of collapsing 100 Mo molecular cores spinning with several values of the ratio of rotational-to-gravitational energies beta=5%-9%. The models include the indirect gravitational potential caused by disc asymmetries. We find that this modifies the barycenter of the disc, causing significant excursions of the central star position, which we term stellar wobbling. The stellar wobbling slows down and protracts the development of gravitational instability in the disc, reducing the number and magnitude of the accretion-driven bursts undergone by the young massive stars, whose properties are in good agreement with that of the burst monitored from the massive protostar M17 MIR. Including stellar wobbling is therefore important for accurate modeling disc structures. Synthetic ALMA interferometric images in the millimeter waveband show that the outcomes of efficient gravitational instability such as spiral arms and gaseous clumps can be detected for as long as the disc is old enough and has already entered the burst mode of accretion.

We use the EAGLE cosmological simulations to study the evolution of the vertical velocity dispersion of cold gas, $\sigma_{z}$, in central disc galaxies and its connection to stellar feedback, gravitational instabilities, cosmological gas accretion and galaxy mergers. To isolate the impact of feedback, we analyse runs that turn off stellar and (or) AGN feedback in addition to a run that includes both. The evolution of $\sigma_z$ and its dependence on stellar mass and star formation rate in EAGLE are in good agreement with observations. Galaxies hosted by haloes of similar virial mass, $\rm M_{200}$, have similar $\sigma_z$ values even in runs where feedback is absent. The prevalence of local instabilities in discs is uncorrelated with $\sigma_z$ at low redshift and becomes only weakly correlated at high redshifts and in galaxies hosted by massive haloes. $\sigma_z$ correlates most strongly with the specific gas accretion rate onto the disc as well as with the degree of misalignment between the inflowing gas and the disc's rotation axis. These correlations are significant across all redshifts and halo masses, with misaligned accretion being the primary driver of high gas turbulence at redshifts $z \lesssim 1$ and for halo masses $\rm M_{200} \lesssim 10^{11.5} M_{\odot}$. Galaxy mergers increase $\sigma_z$ but, because of their rarity, play only a minor role in its evolution. Our results suggest that the turbulence of cold gas in EAGLE discs results from a complex interplay of different physical processes whose relative importance depends on halo mass and redshift.

We present early science results from Deep Investigation of Neutral Gas Origins (DINGO), an HI survey using the Australian Square Kilometre Array Pathfinder (ASKAP). Using ASKAP sub-arrays available during its commissioning phase, DINGO early science data were taken over $\sim$ 60 deg$^{2}$ of the Galaxy And Mass Assembly (GAMA) 23 h region with 35.5 hr integration time. We make direct detections of six known and one new sources at $z < 0.01$. Using HI spectral stacking, we investigate the HI gas content of galaxies at $0.04 < z< 0.09$ for different galaxy colours. The results show that galaxy morphology based on optical colour is strongly linked to HI gas properties. To examine environmental impacts on the HI gas content of galaxies, three sub-samples are made based on the GAMA group catalogue. The average HI mass of group central galaxies is larger than those of satellite and isolated galaxies, but with a lower HI gas fraction. We derive a variety of HI scaling relations for physical properties of our sample, including stellar mass, stellar mass surface density, $NUV-r$ colour, specific star formation rate, and halo mass. We find that the derived HI scaling relations are comparable to other published results, with consistent trends also observed to $\sim$0.5 dex lower limits in stellar mass and stellar surface density. The cosmic HI densities derived from our data are consistent with other published values at similar redshifts. DINGO early science highlights the power of HI spectral stacking techniques with ASKAP.

Although late-type dwarfs and giant stars are substantially different, their flares are thought to originate in similar physical processes and differ only by a scale factor in the energy levels. We study the validity of this approach. We search for characteristics of flares on active giants, which might be statistically different from those on main-sequence stars. We used nearly 4000 flares of 61 giants and 20 stars of other types that were observed with Kepler in long-cadence mode, which is the only suitable database for this comparative study. For every flare, we derived the duration and energy and gathered stellar parameters. Correlations between the flare characteristics and various stellar parameters were investigated. Strong correlations are found between the flare duration and the surface gravity, luminosity, and radii of the stars. Scaled flare shapes appear to be similar on giants and dwarfs with a 30 min cadence. The logarithmic relation of flare energy and duration is steeper for stars with lower surface gravity. Observed flares are longer and more energetic on giants than on dwarfs on average. The generalized linear scaling for the logarithmic relation of flare energy and duration with a universal theoretical slope of $\approx$1/3 should slightly be modified by introducing a dependence on surface gravity.

The atmosphere is one of the most important contamination sources in the ground-based Cosmic Microwave Background (CMB) observations. In this paper, we study three kinds of filters, which are polynomial filter, high-pass filter, and Wiener filter, to investigate their ability for removing atmospheric noise, as well as their impact on the data analysis process through the end-to-end simulations of CMB experiment. We track their performance by analyzing the response of different components of the data, including both signals and noise. In the time domain, the calculation shows that the high-pass filter has the smallest root mean square error and can achieve high filtering efficiency, followed by the Wiener filter and polynomial filter. We then perform map-making with the filtered time ordered data (TOD) to trace the effects from filters on the map domain, and the results show that the polynomial filter gives high noise residual at low frequency, which gives rise to serious leakage to small scales in map domain during the map-making process, while the high-pass filter and Wiener filter do not have such significant leakage. Then we estimate the angular power spectra of residual noise, as well as those of the input signal for comparing the filter effects in the power spectra domain. Finally, we estimate the standard deviation of the filter corrected power spectra to compare the effects from different filters, and the results show that, at low noise level, the three filters give almost comparable standard deviations on the medium and small scales, but at high noise level, the standard deviation of the polynomial filter is significantly larger. These studies can be used for the reduction of atmospheric noise in future ground-based CMB data processing.

We present radial velocity follow-up obtained with ESPRESSO of the M-type star LTT 1445A (TOI-455), for which a transiting planet b with an orbital period of~5.4 days was detected by TESS. We report the discovery of a second transiting planet (LTT 1445A c) and a third non-transiting candidate planet (LTT 1445A d) with orbital periods of 3.12 and 24.30 days, respectively. The host star is the main component of a triple M-dwarf system at a distance of 6.9 pc. We used 84 ESPRESSO high-resolution spectra to determine accurate masses of 2.3$\pm$0.3 $\mathrm{M}_\oplus$ and 1.0$\pm$0.2 $\mathrm{M}_\oplus$ for planets b and c and a minimum mass of 2.7$\pm$0.7 $\mathrm{M}_\oplus$ for planet d. Based on its radius of 1.43$\pm0.09$ $\mathrm{R}_\oplus$ as derived from the TESS observations, LTT 1445A b has a lower density than the Earth and may therefore hold a sizeable atmosphere, which makes it a prime target for the James Webb Space Telescope. We used a Bayesian inference approach with the nested sampling algorithm and a set of models to test the robustness of the retrieved physical values of the system. There is a probability of 85$\%$ that the transit of planet c is grazing, which results in a retrieved radius with large uncertainties at 1.60$^{+0.67}_{-0.34}$ $\mathrm{R}_\oplus$. LTT 1445A d orbits the inner boundary of the habitable zone of its host star and could be a prime target for the James Webb Space Telescope.

Throughout the passage within the Galactic halo, high-velocity clouds (HVCs) sweep up ambient magnetic fields and form stretched and draped configurations of magnetic fields around them. This magnetized layer at the cloud-halo interface can suppress the mixing between the cloud and the surrounding halo gas. For simplicity, many earlier numerical studies adopt spherically symmetric uniform-density clouds as initial conditions. However, observations of HVCs indicate that the clouds are clumpy and turbulent. In this paper, we perform 3D magnetohydrodynamic simulations to study the evolution of clouds with more a realistic initial density distribution. In the parameter regime of our models, simple spherical clouds only grow from the onset of the simulations, while clumpy clouds lose cold gas at early times and then start to grow unless they have very low metallicity. Some massive clumps composing a clumpy cloud share broad similarities with single uniform-density clouds in that they have clear head-tail morphology and magnetic fields draped around themselves. Such similarities are marginal among clumps in the wake of the cloud. With magnetic fields present, the growth of hydrodynamic instabilities is suppressed along the direction of the cloud's motion. Efficient radiative cooling leads to a compact cloud with rich clumps and filaments condensed from the surrounding gas. The slope of the initial density power spectrum is closely related to the density contrast within a clumpy cloud. It, therefore, determines how fast the material is pushed out of the initial cloud and mixed with the halo gas.

Standard sirens -- gravitational wave (GW) sources with an electromagnetic (EM) counterpart -- can be used to measure the Hubble constant directly which should help to ease the existing Hubble tension. However, if the source is moving, a relativistic redshift affects the redshift of the EM counterpart and the apparent distance of the GW source, and thus it needs to be corrected to obtain accurate measurements. We study the effect of velocity on GWs for a source in an expanding universe showing that the total redshift of the wave is equal to the product of the relativistic redshift and the cosmological redshift. We, further, find that a motion of the source changes its apparent distance by a factor $(1+z_{\rm rel})^2$ in contrast to a linear factor for the cosmological redshift. We discuss that the additional factor for the relativistic redshift is a consequence of a velocity-dependent amplitude for GWs. We consider the effect of the velocity on the chirp mass and the apparent distance of the source an observer would infer when ignoring the velocity. We find that for different astrophysical scenarios the error in the chirp mass can range between 0.1\,\% and 7\,\% while the error in the apparent distance can be between 0.25\,\% and 15\,\%. Furthermore, we consider the error introduced in the measurement of the Hubble constant using standard sirens for two cases: (i) when the effect of velocity on the redshift of the EM counterpart is considered but not on the apparent distance obtained from GWs and (ii) when the effect of the velocity is ignored completely. We find that in the first case the error can reach 1\,\% for a source moving due to the peculiar velocity of its host galaxy and that in the second case the error can be more than 5\,\% for a source at the distance of GW150914 with the same velocity.

The race for the most distant object in the Universe has been played by long-duration gamma-ray bursts (GRBs), star-forming galaxies and quasars. GRBs took a temporary lead with the discovery of GRB 090423 at a redshift z=8.2, but now the record-holder is the galaxy GN-z11 at z=11.0. Despite this record, galaxies and quasars are very faint (GN-z11 has a magnitude H=26), hampering the study of the physical properties of the primordial Universe. On the other hand, GRB afterglows are brighter by a factor of >100, with the drawback of lasting only for 1-2 days. Here we describe a novel approach to the discovery of high-redshift (z>6) GRBs, exploiting their near-infrared (nIR) emission properties. Soon after the bright, high-energy prompt phase, a GRB is accompanied by an afterglow. The afterglows of high-redshift GRBs are naturally absorbed, like any other source, at optical wavelengths by Hydrogen along the line of sight in the intergalactic medium (Lyman-alpha absorption). We propose to take advantage of the deep monitoring of the sky by the Vera Rubin Observatory, to simultaneously observe exactly the same fields with a new, dedicated nIR facility. By comparing the two streams of transients, one can pinpoint transients detected in the nIR band and not in the optical band. These fast transients detected only in the nIR and with an AB colour index r-H>3.5 are high-redshift GRBs, with a low contamination rate. Thanks to the depth reached by the Rubin observations, interlopers can be identified, allowing us to discover ~11 GRBs at z>6 per year and ~3 GRBs per year at z>10. This turns out to be one of the most effective probes of the high-redshift Universe.

We present ALMA band 6 observations of the luminous blue variable Eta Car, obtained within the ALMAGAL program. We report SiO J=5-4, SiS J=12-11 and SiN N=5-4 emission in the equatorial region of the Homunculus nebula, constituting the first detection of silicon- and sulphur-bearing molecules in the outskirts of a highly evolved, early-type massive star. SiO, SiS and SiN trace a clumpy equatorial ring that surrounds the central binary at a projected distance of 2 arcsec, delineating the inner rims of the butterfly-shaped dusty region. The formation of silicon-bearing compounds is presumably related to the continuous recycling of dust due to the variable wind regime of Eta Car, that destroys grains and releases silicon back to gas phase. We discuss possible formation routes for the observed species, contextualizing them within the current molecular inventory of Eta Car. We find that the SiO and SiS fractional abundances in localised clumps of the ring, $6.7\times10^{-9}$ and $1.2\times10^{-8}$ respectively, are exceptionally lower than those measured in C- and O-rich AGB stars and cool supergiants; while the higher SiN abundance, $3.6\times10^{-8}$, evidences the nitrogen-rich chemistry of the ejecta. These abundances must be regarded as strict upper limits, since the distribution of H2 in the Homunculus is unknown. In any case, these findings shed new light onto the peculiar molecular ecosystem of Eta Car, and establish its surroundings as a new laboratory to investigate the lifecycle of silicate dust in extreme astrophysical conditions.

Studies of solar twins have key impacts on the astronomical community, but only $\sim$100--200 nearby solar twins ($<$ 1 kpc) have been reliably identified over the last few decades. The aim of our survey (SDST) is to identify $\sim$150--200 distant solar twins and analogues (up to $\lesssim$ 4 kpc) closer to the Galactic Centre. We took advantage of the precise Gaia and Skymapper surveys to select Sun-like candidates in a 2-degree field, which were observed with the HERMES spectrograph on the Anglo-Australian Telescope. We successfully built up the required signal-to-noise ratio (25-per-pixel in the HERMES red band) for most targets as faint as Gaia G of 17.4 mag. The stellar photometric/astrometric parameters (e.g., \teff, \logg, mass) of our candidates are derived in this paper, while the spectroscopic parameters will be presented in the third paper in this SDST series. The selection success rate - the fraction of targets which belong to solar twins or analogues - was estimated from simulated survey data and the Besan\c{c}on stellar population model, and compared with the actual success rate of the survey. We find that expected and actual success rates agree well, indicating that the numbers of solar twins and analogues we discover in SDST are consistent with expectations, affirming the survey approach. These distant solar analogues are prime targets for testing for any variation in the strength of electromagnetism in regions of higher dark matter density, and can make additional contributions to our understanding of, e.g., Galactic chemical evolution in the inner Milky Way.

Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $\gamma$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process regression. We find that the jet particle and magnetic field energy densities are in equipartition ($k_\textrm{r} = 1.08\pm0.18$). The jet apex is located $z_\textrm{91.5 GHz}= 22 - 645$ $R_\textrm{s}$ ($2 - 20 \times 10^{-3}$ pc) upstream of the 3 mm radio core; at that position, the magnetic field amplitude is $B_\textrm{core}^\textrm{91.5 GHz}= 3 - 10$ G. Our results are in good agreement with earlier studies, which utilised very-long-baseline interferometry. Furthermore, we investigate the temporal relation between the ejection of radio and $\gamma$-ray flares. Our results are in favour of the $\gamma$-ray emission being associated with the radio emission. We are able to tentatively connect the ejection of features identified at 43 and 86 GHz to prominent $\gamma$-ray flares. Finally, we compute the multiplicity parameter $\lambda$ and the Michel magnetisation $\sigma_\textrm{M}$ and find that they are consistent with a jet launched by the Blandford & Znajek 1977 mechanism.

In the study of radiation mechanisms of AGNs' jets, hadronuclear (pp) interactions are commonly neglected, because the number density of cold protons in the jet is considered insufficient. Very recently, our previous work proves that pp interactions in the low TeV luminosity AGNs, which have potential to generate detectable very-high-energy (VHE) emission, could be important. Based on this, the one-zone pp model is employed to study low-state quasi-simultaneous spectral energy distributions of a sample of low TeV luminosity AGNs in this work. Our modeling results show that the gamma-ray generated in pp interactions can explain the observed TeV spectra and has contribution to higher energy band that could be detected by the Large High Altitude Air Shower Observatory (LHAASO). In the sample, we suggest that M 87, Mrk 421 and Mrk 501 are the most likely objects to be detected by LHAASO in the near future. Other possible origins of VHE emission are also briefly discussed.

We continue our exploration of the wiggly generalisation of the Velocity-Dependent One Scale Model for cosmic strings, through the study of its allowed asymptotic scaling solutions. We extend the work of a previous paper [Almeida $\&$ Martins, Phys. Rev. D 104 (2021) 043524] by considering the more comprehensive case of a time-varying coarse-graining scale for the string wiggles. The modeling of the evolution of the network therefore relies on three main mechanisms: Hubble expansion, energy transfer mechanisms (e.g., the production of loops and wiggles) and the choice of the scale at which wiggles are coarse-grained. We analyse the role of each of them on the overall behaviour of the network, and thus in the allowed scaling solutions. In Minkowski space, we find that linear scaling, previously observed in numerical simulations without expansion, is not possible with a changing averaging scale. For expanding universes, we find that the three broad classes of scaling solutions -- with the wiggliness disappearing, reaching scaling, or growing -- still exist but are differently impacted by the time evolution of the coarse-graining scale. Nambu-Goto type solutions (without wiggles) are unaffected, growing wiggliness solutions are trivially generalized, while for solutions where wiggliness reaches scaling the expansion rate for which the solution exists is decreased with respect to the one for a fixed coarse-graining scale. Finally, we also show that the inclusion of a time-varying coarse-graining scale allows, in principle, for additional scaling solutions which, although mathematically valid, are not physical. Overall, our mapping of the landscape of the allowed scaling solutions of the wiggly Velocity-Dependent One Scale Model paves the way for the detailed testing of the model, to be done by forthcoming high-resolution field theory and Nambu-Goto simulations.

The extreme space weather conditions resulting from high energetic events likes solar flares and Coronal Mass Ejections (CMEs) demand for reliable space weather forecasting. The magnetic flux tubes while rising through the convection zone gets twisted by the turbulent plasma flows, energizing the system and resulting in flares. We investigate the relationship between the subsurface plasma flows associated with flaring active regions and their surface magnetic flux and current helicity. The near-surface horizontal velocities derived from the ring-diagram analysis of active region patches using Global Oscillation Network Group (GONG) Doppler velocity measurements are used to compute the fluid dynamics descriptors like vertical divergence, vorticity and kinetic helicity used in this work. The flaring active regions are observed to have large value of vertical vorticity and kinetic helicity. Also, the horizontal flow divergence, vorticity, flux, kinetic and current helicities are observed to be significantly correlated and evolve in phase with each other. We observe that the integrated values of the above flow and magnetic parameters observed one day prior to the flare are significantly correlated with the integrated flare intensity of the active region. Hence, we show that strong vorticity/kinetic helicities lead to larger active region twisting, presumably generating high-intensity flares.

Turbulent plasma motion is common in the universe, and invoked in solar flares to drive effective acceleration leading to high energy electrons. Unresolved mass motions are frequently detected in flares from extreme ultraviolet (EUV) observations, which are often regarded as turbulence. However, how this plasma turbulence forms during the flare is still largely a mystery. Here we successfully reproduce observed turbulence in our 3D magnetohydrodynamic simulation where the magnetic reconnection process is included. The turbulence forms as a result of an intricate non-linear interaction between the reconnection outflows and the magnetic arcades below the reconnection site, in which the shear-flow driven Kelvin-Helmholtz Instability (KHI) plays a key role for generating turbulent vortices. The turbulence is produced above high density flare loops, and then propagates to chromospheric footpoints along the magnetic field as Alfvenic perturbations. High turbulent velocities above 200 km s^-1 can be found around the termination shock, while the low atmosphere reaches turbulent velocities of 10 km s^-1 at a layer where the number density is about 10^11 cm^-3. The turbulent region with maximum non-thermal velocity coincides with the region where the observed high-energy electrons are concentrated, demonstrating the potential role of turbulence in acceleration. Synthetic views in EUV and fitted Hinode-EIS spectra show excellent agreement with observational results. An energy analysis demonstrates that more than 10% of the reconnection downflow kinetic energy can be converted to turbulent energy via KHI.

In this paper we consider the effects of adding curvature in extended cosmologies involving a free-to-vary neutrino sector and different parametrizations of Dark Energy (DE). We make use of the Planck 2018 cosmic microwave background temperature and polarization data, Baryon Acoustic Oscillations and Pantheon type Ia Supernovae data. Our main result is that a non-flat Universe cannot be discarded in light of the current astronomical data, because we find an indication for a closed Universe in most of the DE cosmologies explored in this work. On the other hand, forcing the Universe to be flat can significantly bias the constraints on the equation of state of the DE component and its dynamical nature.

The unprecedented angular resolution and sensitivity of ALMA makes it possible to unveil disk populations in distant ($>$2 kpc), embedded young cluster environments. We have conducted an observation towards the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15 $\rm \mu Jy~beam^{-1}$, we detect a total of 38 protostellar disks. These disks have dust masses ranging from about 53 to 1825 $M_\oplus$, assuming a dust temperature of 20 K. This sample is not closely associated with previously identified dense cores, as would be expected for disks around Class 0 protostars. Thus, we expect our sample, being flux limited, to be mainly composed of Class I/flat-spectrum source disks, since these are typically more massive than Class II disks. Furthermore, we find that the distributions of disk masses and radii are statistically indistinguishable with those of the Class I/flat-spectrum objects in the Orion molecular cloud, indicating similar processes are operating in G286.21+0.16 to regulate disk formation and evolution. The cluster center appears to host a massive protostellar system composed of three sources within 1200 au, including a potential binary with 600 au projected separation. Relative to this center, there is no evidence for widespread mass segregation in the disk population. We do find a tentative trend of increasing disk radius versus distance from the cluster center, which may point to the influence of dynamical interactions being stronger in the central regions.

Recent studies have revealed a strong relation between sample-averaged black-hole (BH) accretion rate (BHAR) and star formation rate (SFR) among bulge-dominated galaxies, i.e., "lockstep" BH-bulge growth, in the distant universe. This relation might be closely related to the BH-bulge mass correlation observed in the local universe. To understand further BH-bulge coevolution, we present ALMA CO(2-1) or CO(3-2) observations of 7 star-forming bulge-dominated galaxies at z=0.5-2.5. Using the ALMA data, we detect significant ($>3\sigma$) CO emission from 4 objects. For our sample of 7 galaxies, we measure (or constrain with upper limits) their CO line fluxes and estimate molecular gas masses ($M_{gas}$). We also estimate their stellar masses ($M_{star}$) and SFRs by modelling their spectral energy distributions (SEDs). Using these physical properties, we derive the gas-depletion timescales ($t_{dep} = M_{gas}/SFR$) and compare them with the bulge/BH growth timescales ($t_{grow} = M_{star}/SFR \sim M_{BH}/BHAR$). Our sample generally has $t_{dep}$ shorter than $t_{grow}$ by a median factor of $\gtrsim 4$, indicating that the cold gas will be depleted before significant bulge/BH growth takes place. This result suggests that the BH-bulge lockstep growth is mainly responsible for maintaining their mass relation, not creating it. We note that our sample is small and limited to $z<2.5$; JWST and ALMA will be able to probe to higher redshifts in the near future.

A significant number of stellar binary black hole (sBBH) mergers may be lensed and detected by the third generation gravitational wave (GW) detectors. Their lensed host galaxies may be detectable, which thus helps to accurately localize these sources and provide a new approach to study the origin of sBBHs. In this paper, we investigate the detectability of the lensed host galaxies for the lensed sBBH mergers. We find that the detection fraction of the host galaxies to the lensed GW events can be significantly different for a survey with a given limiting magnitude if sBBHs are produced by different mechanisms, such as the evolution of massive binary stars, the dynamical interactions in dense star clusters, and that assisted by active galactic nuclei or massive black holes. Furthermore, we illustrate that the statistical spatial distribution of those lensed sBBHs in its hosts resulting from different sBBH formation channels can be different from each other. Therefore, with the third generation GW detectors and future large scale galaxy surveys, it is possible to independently constrain the sBBH origin via the detection fraction of those lensed events with identifiable lensing host signatures and/or even constrain the contribution fractions from different sBBH formation mechanisms

We present measurements of spot properties on 31 young stellar objects, based on multi-band data from the HOYS (Hunting Outbursting Young Stars) project. On average the analysis for each object is based on 270 data points during 80 days in at least 3 bands. All the young low-mass stars in our sample show periodic photometric variations. We determine spot temperatures and coverage by comparing the measured photometric amplitudes in optical bands with simulated amplitudes based on atmosphere models, including a complete error propagation. 21 objects in our sample feature cool spots, with spot temperatures 500 - 2500 K below the stellar effective temperature ($T_{\rm eff}$), and a coverage of 0.05 - 0.4. Six more have hot spots, with temperatures up to 3000 K above $T_{\rm eff}$ and coverage below 0.15. The remaining four stars have ambiguous solutions or are AA Tau-type contaminants. All of the stars with large spots (i.e. high coverage $>0.1$) are relatively cool with $T_{\rm eff} < 4500$ K, which could be a result of having deeper convection zones. Apart from that, spot properties show no significant trends with rotation period, infrared excess, or stellar properties. Most notably, we find hot spots in stars that do not show $K-W2$ infrared excess, indicating the possibility of accretion across an inner disk cavity or the presence of plage.

Pulsars are rotating neutron stars that emit periodic electromagnetic radiation. While pulsars generally slow down as they lose energy, some also experience glitches: spontaneous increases of their rotational frequency. According to several models, these glitches can also lead to the emission of long-duration transient gravitational waves (GWs). We present detection prospects for such signals by comparing indirect energy upper limits on GW strain for known glitches with the sensitivity of current and future ground-based GW detectors. We first consider generic constraints based on the glitch size and find that realistic matched-filter searches in the fourth LIGO-Virgo-KAGRA observing run (O4) could make a detection, or set constraints below these indirect upper limits, for equivalents of 36 out of 726 previously observed glitches, and 74 in the O5 run. With the third-generation Einstein Telescope or Cosmic Explorer, 35-40% of glitches would be accessible. When specialising to a scenario where transient mountains produce the post-glitch GW emission, following Yim & Jones (2020), the indirect upper limits are stricter. Out of the smaller set of 119 glitches with measured healing parameter, as needed for predictions under that model, only 6 glitches would have been within reach for O4 and 14 for O5, with a similar percentage as before with third generation detectors. We also discuss how this model matches the observed glitch population.

Quasi-periodic oscillation (QPO) signals are discovered in some fast radio bursts (FRBs) such as FRB 20191221A, as well as in the X-ray burst associated with the galactic FRB from SGR 1935+2154. We revisit the intermediate-field FRB model where the radio waves are generated as fast-magnetosonic waves through magnetic reconnection near the light cylinder. The current sheet in the magnetar wind is compressed by a low frequency pulse emitted from the inner magnetosphere to trigger magnetic reconnection. By incorporating the wave dynamics of the magnetosphere, we demonstrate how the FRB frequency, the single pulse width, and luminosity are determined by the period, magnetic field, QPO frequency and quake energetics of the magnetar. We find that this model can naturally and self-consistently interpret the X-ray/radio event from SGR 1935+2154 and the QPO in FRB 20191221A. It can also explain the observed wide energy range of repeating FRBs in a narrow bandwidth.

During the inspiralling of a white dwarf (WD) into an intermediate-mass black hole (~10^{2-5} M_sun), both gravitational waves (GWs) and electromagnetic (EM) radiation are emitted. Once the eccentric orbit's pericenter radius approaches the tidal radius, the WD would be tidally stripped upon each pericenter passage. The accretion of these stripped mass would produce EM radiation. It is suspected that the recently discovered new types of transients, namely the quasi-periodic eruptions and the fast ultraluminous X-ray bursts, might originate from such systems. Modeling these flares requires a prediction of the amount of stripped mass from the WD and the details of the mass supply to the accretion disk. We run hydrodynamical simulations to study the orbital parameter dependence of the stripped mass. We find that our results match the analytical estimate that the stripped mass is proportional to a 5/2-th power of the excess depth by which the WD overfills its instantaneous Roche lobe at the pericenter. The corresponding fallback rate of the stripped mass is calculated, which may be useful in interpreting the individual flaring light curve in candidate EM sources. We further calculate the long-term mass-loss evolution of a WD during its inspiral and the detectability of the GW and EM signals. The EM signal from the mass-loss stage can be easily detected: the limiting distance is ~320(M_h/10^4 M_sun) Mpc for Einstein Probe. The GW signal, for the space-borne detectors such as Laser Interferometer Space Antenna or TianQin, can be detected only within the Local Supercluster (~33 Mpc).

Parity violation is expected to generate an asymmetry between the amplitude of left and right-handed gravitational-wave modes which leads to a circularly polarized stochastic gravitational-wave background (SGWB). Due to the three independent baselines in the LIGO-Virgo network, we focus on the amplitude difference in strain power characterized by Stokes' parameters and do maximum-likelihood estimation to constrain the polarization degree of SGWB. Our results indicate that there is no evidence for the circularly polarized SGWB in the data. Furthermore, by modeling the SGWB as a power-law spectrum, we place upper limit on the normalized energy density $\Omega_\text{gw}(25\,\text{Hz})<5.3\times10^{-9}$ at $95\%$ confidence level after marginalizing over the polarization degree and spectral index.

We describe pynucastro 2.0, an open source library for interactively creating and exploring astrophysical nuclear reaction networks. We demonstrate new methods for approximating rates and using detailed balance to create reverse rates, show how to build networks and determine whether they are appropriate for a particular science application, and discuss the changes made to the library over the past few years. Finally, we demonstrate the validity of the networks produced and share how we use pynucastro networks in simulation codes.

Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to investigate the sensitivity of the central entropy and the shape of the profiles to changes in the sub-grid model applied to a suite of zoom-in cosmological simulations of a group of mass $M_{500} = 8.8 \times 10^{12}~{\rm M}_\odot$ and a cluster of mass $2.9 \times 10^{14}~{\rm M}_\odot$. Using our reference model, calibrated to match the stellar mass function of field galaxies, we confirm that our simulated groups and clusters contain hot gas with too high entropy in their cores. Additional simulations run without artificial conduction, metal cooling or AGN feedback produce lower entropy levels but still fail to reproduce observed profiles. Conversely, the two objects run without supernova feedback show a significant entropy increase which can be attributed to excessive cooling and star formation. Varying the AGN heating temperature does not greatly affect the profile shape, but only the overall normalisation. Finally, we compared runs with four AGN heating schemes and obtained similar profiles, with the exception of bipolar AGN heating, which produces a higher and more uniform entropy distribution. Our study leaves open the question of whether the entropy core problem in simulations, and particularly the lack of power-law cool-core profiles, arise from incorrect physical assumptions, missing physical processes, or insufficient numerical resolution.

We report the discovery of one possible neutron star binary ($P_{\rm orb} =$ 0.8666 day) by using the LAMOST low-resolution spectroscopic data. The visible companion is a late A-type dwarf ($T_{\rm eff} = 7900 \pm 200$ K; log$g$ $=$ 4.3$\pm$0.2; $M =$ 1.7$\pm$0.1 M$_{\odot}$; $R\ =\ 1.7\pm0.2$ R$_{\odot}$), at a distance of 1.11$\pm0.03$ kpc. No double-lined feature can be seen from the GTC/HORuS high-resolution spectra, thus the radial velocity variation indicates an invisible object hiding in the binary. The system's optical light curves show clear ellipsoidal variability, suggesting that the visible companion is tidal distorted. By fitting the multi-band light curves with the ELC and WD codes, we constrain the mass of the invisible star to be 1.1--1.3 M$_{\odot}$. Spectral disentangling shows no additional component with optical absorption spectra, supporting the system contains one compact object. No X-ray or UV emission are detected in the ROSAT archive observations. Therefore, we suspect the invisible object is more likely a neutron star rather than a white dwarf. Our finding suggests the ability of LAMOST spectroscopic survey to discover X-ray quiescent compact objects.

Geometry and dynamical structure of emission regions in accreting pulsars are shaped by the interplay between gravity, radiation, and strong magnetic field, which significantly affects the opacities of a plasma and radiative pressure under such extreme conditions. Quantitative consideration of magnetic plasma opacities is, therefore, an essential ingredient of any self-consistent modeling of emission region structure of X-ray pulsars. We present results of computations of the Rosseland and Planck mean opacities of a strongly magnetized plasma with a simple chemical composition,namely the solar hydrogen/helium mix. We consider all relevant specific opacities of the magnetized plasma including vacuum polarization effect and contribution of electron-positron pairs where the pair number density is computed in the thermodynamic equilibrium approximation. The magnetic Planck mean opacity determines the radiative cooling of an optically thin strongly magnetized plasma. It is by factor of three smaller than non-magnetic Planck opacity at $k_{\rm B}T < 0.1\,E_{\rm cyc}$ and increases by a factor of $10^2 - 10^4$ at $k_{\rm B}T > 0.3\,E_{\rm cyc}$ due to cyclotron thermal processes. We propose a simple approximate expression which has sufficient accuracy for the magnetic Planck opacity description. We provide the Rosseland opacity in a tabular form computed in the temperature range 1 - 300 keV, magnetic field range $3 \times 10^{10} - 10^{15}$ G, and a broad range of plasma densities. We demonstrate that the scattering on the electron-positron pairs increases the Rosseland opacity drastically at temperatures >50 keV in the case of mass densities typical for accretion channel in X-ray pulsars.

Clouds are ubiquitous\, -- \,they arise for every solar system planet that possesses an atmosphere and have also been suggested as a leading mechanism for obscuring spectral features in exoplanet observations. As exoplanet observations continue to improve, there is a need for efficient and general planetary climate models that appropriately handle the possible cloudy atmospheric environments that arise on these worlds. We generate a new 1D radiative-convective terrestrial planet climate model that self-consistently handles patchy clouds through a parameterized microphysical treatment of condensation and sedimentation processes. Our model is general enough to recreate Earth's atmospheric radiative environment without over-parameterization, while also maintaining a simple implementation that is applicable to a wide range of atmospheric compositions and physical planetary properties. We first validate this new 1D patchy cloud radiative-convective climate model by comparing it to Earth thermal structure data and to existing climate and radiative transfer tools. We produce partially-clouded Earth-like climates with cloud structures that are representative of deep tropospheric convection and are adequate 1D representations of clouds within rocky planet atmospheres. After validation against Earth, we then use our partially clouded climate model and explore the potential climates of super-Earth exoplanets with secondary nitrogen-dominated atmospheres which we assume are abiotic. We also couple the partially clouded climate model to a full-physics, line-by-line radiative transfer model and generate high-resolution spectra of simulated climates. These self-consistent climate-to-spectral models bridge the gap between climate modeling efforts and observational studies of rocky worlds.

With the increasing number of planets discovered by TESS, the atmospheric characterization of small exoplanets is accelerating. L98-59 is a M-dwarf hosting a multi-planet system, and so far, four small planets have been confirmed. The innermost planet b is $\sim15\%$ smaller and $\sim60\%$ lighter than Earth, and should thus have a predominantly rocky composition. The Hubble Space Telescope observed five primary transits of L98-59b in $1.1-1.7\ \mu$m, and here we report the data analysis and the resulting transmission spectrum of the planet. We measure the transit depths for each of the five transits and, by combination, we obtain a transmission spectrum with an overall precision of $\sim20$ ppm in for each of the 18 spectrophotometric channels. With this level of precision, the transmission spectrum does not show significant modulation, and is thus consistent with a planet without any atmosphere or a planet having an atmosphere and high-altitude clouds or haze. The scenarios involving an aerosol-free, H$_2$-dominated atmosphere with H$_2$O or CH$_4$ are inconsistent with the data. The transmission spectrum also disfavors, but does not rules out, an H$_2$O-dominated atmosphere without clouds. A spectral retrieval process suggests that an H$_2$-dominated atmosphere with HCN and clouds or haze may be the preferred solution, but this indication is non-conclusive. Future James Webb Space Telescope observations may find out the nature of the planet among the remaining viable scenarios.

Gamma-ray flares of blazars may be accompanied by high-energy neutrinos due to interactions of high-energy cosmic rays in the jet with photons, as suggested by the detection of the high-energy neutrino IceCube-170922A during a major gamma-ray flare from blazar TXS 0506+056 at the $\sim3\sigma$ significance level. In this work, we present a statistical study of gamma-ray emission from blazars to constrain the contribution of gamma-ray flares to their neutrino output. We construct weekly binned light curves for 145 gamma-ray bright blazars in the {\it Fermi} Large Area Telescope (LAT) Monitored Source List. We derive the fraction of time spent in the flaring state (flare duty cycle) and the fraction of energy released during each flare from the light curves with a Bayesian blocks algorithm. We find that blazars with lower flare duty cycles and energy fractions are more numerous among our sample. We identify no significant differences between blazar sub-classes in terms of flaring activity. Then using a general scaling relation for the neutrino and gamma-ray luminosities, $L_{\nu} \propto (L_{\gamma})^{\gamma}$ with a weighting exponent of ${\gamma} = 1.0 - 2.0$, normalized to the quiescent gamma-ray or X-ray flux of each blazar, we evaluate the neutrino energy flux of each gamma-ray flare. The gamma-ray flare distribution indicates that blazar neutrino emission may be dominated by flares for $\gamma\gtrsim1.5$. The neutrino energy fluxes for one-week and 10-year bins are compared with the declination-dependent IceCube sensitivity to constrain the standard neutrino emission models for gamma-ray flares. Finally, we present the upper-limit contribution of blazar gamma-ray flares to the isotropic diffuse neutrino flux.

High-redshift observations are often biased towards massive and bright galaxies that are not necessarily representative of the full population. In order to accurately study galaxy evolution and mass assembly at these redshifts, observations of ``normal'' main sequence galaxies are required. Here we present Atacama Large Millimeter/Submillimeter Array (ALMA) 0.3" resolution observations of the [CII] emission line at 158$\mu$m of HZ7, a main sequence galaxy at $z=5.25$. Comparing to archival rest-frame UV observations taken by the Hubble Space Telescope (HST), we find strong evidence of the existence of extended [CII] emission, which we estimate to be twice the size of the rest-frame UV emission, yielding one of the first high-redshift objects where a clear signature of a [CII] ``Halo'' has been detected to date. For a matched S\'ersic profile with n = 1, we measured a [CII] effective radius of $0.50\pm 0.04$" (3.07$\pm 0.25$ kpc) and an average rest-frame UV effective radius of $0.2\pm0.04$" ($1.48\pm0.16$ kpc). The [CII] morphology and kinematics of the system suggest a merging event resulting in a non rotating disk system. This event could be responsible for the extended [CII] emission. Alternatively, some potential obscured emission could also explain the [CII] to UV size ratio. These results contribute to the growing consensus with respect to the existence of extended [CII] emission around galaxies.

A tantalizing enigma in extragalactic astronomy concerns the chronology and driving mechanisms of the build-up of late-type galaxies (LTGs). The standard scenario envisages two formation routes, with classical bulges (CBs) assembling first in a quick quasi-monolithic episode followed by gradual disk assembly, and pseudo-bulges (PBs) forming over the Gyr-long secular evolution of LTGs. The expectation is, therefore, the segregation of present-day LTG bulges into two distinct groups. Here we analyse the star formation histories (SFHs) of bulges and disks for 135 LTGs from the CALIFA survey covering the relevant range in LTG mass. In addition, their physical properties were contrasted with predictions from evolutionary synthesis models, adopting exponentially declining SFHs, with an e-folding time 0.1 < $\tau$ < 20 Gyr. Analysis of the SFHs of ~ half-million spaxels consistently reveals that the main properties of bulges and disks show a continuous distribution across total stellar mass. Moreover, the $\tau$ in high-mass LTGs radially increases, suggesting that these grow in an inside-out fashion, while lower-mass LTGs display roughly the same $\tau$ throughout their entire radial extent. Evolutionary synthesis predictions are consistent with observations. Finally, bulges and disks of higher mass LTGs exhibit shorter formation timescales as compared to their lower mass counterparts. Collectively, the obtained results evince a coherent and unified picture for the formation and evolution of LTGs, in which PBs and CBs denote extremities of a continuous sequence. This analysis is consistent with the framework where bulges are assembled with their parent disks by gradual inside-out growth, at a pace that is regulated by the depth of the galactic potential. In accordance is the utter absence of bimodal correlations, as expected if CBs and PBs were to emerge from two distinct formation routes.

This paper systematically studies the relation between metallicity and mass loss of massive stars. We perform one-dimensional stellar evolution simulations and build a grid of $\sim$2000 models with initial masses ranging between 11 and 60 $M_{\odot}$ and absolute metallicities $Z$ between 0.00001 and 0.02. Steady-state winds, comprising hot main-sequence winds and cool supergiant winds, are the main drivers of the mass loss of massive stars in our models. We calculate the total mass loss over the stellar lifetime for each model. Our results reveal the existence of a critical metallicity $Z_{\rm{c}}$ at $Z \sim 10^{-3}$, where the mass loss exhibits a dramatic jump. If $Z>Z_{\rm{c}}$, massive stars tend to evolve into cool supergiants, and a robust cool wind is operational. In contrast, if $Z<Z_{\rm{c}}$, massive stars usually remain as blue supergiants, wherein the cool wind is not activated and the mass loss is generally weak. Moreover, we calculate the wind feedback in a $10^5$ $M_{\odot}$ star cluster with the Salpeter initial mass function. The kinetic energy released by winds does not exhibit any significant transition at $Z_{\rm{c}}$ because the wind velocity of a cool supergiant wind is low and contributes little to the kinetic energy. The effects of critical metallicity provide implications for the fates of metal-poor stars in the early universe.

Using measurements of the [O III], H$\alpha$ and [N II] emission line fluxes originating in the cool (T $\sim10^4$ K) gas that populates the halos of massive early-type galaxies with stellar mass greater than $10^{10.4}$ M$_\odot$, we explore the recent conjecture that active galactic nucleus (AGN) activity preferentially removes the circumgalactic medium (CGM) along the polar (minor-axis) direction. We find deficits in the mean emission line flux of [O III] and H$\alpha$ (65 and 43%, respectively) along the polar vs. planar directions, although due to the large uncertainties in these difficult measurements the results are of marginal statistical significance (1.5$\sigma$). More robustly (97 to 99.9% confidence depending on the statistical test), diagnostic line ratios show stronger AGN ionization signatures along the polar direction at small radii than at other angles or radii. Our results are consistent with the conjecture of an anisotropic CGM in massive, early type galaxies, suggested on independent grounds, that is tied to AGN activity and begin to show the potential of CGM mapping using emission lines.

LUX-ZEPLIN (LZ) collaboration has achieved the strongest constraint on weak-scale dark matter (DM)-nucleon spin-independent (SI) scattering cross section in a large region of parameter space. In this paper, we take a complementary approach and study the prospect of detecting cosmic-ray boosted sub-GeV DM in LZ. In the absence of a signal for DM, we improve upon the previous constraints by a factor of $\sim 2$ using the LZ result for some regions of the parameter space. We also show that upcoming XENONnT and future Darwin experiments will be sensitive to cross sections smaller by factors of $\sim 3$ and $\sim 10$ compared to the current LZ limit, respectively.

We derive purely gravitational constraints on dark matter and cosmic neutrino profiles in the solar system using asteroid (101955) Bennu. We focus on Bennu because of its extensive tracking data and high-fidelity trajectory modeling resulting from the OSIRIS-REx mission. We find that the local density of dark matter is bound by $\rho_{\rm DM}\lesssim 3.3\times 10^{-15}\;\rm kg/m^3 \simeq 6\times10^6\,\bar{\rho}_{\rm DM}$, in the vicinity of $\sim 1.1$ au (where $\bar{\rho}_{\rm DM}\simeq 0.3\;\rm GeV/cm^3$). We show that high-precision tracking data of solar system objects can constrain cosmic neutrino overdensities relative to the Standard Model prediction $\bar{n}_{\nu}$, at the level of $\eta\equiv n_\nu/\bar{n}_{\nu}\lesssim 1.7 \times 10^{11}(0.1 \;{\rm eV}/m_\nu)$ (Saturn), comparable to the existing bounds from KATRIN and other previous laboratory experiments (with $m_\nu$ the neutrino mass). These local bounds have interesting implications for existing and future direct-detection experiments. Our constraints apply to all dark matter candidates but are particularly meaningful for scenarios including solar halos, stellar basins, and axion miniclusters, which predict or allow overdensities in the solar system. Furthermore, introducing a DM-SM long-range fifth force with a strength $\tilde{\alpha}_D$ times stronger than gravity, Bennu can set a constraint on $\rho_{\rm DM}\lesssim \bar{\rho}_{\rm DM}\left(6 \times 10^6/\tilde{\alpha}_D\right)$. These constraints can be improved in the future as the accuracy of tracking data improves, observational arcs increase, and more missions visit asteroids.

We study the prospects in the search of dark matter offered by the newly selected NASA MeV mission COSI (Compton Spectrometer and Imager). This instrument is designed and optimized to detect spectral lines, and we show it offers an exquisite possibility to detect dark matter directly decaying or annihilating into monochromatic gamma-rays. This is the case, for example, for axion-like particles (ALPs) which undergo decay into two photons. Furthermore, we show that COSI can lead to important progress in the quest for primordial black holes (PBHs) dark matter, through measurements of the 511 keV line from the positrons produced via Hawking evaporation. We also outline opportunities for the search of continuum signals, such as those expected from sub-GeV dark matter annihilation/decay into leptons and PBH evaporation into photons. We find that also in this case COSI can lead to improvements of current bounds.

We present single quasiparticle devices as new dark matter (DM) detectors. The threshold of these devices is set by the cooper pair binding energy, and is therefore so low that they can detect DM as light as about an MeV incoming from the Galactic halo, as well as the low-velocity thermalized DM component potentially present in the Earth. Using existing power measurements with these new devices, as well as power measurements with SuperCDMS-CPD, we set new constraints on the DM scattering cross section for DM masses from about 1 MeV to 10 GeV, down to about $10^{-34}-10^{-26}$ cm$^2$ for spin-independent interactions. We outline future directions to improve sensitivity to both halo DM and a thermalized DM population in the Earth using power deposition in quantum devices.

Gravitational-wave (GW) detectors that monitor fluctuations in the separation between inertial test masses (TMs) are sensitive to new forces acting on those TMs. Ultralight dark-photon dark matter (DPDM) coupled to $U(1)_B$ or $U(1)_{B-L}$ charges supplies one such force that oscillates with a frequency set by the DPDM mass. GW detectors operating in different frequency bands are thus sensitive to different DPDM mass ranges. A recent GW detection proposal based on monitoring the separation of certain asteroids in the inner Solar System would have sensitivity to $\mu$Hz frequencies [arXiv:2112.11431]. In this paper, we show how that proposal would also enable access to new parameter space for DPDM coupled to $B$ [respectively, $B-L$] charges in the mass range $4\ [8] \times 10^{-21} \text{eV} \lesssim m_{\text{DM}} \lesssim 2 \times 10^{-19} \text{eV}$, with peak sensitivities about 3 [2] orders of magnitude beyond current best limits on $\varepsilon_B$ [$\varepsilon_{B-L}$] at $m_{\text{DM}} \sim 2 \times 10^{-19} \text{eV}$. Sensitivity could be extended up to $m_{\text{DM}} \sim 3 \times 10^{-18} \text{eV}$ only if noise issues associated with asteroid rotational motion could be overcome.

Circumambient and galactic-scale environments are intermittently present around black holes, especially those residing in active galactic nuclei. As supermassive black holes impart energy on their host galaxy, so the galactic environment affects the geodesic dynamics of solar-mass objects around supermassive black holes and subsequently the gravitational waves emitted from such non-vacuum extreme-mass-ratio binaries. Only recently an exact general-relativistic solution has been found that describes a Schwarzschild black hole immersed in a dark matter halo profile of the Hernquist type. We perform an extensive geodesic analysis of test particles delving in such non-vacuum spacetimes and compare our results with those obtained in vacuum Schwarzschild spacetime, as well as their dominant gravitational-wave emission. Our findings indicate that the radial and polar oscillation frequency ratios, which indicate resonances, descend deeper into the extreme gravity regime as the compactness of the halo increases. This translates to a gravitational redshift of non-vacuum geodesics and their resulting waveforms with respect to the vacuum ones; a phenomenon which has also been observed for ringdown signals in these setups. For compact environments, we find that the apsidal precession of orbits is strongly affected due to the gravitational pull of dark matter; the orbit's axis can rotate in the opposite direction as that of the orbital motion, leading to a retrograde precession drift that depends on the halo's mass, as opposed to the typical prograde precession transpiring in vacuum and galactic-scale environments. Gravitational waves in retrograde-to-prograde orbital alterations demonstrate transient frequency phenomena around a critical non-precessing turning point, thus they may serve as a `smoking gun' for the presence of dense dark matter environments around supermassive black holes.

The distance-inclination degeneracy limits gravitational-wave parameter estimation of compact binary mergers. Although the degeneracy can be partially broken by including higher-order modes or precession, these effects are suppressed in binary neutron stars. In this work we implement a new parameterization of the tidal effects in the binary neutron star waveform, exploiting the binary Love relations, that breaks the distance-inclination degeneracy. The binary Love relations prescribe the tidal deformability of a neutron star as a function of its source-frame mass in an equation-of-state insensitive way, and thus allows direct measurement of the redshift of the source. If the cosmological parameters are assumed to be known, the redshift can be converted to a luminosity distance, and the distance-inclination degeneracy can thus be broken. We implement this new approach, studying a range of binary neutron-star observing scenarios using Bayesian parameter estimation on synthetic data. In the era of the third generation detectors, for observations with signal-to-noise ratios ranging from 6 to 167, we forecast up to a $\sim70\%$ decrease in the $90\%$ credible interval of the distance and inclination, and up to a $\sim50\%$ decrease in that of the source-frame component masses. For edge-on systems, our approach can result in moderate ($\sim50\%$) improvement in the measurements of distance and inclination for binaries with signal-to-noise ratio as low as 10. This prescription can be used to better infer the source-frame masses, and hence refine population properties of neutron stars, such as their maximum mass, impacting nuclear astrophysics. When combined with the search for electromagnetic counterpart observations, the work presented here can be used to put improved bounds on the opening angle of jets from binary neutron star mergers.

Using the horizontal neutral wind observations from the MIGHTI instrument onboard NASA's ICON (Ionospheric Connection Explorer) spacecraft with continuous coverage, we determine the climatology of the mean zonal and meridional winds and the associated mean circulation at low- to middle latitudes ($10^\circ$S-40$^{\circ}$N) for Northern Hemisphere {summer} solstice conditions between 90 km and 200 km altitudes, specifically on 20 June 2020 solstice as well as for a one-month period from 8 June-7 July 2020 {and for Northern winter season from 16 December 2019-31 January 2020, which spans a 47-day period, providing full local time coverage}. The data are averaged within appropriate altitude, longitude, latitude, solar zenith angle, and local time bins to produce mean wind distributions. The geographical distributions and local time variations of the mean horizontal circulation are evaluated. The instantaneous horizontal winds exhibit a significant degree of spatiotemporal variability often exceeding $\pm 150 $ m s$^{-1}$. The daily averaged zonal mean winds demonstrate day-to-day variability. Eastward zonal winds and northward (winter-to-summer) meridional winds are prevalent in the lower thermosphere, which provides indirect observational evidence of the eastward momentum deposition by small-scale gravity waves. The mean neutral winds and circulation exhibit smaller scale structures in the lower thermosphere (90-120 km), while they are more homogeneous in the upper thermosphere, indicating the increasingly dissipative nature of the thermosphere. The mean wind and circulation patterns inferred from ICON/MIGHTI measurements can be used to constrain and validate general circulation models, as well as input for numerical wave models.

Bulk viscosity in cold dark matter is an appealing feature that introduces distinctive phenomenological effects in the cosmological setting as compared to the $\Lambda$CDM model. Under this view, we propose a general parametrization of the bulk viscosity of the form $\xi\sim H^{1-2s} \rho_{m}^{s}$, that covers intriguingly some well-known cases in the Eckart's theory. Some advantages of this novel parametrization are: first, it allows to write the resulting equations of cosmological evolution in the form of an autonomous system for any value of $s$, so a general treatment of the fixed points and stability can be done, and second, the bulk viscosity effect is consistently handled so that it naturally turns off when matter density vanishes. As a main result we find, based on detailed dynamical system analysis, one-parameter family of de-Sitter-like asymptotic solutions with non-vanishing bulk viscosity coefficient during different cosmological periods. Numerical computations are performed jointly along with analytical phase space analysis in order to assess more quantitatively the bulk viscosity effect on the cosmological background evolution. Finally, as a first contact with observation we derive constraints on the free parameters of some bulk viscosity models with specific $s$-exponents from Supernovae Ia and observations of the Hubble parameter, by performing a Bayesian statistical analysis thought the Markov Chain Monte Carlo method.

As a supernova shock expands into space, it may collide with dark matter particles, scattering them up to velocities more than an order of magnitude larger than typical dark matter velocities in the Milky Way. If a supernova remnant is close enough to Earth, and the appropriate age, this flux of high-velocity dark matter could be detectable in direct detection experiments, particularly if the dark matter interacts via a velocity-dependent operator. This could make it easier to detect light dark matter that would otherwise have too little energy to be detected. We show that the Monogem Ring supernova remnant is both close enough and the correct age to produce such a flux, and thus we produce novel direct detection constraints and sensitivities for future experiments.

If dark energy is dynamical due to the evolution of a scalar field, then in general it is expected that the scalar is coupled to matter. While couplings to the standard model particles are highly constrained by local experiments, bounds on couplings to dark matter (DM) are only obtained from cosmological observations and they are consequently weaker. It has recently be pointed out that the coupling itself can become non-zero only at the time of dark energy domination, due to the evolution of dark energy itself, leading to a violation of the equivalence principle (EP) in the dark sector at late times. In this paper we study a specific model and show that such late-time violations of the EP in the DM sector are not strongly constrained by the evolution of the cosmological background and by observables in the linear regime (e.g. from the cosmic microwave background radiation). A study of perturbations in non-linear regime is necessary to constrain late--time violations of the equivalence principle much more strongly.

A modification to the Heisenberg uncertainty principle is called the generalized uncertainty principle (GUP), which emerged due to the introduction of a minimum measurable length, common among phenomenological approaches to quantum gravity. An approach to GUP called linear GUP (LGUP) has recently been developed that satisfies both the minimum measurable length and the maximum measurable momentum, resulting to a phase space volume proportional to the first-order momentum $(1 - \alpha p)^{-4} d^3x d^3p$, where $\alpha$ is the still-unestablished GUP parameter. In this study, we explore the mass-radius relations of LGUP-modified white dwarfs, and provide them with radial perturbations to investigate the dynamical instability arising from the oscillations. We find from the mass-radius relations that LGUP results to a white dwarf with a lower maximum mass, and this effect gets more apparent with larger the values of $\alpha$. We also observe that the mass of the white dwarf corresponding to the vanishing of the square of the fundamental frequency $\omega_0$ is the maximum mass the white dwarf can have in the mass-radius relations. The dynamical instability analysis also shows that instability sets in for all values of the GUP parameters $\alpha$, and at lower central densities $\rho_c$ (corresponding to lower maximum masses) for increasing $\alpha$, which verifies the results obtained from the mass-radius relations plots. Finally, we note that the mass limit is preserved for LGUP-modified white dwarfs, indicating that LGUP supports gravitational collapse of the compact object.

We study the mass-radius relations of finite temperature white dwarfs modified by the quadratic generalized uncertainty principle (QGUP), a prediction that arises from quantum gravity phenomenology. This QGUP approach extends the Heisenberg uncertainty principle by a quadratic term in momenta, which then modifies the phase space volume in the Chandrasekhar equation of state (EoS). This EoS was first calculated by treating the GUP parameter $\beta$ as perturbative. This perturbative EoS exhibits the expected thermal deviation for low pressures, while showing conflicting behaviors in the high pressure regime dependent on the sign of the $j$th order of approximation, $(\mathcal{O}(\beta^j))$. To explore the effects of QGUP further, we proceed with a full numerical simulation, and showed that in general, finite temperatures cause the EoS at low pressures to soften, while QGUP stiffens the EOS at high pressures. This modified EoS was then applied to the Tolman-Oppenheimer-Volkoff equations and its classical approximation to obtain the modified mass-radius relations for general relativistic and Newtonian white dwarfs. The relations for both cases were found to exhibit the expected thermal deviations at small masses, where low-mass white dwarfs are shifted to the high-mass regime at large radii, while high-mass white dwarfs acquire larger masses, beyond the Chandrasekhar limit. Additionally, we find that for sufficiently large values of the GUP parameter and temperature, we obtain mass-radius relations that are completely removed from the ideal case, as high-mass deviations due to GUP and low-mass deviations due to temperature are no longer mutually exclusive.

We study an axion-like particle (ALP) that experiences the first-order phase transition (FOPT) with respect to its mass or potential minimum. This can be realized if the ALP obtains a potential from non-perturbative effects of SU($N$) gauge theory that is confined via the FOPT, or if the ALP is trapped in a false vacuum at high temperatures until it starts to oscillate about the true minimum. The resulting ALP abundance is significantly enhanced compared to the standard misalignment mechanism, explaining dark matter in a broader parameter space that is accessible by e.g. IAXO and DM-radio. Furthermore, the viable parameter space includes a region of the mass $m_a \simeq 10^{-7} - 10^{-8}$ eV and the ALP-photon coupling $g_{a \gamma \gamma} \simeq 10^{-11} {\rm GeV}^{-1}$ that can explain the recent observation of very high energy photons from GRB221009A via axion-photon oscillations. If the ALP in this region explains dark matter, then the ALP has necessarily experienced a first-order phase transition.