Papers for Thursday, Nov 30 2023

Two unresolved questions at galaxy centers, namely the formation of the nuclear star cluster (NSC) and the origin of the gamma-ray excess in the Milky Way (MW) and Andromeda (M31), are both related to the formation and evolution of globular clusters (GCs). They migrate towards the galaxy center due to dynamical friction, and get tidally disrupted to release the stellar mass content including millisecond pulsars (MSPs), which contribute to the NSC and gamma-ray excess. In this study, we propose a semi-analytical model of GC formation and evolution that utilizes the Illustris cosmological simulation to accurately capture the formation epochs of GCs and simulate their subsequent evolution. Our analysis confirms that our GC properties at z=0 are consistent with observations, and our model naturally explains the formation of a massive NSC in a galaxy similar to the MW and M31. We also find a remarkable similarity in our model prediction with the gamma-ray excess signal in the MW. However, our predictions fall short by approximately an order of magnitude in M31, indicating distinct origins for the two gamma-ray excesses. Meanwhile, we utilize the catalog of Illustris halos to investigate the influence of galaxy assembly history. We find that the earlier a galaxy is assembled, the heavier and spatially more concentrated its GC system behaves at z=0. This results in a larger NSC mass and brighter gamma-ray emission from deposited MSPs

In the past decade the study of exoplanet atmospheres at high-spectral resolution, via transmission/emission spectroscopy and cross-correlation techniques for atomic/molecular mapping, has become a powerful and consolidated methodology. The current limitation is the signal-to-noise ratio during a planetary transit. This limitation will be overcome by ANDES, an optical and near-infrared high-resolution spectrograph for the ELT. ANDES will be a powerful transformational instrument for exoplanet science. It will enable the study of giant planet atmospheres, allowing not only an exquisite determination of atmospheric composition, but also the study of isotopic compositions, dynamics and weather patterns, mapping the planetary atmospheres and probing atmospheric formation and evolution models. The unprecedented angular resolution of ANDES, will also allow us to explore the initial conditions in which planets form in proto-planetary disks. The main science case of ANDES, however, is the study of small, rocky exoplanet atmospheres, including the potential for biomarker detections, and the ability to reach this science case is driving its instrumental design. Here we discuss our simulations and the observing strategies to achieve this specific science goal. Since ANDES will be operational at the same time as NASA's JWST and ESA's ARIEL missions, it will provide enormous synergies in the characterization of planetary atmospheres at high and low spectral resolution. Moreover, ANDES will be able to probe for the first time the atmospheres of several giant and small planets in reflected light. In particular, we show how ANDES will be able to unlock the reflected light atmospheric signal of a golden sample of nearby non-transiting habitable zone earth-sized planets within a few tenths of nights, a scientific objective that no other currently approved astronomical facility will be able to reach.

Strong gravitational lensing is a powerful probe of the distribution of matter on sub-kpc scales. It can be used to test the existence of completely dark subhalos surrounding galaxies, as predicted by the standard cold dark matter model, or to test alternative dark matter models. The constraining power of the method depends strongly on photometric and astrometric precision and accuracy. We simulate and quantify the capabilities of upcoming adaptive optics systems and advanced instruments on ground-based telescopes, focusing as an illustration on the Keck Telescope (OSIRIS + KAPA, LIGER + KAPA) and the Thirty Meter Telescope (TMT; IRIS + NFIRAOS). We show that these new systems will achieve dramatic improvements over current ones in both photometric and astrometric precision. Narrow line flux ratio errors below $2\%$, and submilliarcsecond astrometric precision will be attainable for typical quadruply imaged quasars. With TMT, the exposure times required will be of order a few minutes per system, enabling the follow-up of 100-1000 systems expected to be discovered by the Rubin, Euclid, and Roman Telescopes.

We introduce a diffusion-based generative model to describe the distribution of galaxies in our Universe directly as a collection of points in 3-D space (coordinates) optionally with associated attributes (e.g., velocities and masses), without resorting to binning or voxelization. The custom diffusion model can be used both for emulation, reproducing essential summary statistics of the galaxy distribution, as well as inference, by computing the conditional likelihood of a galaxy field. We demonstrate a first application to massive dark matter haloes in the Quijote simulation suite. This approach can be extended to enable a comprehensive analysis of cosmological data, circumventing limitations inherent to summary statistic -- as well as neural simulation-based inference methods.

Astronomical transients, such as supernovae and other rare stellar explosions, have been instrumental in some of the most significant discoveries in astronomy. New astronomical sky surveys will soon record unprecedented numbers of transients as sparsely and irregularly sampled multivariate time series. To improve our understanding of the physical mechanisms of transients and their progenitor systems, early-time measurements are necessary. Prioritizing the follow-up of transients based on their age along with their class is crucial for new surveys. To meet this demand, we present the first method of predicting the age of transients in real-time from multi-wavelength time-series observations. We build a Bayesian probabilistic recurrent neural network. Our method can accurately predict the age of a transient with robust uncertainties as soon as it is initially triggered by a survey telescope. This work will be essential for the advancement of our understanding of the numerous young transients being detected by ongoing and upcoming astronomical surveys.

Parsec-scale separation supermassive black hole binaries in the centre of gas-rich galaxy merger remnants could be surrounded by massive circumbinary discs (CBDs). Black hole mass and spin evolution during the gas-rich binary inspiral are crucial in determining the direction and power of relativistic jets that radio observations with LOFAR and SKAO will probe, and for predicting gravitational wave (GW) emission that IPTA and LISA will measure. We present 3D hydrodynamic simulations capturing gas-rich, self-gravitating CBDs around a $2\times 10^6$M$_{\odot}$ supermassive black hole binary, that probe different mass ratios, eccentricities and inclinations. We employ a sub-grid Shakura-Sunyaev accretion disc to self-consistently model black hole mass and spin evolution together with super-Lagrangian refinement techniques to resolve gas flows, streams and mini-discs within the cavity, which play a fundamental role in torquing and feeding the binary. We find that higher mass ratio and eccentric binaries result in larger cavities, while retrograde binaries result in smaller cavities. All of the simulated binaries are expected to shrink with net gravitational torques being negative. Unlike previous simulations, we do not find preferential accretion onto the secondary black hole. This implies smaller chirp masses at coalescence and hence a weaker GW background. Critically this means that spin-alignment is faster than the binary inspiral timescale even for low mass ratios. However, we find that mini-disc and hence spin alignment is not guaranteed in initially misaligned systems, potentially leading to a significant fraction of recoiled remnants displaced from their host galaxies if chaotic accretion is the dominant feeding channel.

The third data release of Gaia introduced a large catalog of astrometric binaries, out of which about 3,200 are likely main-sequence stars with a white-dwarf (WD) companion. These binaries are typically found with orbital separations of ~1AU, a separation range that was largely unexplored due to observational challenges. Such systems are likely to have undergone a phase of stable mass transfer while the WD progenitor was on the asymptotic giant branch. Here we study the WD mass distribution of a volume-complete sample of binaries with K/M-dwarf primaries and orbital separations ~1AU. We find that the number of massive WDs relative to the total number of WDs in these systems is smaller by an order of magnitude compared to their occurrence among single WDs in the field. One possible reason can be an implicit selection of the WD mass range if these are indeed post-stable-mass-transfer systems. Another reason can be the lack of merger products in our sample compared to the field, due to the relatively tight orbital separations of these systems.

We employ a data-driven approach to investigate the rigidity and spatial dependence of the diffusion of cosmic rays in the turbulent magnetic field of the Milky Way. Our analysis combines data sets from the experiments Voyager, AMS-02, CALET, and DAMPE for a range of cosmic ray nuclei from protons to oxygen. Our findings favor models with a smooth behavior in the diffusion coefficient, indicating a good qualitative agreement with the predictions of self-generated magnetic turbulence models. Instead, the current cosmic-ray data do not exhibit a clear preference for or against inhomogeneous diffusion, which is also a prediction of these models. Future progress might be possible by combining cosmic-ray data with gamma rays or radio observations, enabling a more comprehensive exploration.

V510 Pup (IRAS 08005-2356) is a binary post-AGB system with a fast molecular outflow that has been noted for its puzzling mixture of carbon- and oxygen-rich features in the optical and infrared. To explore this chemical dichotomy and relate it to the kinematics of the source, we present an ACA spectral line survey detailing fourteen newly detected molecules in this pre-planetary nebula. The simultaneous presence of CN/C2H/HC3N and SO/SO2 support the previous conclusion of mixed chemistry, and their line profiles indicate that the C- and O-rich material trace distinct velocity structures in the outflow. This evidence suggests that V510 Pup could harbor a dense O-rich central waist from an earlier stage of evolution, which persisted after a fast C-rich molecular outflow formed. By studying the gas phase composition of this unique source, we aim to reveal new insights into the interplay between dynamics and chemistry in rapidly evolving post-AGB systems.

The recurrent nova RS Ophiuchi (RS Oph) underwent its most recent eruption on 8 August 2021 and became the first nova to produce both detectable GeV and TeV emission. We used extensive X-ray monitoring with the Neutron Star Interior Composition Explorer Mission (NICER) to model the X-ray spectrum and probe the shock conditions throughout the 2021 eruption. The rapidly evolving NICER spectra consisted of both line and continuum emission that could not be accounted for using a single-temperature collisional equilibrium plasma model with an absorber that fully covered the source. We successfully modelled the NICER spectrum as a non-equilibrium ionization collisional plasma with partial-covering absorption. The temperature of the the non-equilibrium plasma show a peak on Day 5 with a kT of approximately 24 keV. The increase in temperature during the first five days could have been due to increasing contribution to the X-ray emission from material behind fast polar shocks or a decrease is the amount of energy being drained from shocks into particle acceleration during that time period. The absorption showed a change from fully covering the source to having a covering fraction of roughly 0.4, suggesting a geometrical evolution of the shock region within the complex global distribution of the circumstellar material. These findings show the evidence of the ejecta interacting with some dense equatorial shell initially and with less dense material in the bipolar regions at later times during the eruption.

The recent detection of the live isotopes $^{60}{\rm Fe}$ and $^{244}{\rm Pu}$ in deep ocean sediments dating back to the past 3-4 Myr poses a serious challenge to the identification of their production site(s). While $^{60}{\rm Fe}$ is usually attributed to standard supernovae, actinides are r-process nucleosynthesis yields, which are believed to be synthesized in rare events, such as special classes of supernovae or binary mergers involving at least one neutron star. Previous works concluded that a single binary neutron star merger cannot explain the observed isotopic ratio. In this work, we consider a set of numerical simulations of binary neutron star mergers producing long-lived massive remnants expelling both dynamical and spiral-wave wind ejecta. The latter, due to a stronger neutrino irradiation, produce also iron-group elements. Assuming that large scale mixing is inefficient before the fading of the kilonova remnant and that the spiral-wave wind is sustained over a 100-200 ms timescale, the ejecta emitted at mid-high latitudes provide a $^{244}{\rm Pu}$ over $^{60}{\rm Fe}$ ratio compatible with observations. The merger could have happened 110-200 pc away from the Earth, and between 3.5 and 4.5 Myr ago. We also compute expected isotopic ratios for 8 other live radioactive nuclides showing that the proposed binary neutron star merger scenario is distinguishable from other scenarios proposed in the literature.

