Papers for Friday, Dec 09 2022

Over the last decade, exoplanetary transmission spectra have yielded an unprecedented understanding about the physical and chemical nature of planets outside our solar system. Physical and chemical knowledge is mainly extracted via fitting competing models to spectroscopic data, based on some goodness-of-fit metric. However, current employed metrics shed little light on how exactly a given model is failing at the individual data point level and where it could be improved. As the quality of our data and complexity of our models increases, there is an urgent need to better understand which observations are driving our model interpretations. Here we present the application of Bayesian leave-one-out cross-validation to assess the performance of exoplanet atmospheric models and compute the expected log pointwise predictive density (elpd$_\text{LOO}$). elpd$_\text{LOO}$ estimates the out-of-sample predictive accuracy of an atmospheric model at data point resolution providing interpretable model criticism. We introduce and demonstrate this method on synthetic HST transmission spectra of a hot Jupiter. We apply elpd$_\text{LOO}$ to interpret current observations of HAT-P-41b and assess the reliability of recent inferences of H$^-$ in its atmosphere. We find that previous detections of H$^{-}$ are dependent solely on a single data point. This new metric for exoplanetary retrievals complements and expands our repertoire of tools to better understand the limits of our models and data. elpd$_\text{LOO}$ provides the means to interrogate models at the single data point level, a prerequisite for robustly interpreting the imminent wealth of spectroscopic information coming from JWST.

We present the first cosmological constraints derived from the analysis of the void size function. This work relies on the final BOSS DR12 data set, a large spectroscopic galaxy catalog, ideal for the identification of cosmic voids. We extract a sample of voids from the distribution of galaxies and we apply a cleaning procedure aimed at reaching high levels of purity and completeness. We model the void size function by means of an extension of the popular volume-conserving model, based on two additional nuisance parameters. Relying on mock catalogs specifically designed to reproduce the BOSS DR12 galaxy sample, we calibrate the extended size function model parameters and validate the methodology. We then apply a Bayesian analysis to constrain the $\Lambda$CDM model and one of its simplest extensions, featuring a constant dark energy equation of state parameter, $w$. Following a conservative approach, we put constraints on the total matter density parameter and the amplitude of density fluctuations, finding $\Omega_{\rm m}=0.29^{+0.07}_{-0.06}$ and $\sigma_8=0.80^{+0.09}_{-0.08}$. Testing the alternative scenario, we derive $w=-1.1\pm 0.3$, in agreement with the $\Lambda$CDM model. These results are independent and complementary to those derived from standard cosmological probes, opening up new ways to identify the origin of potential tensions in the current cosmological paradigm.

Accurately accounting for spectral structure in spectrometer data induced by instrumental chromaticity on scales relevant for detection of the 21-cm signal is among the most significant challenges in global 21-cm signal analysis. In the publicly available EDGES low-band data set (Bowman et al. 2018), this complicating structure is suppressed using beam-factor based chromaticity correction (BFCC), which works by dividing the data by a sky-map-weighted model of the spectral structure of the instrument beam. Several analyses of this data have employed models that start with the assumption that this correction is complete; however, while BFCC mitigates the impact of instrumental chromaticity on the data, given realistic assumptions regarding the spectral structure of the foregrounds, the correction is only partial, which complicates the interpretation of fits to the data with intrinsic sky models. In this paper, we derive a BFCC data model from an analytic treatment of BFCC and demonstrate using simulated observations that the BFCC data model enables unbiased recovery of a simulated global 21-cm signal from beam factor chromaticity corrected data.

Diffuse interstellar bands (DIBs) are broad absorption features associated with interstellar dust and can serve as chemical and kinematic tracers. Conventional measurements of DIBs in stellar spectra are complicated by residuals between observations and best-fit stellar models. To overcome this, we simultaneously model the spectrum as a combination of stellar, dust, and residual components, with full posteriors on the joint distribution of the components. This decomposition is obtained by modeling each component as a draw from a high-dimensional Gaussian distribution in the data-space (the observed spectrum) -- a method we call "Marginalized Analytic Data-space Gaussian Inference for Component Separation" (MADGICS). We use a data-driven prior for the stellar component, which avoids missing stellar features not included in synthetic line lists. This technique provides statistically rigorous uncertainties and detection thresholds, which are required to work in the low signal-to-noise regime that is commonplace for dusty lines of sight. We reprocess all public Gaia DR3 RVS spectra and present an improved 8621 \r{A} DIB catalog, free of detectable stellar line contamination. We constrain the rest-frame wavelength to $8623.14 \pm 0.087$ \r{A} (vacuum), find no significant evidence for DIBs in the Local Bubble from the $1/6^{\rm{th}}$ of RVS spectra that are public, and show unprecedented correlation with kinematic substructure in Galactic CO maps. We validate the catalog, its reported uncertainties, and biases using synthetic injection tests. We believe MADGICS provides a viable path forward for large-scale spectral line measurements in the presence of complex spectral contamination.

Ram pressure stripping (RPS) is an important process to affect the evolution of cluster galaxies and their surrounding environment. We present a large MUSE mosaic for ESO 137-001 and its stripped tails, and studied the detailed distributions and kinematics of the ionized gas and stars. The warm, ionized gas is detected to at least 87 kpc from the galaxy and splits into three tails. There is a clear velocity gradient roughly perpendicular to the stripping direction, which decreases along the tails and disappears beyond $\sim45$ kpc downstream. The velocity dispersion of the ionized gas increases to $\sim80$ km s$^{-1}$ at $\sim20$ kpc downstream and stays flat beyond. The stars in the galaxy disk present a regular rotation motion, while the ionized gas is already disturbed by the ram pressure. Based on the observed velocity gradient, we construct the velocity model for the residual galactic rotation in the tails and discuss the origin and implication of its fading with distance. By comparing with theoretical studies, we interpreted the increased velocity dispersion as the result of the oscillations induced by the gas flows in the galaxy wake, which may imply an enhanced degree of turbulence there. We also compare the kinematic properties of the ionized gas and molecular gas from ALMA, which shows they are co-moving and kinematically mixed through the tails. Our study demonstrates the great potential of spatially resolved spectroscopy in probing the detailed kinematic properties of the stripped gas, which can provide important information for future simulations of RPS.

In the Jupyter ecosystem, data visualization is usually done with "widgets" created as notebook cell outputs. While this mechanism works well in some circumstances, it is not well-suited to presenting interfaces that are long-lived, interactive, and visually rich. Unlike the traditional Jupyter notebook system, the newer JupyterLab application provides a sophisticated extension infrastructure that raises new design possibilities. Here we present a novel user experience (UX) for interactive data visualization in JupyterLab that is based on an "app" that runs alongside the user's notebooks, rather than widgets that are bound inside them. We have implemented this UX for the AAS WorldWide Telescope (WWT) visualization tool. JupyterLab's messaging APIs allow the app to smoothly exchange data with multiple computational kernels, allowing users to accomplish tasks that are not possible using the widget framework. A new Jupyter server extension allows the frontend to request data from kernels asynchronously over HTTP, enabling interactive exploration of gigapixel-scale imagery in WWT. While we have developed this UX for WWT, the overall design and the server extension are portable to other applications and have the potential to unlock a variety of new user activities that aren't currently possible in "science platform" interfaces.

V2756 Sgr is a long-known S-type symbiotic binary that was recently detected in $\sim$ 1 mag eclipse-like fading by the Gaia satellite. This behavior was reported as a Gaia Science Alert under the designation Gaia22eor. V2756 Sgr has not been reported as an eclipsing symbiotic system. In this contribution, we have investigated the recent light curves of this target obtained by Gaia and ASAS-SN survey and supplemented these data by the photometry from the ASAS survey and DASCH archive. In addition, low-resolution BP/RP spectra from Gaia were examined. Based on the presented analysis, we conclude that the observed fading is indeed an eclipse of the hot component of V2756 Sgr by the cool giant. The data also allowed us to refine the orbital period to 725 $\pm$ 3 days.

We present timing solutions for 12 pulsars discovered in the Green Bank North Celestial Cap (GBNCC) 350 MHz pulsar survey, including six millisecond pulsars (MSPs), a double neutron star (DNS) system, and a pulsar orbiting a massive white dwarf companion. Timing solutions presented here include 350 and 820 MHz Green Bank Telescope data from initial confirmation and follow-up as well as a dedicated timing campaign spanning one year. PSR J1122$-$3546 is an isolated MSP, PSRs J1221$-$0633 and J1317$-$0157 are MSPs in black widow systems and regularly exhibit eclipses, and PSRs J2022+2534 and J2039$-$3616 are MSPs that can be timed with high precision and have been included in pulsar timing array experiments seeking to detect low-frequency gravitational waves. PSRs J1221$-$0633 and J2039$-$3616 have Fermi Large Area Telescope $\gamma$-ray counterparts and also exhibit significant $\gamma$-ray pulsations. We measure proper motion for three of the MSPs in this sample and estimate their space velocities, which are typical compared to those of other MSPs. We have detected the advance of periastron for PSR J1018$-$1523 and therefore measure the total mass of the double neutron star system, $m_{\rm tot}=2.3\pm0.3$ M$_{\odot}$. Long-term pulsar timing with data spanning more than one year is critical for classifying recycled pulsars, carrying out detailed astrometry studies, and shedding light on the wealth of information in these systems post-discovery.

The two primary observable quantities of an exoplanet--its mass and radius--alone are not sufficient to probe a rocky exoplanet's interior composition and mineralogy. To overcome this, host-star abundances of the primary planet-building elements (Mg, Si, Fe) are typically used as a proxy for the planet's bulk composition. The majority of small exoplanet hosts, however, do not have available abundance data. Here we present the open-source ExoPlex mass-radius-composition solver. Unlike previous open-source mass-radius solvers, ExoPlex calculates the core chemistry and equilibrium mantle mineralogy for a bulk composition, including effects of mantle FeO content, core light elements and surface water/ice. We utilize ExoPlex to calculate the planetary radii, surface gravities and bulk densities for 10$^6$ model planets up to 2 R$_\oplus$ across these geochemistries, adopting the distribution of FGK stellar abundances to estimate of the range of bulk exoplanet compositions. We outline the $99.7\%$ distribution of radii, surface gravity and bulk densities that define planets as "nominally rocky." Planets outside this range require compositions outside those expected from stellar abundance data, likely making them either Fe-enriched super-Mercuries, or volatile-enriched mini-Neptunes. We apply our classification scheme to a sample of 85 well-resolved exoplanets without available host-star abundances. We estimate only 9 planets are within the "nominally rocky planet zone" at $>70\%$ confidence, while $\sim20\%$ and $\sim30\%$ of this sample can be reasonably classified as super-Mercuries or volatile-rich, respectively. Our results provide observers with a self-consistent way to broadly classify a planet as likely rocky, Mercury-like or volatile-enriched, using mass and radius measurements alone.

