Papers for Monday, Oct 31 2022

Turbulent radiative mixing layers (TRMLs) play an important role in many astrophysical contexts where cool ($\lesssim 10^4$ K) clouds interact with hot flows (e.g., galactic winds, high velocity clouds, infalling satellites in halos and clusters). The fate of these clouds (as well as many of their observable properties) is dictated by the competition between turbulence and radiative cooling; however, turbulence in these multiphase flows remains poorly understood. We have investigated the emergent turbulence arising in the interaction between clouds and supersonic winds in hydrodynamic ENZO-E simulations. In order to obtain robust results, we employed multiple metrics to characterize the turbulent velocity, $v_{\rm turb}$. We find four primary results, when cooling is sufficient for cloud survival. First, $v_{\rm turb}$ manifests clear temperature dependence. Initially, $v_{\rm turb}$ roughly matches the scaling of sound speed on temperature. In gas hotter than the temperature where cooling peaks, this dependence weakens with time until $v_{\rm turb}$ is constant. Second, the relative velocity between the cloud and wind initially drives rapid growth of $v_{\rm turb}$. As it drops (from entrainment), $v_{\rm turb}$ starts to decay before it stabilizes at roughly half its maximum. At late times cooling flows appear to support turbulence. Third, the magnitude of $v_{\rm turb}$ scales with the ratio between the hot phase sound crossing time and the minimum cooling time. Finally, we find tentative evidence for a length-scale associated with resolving turbulence. Under-resolving this scale may cause violent shattering and affect the cloud's large-scale morphological properties.

Alongside the recent increase in discoveries of tidal disruption events (TDEs) have come an increasing number of ambiguous nuclear transients (ANTs). These ANTs are characterized by hot blackbody-like UV/optical spectral energy distributions (SEDs) and smooth photometric evolution, often with hard powerlaw-like X-ray emission. ANTs are likely exotic TDEs or smooth flares originating in otherwise narrow-line active galactic nuclei (AGNs). While their emission in the UV/optical and X-ray has been relatively well-explored, their infrared (IR) emission has not been studied in detail. Here we use the NEOWISE mission and its low-cadence mapping of the entire sky to study mid-infrared dust reprocessing echoes of ANTs. We study 19 ANTs, finding significant MIR flares in 18 objects for which we can estimate an IR luminosity and temperature evolution. The dust reprocessing echoes show a wide range in IR luminosities ($\sim 10^{42} - 10^{45}$ erg s$^{-1}$) with blackbody temperatures largely consistent with sublimation temperature of graphite grains. Excluding the two sources possibly associated with luminous supernovae (ASASSN-15lh and ASASSN-17jz), the dust covering fractions (f$_c$) for detected IR flares lie between 0.05 and 0.91, with a mean of f$_c$ = 0.29 for all ANTs (including limits) and f$_c$ = $0.38 \pm 0.04$ for detections. These covering fractions are much higher than optically-selected TDEs and similar to AGNs. We interpret the high covering fractions in ANT host galaxies as evidence for the presence of a dusty torus.

We perform magnetohydrodynamic simulations of accreting, equal-mass binary black holes in full general relativity focusing on the effect of spin and minidisks on the accretion rate and Poynting luminosity variability. We report on the structure of the minidisks and periodicities in the mass of the minidisks, mass accretion rates, and Poynting luminosity. The accretion rate exhibits a quasi-periodic behavior related to the orbital frequency of the binary in all systems that we study, but the amplitude of this modulation is dependent on the existence of persistent minidisks. In particular, systems that are found to produce persistent minidisks have a much weaker modulation of the mass accretion rate, indicating that minidisks can increase the inflow time of matter onto the black holes, and dampen out the quasi-periodic behavior. This finding has potential consequences for binaries at greater separations where minidisks can be much larger and may dampen out the periodicities significantly.

The ultraviolet (UV) emission of stellar flares may have a pivotal role in the habitability of rocky exoplanets around low-mass stars. Previous studies have used white-light observations to calibrate empirical models which describe the optical and UV flare emission. However, the accuracy of the UV predictions of models have previously not been tested. We combined TESS optical and GALEX UV observations to test the UV predictions of empirical flare models calibrated using optical flare rates of M stars. We find that the canonical 9000 K blackbody model used by flare studies underestimates the GALEX NUV energies of field age M stars by up to a factor of 6.5$\pm$0.7 and the GALEX FUV energies of fully convective field age M stars by 30.6$\pm$10.0. We calculated energy correction factors that can be used to bring the UV predictions of flare models closer in line with observations. We calculated pseudo-continuum flare temperatures that describe both the white-light and GALEX NUV emission. We measured a temperature of 10,700 K for flares from fully convective M stars after accounting for the contribution from UV line emission. We also applied our correction factors to the results of previous studies of the role of flares in abiogenesis. Our results show that M stars do not need to be as active as previously thought in order to provide the NUV flux required for prebiotic chemistry, however we note that flares will also provide more FUV flux than previously modelled.

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will discover an unprecedented number of supernovae (SNe), making spectroscopic classification for all the events infeasible. LSST will thus rely on photometric classification, whose accuracy depends on the not-yet-finalized LSST observing strategy. In this work, we analyze the impact of cadence choices on classification performance using simulated multi-band light curves. First, we simulate SNe with an LSST baseline cadence, a non-rolling cadence, and a presto-color cadence which observes each sky location three times per night instead of twice. Each simulated dataset includes a spectroscopically-confirmed training set, which we augment to be representative of the test set as part of the classification pipeline. Then, we use the photometric transient classification library snmachine to build classifiers. We find that the active region of the rolling cadence used in the baseline observing strategy yields a 25% improvement in classification performance relative to the background region. This improvement in performance in the actively-rolling region is also associated with an increase of up to a factor of 2.7 in the number of cosmologically-useful Type Ia supernovae relative to the background region. However, adding a third visit per night as implemented in presto-color degrades classification performance due to more irregularly sampled light curves. Overall, our results establish desiderata on the observing cadence related to classification of full SNe light curves, which in turn impacts photometric SNe cosmology with LSST.

In this Letter, we show that the nonthermal nature of dark matter freeze-in production leads to large, totally correlated dark matter-photon isocurvature perturbations, which are imprinted in anisotropies of the cosmic microwave background (CMB). Isocurvature is typically expected from inflationary physics, but the isocurvature from freeze-in arises post inflation. We compute the freeze-in of millicharged dark matter, generated from electron-positron annihilations in the early Universe. We find that current CMB observations from \textit{Planck} exclude this scenario for dark matter masses between 1 MeV and 10 GeV at more than $2\sigma$, whereas upcoming CMB experiments will have the sensitivity to reach at least the $4\sigma$ level. We anticipate any scenario in which dark matter is nonthermally produced to generically give rise to isocurvature. Our work opens a new avenue for exploring fundamental dark matter physics through its impact on cosmological observables.

We present the results of a coordinated campaign to simultaneously observe the M star binary Ross 733 simultaneously in the optical and near-ultraviolet (NUV) with TESS and Swift respectively. We observed two flares in the Swift NUV light curve. One of these was decay phase of a flare that was also detected with TESS and the other was only detected in the NUV. We used the TESS light curve to measure the white-light flare rate of Ross 733, and calculate that the system flares with an energy of $10^{33}$ erg once every 1.5 days. We used our simultaneous observations to measure a pseudo-continuum temperature of $7340^{+810}_{-900}$K during the flare decay. We also used our observations to test the NUV predictions of the 9000 K blackbody flare model, and find that it underestimates number of flares we detect in our Swift NUV light curve. We discuss the reasons for this and attribute it to the unaccounted contributions from emission lines and continuum temperatures above 9000 K. We discuss how additional observations are required to break the degeneracy between the two in future multi-wavelength flare campaigns.