Isocurvature fluctuations, where the relative number density of particle species spatially varies, can be generated from initially adiabatic or curvature fluctuations if the various species fall out of or were never in thermal equilibrium. The freezing of the thermal relic dark matter abundance is one such case, but for modes that are still outside the horizon the amplitude is highly suppressed and originates from the small change in the local expansion rate due to the local space curvature produced by the curvature fluctuation. We establish a simple separate-universe method for calculating this generation that applies to both freeze-in and freeze-out models, identify three critical epochs for this process, and give general scaling behaviors for the amplitude in each case: the freezing epoch, the kinetic decoupling epoch and matter-radiation equality. Freeze-out models are typically dominated by spatially modulated annihilation from the latter epochs and can generate much larger isocurvature fluctuations compared with typical freeze-in models, albeit still very small and observationally allowed by cosmic microwave background measurements. We illustrate these results with concrete models where the dark matter interactions are vector or scalar mediated.

We make a comparison of deep SCUBA-2 450 $\mu$m and 850 $\mu$m imaging on the massive lensing cluster field Abell 2744 with Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm data. Our primary goal is to assess how effective the wider-field SCUBA-2 sample, in combination with red JWST priors, is for finding faint dusty star-forming galaxies (DSFGs) compared to the much more expensive mosaicked ALMA observations. We cross-match our previously reported direct ($>5\sigma$) SCUBA-2 sample and red JWST NIRCam prior-selected ($>3\sigma$) SCUBA-2 sample to direct ALMA sources from the DUALZ survey. We find that roughly 95% are confirmed by ALMA. The red priors also allow us to probe deeper in the ALMA image. Next, by measuring the 450 $\mu$m and 850 $\mu$m properties of the full ALMA sample, we show that 46/69 of the ALMA sources are detected at 850 $\mu$m and 24/69 are detected at 450 $\mu$m in the SCUBA-2 images, with a total detection fraction of nearly 75%. All of the robust ($>5\sigma$) ALMA sources that are not detected in at least one SCUBA-2 band lie at 1.2 mm fluxes $\lesssim$ 0.6 mJy and are undetected primarily due to the higher SCUBA-2 flux limits. We also find that the SCUBA-2 detection fraction drops slightly beyond $z=3$, which we attribute to the increasing 1.2 mm to 850 $\mu$m and 1.2 mm to 450 $\mu$m flux ratios combined with the ALMA selection. The results emphasize the power of combining SCUBA-2 data with JWST colors to map the faint DSFG population.

We report the ALMA Band 7 observations of 86 Herschel sources that likely contain gravitationally-lensed galaxies. These sources are selected with relatively faint 500 $\mu$m flux densities between 15 to 85 mJy in an effort to characterize the effect of lensing across the entire million-source Herschel catalogue. These lensed candidates were identified by their close proximity to bright galaxies in the near-infrared VISTA Kilo-Degree Infrared Galaxy Survey (VIKING) survey. Our high-resolution observations (0.15 arcsec) confirm 47 per cent of the initial candidates as gravitational lenses, while lensing cannot be excluded across the remaining sample. We find average lensing masses (log M/M$_{\odot}$ = 12.9 $\pm$ 0.5) in line with previous experiments, although direct observations might struggle to identify the most massive foreground lenses across the remaining 53 per cent of the sample, particularly for lenses with larger Einstein radii. Our observations confirm previous indications that more lenses exist at low flux densities than expected from strong galaxy-galaxy lensing models alone, where the excess is likely due to additional contributions of cluster lenses and weak lensing. If we apply our method across the total 660 sqr. deg. H-ATLAS field, it would allow us to robustly identify 3000 gravitational lenses across the 660 square degree Herschel ATLAS fields.

Cosmological distances are fundamental observables in cosmology. The luminosity ($D_L$), angular diameter ($D_A$) and gravitational wave ($D_{\rm GW}$) distances are all trivially related in General Relativity assuming no significant absorption of photons in the extragalactic medium, also known as cosmic opacity. Supernovae have long been the main cosmological standard candle for the past decades, but bright standard sirens are now a proven alternative, with the advantage of not requiring calibration with other astrophysical sources. Moreover, they can also measure deviations from modified gravity since they can provide evidence for a discrepancy between $D_L$ and $D_{\rm GW}$. However, both gravitational and cosmological parameters are degenerate in the Hubble diagram, making it hard to properly detect beyond standard model physics. Finally, recently a model-independent method was proposed to infer angular diameter distances from large-scale structure which is independent of both early universe and dark energy physics. In this paper we propose a tripartite test of the ratios of these three distances with minimal amount of assumptions regarding cosmology, the early universe, cosmic opacity and modified gravity. We proceed to forecast this test with a combination of uncalibrated LSST and Roman supernovae, Einstein Telescope bright sirens and a joint DESI-like + Euclid-like galaxy survey. We find that even in this very model-independent approach we will be able to detect, in each of many redshift bins, percent-level deviations in these ratios of distances, allowing for very precise consistency checks of $\Lambda$CDM and standard physics.

Observations from wide-field quasar surveys indicate that the quasar auto-correlation length increases dramatically from $z\approx2.5$ to $z\approx4$. This large clustering amplitude at $z\approx4$ has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter halos residing in the most extreme environments at that redshift. In this work, we present a model that simultaneously reproduces both the observed quasar auto-correlation and quasar luminosity functions. The spatial distribution of halos and their relative abundance are obtained via a novel method that computes the halo mass and halo cross-correlation functions by combining multiple large-volume dark-matter-only cosmological simulations with different box sizes and resolutions. Armed with these halo properties, our model exploits the conditional luminosity function framework to describe the stochastic relationship between quasar luminosity, $L$, and halo mass, $M$. Assuming a simple power-law relation $L\propto M^\gamma$ with log-normal scatter, $\sigma$, we are able to reproduce observations at $z\sim 4$ and find that: (a) the quasar luminosity-halo mass relation is highly non-linear ($\gamma\gtrsim2$), with very little scatter ($\sigma \lesssim 0.3$ dex); (b) luminous quasars ($\log_{10} L/{\rm erg s}^{-1} \gtrsim 46.5-47$) are hosted by halos with mass $\log_{10} M/{\rm M}_\odot\gtrsim 13-13.5$; and (c) the implied duty cycle for quasar activity approaches unity ($\varepsilon_{\rm DC}\approx10-60\%$). We also consider observations at $z\approx2.5$ and find that the quasar luminosity-halo mass relation evolves significantly with cosmic time, implying a rapid change in quasar host halo masses and duty cycles, which in turn suggests concurrent evolution in black hole scaling relations and/or accretion efficiency.

The zero-age main sequence (ZAMS) is a critical phase for stellar angular momentum evolution, as stars transition from contraction-dominated spin-up to magnetic wind-dominated spin-down. We present the first robust observational constraints on rotation for FGK stars at $\approx40$ Myr. We have analyzed TESS light curves for 1410 members of five young open clusters with ages between 25-55 Myr: IC 2391, IC 2602, NGC 2451A, NGC 2547, and Collinder 135. In total, we measure 868 rotation periods, including 96 new, high-quality periods for stars around 1 ${M_{\odot}}$. This is an increase of ten times the existing literature sample at the ZAMS. We then use the $\tau^2$ method to compare our data to models for stellar angular momentum evolution. Although the ages derived from these rotation models do not match isochronal ages, we show these observations can clearly discriminate between different models for stellar wind torques. Finally, $\tau^2$ fits indicate that magnetic braking and/or internal angular momentum transport significantly impact rotational evolution even on the pre-main sequence.

FU Orionis-type objects(FUors) are embedded protostars that undergo episodes of high accretion, potentially indicating a widespread but poorly understood phase in the formation of low-mass stars. Gaining a better understanding of the influence exerted by these outbursts on the evolution of the surrounding protoplanetary disc may hold significant implications for the process of planet formation and the evolution of disc chemistry. The heating due to outbursts of high accretion in FUors pushes the snowlines of key volatiles farther out in the disc, so they become easier to observe and study. Among the known FUors, V883 Ori is of particular interest. V883 Ori was the first FUor to show indirect evidence of a resolvable snowline beyond 40 au. By introducing a radial-dependent model of this source including viscous heating, we show that active heating is needed to reproduce the steep thermal profile of dust in the inner disc of V883 Ori. Our disc modeling combines the effect of stellar irradiation and the influence on the disc shape caused by the outburst of accretion. The accuracy of our model is tested by comparing synthetic ALMA images with continuum observations of V883 Ori, showing that the model successfully reproduces the 1.3 mm emission of V883 Ori at high spatial resolution. Our final predictions underline the importance of viscous heating as a predominant heat source for this type of object, changing the physical conditions (shape and temperature) of the disc, and influencing its evolution.

We present an atlas and follow-up spectroscopic observations of 87 thin stream-like structures detected with the STREAMFINDER algorithm in Gaia DR3, of which 29 are new discoveries. Here we focus on using these streams to refine mass models of the Galaxy. Fits with a double power law halo with the outer power law slope set to $-\beta_h=3$ yield an inner power law slope $-\gamma_h=0.97^{+0.17}_{-0.21}$, a scale radius of $r_{0, h}=14.7^{+4.7}_{-1.0}$ kpc, a halo density flattening $q_{m, h}=0.75\pm0.03$, and a local dark matter density of $\rho_{h, \odot}=0.0114\pm0.0007 {\rm M_\odot pc^{-3}}$. Freeing $\beta$ yields $\beta=2.53^{+0.42}_{-0.16}$, but this value is heavily influenced by our chosen virial mass limit. The stellar disks are found to have a combined mass of $4.20^{+0.44}_{-0.53}\times10^{10} {\rm M_\odot}$, with the thick disk contributing $12.4\pm0.7$\% to the local stellar surface density. The scale length of the thin and thick disks are $2.17^{+0.18}_{-0.08}$ kpc and $1.62^{+0.72}_{-0.13}$ kpc, respectively, while their scale heights are $0.347^{+0.007}_{-0.010}$ kpc and $0.86^{+0.03}_{-0.02}$ kpc, respectively. The virial mass of the favored model is $M_{200}=1.09^{+0.19}_{-0.14}\times 10^{12} {\rm M_\odot}$, while the mass inside of 50 kpc is $M_{R<50}=0.46\pm0.03\times 10^{12} {\rm M_\odot}$. We introduce the Large Magellanic Cloud (LMC) into the derived potential models, and fit the "Orphan" stream therein, finding a mass for the LMC that is consistent with recent estimates. Some highlights of the atlas include the nearby trailing arm of $\omega$-Cen, and a nearby very metal-poor stream that was once a satellite of the Sagittarius dwarf galaxy. Finally, we unambiguously detect a hot component around the GD-1 stream, consistent with it having been tidally pre-processed within its own DM subhalo.

