Papers for Tuesday, May 16 2023

The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use SPH simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L$_{*}$ galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume, Romulus25, with stellar masses between 3 $\times$ 10$^{9}$ - 3 $\times$ 10$^{11}$ M$_{\odot}$. We measure the fraction of metals remaining in the ISM and CGM of each galaxy, and calculate the expected mass of its SMBH based on the $M-\sigma$ relation. The deviation of each SMBH from its expected mass, $\Delta M_{BH}$ is compared to the potential of its host via $\sigma$. We find that SMBHs with accreted mass above the empirical $M-\sigma$ relation are about 15\% more effective at removing metals from the ISM than under-massive SMBHs in star forming galaxies. Over-massive SMBHs suppress the overall star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little evidence for the evacuation of gas from their halos, in contrast with other simulations. Finally, we predict that C IV column densities in the CGM of L$_{*}$ galaxies may depend on host galaxy SMBH mass. Our results show that the scatter in the low mass end of $M-\sigma$ relation may indicate how effective a SMBH is at the local redistribution of mass in its host galaxy.

Long-term high-cadence measurements of stellar spectral variability are fundamental to better understand stellar atmospheric properties and stellar magnetism. These, in turn, are fundamental for the detectability of exoplanets as well as the characterization of their atmospheres and habitability. The Sun, viewed as a star via disk-integrated observations, offers a means of exploring such measurements while also offering the spatially resolved observations that are necessary to discern the causes of observed spectral variations. High-spectral resolution observations of the solar spectrum are fundamental for a variety of Earth-system studies, including climate influences, renewable energies, and biology. The Integrated Sunlight Spectrometer at SOLIS, has been acquiring daily high-spectral resolution Sun-as-a-star measurements since 2006.More recently, a few ground-based telescopes with the capability of monitoring the solar visible spectrum at high spectral resolution have been deployed (e.g. PEPSI, HARPS, NEID). However, the main scientific goal of these instruments is to detect exo-planets, and solar observations are acquired mainly as a reference. Consequently, their technical requirements are not ideal to monitor solar variations with high photometric stability, especially over solar-cycle temporal scales.The goal of this white paper is to emphasize the scientific return and explore the technical requirements of a network of ground-based spectrographs devoted to long-term monitoring of disk-integrated solar-spectral variability with high spectral resolution and high photometric stability, in conjunction with disk-resolved observations in selected spectral lines,to complement planet-hunter measurements and stellar-variability studies. The proposed network of instruments offers the opportunity for a larger variety of multidisciplinary studies.

The gravitational wave (GW) and neutrino signals from core-collapse supernovae (CCSNe) are expected to carry pronounced imprints of the standing accretion shock instability (SASI). We investigate whether the correlation between the SASI signatures in the GW and neutrino signals could be exploited to enhance the detection efficiency of GWs. We rely on a benchmark full-scale three-dimensional CCSN simulation with zero-age main sequence mass of $27\ M_\odot$. Two search strategies are explored: 1.~the inference of the SASI frequency range and/or time window from the neutrino event rate detectable at the IceCube Neutrino Observatory; 2.~the use of the neutrino event rate to build a matched filter template. We find that incorporating information from the SASI modulations of the IceCube neutrino event rate can increase the detection efficiency compared to standard GW excess energy searches up to $30\%$ for nearby CCSNe. However, we do not find significant improvements in the overall GW detection efficiency for CCSNe more distant than $1.5$~kpc. We demonstrate that the matched filter approach performs better than the unmodeled search method, which relies on a frequency bandpass inferred from the neutrino signal. The improved detection efficiency provided by our matched filter method calls for additional work to outline the best strategy for the first GW detection from CCSNe.

Ultracompact X-ray binaries (UCXBs) are a distinctive but elusive family of low-mass X-ray binaries (LMXBs) characterised by their tight orbits and degenerate donor stars. Here we present UltraCompCAT, the first online and comprehensive catalogue of UCXBs. The initial version of UltraCompCAT comprises 49 sources, including 20 'confirmed' UCXBs (those with a measured orbital period shorter than 80 min) and 25 systems that we label as 'candidate' based on their multi-wavelength phenomenology. For completeness, we also include four LMXBs with orbital periods in the range of 80 to 120 min, since they might be related (e.g. close progenitors) or even part of the UCXB population that evolved towards longer periods. We discuss the orbital period and Galactic distribution of the catalogue's sample. We provide evidence for the presence of at least two separate groups of UCXBs. One formed by persistent systems with orbital periods shorter than 30 min and a second group of transient objects (70 per cent) with periods in the range of 40 to 60 min. We show that the former group is dominated by sources formed in globular clusters, while the latter accounts for the (known) UCXB population in the Galactic field. We discuss the possible evolutionary channels for both groups.

Identifying the accelerators of Galactic cosmic ray protons (CRs) with energies up to a few PeV ($10^{15}$ eV) remains a theoretical and observational challenge. Supernova remnants (SNRs) represent strong candidates, as they provide sufficient energetics to reproduce the CR flux observed at Earth. However, it remains unclear whether they can accelerate particles to PeV energies, particularly after the very early stages of their evolution. This uncertainty has prompted searches for other source classes and necessitates comprehensive theoretical modeling of the maximum proton energy, $E_{\rm max}$, accelerated by an arbitrary shock. While analytic estimates of $E_{\rm max}$ have been put forward in the literature, they do not fully account for the complex interplay between particle acceleration, magnetic field amplification, and shock evolution. This paper uses a multi-zone, semi-analytic model of particle acceleration based on kinetic simulations to place constraints on $E_{\rm max}$ for a wide range of astrophysical shocks. In particular, we develop relationships between $E_{\rm max}$, shock velocity, size, and ambient medium. We find that SNRs can only accelerate PeV particles under a select set of circumstances, namely, if the shock velocity exceeds $\sim 10^4$ km s$^{-1}$ and escaping particles drive magnetic field amplification. However, older, slower SNRs may still produce observational signatures of PeV particles due to populations accelerated when the shock was younger. Our results serve as a reference for modelers seeking to quickly produce a self-consistent estimate of the maximum energy accelerated by an arbitrary astrophysical shock.

We present the analysis of the rest frame ultraviolet and optical spectra of 30 bright blue quasars at $z\sim3$, selected to examine the suitability of AGN as cosmological probes. In our previous works, we found an unexpectedly high fraction ($\approx 25 \%$) of X-ray weak quasars in the sample. The latter sources also display a flatter UV continuum and a broader and fainter CIV profile in the archival UV data with respect to their X-ray normal counterparts. Here we present new observations with the LBT in both the $zJ$ (rest-frame $\simeq$2300-3100 $\rm \mathring{A}$) and the $K_S$ ($\simeq$4750-5350 $\rm \mathring{A}$) bands. We estimated black hole masses ($M_{\rm BH}$) and Eddington ratios ($\lambda_{\rm Edd}$) from the from the H$\beta$ and MgII emission lines, finding that our $z\sim3$ quasars are on average highly accreting ($\langle \lambda_{\rm Edd} \rangle\simeq 1.2$ and $\langle M_{\rm BH} \rangle\simeq 10^{9.7}M_\odot$), with no difference in $\lambda_{\rm Edd}$ or $M_{\rm BH}$ between X-ray weak and X-ray normal quasars. From the $zJ$ spectra, we derive flux and equivalent width of MgII and FeII, finding that X-ray weak quasars display higher FeII/MgII ratios with respect to typical quasars. FeII/MgII ratios of X-ray normal quasars are instead consistent with other estimates up to $z\simeq6.5$, corroborating the idea of already chemically mature BLRs at early cosmic time. From the $K_S$ spectra, we find that all the X-ray weak quasars present generally weaker [OIII] emission (EW<10 $\rm \mathring{A}$) than the normal ones. The sample as a whole, however, abides by the known X-ray/[OIII] luminosity correlation, hence the different [OIII] properties are likely due to an intrinsically weaker [OIII] emission in X-ray weak objects, associated to the shape of the spectral energy distribution. We interpret these results in the framework of accretion-disc winds.

The first stars were likely more massive than those forming today and thus rapidly evolved, exploding as supernovae and enriching the surrounding gas with their chemical products. In the Local Group, the chemical signature of the first stars has been identified in the so-called Carbon-Enhanced Metal-Poor stars (CEMP-no). On the contrary, a similar C-excess was not found in dense neutral gas traced by high-redshift absorption systems. Here we discuss the recent discovery of three C-enhanced very metal-poor ([Fe/H]< -2) optically thick absorbers at redshift z ~ 3-4, reported by (Saccardi et al. 2023). We show that these absorbers are extra-galactic tracers of the chemical signatures of the first stars, analogous to the CEMP-no stars observed in the Galactic halo and ultra-faint dwarf galaxies. Furthermore, by comparing observations with model predictions we demonstrate that these systems have most likely been imprinted by first stars exploding as low-energy supernovae, which provided > 50% of the metals in these absorbers

We present photometric and spectroscopic data for the nearby Type I supernova (SN Ia) 2019eix (originally classified as a SN Ic), from its discovery day up to 100 days after maximum brightness. Before maximum light SN 2019eix resembles a typical SN Ic, albeit lacking the usual \ion{O}{1} feature. Its lightcurve is similar to the typical SN Ic with decline rates of ($\Delta M_{15,V}= 0.84$) and absolute magnitude of $M_{V}= -18.35$. However, after maximum light this SN has unusual spectroscopic features, a large degree of line blending, significant line blanketing in the blue ($\lambda < 5000$\AA), and strong Ca II absorption features during and after peak brightness. These unusual spectral features are similar to models of sub-luminous thermonuclear explosions, specifically double-detonation models of SNe Ia. Photometrically SN 2019eix appears to be somewhat brighter with slower decline rates than other double detonation candidates. We modeled the spectra using the radiative transfer code TARDIS using SN 1994I (a SN Ic) as a base model to see whether we could reproduce the unusual features of SN 2019eix and found them to be consistent with the exception of the \ion{O}{1} feature. We also compared SN 2019eix with double detonation models and found them to match the observations of SN 2019eix best, but failed to reproduce its full photometric and spectroscopic evolution.