The merger of two neutron stars (NSs) produces the emission of gravitational waves, the formation of a compact object surrounded by a dense and magnetized environment, and the launching of a collimated and relativistic jet, which will eventually produce a short gamma-ray burst (SGRB). The interaction of the jet with the environment has been shown to play a major role in shaping the structure of the outflow that eventually powers the gamma-ray emission. In this paper, we present a set of 2.5 dimensional RMHD simulations that follow the evolution of a relativistic non-magnetized jet through a medium with different magnetization levels, as produced after the merger of two NSs. We find that the predominant consequence of a magnetized ambient medium is that of suppressing instabilities within the jet, and preventing the formation of a series of collimation shocks. One implication of this is that internal shocks lose efficiency, causing bursts with low-luminosity prompt emission. On the other hand, the jet-head velocity and the induced magnetization within the jet are fairly independent from the magnetization of the ambient medium. Future numerical studies with a larger domain are necessary to obtain light curves and spectra in order to better understand the role of magnetized media.

Given their extremely faint apparent brightness, the nature of the first galaxies and how they reionized the Universe's gas are not yet understood. Here we report the discovery, in James Webb Space Telescope (JWST) imaging, of a highly magnified, low-mass (log(M_*/M_sol)=7.70^{+0.11}_{-0.09}) galaxy visible when the Universe was only 510 Myr old, and follow-up prism spectroscopy of the galaxy extending from Lyman alpha to [O III] 5007 in its rest frame. Our JWST spectrum provides [O III] 5007 and H beta detections with a respective signal-to-noise ratio (S/N) of 33 and 7, as well as four additional lines with S/N > 2. These emission lines yield a redshift of z=9.51 and star-formation rate of 0.26^{+0.07}_{-0.05} M_sol per year. The galaxy's large inferred value of [O III]/[O II] = 16 +- 6 suggests that this galaxy has an escape fraction of ionizing radiation larger than 10%, indicating that a population of similar objects could contribute substantially to the reionization budget. Using multiple techniques, we infer a gas oxygen abundance of 12 + log(O/H) = 7.47 +- 0.10 dex, consistent within 2 sigma of the mass-metallicity relation observed for dwarf galaxies in the local Universe.

The three-dimensional characterization of magnetic flux-ropes observed in the heliosphere has been a challenging task for decades. This is mainly due to the limitation to infer the 3D global topology and the physical properties from the 1D time series from any spacecraft. To advance our understanding of magnetic flux-ropes whose configuration departs from the typical stiff geometries, here we present the analytical solution for a 3D flux-rope model with an arbitrary cross-section and a toroidal global shape. This constitutes the next level of complexity following the elliptic-cylindrical (EC) geometry. The mathematical framework was established by Nieves-Chinchilla et al. (2018) ApJ, with the EC flux-rope model that describes the magnetic topology with elliptical cross-section as a first approach to changes in the cross-section. In the distorted-toroidal flux rope model, the cross-section is described by a general function. The model is completely described by a non-orthogonal geometry and the Maxwell equations can be consistently solved to obtain the magnetic field and relevant physical quantities. As a proof of concept, this model is generalized in terms of the radial dependence of current density components. The last part of this paper is dedicated to a specific function, $F(\varphi)=\delta(1-\lambda\cos\varphi)$, to illustrate possibilities of the model. This model paves the way to investigate complex distortions of the magnetic structures in the solar wind. Future investigations will in-depth explore these distortions by analyzing specific events, the implications in the physical quantities, such as magnetic fluxes, heliciy or energy, and evaluating the force balance with the ambient solar wind that allows such distortions.

To pinpoint the peak location of the synchrotron total intensity emission in a spiral arm, we use a map of the spiralarm locations (from the observed arm tangent). Thus In a typical spiral arm in Galactic Quadrant I, we find the peak of the synchrotron radiation to be located about 220 +/-40 pc away from the inner arm edge (hot dust lane) inside the spiral arm. While most of the galactic disk has a clockwise largescale magnetic field, we make a statistical analysis to delimitate more precisely the smaller reverse annulus wiith a counterclockwise galactic magnetic field. We find an annulus width of 2.1 +/-0.3 kpc (measured along the Galactic radius), located from 5.5 to 7.6 kpc from the Galactic Center). The annulus does not overlay with a single spiral arm -- it encompasses segments of two different spiral arms. Using a recent delineation of the position of spiral arms, the field-reversed annulus is seen to encompass the Crux-Centaurus arm (in Galactic Quadrant IV) and the Sagittarius arm (in Galactic Quadrant I).

The motivation to study experimentally CO ice under mimicked interstellar conditions is supported by the large CO gas abundances and ubiquitous presence of CO in icy grain mantles. Upon irradiation in its pure ice form, this highly stable species presents a limited ion and photon-induced chemistry, and an efficient non-thermal desorption. Using infrared spectroscopy, single laser interference, and quadrupole mass spectrometry during CO ice deposition, the CO ice density was estimated as a function of deposition temperature. Only minor variations in the density were found. The proposed methodology can be used to obtain the density of other ice components at various deposition temperatures provided that this value of the density is known for one of these temperatures, which is typically the temperature corresponding to the crystalline form. The apparent tendency of the CO ice density to decrease at deposition temperatures below 14 K is in line with recently published colorimetric measurements. This work allowed to revisit the value of the infrared band strength needed for calculation of the CO ice column density in infrared observations, $8.7 \times 10^{-18} ~ {\rm cm ~ molecule}^{-1}$ at 20 K deposition temperature.

Multiply imaged time-variable sources can be used to measure absolute distances as a function of redshifts and thus determine cosmological parameters, chiefly the Hubble Constant H$_0$. In the two decades up to 2020, through a number of observational and conceptual breakthroughs, this so-called time-delay cosmography has reached a precision sufficient to be an important independent voice in the current ``Hubble tension'' debate between early- and late-universe determinations of H$_0$. The 2020s promise to deliver major advances in time-delay cosmography, owing to the large number of lenses to be discovered by new and upcoming surveys and the vastly improved capabilities for follow-up and analysis. In this review -- after a brief summary of the foundations of the method and recent advances -- we outline the opportunities for the decade and the challenges that will need to be overcome in order to meet the goal of the determination of H$_0$ from time-delay cosmography with 1\% precision and accuracy.

The Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect plays an important role in the rotational properties and evolution of asteroids. While the YORP effect induced by the macroscopic shape of the asteroid and by the presence of surface boulders has been well studied, no investigation has been performed yet regarding how craters with given properties influence this effect. We introduce and estimate the crater-induced YORP effect (CYORP), which arises from the concave structure of the crater, to investigate the magnitude of the resulting torques as a function of varying properties of the crater and the asteroid by a semi-analytical method. The CYORP torque can be comparable to the normal YORP torque when the size of the crater is about one-tenth of the size of the asteroid, or equivalently when the crater/roughness covers one-tenth of the asteroid's surface. Although the torque decreases with the crater size, the combined contribution of all small craters can become non-negligible due to their large number when the commonly used power-law crater size distribution is considered. The CYORP torque of small concave structures, usually considered as surface roughness, is essential to the accurate calculation of the complete YORP torque. Under the CYORP effect that is produced by collisions, asteroids go through a random walk in spin rate and obliquity, with a YORP reset timescale typically of 0.4 Myr. This has strong implications for the rotational evolution and orbital evolution of asteroids. Craters and roughness on asteroid surfaces can influence the YORP torques and therefore the rotational properties and evolution of asteroids. We suggest that the CYORP effect should be considered in the future investigation of the YORP effect on asteroids.