Central Compact Objects (CCOs), neutron stars found near the centre of some Supernova Remnants (SNRs), have been almost exclusively studied in X-rays and are thought to lack the wind nebulae typically seen around young, rotation-powered pulsars. We present the first, spatially-resolved, morphological and spectroscopic study of the optical nebula observed at the location of CXOU J085201.4-461753, the CCO in the heart of the Vela Junior SNR. It is currently the only Galactic CCO with a spatially coincident nebula detected at optical wavelengths, whose exact nature remains uncertain. New MUSE integral field spectroscopy data confirm that the nebula, shaped like a smooth blob extending 8" in diameter, is dominated by [N II]$\lambda\lambda$6548,6583 emission. The data reveals a distinct and previously unobserved morphology of the H$\alpha$ emission, exhibiting an arc-like shape reminiscent of a bow shock nebula. We observe a significantly strong [N II] emission relative to H$\alpha$, with the [N II]$\lambda\lambda$6548,6583 up to 34 times the intensity of the H$\alpha$ emission within the optical nebula environment. Notably, the [N II] and H$\alpha$ structures are not spatially coincident, with the [N II] nebula concentrated to the south of the CCO and delimited by the H$\alpha$ arc-like structure. We detect additional emission in [N I], He I, [S II], [Ar III], [Fe II], and [S III]. We discuss our findings in the light of a photoionization or Wolf-Rayet nebula, pointing to a very massive progenitor and further suggesting that very massive stars do not necessarily make black holes.

The canonical picture of star formation involves disk-mediated accretion, with Keplerian accretion disks and associated bipolar jets primarily observed in nearby, low-mass young stellar objects (YSOs). Recently, rotating gaseous structures and Keplerian disks have been detected around a number of massive (M > 8 solar masses) YSOs (MYSOs) including several disk-jet systems. All of the known MYSO systems are located in the Milky Way, and all are embedded in their natal material. Here we report the detection of a rotating gaseous structure around an extragalactic MYSO in the Large Magellanic Cloud. The gas motions show radial flow of material falling from larger scales onto a central disk-like structure, the latter exhibiting signs of Keplerian rotation, i.e., a rotating toroid feeding an accretion disk and thus the growth of the central star. The system is in almost all aspects comparable to Milky Way high-mass young stellar objects accreting gas via a Keplerian disk. The key difference between this source and its Galactic counterparts is that it is optically revealed, rather than being deeply embedded in its natal material as is expected of such a young massive star. We suggest that this is the consequence of the star having formed in a low-metallicity and low-dust content environment, thus providing important constraints for models of the formation and evolution of massive stars and their circumstellar disks.

Since 2008, the Interstellar Boundary Explorer (IBEX) satellite has been gathering data on heliospheric energetic neutral atoms (ENAs) while being exposed to various sources of background noise, such as cosmic rays and solar energetic particles. The IBEX mission initially released only a qualified triple-coincidence (qABC) data product, which was designed to provide observations of ENAs free of background contamination. Further measurements revealed that the qABC data was in fact susceptible to contamination, having relatively low ENA counts and high background rates. Recently, the mission team considered releasing a certain qualified double-coincidence (qBC) data product, which has roughly twice the detection rate of the qABC data product. This paper presents a simulation-based validation of the new qBC data product against the already-released qABC data product. The results show that the qBCs can plausibly be said to share the same signal rate as the qABCs up to an average absolute deviation of 3.6%. Visual diagnostics at an orbit, map, and full mission level provide additional confirmation of signal rate coherence across data products. These approaches are generalizable to other scenarios in which one wishes to test whether multiple observations could plausibly be generated by some underlying shared signal.

Three-minute oscillations are a common phenomenon in the solar chromosphere above a sunspot. Oscillations can be affected by the energy release process related to solar flares. In this paper, we report on an enhanced oscillation in flare event SOL2012-07-05T21:42 with a period of around three minutes, that occurred at the location of a flare ribbon at a sunspot umbra-penumbra boundary, and was observed both in chromo-spheric and coronal passbands. An analysis of this oscillation was carried out using simultaneous ground-based observations from the Goode Solar Telescope (GST) at the Big Bear Solar Observatory (BBSO) and space-based observations from the Solar Dynamics Observatory (SDO). A frequency shift was observed before and after the flare, with the running penumbral wave that was present with a period of about 200 s before the flare co-existing with a strengthened oscillation with a period of 180 s at the same locations after the flare. We also found a phase difference between different passbands, with the oscillation occurring from high-temperature to low-temperature passbands. Theoretically, the change in frequency is strongly dependent on the variation of the inclination of the magnetic field and the chromospheric temperature. Following an analysis of the properties of the region, we find the frequency change is caused by the slight decrease of the magnetic inclination angle to the local vertical. In addition, we suggest that the enhanced three-minute oscillation is related to the additional heating, maybe due to the downflow, during the EUV late phase of the flare.

One of the many exciting revelations of the New Horizons flyby of Pluto was the observation of global haze layers at altitudes as high as 200 km in the visible wavelengths. This haze is produced in the upper atmosphere through photochemical processes, similar to the processes in Titan's atmosphere. As the haze particles grow in size and descend to the lower atmosphere, they coagulate and interact with the gases in the atmosphere through condensation and sticking processes that serve as temporary and permanent loss processes. New Horizons observations confirm studies of Titan haze analogs suggesting that photochemically produced haze particles harden as they grow in size. We outline in this chapter what is known about the photochemical processes that lead to haze production and outline feedback processes resulting from the presence of haze in the atmosphere, connect this to the evolution of Pluto's atmosphere, and discuss open questions that need to be addressed in future work.

We study preheating via kinetic couplings after dilaton-axion $\alpha$-attractor inflation. We focus on E-model $\alpha$-attractor driven inflation where the inflaton is kinetically coupled to an ultralight axion. In this class of models, the kinetic coupling is related to the form of the potential, and once the amplitude of the scalar curvature spectrum as well as the tensor-to-scalar ratio are specified, the model has no free parameters. We find that kinetic preheating can be extremely efficient, with stronger preheating occurring at parameter values corresponding to smaller values of the tensor-to-scalar ratio. Preheating becomes extremely efficient below $r \lesssim 1.6\times 10^{-5}$.

Upcoming surveys are predicted to discover galaxy-scale strong lenses on the order of $10^5$, making deep learning methods necessary in lensing data analysis. Currently, there is insufficient real lensing data to train deep learning algorithms, but the alternative of training only on simulated data results in poor performance on real data. Domain Adaptation may be able to bridge the gap between simulated and real datasets. We utilize domain adaptation for the estimation of Einstein radius ($\Theta_E$) in simulated galaxy-scale gravitational lensing images with different levels of observational realism. We evaluate two domain adaptation techniques - Domain Adversarial Neural Networks (DANN) and Maximum Mean Discrepancy (MMD). We train on a source domain of simulated lenses and apply it to a target domain of lenses simulated to emulate noise conditions in the Dark Energy Survey (DES). We show that both domain adaptation techniques can significantly improve the model performance on the more complex target domain dataset. This work is the first application of domain adaptation for a regression task in strong lensing imaging analysis. Our results show the potential of using domain adaptation to perform analysis of future survey data with a deep neural network trained on simulated data.

The Distributed Electronic Cosmic-ray Observatory (DECO) is a cell phone app that uses a cell phone camera image sensor to detect cosmic-ray particles and particles from radioactive decay. Images recorded by DECO are classified by a convolutional neural network (CNN) according to their morphology. In this project, we develop a GEANT4-derived simulation of particle interactions inside the CMOS sensor using the Allpix$^2$ modular framework. We simulate muons, electrons, and photons with energy range 10 keV to 100 GeV, and their deposited energy agrees well with expectations. Simulated events are recorded and processed in a similar way as data images taken by DECO, and the result shows both similar image morphology with data events and good quantitative data-Monte Carlo agreement.

In the previous work of Xu & Peng (2021), we investigated the structural and environmental dependence on quenching in the nearby universe. In this work we extend our investigations to higher redshifts by combining galaxies from SDSS and ZFOURGE surveys. In low density, we find a characteristic $\Sigma_{1\ kpc}$ above which the quenching is initiated as indicated by their population-averaged color. $\Sigma^{crit}_{1\ kpc}$ shows only weakly mass-dependency at all redshifts, which suggests that the internal quenching process is more related to the physics that acts in the central region of galaxies. In high density, $\Sigma^{crit}_{1\ kpc}$ for galaxies at $z > 1$ is almost indistinguishable with their low-density counterparts. At $z < 1$, $\Sigma^{crit}_{1\ kpc}$ for low-mass galaxies becomes progressively strongly mass-dependent, which is due to the increasingly stronger environmental effects at lower redshifts. $\Sigma^{crit}_{1\ kpc}$ in low density shows strong redshift evolution with $\sim 1$ dex decrement from $z = 2.5$ to $z = 0$. It is likely due to that at a given stellar mass, the host halo is on average more massive and gas-rich at higher redshifts, hence a higher level of integrated energy from more massive black hole is required to quench. As the halo evolves from cold to hot accretion phase at lower redshifts, the gas is shock-heated and becomes more vulnerable to AGN feedback processes, as predicted by theory. Meanwhile, angular momentum quenching also becomes more effective at low redshifts, which complements a lower level of integrated energy from black hole to quench.

It has been suggested that small chemical anomalies observed in planet-hosting wide binary systems could be due to planet signatures, where the role of the planetary mass is still unknown. We search for a possible planet signature by analyzing the Tc trends in the remarkable binary system HD196067-HD196068. At the moment, only HD196067 is known to host a planet which is near the brown dwarf regime. We take advantage of the strong physical similarity between both stars, which is crucial to achieving the highest possible precision in stellar parameters and elemental chemical abundances. This system gives us a unique opportunity to explore if a possible depletion of refractories in a binary system could be inhibited by the presence of a massive planet. We performed a line-by-line chemical differential study, employing the non-solar-scaled opacities, in order to reach the highest precision in the calculations. After differentially comparing both stars, HD196067 displays a clear deficiency in refractory elements in the Tc plane, a lower iron content (0.051 dex) and also a lower Li I content (0.14 dex) than its companion. In addition, the differential abundances reveal a Tc trend. These targets represent the first cases of an abundance difference around a binary system hosting a super-Jupiter. Although we explored several scenarios to explain the chemical anomalies, none of them can be entirely ruled out. Additional monitoring of the system as well as studies of larger sample of wide binary systems hosting massive planets, are needed to better understand the chemical abundance trend observed in HD196067-68.

We study line driven stellar winds using time-dependent radiation hydrodynamics where the continuum radiation couples to the gas via either a scattering or absorption opacity and there is an additional radiation force due to spectral lines that we model in the Sobolev approximation. We find that in winds with scattering opacities, instabilties tend to be suppressed and the wind reaches a steady state. Winds with absorption opacities are unstable and remain clumpy at late times. Clumps persist because they are continually regenerated in the subcritical part of the flow. Azimuthal gradients in the radial velocity distribution cause a drop in the radial radiation force and provide a mechanism for generating clumps. These clumps form on super-Sobolev scales, but at late times become Sobolev-length sized indicating that our radiation transfer model is breaking down. Inferring the clump distribution at late times therefore requires radiation-hydrodynamic modeling below the Sobolev scale.

