Papers for Monday, Aug 28 2023

If $z > 6$ quasars reside in rare, massive haloes, $\Lambda$CDM cosmology predicts they should be surrounded by an anomalously high number of bright companion galaxies. In this paper, I show that these companion galaxies should also move unusually fast. Using a new suite of cosmological, `zoom-in' hydrodynamic simulations, I here present predictions for the velocity distribution of quasar companion galaxies and its variation with quasar host halo mass at $z \, = \, 6$. Satellites accelerate as they approach the quasar host galaxy, producing a line-of-sight velocity profile that broadens as the distance to the quasar host galaxy decreases. This increase in velocity dispersion is particularly pronounced if the host halo mass is $\gtrsim 5 \times 10^{12} \, \rm M_\odot$. In this case, typical line-of-sight speeds rise to $\approx 500 \, \rm km \, s^{-1}$ at projected radii $\sim 10 \, \rm kpc$. For about $10\%$ of satellites, they should exceed $800 \, \rm km \, s^{-1}$, with $\approx 5\%$ of companions reaching line-of-sight speeds $\sim 1000 \, \rm km \, s^{-1}$. For lower host halo masses $\approx 5 \times 10^{11} - 10^{12} \, \rm M_\odot$, the velocity profile of companion galaxies is significantly flatter. In this case, typical line-of-sight velocities are $\approx 250 \, \rm km \, s^{-1}$ and do not exceed $\approx 500 \, \rm km \, s^{-1}$. A comparison with existing ALMA, JWST and MUSE line-of-sight velocity measurements reveals that observed $z > 6$ quasar companions closely follow the velocity distribution expected for a host halo with mass $\gtrsim 5 \times 10^{12} \, \rm M_\odot$, ruling out a light host halo. Finally, through an estimate of UV and [OIII] luminosity functions, I show that the velocity distribution more reliably discriminates between halo mass than companion number counts, which are strongly affected by cosmic variance.

We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 19.9 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals a hot (T >= 40,000 K) continuum and carbon features observed ~1 day after discovery, SN Ic-like photospheric-phase spectra, and strong forbidden calcium emission starting 38 days after discovery. SN 2022oqm has a relatively high peak luminosity (MB = -17 mag) for CaRTs, making it an outlier in the population. We determine that three power sources are necessary to explain SN 2022oqm's light curve, with each power source corresponding to a distinct peak in its light curve. The first peak of the light curve is powered by an expanding blackbody with a power law luminosity, consistent with shock cooling by circumstellar material. Subsequent peaks are powered by a double radioactive decay model, consistent with two separate sources of photons diffusing through an optically thick ejecta. From the optical light curve, we derive an ejecta mass and 56Ni mass of ~0.89 solar masses and ~0.09 solar masses, respectively. Detailed spectroscopic modeling reveals ejecta that is dominated by intermediate-mass elements, with signs that Fe-peak elements have been well-mixed. We discuss several physical origins for SN 2022oqm and favor a white dwarf progenitor model. The inferred ejecta mass points to a surprisingly massive white dwarf, challenging models of CaRT progenitors.

The first stars are expected to form through molecular-hydrogen (H$_2$) cooling, a channel that is especially sensitive to the thermal and ionization state of gas, and can thus act as a probe of exotic energy injection from decaying or annihilating dark matter (DM). Here, we use a toy halo model to study the impact of DM-sourced energy injection on the H$_2$ content of the first galaxies, and thus estimate the threshold mass required for a halo to form stars at high redshifts. We find that currently allowed DM models can significantly change this threshold, producing both positive and negative feedback. In some scenarios, the extra heating of the gas raises the halo mass required for collapse, whereas in others, energy injection lowers the threshold by increasing the free-electron fraction and catalyzing H$_2$ formation. The direction of the effect can be redshift-dependent. We also bracket the uncertainties from self-shielding of halos from Lyman-Werner radiation. Hence, exotic energy injection can both delay and accelerate the onset of star formation; we show how this can impact the timing of 21cm signals at cosmic dawn. We encourage detailed simulation follow-ups in the most promising regions of parameter space identified in this work.

We present the results obtained from follow-up observations of the MUSE Ultra Deep Field (MUDF) at X-ray energies with XMM-Newton. The MUDF is centred on a unique field with two bright, physically associated quasars at $z\simeq3.23$, separated by $\sim$500 kpc in projection. Both quasars are embedded within extended Ly$\alpha$ nebulae ($\gtrsim 100~\rm kpc$ at a surface brightness flux level of $\approx 6\times 10^{-19} \rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$), whose elongated morphology is suggestive of an extended filament connecting the quasar haloes. The new X-ray observations presented here allow us to characterise the physical properties (e.g. X-ray slope, luminosities, gas column densities) in the innermost region of the MUDF quasars. We find that both quasars are X-ray underluminous compared to objects at similar ultraviolet luminosities. Based on our X-ray spectral analysis, absorbing columns of $N_H(z)\gtrsim$ 10$^{23}$ cm$^{-2}$ appear unlikely, therefore such a weakness is possibly intrinsic. When also including literature data, we do not observe any detectable trend between the area of the nebulae and nuclear luminosities at both the rest-frame 2 keV and 2500 $\rm \mathring{A}$. The area is also not correlated with the X-ray photon index nor with the integrated band flux in the hard band (2$-$10 keV). We also do not find any trend between the extended Ly$\alpha$ emission of the nebulae and the nuclear X-ray luminosity. Finally, the properties of the MUDF quasars' nebulae are consistent with the observed relation between the Ly$\alpha$ integrated luminosity of the nebulae and their area. Our results suggest that the quasar ionization power is not a strong driver of the morphology and size of the nebulae.

Several types of energetic supernovae, such as superluminous supernovae (SLSNe) and broad-line Ic supernovae (Ic-BL SNe), could be powered by the spin-down of a rapidly rotating magnetar. Currently, most models used to infer the parameters for potential magnetar-driven supernovae make several unsuitable assumptions that likely bias the estimated parameters. In this work, we present a new model for magnetar-driven supernovae that relaxes several of these assumptions and an inference workflow that enables accurate estimation of parameters from lightcurves of magnetar-driven supernovae. In particular, in this model, we include the dynamical evolution of the ejecta, coupling it to the energy injected by the magnetar itself while also allowing for non-dipole spin down. We show that the model can reproduce SLSN and Ic-BL SN light curves consistent with the parameter space from computationally expensive numerical models. We also show the results of parameter inference on four well-known example supernovae, demonstrating the model's effectiveness at capturing the considerable diversity in magnetar-driven supernova lightcurves. The model fits each light curve well and recovers parameters broadly consistent with previous works. This model will allow us to explore the full diversity of magnetar-driven supernovae under one theoretical framework, more accurately characterize these supernovae from only photometric data, and make more accurate predictions of future multiwavelength emission to test the magnetar-driven scenario better.

Current gravitational-wave data from stellar-mass black-hole binary mergers suggest a correlation between the binary mass ratio $q$ and the effective spin $\chi_\mathrm{eff}$: more unequal-mass binaries consistently show larger and positive values of the effective spin. Multiple generations of black-hole mergers in dense astrophysical environments may provide a way to form unequal-mass systems, but they cannot explain the observed correlation on their own. We show that the symmetry of the astrophysical environment is a crucial feature to shed light on this otherwise puzzling piece of observational evidence. We present a toy model that reproduces, at least qualitatively, the observed correlation. The model relies on axisymmetric, disk-like environments where binaries participating in hierarchical mergers share a preferential direction. Migration traps in AGN disks are a prime candidate for this setup, hinting at the exciting possibility of constraining their occurrence with gravitational-wave data.

We present a systematic analysis of the X-ray emission of a sample of 17 optically selected, X-ray-detected tidal disruption events (TDEs) discovered between 2014 and 2021. The X-ray light curves show a diverse range of temporal behaviors, with most sources not following the expected power-law decline. The X-ray spectra are mostly extremely soft and consistent with thermal emission from the inner region of an accretion disk that cools as the accretion rate decreases. Three sources show the formation of a hard X-ray corona around 200 days after the UV/optical peak. The shape of the spectral energy distribution, traced by the ratio ($L_{\rm BB}/L_{\rm X}$) between the UV/optical and X-ray luminosities, shows a wide range $L_{\rm BB}/L_{\rm X} \in (0.5, 3000)$ at early-times, and converges to disk-like values $L_{\rm BB}/L_{\rm X} \in (0.5, 10)$ at late-times. The evolution of the derived physical parameters favors a decrease in the optical depth of a reprocessing layer instead of delayed disk formation to explain the late-time X-ray brightening found in several sources. We estimate the fraction of optically discovered TDEs with $L_{\rm X}\geq 10^{42}$ erg s$^{-1}$ to be at least $40\%$, and find that the X-ray loudness is independent of black hole mass. We combine our sample with those from X-ray surveys to construct an X-ray luminosity function, best fitted by a broken power-law with a brake at $\sim 10^{44}$ erg s$^{-1}$. We show that there is no dichotomy between optically and X-ray selected TDEs; instead, there is a continuum of early time $L_{\rm BB}/L_{\rm X}$, at least as wide as $L_{\rm BB}/L_{\rm X} \in (0.1, 3000)$, with optical/X-ray surveys selecting preferentially, but not exclusively, from the higher/lower end of the distribution. Our findings are consistent with an orientation-dependent and time-evolving reprocessing layer, and support viewing-angle unification models.

Populations of stellar-mass black holes (BHs) in globular clusters (GCs) influence their dynamical evolution and have important implications on one of the main formation channels for gravitational wave sources. Inferring the size of these populations remains difficult, however. In this work, multimass models of 34 Milky Way GCs, first presented in Dickson et al., are used to explore the present-day BH populations of a large sample of clusters. Direct constraints on both the total and visible mass components provided by several observables allow these models to accurately determine the distribution of the dark mass (including BHs) within clusters, as we demonstrate in a proof-of-concept fitting of the models to mock observations extracted from Monte Carlo cluster models. New constraints on the BH population retained to the present-day in each cluster are inferred from our models. We find that BH mass fractions ranging from 0 to 1 per cent of the total mass are typically required to explain the observations, except for Omega Cen, for which we infer a mass fraction of 5 per cent, in agreement with previous works. Relationships between the dark remnant populations and other cluster parameters are examined, demonstrating a clear anti-correlation between the amount of BHs and mass segregation between visible stars, as well as a correlation between remnant mass fractions and the dynamical age of clusters. Our inferred BH populations are in good agreement overall with other recent studies using different methodologies, but with notable discrepancies for individual clusters.

The demographics of young exoplanets can shed light onto their formation and evolution processes. Exoplanet properties are derived from the properties of their host stars. As such, it is important to accurately characterize the host stars since any systematic biases in their derivation can negatively impact the derivation of planetary properties. Here, we present a uniform catalog of photometrically-derived stellar effective temperatures, luminosities, radii, and masses for 4,865 young (<1 Gyr) stars in 31 nearby clusters and moving groups within 200 pc. We compared our photometrically-derived properties to a subset of those derived from spectra, and found them to be in good agreement. We also investigated the effect of stellar properties on the detection efficiency of transiting short-period young planets with TESS as calculated in Fernandes et al. 2022, and found an overall increase in the detection efficiency when the new photometrically derived properties were taken into account. Most notably, there is a 1.5 times increase in the detection efficiencies for sub-Neptunes/Neptunes (1.8-6 Re) implying that, for our sample of young stars, better characterization of host star properties can lead to the recovery of more small transiting planets. Our homogeneously derived catalog of updated stellar properties, along with a larger unbiased stellar sample and more detections of young planets, will be a crucial input to the accurate estimation of the occurrence rates of young short-period planets.

Coronal holes are areas on the Sun with open magnetic field lines. They are a source region of the solar wind, but how the wind emerges from coronal holes is not known. We observed a coronal hole using the Extreme Ultraviolet Imager on the Solar Orbiter spacecraft. We identified jets on scales of a few hundred kilometers, which last 20 to 100 seconds and reach speeds of ~100 kilometers per second. The jets are powered by magnetic reconnection and have kinetic energy in the picoflare range. They are intermittent but widespread within the observed coronal hole. We suggest that such picoflare jets could produce enough high-temperature plasma to sustain the solar wind and that the wind emerges from coronal holes as a highly intermittent outflow at small scales.

Observed scatter in the Lyman-alpha opacity of quasar sightlines at $z<6$ has motivated measurements of the correlation between Ly$\alpha$ opacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Ly$\alpha$ troughs at $z\sim5.7$ with a deficit of Lyman-$\alpha$ emitting galaxies (LAEs). In this work, we survey two of the most highly transmissive lines of sight at this redshift, towards the $z=6.02$ quasar SDSS J1306+0356 and the $z=6.17$ quasar PSO J359-06. We find that both fields are underdense in LAEs within 10 $h^{-1}$ Mpc of the quasar sightline, somewhat less extensive than underdensities associated with Ly$\alpha$ troughs. We combine our observations with three additional fields from the literature, and find that while fields with extreme opacities are generally underdense, moderate opacities span a wider density range. The results at high opacities are consistent with models that invoke UV background fluctuations and/or late reionization to explain the observed scatter in IGM Ly$\alpha$ opacities. There is tension at low opacities, however, as the models tend to associate lower IGM Ly$\alpha$ opacities with higher densities. Although the number of fields surveyed is still small, the low-opacity results may support a scenario in which the ionizing background in low-density regions increases more rapidly than some models suggest after becoming ionized. Elevated gas temperatures from recent reionization may also be making these regions more transparent.

The upcoming Comet Interceptor mission involves a parking phase around the Sun-Earth L2 point before transferring to intercept the orbit of a long period comet, interstellar object or a back-up target in the form of a short-period comet. The target is not certain to be known before the launch in 2029. During the parking phase there may thus arise a scenario wherein a decision needs to be taken of whether to go for a particular comet or whether to discard that option in the hope that a better target will appear within a reasonable time frame later on. We present an expectation value-based formalism that could aid in the associated decision making provided that outlined requirements for its implementation exist.

We present a detailed analysis of nearly two decades of optical/UV and X-ray data to study the multi-wavelength pre-explosion properties and post-explosion X-ray properties of nearby SN2023ixf located in M101. We find no evidence of precursor activity in the optical to UV down to a luminosity of $\lesssim 7\times10^{4}\, \rm L_{\odot}$, while X-ray observations covering nearly 18 years prior to explosion show no evidence of luminous precursor X-ray emission down to an absorbed 0.3 - 10.0 keV X-ray luminosity of $\sim6\times10^{36}$ erg s$^{-1}$. Extensive Swift observations taken post-explosion did not detect soft X-ray emission from SN2023ixf within the first $\sim$3.3 days after first light, which suggests a mass-loss rate for the progenitor of $\lesssim5\times10^{-4}\,\rm M_{\odot}$ yr$^{-1}$ or a radius of $\lesssim4\times10^{15}$ cm for the circumstellar material. Our analysis also suggests that if the progenitor underwent a mass-loss episode, this had to occur $>$ 0.5 - 1.5 years prior to explosion, consistent with previous estimates. Swift detected soft X-rays from SN2023ixf $\sim4.25$ days after first light, and it rose to a peak luminosity of $\sim10^{39}$ erg s$^{-1}$ after 10 days and has maintained this luminosity for nearly 50 days post first light. This peak luminosity is lower than expected, given the evidence that SN2023ixf is interacting with dense material. However, this might be a natural consequence of an asymmetric circumstellar medium. X-ray spectra derived from merging all Swift observations over the first 50 days are best described by a two-component bremsstrahlung model consisting of a heavily absorbed and hotter component similar to that found using NuSTAR, and a less-absorbed, cooler component. We suggest that this soft component arises from cooling of the forward shock similar to that found in Type IIn SN2010jl.

Hydrogen loss to space is a key control on the evolution of the Martian atmosphere and the desiccation of the red planet. Thermal escape is thought to be the dominant loss process, but both forward modeling studies and remote sensing observations have indicated the presence of a second, higher-temperature "nonthermal" or "hot" hydrogen component, some fraction of which also escapes. Exothermic reactions and charge/momentum exchange processes produce hydrogen atoms with energy above the escape energy, but H loss via many of these mechanisms has never been studied, and the relative importance of thermal and nonthermal escape at Mars remains uncertain. Here we estimate hydrogen escape fluxes via 47 mechanisms, using newly-developed escape probability profiles. We find that HCO$^+$ dissociative recombination is the most important of the mechanisms, accounting for 30-50% of the nonthermal escape. The reaction CO$_2^+$ + H$_2$ is also important, producing roughly as much escaping H as momentum exchange between hot O and H. Total nonthermal escape from the mechanisms considered amounts to 39% (27%) of thermal escape, for low (high) solar activity. Our escape probability profiles are applicable to any thermospheric hot H production mechanism and can be used to explore seasonal and longer-term variations, allowing for a deeper understanding of desiccation drivers over various timescales. We highlight the most important mechanisms and suggest that some may be important at Venus, where nonthermal escape dominates and much of the literature centers on charge exchange reactions, which do not result in significant escape in this study.

Recent numerical simulations have revealed that dust clumping and planetesimal formation likely proceed in ring-like disc substructures, where dust gets trapped in weakly turbulent pressure maxima. The streaming instability has difficulty operating in such rings with external turbulence and no pressure gradient. To explore potential paths to planetesimal formation in this context, we analyse the stability of turbulent dust-trapping ring under the shearing sheet framework. We self-consistently establish the pressure maximum and the dust ring in equilibrium, the former via a balance of external forcing versus viscosity and the latter via dust drift versus turbulent diffusion. We find two types of $\gtrsim H$-scale instabilities ($H$ being the pressure scale height), which we term the dusty Rossby wave instability (DRWI). Type I is generalised from the standard RWI, which is stationary at the pressure maximum and dominates in relatively sharp pressure bumps. Type II is a newly identified travelling mode that requires the presence of dust. It can operate in relatively mild bumps, including many that are stable to the standard RWI, and its growth rate is largely determined by the equilibrium gas and dust density gradients. We further conduct two-fluid simulations that verify the two types of the DRWI. While Type I leads strong to dust concentration into a large gas vortex similar to the standard RWI, the dust ring is preserved in Type II, and meanwhile exhibiting additional clumping within the ring. The DRWI suggests a promising path towards formation of planetesimals/planetary embryos and azimuthally asymmetric dust structure from turbulent dust-trapping rings.

For the first time, systematic studies of dwarf galaxies are being conducted throughout the Local Volume, including Centaurus A (NGC 5128), which is the nearest elliptical galaxy. Given Centaurus As mass (roughly ten times that of the Milky Way), AGN activity, and recent major mergers, investigating these dwarfs and their star formation physics is imperative. However, simulating the faintest dwarfs in a massive galaxy like Centaurus A with sufficient resolution in a hydrodynamic simulation is computationally expensive and currently unfeasible. In this study, we seek to reproduce Centaurus A dwarfs using the same star formation physics as the Milky Way. We employ the semi-analytic model Galacticus to model dwarfs within a 600 kpc region. Utilizing astrophysical prescriptions and parameters matching the Milky Way satellites, we explore predictions for various properties and star formation histories (SFHs) to investigate environmental effects. We also reproduce cumulative luminosity and luminosity metallicity relations consistent with observations for the overall Centaurus A satellite population, while predicting half-light radii, velocity dispersion, and SFHs for the dwarf galaxies in Centaurus A. The agreement between our predicted SFHs for Centaurus A dwarfs and those of the Milky Way implies the presence of universal processes governing star formation in these galaxies. Overall, our findings shed light on the star formation physics of dwarf galaxies in the Centaurus A system, revealing insights into their properties and dependence on the host environment.

The search for extraterrestrial intelligence (SETI) is a field that has long been within the domain of traditional signal processing techniques. However, with the advent of powerful generative AI models, such as GPT-3, we are now able to explore new ways of analyzing SETI data and potentially uncover previously hidden signals. In this work, we present a novel approach for using generative AI to analyze SETI data, with focus on data processing and machine learning techniques. Our proposed method uses a combination of deep learning and generative models to analyze radio telescope data, with the goal of identifying potential signals from extraterrestrial civilizations. We also discuss the challenges and limitations of using generative AI in SETI, as well as potential future directions for this research. Our findings suggest that generative AI has the potential to significantly improve the efficiency and effectiveness of the search for extraterrestrial intelligence, and we encourage further exploration of this approach in the SETI community. (Disclosure: For the purpose of demonstration, the abstract and title were generated by ChatGPT and slightly modified by the lead author.

The existence of a stochastic gravitational wave background is indicated by the recent pulsar timing array (PTA) experiments. We study the enhanced production of second-order gravitational waves from the scalar perturbations when the universe experiences a transition from the early matter-dominated era to the radiation-dominated era due to Q-ball decay. We extend the analysis in previous work by including the frequency range where density perturbations go non-linear and find that the resultant gravitational wave spectrum can be consistent with that favored by the recent PTA experiment results.

Coronal mass ejections (CMEs) are explosive plasma phenomena prevalently occurring on the Sun and probably on other magnetically active stars. However, how their pre-eruptive configuration evolves toward the main explosion remains elusive. Here, based on comprehensive observations of a long-duration precursor in an event on 2012 March 13, we determine that the heating and slow rise of the pre-eruptive hot magnetic flux rope (MFR) are achieved through a precursor reconnection located above cusp-shaped high-temperature precursor loops. It is observed that the hot MFR threads are built up continually with their middle initially showing an "M" shape and then being separated from the cusp of precursor loops, causing the slow rise of the entire MFR. The slow rise in combination with thermal-dominated hard X-ray source concentrated at the top of the precursor loops shows that the precursor reconnection is much weaker than the flare reconnection of the main eruption. We also perform a three-dimensional magnetohydrodynamics simulation that reproduces the early evolution of the MFR transiting from the slow to fast rise. It is also disclosed that it is the magnetic tension force pertinent to "M"-shaped threads that drives the slow rise, which, however, evolves into a magnetic pressure gradient dominated regime responsible for the rapid-acceleration eruption.

We use an analytic framework to calculate the evolution of binary orbits under a physically-motivated model that accounts for angular momentum loss associated with winds from an accretion disk around the compact objected accretor. Our prescription considers wind mass ejection from the surface of an accretion disk, accounting for a radial mass-loss dependence across the disk surface. We compare this to the standard prescription of angular momentum loss associated with isotropic mass loss from the vicinity of the accretor. The angular momentum loss from a disk-wind is always larger. For mass ratios, $q$, between $2$--$10$, angular momentum loss via a disk wind is $\simeq3$--$40$ times greater than the standard prescription. For the majority of mass ratios and disk properties, accounting for the disk wind can result in considerably smaller orbital separations compared to the standard formalism; the differences being $\simeq 60\%$ depending on how long the effect is integrated for. We conclude that it is important to consider the effects of angular momentum loss from a disk wind when evolving binary orbits.

Based on the Seventh Data Release (DR7) quasar catalog from the Sloan Digital Sky Survey, we investigate the variability of optical quasars in W1, W2, W3 and W4 bands of the Wide-field Infrared Survey Explorer (WISE) and the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE). Adopting the structure function method, we calculate the structure function ($\rm\delta t$=1 yr) which shows no obvious correlations with the bolometric luminosity, the black hole mass and the Eddington ratio. The ensemble structure functions in W1 and W2 bands show that the SF slopes are steeper than those in previous studies which may be caused by different cadence and observational epoch number. We further investigate the relation of variability amplitude $\sigma_m$ between mid-infrared band and optical band, but no obvious correlation is found. No correlation is found between W1-W2 and g-r color. We think the mid-infrared emission of quasars may be smoothed out by the extended dust distribution, thus leading to no obvious correlation. For the radio-loud quasar sub-sample, we further analyze the relation between the variability amplitude in the mid-infrared band and the radio luminosity at 6 cm, but no obvious correlations are found, which indicate the mid-infrared emission contributed from the synchrotron radiation of the relativistic jet is very weak.

We have studied the evolution of HuBi 1-like planetary nebulae, considering several stages of mass injection. We have carried out numerical ionization+1D hydrodynamics+atomic/ionic rate models with our code Coral 1D to reproduce planetary nebulae that present multiple shells produced by different ejection events around the ionizing source. Furthermore, we are interested in comparing numerical simulations with H$\alpha$ and [NII]$\lambda$6584 emission structures and the position-velocity diagrams observed in HuBi 1. This object also has a phase where it has drastically decreased the injection of ionized photons ejected from the source. The result of these different stages of ejection is a nebula with intense [NII] line emission in the inner part of the planetary nebula and an extended HII recombination line emission around the central zone. The model for HuBi 1 shows the capability of our code to explain the hydrodynamical and photoionization evolution in ionization nebulae. This is our first step with a 1D code to study these two physical phenomena at the same time.

Hot exozodiacal dust (HEZD) found around main-sequence stars through interferometric observations in the photometric bands H to L is located close to the dust sublimation radius, potentially at orbital radii comparable to those of close-in exoplanets. Consequently, HEZD has a potential influence on the analysis of the scattered-light polarization of close-in exoplanets and vice versa. We analyze the impact of HEZD on the polarimetric characterization of close-in exoplanets. This study is motivated in particular by the recently proven feasibility of exoplanet polarimetry. Applying the 3D Monte Carlo radiative transfer code POLARIS in an extended and optimized version for radiative transfer in exoplanetary atmospheres and an analytical tool for modeling the HEZD, we simulated and compared the polarization characteristics of the wavelength-dependent scattered-light polarization of HEZD and close-in exoplanets. The varied parameters are the planetary phase angle ($0^\circ-180^\circ$), the dust grain radius ($0.02\ \mu$m $- \ 10\ \mu$m), the HEZD mass ($10^{-10}$$\rm{M}_{\oplus}$ $-\ 10^{-8}$$\rm{M}_{\oplus}$), the orbital inclination ($0^\circ-90^\circ$), the composition of the planetary atmosphere (Mie and Rayleigh scattering atmosphere), the orbital radius of the HEZD ($0.02$ au $-\ 0.4$ au), and the planetary orbital radius ($0.01$ au $-\ 0.05$ au). The dust grain radius has the strongest influence on the polarimetric analysis due to its significant impact on the wavelength-dependent polarization characteristics and the total order of magnitude of the scattered-light polarization. In certain scenarios, the scattered-light polarization of the HEZD even exceeds that of the close-in exoplanet.

We investigate the spectrum of gravitational waves arising from primordial inflation in the presence of a string-theoretical higher curvature correction, specifically, the Gauss-Bonnet coupling term for the inflaton (modulus) field. We show that if the modulus field exhibits a wall-crossing like behavior in the moduli space, there can be a period of Gauss-Bonnet coupling term domination during the usual slow-roll. This phenomenon is potentially detectable as the gravitational wave spectrum exhibits a characteristic peak caused by the brief domination of the Gauss-Bonnet coupling term. We explore the possibility of measuring such gravitational waves with pulsar timing array experiments such as NANOGrav, and future space-borne interferometers such as LISA, DECIGO, and Taiji.

A fundamental question of extra-galactic astronomy that is yet to be fully understood, concerns the evolution of the star formation rate (SFR) and supermassive black hole (SMBH) activity with cosmic time, as well as their interplay and how it impacts galaxy evolution. A primary focus that could shed more light on these questions is the study of merging systems, comprising highly star-forming galaxies (SFGs) and active galactic nuclei (AGN) at the earliest stages of galactic formation. However, it is essential to explore complementary selection methods across multiple wavelengths. The primary objective of this study is to conduct a comprehensive analysis of a sample of high-redshift ($z>3$) far-infrared (far-IR) and radio-emitting galaxies in the highest possible spatial resolution. In order to select the galactic population of our interest, we selected galaxies that present relatively compact radio morphologies at 1.4 GHz as well as a far-IR spectrum that peaks in flux at $\lambda \geq 350 \, \mu m$. For these selection criteria, we used the COSMOS and ECDF-S fields, which provide high spectral and spatial resolution at a multi-wavelength scale. We derived a sample of eight galaxies that were identified either photometrically or spectroscopically at $z>3$ from literature studies and by our team. A thorough investigation of available optical, near-IR, and millimetre (mm) imaging reveals a possible merging scenario in five out of eight cases in our sample. Additionally, available multi-wavelength photometry strongly suggests active star formation at the $10^3 \, M_{\odot}/yr$ level in massive systems co-hosting an active SMBH. Comparison of these results with previous studies, suggests that our selection method preferentially identifies galaxies hosting an active SMBH, as well as a strong SFG component, resulting in high SFR and IR luminosity.

Neutrino Events Reconstruction has always been crucial for IceCube Neutrino Observatory. In the Kaggle competition "IceCube -- Neutrinos in Deep Ice", many solutions use Transformer. We present ISeeCube, a pure Transformer model based on torchscale (the backbone of BEiT-3). Our model has the potential to reach state-of-the-art. By using torchscale, the lines of code drop sharply by about 80% and a lot of new methods can be tested by simply adjusting configs. Also, a relatively primitive loss function Mean Squared Error (MSE) can work really well. Since the model is simple enough, it also has the potential to be used for more purposes such as energy reconstruction, and many new methods such as combining it with GraphNeT can be tested more easily. The code and pretrained models are available at https://github.com/ChenLi2049/ISeeCube

We present measurements of stellar populations properties of a z = 9.1 gravitationally lensed galaxy MACS1149-JD1 using deep JWST NIRISS slitless spectroscopy as well as NIRISS and NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The galaxy is split into four components. Three magnified (${\mu}$ ~ 17) star-forming components are unresolved, giving intrinsic sizes of < 50pc. In addition, the underlying extended component contains the bulk of the stellar mass, formed the majority of its stars ~ 50Myr earlier than the other three components and is not the site of the most active star formation currently. The NIRISS and NIRCam resolved photometry does not confirm a strong Balmer break previously seen in Spitzer. The NIRISS grism spectrum has been extracted for the entire galaxy and shows a clear continuum and Lyman-break, with no Lyman-${\alpha}$ detected.

In this introductory chapter of the Special Issue entitled `The Structure and Evolution of Stars', we highlight the recent major progress made in our understanding in the physics that governs stellar interiors. In so doing, we combine insight from observations, 1D evolutionary modelling and 2+3D rotating (magneto)hydrodynamical simulations. Therefore, a complete and compelling picture of the necessary ingredients in state-of-the-art stellar structure theory and areas in which improvements still need to be made are contextualised. Additionally, the over-arching perspective that links all the themes of subsequent chapters is presented.

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background (CMB) over ~75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the large angular scale CMB polarization to constrain the tensor-to-scalar ratio and the optical depth to last scattering. This paper presents the optical characterization of the 90GHz telescope, which has been observing since July 2018. Observations of the Moon establish the pointing while dedicated observations of Jupiter are used for beam calibration. The standard deviations of the pointing error in azimuth, elevation, and boresight angle are 1.3, 2.1, and 2.0 arcminutes, respectively, over the first 3 years of observations. This corresponds to a pointing uncertainty ~7% of the beam's full width at half maximum (FWHM). The effective azimuthally-symmetrized 1D beam estimated at 90 GHz from per detector intensity beam maps has a FWHM of 0.614+/-0.003 deg and a solid angle of 136.3+/-0.6(stats.)+/-1.1(sys.) usr integrated to a radius of 4 deg. The corresponding beam window function drops to b_ell^2 = 0.92, 0.70, 0.14 at ell = 30, 100, 300, respectively, with relative uncertainties < 2% for ell < 200. Far-sidelobes are studied using detector-centered intensity maps of the Moon and measured to be at a level of 10^-3 or below relative to the peak. The polarization angle of Tau A estimated from preliminary survey maps is 149.6+/-0.2(stats.) deg in equatorial coordinates consistent with prior measurements. Instrumental temperature-to-polarization (T-to-P) leakage is measured at a 95% confidence upper limit of (1.7+/-0.1) x 10^-3 in single detector demodulated data using observations of Jupiter and the Moon. Using pair-differenced demodulated data, a 95% confidence upper limit of 3.6 x 10^-4 is obtained on the T-to-P leakage.

