Papers for Monday, Feb 05 2024

The two known planets in the planetary system of Teegarden's Star are among the most Earth-like exoplanets currently known. Revisiting this nearby planetary system with two planets in the habitable zone aims at a more complete census of planets around very low-mass stars. A significant number of new radial velocity measurements from CARMENES, ESPRESSO, MAROON-X, and HPF, as well as photometry from TESS motivated a deeper search for additional planets. We confirm and refine the orbital parameters of the two know planets Teegarden's Star b and c. We also report the detection of a third planet d with an orbital period of 26.13+-0.04 d and a minimum mass of 0.82+-0.17 M_Earth. A signal at 96 d is attributed to the stellar rotation period. The interpretation of a signal at 172 d remains open. The TESS data exclude transiting short-period planets down to about half an Earth radius. We compare the planetary system architecture of very low-mass stars. In the currently known configuration, the planetary system of Teegarden's star is dynamically quite different from that of TRAPPIST-1, which is more compact, but dynamically similar to others such as GJ 1002.

The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above $10^{14}\rm{M_\odot}$, we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of $\sim 1\rm{kpc}$ in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of $1.48\rm{kpc}$. We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using the the "light traces mass" approach, in this context meaning gas is used as a tracer for the potential, can overestimate offsets by factors of $\sim10$ and $\sim30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section.

Aims. Recent observational studies suggest that feedback from active galactic nuclei (AGNs) may play an important role in the formation and evolution of low-mass dwarf galaxies, an issue that has received little attention from a theoretical perspective. Methods. We investigate this using two sets of 12 cosmological magneto-hydrodynamical simulations of the formation of dwarf galaxies: one set using a version of the AURIGA galaxy formation physics model including AGN feedback and a parallel set with AGN feedback turned off. Results. We show that the full-physics AGN runs satisfactorily reproduce several scaling relations, including the M_ BH-M_*, M_ BH-sigma_* and the baryonic Tully-Fisher relation. We find that the global star formation (SF) of galaxies run with AGN is reduced compared to the one in which AGN has been turned off, suggesting that this type of feedback is a viable way of suppressing SF in low-mass dwarfs, even though none of our galaxies is completely quenched by z=0. Furthermore, we found a tight correlation between the median SF rates and the black-hole-to-stellar mass ratio (M_BH/M_star) in our simulated dwarfs. SF is suppressed due to gas heating in the vicinity of the AGN: less HI gas is available in AGN runs, though the total amount of gas is preserved across the two settings within each galaxy. This indicates that the main effect of AGN feedback in our dwarfs is to heat up and push the gas away from the galaxy's centre rather than expelling it completely. Finally, we show that the two galaxies harbouring the largest BHs have suffered a considerable (up to 65%) reduction in their central dark matter density, pinpointing the role of AGNs in determining the final dark matter mass distribution within dwarf galaxies. This pilot paper highlights the importance of modelling AGN feedback at the lowest mass scales and the impact this can have on dwarf galaxy evolution.

We present a simple method for assessing the predictive performance of high-dimensional models directly in data space when only samples are available. Our approach is to compare the quantiles of observables predicted by a model to those of the observables themselves. In cases where the dimensionality of the observables is large (e.g. multiband galaxy photometry), we advocate that the comparison is made after projection onto a set of principal axes to reduce the dimensionality. We demonstrate our method on a series of two-dimensional examples. We then apply it to results from a state-of-the-art generative model for galaxy photometry (pop-cosmos) that generates predictions of colors and magnitudes by forward simulating from a 16-dimensional distribution of physical parameters represented by a score-based diffusion model. We validate the predictive performance of this model directly in a space of nine broadband colors. Although motivated by this specific example, the techniques we present will be broadly useful for evaluating the performance of flexible, non-parametric population models of this kind, and can be readily applied to any setting where two sets of samples are to be compared.

The Mg II h&k emission lines (2803, 2796 A) are a useful tool for understanding stellar chromospheres and transition regions due to their intrinsic brightness, relatively low interstellar medium (ISM) absorption interference, and abundance of archival spectra available. Similar to other optically thick chromospheric emission lines such as H I Lya, Mg II emissions commonly present with a self-reversed line core, the depth and shape of which vary from star to star. We explore the relationship between self-reversal and the stellar atmosphere by investigating the extent that fundamental stellar parameters affect self-reversal. We present a search for correlations between photospheric parameters such as effective temperature, surface gravity, and metallicity with the Mg II k self-reversal depth for a group of 135 FGKM main sequence stars with high-resolution near-ultraviolet spectra from the Hubble Space Telescope. We modeled the observed MgII k line profiles to correct for ISM attenuation and recover the depth of emission line's self-reversal in relation to the intensity of the line. We used the PHOENIX atmosphere code to homogeneously determine the stellar parameters by computing a suite of stellar atmosphere models that include a chromosphere and transition region, and using archival photometry to guide the models of each star. We quantify the sensitivity of the visible and near-infrared photometry to chromospheric and photospheric parameters. We find weak trends between Mg II k self-reversal depth and age, rotation period, Mg II luminosity, temperature, and mass. All stars in our sample older than ~2 Gyr or rotating slower than ~10 days exhibit self-reversal.

We present pop-cosmos: a comprehensive model characterizing the galaxy population, calibrated to $140,938$ galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in $26$ bands from the ultra-violet to the infra-red. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (parameterized by a flexible score-based diffusion model); a state-of-the-art stellar population synthesis (SPS) model connecting galaxies' instrinsic properties to their photometry; and a data-model for the observation, calibration and selection processes. By minimizing the optimal transport distance between synthetic and real data we are able to jointly fit the population- and data-models, leading to robustly calibrated population-level inferences that account for parameter degeneracies, photometric noise and calibration, and selection effects. We present a number of key predictions from our model of interest for cosmology and galaxy evolution, including the mass function and redshift distribution; the mass-metallicity-redshift and fundamental metallicity relations; the star-forming sequence; the relation between dust attenuation and stellar mass, star formation rate and attenuation-law index; and the relation between gas-ionization and star formation. Our model encodes a comprehensive picture of galaxy evolution that faithfully predicts galaxy colors across a broad redshift ($z<4$) and wavelength range.

Using the 3D smoothed particle hydrodynamics code PHANTOM, we investigate the evolution of the orbital properties of massive black hole binaries embedded in massive discs where gravitational instabilities (GIs) triggered by the disc self-gravity are the only source of angular momentum transport. In particular, we investigate the evolution of binaries with different initial eccentricities $e_0=0.05,\,0.5,\,0.8$ and mass ratios $q=0.1,\,0.3,\,0.9$. Our simulations suggest that there might not be a unique value of critical eccentricity. We find initially more eccentric binaries to tend to higher asymptotic eccentricity values than more circular ones. This implies that there is a range of critical eccentricity values, that depends on the initial condition of the system. In particular, we find the width of this range to be narrower for more unequal binaries. We furthermore measure preferential accretion onto our binaries, finding more accretion onto the primary only for mass ratio $q=0.3$ and eccentricity $e=0.8$. We discuss how this might have implications for the amplitude of the gravitational wave background detected by Pulsar Timing Arrays (PTA) experiments. We finally measure the corresponding value of the viscosity parameter $\alpha$ in our simulations and discuss how this depends on the binary properties.

