Papers for Monday, Oct 10 2022

Many gravitational wave (GW) sources in the LISA band are expected to have non-negligible eccentricity. Furthermore, many of them can undergo acceleration because they reside in the presence of a tertiary. Here we develop analytical and numerical methods to quantify how the compact binary's eccentricity enhances the detection of its peculiar acceleration. We show that the general relativistic precession pattern can disentangle the binary's acceleration-induced frequency shift from the chirp-mass-induced frequency shift in GW template fitting, thus relaxing the signal-to-noise ratio requirement for distinguishing the acceleration by a factor of $10\sim100$. Moreover, by adopting the GW templates of the accelerating eccentric compact binaries, we can enhance the acceleration measurement accuracy by a factor of $\sim100$, compared to the zero-eccentricity case, and detect the source's acceleration even if it does not change during the observational time. For example, a stellar-mass binary black hole (BBH) with moderate eccentricity in the LISA band yields an error of the acceleration measurement $\sim10^{-7}m\cdot s^{-2}$ for $\rm{SNR}=20$ and observational time of $4$ yrs. In this example, we can measure the BBHs' peculiar acceleration even when it is $\sim1\rm pc$ away from a $4\times 10^{6}\rm M_{\odot}$ SMBH. Our results highlight the importance of eccentricity to the LISA-band sources and show the necessity of developing GW templates for accelerating eccentric compact binaries.

Primordial black holes (PBHs) may form from the collapse of matter overdensities shortly after the Big Bang. One may identify their existence by observing gravitational wave (GW) emissions from merging PBH binaries at high redshifts $z\gtrsim 30$, where astrophysical binary black holes (BBHs) are unlikely to merge. The next-generation ground-based GW detectors, Cosmic Explorer and Einstein Telescope, will be able to observe BBHs with total masses of $\mathcal{O}(10-100)~M_{\odot}$ at such redshifts. This paper serves as a companion paper of arXiv:2108.07276, focusing on the effect of higher-order modes (HoMs) in the waveform modeling, which may be detectable for these high redshift BBHs, on the estimation of source parameters. We perform Bayesian parameter estimation to obtain the measurement uncertainties with and without HoM modeling in the waveform for sources with different total masses, mass ratios, orbital inclinations and redshifts observed by a network of next-generation GW detectors. We show that including HoMs in the waveform model reduces the uncertainties of redshifts and masses by up to a factor of two, depending on the exact source parameters. We then discuss the implications for identifying PBHs with the improved single-event measurements, and expand the investigation of the model dependence of the relative abundance between the BBH mergers originating from the first stars and the primordial BBH mergers as shown in arXiv:2108.07276.

Despite their vital importance to understand galaxy evolution and our own Galactic ecosystem, our knowledge of the physical properties of the hot phase of the Milky Way is still inadequate. However, sensitive SRG/eROSITA large area surveys are now providing us with the long sought-after data needed to mend this state of affairs. We present the properties of the soft X-ray emission as observed by eROSITA in the eFEDS field. We measure the temperature and metal abundance of the hot circum-Galactic medium (CGM) to be within $kT_{CGM}=0.153-0.178$ keV and $Z_{CGM}=0.052-0.072$ $Z_\odot$, depending on the contribution of solar wind charge exchange (SWCX). Slightly larger CGM abundances $Z_{CGM}=0.05-0.10$ $Z_\odot$ are possible, considering the uncertain extrapolation of the extragalactic Cosmic X-ray background (CXB) emission below $\sim1$ keV. To recover CGM abundances as large as $Z_{CGM}=0.3$ $Z_\odot$, it must be postulated the presence of an additional component, likely associated with the warm-hot intergalactic medium, providing $\sim15-20$% of the flux in the soft X-ray band. The emission in the soft band is dominated by the CGM, with contributions from the CXB and the local hot bubble. Moreover, the eROSITA data require the presence of an additional component associated with the elusive Galactic corona plus a possible contribution from unresolved M dwarf stars. This component has a temperature of $kT\sim0.4-0.7$ keV and it might be out of thermal equilibrium. It contributes $\sim9$% to the total emission in the 0.6--2 keV band, therefore it is a likely candidate to produce part of the unresolved CXB flux observed in X-ray ultra-deep fields. We also observe a significant contribution to the soft X-ray flux due to SWCX, during periods characterised by stronger solar wind activity, and causing the largest uncertainty on the determination of the CGM temperature.

Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts revealed that we may be missing a large amount of star formation taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically-dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag), combined with a de-blending technique, is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identify 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 are neglected by previous LBGs and H-dropout selection techniques. After performing stacking analysis, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift.

Non-Gaussian likelihoods, ubiquitous throughout cosmology, are a direct consequence of nonlinearities in the physical model. Their treatment requires Monte-Carlo Markov-chain or more advanced sampling methods for the determination of confidence contours. As an alternative, we construct canonical partition functions as Laplace-transforms of the Bayesian evidence, from which MCMC-methods would sample microstates. Cumulants of order $n$ of the posterior distribution follow by direct $n$-fold differentiation of the logarithmic partition function, recovering the classic Fisher-matrix formalism at second order. We connect this approach for weakly non-Gaussianities to the DALI- and Gram-Charlier expansions and demonstrate the validity with a supernova-likelihood on the cosmological parameters $\Omega_m$ and $w$. We comment on extensions of the canonical partition function to include kinetic energies in order to bridge to Hamilton Monte-Carlo sampling, and on ensemble Markov-chain methods, as they would result from transitioning to macrocanonical partition functions depending on a chemical potential. Lastly we demonstrate the relationship of the partition function approach to the Cram\'er-Rao boundary and to information entropies.

Indirect detection opens a unique window for probing thermal dark matter (DM): the same annihilation process that determined the relic abundance in the early Universe drives the present day astrophysical signal. While TeV-scale particles weakly coupled to the Standard Model face undoubted challenges from decades of null searches, the scenario remains compelling, and simple realizations such as Higgsino DM remain largely unexplored. The fate of such scenarios could be determined by gamma-ray observations of the centre of the Milky Way with Imaging Atmospheric Cherenkov Telescopes (IACTs). We consider the ultimate sensitivity of current IACTs to a broad range of TeV-scale DM candidates - including specific ones such as the Wino, Higgsino, and Quintuplet. To do so, we use realistic mock H.E.S.S.-like observations of the inner Milky Way halo, and provide a careful assessment of the impact of recent Milky Way mass modeling, instrumental and astrophysical background uncertainties in the Galactic Center region, and the theoretical uncertainty on the predicted signal. We find that the dominant systematic for IACT searches in the inner Galaxy is the unknown distribution of DM in that region, however, beyond this the searches are currently statistically dominated indicating a continued benefit from more observations. For two-body final states at $1~{\rm TeV}$, we find a H.E.S.S.-like observatory is sensitive to $\langle \sigma v \rangle \sim 3 \times 10^{-26}-4 \times 10^{-25}~{\rm cm}^3{\rm s}^{-1}$, except for neutrino final states, although we find results competitive with ANTARES. In addition, the thermal masses for the Wino and Quintuplet can be probed; the Higgsino continues to be out of reach by at least a factor of a few. Our conclusions are also directly relevant to the next generation Cherenkov Telescope Array, which remains well positioned to be the discovery instrument for thermal DM.

We present a publicly available software package developed for exploring apodized pupil Lyot coronagraph (APLC) solutions for various telescope architectures. In particular, the package optimizes the apodizer component of the APLC for a given focal-plane mask and Lyot stop geometry to meet a set of constraints (contrast, bandwidth etc.) on the coronagraph intensity in a given focal-plane region (i.e. dark zone). The package combines a high-contrast imaging simulation package HCIPy with a third-party mathematical optimizer (Gurobi) to compute the linearly optimized binary mask that maximizes transmission. We provide examples of the application of this toolkit to several different telescope geometries, including the Gemini Planet Imager (GPI) and the High-contrast imager for Complex Aperture Telescopes (HiCAT) testbed. Finally, we summarize the results of a preliminary design survey for the case of a 6~m aperture off-axis space telescope, as recommended by the 2020 NASA Decadal Survey, exploring APLC solutions for different segment sizes. We then use the Pair-based Analytical model for Segmented Telescope Imaging from Space (PASTIS) to perform a segmented wavefront error tolerancing analysis on these solutions.