Transit and radial-velocity surveys over the past two decades have uncovered a significant population of short-period exoplanets. Among them are hot Jupiters, which are gas giant planets with orbital periods of a few days and found in 0.1-1% of Sun-like stars. Hot Jupiters are expected to be engulfed during their host star's radial expansion on the red giant branch. Planetary engulfment has been studied extensively as it may account for observed rapidly rotating and chemically enriched giant stars. We perform 3D hydrodynamical simulations of hot Jupiter engulfment by a 1 solar mass, 4 solar radii early red giant. Our "global" simulations simultaneously resolve the stellar envelope and planetary structure, modelling the hot Jupiter as a polytropic gas sphere. We find that approximately 90% of the hot Jupiter's mass is ablated in the convective part of the giant envelope, which would enhance the surface lithium abundance by 0.1 dex. The hot Jupiter is disrupted by a combination of ram pressure and tidal forces near the base of the convective envelope, with the deepest material penetrating to the radiative zone. The star experiences modest spin-up (~1 km/s), although engulfing a more massive companion could produce a rapidly rotating giant. Drag heating near the surface could exceed the unperturbed stellar luminosity and power an optical transient. For the amount of unbound ejecta recorded in the simulation, H-recombination could also power a transient that is around ten times the pre-engulfment luminosity, for several days.

Gamma-ray bursts (GRBs) have been considered as potential very high-energy photon emitters due to the large amount of energy released as well as the strong magnetic fields involved. However, the detection of TeV photons is not expected from bursts beyond a redshift of $z\gtrsim 0.1$ due to the attenuation caused by the interaction with photons from the extragalactic background light (EBL). The observation of an 18-TeV photon from GRB 221009A (z=0.15) last October 9th, 2022 has challenged what we know about the TeV-emission mechanisms and the extragalactic background. Recent explanations exploiting candidates of dark matter have started to appear. In this paper, we discuss the required conditions and limitations within the most plausible scenario, synchrotron-self Compton (SSC) radiation in the GRB afterglow, to interpret the one 18-TeV photon observation besides the EBL. To avoid the Klein-Nishina effect, we find an improbable value of the microphysical magnetic parameter below $10^{-6}$ for a circumburst medium $>{\rm cm^{-3}}$ (expected in the collapsar scenario). Therefore, we explore possible scenarios in terms of ALPs and dark photon mechanisms to interpret this highly-energetic photon and discuss the implications in the GRB energetics. We find the ALPs and dark photon scenarios explain the 18 TeV photon but not the 251 TeV photon.

Using multifrequency observations, from radio to gamma-rays of the Blazar Mrk 501, we constructed their corresponding light curves, and built its periodograms using RobPer and Lomb-Scargle algorithms. Long-term variability was also studied using the Power Density Spectrum, the Detrended Function Analysis using the software VARTOOLS. The result of these techniques showed an achromatic periodicity <~ 229 d. This combined with the result of a pink colour noise in the spectra lead us to propose that the periodicity was produced by a secondary eclipsing supermassive binary black hole orbiting the primary one locked inside the central engine of Mrk 501. We built a relativistic eclipsing model of this phenomenon using Jacobi elliptical functions finding a periodic relativistic eclipse occurring every ~225 d in all the studied wavebands.

The origin of Ultra High Energy Cosmic Rays is a 60-year old mystery. We show that with more events at the highest energies (above 150~EeV) it may be possible to limit the character of the sources and learn about the intervening magnetic fields. Individual sources become more prominent, relative to the background, as the horizon diminishes. An event-by-event, composition-dependent observatory would allow a ``tomography'' of the sources as different mass and energy groups probe different GZK horizons. A major goal here is to provide a methodology to distinguish between steady and transient or highly variable sources. Using recent Galactic magnetic field models, we calculate ``treasure'' sky maps to identify the most promising directions for detecting Extreme Energy Cosmic Rays (EECR) doublets, events that are close in arrival time and direction. On this basis, we predict the incidence of doublets as a function of the nature of the source host galaxy. Based on the asymmetry in the distribution of time delays, we show that observation of doublets might distinguish source models. In particular the Telescope Array hotspot could exhibit temporal variability as it is in a ``magnetic window'' of small time delays. These considerations could improve the use of data with existing facilities and the planning of future ones such as Global Cosmic Ray Observatory - GCOS.

The Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo Interferometer Collaborations have now detected all three classes of compact binary mergers: binary black hole (BBH), binary neutron star (BNS), and neutron star-black hole (NSBH). For coalescences involving neutron stars, the simultaneous observation of gravitational and electromagnetic radiation produced by an event, has broader potential to enhance our understanding of these events, and also to probe the equation of state (EOS) of dense matter. However, electromagnetic follow-up to gravitational wave (GW) events requires rapid real-time detection and classification of GW signals, and conventional detection approaches are computationally prohibitive for the anticipated rate of detection of next-generation GW detectors. In this work, we present the first deep learning based results of classification of GW signals from NSBH mergers in \textit{real} LIGO data. We show for the first time that a deep neural network can successfully distinguish all three classes of compact binary mergers and separate them from detector noise. Specifically, we train a convolutional neural network (CNN) on $\sim 500,000$ data samples of real LIGO noise with injected BBH, BNS, and NSBH GW signals, and we show that our network has high sensitivity and accuracy. Most importantly, we successfully recover the two confirmed NSBH events to-date (GW200105 and GW200115) and the two confirmed BNS mergers to-date (GW170817 and GW190425), together with $\approx 90\%$ of all BBH candidate events from the third Gravitational Wave Transient Catalog, GWTC-3. These results are an important step towards low-latency real-time GW detection, enabling multi-messenger astronomy.

A numerical detection of the radius-dependent spin transition of dark matter halos is reported. Analyzing the data from the IllustrisTNG simulations, we measure the halo spin vectors at several inner radii within the virial boundaries and investigate their orientations in the principal frames of the tidal and velocity shear fields, called the Tweb and Vweb, respectively. The halo spin vectors in the high-mass section exhibit a transition from the Tweb intermediate to major principal axes as they are measured at more inner radii, which holds for both of the dark matter and baryonic components. The radius threshold at which the transition occurs depends on the smoothing scale, $R_{f}$, becoming larger as $R_{f}$ decreases. For the case of the Vweb, the occurrence of the radius-dependent spin transition is witnessed only when $R_{f}\ge 1\, h^{-1}$Mpc. Repeating the same analysis but with the vorticity vectors, we reveal a critical difference from the spins. The vorticity vectors are always perpendicular to the Tweb (Vweb) major principal axes, regardless of $R_{f}$, which indicates that the halo inner spins are not strongly affected by the generation of vorticity. It is also shown that the halo spins, as well as the Tweb (Vweb) principal axes, have more directional coherence over a wide range of radial distances in the regions where the vorticity vectors have higher magnitudes. The physical interpretations and implications of our results are discussed.

As the number of confirmed exoplanets continues to grow, there is an increased push to spectrally characterize them to determine their atmospheric composition, formation paths, rotation rates, and habitability. However, there is a large population of known exoplanets that either do not transit their star or have been detected via the radial velocity (RV) method at very small angular separations such that they are inaccessible to traditional coronagraph systems. Vortex Fiber Nulling (VFN) is a new single-aperture interferometric technique that uses the entire telescope pupil to bridge the gap between traditional coronagraphy and RV or Transit methods by enabling the direct observation and spectral characterization of targets at and within the diffraction limit. By combining a vortex mask with a single mode fiber, the on-axis starlight is rejected while the off-axis planet light is coupled and efficiently routed to a radiometer or spectrograph for analysis. We have demonstrated VFN in the lab monochromatically in the past. In this paper we present a polychromatic validation of VFN with nulls of $<10^{-4}$ across 15% bandwidth light. We also provide an update on deployment plans and predicted yield estimates for the VFN mode of the Keck Planet Imager and Characterizer (KPIC) instrument. Using PSISIM, a simulation package developed in cooperation with several groups, we assess KPIC VFN's ability to detect and characterize different types of targets including planet candidates around promising young-moving-group stars as well as known exoplanets detected via the RV method. The KPIC VFN on-sky demonstration will pave the road to deployment on future instruments such as Keck-HISPEC and TMT-MODHIS where it could provide high-resolution spectra of sub-Jupiter mass planets down to 5 milliarcseconds from their star.