We study the intermediate tachyon inflation in an anisotropic background. By using the Friedmann equations obtained in the anisotropic geometry, we obtain the slow-roll parameters in the tachyon model. The presence of the anisotropic effects in the slow-roll parameters changes the perturbation parameters in our setup which may change its observational viability. To check this, we perform a numerical analysis and test the results with Planck2018 TT, TE, EE +lowE+lensing+BK14(18)+BAO data. We show that the intermediate anisotropic inflation in some ranges of the anisotropic and intermediate parameters is observationally viable. We also show that the equilateral amplitude of the non-gaussianity in our model is of the order of $10^{-2}-10^{-1}$. By studying the reheating process in our setup, we find that it is possible to have instantaneous reheating in this model. We also find that the temperature during the reheating in our setup is consistent with Big-Bang nucleosynthesis.

We study the spatially resolved outflow properties of CSWA13, an intermediate mass ($M_*=10^{9}~\mathrm{M}_{\odot}$), gravitationally lensed star-forming galaxy at $z=1.87$. We use Keck/KCWI to map outflows in multiple rest-frame ultraviolet ISM absorption lines, along with fluorescent Si II$^*$ emission, and nebular emission from C III] tracing the local systemic velocity. The spatial structure of outflow velocity mirrors that of the nebular kinematics, which we interpret to be a signature of a young galactic wind that is pressurizing the ISM of the galaxy but is yet to burst out. From the radial extent of Si II$^*$ emission, we estimate that the outflow is largely encapsulated within $3.5$ kpc. We explore the geometry (e.g., patchiness) of the outflow by measuring the covering fraction at different velocities, finding that the maximum covering fraction is at velocities $v\simeq-150$ km$\,$s$^{-1}$. Using the outflow velocity ($v_{out}$), radius ($R$), column density ($N$), and solid angle ($\Omega$) based on the covering fraction, we measure the mass loss rate $\log\dot{m}_{out}/(\mathrm{M}_{\odot}\text{yr}^{-1}) = 1.73\pm0.23$ and mass loading factor $\log\eta = 0.04\pm0.34$ for the low-ionization outflowing gas in this galaxy. These values are relatively large and the bulk of the outflowing gas is moving with speeds less than the escape velocity of the galaxy halo, suggesting that the majority of outflowing mass will remain in the circumgalactic medium and/or recycle back into the galaxy. The results support a picture of high outflow rates transporting mass and metals into the inner circumgalactic medium, providing the gas reservoir for future star formation.

The field of sonification, using of non-speech audio for data analysis, is already established in space sciences. Meetings like "The audible Universe" focus on sonification tools applied in astronomy to represent complex data like nebulae and galaxies. Besides, little is said about the translation of images into sound, this challenge that seeks to represent data in 2 or 3 dimensions through a one-dimensional technique. The aforementioned leads to a total lack of consensus regarding the sound parameters to be used and how these are interpreted by people. This work seeks to delve deeper into the existing tools for image sonification, analyzing whether their objective is only outreach or includes the possibility of research. A new proposal is presented, maintaining sonoUno's software focus on research, pointing out the need for reliable techniques that integrate functional diversity people with an active role on research.

The many unique properties of galaxies are shaped by physical processes that affect different components of the galaxy - like the bulges and discs - in different ways, and leave characteristic imprints on the light and spectra of these components. Disentangling their spectra can reveal vital clues that can be traced back in time to understand how galaxies, and their components, form and evolve throughout their lifetimes. With BUDDI, we have decomposed the IFU datacubes in SDSS-MaNGA DR17 into a S\'ersic bulge component and an exponential disc component and extracted their clean bulge and disc spectra. BUDDI-MaNGA is the first and largest statistical sample of such decomposed spectra of 1452 galaxies covering morphologies from ellipticals to late-type spirals. We derived stellar masses of the individual components with SED fitting using BAGPIPES and estimated their mean mass-weighted stellar metallicities and stellar ages using pPXF. With this information in place, we reconstructed the mass assembly histories of the bulges and discs of the 968 spiral galaxies (Sa-Sm Types) in this sample to look for systematic trends with respect to stellar mass and morphology. Our results show a clear downsizing effect especially in the bulges, with more massive components assembling earlier and faster than the less massive ones. Additionally, on comparing the stellar populations of the bulges and discs in these galaxies, we find that a majority of the bulges host more metal-rich and older stars than their disc counterparts. Nevertheless, we also find that there exists a non-negligible fraction of the spiral galaxy population in our sample with bulges that are younger and more metal-rich than their discs. We interpret these results, taking into account how their formation histories and current stellar populations depend on stellar mass and morphology.

The Galileo spacecraft had distant encounters with Earth in 1990 and 1992. Limited Solid State Imager (SSI) data acquired during these encounters has been previously presented, but the majority of the data from these Earth flybys have not been presented in the literature. Observations of Earth taken from afar are both rare and directly relevant to the development of any future exo-Earth direct imaging mission. Here we present a pipeline that vets, calibrates, and measures the disk-integrated brightness of the Earth, in multiple filters, from the complete SSI data sets from both the 1990 and 1992 Galileo flybys. The result is over 1500 usable photometric measurements for Earth as an analog for an exoplanet. The 1990 data set includes full rotational lightcurves in six bandpasses spanning the optical range. The 1992 data set is more limited, with lightcurves only spanning 14 hr. Time-averaged photometry for both encounters is presented while variability and color are discussed relative to findings from NASA's EPOXI mission (which also provided photometric lighturves for Earth). The new Galileo/SSI data are used to further validate the Virtual Planetary Laboratory 3D spectral Earth model, which often serves as a stand-in for true disk-integrated observations of our planet. The revived Galileo/SSI data for Earth is a testament to the ability of NASA's Planetary Data System to maintain data over decades-long timescales. The disk-integrated products derived from these data add to a very short list of calibrated and published whole-disk observations of the Pale Blue Dot.