We followed-up with ESPRESSO the K0V star HIP 29442 (TOI-469), already known to host a validated sub-Neptune companion TOI-469.01. We aim to verify the planetary nature of TOI-469.01. We modelled radial velocity and photometric time series to measure the dynamical mass, radius, and ephemeris, and to characterise the internal structure and composition of TOI-469.01. We confirmed the planetary nature of TOI-469.01. Thanks to ESPRESSO we discovered two additional close-in companions. We also detected their low signal-to-noise transit signals in the TESS light curve. HIP 29442 is a compact multi-planet system, and the three planets have orbital periods $P_{\rm orb, b}=13.63083\pm0.00003$, $P_{\rm orb, c}=3.53796\pm0.00003$, and $P_{\rm orb, d}=6.42975^{+0.00009}_{-0.00010}$ days, and we measured their masses with high precision: $m_{\rm p, b}=9.6\pm0.8~M_{\oplus}$, $m_{\rm p, c}=4.5\pm0.3~M_{\oplus}$, and $m_{\rm p, d}=5.1\pm0.4~M_{\oplus}$. We measured radii and bulk densities of all the planets (the 3$\sigma$ confidence intervals are shown in parenthesis): $R_{\rm p, b}=3.48^{+0.07 (+0.19)}_{-0.08 (-0.28)} ~R_{\oplus}$ and $\rho_{\rm p, b}=1.3\pm0.2 (0.3) g~cm^{-3}$; $R_{\rm p, c}=1.58^{+0.10 (+0.30)}_{-0.11 (-0.34)}~R_{\oplus}$ and $\rho_{\rm p, c}=6.3^{+1.7 (+6.0)}_{-1.3 (-2.7)} g~cm^{-3}$; $R_{\rm p, d}=1.37\pm0.11^{(+0.32)}_{(-0.43)}~R_{\oplus}$ and $\rho_{\rm p, d}=11.0^{+3.4 (+21.0)}_{-2.4 (-6.3)} g~cm^{-3}$. We used the more conservative 3$\sigma$ confidence intervals for the radii as input to the interior structure modelling. We find that HIP 29442 $b$ appears as a typical sub-Neptune, likely surrounded by a gas layer of pure H-He with a mass of $0.27^{+0.24}_{-0.17} M_{\oplus}$ and a thickness of $1.4\pm0.5 R_{\oplus}$. For the innermost companions HIP 29442 $c$ HIP 29442 $d$, the model supports an Earth-like composition.