The Keck Planet Imager and Characterizer (KPIC) is a series of upgrades for the Keck II Adaptive Optics (AO) system and the NIRSPEC spectrograph to enable diffraction-limited, high-resolution ($R>30,000$) spectroscopy of exoplanets and low-mass companions in the K and L bands. Phase I consisted of single-mode fiber injection/extraction units (FIU/FEU) used in conjunction with an H-band pyramid wavefront sensor. Phase II, deployed and commissioned in 2022, adds a 1000-actuator deformable mirror, beam-shaping optics, a vortex coronagraph, and other upgrades to the FIU/FEU. The use of single-mode fibers provides a gain in stellar rejection, a substantial reduction in sky background, and an extremely stable line-spread function on the spectrograph. In this paper we present the results of extensive system-level laboratory testing and characterization showing the instrument's Phase II throughput, stability, repeatability, and other key performance metrics prior to delivery and during installation at Keck. We also demonstrate the capabilities of the various observing modes enabled by the new system modules using internal test light sources. Finally, we show preliminary results of on-sky tests performed in the first few months of Phase II commissioning along with the next steps for the instrument. Once commissioning of Phase II is complete, KPIC will continue to characterize exoplanets at an unprecedented spectral resolution, thereby growing its already successful track record of 23 detected exoplanets and brown dwarfs from Phase I. Using the new vortex fiber nulling (VFN) mode, Phase II will also be able to search for exoplanets at small angular separations less than 45 milliarcseconds which conventional coronagraphs cannot reach.

We have obtained accurate dust reddening from far-ultraviolet (UV) to the mid-infrared (IR) for up to 5 million stars by the star-pair algorithm based on LAMOST stellar parameters along with GALEX, Pan-STARRS 1, Gaia, SDSS, 2MASS, and WISE photometric data. The typical errors are between 0.01-0.03 mag for most colors. We derived the empirical reddening coefficients for 21 colors both in the traditional (single valued) way and as a function of Teff and E(B-V) by using the largest samples of accurate reddening measurements, together with the extinction values from Schlegel et al. The corresponding extinction coefficients have also been obtained. The results are compared with model predictions and generally in good agreement. Comparisons with measurements in the literature show that the Teff- and E(B-V)-dependent coefficients explain the discrepancies between different measurements naturally, i.e., using sample stars of different temperatures and reddening. Our coefficients are mostly valid in the extinction range of 0-0.5 mag and the temperature range of 4000-10000 K. We recommend that the new Teff and E(B-V) dependent reddening and extinction coefficients should be used in the future. A Python package is also provided for the usage of the coefficient

In our previous work [Tamburini and Licata (2017)] we discussed the hypothesis that the ultrafast periodic spectral modulations with frequency $f_S \simeq 0.61$ THz found by Borra and Trottier (2016) in $236$ main sequence stars from a sample of $2.5$ million spectra of galactic halo stars of the Sloan Digital Sky Survey were due to axion-like dark matter piled up in the center of these stars. These temporary matter/dark matter structures are characterized by a spacetime geometry rapidly oscillating at frequencies that depend on the axion mass $m_a$ [Brito {\it et~al.} (2015); Brito {\it et~al.} (2016)]. Borra (2013) found two additional frequencies, $ f_{1,G} \simeq 9.5$ THz and $f_{2,G} \simeq 8.9$ THz, in the Sloan dataset of galaxies, redshifted by the cosmological expansion and, for any redshift value, $f_S + f_{2,G} = f_{1,G}$ is found. Hippke (2019) showed that $f_{2,G}$ is spurious and introduced by the data analysis procedure due to the nonrandom spacing of the spectral absorption lines. This was not even found by Isaacson (2019) when re-observed four of these stars with different instrumentation and data reduction procedure. Interestingly, Hippke found $f_S$ in the solar spectrum but not in the Kurucz (2005) artificial solar spectrum whilst its spectral power estimated by Isaacson resulted below the accepted error $(1\%)$. From these results, we discuss the validity of our ansatz by analyzing the common features present in all the spectra. In the worst case, if all the three frequencies are not real the oscillating axion core models is not valid. Assuming, instead, the validity of $f_S$ from the results from the analysis of the solar spectra, those oscillating modes may be transient modes favoring the axion hypothesis in the mass range $(10^{- 3} < m_a < 2.4 \times 10^{3})~ \mathrm{\mu eV}$, also according to the recent limits from the gamma ray burst GRB221009A.

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

We perform a homogeneous analysis of an unprecedented set of spatially resolved scaling relations (SRs) between ISM components and other properties in the range of scales 0.3-3.4 kpc. We also study some ratios: dust-to-stellar, dust-to-gas, and dust-to-metal. We use a sample of 18 large, spiral, face-on DustPedia galaxies. All the SRs are moderate/strong correlations except the dust-HI SR that does not exist or is weak for most galaxies. The SRs do not have a universal form but each galaxy is characterized by distinct correlations, affected by local processes and galaxy peculiarities. The SRs hold starting from 0.3 kpc, and if a breaking down scale exists it is < 0.3 kpc. By evaluating all galaxies at 3.4 kpc, differences due to peculiarities of individual galaxies are cancelled out and the corresponding SRs are consistent with those of whole galaxies. By comparing subgalactic and global scales, the most striking result emerges from the SRs involving ISM components: the dust-total gas SR is a good correlation at all scales, while the dust-H2 and dust-HI SRs are good correlations at subkpc/kpc and total scales, respectively. For the other explored SRs, there is a good agreement between small and global scales and this may support the picture where the main physical processes regulating the properties and evolution of galaxies occur locally. Our results are consistent with the hypothesis of self-regulation of the SF process. The analysis of subgalactic ratios shows that they are consistent with those derived for whole galaxies, from low to high z, supporting the idea that also these ratios could be set by local processes. Our results highlight the heterogeneity of galaxy properties and the importance of resolved studies on local galaxies in the context of galaxy evolution. They also provide observational constraints to theoretical models and updated references for high-z studies.

Simulation of the cosmic clustering of massive neutrinos is a daunting task, due both to their large velocity dispersion and to their weak clustering power becoming swamped by Poisson shot noise. We present a new approach, the multi-fluid hybrid-neutrino simulation, which partitions the neutrino population into multiple flows, each of which is characterised by its initial momentum and treated as a separate fluid. These fluid flows respond initially linearly to nonlinear perturbations in the cold matter, but slowest flows are later converted to a particle realisation should their clustering power exceed some threshold. After outlining the multi-fluid description of neutrinos, we study the conversion of the individual flows into particles, in order to quantify transient errors, as well as to determine a set of criteria for particle conversion. Assembling our results into a total neutrino power spectrum, we demonstrate that our multi-fluid hybrid-neutrino simulation is convergent to $<3\%$ if conversion happens at $z=19$ and agrees with more expensive simulations in the literature for neutrino fractions as high as $\Omega_\nu h^2 = 0.005$. Moreover, our hybrid-neutrino approach retains fine-grained information about the neutrinos' momentum distribution. However, the momentum resolution is currently limited by free-streaming transients excited by missing information in the neutrino particle initialisation procedure, which restricts the particle conversion to z $\gtrsim 19$ if percent-level resolution is desired.