It has been long believed that oblique and quasi-perpendicular configurations in supernova remnants (SNRs) were inefficient at injecting ions into diffusive shock acceleration (DSA), and that the highest energy Galactic cosmic rays (CRs) must come from parallel or quasi-parallel shocks. However, recent 3D kinetic simulations have shown that high-obliquity shocks can successfully energize particles and produce amplified magnetic fields in the upstream. We aim to investigate the maximum energy particles it is possible to produce in oblique and quasi-perpendicular shocks and whether they are capable of triggering the non-resonant hybrid instability (NRHI). We present a novel setup for hybrid simulations of non-relativistic shocks that use a "faux shock" boundary condition instead of a real shock to significantly reduce the computational cost and that can be run for long enough to study the late-time behaviors of these systems. Our results show that it may be possible for oblique and quasi-perpendicular shocks to transition from early periods of shock drift acceleration (SDA) into DSA at later times, giving particles a brief period of rapid acceleration followed by a long-duration, self-sustaining period of slower energy growth. Furthermore, we find evidence that the NRHI is triggered in the upstream at late times. Oblique and quasi-perpendicular shocks may be an important contributor to high energy CRs, potentially even responsible for the knee in the CR energy spectrum.

We report the serendipitous discovery of a large edge-on protoplanetary disk associated with the infrared source IRAS 23077+6707. The disk's apparent size in the Pan-STARRS (PS1) images is ~11", making this one of the largest known disks on the sky. It is likely a young system, still surrounded by the envelope which is very faint but still visible in the PS1 images in the northern part (alternatively this structure could be filaments from the disk itself). We use the PS1 magnitudes and other available photometric data to construct the spectral energy distribution (SED) of the disk. An optical spectrum indicates that the obscured star is hot, most likely late A. We adopt a distance of 300 pc for this object based on Gaia DR3 extinctions. We model the system using the HOCHUNK3D radiative transfer software and find that the system is consistent with a hot star of effective temperature 8000 K surrounded by a disk of size 1650 AU and mass 0.2 M_solar at inclination 82 degrees.

Globular clusters (GCs), as old as our Galaxy, constantly lose their members to the field as they cross through the Milky Way (MW). These GC escaped stars (or escapees) are suggested to contribute significantly to the MW halo. If a star has left the host GC a long time ago, chemical finger prints, e.g., N enrichment, may reveal its origin. In this work, we aim to establish dynamical connections between N-rich field stars recently identified by LAMOST and the existing MW GCs. By constructing the full action distribution, and combining with metallicity, we found 29 potential GC progenitors for 15 N-rich field stars. Particularly, some of them may be related to MW accretion events. On the other hand, if a star recently left its host GC via tidal evaporation, it still maintain the kinematic properties of the cluster. Here we identify extra-tidal candidates based on their spatial locations, proper motions (PMs), and their position on color-magnitude-diagrams (CMDs). We successfully identified more than 1600 extra-tidal candidates in the vicinity of six Gaia-Enceladus (GE)-related GCs (i.e., NGC 1851, NGC 1904, NGC 6205, NGC 6341, NGC 6779, NGC 7089). The density map of the extra-tidal candidates is confirmed to be an efficient way to find extra-tidal structures. The possible two density peaks at opposite directions of the inner boundary is a good indicator for long stellar stream. Among 95 extra-tidal candidates with spectroscopic radial velocities and metallicity, 54 of them are confirmed to be GC escaped stars, as they share similar properties as host GCs. These extra-tidal candidates are ideal targets for follow-up spectroscopic observation, as it greatly improves the scientific outcome. Once statistically significant number of spectroscopic radial velocities and metallicities are available, the GC dynamical evolution (e.g., mass loss, rotation) can be carefully investigated.

A scale-independent energy-momentum squared gravity (EMSG) allows different gravitational couplings for different types of sources and has been proven to have interesting implications in cosmology. In this paper, the Big Bang Nucleosynthesis (BBN) formalism and the latest observational constraints are being used in order to extract constraints on this class of modified gravity models. The model has been constrained using the light element reaction rates. Using the tight constraint from BBN on the correction term in the Friedman equation due to EMSG, we find a significant deviation in the matter power spectrum of the cosmic microwave background (CMB) even for a very small allowed value of the modification parameter, thus pointing out the sensitivity of the matter power spectrum at the small scales.

Context. One of the most fundamental scaling relations in galaxies is observed between metallicity and stellar mass -- the mass-metallicity relation (MZR) -- although recently a stronger dependence of the gas-phase metallicity with the galactic gravitational potential ($\Phi \rm ZR$) has been reported. Further dependences of metallicity on other galaxy properties have been revealed, with the star formation rate (SFR) being one of the most studied and debated secondary parameters in the relation (the so-called fundamental metallicity relation). Aims. In this work we explore the dependence of the gas-phase metallicity residuals from the MZR and $\Phi \rm ZR$ on different galaxy properties in the search for the most fundamental scaling relation in galaxies. Methods. We applied a random forest regressor algorithm on a sample of 3430 nearby star-forming galaxies from the SDSS-IV MaNGA survey. Using this technique, we explored the effect of 147 additional parameters on the global oxygen abundance residuals obtained after subtracting the MZR. Alternatively, we followed a similar approach with the metallicity residuals from the $\Phi \rm ZR$. Results. The stellar metallicity of the galaxy is revealed as the secondary parameter in both the MZR and the $\Phi \rm ZR$, ahead of the SFR. This parameter reduces the scatter in the relations $\sim 10-15\%$. We find the 3D relation between gravitational potential, gas metallicity, and stellar metallicity to be the most fundamental metallicity relation observed in galaxies.

We report the test results of several independent foreground-cleaning pipelines used in the Ali CMB Polarization Telescope experiment (AliCPT-1), a high-altitude CMB imager in the Northern hemisphere with thousands of detectors dedicated to the search for a primordial CMB polarization $B$-mode signature. Based on simulated data from 4 detector modules and a single season of observation, which we refer to as DC1 data, we employ different and independent pipelines to examine the robustness and effectiveness of the estimates on foreground parameters and the primordial $B$-mode detection. The foreground-cleaning strategies used in the pipelines include the parametric method of template fitting (TF) and the non-parametric methods of the constrained internal linear combination (cILC), the analytical blind separation (ABS), and the generalized least squares (GLS). We examine the impact of possible foreground residuals on the estimate of the CMB tensor-to-scalar ratio ($r$) for each pipeline by changing the contamination components in the simulated maps and varying the foreground models and sky patches for various tests. According to the DC1 data with the simulation input value $r_{\rm true}=0.023$, the foreground residual contamination levels in the TF/ABS/cILC/GLS pipelines are well within the corresponding statistical errors at the $2\sigma$ level. For a selected patch with relatively stronger foreground contamination, all of the proposed pipelines perform robustly in testing. Furthermore, by utilizing the tension estimator, which helps identify significant residual foreground contamination in the detection of the primordial $B$-mode signal by quantifying the discrepancy between various $r$ measurements, we conclude that the presence of small foreground residuals does not lead to any significant inconsistency in the estimation of $r$.