Superthin galaxies are bulgeless low surface brightness galaxies with unusually high major-to-minor axes ratio of the stellar disc, i.e.,$10<a/b<20$. We present Giant Metrewave Radio Telescope (GMRT) \HI{} 21cm radio-synthesis observations of FGC 2366, the thinnest galaxy known with $a/b=21.6$. Employing the 3-D tilted-ring modelling using Fully Automated TiRiFiC (FAT), we determine the structure and kinematics of the \HI{} gas disc, obtaining an asymptotic rotational velocity equal to 100 \kms and a total \HI{} mass equal to 10$^9 M_{\odot}$. Using $z$-band stellar photometry, we obtain a central surface brightness of 22.8 mag ${\rm{arcsec}}^{-2}$, a disc scale length of 2.6 kpc, and a scaleheight of 260 pc. Next, we determine the dark matter density profile by constructing a mass model and find that an NFW dark matter halo best fits the steeply-rising rotation curve. With the above mass inventory in place, we finally construct the dynamical model of the stellar disc of FGC 2366 using the stellar dynamical code "AGAMA". To identify the key physical mechanisms responsible for the superthin vertical structure, we carry out a Principal Component Analysis of the data corresponding to all the relevant dynamical parameters and $a/b$ for a sample of superthin and extremely thin galaxies studied so far. We note that the first two principal components explain 80$\%$ of the variation in the data, and the significant contribution is from the compactness of the mass distribution, which is fundamentally responsible for the existence of superthin stellar discs.