With MaNGA integral field spectroscopy, we present a resolved analysis of star formation for 29 jellyfish galaxies in nearby clusters, identified from radio continuum imaging taken by the Low Frequency Array. Simulations predict enhanced star formation on the "leading half" of galaxies undergoing ram pressure stripping, and in this work we report observational evidence for this elevated star formation. The dividing line (through the galaxy center) that maximizes this star formation enhancement is systematically tied to the observed direction of the ram pressure stripped tail, suggesting a physical connection between ram pressure and this star formation enhancement. We also present a case study on the distribution of molecular gas in one jellyfish galaxy from our sample, IC3949, using ALMA CO J=1-0, HCN J=1-0, and HCO$^+$ J=1-0 observations from the ALMaQUEST survey. The $\mathrm{H_2}$ depletion time (as traced by CO) in IC3949 ranges from $\sim\!1\,\mathrm{Gyr}$ in the outskirts of the molecular gas disk to $\sim\!11\,\mathrm{Gyr}$ near the galaxy center. IC3949 shows a clear region of enhanced star formation on the leading half of the galaxy where the average depletion time is $\sim\!2.7\,\mathrm{Gyr}$, in line with the median value for the galaxy on the whole. Dense gas tracers, HCN and HCO$^+$, are only detected at the galaxy center and on the leading half of IC3949. Our results favour a scenario in which ram pressure compresses the interstellar medium, promoting the formation of molecular gas that in turn fuels a localized increase of star formation.

Non-linear cosmic structures contain valuable information on the expansion history of the background space-time, the nature of dark matter, and the gravitational interaction. The recently developed kinetic field theory of cosmic structure formation (KFT) allows to accurately calculate the non-linear power spectrum of cosmic density fluctuations up to wave numbers of $k\lesssim10\,h\,\mathrm{Mpc}^{-1}$ at redshift zero. Cosmology and gravity enter this calculation via two functions, viz. the background expansion function and possibly a time-dependent modification of the gravitational coupling strength. The success of the cosmological standard model based on general relativity suggests that cosmological models in generalized theories of gravity should have observable effects differing only weakly from those in standard cosmology. Based on this assumption, we derive the functional, first-order Taylor expansion of the non-linear power spectrum of cosmic density fluctuations obtained from the mean-field approximation in KFT in terms of the expansion function and the gravitational coupling strength. This allows us to study non-linear power spectra expected in large classes of generalized gravity theories. To give one example, we apply our formalism to generalized Proca theories.

Velocity dispersion of the massive neutrinos presents a daunting challenge for non-linear cosmological perturbation theory. We consider the neutrino population as a collection of non-linear fluids, each with uniform initial momentum, through an extension of the Time Renormalization Group perturbation theory. Employing recently-developed Fast Fourier Transform techniques, we accelerate our non-linear perturbation theory by more than two orders of magnitude, making it quick enough for practical use. After verifying that the neutrino mode-coupling integrals and power spectra converge, we show that our perturbation theory agrees with N-body neutrino simulations to within $10\%$ for neutrino fractions $\Omega_{\nu,0} h^2 \leq 0.005$ up to wave numbers of $k = 1~h/$Mpc, an accuracy consistent with $\leq 2.5\%$ errors in the neutrino mass determination. Non-linear growth represents a $>10\%$ correction to the neutrino power spectrum even for density fractions as low as $\Omega_{\nu,0} h^2 = 0.001$, demonstrating the limits of linear theory for accurate neutrino power spectrum predictions.

We propose a new approach for large-scale high-dynamic range computational imaging. Deep Neural Networks (DNNs) trained end-to-end can solve linear inverse imaging problems almost instantaneously. While unfolded architectures provide necessary robustness to variations of the measurement setting, embedding large-scale measurement operators in DNN architectures is impractical. Alternative Plug-and-Play (PnP) approaches, where the denoising DNNs are blind to the measurement setting, have proven effective to address scalability and high-dynamic range challenges, but rely on highly iterative algorithms. We propose a residual DNN series approach, where the reconstructed image is built as a sum of residual images progressively increasing the dynamic range, and estimated iteratively by DNNs taking the back-projected data residual of the previous iteration as input. We demonstrate on simulations for radio-astronomical imaging that a series of only few terms provides a high-dynamic range reconstruction of similar quality to state-of-the-art PnP approaches, at a fraction of the cost.