Directional wavelet dictionaries are hierarchical representations which efficiently capture and segment information across scale, location and orientation. Such representations demonstrate a particular affinity to physical signals, which often exhibit highly anisotropic, localised multiscale structure. Many physically important signals are observed over spherical domains, such as the celestial sky in cosmology. Leveraging recent advances in computational harmonic analysis, we design new highly distributable and automatically differentiable directional wavelet transforms on the $2$-dimensional sphere $\mathbb{S}^2$ and $3$-dimensional ball $\mathbb{B}^3 = \mathbb{R}^+ \times \mathbb{S}^2$ (the space formed by augmenting the sphere with the radial half-line). We observe up to a $300$-fold and $21800$-fold acceleration for signals on the sphere and ball, respectively, compared to existing software, whilst maintaining 64-bit machine precision. Not only do these algorithms dramatically accelerate existing spherical wavelet transforms, the gradient information afforded by automatic differentiation unlocks many data-driven analysis techniques previously not possible for these spaces. We publicly release both S2WAV and S2BALL, open-sourced JAX libraries for our transforms that are automatically differentiable and readily deployable both on and over clusters of hardware accelerators (e.g. GPUs & TPUs).

We have re-observed $\rm\sim$40 low-inclination, star-forming galaxies from the MaNGA survey ($\upsigma\sim65$~\kms) at $\sim$6.5 times higher spectral resolution ($\upsigma\sim10$~\kms) using the HexPak integral field unit on the WIYN 3.5m telescope. The aim of these observations is to calibrate MaNGA's instrumental resolution and to characterize turbulence in the warm interstellar medium and ionized galactic outflows. Here we report the results for the H$\rm\upalpha$ region observations as they pertain to the calibration of MaNGA's spectral resolution. Remarkably, we find that the previously-reported MaNGA line-spread-function (LSF) Gaussian width is systematically underestimated by only 1\%. The LSF increase modestly reduces the characteristic dispersion of HII regions-dominated spectra sampled at 1-2 kpc spatial scales from 23 to 20 km s$^{-1}$ in our sample, or a 25\% decrease in the random-motion kinetic energy. This commensurately lowers the dispersion zeropoint in the relation between line-width and star-formation rate surface-density in galaxies sampled on the same spatial scale. This modest zero-point shift does not appear to alter the power-law slope in the relation between line-width and star-formation rate surface-density. We also show that adopting a scheme whereby corrected line-widths are computed as the square root of the median of the difference in the squared measured line width and the squared LSF Gaussian avoids biases and allows for lower SNR data to be used reliably.

We present our past and current long-term monitoring program of water masers in the circumstellar envelopes of evolved stars, augmented by occasional interferometric observations. Using as example the Mira-variable U Her, we identify three types of variability: periodic (following the optical variation), long-term (years-decades) and short-term irregular (weeks-months). We show there are regions in the maser shell where excitation conditions are favourable, which remain stable for many years. Lifetimes of maser clouds in the wind-acceleration zone are of the order of up to a few years. Much longer lifetimes are found for the peculiar case of a maser cloud outside that zone (as in RT Vir), or in some cases where the motion of spectral features can be followed for the entire 2 decade monitoring period (as in red supergiant VX Sgr).

Many quenched galaxies discovered in the early Universe by \textit{JWST} raise fundemental question s on when and how these galaxies became quiescent. Making use of the latest version of the semi-analytic model GAEA that provides good agreement with the observed quenched fractions up to $z\sim 3$, we make predictions for the expected fractions of quiescent galaxies up to $z\sim 7$ and analyze the main quenching mechanism. We find that in a simulated box of $685~{\rm Mpc}$ on a side, the first quenched massive ($M_{\star} \sim 10^{11} {\rm M}_{\odot}$), Milky Way mass, and low mass ($M_{\star} \sim 10^{9.5} {\rm M}_{\odot}$ ) galaxies appear at $z\sim 4.5$, $z\sim 6.2$, and before $z = 7$. Most quenched galaxies identified at early redshifts remain quenched for more than 1 Gyr. Independently of galaxy stellar mass, the dominant quenching mechanism at high redshift is accretion disk feedback (quasar winds) from a central massive black hole, which is triggered by mergers in massive and Milky Way mass galaxies, and by disk instabilities in low-mass galaxies. Environmental stripping become increasingly more important at lower redshift.

This work aims to investigate the behaviour of the lithium abundance in stars with and without detected planets. Our study is based on a sample of 1332 FGK main-sequence stars with measured lithium abundances, for 257 of which planets were detected. Our method reviews the sample statistics and is addressed specifically to the influence of tides and orbital decay, with special attention to planets on close orbits, whose stellar rotational velocity is higher than the orbital period of the planet. In this case, tidal effects are much more pronounced. The analysis also covers the orbital decay on a short timescale, with planets spiralling into their parent star. Furthermore, the sample allows us to study the relation between the presence of planets and the physical properties of their host stars, such as the chromospheric activity, metallicity, and lithium abundance. In the case of a strong tidal influence, we cannot infer from any of the studies described that the behaviour of Li differs between stars that host planets and those that do not. Our sample includes stars with super-solar metallicity ([Fe/H]>0.15 dex) and a low lithium abundance (A(Li) <1.0 dex). This enabled us to analyse scenarios of the origin and existence of these stars. Considering the possible explanation of the F dip, we show that it is not a plausible scenario. Our analysis is based on a kinematic study and concludes that the possible time that elapsed in the travel from their birth places in the central regions of the Galaxy to their current positions in the solar neighbourhood is not enough to explain the high lithium depletion. It is remarkable that those of our high-metallicity low-lithium stars with the greatest eccentricity (e>0.2) are closest to the Galactic centre. A dedicated study of a set of high-metallicity low-Li stars is needed to test the migration-depletion scenario.

We report the first detection in space of a complete sample of nine doubly substituted isotopologues of HCCCN towards the cyanopolyyne peak of TMC-1 using observations of the QUIJOTE line survey taken with the Yebes 40 m telescope. We detected D13CCCN, DC13CCN, DCC13CN, DCCC15N, H13C13CCN, H13CC13CN, HC13C13CN, HCC13C15N, and HC13CC15N through their J=4-3 and J=5-4 lines in the 7 mm window. In addition, we present an extensive analysis of the emission of HCCCN and its singly substituted isotopologues through a large velocity gradient model of the lines detected at 7 mm and 3 mm using the Yebes 40 m and the IRAM 30 m telescopes, respectively. The derived column densities for all the isotopologues are consistent in the two spectral bands for an H2 volume density of 1e4 cm-3 and a kinetic temperature of 10 K. Whereas we observed a 13C fractionation for HCC13CN and other double isotopologues with a 13C atom adjacent to the nitrogen atom, we derived similar C/13C abundance ratios for the three 13C substituted species of DCCCN. This suggests additional chemical discrimination for deuterated isotopologues of HCCCN. Finally, we present the spatial distribution of the J=4-3 and J=5-4 lines from the singly substituted species observed with the Yebes 40 m telescope. The emission peak of the spatial distribution of DCCCN appears to be displaced by 40'' with respect to that of HCCCN and the 13C and 15N isotopologues. In addition to a different formation route for the deuterated species, we could also expect that this differentiation owing to the deuterium fractionation is more efficient at low temperatures, and therefore, that deuterated species trace a colder region of the cloud.