We use the 30 simulations of the Auriga Project to estimate the temporal dependency of the inflow, outflow and net accretion rates onto the discs of Milky Way-like galaxies. The net accretion rates are found to be similar for all galaxies at early times, increasing rapidly up to $\sim 10~\mathrm{M}_\odot \, \mathrm{yr}^{-1}$. After $\sim 6~\mathrm{Gyr}$ of evolution, however, the net accretion rates are diverse: in most galaxies, these exhibit an exponential-like decay, but some systems instead present increasing or approximately constant levels up to the present time. An exponential fit to the net accretion rates averaged over the MW analogues yields typical decay time-scale of $7.2~\mathrm{Gyr}$. The analysis of the time-evolution of the inflow and outflow rates, and their relation to the star formation rate (SFR) in the discs, confirms the close connection between these quantities. First, the inflow$/$outflow ratio stays approximately constant, with typical values of $\dot{M}_\mathrm{out}/ \dot{M}_\mathrm{in} \sim 0.75$, indicating that the gas mass involved in outflows is of the order of 25% lower compared to that involved in inflows. A similar behaviour is found for the SFR$/$inflow rate ratio, with typical values between 0.1 and 0.3, and for the outflow rate$/$SFR which varies in the range $3.5$--$5.5$. Our results show that continuous inflow is key to the SFR levels in disc galaxies, and that the star formation activity and the subsequent feedback in the discs is able to produce mass-loaded galaxy winds in the disc-halo interface.

The Q&U Bolometric Interferometer for Cosmology (QUBIC) is a novel kind of polarimeter optimized for the measurement of the B-mode polarization of the Cosmic Microwave Back-ground (CMB), which is one of the major challenges of observational cosmology. The signal is expected to be of the order of a few tens of nK, prone to instrumental systematic effects and polluted by various astrophysical foregrounds which can only be controlled through multichroic observations. QUBIC is designed to address these observational issues with a novel approach that combines the advantages of interferometry in terms of control of instrumental systematics with those of bolometric detectors in terms of wide-band, background-limited sensitivity.

The composition the atmosphere of Venus results from the integration of many processes entering into play over the entire geological history of the planet. Determining the elemental abundances and isotopic ratios of noble gases (He, Ne, Ar, Kr, Xe) and stable isotopes (H, C, N, O, S) in the Venus atmosphere is a high priority scientific target since it could open a window on the origin and early evolution of the entire planet. This chapter provides an overview of the existing dataset on noble gases and stable isotopes in the Venus atmosphere. The current state of knowledge on the origin and early and long-term evolution of the Venus atmosphere deduced from this dataset is summarized. A list of persistent and new unsolved scientific questions stemming from recent studies of planetary atmospheres (Venus, Earth and Mars) are described. Important mission requirements pertaining to the measurement of volatile elements in the atmosphere of Venus as well as potential technical difficulties are outlined.

I derive the overall best distances for all 402 known galactic novae, and I collect their many properties. The centrepiece is the 74 novae with accurate parallaxes from the new Gaia data release. For the needed priors, I have collected 171 distances based on old methods (including expansion parallaxes and extinction distances). Further, I have collected the V-magnitudes at peak and the extinction measures, so as to produce absolute magnitudes at peak and then derive a crude distance as a prior. Further, I have recognized that 41 per cent of the known novae are concentrated in the bulge, with 68 per cent of these <5.4 degrees from the galactic centre, so the 165 bulge novae must have distances of 8000+-750 parsecs. Putting this all together, I have derived distances to all 402 novae, of which 220 have distances to an accuracy of better than 30 per cent. I find that the disc novae have an exponential scale height of 140+-10 pc. The average peak absolute V-magnitude is -7.45, with an RMS scatter of 1.33 mag. These peak luminosities are significantly correlated with the decline rate (t_3 in days) as M_V,peak = -7.6 + 1.5 * Log[t_3 / 30]. The huge scatter about this relation masks the correlation in many smaller datasets, and makes this relation useless for physical models. The bulge novae are indistinguishable from the disc novae in all properties, except that the novae with red giant companion stars have a strong preference for residing in the bulge population.

We present measurements of the polarization of X-rays in the 2-8 keV band from the pulsar in the ultracompact low mass X-ray binary 4U1626-67 using data from the Imaging X-ray Polarimetry Explorer (IXPE). The 7.66 s pulsations were clearly detected throughout the IXPE observations as well as in the NICER soft X-ray observations, which we use as the basis for our timing analysis and to constrain the spectral shape over 0.4-10 keV energy band. Chandra HETGS high-resolution X-ray spectra were also obtained near the times of the IXPE observations for firm spectral modeling. We find an upper limit on the pulse-averaged linear polarization of <4% (at 95% confidence). Similarly, there was no significant detection of polarized flux in pulse phase intervals when subdividing the bandpass by energy. However, spectropolarimetric modeling over the full bandpass in pulse phase intervals provide a marginal detection of polarization of the power-law spectral component at the 4.8 +/- 2.3% level (90% confidence). We discuss the implications concerning the accretion geometry onto the pulsar, favoring two-component models of the pulsed emission.

There is untapped cosmological information in galaxy redshift surveys in the non-linear regime. In this work, we use the AEMULUS suite of cosmological $N$-body simulations to construct Gaussian process emulators of galaxy clustering statistics at small scales ($0.1-50 \: h^{-1}\,\mathrm{Mpc}$) in order to constrain cosmological and galaxy bias parameters. In addition to standard statistics -- the projected correlation function $w_\mathrm{p}(r_\mathrm{p})$, the redshift-space monopole of the correlation function $\xi_0(s)$, and the quadrupole $\xi_2(s)$ -- we emulate statistics that include information about the local environment, namely the underdensity probability function $P_\mathrm{U}(s)$ and the density-marked correlation function $M(s)$. This extends the model of AEMULUS III for redshift-space distortions by including new statistics sensitive to galaxy assembly bias. In recovery tests, we find that the beyond-standard statistics significantly increase the constraining power on cosmological parameters of interest: including $P_\mathrm{U}(s)$ and $M(s)$ improves the precision of our constraints on $\sigma_8$ by 33%, $\Omega_m$ by 28%, and the growth of structure parameter, $f \sigma_8$, by 18% compared to standard statistics. We additionally find that scales below $4 \: h^{-1}\,\mathrm{Mpc}$ contain as much information as larger scales. The density-sensitive statistics also contribute to constraining halo occupation distribution parameters and a flexible environment-dependent assembly bias model, which is important for extracting the small-scale cosmological information as well as understanding the galaxy-halo connection. This analysis demonstrates the potential of emulating beyond-standard clustering statistics at small scales to constrain the growth of structure as a test of cosmic acceleration. Our emulator is publicly available at https://github.com/kstoreyf/aemulator.

We present multiwavelength observations of supernova (SN) 2017hcc with the Chandra X-ray telescope and the X-ray telescope onboard Swift (Swift-XRT) in X-ray bands, with the Spitzer and the TripleSpec spectrometer in near-infrared (IR) and mid-IR bands and with the Karl G. Jansky Very Large Array (VLA) for radio bands. The X-ray observations cover a period of 29 to 1310 days, with the first X-ray detection on day 727 with the Chandra. The SN was subsequently detected in the VLA radio bands from day 1000 onwards. While the radio data are sparse, synchrotron-self absorption is clearly ruled out as the radio absorption mechanism. The near- and the mid-IR observations showed that late time IR emission dominates the spectral energy distribution. The early properties of \snhcc\ are consistent with shock breakout into a dense mass-loss region, with $\dot M \sim 0.1 M_\odot \rm yr^{-1}$ for a decade. At few 100 days, the mass-loss rate declined to $\sim 0.02 M_\odot \rm yr^{-1}$, as determined from the dominant IR luminosity. In addition, radio data also allowed us to calculate a mass-loss rate at around day 1000, which is two orders of magnitude smaller than the mass-loss rate estimates around the bolometric peak. These values indicate that the SN progenitor underwent an enhanced mass-loss event a decade before the explosion. The high ratio of IR to X-ray luminosity is not expected in simple models and is possible evidence for an asymmetric circumstellar region.