The James Webb Space Telescope (JWST) has revealed extremely distant galaxies at unprecedentedly early cosmic epochs from its deep imaging using the technique of photometric redshift estimation, with its subsequent spectroscopy confirming their redshifts unambiguously, demonstrating the ability of JWST to probe the earliest galaxies, one of its major scientific goals. However, as larger samples continue to be followed up spectroscopically, it has become apparent that nearly all photometric redshifts at these epochs are biased high with confidence >>99%, for as yet unclear reasons. Here we show that this is the same statistical effect that was predicted in different contexts by Sir Arthur Eddington in 1913, in that there exist more lower redshift galaxies to be scattered upwards than the reverse. The bias depends on the shape of the intrinsic redshift distribution, but as an approximate heuristic, all ultra-high photometric redshift estimates must be corrected downwards by up to one standard deviation.

Context. The metal mass fraction of the Sun Z is a key constraint in solar modelling, but its value is still under debate. The standard solar chemical composition of the late 2000s have the ratio of metals to hydrogen Z/X = 0.0181, with a small increase to 0.0187 in 2021, as inferred from 3D non-LTE spectroscopy. However, more recent work on a horizontally and temporally averaged <3D> model claim Z/X = 0.0225, consistent with the high values of twenty-five years ago based on 1D LTE spectroscopy. Aims. We aim to determine a precise and robust value of the solar metal mass fraction from helioseismic inversions, thus providing independent constraints from spectroscopic methods. Methods. We devise a detailed seismic reconstruction technique of the solar envelope, combining multiple inversions and equations of state to accurately and precisely determine the metal mass fraction value. Results. We show that a low value of the solar metal mass fraction corresponding to Z/X = 0.0187 is favoured by helioseismic constraints and that a higher metal mass fraction corresponding to Z/X = 0.0225 are strongly rejected by helioseismic data. Conclusions. We conclude that direct measurement of the metal mass fraction in the solar envelope favours a low metallicity, in line with the 3D non-LTE spectroscopic determination of 2021. A high metal mass fraction as measured using a <3D> model in 2022 is disfavoured by helioseismology for all modern equations of state used to model the solar convective envelope.