Ionization is important for magnetohydrodynamics and chemistry in protoplanetary disks but known ionization sources are often weak along the midplane. We present, for the first time, data from a laboratory experiment, where we measure ions from colliding mm-basalt grains emitted into the surrounding gas phase. This positive detection implies that very basic collisions in early phases of planet formation are sources of ionization. The midplane of protoplanetary disks might be ionized despite the lack of intense radiation sources.

The accretion and ejection of mass in pre-main sequence (PMS) stars are key processes in stellar evolution as they shape the stellar angular momentum transport necessary for the stars' stability. Magnetospheric accretion onto classical T Tauri stars and low-mass PMS stars has been widely studied in the single-star case. This process can not be directly transferred to PMS binary systems, as tidal and gravitation effects, and/or accretion from a circumbinary disc (with variable separation of the components in the case of eccentric orbits) are in place. This work examines the accretion process of two PMS eccentric binaries, DQ Tau and AK Sco, using high-resolution spectropolarimetric time series. We investigate how magnetospheric accretion can be applied to these systems by studying the accretion-related emission lines and the magnetic field of each system. We discover that both systems are showing signs of magnetospheric accretion, despite their slightly different configurations, and the weak magnetic field of AK Sco. Furthermore, the magnetic topology of DQ Tau A shows a change relative to the previous orbital cycle studied: previously dominated by the poloidal component, it is now dominated by the toroidal component. We also report an increase of the component's accretion and the absence of an accretion burst at the apastron, suggesting that the component's magnetic variation might be the cause of the inter-cycle variations of the system's accretion. We conclude on the presence of magnetospheric accretion for both systems, together with gravitational effects, especially for AK Sco, composed of more massive components.

Gamma-ray bursts (GRBs) are particle acceleration sites that can emit photons in the very high-energy (VHE) domain through non-thermal processes. From 2004 until 2018, the current generation of Imaging Atmospheric Cherenkov Telescopes (IACTs) did not detect any GRB in the VHE domain. However, from 2018 to 2020, five detections have been reported. In this work, we try to solve the case of the missing VHE GBRs prior to 2018. We aim to identify GRBs that might have eluded VHE detection in the past years by the H.E.S.S., MAGIC, and VERITAS IACTs. To do so, we study GRBs with known redshift detected by \emph{Swift} from 2004 until June 2022. We first identify all GRBs that could have been observed by these IACTs since 2004, considering observation conditions and visibility constraints. We assume a relation between the X-rays and the VHE gamma rays based on the VHE GRBs detected to date and combine this with the redshift measurements, instrument response information, and observation conditions to predict the observed VHE gamma-ray flux from the \emph{Swift}-XRT measurements. We report findings on 12 bright low-redshift GRBs that we identify as most likely to have been detected in the VHE domain by current IACTs. The rate of IACT-detectable GRBs with ideal observation conditions is $<$1 VHE GRB per year with the current configuration. With its lower energy threshold and higher sensitivity, this rate increases to $\sim$4 VHE GRBs per year with CTA.

Context. The model of disc fragmentation due to gravitational instabilities offers an alternate formation mechanism for gas giant planets, especially those on wide orbits. Aims. Our goal is to determine the 3D structure of disc-instability protoplanets and to examine how this relates to the thermal physics of the fragmentation process. Methods. We modelled the fragmentation of gravitationally unstable discs using the SPH code PHANTOM, and followed the evolution of the protoplanets formed through the first and second-hydrostatic core phases (up to densities 1e-3 g/cm3). Results. We find that the 3D structure of disc-instability protoplanets is affected by the disc environment and the formation history of each protoplanet (e.g. interactions with spiral arms, mergers). The large majority of the protoplanets that form in the simulations are oblate spheroids rather than spherical, and they accrete faster from their poles. Conclusions. The 3D structure of disc-instability protoplanets is expected to affect their observed properties and should be taken into account when interpreting observations of protoplanets embedded in their parent discs.

Galactic winds probe how stellar feedback regulates the mass and metallicity of galaxies. Observations show that galactic winds are multiphase and magnetised. In the local Universe, the dense phase is traced by emission and absorption lines, which reveal the presence of fast-moving clouds embedded in hot streams. Simulations tell us that magnetic fields can shield such clouds and delay their disruption, but there is little discussed on their observational effects. Using 3D MHD simulations, we study the influence of two orientations of the magnetic field (aligned and transverse) on the cloud morphology, temperature and density structure, mixing fraction, ion kinematics, column densities, and absorption spectra. We study supersonic wind-cloud systems with radiative processes, and develop a framework to extract ion column density maps and synthetic absorption spectra. The framework relies on studying ion populations and creating down-the-barrel spectra via an interface that links our PLUTO simulations to TRIDENT using YT, CLOUDY, and STARBURST99. We find that the transverse magnetic field makes the cloud asymmetric, shields and protects dense cold gas, and reduces mixing fractions compared to the aligned case. Ions can reach higher velocities in the transverse field case. The imprints of the initial orientation of the field on the synthetic spectra are: in the cold phase we find no signature of C ii and Si ii when the field is aligned, in the intermediate phase traced by C iv and Si iv we find broader lines in the transverse case, and in the warm phase we find deeper lines for O vi and N v in the aligned case, but they are less sensitive to the field orientation. Magnetic fields significantly affect the absorption spectra of cold clouds. Intermediate ions are the most sensitive to the magnetic field orientation and can potentially yield information about magnetic field topology.

We compiled a list of more than 3500 eclipsing binaries located in and near the Northern Continuous Viewing Zone (NCVZ) of the TESS space telescope that have a sufficient amount of TESS photometry to search for additional hidden components in these systems. We obtained the TESS light curves of all targets in an automated way applying convolution-aided differential photometry on the TESS Full-Frame Images from all available sectors up to Sector 60. Using a new self-developed Python GUI, we visually vetted all of these light curves, determined the eclipsing periods of the objects and calculated their eclipse timing variations (ETVs). The ETV curves were used in order to search for nonlinear variations that could be attributed to a light travel time effect (LTTE) or dynamical perturbations caused by additional components in these systems. We pre-selected 351 such candidates and tried to model their ETVs with the analytic formulae of pure LTTE or the combination of LTTE and dynamical perturbations. In total we could fit a model solution for the ETVs of 135 hierarchical triple candidates in which 10 systems were already known in the literature and the remainder of the 125 systems are new discoveries. Among these systems, there are some more noteworthy ones, such as five tight triples very close to their dynamical stability limit with a period ratio of less than 20 and three newly discovered triply eclipsing triples. We point out that dynamical perturbations are occurring in GZ Dra, which turns out to be a triple. We also made a comparison of the distributions of some orbital parameters coming from our solutions with those from the Kepler sample derived by Borkovits et al. (2016). Finally, we checked the correlations between the available parameters for systems that have Gaia Non-Single Star orbital solutions with those from our ETV solutions. (Abridged)