We present novel techniques and methodology for unresolved photometric characterization of low-Earth Orbit (LEO) satellites. With the Pomenis LEO Satellite Photometric Survey our team has made over 14,000 observations of Starlink and OneWeb satellites to measure their apparent brightness. From the apparent brightness of each satellite, we calculate a new metric: the effective albedo, which quantifies the specularity of the reflecting satellite. Unlike stellar magnitude units, the effective albedo accounts for apparent range and phase angle and enables direct comparison of different satellites. Mapping the effective albedo from multiple observations across the sky produces an all-sky photometric signature which is distinct for each population of satellites, including the various sub-models of Starlink satellites. Space Situational Awareness (SSA) practitioners can use all-sky photometric signatures to differentiate populations of satellites, compare their reflection characteristics, identify unknown satellites, and find anomalous members. To test the efficacy of all-sky signatures for satellite identification, we applied a machine learning classifier algorithm which correctly identified the majority of satellites based solely on the effective albedo metric and with as few as one observation per individual satellite. Our new method of LEO satellite photometric characterization requires no prior knowledge of the satellite's properties and is readily scalable to large numbers of satellites such as those expected with developing communications mega-constellations.

The classical globular clusters found in all galaxy types have half-light radii of $r_{\rm h} \sim$ 2-4 pc, which have been tied to formation in the dense cores of giant molecular clouds. Some old star clusters have larger sizes, and it is unclear if these represent a fundamentally different mode of low-density star cluster formation. We report the discovery of a rare, young "faint fuzzy" star cluster, NGC 247-SC1, on the outskirts of the low-mass spiral galaxy NGC 247 in the nearby Sculptor group, and measure its radial velocity using Keck spectroscopy. We use Hubble Space Telescope imaging to measure the cluster half-light radius of $r_{\rm h} \simeq 12$ pc and a luminosity of $L_V \simeq 4\times10^5 \mathrm{L}_\odot$. We produce a colour-magnitude diagram of cluster stars and compare to theoretical isochrones, finding an age of $\simeq$ 300 Myr, a metallicity of [$Z$/H] $\sim -0.6$ and an inferred mass of $M_\star \simeq 9\times10^4 \mathrm{M}_\odot$. The narrow width of blue-loop star magnitudes implies an age spread of $\lesssim$ 50 Myr, while no old red-giant branch stars are found, so SC1 is consistent with hosting a single stellar population, modulo several unexplained bright "red straggler" stars. SC1 appears to be surrounded by tidal debris, at the end of a $\sim$ 2 kpc long stellar filament that also hosts two low-mass, low-density clusters of a similar age. We explore a link between the formation of these unusual clusters and an external perturbation of their host galaxy, illuminating a possible channel by which some clusters are born with large sizes.

When the extent of protoplanetary melting approached magma ocean (MO)-like conditions, alloy melts efficiently segregated from the silicates to form metallic cores. The nature of the MO of a differentiating protoplanet, i.e., internal or external MO (IMO or EMO), not only determines the abundances of life-essential volatiles like nitrogen (N) and carbon (C) in its core and mantle reservoirs but also the timing and mechanism of volatile loss. Whether the earliest formed protoplanets had IMOs or EMOs is, however, poorly understood. Here we model equilibrium N and C partitioning between alloy and silicate melts in the absence (IMO) or presence (EMO) of vapor degassed atmospheres. Bulk N and C inventories of the protoplanets during core formation are constrained for IMOs and EMOs by comparing the predicted N and C abundances in the alloy melts from both scenarios with N and C concentrations in the parent cores of magmatic iron meteorites. Our results show that in comparison to EMOs, protoplanets having IMOs satisfy N and C contents of the parent cores with substantially lower amounts of bulk N and C present in the parent body during core formation. As the required bulk N and C contents for IMOs and EMOs are in the sub-chondritic and chondritic range, respectively, N and C fractionation models alone cannot be used to distinguish the prevalence of these two end-member differentiation regimes. A comparison of N and C abundances in chondrites with their peak metamorphic temperatures suggests that protoplanetary interiors could lose a substantial portion of their N and C inventories with increasing degrees of thermal metamorphism.

The remnants of binary black hole mergers can be given recoil kick velocities up to $\unit[5,000]{km\,s^{-1}}$ due to anisotropic emission of gravitational waves. E1821+643 is a recoiling supermassive black hole moving at $\sim \unit[2,100]{km\,s^{-1}}$ along the line-of-sight relative to its host galaxy. This suggests a recoil kick of $\sim \unit[2,240]{km\,s^{-1}}$. Such a kick is powerful enough to eject E1821+643 from its $M_\text{gal}\sim2 \times \unit[10^{12}]{M_\odot}$ host galaxy. In this work, we address the question: what are the likely properties of the progenitor binary that formed E1821+643? Using astrophysically motivated priors, we infer that E1821+643 was likely formed from a binary black hole system with masses of $m_1\sim 1.9^{+5.0}_{-3.8}\times \unit[10^9]{M_\odot}$, $m_2\sim 8.1^{+3.9}_{-3.2} \times \unit[10^8]{M_\odot}$ (90\% credible intervals). The black holes in this binary were likely to be spinning rapidly with dimensionless spin magnitudes of ${\chi}_1 = 0.87^{+0.11}_{-0.26}$, ${\chi}_2 = 0.77^{+0.19}_{-0.37}$. Such a high recoil velocity is impossible for spins aligned to the orbital angular momentum axis. This suggest that E1821+643 merged in hot gas, which is thought to provide an environment where spin alignment from accretion proceeds slowly relative to the merger timescale. We infer that E1821+643 is likely to be rapidly rotating with dimensionless spin ${\chi} = 0.92\pm0.04$. A $\unit[2.6\times10^9]{M_\odot}$ black hole, recoiling from a gas-rich environment at $v\sim \unit[2,240]{km\,s^{-1}}$ is likely to persist as an active galactic nuclei for $\sim \unit[860]{Myr}$, in which time it traverses $\sim \unit[2]{Mpc}$.

Titan's northern high latitudes host many large hydrocarbon lakes. Like water lakes on Earth, Titan's lakes are constantly subject to evaporation. This process strongly affects the atmospheric methane abundance, the atmospheric temperature, the lake mixed layer temperature, and the local wind circulation. In this work we use a 2D atmospheric mesoscale model coupled to a slab lake model to investigate the effect of solar and infrared radiation on the exchange of energy and methane between Titan's lakes and atmosphere. The magnitude of solar radiation reaching the surface of Titan through its thick atmosphere is only a few W/m2. However, we find that this small energy input is important and is comparable in absolute magnitude to the latent and sensible heat fluxes, as suggested in the prior study by S. Rafkin and A. Soto (Icarus, 2020). The implementation of a gray radiative scheme in the model confirms the importance of radiation when studying lakes at the surface of Titan. Solar and infrared radiation change the energy balance of the system leading to an enhancement of the methane evaporation rate, an increase of the equilibrium lake temperature almost completely determined by its environment (humidity, insolation, and background wind), and a strengthening of the local sea breeze, which undergoes diurnal variations. The sea breeze efficiently transports methane vapor horizontally, from the lake to the land, and vertically due to rising motion along the sea breeze front and due to radiation-induced turbulence over the land.

Extended and delayed emission around distant TeV sources induced by the effects of propagation of gamma rays through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). We search for delayed GeV emission from the hard-spectrum TeV blazar 1ES 0229+200 with the goal to detect or constrain the IGMF-dependent secondary flux generated during the propagation of TeV gamma rays through the intergalactic medium. We analyze the most recent MAGIC observations over a 5 year time span and complement them with historic data of the H.E.S.S. and VERITAS telescopes along with a 12-year long exposure of the Fermi/LAT telescope. We use them to trace source evolution in the GeV-TeV band over one-and-a-half decade in time. We use Monte Carlo simulations to predict the delayed secondary gamma-ray flux, modulated by the source variability, as revealed by TeV-band observations. We then compare these predictions for various assumed IGMF strengths to all available measurements of the gamma-ray flux evolution. We find that the source flux in the energy range above 200 GeV experiences variations around its average on the 14 years time span of observations. No evidence for the flux variability is found in 1-100 GeV energy range accessible to Fermi/LAT. Non-detection of variability due to delayed emission from electromagnetic cascade developing in the intergalactic medium imposes a lower bound of B>1.8e-17 G for long correlation length IGMF and B>1e-14 G for an IGMF of the cosmological origin. Though weaker than the one previously derived from the analysis of Fermi/LAT data, this bound is more robust, being based on a conservative intrinsic source spectrum estimate and accounting for the details of source variability in the TeV energy band. We discuss implications of this bound for cosmological magnetic fields which might explain the baryon asymmetry of the Universe.