Cen A hosts the closest active galactic nucleus to the Milky Way, which makes it an ideal target for investigating the dynamical processes in the vicinity of accreting supermassive black holes. In this paper, we present 14 Chandra HETGS spectra of the nucleus of Cen A that were observed throughout 2022. We compared them with each other, and contrasted them against the two previous Chandra HETGS spectra from 2001. This enabled an investigation into the spectral changes occurring on timescales of months and 21 years. All Chandra spectra could be well fitted by an absorbed power law with a strong and narrow Fe K$\alpha$ line, a leaked power law feature at low energies, and Si and S K$\alpha$ lines that could not be associated with the central engine. The flux of the continuum varied by a factor of $2.74\pm0.05$ over the course of the observations, whereas the Fe line only varied by $18.8\pm8.8\%$. The photon index increased over 21 years, and the Hydrogen column density varied significantly within a few months as well. The Fe K$\alpha$ line was found at a lower energy than expected from the Cen A redshift, amounting to an excess velocity of $326^{+84}_{-94}~\mathrm{km}~\mathrm{s}^{-1}$ relative to Cen A. We investigated warped accretion disks, bulk motion, and outflows as possible explanations of this shift. The spectra also featured ionized absorption lines from Fe XXV and Fe XXVI, describing a variable inflow.

The common envelope evolution (CEE) is vital in forming short orbital period compact binaries. It covers many objects, such as double compact merging binaries, type Ia supernovae progenitors, binary pulsars, and X-ray binaries. Knowledge about the common envelope (CE) eject efficiency still needs to be improved, though progress has been made recently. Short orbital period hot subdwarf B star plus white dwarf binaries are the most straightforward samples to constrain CEE physics. We apply the known orbital period-white dwarf relation to constrain the sdB progenitor of seven sdB+WD binaries with a known inclination angle. The average value of the CE efficiency parameter is 0.32, which is consistent with previous studies. However, the CE efficiency might not be a constant but is a function of the initial mass ratio based on well-constrained sdB progenitor mass and evolutionary stage. Our results can be used as physical inputs for binary population synthesis simulations on related objects. A similar method can also be applied to study other short orbital period WD binaries.

We present 1-5um spectroscopy of the coldest known brown dwarf, WISE J085510.83-071442.5 (WISE 0855), performed with the Near-Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST). NIRSpec has dramatically improved the measurement of spectral energy distribution of WISE 0855 in terms of wavelength coverage, signal-to-noise ratios, and spectral resolution. We have performed preliminary modeling of the NIRSpec data using the ATMO 2020 models of cloudless atmospheres, arriving at a best fitting model that has T_eff=285 K. That temperature is ~20 K higher than the value derived by combining our luminosity estimate with evolutionary models (i.e., the radius in the model fit to the SED is somewhat smaller than expected from evolutionary models). Through comparisons to the model spectra, we detect absorption in the fundamental band of CO, which is consistent with an earlier detection in a ground-based spectrum and indicates the presence of vertical mixing. Although PH_3 is expected in Y dwarfs that experience vertical mixing, it is not detected in WISE 0855. Previous ground-based M-band spectroscopy of WISE 0855 has been cited for evidence of H_2O ice clouds, but we find that the NIRSpec data in that wavelength range are matched well by our cloudless model. Thus, clear evidence of H_2O ice clouds in WISE 0855 has not been identified yet, but it may still be present in the NIRSpec data. The physical properties of WISE 0855, including the presence of H_2O clouds, can be better constrained by more detailed fitting with both cloudless and cloudy models and the incorporation of unpublished 5-28um data from the Mid-infrared Instrument on JWST.

This study aims to improve the photometric redshifts (photo-$z$s) of galaxies by integrating two contemporary methods: template-fitting and machine learning. Finding the synergy between these two methods was not a high priority in the past, but now that our computer processing power and observational accuracy have increased, we deem it worth investigating. We compared two methods to improve galaxy photometric redshift estimations by using the algorithms ANNz2 and BPz on different photometric and spectroscopic samples from the Sloan Digital Sky Survey (SDSS). We find that the photometric redshift performance of ANNz2 (machine learning) is better than that of BPz (galactic templates), and with the utilisation of the merging technique we introduced, we see that there is an improvement in photo-$z$ when the two strategies are consolidated, providing improvements in $\sigma_{RMS}$ and $\sigma_{68}$ up to [0.0265, 0.0222] in the LRG sample and [0.0471, 0.0471] in the Stripe-82 Sample. This simple demonstration can be used for photo-$z$s of galaxies in fainter and deeper sky surveys, and future work is required to prove its viability in these samples.

The origin and evolution of Saturn's rings is critical to understanding the Saturnian system as a whole. Here, we discuss the physical and chemical composition of the rings, as a foundation for evolutionary models described in subsequent chapters. We review the physical characteristics of the main rings, and summarize current constraints on their chemical composition. Radial trends are observed in temperature and to a limited extent in particle size distribution, with the C ring exhibiting higher temperatures and a larger population of small particles. The C ring also shows evidence for the greatest abundance of silicate material, perhaps indicative of formation from a rocky body. The C ring and Cassini Division have lower optical depths than the A and B rings, which contributes to the higher abundance of the exogenous neutral absorber in these regions. Overall, the main ring composition is strongly dominated by water ice, with minor silicate, UV absorber, and neutral absorber components. Sampling of the innermost D ring during Cassini's Grand Finale provides a new set of in situ constraints on the ring composition, and we explore ongoing work to understand the linkages between the main rings and the D ring. The D ring material is organic- and silicate-rich and water-poor relative to the main rings, with a large population of small grains. This composition may be explained in part by volatile losses in the D ring, and current constraints suggest some degree of fractionation rather than sampling of the bulk D ring material.

Investigating the temperature and density structures of gas in massive protoclusters is crucial for understanding the chemical properties therein. In this study, we present observations of the continuum and thioformaldehyde (H2CS) lines at 345 GHz of 11 massive protoclusters using the Atacama Large Millimeter/submillimeter Array (ALMA) telescope. High spatial resolution and sensitivity observations have detected 145 continuum cores from the 11 sources. H2CS line transitions are observed in 72 out of 145 cores, including line-rich cores, warm cores and cold cores. The H2 column densities of the 72 cores are estimated from the continuum emission which are larger than the density threshold value for star formation, suggesting that H2CS can be widely distributed in star-forming cores with different physical environments. Rotation temperature and column density of H2CS are derived by use of the XCLASS software. The results show the H2CS abundances increase as temperature rises and higher gas temperatures are usually associated with higher H2CS column densities. The abundances of H2CS are positively correlated with its column density, suggesting that the H2CS abundances are enhanced from cold cores, warm cores to line-rich cores in star forming regions.

We analyze the photometric data in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) over $\sim 1,200$ deg$^{2}$ to uncover new halo substructures beyond the distance, $D_{\odot}\sim$ 30 kpc, from the Sun. For this purpose, we develop an isochrone filter for an old, metal-poor stellar system to extract the faint main-sequence stars at a range of distances. With this method, we detect, not only the previously discovered substructures such as the Orphan Stream, but also the new overdensity toward Bo\"otes at 60 $\lesssim\,D_{\odot}\,\lesssim$ 100 kpc and the new stream-like feature toward Pisces at around $D_{\odot}\sim$ 60 kpc. It has been suggested that a small-scale overdensity exists in this direction of Pisces (the so-called Pisces Overdensity), but our results show that the overdensity is widely spread with a tidally elongated feature. Combining our results with the ongoing Hyper Suprime-Cam narrow-band survey and the near-future spectroscopic survey with Prime Focus Spectrograph (PFS) will allow us to place strong constraints on the origin of these halo substructures.

Here it is shown 1) how isocurvature inhomogeneities correlated on causally disconnected (super-horizon) scales are generated from curvature inhomogeneities which are known to be correlated on these scales 2) that super-horizon isocurvature generation is nearly inevitable for non-equilibrium chemical processes 3) that the amplitude of the compositional isocurvature correlations a) can be large for production of rare objects, b) falls off rapidly with separation c) falls off at scales below the horizon when these modes are generated. These two fall-offs results in an "isocurvature bump" in the power spectrum. Isocurvature generation is illustrated by the process of dark matter freeze-in, computed here with both separate universe modelling and linear perturbation theory. For freeze-in the most prominent isocurvature modes are inhomogeneities in the ratio of dark matter to standard model matter. Much smaller inhomogeneities in the ratio of baryons to standard model entropy are also produced. Previous constraints on freeze-in from Ly-$\alpha$ clouds limit the bump enhancement to $\lesssim10\%$ on comoving scales $\lesssim1\,$Mpc. Current observations are not sensitive to the isocurvature modes generated in viable freeze-in models. Results are obtained using a somewhat novel framework to describe cosmological inhomogeneities.

To understand the fundamental parameters of galaxy evolution, we investigated the minimum set of parameters that explain the observed galaxy spectra in the local Universe. We identified four latent variables that efficiently represent the diversity of high-dimensional galaxy spectral energy distributions (SEDs) observed by the Sloan Digital Sky Survey. Additionally, we constructed meaningful latent representation using conditional variational autoencoders trained with different permutations of galaxy physical properties, which helped us quantify the information that these traditionally used properties have on the reconstruction of galaxy spectra. The four parameters suggest a view that complex SED population models with a very large number of parameters will be difficult to constrain even with spectroscopic galaxy data. Through an Explainable AI (XAI) method, we found that the region below 5000\textup{\AA} and prominent emission lines ([O II], [O III], and H$\alpha$) are particularly informative for predicting the latent variables. Our findings suggest that these latent variables provide a more efficient and fundamental representation of galaxy spectra than conventionally considered galaxy physical properties.

We report physical properties of the brightest ($S_{870\,\mu \rm m}=12.4$-$19.2\,$mJy) and not strongly lensed 18 870$\,\mu$m selected dusty star-forming galaxies (DSFGs), also known as submillimeter galaxies (SMGs), in the COSMOS field. This sample is part of an ALMA band$\,$3 spectroscopic survey (AS2COSPEC), and spectroscopic redshifts are measured in 17 of them at $z=2$-$5$. We perform spectral energy distribution analyses and deduce a median total infrared luminosity of $L_{\rm IR}=(1.3\pm0.1)\times10^{13}\,L_{\odot}$, infrared-based star-formation rate of ${\rm SFR}_{\rm IR}=1390\pm150~M_{\odot}\,\rm yr^{-1}$, stellar mass of $M_\ast=(1.4\pm0.6)\times10^{11}\,M_\odot$, dust mass of $M_{\rm dust}=(3.7\pm0.5)\times10^9\,M_\odot$, and molecular gas mass of $M_{\rm gas}= (\alpha_{\rm CO}/0.8)(1.2\pm0.1)\times10^{11}\,M_\odot$, suggesting that they are one of the most massive, ISM-enriched, and actively star-forming systems at $z=2$-$5$. In addition, compared to less massive and less active galaxies at similar epochs, SMGs have comparable gas fractions; however, they have much shorter depletion time, possibly caused by more active dynamical interactions. We determine a median dust emissivity index of $\beta=2.1\pm0.1$ for our sample, and by combining our results with those from other DSFG samples, we find no correlation of $\beta$ with redshift or infrared luminosity, indicating similar dust grain compositions across cosmic time for infrared luminous galaxies. We also find that AS2COSPEC SMGs have one of the highest dust-to-stellar mass ratios, with a median of $0.02\pm0.01$, significantly higher than model predictions, possibly due to too strong of a AGN feedback implemented in the model. Finally, our complete and uniform survey enables us to put constraints on the most massive end of the dust and molecular gas mass functions.