The long-term relaxation of rotating globular clusters is investigated through an extension of the orbit averaged Chandrasekhar non-resonant formalism. A comparison is made with the long-term evolution of the distribution function in action space, measured from averages of sets of $N$-body simulations up to core collapse. The impact of rotation on in-plane relaxation is found to be weak. In addition, we observe a clear match between theoretical predictions and $N$-body measurements. For the class of rotating models considered, we find no strong gravo-gyro catastrophe accelerating core collapse. Both kinetic theory and simulations predict a reshuffling of orbital inclinations from overpopulated regions to underpopulated ones. This trend accelerates as the amount of rotation is increased. Yet, for orbits closer to the rotational plane, the non-resonant prediction does not reproduce numerical measurements. We argue that this mismatch stems from these orbits' coherent interactions, which are not captured by the non-resonant formalism that only addresses local deflections.

To disentangle the different competing physical processes at play in Galactic evolution, a detailed chrono chemicalkinematical, and dynamical characterisation of the disc bimodality is necessary, including high number statistics. Here we make use of an extremely precise subsample of the Gaia DR3 GSP-Spec catalogue of stellar chemophysical parameters. The selected database is composed of 408 800 stars with a median uncertainty of 10 K, 0.03 and 0.01 dex in Teff , log(g) and [M/H], respectively. The stellar parameter precision allows to break the age-metallicity degeneracy of disc stars. For the first time, the disc bimodality in the Kiel diagramme of giant stars is observed, getting rid of interstellar absortion issues. This bimodality produces double Red Giant Branch sequences and Red Clump features for mono-metallicity populations. A comparison with BaSTI isochrones allows to demonstrate that an age gap is needed to explain the evolutionary sequences separation, in agreement with previous age-metallicity relations obtained using Main-Sequence Turn Off stars. A bimodal distribution in the stellar mass-[alpha/Fe] plane is observed at constant metallicity. Finally, a selection of stars with [M/H]=0.45 \pm 0.03 dex shows that the most metal-rich population in the Milky Way disc presents an important proportion of stars with ages in the range 5-13 Gyr. This old, extremely metal-rich population is possibly a mix of migrated stars from the internal Galactic regions, and old disc stars formed before the last major merger of the Milky Way. The Gaia GSP-Spec Kiel diagrammes of disc mono-abundance stellar populations reveal a complex, non linear age-metallicity relation crafted by internal and external processes of Galactic evolution. Their detailed analysis opens new opportunities to reconstruct the puzzle of the Milky Way disc bimodality.

We present the second release of the Meta-catalogue of X-Ray detected Clusters of galaxies (hereafter MCXC-II). The MCXC-II has been compiled from publicly available ROSAT All Sky Survey-based (NORAS, REFLEX, BCS, SGP, NEP, MACS, CIZA, RXGCC) and serendipitous (160SD, 400SD, SHARC, WARPS, and EMSS) cluster catalogues. Redshifts were systematically checked and updated if necessary, and additional redshift information (type and origin) added. The X-ray data were standardised to an overdensity of 500 using a new procedure based on the use of the original flux and aperture measurements available in the input catalogues. MCXC-II contains 2221 entries, essentially completing the census of ROSAT cluster detections by including objects from the REFLEX-II and RXGCC surveys, in addition to providing a complete and fully-homogenised sub-catalogue of all published MACS clusters. Duplicate entries from overlaps between the survey areas of the individual input catalogues were carefully handled. For each cluster the MCXC-II provides three identifiers, a redshift, coordinates, membership in the original catalogue, and standardised [0.1-2.4] keV band luminosity $L_{500}$, total mass $M_{500}$, and radius $R_{500}$. Uncertainties on $L_{500}$ were computed from the flux errors in the original catalogues. MCXC-II additionally furnishes information on overlaps between the input catalogues, gives the luminosity and its uncertainty when measurements from different surveys are available, and provides notes on individual objects.

Context. The galactic winds of starburst galaxies (SBGs) give rise to remarkable structures on kiloparsec scales. However, the evolution and shape of these giant wind bubbles, as well as the properties of the shocks they develop, are not yet fully understood. Aims. We aim to understand what shapes the galactic winds of SBGs, with a particular focus on the role of large-scale magnetic fields in the dynamical evolution of galactic wind-inflated bubbles. In addition, we aim to explore where the conditions for efficient particle acceleration are met in these systems. Methods. We performed magnetohydrodynamic simulations with the AMRVAC code (Adaptive Mesh Refinement Versatile Advection Code) with various configurations of the galactic medium density profile and magnetization. Results. We observe that the large-scale magnetic field, in which galactic winds expand, can impact the structure and evolution of inflated bubbles. However, the typical structures observed in starburst galaxies, such as M82, cannot be solely explained by the magnetic field structures that have been considered. This highlights the importance of other factors, such as the galactic disk, in shaping the galactic bubble. Furthermore, in all the magnetized cases we investigated, the forward wave resulting from the expanding bubbles only results in compression waves, whereas the wind termination shock features high Mach numbers, making it a promising site for diffusive shock acceleration up to $\sim 10^{2}$ PeV. The synthetic X-ray images generated from our models reveal an envelope surrounding the bubbles that extends up to 2 kpc, which could correspond to the polarized emission observed from planar geometry in M82, as well as a large structure inside the bubble corresponding to the shocked galactic wind.

Due to universality and attractor properties, $\alpha$-attractor quintessential inflation establishes direct relations between inflationary observables such as the scalar tilt $n_s$ and the tensor-to-scalar ratio $r$, and late-time dark energy equation of state parameters $w_0$ and $w_a$. In this work, we examine three different physically motivated regimes, considering complete freedom in the parameter $\alpha$, models inspired by supergravity where $\alpha$ takes on values up to $\alpha=7/3$, and Starobinsky inflation ($\alpha=1$). We investigate the consistency and constraints imposed by Cosmic Microwave Background measurements from the Planck satellite, B-mode polarization data from the BICEP/Keck collaboration, and low-redshift observations. Additionally, we consider small-scale CMB measurements released by the Atacama Cosmology Telescope, which give results approaching the Harrison-Zel'dovich spectrum ($n_s \approx 1$). Here $\alpha$-attractors lead to an improved fit over $\Lambda$CDM. For the large-scale CMB measurements, $\alpha\gtrsim2$ models can provide equally good fits as $\Lambda$CDM.