Due to the expansion of our Universe, the redshift of distant objects changes with time. Although the amplitude of this redshift drift is small, it will be measurable with a decade-long campaigns on the next generation of telescopes. Here we present an alternative view of the redshift drift which captures the expansion of the universe in single epoch observations of the multiple images of gravitationally lensed sources. Considering a sufficiently massive lens, with an associated time delay of order decades, simultaneous photons arriving at a detector would have been emitted decades earlier in one image compared to another, leading to an instantaneous redshift difference between the images. We also investigate the effect of peculiar velocities on the redshift difference in the observed images. Whilst still requiring the observational power of the next generation of telescopes and instruments, the advantage of such a single epoch detection over other redshift drift measurements is that it will be less susceptible to systematic effects that result from requiring instrument stability over decade-long campaigns.

Early dark energy models have attracted attention in the context of the recent problem of the Hubble tension. Here we extend these models by taking into account the new density fluctuations generated by the dark energy decays around the recombination phase. We solve the evolution of the density perturbations in dark energy fluid generated at the phase transition of early dark energy as isocurvature perturbations. Assuming that the isocurvature mode is characterized by a power-law power spectrum and is uncorrelated with the standard adiabatic mode, we calculate the CMB angular power spectra and compare them to the Planck data using the Markov-Chain Monte Carlo method. As a result, we obtained the zero-consistent values of the EDE parameters and $H_0=67.47^{+1.01}_{-1.00}~\mathrm{km} \, \mathrm{s}^{-1} \mathrm{Mpc}^{-1}$ at $68 \%$ CL. This $H_0$ value is almost the same as the Planck+$\Lambda$CDM value, $H_0=67.36 \pm 0.54~\mathrm{km} \, \mathrm{s}^{-1} \mathrm{Mpc}^{-1}$, and there is still a $\sim 3 \sigma$ tension between the CMB and Type Ia supernovae observation. Moreover, the amplitude of the spectra induced by the phase transition of the EDE is constrained to be less than one percent of that of the adiabatic mode. This is so small that such non-standard fluctuations cannot appear in the CMB angular spectra. In conclusion, the significant contributions of the EDE to the background and perturbation are excluded by the analysis using the Planck data only, although there is a need for further analysis including the data from the large scale structure of the universe.

The small semi-major axes of Hot Jupiters lead to high atmospheric temperatures of up to several thousand Kelvin. Under these conditions, thermally ionised metals provide a rich source of charged particles and thus build up a sizeable electrical conductivity. Subsequent electromagnetic effects, such as the induction of electric currents, Ohmic heating, magnetic drag, or the weakening of zonal winds have thus far been considered mainly in the framework of a linear, steady-state model of induction. For Hot Jupiters with an equilibrium temperature $T_{eq} > 1500$ K, the induction of atmospheric magnetic fields is a runaway process that can only be stopped by non-linear feedback. For example, the back-reaction of the magnetic field onto the flow via the Lorentz force or the occurrence of magnetic instabilities. Moreover, we discuss the possibility of self-excited atmospheric dynamos. Our results suggest that the induced atmospheric magnetic fields and electric currents become independent of the electrical conductivity and the internal field, but instead are limited by the planetary rotation rate and wind speed. As an explicit example, we characterise the induction process for the hottest exoplanet, KELT-9b by calculating the electrical conductivity along atmospheric $P-T$-profiles for the day- and nightside. Despite the temperature varying between 3000 K and 4500 K, the resulting electrical conductivity attains an elevated value of roughly 1 S/m throughout the atmosphere. The induced magnetic fields are predominately horizontal and might reach up to a saturation field strength of 400 mT, exceeding the internal field by two orders of magnitude.

The focus of this study is to reveal the reason behind a scale problem detected around 1849 in the historical version of the International Sunspot Number Series, i.e. version 1 (Leussu et al, Astronomy and Astrophysics, 559, A28, 2013; Friedli, Solar Phys.291, 2505, 2016). From 1826 to 1848 Heinrich Schwabe's observations were considered primary by Rudolf Wolf, and a shift of primary observer from Schwabe to Wolf in 1849 seems to have led to an inconsistency in the Sunspot Number series. In this study we benefited from various datasets, the most important being Schwabe's raw counts from the Mittheilungen (Prof. Wolf's Journals) that have been digitised at the Royal Observatory of Belgium between 2017 and 2019. We provide a robust quantification of the detected problem by using classic algebraic calculations but also different methods such as a method inspired by Lockwood et al (Journal of Geophysical Research (Space Physics), 119(7), 5172, 2014), hence assigning a modern k-factor to Schwabe's observations before 1849. We also assess the implications of this 1849 inconsistency on the International Sunspot Number series (Versions 1 and 2) before and after 1849.

Spectroscopic observations in X-ray wavelengths provide excellent diagnostics of the temperature distribution in solar flare plasma. The Solar X-ray Monitor (XSM) onboard the Chandrayaan-2 mission provides broad-band disk integrated soft X-ray solar spectral measurements in the energy range of 1-15 keV with high spectral resolution and time cadence. In this study, we analyse X-ray spectra of three representative GOES C-class flares obtained with the XSM to investigate the evolution of various plasma parameters during the course of the flares. Using the soft X-ray spectra consisting of the continuum and well-resolved line complexes of major elements like Mg, Si, and Fe, we investigate the validity of the isothermal and multi-thermal assumptions on the high temperature components of the flaring plasma. We show that the soft X-ray spectra during the impulsive phase of the high intensity flares are inconsistent with isothermal models and are best fitted with double peaked differential emission measure distributions where the temperature of the hotter component rises faster than that of the cooler component. The two distinct temperature components observed in DEM models during the impulsive phase of the flares suggest the presence of the directly heated plasma in the corona and evaporated plasma from the chromospheric footpoints. We also find that the abundances of low FIP elements Mg, Si, and Fe reduces from near coronal to near photospheric values during the rising phase of the flare and recovers back to coronal values during decay phase, which is also consistent with the chromospheric evaporation scenario.

The absolute position of Sgr A*, the compact radio source at the center of the Milky Way, had been uncertain by several tens of milliarcseconds. Here we report improved astrometric measurements of the absolute position and proper motion of Sgr A*. Three epochs of phase-referencing observations were conducted with the VLBA for Sgr A* at 22 and 43 GHz in 2019 and 2020. Using extragalactic radio sources with sub-milliarcsecond accurate positions as reference, we determined the absolute position of Sgr A* at a reference epoch 2020.0 to be at $\alpha$(J2000) = $17^{\rm h} 45^{\rm m}40.^{\rm s}032863~\pm~0.^{\rm s}000016$ and $\delta$(J2000) = $-29^{\circ} 00^{\prime} 28.^{''}24260~\pm~0.^{''}00047$, with an updated proper motion $-3.152~\pm~0.011$ and $-5.586~\pm~0.006$ mas yr$^{-1}$ in the easterly and northerly directions, respectively.

The sunlight reflected from the Moon during a total lunar eclipse has been transmitted through the Earth's atmosphere on the way to the Moon. The combination of multiple scattering and inhomogeneous atmospheric characteristics during that transmission can potentially polarize that light. A similar (although much smaller) effect should also be observable from the atmosphere of a transiting exoplanet. We present the results of polarization observations during the first 15 minutes of totality of the lunar eclipse of 2022 May 16. We find degrees of polarization of 2.1 +/- 0.4 per cent in B, 1.2 +/- 0.3 per cent in V, 0.5 +/- 0.2 per cent in R and 0.2 +/- 0.2 per cent in I. Our polarization values lie in the middle of the range of those reported for previous eclipses, providing further evidence that the induced polarization can change from event to event. We found no significant polarization difference (<0.02 per cent) between a region of dark Mare and nearby bright uplands or between the lunar limb and regions closer to the disk centre due to the different angle of incidence. This further strengthens the interpretation of the polarization's origin being due to scattering in the Earth's atmosphere rather than by the lunar regolith.

Whereas repulsive gravity was considered as a fringe concept until the mid-1990's, the growingexperimental evidence since this epoch for repulsive gravity, in what is now called Dark Energy,for lack of a better understanding of its nature, has led to a vast literature in order to attemptto characterize this repulsive component, and notably its equation of state. In the following, Iwill show that we can use cosmology to test the hypothesis that antimatter is at the origin ofrepulsive gravity, may play the role of a Dark Energy component and, more surprisingly, maymimic the presence of Dark Matter, and justify the MOND phenomenology. More directly,three experiments, AEgIS, ALPHA-g and Gbar, are attempting to measure the action ofgravitation on cold atoms of antihydrogen at CERN in a near future. Finally, I note thatCP violation might be explained by antigravity and I briefly recall the motivations for thisassertion.

We report the results of timing observations of PSR J1952+2630, a 20.7 ms pulsar in orbit with a massive white dwarf companion. With the increased timing baseline, we obtain improved estimates for astrometric, spin, and binary parameters for this system. We get an improvement of an order of magnitude on the proper motion, and, for the first time, we detect three post-Keplerian parameters in this system: the advance of periastron, the orbital decay, and the Shapiro delay. We constrain the pulsar mass to 1.20$^{+0.28}_{-0.29}\rm M_{\odot}$ and the mass of its companion to 0.97$^{+0.16}_{-0.13}\rm M_{\odot}$. The current value of $\dot{P}_{\rm b}$ is consistent with GR expectation for the masses obtained using $\dot{\omega}$ and $h_3$. The excess represents a limit on the emission of dipolar GWs from this system. This results in a limit on the difference in effective scalar couplings for the pulsar and companion (predicted by scalar-tensor theories of gravity; STTs) of $|\alpha_{\rm p}-\alpha_{\rm c}| < 4.8 \times 10^{-3}$, which does not yield a competitive test for STTs. However, our simulations of future campaigns of this system show that by 2032, the precision of $\dot{P}_{\rm b}$ and $\dot{\omega}$ will allow for much more precise masses and much tighter constraints on the orbital decay contribution from dipolar GWs, resulting in $|\alpha_{\rm p}-\alpha_{\rm c}|<1.3 \times 10^{-3}$. We also present the constraints this system will place on the $\{\alpha_0,\beta_0\}$ parameters of DEF gravity by 2032. They are comparable to those of PSR J1738+0333. Unlike PSR J1738+0333, PSR J1952+2630 will not be limited in its mass measurement and has the potential to place even more restrictive limits on DEF gravity in the future. Further improvements to this test will likely be limited by uncertainties in the kinematic contributions to $\dot{P}_{\rm b}$ due to lack of precise distance measurements.

Scalar perturbations in the early Universe create over-dense regions that can collapse into primordial black holes (PBH). This process emits scalar-induced gravitational waves (SIGW) that behaves like an extra radiation component and contributes to the relativistic degrees of freedom ($N_{\rm{eff}}$). We show that $N_{\rm{eff}}$ limits from cosmic microwave background (CMB) give promising sensitivities on both the abundance of PBHs and the primordial curvature perturbation ($\mathcal{P}_{\mathcal{R}}(k)$) at small scales. We show that {\it Planck} and ACTPol data can exclude supermassive PBHs with peak mass $M_{\bullet} \in [3 \times 10^{5}, 5 \times 10^{10}] {\rm{M}}_{\odot}$ as the major component of dark matter, depending on the shape of the PBHs mass distribution. Future CMB-S4 mission is capable of broadening this limit to a vast PBH mass window of $M_{\bullet} \in [8 \times 10^{-5}, 5 \times 10^{10}] {\rm{M}}_{\odot}$, covering sub-stellar masses. These limits correspond to the enhanced sensitivity of $\mathcal{P}_{\mathcal{R}}(k)$ on scales of $k \in [10^1, 10^{22}]\ \rm{Mpc^{-1}}$, which is much smaller than those scales probed by direct perturbation power spectra (CMB and large-scale structure).

We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14-195keV; Swift/BAT) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy and two-phase torus models and clumpy disc+wind models. We find that the smooth, clumpy and two-phase torus models (i.e. without including the polar dusty wind component) and disc+wind models provide best fits for a comparable number of galaxies, 8/24 (33.3%) and 9/24 (37.5%), respectively. We find that the best-fitted models depend on the hydrogen column density (NH), which is related to the X-ray (unobscured/obscured) and/or optical (Sy1/Sy2) classification. In particular, smooth, clumpy and two-phase torus models best reproduce the infrared (IR) emission of AGN with relatively high hydrogen column density (median value of log (NH)=23.5+-0.8; i.e. Sy2). However, clumpy disc+wind models provide the best fits to the nuclear IR spectral energy distributions (SEDs) of Sy1/1.8/1.9 (median value of log (NH)=21.0+-1.0), specifically in the near-IR (NIR) range. The success of the disc+wind models in fitting the NIR emission of Sy1 galaxies is due to the combination of adding large graphite grains to the dust composition and self-obscuration effects caused by the wind at intermediate inclinations. In general, we find that the Seyfert galaxies having unfavourable (favourable) conditions, i.e. nuclear hydrogen column density and Eddington ratio, for launching IR dusty polar outflows are best-fitted with smooth, clumpy and two-phase torus (disk+wind) models confirming the predictions from simulations. Therefore, our results indicate that the nature of the inner dusty structure in AGN depend on the intrinsic AGN properties.