In this Research Highlights, we summarize 31 contributions provided during the Workshop \textit{Multifrequency Behaviour of High Energy Cosmic Sources - XIV}, held in Palermo (Italy) from the 12th to the 17th of June 2023. We will start with the most recent discoveries in the field of gravitational waves (GWs). We will connect this topic to the contributions of Gamma-Ray Bursts (GRBs) associated with GWs and with the Kilonovae (KNe) hunting and, more in general, on GRBs. Continuing on high-energy astrophysics objects, we will delve into Active Galactic Nuclei (AGNs), neutrino astronomy and the study of the primordial universe, both from the space telescopes' observation and from the very recent proposals in terms of cosmological models. From the faraway universe, we will move to the more local scales and discuss the recent observations in Supernova Remnants (SNRs), massive star binaries, globular cluster dynamics, and exoplanets observed by Kepler.

FU Ori-type stars are young stellar objects (YSOs) experiencing luminosity outbursts by a few orders of magnitude, which last for $\sim$$10^2$ years. A dozen of FUors are known up to date, but many more currently quiescent YSOs could have experienced such outbursts in the last $\sim$$10^3$ years. To find observational signatures of possible past outbursts, we utilise ANDES, RADMC-3D code as well as CASA ALMA simulator to model the impact of the outburst on the physical and chemical structure of typical FU Ori systems and how it translates to the molecular lines' fluxes. We identify several combinations of molecular lines that may trace past FU Ori objects both with and without envelopes. The most promising outburst tracers from an observational perspective are the molecular flux combinations of the N$_{2}$H$^{+}$ $J=3-2$, C$^{18}$O $J = 2-1$, H$_2$CO $(J_{\rm K_a, K_c}) = 4_{04}-3_{03}$, and HCN $J = 3-2$ lines. We analyse the processes leading to molecular flux changes and show that they are linked with either thermal desorption or enhanced chemical reactions in the molecular layer. Using observed CO, HCN, N$_2$H$^+$ and H$_2$CO line fluxes from the literature, we identify ten nearby disc systems that might have undergone FU Ori outbursts in the past $\sim$$10^3$ years: [MGM2012] 556, [MGM2012] 371 and [MGM2012] 907 YSOs in L1641, Class II protoplanetary discs around CI Tau, AS 209 and IM Lup and transitional discs DM Tau, GM Aur, LkCa 15 and J1640-2130.

Since 2015 the advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) has detected a large number of gravitational wave events, originating from both binary neutron stars and binary black hole (BBH) mergers. In light of these detections, we simulate the dynamics of ambient test particles in the gravitational potential well of a BBH system close to its inspiral phase with the goal of simulating the associated electromagnetic radiation and resulting spectral energy distribution of such a BBH system. This could shed light on possible detection ranges of electromagnetic counterparts to BBH mergers. The potentials are numerically calculated using finite difference methods, under the assumption of non-rotating black holes with the post-Newtonian Paczynski-Wiita potential approximation in tandem with retarded time concepts analogous to electrodynamics. We find that the frequencies of potential electromagnetic radiation produced by these systems (possibly reaching Earth), range between a few $\text{kHz}$ to a few $100 \text{MHz}$.

Interferometric measurements of the 21cm signal are a prime example of the data-driven era in astrophysics we are entering with current and upcoming experiments. We showcase the use of deep networks that are tailored for the structure of 3D tomographic 21cm light-cones to firstly detect and characterise HI sources and to secondly directly infer global astrophysical and cosmological model parameters. We compare different architectures and highlight how 3D CNN architectures that mirror the data structure are the best-performing model.

Multi-epoch VLBI observations measure 3D water maser motions in protostellar outflows, enabling analysis of inclination and velocity. However, these analyses assume that water masers and shock surfaces within outflows are co-propagating. We compared VLBI data on maser-traced bowshocks in high-mass protostar AFGL 5142-MM1, from seven epochs of archival data from the VLBI Exploration of Radio Astrometry (VERA), obtained from April 2014 to May 2015, and our newly-conducted data from the KVN and VERA Array (KaVA), obtained in March 2016. We find an inconsistency between the expected displacement of the bowshocks and the motions of individual masers. The separation between two opposing bowshocks in AFGL 5142-MM1 was determined to be $337.17\pm0.07~\rm{mas}$ in the KaVA data, which is less than an expected value of $342.1\pm0.7~\rm{mas}$ based on extrapolation of the proper motions of individual maser features measured by VERA. Our measurements imply that the bowshock propagates at a velocity of $24\pm3~\rm{km~s^{-1}}$, while the individual masing gas clumps move at an average velocity of $55\pm5~\rm{km~s^{-1}}$, i.e. the water masers are moving in the outflow direction at double the speed at which the bowshocks are propagating. Our results emphasise that investigations of individual maser features are best approached using short-term high-cadence VLBI monitoring, while long-term monitoring on timescales comparable to the lifetimes of maser features, are better suited to tracing the overall evolution of shock surfaces. Observers should be aware that masers and shock surfaces can move relative to each other, and that this can affect the interpretation of protostellar outflows.

We report a failed solar filament eruption that involves external magnetic reconnection in a quadrupolar magnetic configuration. The evolution exhibits three kinematic evolution phases: a slow-rise phase, an acceleration phase, and a deceleration phase. In the early slow rise, extreme-ultraviolet (EUV) brightenings appear at the expected null point above the filament and are connected to the outer polarities by the hot loops, indicating the occurrence of a breakout reconnection. Subsequently, the filament is accelerated outward, accompanied by the formation of low-lying high-temperature post-flare loops ($>$ 15 MK), complying with the standard flare model. However, after 2--3 minutes, the erupting filament starts to decelerate and is finally confined in the corona. The important finding is that the confinement is closely related to an external reconnection as evidenced by the formation of high-lying large-scale hot loops ($>$ 10 MK) with their brightened footpoints at the outer polarities, the filament fragmentation and subsequent falling along the newly formed large-scale loops, as well as a hard X-ray source close to one of the outer footpoint brightenings. We propose that, even though the initial breakout reconnection and subsequent flare reconnection commence and accelerate the filament eruption, the following external reconnection between the erupting flux rope and overlying field, as driven by the upward filament eruption, makes the eruption finally failed, as validated by the numerical simulation of a failed flux rope eruption.

The origin of the radio synchrotron background (RSB) is currently unknown. Its understanding might have profound implications in fundamental physics or might reveal a new class of radio emitters. In this work, we consider the scenario in which the RSB is due to extragalactic radio sources and measure the angular cross-correlation of LOFAR images of the diffuse radio sky with matter tracers at different redshifts, provided by galaxy catalogs and CMB lensing. We compare these measured cross-correlations to those expected for models of RSB sources. We find that low-redshift populations of discrete sources are excluded by the data, while higher redshift explanations are compatible with available observations. We also conclude that at least 20\% of the RSB surface brightness level must originate from populations tracing the large-scale distribution of matter in the universe, indicating that at least this fraction of the RSB is of extragalactic origin. Future measurements of the correlation between the RSB and tracers of high-redshift sources will be crucial to constraining the source population of the RSB.

We propose that the recently observed diffuse neutrinos by IceCube with energies above 1 PeV might have originated from Sagittarius $\mathrm{A}^{\star}$ located in the galactic disk. This implies that the astrophysical settings of Sagittarius $\mathrm{A}^{\star}$ need to accelerate hadronic cosmic rays to energies of $\sim 100$ PeV or more. Then, the hadronic emission scenario argues that this galactic neutrino source is also a PeV gamma-ray source. Recent observation of galactic diffuse PeV gamma-rays with energies $\sim 1$ PeV by Large High Altitude Air Shower Observatory has also advocated this conjecture. In the present paper, we demonstrate that if protons are accelerated to energies of $\sim 100$ PeV or more as reported by Osmanov {\it et al.} (Astrophys. J. {\bf 835} 164:2017) in Sagittarius $\mathrm{A}^{\star}$ environment, then they might generate PeV neutrinos and gamma rays through cosmic rays-gas/interstellar matter ({\it e.g.} $pp$) interactions. We estimate theoretically the diffuse neutrino flux due to back-to-back charged pion decays and the accompanying gamma-ray flux from neutral pion decays. These results suggest that a fraction ($\simeq 1\%$) of the PeV diffuse neutrino flux observed by IceCube can be explained by the neutrino emission from Sagittarius $\mathrm{A}^{\star}$. Upcoming IceCube Gen2 and Cherenkov telescope array could be able to test our scenario for PeV neutrino and gamma-ray emissions from the only known galactic Pevatron Sagittarius $\mathrm{A}^{\star}$ with CR energies more than $100$ PeV.

We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z=4-6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and six galaxies from archival data, we have individually detected rest-frame 158$\mu$m dust continuum emission from 19 galaxies, nine of which are reported for the first time. The derived far-infrared luminosities are in the range $\log_{10} L_{\rm IR}\,[L_{\odot}]=$10.9-12.4, an order of magnitude lower than previously detected massive dusty star-forming galaxies (DSFGs). The average relationship between the fraction of dust-obscured star formation ($f_{\rm obs}$) and the stellar mass is consistent with previous results at z=4-6 in a mass range of $\log_{10} M_{\ast}\,[M_{\odot}]\sim$9.5-11.0 and show potential evolution from z=6-9. The individual $f_{\rm obs}$ exhibits a significant diversity, and it shows a correlation with the spatial offset between the dust and the UV continuum, suggesting the inhomogeneous dust reddening may cause the source-to-source scatter in $f_{\rm obs}$. The effective radii of the dust emission are on average $\sim$1.5 kpc and are $\sim2$ times more extended than the rest-frame UV. The infrared surface densities of these galaxies ($\Sigma_{\rm IR}\sim2.0\times10^{10}\,L_{\odot}\,{\rm kpc}^{-2}$) are one order of magnitude lower than those of DSFGs that host compact central starbursts. On the basis of the comparable contribution of dust-obscured and dust-unobscured star formation along with their similar spatial extent, we suggest that typical star-forming galaxies at z=4-6 form stars throughout the entirety of their disks.