Deep radio surveys of extragalactic legacy fields trace a large range of spatial and brightness temperature sensitivity scales, and therefore have differing biases to radio-emitting physical components within galaxies. This is particularly true of radio surveys performed at less than 1 arcsec angular resolutions, and so robust comparisons are necessary to better understand the biases present in each survey. We present a multi-resolution and multi-wavelength analysis of the sources detected in a new Very Long Baseline Array (VLBA) survey of the CANDELS GOODS-North field. For the 24 VLBA-selected sources described in Paper I, we augment the VLBA data with EVN data, ~0.1-1 arcsecond angular resolution data provided by VLA and e-MERLIN. This sample includes new AGN detected in this field, thanks to a new source extraction technique that adopts priors from ancillary multi-wavelength data. The high brightness temperatures of these sources (Tb > 10^6 K) confirm AGN cores, that would often be missed or ambiguous in lower-resolution radio data of the same sources. Furthermore, only 15 sources are identified as 'radiative' AGN based on available X-ray and infrared constraints. By combining VLA and VLBA measurements, we find evidence that the majority of the extended radio emission is also AGN dominated, with only 3 sources with evidence for extended potentially star-formation dominated radio emission. We demonstrate the importance of wide-field multi-resolution (arcsecond-milliarcsecond) coverage of the faint radio source population, for a complete picture of the multi-scale processes within these galaxies.

Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ formation, stellar radiation, sink particles, and turbulence. We find that, at early times in the cloud evolution, H$_2$ formation rapidly outpaces dissociation and molecular clouds build their mass from the atomic reservoir in their environment. Rapid H$_2$ formation is also associated with a higher early star formation rate. For the clouds studied here, H$_2$ is strongly out of chemical equilibrium during the early stages of cloud formation but settles into a bursty chemical steady-state about 2 Myrs after the first stars form. At the latest stage of cloud evolution, dissociation outweighs formation and the clouds enter a dispersal phase. We discuss how theories for the molecular cloud lifecycle and the star formation efficiency may be distinguished with observational measurements of H$_2$ fluorescence with a space-based high-resolution FUV spectrometer, such as the proposed Hyperion and Eos NASA Explorer missions. Such missions would enable measurements of the H$_2$ dissociation and formation rates, which we demonstrate can be connected to different phases in a molecular cloud's star-forming life, including cloud building, rapidly star-forming, H$_2$ chemical equilibrium, and cloud destruction.

Self-interacting dark matter (SIDM) theory predicts that dark matter halos experience core-collapse, a process where the halo's inner region rapidly increases in density and decreases in size. The N-body simulations used to study this process can suffer from numerical errors when simulation parameters are selected incorrectly. Optimal choices for simulation parameters are well studied for cold dark matter (CDM), but are not deeply understood when self-interactions are included. In order to perform reliable N-body simulations and model core-collapse accurately we must understand the potential numerical errors, how to diagnose them, and what parameter selections must be made to reduce them. We use the $\texttt{Arepo}$ N-body code to perform convergence tests of core-collapsing SIDM halos across a range of halo concentrations and SIDM cross-sections, and quantify potential numerical issues related to mass resolution, timestep size, and gravitational softening length. Our tests discover that halos with fewer than $10^5$ simulation particles, a resolution typically not met by subhalos in N-body simulations, suffer from significant discreteness noise that leads to variation and extreme outliers in the collapse rate. At our lowest resolution of $N=10^4$ particles, this collapse time variation can reach as high as 20%. At this low resolution we also find a bias in collapse times and a small number of extreme outliers. Additionally, we find that simulations which run far beyond the age of the Universe, which have been used to calibrate SIDM gravothermal fluid models in previous work, have a sensitivity to the timestep size that is not present in shorter simulations or simulations using only CDM. Our work shows that choices of simulation parameters that yield converged results for some halo masses and SIDM models do not necessarily yield convergence for others.

The asymmetry of line profiles, i.e., the secondary component, is crucial to understanding the energy release of coronal holes (CH), quiet sun (QS), and bright points (BPs). We investigate the asymmetry of Si IV 1393.75 {\AA} of the transition-region (TR) line recorded by Interface Region Imaging Spectrometer (IRIS) and co-spatial-temporal Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) data onboard Solar Dynamics Observatory (SDO) for three time series on 26 April 2015, 24 July 2014, 26 July 2014. Most asymmetric profiles are in the complex magnetic field regions of the networks. The asymmetric profiles are fitted with both single and double Gaussian models. The mean value of Doppler velocity of the second component is almost zero (with a significant standard deviation) in QS/CH, which may indicate that the physical process to trigger the secondary Gaussian originates at the formation height of Si IV. While the mean Doppler velocity from secondary Gaussian in BPs is around +4.0 km/s (redshifted). The non-thermal velocities of the secondary Gaussian in all three regions are slightly higher than the single Gaussian. The statistical investigation leads to the prevalence of blueshifted secondary components in QS/CH. However, secondary Gaussian in the BPs redshifted, i.e., the BPs redshift behavior could be interpreted due to the site of reconnection located above the formation height of the Si IV line. The peak intensity of the second component for all three regions is likely to follow a power law that is a signature of the small-scale flaring-like trigger mechanism.

The kinematics and chemical composition of stellar populations of different ages provide crucial information about the evolution of a galaxy. We aim to provide such information for IC 1613, an isolated, gas-rich, star-forming dwarf galaxy in the Local Group. We present here the results of a new spectroscopic study performed with MUSE, an integral-field spectrograph on the Very Large Telescope. We extracted from the data cubes more than 2000 sources from which we separated stellar objects for further spectroscopic analysis. The quality of the data set allowed us to obtain accurate classifications and line-of-sight velocities for about 800 stars. Our sample includes not only Red Giant Branch (RGB) and Main Sequence (MS) stars, but also a number of probable Be and C stars. We also obtained reliable metallicities for about 300 RGB stars. The kinematic analysis revealed for the first time the presence of stellar rotation with high significance. We found a general agreement with the velocity field of the neutral gas component, although the stars showed on average a larger velocity dispersion and slower rotation due to the asymmetric drift. When examining the kinematics of the different stellar components, MS stars appear to closely follow that of the gas, and the velocity dispersion seems to increase towards older stars. Chemical analysis of the RGB stars revealed mean properties comparable to those of other Local Group dwarf galaxies. Our work represents a step forward in understanding the internal processes governing the dynamical evolution of a low-mass galaxy.

A number of modified gravity theories (e.g., $f(R)$-gravity) lead to a Yukawa-like metric in the weak field limit which can be described by two Yukawa parameters, i.e., the strength $\kappa$ and the length scale $\lambda$. The S-stars, orbiting around the supermassive black hole in the Galactic Center, are unique probes to test these gravity theories in relatively strong gravitational field. The Newtonian Yukawa gravity potential or a simple approximation to the Yukawa metric was usually adopted in previous studies when using the orbital motion of S-stars to constrain such modified gravity theories, which may be not sufficiently accurate considering recent and future high resolution observations. In this paper, we first derive the Post-Newtonian (PN) Yukawa motion equation at the 2PN order, and then investigate the high order effects on the orbital motions by comparison with those from the Newtonian Yukawa gravity potential. We further obtain constraints on $\kappa$ by using the observations on the orbital motions of several S-stars (i.e., S2, S38, and S55). Our results show that the current observations of these stars are compatible with the General Relativity and $\kappa$ can be constrained to $|\kappa|<0.01$ with $95\%$ confidence if $\lambda\in(100,250)$\,AU. We also estimate the possible improvements (about an order of magnitude or more) to the constraints by future higher resolution observations and the inclusion of closer S-stars, such as S4716.