In very-high-energy (VHE) gamma-ray astronomy, the community is converging towards the use of a common open data format, called "Data formats for Gamma-ray Astronomy", for the high-level data products. This format is in use for ground-based TeV observatories like H.E.S.S., MAGIC or HAWC, some of whom plan to openly release high-level data products. These efforts are parallel to the development and use of open analysis software such as the Gammapy package. This open initiative has shown that it is possible to define common standards even without governance. With the advent of open VHE observatories (e.g. CTAO, KM3NeT) and an increase in both multi-wavelength and multi-messenger studies, such standards should evolve to support all of VHE multi-messenger astrophysics. For these reasons, a new initiative has been created to specify formats of high-level data from very and ultra high energy gamma-ray facilities and from VHE neutrino detectors. It also aims to better respect the FAIR principles and the IVOA recommendations.This communication will present the Very-high-energy Open Data Format (VODF) project that has been established by eleven VHE astroparticle facilities. Its structure, its organisation and its goal will be presented. Anchored in Open Science, our goal is to solicit comments and future contributions from the VHE astrophysics community.

Since its start in 2014, the lightweight open source Python library Gammapy has come a long way to become a popular data analysis package for high-energy astrophysics. Selected as the official CTAO Science Analysis tool, it is also an approved analysis software within the H.E.S.S. and MAGIC collaborations. The first long-term version, Gammapy v1.0 was released on late 2022. It is compliant with several well-established data conventions in high-energy astrophysics, and provides serialised data products that are interoperable with other software. Event lists and instrument response functions curated within the same format from various instruments can be reduced to data binned in energy, time or spatial coordinates. Thereafter, the flux and morphology of one or more gamma-ray sources can be estimated using Poisson maximum likelihood fitting and assuming a variety of spectral, temporal and spatial models. Flux points, likelihood profiles and light curves extractions are supported. Complex user defined likelihoods and models can also be implemented. In this contribution, we will highlight the main features of Gammapy v1.0, including data reduction and analysis examples from different space and ground-based instruments, applications of various background rejection techniques, and a simultaneous fitting across multiple instruments with astrophysical models. We will also present our plans for the future, showcasing new features such as the support of different event types, unbinned likelihood analysis, spectral unfolding and transient source detections. In addition to an improved API with distributed computing for scalable analysis, enhanced support for all-sky instruments like Fermi-LAT and HAWC is foreseen.