Understanding differences between sub-stellar spectral data and models has proven to be a major challenge, especially for self-consistent model grids that are necessary for a thorough investigation of brown dwarf atmospheres. Using the supervised machine learning method of the random forest, we study the information content of 14 previously published model grids of brown dwarfs (from 1997 to 2021). The random forest method allows us to analyze the predictive power of these model grids, as well as interpret data within the framework of Approximate Bayesian Computation (ABC). Our curated dataset includes 3 benchmark brown dwarfs (Gl 570D, {\epsilon} Indi Ba and Bb) as well as a sample of 19 L and T dwarfs; this sample was previously analyzed in Lueber et al. (2022) using traditional Bayesian methods (nested sampling). We find that the effective temperature of a brown dwarf can be robustly predicted independent of the model grid chosen for the interpretation. However, inference of the surface gravity is model-dependent. Specifically, the BT-Settl, Sonora Bobcat and Sonora Cholla model grids tend to predict logg ~3-4 (cgs units) even after data blueward of 1.2 {\mu}m have been disregarded to mitigate for our incomplete knowledge of the shapes of alkali lines. Two major, longstanding challenges associated with understanding the influence of clouds in brown dwarf atmospheres remain: our inability to model them from first principles and also to robustly validate these models.

Using the GALEX archive, we have discovered extended structures around ten asymptotic giant branch (AGB) stars (out of a total 92 searched) emitting in the far-ultraviolet (FUV) band. In all but one, we find the typical morphology expected for a spherical wind moving relative to, and interacting with the ISM to produce an astrosphere. The exception is V\,Hya whose mass-ejection is known to be highly aspherical, where we find evidence of its large parabolic outflows interacting with the ISM, and its collimated, extreme velocity outflows interacting with the circumstellar medium. For 8 objects with relatively large proper motions, we find (as expected) that the termination-shock region lies in a hemisphere that contains the proper motion vector. Radial intensity cuts for each source have been used to locate the termination shock and the astropause's outer edge. In a few objects, the cuts also reveal faint emission just outside the astropause that likely arises in shocked ISM material. We have used these data, together with published mass-loss rates and wind expansion velocities, to determine the total mass lost and duration for each source -- we find that the duration of and total mass in the shocked wind are significantly larger than their corresponding values for the unshocked wind. The combination of FUV and far-IR data on AGB astrospheres, provides a unique database for theoretical studies (numerical simulations) of wind-ISM interactions. We show that a Cyclical Spatial Heterodyne Spectrometer on a small space-based telescope, can provide high-resolution spectra of astrospheres to confirm the emission mechanism.

We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino-matter interactions and is well-suited to studying the resulting nucleosynthesis and kilonova emission. We find that the ejecta in our simulations under-produce $r$-process abundances beyond the second $r$-process peak. For sufficiently long-lived remnants, these outflows \textit{alone} can produce blue kilonovae, including the blue kilonova component observed for AT2017gfo.

Our ability to interpret observations of galaxies and trace their stellar, gas, and dust content over cosmic time critically relies on our understanding of how the dust abundance and properties vary with environment. Here, we compute the dust surface density across cosmic times to put novel constraints on simulations of the build-up of dust. We provide observational estimates of the dust surface density consistently measured through depletion methods across a wide range of environments, going from the Milky Way up to $z=5.5$ galaxies. These conservative measurements provide complementary estimates to extinction-based observations. In addition, we introduce the dust surface density distribution function -- in analogy with the cold gas column density distribution functions. We fit a power law of the form: $\log f( \Sigma_{\rm Dust})=-1.92 \times \log \Sigma_{\rm Dust} - 3.65$ which proves slightly steeper than for neutral gas and metal absorbers. This observed relation, which can be computed by simulations predicting resolved dust mass functions through 2D projection, provides new constraints on modern dust models.

We present the dayside thermal emission spectrum of WASP-77Ab from 2.8 -- 5.2 $\mu$m as observed with the NIRSpec instrument on the James Webb Space Telescope (JWST). WASP-77Ab was previously found to have a sub-solar metallicity and a solar carbon-to-oxygen (C/O) ratio from H$_2$O and CO absorption lines detected using high-resolution spectroscopy. By performing atmospheric retrievals on the JWST spectrum assuming chemical equilibrium, we find a sub-solar metallicity [M/H]=$-0.91^{+0.24}_{-0.16}$ and C/O ratio $0.36^{+0.10}_{-0.09}$. We identify H$_2$O and CO and constrain their abundances, and we find no CO$_2$ in the spectrum. The JWST and high-resolution spectroscopy results agree within $\sim1\sigma$ for the metallicity and within 1.8$\sigma$ for the C/O ratio. However, our results fit less well in the picture painted by the shorter wavelength spectrum measured by HST WFC3. Comparing the JWST thermal emission spectra of WASP-77Ab and HD 149026b shows that both hot Jupiters have nearly identical brightness temperatures in the near-infrared, but distinctly different atmospheric compositions. Our results reaffirm high-resolution spectroscopy as a powerful and reliable method to measure molecular abundances. Our results also highlight the incredible diversity of hot Jupiter atmospheric compositions.

The emission from shock breakouts (SBOs) represents the earliest electromagnetic (EM) signal emitted by cataclysmic events involving the formation or the merger of neutron stars (NSs). As such, SBOs carry unique information on the structure of their progenitors and on the explosion energy. The characteristic~SBO emission is expected in the UV range, and its detection is one of the key targets of~the ULTRASAT satellite. Among SBO sources, we focus on a specific class involving the formation of fast spinning magnetars in the core-collapse (CC) of massive stars. Fast spinning magnetars are expected to produce a specific signature in the early UV supernova light curve, powered by the extra spin energy quickly released by the NS. Moreover, they are considered as optimal candidates for the emission of long-transient gravitational wave (GW) signals, the detection of which requires early EM triggers to boost the sensitivity of dedicated GW search pipelines. We calculate early supernova UV light curves in the presence of a magnetar central engine, as a function of the explosion energy, ejecta mass and magnetar parameters. We then estimate the ULTRASAT detection horizon (z < 0.15) as a function of the same physical parameters, and the overall expected detection rate finding that magnetar-powered SBOs may represent up to 1/5 of the total events detected by ULTRASAT. Moreover, at the expected sensitivity of the LIGO/Virgo/Kagra O5 science run, one such event occurring within 5 Mpc will providean ideal trigger for a GW long transient search. Future GW detectors like the Einstein Telescope will push the horizon for joint EM-GW detections to 35-40 Mpc.

The Next Generation Arecibo Telescope (NGAT) was a concept presented in a white paper Roshi et al. (2021) developed by members of the Arecibo staff and user community immediately after the collapse of the 305 m legacy telescope. A phased array of small parabolic antennas placed on a tiltable plate-like structure forms the basis of the NGAT concept. The phased array would function both as a transmitter and as a receiver. This envisioned state of the art instrument would offer capabilities for three research fields, viz. radio astronomy, planetary and space & atmospheric sciences. The proposed structure could be a single plate or a set of closely spaced segments, and in either case it would have an equivalent collecting area of a parabolic dish of size 300 m. In this study we investigate the feasibility of realizing the structure. Our analysis shows that, although a single structure ~300 m in size is achievable, a scientifically competitive instrument 130 to 175 m in size can be developed in a more cost effective manner. We then present an antenna configuration consisting of one hundred and two 13 m diameter dishes. The diameter of an equivalent collecting area single dish would be ~130 m, and the size of the structure would be ~146 m. The weight of the structure is estimated to be 4300 tons which would be 53% of the weight of the Green Bank Telescope. We refer to this configuration as NGAT-130. We present the performance of the NGAT-130 and show that it surpasses all other radar and single dish facilities. Finally, we briefly discuss its competitiveness for radio astronomy, planetary and space & atmospheric science applications.

Using deep near-infrared Keck/MOSFIRE observations, we analyze the rest-optical spectra of eight star-forming galaxies in the COSMOS and GOODS-N fields. We reach integration times of $\sim$10 hours in the deepest bands, pushing the limits on current ground-based observational capabilities. The targets fall into two redshift bins -- 5 galaxies at $z \sim 1.7$ and 3 at $z \sim 2.5$ -- and were selected as likely to yield significant auroral-line detections. Even with long integration times, detection of the auroral lines remains challenging. We stack the spectra together into subsets based on redshift, improving the signal-to-noise ratio on the [O III] $\lambda 4364$ auroral emission line and, in turn, enabling a direct measurement of the oxygen abundance for each stack. We compare these measurements to commonly-employed strong-line ratios alongside measurements from the literature. We find that the stacks fall within the distribution of $z>1$ literature measurements, but a larger sample size is needed to robustly constrain the relationships between strong-line ratios and oxygen abundance at high redshift. We additionally report detections of [O I] $\lambda6302$ for eight individual galaxies and composite spectra of 21 targets in the MOSFIRE pointings. We plot their line ratios on the [O III] $\lambda 5008$/H$\beta$ vs. [O I] $\lambda 6302$/H$\alpha$ diagnostic BPT diagram, comparing our targets to local galaxies and H II regions. We find that the [O I]/H$\alpha$ ratios in our sample of galaxies are consistent with being produced in gas ionized by $\alpha$-enhanced massive stars, as has been previously inferred for rapidly-forming galaxies at early cosmic times.

We present new $^{13}$CO(1-0), C$^{18}$O(1-0), HCO$^{+}$(1-0) and H$^{13}$CO$^{+}$(1-0) maps from the IRAM 30m telescope, and a spectrally-resolved [CII] 158 $\mu$m map observed with the SOFIA telescope towards the massive DR21 cloud. This traces the kinematics from low- to high-density gas in the cloud which allows to constrain the formation scenario of the high-mass star forming DR21 ridge. The molecular line data reveals that the sub-filaments are systematically redshifted relative to the dense ridge. We demonstrate that [CII] unveils the surrounding CO-poor gas of the dense filaments in the DR21 cloud. We also show that this surrounding gas is organized in a flattened cloud with curved redshifted dynamics perpendicular to the ridge. The sub-filaments thus form in this curved and flattened mass reservoir. A virial analysis of the different lines indicates that self-gravity should drive the evolution of the ridge and surrounding cloud. Combining all results we propose that bending of the magnetic field, due to the interaction with a mostly atomic colliding cloud, explains the velocity field and resulting mass accretion on the ridge. This is remarkably similar to what was found for at least two nearby low-mass filaments. We tentatively propose that this scenario might be a widespread mechanism to initiate star formation in the Milky Way. However, in contrast to low-mass clouds, gravitational collapse plays a role on the pc scale of the DR21 ridge because of the higher density. This allows more effective mass collection at the centers of collapse and should facilitate massive cluster formation.