We report on the development of a novel pixel charge readout system, Grid Activated Multi-scale pixel readout (GAMPix), which is under development for use in the GammaTPC gamma ray instrument concept. GammaTPC is being developed to optimize the use of liquid argon time projection chamber technology for gamma ray astrophysics, for which a fine grained low power charge readout is essential. GAMPix uses a new architecture with coarse and fine scale instrumented electrodes to solve the twin problems of loss of measured charge after diffusion, and high readout power. Fundamentally, it enables low noise and ultra low power charge readout at the spatial scale limited by diffusion in a time projection chamber, and has other possibly applications, including future DUNE modules.

We search for indirect signals of $\mathscr{O}$(keV) dark matter annihilating or decaying into $\mathscr{O}$(eV) dark photons. These dark photons will be highly boosted and have decay lengths larger than the Milky Way, and can be absorbed by neutrino or dark matter experiments at a rate dependent on the photon-dark photon kinetic mixing parameter and the optical properties of the experiment. We show that current experiments can not probe new parameter space, but future large-scale gaseous detectors with low backgrounds (i.e. CYGNUS, NEXT, PANDAX-III) may be sensitive to this signal when the annihilation cross section is especially large.

Extragalactic and galactic cosmic rays scatter with the cosmic neutrino background during propagation to Earth, yielding a flux of relic neutrinos boosted to larger energies. If an overdensity of relic neutrinos is present in galaxies, and neutrinos are massive enough, this flux might be detectable by high-energy neutrino experiments. For a lightest neutrino of mass $m_{\nu} \sim 0.1$ eV, we find an upper limit on the local relic neutrino overdensity of $\sim 10^{13}$ and an upper limit on the relic neutrino overdensity at TXS 0506+056 of $\sim 10^{10}$. Future experiments like GRAND or IceCube-Gen2 could improve these bounds by orders of magnitude.

We present results from the first experimental demonstration of a tunable thin-shell axion haloscope. This novel geometry decouples the overall volume of the cavity-based resonator from its resonant frequency, thereby evading the steep sensitivity degradation at high-frequencies. An aluminum $2.6$ L ($41$ $\lambda^3$) prototype which tunes from $7.1$ to $8.0$ GHz was fabricated and characterized at room temperature. An axion-sensitive, straightforwardly tunable $\mathrm{TM}$$_{010}$ mode is clearly identified with a room temperature quality factor, $Q$, of $\sim$$5,000$. The on-resonance $E$-field distribution is mapped and found to agree with numerical calculations. Anticipating future cryogenic operation, we develop an alignment protocol relying only on rf measurements of the cavity, maintaining a form factor of $0.57$ across the full tuning range. These measurements demonstrate the feasibility of cavity-based haloscopes with operating volume $V\gg\lambda^3$. We discuss plans for future development and the parameters required for a thin-shell haloscope exploring the post-inflationary axion parameter space ($\sim$$4$ to $\sim$$30$ GHz) at DFSZ sensitivity.

Following the 2019 release by the Event Horizon Telescope Collaboration of the first pictures of a supermassive black hole, there has been an explosion of interest in black hole images, their theoretical interpretation, and their potential use in tests of general relativity. The literature on the subject has now become so vast that an introductory guide for newcomers would appear welcome. Here, we aim to provide an accessible entry point to this growing field, with a particular focus on the black hole "photon ring": the bright, narrow ring of light that dominates images of a black hole and belongs to the black hole itself, rather than to its surrounding plasma. Far from an exhaustive review, this beginner's guide offers a pedagogical review of the key basic concepts and a brief summary of some results at the research frontier.

The CSC (cryogenic scintillating calorimeter) technology devoted to rare event searches is reaching the sensitivity level required for the hunt of dark matter-electron scatterings. Dark matter-electron interactions in scintillating targets are expected to stimulate the emission of single photons, each of energy equal to the target electronic band gap. The electronic band gap in scintillators like NaI/GaAs is of O(eV). The search for this signal can be done by an array of cryogenic light detectors with eV/sub-eV energy resolution. In this work, we describe the detection principle, the detector response and the envisioned detector design to search for dark matter interacting with electrons via the measurement of the scintillation light at millikelvin. First sensitivity projections are provided, which show the potential of this research.

The Planetary Instrument for X-ray Lithochemistry (PIXL) is a rasterable focused-beam X-ray fluorescence (XRF) spectrometer mounted on the arm of National Aeronautics and Space Administration's (NASA) Mars 2020 Perseverance rover. To ensure that PIXL would be capable of performing accurate in-flight compositional analysis of martian targets, in situ, an elemental calibration was performed pre-flight on the PIXL flight instrument in a simulated martian environment. The details of this calibration, and implications for measuring unknown materials on Mars are the subjects of this paper. The major goals of this calibration were both to align the spectrometer to perform accurate elemental analysis and, to derive a matrix of uncertainties that are applied to XRF measurements of all elements in unknown materials. A small set of pure element and pure compound targets and geologically relevant reference materials were measured with the flight hardware in a simulated martian environment. Elemental calibration and quantifications were carried out using PIXL's XRF quantification software (PIQUANT). Uncertainties generated were implemented into the PIQUANT software version employed by the PIXL's data visualization software (PIXLISE). We outline in this work, a list of factors that impact micro-XRF accuracy, the methodology and steps involved in the calibration, details on the fabrication of the uncertainty matrix, instructions on the use and interpretations of the uncertainties applied to unknowns and an assessment on the limitations and areas open to future improvement as part of subsequent calibration efforts.

We derive the equations of motion of relativistic magnetohydrodynamics, as well as microscopic expressions for all of its transport coefficients, from the Boltzmann equation using the method of moments. In contrast to reference Phys. Rev. D 98(7) 2018, where a single component gas was considered, we perform our derivation for a locally neutral fluid composed of two massless particle species with opposite charges. We demonstrate that the magnetohydrodynamical equations of motion become dramatically different for this more realistic system. The shear-stress tensor no longer obeys a single differential equation; it breaks into three non-degenerate components with respect to the magnetic field, each evolving according to different dynamical equations. For large magnetic fields, we further show that the solution of this theory displays oscillatory behaviour that can no longer be described by an Israel-Stewart-like theory. Finally, we investigate the derived equations in a Bjorken flow scenario.