PSR J1101-6101 is an energetic young pulsar which powers the remarkable Lighthouse pulsar wind nebula (PWN). The pulsar belongs to the rare type of radio- and gamma-ray-quiet pulsars which are bright in hard X-rays. Moreover, the Lighthouse PWN is remarkable for its misaligned outflow (which gave rise to the PWN's nickname). Also known as "pulsar filaments", these collimated parsec-scale X-ray structures have been recently discovered in the vicinity of a handful of fast-moving pulsars, and appear unaffected by the ram pressure which confines pulsar tails. We report on NuSTAR observations of PSR J1101-6101 and its misaligned outflow -- the first observation of such a structure above ~10 keV. We detect the outflow up to 25 keV, spatially resolve its spectral evolution with distance from the pulsar, find unambiguous evidence of spectral cooling with distance from the pulsar, and infer physical properties of the particles and magnetic field in the outflow. We also reanalzye archival Chandra data and discuss the outflow's small-scale structure. We detect pulsations from PSR J1101-6101 up to 20 keV, present the X-ray pulse profile, confirm its period derivative, and perform phase-resolved spectroscopy. Lastly, we discuss the X-ray source 2CXO J110158.4-605649 = 2XMM J110158.5--605651 (a serendipitously observed blazar) and suggest it may be the X-ray counterpart to the GeV source 4FGL J1102.0--6054.

A planet's orbital alignment places important constraints on how a planet formed and consequently evolved. The dominant formation pathway of ultra-short period planets ($P<1$ day) is particularly mysterious as such planets most likely formed further out, and it is not well understood what drove their migration inwards to their current positions. Measuring the orbital alignment is difficult for smaller super-Earth/sub-Neptune planets, which give rise to smaller amplitude signals. Here we present radial velocities across two transits of 55 Cancri e, an ultra-short period Super-Earth, observed with the Extreme Precision Spectrograph (EXPRES). Using the classical Rossiter-McLaughlin (RM) method, we measure 55 Cnc e's sky-projected stellar spin-orbit alignment (i.e., the projected angle between the planet's orbital axis and its host star's spin axis) to be $\lambda=10\substack{+17\\ -20}^{\circ}$ with an unprojected angle of $\psi=23\substack{+14\\ -12}^{\circ}$. The best-fit RM model to the EXPRES data has a radial velocity semi-amplitude of just $0.41\substack{+0.09\\ -0.10} m s^{-1}$. The spin-orbit alignment of 55 Cnc e favors dynamically gentle migration theories for ultra-short period planets, namely tidal dissipation through low-eccentricity planet-planet interactions and/or planetary obliquity tides.

We explore a possible scenario of the explosion as a result of core collapses of rotating massive stars that leave a black hole by performing a radiation-viscous-hydrodynamics simulation in numerical relativity. We take moderately and rapidly rotating compact pre-collapse stellar models derived in stellar evolution calculations as the initial conditions. We find that the viscous heating in the disk formed around the central black hole powers an outflow. For rapidly rotating models, the explosion energy is $\gtrsim 3\times10^{51}$ erg, which is comparable to or larger than that of typical stripped-envelope supernovae, indicating that a fraction of such supernovae may be explosions powered by black-hole accretion disks. The explosion energy is still increasing at the end of the simulations with a rate of $>10^{50}$ erg/s, and thus, it may reach $\sim10^{52}$ erg. The nucleosynthesis calculation shows that the mass of $^{56}$Ni amounts to $\gtrsim 0.1M_\odot$, which, together with the high explosion energy, satisfies the required amount for broad-lined type Ic supernovae. The moderately rotating models predict small ejecta mass of order $0.1M_\odot$ and explosion energy of $\lesssim 10^{51}$ erg. Due to the small ejecta mass, these models may predict a short-timescale transient with the rise time 3$-$5 d. It can lead to a bright ($\sim10^{44}$ erg/s) transient like superluminous supernovae in the presence of dense massive circum-stellar medium. Irrespective of the models, the lowest value of the electron fraction of the ejecta is $\gtrsim 0.4$, and thus, the synthesis of heavy $r$-process elements is not found in our calculation.

Aging GONG second generation cameras (Silicon Mountain Design(TM) cameras) were planned to be replaced after their long service of more than a decade. This prompted a market-wide search for a potential replacement detector to meet the GONG science requirements. This report provides some history of the search process, a comparison between CMOS and CCD type sensors and then a quantitative evaluation of potential candidates to arrive at final selection. Further, a feasibility study of the selected sensor for adaptation to GONG optical system was done and sensor characteristics were independently verified in the laboratory. This technical report gives description of these studies and tests.

Fast radio bursts (FRBs), millisecond-duration bursts prevailing in the radio sky, are the latest big puzzle in the universe and have been a subject of intense observational and theoretical investigations in recent years. The rapid accumulation of the observational data has painted the following sketch about the physical origin of FRBs: They predominantly originate from cosmological distances so that their sources produce the most extreme coherent radio emission in the universe; at least some, probably most, FRBs are repeating sources that do not invoke cataclysmic events; and at least some FRBs are produced by magnetars, neutron stars with the strongest magnetic fields in the universe. Many open questions regarding the physical origin(s) and mechanism(s) of FRBs remain. This article reviews the phenomenology and possible underlying physics of FRBs. Topics include: a summary of the observational data, basic plasma physics, general constraints on FRB models from the data, radiation mechanisms, source and environment models, propagation effects, as well as FRBs as cosmological probes. Current pressing problems and future prospects are also discussed.

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.

In spite of the rich phenomenology of the polarization properties of radio pulsars, the rotating vector model (RVM) created 50 years ago remains the best method to determine the beam geometry of a pulsar. We apply the RVM to a sample of 854 radio pulsars observed with the MeerKAT telescope in order to draw conclusions about the population of pulsars as a whole. The main results are that (i) the geometrical interpretation of the position angle traverse is valid in the majority of the population, (ii) the pulsars for which the RVM fails tend to have a high fraction of circular polarization compared to linear polarization, (iii) emission heights obtained through both geometrical and relativistic methods show that the majority of pulsars must have emission heights less than 1000~km independent of spin period, (iv) orthogonal mode jumps are seen in the position angle traverse in about one third of the population. All these results are weakly dependent on the pulsar spin-down energy.

Observations of center-to-limb variations (CLV) of spectral lines and continua provide a good test for the accuracy of models of solar and stellar atmospheric structure and spectral-line formation. They are also widely used to constrain elemental abundances and are becoming increasingly important in atmospheric studies of exoplanets. However, only a few such data sets exist for chromospheric lines. We aim to create a set of standard profiles by means of mosaics made with the CRISP and CHROMIS instruments of the Swedish 1-m Solar Telescope (SST), as well as to explore the robustness of said profiles obtained using this method. For each spectral line we use a mosaic that ranges from the center to the limb. Each of these mosaics are averaged down to 50 individual spectral profiles, spaced by 0.02 in the $\mu$ scale. These profiles are corrected for p-mode oscillations, and their line parameters (equivalent width, line shift, full-width at half-maximum, and line depth) are then compared against literature values where possible. We present a set of 50 average profiles that are spaced equidistantly along the cosine of the heliocentric angle ($\mu$) by steps of 0.02 for five continuum points between 4001 and 7772 \AA, as well as ten of the most commonly observed spectral lines at the SST (Ca II H & K, H$\beta$, Mg I 5173 \AA, C I 5380 \AA, Fe I 6173 \AA, Fe I 6301 \AA, H$\alpha$, O I 7772 \AA, and Ca II 8542 \AA). Center-to-limb variation of line profiles and continua are shared in the CDS as machine-readable tables; providing a quantitative constraint on theoretical models that aim to model stellar atmospheres.

The morphological classification of galaxies is considered a relevant issue and can be approached from different points of view. The increasing growth in the size and accuracy of astronomical data sets brings with it the need for the use of automatic methods to perform these classifications. The aim of this work is to propose and evaluate a method for automatic unsupervised classification of kinematic morphologies of galaxies that yields a meaningful clustering and captures the variations of the fundamental properties of galaxies. We obtain kinematic maps for a sample of 2064 galaxies from the largest simulation of the EAGLE project that mimics integral field spectroscopy (IFS) images. These maps are the input of a dimensionality reduction algorithm followed by a clustering algorithm. We analyse the variation of physical and observational parameters among the clusters obtained from the application of this procedure to different inputs. The inputs studied in this paper are (a) line-of-sight velocity maps for the whole sample of galaxies observed at fixed inclinations, (b) line-of-sight velocity, dispersion, and flux maps together for the whole sample of galaxies observed at fixed inclinations, (c) line-of-sight velocity, dispersion, and flux maps together for two separate subsamples of edge-on galaxies with similar amount of rotation, and (d) line-of-sight velocity, dispersion, and flux maps together for galaxies from different observation angles mixed. The application of the method to solely line-of-sight velocity maps achieves a clear division between slow rotators (SRs) and fast rotators (FRs) and can differentiate rotation orientation. By adding the dispersion and flux information at the input, low rotation edge-on galaxies are separated according to their shapes. Abridged.

Coronal mass ejections (CMEs) result from eruptions of magnetic flux ropes (MFRs) and can possess a three-part structure in white-light coronagraphs, including a bright front, dark cavity and bright core. In the traditional opinion, the bright front forms due to the plasma pileup along the MFR border, the cavity represents the cross section of the MFR, and the bright core corresponds to the erupted prominence. However, this explanation on the nature of the three-part structure is being challenged. In this paper, we report an intriguing event occurred on 2014 June 14 that was recorded by multiple space- and ground-based instruments seamlessly, clearly showing that the CME front originates from the plasma pileup along the magnetic arcades overlying the MFR, and the core corresponds to a hot-channel MFR. Thus the dark cavity is not an MFR, instead it is a low-density zone between the CME front and a trailing MFR. These observations are consistent with a new explanation on the CME structure. If the new explanation is correct, most (if not all) CMEs should exhibit the three-part appearance in their early eruption stage. To examine this prediction, we make a survey study of all CMEs in 2011 and find that all limb events have the three-part feature in the low corona, regardless of their appearances in the high corona. Our studies suggest that the three-part structure is the intrinsic structure of CMEs, which has fundamental importance for understanding CMEs.

We use high-resolution zoom-in cosmological simulations to model outflow triggered by radiation and thermal drivers around the central mass accumulation during direct collapse within the dark matter (DM) halo. The maximal resolution is $8\times 10^{-5}$ pc, and no restrictions are put on the geometry of the inflow/outflow. The central mass is considered {\it prior} to the formation of the supermassive black hole seed at redshift of $z\sim 15.9$, and can constitute either a supermassive star (SMS) of $\sim 10^5\,M_\odot$ surrounded by growing accretion disk or a self-gravitating disk. The radiation transfer is modeled using the ray-tracing algorithm. Due to the high accretion rate of $\sim 1\,M_\odot\,{\rm yr^{-1}}$ determined by the DM halo, accretion is supercritical, resulting in supercritical luminosity which affects the inflow rate, with the duty cycle of $\sim 0.9$. We observe a fast development of hot cavities which quickly extend into polar funnels and expanding dense shells. Within the funnels, fast winds, $\sim 10^3\,{\rm km\,s^{-1}}$, are mass-loaded by the accreting gas. We follow the expanding shells to $\sim 1$ pc, when the shell velocity remains substantially, $\sim 5$ times, above the escape speed. The ionization cones formed by the central UV/X-ray flux extend to the halo virial radius, $R_{\rm h}$. Extrapolating the outflow properties shows that the halo material outside the shell will have difficulty to stop it. We therefore conclude that the expanding wind-driven shell will breakout of the central parsec and will reach $\sim R_{\rm h}$. Finally, the anisotropic accretion flow on sub-parsec scales will attenuate the soft X-rays effect on the H$_2$. Hence, formation of funnels and powerful outflows around, e.g., SMS, can have interesting observational corollaries, to be addressed elsewhere.

Fast radio bursts (FRBs) are millisecond-duration transients with large dispersion measures. The origin of FRBs is still mysterious. One of the methods to comprehend FRB origin is to probe the physical environments of FRB host galaxies. Mapping molecular-gas kinematics in FRB host galaxies is critical because it results in star formation that is likely connected to the birth of FRB progenitors. However, most previous works of FRB host galaxies have focused on its stellar component. Therefore, we, for the first time, report the molecular gas kinematics in the host galaxy of the non-repeating FRB 180924B at $z= 0.3216$. Two velocity components of the CO (3-2) emission line are detected in its host galaxy with the Atacama Large Millimeter/submillimeter Array (ALMA): the peak of one component ($-155.40$ km s$^{-1}$) is near the centre of the host galaxy, and another ($-7.76$ km s$^{-1}$) is near the FRB position. The CO (3-2) spectrum shows asymmetric profiles with A$_{\rm peak}$ $=2.03\pm 0.39$, where A$_{\rm peak}$ is the peak flux density ratio between the two velocity components. The CO (3-2) velocity map also indicates an asymmetric velocity gradient from $-180$ km s$^{-1}$ to 8 km s$^{-1}$. These results indicate a disturbed kinetic structure of molecular gas in the host galaxy. Such disturbed kinetic structures are reported for repeating FRB host galaxies using HI emission lines in previous works. Our finding indicates that non-repeating and repeating FRBs could commonly appear in disturbed kinetic environments, suggesting a possible link between the gas kinematics and FRB progenitors.

We report on the identification of a new $\gamma$-ray emitting narrow-line Seyfert 1 galaxy ($\gamma$-NLS1), SDSS J095909.51+460014.3 (hereinafter J0959+4600, $z$ = 0.399), by establishing its association with a $\gamma$-ray source 4FGL 0959.6+4606, although its low-energy counterpart was suggested to be a radio galaxy 2MASX J09591976+4603515 (hereinafter J0959+4603). \emph{WISE} long-term light curves of these two sources reveal diverse infrared variability patterns. Violent infrared variations of J0959+4600 with an amplitude up to one order of magnitude has been detected, while variability is mild for the other one. More importantly, the peaking time of the infrared flare of J0959+4603 is coincident with the time of brightening of 4FGL 0959.6+4606. At the same time, the flux level of J0959+4603 maintains in the quiescent state. A specific analysis of 15-month Fermi-LAT data, aiming at the high $\gamma$-ray flux state, yields a significant source (TS =41). The corresponding localization analysis suggests that only J0959+4600 falls into the uncertainty area, further supporting the updated association relationship. A broadband spectral energy distribution of J0959+4600 has been drawn and well described by the classic single-zone homogeneous leptonic jet model. The jet properties of J0959+4600 are investigated and found to be comparable with other $\gamma$-NLS1s.

LHAASO J1908+0621 has been recently detected as a source emitting $\gamma$-rays with energies above 100 TeV, and the multiband observations show that a break around 1 TeV appears in the $\gamma$-ray spectrum. We have reanalyzed the GeV $\gamma$-ray properties for the 100 TeV source using 14 years of data recorded by the Fermi Large Area Telescope (Fermi-LAT). The spectrum in the energy range range 30-500 GeV has an index of 1.50 $\pm$ 0.26, which is much smaller than that detected in the TeV $\gamma$-rays. Additionally, the radiation properties of this source are investigated based on a one-zone time-dependent model. In the model, LHAASO J1908+0621 is associated with a pulsar wind nebula (PWN) powered by the pulsar PSR J1907$+$0602. High-energy particles composed of electrons and positrons are injected into the nebula. Multiband nonthermal emission is produced via synchrotron radiation and inverse Compton scattering (ICS). Taking the effect of radiative energy losses and adiabatic cooling into account, the spectral energy distribution from the model with a broken power-law for the distribution of the injected particles can explain the detected fluxes in the $\gamma$-ray bands. The results support that LHAASO J1908+0621 originates from the PWN powered by PSR J1907$+$0602, and the $\gamma$-rays with energy above 100 TeV are produced by the electrons/positrons in the nebula via ICS.

The prototype dwarf nova SS Cyg unexpectedly exhibited an anomalous event in its light curve in the early few months of 2021 in which regular dwarf nova-type outbursts stopped, but small-amplitude fluctuations occurred only. Inspired by this event, we have performed numerical simulations of light curves of SS Cyg by varying mass transfer rates and varying viscosity parameters in the cool disk. We have also studied the effect of gas-stream overflows beyond the outer disk edge in the light curve simulations. We have confirmed that the enhanced mass transfer is unlikely responsible for the 2021 anomalous event and its forerunner. We have found that the enhancement of the viscosity in the disk may reproduce the forerunner of that event but may not be enough to explain the 2021 anomalous event, although the latter result might be particular to our thermal equilibrium curve used. Within our simulations, a model of the gas stream overflow with a slightly higher mass transfer rate than that of our standard model reproduces light curves similar to the 2021 anomalous event. We suggest that the gas-stream overflow is necessary to reproduce that event. The gas-stream overflow may also be responsible for the abnormally high X-ray flux during the normal quiescent state in SS Cyg.

The discrepancy between the mass of galaxies and their rotational velocity is one of the most puzzling scientific phenomena. Despite over a century of research, this phenomenon is not fully understood. Common explanations include dark matter and MOND, among other theories. Here we report on another observation that shows tension between the physics of galaxy rotation and its rotational velocity. We compare the brightness of galaxies, and find that galaxies that spin in the same direction as the Milky Way have different brightness than galaxies that spin in the opposite direction. While such difference in brightness is expected due to Doppler shift, it is expected to be subtle. The results show that the difference in brightness is large enough to be detected by Earth-based telescopes. That observed difference corresponds to physical properties of galaxies with far greater rotational velocity than the rotational velocity of the Milky Way. The difference is consistent in both the Northern galactic pole and the Southern galactic pole, and is not observed in parts of the sky that are perpendicular to the galactic pole. The differences are observed by several different instruments such DECam, SDSS, Pan-STARRS, and HST. The observation is also consistent across annotation methods, including different computer-based methods, manual annotation, or crowdsourcing annotations through Galaxy Zoo, all show similar results. Another possible explanation to the observation is parity violation in the large-scale structure, such that the magnitude of the parity violation was stronger in the earlier Universe. It can also be linked to other anomalies such as the Ho tension. Analysis of Ho using Ia supernovae shows smaller Ho tension when the spin directions of the host galaxies are consistent, although these results are based on a small number of supernovae, and may not be statistically significant.

Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope that show evidence of Shapiro delay: PSRs~J0101$-$6422, J1101$-$6424, J1125$-$6014, J1514$-$4946, J1614$-$2230, J1732$-$5049, and J1909$-$3744. Evidence for Shapiro delay was found in all of the systems, and for three the orientations and data quality enabled strong constraints on their orbital inclinations and component masses. For PSRs~J1125$-$6014, J1614$-$2230 and J1909$-$3744, we determined pulsar masses to be $M_{\rm p} = 1.68\pm 0.17 \, {\rm M_{\odot}}$, $1.94\pm 0.03 \, {\rm M_{\odot}}$ and $1.45 \pm 0.03 \, {\rm M_{\odot}}$, and companion masses to be $M_{\rm c} = 0.33\pm 0.02 \, {\rm M_{\odot}}$, $0.495\pm 0.005 \, {\rm M_{\odot}}$ and $0.205 \pm 0.003 \, {\rm M_{\odot}}$, respectively. This provides the first independent confirmation of PSR~J1614$-$2230's mass, one of the highest known. The Shapiro delays measured for PSRs~J0101$-$6422, J1101$-$6424, J1514$-$4946, and J1732$-$5049 were only weak, and could not provide interesting component mass limits. Despite a large number of millisecond pulsars being routinely timed, relatively few have accurate masses via Shapiro delays. We use simulations to show that this is expected, and provide a formula for observers to assess how accurately a pulsar mass can be determined. We also discuss the observed correlation between pulsar companion masses and spin period, and the anti-correlation between recycled pulsar mass and their companion masses.

Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt through either dynamical or accretion signatures. Gravitational waves (GWs) represent an unparalleled opportunity to survey the sky and detect mergers of IMBHs, making it possible for the first time to constrain their formation, growth, and merger history across cosmic time. While the current network LIGO-Virgo-KAGRA is significantly limited in detecting mergers of IMBH binaries, the next generation of ground-based observatories and space-based missions promise to shed light on the IMBH population through the detection of several events per year. Here, we asses this possibility by determining the optimal network of next-generation of GW observatories to reconstruct the IMBH merger history across cosmic time. We show that Voyager, the Einstein Telescope, and Cosmic Explorer will be able to constrain the distribution of the primary masses of merging IMBHs up to $\sim 10^3\ M_\odot$ and with mass ratio $\gtrsim 0.1$, while LISA will complementary do so at higher mass and smaller mass ratios. Therefore, a network of next-generation ground-based and space-based observatories will potentially reconstruct the merger history of IMBHs. Moreover, IMBHs with masses $\lesssim 5\times 10^3\ M_\odot$ could be observed in multiband up to a redshift of $z\approx 4$, ushering in a new of era GW astronomy.

In this paper, we study the power-law $f(T)$ model using Hubble diagrams of type Ia supernovae (SNIa), quasars (QSOs), Gamma Ray Bursts (GRBs) and the measurements from baryonic acoustic oscillations (BAO) in the framework of the cosmographic method. Using mock data for SNIa, QSOs and GRBs generated based on the power-law $f(T)$ model, we show whether different cosmographic methods are suitable to reconstruct the distance modulus or not. In particular, we investigate the rational PADE polynomials $(3,2)$ and $(2,2)$ in addition to the fourth- and fifth- order Taylor series. We show that PADE $(3,2)$ is the best approximation that can be used in the cosmographic method to reconstruct the distance modulus at both low and high redshifts. In the context of PADE $(3,2)$ cosmographic method, we show that the power-law $f(T)$ model is well consistent with the real observational data from the Hubble diagrams of SNIa, QSOs and GRBs. Moreover, we find that the combination of the Hubble diagram of SNIa and the BAO observation leads to better consistency between the model-independent cosmographic method and the power-law $f(T)$ model. Finally, our observational constraints on the parameter of the effective equation of state of DE, described by the power-law $f(T)$ model, show the phantom-like behavior, especially when the BAO observations are included in our analysis.

The ability to test and constrain theories of cosmic inflation will advance substantially over the next decade. Key data sources include cosmic microwave background (CMB) measurements and observations of the distribution of matter at low-redshift from optical, near-infrared, and 21cm intensity surveys. A positive detection of a CMB B-mode consistent with a primordial stochastic gravitational wave background (SGWB) is widely viewed as a smoking gun for an inflationary phase. Still, a null result does not exclude inflation. However, in a significant class of inflationary scenarios, a low SGWB amplitude is correlated with a more significant running, $\alpha_s$, in the primordial density perturbations than is seen with the simplest inflationary potentials. With this motivation, we forecast the precision with which the spectral index $n_{\rm{s}}$ and $\alpha_{\rm{s}}$ can be constrained by currently envisaged observations, including CMB (Simons Observatory, CMB-S4 and \textit{LiteBIRD}), optical/near infra-red (DESI and SPHEREx), and 21cm intensity mapping (Tianlai and CHIME) surveys. We identify optimal combinations of datasets for constraining the running and show that they may yield additional and informative constraints on the overall inflationary parameter space if the SGWB remains undetected.

Using mock data for the Hubble diagrams of type Ia supernovae (SNIa) and quasars (QSOs) generated based on the standard model of cosmology, and using the least-squares method based on the Markov-Chain-Monte-Carlo (MCMC) algorithm, we first put constraints on the cosmographic parameters in the context of the various model-independent cosmographic methods reconstructed from the Taylor $4^{th}$ and $5^{th}$ order expansions and the Pade (2,2) and (3,2) polynomials of the Hubble parameter, respectively. We then reconstruct the distance modulus in the framework of cosmographic methods and calculate the percentage difference between the distance modulus of the cosmographic methods and that of the standard model. The percentage difference is minimized when the Pade approximation is used which means that the Pade cosmographic method is sufficiently suitable for reconstructing the distance modulus even at high-redshifts. In the next step, using the real observational data for the Hubble diagrams of SNIa, QSOs, gamma-ray-bursts (GRBs), and observations from baryon acoustic oscillations (BAO) in two sets of the low-redshift combination (SNIa+QSOs+GRBs+BAO) embracing the redshift range of $0.01<z<2.26$ and the high-redshift combination (SNIa+QSOs+GRBs) which covers a redshift range of $0.01< z < 5.5$, we put observational constraints on the cosmographic parameters of the Pade cosmography and also the standard model. Our analysis indicates that Pade cosmographic approaches do not reveal any cosmographic tension between the standard model and the observational data. We also confirm this result, using the statistical AIC criteria. Finally, we put the cosmographic method in the redshift-bin data and find a larger value of $\Omega_{m0}$ extracted from $s_0$ parameter compared with those of the $q_0$ parameter and Planck-$\Lambda$CDM values.

The XL-Calibur balloon-borne hard X-ray polarimetry mission comprises a Compton-scattering polarimeter placed at the focal point of an X-ray mirror. The polarimeter is housed within a BGO anticoincidence shield, which is needed to mitigate the considerable background radiation present at the observation altitude of ~40 km. This paper details the design, construction and testing of the anticoincidence shield, as well as the performance measured during the week-long maiden flight from Esrange Space Centre to the Canadian Northwest Territories in July 2022. The in-flight performance of the shield followed design expectations, with a veto threshold <100 keV and a measured background rate of ~0.5 Hz (20-40 keV). This is compatible with the scientific goals of the mission, where %-level minimum detectable polarisation is sought for a Hz-level source rate.

The thesis concerns the analysis of the Active Galactic Nuclei. These are galaxies with an active core. The most luminous type of Active Galactic Nuclei is Quasar. It contains the supermassive black hole at the center. One of the least known subtype of Quasars are: Weak emission-Line Quasars. Their recognizable feature are weak emission-lines. The primary goal of PhD thesis is to evaluate the global parameters such as: the black hole mass, the accretion rate, spin of the black hole, and the inclination of weak emission-line quasars based on the continuum fit method. This method apart from the literature black hole masses estimation methods does not depend on the observed Full Width at Half Maximum of emission line, which could be biased due to the weakness or lack of the emission lines in these quasars. Using the Spectral Energy Distribution of quasars, I have fitted the geometrically thin and optically thick accretion disk model described by Novikov \& Thorne equations. I have obtained the model of the continuum of the accretion disk for the 10 weak emission-line quasars. The second project concerned the description of abnormal, deep absorption of SDSS J110511.15+530806.5 quasar. I checked the correctness of the thesis posed that corona and warm skin concept above/around an accretion disk explain this phenomenon.

Diagnosing the spatial-temporal pattern of magnetic flux on the Sun is vital for understanding the origin of solar magnetism and activity. Here, we report a new form of flux appearance, magnetic outbreak, using observations with an extremely high spatial resolution of 0.16 arcsec from the 1.6-m Goode Solar Telescope (GST) at the Big Bear Solar Observatory. Magnetic outbreak refers to an early growth of unipolar magnetic flux and its later explosion into fragments, in association with plasma upflow and exploding granulations; each individual fragment has flux of 10$^{16}$-10$^{17}$ Mx, moving apart with velocity of 0.5-2.2 km/s. The magnetic outbreak takes place in the hecto-Gauss region of pore moats. In this study, we identify six events of magnetic outbreak during 6-hour observations over an approximate 40$\times$40 arcsec$^{2}$ field of view. The newly discovered magnetic outbreak might be the first evidence of the long-anticipated convective blowup.

We explore the processes of repetitive build-up and explosive release of magnetic energy together with the formation of magnetic flux ropes that eventually resulted into three homologous eruptive flares of successively increasing intensities (i.e., M2.0, M2.6, and X1.0). The flares originated from NOAA active region 12017 during 2014 March 28-29. EUV observations and magnetogram measurements together with coronal magnetic field modeling suggest that the flares were triggered by the eruption of flux ropes embedded by a densely packed system of loops within a small part of the active region. In X-rays, the first and second events show similar evolution with single, compact sources, while the third event exhibits multiple emission centroids with a set of strong non-thermal conjugate sources at 50-100 keV during the HXR peak. The photospheric magnetic field over an interval of approximately 44 hr encompassing the three flares undergoes important phases of emergence and cancellation processes together with significant changes near the polarity inversion lines within the flaring region. Our observations point toward the tether-cutting mechanism as the plausible triggering process of the eruptions. Between the second and third event, we observe a prominent phase of flux emergence which temporally correlates with the build-up phase of free magnetic energy in the active region corona. In conclusion, our analysis reveals an efficient coupling between the rapidly evolving photospheric and coronal magnetic fields in the active region that led to a continued phase of the build-up of free energy, resulting into the homologous flares of successively increasing intensities.

We study the acceleration of ultra-high-energy cosmic rays (UHECRs) at FR-II radio galaxies by performing Monte Carlo simulations for the transport, scattering, and energy change of the CR particles injected into the time-evolving jet flows that are realized through relativistic hydrodynamic (RHD) simulations. Toward that end, we adopt physically motivated models for the magnetic field and particle scattering. By identifying the primary acceleration process among diffusive shock acceleration (DSA), turbulent shear acceleration (TSA), and relativistic shear acceleration (RSA), we find that CRs of $E\lesssim1$ EeV gain energy mainly through DSA in the jet-spine flow and the backflow containing many shocks and turbulence. After they attain $E\gtrsim$ a few EeV, CRs are energized mostly via RSA at the jet-backflow interface, reaching energies well above $10^{20}$ eV. TSA makes a relatively minor contribution. The time-asymptotic energy spectrum of escaping particles is primarily governed by the jet power, shifting to higher energies at more powerful jets. The UHECR spectrum fits well to the double-power-law form, whose break energy, $E_{\rm break}$, corresponds to the size-limited maximum energy. It is close to $d\mathcal{N}/dE\propto E^{-0.5}$ below $E_{\rm break}$, while it follows $d\mathcal{N}/dE\propto E^{-2.6}$ above $E_{\rm break}$, decreasing more gradually than the exponential. The power-law slope of the high-energy end is determined by the energy boosts via non-gradual shear acceleration across the jet-backflow interface and the confinement by the elongated cocoon. We conclude that giant radio galaxies could be major contributors to the observed UHECRs.

We analyze features of the Fe K spectrum of the high mass X-ray binary Cygnus X-3. The spectrum was obtained with the Chandra High Energy Transmission Grating Spectrometer in the third diffraction order. The increased energy resolution of the third order enables us to fully resolve the Fe XV He$\alpha$ complex, the Fe XVI Ly$\alpha$ lines and features of the radiative recombination continua. The emission line spectrum shows the expected features of photoionization equilibrium, excited in the dense stellar wind of the companion star. We detect discrete emission from innershell transitions, in addition to absorption likely due to multiple unresolved transitions in lower ionization states. The emission line intensity ratios observed in the range of the spectrum occupied by the Fe XV $n=1-2$ forbidden and intercombination lines suggest that there is a substantial contribution from resonantly scattered innershell emission from the Li- and Be-like ionization states. The Fe XV forbidden and intercombination lines arise in the ionization zone closest to the compact object, and since they are not subject to radiative transfer effects, we can use them in principle to constrain the radial velocity amplitude of the compact object. We infer that the results indicate a compact object mass of the order of the mass of the Wolf-Rayet companion star, but we note that the presence of resonantly scattered radiation from Li-like ions may complicate the interpretation of the He-like emission spectrum.

The high-mass star-forming region NGC6334I-MM1 underwent an energetic accretion event in January 2015. We report the large-scale ($10 - 100$ AU) and small-scale ($\sim 1$ AU) changes in spatial and velocity structures of 22 GHz water masers as observed with VERA before and during the accretion burst. The masers in the northern bow-shock CM2-W2 brightened, and better traced a bow structure during the burst. In the southern regions, there was both activation and disappearance of associations before and during the burst. We measured the amplitudes, central velocities and FWHMs of about 20 features in each epoch. We found that the linear scale of the brightest feature in CM2-W2 grew from 0.6 AU before the burst to 1.4 AU after the burst, possibly indicating that a larger volume of gas was able to sustain masing action as a consequence of the accretion burst. This feature also had a rapid (0.2 yr) brightness increase by a factor of four, which has been previously reported in long-term single-dish monitoring. We propose that the water maser flare could be explained by an increase of the collisional pump rate due to radiative heating of H$_2$ by increased high energy radiation (UV or X-ray) from the inner protostellar core. We also describe the spot and spectral method of maser proper motion calculations. We argue that for high spectral resolution observations the spectral method is more robust for calculating proper motions than the spot method.

We present the $z\!\approx\!6$ type-1 quasar luminosity function (QLF) based on the Pan-STARRS1 (PS1) quasar survey. The PS1 sample includes 125 quasars at $z\approx5.7-6.2$ with $-28\lesssim M_{1450}\lesssim-25$. Complemented by 48 fainter quasars from the SHELLQs survey, we evaluate the $z\approx6$ QLF over $-28\lesssim M_{1450}\lesssim-22$. Adopting a double power law with an exponential evolution of the quasar density ($\Phi(z)\propto10^{k(z-6)}$; $k=-0.7$), we use a maximum likelihood method to model our data. We find a break magnitude of $M^*=-26.38_{-0.60}^{+0.79}\,\text{mag}$, a faint end slope of $\alpha=-1.70_{-0.19}^{+0.29}$, and a steep bright end slope with $\beta=-3.84_{-1.21}^{+0.63}$. % Based on our new QLF model we determine the quasar comoving spatial density at $z\!\approx\!6$ to be $n( M_{1450}<-26)=1.16_{-0.12}^{+0.13}\,\text{cGpc}^{-3}$. In comparison with the literature, we find the quasar density to evolve with a constant value of $k\approx-0.7$ from $z\approx7$ to $z\approx4$. % Additionally, we derive an ionizing emissivity of $\epsilon_{912}(z=6) =7.23_{-1.02}^{+1.65}\times 10^{22}\,\text{erg}\,\text{s}^{-1}\text{Hz}^{-1}\text{cMpc}^{-3}$ based on the QLF measurement. Given standard assumptions and the recent measurement of the mean free path of Becker et al. (2021) at $z\approx6$ we calculate an HI photoionizing rate of $\Gamma_{\text{HI}}(z{=}6) \approx 6 \times10^{-16}\,\text{s}^{-1}$, strongly disfavoring a dominant role of quasars in hydrogen reionization.

BL Lacertae had undergone a series of historical high flux activity over a year, from August 2020 in the optical to VHE $\gamma$-rays. In this paper, we report on optical flux and spectral variability of the first historical maxima outburst event during October-November in g, r and i bands with the 1.26m telescope at Xinglong observatory, China. We detected significant intranight variations with amplitude rising up to $\sim 30$%, when the fastest variability timescale is found to be a few tens of minutes, giving an emitting region size of the order $10^{-3}$ pc, which corresponds to $\sim 100$ Schwarzschild radius of the central black hole, likely coming from some jet mini-structures. Unlike on intranight timescale, a clear frequency dependent pattern along with symmetric timescales ($\sim$ 11d) of flux variation are detected on long timescale. The spectral evolution was predominated by flattening of the spectra with increasing brightness i.e., a bluer-when-brighter trend in 96% of the cases. On the night before the outburst peak, the color indices clustered in two distinct branches in color--magnitude diagram within a period of $\sim$ 6 hours that is connected to a hard-soft-hard spectral evolution trend extracted from time-resolved spectra. Such trend has never seen in BL Lac or any other blazars before to the best of our knowledge. The results obtained in this study can be explained in the context of shock induced particle acceleration or magnetic re-connection in the jet where turbulent processes most likely resulted the asymmetric flux variation on nightly timescale.

Hydrodynamical interactions between binaries and circumbinary disks (CBDs) play an important role in a variety of astrophysical systems, from young stellar binaries to supermassive black hole binaries. Previous simulations of CBDs have mostly employed locally isothermal equation of state. We carry out two-dimensional viscous hydrodynamic simulations of CBDs around equal-mass, circular binaries, treating the gas thermodynamics by thermal relaxation towards equilibrium temperature (the constant-$\beta$ cooling ansatz, where $\beta$ is the cooling time in units of the local Keplerian time). As an initial study, we use the grid-based code Athena++ on a polar grid, covering an extended disk outside the binary co-orbital region. We find that with a longer cooling time, the accretion variability is gradually suppressed, and the morphology of the CBD becomes more symmetric. The disk also shows evidence of hysterisis behavior depending on the initial conditions. Gas cooling also affects the rate of angular momentum transfer between the binary and the CBD, where given our adopted disk thickness and viscosity ($H/r\sim 0.1$ and $\alpha\sim 0.1$), the binary orbit expands while undergoing accretion for most $\beta$ values between 0 and 4.0 except over a narrow range of intermediate $\beta$ values. The validity of using polar grid excising the central domain is also discussed.

Using data collected at the Pierre Auger Observatory, we search for signatures of instanton-induced processes that would provide evidence of super-heavy particles decaying in the Galactic halo. Such particles could have been produced sufficiently during the post-inflationary epoch to match the relic abundance of dark matter inferred today. The non-observation of these signatures allows us to probe the instanton strength and to derive a bound, the best ever obtained from instanton-mediated processes, on the reduced coupling constant of gauge interactions in the dark sector: $\alpha_X \lesssim 0.09$, for $10^{9} \lesssim M_X/{\rm GeV} < 10^{19}$.

The transport of charged particles in various astrophysical environments permeated by magnetic fields is described in terms of a diffusion process, which relies on diffusion-tensor parameters generally inferred from Monte-Carlo simulations. Based on a red-noise approximation to model the two-point correlation function of the magnetic field experienced by charged particles between two successive times, the diffusion-tensor coefficients were previously derived in the case of pure turbulence. In this contribution to ECRS2022, the derivation is extended to the case of a mean field on top of the turbulence. The results are applicable to a variety of astrophysical environments in regimes where the Larmor radius of the particles is resonant with the power spectrum of the turbulence wavelength (gyro-resonant regime), or where the Larmor radius is greater than the largest turbulence wavelength (high-rigidity regime).

Planetesimal accretion is a key source for heavy-element enrichment in giant planets. It has been suggested that Jupiter's enriched envelope is a result of planetesimal accretion during its growth assuming it formed in a massive planetesimal disk. In this study, we simulate Jupiter's formation in this scenario. We assume in-situ formation and perform N-body simulations to infer the solid accretion rate. We find that tens-Earth masses of planetesimals can be captured by proto-Jupiter during the rapid gas accretion phase. However, if several embryos are formed near Jupiter's core, which is an expected outcome in the case of a massive planetesimal disk, scattering from the embryos increases the eccentricity and inclination of planetesimals and therefore significantly reduces the accretion efficiency. We also compare our results with published semi-analytical models and show that these models cannot reproduce the N-body simulations especially when the planetesimal disk has a large eccentricity and inclination. We show that when the dynamical evolution of planetesimals is carefully modelled, the total mass of captured planetesimals $M_\mathrm{cap,tot}$ is $2 M_\oplus \lesssim M_\mathrm{cap,tot}\lesssim 18 M_\oplus$. The metallicity of Jupiter's envelope can be explained by the planetesimal accretion in our massive disk model despite the low accretion efficiency coming from the high eccentricity and inclination of planetesimals. Our study demonstrates the importance of detailed modelling of planetesimal accretion during the planetary growth and its implications to the heavy-element mass in gaseous planets.

In this study, we expand the discussion of extrasolar planetary research by proposing a new approach to detecting extrasolar moons using Kepler (K2) light curves. We shaped this investigation by comparing transit light curve data of ~5000 main-sequence stars cataloged in the NASA Exoplanet Archive, the Radio Galaxy Zoo Exoplanet Explorers and the Exoplanet Follow Up Observing Program, which served as the repository to collect and analyze supplementary K2 data. By examining K2 light curves, various characteristics related to transits were modeled, which we then compared with confirmed extrasolar planets, variable stars such as eclipsing binaries and noise or gaps in the data. For illustration, perturbations in the timing of two separate transits for 2MASS J08251369+1425306 (Rs ~ 0.346) produced two characteristic decreases in the photon flux followed by two increases, inferring the presence of an extrasolar planet (Rp ~ 0.0129) and its companion, an extrasolar moon (Rm ~ 0.0048). To test our hypothesis, we scrutinized various competing assumptions to resolve the source of the duo photon flux, namely the Rossiter McLaughlin effect, limb darkening, a multiplanet system or a planet occulting one or more sunspots on the surface of the star. However, we uncovered no analogous light curve in K2 to fit the competing hypotheses and account for the proposed extrasolar companion.

We model the early stages of planet formation in the Solar System, including continual planetesimal formation, and planetesimal and pebble accretion onto planetary embryos in an evolving disk driven by a disk wind. The aim is to constrain aspects of planet formation that have large uncertainties by matching key characteristics of the Solar System. The model produces a good fit to these characteristics for a narrow range of parameter space. Planetary growth beyond the ice line is dominated by pebble accretion. Planetesimal accretion is more important inside the ice line. Pebble accretion inside the ice line is slowed by higher temperatures, partial removal of inflowing pebbles by planetesimal formation and pebble accretion further out in the disk, and increased radial velocities due to gas advection. The terrestrial planets are prevented from accreting much water ice because embryos beyond the ice line reach the pebble isolation mass before the ice line enters the terrestrial-planet region. When only pebble accretion is considered, embryos typically remain near their initial mass or grow to the pebble-isolation mass. Adding planetesimal accretion allows Mars-sized objects to form inside the ice line, and allows giant-planet cores to form over a wider region beyond the ice line. In the region occupied by Mercury, pebble Stokes numbers are small. This delays the formation of embryos and stunts their growth, so that only low-mass planets can form here.

The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angular power spectrum of galaxy number counts and cosmic shear data without invoking the so-called 'Limber approximation', in the context of the Rubin Observatory Legacy Survey of Space and Time (LSST). We quantify the performance, in terms of accuracy and speed, of three non-Limber implementations: ${\tt FKEM (CosmoLike)}$, ${\tt Levin}$, and ${\tt matter}$, themselves based on different integration schemes and approximations. We find that in the challenge's fiducial 3x2pt LSST Year 10 scenario, ${\tt FKEM (CosmoLike)}$ produces the fastest run time within the required accuracy by a considerable margin, positioning it favourably for use in Bayesian parameter inference. This method, however, requires further development and testing to extend its use to certain analysis scenarios, particularly those involving a scale-dependent growth rate. For this and other reasons discussed herein, alternative approaches such as ${\tt matter}$ and ${\tt Levin}$ may be necessary for a full exploration of parameter space. We also find that the usual first-order Limber approximation is insufficiently accurate for LSST Year 10 3x2pt analysis on $\ell=200-1000$, whereas invoking the second-order Limber approximation on these scales (with a full non-Limber method at smaller $\ell$) does suffice.

The gravitational coupling between large- and small-scale density perturbations leads to anisotropic distortions to the local small-scale matter fluctuations. Such local anisotropic distortions can be used to reconstruct the large-scale matter distribution, known as tidal reconstruction. In this paper, we apply the tidal reconstruction methods to simulated galaxies in redshift space. We find that redshift-space distortions lead to anisotropic reconstruction results. While the reconstructed radial modes are more noisy mainly due to the small-scale velocity dispersion, the transverse modes are still reconstructed with high fidelity, well correlated with the original large-scale density modes. The bias of the reconstructed field at large scales shows a simple angular dependence which can be described by a form similar to the linear redshift-space distortion. The noise power spectrum is nearly isotropic and scale-independent on large scales. This makes the reconstructed tides field an ideal tracer for cosmic variance cancellation and multi-tracer analysis and has profound implications for future 21~cm intensity mapping surveys.

The parameters of cosmic ray electrons (CRe) and the strength of magnetic field are crucial for studying the particle acceleration of shocks. The recent discovery of eROSITA bubbles suggests that the well-known NPS/Loop I should be a 10-kpc sized relic in the Galactic halo instead of a small local structure near the Sun. By deriving the energy density of CRe and magnetic field strength accounting for the NPS, unprecedentedly precise parameter constraints on particle acceleration for weak shocks in the halo medium can be provided. The parameters of CRe and magnetic field can be derived independently by modeling the gamma-ray and the radio data of the NPS via inverse Compton scattering and synchrotron emission, respectively. Our main results are: (1), the energy density of CRe is $(2-5)\times 10^{-14}$ erg cm-3, and the spectral index is $p\simeq 2.0\pm 0.1$ below the cooling break energy of $\sim 4$ GeV; (2), the magnetic field strength is 4 $\mu$G; (3), the shock acceleration efficiency of CRe is close to 1%. Given the Mach number of 1.5, the high acceleration efficiency and flat spectrum of CRe suggest that re-acceleration of preexisting relativistic electrons should exist in the NPS. In addition, the cooling break of $\sim$ 4 GeV suggests that the cooling timescale is $10^7$ yr, in agreement with the age of the eROSITA bubbles.

We report the discovery of a third planet transiting the star TOI-1260, previously known to host two transiting sub-Neptune planets with orbital periods of 3.127 and 7.493 days, respectively. The nature of the third transiting planet with a 16.6-day orbit is supported by ground-based follow-up observations, including time-series photometry, high-angular resolution images, spectroscopy, and archival imagery. Precise photometric monitoring with CHEOPS allows to improve the constraints on the parameters of the system, improving our knowledge on their composition. The improved radii of TOI-1260b, TOI-1260c are $2.36 \pm 0.06 \rm R_{\oplus}$, $2.82 \pm 0.08 \rm R_{\oplus}$, respectively while the newly discovered third planet has a radius of $3.09 \pm 0.09 \rm R_{\oplus}$. The radius uncertainties are in the range of 3\%, allowing a precise interpretation of the interior structure of the three planets. Our planet interior composition model suggests that all three planets in the TOI-1260 system contains some fraction of gas. The innermost planet TOI-1260b has most likely lost all of its primordial hydrogen-dominated envelope. Planets c and d were also likely to have experienced significant loss of atmospheric through escape, but to a lesser extent compared to planet b.

Within the model of self-gravitating Bose--Einstein condensate (BEC) dark matter (DM) it is argued that the axion-like self-interaction of ultralight bosons provides the existence of rarefied and dense phases, which are predicted earlier on the base of the models with polynomial-like self-interactions. Associating the very short scattering length in BEC DM with the predominant participating composites of few DM particles, we attempt to form a dimer of two particles at a quantum mechanical level, using a smooth $\mu$-deformation of the axion cosine-like potential and replacing the field-dependent argument with the distance between particles. Part of the obtained results concerns potential two-particle scattering with $\mu$-deformed interaction, and they allow us to focus on a special option with unique values of the deformation parameter $\mu=1$ and the coupling constant. In this case of the potential with an infinite scattering length, we get a rather simple solution for the dimer in the ground state. We involve two-channel scattering and Feshbach resonance to describe the formation of a dimer in space. Specifying the parameters of interactions, we reveal a long-lived resonance that occurs when a pair of particles jumps between the open and closed scattering channels with close energy values. This indicates the possibility of participation of such dimers in forming BEC DM halo of galaxies.

We present a new parametric lens model for the massive galaxy cluster Abell 2744 based on the new ultra-deep James Webb Space Telescope (JWST) imaging taken in the framework of the UNCOVER program. These observations constitute the deepest JWST images of a lensing cluster to date, adding to the recent JWST ERS and DDT data taken for this field. The wide field-of-view of UNCOVER ($\sim45$ arcmin$^2$) extends beyond the cluster's well-studied central core and reveals a spectacular wealth of prominent lensed features around two massive cluster sub-structures in the north and north-west, where no multiple images were previously known. The 75 newly uncovered multiple images and candidates of 16 sources allow us, for the first time, to constrain the lensing properties and total mass distribution around these extended cluster structures using strong lensing (SL). Our model yields an effective Einstein radius of $\theta_{E,\mathrm{main}}\simeq23''$ for the main cluster core (for $z_{\mathrm{s}}=2$), enclosing a mass of $M(\theta<\theta_{E,\mathrm{main}})\simeq7.7\times10^{13}$ M$_{\odot}$, and $\theta_{E,\mathrm{NW}}\simeq13''$ for the newly discovered north-western SL structure enclosing $M(\theta<\theta_{E,\mathrm{NW}})\simeq2.2\times10^{13}$ M$_{\odot}$. The northern clump is somewhat less massive with $\theta_{E,\mathrm{N}}\simeq7''$ enclosing $M(\theta<\theta_{E,\mathrm{N}})\simeq8\times10^{12}$ M$_{\odot}$. We find the northern sub-structures of Abell 2744 to broadly agree with the findings from weak lensing (WL) and align with the filamentary structure found by these previous studies. Our model in particular reveals a large area of high magnifications between the various cluster structures, which will be paramount for lensed galaxy studies in the UNCOVER field. The model is made publicly available to accompany the first UNCOVER data release.

The Hayabusa2 mission impact experiment on asteroid Ryugu created an unexpectedly large crater. The associated regime of low-gravity, low-strength cratering remained largely unexplored so far, because these impact conditions cannot be re-created in laboratory experiments on Earth. Here we show that the target cohesion may be very low and the impact probably occurred in the transitional cratering regime, between strength and gravity. For such conditions, our numerical simulations are able to reproduce the outcome of the impact on Ryugu, including the effects of boulders originally located near the impact point. Consistent with most recent analysis of Ryugu and Bennu, cratering scaling-laws derived from our results suggest that surfaces of small asteroids must be very young. However, our results also show that the cratering efficiency can be strongly affected by the presence of a very small amount of cohesion. Consequently, the varying ages of different geological surface units on Ryugu may be due to the influence of cohesion.

We present the results of a set of radiation magnetohydrodynamic simulations of turbulent molecular clouds in which we vary the initial strength of the magnetic field within a range ($1 \lesssim \mu \lesssim 5$) consistent with observations of local giant molecular clouds (GMCs). We find that as we increase the strength of the magnetic field, star formation transitions from unimodal (the baseline case, $\mu=5$, with a single burst of star formation and Salpeter IMF) to bimodal. This effect is clearest in the most strongly magnetized GMCs ($\mu=1$): a first burst of star formation with duration, intensity and IMF comparable to the baseline case, is followed by a second star formation episode in which only low-mass stars are formed. Overall, due to the second burst of star formation, the strongly magnetized case results in a longer star formation period and a higher efficiency of star formation. The second burst is produced by gas that is not expelled by radiative feedback, instead remaining trapped in the GMC by the large-scale B-field, producing a nearly one-dimensional flow of gas along the field lines. The trapped gas has a turbulent and magnetic topology that differs from that of the first phase and strongly suppresses gas accretion onto protostellar cores, reducing their masses. We speculate that this star formation bimodality may be an important ingredient to understand the origin of multiple stellar populations observed in massive globular clusters.

Superhumps are among the abundant variable phenomena observed in the light curves of cataclysmic variables (CVs). They come in two flavours as positive and negative superhumps, distinguished by periods slightly longer or shorter, respectively, than the orbital periods of these interacting binary systems. Positive superhumps are ubiquitous in superoutbursting short period dwarf novae of the SU~UMa type but are less common in longer period systems with accretion disks in a permanent bright state such as novalike variables and most old novae. Negative superhumps do not seem not to have a preference for a particular type of CV. Here, I take advantage of the long high cadence light curves provided by TESS for huge number of stars, selecting all old novae and novalike variables with past reported superhumps for which TESS light curves are available and have not yet been analysed in previous publications in order to study their superhump behaviour. In combination with information taken from the literature the results enable to compile the most complete census of superhumps in these stars so far. As a corollary, for the eclipsing systems in the present sample of objects eclipse epochs derived from the TESS light curves and in some cases from archival light curves are listed and used to update orbital ephemeris and to discuss period changes.

We present analysis of observational data from the Swift Burst Analyser for a sample of 15 short gamma-ray bursts with extended emission (SGRBEEs) which have been processed such that error propagation from Swift's count-rate-to-flux conversion factor is applied to the flux measurements. We apply this propagation to data presented by the Burst Analyser at 0.3-10 keV and also at 15-50 keV, and identify clear differences in the morphologies of the light-curves in the different bands. In performing this analysis with data presented at both 0.3-10 keV, at 15-50 keV, and also at a combination of both bands, we highlight the impact of extrapolating data from their native bandpasses on the light-curve. We then test these data by fitting to them a magnetar-powered model for SGRBEEs, and show that while the model is consistent with the data in both bands, the model's derived physical parameters are generally very loosely constrained when this error propagation is included and are inconsistent across the two bands. In this way, we highlight the importance of the Swift data processing methodology to the details of physical model fits to SGRBEEs.

Breakdown of rotational invariance of the primordial power spectrum manifests in the statistical anisotropy of the observed Cosmic Microwave Background (CMB) radiation. Hemispherical power asymmetry in the CMB may be caused due to a dipolar modulation, indicating the presence of a preferred direction. Appropriately re-scaled local variance maps of the CMB temperature anisotropy data effectively encapsulate this dipolar pattern. As a first-of-its-kind method, we train Artificial Neural Networks (ANNs) with such local variances as input features to distinguish statistically isotropic CMB maps from dipole modulated ones. Our trained ANNs are able to predict components of the amplitude times the unit vector of the preferred direction for mixed sets of modulated and unmodulated maps, with goodness of fit ($R^2$) scores $>97\%$ for full sky, and $>96\%$ for partial sky coverage. On all observed foreground-cleaned CMB maps, the ANNs detect the dipolar modulation signal with overall consistent values of amplitudes and directions. This detection is significant at $97.21\%-99.38\%$ C.L. for all full sky maps, and at $98.34\%-100\%$ C.L. for all partial sky maps. Robustness of the signal holds across full and partial skies, various foreground cleaning methods, inpainting algorithms, instruments and all the different periods of observation for Planck and WMAP satellites. The significant and robust detection of the signal, in addition to the consistency of values of amplitude and directions, as found independent of any pre-existing methods, further mitigates the criticisms of look-elsewhere effects and a posteriori inferences for the preferred dipole direction in the CMB.

The identification of bright quasars at z>6 enables detailed studies of supermassive black holes, massive galaxies, structure formation, and the state of the intergalactic medium within the first billion years after the Big Bang. We present the spectroscopic confirmation of 55 quasars at redshifts 5.6<z<6.5 and UV magnitudes -24.5<M1450<-28.5 identified in the optical Pan-STARRS1 and near-IR VIKING surveys (48 and 7, respectively). Five of these quasars have been independently discovered in other studies. The quasar sample shows an extensive range of physical properties, including 17 objects with weak emission lines, ten broad absorption line quasars, and five with strong radio emission (radio-loud quasars). There are also a few notable sources in the sample, including a blazar candidate at z=6.23, a likely gravitationally lensed quasar at z=6.41, and a z=5.84 quasar in the outskirts of the nearby (D~3 Mpc) spiral galaxy M81. The blazar candidate remains undetected in NOEMA observations of the [CII] and underlying emission, implying a star-formation rate <30-70 Msun/yr. A significant fraction of the quasars presented here lies at the foundation of the first measurement of the z~6 quasar luminosity function from Pan-STARRS1 (introduced in a companion paper). The quasars presented here will enable further studies of the high-redshift quasar population with current and future facilities.

4FGL J1015.5-6030 is an unidentified Fermi-LAT source hosting a bright, extended X-ray source whose X-ray spectrum is consistent with that of a young pulsar, yet no pulsations have been found. Here we report on XMM-Newton timing and Chandra imaging observations of the X-ray counterpart of 4FGL J1015.5-6030. We find no significant periodicity from the source and place a 3$\sigma$ upper-limit on its pulsed fraction of 34$\%$. The Chandra observations resolve the point source from the extended emission. We find that the point source's spectrum is well fit by a blackbody model, with temperature $kT=0.205\pm0.009$ keV, plus a weak power-law component, which is consistent with a thermally emitting neutron star with a magnetospheric component. The extended emission spans angular scales of a few arcseconds up to about 30$''$ from the point source and its spectrum is well fit by a power-law model with a photon index $\Gamma=1.70\pm0.05$. The extended emission's spectrum and 0.5-10 keV luminosity of 4$\times10^{32}$ erg s$^{-1}$ (at a plausible distance of 2 kpc) are consistent with that of a pulsar wind nebula. Based on a comparison to other GeV and X-ray pulsars, we find that this putative pulsar is likely a middle-aged (i.e., $\tau\sim 0.1$--1 Myr) radio-quiet pulsar with $\dot{E}\sim10^{34}-10^{35}$ erg s$^{-1}$.

The Cosmic Microwave Background (CMB) primordial B-modes signal is predicted to be much lower than the polarized Galactic emission (foregrounds) in any region of the sky pointing to the need for sophisticated component separation methods. Among them, the blind Needlet-ILC (NILC) has great relevance given our current poor knowledge of the B-modes foregrounds. However the expected level of spatial variability of the foreground spectral properties complicates the NILC subtraction of the Galactic contamination. In order to reach the ambitious targets of future CMB experiments, we therefore propose a novel extension of the NILC approach, the Multi-Clustering NILC (MC-NILC), which performs NILC variance minimization on separate regions of the sky (clusters) properly chosen to have similar spectral properties of the B-modes foregrounds emission. Clusters are identified thresholding the ratio of B-modes maps at two separate frequencies which is used as tracer of the spatial distribution of the spectral indices of the Galactic emission in B modes. We consider ratios either of simulated foregrounds-only B modes (ideal case) or of cleaned templates of Galactic emission obtained from realistic simulations. In this work we present an application of MC-NILC to the future LiteBIRD satellite, which targets the observation of both reionization and recombination peaks of the primordial B-modes angular power spectrum with a total error on the tensor-to-scalar ratio $\delta r < 0.001$. We show that MC-NILC provides a CMB solution with residual foregrounds and noise contamination that is significantly reduced with respect to NILC and lower than the primordial signal targeted by LiteBIRD at all angular scales for the ideal case and at the reionization peak for a realistic ratio. Thus, MC-NILC will represent a powerful method to mitigate B-modes foregrounds for future CMB polarization experiments.

Cosmic hydrogen reionization and cosmic production of first metals are major phase transitions of the Universe occurring during the first billion years after the Big Bang, but still poorly explored observationally. Using the JWST NIRSpec prism spectroscopy, we report the discovery of a sub-$L_{\ast}$ galaxy at $z_{\rm spec}=8.1623_{-0.0008}^{+0.0007}$, dubbed RXJ2129-z8HeII, via the detection of a series of strong rest-frame UV/optical nebular emission lines and the clear Lyman break. A strong He II $\lambda$1640 emission is present, the highest redshift He II line currently known. Its high rest-frame equivalent width (EW $=19.4\pm3.2$ Angstrom) and extreme flux ratios with respect to UV metal lines and Balmer lines raise the possibility that part of RXJ2129-z8HeII's stellar populations could be Pop III-like. RXJ2129-z8HeII also shows a pronounced UV continuum with an extremely steep (i.e. blue) spectral slope of $\beta=-2.50\pm0.08$, the steepest amongst all spectroscopically confirmed galaxies at $z\gtrsim7$, in support of its very hard ionizing spectrum that could lead to a significant leakage of its ionizing flux. Therefore, RXJ2129-z8HeII is representative of the key galaxy population driving the cosmic reionization. To date, this is also the most compelling case where trace Pop III stars might coexist with more metal-enriched stars.

Surveys with James Webb Space Telescope (JWST) have discovered candidate galaxies in the first 400 Myr of cosmic time. The properties of these distant galaxies provide initial conditions for understanding early galaxy formation and cosmic reionisation. Preliminary indications have suggested these candidate galaxies may be more massive and abundant than previously thought. However, without spectroscopic confirmation of their distances to constrain their intrinsic brightnesses, their inferred properties remain uncertain. Here we report on four galaxies located in the JWST Advanced Deep Extragalactic Survey (JADES) Near-Infrared Camera (NIRCam) imaging with photometric redshifts $z\sim10-13$ subsequently confirmed by JADES JWST Near- Infrared Spectrograph (NIRSpec) observations. These galaxies include the first redshift $z>12$ systems both discovered and spectroscopically confirmed by JWST. Using stellar population modelling, we find the galaxies typically contain a hundred million solar masses in stars, in stellar populations that are less than one hundred million years old. The moderate star formation rates and compact sizes suggest elevated star formation rate surface densities, a key indicator of their formation pathways. Taken together, these measurements show that the first galaxies contributing to cosmic reionisation formed rapidly and with intense internal radiation fields.

Third-order weak lensing statistics are a promising tool for cosmological analyses since they extract cosmological information in the non-Gaussianity of the cosmic large-scale structure. However, such analyses require precise and accurate models for the covariance. In this second paper of a series on third-order weak lensing statistics, we derive and validate an analytic model for the covariance of the third-order aperture statistics $\langle M_\mathrm{ap}^3\rangle$. We derive the covariance model from a real-space estimator for $\langle M_\mathrm{ap}^3\rangle$. We validate the model by comparing it to estimates from simulated Gaussian random fields (GRF) and two sets of N-body simulations. Finally, we perform mock cosmological analyses with the model covariance and the simulation estimate to compare the resulting parameter constraints. We find good agreement between the model and the simulations, both for the GRF and the $N$-body simulations. The figure-of-merit in the $S_8$-$\Omega_\mathrm{m}$ plane from our covariance model is within 3\% of the one obtained from the simulated covariances. We also show that our model, which is based on an estimator using convergence maps, can be used to obtain upper and lower bounds for the covariance of an estimator based on three-point shear correlation functions. This second estimator is required for realistic survey data. In our derivation, we find that the $\langle M_\mathrm{ap}^3\rangle$ covariance cannot be obtained from the bispectrum covariance and that it includes several `finite-field terms' that do not scale with the inverse survey area. Our covariance model is sufficiently accurate for analysing stage III surveys. Covariances for statistics in Fourier space cannot always be straightforwardly converted into covariance for real-space statistics. The modelling code is available at https://github.com/sheydenreich/threepoint/releases/ .

We develop a novel general formalism allowing us to obtain values of the Hubble constant in agreement with late-time observables without degrading the fit to Cosmic Microwave Background data, considering perturbative modifications around a fiducial $\Lambda$CDM cosmology. Taking as proof-of-principle the case of a time-varying electron mass and fine structure constant, we demonstrate that a modified recombination can solve the Hubble tension and lower $S_8$ to match weak lensing measurements. Once baryonic acoustic oscillation and uncalibrated supernovae data are included, however, it is not possible to fully solve the tension with perturbative modifications to recombination.

Many cosmological phenomena lead to the production of primordial black holes in the early Universe. These phenomena often create a population of black holes with extended mass and spin distributions. As these black holes evaporate via Hawking radiation, they can modify various cosmological observables, lead to the production of dark matter, modify the number of effective relativistic degrees of freedom, $N_{\rm eff}$, source a stochastic gravitational wave background and alter the dynamics of baryogenesis. We consider the Hawking evaporation of primordial black holes that feature non-trivial mass and spin distributions in the early Universe. We demonstrate that the shape of such a distribution can strongly affect most of the aforementioned cosmological observables. We outline the numerical machinery we use to undertake this task. We also release a new version of FRISBHEE that handles the evaporation of primordial black holes with an arbitrary mass and spin distribution throughout cosmic history.

We propose a novel dark matter detection method utilizing the excitation of superconducting transmon qubits. Assuming the hidden photon dark matter of a mass of $O(10)\ \mu{\rm eV}$, the classical wave-matter oscillation induces an effective ac electric field via the small kinetic mixing with the ordinary photon. This serves as a coherent drive field for a qubit when it is resonant, evolving it from the ground state towards the first-excited state. We evaluate the rate of such evolution and observable excitations in the measurements, as well as the search sensitivity to the hidden photon dark matter. For a selected mass, one can reach $\epsilon \sim 10^{-12}-10^{-14}$ (where $\epsilon$ is the kinetic mixing parameter of the hidden photon) with a single standard transmon qubit. A simple extension to the frequency-tunable SQUID-based transmon enables the mass scan to cover the whole $4-40\ \mu{\rm eV}$ ($1-10$ GHz) range within a reasonable length of run time. The sensitivity scalability along the number of the qubits also makes it a promising platform in accord to the rapid evolution of the superconducting quantum computer technology.

We study the imprint of light scalar fields on gravitational waves from extreme mass ratio inspirals -- binary systems with a very large mass asymmetry. We first show that, to leading order in the mass ratio, any effects of the scalar on the waveform are captured fully by two parameters: the mass of the scalar and the scalar charge of the secondary compact object. We then use this theory-agnostic framework to show that the future observations by LISA will be able to simultaneously measure both of these parameters with enough accuracy to detect ultra-light scalars.

One of the prominent decay modes of heavy nuclei which are produced in astrophysical environments at temperatures of the order of $10^9$ K is the $\alpha$ ($^4$He) decay. Thermally enhanced $\alpha$ decay rates are evaluated within the standard scheme of a tunneling decay where the $\alpha$ particle tunnels through the potential barrier formed by its interaction with the daughter nucleus. Following the observation that there exist several excited levels with the possibility of an $\alpha$ decay when the daughter nucleus is at a shell closure, we focus in particular on decays producing daughter nuclei with the neutron number, N = 126. Within a statistical approach we find that the half-lives, $t_{1/2}(T)$, for temperatures ranging from $T$ = 0 to 2.4 GK can decrease by 1 - 2 orders of magnitude with the exception of the decay of $^{212}$Po which decays to the doubly magic daughter $^{208}$Pb, where $t_{1/2}(T)$ decreases by 5 orders of magnitude. The effect of these thermally enhanced $\alpha$ decays on the $r$-process nucleosynthesis can be significant in view of the mass build up at the waiting point nuclei with closed neutron shells.

Detecting gravitational waves with a frequency higher than 10 kHz requires a new idea. In previous papers, we proposed magnon gravitational wave detectors and gave the first limit on GHz gravitational waves by reinterpreting the existing data of axion dark matter experiments. In this paper, we show that the sensitivity can be improved by constructing the detector specific to gravitational waves. In particular, we employ an infinite sum of terms in the expansion of Fermi normal coordinates to probe gravitational waves with a wavelength comparable to the detector size. As a consequence, we obtain the sensitivity around $h_c \sim 10^{-20}$.

Recently a series of analyses on the flight time of cosmic photons and neutrinos suggests that the speed of light \emph{in vacuo} takes the energy-dependent form $v(E)\simeq 1-E/E_{\text{LIV}}^{\gamma }$ with $E_{\text{LIV}}^{\gamma }\approx 3.6\times 10^{17}~\text{GeV}$, and meanwhile the speed of neutrinos is proposed to be $v(E)\simeq 1\pm E/E_{\text{LIV}}^{\nu }$ with $E_{\text{LIV}}^{\nu }\approx 6.5\times 10^{17}~\text{GeV}$ and $\pm {}$ representing the helicity dependence. This novel picture immediately urges us to provide a satisfactory theoretical explanation. Among all the attempts to predict the speed variations from quantum gravity, we find that loop quantum gravity can serve as a good candidate for explaining the aforementioned picture consistently.

The (re)introduction of $\Lambda$ into cosmology has spurred debates that touch on central questions in philosophy of science, as well as the foundations of general relativity and particle physics. We provide a systematic assessment of the often implicit philosophical assumptions guiding the methodology of precision cosmology in relation to dark energy. We start by briefly introducing a recent account of scientific progress in terms of risky and constrained lines of inquiry. This allows us to contrast aspects of $\Lambda$ that make it relevantly different from other theoretical entities in science, such as its remoteness from direct observation or manipulability. We lay out a classification for possible ways to explain apparent accelerated expansion but conclude that these conceptually clear distinctions may blur heavily in practice. Finally, we consider the important role played in cosmology by critical tests of background assumptions, approximation techniques, and core principles, arguing that the weak anthropic principle fits into this category. We argue that some core typicality assumptions -- like the Copernican principle and the cosmological principle -- are necessary though not provable, while others -- like the strong anthropic principle and appeals to naturalness or probability in the multiverse -- are not similarly justifiable.

One of the recent attempts to address the Hubble and $\sigma_8$ tensions is to consider the Universe started out not as a de Sitter-like spacetime, but rather anti-de Sitter-like. That is, the Universe underwent an "AdS-to-dS" transition at some point. We study the possibility that there are two dark energy fluids, one of which gave rise to the anti-de Sitter-like early Universe. The interaction is modeled by the Lotka-Volterra equations, commonly used in population biology. We consider "competition" models that are further classified as "unfair competition" and "fair competition". The former involves a quintessence in competition with a phantom, and the second involves two phantom fluids. Surprisingly, even in the latter scenario it is possible for the overall dark energy to cross the phantom divide. The latter model also allows a constant $w$ "AdS-to-dS" transition, thus serving as a counter-example to the claim that such a dark energy must possess a singular equation of state. We also consider a "conversion" model in which a phantom fluid still manages to achieve "AdS-to-dS" transition even if it is being converted into a negative energy density quintessence. In these models, the energy density of the late time effective dark energy is related to the coefficient of the quadratic self-interaction term of the fluids, which is analogous to the resource capacity in population biology.

Warm natural inflation is studied for the case of the original cosine potential. The radiation bath during inflation induces a dissipation (friction) rate in the equation of motion for the inflaton field, which can potentially reduce the field excursion needed for an observationally viable period of inflation. We examine if the dissipation thus provides a mechanism to avoid the large decay constant $f \gtrsim M_{\mathrm{pl}}$ of cold cosine natural inflation. Whereas temperature independent dissipation has previously been shown to alleviate the need for a trans-Planckian decay constant $f$, we illustrate here the difficulties of accommodating a significantly sub-Planckian decay constant ($f<10^{-1}M_{\mathrm{pl}}$) in the case of the following temperature dependent dissipation rates, $\Gamma \propto T^c$, with $c=\{1,3\}$. Such dissipation rates represent physically well-motivated constructions in the literature. For each model, we map its location in the $r$-$n_s$ plane and compare with Cosmic Microwave Background data. For $c=1 \, (c=3)$, we find that agreement with CMB data requires that dissipation be in the weak (moderate) regime and that the minimum allowed value of the decay constant in the potential is $f_{\rm min} = 0.3 \, (0.8)\,M_{\mathrm{pl}}$ respectively.