The maximum energy of electrons accelerated by supernova remnants (SNR) is typically limited by radiative losses. In this scenario, the synchrotron cooling time scale is equal to the acceleration time scale. On the other hand, the low propagation speed of a shock in a dense medium is expected to result in an extended acceleration time scale, thus inducing a decrease in the maximum electron energy for a given SNR age and in the X-ray nonthermal flux. The young Kepler's SNR shows an enhanced efficiency of the acceleration process, which is close to the Bohm limit in the north of its shell, where the shock is slowed down by a dense circumstellar medium. Conversely, in the south, where no interaction with a dense medium is evident and the shock speed is high, the acceleration proceeds with a higher Bohm factor. To investigate this scenario, we studied the temporal evolution of the non-thermal emission, taking advantage of two Chandra X-ray observations of Kepler's SNR (performed in 2006 and 2014). We analyzed the spectra of different filaments both in the north and south of the shell, and measured their proper motion. We found a region with low shock velocity where we measured a significant decrease in flux from 2006 to 2014. This could be the first evidence of fading synchrotron emission in Kepler's SNR. This result suggests that under a certain threshold of shock speed the acceleration process could exit the loss-limited regime.

Outflowing wind as one type of AGN feedback, which involves noncollimated ionized winds prevalent in Seyfert-1 AGNs, impacts their host galaxy by carrying kinetic energy outwards. However, the distance of the outflowing wind is poorly constrained due to a lack of direct imaging observations, which limits our understanding of their kinetic power and therefore makes the impact on the local galactic environment unclear. One potential approach involves a determination of the density of the ionized plasma, making it possible to derive the distance using the ionization parameter {\xi}, which can be measured based on the ionization state. Here, by applying a new time-dependent photoionization model, tpho, in SPEX, we define a new approach, the tpho-delay method, to calculate/predict a detectable density range for warm absorbers of NGC 3783. The tpho model solves self-consistently the time-dependent ionic concentrations, which enables us to study delayed states of the plasma in detail. We show that it is crucial to model the non-equilibrium effects accurately for the delayed phase, where the non-equilibrium and equilibrium models diverge significantly. Finally, we calculate the crossing time to consider the effect of the transverse motion of the outflow on the intrinsic luminosity variation. It is expected that future spectroscopic observations with more sensitive instruments will provide more accurate constraints on the outflow density, and thereby on the feedback energetics.