Environment has long been known to have significant impact on the evolution of galaxies, but here we seek to quantify the subtler differences that might be found in disk galaxies, depending on whether they are isolated, the most massive galaxy in a group (centrals), or a lesser member (satellites). The MaNGA survey allows us to define a large mass-matched sample of 574 galaxies with high-quality integrated spectra in each category. Initial examination of their spectral indices indicates significant differences, particularly in low-mass galaxies. Semi-analytic spectral fitting of a full chemical evolution model to these spectra confirms these differences, with low-mass satellites having a shorter period of star formation and chemical enrichment typical of a closed box, while central galaxies have more extended histories, with evidence of on-going gas accretion over their lifetimes. The derived parameters for gas infall timescale and wind strength suggest that low-mass satellite galaxies have their hot halos of gas effectively removed, while central galaxies retain a larger fraction of gas than isolated galaxies due to the deeper group potential well in which they sit. S0 galaxies form a distinct subset within the sample, particularly at higher masses, but do not bias the inferred lower-mass environmental impact significantly. The consistent picture that emerges underlines the wealth of archaeological information that can be extracted from high-quality spectral data using techniques like semi-analytic spectral fitting.

Various bright structures abound in the chromosphere playing an essential role in the dynamics and evolution therein. Tentatively identifying the wave characteristics in the outer solar atmosphere helps to understand this layer better. One of the most significant aspects of these characteristics is the wave phase speed (PS), which is a dominant contribution to solar coronal heating and Energy distribution in the Sun's atmosphere layers. To obtain energy flux (EF), it is necessary to calculate the filling factor (FF) and the PS. In this study, the FF was determined by tracking the size and intensity of the IRIS bright points (BPs). To estimate an accurate PS and EF, it is necessary to know the chromosphere and transition region (TR) thickness and the phase difference between the two desired levels. chromosphere and TR thickness cannot be measured directly on the disc; This study is performed using spectral data and calibrated based on Doppler velocities. As a result, the PSs in AR and CH, as well as for IRIS BPs have been calculated using the cross-power wavelet transform of Doppler velocities. Consequently, about CH, the PS mean values are from 40 to 180 km/s at network and from 30 to 140 km/s at internetwork; And about AR, are from 80 to 540 km/s at network and 70 to 220 km/s at internetwork. Finally, the EF for the IRIS BPs has been calculated in three different frequencies. The results indicate that the network BPs have an influential role in heating the higher layers while in the internetwork BPs, most of the energy returns to the lower layers.

One of the most fundamental baryonic matter components of galaxies is the neutral atomic hydrogen (H$\,{\rm \scriptsize I}$). At low redshifts, this component can be traced directly through the 21-cm transition, but to infer H$\,{\rm \scriptsize I}$ gas content of the most distant galaxies, a viable tracer is needed. We here investigate the fidelity of the fine structure transition of the ($^2P_{3/2} - ^2P_{1/3}$) transition of singly-ionized carbon [C$\,{\rm \scriptsize II}$] at $158\,\mu$m as a proxy for H$\,{\rm \scriptsize I}$ in a set simulated galaxies at $z\approx 6$, following the work by Heintz et al. (2021). We select 11,125 star-forming galaxies from the SIMBA simulations, with far-infrared line emissions post-processed and modeled within the SIGAME framework. We find a strong connection between [C$\,{\rm \scriptsize II}$] and H$\,{\rm \scriptsize I}$, with the relation between this [C$\,{\rm \scriptsize II}$]-to-H$\,{\rm \scriptsize I}$ relation ($\beta_{\rm [C\,{\rm \scriptsize II}]}$) being anti-correlated with the gas-phase metallicity of the simulated galaxies. We further use these simulations to make predictions for the total baryonic matter content of galaxies at $z\approx 6$, and specifically the HI gas mass fraction. We find mean values of $M_{\rm HI}/M_\star = 1.4$, and $M_{\rm HI}/M_{\rm bar,tot} = 0.45$. These results provide strong evidence for H$\,{\rm \scriptsize I}$ being the dominant baryonic matter component by mass in galaxies at $z\approx 6$.