A tidal interaction between a star and a close-in exoplanet leads to shrinkage of the planetary orbit and eventual tidal disruption of the planet. Measuring the shrinkage of the orbits will allow for the tidal quality parameter of the star ($Q'_\star$) to be measured, which is an important parameter to obtain information about stellar interiors. We analyse data from TESS for two targets known to host close-in hot Jupiters, WASP-18 and WASP-19, to measure the current limits on orbital period variation and provide new constrains on $Q'_\star$. We modelled the transit shape using all the available TESS observations and fitted the individual transit times of each transit. We used previously published transit times together with our results to fit two models, a constant period model, and a quadratic orbital decay model, MCMC algorithms. We find period change rates of $(-0.11\pm0.21)\times10^{-10}$ for WASP-18b and $(-0.35\pm0.22)\times10^{-10}$ for WASP-19b and we do not find significant evidence of orbital decay in these targets. We obtain new lower limits for $Q'_\star$ of $(1.42\pm0.34)\times10^7$ in WASP-18 and $(1.26\pm0.10)\times10^6$ in WASP-19, corresponding to upper limits of the orbital decay rate of $-0.45\times10^{-10}$ and $-0.71\times10^{-10}$, respectively, with a 95% confidence level. We compare our results with other relevant targets for tidal decay studies. We find that the orbital decay rate in both WASP-18b and WASP-19b appears to be smaller than the measured orbital decay of WASP-12b. We show that the minimum value of $Q'_\star$ in WASP-18 is two orders of magnitude higher than that of WASP-12, while WASP-19 has a minimum value one order of magnitude higher, which is consistent with other similar targets. Further observations are required to constrain the orbital decay of WASP-18 and WASP-19.

Aims. We present an implementation of an algorithm for 3D time-dependent Monte Carlo radiative transfer. It allows one to simulate temperature distributions as well as images and spectral energy distributions of the scattered light and thermal reemission radiation for variable illuminating and heating sources embedded in dust distributions, such as circumstellar disks and dust shells on time scales up to weeks. Methods. We extended the publicly available 3D Monte Carlo radiative transfer code POLARIS with efficient methods for the simulation of temperature distributions, scattering, and thermal reemission of dust distributions illuminated by temporally variable radiation sources. The influence of the chosen temporal step width and the number of photon packages per time step as key parameters for a given configuration is shown by simulating the temperature distribution in a spherical envelope around an embedded central star. The effect of the optical depth on the temperature simulation is discussed for the spherical envelope as well as for a model of a circumstellar disk with an embedded star. Finally, we present simulations of an outburst of a star surrounded by a circumstellar disk. Results. The presented algorithm for time-dependent 3D continuum Monte Carlo radiative transfer is a valuable basis for preparatory studies as well as for the analysis of continuum observations of the dusty environment around variable sources, such as accreting young stellar objects. In particular, the combined study of light echos in the optical and near-infrared wavelength range and the corresponding time-dependent thermal reemission observables of variable, for example outbursting sources, becomes possible on all involved spatial scales.

Primordial black holes (PBHs) can be not only cold dark matter candidates but also progenitors of binary black holes observed by LIGO-Virgo-KAGRA (LVK) Collaboration. The PBH mass can be shifted to the heavy distribution if multi-merger processes occur. In this work, we constrain the merger history of PBH binaries using the gravitational wave events from the third Gravitational-Wave Transient Catalog (GWTC-3). Considering four commonly used PBH mass functions, namely the log-normal, power-law, broken power-law, and critical collapse forms, we find that the multi-merger processes make a subdominant contribution to the total merger rate. Therefore, the effect of merger history can be safely ignored when estimating the merger rate of PBH binaries. We also find that GWTC-3 is best fitted by the log-normal form among the four PBH mass functions and confirm that the stellar-mass PBHs cannot dominate cold dark matter.

High-mass gamma-ray binaries are powerful nonthermal galactic sources, some of them hosting a pulsar whose relativistic wind interacts with a likely inhomogeneous stellar wind. So far, modeling these sources including stellar wind inhomogeneities has been done using either simple analytical approaches or heavy numerical simulations, none of which allow for an exploration of the parameter space that is both reasonably realistic and general. Applying different semi-analytical tools together, we study the dynamics and high-energy radiation of a pulsar wind colliding with a stellar wind with different degrees of inhomogeneity to assess the related observable effects. We computed the arrival of clumps to the pulsar wind-stellar wind interaction structure using a Monte Carlo method and a phenomenological clumpy-wind model. The dynamics of the clumps that reach deep into the pulsar wind zone was computed using a semi-analytical approach. This approach allows for the characterization of the evolution of the shocked pulsar wind region in times much shorter than the orbital period. With this three-dimensional information about the emitter, we applied analytical adiabatic and radiative models to compute the variable high-energy emission produced on binary scales. An inhomogeneous stellar wind induces stochastic hour-timescale variations in the geometry of the two-wind interaction structure on binary scales. Depending on the degree of stellar wind inhomogeneity, 10-100% level hour-scale variability in the X-rays and gamma rays is predicted, with the largest variations occurring roughly once per orbit. Our results, based on a comprehensive approach, show that present X-ray and future very-high-energy instrumentation can allow us to trace the impact of a clumpy stellar wind on the shocked pulsar wind emission in a gamma-ray binary.

Galaxy populations are known to exhibit a strong colour bimodality, corresponding to blue star-forming and red quiescent subpopulations. The relative abundance of the two populations has been found to vary with stellar mass and environment. In this paper, we explore the effect of environment considering different types of measurements. We choose a sample of 49, 911 galaxie with $0.05 < z < 0.18$ from the Galaxy And Mass Assembly survey. We study the dependence of the fraction of red galaxies on different measures of the local environment as well as the large-scale "geometric" environment arXiv:1412.2141 [astro-ph.CO] defined by density gradients in the surrounding cosmic web. We find that the red galaxy fraction varies with the environment independently of stellar mass. The local environmental measures exhibit a larger variation in red fraction from lower to higher density, as compared to the geometric environment. By comparing the different environmental densities pairwise, we show that no density measurement fully explains the observed environmental red fraction variation, suggesting the different densities contain different information. We test whether the local environmental measures, when combined together, can explain all the observed environmental red fraction variation. The geometric environment has a small residual effect, and this effect is larger for voids than any other type of geometric environment. This could provide a test of the physics applied to cosmological-scale galaxy evolution simulations as it combines large-scale effects with local environmental impact.

The discovery of radio pulsars (PSRs) around the supermassive black hole (SMBH) in our Galactic Center (GC), Sagittarius A* (Sgr A*), will have significant implications for tests of gravity. In this paper, we predict restrictions on the parameters of the Yukawa gravity by timing a pulsar around Sgr A* with a variety of orbital parameters. Based on a realistic timing accuracy of the times of arrival (TOAs), $\sigma_{\rm TOA}=100\,\mu{\rm s}$, and using a number of 960 TOAs in a 20-yr observation, our numerical simulations show that the PSR-SMBH system will improve current tests of the Yukawa gravity when the range of the Yukawa interaction varies between $10^{1}$-$10^{4}\,{\rm AU}$, and it can limit the graviton mass to be $m_g \lesssim 10^{-24}\,{\rm eV}/c^2$.

Red supergiant stars (RSGs, Minit = 10-40Msun) are known to eject large amounts of material, as much as half of their initial mass during this evolutionary phase. However, the processes powering the mass ejection in low- and intermediate-mass stars do not work for RSGs and the mechanism that drives the ejection remains unknown. Different mechanisms have been proposed as responsible for this mass ejection but so far little is known about the actual processes taking place in these objects. Here we present high angular resolution interferometric ALMA maps of VY CMa continuum and molecular emission, which resolve the structure of the ejecta with unprecedented detail. The study of the molecular emission from the ejecta around evolved stars has been shown to be an essential tool in determining the characteristics of the mass loss ejections. Our aim is thus to use the information provided by these observations to understand the ejections undergone by VY CMa and to determine their possible origins. We inspected the kinematics of molecular emission observed. We obtained position-velocity diagrams and reconstructed the 3D structure of the gas traced by the different species. It allowed us to study the morphology and kinematics of the gas traced by the different species surrounding VY CMa. Two types of ejecta are clearly observed: extended, irregular, and vast ejecta surrounding the star that are carved by localized fast outflows. The structure of the outflows is found to be particularly flat. We present a 3D reconstruction of these outflows and proof of the carving. This indicates that two different mass loss processes take place in this massive star. We tentatively propose the physical cause for the formation of both types of structures. These results provide essential information on the mass loss processes of RSGs and thus of their further evolution.

The presence of debris in Earth's orbit poses a significant risk to human activity in outer space. This debris population continues to grow due to ground launches, loss of external parts from space ships, and uncontrollable collisions between objects. A computationally feasible continuum model for the growth of the debris population and its spatial distribution is therefore critical. Here we propose a diffusion-collision model for the evolution of debris density in Low-Earth Orbit (LEO) and its dependence on ground-launch policy. We parametrize this model and test it against data from publicly available object catalogs to examine timescales for uncontrolled growth. Finally, we consider sensible launch policies and cleanup strategies and how they reduce the future risk of collisions with active satellites or space ships.

A radio wave absorber induces technology innovation in the fields of astrophysics, industrial application for communication, and security. We recently developed and characterized an absorber for millimeter wavelengths. The absorber consists of epoxy, carbon black, and expanded polystyrene beads. The typical diameter of polystyrene beads corresponds to the scale of Mie scattering for the multi-scattering of photons in the absorber. The wide directivity of Mie scattering increases the effective length of the mean free path in the absorber. By applying this effect, we succeeded in improving the performance of the absorber. We characterized the reflectance and transmittance of the absorber in millimeter length and also measured the transmittance in sub-millimeter wavelength.

Weak gravitational lensing by the intervening large-scale structure (LSS) of the Universe is the leading non-linear effect on the anisotropies of the cosmic microwave background (CMB). The integrated line-of-sight mass that causes the distortion -- known as lensing convergence -- can be reconstructed from the lensed temperature and polarization anisotropies via estimators quadratic in the CMB modes, and its power spectrum has been measured from multiple CMB experiments. Sourced by the non-linear evolution of structure, the bispectrum of the lensing convergence provides additional information on late-time cosmological evolution complementary to the power spectrum. However, when trying to estimate the summary statistics of the reconstructed lensing convergence, a number of noise-biases are introduced, as previous studies have shown for the power spectrum. Here, we explore for the first time the noise-biases in measuring the bispectrum of the reconstructed lensing convergence. We compute the leading noise-biases in the flat-sky limit and compare our analytical results against simulations, finding excellent agreement. Our results are critical for future attempts to reconstruct the lensing convergence bispectrum with real CMB data.

As the volume and quality of modern galaxy surveys increase, so does the difficulty of measuring the cosmological signal imprinted in galaxy shapes. Weak gravitational lensing sourced by the most massive structures in the Universe generates a slight shearing of galaxy morphologies called cosmic shear, key probe for cosmological models. Modern techniques of shear estimation based on statistics of ellipticity measurements suffer from the fact that the ellipticity is not a well-defined quantity for arbitrary galaxy light profiles, biasing the shear estimation. We show that a hybrid physical and deep learning Hierarchical Bayesian Model, where a generative model captures the galaxy morphology, enables us to recover an unbiased estimate of the shear on realistic galaxies, thus solving the model bias.

The observability of Lyman-alpha emitting galaxies (LAEs) during the Epoch of Reionization can provide a sensitive probe of the evolving neutral hydrogen gas distribution, thus setting valuable constraints to distinguish different reionization models. In this study, we utilize the new THESAN suite of large-volume (95.5 cMpc) cosmological radiation-hydrodynamic simulations to directly model the Ly$\alpha$ emission from individual galaxies and the subsequent transmission through the intergalactic medium. THESAN combines the AREPO-RT radiation-hydrodynamic solver with the IllustrisTNG galaxy formation model and includes high- and medium-resolution simulations designed to investigate the impacts of halo-mass-dependent escape fractions, alternative dark matter models, and numerical convergence. We find important differences in the Ly$\alpha$ transmission based on reionization history, bubble morphology, frequency offset from line centre, and galaxy brightness. For a given global neutral fraction, Ly$\alpha$ transmission reduces slightly when low mass haloes dominate reionization over high mass haloes. Furthermore, the variation across sightlines for a single galaxy is greater than the variation across all galaxies. This collectively affects the visibility of LAEs, directly impacting observed Ly$\alpha$ luminosity functions (LFs). We employ Gaussian Process Regression using SWIFTEmulator to rapidly constrain an empirical model for dust escape fractions and emergent spectral line profiles to match observed LFs. We find that dust strongly impacts the Ly$\alpha$ transmission and covering fractions of $M_{UV} < -19$ galaxies in $M_{vir} > 10^{11} {\rm M}_{\odot}$ haloes, such that the dominant mode of removing Ly$\alpha$ photons in non-LAEs changes from low IGM transmission to high dust absorption around $z \sim 7$.

We present the discovery of a new optical/X-ray source likely associated with the Fermi $\gamma$-ray source 4FGL J1408.6-2917. Its high-amplitude periodic optical variability, large spectroscopic radial velocity semi-amplitude, evidence for optical emission lines and flaring, and X-ray properties together imply the source is probably a new black widow millisecond pulsar binary. We compile the properties of the 41 confirmed and suspected field black widows, finding a median secondary mass of $0.027\pm0.003\,M_{\odot}$. Considered jointly with the more massive redback millisecond pulsar binaries, we find that the "spider" companion mass distribution remains strongly bimodal, with essentially zero systems having companion masses between $\sim0.07-0.1\,M_{\odot}$. X-ray emission from black widows is typically softer and less luminous than in redbacks, consistent with less efficient particle acceleration in the intrabinary shock in black widows, excepting a few systems that appear to have more efficient "redback-like" shocks. Together black widows and redbacks dominate the census of the fastest-spinning field millisecond pulsars in binaries with known companion types, making up $\gtrsim$80% of systems with $P_{\rm{spin}}<2\,\rm{ms}$. Similar to redbacks, the neutron star masses in black widows appear on average significantly larger than the canonical $1.4\,M_{\odot}$, and many of the highest-mass neutron stars claimed to date are black widows with $M_{\rm{NS}}\gtrsim2.1\,M_{\odot}$. Both of these observations are consistent with an evolutionary picture where spider millisecond pulsars emerge from short orbital period progenitors that had a lengthy period of mass transfer initiated while the companion was on the main sequence, leading to fast spins and high masses.

Early dark energy (EDE) offers a solution to the so-called Hubble tension. Recently, it was shown that the constraints on EDE using Markov Chain Monte Carlo are affected by prior volume effects. The goal of this paper is to present constraints on the fraction of EDE, $f_\mathrm{EDE}$, and the Hubble parameter, $H_0$, which are not subject to prior volume effects. We conduct a frequentist profile likelihood analysis considering Planck cosmic microwave background, BOSS full-shape galaxy clustering, DES weak lensing, and SH0ES supernova data. Contrary to previous findings, we find that $H_0$ for the EDE model is in statistical agreement with the SH0ES direct measurement at $\leq 1.7\,\sigma$ for all data sets. For our baseline data set (Planck + BOSS), we obtain $f_\mathrm{EDE} = 0.087\pm 0.037$ and $H_0 = 70.57 \pm 1.36\, \mathrm{km/s/Mpc}$ at $68\%$ confidence limit. We conclude that EDE is a viable solution to the Hubble tension.

The clustering signals of galaxy clusters are known to be powerful tools for self-calibrating the mass-observable relation and are complementary to cluster abundance and lensing. In this work, we explore the possibility of combining three correlation functions -- cluster lensing, the cluster-galaxy cross-correlation function, and the galaxy auto-correlation function -- to self-calibrate optical cluster selection bias, the boosted clustering and lensing signals in a richness-selected sample mainly caused by projection effects. We develop mock catalogues of redMaGiC-like galaxies and redMaPPer-like clusters by applying Halo Occupation Distribution (HOD) models to N-body simulations and using counts-in-cylinders around massive haloes as a richness proxy. In addition to the previously known small-scale boost in projected correlation functions, we find that the projection effects also significantly boost 3D correlation functions out to scales 100 $h^{-1} \mathrm{Mpc}$. We perform a likelihood analysis assuming survey conditions similar to that of the Dark Energy Survey (DES) and show that the selection bias can be self-consistently constrained at the 10% level. We discuss strategies for applying this approach to real data. We expect that expanding the analysis to smaller scales and using deeper lensing data would further improve the constraints on cluster selection bias.

Thermodynamically induced length fluctuations of high-reflectivity mirror coatings put a fundamental limit on sensitivity and stability of precision optical interferometers like gravitational wave detectors and ultra-stable lasers. The main contribution - Brownian thermal noise - is related to the mechanical loss of the coating material. Owing to their low mechanical losses, Al\textsubscript{0.92}Ga\textsubscript{0.08}As/GaAs crystalline mirror coatings are expected to reduce this limit. At room temperature they have demonstrated lower Brownian thermal noise than with conventional amorphous coatings. %However, no detailed study on the noise constituents from these coatings in optical interferometers has been conducted. We present a detailed study on the spatial and temporal noise properties of such coatings by using them in two independent cryogenic silicon optical Fabry-Perot resonators operated at 4 K, 16 K and 124 K. We confirm the expected low Brownian thermal noise, but also discover two new noise sources that exceed the Brownian noise: birefringent noise that can be canceled via polarization averaging and global excess noise (10 dB above Brownian noise). These new noise contributions are a barrier to improving ultra-stable lasers and the related performance of atomic clocks, and potentially limit the sensitivity of third-generation gravitational wave detectors. Hence, they must be considered carefully in precision interferometry experiments using similar coatings based on semiconductor materials.

The accumulation of certain types of dark matter particles in neutron star cores due to accretion over long timescales can lead to the formation of a mini black hole. In this scenario, the neutron star is destabilized and implodes to form a black hole without significantly increasing its mass. When this process occurs in neutron stars in coalescing binaries, one or both stars might be converted to a black hole before they merge. Thus, in the mass range of $\sim \mbox{1--2}\, M_\odot,$ the Universe might contain three distinct populations of compact binaries: one containing only neutron stars, the second population of only black holes, and a third, mixed population consisting of a neutron star and a black hole. However, it is unlikely to have a mixed population as the various timescales allow for both neutron stars to remain or collapse within a short timescale. In this paper, we explore the capability of future gravitational-wave detector networks, including upgrades of Advanced LIGO and Virgo, and new facilities such as the Cosmic Explorer and Einstein Telescope (XG network), to discriminate between different populations by measuring the effective tidal deformability of the binary, which is zero for binary black holes but nonzero for binary neutron stars. Furthermore, we show that observing the relative abundances of the different populations can be used to infer the timescale for neutron stars to implode into black holes, and in turn, provide constraints on the particle nature of dark matter. The XG network will infer the implosion timescale to within an accuracy of 0.01 Gyr at 90% credible interval and determine the dark matter mass and interaction cross section to within a factor of 2 GeV and 10 cm$^{-2}$, respectively.

The goal of this paper is to introduce a novel likelihood-based inferential framework for axion haloscopes which is valid under the commonly applied "rescanning" protocol. The proposed method enjoys short data acquisition times and a simple tuning of the detector configuration. Local statistical significance and power are computed analytically, avoiding the need of burdensome simulations. Adequate corrections for the look-elsewhere effect are also discussed. The performance of our inferential strategy is compared with that of a simple method which exploits the geometric probability of rescan. Finally, we exemplify the method with an application to a HAYSTAC type axion haloscope.

The detection of gravitational waves from coalescences of binary compact stars by current interferometry experiments has opened up a new era of gravitational-wave astrophysics and cosmology. The search for stochastic gravitational-wave background is underway by correlating signals from a pair of detectors in the detector network formed by the LIGO, Virgo, and KAGRA. In a previous work, we have developed a method based on spherical harmonic expansion to calculate the overlap reduction functions of the LIGO-Virgo-KAGRA network for a polarized stochastic gravitational-wave background. In this work, we will apply the method to calculate the overlap reduction functions of third-generation detectors such as a ground-based network linking the Einstein Telescope and the Cosmic Explorer, and the LISA-Taiji joint space mission.

Space-time parity can solve the strong CP problem and introduces a spontaneously broken $SU(2)_R$ gauge symmetry. We investigate the possibility of baryogenesis from a first-order $SU(2)_R$ phase transition similar to electroweak baryogenesis. We consider a model with the minimal Higgs content, for which the strong CP problem is indeed solved without introducing extra symmetry beyond parity. Although the parity symmetry seems to forbid the $SU(2)_R$ anomaly of the $B-L$ symmetry, the structure of the fermion masses can allow for the $SU(2)_R$ sphaleron process to produce non-zero $B-L$ asymmetry of Standard Model particles so that the wash out by the $SU(2)_L$ sphaleron process is avoided. The setup predicts a new hyper-charged fermion whose mass is correlated with the $SU(2)_R$ symmetry breaking scale and hence with the $SU(2)_R$ gauge boson mass, and depending on the origin of CP violation, with an electron electric dipole moment. In a setup where CP violation and the first-order phase transition are assisted by a singlet scalar field, the singlet can be searched for at future colliders.

We generalize Einstein's General Relativity (GR) by assuming that all matter (including macro-objects) has quantum effects. An appropriate theory to fulfill this task is Gauge Theory Gravity (GTG) developed by the Cambridge group. GTG is a ``spin-torsion" theory, according to which, gravitational effects are described by a pair of gauge fields defined over a flat Minkowski background spacetime. The matter content is completely described by the Dirac spinor field, and the quantum effects of matter are identified as the spin tensor derived from the spinor field. The existence of the spin of matter results in the torsion field defined over spacetime. Torsion field plays the role of Bohmian quantum potential which turns out to be a kind of repulsive force as opposed to the gravitational potential which is attractive. The equivalence principle remains and essential in this theory so that GR is relegated to a locally approximate theory wherein the quantum effects (torsion) are negligible. As a toy model, we assume that the macro matter content can be described by the covariant Dirac equation and apply this theory to the simplest radially symmetric and static gravitational systems. Consequently, by virtue of the cosmological principle, we are led to a static universe model in which the Hubble redshifts arise from the torsion fields.

Gravitational wave memory effects arise from non-oscillatory components of gravitational wave signals, and they are predictions of general relativity in the nonlinear regime that have close connections to the asymptotic properties of isolated gravitating systems. There are many types of memory effects that have been studied in the literature. In this paper we focus on the "displacement" and "spin" memories, which are expected to be the largest of these effects from sources such as the binary black hole mergers which have already been detected by LIGO and Virgo. The displacement memory is a change in the relative separation of two initially comoving observers due to a burst of gravitational waves, whereas the spin memory is a portion of the change in relative separation of observers with initial relative velocity. As both of these effects are small, LIGO, Virgo, and KAGRA can only detect memory effects from individual events that are much louder (and thus rarer) than those that have been detected so far. By combining data from multiple events, however, these effects could be detected in a population of binary mergers. In this paper, we present new forecasts for how long current and future detectors will need to operate in order to measure these effects from populations of binary black hole systems that are consistent with the populations inferred from the detections from LIGO and Virgo's first three observing runs. We find that a second-generation detector network of LIGO, Virgo, and KAGRA operating at the O4 ("design") sensitivity for 1.5 years and then operating at the O5 ("plus") sensitivity for an additional year can detect the displacement memory. For Cosmic Explorer, we find that displacement memory could be detected for individual loud events, and that the spin memory could be detected in a population after 5 years of observation time.

We derive soft theorems for theories in which time symmetries -- symmetries that involve the transformation of time, an example of which are Lorentz boosts -- are spontaneously broken. The soft theorems involve unequal-time correlation functions with the insertion of a soft Goldstone in the far past. Explicit checks are provided for several examples, including the effective theory of a relativistic superfluid and the effective field theory of inflation. We discuss how in certain cases these unequal-time identities capture information at the level of observables that cannot be seen purely in terms of equal-time correlators of the field alone. We also discuss when it is possible to phrase these soft theorems as identities involving equal-time correlators.

We present a parameter estimation framework for gravitational wave (GW) signals that brings together several ideas to accelerate the inference process. First, we use the relative binning algorithm to evaluate the signal-to-noise-ratio timeseries in each detector for a given choice of intrinsic parameters. Second, we decouple the estimation of the intrinsic parameters (such as masses and spins of the components) from that of the extrinsic parameters (such as distance, orientation, and sky location) that describe a binary compact object coalescence. We achieve this by semi-analytically marginalizing the posterior distribution over extrinsic parameters without repeatedly evaluating the waveform for a fixed set of intrinsic parameters. Finally, we augment samples of intrinsic parameters with extrinsic parameters drawn from their appropriate conditional distributions. We implement the method for binaries with aligned spins, restricted to the quadrupole mode of the signal. Using simulated GW signals, we demonstrate that the method produces full eleven-dimensional posteriors that match those from standard Bayesian inference. Our framework takes only ~200 seconds to analyze a typical binary-black-hole signal and ~250 seconds to analyze a typical binary-neutron-star signal using one computing core. Such real-time and accurate estimation of the binary source properties will greatly aid the interpretation of triggers from gravitational wave searches, as well as the search for possible electromagnetic counterparts. We make the framework publicly available via the GW inference package cogwheel.

We study the possible gravitational wave signal and the viability of baryogenesis arising from the electroweak phase transition in an extension of the Standard Model (SM) by a scalar singlet field without a $\mathbb{Z}_2$ symmetry. We first analyze the velocity of the expanding true-vacuum bubbles during the phase transition, confirming our previous finding in the unbroken $\mathbb{Z}_2$ symmetry scenario, where the bubble wall velocity can be computed from first principles only for weak transitions with strength parameters $\alpha \lesssim 0.05$, and the Chapman-Jouguet velocity defines the maximum velocity for which the wall is stopped by the friction from the plasma. We further provide an analytical approximation to the wall velocity in the general scalar singlet scenario without $\mathbb{Z}_2$ symmetry and test it against the results of a detailed calculation, finding good agreement. We show that in the singlet scenario with a spontaneously broken $\mathbb{Z}_2$ symmetry, the phase transition is always weak and we see no hope for baryogenesis. In contrast, in the case with explicit $\mathbb{Z}_2$ breaking there is a region of the parameter space producing a promising baryon yield in the presence of CP violating interactions via an effective operator involving the singlet scalar and the SM top quarks. Yet, we find that this region yields unobservable gravitational waves. Finally, we show that the promising region for baryogenesis in this model may be fully tested by direct searches for singlet-like scalars in di-boson final states at the HL-LHC, combined with present and future measurements of the electron electric dipole moment.