Transiting planets orbiting M dwarfs provide the best opportunity to study the atmospheres of rocky planets with current facilities. As JWST enters its second year of science operations, an important initial endeavor is to determine whether these rocky planets have atmospheres at all. M dwarf host stars are thought to pose a major threat to planetary atmospheres due to their high magnetic activity over several billion-year timescales, and might completely strip atmospheres. Several Cycle 1 and 2 GO and GTO programs are testing this hypothesis, observing a series of rocky planets to determine whether they retained their atmospheres. A key case-study is TRAPPIST-1c, which receives almost the same bolometric flux as Venus. We might, therefore, expect TRAPPIST-1c to possess a thick, $\mathrm{CO}_2$-dominated atmosphere. Instead, Zieba et al. (2023) show that TRAPPIST-1c has little to no CO$_2$ in its atmosphere. To interpret these results, we run coupled time-dependent simulations of planetary outgassing and atmospheric escape to model the evolution of TRAPPIST-1c's atmosphere. We find that the stellar wind stripping that is expected to occur on TRAPPIST-1c over its lifetime can only remove up to $\sim 16$ bar of $\mathrm{CO}_2$, less than the modern $\mathrm{CO}_2$ inventory of either Earth or Venus. Therefore, we infer that TRAPPIST-1c either formed volatile-poor, as compared to Earth and Venus, or lost a substantial amount of $\mathrm{CO}_2$ during an early phase of hydrodynamic hydrogen escape. Finally, we scale our results for the other TRAPPIST-1 planets, finding that the more distant TRAPPIST-1 planets may readily retain atmospheres.

The FRiM fractal operator belongs to a family of operators, called ASAP, defined by an ordered selection of nearest neighbors. This generalization provides means to improve upon the good properties of FRiM. We propose a fast algorithm to build an ASAP operator mimicking the fractal structure of FRiM for pupils of any size and geometry and to learn the sparse coefficients from empirical data. We empirically show the good approximation by ASAP of correlated statistics and the benefits of ASAP for solving phase restoration problems.

I present a review of how late observations of supernovae, of the nebular phase, and much later of supernova remnants (SNRs), and their analysis in 2023 made progress towards breakthroughs in supporting the jittering jets explosion mechanism (JJEM) for core-collapse supernovae (CCSNe) and in introducing the group of lonely white dwarf (WD) scenarios for type Ia supernovae (SNe Ia). The new analyses of CCSN remnants (CCSNRs) reveal point-symmetric morphologies in a way unnoticed before in three CCSNRs. Comparison to multipolar planetary nebulae that are shaped by jets suggests that jets exploded these CCSNe, as predicted by the JJEM, but incompatible with the prediction of the delayed neutrino explosion mechanism. The spherical morphology of the ejecta Pa 30 of the historical type Iax supernova (SN Iax) of 1181 AD, which studies in 2023 revealed, is mostly compatible with the explosion of a lonely WD. Namely, at the explosion time, there is only a WD, without any close companion, although the WD was formed via a close binary interaction, i.e., binary merger. An identification of point-symmetry in SNR G1.9+0.3, a normal SN Ia and the youngest SN in the Galaxy, suggests an SN explosion of a lonely WD inside a planetary nebula (an SNIP). The group of lonely WD scenarios includes the core degenerate scenario and the double degenerate scenario with a merger to explosion delay (MED) time. SN Ia explosions of lonely WDs are common, and might actually account for most (or even all) normal SNe Ia.

Gravitational waves (GWs) induced by primordial fluctuations can be affected by the modification of the sound speed $c^2_{\rm s}$ and the equation of state parameter $w$ once the curvature fluctuations reenter the cosmological horizon. We consider a hypothetical softening of $w$ and $c^2_{\rm s}$ caused by a smooth crossover beyond Standard Model theories, for what we numerically compute the secondary induced GW considering the case of a flat scale-invariant power spectrum. We find that if the amplitude of the power spectrum is sufficiently large, the characteristic feature of the GW signal caused by the smooth crossover can be detected by future spaced-based gravitational wave interferometers and differentiated from the pure radiation case. At the same time, depending on the mass scale where the crossover occurs, such a scenario can have compatibility with primordial black holes being all the dark matter when $\mathcal{A} \sim \mathcal{O}(10^{-3})$, with a mass function very sharply peaked around the horizon mass scale of the minimum of the sound speed.

We have characterized 26 Spitzer/IRAC-selected sources from the SMUVS program that are undetected in the UltraVISTA DR5 H- and/or Ks-band images, covering 94 square arcmin of the COSMOS field which have deep multi-wavelength JWST photometry. We analyzed the JWST/NIRCam imaging from the PRIMER survey and ancillary HST data to reveal the properties of these galaxies from spectral energy distribution (SED) fitting. We find that the majority of these galaxies are detected by NIRCam at <2 micron, with only four remaining as near-infrared dropouts in the deeper JWST images. Our results indicate that the UltraVISTA dropouts candidates are primarily located at z>3 and are characterized by high dust extinctions, with a typical colour excess E(B-V) = 0.5 pm 0.3 and stellar mass log(M*/Msun) = 9.5 pm 1.0. Remarkably, ~75% of these sources show a flux enhancement between the observed photometry and modelled continuum SED that can be attributed to Halpha emission in the corresponding NIRCam bands. The derived (Halpha+ N[II] + S[II]) rest-frame equivalent widths and Halpha star formation rates (SFRs) span values ~100-2200 A and ~5-375 Msun/yr, respectively. The locations of these sources on the SFR-M* plane indicates that 35% of them are starbursts, 40% are main-sequence galaxies and the remaining 25% are located in the star-formation valley. Our sample includes one AGN and three sub-millimeter sources, as revealed from ancillary X-ray and sub-mm photometry. The high dust extinctions combined with the flux boosting from Halpha emission explain why these sources are relatively bright Spitzer galaxies and yet unidentified in the ultra-deep UltraVISTA near-infrared images.

Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial conditions and environment. Those in resonance (with their orbital periods related by a ratio of small integers) are particularly valuable because they imply a system architecture practically unchanged since its birth. Here, we present the observations of six transiting planets around the bright nearby star HD 110067. We find that the planets follow a chain of resonant orbits. A dynamical study of the innermost planet triplet allowed the prediction and later confirmation of the orbits of the rest of the planets in the system. The six planets are found to be sub-Neptunes with radii ranging from 1.94 to 2.85 Re. Three of the planets have measured masses, yielding low bulk densities that suggest the presence of large hydrogen-dominated atmospheres.

We have taken advantage of the implementation of an adaptive optics system on the Themis solar telescope to implement innovative strategies based on an inverse problem formulation for the control loop. Such an approach encompassing the whole system implies the estimation of the pixel variances of the Shack-Hartmann wavefront sensor, a novel real-time method to extract the wavefront slopes as well as their associated noise covariance, and the computation of pseudo-open loop data. The optimal commands are computed by iteratively solving a regularized inverse problem with spatio-temporal constraints including Kolmogorov statistics. The latency of the dedicated real-time control software with conventional CPU is shorter than 300 $\mu$s from the acquisition of the raw 400 x 400 pixel wavefront sensor image to the sending of the commands.

The Parker Solar Probe (PSP) provides us the unprecedentedly close approach observation to the Sun, and hence the possibility of directly understanding the "elementary process" which occurs in the kinetic scale of particles collective interactioin in solar coronal plasmas. We reported a kind of weak solar radio bursts (SRBs), which are detected by PSP when it passed a low-density magnetic channel during its second encounter phase. These weak SRBs have low starting frequecny $\sim 20$ MHz and narrow frequency range from a few tens MHz to a few hundres kHz. Their dynamic spectra display a strongly evolving feature of the intermediate relative drift rate decreasing rapidly from above 0.01/s to below 0.01/s. Analyses based on common empirical models of solar coronal plasmas indicate that these weak SRBs originate from the heliocentric distance $\sim 1.1-6.1~R_S$ (the solar radius), a typical solar wind acceleration region with a low-$\beta$ plasma, and indicate that their soruces have a typic motion velociy $\sim v_A$ (Alfv\'en velocity) obviously lower than that of fast electrons required by effectively exciting SRBs. We propose that solitary kinetic Alfv\'en waves with kinetic scales can be responsible for the generation of these small-scalevweak SRBs, called solitary wave radiation (SWR).

The Planetary Nebula Luminosity Function (PNLF) remains an important extragalactic distance indicator despite a still limited understanding of its most important feature - the bright cut-off. External galaxies benefit from consistent distance and extinction, which makes determining the PNLF easier but detailed study of individual objects much more difficult. Now, the advent of parallaxes from the Gaia mission has dramatically improved distance estimates to planetary nebulae (PNe) in the Milky Way. We have acquired ground-based narrowband imagery and measured the [OIII] fluxes for a volume-limited sample of hundreds of PNe whose best distance estimates from Gaia parallaxes and statistical methods place them within 3 kpc of the Sun. We present the first results of our study, comparing the local PNLF to other galaxies with different formation histories, and discussing how the brightness of the PNe relates to the evolutionary state of their central stars and the properties of the nebula.

In this paper, we compare the scalar field dynamics in axion-like and power law potentials for both positive and negative values of the exponents. We find that, for positive exponents, both the potentials exhibit similar scalar field dynamics and it can be difficult to distinguish them at least at the background level. Even though the potentials are oscillatory in nature for positive exponents scaling solutions can be achieved for larger values of the exponent for which the dynamics can be different during early times. Because of the presence of this scaling nature there is a turnaround in the values of the scalar field equation of state as we increase the values of the exponent in both the potentials. This indicates the deviation from the oscillatory behaviour for the larger values of the exponent. For negative values of the exponent, the dynamics of the scalar field is distinguishable and axion-like potential can give rise to cosmologically viable tracker solutions unlike the power law potentials. So, while for positive exponents we may not distinguish the two potentials for negative exponents the dynamics of the scalar field is distinguishable.

If dark matter (DM) consists of primordial black holes (PBHs) and particles simultaneously, PBHs are generically embedded within particle DM halos. Such "dressed PBHs" (dPBHs) are not subject to typical PBH constraints and can explain the DM abundance in the mass range $10^{-1} \sim 10^2 M_\odot$. We show that diffractive lensing of chirping gravitational waves (GWs) from binary mergers can not only discover, but can also identify dPBH lenses and discriminate them from bare PBHs on the event-by-event basis, with potential to uniquely establish the co-existence of subdominant PBHs and particle DM.

Human biology has a preference for left-handed chiral molecules and an outstanding question is if this is imposed through astrophysical origins. We aim to evaluate the known information about chiral molecules within astrophysical and astrochemical databases, evaluate chemical modeling accuracy, and use high-level CCSD(T) calculations to characterize propylene oxide and other oxirane variants. By comparing these computational values with past laboratory experiments, we find a 99.9% similarity. We also have put together a new database dedicated to chiral molecules and variants of chiral molecules to assist in answering this question.

We use high-resolution data from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) project to investigate the connection between circumnuclear gas reservoirs and nuclear activity in a sample of nearby galaxies. Our sample spans a wide range of nuclear activity types including radio galaxies, Seyfert galaxies, low-luminosity active galactic nuclei (AGN) and inactive galaxies. We use measurements of nuclear millimetre continuum emission along with other archival tracers of AGN accretion/activity to investigate previous claims that at, circumnuclear scales (<100 pc), these should correlate with the mass of the cold molecular gas. We find that the molecular gas mass does not correlate with any tracer of nuclear activity. This suggests the level of nuclear activity cannot solely be regulated by the amount of cold gas around the supermassive black hole (SMBH). This indicates that AGN fuelling, that drives gas from the large scale galaxy to the nuclear regions, is not a ubiquitous process and may vary between AGN type, with timescale variations likely to be very important. By studying the structure of the central molecular gas reservoirs, we find our galaxies have a range of nuclear molecular gas concentrations. This could indicate that some of our galaxies may have had their circumnuclear regions impacted by AGN feedback, even though they currently have low nuclear activity. On the other hand, the nuclear molecular gas concentrations in our galaxies could instead be set by secular processes.