Self-interacting dark matter (SIDM) is an alternative to the standard collisionless cold dark matter model (CDM), allowing for interactions between the dark matter particles through the introduction of a self-scattering cross-section. However, the observable effects between these two scenarios are hard to detect. In this work we present a detailed analysis of an application of galaxy-galaxy lensing to measure with high precision the shapes of cluster halos and how this approach can be used to obtain information regarding the nature of the dark matter particle. Using two sets of simulated data, SIDM and CDM simulations, we compute stacked shear maps centred on several subsets of halos with masses $\gtrsim 10^{13.5} M_\odot$. From these maps, we obtain the quadrupole profiles related to the mean elongation of the particle distribution from which the shape parameters are derived. Accounting for a radial shape variation, this technique provides an enhancement of the observed differences between the simulated data-sets which successfully differentiate between the two cosmological scenarios. In particular, we obtain a larger slope of the shape-radial relation for the halos identified in the SIDM simulation, which are rounder towards the centre. Also, as approaching to the mean virial radius, the projected semi-axis ratios converge to similar values than in the CDM simulation. Moreover, we account for the impact of the neighbouring mass, where a more strongly elongated distributions is found for the halos in the SIDM simulation, indicating that under dark matter self interaction, the large scale structure imprints a more coherent accretion process.

The Extreme Universe Space Observatory on a Super Pressure Balloon II (EUSO-SPB2) launched from Wanaka, New Zealand on May 13th 2023. Consisting of two optical telescopes, EUSO-SPB2 aimed to search for very high energy neutrinos (E>PeV) via Cherenkov radiation, and ultra high energy cosmic rays (UHECRs, E>EeV) via ultraviolet fluorescence. Building on the EUSO-Balloon (2014) and EUSO-SPB1 (2017) missions, the Fluorescence Telescope (FT) comprises 108 multi-anode photomultiplier tubes at the focus of a one meter entrance diameter Schmidt telescope. The FT pointed down at the atmosphere below the SPB's altitude of 33 km. The mission duration was planned to reach up to 100 days. Prior to flight, the instrument was extensively tested in the laboratory and in the field. These measurements, combined with simulations led to an expected peak energy sensitivity around 3 EeV. Combined with a three-times-larger field-of-view than previous EUSO balloon missions, this resulted in an expected observation rate of one UHECR shower per ten hours of observation. The FT was expected to perform the first measurement of UHECRs via fluorescence from sub-orbital space, but was unable to, due to a shortened flight. Nonetheless, these observations of EUSO-SPB2 serve as a stepping stone to future satellite-based missions, such as the Probe of Extreme MultiMessenger Astrophysics (POEMMA), with enormous exposure to the highest energy cosmic rays with all sky coverage. In this contribution we will discuss the performance of the FT in flight as well as preliminary results.

The Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2) is the most advanced balloon mission undertaken by the JEM-EUSO collaboration. EUSO-SPB2 is built on the experience of previous stratosphere missions, EUSO-Balloon and EUSO-SPB, and of the Mini-EUSO space mission currently active onboard the International Space Station. EUSO- SPB2 is equipped with two instruments: a fluorescence telescope aimed at registering ultra-high energy cosmic rays (UHECRs) with an energy above 2 EeV and a Cherenkov telescope built to measure direct Cherenkov emission from cosmic rays with energies above 1 PeV. The EUSO-SPB2 mission will provide pioneering observations on the path towards a space-based multi-messenger observatory. As such, a special attention was paid to the development of triggers and other software aimed at comprehensive data analysis. A whole number of methods based on machine learning (ML) and neural networks was developed during the construction of the experiment and a few others are under active development. Here we provide a brief review of the ML-based methods already implemented in the instrument and the ground software and report preliminary results on the ML-based reconstruction of UHECR parameters for the fluorescence telescope.

We present the first results from CECILIA, a Cycle 1 JWST NIRSpec/MSA program that uses ultra-deep ~30 hour G235M/F170LP observations to target multiple electron temperature-sensitive auroral lines in the spectra of 33 galaxies at z~1-3. Using a subset of 23 galaxies, we construct two ~600 object-hour composite spectra, both with and without the stellar continuum, and use these to investigate the characteristic rest-optical (5700-8500 Angstrom) spectrum of star-forming galaxies at the peak epoch of cosmic star formation. Emission lines of eight different elements (H, He, N, O, Si, S, Ar, and Ni) are detected, with most of these features observed to be <3% the strength of H-alpha. We report the characteristic strength of three auroral lines ([NII]5756, [SIII]6313, and [OII]7322,7332), as well as other semi-strong and faint emission lines, including forbidden [NiII]7380,7414 and the OI 8449 recombination line, some of which have never before been observed outside of the local universe. Using these measurements, we find T_e[NII]=13630+/-2540$ K, representing the first measurement of electron temperature using [NII] in the high-redshift universe. We also see evidence for broad line emission with a FWHM of ~544 km/s; the broad component of H-alpha is 6.01-28.31% the strength of the narrow component and likely arises from star-formation driven outflows. Finally, we briefly comment on the feasibility of obtaining large samples of faint emission lines using JWST in the future.

M35 is a young open cluster and home to an extensive binary population. Using Gaia DR3, Pan-STARRS, and 2MASS photometry with the Bayesian statistical software, BASE-9, we derive precise cluster parameters, identify single and binary cluster members, and extract their masses. We identify 571 binaries down to Gaia G = 20.3 and a lower-limit on the binary frequency of f_b = 0.41 +/- 0.02. We extend the binary demographics by many magnitudes faint-ward of previous (radial-velocity) studies of this cluster and further away from the cluster center (1.78-degrees, roughly 10 core radii). We find the binary stars to be more centrally concentrated than the single stars in the cluster. Furthermore, we find strong evidence for mass segregation within the binary population itself, with progressively more massive binary samples becoming more and more centrally concentrated. For the single stars, we find weaker evidence for mass segregation; only the most massive single stars (> 2.5MSun) appear more centrally concentrated. Given the cluster age of ~200 Myr, and our derived half-mass relaxation time for the cluster of 230 +/- 84 Myr, we estimate ~47% of the binary stars and ~12% of single stars in the cluster have had time to become dynamically mass segregated. Importantly, when we investigate only stars with mass segregation timescales greater than the cluster age, we still find the binaries to be more centrally concentrated than the singles, suggesting the binaries may have formed with a primordially different spatial distribution than the single stars.

We revisit a question asked by Dyson: "Is a graviton detectable?" We demonstrate that in both Dyson's original sense and in a more modern measurement-theoretic sense, it is possible to construct a detector sensitive to single gravitons, and in fact a variety of existing and near-term gravitational wave detectors can achieve this. However, while such a signal would be consistent with the quantization of the gravitational field, we draw on results from quantum optics to show how the same signal could just as well be explained via classical gravitational waves. We outline the kind of measurements that would be needed to demonstrate quantization of gravitational radiation and explain why these are substantially more difficult than simply counting graviton clicks or observing gravitational noise in an interferometer, and likely impossible to perform in practice.

Phantom scalar theories are widely considered in cosmology, but rarely at the quantum level, where they give rise to negative-energy ghost particles. These cause decay of the vacuum into gravitons and photons, violating observational gamma-ray limits unless the ghosts are effective degrees of freedom with a cutoff $\Lambda$ at the few-MeV scale. We update the constraints on this scale, finding that $\Lambda \lesssim 19$ MeV. We further explore the possible coupling of ghosts to a light, possibly massless, hidden sector particle, such as a sterile neutrino. Vacuum decays can then cause the dark matter density of the universe to grow at late times. The combined phantom plus dark matter fluid has an effective equation of state $w < -1$, and functions as a new source of dark energy. We derive constraints from cosmological observables on the rate of vacuum decay into such a phantom fluid. We find a mild preference for the ghost model over the standard cosmological one, and a modest amelioration of the Hubble and $S_8$ tensions.