We simultaneously measured the spatially-resolved CO-to-H$_{2}$ conversion factor ($\alpha_\mathrm{CO}$) and dust-to-gas ratio (DGR) in nearby galaxies on a kiloparsec scale. In this study, we used $^{12}$CO($J=1-0$) data obtained by the Nobeyama 45-m radio telescope with HI and dust mass surface densities. We obtained the values of global $\alpha_\mathrm{CO}$ and DGR in 22 nearby spiral galaxies, with averages of $2.66 \pm 1.36\ M_\odot\ \mathrm{pc}^{-2}\ (\mathrm{K\ km\ s^{-1}})^{-1}$ and $0.0052 \pm 0.0026$, respectively. Furthermore, the radial variations of $\alpha_\mathrm{CO}$ and DGR in four barred spiral galaxies (IC 342, NGC 3627, NGC 5236, and NGC 6946) were obtained by dividing them into the inner and outer regions with a boundary of $0.2R_{25}$, where $R_{25}$ is the isophotal radius at 25 mag arcsec$^{-2}$ in the $B$ band. The averages of $\alpha_\mathrm{CO}$ and DGR in the inner region ($\leq 0.2R_{25}$) are $0.36 \pm 0.08\ M_\odot\ \mathrm{pc}^{-2}\ (\mathrm{K\ km\ s^{-1}})^{-1}$ and $0.0199 \pm 0.0058$, while those in the outer region ($> 0.2R_{25}$) are $1.49 \pm 0.76\ M_\odot\ \mathrm{pc}^{-2}\ (\mathrm{K\ km\ s^{-1}})^{-1}$ and $0.0084 \pm 0.0037$, respectively. The value of $\alpha_\mathrm{CO}$ in the outer region is 2.3 to 5.3 times larger than that of the inner region. When separated into the inner and outer regions, we find that $\alpha_\mathrm{CO}$ and DGR correlate with the metallicity and the star formation rate surface density. The value of $\alpha_\mathrm{CO}$ derived in this study tends to be smaller than those obtained in previous studies for the Milky Way and nearby star-forming galaxies. This fact can be attributed to our measurements being biased toward the inner region; we measured $\alpha_\mathrm{CO}$ at 0.85 and 0.76 times smaller in radius than the previous works for nearby star-forming galaxies and the Milky Way, respectively.

We present the surface magnetic field conditions of the brightest pulsating RV Tauri star, R Sct. Our investigation is based on the longest spectropolarimetric survey ever performed on this variable star. The analysis of high resolution spectra and circular polarization data give sharp information on the dynamics of the atmosphere and the surface magnetism, respectively. Our analysis shows that surface magnetic field can be detected at different phases along a pulsating cycle, and that it may be related to the presence of a radiative shock wave periodically emerging out of the photosphere and propagating throughout the stellar atmosphere.

We study a Jupiter-mass planet formation for the first time in radiative magneto-hydrodynamics (MHD) simulations and compare it with pure hydrodynamical simulations, as well as to different isothermal configurations. We found that the meridional circulation is the same in every setup. The planetary spiral wakes drive a vertical stirring inside the protoplanetary disc and the encounter with these shock fronts also helps in delivering gas vertically onto the Hill-sphere. The accretion dynamics are unchanged: the planet accretes vertically, and there is outflow in the midplane regions inside the Hill-sphere. We determined the effective $\alpha$-viscosity generated in the disc by the various angular momentum loss mechanisms, which showed that magnetic fields produce high turbulence in the ideal MHD limit, that grows from $\alpha \sim 10^{-2.5}$ up to $\sim 10^{-1.5}$ after the planet spirals develop. In the HD simulations, the planetary spirals contribute to $\alpha \sim 10^{-3}$, making this a very important angular momentum transport mechanism. Due to the various $\alpha$ values in the different setups, the gap opening is different in each case. In the radiative MHD setups, the high turbulent viscosity prevents gap opening, leading to a higher Hill mass, and no clear dust trapping regions. While the Hill accretion rate is $10^{-6} \rm{M_{Jup}/yr}$ in all setups, the accretion variability is orders of magnitude higher in radiative runs than in isothermal ones. Finally, with higher-resolution runs, the magneto-rotational instability started to be resolved, changing the effective viscosity and increasing the heating in the disc.

The growing range of automated algorithms for the identification of molecular clouds and clumps in large observational datasets has prompted the need for the direct comparison of these procedures. However, these methods are complex and testing for biases is often problematic: only a few of them have been applied to the same data set or calibrated against a common standard. We compare the Fellwalker method, a widely used watershed algorithm, to the more recent Spectral Clustering for Interstellar Molecular Emission Segmentation (SCIMES). SCIMES overcomes sensitivity and resolution biases that plague many friends-of-friends algorithms by recasting cloud segmentation as a clustering problem. Considering the \ce{^{13}CO}/\ce{C^{18}O} ($J = 3 - 2$) Heterodyne Inner Milky Way Plane Survey (CHIMPS) and the CO High-Resolution Survey (COHRS), we investigate how these two different approaches influence the final cloud decomposition. Although the two methods produce largely similar statistical results over the CHIMPS dataset, FW appears prone to over-segmentation, especially in crowded fields where gas envelopes around dense cores are identified as adjacent, distinct objects. FW catalogue also includes a number of fragmented clouds that appear as different objects in a line-of-sight projection. In addition, cross-correlating the physical properties of individual sources between catalogues is complicated by different definitions, numerical implementations, and design choices within each method, which make it very difficult to establish a one-to-one correspondence between the sources.

We aim to unravel the interplay between bars, star formation (SF), and active galactic nuclei (AGNs) in barred galaxies. To this end, we utilize the SDSS DR12 to select a sample of nearby (0.02 < z < 0.06) disk galaxies that are suitable for bar examination ($M_r < -20.12$ and inclination $\lesssim$ 53$^{\circ}$). We identify 3662 barred galaxies and measure the length and axis ratio of each bar. We invent new bar parameters that mitigate the stellar and bulge mass biases and show, for the first time, that the evolution of non-AGN and AGN-hosting barred galaxies should be tracked using different bar parameters; the bar length for non-AGN galaxies and the bar axis ratio for AGN-hosting galaxies. Our analysis confirms that barred galaxies have a higher specific SF rate than unbarred control galaxies. Moreover, we find a positive correlation of bar length with both the SF enhancement and the centrally star-forming galaxy fraction, indicating the interconnectivity of bars and SF through the bar-driven gas inflow. We also find that while the AGN fraction of barred galaxies is the same as that of the unbarred control sample, galaxies hosting more massive black holes (BHs) have rounder (i.e., higher axis ratio) bars, implying that the bar is not a cause of AGN activity; rather, AGNs appear to regulate bars. Our findings corroborate theoretical predictions that bars in non-AGN galaxies grow in length, and bars in AGN-hosting galaxies become rounder as BHs grow and eventually get destroyed.

Disk winds are thought to play a critical role in the evolution and dispersal of protoplanetary disks. The primary diagnostic of this physics is emission from the wind, especially in the low-velocity component of the [O I] $\lambda6300$ line. However, the interpretation of the line is usually based on spectroscopy alone, which leads to confusion between magnetohydrodynamic winds and photoevaporative winds. Here, we report that in high-resolution VLT/MUSE spectral mapping of TW~Hya, 80 % of the [O ] emission is confined to within 1 AU radially from the star. A generic model of a magnetothermal wind produces [O I] emission at the base of the wind that broadly matches the flux and the observed spatial and spectral profiles. The emission at large radii is much fainter that predicted from models of photoevaporation, perhaps because the magnetothermal wind partially shields the outer disk from energetic radiation from the central star. This result calls into question the assumed importance of photoevaporation in disk dispersal predicted by models.