Binary systems containing a compact object may exhibit periodic brightening episodes due to gravitational lensing as the compact object transits the companion star. Such "self-lensing" signatures have been detected before for white dwarf binaries. We attempt to use these signatures to identify detached stellar-mass neutron star and black hole binaries using data from the Zwicky Transient Facility (ZTF). We present a systematic search for self-lensing signals in Galactic binaries from a subset of high-cadence ZTF data taken in 2018. We identify 19 plausible candidates from the search, although because each candidate is observed to only brighten once, other origins such as stellar flares are more likely. We discuss prospects for more comprehensive future searches of the ZTF data.

S 308 is an X-ray emitting bubble that surrounds the Wolf-Rayet star WR6. The structure shines in the optical as well and is thus known as the Dolphin Nebula. Due to its large angular extent, it has been covered at only 90% with past XMM-Newton observations. Thanks to the unique dataset provided by the all-sky survey performed in X-rays by SRG/eROSITA, we can show for the first time the image of the bubble in its entire extent in this band, together with its spectral characterization. Moreover, we have tried to apply the same procedure for other wind-blown bubbles detected in the optical/IR and we searched for X-ray extended emission around them. We first analyzed the diffuse emission of S308, providing a detailed spectral analysis. We then considered a sample of 22 optical/IR selected wind-blown bubbles from a previous study based on WISE data, providing an estimate of the X-ray flux for the first time. We obtained the best fit for S308 with a two-temperature non-equilibrium plasma model (kT$_{1}=0.8_{-0.3}^{+0.8}$ keV and kT$_{2}=2_{-1}^{+3}$ keV) showing super-solar N abundance and low absorption. We did not detect any of the 22 optical/IR emitting bubbles in X-rays, but using our best fit model, we estimated the 3$\sigma$ flux upper limits for each bubble. We demonstrate the new possibility offered by SRG/eROSITA to study known wind-blown bubbles and look for other ones. A two-temperature plasma description seems to fit the data quite well for S308. Since all of the 22 bubbles studied still remain undetected by SRG/eROSITA, it is very likely that absorption effects and spatial compactness are responsible for the challenges standing in the way of detecting these bubbles in soft X-rays.

We derive distances and masses of stars from the Sloan Digital Sky Survey (SDSS) Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 17 (DR17) using simple neural networks. Training data for distances comes from Gaia EDR3, supplemented by literature distances for star clusters. For masses, the network is trained using asteroseismic masses for evolved stars and isochrone masses for main sequence stars. The models are trained on effective temperature, surface gravity, metallicity and carbon and nitrogen abundances. We found that our distance predictions have median fractional errors that range from $\approx 20\%$ at low log g and $\approx 10\%$ at higher log g with a standard deviation of $\approx 11\%$. The mass predictions have a standard deviation of $\pm 12\%$. Using the masses, we derive ages for evolved stars based on the correspondence between mass and age for giant stars given by isochrones. The results are compiled into a Value Added Catalog (VAC) called DistMass that contains distances and masses for 733901 independent spectra, plus ages for 396548 evolved stars.

We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 \lya absorption lines fitted from 12 HST STIS quasar spectra, with a total pathlength of \Delta z=2.1. We employ the machine-learning-based inference method that uses joint b-N distributions obtained from \lyaf decomposition. Our results show that the HI photoionization rates, \Gamma, are in good agreement with the recent UV background synthesis models, with \log (\Gamma/s^{-1})={-11.79}^{0.18}_{-0.15}, -11.98}^{0.09}_{-0.09}, and {-12.32}^{0.10}_{-0.12} at z=1.4, 1.2, and 1 respectively. We obtain the IGM temperature at the mean density, T_0, and the adiabatic index, \gamma, as [\log (T_0/K), \gamma]= [{4.13}^{+0.12}_{-0.10}, {1.34}^{+0.10}_{-0.15}], [{3.79}^{+0.11}_{-0.11}, {1.70}^{+0.09}_{-0.09}] and [{4.12}^{+0.15}_{-0.25}, {1.34}^{+0.21}_{-0.26}] at z=1.4, 1.2 and 1 respectively. Our measurements of T_0 at z=1.4 and 1.2 are consistent with the expected trend from z<3 temperature measurements as well as theoretical expectations that, in the absence of any non-standard heating, the IGM should cool down after HeII reionization. Whereas, our T_0 measurements at z=1 show unexpectedly high IGM temperature. However, because of the relatively large uncertainty in these measurements of the order of \Delta T_0~5000 K, mostly emanating from the limited redshift path length of available data in these bins, we can not definitively conclude whether the IGM cools down at z<1.5. Lastly, we generate a mock dataset to test the constraining power of future measurement with larger datasets. The results demonstrate that, with redshift pathlength \Delta z \sim 2 for each redshift bin, three times the current dataset, we can constrain the T_0 of IGM within 1500K. Such precision would be sufficient to conclusively constrain the history of IGM thermal evolution at z < 1.5.

The clustering of gravitational waves in luminosity distance space is emerging as a promising probe of the growth of structure. Just like for galaxies, its osbervation is subject to a number of relativistic corrections that affect the measured signal and need to be accounted for when fitting theoretical models to the data. We derive the full expression for the number count of gravitational waves in luminosity distance space, including all relativistic corrections, in LCDM and in scalar-tensor theories with luminal propagation of tensors. We investigate the importance of each relativistic effect and the detectability of the total signal by current and planned GW detectors. We consider also supernovae in luminosity distance space, highlighting the differences with gravitational waves in the case of scalar-tensor theories. We carry out a thorough comparison among the number count of gravitational waves and supernovae in luminosity distance space, and that of galaxies in redshift space. We show how the relativistic corrections contain useful complementary information on the growth of perturbations and on the underlying theory of gravity, highlighting the synergy with other cosmological probes.

About a dozen exoplanetary systems have been discovered with three or more planets participating in a sequence of mean-motion resonances. The unique and complex architectures of these so-called "resonant chains" motivate efforts to characterize their planets holistically. In this work, we perform a comprehensive exploration of the spin-axis dynamics of planets in resonant chains. Planetary spin states are closely linked with atmospheric dynamics and habitability and are thus especially relevant to resonant chains like TRAPPIST-1, which hosts several temperate planets. Considering a set of observed resonant chains, we calculate the equilibrium states of the planetary axial tilts ("obliquities"). We show that high obliquity states exist for $\sim60\%$ of planets in our sample, and many of these states can be stable in the presence of tidal dissipation. Using case studies of two observed systems (Kepler-223 and TOI-1136), we demonstrate how these high obliquity states could have been attained during the initial epoch of disk-driven orbital migration that established the resonant orbital architectures. We show that the TRAPPIST-1 planets most likely have zero obliquities, with the possible exception of planet d. Overall, our results highlight that both the orbital and spin states of resonant chains are valuable relics of the early stages of planet formation and evolution.

The purpose of FFREE - the new optical bench devoted to experiments on high-contrast imaging at LAOG - consists in the validation of algorithms based on off-line calibration techniques and adaptive optics (AO) respectively for the wavefront measurement and its compensation. The aim is the rejection of the static speckles pattern arising in a focal plane after a diffraction suppression system (based on apodization or coronagraphy) by wavefront pre-compensation. To this aim, FFREE has been optimized to minimize Fresnel propagation over a large near infrared (NIR) bandwidth in a way allowing efficient rejection up to the AO control radius, it stands then as a demonstrator for the future implementation of the optics that will be common to the scientific instrumentation installed on EPICS.

We compute bias, variance, and approximate confidence intervals for the efficiency of a random selection process under various special conditions that occur in practical data analysis. We consider the following cases: a) the number of trials is not constant but drawn from a Poisson distribution, b) the samples are weighted, c) the numbers of successes and failures have a variance which exceeds that of a Poisson process, which is the case, for example, when these numbers are obtained from a fit to mixture of signal and background events. Generalized Wilson intervals based on these variances are computed, and their coverage probability is studied. The efficiency estimators are unbiased in all considered cases, except when the samples are weighted. The standard Wilson interval is also suitable for case a). For most of the other cases, generalized Wilson intervals can be computed with closed-form expressions.

Barlow and Beeston presented an exact likelihood for the problem of fitting a composite model consisting of binned templates obtained from Monte-Carlo simulation which are fitted to equally binned data. Solving the exact likelihood is technically challenging, and therefore Conway proposed an approximate likelihood to address these challenges. In this paper, a new approximate likelihood is derived from the exact Barlow-Beeston one. The new approximate likelihood and Conway's likelihood are generalized to problems of fitting weighted data with weighted templates. The performance of estimates obtained with all three likelihoods is studied on two toy examples: a simple one and a challenging one. The performance of the approximate likelihoods is comparable to the exact Barlow-Beeston likelihood, while the performance in fits with weighted templates is better. The approximate likelihoods evaluate faster than the Barlow-Beeston one when the number of bins is large.

The quantum extension of the Kruskal spacetime indicates the existence of a companion black hole in the universe earlier than ours. It is shown that the radiations from the companion black hole can enter its horizon, pass through the deep Planck region, and show up from the white hole in our universe. These radiations inlay extra bright rings in the image of the black hole in our universe, and some of these rings appear distinctly in the shadow region. Therefore, the image of the black hole observed by us encodes the information of quantum gravity. The positions and widths of the bright rings are predicted precisely. The predictive values for supermassive black holes are universal for a quite general class of quantum-modified spacetimes with the phenomenon of black hole to white hole transition. Thus, our result opens a new experimental window to test this phenomenon predicted by quantum gravity.

Wide binary stars are used to test the modified gravity called Scalar-Tensor-Vector Gravity or MOG. This theory is based on the additional gravitational degrees of freedom, the scalar field $G=G_N(1+\alpha)$, where $G_N$ is Newton's constant, and the massive (spin-1 graviton) vector field $\phi_\mu$. The wide binaries have separations of 2-30 kAU. The MOG acceleration law, derived from the MOG field equations and equations of motion of a massive test particle for weak gravitational fields, depends on the enhanced gravitational constant $G=G_N(1+\alpha)$ and the effective running mass $\mu$. The magnitude of $\alpha$ depends on the physical length scale or averaging scale $\ell$ of the system. The modified MOG acceleration law for weak gravitational fields predicts that for the solar system and for the wide binary star systems gravitational dynamics follows Newton's law.