We propose a Bayesian meta-analysis to infer the current expansion rate of the Universe, called the Hubble constant ($H_0$), via time delay cosmography. Inputs of the meta-analysis are estimates of two properties for each pair of gravitationally lensed images; time delay and Fermat potential difference estimates with their standard errors. A meta-analysis can be appealing in practice because obtaining each estimate from even a single lens system involves substantial human efforts, and thus estimates are often separately obtained and published. This work focuses on combining these estimates from independent studies to infer $H_0$ in a robust manner. For this purpose, we adopt Student's $t$ error for the inputs of the meta-analysis. We investigate properties of the resulting $H_0$ estimate via two simulation studies with realistic imaging data. It turns out that the meta-analysis can infer $H_0$ with sub-percent bias and about 1 percent level of coefficient of variation, even when 30 percent of inputs are manipulated to be outliers. We also apply the meta-analysis to three gravitationally lensed systems, and estimate $H_0$ by $75.632 \pm 6.918$ (km/second/Mpc), which covers a wide range of $H_0$ estimates obtained under different physical processes. An R package, h0, is publicly available for fitting the proposed meta-analysis.

While first order phase transitions (FOPTs) have been extensively studied as promising cosmological sources of gravitational waves, the phenomenon of particle production from the dynamics of the background field during FOPTs has received relatively little attention in the literature, where it has only been studied with semi-analytic estimates in some simplified settings. This paper provides improved numerical studies of this effect in more realistic frameworks, revealing important qualitative details that have been missed in the literature. We also provide easy to use analytic formulae that can be used to calculate particle production in generic FOPT setups.

Joint observations of gravitational waves and electromagnetic counterparts will answer questions about cosmology, gamma-ray bursts, and the behaviour of matter at supranuclear densities. The addition of a Southern-Hemisphere gravitational-wave observatory to proposed global networks creates a longer baseline, which is beneficial for sky localisation. We analyse how an observatory in Australia can enhance the multi-messenger astronomy capabilities of future networks. We estimate the number of binary neutron star mergers with joint observations of gravitational waves and kilonova counterparts detectable by the Vera C. Rubin Observatory. First, we consider a network of upgrades to current observatories. Adding an Australian observatory to a three-observatory network (comprising two observatories in the USA and one in Europe) boosts the rate of joint observations from $2.5^{+4.5}_{-2.0}$ per year to $5.6^{+10}_{-4.5}$ per year (a factor of two improvement). Then, we consider a network of next-generation observatories. Adding a $20$ km Australian observatory to a global network of a Cosmic Explorer $40$ km in the USA and an Einstein Telescope in Europe only marginally increases the rate from $40^{+71}_{-32}$ per year to $44^{+79}_{-35}$ per year (a factor of 1.1 improvement). The addition of an Australian observatory, however, ensures that at least two observatories are online far more often. When the Cosmic Explorer $40$ km is offline for a major upgrade, the Australian observatory increases the joint observation rate from $0.5^{+0.8}_{-0.4}$ per year to $38^{+68}_{-30}$ per year (a factor of 82 improvement). When the Einstein Telescope is offline, the joint observation rate increases from $0.2^{+0.3}_{-0.1}$ per year to $19^{+34}_{-15}$ per year (a factor of 113 improvement). We sketch out the broader science case for a Southern-Hemisphere gravitational-wave observatory.

We discuss halo-independent constraints on the Inelastic Dark Matter scenario, in which a Weakly Interaction Massive Particle (WIMP) state $\chi$ with mass $m_\chi$ interacts with nuclear targets by upscattering to a heavier state $\chi^{\prime}$ with mass $m_\chi+\delta$. In order to do so we adopt the single-stream method, that exploits the complementarity of Direct Detection (DD) and Capture in the Sun to extend the experimental sensitivity to the full range of incoming WIMP speeds. We show that a non-vanishing mass splitting $\delta$ modifies such range, and that for particular combinations of $m_\chi$ and $\delta$ the complementarity between the two detection techniques required by the method is lost. Specifically, assuming for the escape velocity in our Galaxy $u_{esc}$ the reference value $u_{esc}^{ref}$ = 560 km/s a halo-independent bound is possible when $\delta\lesssim$ 510 keV for a Spin-Independent interaction and when $\delta\lesssim$ 245 keV for a Spin-Dependent interaction (with the Spin-Independent value slightly reduced to $\delta\lesssim$ 490 keV when $u_{esc}>u_{esc}^{ref}$). We find that in the low-mass regime the bound from capture in the Sun is always more constraining than that for DD and is sufficient alone to provide a halo-independent constraint, while for large WIMP masses the halo-independent bound is given by a combination of capture in the Sun and DD. In this latter case, while for increasing values of $\delta$ the sensitivity range of initial speeds of the WIMP is reduced for both DD and capture in the Sun, such effect is more pronounced for DD than for capture. We also find that, for $u_{esc}$ = $u_{esc}^{ref}$, unless the mass of the target used in DD is larger than about four times that of the target driving capture in the Sun, DD does not play any role in the determination of the maximal value of $\delta$ for which a halo-independent bound is possible.

A series of pulsar timing arrays (PTAs) recently observed gravitational waves at the nanohertz frequencies. Motivated by this remarkable result, we present a novel class of Pati-Salam models that give rise to a network of metastable cosmic strings, offering a plausible explanation for the observed PTA data. Besides, we introduce a hybrid inflationary scenario to eliminate magnetic monopoles that arise during the spontaneous symmetry breaking of the Pati-Salam gauge group to the Standard Model. The resulting scalar spectral index is compatible with Planck data, and the tensor-to-scalar ratio is anticipated to be extremely small. Moreover, we incorporate a non-thermal leptogenesis to generate the required baryon asymmetry in our framework. Finally, the gravitational wave spectra generated by the metastable cosmic strings not only correspond to signals observed in recent PTAs, including NANOGrav, but are also within the exploration capacity of both present and future ground-based and space-based experiments.

This work presents a large-scale simulation study investigating the deployment and operation of distributed swarms of CubeSats for interplanetary missions to small celestial bodies. Utilizing Taylor numerical integration and advanced collision detection techniques, we explore the potential of large CubeSat swarms in capturing gravity signals and reconstructing the internal mass distribution of a small celestial body while minimizing risks and Delta V budget. Our results offer insight into the applicability of this approach for future deep space exploration missions.

The possibility that the vacuum energy density (VED), $\rho_{\rm vac}$, could be time dependent in the expanding Universe is intuitively more reasonable than just a rigid cosmological constant for the entire cosmic history. The dynamics of $\rho_{\rm vac}=\rho_{\rm vac}(H)$ as a function of the Hubble rate, $H(t)$, most likely contributes to alleviate cosmological problems and tensions, having also implications on the so-called fundamental `constants' of Nature, which should be slowly drifting with the cosmic expansion owing to the fluctuations of the quantum vacuum. This includes the gravitational `constant' $G$, but also the gauge and Yukawa couplings as well as the particle masses themselves (both of dark matter and baryonic matter). The subtle exchange of energy involved is the basis for the ``micro and macro connection''. Herein, I discuss not only this connection as a possibility but show that it is in fact a generic prediction of QFT in cosmological spacetime which is fully compatible with general covariance. This fact has not been pointed out until recently when an appropriate renormalization framework for the VED has been found which is free from the usual conundrums associated with the cosmological constant problem.