Exciting clues to isotropic cosmic birefringence have recently been detected in the $EB$ cross-power spectra of the polarization data of the cosmic microwave background (CMB). Early Dark Energy (EDE) models with a pseudoscalar field coupled to photons via a Chern-Simons term can be used to explain this phenomenon, and can also potentially be used to simultaneously resolve the $H_0$ tension. In this work we incorporate an early dark energy scalar field, including a Chern-Simons coupling, into an existing Boltzmann solver and numerically recover the $EB$ cross-power spectrum for two models in the literature; the $\alpha$-attractor, and the Rock `n' Roll field. We find that both the models fit the $EB$ spectra, and the $EB$ spectra alone do not possess sufficient constraining power to distinguish the two models based on current data.

Galaxy surveys map the three-dimensional distribution of matter in the Universe, encoding information about both the primordial cosmos and its subsequent evolution. By comparing the angular and physical scales of features in the galaxy distribution, we can compute the physical distance to the sample, and thus extract the Hubble parameter, $H_0$. In this chapter, we discuss how this is performed in practice, introducing two key ``standard rulers''. The first, the sound horizon at recombination, leads to baryon acoustic oscillations, and, by combining with external data from the CMB or Big Bang Nucleosynthesis, leads to a competitive $H_0$ constraint. Information can also be extracted from the physical scale of the horizon at matter-radiation equality; though somewhat less constraining, this depends on very different physics and is an important validation test of the physical model. We discuss how both such constraints can be derived (using ``template'' and ``full-shape'' methodologies), and present a number of recent constraints from the literature, some of which are comparable in precision to (and independent from) Planck. Finally, we discuss future prospects for improving these constraints in the future.

The recent detection of the Geminga PWN by HAWC in the multi-TeV band allows us to infer precious information about the transport of pairs in the immediate surroundings of the pulsar and on the spectrum of pairs contributed by a Geminga-like pulsar to the spectrum of pairs in the cosmic radiation. Moreover, this detection allows us to address the issue of how typical are the so-called TeV halos associated to PWNe. Our calculations confirm the need to have suppressed diffusion in a region of at least $20-50\,$pc around the pulsar, and are used here to infer precious constraints on the spectrum of pairs accelerated at the termination shock: more specifically, we discuss the conditions under which such a spectrum is consistent with that typically expected in a PWN and how it gets modified once it escapes the halo. Finally, we discuss the implications of the existence of a TeV halo around Geminga in terms of acceleration of protons in the pulsar environment, a topic of profound relevance for the whole field of particle acceleration and physics of pulsars.

Cosmology constraints serve as a crucial criterion in discriminating cosmological models. The traditional combined method to constrain the cosmological parameters designates the corresponding theoretical value and observational data as functions of redshift, however, sometimes the redshift cannot be measured directly, or the measurement error is large, or the definition of redshift is controversial. In this paper, we propose a novel joint method to constrain parameters that eliminates the redshift $z$ and makes full use of the multiple observables $\left\lbrace \mathcal{F}_{1,\mathrm{obs}},\mathcal{F}_{2,\mathrm{obs}},\cdots,\mathcal{F}_{M,\mathrm{obs}}\right\rbrace$ spanning in $M$-dimensional joint observables space. Considering the generality of the mathematical form of the cosmological models and the guidance from low to high dimensions, we firstly validate our method in a three-dimensional joint observables space spanned by $H(z)$, $f\sigma_{8}(z)$ and $D_{A}(z)$, where the three coordinates can be considered redshift-free measurements of the same celestial body (or shared-redshift data reconstructed model independently). Our results are consistent with the traditional combined method but with lower errors, yielding $H_0=68.7\pm0.1\mathrm{~km} \mathrm{~s}^{-1}\mathrm{~Mpc}^{-1}$, $\Omega_{m0}=0.289\pm0.003$, $\sigma_{8}=0.82\pm0.01$ and showing alleviated parametric degeneracies to some extent. In principle, our joint constraint method allows an extended form keeping the redshift information as an independent coordinate and can also be readily degraded to the form of a traditional combined method to constrain parameters.

Tidal disruption events are common in the Universe, which may occur in various compact star systems and could account for many astrophysical phenomena. Depending on the separation between the central compact star and its companion, either a full disruption or a partial disruption may occur. The partial disruption of a rocky planet around a neutron star can produce kilometer-sized clumps, but the main portion of the planet can survive. The dynamical evolution of these clumps is still poorly understood. In this study, the characteristics of partial disruption of a rocky planet in a highly elliptical orbit around a neutron star is investigated. The periastron of the planet is assumed to be very close to the neutron star so that it would be partially disrupted by tidal force every time it passes through the periastron. It is found that the fragments generated in the process will change their orbits on a time scale of a few orbital periods due to the combined influence of the neutron star and the remnant planet, and will finally collide with the central neutron star. Possible outcomes of the collisions are discussed.

We investigate the environmental dependence of the gas-phase metallicity for galaxies at $z=0$ to $z\gtrsim 2$ and the underlying physical mechanisms driving this dependence using state-of-the-art cosmological hydrodynamical simulations. We find that, at fixed stellar mass, central galaxies in massive halos have lower gas-phase metallicity than those in low-mass halos. On the contrary, satellite galaxies residing in more massive halos are more metal-rich. The combined effect is that massive galaxies are more metal-poor in massive halos, and low-mass galaxies are more metal-rich in massive halos. By inspecting the environmental dependence of other galaxy properties, we identify that the accretion of low-metallicity gas is responsible for the environmental dependence of central galaxies at high $z$, whereas the AGN feedback processes play a crucial role at low $z$. For satellite galaxies, we find that both the suppression of gas accretion and the stripping of existing gas are responsible for their environmental dependence, with negligible effect from the AGN feedback. Finally, we show that the difference of gas-phase metallicity as a function of stellar mass between protocluster and field galaxies agrees with recent observational results, for example from the MAMMOTH-Grism survey.

The linear stability analysis of a stratified rotating fluid (see paper I) showed that disks with a baroclinic stratification under the influence of thermal relaxation will become unstable to thermal instabilities. One instability is the Goldreich-Schubert-Fricke instability (GSF), which is the local version of the Vertical Shear Instability (VSI) and the other is a thermal overstability, the Convective Overstability (COS). In the present paper we reproduce the analytic predicted growth rates for both instabilities in numerical experiments of small axisymmetric sections of vertically isothermal disks with a radial temperature gradient, especially for cooling times longer than the critical cooling time for VSI. In this cooling time regime our simulations reveal the simultaneous and independent growth of both modes: COS and GSF. We consistently observe that GSF modes exhibit a faster growth rate compared to COS modes. Near the midplane, GSF modes eventually stop growing, while COS modes continue to grow and ultimately dominate the flow pattern. Away from the midplane, we find GSF modes to saturate, when bands of constant angular momentum have formed. In these bands we observe the formation and growth of eddies driven by the baroclinic term, further enhancing the velocity perturbations. In geophysics this effect is known as horizontal convection or sea-breeze instability. Three-dimensional simulations will have to show whether similar effects will occur when axisymmetry is not enforced. Our local simulations help to reveal the numerical resolution requirements to observe thermal instabilities in global simulations of disks around young stars.

Aims: We aim at spatially and spectrally resolving the innermost scale of the young stellar object CI Tau to constrain the inner disk properties and better understand the magnetospheric accretion phenomenon. Methods: The high sensitivity offered by the combination of the four 8-m telescopes of the VLTI allied with the spectral resolution of the K-band beam combiner GRAVITY offers a unique capability to probe the sub-au scale of the CI Tau system, tracing both dust and gas emission regions. We develop a geometrical model to fit the interferometric observables and constrain the physical properties of the inner dusty disk. The continuum-corrected pure line visibilities have been used to estimate the size of the Br$\gamma$ emitting region. Results: From the K-band continuum study, we report an highly inclined resolved inner dusty disk, with an inner edge located at a distance of $21\pm2\,R_\star$ from the central star, which is significantly larger than the dust sublimation radius (R$_{sub}= 4.3$ to $8.6\,R_\star$). The inner disk appears misaligned compared to the outer disk observed by ALMA and the non-zero closure phase indicates the presence of a bright asymmetry on the south-west side. From the differential visibilities across the Br$\gamma$ line, we resolve the line emitting region, and measure a size of $4.8^{+0.8}_{-1.0}$ $R_\star$. Conclusions: The extended inner disk edge compared to the dust sublimation radius is consistent with the claim of an inner planet, CI Tau b, orbiting close-in. The inner-outer disk misalignment may be induced by gravitational torques or magnetic warping. The size of the Br$\gamma$ emitting region is consistent with the magnetospheric accretion process. Assuming it corresponds to the magnetospheric radius, it is significantly smaller than the co-rotation radius, which suggests an unstable accretion regime that is consistent with CI Tau being a burster.

A sample of X-ray detected reverberation-mapped quasars provides a unique opportunity to compare cosmological constraints inferred using two well-established relations - the X-ray/UV luminosity ($L_{X}-L_{UV}$) relation and the broad-line region radius-UV monochromatic luminosity ($R-L$) relation. $L_{X}-L_{UV}$ and $R-L$ luminosity distances to the same quasars exhibit a distribution of their differences that is generally positively skewed for the six cosmological models we consider. This behaviour can be interpreted qualitatively to arise as a result of the dust extinction of UV/X-ray quasar emission. We show that the extinction always contributes to the non-zero difference between $L_{X}-L_{UV}$-based and $R-L$-based luminosity distances and we derive a linear relationship between the X-ray/UV colour index $E_{X-UV}$ and the median/mean value of the luminosity-distance difference, which also depends on the value of the $L_{X}-L_{UV}$ relation slope. Taking into account the prevailing positive values of the luminosity-distance difference median, we estimate an average X-ray/UV colour index of $\overline{E}_{X-UV}=0.089 \pm 0.019$ mag, while the value based on the positive mean values of the difference is $\overline{E}_{X-UV}=0.050\pm 0.013$ mag. We demonstrate that this amount of extinction is typical for the majority of quasars since it originates in the circumnuclear and interstellar media of host galaxies. It can only be slightly alleviated by the standard hard X-ray and far-UV extinction cuts used by Lusso et al. (2020). Consequently, the $L_{X}-L_{UV}$ relation QSO data compilation of Lusso et al. (2020) cannot be used for cosmological purposes.

We couple the DELPHI framework for galaxy formation with a model for the escape of ionizing photons to study both its variability with galaxy assembly and the resulting key reionization sources. In this model, leakage either occurs through a fully ionized gas distribution (ionization bounded) or additionally through channels cleared of gas by supernova explosions (ionization bounded + holes). The escape fraction is therefore governed by a combination of the density and star formation rate. Having calibrated our star formation efficiencies to match high-$z$ observables, we find the central gas density to regulate the boundary between high ($>0.70$) and low ($<0.06$) escape fractions. As galaxies become denser at higher redshifts, this boundary shifts from $M_{h}\simeq 10^{9.5}\mathrm{M_{\odot}}$ at $z\sim 5$ to $M_{h}\simeq 10^{7.8}\mathrm{M_{\odot}}$ at $z\sim 15$. While leakage is entirely governed through holes above this mass range, it is not affecting general trends for lower masses. We find the co-evolution of galaxy assembly and the degree of leakage to be mass and redshift dependent, driven by an increasing fraction of $f_{\mathrm{esc}}<0.06$ galaxies at increasing mass and redshift. The variability in the escape of ionizing photons is driven by the underlying variations in our dark matter assembly histories. Galaxies with $M_h < 10^{7.9} ~ (10^{8.9})M_{\odot}$ provide half of the escaping ionizing emissivity by $z \sim 10 ~ (5)$ in the ionization bounded model. On the other hand, galaxies that purely leak through holes contribute $6$ $(13)\%$ at $z\sim 5$ $(15)$. We end by exploring the impact of two reionization feedback scenarios, in which we suppress the gas content of galaxies with $T_{\mathrm{vir}}<20000\mathrm{K}$ and $v_{c}<30\mathrm{kms^{-1}}$ residing in ionized regions.

This paper deals with the problem of constructing a flight scheme to Venus, in which a spacecraft flying to the planet after a gravity assist maneuver and transition to a resonant orbit in order to re-encounter with Venus, makes a passage of a minor celestial body. The 117 candidate asteroids from the NASA JPL catalogue, whose diameter exceeds 1 km, were selected. The flight trajectories which meet the criteria of impulse-free both flyby Venus and asteroid, and the subsequent landing on the surface of Venus were found within the interval of launch dates from 2029 to 2050. The trajectory of the spacecraft flight from the Earth to Venus including flyby of Venus and asteroids with a subsequent landing on the surface of Venus was analyzed.

In Kato et al. (2019, arXiv:1909.00910), I reported on a double outburst and rebrightenings in 2018 in V544 Her. Such a phenomenon is usually observed in WZ Sge stars which evolved after the period bounce and the colors of V544 Her in quiescence apparently exclude this possibility. Although this phenomenon was considered to be rare, I detected almost exactly the same one in 2021 using ZTF, ATLAS and ASAS-SN public data. I also detected a phenomenon very similar to this in ASASSN-19yt in 2022. The same object showed a different type of outburst in 2019 whose morphology looked like that of an SS Cyg star. If ASASSN-19yt is an SU UMa star, the morphology of the 2019 outburst would challenge our knowledge in SU UMa stars. If this object, or V544 Her, is an SS Cyg star, what causes a double outburst and rebrightenings would become an unsolved problem in dwarf novae.

We here report a probable detection of a stellar coronal mass ejection (CME) in active M dwarf KIC 8093473 by performing an analysis on its time resolved X-ray spectra observed by XMM-Newton satellite. Compared to the value at quiescent state and the interstellar one, our spectral modeling returns a marginal (and probably evolving) excess of hydrogen column density in the flare state at a significance level of 1$\sigma$, which can be understood by an additional absorption due to a flare-associate CME. The CME mass is then estimated to be $\sim7\times10^{18}-2\times10^{20}$ g according to the ice cream cone model.

In this paper, we carry out a detailed analysis of the M1.6 class eruptive flare occurring in NOAA active region 13078 on 2022 August 19. The flare is associated with a fast coronal mass ejection (CME) propagating in the southwest direction with an apparent speed of $\sim$926 km s$^{-1}$. Meanwhile, a shock wave is driven by the CME at the flank. The eruption of CME generates an extreme-ultraviolet (EUV) wave expanding outward from the flare site with an apparent speed of $\geq$200 km s$^{-1}$. As the EUV wave propagates eastward, it encounters and interacts with the low-lying adjacent coronal loops (ACLs), which are composed of two loops. The compression of EUV wave results in contraction, expansion, and transverse vertical oscillations of ACLs. The commencements of contraction are sequential from western to eastern footpoints and the contraction lasts for $\sim$15 minutes. The speeds of contraction lie in the range of 13$-$40 km s$^{-1}$ in 171 {\AA} and 8$-$54 km s$^{-1}$ in 193 {\AA}. A long, gradual expansion follows the contraction at lower speeds. Concurrent vertical oscillations are superposed on contraction and expansion of ACLs. The oscillations last for 2$-$9 cycles and the amplitudes are $\leq$4 Mm. The periods are between 3 to 12 minutes with an average value of 6.7 minutes. The results show rich dynamics of coronal loops.

Gravitational wave observations indicate the existence of merging black holes (BHs) with high spin ($a\gtrsim0.3$), whose formation pathways are still an open question. A possible way to form those binaries is through the tidal spin-up of a Wolf-Rayet (WR) star by its BH companion. In this work, we investigate this scenario by directly calculating the tidal excitation of oscillation modes in WR star models, determining the tidal spin-up rate, and integrating the coupled spin-orbit evolution for WR-BH binaries. We find that for short-period orbits and massive WR stars, the tidal interaction is mostly contributed by standing gravity modes, in contrast to Zahn's model of travelling waves which is frequently assumed in the literature. The standing modes are less efficiently damped than traveling waves, meaning that prior estimates of tidal spin-up may be overestimated. We show that tidal synchronization is rarely reached in WR-BH binaries, and the resulting BH spins have $a \lesssim 0.4$ for all but the shortest period ($P_{\rm orb} \! \lesssim 0.5 \, {\rm d}$) binaries. Tidal spin-up in lower-mass systems is more efficient, providing an anti-correlation between the mass and spin of the BHs, which could be tested in future gravitational wave data. Non-linear damping processes are poorly understood but may allow for more efficient tidal spin-up. We also discuss a new class of gravito-thermal modes that appear in our calculations.

French astronomer Honor\'e Flaugergues compiled astronomical observations in a series of hand-written notebooks for 1782$\unicode{x2013}$1830, which are preserved at Paris Observatory. We reviewed these manuscripts and encoded the records that contain sunspot measurements into a numerical table for further analysis. All measurements are timings and we found three types of measurements allowing the reconstruction of heliographic coordinates. In the first case, the Sun and sunspots cross vertical and horizontal wires, in the second case, one vertical and two mirror-symmetric oblique wires, and in the third case, a rhombus-shaped set of wires. Additionally, timings of two solar eclipses also provided a few sunspot coordinates. As a result, we present the time--latitude (butterfly) diagram of the reconstructed sunspot coordinates, which covers the period of the Dalton Minimum and confirms consistency with those of Derfflinger and Prantner. We identify four solar cycles in this diagram and discuss the observed peculiarities as well as the data reliability.

On December 2021, a new camera box for two-colour simultaneous visible photometry was successfully installed on the ASTEP telescope at the Concordia station in Antarctica. The new focal box offers increased capabilities for the ASTEP+ project. The opto-mechanical design of the camera was described in a previous paper. Here, we focus on the laboratory tests of each of the two cameras, the low-temperature behaviour of the focal box in a thermal chamber, the on-site installation and alignment of the new focal box on the telescope, the measurement of the turbulence in the tube and the operation of the telescope equipped with the new focal box. We also describe the data acquisition and the telescope guiding procedure and provide a first assessment of the performances reached during the first part of the 2022 observation campaign. Observations of the WASP19 field, already observed previously with ASTEP, demonstrates an improvement of the SNR by a factor 1.7, coherent with an increased number of photon by a factor of 3. The throughput of the two cameras is assessed both by calculation of the characteristics of the optics and quantum efficiency of the cameras, and by direct observations on the sky. We find that the ASTEP+ two-colour transmission curves (with a dichroic separating the fluxes at 690nm) are similar to those of GAIA in the blue and red channels, but with a lower transmission in the ASTEP+ red channel leading to a 1.5 magnitude higher B-R value compared to the GAIA B-R value. With this new setting, the ASTEP+ telescope will ensure the follow-up and the characterization of a large number of exoplanetary transits in the coming years in view of the future space missions JWST and Ariel.

Blue stragglers are anomalously luminous core hydrogen-burning stars formed through mass-transfer in binary/triple systems and stellar collisions. Their physical and evolutionary properties are largely unknown and unconstrained. Here we analyze 320 high-resolution spectra of blue stragglers collected in eight galactic globular clusters with different structural characteristics and show evidence that the fraction of fast rotating blue stragglers (with rotational velocities larger than 40 km/s) increases for decreasing central density of the host system. This trend suggests that fast spinning blue stragglers prefer low-density environments and promises to open an unexplored route towards understanding the evolutionary processes of these stars. Since large rotation rates are expected in the early stages of both formation channels, our results provide direct evidence for recent blue straggler formation activity in low-density environments and put strong constraints on the timescale of the collisional blue straggler slow-down processes.

Previous studies have identified two emission modes in PSR B1859+07: a normal mode that has three prominent components in the average profile, with the trailing one being the brightest, and an anomalous mode (i.e. the A mode) where emissions seem to be shifted to an earlier phase. Within the normal mode, further analysis has revealed the presence of two sub-modes, i.e. the cW mode and cB mode, where the central component can appear either weak or bright. As for the anomalous mode, a new bright component emerges in the advanced phase while the bright trailing component in the normal mode disappears. New observations of PSR B1859+07 by using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) have revealed the existence of a previously unknown emission mode, dubbed as the Af mode. In this mode, all emission components seen in the normal and anomalous modes are detected. Notably, the mean polarization profiles of both the A and Af modes exhibit an orthogonal polarization angle jump in the bright leading component. The polarization angles for the central component in the original normal mode follow two distinct orthogonal polarization modes in the A and Af modes respectively. The polarization angles for the trailing component show almost the same but a small systematic shift in the A and Af modes, roughly following the values for the cW and cB modes. Those polarization features of this newly detected emission mode imply that the anomalous mode A of PSR B1859+07 is not a result of ``phase shift" or ``swooshes" of normal components, but simply a result of the varying intensities of different profile components. Additionally, subpulse drifting has been detected in the leading component of the Af mode.

The Andromeda Galaxy (M31) is the Local Group galaxy that is most similar to the Milky Way (MW). The similarities between the two galaxies make M31 useful for studying integrated properties common to spiral galaxies. We use the data from the recent QUIJOTE-MFI Wide Survey, together with new raster observations focused on M31, to study its integrated emission. The addition of raster data improves the sensitivity of QUIJOTE-MFI maps by a factor greater than 3. Our main interest is to confirm if anomalous microwave emission (AME) is present in M31, as previous studies have suggested. To do so, we built the integrated spectral energy distribution of M31 between 0.408 and 3000 GHz. We then performed a component separation analysis taking into account synchrotron, free-free, AME and thermal dust components. AME in M31 is modelled as a log-normal distribution with maximum amplitude, $A_{\rm AME}$, equal to $1.06\pm0.30$ Jy. It peaks at $\nu_{\rm AME}=17.28\pm3.08$ GHz with a width of $W_{\rm AME}=0.57\pm0.15$. Both the Akaike and Bayesian Information Criteria find the model without AME to be less than 1 % as probable as the one taking AME into consideration, thus strongly favouring the presence of AME in M31. We find that the AME emissivity in M31 is $\epsilon_{\rm AME}^{\rm 28.4\,GHz}=9.1\pm2.9$ $\mu$K/(MJy/sr), similar to that computed for the MW. We also provide the first upper limits for the AME polarization fraction in an extragalactic object. M31 remains the only galaxy where an AME measurement has been made of its integrated spectrum.

We compute the structure of a Newtonian, multi-ion radiation-mediated shock for different compositions anticipated in various stellar explosions, including supernovae, gamma-ray bursts, and binary neutron star mergers, using a multi-fluid RMS model that incorporates a self-consistent treatment of electrostatic coupling between the different plasma constituents. We find a significant velocity separation between ions having different charge-to-mass ratios in the immediate shock downstream and demonstrate that in fast enough shocks ion-ion collisions can trigger fusion and fission events at a relatively large rate. Our analysis does not take into account potential kinetic effects, specifically, anomalous coupling through plasma microturbulence, that can significantly reduce the velocity spread downstream, below the activation energy for nuclear reactions. A rough estimate of the scale separation in RMS suggests that for shocks propagating in BNS merger ejecta, the anomalous coupling length may exceed the radiation length, allowing a considerable composition change behind the shock via inelastic collisions of $\alpha$ particles with heavy elements at shock velocities $\beta_u\gtrsim0.2$. Moreover, a sufficient abundance of free neutrons upstream of the shock can also trigger fission through neutron capture reactions downstream. The resultant change in the composition profile may affect the properties of the early kilonova emission. The implications for other exploding systems are also briefly discussed.

The 21-cm signal from the dark ages provides a potential new probe of fundamental cosmology. While exotic physics could be discovered, here we quantify the expected benefits within the standard cosmology. A measurement of the global (sky-averaged) 21-cm signal to the precision of thermal noise from a 1,000 hour integration would yield a $5.5\%$ measurement of a combination of cosmological parameters. A 10,000 hour integration would improve this to $1.8\%$, and constrain the cosmic Helium fraction as well as Planck. Precision cosmology with 21-cm fluctuations requires a collecting area of $10\,{\rm km}^2$ (which corresponds to 400,000 stations), which with a 1,000 hour integration would exceed the same global case. Enhancing the collecting area or integration time $\times$10 would yield a $0.5\%$ parameter combination, a Helium measurement five times better than Planck, and a constraint on the neutrino mass as good as Planck. Our analysis sets a baseline for upcoming lunar and space-based dark ages experiments.

Transmission spectroscopy is one of the most successful methods of learning about exoplanet atmospheres. The process of retrievals using transmission spectroscopy consists of creating numerous forward models and comparing them to observations to solve the inverse problem of constraining the atmospheric properties of exoplanets. We explore the impact of one simplifying assumption commonly employed by forward models of transiting exoplanets: namely that the planet can be treated as an isolated, non-rotating spherical body. The centrifugal acceleration due to a planet's rotation opposes the gravitational pull on a planet's atmosphere and increases its scale height. Conventional forward models used for retrievals generally do not include this effect. We find that atmospheric retrievals produce significantly different results for close-in planets with low gravity when this assumption is removed, e.g., differences between true and retrieved values of gas abundances greater than 1$\sigma$ for a simulated planet analogous to WASP-19 b. We recommend that the correction to the atmospheric scale height due to this effect be taken into account for the analysis of high precision transmission spectra of exoplanets in the future, most immediately JWST Cycle 1 targets WASP-19 b and WASP-121 b.

The Hunga Tonga-Hunga Ha'apai volcano erupted on 15 January 2022 with an energy equivalent to around 61 megatons of TNT. The explosion was bigger than any other volcanic eruption so far in the 21st century. Huge quantities of particles, including dust and water vapour, were released into the atmosphere. We present the results of a preliminary study of the effects of the explosion on observations taken at Paranal Observatory using a range of instruments. These effects were not immediately transitory in nature, and a year later stunning sunsets are still being seen at Paranal.

In this paper we report on the follow-up of five potential exoplanets detected with Gaia astrometry and provide an overview of what is currently known about the nature of the entire Gaia astrometric exoplanet candidate sample, 72 systems in total. We discuss the primary false-positive scenario for astrometric planet detections: binary systems with alike components that produce small photocenter motions, mimicking exoplanets. These false positives can be identified as double-lined SB2 binaries through analysis of high resolution spectra. Doing so we find that three systems, Gaia DR3 1916454200349735680, Gaia DR3 2052469973468984192, and Gaia DR3 5122670101678217728 are indeed near equal mass double star systems rather than exoplanetary systems. The spectra of the other two analyzed systems, HD 40503 and HIP 66074, are consistent with the exoplanet scenario in that no second set of lines can be found in the time series of publicly available high resolution spectra. However, their Gaia astrometric solutions imply radial-velocity semi-amplitudes $\sim$\,3 (HD 40503) and $\sim$\,15 (HIP 66074) larger than what was observed with ground based spectrographs. The Gaia astrometry orbital solutions and ground-based radial-velocity measurements exhibit inconsistencies in six out of a total of 12 exoplanet candidate systems where such data are available, primarily due to substantial differences between observed ground-based radial-velocity semi-amplitudes and those implied by the Gaia orbits. We investigated various hypotheses as to why this might be the case, and though we found no clear perpetrator, we note that a mismatch in orbital inclination offers the most straightforward explanation.

There are 214 X-ray point-sources ($L_{\rm X}>10^{35} \mathrm{erg/s}$) identified as X-ray binaries (XRBs) in the nearby spiral galaxy M83. Since XRBs are powered by accretion onto a neutron star or a black hole from a companion/donor star these systems are promising progenitors of merging double compact objects (DCOs): black hole - black hole (BH-BH), black hole - neutron star (BH-NS), or neutron star - neutron star (NS-NS) systems. The connection (i.e. XRBs evolving into DCOs) may provide some hints to the yet unanswered question: what is the origin of the LIGO/Virgo/KAGRA mergers? Available observations do not allow to determine what will be the final fate of the XRBs observed in M83. Yet, we can use evolutionary model of isolated binaries to reproduce the population of XRBs in M83 by matching model XRBs numbers/types/luminosities to observations. Knowing the detailed properties of M83 model XRBs (donor/accretor masses, their evolutionary ages and orbits) we follow their evolution to the death of donor stars to check whether any merging DCOs are formed. Although all merging DCOs in our isolated binary evolution model go through the XRB phase (defined as reaching X-ray luminosity from RLOF/wind accretion onto NS/BH above $10^{35}$ erg/s), only very few XRBs evolve to form merging (in Hubble time) DCOs. For M83 with its solar-like metallicity stars and continiuous star-formation we find that only $\sim 1-2\%$ of model XRBs evolve into merging DCOs depending on the adopted evolutionary physics. This is caused by (i) merger of donor star with compact object during common envelope phase, (ii) binary disruption at the supernova explosion of donor star, (iii) formation of a DCO on a wide orbit (merger time longer than Hubble time).

The extremely high brightness temperatures of fast radio bursts (FRBs) imply that the radiation process must be coherent, but the radiation mechanism is still unknown. The observed properties of narrow spectra and polarization distributions could be used to constrain the radiation mechanism of FRBs. In this work, we discuss the implications of the spectra and polarizations of FRBs from the perspective of intrinsic radiation mechanisms. We first analyze the observed relative spectral bandwidth of radio bursts from an FRB repeater. Furthermore, we generally discuss the properties of the spectra and polarization of the radiation mechanisms involving the relativistic particle's perpendicular acceleration, which depends on the relation between the particle's deflection angle $\psi$ and the radiation beaming angle $1/\gamma$. We find that: (1) If the narrow spectra of FRBs are attributed to the intrinsic radiation mechanism of a single particle, the condition of $\gamma\psi\ll1$ would be necessary, in which scenario, the observed number fraction between linearly and circularly polarized bursts of some FRB repeaters might be due to the propagation effects; (2) Coherent process by multiple particles with some special distributions can lead to a narrow spectrum even for the scenario with $\gamma\psi\gg1$; (3) If the observed number fraction between linearly and circularly polarized bursts is attributed to the radiation mechanism with $\gamma\psi\gg1$, the cumulative distributions of the linear and circular polarization degrees would mainly depend on the particle's beaming distribution.

Imaging Air Cherenkov Telescopes (IACTs) are essential to ground-based observations of gamma rays in the GeV to TeV regime. One particular challenge of ground-based gamma-ray astronomy is an effective rejection of the hadronic background. We propose a new deep-learning-based algorithm for classifying images measured using single or multiple Imaging Air Cherenkov Telescopes. We interpret the detected images as a collection of triggered sensors that can be represented by graphs and analyzed by graph convolutional networks. For images cleaned of the light from the night sky, this allows for an efficient algorithm design that bypasses the challenge of sparse images in deep learning approaches based on computer vision techniques such as convolutional neural networks. We investigate different graph network architectures and find a promising performance with improvements to previous machine-learning and deep-learning-based methods.

The detailed simulation of extended air showers (EAS) and their emission of Cherenkov and fluorescence light requires increasing computation time and storage volume with increasing energy of the primary particle. Given these limitations, it is currently challenging to optimize configurations of imaging air Cherenkov telescopes at photon energies beyond approximately 100 TeV. Additionally, the existing simulation frameworks are not capable of capturing the interplay of Cherenkov and fluorescence light emission at large zenith angle distances ($\gtrsim 70^\circ$), where the collection area of Cherenkov telescopes considerably increases. Here, we present EASpy, a framework for the simulation of EAS at large zenith angles using parametrizations for electron-positron distributions. Our proposed approach for the emission of fluorescence and Cherenkov light and the subsequent imaging of these components by Imaging Atmospheric Cherenkov Telescopes (IACTs) aims to provide flexibility and accuracy while at the same time it reduces the computation time considerably compared to full Monte Carlo simulations. We find excellent agreement of the resulting Cherenkov images when comparing results obtained from EASpy with the de-facto standard simulation tool CORSIKA and sim_telarray. In the process of verifying our approach, we have found that air shower images appear wider and longer with increasing impact distance at large zenith angles, an effect that has previously not been noted. We also investigate the distribution of light on the ground for fluorescence and Cherenkov emission and highlight their key differences to distributions at moderate zenith angles.

Coronal mass ejections (CMEs) are violent ejections of magnetized plasma from the Sun, which can trigger geomagnetic storms, endanger satellite operations and destroy electrical infrastructures on the Earth. After systematically searching Sun-as-a-star spectra observed by the Extreme-ultraviolet Variability Experiment (EVE) onboard the Solar Dynamics Observatory (SDO) from May 2010 to May 2022, we identified eight CMEs associated with flares and filament eruptions by analyzing the blue-wing asymmetry of the O III 52.58 nm line profiles. Combined with images simultaneously taken by the 30.4 nm channel of the Atmospheric Imaging Assembly onboard SDO, the full velocity and propagation direction for each of the eight CMEs are derived. We find a strong correlation between geomagnetic indices (Kp and Dst) and the angle between the CME propagation direction and the Sun-Earth line, suggesting that Sun-as-a-star spectroscopic observations at EUV wavelengths can potentially help to improve the prediction accuracy of the geoeffectiveness of CMEs. Moreover, an analysis of synthesized long-exposure Sun-as-a-star spectra implies that it is possible to detect CMEs from other stars through blue-wing asymmetries or blueshifts of spectral lines.

We aim at searching for exoplanets on the whole ESO/VLT-SPHERE archive with improved and unsupervised data analysis algorithm that could allow to detect massive giant planets at 5 au. To prepare, test and optimize our approach, we gathered a sample of twenty four solar-type stars observed with SPHERE using angular and spectral differential imaging modes. We use PACO, a new generation algorithm recently developed, that has been shown to outperform classical methods. We also improve the SPHERE pre-reduction pipeline, and optimize the outputs of PACO to enhance the detection performance. We develop custom built spectral prior libraries to optimize the detection capability of the ASDI mode for both IRDIS and IFS. Compared to previous works conducted with more classical algorithms than PACO, the contrast limits we derived are more reliable and significantly better, especially at short angular separations where a gain by a factor ten is obtained between 0.2 and 0.5 arcsec. Under good observing conditions, planets down to 5 MJup, orbiting at 5 au could be detected around stars within 60 parsec. We identified two exoplanet candidates that require follow-up to test for common proper motion. In this work, we demonstrated on a small sample the benefits of PACO in terms of achievable contrast and of control of the confidence levels. Besides, we have developed custom tools to take full benefits of this algorithm and to quantity the total error budget on the estimated astrometry and photometry. This work paves the way towards an end-to-end, homogeneous, and unsupervised massive re-reduction of archival direct imaging surveys in the quest of new exoJupiters.

The solar corona produces coronal rain, hundreds of times colder and denser material than the surroundings. Coronal rain is known to be deeply linked to coronal heating, but its origin, dynamics, and morphology are still not well understood. The leading theory for its origin is thermal instability (TI) occurring in coronal loops in a state of thermal non-equilibrium (TNE), the TNE-TI scenario. Under steady heating conditions, TNE-TI repeats in cycles, leading to long-period EUV intensity pulsations and periodic coronal rain. In this study, we investigate coronal rain on the large spatial scales of an active region (AR) and over the long temporal scales of EUV intensity pulsations to elucidate its distribution at such scales. We conduct a statistical study of coronal rain observed over an AR off-limb with IRIS and SDO imaging data, spanning chromospheric to transition region (TR) temperatures. The rain is widespread across the AR, irrespective of the loop inclination, and with minimal variation over the 5.45-hour duration of the observation. Most rain has a downward ($87.5\%$) trajectory; however, upward motions ($12.5\%$) are also ubiquitous. The rain dynamics are similar over the observed temperature range, suggesting that the TR and chromospheric emission are co-located on average. The average clump widths and lengths are similar in the SJI channels and wider in the AIA 304 channel. We find ubiquitous long-period EUV intensity pulsations in the AR. Short-term periodicity is found (16 min) linked to the rain appearance, which constitutes a challenge to explain under the TNE-TI scenario.

We study how future Type-Ia supernovae (SNIa) standard candles detected by the Vera C. Rubin Observatory (LSST) can constrain some cosmological models. We use a realistic three-year SNIa simulated dataset generated by the LSST Dark Energy Science Collaboration (DESC) Time Domain pipeline, which includes a mix of spectroscopic and photometrically identified candidates. We combine this data with Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillation (BAO) measurements to estimate the dark energy model parameters for two models- the baseline $\Lambda$CDM and Chevallier-Polarski-Linder (CPL) dark energy parametrization. We compare them with the current constraints obtained from joint analysis of the latest real data from the Pantheon SNIa compilation, CMB from Planck 2018 and BAO. Our analysis finds tighter constraints on the model parameters along with a significant reduction of correlation between $H_0$ and $\sigma_8$. We find that LSST is expected to significantly improve upon the existing SNIa data in the critical analysis of cosmological models.

In this work, we consider four different galaxy populations and two distinct global environments in the local Universe (z $\leq 0.11$) to investigate the evolution of transitional galaxies (such as star-forming spheroids and passive discs) across different environments. Our sample is composed of 3,899 galaxies within the R$_{200}$ radius of 231 clusters and 11,460 field galaxies. We also investigate the impact of the cluster's dynamic state, as well as the galaxy's location in the projected phase space diagram (PPS). We found that although the cluster environment as a whole influences galaxy evolution, the cluster dynamical state does not. Furthermore, star-forming galaxies represent recent cluster arrivals in comparison to passive galaxies (especially in the case of early-types). Among the ETGs, we find that the D$_n(4000)$ and H$_\delta$ parameters indicate a smooth transition between the subpopulations. In particular, for the SF-ETGs, we detect a significant difference between field and cluster galaxies, as a function of stellar mass, for objects with Log $M_*$/M$_{\odot} > 10.5$. Analyzing the color gradient, the results point toward a picture where field galaxies are more likely to follow the monolithic scenario, while the cluster galaxies the hierarchical scenario. In particular, if we split the ETGs into lenticulars and ellipticals, we find that the steeper color gradients are more common for the lenticulars. Finally, our results indicate the need for galaxy pre-processing in smaller groups, before entering clusters.

In the context of a recently proposed contact-binary model of the Kreutz system, all its members are products of the process of cascading fragmentation of the two lobes of the parent, Aristotle's comet of 372 BC. This process presumably began with the lobes' separation from each other near aphelion. However, not every object in a Kreutz-like orbit is a Kreutz sungrazer. Any surviving sungrazer that had split off from the progenitor before the lobes separated, as well as its surviving fragments born in any subsequent tidal or nontidal event, are by definition not members of the Kreutz system. Yet, as parts of the same progenitor, they belong -- as do all Kreutz sungrazers -- to a broader assemblage of related objects, which I refer to as a super-Kreutz system. After estimating the ratio of the number of super-Kreutz members to nonmembers among potential historical sungrazers, I generate representative extended pedigree charts for both the Kreutz system and super-Kreutz system. While the fragmentation paths and relationships among the individual sungrazers or potential sungrazers in the two charts are (with at most a few exceptions) arbitrary, the purpose of the exercise is to suggest that the Kreutz system proper could in effect represent an ultimate deagglomeration stage of the super-Kreutz system.

The orbits of the main satellites of Uranus are expected to slowly drift away owing to tides raised in the planet. As a result, the 5/3 mean motion resonance between Ariel and Umbriel was likely encountered in the past. Previous studies have shown that, in order to prevent entrapment in this resonance, the eccentricities of the satellites must be larger than $\sim 0.01$ at the epoch, which is hard to explain. On the other hand, if the satellites experience some temporary capture and then escape, the inclinations rise to high values that are not observed today. We have revisited this problem both analytically and numerically focussing on the inclination, using a secular two-satellite model with circular orbits. We show that if the inclination of Umbriel was around $0.15^{\circ}$ at the time of the 5/3 resonance encounter, capture can be avoided in about $60\%$ of the cases. Moreover, after the resonance crossing, the inclination of Umbriel drops to a mean value around $0.08^{\circ}$, which is close to the presently observed one. The final inclination of Ariel is distributed between $0.01^{\circ}$ and $0.25^{\circ}$ with a nearly equal probability, which includes the present mean value of $0.02^{\circ}$.

The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements.

We investigate the properties of dark energy halos in models with a nonminimal coupling in the dark sector. We show, using a quasistatic approximation, that a coupling of the mass of dark matter particles to a standard quintessence scalar field $\phi$ generally leads to the formation of dark energy concentrations in and around compact dark matter objects. These are associated with regions where scalar field gradients are large and the dark energy equation of state parameter is close to $-1/3$. We find that the energy and radius of a dark energy halo are approximately given by $E_{\rm halo} \sim \boldsymbol{\beta}^2 \varphi \, m$ and $r_{\rm halo} \sim \sqrt{\boldsymbol{\beta} \,\varphi ({R}/{H})}$, where $\varphi=Gm/(R c^2)$, $m$ and $R$ are, respectively, the mass and radius of the associated dark matter object, $\boldsymbol{\beta} = -d \ln m/d \phi$ is the nonminimal coupling strength parameter, $H$ is the Hubble parameter, $G$ is the gravitational constant, and $c$ is the speed of light in vacuum. We further show that current observational limits on $\boldsymbol{\beta}$ over a wide redshift range lead to stringent constraints on $E_{\rm halo}/m$ and, therefore, on the impact of dark energy halos on the value of the dark energy equation of state parameter. We also briefly comment on potential backreaction effects that may be associated with the breakdown of the quasistatic approximation and determine the regions of parameter space where such a breakdown might be expected to occur.

Understanding the star formation rate (SFR) variability and how it depends on physical properties of galaxies is important for developing and testing the theory of galaxy formation. We investigate how statistical measurements of the extragalactic background light (EBL) can shed light on this topic and complement traditional methods based on observations of individual galaxies. Using semi-empirical models of galaxy evolution and SFR indicators sensitive to different star formation timescales (e.g., H$\alpha$ and UV continuum luminosities), we show that the SFR variability, quantified by the joint probability distribution of the SFR indicators (i.e., the bivariate conditional luminosity function), can be characterized as a function of galaxy mass and redshift through the cross-correlation between deep, near-infrared maps of the EBL and galaxy distributions. As an example, we consider combining upcoming SPHEREx maps of the EBL with galaxy samples from Rubin/LSST. We demonstrate that their cross-correlation over a sky fraction of $f_\mathrm{sky}\sim0.5$ can constrain the joint SFR indicator distribution at high significance up to $z\sim2.5$ for mass-complete samples of galaxies down to $M_{*}\sim10^9\,M_{\odot}$. These constraints not only allow models of different SFR variability to be distinguished, but also provide unique opportunities to investigate physical mechanisms that require large number statistics such as environmental effects. The cross-correlations investigated illustrate the power of combining cosmological surveys to extract information inaccessible from each data set alone, while the large galaxy populations probed capture ensemble-averaged properties beyond the reach of targeted observations towards individual galaxies.

We describe a new mechanism that gives rise to dissipation during cosmic inflation. In the simplest implementation, the mechanism requires the presence of a massive scalar field with a softly-broken global $U(1)$ symmetry, along with the inflaton field. Particle production in this scenario takes place on parametrically sub-horizon scales, at variance with the case of dissipation into gauge fields. Consequently, the backreaction of the produced particles on the inflationary dynamics can be treated in a \textit{local} manner, allowing us to compute their effects analytically. We determine the parametric dependence of the power spectrum which deviates from the usual slow-roll expression. Non-Gaussianities are always sizeable whenever perturbations are generated by the noise induced by dissipation: $f_{\rm NL}^{\rm eq} \gtrsim {O}(10)$.

We study the impact of the ambient fluid on the evolution of collapsing false vacuum bubbles by simulating the dynamics of a coupled bubble-particle system. A significant increase in the mass of the particles across the bubble wall leads to a buildup of those particles inside the false vacuum bubble. We show that the backreaction of the particles on the bubble slows or even reverses the collapse. Consequently, if the particles in the true vacuum become heavier than in the false vacuum, the particle-wall interactions always decrease the compactness that the false vacuum bubbles can reach making their collapse to black holes less likely.

New light scalar degrees of freedom may alleviate the dark matter and dark energy problems, but if coupled to matter, they generally mediate a fifth force. In order for this fifth force to be consistent with existing constraints, it must be suppressed close to matter sources, e.g. through a non-linear screening mechanism. In this work, we investigate the non-relativistic two-body problem in shift-symmetric scalar-tensor theories that exhibit kinetic screening ($k$-mouflage), both numerically and analytically. We develop an approximate scheme, based on a Hodge-Helmholtz decomposition of the Noether current associated to the shift symmetry, allowing for a qualitative insight into the dynamics and yielding results in good agreement with the numerical ones in most of the parameter space. We apply the formalism to polynomial $k$-essence and to Dirac-Born-Infeld (DBI) type theories, as well as to theories that develop ``anti-screening''. In the deep nonlinear regime, we find that the fifth force is screened slightly more efficiently in equal-mass systems than in extreme mass-ratio ones. However, we find that systems with comparable masses also exhibit regions where the screening is ineffective. These descreened spheroidal regions (bubbles) could in principle be probed in the solar system with sufficiently precise space accelerometers.

One of the aims of Aether Scalar Tensor Theory (AeST) is to reproduce the successes of Modified Newtonian Dynamics (MOND) on galactic scales. Indeed, the quasi-static limit of AeST can achieve precisely this, assuming that the vector field $\vec{A}$ vanishes. However, this assumption of a vanishing vector field is often inconsistent. Here, we show how to correctly take into account the vector field and find that the quasi-static limit depends on a model parameter $m_\times$. In the limit $m_\times \to 0$, one recovers the quasi-static limit with a vanishing vector field. In particular, one finds a two-field version of MOND. In the opposite limit, $m_\times \to \infty$, one finds a single-field version of MOND. We show that, in practice, much of the phenomenology of the quasi-static limit depends only very little on the value of $m_\times$. Still, for some observational tests, such as those involving wide binaries, $m_\times$ has percent-level effects that may be important.

The introduction of nonlinear nanophotonic devices to the field of optical frequency comb metrology has enabled new opportunities for low-power and chip-integrated clocks, high-precision frequency synthesis, and broad bandwidth spectroscopy. However, most of these advances remain constrained to the near-infrared region of the spectrum, which has restricted the integration of frequency combs with numerous quantum and atomic systems in the ultraviolet and visible. Here, we overcome this shortcoming with the introduction of multi-segment nanophotonic thin-film lithium niobate (LN) waveguides that combine engineered dispersion and chirped quasi-phase matching for efficient supercontinuum generation via the combination of $\chi^{(2)}$ and $\chi^{(3)}$ nonlinearities. With only 90 pJ of pulse energy at 1550 nm, we achieve gap-free frequency comb coverage spanning 330 to 2400 nm. The conversion efficiency from the near-infrared pump to the UV-Visible region of 350-550 nm is nearly 20%. Harmonic generation via the $\chi^{(2)}$ nonlinearity in the same waveguide directly yields the carrier-envelope offset frequency and a means to verify the comb coherence at wavelengths as short as 350 nm. Our results provide an integrated photonics approach to create visible and UV frequency combs that will impact precision spectroscopy, quantum information processing, and optical clock applications in this important spectral window.

Neutrinos as almost massless particles could mediate long-range forces, known as neutrino forces. In this talk, I will introduce some theoretical aspects of neutrino forces, including why the potential of a neutrino force has the $1/r^{5}$ form and how it may vary under different circumstances. Experimental probes and possible implications for cosmology are also briefly discussed.

It is well-known that the primordial scalar curvature and tensor perturbations, $\zeta$ and $\gamma_{ij}$, are conserved on super-horizon scales in minimal inflation models. However, their wave functional has a rapidly oscillating phase which is slow-roll unsuppressed, as can be seen either from boundary (total-derivative) terms of cosmological perturbations, or the WKB approximation of the Wheeler-DeWitt equation. Such an oscillatory phase involves gravitational non-linearity between scalar and tensor perturbations. By tracing out unobserved modes, the oscillatory phase causes faster decoherence of primordial gravitons compared to those by bulk interactions. Our results put a stronger lower bound of decoherence effect to the recent proposals probing squeezed primordial gravitons.

The effect of proportional electroluminescence (EL) is used to record the primary ionization signal (S2) in the gas phase of two-phase argon detectors for dark matter particle (WIMP) searches and low-energy neutrino experiments. Our previous studies of EL time properties revealed the presence of two unusual slow components in S2 signal of two-phase argon detector, with time constants of about 4-5 $\mu$s and 50 $\mu$s. The puzzle of slow components is that their time constants and contributions to the overall signal increase with electric field (starting from a certain threshold), which cannot be explained by any of the known mechanisms of photon and electron emission in two-phase media. There are indications that these slow components result from delayed electrons, temporarily trapped during their drift in the EL gap on metastable negative argon ions of yet unknown nature. In this work, this hypothesis is convincingly confirmed by studying the time properties of electroluminescence in a Thick Gas Electron Multiplier (THGEM) coupled to the EL gap of two-phase argon detector. In particular, an unusual slow component in EL signal, similar to that observed in the EL gap, was observed in THGEM itself. In addition, with the help of THGEM operated in electron multiplication mode, the slow component was observed directly in the charge signal, unambiguously confirming the effect of trapped electrons in S2 signal. These results will help to unravel the puzzle of slow components in two-phase argon detectors and thus to understand the background in low-mass WIMP searches.

Ultralight boson fields, with a mass around $10^{-23}\text{eV}$, are promising candidates for the elusive cosmological dark matter. These fields induce a periodic oscillation of the spacetime metric in the nanohertz frequency band, which is detectable by pulsar timing arrays. In this paper, we investigate the gravitational effect of ultralight tensor dark matter on the arrival time of radio pulses from pulsars. We find that the pulsar timing signal caused by tensor dark matter exhibits a different angular dependence than that by scalar and vector dark matter, making it possible to distinguish the ultralight dark matter signal with different spins. Combining the gravitational effect and the coupling effect of ultralight tensor dark matter with standard model matter provides a complementary way to constrain the coupling parameter $\alpha$. We estimate $\alpha \lesssim 10^{-6}\sim 10^{-5}$ in the mass range $m<5\times 10^{-23}\mathrm{eV}$ with current pulsar timing array.

We present a scalar-driven sterile neutrino production model where the interaction with the ultralight scalar field modifies the oscillation production of sterile neutrinos in the early universe. The model effectively suppresses the production of sterile neutrinos at low temperatures due to the heavy scalar mass, resulting in a colder matter power spectrum that avoids constraints from small-scale structure observations. In this model, the dominant dark matter relic is from sterile neutrinos, with only a small fraction originating from the ultralight scalar. Furthermore, the model predicts a detectable X/Gamma-ray flux proportional to the cubic density of local sterile neutrinos for a light scalar mass due to the light scalar coupling tosterile neutrinos. This distinguishes our model from normal decaying dark matter, which has a linear dependence on the density. In addition, the model predicts a potential low-energy monochromatic neutrino signal that can be detectable by future neutrino telescopes.

We utilized the duality between massive vector and massive Kalb-Ramond fields to derive an effective action for Abelian-Higgs cosmic strings. This enabled us to determine the classically renormalized string tension and facilitate calculations for back-reaction effects. Additionally, we derived a comprehensive expression for the energy flux of radiation emitted by Abelian-Higgs cosmic strings. Applying this equation to a cuspless loop, we obtained that the loop lifetime is proportional to the square of the loop length, which is in agreement with field-theory simulations.

The observation of transient gravitational waves is hindered by the presence of transient noise, colloquially referred to as glitches. These glitches can often be misidentified as gravitational waves by searches for unmodeled transients using the excess-power type of methods and sometimes even excite template waveforms for compact binary coalescences while using matched filter techniques. They thus create a significant background in the searches. This background is more critical in getting identified promptly and efficiently within the context of real-time searches for gravitational-wave transients. Such searches are the ones that have enabled multi-messenger astrophysics with the start of the Advanced LIGO and Advanced Virgo data taking in 2015 and they will continue to enable the field for further discoveries. With this work we propose and demonstrate the use of a signal-based test that quantifies the fidelity of the time-frequency decomposition of the putative signal based on first principles on how astrophysical transients are expected to be registered in the detectors and empirically measuring the instrumental noise. It is based on the Q-transform and a measure of the occupancy of the corresponding time-frequency pixels over select time-frequency volumes; we call it ``QoQ''. Our method shows a 40% reduction in the number of retraction of public alerts that were issued by the LIGO-Virgo-KAGRA collaborations during the third observing run with negligible loss in sensitivity. Receiver Operator Characteristic measurements suggest the method can be used in online and offline searches for transients, reducing their background significantly.

The investigation of the phase state of dense matter is hindered by complications of first-principle nonperturbative quantum chromodynamics. By performing the first consistent general-relativistic calculations of tidal-excited g-mode of neutron stars with a first-order strong interaction phase transition in the high-density core, we demonstrate that gravitational wave signal during binary neutron star inspiral probes their innermost hadron-quark transition and provides potent constraints from present and future gravitational-wave detectors.

We investigate the propagation of certain non-plane wave solutions to Maxwell's equations in both flat and curved spacetimes. We find that the effective signal velocity associated to such solutions need not be $c$ and that the signal need not propagate along null geodesics; indeed, more than this, we find that the information encoded in the signals associated with such solutions can be substantially non-local. Having established these results, we then turn to their conceptual-philosophical-foundational significance -- which, in brief, we take to be the following: (i) one should not assume that all electromagnetic waves generated in the cosmos are localised plane wave packages; thus, (ii) one cannot assume that signals reaching us from the cosmos arrive with a particular velocity (namely, $c$), and that such signals encode local information regarding their sources; therefore (iii) astrophysicists and cosmologists should be wary about making such assumptions in their inferences from obtained data -- for to do so may lead to incorrect inferences regarding the nature of our universe.

We study the theoretical possibility of $^3P_0$ neutron superfluid in dilute spin-polarized neutron matter, which may be relevant to the crust region of a magnetized neutron star. In such a dilute regime where the neutron Fermi energy is less than 1 MeV, the $^1S_0$ neutron superfluid can be exhausted by a strong magnetic field of the compact star. At low-energy limit relevant for dilute neutron matter, the $^3P_0$ interaction is stronger than the $^3P_2$ one, which is believed to induce the triplet superfluid in the core. We present the ground-state phase diagram of dilute neutron matter with respect to the magnetic field and numerically estimate the critical temperature of $^3P_0$ neutron superfluid, which exceeds $10^7$~K.