The probability distribution of density in isothermal, supersonic, turbulent gas is approximately lognormal. This behaviour can be traced back to the shock waves travelling through the medium, which randomly adjust the density by a random factor of the local sonic Mach number squared. Provided a certain parcel of gas experiences a large number of shocks, due to the central limit theorem, the resulting distribution for density is lognormal. We explore a model in which parcels of gas undergo finite number of shocks before relaxing to the ambient density, causing the distribution for density to deviate from a lognormal. We confront this model with numerical simulations with various r.m.s. Mach numbers ranging from subsonic as low as 0.1 to supersonic at 25. We find that the fits to the finite formula are an order of magnitude better than a lognormal. The model naturally extends even to subsonic flows, where no shocks exist.

Thermonuclear (type-I) bursts arise from unstable ignition of accumulated fuel on the surface of neutron stars in low-mass X-ray binaries. Measurements of burst properties in principle enable observers to infer the properties of the host neutron star and mass donors, but a number of confounding astrophysical effects contribute to systematic uncertainties. Here we describe some commonly-used approaches for determining system parameters, including composition of the burst fuel, and introduce a new suite of software tools, concord, intended to fully account for astrophysical uncertainties. Comparison of observed burst properties with the predictions of numerical models is a complementary method of constraining host properties, and the tools presented here are intended to make comprehensive model-observation comparisons straightforward. When combined with the extensive samples of burst observations accumulated by X-ray observatories, these software tools will provide a step-change in the amount of information that can be inferred about typical burst sources.

We report the properties of more than 600 bursts (including cluster-bursts) detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during an extremely active episode on UTC September 25-28, 2021, in a series of four papers. The observations were carried out in the band of 1.0 - 1.5 GHz by using the center beam of the L-band 19-beam receiver. We monitored the source in sixteen 1-hour sessions and one 3-hour session spanning 23 days. All the bursts were detected during the first four days. In this first paper of the series, we perform a detailed morphological study of 624 bursts using the 2-dimensional frequency-time ``waterfall'' plots, with a burst (or cluster-burst) defined as an emission episode during which the adjacent emission peaks have a separation shorter than 400 ms. The duration of a burst is therefore always longer than 1 ms, with the longest up to more than 120 ms. The emission spectra of the sub-bursts are typically narrow within the observing band with a characteristic width of $\sim$277 MHz. The center frequency distribution has a dominant peak at about 1091.9 MHz and a secondary weak peak around 1327.9 MHz. Most bursts show a frequency-downward-drifting pattern. Based on the drifting patterns, we classify the bursts into five main categories: downward drifting (263) bursts, upward drifting (3) bursts, complex (203), no drifting (35) bursts, and no evidence for drifting (121) bursts. Subtypes are introduced based on the emission frequency range in the band (low, middle, high and wide) as well as the number of components in one burst (1, 2, or multiple). We measured a varying scintillation bandwidth from about 0.5 MHz at 1.0 GHz to 1.4 MHz at 1.5 GHz with a spectral index of 3.0.

As the third paper in the multiple-part series, we report the statistical properties of radio bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio telescope (FAST) during an extremely active episode between the 25th and the 28th of September 2021 (UT). We focus on the polarisation properties of 536 bright bursts with $\mathrm{S/N}>50$. We found that the Faraday rotation measures (RMs) monotonically dropped from $-579 \ {\rm rad \ m^{-2}}$ to $-605 \ {\rm rad \ m^{-2}}$ in the 4-day window. The RM values were compatible with the values ($-300$ to $-900\ {\rm rad \ m^{-2}}$ ) reported 4 month ago (Xu et al. 2022). However, the RM evolution rate in the current observation window was at least an order of magnitude smaller than the one ($\sim 500\ {\rm rad \ m^{-2}\, day^{-1}}$) previously reported during the rapid RM-variation phase, but is still higher than the one ($\le 1\ {\rm rad \ m^{-2} day^{-1}}$ ) during the later RM no-evolution phase. The bursts of FRB 20201124A were highly polarised with the total degree of polarisation (circular plus linear) greater than 90% for more than 90\% of all bursts. The distribution of linear polarisation position angles (PAs), degree of linear polarisation ($L/I$), and degree of circular polarisation ($V/I$) can be characterised with unimodal distribution functions. During the observation window, the distributions became wider with time, i.e. with larger scatter, but the centroids of the distribution functions remained nearly constant. For individual bursts, significant PA variations (confidence level 5-$\sigma$) were observed in 33% of all bursts. The polarisation of single pulses seems to follow certain complex trajectories on the Poincar\'e sphere, which may shed light on the radiation mechanism at the source or the plasma properties along the path of FRB propagation.

We report the properties of more than 800 bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio telescope (FAST) during an extremely active episode on UTC September 25th-28th, 2021 in a series of four papers. In this fourth paper of the series, we present a systematic search of the spin period and linear acceleration of the source object from both 996 individual pulse peaks and the dedispersed time series. No credible spin period was found from this data set. We rule out the presence of significant periodicity in the range between 1 ms to 100 s with a pulse duty cycle $< 0.49\pm0.08$ (when the profile is defined by a von-Mises function, not a boxcar function) and linear acceleration up to $300$ m s$^{-2}$ in each of the four one-hour observing sessions, and up to $0.6$ m s$^{-2}$ in all 4 days. These searches contest theoretical scenarios involving a 1 ms to 100 s isolated magnetar/pulsar with surface magnetic field $<10^{15}$ G and a small duty cycle (such as in a polar-cap emission mode) or a pulsar with a companion star or black hole up to 100 M$_{\rm \odot}$ and $P_b>10$ hours. We also perform a periodicity search of the fine structures and identify 53 unrelated millisecond-timescale "periods" in multi-components with the highest significance of 3.9 $\sigma$. The "periods" recovered from the fine structures are neither consistent nor harmonically related. Thus they are not likely to come from a spin period. We caution against claiming spin periodicity with significance below $\sim$ 4 $\sigma$ with multi-components from one-off FRBs. We discuss the implications of our results and the possible connections between FRB multi-components and pulsar micro-structures.

Knowledge of the three-dimensional structure of Galactic molecular clouds is important for understanding how clouds are affected by processes such as turbulence and magnetic fields and how this structure effects star formation within them. Great progress has been made in this field with the arrival of the Gaia mission, which provides accurate distances to $\sim10^{9}$ stars. Combining these distances with extinctions inferred from optical-IR, we recover the three-dimensional structure of 16 Galactic molecular cloud complexes at $\sim1$pc resolution using our novel three-dimensional dust mapping algorithm \texttt{Dustribution}. Using \texttt{astrodendro} we derive a catalogue of physical parameters for each complex. We recover structures with aspect ratios between 1 and 11, i.e.\ everything from near-spherical to very elongated shapes. We find a large variation in cloud environments that is not apparent when studying them in two-dimensions. For example, the nearby California and Orion A clouds look similar on-sky, but we find California to be more sheet-like, and massive, which could explain their different star-formation rates. In Carina, our most distant complex, we observe evidence for dust sputtering, which explains its measured low dust mass. By calculating the total mass of these individual clouds, we demonstrate that it is necessary to define cloud boundaries in three-dimensions in order to obtain an accurate mass; simply integrating the extinction overestimates masses. We find that Larson's relationship on mass vs radius holds true whether you assume a spherical shape for the cloud or take their true extents.

The observational data of primordial black holes and scalar-induced gravitational waves can constrain the primordial curvature perturbation at small scales. We parameterize the primordial curvature perturbation by a broken power law form and find that it is consistent with many inflation models that can produce primordial black holes, such as nonminimal derivative coupling inflation, scalar-tensor inflation, Gauss-Bonnet inflation, and K/G inflation. The constraints from primordial black holes on the primordial curvature power spectrum with the broken power law form are obtained, where the fraction of primordial black holes in dark matter is calculated by the peak theory. Both the real-space top-hat and the Gauss window functions are considered. The constraints on the amplitude of primordial curvature perturbation with Gauss window function are around three times larger than those with real-space top-hat window function. The constraints on the primordial curvature perturbation from the NANOGrav 12.5yrs data sets are displayed, where the NANOGrav signals are assumed as the scalar-induced gravitational waves, and only the first five frequency bins are used.