The gravitational wave (GW) spectrum at frequencies above a kHz is a largely unexplored frontier. We show that detectors with sensitivity to single-phonon excitations in crystal targets can search for GWs with frequencies, $\mathrm{THz} \lesssim f \lesssim 100 \, \mathrm{THz}$, corresponding to the range of optical phonon energies, $\mathrm{meV} \lesssim \omega \lesssim 100 \, \mathrm{meV}$. Such detectors are already being built to search for light dark matter (DM), and therefore sensitivity to high-frequency GWs will be achieved as a byproduct. We begin by deriving the absorption rate of a general GW signal into single phonons. We then focus on carefully defining the detector sensitivity to monochromatic and chirp signals, and compute the detector sensitivity for many proposed light DM detection targets. The detector sensitivity is then compared to the signal strength of candidate high-frequency GW sources, e.g., superradiant annihilation and black hole inspiral, as well as other recent detector proposals in the $\mathrm{MHz} \lesssim f \lesssim 100 \, \mathrm{THz}$ frequency range. With a judicious choice of target materials, a collection of detectors could optimistically achieve sensitivities to monochromatic signals with $h_0 \sim 10^{-23} - 10^{-25}$ over $\mathrm{THz} \lesssim f \lesssim 100 \, \mathrm{THz}$.

Stable dark matter particles may arise as pseudo-Goldstone bosons from the confinement of dark quarks interacting via a non-Abelian gauge force. Their relic abundance is determined not by annihilations into visible particles but by dark pion number-changing processes within the dark sector, such as $3 \pi_D \to 2 \pi_D$. However, if the dark vector mesons $\rho_D$ are light enough for $3 \pi_D \to \pi_D \rho_D$ annihilations to be kinematically allowed, this process dominates and significantly delays freeze-out. As a result, the preferred dark matter mass scale increases and bounds from the Bullet Cluster can be evaded.

We describe a simple and efficient Python code to perform Bayesian forecasting for gravitational waves (GW) produced by Extreme-Mass-Ratio-Inspiral systems (EMRIs). The code runs on GPUs for an efficient parallelised computation of thousands of waveforms and sampling of the posterior through a Markov-Chain-Monte-Carlo (MCMC) algorithm. EMRI_MC generates EMRI waveforms based on the so--called kludge scheme, and propagates it to the observer accounting for cosmological effects in the observed waveform due to modified gravity/dark energy. Extending the code to more accurate schemes for the generation of the waveform is straightforward. Despite the known limitations of the kludge formalism, we believe that the code can provide a helpful resource for the community working on forecasts for interferometry missions in the milli-Hz scale, predominantly, the satellite-mission LISA.

The direct detection of gravitational waves by LIGO has heralded a new era for astronomy and physics. Typically the gravitational waves observed by LIGO are dominated by noise. In this work we use Deep Convolutional Neural Networks (specifically U-Nets) to filter a clean signal from noisy data. We present two realizations of U-Net filters, the Noise2Clean U-Net filter which is trained using noisy and clean realizations of the same signal, as well as Noise2Noise U-Net which is trained on two separate noisy realization of the same signal. We find that the U-Nets successfully filter signal from noise. We also benchmark the performance of U-Nets by using them to detect the binary presence or absence of gravitational wave signals in data.

Thermal noise of the dielectric mirror coatings can limit laser-optical high-precision measurements. Coatings made of amorphous silicon and silicon nitride could provide a remedy for both gravitational-wave detectors and optical clocks. However, the absorption spectra of these materials require laser wavelengths around 2 $\mu$m. For GW detectors, ultra-stable laser light of tens or hundreds of watts is needed. Here, we report the production of nearly 30 W of ultra-stable laser light at 2128 nm by frequency conversion of 1064 nm light from a master oscillator power amplifier system. We achieve an external conversion efficiency of (67.5 $\pm$ 0.5) % via optical parametric oscillation and a relative power noise in the range of $10^{-6}$/$\sqrt{\text{Hz}}$ at 100 Hz, which is almost as low as that of the input light and underlines the potential of our approach.

pynucastro is a python library that provides visualization and analyze techniques to classify, construct, and evaluate nuclear reaction rates and networks. It provides tools that allow users to determine the importance of each rate in the network, based on a specified list of thermodynamic properties. Additionally, pynucastro can output a network in C++ or python for use in simulation codes, include the AMReX-Astrophysics simulation suite. We describe the changes in pynucastro since the last major release, including new capabilities that allow users to generate reduced networks and thermodynamic tables for conditions in nuclear statistical equilibrium.

A consequence of QCD axion dark matter being born after inflation is the emergence of ultra-small-scale substructures known as miniclusters. Although miniclusters merge to form minihalos, this intrinsic granularity is expected to remain imprinted on small scales in our galaxy, leading to potentially damning consequences for the campaign to detect axions directly on Earth. This picture, however, is modified when one takes into account the fact that encounters with stars will tidally strip mass from the miniclusters, creating pc-long tidal streams that act to refill the dark matter distribution. Here we ask whether or not this stripping rescues experimental prospects from the worst-case scenario in which the majority of axions remain bound up in unobservably small miniclusters. We find that the density sampled by the Earth on mpc-scales will be, on average, around 70-90% of the known local DM density, and at a typical point in the solar neighbourhood, we expect most of the dark matter to be comprised of debris from $\mathcal{O}(10^2$-$10^3)$ overlapping streams. If haloscopes can measure the axion signal with high-enough frequency resolution, then these streams are revealed in the form of an intrinsically spiky lineshape, in stark contrast with the standard assumption of a smooth, featureless Maxwellian distribution -- a unique prediction that constitutes a way for experiments to distinguish between pre and post-inflationary axion cosmologies.

Mid-frequency band gravitational-wave detectors will be complementary for the existing Earth-based detectors (sensitive above 10 Hz or so) and the future space-based detectors such as LISA, which will be sensitive below around 10 mHz. A ground-based superconducting omnidirectional gravitational radiation observatory (SOGRO) has recently been proposed along with several design variations for the frequency band of 0.1 to 10 Hz. For three conceptual designs of SOGRO (e.g., pSOGRO, SOGRO and aSOGRO), we examine their multi-channel natures, sensitivities and science cases. One of the key characteristics of the SOGRO concept is its six detection channels. The response functions of each channel are calculated for all possible gravitational wave polarizations including scalar and vector modes. Combining these response functions, we also confirm the omnidirectional nature of SOGRO. Hence, even a single SOGRO detector will be able to determine the position of a source and polarizations of gravitational waves, if detected. Taking into account SOGRO's sensitivity and technical requirements, two main targets are most plausible: gravitational waves from compact binaries and stochastic backgrounds. Based on assumptions we consider in this work, detection rates for intermediate-mass binary black holes (in the mass range of hundreds up to $10^{4}$ $M_\odot$) are expected to be $0.0014-2.5 \,\, {\rm yr}^{-1}$. In order to detect stochastic gravitational wave background, multiple detectors are required. Two aSOGRO detector networks may be able to put limits on the stochastic background beyond the indirect limit from cosmological observations.

A long-standing problem within the study of cosmic inflation consists in fully reconciling the stochastic approach with perturbation theory. A complete connection between both formalisms has remained elusive even in the simple case of a single scalar field with self interactions determined by an arbitrary potential, in a fixed de Sitter background with a constant expansion rate. Using perturbation theory, we offer an exact calculation of the one-point probability density function for primordial fluctuations, valid to first order in the potential. We examine under which conditions our solution respects the Fokker-Planck equation encountered within the stochastic approach. We identify discrepancies and elucidate their origins, allowing us to shed light on the validity of the stochastic formalism.

It has recently been shown that a significant slowdown of many stars can be attributed to the emergence of a strong magnetic field within the radiative region, where heat is transferred through radiation in a stably stratified layer. Here, we describe how this transition can be understood as a subcritical bifurcation to small-scale turbulence in linearly stable flows. The turbulence is sustained by a nonlinear mean-field dynamo and can be observed down to relatively small differential rotation, arbitrarily far from the linear onset of any hydrodynamic instability. In this regime, turbulent fluctuations provide diffusivity-free transfer of angular momentum that increases the transport generated by the magnetic field triggering the turbulence. Finally, we present a simple nonlinear model that captures this scenario and can be used as a general description of the transition to turbulence in astrophysical flows, as long as it involves a competition between a large-scale dynamo, and a small-scale magnetic instability.

Deviations from General Relativity can alter the quasi-normal mode (QNM) ringdown of perturbed black holes. It is known that a shift-symmetric (hence massless) scalar can only introduce black hole hair if it couples to the Gauss-Bonnet invariant, in which case the scalar charge is fixed with respect to the black hole mass and controlled by the strength of that coupling. The charge per unit mass decreases with the mass and can, therefore, be used as a perturbative parameter for black holes that are sufficiently large with respect to the scale suppressing the deviation from General Relativity or the Standard model. We construct an effective field theory scheme for QNMs using this perturbative parameter to capture deviations from Kerr for both the background and the perturbations. We demonstrate that up to second order in the charge per unit mass, QNMs can be calculated by solving standard linearised perturbation equations for the Kerr metric with sources depending on solutions of the same equations up to first order. It follows that corrections to the QNM frequencies are heavily suppressed for sufficiently massive black holes, meaning that LISA is very unlikely to detect any evidence of scalar hair in ringdown signals.

It was recently proposed that five-dimensional inflation can relate the causal size of the observable universe to the present weakness of gravitational interactions by blowing up an extra compact dimension from the microscopic fundamental length of gravity to a large size in the micron range, as required in the Dark Dimension proposal. Here, we compute the power spectrum of all primordial fluctuations emerging from a 5-dimensional inflaton in a slow-roll region of its potential, showing an interesting change of behaviour at large scales corresponding to angles larger than about 10 degrees in the sky.

We show that the maximal frequency of cosmic gravitons must not exceed the THz domain. From a classical viewpoint, both in conventional inflationary scenarios and in bouncing models the largest frequency of the spectrum overshoots the MHz band even if its specific signature is model dependent. According to a quantum mechanical perspective the maximal frequency is instead associated with the range of energies where a single pair of gravitons with opposite (comoving) three-momenta is produced. The upper limit on the largest frequency determines the minimal chirp amplitude [typically ${\mathcal O}(10^{-32})$] required for a direct detection of a cosmic signal in the THz band. Below this limiting frequency the minimal chirp amplitude can be enhanced so that the optimal range ultimately depends on the physical properties of the diffuse backgrounds. In case a hypothetical instrument (at present just a figment of a hopeful imagination) would reach chirp amplitudes down to ${\mathcal O}(10^{-30})$ in the MHz or GHz bands, the Bose-Einstein correlations could be used to probe the properties of cosmic gravitons and their super-Poissonian statistics.

In this Letter, we propose to detect the interaction of a hypothetical coherently evolving cosmological scalar field with an orbital network of quantum sensors, focusing on the GPS satellite network as a test example. Cosmological scenarios, such as a scalar-tensor theory for dark energy or the axi-Higgs model, suggest that such a field may exist. As this field would be (approximately) at rest in the CMB frame, it would exhibit a dipole as a result of the movement of our terrestrial observers relative to the CMB. While the current sensitivity of the GPS network is insufficient to detect a cosmological dipole, future networks of quantum sensors on heliocentric orbits, using state-of-the-art atomic clocks, can reach and exceed this requirement.