The Galactic bulge at latitude $4<|b|$(deg.)$<10$ was claimed to show an X-shape, which means that stellar density distributions along the line of sight have a double peak. However, this double peak is only observed with the red-clump population, and doubt has been cast on its use as a perfect standard candle.\ As such, a boxy bulge without an X-shape is not discarded. We aim to constrain the shape of the bulge making use of a different population: Mira variables from the new Optical Gravitational Lensing Experiment data release, OGLE-IV, with an average age of 9 Gyr. We analysed an area of the bulge far from the plane, where we fitted the density of the Miras with boxy bulge and X-shaped bulge models and calculated the probability of each model. We find that the probability of a boxy bulge fitting the data is $p=0.19$, whereas the probability for the X-shaped bulge is only $p=2.85 \cdot 10^{-6}$ (equivalent to a tension of the model with the data of a 4.7$\sigma $ level). Therefore, the boxy bulge model seems to be more appropriate for describing the Galactic bulge, although we cannot exclude any model with complete certainty.

We report the properties of more than 800 bursts detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during an extremely active episode on UTC September 25-28, 2021 in a series of four papers. In this second paper of the series, we mainly focus on the energy distribution of the detected bursts. The event rate initially increased exponentially but the source activity stopped within 24 hours after the 4th day. The detection of 542 bursts in one hour during the fourth day marked the highest event rate detected from one single FRB source so far. The bursts have complex structures in the time-frequency space. We find a double-peak distribution of the waiting time, which can be modeled with two log-normal functions peaking at 51.22 ms and 10.05 s, respectively. Compared with the emission from a previous active episode of the source detected with FAST, the second distribution peak time is smaller, suggesting that this peak is defined by the activity level of the source. We calculate the isotropic energy of the bursts using both a partial bandwidth and a full bandwidth and find that the energy distribution is not significantly changed. We find that an exponentially connected broken-power-law function can fit the cumulative burst energy distribution well, with the lower and higher-energy indices being $-1.22\pm0.01$ and $-4.27\pm0.23$, respectively. Assuming a radio radiative efficiency of $\eta_r = 10^{-4}$, the total isotropic energy of the bursts released during the four days when the source was active is already $3.9\times10^{46}$ erg, exceeding $\sim 23\%$ of the available magnetar dipolar magnetic energy. This challenges the magnetar models invoking an inefficient radio emission (e.g. synchrotron maser models).

We present the development of a Skipper Charge-Coupled Device (CCD) focal plane prototype for the SOAR Telescope Integral Field Spectrograph (SIFS). This mosaic focal plane consists of four 6k $\times$ 1k, 15 $\mu$m pixel Skipper CCDs mounted inside a vacuum dewar. We describe the process of packaging the CCDs so that they can be easily tested, transported, and installed in a mosaic focal plane. We characterize the performance of $\sim 650 \mu$m thick, fully-depleted engineering-grade Skipper CCDs in preparation for performing similar characterization tests on science-grade Skipper CCDs which will be thinned to 250$\mu$m and backside processed with an antireflective coating. We achieve a single-sample readout noise of $4.5 e^{-} rms/pix$ for the best performing amplifiers and sub-electron resolution (photon counting capabilities) with readout noise $\sigma \sim 0.16 e^{-} rms/pix$ from 800 measurements of the charge in each pixel. We describe the design and construction of the Skipper CCD focal plane and provide details about the synchronized readout electronics system that will be implemented to simultaneously read 16 amplifiers from the four Skipper CCDs (4-amplifiers per detector). Finally, we outline future plans for laboratory testing, installation, commissioning, and science verification of our Skipper CCD focal plane.

Galaxy Zoo: Clump Scout is a web-based citizen science project designed to identify and spatially locate giant star forming clumps in galaxies that were imaged by the Sloan Digital Sky Survey Legacy Survey. We present a statistically driven software framework that is designed to aggregate two-dimensional annotations of clump locations provided by multiple independent Galaxy Zoo: Clump Scout volunteers and generate a consensus label that identifies the locations of probable clumps within each galaxy. The statistical model our framework is based on allows us to assign false-positive probabilities to each of the clumps we identify, to estimate the skill levels of each of the volunteers who contribute to Galaxy Zoo: Clump Scout and also to quantitatively assess the reliability of the consensus labels that are derived for each subject. We apply our framework to a dataset containing 3,561,454 two-dimensional points, which constitute 1,739,259 annotations of 85,286 distinct subjects provided by 20,999 volunteers. Using this dataset, we identify 128,100 potential clumps distributed among 44,126 galaxies. This dataset can be used to study the prevalence and demographics of giant star forming clumps in low-redshift galaxies. The code for our aggregation software framework is publicly available at: https://github.com/ou-astrophysics/BoxAggregator

Dilatons (and moduli) couple to the masses and coupling constants of ordinary matter, and these quantities are fixed by the local value of the dilaton field. If, in addition, the dilaton with mass $m_\phi$ contributes to the cosmic dark matter density, then such quantities oscillate in time at the dilaton Compton frequency. We show how these oscillations lead to broadening and shifting of the Voigt profile of the Ly$\alpha$ forest, in a manner that is correlated with the local dark matter density. We further show how tomographic methods allow the effect to be reconstructed by observing the Ly$\alpha$ forest spectrum of distant quasars. We then simulate a large number of quasar lines of sight using the lognormal density field, and forecast the ability of future astronomical surveys to measure this effect. We find that in the ultra low mass range $10^{-32}\text{ eV}\leq m_\phi\leq 10^{-28}\text{ eV}$ upcoming observations can improve over existing limits to the dilaton electron mass and fine structure constant couplings set by fifth force searches by up to five orders of magnitude. Our projected limits apply assuming that the ultralight dilaton makes up a few percent of the dark matter density, consistent with upper limits set by the cosmic microwave background anisotropies.

Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process in which the window functions play a crucial role. In a delay-based analysis, the mapping between instrumental space, made of per-baseline delay spectra, and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand the power spectrum measured by an interferometer, we assess the impact of instrument characteristics and analysis choices on the estimator by deriving its exact window functions, outside of the delay approximation. Focusing on HERA as a case study, we find that observations made with long baselines tend to correspond to enhanced low-k tails of the window functions, which facilitate foreground leakage outside the wedge, whilst the choice of bandwidth and frequency taper can help narrow them down. With the help of simple test cases and more realistic visibility simulations, we show that, apart from tracing mode mixing, the window functions can accurately reconstruct the power spectrum estimator of simulated visibilities. We note that the window functions depend strongly on the chromaticity of the beam, and less on its spatial structure - a Gaussian approximation, ignoring side lobes, is sufficient. Finally, we investigate the potential of asymmetric window functions, down-weighting the contribution of low-k power to avoid foreground leakage. The window functions presented in this work correspond to the latest HERA upper limits for the full Phase I data. They allow an accurate reconstruction of the power spectrum measured by the instrument and can be used in future analyses to confront theoretical models and data directly in cylindrical space.

Scattered light imaging studies have detected nearly two dozen spiral arm systems in circumstellar disks, yet the formation mechanisms for most of them are still under debate. Although existing studies can use motion measurements to distinguish leading mechanisms such as planet-disk interaction and disk self-gravity, close-in stellar flybys can induce short-lived spirals and even excite arm-driving planets into highly eccentric orbits. With unprecedented stellar location and proper motion measurements from Gaia DR3, here we study for known spiral arm systems their flyby history with their stellar neighbours by formulating an analytical on-sky flyby framework. For stellar neighbors currently located within 10 pc from the spiral hosts, we restrict the flyby time to be within the past $10^4$ yr and the flyby distance to be within $10$ times the disk extent in scattered light. Among a total of $12570$ neighbors that are identified in Gaia DR3 for $20$ spiral systems, we do not identify credible flyby candidates for isolated systems. Our analysis suggests that close-in recent flyby is not the dominant formation mechanism for isolated spiral systems in scattered light.

We construct globally hyperbolic spacetimes such that each slice $\{t=t_0\}$ of the universal time $t$ is a model space of constant curvature $k(t_0)$ which may not only vary with $t_0\in\mathbb{R}$ but also change its sign. The metric is smooth and slightly different to FLRW spacetimes, namely, $g=-dt^2+dr^2+ S_{k(t)}^2(r) g_{\mathbb{S}^{n-1}}$, where $g_{\mathbb{S}^{n-1}}$ is the metric of the standard sphere, $S_{k(t)}(r)=\sin(\sqrt{k(t)}\, r)/\sqrt{k(t)}$ when $k(t)\geq 0$ and $S_{k(t)}(r)=\sinh(\sqrt{-k(t)}\, r)/\sqrt{-k(t)}$ when $k(t)\leq 0$. In the open case, the $t$-slices are (non-compact) Cauchy hypersurfaces of curvature $k(t)\leq 0$, thus homeomorphic to $\mathbb{R}^n$; a typical example is $k(t)=-t^2$ (i.e., $S_{k(t)}(r)=\sinh(tr)/t$). In the closed case, $k(t)>0$ somewhere, a slight extension of the class shows how the topology of the $t$-slices changes. This makes at least one comoving observer to disappear in finite time $t$ showing some similarities with an inflationary expansion. Anyway, the spacetime is foliated by Cauchy hypersurfaces homeomorphic to spheres, not all of them $t$-slices.

Neutron star mergers provide a unique laboratory for the study of strong-field gravity coupled to Quantum Chromodynamics in extreme conditions. The frequencies and amplitudes of the resulting gravitational waves encode invaluable information about the merger. Simulations to date have shown that these frequencies lie in the kilo-Hertz range. They have also shown that, if Quantum Chromodynamics possesses a first-order phase transition at high baryon density, then this is likely to be accessed during the merger dynamics. Here we show that this would result in the nucleation of superheated and/or supercompressed bubbles whose subsequent dynamics would produce gravitational waves in the Mega-Hertz range. We estimate the amplitude of this signal and show that it may fall within the expected sensitivity of future superconducting radio-frequency cavity detectors for mergers at distances up to tens of Mega-parsecs.

We summarize recent progress in the development, application, and understanding of effective field theories and highlight promising directions for future research. This Report is prepared as the TF02 "Effective Field Theory" topical group summary for the Theory Frontier as part of the Snowmass 2021 process.

Many regions across the globe broke their surface temperature records in recent years, further sparking concerns about the impending arrival of "tipping points" later in the 21st century. This study analyzes observed global surface temperature trends in three target latitudinal regions: the Arctic Circle, the Tropics, and the Antarctic Circle. We show that global warming is accelerating unevenly across the planet, with the Arctic warming at approximately three times the average rate of our world. We further analyzed the reliability of latitude-dependent surface temperature simulations from a suite of Coupled Model Intercomparison Project Phase 6 models and their multi-model mean. We found that GISS-E2-1-G and FGOALS-g3 were the best-performing models based on their statistical abilities to reproduce observational, latitude-dependent data. Surface temperatures were projected from ensemble simulations of the Shared Socioeconomic Pathway 2-4.5 (SSP2-4.5). We estimate when the climate will warm by 1.5, 2.0, and 2.5 degrees C relative to the preindustrial period, globally and regionally. GISS-E2-1-G projects that global surface temperature anomalies would reach 1.5, 2.0, and 2.5 degrees C in 2024 (+/-1.34), 2039 (+/-2.83), and 2057 (+/-5.03) respectively, while FGOALS-g3 predicts these "tipping points" would arrive in 2024 (+/-2.50), 2054 (+/-7.90), and 2087 (+/-10.55) respectively. Our results reaffirm a dramatic, upward trend in projected climate warming acceleration, with upward concavity in 21st century projections of the Arctic, which could lead to catastrophic consequences across the Earth. Further studies are necessary to determine the most efficient solutions to reduce global warming acceleration and maintain a low SSP, both globally and regionally.

We investigate properties of the conserved charge in general relativity, recently proposed by one of the present authors with his collaborators, in the inflation era, the matter dominated era and the radiation dominated era of the expanding Universe. We show that the conserved charge becomes the Bekenstein-Hawking entropy in the inflation era, and it becomes the matter entropy and the radiation entropy in the matter and radiation dominated eras, respectively, while the charge itself is always conserved. These properties are qualitatively confirmed by a numerical analysis of a model with a scalar field and radiations. Results in this paper provide more evidences on the interpretation that the conserved charge in general relativity corresponds to entropy.

Solar Orbiter provides dust detection capability in inner heliosphere, but estimating physical properties of detected dust from the collected data is far from straightforward. First, a physical model for dust collection considering a Poisson process is formulated. Second, it is shown that dust on hyperbolic orbits is responsible for the majority of dust detections with Solar Orbiter's Radio and Plasma Waves (SolO/RPW). Third, the model for dust counts is fitted to SolO/RPW data and parameters of the dust are inferred, namely: radial velocity, hyperbolic meteoroids predominance, and solar radiation pressure to gravity ratio as well as uncertainties of these. Non-parametric model fitting is used to get the difference between inbound and outbound detection rate and dust radial velocity is thus estimated. A hierarchical Bayesian model is formulated and applied to available SolO/RPW data. The model uses the methodology of Integrated Nested Laplace Approximation, estimating parameters of dust and their uncertainties. SolO/RPW dust observations can be modelled as a Poisson process in a Bayesian framework and observations up to this date are consistent with the hyperbolic dust model with an additional background component. Analysis suggests a radial velocity of the hyperbolic component around $(63 \pm 7) \mathrm{km/s}$ with the predominance of hyperbolic dust about $(78 \pm 4) \%$. The results are consistent with hyperbolic meteoroids originating between $0.02 \mathrm{AU}$ and $0.1 \mathrm{AU}$ and showing substantial deceleration, which implies effective solar radiation pressure to gravity ratio $\gtrsim 0.5$. The flux of hyperbolic component at $1 \mathrm{AU}$ is found to be $(1.1 \pm 0.2) \times 10^{-4} \mathrm{m^{-2}s^{-1}}$ and the flux of background component at $1 \mathrm{AU}$ is found to be $(5.4 \pm 1.5) \times 10^{-5} \mathrm{m^{-2}s^{-1}}$.

We investigate the field equations of the conformally invariant models of gravity with curvature-matter coupling, constructed in Weyl geometry, by using the Palatini formalism. We consider the case in which the Lagrangian is given by the sum of the square of the Weyl scalar, of the strength of the field associated to the Weyl vector, and a conformally invariant geometry-matter coupling term, constructed from the matter Lagrangian and the Weyl scalar. After substituting the Weyl scalar in terms of its Riemannian counterpart, the quadratic action is defined in Riemann geometry, and involves a nonminimal coupling between the Ricci scalar and the matter Lagrangian. For the sake of generality, a more general Lagrangian, in which the Weyl vector is nonminmally coupled with an arbitrary function of the Ricci scalar, is also considered. By varying the action independently with respect to the metric and the connection, the independent connection can be expressed as the Levi-Civita connection of an auxiliary, Ricci scalar and Weyl vector dependent metric, which is related to the physical metric by means of a conformal transformation. The field equations are obtained in both the metric and the Palatini formulations. The cosmological implications of the Palatini field equations are investigated for three distinct models corresponding to different forms of the coupling functions. A comparison with the standard $\Lambda$CDM model is also performed, and we find that the Palatini type cosmological models can give an acceptable description of the observations.

Leaver's method has been the standard for computing the quasinormal mode (QNM) frequencies for a Kerr black hole (BH) for a few decades. We start with a spectral variant of Leaver's method introduced by Cook and Zalutskiy (arXiv: 1410.7698) and propose improvements in the form of computing the necessary derivatives analytically, rather than by numerical finite differencing. We also incorporate this derivative information into qnm, a Python package which finds the QNM frequencies via the spectral variant of Leaver's method. We confine ourselves to first derivatives only.

We discuss some testable predictions of a non-supersymmetric $SO(10)$ model supplemented by a Peccei-Quinn symmetry. We utilize a symmetry breaking pattern of $SO(10)$ that yields unification of the Standard Model gauge couplings, with the unification scale also linked to inflation driven by an $SO(10)$ singlet scalar field with a Coleman-Weinberg potential. Proton decay mediated by the superheavy gauge bosons may be observable at the proposed Hyper-Kamiokande experiment. Due to an unbroken $Z_2$ gauge symmetry from $SO(10)$, the model predicts the presence of a stable intermediate mass fermion which, together with the axion, provides the desired relic abundance of dark matter. The model also predicts the presence of intermediate scale topologically stable monopoles and strings that survive inflation. The monopoles may be present in the Universe at an observable level. We estimate the stochastic gravitational wave background emitted by the strings and show that it should be testable in a number of planned and proposed space and land based experiments. Finally, we show how the observed baryon asymmetry in the Universe is realized via non-thermal leptogenesis.